The search was performed with the function FindRoot in Mathematica ... obtain the pdf( - ), which then was convolved with the derivative of a Gaussian in.
Auxiliary Material for On the magnetocrystalline anisotropy of greigite (Fe3S4) Michael Winklhofer1, Liao Chang2,3, and Stephan H. K. Eder1 1
Department of Earth and Environmental Sciences, Ludwig-Maximilians-University Munich, Theresienstr. 41, D-80333 Munich, Germany
2
Research School of Earth Sciences, The Australian National University, Canberra, ACT 0200, Australia 3
Paleomagnetic Laboratory βFort Hoofddijkβ, Department of Earth Sciences, Utrecht University, Budapestlaan 17, 3584 CD Utrecht, The Netherlands
Geochemistry, Geophysics, Geosystems, 2014 Introduction The auxiliary material consists of an Appendix (Sections Ax) and supplementary slides. Figures A1, A2 and Tables A1-A3 are included in this README. Structure of Appendix: A1. Higher order cubic magnetocrystalline anisotropy terms a) Series representations of πΉ!" and definitions of πΎ! , πΎ! , πΎ! Β b) Easy axes orientation along symmetry axes and non-major cubic axes (Figure A1a,b) c) Temperature dependence of the πΎ!
A2. Computational details of the fitting procedure a) Calculation of resonance fields b) Forward calculation of powder spectra c) Fitting of FMR powder spectra d) Performance
A3. Results from spectral fitting a) Sensitivity study: Role of uncertainty in g-factor (Table A1a,b) b) Fits for X-band data of cement-embedded samples recorded at 290 K (Table A2a,b) c) Comparison between X- and Q-band FMR data for sample SYN519 (Figure A2 and Table A3)
Auxiliary slides βTEM_EDβ (Fig_SI_Revisiting_ED_YamaguchiWada.pdf) The supplement shows our re-analysis of the electron diffraction pattern that Yamaguchi and Wada [1970] considered as evidence of a easy axes orientation in greigite. We suggest that the reflections in question are not due to (400), but (440). Either way, the incompleteness of the original diffractograms does not allow for a texture analysis.
Appendix A1 Higher order cubic magnetocrystalline anisotropy terms a) Series representations of πππ§ and definitions of π²π , π²π , π²π Β The cubic magnetocrystalline anisotropy energy (MAE) density is usually expanded in even powers of directional cosines πΌ! of the magnetization (i=1,2,3), πΉ!" πΌ! , πΌ! , πΌ! ; πΎ! , πΎ! , πΎ! , β¦ = πΎ! + πΎ! πΌ!! πΌ!! + πΌ!! πΌ!! + πΌ!! πΌ!! + +πΎ! πΌ!! πΌ!! πΌ!! + πΎ! πΌ!! πΌ!! + πΌ!! πΌ!! + πΌ!! πΌ!! + β― ,
(A1) Β
where πΎ! , πΎ! , and Β πΎ! are referred to as first-, second-, and third-order magnetocrystalline ansiotropy constants, respectively. Importantly, πΉ!" is the Helmholtz free energy, which is measured at constant temperature and volume (strain free), which should hold reasonably well for particles embedded in a mechanically strong matrix. The following representation for the MAE is also found in the literature, πΉ!" β²(πΌ! , πΌ! , πΌ! ; πΎ! , πΎ!! , πΎ! ) = πΎ! + πΎ! πΌ!! πΌ!! + πΌ!! πΌ!! + πΌ!! πΌ!! + +πΎ!! πΌ!! πΌ!! πΌ!! + πΎ! πΌ!! πΌ!! + πΌ!! πΌ!! + πΌ!! πΌ!! ! ,
(A2) Β
which because of πΌ!! + πΌ!! + πΌ!! = 1 Β is equal to πΉ!" (πΌ! , πΌ! , πΌ! ; πΎ! , πΎ!! , πΎ! ) = πΎ! + πΎ! πΌ!! πΌ!! + πΌ!! πΌ!! + πΌ!! πΌ!! + + πΎ!! + 2πΎ! πΌ!! πΌ!! πΌ!! + πΎ! πΌ!! πΌ!! + πΌ!! πΌ!! + πΌ!! πΌ!! ,
(A3)
so that πΎ! of (A1) corresponds to πΎ!! + 2πΎ! of (A2). This mixing of coefficients for the πΌ!! πΌ!! πΌ!! term in Eq. (A3) is confusing and can be avoided by using (A1). b) Easy axes orientation along symmetry axes and non-major cubic axes Depending on the relative magnitude of πΎ! , πΎ! , πΎ! , any of the crystallographic axis of high symmetry ( β four fold, - two fold, and three fold) can be an easy, intermediate, or hard axis (see Fig. A1a for πΎ! >0 and Fig. S1b for πΎ! 0 is the minimum of the MAE, while it is the maximum for πΎ! 0), then there is a region in the πΎ! /πΎ! , πΎ! /πΎ! plane (above the red line in Fig. A1b), where the easy axis lies between a and an adjacent axis (or between a and axis, black field). Such a more general (or ) orientation of the easy axis was indeed observed in some cubic Laves compounds RFe2 (R=Ce, Lu) [Atzmony & Dariel, 1974].
Figure A1a βPhase diagramβ showing the energetic ordering of symmetry axes (easy axis < intermediate axis < hard axis) for K 1 > 0 with the magnetocrystalline anisotropy energy given by (A1). Each one of the symmetry axes can be an easy, intermediate, or hard axis.
Figure A1b As in Fig A1a, now for K 1 < 0. As opposed to the case K1 > 0, the easy axis now need not coincide with a major symmetry axis. Above the red line, the easy axis lies between a axis and axis and has a [uuw] orientation. In the black field, the easy axis is between and . The relative ordering of the symmetry axes in terms of energy is not affected by the tilt. Below the red line, the (left) or axes define the easy axes of magnetization.
c) Temperature dependence of the π²π Although representation (A1) of the MAE is compact and useful for computational purposes, it is not the most appropriate representation for analysis of the temperature dependence of the MAE in terms of the Akulov-Zener-law [see Callen and Callen, 1966, J. Phys. Chem. Solids, 27, 1271-85] which states that π! (π) βπ(βπ!!)/! π
βπ (π) Β = , (A4) π
βπ (0) Β π! (0) Β where the π
βπ are magnetocrystalline anisotropy coefficients, which are related to an expansion of the MAE in terms of spherical harmonics. Including π
βπ coefficients up to order βπ = 6, Zener [1954, Phys. Rev. 96, 1335-37;] gave the following expression for a βzero-meanβ cubic MAE ! πΉ!" in terms of directional cosines ! πΉ!" πΌ! , πΌ! , πΌ! ; π
! , π
! = π
! πΌ!! πΌ!! + πΌ!! πΌ!! + πΌ!! πΌ!! β
+π
! πΌ!! πΌ!! πΌ!! β where π
! = πΎ! +
! !!
1 2 πΌ!! πΌ!! + πΌ!! πΌ!! + πΌ!! πΌ!! + , 11 231
1 +β― 5 (A5)
πΎ! and π
! = πΎ! . The numerical constants in (A3) (1/5 and -2/231) are
! chosen such that the average of πΉ!" πΌ! , πΌ! , πΌ! ; π
! , π
! over the surface of the unit sphere ! ! ! ! vanishes, while the average of πΉ!" πΌ! , πΌ! , πΌ! ; πΎ! , πΎ! is πΎ! + πΎ! . Note that β β
β
! !"#
!
!"#
!"#
!
! !!
=
.
When πΎ! cannot be neglected in (A1), then πΎ! will contribute to π
! and π
! and equation (A5) needs to be expanded to a higher order. The relationships between the π
βπ and πΎ! can be obtained with the procedure detailed in Callen and Callen [1960, J. Phys. Chem. Solids, 16, 310-328, Appendix B therein], which they carried out to order βπ = 6. When including πΎ! (i.e., going to order βπ = 8), we obtain the following relationships: 1 24 1 18 πΎ! + πΎ! , Β Β πΎ! = π
! β π
! β π
, 11 143 11 65 ! 6 6 Β (A6) Β Β Β Β π
! = πΎ! β πΎ! , πΎ! = π
! + π
! , 5 5 π
! = πΎ! , πΎ! = π
! ,
π
! = πΎ! +
!
and π
! = πΎ! + 21πΎ! + πΎ! + 3πΎ! . Therefore, if πΎ! is similar in magnitude to πΎ! , it has a !"# stronger influence on the anisotropy than πΎ! has. Hence the temperature dependence of πΎ! and πΎ! predicted by (A4) then is: 1 24 π! π πΎ! π! + πΎ! π! π! π! 11 143 !" 6 18 π! π π! π β πΎ! π! β Β πΎ! T! π! π! π! π! 5 65
!"
!"
!"
πΎ! (π) Β = πΎ! π! + β
1 Β πΎ π 11 ! !
and 6 πΎ! π = πΎ! π! β πΎ! π! 5
π! π π! π!
6 + Β πΎ! T! 5
π! π π! π!
+ !"
, (A7)
,
(A8)
For greigite, Chang et al. [2008] reported a 7-8% decrease of M! from π! =5 K to room temperature.
A2: Computational details of the fitting procedure a) Calculation of resonance fields: The resonance field π΅!"# (π! , π! ) for a given set of anisotropy parameters and for a given orientation of the external field vector with respect to the [001] axis of the crystal was obtained by simultaneously solving (numerically) the resonance condition (Eq. 1 in paper) and the equilibrium condition (ππΉ/ππ = 0, ππΉ/ππ = 0) for π β π!" , π β π!" , and π΅ β π΅!"# . The search was performed with the function FindRoot in Mathematica (Wolfram Research). The initial guesses for π!" and π!" were π! and π! , respectively, and the search range for π΅!"# Β was restricted to the range defined by the minima and !!!!! !!!"! !!""! , π΅!"# , π΅!"# } (see Eq. 3 in paper), except for maxima of the set SB = {π΅!"# or easy axis orientations. As initial guess for π΅!"# we chose the value from the analytical first-order approximation derived by SchlΓΆman [1958], provided that value was inside the search range. In case it was not, the intermediate value of SB was used. b) Forward calculation of FMR powder spectra: π΅!"# (π! , π! ) was first computed according to A2.a) for 100 points in the field ! ! orientation range π! = 0, , π! = [0.01, ]. The π΅!"# (π! , π! ) surface was then !
!
interpolated on a finer grid (100 x 100 points) equidistant in π! and cos π! , using an interpolation order of 3. The resulting π΅!"# π! , cos π! dataset was histogrammed to obtain the pdf(π΅!"# ), which then was convolved with the derivative of a Gaussian in order to obtain a derivative absorption spectrum. The convolution was performed in Fourier space, using the discrete Fourier transform. c) Fitting of FMR powder spectra The best-fit model parameters K1, K2, K3, g, Ξπ΅ were obtained by minimizing the sum of squared differences between experimental and model spectra, using the NMinimize function in Mathematica with the NelderMead method (a simplex scheme). Both experimental and simulated spectra were normalized in the same way [ (max-min)/2 norm] so as to avoid another fit parameter. Control minimizations were performed with the computation-time intensive SimulatedAnnealing method to check if there were other parameter combinations in the global search space which had not been found with the simplex method under the given initial guesses.
d) Performance: The simulation of a derivative absorption spectrum for a given anisotropy model takes about a second on a MacBook Air. The fitting of an experimental spectrum may take several hours (for the full set of five fitting parameters) depending on the distance between initial guesses and final values, but can be reduced to about 10 - 40 minutes when two of the five model parameters have prescribed values (i.e., K3=0, g=2.08). Of course, the run time depends on the requested precision and on how close the initial guesses are to the final guess.
A3. Results from spectral fitting a) Sensitivity study: Role of uncertainties in g-factor The intrinsic g value (βg-factorβ) of greigite is not known, but likely to be in the range 2.05 to 2.1 (see paper). To find out how uncertainties in g affect the fitted value of K1, we forward calculated FMR absorption derivative spectra for a given set of anisotropy parameters K1, K2, g and effective intrinsic Gaussian linewidth Ξπ΅ and then fitted the spectra with prescribed values of g. The results in Table A1a,b show that overestimation of the g value leads to an underestimation of the magnitude of K1 and vice versa. The fit error parameter s depends quadratically on (gest β gtrue), where gtrue and gest are the true and estimated values of g. g πΎ!!"#
2.04
2.05
2.06
2.07
2.08
2.09
2.1
2.12
2.14
Β Β
-2.987
-2.916
-2.847
-2.773
-2.701
-2.630
-2.560
-2.405
-2.249
Β Β
-0.578
-0.586
-0.582
-0.586
-0.583
-0.582
-0.583
-0.595
-0.590
Ξπ΅ [mT]
57.6
58.5
58.8
59.4
60.0
60.7
61.3
62.9
64.4
Ξ΄π΅!" [mT]
4.65
5.28
6.01
6.84
7.67
8.48
9.25
11.17
13.37
s
0.020
0.01
0.005
0.001
0. Ξ΄π΅!" = π΅!"# β π΅!"
ii) for K1 0 Β models for sample SYN519 is twice as large as that for sample SYN702 (compare Table 1 in paper with Table A3 here), which is readily explained by the different π factors (1.93 vs. 2.08). As a consequence, the anisotropy field is twice as large, too. Even so, the energy barriers ΞπΉ!" remain at about 2.5 kJ/m3, irrespective of the g factor. As far as the πΎ! < 0 models are concerned, the magnitude of πΎ! in the models for sample SYN519 is between four and five times larger compared to the models for sample SYN702 (compare Table 1 and Table A3), and most strikingly, πΎ! now has large positive values instead of small negative values. Yet, as with the πΎ! > 0 models for SYN519, the energy barriers Β βπΉ!" calculated for the πΎ! < 0 models do not differ much from those listed in Table 1. The large negative πΎ! values obtained are related to the high field side of the FMR derivative spectrum (Figure A2), which at Q-band frequencies has a long negative tail that extends beyond the maximum field of 2 T. This tail is even more pronounced for sister sample SYN706 from Chang et al. [2011].
Model
πΎ! Β [kJ/m3]
πΎ! πΎ!
πΎ! πΎ!
g
a519
β68.9
β2.2
β1.3
2.08
b519
β72.8
β2.4
β1.2
2.05
c519
β83.3
β2.2
β0.3
1.89
d519
β74.6
β1.5
β0.3
1.90
e519
β69.4
β1.8
0
1.90
f519
β33.3
β0.4
β0.2
2.04
519
β80.4
β5.7
β2.4
1.92
a519
+16.3
+3.5
0
2.08
b519
+16.8
+4.1
β0.8
2.04
c519
+16.7
+5.6
+0.7
2.04
d519
+18.5
+4.5
β0.3
2.04
f519
+20.6
+9.8
β2.1
1.93
f [GHz] 9.4 (34.0) (9.4) 34.0 9.4 34.0 9.4 34.0 9.4 34.0 (9.4) 34.0 9.4 34.0 (9.4) 34.0 (9.4) 34.0 (9.4) 34.0 (9.4) 34.0 9.4 34.0
Ξπ΅ [mT] 73 (95) (71) 112 77 91 66 105 81 103 (70) 100 61 92 (68) 80 (63) 88 (75) 90 (67) 89 65 88
fit categ. 1 3 2 2 3 2 2 1 2 2 3 3 1 2 3 2 4 3 4 2 3 2 2 2
π΅!" [mT]
ΞπΉ!" [kJ/m3]
104
2.5
82
1.9
116
0.6
205
2.6
150
1.1
160
2.5
208
1.0
135
4.0
140
3.3
241
4.9
154
4.2
171
2.4
Table A3 Overview of anisotropy models A! = πΎ! , πΎ! , πΎ! ; π ! found for greigite powder SYN519 dispersed in eicosane. Symbols as in Table 1 of paper. The column βfit categoryβ (1 = very good) is a qualitative rank for the goodness of fit for a given solution at X-band (9.40 GHz) and Q-band (34.0 GHz). For illustration, see Figure A2. For most of the cases, experimental data were only fitted in one band, so that the values for ΞB in the other band indicate approximate values for ΞB and are therefore set in parentheses.
Figure A2 (Following three pages) X- and Q-band data for greigite powder SYN519 dispersed in eicosane (solid line), with simulated spectra (red dots) for the models listed in Table A3. All spectra are normalized.
aS519 Γ» Q-band
1.0
1.0
0.5
0.5
dc''ΓͺdB
dc''ΓͺdB
aS519 Γ» X-band
0.0
-0.5
-0.5 -1.0 0
0.0
200
B @mTD 400
600
-1.0 600
800
bS519 Γ» Q-band
1.0
1.0
0.5
0.5
dc''ΓͺdB
dc''ΓͺdB
bS519 Γ» X-band
0.0
0.0
-0.5
-0.5
-1.0 0
-1.0 600
200
B @mTD 400
600
800
1.0
1.0
0.5
0.5
0.0
0.0 -0.5
-0.5 200
B @mTD 400
600
-1.0 600
800
dS519 Γ» X-band 1.0
0.5
dc''ΓͺdB
dc''ΓͺdB
B @mTD
800 1000 1200 1400 1600 1800 2000
dS519 Γ» Q-band
1.0
0.0
0.5 0.0 -0.5
-0.5 -1.0 0
B @mTD
800 1000 1200 1400 1600 1800 2000
cS519 Γ» Q-band
dc''ΓͺdB
dc''ΓͺdB
cS519 Γ» X-band
-1.0 0
B @mTD
800 1000 1200 1400 1600 1800 2000
200
B @mTD 400
600
800
-1.0 600 800 1000 1200 1400 1600 1800 2000
B @mTD
eS519 Γ» Q-band
1.0
1.0
0.5
0.5
dc''ΓͺdB
dc''ΓͺdB
eS519 Γ» X-band
0.0
-0.5
-0.5
B @mTD fS519 Γ» X-band
200
400
600
-1.0 600 800 1000 1200 1400 1600 1800 2000
B @mTD fS519 Γ» Q-band
800
1.0
1.0
0.5
0.5
dc''ΓͺdB
dc''ΓͺdB
-1.0 0
0.0
0.0
0.0 -0.5
-1.0 0
-1.0 600 800 1000 1200 1400 1600 1800 2000
B @mTD S519 Γ» X-band
200
400
600
B @mTD S519 Γ» Q-band
800
1.0
1.0
0.5
0.5
dc''ΓͺdB
dc''ΓͺdB
-0.5
0.0
-0.5
-0.5
B @mTD aS519 Γ» X-band
200
400
600
-1.0 600 800 1000 1200 1400 1600 1800 2000
B @mTD aS519 Γ» Q-band
800
1.0
1.0
0.5
0.5
dc''ΓͺdB
dc''ΓͺdB
-1.0 0
0.0
0.0
0.0
-0.5
-0.5
-1.0 0
-1.0 600 800 1000 1200 1400 1600 1800 2000
200
B @mTD 400
600
800
B @mTD
bS519 Γ» Q-band
1.0
1.0
0.5
0.5
dc''ΓͺdB
dc''ΓͺdB
bS519 Γ» X-band
0.0
0.0 -0.5
-0.5 -1.0 0
200
B @mTD 400
600
-1.0 600
800
1.0
1.0
0.5
0.5
dc''ΓͺdB
dc''ΓͺdB
cS519 Γ» X-band
0.0
0.0 -0.5
-0.5 -1.0 0
200
B @mTD 400
600
800
-1.0 600
1.0
1.0
0.5
0.5
dc''ΓͺdB
dc''ΓͺdB
dS519 Γ» X-band
0.0
-0.5
-1.0 0
-1.0 600
B @mTD eS519 Γ» X-band
200
400
600
800
0.5
0.5
dc''ΓͺdB
1.0
0.0
0.0
-0.5
-0.5
-1.0 0
-1.0 600
B @mTD 400
600
B @mTD
800 1000 1200 1400 1600 1800 2000
eS519 Γ» Q-band
1.0
200
B @mTD dS519 Γ» Q-band
800 1000 1200 1400 1600 1800 2000
0.0
-0.5
dc''ΓͺdB
B @mTD cS519 Γ» Q-band
800 1000 1200 1400 1600 1800 2000
800
B @mTD
800 1000 1200 1400 1600 1800 2000