Abstract:Many popular methods for the calculation of chemical potentials rely on the insertion of test particles into the target system. In the case of liquids and liquid mixtures, this procedure increases in difficulty upon increasing density or concentration, and the use of sophisticated enhanced sampling techniques becomes inevitable. In this work we propose an alternative strategy, spatially resolved thermodynamic integration, or SPARTIAN for short. Here, molecules are described with atomistic resolution in a simulation subregion, and as ideal gas particles in a larger reservoir. All molecules are free to diffuse between subdomains adapting their resolution on the fly. To enforce a uniform density profile across the simulation box, a single-molecule external potential is computed, applied, and identified with the difference in chemical potential between the two resolutions. Since the reservoir is represented as an ideal gas bath, this difference exactly amounts to the excess chemical potential of the target system. The present approach surpasses the high density/concentration limitation of particle insertion methods because the ideal gas molecules entering the target system region spontaneously adapt to the local environment. The ideal gas representation contributes negligibly to the computational cost of the simulation, thus allowing one to make use of large reservoirs at minimal expenses. The method has been validated by computing excess chemical potentials for pure Lennard-Jones liquids and mixtures, SPC and SPC/E liquid water, and aqueous solutions of sodium chloride. The reported results well reproduce literature data for these systems.

Extracting free energies from local composition fluctuations in solids: A frequency distribution analysis of simulated atom probe data obtained by TAPSim

Quantifying free energies through local composition fluctuations in binary solid solutions: A proof-of-concept study using atomistic embedded-atom simulations

An Effective Method to Estimate Composition Amplitude of Spinodal Decomposition for Atom Probe Tomography Validated by Phase Field Simulations

Parameter free quantitative analysis of atom probe data by correlation functions: Application to the precipitation in Al-Zn-Mg-Cu

Efficient Approaches To Solutions Of Partition Function For Condensed Matters

Calculations of excess free energies of precipitates via direct thermodynamic integration across phase boundaries

Composition-dependent interatomic potentials: A systematic approach to modelling multicomponent alloys

Uncertainty Quantification of DFT-predicted Finite Temperature Thermodynamic Properties within the Debye Model

Automated free energy calculation from atomistic simulations

Local temperature measurement in molecular dynamics simulations with rigid constraints

Understanding atom probe's analytical performance for iron oxides using correlation histograms and ab initio calculations

A novel method for calculating relative free energy of similar molecules in two environments

First-principles modeling of temperature and concentration dependent solubility in the phase separating Fe$_x$Cu$_{1-x}$ alloy system

Spatially Resolved Thermodynamic Integration: An Efficient Method to Compute Chemical Potentials of Dense Fluids

Solution to the key problem of statistical physics -- calculations of partition function of many-body systems

Uncertainties in interfacial excess calculations from atom probe tomography data

Scalable Ab Initio Electronic Structure Methods with Near Chemical Accuracy for Main Group Chemistry

Atomic Decompositions of Periodic Electronic-Structure Simulations

Crystallographic dependence of Field Evaporation Energy Barrier in metals using Field Evaporation Energy Loss Spectroscopy mapping

Local Crystallography Analysis for Atomically Resolved Scanning Tunneling Microscopy Images