Perturbations in vibrational diatomic spectra: Factorization of the molecular wave function

On the numerical solution of the exact factorization equations

The exact factorization (EF) approach to coupled electron-ion dynamics recasts the time-dependent molecular Schrödinger equation as two coupled equations, one for the nuclear wavefunction and one for the conditional electronic wavefunction. The potentials appearing in these equations have provided insight into non-adiabatic processes, and new practical non-adiabatic dynamics methods have been formulated starting from these equations. Here, we provide a first demonstration of a self-consistent solution of the exact equations, with a preliminary analysis of their stability and convergence properties. The equations have an unprecedented mathematical form, involving a Hamiltonian outside the class of Hermitian Hamiltonians usually encountered in time-propagation, and so the usual numerical methods for time-dependent Schrödinger fail when applied in a straightforward way to the EF equations. We find an approach that enables stable propagation long enough to witness non-adiabatic behavior in a model system before non-trivial instabilities take over. Implications for the development and analysis of EF-based methods are discussed.

Nodeless vibrational amplitudes and quantum nonadiabatic dynamics in the nested funnel for a pseudo Jahn-Teller molecule or homodimer

The nonadiabatic states and dynamics are investigated for a linear vibronic coupling Hamiltonian with a static electronic splitting and weak off-diagonal Jahn-Teller coupling through a single vibration with a vibrational-electronic resonance. With a transformation of the electronic basis, this Hamiltonian is also applicable to the anti-correlated vibration in a symmetric homodimer with marginally strong constant off-diagonal coupling, where the non-adiabatic states and dynamics model electronic excitation energy transfer or self-exchange electron transfer. For parameters modeling a free-base naphthalocyanine, the nonadiabatic couplings are deeply quantum mechanical and depend on wavepacket width; scalar couplings are as important as the derivative couplings that are usually interpreted to depend on vibrational velocity in semiclassical curve crossing or surface hopping theories. A colored visualization scheme that fully characterizes the non-adiabatic states using the exact factorization is developed. The nonadiabatic states in this nested funnel have nodeless vibrational factors with strongly avoided zeroes in their vibrational probability densities. Vibronic dynamics are visualized through the vibrational coordinate dependent density of the time-dependent dipole moment in free induction decay. Vibrational motion is amplified by the nonadiabatic couplings, with asymmetric and anisotropic motions that depend upon the excitation polarization in the molecular frame and can be reversed by a change in polarization. This generates a vibrational quantum beat anisotropy in excess of 2/5. The amplitude of vibrational motion can be larger than that on the uncoupled potentials, and the electronic population transfer is maximized within one vibrational period. Most of these dynamics are missed by the adiabatic approximation, and some electronic and vibrational motions are completely suppressed by the Condon approximation of a coordinate-independent transition dipole between adiabatic states. For all initial conditions investigated, the initial nonadiabatic electronic motion is driven towards the lower adiabatic state, and criteria for this directed motion are discussed.

On the Dynamics through a Conical Intersection

Conical intersections represent critical topological features of potential energy surfaces and open ultrafast nonradiative deactivation channels for photoexcited molecules. In the following, we investigate how this funneling picture is transposed in the eyes of the exact factorization formalism for a 2D model system. The exact factorization of the total molecular wave function leads to the fundamental concept of time-dependent potential energy surface and time-dependent vector potential, whose behavior during a dynamics through a conical intersection has up to now remained unexplored. Despite the fact that these quantities might be viewed as time-dependent generalizations of the adiabatic potential energy surfaces and the nonadiabatic coupling vectors, characteristic quantities appearing in the Born-Oppenheimer framework, we observe that they do not exhibit particular topological features in the region of conical intersection but still reflect the complex dynamics of the nuclear wavepacket.

Exact Potential Driving the Electron Dynamics in Enhanced Ionization of H(2)(+)

It was recently shown that the exact factorization of the electron-nuclear wave function allows the construction of a Schrödinger equation for the electronic system, in which the potential contains exactly the effect of coupling to the nuclear degrees of freedom and any external fields. Here we study the exact potential acting on the electron in charge-resonance enhanced ionization in a model one-dimensional H(2)(+) molecule. We show there can be significant differences between the exact potential and that used in the traditional quasistatic analyses, arising from nonadiabatic coupling to the nuclear system, and that these are crucial to include for accurate simulations of time-resolved ionization dynamics and predictions of the ionization yield.