Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Numerical study of the nonequilibrium dynamics of 1-D electron-phonon systems using a local basis optimization

Christoph Brockt

Published Jan 1, 2018• 2 citations
Publisher

Abstract

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.

DOI
10.15488/3477

Related papers

Matrix-product-state method with a dynamical local basis optimization for bosonic systems out of equilibrium

Christoph Brockt, Florian Dorfner, Lev Vidmar +2 more

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.

arXiv:1508.00694•Physical Review B•Dec 9, 2015•67 citations

Transport and dynamical properties of electron-phonon coupled systems

D. J. Jansen

In this thesis, I investigate the transport and dynamical properties of electron-phonon coupled systems using numerical methods. I focus on the one-dimensional Holstein model, which incorporates both polaron and charge density wave (CDW) physics, and use matrix-product-state-based algorithms. In the first part, I present work where we combine purification with local basis optimization to study the finite-temperature properties of the model. This gives access to polaron spectral functions, polaron and bipolaron optical conductivity, and the energy-transport coefficient at finite filling. The second part focuses on the Holstein model coupled to two tight-binding leads. There, we investigate how current is transported through the structure and how CDW states break down. In the last part, we analyze the spread of an electron in a lattice and the decay of CDWs. We aim to understand if a full quantum mechanical treatment of the phonons is necessary or if trajectory-based algorithms are applicable.

•Jan 1, 2023•0 citations

Scattering of an electronic wave packet by a one-dimensional electron-phonon-coupled structure

Christoph Brockt, Eric Jeckelmann

We investigate the scattering of an electron by phonons in a small structure between two one-dimensional tight-binding leads. This model mimics the quantum electron transport through atomic wires or molecular junctions coupled to metallic leads. The electron-phonon-coupled structure is represented by the Holstein model. We observe permanent energy transfer from the electron to the phonon system (dissipation), transient self-trapping of the electron in the electron-phonon-coupled structure (due to polaron formation and multiple reflections at the structure edges), and transmission resonances that depend strongly on the strength of the electron-phonon coupling and the adiabaticity ratio. A recently developed TEBD algorithm, optimized for bosonic degrees of freedom, is used to simulate the quantum dynamics of a wave packet launched against the electron-phonon-coupled structure. Exact results are calculated for a single electron-phonon site using scattering theory and analytical approximations are obtained for limiting cases.

arXiv:1612.04665•Physical review. B./Physical review. B•Feb 23, 2017•9 citations

Charge-density-wave melting in the one-dimensional Holstein model

Jan Stolpp, J. Herbrych, Florian Dorfner +2 more

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 ${M}_{\mathrm{ph}}=40$ bare phonons per site by working with $O(10)$ optimal phonon modes.

arXiv:1911.01718•Physical review. B./Physical review. B•Jan 17, 2020•44 citations

Nonequilibrium electron-lattice dynamics

Sharareh Sayyad

Optimizing various properties of materials have been under intensive focus for many decades. In recent years, the exposure of ultrafast techniques enhance the possibility of seizing this goal. As by employing these techniques, one can expose the system under strong short fields, and investigate the transient response of various low-lying excitations. Such studies will strengthen our understanding of interplays among various degrees of freedom in correlated materials, and thereby facilitate the chance to manipulate different properties of these systems. Of our particular attention is the electron-lattice dynamics, where modulating a particular vibrational mode of the crystal lattice, will lead to the emergence of new quasi stable'' lattice structures as well as novel transient phenomena. Aside from these experimental prospects, the above mentioned goals were also the subject of an extensive theoretical studies. In this thesis, we investigated the in and out of equilibrium dynamics of electronic systems which are interacting with lattice degrees of freedom. We employed the two-time Keldysh Green's function within the framework of the dynamical mean field theory to conduct our researches on this subject. In our first attempt, we studied the dynamical response of a one-electron system under a linear electron-lattice interaction. We explored the real-time formation of a bound quasiparticle, known as a polaron, under various parameter regimes of the model. We presented that when the carrier is suddenly coupled to the lattice degrees of freedom, the long-time response of the system in the strong coupling limit exhibits both bound, polaronic, and unbound, delocalized, states. We characterized the nature of such a mixed state using the adiabatic picture of the model. We also showed that the formed polaron is dressed by excited phononic states, in contrast to its equilibrium counterpart, where the ground-state phononic states are bounding the electron. We also addressed the response of a many-body electronic system, weakly coupled to a bosonic bath, under various time dependent excitations, which brings the system to a bad-metallic phase. We employed a model which treated the Coulomb repulsive interaction locally. We presented that the relaxation of this excited state varies by the strength of the electron-electron interaction. We demonstrated that, for both insulating and metallic regimes of the system, a slow relaxation dynamics can be observed close to the phase transition from the metallic to insulating regimes of the system. Our analysis revealed that in the metallic side, this slow evolution is rooted in the spin-related physics of the problem. However, in the insulating regimes of the model, our results showed that presence of a gap in the spectral density of the insulating system prohibits the relaxation of the excited particles. We also studied the associated timescale for the adiabatic dynamics near the Mott transition. Our results exhibited that there is a nonadiabatic window'' in the metallic side of the phase transition, in which the adiabatic timescale is extremely long. We pointed out that the presence of such a window is intertwined with the slow relaxation timescale of the spin-related dynamics close to the Mott transition. We also captured an induced metallic phase near metal-to-insulator transition out of equilibrium. We initiated this study using the slave-rotor impurity solver, which is been coupled to the self-energies obtained from the first-order Luttinger-Ward functional for the Phonon degrees of freedom. We presented the emergence of the larger spectral density at small frequencies both in metallic and insulating phases close to the phase transition, which is surrounded by pseudo gaps. We explored the physics of such spectral densities as well as their dependencies on the phononic parameters. We postulated that this response of the system is due to the evolution of the system, in a non-thermal mechanism, towards a temperature which is less than the initial temperature of the system, before coupling to the bath. We also pointed out that the associated very slow relaxation dynamics is the result of the formation of psuedo-gaps around the quasiparticle peak.

•Jan 17, 2017•0 citations

Hybrid quantum-classical matrix-product state and Lanczos methods for electron-phonon systems with strong electronic correlations: Application to disordered systems coupled to Einstein phonons

Heiko Georg Menzler, Suman Mondal, Fabian Heidrich-Meisner

We present two quantum-classical hybrid methods for simulating the time-dependence of electron-phonon systems that treat electronic correlations numerically exactly and optical-phonon degrees of freedom classically. These are a time-dependent Lanczos and a matrix-product state method, each combined with the multi-trajectory Ehrenfest approach. Due to the approximations, reliable results are expected for the adiabatic regime of small phonon frequencies. We discuss the convergence properties of both methods for a system of interacting spinless fermions in one dimension and provide a benchmark for the Holstein chain. As a first application, we study the decay of charge density wave order in a system of interacting spinless fermions coupled to Einstein oscillators and in the presence of quenched disorder. We investigate the dependence of the relaxation dynamics on the electron-phonon coupling strength and provide numerical evidence that the coupling of strongly disordered systems to classical oscillators leads to delocalization, thus destabilizing the (finite-size) many-body localization in this system.

cond-mat.str-elCross-listedarXiv:2512.10899•Physical review. B./Physical review. B•Dec 12, 2025•3 citations
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.
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 thesis, I investigate the transport and dynamical properties of electron-phonon coupled systems using numerical methods. I focus on the one-dimensional Holstein model, which incorporates both polaron and charge density wave (CDW) physics, and use matrix-product-state-based algorithms. In the first part, I present work where we combine purification with local basis optimization to study the finite-temperature properties of the model. This gives access to polaron spectral functions, polaron and bipolaron optical conductivity, and the energy-transport coefficient at finite filling. The second part focuses on the Holstein model coupled to two tight-binding leads. There, we investigate how current is transported through the structure and how CDW states break down. In the last part, we analyze the spread of an electron in a lattice and the decay of CDWs. We aim to understand if a full quantum mechanical treatment of the phonons is necessary or if trajectory-based algorithms are applicable.
We investigate the scattering of an electron by phonons in a small structure between two one-dimensional tight-binding leads. This model mimics the quantum electron transport through atomic wires or molecular junctions coupled to metallic leads. The electron-phonon-coupled structure is represented by the Holstein model. We observe permanent energy transfer from the electron to the phonon system (dissipation), transient self-trapping of the electron in the electron-phonon-coupled structure (due to polaron formation and multiple reflections at the structure edges), and transmission resonances that depend strongly on the strength of the electron-phonon coupling and the adiabaticity ratio. A recently developed TEBD algorithm, optimized for bosonic degrees of freedom, is used to simulate the quantum dynamics of a wave packet launched against the electron-phonon-coupled structure. Exact results are calculated for a single electron-phonon site using scattering theory and analytical approximations are obtained for limiting cases.
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{M}_{\mathrm{ph}}=40Mph​=40 bare phonons per site by working with O(10)O(10)O(10) optimal phonon modes.
Optimizing various properties of materials have been under intensive focus for many decades. In recent years, the exposure of ultrafast techniques enhance the possibility of seizing this goal. As by employing these techniques, one can expose the system under strong short fields, and investigate the transient response of various low-lying excitations. Such studies will strengthen our understanding of interplays among various degrees of freedom in correlated materials, and thereby facilitate the chance to manipulate different properties of these systems. Of our particular attention is the electron-lattice dynamics, where modulating a particular vibrational mode of the crystal lattice, will lead to the emergence of new quasi stable'' lattice structures as well as novel transient phenomena. Aside from these experimental prospects, the above mentioned goals were also the subject of an extensive theoretical studies. In this thesis, we investigated the in and out of equilibrium dynamics of electronic systems which are interacting with lattice degrees of freedom. We employed the two-time Keldysh Green's function within the framework of the dynamical mean field theory to conduct our researches on this subject. In our first attempt, we studied the dynamical response of a one-electron system under a linear electron-lattice interaction. We explored the real-time formation of a bound quasiparticle, known as a polaron, under various parameter regimes of the model. We presented that when the carrier is suddenly coupled to the lattice degrees of freedom, the long-time response of the system in the strong coupling limit exhibits both bound, polaronic, and unbound, delocalized, states. We characterized the nature of such a mixed state using the adiabatic picture of the model. We also showed that the formed polaron is dressed by excited phononic states, in contrast to its equilibrium counterpart, where the ground-state phononic states are bounding the electron. We also addressed the response of a many-body electronic system, weakly coupled to a bosonic bath, under various time dependent excitations, which brings the system to a bad-metallic phase. We employed a model which treated the Coulomb repulsive interaction locally. We presented that the relaxation of this excited state varies by the strength of the electron-electron interaction. We demonstrated that, for both insulating and metallic regimes of the system, a slow relaxation dynamics can be observed close to the phase transition from the metallic to insulating regimes of the system. Our analysis revealed that in the metallic side, this slow evolution is rooted in the spin-related physics of the problem. However, in the insulating regimes of the model, our results showed that presence of a gap in the spectral density of the insulating system prohibits the relaxation of the excited particles. We also studied the associated timescale for the adiabatic dynamics near the Mott transition. Our results exhibited that there is a nonadiabatic window'' in the metallic side of the phase transition, in which the adiabatic timescale is extremely long. We pointed out that the presence of such a window is intertwined with the slow relaxation timescale of the spin-related dynamics close to the Mott transition. We also captured an induced metallic phase near metal-to-insulator transition out of equilibrium. We initiated this study using the slave-rotor impurity solver, which is been coupled to the self-energies obtained from the first-order Luttinger-Ward functional for the Phonon degrees of freedom. We presented the emergence of the larger spectral density at small frequencies both in metallic and insulating phases close to the phase transition, which is surrounded by pseudo gaps. We explored the physics of such spectral densities as well as their dependencies on the phononic parameters. We postulated that this response of the system is due to the evolution of the system, in a non-thermal mechanism, towards a temperature which is less than the initial temperature of the system, before coupling to the bath. We also pointed out that the associated very slow relaxation dynamics is the result of the formation of psuedo-gaps around the quasiparticle peak.
We present two quantum-classical hybrid methods for simulating the time-dependence of electron-phonon systems that treat electronic correlations numerically exactly and optical-phonon degrees of freedom classically. These are a time-dependent Lanczos and a matrix-product state method, each combined with the multi-trajectory Ehrenfest approach. Due to the approximations, reliable results are expected for the adiabatic regime of small phonon frequencies. We discuss the convergence properties of both methods for a system of interacting spinless fermions in one dimension and provide a benchmark for the Holstein chain. As a first application, we study the decay of charge density wave order in a system of interacting spinless fermions coupled to Einstein oscillators and in the presence of quenched disorder. We investigate the dependence of the relaxation dynamics on the electron-phonon coupling strength and provide numerical evidence that the coupling of strongly disordered systems to classical oscillators leads to delocalization, thus destabilizing the (finite-size) many-body localization in this system.