Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2412.15432

Long-range exchange coupling in a magnetic multilayer system

Leonardo Cordeiro de Souza, Rafael Alves Dumer, Maurício Godoy

Published Dec 19, 2024• 0 citations
arXivHTMLPDF

Abstract

In this study, we explored the magnetic coupling in a multilayer system consisting of thin layers separated by a distance $d$. We have employed Monte Carlo (MC) simulations to calculate the thermodynamic quantities such as the magnetization per spin $m_{L}^μ$, magnetic susceptibility $χ_{L}^μ$, and the reduced fourth-order Binder cumulant $U_{L}^μ$ as a function of temperature $T$ and for several values of lattice size $L$. These quantities were obtained for each layer ($μ=l)$, for the bulk system ($μ=b)$, as a function of interlayer distance $d$ and the parameter $α$ that defines the scale length of the exponential decay of the interaction. Furthermore, we applied the finite-size scaling theory to calculate the critical exponents. Our results reveal that the system exhibits 2D Ising exponents when the layers are sufficiently separated. On the other hand, in the compact limit, where $d$ equals the distance between two adjacent sites within the same layer, our bulk results show that the system exhibits 3D Ising critical exponents, provided that the number of layers $L_{z}$ increases proportionally to the layer sizes. Even with the separation increasing up to $d=2.0$, the layers are so correlated that the set of critical exponents retains the values of 3D critical exponents. However, when the number of layers $L_{z}$ remains fixed, even in the compact limit with periodic boundary conditions, only the exponent $β$ aligns closely with the predicted literature values, on the other hand, the other exponents show significant deviations.

DOI
10.48550/arxiv.2412.15432

Related papers

Finite-size effects in layered magnetic systems

Dragi Karevski, Malte Henkel

Thermal and magnetic effects in a system consisting of thin layers of coupled Ising spins with S=1/2 and S=1 are considered. The specific heat and the correlation length display maxima at two different temperatures. It is discussed in what sense these maxima can be interpreted as a finite-size rounding of a thermodynamic singularity associated with a phase transition. The connection with ordinary, extraordinary, and special surface phase transitions is made. In two dimensions (2D), the surface critical exponents are calculated from conformal invariance. The bulk and surface finite-size scaling of the order-parameter profiles at the transition points is discussed. In 2D, an exact scaling function for the profiles is suggested through conformal invariance arguments for the (extra)ordinary transition.

arXiv:cond-mat/9610079•Physical review. B, Condensed matter•Mar 1, 1997•22 citations

Is the hyperscaling relation violated below the upper critical dimension in some particular cases?

H. T. Diep, Van Thanh Ngo

In this review, we show our results with new interpretation on the critical exponents of thin films obtained by high-performance multi-histogram Monte Carlo simulations. The film thickness $N_z$ consists of a few layers up to a dozen of layers in the $z$ direction. The free boundary condition is applied in this direction while in the $xy$ plane periodic boundary conditions are used. Large $xy$ plane sizes are used for finite-size scaling. The Ising model is studied with nearest-neighbor (NN) interaction. When $N_z=1$, namely the two-dimensional (2D) system, we find the critical exponents given by the renormalization group. While, for $N_z>1$, the critical exponents calculated with the high-precision multi-histogram technique show that they deviate slightly but systematically from the 2D values. If we use these values of critical exponents in the hyperscaling relation with $d=2$, then the hyperscaling relation is violated. However, if we use the hyperscaling relation and the critical exponents obtained for $N_z>1$ to calculate the dimension of the system, we find the system dimension slightly larger than 2. This can be viewed as an "effective" dimension. More discussion is given in the paper. We also show the cross-over between the first- and second-order transition while varying the film thickness in an antiferromagnetic FCC Ising frustrated thin film. In addition, we will show evidence that when a 2D system has two order parameters of different symmetries with a single transition, the critical exponents are new, suggesting a universality class of coupled two-symmetry breakings. In this case, the 2D hyperscaling does not hold. Another case is the 3D Ising model coupled to the lattice vibration: the critical exponents deviate from the 3D Ising ones, the results suggest the violation of the hyperscaling.

arXiv:2507.02159•Jul 2, 2025•0 citations

Transfer matrix and Monte Carlo tests of critical exponents in Ising model

J. Kaupužs

The corrections to finite-size scaling in the critical two-point correlation function G(r) of 2D Ising model on a square lattice have been studied numerically by means of exact transfer-matrix algorithms. The systems have been considered, including up to 800 spins. The calculation of G(r) at a distance r equal to the half of the system size L shows the existence of an amplitude correction proportional to 1/L^2. A nontrivial correction ~1/L^0.25 of a very small magnitude also has been detected, as it can be expected from our recently developed GFD (grouping of Feynman diagrams) theory. Monte Carlo simulations of the squared magnetization of 3D Ising model have been performed by Wolff's algorithm in the range of the reduced temperatures t =< 0.000086 and system sizes L =< 410. The effective critical exponent beta_eff(t) tends to increase above the currently accepted numerical values. The critical coupling K_c=0.22165386(51) has been extracted from the Binder cumulant data within 96 =< L =< 384. The critical exponent 1/nu, estimated from the finite-size scaling of the derivatives of the Binder cumulant, tends to decrease slightly below the RG value 1.587 for the largest system sizes. The finite-size scaling of accurately simulated maximal values of the specific heat C_v in 3D Ising model confirms a logarithmic rather than power-like critical singularity of C_v.

arXiv:cond-mat/0405197•May 10, 2004•1 citation

Magnetic properties of Ising thin films with cubic lattices

Yongyut Laosiritaworn, J. Poulter, Julie B. Staunton

We have used Monte Carlo simulations and mean-field analysis to observe the magnetic behavior of Ising thin films with cubic lattice structures as a function of temperature and thickness, especially in the critical region. Magnetization and magnetic susceptibility, including layer variation, are investigated. We find that the magnetic behavior changes from two-dimensional to three-dimensional character with increasing film thickness. Both the crossover of the critical temperature from a two-dimensional to a bulk value and the shift exponent are observed. Nevertheless, with support from a scaling function, the simulations show that the effective critical exponents for films with large enough layer extents only vary a little from their two-dimensional values. This, in particular, provides an indication of two-dimensional universality in the thin films.

arXiv:cond-mat/0307291•Physical Review B•Sep 22, 2004•53 citations

Surface effects on magnetic phase transitions

Kurt Binder, Pierre C. Hohenberg

The surface scaling theory previously presented by the authors is developed further, and derived heuristically from a cluster model. Monte Carlo calculations are carried out to obtain the spatial and temperature dependence of the magnetization in Ising and Heisenberg systems with free surfaces. The exponent ${\ensuremath{\beta}}_{1}$ of the (surface) layer magnetization is shown to agree with the scaling value (${\ensuremath{\beta}}_{1}\ensuremath{\approx}\frac{2}{3}$) previously derived. In the Heisenberg system, the results at low temperature agree with a spin-wave calculation by Mills and Maradudin. Ising models with modified exchange ${J}_{s}=J(1+\ensuremath{\Delta})\ensuremath{\ne}J$ on the surface are considered, both in mean-field theory and by means of high-temperature-series expansions. The critical value ${\ensuremath{\Delta}}_{c}$ for surface ordering is found from the series to be 0.6, compared to the mean-field value of 0.25. For $\ensuremath{\Delta}>{\ensuremath{\Delta}}_{c}$ there is a temperature region in which the surface behaves like a bulk two-dimensional Ising model near its phase transition. The critical exponents experience a crossover at $\ensuremath{\Delta}={\ensuremath{\Delta}}_{c}$, which is reflected in poorly behaved series, and effective exponents differing from the true ones for $\ensuremath{\Delta}\ensuremath{\lesssim}{\ensuremath{\Delta}}_{c}$. In the case of weakened surface exchange ($0<{J}_{s}<J$), the layer magnetization is shown to fit a linear temperature dependence over a large temperature range below ${T}_{c}$, thus providing a possible explanation for previous experiments. For sufficiently strong negative ${J}_{s}$, mean-field theory predicts that the surface will order antiferromagnetically while the bulk is ferromagnetic.

Physical review. B, Solid state•Mar 1, 1974•566 citations

Glassy phases in a layered ising model with random interlayer exchange

P. N. Timonin

Thermodynamics of a layered Ising model with infinite-range ferromagnetic intralayer interaction and random nearest-neighbor interlayer coupling is considered. A detailed analysis of the model with vanishing average interlayer coupling is presented. The Gibbs free energy is found in the critical region, and the existence of many metastable states is demonstrated. Thermodynamic parameters of the system are found for periodic states. As the mean square interlayer coupling increases, the equilibrium state of the system undergoes an infinite sequence of first-order phase transitions, the number of magnetic planes and the distance between them change discontinuously, and so do both bulk magnetization and magnetic susceptibility.

Journal of Experimental and Theoretical Physics•Mar 1, 2006•0 citations
In this study, we explored the magnetic coupling in a multilayer system consisting of thin layers separated by a distance ddd. We have employed Monte Carlo (MC) simulations to calculate the thermodynamic quantities such as the magnetization per spin mLμm_{L}^μmLμ​, magnetic susceptibility χLμχ_{L}^μχLμ​, and the reduced fourth-order Binder cumulant ULμU_{L}^μULμ​ as a function of temperature TTT and for several values of lattice size LLL. These quantities were obtained for each layer (μ=l)μ=l)μ=l), for the bulk system (μ=b)μ=b)μ=b), as a function of interlayer distance ddd and the parameter ααα that defines the scale length of the exponential decay of the interaction. Furthermore, we applied the finite-size scaling theory to calculate the critical exponents. Our results reveal that the system exhibits 2D Ising exponents when the layers are sufficiently separated. On the other hand, in the compact limit, where ddd equals the distance between two adjacent sites within the same layer, our bulk results show that the system exhibits 3D Ising critical exponents, provided that the number of layers LzL_{z}Lz​ increases proportionally to the layer sizes. Even with the separation increasing up to d=2.0d=2.0d=2.0, the layers are so correlated that the set of critical exponents retains the values of 3D critical exponents. However, when the number of layers LzL_{z}Lz​ remains fixed, even in the compact limit with periodic boundary conditions, only the exponent βββ aligns closely with the predicted literature values, on the other hand, the other exponents show significant deviations.
Thermal and magnetic effects in a system consisting of thin layers of coupled Ising spins with S=1/2 and S=1 are considered. The specific heat and the correlation length display maxima at two different temperatures. It is discussed in what sense these maxima can be interpreted as a finite-size rounding of a thermodynamic singularity associated with a phase transition. The connection with ordinary, extraordinary, and special surface phase transitions is made. In two dimensions (2D), the surface critical exponents are calculated from conformal invariance. The bulk and surface finite-size scaling of the order-parameter profiles at the transition points is discussed. In 2D, an exact scaling function for the profiles is suggested through conformal invariance arguments for the (extra)ordinary transition.
In this review, we show our results with new interpretation on the critical exponents of thin films obtained by high-performance multi-histogram Monte Carlo simulations. The film thickness NzN_zNz​ consists of a few layers up to a dozen of layers in the zzz direction. The free boundary condition is applied in this direction while in the xyxyxy plane periodic boundary conditions are used. Large xyxyxy plane sizes are used for finite-size scaling. The Ising model is studied with nearest-neighbor (NN) interaction. When Nz=1N_z=1Nz​=1, namely the two-dimensional (2D) system, we find the critical exponents given by the renormalization group. While, for Nz>1N_z>1Nz​>1, the critical exponents calculated with the high-precision multi-histogram technique show that they deviate slightly but systematically from the 2D values. If we use these values of critical exponents in the hyperscaling relation with d=2d=2d=2, then the hyperscaling relation is violated. However, if we use the hyperscaling relation and the critical exponents obtained for Nz>1N_z>1Nz​>1 to calculate the dimension of the system, we find the system dimension slightly larger than 2. This can be viewed as an "effective" dimension. More discussion is given in the paper. We also show the cross-over between the first- and second-order transition while varying the film thickness in an antiferromagnetic FCC Ising frustrated thin film. In addition, we will show evidence that when a 2D system has two order parameters of different symmetries with a single transition, the critical exponents are new, suggesting a universality class of coupled two-symmetry breakings. In this case, the 2D hyperscaling does not hold. Another case is the 3D Ising model coupled to the lattice vibration: the critical exponents deviate from the 3D Ising ones, the results suggest the violation of the hyperscaling.
The corrections to finite-size scaling in the critical two-point correlation function G(r) of 2D Ising model on a square lattice have been studied numerically by means of exact transfer-matrix algorithms. The systems have been considered, including up to 800 spins. The calculation of G(r) at a distance r equal to the half of the system size L shows the existence of an amplitude correction proportional to 1/L^2. A nontrivial correction ~1/L^0.25 of a very small magnitude also has been detected, as it can be expected from our recently developed GFD (grouping of Feynman diagrams) theory. Monte Carlo simulations of the squared magnetization of 3D Ising model have been performed by Wolff's algorithm in the range of the reduced temperatures t =< 0.000086 and system sizes L =< 410. The effective critical exponent beta_eff(t) tends to increase above the currently accepted numerical values. The critical coupling K_c=0.22165386(51) has been extracted from the Binder cumulant data within 96 =< L =< 384. The critical exponent 1/nu, estimated from the finite-size scaling of the derivatives of the Binder cumulant, tends to decrease slightly below the RG value 1.587 for the largest system sizes. The finite-size scaling of accurately simulated maximal values of the specific heat C_v in 3D Ising model confirms a logarithmic rather than power-like critical singularity of C_v.
We have used Monte Carlo simulations and mean-field analysis to observe the magnetic behavior of Ising thin films with cubic lattice structures as a function of temperature and thickness, especially in the critical region. Magnetization and magnetic susceptibility, including layer variation, are investigated. We find that the magnetic behavior changes from two-dimensional to three-dimensional character with increasing film thickness. Both the crossover of the critical temperature from a two-dimensional to a bulk value and the shift exponent are observed. Nevertheless, with support from a scaling function, the simulations show that the effective critical exponents for films with large enough layer extents only vary a little from their two-dimensional values. This, in particular, provides an indication of two-dimensional universality in the thin films.
The surface scaling theory previously presented by the authors is developed further, and derived heuristically from a cluster model. Monte Carlo calculations are carried out to obtain the spatial and temperature dependence of the magnetization in Ising and Heisenberg systems with free surfaces. The exponent ${\ensuremath{\beta}}_{1}$ of the (surface) layer magnetization is shown to agree with the scaling value (${\ensuremath{\beta}}_{1}\ensuremath{\approx}\frac{2}{3}$) previously derived. In the Heisenberg system, the results at low temperature agree with a spin-wave calculation by Mills and Maradudin. Ising models with modified exchange ${J}_{s}=J(1+\ensuremath{\Delta})\ensuremath{\ne}J$ on the surface are considered, both in mean-field theory and by means of high-temperature-series expansions. The critical value ${\ensuremath{\Delta}}_{c}$ for surface ordering is found from the series to be 0.6, compared to the mean-field value of 0.25. For $\ensuremath{\Delta}>{\ensuremath{\Delta}}_{c}$ there is a temperature region in which the surface behaves like a bulk two-dimensional Ising model near its phase transition. The critical exponents experience a crossover at $\ensuremath{\Delta}={\ensuremath{\Delta}}_{c}$, which is reflected in poorly behaved series, and effective exponents differing from the true ones for $\ensuremath{\Delta}\ensuremath{\lesssim}{\ensuremath{\Delta}}_{c}$. In the case of weakened surface exchange (0<Js<J0<{J}_{s}<J0<Js​<J), the layer magnetization is shown to fit a linear temperature dependence over a large temperature range below Tc{T}_{c}Tc​, thus providing a possible explanation for previous experiments. For sufficiently strong negative Js{J}_{s}Js​, mean-field theory predicts that the surface will order antiferromagnetically while the bulk is ferromagnetic.
Thermodynamics of a layered Ising model with infinite-range ferromagnetic intralayer interaction and random nearest-neighbor interlayer coupling is considered. A detailed analysis of the model with vanishing average interlayer coupling is presented. The Gibbs free energy is found in the critical region, and the existence of many metastable states is demonstrated. Thermodynamic parameters of the system are found for periodic states. As the mean square interlayer coupling increases, the equilibrium state of the system undergoes an infinite sequence of first-order phase transitions, the number of magnetic planes and the distance between them change discontinuously, and so do both bulk magnetization and magnetic susceptibility.