Abstract
The classical spin system consisting of three spins with Heisenberg interaction is an example of a completely integrable mechanical system. In this paper, we explicitly calculate thermodynamic quantities such as density of states, specific heat, susceptibility and spin autocorrelation functions. These calculations are performed (semi-)analytically and shown to agree with corresponding Monte Carlo simulations. It is shown that the thermodynamic functions depend qualitatively on the character of the system in terms of its frustration, especially w. r. t. their low temperature limit. For the long-time autocorrelation function, we find, for certain values of the coupling constants, a decay to constant values in the form of an 1/t damped harmonic oscillation and propose a theoretical explanation.
-
Author contributions: All the authors have accepted responsibility for the entire content of this submitted manuscript and approved submission.
-
Research funding: None declared.
-
Conflict of interest statement: The authors declare no conflicts of interest regarding this article.
Appendix A: Details of the explicit time evolution
Instead of u we will use the variable x given by
where the constants x0 and g will be determined later such that the Weierstrass differential equation [[30], 23.3.10] is obtained. Also v and w can be linearly expressed in terms of the variable x in the form:
We consider the time derivative of x:
By substituting (A2), (A3) and (A1), the square of (A5) can be written as a 3rd order polynomial Π(x). g and x0 will be chosen such that the cubic term of Π(x) reads 4x3 and the quadratic term of Π(x) vanishes and hence
This is achieved by setting
and
The explicit form of the coefficients g2 and g3 is more complicated:
and
where
For statistical considerations parts of the phase space with zero measure can be neglected (but note that the dos may diverge for states with aperiodic motion according to (48)). Hence we can restrict ourselves to the “generic case” where certain exceptions are excluded, see [11]. In this generic case the polynomial Π(x) will have three real simple roots x1 < x2 < x3 satisfying x1 + x2 + x3 = 0 and Π(x) > 0 for x1 < x < x2. The explicit form of the roots is known but of overwhelming complexity if expressed in terms of the physical parameters ɛ, σ, J1, J2, J3.
It follows [11] that in the generic case (A5) has the solution
with the above-mentioned parameters g2, g3 and the imaginary parameter t0 can be expressed through an elliptical integral:
see [[30], 23.6.34-35] and [[40], 17.4.61 ff]. Moreover, this solution will be
Hence, for given coupling constants J1, J2, J3, the period
Appendix B: Distribution of random variables with a saddle point
We will prove the following proposition which is tailored to its application in Section 7.2.2 and not formulated as general as possible:
Proposition 1
Let
for all
Then ρ
f
has a logarithmic singuarity at
Proof
Without loss of generality we may assume
where H(ϵ) is the hyperbolic region
see Figure 16. By assumption, φ(0, 0) > 0 and, since φ is continuous, we may choose R > 0 so small such that
This implies
where
see Figure 16. If u↦ρ f (u) would be continuous in a neighborhood of u = 0 then it would follow that
which is a contradiction due to the divergence of −log ϵ. Hence ρ f (0) is divergent and the singularity is, at least, of logarithmic order.
The singularity is exactly of logarithmic order since the contribution to ρ
f
(0) from other possible zeroes

Plot of the hyperbolic domain H(ϵ) given by |x y| ≤ ϵ and |x|, |y| ≤ R. The area of the intersection of H(ϵ) with the positive quadrant is given by
Appendix C: Asymptotic expansions of Fourier integrals involving logarithmic singularities
We consider the case of a distribution function ρ(ω) with a logarithmic singularity at ω = ω0 and will investigate the decay of the corresponding Fourier transform
where ϕ(ω) is N times continuously differentiable for γ1 < ω < γ2 and γ1 < ω0 < γ2. For our purposes we may assume that ϕ(ω) is a real function. Then we consider the Fourier integral
and the asymptotic expansion of
We will utilize this solution to obtain the asymptotic expansion for the two-sided Fourier integral (C2). For this purpose, we will quote the corresponding theorem 4 of [41] in full detail, with slight modifications according to our notation.
Proposition 2
(Erdélyi). Under the preceding assumptions on ϕ(ω) we have
for t → +∞, where ψ(z) denotes the logarithmic derivative of Γ(z).
Let us denote the complementary integral of (C3) by
such that
where we have used that ϕ(ω) is real. This form of
we thus obtain from (C4) and the replacement ω0↦ −ω0
This entails
and, finally,
We note that the terms containing log t cancel and the leading term corresponding to n = 0 in the asymptotic expansion (C13) is proportional to 1/t.
References
[1] E. H. Lieb, “The classical limit of quantum spin systems,” Commun. Math. Phys., vol. 31, pp. 327–340, 1973. https://doi.org/10.1007/bf01646493.Suche in Google Scholar
[2] J. Fröhlich, A. Knowles, and E. Lenzmann, “Semi-classical dynamics in quantum spin systems,” Lett. Math. Phys., vol. 82, pp. 275–296, 2007. https://doi.org/10.1007/s11005-007-0202-y.Suche in Google Scholar
[3] P. Pyykkö, “Magically magnetic gadolinium,” Nat. Chem., vol. 7, p. 680, 2015. https://doi.org/10.1038/nchem.2287.Suche in Google Scholar PubMed
[4] K. B. Ghiassi, M. M. Olmstead, and A. L. Balch, “Gadolinium-containing endohedral fullerenes: structures and function as magnetic resonance imaging (MRI) agents,” Dalton Trans., vol. 43, pp. 7346–7358, 2014. https://doi.org/10.1039/c3dt53517g.Suche in Google Scholar PubMed
[5] L. Qin, G.-J. Zhou, Y.-Z. Yu, et al.., “Topological self-assembly of highly-symmetric lanthanide clusters: a magnetic study of exchange- coupling “fingerprints” in giant gadolinium(III) cages,” J. Am. Chem. Soc., vol. 45, pp. 16405–16411, 2017. https://doi.org/10.1021/jacs.7b09996.Suche in Google Scholar PubMed
[6] L. Qin, H.-L. Zhang, Y.-Q. Zhai, et al.., “A giant spin molecule with ninety-six parallel unpaired electrons,” iScience, vol. 24, no. 4, p. 102350, 2021. https://doi.org/10.1016/j.isci.2021.102350.Suche in Google Scholar PubMed PubMed Central
[7] H.-J. Schmidt, A. Lohmann, and J. Richter, “Eighth-order high-temperature expansion for general Heisenberg Hamiltonians,” Phys. Rev. B, vol. 84, p. 104443, 2011. https://doi.org/10.1103/physrevb.84.104443.Suche in Google Scholar
[8] G. M. Wysin, Magnetic Excitations and Geometric Confinement, Bristol, England, IOP Publishing, 2015.10.1088/978-0-7503-1074-1Suche in Google Scholar
[9] A. K. Boudalis, G. Rogez, and P. Turek, “Determination of the distributions of the spin-Hamiltonian parameters in spin triangles: a combined magnetic susceptometry and electron paramagnetic resonance spectroscopic study of the highly symmetric [Cr3O(PhCOO)6(py)3](ClO4)·0.5py,” Inorg. Chem., vol. 57, pp. 13259–13269, 2018. https://doi.org/10.1021/acs.inorgchem.8b01764.Suche in Google Scholar PubMed
[10] V. I. Arnol’d, Mathematical Methods of Classical Mechanics, Berlin, Springer, 1978.10.1007/978-1-4757-1693-1Suche in Google Scholar
[11] H.-J. Schmidt, The Classical Spin Triangle as an Integrable System. Preprint math-ph:2109.04841, 2021.Suche in Google Scholar
[12] H.-J. Schmidt, “The general spin triangle,” Int. J. Mod. Phys. B, vol. 278, no. 16, p. 1350064, 2013. https://doi.org/10.1142/s0217979213500641.Suche in Google Scholar
[13] H.-J. Schmidt, C. Schröder, E. Hägele, and M. Luban, “Dynamics and thermodynamics of a pair of interacting dipoles,” J. Phys. Math. Theor., vol. 48, p. 185002, 2015. https://doi.org/10.1088/1751-8113/48/18/185002.Suche in Google Scholar
[14] J. H. Luscombe, M. Luban, and F. Borsa, “Classical Heisenberg model of magnetic molecular ring clusters: accurate approximants for correlation functions and susceptibility,” J. Chem. Phys., vol. 108, no. 17, pp. 7266–7273, 1998. https://doi.org/10.1063/1.476144.Suche in Google Scholar
[15] O. Ciftja, M. Luban, M. Auslender, and J. H. Luscombe, “Equation of state and spin-correlation functions of ultrasmall classical Heisenberg magnets,” Phys. Rev. B, vol. 60, no. 14, pp. 10122–10133, 1999. https://doi.org/10.1103/physrevb.60.10122.Suche in Google Scholar
[16] D. Mentrup, J. Schnack, and M. Luban, “Spin dynamics of quantum and classical Heisenberg dimers,” Physica A, vol. 272, pp. 153–161, 1999. https://doi.org/10.1016/s0378-4371(99)00239-3.Suche in Google Scholar
[17] D. Mentrup, H.-J. Schmidt, J. Schnack, and M. Luban, “Transition from quantum to classical Heisenberg trimers: thermodynamics and time correlation functions,” Physica A, vol. 278, pp. 214–221, 2000. https://doi.org/10.1016/s0378-4371(99)00571-3.Suche in Google Scholar
[18] O. Ciftja, “Spin correlation functions of some frustrated ultra-small classical Heisenberg clusters,” Physica A, vol. 268, pp. 541–557, 2000. https://doi.org/10.1016/s0378-4371(00)00326-5.Suche in Google Scholar
[19] M. Ameduri and R. A. Klemm, “Time correlation functions of three classical Heisenberg spins on an isosceles triangle and on a chain,” Phys. Rev. B, vol. 66, p. 224404, 2002. https://doi.org/10.1103/physrevb.66.224404.Suche in Google Scholar
[20] O. Ciftja, “Spin dynamics of an ultra-small nanoscale molecular magnet,” Nanoscale Res. Lett., vol. 2, pp. 168–174, 2007. https://doi.org/10.1007/s11671-007-9049-5.Suche in Google Scholar
[21] O. Ciftja, “Exact classical spin dynamics of high spin nanoscale molecular magnetic clusters,” Phys. Lett. A, vol. 424, p. 127826, 2022. https://doi.org/10.1016/j.physleta.2021.127826.Suche in Google Scholar
[22] F. Borsa and M. Mali, “Experimental study of high-temperature spin dynamics in one-dimensional Heisenberg paramagnets,” J. Phys. B, vol. 8, no. 5, pp. 2215–2219, 1974. https://doi.org/10.1103/physrevb.9.2215.Suche in Google Scholar
[23] H. Takano and S. Miyashita, “Relaxation of the spin autocorrelation function in the kinetic ising model with bond dilution,” J. Phys. Soc. Jpn, vol. 58, pp. 3871–3874, 1989. https://doi.org/10.1143/jpsj.58.3871.Suche in Google Scholar
[24] A. Lascialfari, Z. H. Jang, F. Borsa, D. Gatteschi, A. Cornia, “Comparison of the spin dynamics in different types of molecular magnetic rings from 1H NMR,” J. Appl. Phys., vol. 83, no. 11, pp. 6946–6948, 1989.10.1063/1.367564Suche in Google Scholar
[25] M. L. Baker, T. Guidi, S. Carretta, et al.., “Spin dynamics of molecular nanomagnets unravelled at atomic scale by four-dimensional inelastic neutron scattering,” Nat. Phys., vol. 8, no. 12, pp. 906–911, 2012. https://doi.org/10.1038/nphys2431.Suche in Google Scholar
[26] R. A. Klemm and M. Luban, “Zero-field time correlation functions of four classical Heisenberg spins on a ring,” Phys. Rev. B, vol. 64, p. 104424, 2001. https://doi.org/10.1103/physrevb.64.104424.Suche in Google Scholar
[27] H.-J. Schmidt and M. Luban, “Classical ground states of symmetric Heisenberg spin systems,” J. Phys. Math. Gen., vol. 36, no. 23, p. 6351, 2003. https://doi.org/10.1088/0305-4470/36/23/306.Suche in Google Scholar
[28] H.-J. Schmidt, Theory of Ground States For Classical Heisenberg Spin Systems I. Preprint cond-mat:1701.02489v2, 2017.Suche in Google Scholar
[29] H.-J. Schmidt, Theory of Ground States For Classical Heisenberg Spin Systems III. Preprint cond-mat:1707.06512, 2017.Suche in Google Scholar
[30] F. W. J. Olver, A. B. Olde Daalhuis, D. W. Lozier, et al.., Eds. NIST Digital Library of Mathematical Functions. 2021. Available at: http://dlmf.nist.gov/. Release 1.1.1 of 2021-03-15.Suche in Google Scholar
[31] J. Schnack, “Effects of frustration on magnetic molecules: a survey from Olivier Kahn until today,” Dalton Trans., vol. 39, pp. 4677–4686, 2010. https://doi.org/10.1039/b925358k.Suche in Google Scholar PubMed
[32] F. Wang and D. P. Landau, “Efficient, multiple-range random walk algorithm to calculate the density of states,” Phys. Rev. Lett., vol. 86, pp. 2050–2053, 2001. https://doi.org/10.1103/physrevlett.86.2050.Suche in Google Scholar
[33] T. Nakanishi and S. Yamamoto, “Intrinsic double-peak structure of the specific heat in low-dimensional quantum ferrimagnets,” Phys. Rev. B, vol. 65, p. 214418, 2002. https://doi.org/10.1103/physrevb.65.214418.Suche in Google Scholar
[34] M. de Souza, R. Paupitz, A. Seridonio, and R. E. Lagos, “Specific heat anomalies in solids described by a multilevel model,” Braz. J. Phys., vol. 46, pp. 206–212, 2016. https://doi.org/10.1007/s13538-016-0404-9.Suche in Google Scholar
[35] K. Karl’ová, J. Strečka, and T. Madaras, “The Schottky-type specific heat as an indicator of relative degeneracy between ground and first-excited states: the case study of regular Ising polyhedra,” Phys. B Condens. Matter, vol. 488, pp. 49–56, 2016.10.1016/j.physb.2016.01.033Suche in Google Scholar
[36] N. Wiener, “Generalized harmonic analysis,” Acta Math., vol. 55, pp. 117–258, 1930. https://doi.org/10.1007/bf02546511.Suche in Google Scholar
[37] A. Khintchine, “Korrelationstheorie der stationären stochastischen Prozesse,” Math. Ann., vol. 109, no. 1, pp. 604–615, 1934. https://doi.org/10.1007/bf01449156.Suche in Google Scholar
[38] M. Luban and J. H. Luscombe, “Equilibrium time correlation functions and the dynamics of fluctuations,” Am. J. Phys., vol. 67, pp. 1161–1169, 1999. https://doi.org/10.1119/1.19102.Suche in Google Scholar
[39] N. Chernov, “Decay of correlations,” Scholarpedia, vol. 3, no. 4, p. 4862, 2008. https://doi.org/10.4249/scholarpedia.4862.Suche in Google Scholar
[40] M. Abramowitz and I. A. Stegun, Eds. Handbook of Mathematical Functions, New York, Dover, 1972.Suche in Google Scholar
[41] A. Erdélyi, “Asymptotic expansions of Fourier integrals involving logarithmic singularities,” J. Soc. Ind. Appl. Math., vol. 4, no. 1, pp. 38–47, 1956. https://doi.org/10.1137/0104003.Suche in Google Scholar
© 2022 Walter de Gruyter GmbH, Berlin/Boston
Artikel in diesem Heft
- Frontmatter
- Dynamical Systems & Nonlinear Phenomena
- Qualitative behavior of a discrete predator–prey system under fear effects
- Nonlinear behaviour of ion acoustic shock waves in a two-electron temperature nonthermal complex plasma
- Gravitation & Cosmology
- Study of baryogenesis in the framework of Hořava–Lifshitz cosmology with Starobinsky potential
- Solid State Physics & Materials Science
- Role of graphene-oxide and reduced-graphene-oxide on the performance of lead-free double perovskite solar cell
- Thermodynamics & Statistical Physics
- Thermodynamics of the classical spin triangle
Artikel in diesem Heft
- Frontmatter
- Dynamical Systems & Nonlinear Phenomena
- Qualitative behavior of a discrete predator–prey system under fear effects
- Nonlinear behaviour of ion acoustic shock waves in a two-electron temperature nonthermal complex plasma
- Gravitation & Cosmology
- Study of baryogenesis in the framework of Hořava–Lifshitz cosmology with Starobinsky potential
- Solid State Physics & Materials Science
- Role of graphene-oxide and reduced-graphene-oxide on the performance of lead-free double perovskite solar cell
- Thermodynamics & Statistical Physics
- Thermodynamics of the classical spin triangle