Search

The kinetics of many chemical reactions can be readily explained with a statistical approach, for example, using a form of transition state theory and comparing calculated Gibbs energies along the reaction coordinate(s). However, there are cases where this approach fails, notably when the vibrational relaxation of the molecule to its statistical equilibrium occurs on the same time scale as the reaction dynamics, whether it is caused by slow relaxation, a fast reaction, or both. These nonstatistical phenomena are then often explored computationally using (quasi)classical ab initio molecular dynamics by calculating a large number of trajectories while being prone to issues such as zero-point energy leakage. On the other side of the field, we see resource-intensive quantum dynamics simulations, which significantly limit the size of explorable systems. We find that using a Fermi's golden rule type of model for vibrational relaxation, based on anharmonic coupling constants, we can extract the same qualitative information while giving insights into how to enhance (or destroy) the bottlenecks causing the phenomena. We present this model as a middle ground for exploring complex nonstatistical behavior, capable of treating medium-sized organic molecules or biologically relevant fragments. We also cover the challenges involved, in particular quantifying the excess energy in terms of vibrational modes. Relying on readily available electronic structure methods and providing results in a simple master equation form, this model shows promise as a screening tool for opportunities in mode-selective chemistry without external control.

A combined computational and experimental study reveals that ortho -, meta - and para -aminobiphenyl isomers undergo distinctly different photochemical reactions involving proton transfer. Deuterium exchange experiments show that the ortho -isomer undergoes a facile photoprotonation at a carbon atom via excited-state intramolecular proton transfer (ESIPT). The meta -isomer undergoes water-assisted excited-state proton transfer (ESPT) and a photoredox reaction via proton-coupled electron transfer (PCET). The para -isomer undergoes a water-assisted ESPT reaction. All three reactions take place in the singlet excited-state, except for the photoredox process of the meta -isomer, which involves a triplet excited-state. Computations illustrate the important role of excited-state antiaromaticity relief in these photoreactions., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

The symmetric bissilyl-dione 3 reveals two well-separated n → π* absorption bands at λ max = 637 nm ( ε = 140 mol -1 dm 3 cm -1 ) and 317 nm ( ε = 2460 mol -1 dm 3 cm -1 ). Whereas excitation of 3 at λ = 360/365 nm affords an isolable siloxyketene 4 in excellent yields, irradiation at λ = 590/630 nm leads to the stereo-selective and quantitative formation of the siloxyrane 5. These remarkable wavelength-dependent rearrangements are based on the electronic and steric properties provided by the hypersilyl groups. While the siloxyketene 4 is formed via a hitherto unknown 1,3-hypersilyl migration via the population of a second excited singlet state (S 2 , λ max = 317 nm, a rare case of anti-Kasha reactivity), the siloxyrane 5 emerges from the first excited triplet state (T 1 via S 1 λ max = 637 nm). These distinct reaction pathways can be traced back to specific energy differences between the S 2 , S 1 and T 1 , an electronic consequence of the bissilyl substited α-dione (the "pearl"). The hypersilyl groups act as protective ''oyster shell", which are responsible for the clean formation of 4 and 5 basically omitting side products. We describe novel synthetic pathways to achieve hypersilyl substitution (3) and report an in-depth investigation of the photorearrangements of 3 using UV/vis, in situ IR, NMR spectroscopy and theoretical calculations., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

Challenging the basis of our chemical intuition, recent experimental evidence reveals the presence of a new type of intrinsic fluorescence in biomolecules that exists even in the absence of aromatic or electronically conjugated chemical compounds. The origin of this phenomenon has remained elusive so far. In the present study, we identify a mechanism underlying this new type of fluorescence in different biological aggregates. By employing non-adiabatic ab initio molecular dynamics simulations combined with a data-driven approach, we characterize the typical ultrafast non-radiative relaxation pathways active in non-fluorescent peptides. We show that the key vibrational mode for the non-radiative decay towards the ground state is the carbonyl elongation. Non-aromatic fluorescence appears to emerge from blocking this mode with strong local interactions such as hydrogen bonds. While we cannot rule out the existence of alternative non-aromatic fluorescence mechanisms in other systems, we demonstrate that this carbonyl-lock mechanism for trapping the excited state leads to the fluorescence yield increase observed experimentally, and set the stage for design principles to realize novel non-invasive biocompatible probes with applications in bioimaging, sensing, and biophotonics., (© 2023. The Author(s).)

This work shows how the N 1s photoemission (PE) spectrum of self-associated melamine molecules in aqueous solution has been successfully rationalized using an integrated computational approach encompassing classical metadynamics simulations and quantum calculations based on density functional theory (DFT). The first approach allowed us to describe interacting melamine molecules in explicit waters and to identify dimeric configurations based on π-π and/or H-bonding interactions. Then, N 1s binding energies (BEs) and PE spectra were computed at the DFT level for all structures both in the gas phase and in an implicit solvent. While pure π-stacked dimers show gas-phase PE spectra almost identical to that of the monomer, those of the H-bonded dimers are sensibly affected by NH···NH or NH···NC interactions. Interestingly, the solvation suppresses all of the non-equivalences due to the H-bonds yielding similar PE spectra for all dimers, matching very well our measurements.

Excited state intramolecular proton transfer (ESIPT) has been documented from an amino NH 2 group to a carbon atom of an adjacent aromatic ring. This finding changes the paradigm, as hitherto such processes have not been considered as plausible due to slow protonation of carbon and low (photo)acidity of the NH 2 group. The ESIPT was studied by irradiation of 2-(2-aminophenyl)naphthalene in CH 3 CN-D 2 O, whereupon regiospecific incorporation of deuterium takes place at the naphthalene position 1, with a quantum yield of Φ = 0.11. A synergy of experimental and computational investigations completely unraveled the mechanism of this important photochemical reaction. Upon excitation to the photoreactive S 2 (L a ) state, a favorable redistribution of charge sets the stage for ESIPT to the carbon atom in naphthalene position 1. H 2 O molecules are needed, as they increase the excitation energy and oscillator strength for the population of the S 2 (L a ) state. The gain in energy is used to surmount a small energy barrier on the pathway from the Franck-Condon geometry to the conical intersection with the S 0 , delivering aza-quinone methide.

UV irradiation of RNA leads to the formation of intra- and inter-strand crosslinks of cyclobutane type. Despite the importance of this reaction, relatively little is known about how the mutual orientation of the two bases affects the outcome of the reaction. Here we report a comparative nonadiabatic molecular dynamics study of face-to-back (F2B) and face-to-face (F2F) stacked uracil-water clusters. The computations were performed using the second-order algebraic-diagrammatic-construction (ADC(2)) method. We found that F2B stacked uracil-water clusters either relax non-reactively to the ground state by an ethylenic twist around the CC bond or remain in the lowest nπ* state in which the two bases gradually move away from each other. This finding is consistent with the low propensity for the formation of intra-strand cyclobutane dimers between adjacent RNA bases. On the contrary, in F2F stacked uracil-water clusters, in addition to non-reactive deactivation, we found a pro-reactive deactivation pathway, which may lead to the formation of cyclobutane uracil dimers in the electronic ground state. On a qualitative level, the observed photodynamics of F2F stacked uracil-water clusters explains the greater propensity of RNA to form inter-strand cyclobutane-type crosslinks.

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'.

Photoelectron angular distribution (PAD) in the laboratory frame for randomly oriented molecules is typically described by a single anisotropy parameter, the so-called asymmetry parameter. However, especially from a theoretical perspective, it is more natural to consider molecular photoionization by using a molecular frame. The molecular frame PADs (MFPADs) may be used to extract information about the electronic structure of the system studied. In the last decade, significant experimental efforts have been directed to MFPAD measurements. MFPADs are highly characterizing signatures of the final ionic states. In particular, they are very sensitive to the nature of the final state, which is embodied in the corresponding Dyson orbital. In our previous work on acetylacetone, a prototype system for studying intra-molecular hydrogen bond interactions, we followed the dynamics of the excited states involved in the photoexcitation-deexcitation process of this molecule. It remains to be explored the possibility of discriminating between different excited states through the MFPAD profiles. The calculation of MFPADs to differentiate excited states can pave the way to the possibility of a clear discrimination for all the cases where the recognition of excited states is otherwise intricate.

Recently, so-called "nontraditional intrinsic luminescence" has been reported in several macromolecular systems. Although DABCO (1,4-diazabicyclo[2.2.2]octane) is the first system in which the effect was observed, a thorough analysis of the optical properties of the molecule, which would reveal the origin of this mysterious effect, is still pending. We perform an advanced post-Hartree-Fock treatment of the low-lying electronic states of this molecule, which need to be described with care because of their pronounced Rydberg character. We take a deeper look into the low-lying electronic transitions of DABCO targeting the explanation of the complex vibronic structures of its absorption and fluorescence spectra. Two electronic states, the 1 E'(n + 3p xy ) and 1 A 2 (n ″ (n + 3p z ) states, contribute to the absorption spectrum in the 39000-46000 cm -1 spectral range. We also reveal the spectroscopic signature of the 1 (n 2 ″ (n + 3p z (n 1 A 1 ' (n + 3s) state and the electronic ground state. The spectrum, which can be interpreted in terms of populating a few vibrational normal modes, is shifted toward visible wavelengths mostly due to the forced interaction of the lone pair electrons of the two nitrogen atoms. Our work on DABCO may help to rationalize the luminescence of more complex systems containing tertiary amine groups.

Two-dimensional (2D) electronic spectroscopy is a powerful nonlinear technique which provides spectroscopic information on two frequency axes as well as dynamical information as a function of the so-called waiting time. Herein, an ab initio theoretical framework for the simulation of electronic 2D spectra has been developed. The method is based on the classical approximation to the doorway-window representation of three-pulse photon-echo signals and the description of nuclear motion by classical trajectories. Nonadiabatic effects are taken into account by a trajectory surface-hopping algorithm. 2D electronic spectra were simulated with ab initio on-the-fly trajectory calculations using the ADC(2) electronic-structure method for the pyrazine molecule, which is a benchmark system for ultrafast radiationless decay through conical intersections. It is demonstrated that 2D spectroscopy with subfemtosecond UV pulses can provide unprecedented detailed information on the ultrafast photodynamics of polyatomic molecules.

A computational protocol for simulating time-resolved photoelectron signals of medium-sized molecules is presented. The procedure is based on a trajectory surface-hopping description of the excited-state dynamics and a combined Dyson orbital and multicenter B-spline approach for the computation of cross sections and asymmetry parameters. The accuracy of the procedure has been illustrated for the case of ultrafast internal conversion of gas-phase pyrazine excited to the 1 B 2 u (ππ*) state. The simulated spectra and the asymmetry map are compared to the experimental data, and a very good agreement was obtained without applying any energy-dependent rescaling or broadening. An interesting side result of this work is the finding that the signature of the 1 A u ( n π*) state is indistinguishable from that of the 1 B 3 u ( n π*) state in the time-resolved photoelectron spectrum. By locating four symmetrically equivalent minima on the lowest-excited (S 1 ) adiabatic potential energy surface of pyrazine, we revealed the strong vibronic coupling of the 1 A u ( n π*) and 1 B 3 u ( n π*) states near the S 1 ← S 0 band origin.

Despite many studies, the mechanisms of nonradiative relaxation of uracil in the gas phase and in aqueous solution are still not fully resolved. Here we combine theoretical UV absorption spectroscopy with nonadiabatic dynamics simulations to identify the photophysical mechanisms that can give rise to experimentally observed decay time constants. We first compute and theoretically assign the electronic spectra of uracil using the second-order algebraic-diagrammatic-construction (ADC(2)) method. The obtained electronic states, their energy differences and state-specific solvation effects are the prerequisites for understanding the photodynamics. We then use nonadiabatic trajectory-surface-hopping dynamics simulations to investigate the photoinduced dynamics of uracil and uracil-water clusters. In contrast to previous studies, we found that a single mechanism - the ethylenic twist around the C[double bond, length as m-dash]C bond - is responsible for the ultrafast component of the nonradiative decay, both in the gas phase and in solution. Very good agreement with the experimentally determined ultrashort decay time constants is obtained.

Amyloids have unique structural, chemical, and optical properties. Although much theoretical effort has been directed toward understanding amyloid nucleation, the understanding of their optical properties has remained rather limited. In particular, the photophysical mechanisms leading to near-UV excitation and characteristic blue-green luminescence in amyloid systems devoid of aromatic amino acids have not been resolved. We use ab initio static calculations and nonadiabatic dynamics simulations to study the excited electronic states of model amyloid-like peptides. We show that their photophysics is essentially governed by the multitude of nπ* states with excitation localized on the amide groups. The strong stabilization of the nπ* states with respect to the amide group deplanarization and the concomitant increase of the oscillator strength make excitation in the near-UV possible. With respect to emission, our dynamics simulations revealed that the amyloid cross β arrangement effectively hinders the nonradiative relaxation channels usually operative in proteins. Finally, we show that after relaxation of the photoexcited peptides toward the minimum of the different nπ* states, fluorescence takes place in the visible (green) part of the electromagnetic spectrum.

Although diazoalkanes are important carbene precursors in organic synthesis, a comprehensive mechanism of photochemical formation of carbenes from diazoalkanes has not been proposed. Synergies of experiments and computations demonstrate the involvement of higher excited singlet states in the photochemistry of diazoalkanes. In all investigated diazoalkanes, excitation to S 1 results in nonreactive internal conversion to S 0 . On the contrary, excitation to higher-lying singlet states (S n , n > 1) drives the reaction toward a different segment of the S 1 /S 0 conical intersection seam and results in nitrogen elimination and formation of carbenes.

The theoretical assignment of electronic spectra of polyatomic molecules is a challenging problem that requires the specification of the character of a large number of electronic states. We propose a procedure for automatically determining the character of electronic transitions and apply it to the study of UV spectra of DNA bases in the gas phase and in the aqueous environment. The procedure is based on the computation of electronic wave function overlaps and accounts for an extensive sampling of nuclear geometries. Novelties of this work are the theoretical assignment of the electronic spectra of DNA bases up to 190 nm and a state specific analysis of solvation effects. By accounting for different effects contributing to the total solvent shift we obtained a good agreement between the computed and experimental spectra. Effects of vibrational averaging, temperature and solvent-induced structural changes shift excitation energies to lower values. Solvent-solute electrostatic interactions are state specific and strongly destabilize nRyd states, and to lesser extent nπ* and πRyd states. Altogether, this results in the stabilization of ππ* states and destabilization of nπ*, πRyd and nRyd states in solution.

Two algorithms for calculating overlaps between CIS (or TDDFT) type excited state wave functions are presented, one based on an expansion of overlap determinants into level 2 minors (OL2M) and the other based on an expansion of the wave functions into natural transition orbitals (ONTO). Both algorithms are significantly faster than previously available algorithms, with the ONTO algorithm reducing the cost of a single overlap element calculation by a factor of the square of the number of occupied orbitals in the system. The algorithm exhibits orders of magnitude faster calculations for large systems and significantly increases the size of systems for which TDDFT based nonadiabatic dynamics simulations can be performed. The OL2M algorithm is substantially slower for a single overlap matrix element but becomes preferred when overlaps between large numbers of states are required. Additionally, we test the accuracy of approximate overlaps calculated using truncated wave functions and show that truncation can lead to large errors in the overlaps. Lastly, we provide examples of applications for wave function overlaps outside the context of nonadiabatic dynamics.

Trajectory surface hopping (TSH) methods have been widely used to study photoinduced nonadiabatic processes. In the present study, nonadiabatic dynamics simulations with the widely used Tully's fewest switches surface hopping (FSSH) algorithm and a Landau-Zener-type TSH (LZSH) algorithm have been performed for the internal conversion dynamics of pyrazine. The accuracy of the two TSH algorithms has been critically evaluated by a direct comparison with exact quantum dynamics calculations for a model of pyrazine. The model comprises the three lowest excited electronic states (B 3u (nπ * ), A 1u (nπ * ), and B 2u (ππ * )) and the nine most relevant vibrational degrees of freedom. Considering photoexcitation to the diabatic B 2u (ππ * ) state, we examined the time-dependent diabatic and adiabatic electronic population dynamics. It is found that the diabatic populations obtained with both TSH methods are in good agreement with the exact quantum results. Fast population oscillations between the B 3u (nπ * ) and A 1u (nπ * ) states, which reflect nonadiabatic electronic transitions driven by coherent dynamics in the normal mode Q 8a , are qualitatively reproduced by both TSH methods. In addition to the model study, the TSH methods have been interfaced with the second-order algebraic diagrammatic construction ab initio electronic-structure method to perform full-dimensional on-the-fly nonadiabatic dynamics simulations for pyrazine. It is found that the electronic population dynamics obtained with the LZSH method is in excellent agreement with that obtained by the FSSH method using a local diabatization algorithm. Moreover, the electronic populations of the full-dimensional on-the-fly calculations are in excellent agreement with the populations of the three-state nine-mode model, which confirms that the internal conversion dynamics of pyrazine is accurately represented by this reduced-dimensional model on the time scale under consideration (200 fs). The original FSSH method, in which the electronic wave function is propagated in the adiabatic representation, yields less accurate results. The oscillations in the populations of the diabatic B 3u (nπ * ) and A 1u (nπ * ) states driven by the mode Q 8a are also observed in the full-dimensional dynamics simulations.

The computational analysis of the isomer- and conformer-dependent photochemistry of 1- and 2-naphthols and their microsolvated water clusters is motivated by their very different excited state reactivities. We present evidence that 1- and 2-naphthol follow distinct excited state deactivation pathways. The deactivation of 2-naphthols, 2-naphthol water clusters, as well as of the anti conformer of 1-naphthol is mediated by the optically dark 1 πσ* state. The dynamics of the 1 πσ* surface leads to the homolytic cleavage of the OH bond. On the contrary, the excited state deactivation of syn 1-naphthol and 1-naphthol water clusters follows an uncommon reaction pathway. Upon excitation to the bright 1 ππ*(L a ) state, a highly specific excited state hydrogen transfer (ESHT) to carbon atoms C8 and C5 takes place, yielding 1,8- and 1,5-naphthoquinone methides. The ESHT pathway arises from the intrinsic electronic properties of the 1 ππ*(L a ) state of 1-naphthols., (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

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.