On the origin of the shift between vertical excitation and band maximum in molecular photoabsorption

Evaluating the interactions between vibrational modes and electronic transitions using frontier orbital energy derivatives

Vibrations affect molecular optoelectronic properties, even at zero kelvin. Accounting for these effects using computational modelling is costly, as it requires many calculations at geometries distorted from equilibrium. Here, we propose a low-cost method for identifying vibrations most strongly coupled to the electronic structure, based on using orbital energy derivatives as a diagnostic.

A DFT/MRCI Hamiltonian parameterized using only ab initio data: I. valence excited states

A new combined density functional theory and multi-reference configuration interaction (DFT/MRCI) Hamiltonian parameterized solely using the benchmark ab initio vertical excitation energies obtained from the QUEST databases is presented. This new formulation differs from all previous versions of the method in that the choice of the underlying exchange-correlation (XC) functional employed to construct the one-particle (orbital) basis is considered, and a new XC functional, QTP17, is chosen for its ability to generate a balanced description of core and valence vertical excitation energies. The ability of the new DFT/MRCI Hamiltonian, termed QE8, to furnish accurate excitation energies is confirmed using benchmark quantum chemistry computations, and a mean absolute error of 0.16 eV is determined for the wide range of electronic excitations included in the validation dataset. In particular, the QE8 Hamiltonian dramatically improves the performance of DFT/MRCI for doubly excited states. The performance of fast approximate DFT/MRCI methods, p-DFT/MRCI and DFT/MRCI(2), is also evaluated using the QE8 Hamiltonian, and they are found to yield excitation energies in quantitative agreement with the parent DFT/MRCI method, with the two methods exhibiting a mean difference of 0.01 eV with respect to DFT/MRCI over the entire benchmark set.

2-in-1 Phase Space Sampling for Calculating the Absorption Spectrum of the Hydrated Electron

The investigation of vibrational effects on absorption spectrum calculations often employs Wigner sampling or thermal sampling. While Wigner sampling incorporates zero-point energy, it may not be suitable for flexible systems. Thermal sampling is applicable to anharmonic systems yet treats nuclei classically. The application of generalized smoothed trajectory analysis (GSTA) as a postprocessing method allows for the incorporation of nuclear quantum effects (NQEs), combining the advantages of both sampling methods. We demonstrate this approach in computing the absorption spectrum of a hydrated electron. Theoretical exploration of the hydrated electron and its embryonic forms, such as water cluster anions, poses a significant challenge due to the diffusivity of the excess electron and the continuous motion of water molecules. In many previous studies, the wave nature of atomic nuclei is often neglected, despite the substantial impact of NQEs on thermodynamic and spectroscopic properties, particularly for hydrogen atoms. In our studies, we examine these NQEs for the excess electrons in various water systems. We obtained structures from mixed classical-quantum simulations for water cluster anions and the hydrated electron by incorporating the quantum effects of atomic nuclei with the filtration of the classical trajectories. Absorption spectra were determined at different theoretical levels. Our results indicate significant NQEs, red shift, and broadening of the spectra for hydrated electron systems. This study demonstrates the applicability of GSTA to complex systems, providing insights into NQEs on energetic and structural properties.

Excited state properties of an A-D-A non-fullerene electron acceptor: a LC-TD-DFTB study

Understanding charge transfer processes is essential to estimate the performance of organic photovoltaic technologies. Although experimental production is on the rise, predictability strongly relies on theoretical modeling, which is limited to the size of semiconductors. As a computationally favorable approach, we benchmarked the long-range corrected (LC) time-dependent (TD) formulation of the semi-empirical density functional-based tight-binding method (DFTB) for three polycyclic aromatic hydrocarbons (PAHs) and studied the DTP-IC-4Ph molecule, a PAH-based non-fullerene electron acceptor (NFA) with an A-D-A backbone structure. After a thorough investigation into the long-range parameter (ω) tuning for naphthalene, anthracene and pyrene, the excitation energies, oscillator strengths and Natural Transition Orbitals (NTOs) were compared with the standard ωB97X-D/6-31G(d,p) level of theory and the ADC2/6-31G(d,p) multiconfigurational method. We estimated mobility-related properties of the NFA and considered 1000 thermally accessible configurations to qualitatively reproduce the experimental absorption profile and investigate the energetic disorder. Finally, we conducted a fragment-based analysis using the one-electron transition density matrix (1TDM) to determine the character of the excited states and investigate the effect of side chains on exciton formation. Our results are sensitive to the level of theory and highly dependent on the long-range parameter but suggest that the presence of alkyl chains promotes a higher average charge delocalization and allows for additional hopping mechanisms, favoring the charge transfer dynamics.

Photodissociation and formation of an ion-pair in CH2 FCl (HCFC-31)

Although CH2 FCl (HCFC-31) recently became of great atmospheric importance, studies concerning its excited states are almost nonexistent. Several excited singlet states were studied (valence nσ* and Rydberg n3s, n3p, σ3s, and σ3p) through highly correlated multireference configuration interaction with singles and doubles, including extensivity correction. Comparison with the states of CH3 Cl indicates a strong influence of the F atom. Potential energy curves suggest formation of an electrostatically bound complex that relaxes to a hydrogen-bonded contact ion-pair (HBCIP) which can decay yielding CH2 F + Cl or to the ground state minimum of CH2 FCl. The HBCIP has a dipole moment of 9.57 D, a CI wavefunction described as 0.65ionic + 0.20biradical and it is strongly bonded by 4.72 eV. Its H bond has characteristics of moderate and strong H bonds. The simulated absorption spectrum confirms the nσ* assignment for the first and suggests the n3s + n3pσ assignment for the second band.

Role of Singles Amplitudes in ADC(2) and CC2 for Low-Lying Electronically Excited States

The closely related second-order methods CC2 and ADC(2) usually perform very similarly for single excitations of organic molecules. However, as rationalized in this work, significant deviations between these two methods can arise if the ground state and a low-lying singly excited state arise from a strong coupling between their leading configurations. Such a configuration mixing is partially accounted for in CC2 through the ground-state singles amplitudes but is omitted in ADC(2). This can cause unusual deviations between the results obtained with these methods. In this work, we study how severe this effect can become at the example of two solvatochromic dyes: the negatively solvatochromic betaine dye N1-tBu and the positively solvatochromic bithiophene P1. These two dyes allow one to study the limits of both small and somewhat larger excitation energies and configuration mixing by tuning the S0 → S1 transition energy through the polarity of the environment. Higher-level calculations at the CC3 level provide information on the accuracy of ADC(2) and CC2 in these cases. The most extreme deviation between ADC(2) and CC2 is found for N1-tBu in vacuum, where the ADC(2) result is 0.45 eV below that of CC2. In this case, the methodical error of CC2 with respect to CC3 is only 0.05 eV. With increasing excitation energy in polar solvents, the CC2-ADC(2) deviation decreases and reaches a value of only 0.15 eV. For P1, which has larger excitation energies, these effects are reversed due to the opposite solvatochromism but also smaller in magnitude: the deviation increases from 0.08 eV in vacuum to 0.16 eV in the so-called conductor limit of the continuum solvation model. Although for these two dyes larger deviations are observed for smaller excitation energies, the extent of configuration mixing does not generally correlate with only the size of excitation energy. For example, s-triazine (0.15 eV), formamide (0.19 eV), and formaldehyde (0.23 eV) also show large deviations between CC2 and ADC(2) despite their much higher excitation energies compared to those of N1-tBu and P1.

Benchmarks for transition metal spin-state energetics: why and how to employ experimental reference data?

Accurate prediction of energy differences between alternative spin states of transition metal complexes is essential in computational (bio)inorganic chemistry-for example, in characterization of spin crossover materials and in the theoretical modeling of open-shell reaction mechanisms-but it remains one of the most compelling problems for quantum chemistry methods. A part of this challenge is to obtain reliable reference data for benchmark studies, as even the highest-level applicable methods are known to give divergent results. This Perspective discusses two possible approaches to method benchmarking for spin-state energetics: using either theoretically computed or experiment-derived reference data. With the focus on the latter approach, an extensive general review is provided for the available experimental data of spin-state energetics and their interpretations in the context of benchmark studies, targeting the possibility of back-correcting the vibrational effects and the influence of solvents or crystalline environments. With a growing amount of experience, these effects can be now not only qualitatively understood, but also quantitatively modeled, providing the way to derive nearly chemically accurate estimates of the electronic spin-state gaps to be used as benchmarks and advancing our understanding of the phenomena related to spin states in condensed phases.

First hyperpolarizability of the di-8-ANEPPS and DR1 nonlinear optical chromophores in solution. An experimental and multi-scale theoretical chemistry study

The solvent effects on the linear and second-order nonlinear optical properties of an aminonaphtylethenylpyridinium (ANEP) dye are investigated by combining experimental and theoretical chemistry methods. On the one hand, deep near infrared (NIR) hyper-Rayleigh scattering (HRS) measurements (1840-1950 nm) are performed on solutions of di-8-ANEPPS in deuterated chloroform, dimethylformamide, and dimethylsulfoxide to determine their first hyperpolarizablity (βHRS). For the first time, these HRS experiments are carried out in the picosecond regime in the deep NIR with very moderate (≤3 mW) average input power, providing a good signal-to-noise ratio and avoiding solvent thermal effects. Moreover, the frequency dispersion of βHRS is investigated for Disperse Red 1 (DR1), a dye commonly used as HRS external reference. On the other hand, these are compared with computational chemistry results obtained by using a sequential molecular dynamics (MD) then quantum mechanics (QM) approach. The MD method allows accounting for the dynamical nature of the molecular structures. Then, the QM part is based on TDDFT/M06-2X/6-311+G* calculations using solvation models ranging from continuum to discrete ones. Measurements report a decrease of the βHRS of di-8-ANEPPS in more polar solvents and these effects are reproduced by the different solvation models. For di-8-ANEPPS and DR1, comparisons show that the use of a hybrid solvation model, combining the description of the solvent molecules around the probe by point charges with a continuum model, already achieves quasi quantitative agreement with experiment. These results are further improved by using a polarizable embedding that includes the atomic polarizabilities in the solvent description.

Tautomeric equilibrium and spectroscopic properties of 8-azaguanine revealed by quantum chemistry methods

8-azaguanine is a triazolopyrimidine nucleobase analog possessing potent antibacterial and antitumor activities, and it has been implicated as a lead molecule in cancer and malaria therapy. Its intrinsic fluorescence properties can be utilized for monitoring its interactions with biological polymers like proteins or nucleic acids. In order to better understand these interactions, it is important to know the tautomeric equilibrium of this compound. In this work, the tautomeric equilibrium of all natural neutral and anionic compound forms (except highly improbable imino-enol tautomers) as well as their methyl derivatives and ribosides was revealed by quantum chemistry methods. It was shown that, as expected, tautomers protonated at positions 1 and 9 dominate neutral forms both in gas phase and in aqueous solution. 8-azaguanines methylated at any position of the triazole ring are protonated at position 1. The computed vertical absorption and emission energies are in very good agreement with the experimental data. They confirm the validity of the assumption that replacing the proton with the methyl group does not significantly change the positions of absorption and fluorescence peaks.

Efficient Protocol for Computing MCD Spectra in a Broad Frequency Range Combining Resonant and Damped CC2 Quadratic Response Theory

Coupled cluster response theory offers a path to high-accuracy calculations of spectroscopic properties, such as magnetic circular dichroism (MCD). However, divergence or slow convergence issues are often encountered for electronic transitions in high-energy regions with a high density of states. This is here addressed for MCD by an implementation of damped quadratic response theory for resolution-of-identity coupled cluster singles-and-approximate-doubles (RI-CC2), along with an implementation of the MCD A term from resonant response theory. Combined, damped and resonant response theory calculations provide an efficient strategy to obtain MCD spectra over a broad frequency range and for systems that include highly symmetric molecules with degenerate excited states. The protocol is illustrated by application to zinc tetrabenzoporphyrin in the energy region of 2-8 eV and comparison to experimental data. Timings are reported for the resonant and damped approaches, showing that a greater part of the calculation time is consumed by the construction of the building blocks for the final MCD ellipticity. A recommendation on how to use the procedure is outlined.

Theoretical investigation on the functional group modulation of UV-Vis absorption profiles of triphenylamine derivatives

Understanding the functional group modulation of electronic structure and excitation is pivotal to the design of organic small molecules (OSMs) for photoelectric applications. In this study, we employed density functional theory (DFT) and time-dependent DFT (TDDFT) calculations to explore the unique absorption character of four triphenylamine photosensitizers. The various conformations were investigated given the multiple single bonds in the compounds, and the resemblance in the electronic structure of different conformations is affirmed because the coplanarity and consequent long-range conjugation is maintained regardless of the orientation of the flexible blocks. Six functionals were evaluated, and MN15 was found to successfully reproduce the intense secondary absorption peak for the double 3,4-ethylenedioxythiophene (EDOT) modified sensitizer over B3LYP, PBE0, M062X, CAM-B3LYP, and ωB97XD. The introduction of EDOT gives rise to a new excited state S4, which is a local excitation constrained in the EDOT substituent triphenylamine block. This new excited state S4, in combination with inherent S2 and S3 derived from prototype molecule TPA-Pyc, jointly contributes to the hump of the secondary absorption peak of ETE-Pyc and finally affects the light-harvesting ability of the dye-sensitized TiO2 photoanode. The current findings provide guidance toward the rational design of OSMs with good light-harvest ability.

Simulation of the VUV Absorption Spectra of Oxygenates and Hydrocarbons: A Joint Theoretical-Experimental Study

Vacuum UV absorption spectroscopy is regularly used to provide unambiguous identification of a target species, insight into the electronic structure of molecules, and quantitative species concentrations. As molecules of interest have become more complex, theoretical spectra have been used in tandem with laboratory spectroscopic analysis or as a replacement when experimental data is unavailable. However, it is difficult to determine which theoretical methodologies can best simulate experiment. This study examined the performance of EOM-CCSD and 10 TD-DFT functionals (B3LYP, BH&HLYP, BMK, CAM-B3LYP, HSE, M06-2X, M11, PBE0, ωB97X-D, and X3LYP) to produce reliable vacuum UV absorption spectra for 19 small oxygenates and hydrocarbons using vertical excitation energies. The simulated spectra were analyzed against experiment using both a qualitative analysis and quantitative metrics, including cosine similarity, relative integral change, mean signed error, and mean absolute error. Based on our ranking system, it was determined that M06-2X was consistently the top performing TD-DFT method with BMK, CAM-B3LYP, and ωB97X-D also producing reliable spectra for these small combustion species.

Electronic Spectroscopy and Dissociation Dynamics of Vinyl-Substituted Criegee Intermediates: 2-Butenal Oxide and Comparison with Methyl Vinyl Ketone Oxide and Methacrolein Oxide Isomers

The 2-butenal oxide Criegee intermediate [(CH₃CH═CH)CHOO], an isomer of the four-carbon unsaturated Criegee intermediates derived from isoprene ozonolysis, is characterized on its first π* ← π electronic transition and by the resultant dissociation dynamics to O (¹D) + 2-butenal [(CH₃CH═CH)CHO] products. The electronic spectrum of 2-butenal oxide under jet-cooled conditions is observed to be broad and unstructured with peak absorption at 373 nm, spanning to half maxima at 320 and 420 nm, and in good accord with the computed vertical excitation energies and absorption spectra obtained for its lowest energy conformers. The distribution of total kinetic energy released to products is ascertained through velocity map imaging of the O (¹D) products. About half of the available energy, deduced from the theoretically computed asymptotic energy, is accommodated as internal excitation of the 2-butenal fragment. A reduced impulsive model is introduced to interpret the photodissociation dynamics, which accounts for the geometric changes between 2-butenal oxide and the 2-butenal fragment, and vibrational activation of associated modes in the 2-butenal product. Application of the reduced impulsive model to the photodissociation of isomeric methyl vinyl ketone oxide reveals greater internal activation of the methyl vinyl ketone product arising from methyl internal rotation and rock, which is distinctly different from the dissociation dynamics of 2-butenal oxide or methacrolein oxide.

Revealing and Tuning the Photophysics of C=N Containing Photothermal Molecules: Excited State Dynamics Simulations

Molecular photothermal conversion materials are recently attracting increasing attention for phototherapy applications. Herein we investigate the excitation and de-excitation processes of a photothermal molecule (C1TI) that is among the recently developed class of small-molecule-based photothermal imines with superb photothermal conversion efficiencies (PTCEs) up to 90% and a molecule (M2) that is constructed by replacing the amino group of C1TI with an H atom, via excited-state dynamics simulations based on the time-dependent density functional theory (TD-DFT). The simulations reveal fast (<150 fs of average time) nonradiative decays of the lowest excited singlet (S1) state to a conical intersection (CI) with the ground (S0) state in high yields (C1TI: 93.9% and M2: 87.1%). The fast decays, driven by C=N bond rotation to a perpendicular structural configuration, are found to be barrierless. The slight structural difference between C1TI and M2 leads to drastically different S0-S1 energy surfaces, especially M2 features a relatively much lower CI (0.8 eV in energy) and much more decay energy (1.0 eV) to approach the CI. This work provides insights into the de-excitation mechanisms and the performance tuning of C=N enabled photothermal materials.

Pre-Dewar structure modulates protonated azaindole photodynamics

Recent experimental work revealed that the lifetime of the S₃ state of protonated 7-azaindole is about ten times longer than that of protonated 6-azaindole. We simulated the nonradiative decay pathways of these molecules using trajectory surface hopping dynamics after photoexcitation into S₃ to elucidate the reason for this difference. Both isomers mainly follow a common ππ* relaxation pathway involving multiple state crossings while coming down from S₃ to S₁ in the subpicosecond time scale. However, the simulations reveal that the excited-state topographies are such that while the 6-isomer can easily access the region of nonadiabatic transitions, the internal conversion of the 7-isomer is delayed by a pre-Dewar bond formation with a boat conformation.

Calculation of Metallocene Ionization Potentials via Auxiliary Field Quantum Monte Carlo: Toward Benchmark Quantum Chemistry for Transition Metals

The accurate ab initio prediction of ionization energies is essential to understanding the electrochemistry of transition metal complexes in both materials science and biological applications. However, such predictions have been complicated by the scarcity of gas phase experimental data, the relatively large size of the relevant molecules, and the presence of strong electron correlation effects. In this work, we apply all-electron phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) utilizing multideterminant trial wave functions to six metallocene complexes to compare the computed adiabatic and vertical ionization energies with experimental results. We find that ph-AFQMC yields mean absolute errors (MAEs) of 1.69 ± 1.02 kcal/mol for the adiabatic energies and 2.85 ± 1.13 kcal/mol for the vertical energies. We also carry out density functional theory (DFT) calculations using a variety of functionals, which yields MAEs of 3.62-6.98 kcal/mol and 3.31-9.88 kcal/mol, as well as one variant of localized coupled cluster calculations (DLPNO-CCSD(T₀) with moderate PNO cutoffs), which has MAEs of 4.96 and 6.08 kcal/mol, respectively. We also test the reliability of DLPNO-CCSD(T₀) and DFT on acetylacetonate (acac) complexes for adiabatic energies measured in the same manner experimentally, and we find higher MAEs, ranging from 4.56 to 10.99 kcal/mol (with a different ordering) for DFT and 6.97 kcal/mol for DLPNO-CCSD(T₀). Finally, by utilizing experimental solvation energies, we show that accurate reduction potentials in solution for the metallocene series can be obtained from the AFQMC gas phase results.

Cyclo[18]carbon including zero-point motion: ground state, first singlet and triplet excitations, and hole transfer

Recent synthesis of cyclo[18]carbon has spurred increasing interest in carbon rings. We focus on a comparative inspection of ground and excited states, as well as of hole transfer properties of cumulenic and polyynic cyclo[18]carbon via Density Functional Theory (DFT), time-dependent DFT (TD-DFT) and real-time time-dependent DFT (RT-TDDFT). Zero-point vibrations are also accounted for, using a Monte Carlo sampling technique and a less exact, yet mode-resolved, quadratic approximation. The inclusion of zero-point vibrations leads to a red-shift on the HOMO-LUMO gap and the first singlet and triplet excitation energies of both conformations, correcting the values of the 'static' configurations by 9% to 24%. Next, we oxidize the molecule, creating a hole at one carbon atom. Hole transfer along polyynic cyclo[18]carbon is decreased in magnitude compared to its cumulenic counterpart and lacks the symmetric features the latter displays. Contributions by each mode to energy changes and hole transfer between diametrically opposed atoms vary, with specific bond-stretching modes being dominant.

Excited state absorption of DNA bases in the gas phase and in chloroform solution: a comparative quantum mechanical study

We study the excited state absorption (ESA) properties of the four DNA bases (thymine, cytosine, adenine, and guanine) by different single reference quantum mechanical methods, namely, equation of motion coupled cluster singles and doubles (EOM-CCSD), singles, doubles and perturbative triples (EOM-CC3), and time-dependent density functional theory (TD-DFT), with the long-range corrected CAM-B3LYP functional. Preliminary results at the Tamm-Dancoff (TDA) CAM-B3LYP level using the maximum overlap method (MOM) are reported for thymine. In the gas phase, the three methods predict similar One Photon Absorption (OPA) spectra, which are consistent with the experimental results and with the most accurate computational studies available in the literature. The ESA spectra are then computed for the ππ* states (one for pyrimidine, two for purines) associated with the lowest-energy absorption band, and for the close-lying nπ* state. The EOM-CC3, EOM-CCSD and CAM-B3LYP methods provide similar ESA spectral patterns, which are also in qualitative agreement with literature RASPT2 results. Once validated in the gas phase, TD-CAM-B3LYP has been used to compute the ESA in chloroform, including solvent effects by the polarizable continuum model (PCM). The predicted OPA and ESA spectra in chloroform are very similar to those in the gas phase, most of the bands shifting by less than 0.1 eV, with a small increase of the intensities and a moderate destabilization of the nπ* state. Finally, ESA spectra have been computed from the minima of the lowest energy ππ* state, and found in line with the available experimental transient absorption spectra of the nucleosides in solution, providing further validation of our computational approach.

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.

Assessment of time-dependent density functionals for the electronic excitation energies of organic dyes used in DSSCs

The absorption spectra modeled as the vertical excitation energies of 13 dye sensitizers used in dye-sensitized solar cells (DSSCs) are benchmarked by means of time-dependent (TD)-DFT, using 36 functionals from different DFT rungs.

Active Learning Configuration Interaction for Excited-State Calculations of Polycyclic Aromatic Hydrocarbons

We present the active learning configuration interaction (ALCI) method for multiconfigurational calculations based on large active spaces. ALCI leverages the use of an active learning procedure to find important electronic configurations among the full configurational space generated within an active space. We tested it for the calculation of singlet-singlet excited states of acenes and pyrene using different machine learning algorithms. The ALCI method yields excitation energies within 0.2-0.3 eV from those obtained by traditional complete active-space configuration interaction (CASCI) calculations (affordable for active spaces up to 16 electrons in 16 orbitals) by including only a small fraction of the CASCI configuration space in the calculations. For larger active spaces (we tested up to 26 electrons in 26 orbitals), not affordable with traditional CI methods, ALCI captures the trends of experimental excitation energies. Overall, ALCI provides satisfactory approximations to large active-space wave functions with up to 10 orders of magnitude fewer determinants for the systems presented here. These ALCI wave functions are promising and affordable starting points for the subsequent second-order perturbation theory or pair-density functional theory calculations.

Multiscale Modeling of Electronic Spectra Including Nuclear Quantum Effects

Theoretical prediction of electronic absorption spectra without input from experiments is no easy feat, as it requires addressing all of the factors that affect line shapes. In practice, however, the methodologies are limited to treat these ingredients only to a certain extent. Here, we present a multiscale protocol that addresses the temperature, solvent, and nuclear quantum effects as well as anharmonicity and the reconstruction of the final spectra from individual transitions. First, quantum mechanics/molecular mechanics (QM/MM) molecular dynamics is conducted to obtain trajectories of solute-solvent configurations, from which the corresponding quantum-corrected ensembles are generated through the generalized smoothed trajectory analysis (GSTA). The optical spectra of the ensembles are then produced by calculating vertical transitions using time-dependent density-functional theory (TDDFT) with implicit solvation. To obtain the final spectral shapes, the stick spectra from TDDFT are convoluted with Gaussian kernels where the half-widths are determined by a statistically motivated strategy. We have tested our method by calculating the UV-vis spectra of a recently discovered acridine photocatalyst in two redox states. Vibronic progressions and broadenings due to the finite lifetime of the excited states are not included in the methodology yet. Nuclear quantization affects the relative peak intensities and widths, which is necessary to reproduce the experimental spectrum. We have also found that using only the optimized geometry of each molecule works surprisingly well if a proper empirical broadening factor is applied. This is explained by the rigidity of the conjugated chromophore moieties of the selected molecules, which are mainly responsible for the excitations in the spectra. In contrast, we have also shown that other parts of the molecules are flexible enough to feature anharmonicities that impair the use of other techniques such as Wigner sampling.

The importance of finite temperature and vibrational sampling in the absorption spectrum of a nitro-functionalized Ru(ii) water oxidation catalyst

Consideration of finite temperature and vibrational motion can be an essential component for accurate simulations of absorption spectra. Here we use finite-temperature Wigner phase-space sampling to investigate the intense absorption of the water oxidation catalyst Ru(dppip-NO2) in the visible (vis) region. The influence of vibrational and torsional motions as well as temperature effects are addressed for the different protonation forms of the pH-sensitive dppip-NO2 ligand of the catalyst. Excitations to the nitrophenyl group and π-system of dppip-NO2, which characterize the absorption band in the equilibrium spectra, experience energy shifts and a significant decrease in oscillator strength when nuclear motion is considered. The importance of excitations to the nitrophenyl group for the vis band is reduced in the spectra computed from the 300 K ensembles, which feature broad distributions of the corresponding dihedral angles. The effects of vibrational sampling on the absorption spectra may be attributed to nitrophenyl and, in particular, to NO2 torsional motions. We expect finite temperature and vibrational sampling to be important for simulating the absorption spectra of other transition metal complexes with flexible ligands or nitro-aromatic motifs.

The Quest to Simulate Excited-State Dynamics of Transition Metal Complexes

This Perspective describes current computational efforts in the field of simulating photodynamics of transition metal complexes. We present the typical workflows and feature the strengths and limitations of the different contemporary approaches. From electronic structure methods suitable to describe transition metal complexes to approaches able to simulate their nuclear dynamics under the effect of light, we give particular attention to build a bridge between theory and experiment by critically discussing the different models commonly adopted in the interpretation of spectroscopic experiments and the simulation of particular observables. Thereby, we review all the studies of excited-state dynamics on transition metal complexes, both in gas phase and in solution from reduced to full dimensionality.

Systematic improvement of molecular excited state calculations by inclusion of nuclear quantum motion: A mode-resolved picture and the effect of molecular size

The energies of molecular excited states arise as solutions to the electronic Schrödinger equation and are often compared to experiment. At the same time, nuclear quantum motion is known to be important and to induce a redshift of excited state energies. However, it is thus far unclear whether incorporating nuclear quantum motion in molecular excited state calculations leads to a systematic improvement of their predictive accuracy, making further investigation necessary. Here, we present such an investigation by employing two first-principles methods for capturing the effect of quantum fluctuations on excited state energies, which we apply to the Thiel set of organic molecules. We show that accounting for zero-point motion leads to much improved agreement with experiment, compared to "static" calculations that only account for electronic effects, and the magnitude of the redshift can become as large as 1.36 eV. Moreover, we show that the effect of nuclear quantum motion on excited state energies largely depends on the molecular size, with smaller molecules exhibiting larger redshifts. Our methodology also makes it possible to analyze the contribution of individual vibrational normal modes to the redshift of excited state energies, and in several molecules, we identify a limited number of modes dominating this effect. Overall, our study provides a foundation for systematically quantifying the shift of excited state energies due to nuclear quantum motion and for understanding this effect at a microscopic level.

Excited States of Bromopyrimidines Probed by VUV Photoabsorption Spectroscopy and Theoretical Calculations

We report absolute photoabsorption cross sections for gas-phase 2- and 5-bromopyrimidine in the 3.7-10.8 eV energy range, in a joint theoretical and experimental study. The measurements were carried out using high-resolution vacuum ultraviolet synchrotron radiation, with quantum chemical calculations performed through the nuclear ensemble approach in combination with time-dependent density functional theory, along with additional Franck-Condon Herzberg-Teller calculations for the first absorption band (3.7-4.6 eV). The cross sections of both bromopyrimidines are very similar below 7.3 eV, deviating more substantially from each other at higher energies. In the 7.3-9.0 eV range where the maximum cross-section is found, a single and broad band is observed for 5-bromopyrimidine, while more discernible features appear in the case of 2-bromopyrimidine. Several π* ← π transitions account for the most intense bands, while weaker ones are assigned to transitions involving the nitrogen and bromine lone pairs, the antibonding σ*Br orbital, and the lower-lying Rydberg states. A detailed comparison with the available photo-absorption data of bromobenzene is also reported. We have found significant differences regarding the main absorption band, which is more peaked in bromobenzene, becoming broader and shifting to higher energies in both bromopyrimidines. In addition, there is a significant suppression of vibrational structures and of Rydberg states in the pair of isomers, most noticeably for 2-bromopyrimidine.

MLatom 2: An Integrative Platform for Atomistic Machine Learning

Atomistic machine learning (AML) simulations are used in chemistry at an ever-increasing pace. A large number of AML models has been developed, but their implementations are scattered among different packages, each with its own conventions for input and output. Thus, here we give an overview of our MLatom 2 software package, which provides an integrative platform for a wide variety of AML simulations by implementing from scratch and interfacing existing software for a range of state-of-the-art models. These include kernel method-based model types such as KREG (native implementation), sGDML, and GAP-SOAP as well as neural-network-based model types such as ANI, DeepPot-SE, and PhysNet. The theoretical foundations behind these methods are overviewed too. The modular structure of MLatom allows for easy extension to more AML model types. MLatom 2 also has many other capabilities useful for AML simulations, such as the support of custom descriptors, farthest-point and structure-based sampling, hyperparameter optimization, model evaluation, and automatic learning curve generation. It can also be used for such multi-step tasks as Δ-learning, self-correction approaches, and absorption spectrum simulation within the machine-learning nuclear-ensemble approach. Several of these MLatom 2 capabilities are showcased in application examples.

Molecular excited states through a machine learning lens

Theoretical simulations of electronic excitations and associated processes in molecules are indispensable for fundamental research and technological innovations. However, such simulations are notoriously challenging to perform with quantum mechanical methods. Advances in machine learning open many new avenues for assisting molecular excited-state simulations. In this Review, we track such progress, assess the current state of the art and highlight the critical issues to solve in the future. We overview a broad range of machine learning applications in excited-state research, which include the prediction of molecular properties, improvements of quantum mechanical methods for the calculations of excited-state properties and the search for new materials. Machine learning approaches can help us understand hidden factors that influence photo-processes, leading to a better control of such processes and new rules for the design of materials for optoelectronic applications.

Describing ground and excited state potential energy surfaces for molecular photoswitches using coupled cluster models

In this article, we use two extensively studied systems, a retinal model system and azobenzene, to explore the use of coupled cluster models for describing ground and singlet excited state potential energy surfaces of photoswitchable systems. While not being suitable for describing nuclear dynamics of photoisomerization, coupled cluster models have useful attributes, such as the inclusion of dynamical correlation, their black box nature, and the systematic improvement offered by truncation level. Results for the studied systems show that when triple excitations (here through the CC3 model) are included, ground and excited state potential energy surfaces for isomerization paths may reliably be generated, also for states of doubly excited character. For ground state equilibrium cis- and trans-azobenzene, the molecular geometry and basis set is seen to significantly impact the vertical excitation energies for the two lowest excited states. Efficient implementations of coupled cluster models can therefore constitute valuable tools for investigating photoswitchable systems and can be used for preliminary black box studies to gather information before more complicated excited state dynamics approaches are pursued.

Mechanistic insight into the fluorescence activity of forensic fingerprinting reagents

Fingerprint detection is still the primary investigative technique for deciphering criminal inquiries and identifying individuals. The main forensic fingerprinting reagents (FFRs) currently in use can require multiple treatment steps to produce fingerprints of sufficient quality. Therefore, the development of new, more effective FFRs that require minimal chemical treatment is of great interest in forensic chemistry. In this work, prudently crafted density functional theory and time-dependent density functional theory calculations are utilized to derive mechanistic insight into the optical activity of the non-fluorescent product of ninhydrin, diketohydrindylidenediketohydrindamine (DYDA), and fluorescent product of DFO (1,8-diazafluoren-9-one). We investigate various protonation sites to gain an understanding of isomeric preference in the solid-state material. A relaxed scan of a single torsion angle rotation in the S₁ minimized geometry of the O-protonated DYDA isomer suggests a conical intersection upon ∼10° rotation. We show that the absence of a rigid hydrogen-bonded network in the crystal structure of DYDA supports the hypothesis of torsion rotation, which leads de-excitation to occur readily. Conversely, for the fluorescent DFO product, our calculations support an avoided crossing suggestive of a non-radiative mechanism when the torsion angle is rotated by about ∼100°. This mechanistic insight concurs with experimental observations of fluorescence activity in DFO and may aid the photophysical understanding of poorly visualized fingerprints due to weak fluorescence. We show that identifying suggestive avoided crossings via the method described here can be used to initialize thoughts toward the computational design of FFRs.

Perfluoro effect on the electronic excited states of para-benzoquinone revealed by experiment and theory

We report a comprehensive study on the electronic excited states of tetrafluoro-1,4-benzoquinone, through high-resolution vacuum ultraviolet photoabsorption spectroscopy and time-dependent density functional theory calculations performed within the nuclear ensemble approach. Absolute cross section values were experimentally determined in the 3.8-10.8 eV energy range. The present experimental results represent the highest resolution data yet reported for this molecule and reveal previously unresolved spectral structures. The interpretation of the results was made in close comparison with the available data for para-benzoquinone [Jones et al., J. Chem. Phys., 2017, 146, 184303]. While the dominant absorption features for both molecules arise from analogous π* ← π transitions, some remarkable differences have been identified. The perfluoro effect manifests in different ways: shifts in band positions and cross sections, appearance of features associated with excitations to σCF* orbitals, and spectrum broadening by quenching of either vibrational or Rydberg progressions. The level of agreement between experiment and theory is very satisfactory, yet that required the inclusion of nuclear quantum effects in the calculations. We have also discussed the role of temperature on the absorption spectrum, as well as the involvement of core-excited resonances in promoting dissociative electron attachment reactions in the 3-5 eV range.

Spin-state energetics of metallocenes: How do best wave function and density functional theory results compare with the experimental data?

We benchmark the accuracy of quantum-chemical methods, including wave function theory methods [coupled cluster theory at the CCSD(T) level, multiconfigurational perturbation-theory (CASPT2, NEVPT2) and internally contracted multireference configuration interaction (MRCI)] and 30 density functional theory (DFT) approximations, in reproducing the spin-state splittings of metallocenes. The reference values of the electronic energy differences are derived from the experimental spin-crossover enthalpy for manganocene and the spectral data of singlet-triplet transitions for ruthenocene, ferrocene, and cobaltocenium. For ferrocene and cobaltocenium we revise the previous experimental interpretations regarding the lowest triplet energy; our argument is based on the comparison with the lowest singlet excitation energy and herein reported, carefully determined absorption spectrum of ferrocene. When deriving vertical energies from the experimental band maxima, we go beyond the routine vertical energy approximation by introducing vibronic corrections based on simulated vibrational envelopes. The benchmarking result confirms the high accuracy of the CCSD(T) method (in particular, for UCCSD(T) based on Hartree-Fock orbitals we find for our dataset: maximum error 0.12 eV, weighted mean absolute error 0.07 eV, weighted mean signed error 0.01 eV). The high accuracy of the single-reference method is corroborated by the analysis of a multiconfigurational character of the complete active space wave function for the triplet state of ferrocene. On the DFT side, our results confirm the non-universality problem with approximate functionals. The present study is an important step toward establishing an extensive and representative benchmark set of experiment-derived spin-state energetics for transition metal complexes.

Accurate Computation of the Absorption Spectrum of Chlorophyll a with Pair Natural Orbital Coupled Cluster Methods

The ability to accurately compute low-energy excited states of chlorophylls is critically important for understanding the vital roles they play in light harvesting, energy transfer, and photosynthetic charge separation. The challenge for quantum chemical methods arises both from the intrinsic complexity of the electronic structure problem and, in the case of biological models, from the need to account for protein–pigment interactions. In this work, we report electronic structure calculations of unprecedented accuracy for the low-energy excited states in the Q and B bands of chlorophyll a. This is achieved by using the newly developed domain-based local pair natural orbital (DLPNO) implementation of the similarity transformed equation of motion coupled cluster theory with single and double excitations (STEOM-CCSD) in combination with sufficiently large and flexible basis sets. The results of our DLPNO–STEOM-CCSD calculations are compared with more approximate approaches. The results demonstrate that, in contrast to time-dependent density functional theory, the DLPNO–STEOM-CCSD method provides a balanced performance for both absorption bands. In addition to vertical excitation energies, we have calculated the vibronic spectrum for the Q and B bands through a combination of DLPNO–STEOM-CCSD and ground-state density functional theory frequency calculations. These results serve as a basis for comparison with gas-phase experiments.

Intersystem Crossing and Triplet Dynamics in an Iron(II) N-Heterocyclic Carbene Photosensitizer

The electronic excited states of the iron(II) complex [Fe^II(tpy)(pyz-NHC)]²⁺ [tpy = 2,2′:6′,2″-terpyridine; pyz-NHC = 1,1′-bis(2,6-diisopropylphenyl)pyrazinyldiimidazolium-2,2′-diylidene] and their relaxation pathways have been theoretically investigated. To this purpose, trajectory surface-hopping simulations within a linear vibronic coupling model including a 244-dimensional potential energy surface (PES) with 20 singlet and 20 triplet coupled states have been used. The simulations show that, after excitation to the lowest-energy absorption band of predominant metal-to-ligand charge-transfer character involving the tpy ligand, almost 80% of the population undergoes intersystem crossing to the triplet manifold in about 50 fs, while the remaining 20% decays through internal conversion to the electronic ground state in about 300 fs. The population transferred to the triplet states is found to deactivate into two different regions of the PESs, one where the static dipole moment is small and shows increased metal-centered character and another with a large static dipole moment, where the electron density is transferred from the tpy to pyz-NHC ligand. Coherent oscillations of 400 fs are observed between these two sets of triplet populations, until the mixture equilibrates to a ratio of 60:40. Finally, the importance of selecting suitable normal modes is highlighted—a choice that can be far from straightforward in transition-metal complexes with hundreds of degrees of freedom.

Trajectory surface-hopping simulations with a linear vibronic coupling model reveal the competition of major intersystem crossing versus minor internal conversion dynamics in an iron(II) N-heterocyclic carbene dye. The triplet population bifurcates into two regions of the potential energy surfaces, characterized by small and large static dipole moments due to different electronic character and showing coherent oscillations of 400 fs until both triplet populations coexist in a mixture of 60:40.

Machine Learning for Absorption Cross Sections

We present a machine learning (ML) method to accelerate the nuclear ensemble approach (NEA) for computing absorption cross sections. ML-NEA is used to calculate cross sections on vast ensembles of nuclear geometries to reduce the error due to insufficient statistical sampling. The electronic properties-excitation energies and oscillator strengths-are calculated with a reference electronic structure method only for a relatively few points in the ensemble. The KREG model (kernel-ridge-regression-based ML combined with the RE descriptor) as implemented in MLatom is used to predict these properties for the remaining tens of thousands of points in the ensemble without incurring much of additional computational cost. We demonstrate for two examples, benzene and a 9-dicyanomethylene derivative of acridine, that ML-NEA can produce statistically converged cross sections even for very challenging cases and even with as few as several hundreds of training points.

Simulation of Excitation by Sunlight in Mixed Quantum-Classical Dynamics

This paper proposes a method to simulate nonadiabatic dynamics initiated by thermal light, including solar radiation, in the frame of mixed quantum-classical (MQC) methods, like surface hopping. The method is based on the Chenu-Brumer approach, which treats the thermal radiation as an ensemble of coherent pulses. It is composed of three steps: (1) sampling initial conditions from a broad blackbody spectrum, (2) dynamics propagation using conventional methods, and (3) ensemble averaging considering the field and realization time of the pulses. The application of MQC dynamics with pulse ensembles (MQC-PE) to a model system of nucleic acid photophysics showed the emergence of a steady excited-state population. In another test case, modeling retinal photophysics, MQC-PE predicted that although the photoisomerization occurs within 200 fs, it may take tens of microseconds of continuous solar irradiation to photoactivate a single retinal. Such emergent long timescales may impact our understanding of biological and technological phenomena occurring under solar radiation.

Simulation of Excitation by Sunlight in Mixed Quantum-Classical Dynamics

This paper proposes a method to simulate nonadiabatic dynamics initiated by thermal light, including solar radiation, in the frame of mixed quantum-classical (MQC) methods, like surface hopping. The method is based on the Chenu–Brumer approach, which treats the thermal radiation as an ensemble of coherent pulses. It is composed of three steps: (1) sampling initial conditions from a broad blackbody spectrum, (2) dynamics propagation using conventional methods, and (3) ensemble averaging considering the field and realization time of the pulses. The application of MQC dynamics with pulse ensembles (MQC-PE) to a model system of nucleic acid photophysics showed the emergence of a steady excited-state population. In another test case, modeling retinal photophysics, MQC-PE predicted that although the photoisomerization occurs within 200 fs, it may take tens of microseconds of continuous solar irradiation to photoactivate a single retinal. Such emergent long timescales may impact our understanding of biological and technological phenomena occurring under solar radiation.

The effect of N-heterocyclic carbene units on the absorption spectra of Fe(ii) complexes: a challenge for theory

The absorption spectra of five Fe(ii) homoleptic and heteroleptic complexes containing strong sigma-donating N-heterocyclic carbene (NHC) and polypyridyl ligands have been theoretically characterized using a tuned range-separation functional.

Synthesis, photophysical and electrochemical properties of novel and highly fluorescent difluoroboron flavanone β-diketonate complexes

A convenient method for the synthesis of flavanone difluoroboron complexes with strong absorption and fluorescence properties is described.