In a recent Letter [Phys. Rev. Lett. 88, 256403(2002), cond-mat/0109158] Cazalilla and Marston proposed a time-dependent density- matrix renormalization group (TdDMRG) algorithm for the accurate evaluation of out-of-equilibrium properties of quantum many-body systems. For a point contact junction between two Luttinger liquids, a current oscillation develops after initial transient in the insulating regime. Here we would like to point out that (a) the observed oscillation is an artifact of the method; (b) the TdDMRG can be significantly improved by incorporating the non-equilibrium evolution of the goundstate into the density matrix.
In a recent Letter [Phys. Rev. Lett. 88, 256403(2002), cond-mat/0109158] Cazalilla and Marston proposed a time-dependent density- matrix renormalization group (TdDMRG) algorithm for the accurate evaluation of out-of-equilibrium properties of quantum many-body systems. For a point contact junction between two Luttinger liquids, a current oscillation develops after initial transient in the insulating regime. Here we would like to point out that (a) the observed oscillation is an artifact of the method; (b) the TdDMRG can be significantly improved by incorporating the non-equilibrium evolution of the goundstate into the density matrix.
The density-matrix renormalization-group algorithm is extended to treat time-dependent problems. The method provides a systematic and robust tool to explore out-of-equilibrium phenomena in quantum many-body systems. We illustrate the method by showing that attractive interactions enhance the tunneling current between two Luttinger liquids, whereas repulsive interactions suppress it, in qualitative agreement with analytical predictions. Enhancement of the transport current through a quantum dot in the Kondo regime is also exhibited.
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 develop an alternative time-dependent numerical renormalization group (TDNRG) formalism for multiple quenches and implement it to study the response of a quantum impurity system to a general pulse. Within this approach, we reduce the contribution of the NRG approximation to numerical errors in the time evolution of observables by a formulation that avoids the use of the generalized overlap matrix elements in our previous multiple-quench TDNRG formalism [Nghiem et al., Phys. Rev. B 89, 075118 (2014); Phys. Rev. B 90, 035129 (2014)]. We demonstrate that the formalism yields a smaller cumulative error in the trace of the projected density matrix as a function of time and a smaller discontinuity of local observables between quenches than in our previous approach. Moreover, by increasing the switch-on time, the time between the first and last quench of the discretized pulse, the long-time limit of observables systematically converges to its expected value in the final state, i.e., the more adiabatic the switching, the more accurately is the long-time limit recovered. The present formalism can be straightforwardly extended to infinite switch-on times. We show that this yields highly accurate results for the long-time limit of both thermodynamic observables and spectral functions, and overcomes the significant errors within the single quench formalism [Anders et al., Phys. Rev. Lett. 95, 196801 (2005); Nghiem et al., Phys. Rev. Lett. 119, 156601 (2017)]. This improvement provides a first step towards an accurate description of nonequilibrium steady states of quantum impurity systems, e.g., within the scattering states NRG approach [Anders, Phys. Rev. Lett. 101, 066804 (2008)].
We present a detailed comparison of three different methods designed to tackle non-equilibrium quantum transport, namely the functional renormalization group (fRG), the time-dependent density matrix renormalization group (tDMRG) and the iterative summation of real-time path integrals (ISPI). For the non-equilibrium single-impurity Anderson model (including a Zeeman term at the impurity site), we demonstrate that the three methods are in quantitative agreement over a wide range of parameters at the particle–hole symmetric point as well as in the mixed-valence regime. We further compare these techniques with two quantum Monte Carlo approaches and the time-dependent numerical renormalization group method.
In this thesis, strongly correlated quantum many-body systems in equilibrium and in out-of-equilibrium situations are investigated. This is done by applying and developing well established numerical methods. The focus of the thesis lies in the development and application of the density matrix renormalization group method (DMRG) to quantum many-body systems out of equilbrium. In this thesis, in addition to the DMRG, methods for the exact diagonalization of the Hamiltonian of the system, like the Lanczos- or the Jacobi-Davidson method, are treated. An extension of the Lanczos method makes it possible to treat the time evolution of strongly correlated quantum systems with an accuracy comparable to machine precision. This method is the basis for a possible extension of the DMRG for the treatment of systems out of equilibrium, the so-called adaptive time-dependent DMRG ("adaptive t-DMRG"). A second variant of the adaptive t-DMRG uses the Suzuki-Trotter decomposition of the time-evolution operator. An error analysis demonstrates that both methods, for a suitable choice of control parameters, have errors < 1% at the end of the time evolution in complicated quantities like, e.g., the momentum distribution function. We apply these numerical methods to investigate the quantum critical behavior of a variant of the Hubbard model and to treat two non-equilibrium situations of current interest. Extensive use of the DMRG makes it possible to clarify the quantum-critical behavior of the ionic Hubbard model; in particular, the numerical results demonstrate quite clearly that the correct scenario has two critical points; at the first critical point, only the charge degrees of freedom are critical, while at the second one only the spin degrees of freedom are critical. Quite surprisingly, the Mott-insulator phase in the strong-coupling regime shows a divergent electrical susceptibility. Next, the dynamics of a system of so-called soft-core bosons is treated. The particles initially are trapped in a deep box-potential and the dynamics is investigated after releasing them from this trapping potential. Similar to exact results obtained for so-called hard-core bosons, (quasi-)coherent matter waves emerge, demonstrating the possibility of realizing an atom laser in experiments on optical lattices. We find that the wave vector of the emerging matter wave can be tuned by changing the strength of the interaction between the particles. The second non-equilibrium situation is a so-called "quantum-quench" of a system of strongly correlated fermions, i.e., the evolution of the system after suddenly changing an intrinsic parameter like, e.g., the interaction strength between the particles, is investigated. In this case, we are mainly interested in the long-time behavior; in particular, the question arises whether, relying on Boltzmann's ergodic hypothesis, relaxation to a thermal state is obtained. We find that, in general, non-thermal final states are reached, which can be described by generalized Gibbs-Boltzmann ensembles.
In a recent paper [Phys. Rev. B 59, 9699 (1999)], Chandross and Hicks claim to present a new density-matrix renormalization group (DMRG) method for dealing with excited-states of quantum lattice models. The proposed improvement to the DMRG---the inclusion of excited-state wave functions in addition to the ground state in the density matrix when calculating excitations---is in fact standard practice, is clearly stated in White's original papers, and has been used repeatedly by many groups to study excited states. The authors apply the method to the extended, dimerized Hubbard model for conjugated polymers. The criteria for determining whether states are bound or not are assessed. The authors claim that their results show that the optically important ‘‘1Bu′′ state is bound (excitonic), in contrast to a previous study. However, the discussion is qualitative, and the authors arrive at conclusions on the basis of results for one lattice size only. We show that when the criterion of Chandross and Hicks is developed into a quantitative definition of particle-hole separation, with the finite-size dependence analyzed, the implication is that the 1Bu state is unbound, at least in the sense of the density-density correlation function, in keeping with the conclusions of a previous study.