Electron-nuclear entanglement in the time-dependent molecular wavefunction

Generalized Semiclassical Ehrenfest Method: A Route to Wave Function-Free Photochemistry and Nonadiabatic Dynamics with Only Potential Energies and Gradients

We reconsider recent methods by which direct dynamics calculations of electronically nonadiabatic processes can be carried out while requiring only adiabatic potential energies and their gradients. We show that these methods can be understood in terms of a new generalization of the well-known semiclassical Ehrenfest method. This is convenient because it eliminates the need to evaluate electronic wave functions and their matrix elements along the mixed quantum-classical trajectories. The new approximations and procedures enabling this advance are the curvature-driven approximation to the time-derivative coupling, the generalized semiclassical Ehrenfest method, and a new gradient correction scheme called the time-derivative matrix (TDM) scheme. When spin-orbit coupling is present, one can carry out dynamics calculations in the fully adiabatic basis using potential energies and gradients calculated without spin-orbit coupling plus the spin-orbit coupling matrix elements. Even when spin-orbit coupling is neglected, the method is useful because it allows calculations by electronic structure methods for which nonadiabatic coupling vectors are unavailable. In order to place the new considerations in context, the article starts out with a review of background material on trajectory surface hopping, the semiclassical Ehrenfest scheme, and methods for incorporating decoherence. We consider both internal conversion and intersystem crossing. We also review several examples from our group of successful applications of the curvature-driven approximation.

Differential Shannon entropies and correlation measures for Born-Oppenheimer electron-nuclear dynamics: numerical results and their analytical interpretation

We study the Born-Oppenheimer dynamics within a model for a coupled electron-nuclear motion. Differential Shannon entropies are calculated from the time-dependent probability densities of the combined system and, using single particle densities, entropies for the electronic and nuclear degrees of freedom are derived. These functions provide information on details of the wave packet motion. From the entropies, we determine the mutual information which characterizes particle correlations. This quantity is compared to other measures of electron-nuclear entanglement. Numerical results are interpreted within an analytically solvable approach, and it is documented how these functions depend on properties of the Born-Oppenheimer wave function and, in particular, how dynamical effects like wave packet focusing and dispersion influence the correlation between the particles.

Surface modeling of photocatalytic materials for water splitting

The photocatalyst surface is central to photocatalytic reactions. However, it has been a challenge to explicitly understand both the surface configuration and the structure-dependent photocatalytic properties at the atomic level. First-principles density functional theory (DFT) calculations provide a versatile method that makes up for the lack of experimental surface studies. In DFT calculations, the initial surface model greatly affects the accuracy of the calculation results. Consequently, establishing a more realistic and more reliable material surface models is undoubtedly the first step and the most important link in theoretical calculations. The aim of this Perspective is to provide a general understanding of the methods for the surface modeling of photocatalytic materials in recent years. We begin with a discussion of the basic theories applied in photocatalytic surface research, followed by an explanation of the importance of surface modeling in photocatalysis. We then elaborate on the advantages and disadvantages of the basic surface model and briefly describe the latest surface modeling methods. Finally, we evaluate the rationality of current surface modeling methods. We summarize this Perspective by prospecting the developing directions of photocatalytic surface research in the future. It is believed that a reasonable surface model should be verified by both experimental characterization and theoretical computation with negative feedback.

Quantum-classical nonadiabatic dynamics of Floquet driven systems

We develop a trajectory-based approach for excited-state molecular dynamics simulations of systems subject to an external periodic drive. We combine the exact-factorization formalism, allowing us to treat electron-nuclear systems in nonadiabatic regimes, with the Floquet formalism for time-periodic processes. The theory is developed starting with the molecular time-dependent Schrödinger equation with the inclusion of an external periodic drive that couples to the system dipole moment. With the support of the Floquet formalism, quantum dynamics is approximated by combining classical-like, trajectory-based, nuclear evolution with electronic dynamics represented in the Floquet basis. The resulting algorithm, which is an extension of the coupled-trajectory mixed quantum-classical scheme for periodically driven systems, is applied to a model study, exactly solvable, with different field intensities.

Relaxation dynamics through a conical intersection: Quantum and quantum-classical studies

We study the relaxation process through a conical intersection of a photo-excited retinal chromophore model. The analysis is based on a two-electronic-state two-dimensional Hamiltonian developed by Hahn and Stock [J. Phys. Chem. B 104 1146 (2000)] to reproduce, with a minimal model, the main features of the 11-cis to all-trans isomerization of the retinal of rhodopsin. In particular, we focus on the performance of various trajectory-based schemes to nonadiabatic dynamics, and we compare quantum-classical results to the numerically exact quantum vibronic wavepacket dynamics. The purpose of this work is to investigate, by analyzing electronic and nuclear observables, how the sampling of initial conditions for the trajectories affects the subsequent dynamics.

Quantum Trajectories for the Dynamics in the Exact Factorization Framework: A Proof-of-Principle Test

In the framework of the exact factorization of the time-dependent electron-nuclear wave function, we investigate the possibility of solving the nuclear time-dependent Schrödinger equation based on trajectories. The nuclear equation is separated in a Hamilton-Jacobi equation for the phase of the wave function, and a continuity equation for its (squared) modulus. For illustrative adiabatic and nonadiabatic one-dimensional models, we implement a procedure to follow the evolution of the nuclear density along the characteristics of the Hamilton-Jacobi equation. Those characteristics are referred to as quantum trajectories, since they are generated via ordinary differential equations similar to Hamilton's equations, but including the so-called quantum potential, and they can be used to reconstruct exactly the quantum-mechanical nuclear wave function, provided infinite initial conditions are propagated in time.