We consider the application of the original Meyer-Miller (MM) Hamiltonian to mapping fermionic quantum dynamics to classical equations of motion. Non-interacting fermionic and bosonic systems share the same one-body density dynamics when evolving fro
m the same initial many-body state. The MM classical mapping is exact for non-interacting bosons, and therefore it yields the exact time-dependent one-body density for non-interacting fermions as well. Starting from this observation, the MM mapping is compared to different mappings specific for fermionic systems, namely the spin mapping (SM) with and without including a Jordan-Wigner transformation, and the Li-Miller mapping (LMM). For non-interacting systems, the inclusion of fermionic anti-symmetry through the Jordan-Wigner transform does not lead to any improvement in the performance of the mappings and instead it worsens the classical description. For an interacting impurity model and for models of excitonic energy transfer, the MM and LMM mappings perform similarly, and in some cases the former outperforms the latter when compared to a full quantum description. The classical mappings are able to capture interference effects, both constructive and destructive, that originate from equivalent energy transfer pathways in the models.
Nonadiabatic molecular dynamics occur in a wide range of chemical reactions and femtochemistry experiments involving electronically excited states. These dynamics are hard to treat numerically as the systems complexity increases and it is thus desira
ble to have accurate yet affordable methods for their simulation. Here, we introduce a linearized semiclassical method, the generalized discrete truncated Wigner approximation (GDTWA), which is well-established in the context of quantum spin lattice systems, into the arena of chemical nonadiabatic systems. In contrast to traditional continuous mapping approaches, e.g. the Meyer-Miller-Stock-Thoss and the spin mappings, GDTWA samples the electron degrees of freedom in a discrete phase space, and thus forbids an unphysical unbounded growth of electronic state populations. The discrete sampling also accounts for an effective reduced but non-vanishing zero-point energy without an explicit parameter, which makes it possible to treat the identity operator and other operators on an equal footing. As numerical benchmarks on two Linear Vibronic Coupling models show, GDTWA has a satisfactory accuracy in a wide parameter regime, independently of whether the dynamics is dominated by relaxation or by coherent interactions. Our results suggest that the method can be very adequate to treat challenging nonadiabatic dynamics problems in chemistry and related fields.
We study the generation of electronic ring currents in the presence of nonadiabatic coupling using circularly polarized light. For this, we introduce a solvable model consisting of an electron and a nucleus rotating around a common center and subject
to their mutual Coulomb interaction. The simplicity of the model brings to the forefront the non-trivial properties of electronic ring currents in the presence of coupling to the nuclear coordinates and enables the characterization of various limiting situations transparently. Employing this model, we show that vibronic coupling effects play a crucial role even when a single $E$ degenerate eigenstate of the system supports the current. The maximum current of a degenerate eigenstate depends on the strength of the nonadiabatic interactions. In the limit of large nuclear to electronic masses, in which the Born-Oppenheimer approximation becomes exact, constant ring currents and time-averaged oscillatory currents necessarily vanish.
The effect of nuclear dynamics and conical intersections on electronic coherences is investigated employing a two-state, two-mode linear vibronic coupling model. Exact quantum dynamical calculations are performed using the multi-configuration time-de
pendent Hartree method (MCTDH). It is found that the presence of a non-adiabatic coupling close to the Franck-Condon point can preserve electronic coherence to some extent. Additionally, the possibility of steering the nuclear wavepackets by imprinting a relative phase between the electronic states during the photoionization process is discussed. It is found that the steering of nuclear wavepackets is possible given that a coherent electronic wavepacket embodying the phase difference passes through a conical intersection. A conical intersection close to the Franck-Condon point is thus a necessary prerequisite for control, providing a clear path towards attochemistry.
Nonequilibrium dynamical mean-field theory (DMFT) solves correlated lattice models by obtaining their local correlation functions from an effective model consisting of a single impurity in a self-consistently determined bath. The recently developed m
apping of this impurity problem from the Keldysh time contour onto a time-dependent single-impurity Anderson model (SIAM) [C. Gramsch et al., Phys. Rev. B 88, 235106 (2013)] allows one to use wave function-based methods in the context of nonequilibrium DMFT. Within this mapping, long times in the DMFT simulation become accessible by an increasing number of bath orbitals, which requires efficient representations of the time-dependent SIAM wave function. These can be achieved by the multiconfiguration time-dependent Hartree (MCTDH) method and its multi-layer extensions. We find that MCTDH outperforms exact diagonalization for large baths in which the latter approach is still within reach and allows for the calculation of SIAMs beyond the system size accessible by exact diagonalization. Moreover, we illustrate the computation of the self-consistent two-time impurity Greens function within the MCTDH second quantization representation.
We develop the multi-layer multi-configuration time-dependent Hartree method for bosons (ML-MCTDHB), a variational numerically exact ab-initio method for studying the quantum dynamics and stationary properties of bosonic systems. ML-MCTDHB takes adva
ntage of the permutation symmetry of identical bosons, which allows for investigations of the quantum dynamics from few to many-body systems. Moreover, the multi-layer feature enables ML-MCTDHB to describe mixed bosonic systems consisting of arbitrary many species. Multi-dimensional as well as mixed-dimensional systems can be accurately and efficiently simulated via the multi-layer expansion scheme. We provide a detailed account of the underlying theory and the corresponding implementation. We also demonstrate the superior performance by applying the method to the tunneling dynamics of bosonic ensembles in a one-dimensional double well potential, where a single-species bosonic ensemble of various correlation strengths and a weakly interacting two-species bosonic ensemble are considered.
