• Spontaneous discrete symmetry breaking can be described in the framework of Projected Entangled Pair States (PEPS) by linearly superposing local tensors belonging to two (or more) symmetry classes of tensors. This is illustrated in the case of a frustrated spin-1 Heisenberg model on the square lattice, which hosts a nematic spin liquid spontaneously breaking lattice $\pi/2$-rotation symmetry. A superposition of SU(2)-symmetric PEPS tensors belonging to two irreducible representations of the lattice point group is shown to capture accurately the properties of the nematic phase, as shown from a comparison to Exact Diagonalisations and Density Matrix Renormalization Group results.
  • We use numerical techniques to study dynamical properties at finite temperature ($T$) of the Heisenberg spin chain with random exchange couplings, which realizes the random singlet (RS) fixed point in the low-energy limit. Specifically, we study the dynamic spin structure factor $S(q,\omega)$, which can be probed directly by inelastic neutron scattering experiments and, in the limit of small $\omega$, in nuclear magnetic resonance (NMR) experiments through the spin-lattice relaxation rate $1/T_1$. Our work combines three complementary methods: exact diagonalization, matrix-product-state algorithms, and stochastic analytic continuation of quantum Monte Carlo results in imaginary time. Unlike the uniform system, whose low-energy excitations at low $T$ are restricted to $q$ close to $0$ and $\pi$, our study reveals a continuous narrow band of low-energy excitations in $S(q,\omega)$, extending throughout the Brillouin zone. Close to $q=\pi$, the scaling properties of these excitations are well captured by the RS theory, but we also see disagreements with some aspects of the predicted $q$-dependence further away from $q=\pi$. Furthermore we find spin diffusion effects close to $q=0$ that are not contained within the RS theory but give non-negligible contributions to the mean $1/T_1$. To compare with NMR experiments, we consider the distribution of the local $1/T_1$ values, which is broad, approximately described by a stretched exponential. The mean value first decreases with $T$, but starts to increase and diverge below a crossover temperature. Although a similar divergent behavior has been found for the static uniform susceptibility, this divergent behavior of $1/T_1$ has never been seen in experiments. Our results show that the divergence of the mean $1/T_1$ is due to rare events in the disordered chains and is concealed in experiments, where the typical $1/T_1$ value is accessed.
  • Motivated by near-term experiments with ultracold alkaline-earth atoms confined to optical lattices, we establish numerically and analytically the phase diagram of two-leg SU($N$) spin ladders. Two-leg ladders provide a rich and highly non-trivial extension of the single chain case on the way towards the relatively little explored two dimensional situation. Focusing on the experimentally relevant limit of one fermion per site, antiferromagnetic exchange interactions, and $2\leq N \leq 6$, we show that the phase diagrams as a function of the interchain (rung) to intrachain (leg) coupling ratio $J_\perp/J_\Vert$ strongly differ for even vs. odd $N$. For even $N=4$ and 6, we demonstrate that the phase diagram consists of a single valence bond crystal (VBC) with a spatial period of $N/2$ rungs. For odd $N=3$ and 5, we find surprisingly rich phase diagrams exhibiting three distinct phases. For weak rung coupling, we obtain a VBC with a spatial period of $N$ rungs, whereas for strong coupling we obtain a critical phase related to the case of a single chain. In addition, we encounter intermediate phases for odd $N$, albeit of a different nature for $N=3$ as compared to $N=5$. For $N=3$, we find a novel gapless intermediate phase with $J_\perp$-dependent incommensurate spatial fluctuations in a sizeable region of the phase diagram. For $N=5$, there are strong indications for a narrow potentially gapped intermediate phase, whose nature is not entirely clear.
  • We study the dynamic spin structure factor of the spin-$1/2$ square-lattice Heisenberg antiferromagnet and of the $J$-$Q$ model (with 4-spin interactions $Q$ and Heisenberg exchange $J$). Using an improved method for stochastic analytic continuation of imaginary-time correlation functions computed with QMC simulations, we can treat the sharp ($\delta$-function) contribution from spinwaves (magnons) and a continuum at higher energy. The results for the Heisenberg model agree with neutron scattering experiments on Cu(DCOO)$_2$$\cdot$4D$_2$O, where a broad spectral-weight continuum at $q=(\pi,0)$ was interpreted as deconfined spinons. Our results at $(\pi,0)$ show a similar reduction of the magnon weight and a large continuum, while the continuum is much smaller at $q=(\pi/2,\pi/2)$ (as also seen experimentally). Turning on $Q$, we observe a rapid reduction of the $(\pi,0)$ magnon weight to zero, well before the deconfined quantum phase transition into a spontaneously dimerized state. We re-interpret the picture of deconfined spinons at $(\pi,0)$ in the experiments as nearly deconfined spinons---a precursor to deconfined quantum criticality. To further elucidate the picture of a fragile $(\pi,0)$-magnon in the Heisenberg model and its depletion in the $J$-$Q$ model, we introduce an effective model in which a magnon can split into two spinons that do not separate but fluctuate in and out of the magnon space (in analogy with the resonance between a photon and a particle-hole pair in the exciton-polariton problem). The model reproduces the $(\pi,0)$ and $(\pi/2,\pi/2)$ features of the Heisenberg model. It can also account for the rapid loss of the $(\pi,0)$ magnon with increasing $Q$ and a remarkable persistence of a large magnon pole at $q=(\pi/2,\pi/2)$ even at the deconfined critical point.
  • We study the interplay between disorder and interactions for emergent bosonic degrees of freedom induced by an external magnetic field in the Br-doped spin-gapped antiferromagnetic material Ni(Cl$_{1-x}$Br$_x$)$_2$-4SC(NH$_2$)$_2$ (DTNX). Building on nuclear magnetic resonance experiments at high magnetic field [A. Orlova et al., Phys. Rev. Lett. 118, 067203 (2017)], we describe the localization of isolated impurity states, providing a realistic theoretical modeling for DTNX. Going beyond single impurity localization we use quantum Monte Carlo simulations to explore many-body effects from which pairwise effective interactions lead to a (impurity-induced) BEC revival [M. Dupont, S. Capponi, and N. Laflorencie, Phys. Rev. Lett. 118, 067204 (2017)]. We further address the question of the existence of a many-body localized Bose-glass (BG) phase in DTNX, which is found to compete with a series of a new kind of BEC regimes made out of the multi-impurity states. The global magnetic field-temperature phase diagram of DTNX reveals a very rich structure for low impurity concentration, with consecutive disorder-induced BEC mini-domes separated by intervening many-body localized BG regimes. Upon increasing the impurity level, multiple mini-BEC phases start to overlap, while intermediate BG regions vanish.
  • Using large-scale simulations based on matrix product state and quantum Monte Carlo techniques, we study the superfluid to Bose glass-transition for one-dimensional attractive hard-core bosons at zero temperature, across the full regime from weak to strong disorder. As a function of interaction and disorder strength, we identify a Berezinskii-Kosterlitz-Thouless critical line with two different regimes. At small attraction where critical disorder is weak compared to the bandwidth, the critical Luttinger parameter $K_c$ takes its universal Giamarchi-Schulz value $K_{c}=3/2$. Conversely, a non-universal $K_c>3/2$ emerges for stronger attraction where weak-link physics is relevant. In this strong disorder regime, the transition is characterized by self-similar power-law distributed weak links with a continuously varying characteristic exponent $\alpha$.
  • Building on recent NMR experiments [Phys. Rev. Lett. 118, 067203], we theoretically investigate the high magnetic field regime of the disordered quasi-one-dimensional $S=1$ antiferromagnetic material Ni(Cl$_{1-x}$Br$_x$)$_2$-4SC(NH$_2$)$_2$. The interplay between disorder, chemically controlled by Br-doping, interactions, and the external magnetic field, leads to a very rich phase diagram. Beyond the well-known antiferromagnetically ordered regime, analog of a Bose condensate of magnons, which disappears when $H\ge 12.3$ T, we unveil a resurgence of phase coherence at higher field $H\sim 13.6$ T, induced by the doping. Interchain couplings stabilize finite temperature long-range order whose extension in the field - temperature space is governed by the concentration of impurities $x$. Such a "mini-condensation" contrasts with previously reported Bose-glass physics in the same regime, and should be accessible to experiments.
  • By measuring the nuclear magnetic resonance (NMR) 1/T_1 relaxation rate in the Br (bond) doped DTN compound, Ni(Cl(1-x)Br_x)2-4SC(NH2)2 (DTNX), we show that the low-energy spin dynamics of its high magnetic field "Bose-glass" regime is dominated by a strong peak of spin fluctuations found at the nearly doping-independent position H* = 13.6 T. From its temperature and field dependence we conclude that this corresponds to a level crossing of the energy levels related to the doping-induced impurity states. Observation of the local NMR signal from the spin adjacent to the doped Br allowed us to fully characterize this impurity state. We have thus quantified a microscopic theoretical model that paves the way to better understanding of the Bose-glass physics in DTNX, as revealed in the related theoretical study [M. Dupont, S. Capponi, and N. Laflorencie, Phys. Rev. Lett. 118, 067204 (2017), arXiv:1610.05136].
  • Nuclear magnetic resonance (NMR) study of the high magnetic field (H) part of the Bose-Einstein condensed (BEC) phase of the quasi-onedimensional (quasi-1D) antiferromagnetic quantum spin-chain compound NiCl2-4SC(NH2)2 (DTN) was performed. We precisely determined the phase boundary, Tc(H), down to 40 mK; the critical boson density, n_c(Tc); and the absolute value of the BEC order parameter S_perp at very low temperature (T = 0.12 K). All results are accurately reproduced by numerical quantum Monte Carlo simulations of a realistic three-dimensional (3D) model Hamiltonian. Approximate analytical predictions based on the 1D Tomonaga-Luttinger liquid description are found to be precise for Tc(H), but less so for S_perp(H), which is more sensitive to the strength of 3D couplings, in particular close to the critical field. A mean-field treatment, based on the Hartree-Fock-Popov description, is found to be valid only up to n_c = 4% (T < 0.3 K), while for higher n_c boson interactions appear to modify the density of states.
  • We present numerical evidence that the spin-1/2 Heisenberg model on the two-dimensional checkerboard lattice exhibits several magnetization plateaux for m=0, 1/4, 1/2 and 3/4, where m is the magnetization normalized by its saturation value. These incompressible states correspond to somehow similar valence-bond crystal phases that break lattice symmetries, though they are different from the already established plaquette phase for m=0. Our results are based on Exact Diagonalization as well as Density Matrix Renormalization Group large-scale simulations, and interpreted in terms of simple trial wavefunctions.
  • Fermions hopping on a hexagonal lattice represent one of the most active research field in condensed matter since the discovery of graphene in 2004 and its numerous applications. Another exciting aspect of the interplay between geometry and quantum mechanical effects is given by the Haldane model (F. D. M. Haldane, Phys. Rev. Lett. 61, 2015 (1988)) where spinless fermions experiencing a certain flux pattern on the honeycomb lattice leads to the stabilization of a topological phase of matter, distinct from a Mott insulator and dubbed Chern insulator nowadays. In this context, it is crucial to understand the role of interactions and this review will describe recent results that have been obtained for a minimal model, namely spinless fermions fermions with nearest and next-nearest neighbour density-density interactions on the honeycomb lattice at half-filling. Topics addressed include an introduction of the minimal model and a discussion of the possible instabilities of the Dirac semimetal, a presentation of various theroretical and numerical approaches, and a summary of the results with a particular emphasis on the stability or not of some exotic quantum phases such as charge ordered ones (similar to Wigner crystals) or spontaneous Chern insulator phase.
  • We study the quantum phase diagram of a system of hard-core bosons on the Kagome lattice with nearest-neighbor repulsive interactions, for arbitrary densities, by means of the hierarchical mean field theory and exact diagonalization techniques. This system is isomorphic to the spin S=1/2 XXZ model in presence of an external magnetic field, a paradigmatic example of frustrated quantum magnetism. In the non-frustrated regime, we find two crystal phases at densities 1/3 and 2/3 that melt into a superfluid phase when increasing the hopping amplitude, in semi-quantitative agreement with quantum Monte Carlo computations. In the frustrated regime and away from half-filling, we find a series of plateaux with densities commensurate with powers of 1/3. The broader density plateaux (at densities 1/3 and 2/3) are remnants of the classical degeneracy in the Ising limit. For densities near half-filling, this staircase of crystal phases melts into a superfluid, which displays finite chiral currents when computed with clusters having an odd number of sites. Both the staircase of crystal phases and the superfluid phase prevail in the non-interacting limit, suggesting that the lowest dispersionless single-particle band may be at the root of this phenomenon.
  • We present results of numerical simulations performed on one-dimensional spin chains in order to extract the so-called relaxation rate $1/T_1$ accessible through NMR experiments. Building on numerical tensor network methods using the Matrix Product States (MPS) formalism, we can follow the non-trivial crossover occurring in critical chains between the high-temperature diffusive classical regime and the low-temperature response described by the Tomonaga-Luttinger liquid (TLL) theory, for which analytical expressions are known. In order to compare analytics and numerics, we focus on a generic spin-$1/2$ XXZ chain which is a paradigm of gapless TLL, as well as a more realistic spin-$1$ anisotropic chain, modelling the DTN material, which can be either in a trivial gapped phase or in a TLL regime induced by an external magnetic field. Thus, by monitoring the finite temperature crossover, we provide quantitative limits on the range of validity of TLL theory, that will be useful when interpreting experiments on quasi one-dimensional materials.
  • Spontaneous symmetry breaking is deeply related to dimensionality of system. The Neel order going with spontaneous breaking of $U(1)$ symmetry is safely allowed at any temperature for three-dimensional systems but allowed only at zero temperature for purely two-dimensional systems. We closely investigate how smoothly the ordering process of the three-dimensional system is modulated into that of the two-dimensional one with reduction of dimensionality, considering spatially anisotropic quantum antiferromagnets. We first show that the N\'eel temperature is kept finite even in the two-dimensional limit although the N\'eel order is greatly suppressed for low-dimensionality. This feature of the N\'eel temperature is highly nontrivial, which dictates how the order parameter is squashed under the reduction of dimensionality. Next we investigate this dimensional modulation of the order parameter. We develop our argument taking as example a coupled spin-ladder system relevant for experimental studies. The ordering process is investigated multidirectionally using theoretical techniques of a mean-field method combined with analytical (exact solutions of quantum field theories) or numerial (density-matrix renormalization-group) method, a variational method, a renormalization-group study, linear spin-wave theory, and quantum Monte-Carlo simulation. We show that these methods independent of each other lead to the same conclusion about the dimensional modulation.
  • We consider the spin-1/2 hamiltonians proposed by Cano and Fendley [J. Cano and P. Fendley, Phys. Rev. Lett. 105, 067205 (2010)] which were built to promote the well-known Rokshar-Kivelson (RK) point of quantum dimer models to spin-1/2 wavefunctions. We first show that these models, besides the exact degeneracy of RK point, support gapless spinless excitations as well as a spin gap in the thermodynamic limit, signatures of an unusual spin liquid. We then extend the original construction to create a continuous family of SU(2) invariant spin models that reproduces the phase diagram of the quantum dimer model, and in particular show explicit evidences for existence of columnar and staggered phases. The original models thus appear as multicritical points in an extended phase diagram. Our results are based on the use of a combination of numerical exact simulations and analytical mapping to effective generalized quantum dimer models.
  • We investigate the phase diagram of spinless fermions with nearest and next-nearest neighbour density-density interactions on the honeycomb lattice at half-filling. Using Exact Diagonalization techniques of the full Hamiltonian and constrained subspaces, combined with a careful choice of finite-size clusters, we determine the different charge orderings that occur for large interactions. In this regime we find a two-sublattice N\'eel-like state, a charge modulated state with a tripling of the unit cell, a zig-zag phase and a novel charge ordered states with a 12 site unit cells we call N\'eel domain wall crystal, as well as a region of phase separation for attractive interactions. A sizeable region of the phase diagram is classically degenerate, but it remains unclear whether an order-by-disorder mechanism will lift the degeneracy. For intermediate repulsion we find evidence for a Kekul\'e or plaquette bond-order wave phase. We also investigate the possibility of a spontaneous Chern insulator phase (dubbed topological Mott insulator), as previously put forward by several mean-field studies. Although we are unable to detect convincing evidence for this phase based on energy spectra and order parameters, we find an enhancement of current-current correlations with the expected spatial structure compared to the non-interacting situation. While for the studied $t{-}V_1{-}V_2$ model the phase transition to the putative topological Mott insulator is preempted by the phase transitions to the various ordered states, our findings might hint at the possibility for a topological Mott insulator in an enlarged Hamiltonian parameter space, where the competing phases are suppressed.
  • We investigate the low-energy properties as well as quantum and thermal phase transitions of the Baxter-Wu model in a transverse magnetic field. Our study relies on stochastic series expansion quantum Monte Carlo and on series expansions about the low- and high-field limits at zero temperature using the quantum finite-lattice method on the triangular lattice. The phase boundary consists of a second-order critical line in the 4-state Potts model universality class starting from the pure Baxter-Wu limit meeting a first-order line connected to the zero-temperature transition point ($h\approx2.4$, $T=0$). Both lines merge at a tricritical point approximatively located at ($h\approx 2.3J$, $T\approx J$).
  • We clarify the existence of several magnetization plateaux for the kagome $S=1/2$ antiferromagnetic Heisenberg model in a magnetic field. Using approximate or exact localized magnon eigenstates, we are able to describe in a similar manner the plateau states that occur for magnetization per site $m=1/3$, $5/9$, and $7/9$ of the saturation value. These results are confirmed using large-scale Exact Diagonalization on lattices up to 63 sites.
  • The Kagome Heisenberg antiferromagnet is mapped onto an effective Hamiltonian on the star superlattice by Contractor Renormalization. Comparison of ground state energies on large lattices to Density Matrix Renormalization Group justifies truncation of effective interactions at range 3. Within our accuracy, magnetic and translational symmetries are not broken (i.e. a spin liquid ground state). However, we discover doublet spectral degeneracies which signal the onset of p6 - chirality symmetry breaking. This is understood by simple mean field analysis. Experimentally, the p6 chiral order parameter should split the optical phonons degeneracy near the zone center. Addition of weak next to nearest neighbor coupling is discussed.
  • In a recent letter [D. Poletti et al., EPL 93, 37008 (2011)] a model of attractive spinless fermions on the honeycomb lattice at half filling has been studied by mean-field theory, where distinct homogenous phases at rather large attraction strength $V>3.36$, separated by (topological) phase transitions, have been predicted. In this comment we argue that without additional interactions the ground states in these phases are not stable against phase separation. We determine the onset of phase separation at half filling $V_{ps}\approx 1.7$ by means of infinite projected entangled-pair states (iPEPS) and exact diagonalization.
  • The phase diagram of one-dimensional earth-alkaline fermionic atoms and ytterbium 171 atoms is investigated by means of a low-energy approach and density-matrix renormalization group calculations. For incommensurate filling, four gapless phases with a spin gap are found and consist of two superconducting instabilities and two coexisting bond and charge density-waves instabilities. In the half-filled case, seven Mott-insulating phases arise with the emergence of four non-degenerate phases with exotic hidden orderings.
  • We extend the definition of the recently introduced valence bond entanglement entropy to arbitrary SU(2) wave functions of S=1/2 spin systems. Thanks to a reformulation of this entanglement measure in terms of a projection, we are able to compute it with various numerical techniques for frustrated spin models. We provide extensive numerical data for the one-dimensional J1-J2 spin chain where we are able to locate the quantum phase transition by using the scaling of this entropy with the block size. We also systematically compare with the scaling of the von Neumann entanglement entropy. We finally underline that the valence-bond entropy definition does depend on the choice of bipartition so that, for frustrated models, a "good" bipartition should be chosen, for instance according to the Marshall sign.
  • The behavior of the ground-state fidelity susceptibility in the vicinity of a quantum critical point is investigated. We derive scaling relations describing its singular behavior in the quantum critical regime. Unlike it has been found in previous studies, these relations are solely expressed in terms of conventional critical exponents. We also describe in detail a quantum Monte Carlo scheme that allows for the evaluation of the fidelity susceptibility for a large class of many-body systems and apply it in the study of the quantum phase transition for the transverse-field Ising model on the square lattice. Finite size analysis applied to the so obtained numerical results confirm the validity of our scaling relations. Furthermore, we analyze the properties of a closely related quantity, the ground-state energy's second derivative, that can be numerically evaluated in a particularly efficient way. The usefulness of both quantities as alternative indicators of quantum criticality is examined.
  • A quantum phase transition is typically induced by tuning an external parameter that appears as a coupling constant in the Hamiltonian. Another route is to vary the global symmetry of the system, generalizing, e.g., SU(2) to SU(N). In that case, however, the discrete nature of the control parameter prevents one from identifying and characterizing the transition. We show how this limitation can be overcome for the SU(N) Heisenberg model with the help of a singlet projector algorithm that can treat N continuously. On the square lattice, we find a direct, continuous phase transition between Neel-ordered and crystalline bond-ordered phases at Nc=4.57(5) with critical exponents z=1 and beta/nu=0.81(3).
  • When a system undergoes a quantum phase transition, the ground-state wave-function shows a change of nature, which can be monitored using the fidelity concept. We introduce two Quantum Monte Carlo schemes that allow the computation of fidelity and its susceptibility for large interacting many-body systems. These methods are illustrated on a two-dimensional Heisenberg model, where fidelity estimators show marked behaviours at two successive quantum phase transitions. We also develop a scaling theory which relates the divergence of the fidelity susceptibility to the critical exponent of the correlation length. A good agreement is found with the numerical results.