An official website of the United States government
Predicting the asymmetric structure and dynamics of solvated hydroxide and hydronium in water from ab initio molecular dynamics (AIMD) has been a challenging task. The difficulty mainly comes from a lack of accurate and efficient exchange–correlation functional in elucidating the amphiphilic nature and the ubiquitous proton transfer behaviors of the two ions. By adopting the strongly constrained and appropriately normed (SCAN) meta-generalized gradient approximation functional in AIMD simulations, we systematically examine the amphiphilic properties, the solvation structures, the electronic structures, and the dynamic properties of the two water ions. In particular, we compare these results to those predicted by the PBE0-TS functional, which is an accurate yet computationally more expensive exchange–correlation functional. We demonstrate that the general-purpose SCAN functional provides a reliable choice for describing the two water ions. Specifically, in the SCAN picture of water ions, the appearance of the fourth and fifth hydrogen bonds near hydroxide stabilizes the pot-like shape solvation structure and suppresses the structural diffusion, while the hydronium stably donates three hydrogen bonds to its neighbors. We apply a detailed analysis of the proton transfer mechanism of the two ions and find the two ions exhibit substantially different proton transfer patterns. Our AIMD simulations indicate that hydroxide diffuses more slowly than hydronium in water, which is consistent with the experimental results.
more »
« less
Thermodynamic and Transport Properties of Tetrabutylphosphonium Hydroxide and Tetrabutylphosphonium Chloride–Water Mixtures via Molecular Dynamics Simulation
Thermodynamic, structural, and transport properties of tetrabutylphosphonium hydroxide (TBPH) and tetrabutylphosphonium chloride (TBPCl)–water mixtures have been investigated using all-atom molecular dynamics simulations in response to recent experimental work showing the TBPH–water mixtures capability as a cellulose solvent. Multiple transitional states exist for the water—ionic liquid (IL) mixture between 70 and 100 mol% water, which corresponds to a significant increase in water hydrogen bonds. The key transitional region, from 85 to 92.5 mol% water, which coincides with the mixture’s maximum cellulose solubility, reveals small and distinct water veins with cage structures formed by the TBP+ ions, while the hydroxide and chloride ions have moved away from the P atom of TBP+ and are strongly hydrogen bonded to the water. The maximum cellulose solubility of the TBPH–water solution at approximately 91.1 mol% water, appears correlated with the destruction of the TBP’s interlocking structure in the simulations, allowing the formation of water veins and channeling structures throughout the system, as well as changing from a subdiffusive to a near-normal diffusive regime, increasing the probability of the IL’s interaction with the cellulose polymer. A comparison is made between the solution properties of TBPH and TBPCl with those of alkylimidazolium-based ILs, for which water appears to act as anti-solvent rather than a co-solvent.
more »
« less
- Award ID(s):
- 1726534
- PAR ID:
- 10196964
- Date Published:
- Journal Name:
- Polymers
- Volume:
- 12
- Issue:
- 1
- ISSN:
- 2073-4360
- Page Range / eLocation ID:
- 249
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Deep eutectic solvents (DES) are compounds of a hydrogen bond donor (HBD) and a hydrogen bond acceptor (HBA) that contain a depressed melting point compared to their individual constituents. DES have been studied for their use as carbon capture media and biogas upgrading. However, contaminants’ presence in biogas might affect the carbon capture by DES. In this study, conductor-like screening model for real solvents (COSMO-RS) was used to determine the effect of temperature, pressure, and selective contaminants on five DES’ namely, choline chloride-urea, choline chloride-ethylene glycol, tetra butyl ammonium chloride-ethylene glycol, tetra butyl ammonium bromide-decanoic acid, and tetra octyl ammonium chloride-decanoic acid. Impurities studied in this paper are hydrogen sulfide, ammonia, water, nitrogen, octamethyltrisiloxane, and decamethylcyclopentasiloxane. At infinite dilution, CO2 solubility dependence upon temperature in each DES was examined by means of Henry’s Law constants. Next, the systems were modeled from infinite dilution to equilibrium using the modified Raoults’ Law, where CO2 solubility dependence upon pressure was examined. Finally, solubility of CO2 and CH4 in the various DES were explored with the presence of varying mole percent of selective contaminants. Among the parameters studied, it was found that the HBD of the solvent is the most determinant factor for the effectiveness of CO2 solubility. Other factors affecting the solubility are alkyl chain length of the HBA, the associated halogen, and the resulting polarity of the DES. It was also found that choline chloride-urea is the most selective to CO2, but has the lowest CO2 solubility, and is the most polar among other solvents. On the other hand, tetraoctylammonium chloride-decanoic acid is the least selective, has the highest maximum CO2 solubility, is the least polar, and is the least affected by its environment.more » « less
-
D-Mannosamine hydrochloride (2-amino-2-deoxy-D-mannose hydrochloride), C 6 H 14 NO 5 + ·Cl − , (I), crystallized from a methanol/ethyl acetate/ n -hexane solvent mixture at room temperature in a 4 C 1 chair conformation that is slightly distorted towards the C3,O5 B form. A comparison of the structural parameters of (I) with the corresponding parameters in α-D-glucosamine hydrochloride, (II), and β-D-galactosamine hydrochloride, (III)/(III′), was undertaken to evaluate the effects of ionic hydrogen bonding on structural properties. Three types of ionic hydrogen bonds are present in the crystals of (I)–(III)/(III′), i.e. N + —H...O, N + —H...Cl − , and O—H...Cl − . The exocyclic structural parameters in (I), (II), and (III)/(III′) appear to be most influenced by this bonding, especially the exocyclic hydroxy groups, which adopt eclipsed conformations enabled by ionic hydrogen bonding to the chloride anion. Anomeric disorder was observed in crystals of (I), with an α:β ratio of 37:63. However, anomeric configuration appears to exert minimal structural effects; that is, bond lengths, bond angles, and torsion angles are essentially identical in both anomers. The observed disorder at the anomeric C atom of (I) appears to be caused by the presence of the chloride anion and atom O3 or O4 in proximal voids, which provide opportunities for hydrogen bonding to atom O1 in both axial and equatorial orientations.more » « less
-
The transport of excess protons and hydroxide ions in water underlies numerous important chemical and biological processes. Accurately simulating the associated transport mechanisms ideally requires utilizing ab initio molecular dynamics simulations to model the bond breaking and formation involved in proton transfer and path-integral simulations to model the nuclear quantum effects relevant to light hydrogen atoms. These requirements result in a prohibitive computational cost, especially at the time and length scales needed to converge proton transport properties. Here, we present machine-learned potentials (MLPs) that can model both excess protons and hydroxide ions at the generalized gradient approximation and hybrid density functional theory levels of accuracy and use them to perform multiple nanoseconds of both classical and path-integral proton defect simulations at a fraction of the cost of the corresponding ab initio simulations. We show that the MLPs are able to reproduce ab initio trends and converge properties such as the diffusion coefficients of both excess protons and hydroxide ions. We use our multi-nanosecond simulations, which allow us to monitor large numbers of proton transfer events, to analyze the role of hypercoordination in the transport mechanism of the hydroxide ion and provide further evidence for the asymmetry in diffusion between excess protons and hydroxide ions.more » « less
-
Abstract The electrochemical stability window of water is known to vary with the type and concentration of dissolved salts. However, the underlying influence of ions on the thermodynamic stability of aqueous solutions has not been fully understood. Here, we investigated the electrolytic behaviors of aqueous electrolytes as a function of different ions. Our findings indicate that ions with high ionic potentials, i.e., charge density, promote the formation of their respective hydration structures, enhancing electrolytic reactions via an inductive effect, particularly for small cations. Conversely, ions with lower ionic potentials increase the proportion of free water molecules—those not engaged in hydration shells or hydrogen‐bonding networks—leading to greater electrolytic stability. Furthermore, we observe that the chemical environment created by bulky ions with lower ionic potentials impedes electrolytic reactions by frustrating the solvation of protons and hydroxide ions, the products of oxygen evolution reaction (OER) and hydrogen evolution reaction (HER), respectively. We found that the solvation of protons plays a more substantial role than that of hydroxide, which explains a greater shift for OER than for HER, a puzzle that cannot be rationalized by the notion of varying O−H bond strengths of water. These insights will help the design of aqueous systems.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.3390/polym12010249
