Accurate time propagation for the Schrodinger equation with an explicitly time-dependent Hamiltonian

Simple and General Unitarity Conserving Numerical Real-Time Propagators of the Time-Dependent Schrödinger Equation Based on Magnus Expansion

Magnus expansion (ME) provides a general way to expand the real-time propagator of a time-dependent Hamiltonian within the exponential such that the unitarity is satisfied at any order. We use this property and explicit integration of Lagrange interpolation formulas for the time-dependent Hamiltonian within each time interval and derive approximations that preserve unitarity for the differential time evolution operators of general time-dependent Hamiltonians. The resulting second-order approximation is the same as using the average of Hamiltonians for two end points of time. We identify three fourth-order approximations involving commutators of Hamiltonians at different times and also derive a sixth-order expression. A test of these approximations along with other available expressions for a two-state time-dependent Hamiltonian with sinusoidal time dependences provides information on the relative performance of these approximations and suggests that the derived expressions can serve as useful numerical tools for time evolution in time-resolved spectroscopy, quantum control, quantum sensing, real-time ab initio quantum dynamics, and open system quantum dynamics.

Nonadiabatic wave packet dynamics and predissociation resonances in sodium hydride

Vibrational wave packet dynamics provides an opportunity to explore the energy landscape and the population transfer between nonadiabatically coupled excited electronic states. Here the coupled nonadiabatic dynamics of the C¹Σ⁺ and D¹Σ⁺ states of sodium hydride (NaH) in the gas phase in the adiabatic picture is studied, using a sequence of ultra-fast laser pulses in the femtosecond region. Emergence of different population dynamics and dissociation probabilities is shown by carefully choosing the pulse wavelength, duration and time-shift between the pulses, exciting the molecule from the ground X¹Σ⁺ state via the immediate A¹Σ⁺ state. Quantum dynamics simulations were performed in the adiabatic picture, avoiding the adiabatic to diabatic transformation. Predissociation resonances, i.e. vibrational states with finite lifetimes, arise due to nonadiabatic couplings between bound and continuum states. Here accurate resonance energies and widths are computed providing further insight into the dissociation dynamics.

Equation-of-Motion Methods for the Calculation of Femtosecond Time-Resolved 4-Wave-Mixing and N-Wave-Mixing Signals

Femtosecond nonlinear spectroscopy is the main tool for the time-resolved detection of photophysical and photochemical processes. Since most systems of chemical interest are rather complex, theoretical support is indispensable for the extraction of the intrinsic system dynamics from the detected spectroscopic responses. There exist two alternative theoretical formalisms for the calculation of spectroscopic signals, the nonlinear response-function (NRF) approach and the spectroscopic equation-of-motion (EOM) approach. In the NRF formalism, the system-field interaction is assumed to be sufficiently weak and is treated in lowest-order perturbation theory for each laser pulse interacting with the sample. The conceptual alternative to the NRF method is the extraction of the spectroscopic signals from the solutions of quantum mechanical, semiclassical, or quasiclassical EOMs which govern the time evolution of the material system interacting with the radiation field of the laser pulses. The NRF formalism and its applications to a broad range of material systems and spectroscopic signals have been comprehensively reviewed in the literature. This article provides a detailed review of the suite of EOM methods, including applications to 4-wave-mixing and N-wave-mixing signals detected with weak or strong fields. Under certain circumstances, the spectroscopic EOM methods may be more efficient than the NRF method for the computation of various nonlinear spectroscopic signals.

Nano-chevron quantum dot for spin-qubit applications

We study the theoretical properties of a parabolic hBN/MoS₂/hBN heterostructure quantum dot potential generated via electrostatic gates and its interaction with a cobalt nano chevron. We demonstrate that such an example system can undergo electric dipole spin resonance for a single electron isolated to the K' valley within the MoS₂ monolayer, and such a system can achieve pi-rotation times of approximately 5.5 ns under the influence of a 20.89 GHz driving field. Our proposed system requires operating conditions easily achievable with current experimental methods and would allow for the all-electrical control of a spin-qubit within an MoS₂ device. Our results show that such a system is experimentally feasible and would have comparable properties to that of more traditional silicon based spin-qubits. Furthermore, the design of the device can be applied to other material systems beyond MoS₂ and cobalt. In theory, the proposed structure could make use of any 2D material that experiences strong proximity exchange interactions with other magnetic materials, which makes our proposed design highly general.

Simulating Electron Dynamics of Complex Molecules with Time-Dependent Complete Active Space Configuration Interaction

Time-dependent electronic structure methods are growing in popularity as tools for modeling ultrafast and/or nonlinear processes, for computing spectra, and as the electronic structure component of mean-field molecular dynamics simulations. Time-dependent configuration interaction (TD-CI) offers several advantages over the widely used real-time time-dependent density functional theory: namely, that it correctly models Rabi oscillations; it offers a spin-pure description of open-shell systems; and a hierarchy of TD-CI methods can be defined that systematically approach the exact solution of the time-dependent Schrodinger equation (TDSE). In this work, we present a novel TD-CI approach that extends TD-CI to large complete active-space configuration expansions. Such extension is enabled by use of a direct configuration interaction approach that eliminates the need to explicitly build, store, or diagonalize the Hamiltonian matrix. Graphics processing unit (GPU) acceleration enables fast solution of the TDSE even for large active spaces-up to 12 electrons in 12 orbitals (853776 determinants) in this work. A symplectic split operator propagator yields long-time norm conservation. We demonstrate the applicability of our approach by computing the response of a large molecule with a strongly correlated ground state, decacene (C₄₂H₂₄), to various pulses (δ-function, transform limited, chirped). Our simulations predict that chirped pulses can be used to induce dipole-forbidden transitions. Simulations of decacene using the 6-31G(d) basis set and a 12 electrons/12 orbitals active space took 20.1 h to propagate for 100 fs with a 1 attosecond time step on a single NVIDIA K40 GPU. Convergence with respect to time step is found to depend on the property being computed and the chosen active space.

Exponential propagators for the Schrödinger equation with a time-dependent potential

We consider the numerical integration of the Schrödinger equation with a time-dependent Hamiltonian given as the sum of the kinetic energy and a time-dependent potential. Commutator-free (CF) propagators are exponential propagators that have shown to be highly efficient for general time-dependent Hamiltonians. We propose new CF propagators that are tailored for Hamiltonians of the said structure, showing a considerably improved performance. We obtain new fourth- and sixth-order CF propagators as well as a novel sixth-order propagator that incorporates a double commutator that only depends on coordinates, so this term can be considered as cost-free. The algorithms require the computation of the action of exponentials on a vector similar to the well-known exponential midpoint propagator, and this is carried out using the Lanczos method. We illustrate the performance of the new methods on several numerical examples.

Optimal control of orientation and entanglement for two dipole–dipole coupled quantum planar rotors

Optimal control simulations are performed for orientation and entanglement of two dipole–dipole coupled identical quantum rotors.

Ultrafast dynamics induced by the interaction of molecules with electromagnetic fields: Several quantum, semiclassical, and classical approaches

Several strategies for simulating the ultrafast dynamics of molecules induced by interactions with electromagnetic fields are presented. After a brief overview of the theory of molecule-field interaction, we present several representative examples of quantum, semiclassical, and classical approaches to describe the ultrafast molecular dynamics, including the multiconfiguration time-dependent Hartree method, Bohmian dynamics, local control theory, semiclassical thawed Gaussian approximation, phase averaging, dephasing representation, molecular mechanics with proton transfer, and multipolar force fields. In addition to the general overview, some focus is given to the description of nuclear quantum effects and to the direct dynamics, in which the ab initio energies and forces acting on the nuclei are evaluated on the fly. Several practical applications, performed within the framework of the Swiss National Center of Competence in Research "Molecular Ultrafast Science and Technology," are presented: These include Bohmian dynamics description of the collision of H with H2, local control theory applied to the photoinduced ultrafast intramolecular proton transfer, semiclassical evaluation of vibrationally resolved electronic absorption, emission, photoelectron, and time-resolved stimulated emission spectra, infrared spectroscopy of H-bonding systems, and multipolar force fields applications in the condensed phase.

Efficient methods for linear Schrödinger equation in the semiclassical regime with time-dependent potential

We build efficient and unitary (hence stable) methods for the solution of the linear time-dependent Schrödinger equation with explicitly time-dependent potentials in a semiclassical regime. The Magnus–Zassenhaus schemes presented here are based on a combination of the Zassenhaus decomposition (Bader et al. 2014 Found. Comput. Math. 14 , 689–720. ( doi:10.1007/s10208-013-9182-8 )) with the Magnus expansion of the time-dependent Hamiltonian. We conclude with numerical experiments.

Numerical solutions of the Schrödinger equation with source terms or time-dependent potentials

We develop an approach to solving numerically the time-dependent Schrödinger equation when it includes source terms and time-dependent potentials. The approach is based on the generalized Crank-Nicolson method supplemented with an Euler-MacLaurin expansion for the time-integrated nonhomogeneous term. By comparing the numerical results with exact solutions of analytically solvable models, we find that the method leads to precision comparable to that of the generalized Crank-Nicolson method applied to homogeneous equations. Furthermore, the systematic increase in precision generally permits making estimates of the error.

Solution of the time-dependent Schrödinger equation using time-dependent basis functions

The time-dependent Schrödinger equation is solved by using an explicitly time-dependent basis. This approach allows efficient reflection-free time propagation of the wave function. The applicability of the method is illustrated by solving various time-dependent problems including the calculation of the above threshold ionization of a model atom and the optical absorption spectrum of a sodium dimer.