AFOSR-949 THE UNIVERSITY OF MICHIGAN DEPARTMENT OF ELECTRICAL ENGINEERING Electromagnetic Materials Laboratory Technical Note STUDY PROGRAM FOR A MAGNETIC TIMER UTILIZING THE PROPAGATION OF DOMAIN WALLS W. W. Raymond Do Mo Grimes ORA Project 04316 under contract witho AIR FORCE OFFICE OF SCIENTIFIC RESEARCH AIR RESEARCH AND DEVELOPMENT COMMAND CONTRACT NOo AF 49(638)-986 WASHINGTON, D oC administered through' OFFICE OF RESEARCH ADMINISTRATION ANN ARBOR July 1961

FOREWORD The following is the final report on a study program undertaken to investigate the possible use of controlled domain wall propagation in a magnetic material as a continuously variable time delay device. This program was jointly financed through DOFL Contract DA-99-186-502-ORD-945, and AFOSR Contract 49(638)-986. This report is jointly being used as a final report to DOFL and as a technical report to AFOSR. The equation of domain wall propagation is developed in terms of the fundamental parameters of the magnetic material. Conditions required for precisely controlled wall motion are discussed in terms of the shape and physical characteristics of the media. Those external parameters which influence propagation are elaborated on in view of how they may be controlled to optimize the shape and velocity of the propagating magnetic discontinuity. Available data from the literature on the range in domain wall velocity in various materials are discussed and are compiled in tabulated form. Included are some results of experimental data obtained in this laboratory during the period of this contract. Finally, several geometries are suggested for a suitable time delay system. In view of the high domain wall velocities, none of these systems are of use for time delays greater than several seconds. iii

TABLE OF CONTENTS Page LIST OF TABLES vii LIST OF FIGURES ix INTRODUCTION 1 PROPERTIES OF THE MAGNETIC MEDIA WHICH AFFECT WALL MOTION 2 AXIAL PROPAGATION OF A DOMAIN WALL IN A FERROMAGNETIC CYLINDER 8 SHAPE OF THE REVERSED NUCLEUS AND MOVING WALL 10 INDUCED VOLTAGE AND EDDY CURRENT DAMPING 12 WALL DAMPING DUE TO RELAXATION LOSS 15 DESIGN OF A MAGNETIC TIME DELAY SYSTEM 22 APPENDIX: THE 180~ MAGNETIC DOMAIN WALL 26 REFERENCE 47 v

LIST OF TABLES Table Page I. Data on Domain Wall Motion in Permalloys at Low Field 33 II. Data on Domain Wall Motion in,48 Ni - Fe Tapes (Ref. 12) 34 III. Data on Domain Wall Motion in Single Crystals 35 IV. Physical Data for.72 Ni -.28 Fe Permalloy 36 vii

LIST OF FIGURES Figure Page 1l Nucleation of a reversed domain in a uniformly magnetized cylinder. 37 2. Domain wall shape used for calculation of eddy current damping. 37 53 System for measurement of wall velocity, 38 4. Output voltage of search coils wound on the permalloy wire. 39 5. Axial wall velocity vs. applied field for a.020 cm permalloy wire under tension. 40 6. Wall velocity vs. applied field for a.020 cm permalloy wire in an eddy current damping tube. 41 7~ Hysteresis of Ho and R with change in tension in a.020 cm permalloy wire. 42 8. Axial wall velocity vs. applied field for a.0114 cm permalloy wire under tension. 43 9. Hysteresis of Ho and R with change in tension in a *014 cm annealed permalloy wire. 44 10. Hysteresis of Ho and R with change in tension in a.0114 cm permalloy wire after plastic stretch. 45 11. Motion of domain walls nucleated at the ends of main solenoid. 46 ix

INTRODUCTION The first concept of a magnetic time delay device would embody the main features of the classical experiments conducted by Sixtus and Tonks. A long thin wire or-tape of homogeneous magnetic material is placed in a long solenoid. The magnetization of the wire or tape is completely aligned axially by the solenoid field,.Hs. When the main field is removed the magnetization must remain aligned, which requires the wire to have a longitudinal anisotropy. The main field is now reduced to zero, reversed, and increased to a value H, less than the axial anisotropy field 2K/M so that the original direction of magnetization is unchanged. A local pulse field Hp in the direction of the reversed main field is now applied near one end of the solenoid, If the pulse field plus the solenoid field exceeds a critical value Ho a nucleus of reversed magnatization is formed. The pressure exerted by the main field on the domain wall of the nucleus causes it to grow. When the wall reaches the wire surface, it splits into two mirror halves, each propagating in opposite direction to the ends of the wire, thus progressively reversing the magnetization in the wire. One or more small sensing coils are would on the wire. The time delay is the elapsed time from the start of the pulse field until the domain passes through the sensing coils. The delay is adjusted by either changing the position of the sensing coil relative to the pulse coil, or changing the main solenoid field or both. 1

PROPERTIES OF THE MAGNETIC MEDIA WHICH AFFECT WALL MOTION The range of wall velocity is controlled by the solenoid field, however the range is rather limited (see Figures 5 and 8). For low applied fields the dependency of velocity on field is linear and may be expressed as: v = R(H - Ho) (1) where R is the wall mobility constant, to be developed, H is the applied field, and Ho is the critical nucleating field. The difference (H-Ho), exerts a pressure to move the wall, or in equivalence supplies the energy to all losses associated with wall motion. There is-both a.'lower and upper limit to the velocity for coherent wall motion. Below the lower limit the wall shape and area vary and it progresses in erratic jumps dependent on the local variation in structure of the magnetic media. Above the upper limit where the applied field now approaches the magnitude of the anisotropy field, 2K/M, nucleation may occur at random along the media, hence several walls may propagate simultaneously. If the field is further increased a state of non-uniform rotation of the magnetization takes place, this is the lower limit for fast switching devices. If the applied field exceeds the anisotropy field, reversal takes place by an almost simultaneous rotation of magnetization.2 In order to adequately discuss the factors which enter the wall mobility equation, a short summary is given of the energy and loss terms relevant to 2

the phenomena of wall motion. An appendix is included in which a brief development is given of the wall parameters for an ideal 180~ wall. See the review articles given in Reference 3 for a more extensive development. In the following energy density equations, constant terms which do not enter the wall equations are omitted. The exchange energy density is given by: f e x - )2 A[( 2 fex = A[(V 1) + (V 2) + (V a ] (2) where A is the exchange constant and the alphas are the direction cosines of the magnetization vector with respect to the principal axes of the cubic crystal. Equation (2) gives the energy per unit volume of interaction between an electron spin and its nearest neighbors which contribute to the magnetization. The energy is isotropic and has a minimum value when the spins are completely aligned. Values of A reported in the literature vary widely in different materials, however they are of the same order of magnitude, An average value of 2 x 10-6 ergs/cm3 will be used here. For a cubic crystal the crystalline anisotropy energy density is: fan = Kl[Qo|a + o + G + K2o] + Ka2[Cao]9 (5) which measures the energy required to rotate the magnetization from an easy to hard direction in the unstrained crystal. Positive anisotropy indicates the cube edges are the direction of easy magnetization. The anisotropy energy is extremely sensitive to temperature, changing at approximately the tenth of higher power of temperature, increasing as the temperature is lowered, 5

It also may vary as much as an order of magnitude depending on the state of order in the material. For permalloys at room temperature K2 is about 6 percent of the magnitude of K1 and may be neglected. Values of K1 range from 12 x 103 ergs/cm3 for 40 percent nickel content to zero at 74 percent nickel (Ni3Fe) in the quenched state. In permalloys which have been slowly cooled K1 varies from 10 x 10-3 at 40 percent nickel, zero at 63 percent nickel to -25 x 103 at 74 percent nickel. The range for furnace cooled materials lies in between the above values. Uniaxial anisotropy may be induced by annealing in a magnetic field. The induced anisotropy will vary from zero to 4 x 10s ergs/cm3 dependent on the cooling rate. K1 is taken to be -7 x 103 ergs/cm3 for the 1.72 Ni -.28 Fe material investigated here. The magnetoelastic or magnetostrictive energy density for the isotropic case, Xloo = 111 = X is' fme = 3/2 X T sin2 P (4) where X is the magnetostrictive strain constant, T the applied tensile stress, and P is the angle between the magnetization vector and the direction of applied stress. The magnetostrictive energy arises from the interaction between the magnetization and the mechanical strain of the crystal latticeo The energy is stored in the deformation strain of the lattice. Positive X indicates the magnetization will have an easy direction in the direction of applied tension or perpendicular to an applied compressive stresso Both X100oo and 111 are positive in the range from about 4 percent nickel to slightly less than 80 percent nickel and is independent of heat treatment 4

with the exception of the 68 to 80 percent range. For.72 Ni -.28 Fe, quenched state, the magnetostriction is nearly isotropic and. is taken to be about 15 x 106. Crystal anisotropy energy and magnetostrictive energy are closely related, the source of the latter being somewhat obscure. It is noted that a lattice of high symmetry has low anisotropy, and of low symmetry, high anisotropy as is evident from comparison of the respective values for iron and cobalt. The magnetostatic energy density: f = - H H (5) gives the interaction of the magnetization as a dipole with the applied field. H is taken as the effective field interior to the magnetic media, one component of which is the applied field. The energy of the magnetization in its own field is: fs - 2Hd M (6) where Hd is the demagnetizing field which arises when V * M ~ 0. The demagnetizing field is a function of the shape of the magnetic media and the change in M through the media. Theoretically it may be calculated from: Hd = 1- (M V) V dv (7) vol where Hd(r') and M(r), however the solution depends on a knowledge of M(r)o IThe demagnetizing field may be most conveniently stated in terms of the de5

magnetizing factor Hd = - N M (8) LOSS MECHANISMS When a field is applied to a magnetic media, the magnetization changes in direction by rotation about the field. The loss accompanying the rotation is attributed to spin-spin and spin-lattice interactions of a complex nature. With loss the motion of precession about the applied field is damped and the magnetization motion decays exponentially toward the direction of the field. 6 The process is given on a phenomenological basis by: -/a M = Mx -- (9 at 7M aty where = oe, ge Ym = mg o ge is the gyromagnetic ratio 2m 2m g = 2, is the Lande factor a= dimensionless loss factor. This form, suggested by Gilbert, gives a damping term proportional to the time rate of change of magnetization and results in a minimum time of reversal for a finite value of the damping constant. Due to the inertia of the spin system, the motion of magnetization in a moving wall will result in kinetic energy term proportional to the wall velocity squared. The inertia of the spin system summed through the moving wall is expressed as an equivalent wall mass, see appendix. 6

The change of magnetization with time in the moving wall induces an emf, and as a result eddy currents flow. The resultant power loss is a function of the shape of the domains and the magnetic material, conductivity, and wall velocity. The eddy current distribution is a boundary value problem, solved as follows: V2 J = 0 (10) everwhere except within the wall. The boundary conditions are: n- = 0 (11) at the conducting surface and + n xJ = 2poaMv (12) at the wall surface. n is the unit outward normal at the media surface and at the domain wall respectively. Both the relaxation damping of the magnetization and eddy current loss affect the wall mobility. Increase of a, the damping factor, and a the conductivity decrease the wall mobility, thus for low wall velocity in a given field it is desirable to increase damping and conductivity. Controlled domain wall motion is only observed in on permalloy tapes above one eighth mil in thickness. In ferrites and very thin nickel-iron films or tapes relaxation damping is dominant and eddy current damping negligible.7 7

AXIAL PROPAGATION OF A DOMAIN WALL IN A FERROMAGNETIC CYLINDER Since Sixtus and Tonks first studied the propagation of domain walls in permalloy wires, at least six other groups have extended the experiments (see References, 1, 8, 9, 10, 11, and 12). Extensive experimental data have been compiled on various compositions in the ranges of nickel composition from lO10 to 20% and 40%0 to 80%. The data show variation of velocity with composition, heat treatment, cold work, tension and torsion; however little or no theoretical development was given. The first order analysis is due to Doringl3 and a brief review is given here with some modification. The theoretical problem of the reversed domain nucleation, growth, and propagation may be approached rigorously by minimizing the total free energy of the magnetic configuration, with consideration of the loss terms. A variational solution would then yield the dynamic equations of motion of the magnetization. These equations may then be solved subject to the boundary conditions imposed at the media surface and through the moving region 14 of magnetization reversal, as suggested by W. F. Brown. The problem, however, is greatly complicated by two factors, namely: the demagnetizing field, Equation (7) and the retarding field produced by eddy currents. Both are a function of the magnetization spatial configuration in the zone of reversal and the magnetization configuration is dependent of these fields. The resultant equations are non-linear partial differential equations in space and time. No intensive effort has been made to solve these equations within the limited objective of this project, however they are of importance and a rigorous solution is required for a complete understanding of the problem. 8

A simplified approach requires a start from the experimentally determined shape of a reversed magnetic nucleus, and the assumption of an invariant shape as the reversed region grows to the media surface, thereafter propagating along the structure. Sixtus and Tonks measured the shape of reversed magnetic region "frozen" in wires. When the nucleus region did not extend to the wire surface, it formed a prolate spheroid with the long axis parallel to the wire axis. Doring 3 derived the critical ratio of major to minor axis for the nucleus which agreed well with experimental results. The assumption is made that the reversal nucleus exists as a uniformly magnetized prolate spheroid, separated from the surrounding magnetic media by a domain wallo Calculations from the integrated search coil voltage waveform show the shape of the moving discontinuity to be roughly one half of a spheroid with semiminor axis equal to the wire radius; it is to be noted here that other shapes may be made to fit equally well but do not simplify the problem of the demagnetizing field. The ratio of semi-major axis to semi-minor axis is normally from 800 to 2,000. The lower figure applies to cold drawn permalloy wire, and is almost independent of wire diameter Hydrogen annealed wires have an increase of the ratio from 900 to 2,000 as the wire diameter is decreased, as determined by the affect of the eddy current retarding field on the shape of the wall. 9

SHAPE OF THE REVERSED NUCLEUS AND MOVING WALL Assume the magnetic media to be a long cylinder of radius ao The cylinder is initially fully magnetized to the right, see Figure 1. A field is applied axially to the lefto At the center of the cylinder a domain will nucleate, which before it grows has an inner region uniformly magnetized to the left. The energy of the nucleus is U = 2oM(H - Hd) dv - aw ds (13) where reference is made to Equations (5) and (6) for the magnetostatic energy terms, and to the appendix Equation (43) for the wall energy term. If the nucleus is to grow the applied field must exceed some critical value Ho, if less the nucleus will disappear, hence d U ) 2ioMHo dv (14) The demagnetizing energy term may be rewritten taking Hx = -NxMx, Hz = -NzMz, where Nx and Nz are the respective demagnetizing factors in the directions indicated. The energy density due to the demagnetizing field is f = | o(M x H)' d9 = - 2toM Nz (15) The reversed nucleus is taken to be a rotational spheroid with semiminor axis c, semi-major axis b, with b > > c. The demagnetizing factor for this shape is: = 1 [ln(2k) - 1] (16) k10

Where k = b/c > 500, thus (16) is a very good approximation for Nz~ The total energy for the nucleus may be approximated now by: U = 8i c2boHM - awl2 cb - 8t oM2 c [ln(2k) - 1] (17) 3 3 b The critical length of the nucleus with respect to a given radius is found by differentiating the total energy with respect to b, and then placing this result in (14) using the equal sign. From (14): au= 8it oHoMc2 (18) ab 3 The result of this operation after some rearrangement is: ln(2k) - 2 c - 8_oM(H - Ho) (19) k2 8t8~oM2 If H exceeds Ho for a given radius c, the nucleus will continue to growo When c = a, the radius of the cylinder, the surface of the nucleus on merging with the boundary splits into two sections of length b. Each section propagates in opposite directions at a velocity determined by (H-Ho). The shape of the propagating domain wall is now considered to be one symmetric half of a prolate spheroid. The applied field less the critical field must now supply the eddy current and damping losses. The energy integral is taken over the volume containing the region of the propagating wall, since the magnetization in the remainder of the cylinder is stationary. The true shape of the wall in motion has a greater slope at the cylinder boundary and less slope along the lateral surface than is given by the above assumed shape 11

Nucleation, in reality, may occur anywhere in the volume of the media inside the pulse coil. Two or more nuclei may form, not necessarily simultaneously, grow, merge, then propagate along the wireo Since the pulse coil has a short rise time, eddy currents induced by the pulse field delay the penetration of the field into the wire, thus it is more probable that a nucleus be formed near the surface rather than at the center of the wire. The nucleus would spread over the wire surface and then penetrate to the center, where it would then break into two sections moving in opposite directions with the tip at the center trailingo The general shape of the wall is approximately the same whether the tip leads or trails in motion. Either type of motion is equally probable. INDUCED VOLTAGE AND EDDY CURRENT DAMPING The solution for the eddy current loss would proceed as outlined in Equation (10) through (12), however the circular symmetry of the domain wall allows a more direct approach. Since the length to radius ratio of the wall is greater than 500, the wall is taken to have a conic shape as shown in Figure 20 The equation of the domain wall is now simply x = k z, k = b/a and the coordinates are moving with the wall at velocity v. The voltage induced in any path circularly concentric with the wire center is: 4xpoMzv V = -4t for r > x k2 V = 0 for r < x (20) 12

Equation (20) gives the voltage per turn induced in any coil wound about the wire The voltage gradient is 9 directed, and constant along any circular path of radius r, from (20) 2iioMvz E;2 r >x (21) k2 r The eddy current power loss is found by integrating aEn through all the volume of conducting media outside the domain wall and over the length 0 < z < b. Inside the wire Pe = 2c E rdrdz (22) o0 X which reduces to i A ioMv o3 Pe- 2) 18 (23) The power loss ina surrounding conducting tube of inside radius al, outer radius a2, and conductivity o2 is from (22) changing the limits on x from al to a2: P Lov)2b3 C n ) (24) The power input of the applied field is: P = 2oMHeta2v (25) 13

Equating the sum of (23) and (24) to (25) gives, after simplification: v = 9bHe (26) 4pMa2 a[l + 3 a i (2n ) or the mobility due to eddy currents alone is: 9b Re =b (27) 440Ma2a[l + 5 In ( ] Where the second term in the denominator of (27) gives the contribution of an eddy current damping tube to decrease the mobility. Assuming the wall to be one half of a rotational ellipsoid, the equation for the voltage induced in a circular path is unchanged, however instead of Equation (27), the mobility becomes Re — 9b (28) [IoMa a[l +.52 - in (j From a comparison of (27) and (28) it is evident that the ellipsoid wall results in a mobility of only 17 percent of that for the conic wall,. The exact mobility cannot be determined until a more proper equation for the wall shape is found. The formal significance of these equations, since they are dimensionally correct, is the manner in which conductivity, magnetization, and wire size affect mobility. The eddy current mobility may be expressed more generally; Re 2 b 4 a (29) 14,3 G2 1_ 40lMa af[l + l In -)] 14

where: f =.111, conic wall f =.64, half rotational ellipsoid wall WALL DAMPING DUE TO RELAXATION LOSS The total loss due to spin relaxation processes cannot be given with any accuracy until a better model of the shape of moving wall is developed. This loss is dependent on knowledge of the spin angular variation through the wall which is unknown in detail. However the power loss per unit of wall area should be of the same form as for a plane wall, thus the total loss will be given to a first order only in terms of a dimensionless constant proportional to the cross sectional area of the wire, inversely proportional to the area of the domain wall, let this constant be q, then in terms of Equation (49) Rr = q5Y ) (30) This at least shows the factors which control relaxation damping mobility. GENERAL EQUATION FOR WALL VELOCITY From the definition of the wall mobility, as derived from the force equation, see appendix 2,uoM v = (H - Ho) (31) 15

where H-Ho is the field which supplies the loss as the wall moves. The total mobility term, considering eddy current, and relaxation is determined then from the reciprocal relation; since Br and Be are additive, or R- = Re-1 + Rr- (32) The final wall equation becomes v = R[H - Ho] (33) where Re and Rr are given by Equations (29) and (30) respectively~ As an example of the relative contribution of eddy current and relaxation damping to the mobility, data were taken from a run on a o72 Ni - o28 Fe wire o0114 cm diameter, see Figure 8o The experimental data are as follows: T = 58 7 kg/mm2 R = 19,500 cm/sec-oe or R = 2o45 mks b = 4085 cm k = 850 Using the data given in Table IV and Equation (29), the eddy current mobility is: Re = 9o37 mks, conic wall Re = lo62 mks, ellipse wall 16

From (32), the relaxation mobilities are: Rr = 3 32 mks, conic wall Rr = negative value, ellipse wall The conic wall gives a realizable relaxation term; that for the elliptic wall is not realizable, hence the wall shape is probably intermediate between these two extremeso The anisotropy constant due to applied tension is 2.58 x 105 ergs/cm3 or at least 100 times the crystalline anisotropy energy. The wall which parameter b = 2.8 x 10-6 cm or several hundred lattice constants. FACTORS CONTROLLING WALL MOTION The material parameters which determine the wall velocity are evident from an inspection of the mobility Equations (29) and (30); however, since these factors add reciprocally to give the net inverse mobility, it is more convenient to look at the gross damping factor B, from Equation (31): B oma (2Gt0M)2 2a2 of Kl 2 qb (1 + c) b af I 1 where: f = o111, conic wall f =.64, elliptic wall All material parameters in (34) increase with a decrease of temperature. The axial wall length increases with decrease in temperature as seen from Equation (19)o The combined effect is to increase B, hence decrease the 17

mobility by lowering the temperature. The critical field Ho also increases as temperature drops. The mobility in a nickel ferrite crystal at room temperature is 150 times the value at liquid nitrogen,17 the change in the range of room temperature is less than 10 percent. The mobility increases from 42 cm/sec-oe at 100C to 58 at 950C in.60 Ni -.40 Fe permalloy, see Tables I and III. Since both mobility and critical field are temperature sensitive, temperature control is a requirement for constant delay time. The field difference for low field wall motion is normally less than an oersted, since stray fields may easily exceed this value, shielding is necessary. The final important requirement for uniform wall motion is a high uniaxial anisotropy in the direction of wall motion. If the magnetization is not aligned, wall shape and velocity have high local variations, thus time delays are not reproducible and pulse voltages are variableo Uniaxial anisotropy is brought about in three ways: (a) High percent reduction in cross section by drawing or rolling. The residual stress due to plastic elongation in the roll direction produces a net uniaxial anisotropy in this direction with positive magnetostriction, transverse if the magnetostriction is negative. The alignment varies through the cross section, being very poor near the lateral surfaceso 2 (b) Alignment of the easy axis of the crystallites in the axial direction by anneal in a magnetic fieldo5 Uniaxial alignment of magnetization is low in soft cubic magnetic materialso 18

(c) Tension in the axial direction for materials with positive magnetostrictiono Compression or combined tension and torsion with negative magnetostriction. Optimum stress is just below the elastic limit22 No single method produces an ideal alignment. The combination of the above three techniques in the order listed produces the best condition for uniform wall motion, however the mobility is extremely high -- up to 78,000 cm/sec-oe. A permalloy baked for about an hour in hydrogen, quenched in air, then stressed to just below the elastic limit yielded a sample with uniform wall motion, and the mobility was reduced to 14,000 cm/sec-oe. The addition of an eddy current damping tube to this sample further reduced the mobility to 7,500 cm/sec-oe. In comparison a cold drawn sample of the same material, when placed in the same eddy current tube, resulted in a reduction in mobility from 13,500 to 3,500 cm/sec-oe, however the wall motion is erratic. The above data were taken with a.60 Ni -.40 Fe wire,.005 cm radius. The following discussion refers to data obtained during the period of this contract. Figures 4 through 10 show results obtained on.72 Ni -.28 Fe wire using the experimental setup shown in Figure 3. The propagating field was measured with a d-c digital voltmeter to within ~53 percent. The pulse and saturation fields were measured on an oscilloscope to within ~10 percent. All fields have been corrected for the ambient field which averaged from.45 to.50 oersted along the solenoid axis. The wire was obtained from the Driver-Harris Company, Harrison, No J. After cold reduction to the required size, the wire is baked in a hydrogen 19

atmosphere, then air quenched to room temperature. The samples were unfortunately subjected to sharp bends on shipment, which resulted in erratic wall motion at low tension. After a slight plastic set wall motion became very uniform over the linear range of propagation. A comparison of Figures 5 and 6 show the effect of added eddy current damping on wall mobility. The.020 diameter wire was placed in a.320 by.158 cm diameter brass tube which decreased the mobility by over 50 percent for the same tension. These two graphs are not a good means of comparing the change in driving field with the addition of the damping tubeo The data for Figure 5 were taken subsequent to the removal of the conducting tube, after which the wire was given an extreme plastic set, hence the large increase in Ho. The effect of tension on mobility and Ho for an initially unstrained sample is illustrated in Figure 9a. Both mobility and Ho reach a minimum just below the elastic limit. After plastic set Ho increases but R remains less than the initial value due to the residual strain. In Figures 7 and 10 are shown the hysteresis in mobility and Ho with change in tension after plastic set. These curves are reproducible in samples subjected to the same degree of plastic stretch, A comparison of Figures 5 and 8, or 7 and 10, show the increase in wall mobility with decrease in wire radius. The ratio of wire areas is 3518, the ratio of mobilities at 50 kg/mm2 tension is 0.80, which compares favorably with the respective ratios predicted from Equation (34). Figure 11 shows the motion of walls caused by reverse domains nucleated 20

at the ends of a wire extending beyond the drive field solenoido End nucleation may be minimized by termination within the region of uniform drive field~ The termination should be designed so as not to increase the tension stress over that in the wire proper, or introduce torsional stress. Tension induced anisotropy did not give complete alignment of the magnetizationo The B-H loops taken for the wire under tension are square, with a maximum magnetization of 7,200 gauss at all tensions above 2.1 kg/mm20 The maximum drive field was 6 oersted, the coercive field was 2.8 oe. The coercive force shows little or no relation to Ho, the latter being quite sensitive to strain. The saturation magnetization at room temperature is 12,000 gauss, or only 60q of the magnetization is axially reversed by the wall motion. Calculations of eddy current damping and induced voltage in sense coils should use the switched magnetization AM rather than magnetization M. Two permalloy tapes were tested but no uniform axial wall motion was found. The tapes were 1/8 mil by 1/16 inch in cross section, and were cold worked with no subsequent annealo The.50 Ni -.50 Fe tape under high pulsed field exhibited random nucleation and wall motion across the tape.24 The 4-79 Mo permalloy tape broke before sufficient tension was applied. The latter material has low magnetostriction and conductivity, neither conducive to uniform axial wall motiono Both tapes had extremely rough edges due to slitting. The edge roughness probably caused local regions with a closure domain structure, which was not removed by the saturation field. These regions are nucleation points 21

for walls which move across the tape width when the drive field is applied, A non-uniform strain distribution was quite noticeable across the tape width which.was probably due to residual strain induced during cold rolling. This would be removed by heat treatment. Non-uniformity of axial wall motion appears to exist even in annealed tapes as evidenced by the mobility data in Table II for.48 Ni - Fe permalloy. For the same cross sectional area the wall mobility in a thin tape will be considerably higher than in the wire since the tape geometry lowers the eddy current contribution to the damping. DESIGN OF A MAGNETIC TIME DELAY SYSTEM In summary, the following criteria should be considered: (a) The magnetic media should be in a long filament geometry, probably circular in cross section rather than rectangular for the reasons noted above. The long, small cross section configuration gives a minimum demagnetizing factor in the preferred axial direction of magnetization, (b) The material, since it is polycrystalline, must have high induced axial anisotropy. Crystal anisotropy should be low and magnetostriction high, however the two are uniquely related, so a compromise is necessary. The.48 to.79 percent nickel content permalloys show a decrease in both quantities as the nickel content is increased. Probably a composition near 60o% nickel is best since the magnetostriction is nearly isotropic, and the crystal anisotropy very low in the ordered state. 22

(c) The wire should be stressed to a point just below the elastic limit, and the structure designed to maintain constant stress. (d) Reverse saturation and driving fields must be very uniform along the axis of propagation~ (e) The magnetic wire must be well shielded from stray fields. (f) A conductive eddy current envelope may be used to lower wall mobility, (g) Some means of temperature control is required to maintain constant mobility. For a one shot time delay d-c supplies may be used. The wire must be initially saturated opposite to the direction of wall motion. To initiate the time delay, the saturation current circuit is opened and simultaneously the reverse domain coil circuit and drive field circuit are energized. Repetitive time delays can only be generated by pulsing the saturation solenoid and nucleation coil. The maximum time delay is small. With reference to Figure 6, the lowest mobility is 4,200 cm/sec-oe, and the lowest velocity is slightly below 2,000 cm/sec. This corresponds to a delay of 50 milliseconds per meter. Long delays can only be realized by use of extremely long lengths of wireo Long time delays in a small package seem improbable in the light of the above restrictions, The following structures are suggested for a magnetic time delay device: (a) The twister geometry appears favorable subject to the requirement of precisely controlled tension along the magnetic wire, This structure 23

simplifies the problem of obtaining a uniform driving field along the wire. A model has been built here, but not tested. The scale is large: a 1/8 inch center conductor wrapped with a helix of.020 permalloy wire. The objective here is to determine whether the stressed wire relaxes with time, if so, to what extent the wall velocity is changed. (b) The magnetic wire is wound under tension helixally over the surface of a thin walled rigid plastic tubeo The saturation and drive field windings are wound as toroidal coils over the magnet wire. For a one shot delay the current sources may be packed inside the tube. (c) An intriguing configuration should be mentionedo A thin uniform layer of permalloy, having negative magnetostriction is plated or bonded on 25 a conducting wire. Permalloy with a nickel content exceeding 85 percent is suitable. When the wire is placed under tension the easy direction of magnetization is circumferential. The saturation and driving fields are produced by currents flowing through the conducting core. The nucleation coil is wound over the surface of the magnetic layer in a skewed manner so as to create a local field reversing the magnetization in the circumferential direction. The circular wall now moves down the structure at a velocity determined by the driving field. This is an ideal geometry since the demagnetizing field in the easy direction is zero. The principal drawback appears to be obtaining uniform properties in the magnetic layer. (d) Another device based on a different principal may prove useful. This is the thin magnetic film shift register developed by Broadbent and 26 McLung. A strip of thin film or tape is placed between axially displaced 24

pump windings. These windings are pulsed alternately to form a pump field which progressively travels along the magnetic strip. The pulse frequency controls the field velocity. The strip is first saturated in the direction opposite to pump field motion. The saturation field is removed and a reversed domain nucleated near one end of the strip. The pump field now forces the reversed domain to progressively jump along the strip. The length of time delay is only limited by the gradual degeneration in shape of the reversed domain as it moves along the tape. Several other structures may be suggested, however the requirements of constant tension and uniform drive field immediately narrows the field to the four described. Structures (b) and (c) should be the basis of further development work necessary to realize a compact magnetic time delay deviceo 25

APPENDIX THE 180~ MAGNETIC DOMAIN WALL Expressing Equation (2) through (6) for the free energy density terms in spherical coordinates: fex = A[(V@)2 + sin2 G (V])2] (35) fan K1 [sin2 20 + sin4 ~ sin2 20] + Km sin2 m (36) fan - ~fh = - oM ~ (H - Rd) (37) The angle 9 gives the orientation of M with respect to a z axis, and M sin 9 subtends the angle 0 with the x axis in a right handed xyz coordinate systemEquation (36) gives the sum of magnetocrystalline and strain induced anisotropy energy, where Km = 23 T, and P is the angle between the applied tension stress and the magnetization. With no field applied the magnetization exists in a state of equilibrium which may be found by minimizing the total free energy of the ferromagnetic body. From experimental observation it is found that the static state of magnetization exists as uniformly magnetized domains oriented parallel to the direction of least anisotropy energy. Near crystal surfaces or nonmagnetic inclusions, closure domains form to minimize the self magnetostatic energy, hence the gross domain configuration is determined by the size and shape of the magnetic body and size, shape, and spatial orientation of the 26

crystals which compose the body. The regions between adjacent uniformly magnetized domains -- the domain walls, are relatively thin in comparison to the thickness of the bounding domains. The walls seek an orientation to minimize the anisotropy, exchange and magnetostatic self energy due to the spatial variation of magnetization separating the uniformly magnetized regions. An exact minimum energy configuration is not theoretically determinable by simply assuming domain walls and drawing in the proper closure domains, but at least this process leads to one of several realizable configurations. The possible wall orientations which may exist in a single crystal is determined by the anisotropy energy,27 only the simple 180~ plane wall is considered here. Consider an infinite wall in the x-y plane, with a tensional stress applied in the x direction. The adjacent domains are aligned with M = +Mx, z =- o, M =-M x) z = + coo The applied field is zero and only a 0 variation is allowed, hence the demagnetizing field normal to the wall is zero. The free energy per unit area in the volume of a cylinder normal to the wall is: 00 aw = fdz (38) where the energy density f is given by the sum of Equations (35) and (36)~ The state of minimum energy is found by applying the variational operator 1d 5f j = o (39) 27

to the total free energy density, subject to the condition that the variation is zero at z = + oo The result is: 2A - fan = 0(40) az2 as which states that the torques exerted by the anisotropy and exchange fields on the magnetization are equal and opposite through the wall. Multiplying by d0 and integrating gives: X = fan 1/ (41) 6z lA Solving (41) for the wall thickness: z = lIn tan (42) z = n K112 K 2 j1/ + K+ + sec2 A 1/2 where 6 = - is the wall thickness parameter. Mathematically the \Ki+KmJ wall has an infinite thickness, however most of the angular variation takes place in a width of 10 6 normal to the wall. The wall energy per unit area is found by inserting (41) and (42) in Equation (38) and integrating: 1/2 1 Km 1/2 aw = 2[A(K1 + Km)] 1+ K1 tanh + K- (43) iw YnKm 172 K1 ) Let an external field be applied in the x direction parallel to the wall in an x-y plane. If the field is much less than the anisotropy field, 2K/M, the magnetization within the domains remain respectively parallel and 28

antiparallel to the field. In the intervening wall, however, the magnetization is at an angle 0 with respect to the field and experience a torque which causes it to start a precession about the applied field. There results a component of magnetization normal to the wall surfaceo Since V B = 0 28 always, a demagnetizing field He is formed perpendicular to the wall which is much greater than the applied field. The magnetization within the wall, thus perturbed by the low applied, field now precesses about the field He at an angular velocity yHeo The wall region moves parallel to its normal in the direction of the domain aligned antiparallel to the field, The spin angular distribution within the moving wall is essentially the same as within the wall at rest, hence 0= f(z - vt) and _ = - v_ = yHe at az or He = ) (44) The energy density within the wall associated with He is 2 which is proportional to the square of the wall velocity. This term is a measure of the additional energy stored in the inertia of the spin system of the moving wall volume, and is conveniently expressed as the kinetic energy of 29 an equivalent wall mass9 per unit area: mY2 ~H dv. 2 2 2 29

Integration over a cylinder of unit area and infinite length normal to the wall, using Equation (44) gives: m 2= 2 do (45) Substitution of Equation (41) into (45) gives the wall mass: ~~~~Km. 1/2+ m = 2L t-1 + i 12 ta1h- 1 (+ 46) K K For either K1 = 0 or Km = 0, Equation (46) simplifies to: m = - (46a) y76 The power loss per unit of wall area due to relaxation damping is Bv equivalent to a gross viscous damping loss B- From Equation (9), making use of the fact that He is much greater than the applied field, the power loss per unit volume in the wall is: ~H =M = sin (47) o? +1 Integrating over a cylinder of unit area, and infinite length normal to the wall gives the power input per unit area necessary to supply the losso From Equation (47) and (44): 2oMrv = v sin2 ( ) (48) Substitution of Equation (41) into (48) after integration and rearrangement gives the gross viscous damping factor: 50

-7(c0+l)5 2 K 4 12 and. for Km ~ K1 to: 2OQM 4 B o(49b) 7(c?+1)6 3 General Equation of Wall Motion: For a plane 180~ domain wall with the field applied parallel to the wall, the differential equation of motion is: Mz + Bz + Osz = 2ktoMH (50) The (az) term is the force which restrains the wall to the initial static equilibrium position0 The term is assumed linear and lossless for motion of the wall over short distances~ For wall motion over a large distance the term varies with the wall position and represents an irreversible loss, due to lattice inperfections and non-magnetic incliddionsO The loss involved in this term has not been placed on a satisfactory quantitative basis and is not accounted for here, (See Reference 3a page 549 for additional detail ) With a constant field applied, since the wall mass is small, the wall almost immediately reaches a constant velocity~ The equation of motion may be more conveniently written as: 31

2toMH cz V =- > -> B B or finally, for low fields v = R[H - H] (51) where R = ~M (52) B Ho = < > (53) B R is the wall mobility, which is a constant for applied fields much less than the effective anisotropy field 2K/M. Ho is an average measure of the domain wall stiffness a over the distance it moves. Both terms, R and Ho, show minor variation as noted by repeated measurement of wall velocity vso applied field under identical conditions of wall propagation. 32

TABLE I DATA ON DOMAIN WALL MOTION IN PERMALLOYS AT LOW FIELD Tension Velocity Range H x 10-3 H Diameter Temp R Composition State 213og/m e 3._-'kg/mm2 cm/sec x 10-3 cm/sec-oe oe cm Ref..14 Ni - Fe Cold drawn 38 2.5 - 10 25 5.4.038 20 - 30 1 46 2,5 - 20 28 5.0 62 350 - 15 32 4.6 92 35.0 - 15 20 3.4.10 Ni - Fe Cold drawn 30 3-0 - 9.0,19 3-7.0334 20 - 30 9 42 350 - 12 19 3.3 o50 Ni -.50 Fe Annealed 4.8 135.0 - 24 47 0.01 20 - 30 10 8,8 8.0 - 24 42.09 12.8 5.0 - 20 42.17.60 Ni -.40 Fe Annealed. 4.8 355 - 10 28.19.01 20 - 30 10 8,8 355 - 10 31.23 12,8 535 - 10 30.31.60 Ni -.40 Fe Annealed. - 14.4.11.005 - 180 10 42,2.05 10 580o.05 95 ~60 Ni -.40 Fe Cold drawn 30 80 - 12 -.01 20 30 10.60 Ni - o40 Fe Cold drawn - 4,2 3536.o005 180 10 11o,9 244 10 16.2 2.22 95

TABLE II DATA ON DOMAIN WALL MOTION IN.48 Ni - Fe TAPES (Ref, 12) T W Ho b R x 10-3 State State mil mil oe cm cm/sec-oe Cold rolled 1 100 4.1 20 + 3 17 ~ 1 Cold rolled 1 40 4.5 9 + 2 22 + 1 Cold rolled 1/8 100 6,1 7 + 2 9 ~ 3 Annealed 1 100 03.5 50 ~ 20 140 + 20 Annealed 1 40 0.48 20 + 3 170 + 20 Annealed 1/8 100 0.31 16 ~ 5 120 + 30 54

TABLE III DATA ON DOMAIN WALL MOTION IN SINGLE CRYSTALS Ho R Velocity Range Tempo Composition e cm secToe mWall Re'............oe. cm/sec-oe cm/sec, K......04 Si - Fe 3.8 - 85 (100) 15 4,7 - 180 5.4 - 290 Fe304.06 1,900 100 - 850 - (110) 16 (NiO) 91(FeO) 09Fe203 1.5 20,000 (110) 6 (NiO)75 (FeO) 25Fe20.075 26,160 0 - 6,000 201 (110) 17 Mn 4Fe1. 604 1.4 94 0 - 500 77 (112) 19.015 805 160 - 370 77 (110) 5o4p x 10l Fe Whisker 0+ 2,500 100,000 - 200,000 77 (100) 20 0+ 7000, 0000 100,000 - 400oo00 200 0+ 8,000 100,000 - 500,000 298

TABLE IV PHYSICAL DATA FOR.72 Ni -.28 Fe PERMALLOY Values are for the temperature range of 25 - 30~C Magnetostriction constant: (assumed to be isotropic) X' 15 x 10-6 Crystal anisotropy energy: ergs/cm3 KI = 2 x 103 quenched K1 = -7 x 103 furnace cooled K1 = -20 x 103 slowly cooled K2/K1 < o06 Exchange constant: ergs/cm A 2 x 10-6 reported values range from 1.6 x 10-6 to 3.3 x 106 29 Relaxation damping constant: dimensionless 0.3 < a < 9 dependent on frequency temperature, and technique of measurement. Magnetization: gauss saturation:!tMs = 12,000 from B-H curves 4tM = 7,200 Conductivity: mhos/cm permalloy, a = 5.02 x 104 brass tube, a = 14.9 x 10 56

I a~ z M & z H Fig. 1. Nucleation of a reversed domain in a uniformly magnetized cylinder. x H v F _ 1- -— Z -z I b Fig. 2. Domain wall shape used for calculation of eddy current damping. 37

SYNC <1 1<i _ -f',SQUARE WAVE PULSE GENERATOR -o2 COIL Li Hs - 3 SEARCH Lp__ *-2, J II COILS HC lo WITC Fig. 1. Sysm fr HS DIGITAL 8VOLTMETER PROPAGATION SOLENOID -- SATURATION SOLENOID 5700 TURNS ~J34 5700 TURNS *34 TENSION ~ —-lo -- 0 ol TO NUCLEATION ~.. --- o — 2 PULSE COIL SQUARE WAVE Hp GENERATOR FROM HS ~ —— ~ SEARCH 3 0 —' ~Hpc__ - 0~ —-- ~COILS 4 o ELECTRONIC CRO AMPLIFIER CRO SWITCH Fig. 3. System for measurement of wall velocity. 38

Scale:.23 mv/cm x 1 ms/cm Scale:.5 mv/cm x 1 ms/cm'IO (a).0203 am diam. in a brass tbe (b).0114 cm diam. Fig. 4. Output voltage of search coils wound on the permalloy wire.

! - I I Annealed.72 Ni.28 Fe 12 After Plastic Set at 71. kg/mm2.0203 cm Diameter T: kg/mm2 R x 103 cm/sec-oe. 0 30.4 17.9 // I0 A 36.8 15.2 oD 45.8 12.5 V 55. 1 11.4 + + 67.5 10.8 - 8-V 55.1 12.0/ 7 So U 45.8 14.0 cn A 36.8 15.3 E 30.4 20.4.6 0 /.. 4' /./ /4 - 2 / /o o/,/ / / / / 1.2 1.4 1.6 1.8 2.0 2.2 2.4 2.6 2.8 H IN OERSTEDS Fig. 5. Axial wall velocity vs. applied field for a.020 cm permalloy wire under tension.

8..72 Ni-.28 Fe Wire in a Brass Tube After Plastic Set at 61.3 kg/mm2 Diameter in cms: Wire =.0203 Tube=.320 0.0. by.158 LD. 1I to 6. T: kg/mm Rx 10: cm/sec-oe. Io 6. —o 0 27.8 7.81 x A 42.8 4.96 o D 49.0 4.53 -l V 55. 1 4.22 4.-0 61.3 4.18 0U ITension Increased z 2. ~/ / // / oL-/ /V _ 4,i 1.8 2___/.6 1.0 1.4 1.8 2.2 2.6 H IN OERSTEDS Fig. 6. Wall velocity vs. applied field for a.020 cm permalloy wire in an eddy current damping tube. 41

2. 5 2.5 ~Annealed.72 Ni-.28 Fe After Plastic Set At 71 kg/mm2.0203 cm Diameter 2.3 0 Ho vsT 0 R vsT 2.1 \ 20 Cn \ \W \\ C C,) 03 b0 ~~~~~~~C3l) w~~~~~~~~~~~~~~~~~~~~~ peUmalloy wire 1 9 —--—. —5- ^ --------- 184m42 ~ \:: a 0 0 z w wie 1.3 -- 12 TENSION: kg/mm2 Fig. 7. Hysteresis of Ho and. R with change in tension in a.020 cm permalloy wire. 42

12 __ __ _ Annealed.72 Ni-.28 Fe Wire;.0114 cm Diameter, After Plastic Set at 58.7 kg/mm2 10- T:kg/mm2 RxlO-3 m / cm/sec-oe 0 19.6 24.4 to 0 39.1 20.2 o0 58.7 19.5 X 8_ 39.1 20.8 W 0 19.6 24.4 0) E A ~F / A w v v 0 _6 wu 4 1.2 1.2 1.4 1.6 1.8 2.0 2.2 2.4 H IN OERSTEDS Fig. 8. Axial wall velocity vs. applied field for a.0114 cm permalloy wire under tension.

72 Ni -.28 Fe 2.2 A\ Annealed 34.0114 cm Diameter -\- R vsT On — 0 — Ho vsT a oj \ \ 5 0 1- \ E o 1.0 22 \ rr15 25 35 45 55 TENSION: kg/mm2 Fig. 9. Hysteresis of Ho and R with change in tension in a.0114 cm annealed permalloy wire. 44

.72 Ni-.28 Fe 1 C.8____ iAfter Plastic Set 1.8 30.0114 cm Diameter )n — ) — R vsT o c3 Ho vsT _, D: X w 6 18 I 1.0 22 E.6...18 15 25 35 45 55 TENSION: kg/mm2 Fig. 10. Hysteresis of Ho and R with change in tension in a.0114 cm permalloy wire after plastic stretch. 45

C) o Motion is u' 0 Irregular Cn l \.72 Ni -.28 Fe I_ J I \I | After Plastic Set Et~_^,-l o~.)~ ~.0114 cm Diameter 10' 1 I ~ / — __\ —__H = 1.26 oe W\ Hp= 0 SU //) \Hs 11.5 < Dir. of Dir. of T 58.7 g/mm a. Motion Motion <[ 5 v = 4,500. L= 4,240 ~ 3 cm/sec cm/sec 0 w I — 0 20 40 60 80 100 DISTANCE ALONG WIRE: cm Fig. 11. Motion of domain walls nucleated at the ends of main solenoid.

REFERENCE 1. Sixtus, K. J. and Tonks, L., "Propagation of Large Barkhausen Discontinuities," Phys. Rev. 37, 930 (1931). "Further Experiments on the Propagation of Large Barkhausen Discontinuities," Physo Rev. 39, 357 (1932a).'Propagation of Large Barkhausen Discontinuities II," Phys. Rev. 42, 419 (1932b), "Propagation of Large Barkhausen Discontinuities III, Effect of Circular Field with Torsion," Phys. Rev. 43, 70 (1933a). "Propagation of Large Barkhausen Discontinuities IV, Regions of Reversed Magnetization," Phys. Rev, 43, 931 (1933b), 2. Gyorgy, E. M., "Flux Reversal in Soft Ferromagnetics," J. Appl, Phys. 31, 110o-s (1960), 3a. Kittel, C., and Gait, Jo Ko, "Ferromagnetic Domain Theory," Solid State Physics 3, 439 (1956), 3b. Bozorth, Ro M., "Ferromagnetism," Do Van Nostrand Co. (1951). 3c. Stewart, K. Ho, "Ferromagnetic Domains," Cambridge University Press (1954)o d-. Bozorth, R. Mo, and Walker, Jo Go, "Magnetic Crystal Anisotropy and Magnetostriction in Iron-Nickel Alloys," Physo Rev. 89. 624 (1953)o 5. Ferguson, Eo To, "Uniaxial Anisotropy Induced in Fe - Ni Alloys by Magnetic Anneal," J. Appl. Phys. 29, 252 (1958)o 6. Gilbert, T. L., and Kelley, J. Mo, "Anomalous Rotational Damping in Ferromagnetic Sheets," Proc. Pitt. Conf. on Mag. and Mag. Mato, AIEE, T-78, 253 (1955)o 7- Menyuk, N., "Magnetic Materials for Digital Computer Components II. Magnetic Characteristics of Ultra-Thin Mo-Permalloy Cores," Jo Appl. PhySo 26, 692 (1955). 8. Haake, Ho, "Uber die Keime und die Ausbreitung der Ummagnetisierung bei grossen Barkhausen-Sprungen," Z. f, Phys. 113, 218 (1939). 9. Reinhart, Ro E., "Large Barkhausen Discontinuities and their Propagation in Ni - Fe Alloys," Part I, Phys. Rev. 45, 420 (1934); Part II, Phys. Rev. 45, 483 (1934)o 10. Dijkstra, L. Jo, and Snoek, Jo L., "On the Propagation of Large Barkhausen Discontinuities in Ni - Fe Alloys," Philips Res. Repo 4, 334 (1949)o 47

11. Ogawa, S., "Uber der grossen Barkhausensprung. Dritter Teil. Anwachsen und Verschwinden des Ummagnetisierungskeimes," Sci Rep* Tohoku Univo 1, 53 (1949). 12. Cole, Ro W,, "Motion of Ferromagnetic Domain Walls in.48 Ni - Fe Tapes," J. Appl. Phys. 27, 1104 (1956). 135 Doring, W., "Uber das Anwachsen der Ummagnetisierungskeime bei grossen Barkhausen - Sprungen," Z. f. Physo 108, 137 (1937). 14. Brown, W. F., Jr., "Micromagnetics, Domains and Resonance," J* Applo Phys. 30, 62-S (1959). 15. Williams, H. J., and Schockley, W., and Kittel, C., "Studies of the Propagation Velocity of a Ferromagnetic Domain Boundary," Phys. Rev, 80, 1090 (1950). 16, Galt, J. K., Andrus, J., and Hopper, H. G., "Motion of Domain Walls in Ferrite Crystals," Rev. Mod. Phys. 25, 93 (1953). 17. Galt, J. K., "Motion of Individual Domain Walls in a Nickel-Iron Ferrite," Bell Sys. Tech. J. 33, 1023 (1954). 18, Gait, J. K., "Motion of a Ferromagnetic Domain Wall in Fe304," Phys. Rev. 85, 664 (1952). 19< Dillon, J. F., and Earl, H. E., "Domain Wall Motion and Ferrimagnetic Resonance in a Manganese Ferrite," J, Appl. Phys. 30, 202 (1959). 20, DeBlois, R. W., "Magnetization Processes: Reversals and Losses," J. Appl. Phys. 29, 459 (1958). 21, Olmen, R, W., and Mitchell, E. M,, "Slow Domain Wall Motion in Homogeneous Vacuum Deposited Iron-Nickel Films," J. Appl. Phys. 20, 258-S (1959). 22, Neurath, P. W., "The Effect of Plastic and Elastic Stresses on the Losses and the Domain Configurations of Grain Oriented 3% Si - Fe," Proc.. Pitt. Conf. on Mag. and Mag. Mat., AIEE, T-78, 92 (1955). 235 Chikazumi, S., "Magnetic Anisotropy Induced by Magnetic Annealing and by Cold Working of Ni3Fe Crystals," J. Appl. Phys. 29, 346 (1958)o 24, Rodbell D. S., and Bean, C. P., "Influence of Pulsed Magnetic Fields on the Reversal of Magnetization in Square Loop Metallic Tapes," Jo Appl. Phys. 26, 1318 (1955)o 25. Long, T. Ro, "Electrodeposited Memory Elements for a Nondestructive Memory," tJ Appl. Phys. 31, 123-S (1960). 48

26, Broadbent, K. D., and McClung, F. Jo, "A Thin Magnetic Film Shift Register," Hughes Res. Lab., Res. Rept. 136 (1960). 27. Lilley, Bo, "Energies and Widths of Domain Boundaries in Ferromagnetics," Phil. Mag. 41, 792 (1950). 28o Becker, R., "Dynamics of the Block Wall and Permeability of High Frequencies," J. Physo Rad. 12, 332 (1951). 29. Doring, W,, "Uber die Tragheit der Wande Zwischen Weisschen Bezirken," 4 f Naturo 3A, 373 (1948). 30~ Weertman, J. R., and Rado, Go To "Low Temperature Spin-Wave Resonance at 3,000 and 4,000 mc in a Permalloy Having Nearly Zero Magnetocrystalline Anisotropy," J. Appl. Phys. 29, 328 (1958)o 49

DISTRIBUTION LIST (One copy unless otherwise noted) Air Force Office of Scientific Wright Air Development Division Research 3 (WWAD) Washington 25, D.C. Wright-Patterson Air Force Base Attn: Solid State Sciences Division Ohio Air Force Office of Scientific Wright Air Development Division Research 2 (Metals and Ceramics Laboratory) Attn: Technical Library (SRGL) Materials Central Wright-Patterson Air Force Base Armed Services Technical Ohio Information Agency (TIPCR) 10 Arlington Hall Station Wright Air Development Division Arlington 12, Virginia (Physics Laboratory) Materials Central Commander 2 Wright-Patterson Air Force Base Army Rocket & Guided Missile Agency Ohio Redstone Arsenal, Alabama Attn: ORDXR-OTL Wright Air Development Division (Materials Information Branch) RAND Corporation 2 Materials Central 1700 Main Street Wright-Patterson Air Force Base Santa Monica, California Ohio Air Research and Development Command Institute of Technology (AU) Library (RDRS) MCLI-LIB, Building 125, Area B Andrews Air Force Base Wright-Patterson Air Force Base Washington 25, D.Co Ohio EOAFRD ARL (Metallurgy) Air Research and Development Command AFRD 47 Cantersteen Wright-Patterson Air Force Base Brussels, Belgium Ohio Hq., USAF (AFDRT) ARL (Physics, Solid State) Washington 25, D.C. AFRD Wright-Patterson Air Force Base ARL (Technical Library) Ohio Building 450 Wright-Patterson Air Force Base AFOSR (SRLTL) Ohio Holloman Air Force Base, New Mexico

DISTRIBUTION LIST (Continued) AFCRL (CRRELA) Director L. Go Hanscom Field U. S. Naval Research Laboratory Bedford, Massachusetts Washington 25, D.C. Attn: Library AFFTC (FTCTL) Edwards Air Force Base, California Department of the Navy Office of Naval Research AEDC (AEOIM) Washington 25, D.C. Arnold Air Force Station Attn: Code 423 Tennessee Attn: Code 421 AFSWC (SWOI) Officer in Charge Kirtland Air Force Base Office of Naval Research New Mexico Navy No. 100 Fleet Post Office Office of the Chief of Research New York, New York and Development Department of the Army Commanding Officer Attn: Scientific Information Naval Radiological Defense Laboratory Washington 25, D.C. San Francisco Naval Shipyard San Francisco 24, California Army Research Office Box CM, Duke Station Dr. D. F. Bleil Durham, North Carolina Associate Technical Director Attn: CRD-AA-IP for Research U. S. Naval Ordnance Laboratory Commanding Officer White Oak, Silver Spring, Maryland Ordnance Materials Research Office Watertown Arsenal National Aeronautics & Space Agency Watertown 72, Massachusetts 1520 H Street, NoW. Attn: PS&C Division Washington 25, D.Co Commanding Officer High Speed Flight Station (NASA) Watertown Arsenal Edwards Air Force Base, California Watertown 72, Massachusetts Attn: Technical Library Attn~ Technical Reports Section Langley Research Center (NASA) Commander Langley Air Force Base, Virginia Signal Corps Engineering Laboratory Attn- Technical Library Fort Monmouth, New Jersey Attn~ SIGFM/EL-RPO U. S. At;omic Energy Commission Library Branch Ames Research Center (NASA) Technical Information Division, ORE Moffett Field, California Po 0. Box E Attn: Technical Library Oak Ridge, Tennessee

DISTRIBUTION LIST (Concluded) Lewis Research Center (NASA) Knolls Atomic Power Laboratory 21000 Brookpark Road P. 0. Box 1072 Cleveland 35, Ohio Schenectady, New York Attn: Technical Library Attn: Document Librarian Wallops Station (NASA) National Science Foundation Wallops Island, Virginia 1901 Constitution Avenue, N.W. Attn: Technical Library Washington 25, D.C. Division of Research National Bureau of Standards Library U. S. Atomic Energy Commission Room 203, Northwest Building Division Office Washington 25, D.C. Washington 25, D.C. Director Major John Radcliffe Office of Technical Services ANP Office Department of Commerce U. S. Atomic Energy Commission Technical Reports Branch Washington 25, D.C. Washington 25, D.C. Oak Ridge National Laboratory Chairman Post Office Box P Canadian Joint Staff (DRB/DSIS) Oak Ridge, Tennessee 2450 Massachusetts Avenue, N.W. Attn: Central Files Washington, D.C. Brookhaven National Laboratory Defense Research Member Upton, Long Island, New York Canadian Joint Staff Attn: Research Laboratory Director of Engineering Research Defense Research Board Argonne National Laboratory Ottawa, Canada 9700 S. Cass Avenue Attn: Mr. H. C. Oatway Argonne, Illinois Attn: Librarian Institute of the.Aeronautical Sciences 2 East 64th Street Document Custodian New York 21, New York Los Alamos Scientific Laboratory Attn: Librarian P. 0. Box 1663 Los Alamos, New Mexico Ames Laboratory Iowa State College P. O. Box 14A, Station A Ames, Iowa

UNIVERSITY OF MICHIGAN 3 901 0369 4411111 3 9015 03695 6244