r-srura- APPLICATION OF A HIGH PRESSURE MBROBOMB T0 HYDROSTATK} PRESSURIZATION 0F MAGNET” E Thesis Io: tho Doom of 81.8. MICHIGAN STATE UNIVERSITY JOHN KELLY ”AH“! 1977 .t."' ABSTRACT APPLICATION OF A HIGH PRESSURE MICROBOMB TO HYDROSTATIC PRESSURIZATION OF MAGNETITE By John Kelly Maher This thesis presents a study of hydrostatic pressure effects on the remanent magnetization of magnetite. Two types of samples were used; a single crystal of magnetite, and a synthetic "rock" composed of a resin matrix and a powder of synthetic magnetite. The samples were initially saturated in a 2 kilogauss field, then pressurized to 4.0 kilobars in an oil-filled non- magnetic chamber. Isothermal remanent magnetization was measured during pressurization by means of a ballistic magnetometer. Pressure was increased and decreased in a continuous cycle during the tests. During the first pressurization, the sample magnetization dropped to about 70% of the initial saturation magnetization with a pressure of 1.5 kilobars, and continued to drOp at a slower rate to about 55% of the initial saturation with a pressure of 3.0 kilobars. Pressure was reduced to 2.0 kilobars with no further change in magnetization, then increased to 4.0 kilobars with another 5% drop in magnetization to about 50% of the initial value. The thesis also presents some alterations made in a high pressure microbomb used by earlier researchers (Car- michael, et. al., 1968). The paper discusses both the development and use of the bomb and the complimentary apparatus required for pressurization and magnetic measurements. The latter includes a ballistic magnetometer specifically designed by R.S. Carmichael for use with this system. ACKNOWLEDGEMENTS I would like to thank my advisor, Dr. Robert 8. Car- michael, for his advise and support during this work. He not only designed and supplied a good deal of the equipment used, but also furnished many of the ideas for this work. Dr. Ted Vinson of the Civil Engineering Department supplied the hydraulic press and other items of support. Dr. Jerry Cowen of the Physics Department loaned a digital voltmeter and integrator for the ballistic magnetometer. I would like to extend special thanks to my colleague Joon Kim, who assisted in much of the experimental work, and frequently put up with my unorthodox working hours. The figures were produced by Yvonne Plaisance and Jake Keeser of Chevron Oil Company. ii APPLICATION OF A HIGH PRESSURE MICROBOMB TO HYDROSTATIC PRESSURIZATION OF MAGNETITE A THESIS SUBMITTED TO THE DEPARTMENT OF GEOLOGY OF MICHIGAN STATE UNIVERSITY IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF MASTER OF SCIENCE By JOHN KELLY MAHER JUNE, 1976 TABLE OF CONTENTS List of Figures I. II. III. IV. VI. Introduction , Theory A. Domain Energies 1. Exchange energy 2. Magnetocrystalline anisotropy energy 3. Domain wall energy 4. Magnetostrictive anisotropy energy 5. Magnetostatic energy 6. Field energy B. Grain Size Effects Equipment and Experimental Techniques A. Previous work B. Pressure microbomb C. Operation D. Ballistic Magnetometer Sample Specifications and Experimental Results A. General prOperties of magnetite B. Samples C. Magnetic results Interpretation of Results Conclusion Appendix A - Bomb Calibration Data Bibliography iv 14 14 14 19 22 29 29 30 38 40 42 43 I-l I-2 I-3 II-l 11-2 11-3 11-4 II-S 11-6 11-7 III-1 III-2 III-3 III-4 LIST OF FIGURES Possible domain structures Fringing field at a plane surface Variation in rotation with distance from domain wall center Non-magnetic pressure chamber and ballistic magnetometer Piston assemblies A typical calibration of pressure microbomb us1ng NH4F Hydraulic press force vs internal confining pressure inside bomb Sense coil field Mutual cancellation of solenoid and lab fields Schematic diagram of equipment set-up Magnetization vs pressure - sample 7 Magnetization vs pressure - sample PH-3 Magnetization vs pressure - sample PH-3 Magnetization vs pressure - sample PH-3 12 12 15 18 21 23 24 26 27 34 35 36 37 I. INTRODUCTION To date, laboratory research on magnetic properties of earth materials under simulated crustal conditions has cen- tered on high temperature experiments and tests involving application of uniaxial pressure. Little work has been done concerning properties of minerals and rocks under realistic high hydrostatic pressure. The latter can be an important factor in the behavior of physical parameters of rocks and minerals. For example, Stesky and Brace (1973) have shown that the electrical conductivity of oceanic basalts varies with increased confining pressure, and Carmichael (1969) has shown a variation of coercive force in magnetite with hydro- static pressure. Breiner and Kovach (1966) used pressure- induced changes in magnetization of rock to try to forecast seismic events along the San Andreas fault. Other recent studies of applications of geopiezomagnetism to earthquake prediction include those by Golovkov (1969), Rikitake (1968) andJohnston (1975). Davis and Stacey (1971) noted that a local magnetic anomally was produced in the area of a dam, as a result of crustal loading when the reservoir was filled. Since most crustal rocks are, or have been, subjected to confining (hydrostatic) pressures of at least several kilobars, it is important in the study of geopiezomagnetism to understand the effects of confining pressure on magnetic materials. Both transient and permanent effects of pressure need to be studied to understand the true nature of magnetic remanence in rocks, and how it changes in the lithosphere. Specifically, this paper is intended to study the nature of the variation of the magnetization of magnetite within the pressure range of one atmosphere to five kilobars. Isothermal remanent magnetization was monitored as hydro- static pressure on the magnetite samples was increased from one atmosphere to about three kilobars, then decreased to two kilobars, increased again to four kilobars and finally released to about one kilobar (the pressure held within the pressure bomb by the frictional force of the packing mater- ial). A significant drop was noticed in the magnetization with the first pressurization to two kilobars, J dropping to about 70% of the initial saturation. Little change in J was noticed beyond 2.0 kilobars. The samples used were a magnetite crystal and a synthe- tic magnetite powder. The powder was mixed with a resin matrix to simulate a rock with a high magnetite content. A fundamental part of this type of study is the develop- ment of technology capable of both simulating a high pres- sure environment and allowing for convenient measurement of the magnetic properties of the sample. This work is intended to contribute something both to the technology of current high pressure experimentation and to the understanding of rock magnetism. II. THEORY Domain theory describes the magnetization of individual mineral grains; the magnetization of such grains in rock governs the magnetic behavior of the rock as a whole and is therefore of prime importance in the study of rock magnetism. The limit of how well the magnetic properties of minerals are understood is the limit to how well paleomagnetism may be understood and relied upon. Grain size is a prime factor in the magnetization of mineral grains and as such is an important factor in how long a rock will hold its natural remanent magnetization, and to what extent changing conditions affect that remanence. The effect of grain size may be better understood by viewing the energies involved in rock magnetization. Both the energies and sizes of grains are interrelated with the structure of domains within ferromagnetic mineral grains. A number of good reviews of the subject may be found in the literature (Neel, 1955; Stacey, 1963; Kittel, 1949, Chikazumi, 1964; Stacey and Banerjee, 1974). DOMAIN ENERGIES The purpose of domains in ferromagnetic minerals is to minimize the total magnetic energy of the particular speci- men. Since ferrites possess spontaneous magnetization below the Curie temperature, magnetic fields will be set up within the material, and these will extend beyond the sample sur— face. This increases the magnetic energy of the specimen, causing it to be in an unstable state-~hence the formation of domains. The total magnetic energy of a ferrite grain may be divided into six components, which will be briefly described. 1. Exchange Energy, E This is produced by the interaction x. of electron spins of adjacent atoms, and was derived initial- ly by Heisenberg. This energy can be minimized when the magnetic moments of adjacent spins are parallel. If a domain wall between two antiparallel domains is considered to be of infinitessimal thickness, then the angle between the spin axes will be 180°, with a maximum exchange energy. But by making the domain wall relatively thick, with many contained spin axes, then the angle between any two spin axes can be reduced, thus reducing the exchange energy density. 2. Magnetocrystalline Anisotropy Energy, Ek- Sometimes called the anisotropy energy or crystal magnetic anisotropy energy (Chikazumi, 1964), this energy is a measure of the preference of a ferromagnetic crystal for spontaneous magne- tization along particular crystallographic axes (Kittel, 1949; Chikazumi, 1964). The direction of preferred magnetization is called the direction of "easy” magnetization; the most difficult directions are called the hard directions. The excess energy required to magnetize a grain in the hard di- rection, as opposed to the easy direction, is the anisotropy energy, Ek (Kittel, 1949). '5 Expressions for Ek, in terms of the direction cosines of the internal magnetization, have been worked out in most discussions (Chikazumi, 1964; Kittel, 1949; Stacey, 1963). For a cubic crystal system, two anisotropy constants are used, K1 and K2; the value of Ek for a cubic crystal, such as given below (01, 62, 03 are the direction cosines of the magnetization vector with respect to the three cube edges): 2 2 2 2 2 2 + “2 “3 I 0‘3 0‘1) I K2 0‘1 “2 “3 Both anisotropy constants are temperature-dependent, and can vary enough to change sign of the net Ek' K1 is generally more important than K2 (Stacey, 1963). In magnetite, the easy axes are the axes, and the anisotrOpy constants are negative, so we are primarily con- cerned with a magnetization inclined at some angle 0 to the axis. To a second-order function in sin 0 (Stacey and Banerjee, 1974) the above equation reduces to K . E.= e +3.9) - as! +212) ., the second term dropping out as the magnetization aligns itself with the easy axis. For pure magnetite at 290°K, K1 = -1.36 x 105 ergs/cc and K2 = -0.44 x 105 ergs/cc, making at -1.80 x 105 E ergs/cc, lowest along <111>. k! A more gradual rotation of spins in domain walls, while decreasing the exchange energy, increases Ek by forcing some spins out of the easy direction of magnetization. However, ‘while EX is affected by this movement of spins within the crystal, the exchange energy itself does not produce any anis- otropic effects. The origin of magnetocrystalline anisotropy is due rather to atomic interactions: "the spin interacts with the orbital motion (of the electrons) by means of the spin—orbit coupling and the orbital motion in turn interacts with the crystal structure by means of the electrostatic fields and overlapping wave functions associated with neigh- boring atoms in the lattice" (Kittel, 1949). 3. Domain Wall Energy, Ew. This is a function of the mini- mized sum of the exchange and anisotropy energies (Ex + Ek). In a uniaxial model the exchange energy per unit area of wall is given by (2) EX = fA(6¢/6X)2dx (Craik and Tebble, 1965) where A is the exchange constant and 4 is the angle of spin rotation as a function of x (see Fig. I-3). In the same model, Ek per unit area is given by = ' 2 (3) ER le Sln (¢)dx where K is the anisotropy coefficient. The energy of the wall using the above values for Ex and Ek is then given by (4) Ew = 2 (AK1)% (Stacey, 1963) this being approximately 1 erg/cmz, while letting A equal the exchange energy for spin S and lattice spacing "a" -- A = J is the exchange integral. 7 4. Magnetostrictive Anisotropy Energy, EA . When a crystal o is stressed, the domains and walls within it are physically moved. This movement causes subsequent changes in the direc~ tions of the domain magnetization vectors, and a change in the net magnetization of the crystal. The energy associated with this change in magnetization is the magnetostrictive anisotropy energy (also referred to as the magnetoelastic energy or the magnetostrictive strain energy). BAG is defined to be zero for an unstrained lattice (Kittel, 1949). It should be kept clear that when the magnetization induces a physical change in the shape of a specimen the effect is termed magnetostriction; when applied stress on the specimen induces a change in magnetization, the effect is termed in- verse magnetostriction, or piezomagnetism. When a ferromagnetic body is placed in a magnetic field, interaction between the atomic magnetic moments causes a change in bond lengths, since E10 is a function of the inter- atomic distance 3. For a cubic lattice, EA may be expressed o in terms of the lattice strain tensor components (eij) and the direction cosines of the domain magnetization (ai) as (Chikazumi, 1964) for the one-dimensional case: (5) Bio = B1<§xx(a - l/3)+ eyy(a - 1/3) + ezz(a - 1/3)) + Bz(exya1a2 + eyzaZa3 + ezxa3al) _ 91L _ where B1 — N(dr ro , B2 - ZNL and N = demagnetization factor, L = specimen length. 8 If a stress 0 is applied to a ferromagnetic material, the energy due to that stress may be expressed by giving the strain tensor components of the magnetoelastic energy in terms of the direction cosines gi, the stress, and the elas- tic constants SijI as, _ p 2 2 2 (6) exx — o(sllg1 + 512(g2 + g3)).... By substituting (6) and (5), we arrive at a new expres- sion for E : Ao _ _ 2 2 2 2 2 2 _ (7) Exo ‘ B10(511 S12)(31g1 I 3282 I 3333 1/3) I B2°544(31a2g1g2 I azasgzgs I a3alg3g1) If, for example, the magnetization of the domain is parallel to the (111) face, Elo may be expressed in terms of 0, the angle between the direction of tension in the crystal and the (111) face: (8) cos(0) = (1/3)2 (g1 + g2 + g3) and 3 (9) BAG = + 3 Alllocosz(0) The same holds for other domain directions, with the substitution of the correct value for A, the saturation magnetization strain (Chikazumi, 1964). The dependence of the magnetostrictive energy on the domain configuration becomes complex in the threeédimensional real case. However, a two-dimensional example computed for 9 the four domain case of Fig. I-1.c (Stacey and Banerjee, 1974) serves to give an order of magnitude value for EA 3 S o in magne- tite of 6 x 10 ergs/cc, compared to -2 x 10 ergs/cc for Bk. 5. Magnetostatic Energy, Em. When free poles exist at the surface of a grain, they give rise to the magnetostatic energy, given by (10) Em = %NJ: (for an ellipsoid) (Stacey, 1963) where JS is the saturation magnetization parallel to the long axis of the ellipsoid and N is the demagnetization factor. The strength of Em is dependent upon the shape of the sample and the number and strength of free surface poles. 6. Field Energy, Eh. A grain may also acquire a field energy when in a magnetic field H, with magnetization l of the sample. The energy per unit volume is given as (11) Eh = HJ cosO (0 = angle between J a H) (Stacey, 1963). This may cause the enlargement of domains oriented in the general direction of H, while opposing domains will be re- duced. GRAIN SIZE EFFECTS The importance of grain size in the study of domains is evidenced by the several properties dependent upon it: such properties include coercive sorce, susceptibility, and transi- tions from the superparamagnetic state to single domain and multidomain states. 10 Consider first a spherical isotropic grain of sufficient size to contain only one domain. The energy of the grain is then given by the magnetostatic energy only (Nagata, 1961): (12) Em = l/zNJ:(4/3)1TR3 , where N = “24/3. As grain size increases, a point is reached where the inclusion of a domain wall will decrease the energy of the grain, by decreasing the magnetostatic energy by a greater amount than is gained by the addition of the wall energy. The new grain energy is given by the sum of BW and Em, the new total energy now being reduced by 50%: 2 (13) E = %(%J:R3)(4n/3)2 + nR a total w where aw is the wall energy per unit area. The theoretical critical radius, R2, may be found by equating the total energies of the single and multidomain states. In this instance, the value for R2 is (Nagata, 1961) _ 2 (14) R2 - 9aw/41TJS Single domain particles range in size from 0.1 microns in diameter to 0.5 microns. From 0.5 to 20.0 microns grains undergo a transition from the single domain state to the multidomain state. Grains in this range are termed pseudo- single domain particles. Grains larger than 20.0 microns zare multidomain only. (These values are for magnitite). Single domain particles have coercive forces that are ggreater than either SPM or MD particles. Neel (1955) states 11 that this occurs because wall displacements, necessary in the movement of a multidomain particle through the hysteresis cycle, require a release of energy, while this is not a fac- tor with single domain grains. SPM particles have smaller values of HC because below a certain volume, 1, the thermal agitation energy kT becomes large with respect to the energy terms dependent upon X and the anisotropy constant E, such that kT exceeds Kv Sin2 (0). This perturbs the precession of the atomic moments, altering O and causing the remanence to decay with time, according to Jr - Joexp(-t/t0), to being the relaxation time (Neel, 1955). POSSIBLE DOMAIN STRUCTURES (After Stacey, 1963) FRINGING FIELD AT A PLANE B (After Stacey, 1963) 13 180 1 I l 150v / ¢ 90 FIGURE l-3 ABOVE IS THE VARIATION IN ROTATION WITH DISTANCE FROM THE DOMAIN WALL CENTER. BELOW IS THE MAGNETIZATION DIRECTION WITHIN THE DOMAIN WALL. (After Craik and TobeOI. III. EQUIPMENT AND EXPERIMENTAL TECHNIQUES A large part of this study involved the development of new high-pressure apparatus and subsequent modifications of this equipment. The following is a discussion of the equip- ment and techniques used. PREVIOUS WORK The study of geophysical properties of earth materials under hydrostatic pressure has previously been largely re- stricted to the pressure range of less than two kilobars, with the exception of some velocity studies and some work in electrical conductivity (Stesky and Brace, 1973). Experi- ments in solid-state physics have gone into hundreds of kilobars, but these studies usually begin at approximately 10 kilobars, and are not done on earth materials. The pur- pose of the present study is to partially fill the 2-—10 Kb gap. Until recently, research into the magnetic properties of minerals has generally been restricted to tests under uniaxial pressures. However, effects of high hydrostatic pressure on the Curie temperature of magnetite (Schult, 1970), and on coercive force and saturation remanence (Carmichael, 1969) have been observed. PRESSURE MICROBOMB The pressure vessel used in this study is essentially that described in previous studies (Carmichael, et al, 1968). It is a piston apparatus (Fig. II-l) of beryllium copper 14 15 NO N- MAGNETIC PRESSURE CHAMBER AND BALLISTIC MAGNETO METER PISTON SENSE COIL / LOCKING CAP I“ VT >- m J J O D: *— SAMPLE I I‘ "I I I | SOLENOID __I |__ _I . N \ _ ‘Q X ‘ 3 T l' I I ‘K PISTON SUPPORT CYLJNDER FIG. ”-1 16 alloy (1.81% Be by weight) and is capable of sustaining internal pressures to 11 kilobars. The bomb is diamagnetic, with a susceptibility of -0.6x10.7 emu/gm at 20°C, so magne- tic changes due to the bomb with variations in applied fields are negligible compared to changes in magnetite and other samples. Aside from its non-ferromagnetic material, the bomb has a number of other unique attributes. It is capable of pro- ducing both hydrostatic pressure and a uniaxial pressure on the sample simultaneously, and needs only a simple uniaxial press to achieve this. The bomb is small enough that magne- tic measurements may be made with a small solenoid during pressurization, and it is easily portable. The pressure may be held in the bomb for long periods of time by means of the locking caps, which may be screwed down against the mushroom pistons on either end, thus making possible the study of long-term pressure effects. The present version of the bomb has a larger sample space than the original, allow- ing the study of samples up to 1.2 cm in diameter and 2.5 cm in length. After machining, the bomb used here was heat treated at 3160C for four hours. After heat treating, the parts were work hardened by the method outlined by Paul, et al (1959). All materials were work hardened up to 5.5 Kb. Following several pressure runs to four kilobars, the large piston's mushroom cap sheared around the base. The outer ring of the cap had been supported by the washers and l7 shoulder of the main piston, while the inner portion of the cap was forced down into the recess in the outer piston by the confining pressure (see Fig. 11-2). This was corrected by modifying the piston and cap: the flanged ends of each cap were increased in thickness, and the length of the neck of the cap was increased. The modification corrected the problem with the mushroom caps, but the stress previously absorbed by the caps was now transferred to the shoulders of the main pistons. No ill effects were noted at pressures below four kilobars, but at pressures between 4.0 and 5.5 kilobars the main piston of the large piston assembly began to deform, expanding the portion of the piston between 0.2 and 0.5 centimeters below the bottom of the recess. This expansion made removal of the piston assembly difficult, and greatly increased the friction between the cylinder wall and the piston during pressure runs, thereby adding an uncertainty in the value of the internal confining pressure with increasing applied press force. The piston was remachined to its original diameter, but this weakened it enough to cause shearing in the previously ex- panded area. The large piston was again modified, this time by re- ducing the diameter of the recess and corresponding neck of the mushroom cap. This allowed the stress to be absorbed by a greater area of the piston, and to be concentrated more toward the center of the piston instead of at its perimeter. The final versions of the pistons are shown in Figure II-Z. LARGE PISTON ASSEMBLY PISTON [11 Be— Cu WASHERS (Ill /.<-— LEAD GASKET ‘— RUBBER GASKET , 96-h Hm ALL PIECES HAVE 4— MUSHROOM CAP CIRCULAR CROSS—SECTION SMALL PISTON ASSEMBLY <— PISTON /"D Be—Cu WASHERS \h-EI 84— Pb GASKET ‘3 9"- RUBBER GASKET <—- MUSHROOM CAP SCALE ‘22:: 0 1,200 1" PISTON ASSEMBLIES FIG. “—2 19 One final modification of the previous bomb was made. When extracting the pistons and sample from the bomb after pressurization, a force of several hundred pounds is required to drive out the large piston. Single crystal samples gen- erally cannot sustain such large directed stress without crushing. To prevent destroying the samples after each pres- surization, a brass bushing was made to fit around the sample. Thus, when the large piston is driven out by forcing the small piston through the entire length of the bomb, this force is absorbed by the bushing and not the sample. Small notches were filed in each of the bushings to allow free fluid flow around the sample during pressurization. OPERATION Pressurizing techniques for the current research closely followed tyose outlined for the original bomb (Carmichael, et al, 1968). A fluid pressure-transmitting medium of a 1:1 mixture of kerosene and transformer oil was used. A fluid pressure medium makes for more difficult sealing, but some advantages are gained over solid pressure-transmitting medi- ums: uniform confining pressure is assured, so there is no concern about possible pressure gradients (Paul, et al, 1959) and though the oil mixture congeals at low temperature, there is a low fractional change in volume and therefore little hydrostatic pressure loss (Carmichael, et a1, 1968). The sample was first coated with resin to prevent oil from entering any small cracks and causing the sample to disintegrate under pressure. The sample was then glued to 20 the base plate and inserted into the microbomb, One piston was then driven into the cylinder far enough to allow the corresponding locking cap to be secured. The chamber was then filled with oil, the other piston enserted and the second locking cap secured. Seals around the pistons were provided by two gaskets, one of rubber and one of lead. The rubber gasket expands against the wall of the cylinder upon slight pressurization, and may be tightened by hand to a sufficient degree to pre- vent leakage. As pressure is applied by a press, the lead gasket is compressed between the two Be-Cu washers, flowing out around the piston to form a seal sufficient to withstand the pressures involved (Fig. II-Z). Internal hydrostatic pressure was determined using a sample of NH4F. This crystal undergoes a change of phase at 3.65 Kb, and reduces in volume by about 30%. The sample of NH4F was sealed by an impermeable rubber membrane, and pis- ton displacement was monitored versus applied press pressure. A typical calibration curve Of several tests for this micro- bomb is given in Figure 11-3. The sudden change in the slope of the curve indicates the point where internal pressure has reached 3.65 Kb; the sample has suddenly reduced its volume, and the pistons have moved in rapidly to maintain a constant pressure. The half-width of the hysteresis loop of each pressure cycle is approximately equal to the packing friction between the large piston and the cylinder wall (Paul, et al, 1959). Internal hydrostatic pressure is assumed to be a 21 “3.22 0253 mEOmOmQE wmammmcn. “.0 20.53.5943 4(03: ( T: 9.50:. .3 82x memo“. is m. E m. a. : o. a m A o o _ _ _ _ _ L _ _ _ o I I _ 386 new I oz_m womOu 8mm... 03:42.32... T: 230: 8.. 8N. + memo". 33%? u . 3.. 82x. memo“. is a «N cm 2 2 3 «F o. m o v N o o — L _ _ _ A _ A _ a I P w I N 3 H N V fil )1 m mm 0 4. 8 | .50 02.93.. w m I m> / is. m a U a V m>mao 02.96.23 I m m \ \ H \ \ 3 \ \ \ I \ \ \ I o \ \ \ \ \ [ _ _ L.\ _ _ _ _ _ _ _ _ h 24 I: .21 224.. me do Emmmao :53. 3m... .:oo 32% .20 2: mice “5 :mhzmu Pr momzwm mho2 + mmmmm ) m0hI malhwm .PstESOm ".0 E mow 20....