We study the ground-state phase diagrams and properties of spin-1/2 Heisenberg models on the diamond-like-decorated square and triangular lattices. The diamond-like-decorated square (triangular) lattice is a lattice in which the bonds of a square (triangular) lattice are replaced with diamond units. The diamond unit has two types of antiferromagnetic exchange interactions, and the ratio λ of the strength of the diagonal bond to that of the other four edges determines the ground-state properties. In particular, the macroscopically degenerated tetramer-dimer states, which are equivalent to the dimer covering states of the original lattices, are stabilized for λc < λ < 2, where the value of λc depends on the lattices. To determine the phase diagrams and boundaries λc, we employ the modified spin wave (MSW) method and the quantum Monte Carlo (QMC) method to estimate the ground-state energies of the ferrimagnetic states for λ < λc, where we can consider the mixed spin-1 and spin-1/2 Lieb-lattice and triangular Lieb-lattice Heisenberg antiferromagnets instead, and obtain λc(square)=0.974 and λc(triangular)=0.988. We also calculate the long-range order (LRO) parameters using the MSW and QMC methods and find the scaling relations where the spin reductions of each sublattice are inversely proportional to the number of sublattice sites. We prove these scaling relations by applying an infinitesimal uniform magnetic field. Furthermore, by examining the calculation process in the MSW, we clarify the mathematical structure behind the scaling relations for the sublattice LROs.
The statistics of the ground-state and domain-wall energies for the two-dimensional random-bond Ising model on square lattices with independent, identically distributed bonds having probability
p of ${J}_{ij}=\ensuremath{-}1$ and $(1\ensuremath{-}p)$ of
Jij=+1 are studied. We are able to consider large samples of up to
3202 spins by using sophisticated matching algorithms. We study $L\ifmmode\times\else\texttimes\fi{}L$ systems, but we also consider $L\ifmmode\times\else\texttimes\fi{}M$ samples, for different aspect ratios
R=L∕M. We find that the scaling behavior of the ground-state energy and its sample-to-sample fluctuations inside the spin-glass region $({p}_{c}\ensuremath{\leqslant}p\ensuremath{\leqslant}1\ensuremath{-}{p}_{c})$ are characterized by simple scaling functions. In particular, the fluctuations exhibit a cusplike singularity at
pc. Inside the spin-glass region, the average domain-wall energy converges to a finite nonzero value as the sample size becomes infinite, holding
R fixed. Here, large finite-size effects are visible, which can be explained for all
p by a single exponent $\ensuremath{\omega}\ensuremath{\approx}2∕3$, provided higher-order corrections to scaling are included. Finally, we confirm the validity of aspect-ratio scaling for $R\ensuremath{\rightarrow}0$: the distribution of the domain-wall energies converges to a Gaussian for $R\ensuremath{\rightarrow}0$, although the domain walls of neighboring subsystems of size $L\ifmmode\times\else\texttimes\fi{}L$ are not independent.
We study, on the basis of the general entangled-plaquette variational ansatz, the ground-state (GS) properties of the spin- antiferromagnetic Heisenberg model on the triangular lattice. Our numerical estimates are in good agreement with available exact results and comparable, for large system sizes, to those computed via the best alternative numerical approaches, or by means of variational schemes based on specific (i.e. incorporating problem-dependent terms) trial wave functions (WFs). The extrapolation to the thermodynamic limit of our results for lattices comprising up to N =324 spins yields an upper bound of the GS energy per site (in units of the exchange coupling) of − 0.5458(2) [−0.4074(1) for the XX model], while the estimated infinite-lattice order parameter is 0.3178(5) (i.e. ∼64% of the classical value).
Various lattice geometries and boundary conditions are used to investigate valence-bond-solid (VBS) ordering in the ground state of an
S=1/2 square-lattice quantum spin model---the
J-
Q model, in which four- or six-spin interactions
Q are added to the standard Heisenberg exchange
J. Ground-state results for finite systems (with up to thousands of spins) are obtained using an unbiased projector quantum Monte Carlo method. It is found that great care has to be taken when extrapolating the order parameter to infinite lattice size, in particular, in cylinder geometry. Even though strong VBS order exists in two dimensions, and is established clearly with increasing system size on $L\ifmmode\times\else\texttimes\fi{}L$ lattices (or ${L}_{x}\ifmmode\times\else\texttimes\fi{}{L}_{y}$ lattices with a fixed aspect ratio
Lx/Ly of order 1), only short-range VBS correlations are observed on long cylinders (when ${L}_{x}\ensuremath{\rightarrow}\ensuremath{\infty}$ at fixed
Ly). The correlation length increases with the cylinder width, until long-range order sets in at a ``critical'' width. This width is very large even when the 2D order is relatively strong. For example, for a system in which the order parameter is
70 of the largest possible value,
Ly=8 is required for ordering. Extrapolations of the VBS order parameter based on correlation functions (the square of the order parameter) for small $L\ifmmode\times\else\texttimes\fi{}L$ lattices can also be misleading. For a
20-ordered system, results for
L up to $\ensuremath{\approx}$20 appear to extrapolate clearly to a vanishing order parameter, while for larger lattices the scaling behavior crosses over and extrapolates to a nonzero value (with exponentially small finite-size corrections). The VBS order parameter also exhibits interesting edge effects related the known emergent U(1) symmetry close to a ``deconfined'' critical point, which, if not considered properly, can lead to wrong conclusions for the thermodynamic limit. The observed finite-size behavior for small $L\ifmmode\times\else\texttimes\fi{}L$ lattices and long cylinders is very similar to that predicted for a
Z2 spin liquid. The results therefore raise concerns about recent numerical work claiming
Z2 spin-liquid ground states in 2D frustrated quantum spin systems, in particular, the Heisenberg model with nearest and next-nearest-neighbor couplings. Based on the results presented so far, a VBS state in this system cannot be ruled out.
We consider a one-dimensional lattice of Ising-type variables where the ferromagnetic exchange interaction
J between neighboring sites is frustrated by a long-ranged antiferromagnetic interaction of strength
g between the sites
i and
j, decaying as ${|i\ensuremath{-}j|}^{\ensuremath{-}\ensuremath{\alpha}}$, with $\ensuremath{\alpha}>1$. For $\ensuremath{\alpha}$ smaller than a certain threshold ${\ensuremath{\alpha}}_{0}$, which is larger than 2 and depends on the ratio
J/g, the ground state consists of an ordered sequence of segments with equal length and alternating magnetization. The width of the segments depends on both $\ensuremath{\alpha}$ and the ratio
J/g. Our Monte Carlo study shows that the on-site magnetization vanishes at finite temperatures and finds no indication of any phase transition. Yet, the modulation present in the ground state is recovered at finite temperatures in the two-point correlation function, which oscillates in space with a characteristic spatial period. The latter depends on $\ensuremath{\alpha}$ and
J/g and decreases smoothly from the ground-state value as the temperature is increased. Such an oscillation of the correlation function is exponentially damped over a characteristic spatial scale, the correlation length, which asymptotically diverges roughly as the inverse of the temperature as
T=0 is approached. This suggests that the long-range interaction causes the Ising chain to fall into a universality class consistent with an underlying continuous symmetry. The ${e}^{\ensuremath{\Delta}/T}$ temperature dependence of the correlation length and the uniform ferromagnetic ground state, characteristic of the
g=0 discrete Ising symmetry, are recovered for $\ensuremath{\alpha}>{\ensuremath{\alpha}}_{0}$.
An importance-sampling Monte Carlo method is used to study $N\ifmmode\times\else\texttimes\fi{}N$ Ising square lattices with nearest-neighbor interactions and either free edges or periodic boundary conditions. The internal energy, specific heat, order parameter, susceptibility, and near-neighbor spin-spin correlation functions of the finite lattices are determined as a function of
N and extrapolated to the corresponding infinite-system values. The effect of finite size is greater for free edges in all cases. The results agree well with predictions of finite size scaling theory and the shape functions as well as amplitudes of surface contribution terms are determined.