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

Giovannini, Tommaso; Lafiosca, Piero; Cappelli, Chiara

doi:10.1021/acs.jctc.7b00776

Cited by 69 publications

(119 citation statements)

References 125 publications

Supporting

Mentioning

112

Contrasting

Order By: Relevance

“…The computed larger solvatochromism for QM/FQFµ is thus probably due to the fact that both GS ans ES dipole moments are increased of almost 50% with respect to QM/FQ. Howevef, it is worth noticing that molecular dipole moments are highly affected by non-electrostatic interactions, in particular repulsion energy terms, as recently reported by some of the present authors 32. The inclusion of such contributions can in principle reduce the molecular dipole moments of both GS and ES.…”

supporting

confidence: 76%

“…The interaction between the QM and MM portions is usually described in terms of electrostatic forces, although approaches to include non-electrostatic QM/MM interactions have been proposed recently. [32][33][34] In order to fully capure the physics of solute-…”

mentioning

confidence: 99%

“…As a consequence, non-electrostatic interactions should reasonably shift QM/FQFµ values towards the experimental data. On the other hand, in light of such considerations, the good performance of QM/FQ for pyridine and pyrimidine might be ascribed to error cancellations between the description of electrostatic effects and the lack of explicit consideration of exchange-repulsion effects [32][33][34]. 5 Summary and ConclusionsIn this work, we have extended the fully polarizable QM/FQFµ approach to the evaluation of excited state energies.…”

mentioning

confidence: 99%

See 2 more Smart Citations

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

Giovannini

Riso

Ambrosetti

et al. 2019

The Journal of Chemical Physics

Self Cite

View full text Add to dashboard Cite

Fully polarizable QM/MM approach based on fluctuating charges and fluctuating dipoles, named QM/FQFµ (J. Chem. Theory Comput. 2019, 15 2233-2245), 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 1 arXiv:1906.03852v1 [physics.chem-ph] 10 Jun 2019 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. IntroductionExcited-state phenomena play a crucial role in many application fields, as for instance photocatalysis, optical information storage and solar cells. In the past decades, theoretical modeling of excited-state properties of molecules in the gas-phase has become a widespread strategy of investigation, 1 giving precious information on, for instance, the nature of the electronic excitation, 2-4 nuclei-electron coupling effects 5-7 and excited state electron dynamics. [8][9][10] However, for electronic phenomena taking place in the condensed phase, 11-19 the interplay between the molecule and its environment can substantially alter the electronic response to external electromagnetic fields. Therefore, any accurate modeling of excited states of solvated systems asks for reliable theoretical approaches to include the effects of the environment at all levels of the excitation phenomenon.Most of the currently available approaches to describe the effects of the external environment on molecular properties belong to the class of the so-called focused models; 20-25 the attention is focused on the molecule and the environment is treated a lower level of sophistication as it modifies, but not determines, the molecular response to the external radiation. In order to keep the atomistic description of the environment, thus substantially overcoming well-known and amply used continuum solvent descriptions, 23,25-29 multiscale QM/Molecular Mechanics (MM) approaches have been developed. 30,31 In such models, the molecule (solute) is treated at the QM level, whereas the environment (solvent) is modelled by means of classical MM force fields (FF). The interaction between the QM and MM portions is usually described in terms of electrostatic forces, although approaches to include non-electrostatic QM/MM interactions have been proposed recently. [32][33][34] In order to fully capure the physics of solute-

show abstract

supporting

confidence: 76%

mentioning

confidence: 99%

mentioning

confidence: 99%

See 1 more Smart Citation

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

Giovannini

Riso

Ambrosetti

et al. 2019

The Journal of Chemical Physics

Self Cite

View full text Add to dashboard Cite

show abstract

“…where V QM ½qðr i Þ is the electrostatic potential due to the QM density of charge at the i-th FQ q i placed at r i . Notice that nonelectrostatic interaction terms, which have been recently proposed by some of us, [63] will not be considered in this work. By exploiting a self-consistent field (SCF) description of the QM moiety, the global QM/MM energy functional reads: [45,46,64] E½P; q; k5trhP1…”

Section: T He Ore Tica L M Ode Lmentioning

confidence: 99%

Simulating vertical excitation energies of solvated dyes: From continuum to polarizable discrete modeling

Giovannini

Macchiagodena

Ambrosetti

et al. 2018

Int J of Quantum Chemistry

Self Cite

View full text Add to dashboard Cite

We present a computational study on the spectroscopic properties of UV-Vis absorbing dyes in water solution. We model the solvation environment by using both continuum and discrete models, with and without polarization, to establish how the physical and chemical properties of the solute-solvent interaction may affect the spectroscopic response of aqueous systems. Seven different compounds were chosen, representing different classes of organic molecules. The classical atomistic description of the solvent molecules was enriched with polarization effects treated by means of the fluctuating charges (FQ) model, propagated to the first-order response function of the quantum-mechanical (QM) solute to include its effects withing the modeling of the electronic excitations of the systems. Results obtained with the QM/FQ model were compared with those from continuum solvation models as well as nonpolarizable atomistic models, and then confronted with the experimental values to determine the accuracy that can be expected with each level of theory. Moreover, a thorough structural analysis using molecular dynamics simulations is provided for each system. K E Y W O R D S excitation energies, QM/FQ, QM/MM, QM/PCM, solvent effects, TD-DFT 1 | I N TR O DU C TI O N One-photon absorption spectroscopy within the UV-Visible range is often the most direct and inexpensive analytical tool that can be used to study the electronic properties of a system. Most commonly, such measurements are carried out on solvated samples, with water being a ubiquitous choice.With the gradual increase in the complexity of the systems under investigation, the correct interpretation of experimental data is increasingly reliant on their calculated ab initio counterparts. Many theoretical models based on quantum mechanics (QM), accompanied by their computational implementations, have been presented over the years offering different levels of compromise between the computational cost and the accuracy of the results. [1][2][3] At present, methods based on density functional theory (DFT) and its time-dependent counterpart (TD-DFT) have become the most popular choice for the simulation of absorption spectra of medium-large organic molecular systems thanks to their versatility stemming from the freedom of choice of density functional and basis set, as well as the favorable scaling with system size which allows their application to increasingly large systems. [1,[4][5][6] Many benchmarks studies have been presented elaborating on the merits and limitations of TD-DFT for the simulation of UV-Vis spectroscopy, as well as on the most appropriate choice of functional and basis set combination for different types of system. [6][7][8][9][10][11][12][13][14][15][16] And though many computational studies are carried out on isolated systems, solvent effects should not be neglected for the presence of the solvation environment can significantly alter the electronic absorption properties of a system, both qualitatively and quantitatively. [17][18][19][20][21][22][23][24][25][2...

show abstract

“…One way of classifying them relies on how they account for polarization and polarizability of the MM subsystem. 29 This includes: induced dipoles, 27,[30][31][32][33][34] Drude oscillators, 35 fluctuating charges, 36,37 charge-on-spring models, 38 or atom-centered multipoles. [39][40][41][42][43] Some of the methods have been extended to utilize QM models beyond DFT.…”

Section: Introductionmentioning

confidence: 99%

Reciprocal Polarizable Embedding with a Transferable H2O Potential Function II: Application to (H2O)n Clusters & Liquid Water

Dohn

Jónsson²,

Jónsson³

2019

Preprint

View full text Add to dashboard Cite

The incorporation of polarization in multiscale quantum-mechanics / molecularmechanics (QM/MM) simulations is important for a variety of applications, e.g. chargetransfer reactions. A recently developed formalism based on a density functional theory description of the QM region and a potential energy function for H 2 O molecules that includes quadrupole as well as dipole polarizability of the MM region is used to simulate liquid water and water clusters. Analysis of the energy, atomic forces, MM polarization, and structure is presented. A quantitative assessment of the QM/MM-MM/MM interaction energy differences of all possible QM/MM configurations of (H 2 O) n clusters shows that the interquartile range of the distributions of the QM/MM binding energies are never more than 20 meV/molecule higher or lower than the binding energies produced with either of the single-model results. Comparing these interaction energy differences with the QM/MM induction differences show that they are not systematically caused by the induced MM moments of our polarizable embedding scheme.Optimized hexamer geometries as well as liquid water structure are shown to be improved in comparison with results obtained using point-charge based embedding models neglecting polarization.

show abstract

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

Cited by 69 publications

References 125 publications

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

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

Simulating vertical excitation energies of solvated dyes: From continuum to polarizable discrete modeling

Reciprocal Polarizable Embedding with a Transferable H2O Potential Function II: Application to (H2O)n Clusters & Liquid Water

Contact Info

Product

Resources

About