Nadia Rega

Energies, structures, and electronic properties of molecules in solution with the C-PCM solvation model

The conductor‐like solvation model, as developed in the framework of the polarizable continuum model (PCM), has been reformulated and newly implemented in order to compute energies, geometric structures, harmonic frequencies, and electronic properties in solution for any chemical system that can be studied in vacuo. Particular attention is devoted to large systems requiring suitable iterative algorithms to compute the solvation charges: the fast multipole method (FMM) has been extensively used to ensure a linear scaling of the computational times with the size of the solute. A number of test applications are presented to evaluate the performances of the method.

New developments in the polarizable continuum model for quantum mechanical and classical calculations on molecules in solution

The polarizable continuum model (PCM), used for the calculation of molecular energies, structures, and properties in liquid solution has been deeply revised, in order to extend its range of applications and to improve its accuracy. The main changes effect the definition of solute cavities, of solvation charges and of the PCM operator added to the molecular Hamiltonian, as well as the calculation of energy gradients, to be used in geometry optimizations. The procedure can be equally applied to quantum mechanical and to classical calculations; as shown also with a number of numerical tests, this PCM formulation is very efficient and reliable. It can also be applied to very large solutes, since all the bottlenecks have been eliminated to obtain a procedure whose time and memory requirements scale linearly with solute size. The present procedure can be used to compute solvent effects at a number of different levels of theory on almost all the chemical systems which can be studied in vacuo.

Polarizable dielectric model of solvation with inclusion of charge penetration effects

An approximate method, recently proposed to include in continuum solvation models the effects of electronic charge lying outside the solute cavity, has been adapted and implemented in the framework of the polarizable continuum model (PCM). This formulation exploits all the features already developed for the other PCM versions; it provides molecular free energies, gradients and second derivatives with respect to nuclear coordinates. The performances of this method have been tested in comparison with other PCM versions, in particular, we examined the reliability of this technique to reproduce actual volume charge distribution effects, compared to traditional procedures based on Gauss’ Law.

Vibrational computations beyond the harmonic approximation: Performances of the B3LYP density functional for semirigid molecules

The performances of the B3LYP density functional in the computation of harmonic and anharmonic frequencies were tested using 14 standard basis sets of double and triple zeta quality for a set of semirigid molecules containing from 4 to 12 atoms. The quality of the results is assessed by comparison with the most reliable computations available in the literature. The study reveals that the relatively cheap 6‐31+G(d,p) basis set performs a very good job for harmonic frequency calculations and that B3LYP anharmonicities are in close agreement with the reference values irrespective of the basis set used. On these grounds “hybrid force fields” are proposed to achieve the best compromise between computer time and quality of the results.

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

Interpretation of structural properties and dynamic behavior of molecules in solution is of fundamental importance to understand their stability, chemical reactivity, and catalytic action. While information can be gained, in principle, by a variety of spectroscopic techniques, the interpretation of the rich indirect information that can be inferred from the analysis of experimental spectra is seldom straightforward because of the subtle interplay of several different effects, whose specific role is not easy to separate and evaluate. In such a complex scenario, theoretical studies can be very helpful at two different levels: (i) supporting and complementing experimental results to determine the structure of the target molecule starting from its spectral properties; (ii) dissecting and evaluating the role of different effects in determining the observed spectroscopic properties. This is the reason why computational spectroscopy is rapidly evolving from a highly specialized research field into a versatile and widespread tool for the assignment of experimental spectra and their interpretation in terms of chemical physical effects. In such a situation, it becomes important that both computationally and experimentally oriented chemists are aware that new methodological advances and integrated computational strategies are available, providing reliable estimates of fundamental spectral parameters not only for relatively small molecules in the gas phase but also for large and flexible molecules in condensed phases. In this Account, we review the most significant methodological contributions from our research group in this field, and by exploiting some recent results of their application to the computation of IR, UV-vis, NMR, and EPR spectral parameters, we discuss the microscopic mechanisms underlying solvent and vibrational effects on the spectral parameters. After reporting some recent achievements for the study of excited states by first principle quantum mechanical approaches, we focus on the treatment of environmental effects by means of mixed discrete-continuum solvent models and on effective methods for computing vibronic contributions to the spectra. We then discuss some new developments, mainly based on time-dependent approaches, allowing us to go beyond the determination of spectroscopic parameters toward the simulation of line widths and shapes. Although further developments are surely needed to improve the accuracy and effectiveness of several items in the proposed approach, we try to show that the first important steps toward a direct comparison between the results obtained in vitro and those obtained in silico have been made, making easier fruitful crossovers among experiments, computations and theoretical models, which would be decisive for a deeper understanding of the spectral behavior associated with complex systems and processes.

Development and validation of reliable quantum mechanical approaches for the study of free radicals in solution

The hyperfine parameters of a number of representative free radicals have been computed by post‐Hartree–Fock and density functional approaches including averaging effects from large amplitude vibrations and solvent effects through a recent implementation of the polarizable continuum model. Our results show that fully ab initio hyperfine splittings are accurate enough to back the interpretation of experimental data and to allow an unbiased judgement of the role played by electronic, vibrational, and environmental effects in determining the observed value. The very good results obtained by a density functional approach including some Hartree–Fock exchange both for intrinsic values and for solvent shifts pave the route for the investigation of large biologically significant radicals in their natural aqueous medium.

A hybrid explicit/implicit solvation method for first-principle molecular dynamics simulations.

In this work, we present a hybrid explicit/implicit solvation model, well suited for first-principles molecular dynamics simulations of solute-solvent systems. An effective procedure is presented that allows to reliably model a solute with a few explicit solvation shells, ensuring solvent bulk behavior at the boundary with the continuum. Such an approach is integrated with high-level ab initio methods using localized basis functions to perform first-principles or mixed quantum mechanics/molecular mechanics simulations within the extended-Lagrangian formalism. A careful validation of the model along with illustrative applications to solutions of acetone and glycine radical are presented, considering two solvents of different polarity, namely, water and chloroform. Results show that the present model describes dynamical and solvent effects with an accuracy at least comparable to that of conventional approaches based on periodic boundary conditions.

Fluorescence Lifetimes and Quantum Yields of Rhodamine Derivatives: New Insights from Theory and Experiment

Although lifetimes and quantum yields of widely used fluorophores are often largely characterized, a systematic approach providing a rationale of their photophysical behavior on a quantitative basis is still a challenging goal. Here we combine methods rooted in the time-dependent density functional theory and fluorescence lifetime imaging microscopy to accurately determine and analyze fluorescence signatures (lifetime, quantum yield, and band peaks) of several commonly used rhodamine and pyronin dyes. We show that the radiative lifetime of rhodamines can be correlated to the charge transfer from the phenyl toward the xanthene moiety occurring upon the S(0) ← S(1) de-excitation, and to the xanthene/phenyl relative orientation assumed in the S(1) minimum structure, which in turn is variable upon the amino and the phenyl substituents. These findings encourage the synergy of experiment and theory as unique tool to design finely tuned fluorescent probes, such those conceived for modern optical sensors.

A quantum mechanical/molecular dynamics/mean field study of acrolein in aqueous solution: Analysis of H bonding and bulk effects on spectroscopic properties

A novel molecular dynamics methodology recently proposed by our group [Rega et al., Chem. Phys. Lett. 422, 367 (2006)], which is based on an integrated hybrid potential rooted in high level quantum mechanical methods using localized basis functions and nonperiodic boundary conditions, has been applied to study acrolein in aqueous solution. The solute structural rearrangement and its hydrogen-bonding pattern due to the interactions with water have been analyzed in some detail. Moreover, the solvent effects on the UV n-->pi* vertical transition and on the NMR 13C and 17O shielding constants of acrolein have been investigated theoretically by performing a posteriori quantum mechanical calculations on a statistically significant number of snapshots extracted from both gas-phase and aqueous solution simulations. Results show that such effective computational strategy can be successfully used to improve our understanding, at atomic level, of important spectroscopic observables.

Exploring the Metric of Excited State Proton Transfer Reactions

The excited state proton transfer (ESPT) reaction taking place between 7-hydroxy-4-(trifluorometyl)coumarin and 1-methylimidazole, recently experimentally characterized, has been here considered as a case study to illustrate the possibility of using theoretical approaches rooted in density functional theory (DFT) and time-dependent DFT (TD-DFT) for the description of complex reactions occurring at the excited state. In particular, beside identifying all stable species occurring at the ground and excited state during the ESPT reaction, a quantitative characterization of their photophysical properties, such as absorption and emission, is obtained by properly including solvent effects. More interestingly, a computational protocol enabling one to locate possible reaction pathways for the ESPT is here proposed. This protocol is based on the use of density based indices purposely developed to characterize the properties of vertical and relaxed excited states which allow one to discriminate the most favorable reaction paths on potential energy surfaces that are in the case of ESPT intrinsically very flat and difficult to characterize based on sole energy criteria, thus opening a new scenario for the description of photoinduced proton transfer reactions.