We investigate the propagation of density-wave packets in a Bose-Hubbard model using the adaptive time-dependent density-matrix renormalization group method. We discuss the decay of the amplitude with time and the dependence of the velocity on density, interaction strength, and the height of the perturbation in a numerically exact way, covering arbitrary interactions and amplitudes of the perturbation. In addition, we investigate the effect of self-steepening due to the amplitude dependence of the velocity and discuss the possibilities for an experimental detection of the moving wave packet in time-of-flight pictures. By comparing the sound velocity to theoretical predictions, we determine the limits of a Gross-Pitaevskii- or Bogoliubov-type description and the regime where repulsive one-dimensional Bose gases exhibit fermionic behavior.
We investigate the propagation of density-wave packets in a Bose-Hubbard model using the adaptive time-dependent density-matrix renormalization group method. We discuss the decay of the amplitude with time and the dependence of the velocity on density, interaction strength, and the height of the perturbation in a numerically exact way, covering arbitrary interactions and amplitudes of the perturbation. In addition, we investigate the effect of self-steepening due to the amplitude dependence of the velocity and discuss the possibilities for an experimental detection of the moving wave packet in time-of-flight pictures. By comparing the sound velocity to theoretical predictions, we determine the limits of a Gross-Pitaevskii- or Bogoliubov-type description and the regime where repulsive one-dimensional Bose gases exhibit fermionic behavior.
In this work we develop a new numerical method, the adaptive time-dependent density-matrix renormalization-group (adaptive t-DMRG), which turns out to be very well suited to investigate time-dependent phenomena in one-dimensional strongly correlated systems. The method is based on the original density-matrix renormalization group method by White and the time evolving block-decimation procedure by Vidal. We further show the applicability of the new method to different physical systems: bosonic, fermionic, and spin chains. We first discuss the consequences of the presence of an external trapping potential on the static properties of trapped cold bosonic atoms in an optical lattice using the Bose-Hubbard model. It is shown that a properly rescaled one-particle density matrix characterizes superfluid versus insulating states just as in the homogeneous system. We note that the superfluid to Mott-insulating transition is seen most directly in the half width of the interference peak. Then we turn to the evolution of density perturbations in these systems. In particular, we discuss the dependence of the velocity and the decay of the amplitude on density, interaction strength, and the extent and height of the perturbation in a numerically exact way, covering a wide range of interaction and perturbation strengths. By comparing our results for the sound velocity to theoretical predictions, we determine the limits of a Gross-Pitaevskii or Bogoliubov type description and the regime where repulsive one-dimensional Bose gases exhibit fermionic behaviour. As a second application of the adaptive t-DMRG, we investigate the phenomenon of spin-charge separation in the Hubbard model and propose an experimental setup for its observation in cold Fermi gases. We show the robustness of this separation beyond the low-energy regime by studying the time evolution of density wave packets of finite strength and at length scales down to a few lattice spacings. A striking signature of spin-charge separation is found in 1D cold Fermi gases in a harmonic trap using the different propagation properties of the liquid and Mott-insulating phases. The third application relates to the area of spintronics. A simplified model for the magnetization transport between two coupled reservoirs with opposite spin polarization is studied. We investigate whether a simple long time limit for the spin transport exists and if so, how it depends on the properties of the system such as integrability and criticality. Time-scales accessible to us are of the order of 100 units of time measured in 1/J while maintaining insignificant error in the observables.
We study the propagation of sound waves in a Bose-Einstein condensate trapped in a one-dimensional optical lattice. We find that the velocity of the propagation of sound wave packets decreases with increasing optical lattice depth, as predicted by the Bogoliubov theory. The strong interplay between nonlinearities and the periodicity of the external potential generates phenomena that are not present in the uniform case. Shock waves, for instance, can propagate slower than sound waves, due to the negative curvature of the dispersion relation. Moreover, nonlinear corrections to the Bogoliubov theory appear to be important even with very small density perturbations, inducing a saturation of the amplitude of the sound signal.
We propose a setup in which Andreev-like reflections predicted for 1D transport systems could be observed time dependently using cold atoms in a 1D optical lattice. Using time-dependent density matrix renormalization group methods we analyze the wave packet dynamics as a density excitation propagates across a boundary in the interaction strength. These phenomena exhibit good correspondence with predictions from Luttinger liquid models and could be observed in current experiments in the context of the Bose-Hubbard model.
We analyze the time evolution of the Bose-Hubbard model after a sudden quantum quench to a weakly interacting regime. Specifically, motivated by a recent experiment at Kyoto University, we numerically simulate redistribution of the kinetic and on-site interaction energies at an early time, which was observed in nonequilibrium dynamics of ultracold Bose gases in a cubic optical lattice starting with a singly occupied Mott-insulator state. In order to compute the short-time dynamics corresponding to the experimental situation, we apply the truncated-Wigner approximation to the Bose-Hubbard model on a cubic lattice. We show that our semiclassical approach quantitatively reproduces the fast redistribution dynamics. We further analyze spatial spreading of density-density correlations at equal time in the Bose-Hubbard model on a square lattice with a large filling factor. When the system is initially prepared in a coherent state, we find that a propagation velocity of the correlation wave packet in the correlation function strongly depends on the final interaction strength, and it is bounded by twice the maximum group velocity of the elementary excitations. In contrast, when the system is initially in a Mott-insulator state, the propagation velocity of the wave packet is approximately independent of the final interaction strength.
Ultracold bosonic gases in optical lattices are strongly correlated quantum systems similar to solids. The strong correlation between the electrons in a solid on the one hand, and the bosonic atoms in optical lattice on the other, exhibit various quantum phenomena like insulation, conductivity, localization of electrons and atoms, respectively. Controlled by the intensity of the lattice laser, the ultracold bosonic gas can be transferred from a regime with superfluid character for shallow lattices into a regime of strong correlations, the Mott insulator. As an additional external parameter besides the lattice depth, one can generate spatial inhomogeneities by superimposing an additional standing wave (so-called two-color superlattices), which gives rise to localization effects or the formation of a Bose-glass phase. In the present work, numerical simulations are employed in order to investigate characteristic signatures of the quantum phases in the low-energy excitation spectrum of one-dimensional systems. We simulate temporal small amplitude modulations of the optical lattice in analogy to experiments, and evaluate the response of the system from the time-evolved initial state. The lattice systems are described in the framework of the Bose-Hubbard model. For the evaluation of the time-evolved state, we employ several numerical methods. We analyze systems of small size (6 particles on 6 sites) using an exact time-evolution by integration of the time-dependent Schrödinger equation. The formulation of an importance truncation scheme enables us to retain only the relevant components of the model space in the strongly correlated regime and, thus, allows for the investigation of systems with 10 particles on 10 sites using exact time-evolution. Based on this method, we present results of the Mott-insulating regime as well as for the Bose-glass phase. Furthermore, we employ particle-hole methods, which allow for the treatment of systems with experimental lattice sizes and particle numbers. Starting from the equation of motion method we adapt the Tamm-Dancoff approximation as well as the random-phase approximation for the occupation number representation of the Bose-Hubbard model. We present results of simulations of up to 50 particles on 50 sites and discuss the impact of the lattice depth on the low-energy excitations (U-resonance). Moreover, the impact of a two-color superlattice and the variation of its amplitude is investigated.
We analyze the quench dynamics of a one-dimensional bosonic Mott insulator and focus on the time evolution of density correlations. For these we identify a pronounced propagation front, the velocity of which, once correctly extrapolated at large distances, can serve as a quantitative characteristic of the many-body Hamiltonian. In particular, the velocity allows the weakly interacting regime, which is qualitatively well described by free bosons, to be distinguished from the strongly interacting one, in which pairs of distinct quasiparticles dominate the dynamics. In order to describe the latter case analytically, we introduce a general approximation to solve the Bose-Hubbard Hamiltonian based on the Jordan-Wigner fermionization of auxiliary particles. This approach can also be used to determine the ground-state properties. As a complement to the fermionization approach, we derive explicitly the time-dependent many-body state in the noninteracting limit and compare our results to numerical simulations in the whole range of interactions of the Bose-Hubbard model.