We present a method for simulating the time evolution of one-dimensional correlated electron-phonon systems which combines the time-evolving block decimation algorithm with a dynamical optimization of the local basis. This approach can reduce the computational cost by orders of magnitude when boson fluctuations are large. The method is demonstrated on the nonequilibrium Holstein polaron by comparison with exact simulations in a limited functional space and on the scattering of an electronic wave packet by local phonon modes. Our study of the scattering problem reveals a rich physics including transient self-trapping and dissipation.
We present a method for simulating the time evolution of one-dimensional correlated electron-phonon systems which combines the time-evolving block decimation algorithm with a dynamical optimization of the local basis. This approach can reduce the computational cost by orders of magnitude when boson fluctuations are large. The method is demonstrated on the nonequilibrium Holstein polaron by comparison with exact simulations in a limited functional space and on the scattering of an electronic wave packet by local phonon modes. Our study of the scattering problem reveals a rich physics including transient self-trapping and dissipation.
In this dissertation a newly developed numerical method is presented, which is optimized for the time evolution of one-dimensional lattice systems with large local Hilbert spaces. This method extends the time-evolving block-decimation (TEBD) to include a local basis optimization (LBO), which has already been successfully combined with ground state methods. The algorithm is based on matrix product states (MPS), which can represent the quantum state of a one-dimensional chain in most cases with a number of parameters, that is not exponentially increasing with the chain length. The LBO causes a reduction of the simulation times that is linear in the bond dimension of the MPS. To demonstrate the advantages of this method, we apply the TEBD-LBO to electron-phonon (e-p) systems. In this thesis these are described by the Holstein model, which goes beyond semi-classical approximations and covers the full quantum statistics of the phonons. The understanding of the nonequilibrium dynamics of charge carriers coupled to lattice vibrations is of great importance for research areas like transport through quasi one-dimensional conductors, photo-generated phase transitions and time-resolved spectroscopy. First, the energy transfer from a highly excited electron to the phononic degrees of freedom on a small chain is studied. In the various parameter regimes different types of relaxation occur. In any case, after a certain time the system reaches a state, where on average no energy is exchanged between the electron and phonons. This can either mean a constant kinetic energy or oscillations with a constant amplitude and frequency. Next, long perfectly conducting leads without coupling to phonons were attached on both sides of the small chain. In the left lead an electron with density distribution in the shape of a Gaussian wave packet is injected with momentum towards the e-p coupled structure in the middle. This structure acts as an impurity in an otherwise perfectly conducting chain. The investigation shows resonance effects in the transmission and reflection at this impurity. Further, the electron can transfer a part of its energy permanently to the phonons, which results in a reduction of the velocity. Finally, two mechanisms are presented that lead to self-trapping of the electron on the e-p coupled structure.
Nonequilibrium dynamics of the open chain Holstein-Hubbard model is studied using the linear time-scaling GKBA+ODE scheme developed in Pavlyukh et al. [Phys. Rev. B 105, 125134 (2022)]. We focus on the set of parameters relevant for photovoltaic materials, i.e., a pair of electrons interacting with phonons at the crossover between the adiabatic and antiadiabatic regimes and at moderately large electron-electron interaction. By comparing with exact solutions for two corner cases, we demonstrate the accuracy of the T matrix (in the pp channel) and the second-order Fan (GD) approximations for the treatment of electronic ($e\text{\ensuremath{-}}e$) and electron-phonon (e-ph) correlations, respectively. The feedback of electrons on phonons is consistently included and is shown to be mandatory for the total energy conservation. When two interactions are simultaneously present, our simulations offer a glimpse into the dynamics of doublons and polarons unveiling the formation, propagation and decay of these quasiparticles, energy redistribution between them and self-trapping of electrons.
In this thesis, different wave-function based numerical methods are used to study many-body systems with bosonic degrees of freedom in 1d. The main focus is on the Holstein model of spinless fermions, where ground-state phases, nonequilibrium dynamics, and thermalization are investigated. Using matrix-product-state methods for the investigation of such models with electron-phonon coupling poses special numerical challenges. Local basis optimization, projected purification, and subspace expansion are used to overcome these challenges.
Nonequilibrium dynamics of the open chain Holstein-Hubbard model is studied using the linear time-scaling GKBA+ODE scheme developed in the preceeding paper. We focus on the set of parameters relevant for photovoltaic materials, i.\,e., a pair of electrons interacting with phonons at the cross-over between the adiabatic and anti-adiabatic regimes and at moderately large electron-electron interaction. By comparing with exact solutions for two corner cases, we demonstrate the accuracy of the T-matrix (in the pp channel) and the second-order Fan (GD) approximations for the treatment of electronic (e-e) and electron-phonon (e-ph) correlations, respectively. The feedback of electron on phonons is consistently included and is shown to be mandatory for the total energy conservation. When two interactions are simultaneously present, our simulations offer a glimpse into the dynamics of doublons and polarons unveiling the formation, propagation and decay of these quasiparticles, energy redistribution between them and self-trapping of electrons.
We study the Holstein model of spinless fermions, which at half filling exhibits a quantum phase transition from a metallic Tomonaga-Luttinger liquid phase to an insulating charge-density-wave (CDW) phase at a critical electron-phonon coupling strength. In our work, we focus on the real-time evolution starting from two different types of initial states that are CDW ordered: (i) ideal CDW states with and without additional phonons in the system and (ii) correlated ground states in the CDW phase. We identify the mechanism for CDW melting in the ensuing real-time dynamics and show that it strongly depends on the type of initial state. We focus on the far-from-equilibrium regime and emphasize the role of electron-phonon coupling rather than dominant electronic correlations, thus complementing a previous study of photoinduced CDW melting [H. Hashimoto and S. Ishihara, Phys. Rev. B 96, 035154 (2017)]. The numerical simulations are performed by means of matrix-product-state based methods with a local basis optimization (LBO). Within these techniques, one rotates the local (bosonic) Hilbert spaces adaptively into an optimized basis that can then be truncated while still maintaining a high precision. In this work, we extend the time-evolving block decimation (TEBD) algorithm with LBO, previously applied to single-polaron dynamics, to a half-filled system. We demonstrate that in some parameter regimes, a conventional TEBD method without LBO would fail. Furthermore, we introduce and use a ground-state density-matrix renormalization group method for electron-phonon systems using local basis optimization. In our examples, we account for up to Mph=40 bare phonons per site by working with O(10) optimal phonon modes.
We adapt the time-evolving block decimation (TEBD) algorithm, originally devised to simulate the dynamics of one-dimensional quantum systems, to simulate the time evolution of nonequilibrium stochastic systems. We describe this method in detail; a system's probability distribution is represented by a matrix product state (MPS) of finite dimension and then its time evolution is efficiently simulated by repeatedly updating and approximately refactorizing this representation. We examine the use of MPS as an approximation method, looking at parallels between the interpretations of applying it to quantum state vectors and probability distributions. In the context of stochastic systems we consider two types of factorization for use in the TEBD algorithm: non-negative matrix factorization (NMF), which ensures that the approximate probability distribution is manifestly non-negative, and the singular value decomposition (SVD). Comparing these factorizations, we find the accuracy of the SVD to be substantially greater than current NMF algorithms. We then apply TEBD to simulate the totally asymmetric simple exclusion process (TASEP) for systems of up to hundreds of lattice sites in size. Using exact analytic results for the TASEP steady state, we find that TEBD reproduces this state such that the error in calculating expectation values can be made negligible even when severely compressing the description of the system by restricting the dimension of the MPS to be very small. Out of the steady state we show for specific observables that expectation values converge as the dimension of the MPS is increased to a moderate size.