Direct nonadiabatic quantum dynamics simulations of the photodissociation of phenol

Modeling Photodissociation: Quantum Dynamics Simulations of Methanol

A comprehensive computational study of the gas-phase photodissociation dynamics of methanol is presented. Using a multiconfigurational active space based method (RASSCF) to obtain multidimensional potential energy surfaces (PESs) on-the-fly, direct quantum dynamics simulations were run using the variational multi-configurational Gaussian method (DD-vMCG). Different initial excitation energies were simulated to investigate the dependence of the branching ratios on the electronic state being populated. A detailed mechanistic explanation is provided for the observed differences with respect to the excitation energy. Population of the lowest lying excited state of methanol leads to rapid hydroxyl hydrogen loss as the main dissociation channel. This is rationalized by the strongly dissociative nature of the PES cut along the O-H stretching coordinate, confirmed by the broad feature in the absorption spectrum. In contrast, more energetic excitations lead mainly to C-O bond breaking. Again, analysis of the diabatic surfaces offers a clear explanation in terms of the nature of the electronic states involved and the coupling between them. The type of calculations presented, as well as the subsequent analysis of the results, should be seen as a general workflow for the modeling of photochemical reactions.

On the multiphoton ionisation photoelectron spectra of phenol

The phenol molecule is a prototype for non-adiabatic dynamics and the excited-state photochemistry of biomolecules. In this article, we report a joint theoretical and experimental investigation on the resonance enhanced multiphoton ionisation photoelectron (REMPI) spectra of the two lowest ionisation bands of phenol. The focus is on the theoretical interpretation of the measured spectra using quantum dynamics simulations. These were performed by numerically solving the time-dependent Schrödinger equation using the multi-layer variant of the multiconfiguration time-dependent Hartree algorithm together with a vibronic coupling Hamiltonian model. The ionising laser pulse is modelled explicitly within the ionisation continuum model to simulate experimental femtosecond 1+1 REMPI photoelectron spectra. These measured spectra are sensitive to very short lived electronically excited states, providing a rigorous benchmark for our theoretical methods. The match between experiment and theory allows for an interpretation of the features of the spectra at different wavelengths and shows that there are features due to both 'direct' and 'indirect' ionisation, resulting from non-resonant and resonant excitation by the pump pulse.

Quantum Dynamics of Photodissociation: Recent Advances and Challenges

Recent advances in constructing accurate potential energy surfaces and nonadiabatic couplings from high-level ab initio data have revealed detailed potential landscapes in not only the ground electronic state but also excited ones. They enabled quantitatively accurate characterization of photoexcited reactive systems using quantum mechanical methods. In this Perspective, we survey the recent progress in quantum mechanical studies of adiabatic and nonadiabatic photodissociation dynamics, focusing on initial state control and product energy disposal. These new insights helped to understand quantum effects in small prototypical systems, and the results serve as benchmarks for developing more approximate theoretical methods.

What Controls the Quality of Photodynamical Simulations? Electronic Structure Versus Nonadiabatic Algorithm

The field of nonadiabatic dynamics has matured over the last decade with a range of algorithms and electronic structure methods available at the moment. While the community currently focuses more on developing and benchmarking new nonadiabatic dynamics algorithms, the underlying electronic structure controls the outcome of nonadiabatic simulations. Yet, the electronic-structure sensitivity analysis is typically neglected. In this work, we present a sensitivity analysis of the nonadiabatic dynamics of cyclopropanone to electronic structure methods and nonadiabatic dynamics algorithms. In particular, we compare wave function-based CASSCF, FOMO-CASCI, MS- and XMS-CASPT2, density-functional REKS, and semiempirical MRCI-OM3 electronic structure methods with the Landau-Zener surface hopping, fewest switches surface hopping, and ab initio multiple spawning with informed stochastic selection algorithms. The results clearly demonstrate that the electronic structure choice significantly influences the accuracy of nonadiabatic dynamics for cyclopropanone even when the potential energy surfaces exhibit qualitative and quantitative similarities. Thus, selecting the electronic structure solely on the basis of the mapping of potential energy surfaces can be misleading. Conversely, we observe no discernible differences in the performance of the nonadiabatic dynamics algorithms across the various methods. Based on the above results, we discuss the present-day practice in computational photodynamics.

Mode-dependent H atom tunneling dynamics of the S1 phenol is resolved by the simple topographic view of the potential energy surfaces along the conical intersection seam

Mode-dependent H atom tunneling dynamics of the O-H bond predissociation of the S₁ phenol has been theoretically analyzed. As the tunneling is governed by the complicated multi-dimensional potential energy surfaces that are dynamically shaped by the upper-lying S₁(ππ*)/S₂(πσ*) conical intersection, the mode-specific tunneling dynamics of phenol (S₁) has been quite formidable to be understood. Herein, we have examined the topography of the potential energy surface along the particular S₁ vibrational mode of interest at the nuclear configurations of the S₁ minimum and S₁/S₂ conical intersection. The effective adiabatic tunneling barrier experienced by the reactive flux at the particular S₁ vibrational mode excitation is then uniquely determined by the topographic shape of the potential energy surface extended along the conical intersection seam coordinate associated with the particular vibrational mode. The resultant multi-dimensional coupling of the specific vibrational mode to the tunneling coordinate is then reflected in the mode-dependent tunneling rate as well as nonadiabatic transition probability. Remarkably, the mode-specific experimental result of the S₁ phenol tunneling reaction [K. C. Woo and S. K. Kim, J. Phys. Chem. A 123, 1529-1537 (2019)] (in terms of the qualitative and relative mode-dependent dynamic behavior) could be well rationalized by semi-classical calculations based on the mode-specific topography of the effective tunneling barrier, providing the clear conceptual insight that the skewed potential energy surfaces along the conical intersection seam (strongly or weakly coupled to the tunneling reaction coordinate) may dictate the tunneling dynamics in the proximity of the conical intersection.

Mixed-quantum-classical or fully-quantized dynamics? A unified code to compare methods

Three methods for non-adiabatic dynamics are compared to highlight their capabilities. Multi-configurational time-dependent Hartree is a full grid-based solution to the time-dependent Schrödinger equation, variational multi-configurational Gaussian (vMCG) uses a less flexible but unrestricted Gaussian wavepacket basis, and trajectory surface hopping (TSH) replaces the nuclear wavepacket with a swarm of classical trajectories. Calculations with all methods using a model Hamiltonian were performed. The vMCG and TSH were also then run in a direct dynamics mode, with the potential energy surfaces calculated on-the-fly using quantum chemistry calculations. All dynamics calculations used the Quantics package, with the TSH calculations using a new interface to a surface hopping code. A novel approach to calculate adiabatic populations from grid-based quantum dynamics using a time-dependent discrete variable representation is presented, allowing a proper comparison of methods. This article is part of the theme issue 'Chemistry without the Born-Oppenheimer approximation'.

Diabatic States of Molecules

Quantitative simulations of electronically nonadiabatic molecular processes require both accurate dynamics algorithms and accurate electronic structure information. Direct semiclassical nonadiabatic dynamics is expensive due to the high cost of electronic structure calculations, and hence it is limited to small systems, limited ensemble averaging, ultrafast processes, and/or electronic structure methods that are only semiquantitatively accurate. The cost of dynamics calculations can be made manageable if analytic fits are made to the electronic structure data, and such fits are most conveniently carried out in a diabatic representation because the surfaces are smooth and the couplings between states are smooth scalar functions. Diabatic representations, unlike the adiabatic ones produced by most electronic structure methods, are not unique, and finding suitable diabatic representations often involves time-consuming nonsystematic diabatization steps. The biggest drawback of using diabatic bases is that it can require large amounts of effort to perform a globally consistent diabatization, and one of our goals has been to develop methods to do this efficiently and automatically. In this Feature Article, we introduce the mathematical framework of diabatic representations, and we discuss diabatization methods, including adiabatic-to-diabatic transformations and recent progress toward the goal of automatization.

Semiclassical Trajectory Studies of Reactive and Nonreactive Scattering of OH(A 2S+) by H2 Based on an Improved Full-Dimensional Ab Initio Diabatic Potential Energy Matrix

We present a new full-dimensional diabatic potential energy matrix (DPEM) for electronically nonadiabatic collisions of OH( A 2 Σ + ) with H 2 , and we calculate the probabilities of electronically adiabatic inelastic collisions, nonreactive quenching, and reactive quenching to form H 2 O + H. The DPEM was fitted using a many-body expansion with permutationally invariant polynomials in bond-order functions to represent the many-body part. The dynamics calculations were carried out with the fewest-switches with time uncertainty and stochastic decoherence (FSTU/SD) semiclassical trajectory method. We present results both for head-on collisions (impact parameter b equal to zero) and for a full range of impact parameters. The results are compared to experiment and to earlier FSTU/SD and quantum dynamics calculations with a previously published DPEM. The various theoretical results all agree that nonreactive quenching dominates reactive quenching, but there are quantitative differences between the two DPEMs and between the b = 0 results and the all- b results, especially for the probability of reactive quenching.