# Quadrupole Excitation in Tunnel Splitting Oscillation in NanoparticleMn12

Authors: Yousef Yousefi; Khikmat Kh. Muminov
Journal: Advances in Condensed Matter Physics
 (2012)

---

## Abstract

We analyze the interference between tunneling paths that occur
for a spin system with special Hamiltonian both for dipole and quadrupole
excitations. Using an instanton approach, we find that as the strength
of the second-order transverse anisotropy is increased, the tunnel splitting for both excitations is modulated, with zeros occurring periodically and the number of quenching points for quadrupole excitation
decreasing. This effect results from the interference of four tunneling
paths connecting easy-axis spin orientations and occurs in the absence
of any magnetic field.


---

## Body

## 1. Introduction

Berry phase effects play an important role in spin dynamics. Tunneling of a spin (or magnetic particle) between degenerate orientations can be calculated by such Berry phase effects via the interference between tunneling paths, with quenching of the tunnel splitting occurring when tunneling paths destructively interfere [1, 2].The first conclusive experimental evidence of tunneling in molecular magnet Mn12O12(CH3COO)16(H2O)4 called Mn12 was provided by Friedman [3]. He applies a magnetic field along the easyaxis of the crystal at a variety of temperatures and found a curve punctuated by sharp steps at regular intervals in the parameter λ. Also this observation that is confirmed by Hernandez et al. [4] and Thomas et al. [5] shows a series of steps in the hysteresis loops in Mn12 in low temperature. The steps observed in the hysteresis loops at nearly equal intervals of magnetic field are due to enhanced relaxation of the magnetization at the resonant fields when levels on opposite side of the anisotropy barrier coincide in energy (see Figure 1).The first curve was obtained from numerical calculation of the tunnel splitting for Mn12 and the second curve is the magnetization plot of this magnetic molecule in low temperature.
(a)(b)The common feature of molecular magnets responsible for their interesting behavior is the strong spin-spin coupling between metallic ions in the core of the molecule. In this magnet molecule, there are 4 Mn4+ atoms (spin 3/2) surrounded by 8 Mn3+ atoms (spin 2), and they retain a relative orientation such as to give the whole molecule a spin of 10 (= 8(2)−4(3/2)) [3].The Mn12 single-molecule magnet has a fourfold transverse magnetic anisotropy and displays resonant tunneling between two easy-axis orientations [3]. Here we consider the interference that occurs between tunneling paths in such a system and how that interference can be modified by the presence of a second-order transverse anisotropy perturbation and multipole excitations. Such a perturbation could potentially be induced by the application of uniaxial pressure to a sample of Mn12. We find that for both excitations the tunnel splitting is periodically quenched as a function of the strength of the perturbation. This interference effect takes place in the absence of any magnetic field. In the present study, the fourth-order anisotropy term is the primary transverse anisotropy in the problem, producing four interfering paths, but the interference is modulated by the strength of the second-order anisotropy, leading to periodic quenches as that term is varied. In particular, we consider a spin governed by the Hamiltonian(1)Ĥ=Ĵ2-Ĵz2+k(Ĵ+4+Ĵ-4)+λ(Ĵx2-Ĵy2),
where the Ĵs are standard spin operators. We restrict the values of k and λ such that the z-axis is the spin easy axis. In the above Hamiltonian, the leading term -Ĵz2 established the double well potential and the easy axis and the term fourth-order transverse anisotropy k(Ĵ+4+Ĵ-4) are responsible for 4-fold rotational symmetries. The third term is a second-order transverse anisotropy that is present in many low-symmetry SMMs.In this paper, firstly we calculated periodically depended of spin tunneling for Mn12 in SU(2) group, in other word, we considered only dipole excitation in Hamiltonian. Due to the symmetry of spin operators in Hamiltonian, for being more accurate, we have to consider other multipole excitation. Then, we consider both dipole and quadrupole excitations.


## 2. Theory and Calculations

The instanton method is an efficient way of calculating tunnel splitting, both for particles [6] and for spin [7]. It is based on evaluating the path integral for a certain propagator in the steepest-descent approximation and is designed to be asymptotically correct in the semiclassical limit (J→∞ or ℏ→0). Instantons are classical paths that run between degenerate classical minima of the energy. By classical we mean that the path obeys the principle of least action and satisfies energy conservation. However, a path along which energy is conserved can not run between two minima and still have real coordinates and momenta. Hence, one must enlarge the notion of a classical path and allow the coordinates and/or momenta to become complex.In the continuous approach the tunnel splitting can be computed as the functional integral [8](2)∫exp⁡(-S(ẋ,x,t))dμ(x),
where dμ(x)=∏i=1Ndxi. In Su(2) group N=2 and in Su(3) group N=4, and N is the number of degree of freedoms, and the action is given by(3)S(ẋ,x,t)=SK+SD+SB.The first term,SK, is kinetic term which has the properties of a Berry phase, which can give rise to interference between different trajectories. The second term, SD, is dynamical term and the final term is boundary term this term depends explicitly on the boundary values of the path.When the functional action is specified, the instanton recipe for calculating the tunnel splitting is as follows. Let there be a number of instantons, that is, least action paths, labeled byk, and let the actions for these various paths be Skcl. The tunneling amplitude is given by (4)Δ=∑kDkexp⁡(-Skcl⁡).The prefactorDk results from integrating the Gaussian fluctuations around the kth instanton.For nanoparticle Mn12, the tunneling amplitude calculated and obtained the following relation [3]:(5)Δ≈4e-SRcos⁡⁡(Jπ)cos⁡⁡(JSI),
where(6)SI=Re(∫π/23π/2P(ϕ)dϕ),
and P(ϕ) is solution of classical enrgy equation. SR is the real part of the Euclidean action and the same for all instantons. The factor cos(Jπ) ensures that the tunnel splitting is zero for half-integer spin.The classical energy is (thex-axis is chosen as the azimuth because the-path will not encounter it and ϕ is measured from the z axis)(7)Ecl⁡(θ,ϕ)=J2(1-sin⁡2⁡θsin⁡2ϕ)+λJ2(cos⁡2⁡θ-sin⁡2⁡θcos⁡2⁡ϕ)+2kJ4(cos⁡4⁡θ+sin⁡4⁡θcos⁡4⁡ϕ-6cos⁡2⁡θsin⁡2⁡θcos⁡2⁡ϕ).The minimum of energy obtained at(8)θ=π2ϕ=π2,3π2.The minima of enregy are(9)Ecl⁡(π2,π2)=0Ecl⁡(π2,3π2)=0.Making the substitutionP=cosθ, the instanton must satisfy the constraint(10)J2(1-(1-P2)sin⁡2⁡ϕ)+λJ2(P2-(1-P2)cos⁡2⁡ϕ)+2kJ4(P4+(1-P2)2cos⁡4⁡ϕ-6P2(1-P2)cos⁡2⁡ϕ)=0.The four solutions for the above equation are:(11)P(ϕ)=±A(ϕ)±iB(ϕ)C(ϕ),
where(12)A(ϕ)=2acos⁡2⁡ϕ(7+cos⁡⁡2ϕ)-λ(1+cos⁡2⁡)-sin⁡2⁡ϕB(ϕ)=acos⁡2⁡ϕ(1+k-λ+kcos⁡2ϕ)×(35+28cos⁡⁡2ϕ+cos⁡⁡4ϕ)-(λ(1+cos⁡2⁡ϕ)-2acos⁡2⁡ϕ(7+cos⁡2ϕ)+sin⁡2⁡ϕ)2C(ϕ)=4a(cos⁡4⁡ϕ+6cos⁡2⁡ϕ+1),
and a=J2k.Solutions (11), despite conserving energy, fail to meet the boundary conditions P(π/2)=P(3π/2)=0. These solutions gives rise to what Garg called boundary jump instantons [9]. As a final note, we consider only the following value of P(ϕ):(13)P(ϕ)=A(ϕ)+iB(ϕ)C(ϕ).Substituting above relation in (6) and using from (5), splitting of level energy calculated. Up this section, we have just considered the dipole excitation of spin system. Because we will analyze the periodical dependence of spin tunneling, by relation (5), plotted function Cos(JSI) versus λ.If quadrupole excitation is considered, the classical energy is changed in the following form:(14)Ecl⁡(θ,ϕ)=J2(1-sin⁡2⁡θcos⁡2⁡2gsin⁡2⁡ϕ)+λJ2cos⁡2⁡2g(cos⁡2⁡θ-sin⁡2⁡θcos⁡2⁡ϕ)+2kJ4cos⁡4⁡2g(cos⁡4⁡θ+sin⁡4⁡θcos⁡2⁡ϕ-6cos⁡2θsin⁡2⁡θcos⁡2⁡ϕ).The minimum of energy obtained at(15)θ=π2,g=π4,ϕ=π2,3π2.Also minima of enrgy are(16)Ecl⁡(π2,π2,π4)=0Ecl⁡(π2,3π2,π4)=0.Making the substitutionP=Cosθ, the instanton must satisfy the constraint:(17)J2(1-(1-P2)cos⁡2⁡gsin⁡2⁡ϕ)+λJ2cos⁡2⁡2g(P2-(1-P2)cos⁡2⁡ϕ)+2kJ4cos⁡4⁡2g(P4+(1-P2)2cos⁡4⁡ϕ-6P2(1-P2)cos⁡2⁡ϕ)=0.Similar to previous section, solution that has the following form is considered:(18)P(ϕ,g)=A(ϕ,g)+iB(ϕ,g)C(ϕ,g),
where (19)A(ϕ,g)=2acos⁡2⁡ϕsec2⁡2g(7+cos⁡⁡2ϕ)-12(1+3λ+(λ-1)cos⁡⁡2ϕ)sec4⁡2gB(ϕ,g)=8a(1+6cos⁡2⁡ϕ+cos⁡4⁡ϕ)×(1+2kcos⁡4⁡2gcos⁡4ϕ-cos⁡22g(λcos⁡2ϕ+sin⁡2ϕ))-(λ+sin⁡2⁡ϕ+cos⁡2⁡ϕ(λ-2acos⁡2⁡2g(7+cos⁡⁡2ϕ)))2C(ϕ,g)=4a(1+6cos⁡2⁡ϕ+cos⁡4ϕ)cos⁡2⁡2g.Substituting relation (19) in (6) and using (5), splitting of energy level was calculated. In this section, we consider both dipole and quadrupole excitations of spin system. Like the previous section, plotted function cos(JSI) versus λ. In this plot, a=10-4 and k=10-6 are considered.As seen from Figures2 and 3, when we added the quadrapole excitation to Hamilton, the number of quenching point reduced from 5 to 4, also by using Figure 1 we see that the number of steps is 4; that is equal to the number of quenching points and then these steps in magnetization graphic are quenching point is in spin tunneling phenomena.Figure 2
Red curve isRe(P)/π and blue curve is cos(JSI) versus λ obtained by using (11). In this graph, only dipole excitation was considered.Figure 3
Red curve isRe(P)/π and blue curve is cos(JSI) versus λ obtained using (16). In this graph both dipole and quadrupole excitations are considered.As mentioned in the introduction, there are step changes in the experimental hysteresis loop of this crystal that have been explained considering only dipolar interaction between molecules. For a more accurate description of the experimental results, we found that we have to consider both dipolar and quadrupole excitations.


## 3. Discussion

In this paper, we consider the spin tunneling phenomena to parameterλ in both groups SU(2) and SU(3). Graph of function Re(P)versus λ in SU(2) group is a straight line started from 5 and with increasing λ up to 0.04, this straight line decreases and in point λ=0.04 the value is zero and remains in this value. But above function graph in SU(3) group, started less than 4, is not straight line and in point λ=0.04 the value is zero and remains in this value. The number of quenching points in SU(2) group (dipole excitation) is 5, these points decreases to 4, when quadrupole excitation is considered.



---

*Source: 530765-2012-05-16.xml*