Heat jet approach for finite temperature atomic simulations of triangular lattice

2017 ◽  
Vol 59 (5) ◽  
pp. 843-859 ◽  
Author(s):  
Baiyili Liu ◽  
Shaoqiang Tang ◽  
Jun Chen
2001 ◽  
Vol 79 (11-12) ◽  
pp. 1459-1461 ◽  
Author(s):  
W Stephan ◽  
B W Southern

The phase diagram of the classical anisotropic (XXZ) Heisenberg model on the two-dimensional triangular lattice is investigated using Monte Carlo methods. In the easy-axis limit, two finite-temperature vortex-unbinding transitions have been observed. In the easy-plane limit, there also appear to be two distinct finite-temperature phase transitions that are very close in temperature. The upper transition corresponds to an Ising-like chirality ordering and the lower temperature transition corresponds to a Kosterlitz–Thouless vortex-unbinding transition. These phase-transition lines all meet at the Heisenberg point and provide strong evidence that the isotropic model undergoes a novel finite-temperature phase transition. PACS Nos.: 75.10Hk, 75.40Mg


2013 ◽  
Vol 91 (7) ◽  
pp. 542-547 ◽  
Author(s):  
Solomon A. Owerre

We present the linear spin wave theory calculation of the superfluid phase of a hard-core boson J-K model with nearest neighbour exchange J and four-particle ring-exchange K at half filling on the triangular lattice, as well as the phase diagrams of the system at zero and finite temperatures. A similar analysis has been done on a square lattice (Schaffer et al. Phys. Rev. B, 80, 014503 (2009)). We find similar behaviour to that of a square lattice but with different spin wave values of the thermodynamic quantities. We also find that the pure J model (XY model), which has a well-known uniform superfluid phase with an ordered parameter [Formula: see text] at zero temperature is quickly destroyed by the inclusion of negative-K ring-exchange interactions, favouring a state with a (4π/3, 0) ordering wavevector. We further study the behaviour of the finite-temperature Kosterlitz–Thouless phase transition (TKT) in the uniform superfluid phase, by forcing the universal quantum jump condition on the finite-temperature spin wave superfluid density. We find that for K < 0, the phase boundary monotonically decreases to T = 0 at K/J = −4/3, where a phase transition is expected and TKT decreases rapidly, while for positive K, TKT reaches a maximum at some K ≠ 0. It has been shown on a square lattice using quantum Monte Carlo (QMC) simulations that for small K > 0 away from the XY point, the zero-temperature spin stiffness value of the XY model is decreased (Melko and Sandvik. Ann. Phys. 321, 1651 (2006)). Our result seems to agree with this trend found in QMC simulations for two-dimensional systems.


2015 ◽  
Vol 18 (5) ◽  
pp. 1445-1460 ◽  
Author(s):  
Shaoqiang Tang ◽  
Baiyili Liu

AbstractIn this paper, we propose a heat jet approach for atomic simulations at finite temperature. Thermal fluctuations are injected into an atomic subsystem from its boundaries, without modifying the governing equations for the interior domain. More precisely, we design a two way local boundary condition, and take the incoming part of a phonon representation for thermal fluctuation input. In this way, nonthermal wave propagation simulations are effectively performed at finite temperature. We further apply this approach to nonlinear chains with the Morse potential. Chains with model parameters fitted to carbon and gold are simulated at room temperature with fidelity.


2019 ◽  
Vol 21 (7) ◽  
pp. 073015 ◽  
Author(s):  
Shengjie Jin ◽  
Xinxin Guo ◽  
Peng Peng ◽  
Xuzong Chen ◽  
Xiaopeng Li ◽  
...  

2006 ◽  
Vol 17 (04) ◽  
pp. 591-600
Author(s):  
KWANGHOON CHUNG ◽  
MOOKYUNG CHEON ◽  
IKSOO CHANG

The critical coarsening dynamics of the spin S =1/2, 3/2 antiferromagnetic Ising model on a triangular lattice is studied by the dynamic Monte Carlo simulation using a heat bath algorithm. The triangular antiferromagnetic Ising (TAI) model possesses an intrinsic geometrical frustration and a large degeneracy of ground state which may affect the equilibrium and non-equilibrium critical behaviors. The S =1/2 TAI has no phase transition at a finite temperature, but it was suggested that the S =3/2 TAI has the Kosterlitz–Thouless (KT)-type phase transition at a finite temperature. We employ a finite size scaling approach for the correlation function from the short-time dynamics of the S =1/2, 3/2 TAI to calculate the values of the static critical exponent η and the dynamic exponent z. For the S =1/2 TAI, our dynamic scaling analysis provides η =0.498±0.006 and z =2.278±0.020 at T =0, agreeing with the previous results. For the S =3/2 TAI, after identifying a KT-transition temperature TKT =0.51±0.01, we find the temperature-dependent η ranging from 0.301±0.008 at T =0.51 to 0.224±0.016 at T =0 along the KT-line whereas the value of z =2.20±0.06 is constant for T≤TKT. It is shown that the spin S =3/2 TAI model and the two-dimensional XY model, sharing the KT-type phase transition, exhibit similar static critical and coarsening dynamics behavior. For both the S =1/2, 3/2 TAI, the value of z (η) is compatible with (larger than) that of the Ising model at Tc and the XY model for T≤TKT in two-dimension. Our results imply that although the quasi-long-range order disappears with ηXY =0 in the two-dimensional XY model at T =0, the S =3/2 TAI still maintains it with η TAI =0.224 due to the effect of a frustration and a high degeneracy of ground state.


Sign in / Sign up

Export Citation Format

Share Document