We present for the first time a quantum mechanics/molecular mechanics scheme which combines quantum Monte Carlo with the reaction field of classical polarizable dipoles (QMC/MMpol). In our approach, the optimal dipoles are self-consistently generated at the variational Monte Carlo level and then used to include environmental effects in diffusion Monte Carlo. We investigate the performance of this hybrid model in describing the vertical excitation energies of prototypical small molecules solvated in water, namely, methylenecyclopropene and s-trans acrolein. Two polarization regimes are explored where either the dipoles are optimized with respect to the ground-state solute density (polGS) or different sets of dipoles are separately brought to equilibrium with the states involved in the electronic transition (polSS). By comparing with reference supermolecular calculations where both solute and solvent are treated quantum mechanically, we find that the inclusion of the response of the environment to the excitation of the solute leads to superior results than the use of a frozen environment (point charges or polGS), in particular, when the solute-solvent coupling is dominated by electrostatic effects which are well recovered in the polSS condition. QMC/MMpol represents therefore a robust scheme to treat important environmental effects beyond static point charges, combining the accuracy of QMC with the simplicity of a classical approach.

Introducing QMC/MMpol: Quantum Monte Carlo in Polarizable Force Fields for Excited States

FLORIS, FRANCA MARIA;AMOVILLI, CLAUDIO;MENNUCCI, BENEDETTA;
2016-01-01

Abstract

We present for the first time a quantum mechanics/molecular mechanics scheme which combines quantum Monte Carlo with the reaction field of classical polarizable dipoles (QMC/MMpol). In our approach, the optimal dipoles are self-consistently generated at the variational Monte Carlo level and then used to include environmental effects in diffusion Monte Carlo. We investigate the performance of this hybrid model in describing the vertical excitation energies of prototypical small molecules solvated in water, namely, methylenecyclopropene and s-trans acrolein. Two polarization regimes are explored where either the dipoles are optimized with respect to the ground-state solute density (polGS) or different sets of dipoles are separately brought to equilibrium with the states involved in the electronic transition (polSS). By comparing with reference supermolecular calculations where both solute and solvent are treated quantum mechanically, we find that the inclusion of the response of the environment to the excitation of the solute leads to superior results than the use of a frozen environment (point charges or polGS), in particular, when the solute-solvent coupling is dominated by electrostatic effects which are well recovered in the polSS condition. QMC/MMpol represents therefore a robust scheme to treat important environmental effects beyond static point charges, combining the accuracy of QMC with the simplicity of a classical approach.
2016
Guareschi, Riccardo; Zulfikri, Habiburrahman; Daday, Csaba; Floris, FRANCA MARIA; Amovilli, Claudio; Mennucci, Benedetta; Filippi, Claudia
File in questo prodotto:
File Dimensione Formato  
acs.jctc.6b00044.pdf

solo utenti autorizzati

Tipologia: Versione finale editoriale
Licenza: NON PUBBLICO - Accesso privato/ristretto
Dimensione 688.84 kB
Formato Adobe PDF
688.84 kB Adobe PDF   Visualizza/Apri   Richiedi una copia
paper_qmmmpol_GUARESCHI-2.pdf

accesso aperto

Tipologia: Documento in Post-print
Licenza: Tutti i diritti riservati (All rights reserved)
Dimensione 538.63 kB
Formato Adobe PDF
538.63 kB Adobe PDF Visualizza/Apri

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11568/795787
Citazioni
  • ???jsp.display-item.citation.pmc??? 3
  • Scopus 27
  • ???jsp.display-item.citation.isi??? 25
social impact