Consistent Density Functional Theory-Based Description of Ion Hydration Through Density-Corrected Many-Body Representations

Francesco Paesani,Etienne Palos,Alessandro Caruso
DOI: https://doi.org/10.26434/chemrxiv-2023-vp6ns-v3
2023-11-02
Abstract:Delocalization error constrains the accuracy of density functional theory (DFT) in describing molecular interactions in ion–water systems. Using Na+ and Cl− in water as model systems, we calculate the effects of delocalization error in the SCAN functional for describing ion–water and water–water interactions in hydrated ions, and demonstrate that density-corrected SCAN (DC-SCAN) predicts n-body and interaction energies with an accuracy approaching coupled cluster theory. The performance of DC-SCAN is size-consistent, maintaining an accurate description of molecular in- teractions well beyond the first solvation shell. Molecular dynamics simulations at ambient conditions with MB-SCAN(DC) potentials, derived from the many-body expansion, predict the solvation structure of Na+ and Cl− in quantitative agreement with reference data, while simultaneously reproducing the structure of liquid water. Beyond rationalizing the accuracy of density-corrected models of ion hydration, our findings suggest that our unified density-corrected many-body formalism holds great promise for efficient DFT-based simulations of condensed-phase systems with chemical accuracy.
Chemistry
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is the influence of the delocalization error in density functional theory (DFT) on the description of molecular interactions in ion - water systems, especially in the application of the hydration process of Na+ and Cl− in water. Specifically, by using the density - corrected SCAN functional (DC - SCAN), the paper aims to improve the prediction accuracy of ion - water and water - water interactions to a level close to that of the coupled - cluster theory and maintain size consistency, so that it can accurately describe molecular interactions beyond the first solvent shell. In addition, the paper also explores the ability of the MB - SCAN(DC) potential energy based on the many - body expansion to predict the hydration structures of Na+ and Cl− under normal temperature conditions while reproducing the structure of liquid water. The core objectives of the paper are: 1. **Evaluate the delocalization error in the SCAN functional**: Analyze the delocalization error in the SCAN functional when describing ion - water and water - water interactions by calculating the interaction energies of Na+ and Cl− in water. 2. **Propose a density - correction method**: Introduce the density - corrected SCAN functional (DC - SCAN) to reduce the prediction deviation caused by the delocalization error. 3. **Verify the accuracy of DC - SCAN**: Through molecular dynamics simulations, verify the prediction ability of DC - SCAN for the hydration structures of Na+ and Cl− and the structure of liquid water under normal temperature conditions, and ensure its consistency with the reference data. 4. **Explore the application of many - body formalism**: Propose a unified density - corrected many - body formalism for efficient and accurate DFT - based condensed - phase system simulations. In summary, this research aims to provide a more accurate and reliable method for simulating and understanding the hydration process of ions in water by improving the DFT method, which is of great significance for the research of natural and industrial chemical phenomena such as biomolecule stabilization, catalysis and transport processes, and electrochemical processes.