CPL of Mellein and Related Natural Compounds: Analysis of the ESIPT Phenomenon.

(R)-(-)-Mellein, (3R,4R)-4-hydroxymellein and (3R,4S)-4-hydroxymellein obtained from fungi, i.e. from Diplodia globulosa, were investigated as a class of natural products presenting ESIPT (excited state intramolecular proton transfer) phenomenon, through fluorescence and CPL (circularly polarized luminescence). The study was preceded by the assessment of the absolute configuration through ECD and VCD (electronic and vibrational circular dichroism) spectroscopies in addition to NMR spectra. It is found that ESIPT takes place in these systems very rapidly, and no dual fluorescence has been observed. The experimental study is backed up by TD-DFT calculations of ECD and CPL spectra, plus MD dynamics to follow proton transfer in the excited state and careful analysis of the puckering dynamics of the lactone ring. Deprotonated forms of the three compounds were also investigated by the same chiroptical experimental and theoretical methods, showing how one can find in natural compounds not only biological activity but also biologically compatible sensing probes.

Harmonic and anharmonic vibrational computations for biomolecular building blocks: Benchmarking DFT and basis sets by theoretical and experimental IR spectrum of glycine conformers.

Advanced vibrational spectroscopic experiments have reached a level of sophistication that can only be matched by numerical simulations in order to provide an unequivocal analysis, a crucial step to understand the structure-function relationship of biomolecules. While density functional theory (DFT) has become the standard method when targeting medium-size or larger systems, the problem of its reliability and accuracy are well-known and have been abundantly documented. To establish a reliable computational protocol, especially when accuracy is critical, a tailored benchmark is usually required. This is generally done over a short list of known candidates, with the basis set often fixed a priori. In this work, we present a systematic study of the performance of DFT-based hybrid and double-hybrid functionals in the prediction of vibrational energies and infrared intensities at the harmonic level and beyond, considering anharmonic effects through vibrational perturbation theory at the second order. The study is performed for the six-lowest energy glycine conformers, utilizing available "state-of-the-art" accurate theoretical and experimental data as reference. Focusing on the most intense fundamental vibrations in the mid-infrared range of glycine conformers, the role of the basis sets is also investigated considering the balance between computational cost and accuracy. Targeting larger systems, a broad range of hybrid schemes with different computational costs is also tested.

Mid-IR and CH stretching vibrational circular dichroism spectroscopy to distinguish various sources of chirality: The case of quinophaneoxazoline based ruthenium(II) complexes.

Five diastereomers of ruthenium(II) complexes based on quinolinophaneoxazoline ligands were investigated by vibrational circular dichroism (VCD) in the mid-IR and CH stretching regions. Diastereomers differ in three sources of chirality: the planar chirality of the quinolinophane moiety, the central chirality of an asymmetric carbon atom of the oxazoline ring, and the chirality of the ruthenium atom. VCD, allied to DFT calculations, has been found to be effective in disentangling the various forms of chirality. In particular, a VCD band is identified in the CH stretching region directly connected to the chirality of the metal. The analysis of the calculated VCD spectra is carried out by partitioning the complexes into fragments. The anharmonic analysis is also performed with a recently proposed reduced-dimensionality approach: such treatment is particularly important when examining spectroscopic regions highly perturbed by resonances, like the CH stretching region.

Scaling-up VPT2: A feasible route to include anharmonic correction on large molecules.

Vibrational analysis plays a crucial role in the investigation of molecular systems. Though methodologies like second-order vibrational perturbation theory (VPT2) have paved the way to more accurate simulations, the computational cost remains a difficult barrier to overcome when the molecular size increases. Building upon recent advances in the identification of resonances, we propose an approach making anharmonic simulations possible for large-size systems, typically unreachable by standard means. This relies on the fact that, often, only portions of the whole spectra are of actual interest. Therefore, the anharmonic corrections can be included selectively on subsets of normal modes directly related to the regions of interest. Starting from the VPT2 equations, we evaluate rigorously and systematically the impact of the truncated anharmonic treatment onto simulations. The limit and feasibility of the reduced-dimensionality approach are detailed, starting on a smaller model system. The methodology is then challenged on the IR absorption and vibrational circular dichroism spectra of an organometallic complex in three different spectral ranges.

Combination of Resonance and Non-Resonance Chiral Raman Scattering in a Cobalt(III) Complex.

Resonance Raman optical activity (RROA) spectra with high sensitivity reveal details on molecular structure, chirality, and excited electronic properties. Despite the difficulty of the measurements, the recorded data for the Co(III) complex with S,S-N,N-ethylenediaminedisuccinic acid are of exceptional quality and, coupled with the theory, spectacularly document the molecular behavior in resonance. This includes a huge enhancement of the chiral scattering, contribution of the antisymmetric polarizabilities to the signal, and the Herzberg-Teller effect significantly shaping the spectra. The chiral component is by about one order of magnitude bigger than for an analogous aluminum complex. The band assignment and intensity profile were confirmed by simulations based on density functional and vibronic theories. The resonance was attributed to the S0 →S3 transition, with the strongest signal enhancement of Raman and ROA spectral bands below about 800 cm-1 . For higher wavenumbers, other excited electronic states contribute to the scattering in a less resonant way. RROA spectroscopy thus appears as a unique tool to study the structure and electronic states of absorbing molecules in analytical chemistry, biology, and material science.

Matrix-isolation and cryosolution-VCD spectra of α-pinene as benchmark for anharmonic vibrational spectra calculations.

The inclusion of anharmonicity in vibrational spectral analyis remains associated to small molecular systems with up to a dozen of atoms, with half a dozen of non-hydrogen atoms, typically thesize of propylene oxide. One may see two reasons for this: first of all, larger systems are often thought to be computationally too demanding (high computational costs) for a full anharmonic vibrational analysis. Second, the identification of resonances and their correction is often considered something only expert theoreticians could address because of the lack of unequivocal criteria. In this contribution, we illustrate that resonances can indeed become a complex problem, which can be handled almost transparently thanks to recent advances in vibrational perturbation theory (VPT2). The study also emphasizes the importance and the central role played by experiment in benchmarking novel theoretical approaches. In fact, we herein provide the currently highest resolution VCD spectra available for α- and ß-pinene obtained under matrix-isolation conditions and in liquid Xenon as solvent. They are interpreted by VPT2 with novel tests for the identification of resonances. Hence, the study demonstrates the mutual stimulation of advances in both experimental techniques and computational models.

Perturb-Then-Diagonalize Vibrational Engine Exploiting Curvilinear Internal Coordinates.

The present paper is devoted to the implementation and validation of a second-order perturbative approach to anharmonic vibrations, followed by variational treatment of strong couplings (GVPT2) based on curvilinear internal coordinates. The main difference with respect to the customary Cartesian-based formulation is that the kinetic energy operator is no longer diagonal, and has to be expanded as well, leading to additional terms which have to be taken into proper account. It is, however, possible to recast all the equations as well-defined generalizations of the corresponding Cartesian-based counterparts, thus achieving a remarkable simplification of the new implementation. Particular attention is paid to the treatment of Fermi resonances with significant number of test cases analyzed fully, validating the new implementation. The results obtained in this work confirm that curvilinear coordinates strongly reduce the strength of inter-mode couplings compared to their Cartesian counterparts. This increases the reliability of low-order perturbative treatments for semi-rigid molecules and paves the way toward the reliable representation of more flexible molecules where small- and large-amplitude motions can be safely decoupled and treated at different levels of theory.

An Effective and Automated Processing of Resonances in Vibrational Perturbation Theory Applied to Spectroscopy.

The broader availability of cost-effective methodologies like second-order vibrational perturbational theory (VPT2), also in general-purpose quantum chemical programs, has made the inclusion of anharmonic effects in vibrational calculations easier, paving the way to more accurate simulations. Combined with modern computing hardware, VPT2 can be used on relatively complex molecular systems with dozen of atoms. However, the problem of resonances and their corrections remains a critical pitfall of perturbative methods. Recent works have highlighted the sensitivity of band intensities to even subtle resonance effects, underlying the importance of a correct treatment to predict accurate spectral bandshapes. This aspect is even more critical with chiroptical spectroscopies whose signal is weak. This has motivated the present work in exploring robust methods and criteria to identify resonances not only in energy calculations but also on the transition moments. To study their performance, three molecules of representative sizes ranging from ten to several dozens of atoms were chosen. The impact of resonances, as well as the accuracy achievable once they are properly treated, is illustrated by the changes in spectral bandshapes, including chiroptical spectroscopies.

Scissor-like Face to Face π-π Stacking: A Surprising Preference Induced by the Isocyano Group in the Self-Assembled Dimer of Phenyl Isocyanide.

Phenyl isocyanide has been chosen as a prototype to probe the π-π interaction modulated by the -NC group, which has a chameleonic nature with two main resonance forms showing a triple bond and being carbenoid. The rotational spectroscopic investigation complemented with theoretical analyses indicates that the phenyl isocyanide dimer has a scissor-like configuration controlled by dispersive forces along with the formation of π-π stacking. This is the first rotational spectroscopic evidence, to the best of our knowledge, that the mono-substitution by an -NC group on benzene can activate the meta position in forming noncovalent interactions. This work also provides experimental evidence on the importance of substituent effects in modulating π-π stacked structures, as well as practical proof of a biased interaction behavior of isocyanide-substituted aromatic molecules.

Anharmonic Vibrational Raman Optical Activity of Methyloxirane: Theory and Experiment Pushed to the Limits.

Combining Raman scattering and Raman optical activity (ROA) with computer simulations reveals fine structural and physicochemical properties of chiral molecules. Traditionally, the region of interest comprised fundamental transitions within 200-1800 cm-1. Only recently, nonfundamental bands could be observed as well. However, theoretical tools able to match the observed spectral features and thus assist their assignment are rather scarce. In this work, we present an accurate and simple protocol based on a three-quanta anharmonic perturbative approach that is fully fit to interpret the observed signals of methyloxirane within 150-4500 cm-1. An unprecedented agreement even for the low-intensity combination and overtone transitions has been achieved, showing that anharmonic Raman and ROA spectroscopies can be valuable tools to understand vibrations of chiral molecules or to calibrate computational models.

Anharmonic Aspects in Vibrational Circular Dichroism Spectra from 900 to 9000 cm-1 for Methyloxirane and Methylthiirane.

Vibrational circular dichroism (VCD) spectra and the corresponding IR spectra of the chiral isomers of methyloxirane and of methylthiirane have been reinvestigated, both experimentally and theoretically, with particular attention to accounting for anharmonic corrections, as calculated by the GVPT2 approach. De novo recorded VCD spectra in the near IR (NIR) range regarding CH-stretching overtone transitions, together with the corresponding NIR absorption spectra, were also considered and accounted for, both with the GVPT2 and with the local mode approaches. Comparison of the two methods has permitted us to better describe the nature of active "anharmonic" modes in the two molecules and the role of mechanical and electrical anharmonicity in determining the intensities of VCD and IR/NIR data. Finally, two nonstandard IR/NIR regions have been investigated: the first one about ≈2000 cm-1, involving mostly two-quanta bending mode transitions, the second one between 7000 and 7500 cm-1 involving three-quanta transitions containing CH-stretching overtones and HCC/HCH bending modes.

Accounting for molecular flexibility in photoionization: case of tert-butyl hydroperoxide.

tert-Butyl hydroperoxide (tBuOOH) is a common intermediate in the oxidation of organic compounds that needs to be accurately quantified in complex gas mixtures for the development of chemical kinetic models of low temperature combustion. This work presents a combined theoretical and experimental investigation on the synchrotron-based VUV single photon ionization of gas-phase tBuOOH in the 9.0 - 11.0 eV energy range, including dissociative ionization processes. Computations consist of the determination of the structures, vibrational frequencies and the energetics of neutral and ionic tBuOOH. The Franck-Condon spectrum for the tBuOOH+ (X+) + e- ← tBuOOH (X) + hν transition is computed, where special treatment is undertaken because of the flexibility of tBuOOH, in particular regarding the OOH group. Through comparison of the experimental mass-selected threshold photoelectron spectra with explicitly correlated coupled cluster calculations and Franck-Condon simulations that account for the flexibility of the molecule, an estimation of the ionization energy is given. The appearance energy of the only fragment observed within the above-mentioned energy range, identified as the tert-butyl C4H9+, is also reported. Finally, the signal branching ratio between the parent and the fragment ions is provided as a function of photon energy, essential to quantify tBuOOH in gas-phase oxidation/combustion experiments via advanced mass spectrometry techniques.

Modulation of π character upon complexation captured by molecular rotation spectra.

Two configurations of the furan-CF3Cl complex have been observed by high-resolution rotational spectroscopy. One is characterized by a dominant Cl lone pair⋯π*aromatic interaction and the other is stabilized by a C-Cl···π-CC- halogen bond. This is the first rotational spectroscopic evidence, to the best of our knowledge, that shows how a complexation with a partner like CF3Cl (the weak lone pair belt of Cl, to be more specific) can modulate both the aromatic π* and diene π characters of a heteroaromatic molecule in the formation of non-covalent interactions. The results emphasize the partner molecules' role in modulating the π electron structure and will not only deepen our understanding on non-covalent interactions but also lead to better designs of heteroaromatic-based drugs and materials.

Accurate Quantum Chemical Spectroscopic Characterization of Glycolic Acid: A Route Toward its Astrophysical Detection.

The first step to shed light on the abiotic synthesis of biochemical building blocks, and their further evolution toward biological systems, is the detection of the relevant species in astronomical environments, including earthlike planets. To this end, the species of interest need to be accurately characterized from structural, energetic, and spectroscopic viewpoints. This task is particularly challenging when dealing with flexible systems, whose spectroscopic signature is ruled by the interplay of small- and large-amplitude motions (SAMs and LAMs, respectively) and is further tuned by the conformational equilibrium. In such instances, quantum chemical (QC) calculations represent an invaluable tool for assisting the interpretation of laboratory measurements or even observations. In the present work, the role of QC results is illustrated with reference to glycolic acid (CH2OHCOOH), a molecule involved in photosynthesis and plant respiration and a precursor of oxalate in humans, which has been detected in the Murchison meteorite but not yet in the interstellar medium or in planetary atmospheres. In particular, the equilibrium structure of the lowest-energy conformer is derived by employing the so-called semiexperimental approach. Then, accurate yet cost-effective QC calculations relying on composite post-Hartree-Fock schemes and hybrid coupled-cluster/density functional theory approaches are used to predict the structural and ro-vibrational spectroscopic properties of the different conformers within the framework of the second-order vibrational perturbation theory. A purposely tailored discrete variable representation anharmonic approach is used to treat the LAMs related to internal rotations. The computed spectroscopic data, particularly those in the infrared region, complement the available experimental investigations, thus enhancing the possibility of an astronomical detection of this molecule.

General Perturb-Then-Diagonalize Model for the Vibrational Frequencies and Intensities of Molecules Belonging to Abelian and Non-Abelian Symmetry Groups.

In this paper, we show that the standard second-order vibrational perturbation theory (VPT2) for Abelian groups can be used also for non-Abelian groups without employing specific equations for two- or threefold degenerate vibrations but rather handling in the proper way all the degeneracy issues and deriving the peculiar spectroscopic signatures of non-Abelian groups (e.g., l-doubling) by a posteriori transformations of the eigenfunctions. Comparison with the results of previous conventional implementations shows a perfect agreement for the vibrational energies of linear and symmetric tops, thus paving the route to the transparent extension of the equations already available for asymmetric tops to the energies of spherical tops and the infrared and Raman intensities of molecules belonging to non-Abelian symmetry groups. The whole procedure has been implemented in our general engine for vibro-rotational computations beyond the rigid rotor/harmonic oscillator model and has been validated on a number of test cases.

Unusual binary aggregates of perylene bisimide revealed by their electronic transitions in helium nanodroplets and DFT calculations.

The S1 ← S0 electronic transition of perylene bisimide (PBI) and its binary aggregates were investigated using a combination of helium nanodroplet isolation spectroscopy and computational methods. First, well-resolved vibronic bands of the PBI monomer obtained under the superfluid helium nanodroplet environment were compared to simulated vibronic spectra with anharmonic corrections of the band positions. Second, about ten sets of weaker vibronic bands were observed, which show similar vibronic patterns as that of the PBI monomer and have their band origins red-shifted by about 8 to 218 cm-1. Experimental Poisson curve analyses, performed at the origins of these new sets of bands and the PBI monomer, indicate that the carriers of these weaker red-shifted vibronic bands are binary adducts of PBI. Three types of PBI dimer structures where the electronic transition dipole moments of the two subunits are perpendicular to each other were proposed as possible carriers of these red-shifted vibronic patterns. Extensive vibronic simulations were carried out in a multi-step procedure with TD-DFT, vertical Hessian, and finally adiabatic Hessian approaches. Small red-shifted band origins and very similar vibronic patterns to that of the monomer were predicted for unusual, T-shaped, type I dimer structures and are in close agreement with the experimental data. The combined experimental and theoretical results indicate that the helium nanodroplet environment enables the formation of these unusual T-shaped dimers and stabilizes them.

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.

Unraveling the internal conversion process within the Q-bands of a chlorophyll-like-system through surface-hopping molecular dynamics simulations.

The non-radiative relaxation process within the Q-bands of chlorophylls represents a crucial preliminary step during the photosynthetic mechanism. Despite several experimental and theoretical efforts performed in order to clarify the complex dynamics characterizing this stage, a complete understanding of this mechanism is still far to be reached. In this study, non-adiabatic excited-state molecular dynamic simulations have been performed to model the non-radiative process within the Q-bands for a model system of chlorophylls. This system has been considered in the gas phase and then, to have a more representative picture of the environment, with implicit and mixed implicit-explicit solvation models. In the first part of this analysis, absorption spectra have been simulated for each model in order to guide the setup for the non-adiabatic excited-state molecular dynamic simulations. Then, non-adiabatic excited-state molecular dynamic simulations have been performed on a large set of independent trajectories and the population of the Qx and Qy states has been computed as the average of all the trajectories, estimating the rate constant for the process. Finally, with the aim of investigating the possible role played by the solvent in the Qx-Qy crossing mechanism, an essential dynamic analysis has been performed on the generated data, allowing one to find the most important motions during the simulated dynamics.

Theoretical Investigation of the Circularly Polarized Luminescence of a Chiral Boron Dipyrromethene (BODIPY) Dye.

Over the last decade, molecules capable of emitting circularly polarized light have attracted growing attention for potential technological and biological applications. The efficiency of such devices depend on multiple parameters, in particular the magnitude and wavelength of the peak of emitted light, and also on the dissymmetry factor for chiral applications. In light of these considerations, molecular systems with tunable optical properties, preferably in the visible spectral region, are particularly appealing. This is the case of boron dipyrromethene (BODIPY) dyes, which exhibit large molecular absorption coefficients, have high fluorescence yields, are very stable, both thermally and photochemically, and can be easily functionalized. The latter property has been extensively exploited in the literature to produce chromophores with a wide range of optical properties. Nevertheless, only a few chiral BODIPYs have been synthetized and investigated so far. Using a recently reported axially chiral BODIPY derivative where an axially chiral BINOL unit has been attached to the chromophore unit, we present a comprehensive computational protocol to predict and interpret the one-photon absorption and emission spectra, together with their chiroptical counterparts. From the physico-chemical properties of this molecule, it will be possible to understand the origin of the circularly polarized luminescence better, thus helping to fine-tune the properties of interest. The sensitivity of such processes require accurate results, which can be achieved through a proper account of the vibrational structure in optical spectra. Methodologies to compute vibrationally-resolved electronic spectra can now be applied on relatively large chromophores, such as BODIPYs, but require more extensive computational protocols. For this reason, particular attention is paid in the description of the different steps of the protocol, and the potential pitfalls. Finally, we show how, by means of appropriate tools and approaches, data from intermediate steps of the simulation of the final spectra can be used to obtain further insights into the properties of the molecular system under investigation and the origin of the visible bands.

Role of specific solute-solvent interactions on the photophysical properties of distyryl substituted BODIPY derivatives.

Computational spectroscopy is now a valuable tool to better understand experimental spectroscopic data. Recent advancements in computational methods allow the study of electronic excited states of medium-large molecular systems with greater accuracy, simulating spectra directly comparable with experiments. Besides the system of interest, proper modelling of the environment effects is essential to reach a sufficient level of accuracy. This is especially important for solvents with high polarity or those able to establish specific interactions with the solute. One of the most straightforward and fastest ways to deal with solvent effects is the use of implicit models, like the polarizable continuum model (PCM). However, these models cannot describe specific solute-solvent interactions. In this case, one possible solution is the adoption of a mixed model that includes the solute and few solvent molecules within a PCM cavity. With the aim of investigating the role played by specific solute-solvent interactions, the behaviour of styryl substituted BODIPY systems in methanol has been studied at a theoretical level, considering both implicit and mixed implicit-explicit models. In the first part of our analysis, vibrationally resolved electronic spectra for absorption and emission processes have been simulated, using several sets of coordinates to represent normal modes, and compared with experiments. Then, to verify if specific solute-solvent interactions play a crucial role during the ultrafast intraband relaxation processes, 2D electronic spectra were also simulated, which could provide valuable information on the complex ultrafast dynamics.