Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Hydroxyls”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Simulation insights into wetting properties of hydrogen-brine-clay for hydrogen geo-storage

Hydrogen geo-storage is attracting substantial interdisciplinary interest as a cost-effective and sustainable option for medium- and long-term storage. Hydrogen can be stored underground in diverse formations, including aquifers, salt caverns, and depleted oil and gas reservoirs. The wetting dynamics of the hydrogen-brine-rock system are critical for assessing both structural and residual storage capacities, and ensuring containment safety. Through molecular dynamics simulations, we explore how varying concentrations of cushion gases (CO 2 or CH 4 ) influence the wetting properties of hydrogen-brine-clay systems under geological conditions (15 MPa and 333 K). We employed models of talc and the hydroxylated basal face of kaolinite (kaoOH) as clay substrates. Our findings reveal that the effect of cushion gases on hydrogen-brine-clay wettability is strongly dependent on the clay-brine interactions. Notably, CO 2 and CH 4 reduce the water wettability of talc in hydrogen-brine-talc systems, while exerting no influence on the wettability of hydrogen-brine-kaoOH systems. Detailed analysis of free energy of cavity formation near clay surfaces, clay-brine interfacial tensions, and the Willard-Chandler surface for gas-brine interfaces elucidate the molecular mechanisms underlying wettability changes. Our simulations identify empirical correlations between wetting properties and the average free energy required to perturb a flat interface when clay-brine interactions are less dominant. Here, our thorough thermodynamic analysis of rock-fluid and fluid-fluid interactions, aligning with key experimental observations, underscores the utility of simulated interfacial properties in refining contact angle measurements and predicting experimentally relevant properties. These insights significantly enhance the assessment of gas geo-storage potential. Prospectively, the approaches and findings obtained from this study could form a basis for more advanced multiscale simulations that consider a range of geological and operational variables, potentially guiding the development and improvement of geo-storage systems in general, with a particular focus on hydrogen storage.

25 ENERGY STORAGE↗
DOE OSTI2020

Controllable depolymerization of lignin using carbocatalyst graphene oxide under mild conditions

Many efforts have been devoted to depolymerizing lignin to high-value aromatic monomers. However, a controllable and complete lignin inter-unit bond cleavage remains difficult because of the randomness in lignin structure and low selectivity in catalysis. Herein, we demonstrate that a sole carbocatalyst graphene oxide (GO) can selectively oxidize aryl- and alkyl-hydroxyl groups of lignin without the presence of any additional oxidizers, leading to lignin degradation under mild conditions. The results from lignin dimers, lignin oligomers, and organosolv lignin indicate that the depolymerization is catalyzed via an oxidative mechanism through an enol ether intermediate, a stable intermediate that has been chemically isolated here. Quantum chemistry calculations at the density functional theory level indicate that the phenolic –OH group in lignin facilitates the formation of radical cations and thus prohibits ring opening reactions, leading to the complete conversion of lignin dimers to guaiacol through the enol ether intermediate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Oxidation rates and redox stabilization of ferrous iron in trioctahedral smectites

Iron(II)-bearing trioctahedral smectites (saponites) form during anoxic alteration of basaltic rock. They are predicted to have been widespread on the early Earth and are observed in the oceanic subsurface today. Smectite structures, including the occupancy of sites in the octahedral sheet, affect iron redox behavior but the rates and products of trioctahedral smectite oxidation have been largely unexplored to date. In this study we synthesized two Fe(II)-bearing trioctahedral smectites, one moderate (22 wt% Fe) and one high (27 wt% Fe) in iron content. We then examined the rate, extent, and products of their oxidation by dissolved oxygen, nitrite, and hydrogen peroxide. Dissolved oxygen caused partial oxidation of Fe(II) in the smectites with 14 to 43% of Fe(II) unoxidized after 20 to 30 days of exposure. The rate and extent of oxidation correlated with the dissolved oxygen concentration and the Fe(II) content of the clay. The incomplete oxidation in these experiments is consistent with the mixed-valent trioctahedral smectites observed in oxidized natural samples but contrasts with the complete reoxidation by oxygen shown by chemically- or microbially-reduced dioctahedral smectites. Oxidation of structural Fe(II) by 5 mmol L -1 nitrite was negligible for the moderate-iron smectite and yielded only ~17% oxidation after 54 days of reaction for the high-iron smectite. Hydrogen peroxide caused rapid and near-complete oxidation of both clays. Powder X-ray diffraction, variable-temperature Mössbauer spectroscopy, and extended X-ray absorption fine structure spectroscopy together detected no crystalline or short-range-ordered secondary phases and show that oxidized iron remained in the trioctahedral smectite structure. The recalcitrant Fe(II) pool in oxidized trioctahedral smectites exists in less distorted sites than Fe(II) in the initial clays. Its unreactive nature at prolonged reaction times indicates an elevated redox potential generated by the local coordination environment. Slower oxidation rates create a larger recalcitrant Fe(II) pool, suggesting kinetic competition between oxidation and a process involved in redox stabilization, such as electron exchange between octahedral iron sites or deprotonation of hydroxyl groups in the structure. The resistance to complete oxidation of trioctahedral ferrous smectites and their full retention of iron demonstrates that transitions from anoxic to oxic conditions generate mixed-valence smectites rather than a mixture of new phases. Identifying the diagenetic products of mixed-valent trioctahedral smectites may provide an indicator in the rock record of past redox cycling. In conclusion, substantial portions of structural Fe(II) in trioctahedral smectites display slow abiotic oxidation kinetics and represent potential electron donors for both microaerophilic iron oxidizing and nitrate-reducing, iron-oxidizing microorganisms in altered mafic rocks and related settings.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

An engineered lactate oxidase based electrochemical sensor for continuous detection of biomarker lactic acid in human sweat and serum

Lactate levels in humans reveal intensity and duration of exertion and provide a critical readout for the severity of life-threatening illnesses such as pediatric sepsis. Using the lactate oxidase enzyme (Lox) from Aerococcus viridians, we demonstrated its functionality for lactate electrochemical sensing in physiological fluids in a lab setting. The structure and dynamics of LOx were validated by crystallography, X-ray scattering, and hydroxyl radical protein footprinting. This provided a validated protein template for understanding and designing an enzyme-based electrochemical sensing elements. Using this template, LOx enzyme variants were generated and compared. Comparison of the variants demonstrates that one exhibits effective lactate sensing at significantly reduced operating voltages. Additionally, we demonstrate that the four hexahistidine-tags on each enzyme tetramer are sufficient for immobilization to create a durable, functional sensor, with no need for a covalent attachment, enabling self-immobilization and eliminating the need for additional immobilization steps. The functionality of the LOx enzyme variants was verified at physiological lactate concentrations in both human serum (0–4 mM) and artificial sweat (0–100 mM) using 3-electrode setups for analysis of the three variants in parallel. Accuracy of measurement in both artificial sweat and human serum were high. Employing a microfluidic flow cell, we successfully monitored varying lactate levels in physiological fluids continuously over a 2h period. Overall, this optimized LOx enzyme, which self-immobilizes onto gold sensing electrodes, facilitates efficient and reliable lactate detection and continuous monitoring at reduced operating voltages suitable for further development towards commercial use.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Nanolignin-based high internal phase emulsions for efficient protection of curcumin against UV degradation

As an emulsifier, lignin exhibits excellent UV resistance on drug-loaded emulsion systems for drug delivery. However, due to the structural variation and complexity of lignins from various origins, their UV shielding performance varies with the techniques for lignin extraction, which impacts properties and the protection efficiency of lignin-based HIPEs (high internal phase emulsions). In this work, lignin nanoparticles, prepared from three lignin preparations of Eucalyptus, were used in HIPEs delivery systems to protect curcumin from degradation by UV radiation. Structures of the lignin preparations were characterized using 2D HSQC (heteronuclear single-quantum coherence) NMR (nuclear magnetic resonance), 31 P NMR, and GPC (gel permeation chromatography). The residual curcumin level after 36h UV exposure in the nanolignin-based HIPEs was over 72%, much higher than that (< 10% after 24 h UV exposure) in the oil phase without lignin, indicating that the nanolignin-based HIPEs with enhanced UV shielding ability protect curcumin better. Of the three lignin preparations, AL (alkali lignin), with the lowest molecular weight, highest contents of phenolic hydroxyl and carboxyl groups, and highest S/G ratio, displayed the best anti-UV radiation ability and the most uniform nanoparticle size.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Effects of the advanced organosolv pretreatment strategies on structural properties of woody biomass

Plants have evolved their lignocellulosic cell walls through complex structural and chemical mechanisms to protect itself against microbial attacks, which makes native lignocellulosic biomass recalcitrant to enzymatic deconstruction. Pretreatment is a crucial step in the biological conversion of biomass as it can render structural changes in the plant cell wall to reduce the biomass recalcitrance, thus enhancing its sugar release performance. There have been many efforts to develop effective pretreatment technologies to overcome the biomass recalcitrance with a primary focus on the efficient conversion of biomass carbohydrates to liquid fuels, while lignin is significantly underutilized despite its bulk amount and high-value opportunities. In this study, the effects of two recent organosolv pretreatment strategies, co-solvent enhanced lignocellulosic fractionation (CELF) and γ-Valerolactone (GVL) pretreatments, on physicochemical properties of poplar were investigated and compared with the effects of conventional ethanol organosolv pretreatment. Diverse physicochemical properties of biomass including chemical compositions, molecular weights of cellulose and lignin, aromatics and inter-unit linkages of lignin, lignin hydroxyl group contents, cellulose crystallinity, and accessible surface area of cellulose were analyzed before and after pretreatments. The results revealed how each organic solvent pretreatment system affected biomass structural characteristics and recalcitrance.

Meng, Xianzhi↗
DOE OSTI2020

The effect of lignin degradation products on the generation of pseudo-lignin during dilute acid pretreatment

Pseudo-lignin is an insoluble material that tends to deposit on acidic pretreated fiber surface and negatively impacts the biological conversion of biomass to valued added products. Employing a laboratory generated acidic hydrolysate solution from dilute acid pretreatment of bamboo residues, we sought to seek an improved understanding of pseudo-lignin formation from fragmented carbohydrates. In addition, we also performed experiments to investigate how lignin model compounds interact with carbohydrates en route to pseudo-lignin’s formation. The yields of pseudo-lignin generated from xylose (CXPL) were higher than those generated from glucose (CGPL). Gas chromatography-Mass spectrometry and Foloin-Ciocalteau analysis suggested that there were more aromatic compounds in hydrolyzate during CGPL’s formation. The presence of lignin phenolic model compounds had both a positive and negative effect upon quantitative yields of pseudo-lignin. Gel permeation chromatography and NMR analysis ( 31 P and 2D-HSQC NMR) of the pseudo-lignin samples generated in the presence and absence of 4-hydroxybenzoic acid (4-HA), employed as a lignin model compound, revealed drastic differences in molecular weight, changes in hydroxyl group content, and hydrocarbon bonds, suggesting significant differences in the pseudo-lignin reaction pathways when certain lignin model compounds are present. Finally, our findings demonstrate the effects of both hemicellulose carbohydrate makeup and lignin model compound impact on the formation of pseudo-lignin and its chemical properties, suggesting a potential way to minimize pseudo-lignin formation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Conversion of Loblolly pine biomass residues to bio-oil in a two-step process: Fast pyrolysis in the presence of zeolite and catalytic hydrogenation

In this work, Loblolly pine biomass residue was converted to bio-oil in a two-step process, consisting of 1) fast pyrolysis in the presence of zeolite ZSM-5 as a catalyst to produce pyrolysis oil, 2) hydrogenation of pyrolysis oil using formic acid as the hydrogen source in presence of Ru/activated carbon catalyst. Pyrolysis oils were analyzed by 13 C, 31 P and HSQC-NMR and the results revealed that the zeolite-induced catalytic fast pyrolysis process led to effective demethoxylation, producing more catechol and p-hydroxy-phenyl hydroxyl groups in the bio-oils, resulting in a decrease in the methoxyl group content by about 85 % and rich aromatic structures in the pyrolysis oils. The properties of pyrolysis oil with and without zeolite were in the bio-oil range. Hydrogenated pyrolysis oil showed that 79 % of the aromatic protons are eliminated and 87 % of protons are aliphatic in nature, with no oxygen attached to the α-carbon.

09 BIOMASS FUELS↗
DOE OSTI2020

Oxidation of phenolic compounds during autothermal pyrolysis of lignocellulose

Fast pyrolysis is traditionally defined as the rapid decomposition of organic material in the absence of oxygen to produce primarily a liquid product known as bio-oil. However, the introduction of small amounts of oxygen to the process holds prospects of internally generating the energy needed for pyrolysis. The present study investigates the partial oxidation of lignin-derived compounds during pyrolysis, which generates both carbon oxides and aromatic carbonyl compounds. Analysis of lignin derived phenolic compounds was performed to determine if the composition had changed under oxidative conditions. NMR analyses indicates aromatic carbonyls increased under oxidative conditions, with a corresponding decrease in phenolic hydroxyl groups. Finally, model phenolic compounds were pyrolyzed to help understand the role of partial oxidation during autothermal pyrolysis of lignocellulosic biomass.

09 BIOMASS FUELS↗
DOE OSTI2022

Corrosion-resistant polyurethane coatings from structure-homogenized biorefinery lignin through fractionation and oxypropylation

Using lignin for polyurethane (PU) synthesis is challenging because of its heterogeneous structure, i.e., broad molecule weight (MW) distribution and diverse hydroxyl (OH) functionality. This complexity influences the crosslinking reaction between lignin and isocyanate and affects the quality of PU products. Here, we propose a two-step strategy for lignin structure homogenization by lowering its MW and uniforming OH functionality (step 2). Step 1 involves a bioethanol fractionation process, yielding fractionated lignin with one order lower MW than the raw lignin while exhibiting better processability for the subsequent process. Step 2 is performed by an oxyalkylation process using fractionated lignin that converts phenolic and carboxylic OHs to aliphatic OHs. As a result, we obtained structure-homogenized lignin with good compatibility and crosslinking with hexamethylene diisocyanate for PU coating synthesis. Moreover, PU coatings prepared from structure-homogenized lignin display a defect-free microscopic structure and show better corrosion resistance than coatings fabricated from its raw lignin alternative. Overall, the structure homogenization process is a beneficial way of using lignin as biopolyol in PU synthesis toward the application-oriented design of corrosion-resistant lignin-based PU coatings.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Biophysical analysis of the plant-specific GIPC sphingolipids reveals multiple modes of membrane regulation

The plant plasma membrane (PM) is an essential barrier between the cell and the external environment, controlling signal perception and transmission. It consists of an asymmetrical lipid bilayer made up of three different lipid classes: sphingolipids, sterols, and phospholipids. The glycosyl inositol phosphoryl ceramides (GIPCs), representing up to 40% of total sphingolipids, are assumed to be almost exclusively in the outer leaflet of the PM. However, their biological role and properties are poorly defined. In this study, we investigated the role of GIPCs in membrane organization. Because GIPCs are not commercially available, we developed a protocol to extract and isolate GIPC-enriched fractions from eudicots (cauliflower and tobacco) and monocots (leek and rice). Lipidomic analysis confirmed the presence of trihydroxylated long chain bases and 2-hydroxylated very long-chain fatty acids up to 26 carbon atoms. The glycan head groups of the GIPCs from monocots and dicots were analyzed by gas chromatograph-mass spectrometry, revealing different sugar moieties. Multiple biophysics tools, namely Langmuir monolayer, ζ-Potential, light scattering, neutron reflectivity, solid state 2H-NMR, and molecular modeling, were used to investigate the physical properties of the GIPCs, as well as their interaction with free and conjugated phytosterols. We showed that GIPCs increase the thickness and electronegativity of model membranes, interact differentially with the different phytosterols species, and regulate the gel-to-fluid phase transition during temperature variations. These results unveil the multiple roles played by GIPCs in the plant PM.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Primer terminal ribonucleotide alters the active site dynamics of DNA polymerase η and reduces DNA synthesis fidelity

DNA polymerases catalyze DNA synthesis with high efficiency, which is essential for all life. Extensive kinetic and structural efforts have been executed in exploring mechanisms of DNA polymerases, surrounding their kinetic pathway, catalytic mechanisms, and factors that dictate polymerase fidelity. Recent time-resolved crystallography studies on DNA polymerase η (Pol η) and β have revealed essential transient events during the DNA synthesis reaction, such as mechanisms of primer deprotonation, separated roles of the three metal ions, and conformational changes that disfavor incorporation of the incorrect substrate. DNA-embedded ribonucleotides (rNs) are the most common lesion on DNA and a major threat to genome integrity. While kinetics of rN incorporation has been explored and structural studies have revealed that DNA polymerases have a steric gate that destabilizes ribonucleotide triphosphate binding, the mechanism of extension upon rN addition remains poorly characterized. Using steady-state kinetics, static and time-resolved X-ray crystallography with Pol η as a model system, we showed that the extra hydroxyl group on the primer terminus does alter the dynamics of the polymerase active site as well as the catalysis and fidelity of DNA synthesis. During rN extension, Pol η error incorporation efficiency increases significantly across different sequence contexts. Finally, our systematic structural studies suggest that the rN at the primer end improves primer alignment and reduces barriers in C2'-endo to C3'-endo sugar conformational change. Overall, our work provides further mechanistic insights into the effects of rN incorporation on DNA synthesis.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Human cytochrome P450 17A1 structures with metabolites of prostate cancer drug abiraterone reveal substrate-binding plasticity and a second binding site

Abiraterone acetate is a first-line therapy for castration-resistant prostate cancer. This prodrug is deacetylated in vivo to abiraterone, which is a potent and specific inhibitor of cytochrome P450 17A1 (CYP17A1). CYP17A1 performs two sequential steps that are required for the biosynthesis of androgens that drive prostate cancer proliferation, analogous to estrogens in breast cancer. Abiraterone can be further metabolized in vivo on the steroid A ring to multiple metabolites that also inhibit CYP17A1. Despite its design as an active-site-directed substrate analog, abiraterone and its metabolites demonstrate mixed competitive/noncompetitive inhibition. To understand their binding, we solved the X-ray structures of CYP17A1 with three primary abiraterone metabolites. Despite different conformations of the steroid A ring and substituents, all three bound in the CYP17A1 active site with the steroid core packed against the I helix and the A ring C3 keto or hydroxyl oxygen forming a hydrogen bond with N202 similar to abiraterone itself. The structure of CYP17A1 with 3-keto, 5α-abiraterone was solved to 2.0 Å, the highest resolution to date for a CYP17A1 complex. This structure had additional electron density near the F/G loop, which is likely a second molecule of the inhibitor and which may explain the noncompetitive inhibition. Mutation of the adjacent Asn52 to Tyr positions its side chain in this space, maintains enzyme activity, and prevents binding of the peripheral ligand. Collectively, our findings provide further insight into abiraterone metabolite binding and CYP17A1 function.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Characterization of lignin-degrading enzyme PmdC, which catalyzes a key step in the synthesis of polymer precursor 2-pyrone-4,6-dicarboxylic acid

Pyrone-2,4-dicarboxylic acid (PDC) is a valuable polymer precursor that can be derived from the microbial degradation of lignin. The key enzyme in the microbial production of PDC is 4-carboxy-2-hydroxymuconate-6-semialdehyde (CHMS) dehydrogenase, which acts on the substrate CHMS. We present the crystal structure of CHMS dehydrogenase (PmdC from Comamonas testosteroni) bound to the cofactor NADP, shedding light on its three-dimensional architecture, and revealing residues responsible for binding NADP. Using a combination of structural homology, molecular docking, and quantum chemistry calculations, we have predicted the binding site of CHMS. Key histidine residues in a conserved sequence are identified as crucial for binding the hydroxyl group of CHMS and facilitating dehydrogenation with NADP. Mutating these histidine residues results in a loss of enzyme activity, leading to a proposed model for the enzyme's mechanism. These findings are expected to help guide efforts in protein and metabolic engineering to enhance PDC yields in biological routes to polymer feedstock synthesis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Catalytic consequences of hydrogen addition events and solvent-adsorbate interactions during guaiacol-H 2 reactions at the H 2 O-Ru(0 0 0 1) interface

Catalytic reactions of biomass-derived phenolics and H 2 occur on transition metal surfaces via competitive C–O cleavage and ring saturation pathways, with both requiring multiple hydrogen addition events before forming their respective rate limiting transition states. These events are markedly affected by solvent chemical identity, with polar protic solvents ionizing hydrogen adatoms (H*) to interfacial protons (H + ) and opening up new catalytic routes. Here, we establish the reaction coordinate space for guaiacol-H 2 reactions on Ru(0 0 0 1) using density functional theory and describe the atomic-scale effect of a polar protic solvent, H 2 O. Coupled H + and H* attack leads to quasi-equilibrated enol and keto intermediates as the precursors for C–O cleavage and ring saturation, respectively. For C–O cleavage, H 2 O solvent enables a lower energy pathway via concomitant transfer of the hydroxyl H + to the methoxy oxygen during C–OCH 3 cleavage, forming a charge separated [Ru(s)–(C 6 H 5 O – )…(H + )…OCH 3 ] transition state and reducing the barrier by up to 0.8 eV as compared to unassisted C–OCH 3 cleavage. For ring saturation, H* attack onto an unsaturated meta carbon is rate limiting with no direct solvent participation, suggesting that protic polar solvents selectively promote the C–O cleavage pathway. Taken together, we show that activating guaiacol for either C–O bond cleavage or ring saturation product formation depends on the reactive hydrogen identity (H* or H + ), enol/keto isomerization equilibrium, and accessibility of the proton assisted Car–OCH 3 cleavage transition state. Here, all such factors are tunable via changes to the solvent or metal identity.

09 BIOMASS FUELS↗
DOE OSTI2021

Effects of water on the kinetics of acetone hydrogenation over Pt and Ru catalysts

We employed an approach combining reaction kinetics measurements at steady state conditions, electronic structure calculations employing density functional theory, and microkinetic modeling for acetone hydrogenation to provide insights into the effects of water on metal catalyst surfaces for the hydrogenation of oxygenates over a wide range of reaction conditions. Elucidation of the repulsive interactions due to adsorbed water molecules at various reaction conditions provides a basis to formulate rate expressions for heterogeneous catalytic processes of biomass oxygenates. Reaction kinetics experiments were carried out at partial pressures of H 2 , acetone, water and helium in the range of 0.51-0.79, 0.02-0.13, 0.08-0.23, 0-0.28 atm, respectively. We show that the addition of water enhances the hydrogenation rate at 353 K and 1 atm on oxophilic metal catalysts such as Ru/C, whereas the same promotional effect of water is not observed for Pt-based catalysts. Microkinetic model predictions for the hydrogenation of acetone on Ru in the absence and presence of water, using enthalpies and entropies obtained from DFT calculations, were in agreement with the experimentally observed reaction orders and activation barriers. The model shows that a water-assisted hydroxypropyl path is expected to be the favored path on Ru with a rate-determining step of H-OH-mediated hydrogenation of C 3 H 6 OH (i.e., the hydroxypropyl intermediate formed by H 2 O-mediated initial hydrogenation of acetone) to produce isopropyl alcohol (IPA). Furthermore, hydrogen, acetone, hydroxypropyl intermediate and hydroxyl species were predicted to be abundant on the Ru surface with a high coverage of nearly 85%. Lastly, the combined studies of computational and experimental catalysis on hydrogenation reactions help to elucidate the mechanistic role of water on metal catalyzed reactions for producing chemical building blocks from biomass-derived oxygenates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Protonation Dynamics of Confined Ethanol–Water Mixtures in H-ZSM-5 from Machine Learning-Driven Metadynamics

Zeolites are indispensable heterogeneous catalysts in industrial chemical processes, valued for their strong Brønsted acidity, well-defined microporous frameworks, and tunable pore structures. Their catalytic activity arises primarily from Brønsted acid sites (BAS), typically present as bridging hydroxyl groups (Si–OH–Al). Under aqueous reaction conditions, these protons interact dynamically with water and alcohol molecules, leading to complex solvation and protonation behavior within confined pores. In this study, we investigate the protonation equilibrium occurring between ethanol and water at the BAS of acidic zeolites under varying hydration levels, i.e., C2H5OH–(H2O)n, n=1–4. Local structure was analyzed through an adaptive-learning global optimization algorithm, while enhanced sampling molecular dynamics simulations with Well-Tempered Metadynamics (WMetaD) and machine learning interatomic potentials (MLPs) provide free-energy surfaces (FES) at variable hydration levels. The results reveal a strong dependence of proton localization on the degree of hydration. At low hydration (1 water molecule), the proton resides predominantly on ethanol; with 2 water molecules, it shifts toward water, and at higher hydration (3 or more water molecules), it becomes extensively delocalized over the water cluster. These findings underscore the critical role of solvation in modulating acid site behavior and suggest that a minimum of three water molecules is necessary to fully stabilize the proton on water within the zeolite framework. This solvation threshold has significant implications for catalytic processes, particularly in biomass conversion reactions where alcohol protonation is a key step in dehydration mechanisms.

machine learning↗
DOE OSTI2024

Variations in proton transfer pathways and energetics on pristine and defect-rich quartz surfaces in water: Insights into the bimodal acidities of quartz

Hypothesis. Understanding the mechanisms of proton transfer on quartz surfaces in water is critical for a range of processes in geochemical, environmental, and materials sciences. The wide range of surface acidities (>9 pKa units) found on the ubiquitous mineral quartz is caused by the structural variations of surface silanol groups. Molecular scale simulations provide essential tools for elucidating the origin of site-specific surface acidities. Simulations. Here, we used density-functional tight-binding-based molecular dynamics combined with rare-event metadynamics simulations to probe the mechanisms of deprotonation reactions from ten representative surface silanol groups found on both pristine and defect-rich quartz (1 0 1) surfaces with Si vacancies. Findings. The results show that deprotonation is a highly dynamic process where both the surface hydroxyls and bridging oxygen atoms serve as the proton acceptors, in addition to water. Deprotonation of embedded silanols through intrasurface proton transfer exhibited lower pKa values with less H-bond participation and higher energy barriers, suggesting a new mechanism to explain the bimodal acidity observed on quartz surface. Defect sites, recently shown to comprise a significant portion of the quartz (1 0 1) surface, diversify the coordination and local H-bonding environments of the surface silanols, changing both the deprotonation pathways and energetics, leading to a wider range of pKa values (2.4 to 11.5) than that observed on pristine quartz surface (10.4 and 12.1).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗