Non-adiabatic dynamics of photoexcited cyclobutanone: Predicting structural measurements from trajectory surface hopping with XMS-CASPT2 simulations.

Over the years, theoretical calculations and scalable computer simulations have complemented ultrafast experiments, as they offer the advantage of overcoming experimental restrictions and having access to the whole dynamics. This synergy between theory and experiment promises to yield a deeper understanding of photochemical processes, offering valuable insights into the behavior of complex systems at the molecular level. However, the ability of theoretical models to predict ultrafast experimental outcomes has remained largely unexplored. In this work, we aim to predict the electron diffraction signals of an upcoming ultrafast photochemical experiment using high-level electronic structure calculations and non-adiabatic dynamics simulations. In particular, we perform trajectory surface hopping with extended multi-state complete active space with second order perturbation simulations for understanding the photodissociation of cyclobutanone (CB) upon excitation at 200 nm. Spin-orbit couplings are considered for investigating the role of triplet states. Our simulations capture the bond cleavage after ultrafast relaxation from the 3s Rydberg state, leading to the formation of the previously observed primary photoproducts: CO + cyclopropane/propene (C3 products), ketene, and ethene (C2 products). The ratio of the C3:C2 products is found to be about 1:1. Within 700 fs, the majority of trajectories transition to their electronic ground state, with a small fraction conserving the initial cyclobutanone ring structure. We found a minimal influence of triplet states during the early stages of the dynamics, with their significance increasing at later times. We simulate MeV-ultrafast electron diffraction (UED) patterns from our trajectory results, linking the observed features with specific photoproducts and the underlying structural dynamics. Our analysis reveals highly intense features in the UED signals corresponding to the photochemical processes of CB. These features offer valuable insights into the experimental monitoring of ring opening dynamics and the formation of C3 and C2 photoproducts.

H2O-HCN complex: A new potential energy surface and intermolecular rovibrational states from rigorous quantum calculations.

In this work the H2O-HCN complex is quantitatively characterized in two ways. First, we report a new rigid-monomer 5D intermolecular potential energy surface (PES) for this complex, calculated using the symmetry-adapted perturbation theory based on density functional theory method. The PES is based on 2833 ab initio points computed employing the aug-cc-pVQZ basis set, utilizing the autoPES code, which provides a site-site analytical fit with the long-range region given by perturbation theory. Next, we present the results of the quantum 5D calculations of the fully coupled intermolecular rovibrational states of the H2O-HCN complex for the total angular momentum J values of 0, 1, and 2, performed on the new PES. These calculations rely on the quantum bound-state methodology developed by us recently and applied to a variety of noncovalently bound binary molecular complexes. The vibrationally averaged ground-state geometry of H2O-HCN determined from the quantum 5D calculations agrees very well with that from the microwave spectroscopic measurements. In addition, the computed ground-state rotational transition frequencies, as well as the B and C rotational constants calculated for the ground state of the complex, are in excellent agreement with the experimental values. The assignment of the calculated intermolecular vibrational states of the H2O-HCN complex is surprisingly challenging. It turns out that only the excitations of the intermolecular stretch mode can be assigned with confidence. The coupling among the angular degrees of freedom (DOFs) of the complex is unusually strong, and as a result most of the excited intermolecular states are unassigned. On the other hand, the coupling of the radial, intermolecular stretch mode and the angular DOFs is weak, allowing straightforward assignment of the excitation of the former.

Different flavors of exact-factorization-based mixed quantum-classical methods for multistate dynamics.

The exact factorization approach has led to the development of new mixed quantum-classical methods for simulating coupled electron-ion dynamics. We compare their performance for dynamics when more than two electronic states are occupied at a given time, and analyze: (1) the use of coupled versus auxiliary trajectories in evaluating the electron-nuclear correlation terms, (2) the approximation of using these terms within surface-hopping and Ehrenfest frameworks, and (3) the relevance of the exact conditions of zero population transfer away from nonadiabatic coupling regions and total energy conservation. Dynamics through the three-state conical intersection in the uracil radical cation as well as polaritonic models in one dimension are studied.

Conformer-Specific Dissociation Dynamics in Dimethyl Methylphosphonate Radical Cation.

The dynamics of the dimethyl methylphosphonate (DMMP) radical cation after production by strong field adiabatic ionization have been investigated. Pump-probe experiments using strong field 1300 nm pulses to adiabatically ionize DMMP and a 800 nm non-ionizing probe induce coherent oscillations of the parent ion yield with a period of about 45 fs. The yields of two fragments, PO2C2H7+ and PO2CH4+, oscillate approximately out of phase with the parent ion, but with a slight phase shift relative to each other. We use electronic structure theory and nonadiabatic surface hopping dynamics to understand the underlying dynamics. The results show that while the cation oscillates on the ground state along the P=O bond stretch coordinate, the probe excites population to higher electronic states that can lead to fragments PO2C2H7+ and PO2CH4+. The computational results combined with the experimental observations indicate that the two conformers of DMMP that are populated under experimental conditions exhibit different dynamics after being excited to the higher electronic states of the cation leading to different dissociation products. These results highlight the potential usefulness of these pump-probe measurements as a tool to study conformer-specific dynamics in molecules of biological interest.

Exact-Factorization-Based Surface Hopping for Multistate Dynamics.

A surface-hopping algorithm recently derived from the exact factorization approach, SHXF [Ha et al. J. Phys. Chem. Lett. 2018, 9, 1097], introduces an additional term in the electronic equation of surface hopping that couples electronic states through the quantum momentum. This term not only provides a first-principles description of decoherence, but here we show it is crucial to accurately capture nonadiabatic dynamics when more than two states are occupied at any given time. Using a vibronic coupling model of the uracil cation, we show that the lack of this term in traditional surface-hopping methods, including those with decoherence corrections, leads to failure to predict the dynamics through a three-state intersection, while SHXF performs similarly to the multiconfiguration time-dependent Hartree quantum dynamics benchmark.

Study of the Decoherence Correction Derived from the Exact Factorization Approach for Nonadiabatic Dynamics.

We present a detailed study of the decoherence correction to surface hopping that was recently derived from the exact factorization approach. Ab initio multiple spawning calculations that use the same initial conditions and the same electronic structure method are used as a reference for three molecules: ethylene, the methaniminium cation, and fulvene, for which nonadiabatic dynamics follows a photoexcitation. A comparison with the Granucci-Persico energy-based decoherence correction and the augmented fewest-switches surface-hopping scheme shows that the three decoherence-corrected methods operate on individual trajectories in a qualitatively different way, but the results averaged over trajectories are similar for these systems.

Impulsive alignment of 4He-CH3I: A theoretical study.

We simulate the non-adiabatic laser alignment of the weakly bound 4He-CH3I complex based on a quantum mechanical wave packet calculation for a model He-CH3I interaction potential. Two different regimes are found depending on the laser intensity. At intensities typical of non-adiabatic alignment experiments, the rotational dynamics resembles that of the isolated molecule. This is attributed to the fact that after the initial prompt alignment peak the complex rapidly dissociates. The subsequent revival pattern is due to the free rotation of the molecule detached from the helium atom. It is superimposed to a flat background corresponding to ∼20% of the wave packet which remains bound, containing lower rotational excitation. At lower intensities, dissociation is avoided but the rotational excitation is not high enough to provide an efficient alignment and a broad non-regular structure is observed. Besides, the interaction of the He atom with the molecule quenches any possible alignment. These interpretations are based on the calculation of different observables related to the rotational motion. We compare our findings with recent experimental and theoretical results of non-adiabatic alignment of linear molecules solvated in helium nanodroplets or weakly interacting with one helium atom.

Helium-induced electronic transitions in photo-excited Ba+-Hen exciplexes.

The possibility for helium-induced electronic transitions in a photo-excited atom is investigated using Ba+ excited to the 6p 2P state as a prototypical example. A diabatization scheme has been designed to obtain the necessary potential energy surfaces and couplings for complexes of Ba+ with an arbitrary number of helium atoms. It involves computing new He-Ba+ electronic wave functions and expanding them in determinants of the non-interacting complex. The 6p 2P ← 6s 2S photodissociation spectrum of He⋯Ba+ calculated with this model shows very weak coupling for a single He atom. However, several electronic relaxation mechanisms are identified, which could potentially explain the expulsion of barium ions from helium nanodroplets observed experimentally upon Ba+ photoexcitation. For instance, an avoided crossing in the ring-shaped He7Ba+ structure is shown to provide an efficient pathway for fine structure relaxation. Symmetry breaking by either helium density fluctuations or vibrations can also induce efficient relaxation in these systems, e.g., bending vibrations in the linear He2Ba+ excimer. The identified relaxation mechanisms can provide insight into helium-induced non-adiabatic transitions observed in other systems.

Quantum wave-packet dynamics in spin-coupled vibronic states.

Extending the Shin-Metiu two-electron Hamiltonian, we construct a new Hamiltonian with effective singlet-triplet couplings. The Born-Oppenheimer electronic potentials and couplings are obtained for different parameters, and the laser-free dynamics is calculated with the full Hamiltonian and in the adiabatic limit. We compare the dynamics of the system using nuclear wave packets for different numbers of Born-Oppenheimer potentials and vibronic wave packets on a full 3-dimensional (two electron coordinates plus one nuclear coordinate) grid. Using strong fields, we show that it is possible to dynamically lock the spin state of the system by decoupling the singlet-triplet transition via a nonresonant dynamic Stark effect in the adiabatic limit. Although a similar spin-locking mechanism is observed in the dynamics of vibronic wave packets, multiphoton ionization cannot be neglected leading to the breakdown of the control scheme.