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.

All-atom molecular dynamic simulations and relative NMR spectra study of weak C-H⋯⋯⋯O contacts in amide-water systems

Different Weak C−H···O Contacts in N-Methylacetamide-Water System: Molecular Dynamics Simulations and NMR Experimental Study

All-atom Simulation and Excess Properties Study on Intermolecular Interactions of Amide–water System

New Insights into a Hydrogen Bond: Hyper-Raman Spectroscopy of DMSO-water Solution.

Structures and Hydrogen Bonding Analysis of N,N-Dimethylformamide and N,N-Dimethylformamide−Water Mixtures by Molecular Dynamics Simulations

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

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

1H NMR studies of DMF-H2O association system

An All-Atom Simulation Study on Intermolecular Interaction of DMSO–water System

H-1 NMR studies of DMF-H2O association system

Molecular Dynamics Simulation for Formamide Aqueous Solution

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)

The Methyl C-H Blueshift in N,N-dimethylformamide-water Mixtures Probed by Two-Dimensional Fourier-transform Infrared Spectroscopy.

Comparative Study of Hydrogen Bonding Interactions Between N-methylacetamide and Methyl Acetate/Ethyl Formate

Automated Fragmentation QM/MM Calculation of Amide Proton Chemical Shifts in Proteins with Explicit Solvent Model.

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

<sup>2</sup>H NMR Studies of Simple Organic Groups Covalently Attached to Ordered Mesoporous Silica

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

Molecular Dynamics Simulation for Dmf Aqueous Solution

Sum-frequency Vibrational Spectroscopic Studies of Water/vapor Interfaces

The intermolecular interactions of methanol diluted in protic and aprotic solvent probed by polarized Raman spectroscopy and HNMR