Simulation of the photodeactivation of formamide in the nO-pi* and pi-pi* states: an ab initio on-the-fly surface-hopping dynamics study

Photodissociation of leucine-enkephalin protonated peptide: an experimental and theoretical perspective

Understanding the competing processes that govern far ultraviolet photodissociation (FUV-PD) of biopolymers such as proteins is a challenge. Here, we report a combined experimental and theoretical investigation of FUV-PD of protonated leucine-enkephalin pentapeptide ([YGGFL + H]+) in the gas-phase. Time-dependent density functional theory (TD-DFT) calculations in combination with experiments and previous results for amino acids and shorter peptides help in rationalizing the evolution of the excited states. The results confirm that fragmentation of [YGGFL + H]+ results mainly from vibrationally excited species in the ground electronic state, populated after internal conversion. We also propose fragmentation mechanisms for specific photo-fragments such as tyrosine side chain loss (with an extra hydrogen) or hydrogen loss. In general, we observe the same mechanisms as for smaller peptides or protonated Tyr and Phe, that are not quenched by the presence of other amino acids. Nevertheless, we also found some differences, as for H loss, in part due to the fact that the charge is solvated by the peptide chain and not only by the COOH terminal group.

A Local Diabatisation Method for Two-State Adiabatic Conical Intersections

A methodology to locally characterize conical intersections (CIs) between two adiabatic electronic states for which no nonadiabatic coupling (NAC) vectors are available is presented. Based on the Hessian and gradient at the CI, the branching space coordinates are identified. The potential energy surface around the CI in the branching space is expressed in the diabatic representation, from which the NAC vectors can be calculated in a wave-function-free, energy-based approach. To demonstrate the universality of the developed methodology, the minimum-energy CI (MECI) between the first (S1) and second (S2) singlet excited states of formamide is investigated at the state-averaged complete active space self-consistent field (SA-CASSCF) and extended multistate complete active space second-order perturbation theory (XMS-CASPT2) levels of theory. In addition, the asymmetrical MECI between the ground state (S0) and S1 of cyclopropanone is evaluated using SA-CASSCF, as well as (ME)CIs between the S1 and S2 states of benzene using SA-CASSCF and time-dependent density functional theory (TDDFT). Finally, a CI between the S1 and S2 excited states of thiophene was analyzed using TDDFT.

NWChem: Recent and Ongoing Developments

This paper summarizes developments in the NWChem computational chemistry suite since the last major release (NWChem 7.0.0). Specifically, we focus on functionality, along with input blocks, that is accessible in the current stable release (NWChem 7.2.0) and in the "master" development branch, interfaces to quantum computing simulators, interfaces to external libraries, the NWChem github repository, and containerization of NWChem executable images. Some ongoing developments that will be available in the near future are also discussed.

The primary photo-dissociation dynamics of aqueous formamide and dimethylformamide

We study the primary dissociation dynamics of aqueous formamide (HCONH2) and dimethylformamide (HCON(CH3)2) induced by photo-excitation at λ = 200 nm. The photolysis is recorded with sub-picosecond time resolution by...

An Ab Initio Multiple Cloning Method for Non-Adiabatic Excited-State Molecular Dynamics in NWChem

The recently developed ab initio multiple cloning (AIMC) approach based on the multiconfigurational Ehrenfest (MCE) method provides a powerful and accurate way of describing the excited-state dynamics of molecular systems. The AIMC method is a controlled approximation to nonadiabatic dynamics with a particular strength in the proper description of decoherence effects because of the branching of vibrational wavepackets at a level crossing. Here, we report a new implementation of the AIMC algorithm in the open source NWChem computational chemistry program. The framework combines linear-response time-dependent density functional theory with Ehrenfest mean-field theory to determine the equations of motion for classical trajectories. The multidimensional wave function is decomposed into a superposition of Gaussian coherent states guided by Ehrenfest trajectories (i.e., MCE approach), which can clone with fully quantum mechanical amplitudes and phases. By using an efficient time-derivative based nonadiabatic coupling approach within the AIMC method, all observables are calculated on-the-fly in the nonadiabatic molecular dynamics process. As a representative example, we apply our implementation to study the ultrafast photoinduced electronic and vibrational energy transfer in a pyridine molecule. The effects of the cloning procedure on electronic and vibrational coherence, relaxation and unidirectional energy transfer are discussed. This new AIMC implementation provides a high-level nonadiabatic molecular dynamics framework for simulating photoexcited dynamics in complex molecular systems and experimentally relevant ultrafast spectroscopic probes, such as nonlinear coherent optical and X-ray signals.

Unimolecular photolysis mechanisms of formamide: the effect of thermal energy

The C–N dissociation is the most favorable in S₁ state, and radical products (1) are precursors in other channels.

Unimolecular decomposition of formamide via direct chemical dynamics simulations

Classical chemical dynamics simulations show that formamide (NH₂CHO) can dissociate via multiple pathways, either by direct dissociations or via intramolecular rearrangements to different isomers followed by dissociation.

Nonradiative Relaxation Mechanisms of UV Excited Phenylalanine Residues: A Comparative Computational Study

The present work is directed toward understanding the mechanisms of excited state deactivation in three neutral model peptides containing the phenylalanine residue. The excited state dynamics of theγL(g+)folded form of N-acetylphenylalaninylamide (NAPA B) and its amide-N-methylated derivative (NAPMA B) is reviewed and compared to the dynamics of the monohydrated structure of NAPA (NAPAH). The goal is to unravel how the environment, and in particular solvation, impacts the photodynamics of peptides. The systems are investigated using reaction path calculations and surface hopping nonadiabatic dynamics based on the coupled cluster doubles (CC2) method and time-dependent density functional theory. The work emphasizes the role that excitation transfer from the phenylππ*to amidenπ*state plays in the deactivation of the three systems and shows how the ease of out-of-plane distortions of the amide group determines the rate of population transfer between the two electronic states. The subsequent dynamics on thenπ*state is barrierless along several pathways and leads to fast deactivation to the ground electronic state.

Photo-stability of peptide-bond aggregates: N-methylformamide dimers

The formation of weakly-bound dimers of N-methylformamide (NMF) and the photochemistry of these dimers after irradiation at 248 nm were explored using matrix-isolation spectroscopy. Calculations were used to characterize the diverse isomers and assign their IR spectra; non-adiabatic dynamics was simulated to understand their photo-deactivation mechanism. The most stable dimers, and , were obtained by trans-trans aggregation (N-HO[double bond, length as m-dash]C interactions) and could be identified in the matrix. The main products formed after irradiation are the trans-cis dimers ( and ), also stabilized by N-HO[double bond, length as m-dash]C interactions. In contrast to the photochemistry of the monomers, no dissociative products were observed after 248 nm irradiation of the dimers. The absence of dissociative products can be explained by a proton-transfer mechanism in the excited state that is faster than the photo-dissociative mechanism. The fact that hydrogen bonding has such a significant effect on the photochemical stability of NMF has important implications to understand the stability of peptide-bonded systems to UV irradiation.

Meteorite-catalyzed syntheses of nucleosides and of other prebiotic compounds from formamide under proton irradiation

Liquid formamide has been irradiated by high-energy proton beams in the presence of powdered meteorites, and the products of the catalyzed resulting syntheses were analyzed by mass spectrometry. Relative to the controls (no radiation, or no formamide, or no catalyst), an extremely rich, variegate, and prebiotically relevant panel of compounds was observed. The meteorites tested were representative of the four major classes: iron, stony iron, chondrites, and achondrites. The products obtained were amino acids, carboxylic acids, nucleobases, sugars, and, most notably, four nucleosides: cytidine, uridine, adenosine, and thymidine. In accordance with theoretical studies, the detection of HCN oligomers suggests the occurrence of mechanisms based on the generation of radical cyanide species (CN·) for the synthesis of nucleobases. Given that many of the compounds obtained are key components of extant organisms, these observations contribute to outline plausible exogenous high-energy-based prebiotic scenarios and their possible boundary conditions, as discussed.

Photodynamic behavior of electronic coupling in a N-methylformamide dimer

The role of the bridging water molecules has been studied during the excited state photodynamics of a N-methylformamide dimer in complex with water molecules employing the complete active space self-consistent field (CASSCF) and CAS perturbation theory (CASPT2) methods.

Guanidine and guanidinium cation in the excited state--theoretical investigation

Diverse ab initio and density-functional-theory methods were used to investigate geometries, energies, and electronic absorption spectra of guanidine and its protonated form, as well as their photo-deactivation processes. It was shown that the guanidine is a weakly absorbing species with the excitation spectrum consisting mostly of transitions to the Rydberg excited states and one valence n-π4 state. The lowest energy band has a maximum at ca. 6.9 eV (∼180 nm). The protonation of guanidine affects its excitation spectrum substantially. A major shift of the Rydberg states to higher energies is clearly visible and strongly absorbing transitions from the ground state to the π3-π4 and π2-π4 states appears at 7.8 eV (∼160 nm). Three low-lying conical intersections (two for guanidine and one for protonated guanidine) between the ground state and the first excited singlet state were located. They are accessible from the Franck-Condon region through amino N-H stretching and out-of-plane deformations in guanidine and protonated guanidine, respectively. The relaxation of the π3-3s Rydberg state via amino N-H bond stretching was hindered by a barrier. The nondissociated conical intersection in protonated guanidine mediates the radiationless deactivation of the compound after excitation into the π3-π4 state. This fact is detrimental for the photostability of guanidine, since its conjugate acid is stable in aqueous solution over a wide pH range and in protein environment, where guanidinium moiety in arginine is expected to be in a protonated form.

Nonadiabatic excited-state molecular dynamics: modeling photophysics in organic conjugated materials

To design functional photoactive materials for a variety of technological applications, researchers need to understand their electronic properties in detail and have ways to control their photoinduced pathways. When excited by photons of light, organic conjugated materials (OCMs) show dynamics that are often characterized by large nonadiabatic (NA) couplings between multiple excited states through a breakdown of the Born-Oppenheimer (BO) approximation. Following photoexcitation, various nonradiative intraband relaxation pathways can lead to a number of complex processes. Therefore, computational simulation of nonadiabatic molecular dynamics is an indispensable tool for understanding complex photoinduced processes such as internal conversion, energy transfer, charge separation, and spatial localization of excitons. Over the years, we have developed a nonadiabatic excited-state molecular dynamics (NA-ESMD) framework that efficiently and accurately describes photoinduced phenomena in extended conjugated molecular systems. We use the fewest-switches surface hopping (FSSH) algorithm to treat quantum transitions among multiple adiabatic excited state potential energy surfaces (PESs). Extended molecular systems often contain hundreds of atoms and involve large densities of excited states that participate in the photoinduced dynamics. We can achieve an accurate description of the multiple excited states using the configuration interaction single (CIS) formalism with a semiempirical model Hamiltonian. Analytical techniques allow the trajectory to be propagated "on the fly" using the complete set of NA coupling terms and remove computational bottlenecks in the evaluation of excited-state gradients and NA couplings. Furthermore, the use of state-specific gradients for propagation of nuclei on the native excited-state PES eliminates the need for simplifications such as the classical path approximation (CPA), which only uses ground-state gradients. Thus, the NA-ESMD methodology offers a computationally tractable route for simulating hundreds of atoms on ~10 ps time scales where multiple coupled excited states are involved. In this Account, we review recent developments in the NA-ESMD modeling of photoinduced dynamics in extended conjugated molecules involving multiple coupled electronic states. We have successfully applied the outlined NA-ESMD framework to study ultrafast conformational planarization in polyfluorenes where the rate of torsional relaxation can be controlled based on the initial excitation. With the addition of the state reassignment algorithm to identify instances of unavoided crossings between noninteracting PESs, NA-ESMD can now be used to study systems in which these so-called trivial unavoided crossings are expected to predominate. We employ this technique to analyze the energy transfer between poly(phenylene vinylene) (PPV) segments where conformational fluctuations give rise to numerous instances of unavoided crossings leading to multiple pathways and complex energy transfer dynamics that cannot be described using a simple Förster model. In addition, we have investigated the mechanism of ultrafast unidirectional energy transfer in dendrimers composed of poly(phenylene ethynylene) (PPE) chromophores and have demonstrated that differential nuclear motion favors downhill energy transfer in dendrimers. The use of native excited-state gradients allows us to observe this feature.

Photodissociation dynamics of propargylene, HCCCH

We report a joint theoretical and experimental study on the photodissociation of the C₃H₂ isomer propargylene, HCCCH, combining velocity map imaging with nonadiabatic surface hopping calculations.

Decomposition pathways of the neutral and protonated formamide in some lower-lying excited states

Unimolecular decompositions of neutral (NH2CHO) and protonated (NH3CHO(+)) formamide, an active precursor of biomolecules in prebiotic chemistry, are investigated in the ground (S0) and first triplet (T1) and singlet (S1) excited states. Different decomposition channels including the homolytic bond dissociations, dehydration, decarbonylation, dehydrogenation, etc., are explored using coupled-cluster theory (CCSD(T)/CBS method) for both S0 and T1 states and RASPT2(18,15)/6-31G(d,p) computations for the S1 state. On S1 and T1 energy surfaces, formamide preferentially follows C-N homolytic bond cleavages forming NH2 + HCO radical pairs. Formation of HCN and HNC from dehydration of neutral and protonated formamide via formimic acid and aminohydroxymethylene isomers has higher energy barriers. A strong stabilization upon triplet excitation of the two latter isomers significantly facilitates the interconversions between isomers, and thus considerably reduces the energy barriers for dehydration pathways. The most probable pathways for HCN and HNC generation are found to be dehydration of formamide in the T1 state. Dehydration pathways from the neutral S1 and protonated T1 forms lead to stable complexes of HCN and HNC with water but are associated with large energy barriers. Overall, in the lower-lying excited states of either neutral or protonated formamide, dehydration is not competitive with homolytic C-N bond cleavages, which finally lead to formation of CO.

Unravelling the details of vitamin D photosynthesis by non-adiabatic molecular dynamics simulations

We investigate the photodynamics of vitamin D derivatives by a fully analytical implementation of the linear response time-dependent density functional theory surface hopping method (LR-TDDFT-SH). Our study elucidates the dynamics of the processes involved in vitamin D formation at the molecular level and with femtosecond resolution. We explain the major experimental findings and provide new insights that cannot directly be obtained from experiments: firstly, we investigate the dynamics of the photoinduced ring-opening of provitamin D (Pro) and cyclohexadiene (CHD) and the subsequent rotational isomerization. In agreement with recent experiments and CC2 calculations, only the bright S(1) state is involved in the ring-opening reaction. Our calculations confirm the experimentally reported 5 : 1 ratio between the excited state lifetimes of Pro and CHD. The longer lifetimes of Pro are attributed to steric constraints of the steroid skeleton and to temperature effects, both emerging directly from our simulations. For CHD and Pro, we present an explanation of the biexponential decay recently reported by Sension and coworkers [Tang et al., J. Phys. Chem., 2011, 134, 104503]: our calculations suggest that the fast and slow components arise from a reactive and an unreactive reaction pathway, respectively. Secondly, we assess the wavelength dependent photochemistry of previtamin D (Pre). Using replica exchange molecular dynamics we sample the Pre conformers present at thermal equilibrium. Based on this ensemble we explain the conformation dependent absorption and the essential features of Pre photochemistry. Consistent with the experiments, we find ring-closure to occur mostly after excitation of the cZc conformers and at lower energies, whereas Z/E isomerization of the central double bond preferably occurs after excitation at higher energies. For the isomerization we provide the first theoretical evidence of the proposed hula-twist mechanism. Our results show that LR-TDDFT-SH is a highly valuable tool for studying the photochemistry of moderately large systems, even though challenges remain in the vicinity of conical intersections.

Photophysics and deactivation pathways of thymine

Combined complete active space perturbation theory (CASPT2) and multireference configuration interaction calculations with single and double excitations (MR-CISD) were performed in order to explore possible deactivation pathways of thymine after photoexcitation. Equilibrium geometries are reported together with a total of eight extremes (minima or maxima) on the crossing seam (MXS), through which such radiationless transitions may occur. Furthermore, conformational analysis allows grouping these conical intersections in five distinct types. Reaction paths were calculated connecting the S1 (1)n pi* minimum with the lowest-energy MXS of each group. Two distinct types of paths were observed, both with features that should delay the internal conversion to the ground state. This is shown to provide a possible explanation for the relatively long excited-state lifetime of thymine.