Protonic Defects in Hydrogen Bonded Liquids: Structure and Dynamics in Ammonia and Comparison with Water

Proton Transport Scenarios in Sulfuric Acid Explored via Ab Initio Molecular Dynamics Simulations

Sulfuric acid (H2SO4), a highly reactive reagent containing intrinsic protonic charge carriers, has been studied via ab initio molecular dynamics simulations. Specifically, we explore the solvation shell structure of the protonic defects, H1SO4- and H3SO4+, as well as the underlying proton transport mechanisms in both the neat and hydrated H2SO4 solutions. Our findings reveal a significant contraction of the dynamic hydrogen-bonded network around the protonic defects, which resembles features seen in water. The simulations provide estimates of the structural relaxation time scales for proton release from both the covalent O-H bonds (∼23 ps) and the hydrogen bonds (∼0.4 ps). In contrast to water, our analysis of the proton transfer scenarios in sulfuric acid reveals correlated events mediated by the formation of longer (up to four) hydrogen-bonded Grotthuss chains.

Facile ammonium oxidation to nitrogen gas in acid wastewater by in situ photogenerated chlorine radicals

The treatment of low-concentration ammonium (e.g., <50 mg L^-1) in highly acidic wastewaters through traditional biological nitrification, physical separation, or chemical stripping remains a huge challenge. Herein, we report that photocatalytic ammonium oxidation using bismuth oxychloride (BiOCl) can successfully occur in Cl^--laden solutions within a pH range of 1.0-6.0. All reactions follow pseudo-zero-order kinetics (with rate constants of 0.27-0.32 mg L^-1 min^-1 at pH 2.0-6.0 and 0.14 mg L^-1 min^-1 at pH 1.0), indicating the saturation of reactive species by the reactants. The interlayer is self-oxidized by the valence band holes (VB h⁺), resulting in the formation of Cl^• and subsequently HClO, which is excited upon UV irradiation to provoke consecutive photoreactions for chlorine radical generation. Compared to the free chlorine, HO^•, Cl^•, and Cl₂^•-, the ClO^• produced using the UV/BiOCl system plays a predominant role in oxidizing ammonium under acidic conditions. BiOCl exhibits good stability because of the compensation of Cl^- from solution and maintains high activity under different conditions (i.e., different cations and co-existing anions, temperatures, and initial substrate concentrations). The successful removal of ammonium from real wastewater using the UV/BiOCl system suggests that this is a promising method for treating diluted ammonium under highly acidic conditions.

Metallization of Shock-Compressed Liquid Ammonia

Ammonia is predicted to be one of the major components in the depths of the ice giant planets Uranus and Neptune. Their dynamics, evolution, and interior structure are insufficiently understood and models rely imperatively on data for equation of state and transport properties. Despite its great significance, the experimentally accessed region of the ammonia phase diagram today is still very limited in pressure and temperature. Here we push the probed regime to unprecedented conditions, up to ∼350  GPa and ∼40 000  K. Along the Hugoniot, the temperature measured as a function of pressure shows a subtle change in slope at ∼7000  K and ∼90  GPa, in agreement with ab initio simulations we have performed. This feature coincides with the gradual transition from a molecular liquid to a plasma state. Additionally, we performed reflectivity measurements, providing the first experimental evidence of electronic conduction in high-pressure ammonia. Shock reflectance continuously rises with pressure above 50 GPa and reaches saturation values above 120 GPa. Corresponding electrical conductivity values are up to 1 order of magnitude higher than in water in the 100 GPa regime, with possible significant contributions of the predicted ammonia-rich layers to the generation of magnetic dynamos in ice giant interiors.

Aqueous Solvation of Ammonia and Ammonium: Probing Hydrogen Bond Motifs with FT-IR and Soft X-ray Spectroscopy

In a multifaceted investigation combining local soft X-ray and vibrational spectroscopic probes with ab initio molecular dynamics simulations, hydrogen-bonding interactions of two key principal amine compounds in aqueous solution, ammonia (NH₃) and ammonium ion (NH₄⁺), are quantitatively assessed in terms of electronic structure, solvation structure, and dynamics. From the X-ray measurements and complementary determination of the IR-active hydrogen stretching and bending modes of NH₃ and NH₄⁺ in aqueous solution, the picture emerges of a comparatively strongly hydrogen-bonded NH₄⁺ ion via N-H donating interactions, whereas NH₃ has a strongly accepting hydrogen bond with one water molecule at the nitrogen lone pair but only weakly N-H donating hydrogen bonds. In contrast to the case of hydrogen bonding among solvent water molecules, we find that energy mismatch between occupied orbitals of both the solutes NH₃ and NH₄⁺ and the surrounding water prevents strong mixing between orbitals upon hydrogen bonding and, thus, inhibits substantial charge transfer between solute and solvent. A close inspection of the calculated unoccupied molecular orbitals, in conjunction with experimentally measured N K-edge absorption spectra, reveals the different nature of the electronic structural effects of these two key principal amine compounds imposed by hydrogen bonding to water, where a pH-dependent excitation energy appears to be an intrinsic property. These results provide a benchmark for hydrogen bonding of other nitrogen-containing acids and bases.

Ab initio molecular dynamics

In this chapter, an introduction to ab initio molecular dynamics (AIMD) has been given. Many of the basic concepts, like the Hellman-Feynman forces, the difference between the Car-Parrinello molecular dynamics and AIMD, have been explained. Also a very versatile AIMD code, the CP2K, has been introduced. On the application, the emphasis was on the aqueous systems and chemical reactions. The biochemical applications have not been discussed in depth.

Water-assisted proton transfer in ferredoxin I

The role of water molecules in assisting proton transfer (PT) is investigated for the proton-pumping protein ferredoxin I (FdI) from Azotobacter vinelandii. It was shown previously that individual water molecules can stabilize between Asp(15) and the buried [3Fe-4S](0) cluster and thus can potentially act as a proton relay in transferring H(+) from the protein to the μ(2) sulfur atom. Here, we generalize molecular mechanics with proton transfer to studying proton transfer reactions in the condensed phase. Both umbrella sampling simulations and electronic structure calculations suggest that the PT Asp(15)-COOH + H(2)O + [3Fe-4S](0) → Asp(15)-COO(-) + H(2)O + [3Fe-4S](0) H(+) is concerted, and no stable intermediate hydronium ion (H(3)O(+)) is expected. The free energy difference of 11.7 kcal/mol for the forward reaction is in good agreement with the experimental value (13.3 kcal/mol). For the reverse reaction (Asp(15)-COO(-) + H(2)O + [3Fe-4S](0)H(+) → Asp(15)-COOH + H(2)O + [3Fe-4S](0)), a larger barrier than for the forward reaction is correctly predicted, but it is quantitatively overestimated (23.1 kcal/mol from simulations versus 14.1 from experiment). Possible reasons for this discrepancy are discussed. Compared with the water-assisted process (ΔE ≈ 10 kcal/mol), water-unassisted proton transfer yields a considerably higher barrier of ΔE ≈ 35 kcal/mol.

A first principles molecular dynamics study of excess electron and lithium atom solvation in water-ammonia mixed clusters: structural, spectral, and dynamical behaviors of [(H2O)5NH3]- and Li(H2O)5NH3 at finite temperature

First principles molecular dynamics simulations are carried out to investigate the solvation of an excess electron and a lithium atom in mixed water-ammonia cluster (H(2)O)(5)NH(3) at a finite temperature of 150 K. Both [(H(2)O)(5)NH(3)](-) and Li(H(2)O)(5)NH(3) clusters are seen to display substantial hydrogen bond dynamics due to thermal motion leading to many different isomeric structures. Also, the structures of these two clusters are found to be very different from each other and also very different from the corresponding neutral cluster without any excess electron or the metal atom. Spontaneous ionization of Li atom occurs in the case of Li(H(2)O)(5)NH(3). The spatial distribution of the singly occupied molecular orbital shows where and how the excess (or free) electron is primarily localized in these clusters. The populations of single acceptor (A), double acceptor (AA), and free (NIL) type water and ammonia molecules are found to be significantly high. The dangling hydrogens of these type of water or ammonia molecules are found to primarily capture the free electron. It is also found that the free electron binding motifs evolve with time due to thermal fluctuations and the vertical detachment energy of [(H(2)O)(5)NH(3)](-) and vertical ionization energy of Li(H(2)O)(5)NH(3) also change with time along the simulation trajectories. Assignments of the observed peaks in the vibrational power spectra are done and we found a one to one correlation between the time-averaged populations of water and ammonia molecules at different H-bonding sites with the various peaks of power spectra. The frequency-time correlation functions of OH stretch vibrational frequencies of these clusters are also calculated and their decay profiles are analyzed in terms of the dynamics of hydrogen bonded and dangling OH modes. It is found that the hydrogen bond lifetimes in these clusters are almost five to six times longer than that of pure liquid water at room temperature.

Excess electron and lithium atom solvation in water clusters at finite temperature: an ab initio molecular dynamics study of the structural, spectral, and dynamical behavior of (H2O)6- and Li(H2O)6

The roles of hydrogen bonds in the solvation of an excess electron and a lithium atom in water hexamer cluster at 150 K have been studied by means of ab initio molecular dynamics simulations. It is found that the hydrogen bonded structures of (H(2)O)(6)(-) and Li(H(2)O)(6) clusters are very different from each other and they dynamically evolve from one conformer to other along their simulation trajectories. The populations of the single acceptor, double acceptor, and free type water molecules are found to be significantly high unlike that in pure water clusters. Free hydrogens of these type of water molecules primarily capture the unbound electron density in these clusters. It is found that the binding motifs of the free electron evolve with time and the vertical detachment energy of (H(2)O)(6)(-) and vertical ionization energy of Li(H(2)O)(6) also change with time. Assignments of the observed peaks in vibrational power spectra are done, and we found direct correlations between the time-averaged population of water molecules in different hydrogen bonding states and the spectral features. The dynamical aspects of these clusters have also been studied through calculations of time correlations of instantaneous stretch frequencies of OH modes which are obtained from the simulation trajectories through a time series analysis.

A multistate empirical valence bond model for solvation and transport simulations of OH- in aqueous solutions

We describe a new multistate empirical valence bond (MS-EVB) model of OH(-) in aqueous solutions. This model is based on the recently proposed "charged ring" parameterization for the intermolecular interaction of hydroxyl ion with water [Ufimtsev, et al., Chem. Phys. Lett., 2007, 442, 128] and is suitable for classical molecular simulations of OH(-) solvation and transport. The model reproduces the hydration structure of OH(-)(aq) in good agreement with experimental data and the results of ab initio molecular dynamics simulations. It also accurately captures the major structural, energetic, and dynamic aspects of the proton transfer processes involving OH(-) (aq). The model predicts an approximately two-fold increase of the OH(-) mobility due to proton exchange reactions.

Pseudo-Jahn-Teller origin of the low barrier hydrogen bond in N(2)H(7) (+)

The microscopic origin and quantum effects of the low barrier hydrogen bond (LBHB) in the proton-bound ammonia dimer cation N(2)H(7) (+) were studied by means of ab initio and density-functional theory (DFT) methods. These results were analyzed in the framework of vibronic theory and compared to those obtained for the Zundel cation H(5)O(2) (+). All geometry optimizations carried out using wavefunction-based methods [Hartree-Fock, second and fourth order Moller-Plesset theory (MP2 and MP4), and quadratic configuration interaction with singles and doubles excitations (QCISD)] lead to an asymmetrical H(3)N-H(+)cdots, three dots, centeredNH(3) conformation (C(3v) symmetry) with a small energy barrier (1.26 kcalmol in MP4 and QCISD calculations) between both equivalent minima. The value of this barrier is underestimated in DFT calculations particularly at the local density approximation level where geometry optimization leads to a symmetric H(3)Ncdots, three dots, centeredH(+)cdots, three dots, centeredNH(3) structure (D(3d) point group). The instability of the symmetric D(3d) structure is shown to originate from the pseudo-Jahn-Teller mixing of the electronic (1)A(1g) ground state with five low lying excited states of A(2u) symmetry through the asymmetric alpha(2u) vibrational mode. A molecular orbital study of the pseudo-Jahn-Teller coupling has allowed us to discuss the origin of the proton displacement and the LBHB formation in terms of the polarization of the NH(3) molecules and the transfer of electronic charge between the proton and the NH(3) units (rebonding). The parallel study of the H(5)O(2) (+) cation, which presents a symmetric single-well structure, allows us to analyze why these similar molecules behave differently with respect to proton transfer. From the vibronic analysis, a unified view of the Rudle-Pimentel three-center four-electron and charge transfer models of LBHBs is given. Finally, the large difference in the N-N distance in the D(3d) and C(3v) configurations of N(2)H(7) (+) indicates a large anharmonic coupling between alpha(2u)-alpha(1g) modes along the proton-transfer dynamics. This issue was explored by solving numerically the vibrational Schrodinger equation corresponding to the bidimensional E[Q(alpha(2u)),Q(alpha(1g))] energy surface calculated at the MP46-311++G(**) level of theory.

Reactive force fields for proton transfer dynamics

A force field-inspired method based on fitted, high-quality multidimensional potential energy surfaces to follow proton transfer (PT) reactions in molecular dynamics simulations is presented. In molecular mechanics with proton transfer (MMPT) a system is partitioned into a region where proton transfer takes place and the remaining degrees of freedom which are treated with a conventional force field. The implementation of the method and applications to specific chemically and biologically relevant scenarios are presented. MMPT is developed in view of two primary areas in mind: to follow the molecular dynamics of proton transfer in the condensed phase on realistic time scales and to adapt the shape (morphing) of the potential energy surface for specific applications. MMPT is applied to PT in protonated ammonia dimer, double proton transfer in 2-pyridone-2-hydroxypyridine, and the first step of PT from a protein side-chain towards a buried [3Fe4S] cluster in ferredoxin I. Specific findings of the work include the fundamental role of the N-N vibration as the gating mode for PT in NH4+...NH3 and the qualitative understanding of PT from the protein to a metastable active-site water molecule in Ferredoxin I.

Investigating the relationship between infrared spectra of shared protons in different chemical environments: a comparison of protonated diglyme and protonated water dimer

The energetics, dynamics, and infrared spectroscopy of the shared proton in different chemical environments is investigated using molecular dynamics simulations. A three-dimensional potential energy surface (PES) suitable for describing proton transfer between an acceptor and a donor oxygen atom is combined with an all-atom force field to carry out reactive molecular dynamics simulations. The construction of the fully dimensional PES is inspired from the established mixed quantum mechanics/molecular mechanics treatment of larger systems. The "morphing potential" method is used to transform the generic PES for proton transfer along an O...H+...O motif into a three-dimensional PES for proton transfer in protonated diglyme. Using molecular dynamics simulations at finite temperature, the gas phase infrared spectra are calculated for both species from the Fourier transform of the dipole moment autocorrelation function. For protonated diglyme the modes involving the H+ motion are strongly mixed with other degrees of freedom. At low temperature, the O...H+...O asymmetric stretching vibration is found at 870 cm-1, whereas for H5O2+ this band is at 724 cm-1. As expected, the vibrational bands of protonated diglyme show no temperature dependence whereas for H5O2+ at T = 100 K the proton transfer mode is found at 830 cm-1, in good agreement with 861 cm-1 from very recent molecular dynamics simulations.

Simulations of Liquid Ammonia Based on the Combined Quantum Mechanical/Molecular Mechanical (QM/MM) Approach

Two combined quantum mechanics/molecular mechanics (QM/MM) molecular dynamics simulations, namely, HF/MM and B3LYP/MM, have been performed to investigate the local structure and dynamics of liquid ammonia. The most interesting region, a sphere containing a central reference molecule and all its nearest surrounding molecules (first coordination shell), was treated by the Hartree-Fock (HF) and hybrid density functional B3LYP methods, whereas the rest of the system was described by the classical pair potentials. On the basis of both HF and B3LYP methods, it is observed that the hydrogen bonding in this peculiar liquid is weak. The structure and dynamics of this liquid are suggested to be determined by the steric packing effects, rather than by the directional hydrogen bonding interactions. Compared to previous empirical as well as Car-Parrinello (CP) molecular dynamics studies, our QM/MM simulations provide detailed information that is in better agreement with experimental data.

On the Quantum Nature of an Excess Proton in Liquid Hydrogen Fluoride

Liquid hydrogen fluoride consists of chains of hydrogen-bonded molecules. The nature of an excess proton in liquid HF, which is of interest not only for its own sake, but also in relation to super-acid chemistry and to its behavior in water, has been studied using computer simulations. The methodology employed is the density-functional-theory-based path-integral Car-Parrinello ab initio molecular dynamics. The excess proton, which formally exists as a H2F+ or a H2F2+ defect in an HF chain, is found to strongly perturb the chain to which it is attached. Moreover, due to large zero-point energy, the charge defect is largely delocalized over several HF molecules.

Ab initiomolecular-dynamics study of liquid formamide

Properties of neat liquid formamide (HCONH2) have been studied by the combination of gradient-corrected density-functional theory, norm-conserving pseudopotentials, and the adaptive finite-element method. The structural and dynamical quantities have been calculated through molecular dynamics simulations under the Born-Oppenheimer approximation. Satisfactory agreement with experimental data was obtained for both intramolecular and intermolecular properties. Our results are also compared with those of the empirical potential functions to clarify their accuracies.

On the influence of semirigid environments on proton transfer along molecular chains

The dynamics of proton transfer along ammonia chains (chemical composition N(x)H(+)(3x+1), x=2, 4, and 6) in a constraining environment is investigated by ab initio molecular dynamics simulations. A carbon nanotube of defined length and diameter is used as an idealized constraining environment such that the ammonia chain is forced to maintain its quasilinear geometry. It is found that, although the energetics of proton transport shows considerable energetic barriers, proton translocation along the wire is possible at finite temperature for all chain lengths studied. The proton transport involves rotational reorientation of the proton-carrying ammonia molecule. High level ab initio calculations (MP2/aug-cc-pVTZ) yield barriers for internal rotation of 9.1 kcal/mol for NH(4) (+)-NH(3) and 11.7 kcal/mol for OH(3) (+)-OH(2), respectively. The infrared spectrum calculated from the dipole-dipole autocorrelation function shows distinct spectral features in the regions (2000-3000 cm(-1)) where the NHN proton transfer mode is expected to absorb. Assigning moderate opposite total charges between 0.002 and 0.2e to the carbon atoms at the end caps of the nanotube leads to a considerable speedup of the proton transfer.

The nature and transport mechanism of hydrated hydroxide ions in aqueous solution

Compared to other ions, protons (H(+)) and hydroxide ions (OH(-)) exhibit anomalously high mobilities in aqueous solutions. On a qualitative level, this behaviour has long been explained by 'structural diffusion' the continuous interconversion between hydration complexes driven by fluctuations in the solvation shell of the hydrated ions. Detailed investigations have led to a clear understanding of the proton transport mechanism at the molecular level. In contrast, hydroxide ion mobility in basic solutions has received far less attention, even though bases and base catalysis play important roles in many organic and biochemical reactions and in the chemical industry. The reason for this may be attributed to the century-old notion that a hydrated OH(-) can be regarded as a water molecule missing a proton, and that the transport mechanism of such a 'proton hole' can be inferred from that of an excess proton by simply reversing hydrogen bond polarities. However, recent studies have identified OH(-) hydration complexes that bear little structural similarity to proton hydration complexes. Here we report the solution structures and transport mechanisms of hydrated hydroxide, which we obtained from first-principles computer simulations that explicitly treat quantum and thermal fluctuations of all nuclei. We find that the transport mechanism, which differs significantly from the proton hole picture, involves an interplay between the previously identified hydration complexes and is strongly influenced by nuclear quantum effects.