Multi-state trajectory approach to non-adiabatic dynamics: General formalism and the active state trajectory approximation

Nonadiabatic Field: A Conceptually Novel Approach for Nonadiabatic Quantum Molecular Dynamics

Reliable trajectory-based nonadiabatic quantum dynamics methods at the atomic/molecular level are critical for the practical understanding and rational design of many important processes in real large/complex systems, where the quantum dynamical behavior of electrons and that of nuclei are coupled. The paper reports latest progress of nonadiabatic field (NaF), a conceptually novel approach for nonadiabatic quantum dynamics with independent trajectories. Substantially different from the mainstreams of Ehrenfest-like dynamics and surface hopping methods, the nuclear force in NaF involves the nonadiabatic force arising from the nonadiabatic coupling between different electronic states, in addition to the adiabatic force contributed by a single adiabatic electronic state. NaF is capable of faithfully describing the interplay between electronic and nuclear motion in a broad regime, which covers where the relevant electronic states keep coupled in a wide range or all the time and where the bifurcation characteristic of nuclear motion is essential. NaF is derived from the exact generalized phase space formulation with coordinate-momentum variables, where constraint phase space (CPS) is employed for discrete electronic-state degrees of freedom (DOFs) and infinite Wigner phase space is used for continuous nuclear DOFs. We propose efficient integrators for the equations of motion of NaF in both adiabatic and diabatic representations. Since the formalism in the CPS formulation is not unique, NaF can in principle be implemented with various phase space representations of the time correlation function (TCF) for the time-dependent property. They are applied to a suite of representative gas-phase and condensed-phase benchmark models where numerically exact results are available for comparison. It is shown that NaF is relatively insensitive to the phase space representation of the electronic TCF and will be a potential tool for practical and reliable simulations of the quantum mechanical behavior of both electronic and nuclear dynamics of nonadiabatic transition processes in real systems.

Two-oscillator mapping modification of the Poisson bracket mapping equation formulation of the quantum-classical Liouville equation

Mapping basis solutions provide efficient ways for simulating mixed quantum-classical (MQC) dynamics in complex systems by matching multiple quantum states of interest to some fictitious physical states. Recently, various MQC methods were devised such that two harmonic oscillators are employed to represent each electronic state, showing improvements over one-oscillator-based methods. Here, we introduce and analyze newly modified mapping approximations of the quantum-classical Liouville equation (QCLE) using two oscillators for each electronic state. We design two separate mapping relations that we can adopt toward simulating dynamics and computing expectation values. Through the process, two MQC methods can be constructed, one of which actually reproduces the population dynamics of the forward and backward trajectory solution of QCLE. By applying the methods to spin-boson systems with a range of parameters, we find out that the choice of mapping relations greatly affects the simulation results. We also show that further improvement is possible through using modified identity operator formulations. Our findings may be helpful in constructing improved MQC methods in the future.

Topology of quantum coherence networks in singlet fission: mapping exciton states into real space and the dislocation induced three dimensional manifolds

An understanding of the global structure of quantum coherence networks in coupled multistate systems is of great importance for the development of emerging quantum technologies such as quantum control and quantum materials design. Here, we study the topology of a quantum coherence network of a typical singlet exciton fission system by mapping the exciton states into crystal structures in real space. The defects in crystals could lead to changes in the topological structures, and also fission dynamics. In particular, we demonstrate that the dislocation induced three dimensional manifold, which differs from its lower dimensional counterparts globally, could generate exotic global structures, such as chiral spirals, and modulate singlet fission substantially. The findings may shed light on the new possibilities of engineering effective structures for fission materials.

A molecular perspective on Tully models for nonadiabatic dynamics

We present a series of standardized molecular tests for nonadiabatic dynamics, reminiscent of the one-dimensional Tully models proposed in 1990.

Initial sampling in symmetrical quasiclassical dynamics based on Li-Miller mapping Hamiltonian

A symmetrical quasiclassical (SQC) dynamics approach based on the Li-Miller (LM) mapping Hamiltonian (SQC-LM) was employed to describe nonadiabatic dynamics. In principle, the different initial sampling procedures may be applied in the SQC-LM dynamics, and the results may be dependent on different initial sampling. We provided various initial sampling approaches and checked their influence. We selected two groups of models including site-exciton models for exciton dynamics and linear vibronic coupling models for conical intersections to test the performance of SQC-LM dynamics with the different initial sampling methods. The results were examined with respect to those of the accurate multiconfigurational time-dependent Hartree (MCTDH) quantum dynamics. For both the models, the SQC-LM method more-or-less gives a reasonable description of the population dynamics, while the influence of the initial sampling approaches on the final results is noticeable. It seems that the suitable initial sampling methods should be determined by the system under study. This indicates that the combination of the SQC-LM method with a suitable sampling approach may be a potential method in the description of nonadiabatic dynamics.

Disparity in Photoexcitation Dynamics between Vertical and Lateral MoS2/WSe2 Heterojunctions: Time-Domain Simulation Emphasizes the Importance of Donor-Acceptor Interaction and Band Alignment

Two-dimensional transition metal dichalcogenides (TMDs) heterojunctions are appealing candidates for optoelectronics and photovoltaics. Using time-domain density functional theory combined with nonadiabatic (NA) molecular dynamics, we show that photoexcitation dynamics exhibit a significant difference in the vertical and lateral MoS₂/WSe₂ heterojunctions arising from the disparity in the donor-acceptor interaction and fundamental band alignment. The obtained electron transfer time scale in the vertical heterojunction shows excellent agreement with experiment. Hole transfer proceeds 1.5 times slower. The electron-hole recombination is 3 orders of magnitude longer than the charge separation, which favors solar cell applications. On the contrary, the lateral heterojunction shows no band offsets steering charge separation. The excited electron is localized at the interface that attracts holes to form an exciton-like state due to Coulomb interaction, suggesting potential applications in light-emitting devices. The coupled electron and hole wave functions increase NA coupling and the coherence time, accelerating electron-hole recombination by a factor of 3 compared with the vertical case. The atomistic studies advance our understanding of the photoinduced charge-phonon dynamics in TMDs heterojunctions.