Next Article in Journal
Ideal Agent System with Triplet States: Model Parameter Identification of Agent–Field Interaction
Next Article in Special Issue
QUBO Problem Formulation of Fragment-Based Protein–Ligand Flexible Docking
Previous Article in Journal
Dynamic Feature Extraction-Based Quadratic Discriminant Analysis for Industrial Process Fault Classification and Diagnosis
Previous Article in Special Issue
The Capabilities of Boltzmann Machines to Detect and Reconstruct Ising System’s Configurations from a Given Temperature
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ising Ladder with Four-Spin Plaquette Interaction in a Transverse Magnetic Field

by
Maria Eugenia S. Nunes
1,
Francisco Welington S. Lima
2 and
Joao A. Plascak
3,4,5,*
1
Departamento de Física, Universidade Federal de Ouro Preto, Campus Universitário Morro do Cruzeiro, ICEB, Ouro Preto 35400-000, MG, Brazil
2
Dietrich Stauffer Computational Physics Lab, Departamento de Física, Universidade Federal do Piauí, Teresina 64049-550, PI, Brazil
3
Departamento de Física, Centro de Ciências Exatas e da Natureza CCEN, Universidade Federal da Paraíba, Cidade Universitária, João Pessoa 58051-970, PB, Brazil
4
Departamento de Física, Universidade Federal de Minas Gerais, C. P. 702, Belo Horizonte 30123-970, MG, Brazil
5
Department of Physics and Astronomy, University of Georgia, Athens, GA 30602, USA
*
Author to whom correspondence should be addressed.
Entropy 2023, 25(12), 1665; https://doi.org/10.3390/e25121665
Submission received: 6 November 2023 / Revised: 30 November 2023 / Accepted: 12 December 2023 / Published: 16 December 2023
(This article belongs to the Special Issue Ising Model: Recent Developments and Exotic Applications II)

Abstract

:
The spin-1/2 quantum transverse Ising model, defined on a ladder structure, with nearest-neighbor and four-spin interaction on a plaquette, was studied by using exact diagonalization on finite ladders together with finite-size-scaling procedures. The quantum phase transition between the ferromagnetic and paramagnetic phases has then been obtained by extrapolating the data to the thermodynamic limit. The critical transverse field decreases as the antiferromagnetic four-spin interaction increases and reaches a multicritical point. However, the exact diagonalization approach was not able to capture the essence of the dimer phase beyond the multicritical transition.

1. Introduction

Low-dimensional quantum spin models have long been the subject of many theoretical and experimental approaches. The interest in such systems mainly derives from two factors. The first is the high experimental appeal because some physical realizations, when viewed by the relevant microscopic quantum interactions, consist, in fact, of single molecules, linear chains, or even ladder-like structures. To cite a few examples among the vast list of experimental results that fall into these scenarios, we have (i) molecular magnetism in Cu5-NIPA [1], V6-like magnetic molecules [2], and multiferroics [3] (these compounds can be considered zero-dimensional systems); (ii) TMMC, CsNiF 3 , and CuCl 2 2NC 5 H 5 systems exhibiting one-dimensional characteristics [4]; (iii) ladder-type structures as in (VO) 2 P 2 O 7 [5], Cu(C 5 H 12 N 2 ) 2 2Cl 4 [6], as well as La 2 CuO 4 and La 6 Ca 8 Cu 24 O 41 cuprates [7,8,9,10]. Of course, the majority of physical realizations indeed have the main quantum interactions along the three crystalline directions. Second, it is a real challenge to theoretically treat quantum models even in low spatial dimensions, since the statistical mechanics of a d-dimensional quantum model (d being the dimension of the lattice) is equivalent to a d + 1 -classical model [11], implying that even one-dimensional quantum systems can undergo a critical phase transition at zero temperature, the so-called quantum phase transition.
Regarding the low-dimensional theoretical models, the Ising, XY, and Heisenberg quantum spin systems have been widely studied in the literature, mostly with nearest-neighbor interactions and also in the presence of transverse fields (see, for example, Refs. [11,12]). In particular, spin-1/2 XY chains, with the transverse Ising case being a projection of the XY Hamiltonian when no interaction in the y-direction is present, have been studied not only for their static phase transition behavior but also from a dynamical perspective by computing the time-dependent correlation functions and the corresponding spectral functions [13,14,15,16]. On the other hand, spin-1/2 anisotropic Heisenberg chains have been treated according to different approaches, including the Bethe ansatz and Jordan–Wigner transformations (for example, see Ref. [17] and references therein). Furthermore, the effects of randomness in the isotropic XY chain [17], which makes the theoretical treatment even more difficult to tackle, have been introduced in the nearest-neighbor exchange and other types of interactions, as well as in the transverse field.
Despite the nearest-neighbor exchange being the relevant interaction in all physical realizations and models cited above, four-spin interactions in a plaquette have been suggested to reproduce the dispersion relation observed in inelastic neutron scattering experiments on cuprates [7,8,9,10] and even to stabilize a chiral spin liquid on a triangular lattice [18]. This new experimental motivation has lead to an increase in theoretical investigations of quantum Ising–Heisenberg-like models defined on a ladder [19,20], especially with four-spin interaction, which has been shown to induce other unusual types of order, such as scalar chirality and intra-rail staggered dimerization (see, for instance, Ref. [21] and references therein).
The spin-1/2 quantum transverse Ising model, which is the simplest model defined on a ladder structure with nearest-neighbor interaction, has thus been widely studied in the literature. In particular, the model including a four-spin interaction on a plaquette has been recently considered by using the density matrix renormalization group (DMRG) approach [21]. (The dynamical behavior has also been recently studied via the recurrence relation method [22]) The ground state phase diagram has been obtained in the transverse field versus four-spin interaction plane. In this case, a ferromagnetic-to-paramagnetic second-order phase transition is observed, and, for sufficiently high antiferromagnetic four-spin couplings, the system presents a dimer phase and a multicritical point. Although the transition for ferromagnetic four-spin couplings is well established, with the two degenerated ferromagnetic phases becoming equal to the paramagnetic phase at the transition point, the situation is less clear for antiferromagnetic four-spin couplings, where the ordered phase is given by a superposition of dimerized states [21]. In view of this fact, it is interesting to investigate this model with different techniques in order to obtain a better picture of the phases along the dimerized transition line. Since the finite-size-scaling approach, based on exact results on finite lattices, has been shown to be one of the most accurate methods to treat such models, it was employed here to obtain the quantum phase transition of the transverse Ising ladder with four-spin interactions. It would also be quite interesting and illustrative to compare the present results of the phase diagram to those obtained for the one-dimensional model with linear four-spin coupling, which have been obtained by using the same exact diagonalization procedure [23].
Thus, based on the above, in the present work, we have revisited the transverse Ising model with four-spin interaction in a ladder structure and used exact diagonalization on finite lattices. The energy gap has been computed, and, from the corresponding finite-size-scaling (FSS) relation, the critical transverse field has been obtained for several values of four-spin interaction. The ferromagnetic–paramagnetic phase transition line has been obtained by extrapolating the data to the thermodynamic limit. The results agree well with those previously obtained from the DMRG [21]. However, contrary to the model with four-spin interactions in one dimension [23], it has been noticed that, in the ladder structure, there are much stronger finite-size effects, and one indeed has to use the larger possible lattices for the extrapolations. In addition, the ferromagnetic–paramagnetic transition line on the ladder does not go to zero at the multiphase point in the classical limit of vanishing transverse field, as does the one-dimensional model.
The remainder of the paper is structured as follows. In the next section, the model is defined, and the configurations in the classical limit are discussed. Section 3 describes the theoretical approach using the exact diagonalization on finite ladders and the corresponding FSS relation for the energy gap and critical transverse field. The results of the transition line are presented in Section 4, and some concluding remarks are addressed in the final section.

2. Model

The Hamiltonian of the spin-1/2 Ising model, defined on a ladder structure, as depicted in Figure 1, with four-spin interaction and in the presence of a transverse field, can be written as
H = J i = 1 L σ 1 , i z σ 1 , i + 1 z + σ 2 , i z σ 2 , i + 1 z J i = 1 L σ 1 , i z σ 2 , i z + J 4 i = 1 L σ 1 , i z σ 1 , i + 1 z σ 2 , i z σ 2 , i + 1 z B i = 1 L ( σ 1 , i x + σ 2 , i x ) ,
where J > 0 is the nearest-neighbor exchange interaction along the rails and in the rungs, J 4 is the four-spin interaction connecting the spins in a plaquette, and B is the transverse external field applied in the x direction. The sums are over the L rungs of a ladder with periodic boundary conditions in the direction of the two side rails denoted, respectively, by 1 and 2. The spin-1/2 operators σ , i α , with = 1 , 2 and α = x , z , are given by the Pauli spin matrices.
Due to the positiveness convention of the four-spin interaction in the Hamiltonian (1), when J 4 < 0 , the plaquette ordering is also ferromagnetic for low values of transverse fields. As a result, the system undergoes a second-order transition from a ferromagnetic phase to a paramagnetic phase at a critical transverse field value B c . This critical transverse field decreases as J 4 increases. As has been shown in Ref. [21], when J 4 > 0 , this transition persists until some value of J 4 (that is B-dependent) where a dimer phase is set up in the ladder.
At B = 0 , one has a classical system with a multiphase point at J 4 = 3 / 2 dividing the classical axis J 4 into two regions: (1) for J 4 < 3 / 2 a double degenerated ferromagnetic phase and (2) for J 4 3 / 2 a ground state that is 2 L + 1 degenerate, with a residual entropy S r = ( L + 1 ) k B ln 2 (here, k B is the Boltzmann constant).
Contrary to the one-dimensional model with four-spin interaction along a straight line [23], the ferromagnetic–paramagnetic transition line does not terminate at the multiphase point at B = 0 . This is a consequence of the ladder structure allowing the appearance of the rung dimerized phase [21].

3. Theoretical Background

We have used finite ladders of length L, with N = 2 L sites, and periodic boundary conditions. For each ladder, the energy gap G L ( B , J 4 ) has been computed. The quantity G L ( B , J 4 ) is given by
G L ( B , J 4 ) = E L 1 E L 0 ,
where E L 1 E L 0 is the energy gap between the first excited state and the ground state, respectively.
G L ( B , J 4 ) is equivalent to the correlation length in thermal systems and satisfies the FSS relation [24]
L G L ( B c L , J 4 ) = L G L ( B c L , J 4 ) ,
for two finite ladders of length L and L , with usually L > L . From the above relation, it is possible to estimate B c L , the critical field for the ladder pair ( L , L ).
For sufficiently large lattices, it is expected that the quantities L G L ( B , J 4 ) , as a function of B for a given value of J 4 , cross at the same point B c , the critical transverse field for the considered value of J 4 . However, in some cases, where the size of the finite lattices is insufficient, residual corrections make the crossing points B c L suffer a systematic shift as L varies. Nevertheless, there is an additional FSS relation for the L dependence on the transition points B c L given by [25]
B c L = B c + a L 1 / ν 1 + b L ω ,
where B c = B c L is the critical transverse field in the thermodynamic limit, ν is the correlation length critical exponent, ω is the correction-to-scaling exponent, and a and b are non-universal constants. In this way, for every chosen value of J 4 and reference ladder L , one obtains the crossing point B c L for the ladder pair ( L , L ) through Equation (3). Next, using Equation (4) for various values of L, the desired extrapolated value of critical transverse field can be computed.
In the present work, the energy gap has been obtained through exact diagonalization of finite ladders with length L in the range 2 L 11 . For the larger lattices, we have employed the Lanczos diagonalization procedure [26].

4. Results and Discussion

In what follows, J = 1 . This value of J can also be interpreted as measuring the four-spin interaction J 4 , as well as the transverse field strength B, in units of J.
As an example of the behavior of energy gap as a function of the transverse field, Figure 2 shows L G L ( B , J 4 ) for several values of L and J 4 = 0 . In this case, we simply have the transverse Ising model on a ladder [27]. It is clear that a small region of unique crossings is only achieved for the larger ladders. The estimate of the critical transverse field, in this case, is revealed by the corresponding arrow.
It is also evident in the figure that for smaller ladders, there is a shift in the values of crossings B c L . We can thus compute B c L by taking a reference lattice L and several values of L > L . From the results of Figure 2, it is possible to obtain a clear crossing by using L = 2 , 3 , and 4. The results obtained are depicted in Figure 3, which allows one to obtain additional extrapolated estimates of B c through fits to Equation (4) for each value of the reference ladder L . The present estimate B c = 1.8322(2) agrees very well with B c = 1.83213 from the DMRG of Ref. [21] using a finite ladder with length L = 28 .
It should be stressed that in using Equation (4), the exponents ν and θ are required. Although from Equation (3) it is possible to compute ν using renormalization group ideas and estimates of B c [23,28] (for instance, from the values of the crossings B c L for the larger ladders), the obtained results are close to the exact ones for this universality class. However, the correction-to-scaling exponent cannot be obtained in this way and, at least, should be treated as an extra adjustable parameter. For this reason, and due to the fact that we are indeed more interested in the location of the transition line, we have resorted to the known values ν = 1 and θ = 2 for the Ising universality class in order to compute the extrapolated critical transverse field.
Similar curves, with similar FSS analysis, are obtained when J 4 > 0 . Figure 4 shows the critical transverse field B c as a function of J 4 obtained from the present method in comparison to some results from the DMRG of Ref. [21]. One can clearly see that, contrary to the very same model in one dimension [23] (with four-spin interactions along a straight line), the transition line definitely does not go to zero as J 4 tends to 3/2, the multiphase point. The agreement with DMRG results is not only apparent in Figure 4 but also in Table 1, which gives a more detailed numerical comparison for some selected values of J 4 .

5. Concluding Remarks

The transverse Ising model, defined on a ladder and with four-spin interactions on a plaquette, has been studied by using the finite-size-scaling approach with exact diagonalization of finite ladders. The method has achieved acceptable results for the ferromagnetic-to-paramagnetic transition line in the transverse field versus four-spin interaction plane. The computed critical transverse field in the thermodynamic limit is also in good agreement with that from the density matrix renormalization group procedure.
Despite the expected efficacy of the FSS along the ferromagnetic–paramagnetic transition line, the model shows, in fact, strong finite-size effects, and only with considerably larger lattices could the results be accurately obtained. The finite-size effects here are, in fact, stronger than those seen in the one-dimensional model [28]. The ferromagnetic–paramagnetic transition line definitely does not meet the multiphase point at zero field, and there is no reentrancy along this transition line whatsoever. This behavior is contrary to the one-dimensional model although in agreement with the DMRG approach providing a dimerized phase for sufficiently high values of antiferromagnetic four-spin interaction.
One unexpected feature is the exact diagonalization on finite lattices method being unable to locate the dimerized phase and the corresponding transition to the ferromagnetic and paramagnetic phases that occurs. For values J 4 > 1.9 , the crossings become more diffuse, and in some regions it is not possible to obtain good fits with the expected scaling relation. More suitable quantities than the energy gap may be necessary to unveil the transition characteristics and the microscopic spin behavior for larger values of four-spin interaction by using the FSS procedure. Accordingly, a detailed characterization of the nature of the multicritical point, whether tricritical, tetracritical, etc. [29], remains to be conducted.

Author Contributions

Conceptualization, J.A.P. and M.E.S.N.; Methodology, F.W.S.L.; Formal analysis, M.E.S.N.; Investigation, M.E.S.N.; Writing—original draft, J.A.P.; Writing—review & editing, J.A.P., M.E.S.N. and F.W.S.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work has been partially funded by the Brazilian agencies CNPq, CAPES, FAPEMIG, and FUNCAP.

Data Availability Statement

Data are contained within the article.

Acknowledgments

The authors would like to thank F. C. Xavier and R. G. Pereira for fruitful discussions and the former for the invaluable help in implementing the Lanczos algorithm in the ladder lattice. Special thanks are due to S.-H. Tsai for the help in using the GACRC facilities of the cluster at the University of Georgia in the USA and G. Weber for the use of the workstations of the Statistical Physics group at UFMG in Brazil. Computational support by the system SGI Altix 1350, the CENAPAD.UNICAMP-USP computational park, Campinas, São Paulo - BRAZIL, is also gratefully acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Nath, R.; Tsirlin, A.A.; Khuntia, P.; Janson, O.; Förster, T.; Padmanabhan, M.; Li, J.; Skourski, Y.; Baenitz, M.; Rosner, H.; et al. Magnetization and spin dynamics of the spin S = 1/2 hourglass nanomagnet Cu5(OH)2 (NIPA) 4·10H2O. Phys. Rev. B 2013, 87, 214417. [Google Scholar] [CrossRef]
  2. Torrico, J.; Plascak, J.A. Ground state and thermodynamic properties of spin-1/2 isosceles Heisenberg triangles for V6-like magnetic molecules. Phys. Rev. E 2020, 102, 062116. [Google Scholar] [CrossRef] [PubMed]
  3. Valentim, A.; Garcia, D.J.; Plascak, J.A. Magnetoelectric coupling on fused azulene oligomers. Phys. Rev. B 2022, 105, 174426. [Google Scholar] [CrossRef]
  4. Mikeska, H.J. Magnetic chains in solids. J. Magn. Magn. Mater. 1979, 13, 35–49. [Google Scholar] [CrossRef]
  5. Johnston, D.C.; Johnson, J.W.; Goshorn, D.P.; Jacobson, A.J. Magnetic susceptibility of (VO)2 P2O7: A one-dimensional spin-1/2 Heisenberg antiferromagnet with a ladder spin configuration and a singlet ground state. Phys. Rev. B 1987, 35, 219. [Google Scholar] [CrossRef] [PubMed]
  6. Hagiwara, M.; Katori, H.A.; Schollwöck, U.; Mikeska, H.J. Specific heat of an S = 1/2 Heisenberg ladder compound Cu2(C5H12N2)2Cl4 under magnetic fields. Phys. Rev. B 2000, 62, 1051. [Google Scholar] [CrossRef]
  7. Brehmer, S.; Mikeska, H.J.; Müller, M.; Nagaosa, N.; Uchida, S. Effects of biquadratic exchange on the spectrum of elementary excitations in spin ladders. Phys. Rev. B 1999, 60, 329. [Google Scholar] [CrossRef]
  8. Matsuda, M.; Katsumata, K.; Eccleston, R.S.; Brehmer, S.; Mikeska, H.J. Magnetic excitations and exchange interactions in the spin-1/2 two-leg ladder compound La6Ca8Cu24O41. Phys. Rev. B 2000, 62, 8903. [Google Scholar] [CrossRef]
  9. Coldea, R.; Hayden, S.M.; Aeppli, G.; Perring, T.G.; Frost, C.D.; Mason, T.E.; Cheong, S.W.; Fisk, Z. Spin waves and electronic interactions in La2CuO4. Phys. Rev. Lett. 2011, 86, 5377. [Google Scholar] [CrossRef]
  10. Larsen, C.B.; Rømer, A.T.; Janas, S.; Treue, F.; Mønsted, B.; Shaik, N.E.; Rønnow, H.M.; Lefmann, K. Exact diagonalization study of the Hubbard-parametrized four-spin ring exchange model on a square lattice. Phys. Rev. B 2019, 99, 054432. [Google Scholar] [CrossRef]
  11. Sachdev, S. Quantum Phase Transitions, 2nd ed.; Cambridge University Press: New York, NY, USA, 2011. [Google Scholar]
  12. Parkinson, J.B.; Farnell, D.J.J. An Introduction to Quantum Spin Systems; Lecture Notes Physic 816; Springer: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
  13. Derzhko, O.; Richter, J. Solvable model of a random spin-XY chain. Phys. Rev. B 1997, 55, 14298. [Google Scholar] [CrossRef]
  14. Derzhko, O.; Krokhmalskii, T. Dynamic structure factor of the spin-1/2 transverse Ising chain. Phys. Rev. B 1997, 56, 11659. [Google Scholar] [CrossRef]
  15. Menchyshyn, O.; Ohanyan, V.; Verkholyak, T.; Krokhmalskii, T.; Derzhko, O. Magnetism-driven ferroelectricity in spin-1/2 X Y chains. Phys. Rev. B 2015, 92, 184427. [Google Scholar] [CrossRef]
  16. Derzhko, O.; Krokhmalskii, T.; Stolze, J. Dynamics of the spin-1/2 isotropic XY chain in a transverse field. J. Phys. A Math. Gen. 2000, 33, 3063. [Google Scholar] [CrossRef]
  17. Verkholyak, T.; Derzhko, O.; Krokhmalskii, T.; Stolze, J. Dynamic properties of quantum spin chains: Simple route to complex behavior. Phys. Rev. B 2007, 76, 144418. [Google Scholar] [CrossRef]
  18. Cookmeyer, T.; Motruk, J.; Moore, J.E. Four-spin terms and the origin of the chiral spin liquid in Mott insulators on the triangular lattice. Phys. Rev. Lett. 2021, 127, 087201. [Google Scholar] [CrossRef] [PubMed]
  19. Derzhko, O.; Krokhmalskii, T.; Richter, J. Emergent Ising degrees of freedom in frustrated two-leg ladder and bilayer s = 1/2 Heisenberg antiferromagnets. Phys. Rev. B 2010, 82, 214412. [Google Scholar] [CrossRef]
  20. Zad, H.A.; Ohanyan, V.; Zoshki, A.; Strečka, J. Magnetization plateaus and enhanced magnetocaloric effect of a spin-1/2 Ising-Heisenberg and Heisenberg double sawtooth ladder with four-spin interaction. Phys. Rev. E 2023, 108, 044132. [Google Scholar]
  21. Xavier, J.C.; Pereira, R.G.; Nunes, M.E.S.; Plascak, J.A. Coexistence of spontaneous dimerization and magnetic order in a transverse-field Ising ladder with four-spin interactions. Phys. Rev. B 2022, 105, 024430. [Google Scholar] [CrossRef]
  22. de Souza, W.L.; de Mello Silva, É.; Martins, P.H. Dynamics of the spin-1/2 Ising two-leg ladder with four-spin plaquette interaction and transverse field. Phys. Rev. E 2020, 101, 042104. [Google Scholar] [CrossRef]
  23. Guimaraes, P.R.C.; Plascak, J.A.; Barreto, F.C.S.; Florencio, J. Quantum phase transitions in the one-dimensional transverse Ising model with second-neighbor interactions. Phys. Rev. B 2002, 66, 064413. [Google Scholar] [CrossRef]
  24. Barber, M.N. Phase Transitions and Critical Phenomena; Domb, C., Lebowitz, J.L., Eds.; Academic Press: London, UK, 1983. [Google Scholar]
  25. Landau, D.P.; Binder, K. A Guide to Monte Carlo Simulation in Statistical Physics, 5th ed.; Cambridge University Press: Cambridge, UK, 2021. [Google Scholar]
  26. Golub, G.H.; van Loan, C.F. Matrix Computations, 3rd ed.; Johns Hopkins University Press: Baltimore, MD, USA, 1996. [Google Scholar]
  27. Xavier, J.C.; Ramos, F.B. Universal finite-size corrections of the entanglement entropy of quantum ladders and the entropic area law. J. Stat. Mech. Theory Exp. 2014, 2014, P10034. [Google Scholar] [CrossRef]
  28. Bonfim, O.F.d.; Florêncio, J. Quantum phase transitions in the transverse one-dimensional Ising model with four-spin interactions. Phys. Rev. B 2006, 74, 134413. [Google Scholar] [CrossRef]
  29. Dias, D.A.; Lima, F.W.S.; Plascak, J.A. Generalized Gibbs Phase Rule and Multicriticality Applied to Magnetic Systems. Entropy 2021, 24, 63. [Google Scholar] [CrossRef]
Figure 1. Sketch of a ladder with rails labeled as 1 and 2. The interaction along the rails and in the rungs is J. The circles represent the four-spin interaction J 4 in each plaquette. σ n , i z , with n = 1 , 2 , are the z-component spin operators. The z axis is taken to be perpendicular to the plane of the figure and the transverse field B is along the x axis.
Figure 1. Sketch of a ladder with rails labeled as 1 and 2. The interaction along the rails and in the rungs is J. The circles represent the four-spin interaction J 4 in each plaquette. σ n , i z , with n = 1 , 2 , are the z-component spin operators. The z axis is taken to be perpendicular to the plane of the figure and the transverse field B is along the x axis.
Entropy 25 01665 g001
Figure 2. The quantity L G L ( B , J 4 ) , defined in the text, as a function of the transverse field B, for several values of L (indicated in the curves) and J 4 = 0 . The arrow indicates the estimate of B c when considering only the larger ladders.
Figure 2. The quantity L G L ( B , J 4 ) , defined in the text, as a function of the transverse field B, for several values of L (indicated in the curves) and J 4 = 0 . The arrow indicates the estimate of B c when considering only the larger ladders.
Entropy 25 01665 g002
Figure 3. Crossing transverse field B c L , as a function of L 1 / ν ( ν = 1 ), for J 4 = 0 and different values of the reference ladder L (indicated in the legend). In the legend we find the extrapolated values of B c , using fits to Equation (4) with θ = 2 , for each curve. Furthermore, indicated in the legend come the values from the crossings of the larger ladders in Figure 2 (cross) and from Ref. [21].
Figure 3. Crossing transverse field B c L , as a function of L 1 / ν ( ν = 1 ), for J 4 = 0 and different values of the reference ladder L (indicated in the legend). In the legend we find the extrapolated values of B c , using fits to Equation (4) with θ = 2 , for each curve. Furthermore, indicated in the legend come the values from the crossings of the larger ladders in Figure 2 (cross) and from Ref. [21].
Entropy 25 01665 g003
Figure 4. Critical values of the transverse field B c as a function of the four-spin interaction J 4 . FM and PM stand for the ferromagnetic and paramagnetic phases, respectively. The full line represents the present results and the circles the values from the DMRG of Ref. [21]. The diamond is located at the multiphase point at B c = 0 .
Figure 4. Critical values of the transverse field B c as a function of the four-spin interaction J 4 . FM and PM stand for the ferromagnetic and paramagnetic phases, respectively. The full line represents the present results and the circles the values from the DMRG of Ref. [21]. The diamond is located at the multiphase point at B c = 0 .
Entropy 25 01665 g004
Table 1. Critical values of the transverse field B c in the thermodynamic limit for some values of the four-spin interaction J 4 . The first row gives the results from the present work and the second row those from the DMRG calculation [21] with L = 28 .
Table 1. Critical values of the transverse field B c in the thermodynamic limit for some values of the four-spin interaction J 4 . The first row gives the results from the present work and the second row those from the DMRG calculation [21] with L = 28 .
J 4 00.511.51.71.85
B c 1.8322(2)1.5726(2)1.2938(2)0.9813(3)0.834(3)0.697(2)this work
1.83213  1.57226  1.29334  0.98041  0.83579 0.70510 DMR [21]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Nunes, M.E.S.; Lima, F.W.S.; Plascak, J.A. Ising Ladder with Four-Spin Plaquette Interaction in a Transverse Magnetic Field. Entropy 2023, 25, 1665. https://doi.org/10.3390/e25121665

AMA Style

Nunes MES, Lima FWS, Plascak JA. Ising Ladder with Four-Spin Plaquette Interaction in a Transverse Magnetic Field. Entropy. 2023; 25(12):1665. https://doi.org/10.3390/e25121665

Chicago/Turabian Style

Nunes, Maria Eugenia S., Francisco Welington S. Lima, and Joao A. Plascak. 2023. "Ising Ladder with Four-Spin Plaquette Interaction in a Transverse Magnetic Field" Entropy 25, no. 12: 1665. https://doi.org/10.3390/e25121665

APA Style

Nunes, M. E. S., Lima, F. W. S., & Plascak, J. A. (2023). Ising Ladder with Four-Spin Plaquette Interaction in a Transverse Magnetic Field. Entropy, 25(12), 1665. https://doi.org/10.3390/e25121665

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop