The magnetic properties of the two-dimensional S = 1 2 (quantum) antiferromagnetic Heisenberg model on a honeycomb lattice with and without interlayer coupling are studied by means of a continuous Euclidean time Quantum-Monte-Carlo algorithm.The internal energy, the magnetic susceptibility and the staggered magnetization are determined in the full temperature range.For the two-dimensional system the groundstate energy/bond is found to be E hc 0 = -0.36303(13),and the zero temperature staggered magnetization m st = 0.2681(8).For coupled planes of honeycomb systems a phase transition from an ordered phase to a disordered phase is found at T /J = 0.695(10).
The magnetic properties of the two-dimensional S = 1 2 (quantum) antiferromagnetic Heisenberg model on a honeycomb lattice with and without interlayer coupling are studied by means of a continuous Euclidean time Quantum-Monte-Carlo algorithm.The internal energy, the magnetic susceptibility and the staggered magnetization are determined in the full temperature range.For the two-dimensional system the groundstate energy/bond is found to be E hc 0 = -0.36303(13),and the zero temperature staggered magnetization m st = 0.2681(8).For coupled planes of honeycomb systems a phase transition from an ordered phase to a disordered phase is found at T /J = 0.695(10).
The continuous imaginary-time quantum Monte Carlo method with the worm update algorithm is applied to explore the ground-state properties of the spin-1/2 Heisenberg model with antiferromagnetic (AF) coupling J>0 and ferromagnetic (F) coupling J^{'}<0 along zigzag and armchair directions, respectively, on honeycomb lattice. It is found that by enhancing the F coupling J^{'} between zigzag AF chains, the system is smoothly crossover from one-dimensional zigzag spin chains to a two-dimensional magnetic ordered state. In absence of an external field, the system is in a stripe-ordered phase. In the presence of uniform and staggered fields, the uniform and staggered out-of-plane magnetizations appear while the stripe order remains in the xy plane, and a second-order quantum phase transition (QPT) at a critical staggered field is observed. The critical exponents of correlation length for QPTs induced by a staggered field for the cases with J>0, J^{'}<0 and J<0, J^{'}>0 are obtained to be ν=0.70046(1) and 0.7086(3), respectively, indicating that both cases belong to O(3) universality. The corresponding dynamic and susceptibility exponent z and γ/ν are fitted to be 1.006572(9), 1.9412(2) and 1.004615(8), 1.96121(9) for the two cases, respectively. The scaling behavior in a staggered field is analyzed, and the ground-state phase diagrams in the plane of coupling ratio and staggered field are presented for two cases. The temperature dependence of susceptibility and specific heat of both systems in external magnetic fields is also discussed. A Kosterlitz-Thouless phase transition is found for the present system in a uniform field.
We consider the spin-1/2 antiferromagnetic Heisenberg model on a bilayer honeycomb lattice including interlayer frustration in the presence of an external magnetic field. In the vicinity of the saturation field, we map the low-energy states of this quantum system onto the spatial configurations of hard hexagons on a honeycomb lattice. As a result, we can construct effective classical models (lattice-gas as well as Ising models) on the honeycomb lattice to calculate the properties of the frustrated quantum Heisenberg spin system in the low-temperature regime. We perform classical Monte Carlo simulations for a hard-hexagon model and adopt known results for an Ising model to discuss the finite-temperature order-disorder phase transition that is driven by a magnetic field at low temperatures. We also discuss an effective-model description around the ideal frustration case and find indications for a spin-flop-like transition in the considered isotropic spin model.
We study the zero-temperature quantum phase diagram of a spin-1 Heisenberg antiferromagnet on the honeycomb lattice with both nearest-neighbor exchange coupling J1>0 and frustrating next-nearest-neighbor coupling ${J}_{2}\ensuremath{\equiv}\ensuremath{\kappa}{J}_{1}>0$, using the coupled cluster method implemented to high orders of approximation, and based on model states with different forms of classical magnetic order. For each we calculate directly in the bulk thermodynamic limit both ground-state low-energy parameters (including the energy per spin, magnetic order parameter, spin stiffness coefficient, and zero-field uniform transverse magnetic susceptibility) and their generalized susceptibilities to various forms of valence-bond crystalline (VBC) order, as well as the energy gap to the lowest-lying spin-triplet excitation. In the range $0<\ensuremath{\kappa}<1$ we find evidence for four distinct phases. Two of these are quasiclassical phases with antiferromagnetic long-range order, one with two-sublattice N\'eel order for $\ensuremath{\kappa}<{\ensuremath{\kappa}}_{{c}_{1}}=0.250(5)$, and another with four-sublattice N\'eel-II order for $\ensuremath{\kappa}>{\ensuremath{\kappa}}_{{c}_{2}}=0.340(5)$. Two different paramagnetic phases are found to exist in the intermediate region. Over the range ${\ensuremath{\kappa}}_{{c}_{1}}<\ensuremath{\kappa}<{\ensuremath{\kappa}}_{c}^{i}=0.305(5)$ we find a gapless phase with no discernible magnetic order, which is a strong candidate for being a quantum spin liquid, while over the range ${\ensuremath{\kappa}}_{c}^{i}<\ensuremath{\kappa}<{\ensuremath{\kappa}}_{{c}_{2}}$ we find a gapped phase, which is most likely a lattice nematic with staggered dimer VBC order that breaks the lattice rotational symmetry.
The quantum phase transition, scaling behaviors, and thermodynamics in the spin-1/2 quantum Heisenberg model with antiferromagnetic coupling J>0 in the armchair direction and ferromagnetic interaction J^{'}<0 in the zigzag direction on a honeycomb lattice are systematically studied using the continuous-time quantum Monte Carlo method. By calculating the Binder ratio Q_{2} and spin stiffness ρ in two directions for various coupling ratios α=J^{'}/J under different lattice sizes, we found that a quantum phase transition from the dimerized phase to the stripe phase occurs at the quantum critical point α_{c}=-0.93. Through the finite-size scaling analysis on Q_{2}, ρ_{x}, and ρ_{y}, we determined the critical exponent related to the correlation length ν to be 0.7212(8), implying that this transition falls into a classical Heisenberg O(3) universality. A zero magnetization plateau is observed in the dimerized phase, whose width decreases with increasing α. A phase diagram in the coupling ratio α-magnetic field h plane is obtained, where four phases, including dimerized, stripe, canted stripe, and polarized, are identified. It is also unveiled that the temperature dependence of the specific heat C(T) for different α's intersects precisely at one point, similar to that of liquid ^{3}He under different pressures and several magnetic compounds under various magnetic fields. The scaling behaviors of Q_{2}, ρ, and C(T) are carefully analyzed. The susceptibility is compared with the experimental data to give the magnetic parameters of both compounds.
We use a combination of analytical and numerical techniques to study the phase diagram of the frustrated Heisenberg model on the bilayer honeycomb lattice. Using the Schwinger-boson description of the spin operators followed by a mean-field decoupling, the magnetic phase diagram is studied as a function of the frustration coupling J2 and the interlayer coupling ${J}_{\ensuremath{\perp}}$. The presence of both magnetically ordered and disordered phases is investigated by means of the evaluation of ground-state energy, spin gap, local magnetization, and spin-spin correlations. We observe a phase with a spin gap and short-range N\'eel correlations that survives for nonzero next-nearest-neighbor interaction and interlayer coupling. Furthermore, we detect signatures of a reentrant behavior in the melting of the N\'eel phase and symmetry restoring when the system undergoes a transition from an on-layer nematic valence-bond crystal phase to an interlayer valence-bond crystal phase. We complement our work with exact diagonalization on small clusters and dimer-series expansion calculations, together with a linear spin-wave approach to study the phase diagram as a function of the spin S, the frustration, and the interlayer couplings.
The ground-state energy, staggered magnetization, and correlation function are computed for the two-dimensional spin-(1/2 quantum antiferromagnetic Heisenberg model (believed to be relevant to the high-temperature superconducting materials without doping) on lattices up to 24\ifmmode\times\else\texttimes\fi{}24 by Monte Carlo methods extrapolated to zero temperature. A more accurate value for the ground-state energy is determined by the method of Neumann and Ulam by extrapolating from lattices up to 12\ifmmode\times\else\texttimes\fi{}12.