Monte Carlo determination of the low-energy constants for a two-dimensional spin-1 Heisenberg model with spatial anisotropy

The low-energy constants, namely the spin stiffness ρ _s , the staggered magnetization density ? _s per area, and the spinwave velocity c of the two-dimensional (2D) spin-1 Heisenberg model on the square and rectangular lattices are determined using the first principles Monte Carlo method. In particular, the studied models have different antiferromagnetic couplings J ₁ and J ₂ in the spatial 1- and 2-directions, respectively. For each considered J ₂∕J ₁, the aspect ratio of the corresponding linear box sizes L ₂∕L ₁ used in the simulations is adjusted so that the squares of the two spatial winding numbers take the same values. In addition, the relevant finite-volume and -temperature predictions from magnon chiral perturbation theory are employed in extracting the numerical values of these low-energy constants. Our results of ρ _s1 are in quantitative agreement with those obtained by the series expansion method over a broad range of J ₂∕J ₁. This in turn provides convincing numerical evidence for the quantitative correctness of our approach. The ? _s and c presented here for the spatially anisotropic models are new and can be used as benchmarks for future related studies.     

BCS-BEC crossover on the two-dimensional honeycomb lattice

The attractive Hubbard model on the honeycomb lattice exhibits, at half filling, a quantum critical point between a semimetal with massless Dirac fermions and an s-wave superconductor (SC). We study the BCS-BEC crossover in this model away from half filling at zero temperature and show that the appropriately defined crossover line (in the interaction-density plane) passes through the quantum critical point at half filling. For a range of densities around half filling, the "underlying Fermi surface" of the SC, defined as the momentum space locus of minimum energy quasiparticle excitations, encloses an area which changes nonmonotonically with interaction. We also study fluctuations in the SC and the semimetal, and show the emergence of an undamped Leggett mode deep in the SC. Finally, we consider possible implications for ultracold atoms in optical lattices and the high temperature SCs.     

Chiral extrapolation and determination of low-energy constants from lattice data

We propose analytic approximations of chiral SU(3)$\mathit{SU}(3)$ amplitudes for the extrapolation of lattice data to the physical meson masses. The method allows the determination of NNLO low-energy constants in a controllable fashion. We test the approach with recent lattice data for the ratio FK/Fπ$F_K/F_{\pi }$ of meson decay constants.     

Spin excitations and thermodynamics of the antiferromagnetic Heisenberg model on the layered honeycomb lattice

We present a spin-rotation-invariant Green-function theory for the dynamic spin susceptibility in the spin-1/2 antiferromagnetic Heisenberg model on a stacked honeycomb lattice. Employing a generalized mean-field approximation for arbitrary temperatures, the thermodynamic quantities (two-spin correlation functions, internal energy, magnetic susceptibility, staggered magnetization, Néel temperature, correlation length) and the spin-excitation spectrum are calculated by solving a coupled system of self-consistency equations for the correlation functions. The temperature dependence of the magnetic (uniform static) susceptibility is ascribed to antiferromagnetic short-range order. The Néel temperature is calculated for arbitrary interlayer couplings. Our results are in a good agreement with numerical computations for finite clusters and with available experimental data on the β-Cu₂V₂O₂ compound.     

The frustrated Heisenberg antiferromagnet on the honeycomb lattice: J1-J2 model

We study the ground-state phase diagram of the frustrated spin-[Formula: see text] antiferromagnet with J(2) = xJ(1) > 0 (J(1) > 0) on the honeycomb lattice, using the coupled-cluster method. We present results for the ground-state energy, magnetic order parameter and plaquette valence-bond crystal (PVBC) susceptibility. We find a paramagnetic PVBC phase for x(c(1)) < x < x(c(2)), where x(c(1)) ≈ 0.207 ± 0.003 and x(c(2)) ≈ 0.385 ± 0.010. The transition at x(c(1)) to the Néel phase seems to be a continuous deconfined transition (although we cannot exclude a very narrow intermediate phase in the range 0.21 ? x ? 0.24), while that at x(c(2)) is of first-order type to another quasiclassical antiferromagnetic phase that occurs in the classical version of the model only at the isolated and highly degenerate critical point [Formula: see text]. The spiral phases that are present classically for all values x > 1/6 are absent for all x ? 1.     

Exotic entanglement scaling of Heisenberg antiferromagnet on honeycomb lattice

The scaling behaviors of entanglement entropy (EE) against dimension cut-off of density matrix renormalization group (DMRG) in an anisotropic Heisenberg model on honeycomb lattice are investigated. In the gapped dimer phase, the entanglement spectrum (ES) exhibits large gaps and the EE shows an unexpected linear scaling before convergence. In contrast in the gapless Néel phase, the ES decays in a much smoother way, and the EE scales logarithmically. Our calculations show that the linear scaling in the dimer phase originates from one dominant Schmidt number plus n (nearly) degenerate Schmidt numbers that are much smaller than the dominant one. The non-trivial entanglement-scaling properties of the dimer and Néel phases could potentially be used for their detections.     

Blume-Emery-Griffiths model on the honeycomb lattice

Exact results are obtained for a spin-1 system on the honeycomb lattice with the Blume-Emery-Griffiths Hamiltonian –/kT =J _i,j S _i S _j +Ki,jS _i ² _j ² – _i S _i ² +HS _i subject to the constraintK=–ln coshJ. ForJ>0, the system behaves like a spin-1/2 Ising ferromagnet with the free energy analytic everywhere except at the first-order phase boundaryH=0, tanhJ<(2+e )/ . Derivatives of the free energy across this boundary are discontinuous and we obtain the exact expression for the spontaneous magnetization. ForJ<0, the system can be transcribed into an antiferromagnetic spin-1/2 Ising model in a real magnetic field, and from this equivalence portions of the exact phase boundary are determined.     

Ising transition in the two-dimensional quantum J1-J2 Heisenberg model

We study the thermodynamics of the spin-S two-dimensional quantum Heisenberg antiferromagnet on the square lattice with nearest (J1) and next-nearest (J2) neighbor couplings in its collinear phase (J(2)/J(1)>0.5), using the pure-quantum self-consistent harmonic approximation. Our results show the persistence of a finite-temperature Ising phase transition for every value of the spin, provided that the ratio J(2)/J(1) is greater than a critical value corresponding to the onset of collinear long-range order at zero temperature. We also calculate the spin and temperature dependence of the collinear susceptibility and correlation length, and we discuss our results in light of the experiments on Li2VOSiO4 and related compounds.     

Evidence for deconfined quantum criticality in a two-dimensional Heisenberg model with four-spin interactions

Using ground-state projector quantum Monte Carlo simulations in the valence-bond basis, it is demonstrated that nonfrustrating four-spin interactions can destroy the Néel order of the two-dimensional S=1/2 Heisenberg antiferromagnet and drive it into a valence-bond solid (VBS) phase. Results for spin and dimer correlations are consistent with a single continuous transition, and all data exhibit finite-size scaling with a single set of exponents, z=1, nu=0.78+/-0.03, and eta=0.26+/-0.03. The unusually large eta and an emergent U(1) symmetry, detected using VBS order parameter histograms, provide strong evidence for a deconfined quantum critical point.     

Universal effective coupling constants for the generalized Heisenberg model

The aim of this study is to find universal critical values of the effective dimensionless coupling constant g ₆ and refined universal values g ₄ for Heisenberg ferromagnets with n-component order parameters. These constants appear in the equation of state and determine the nonlinear susceptibilities χ ₄ and χ ₆ in the critical region. Calculations are made of the first three terms of the expansion of g ₆ in powers of g ₄ in the limits of O(n) symmetry three-dimensional λϕ ⁴ theory, the resultant series is resummed by the Padé-Borel method, and then by substituting the fixed point coordinates g ₄ ^* in the resultant expression, numerical values of g ₆ ^* are obtained for different n. These numbers g ₄ ^* for n>3 were determined from a six-loop expansion for the β-function resummed using the Padé-Borel-Leroy technique. An analysis of the accuracy of these g ₆ ^* values showed that they may differ from the true values by no more than 1.6%. These values of g ₆ ^* were compared with those obtained by the 1/n expansion method which allowed the level of accuracy of this method to be assessed. Fiz. Tverd. Tela (St. Petersburg) 40, 1284–1290 (July 1998)     

Half-filled Kondo lattice on the honeycomb lattice

The unique linear density of state around the Dirac points for the honeycomb lattice brings much novel features in strongly correlated models. Here we study the ground-state phase diagram of the Kondo lattice model on the honeycomb lattice at half-filling by using an extended mean-field theory. By treating magnetic interaction and Kondo screening on an equal footing, it is found that besides a trivial discontinuous first-order quantum phase transition between well-defined Kondo insulator and antiferromagnetic insulating state, there can exist a wide coexistence region with both Kondo screening and antiferromagnetic orders in the intermediate coupling regime. In addition, the stability of Kondo insulator requires a minimum strength of the Kondo coupling. These features are attributed to the linear density of state, which are absent in the square lattice. Furthermore, fluctuation effect beyond the mean-field decoupling is analyzed and the corresponding antiferromagnetic spin-density-wave transition falls into the O(3) universal class. Comparatively, we also discuss the Kondo necklace and the Kane-Mele-Kondo (KMK) lattice models on the same lattice. Interestingly, it is found that the topological insulating state is unstable to the usual antiferromagnetic ordered states at half-filling for the KMK model. The present work may be helpful for further study on the interplay between conduction electrons and the densely localized spins on the honeycomb lattice.     

An investigation of the quantum J 1 - J 2 - J 3 model on the honeycomb lattice

We have investigated the quantum J ₁ - J ₂ - J ₃ model on the honeycomb lattice with exact diagonalizations and linear spin-wave calculations for selected values of J ₂ / J ₁ , J ₃ / J ₁ and antiferromagnetic (J ₁ > 0) or ferromagnetic (J ₁ < 0) nearest neighbor interactions. We found a variety of quantum effects: “order by disorder" selection of a Néel ordered ground-state, good candidates for non-classical ground-states with dimer long range order or spin-liquid like. The purely antiferromagnetic Heisenberg model is confirmed to be Néel ordered. Comparing these results with those observed on the square and triangular lattices, we enumerate some conjectures on the nature of the quantum phases in the isotropic models. Received 17 November 2000 and Received in final form 21 January 2001     

Physical realization of topological excitations in quantum Heisenberg ferromagnet on lattice

The physical spin configurations corresponding to the topological excitations expected to be present in the XY limit of a purely quantum spin \frac 12\frac {1}{2} Heisenberg ferromagnet, are investigated on a two dimensional square lattice. Quantum vortices (anti-vortices) are constructed by making use of the coherent spin field components from meronic (anti-meronic) configurations in the limiting case of very strong XY anisotropy. The equations for pseudo-time evolution of coherent spin fields used in this analysis, are obtained by applying variational principle on the quantum Euclidean action corresponding to the Heisenberg ferromagnet on lattice. The important role of the associated topological-like term extracted from the Wess-Zumino-like contribution, is highlighted in our procedure. It is shown that this term can identify a large class of vortices (anti-vortices). In particular the excitations having odd topological charges form this class. Our present formalism is markedly different from the conventional approach for the construction of quantum vortices (anti-vortices).