From ab initio quantum chemistry to molecular dynamics: The delicate case of hydrogen bonding in ammonia

Abstract

The ammonia dimer (NH3)2 has been investigated using high--level ab initio quantum chemistry methods and density functional theory (DFT). The structure and energetics of important isomers is obtained to unprecedented accuracy without resorting to experiment. The global minimum of eclipsed C_s symmetry is characterized by a significantly bent hydrogen bond which deviates from linearity by about 20 degrees. In addition, the so-called cyclic C_{2h} structure is extremely close in energy on an overall flat potential energy surface. It is demonstrated that none of the currently available (GGA, meta--GGA, and hybrid) density functionals satisfactorily describe the structure and relative energies of this nonlinear hydrogen bond. We present a novel density functional, HCTH/407+, designed to describe this sort of hydrogen bond quantitatively on the level of the dimer, contrary to e.g. the widely used BLYP functional. This improved functional is employed in Car-Parrinello ab initio molecular dynamics simulations of liquid ammonia to judge its performance in describing the associated liquid. Both the HCTH/407+ and BLYP functionals describe the properties of the liquid well as judged by analysis of radial distribution functions, hydrogen bonding structure and dynamics, translational diffusion, and orientational relaxation processes. It is demonstrated that the solvation shell of the ammonia molecule in the liquid phase is dominated by steric packing effects and not so much by directional hydrogen bonding interactions. In addition, the propensity of ammonia molecules to form bifurcated and multifurcated hydrogen bonds in the liquid phase is found to be negligibly small.

Journal of Chemical Physics, in press (305335JCP)

References in corpus (5)

• Separable Dual Space Gaussian Pseudo-potentials
• Towards standard methods for benchmark quality ab initio thermochemistry --- W1 and W2 theory
• Assessment of W1 and W2 theories for the computation of electron affinities, ionization potentials, heats of formation, and proton affinities
• Hydrogen Bond Dynamics Near A Micellar Surface: Origin of the Universal Slow Relaxation at Complex Aqueous Interfaces
• Fast and slow dynamics of hydrogen bonds in liquid water

Cited by in corpus (8)

• Development of Novel Density Functionals for Thermochemical Kinetics
• Coupled cluster benchmarks of water monomers and dimers extracted from DFT liquid water: the importance of monomer deformations
• Ab-initio cluster approach for high-harmonic generation in liquids
• Molecular reorientation in hydrogen-bonding liquids: through algebraic $\sim t^{-3/2}$ relaxation toward exponential decay
• Benzene radical anion in the context of the Birch reduction: when solvation is the key
• Ionization of Ammonia Nanoices With Adsorbed Methanol Molecules
• Machine learning for accuracy in density functional approximations
• Enhanced sampling of robust molecular datasets with uncertainty-based collective variables