Effective fully polarizable QM/MM approaches to compute Raman and Raman Optical Activity spectra in aqueous solution

Raman and Raman Optical Activity (ROA) signals are amply affected by solvent effects, especially in the presence of strongly solute-solvent interactions such as Hydrogen Bonding (HB). In this work, we extend the fully atomistic polarizable Quantum Mechanics/Molecular Mechanics approach, based on the Fluctuating Charges and Fluctuating Dipoles force field to the calculation of Raman and ROA spectra. Such an approach is able to accurately describe specific HB interactions, by also accounting for anisotropic contributions due to the inclusion of fluctuating dipoles. To highlight the potentiality of the novel approach, Raman and ROA spectra of L-Serine and L-Cysteine dissolved in aqueous solution are computed and compared both with alternative theoretical approaches and experimental measurements.

Multiscale Frozen Density Embedding/Molecular Mechanics Approach for Simulating Magnetic Response Properties of Solvated Systems

We present a three-layer hybrid quantum mechanical/quantum embedding/molecular mechanics approach for calculating nuclear magnetic resonance (NMR) shieldings and J-couplings of molecular systems in solution. The model is based on the frozen density embedding (FDE) and polarizable fluctuating charges (FQ) and fluctuating dipoles (FQFμ) force fields and permits the accurate ab initio description of short-range nonelectrostatic interactions by means of the FDE shell and cost-effective treatment of long-range electrostatic interactions through the polarizable force field FQ(Fμ). Our approach's accuracy and potential are demonstrated by studying NMR spectra of Brooker's merocyanine in aqueous and nonaqueous solutions.

QM/Classical Modeling of Surface Enhanced Raman Scattering Based on Atomistic Electromagnetic Models

We present quantum mechanics (QM)/frequency dependent fluctuating charge (QM/ωFQ) and fluctuating dipoles (QM/ωFQFμ) multiscale approaches to model surface-enhanced Raman scattering spectra of molecular systems adsorbed on plasmonic nanostructures. The methods are based on a QM/classical partitioning of the system, where the plasmonic substrate is treated by means of the atomistic electromagnetic models ωFQ and ωFQFμ, which are able to describe in a unique fashion and at the same level of accuracy the plasmonic properties of noble metal nanostructures and graphene-based materials. Such methods are based on classical physics, i.e. Drude conduction theory, classical electrodynamics, and atomistic polarizability to account for interband transitions, by also including an ad-hoc phenomenological correction to describe quantum tunneling. QM/ωFQ and QM/ωFQFμ are thus applied to selected test cases, for which computed results are compared with available experiments, showing the robustness and reliability of both approaches.

Correction to "Going Beyond the Limits of Classical Atomistic Modeling of Plasmonic Nanostructures"

[This corrects the article DOI: 10.1021/acs.jpcc.1c04716.].

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.

Do We Really Need Quantum Mechanics to Describe Plasmonic Properties of Metal Nanostructures?

Optical properties of metal nanostructures are the basis of several scientific and technological applications. When the nanostructure characteristic size is of the order of few nm or less, it is generally accepted that only a description that explicitly describes electrons by quantum mechanics can reproduce faithfully its optical response. For example, the plasmon resonance shift upon shrinking the nanostructure size (red-shift for simple metals, blue-shift for d-metals such as gold and silver) is universally accepted to originate from the quantum nature of the system. Here we show instead that an atomistic approach based on classical physics, ωFQFμ (frequency dependent fluctuating charges and fluctuating dipoles), is able to reproduce all the typical "quantum" size effects, such as the sign and the magnitude of the plasmon shift, the progressive loss of the plasmon resonance for gold, the atomistically detailed features in the induced electron density, and the non local effects in the nanoparticle response. To support our findings, we compare the ωFQFμ results for Ag and Au with literature time-dependent DFT simulations, showing the capability of fully classical physics to reproduce these TDDFT results. Only electron tunneling between nanostructures emerges as a genuine quantum mechanical effect, that we had to include in the model by an ad hoc term.

Absorption Properties of Large Complex Molecular Systems: The DFTB/Fluctuating Charge Approach

We report on the first formulation of a novel polarizable QM/MM approach, where the density functional tight binding (DFTB) is coupled to the fluctuating charge (FQ) force field. The resulting method (DFTB/FQ) is then extended to the linear response within the TD-DFTB framework and challenged to study absorption spectra of large condensed-phase systems.

Going Beyond the Limits of Classical Atomistic Modeling of Plasmonic Nanostructures

Theoretical modeling of plasmonic phenomena is of fundamental importance for rationalizing experimental measurements. Despite the great success of classical continuum modeling, recent technological advances allowing for the fabrication of structures defined at the atomic level require to be modeled through atomistic approaches. From a computational point of view, the latter approaches are generally associated with high computational costs, which have substantially hampered their extensive use. In this work, we report on a computationally fast formulation of a classical, fully atomistic approach, able to accurately describe both metal nanoparticles and graphene-like nanostructures composed of roughly 1 million atoms and characterized by structural defects.

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.

Effective computational route towards vibrational optical activity spectra of chiral molecules in aqueous solution

A polarizable QM/MM approach to accurately compute the Vibrational Optical Activity (VOA) spectra of chiral systems is proposed and applied to aqueous solutions of (l)-methyl lactate and (S)-glycidol.

A General Route to Include Pauli Repulsion and Quantum Dispersion Effects in QM/MM Approaches

A methodology to account for nonelectrostatic interactions in Quantum Mechanical (QM)/Molecular Mechanics (MM) approaches is developed. Formulations for Pauli repulsion and dispersion energy, explicitly depending on the QM density, are derived. Such expressions are based on the definition of an auxiliary density on the MM portion and the Tkatchenko-Scheffler (TS) approach, respectively. The developed method is general enough to be applied to any QM/MM method and partition, provided an accurate tuning of a small number of parameters is obtained. The coupling of the method with both nonpolarizable and fully polarizable QM/fluctuating charge (FQ) approaches is reported and applied. A suitable parametrization for the aqueous solution, so that its most representative features are well reproduced, is outlined. Then, the obtained parametrization and method are applied to calculate the nonelectrostatic (repulsion and dispersion) interaction energy of nicotine in aqueous solution.