Effective Inclusion of Mechanical and Electrical Anharmonicity in Excited Electronic States: VPT2-TDDFT Route

Assessing the accuracy of TD-DFT excited-state geometries through optimal tuning with GW energy levels

We study the accuracy of excited state (ES) geometries using optimally tuned LC-PBE functionals with tuning based on GW quasiparticle energies. We compare the results obtained with the PBE, PBE0, non-tuned, and tuned LC-PBE functionals with available high-level CC reference values as well as experimental data. First, we compare ES geometrical parameters obtained for three different types of systems: molecules composed of a few atoms, 4-(dimethylamino)benzonitrile (DMABN), and conjugated dyes. To this end, we used wave-function results as benchmarks. Next, we evaluate the accuracy of the theoretically simulated spectra as compared to the experimental ones for five large dyes. Our results show that, besides small compact molecules for which tuning LC-PBE does not allow obtaining geometries more accurate than those computed with standard functionals, tuned range-separated functionals are clearly to be favored, not only for ES geometries but also for 0-0 energies, band shapes, and intensities for absorption and emission spectra. In particular, the results indicate that GW-tuned LC-PBE functionals provide improved matching with experimental spectra as compared to conventionally tuned functionals. It is an open question whether TD-DFT with GW-tuned functionals can qualitatively mimic the actual many-body Bethe-Salpeter (BSE/GW) formalism for which analytic ionic gradients remain to be developed.

Unbiased Comparison between Theoretical and Experimental Molecular Structures and Properties: Toward an Accurate Reduced-Cost Evaluation of Vibrational Contributions

The tremendous development of hardware and software is constantly increasing the role of quantum chemical (QC) computations in the assignment and interpretation of experimental results. However, an unbiased comparison between theory and experiment requires the proper account of vibrational averaging effects. In particular, high-resolution spectra in the gas phase are now available for molecules containing up to about 50 atoms, which are too large for a brute-force approach with the available QC methods of sufficient accuracy. In the present paper, we introduce hybrid approaches, which allow the accurate evaluation of vibrational averaging effects for molecules of this size beyond the harmonic approximation, with special attention being devoted to rotational constants. After the validation of new tools for relatively small molecules, the β-estradiol hormone and a prototypical molecular motor have been considered to witness the feasibility of accurate computations for large molecules.

Reassignment of the Photoelectron Spectrum of Methylketene Using a Hybrid Model of Harmonic and Anharmonic Oscillators to Compute Franck-Condon Factors

We constructed a hybrid model of harmonic and anharmonic oscillators to compute Franck-Condon factors and interpret the photoelectron spectrum of methylketene. The equilibrium structures of methylketene and its cation were optimized, and then, the harmonic and anharmonic vibrational frequencies were computed using the B3LYP, PBE0, APFD, and ωB97XD approaches of the density functional theory. The photoelectron spectrum of methylketene was simulated by computing the Franck-Condon factors with both the harmonic and hybrid models. The adiabatic ionization energy of methylketene was computed by using the CCSD(T) approach extrapolating to the complete basis set limit. The simulated photoelectron spectra are consistent with those from the experiment for both the harmonic and hybrid models. However, the error in band positions is reduced by using the hybrid model. The computed adiabatic ionization energies of methylketene are in agreement with the experiment, with the smallest error being 0.017 eV. Our interpretation based on the theoretical spectrum led to the reassignment of the experimental photoelectron spectrum of methylketene.

Rational Design of Asymmetric Polymethines to Attain NIR(II) Bioimaging at >1100 nm

Organic molecules having emission in the NIR(II) region are emergent and receiving enormous attention. Unfortunately, attaining accountable organic emission intensity around the NIR(II) region is hampered by the dominant internal conversion operated by the energy gap law, where the emission energy gap and the associated internal reorganization energy λ_int play key roles. Up to the current stage, the majority of the reported organic NIR(II) emitters belong to those polymethines terminated by two symmetric chromophores. Such a design has proved to have a small λ_int that greatly suppresses the internal conversion. However, the imposition of symmetric chromophores is stringent, limiting further development of organic NIR(II) dyes in diversity and versatility. Here, we propose a new concept where as far as the emissive state of the any asymmetric polymethines contains more or less equally transition density between two terminated chromophores, λ_int can be as small as that of the symmetric polymethines. To prove the concept, we synthesize a series of new polymethines terminated by xanthen-9-yl-benzoic acid and 2,4-diphenylthiopyrylium derivatives, yielding AJBF1112 and AEBF1119 that reveal emission peak wavelength at 1112 and 1119 nm, respectively. The quantum yield is higher than all synthesized symmetric polymethines of 2,4-diphenylthiopyrylium derivatives (SC1162, 1182, 1185, and 1230) in this study. λ_int were calculated to be as small as 6.2 and 7.3 kcal/mol for AJBF1112 and AEBF1119, respectively, proving the concept. AEBF1119 was further prepared as a polymer dot to demonstrate its in vitro specific cellular imaging and in vivo tumor/bone targeting in the NIR(II) region.

Applicability of the Thawed Gaussian Wavepacket Dynamics to the Calculation of Vibronic Spectra of Molecules with Double-Well Potential Energy Surfaces

Simulating vibrationally resolved electronic spectra of anharmonic systems, especially those involving double-well potential energy surfaces, often requires expensive quantum dynamics methods. Here, we explore the applicability and limitations of the recently proposed single-Hessian thawed Gaussian approximation for the simulation of spectra of systems with double-well potentials, including 1,2,4,5-tetrafluorobenzene, ammonia, phosphine, and arsine. This semiclassical wavepacket approach is shown to be more robust and to provide more accurate spectra than the conventional harmonic approximation. Specifically, we identify two cases in which the Gaussian wavepacket method is especially useful due to the breakdown of the harmonic approximation: (i) when the nuclear wavepacket is initially at the top of the potential barrier but delocalized over both wells, e.g., along a low-frequency mode, and (ii) when the wavepacket has enough energy to classically go over the low potential energy barrier connecting the two wells. The method is efficient and requires only a single classical ab initio molecular dynamics trajectory, in addition to the data required to compute the harmonic spectra. We also present an improved algorithm for computing the wavepacket autocorrelation function, which guarantees that the evaluated correlation function is continuous for arbitrary size of the time step.

Sum-frequency vibrational spectroscopy of methanol at interfaces due to Fermi resonance

We present a theoretical method of studying sum-frequency vibrational spectroscopy for the CH3 group of methanol at interfaces due to Fermi resonance, which provides a novel and untraditional point of view with respect to traditional approaches.

Spectroscopic Manifestations of Indirect Vibrational State Mixing: Novel Anharmonic Effects on a Prereactive H Atom Transfer Surface

The NH stretch region of the IR spectrum of methyl anthranilate is modeled in the S₁ state to understand the connection between the absence of this fundamental in the fluorescence-dip infrared spectra of Blodgett et al. [Phys. Chem. Chem. Phys. 2020, 22, 14077] and its relevance to the H atom dislocation that occurs upon electronic excitation. A set of coordinates are chosen that highlight the role of certain low-frequency modes. A Hamiltonian is developed in which a large-amplitude two-dimensional surface describing the H-bonded H atom is linearly and quadratically coupled to the remaining degrees of freedom which are treated at the harmonic level. The surface is calculated within the time-dependent density functional theory framework by using the B3LYP/6-311++(d, p) level of theory with dispersion. Our spectral results show that indirect couplings lead to massive intensity sharing over hundreds of wavenumbers. This sharing is predicted to be dramatically reduced upon deuteration. The spectral broadening mechanism is found to involve off-resonant doorway states that are themselves strongly coupled to states nearly degenerate with the NH stretch fundamental and represents a complementary mechanism to previous explanations based on Fermi resonance or the presence of Franck-Condon like combination bands with low-frequency motions. Consistent with the spectra predictions, time-dependent calculations show that if the NH stretch fundamental were excited with an ultrafast laser, it would decay within 40 fs. The competition between H atom dislocation and vibrational relaxation is discussed.

Interplay of stereo-electronic, vibronic and environmental effects in tuning the chiroptical properties of an Ir(III) cyclometalated N-heterocyclic carbene

Chiroptical spectra are among the most suitable techniques for investigating the ground and excited electronic states of chiral systems, but their interpretation is not straightforward and strongly benefits from quantum chemical simulations, provided that the employed computational model is sufficiently accurate and deals properly with stereo-electronic, vibrational averaging and environmental effects. Since the synergy among all these effects is only rarely accounted for, especially for large and flexible organometallic systems, the main aim of this contribution is to illustrate the latest developments of computational approaches rooted into the density functional theory for describing stereo-electronic effects and complemented by effective techniques to deal with vibrational modulation effects and solvatochromic shifts. In this connection, chiral iridium complexes offer an especially suitable case study in view of their bright phosphorescence, which is particularly significant for building effective light emitting diodes (OLEDs) and biomarkers and can be finely tuned by the nature of the metal ligands. For instance, a recently synthesized family of cycloiridiated complexes, KC and KD, bearing a pentahelicenic N-heterocyclic carbene (KB), has shown an enhanced long-lasting, bright phosphorescence. Deeper insights into the still unclear nature and origin of the enhancement could be gained by the interpretation of the chiroptical spectra, which is quite challenging in view of the presence of two sources of chirality, the chiral center on Ir and the chiral axis related to the helicene ligand, in addition to the relativistic effects related to the presence of the Ir center. At the same time, the large dimensions of KC and KD hamper the use of the most sophisticated (but prohibitively expensive) computational models, so that more approximate approaches must be validated on a suitable model compound. To this end, after optimizing the computational scheme on a model system devoid of the helicene moiety (KA), we have performed a comprehensive investigation of the KC and KD spectra, whose interpretation is further aided by novel graphical tools. The discussion and analysis of the results will not be focused on the theoretical background, but, rather, on practical details (specific functional, basis set, vibronic model, solvent regime) with the aim of providing general guidelines for the use of last-generation computational spectroscopy tools also by non-specialists.

The missing NH stretch fundamental in S1 methyl anthranilate: IR-UV double resonance experiments and local mode theory

The infrared spectra of jet-cooled methyl anthranilate (MA) and the MA-H₂O complex are reported in both S₀ and S₁ states, recorded using fluorescence-dip infrared (FDIR) spectroscopy under jet-cooled conditions. Using a combination of local mode CH stretch modeling and scaled harmonic vibrational character, a near-complete assignment of the infrared spectra is possible over the 1400-3700 cm^-1 region. While the NH stretch fundamentals are easily observed in the S₀ spectrum, in the S₁ state, the hydrogen bonded NH stretch shift is not readily apparent. Scaled harmonic calculations predict this fundamental at just below 2900 cm^-1 with an intensity around 400 km mol^-1. However, the experimental spectrum shows no evidence of this transition. A local mode theory is developed in which the NH stretch vibration is treated adiabatically. Minimizing the energy of the corresponding stretch state with one quantum of excitation leads to a dislocation of the H atom where there is equal sharing between N and O atoms. The sharing occurs as a result of significant molecular arrangement due to strong coupling of this NH stretch to other internal degrees of freedom and in particular to the contiguous HNC bend. A two-dimensional model of the coupling between the NH stretch and this bend highlights important nonlinear effects that are not captured by low order vibrational perturbation theory. In particular, the model predicts a dramatic dilution of the NH stretch oscillator strength over many transitions spread over more than 1000 cm^-1, making it difficult to observe experimentally.

Vibrations of the guanine-cytosine pair in chloroform: an anharmonic computational study

We compute at the anharmonic level the vibrational spectra of the Watson-Crick dimer formed by guanosine (G) and cytidine (C) in chloroform, together with those of G, C and the most populated GG dimer. The spectra for deuterated and partially deuterated GC are also computed. We use DFT calculations, with B3LYP and CAM-B3LYP as reference functionals. Solvent effects from chloroform are included via the Polarizable Continuum Model (PCM), and by performing tests on models including up two chloroform molecules. Both B3LYP and CAM-B3LYP calculations reproduce the shape of the experimental spectra well in the fingerprint region (1500-1700 cm-1) and in the N-H stretching region (2800-3600 cm-1), with B3LYP providing better quantitative agreement with experiments. According to our calculations, the N-H amido streching mode of G falls at ∼2900 cm-1, while the N-H amino of G and C falls at ∼3100 cm-1 when hydrogen-bonded, or ∼3500 cm-1 when free. Overtone and combination bands strongly contribute to the absorption band at ∼3300 cm-1. Inclusion of bulk solvent effects significantly increases the accuracy of the computed spectra, while solute-solvent interactions have a smaller, though still noticeable, effect. Some key aspects of the anharmonic treatment of strongly vibrationally coupled supermolecular systems and the related methodological issues are also discussed.

Adiabatic-Molecular Dynamics Generalized Vertical Hessian Approach: A Mixed Quantum Classical Method To Compute Electronic Spectra of Flexible Molecules in the Condensed Phase

We present a general mixed quantum classical method that couples classical molecular dynamics (MD) and vibronic models to compute the shape of electronic spectra of flexible molecules in the condensed phase without, in principle, any phenomenological broadening. It is based on a partition of the nuclear motions of the solute + solvent system in "soft" and "stiff" vibrational modes and an adiabatic hypothesis that assumes that stiff modes are much faster than soft ones. In this framework, the spectrum is rigorously expressed as a conformational integral of quantum vibronic spectra along the stiff coordinates only. Soft modes enter at the classical level through the conformational distribution that is sampled with classical MD runs. In each configuration, reduced-dimensionality quadratic Hamiltonians are built in the space of the stiff coordinates only, thanks to a generalization of the Vertical Hessian harmonic model and an iterative application of projectors in internal coordinates to remove soft modes. Quantum vibronic spectra, specific for each sampled configuration of the soft coordinates, are then computed at the desired temperature with efficient time-dependent techniques, and the global spectrum simply arises from their average. For consistency of the whole procedure, classical MD runs are performed with quantum-mechanically derived force fields, parameterized at the same level of theory selected for generating the quadratic Hamiltonians along the stiff coordinates. Application to N-methyl-6-oxyquinolinium betaine in water, dithiophene in ethanol, and cyanidine in water is presented to show the performance of the method.

Anharmonic excited state frequencies of para-difluorobenzene, toluene and catechol using analytic RI-CC2 second derivatives

Analytic second nuclear derivatives for excited electronic state energies have been implemented for the resolution-of-the-identity accelerated CC2, CIS(D_∞) and ADC(2) models. Our efficient implementation with O(N²) memory demands enables the treatment of medium sized molecules with large basis sets and high numerical precision and thereby paves the way for semi-numerical evaluation of the higher-order derivatives required for anharmonic corrections to excited state vibrational frequencies. We compare CC2 harmonic and anharmonic excited state frequencies with experimental values for para-difluorobenzene, toluene and catechol. Basis set problems occur for out-of-plane bending vibrations due to intramolecular basis set superposition error. For non-planar molecules and in plane modes of planar molecules, the agreement between theory and experiment is better than 30 cm^-1 on average and we reassign a number of experimental bands on the basis of the ab initio predictions.

Single-Hessian thawed Gaussian approximation

To alleviate the computational cost associated with on-the-fly ab initio semiclassical calculations of molecular spectra, we propose the single-Hessian thawed Gaussian approximation in which the Hessian of the potential energy at all points along an anharmonic classical trajectory is approximated by a constant matrix. The spectra obtained with this approximation are compared with the exact quantum spectra of a one-dimensional Morse potential and with the experimental spectra of ammonia and quinquethiophene. In all cases, the single-Hessian version performs almost as well as the much more expensive on-the-fly ab initio thawed Gaussian approximation and significantly better than the global harmonic schemes. Remarkably, unlike the thawed Gaussian approximation, the proposed method conserves energy exactly, despite the time dependence of the corresponding effective Hamiltonian, and, in addition, can be mapped to a higher-dimensional time-independent classical Hamiltonian system. We also provide a detailed comparison with several related approximations used for accelerating prefactor calculations in semiclassical simulations.

On the simulation of vibrationally resolved electronic spectra of medium-size molecules: the case of styryl substituted BODIPYs

BODIPY dyes are used in a variety of applications because of their peculiar spectroscopic and photo-physical properties that vary depending on the stereochemistry of the functional groups attached to the boron-dipyrromethene core structure. In this work, we have applied several computational methods, adapted for semi-rigid molecules based on the Franck-Condon principle, for the study of the optical properties of BODIPY systems and for the understanding of the influence of functional groups on their spectroscopic features. We have analyzed the electronic spectra of two styryl substituted BODIPY molecules of technological interest, properly taking into account the vibronic contribution. For comparison with recently recorded experimental data in methanol, the vibrationally resolved electronic spectra of these systems were computed using both Time-Independent (TI) and Time-Dependent (TD) formalisms. The first step toward the analysis of optical properties of the styryl modified BODIPYs was a benchmark of several density functionals, to select the most appropriate one. We have found that all benchmarked functionals provide good results in terms of band shape but some of them show strong discrepancies in terms of band position. Beyond the issue of the electronic structure calculation method, different levels of sophistication can be adopted for the calculation of vibronic transitions. In this study, the effect of mode couplings and the influence of the Herzberg-Teller terms on the theoretical spectra has been investigated. It has been found that all levels of theory considered give reproducible results for the investigated systems: band positions and shapes are similar at all levels and little improvements have been found in terms of band shape with the inclusion of Herzberg-Teller effect. Inclusion of temperature effects proved to be challenging due to the important impact of large amplitude motions. Better agreement can be achieved by adopting a suitable set of coordinates coupled with a reduced-dimensionality scheme.

Time-Dependent Formulation of Resonance Raman Optical Activity Spectroscopy

In this work, we extend the theoretical framework recently developed for the simulation of resonance Raman (RR) spectra of medium-to-large sized systems to its chiral counterpart, namely, resonance Raman optical activity (RROA). The theory is based on a time-dependent (TD) formulation, with the transition tensors obtained as half-Fourier transforms of the appropriate cross-correlation functions. The implementation has been kept as general as possible, supporting adiabatic and vertical models for the PES representation, both in Cartesian and internal coordinates, with the possible inclusion of Herzberg-Teller (HT) effects. Thanks to the integration of this TD-RROA procedure within a general-purpose quantum-chemistry program, both solvation and leading anharmonicity effects can be included in an effective way. The implementation is validated on one of the smallest chiral molecule (methyloxirane). Practical applications are illustrated with three medium-size organic molecules (naproxen-OCD₃, quinidine and 2-Br-hexahelicene), whose simulated spectra are compared to the corresponding experimental data.

High-Resolution Absorption and Electronic Circular Dichroism Spectra of (R)-(+)-1-Phenylethanol. Confident Interpretation Based on the Synergy between Experiments and Computations

Using density functional theory and its time-dependent extension for excited states, the S₀ →S₁ high-resolution vibronic absorption and electronic circular dichroism spectra of (R)-(+)-1-phenylethanol are computed and compared to experimental spectra measured in jet-cooled conditions in the region within 1000 cm^-1 of the 0-0 transition. The agreement between theory and computation is satisfactory and allows a confident assignment of several experimental bands in terms of fundamentals of different modes. Cases are documented for which the analysis of optical anisotropy factors, owing to their signed nature, remarkably enhances the possibility of a robust assignment of the experimental absorption bands. Computational analysis shows that the experimental spectra are dominated by Herzberg-Teller contributions and that the electronic circular dichroism spectrum and the anisotropy factors are also strongly modulated by the effect of Duschinsky mixings.

Theoretical studies of atmospheric molecular complexes interacting with NIR to UV light

Theoretical structural and spectroscopic data for weakly bonded atmospheric complexes of formaldehyde interacting with Earth’s electromagnetic spectrum.

Accurate Excited-State Geometries: A CASPT2 and Coupled-Cluster Reference Database for Small Molecules

We present an investigation of the excited-state structural parameters determined for a large set of small compounds with the dual goals of defining reference values for further works and assessing the quality of the geometries obtained with relatively cheap computational approaches. In the first stage, we compare the excited-state geometries obtained with ADC(2), CC2, CCSD, CCSDR(3), CC3, and CASPT2 and large atomic basis sets. It is found that CASPT2 and CC3 results are generally in very good agreement with one another (typical differences of ca. 3 × 10^–3 Å) when all electrons are correlated and when the aug-cc-pVTZ atomic basis set is employed with both methods. In a second stage, a statistical analysis reveals that, on the one hand, the excited-state (ES) bond lengths are much more sensitive to the selected level of theory than their ground-state (GS) counterparts and, on the other hand, that CCSDR(3) is probably the most cost-effective method delivering accurate structures. Indeed, CCSD tends to provide too compact multiple bond lengths on an almost systematic basis, whereas both CC2 and ADC(2) tend to exaggerate these bond distances, with more erratic error patterns, especially for the latter method. The deviations are particularly marked for the polarized CO and CN bonds, as well as for the puckering angle in formaldehyde homologues. In the last part of this contribution, we provide a series of CCSDR(3) GS and ES geometries of medium-sized molecules to be used as references in further investigations.