By using a simulated annealing approach, Monte Carlo and molecular-dynamics techniques we have studied static and dynamic behavior of the classical two-dimensional anisotropic Heisenberg model. We have obtained numerically that the vortex developed in such a model exhibit two different behaviors depending if the value of the anisotropy $lambda$ lies below or above a critical value $lambda_c$ . The in-plane and out-of-plane correlation functions ($S^{xx}$ and $S^{zz}$) were obtained numerically for $lambda < lambda_c$ and $lambda > lambda_c$ . We found that the out-of-plane dynamical correlation function exhibits a central peak for $lambda > lambda_c$ but not for $lambda < lambda_c$ at temperatures above $T_{BKT}$ .
Ultrathin magnetic films can be modeled as an anisotropic Heisenberg model with long-range dipolar interactions. It is believed that the phase diagram presents three phases: An ordered ferromagnetic phase I, a phase characterized by a change from out-of-plane to in-plane in the magnetization II, and a high-temperature paramagnetic phase III. It is claimed that the border lines from phase I to III and II to III are of second order and from I to II is first order. In the present work we have performed a very careful Monte Carlo simulation of the model. Our results strongly support that the line separating phases II and III is of the BKT type.
We study classical and quantum Heisenberg antiferromagnets with exchange anisotropy of XXZ-type and crystal field single-ion terms of quadratic and cubic form in a field. The magnets display a variety of phases, including the spin-flop (or, in the quantum case, spin-liquid) and biconical (corresponding, in the quantum lattice gas description, to supersolid) phases. Applying ground-state considerations, Monte Carlo and density matrix renormalization group methods, the impact of quantum effects and lattice dimension is analysed. Interesting critical and multicritical behaviour may occur at quantum and thermal phase transitions.
The classical Heisenberg model in two spatial dimensions constitutes one of the most paradigmatic spin models, taking an important role in statistical and condensed matter physics to understand magnetism. Still, despite its paradigmatic character and the widely accepted ban of a (continuous) spontaneous symmetry breaking, controversies remain whether the model exhibits a phase transition at finite temperature. Importantly, the model can be interpreted as a lattice discretization of the $O(3)$ non-linear sigma model in $1+1$ dimensions, one of the simplest quantum field theories encompassing crucial features of celebrated higher-dimensional ones (like quantum chromodynamics in $3+1$ dimensions), namely the phenomenon of asymptotic freedom. This should also exclude finite-temperature transitions, but lattice effects might play a significant role in correcting the mainstream picture. In this work, we make use of state-of-the-art tensor network approaches, representing the classical partition function in the thermodynamic limit over a large range of temperatures, to comprehensively explore the correlation structure for Gibbs states. By implementing an $SU(2)$ symmetry in our two-dimensional tensor network contraction scheme, we are able to handle very large effective bond dimensions of the environment up to $chi_E^text{eff} sim 1500$, a feature that is crucial in detecting phase transitions. With decreasing temperatures, we find a rapidly diverging correlation length, whose behaviour is apparently compatible with the two main contradictory hypotheses known in the literature, namely a finite-$T$ transition and asymptotic freedom, though with a slight preference for the second.
By considering the appropriate finite-size effect, we explain the connection between Monte Carlo simulations of two-dimensional anisotropic Heisenberg antiferromagnet in a field and the early renormalization group calculation for the bicritical point in $2+epsilon$ dimensions. We found that the long length scale physics of the Monte Carlo simulations is indeed captured by the anisotropic nonlinear $sigma$ model. Our Monte Carlo data and analysis confirm that the bicritical point in two dimensions is Heisenberg-like and occurs at T=0, therefore the uncertainty in the phase diagram of this model is removed.
The existence of Neel order in the S=1/2 Heisenberg model on the square lattice at T=0 is shown using inequalities set up by Kennedy, Lieb and Shastry in combination with high precision Quantum Monte Carlo data.