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

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 Actinic Photochemistry of 2-Hydroperoxypropanal.

The photochemical reactions triggered by the sunlight absorption of transient volatile organic compounds in the troposphere are notoriously difficult to characterize experimentally due to the unstable and short-lived nature of these organic molecules. Some members of this family of compounds are likely to exhibit a rich photochemistry given the diversity of functional groups they can bear. Even more interesting is the photochemical fate of volatile organic compounds bearing more than one functional group that can absorb light─this is the case, for example, of α-hydroperoxycarbonyls, which are formed during the oxidation of isoprene. Experimental observables characterizing the photochemistry of these molecules like photoabsorption cross-sections or photolysis quantum yields are currently missing, and we propose here to leverage a recently developed computational protocol to predict in silico the photochemical fate of 2-hydroperoxypropanal (2-HPP) in the actinic region. We combine different levels of electronic structure methods─SCS-ADC(2) and XMS-CASPT2─with the nuclear ensemble approach and trajectory surface hopping to understand the mechanistic details of the possible nonradiative processes of 2-HPP. In particular, we predict the photoabsorption cross-section and the wavelength-dependent quantum yields for the observed photolytic pathways and combine them to determine in silico photolysis rate constants. The limitations of our protocol and possible future improvements are discussed.

Calculating Photoabsorption Cross-Sections for Atmospheric Volatile Organic Compounds.

Characterizing the photochemical reactivity of transient volatile organic compounds (VOCs) in our atmosphere begins with a proper understanding of their abilities to absorb sunlight. Unfortunately, the photoabsorption cross-sections for a large number of transient VOCs remain unavailable experimentally due to their short lifetime or high reactivity. While structure-activity relationships (SARs) have been successfully employed to estimate the unknown photoabsorption cross-sections of VOCs, computational photochemistry offers another promising strategy to predict not only the vertical electronic transitions of a given molecule but also the width and shape of the bands forming its absorption spectrum. In this work, we focus on the use of the nuclear ensemble approach (NEA) to determine the photoabsorption cross-section of four exemplary VOCs, namely, acrolein, methylhydroperoxide, 2-hydroperoxy-propanal, and (microsolvated) pyruvic acid. More specifically, we analyze the influence that different strategies for sampling the ground-state nuclear density-Wigner sampling and ab initio molecular dynamics with a quantum thermostat-can have on the simulated absorption spectra. We highlight the potential shortcomings of using uncoupled harmonic modes within Wigner sampling of nuclear density to describe flexible or microsolvated VOCs and some limitations of SARs for multichromophoric VOCs. Our results suggest that the NEA could constitute a powerful tool for the atmospheric community to predict the photoabsorption cross-section for transient VOCs.

Caveat when using ADC(2) for studying the photochemistry of carbonyl-containing molecules.

Several electronic-structure methods are available to study the photochemistry and photophysics of organic molecules. Among them, ADC(2) stands as a sweet spot between computational efficiency and accuracy. As a result, ADC(2) has recently seen its number of applications booming, in particular to unravel the deactivation pathways and photodynamics of organic molecules. Despite this growing success, we demonstrate here that care has to be taken when studying the nonradiative pathways of carbonyl-containing molecules, as ADC(2) appears to suffer from a systematic flaw.

On the Theoretical Determination of Photolysis Properties for Atmospheric Volatile Organic Compounds.

Volatile organic compounds (VOCs) are ubiquitous atmospheric molecules that generate a complex network of chemical reactions in the troposphere, often triggered by the absorption of sunlight. Understanding the VOC composition of the atmosphere relies on our ability to characterize all of their possible reaction pathways. When considering reactions of (transient) VOCs with sunlight, the availability of photolysis rate constants, utilized in general atmospheric models, is often out of experimental reach due to the unstable nature of these molecules. Here, we show how recent advances in computational photochemistry allow us to calculate in silico the different ingredients of a photolysis rate constant, namely, the photoabsorption cross-section and wavelength-dependent quantum yields. The rich photochemistry of tert-butyl hydroperoxide, for which experimental data are available, is employed to test our protocol and highlight the strengths and weaknesses of different levels of electronic structure and nonadiabatic molecular dynamics to study the photochemistry of (transient) VOCs.

Direct Observation of Aggregation-Induced Emission Mechanism.

The mechanism of aggregation-induced emission, which overcomes the common aggregation-caused quenching problem in organic optoelectronics, is revealed by monitoring the real time structural evolution and dynamics of electronic excited state with frequency and polarization resolved ultrafast UV/IR spectroscopy and theoretical calculations. The formation of Woodward-Hoffmann cyclic intermediates upon ultraviolet excitation is observed in dilute solutions of tetraphenylethylene and its derivatives but not in their respective solid. The ultrafast cyclization provides an efficient nonradiative relaxation pathway through crossing a conical intersection. Without such a reaction mechanism, the electronic excitation is preserved in the molecular solids and the molecule fluoresces efficiently, aided by the very slow intermolecular charge and energy transfers due to the well separated molecular packing arrangement. The mechanisms can be general for tuning the properties of chromophores in different phases for various important applications.

Semiclassical Approach to Photophysics Beyond Kasha's Rule and Vibronic Spectroscopy Beyond the Condon Approximation. The Case of Azulene.

Azulene is a prototypical molecule with an anomalous fluorescence from the second excited electronic state, thus violating Kasha's rule, and with an emission spectrum that cannot be understood within the Condon approximation. To better understand the photophysics and spectroscopy of azulene and other nonconventional molecules, we developed a systematic, general, and efficient computational approach combining the semiclassical dynamics of nuclei with ab initio electronic structure. First, to analyze the nonadiabatic effects, we complement the standard population dynamics by a rigorous measure of adiabaticity, estimated with the multiple-surface dephasing representation. Second, we propose a new semiclassical method for simulating non-Condon spectra, which combines the extended thawed Gaussian approximation with the efficient single-Hessian approach. S1 ← S0 and S2 ← S0 absorption and S2 → S0 emission spectra of azulene, recorded in a new set of experiments, agree very well with our calculations. We find that accuracy of the evaluated spectra requires the treatment of anharmonicity, Herzberg-Teller, and mode-mixing effects.

Modulation of charge transfer by N-alkylation to control photoluminescence energy and quantum yield.

Charge transfer in organic fluorophores is a fundamental photophysical process that can be either beneficial, e.g., facilitating thermally activated delayed fluorescence, or detrimental, e.g., mediating emission quenching. N-Alkylation is shown to provide straightforward synthetic control of the charge transfer, emission energy and quantum yield of amine chromophores. We demonstrate this concept using quinine as a model. N-Alkylation causes changes in its emission that mirror those caused by changes in pH (i.e., protonation). Unlike protonation, however, alkylation of quinine's two N sites is performed in a stepwise manner to give kinetically stable species. This kinetic stability allows us to isolate and characterize an N-alkylated analogue of an 'unnatural' protonation state that is quaternized selectively at the less basic site, which is inaccessible using acid. These materials expose (i) the through-space charge-transfer excited state of quinine and (ii) the associated loss pathway, while (iii) developing a simple salt that outperforms quinine sulfate as a quantum yield standard. This N-alkylation approach can be applied broadly in the discovery of emissive materials by tuning charge-transfer states.

Mechanisms of fluorescence quenching in prototypical aggregation-induced emission systems: excited state dynamics with TD-DFTB.

A recent implementation of time-dependent tight-binding density functional theory is employed in excited state molecular dynamics for the investigation of the fluorescence quenching mechanism in 3 prototypical aggregation-induced emission systems. An assessment of the accuracy of the electronic structure method is done by comparison with previous theoretical work while dynamics simulations were extended to the condensed phase to obtain excited state lifetimes comparable to experiment. A thorough investigation is done on tetraphenylethylene in order to resolve the on-going debate on the role of specific deactivation mechanisms. Both gas phase and solvent dynamics were computed for fulvene and silole derivatives.

Infrared Spectroscopy as a Probe of Electronic Energy Transfer.

We have combined electronic and vibrational spectroscopy in a cryogenic ion trap to produce highly resolved, conformer-selective spectra for the ground and excited states of a peptide containing two chromophores. These spectra permit us to determine the precise three-dimensional structure of the peptide and give insight into the migration of the electronic excitation from phenylalanine to tyrosine because changes in the excited-state infrared spectra are sensitive to localization of the electronic energy in each chromophore. The well-controlled experimental conditions make this result a stringent test for theoretical methods dealing with electronic energy transfer.

Photochemistry of 1- and 2-Naphthols and Their Water Clusters: The Role of 1 ππ*(La ) Mediated Hydrogen Transfer to Carbon Atoms.

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 ππ*(La ) 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 ππ*(La ) state of 1-naphthols.

Enhancing the power conversion efficiency of dye-sensitized solar cells via molecular plasmon-like excitations.

We introduce a tactic for employing molecular plasmon-like excitations to enhance solar-to-electric power conversion efficiency of dye-sensitized solar cells. We offer general design principles of dimeric dyes, in which a strong plasmonic interaction between two π-conjugated moieties is promoted. The π-stacked conformations of these dimeric dyes result in a desirable broadened absorption and a longer absorption onset wavelength.

Rationalizing fluorescence quenching in meso-BODIPY dyes.

meso-Substituted boron-dipyrromethene (BODIPY) dyes are a puzzling class of molecules, which feature contrasting emissive behaviors. The full mechanistic picture for these distinctive properties is still missing. Using static and dynamic excited state computations we unravel the key reasons behind these divergences.

Low-Lying ππ* States of Heteroaromatic Molecules: A Challenge for Excited State Methods.

The description of low-lying ππ* states of linear acenes by standard electronic structure methods is known to be challenging. Here, we broaden the framework of this problem by considering a set of fused heteroaromatic rings and demonstrate that standard electronic structure methods do not provide a balanced description of the two (typically) lowest singlet state (La and Lb) excitations. While the Lb state is highly sensitive to correlation effects, La suffers from the same drawbacks as charge transfer excitations. We show that the comparison between CIS/CIS(D) can serve as a diagnostic for detecting the two problematic excited states. Standard TD-DFT and even its spin-flip variant lead to inaccurate excitation energies and interstate gaps, with only a double hybrid functional performing somewhat better. The complication inherent to a balanced description of these states is so important that even CC2 and ADC(2) do not necessarily match the ADC(3) reference.

How does tetraphenylethylene relax from its excited states?

Tetraphenylethylene is a prototypical example of a molecule displaying aggregation-induced emission. Despite many studies on the optical properties of TPE and its derivatives, the origin of the non-emissive behavior in the gas phase or in dilute solutions has yet to be unravelled. Here, we identify the ultrafast deactivation mechanisms responsible for the fluorescence quenching in isolated TPE.

Intramolecular symmetry-adapted perturbation theory with a single-determinant wavefunction.

We introduce an intramolecular energy decomposition scheme for analyzing non-covalent interactions within molecules in the spirit of symmetry-adapted perturbation theory (SAPT). The proposed intra-SAPT approach is based upon the Chemical Hamiltonian of Mayer [Int. J. Quantum Chem. 23(2), 341-363 (1983)] and the recently introduced zeroth-order wavefunction [J. F. Gonthier and C. Corminboeuf, J. Chem. Phys. 140(15), 154107 (2014)]. The scheme decomposes the interaction energy between weakly bound fragments located within the same molecule into physically meaningful components, i.e., electrostatic-exchange, induction, and dispersion. Here, we discuss the key steps of the approach and demonstrate that a single-determinant wavefunction can already deliver a detailed and insightful description of a wide range of intramolecular non-covalent phenomena such as hydrogen bonds, dihydrogen contacts, and π - π stacking interactions. Intra-SAPT is also used to shed the light on competing intra- and intermolecular interactions.

Qualitatively Incorrect Features in the TDDFT Spectrum of Thiophene-Based Compounds.

Ab initio molecular electronic structure computations of thiophene-based compounds constitute an active field of research prompted by the growing interest in low-cost materials for organic electronic devices. In particular, the modeling of electronically excited states and other time-dependent phenomena has moved toward the description of more realistic albeit challenging systems. We demonstrate that due to its underlying approximations, time-dependent density functional theory predicts results that are qualitatively incorrect for thiophene and thienoacenes, although not for oligothiophene chains. The failure includes spurious state inversion and excitation characters, wrong distribution of oscillator strengths and erroneous potential energy surfaces. We briefly analyze possible origins of this behavior and identify alternative methods that alleviate these problems.

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.

Excited state dynamics of thiophene and bithiophene: new insights into theoretically challenging systems.

The computational elucidation and proper description of the ultrafast deactivation mechanisms of simple organic electronic units, such as thiophene and its oligomers, is as challenging as it is contentious. A comprehensive excited state dynamics analysis of these systems utilizing reliable electronic structure approaches is currently lacking, with earlier pictures of the photochemistry of these systems being conceived based upon high-level static computations or lower level dynamic trajectories. Here a detailed surface hopping molecular dynamics of thiophene and bithiophene using the algebraic diagrammatic construction to second order (ADC(2)) method is presented. Our findings illustrate that ring puckering plays an important role in thiophene photochemistry and that the photostability increases when going upon dimerization into bithiophene.

Adjusting the Local Arrangement of π-Stacked Oligothiophenes through Hydrogen Bonds: A Viable Route to Promote Charge Transfer.

We show that substituting quaterthiophene cores with strong H-bond aggregators, such as urea groups, provides an efficient way to adjust the mutual in-plane displacements of the semiconducting units and promote charge transfer. Our 2-D structure-property mapping reveals that the insertion of substituents induces up to 2.0 Å longitudinal and transversal displacements between the π-conjugated moieties. Some of these relative displacements lead to improved cofacial orbital overlaps that are otherwise inaccessible due to Pauli repulsion. Our results also emphasize that the fine-tuning of in-plane displacements is more effective than achieving "tighter" packing to promote charge-transfer properties.