Kohn-Sham fragment energy decomposition analysis.

We introduce the concept of Kohn-Sham fragment localized molecular orbitals (KS-FLMOs), which are Kohn-Sham molecular orbitals (MOs) localized in specific fragments constituting a generic molecular system. In detail, we minimize the local electronic energies of various fragments, while maximizing the repulsion between them, resulting in the effective localization of the MOs. We use the developed KS-FLMOs to propose a novel energy decomposition analysis, which we name Kohn-Sham fragment energy decomposition analysis, which allows for rationalizing the main non-covalent interactions occurring in interacting systems both in vacuo and in solution, providing physical insights into non-covalent interactions. The method is validated against state-of-the-art energy decomposition analysis techniques and with high-level calculations.

Modeling UV/Vis Absorption Spectra of Food Colorants in Solution: Anthocyanins and Curcumin as Case Studies.

We present a comprehensive computational study of UV/Vis absorption spectra of significant food colorants, specifically anthocyanins and curcumin tautomers, dissolved in polar protic solvents, namely water and ethanol. The absorption spectra are simulated using two fully polarizable quantum mechanical (QM)/molecular mechanics (MM) models based on the fluctuating charge (FQ) and fluctuating charge and dipoles (FQFµ) force fields. To accurately capture the dynamical aspects of the solvation phenomenon, atomistic approaches are combined with configurational sampling obtained through classical molecular dynamics (MD) simulations. The calculated QM/FQ and QM/FQFµ spectra are then compared with experiments. Our findings demonstrate that a precise reproduction of the UV/Vis spectra of the studied pigments can be achieved by adequately accounting for configurational sampling, polarization effects, and hydrogen bonding interactions.

UV-Resonance Raman Spectra of Systems in Complex Environments: A Multiscale Modeling Applied to Doxorubicin Intercalated into DNA.

UV-Resonance Raman (RR) spectroscopy is a valuable tool to study the binding of drugs to biomolecular receptors. The extraction of information at the molecular level from experimental RR spectra is made much easier and more complete thanks to the use of computational approaches, specifically tuned to deal with the complexity of the supramolecular system. In this paper, we propose a protocol to simulate RR spectra of complex systems at different levels of sophistication, by exploiting a quantum mechanics/molecular mechanics (QM/MM) approach. The approach is challenged to investigate RR spectra of a widely used chemotherapy drug, doxorubicin (DOX) intercalated into a DNA double strand. The computed results show good agreement with experimental data, thus confirming the reliability of the computational protocol.

Computational Insights into the Adsorption of Ligands on Gold Nanosurfaces.

We study the adsorption process of model peptides, nucleobases, and selected standard ligands on gold through the development of a computational protocol based on fully atomistic classical molecular dynamics (MD) simulations combined with umbrella sampling techniques. The specific features of the interface components, namely, the molecule, the metallic substrate, and the solvent, are taken into account through different combinations of force fields (FFs), which are found to strongly affect the results, especially changing absolute and relative adsorption free energies and trends. Overall, noncovalent interactions drive the process along the adsorption pathways. Our findings also show that a suitable choice of the FF combinations can shed light on the affinity, position, orientation, and dynamic fluctuations of the target molecule with respect to the surface. The proposed protocol may help the understanding of the adsorption process at the microscopic level and may drive the in-silico design of biosensors for detection purposes.

Linear response properties of solvated systems: a computational study.

We present a computational study of static and dynamic linear polarizabilities in solution. We use different theoretical approaches to describe solvent effects, ranging from quantum mechanics/molecular mechanics (QM/MM) to quantum embedding approaches. In particular, we consider non-polarizable and polarizable QM/MM methods, the latter based on the fluctuating charge (FQ) force field. In addition, we use a quantum embedding method defined in the context of multilevel Hartree-Fock (MLHF), where the system is divided into active and inactive regions, and combine it with a third layer described by means of the FQ model. The multiscale approaches are then used as reference wave functions for equation-of-motion coupled cluster (EOM-CC) response properties, allowing for the account of electron correlation. The developed models are applied to the calculation of linear response properties of two organic moieties-namely, para-nitroaniline and benzonitrile-in non-aqueous solvents-1,4-dioxane, acetonitrile, and tetrahydrofuran. The computed polarizabilities are then discussed in terms of the physico-chemical solute-solvent interactions described by each method (electrostatic, polarization and Pauli repulsion), and finally compared with the available experimental references.

Assessing the quality of QM/MM approaches to describe vacuo-to-water solvatochromic shifts.

The performance of different quantum mechanics/molecular mechanics embedding models to compute vacuo-to-water solvatochromic shifts is investigated. In particular, both nonpolarizable and polarizable approaches are analyzed and computed results are compared to reference experimental data. We show that none of the approaches outperform the others and that errors strongly depend on the nature of the molecular transition to be described. Thus, we prove that the best choice of embedding model highly depends on the molecular system and that the use of a specific approach as a black box can lead to significant errors and, sometimes, totally wrong predictions.

Strong Coupling between Localized Surface Plasmons and Molecules by Coupled Cluster Theory.

Plasmonic nanocavities enable the confinement of molecules and electromagnetic fields within nanometric volumes. As a consequence, the molecules experience a remarkably strong interaction with the electromagnetic field to such an extent that the quantum states of the system become hybrids between light and matter: polaritons. Here, we present a nonperturbative method to simulate the emerging properties of such polaritons: it combines a high-level quantum chemical description of the molecule with a quantized description of the localized surface plasmons in the nanocavity. We apply the method to molecules of realistic complexity in a typical plasmonic nanocavity, featuring also a subnanometric asperity (picocavity). Our results disclose the effects of the mutual polarization and correlation of plasmons and molecular excitations, disregarded so far. They also quantify to what extent the molecular charge density can be manipulated by nanocavities and stand as benchmarks to guide the development of methods for molecular polaritonics.

Ring Vibrations to Sense Anionic Ibuprofen in Aqueous Solution as Revealed by Resonance Raman.

We unravel the potentialities of resonance Raman spectroscopy to detect ibuprofen in diluted aqueous solutions. In particular, we exploit a fully polarizable quantum mechanics/molecular mechanics (QM/MM) methodology based on fluctuating charges coupled to molecular dynamics (MD) in order to take into account the dynamical aspects of the solvation phenomenon. Our findings, which are discussed in light of a natural bond orbital (NBO) analysis, reveal that a selective enhancement of the Raman signal due to the normal mode associated with the C-C stretching in the ring, νC=C, can be achieved by properly tuning the incident wavelength, thus facilitating the recognition of ibuprofen in water samples.

An efficient and robust procedure to calculate absorption spectra of aqueous charged species applied to NO2.

Accurate calculation of absorption spectra of aqueous NO2- requires rigorously sampling the quantum potential energy surfaces for microsolvation of NO2- with at least five explicit water molecules and embedding the resulting clusters in a continuum solvent accounting for the statistical weighted contributions of individual isomers. This method, which we address as ASCEC + PCM, introduces several desired features when compared against MD simulations derived QM/MM spectra: comparatively fewer explicit solvent molecules to be treated with expensive QM methods, the identification of equilibrium structures in the quantum PES to be used in further vibrational spectroscopy, and the unequivocal identification of cluster orbitals undergoing electronic transitions and charge transfer that originate the spectral bands.

Combining multilevel Hartree-Fock and multilevel coupled cluster approaches with molecular mechanics: a study of electronic excitations in solutions.

We investigate the coupling of different quantum-embedding approaches with a third molecular-mechanics layer, which can be either polarizable or non-polarizable. In particular, such a coupling is discussed for the multilevel families of methods, in which the system is divided into an active and an inactive orbital space. The computational cost of the resulting three-layer approaches is reduced by treating the long-range interactions at the classical level. The developed methods are tested by the calculation of excitation energies of molecular systems in aqueous solution, for which an atomistic description of the environment is crucial to correctly reproduce the specific solute-solvent interactions, such as hydrogen bonding. In particular, we present the results obtained for three different moieties - acrolein, pyridine and para-nitroaniline - showing that an almost perfect agreement with experimental data can be achieved when the relevant physico-chemical interactions are included in the modeling of the condensed phase.

A polarizable three-layer frozen density embedding/molecular mechanics approach.

We present a novel multilayer polarizable embedding approach in which the system is divided into three portions, two of which are treated using density functional theory and their interaction is based on frozen density embedding (FDE) theory, and both also mutually interact with a polarizable classical layer described using an atomistic model based on fluctuating charges (FQ). The efficacy of the model is demonstrated by extending the formalism to linear response properties and applying it to the simulation of the excitation energies of organic molecules in aqueous solution, where the solute and the first solvation shell are treated using FDE, while the rest of the solvent is modeled using FQ charges.

Molecular spectroscopy of aqueous solutions: a theoretical perspective.

Computational spectroscopy is an invaluable tool to both accurately reproduce the spectra of molecular systems and provide a rationalization for the underlying physics. However, the inherent difficulty to accurately model systems in aqueous solutions, owing to water's high polarity and ability to form hydrogen bonds, has severely hampered the development of the field. In this tutorial review we present a technique developed and tested in recent years based on a fully atomistic and polarizable classical modeling of water coupled with a quantum mechanical description of the solute. Thanks to its unparalleled accuracy and versatility, this method can change the perspective of computational and experimental chemists alike.

Theory and algorithms for chiroptical properties and spectroscopies of aqueous systems.

Chiroptical properties and spectroscopies are valuable tools to study chiral molecules and assign absolute configurations. The spectra that result from chiroptical measurements may be very rich and complex, and hide much of their information content. For this reason, the interplay between experiments and calculations is especially useful, provided that all relevant physico-chemical interactions that are present in the experimental sample are accurately modelled. The inherent difficulty associated to the calculation of chiral signals of systems in aqueous solutions requires the development of specific tools, able to account for the peculiarities of water-solute interactions, and especially its ability to form hydrogen bonds. In this perspective we discuss a multiscale approach, which we have developed and challenged to model the most used chiroptical techniques.

Absorption spectra of xanthines in aqueous solution: a computational study.

We present a detailed computational analysis of the UV/Vis spectra of caffeine, paraxanthine and theophylline in aqueous solution. A hierarchy of solvation approaches for modeling the aqueous environment have been tested, ranging from the continuum model to the non-polarizable and polarizable quantum mechanical (QM)/molecular mechanics (MM) models, with and without the explicit inclusion of water molecules in the QM portion. The computed results are directly compared with the experimental data, thus highlighting the role of electrostatic, polarization and hydrogen boding solute-solvent interactions.

eT 1.0: An open source electronic structure program with emphasis on coupled cluster and multilevel methods.

The eT program is an open source electronic structure package with emphasis on coupled cluster and multilevel methods. It includes efficient spin adapted implementations of ground and excited singlet states, as well as equation of motion oscillator strengths, for CCS, CC2, CCSD, and CC3. Furthermore, eT provides unique capabilities such as multilevel Hartree-Fock and multilevel CC2, real-time propagation for CCS and CCSD, and efficient CC3 oscillator strengths. With a coupled cluster code based on an efficient Cholesky decomposition algorithm for the electronic repulsion integrals, eT has similar advantages as codes using density fitting, but with strict error control. Here, we present the main features of the program and demonstrate its performance through example calculations. Because of its availability, performance, and unique capabilities, we expect eT to become a valuable resource to the electronic structure community.

Simulating Absorption Spectra of Flavonoids in Aqueous Solution: A Polarizable QM/MM Study.

We present a detailed computational study of the UV/Vis spectra of four relevant flavonoids in aqueous solution, namely luteolin, kaempferol, quercetin, and myricetin. The absorption spectra are simulated by exploiting a fully polarizable quantum mechanical (QM)/molecular mechanics (MM) model, based on the fluctuating charge (FQ) force field. Such a model is coupled with configurational sampling obtained by performing classical molecular dynamics (MD) simulations. The calculated QM/FQ spectra are compared with the experiments. We show that an accurate reproduction of the UV/Vis spectra of the selected flavonoids can be obtained by appropriately taking into account the role of configurational sampling, polarization, and hydrogen bonding interactions.

Interplay between conformational and solvent effects in UV-visible absorption spectra: curcumin tautomers as a case study.

We present a combined theoretical and experimental study on the UV-vis spectra of enol-keto (EK) and keto-keto (KK) tautomeric forms of curcumin dissolved in aqueous solution. Solvent effects have been investigated by resorting to the implicit polarizable continuum model (QM/PCM) and non-polarizable and fully polarizable QM/MM approaches, the latter based on the fluctuating charges (FQ) force-field. In particular, all methods are challenged to rationalize the contribution of conformational, electrostatic and polarization effects in the calculation of the vertical excitation spectra of curcumin tautomers. The obtained results highlight that for both tautomers specific solute-solvent hydrogen-bond interactions play a minor role with respect to conformational and electrostatic effects.

A combined experimental and theoretical study of optical rotatory dispersion for (R)-glycidyl methyl ether in aqueous solution.

The dispersive optical activity for aqueous solutions of non-rigid (R)-glycidyl methyl ether (R-GME) has been explored synergistically from experimental and theoretical perspectives. Density functional theory analyses performed with the polarizable continuum model for implicit solvation identified nine low-lying stable conformers that are interconverted by rotation about two large-amplitude torsional coordinates. The antagonistic chiroptical signatures predicted for these structural isomers were averaged under a Boltzmann-weighting ansatz to estimate the behavior expected for a thermally equilibrated ensemble. This led to optical rotatory dispersion profiles that reproduced the overall shape of observations but failed to achieve uniform agreement with measured specific-rotation values even when anharmonic vibrational corrections were applied. A mixed QM/FQ paradigm, whereby quantum-mechanical (QM) calculations of optical activity were combined with classical molecular dynamics simulations of explicit solvation that included mutual-polarization effects by means of fluctuating charges (FQ), was enlisted to elucidate the microsolvation environment and gauge its impact upon conformer distributions and response properties. Although quantitative accord with experiments remained elusive, this approach revealed strong variations in the magnitude and sign of rotatory powers for R-GME as the configuration of surrounding water molecules evolved, thereby highlighting the inherently dynamical nature of the solvated chiroptical response, calling into question the validity of "static" descriptions based on the presumption of distinct energy minima, and giving insight into the inherent complexity posed by the modeling of such properties for solvated systems.

Effective yet reliable computation of hyperfine coupling constants in solution by a QM/MM approach: Interplay between electrostatics and non-electrostatic effects.

In this paper, we have extended to the calculation of hyperfine coupling constants, the model recently proposed by some of the present authors [Giovannini et al., J. Chem. Theory Comput. 13, 4854-4870 (2017)] to include Pauli repulsion and dispersion effects in Quantum Mechanical/Molecular Mechanics (QM/MM) approaches. The peculiarity of the proposed approach stands in the fact that repulsion/dispersion contributions are explicitly introduced in the QM Hamiltonian. Therefore, such terms not only enter the evaluation of energetic properties but also propagate to molecular properties and spectra. A novel parametrization of the electrostatic fluctuating charge force field has been developed, thus allowing a quantitative reproduction of reference QM interaction energies. Such a parametrization has been then tested against the prediction of EPR parameters of prototypical nitroxide radicals in aqueous solutions.

Electronic transitions for a fully polarizable QM/MM approach based on fluctuating charges and fluctuating dipoles: Linear and corrected linear response regimes.

The fully polarizable Quantum Mechanics/Molecular Mechanics (QM/MM) approach based on fluctuating charges and fluctuating dipoles, named QM/FQFµ [T. Giovannini et al., J. Chem. Theory Comput. 15, 2233 (2019)], is extended to the calculation of vertical excitation energies of solvated molecular systems. Excitation energies are defined within two different solvation regimes, i.e., linear response (LR), where the response of the MM portion is adjusted to the QM transition density, and corrected-Linear Response (cLR) in which the MM response is adjusted to the relaxed QM density, thus being able to account for charge equilibration in the excited state. The model, which is specified in terms of three physical parameters (electronegativity, chemical hardness, and polarizability) is applied to vacuo-to-water solvatochromic shifts of aqueous solutions of para-nitroaniline, pyridine, and pyrimidine. The results show a good agreement with their experimental counterparts, thus highlighting the potentialities of this approach.