That Little Extra Kick: Nonadiabatic Effects in Acetaldehyde Photodissociation

Ultrafast photochemistry and electron diffraction for cyclobutanone in the S2 state: Surface hopping with time-dependent density functional theory

We simulate the photodynamics of gas-phase cyclobutanone excited to the S2 state using fewest switches surface hopping (FSSH) dynamics powered by time-dependent density functional theory (TDDFT). We predict a total photoproduct yield of 8%, with a C3:C2 product ratio of 0 trajectories to 8 trajectories. One primary S2 → S1 conical intersection is identified involving the compression of an α-carbon-carbon-hydrogen bond angle. Excited state lifetimes computed with respect to electronic state populations were found to be 3.96 ps (S2 → S1) and 498 fs (S1 → S0). We also generate time-resolved difference pair distribution functions (ΔPDFs) from our TDDFT-FSSH dynamics results in order to generate direct comparisons with ultrafast electron diffraction experiment observables. Global and target analysis of time-resolved ΔPDFs produced a distinct set of lifetimes: (i) a 0.548 ps decay and (ii) a 1.69 ps decay, both resembling the S2 minimum, as well as (iii) a long decay that resembles the S1 minimum geometry and the fully separated C2 products. Finally, we contextualize our results by considering the impact of the most likely sources of significant errors.

A high-level ab initio study of the photodissociation of acetaldehyde

Acetaldehyde is a very relevant atmospheric species whose photodissociation has been extensively studied in the first absorption band both experimentally and theoretically. Very few works have been reported on acetaldehyde photodissociation at higher excitation energies. In this work, the photodissociation dynamics of acetaldehyde is investigated by means of high-level multireference configuration interaction ab initio calculations. Five different fragmentation pathways of acetaldehyde are explored by calculating the potential-energy curves of the ground and several excited electronic states along the corresponding dissociating bond distances. The excitation energy range covered in the study is up to 10 eV, nearly the ionization energy of acetaldehyde. We intend to rationalize the available experimental results and, in particular, to elucidate why some of the studied fragmentation pathways are experimentally observed in the different excitation energy regions and some others are not. Based on the shape of the calculated potential curves, we are able to explain the main findings of the available experiments, also suggesting possible dynamical dissociation mechanisms in the different energy regions. Thus, the reported potential curves are envisioned as a useful tool to interpret the currently available experiments as well as future ones on acetaldehyde photodissociation at excitation wavelengths in the range studied here.

Parametrically Managed Activation Functions for Improved Global Potential Energy Surfaces for Six Coupled 5A' States and Fourteen Coupled 3A' States of O + O2

We report new potential energy surfaces for six coupled 5A' states and 14 coupled 3A' states of O3. The new surfaces are created by parametrically managed diabatization by deep neural network (PM-DDNN). The PM-DDNN method uses calculated adiabatic potential energy surfaces to discover and fit an underlying adiabatic-equivalent set of diabatic surfaces and their couplings and obtains the fit to the adiabatic surfaces by diagonalization of the diabatic potential energy matrix (DPEM). The procedure yields the adiabatic surfaces and their gradients, as well as the DPEM and its gradient. If desired one can also compute the nonadiabatic coupling due to the transformation. The present work improves on previous work by using a new coordinate to guide the decay of the neural network contribution to the many-body fit to the whole DPEM. The main objective was to obtain smoother potentials than the previous ones with better suitability for dynamics calculations, and this was achieved. Furthermore, we obtained suitably small deviations from the input reference data. For the six coupled 5A' surfaces, the 60,366 data below 10 eV are fit with a mean unsigned error (MUE) of 49 meV, and for the 14 coupled 3A' surfaces, the 76,733 data below 10 eV are fit with an MUE of 28 meV. The data below 5 eV fit even more accurately with MUEs of 37 meV (5A') and 20 meV (3A').

Molecular Dynamics Study of the Photodegradation of Polymeric Chains

The development of reusable polymeric materials inspires an attempt to combine renewable biomass with upcycling to form a biorenewable closed system. It has been reported that 2,5-furandicarboxylic acid (FDCA) can be recovered for recycling when incorporated as monomers into photodegradable polymeric systems. Here, we conduct density functional theory (DFT) studies with periodic boundary conditions on microscopic structures involved in the photodegradation of polymeric chains incorporating FDCA and 2-nitro-1,3-benzenedimethanol. The photodegradation process of polymeric chains is studied using time-dependent excited-state molecular dynamics (TDESMD) in vacuum and aqueous environments. Changes in the photophysical properties for reaction intermediates are characterized by ground-state observables. The distribution of reaction intermediates and products is obtained from TDESMD trajectories using cheminformatics techniques. Results show that a higher degree of polymeric chain degradation is achieved in the vacuum environment. Additionally, one finds that the FDCA molecule is recoverable in the aqueous environment, in qualitative agreement with experimental findings.

Hot Carrier Dynamics at Ligated Silicon(111) Surfaces: A Computational Study

We provide a case-study for thermal grafting of benzenediazonium bromide onto a hydrogenated Si(111) surface using ab initio molecular dynamics (AIMD) calculations. A sequence of reaction steps is identified in the AIMD trajectory, including the loss of N₂ from the diazonium salt, proton transfer from the surface to the bromide ion that eliminates HBr, and deposition of the phenyl group onto the surface. We next assess the influence of the phenyl groups on photophysics of hydrogen-terminated Si(111) slabs. The nonadiabatic couplings necessary for a description of the excited-state dynamics are calculated by combining ab initio electronic structures and reduced density matrix formalism with Redfield theory. The phenyl-terminated slab shows reduced nonradiative relaxation and recombination rates of hot charge carriers in comparison with the hydrogen-terminated slab. Altogether, our results provide atomistic insights revealing that (i) the diazonium salt thermally decomposes at the surface allowing the formation of covalently bonded phenyl group, and (ii) the coverage of phenyl groups on the surface slows down charge carrier cooling driven by electron-phonon interactions, which increases photoluminescence efficiency at the near-infrared spectral region.

TURBOMOLE: Modular program suite for ab initio quantum-chemical and condensed-matter simulations

TURBOMOLE is a collaborative, multi-national software development project aiming to provide highly efficient and stable computational tools for quantum chemical simulations of molecules, clusters, periodic systems, and solutions. The TURBOMOLE software suite is optimized for widely available, inexpensive, and resource-efficient hardware such as multi-core workstations and small computer clusters. TURBOMOLE specializes in electronic structure methods with outstanding accuracy-cost ratio, such as density functional theory including local hybrids and the random phase approximation (RPA), GW-Bethe-Salpeter methods, second-order Møller-Plesset theory, and explicitly correlated coupled-cluster methods. TURBOMOLE is based on Gaussian basis sets and has been pivotal for the development of many fast and low-scaling algorithms in the past three decades, such as integral-direct methods, fast multipole methods, the resolution-of-the-identity approximation, imaginary frequency integration, Laplace transform, and pair natural orbital methods. This review focuses on recent additions to TURBOMOLE's functionality, including excited-state methods, RPA and Green's function methods, relativistic approaches, high-order molecular properties, solvation effects, and periodic systems. A variety of illustrative applications along with accuracy and timing data are discussed. Moreover, available interfaces to users as well as other software are summarized. TURBOMOLE's current licensing, distribution, and support model are discussed, and an overview of TURBOMOLE's development workflow is provided. Challenges such as communication and outreach, software infrastructure, and funding are highlighted.

Real-time tracking of the entangled pathways in the multichannel photodissociation of acetaldehyde

The roaming mechanism, an unconventional reaction path, was discovered more than a decade ago in the studies of formaldehyde photodissociation, H₂CO → H₂ + CO. Since then, observations of roaming have been claimed in numerous photochemical processes. A closer examination of the presented data, however, revealed that evidence for roaming is not always unequivocal, and some of the conclusions could be misleading. We report here an in-depth, joint experimental and theoretical study of the title reaction. By tracking the time-evolution of the pair-correlated product state distributions, we decipher the competing, interwoven reaction pathways that lead to the radical (CH₃ + HCO) and molecular (CH₄ + CO) products. Possible roaming pathways are then elucidated and a more precise descriptor of the phenomenon is delineated.

Non-adiabatic Excited-State Molecular Dynamics: Theory and Applications for Modeling Photophysics in Extended Molecular Materials

Optically active molecular materials, such as organic conjugated polymers and biological systems, are characterized by strong coupling between electronic and vibrational degrees of freedom. Typically, simulations must go beyond the Born-Oppenheimer approximation to account for non-adiabatic coupling between excited states. Indeed, non-adiabatic dynamics is commonly associated with exciton dynamics and photophysics involving charge and energy transfer, as well as exciton dissociation and charge recombination. Understanding the photoinduced dynamics in such materials is vital to providing an accurate description of exciton formation, evolution, and decay. This interdisciplinary field has matured significantly over the past decades. Formulation of new theoretical frameworks, development of more efficient and accurate computational algorithms, and evolution of high-performance computer hardware has extended these simulations to very large molecular systems with hundreds of atoms, including numerous studies of organic semiconductors and biomolecules. In this Review, we will describe recent theoretical advances including treatment of electronic decoherence in surface-hopping methods, the role of solvent effects, trivial unavoided crossings, analysis of data based on transition densities, and efficient computational implementations of these numerical methods. We also emphasize newly developed semiclassical approaches, based on the Gaussian approximation, which retain phase and width information to account for significant decoherence and interference effects while maintaining the high efficiency of surface-hopping approaches. The above developments have been employed to successfully describe photophysics in a variety of molecular materials.

Multistate hybrid time-dependent density functional theory with surface hopping accurately captures ultrafast thymine photodeactivation

We report an efficient analytical implementation of first-order nonadiabatic derivative couplings between arbitrary Born–Oppenheimer states in the hybrid time-dependent density functional theory (TDDFT) framework using atom-centered basis functions.

Tuning the photoreactivity of Z-hexatriene photoswitches by substituents - a non-adiabatic molecular dynamics study

To understand how substituents can be used to increase the quantum yield of photochemical electrocyclic ring-closing of the Z-hexa-1,3,5-triene (HT) photoswitch forming cyclohexadiene (CHD), we investigate the S1 photo dynamics of HT and its derivatives 2,5-dimethyl-HT (DMHT), 2-isopropyl-5-methyl-HT (2,5-IMHT), 1-isopropyl-4-methyl-HT (1,4-IMHT), and 2,5-diisopropyl-HT (DIHT) using time-dependent density functional theory surface hopping dynamics. We report detailed photoproduct distributions, formation mechanisms, branching ratios, and wavelength-dependent product quantum yields. Most products have been confirmed experimentally and include all-trans HT derivatives, cyclopropanes, cyclobutenes, cyclopentene, cyclohexadienes, and bicyclic compounds. Regarding CHD formation, we find that for the 2,5-substituted derivatives DMHT, 2,5-IMHT, and DIHT, the branching ratios increase with increasing size of the substituents. In contrast the branching ratios of the E/Z-isomerization decrease with increasing size of the substituents. Due to steric interactions, increasing the size of the substituents increases the amount of gZg rotamers in the ground state, which are prone to CHD formation and have lower E/Z-isomerization probability. Furthermore, we find [1,4], [1,5], and [1,6]-sigmatropic hydrogen shift reactions occurring at large percentages (5% to 15%); for sterical reasons these reactions stem from tZg conformers. DIHT shows the lowest percentage of side product formation among the 2,5-substituted molecules and highest CHD branching ratio; its CHD quantum yield can be increased up to more than 64%, by excitation of DIHT on the red tail of its absorption spectrum, whereas the Z/E-isomerization is reduced below 5% and side reactions practically vanish. This makes DIHT the best candidate for applications in molecular switches.

First-Principles Prediction of Wavelength-Dependent Product Quantum Yields

We present a method to predict wavelength-dependent product quantum yields (PQYs) for photochemical reactions and applied it to Z/E-isomerization and several ring-closing reactions of Z-2,5-dimethyl-1,3,5-hexatriene and truncated previtamin D. Using branching ratios from surface hopping molecular dynamics, individual trajectories are correlated with the absorption spectra of their initial structures. The wavelength-dependent PQYs are computed by dividing the average spectrum of the initial structures of the product-forming trajectories by the average spectrum of all initial structures. Accurate absorption spectra are calculated using the correlated ADC(2) method with an implicit solvent. Calculations reproduce the experimentally found trend of increasing six-ring formation and decreasing Z/E-isomerization on the red side of the spectrum. Over all seven reactions studied, the mean absolute error (MAE) between experimental and calculated PQYs (MAE) amounts to 8.1%, and the largest MAE is 18.6%. For four reactions, predicted values agree quantitatively with experiments within 5.6%.

Unraveling Photodimerization of Cyclohexasilane from Molecular Dynamics Studies

Photoinduced reactions of a pair of cyclohexasilane (CHS) monomers are explored by time-dependent excited-state molecular dynamics (TDESMD) calculations. In TDESMD trajectories, one observes vivid reaction events including dimerization and fragmentation. A general reaction pathway is identified as (i) ring-opening formation of a dimer, (ii) rearrangement induced by bond breaking, and (iii) decomposition through the elimination of small fragments. The identified pathway supports the chemistry proposed for the fabrication of silicon-based materials using CHS as a precursor. In addition, we find dimers have smaller HOMO-LUMO gaps and exhibit a red shift and line-width broadening in the computed photoluminescence spectra compared with a pair of CHS monomers.

Photofragmentation of Tetranitromethane: Spin-Unrestricted Time-Dependent Excited-State Molecular Dynamics

In this study, the photofragmentation dynamics of tetranitromethane (TNM) is explored by a spin-unrestricted time-dependent excited-state molecular dynamics (u-TDESMD) algorithm based on Rabi oscillations and principles similar to trajectory surface hopping, with a midintensity field approximation. The leading order process is represented by the molecule undergoing cyclic excitations and de-excitations. During excitation cycles, the nuclear kinetic energy is accumulated to overcome the dissociation barriers in the reactant and a sequence of intermediates. The dissociation pathway includes the ejection of NO₂ groups followed by the formation of NO and CO. The simulated mass spectra at the ab initio level, based on the bond length in possible fragments, are extracted from simulation trajectories. The recently developed methodology has the potential to model and monitor photoreactions with open-shell intermediates and radicals.

Mechanism of photocatalytic water oxidation on small TiO2 nanoparticles

We present the first unconstrained nonadiabatic molecular dynamics (NAMD) simulations of photocatalytic water oxidation by small hydrated TiO2 nanoparticles using Tully surface hopping and time-dependent density functional theory. The results indicate that ultrafast electron-proton transfer from physisorbed water to the photohole initiates the photo-oxidation on the S1 potential energy surface. The new mechanism readily explains the observation of mobile hydroxyl radicals in recent experiments. Two key driving forces for the photo-oxidation reaction are identified: localization of the electron-hole pair and stabilization of the photohole by hydrogen bonding interaction. Our findings illustrate the scope of recent advances in NAMD methods and emphasize the importance of explicit simulation of electronic excitations.

The role of tachysterol in vitamin D photosynthesis – a non-adiabatic molecular dynamics study

To investigate the role of tachysterol in the regulation of vitamin D photosynthesis, we studied its absorption properties and photodynamics by ab initio methods and non-adiabatic molecular dynamics.