We consider the quantum evolution and relaxation of an electron or hole which is coupled to a set of bath modes. In most applications the bath modes would be the vibronic coordinates but the model considered applies to any type of dynamic boson environment. The method is developed specifically for the problem of dynamic polaron formation in small nonperiodic systems. It can describe a broad group of experimental situations, including in particular electron localization in organics and polymeric materials and devices. The immediate bath is allowed to dissipate energy to a secondary bath. The bath obeys classical dynamics which puts some restriction on the range of validity of this approach. Using the density matrix formalism on a tight binding model consisting of a linear chain coupled to vibronic coordinates, we demonstrate in real time how the interaction with a dissipative bath makes the initial quantum distribution reach a steady-state population. This calculation is based on the Ehrenfest dynamics approximation. As an example we consider coupling at a single impurity site and find that for given parameters (bath coupling, site energy, and relaxation rate), the particle becomes dynamically localized in space on a particular time scale. This localized particle can be called a polaron. We define a population formation time in the same way as done in the experimental measurement. This formation time is studied as a function of the coupling strength, bandwidth, and energy dissipation rate. Energy dissipation plays a crucial role in the spatial localization process. The formation time shortens as the electron-vibration coupling increases, and as the intersite tunneling increases, but lengthens with impurity trap depth. Polaron formation is suppressed for sufficiently wide electronic bands.
We consider the quantum evolution and relaxation of an electron or hole which is coupled to a set of bath modes. In most applications the bath modes would be the vibronic coordinates but the model considered applies to any type of dynamic boson environment. The method is developed specifically for the problem of dynamic polaron formation in small nonperiodic systems. It can describe a broad group of experimental situations, including in particular electron localization in organics and polymeric materials and devices. The immediate bath is allowed to dissipate energy to a secondary bath. The bath obeys classical dynamics which puts some restriction on the range of validity of this approach. Using the density matrix formalism on a tight binding model consisting of a linear chain coupled to vibronic coordinates, we demonstrate in real time how the interaction with a dissipative bath makes the initial quantum distribution reach a steady-state population. This calculation is based on the Ehrenfest dynamics approximation. As an example we consider coupling at a single impurity site and find that for given parameters (bath coupling, site energy, and relaxation rate), the particle becomes dynamically localized in space on a particular time scale. This localized particle can be called a polaron. We define a population formation time in the same way as done in the experimental measurement. This formation time is studied as a function of the coupling strength, bandwidth, and energy dissipation rate. Energy dissipation plays a crucial role in the spatial localization process. The formation time shortens as the electron-vibration coupling increases, and as the intersite tunneling increases, but lengthens with impurity trap depth. Polaron formation is suppressed for sufficiently wide electronic bands.
We use a one-dimensional tight binding model with an impurity site characterized by electron-vibration coupling, to describe electron transfer and localization at zero temperature, aiming to examine the process of polaron formation in this system. In particular we focus on comparing a semiclassical approach that describes nuclear motion in this many vibronic-states system on the Ehrenfest dynamics level to a numerically exact fully quantum calculation based on the Bonca-Trugman method [J. Bonča and S. A. Trugman, Phys. Rev. Lett. 75, 2566 (1995)]. In both approaches, thermal relaxation in the nuclear subspace is implemented in equivalent approximate ways: In the Ehrenfest calculation the uncoupled (to the electronic subsystem) motion of the classical (harmonic) oscillator is simply damped as would be implied by coupling to a Markovian zero temperature bath. In the quantum calculation, thermal relaxation is implemented by augmenting the Liouville equation for the oscillator density matrix with kinetic terms that account for the same relaxation. In both cases we calculate the probability to trap the electron by forming a polaron and the probability that it escapes to infinity. Comparing these calculations, we find that while both result in similar long time yields for these processes, the Ehrenfest-dynamics based calculation fails to account for the correct time scale for the polaron formation. This failure results, as usual, from the fact that at the early stage of polaron formation the classical nuclear dynamics takes place on an unphysical average potential surface that reflects the distributed electronic population in the system, while the quantum calculation accounts fully for correlations between the electronic and vibrational subsystems.
We present results for the time evolution of a one-dimensional system consisting of an electron, described by a tight-binding Hamiltonian and a harmonic latttice, coupled by a deformation-type potential. We solve numerically the nonlinear system of equations of motion for this model in order to study the effects of varying the electronic effective mass for several initial conditions and coupling strengths. Different types of localized and extended states are formed with features that are absent from the traditional polaronic states and depend very strongly on the initial electronic configuration and effective mass in a very often unexpected manner. We find that, in general, an increase of the initial electronic energy decreases the ability of the system to form localized states. However, a large effective mass favors localized polaron formation for initially localized electrons, but this is not always the case for initially extended electronic states. In the latter case, increasing the effective mass of an electron initially close to the bottom of the band makes localization more difficult, while for an initially highly excited electron, localized polaron formation is possible only when the electronic effective mass and the atomic masses of the lattice become of the same order. Finally, for a small parameter range, we find an impressive recurrence, a periodic and a complete exchange between the electronic and vibrational degrees of freedom of a small part of the initial electronic energy.
We investigate the dissipative dynamics of multilevel quantum systems strongly coupled either to a lossy cavity mode or directly to a bosonic environment. By deriving spectral localization bounds, we establish conditions under which strong system-bath coupling gives rise to a hierarchy of population relaxation timescales. Our approach builds on the reaction-coordinate polaron-transform framework. By mapping the original strong-coupling problem onto an effective weakly dissipative model, we analyze the spectrum of the resulting Liouvillian superoperator through localization bounds. For the generalized V model, we find that strong system-bath coupling gives rise to a bright-dark structure in the effective system-bath coupling operator: a single collective mode remains strongly coupled to the environment, while the remaining modes become progressively dark. Consequently, the dynamics separate into fast and slow sectors and, at finite coupling strengths, develop a hierarchy of population relaxation timescales. Numerical simulations based on both secular and non-secular quantum master equations corroborate the emergence of timescale separation and the pronounced slowing down of dissipative dynamics at strong coupling. Our results reveal a general mechanism underlying anomalously slow relaxation in strongly coupled open quantum systems and provide a route for engineering long-lived states through system-environment interactions.
We investigate the effect of inter-mode coupling on the vibrational relaxation dynamics of molecules in weak dissipative environments. The simulations are performed within the reduced density matrix formalism in the Markovian regime, assuming a Lindblad form for the system-bath interaction. The prototypical two-dimensional model system representing two CO molecules approaching a Cu(100) surface is adapted from an ab initio potential, while the diatom-diatom vibrational coupling strength is systematically varied. In the weak system-bath coupling limit and at low temperatures, only first order non-adiabatic uni-modal coupling terms contribute to surface-mediated vibrational relaxation. Since dissipative dynamics is non-unitary, the choice of representation will affect the evolution of the reduced density matrix. Two alternative representations for computing the relaxation rates and the associated operators are thus compared: the fully coupled spectral basis, and a factorizable ansatz. The former is well-established and serves as a benchmark for the solution of Liouville-von Neumann equation. In the latter, a contracted grid basis of potential-optimized discrete variable representation is tailored to incorporate most of the inter-mode coupling, while the Lindblad operators are represented as tensor products of one-dimensional operators, for consistency. This procedure results in a marked reduction of the grid size and in a much more advantageous scaling of the computational cost with respect to the increase of the dimensionality of the system. The factorizable method is found to provide an accurate description of the dissipative quantum dynamics of the model system, specifically of the time evolution of the state populations and of the probability density distribution of the molecular wave packet. The influence of intra-molecular vibrational energy redistribution appears to be properly taken into account by the new model on the whole range of coupling strengths. It demontrates that most of the mode mixing during relaxation is due to the potential part of the Hamiltonian and not to the coupling among relaxation operators.
We study the quantum dynamics of small polaron formation and polaron transport through finite quantum structures in the framework of the one-dimensional Holstein model with site-dependent potentials and interactions. Combining Lanczos diagonalization with Chebyshev moment expansion of the time evolution operator, we determine how different initial states, representing stationary ground states or injected wave packets, after an electron-phonon interaction quench, develop in real space and time. Thereby, the full quantum nature and dynamics of electrons and phonons is preserved. We find that the decay out of the initial state sensitively depends on the energy and momentum of the incoming particle, the electron-phonon coupling strength, and the phonon frequency, whereupon bound polaron-phonon excited states may emerge in the strong-coupling regime. The tunneling of a Holstein polaron through a quantum wall or dot is generally accompanied by strong phonon number fluctuations due to phonon emission and reabsorption processes.
We propose an alternative concept for the dynamics of a quantum bath, the Chebyshev space, and a method based on this concept, the Chebyshev space method. The Chebyshev space is an abstract vector space that exactly represents the fermionic or bosonic bath degrees of freedom, without a discretization of the bath density of states. Relying on Chebyshev expansions, the Chebyshev space representation of a bath has very favorable properties with respect to extremely precise and efficient calculations of ground state properties, static and dynamical correlations, and time evolution for a great variety of quantum systems. The aim of the present work is to introduce the Chebyshev space in detail and to demonstrate the capabilities of the Chebyshev space method. Although the central idea is derived in full generality, the focus is on model systems coupled to fermionic baths. In particular, we address quantum impurity problems, such as an impurity in a host or a bosonic impurity with a static barrier, and the motion of a wave packet on a chain coupled to leads. For the bosonic impurity, the phase transition from a delocalized electron to a localized polaron in arbitrary dimension is detected. For the wave packet on a chain, we show how the Chebyshev space method implements different boundary conditions, including transparent boundary conditions replacing infinite leads. Furthermore, the self-consistent solution of the Holstein model in infinite dimension is calculated. With the examples, we demonstrate how highly accurate results for system energies, correlation and spectral functions, and time dependence of observables are obtained with modest computational effort.