Converging toward a practical solution of the Holstein molecular crystal model

Is there a finite mobility for the one vibrational mode Holstein model? Implications from real time simulations

The question of whether there exists a finite mobility in the standard Holstein model with one vibrational mode on each site remains unclear. In this Communication, we approach this problem by employing the hierarchical equation of motion method to simulate model systems where the vibrational modes are dissipative. It is found that, as the friction becomes smaller, the charge carrier mobility increases significantly and a friction-free limit cannot be obtained. The current autocorrelation functions are also calculated for the friction-free Holstein model, and converged results cannot be obtained with an increase in the number of sites. Based on these observations, we conclude that a finite mobility cannot be defined for the standard Holstein model in the parameter regime explored in this work.

Spectral Functions of the Holstein Polaron: Exact and Approximate Solutions

It is generally accepted that the dynamical mean field theory gives a good solution of the Holstein model, but only in dimensions greater than two. Here, we show that this theory, which becomes exact in the weak coupling and in the atomic limit, provides an excellent, numerically cheap, approximate solution for the spectral function of the Holstein model in the whole range of parameters, even in one dimension. To establish this, we make a detailed comparison with the spectral functions that we obtain using the newly developed momentum-space numerically exact hierarchical equations of motion method, which yields electronic correlation functions directly in real time. We crosscheck these conclusions with our path integral quantum Monte Carlo and exact diagonalization results, as well as with the available numerically exact results from the literature.

Probability-current analysis of energy transport in open quantum systems

We introduce a probability-current analysis of excitation energy transfer between states of an open quantum system. Expressing the energy transfer through currents of excitation probability between the states in a site representation enables us to gain key insights into the energy transfer dynamics. In particular, the analysis yields direct identification of the pathways of energy transport in large networks of sites and quantifies their relative weights, as well as the respective contributions of unitary dynamics, coherence, dephasing, and relaxation and dissipation processes to the energy transfer. It thus provides much more information than studying only excitation probabilities of the states as a function of time. Our analysis is general and can be readily applied to a broad range of dynamical descriptions of open quantum system dynamics with coupling to non-Markovian or Markovian environments.

Mobility of Holstein polaron at finite temperature: an unbiased approach

We present the first unbiased results for the mobility μ of a one-dimensional Holstein polaron obtained by numerical analytic continuation combined with diagrammatic and worldline Monte Carlo methods in the thermodynamic limit. We have identified for the first time several distinct regimes in the λ-T plane including a band conduction region, incoherent metallic region, an activated hopping region, and a high-temperature saturation region. We observe that although mobilities and mean free paths at different values of λ differ by many orders of magnitude at small temperatures, their values at T larger than the bandwidth become very close to each other.

Polaronic discontinuities induced by off-diagonal coupling

In this paper, we study a form of the Holstein molecular crystal model in which the influence of lattice vibrations on the transfers of electronic excitations between neighboring sites (off-diagonal coupling) is taken into account. Using the Toyozawa Ansatz and the Lanczos algorithm, the Holstein Hamiltonian with two types of off-diagonal coupling is studied focusing on a number of analyticity issues in the ground state. For finite-sized lattices and antisymmetric coupling, a sequence of discontinuities are found in the polaron energy dispersion, the size of the ground-state phonon cloud, and the linearized von Neumann entropy used to quantify the quantum entanglement between the exciton and the phonons in the ground state. Such behavior is accompanied by a shift of the ground-state crystal momentum from zero to nonzero values as the coupling strength is increased. In the thermodynamic limit, all discontinuities associated with antisymmetric coupling vanish except the one corresponding to the initial departure of the ground-state wavevector from the Brillouin zone center. For the case of symmetric off-diagonal coupling, a smooth crossover is found to exist in all parameters regimes.

Study of the one-dimensional Holstein model with next-nearest-neighbor hopping

We present the effect of next-nearest-neighbor hopping on Holstein polarons. The energy and wavefunction of the ground state and first excited state are obtained by using a conjugate-gradient technique on the variational Hilbert space of the electron-phonon (e-ph) system. We increase the size of the variational space systematically to obtain high accuracy. With the converged wavefunction at our disposal we obtain various physical quantities and correlation functions to get a holistic insight into the physics of e-ph interaction. Our study reveals that the effect of next-nearest-neighbor hopping is most prominent at the physically interesting crossover regime of the e-ph coupling and in the adiabatic regime. We have also performed analytical strong coupling second order perturbation with a Lang-Firsov (LF) phonon basis to account for some features in the strong coupling regime.

Influence of complex exciton-phonon coupling on optical absorption and energy transfer of quantum aggregates

We present a theory that efficiently describes the quantum dynamics of an electronic excitation that is coupled to a continuous, highly structured phonon environment. Based on a stochastic approach to non-Markovian open quantum systems, we develop a dynamical framework that allows us to handle realistic systems where a fully quantum treatment is desired yet the usual approximation schemes fail. The capability of the method is demonstrated by calculating spectra and energy transfer dynamics of mesoscopic molecular aggregates, elucidating the transition from fully coherent to incoherent transfer.