Dataset Information

Effect of Water Content on Thermodynamic Properties of Compressed Hydrogen.

ABSTRACT: Force field-based molecular simulations were used to calculate thermal expansivities, heat capacities, and Joule-Thomson coefficients of binary (standard) hydrogen-water mixtures for temperatures between 366.15 and 423.15 K and pressures between 50 and 1000 bar. The mole fraction of water in saturated hydrogen-water mixtures in the gas phase ranges from 0.004 to 0.138. The same properties were calculated for pure hydrogen at 323.15 K and pressures between 100 and 1000 bar. Simulations were performed using the TIP3P and a modified TIP4P force field for water and the Marx, Vrabec, Cracknell, Buch, and Hirschfelder force fields for hydrogen. The vapor-liquid equilibria of hydrogen-water mixtures were calculated along the melting line of ice Ih, corresponding to temperatures between 264.21 and 272.4 K, using the TIP3P force field for water and the Marx force field for hydrogen. In this temperature range, the solubilities and the chemical potentials of hydrogen and water were obtained. Based on the computed solubility data of hydrogen in water, the freezing-point depression of water was computed ranging from 264.21 to 272.4 K. The modified TIP4P and Marx force fields were used to improve the solubility calculations of hydrogen-water mixtures reported in our previous study [Rahbari A.;J. Chem. Eng. Data2019, 64, 4103-4115] for temperatures between 323 and 423 K and pressures ranging from 100 to 1000 bar. The chemical potentials of ice Ih were calculated as a function of pressure between 100 and 1000 bar, along the melting line for temperatures between 264.21 and 272.4 K, using the IAPWS equation of state for ice Ih. We show that at low pressures, the presence of water has a large effect on the thermodynamic properties of compressed hydrogen. Our conclusions may have consequences for the energetics of a hydrogen refueling station using electrochemical hydrogen compressors.

SUBMITTER: Rahbari A

PROVIDER: S-EPMC8154567 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Protein-protein interactions (PPIs) have been identified as a vital regulator of cellular pathways and networks. However, the determinants that control binding affinity and specificity at protein surfaces are incompletely characterized and thus unable to be exploited for the purpose of developing PPI inhibitors to control cellular pathways in disease states. One of the key factors in intermolecular interactions that remains poorly understood is the role of water molecules and in particular the importance of solvent entropy. This factor is expected to be particularly important at protein surfaces, and the release of water molecules from hydrophobic regions is one of the most important drivers of PPIs. In this work, we have studied the protein surface of a mutant of the protein RadA to quantify the thermodynamics of surface water molecules. RadA and its human homologue RAD51 function as recombinases in the process of homologous recombination. RadA binds to itself to form oligomeric structures and thus contains a well-characterized protein-protein binding surface. Similarly, RAD51 binds either to itself to form oligomers or to the protein BRCA2 to form filaments. X-ray crystallography has determined that the same interface functions in both interactions. Work in our group has generated a partially humanized mutant of RadA, termed MAYM, which has been crystallized in the apo form. We studied this apo form of MAYM using a combination of molecular dynamics (MD) simulations and inhomogeneous fluid solvation theory (IFST). The method locates a number of the hydration sites observed in the crystal structure and locates hydrophobic sites where hydrophobic species are known to bind experimentally. The simulations also highlight the importance of the restraints placed on the protein in determining the results. Finally, the results identify a correlation between the predicted entropy of water molecules at a given site and the solvent-accessible surface area and suggest that correlations between water molecules only need to be considered for water molecules separated by less than 3.2 Å. The combination of MD and IFST has been used previously to study PPIs and represents one of the few existing methods to quantify solvent thermodynamics. This is a vital aspect of molecular recognition and one which we believe must be developed.

Project description:BACKGROUND:Gram-negative bacteria use periplasmic-binding proteins (bPBP) to transport nutrients through the periplasm. Despite immense diversity within the recognized substrates, all members of the family share a common fold that includes two domains that are separated by a conserved hinge. The hinge allows the protein to cycle between open (apo) and closed (ligated) conformations. Conformational changes within the proteins depend on a complex interplay of mechanical and thermodynamic response, which is manifested as an increase in thermal stability and decrease of flexibility upon ligand binding. RESULTS:We use a distance constraint model (DCM) to quantify the give and take between thermodynamic stability and mechanical flexibility across the bPBP family. Quantitative stability/flexibility relationships (QSFR) are readily evaluated because the DCM links mechanical and thermodynamic properties. We have previously demonstrated that QSFR is moderately conserved across a mesophilic/thermophilic RNase H pair, whereas the observed variance indicated that different enthalpy-entropy mechanisms allow similar mechanical response at their respective melting temperatures. Our predictions of heat capacity and free energy show marked diversity across the bPBP family. While backbone flexibility metrics are mostly conserved, cooperativity correlation (long-range couplings) also demonstrate considerable amount of variation. Upon ligand removal, heat capacity, melting point, and mechanical rigidity are, as expected, lowered. Nevertheless, significant differences are found in molecular cooperativity correlations that can be explained by the detailed nature of the hydrogen bond network. CONCLUSION:Non-trivial mechanical and thermodynamic variation across the family is explained by differences within the underlying H-bond networks. The mechanism is simple; variation within the H-bond networks result in altered mechanical linkage properties that directly affect intrinsic flexibility. Moreover, varying numbers of H-bonds and their strengths control the likelihood for energetic fluctuations as H-bonds break and reform, thus directly affecting thermodynamic properties. Consequently, these results demonstrate how unexpected large differences, especially within cooperativity correlation, emerge from subtle differences within the underlying H-bond network. This inference is consistent with well-known results that show allosteric response within a family generally varies significantly. Identifying the hydrogen bond network as a critical determining factor for these large variances may lead to new methods that can predict such effects.

Project description:We perform path-integral molecular dynamics (PIMD) simulations of H2O and D2O using the q-TIP4P/F model. Simulations are performed at P = 1 bar and over a wide range of temperatures that include the equilibrium (T≥ 273 K) and supercooled (210 ≤T < 273 K) liquid states of water. The densities of both H2O and D2O calculated from PIMD simulations are in excellent agreement with experiments in the equilibrium and supercooled regimes. We also evaluate important thermodynamic response functions, specifically, the thermal expansion coefficient αP(T), isothermal compressibility κT(T), isobaric heat capacity CP(T), and static dielectric constant ε(T). While these properties are in excellent [αP(T) and κT(T)] or semi-quantitative agreement [CP(T) and ε(T)] with experiments in the equilibrium regime, they are increasingly underestimated upon further cooling. It follows that the inclusion of nuclear quantum effects in PIMD simulations of (q-TIP4P/F) water is not sufficient to reproduce the anomalous large fluctuations in density, entropy, and electric dipole moment characteristic of supercooled water. It has been hypothesized that water may exhibit a liquid-liquid critical point (LLCP) in the supercooled regime at P > 1 bar and that such a LLCP generates a maximum in CP(T) and κT(T) at 1 bar. Consistent with this hypothesis and in particular, with experiments, we find a maximum in the κT(T) of q-TIP4P/F light and heavy water at T≈ 230-235 K. No maximum in CP(T) could be detected down to T≥ 210 K. We also calculate the diffusion coefficient D(T) of H2O and D2O using the ring-polymer molecular dynamics (RPMD) technique and find that computer simulations are in remarkable good agreement with experiments at all temperatures studied. The results from RPMD/PIMD simulations are also compared with the corresponding results obtained from classical MD simulations of q-TIP4P/F water where atoms are represented by single interacting sites. Surprisingly, we find minor differences in most of the properties studied, with CP(T), D(T), and structural properties being the only (expected) exceptions.

Project description:One of the important parameters in water management of proton exchange membranes is the electro-osmotic drag (EOD) coefficient of water. The value of the EOD coefficient is difficult to justify, and available literature data on this for Nafion membranes show scattering from in experiments and simulations. Here, we use a classical all-atom model to compute the EOD coefficient and thermodynamic properties of water from molecular dynamics simulations for temperatures between 330 and 420 K, and for different water contents between λ = 5 and λ = 20. λ is the ratio between the moles of water molecules to the moles of sulfonic acid sites. This classical model does not capture the Grotthuss mechanism; however, it is shown that it can predict the EOD coefficient within the range of experimental values for λ = 5 where the vehicular mechanism dominates proton transfer. For λ > 5, the Grotthuss mechanism becomes dominant. To obtain the EOD coefficient, average velocities of water and ions are computed by imposing different electric fields to the system. Our results show that the velocities of water and hydronium scale linearly with the electric field, resulting in a constant ratio of ca. 0.4 within the error bars. We find that the EOD coefficient of water linearly increases from 2 at λ = 5 to 8 at λ = 20 and the results are not sensitive to temperature. The EOD coefficient at λ = 5 is within the range of experimental values, confirming that the model can capture the vehicular transport of protons well. At λ = 20, due to the absence of proton hopping in the model, the EOD coefficient is overestimated by a factor of 3 compared to experimental values. To analyze the interactions between water and Nafion, the partial molar enthalpies and partial molar volumes of water are computed from molecular dynamics simulations. At different water uptakes, multiple linear regression is used on raw simulation data within a narrow composition range of water inside the Nafion membrane. The partial molar volumes and partial molar excess enthalpies of water asymptotically approach the molar volumes and molar excess enthalpies of pure water for water uptakes above 5. This confirms the model can capture the bulklike behavior of water in the Nafion at high water uptakes.

Dataset Information

Effect of Water Content on Thermodynamic Properties of Compressed Hydrogen.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets