Dimer-dimer Correlation Functions Research Articles

Pinwheel valence bond crystal ground state of the spin- 12 Heisenberg antiferromagnet on the shuriken lattice

We investigate the nature of the ground-state of the spin-$\frac{1}{2}$ Heisenberg antiferromagnet on the $shuriken$ lattice by complementary state-of-the-art numerical techniques, such as variational Monte Carlo (VMC) with versatile Gutzwiller-projected Jastrow wave functions, unconstrained multi-variable variational Monte Carlo (mVMC), and pseudo-fermion/Majorana functional renormalization group (PF/PM-FRG) methods. We establish the presence of a quantum paramagnetic ground state and investigate its nature, by classifying symmetric and chiral quantum spin liquids, and inspecting their instabilities towards competing valence-bond-crystal (VBC) orders. Our VMC analysis reveals that a VBC with a pinwheel structure emerges as the lowest-energy variational ground state, and it is obtained as an instability of the U(1) Dirac spin liquid. Analogous conclusions are drawn from mVMC calculations employing accurate BCS pairing states supplemented by symmetry projectors, which confirm the presence of pinwheel VBC order by a thorough analysis of dimer-dimer correlation functions. Our work highlights the nontrivial role of quantum fluctuations via the Gutzwiller projector in resolving the subtle interplay between competing orders.

Topological sectors, dimer correlations, and monomers from the transfer-matrix solution of the dimer model

We solve the classical square-lattice dimer model with periodic boundaries and in the presence of a field t that couples to the (vector) flux, by diagonalizing a modified version of Lieb's transfer matrix. After deriving the torus partition function in the thermodynamic limit, we show how the configuration space divides into topological sectors corresponding to distinct values of the flux. Additionally, we demonstrate in general that expectation values are t independent at leading order, and obtain explicit expressions for dimer occupation numbers, dimer-dimer correlation functions, and the monomer distribution function. The last of these is expressed as a Toeplitz determinant, whose asymptotic behavior for large monomer separation is tractable using the Fisher-Hartwig conjecture. Our results reproduce those previously obtained using Pfaffian techniques.

Mean-field theory of interacting triplons in a two-dimensional valence-bond solid: Stability and properties of many-triplon states

We study a system of interacting triplons (the elementary excitations of a valence-bond solid) described by an effective interacting boson model derived within the bond-operator formalism. In particular, we consider the square lattice spin-1/2 $J_1$-$J_2$ antiferromagnetic Heisenberg model, focus on the intermediate parameter region, where a quantum paramagnetic phase sets in, and consider the columnar valence-bond solid phase. Within the bond-operator theory, the Heisenberg model is mapped into an effective boson model in terms of triplet operators $t$. The effective boson model is studied at the harmonic approximation and the energy of the triplons and the expansion of the triplon operators $b$ in terms of the triplet operators $t$ are determined. Such an expansion allows us to performed a second mapping, and therefore, determine an effective interacting boson model in terms of the triplon operators $b$. We then consider systems with a fixed number of triplons and determined the ground-state energy and the spectrum of elementary excitations within a mean-field approximation. We show that many-triplon states are stable, the lowest-energy ones are constituted by a small number of triplons, and the excitation gaps are finite. For $J_2=0.48 J_1$ and $J_2=0.52 J_1$, we also calculate spin-spin and dimer-dimer correlation functions, dimer order parameters, and the bipartite von Neumann entanglement entropy within our mean-field formalism in order to determine the properties of the many-triplon state as a function of the triplon number. We find that the spin and the dimer correlations decay exponentially and that the entanglement entropy obeys an area law, regardless the triplon number. Moreover, only for $J_2=0.48 J_1$, the spin correlations indicate that the many-triplon states with large triplon number might display a more homogeneous singlet pattern than the columnar valence-bond solid.

Quantum criticality on a chiral ladder: An SU(2) infinite density matrix renormalization group study

In this paper we study the ground state properties of a ladder Hamiltonian with chiral $SU(2)$-invariant spin interactions, a possible first step towards the construction of truly two dimensional non-trivial systems with chiral properties starting from quasi-one dimensional ones. Our analysis uses a recent implementation by us of $SU(2)$ symmetry in tensor network algorithms, specifically for infinite Density Matrix Renormalization Group (iDMRG). After a preliminary analysis with Kadanoff coarse-graining and exact diagonalization for a small-size system, we discuss its bosonization and recap the continuum limit of the model to show that it corresponds to a conformal field theory, in agreement with our numerical findings. In particular, the scaling of the entanglement entropy as well as finite-entanglement scaling data show that the ground state properties match those of the universality class of a $c = 1$ conformal field theory (CFT) in $(1+1)$ dimensions. We also study the algebraic decay of spin-spin and dimer-dimer correlation functions, as well as the algebraic convergence of the ground state energy with the bond dimension, and the entanglement spectrum of half an infinite chain. Our results for the entanglement spectrum are remarkably similar to those of the spin-$1/2$ Heisenberg chain, which we take as a strong indication that both systems are described by the same CFT at low energies, i.e., an $SU(2)_1$ Wess-Zumino-Witten theory. Moreover, we explain in detail how to construct Matrix Product Operators for $SU(2)$-invariant three-spin interactions, something that had not been addressed with sufficient depth in the literature.

Non-Abelian chiral spin liquid in a quantum antiferromagnet revealed by an iPEPS study

Abelian and non-Abelian topological phases exhibiting protected chiral edge modes are ubiquitous in the realm of the Fractional Quantum Hall (FQH) effect. Here, we investigate a spin-1 Hamiltonian on the square lattice which could, potentially, host the spin liquid analog of the (bosonic) non-Abelian Moore-Read FQH state, as suggested by Exact Diagonalisation of small clusters. Using families of fully SU(2)-spin symmetric and translationally invariant chiral Projected Entangled Pair States (PEPS), variational energy optimization is performed using infinite-PEPS methods, providing good agreement with Density Matrix Renormalisation Group (DMRG) results. A careful analysis of the bulk spin-spin and dimer-dimer correlation functions in the optimized spin liquid suggests that they exhibit long-range "gossamer tails". We argue these tails are finite-$D$ artifacts of the chiral PEPS, which become irrelevant when the PEPS bond dimension $D$ is increased. From the investigation of the entanglement spectrum, we observe sharply defined chiral edge modes following the prediction of the SU(2)$_2$ Wess-Zumino-Witten theory and exhibiting a conformal field theory (CFT) central charge $c=3/2$, as expected for a Moore-Read chiral spin liquid. We conclude that the PEPS formalism offers an unbiased and efficient method to investigate non-Abelian chiral spin liquids in quantum antiferromagnets.

Height fluctuations in interacting dimers

We consider a non-integrable model for interacting dimers on the two-dimensional square lattice. Configurations are perfect matchings of $\mathbb Z^2$, i.e. subsets of edges such that each vertex is covered exactly once ("close-packing" condition). Dimer configurations are in bijection with discrete height functions, defined on faces $\boldsymbol{\xi}$ of $\mathbb Z^2$. The non-interacting model is "integrable" and solvable via Kasteleyn theory; it is known that all the moments of the height difference $h_{\boldsymbol{\xi}}-h_{\boldsymbol{\eta}}$ converge to those of the massless Gaussian Free Field (GFF), asymptotically as $|{\boldsymbol{\xi}}-{\boldsymbol{\eta}}|\to \infty$. We prove that the same holds for small non-zero interactions, as was conjectured in the theoretical physics literature. Remarkably, dimer-dimer correlation functions are instead not universal and decay with a critical exponent that depends on the interaction strength. Our proof is based on an exact representation of the model in terms of lattice interacting fermions, which are studied by constructive field theory methods. In the fermionic language, the height difference $h_{\boldsymbol{\xi}}-h_{\boldsymbol{\eta}}$ takes the form of a non-local operator, consisting of a sum of monomials along an {\it arbitrary} path connecting $\boldsymbol{\xi}$ and $\boldsymbol{\eta}$. As in the non-interacting case, this path-independence plays a crucial role in the proof.

Quantum Spin Liquid in Spin 1/2J1–J2Heisenberg Model on Square Lattice: Many-Variable Variational Monte Carlo Study Combined with Quantum-Number Projections

The nature of quantum spin liquids is studied for the spin-$1/2$ antiferromagnetic Heisenberg model on a square lattice containing exchange interactions between nearest-neighbor sites, $J_1$, and those between next-nearest-neighbor sites, $J_2$. We perform variational Monte Carlo simulations together with the quantum-number-projection technique and clarify the phase diagram in the ground state together with its excitation spectra. We obtain the nonmagnetic phase in the region $0.4< J_2/J_1\le 0.6$ sandwiched by the staggered and stripe antiferromagnetic (AF) phases. Our direct calculations of the spin gap support the notion that the triplet excitation from the singlet ground state is gapless in the range of $0.4 < J_2/J_1 \le 0.5$, while the gapped valence-bond-crystal (VBC) phase is stabilized for $0.5 < J_2/J_1 \le 0.6$. The VBC order is likely to have the columnar symmetry with a spontaneous symmetry breaking of the $C_{4v}$ symmetry. The power-law behaviors of the spin-spin and dimer-dimer correlation functions in the gapless region are consistent with the emergence of the algebraic quantum-spin-liquid phase (critical phase). The exponent of the spin correlation $\langle S(0)S(r)\rangle \propto 1/r^{z+\eta}$ at a long distance $r$ appears to increase from $z+\eta\sim 1$ at $J_2/J_1\sim0.4$ toward the continuous transition to the VBC phase at $J_1/J_1\sim0.5$. Our results, however, do not fully exclude the possibility of a direct quantum transition between the staggered AF and VBC phases with a wide critical region and deconfined criticality.

Featureless quantum spin liquid,13-magnetization plateau state, and exotic thermodynamic properties of the spin-12frustrated Heisenberg antiferromagnet on an infinite Husimi lattice

By utilizing tensor-network-based methods, we investigate the zero- and finite-temperature properties of the spin-1/2 Heisenberg antiferromagnetic (HAF) model on an infinite Husimi lattice that contains 3/2 sites per triangle. The ground state of this model is found to possess vanishing local magnetization and is featureless; the spin-spin and dimer-dimer correlation functions are verified to decay exponentially; and its ground-state energy per site is determined to be $e_0=-0.4343(1)$, which is very close to that [$e_0=-0.4386(5)$] of the intriguing kagome HAF model. The magnetization curve shows the absence of a zero-magnetization plateau, implying a gapless excitation. A $1/3$-magnetization plateau with spin up-up-down state is observed, which is selected and stabilized by quantum fluctuations. A ground state phase diagram under magnetic fields is presented. Moreover, both magnetic susceptibility and the specific heat are studied, whose low-temperature behaviors reinforce the conclusion that the HAF model on the infinite Husimi lattice owns a gapless and featureless spin liquid ground state.

Plaquette valence-bond crystal in the frustrated Heisenberg quantum antiferromagnet on the square lattice

Using both exact diagonalizations and diagonalizations in a subset of short-range valence bond singlets, we address the nature of the groundstate of the Heisenberg spin-1/2 antiferromagnet on the square lattice with competing next-nearest and next-next-nearest neighbor antiferromagnetic couplings (J1-J2-J3 model). A detailed comparison of the two approaches reveals a region along the line (J2+J3)/J1=1/2, where the description in terms of nearest-neighbor singlet coverings is excellent, therefore providing evidence for a magnetically disordered region. Furthermore a careful analysis of dimer-dimer correlation functions, dimer structure factors and plaquette-plaquette correlation functions provides striking evidence for the presence of a plaquette valence bond solid order in part of the magnetically disordered region.

Zero-temperature properties of the quantum dimer model on the triangular lattice

Using exact diagonalizations and Green's function Monte Carlo simulations, we have studied the zero-temperature properties of the quantum dimer model on the triangular lattice on clusters with up to 588 sites. A detailed comparison of the properties in different topological sectors as a function of the cluster size and for different cluster shapes has allowed us to identify different phases, to show explicitly the presence of topological degeneracy in a phase close to the Rokhsar-Kivelson point, and to understand finite-size effects inside this phase. The nature of the various phases has been further investigated by calculating dimer-dimer correlation functions. The present results confirm and complement the phase diagram proposed by Moessner and Sondhi on the basis of finite-temperature simulations [Phys. Rev. Lett. {\bf 86}, 1881 (2001)].

Dynamical dimer-dimer correlation functions from exact diagonalization

A regularization method is presented to deduce dynamic correlation functions from exact diagonalization calculations. It is applied to dimer-dimer correlation functions in quantum spin chains relevant for the description of spin-Peierls systems. Exact results for the $\mathrm{XY}$ model are presented. The analysis draws into doubt that the dimer-dimer correlation functions show the same scale invariance as spin-spin correlation functions. The results are applied to describe the quasielastic scattering in ${\mathrm{CuGeO}}_{3}$ and the hardening of the Peierls-active phonons.

Dirac structure, RVB, and Goldstone modes in thekagoméantiferromagnet

We show that there exists a long-range RVB state for the kagome lattice spin-1/2 Heisenberg antiferromagnet for which the spinons have a massless Dirac spectrum. By considering various perturbations of the RVB state which give mass to the fermions by breaking a symmetry, we are able to describe a wide-ranging class of known states on the kagome lattice, including spin-Peierls solid and chiral spin liquid states. Using an RG treatment of fluctuations about the RVB state, we propose yet a different symmetry breaking pattern and show how collective excitations about this state account for the gapless singlet modes seen experimentally and numerically. We make further comparison with numerics for Chern numbers, dimer-dimer correlation functions, the triplet gap, and other quantities. To accomplish these calculations, we propose a variant of the SU(N) theory which enables us to include many of the effects of Gutzwiller projection at the mean-field level.

Test of finite-size scaling predictions in three dimensions

A constrained monomer-dimer (CMD) model on decorated hypercubic lattices, which is exactly solvable in any dimensionality d, is studied. The model is critical in the limit of maximal packing of the dimer configurations allowed by the constraint. For d=2 it is equivalent to the rooted-tree model on the square lattice considered by Duplantier and David (1988). The validity of two finite-size scaling predictions in higher dimensionalities is studied: (1) the existence of logarithmic corrections in the free energy due to corners, and (ii) the amplitude-exponent relation for the pair correlation function. The logarithmic finite-size corrections are obtained for arbitrary d, for boundary conditions which are periodic for d'>or=0 dimensions and free in the remaining d-d' dimensions. Dimer-dimer correlation functions are studied at d=3 for a system finite in two dimensions and infinite in the third. It is shown that the amplitude-exponent relation for the CMD model holds for d=2 and breaks down for d=3.