Bilinear-biquadratic Heisenberg Research Articles

Spin transport in the bilinear–biquadratic Heisenberg model on square lattice

We analyze the spin coefficients in the spin-1 bilinear–biquadratic Heisenberg model (BBQ) using the SU(3) Schwinger boson approach. The calculations were performed for the model in the ferro- and antiferromagnetic phases, where the AC and DC spin conductivities where computed and the effect of the biquadratic coupling Jbq2 and next-nearest neighboring exchange interactions J2 on spin conductivity analyzed. The aim is to verify the influence of phase transition between the different phases of the system on regular part of the spin conductivity σreg(ω) (AC contribution) and Drude’s weight, that represents the DC contribution. We get that the highest peaks in the spin conductivity in the ferro- and antiferromagnetic phases of the model depends strongly on the values of the exchange coupling constants J2 and Jbq2, as expected.

Phase diagrams and excitations of anisotropic S=1 quantum magnets on the triangular lattice

The $S=1$ bilinear-biquadratic Heisenberg exchange model on the triangular lattice with a single-ion anisotropy has previously been shown to host a number of exotic magnetic and nematic orders [Moreno-Cardoner $\textit{et al.}$, Phys. Rev. B $\textbf{90}$, 144409 (2014)], including an extensive region of "supersolid" order. In this work, we amend the model by an XXZ anisotropy in the exchange interactions. Tuning to the limit of an exactly solvable $S=1$ generalized Ising-/Blume-Capel-type model provides a controlled limit to access phases at finite transverse exchange. Notably, we find an additional macroscopically degenerate region in the phase diagram and study its fate under perturbation theory. We further map out phase diagrams as a function of the XXZ anisotropy parameter, ratio of bilinear and biquadratic interactions and single-ion anisotropy, and compute corrections to the total ordered moment in various phases using systematically constructed linear flavor-wave theory. We also present linear flavor-wave spectra of various states, finding that the lowest-energy band in three-sublattice generalized (i.e. with $S^z=\pm1,0$) Ising/Blume-Capel states, stabilized by strong exchange anisotropies, is remarkably flat, opening up the way to flat-band engineering of magnetic excitation via stabilizing non-trivial Ising-ordered ground states.

Dimensional enhancements in a quantum battery with imperfections

Power storage devices are shown to increase their efficiency if they are designed by using quantum systems. We show that the average power output of a quantum battery based on a quantum interacting spin model, charged via a local magnetic field, can be enhanced with the increase of spin quantum number. In particular, we demonstrate such increment in the power output when the initial state of the battery is prepared as the ground or canonical equilibrium state of the spin-j XY model and the bilinear-biquadratic spin-j Heisenberg chain (BBH) in presence of the transverse magnetic field. Interestingly, we observe that in the case of the XY model, a trade-off relation exists between the range of interactions in which the power increases and the dimension while for the BBH model, the improvements depend on the phase in which the initial state is prepared. Moreover, we exhibit that such dimensional advantages persist even when the battery-Hamiltonian has some defects or when the initial battery-state is prepared at finite temperature.

On a Class of Orthogonal-Invariant Quantum Spin Systems on the Complete Graph

AbstractWe study a two-parameter family of quantum spin systems on the complete graph, which is the most general model invariant under the complex orthogonal group. In spin $S=\frac {1}{2}$ it is equivalent to the XXZ model, and in spin $S=1$ to the bilinear-biquadratic Heisenberg model. The paper is motivated by the work of Björnberg, whose model is invariant under the (larger) complex general linear group. In spin $S=\frac {1}{2}$ and $S=1$ we give an explicit formula for the free energy for all values of the two parameters, and for spin $S>1$ for when one of the parameters is non-negative. This allows us to draw phase diagrams and determine critical temperatures. For spin $S=\frac {1}{2}$ and $S=1$, we give the left and right derivatives as the strength parameter of a certain magnetisation term tends to zero, and we give a formula for a certain total spin observable, and heuristics for the set of extremal Gibbs states in several regions of the phase diagrams, in the style of a recent paper of Björnberg, Fröhlich, and Ueltschi. The key technical tool is expressing the partition function in terms of the irreducible characters of the symmetric group and the Brauer algebra. The parameters considered include, and go beyond, those for which the systems have probabilistic representations as interchange processes.

Bilinear-biquadratic spin-1 model in the Haldane and dimerized phases

We study the low-lying spectrum of the bilinear-biquadratic Heisenberg model in the dimerized and Haldane phases using a tensor renormalization method. At the critical point $\ensuremath{\theta}=\ensuremath{-}\ensuremath{\pi}/4$, the finite-size spectrum predicted by the Wess-Zumino-Witten (WZW) model can only partly be confirmed. We find a singlet-singlet gap which does not fit into the WZW systematics. The results obtained are compared to Bethe ansatz, exact diagonalization, and density-matrix renormalization group calculations for specific parameters.

Universal Qudit Hamiltonians

A family of quantum Hamiltonians is said to be universal if any other finite-dimensional Hamiltonian can be approximately encoded within the low-energy space of a Hamiltonian from that family. If the encoding is efficient, universal families of Hamiltonians can be used as universal analogue quantum simulators and universal quantum computers, and the problem of approximately determining the ground-state energy of a Hamiltonian from a universal family is QMA-complete. One natural way to categorise Hamiltonians into families is in terms of the interactions they are built from. Here we prove universality of some important classes of interactions on qudits (d-level systems):We completely characterise the k-qudit interactions which are universal, if augmented with arbitrary Hermitian 1-local terms.We find that, for all k geqslant 2 and all local dimensions d geqslant 2, almost all such interactions are universal aside from a simple stoquastic class.We prove universality of generalisations of the Heisenberg model that are ubiquitous in condensed-matter physics, even if free 1-local terms are not provided. We show that the SU(d) and SU(2) Heisenberg interactions are universal for all local dimensions d geqslant 2 (spin geqslant 1/2), implying that a quantum variant of the Max-d-Cut problem is QMA-complete. We also show that for d=3 all bilinear-biquadratic Heisenberg interactions are universal. One example is the general AKLT model.We prove universality of any interaction proportional to the projector onto a pure entangled state.

Towards a quantum Monte Carlo for lattice systems

In this work we build the foundations of a quantum Monte Carlo as a stochastic numerical method to solve lattice many-body quantum systems with nearest-neighbor interactions at most. As motivation, we briefly describe the bilinear-biquadratic Heisenberg model with an external field, for spin-1 particles, as an effective Hamiltonian of the Bose-Hubbard model with an external quadratic Zeeman field in the Mott insulator phase at unit filling. Then, we discuss how to implement the world line Monte Carlo with local updates to circumvent the difficulties that arise on these type of systems by mapping the quantum partition function into the one of an effective classical model, in one additional dimension, given by the imaginary time evolution of the system. Such a mapping is performed by means of the Suzuki-Trotter decomposition, which transforms the original partition function into a summation of weights given by the classical configurations. Later, we present a set of observables that can be measured through this method and show how to use a Metropolis update scheme to accomplish the measurements. At last, we present the maximization of the configuration weights for three parameter sets as the first and relevant step to perform future measurements.

Anisotropic magnetic excitations of a frustrated bilinear-biquadratic spin model: Implications for spin waves of detwinned iron pnictides

Elucidating the nature of spin excitations is important to understanding the mechanism of superconductivity in the iron pnictides. Motivated by recent inelastic neutron scattering measurements in the nearly 100% detwinned BaFe$_{2}$As$_{2}$, we study the spin dynamics of an $S=1$ frustrated bilinear-biquadratic Heisenberg model in the antiferromagnetic phase with wavevector $(\pi,0)$. The biquadratic interactions are treated in a dynamical way using a flavor-wave theory in an $SU(3)$ representation. Besides the magnon (dipolar) excitations, the biquadratic interactions give rise to quadrupolar excitations at high energies. We find that the quadrupolar wave significantly influences, in an energy dependent way, the anisotropy between the spin excitation spectra along the $(\pi,0)$ and $(0,\pi)$ directions in the wave vector space. Our theoretical results capture the essential behavior of the spin dynamics measured in the antiferromagnetic phase of the detwinned BaFe$_2$As$_2$. More generally, our results underscore the importance of electron correlation effects for the microscopic physics of the iron pnictides.

Density matrix renormalization group study of nematicity in two dimensions: Application to a spin-1 bilinear-biquadratic model on the square lattice

Nematic order is an exotic property observed in several strongly correlated systems, such as the iron-based superconductors. Using large-scale density matrix renormalization group (DMRG) techniques, we study at zero-temperature the nematic spin liquid that competes with spin dipolar and quadrupolar orders. We use these nematic orders to characterize different quantum phases and quantum phase transitions. More specifically, we study a spin-$1$ bilinear-biquadratic Heisenberg model on the square lattice with couplings beyond nearest neighbors. We focus on parameter regions around the highly symmetric $SU(3)$ point where the bilinear and biquadratic interactions are equal. With growing further-neighbor biquadratic interactions, we identify different spin dipolar and quadrupolar orders. We find that the DMRG results on cylindrical geometries correctly detect nematicity in different quantum states and accurately characterize the quantum phase transitions among them. Therefore, spin-driven nematicity -- here defined as the spontaneous breaking of the lattice invariance under a 90$^o$ rotation -- is an order parameter which can be studied directly in DMRG calculations in two dimensions in different quantum states.

Broken mirror symmetry, incommensurate spin correlations, and B2g nematic order in iron pnictides

Motivated by growing indications for a distinct form of nematic correlations in the extremely hole doped iron pnictide compounds $A$Fe$_2$As$_2$ ($A$=K,Rb,Cs), we consider spin-driven nematic order in the general case of incommensurate magnetic fluctuations. We classify the nematic order parameters by broken mirror symmetries of the tetragonal $D_{4h}$ point group, and use this scheme to construct a general Ginzburg-Landau theory that links the nematic order to spatial pattern of magnetic fluctuations. Our analysis points to incommensurate magnetic fluctuations of wavevector $(q,q)$ as underlying a $B_{2g}$ nematic order in $A$Fe$_2$As$_2$. We substantiate this idea by microscopic calculations of the nematic order based on 3-sublattice $(2\pi/3,2\pi/3)$ spin correlations in an extended bilinear-biquadratic Heisenberg model. We summarize the existing evidence in support of the proposed mechanism, suggest ways to further test it in future experiments, and discuss the implications of our results for iron-based high temperature superconductivity.

Universal tensor-network algorithm for any infinite lattice

We present a general graph-based Projected Entangled-Pair State (gPEPS) algorithm to approximate ground states of nearest-neighbor local Hamiltonians on any lattice or graph of infinite size. By introducing the structural-matrix which codifies the details of tensor networks on any graphs in any dimension $d$, we are able to produce a code that can be essentially launched to simulate any lattice. We further introduce an optimized algorithm to compute simple tensor updates as well as expectation values and correlators with a mean-field-like effective environments. Though not being variational, this strategy allows to cope with PEPS of very large bond dimension (e.g., $D=100$), and produces remarkably accurate results in the thermodynamic limit in many situations, and specially when the correlation length is small and the connectivity of the lattice is large. We prove the validity of our approach by benchmarking the algorithm against known results for several models, i.e., the antiferromagnetic Heisenberg model on a chain, star and cubic lattices, the hardcore Bose-Hubbard model on square lattice, the ferromagnetic Heisenberg model in a field on the pyrochlore lattice, as well as the $3$-state quantum Potts model in field on the kagome lattice and the spin-$1$ bilinear-biquadratic Heisenberg model on the triangular lattice. We further demonstrate the performance of gPEPS by studying the quantum phase transition of the $2d$ quantum Ising model in transverse magnetic field on the square lattice, and the phase diagram of the Kitaev-Heisenberg model on the hyperhoneycomb lattice. Our results are in excellent agreement with previous studies.

Ground-state study of the spin-1 bilinear-biquadratic Heisenberg model on the triangular lattice using tensor networks

Making use of infinite projected entangled pair states, we investigate the ground state phase diagram of the nearest-neighbor spin-1 bilinear-biquadratic Heisenberg model on the triangular lattice. In agreement with previous studies, we find the ferromagnetic, 120 degree magnetically ordered, ferroquadrupolar and antiferroquadrupolar phases, and confirm that all corresponding phase transitions are first order. Moreover, we provide an accurate estimate of the location of the ferroquadrupolar to 120 degree magnetically ordered phase transition, thereby fully establishing the phase diagram. Also, we do not encounter any signs of the existence of a quantum paramagnetic phase. In particular, contrary to the equivalent square lattice model, we demonstrate that on the triangular lattice the one-dimensional Haldane phase does not reach all the way up to the two-dimensional limit. Finally, we investigate the possibility of an intermediate partially-magnetic partially-quadrupolar phase close to $\theta = \pi/2$, and we show that, also contrary to the square lattice case, this phase is not present on the triangular lattice.

Dimerization in the anisotropic bilinear-biquadratic Heisenberg model

Two recent publications report different boundaries for the dimerized phase of the bilinear-biquadratic spin-1 Heisenberg model with quadratic Zeeman anisotropy (see Refs. [9] and [10]). We numerically address these discrepancies for the purely biquadratic model with quadratic Zeeman term using methods based on matrix product states (MPS). We extrapolate the results beyond the purely biquadratic model in order to determine the phase boundaries of the dimerized phase and propose a modified phase diagram. Our MPS implementation incorporates SU(2) symmetry explicitly and we discuss its implementation into the MPS representation in the presence of the Zeeman anisotropy.

A tensor network study of the complete ground state phase diagram of the spin-1 bilinear-biquadratic Heisenberg model on the square lattice

Using infinite projected entangled pair states, we study the ground state phase diagram of the spin-1 bilinear-biquadratic Heisenberg model on the square lattice directly in the thermodynamic limit. We find an unexpected partially nematic partially magnetic phase in between the antiferroquadrupolar and ferromagnetic regions. Furthermore, we describe all observed phases and discuss the nature of the phase transitions involved.

Dimerization in ultracold spinor gases with Zeeman splitting

Two recent publications report different boundaries for the dimerized phase of the bilinear-biquadratic spin-1 Heisenberg model with quadratic Zeeman effect. We address these discrepancies for the biquadratic model with quadratic Zeeman term and explain the differences. Based on our numerical results the phase boundaries of the dimerized phase are determined.

Bilinear-biquadratic spin-1 rings: an SU(2)-symmetric MPS algorithm for periodic boundary conditions

An efficient algorithm for SU(2) symmetric matrix product states with periodic boundary conditions is proposed and implemented. It is applied to a study of the spectrum and correlation properties of the spin-1 bilinear-biquadratic Heisenberg model. We characterize the various phases of this model by the lowest states of the spectrum with angular momentum for systems of up to 100 spins. Furthermore, we provide precision results for the dimerization correlator as well as the string correlator.

Emergent Haldane phase in the S=1 bilinear-biquadratic Heisenberg model on the square lattice

Infinite projected entangled pair states simulations of the $S=1$ bilinear-biquadratic Heisenberg model on the square lattice reveal an emergent Haldane phase in between the previously predicted antiferromagnetic and 3-sublattice 120$^\circ$ magnetically ordered phases. This intermediate phase preserves SU(2) spin and translational symmetry but breaks lattice rotational symmetry, and it can be adiabatically connected to the Haldane phase of decoupled $S=1$ chains. Our results contradict previous studies which found a direct transition between the two magnetically ordered states.

Color ice states, weathervane modes, and order by disorder in the bilinear-biquadratic pyrochlore Heisenberg antiferromagnet

We study the pyrochlore Heisenberg antiferromagnet with additional positive biquadratic interaction in the semiclassical limit. The classical ground state manifold of the model contains an extensively large family of non-coplanar spin states known as "color ice states". Starting from a color ice state, a subset of spins may rotate collectively at no energy cost. Such excitation may be viewed in this three-dimensional system as a "membrane-like" analog of the well-known weathervane modes in the classical kagome Heisenberg antiferromagnet. We investigate the weathervane modes in detail and elucidate their physical properties. Furthermore, we study the order by disorder phenomenon in this model, focusing on the role of harmonic fluctuations. Our computationally-limited phase space search suggests that quantum fluctuations select three different states as the magnitude of the biquadratic interaction increases relative to the bilinear interaction, implying a sequence of phase transitions solely driven by fluctuations.

Spin Ferroquadrupolar Order in the Nematic Phase of FeSe.

We provide evidence that spin ferroquadrupolar (FQ) order is the likely ground state in the nonmagnetic nematic phase of stoichiometric FeSe. By studying the variational mean-field phase diagram of a bilinear-biquadratic Heisenberg model up to the 2nd nearest neighbor, we find the FQ phase in close proximity to the columnar antiferromagnet commonly realized in iron-based superconductors; the stability of the FQ phase is further verified by the density matrix renormalization group. The dynamical spin structure factor in the FQ state is calculated with flavor-wave theory, which yields a qualitatively consistent result with inelastic neutron scattering experiments on FeSe at both low and high energies. We verify that FQ can coexist with C_{4} breaking environments in the mean-field calculation, and further discuss the possibility that quantum fluctuations in FQ act as a source of nematicity.

Existence of Néel Order in the S=1 Bilinear-Biquadratic Heisenberg Model via Random Loops

We consider the general spin-1 SU(2) invariant Heisenberg model with a two-body interaction. A random loop model is introduced and relations to quantum spin systems is proved. Using this relation it is shown that for dimensions 3 and above N\'eel order occurs for a large range of values of the relative strength of the bilinear ($-J_1$) and biquadratic ($-J_2$) interaction terms. The proof uses the method of reflection positivity and infrared bounds. Links between spin correlations and loop correlations are proved.