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

Marcos Rigol

Marcos Rigol appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.quant-gascond-mat.stat-mechcond-mat.str-elquant-phcond-mat.dis-nncond-mat.mes-hallcond-mat.mtrl-scihep-thphysics.atom-phphysics.comp-ph

Catalog footprint

What is connected

77works
10topics
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

77 published item(s)

preprint2022arXiv

Statistical properties of the off-diagonal matrix elements of observables in eigenstates of integrable systems

We study the statistical properties of the off-diagonal matrix elements of observables in the energy eigenstates of integrable quantum systems. They have been found to be dense in the spin-1/2 XXZ chain, while they are sparse in noninteracting systems. We focus on the quasimomentum occupation of hard-core bosons in one dimension, and show that the distributions of the off-diagonal matrix elements are well described by generalized Gamma distributions, in both the presence and absence of translational invariance but not in the presence of localization. We also show that the results obtained for the off-diagonal matrix elements of observables in the spin-1/2 XXZ model are well described by a generalized Gamma distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Volume-law entanglement entropy of typical pure quantum states

The entanglement entropy of subsystems of typical eigenstates of quantum many-body Hamiltonians has been recently conjectured to be a diagnostic of quantum chaos and integrability. In quantum chaotic systems it has been found to behave as in typical pure states, while in integrable systems it has been found to behave as in typical pure Gaussian states. In this tutorial, we provide a pedagogical introduction to known results about the entanglement entropy of subsystems of typical pure states and of typical pure Gaussian states. They both exhibit a leading term that scales with the volume of the subsystem, when smaller than one half of the volume of the system, but the prefactor of the volume law is fundamentally different. It is constant (and maximal) for typical pure states, and it depends on the ratio between the volume of the subsystem and of the entire system for typical pure Gaussian states. Since particle-number conservation plays an important role in many physical Hamiltonians, we discuss its effect on typical pure states and on typical pure Gaussian states. We prove that while the behavior of the leading volume-law terms does not change qualitatively, the nature of the subleading terms can change. In particular, subleading corrections can appear that depend on the square root of the volume of the subsystem. We unveil the origin of those corrections. Finally, we discuss the connection between the entanglement entropy of typical pure states and analytical results obtained in the context of random matrix theory, as well as numerical results obtained for physical Hamiltonians.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Universal Landau-Zener regimes in dynamical topological phase transitions

In finite systems driven unitarily across topological phase transitions, the Chern number and the Bott index have been found to exhibit different behaviors depending on the boundary conditions and on the commensurability of the lattice. For periodic boundary conditions, the Chern number does not change for finite commensurate lattices (or in the thermodynamic limit). On the other hand, the Chern number can change for incommensurate lattices with periodic boundary conditions and the Bott index can change for lattices with open boundary conditions. Here we show that the scalings of the fields at which those two indices change exhibit Landau-Zener and near-adiabatic regimes depending on the speed at which the strength of the drive is ramped up and on the system size. Those regimes are preceded by a regime in which the topological indices do not change. The latter is the only regime that, for nonvanishing ramp speeds, survives in the thermodynamic limit. We then show that the dc Hall response can be used to detect topological phase transitions independently of the behavior of the topological indices.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Eigenstate Thermalization in a Locally Perturbed Integrable System

Eigenstate thermalization is widely accepted as the mechanism behind thermalization in generic isolated quantum systems. Using the example of a single magnetic defect embedded in the integrable spin-1/2 $XXZ$ chain, we show that locally perturbing an integrable system can give rise to eigenstate thermalization. Unique to such setups is the fact that thermodynamic and transport properties of the unperturbed integrable chain emerge in properties of the eigenstates of the perturbed (nonintegrable) one. Specifically, we show that the diagonal matrix elements of observables in the perturbed eigenstates follow the microcanonical predictions for the integrable model, and that the ballistic character of spin transport in the integrable model is manifest in the behavior of the off-diagonal matrix elements of the current operator in the perturbed eigenstates.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Low-frequency behavior of off-diagonal matrix elements in the integrable XXZ chain and in a locally perturbed quantum-chaotic XXZ chain

We study the matrix elements of local operators in the eigenstates of the integrable XXZ chain and of the quantum-chaotic model obtained by locally perturbing the XXZ chain with a magnetic impurity. We show that, at frequencies that are polynomially small in the system size, the behavior of the variances of the off-diagonal matrix elements can be starkly different depending on the operator. In the integrable model we find that, as the frequency $ω\rightarrow0$, the variances are either nonvanishing (generic behavior) or vanishing (for a special class of operators). In the quantum-chaotic model, on the other hand, we find the variances to be nonvanishing as $ω\rightarrow0$ and to indicate diffusive dynamics. We highlight which properties of the matrix elements of local operators are different between the integrable and quantum-chaotic models independently of the specific operator selected.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Entanglement and matrix elements of observables in interacting integrable systems

We study the bipartite von Neumann entanglement entropy and matrix elements of local operators in the eigenstates of an interacting integrable Hamiltonian (the paradigmatic spin-1/2 XXZ chain), and we contrast their behavior with that of quantum chaotic systems. We find that the leading term of the average (over all eigenstates in the zero magnetization sector) eigenstate entanglement entropy has a volume-law coefficient that is smaller than the universal (maximal entanglement) one in quantum chaotic systems. This establishes the entanglement entropy as a powerful measure to distinguish integrable models from generic ones. Remarkably, our numerical results suggest that the volume-law coefficient of the average entanglement entropy of eigenstates of the spin-1/2 XXZ Hamiltonian is very close to, or the same as, the one for translationally invariant quadratic fermionic models. We also study matrix elements of local operators in the eigenstates of the spin-1/2 XXZ Hamiltonian at the center of the spectrum. For the diagonal matrix elements, we show evidence that the support does not vanish with increasing system size, while the average eigenstate-to-eigenstate fluctuations vanish in a power-law fashion. For the off-diagonal matrix elements, we show that they follow a distribution that is close to (but not quite) log-normal, and that their variance is a well-defined function of $ω=E_α-E_β$ ($\{E_α\}$ are the eigenenergies) proportional to $1/D$, where $D$ is the Hilbert space dimension.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Observation of Dynamical Fermionization

We observe dynamical fermionization, where the momentum distribution of a Tonks-Girardeau (T-G) gas of strongly interacting bosons in 1D evolves from bosonic to fermionic after its axial confinement is removed. The asymptotic momentum distribution after expansion in 1D is the distribution of rapidities, which are the conserved quantities associated with many-body integrable systems. Rapidities have not previously been measured in any interacting many-body quantum system. Our measurements agree well with T-G gas theory. We also study momentum evolution after the trap depth is suddenly changed to a new non-zero value. We observe the predicted bosonic-fermionic oscillations and see deviations from the theory outside of the T-G gas limit.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dynamical preparation of Floquet Chern insulators

Realizing topological insulators is of great current interest because of their remarkable properties and possible future applications. There are recent proposals, based on Floquet analyses, that one can generate topologically nontrivial insulators by periodically driving topologically trivial ones. Here we address what happens if one follows the dynamics in such systems. Specifically, we present an exact study of the time evolution of a graphene-like system subjected to a circularly polarized electric field. We prove that, for infinite (translationally invariant) systems, the Chern number is conserved under unitary evolution. For systems with boundaries, on the other hand, we show that a properly defined topological invariant, the Bott index, can change. Hence, it should be possible to experimentally prepare topological states starting from non-topological ones. We show that the chirality of the edge current in such systems can be controlled by adjusting the filling.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Eigenstate thermalization in the two-dimensional transverse field Ising model

We study the onset of eigenstate thermalization in the two-dimensional transverse field Ising model (2D-TFIM) in the square lattice. We consider two non-equivalent Hamiltonians: the ferromagnetic 2D-TFIM and the antiferromagnetic 2D-TFIM in the presence of a uniform longitudinal field. We use full exact diagonalization to examine the behavior of quantum chaos indicators and of the diagonal matrix elements of operators of interest in the eigenstates of the Hamiltonian. A finite size scaling analysis reveals that quantum chaos and eigenstate thermalization occur in those systems whenever the fields are nonvanishing and not too large.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

From Quantum Chaos and Eigenstate Thermalization to Statistical Mechanics and Thermodynamics

This review gives a pedagogical introduction to the eigenstate thermalization hypothesis (ETH), its basis, and its implications to statistical mechanics and thermodynamics. In the first part, ETH is introduced as a natural extension of ideas from quantum chaos and random matrix theory (RMT). To this end, we present a brief overview of classical and quantum chaos, as well as RMT and some of its most important predictions. The latter include the statistics of energy levels, eigenstate components, and matrix elements of observables. Building on these, we introduce the ETH and show that it allows one to describe thermalization in isolated chaotic systems without invoking the notion of an external bath. We examine numerical evidence of eigenstate thermalization from studies of many-body lattice systems. We also introduce the concept of a quench as a means of taking isolated systems out of equilibrium, and discuss results of numerical experiments on quantum quenches. The second part of the review explores the implications of quantum chaos and ETH to thermodynamics. Basic thermodynamic relations are derived, including the second law of thermodynamics, the fundamental thermodynamic relation, fluctuation theorems, and the Einstein and Onsager relations. In particular, it is shown that quantum chaos allows one to prove these relations for individual Hamiltonian eigenstates and thus extend them to arbitrary stationary statistical ensembles. We then show how one can use these relations to obtain nontrivial universal energy distributions in continuously driven systems. At the end of the review, we briefly discuss the relaxation dynamics and description after relaxation of integrable quantum systems, for which ETH is violated. We introduce the concept of the generalized Gibbs ensemble, and discuss its connection with ideas of prethermalization in weakly interacting systems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Fundamental Asymmetry in Quenches Between Integrable and Nonintegrable Systems

We study quantum quenches between integrable and nonintegrable hard-core boson models in the thermodynamic limit with numerical linked cluster expansions. We show that while quenches in which the initial state is a thermal equilibrium state of an integrable model and the final Hamiltonian is nonintegrable (quantum chaotic) lead to thermalization, the reverse is not true. While this might appear counterintuitive given the fact that the eigenstates of both Hamiltonians are related by a unitary transformation, we argue that it is generic. Hence, the lack of thermalization of integrable systems is robust against quenches starting from stationary states of nonintegrable ones. Nonintegrable systems thermalize independently of the nature of the initial Hamiltonian.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Generalized Gibbs ensemble in integrable lattice models

The generalized Gibbs ensemble (GGE) was introduced ten years ago to describe observables in isolated integrable quantum systems after equilibration. Since then, the GGE has been demonstrated to be a powerful tool to predict the outcome of the relaxation dynamics of few-body observables in a variety of integrable models, a process we call generalized thermalization. This review discusses several fundamental aspects of the GGE and generalized thermalization in integrable systems. In particular, we focus on questions such as: which observables equilibrate to the GGE predictions and who should play the role of the bath; what conserved quantities can be used to construct the GGE; what are the differences between generalized thermalization in noninteracting systems and in interacting systems mappable to noninteracting ones; why is it that the GGE works when traditional ensembles of statistical mechanics fail. Despite a lot of interest in these questions in recent years, no definite answers have been given. We review results for the XX model and for the transverse field Ising model. For the latter model, we also report original results and show that the GGE describes spin-spin correlations over the entire system. This makes apparent that there is no need to trace out a part of the system in real space for equilibration to occur and for the GGE to apply. In the past, a spectral decomposition of the weights of various statistical ensembles revealed that generalized eigenstate thermalization occurs in the XX model (hard-core bosons). Namely, eigenstates of the Hamiltonian with similar distributions of conserved quantities have similar expectation values of few-spin observables. Here we show that generalized eigenstate thermalization also occurs in the transverse field Ising model.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Multiband Effects and the Bose-Hubbard Model in One-Dimensional Lattices

We study phase diagrams of one-dimensional bosons with contact interactions in the presence of a lattice. We use the worm algorithm in continuous space and focus on the incommensurate superfluid Mott-insulator transition. Our results are compared to those from the one-band Bose-Hubbard model. When Wannier states are used to determine the Bose-Hubbard model parameters, the comparison unveils an apparent breakdown of the one-band description for strong interactions, even for the Mott-insulating state with an average of one particle per site ($n=1$) in deep lattices. We introduce an inverse confined scattering analysis to obtain the ratio $U/J$, with which the Bose-Hubbard model provides correct results for strong interactions, deep lattices, and $n=1$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Observation of Spatial Charge and Spin Correlations in the 2D Fermi-Hubbard Model

Strong electron correlations lie at the origin of transformative phenomena such as colossal magneto-resistance and high-temperature superconductivity. Already near room temperature, doped copper oxide materials display remarkable features such as a pseudo-gap and a "strange metal" phase with unusual transport properties. The essence of this physics is believed to be captured by the Fermi-Hubbard model of repulsively interacting, itinerant fermions on a lattice. Here we report on the site-resolved observation of charge and spin correlations in the two-dimensional (2D) Fermi-Hubbard model realized with ultracold atoms. Antiferromagnetic spin correlations are maximal at half-filling and weaken monotonically upon doping. Correlations between singly charged sites are negative at large doping, revealing the Pauli and correlation hole\textemdash a suppressed probability of finding two fermions near each other. However, as the doping is reduced below a critical value, correlations between such local magnetic moments become positive, signaling strong bunching of doublons and holes. Excellent agreement with numerical linked-cluster expansion (NLCE) and determinantal quantum Monte Carlo (DQMC) calculations is found. Positive non-local moment correlations directly imply potential energy fluctuations due to doublon-hole pairs, which should play an important role for transport in the Fermi-Hubbard model.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Universal scaling of density and momentum distributions in Lieb-Liniger gases

We present an exact numerical study of the scaling of density and momentum distribution functions of harmonically trapped one-dimensional bosons with repulsive contact interactions at zero and finite temperatures. We use path integral quantum Monte Carlo with worm updates in our calculations at finite interaction strengths, and the Bose-Fermi mapping in the Tonks-Girardeau regime. We discuss the homogeneous case and, within the local density approximation, use it to motivate the scaling in the presence of a harmonic trap. For the momentum distribution function, we pay special attention to the high momentum tails and their $k^{-4}$ asymptotic behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Many-body localization and thermalization in disordered Hubbard chains

We study the many-body localization transition in one-dimensional Hubbard chains using exact diagonalization and quantum chaos indicators. We also study dynamics in the delocalized (ergodic) and localized phases and discuss thermalization and eigenstate thermalization, or the lack thereof, in such systems. Consistently within the indicators and observables studied, we find that ergodicity is very robust against disorder, namely, even in the presence of weak Hubbard interactions the disorder strength needed for the system to localize is large. We show that this robustness might be hidden by finite size effects in experiments with ultracold fermions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Mott-insulator-to-superconductor transition in a two-dimensional superlattice

We use quantum Monte Carlo and exact diagonalization calculations to study the Mott-insulator to superconductor quantum phase transition in a two-dimensional fermionic Hubbard model with attractive interactions in the presence of a superlattice potential. The model introduced offers unique possibilities to study such transitions in optical lattice experiments. We show that, in regimes with moderate to strong interactions, the transition belongs to the 3D-XY universality class. We also explore the character of the lowest energy charge excitations in the insulating and superconducting phases and show that they can be fermionic or bosonic depending on the parameters chosen.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Observation of coherent quench dynamics in a metallic many-body state of fermionic atoms

Quantum simulation with ultracold atoms has become a powerful technique to gain insight into interacting many-body systems. In particular, the possibility to study nonequilibrium dynamics offers a unique pathway to understand correlations and excitations in strongly interacting quantum matter. So far, coherent nonequilibrium dynamics has exclusively been observed in ultracold many-body systems of bosonic atoms. Here we report on the observation of coherent quench dynamics of fermionic atoms. A metallic state of ultracold spin-polarised fermions is prepared along with a Bose-Einstein condensate in a shallow three-dimensional optical lattice. After a quench that suppresses tunnelling between lattice sites for both the fermions and the bosons, we observe long-lived coherent oscillations in the fermionic momentum distribution, with a period that is determined solely by the Fermi-Bose interaction energy. Our results show that coherent quench dynamics can serve as a sensitive probe for correlations in delocalised fermionic quantum states and for quantum metrology.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Quantum distillation and confinement of vacancies in a doublon sea

Ultracold atomic gases have revolutionized the study of non-equilibrium dynamics in quantum many-body systems. Many counterintuitive non-equilibrium effects have been observed, such as suppressed thermalization in a one-dimensional (1D) gas, the formation of repulsive self-bound dimers, and identical behaviors for attractive and repulsive interactions. Here, we observe the expansion of a bundle of ultracold 1D Bose gases in a flat-bottomed optical lattice potential. By combining in situ measurements with photoassociation, we follow the spatial dynamics of singly, doubly, and triply occupied lattice sites. The system sheds interaction energy by dissolving some doublons and triplons. Some singlons quantum distill out of the doublon center, while others remain confined. Our Gutzwiller mean-field model captures these experimental features in a physically clear way. These experiments might be used to study thermalization in systems with particle losses or the evolution of quantum entanglement, or if applied to fermions, to prepare very low entropy states.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Quantum quenches and many-body localization in the thermodynamic limit

We use thermalization indicators and numerical linked cluster expansions to probe the onset of many-body localization in a disordered one-dimensional hard-core boson model in the thermodynamic limit. We show that after equilibration following a quench from a delocalized state, the momentum distribution indicates a freezing of one-particle correlations at higher values than in thermal equilibrium. The position of the delocalization to localization transition, identified by the breakdown of thermalization with increasing disorder strength, is found to be consistent with the value from the level statistics obtained via full exact diagonalization of finite chains. Our results strongly support the existence of a many-body localized phase in the thermodynamic limit.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Small quench dynamics as a probe for trapped ultracold atoms

Finite systems of bosons and/or fermions described by the Hubbard model can be realized using ultracold atoms confined in optical lattices. The ground states of these systems often exhibit a coexistence of compressible superfluid and incompressible Mott insulating regimes. We analyze such systems by studying the out-of-equilibrium dynamics following a weak sudden quench of the trapping potential. In particular, we show how the temporal variance of the site occupations reveals the location of spatial boundaries between compressible and incompressible regions. The feasibility of this approach is demonstrated for several models using numerical simulations. We first consider integrable systems, hard-core bosons (spinless fermions) confined by a harmonic potential, where space separated Mott and superfluid phases coexist. Then, we analyze a nonintegrable system, a $J-V-V'$ model with coexisting charge density wave and superfluid phases. We find that the temporal variance of the site occupations is a more effective measure than other standard indicators of phase boundaries such as a local compressibility. Based on these examples, we argue that analyzing temporal fluctuations is a valuable experimental tool for exploring phase boundaries in trapped atom systems.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Thermodynamics of two-dimensional spin models with bimodal random-bond disorder

We use numerical linked cluster expansions to study thermodynamic properties of the two-dimensional spin-1/2 Ising, XY, and Heisenberg models with bimodal random-bond disorder on the square and honeycomb lattices. In all cases, the nearest-neighbor coupling between the spins takes values $\pm J$ with equal probability. We obtain the disorder averaged (over all disorder configurations) energy, entropy, specific heat, and uniform magnetic susceptibility in each case. These results are compared with the corresponding ones in the clean models. Analytic expressions are obtained for low orders in the expansion of these thermodynamic quantities in inverse temperature.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Coherent quench dynamics in the one-dimensional Fermi-Hubbard model

Recently, it has been shown that the momentum distribution of a metallic state of fermionic atoms in a lattice Fermi-Bose mixture exhibits coherent oscillations after a global quench that suppresses tunneling. The oscillation period is determined by the Fermi-Bose interaction strength. Here we show that similar dynamics occurs in the fermionic Hubbard model when we quench a noninteracting metallic state by introducing a Hubbard interaction and suppressing tunneling. The period is determined primarily by the interaction strength. Conversely, we show that one can accurately determine the Hubbard interaction strength from the oscillation period, taking into account corrections from any small residual tunneling present in the final Hamiltonian. Such residual tunneling shortens the period and damps the oscillations, the latter being visible in the Fermi-Bose experiment.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Driven dipole oscillations and the lowest energy excitations of strongly interacting lattice bosons in a harmonic trap

We show that the analysis of the time evolution of the occupation of site and momentum modes of harmonically trapped lattice hard-core bosons, under driven dipole oscillations, allows one to determine the energy of the lowest one-particle excitations of the system in equilibrium. The analytic solution of a single particle in the absence of a lattice is used to identify which function of those time-dependent observables is best fit for the analysis, as well as to relate the dynamic response of the system to its single-particle spectrum. In the presence of the lattice and of multiple particles, a much richer and informative dynamical response is observed under the drive.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Linked-cluster expansion for the Green's function of the infinite-U Hubbard model

We implement a highly efficient strong-coupling expansion for the Green's function of the Hubbard model. In the limit of extreme correlations, where the onsite interaction is infinite, the evaluation of diagrams simplifies dramatically enabling us to carry out the expansion to the eighth order in powers of the hopping amplitude. We compute the finite-temperature Green's function analytically in the momentum and Matsubara frequency space as a function of the electron density. Employing Pade approximations, we study the equation of state, Kelvin thermopower, momentum distribution function, quasiparticle fraction, and quasiparticle lifetime of the system at temperatures lower than, or of the order of, the hopping amplitude. We also discuss several different approaches for obtaining the spectral functions through analytic continuation of the imaginary frequency Green's function, and show results for the system near half filling. We benchmark our results for the equation of state against those obtained from a numerical linked-cluster expansion carried out to the eleventh order.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Nature of the phases in the frustrated XY model on the honeycomb lattice

We study the phase diagram of the frustrated XY model on the honeycomb lattice by using accurate correlated wave functions and variational Monte Carlo simulations. Our results suggest that a spin-liquid state is energetically favorable in the region of intermediate frustration, intervening between two magnetically ordered phases. The latter ones are represented by classically ordered states supplemented with a long-range Jastrow factor, which includes relevant correlations and dramatically improves the description provided by the purely classical solution of the model. The construction of the spin-liquid state is based on a decomposition of the underlying bosonic particles in terms of spin-1/2 fermions (partons), with a Gutzwiller projection enforcing no single occupancy, as well as a long-range Jastrow factor.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Nonequilibrium dynamics of one-dimensional hard-core anyons following a quench: Complete relaxation of one-body observables

We demonstrate the role of interactions in driving the relaxation of an isolated integrable quantum system following a sudden quench. We consider a family of integrable hard-core lattice anyon models that continuously interpolates between noninteracting spinless fermions and strongly interacting hard-core bosons. A generalized Jordan-Wigner transformation maps the entire family to noninteracting fermions. We find that, aside from the singular free-fermion limit, the entire single-particle density matrix and therefore all one-body observables relax to the predictions of the generalized Gibbs ensemble (GGE). This demonstrates that, in the presence of interactions, correlations between particles in the many-body wave function provide the effective dissipation required to drive relaxation of all one-body observables to the GGE. This relaxation does not depend on translational invariance, or the tracing out of any spatial domain of the system.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Quantum Quenches in the Thermodynamic Limit

We introduce a linked-cluster based computational approach that allows one to study quantum quenches in lattice systems in the thermodynamic limit. This approach is used to study quenches in one-dimensional lattices. We provide evidence that, in the thermodynamic limit, thermalization occurs in the nonintegrable regime but fails at integrability. A phase transitionlike behavior separates the two regimes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Quantum quenches in the thermodynamic limit. II. Initial ground states

A numerical linked-cluster algorithm was recently introduced to study quantum quenches in the thermodynamic limit starting from thermal initial states [M. Rigol, Phys. Rev. Lett. 112, 170601 (2014)]. Here, we tailor that algorithm to quenches starting from ground states. In particular, we study quenches from the ground state of the antiferromagnetic Ising model to the XXZ chain. Our results for spin correlations are shown to be in excellent agreement with recent analytical calculations based on the quench action method. We also show that they are different from the correlations in thermal equilibrium, which confirms the expectation that thermalization does not occur in general in integrable models even if they cannot be mapped to noninteracting ones.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Quenching the Anisotropic Heisenberg Chain: Exact Solution and Generalized Gibbs Ensemble Predictions

We study quenches in integrable spin-1/2 chains in which we evolve the ground state of the antiferromagnetic Ising model with the anisotropic Heisenberg Hamiltonian. For this nontrivially interacting situation, an application of the first-principles-based quench action method allows us to give an exact description of the postquench steady state in the thermodynamic limit. We show that a generalized Gibbs ensemble, implemented using all known local conserved charges, fails to reproduce the exact quench action steady state and to correctly predict postquench equilibrium expectation values of physical observables. This is supported by numerical linked-cluster calculations within the diagonal ensemble in the thermodynamic limit.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Self trapping in the two-dimensional Bose-Hubbard model

We study the expansion of harmonically trapped bosons in a two-dimensional lattice after suddenly turning off the confining potential. We show that, in the presence of multiple occupancies per lattice site and strong interactions, the system exhibits a clear dynamical separation into slowly and rapidly expanding clouds. We discuss how this effect can be understood within a simple picture by invoking doublons and Bose enhancement. This picture is corroborated by an analysis of the momentum distribution function in the regions with slowly and rapidly expanding bosons.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Spiral antiferromagnets beyond the spin-wave approximation: frustrated $XY$ and Heisenberg models in the honeycomb lattice

We examine the stability of classical states with a generic incommensurate spiral order against quantum fluctuations. Specifically, we focus on the frustrated spin-1/2 $XY$ and Heisenberg models on the honeycomb lattice with nearest-neighbor $J_1$ and next-nearest-neighbor $J_2$ antiferromagnetic couplings. Our variational approach is based on the Jastrow wave functions, which include quantum correlations on top of classical spin waves. We perform a systematic optimization of wave vectors and Jastrow pseudo-potentials within this class of variational states and find that quantum fluctuations favor collinear states over generic coplanar spirals. The Néel state with ${\bf Q}=(0,0)$ extends its stability well beyond the classical value $J_2/J_1=1/6$. Most importantly, the collinear states with ${\bf Q}=(0,2π/\sqrt{3})$ (and the two symmetry-related states) are found to be stable in a large regime with intermediate frustration, while at the classical level they are limited to the point $J_2/J_1=0.5$. For large frustration, the $120^\circ$ state is stabilized for finite values of $J_2/J_1$ in both models.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Dynamics and thermalization in correlated one-dimensional lattice systems

We review exact approaches and recent results related to the relaxation dynamics and description after relaxation of various one-dimensional lattice systems of hard-core bosons after a sudden quench. We first analyze the integrable case, where the combination of analytical insights and computational techniques enable one to study large system sizes. Thermalization does not occur in this regime. However, after relaxation, observables can be described by a generalization of the Gibbs ensemble. We then utilize full exact diagonalization to study what happens as integrability is broken. We show that thermalization does occur in finite nonintegrable systems provided they are sufficiently far away from the integrable point. We argue that the onset of thermalization can be understood in terms of the eigenstate thermalization hypothesis.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Electronic spectral properties of the two-dimensional infinite-U Hubbard model

A strong-coupling series expansion for the Green's function and the extremely-correlated Fermi liquid (ECFL) theory are used to calculate the moments of the electronic spectral functions of the infinite-U Hubbard model. Results from these two complementary methods agree very well at both, low densities, where the ECFL solution is the most accurate, and at high to intermediate temperatures, where the series converge. We find that a modified first moment, which underestimates the contributions from the occupied states and is accessible in the series through the time-dependent Green's function, best describes the quasiparticle peak location in the strongly-correlated regime. This is examined by the ECFL results at low temperatures, where it is shown that the spectral function is largely skewed towards the occupied states.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Finite-temperature properties of strongly correlated fermions in the honeycomb lattice

We study finite-temperature properties of strongly interacting fermions in the honeycomb lattice using numerical linked-cluster expansions and determinantal quantum Monte Carlo simulations. We analyze a number of thermodynamic quantities, including the entropy, the specific heat, uniform and staggered spin susceptibilities, short-range spin correlations, and the double occupancy at and away from half filling. We examine the viability of adiabatic cooling by increasing the interaction strength for homogeneous as well as for trapped systems. For the homogeneous case, this process is found to be more efficient at finite doping than at half filling. That, in turn, leads to an efficient adiabatic cooling in the presence of a trap which, starting with even relatively high entropies, can drive the system to have a Mott insulating phase with substantial antiferromagnetic correlations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Fluctuation-Dissipation Theorem in an Isolated System of Quantum Dipolar Bosons after a Quench

We examine the validity of fluctuation-dissipation relations in isolated quantum systems taken out of equilibrium by a sudden quench. We focus on the dynamics of trapped hard-core bosons in one-dimensional lattices with dipolar interactions whose strength is changed during the quench. We find that fluctuation-dissipation relations hold if the system is nonintegrable after the quench. They also hold if the system is integrable after the quench if the initial state is an equilibrium state of a nonintegrable Hamiltonian. However, they fail if the system is integrable both before and after quenching.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Initial state dependence of the quench dynamics in integrable quantum systems. III. Chaotic states

We study sudden quantum quenches in which the initial states are selected to be either eigenstates of an integrable Hamiltonian that is nonmappable to a noninteracting one or a nonintegrable Hamiltonian, while the Hamiltonian after the quench is always integrable and mappable to a noninteracting one. By studying weighted energy densities and entropies, we show that quenches starting from nonintegrable (chaotic) eigenstates lead to an "ergodic" sampling of the eigenstates of the final Hamiltonian, while those starting from the integrable eigenstates do not (or at least it is not apparent for the system sizes accessible to us). This goes in parallel with the fact that the distribution of conserved quantities in the initial states is thermal in the nonintegrable cases and nonthermal in the integrable ones, and means that, in general, thermalization occurs in integrable systems when the quench starts form an eigenstate of a nonintegrable Hamiltonian (away from the edges of the spectrum), while it fails (or requires larger system sizes than those studied here to become apparent) for quenches starting at integrable points. We test those conclusions by studying the momentum distribution function of hard-core bosons after a quench.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Localization and the effects of symmetries in the thermalization properties of one-dimensional quantum systems

We study how the proximity to an integrable point or to localization as one approaches the atomic limit, as well as the mixing of symmetries in the chaotic domain, may affect the onset of thermalization in finite one-dimensional systems. We consider systems of hard-core bosons at half-filling with nearest neighbor hopping and interaction, and next-nearest neighbor interaction. The latter breaks integrability and induces a ground-state superfluid to insulator transition. By full exact diagonalization, we study chaos indicators and few-body observables. We show that when different symmetry sectors are mixed, chaos indicators associated with the eigenvectors, contrary to those related to the eigenvalues, capture the onset of chaos. The results for the complexity of the eigenvectors and for the expectation values of few-body observables confirm the validity of the eigenstate thermalization hypothesis in the chaotic regime, and therefore the occurrence of thermalization. We also study the properties of the off-diagonal matrix elements of few-body observables in relation to the transition from integrability to chaos and from chaos to localization.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Quantum quenches and thermalization in one-dimensional fermionic systems

We study the dynamics and thermalization of strongly correlated fermions in finite one-dimensional lattices after a quantum quench. Our calculations are performed using exact diagonalization. We focus on one- and two-body observables such as the momentum distribution function [n(k)] and the density-density structure factor [N(k)], respectively, and study the effects of approaching an integrable point. We show that while the relaxation dynamics and thermalization of N(k) for fermions is very similar to the one of hardcore bosons, the behavior of n(k) is distinctively different. The latter observable exhibits a slower relaxation dynamics in fermionic systems. We identify the origin of this behavior, which is related to the off-diagonal matrix elements of n(k) in the basis of the eigenstates of the Hamiltonian. More generally, we find that thermalization occurs far away from integrability and that it breaks down as one approaches the integrable point.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Scaling of the gap, fidelity susceptibility, and Bloch oscillations across the superfluid to Mott insulator transition in the one-dimensional Bose-Hubbard model

We investigate the interaction-induced superfluid-to-Mott insulator transition in the one-dimensional Bose-Hubbard model (BHM) for fillings $n=1$, $n=2$, and $n=3$ by studying the single-particle gap, the fidelity susceptibility, and the amplitude of Bloch oscillations via density-matrix renormalization-group methods. We apply a generic scaling procedure for the gap, which allows us to determine the critical points with very high accuracy. We also study how the fidelity susceptibility behaves across the phase transition. Furthermore, we show that in the BHM, and in a system of spinless fermions, the amplitude of Bloch oscillations after a tilt of the lattice vanishes at the critical points. This indicates that Bloch oscillations can serve as a tool to detect the transition point in ongoing experiments with ultracold gases.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Single-particle and many-body analyses of a quasiperiodic integrable system after a quench

In general, isolated integrable quantum systems have been found to relax to an apparent equilibrium state in which the expectation values of few-body observables are described by the generalized Gibbs ensemble. However, recent work has shown that relaxation to such a generalized statistical ensemble can be precluded by localization in a quasiperiodic lattice system. Here we undertake complementary single-particle and many-body analyses of noninteracting spinless fermions and hard-core bosons within the Aubry-Andre model to gain insight into this phenomenon. Our investigations span both the localized and delocalized regimes of the quasiperiodic system, as well as the critical point separating the two. Considering first the case of spinless fermions, we study the dynamics of the momentum distribution function and characterize the effects of real-space and momentum-space localization on the relevant single-particle wave functions and correlation functions. We show that although some observables do not relax in the delocalized and localized regimes, the observables that do relax in these regimes do so in a manner consistent with a recently proposed Gaussian equilibration scenario, whereas relaxation at the critical point has a more exotic character. We also construct various statistical ensembles from the many-body eigenstates of the fermionic and bosonic Hamiltonians and study the effect of localization on their properties.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Short Introduction to Numerical Linked-Cluster Expansions

We provide a pedagogical introduction to numerical linked-cluster expansions (NLCEs). We sketch the algorithm for generic Hamiltonians that only connect nearest-neighbor sites in a finite cluster with open boundary conditions. We then compare results for a specific model, the Heisenberg model, in each order of the NLCE with the ones for the finite cluster calculated directly by means of full exact diagonalization. We discuss how to reduce the computational cost of the NLCE calculations by taking into account symmetries and topologies of the linked clusters. Finally, we generalize the algorithm to the thermodynamic limit, and discuss several numerical resummation techniques that can be used to accelerate the convergence of the series.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Alternatives to Eigenstate Thermalization

An isolated quantum many-body system in an initial pure state will come to thermal equilibrium if it satisfies the eigenstate thermalization hypothesis (ETH). We consider alternatives to ETH that have been proposed. We first show that von Neumann's quantum ergodic theorem relies on an assumption that is essentially equivalent to ETH. We also investigate whether, following a sudden quench, special classes of pure states can lead to thermal behavior in systems that do not obey ETH, namely, integrable systems. We find examples of this, but only for initial states that obeyed ETH before the quench.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

An Exactly Solvable Model for the Integrability-Chaos Transition in Rough Quantum Billiards

A central question of dynamics, largely open in the quantum case, is to what extent it erases a system's memory of its initial properties. Here we present a simple statistically solvable quantum model describing this memory loss across an integrability-chaos transition under a perturbation obeying no selection rules. From the perspective of quantum localization-delocalization on the lattice of quantum numbers, we are dealing with a situation where every lattice site is coupled to every other site with the same strength, on average. The model also rigorously justifies a similar set of relationships recently proposed in the context of two short-range-interacting ultracold atoms in a harmonic waveguide. Application of our model to an ensemble of uncorrelated impurities on a rectangular lattice gives good agreement with ab initio numerics.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Effect of particle statistics in strongly correlated two-dimensional Hubbard models

We study the onset of particle statistics effects as the temperature is lowered in strongly correlated two-dimensional Hubbard models. We utilize numerical linked-cluster expansions and focus on the properties of interacting lattice fermions and two-component hard-core bosons. In the weak-coupling regime, where the ground state of the bosonic system is a superfluid, the thermodynamic properties of the two systems at half filling exhibit very large differences even at high temperatures. In the strong-coupling regime, where the low-temperature behavior is governed by a Mott insulator for either particle statistics, the agreement between the thermodynamic properties of both systems extends to regions where the antiferromagnetic (iso)spin correlations are exponentially large. We analyze how particle statistics affects adiabatic cooling in those systems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Entropy of isolated quantum systems after a quench

A diagonal entropy, which depends only on the diagonal elements of the system's density matrix in the energy representation, has been recently introduced as the proper definition of thermodynamic entropy in out-of-equilibrium quantum systems. We study this quantity after an interaction quench in lattice hard-core bosons and spinless fermions, and after a local chemical potential quench in a system of hard-core bosons in a superlattice potential. The former systems have a chaotic regime, where the diagonal entropy becomes equivalent to the equilibrium microcanonical entropy, coinciding with the onset of thermalization. The latter system is integrable. We show that its diagonal entropy is additive and different from the entropy of a generalized Gibbs ensemble, which has been introduced to account for the effects of conserved quantities at integrability.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Expansion of Bose-Hubbard Mott insulators in optical lattices

We study the expansion of bosonic Mott insulators in the presence of an optical lattice after switching off a confining potential. We use the Gutzwiller mean-field approximation and consider two different setups. In the first one, the expansion is restricted to one direction. We show that this leads to the emergence of two condensates with well defined momenta, and argue that such a construct can be used to create atom lasers in optical lattices. In the second setup, we study Mott insulators that are allowed to expand in all directions in the lattice. In this case, a simple condensate is seen to develop within the mean-field approximation. However, its constituent bosons are found to populate many nonzero momentum modes. An analytic understanding of both phenomena in terms of the exact dispersion relation in the hard-core limit is presented.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Initial state dependence of the quench dynamics in integrable quantum systems. II. Thermal states

We study properties of isolated integrable quantum systems after a sudden quench starting from thermal states. We show that, even if the system is initially in thermal equilibrium at finite temperature, the diagonal entropy after a quench remains a fraction of the entropy in the generalized ensembles introduced to describe integrable systems after relaxation. The latter is also, in general, different from the entropy in thermal equilibrium. Furthermore, we examine the difference between the distribution of conserved quantities in the thermal and generalized ensembles after a quench and show that they are also, in general, different from each other. This explains why these systems fail to thermalize. A finite size scaling analysis is presented for each quantity, which allows us making predictions for thermodynamically large lattice sizes.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Noise Correlation Scalings: Revisiting the Quantum Phase Transition in Incommensurate Lattices with Hard-Core Bosons

Finite size scalings of the momentum distribution and noise correlations are performed to study Mott insulator, Bose glass, and superfluid quantum phases in hard-core bosons (HCBs) subjected to quasi-periodic disorder. The exponents of the correlation functions at the Superfluid to Bose glass (SF-BG) transition are found to be approximately one half of the ones that characterizes the superfluid phase. The derivatives of the peak intensities of the correlation functions with respect to quasiperiodic disorder are shown to diverge at the SF-BG critical point. This behavior does not occur in the corresponding free fermion system, which also exhibits an Anderson-like transition at the same critical point, and thus provides a unique experimental tool to locate the phase transition in interacting bosonic systems. We also report on the absence of primary sublattice peaks in the momentum distribution of the superfluid phase for special fillings.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Numerical study of the thermodynamics of clinoatacamite

We study the thermodynamic properties of the clinoatacamite compound, Cu_2(OH)_3Cl, by considering several approximate models. They include the Heisenberg model on (i) the uniform pyrochlore lattice, (ii) a very anisotropic pyrochlore lattice, and (iii) a kagome lattice weakly coupled to spins that sit on a triangular lattice. We utilize the exact diagonalization of small clusters with periodic boundary conditions and implement a numerical linked-cluster expansion approach for quantum lattice models with reduced symmetries, which allows us to solve model (iii) in the thermodynamic limit. We find a very good agreement between the experimental uniform susceptibility and the numerical results for models (ii) and (iii), which suggests a weak ferromagnetic coupling between the kagome and triangular layers in clinoatacamite. We also study thermodynamic properties in a geometrical transition between a planar pyrochlore lattice and the kagome lattice.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Quantum quenches in disordered systems: Approach to thermal equilibrium without a typical relaxation time

We study spectral properties and the dynamics after a quench of one-dimensional spinless fermions with short-range interactions and long-range random hopping. We show that a sufficiently fast decay of the hopping term promotes localization effects at finite temperature, which prevents thermalization even if the classical motion is chaotic. For slower decays, we find that thermalization does occur. However, within this model, the latter regime falls in an unexpected universality class, namely, observables exhibit a power-law (as opposed to an exponential) approach to their thermal expectation values.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Quenches in a quasi-disordered integrable lattice system: Dynamics and statistical description of observables after relaxation

We study the dynamics and the resulting state after relaxation in a quasi-disordered integrable lattice system after a sudden quench. Specifically, we consider hard-core bosons in an isolated one-dimensional geometry in the presence of a quasi-periodic potential whose strength is abruptly changed to take the system out of equilibrium. In the delocalized regime, we find that the relaxation dynamics of one-body observables, such as the density, the momentum distribution function, and the occupation of the natural orbitals, follow, to a good approximation, power laws. In that regime, we also show that the observables after relaxation can be described by the generalized Gibbs ensemble, while such a description fails for the momentum distribution and the natural orbital occupations in the presence of localization. At the critical point, the relaxation dynamics is found to be slower than in the delocalized phase.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Short-Range Correlations and Cooling of Ultracold Fermions in the Honeycomb Lattice

We use determinantal quantum Monte Carlo simulations and numerical linked-cluster expansions to study thermodynamic properties and short-range spin correlations of fermions in the honeycomb lattice. We find that, at half filling and finite temperatures, nearest-neighbor spin correlations can be stronger in this lattice than in the square lattice, even in regimes where the ground state in the former is a semimetal or a spin liquid. The honeycomb lattice also exhibits a more pronounced anomalous region in the double occupancy that leads to stronger adiabatic cooling than in the square lattice. We discuss the implications of these findings for optical lattice experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Superfluid to normal phase transition in strongly correlated bosons in two and three dimensions

Using quantum Monte Carlo simulations, we investigate the finite-temperature phase diagram of hard-core bosons (XY model) in two- and three-dimensional lattices. To determine the phase boundaries, we perform a finite-size-scaling analysis of the condensate fraction and/or the superfluid stiffness. We then discuss how these phase diagrams can be measured in experiments with trapped ultracold gases, where the systems are inhomogeneous. For that, we introduce a method based on the measurement of the zero-momentum occupation, which is adequate for experiments dealing with both homogeneous and trapped systems, and compare it with previously proposed approaches.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Weak and strong typicality in quantum systems

We study the properties of mixed states obtained from eigenstates of many-body lattice Hamiltonians after tracing out part of the lattice. Two scenarios emerge for generic systems: (i) the diagonal entropy becomes equivalent to the thermodynamic entropy when a few sites are traced out (weak typicality); and (ii) the von Neumann (entanglement) entropy becomes equivalent to the thermodynamic entropy when a large fraction of the lattice is traced out (strong typicality). Remarkably, the results for few-body observables obtained with the reduced, diagonal, and canonical density matrices are very similar to each other, no matter which fraction of the lattice is traced out. Hence, for all physical quantities studied here, the results in the diagonal ensemble match the thermal predictions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Analytical and numerical study of trapped strongly correlated bosons in two- and three-dimensional lattices

We study the ground-state properties of trapped inhomogeneous systems of hardcore bosons in two- and three-dimensional lattices. We obtain our results both numerically, using quantum Monte Carlo techniques, and via several analytical approximation schemes, such as the Gutzwiller-mean-field approach, a cluster-mean-field method and a spin-wave analysis which takes quantum fluctuations into account. We first study the homogeneous case, for which simple analytical expressions are obtained for all observables of interest, and compare the results with the numerical ones. We obtain the equation of state of the system along with other thermodynamic properties such as the free energy, kinetic energy, superfluid density, condensate fraction and compressibility. In the presence of a trap, superfluid and insulating domains coexist in the system. We show that the spin-wave-based method reproduces the quantum Monte-Carlo results for global as well as for local quantities with a high degree of accuracy. We also discuss the validity of the local density approximation in those systems. Our analysis can be used to describe bosons in optical lattices where the onsite interaction U is much larger than the hopping amplitude t.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Collapse and revival oscillations as a probe for the tunneling amplitude in an ultra-cold Bose gas

We present a theoretical study of the quantum corrections to the revival time due to finite tunneling in the collapse and revival of matter wave interference after a quantum quench. We study hard-core bosons in a superlattice potential and the Bose-Hubbard model by means of exact numerical approaches and mean-field theory. We consider systems without and with a trapping potential present. We show that the quantum corrections to the revival time can be used to accurately determine the value of the hopping parameter in experiments with ultracold bosons in optical lattices.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Generalized Thermalization in an Integrable Lattice System

After a quench, observables in an integrable system may not relax to the standard thermal values, but can relax to the ones predicted by the generalized Gibbs ensemble (GGE) [M. Rigol et al., Phys. Rev. Lett. 98, 050405 (2007)]. The GGE has been shown to accurately describe observables in various one-dimensional integrable systems, but the origin of its success is not fully understood. Here we introduce a microcanonical version of the GGE and provide a justification of the GGE based on a generalized interpretation of the eigenstate thermalization hypothesis, which was previously introduced to explain thermalization of nonintegrable systems. We study relaxation after a quench of one-dimensional hard-core bosons in an optical lattice. Exact numerical calculations for up to 10 particles on 50 lattice sites (~10^10 eigenstates) validate our approach.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Initial state dependence of the quench dynamics in integrable quantum systems

We identify and study classes of initial states in integrable quantum systems that, after the relaxation dynamics following a sudden quench, lead to near-thermal expectation values of few-body observables. In the systems considered here, those states are found to be insulating ground states of lattice hard-core boson Hamiltonians. We show that, as a suitable parameter in the initial Hamiltonian is changed, those states become closer to Fock states (products of single site states) as the outcome of the relaxation dynamics becomes closer to the thermal prediction. At the same time, the energy density approaches a Gaussian. Furthermore, the entropy associated with the generalized canonical and generalized grand-canonical ensembles, introduced to describe observables in integrable systems after relaxation, approaches that of the conventional canonical and grand-canonical ensembles. We argue that those classes of initial states are special because a control parameter allows one to tune the distribution of conserved quantities to approach the one in thermal equilibrium. This helps in understanding the approach of all the quantities studied to their thermal expectation values. However, a finite-size scaling analysis shows that this behavior should not be confused with thermalization as understood for nonintegrable systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Kaleidoscope of exotic quantum phases in a frustrated XY model

The existence of quantum spin liquids was first conjectured by Pomeranchuk some 70 years ago, who argued that frustration in simple antiferromagnetic theories could result in a Fermi-liquid-like state for spinon excitations. Here we show that a simple quantum spin model on a honeycomb lattice hosts the long sought for Bose metal with a clearly identifiable Bose surface. The complete phase diagram of the model is determined via exact diagonalization and is shown to include four distinct phases separated by three quantum phase transitions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Magnetism and pairing of two-dimensional trapped fermions

The emergence of local phases in a trapped two-component Fermi gas in an optical lattice is studied using quantum Monte Carlo simulations. We treat temperatures that are comparable or lower than those presently achievable in experiments and large enough systems that both magnetic and paired phases can be detected by inspection of the behavior of suitable short-range correlations. We use the latter to suggest the interaction strength and temperature range at which experimental observation of incipient magnetism and d-wave pairing are more likely and evaluate the relation between entropy and temperature in two-dimensional confined fermionic systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Phase diagram of the half-filled one-dimensional t-V-V' model

We study the phase diagram of spinless fermions with nearest and next-nearest-neighbor interactions in one dimension utilizing the (finite-size) density-matrix renormalization group (DMRG) method. The competition between nearest and next-nearest-neighbor interactions and nearest-neighbor hopping generates four phases in this model: two charge-density-wave insulators, a Luttinger liquid phase, and a bond-order phase. We use finite-size scaling of the gap and various structure factors to determine the phase diagram.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Scaling of noise correlations in one-dimensional-lattice-hard-core-boson systems

Noise correlations are studied for systems of hard-core bosons in one-dimensional lattices. We use an exact numerical approach based on the Bose-Fermi mapping and properties of Slater determinants. We focus on the scaling of the noise correlations with system size in superfluid and insulating phases, which are generated in the homogeneous lattice, with period-two superlattices and with uniformly distributed random diagonal disorder. For the superfluid phases, the leading contribution is shown to exhibit a density-independent scaling proportional to the system size, while the first subleading term exhibits a density-dependent power-law exponent.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Thermodynamics and phase transitions for the Heisenberg model on the pinwheel distorted kagome lattice

We study the Heisenberg model on the pinwheel distorted kagome lattice as observed in the material Rb_2Cu_3SnF_12. Experimentally relevant thermodynamic properties at finite temperatures are computed utilizing numerical linked-cluster expansions. We also develop a Lanczos-based, zero-temperature, numerical linked cluster expansion to study the approach of the pinwheel distorted lattice to the uniform kagome-lattice Heisenberg model. We find strong evidence for a phase transition before the uniform limit is reached, implying that the ground state of the kagome-lattice Heisenberg model is likely not pinwheel dimerized and is stable to finite pinwheel-dimerizing perturbations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Thermodynamics of strongly interacting fermions in two-dimensional optical lattices

We study finite-temperature properties of strongly correlated fermions in two-dimensional optical lattices by means of numerical linked cluster expansions, a computational technique that allows one to obtain exact results in the thermodynamic limit. We focus our analysis on the strongly interacting regime, where the on-site repulsion is of the order of or greater than the band width. We compute the equation of state, double occupancy, entropy, uniform susceptibility, and spin correlations for temperatures that are similar to or below the ones achieved in current optical lattice experiments. We provide a quantitative analysis of adiabatic cooling of trapped fermions in two dimensions, by means of both flattening the trapping potential and increasing the interaction strength.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Thermodynamics of the Antiferromagnetic Heisenberg Model on the Checkerboard Lattice

Employing numerical linked-cluster expansions (NLCEs) along with exact diagonalizations of finite clusters with periodic boundary condition, we study the energy, specific heat, entropy, and various susceptibilities of the antiferromagnetic Heisenberg model on the checkerboard lattice. NLCEs, combined with extrapolation techniques, allow us to access temperatures much lower than those accessible to exact diagonalization and other series expansions. We find that the high-temperature peak in specific heat decreases as the frustration increases, consistent with the large amount of unquenched entropy in the region around maximum classical frustration, where the nearest-neighbor and next-nearest neighbor exchange interactions (J and J', respectively) have the same strength, and with the formation of a second peak at lower temperatures. The staggered susceptibility shows a change of character when J' increases beyond 0.75J, implying the disappearance of the long-range antiferromagnetic order at zero temperature. For J'=4J, in the limit of weakly coupled crossed chains, we find large susceptibilities for stripe and Neel order with Q=(pi/2,pi/2) at low temperatures with antiferromagnetic correlations along the chains. Other magnetic and bond orderings, such as a plaquette valence-bond solid and a crossed-dimer order suggested by previous studies, have also been investigated.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Topological Phase Transitions for Interacting Finite Systems

In this paper, we investigate signatures of topological phase transitions in interacting systems. We show that the key signature is the existence of a topologically protected level crossing, which is robust and sharply defines the topological transition, even in finite-size systems. Spatial symmetries are argued to play a fundamental role in the selection of the boundary conditions to be used to locate topological transitions in finite systems. We discuss the theoretical implications of this result, and utilize exact diagonalization to demonstrate its manifestations in the Haldane-Fermi-Hubbard model. Our findings provide an efficient way to detect topological transitions in experiments and in numerical calculations that cannot access the ground-state wave function.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Comment on "Quenches in quantum many-body systems: One-dimensional Bose-Hubbard model reexamined" [arXiv:0810.3720]

In a recent paper Roux [Phys. Rev. A 79, 021608(R) (2009), arXiv:0810.3720] argued that thermalization in a Bose-Hubbard system, after a quench, follows from the approximate Boltzmann distribution of the overlap between the initial state and the eigenstates of the final Hamiltonian. We show here that the distribution of the overlaps is in general not related to the canonical (or microcanonical) distribution and, hence, it cannot explain why thermalization occurs in quantum systems.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Interaction effects and quantum phase transitions in topological insulators

We study strong correlation effects in topological insulators via the Lanczos algorithm, which we utilize to calculate the exact many-particle ground-state wave function and its topological properties. We analyze the simple, noninteracting Haldane model on a honeycomb lattice with known topological properties and demonstrate that these properties are already evident in small clusters. Next, we consider interacting fermions by introducing repulsive nearest-neighbor interactions. A first-order quantum phase transition was discovered at finite interaction strength between the topological band insulator and a topologically trivial Mott insulating phase by use of the fidelity metric and the charge-density-wave structure factor. We construct the phase diagram at $T = 0$ as a function of the interaction strength and the complex phase for the next-nearest-neighbor hoppings. Finally, we consider the Haldane model with interacting hard-core bosons, where no evidence for a topological phase is observed. An important general conclusion of our work is that despite the intrinsic nonlocality of topological phases their key topological properties manifest themselves already in small systems and therefore can be studied numerically via exact diagonalization and observed experimentally, e.g., with trapped ions and cold atoms in optical lattices.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Onset of quantum chaos in one-dimensional bosonic and fermionic systems and its relation to thermalization

By means of full exact diagonalization, we study level statistics and the structure of the eigenvectors of one-dimensional gapless bosonic and fermionic systems across the transition from integrability to quantum chaos. These systems are integrable in the presence of only nearest-neighbor terms, whereas the addition of next-nearest neighbor hopping and interaction may lead to the onset of chaos. We show that the strength of the next-nearest neighbor terms required to observe clear signatures of nonintegrability is inversely proportional to the system size. Interestingly, the transition to chaos is also seen to depend on particle statistics, with bosons responding first to the integrability breaking terms. In addition, we discuss the use of delocalization measures as main indicators for the crossover from integrability to chaos and the consequent viability of quantum thermalization in isolated systems.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Phase diagram of the hardcore Bose-Hubbard model on a checkerboard superlattice

We obtain the complete phase diagram of the hardcore Bose-Hubbard model in the presence of a period-two superlattice in two and three dimensions. First we acquire the phase boundaries between the superfluid phase and the `trivial' insulating phases of the model (the completely-empty and completely-filled lattices) analytically. Next, the boundary between the superfluid phase and the half-filled Mott-insulating phase is obtained numerically, using the stochastic series expansion (SSE) algorithm followed by finite-size scaling. We also compare our numerical results against the predictions of several approximation schemes, including two mean-field approaches and a fourth-order strong-coupling expansion (SCE), where we show that the latter method in particular is successful in producing an accurate picture of the phase diagram. Finally, we examine the extent to which several approximation schemes, such as the random phase approximation and the strong-coupling expansion, give an accurate description of the momentum distribution of the bosons inside the insulating phases.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Quantum chaos and thermalization in gapped systems

We investigate the onset of thermalization and quantum chaos in finite one-dimensional gapped systems of hard-core bosons. Integrability in these systems is broken by next-nearest-neighbor repulsive interactions, which also generate a superfluid to insulator transition. By employing full exact diagonalization, we study chaos indicators and few-body observables. We show that with increasing system size, chaotic behavior is seen over a broader range of parameters and, in particular, deeper into the insulating phase. Concomitantly, we observe that, as the system size increases, the eigenstate thermalization hypothesis extends its range of validity inside the insulating phase and is accompanied by the thermalization of the system.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Strongly Interacting Atom Lasers in Three Dimensional Optical Lattices

We show that the dynamical melting of a Mott insulator in a three-dimensional lattice leads to condensation at nonzero momenta, a phenomenon that can be used to generate strongly interacting atom lasers in optical lattices. For infinite onsite repulsion, the case considered here, the momenta at which bosons condense is determined analytically and found to have a simple dependence on the hopping amplitudes. The occupation of the condensates is shown to scale linearly with the total number of atoms in the initial Mott insulator. Our results are obtained using a Gutzwiller-type mean-field approach, gauged against exact diagonalization solutions of small systems.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Kagome lattice antiferromagnets and Dzyaloshinsky-Moriya interactions

We study the consequences of in-plane (D_p) and out-of-plane (D_z) Dzyaloshinsky-Moriya (DM) interactions on the thermodynamic properties of spin-1/2} Heisenberg model on the kagome lattice using numerical linked cluster expansions and exact diagonalization, and contrast them with those of other perturbations such as exchange anisotropy and dilution. We find that different combinations of the DM anisotropies lead to a wide variety of thermodynamic behavior, which are quite distinct from those of most other perturbations. We argue that the sudden upturn seen experimentally in the susceptibility of the material ZnCu3(OH)6Cl2 can be understood in terms of Dzyaloshinsky-Moriya anisotropies with D_p>|D_z|. We also show that the measured specific heat of the material puts further constraints on the allowed DM parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Bethe-Ansatz density-functional theory of ultracold repulsive fermions in one-dimensional optical lattices

We present an extensive numerical study of ground-state properties of confined repulsively interacting fermions on one-dimensional optical lattices. Detailed predictions for the atom-density profiles are obtained from parallel Kohn-Sham density-functional calculations and quantum Monte Carlo simulations. The density-functional calculations employ a Bethe-Ansatz-based local-density approximation for the correlation energy, which accounts for Luttinger-liquid and Mott-insulator physics. Semi-analytical and fully numerical formulations of this approximation are compared with each other and with a cruder Thomas-Fermi-like local-density approximation for the total energy. Precise quantum Monte Carlo simulations are used to assess the reliability of the various local-density approximations, and in conjunction with these allow to obtain a detailed microscopic picture of the consequences of the interplay between particle-particle interactions and confinement in one-dimensional systems of strongly correlated fermions.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Confinement control by optical lattices

It is shown that the interplay of a confining potential with a periodic potential leads for free particles to states spatially confined on a fraction of the total extension of the system. A more complex `slicing' of the system can be achieved by increasing the period of the lattice potential. These results are especially relevant for fermionic systems, where interaction effects are in general strongly reduced for a single species at low temperatures.

0 citations0 reviewsNo code linkedOpen access