CC BY
16
Calculations of the onset temperature for tunneling in multispin systems
S.M. Vlasov1'2, P. F. Bessarab1'2, V.M. Uzdin1'3, H. Jonsson2'4 xITMO University, Kronverkskiy, 49, St. Petersburg, 197101, Russia 2 Science Institute and Faculty of Physical Sciences, Univ. of Iceland, 107 Reykjavik, Iceland 3St. Petersburg State University, St. Petersburg, 198504, Russia 4Department of Applied Physics, Aalto University, Espoo, FIN-00076, Finland
v_uzdin@mail.ru
PACS 75.45.+j , 82.20.Xr DOI 10.17586/2220-8054-2017-8-4-454-461
Transitions between magnetic states of a system coupled to a heat bath can occur by exceeding the energy barrier, but as temperature is lowered quantum mechanical tunneling through the barrier becomes the dominant transition mechanism. A method is presented for estimating the onset temperature for tunneling in a system with an arbitrary number of spins using the second derivatives of the energy with respect to the orientation of the magnetic vectors at the first order saddle point on the energy surface characterizing the over-the-barrier mechanism. An application to a monomer and a dimer of molecular magnets containing a Mn4 group is presented and the result found to be in excellent agreement with reported experimental measurements.
Keywords: magnetic transitions, tunneling, rate theory, molecular magnet.
Received: 5 August 2017 Revised: 15 August 2017
1. Introduction
An assessment of the stability of magnetic states is an important problem in the theory of magnetism [1,2]. Magnetic nanosystems are candidates for spintronics and high density storage devices where a critical issue is the lifetime of prepared magnetic states. A given magnetic state of a nanosystem coupled to a heat bath can have a finite lifetime because of thermally activated transitions involving leaps over an energy barrier separating it from other magnetic states of the system. But, at low enough temperature, quantum mechanical tunneling through the barrier will become the dominant transition mechanism. It is important to have a tool to estimate the temperature at which the tunneling mechanism becomes predominant, thereby reducing the lifetime from what would be expected based on the over-the-barrier mechanism.
Quantum tunneling in spin systems has been the topic of many theoretical [3] and experimental studies [4-7] over the past few decades. At high enough temperature, quantum effects are not significant and transitions are due to thermal activation described by the Arrhenius law. However, in the limit of zero temperature, only tunneling out of the lowest energy level of the initial state can occur. In an intermediate temperature range, tunneling can be thermally assisted, i.e. occurring from thermally activated levels of the initial state. The crossover from over-the-barrier to quantum tunneling can, in such cases, be relatively smooth. The shape of the energy surface affects how sharp this crossover is [8].
Several theoretical studies of the crossover in simple spin models containing one [9-12] or two [13] spins have been carried out using the effective potential method [14]. The effective potential method uses an exact mapping of a spin system onto a particle system with an effective potential and mass that can be coordinate dependent. However, this method is not general and can only be applied to simple model systems. High-order anisotropy terms, such as (S+ + S—) and (S+ + S—), have not been incorporated in the mapping approach [15], but can significantly affect the tunneling rate [5,16]. So far, systems with high-order anisotropy have only been studied numerically by direct diagonalization of the Hamiltonian [17,18].
Here, a method for calculating the onset temperature of thermally assisted tunneling in a multispin system is presented and applied to two systems: a monomer and a dimer of molecular magnets containing Mn4 groups. The method involves an analysis of the shape of the energy surface in the vicinity of the first order saddle point on the energy surface representing the transition state of the over-the-barrier mechanism. The present article extends a previous report on systems consisting of a single magnetic vector (a single spin, or multiple spins within the macro-spin approximation) [19]. With the present formulation, the onset temperature for tunneling in systems consisting of an arbitrary number of spins can be estimated.
The article is organized as follows: In the next section, the methodology for estimating the lifetime of a magnetic state due to over-the-barrier transitions is briefly reviewed. Then, the expression for the onset temperature for quantum mechanical tunneling is derived in section 3. Applications are presented in section 4, first to a monomer and then to a dimer of molecular magnets. The article concludes with a summary in section 5.
2. Over-the-barrier transitions
A system of N spins is described by the orientation of each of the spins, w = {61, 62,... ^n, }, where 6 and ^ are spherical polar coordinates. The energy as a function of these variables, U(w), represents an energy surface and the magnetic states of the system correspond to local minima on this surface. The mechanism of a transition from one state to another involving a jump over an energy barrier is characterized by a minimum energy path (MEP) connecting the corresponding minima. At every point on an MEP, the gradient of the energy point is along the path. The MEP is a transition path with maximal statistical weight, assuming a Boltzmann distribution is established and maintained in the initial state. The MEP can be found using the geodesic nudged elastic band method where the path is discretized by generating a set of replicas of the system and an optimization algorithm is used to bring the replicas from some initial location to an MEP [20].
The point of highest energy along the MEP, wt, represents a first order saddle point on the energy surface, a point where the Hessian has one and only one negative eigenvalue. Within harmonic transition state theory (HTST), the transition state is taken to be a hyperplane going through the saddle point with normal pointing along the direction of the unstable mode (the eigenvector corresponding to the negative eigenvalue of the Hessian) [21,22]. An estimate of the activation energy for the transition is then obtained as AE = U(w^) — U(wm), where U(w^) is the energy at the first order saddle point and U(wm) is the energy of the initial state minimum.
The rate of thermally activated transitions involving jumps over the energy barrier can be estimated using HTST for magnetic systems [21]
rHTST = VQe-AE/kB T, (1)
where v0 is a pre-exponential factor that depends on the Hessian evaluated at the saddle point and at the minimum of the energy surface. Calculations of transition rates in magnetic systems using this approach have been carried out in studies of, for example, remagnetization in small islands adsorbed on a solid surface [23], analysis of temperature dependence of hysteresis loops [24], annihilation of magnetic skyrmions [25,26] and interaction of a magnetic tip with a solid surface [27].
Quantum mechanical effects can influence the transition rate even at temperature above the onset temperature for tunneling. This occurs because of the zero point energy (ZPE) which can raise the energy of the initial state and the transition state to a different extent. Within the harmonic approximation, the correction is obtained by adding the ZPE for each vibrational (magnon) mode to the initial state energy and the saddle point energy. As there is one fewer vibrational mode at the transition state than at the initial state, this effect tends to reduce the energy barrier and lower the activation energy for the transition, analogous to the so-called kinetic isotope effect in atomic rearrangements (chemical reactions and atom diffusion). The ZPE corrected activation energy is
AEzpec = (u(wt) + 2 NC — (u(wm) + 2 £ , (2)
where v^ and vm i are the vibrational frequencies [28,29] calculated at the saddle point and at the minimum, respectively. As can be seen from the molecular magnet example below, the ZPE correction can significantly change the acitvation energy and thereby the slope in the Arrhenius graf.
3. Onset temperature for tunneling
Tunneling through an energy barrier can take place because of quantum mechanical delocalization which enables the system to avoid the top of the energy barrier. In a thermalized system, the effect of delocalization can be described in a convenient way using statistical Feynman Path Integrals. The instanton technique can be used to identify when quantum delocalization becomes large enough to make tunneling the dominant transition mechanism. In instanton theory the transition rate is written in terms of the imaginary-time (t = it) action, S, and a stationary phase approximation [30-32]. The transition rate constant can be expressed as
r = D exp(-S[w(t)]inst/h), (3)
where w(t) is a path corresponding to a periodic trajectory with period t = ¡3h satisfying w(0) = w(p%), and the pre-factor D takes into account the probability of paths deviating from w(t). This special path, often referred to
as the instanton, corresponds to a stationary point of the action, a saddle point on the action surface where the path is extended, corresponding to quantum delocalization [32-34].
In our previous work [19], we considered systems with a single spin. Here, a generalization to many-spin systems is described. The Hilbert space of a many-spin system is a tensor product of the Hilbert space of the individual spins. The Euclidean (imaginary-time) action for a system with N spins of length s is defined by [35-37]
/3/2 N /3/2
S (H) = -is J ¿(1 - cos dr + J U (w(r))dr
(4)
-//2 j 1 -//2 where the first term is the sum of the Berry phase of the individual spins and H = w(r) represents a closed trajectory with period p.
From here on we use atomic units with h = 1, = 1/2 and the mass and charge of an electron have unit value, e = 1 and me = 1. A scaled gyromagnetic ratio is defined as g = g/2.
In order to find the onset temperature for tunneling, Tc, the action is expanded to second order in the vicinity of the first order saddle point on the energy surface, H = + eSH,
S(H) = pU (wt) + SS +1 S2S, (5)
2
/3/2 N /3/2
S2S = -2is J E (ôOjS4>j sin j dr + J SH H SHTdr,
-3/2 j = 1
(6)
-3/2
where H is the Hessian matrix of the energy surface U(w) at the saddle point wt. As e ^ 0, S2S becomes a quadratic form of the Hessian, H, which has one and only one real negative eigenvalue.
The task is to find the temperature at which ^ becomes an instanton, i.e. a saddle point at the action surface with quantum delocalization. At that point, a zero mode appears corresponding to displacement along the path and thus constant S. This signals the transition from thermally activated jumps to quantum tunneling [32,33].
Since the instanton is a closed path, S6k and S^k can be expanded in a Fourier series and S2S in eqn. (6) rewritten as
N
2 S2S^ ,^) = p ¿IJ2
n=0 \fc=0
2ns sin 6
p
k n(^6f - 6kn)
+
E
a6k 6n + b{4>n 6j + C 6j) + c4>n4>n
(7)
The matrix representing the quadratic form of the action is infinite and has a block diagonal form
/70 \
G =
Yi
(8)
.
where
Yn = H +
0 • • 0 -kn •• • 0
0 . . 0 0 . .0
0 • •0 0 • • kn • kN
kn • •0 0 • • • 0
0. .. 0 0 . .0
0 • kn kN 0 • • • 0
(9)
. / where = 2nnTs sin j
An expression for the onset temperature for tunneling can be obtained by analyzing the determinant of G:
det(G) = H det(Yn).
(10)
n=0
Y
n
Consider first 70 which is the Hessian of the energy surface at the first order saddle point. There, 70 contains one negative eigenvalue so det(Y0) < 0. However, the other blocks, y„ (for n > 0), depend on temperature and at T > Tc are positive definite
N
det(7n) - (2nsn)2N sin2 0t > 0,
(11)
j=i
and det(G) < 0. As the temperature is lowered to Tc, a zero of G appears signaling a stationary point involving quantum delocalization. This must arise from 71 since the temperature in each y„ is scaled by n and blocks with n > 1 thus have zeros at even lower temperature. The onset temperature for tunneling is, therefore, found by solving
det(7i(Tc))=0. (12)
This can be easily done numerically for multi-spin systems. All that is required for the calculation are the second derivatives of the energy at the first order saddle point corresponding to the over-the-barrier mechanism. For a single-spin system, an analytical expression for the onset temperature can be obtained
\/b2 — ac
Tc =
where
d2U (0t,4>t)
362
2nskB sin 0t ' d2U (0t,4>t)
et,$t
34>2
b =
d2U (0t,4>t)
et,0t
3030
0t,0t
The pair of eigenvalues of the 71 block are
a + c J(a - c)2 + 4b2 - 4k2
(13)
(14)
(15)
2 2
In addition to the crossover temperature, we can define the temperature T0 at which the two eigenvalues become complex conjugates
■s/(a — c)2 + 4b2
To =
(16)
4nskB sin 01
At this temperature, both eigenvalues are equal to (a + c)/2 and could be negative. At higher temperature, the eigenvalues are complex. In the temperature range between Tc and T0 the two eigenvalues are real, either both are positive or both are negative.
a
c
4. Applications. Molecular magnets
Previously, we have presented calculations and compared with experimental results for the single molecular magnet containing a Mn4 group [19]. We briefly review those calculations here before discussing the dimer. The molecular magnet has a total spin of s = 9/2 and can be described with the following Hamiltonian
# = -DS2 — BS4 — gm0S • H + Htrans, (17)
where D and B are axial anisotropy constants, B is an applied magnetic field and Htrans gives the in-plane anisotropy which is small compared with D. The values of the parameters have been determined from experimental measurements (not involving tunneling) [38]: D = 0.74 K, B = —3.8 mK and g = 2. The transverse term is chosen to be fourth-order E(S+ + S- with E = 9.8 x 10-4 K [38].
For the single molecular magnet it is possible to give an analytical formula for the onset temperature. Expanding (12) with the Mn4 system Hamiltonian (17) at the saddle point gives
Tc =4sVDE — 4E 2s2, (18)
n
and Tc = 0.61 K. This result is in close agreement with the reported experimental results giving 0.67 K [38]. The dependence of eigenvalues on temperature is shown in Fig. 1.
The calculated transition rate using HTST (applicable for T > Tc) and the predicted onset temperature for tunneling in the monomer are shown in Fig. 2, along with experimental data. At temperature below Tc, the tunneling rate is expected to be similar to the HTST rate at Tc. The agreement between the calculations and measurements is good.
temperature, [K]
Fig. 1. Eigenvalues of the y1 matrix block for the monomer (eq. (17)) Mn4 molecular magnet. Solid line: real part. Dashed line: imaginary part. At temperature above T0 the eigenvalues are complex conjugates and give a positive contribution into the determinant. At temperature below the Tc, y1 has only one negative eigenvalue, therefore the determinant changes sign at Tc. In the temperature interval Tc <T <T0 y1 has two negative eigenvalues
Fig. 2. The calculated lifetime for a monomer and a dimer of molecular magnets containing Mn4 groups. Results of calculations using harmonic transition state theory with zero point energy correction for the monomer (solid gray line) and the dimer (solid black line) are shown. Calculated lifetime for the dimer without ZPE correction (dashed line) is also shown for comparison. Calculated values for the onset temperature for tunneling calculated from (18) using the Hamiltonian in (17) for the monomer and (19) for the dimer (dashed vertical lines) with parameters taken from experimental measurements are also shown. Experimentally measured transition rate for the monomer [38] is shown with gray triangles and for the dimer [40] with black triangles. Below Tc, the tunneling rate is expected to be similar to the HTST rate at Tc. Excellent agreement is obtained between the calculated and measured results. Inset: The determinant of the y1 matrix block for the dimer, showing how it becomes negative at Tc
We turn now to calculations for a dimer of Mn4 molecular magnets and compare with experimental measurements on [Mn4O3Cl4(O2CEt)3(py)3]2. Each unit of the dimer is described with a magnetic moment of s = 9/2, and the Hamiltonian of the system is taken to be
# = Hi + HH2 + JSi • S2, (19)
where both H1 and H2 have the same form as (17) with the parameters taken from [39] and J = 0.12 K is isotropic superexchange [40]. The transverse term is chosen to include both 4th- and 2nd-order anisotropies: C(S+ + S—) and E(S+ + S—), with C = 0.032 K and E = -7.5 x 10—5 K, in order to match the experimental results of the rate above Tc. The energy surface and a minimum energy path for a transition are shown in Fig. 3.
0 nl 4 n/ 2 3n/4 n
<t>8, trad]
Fig. 3. Energy surface for a dimer molecular magnets containing Mn4 group when dA = dB = n/2. The minimum energy path for a remagnetization transition is shown with black dots. The insets show the orientation of the 9/2 spins of each Mn4 unit at the initial state (left), first order saddle point (middle), and product state (right).
The calculated crossover temperature for the dimer is Tc = 0.40 K in close agreement with experimental measurements [40]. Due to the weak exchange interaction every spin flips almost independently.
in order to see the effect of coupling between the two units more clearly, we also calculated and analyzed a system with larger exchange interaction, taken to be of the same order as the anisotropy, J = 0.6 K. The calculated onset temperature for the tunneling is found to be Tc = 0.35 K. Due to the strong exchange interaction between spins eigenvalues of y1 behave differently, but we still can see the temperature at which one of them crosses zero (see Fig. 4).
The result of calculation of the HTST transition rate and onset temperature for tunneling of the dimer (19) using experimentally determined parameter values is shown in Fig. 2. The ZPE correction changes the slope in the Arrhenius graph, bringing the calculated HTST results closer to the experimental values.
5. Summary
A method is presented for finding the onset temperature of thermally activated tunneling in magnetic systems that are characterized by spin vectors with orientation prescribed by continuous angular variables. A study is made for the Mn4 molecular magnets where both the high temperature jump rate and the onset temperature for tunneling are calculated and compared with experimental data. Excellent agreement is obtained. The rate calculated above the onset temperature for tunneling includes zero point energy correction. The onset temperature for tunneling is found by analyzing eigenvalues of the y1 matrix block at the first order saddle point on the energy surface.
Previously, the onset temperature for tunneling in single-spin or simple (without higher-order anisotropy) two-spin systems can be assessed by other methods, such as the effective potential method or numerical methods involving direct diagonalization of the Hamiltonian, but the method presented here is applicable to systems with an arbitrary number of spins. The method presented here, makes it possible to estimate the onset temperature for tunneling in a magnetic system described by spin vectors as long as the second derivatives of the energy with respect to the angles describing the orientation of vectors can be evaluated at the saddle point on the energy surface.
0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7
temperature, [K]
Fig. 4. Eigenvalues of the 71 matrix block of the dimer (eq. (19)) when the exchange parameter is chosen to be J = 0.6 K. Solid line: real part. Dashed line: imaginary part. The eigenvalues depend on temperature in a similar way as is shown for the monomer in Fig. 1. A pair of eigenvalues becomes real and remains positive, while a second pair of eigenvalues changes sign at Tc from negative to positive as temperature is lowered
The most peculiar characteristic of magnetic systems is that the onset temperature for tunneling depends not only on second derivatives along the unstable mode at the first order saddle point, as in particle systems, but on the curvature of the energy surface at the first order saddle point. This is due to the inseparability of the spin Hamiltonian.
Acknowledgements
This work was supported by the Icelandic Research Fund and the Academy of Finland (grant 278260) and partially financially supported by the Government of the Russian Federation (grant 074-U01).
References
[1] Brown, Jr. W.F. Thermal Fluctuations of Fine Ferromagnetic Particles. IEEE Transactions on Magnetics, 1979, MAG-15, P. 1196-1208.
[2] Braun H.-B. Kramers's rate theory, broken symmetries and magnetization reversal. J. Appl. Physics, 1994, 76, P. 6310-6315.
[3] Chudnovsky E.M., Tejada J. Macroscopic Quantum Tunneling of the Magnetic Moment. Cambridge University Press, New York.: 1998.
[4] Bokacheva L., Kent A.D. and Walters M.A. Crossover between Thermally Assisted and Pure Quantum Tunneling in Molecular Magnet Mni2-Acetate. Phys. Rev. Lett., 2000, 85, P. 4803-4806.
[5] Mertes K.M., Zhong Y., et al. Abrupt crossover between thermally activated relaxation and quantum tunneling in a molecular magnet. Europhys. Lett., 2001, 55, P. 874-879.
[6] Wernsdorfer W., Murugesu M. and Christou G. Resonant Tunneling in Truly Axial Symmetry Mni2 Single-Molecule Magnets: Sharp Crossover between Thermally Assisted and Pure Quantum Tunneling. Phys. Rev Lett., 2006, 96, 057208 (4 pp.)
[7] Loth S., Baumann S., et al. Bistability in atomic-scale antiferromagnets. Science, 2012, 335, P. 196-199.
[8] Chudnovsky E.M. Phase transitions in the problem of the decay of a metastable state. Phys. Rev. A, 1992, 46, P. 8011-8014.
[9] Chudnovsky E. M., Garanin D. A. First- and Second-Order Transitions between Quantum and Classical Regimes for the Escape Rate of a Spin System. Phys. Rev. Lett., 1997, 79, P. 4469-4472.
[10] Park C.-S., Yoo S.-K. and Yoon D.-H. Quantum-classical crossover of the escape rate in a biaxial spin system with an arbitrarily directed magnetic field. Phys. Rev. B, 2000, 61, P. 11618-11624.
[11] Kim G.-H. and Chudnovsky E.M. Escape-rate crossover between quantum and classical regimes in molecular magnets: A diagonalization approach. Europhys. Lett., 2000, 52, P. 681-687.
[12] Garanin D.A. and Chudnovsky E.M. Quantum statistical metastability for a finite spin. Phys. Rev. B, 2000, 63, 024418 (7 pp.)
[13] Owerre S. A. and Paranjape M. B. Phase transition between quantum and classical regimes for the escape rate of dimeric molecular nanomagnets in a staggered magnetic field. Phys. Lett. A, 2014, 378, P. 1407-1412.
[14] Ulyanov V.V. and Zaslavskii O.B. New methods in the theory of quantum spin systems. Phys. Rep., 1992, 216, P. 179-251.
[15] Owerre S. A., Paranjape M. B. Macroscopic quantum tunneling and quantum-classical phase transitions of the escape rate in large spin systems. Phys. Rep., 2015, 546, P. 1-60.
[16] Wernsdorfer W. and Sessoli R. Quantum Phase Interference and Parity Effects in Magnetic Molecular Clusters. Science, 1999, 284, P. 133-135.
[17] Kang D.H. and Kim G.-H. Theoretical study of abrupt or gradual crossover of the escape rate in single-molecule magnets. Phys. Rev. B, 2006, 74, 184418 (5 pp.)
[18] Kim G.-H., Kang D.H. and Shin M.C. Quantum-classical crossover of the escape rate in the biaxial nanomagnets with a higher order symmetry. Eur. Phys. J. B, 2011, 83, P. 63-67.
[19] Vlasov S., Bessarab P.F., Uzdin V.M. and Jonsson H. Classical to quantum mechanical tunneling mechanism crossover in thermal transitions between magnetic states. Faraday Discuss., 2016, 195, P. 93-109.
[20] Bessarab P.F., Uzdin V.M. and Jonsson H. Method for finding mechanism and activation energy of magnetic transitions, applied to skyrmion and antivortex annihilation. Comp. Phys. Commun., 2015, 196, P. 335-347.
[21] Bessarab P.F., Uzdin V.M. and Jonsson H. Harmonic Transition State Theory of Thermal Spin Transitions. Phys. Rev. B, 2012, 85, 184409 (4 pp.)
[22] Bessarab P.F., Uzdin V.M. and Jonsson H. Potential Energy Surfaces and Rates of Spin Transitions. Z. Phys. Chem., 2013, 227, P. 15431557
[23] Bessarab P.F., Uzdin V.M. and Jonsson H. Size and Shape Dependence of Thermal Spin Transitions in Nanoislands. Phys. Rev. Lett., 2013, 110, 020604 (5 pp.)
[24] Moskalenko M. , Bessarab P.F., Uzdin V.M. and Jonsson H. Qualitative Insight and Quantitative Analysis of the Effect of Temperature on the Coercivity of a Magnetic System. AIP Advances, 2016, 6, 025213 (8 pp.)
[25] Lobanov I., Jonsson H. and Uzdin V. M. Mechanism and activation energy of magnetic skyrmion annihilation obtained from minimum energy path calculations. Phys. Rev. B, 2016, 94, 174418 (7 pp.)
[26] Bessarab P.F. Comment on "Path to collapse for an isolated Neel skyrmion". Phys. Rev. B, 95, 2017, 136401 (2 pp.)
[27] Ivanov A. , Bessarab P.F., Uzdin V.M. and Jonsson H. Magnetic exchange force microscopy: Theoretical analysis of induced magnetization reversals. Nanoscale, 2017, Accepted Manuscript, DOI: 10.1039/C7NR04036A
[28] Smit J., Beljers H. Ferromagnetic resonance absorption in BaFei2Oig, a highly anisotropic crystal. Philips Res. Rep., 1955,10, P. 113-130.
[29] Suhl H. Ferromagnetic resonance in nickel ferrite between one and two kilomegacycles. Phys. Rev., 1955, 97, P. 555-557.
[30] Miller W.H. Semiclassical limit of quantum mechanical transition state theory for nonseparable systems. J. Chem. Phys., 1975, 62, P. 1899-1906.
[31] Callen C., Coleman S. Fate of the false vacuum. II. First quantum corrections. Phys. Rev. D, 1977, 16, P. 1762-1768
[32] Benderskii V. A., Goldanskii V. I. and Makarov D. E. Quantum dynamics in low-temperature chemistry. Phys. Rep., 1993, 233, P. 195-339.
[33] Mills G., Schenter G. K., Makarov D. E. and Jonsson H. Generalized path integral based quantum transition state theory. Chem. Phys. Letters, 1997, 278, P. 91-96.
[34] Mills G., Schenter G. K., Makarov D. E. and Jonsson H. RAW Quantum transition state theory. proceeding of "Classical and Quantum Dynamics in Condensed Phase Simulations", Villa Marigola, Lerici, Italy, 1997, 405 pp.
[35] KIauder J. Path integrals and stationary-phase approximations. Phys. Rev. D, 1979, 19, P. 2349-2356.
[36] Kochetov E. SU(2) coherentstate path integral. J. Math. Phys., 1995, 36, P. 4667-4679.
[37] Fradkin E. Field Theories of Condensed Matter Physics. Cambridge University Press, New York.: 2013, 190 pp.
[38] Aubin S. M. J. et. al. Resonant Magnetization Tunneling in the Trigonal Pyramidal MnIVMn^11 Complex [Mn4O3Cl(O2CCH3)3(dbm)3]. J. Am. Chem. Sac., 1998, 120, P. 4991-5004.
[39] Edwards R. S., Hill S., et. al. A comparative high frequency EPR study of monomeric and dimeric Mn4 single-molecule magnets, Polyhedron, 2003, 22, P. 1911-1916.
[40] Wernsdorfer W., Aliaga-Alcalde N., Hendrickson D. N., Christou G. Exchange-biased quantum tunnelling in a supramolecular dimer of single-molecule magnets. Nature, 2002, 416, P. 406-409.
