
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 spin1 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 lowenergy 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 spinlattice relaxation rate $1/T_1$. Our work combines three
complementary methods: exact diagonalization, matrixproductstate algorithms,
and stochastic analytic continuation of quantum Monte Carlo results in
imaginary time. Unlike the uniform system, whose lowenergy excitations at low
$T$ are restricted to $q$ close to $0$ and $\pi$, our study reveals a
continuous narrow band of lowenergy 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 nonnegligible 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 nearterm experiments with ultracold alkalineearth atoms
confined to optical lattices, we establish numerically and analytically the
phase diagram of twoleg SU($N$) spin ladders. Twoleg ladders provide a rich
and highly nontrivial 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$ squarelattice
Heisenberg antiferromagnet and of the $J$$Q$ model (with 4spin interactions
$Q$ and Heisenberg exchange $J$). Using an improved method for stochastic
analytic continuation of imaginarytime 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 spectralweight 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 reinterpret the picture of deconfined spinons at $(\pi,0)$
in the experiments as nearly deconfined spinonsa 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 particlehole pair in the
excitonpolariton 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 Brdoped
spingapped antiferromagnetic material Ni(Cl$_{1x}$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 manybody effects from which pairwise
effective interactions lead to a (impurityinduced) 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 manybody localized Boseglass (BG)
phase in DTNX, which is found to compete with a series of a new kind of BEC
regimes made out of the multiimpurity states. The global magnetic
fieldtemperature phase diagram of DTNX reveals a very rich structure for low
impurity concentration, with consecutive disorderinduced BEC minidomes
separated by intervening manybody localized BG regimes. Upon increasing the
impurity level, multiple miniBEC phases start to overlap, while intermediate
BG regions vanish.

Using largescale simulations based on matrix product state and quantum Monte
Carlo techniques, we study the superfluid to Bose glasstransition for
onedimensional attractive hardcore bosons at zero temperature, across the
full regime from weak to strong disorder. As a function of interaction and
disorder strength, we identify a BerezinskiiKosterlitzThouless 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 GiamarchiSchulz value $K_{c}=3/2$. Conversely, a nonuniversal
$K_c>3/2$ emerges for stronger attraction where weaklink physics is relevant.
In this strong disorder regime, the transition is characterized by selfsimilar
powerlaw 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
quasionedimensional $S=1$ antiferromagnetic material
Ni(Cl$_{1x}$Br$_x$)$_2$4SC(NH$_2$)$_2$. The interplay between disorder,
chemically controlled by Brdoping, interactions, and the external magnetic
field, leads to a very rich phase diagram. Beyond the wellknown
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 longrange order whose extension in the field 
temperature space is governed by the concentration of impurities $x$. Such a
"minicondensation" contrasts with previously reported Boseglass 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(1x)Br_x)24SC(NH2)2 (DTNX), we show
that the lowenergy spin dynamics of its high magnetic field "Boseglass"
regime is dominated by a strong peak of spin fluctuations found at the nearly
dopingindependent 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 dopinginduced 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 Boseglass
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 BoseEinstein condensed (BEC) phase of the quasionedimensional (quasi1D)
antiferromagnetic quantum spinchain compound NiCl24SC(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 threedimensional (3D) model Hamiltonian. Approximate analytical
predictions based on the 1D TomonagaLuttinger 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
meanfield treatment, based on the HartreeFockPopov 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 spin1/2 Heisenberg model on the
twodimensional 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
valencebond 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
largescale 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 nextnearest neighbour densitydensity interactions on the honeycomb
lattice at halffilling.
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 hardcore bosons on the
Kagome lattice with nearestneighbor 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 nonfrustrated 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 semiquantitative agreement with quantum
Monte Carlo computations. In the frustrated regime and away from halffilling,
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 halffilling, 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
noninteracting limit, suggesting that the lowest dispersionless
singleparticle band may be at the root of this phenomenon.

We present results of numerical simulations performed on onedimensional spin
chains in order to extract the socalled 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 nontrivial crossover
occurring in critical chains between the hightemperature diffusive classical
regime and the lowtemperature response described by the TomonagaLuttinger
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 onedimensional 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 threedimensional systems but allowed only at
zero temperature for purely twodimensional systems. We closely investigate how
smoothly the ordering process of the threedimensional system is modulated into
that of the twodimensional one with reduction of dimensionality, considering
spatially anisotropic quantum antiferromagnets. We first show that the N\'eel
temperature is kept finite even in the twodimensional limit although the
N\'eel order is greatly suppressed for lowdimensionality. 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 spinladder system relevant for experimental studies. The
ordering process is investigated multidirectionally using theoretical
techniques of a meanfield method combined with analytical (exact solutions of
quantum field theories) or numerial (densitymatrix renormalizationgroup)
method, a variational method, a renormalizationgroup study, linear spinwave
theory, and quantum MonteCarlo simulation. We show that these methods
independent of each other lead to the same conclusion about the dimensional
modulation.

We consider the spin1/2 hamiltonians proposed by Cano and Fendley [J. Cano
and P. Fendley, Phys. Rev. Lett. 105, 067205 (2010)] which were built to
promote the wellknown RoksharKivelson (RK) point of quantum dimer models to
spin1/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
nextnearest neighbour densitydensity interactions on the honeycomb lattice at
halffilling. Using Exact Diagonalization techniques of the full Hamiltonian
and constrained subspaces, combined with a careful choice of finitesize
clusters, we determine the different charge orderings that occur for large
interactions. In this regime we find a twosublattice N\'eellike state, a
charge modulated state with a tripling of the unit cell, a zigzag 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 orderbydisorder mechanism will lift the
degeneracy. For intermediate repulsion we find evidence for a Kekul\'e or
plaquette bondorder wave phase. We also investigate the possibility of a
spontaneous Chern insulator phase (dubbed topological Mott insulator), as
previously put forward by several meanfield studies. Although we are unable to
detect convincing evidence for this phase based on energy spectra and order
parameters, we find an enhancement of currentcurrent correlations with the
expected spatial structure compared to the noninteracting 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 lowenergy properties as well as quantum and thermal phase
transitions of the BaxterWu model in a transverse magnetic field. Our study
relies on stochastic series expansion quantum Monte Carlo and on series
expansions about the low and highfield limits at zero temperature using the
quantum finitelattice method on the triangular lattice. The phase boundary
consists of a secondorder critical line in the 4state Potts model
universality class starting from the pure BaxterWu limit meeting a firstorder
line connected to the zerotemperature 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 largescale 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 meanfield 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 entangledpair states
(iPEPS) and exact diagonalization.

The phase diagram of onedimensional earthalkaline fermionic atoms and
ytterbium 171 atoms is investigated by means of a lowenergy approach and
densitymatrix 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 densitywaves
instabilities. In the halffilled case, seven Mottinsulating phases arise with
the emergence of four nondegenerate 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 onedimensional J1J2 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 valencebond 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 groundstate 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 manybody systems and apply it in the study of the quantum
phase transition for the transversefield 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 groundstate 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 Neelordered and
crystalline bondordered 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 groundstate
wavefunction 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 manybody
systems. These methods are illustrated on a twodimensional 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.