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.

Ammonia-water cation and ammonia dimer cation

Ab Initio Investigation Of Structure, Stability, Thermal Behavior, Bonding, And Infrared Spectra Of Ionized Water Cluster (H2o)(6)(+) (Vol 145, 154307, 2016)

Ab Initio Investigation of the Lower Energy Candidate Structures for (h2o)(5)+ Water Cluster

Ab Initio Investigation of Possible Candidate Structures and Properties of Water Cluster (H2O)7+ Via Particle Swarm Optimization Method

Possible Low-Energy Isomers of OH (H2O)4 (n = 0, ±1) Clusters Via the Particle Swarm Optimization Algorithm: an Ab Initio Study

Ab Initio Study of Cationic Water Cluster (H2O)(9)(+) Via Particle Swarm Optimization Algorithm

Ab initio molecular orbital calculations of the infrared spectra of hydrogen-bonded complexes of water, ammonia and hydroxylamine

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

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

A Comparison Of Ab Initio Calculation Methods On The Structure Of Nh3 Center Dot Center Dot Center Dot H2o Hydrogen-Bonded Cluster

Dynamics and structural changes of small water clusters on ionization

Understanding the Hydrogen-Bonded Clusters of Ammonia (NH 3 ) n ( n = 3–6): Insights from the Electronic Structure Theory

Theoretical Investigation On The Interactions Between Hrncch And X( X=N-2, O-2, H2o, Nh3)

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

Structures and Interactions in N-methylacetamide-water Mixtures Studied by IR Spectra and Density Functional Theory

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

Study on Uv, Ir and Nmr Spectra of Double Hydrogen-Bonded Complexes

Theoretical Study of Binding Interactions and Vibrational Raman Spectra of Water in Hydrogen-Bonded Anionic Complexes: (H2o)n- (n = 2 and 3), H2O···X-(X = F, Cl, Br, and I), and H2O···M- (M = Cu, Ag, and Au)

Systematic Theoretical Study Of Structures And Bondings Of The Charge-Transfer Complexes Of Ammonia With Hx, Xy, And X(2) (X And Y Are Halogens)

X-ray Absorption Spectra for Aqueous Ammonia and Ammonium: Quantum Mechanical versus Molecular Mechanical Embedding Schemes

Structures, spectroscopy, binding and clustering energies of the hydrated copper dication clusters