Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydration”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18

Anhydrous vs Hydrated f-Element Acetate Polymers Dictated by the Stoichiometry of Protic Acidic/Basic Azole Mixtures

Continuing our investigations of ionic liquid (IL) based routes to a library of f-element/soft donor complexes which could be studied crystallographically, we have explored the dissolution of f-element salts in protic imidazole-based ILs containing only soft donors at high temperatures to drive off volatiles, including water and carboxylic or mineral acids. Here we present our results, reacting acidic and basic azoles in 1:3 or 1:1 stoichiometric compositions at elevated temperature, followed by saturation with Nd(OAc) 3 ·xH 2 O or Ce(OAc) 3 ·xH 2 O, which led to 13 new metal–acetate polymeric complexes identified by single-crystal X-ray diffraction. We found that the diversity in coordination modes of the simple acetate ligand that interfere with substitution of the softer N donors led to several readily crystallizable complexes forming two distinct groups with respect to f-element interaction with the ionic liquid precursors. When the acidic/basic azole ratio was 1:3, acetate and a neutral basic azole were found to be coordinated to the metal centers but no water, although in one case (2) water was observed in the secondary coordination sphere: [Ce(μ 2 -OAc) 3 (C 1 im)] n (1, C 1 im = 1-methylimidazole), [Nd(μ 2 -(OAc) 3 (C 1 im)] n ·nH 2 O (2), [Ce(μ 2 -OAc) 3 (C 2 im)] n (3, C 2 im = 1-ethylimidazole), [Ln(μ 2 -OAc) 3 DMF] n (Ln = Nd (4), Ce (5); dimethylformamide (DMF) was substituted for the azole mixture), and [Nd(μ 2 -OAc) 3 (C 4 im)] n (6, C 4 im = 1-butylimidazole). However, when the stoichiometric ratio was 1:1, water was always observed coordinated to the metal ions with the acidic azole included in the structure as a solvate or cocrystal, despite a higher reaction temperature: [Nd(μ 2 -OAc) 3 (OH 2 )] n ·n(1,2,3-Taz) (7, 1,2,3-Taz = 1,2,3-triazole), [Ln(μ 2 -OAc) 3 (OH 2 )] n ·n(4,5-DCim) (Ln = Nd (8), Ce (9), 4,5-DCim = 4,5-dicyanoimidazole), [Ln(μ 2 -OAc) 3 (OH 2 )] n ·n(3,5-diNH 2 -1,2,4-Taz) (Ln = Nd (10), Ce (11), 3,5-diNH 2 -1,2,4-Taz = 3,5-diamino-1,2,4-triazole), [Ce(μ 2 -OAc) 3 (OH 2 )] n ·n(3-NH 2 -1,2,4-Taz) (12, 3-NH 2 -1,2,4-Taz = 3-amino-1,2,4-triazole), and [Nd(μ 2 -OAc) 3 (OH 2 )] n ·n(5-NH 2 -Tz) (13, 5-NH2-Tz = 5-aminotetrazole). All of the compounds retain the Ln:OAc– ratio of 1:3 and form 1D polymeric chains; however, they exhibit a variety of coordination modes affecting the degree of chain condensation. The isolation of both hydrated and anhydrous products revealed different abilities of the investigated soft N-donors to compete with O-donors finding their place in the coordination sphere of the lanthanide or in the crystal lattice.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Breakthrough Zn(II) Catalyst for Direct Air Capture Employing CO 2 Hydration

Direct air capture (DAC) represents a vital technology for atmospheric CO₂ remediation, but few studies have tested catalysts at dilute atmospheric CO₂ concentrations. Inspired by the carbonic anhydrase metalloenzyme, we report a catalytic DAC strategy employing robust zinc(II) enzyme mimics that enable efficient CO₂ sequestration pathways. A catalyst-mediated CO₂ hydration cycle in aqueous sorbents facilitates accelerated capture from dilute atmospheric air, thereby addressing the kinetic limitations observed in carbonate-based systems. Our developed complexes [ZnC1] and [ZnC2] enhance capture rates up to two-fold at millimolar concentrations and improve CO₂ mass transfer by 40-60% in 1 M K₂CO₃ sorbent under ambient conditions. These bench-stable, earth-abundant zinc catalysts operate effectively under dilute CO₂ concentrations, overcoming the kinetic limitations of conventional carbonate-based sorbents. As a result, mechanistic studies support a biomimetic catalytic cycle that facilitates rapid CO₂ conversion, demonstrating that catalyst-assisted DAC can enable energy-efficient, scalable carbon capture technologies.

Atmospheric chemistry↗

Water Defect Stabilizes the Bi 3+ Lone-Pair Electronic State Leading to an Unusual Aqueous Hydration Structure

In this work, the aqueous hydration structure of Bi 3+ is probed using a combination of extended X-ray absorption fine structure (EXAFS) spectroscopy and density functional theory (DFT) simulations of ion-water clusters and of condensed-phase solutions. Anomalous features in the EXAFS spectra are found to be associated with a highly asymmetric first-solvent water shell. The aqueous chemistry and structure of the Bi 3+ ion is dramatically controlled by the water stabilization of a lone-pair electronic state involving the mixed 6s and 6p orbitals. This leads to a distinct multi-modal distribution of waters in the first shell that are separated by about 0.2Å. The lone-pair structure is stabilized by a collective response of multiple waters that are localized near the lone-pair anti-bonding site. The findings indicate that the lone-pair stereochemistry of aqueous Bi 3+ plays a major role in the binding of water and ligands in aqueous solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Affinity for OH – Produces Four-Coordinated Zn 2+ Impurities in Hydrated Amorphous Calcium Carbonate

Using ab initio based molecular dynamics and electronic structure calculations, we show that Zn impurities in hydrated amorphous calcium carbonate (ACC) have a much lower coordination number than other divalent impurities due to covalent interactions between the 3d Zn shell and the oxygen atoms of the carbonate and water groups. Further, the local structure around Zn in ACC, including the predicted low coordination number, is confirmed by X-ray absorption spectroscopy of synthetic Zn-bearing ACC. The strong Zn–O chemical interaction leads to substantial water dissociation and slightly disrupts the hydrogen bonding network. Implications of Zn 2+ incorporation for ACC stability are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The H•/H– Redox Couple and Absolute Hydration Energy of H–

A supermolecule-continuum approach with water clusters up to n = 16 H2O molecules has been used to predict the absolute hydration free energies at 298 K ($\Delta$G hyd ) of both hydrogen (H•) and hydride (H - ) to be 4.6 and -78 ± 2 kcal/mol, respectively. These values are combined with a high accuracy prediction of the gas phase electron affinity ($\Delta$G gas,298K = -16.9 kcal/mol) to determine the aqueous electron affinity of H• of 99.5 ± 2 kcal/mol, which yields a reduction potential for H• vs SHE of -0.03 ± 0.15 V. This value is in agreement within 0.2 V with most estimates obtained using a wide variety of approaches. The results of this study can be used to improve the absolute hydricity scale in water which provides additional insights into how a putative hydride interacts with solvent, but do not change the ability to predict the relative reactivity of two species using relative hydricity scales.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

QM Investigation of Rare Earth Ion Interactions with First Hydration Shell Waters and Protein-Based Coordination Models

Here, conventional methods for extracting rare earth metals (REMs) from mined mineral ores are inefficient, expensive, and environmentally damaging. Recent discovery of lanmodulin (LanM), a protein that coordinates REMs with high-affinity and selectivity over competing ions, provides inspiration for new REM refinement methods. Here, we used quantum mechanical (QM) methods to investigate trivalent lanthanide cation (Ln 3+ ) interactions with coordination systems representing bulk solvent water and protein binding sites. Energy decomposition analysis (EDA) showed differences in the energetic components of Ln 3+ interaction with representatives of solvent (water, H 2 O) and protein binding sites (acetate, CH 3 COO – ), highlighting the importance of accurate description of electrostatics and polarization in computational modeling of REM interactions with biological and bioinspired molecules. Relative binding free energies were obtained for Ln 3+ with coordination complexes originating from binding sites in PDB structures of a lanthanum binding peptide (PDB entry 7CCO) and LanM, with explicit consideration of the first hydration shell waters, according to quasi-chemical theory (QCT). Beyond the first shell, the bulk solvent environment was represented with an implicit continuum model. Ln 3+ interactions with (H 2 O) 9 and both binding site models became more favorable, moving down the periodic series. This trend was more pronounced with the protein binding site models than with water, resulting in affinity increasing with periodic number, except for the last REM, Lu 3+ , which bound less favorably than the preceding element, Yb 3+ . Using the truncated 7CCO binding site model, the magnitude and trend of the experimental Ln 3+ relative binding free energies for the whole 7CCO peptide were reproduced. Conversely, the previously reported experimental data for LanM show a preference for the earlier lanthanides; this is likely due to longer-range interactions and cooperative effects, which are not represented by the reduced models. Using the truncated 7CCO binding site model, the magnitude and trend of the experimental Ln 3+ relative binding free energies for the whole 7CCO peptide were reproduced. In contrast to the previously reported experimental data for LanM, the peptide preferentially binds the earlier lanthanides. This difference likely arises due to longer-range interactions and cooperative effects not represented by the peptide. Further investigation of Ln 3+ interactions with whole proteins using polarizable molecular mechanics models with explicit solvent is warranted to understand the influence of longer-ranged interactions, cooperativity, and bulk solvent. Nevertheless, the present work provides new insights into Ln 3+ interactions with biomolecules and presents an effective computational platform for designing specific single-site REM binding peptides more efficiently.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Intermediate Species in the Crystallization of Sodium Aluminate Hydroxy Hydrates

In highly alkaline “water-in-salt” Na 2 O/Al 2 O 3 /H 2 O solutions where the monomeric Al(OH) 4 - anion dominates, isolation of transitional species that seed crystallization of sodium aluminate salt hydrates has been challenging. For example, discrimination of dimeric [e.g., Al 2 O(OH) 6 2- ] from species via 27 Al nuclear magnetic resonance (NMR) spectroscopy is limited via fast interconversion through hydrolysis and condensation reactions. Despite this, using magic-angle spinning NMR ( 27 Al MAS NMR) at high magnetic field strengths (14.1 and 20.0 T) on crystals of nonasodium bis(hexahydroxyaluminate) trihydroxide hexahydrate [NSA, Na 9 [Al(OH) 6 ] 2 ·3(OH)·6(H 2 O)], we observed a 27 Al resonance consistent with Na 2 Al 2 O(OH) 6 . This tetrahedral dimeric species was also identified in the mother liquor by Raman spectroscopy while remaining expectedly unresolved by high-field liquid-state 27 Al NMR. Its substantial abundance in the mother liquor and as either an interstitial or an adsorbate on precipitated solids suggest that it plays a crucial role in NSA crystallization. Here, our detailed characterization, which includes single-pulse 27 Al MAS NMR, nutation 27 Al MAS NMR, and multiple quantum 27 Al MAS NMR, provides new molecular insights into how aluminate transitions from tetrahedral in solution to octahedral in solid. This understanding could ultimately lead to better predictions of aluminum phase change behavior in caustic sodium hydroxide with wide-ranging technological and environmental implications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Inelastic Neutron Spectra of Uranium Tetrafluoride Hydrate, UF 4 (H 2 O) 2.5

Uranium tetrafluoride hydrate (UFH) is formed by immersing anhydrous UF4 under water for 12 h. UFH is therefore clearly a chemical species of environmental concern, as anhydrous UF4 is an intermediate uranium form in the nuclear fuel cycle. We use inelastic neutron scattering (INS) to probe the full vibrational spectra of UFH and its deuterated analogue in an effort to improve the fundamental understanding of its vibrational spectra. Coupled with density functional theory (DFT) calculations, the first for this compound, and full spectral modeling, we generate the complete vibrational spectra of UFH and compare them to prior optical spectroscopic results. In particular, the combination of DFT with INS allows us to identify multiple distinct chemical environments in the water bending and OH stretching regions. Whereas the water molecules directly bound to the U atoms execute OH stretching around 3600 cm –1 , a second class of H-bonded waters vibrate below 3000 cm –1 , an indicator of strong H bonding. In addition, a class of librational water modes are observed between 400 and 900 cm–1, which themselves can be separated in energy according to their chemical environments. Furthermore, measurements presented herein directly assist in the assignment of certain spectral features in the infrared spectrum and will inform future investigations of UFH for environmental or forensics purposes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Site-Selective Atomic Layer Deposition on Rutile TiO 2 : Selective Hydration as a Route to Target Point Defects

Routes to area-and especially site-selective atomic layer deposition (ALD) remain an enticing challenge in precision surface science, despite the potentially game-changing capability for many energy applications. An unparalleled level of surface reaction control is required to direct ALD to select sites on the same nominal material, for example, targeted growth on distinct phases, facets, step-edges, and/or defects. However, as a sequential surface synthesis method, ALD is uniquely suited to these challenges, including the possibility of selective deposition at defective surface atom arrangements. Here, we computationally identify conditions for site-selective ALD through hydration of surface defects, including oxygen vacancies and titanium interstitials on low-index rutile TiO 2 facets. First-principles computation is used to predict, as a function of temperature, the hydroxylation of defects that are targeted by proton-exchange-mediated ALD processes. In situ ellipsometric measurements of ALD Al 2 O 3 nucleation on TiO 2 (110) single crystals prepared with and without abundant oxygen vacancies demonstrate striking contrast, corroborating computational predictions and revealing a mechanistically clear path to site-selective ALD.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Guest-Host Interactions in Clathrate Hydrates: Benchmark MP2 and CCSD(T)/CBS Binding Energies of CH4, CO2 and H2S in (H2O)20 Cages

We present benchmark binding energies of naturally occurring gas molecules CH4, CO2, and H2S in the small cage, namely the pentagonal dodecahedron (512) (H2O)20, which is one of the constituent cages of the 3 major lattices (structures I, II and H) of clathrate hydrates. These weak interactions require higher levels of electron correlation and converge slowly with increasing basis set to the Complete Basis Set (CBS) limit, necessitating the use of large basis sets up to the augcc- pV5Z and subsequent correction for Basis Set Superposition Error (BSSE). For the host hollow (H2O)20 cages we have identified a most stable isomer with binding energy of -200.8 ± 2.1 kcal/mol at the CCSD(T)/CBS limit (-199.2 ± 0.5 kcal/mol at the MP2/CBS limit). Additionally, we report converged second order Moller-Plesset (MP2) CBS binding energies for the encapsulation of guests in the (H2O)20 cage of -4.3 ± 0.1 for CH4@(H2O)20, -6.6 ± 0.1 for CO2@(H2O)20 and -8.5 ± 0.1 kcal/mol for H2S@(H2O)20, respectively. For CH4@(H2O)20, exhibiting the weakest encapsulation affinity among the three, we report CCSD(T)/aug-cc-pVTZ binding energies and, based on them, a CCSD(T)/CBS estimate of -4.75 ± 0.1 kcal/mol. To the best of our knowledge, the CCSD(T)/aug-cc-pVTZ calculation for CH4@(H2O)20 is the largest one reported to date (168 valence electrons, 1978 basis functions and the correlation of 84 doubly occupied and 1873 virtual orbitals) and required a scalable implementation of the (T) module on 6144 nodes (350208 cores) of the “Cori” supercomputer at the National Energy Research Supercomputing Center (NERSC) for a total execution time of 195 minutes (for the (T) part). These efficient scalable implementations of highly correlated methods offer the capability to obtain long-lasting benchmarks of intermolecular interactions in complex systems. They also provide a path towards parametrizing classical potentials needed to study the dynamical and transport properties in these complex systems as well as assess the accuracy of lower scaling electronic structure methods such as Density Functional Theory (DFT) and MP2 including its spin-biased variants.

Heindel, Joseph↗

Near-Quantitative Predictions of the First-Shell Coordination Structure of Hydrated First-Row Transition Metal Ions Using K-Edge X-ray Absorption Near-Edge Spectroscopy

Understanding the solvation structure of transition metal ions is important due to its broad impact on the kinetics, stability, and reactivity in a variety of applications in geochemistry, biochemistry, energy storage, and environmental chemistry. Using water as the common ligand, we study the X-ray absorption pre-edge and near-edge spectra at the K-edge of a near-complete series of hydrated first-row transition metal ions with d-orbital occupancy ranging from d 2 to d 10 . Starting with optimized structures that were derived from an explicit solvation treatment at the density functional theory (DFT) level, we then compute the pre-edge X-ray absorption spectra at the metal ion K-edge at the time dependent density functional theory (TDDFT) and restricted active-space second-order perturbation theory (RASPT2) levels of theory. TDDFT calculations provide accurate results for spectra that are dominated by single excitations, while significant improvements were obtained with RASPT2 calculations that correctly distinguish between singly and doubly excited states, where relevant, with quantitative accuracy compared with experiment. We analyze and assign the dierent pre-edge features for each metal ion in order to reveal the impact of the variations in the d orbital occupancy on the rst-shell coordination environment. For completeness, we also report the lowest energy ligand-field d-d transitions for all the transition metal aqua ions considered in this study using complete active space second order perturbation theory (CASPT2).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Birth of the Hydrated Electron via Charge-Transfer-to-Solvent Excitation of Aqueous Iodide

A primary means to generate hydrated electrons in laboratory experiments is excitation to the charge-transfer-to-solvent (CTTS) state of a solute such as I – (aq), but this initial step in the genesis of e – (aq) has never been simulated directly using ab initio molecular dynamics. Here, we report the first such simulations, combining ground- and excited-state simulations of I – (aq) with a detailed analysis of fluctuations in the Coulomb potential experienced by the nascent solvated electron. What emerges is a two-step picture of the evolution of e – (aq) starting from the CTTS state: I – (aq) + hν → I –* (aq) → I • (aq) + e – (aq). Notably, the equilibrated ground state of e – (aq) evolves from I –* (aq) without any nonadiabatic transitions, simply as a result of solvent reorganization. The methodology used here should be applicable to other photochemical electron transfer processes in solution, an important class of problems directly relevant to photocatalysis and energy transfer.

Carter-Fenk, Kevin↗

Ionic Associations and Hydration in the Electrical Double Layer of Water-in-Salt Electrolytes

Water-in-Salt-Electrolytes (WiSEs) are an exciting class of concentrated electrolytes finding applications in energy storage devices because of their expanded electrochemical stability window, good conductivity and cation transference number, and fire-extinguishing properties. These distinct properties are thought to originate from the presence of an anion-dominated ionic network and interpenetrating water channels for cation transport, which indicates that associations in WiSEs are crucial to understanding their properties. Currently, associations have mainly been investigated in the bulk, while little attention has been given to the electrolyte structure near electrified interfaces. Here, we develop a theory for the electrical double layer (EDL) of WiSEs, where we consistently account for the thermoreversible associations of species into Cayley tree aggregates. The theory predicts an asymmetric structure of the EDL. At negative voltages, hydrated Li + dominates, and cluster aggregation is initially slightly enhanced before disintegration at larger voltages. At positive voltages, when compared to the bulk, clusters are strictly diminished. Performing atomistic molecular dynamics (MD) simulations of the EDL of WiSE provides EDL data for validation and bulk data for parametrization of our theory. Validating the predictions of our theory against MD showed good qualitative agreement. Furthermore, we performed electrochemical impedance measurements to determine the differential capacitance of the studied LiTFSI WiSE and also found reasonable agreement with our theory. Overall, the developed approach can be used to investigate ionic aggregation and solvation effects in the EDL, which, among other properties, can be used to understand the precursors for solid-electrolyte interphase formation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Pathways of Photocatalytic Oxidation of Formic Acid on Dry and Hydrated Anatase TiO 2 Surfaces

The photocatalytic oxidation of formic acid (FA), which is one of the most abundant volatile organic compounds, is a promising air remediation technology inspired by nature. However, the detailed mechanism of this photocatalytic reaction on the surface of TiO 2 , a typical photocatalyst, is not yet well-understood. In this work, we present a computational mechanistic study of the thermal vs photocatalytic oxidation of FA on dry and hydrated anatase TiO 2 (101) surfaces, based on periodic hybrid density functional theory (DFT) calculations, in which the photooxidation is treated as an excited-state process in a constrained triplet spin state. We first compare the adsorption modes of FA on the anatase (101) surface in the ground and excited states, followed by identification of the corresponding reaction intermediates that lead to the formation of CO 2 . We unveil the pivotal role of photogenerated holes localized at surface under-coordinated oxygen sites in mediating the C−H bond cleavage, thereby promoting CO 2 formation through a highly stable intermediate and an exergonic reaction step. Further investigation of the effect of coadsorbed water molecules shows that hydrogen bonding with water stabilizes FA in a monodentate configuration. This is favored over the unreactive bidentate structure that is the most stable under dry conditions, thus providing insight into the experimentally observed increase of the reaction rate in the presence of water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Pathways to High-Performance Salt Hydrate Thermochemical Energy Storage Materials and Systems

Thermochemical materials (TCMs) based on salt hydrates are promising for thermal energy storage as they combine high energy densities with low reaction temperatures. However, their adoption is hindered by poor structural integrity and degradation under hygrothermal cycling. Storage performance is governed not only by the chemical reaction, but also by the coupled thermo-chemo-mechanical behavior that evolves with cycling. Understanding and controlling this coupling across length scales (material-to-reactor) is necessary to improve TCM stability and lifetime. In this perspective, we discuss the shortcomings of current characterization approaches and emphasize the need for measuring transport properties and structural transformations using in situ techniques that capture the dynamic evolution of these materials. We also outline opportunities for multiscale modeling frameworks that link thermodynamics and mechanics, enabling predictive evaluation of composite architectures designed for cycling stability. We conclude by identifying research questions that must be addressed to transform TCMs into viable energy storage technologies.

Composites↗

Leveraging Hydration Forces for Size-Specific Nanoparticle Enrichment with a Redox-Responsive Silica-Binding Elastin-Like Polypeptide

Elastin-like polypeptides (ELPs) are low-complexity proteins that coacervate above a characteristic lower critical solution temperature (LCST). While the thermoresponsiveness of ELPs has been widely exploited in the biomedical and biomaterials fields, their ability to mediate nanoparticle assembly below their transition temperature remains largely unexplored. Here, we show that unmodified ELPs induce the reversible flocculation of silica nanoparticles (SiNPs) by forming backbone hydrogen bonds with surface silanols. Interparticle bridging is modulated by ELP length and concentration and by the presence of N- and C-terminal anchoring groups such as a cysteine residue and a Car9 silica-binding peptide. Using a redox-responsive fusion protein consisting of disulfide-bonded ELP domains terminated by Car9 segments, we stabilize 20 nm SiNPs under oxidizing conditions while triggering particle flocculation upon addition of reductant. We find that SiNP sedimentation under reducing conditions exhibits a sharp dependency on particle size that arises from the curvature-dependent structure of surface silanols. While the isolated silanols of SiNPs smaller than 30 nm are efficiently engaged by the ELP domains of Car9-anchored proteins, repulsion forces associated with the presence of a layer of molecular water together with increased electrostatic repulsion preclude efficient engagement of H-bonded silanols displayed on the surface of SiNPs larger than 60 nm. We harness these findings to selectively enrich SiNPs based on size and expand the concept to titania (TiO2) by demonstrating that rutile nanoparticles can be stabilized or sedimented with solid-binding ELPs by adjusting the solution pH to promote or discourage the formation of a hydration layer. These strategies should prove broadly useful for the separation of other oxides and their polymorphs and provide a tunable strategy for nanoparticle assembly and bioinspired colloidal design.

ELP↗

Solvent Dynamics Are Critical to Understanding Carbon Dioxide Dissolution and Hydration in Water

Simulations of carbon dioxide (CO 2 ) in water may aid in understanding the impact of its accumulation in aquatic environments and help advance technologies for carbon capture and utilization (via, e.g., mineralization). Quantum mechanical (QM) simulations based on static molecular models with polarizable continuum solvation poorly reproduce the energetics of CO 2 hydration to form carbonic acid in water, independent of the level of QM theory employed. Only with density-functional-theory-based molecular dynamics and rare-event sampling, followed by energy corrections based on embedded correlated wavefunction theory (in conjunction with density functional embedding theory), can a close agreement between theory and experiment be achieved. Such multilevel simulations can serve as benchmarks for simpler, less costly models, giving insight into potential errors of the latter. The strong influence of sampling/averaging over dynamical solvent configurations on the energetics stems from the difference in polarity of both the transition state and product (both polar) versus the reactant (nonpolar). When a solute undergoes a change in polarity during reaction, affecting its interaction with the solvent, careful assessment of the energetic contribution of the solvent response to this change is critical. Here, we show that static models (without structural sampling) that incorporate three explicit water molecules can yield far superior results than models with more explicit water molecules because fewer water molecules yield less configurational artifacts. Static models intelligently incorporating both explicit (molecules directly participating in the reaction) and implicit solvation, along with a proper QM theory, e.g., CCSD(T) for closed-shell systems, can close the accuracy gap between static and dynamic models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Switching of Protonation Sites in Hydrated Nicotine via a Grotthuss Mechanism

The switching of the protonation sites in hydrated nicotine, probed by experimental infrared (IR) spectroscopy and theoretical ab initio calculations, is facilitated via a Grotthuss instead of a bimolecular proton transfer (vehicle) mechanism at the experimental temperature (T = 130 K) as unambiguously confirmed by experiments with deuterated water. In contrast, the bimolecular vehicle mechanism is preferred at higher temperatures (T = 300 K) as determined by theory. Further, the Grotthuss mechanism for the concerted proton transfer results in the production of nicotine’s bioactive and addictive pyrrolidine-protonated (Pyrro-H + ) protomer with just 5 water molecules. Theoretical analysis suggests that the concerted proton transfer occurs via hydrogen-bonded bridges consisting of a 3 water molecule “core” that connects the pyridine protonated (Pyri-H + ) with the pyrrolidine-protonated (Pyrro-H + ) protomers. Additional water molecules attached as acceptors to the hydrogen-bonded “core” bridge result in lowering the reaction barrier of the concerted proton transfer down to less than 6 kcal/mol, which is consistent with the experimental observations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗