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

Ian P. McCulloch

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

ResearcherUnclaimed source record
cond-mat.str-elquant-phcond-mat.quant-gascond-mat.stat-mechcond-mat.mes-hallhep-latphysics.comp-ph

Catalog footprint

What is connected

30works
7topics
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

30 published item(s)

preprint2023arXiv

Disorder-Free Localization in $2+1$D Lattice Gauge Theories with Dynamical Matter

Disorder-free localization (DFL) has been established as a mechanism of strong ergodicity breaking in $1+1$D lattice gauge theories (LGTs) with dynamical matter for quenches starting in homogeneous initial states that span an extensive number of gauge superselection sectors. Nevertheless, the fate of DFL in $2+1$D in the presence of dynamical matter has hitherto remained an open question of great interest in light of the instability of quenched-disorder many-body localization in higher spatial dimensions. Using infinite matrix product state calculations, we show that DFL survives in $2+1$D LGTs, albeit it is generally less pronounced than in $1+1$D, and highly depends on the matter configuration of the initial state. Through suitable matter configurations, we are able to relate and compare the $1+1$D and $2+1$D cases, showing that the main ingredient for the strength of DFL in our setup is the \textit{propagation directionality} of matter. Our results suggest that, generically, DFL is weakened with increasing spatial dimension, although it can be made independent of the latter by minimizing the propagation directionality of matter in the initial state.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Dynamical phase transitions in the two-dimensional transverse-field Ising model

We investigate two separate notions of dynamical phase transitions in the two-dimensional nearest-neighbor transverse-field Ising model on a square lattice using matrix product states and a new \textit{hybrid} infinite time-evolving block decimation algorithm, where the model is implemented on an infinitely long cylinder with a finite diameter along which periodic boundary conditions are employed. Starting in an ordered initial state, our numerical results suggest that quenches below the dynamical critical point give rise to a ferromagnetic long-time steady state with the Loschmidt return rate exhibiting \textit{anomalous} cusps even when the order parameter never crosses zero. Within the accessible timescales of our numerics, quenches above the dynamical critical point suggest a paramagnetic long-time steady state with the return rate exhibiting \textit{regular} cusps connected to zero crossings of the order parameter. Additionally, our simulations indicate that quenching slightly above the dynamical critical point leads to a coexistence region where both anomalous and regular cusps appear in the return rate. Quenches from the disordered phase further confirm our main conclusions. Our work supports the recent finding that anomalous cusps arise only when local spin excitations are the energetically dominant quasiparticles. Our results are accessible in modern Rydberg experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Robustness of gauge-invariant dynamics against defects in ultracold-atom gauge theories

Recent years have seen strong progress in quantum simulation of gauge-theory dynamics using ultracold-atom experiments. A principal challenge in these efforts is the certification of gauge invariance, which has recently been realized in [B.~Yang et al., arXiv:2003.08945]. One major but poorly investigated experimental source of gauge-invariance violation is an imperfect preparation of the initial state. Using the time-dependent density-matrix renormalization group, we analyze the robustness of gauge-invariant dynamics against potential preparation defects in the above ultracold-atom implementation of a $\mathrm{U}(1)$ gauge theory. We find defects related to an erroneous initialization of matter fields to be innocuous, as the associated gauge-invariance violation remains strongly localized throughout the time evolution. A defect due to faulty initialization of the gauge field leads to a mild proliferation of the associated violation. Furthermore, we characterize the influence of immobile and mobile defects by monitoring the spread of entanglement entropy. Overall, our results indicate that the aforementioned experimental realization exhibits a high level of fidelity in the gauge invariance of its dynamics at all evolution times. Our work provides strong evidence that ultracold-atom setups can serve as an extremely reliable framework for the quantum simulation of gauge-theory dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Cumulants and Scaling Functions of Infinite Matrix Product States

The order parameter cumulants of infinite matrix product ground states are evaluated across a quantum phase transition. A scheme using the Binder cumulant, finite-entanglement scaling and scaling functions to obtain the critical point and exponents of the correlation length and cumulants is presented. Analogous to the scaling relations that relate the exponents of various thermodynamic quantities, a cumulant exponent relation is derived and used to check the consistency and relationship between the cumulant exponents. This scheme gives a numerically economical way of accurately obtaining the critical exponents. Examples of this scheme are shown for four one-dimensional models - the transverse field Ising model, the topological Kondo insulator, the S = 1 Heisenberg chain with single-ion anisotropy and the Bose-Hubbard model. A two-dimensional model is also exemplified in the square lattice transverse field Ising model on an infinite cylinder. These exemplary systems portray a variety of local and string order parameters as well as phase transition classes that can be studied with the scaling functions and infinite matrix product states.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Hybrid infinite time-evolving block decimation algorithm for long-range multi-dimensional quantum many-body systems

In recent years, the infinite time-evolution block decimation (iTEBD) method has been demonstrated to be one of the most efficient and powerful numerical schemes for time-evolution in one-dimensional quantum many-body systems. However, a major shortcoming of the method, along with other state-of-the-art algorithms for many-body dynamics, has been their restriction to one spatial dimension. We present an algorithm based on a \textit{hybrid} extension of iTEBD where finite blocks of a chain are first locally time-evolved before an iTEBD-like method combines these processes globally. This in turn permits simulating the dynamics of many-body systems in the thermodynamic limit in $d\geq1$ dimensions including in the presence of long-range interactions. Our work paves the way for simulating the dynamics of many-body phenomena that occur exclusively in higher dimensions, and whose numerical treatments have hitherto been limited to exact diagonalization of small systems, which fundamentally limits a proper investigation of dynamical criticality. We expect the algorithm presented here to be of significant importance to validating and guiding investigations in state-of-the-art ion-trap and ultracold-atom experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Optical conductivity of the Hubbard chain away from half filling

We consider the optical conductivity $σ_1(ω)$ in the metallic phase of the one-dimensional Hubbard model. Our results focus on the vicinity of half filling and the frequency regime around the optical gap in the Mott insulating phase. By means of a density-matrix renormalization group implementation of the correction-vector approach, $σ_1(ω)$ is computed for a range of interaction strengths and dopings. We identify an energy scale $E_{\rm opt}$ above which the optical conductivity shows a rapid increase. We then use a mobile impurity model in combination with exact results to determine the behavior of $σ_1(ω)$ for frequencies just above $E_{\rm opt}$ which is in agreement with our numerical data. As a main result, we find that this onset behavior is not described by a power law.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Strictly Single-Site DMRG Algorithm with Subspace Expansion

We introduce a strictly single-site DMRG algorithm based on the subspace expansion of the Alternating Minimal Energy (AMEn) method. The proposed new MPS basis enrichment method is sufficient to avoid local minima during the optimisation, similarly to the density matrix perturbation method, but computationally cheaper. Each application of $\hat H$ to $|Ψ\rangle$ in the central eigensolver is reduced in cost for a speed-up of $\approx (d + 1)/2$, with $d$ the physical site dimension. Further speed-ups result from cheaper auxiliary calculations and an often greatly improved convergence behaviour. Runtime to convergence improves by up to a factor of 2.5 on the Fermi-Hubbard model compared to the previous single-site method and by up to a factor of 3.9 compared to two-site DMRG. The method is compatible with real-space parallelisation and non-abelian symmetries.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Fast convergence of imaginary time evolution tensor network algorithms by recycling the environment

We propose an environment recycling scheme to speed up a class of tensor network algorithms that produce an approximation to the ground state of a local Hamiltonian by simulating an evolution in imaginary time. Specifically, we consider the time-evolving block decimation (TEBD) algorithm applied to infinite systems in 1D and 2D, where the ground state is encoded, respectively, in a matrix product state (MPS) and in a projected entangled-pair state (PEPS). An important ingredient of the TEBD algorithm (and a main computational bottleneck, especially with PEPS in 2D) is the computation of the so-called environment, which is used to determine how to optimally truncate the bond indices of the tensor network so that their dimension is kept constant. In current algorithms, the environment is computed at each step of the imaginary time evolution, to account for the changes that the time evolution introduces in the many-body state represented by the tensor network. Our key insight is that close to convergence, most of the changes in the environment are due to a change in the choice of gauge in the bond indices of the tensor network, and not in the many-body state. Indeed, a consistent choice of gauge in the bond indices confirms that the environment is essentially the same over many time steps and can thus be re-used, leading to very substantial computational savings. We demonstrate the resulting approach in 1D and 2D by computing the ground state of the quantum Ising model in a transverse magnetic field.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Imaginary-time matrix product state impurity solver for dynamical mean-field theory

We present a new impurity solver for dynamical mean-field theory based on imaginary-time evolution of matrix product states. This converges the self-consistency loop on the imaginary-frequency axis and obtains real-frequency information in a final real-time evolution. Relative to computations on the real-frequency axis, required bath sizes are much smaller and less entanglement is generated, so much larger systems can be studied. The power of the method is demonstrated by solutions of a three band model in the single and two-site dynamical mean-field approximation. Technical issues are discussed, including details of the method, efficiency as compared to other matrix product state based impurity solvers, bath construction and its relation to real-frequency computations and the analytic continuation problem of quantum Monte Carlo, the choice of basis in dynamical cluster approximation, and perspectives for off-diagonal hybridization functions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Non-thermal melting of Neel order in the Hubbard model

We study the unitary time evolution of antiferromagnetic order in the Hubbard model after a quench starting from the perfect Néel state. In this setup, which is well suited for experiments with cold atoms, one can distinguish fundamentally different pathways for melting of long-range order at weak and strong interaction. In the Mott insulating regime, melting of long-range order occurs due to the ultra-fast transfer of energy from charge excitations to the spin background, while local magnetic moments and their exchange coupling persist during the process. The latter can be demonstrated by a local spin-precession experiment. At weak interaction, local moments decay along with the long-range order. The dynamics is governed by residual quasiparticles, which are reflected in oscillations of the off-diagonal components of the momentum distribution. Such oscillations provide an alternative route to study the prethermalization phenomenon and its influence on the dynamics away from the integrable (noninteracting) limit. The Hubbard model is solved within nonequilibrium dynamical mean-field theory, using the density matrix-renormalization group as an impurity solver.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Phase diagram of the J1-J2 Heisenberg model on the kagome lattice

We perform an extensive density matrix renormalization group (DMRG) study of the ground-state phase diagram of the spin-1/2 J_1-J_2 Heisenberg model on the kagome lattice. We focus on the region of the phase diagram around the kagome Heisenberg antiferromagnet, i.e., at J_2=0. We investigate the static spin structure factor, the magnetic correlation lengths, and the spin gaps. Our results are consistent with the absence of magnetic order in a narrow region around J_2\approx 0, although strong finite-size effects do not allow us to accurately determine the phase boundaries. This result is in agreement with the presence of an extended spin-liquid region, as it has been proposed recently. Outside the disordered region, we find that for ferromagnetic and antiferromagnetic J_2 the ground state displays signatures of the magnetic order of the \sqrt{3}\times\sqrt{3} and the q=0 type, respectively. Finally, we focus on the structure of the entanglement spectrum (ES) in the q=0 ordered phase. We discuss the importance of the choice of the bipartition on the finite-size structure of the ES.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Spectral functions and time evolution from the Chebyshev recursion

We link linear prediction of Chebyshev and Fourier expansions to analytic continuation. We push the resolution in the Chebyshev-based computation of $T=0$ many-body spectral functions to a much higher precision by deriving a modified Chebyshev series expansion that allows to reduce the expansion order by a factor $\sim\frac{1}{6}$. We show that in a certain limit the Chebyshev technique becomes equivalent to computing spectral functions via time evolution and subsequent Fourier transform. This introduces a novel recursive time evolution algorithm that instead of the group operator $e^{-iHt}$ only involves the action of the generator $H$. For quantum impurity problems, we introduce an adapted discretization scheme for the bath spectral function. We discuss the relevance of these results for matrix product state (MPS) based DMRG-type algorithms, and their use within dynamical mean-field theory (DMFT). We present strong evidence that the Chebyshev recursion extracts less spectral information from $H$ than time evolution algorithms when fixing a given amount of created entanglement.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Chebyshev Matrix Product State Impurity Solver for the Dynamical Mean-Field Theory

We compute the spectral functions for the two-site dynamical cluster theory and for the two-orbital dynamical mean-field theory in the density-matrix renormalization group (DMRG) framework using Chebyshev expansions represented with matrix product states (MPS). We obtain quantitatively precise results at modest computational effort through technical improvements regarding the truncation scheme and the Chebyshev rescaling procedure. We furthermore establish the relation of the Chebyshev iteration to real-time evolution, and discuss technical aspects as computation time and implementation in detail.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Comment on "Phase separation in a two-species Bose mixture"

In an article in 2007, Mishra, Pai, and Das [Phys. Rev. A 76, 013604 (2007)] investigated the two-component Bose-Hubbard model using the numerical DMRG procedure. In the regime of inter-species repulsion $U^{ab}$ larger than the intra-species repulsion $U$, they found a transition from a uniform miscible phase to phase-separation occurring at a finite value of $U$ , e.g., at around $U = 1.3$ for $Δ= U^{ab}/U = 1.05$ and $ρ_{a} = ρ_{b} = 1/2$. In this comment, we show that this result is not correct and in fact the two-component Bose-Hubbard model is unstable to phase-separation for any $U^{ab} > U > 0$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Domain-wall melting in ultracold boson systems with holes and spin-flip defects

Quantum magnetism is a fundamental phenomenon of nature. As of late, it has garnered a lot of interest because experiments with ultracold atomic gases in optical lattices could be used as a simulator for phenomena of magnetic systems. A paradigmatic example is the time evolution of a domain-wall state of a spin-1/2 Heisenberg chain, the so-called domain-wall melting. The model can be implemented by having two species of bosonic atoms with unity filling and strong on-site repulsion U in an optical lattice. In this paper, we study the domain-wall melting in such a setup on the basis of the time-dependent density matrix renormalization group (tDMRG). We are particularly interested in the effects of defects that originate from an imperfect preparation of the initial state. Typical defects are holes (empty sites) and flipped spins. We show that the dominating effects of holes on observables like the spatially resolved magnetization can be taken account of by a linear combination of spatially shifted observables from the clean case. For sufficiently large U, further effects due to holes become negligible. In contrast, the effects of spin flips are more severe as their dynamics occur on the same time scale as that of the domain-wall melting itself. It is hence advisable to avoid preparation schemes that are based on spin-flips.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Quantum magnetism of bosons with synthetic gauge fields in one-dimensional optical lattices: a Density Matrix Renormalization Group study

In this paper, we provide a comprehensive study of the quantum magnetism in the Mott insulating phases of the 1D Bose-Hubbard model with abelian or non-abelian synthetic gauge fields, using the Density Matrix Renormalization Group (DMRG) method. We focus on the interplay between the synthetic gauge field and the asymmetry of the interactions, which give rise to a very general effective magnetic model: a XYZ model with various Dzyaloshinskii-Moriya (DM) interactions. The properties of the different quantum magnetic phases and phases transitions of this model are investigated.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Solving nonequilibrium dynamical mean-field theory using matrix product states

We solve the nonequilibrium dynamical mean-field theory (DMFT) using matrix product states (MPS). This allows us to treat much larger bath sizes and by that reach substantially longer times (factor $\sim$ 2 -- 3) than with exact diagonalization. We show that the star geometry of the underlying impurity problem can have substantially better entanglement properties than the previously favoured chain geometry. This has immense consequences for the efficiency of an MPS-based description of general impurity problems: in the case of equilibrium DMFT, it leads to an orders-of-magnitude speedup. We introduce an approximation for the two-time hybridization function that uses time-translational invariance, which can be observed after a certain relaxation time after a quench to a time-independent Hamiltonian.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Miscible-Immiscible Quantum Phase Transition in Coupled Two-Component Bose-Einstein Condensates in 1D Optical Lattices

Using numerical techniques, we study the miscible-immiscible quantum phase transition in a linearly coupled binary Bose-Hubbard model Hamiltonian that can describe low-energy properties of a two-component Bose-Einstein condensate in optical lattices. With the quantum many-body ground state obtained from density matrix renormalization group algorithm, we calculate the characteristic physical quantities of the phase transition controlled by the linear coupling between two components. Furthermore we calculate the Binder cumulant to determine the critical point and draw the phase diagram. The strong-coupling expansion shows that in the Mott insulator regime the model Hamiltonian can be mapped to a spin 1/2 XXZ model with a transverse magnetic field.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Creation and dynamics of remote spin-entangled pairs in the expansion of strongly correlated fermions in an optical lattice

We consider the nonequilibrium dynamics of an interacting spin-1/2 fermion gas in a one-dimensional optical lattice after switching off the confining potential. In particular, we study the creation and the time evolution of spatially separated, spin-entangled fermionic pairs. The time-dependent density-matrix renormalization group is used to simulate the time evolution and evaluate the two-site spin correlation functions, from which the concurrence is calculated. We find that the typical distance between entangled fermions depends crucially on the onsite interaction strength, and that a time-dependent modulation of the tunneling amplitude can enhance the production of spin-entanglement. Moreover, we discuss the prospects of experimentally observing these phenomena using spin-dependent single-site detection.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Expansion dynamics of interacting bosons in homogeneous lattices in one and two dimensions

We experimentally and numerically investigate the expansion of initially localized ultracold bosons in homogeneous one- and two-dimensional optical lattices. We find that both dimensionality and interaction strength crucially influence these non-equilibrium dynamics. While the atoms expand ballistically in all integrable limits, deviations from these limits dramatically suppress the expansion and lead to the appearance of almost bimodal cloud shapes, indicating diffusive dynamics in the center surrounded by ballistic wings. For strongly interacting bosons, we observe a dimensional crossover of the dynamics from ballistic in the one-dimensional hard-core case to diffusive in two dimensions, as well as a similar crossover when higher occupancies are introduced into the system.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Dynamical windows for real-time evolution with matrix product states

We propose the use of a dynamical window to investigate the real-time evolution of quantum many-body systems in a one-dimensional lattice. In a recent paper [H. Phien et al, arxiv:????.????], we introduced infinite boundary conditions (IBC) in order to investigate real-time evolution of an infinite system under a local perturbation. This was accomplished by restricting the update of the tensors in the matrix product state to a finite window, with left and right boundaries held at fixed positions. Here we consider instead the use of a dynamical window, namely a window where the positions of left and right boundaries are allowed to change in time. In this way, all simulation efforts can be devoted to the space-time region of interest, which leads to a remarkable reduction in computational costs. For illustrative purposes, we consider two applications in the context of the spin-1 antiferromagnetic Heisenberg model in an infinite spin chain: one is an expanding window, with boundaries that are adjusted to capture the expansion in time of a local perturbation of the system; the other is a moving window of fixed size, where the position of the window follows the front of a propagating wave.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Infinite boundary conditions for matrix product state calculations

We propose a formalism to study dynamical properties of a quantum many-body system in the thermodynamic limit by studying a finite system with infinite boundary conditions (IBC) where both finite size effects and boundary effects have been eliminated. For one-dimensional systems, infinite boundary conditions are obtained by attaching two boundary sites to a finite system, where each of these two sites effectively represents a semi-infinite extension of the system. One can then use standard finite-size matrix product state techniques to study a region of the system while avoiding many of the complications normally associated with finite-size calculations such as boundary Friedel oscillations. We illustrate the technique with an example of time evolution of a local perturbation applied to an infinite (translationally invariant) ground state, and use this to calculate the spectral function of the S=1 Heisenberg spin chain. This approach is more efficient and more accurate than conventional simulations based on finite-size matrix product state and density-matrix renormalization-group approaches.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Nature of the Spin Liquid Ground State of the S=1/2 Kagome Heisenberg Model

We perform a density-matrix renormalization group (DMRG) study of the S=1/2 Heisenberg antiferromagnet on the kagome lattice to identify the conjectured spin liquid ground state. Exploiting SU(2) spin symmetry, which allows us to keep up to 16,000 DMRG states, we consider cylinders with circumferences up to 17 lattice spacings and find a spin liquid ground state with an estimated per site energy of -0.4386(5), a spin gap of 0.13(1), very short-range decay in spin, dimer and chiral correlation functions and finite topological entanglement gamma consistent with gamma=ld 2, ruling out gapless, chiral or non-topological spin liquids. At the same time, DMRG results provide strong evidence for a gapped topological Z2 spin liquid.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Stroboscopic observation of quantum many-body dynamics

Recent experiments have demonstrated single-site resolved observation of cold atoms in optical lattices. Thus, in the future it may be possible to take repeated snapshots of an interacting quantum many-body system during the course of its evolution. Here we address the impact of the resulting Quantum (anti-)Zeno physics on the many-body dynamics. We use the time-dependent density-matrix renormalization group to obtain the time evolution of the full many-body wave function, which is then periodically projected in order to simulate realizations of stroboscopic measurements. For the example of a one-dimensional lattice of spin-polarized fermions with nearest-neighbor interactions, we find regimes for which many-particle configurations are stabilized and destabilized depending on the interaction strength and the time between observations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Chebyshev matrix product state approach for spectral functions

We show that recursively generated Chebyshev expansions offer numerically efficient representations for calculating zero-temperature spectral functions of one-dimensional lattice models using matrix product state (MPS) methods. The main features of this Chebychev matrix product state (CheMPS) approach are: (i) it achieves uniform resolution over the spectral function's entire spectral width; (ii) it can exploit the fact that the latter can be much smaller than the model's many-body bandwidth; (iii) it offers a well-controlled broadening scheme; (iv) it is based on a succession of Chebychev vectors |t_n>, (v) whose entanglement entropies were found to remain bounded with increasing recursion order n for all cases analyzed here; (vi) it distributes the total entanglement entropy that accumulates with increasing n over the set of Chebyshev vectors |t_n>. We present zero-temperature CheMPS results for the structure factor of spin-1/2 antiferromagnetic Heisenberg chains and perform a detailed finite-size analysis. Making comparisons to three benchmark methods, we find that CheMPS (1) yields results comparable in quality to those of correction vector DMRG, at dramatically reduced numerical cost; (2) agrees well with Bethe Ansatz results for an infinite system, within the limitations expected for numerics on finite systems; (3) can also be applied in the time domain, where it has potential to serve as a viable alternative to time-dependent DMRG (in particular at finite temperatures). Finally, we present a detailed error analysis of CheMPS for the case of the noninteracting resonant level model.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Landau-Zener sweeps and sudden quenches in coupled Bose-Hubbard chains

We simulate numerically the dynamics of strongly correlated bosons in a two-leg ladder subject to a time-dependent energy bias between the two chains. When all atoms are initially in the leg with higher energy, we find a drastic reduction of the inter-chain particle transfer for slow linear sweeps, in quantitative agreement with recent experiments. This effect is preceded by a rapid broadening of the quasi-momentum distribution of atoms, signaling the presence of a bath of low-energy excitations in the chains. We further investigate the scenario of quantum quenches to fixed values of the energy bias. We find that for large enough density the momentum distribution relaxes to that of an equilibrium thermal state with the same energy.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Probing the relaxation towards equilibrium in an isolated strongly correlated 1D Bose gas

The problem of how complex quantum systems eventually come to rest lies at the heart of statistical mechanics. The maximum entropy principle put forward in 1957 by E. T. Jaynes suggests what quantum states one should expect in equilibrium but does not hint as to how closed quantum many-body systems dynamically equilibrate. A number of theoretical and numerical studies accumulate evidence that under specific conditions quantum many-body models can relax to a situation that locally or with respect to certain observables appears as if the entire system had relaxed to a maximum entropy state. In this work, we report the experimental observation of the non-equilibrium dynamics of a density wave of ultracold bosonic atoms in an optical lattice in the regime of strong correlations. Using an optical superlattice, we are able to prepare the system in a well-known initial state with high fidelity. We then follow the dynamical evolution of the system in terms of quasi-local densities, currents, and coherences. Numerical studies based on the time-dependent density-matrix renormalization group method are in an excellent quantitative agreement with the experimental data. For very long times, all three local observables show a fast relaxation to equilibrium values compatible with those expected for a global maximum entropy state. We find this relaxation of the quasi-local densities and currents to initially follow a power-law with an exponent being significantly larger than for free or hardcore bosons. For intermediate times the system fulfills the promise of being a dynamical quantum simulator, in that the controlled dynamics runs for longer times than present classical algorithms based on matrix product states can efficiently keep track of.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Real-time energy dynamics in spin-1/2 Heisenberg chains

We study the real-time dynamics of the local energy density in the spin-1/2 XXZ chain starting from initial states with an inhomogeneous profile of bond energies. Numerical simulations of the dynamics of the initial states are carried out using the adaptive time-dependent density matrix renormalization group method. We analyze the time dependence of the spatial variance associated with the local energy density to classify the dynamics as either ballistic or diffusive. Our results are consistent with ballistic behavior both in the massless and the massive phase. We also study the same problem within Luttinger Liquid theory and obtain that energy wave-packets propagate with the sound velocity. We recover this behavior in our numerical simulations in the limit of very weakly perturbed initial states.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Coulomb interaction effects and electron spin relaxation in the 1d Kondo lattice model

We study the effects of the Coulomb interaction in the one dimensional Kondo lattice model on the phase diagram, the static magnetic susceptibility and electron spin relaxation. We show that onsite Coulomb interaction supports ferromagnetic order and nearest neighbor Coulomb interaction drives, depending on the electron filling, either a paramagnetic or ferromagnetic order. Furthermore we calculate electron quasiparticle life times, which can be related to electron spin relaxation and decoherence times, and explain their dependence on the strength of interactions and the electron filling in order to find the sweet spot of parameters where the relaxation time is maximized. We find that effective exchange processes between the electrons dominate the spin relaxation and decoherence rate.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Valence bond entanglement entropy of frustrated spin chains

We extend the definition of the recently introduced valence bond entanglement entropy to arbitrary SU(2) wave functions of S=1/2 spin systems. Thanks to a reformulation of this entanglement measure in terms of a projection, we are able to compute it with various numerical techniques for frustrated spin models. We provide extensive numerical data for the one-dimensional J1-J2 spin chain where we are able to locate the quantum phase transition by using the scaling of this entropy with the block size. We also systematically compare with the scaling of the von Neumann entanglement entropy. We finally underline that the valence-bond entropy definition does depend on the choice of bipartition so that, for frustrated models, a "good" bipartition should be chosen, for instance according to the Marshall sign.

0 citations0 reviewsNo code linkedOpen access