Deciphering the Influence of Ground-State Distributions on the Calculation of Photolysis Observables

Abstract

Nonadiabatic molecular dynamics offers a powerful tool for studying the photochemistry of molecular systems. Key to any nonadiabatic molecular dynamics simulation is the definition of its initial conditions (ICs), ideally representing the initial molecular quantum state of the system of interest. In this work, we provide a detailed analysis of how ICs may influence the calculation of experimental observables by focusing on the photochemistry of methylhydroperoxide (MHP), the simplest and most abundant organic peroxide in our atmosphere. We investigate the outcome of trajectory surface hopping simulations for distinct sets of ICs sampled from different approximate quantum distributions, namely harmonic Wigner functions and ab initio molecular dynamics using a quantum thermostat (QT). Calculating photoabsorption cross-sections, quantum yields, and translational kinetic energy maps from the results of these simulations reveals the significant effect of the ICs, in particular when low-frequency (∼ a few hundred cm–1) normal modes are connected to the photophysics of the molecule. Overall, our results indicate that sampling ICs from ab initio molecular dynamics using a QT is preferable for flexible molecules with photoactive low-frequency modes. From a photochemical perspective, our nonadiabatic dynamics simulations offer an explanation for a low-energy tail observed at high excitation energy in the translational kinetic energy map of MHP.

A Theoretical Perspective on the Photochemistry of Boron–Nitrogen Lewis Adducts

“…The combination of ADC(2) with the aug-cc-pVDZ basis set yields vertical transition energies and oscillator strengths in excellent agreement with EOM-CCSD/aug-cc-pVQZ for the nine lowest excited states investigated (Table S2). These results are further validated by XMS-CASPT2 calculations, even though some oscillator strengths appear weaker with this methoda possible weakness of XMS-CASPT2 already highlighted in earlier work (see, for example, ref ). Based on this first validation of ADC(2)/aug-cc-pVDZ in the FC region, we now extend our benchmark to the investigation of the photodissociation pathways expected for H 3 N–BH 3 .…”

Deconvolution of the X-ray absorption spectrum of trans-1,3-butadiene with resonant Auger spectroscopy

“…38 and 39, and the method has been successfully used in numerous cases. [40][41][42][43][44][45][46] The ground state MD simulations were performed at the B3LYP/6-31+g* level using the quantum thermostat, 47,48 with the temperature set to 300 K. The total length of the simulation was 50 ps and the time step was set to 0.5 fs. Within the NEA approach, a larger set of 400 and a smaller set of 50 geometries were selected from the ground state MD.…”

Advanced theoretical design of light-driven molecular rotary motors: enhancing thermal helix inversion and visible-light activation

“…However, more precise methods such as the quantum thermostat technique have been devised recently, and comparisons between these methods and experimental data are being made. 46 Although the molecular system in the current study includes low-frequency modes, the purpose of the TSH simulations using OM2/MRCI is qualitative rather than quanti-tative, so we used conventional Wigner sampling to set the initial conditions. To enhance accuracy, the empirical decoherence correction proposed by Granucci et al 47 (0.1 a.u.)…”