A reciprocal-space formulation of mixed quantum-classical dynamics

Mixed quantum-classical modeling of exciton-phonon scattering in solids: Application to optical linewidths of monolayer MoS2

We present a mixed quantum-classical framework for the microscopic and non-Markovian modeling of exciton-phonon scattering in solid-state materials and apply it to calculate the optical linewidths of monolayer MoS2. Within this framework, we combine reciprocal-space mixed quantum-classical dynamics with models for the quasiparticle band structure as well as the electron-hole and carrier-phonon interactions, parametrized against ab initio calculations, although noting that a direct interfacing with ab initio calculations is straightforward in principle. We introduce various parameters for truncating the Brillouin zone to select regions of interest. Variations of these parameters allow us to determine linewidths in the limit of asymptotic material sizes. The obtained asymptotic linewidths are found to agree favorably with experimental measurements across a range of temperatures. As such, our framework establishes itself as a promising route toward unraveling the non-Markovian and microscopic principles governing the nonadiabatic dynamics of solids.

Mixed Quantum-Classical Dynamics under Arbitrary Unitary Basis Transformations

A common approach to minimizing the cost of quantum computations is by unitarily transforming a quantum system into a basis that can be optimally truncated. Here, we derive classical equations of motion subjected to similar unitary transformations and propose their integration into mixed quantum-classical dynamics, allowing this class of methods to be applied within arbitrary bases for both the quantum and classical coordinates. To this end, canonical positions and momenta of the classical degrees of freedom are combined into a set of complex-valued coordinates amenable to unitary transformations. We demonstrate the potential of the resulting approach by means of surface hopping calculations of an electronic carrier scattering onto a single impurity in the presence of phonons. Appropriate basis transformations, capturing both the localization of the impurity and the delocalization of higher-energy excitations, are shown to faithfully capture the dynamics within a fraction of the classical and quantum basis sets.

Floquet nonadiabatic mixed quantum-classical dynamics in periodically driven solid systems

In this paper, we introduce the Floquet mean-field dynamics and Floquet surface hopping approaches to study the nonadiabatic dynamics in periodically driven solid systems. We demonstrate that these two approaches can be formulated in both real and reciprocal spaces. Using the two approaches, we are able to simulate the interaction between electronic carriers and phonons under periodic drivings, such as strong light-matter interactions. Employing the Holstein and Peierls models, we show that strong light-matter interactions can effectively modulate the dynamics of electronic population and mobility. Notably, our study demonstrates the feasibility and effectiveness of modeling low-momentum carriers' interactions with phonons using a truncated reciprocal space basis, an approach impractical in real space frameworks. Moreover, we reveal that even with a significant truncation, carrier populations derived from surface hopping maintain greater accuracy compared to those obtained via mean-field dynamics. These results underscore the potential of our proposed methods in advancing the understanding of carrier-phonon interactions in various periodically driven materials.

Overcoming positivity violations for density matrices in surface hopping

Fewest-switches surface hopping (FSSH) has emerged as one of the leading methods for modeling the quantum dynamics of molecular systems. While its original formulation was limited to adiabatic populations, the growing interest in the application of FSSH to coherent phenomena prompts the question of how one should construct a complete density matrix based on FSSH trajectories. A straightforward solution is to define adiabatic coherences based on wavefunction coefficients. In this paper, we demonstrate that inconsistencies introduced in the density matrix through such treatment may lead to a violation of positivity. We furthermore show that a recently proposed coherent generalization of FSSH results in density matrices that satisfy positivity while yielding improved accuracy throughout much (but not all) of parameter space.

Surface hopping dynamics in periodic solid-state materials with a linear vibronic coupling model

We report a surface hopping approach in which the implemented linear vibronic coupling Hamiltonian is constructed and the electronic wavefunction is propagated in the reciprocal space. The parameters of the linear vibronic coupling model, including onsite energies, phonon frequencies, and electron-phonon couplings, are calculated with density-functional theory and density-functional perturbation theory and interpolated in fine sampling points of the Brillouin zone with maximally localized Wannier functions. Using this approach, we studied the relaxation dynamics of the photo-excited hot carrier in a one-dimensional periodic carbon chain. The results show that the completeness of the number of Hilbert space k points and the number of phonon q points plays an important role in the hot carrier relaxation processes. By calculating the relaxation times of hot carriers under different reciprocal space sampling and extrapolating with the stretched-compressed exponential function, the relaxation times of hot electrons and holes in the quasi-continuous energy band are obtained. By considering the feedback effect in the hopping processes and analyzing the time-dependent phonon energy in different normal modes, we found that the long-wave longitudinal optical phonons play a major role in the relaxation dynamics of hot electrons and holes. We, therefore, provided herein an efficient and accurate approach for modeling the photophysical processes in periodic solid-state material systems.

A reciprocal-space formulation of surface hopping

Surface hopping has seen great success in describing molecular phenomena where electronic excitations tend to be localized, but its application to materials with band-like electronic properties has remained limited. Here, we derive a formulation of fewest-switches surface hopping where both the quantum and classical equations of motion are solved entirely in terms of reciprocal-space coordinates. The resulting method is directly compatible with band structure calculations and allows for the efficient description of band-like phenomena by means of a truncation of the Brillouin zone. Using the Holstein and Peierls models as examples, we demonstrate the formal equivalence between real-space and reciprocal-space surface hopping and assess their accuracy against mean-field mixed quantum-classical dynamics and numerically exact results. Having very similar equations of motion, reciprocal-space surface hopping can be straightforwardly incorporated in existing (real-space) surface hopping implementations.