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
- NSF-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
-
Cleavage of aromatic carbon–chlorine bonds is critical for the degradation of toxic industrial compounds. Here, we solved the X-ray crystal structure of chlorothalonil dehalogenase (Chd) from Pseudomonas sp. CTN-3, with 15 of its N-terminal residues truncated (Chd T ), using single-wavelength anomalous dispersion refined to 1.96 Å resolution. Chd has low sequence identity (<15%) compared with all other proteins whose structures are currently available, and to the best of our knowledge, we present the first structure of a Zn(II)-dependent aromatic dehalogenase that does not require a coenzyme. Chd T forms a “head-to-tail” homodimer, formed between two α-helices from each monomer, with three Zn(II)-binding sites, two of which occupy the active sites, whereas the third anchors a structural site at the homodimer interface. The catalytic Zn(II) ions are solvent-accessible via a large hydrophobic (8.5 × 17.8 Å) opening to bulk solvent and two hydrophilic branched channels. Each active-site Zn(II) ion resides in a distorted trigonal bipyramid geometry with His 117 , His 257 , Asp 116 , Asn 216 , and a water/hydroxide as ligands. A conserved His residue, His 114 , is hydrogen-bonded to the Zn(II)-bound water/hydroxide and likely functions as the general acid-base. We examined substrate binding by docking chlorothalonil (2,4,5,6-tetrachloroisophtalonitrile, TPN) into the hydrophobic channel and observed that the most energetically favorable pose includes a TPN orientation that coordinates to the active-site Zn(II) ions via a CN and that maximizes a π–π interaction with Trp 227 . On the basis of these results, along with previously reported kinetics data, we propose a refined catalytic mechanism for Chd-mediated TPN dehalogenation.more » « less
-
Tunneling field effect transistors (TFETs) have gained much interest in the previous decade for use in low power CMOS electronics due to their sub-thermal switching [1]. To date, all TFETs are fabricated as vertical nanowires or fins with long, difficult processes resulting in long learning cycle and incompatibility with modern CMOS processing. Because most TFETs are heterojunction TFETs (HJ-TFETs), the geometry of the device is inherently vertically because dictated by the orientation of the tunneling HJ, achieved by typical epitaxy. Template assisted selective epitaxy was demonstrated for vertical nanowires [2] and horizontally arranged nanorods [3] for III-V on Si integration. In this work, we report results on the area selective and template assisted epitaxial growth of InP, utilizing SiO2 based confined structures on InP substrates, which enables horizontal HJs, that can find application in the next generation of TFET devices. The geometries of the confined structures used are so that only a small area of the InP substrate, dubbed seed, is visible to the growth atmosphere. Growth is initiated selectively only at the seed and then proceeds in the hollow channel towards the source hole. As a result, growth resembles epitaxial lateral overgrowth from a single nucleation point [4], reaping the benefits of defect confinement and, contrary to spontaneous nanowire growth, allows orientation in an arbitrary, template defined direction. Indium phosphide 2-inch (110) wafers are used as the starting substrate. The process flow (Fig.1) consists of two plasma enhanced chemical vapor deposition (PECVD) steps of SiO2, appropriately patterned with electron beam lithography (EBL), around a PECVD amorphous silicon sacrificial layer. The sacrificial layer is ultimately wet etched with XeF2 to form the final, channel like template. Not shown in the schematic in Fig.1 is an additional, ALD deposited, 3 nm thick, alumina layer which prevents plasma damage to the starting substrate and is removed via a final tetramethylammonium hydroxide (TMAH) based wet etch. As-processed wafers were then diced and loaded in a Thomas Swan Horizontal reactor. Successful growth conditions found were 600°C with 4E6 mol/min of group III precursor, a V/III ratio of 400 and 8 lpm of hydrogen as carrier gas. Trimethylindium (TMIn) and tertiarybutylphosphine (TBP) were used as In and P precursors respectively. Top view SEM (Fig.2) confirms growth in the template thanks to sufficient Z-contrast despite the top oxide layer, not removed before imaging. TEM imaging shows a cross section of the confined structure taken at the seed hole (Fig.3). The initial growth interface suggests growth was initiated at the seed hole and atomic order of the InP conforms to the SiO2 template both at the seed and at the growth front. A sharp vertical facet is an encouraging result for the future development of vertical HJ based III-V semiconductor devices.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
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.3390/polym12010249
