BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Frédéric Mila

Frédéric Mila appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.str-elcond-mat.quant-gascond-mat.stat-mechquant-phhep-thcond-mat.dis-nncond-mat.mtrl-scicond-mat.othercond-mat.supr-con

Catalog footprint

What is connected

44works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

44 published item(s)

preprint2024arXiv

Critical line of the triangular Ising antiferromagnet in a field from a $C_3$-symmetric corner transfer matrix algorithm

The corner transfer matrix renormalization group (CTMRG) algorithm has been extensively used to investigate both classical and quantum two-dimensional (2D) lattice models. The convergence of the algorithm can strongly vary from model to model depending on the underlying geometry and symmetries, and the presence of algebraic correlations. An important factor in the convergence of the algorithm is the lattice symmetry, which can be broken due to the necessity of mapping the problem onto the square lattice. We propose a variant of the CTMRG algorithm, designed for models with $C_3$-symmetry, which we apply to the conceptually simple yet numerically challenging problem of the triangular lattice Ising antiferromagnet in a field, at zero and low temperatures. We study how the finite-temperature three-state Potts critical line in this model approaches the ground-state Kosterlitz-Thouless transition driven by a reduced field ($h/T$). In this particular instance, we show that the $C_3$-symmetric CTMRG leads to much more precise results than both existing results from exact diagonalization of transfer matrices and Monte Carlo.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Critical properties of the quantum Ashkin-Teller chain with chiral perturbations

We investigate the nature of the phase transitions in the quantum Ashkin-Teller chain in the presence of chiral perturbations. We locate the Lifshitz line separating a region of direct chiral transitions from the region where the transition is through an intermediate floating phase. Furthermore, we identify a small region in the vicinity of the four-state Potts point where chiral perturbations are irrelevant and where the transition remains conformal. Implications to Rydberg atoms experiments are briefly discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Commensurate-incommensurate transition in the chiral Ashkin-Teller model

We investigate the classical chiral Ashkin-Teller model on a square lattice with the corner transfer matrix renormalisation group (CTMRG) algorithm. We show that the melting of the period-4 phase in the presence of a chiral perturbation takes different forms depending on the coefficient of the four-spin term in the Ashkin-Teller model. Close to the clock limit of two decoupled Ising models, the system undergoes a two-step commensurate-incommensurate transition as soon as the chirality is introduced, with an intermediate critical floating phase bounded by a Kosterlitz-Thouless transition at high temperature and a Pokrovsky-Talapov transition at low temperature. By contrast, close to the four-states Potts model, we argue for the existence of a unique commensurate-incommensurate transition that belongs to the chiral universality class, and for the presence of a Lifshitz point where the ordered, disordered and floating phases meet. Finally, we map the whole phase diagram, which turns out to be in qualitative agreement with the 40 year old prediction by Huse and Fisher.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

From SU(2)$_5$ to SU(2)$_3$ Wess-Zumino-Witten transitions in a frustrated spin-5/2 chain

We investigate the properties of a frustrated spin-5/2 chain with next-nearest neighbor two and three-site interactions, with two questions in mind: the nature of the transition into the dimerized phase induced by the three-site interaction, and the possible presence of a critical floating phase at intermediate values of the next-nearest neighbor interaction. We provide strong evidence that the continuous transition into the dimerized phase, which has been found to be generically in the Wess-Zumino-Witten SU(2)$_{2S}$ universality class up to spin $S=2$, is SU(2)$_5$ only at two isolated points of the phase diagram, and that it is SU(2)$_3$ in between, in agreement with the presence of two relevant operators allowed by symmetry for SU(2)$_5$, and with the conservation of the parity of the level index along the renormalization flow between SU(2)$_k$ theories with different values of $k$. We also find that the dimerization induced by the next-nearest neighbor interaction is a three step process, with first a small partially dimerized phase followed by a broad critical floating phase with incommensurate correlations before the fully dimerized phase is reached. Implications for the iron oxide ${\mathrm{Bi}}_{3}{\mathrm{FeMo}}_{2}{\mathrm{O}}_{12}$ are briefly discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Possibility to detect the bound state of the Heisenberg ferromagnetic chain at intermediate temperature

Motivated by the lack of direct evidence with inelastic neutron scattering of the well documented bound state of Heisenberg ferromagnets, we use the time-dependent Thermal Density Matrix Renormalization Group algorithm to study the temperature dependence of the dynamical spin structure factor of Heisenberg ferromagnetic spin chains. For spin-1/2, we show that the bound state appears as a well defined excitation with significant spectral weight in the temperature range $J/12 \lesssim T \lesssim J/3$, pointing to the possibility of detecting it with inelastic neutron scattering near $k=π$ provided the temperature is neither too low nor too high - at low temperature, the spectral weight only grows as $T^{3/2}$, and at high temperature the bound state peak merges with the two-magnon continuum. For spin-1, the situation is more subtle because the bound state with two neighboring spin flips competes with an anti-bound state with two spin-flips on the same site. As a consequence, the relative spectral weight of the bound state is smaller than for spin-1/2, and a weak resonance due to the anti-bound state appears in the continuum. A clearer signature of the bound state (resp. anti-bound state) can be obtained if a negative (resp. positive) biquadratic interaction is present.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Spinon confinement and deconfinement in spin-1 chains

Motivated by the rich phase diagrams recently unveiled in frustrated spin-1 chains with competing next-nearest neighbour and four-spin interactions, we investigate the nature of the elementary excitations of spin-1 chains in the vicinity of phase transitions using the matrix-product state (MPS) representation of elementary excitations. First, we show that both spinons and magnons naturally arise in SU(2) invariant spin chains when describing ground states and elementary excitations using MPS. Then we investigate the nature of the elementary excitations across the first-oder transition between the Haldane phase and the topologically trivial next-nearest neighbor Haldane phase. We show explicitly that spinons deconfine at the transition, and we calculate the dispersion of these deconfined spinons with MPS. We also show that, immediately away from the transition line, spinons confine, forming dispersive spinon/anti-spinon bound states on both sides of the transition. Finally, we show that deconfined spinons also appear at the transition between the Haldane phase and the spontaneously dimerized phase when this transition is first order.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Thermal Ising transition in the spin-1/2 J1-J2 Heisenberg model

Using an SU(2) invariant finite-temperature tensor network algorithm, we provide strong numerical evidence in favor of an Ising transition in the collinear phase of the spin-1/2 $J_1-J_2$ Heisenberg model on the square lattice. In units of $J_2$, the critical temperature reaches a maximal value of $T_c/J_2\simeq 0.18$ around $J_2/J_1\simeq 1.0$. It is strongly suppressed upon approaching the zero-temperature boundary of the collinear phase $J_2/J_1\simeq 0.6$, and it vanishes as $1/\log(J_2/J_1)$ in the large $J_2/J_1$ limit, as predicted by Chandra, Coleman and Larkin [Phys. Rev. Lett. 64, 88, 1990]. Enforcing the SU(2) symmetry is crucial to avoid the artifact of finite-temperature SU(2) symmetry breaking of U(1) algorithms, opening new perspectives in the investigation of the thermal properties of quantum Heisenberg antiferromagnets.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Phase diagram of the Shastry-Sutherland Compound SrCu2(BO3)2 under extreme combined conditions of field and pressure

Motivated by the intriguing properties of the Shastry-Sutherland compound SrCu2(BO3)2 under pressure, with a still debated intermediate plaquette phase appearing at around 20 kbar and a possible deconfined critical point at higher pressure upon entering the antiferromagnetic phase, we have investigated its high-field properties in this pressure range using tunnel diode oscillator (TDO) measurements. The two main new phases revealed by these measurements are fully consistent with those identified by infinite Projected Entangled Pair states (iPEPS) calculations of the Shastry-Sutherland model, a 1/5 plateau and a 10 x 2 supersolid. Remarkably, these phases are descendants of the full-plaquette phase, the prominent candidate for the intermediate phase of SrCu2(BO3)2. The emerging picture for SrCu2(BO3)2 is shown to be that of a system dominated by a tendency to an orthorhombic distortion at intermediate pressure, an important constraint on any realistic description of the transition into the antiferromagnetic phase.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Quantum Monte Carlo simulations of highly frustrated magnets in a cluster basis: The two-dimensional Shastry-Sutherland model

Quantum Monte Carlo (QMC) simulations constitute nowadays one of the most powerful methods to study strongly correlated quantum systems, provided that no "sign problem" arises. However, many systems of interest, including highly frustrated magnets, suffer from an average sign that is close to zero in standard QMC simulations. Nevertheless, a possible sign problem depends on the simulation basis, and here we demonstrate how a suitable choice of cluster basis can be used to eliminate or at least reduce the sign problem in highly frustrated magnets that were so far inaccessible to efficient QMC simulations. We focus in particular on the application of a two-spin (dimer)-based QMC method to the thermodynamics of the spin-1/2 Shastry-Sutherland model for SrCu$_2$(BO$_3$)$_2$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Floating, critical and dimerized phases in a frustrated spin-3/2 chain

We study spontaneous dimerization and emergent criticality in a spin-3/2 chain with antiferromagnetic nearest-neighbor $J_1$, next-nearest-neighbor $J_2$ and three-site $J_3$ interactions. In the absence of three-site interaction $J_3$, we provide evidence that the model undergoes a remarkable sequence of three phase transitions as a function of $J_2/J_1$, going successively through a critical commesurate phase, a partially dimerized gapped phase, a critical floating phase with quasi-long-range incommensurate order, to end up in a fully dimerized phase at very large $J_2/J_1$. In the field theory language, this implies that the coupling constant of the marginal operator responsible for dimerization changes sign three times. For large enough $J_3$, the fully dimerized phase is stabilized for all $J_2$, and the phase transitions between the critical phases and this phase are both Wess-Zumino-Witten (WZW) SU(2)$_3$ along part of the boundary and turn first order at some point due to the presence of a marginal operator in the WZW SU(2)$_3$ model. By contrast, the transition between the two dimerized phase is always first order, and the phase transitions between the partially dimerized phase and the critical phases are Kosterlitz-Thouless. Finally, we discuss the intriguing spin-1/2 edge states that emerge in the partially dimerized phase for even chains. Unlike their counterparts in the spin-1 chain, they are not confined and disappear upon increasing $J_2$ in favour of a reorganization of the dimerization pattern.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Haldane Gap of the Three-Box Symmetric $\mathrm{SU}(3)$ Chain

Motivated by the recent generalization of the Haldane conjecture to $\mathrm{SU}(3)$ chains [M. Lajkó et al., Nucl. Phys. B924, 508 (2017)] according to which a Haldane gap should be present for symmetric representations if the number of boxes in the Young diagram is a multiple of three, we develop a density matrix renormalization group algorithm based on standard Young tableaus to study the model with three boxes directly in the representations of the global $\mathrm{SU}(3)$ symmetry. We show that there is a finite gap between the singlet and the symmetric $[3\,0\,0]$ sector $Δ_{[3\,0\,0]}/J = 0.040\pm0.006$ where $J$ is the antiferromagnetic Heisenberg coupling, and we argue on the basis of the structure of the low energy states that this is sufficient to conclude that the spectrum is gapped.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Kibble-Zurek exponent and chiral transition of the period-4 phase of Rydberg chains

Chains of Rydberg atoms have emerged as an amazing playground to study quantum physics in 1D. Playing with inter-atomic distances and laser detuning, one can in particular explore the commensurate-incommensurate transition out of charge-density waves through the Kibble-Zurek mechanism, and the possible presence of a chiral transition with dynamical exponent $z>1$. Here we address this problem theoretically with effective blockade models where the short-distance repulsions are replaced by a constraint of no double occupancy. For the period-4 phase, we show there is an Ashkin-Teller transition point with exponent $ν=0.78$ surrounded by a direct chiral transition with a dynamical exponent $z=1.14$ and a Kibble-Zurek exponent $μ=0.4$. For Rydberg atoms with a van der Waals potential, we suggest that the experimental value $μ=0.25$ is due to a chiral transition with $z\simeq 1.9$ and $ν\simeq 0.47$ surrounding an Ashkin-Teller transition close to the 4-state Potts universality.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Magnetic field induced bound states in spin-$\frac{1}{2}$ ladders

Motivated by the intriguing mode splittings in a magnetic field recently observed with inelastic neutron scattering in the spin ladder compound (C$_5$H$_{12}$N)$_2$CuBr$_4$ (BPCB), we investigate the nature of the spin ladder excitations using DMRG and analytical arguments. Starting from the fully frustrated ladder, for which we derive the low-energy spectrum, we show that bound states are generically present close to $k=0$ in the dynamical structure factor of spin ladders above $H_{c 1}$, and that they are characterized by a field-independent binding energy and an intensity that grows with $H-H_{c 1}$. These predictions are shown to explain quantitatively the split modes observed in BPCB.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Reduction of the sign problem near $T=0$ in quantum Monte Carlo simulations

Building on a recent investigation of the Shastry-Sutherland model [S. Wessel et al., Phys. Rev. B 98, 174432 (2018)], we develop a general strategy to eliminate the Monte Carlo sign problem near the zero temperature limit in frustrated quantum spin models. If the Hamiltonian of interest and the sign-problem-free Hamiltonian---obtained by making all off-diagonal elements negative in a given basis---have the same ground state and this state is a member of the computational basis, then the average sign returns to one as the temperature goes to zero. We illustrate this technique by studying the triangular and kagome lattice Heisenberg antiferrromagnet in a magnetic field above saturation, as well as the Heisenberg antiferromagnet on a modified Husimi cactus in the dimer basis. We also provide detailed appendices on using linear programming techniques to automatically generate efficient directed loop updates in quantum Monte Carlo simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Solving frustrated Ising models using tensor networks

Motivated by the recent success of tensor networks to calculate the residual entropy of spin ice and kagome Ising models, we develop a general framework to study frustrated Ising models in terms of infinite tensor networks %, i.e. tensor networks that can be contracted using standard algorithms for infinite systems. This is achieved by reformulating the problem as local rules for configurations on overlapping clusters chosen in such a way that they relieve the frustration, i.e. that the energy can be minimized independently on each cluster. We show that optimizing the choice of clusters, including the weight on shared bonds, is crucial for the contractibility of the tensor networks, and we derive some basic rules and a linear program to implement them. We illustrate the power of the method by computing the residual entropy of a frustrated Ising spin system on the kagome lattice with next-next-nearest neighbour interactions, vastly outperforming Monte Carlo methods in speed and accuracy. The extension to finite-temperature is briefly discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Spin-$\frac{1}{2}$ kagome Heisenberg antiferromagnet with strong breathing anisotropy

We study the zero-temperature phase diagram of the spin-$\frac{1}{2}$ Heisenberg model with breathing anisotropy (i.e., with different coupling strength on the upward and downward triangles) on the kagome lattice. Our study relies on large scale tensor network simulations based on infinite projected entangled-pair state and infinite projected entangled-simplex state methods adapted to the kagome lattice. Our energy analysis suggests that the U(1) algebraic quantum spin-liquid (QSL) ground-state of the isotropic Heisenberg model is stable up to very large breathing anisotropy until it breaks down to a critical lattice-nematic phase that breaks rotational symmetry in real space through a first-order quantum phase transition. Our results also provide further insight into the recent experiment on vanadium oxyfluoride compounds which has been shown to be relevant platforms for realizing QSL in the presence of breathing anisotropy.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

SrCu$_2$(BO$_3$)$_2$ under pressure: a first-principle study

Using density-functional theory (DFT) band-structure calculations, we study the crystal structure, the lattice dynamics, and the magnetic interactions in the Shastry-Sutherland magnet SrCu$_2$(BO$_3$)$_2$ under pressure, concentrating on experimentally relevant pressures up to 4 GPa. All tendencies inside the candidate intermediate phases are thoroughly worked out, including specific predictions for some Raman active phonon modes that could be used to pin down the nature of the intermediate phase.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Generalization of the Haldane conjecture to SU($n$) chains

Recently, SU(3) chains in the symmetric and self-conjugate representations have been studied using field theory techniques. For certain representations, namely rank-$p$ symmetric ones with $p$ not a multiple of 3, it was argued that the ground state exhibits gapless excitations. For the remaining representations considered, a finite energy gap exists above the ground state. In this paper, we extend these results to SU($n$) chains in the symmetric representation. For a rank-$p$ symmetric representation with $n$ and $p$ coprime, we predict gapless excitations above the ground state. If $p$ is a multiple of $n$, we predict a unique ground state with a finite energy gap. Finally, if $p$ and $n$ have a greatest common divisor $1<q<n$, we predict a ground state degeneracy of $n/q$, with a finite energy gap. To arrive at these results, we derive a non-Lorentz invariant flag manifold sigma model description of the SU($n$) chains, and use the renormalization group to show that Lorentz invariance is restored at low energies. We then make use of recently developed anomaly matching conditions for these Lorentz-invariant models. We also review the Lieb-Shultz-Mattis-Affleck theorem, and extend it to SU($n$) models with longer range interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Perturbative approach to quantum tunneling and quantum interferences in spin clusters

Collective tunneling is a ubiquitous phenomenon in finite-size spin clusters that shows up in systems as diverse as molecular magnets or spin clusters adsorbed at surfaces. The problem we explore is to understand how small flipping terms can cooperate to flip a large spin to the opposite direction or a cluster of interacting elementary Ising spins into the time-reversed state. These high order processes will involve at least two channels, a single spin-flip channel due to a transverse field and a two-spin flip channel due to exchange or due to single-ion anisotropies. In the present paper, we show that high-order perturbation theory can be formulated and evaluated with the help of simple recurrence relations, leading to a compact theory of tunnelling in macroscopic spins, in one-dimensional clusters, as well as in small higher-dimensional clusters. This is demonstrated explicitly in the case of the Ising model with a transverse field and transverse exchange, and in the case of macroscopic spins with uniaxial anisotropy. Our approach provides a transparent theory of level crossings, where the tunneling between time reversed configurations vanishes as a function of the external field. The crossings result from destructive quantum interferences between competing flipping channels. Our theory consistently predicts N crossings in chains of N Ising spins, 2S crossings in single spins of magnitude S, and yields explicit analytical formulae for the level crossings of open chains and macroscopic spins. Disorder can be easily implemented in this perturbative formalism. Leading disorder effects can be treated analytically for spin rings. At the smallest transverse field crossing, the suppression of tunneling is most robust with respect to disorder. We briefly discuss the implications of our findings for the use of realistic spin clusters on surfaces to store information.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Chiral spin liquids in triangular lattice SU(N) fermionic Mott insulators with artificial gauge fields

We show that, in the presence of a $π/2$ artificial gauge field per plaquette, Mott insulating phases of ultra-cold fermions with $SU(N)$ symmetry and one particle per site generically possess an extended chiral phase with intrinsic topological order characterized by a multiplet of $N$ low-lying singlet excitations for periodic boundary conditions, and by chiral edge states described by the $SU(N)_1$ Wess-Zumino-Novikov-Witten conformal field theory for open boundary conditions. This has been achieved by extensive exact diagonalizations for $N$ between $3$ and $9$, and by a parton construction based on a set of $N$ Gutzwiller projected fermionic wave-functions with flux $π/N$ per triangular plaquette. Experimental implications are briefly discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Multi-triplet bound states and finite-temperature dynamics in highly frustrated quantum spin ladders

Low-dimensional quantum magnets at finite temperatures present a complex interplay of quantum and thermal fluctuation effects in a restricted phase space. While some information about dynamical response functions is available from theoretical studies of the one-triplet dispersion in unfrustrated chains and ladders, little is known about the finite-temperature dynamics of frustrated systems. Experimentally, inelastic neutron scattering studies of the highly frustrated two-dimensional material SrCu$_2$(BO$_3$)$_2$ show an almost complete destruction of the one-triplet excitation band at a temperature only 1/3 of its gap energy, accompanied by strong scattering intensities for apparent multi-triplet excitations. We investigate these questions in the frustrated spin ladder and present numerical results from exact diagonalization for the dynamical structure factor as a function of temperature. We find anomalously rapid transfer of spectral weight out of the one-triplet band and into both broad and sharp spectral features at a wide range of energies, including below the zero-temperature gap of this excitation. These features are multi-triplet bound states, which develop particularly strongly near the quantum phase transition, fall to particularly low energies there, and persist to all the way to infinite temperature. Our results offer valuable insight into the physics of finite-temperature spectral functions in SrCu$_2$(BO$_3$)$_2$ and many other highly frustrated spin systems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Semiclassical theory of the magnetization process of the triangular lattice Heisenberg model

Motivated by the numerous examples of 1/3 magnetization plateaux in the triangular lattice Heisenberg an- tiferromagnet with spins ranging from 1/2 to 5/2, we revisit the semiclassical calculation of the magnetization curve of that model, with the aim of coming up with a simple method that allows one to calculate the full mag- netization curve, and not just the critical fields of the 1/3 plateau. We show that it is actually possible to calculate the magnetization curve including the first quantum corrections and the appearance of the 1/3 plateau entirely within linear spin-wave theory, with predictions for the critical fields that agree to order 1/S with those derived a long-time ago on the basis of arguments that required to go beyond linear spin-wave theory. This calculation relies on the central observation that there is a kink in the semiclassical energy at the field where the classical ground state is the collinear up-up-down structure, and that this kink gives rise to a locally linear behavior of the energy with the field when all semiclassical ground states are compared to each other for all fields. The magnetization curves calculated in this way for spin 1/2, 1 and 5/2 are shown to be in good agreement with available experimental data.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Stabilization of the chiral phase of the SU($6m$) Heisenberg model on the honeycomb lattice with $m$ particles per site for $m$ larger than 1

We show that, when $N$ is a multiple of 6 ($N=6m$, $m$ integer), the \SU{N} Heisenberg model on the honeycomb lattice with $m$ particles per site has a clear tendency toward chiral order as soon as $m\geq 2$. This conclusion has been reached by a systematic variational Monte Carlo investigation of Gutzwiller projected wave-functions as a function of $m$ between the case of one particle per site ($m=1$), for which the ground state has recently been shown to be in a plaquette singlet state, and the $m\rightarrow \infty$ limit, where a mean-field approach has established that the ground state has chiral order. This demonstrates that the chiral phase can indeed be stabilized for not too large values of $m$, opening the way to its experimental realisations in other lattices.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

SU(6) Heisenberg model on the honeycomb lattice: competition between plaquette and chiral order

We revisit the SU(6) Heisenberg model on the honeycomb lattice, which has been predicted to be a chiral spin liquid by mean-field theory [G. Szirmai et al., Phys. Rev. A 84, 011611 (2011)]. Using exact diagonalizations of finite clusters, infinite projected entangled pair states simulations, and variational Monte Carlo simulations based on Gutzwiller projected wave functions, we provide strong evidence in favour of the competing plaquette state, which was reported to be higher but close by in energy according to mean-field theory. This is further confirmed by the investigation of the model with a ring exchange term, which shows that there is a transition between the plaquette state and the chiral state at a finite value of the ring exchange term.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Variational Monte-Carlo investigation of SU($N$) Heisenberg chains

Motivated by recent experimental progress in the context of ultra-cold multi-color fermionic atoms in optical lattices, we have investigated the properties of the SU($N$) Heisenberg chain with totally antisymmetric irreducible representations, the effective model of Mott phases with $m < N$ particles per site. These models have been studied for arbitrary $N$ and $m$ with non-abelian bosonization [I. Affleck, Nuclear Physics B 265, 409 (1986); 305, 582 (1988)], leading to predictions about the nature of the ground state (gapped or critical) in most but not all cases. Using exact diagonalization and variational Monte-Carlo based on Gutzwiller projected fermionic wave functions, we have been able to verify these predictions for a representative number of cases with $N \leq 10$ and $m \leq N/2$, and we have shown that the opening of a gap is associated to a spontaneous dimerization or trimerization depending on the value of m and N. We have also investigated the marginal cases where abelian bosonization did not lead to any prediction. In these cases, variational Monte-Carlo predicts that the ground state is critical with exponents consistent with conformal field theory.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Quantum dimer model for the spin-1/2 kagome Z2 spin liquid

We revisit the description of the low-energy singlet sector of the spin-1/2 Heisenberg antiferromagnet on kagome in terms of an effective quantum dimer model. With the help of exact diagonalizations of appropriate finite-size clusters, we show that the embedding of a given process in its kagome environment leads to dramatic modifications of the amplitudes of the elementary loop processes, an effect not accessible to the standard approach based on the truncation of the Hamiltonian to the nearest-neighbour valence-bond basis. The resulting parameters are consistent with a Z$_2$ spin liquid rather than with a valence-bond crystal, in agreement with the last density matrix renormalization group results.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Semiclassical evidence of columnar order in the fully frustrated transverse field Ising model on the square lattice

We investigate the zero-temperature phase diagram of the fully frustrated transverse field Ising model on the square lattice both in the classical limit and in the presence of quantum fluctuations. At the classical level (the limit of infinite spin $S$), we find that upon decreasing the transverse field $Γ$ this model exhibits a phase transition from the fully polarized state into an eight-fold degenerate translational symmetry breaking state. This phase can be identified to correspond to plaquette order in the dimer language and remains the lowest-energy state in the entire range of fields below the critical one, $Γ_c$. The eight-fold degenerate solution which corresponds to columnar order in the dimer language is a saddle point of the classical energy. It is degenerate with the plaquette solution at $Γ=0$ and is only slightly higher in energy in the whole interval $0<Γ<Γ_c$. The effect of quantum fluctuations is investigated in the context of a large S expansion both for the plaquette and columnar structures. For this purpose we employ an approximate method allowing to estimate from above the fluctuation-induced correction to the energy of a configuration which at the classical level is a saddle point of the energy, \textit{not} a local minimum. Although the convergence of the $1/S$ expansion in the $Γ/J\rightarrow 0$ limit remains an open question, harmonic quantum fluctuations show a clear tendency to overcome the energy difference between the two states and to change the classical picture favoring the columnar order over the plaquette one in a wide parameter range.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Berry phase investigation of spin-S ladders

We investigate the properties of antiferromagnetic spin-S ladders with the help of local Berry phases defined by imposing a twist on one or a few local bonds. In gapped systems with time reversal symmetry, these Berry phases are quantized, hence able in principle to characterize different phases. In the case of a fully frustrated ladder where the total spin on a rung is a conserved quantity that changes abruptly upon increasing the rung coupling, we show that two Berry phases are relevant to detect such phase transitions: the rung Berry phase defined by imposing a twist on one rung coupling, and the twist Berry phase defined by twisting the boundary conditions along the legs. In the case of non-frustrated ladders, we have followed the fate of both Berry phases when interpolating between standard ladders and dimerized spin chains. A careful investigation of the spin gap and of edge states shows that a change of twist Berry phase is associated to a quantum phase transition at which the bulk gap closes, and at which, with appropriate boundary conditions, edge states appear or disappear, while a change of rung Berry phase is not necessarily associated to a quantum phase transition. The difference is particularly acute for regular ladders, in which the twist Berry phase does not change at all upon increasing the rung coupling from zero to infinity while the rung Berry phase changes 2S times. By analogy with the fully frustrated ladder, these changes are interpreted as cross-overs between domains in which the rungs are in different states of total spin from 0 in the strong rung limit to 2S in the weak rung limit. This interpretation is further supported by the observation that these cross-overs turn into real phase transitions as a function of rung coupling if one rung is strongly ferromagnetic, or equivalently if one rung is replaced by a spin 2S impurity.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Competing states in the SU(3) Heisenberg model on the honeycomb lattice: Plaquette valence-bond crystal versus dimerized color-ordered state

Conflicting predictions have been made for the ground state of the SU(3) Heisenberg model on the honeycomb lattice: Tensor network simulations found a plaquette order [Zhao et al, Phys. Rev. B 85, 134416 (2012)], where singlets are formed on hexagons, while linear flavor-wave theory (LFWT) suggested a dimerized, color ordered state [Lee and Yang, Phys. Rev. B 85, 100402 (2012)]. In this work we show that the former state is the true ground state by a systematic study with infinite projected-entangled pair states (iPEPS), for which the accuracy can be systematically controlled by the so-called bond dimension $D$. Both competing states can be reproduced with iPEPS by using different unit cell sizes. For small $D$ the dimer state has a lower variational energy than the plaquette state, however, for large $D$ it is the latter which becomes energetically favorable. The plaquette formation is also confirmed by exact diagonalizations and variational Monte Carlo studies, according to which both the dimerized and plaquette states are non-chiral flux states.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Competition between three-sublattice order and superfluidity in the quantum 3-state Potts model of ultracold bosons and fermions on a square optical lattice

We study a quantum version of the three-state Potts model that includes as special cases the effective models of bosons and fermions on the square lattice in the Mott insulating limit. It can be viewed as a model of quantum permutations with amplitudes J_parallel and J_perp for identical and different colors, respectively. For J_parallel=J_perp>0, it is equivalent to the SU(3) Heisenberg model, which describes the Mott insulating phase of 3-color fermions, while the parameter range J_perp<min(0,-J_parallel) can be realized in the Mott insulating phase of 3-color bosonic atoms. Using linear flavor wave theory, infinite projected entangled-pair states (iPEPS), and continuous-time quantum Monte-Carlo simulations, we construct the full T=0 phase diagram, and we explore the T>0 properties for J_perp<0. For dominant antiferromagnetic J_parallel interactions, a three-sublattice long-range ordered stripe state is selected out of the ground state manifold of the antiferromagnetic Potts model by quantum fluctuations. Upon increasing |J_perp|, this state is replaced by a uniform superfluid for J_perp<0, and by an exotic three-sublattice superfluid followed by a two-sublattice superfluid for J_perp>0. The transition out of the uniform superfluid (that can be realized with bosons) is shown to be a peculiar type of Kosterlitz-Thouless transition with three types of elementary vortices.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Realization of higher Wess-Zumino-Witten models in spin chains

Building on the generalization of the exactly dimerized Majumdar-Ghosh ground state to arbitrary spin S for the Heisenberg chain with a three-site term (S_{i-1} S_i)(S_i S_{i+1}) + H.c., we use density-matrix renormalization group simulations and exact diagonalizations to determine the nature of the dimerization transition for S=1, 3/2 and 2. The resulting central charge and critical exponent are in good agreement with the SU(2)_{k=2S} Wess-Zumino-Witten values c=3k/(2+k) and η=3/(2+k). Since the 3-site term that induces dimerization appears naturally if exchange interactions are calculated beyond second order, these results suggest that SU(2)_{k>1} Wess-Zumino-Witten models might finally be realized in actual spin chains.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Zero-temperature Monte Carlo study of the non-coplanar phase of the classical bilinear-biquadratic Heisenberg model on the triangular lattice

We investigate the ground-state properties of the highly degenerate non-coplanar phase of the classical bilinear-biquadratic Heisenberg model on the triangular lattice with Monte Carlo simulations. For that purpose, we introduce an Ising pseudospin representation of the ground states, and we use a simple Metropolis algorithm with local updates, as well as a powerful cluster algorithm. At sizes that can be sampled with local updates, the presence of long-range order is surprisingly combined with an algebraic decay of correlations and the complete disordering of the chirality. It is only thanks to the investigation of unusually large systems (containing $\sim 10^8$ spins) with cluster updates that the true asymptotic regime can be reached and that the system can be proven to consist of equivalent (i.e., equally ordered) sublattices. These large-scale simulations also demonstrate that the scalar chirality exhibits long-range order at zero temperature, implying that the system has to undergo a finite-temperature phase transition. Finally, we show that the average distance in the order parameter space, which has the structure of an infinite Cayley tree, remains remarkably small between any pair of points, even in the limit when the real space distance between them tends to infinity.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Anticollinear magnetic order induced by impurities in the frustrated Heisenberg model of pnictides

We present Monte Carlo simulations for a classical antiferromagnetic Heisenberg model with both nearest ($J_1$) and next-nearest ($J_2$) exchange couplings on the square lattice in the presence of non-magnetic impurities. We show that the order-by-disorder entropy selection, associated with the Ising-like phase transition that appears for $J_2/J_1>1/2$ in the pure spin model, is quenched at low temperature due to the presence of non-magnetic impurities. Evidences that a new competing order is stabilized around the impurities, and in turn induces a re-entrance phase transition are reported. Implications for local magnetic measurement of the parent compound of iron pnictides are briefly discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Antiferromagnetic spin-S chains with exactly dimerized ground states

We show that spin S Heisenberg spin chains with an additional three-body interaction of the form (S_{i-1}S_{i})(S_{i}S_{i+1})+h.c. possess fully dimerized ground states if the ratio of the three-body interaction to the bilinear one is equal to 1/(4S(S+1)-2). This result generalizes the Majumdar-Ghosh point of the J_1-J_2 chain, to which the present model reduces for S=1/2. For S=1, we use the density matrix renormalization group method (DMRG) to show that the transition between the Haldane and the dimerized phases is continuous with central charge c=3/2. Finally, we show that such a three-body interaction appears naturally in a strong-coupling expansion of the Hubbard model, and we discuss the consequences for the dimerization of actual antiferromagnetic chains.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Entropy dependence of correlations in one-dimensional SU(N) antiferromagnets

Motivated by the possibility to load multi-color fermionic atoms in optical lattices, we study the entropy dependence of the properties of the one-dimensional antiferromagnetic SU(N) Heisenberg model, the effective model of the SU(N) Hubbard model with one particle per site (filling 1/N). Using continuous-time world line Monte Carlo simulations for N=2 to 5, we show that characteristic short-range correlations develop at low temperature as a precursor of the ground state algebraic correlations. We also calculate the entropy as a function of temperature, and we show that the first sign of short-range order appears at an entropy per particle that increases with N and already reaches 0.8k_B at N=4, in the range of experimentally accessible values.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Evidence of columnar order in the fully frustrated transverse field Ising model on the square lattice

Using extensive classical and quantum Monte Carlo simulations, we investigate the ground-state phase diagram of the fully frustrated transverse field Ising model on the square lattice. We show that pure columnar order develops in the low-field phase above a surprisingly large length scale, below which an effective U(1) symmetry is present. The same conclusion applies to the Quantum Dimer Model with purely kinetic energy, to which the model reduces in the zero-field limit, as well as to the stacked classical version of the model. By contrast, the 2D classical version of the model is shown to develop plaquette order. Semiclassical arguments show that the transition from plaquette to columnar order is a consequence of quantum fluctuations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Semiclassical approach to ground-state properties of hard-core bosons in two dimensions

Motivated by some inconsistencies in the way quantum fluctuations are included beyond the classical treatment of hard-core bosons on a lattice in the recent literature, we revisit the large-S semi-classical approach to hard-core bosons on the square lattice at T=0. First of all, we show that, if one stays at the purely harmonic level, the only correct way to get the 1/S correction to the density is to extract it from the derivative of the ground state energy with respect to the chemical potential, and that to extract it from a calculation of the ground state expectation value of the particle number operator, it is necessary to include 1/\sqrt{S} corrections to the harmonic ground state. Building on this alternative approach to get 1/S corrections, we provide the first semiclassical derivation of the momentum distribution, and we revisit the calculation of the condensate density. The results of these as well as other physically relevant quantities such as the superfluid density are systematically compared to quantum Monte Carlo simulations. This comparison shows that the logarithmic corrections in the dilute Bose gas limit are only captured by the semi-classical approach if the 1/S corrections are properly calculated, and that the semi-classical approach is able to reproduce the 1/k divergence of the momentum distribution at k=0. Finally, the effect of 1/S^2 corrections is briefly discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Spin-wave analysis of the transverse-field Ising model on the checkerboard lattice

The ground state properties of the S=1/2 transverse-field Ising model on the checkerboard lattice are studied using linear spin wave theory. We consider the general case of different couplings between nearest neighbors (J1) and next-to-nearest neighbors (J2). In zero field the system displays a large degeneracy of the ground state, which is exponential in the system size (for J1=J2) or in the system's linear dimensions (for J2>J1). Quantum fluctuations induced by a transverse field are found to be unable to lift this degeneracy in favor of a classically ordered state at the harmonic level. This remarkable fact suggests that a quantum-disordered ground state can be instead promoted when non-linear fluctuations are accounted for, in agreement with existing results for the isotropic case J1=J2. Moreover spin-wave theory shows sizable regions of instability which are further candidates for quantum-disordered behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Three-sublattice order in the SU(3) Heisenberg model on the square and triangular lattice

We present a numerical study of the SU(3) Heisenberg model of three-flavor fermions on the triangular and square lattice by means of the density-matrix renormalization group (DMRG) and infinite projected entangled-pair states (iPEPS). For the triangular lattice we confirm that the ground state has a three-sublattice order with a finite ordered moment which is compatible with the result from linear flavor wave theory (LFWT). The same type of order has recently been predicted also for the square lattice [PRL 105, 265301 (2010)] from LFWT and exact diagonalization. However, for this case the ordered moment cannot be computed based on LFWT due to divergent fluctuations. Our numerical study clearly supports this three-sublattice order, with an ordered moment of m=0.2-0.4 in the thermodynamic limit.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Phase diagram and continuous pair-unbinding transition of the bilinear-biquadratic S=1 Heisenberg chain in a magnetic field

We investigate the properties of the Heisenberg S=1 chain with bilinear and biquadratic interactions in a magnetic field using the Density Matrix Renormalization Group, Bethe ansatz and field theoretical considerations. In a large region of the parameter space, we identify a magnetized ferroquadrupolar Luttinger liquid consisting of a quasi-condensate of bound magnon pairs. This liquid undergoes a continuous pair unbinding transition to a more conventional Luttinger liquid region obtained by polarizing the system above the Haldane gap region. This pair unbinding transition is shown to be in the Ising universality class on top of a Luttinger liquid, leading to an effective central charge 3/2. We also revisit the nature of the partially polarized Luttinger liquid around and above the Uimin-Lai-Sutherland point. Our results confirm that this is a two-component liquid and rule out the formation of a single-component vector chiral phase.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Phase separation versus supersolid behavior in frustrated antiferromagnets

We investigate the competition between spin-supersolidity and phase separation in a frustrated spin-half model of weakly coupled dimers. We start by considering systems of hard-core bosons on the square lattice, onto which the low-energy physics of the herein investigated spin model can be mapped, and devise a criterion for gauging the interplay between supersolid order and domain wall formation based on strong coupling arguments. Effective bosonic models for the spin model are derived via the contractor renormalization (CORE) algorithm and we propose to combine a self-consistent cluster mean-field solution with our criterion for the occurrence of phase separation to derive the phase diagram as a function of frustration and magnetic field. In the limit of strong frustration, the model is shown to be unstable toward phase separation, in contradiction with recently published results. However, a region of stable supersolidity is identified for intermediate frustration, in a parameter range not investigated so far and of possible experimental relevance.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Theory of inelastic light scattering in spin-1 systems: resonant regimes and detection of quadrupolar order

Motivated by the lack of an obvious spectroscopic probe to investigate non-conventional order such as quadrupolar orders in spin S>1/2 systems, we present a theoretical approach to inelastic light scattering for spin-1 quantum magnets in the context of a two-band Hubbard model. In contrast to the S=1/2 case, where the only type of local excited state is a doubly occupied state of energy $U$, several local excited states with occupation up to 4 electrons are present. As a consequence, we show that two distinct resonating scattering regimes can be accessed depending on the incident photon energy. For $\hbarω_{in}\lesssim U$, the standard Loudon-Fleury operator remains the leading term of the expansion as in the spin-1/2 case. For $\hbarω_{in}\lesssim4U$, a second resonant regime is found with a leading term that takes the form of a biquadratic coupling $\sim({\bf S}_{i}\cdot{\bf S}_{j)^{2}$. Consequences for the Raman spectra of S=1 magnets with magnetic or quadrupolar order are discussed. Raman scattering appears to be a powerful probe of quadrupolar order.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Unconventional magnetization plateaus in a Shastry-Sutherland spin tube

Using density matrix renormalization group (DMRG) and perturbative continuous unitary transformations (PCUTs), we study the magnetization process in a magnetic field for all coupling strengths of a quasi-1D version of the 2D Shastry-Sutherland lattice, a frustrated spin tube made of two orthogonal dimer chains. At small inter-dimer coupling, plateaus in the magnetization appear at 1/6, 1/4, 1/3, 3/8, and 1/2. As in 2D, they correspond to a Wigner crystal of triplons. However, close to the boundary of the product singlet phase, plateaus of a new type appear at 1/5 and 3/4. They are stabilized by the localization of {\it bound states} of triplons. Their magnetization profile differs significantly from that of single triplon plateaus and leads to specific NMR signatures. We address the possibility to stabilize such plateaus in further geometries by analyzing small finite clusters using exact diagonalizations and the PCUTs.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

A 2D quantum antiferromagnet with a four-fold degenerate valence-bond-solid ground state

We study the competition between antiferromagnetic order and valence bond solid formation in a two-dimensional frustrated spin-1/2 model. The J1-J2 model on the square lattice is further frustrated by introducing products of three-spin projectors which stabilize four dimer-product states as degenerate ground states. These four states are reminiscent of the dimerized singlet ground state of the Shastry-Sutherland model. Using exact diagonalizations, we study the evolution of the ground state by varying theratio of interactions. For a large range of parameters (J2 > 0.25J1), our model shows a direct transition between the valence-bond-solid phase and the collinear antiferromagnetic phase. For small values of J2, several intermediate phases appear which are also analyzed.

0 citations0 reviewsNo code linkedOpen access