Abstract:The structure of aqueous ammonia solutions is investigated through classical molecular dynamics (MD) and ab initio molecular dynamics (AIMD) simulations. We have preliminarily compared three well-known classical force fields for liquid water (SPC, SPC/E, and TIP4P) in order to identify the most accurate one in reproducing AIMD results obtained at the Generalized Gradient Approximation (GGA) and meta-GGA levels of theory. Liquid ammonia has been simulated by implementing an optimized force field recently developed by Chettiyankandy et al. [Fluid Phase Equilib. 511, 112507 (2020)]. Analysis of the radial distribution functions for different ammonia concentrations reveals that the three water force fields provide comparable estimates of the mixture structure, with the SPC/E performing slightly better. Although a fairly good agreement between MD and AIMD is observed for conditions close to the equimolarity, at lower ammonia concentrations, important discrepancies arise, with classical force fields underestimating the number and strength of H-bonds between water molecules and between water and ammonia moieties. Here, we prove that these drawbacks are rooted in a poor sampling of the configurational space spanned by the hydrogen atoms lying in the H-bonds of H2O⋯H2O and, more critically, H2O⋯NH3 neighbors due to the lack of polarization and charge transfer terms. This way, non-polarizable classical force fields underestimate the proton affinity of the nitrogen atom of ammonia in aqueous solutions, which plays a key role under realistic dilute ammonia conditions. Our results witness the need for developing more suited polarizable models that are able to take into account these effects properly.

Molecular dynamics simulations of a potassium ion and an iodide ion in liquid ammonia

In-depth Investigation of the Iodine-Nitrogen Interaction Mechanism Using Ab Initio and Density Functional Theory Calculations

Study on properties of liquid ammonia via molecular dynamics simulation based on ABEEMσπ polarisable force field

Unravelling the structural and dynamical properties of concentrated aqueous ammonium nitrate solutions: MD simulation studies

Rotational Mechanism of Ammonium Ion in Water and Methanol

The structure of water-ammonia mixtures from classical and ab initio molecular dynamics

Molecular dynamics simulation of potassium along the liquid-vapor coexistence curve

Interactions between cyclic nucleotides and common cations: an ab initio molecular dynamics study

The structure of water–ammonia mixtures from classical and ab initio molecular dynamics

Hydration of NH_4^+ in Water: Bifurcated Hydrogen Bonding Structures and Fast Rotational Dynamics

NMR and computational studies of ammonium ion binding to dibenzo-18-crown-6

Similarities and differences between potassium and ammonium ions in liquid water: a first-principles study

Kinetic and dynamic studies of the NH 2 + + H 2 reaction on a high-level ab initio potential energy surface

Jump Rotational Mechanism of Ammonium Ion in Aqueous Solutions

Hydration of NH$_4^+$ in Water: Bifurcated Hydrogen Bonding Structures and Fast Rotational Dynamics

Properties of Hydrogen Bonds in the Protic Ionic Liquid Ethylammonium Nitrate. DFT versus DFTB Molecular Dynamics

Molecular dynamics simulations of ion clustering and conductivity in NaI/ether solutions. I. Effect of ion charge

A Network Approach to Unravel Correlated Ion Pair Dynamics in Protic Ionic Liquids. The Case of Triethylammonium Nitrate

Ammonia-water cation and ammonia dimer cation

Water rotational jump driven large amplitude molecular motions of nitrate ions in aqueous potassium nitrate solution

Ab initio studies of structures and properties of small potassium clusters