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.
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.
Using the adaptive time-dependent density-matrix renormalization group method for the 1D Hubbard model, the splitting of local perturbations into separate wave packets carrying charge and spin is observed in real time. We show the robustness of this separation beyond the low-energy Luttinger liquid theory by studying the time evolution of single particle excitations and density wave packets. A striking signature of spin-charge separation is found in 1D cold Fermi gases in a harmonic trap at the boundary between liquid and Mott-insulating phases. We give quantitative estimates for an experimental observation of spin-charge separation in an array of atomic wires.
One of the key features of one-dimensional systems is the phenomenon of spin–charge separation. We show this separation of spin and charge in real time in the 1D Hubbard model by the application of the adaptive time-dependent density-matrix renormalization group method. It is found for single particle excitations and density wave packets far beyond the low-energy Luttinger liquid theory. An experimental set-up is proposed on how to observe the spin–charge separation with ultracold Fermi gases in a harmonic trap using the presence of the Mott-insulator region. Quantitative estimates for experimental parameters are given.
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.
The goal of this thesis is to study the transport properties and real-time dynamics of quantum magnets and ultra-cold atomic gases in one spatial dimension using numerical methods. The focus will be on the discussion of diffusive versus ballistic dynamics along with a detailed analysis of characteristic velocities in ballistic regimes. For the simulation of time-dependent density profiles we use the adaptive time-dependent density matrix renormalization group (DMRG). This numerical method allows for the simulation of time-dependent wave functions close to as well as far from equilibrium in a controlled manner. The studies of one-dimensional quantum magnets are partially motivated by the experimental evidence for a highly anisotropic and for insulators comparably high thermal conductivity of certain cuprates. We use linear response theory to study transport coefficients at arbitrary temperatures by diagonalizing small systems exactly and then calculating the current-current correlation functions. As first application we discuss the spin transport in the spin-1/2 Heisenberg chain with anisotropic exchange interactions (XXZ-chain). The second application of exact diagonalization, here in combination with time-dependent DMRG, is a discussion of the transverse components of the current-current correlation function. While usually only a Zeeman field is considered in the theory of transport coefficients, we here investigate the dynamic induced by an additional transverse magnetic field. We find that in this scenario the current-current correlation function exhibits coherent oscillations. In addition a second non-trivial frequency, different from the one expected from the usual Larmor precession, emerges and is studied varying temperature and field. Finally we calculate the frequency-dependent spin and heat conductivity of dimerized spin chains in a magnetic field. Motivated by the recent experimental studies of the phase diagram of C5H12N2CuBr4 we take the dimerized chain as a minimal model that exhibits features of the low-temperature region of the observed phase diagram. As a main result, the spin and heat conductivity obtained from linear response theory are enhanced in the field-induced gapless phase. The last application in the field of one-dimensional quantum magnets is the simulation of time-dependent energy-density wave-packets close to as well as far from equilibrium using the time-dependent density renormalization group. The main results are ballistic energy dynamics independently of how far out-of-equilibrium the initial state is and a detailed understanding of the average expansion velocity. The applications in the field of ultra-cold atomic gases focus on the sudden expansion of an initially trapped gas into an empty optical lattice. This setup was recently realized in an experiment performed by U. Schneider {\it et al.} and discussed in the context of electronic transport in the two-dimensional and the three-dimensional Fermi-Hubbard model. Here we investigate the sudden expansion of three different setups: For the expansion of a spin-balanced cloud of fermions, we identify the ballistic regime, and therein investigate the average expansion velocity of the cloud. As a main result the expansion velocity is determined by a small subset of the initial condition over a wide range of parameters. For instance, the Mott-insulating phase of the Hubbard model is characterized by a constant expansion velocity independently of the strength of the interaction. In the case of spinless bosons, we study the expansion from initial states that have a fixed particle number per lattice site and a certain concentration of defects. We study the expansion velocity as a function of interaction strength and investigate whether the time-dependent momentum distribution functions indicate a dynamical quasi-condensation. The last example is the sudden expansion of a spin-polarized gas of fermions in the presence of attractive interactions. This study is motivated by current effort to experimentally detect the Fulde-Ferrell-Larkin-Ovchinnikov state. Our results for the time-dependent momentum distribution functions and the wave-function of the pair condensate suggest that the signatures of the FFLO state vanish quickly, yet a stationary form of the momentum distribution also emerges fast. The latter is shown to be determined by the initial conditions, which might eventually allow for an indirect detection of the FFLO phase.
In recent years, ultracold bosonic and fermionic atom gases in optical lattices have allowed us to realize experimentally phenomena originally predicted or observed in condensed matter physics with unprecedented control and purity. Among others, proposals have been made to demonstrate a key feature of one dimensional (1D) quantum systems, the so-called spin-charge separation, which has been quite elusive in condensed matter systems. However, as we want to discuss in this paper using the time-dependent density-matrix renormalization group (DMRG), beyond a mere demonstration of spin-charge separation cold atom gases will also shed light on previously inaccessible aspects of this phenomenon. Moreover, the previous proposals for experimental realizations will be revisited taking into account the modelling of an external magnetic field.
This thesis deals with out-of-equilibrium transport phenomena in strongly correlated low-dimensional fermionic systems, with special emphasis on π-conjugated molecular materials. The focus of this work is to study real-time dynamics of spin and charge transport in these systems in order to investigate non-equilibrium transport in single-molecule electronic and spintronic devices.
Chapter 1 describes the electronic structure and dynamics of strongly correlated fermionic systems in general, and in one-dimension, in particular. For this purpose, effective low-energy model Hamiltonians (used in this work) are discussed. Whenever applicable, approximate analytical and numerical methods commonly used in the literature to deal with these model Hamiltonians, are outlined. In the context of one-dimensional strongly correlated fermionic systems, analytical techniques like the Bethe ansatz and bosonization, and numerical procedures like exact diagonalization and DMRG, used for solving finite systems, are discussed in detail.
Chapter 2 provides an overview of the different zero-temperature (T = 0) time-dependent DMRG algorithms, which have been used to study out-of-equilibrium time-dependent phenomena in low-dimensional strongly correlated systems.
In Chapter 3 we employ the time-dependent DMRG algorithm proposed by Luo, Xiang and Wang [Phys. Rev. Lett. 91, 049701 (2003)], to study the role of dimerization and electronic correlations on the dynamics of spin-charge separation. We employ the H¨uckel and Hubbard models for our studies. We have modified the algorithm proposed by Luo et. al to overcome some of its limitations.
Chapter 4 presents a generalized adaptive time-dependent density matrix renormalization group (DMRG) scheme developed by us, called the Double Time Window Targeting (DTWT) technique, which is capable of giving accurate results with lesser computational resources than required by the existing methods. This procedure originates from the amalgamation of the features of pace keeping DMRG algorithm, first proposed by Luo et. al, [Phys.Rev. Lett. 91, 049701 (2003)], and the time-step targeting (TST) algorithm by Feiguin and White [Phys. Rev. B 72, 020404 (2005)].
In chapter 5 we apply the Double Time Window Targeting (DTWT) technique, which was discussed in the previous chapter, for studying real-time quantum dynamics of spin-charge separation in π-conjugated polymers. We employ the Pariser-Parr-Pople (PPP) model which has long-range electron-electron interactions. For investigating real-time dynamics of spin and charge transport, we inject a hole at one end of polyene chains of different lengths and study the temporal evolution of its spin and charge degrees of freedom, using the DTWT td-DMRG algorithm.
Chapter 6 we investigate the effect of terminal substituents on the dynamics of spin and charge transport in donor-acceptor substituted polyenes (D- (CH)x- A) chains, also known as push-pull polyenes. We employ long-range correlated model Hamiltonian for the D- (CH)x- A system and, real-time DMRG dynamics for…