Assessment of Franck-Condon Methods for Computing Vibrationally Broadened UV-vis Absorption Spectra of Flavin Derivatives: Riboflavin, Roseoflavin, and 5-Thioflavin

Electronic Structure Methods for Simulating Flavin's Spectroscopy and Photophysics: Comparison of Multi-reference, TD-DFT, and Single-Reference Wave Function Methods

The use of flavins and flavoproteins in photocatalytic, sensing, and biotechnological applications has led to a growing interest in computationally modeling the excited-state electronic structure and photophysics of flavin. However, there is limited consensus regarding which computational methods are appropriate for modeling flavin's photophysics. We compare the energies of low-lying excited states of flavin computed with time-dependent density functional theory (TD-DFT), equation-of-motion coupled cluster (EOM-EE-CCSD), scaled opposite-spin configuration interaction [SOS-CIS(D)], multiconfiguration pair-density functional theory (MC-PDFT), and several multireference perturbation theory (MR-PT2) methods. In the first part, we focus on excitation energies of the first singlet excited state (S1) of five different redox and protonation states of flavin, with the goal of finding a suitable active space for MR-PT2 calculations. In the second part, we construct two sets of one-dimensional potential energy surfaces connecting the S0 and S1 equilibrium geometries (S0-S1 path) and the S1 (π,π*) and S2 (n,π*) equilibrium geometries (S1-S2 path). The first path therefore follows a Franck-Condon active mode of flavin while the second path maps crossings points between low-lying singlet and triplet states in flavin. We discuss the similarities and differences in the TD-DFT, EOM-EE-CCSD, SOS-CIS(D), MC-PDFT and MR-PT2 energy profiles along these paths. We find that (TD-)DFT methods are suitable for applications such as simulating the spectra of flavins but are inconsistent with several other methods when used for some geometry optimizations and when describing the energetics of dark (n,π*) states. MR-PT2 methods show promise for the simulation of flavin's low-lying excited states, but the selection of orbitals for the active space and the number of roots used for state averaging must be done carefully to avoid artifacts. Some properties, such as the intersystem crossing geometry and energy between the S1 (π,π*) and T2 (n,π*) states, may require additional benchmarking before they can be determined quantitatively.

Red Shift in the Absorption Spectrum of Phototropin LOV1 upon the Formation of a Semiquinone Radical: Reconstructing the Orbital Architecture

Flavin mononucleotide (FMN) is a ubiquitous blue-light pigment due to its ability to drive one- and two-electron transfer reactions. In both light-oxygen-voltage (LOV) domains of phototropin from the green algae Chlamydomonas reinhardtii, FMN is noncovalently bound. In the LOV1 cysteine-to-serine mutant (C57S), light-induced electron transfer from a nearby tryptophan occurs, and a transient spin-correlated radical pair (SCRP) is formed. Within this photocycle, nuclear hyperpolarization is created by the solid-state photochemically induced dynamic nuclear polarization (photo-CIDNP) effect. In a side reaction, a stable protonated semiquinone radical (FMNH·) forms undergoing a significant bathochromic shift of the first electronic transition from 445 to 591 nm. The incorporation of phototropin LOV1-C57S into an amorphous trehalose matrix, stabilizing the radical, allows for application of various magnetic resonance experiments at ambient temperatures, which are combined with quantum-chemical calculations. As a result, the bathochromic shift of the first absorption band is explained by lifting the degeneracy of the molecular orbital energy levels for electrons with alpha and beta spins in FMNH· due to the additional electron.

How Aqueous Solvation Impacts the Frequencies and Intensities of Infrared Absorption Bands in Flavin: The Quest for a Suitable Solvent Model

FTIR spectroscopy accompanied by quantum chemical simulations can reveal important information about molecular structure and intermolecular interactions in the condensed phase. Simulations typically account for the solvent either through cluster quantum mechanical (QM) models, polarizable continuum models (PCM), or hybrid quantum mechanical/molecular mechanical (QM/MM) models. Recently, we studied the effect of aqueous solvent interactions on the vibrational frequencies of lumiflavin, a minimal flavin model, using cluster QM and PCM models. Those models successfully reproduced the relative frequencies of four prominent stretching modes of flavin's isoalloxazine ring in the diagnostic 1450-1750 cm-1 range but poorly reproduced the relative band intensities. Here, we extend our studies on this system and account for solvation through a series of increasingly sophisticated models. Only by combining elements of QM clusters, QM/MM, and PCM approaches do we obtain an improved agreement with the experiment. The study sheds light more generally on factors that can impact the computed frequencies and intensities of IR bands in solution.

Tuning the Quantum Chemical Properties of Flavins via Modification at C8

Flavins are central to countless enzymes but display different reactivities depending on their environments. This is understood to reflect modulation of the flavin electronic structure. To understand changes in orbital natures, energies, and correlation over the ring system, we begin by comparing seven flavin variants differing at C8, exploiting their different electronic spectra to validate quantum chemical calculations. Ground state calculations replicate a Hammett trend and reveal the significance of the flavin π-system. Comparison of higher-level theories establishes CC2 and ACD(2) as methods of choice for characterization of electronic transitions. Charge transfer character and electron correlation prove responsive to the identity of the substituent at C8. Indeed, bond length alternation analysis demonstrates extensive conjugation and delocalization from the C8 position throughout the ring system. Moreover, we succeed in replicating a particularly challenging UV/Vis spectrum by implementing hybrid QM/MM in explicit solvents. Our calculations reveal that the presence of nonbonding lone pairs correlates with the change in the UV/Vis spectrum observed when the 8-methyl is replaced by NH₂, OH, or SH. Thus, our computations offer routes to understanding the spectra of flavins with different modifications. This is a first step toward understanding how the same is accomplished by different binding environments.

Ionic Atmosphere Effect on the Absorption Spectrum of a Flavoprotein: A Reminder to Consider Solution Ions

This study utilizes the FMN-dependent NADH:quinone oxidoreductase from Pseudomonas aeruginosa PAO1 to investigate the effect of introducing an active site negative charge on the flavin absorption spectrum both in the absence and presence of a long-range electrostatic potential coming from solution ions. There were no observed changes in the flavin UV-visible spectrum when an active site tyrosine (Y277) becomes deprotonated in vitro. These results could only be reproduced computationally using average solvent electrostatic configuration (ASEC) QM/MM simulations that include both positive and negative solution ions. The same calculations performed with minimal ions to neutralize the total protein charge predicted that deprotonating Y277 would significantly alter the flavin absorption spectrum. Analyzing the distribution of solution ions indicated that the ions reorganize around the protein surface upon Y277 deprotonation to cancel the effect of the tyrosinate on the flavin absorption spectrum. Additional biochemical experiments were performed to test this hypothesis.

Investigating the biological actions of some Schiff bases using density functional theory study

Schiff base ligands have wide varieties of application in several fields. One of which is the biological actions they possess such as anti-fungal, anti-bacterial, anti-malarial, and anti-viral characteristics. In this study, some synthesized phenylimino-based Schiff bases were investigated using density functional theory (DFT) to unravel their biological descriptors. The gas-phase quantum chemical calculation was done on the Schiff base 3-((E)-(phenylimino)methyl)benzene-1,2-diol and other synthesized analogues to evaluate their reactivity and stability properties including the substituent effect on the basic molecule. The Coulomb-attenuating method (CAM-B3LYP) functional was employed for the theoretical calculations. The Nuclear Magnetic Resonance (NMR), Fourier Transform-Infrared (FT-IR), Ultraviolet/visible spectroscopies calculated agrees with the experimental values. The obtained charge transfer and electronic features provide useful information regarding the active sites for biological application in the compounds.

Spectroscopic evidence for direct flavin-flavin contact in a bifurcating electron transfer flavoprotein

A remarkable charge transfer (CT) band is described in the bifurcating electron transfer flavoprotein (Bf-ETF) from Rhodopseudomonas palustris (RpaETF). RpaETF contains two FADs that play contrasting roles in electron bifurcation. The Bf-FAD accepts electrons pairwise from NADH, directs one to a lower-reduction midpoint potential (E°) carrier, and the other to the higher-E° electron transfer FAD (ET-FAD). Previous work noted that a CT band at 726 nm formed when ET-FAD was reduced and Bf-FAD was oxidized, suggesting that both flavins participate. However, existing crystal structures place them too far apart to interact directly. We present biochemical experiments addressing this conundrum and elucidating the nature of this CT species. We observed that RpaETF missing either FAD lacked the 726 nm band. Site-directed mutagenesis near either FAD produced altered yields of the CT species, supporting involvement of both flavins. The residue substitutions did not alter the absorption maximum of the signal, ruling out contributions from residue orbitals. Instead, we propose that the residue identities modulate the population of a protein conformation that brings the ET-flavin and Bf-flavin into direct contact, explaining the 726 nm band based on a CT complex of reduced ET-FAD and oxidized Bf-FAD. This is corroborated by persistence of the 726 nm species during gentle protein denaturation and simple density functional theory calculations of flavin dimers. Although such a CT complex has been demonstrated for free flavins, this is the first observation of such, to our knowledge, in an enzyme. Thus, Bf-ETFs may optimize electron transfer efficiency by enabling direct flavin-flavin contact.

Femtosecond-to-nanosecond dynamics of flavin mononucleotide monitored by stimulated Raman spectroscopy and simulations

Flavin mononucleotide (FMN) belongs to the large family of flavins, ubiquitous yellow-coloured biological chromophores that contain an isoalloxazine ring system. As a cofactor in flavoproteins, it is found in various enzymes and photosensory receptors, like those featuring the light-oxygen-voltage (LOV) domain. The photocycle of FMN is triggered by blue light and proceeds via a cascade of intermediate states. In this work, we have studied isolated FMN in an aqueous solution in order to elucidate the intrinsic electronic and vibrational changes of the chromophore upon excitation. The ultrafast transitions of excited FMN were monitored through the joint use of femtosecond stimulated Raman spectroscopy (FSRS) and transient absorption spectroscopy encompassing a time window between 0 ps and 6 ns with 50 fs time resolution. Global analysis of the obtained transient visible absorption and transient Raman spectra in combination with extensive quantum chemistry calculations identified unambiguously the singlet and triplet FMN populations and addressed solvent dynamics effects. The good agreement between the experimental and theoretical spectra facilitated the assignment of electronic transitions and vibrations. Our results represent the first steps towards more complex experiments aimed at tracking structural changes of FMN embedded in light-inducible proteins upon photoexcitation.

Nitrofurantoin-melamine monohydrate (cocrystal hydrate): Probing the role of H-bonding on the structure and properties using quantum chemical calculations and vibrational spectroscopy

Cocrystal monohydrate of nitrofurantoin (NF) with melamine (MELA) has been studied as NF is an antibacterial drug used for the treatment of urinary tract infections. The structure of nitrofurantoin-melamine-monohydrate (NF-MELA-H₂O) is characterized by FT-IR and FT-Raman spectroscopy. The energies and vibrational frequencies of the optimized structures calculated using quantum chemical calculations. Supported by normal coordinate analyses and potential energy distributions (PEDs), the complete vibrational assignments recommended for the observed fundamentals of cocrystal hydrate. With the aim of inclusion of all the H-bond interactions, dimer of NF-MELA-H₂O has been studied as only two molecules of cocrystal hydrate are present in the unit cell. By the study of dimeric model consistent assignment of the FT-IR and FT-Raman spectrum obtained. H-bonds are of essential importance in an extensive range of molecular sciences. The vibrational analyses depict existence of H-bonding (O-H⋯N) between water O-H and pyridyl N atom of MELA in both monomer and dimer. To probe the strength and nature of H-bonding in monomer and dimer, topological parameters such as electron density (ρ_BCP), Laplacian of electron density (∇²ρ_BCP), total electron energy density (H_BCP) and H-bond energy (E_HB) at bond critical points (BCP) are evaluated by quantum theory of atoms in molecules (QTAIM). Natural bond orbitals (NBOs) analyses are carried out to study especially the intra and intermolecular H-bonding and their second order stabilization energy (E⁽²⁾). The value of HOMO-LUMO energy band gap for NF-MELA-H₂O (monomer and dimer both) is less than NF, showing more chemical reactivity for NF-MELA-H₂O. Chemical reactivity has been described with the assistance of electronic descriptors. Global electrophilicity index (ω = 7.3992 eV) shows that NF-MELA-H₂O behaves as a strong electrophile than NF. The local reactivity descriptors analyses such as Fukui functions, local softnesses and electrophilicity indices performed to determine the reactive sites within NF-MELA-H₂O. In MEP map of NF-MELA (monomer and dimer) electronegative regions are about NO₂ and C=O group of NF, although the electropositive regions are around NH₂, N-H group and H₂O molecule. Molar refractivity (MR) value of NF-MELA-H₂O (monomer and dimer) lies within the range set by Lipinski's modified rules. This study could set as an example to study the H-bond interactions in pharmaceutical cocrystals.

Spectroscopic, thermodynamic and computational evidence of the locations of the FADs in the nitrogen fixation-associated electron transfer flavoprotein

Flavin-based electron bifurcation allows enzymes to redistribute energy among electrons by coupling endergonic and exergonic electron transfer reactions. Diverse bifurcating enzymes employ a two-flavin electron transfer flavoprotein (ETF) that accepts hydride from NADH at a flavin (the so-called bifurcating FAD, Bf-FAD). The Bf-FAD passes one electron exergonically to a second flavin thereby assuming a reactive semiquinone state able to reduce ferredoxin or flavodoxin semiquinone. The flavin that accepts one electron and passes it on via exergonic electron transfer is known as the electron transfer FAD (ET-FAD) and is believed to correspond to the single FAD present in canonical ETFs, in domain II. The Bf-FAD is believed to be the one that is unique to bifurcating ETFs, bound between domains I and III. This very reasonable model has yet to be challenged experimentally. Herein we used site-directed mutagenesis to disrupt FAD binding to the presumed Bf site between domains I and III, in the Bf-ETF from Rhodopseudomonas palustris (RpaETF). The resulting protein contained only 0.80 ± 0.05 FAD, plus 1.21 ± 0.04 bound AMP as in canonical ETFs. The flavin was not subject to reduction by NADH, confirming absence of Bf-FAD. The retained FAD displayed visible circular dichroism (CD) similar to that of the ET-FAD of RpaETF. Likewise, the mutant underwent two sequential one-electron reductions forming and then consuming anionic semiquinone, reproducing the reactivity of the ET-FAD. These data confirm that the retained FAD in domain II corresponds the ET-FAD. Quantum chemical calculations of the absorbance and CD spectra of each of WT RpaETF's two flavins reproduced the observed differences between their CD and absorbance signatures. The calculations for the flavin bound in domain II agreed better with the spectra of the ET-flavin, and those calculated based on the flavin between domains I and III agreed better with spectra of the Bf-flavin. Thus calculations independently confirm the locations of each flavin. We conclude that the site in domain II harbours the ET-FAD whereas the mutated site between domains I and III is the Bf-FAD site, confirming the accepted model by two different tests.

Structural and Reactivity Analyses of Nitrofurantoin–4-dimethylaminopyridine Salt Using Spectroscopic and Density Functional Theory Calculations

Pharmaceutical salt, nitrofurantoin–4-dimethylaminopyridine (NF-DMAP), along with its native components NF and DMAP are scrutinized by FT-IR and FT-Raman spectroscopy along with density functional theory so that an insight into the H-bond patterns in the respective crystalline lattices can be gained. Two different functionals, B3LYP and wB97X-D, have been used to compare the theoretical results. The FT-IR spectra obtained for NF-DMAP and NF clearly validate the presence of C33–H34⋅⋅⋅O4 and N23–H24⋅⋅⋅N9 hydrogen bonds by shifting in the stretching vibration of –NH and –CH group of DMAP+ towards the lower wavenumber side. To explore the significance of hydrogen bonding, quantum theory of atoms in molecules (QTAIM) has been employed, and the findings suggest that the N23–H24⋅⋅⋅N9 bond is a strong intermolecular hydrogen bond. The decrement in the HOMO-LUMO gap, which is calculated from NF → NF-DMAP, reveals that the active pharmaceutical ingredient is chemically less reactive compared to the salt. The electrophilicity index (ω) profiles for NF and DMAP confirms that NF is acting as electron acceptor while DMAP acts as electron donor. The reactive sites of the salt are plotted by molecular electrostatic potential (MEP) surface and calculated using local reactivity descriptors.

Fluorescence Properties of Flavin Semiquinone Radicals in Nitronate Monooxygenase

Fluorescent cofactors like flavins can be exploited to probe their local environment with spatial and temporal resolution. Although the fluorescence properties of the oxidized and two-electron-reduced states of flavins have been studied extensively, this is not the case for the one-electron-reduced state. Both the neutral and anionic semiquinones have proven particularly challenging to examine, as they are unstable in solution and are transient, short-lived species in many catalytic cycles. Here, we report that the nitronate monooxygenase (NMO) from Pseudomonas aeruginosa PAO1 is capable of stabilizing both semiquinone forms anaerobically for hours, thus enabling us to study their spectroscopy in a constant protein environment. We found that in the active site of NMO, the anionic semiquinone exhibits no fluorescence, whereas the neutral semiquinone radical shows a relatively strong fluorescence, with a behavior that violates the Kasha-Vavilov rule. These fluorescence properties are discussed in the context of time-dependent density functional theory calculations, which reveal low-lying dark states in both systems.

Electronic spectra of flavin in different redox and protonation states: a computational perspective on the effect of the electrostatic environment

This study discusses how UV/vis absorption spectra of flavin in different redox and protonation states are shifted by the nearby electrostatic microenvironment.

Spectroscopic and molecular structure (monomeric and dimeric model) investigation of Febuxostat: A combined experimental and theoretical study

Febuxostat (FXT) is a urate-lowering drug and xanthine oxidase inhibitor which is used for the treatment of hyperuricemia and gout caused by increased levels of uric acid in the blood (hyperuricemia). The present study aims to provide deeper knowledge of the structural, vibrational spectroscopic and physiochemical properties of FXT based on monomeric and dimeric model with the aid of combination of experimental and computational methods. The conformational analysis of form Q has been done to predict the possible structure of unknown form A. Vibrational spectra of form A and Q has been compared to get an idea of hydrogen bonding interactions of form A. A computational study of FXT has been executed at different level (B3LYP, M06-2X, WB97XD) of theory and 6-31 G (d, p) basis set for dimeric model to elucidate the nature of intermolecular hydrogen bond. The red shift observed in the stretching modes of OH, CO groups and blue shift in stretching mode of CN group in experimental as well as in theoretical spectra explains the involvement of these groups in intermolecular hydrogen bonding. NBO analysis shows that change in electron density (ED) in the lone pair orbital to σ* antibonding orbital (LP1 (N39) → σ* (O3-H38)) with maximum value of E(2) energy confirms the presence of hydrogen bond (N39⋯H38-O3) leading to dimer formation. Study of topological parameters was executed for dimer using Bader's atoms in molecules (AIM) theory predicting the partially covalent nature of hydrogen bonds present in the molecule. The study of molecular electrostatic potential surface (MEPS) map ascertains that the CO, CN group are prone to electrophilic attack and OH group is active towards nucleophilic attack. The lower energy band gap and higher value of softness of dimeric model of FXT indicates its more reactivity, polarisability than monomeric model. The local reactivity descriptors predict the order of reactive sites towards electrophilic, nucleophilic and radical attack. An investigation made to determine the ligand protein interaction of FXT through docking with different molecular targets reveals the inhibitive as well as antibacterial nature of FXT.

A computational study of the vibrationally-resolved electronic circular dichroism spectra of single-chain transoid and cisoid oligothiophenes in chiral conformations

We simulate the vibronic profile of the electronic circular dichroism (ECD) spectra of oligothiophenes in cisoid and transoid chiral arrangements. We consider oligomers of different lengths, from two to fifteen units, and investigate extensively how the ECD spectral shapes depend on the inter-ring torsions. In general, the molecular structures we consider are not stationary points of the ground state potential energy surface. Therefore, in order to perform vibronic calculations, we present a new computational protocol able to define reduced-dimensionality models where the effect of the off-equilibrium modes is removed. This is done adopting a description of the vibrational motions in curvilinear internal coordinates, and vertical harmonic models coupled with an iterative application of projectors to define energy Hessians, and therefore effective normal modes, in the space complementary to the one of the off-equilibrium coordinates. Although we consider both Franck-Condon and Herzberg-Teller contributions, the results show that transoid twisted ribbons always give rise to monosignated ECD spectra, while bi-signated and multi-signated spectra are expected for cisoid helices. These findings are explained on the basis of the different transition strengths of the lowest electronic states imparted by the different spatial arrangement, that is almost linear for transoid structures and more globular for cisoid ones. We predicted the chiroptical response of a large number of possible molecular arrangements. These data are employed to critically discuss the experimental ECD of polythiophenes in different experimental conditions, forming either aggregates or host-guest complexes. The method here proposed to perform vibronic calculations in reduced-dimensionality models is of general applicability and its potential interest goes beyond the practical application presented here.

Absorption and luminescence spectroscopy of mass-selected flavin adenine dinucleotide mono-anions

We report the absorption profile of isolated Flavin Adenine Dinucleotide (FAD) mono-anions recorded using photo-induced dissociation action spectroscopy. In this charge state, one of the phosphoric acid groups is deprotonated and the chromophore itself is in its neutral oxidized state. These measurements cover the first four optical transitions of FAD with excitation energies from 2.3 to 6.0 eV (210-550 nm). The S₀ → S₂ transition is strongly blue shifted relative to aqueous solution, supporting the view that this transition has a significant charge-transfer character. The remaining bands are close to their solution-phase positions. This confirms that the large discrepancy between quantum chemical calculations of vertical transition energies and solution-phase band maxima cannot be explained by solvent effects. We also report the luminescence spectrum of FAD mono-anions in vacuo. The gas-phase Stokes shift for S₁ is 3000 cm^-1, which is considerably larger than any previously reported for other molecular ions and consistent with a significant displacement of the ground and excited state potential energy surfaces. Consideration of the vibronic structure is thus essential for simulating the absorption and luminescence spectra of flavins.

Elucidation of the intrinsic optical properties of hydrogen-bonded and protonated flavin chromophores by photodissociation action spectroscopy

The intrinsic optical properties of the flavin chromophore when engaged in hydrogen bonding or being protonated were elucidated by photo-induced action spectroscopy and computations.

Franck-Condon Models for Simulating the Band Shape of Electronic Absorption Spectra

Band shape is an essential ingredient in the simulation of electronic absorption spectra. The excitation of multiple series of vibrational levels during an electronic excitation is a main contributor to band shapes. Here we present two simple models based on the Franck-Condon displaced-harmonic-oscillator model. The models are both derived from the time-dependent formulation of electronic spectroscopy. They assume that the transition dipoles do not depend on geometry and that the potential energy surfaces are locally quadratic; one model is second order in time and is called LQ2, and the other is third order in time and is called LQ3. These models are suitable for simulating the unresolved vibronic band shapes of electronic spectra that involve many vibrational modes. The models are straightforward and can be easily applied to simulate absorption spectra that are composed of many electronic transitions. As compared to carrying out molecular dynamics simulations, they require relatively few electronic structure calculations, and the additional cost for constructing the spectra is negligible. Therefore, the models are suitable for simulating the spectra of complex systems such as transition-metal complexes.

Theoretical Assessment of Excited State Gradients and Resonance Raman Intensities for the Azobenzene Molecule

The ground state geometries and vibrational frequencies as well as the excitation energies and excited state gradients of the S₁(nπ*) and S₂(ππ*) states of trans- and cis-azobenzene are investigated by several DFT methods, namely B3LYP, PBE, M06-2X, CAM-B3LYP, and ωB97X. Excited state properties and in particular gradients are also assessed using the wave function based methods EOM-CCSD and RASPT2/RASSCF. Comparison with experimental data shows that the B3LYP functional gives the most accurate results for the ground state geometry and vibrational frequencies. The analysis of the vertical excitation energies reveals that the RASPT2 approach provides the most accurate excitation energies with deviations of the order of 0.1 eV. Among the TDDFT methods, the CAM-B3LYP functional shows the best performance on the excitation energies. By assessing the excited state gradients with respect to the reference RASPT2 data, the most accurate gradients are obtained with B3LYP, whereas other functionals as well as the EOM-CCSD and RASSCF calculations give less consistent results. Overall, despite the tendency of B3LYP to underestimate the excitation energies, this functional provides the most balanced description of both ground and excited state properties for both isomers of azobenzene in the Franck-Condon region.

Advances in Quantum Mechanochemistry: Electronic Structure Methods and Force Analysis

In quantum mechanochemistry, quantum chemical methods are used to describe molecules under the influence of an external force. The calculation of geometries, energies, transition states, reaction rates, and spectroscopic properties of molecules on the force-modified potential energy surfaces is the key to gain an in-depth understanding of mechanochemical processes at the molecular level. In this review, we present recent advances in the field of quantum mechanochemistry and introduce the quantum chemical methods used to calculate the properties of molecules under an external force. We place special emphasis on quantum chemical force analysis tools, which can be used to identify the mechanochemically relevant degrees of freedom in a deformed molecule, and spotlight selected applications of quantum mechanochemical methods to point out their synergistic relationship with experiments.

Importance of Vibronic Effects in the UV-Vis Spectrum of the 7,7,8,8-Tetracyanoquinodimethane Anion

We present a computational method for simulating vibronic absorption spectra in the ultraviolet-visible (UV-vis) range and apply it to the 7,7,8,8-tetracyanoquinodimethane anion (TCNQ^-), which has been used as a ligand in black absorbers. Gaussian broadening of vertical electronic excitation energies of TCNQ^- from linear-response time-dependent density functional theory produces only one band, which is qualitatively incorrect. Thus, the harmonic vibrational modes of the two lowest doublet states were computed, and the vibronic UV-vis spectrum was simulated using the displaced harmonic oscillator approximation, the frequency-shifted harmonic oscillator approximation, and the full Duschinsky formalism. An efficient real-time generating function method was implemented to avoid the exponential complexity of conventional Franck-Condon approaches to vibronic spectra. The obtained UV-vis spectra for TCNQ^- agree well with experiment; the Duschinsky rotation is found to have only a minor effect on the spectrum. Born-Oppenheimer molecular dynamics simulations combined with calculations of the electronic excitation energies for a large number of molecular structures were also used for simulating the UV-vis spectrum. The Born-Oppenheimer molecular dynamics simulations yield a broadening of the energetically lowest peak in the absorption spectrum, but additional vibrational bands present in the experimental and simulated quantum harmonic oscillator spectra are not observed in the molecular dynamics simulations. Our results underline the importance of vibronic effects for the UV-vis spectrum of TCNQ^-, and they establish an efficient method for obtaining vibronic spectra using a combination of linear-response time-dependent density functional theory and a real-time generating function approach.

Revisiting Vertical Models To Simulate the Line Shape of Electronic Spectra Adopting Cartesian and Internal Coordinates

Vertical models for the simulation of spectroscopic line shapes expand the potential energy surface (PES) of the final state around the equilibrium geometry of the initial state. These models provide, in principle, a better approximation of the region of the band maximum. At variance, adiabatic models expand each PES around its own minimum. In the harmonic approximation, when the minimum energy structures of the two electronic states are connected by large structural displacements, adiabatic models can breakdown and are outperformed by vertical models. However, the practical application of vertical models faces the issues related to the necessity to perform a frequency analysis at a nonstationary point. In this contribution we revisit vertical models in harmonic approximation adopting both Cartesian (x) and valence internal curvilinear coordinates (s). We show that when x coordinates are used, the vibrational analysis at nonstationary points leads to a deficient description of low-frequency modes, for which spurious imaginary frequencies may even appear. This issue is solved when s coordinates are adopted. It is however necessary to account for the second derivative of s with respect to x, which here we compute analytically. We compare the performance of the vertical model in the s-frame with respect to adiabatic models and previously proposed vertical models in x- or Q₁-frame, where Q₁ are the normal coordinates of the initial state computed as combination of Cartesian coordinates. We show that for rigid molecules the vertical approach in the s-frame provides a description of the final state very close to the adiabatic picture. For sizable displacements it is a solid alternative to adiabatic models, and it is not affected by the issues of vertical models in x- and Q₁-frames, which mainly arise when temperature effects are included. In principle the G matrix depends on s, and this creates nonorthogonality problems of the Duschinsky matrix connecting the normal modes of initial and final states in adiabatic approaches. We highlight that such a dependence of G on s is also an issue in vertical models, due to the necessity to approximate the kinetic term in the Hamiltonian when setting up the so-called GF problem. When large structural differences exist between the initial and the final-state minima, the changes in the G matrix can become too large to be disregarded.

General formulation of vibronic spectroscopy in internal coordinates

Our general platform integrating time-independent and time-dependent evaluations of vibronic effects at the harmonic level for different kinds of absorption and emission one-photon, conventional and chiral spectroscopies has been extended to support various sets of internal coordinates. Thanks to the implementation of analytical first and second derivatives of different internal coordinates with respect to cartesian ones, both vertical and adiabatic models are available, with the inclusion of mode mixing and, possibly, Herzberg-Teller contributions. Furthermore, all supported non-redundant sets of coordinates are built from a fully automatized algorithm using only a primitive redundant set derived from a bond order-based molecular topology. Together with conventional stretching, bending, and torsion coordinates, the availability of additional coordinates (including linear and out-of-plane bendings) allows a proper treatment of specific systems, including, for instance, inter-molecular hydrogen bridges. A number of case studies are analysed, showing that cartesian and internal coordinates are nearly equivalent for semi-rigid systems not experiencing significant geometry distortions between initial and final electronic states. At variance, delocalized (possibly weighted) internal coordinates become much more effective than their cartesian counterparts for flexible systems and/or in the presence of significant geometry distortions accompanying electronic transitions.

Theoretical investigation of the broad one-photon absorption line-shape of a flexible symmetric carbazole derivative

The one-photon absorption spectrum of a carbazole derivative has been studied by employing density functional response theory combined with a mixed quantum/classical approach to simulate the spectral shape.

A proposal for a dipole-generated BLUF domain mechanism

The resting and signaling structures of the blue-light sensing using flavin (BLUF) photoreceptor domains are still controversially debated due to differences in the molecular models obtained by crystal and NMR structures. Photocycles for the given preferred structural framework have been established, but a unifying picture combining experiment and theory remains elusive. We summarize present work on the AppA BLUF domain from both experiment and theory. We focus on IR and UV/vis spectra, and to what extent theory was able to reproduce experimental data and predict the structural changes upon formation of the signaling state. We find that the experimental observables can be theoretically reproduced employing any structural model, as long as the orientation of the signaling essential Gln63 and its tautomer state are a choice of the modeler. We also observe that few approaches are comparative, e.g., by considering all structures in the same context. Based on recent experimental findings and a few basic calculations, we suggest the possibility for a BLUF activation mechanism that only relies on electron transfer and its effect on the local electrostatics, not requiring an associated proton transfer. In this regard, we investigate the impact of dispersion correction on the interaction energies arising from weakly bound amino acids.