arXiv:1512.03030 · Physical review. B./Physical review. B

Quantum phase diagram of a frustrated antiferromagnet on the bilayer honeycomb lattice

Hao Zhang, Carlos Alberto Lamas, Marcelo Arlego, Wolfram Brenig

Published Jun 27, 2016• 19 citations

Abstract

We study the spin-1/2 Heisenberg antiferromagnet on a bilayer honeycomb lattice including interlayer frustration. Using a set of complementary approaches, namely, Schwinger bosons, dimer series expansion, bond operators, and exact diagonalization, we map out the quantum phase diagram. Analyzing ground-state energies and elementary excitation spectra, we find four distinct phases, corresponding to three collinear magnetic long-range ordered states, and one quantum disordered interlayer dimer phase. We detail that the latter phase is adiabatically connected to an exact singlet product ground state of the bilayer, which exists along a line of maximum interlayer frustration. The order within the remaining three phases will be clarified.

DOI: 10.1103/physrevb.93.235150

Related papers

Quantum phases in the frustrated Heisenberg model on the bilayer honeycomb lattice

Hao Zhang, Marcelo Arlego, Carlos Alberto Lamas

arXiv:1312.6324Physical Review BJan 6, 201431 citations

Self consistent study of the quantum phases in a frustrated antiferromagnet on the bilayer honeycomb lattice

Arlego, Marcelo José Fabián, Carlos Alberto Lamas, Zhang, H.

Apr 1, 20149 citations

Quantum phase diagram of a frustrated spin-$\frac{1}{2}$ $J_{1}$--$J_{2}$--$J_{1}^{\perp}$ Heisenberg model on an $AA$-stacked bilayer honeycomb lattice

P. H. Y. Li, R. F. Bishop

May 3, 20180 citations

Quantum phase diagram of the bilayer Kitaev-Heisenberg model

Elahe Samimi, Antonia Duft, Patrick Adelhardt +2 more

arXiv:2412.17495Physical review. B./Physical review. BApr 16, 20260 citations

Series expansion studies of the J 1 − J 2 Heisenberg bilayer

Erik Wagner, Wolfram Brenig

Physical review. B./Physical review. BSep 14, 20230 citations

Series expansion studies of the $J_1$-$J_2$-Heisenberg bilayer

Wagner, Erik, Wolfram Brenig

arXiv:2302.11599Feb 22, 20230 citations

arXiv:1512.03030 · Physical review. B./Physical review. B

Quantum phase diagram of a frustrated antiferromagnet on the bilayer honeycomb lattice

Hao Zhang, Carlos Alberto Lamas, Marcelo Arlego, Wolfram Brenig

Published Jun 27, 2016• 19 citations

arXiv HTML PDF

Abstract

DOI: 10.1103/physrevb.93.235150

Related papers

Quantum phases in the frustrated Heisenberg model on the bilayer honeycomb lattice

Hao Zhang, Marcelo Arlego, Carlos Alberto Lamas

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 ${J}_{2}$ 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.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 ${J}_{2}$ 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.

arXiv:1312.6324Physical Review BJan 6, 201431 citations

Self consistent study of the quantum phases in a frustrated antiferromagnet on the bilayer honeycomb lattice

Arlego, Marcelo José Fabián, Carlos Alberto Lamas, Zhang, H.

We study the frustrated Heisenberg model on the bilayer honeycomb lattice. The ground-state energy and spin gap are calculated, using different bosonic representations at mean field level and numerical calculations, to explore different sectors of the phase diagram. In particular we make use of a bond operator formalism and series expansion calculations to study the extent of dimer inter-layer phase. On the other hand we use the Schwinger boson method and exact diagonalization on small systems to analyze the evolution of on-layer phases. In this case we specifically observe a phase that presents a spin gap and short range Neel correlations that survives even in the presence of non-zero next-nearest-neighbor interaction and inter-layer coupling.We study the frustrated Heisenberg model on the bilayer honeycomb lattice. The ground-state energy and spin gap are calculated, using different bosonic representations at mean field level and numerical calculations, to explore different sectors of the phase diagram. In particular we make use of a bond operator formalism and series expansion calculations to study the extent of dimer inter-layer phase. On the other hand we use the Schwinger boson method and exact diagonalization on small systems to analyze the evolution of on-layer phases. In this case we specifically observe a phase that presents a spin gap and short range Neel correlations that survives even in the presence of non-zero next-nearest-neighbor interaction and inter-layer coupling.

Apr 1, 20149 citations

Quantum phase diagram of a frustrated spin- $\frac{1}{2}$ $J_{1}$ -- $J_{2}$ -- $J_{1}^{\perp}$ Heisenberg model on an $AA$ -stacked bilayer honeycomb latticeQuantum phase diagram of a frustrated spin-$\frac{1}{2}$ $J_{1}$--$J_{2}$--$J_{1}^{\perp}$ Heisenberg model on an $AA$-stacked bilayer honeycomb lattice

P. H. Y. Li, R. F. Bishop

The zero-temperature quantum phase diagram of the spin- $\frac{1}{2}$ $J_{1}$ -- $J_{2}$ -- $J_{1}^{\perp}$ model on an $AA$ -stacked bilayer honeycomb lattice is investigated using the coupled cluster method (CCM). The model comprises two monolayers in each of which the spins, residing on honeycomb-lattice sites, interact via both nearest-neighbor (NN) and frustrating next-nearest-neighbor isotropic antiferromagnetic (AFM) Heisenberg exchange iteractions, with respective strengths $J_{1} > 0$ and $J_{2} \equiv \kappa J_{1}>0$ . The two layers are coupled via a comparable Heisenberg exchange interaction between NN interlayer pairs, with a strength $J_{1}^{\perp} \equiv \delta J_{1}$ . The complete phase boundaries of two quasiclassical collinear AFM phases, namely the Neel and Neel-II phases, are calculated in the $\kappa \delta$ half-plane with $\kappa > 0$ . Whereas on each monolayer in the Neel state all NN pairs of spins are antiparallel, in the Neel-II state NN pairs of spins on zigzag chains along one of the three equivalent honeycomb-lattice directions are antiparallel, while NN interchain spins are parallel. We calculate directly in the thermodynamic (infinite-lattice) limit both the magnetic order parameter $M$ and the excitation energy $\Delta$ from the $s^{z}_{T}=0$ ground state to the lowest-lying $|s^{z}_{T}|=1$ excited state (where $s^{z}_{T}$ is the total $z$ component of spin for the system as a whole, and where the collinear ordering lies along the $z$ direction) for both quasiclassical states used (separately) as the CCM model state, on top of which the multispin quantum correlations are then calculated to high orders ( $n \leq 10$ ) in a systematic series of approximations involving $n$ -spin clusters. The sole approximation made is then to extrapolate the sequences of $n$ th-order results for $M$ and $\Delta$ to the exact limit, $n \to \infty$ .The zero-temperature quantum phase diagram of the spin-$\frac{1}{2}$ $J_{1}$--$J_{2}$--$J_{1}^{\perp}$ model on an $AA$-stacked bilayer honeycomb lattice is investigated using the coupled cluster method (CCM). The model comprises two monolayers in each of which the spins, residing on honeycomb-lattice sites, interact via both nearest-neighbor (NN) and frustrating next-nearest-neighbor isotropic antiferromagnetic (AFM) Heisenberg exchange iteractions, with respective strengths $J_{1} > 0$ and $J_{2} \equiv \kappa J_{1}>0$. The two layers are coupled via a comparable Heisenberg exchange interaction between NN interlayer pairs, with a strength $J_{1}^{\perp} \equiv \delta J_{1}$. The complete phase boundaries of two quasiclassical collinear AFM phases, namely the Neel and Neel-II phases, are calculated in the $\kappa \delta$ half-plane with $\kappa > 0$. Whereas on each monolayer in the Neel state all NN pairs of spins are antiparallel, in the Neel-II state NN pairs of spins on zigzag chains along one of the three equivalent honeycomb-lattice directions are antiparallel, while NN interchain spins are parallel. We calculate directly in the thermodynamic (infinite-lattice) limit both the magnetic order parameter $M$ and the excitation energy $\Delta$ from the $s^{z}_{T}=0$ ground state to the lowest-lying $|s^{z}_{T}|=1$ excited state (where $s^{z}_{T}$ is the total $z$ component of spin for the system as a whole, and where the collinear ordering lies along the $z$ direction) for both quasiclassical states used (separately) as the CCM model state, on top of which the multispin quantum correlations are then calculated to high orders ($n \leq 10$) in a systematic series of approximations involving $n$-spin clusters. The sole approximation made is then to extrapolate the sequences of $n$th-order results for $M$ and $\Delta$ to the exact limit, $n \to \infty$.

May 3, 20180 citations

Quantum phase diagram of the bilayer Kitaev-Heisenberg model

Elahe Samimi, Antonia Duft, Patrick Adelhardt +2 more

We study the ground-state phase diagram of the spin- $1/2$ Kitaev-Heisenberg model on the bilayer honeycomb lattice with large-scale tensor network calculations based on the infinite projected entangled pair state technique as well as high-order series expansions. We find that beyond various magnetically ordered phases, including ferromagnetic, zigzag, antiferromagnetic (AFM) and stripy states, two extended quantum spin liquid phases arise in the proximity of the Kitaev limit. While these ordered phases also appear in the monolayer Kitaev-Heisenberg model, our results further show that a valence bond solid state emerges in a relatively narrow range of parameter space between the AFM and stripy phases, which can be adiabatically connected to isolated Heisenberg dimers. Our results highlight the importance of considering interlayer interactions on the emergence of novel quantum phases in the bilayer Kitaev materials.We study the ground-state phase diagram of the spin-$1/2$ Kitaev-Heisenberg model on the bilayer honeycomb lattice with large-scale tensor network calculations based on the infinite projected entangled pair state technique as well as high-order series expansions. We find that beyond various magnetically ordered phases, including ferromagnetic, zigzag, antiferromagnetic (AFM) and stripy states, two extended quantum spin liquid phases arise in the proximity of the Kitaev limit. While these ordered phases also appear in the monolayer Kitaev-Heisenberg model, our results further show that a valence bond solid state emerges in a relatively narrow range of parameter space between the AFM and stripy phases, which can be adiabatically connected to isolated Heisenberg dimers. Our results highlight the importance of considering interlayer interactions on the emergence of novel quantum phases in the bilayer Kitaev materials.

arXiv:2412.17495Physical review. B./Physical review. BApr 16, 20260 citations

Series expansion studies of the J 1 − J 2 Heisenberg bilayer

Erik Wagner, Wolfram Brenig

We study a bilayer of the frustrated ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ Heisenberg model on the square lattice. Starting from the dimer limit at strong interlayer coupling, we perform series expansions using the perturbative continuous unitary transformation, based on the flow equation method, in order to determine the spectrum up to the two-triplon sector. From the one-triplon dispersion, we obtain quantum critical lines for transitions from the dimer phase into either N\'eel or collinear magnetic order. For low to intermediate frustration these transitions are consistent with existing findings, based on the magnetic phases. In the region of strongest frustration, i.e., ${J}_{2}/{J}_{1}\ensuremath{\sim}0.5$, we provide an estimate for the stability of the anticipated single-layer quantum spin-liquids against finite interlayer coupling. In the two-triplon sector, we find a set of well defined (anti)bound states, which can be classified according to total spin and in-plane rotational symmetry. For vanishing frustration, these states agree with the previous series expansion analysis. For ${J}_{2}/{J}_{1}\ensuremath{\gtrsim}0.5$, we provide evidence for a close-by condensation of one-triplon and two-triplon singlet bound states, suggesting that between the dimer and the collinear state additional phases may intervene.We study a bilayer of the frustrated ${J}_{1}\text{\ensuremath{-}}{J}_{2}$ Heisenberg model on the square lattice. Starting from the dimer limit at strong interlayer coupling, we perform series expansions using the perturbative continuous unitary transformation, based on the flow equation method, in order to determine the spectrum up to the two-triplon sector. From the one-triplon dispersion, we obtain quantum critical lines for transitions from the dimer phase into either N\'eel or collinear magnetic order. For low to intermediate frustration these transitions are consistent with existing findings, based on the magnetic phases. In the region of strongest frustration, i.e., ${J}_{2}/{J}_{1}\ensuremath{\sim}0.5$, we provide an estimate for the stability of the anticipated single-layer quantum spin-liquids against finite interlayer coupling. In the two-triplon sector, we find a set of well defined (anti)bound states, which can be classified according to total spin and in-plane rotational symmetry. For vanishing frustration, these states agree with the previous series expansion analysis. For ${J}_{2}/{J}_{1}\ensuremath{\gtrsim}0.5$, we provide evidence for a close-by condensation of one-triplon and two-triplon singlet bound states, suggesting that between the dimer and the collinear state additional phases may intervene.

Physical review. B./Physical review. BSep 14, 20230 citations

Series expansion studies of the $J_1$ - $J_2$ -Heisenberg bilayerSeries expansion studies of the $J_1$-$J_2$-Heisenberg bilayer

Wagner, Erik, Wolfram Brenig

We study a bilayer of the frustrated $J_1$ - $J_2$ Heisenberg-model on the square lattice. Starting from the dimer limit at strong interlayer coupling, we perform series expansions using the perturbative Continuous Unitary Transformation, based on the flow equation method, in order to determine the spectrum up to the two-triplon sector. From the one-triplon dispersion we obtain quantum critical lines for transitions from the dimer phase into either Néel or collinear magnetic order. For low to intermediate frustration these transitions are consistent with existing findings, based on the magnetic phases. In the region of strongest frustration, i.e. $J_2/J_1\sim 0.5$ , we provide an estimate for the stability of the anticipated single-layer quantum spin-liquids against finite interlayer coupling. In the two-triplon sector we find a set of well defined (anti)bound states, which can be classified according to total spin and in-plane rotational symmetry. For vanishing frustration these states agree with previous series expansion analysis. For $J_2/J_1\gtrsim 0.5$ we provide evidence for a close-by condensation of one-triplon and two-triplon singlet bound states, suggesting that between the dimer and the collinear state, additional phases may intervene.We study a bilayer of the frustrated $J_1$-$J_2$ Heisenberg-model on the square lattice. Starting from the dimer limit at strong interlayer coupling, we perform series expansions using the perturbative Continuous Unitary Transformation, based on the flow equation method, in order to determine the spectrum up to the two-triplon sector. From the one-triplon dispersion we obtain quantum critical lines for transitions from the dimer phase into either Néel or collinear magnetic order. For low to intermediate frustration these transitions are consistent with existing findings, based on the magnetic phases. In the region of strongest frustration, i.e. $J_2/J_1\sim 0.5$, we provide an estimate for the stability of the anticipated single-layer quantum spin-liquids against finite interlayer coupling. In the two-triplon sector we find a set of well defined (anti)bound states, which can be classified according to total spin and in-plane rotational symmetry. For vanishing frustration these states agree with previous series expansion analysis. For $J_2/J_1\gtrsim 0.5$ we provide evidence for a close-by condensation of one-triplon and two-triplon singlet bound states, suggesting that between the dimer and the collinear state, additional phases may intervene.

arXiv:2302.11599Feb 22, 20230 citations

{J}_{2}

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.

Quantum phase diagram of a frustrated spin-

\frac{1}{2}

J_{1}

J_{2}

J_{1}^{\perp}

Heisenberg model on an

AA

-stacked bilayer honeycomb lattice

The zero-temperature quantum phase diagram of the spin-

\frac{1}{2}

J_{1}

J_{2}

J_{1}^{\perp}

model on an

AA

-stacked bilayer honeycomb lattice is investigated using the coupled cluster method (CCM). The model comprises two monolayers in each of which the spins, residing on honeycomb-lattice sites, interact via both nearest-neighbor (NN) and frustrating next-nearest-neighbor isotropic antiferromagnetic (AFM) Heisenberg exchange iteractions, with respective strengths

J_{1} > 0

and

J_{2} \equiv \kappa J_{1}>0

. The two layers are coupled via a comparable Heisenberg exchange interaction between NN interlayer pairs, with a strength

J_{1}^{\perp} \equiv \delta J_{1}

. The complete phase boundaries of two quasiclassical collinear AFM phases, namely the Neel and Neel-II phases, are calculated in the

\kappa \delta

half-plane with

\kappa > 0

. Whereas on each monolayer in the Neel state all NN pairs of spins are antiparallel, in the Neel-II state NN pairs of spins on zigzag chains along one of the three equivalent honeycomb-lattice directions are antiparallel, while NN interchain spins are parallel. We calculate directly in the thermodynamic (infinite-lattice) limit both the magnetic order parameter

M

and the excitation energy

\Delta

from the

s^{z}_{T}=0

ground state to the lowest-lying

|s^{z}_{T}|=1

excited state (where

s^{z}_{T}

is the total

z

component of spin for the system as a whole, and where the collinear ordering lies along the

z

direction) for both quasiclassical states used (separately) as the CCM model state, on top of which the multispin quantum correlations are then calculated to high orders (

n \leq 10

) in a systematic series of approximations involving

n

-spin clusters. The sole approximation made is then to extrapolate the sequences of

n

th-order results for

M

and

\Delta

to the exact limit,

n \to \infty

We study the ground-state phase diagram of the spin-

1/2

Kitaev-Heisenberg model on the bilayer honeycomb lattice with large-scale tensor network calculations based on the infinite projected entangled pair state technique as well as high-order series expansions. We find that beyond various magnetically ordered phases, including ferromagnetic, zigzag, antiferromagnetic (AFM) and stripy states, two extended quantum spin liquid phases arise in the proximity of the Kitaev limit. While these ordered phases also appear in the monolayer Kitaev-Heisenberg model, our results further show that a valence bond solid state emerges in a relatively narrow range of parameter space between the AFM and stripy phases, which can be adiabatically connected to isolated Heisenberg dimers. Our results highlight the importance of considering interlayer interactions on the emergence of novel quantum phases in the bilayer Kitaev materials.

Series expansion studies of the

J_1

J_2

-Heisenberg bilayer

We study a bilayer of the frustrated

J_1

J_2

Heisenberg-model on the square lattice. Starting from the dimer limit at strong interlayer coupling, we perform series expansions using the perturbative Continuous Unitary Transformation, based on the flow equation method, in order to determine the spectrum up to the two-triplon sector. From the one-triplon dispersion we obtain quantum critical lines for transitions from the dimer phase into either Néel or collinear magnetic order. For low to intermediate frustration these transitions are consistent with existing findings, based on the magnetic phases. In the region of strongest frustration, i.e.

J_2/J_1\sim 0.5

, we provide an estimate for the stability of the anticipated single-layer quantum spin-liquids against finite interlayer coupling. In the two-triplon sector we find a set of well defined (anti)bound states, which can be classified according to total spin and in-plane rotational symmetry. For vanishing frustration these states agree with previous series expansion analysis. For

J_2/J_1\gtrsim 0.5

we provide evidence for a close-by condensation of one-triplon and two-triplon singlet bound states, suggesting that between the dimer and the collinear state, additional phases may intervene.