Computing chemical potentials of solutions from structure factors

Abstract

The chemical potential of a component in a solution is defined as the free energy change as the amount of the component changes. Computing this fundamental thermodynamic property from atomistic simulations is notoriously difficult, because of the convergence issues in free energy methods and finite size effects. This paper presents the S0 method, which can be used to obtain chemical potentials from static structure factors computed from equilibrium molecular dynamics simulations under the isothermal-isobaric ensemble. The S0 method is demonstrated on the systems of binary Lennard-Jones particles, urea--water mixtures, a NaCl aqueous solution, and a high-pressure carbon-hydrogen mixture.

References in corpus (1)

• Computation of the chemical potential and solubility of amorphous solids

Cited by in corpus (8)

• Diamond formation from hydrocarbon mixtures in planets
• A streamlined molecular-dynamics workflow for computing solubilities of molecular and ionic crystals
• Grand-canonical molecular dynamics simulations powered by a hybrid 4D nonequilibrium MD/MC method: Implementation in LAMMPS and applications to electrolyte solutions
• Theory of Nonequilibrium Multicomponent Coexistence
• Computing chemical potentials of adsorbed or confined fluids
• Combining Molecular Dynamics and Experimental Methods for the Parametrization of Binary Carbonate-Based Electrolytes
• Robust Kirkwood-Buff inversion in complex mixtures via reciprocal-space methods
• Computing ternary liquid phase diagrams: Fe-Cu-Ni