Quantum Chemical Microsolvation via Free Energy-Based Solvation Site Identification.
- 1Institute of Materials Chemistry, TU Wien, Getreidemarkt 9, 1060Vienna, Austria.
Abstract
In quantum chemical descriptions of solutions, solvents are treated either implicitly, accounting only for bulk properties, or explicitly by including solvent molecules in the calculations. Although explicit treatments are more accurate, they are often computationally prohibitive. As a cost-effective alternative, microsolvation models have emerged that include only a limited number of solvent molecules to capture the most relevant solute-solvent interactions. However, determining their optimal number, positions, and orientations remains a challenge. Existing approaches typically either build up solvent clusters in a stepwise, order-dependent manner or select solvent molecules based on geometric criteria without directly quantifying interaction strength. Here, we introduce the Free Energy Based Identification of Solvation Sites (FEBISS) 2, a revised microsolvation protocol based on Grid Inhomogeneous Solvation Theory (GIST) applied to molecular dynamics trajectories. The method identifies thermodynamically favorable solvation sites from Free Energy contributions. It automatically places, orients, and ranks solvent molecules by interaction strength, using a statistically averaged, order-independent orientation procedure that is applicable to arbitrary rigid solvents rather than water alone. This approach provides a physically grounded criterion for selecting the number of solvent molecules in microsolvated clusters. We demonstrate the applicability of the approach across a range of systems, from small organic molecules to transition-metal catalysts, and for various solvents. We show that explicit microsolvation affects the electronic structure even in systems traditionally considered weakly solvent-sensitive, highlighting FEBISS 2 as a practical methodology to set up physically meaningful microsolvation clusters for quantum chemical studies.
Related
- Investigating the Electrochemical Double Layer with Quantum‐Chemical Simulations and Implicit Solvation Models
- Quantitative Modelling of Solvation in 2D Materials via Alchemical Decoupling and Potential Energy Surfaces
- Extended QM/Me: Quantum-Mechanical Explicit Solvation at Solid–Liquid Interfaces
- Quantum-classical solvation hydrodynamics: a Hamiltonian modeling framework
- COMPLAX: A Python Tool Package for Automated Microsolvation of Organic Molecules
BibTeX
@article{magenheim2026quantum,
title = {{Quantum Chemical Microsolvation via Free Energy-Based Solvation Site Identification.}},
author = {Lukas B Magenheim and Maren Podewitz},
journal = {Journal of chemical theory and computation},
year = {2026},
month = sep,
doi = {10.1021/acs.jctc.6c00994},
url = {https://europepmc.org/article/MED/42917476},
note = {CC BY},
}
RIS
TY - JOUR
TI - Quantum Chemical Microsolvation via Free Energy-Based Solvation Site Identification.
AU - Lukas B Magenheim
AU - Maren Podewitz
PY - 2026
DA - 2026/09/01
JO - Journal of chemical theory and computation
DO - 10.1021/acs.jctc.6c00994
UR - https://europepmc.org/article/MED/42917476
AB - In quantum chemical descriptions of solutions, solvents are treated either implicitly, accounting only for bulk properties, or explicitly by including solvent molecules in the calculations. Although explicit treatments are more accurate, they are often computationally prohibitive. As a cost-effective alternative, microsolvation models have emerged that include only a limited number of solvent molecules to capture the most relevant solute-solvent interactions. However, determining their optimal number, positions, and orientations remains a challenge. Existing approaches typically either build up solvent clusters in a stepwise, order-dependent manner or select solvent molecules based on geometric criteria without directly quantifying interaction strength. Here, we introduce the Free Energy Based Identification of Solvation Sites (FEBISS) 2, a revised microsolvation protocol based on Grid Inhomogeneous Solvation Theory (GIST) applied to molecular dynamics trajectories. The method identifies thermodynamically favorable solvation sites from Free Energy contributions. It automatically places, orients, and ranks solvent molecules by interaction strength, using a statistically averaged, order-independent orientation procedure that is applicable to arbitrary rigid solvents rather than water alone. This approach provides a physically grounded criterion for selecting the number of solvent molecules in microsolvated clusters. We demonstrate the applicability of the approach across a range of systems, from small organic molecules to transition-metal catalysts, and for various solvents. We show that explicit microsolvation affects the electronic structure even in systems traditionally considered weakly solvent-sensitive, highlighting FEBISS 2 as a practical methodology to set up physically meaningful microsolvation clusters for quantum chemical studies.
SN - 1549-9618
SN - 1549-9626
ER -