Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “potential energy surface”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Quantitative Measurement of Positive and Negative Ion Species Ejected from a Li–O–H Surface by Hydrogen and Noble Gas Ion Irradiation

We report sputtering yields of Li + , H - , O - , and OH x - ion species from an Li–O–H surface for H, D, He, Ne, and Ar ion irradiation at 45° incidence in the energy range of 30–2000 eV. A Li film was deposited on a stainless steel target using Li evaporators in the LTX-β vessel, using the LTX-β Sample Exposure Probe (SEP), which includes an ultrahigh vacuum suitcase for transferring targets without significant contamination from air exposure. The SEP was used to transfer the Li-coated target from LTX-β to a separate Sample Exposure Station (SES) to perform ion exposure measurements. The SEP was also used for characterization of the Li-coated target utilizing X-ray photoelectron spectroscopy in a different chamber, showing that the lithium film surface was oxidized. Ion exposures were performed using an electron cyclotron resonance plasma source in the SES. Sputtered/ejected species were sampled by a quadrupole mass spectrometer with capabilities for detecting positive and negative ions, and an energy filter for determining the mean kinetic energy of the ejected ion species. All ion irradiations caused Li + ions to be ejected, while causing impurity ions such as H + , H - , O - and OH - to be ejected. Measured ion energies of Li + ions from a Li–O–H surface suggested that the typical sheath potential on the divertor surface can trap sputtered Li + ions, which were previously reported as ~ 60% of total sputtered Li species from Li targets (Allain and Ruzic in Nucl Fusion 42:202, 2002). Hence, our results for the sputtering yields of ejected ion species and their associated ion energies from a Li–O–H surface indicates that lithium sputtering is suppressed and impurity removal is enhanced due to the sheath potential at the divertor surface for fusion reactor applications.

Plasma-material interaction↗
DOE OSTI2020

The Land Sparing, Water Surface Use Efficiency, and Water Surface Transformation of Floating Photovoltaic Solar Energy Installations

Floating photovoltaic solar energy installations (FPVs) represent a new type of water surface use, potentially sparing land needed for agriculture and conservation. However, standardized metrics for the land sparing and resource use efficiencies of FPVs are absent. These metrics are critical to understanding the environmental and ecological impacts that FPVs may potentially exhibit. Here, we compared techno-hydrological and spatial attributes of four FPVs spanning different climatic regimes. Next, we defined and quantified the land sparing and water surface use efficiency (WSUE) of each FPV. Lastly, we coined and calculated the water surface transformation (WST) using generation data at the world’s first FPV (Far Niente Winery, California). The four FPVs spare 59,555 m2 of land and have a mean land sparing ratio of 2.7:1 m2 compared to ground-mounted PVs. Mean direct and total capacity-based WSUE is 94.5 ± 20.1 SD Wm−2 and 35.2 ± 27.4 SD Wm−2, respectively. Direct and total generation-based WST at Far Niente is 9.3 and 13.4 m2 MWh−1 yr−1, respectively; 2.3 times less area than ground-mounted utility-scale PVs. Our results reveal diverse techno-hydrological and spatial attributes of FPVs, the capacity of FPVs to spare land, and the utility of WSUE and WST metrics.

14 SOLAR ENERGY↗
DOE OSTI2026

A bioinspired approach for adaptive solid-solid phase change material coatings with optimized surface features for passive thermal regulation

The necessity to reduce global energy consumption calls for innovative strategies in building thermal management. Passive thermal regulation, particularly through bio-inspired designs, offers a promising avenue by mimicking nature's efficient control of optical properties. This research introduces a novel, climate-responsive coating that integrates optimized bio-inspired surface features with a solid-solid phase change material (SS-PCM) to dynamically manage solar absorptivity without adding additional thickness, enabling both heating and cooling as needed. Drawing on the photonic architectures of the Saharan silver ant and Morpho Didius butterfly, we employed a modeling and multi-objective optimization framework to tailor these surface features. Simulations reveal that surface texture, rather than the intrinsic phase transition of the SS PCM, dominates optical control. Relative to a flat SS PCM coating, optimized isotropic random roughness and broader range features yielded the highest passive heating power increase of about 144 % and 319 % respectively suitable for cold climates. Saharan ant-inspired features enhanced passive cooling for hot climates, achieving a 21.8 % improvement. For moderate climates, Butterfly-wing-inspired surface features provided a balanced enhancement of 19 % for heating and 7 % for cooling. Across all cases, the optimized surface features reduced combined heating and cooling energy demand more effectively than the baseline coating, while preserving material thickness. These findings demonstrate that climate-adaptive, optimized bio-inspired surface features can unlock the full potential of SS PCM coatings, providing a versatile pathway to significant energy savings in buildings and other applications. The methodology establishes a framework for designing next-generation adaptive envelopes that leverage natural photonic principles for high-impact, low-cost thermal regulation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Finite element procedure for thermomechanical and structural integrity analysis of beam intercepting devices subjected to free electron laser

High-energy particles, including photons (x-ray, γ-ray, bremsstrahlung), electrons, and protons, possess the capability to penetrate materials and deposit energy within them. The degree of absorption depends on both the energy and type of particles, as well as the properties of the materials with which they interact. This energy deposition can manifest either at the material's surface or throughout its volume, potentially resulting in various failure modes. The primary aim of this paper is to establish a structured analysis methodology for evaluating the structural integrity of beam-intercepting devices when subjected to high-energy particles. Here, the paper also reviews some of the underlying physics, pertinent to the scope of the thermomechanical analysis, potential failure modes, and introduces verification and validation methodologies. Engineers and researchers can utilize the guidelines presented in this paper to effectively plan the development of beam intercepting devices, thereby ensuring their reliability and performance in the presence of high-energy particle exposure.

Finite element analysis↗
DOE OSTI2022

Floating solar panels on reservoirs impact phytoplankton populations: A modelling experiment

Floating solar photovoltaic (FPV) deployments are increasing globally as the switch to renewable energy intensifies, representing a considerable water surface transformation. FPV installations can potentially impact aquatic ecosystem function, either positively or negatively. However, these impacts are poorly resolved given the challenges of collecting empirical data for field or modelling experiments. In particular, there is limited evidence on the response of phytoplankton to changes in water body thermal dynamics and light climate with FPV. Given the importance of understanding phytoplankton biomass and species composition for managing ecosystem services, we use an uncertainty estimation approach to simulate the effect of FPV coverage and array siting location on a UK reservoir. FPV coverage was modified in 10% increments from a baseline with 0% coverage to 100% coverage for three different FPV array siting locations based on reservoir circulation patterns. Results showed that FPV coverage significantly impacted thermal properties, resulting in highly variable impacts on phytoplankton biomass and species composition. The impacts on phytoplankton were often dependent on array siting location as well as surface coverage. Changes to phytoplankton species composition were offset by the decrease in phytoplankton biomass associated with increasing FPV coverage. We identified that similar phytoplankton biomass reductions could be achieved with less FPV coverage by deploying the FPV array on the water body's faster-flowing area than the central or slower flowing areas. The difference in response dependent on siting location could be used to tailor phytoplankton management in water bodies. Simulation of water body-FPV interactions efficiently using an uncertainty approach is an essential tool to rapidly develop understanding and ultimately inform FPV developers and water body managers looking to minimise negative impacts and maximise co-benefits.

14 SOLAR ENERGY↗
DOE OSTI2023

Investigating Spillover Energy as a Descriptor for Single-Atom Alloy Catalyst Design

The identification of thermodynamic descriptors of catalytic performance is essential for the rational design of heterogeneous catalysts. Here, we investigate how spillover energy, a descriptor quantifying whether intermediates are more stable at the dopant or host metal sites, can be used to design single-atom alloys (SAAs) for formic acid dehydrogenation. Further, using theoretical calculations, we identify NiCu as a SAA with favorable spillover energy and demonstrate that formate intermediates produced after the initial O–H activation are more stable at Ni sites where rate-determining C–H activation occurs. Surface science experiments demonstrated that NiCu(111) SAAs are more reactive than Cu(111) while they still follow the formate reaction pathway. However, reactor studies of silica-supported NiCu SAA nanoparticles showed only a modest improvement over Cu resulting from surface coverage effects. Overall, this study demonstrates the potential of engineering SAAs using spillover energy as a design parameter and highlights the importance of adsorbate–adsorbate interactions under steady-state operation.

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

Ternary PtIrNi Catalysts for Efficient Electrochemical Ammonia Oxidation

Due to high energy density, carbon-free feature, and easiness to be liquefied, ammonia (NH 3 ) has proved to be an effective alternative to hydrogen in low temperature fuel cells via its direct the ammonia oxidation reaction (AOR) for electricity generation. However, the kinetically sluggish AOR has prohibitively hindered the attractive direct ammonia fuel cell (DAFC) applications. Here we report an efficient AOR catalyst, in which ternary PtIrNi alloy nanoparticles well dispersed on a binary composite support consisting of porous silicon dioxide (SiO 2 ) and carboxyl-functionalized carbon nanotube (PtIrNi/SiO 2 -CNT-COOH) through a sonochemical-assisted synthesis strategy. The PtIrNi alloy nanoparticles, with the aid of abundant OHad provided by porous SiO 2 , and the improved electrical conductivity by CNTs, exhibit remarkable catalytic activity for the AOR in alkaline media. It is evidenced by a lower onset potential (~0.40 V vs. RHE) at room temperature, than that of commercial PtIr/C (ca. 0.43 V vs. RHE). Increasing NH 3 concentrations and operation temperatures significantly improve AOR performance of this catalyst. Specifically, AOR activity of the optimal PtIrNi nanoparticle catalyst can be significantly enhanced by elevating the temperature to 80ºC, with a much lower onset potential (~0.32 V vs. RHE), indicating that DAFC can be operated at higher temperature for increased performance. Constant-potential density functional theory (DFT) calculations showed that the Pt-Ir ensembles on {100}-terminated surfaces serve as the active site. Importantly, the introduction of Ni raises the center energy of the density of states projected onto the group d-orbitals of surface sites and thus lowers the theoretical onset potential for *NH 2 dehydrogenation to *NH when compared to Pt and Pt 3 Ir alloy.

25 ENERGY STORAGE↗
DOE OSTI2023

Excitation of surface waves in 3D ion beam neutralization

Neutralization of beams with 2D and 3D geometries by the electrons emitted from an external source is studied using particle-in-cell simulations. Our work reveals that the high-energy electrons excite Trivelpiece–Gould (TG) surface waves in the beams with 3D axisymmetric geometries. These high-energy electrons are generated because of a large amplitude electrostatic solitary wave (ESW) that forms near the electron source and has an electric potential amplitude more than three times the electron thermal energy. Here, we also find that surface wave excitation only happens when the beam radius is large enough at the ion source to attract enough electrons that could form the large amplitude ESW. A comparison of the 3D TG surface wave dispersion relation with an expression for 2D surface waves reveals that they become excited in 3D axisymmetric but not in 2D planar beam because of a higher phase speed requirement in the latter case.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Microkinetic Modeling: A Tool for Rational Catalyst Design

The design of heterogeneous catalysts relies on understanding the fundamental surface kinetics that controls catalyst performance, and microkinetic modeling is a tool that can help the researcher in streamlining the process of catalyst design. Microkinetic modeling is used to identify critical reaction intermediates and rate-determining elementary reactions, thereby providing vital information for designing an improved catalyst. In this review, we summarize general procedures for developing microkinetic models using reaction kinetics parameters obtained from experimental data, theoretical correlations, and quantum chemical calculations. We examine the methods required to ensure the thermodynamic consistency of the microkinetic model. We describe procedures required for parameter adjustments to account for the heterogeneity of the catalyst and the inherent errors in parameter estimation. We discuss the analysis of microkinetic models to determine the rate determining reactions using the degree of rate control and reversibility of each elementary reaction. We introduce incorporation of Brønsted–Evans–Polanyi relations and scaling relations in microkinetic models and the effects of these relations on catalytic performance and formation of volcano curves. We review the analysis of reaction schemes in terms of the maximum rate of elementary reactions, and we outline a procedure to identify kinetically significant transition states and adsorbed intermediates. Here, we explore the application of generalized rate expressions for the prediction of optimal binding energies of important surface intermediates and to estimate the extent of potential rate improvement. We also explore the application of microkinetic modeling in homogeneous catalysis, electro-catalysis, and transient reaction kinetics. We conclude by highlighting the challenges and opportunities in the application of microkinetic modeling for catalyst design.

42 ENGINEERING↗
DOE OSTI2023

Effects of Lower Troposphere Vertical Mixing on Simulated Clouds and Precipitation Over the Amazon During the Wet Season

Planetary boundary layer (PBL) schemes parameterize unresolved turbulent mixing within the PBL and free troposphere (FT). Previous studies reported that precipitation simulation over the Amazon in South America is quite sensitive to PBL schemes and the exact relationship between the turbulent mixing and precipitation processes is, however, not disentangled. In this study, regional climate simulations over the Amazon in January–February 2019 are examined at process level to understand the precipitation sensitivity to PBL scheme. The focus is on two PBL schemes, the Yonsei University (YSU) scheme, and the asymmetric convective model v2 (ACM2) scheme, which show the largest difference in the simulated precipitation. During daytime, while the FT clouds simulated by YSU dissipate, clouds simulated by ACM2 maintain because of enhanced moisture supply due to the enhanced vertical moisture relay transport process: (a) vertical mixing within PBL transports surface moisture to the PBL top, and (b) FT mixing feeds the moisture into the FT cloud deck. Due to the thick cloud deck over Amazon simulated by ACM2, surface radiative heating is reduced and consequently the convective available potential energy is reduced. As a result, precipitation is weaker from ACM2. Two key parameters dictating the vertical mixing are identified, p , an exponent determining boundary layer mixing and λ , a scale dictating FT mixing. Sensitivity simulations with altered p , λ , and other treatments within YSU and ACM2 confirm the precipitation sensitivity. The FT mixing in the presence of clouds appears most critical to explain the sensitivity between YSU and ACM2.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

High-resolution CH stretch spectroscopy of jet-cooled cyclopentyl radical: First insights into equilibrium structure, out-of-plane puckering, and IVR dynamics

First, high-resolution sub-Doppler infrared spectroscopic results for cyclopentyl radical (C 5 H 9 ) are reported on the α-CH stretch fundamental with suppression of spectral congestion achieved by adiabatic cooling to T rot ≈ 19(4) K in a slit jet expansion. Surprisingly, cyclopentyl radical exhibits a rotationally assignable infrared spectrum, despite 3N – 6 = 36 vibrational modes and an upper vibrational state density (ρ ≈ 40–90 #/cm –1 ) in the critical regime (ρ ≈ 100 #/cm –1 ) necessary for onset of intramolecular vibrational relaxation (IVR) dynamics. Such high-resolution data for cyclopentyl radical permit detailed fits to a rigid-rotor asymmetric top Hamiltonian, initial structural information for ground and vibrationally excited states, and opportunities for detailed comparison with theoretical predictions. Specifically, high level ab initio calculations at the coupled-cluster singles, doubles, and perturbative triples (CCSD(T))/ANO0, 1 level are used to calculate an out-of-plane bending potential, which reveals a C 2 symmetry double minimum 1D energy surface over a C 2v transition state. Here, the inversion barrier [V barrier ≈ 3.7(1) kcal/mol] is much larger than the effective moment of inertia for out-of-plane bending, resulting in localization of the cyclopentyl wavefunction near its C 2 symmetry equilibrium geometry and tunneling splittings for the ground state too small (<1 MHz) to be resolved under sub-Doppler slit jet conditions. The persistence of fully resolved high-resolution infrared spectroscopy for such large cyclic polyatomic radicals at high vibrational state densities suggests a “deceleration” of IVR for a cycloalkane ring topology, much as low frequency torsion/methyl rotation degrees of freedom have demonstrated a corresponding “acceleration” of IVR processes in linear hydrocarbons.

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

Predicting Catalyst Surface Stability Under Reaction Conditions Using Deep Reinforcement Learning and Machine Learning Potentials

Catalysts are critical for most large-scale energy intensive chemical transformation processes, such as energy storage, liquid fuel production, and the formation of chemical building blocks. The catalyst composition, structure, and morphology impact the performance under reaction conditions and influences properties like activity and selectivity. Many industrial catalysts offer imperfect activity/selectivity, or contain expensive metals. Further, catalysts can deactivate over time as the morphology changes or harsh reactive environments alter the surface structure. Methods to automatically model and predict the kinetics of how catalyst surfaces will restructure would enable engineers to design around these challenges, improve performance, increase longevity. This project investigated a specific type of machine learning model, deep reinforcement learning, machine learning models to act as surrogates for the physical system, and compared the results of these approaches with a specialized high-throughput experimental synthesis and measurement platform.

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

Realizing a topological diode effect on the surface of a topological Kondo insulator

Introducing the concept of topology into material science has sparked a revolution from classic electronic and optoelectronic devices to topological quantum devices. The latter has potential for transferring energy and information with unprecedented efficiency. Here, we demonstrate a topological diode effect on the surface of a three-dimensional material, SmB 6 , a candidate topological Kondo insulator. The diode effect is evidenced by pronounced rectification and photogalvanic effects under electromagnetic modulation and radiation at radio frequency. Our experimental results and modeling suggest that these prominent effects are intimately tied to the spatially inhomogeneous formation of topological surface states (TSS) at the intermediate temperature. This work provides a manner of breaking the mirror symmetry (in addition to the inversion symmetry), resulting in the formation of pn-junctions between puddles of metallic TSS. Further, this effect paves the way for efficient current rectifiers or energy-harvesting devices working down to radio frequency range at low temperature, which could be extended to high temperatures using other topological insulators with large bulk gap.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Cobalt-free and high-rate stable 5V lithium nickel manganese oxide spinel cathodes enabled via surface oxygen vacancies

Spinel LiNi 0.5 Mn 1.5 O 4 offers both the high-rate, low-cost and safety advantages of LiFePO 4 and the high energy density of LiNiₓMnᵧCo₁₋ₓ₋ᵧO₂ and LiNiₓCoᵧAlzO₂ cathodes. However, the large operating voltage of these materials induces electrolyte oxidation, which degrades the interface and drives Mn dissolution. These reactions are further exacerbated at high rates due to temperature rise. In this study, we discover that ammoniacal treatment followed by annealing introduces a high density of oxygen vacancies in the “near-surface region” of LiNi 0.5 Mn 1.5 O 4 particles. These vacancies release electrons changing the oxidation state of Mn and suppressing its tendency to oxidize the electrolyte. Further, these vacancies enhance the electrode’s electronic conductivity (by ∼3-fold) and Li + diffusivity (by ∼2-fold) greatly improving charge transport, especially when operated at high rates. This results in an across-the-board improvement in self-discharge, specific capacity, energy density, rate capability, coulombic efficiency and cycling stability. When cycled at ∼200 mA g −1 , the capacity fade averaged over 3000 cycles for the surface vacancy-enriched material is ∼0.0167% per cycle compared to an order of magnitude higher fade rate for the baseline material. In conclusion, these findings reveal the potential of targeted surface oxygen vacancy doping to develop cobalt-free and high energy density cathodes that tolerate fast charging and deliver improved cycle life.

Cobalt-free cathodes↗
DOE OSTI2020

Oxidation State and Surface Reconstruction of Cu under CO 2 Reduction Conditions from In Situ X-ray Characterization

The electrochemical CO 2 reduction reaction (CO 2 RR) using Cu-based catalysts holds great potential for producing valuable multi-carbon products from renewable energy. However, the chemical and structural state of Cu catalyst surfaces during the CO 2 RR remains a matter of debate. In this paper, we show the structural evolution of the near-surface region of polycrystalline Cu electrodes under in situ conditions through a combination of grazing incidence X-ray absorption spectroscopy (GIXAS) and X-ray diffraction (GIXRD). The in situ GIXAS reveals that the surface oxide layer is fully reduced to metallic Cu before the onset potential for CO 2 RR, and the catalyst maintains the metallic state across the potentials relevant to the CO 2 RR. We also find a preferential surface reconstruction of the polycrystalline Cu surface toward (100) facets in the presence of CO 2 . Quantitative analysis of the reconstruction profiles reveals that the degree of reconstruction increases with increasingly negative applied potentials, and it persists when the applied potential returns to more positive values. These findings show that the surface of Cu electrocatalysts is dynamic during the CO 2 RR, and emphasize the importance of in situ characterization to understand the surface structure and its role in electrocatalysis.

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

Strain-Mediated Sn Incorporation and Segregation in Compositionally Graded Ge 1– x Sn x Epilayers Grown by MBE at Different Temperatures

We investigated the process of Sn incorporation and surface segregation for compositionally graded Ge 1–x Sn x epilayers grown on high-quality Ge (001) substrates. Here, the growth resulted in pseudomorphic GeSn layers with a ~6% maximal Sn fraction at a constant substrate temperature. The maximal fraction of Sn was increased to 9.0% when the growth temperature was continuously lowered while increasing the Sn flux. The analysis of surface droplets and SIMS profiles of elemental composition give evidence of Sn rejection during the growth, potentially associated with a critical energy of elastic strain. The intentional reduction of the coherent strain by decreasing the Sn flux near the sample surface has been shown to trap a higher fraction of Sn in the Ge 1–x Sn x layer and lower surface segregation. Our results demonstrate that strain relief by misfit dislocations in the compositionally graded layer is inhibited, which leads to Sn segregation. Specifically, the compressive strain in the graded Ge 1–x Sn x epilayer is effectively “zero” near the interface with the Ge substrate and increases up to about –1.5 × 10 –2 near the surface. Thus, although the nucleation of a dislocation may reduce the compressive strain for the top region of the epilayer, it is not beneficial for the bottom region.

42 ENGINEERING↗
DOE OSTI2021

Enabling bifacial thin film devices by developing a back surface field using Cu x AlO y

Bifacial solar cells have the potential to increase the energy yield per unit area over traditional monofacial devices without significant added cost, driving $/kWh costs lower and accelerating the adoption of solar photovoltaics. However, the performance of bifacial thin film solar cells significantly lags that achieved by crystalline silicon cells. Here we incorporate wide bandgap Cu x AlO y as a back buffer layer for CdTe devices and achieve a backside illuminated device with high current density and high fill factor. Moreover, these values remain nearly constant even as the absorber layer thickness changes, indicating that a fully-depleted device is not required for efficient charge collection. We show that this response is indicative of a back surface field, albeit with a persistent high back surface recombination velocity. By managing electron reflection, we achieved a backside illumination conversion efficiency of 7.1% and bifaciality of 0.55 for a 3.3 µm CdTe device and 8.0% and 0.62 for a 2 µm device. Future improvements can be made by identifying and incorporating a passivation material that reduces the back surface recombination velocity.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A Deep Potential model for liquid–vapor equilibrium and cavitation rates of water

Computational studies of liquid water and its phase transition into vapor have traditionally been performed using classical water models. Here, we utilize the Deep Potential methodology—a machine learning approach—to study this ubiquitous phase transition, starting from the phase diagram in the liquid–vapor coexistence regime. The machine learning model is trained on ab initio energies and forces based on the SCAN density functional, which has been previously shown to reproduce solid phases and other properties of water. Here, we compute the surface tension, saturation pressure, and enthalpy of vaporization for a range of temperatures spanning from 300 to 600 K and evaluate the Deep Potential model performance against experimental results and the semiempirical TIP4P/2005 classical model. Moreover, by employing the seeding technique, we evaluate the free energy barrier and nucleation rate at negative pressures for the isotherm of 296.4 K. Further, we find that the nucleation rates obtained from the Deep Potential model deviate from those computed for the TIP4P/2005 water model due to an underestimation in the surface tension from the Deep Potential model. From analysis of the seeding simulations, we also evaluate the Tolman length for the Deep Potential water model, which is (0.091 ± 0.008) nm at 296.4 K. Finally, we identify that water molecules display a preferential orientation in the liquid–vapor interface, in which H atoms tend to point toward the vapor phase to maximize the enthalpic gain of interfacial molecules. We find that this behavior is more pronounced for planar interfaces than for the curved interfaces in bubbles. This work represents the first application of Deep Potential models to the study of liquid–vapor coexistence and water cavitation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗