Quantum mechanical computations and spectroscopy: from small rigid molecules in the gas phase to large flexible molecules in solution

Photophysics of a nucleic acid-protein crosslinking model strongly depends on solvation dynamics: an experimental and theoretical study

We present a combined experimental and theoretical study of the photophysics of 5-benzyluracil (5BU) in methanol, which is a model system for interactions between nucleic acids and proteins. A molecular dynamics study of 5BU in solution through efficient DFT-based hybrid ab initio potentials revealed a remarkable conformational flexibility - allowing the population of two main conformers - as well as specific solute-solvent interactions, which both appear as relevant factors for the observed 5BU optical absorption properties. The simulated absorption spectrum, calculated on such an ensemble, enabled a molecular interpretation of the experimental UV-Vis lowest energy band, which is also involved in the induced photo-reactivity upon irradiation. In particular, the first two excited states (mainly involving the uracil moiety) both contribute to the 5BU lowest energy absorption. Moreover, as a key finding, the nature and brightness of such electronic transitions are strongly influenced by 5BU conformation and the microsolvation of its heteroatoms.

Comprehensive Analysis of Methyl-β-D-ribofuranoside: A Multifaceted Spectroscopic and Theoretical Approach

This study presents a comprehensive analysis of the vibrational spectra of methyl-β-D-ribofuranoside. Employing a combination of inelastic neutron scattering, Raman, and infrared spectroscopy allows for the observation of all modes regardless of the selection rules. The experimental techniques were complemented by density functional theory computational methods using both gas-phase (Gaussian) and solid-state (CRYSTAL, CASTEP) approaches to provide an unambiguous assignment of the defining vibrational features. Two distinct structures of the molecule were identified in the unit cell, differentiated mainly by the orientation of the furanose ring O-H bonds. The low-energy region of the spectrum (<400 cm-1) is dominated by lattice vibrations and functional group rotation, while the midenergy region is dominated by out-of-plane bending motions of the furanose ring (400-900 cm-1) and by C-H bending in the methyl and methylene groups (1400-1600 cm-1). The high-energy region (>2800 cm-1) encompasses the C-H and O-H stretching modes and offers convincing evidence of at least one H-bonding interaction between the two structures of methyl-β-D-ribofuranoside.

Organic compounds for solid state luminescence enhancement/aggregation induced emission: a theoretical perspective

Organic luminophores displaying one or more forms of luminescence enhancement in solid state are extremely promising for the development and performance optimization of functional materials essential to many modern key technologies. Yet, the effort to harness their huge potential is riddled with hurdles that ultimately come down to a limited understanding of the interactions that result in the diverse molecular environments responsible for the macroscopic response. In this context, the benefits of a theoretical framework able to provide mechanistic explanations to observations, supported by quantitative predictions of the phenomenon, are rather apparent. In this perspective, we review some of the established facts and recent developments about the current theoretical understanding of solid-state luminescence enhancement (SLE) with an accent on aggregation-induced emission (AIE). A description of the macroscopic phenomenon and the questions it raises is accompanied by a discussion of the approaches and quantum chemistry methods that are more apt to model these molecular systems with the inclusion of an accurate yet efficient simulation of the local environment. A sketch of a general framework, building from the current available knowledge, is then attempted via the analysis of a few varied SLE/AIE molecular systems from literature. A number of fundamental elements are identified offering the basis for outlining design rules for molecular architectures exhibiting SLE that involve specific structural features with the double role of modulating the optical response of the luminophores and defining the environment they experience in solid state.

An Expedited Route to Optical and Electronic Properties at Finite Temperature via Unsupervised Learning

Electronic properties and absorption spectra are the grounds to investigate molecular electronic states and their interactions with the environment. Modeling and computations are required for the molecular understanding and design strategies of photo-active materials and sensors. However, the interpretation of such properties demands expensive computations and dealing with the interplay of electronic excited states with the conformational freedom of the chromophores in complex matrices (i.e., solvents, biomolecules, crystals) at finite temperature. Computational protocols combining time dependent density functional theory and ab initio molecular dynamics (MD) have become very powerful in this field, although they require still a large number of computations for a detailed reproduction of electronic properties, such as band shapes. Besides the ongoing research in more traditional computational chemistry fields, data analysis and machine learning methods have been increasingly employed as complementary approaches for efficient data exploration, prediction and model development, starting from the data resulting from MD simulations and electronic structure calculations. In this work, dataset reduction capabilities by unsupervised clustering techniques applied to MD trajectories are proposed and tested for the ab initio modeling of electronic absorption spectra of two challenging case studies: a non-covalent charge-transfer dimer and a ruthenium complex in solution at room temperature. The K-medoids clustering technique is applied and is proven to be able to reduce by ∼100 times the total cost of excited state calculations on an MD sampling with no loss in the accuracy and it also provides an easier understanding of the representative structures (medoids) to be analyzed on the molecular scale.

Spatial Decay and Limits of Quantum Solute-Solvent Interactions

Molecular excitations in the liquid-phase environment are renormalized by the surrounding solvent molecules. Herein, we employ the GW approximation to investigate the solvation effects on the ionization energy of phenol in various solvent environments. The electronic effects differ by up to 0.4 eV among the five investigated solvents. This difference depends on both the macroscopic solvent polarizability and the spatial decay of the solvation effects. The latter is probed by separating the electronic subspace and the GW correlation self-energy into fragments. The fragment correlation energy decays with increasing intermolecular distance and vanishes at ∼9 Å, and this pattern is independent of the type of solvent environment. The 9 Å cutoff defines an effective interacting volume within which the ionization energy shift per solvent molecule is proportional to the macroscopic solvent polarizability. Finally, we propose a simple model for computing the ionization energies of molecules in an arbitrary solvent environment.

Understanding Charge Dynamics in Dense Electronic Manifolds in Complex Environments

Photoinduced charge transfer (CT) excited states and their relaxation mechanisms can be highly interdependent on the environment effects and the consequent changes in the electronic density. Providing a molecular interpretation of the ultrafast (subpicosecond) interplay between initial photoexcited states in such dense electronic manifolds in condensed phase is crucial for improving and understanding such phenomena. Real-time time-dependent density functional theory is here the method of choice to observe the charge density, explicitly propagated in an ultrafast time domain, along with all time-dependent properties that can be easily extracted from it. A designed protocol of analysis for real-time electronic dynamics to be applied to time evolving electronic density related properties to characterize both in time and in space CT dynamics in complex systems is here introduced and validated, proposing easy to be read cross-correlation maps. As case studies to test such tools, we present the photoinduced charge-transfer electronic dynamics of 5-benzyluracil, a mimic of nucleic acid/protein interactions, and the metal-to-ligand charge-transfer electronic dynamics in water solution of [Ru(dcbpy)2(NCS)2]4-, dcbpy = (4,4'-dicarboxy-2,2'-bipyridine), or "N34-", a dye sensitizer for solar cells. Electrostatic and explicit ab initio treatment of solvent molecules have been compared in the latter case, revealing the importance of the accurate modeling of mutual solute-solvent polarization on CT kinetics. We observed that explicit quantum mechanical treatment of solvent slowed down the charge carriers mobilities with respect to the gas-phase. When all water molecules were modeled instead as simpler embedded point charges, the electronic dynamics appeared enhanced, with a reduced hole-electron distance and higher mean velocities due to the close fixed charges and an artificially increased polarization effect. Such analysis tools and the presented case studies can help to unveil the influence of the electronic manifold, as well as of the finite temperature-induced structural distortions and the environment on the ultrafast charge motions.

MISPR: an open-source package for high-throughput multiscale molecular simulations

Computational tools provide a unique opportunity to study and design optimal materials by enhancing our ability to comprehend the connections between their atomistic structure and functional properties. However, designing materials with tailored functionalities is complicated due to the necessity to integrate various computational-chemistry software (not necessarily compatible with one another), the heterogeneous nature of the generated data, and the need to explore vast chemical and parameter spaces. The latter is especially important to avoid bias in scattered data points-based models and derive statistical trends only accessible by systematic datasets. Here, we introduce a robust high-throughput multi-scale computational infrastructure coined MISPR (Materials Informatics for Structure-Property Relationships) that seamlessly integrates classical molecular dynamics (MD) simulations with density functional theory (DFT). By enabling high-performance data analytics and coupling between different methods and scales, MISPR addresses critical challenges arising from the needs of automated workflow management and data provenance recording. The major features of MISPR include automated DFT and MD simulations, error handling, derivation of molecular and ensemble properties, and creation of output databases that organize results from individual calculations to enable reproducibility and transparency. In this work, we describe fully automated DFT workflows implemented in MISPR to compute various properties such as nuclear magnetic resonance chemical shift, binding energy, bond dissociation energy, and redox potential with support for multiple methods such as electron transfer and proton-coupled electron transfer reactions. The infrastructure also enables the characterization of large-scale ensemble properties by providing MD workflows that calculate a wide range of structural and dynamical properties in liquid solutions. MISPR employs the methodologies of materials informatics to facilitate understanding and prediction of phenomenological structure-property relationships, which are crucial to designing novel optimal materials for numerous scientific applications and engineering technologies.

Direct observation of the solvent organization and nuclear vibrations of [Ru(dcbpy)2(NCS)2]4-, [dcbpy = (4,4'-dicarboxy-2,2'-bipyridine)], via ab initio molecular dynamics

Environmental effects can drastically influence the optical properties and photoreactivity of molecules, particularly in the presence of polar and/or protic solvents. In this work we investigate a negatively charged Ru(II) complex, [Ru(dcbpy)₂(NCS)₂]^4- [dcbpy = (4,4'-dicarboxy-2,2'-bipyridine)], in water solution, since this system belongs to a broader class of transition-metal compounds undergoing upon photo-excitation rapid and complex charge transfer (CT) dynamics, which can be dictated by structural rearrangement and solvent environment. Ab initio molecular dynamics (AIMD) relying on a hybrid quantum/molecular mechanics scheme is used to probe the equilibrium microsolvation around the metal complex in terms of radial distribution functions of the main solvation sites and solvent effects on the overall equilibrium structure. Then, using our AIMD-based generalized normal mode approach, we investigate how the ligand vibrational spectroscopic features are affected by water solvation, also contributing to the interpretation of experimental Infra-Red spectra. Two solvation sites are found for the ligands: the sulfur and the oxygen sites can interact on average with ∼4 and ∼3 water molecules, respectively, where a stronger interaction of the oxygen sites is highlighted. On average an overall dynamic distortion of the C₂ symmetric gas-phase structure was found to be induced by water solvation. Vibrational analysis reproduced experimental values for ligand symmetric and asymmetric stretchings, linking the observed shifts with respect to the gas-phase to a complex solvent distribution around the system. This is the groundwork for future excited-state nuclear and electronic dynamics to monitor non-equilibrium processes of CT excitation in complex environments, such as exciton migration in photovoltaic technologies.

A Joint Venture of Ab Initio Molecular Dynamics, Coupled Cluster Electronic Structure Methods, and Liquid-State Theory to Compute Accurate Isotropic Hyperfine Constants of Nitroxide Probes in Water

The isotropic hyperfine coupling constant (HFCC, Aiso) of a pH-sensitive spin probe in a solution, HMI (2,2,3,4,5,5-hexamethylimidazolidin-1-oxyl, C9H19N2O) in water, is computed using an ensemble of state-of-the-art computational techniques and is gauged against X-band continuous wave electron paramagnetic resonance (EPR) measurement spectra at room temperature. Fundamentally, the investigation aims to delineate the cutting edge of current first-principles-based calculations of EPR parameters in aqueous solutions based on using rigorous statistical mechanics combined with correlated electronic structure techniques. In particular, the impact of solvation is described by exploiting fully atomistic, RISM integral equation, and implicit solvation approaches as offered by ab initio molecular dynamics (AIMD) of the periodic bulk solution (using the spin-polarized revPBE0-D3 hybrid functional), embedded cluster reference interaction site model integral equation theory (EC-RISM), and polarizable continuum embedding (using CPCM) of microsolvated complexes, respectively. HFCCs are obtained from efficient coupled cluster calculations (using open-shell DLPNO-CCSD theory) as well as from hybrid density functional theory (using revPBE0-D3). Re-solvation of "vertically desolvated" spin probe configuration snapshots by EC-RISM embedding is shown to provide significantly improved results compared to CPCM since only the former captures the inherent structural heterogeneity of the solvent close to the spin probe. The average values of the Aiso parameter obtained based on configurational statistics using explicit water within AIMD and from EC-RISM solvation are found to be satisfactorily close. Using either such explicit or RISM solvation in conjunction with DLPNO-CCSD calculations of the HFCCs provides an average Aiso parameter for HMI in aqueous solution at 300 K and 1 bar that is in good agreement with the experimentally determined one. The developed computational strategy is general in the sense that it can be readily applied to other spin probes of similar molecular complexity, to aqueous solutions beyond ambient conditions, as well as to other solvents in the longer run.

Ab-initio molecular dynamics and hybrid explicit-implicit solvation model for aqueous and nonaqueous solvents: GFP chromophore in water and methanol solution as case study

Solute-solvent interactions are proxies for understanding how the electronic density of a chromophore interacts with the environment in a more exhaustive way. The subtle balance between polarization, electrostatic, and non-bonded interactions need to be accurately described to obtain good agreement between simulations and experiments. First principles approaches providing accurate configurational sampling through molecular dynamics may be a suitable choice to describe solvent effects on solute chemical-physical properties and spectroscopic features, such as optical absorption of dyes. In this context, accurate energy potentials, obtained by hybrid implicit/explicit solvation methods along with employing nonperiodic boundary conditions, are required to represent bulk solvent around a large solute-solvent cluster. In this work, a novel strategy to simulate methanol solutions is proposed combining ab initio molecular dynamics, a hybrid implicit/explicit flexible solvent model, nonperiodic boundary conditions, and time dependent density functional theory. As case study, the robustness of the proposed protocol has been gauged by investigating the microsolvation and electronic absorption of the anionic green fluorescent protein chromophore in methanol and aqueous solution. Satisfactory results are obtained, reproducing the microsolvation layout of the chromophore and, as a consequence, the experimental trends shown by the optical absorption in different solvents.

Deconvolution of conformational exchange from Raman spectra of aqueous RNA nucleosides

Ribonucleic acids (RNAs) are key to the central dogma of molecular biology. While Raman spectroscopy holds great potential for studying RNA conformational dynamics, current computational Raman prediction and assignment methods are limited in terms of system size and inclusion of conformational exchange. Here, a framework is presented that predicts Raman spectra using mixtures of sub-spectra corresponding to major conformers calculated using classical and ab initio molecular dynamics. Experimental optimization allowed purines and pyrimidines to be characterized as predominantly syn and anti, respectively, and ribose into exchange between equivalent south and north populations. These measurements are in excellent agreement with Raman spectroscopy of ribonucleosides, and previous experimental and computational results. This framework provides a measure of ribonucleoside solution populations and conformational exchange in RNA subunits. It complements other experimental techniques and could be extended to other molecules, such as proteins and carbohydrates, enabling biological insights and providing a new analytical tool.

Solvent-Dependent Structures of Natural Products Based on the Combined Use of DFT Calculations and 1H-NMR Chemical Shifts

Detailed solvent and temperature effects on the experimental 1H-NMR chemical shifts of the natural products chrysophanol (1), emodin (2), and physcion (3) are reported for the investigation of hydrogen bonding, solvation and conformation effects in solution. Very small chemical shift of │Δδ│ < 0.3 ppm and temperature coefficients │Δδ/ΔΤ│ ≤ 2.1 ppb/K were observed in DMSO-d6, acetone-d6 and CDCl3 for the C(1)-OH and C(8)-OH groups which demonstrate that they are involved in a strong intramolecular hydrogen bond. On the contrary, large chemical shift differences of 5.23 ppm at 298 K and Δδ/ΔΤ values in the range of -5.3 to -19.1 ppb/K between DMSO-d6 and CDCl3 were observed for the C(3)-OH group which demonstrate that the solvation state of the hydroxyl proton is a key factor in determining the value of the chemical shift. DFT calculated 1H-NMR chemical shifts, using various functionals and basis sets, the conductor-like polarizable continuum model, and discrete solute-solvent hydrogen bond interactions, were found to be in very good agreement with the experimental 1H-NMR chemical shifts even with computationally less demanding level of theory. The 1H-NMR chemical shifts of the OH groups which participate in intramolecular hydrogen bond are dependent on the conformational state of substituents and, thus, can be used as molecular sensors in conformational analysis. When the X-ray structures of chrysophanol (1), emodin (2), and physcion (3) were used as input geometries, the DFT-calculated 1H-NMR chemical shifts were shown to strongly deviate from the experimental chemical shifts and no functional dependence could be obtained. Comparison of the most important intramolecular data of the DFT calculated and the X-ray structures demonstrate significant differences for distances involving hydrogen atoms, most notably the intramolecular hydrogen bond O-H and C-H bond lengths which deviate by 0.152 tο 0.132 Å and 0.133 to 0.100 Å, respectively, in the two structural methods. Further differences were observed in the conformation of -OH, -CH3, and -OCH3 substituents.

Computational prediction of chiroptical properties in structure elucidation of natural products

Covering: up to 2019This review covers the current status of the quantum mechanical prediction of chiroptical properties, such as electronic CD and optical rotation, as needed for stereochemical assignments in new natural products. The reliability of the prediction of chiroptical properties is steadily increasing, with a parallel decrease in the required computational resources. Now, quantum mechanical calculations for a medium-sized natural product can be reliably performed by natural product chemists on a mainstream PC. This review is aimed to guide natural product chemists through the numerous steps involved in such calculations. Through a concise, but comprehensive, discussion of the current computational practice, enriched by a few illustrative examples, this review provides readers with the theoretical background and practical knowledge needed to select the most appropriate parameters for performing the calculations, to anticipate possible problems, and to critically evaluate the reliability of their computational results. Common reasons for mistakes are also discussed; in particular, the importance of the correct evaluation of conformational ensembles of flexible molecules (an aspect often overlooked in current research) is stressed.

Understanding the vibrational spectra of crystalline isoniazid: Raman, IR and INS spectroscopy and solid-state DFT study

This work presents a comprehensive spectroscopic analysis of crystalline isoniazid, one of the main drugs in tuberculosis chemotherapy, using a blend of spectroscopic and computational methods. Mid- and far-infrared, Raman, and inelastic neutron scattering spectroscopies, with contribution of isotopic substitution are combined with discrete and periodic DFT quantum chemical calculations. This combined approach successfully reproduces the whole spectral range, allowing a sound assignment of all the vibrational bands. Previous misassignments have been corrected and several spectral features of isoniazid crystal are reported for the first time. Virtues and limitations of the computational approach (periodic and discrete) are also discussed in light of the present state-of-the-art in the field.

Conformational Equilibrium of Cinchonidine in C6D12 Solution. Alternative NMR/DFT Approach

¹H NMR and ¹³C NMR chemical shifts as well as conformation dependent vicinal ¹H-¹H spin-spin coupling constants for cinchonidine in a dilute C₆D₁₂ solution have been measured. These data have been interpreted in detail exploiting the results of the extensive quantum chemistry calculations of molecular geometry and NMR parameters of the molecule, performed using the density functional theory (DFT) B3LYP/6-311++G(2d,p) polarizable continuum model (PCM) level of theory. The experimental values of NMR parameters for cinchonidine have been reproduced very well in terms of parameters calculated for key conformers of this molecule. Simultaneously, the analysis has provided us with a lot of information on conformational equilibrium of cinchonidine in the investigated solution. These findings remain in general agreement with the conclusions of other works, based on NOESY spectra or other physicochemical data. Thus, a careful quantitative interpretation of easily measurable NMR chemical shifts can be an independent and valuable source of structural information even in such complex cases as cinchonidine in solution.

The importance of nuclear quantum effects in spectral line broadening of optical spectra and electrostatic properties in aromatic chromophores

In this work, we examine the importance of nuclear quantum effects on capturing the line broadening and vibronic structure of optical spectra. We determine the absorption spectra of three aromatic molecules indole, pyridine, and benzene using time dependent density functional theory with several molecular dynamics sampling protocols: force-field based empirical potentials, ab initio simulations, and finally path-integrals for the inclusion of nuclear quantum effects. We show that the absorption spectrum for all these chromophores are similarly broadened in the presence of nuclear quantum effects regardless of the presence of hydrogen bond donor or acceptor groups. We also show that simulations incorporating nuclear quantum effects are able to reproduce the heterogeneous broadening of the absorption spectra even with empirical force fields. The spectral broadening associated with nuclear quantum effects can be accounted for by the broadened distribution of chromophore size as revealed by a particle in the box model. We also highlight the role that nuclear quantum effects have on the underlying electronic structure of aromatic molecules as probed by various electrostatic properties.

Vibrational band assignment of 2-ethyl-1-hexanol

The vibrational structure of 2-ethyl-1-hexanol is of great interest because of its industrial and military applications. However, detailed spectral analysis is challenging due to its flexibility. This paper reports a detailed analysis of the gas and liquid phase vibrational spectra of 2-ethyl-1-hexanol using the Fourier transform infrared spectroscopy and Raman experimental data. By performing a detailed exploration of the conformational space in this work, the theoretical spectra reproduced almost all experimental details observed, and assigned internal valence coordinates to all of the experimentally observed bands in the floppy 2-ethyl-1-hexanol molecule. Relative contributions from the various internal valence coordinates to the experimental vibrational bands are directly compared between the liquid phase Raman band and the gas and liquid phase infrared band.

Hydrogen Atomic Positions of O-H···O Hydrogen Bonds in Solution and in the Solid State: The Synergy of Quantum Chemical Calculations with ¹H-NMR Chemical Shifts and X-ray Diffraction Methods

The exact knowledge of hydrogen atomic positions of O-H···O hydrogen bonds in solution and in the solid state has been a major challenge in structural and physical organic chemistry. The objective of this review article is to summarize recent developments in the refinement of labile hydrogen positions with the use of: (i) density functional theory (DFT) calculations after a structure has been determined by X-ray from single crystals or from powders; (ii) ¹H-NMR chemical shifts as constraints in DFT calculations, and (iii) use of root-mean-square deviation between experimentally determined and DFT calculated ¹H-NMR chemical shifts considering the great sensitivity of ¹H-NMR shielding to hydrogen bonding properties.

Solvation of Uracil and Its Derivatives by DMSO: A DFT-Supported 1H NMR and 13C NMR Study

¹H NMR and ¹³C NMR spectra of uracil, thymine, 5-hydroxymethyluracil, 5,6-dihydrouracil, and 5,6-dihydrothymine in DMSO-d₆ solutions have been measured. Additionally, molecular structures as well as NMR parameters of these compounds and their various solvates have been calculated using DFT B3LYP/6-311++G(2d,p) PCM(DMSO) method. The analysis of the chemical shift data for these compounds has shown that, indeed, in DMSO solutions they occur as equilibrium mixtures of free molecules and solvates in which solute and solvent molecules are joined by NH···O or OH···O hydrogen bonds. The populations of particular species present in the solutions have been estimated. Moreover, it has been found that 5,6-dihydrothymine exists in DMSO solution preferentially in conformation with the methyl group occupying the pseudoequatorial position. This finding is based on the molecular energy calculations and remains in full agreement with the interpretation of NMR data and theoretical calculations of NMR parameters.

Investigation of two- and three-bond carbon-hydrogen coupling constants in cinnamic acid based compounds

Two- and three-bond coupling constants (² J_HC and ³ J_HC ) were determined for a series of 12 substituted cinnamic acids using a selective 2D inphase/antiphase (IPAP)-single quantum multiple bond correlation (HSQMBC) and 1D proton coupled ¹³ C NMR experiments. The coupling constants from two methods were compared and found to give very similar values. The results showed coupling constant values ranging from 1.7 to 9.7 Hz and 1.0 to 9.6 Hz for the IPAP-HSQMBC and the direct ¹³ C NMR experiments, respectively. The experimental values of the coupling constants were compared with discrete density functional theory (DFT) calculated values and were found to be in good agreement for the ³ J_HC . However, the DFT method under estimated the ² J_HC coupling constants. Knowing the limitations of the measurement and calculation of these multibond coupling constants will add confidence to the assignment of conformation or stereochemical aspects of complex molecules like natural products. Copyright © 2016 John Wiley & Sons, Ltd.

Quantum chemical and spectroscopic investigations of 4-Hydroxy-7-methyl-1,8-naphthyridine-3-carboxylic acid

The FTIR (4000–400[Formula: see text]cm[Formula: see text]) and the FT-Raman spectra (4000–50[Formula: see text]cm[Formula: see text]) of 4-Hydroxy-7-methyl-1,8-naphthyridine-3-carboxylic acid are recorded and investigated. The spectra are interpreted using anharmonic frequency computations by VPT2, VSCF and PT2-VSCF methods within DFT/6-311G(d,p) framework. The root mean square (RMS) values indicate that VSCF computed frequencies are in close agreement with the observed frequencies. The combination and overtone bands are also identified in the FTIR spectrum. The intermolecular O-H[Formula: see text]O hydrogen bonding interactions are discussed in the dimer structure of the molecule. The magnitudes of the coupling between pair of modes are also computed. The electronic spectra in water and ethanol solvents are analyzed using TD-B3LYP/6-311[Formula: see text]G(d,p) level of theory. Molecular electrostatic potential (MEP) and HOMO-LUMO analysis are also performed.

(1)H NMR spectra of alcohols in hydrogen bonding solvents: DFT/GIAO calculations of chemical shifts

Proton nuclear magnetic resonance (NMR) shifts of aliphatic alcohols in hydrogen bonding solvents have been computed on the basis of density functional theory by applying the gauge-including atomic orbital method to geometry-optimized alcohol/solvent complexes. The OH proton shifts and hydrogen bond distances for methanol or ethanol complexed with pyridine depend very much on the functional employed and very little on the basis set, provided it is sufficiently large to give the correct quasi-linear hydrogen bond geometry. The CH proton shifts are insensitive to both the functional and the basis set. NMR shifts for all protons in several alcohol/pyridine complexes are calculated at the Perdew, Burke and Ernzerhof PBE0/cc-pVTZ//PBE0/6-311 + G(d,p) level in the gas phase. The results correlate with the shifts for the pyridine-complexed alcohols, determined by analysing data from the NMR titration of alcohols against pyridine. More pragmatically, computed shifts for a wider range of alcohols correlate with experimental shifts in neat pyridine. Shifts for alcohols in dimethylsulfoxide, based on the corresponding complexes in the gas phase, correlate well with the experimental values, but the overall root mean square difference is high (0.23 ppm), shifts for the OH, CHOH and other CH protons being systematically overestimated, by averages of 0.42, 0.21 and 0.06 ppm, respectively. If the computed shifts are corrected accordingly, a very good correlation is obtained with a gradient of 1.00 ± 0.01, an intercept of 0.00 ± 0.02 ppm and a root mean square difference of 0.09 ppm. This is a modest improvement on the result of applying the CHARGE programme to a slightly different set of alcohols. Some alcohol complexes with acetone and acetonitrile were investigated both in the gas phase and in a continuum of the relevant solvent.

Computational Investigation on the Spectroscopic Properties of Thiophene Based Europium β-Diketonate Complexes

The adiabatic transition energies from the lowest triplet states of four Europium tris β-diketonate/phenantroline complexes have been determined in vacuo and in dicholomethane solution by the ΔSCF approach at the density functional theory level, using the PBE1PBE and the CAM-B3LYP hybrid functionals. The calculated adiabatic transition energies have been compared with the experimental 0-0 transitions of each complex determined from phosphorescence spectra of the corresponding Gd(3+) complexes and followed by direct comparison between simulated and experimental spectra line shapes. For compound 1, the Eu(TTA)3Phen system, triplet states other than the lowest one and conformational isomers other than the one present in the crystallographic structure have been considered. In the crystallographic structure, this compound presents three quasi-degenerate low energy triplet states, differing for the TTA ligand where the two unpaired electrons are localized and showing close adiabatic transition energies. For compound 1, the lowest triplet states of the four investigated conformational isomers show similar characteristics and close adiabatic transition energies. On the basis of these results, an investigation of compounds 2-4 (Eu(Br-TTA)3Phen, Eu(DTDK)3Phen, and Eu(MeT-TTA)3) has been performed by considering only the isomer present in the crystallographic structure and only the lowest triplet state of each compound. For compounds 1-3, the energies of the lowest triplet states calculated by both functionals in solution including zero-point energy corrections well reproduce the experimental trends as well as the values of the adiabatic transition energies: CAM-B3LYP, the best performing functional, provides energies of the lowest triplet state with deviations from experiments lower than 1200 cm(-1). Also, the calculated vibrationally resolved phosphorescence spectra and UV-vis absorptions well reproduce the main features of their experimental counterparts. Significant differences between calculated and experimental results are observed for compound 4, for which difficulties in the experimental determination of the triplet state energy were encountered: our results show that the negligible photoluminescence quantum yield of this compound is due to the fact that the energy of the most stable triplet state is significantly lower than that of the resonance level of the Europium ion, and thus the energy transfer process is prevented. These results confirm the reliability of the adopted computational approach in calculating the energy of the lowest triplet state energy of these systems, a key parameter in the design of new ligands for lanthanide complexes presenting large photoluminescence quantum yields.

Solvatochromic Shift of Brooker's Merocyanine: Hartree-Fock Exchange in Time Dependent Density Functional Calculation and Hydrogen Bonding Effect

The Brooker's merocyanine exhibits a large hypsochromic shift from an apolar aprotic solvent to a polar protic solvent. Quantum chemical calculations have been performed to study the solvatochromism, but there remained a discrepancy between the calculated and experimental solvatochromic shifts. In this paper we evaluate quantum mechanically the excitation energies of the Brooker's merocyanine in water, methanol, acetonitrile, and dichloromethane to investigate what are important factors to accurately model the solvatochromism of the dye by using TDDFT in combination with implicit and explicit solvation models including the PCM, PCMSMD, RISM-SCF-SEDD, and mean-field QM/MM. The results severely depend on the density functional, especially on the amount of Hartree-Fock exchange included in the functional. Furthermore, an explicit description of the solute-solvent hydrogen bonds makes a non-negligible contribution to the shift. The experimental large solvatochromic shift can be accurately reproduced by the TDDFT/RISM-SCF-SEDD and mean-field QM/MM calculations with the LC-BOP functional, although the excitation energies in solutions are considerably overestimated. We also estimated the excitation energies and the solvatochromic shift at the SAC-CI/RISM-SCF-SEDD and mean-field QM/MM level, which are in very good agreement with the experimental values. These results indicate that if an explicit solvent model is used, TDDFT calculations using such a long-range corrected functional can accurately model the solvatochromism. However, an ab initio quantum chemical method including sufficient electron correlation effects is required to reproduce not only the solvatochromism but also the excitation energies in solutions.

Dye chemistry with time-dependent density functional theory

In this perspective, we present an overview of the determination of excited-state properties of "real-life" dyes, and notably of their optical absorption and emission spectra, performed during the last decade with time-dependent density functional theory (TD-DFT). We discuss the results obtained with both vertical and adiabatic (vibronic) approximations, choosing relevant examples for several series of dyes. These examples include reproducing absorption wavelengths of numerous families of coloured molecules, understanding the specific band shape of amino-anthraquinones, optimising the properties of dyes used in solar cells, mimicking the fluorescence wavelengths of fluorescent brighteners and BODIPY dyes, studying optically active biomolecules and photo-induced proton transfer, as well as improving the properties of photochromes.

Relation of molecular structure to Franck-Condon bands in the visible-light absorption spectra of symmetric cationic cyanine dyes

A Franck-Condon (FC) model is used to study the solution-phase absorbance spectra of a series of seven symmetric cyanine dyes having between 22 and 77 atoms. Electronic transition energies were obtained from routine visible-light absorbance and fluorescence emission spectra. Harmonic normal modes were computed using density functional theory (DFT) and a polarizable continuum solvent model (PCM), with frequencies corrected using measured mid-infrared spectra. The model predicts the relative energies of the two major vibronic bands to within 5% and 11%, respectively, and also reproduces structure-specific differences in vibronic band shapes. The bands themselves result from excitation of two distinct subsets of normal modes, one with frequencies between 150 and 625cm(-1), and the other between 850 and 1480cm(-1). Vibronic transitions excite symmetric in-plane bending of the polymethine chain, in-plane bends of the polymethine and aromatic C-H bonds, torsions and deformations of N-alkyl substituents, and in the case of the indocyanines, in-plane deformations of the indole rings. For two dyes, the model predicts vibronic coupling into symmetry-breaking torsions associated with trans-cis photoisomerization.

An ultrafast spectroscopic and quantum mechanical investigation of multiple emissions in push–pull pyridinium derivatives bearing different electron donors

An interconversion between rotamers in S₁ was evidenced for a cationic push–pull pyrenyl derivative through a joint femtosecond fluorescence up-conversion and density functional theoretical study.

Tetrathiafulvalene derivatives as cation sensor: density functional theory investigation of the hyper-Rayleigh scattering first hyperpolarizability

The hyper-Rayleigh scattering first hyperpolarizability of a series of extended tetrathiafulvalene (TTF) and TTF derivatives have been theoretically investigated using density functional theory to explore their use as potential cation sensor.

Modeling EPR parameters of nitrogen containing conjugated radical cations

DFT investigation on conjugated radical cations containing¹⁴N nucleus to obtain accurate isotropic hyperfine coupling constants.

Solvation chemical shifts of perylenic antenna molecules from molecular dynamics simulations

Solvation-induced shifts in molecular properties can be realistically simulated by employing a dynamic model with explicit solvent molecules. In this work, (13)C NMR chemical shifts of various candidate antenna molecules for dye-sensitised solar cells have been studied by using density-functional theory calculations both in vacuo and by employing a dynamic solvation model. The solvent effects were investigated using instantaneous molecular dynamics snapshots containing the antenna molecule and surrounding acetonitrile solvent molecules. Such calculations take into account the main mechanisms of solvation-induced chemical shifts. We have analysed the contributions to the solvent shift due to the solvent susceptibility anisotropy, changes in the density of the virtual orbital space and the accessibility of the excited states to the pronouncedly local magnetic hyperfine operator. We present Lorentzian-broadened chemical shift stick spectra in which a comparison of the in vacuo and dynamic-solvation model results is graphically illustrated. The results show that the solvent-accessible atoms at the perimeter of the solute are influenced by the virtual states of the solvent molecules, which are visible to the hyperfine operators of the perimeter nuclei. This enables efficient coupling of the ground state of the solute to the magnetically allowed excited states, resulting in a positive chemical shift contribution of the perimeter nuclei. As a result of solvation, the chemical shift signals of perimeter nuclei are found to be displaced towards larger chemical shift values, whereas the nuclei of the inner region of the solute molecules show the opposite trend. The solvent susceptibility anisotropy is found to cause a small and practically constant contribution.

Mode-specific vibrational relaxation of photoexcited guanosine 5'-monophosphate and its acid form: a femtosecond broadband mid-IR transient absorption and theoretical study

UV-pump/broadband-mid-IR-probe transient absorption (TA) experiments and ab initio quantum mechanical (QM) calculations were used to investigate the photophysics in heavy water of the neutral and acid forms of guanosine 5'-monophosphate (GMP and GMPD(+), respectively). Excited GMP undergoes ultrafast internal conversion (IC) and returns to the electronic ground state in less than one picosecond with a large amount of excess vibrational energy. The spectroscopic signals are dominated by vibrational cooling - a process in which the solute dissipates vibrational energy to the solvent. For neutral GMP, cooling proceeds with a time constant of 3 ps. Following IC, at least some medium-frequency modes such as the carbonyl stretch and an in-plane ring vibration are excited, suggesting that the vibrational energy distribution is non-statistical. This is consistent with predicted structural changes upon passage through the S1/S0 conical intersection. GMPD(+) differs from GMP by a single deuteron at the N7 position, but has a dramatically longer lifetime of 200 ps. Vibrational cooling of the S1 state of GMPD(+) was monitored via several medium-frequency modes that were assigned using QM calculations. These medium-frequency modes are also vibrationally excited in a non-statistical fashion. Excitation of these modes is in line with the change in geometry at the S1 minimum of GMPD(+) predicted by QM calculations. Furthermore, these modes relax at different rates, fully consistent with QM calculations, which predict that excited vibrational states of the carbonyl stretch couple strongly to the D2O solvent and thus deactivate via intermolecular energy transfer (IET). In contrast, the ring stretch couples strongly to other ring modes of the guanine chromophore and appears to decay via intramolecular vibrational energy redistribution (IVR).

New developments of a multifrequency virtual spectrometer: stereo-electronic, dynamical, and environmental effects on chiroptical spectra

Computational spectroscopy has recently evolved from a field reserved for specialists toward a general tool allowing interpretations and analyses of experimental results. However, the current practice of providing tables of transitions for rigid geometries, possibly tuned by phenomenological broadening, is by far too naive. In order to improve this situation in the last few years we have been developing a general, robust, and user-friendly virtual spectrometer (VS) able to complement experimental studies for complex systems in condensed phases. The VS is based on flexible graphical pre- and postprocessing tools interfaced with general number-crunching software. This last tool is rooted in several electronic structure methodologies (DFT, TD-DFT, post-Hartree-Fock), powerful discrete/continuum models for describing environmental effects, and general vibrational and vibronic models. These last topics are the main focus of this work, which sketches our latest developments related to effective inclusion of anharmonic contributions, together with time-independent and/or time-dependent descriptions of vibronic transitions including Franck-Condon, Herzberg-Teller, and Duschinsky effects. Some test cases are described in some detail with the aim of showing the role of different effects in ruling vibrational (VCD) and electronic (ECD, CPL) chiral spectroscopies.

1H and 13C NMR scaling factors for the calculation of chemical shifts in commonly used solvents using density functional theory

Calculation of NMR chemical shifts and coupling constants using quantum mechanical calculations [density functional theory (DFT)], has become a very popular tool for the determination of conformation and the assignment of stereochemistry within a molecule. We present the scaling factors (linear regression parameters) from 10 DFT methods for 10 commonly used NMR solvents using the same set of reference compounds. The results were compared with the corresponding gas-phase calculations to assess the inclusion of the polarizable continuum model for solvent effects.

Fully anharmonic IR and Raman spectra of medium-size molecular systems: accuracy and interpretation

Computation of full infrared (IR) and Raman spectra (including absolute intensities and transition energies) for medium- and large-sized molecular systems beyond the harmonic approximation is one of the most interesting challenges of contemporary computational chemistry. Contrary to common beliefs, low-order perturbation theory is able to deliver results of high accuracy (actually often better than those issuing from current direct dynamics approaches) provided that anharmonic resonances are properly managed. This perspective sketches the recent developments in our research group toward the development of a robust and user-friendly virtual spectrometer rooted in second-order vibrational perturbation theory (VPT2) and usable also by non-specialists essentially as a black-box procedure. Several examples are explicitly worked out in order to illustrate the features of our computational tool together with the most important ongoing developments.