The photodissociation of solvated cyclopropanone and its hydrate explored via non-adiabatic molecular dynamics using ΔSCF.

The decay of cyclopropanone is a typical example of a photodecomposition process. Ethylene and carbon monoxide are formed following the excitation to the first singlet excited state through a symmetrical or asymmetrical pathway. The results obtained with non-adiabatic molecular dynamics (NAMD) using the delta self-consistent field (ΔSCF) method correspond well to previous experimental and multireference theoretical studies carried out in the gas phase. Moreover, this efficient methodology allows NAMD simulations of cyclopropanone in aqueous solution to be performed, which reveal analogue deactivation mechanisms, but a shorter lifetime and reduced photodissociation as compared to the gas-phase. The excited state dynamics of cyclopropanone hydrate, an enzyme inhibitor, in an aqueous environment are reported as well. Cyclopropanone hydrate strongly interacts with the surrounding solvent via the formation of hydrogen bonds. Excitation to the first singlet excited state shows an asymmetric pathway with cyclopropanone hydrate and propionic acid as the main photoproducts.

The ΔSCF method for non-adiabatic dynamics of systems in the liquid phase.

Computational studies of ultrafast photoinduced processes give valuable insights into the photochemical mechanisms of a broad range of compounds. In order to accurately reproduce, interpret, and predict experimental results, which are typically obtained in a condensed phase, it is indispensable to include the condensed phase environment in the computational model. However, most studies are still performed in vacuum due to the high computational cost of state-of-the-art non-adiabatic molecular dynamics (NAMD) simulations. The quantum mechanical/molecular mechanical (QM/MM) solvation method has been a popular model to perform photodynamics in the liquid phase. Nevertheless, the currently used QM/MM embedding techniques cannot sufficiently capture all solute-solvent interactions. In this Perspective, we will discuss the efficient ΔSCF electronic structure method and its applications with respect to the NAMD of solvated compounds, with a particular focus on explicit quantum mechanical solvation. As more research is required for this method to reach its full potential, some challenges and possible directions for future research are presented as well.

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.

Timescales of N-H bond dissociation in pyrrole: a nonadiabatic dynamics study.

The excitation wavelength dependent photodynamics of pyrrole are investigated by nonadiabatic trajectory-surface-hopping dynamics simulations based on time dependent density functional theory (TDDFT) and the algebraic diagrammatic construction method to the second order (ADC(2)). The ADC(2) results confirm that the N-H bond dissociation occurring upon excitation at the origin of the first excited state, S1(πσ*), is driven by tunnelling [Roberts et al., Faraday Discuss., 2013, 163, 95] as a barrier of ΔE = 1780 cm(-1) traps the population in a quasi-bound minimum. Upon excitation to S1(πσ*) in the wavelength range of 236-240 nm, direct dissociation of the N-H bond takes place with a time constant of 28 fs. The computed time constant is in very good agreement with recently reported measurements. Excitation to the lowest B2 state in the 198-202 nm range returns a time constant for N-H fission of 48 fs at the B3LYP/def2-TZVPD level, in perfect agreement with the experiment [Roberts et al. Faraday Discuss., 2013, 163, 95]. For the same wavelength range the ADC(2)/aug-cc-pVDZ decay constant is more than three times longer than the experimentally reported one. The accuracy of the B3LYP/def2-TZVPD dynamics is checked against reference complete-active-space second-order perturbation theory (CASPT2) calculations and explained in terms of correct topography of the ππ* surface and the lack of mixing between the ππ* and the 3px Rydberg states which occurs in the ADC(2) method.

Non-radiative relaxation of UV photoexcited phenylalanine residues: probing the role of conical intersections by chemical substitution.

A conformation-selective photophysics study in phenylalanine model peptides, combining pump-probe gas phase experiments and excited state calculations, highlights for the first time the quenching properties of a primary amide group (through its nπ* excited state) along with the effect of vibrational energy that facilitates access to the conical intersection area.

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.

The Role of Aqueous Solvation on the Intersystem Crossing of Nitrophenols.

The photochemistry of nitrophenols is a source of smog as nitrous acid is formed from their photolysis. Nevertheless, computational studies of the photochemistry of these widespread toxic molecules are scarce. In this work, the initial photodeactivation of ortho-nitrophenol and para-nitrophenol is modeled, both in gas phase and in aqueous solution to simulate atmospheric and aerosol environments. A large number of excited states, six for ortho-nitrophenol and 11 for para-nitrophenol, have been included and were all populated during the decay. Moreover, periodic time-dependent density functional theory (TDDFT) is used for both the explicitly included solvent and the solute. A comparison to periodic QM/MM (TDDFT/MM), with electrostatic embedding, is made, showing notable differences between the decays of solvated nitrophenols simulated with QM/MM and full (TD)DFT. A reduced intersystem crossing in aqueous solution could be observed thanks to the surface hopping approach using explicit, periodic TDDFT solvation including spin-orbit couplings.

Unraveling the mechanisms of nonradiative deactivation in model peptides following photoexcitation of a phenylalanine residue.

The mechanisms of nonradiative deactivation of a phenylalanine residue after near-UV photoexcitation have been investigated in an isolated peptide chain model (N-acetylphenylalaninylamide, NAPA) both experimentally and theoretically. Lifetime measurements at the origin of the first ππ* state of jet-cooled NAPA molecules have shown that (i) among the three most stable conformers of the molecule, the folded conformer NAPA B is ∼50-times shorter lived than the extended major conformer NAPA A and (ii) this lifetime is virtually insensitive to deuteration at the NH(2) and NH sites. Concurrent time-dependent density functional theory (TDDFT) based nonadiabatic dynamics simulations in the full dimensionality, carried out for the NAPA B conformer, provided direct insights on novel classes of ultrafast deactivation mechanisms, proceeding through several conical intersections and leading in fine to the ground state. These mechanisms are found to be triggered either (i) by a stretch of the N(Phe)H bond, which leads to an H-transfer to the ring, or (ii) by specific backbone amide distortions. The potential energy surfaces of the NAPA conformers along these critical pathways have been characterized more accurately using the coupled cluster doubles (CC2) method and shown to exhibit barriers that can be overcome with moderate excess energies. These results analyzed in the light of the experimental findings enabled us to assign the short lifetime of NAPA B conformer to a number of easily accessible exit channels from the initial ππ* surface, most importantly the one involving a transfer of electronic excitation to an nπ* surface, induced by distortions of the backbone peptide bond.

Very efficient generation of quinone methides through excited state intramolecular proton transfer to a carbon atom.

Irradiation of 2-phenyl-1-naphthol (6) in CH(3) CN/D(2) O (3:1) leads to very efficient incorporation of deuterium at the ortho-positions of the adjacent phenyl ring (overall Φ=0.73±0.07), along with minor incorporation at the naphthalene positions 5 and 8. These finding are explained by excited state intramolecular proton transfer (ESIPT) from the phenolic OH group to the corresponding carbon atoms, the main pathway giving rise to quinone methide (QM) 7, which has been characterized by LFP (τ≈20 ns; 460 nm). The ESIPT reaction paths have been explored with the second order approximate coupled cluster (CC2) method. In nonprotic solvents the ESIPT from the naphthol O-H to the ortho-position of the phenyl ring proceeds in a barrierless manner along the (1) L(a) energy surface via a conical intersection with the S(0) state, delivering 7. In aqueous solvent, clusters with H(2) O are formed wherein proton transfer (PT) to solvent and a H(2) O-mediated relay mechanism gives rise to naphtholates and QMs. The results are compared with 2-phenylphenol (3) that also undergoes barrierless ESIPT giving a QM via a conical intersection. However, due to an unfavorable conformation in the ground state, the quantum efficiency for ESIPT of 3 is significantly lower (Φ for D-exchange=0.041). These results show that ESIPT from phenol to carbon need not be an intrinsically inefficient process.

Desmotropy, polymorphism, and solid-state proton transfer: four solid forms of an aromatic o-hydroxy Schiff base.

The Schiff base derived from salicylaldehyde and 2-amino-3-hydroxypyridine affords a diversity of solid forms, two polymorphic pairs of the enol-imino (D1 a and D1 b) and keto-amino (D2 a and D2 b) desmotropes. The isolated phases, identified by IR spectroscopy, X-ray crystallography, and (13)C cross-polarization/magnetic angle spinning (CP/MAS) NMR spectroscopy, display essentially planar molecular conformations characterized by strong intramolecular hydrogen bonds of the O-H⋅⋅⋅N (D1) or N-H⋅⋅⋅O (D2) type. A change in the position of the proton within this O⋅⋅⋅H⋅⋅⋅N system is accompanied by substantially different molecular conformations and, subsequently, by divergent supramolecular architectures. The appearance and interconversion conditions for each of the four phases have been established on the basis of a number of solution and solvent-free experiments, and evaluated against the results of computational studies. Solid phases readily convert into the most stable form (D1 a) upon exposure to methanol vapor, heating, or by mechanical treatment, and these transformations are accompanied by a change in the color of the sample. The course of thermally induced transformations has been monitored in detail by means of temperature-resolved powder X-ray diffraction and infrared spectroscopy. Upon dissolution, all forms equilibrate immediately, as confirmed by NMR and UV/Vis spectroscopy in several solvents, with the equilibrium shifted far towards the enol tautomer. This study reveals the significance of peripheral groups in the stabilization of metastable tautomers in the solid state.

Photoinduced dynamics of formic acid monomers and dimers: the role of the double hydrogen bond.

Nonadiabatic dynamics in the framework of time-dependent density functional theory was used to simulate gas-phase relaxation dynamics of pairs of conformations of formic acid monomers (cis and trans FAM) and dimers (acyclic aFAD and cyclic cFAD). In the early phase of the excited state dynamics, elongation of the C═O bond and pyramidalization of the carbon atom is observed in both FAM and FAD. Subsequently, the photodynamics of FAM is shown to be dominated by fragmentation processes occurring mostly in the excited state and resulting in HCO and OH radicals. In only a few cases does the dissociation take place from the vibrationally excited ground electronic state, whereby CO and H(2)O are the major reaction products. In the dimers, single proton transfer triggers ultrafast relaxation to the ground electronic state. In the single hydrogen bonded dimer about half of the trajectories dissociate into electronically excited monomers, whereas this potentially destructive dissociation is effectively suppressed in the double hydrogen bonded dimer. Upon relaxation to the ground electronic state, separation of FAD into monomers takes place, but without their further fragmentation on the time scale of the simulation. We conclude that the crucial difference between the FAM and FAD photodynamics is that the latter is dominated by nondestructive radiationless deactivation pathways during which a key protective role is assumed by the single (aFAD) or double (cFAD) intermonomer hydrogen bonds.

Spin-Orbit Couplings for Nonadiabatic Molecular Dynamics at the ΔSCF Level.

A procedure for the calculation of spin-orbit coupling (SOC) at the delta self-consistent field (ΔSCF) level of theory is presented. Singlet and triplet excited electronic states obtained with the ΔSCF method are expanded into a linear combination of singly excited Slater determinants composed of ground electronic state Kohn-Sham orbitals. This alleviates the nonorthogonality between excited and ground electronic states and introduces a framework, similar to the auxiliary wave function at the time-dependent density functional theory (TD-DFT) level, for the calculation of observables. The ΔSCF observables of the formaldehyde system were compared to reference TD-DFT values. Our procedure gives all components (energies, gradients, nonadiabatic couplings, and SOC terms) at the ΔSCF level of theory for conducting efficient, full-atomistic nonadiabatic molecular dynamics with intersystem crossing, particularly in condensed phase systems.

Fluxionality of hydrogen ligands in Fe(H)2(H2)(PEtPh2)3.

Extensive computational investigations along with additional quasielastic neutron scattering data were used to obtain a consistent picture of the extensive fluxionality of hydride and dihydrogen ligands in Fe(H)(2)(H(2))(PEtPh(2))(3) over a wide range of temperatures from 1.5 to 320 K. We were able to identify three different regimes in the dynamical processes based on activation energies obtained from line spectral broadening. The rotational tunneling lines (coherent exchange of the two hydrogens of the H(2) ligand) are broadened with increasing temperature by incoherent exchange up to about 80 K at which point they merge into a quasielastic spectrum from 100 K to about 225 K. The effective activation energies for the two regions are 0.14 and 0.1 kcal mol(-1), respectively. A third dynamical process with a higher activation energy of 0.44 kcal mol(-1) dominates above 225 K, which we attribute to a quantum dynamical exchange of dihydrogen and hydride ligands. Our detailed density functional theory (DFT) structural calculations involving the three functionals (B3LYP, TPSS, and wB97XD) provide a good account of the experimental structure and rotational barriers when only the hydrogen ligands are relaxed. Full relaxation of the "gas-phase" molecule, however, appears to occur to a greater degree than what is possible in the crystal structure. The classical dihydrogen-hydride exchange path involves a cis-dihydrogen and tetrahydride structure with energies of 6.49 and 7.38 kcal mol(-1), respectively. Experimental observation of this process with much lower energies would seem to suggest involvement of translational tunneling in addition to the rotational tunneling. Dynamics of this type may be presumed to be important in hydrogen spillover from metal particles, and therefore need to be elucidated in an effort to utilize this phenomenon.

ΔSCF with Subsystem Density Embedding for Efficient Nonadiabatic Molecular Dynamics in Condensed-Phase Systems.

An approach combining subsystem density embedding with the variational delta self-consistent field is presented, which extends current capabilities for excited-electronic-state calculations. It was applied on full-atomic nonadiabatic dynamics of a solvated diimide system, demonstrating that comparable accuracy can be achieved for this system for the investigated configuration space and with a shorter simulation time than the computationally more expensive conventional Kohn-Sham density functional theory-based method. This opens a new pragmatic technique for efficient simulation of nonadiabatic processes in the condensed phase, in particular, for liquids.

Trajectory Surface Hopping Nonadiabatic Molecular Dynamics with Kohn-Sham ΔSCF for Condensed-Phase Systems.

We present an efficient approach for surface hopping-based nonadiabatic dynamics in the condensed phase. For the systems studied, a restricted Kohn-Sham orbital formulation of the delta self-consistent field (ΔSCF) method was used for efficient calculation of excited electronic states. Time-dependent density functional theory (DFT) is applied to aid excited-state SCF convergence and provide guess electronic state densities. Aside from that the Landau-Zener procedure simplifies the surface hopping between electronic states. By utilizing the combined Gaussian and plane waves approach with periodic boundary conditions the method is easily applicable to full atomistic DFT simulations of condensed-phase systems and was used to study the nonradiative deactivation mechanism of photoexcited diimide in water solution.

Mechanisms of photostability in kynurenines: a joint electronic-structure and dynamics study.

Kynurenines are UV filters found in the human ocular lens which protect the retina from radiation damage. We report on ab initio investigations of the photochemistry of the cis and trans conformers of kynurenine and of an intramolecularly hydrogen-bonded conformer of 3-hydroxykynurenine O-ß-D-glucoside. We have explored the excited-state reaction paths for several radiationless excited-state deactivation processes in kynurenines. We show that electron-driven proton-transfer processes mediated by an excited state of charge-transfer character exhibit negligible barriers and that the relevant potential-energy profiles are lower in energy than the lowest absorbing ππ* state. In these proton-transfer processes, a proton moves from one of the amino groups of kynurenine to the keto group. We also report on nonadiabatic trajectory-surface-hopping molecular-dynamics simulations for photoexcited kynurenine. These simulations show that the cis and trans conformers of kynurenine deactivate on a femtosecond-to-picosecond time scale preferably via electron-driven proton transfer from one of the amino groups to the keto group. Cis kynurenine deactivates via a ring-N-H···O═C proton-transfer process. Trans kynurenine tends to undergo trans → cis isomerization before deactivating via the same process. These results suggest that the deactivation process involving the ring-amino group in the cis conformer of kynurenine is the most efficient excited-state deactivation process in kynurenines. The joint electronic-structure calculations and dynamics simulations provide a new level of mechanistic insight into the efficient UV-filtering capacity of kynurenines.

Calculating CD Spectra of Flexible Peptides: An Assessment of TD-DFT Functionals.

Calculations of CD spectra can provide important structural information for peptide systems. Although TD-DFT is an attractive method for these calculations, recent studies have pointed to problems with modeling charge transfer excitations. Motivated by these problems, we benchmarked the performances of CAM-B3LYP, ωB97X-D, M06-2X, B3LYP, and PBE0 against high level ab initio RICC2 calculations for selected peptide structures. Furthermore, we compared the performance of the functionals with the experimentally available data. Our results show that long-range corrected functionals (CAM-B3LYP and ωB97X-D) correlate relatively well with RICC2 calculations, as does the meta-hybrid M06-2X, while the global hybrid functionals (B3LYP and PBE0) exhibit the aforementioned charge-transfer artifacts. On the other hand, PBE0 and even more so M06-2X and B3LYP produce spectra in better agreement with the experimental data. We have clarified this apparent discrepancy by finding that the surplus charge-transfer excitations, exhibited by B3LYP and PBE0, seem to have a negligible contribution to the final spectra, once appropriate structural averaging is performed.