Hydrogen Bond Network Of Water Research Articles

Strategy for Modeling the Infrared Spectra of Ion-Containing Water Drops.

Hydrated ions are ubiquitous in environmental and biological media. Understanding the perturbation exerted by an ion on the water hydrogen bond network is possible in the nanodrop regime by recording vibrational spectra in the O-H bond stretching region. This has been achieved experimentally in recent years by forming gaseous ions containing tens to hundreds of water molecules and recording their infrared photodissociation spectra. In this paper, we demonstrate the capabilities of a modeling strategy based on an extension of the AMOEBA polarizable force field to implement water atomic charge fluctuations along with those of intramolecular structure along the dynamics. This supplementary flexibility of nonbonded interactions improves the description of the hydrogen bond network and, therefore, the spectroscopic response. Finite temperature IR spectra are obtained from molecular dynamics simulations by computing the Fourier transform of the dipole moment autocorrelation function. Simulations of 1-2 ns are required for extensive sampling in order to reproduce the experimental spectra. Furthermore, bands are assigned with the driven molecular dynamics approach. This method package is shown to compare successfully with experimental spectra for 11 ions in water drops containing 36-100 water molecules. In particular, band frequency shifts of the free O-H stretching modes at the cluster surface are well reproduced as a function of both ion charge and drop size.

Shrinking of Rapidly Evaporating Water Microdroplets Reveals their Extreme Supercooling.

The fast evaporative cooling of micrometer-sized water droplets in a vacuum offers the appealing possibility to investigate supercooled water-below the melting point but still a liquid-at temperatures far beyond the state of the art. However, it is challenging to obtain a reliable value of the droplet temperature under such extreme experimental conditions. Here, the observation of morphology-dependent resonances in the Raman scattering from a train of perfectly uniform water droplets allows us to measure the variation in droplet size resulting from evaporative mass losses with an absolute precision of better than 0.2%. This finding proves crucial to an unambiguous determination of the droplet temperature. In particular, we find that a fraction of water droplets with an initial diameter of 6379±12 nm remain liquid down to 230.6±0.6 K. Our results question temperature estimates reported recently for larger supercooled water droplets and provide valuable information on the hydrogen-bond network in liquid water in the hard-to-access deeply supercooled regime.

Persistent homology analysis of ion aggregations and hydrogen-bonding networks.

Despite the great advancement of experimental tools and theoretical models, a quantitative characterization of the microscopic structures of ion aggregates and their associated water hydrogen-bonding networks still remains a challenging problem. In this paper, a newly-invented mathematical method called persistent homology is introduced, for the first time, to quantitatively analyze the intrinsic topological properties of ion aggregation systems and hydrogen-bonding networks. The two most distinguishable properties of persistent homology analysis of assembly systems are as follows. First, it does not require a predefined bond length to construct the ion or hydrogen-bonding network. Persistent homology results are determined by the morphological structure of the data only. Second, it can directly measure the size of circles or holes in ion aggregates and hydrogen-bonding networks. To validate our model, we consider two well-studied systems, i.e., NaCl and KSCN solutions, generated from molecular dynamics simulations. They are believed to represent two morphological types of aggregation, i.e., local clusters and extended ion networks. It has been found that the two aggregation types have distinguishable topological features and can be characterized by our topological model very well. Further, we construct two types of networks, i.e., O-networks and H2O-networks, for analyzing the topological properties of hydrogen-bonding networks. It is found that for both models, KSCN systems demonstrate much more dramatic variations in their local circle structures with a concentration increase. A consistent increase of large-sized local circle structures is observed and the sizes of these circles become more and more diverse. In contrast, NaCl systems show no obvious increase of large-sized circles. Instead a consistent decline of the average size of the circle structures is observed and the sizes of these circles become more and more uniform with a concentration increase. As far as we know, these unique intrinsic topological features in ion aggregation systems have never been pointed out before. More importantly, our models can be directly used to quantitatively analyze the intrinsic topological invariants, including circles, loops, holes, and cavities, of any network-like structures, such as nanomaterials, colloidal systems, biomolecular assemblies, among others. These topological invariants cannot be described by traditional graph and network models.

Development of non-bonded interaction parameters between graphene and water using particle swarm optimization.

New Lennard-Jones parameters have been developed to describe the interactions between atomistic model of graphene, represented by REBO potential, and five commonly used all-atom water models, namely SPC, SPC/E, SPC/Fw, SPC/Fd, and TIP3P/Fs by employing particle swarm optimization (PSO) method. These new parameters were optimized to reproduce the macroscopic contact angle of water on a graphene sheet. The calculated line tension was in the order of 10-11 J/m for the droplets of all water models. Our molecular dynamics simulations indicate the preferential orientation of water molecules near graphene-water interface with one OH bond pointing toward the graphene surface. Detailed analysis of simulation trajectories reveals the presence of water molecules with ≤∼1, ∼2, and ∼4 hydrogen bonds at the surface of air-water interface, graphene-water interface, and bulk region of the water droplet, respectively. Presence of water molecules with ≤∼1 and ∼2 hydrogen bonds suggest the existence of water clusters of different sizes at these interfaces. The trends observed in the libration, bending, and stretching bands of the vibrational spectra are closely associated with these structural features of water. The inhomogeneity in hydrogen bond network of water at the air-water and graphene-water interface is manifested by broadening of the peaks in the libration band for water present at these interfaces. The stretching band for the molecules in water droplet shows a blue shift as compared to the pure bulk water, which conjecture the presence of weaker hydrogen bond network in a droplet. © 2017 Wiley Periodicals, Inc.

Direct Anionic Effect on Water Structure and Indirect Anionic Effect on Peptide Backbone Hydration State Revealed by Thin-Layer Infrared Spectroscopy.

In this work, the anionic effect on water structure and on the peptide backbone and water interaction was investigated directly in aqueous solution using thin-layer transmission infrared spectroscopy. The chaotropic anions were found to weaken the water hydrogen-bonding strength and red shift the HOH bending frequency, while the kosmotropic anions were found to strengthen the water hydrogen-bonding network and blue shift the HOH bending frequency. The kosmotropes, especially F-, blue shift the vibrational frequencies of both amide II and amide III bands of N-methylacetamide (NMA), indicating NMA is in the "salting-in" state; while the chaotropes (Cl-, NO3-, Br-, I-, and SCN-) red shift the frequencies of the two normal modes, indicating NMA is in the "salting-out" state. Furthermore, the changes of the vibrational frequencies of the HOH bending, amide II and III bands were found to generally follow the Hofmeister anionic series. Our results suggest that hydrated anion influences the peptide backbone mainly through the N-H group, but a weak and indirect effect through the amide C═O group also contributes. Thus, these amide modes can be used as vibrational measures of anionic influences on peptide backbone's hydration state. Our work also suggests that deuteration of the amide unit decreases the sensitivity of the amide II and III vibrational modes in this regard.

An "Iceberg" Coating Preserves Bulk Hydration Dynamics in Aqueous PEG Solutions.

Ultrafast picosecond time scale two-dimensional infrared (2D-IR) spectroscopy of a new water-soluble transition metal complex acting as a vibrational probe shows that over a range of concentration and poly(ethylene glycol) (PEG) molecular mass (2000, 8000, and 20000 Da) the time scale of the sensed hydration dynamics differs negligibly from bulk water (D2O). PEG is well-known to establish a highly stable hydration shell because the spacing between adjacent ethereal oxygens nearly matches water's hydrogen-bonding network. Although these first-shell water molecules are likely significantly retarded, they present an interface to subsequent hydration shells and thus diminish the largely entropic perturbation to water's orientational dynamics. In addition to the longer PEGs, a series of concentration-dependent 2D-IR measurements using aqueous PEG-400 show a pronounced hydration slowdown in the vicinity of the critical overlap concentration (c*). Comparison between these dynamical results and previously reported steady-state infrared spectroscopy of aqueous PEG-1000 solutions reveals a strikingly identical dependence on number of water molecules per ethylene oxide monomer, scaled according to the critical overlap concentration.

The Bend+Libration Combination Band Is an Intrinsic, Collective, and Strongly Solute-Dependent Reporter on the Hydrogen Bonding Network of Liquid Water.

Water is an extensively self-associated liquid due to its extensive hydrogen bond (H-bond) forming ability. The resulting H-bonded network fluid exhibits nearly continuous absorption of light from the terahertz to the near-IR region. The relatively weak bend+libration water combination band (centered at 2130 cm-1) has been largely overlooked as a reporter of liquid water's structure and dynamics despite its location in a convenient region of the IR for spectroscopic study. The intermolecular nature of the combination band leads to a unique absorption signal that reports collectively on the rigidity of the H-bonding network in the presence of many different solutes. This study reports comprehensively how the combination band acts as an intrinsic and collective probe in various chemically and biologically relevant solutions, including salts of varying character, denaturants, osmolytes, crowders, and surfactants that form reverse micelles and micelles. While we remark on changes in the line width and intensity of this combination band, we mainly focus on the frequency and how the frequency reports on the collective H-bonding network of liquid water. We also comment on the "association band" moniker often applied to this band and how to evaluate discrete features in this spectral region that sometimes appear in the IR spectra of specific kinds of aqueous samples of organic solutes, especially those with very high solute concentrations, with the conclusion that most of these discrete spectral features come exclusively from the solutes and do not report on the water. Contrasts are drawn throughout this work between the collective and delocalized reporting ability of the combination band and the response of more site-specific vibrations like the much-investigated OD stretch of HDO in H2O: the combination band is a unique reporter of H-bonding structure and dynamics and fundamentally different than any local mode probe. Since this band appears as the spectroscopic "background" for many local-mode reporter groups, we note the possibility of observing both local and collective solvent dynamics at the same time in this spectral region.

Ion aggregation in high salt solutions. VII. The effect of cations on the structures of ion aggregates and water hydrogen-bonding network.

Ions in high salt solutions have a strong propensity to form polydisperse ion aggregates with broad size and shape distributions. In a series of previous comparative investigations using femtosecond IR pump-probe spectroscopy, molecular dynamics simulation, and graph theoretical analysis, we have shown that there exists a morphological difference in the structures of ion aggregates formed in various salt solutions. As salt concentration increases, the ions in high salt solutions form either cluster-like structures excluding water molecules or network-like structures entwined with water hydrogen-bonding networks. Interestingly, such morphological characteristics of the ion aggregates have been found to be in correlation with the solubility limits of salts. An important question that still remains unexplored is why certain salts with different cations have notably different solubility limits in water. Here, carrying out a series of molecular dynamics simulations of aqueous salt solutions and analyzing the distributions and connectivity patterns of ion aggregates with a spectral graph analysis method, we establish the relationship between the salt solubility and the ion aggregate morphology with a special emphasis on the cationic effects on water structures and ion aggregation. We anticipate that the understanding of large scale ion aggregate structures revealed in this study will be critical for elucidating the specific ion effects on the solubility and conformational stability of co-solute molecules such as proteins in water.

Taurine as a water structure breaker and protein stabilizer

The enhancing effect on the water structure has been confirmed for most of the osmolytes exhibiting both stabilizing and destabilizing properties in regard to proteins. The presented work concerns osmolytes, which should be classified as “structure breaking” solutes: taurine and N,N,N-trimethyltaurine (TMT). Here, we combine FTIR spectroscopy, DSC calorimetry and DFT calculations to gain an insight into the interactions between osmolytes and two proteins: lysozyme and ubiquitin. Despite high structural similarity, both osmolytes exert different influence on protein stability: taurine is a stabilizer, TMT is a denaturant. We show also that taurine amino group interacts directly with the side chains of proteins, whereas TMT does not interact with proteins at all. Although two solutes weaken on average the structure of the surrounding water, their hydration spheres are different. Taurine is surrounded by two populations of water molecules: bonded with weak H-bonds around sulfonate group, and strongly bonded around amino group. The strong hydrogen-bonded network of water molecules around the amino group of taurine further improves properties of enhanced protein hydration sphere and stabilizes the native protein form. Direct interactions of this group with surface side chains provide a proper orientation of taurine and prevents the {text{SO}}_{3}^{ - } group from negative influence. The weakened {text{SO}}_{3}^{ - } hydration sphere of TMT breaks up the hydrogen-bonded network of water around the protein and destabilizes it. However, TMT at low concentration stabilize both proteins to a small extent. This effect can be attributed to an actual osmophobic effect which is overcome if the concentration increases.

Proton Diffusion through Bilayer Pores.

The transport of protons through channels in complex environments is important in biology and materials science. In this work, we use multistate empirical valence bond simulations to study proton transport within a well-defined bilayer pore in a lamellar Lβ phase lyotropic liquid crystal (LLC). The LLC is formed from the self-assembly of dicarboxylate gemini surfactants in water, and a bilayer-spanning pore of radius of approximately 3-5 Å results from the uneven partitioning of surfactants between the two leaflets of the lamella. Local proton diffusion within the pore is significantly faster than diffusion at the bilayer surface, which is due to the greater hydrophobicity of the surfactant/water interface within the pore. Proton diffusion proceeds by surface transport along exposed hydrophobic pockets at the surfactant/water interface and depends on the continuity of hydronium-water hydrogen bond networks. At the bilayer surface, there is a reduced fraction of the "Zundel" intermediates that are central to the Grotthuss transport mechanism, whereas the fraction of these species within the bilayer pore is similar to that in bulk water. Our results demonstrate that the chemical nature of the confining interface, in addition to confinement length scale, is an important determiner of local proton transport in nanoconfined aqueous environments.

Multivariate Analysis of Mixed Lipid Aggregate Phase Transitions Monitored Using Raman Spectroscopy.

The phase behavior of aqueous 1,2-dimyristoyl-sn-glycero-3-phosphorylcholine (DMPC)/1,2-dihexanoyl-sn-glycero-3-phosphocholine (DHPC) mixtures between 8.0 ℃ and 41.0 ℃ were monitored using Raman spectroscopy. Temperature-dependent Raman matrices were assembled from series of spectra and subjected to multivariate analysis. The consensus of pseudo-rank estimation results is that seven to eight components account for the temperature-dependent changes observed in the spectra. The spectra and temperature response profiles of the mixture components were resolved by applying a variant of the non-negative matrix factorization (NMF) algorithm described by Lee and Seung (1999). The rotational ambiguity of the data matrix was reduced by augmenting the original temperature-dependent spectral matrix with its cumulative counterpart, i.e., the matrix formed by successive integration of the spectra across the temperature index (columns). Successive rounds of constrained NMF were used to isolate component spectra from a significant fluorescence background. Five major components exhibiting varying degrees of gel and liquid crystalline lipid character were resolved. Hydrogen-bonded water networks exhibiting varying degrees of organization are associated with the lipid components. Spectral parameters were computed to compare the chain conformation, packing, and hydration indicated by the resolved spectra. Based on spectral features and relative amounts of the components observed, four components reflect long chain lipid response. The fifth component could reflect the response of the short chain lipid, DHPC, but there were no definitive spectral features confirming this assignment. A minor component of uncertain assignment that exhibits a striking response to the DMPC pre-transition and chain melting transition also was recovered. While none of the spectra resolved exhibit features unequivocally attributable to a specific aggregate morphology or step in the gelation process, the results are consistent with the evolution of mixed phase bicelles (nanodisks) and small amounts of worm-like DMPC/DHPC aggregates, and perhaps DHPC micelles, at low temperature to suspensions of branched and entangled worm-like aggregates above the DMPC gel phase transition and perforated multi-lamellar aggregates at high temperature.

Signatures of Solvation Thermodynamics in Spectra of Intermolecular Vibrations.

This study explores the thermodynamic and vibrational properties of water in the three-dimensional environment of solvated ions and small molecules using molecular simulations. The spectrum of intermolecular vibrations in liquid solvents provides detailed information on the shape of the local potential energy surface, which in turn determines local thermodynamic properties such as the entropy. Here, we extract this information using a spatially resolved extension of the two-phase thermodynamics method to estimate hydration water entropies based on the local vibrational density of states (3D-2PT). Combined with an analysis of solute–water and water–water interaction energies, this allows us to resolve local contributions to the solvation enthalpy, entropy, and free energy. We use this approach to study effects of ions on their surrounding water hydrogen bond network, its spectrum of intermolecular vibrations, and resulting thermodynamic properties. In the three-dimensional environment of polar and nonpolar functional groups of molecular solutes, we identify distinct hydration water species and classify them by their characteristic vibrational density of states and molecular entropies. In each case, we are able to assign variations in local hydration water entropies to specific changes in the spectrum of intermolecular vibrations. This provides an important link for the thermodynamic interpretation of vibrational spectra that are accessible to far-infrared absorption and Raman spectroscopy experiments. Our analysis provides unique microscopic details regarding the hydration of hydrophobic and hydrophilic functional groups, which enable us to identify interactions and molecular degrees of freedom that determine relevant contributions to the solvation entropy and consequently the free energy.

Configurational Disorder of Water Hydrogen-Bond Network at the Protein Dynamical Transition.

We introduce a novel strategy to quantify the disorder of extended water-water hydrogen-bond (HB) networks sampled in particle-based computer simulations. The method relies on the conformational clustering of the HB connectivity states. We successfully applied it to unveil the fine relationship among the protein dynamical transition in hydrated powder, which marks the activation of protein flexibility at Td ≈ 240 K, and the sudden increase in the configurational disorder of the water HB network enveloping the proteins. Our finding links, in the spirit of the Adam-Gibbs relationship, the diffusivity of protein atoms, as quantified by the hydrogen mean-square displacements, and the thermodynamic solvent configurational entropy.

Perturbation of hydrogen bonding in hydrated pyrrole-2-carboxaldehyde complexes.

The interaction of external water molecules with hydrated pyrrole-2-carboxaldehyde PCL/(H2O) n complexes was investigated. The work was supported by both theoretical [DFT/TD-DFT methods using 6-311G++(d,p) basis set in the ground (S0) and excited (S1, S2, S3)states] and experimental [UV-Vis, FTIR and Raman] verification. The focus of the present work was on the weak intermolecular O-H⋯O, N-H⋯O-H hydrogen bonded interaction (IerHB) between PCL and external water molecules, and the influence of increasing the number of water molecules to form hydrated PCL/(H2O)n complexes. Effects were observed on different vibrational normal modes and on electronic transition levels. A hydrogen-bonded network of water induces a shift to higher energy in certain normal modes of PCL to form stable PCL/(H2O)n complexes by lowering the barrier energy. Potential energy distribution (PED) analysis indicates a significant charge transfer from PCL to water by creating a water bridge. Hydrogen bonding effects account for the substantial red shift and broadness in νNH, νCO vibrational modes. Water rearrangement turns out to be the main driving force for hydrated complex formation. Graphical abstract Stability of PCL/(H2O)4 hydarted complex.

Hydrogen Bond Networks and Hydrophobic Effects in the Amyloid β30-35 Chain in Water: A Molecular Dynamics Study.

We have studied the conformational landscape of the C-terminal fragment of the amyloid protein Aβ30-35 in water using well-tempered metadynamics simulations and found that it resembles an intrinsically disordered protein. The conformational fluctuations of the protein are facilitated by a collective reorganization of both protein and water hydrogen bond networks, combined with electrostatic interactions between termini as well as hydrophobic interactions of the side chains. The stabilization of hydrophobic interactions in one of the conformers involves a collective collapse of the side chains along with a squeeze-out of water sandwiched between them. The charged N- and C-termini play a critical role in stabilizing different types of protein conformations, including those involving contact-ion salt bridges as well as solvent-mediated interactions of the termini and the amide backbone. We have examined this by probing the distribution of directed water wires forming the hydrogen bond network enveloping the polypeptide. Water wires and their fluctuations form an integral part of structural signature of the protein conformation.

The Jones-Ray effect reinterpreted: Surface tension minima of low ionic strength electrolyte solutions are caused by electric field induced water-water correlations

The surface tension of electrolyte solutions exhibits a minimum at millimolar electrolyte concentrations and then rises with increasing concentration. This minimum, known as the Jones-Ray effect, has been hotly debated over the past ∼80years. If not considered as an artifact, it is typically ascribed to a phenomenological rare binding site for ions or ion pairs. Here, we propose an alternative underlying mechanism, namely that the hydrogen bond network of water responds to the collective electrostatic field of ions by increasing its orientational order, supported by recent surface tension measurements of NaCl solutions in H2O and D2O, and second harmonic scattering experiments in combination with ion resonant second harmonic reflection experiments. Recent thermodynamic and purely electrostatic treatments of the surface tension provide support for this interpretation. In addition, concerns related to possible artifacts influencing the measurements are quantified experimentally.

Water-anion hydrogen bonding dynamics: Ultrafast IR experiments and simulations.

Many of water's remarkable properties arise from its tendency to form an intricate and robust hydrogen bond network. Understanding the dynamics that govern this network is fundamental to elucidating the behavior of pure water and water in biological and physical systems. In ultrafast nonlinear infrared experiments, the accessible time scales are limited by water's rapid vibrational relaxation (1.8 ps for dilute HOD in H2O), precluding interrogation of slow hydrogen bond evolution in non-bulk systems. Here, hydrogen bonding dynamics in bulk D2O were studied from the perspective of the much longer lived (36.2 ps) CN stretch mode of selenocyanate (SeCN-) using polarization selective pump-probe (PSPP) experiments, two-dimensional infrared (2D IR) vibrational echo spectroscopy, and molecular dynamics simulations. The simulations make use of the empirical frequency mapping approach, applied to SeCN- for the first time. The PSPP experiments and simulations show that the orientational correlation function decays via fast (2.0 ps) restricted angular diffusion (wobbling-in-a-cone) and complete orientational diffusive randomization (4.5 ps). Spectral diffusion, quantified in terms of the frequency-frequency correlation function, occurs on two time scales. The initial 0.6 ps time scale is attributed to small length and angle fluctuations of the hydrogen bonds between water and SeCN-. The second 1.4 ps measured time scale, identical to that for HOD in bulk D2O, reports on the collective reorganization of the water hydrogen bond network around the anion. The experiments and simulations provide details of the anion-water hydrogen bonding and demonstrate that SeCN- is a reliable vibrational probe of the ultrafast spectroscopy of water.

Quantifying Hydrogen-Bond Populations in Dimethyl Sulfoxide/Water Mixtures.

Dimethyl sulfoxide (DMSO) disrupts the hydrogen-bond networks in water. The widespread use of DMSO as a cosolvent, along with its unusual attributes, have inspired numerous studies. Herein, infrared absorption spectroscopy of the S=O stretching mode combined with molecular dynamics and quantum chemistry models were used to directly quantify DMSO/water hydrogen-bond populations in binary mixtures. Singly H-bonded species are dominant at 10 mol %, due to strong DMSO-water interactions. We found an unexpected increase in non-hydrogen-bonded DMSO near the eutectic point (ca. 35 mol %) which also correlates with several abnormalities in the bulk solution properties. We find evidence for three distinct regimes: 1) strong DMSO-water interactions (<30 mol %), 2) ideal-solution-like (30-90 mol %), and 3) self-interaction, or aggregation, regime (>90 mol %). We propose a "step in" mechanism, which involves hydrogen bonding between water and the DMSO aggregate species.

Real-time monitoring of hydrophobic aggregation reveals a critical role of cooperativity in hydrophobic effect

The hydrophobic interaction drives nonpolar solutes to aggregate in aqueous solution, and hence plays a critical role in many fundamental processes in nature. An important property intrinsic to hydrophobic interaction is its cooperative nature, which is originated from the collective motions of water hydrogen bond networks surrounding hydrophobic solutes. This property is widely believed to enhance the formation of hydrophobic core in proteins. However, cooperativity in hydrophobic interactions has not been successfully characterized by experiments. Here, we quantify cooperativity in hydrophobic interactions by real-time monitoring the aggregation of hydrophobic solute (hexaphenylsilole, HPS) in a microfluidic mixer. We show that association of a HPS molecule to its aggregate in water occurs at sub-microsecond, and the free energy change is −5.8 to −13.6 kcal mol−1. Most strikingly, we discover that cooperativity constitutes up to 40% of this free energy. Our results provide quantitative evidence for the critical role of cooperativity in hydrophobic interactions.

Density Functional Theory and Car-Parrinello Molecular Dynamics Study of the Hydrogen-Producing Mechanism of the Co(dmgBF2)2 and Co(dmgH)2 Cobaloxime Complexes in Acetonitrile-Water Solvent.

The catalytic hydrogen-producing processes of two prototypical cobaloxime catalysts, Co(dmgBF2)2 (dmgBF2 = difluoroboryl-dimethylglyoxime) and Co(dmgH)2 (dmgH = dimethylglyoxime), were studied by density functional theory (DFT) and Car-Parrinello molecular dynamics (CPMD) simulations in the explicit acetonitrile-water solvent. Our study demonstrates the key role of water molecules as shuttles to deliver protons to the cobalt active centers of these catalysts. However, the transfer of protons to the cobalt centers also competes with the diffusion of the proton away from the complex via the hydrogen bond network of water. Protons were found to react with the oxygen of the side group of Co(dmgH)2, while a similar reaction was not observed for Co(dmgBF2)2. This explains the experimentally observed relative instability of Co(dmgH)2 in the acidic medium. The rate-limiting step of the hydrogen-producing process was found to be the first proton transfer to the cobalt center for both cobaloxime complexes. Structural and electron population analysis was carried out to provide insight into the origin of the difference of the proton transfer free-energy barriers of these two cobalt complexes. Our study has contributed to the key microscopic understanding of the hydrogen-producing process by this class of catalysts.