Advanced search
Start date
Betweenand


Toward a numerically efficient description of bulk-solvated anionic states

Full text
Author(s):
Kiataki, Matheus B. ; Coutinho, Kaline ; Varella, Marcio T. do N.
Total Authors: 3
Document type: Journal article
Source: Journal of Chemical Physics; v. 161, n. 3, p. 10-pg., 2024-07-21.
Abstract

We investigate the vertical electron attachment energy (VAE) of 1-methyl-4-nitroimidazole, a model radiosensitizer, employing quantum mechanics/molecular mechanics (QM/MM) and QM/polarized continuum (QM/PCM) solvation models. We considered the solvent-excluded surface (QM/PCM-SES) and Van der Waals (QM/PCM-VDW) cavities within the PCM framework, the electrostatic embedding QM/MM (EE-QM/MM) model, and the self-consistent sequential QM/MM polarizable electrostatic embedding (scPEE-S-QM/MM) model. Due to slow VAE convergence concerning the number of QM solvent molecules, full QM calculations prove inefficient. Ensemble averages in these calculations do not align with VAEs computed for the representative solute-solvent configuration. QM/MM and QM/PCM calculations show agreement with each other for sufficiently large QM regions, although the QM/PCM-VDW model exhibits artifacts linked to the cavity. QM/MM models demonstrate good agreement between ensemble averages and VAEs calculated with the representative configuration. Notably, the VAE computed with the scPEE-S-QM/MM model achieves faster convergence concerning the number of QM water molecules compared to the EE-QM/MM model, attributed to enhanced efficiency from MM charge polarization in the scPEE-S-QM/MM approach. This emphasizes the importance of QM/classical models with accurate solute-solvent and solvent-solvent mutual polarization for obtaining converged VAEs at a reasonable computational cost. The full-QM approach is very inefficient, while the microsolvation model is inaccurate. Computational savings in QM/MM models result from electrostatic embedding and the representative configuration, with the scPEE-S-QM/MM approach emerging as an efficient tool for describing bulk-solvated anions within the QM/MM framework. Its potential extends to improving transient anion state descriptions in biomolecules and radiosensitizers, especially given the frequent employment of microsolvation models. (AU)

FAPESP's process: 23/07835-2 - Computational Developments in Molecular Modeling of Complex Fluids
Grantee:Matheus Bacigalupo Kiataki
Support Opportunities: Scholarships in Brazil - Technical Training Program - Technical Training
FAPESP's process: 20/16155-7 - Molecular models for positronic, biological e materials sciences
Grantee:Márcio Teixeira do Nascimento Varella
Support Opportunities: Regular Research Grants
FAPESP's process: 14/50983-3 - INCT 2014: complex fluids
Grantee:Antonio Martins Figueiredo Neto
Support Opportunities: Research Projects - Thematic Grants
FAPESP's process: 21/09016-3 - Structural and electronic characterization of biogenic molecules emitted by the Amazon flora
Grantee:Kaline Rabelo Coutinho
Support Opportunities: Regular Research Grants