SEARCH · Engineering Papers

Results for “Interfacial tension”

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.

SourceResearch connectionFrom yearThrough year

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

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Synergistic Assembly of Charged Oligomers and Amino Acids at the Air–Water Interface: An Avenue toward Surface-Directed CO 2 Capture

Interfaces are considered a major bottleneck in the capture of CO 2 from air. Efforts to design surfaces to enhance CO 2 capture probabilities are challenging due to the remarkably poor understanding of chemistry and self-assembly taking place at these interfaces. In this report we leverage surface-specific vibrational spectroscopy, Langmuir trough techniques, and simulations to mechanistically elucidate how cationic oligomers can drive surface localization of amino acids (AAs) that serve as CO 2 capture agents speeding up the apparent rate of absorption. We demonstrate how tuning these interfaces provides a means to facilitate CO 2 capture chemistry to occur at the interface, while lowering surface tension and improving transport/reaction probabilities. We show that in the presence of interfacial AA-rich aggregates, one can improve capture probabilities vs that of a bare interface, which holds promise in addressing climate change through the removal of CO 2 via tailored interfaces and associated chemistries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Covalently-Bonded Laminar Assembly of Van der Waals Semiconductors with Polymers: Toward High-Performance Flexible Devices

Van der Waals semiconductors (vdWS) offer superior mechanical and electrical properties and are promising for flexible microelectronics when combined with polymer substrates. However, the self-passivated vdWS surfaces and their weak adhesion to polymers tend to cause interfacial sliding and wrinkling, and thus, are still challenging the reliability of vdWS-based flexible devices. Here, in this study, an effective covalent vdWS–polymer lamination method with high stretch tolerance and excellent electronic performance is reported. Using molybdenum disulfide (MoS 2 ) and polydimethylsiloxane (PDMS) as a case study, gold–chalcogen bonding and mercapto silane bridges are leveraged. The resulting composite structures exhibit more uniform and stronger interfacial adhesion. This enhanced coupling also enables the observation of a theoretically predicted tension-induced band structure transition in MoS 2 . Moreover, no obvious degradation in the devices’ structural and electrical properties is identified after numerous mechanical cycle tests. This high-quality lamination enhances the reliability of vdWS-based flexible microelectronics, accelerating their practical applications in biomedical research and consumer electronics.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Fluctuating hydrodynamics and the Rayleigh–Plateau instability

The Rayleigh–Plateau instability occurs when surface tension makes a fluid column become unstable to small perturbations. At nanometer scales, thermal fluctuations are comparable to interfacial energy densities. Consequently, at these scales, thermal fluctuations play a significant role in the dynamics of the instability. These microscopic effects have previously been investigated numerically using particle-based simulations, such as molecular dynamics (MD), and stochastic partial differential equation–based hydrodynamic models, such as stochastic lubrication theory. In this paper, we present an incompressible fluctuating hydrodynamics model with a diffuse-interface formulation for binary fluid mixtures designed for the study of stochastic interfacial phenomena. An efficient numerical algorithm is outlined and validated in numerical simulations of stable equilibrium interfaces. We present results from simulations of the Rayleigh–Plateau instability for long cylinders pinching into droplets for Ohnesorge numbers of Oh = 0.5 and 5.0. Both stochastic and perturbed deterministic simulations are analyzed and ensemble results show significant differences in the temporal evolution of the minimum radius near pinching. Short cylinders, with lengths less than their circumference, were also investigated. As previously observed in MD simulations, we find that thermal fluctuations cause these to pinch in cases where a perturbed cylinder would be stable deterministically. Finally, we show that the fluctuating hydrodynamics model can be applied to study a broader range of surface tension–driven phenomena.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Phase Separation, Capillarity, and Odd-Surface Flows in Chiral Active Matter

Active phase separations evade canonical thermodynamic descriptions and have thus challenged our understanding of coexistence and interfacial phenomena. Considerable progress has been made towards a nonequilibrium theoretical description of these traditionally thermodynamic concepts. Spatial parity symmetry is conspicuously assumed in much of this progress, despite the ubiquity of chirality in experimentally realized systems. Here, in this Letter, we derive a theory for the phase coexistence and interfacial fluctuations of a system that microscopically violates spatial parity. We find suppression of the phase separation as chirality is increased as well as the development of steady-state currents tangential to the interface dividing the phases. These odd flows are irrelevant to stationary interfacial properties, with stability, capillary fluctuations, and surface area minimization determined entirely by the capillary surface tension. Using large-scale Brownian dynamics simulations, we find excellent agreement with our theoretical scaling predictions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Surface enrichment dictates block copolymer orientation

Orientation of block copolymer (BCP) morphology in thin films is critical to applications as nanostructured coatings. Despite being well-studied, the ability to control BCP orientation across all possible block constituents remains challenging. Here, in this study, we deploy coarse-grained molecular dynamics simulations to study diblock copolymer ordering in thin films, focusing on chain makeup, substrate surface energy, and surface tension disparity between the two constituent blocks. We explore the multi-dimensional parameter space of ordering using a machine-learning approach, where an autonomous loop using a Gaussian process (GP) control algorithm iteratively selects high-value simulations to compute. The GP kernel was engineered to capture known symmetries. The trained GP model serves as both a complete map of system response, and a robust means of extracting material knowledge. We demonstrate that the vertical orientation of BCP phases depends on several counter-balancing energetic contributions, including entropic and enthalpic material enrichment at interfaces, distortion of morphological objects through the film depth, and of course interfacial energies. BCP lamellae are found more resistant to these effects, and thus more robustly form vertical orientations across a broad range of conditions; while BCP cylinders are found to be highly sensitive to surface tension disparity.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2025

Sharp front tracking with geometric interface reconstruction

Here, this paper presents a novel sharp front-tracking method designed to address limitations in classical front-tracking approaches, specifically their reliance on smooth interpolation kernels and extended stencils for coupling the front and fluid mesh. In contrast, the proposed method employs exclusively sharp, localized interpolation and spreading kernels, restricting the coupling to the interfacial fluid cells–those containing the interface/front. This localized coupling is achieved by integrating a divergence-preserving velocity interpolation method with a piecewise parabolic interface calculation (PPIC) and a polyhedron intersection algorithm to compute the indicator function and local interface curvature. Surface tension is computed using the Continuum Surface Force (CSF) method, maintaining consistency with the sharp representation. Additionally, we propose an efficient local roughness smoothing implementation to account for surface mesh undulations, which is easily applicable to any triangulated surface mesh. Building on our previous work, the primary innovation of this study lies in the localization of the coupling for both the indicator function and surface tension calculations. By reducing the interface thickness on the fluid mesh to a single cell, as opposed to the 4–5 cell spans typical in classical methods, the proposed sharp front-tracking method achieves a highly localized and accurate representation of the interface. This sharper representation mitigates parasitic currents and improves force balancing, making it particularly suitable for scenarios where the interface plays a critical role, such as microfluidics, fluid-fluid interactions, and fluid-structure interactions. The proposed method is comprehensively validated and tested on canonical interfacial flow problems, including stationary and translating Laplace equilibria, oscillating droplets, and rising bubbles. The presented results demonstrate that the sharp front-tracking method significantly outperforms the classical approach in terms of accuracy, stability, and computational efficiency. Notably, parasitic currents are reduced by approximately two orders of magnitude and stable results are obtained for parameter ranges where classical front tracking fails to converge.

42 ENGINEERING↗

DOE OSTI2025

Energy absorption of architectured PrintCast interpenetrating composites in tension

Additively manufactured (AM) metal-metal composites consisting of PrintCasted 316 L austenitic stainless-steel lattice structures infiltrated with A356 casting alloy, have recently been developed for use in high energy absorption systems with potential applications ranging from static load bearing to dynamic blast containment structures. This system has a unique mechanical behavior as the volume fraction of lattice increases showing a transition from localized to de-localized failure and dramatic increase in energy absorption capability. In this work, PrintCast A356/316 L composite tensile specimens were produced with lattice volume fractions ranging from 20 % to 50 % to capture the range of this behavior. Finite element simulations support neutron diffraction measurements of stress state. Results illustrate that in tension, the reinforcement material is in tension while the matrix support material is in compression, information offering significant insight into the transition to de-localized failure. Moreover, the simulation results provide further insight into how interfacial bonding (or lack of bonding) affects the energy absorption capabilities of the PrintCast composites.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Are Contact Angle Measurements Useful for Oxide-Coated Liquid Metals?

This work establishes that static contact angles for gallium-based liquid metals have no utility despite the continued and common use of such angles in the literature. In the presence of oxygen, these metals rapidly form a thin (~1–3 nm) surface oxide “skin” that adheres to many surfaces and mechanically impedes its flow. This property is problematic for contact angle measurements, which presume the ability of liquids to flow freely to adopt shapes that minimize the interfacial energy. We show here that advancing angles for a metal are always high (>140°)—even on substrates to which it adheres—because the solid native oxide must rupture in tension to advance the contact line. The advancing angle for the metal depends subtly on the substrate surface chemistry but does not vary strongly with hydrophobicity of the substrate. During receding measurements, the metal droplet initially sags as the liquid withdraws from the “sac” formed by the skin and thus the contact area with the substrate initially increases despite its volumetric recession. Furthermore, the oxide pins at the perimeter of the deflated “sac” on all the surfaces are tested, except for certain rough surfaces. With additional withdrawal of the liquid metal, the pinned angle gets smaller until eventually the oxide “sac” collapses. Thus, static contact angles can be manipulated mechanically from 0° to >140° due to hysteresis and are therefore uninformative. We also provide recommendations and best practices for wetting experiments, which may find use in applications that use these alloys such as soft electronics, composites, and microfluidics.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Probing air-water interfaces of dibutyl phosphoric acid (HDBP) aqueous solutions using vibrational sum frequency generation (vSFG) spectroscopy

Liquid-liquid extraction is a separation technique implemented in a wide variety of areas, achieving particular success in both the nuclear and biomedical fields. In this work, vibrational sum frequency generation spectroscopy (VSFG) and surface tension measurements were used to investigate the adsorption of dibutyl phosphate (DBP) at air-aqueous interfaces to simulate liquid-liquid systems relevant to the Plutonium Uranium Redox Extraction (PUREX) Process. The objective of this work is to establish qualitative relationships between changes in the bulk aqueous phase concentrations of DBP and its concentration and structure at air-liquid interface as probed with VSFG. Nitric acid concentration and solution ionic strength were varied to examine their effect on the interfacial DBP.. Introduction of DBP into neat water resulted in reduction of the VSFG spectral intensity in the dangling O-H region (3680 – 3800 cm -1 ) but large increase in the H-bonded O-H stretch frequency region (3000 – 3500 cm-1) and the appearance of the CH 3 symmetric stretch and CH 3 Fermi resonance peaks at ~ 2880 and 2945 cm -1 , respectively, indicating DBP at the air-water interface. The intensity of the C-H strecth peaks increased as DBP concentration increased from 0.24 to 32 mM, accompanied by a decreasing surface tension values. At fixed DBP concentration, the addition of either or both of HNO 3 and NaNO 3 to an ionic strength of 1 M or 3 M led to significant reduction of the O-H VSFG peaks and enhancement of the C-H peaks. The origins of these experimental observations are attributed to both the increased HDBP molecules partitioning and adsorption to the interface and the protonation of the interfacial DBP- molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Surface Charge at an Oil/Water Interface Promoting the Enhancement and Assignment of Water Combination Bands

Given the different polarities and distinct chemical structures, it is no wonder that oil and water are incompatible neighbors. They can be forced to cooperate if provided sufficient energy or a suitable intermediary. Emulsions are the result of such a cooperation, where an emulsifying agent (the intermediary) in the form of a surfactant or a cocktail of surface stabilizing compounds allows for the ‘mixing’ of oil and water. Surfactant stabilized emulsions that are nanoscale in size (nanoemulsions) have caught the attention of many because of their newly discovered use in the synthesis of nanomaterials, oil recovery, drug delivery, and food science. Unlike smaller emulsions (e.g. microemulsions with 2-20 nm diameters), very little is known about the molecular level properties of nanoemulsions (NEs) whose diameters are in the 100-500 nm range. One thing that is well known is the macroscopic properties of droplets from these two size regimes are dramatically different, often for reasons not yet fully understood. Nanoemulsions do not form spontaneously and require an external shear beyond ordinary procedures used to create micelles. This shear is required to overcome the effects of surface tension and create droplets in the nanoscale regime. Once formed, however, they can be stable for long periods of time, making them very attractive for commercial and pharmaceutical products requiring a long shelf life. What is not well known is the interfacial molecular processes that lead to the formation and stability of nanoemulsions including those stabilized by surfactants and/or polymers. Such knowledge is important in understanding the molecular structure-function relationship that assists in predicting what additives and solvents are necessary for a specific application. Key to advancing nanoemulsion utility is to develop the needed molecular-level understanding of the structure and bonding interactions present at the surface of these soft nanoparticles.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Dynamics of the water-plasma interface in various discharge modes of atmospheric-pressure plasmas

This paper presents an experimental study of the dynamics of the water-plasma interface during the interaction high-voltage (HV) atmospheric pressure discharge with the water surface. Under the influence of the voltage applied to the pin type electrode, discharges are formed at the air layer between the sharp-tip HV electrode and deionized water. Three discharge regimes were identified from synchronized electrical waveforms and imaging: (i) a weak linear regime at 3–10.6 kV, (ii) a branching streamer regime at the maximum applied voltage of 12.6 kV, and (iii) a continuous arc-like regime occurring during the transition to lower voltage and higher current (e.g., U ≈ 3.2 kV, I ≈ 4.07 mA). High-speed shadowgraph images showed a symmetric interfacial cavity whose depth increased nonlinearly with voltage from h 0 ≈ 0.1 mm at 3 kV to h 0 ≈ 2.7 mm at 10.6 kV, reaching a maximum depth of h 0 ≈ 5.9 mm at 12.6 kV, while the cavity disappeared in the continuous-channel regime and was replaced by outward-propagating wave-like motion. The deformation of the water surface during the discharge-water interaction is governed by the balance between electric field forces, surface tension, gravitational forces, and electrohydrodynamic forces, whose relative contributions vary with the applied voltage. This effect was explained by quantitatively estimating the magnitudes of the acting forces and establishing the force balance. Optical emission spectroscopy in the continuous-channel regime indicated air plasma signatures and yielded a gas temperature of approximately 400 K, while prolonged operation (≈ 20 min) increased the water temperature to ~ 70 °C, reduced the water-layer thickness from 6 mm to 3 mm, and decreased pH from 7 to 4. These regime-resolved, quantitative results clarify how the dominant interfacial forcing shifts with discharge mode and provide a mechanistic basis for controlling plasma-water interactions in atmospheric-pressure applications.

atmospheric pressure plasma↗

DOE OSTI2024

Polymer-fiber-reinforced polymers with enhanced interfacial bonding between polypropylene fiber and polyethylene matrix

Self-reinforced composites (SRCs) consist of reinforcing fibers and a base matrix made of the same thermoplastic polymer, offering lightweight, recyclability, and sustainability benefits. However, limited research exists on composites where the reinforcing thermoplastic polymer fibers differ from the base thermoplastic matrix. Here, this study focuses on investigating the mechanical behavior of such composites and exploring different surface modification methods to enhance the fiber/matrix interfacial bonding using polypropylene fibers and a polyethylene matrix as an example. It is shown that surface treatment with a commercial adhesion promoter containing n-butyl acetate significantly improves the interfacial shear strength between polypropylene fibers and the polyethylene matrix, increasing it by 145% compared to other methods investigated. Additionally, increasing the length of the embedded polymer fiber in the matrix leads to a notable increase in specific interfacial energy. Consequently, the thermoplastic polymer-fiber-reinforced polymers (PFRPs) using surface-treated woven polypropylene fabrics and a polyethylene matrix exhibit a 20% higher tensile strength and a 65% higher toughness compared to non-treated PFRPs. This study also shows that specific mechanical properties (normalized by the composite density) of the investigated woven PFRPs are similar to those of non-treated SRCs under uni-axial tension. Particularly, their ductility outperforms carbon-/glass-/aramid-fiber-reinforced polymers by at least 6 times at a same fiber volume fraction. The investigation of such composites and the exploration of surface modification methods present important progress in the field of thermoplastic PFRPs, which serve as a solution for addressing concerns related to recyclability and sustainability.

Fiber pull-out↗

DOE OSTI2025

Probing specific ion effects at air-aqueous dibutyl phosphate interfaces using vibrational sum frequency generation spectroscopy

Molecular properties at air–liquid and liquid–liquid interface hold the key to many processes involving molecular transport across phase boundaries from aerosol formation to carbon cycling and material separation using solvent extraction techniques. Using dibutyl phosphate (DBP) as a representative for partially aqueous soluble surfactants, the specific ion effect (SIE) of the Hofmeister series cations Cs + , Na + , Li + , and Mg 2+ on the partition and interaction between surfactant molecules and water molecules in the air–aqueous interface are investigated using vibrational sum frequency generation spectroscopy and surface tension measurements. In the presence of 1 mM and 1M bulk aqueous phase ionic strength salt concentrations, fundamental qualitative relationships are observed for the salting out of DBP relative to bulk aqueous phase nitrate salt concentrations and the specific cations species. At 1 mM ionic strength, the interfacial charge and hence the interfacial potential modulates the electrostatic interactions; in particular, the counter cations partially screen the negatively charged interface induced by the DBP in a direct Hofmeister order. At 1M ionic strength, the electric field at the interface or interfacial potential is effectively neutralized, and the counter cations promote the partitioning of DBP to the interface depending on their specific interaction with the DBP head group and metal ion hydration properties. The present results lay a foundation to study SIEs of heavier metals on hydrophobic-aqueous DBP interfaces.

Louie, Christina [Washington State University, Pul↗

DOE OSTI2023

Large-scale molecular dynamics simulations of bubble collapse in water: Effects of system size, water model, and nitrogen

Molecular dynamics simulations in the microcanonical ensemble are performed to study the collapse of a bubble in liquid water using the single-site mW and the four-site TIP4P/2005 water models. To study system size effects, simulations for pure water systems are performed using periodically replicated simulation boxes with linear dimensions, L, ranging from 32 to 512 nm with the largest systems containing 8.7 × 10 6 and 4.5 × 10 9 molecules for the TIP4P/2005 and mW water models, respectively. The computationally more efficient mW water model allows us to reach converging behavior when the bubble dynamics results are plotted in reduced units, and the limiting behavior can be obtained through linear extrapolation in L –1 . Qualitative differences are observed between simulations with the mW and TIP4P/2005 water models, but they can be explained by the models’ differences in predicted viscosity and surface tension. Although bubble collapse occurs on time scales of only hundreds of picoseconds, the system sizes used here are sufficiently large to obtain bubble dynamics consistent with the Rayleigh–Plesset equation when using the models’ thermophysical properties as input. For the conditions explored here, extreme heating of the interfacial water molecules near the time of collapse is observed for the larger mW water systems (but the model underpredicts the viscosity), whereas heating is less pronounced for the TIP4P/2005 water systems because its larger viscosity contribution slows the collapse dynamics. The presence of nitrogen within the bubble only starts to affect bubble dynamics near the very end of the initial collapse, leading to an incomplete collapse and strong rebound for the mW water model. Although nitrogen is non-condensable at 300 K, it becomes highly compressed and reaches a liquid-like density near the collapse point. We find that the dissolution of nitrogen is much slower than the movement of the collapsing water front, and the re-expansion of the dense nitrogen droplet gives rise to bubble rebound. The incompatibility of the collapse and dissolution time scales should be considered for continuum-scale modeling of bubble dynamics. Finally, we also confirm that the diffusion coefficient for dissolved nitrogen is insensitive to pressure as the liquid transitions from a compressed to a stretched state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Sensitivity of the strength and toughness of concrete to the properties of the interfacial transition zone

We report Civil infrastructure is made primarily of concrete structures or components and therefore understanding durability and fracture behavior of concrete is of utmost importance. Concrete contains an interfacial transition zone (ITZ), a porous region surrounding the aggregates, that is often considered to be the weakest region in the concrete. The ITZ is poorly characterized and property estimates for the ITZ differ considerably. In this simulation study, representative concrete mesostructures are produced by packing coarse aggregates with realistic geometries into a mortar matrix. A meshless numerical method, peridynamics, is utilized to simulate the mechanical response including fracture under uniaxial compression and tension. The sensitivity of the stiffness and fracture toughness of the samples to the ITZ properties is computed, showing strong relationships between the ITZ properties and the effective modulus and effective yield strength of the concrete. These results provides insight into the influence of the poorly characterized ITZ on the stiffness and strength of concrete. This work showcases the applicability of peridynamics to concrete systems, matching experimental strength and modulus values. Additionally, relationships between the ITZ’s mechanical properties and the overall concrete strength and stiffness are presented to enable future design decisions.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Coupling of Adhesion and Anti-Freezing Properties in Hydrogel Electrolytes for Low-Temperature Aqueous-Based Hybrid Capacitors

Solid-state zinc-ion capacitors are emerging as promising candidates for large-scale energy storage owing to improved safety, mechanical and thermal stability and easy-to-direct stacking. Hydrogel electrolytes are appealing solid-state electrolytes because of eco-friendliness, high conductivity and intrinsic flexibility. However, the electrolyte/electrode interfacial contact and anti-freezing properties of current hydrogel electrolytes are still challenging for practical applications of zinc-ion capacitors. Here, we report a class of hydrogel electrolytes that couple high interfacial adhesion and anti-freezing performance. The synergy of tough hydrogel matrix and chemical anchorage enables a well-adhered interface between hydrogel electrolyte and electrode. Meanwhile, the cooperative solvation of ZnCl 2 and LiCl hybrid salts renders the hydrogel electrolyte high ionic conductivity and mechanical elasticity simultaneously at low temperatures. More significantly, the Zn||carbon nanotubes hybrid capacitor based on this hydrogel electrolyte exhibits low-temperature capacitive performance, delivering high-energy density of 39 Wh kg –1 at –60 °C with capacity retention of 98.7% over 10,000 cycles. With the benefits of the well-adhered electrolyte/electrode interface and the anti-freezing hydrogel electrolyte, the Zn/Li hybrid capacitor is able to accommodate dynamic deformations and function well under 1000 tension cycles even at –60 °C. This work provides a powerful strategy for enabling stable operation of low-temperature zinc-ion capacitors.

25 ENERGY STORAGE↗

DOE OSTI2021

Ultrasonic additive manufacturing of zirconium: Pilot results

Ultrasonic additive manufacturing (UAM) was successfully applied to the Zirconium material system to create a three-dimensional component prototype. The UAM process resulted in grain size refinement and significant localized deformation. Electron backscatter diffraction analysis revealed activated slip systems only applicable at high deformation temperatures, while texture analysis showed a decrease in measurable texture of the UAM build in comparison to the initial Zr-foil. Average tensile strengths in the X (sonotrode travel), Y (vibration), and Z (build) directions averaged 435 MPa, 458 MPa, and 359 MPa, respectively, demonstrating the viability of UAM for Zr-based materials. Furthermore, the delamination of some Z specimens along foil boundaries during loading suggests an interplay between interfacial Ti impurities introduced during welding and the spatial dependence of weld quality for this material system.

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↗