Time-dependent currents of 1D bosons in an optical lattice

We analyse the time-dependence of currents in a 1D Bose gas in an optical lattice. For a 1D system, the stability of currents induced by accelerating the lattice exhibits a broad crossover as a function of the magnitude of the acceleration, and the strength of the inter-particle interactions. This differs markedly from mean-field results in higher dimensions. Using the infinite Time Evolving Block Decimation algorithm, we characterise this crossover by making quantitative predictions for the time-dependent behaviour of the currents and their decay rate. We also compute the time-dependence of quasi-condensate fractions which can be measured directly in experiments. We compare our results to calculations based on phase-slip methods, finding agreement with the scaling as the particle density increases, but with significant deviations near unit filling.Comment: 19 pages, 10 figure

Non-equilibrium dynamics of bosonic atoms in optical lattices: Decoherence of many-body states due to spontaneous emission

We analyze in detail the heating of bosonic atoms in an optical lattice due to incoherent scattering of light from the lasers forming the lattice. Because atoms scattered into higher bands do not thermalize on the timescale of typical experiments, this process cannot be described by the total energy increase in the system alone (which is determined by single-particle effects). The heating instead involves an important interplay between the atomic physics of the heating process and the many-body physics of the state. We characterize the effects on many-body states for various system parameters, where we observe important differences in the heating for strongly and weakly interacting regimes, as well as a strong dependence on the sign of the laser detuning from the excited atomic state. We compute heating rates and changes to characteristic correlation functions based both on perturbation theory calculations, and a time-dependent calculation of the dissipative many-body dynamics. The latter is made possible for 1D systems by combining time-dependent density matrix renormalization group (t-DMRG) methods with quantum trajectory techniques.Comment: 17 pages, 14 figure

Real-time dynamics at finite temperature by DMRG: A path-integral approach

We propose a path-integral variant of the DMRG method to calculate real-time correlation functions at arbitrary finite temperatures. To illustrate the method we study the longitudinal autocorrelation function of the $XXZ$-chain. By comparison with exact results at the free fermion point we show that our method yields accurate results up to a limiting time which is determined by the spectrum of the reduced density matrix.Comment: 5 pages, 4 figure

Time-dependent evolution of two coupled Luttinger liquids

We consider two disconnected Luttinger liquids which are coupled at $t=0$ through chiral density-density interactions. Both for $t<0$ and $t \geq 0$ the system is exactly solvable by means of bosonization and this allows to evaluate analytically the time-dependence of correlation functions. We find that in the long-time limit the critical exponent governing the one-particle correlation function differs from the exponent dictated by the equilibrium ground state of the coupled system. We also discuss how this reflects on some physical quantities which are accessible in real experiments.Comment: 6 pages, 1 eps fig, revised version accepted for publication in Phys. Rev.

Measuring entanglement growth in quench dynamics of bosons in an optical lattice

We discuss a scheme to measure the many-body entanglement growth during quench dynamics with bosonic atoms in optical lattices. By making use of a 1D or 2D setup in which two copies of the same state are prepared, we show how arbitrary order Renyi entropies can be extracted using tunnel-coupling between the copies and measurement of the parity of on-site occupation numbers, as has been performed in recent experiments. We illustrate these ideas for a Superfluid-Mott insulator quench in the Bose-Hubbard model, and also for hard-core bosons, and show that the scheme is robust against imperfections in the measurements.Comment: 4+ pages plus supplementary materia

Variational study of hard-core bosons in a 2-D optical lattice using Projected Entangled Pair States (PEPS)

We have studied the system of hard-core bosons on a 2-D optical lattice using a variational algorithm based on projected entangled-pair states (PEPS). We have investigated the ground state properties of the system as well as the responses of the system to sudden changes in the parameters. We have compared our results to mean field results based on a Gutzwiller ansatz.Comment: 9 pages, 9 figure

Dynamical crystal creation with polar molecules or Rydberg atoms in optical lattices

We investigate the dynamical formation of crystalline states with systems of polar molecules or Rydberg atoms loaded into a deep optical lattice. External fields in these systems can be used to couple the atoms or molecules between two internal states: one that is weakly interacting and one that exhibits a strong dipole-dipole interaction. By appropriate time variation of the external fields we show that it is possible to produce crystalline states of the strongly interacting states with high filling fractions chosen via the parameters of the coupling.We study the coherent dynamics of this process in one dimension (1D) using a modified form of the time-evolving block decimation (TEBD) algorithm, and obtain crystalline states for system sizes and parameters corresponding to realistic experimental configurations. For polar molecules these crystalline states will be long-lived, assisting in a characterization of the state via the measurement of correlation functions. We also show that as the coupling strength increases in the model, the crystalline order is broken. This is characterized in 1D by a change in density-density correlation functions, which decay to a constant in the crystalline regime, but show different regions of exponential and algebraic decay for larger coupling strengths

Real time evolution using the density matrix renormalization group

We describe an extension to the density matrix renormalization group method incorporating real time evolution into the algorithm. Its application to transport problems in systems out of equilibrium and frequency dependent correlation functions is discussed and illustrated in several examples. We simulate a scattering process in a spin chain which generates a spatially non-local entangled wavefunction.Comment: 4 pages, 4 eps figures, some minor corrections in text and Eq.(3

Entanglement growth in quench dynamics with variable range interactions

Studying entanglement growth in quantum dynamics provides both insight into the underlying microscopic processes and information about the complexity of the quantum states, which is related to the efficiency of simulations on classical computers. Recently, experiments with trapped ions, polar molecules, and Rydberg excitations have provided new opportunities to observe dynamics with long-range interactions. We explore nonequilibrium coherent dynamics after a quantum quench in such systems, identifying qualitatively different behavior as the exponent of algebraically decaying spin-spin interactions in a transverse Ising chain is varied. Computing the build-up of bipartite entanglement as well as mutual information between distant spins, we identify linear growth of entanglement entropy corresponding to propagation of quasiparticles for shorter range interactions, with the maximum rate of growth occurring when the Hamiltonian parameters match those for the quantum phase transition. Counter-intuitively, the growth of bipartite entanglement for long-range interactions is only logarithmic for most regimes, i.e., substantially slower than for shorter range interactions. Experiments with trapped ions allow for the realization of this system with a tunable interaction range, and we show that the different phenomena are robust for finite system sizes and in the presence of noise. These results can act as a direct guide for the generation of large-scale entanglement in such experiments, towards a regime where the entanglement growth can render existing classical simulations inefficient.Comment: 17 pages, 7 figure