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

Physical review. B, Condensed matter

Numerical study of the ground-state properties of a frustratedXYmodel

P. Gawiec, Daniel R. Grempel

Published Aug 1, 1991• 21 citations
Publisher

Abstract

We study numerically the ground-state morphology and properties of the low-lying states of an XY model with competing interactions for several values of the concentration x and of the ratio between ferromagnetic and antiferromagnetic exchange constants, \ensuremath{\lambda}. If \ensuremath{\lambda}>1, the ground state is likely to be a spin glass for all nonzero concentrations. In the region \ensuremath{\lambda}\ensuremath{\le}1, the ground state is a disordered canted ferromagnet for low dilutions with a transition to a pure spin-glass phase at a \ensuremath{\lambda}-dependent critical concentration. The analysis of the spin-spin correlation functions gives insight on the mechanism responsible for the degradation of the magnetic order at low dilutions. An unusual law of decay of correlations suggests that the canted ferromagnetic phase that we find at small \ensuremath{\lambda} and x may actually be a phase with power-law decay of correlations.

DOI
10.1103/physrevb.44.2613

Related papers

Ground State Properties of the ±JIsing Model in Two Dimensions

Yukiyasu Ozeki

We investigate the ground state properties of the ± J Ising model with asymmetric probability weight in two dimensions. We propose a numerical transfer matrix method which enables us to calculate directly the ground-state properties such as correlation functions, defect energy and entropy at T =0. We find that the spin glass correlation function decays in a power law like \(\tilde{g}(r)\sim{r}^{-\tilde{\eta}}\) with \(\tilde{\eta}\simeq0.21\) in the range \(0.5\leq{p}\lesssim0.85\), suggesting the weak universality of spin glass-the concentration of + J bond is denoted by p . At the ferromagnetic phase boundary locating around p =0.89, the ferromagnetic correlation function decays as \(g(r)\sim{r}^{-\tilde{\eta}}\) with η≃0.240±0.016, which suggests the weak universality of ferromagnetic critical exponent from the non-random critical point ( p =1, T =2.269···) to the ground sate counterpart. The behavior of correlation functions, and also that of defect energies indicate the existence of an extra-phase, called the random antiphase state, in \(0.85\lesssim{p}\lesssim0.89\).

Journal of the Physical Society of Japan•Oct 15, 1990•61 citations

Ferrimagnetism in a disordered Ising model

Giovanni Paladin, Michele Pasquini, Maurizio Serva

We introduce a one dimensional Ising model with two competing interactions: nearest neighbor random couplings +/-J with equal probability and a positive infinite range coupling Lambda. At low temperature T the model exhibits a first order phase transition between a ferromagnetic state (with magnetization m(1) = 1 at T = 0) and a << ferrimagnetic >> state (with m(2) = 2/3 at T = 0), when the disorder strength J/Lambda is increased. For 5/12 < J/Lambda < 1, a whole spectrum of ferrimagnetic ground states with magnetization m(n) = 2/(n + 1) (n = 2, ..., infinity) is present while for J/lambda > 1 the ground state is given by a trivial one dimensional spin glass with m = 0. The main qualitative features of the model can be described by a simplified annealed model where the random couplings can arrange themselves to minimize free energy with the constraint that the number of positive couplings is fixed by the law of large numbers in the thermodynamic limit. This model is exactly solved at all temperatures and the diagram of phase is calculated.

Journal de Physique I•Nov 1, 1994•7 citations

XY model with antinematic interaction

Milan ˇ Zukovič

We consider the XY model with ferromagnetic (FM) and antinematic (AN) nearest-neighbor interactions on a square lattice for a varying interaction strength ratio. Besides the expected FM and AN quasi-long-range order (QLRO) phases we identify at low temperatures another peculiar canted ferromagnetic (CFM) QLRO phase, resulting from the competition between the collinear FM and noncollinear AN ordering tendencies. In the CFM phase neighboring spins that belong to different sublattices are canted by a nonuniversal (dependent on the interaction strength ratio) angle and the ordering is characterized by a fast-decaying power-law intrasublattice correlation function. Compared to the FM phase, in the CFM phase correlations are significantly diminished by the presence of zero-energy domain walls due to the inherent degeneracy caused by the AN interactions. We present the phase diagram as a function of the interaction strength ratio and discuss the character of the respective phases as well as the transitions between them.

arXiv:1903.05586•Physical review. E•Jun 12, 2019•11 citations

Spin-spin correlation function of the two-dimensional XY model with weak site or bond dilution

O. Kapikranian

The spin-spin correlation function of the two-dimensional (2D) $XY$ model decays as a power law at all temperatures below the Berezinskii-Kosterlitz-Thouless transition point with a temperature-dependent exponent $\ensuremath{\eta}=\ensuremath{\eta}(T/J)$ ($J$ is the ferromagnetic coupling strength). It is known from computer experiments that in the 2D $XY$ model with site or bond dilution this exponent depends on the concentration $p$ of removed sites or bonds as well. Knowing the slope $\ensuremath{\partial}\ensuremath{\eta}/\ensuremath{\partial}p$ at point $p=0$, one can predict the value of the exponent for small dilution concentrations: ${\ensuremath{\eta}(p)\ensuremath{\simeq}\ensuremath{\eta}(0)+p(\ensuremath{\partial}\ensuremath{\eta}/\ensuremath{\partial}p)|}_{p=0}$. As shown in this paper, the spin-wave Hamiltonian allows one to obtain exact results for this slope: ${(\ensuremath{\partial}\ensuremath{\eta}/\ensuremath{\partial}p)|}_{p=0}=T/(2J)+O({(T/J)}^{2})$ and $T/(\ensuremath{\pi}J)+O({(T/J)}^{2})$ for site and for bond dilution, respectively.

arXiv:1111.3939•Physical Review B•Mar 5, 2012•3 citations

Magnetic correlations in two-dimensional spin-glasses

I. Morgenstern, Kurt Binder

By a recursive method numerically exact free energies are calculated for square $L\ifmmode\times\else\texttimes\fi{}L$ Ising lattices, with $6<~L<~18$, for several kinds of frozen-in bond disorder: (i) bonds $\ifmmode\pm\else\textpm\fi{}J$ with various concentrations of negative bonds; (ii) bonds distributed according to a Gaussian distribution. Ground states of these systems are identified, the response to "ordering fields" is studied, and the correlation function ${〈{S}_{0}{S}_{R}〉}_{T}^{2}$ is calculated as a function of temperature for various distances $R$ in the lattice. This correlation is found to decay strongly (presumably exponentially) with increasing $R$ even at temperatures distinctly below the apparent freezing temperature ${T}_{f}$ of previous Monte Carlo simulations; this "freezing transition" is hence unambiguously identified as a nonequilibrium effect. However, the correlation length is found to become long ranged at low temperatures, and it is suggested that a phase transition still occurs at $T=0$; while in the Gaussian model the spin-glass order parameter $q(T=0)=1$, it is found that $q\ensuremath{\equiv}0$ in the $\ifmmode\pm\else\textpm\fi{}J$ model where rather a power-law decay of correlations ${〈{S}_{0}{S}_{R}〉}_{T=0}^{2}$ occurs. Performing Monte Carlo simulations for precisely the same systems, the cooling times necessary to reach the true ground states of the system are identified, as well as the simulation times necessary to reach thermal equilibrium for the correlation functions. These times are found to increase so strongly with $L$ that for systems of macroscopic size the correct thermal equilibrium is probably irrelevant for experimental purposes. Rather a statistical mechanics based on the many long-lived metastable states would be required.

Physical review. B, Condensed matter•Jul 1, 1980•304 citations

Ground-State Magnetization in Disordered Systems: Exchange versus Off-Diagonal Interaction Fluctuations

Philippe Jacquod, A. Douglas Stone

In a Hartree-Fock picture, itinerant ferromagnetism results from a competition between kinetic and exchange energy, with the magnetized state being favored at large interaction strength. In a recent paper [1] we showed that in contrast to this average effect, fluctuations of off-diagonal interaction matrix elements systematically reduce the ground-state magnetization. When the interaction dominates, the occurence of a non-zero ground-state magnetization depends on the ratio λ between average exchange energy and the fluctuation amplitude of the off-diagonal matrix elements, and a nonzero critical value λc is necessary to magnetize the ground state. We extend these results by numerical studies of λ for standard tight-binding models which indicate a regime of intermediate disorder where off-diagonal fluctuations should play an important role for ground-state magnetization. We also emphasize the presence of strong correlations between minimal eigenvalues of different magnetization blocks which further reduce the probability of a nonzero ground-state spin.

arXiv:cond-mat/0003352•physica status solidi (b)•Mar 1, 2000•1 citation
We study numerically the ground-state morphology and properties of the low-lying states of an XY model with competing interactions for several values of the concentration x and of the ratio between ferromagnetic and antiferromagnetic exchange constants, \ensuremath{\lambda}. If \ensuremath{\lambda}>1, the ground state is likely to be a spin glass for all nonzero concentrations. In the region \ensuremath{\lambda}\ensuremath{\le}1, the ground state is a disordered canted ferromagnet for low dilutions with a transition to a pure spin-glass phase at a \ensuremath{\lambda}-dependent critical concentration. The analysis of the spin-spin correlation functions gives insight on the mechanism responsible for the degradation of the magnetic order at low dilutions. An unusual law of decay of correlations suggests that the canted ferromagnetic phase that we find at small \ensuremath{\lambda} and x may actually be a phase with power-law decay of correlations.
We investigate the ground state properties of the ± J Ising model with asymmetric probability weight in two dimensions. We propose a numerical transfer matrix method which enables us to calculate directly the ground-state properties such as correlation functions, defect energy and entropy at T =0. We find that the spin glass correlation function decays in a power law like g~(r)∼r−η~\tilde{g}(r)\sim{r}^{-\tilde{\eta}}g~​(r)∼r−η~​ with η~≃0.21\tilde{\eta}\simeq0.21η~​≃0.21 in the range 0.5≤p≲0.850.5\leq{p}\lesssim0.850.5≤p≲0.85, suggesting the weak universality of spin glass-the concentration of + J bond is denoted by p . At the ferromagnetic phase boundary locating around p =0.89, the ferromagnetic correlation function decays as g(r)∼r−η~g(r)\sim{r}^{-\tilde{\eta}}g(r)∼r−η~​ with η≃0.240±0.016, which suggests the weak universality of ferromagnetic critical exponent from the non-random critical point ( p =1, T =2.269···) to the ground sate counterpart. The behavior of correlation functions, and also that of defect energies indicate the existence of an extra-phase, called the random antiphase state, in 0.85≲p≲0.890.85\lesssim{p}\lesssim0.890.85≲p≲0.89.
We introduce a one dimensional Ising model with two competing interactions: nearest neighbor random couplings +/-J with equal probability and a positive infinite range coupling Lambda. At low temperature T the model exhibits a first order phase transition between a ferromagnetic state (with magnetization m(1) = 1 at T = 0) and a << ferrimagnetic >> state (with m(2) = 2/3 at T = 0), when the disorder strength J/Lambda is increased. For 5/12 < J/Lambda < 1, a whole spectrum of ferrimagnetic ground states with magnetization m(n) = 2/(n + 1) (n = 2, ..., infinity) is present while for J/lambda > 1 the ground state is given by a trivial one dimensional spin glass with m = 0. The main qualitative features of the model can be described by a simplified annealed model where the random couplings can arrange themselves to minimize free energy with the constraint that the number of positive couplings is fixed by the law of large numbers in the thermodynamic limit. This model is exactly solved at all temperatures and the diagram of phase is calculated.
We consider the XY model with ferromagnetic (FM) and antinematic (AN) nearest-neighbor interactions on a square lattice for a varying interaction strength ratio. Besides the expected FM and AN quasi-long-range order (QLRO) phases we identify at low temperatures another peculiar canted ferromagnetic (CFM) QLRO phase, resulting from the competition between the collinear FM and noncollinear AN ordering tendencies. In the CFM phase neighboring spins that belong to different sublattices are canted by a nonuniversal (dependent on the interaction strength ratio) angle and the ordering is characterized by a fast-decaying power-law intrasublattice correlation function. Compared to the FM phase, in the CFM phase correlations are significantly diminished by the presence of zero-energy domain walls due to the inherent degeneracy caused by the AN interactions. We present the phase diagram as a function of the interaction strength ratio and discuss the character of the respective phases as well as the transitions between them.
The spin-spin correlation function of the two-dimensional (2D) XYXYXY model decays as a power law at all temperatures below the Berezinskii-Kosterlitz-Thouless transition point with a temperature-dependent exponent $\ensuremath{\eta}=\ensuremath{\eta}(T/J)$ (JJJ is the ferromagnetic coupling strength). It is known from computer experiments that in the 2D XYXYXY model with site or bond dilution this exponent depends on the concentration ppp of removed sites or bonds as well. Knowing the slope $\ensuremath{\partial}\ensuremath{\eta}/\ensuremath{\partial}p$ at point p=0p=0p=0, one can predict the value of the exponent for small dilution concentrations: ${\ensuremath{\eta}(p)\ensuremath{\simeq}\ensuremath{\eta}(0)+p(\ensuremath{\partial}\ensuremath{\eta}/\ensuremath{\partial}p)|}_{p=0}$. As shown in this paper, the spin-wave Hamiltonian allows one to obtain exact results for this slope: ${(\ensuremath{\partial}\ensuremath{\eta}/\ensuremath{\partial}p)|}_{p=0}=T/(2J)+O({(T/J)}^{2})$ and $T/(\ensuremath{\pi}J)+O({(T/J)}^{2})$ for site and for bond dilution, respectively.
By a recursive method numerically exact free energies are calculated for square $L\ifmmode\times\else\texttimes\fi{}L$ Ising lattices, with 6< L< 186<~L<~186< L< 18, for several kinds of frozen-in bond disorder: (i) bonds $\ifmmode\pm\else\textpm\fi{}J$ with various concentrations of negative bonds; (ii) bonds distributed according to a Gaussian distribution. Ground states of these systems are identified, the response to "ordering fields" is studied, and the correlation function 〈S0SR〉T2{〈{S}_{0}{S}_{R}〉}_{T}^{2}〈S0​SR​〉T2​ is calculated as a function of temperature for various distances RRR in the lattice. This correlation is found to decay strongly (presumably exponentially) with increasing RRR even at temperatures distinctly below the apparent freezing temperature Tf{T}_{f}Tf​ of previous Monte Carlo simulations; this "freezing transition" is hence unambiguously identified as a nonequilibrium effect. However, the correlation length is found to become long ranged at low temperatures, and it is suggested that a phase transition still occurs at T=0T=0T=0; while in the Gaussian model the spin-glass order parameter q(T=0)=1q(T=0)=1q(T=0)=1, it is found that $q\ensuremath{\equiv}0$ in the $\ifmmode\pm\else\textpm\fi{}J$ model where rather a power-law decay of correlations 〈S0SR〉T=02{〈{S}_{0}{S}_{R}〉}_{T=0}^{2}〈S0​SR​〉T=02​ occurs. Performing Monte Carlo simulations for precisely the same systems, the cooling times necessary to reach the true ground states of the system are identified, as well as the simulation times necessary to reach thermal equilibrium for the correlation functions. These times are found to increase so strongly with LLL that for systems of macroscopic size the correct thermal equilibrium is probably irrelevant for experimental purposes. Rather a statistical mechanics based on the many long-lived metastable states would be required.
In a Hartree-Fock picture, itinerant ferromagnetism results from a competition between kinetic and exchange energy, with the magnetized state being favored at large interaction strength. In a recent paper [1] we showed that in contrast to this average effect, fluctuations of off-diagonal interaction matrix elements systematically reduce the ground-state magnetization. When the interaction dominates, the occurence of a non-zero ground-state magnetization depends on the ratio λ between average exchange energy and the fluctuation amplitude of the off-diagonal matrix elements, and a nonzero critical value λc is necessary to magnetize the ground state. We extend these results by numerical studies of λ for standard tight-binding models which indicate a regime of intermediate disorder where off-diagonal fluctuations should play an important role for ground-state magnetization. We also emphasize the presence of strong correlations between minimal eigenvalues of different magnetization blocks which further reduce the probability of a nonzero ground-state spin.