Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “water cluster”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

VUV Photoionization Dynamics and Reactivity of Heterogeneous Water Clusters

This feature focuses on bridging isolated water and bulk water studies. It assembles experiments and theory on water clusters, mainly probed by vacuum ultraviolet (VUV) radiation, and summarizes what the effects are of intermolecular interactions on both the spectroscopy and the VUV-induced processes in water. In particular, it highlights studies of heterogeneous water clusters─those incorporating other molecular species such as naphthalene, methane, formic acid, and glycerol─which serve as model systems to investigate the fundamental roles of water in hydrogen bonding networks, proton transfer, and astrochemical processes. These mixed clusters provide a platform to revisit the stability of protonated water clusters and compile observations of ionization-induced structural rearrangements and fragmentation, especially in systems involving hydroxyl-rich cosolvents. This feature then explores energy transfer mechanisms in molecular clusters, following VUV photoexcitation. It concludes with potential future directions. First, we investigate excited-state dynamics in molecular clusters through direct probing, complementing the synchrotron studies discussed here. Second, exploring gas-phase molecule evaporation from confined spaces and interfaces using advanced spectroscopic techniques sheds light on these ubiquitous, yet currently debated, molecular processes.

Cluster chemistry↗
DOE OSTI2025

Water Clustering Modulates Activity and Enables Hydrogenated Product Formation during Carbon Monoxide Electroreduction in Aprotic Media

Water solvation plays a critical role in a wide range of electrochemical transformations, but its role is often convoluted since water is typically used as both a solvent and a proton source. Here, we experimentally control water speciation and activity using aprotic solvent media during the carbon monoxide reduction reaction (CORR). Remarkably, we show that aprotic solvents that support microheterogeneous water–water clusters lead to significant amounts of CORR products (methane and ethylene) with a maximum ethylene Faradaic efficiency of 22% in acetonitrile (χ H 2 O = 0.2). In contrast, microhomogeneous systems–where water integrates into the solvents’ intermolecular binding network and has lower activity–primarily support the undesired hydrogen evolution reaction (HER). As a result, insights gained expand our understanding of water activity and nonaqueous electrolyte design for other important transformation reactions beyond CO reduction, such as CO2RR and HER.

Electrolytes↗
DOE OSTI2025

Theory of the photomolecular effect

It is well-known that water in liquid and vapor phases exhibits weak visible-light absorption. Recent experiments, however, show that at the liquid-air interface, absorption drastically increases, accelerating evaporation beyond thermal limits by 2–5 times. Strikingly, evaporation peaks at green wavelengths, despite no corresponding absorptance peak. The underlying mechanism of this observation, termed the photomolecular effect, remains puzzling, particularly as water molecules do not exhibit resonance peaks in the visible spectrum. Here, we present a theoretical model explaining the effect. We show that surface-bound water clusters undergo non-resonant photon-driven evaporation, with green light being particularly effective. Crucially, we do not expect green light to couple more strongly than other wavelengths to the molecules at the surface, rather the energy of green light is more effectively used to vaporize water. Our model accounts for why the evaporation peak does not align with absorptance and provides a quantum mechanical explanation of how a single photon can vaporize an entire molecular cluster. This model challenges conventional views on water-light interactions, revealing a fundamentally non-thermal evaporation mechanism. Furthermore, our findings have implications for water purification, energy-efficient drying, and climate modeling, opening new pathways for optimizing evaporation-based technologies.

Interfacial evaporation↗
DOE OSTI2025

Linear-scaling quadruple excitations in local pair natural orbital coupled-cluster theory

Here, we present a fast, asymptotically linear-scaling implementation of the perturbative quadruples energy correction in coupled-cluster theory using local natural orbitals. Our work follows the domain-based local pair natural orbital (DLPNO) approach previously applied to lower levels of excitations in coupled-cluster theory. Our DLPNO-CCSDT(Q) algorithm uses converged doubles and triples amplitudes from a preceding DLPNO-CCSDT computation to compute the quadruples amplitude and energy in the quadruples natural orbital (QNO) basis. We demonstrate the compactness of the QNO space, showing that more than 95% of the (Q) correction can be recovered using relatively loose natural orbital cutoffs, compared to the tighter cutoffs used in pair and triples natural orbitals at lower levels of coupled-cluster theory. We also highlight the accuracy of our algorithm in the computation of relative energies, which yields deviations of sub-kJ mol −1 in relative energy compared to the canonical CCSDT(Q). Timings are conducted on a series of growing linear alkanes (up to 10 carbons and 608 basis functions) and water clusters (up to 49 water molecules and 2842 basis functions) to establish the asymptotic linear-scaling of our DLPNO-(Q) algorithm.

Auxiliary functions↗
DOE OSTI2025

Linear-Scaling Local Natural Orbital-Based Full Triples Treatment in Coupled-Cluster Theory

We present an efficient, asymptotically linear-scaling implementation of the canonically O(N 8 ) coupled-cluster method with singles, doubles, and full triples excitations (CCSDT) method. We apply the domain-based local pair natural orbital (DLPNO) approach for computing CCSDT amplitudes. Our method, called DLPNO–CCSDT, uses the converged coupled-cluster amplitudes from a preceding DLPNO–CCSD(T) computation as a starting point for the solution of the CCSDT equations in the local natural orbital basis. To simplify the working equations, we t1-dress our two-electron integrals and Fock matrices, allowing our equations to take on the form of CCDT. With appropriate parameters, our method can recover more than 99.99% of the total canonical CCSDT correlation energy. In addition, we demonstrate that our method consistently yields sub-kJ mol –1 errors in relative energies when compared to canonical CCSDT, and, likewise, when computing the difference between CCSDT and CCSD(T). Finally, to highlight the low scaling of our algorithm, we present timings on linear alkanes (up to 30 carbons and 730 basis functions) and water clusters (up to 131 water molecules and 3144 basis functions).

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 OSTI2025

Microcarbonation of Naphthalene: An Experimental and Computational Study of Photoionization in Naphthalene-Carbon Dioxide Clusters

The photoionization of naphthalene (N)-carbon dioxide (CO 2 ) clusters was studied using tunable vacuum ultraviolet (VUV) radiation from a synchrotron in the photon range of 8.0 to 13.7 eV, in combination with time-of-flight mass spectrometry. Clusters of monomer, dimer, and trimer naphthalene with CO 2 (N­(CO 2 ) 0–6 , N 2 ­(CO 2 ) 0–3 , N3) were observed. The lowest-energy conformers were obtained via a conformer search, followed by geometry optimizations at the ωB97X-V2/aug-cc-pVTZ (monomer) and ωB97X-V2/aug-cc-pVDZ (dimer) levels of theory. Carbon dioxide was found to preferentially cluster on top of the naphthalene molecule (in an out-of-plane configuration). From the mass spectra, photoionization intensity curves (PICs) were constructed, and appearance energies (AEs) were determined. No substantial trend in AE was observed with increasing size of the naphthalene-carbon dioxide clusters; rather, AE oscillations around the value for pure naphthalene were observed. These AE oscillations are also observed in a recently studied naphthalene-water cluster system, though in this system a slight downward trend in AE for pure naphthalene clusters (N1/2/3/4) was observed. The differences between the two systems are attributed to differing interaction strengths. In conclusion, understanding these differences may aid in determining how photoprocessing can proceed differently depending on the dominant matrix component in interstellar ices.

Wannenmacher, Anna [Lawrence Berkeley National Lab↗
DOE OSTI2025

Synergistic Density Functional Theory and Molecular Dynamics Approach to Elucidate PNIPAM–Water Interaction Mechanisms

This study employs Density Functional Theory (DFT) and Molecular Dynamics (MD) simulations to investigate interactions between water molecules and Poly(N-isopropylacrylamide) (PNIPAM). DFT reveals preferential water binding sites, with enhanced binding energy observed in the linker zone. Quantum Theory of Atoms in Molecules (QTAIM) and electron localization function (ELF) analyses highlight the roles of hydrogen bonding and steric hindrance. MD simulations unveil temperature-dependent hydration dynamics, with structural transitions marked by changes in the radius of gyration (Rg) and the radial distribution function (RDF), aligning with DFT findings. Our work goes beyond prior studies by combining a DFT, QTAIM and MD simulations approach across different PNIPAM monomer-to-30mer structures. It introduces a systematic quantification of pseudo-saturation thresholds and explores water clustering dynamics with structural specificity, which have not been previously reported in the literature. These novel insights establish a more complete molecular-level picture of PNIPAM hydration behavior and temperature responsiveness, emphasizing the importance of amide hydrogen and carbonyl oxygen sites in hydrogen bonding, which weakens above the lower critical solution temperature (LCST), resulting in increased hydrophobicity and paving the way for understanding water sorption mechanisms, offering guidance for future applications such as dehumidification and atmospheric water harvesting.

08 HYDROGEN↗
DOE OSTI2025

Supercooling of Alaskan Beetle Larvae as a Winter Survival Strategy

Insects are able to survive subfreezing temperatures by either limiting ice crystal formation in their bodies or through freeze avoidance. Beetle larvae are able to avoid freezing in winter by dehydrating in the fall months and replacing their body water content with high concentrations of glycerol. This enables the body fluid of the insect to supercool, and even vitrify, recovering unharmed when the temperature warms in the spring. Using nondestructive, high‐energy X‐ray synchrotron diffraction experiments, direct insight into how cryopreservation occurs at the atomic level within the beetle larvae has been obtained. The results shed light on the molecular‐level interactions associated with the mechanism responsible for surviving freezing temperatures. The molecular models of severely dehydrated Alaskan beetle larvae, based on glycerol‐water mixtures, yield a total of 4.2 ± 1.2 intermolecular hydrogen bonds per glycerol molecule at 275 K, in good agreement with existing molecular dynamics simulations. Most importantly, they show that if just over half the body fluid content is water, the water clusters are too small to form ice crystals that cause cellular damage.

Benmore, Chris J. [X‐Ray Science Division Advanced↗
DOE OSTI2026

Using Density-Corrected DFT to Understand Density-Driven and Functional-Dependent Errors in Ab Initio Simulations of the Hydrated Electron

The hydrated electron, an excess electron in liquid water, plays a crucial role in a plethora of chemical processes, motivating extensive research efforts to characterize its structure, dynamics, and reactivity in solution. Recent theoretical approaches to understanding this intriguing object have involved ab initio simulations based on density functional theory (DFT). Although DFT allows for the study of hydrated electron reactivity and quantum mechanical behavior, it is well-known that anionic systems can suffer from significant density-driven errors (DDEs). Density-corrected DFT (DC-DFT) provides a framework to mitigate such errors; the method reduces DDEs by replacing the self-consistent (SC) density associated with a given density functional with the Hartree–Fock (HF) density. Since HF densities tend to be more localized than DFT SC densities, the DC-DFT scheme significantly improves errors in calculations where the SC density is spuriously delocalized. Here, we investigate how the use of density correction affects the calculated properties of the DFT-simulated (PBEh) hydrated electron, a particularly challenging diffuse anionic system to simulate. First, we analyze charge delocalization in a system consisting of a model octahedral hydrated electron water cluster (the so-called Kevan structure) along with a spatially separated sulfur atom. We show that the use of density correction indeed reduces DDEs in comparison to a standard DFT global hybrid functional. We then propagate molecular dynamics trajectories of the hydrated electron using DC-DFT, where we find that DC further localizes electron density in the cavity region, a signature of reduced charge delocalization. Unfortunately, the decreased radius of gyration of the spin density and corresponding tightening of the local solvation structure from density correction causes predicted observables to deviate further from experimental measurements than when density correction is not employed. Here, we argue that DC’s worse agreement with experiment results from the removal of a fortuitous cancellation of errors that is intrinsic to the PBEh functional. This indicates that the difficulties with DFT to simulate hydrated electrons are primarily due to the inherent approximations in DFT rather than to density-driven errors.

Density functional theory↗
DOE OSTI2025

Untangling Sources of Error in the Density-Functional Many-Body Expansion

The many-body expansion provides a framework for data-driven applications of electronic structure theory, including parametrization of classical force fields and machine learning. In this article, we demonstrate that its use significantly amplifies quadrature grid errors when modern density-functional approximations are employed. Standard grids that work well in conventional density-functional calculations result in runaway error accumulation when used with the many-body expansion. At the same time, delocalization error is also exacerbated, leading to exaggerated estimates of nonadditive n-body interactions. This is illustrated for anion–water clusters using the SCAN, r2SCAN, ωB97X-V and ωB97M-V functionals. By employing dense quadrature grids, the inherent self-interaction error is exposed, which can then be mitigated using a variety of other strategies.

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

Dehydration of Methyl Lactate on Alkali Cation-Exchanged Faujasite: Effects of Metal Cation Identity and Water Pressure

Turnover rates for the catalytic dehydration of methyl lactate (ML) over ion-exchanged faujasite (FAU) catalysts depend on the identity of alkali metal cations (Na + , K + , Cs + ) and local solvation effects. Analysis of rate measurements and in situ infrared spectroscopy gives evidence that the reaction involves kinetically relevant dissociation of adsorbed methyl lactate upon alkali metal cations. This process involves concerted methyl transfer to the surface and dissociation of the alkali metal from the framework, which occurs at cationic active sites that remain predominantly unoccupied under relevant conditions (0.5–10 kPa ML, 0.5–15 kPa H 2 O, 563–583 K). Despite the mechanistic similarities, apparent activation enthalpies (ΔH app ‡ ) decrease linearly (47 kJ mol –1 from Na + to Cs + ) with ionization energy and cationic radius, and apparent activation entropies (ΔS app ‡ ) decrease 74 J mol –1 K –1 . These trends reflect electrostatic interactions that stabilize the cations to the anionic sites on the zeolite: stronger association between these charges leads to increasingly endothermic processes to displace the alkali metal to form a cationic methoxy and an intrapore metal lactate intermediate. Water physisorption measurements suggest alkali metal ions bind superstoichiometric quantities of water within FAU pores, and in situ infrared spectra suggest the concerted adsorption of ML requires reorganization of this water. Consequently, these processes introduce entropic gains that partially offset entropy losses associated with ML adsorption. Hence, turnover rates differ only by a factor of 2 among Na-, K-, and Cs-FAU at 573 K (ΔΔG app ‡ = 5 kJ mol –1 ). These findings demonstrate the interplay of alkali metal ions with zeolite active sites and intrapore water clusters for ML dehydration, indicating that these interactions can be leveraged to deliver optimal performance under different reaction conditions.

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

A Transferable Force Field for Predicting Adsorption and Diffusion of Water in Cationic Zeolites with Coupled Cluster Accuracy

We present a transferable force field for water in proton-exchanged, alkali (Li, Na, K, Rb, and Cs) metal-exchanged, and alkaline-earth (Mg, Ca, Sr, and Ba) metal-exchanged zeolites. The fitting methodology is based on adsorbate–adsorbent interaction energies obtained from periodic density functional theory calculations and corrected using the coupled-cluster method applied to small model clusters. To ensure an accurate prediction of both adsorption and diffusion properties of water, sets of configurations that sample both adsorption sites and intracrystalline hopping transition states were used in the fitting. The quality of the force field is assessed for a wide range of zeolites with different topologies and chemical compositions, demonstrating good agreement between theoretical predictions and experimental measurements of water adsorption and diffusion.

Adsorption↗
DOE OSTI2026

Superstructure Optimization for Brine Valorization from Brackish Water Desalination

This poster presents preliminary results from a superstructure optimization framework developed to identify cost-optimal brine valorization configurations for brackish water desalination plants across diverse U.S. regional feed chemistries. The study uses brackish groundwater compositions from Arizona, California, Florida, New Mexico, and Texas. Using Pyomo Generalized Disjunctive Programming (GDP) within the WaterTAP modeling environment, the optimization framework simultaneously evaluates thousands of candidate treatment configurations, spanning nanofiltration, reverse osmosis, and chemical precipitation, to minimize the levelized cost of water (LCOW) while meeting water recovery targets and product recovery constraints. Results across eight representative feed clusters demonstrate water recovery rates of 57–87% and net LCOW values ranging from -$0.032/m³ (net revenue-positive) to $0.80/m. Notably, no single process configuration was optimal across all feed types, underscoring the necessity of feed-specific optimization. Products targeted include calcium carbonate (CaCO₃) at $0.01/kg and sodium chloride (NaCl) at $0.10/kg, both at 95% purity, with product revenues offsetting treatment costs in several scenarios. The work advances NAWI's process systems engineering capabilities for multi-configuration screening.

58 GEOSCIENCES↗
DOE OSTI2025

Nanoscopic Structure of the Interface between Reduced TiO 2– x (110) and Water Vapor

We investigate the interaction between a reduced rutile TiO 2–x (110)-(1×4) surface and water vapor at ambient pressures using atomic force microscopy and X-ray photoelectron spectroscopy (AP-AFM and AP-XPS). Our results reveal that water molecules strongly interact with the reduced surface, leading to hydroxylation and localized clustering of water molecules. In defect-rich regions, AFM tip-induced restructuring causes removal of the topmost surface layer, highlighting the lowered cohesive energy of the surface atoms upon hydroxylation. These findings provide new insights into the water adsorption and restructuring mechanisms on reducible transition metal oxides, relevant for catalytic and environmental applications.

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

An atomic cluster expansion (ACE) potential for water under extreme conditions

We present a machine learning interatomic potential for water designed to capture its complex multiphase behavior, including both molecular and superionic ice phases. The potential is based on the atomic cluster expansion (ACE) formulation and has been parameterized to enable high-fidelity molecular dynamics simulations of water under extreme conditions, for pressures up to 100 GPa and for temperatures between 500 and 6000 K. A diverse range of configurations was generated through ab initio molecular dynamics (AI-MD) simulations, covering insulating and superionic ice phases, liquid water, and dissociated plasma phase. We demonstrate that the H 2 O ACE potential accurately reproduces experimental and DFT predicted isotherms and Hugoniots. Crucially, the potential is able to capture the intricate phase behavior of water, including the transition from molecular fluid to the appropriate solid ice phases, and the superionic ice phases. This work provides a robust interatomic potential that can be used for large-scale, accurate simulations of water under extreme thermodynamic conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Photocatalytic Overall Water Splitting at the Integrated Rh–MoRhO x Cluster Heterostructure on InGaN/GaN Nanowires

The quest for efficient solar-driven water splitting, a promising avenue for clean fuel production, faces challenges due to limited solar energy conversion efficiency. Traditional approaches study the overall water splitting as two spatially separate half reactions on two unrelated sites, hindering full utilization of photogenerated charge and water molecules. To overcome these limitations, an integrated cluster heterostructure catalyst on InGaN/GaN semiconductor nanowires is proposed for the effective utilization of photogenerated charge carriers and water molecules on the same redox localization. By establishing the fast charge extraction kinetics based on InGaN/GaN nanowires, the integration of Rh and MoRhOx clusters on the nanowire surface enables simultaneous and fast hydrogen/oxygen evolution reactions at the cluster heterostructure. Furthermore, the integrated strategy can enhance the charge redistribution across the heterostructure between the two clusters, further optimizing adsorption of reaction intermediates on each cluster for boosted photocatalytic water splitting activity. Consequently, the integrated heterostructure triggers a 40-fold increased hydrogen production efficiency in an artificial leaf system. This study provides valuable insights for the rational design of advanced heterostructured photocatalysts for water splitting and beyond.

GaN nanowire↗