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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2021

Origins of electronic bands in the antiferromagnetic topological insulator MnBi 2 Te 4

Despite the rapid progress in understanding the first intrinsic magnetic topological insulator MnBi 2 Te 4 , its electronic structure remains a topic under debates. In this work, we perform a thorough spectroscopic investigation into the electronic structure of MnBi 2 Te 4 via laser-based angle-resolved photoemission spectroscopy. Through quantitative analysis, we estimate an upper bound of 3 meV for the gap size of the topological surface state. Furthermore, our circular dichroism measurements reveal band chiralities for both the topological surface state and quasi-2D bands, which can be well reproduced in a band hybridization model. A numerical simulation of energy-momentum dispersions based on a four-band model with an additional step potential near the surface provides a promising explanation for the origin of the quasi-2D bands. Our study represents a solid step forward in reconciling the existing controversies in the electronic structure of MnBi 2 Te 4 , and provides an important framework to understand the electronic structures of other relevant topological materials MnBi 2n Te 3n+1 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Enhanced Boundary Layer Height Detection Using Ceilometer, Surface Meteorology, and Radiation Products With a Random Forest Ensemble Method

This study develops and evaluates a Random Forest (RF) model for estimating planetary boundary layer height (PBLH) using 9 years of data from the Atmospheric Radiation Measurement Southern Great Plains (ARM SGP) user facility, with potential application in the NOAA Surface Radiation (SURFRAD) Network. The model integrates ceilometer, surface meteorology, and radiation measurements, and is trained using thermodynamic PBLH estimates derived from radiosondes. This approach aims to bridge gaps between aerosol-based and thermodynamic-based PBLH estimates. The RF model outperformed traditional methods during daytime and better captured transition periods, demonstrating improved accuracy and robustness. At ARM SGP, it showed a substantial reduction in both bias and RMSE, with a bias near zero (−4.9 m) compared with traditional Haar Wavelet (HW) (70.9 m) and Vaisala BL-View software (124.1 m), and an RMSE of 303.2 m, lower than both BL-View (566.9 m) and HW (404.6 m). During daytime hours, RF consistently outperformed both alternatives, maintaining lower bias and RMSE across all periods. At a second evaluation site, RF achieved the lowest overall RMSE (323.7 m), similar to HW (326.4 m) and significantly better than BL-View (738.3 m). However, all models showed reduced accuracy under stable nighttime conditions, limiting the reliability of PBLH estimates. Key predictors for the model included the lifting condensation level height (LCLH), aerosol gradients, and month for seasonal variability. The study underscores the potential of integrating machine learning with multiple data sets such as surface energy and thermodynamic data to advance PBLH estimation.

boundary layer height↗
DOE OSTI2022

Magnetism and magnetotransport in the kagome antiferromagnet Mn 3 Ge

We perform classical Monte Carlo and stochastic Landau-Lifshitz-Gilbert simulations to study the temperature-dependent magnetism of the kagome antiferromagnet Weyl metal Mn3Ge, and we find that a long-range chiral order sets in at a transition temperature well below the Neel temperature (T N ). Based on the crystalline symmetries imposed by the chiral magnetic order, we argue for the presence of multiple isoenergetic Weyl nodes (nodes that are at the same energy and with a congruent Fermi surface around them) near the chemical potential. Using the semiclassical Boltzmann equations, we show that the combined contribution to the net longitudinal magnetoconductance (LMC) and the planar Hall conductance (PHC) from tilted Weyl nodes can lead to signatures that are qualitatively distinct from that of a single pair of Weyl nodes. In particular, we show that magnetic orders with different chiralities can give rise to different periods in LMC and PHC as a function of the in-plane magnetic field direction. As a result, this is ultimately related to differences in the symmetry-imposed constraints on the Weyl nodes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

WBG-Enabled Current-Source Inverters for Integrated PM Motor Drives

This project takes advantage of the special capabilities of wide-bandgap (WBG) power semiconductor devices to develop innovative power electronics in the form of new current-source inverters (CSIs) that offer promising advantages over the dominant voltage-source inverter (VSI) topology. These CSIs, in turn, are extremely well-suited for integration into the housings of permanent magnet synchronous machines (PMSMs) to form integrated motor drives (IMDs). These new IMDs offer great promise for achieving major energy savings in a wide variety of applications that benefit from adjustable-speed control, including air conditioners, laundry appliances, industrial pumps/compressors, electric vehicles, and aerospace drives. WBG devices play a critical enabling role in these new IMDs because of their transformative features including much higher switching frequencies, lower losses, and compatibility with high operating temperatures. When incorporated into new CSI-based IMDs with PM machines, these WBG switches open the door to achieving major increases in power density, drive system efficiency, and fault tolerance, as well as substantial reductions in electromagnetic interference (EMI), manufacturing cost, and temperature-induced failures. The higher operating temperature capability of WBG devices compared to conventional silicon power devices is very appealing in IMD applications because the power electronics is mounted in close proximity to the motors which typically operate at temperatures well above the maximum limits of today’s silicon-based power electronics. This project has succeeded in designing, building, and testing multiple prototype versions of this CSI-IMD that have overcome many technical challenges in order to demonstrate the impressive performance improvements that can be achieved by the WBG-enabled CSI-IMD. The preliminary demo and bench-top versions of the current-source inverter developed during the first two years of the project were critical to laying the technical foundations for the 3 kW prototype CSI-IMD unit that was successfully built and tested during the third year. This prototype CSI-IMD unit was designed to fit within the housing envelope of the original permanent magnet (PM) synchronous machine in order to meet the demanding power density requirements that were set at the beginning of the project. All of the remaining performance objectives set for the prototype CSI-IMD unit including efficiency and electromagnetic interference (EMI) were also met. The last 18 months of the project were devoted to developing further enhancements of the WBG-enabled CSI-IMD technology that better prepare it for commercial production. More specifically, an upgraded version of the prototype CSI-IMD unit was developed that moved the power electronics into the same machine housing chamber as the motor, substantially raising the thermal demands on the power electronics. Tests with five different combinations of motor enclosure types and air cooling configurations were evaluated. Importantly, this work confirmed that the power electronics can deliver its full rated power and still operate well within its maximum temperature limits even for worst-case conditions when the housing is “totally-enclosed” without any openings for air to enter or exit the enclosure, and no blower/fan is provided to blow air over the outside surface of the enclosure. Reaching the full performance and energy-savings potential of this disruptive CSI-IMD motor drive technology is highly consistent with ARPA-E’s stated mission to “enhance the economic and energy security of the United States” while also supporting its commitment to “ensure that the U.S. maintains a technological lead in developing and deploying advanced energy technologies”. Follow-on projects are under way to explore the scalability of WBG-enabled CSI-IMD technology to 100 kW (peak) electric vehicle traction drives and fault-tolerant modular motor drives for future electrified aircraft propulsion applications.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Additive Manufacturing of Thermal Energy Storage Composites with Microencapsulated Phase Change Materials Supported in a Multi-Polymer Matrix

Advanced manufacturing techniques, such as additive manufacturing (AM), that can directly integrate phase change materials (PCMs) have garnered interest in recent years due to their potential for development of highly efficient thermal energy storage architectures. Complex, high surface area geometries embedded with PCMs that are only feasible with AM can improve thermal management with reduced material waste. Our work focuses on developing composite filaments with microencapsulated phase change materials (MEPCM) bound within a single or dual polymer matrix that can be processed through standard filament extruders and additively manufactured using off-the-shelf 3D printers. Polymer powders, rather than polymer pellets, were key to homogenously mixed filaments achieving high MEPCM loadings with no deterioration in thermal energy storage (TES) capability during extrusion. Composite filaments contain upwards of 60 wt% MEPCM and were printed without loss in feature resolution, print speed, or layer adhesion. Storage enthalpies of printed composites range from 100 - 130 kJ/kg, which were within 5% of the theoretical enthalpy based on weight fraction of MEPCM and maintained enthalpies within 1% over 500 thermal cycles. We can reliably manufacture low density, high surface area structures like 15% gyroid infill, along with dense, compact pucks at a 100% concentric infill. Prints were also scalable to a 900 cm3 honeycomb infill heat exchanger model that has an estimated energy storage capacity of 9 Wh.

3D printing↗
DOE OSTI2020

Simulating the Effects of Surface Energy Partitioning on Convective Organization: Case Study and Observations in the US Southern Great Plains

Realistic cloud-resolving simulations were performed to study the effects of surface energy partitioning (surface sensible and latent heat fluxes) on the organization of isolated convection into larger mesoscale convective systems (MCSs) near the US Southern Great Plains. The role of cold pools in mediating surface-convection interactions was explored. Better organized MCSs tended to occur in the experiments with perturbed wetter soil (and more active vegetation), regardless of the effects of soil moisture on the diurnal timing of convective triggering. Wetter soil led to shallower boundary layers and more convective available potential energy than drier soil. The roles of cold pools on convection are lifecycle-stage dependent: A dry surface allows more numerous colliding cold pools, thereby aiding in convective triggering by reducing entrainment in the early stages, and providing gust-front uplift in later stages. However, horizontal propagation of the cold-pool density current can outrun the convective system, creating a slantwise updraft and thus weakening the gust front uplift in later stages. This effect calls into question previous cold pool parameterizations, in which the gust front uplift is mainly proportional to the negative buoyancy of cold air. Finally, both the model and observation show an enhancement of surface latent heat flux during the passage of a gust front at night, suggesting that a positive feedback between the surface and convection helps MCSs to persist into the nighttime.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Processes at lithium-hydride/deuteride surfaces upon low energy impact of H/D

Sputtering, reflection, and retention processes at amorphous and crystalline lithium hydride surfaces due to impact of low energy (1–100 eV) hydrogen and deuterium atoms over the range of 0 o −85 o angle of incidence at 300 K surface temperature were investigated by atomistic computational methods. Classical molecular dynamics simulations were performed with improved reactive bond-order force field (ReaxFF) potentials that include long-range polarization effects. In addition to probabilities of surface processes, the energy and angular spectra of ejected particles were obtained. Comparison of these results with those previously obtained on pristine lithium surfaces indicates the importance of saturation of the Li surface and near-surface region with hydrogen. We show that such saturation, which is typical in both laboratory and fusion device experiments with lithium coating of the plasma-facing surfaces, significantly changes the surface processes with hydrogen irradiation in the understudied low-energy region of impact energies.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Clean Energy Technology Applications on US Mine Land: Technical Analysis

As the United States transitions toward a clean energy economy, an opportunity exists for redeveloping the more than 17,000 mine land sites located across the nation with clean energy technologies, which have a combined potential for generating more than 85 GW of clean electricity. This report provides an overview of the potential of demonstrating and deploying clean energy projects on current and former mine land. Clean energy project refers to a project that demonstrates one or more of the following technologies: solar; microgrids; geothermal; direct air capture; fossil-fueled electricity generation with carbon capture, utilization, and sequestration; energy storage, including pumped storage hydropower and compressed air energy storage; and advanced nuclear technologies. The report discusses the following technologies and their potential for creating jobs and generating tax revenue that would result in direct and indirect benefits to the local economy: Solar photovoltaics (PV) is being developed on current and former mine land in various parts of the world, including the United States. This approach is attractive because it requires limited infrastructure investment and would utilize the bare surfaces of mines and tailing ponds. Solar resource availability may be greater in the southern regions, including the Interior and Appalachian Basins and the southwestern United States. However, since some mine land sites include areas of significant change in elevation, the deployment of PV on mine land may require sophisticated planning to account for shading and irradiance, or may require regrading of the areas. PV does not create significant environmental risks and generally does not face public resistance; Geothermal systems are often spatially and genetically associated with ore deposits, and in some cases, they have been discovered while in search for epithermal mineral resources. Numerous diverse geothermal applications have been employed at mine land around the world, including power generation, mineral extraction from geothermal brines, process heating, direct use for other mining operations, and direct use for non-mining operations and subsurface energy storage, including geothermal heat pumps. Case studies highlighting these applications provide key lessons relating to identifying drivers and barriers to geothermal resource deployment and can be used to create screening tools for identifying the types and locations of mine land most amenable to utilizing geothermal resources; Carbon capture, utilization, and sequestration technologies include direct air capture (DAC) and enhanced weathering. DAC technologies include air contactors, regeneration systems, and CO 2 compression systems. Captured CO 2 can be converted to valuable feedstocks or possibly injected into abandoned subsurface mines where it would be absorbed by alkaline rock waste and mine tailings or by the porous minerals along the walls of the mine. DAC systems can be coupled with energy sources such as wind, solar, grid, or geothermal. Many DAC systems require a source of water or steam; however, some are expected to be net producers of water. Local impacts of DAC systems are expected to be low, and are related to land footprint, material disposal, and upstream impacts of energy and material production; Compressed air energy storage is an established energy storage technology in salt caverns. It has the potential for implementation in underground mines by pressurizing and storing a large amount of air using electrical compressors when excess electricity is available. When a need for discharge emerges, the air is used to spin turbines and produce the necessary volume of electricity. Abandoned or unused mine openings, including shafts, adits, access tunnels, and mined workings of any orientation, offer potential for vast amounts of compressed air energy storage if the site characteristics meet operational requirements; Pumped hydropower storage can be implemented in surface and subsurface mines. In surface mine applications, both reservoirs may be located in a mine pit or artificial reservoirs made of excavated materials. In subsurface mines, the lower reservoir may be implemented by waterproofing and flooding mine shafts and tunnels. The water is then pumped from the lower reservoir to the upper reservoir during periods of low load and high production, and it is discharged through the turbines during periods of peak demand. The potential environmental damages associated with acidity of mine water or the presence of toxic chemicals incentivizes the development of closed-loop technologies, in which water circulates inside the pumped hydropower facility without being discharged into the external water basins; Advanced nuclear energy technologies include small modular reactors, which can be deployed locally to produce electricity and heat. Such units require seismic stability and a supply of cooling water, but population constraints may exist in some areas. Therefore, remote mine land could represent an optimal location for siting advanced nuclear energy technologies.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

Surface Functionalized Barium Titanate Nanoparticles: A Combined Experimental and Computational Study

Barium titanate (BTO) nanoparticles show great potential for use in electrostatic capacitors with high energy density. This includes both polymer composite and sintered capacitors. However, questions about the nanoparticles' size distribution, amount of agglomeration, and surface ligand effect on performance properties remain. Reducing particle agglomeration is a crucial step to understanding the properties of nanoscale particles, as agglomeration has significant effects on the composite dielectric constant. BTO surface functionalization using phosphonic acids is known reduce BTO nanoparticle agglomeration. We explore solution synthesized 10 nm BTO particles with tert-butylphosphonic acid ligands. Recent methods to quantifying agglomeration using an epoxy matrix before imaging shows that tert-butylphosphonic acid ligands reduce BTO agglomeration by 33%. Thermometric, spectroscopic, and computational methods provide confirmation of ligand binding and provide evidence of multiple ligand binding modes on the BTO particle surface.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Vibrational sum frequency generation spectroscopy reveals the inertness of chromium oxide (001) surfaces

Nanoengineered metal oxides such as Cr(III)-oxide (chromia) films have diverse potential applications in corrosion inhibition, remediation, energy generation, catalysis, data storage, and biological and environmental systems. Concerns about material degradation or oxidation to toxic chromate necessitate an understanding of chromia/aqueous interfaces, beginning with their hydroxylation and hydration behavior. Vibrational sum-frequency generation spectroscopy (vSFG) provides specific molecular-level information about water at the oxide/aqueous junction with high surface selectivity. To overcome the strong absorber problem typical of certain metal oxides in the UV–visible range during nonlinear optical studies, we employed molecular beam epitaxy to deposit transparent, epitaxial nanofilms of Cr 2 O 3 with (001) crystalline orientation on sapphire (Al 2 O 3 (001)) substrates, as confirmed by atomic force microscopy and X-ray diffraction. vSFG spectra of the air and water interfaces of the Cr 2 O 3 (001) films reveal hydroxyl features corresponding to both dissociated and molecular water on the surface. In contrast to the dangling hydroxyls found on the bare Al 2 O 3 (001) substrate, the hydroxyl groups on the deposited Cr 2 O 3 (001) nanofilm do not readily undergo isotopic H/D exchange when exposed to varying forms of D 2 O under ambient conditions. When considering the chemistries of the corresponding trivalent cations in aqueous solution, the finding is at variance with their similar acidities and proton exchange dynamics but consistent with markedly slower inner-sphere water ligand exchange of hexaquo Cr 3+ . This finding challenges the idea that proton and water exchange at the water/Cr₂O₃(001) interface are solely correlated, driven by the strength of metal-bridging oxygen bonds and surface hydroxyl distribution, without direct causation.

Boamah-Agyemang, Mavis D. [Pacific Northwest Natio↗
DOE OSTI2022

High-throughput screening to predict highly active dual-atom catalysts for electrocatalytic reduction of nitrate to ammonia

Ammonia is an essential chemical owing to its importance in fertilizer production and other industrial applications. Electrocatalytic nitrate reduction to ammonia (NO 3 RR) holds great promise for low-temperature ammonia production while simultaneously addressing nitrate-based environmental concerns. To provide the mechanistic understanding needed to design an effective electrocatalyst, we systematically investigated the catalytic performance of metal-based dual-atom catalysts (DACs) anchored on two-dimensional (2D) expanded phthalocyanine (Pc) for NO 3 RR. We found that NO 3 RR can efficiently produce ammonia on Cr 2 -Pc, V 2 -Pc, Ti 2 -Pc, and Mn 2 -Pc surfaces with low limiting potentials of – 0.02, – 0.25, – 0.34, and – 0.41 V RHE , respectively. Moreover, using the free energy difference of *NO 3 - and *H as a descriptor, we found that the hydrogen evolution reaction is significantly suppressed on the DAC surface due to an ensemble effect in which the two metal atoms cooperate to selectively form ammonia. We performed high-throughput screening to develop an efficient metal-based DAC for NO 3 - reduction, followed by a mechanistic study to elucidate the NO3RR pathway on the DAC. Finally, this work provides design information for advancing sustainable ammonia synthesis under ambient conditions.

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

Moisture Performance of Unvented Attics With Vapor Diffusion Ports and Buried Ducts in Hot, Humid Climates

Energy efficiency measures, such as cool roofs, radiant barriers, interior radiative control coatings, and buried ducts are increasing in popularity and are promoted by energy codes because of their energy-saving potential. However, these strategies can also pose moisture risks in attics by lowering surface temperatures and increasing condensation potential and moisture accumulation. Of particular concern in hot-humid climates is dripping condensation on cold air-conditioning ductwork in the summer - commonly referred to as duct "sweating" - which threatens the attic floor with conditions conducive to mold growth and rot. One strategy to mitigate these moisture issues is to wrap ductwork in thicker duct-wrap insulation with an integrated exterior vapor barrier, but thick duct wrap can be difficult to come by, expensive, and unwieldy to work with. This study explores an alternative strategy of reducing moisture issues while embracing energy efficiency by using unvented attics with vapor diffusion ports and buried ductwork in hot-humid climates. Vapor diffusion ports have been studied so far in a wide range of U.S. climates, mostly in the context of conditioned attics. In this study, the strategy is implemented in the novel context of hot-humid climates with ductwork sitting atop blown-in attic floor insulation in unconditioned attics. Using a combination of field experiments and hygrothermal modeling, the findings of this project indicate that an unvented attic with vapor diffusion ports and buried ducts may be a key part of a successful low-cost method for reducing the attic moisture load by venting excess moisture out of the attic, keeping duct-jacket surfaces above dew point temperature, and keeping the roof deck safe from winter moisture accumulation.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Conducting thermal energy to the membrane/water interface for the enhanced desalination of hypersaline brines using membrane distillation

Membrane distillation (MD) is a membrane-based thermal desalination process capable of treating hypersaline brines. Standard MD systems rely on preheating the feed to drive the desalination process. However, relying on the feed to carry thermal energy is limited by a decline of the thermal driving force as the water moves across the membrane, and temperature polarization. In contrast, supplying heat directly into the feed channel, either through the membrane or other channel surfaces, has the potential of minimizing temperature polarization, increasing single-pass water recoveries, and decreasing the number of heat exchangers in the system. When solar thermal energy can be utilized, particularly if the solar heat is optimally delivered to enhance water evaporation and process performance, MD processes can potentially be improved in terms of energy efficiency, environmental sustainability, or operating costs. Here we describe an MD process using layered composite membranes that include a high-thermal-conductivity layer for supplying heat directly to the membrane-water interface and the flow channel. Here, the MD system showed stable performance with water flux up to 9 L/m 2 /hr, and salt rejection >99.9% over hours of desalinating hypersaline feed (100 g/L NaCl). In addition to bench-scale system, we developed a computational fluid dynamics model that successfully described the transport phenomena in the system.

42 ENGINEERING↗
DOE OSTI2024

Improving the reliability of machine learned potentials for modeling inhomogeneous liquids

The atomic-scale response of inhomogeneous fluids at interfaces and surrounding solute particles plays a critical role in governing chemical, electrochemical, and biological processes. Classical molecular dynamics simulations have been applied extensively to simulate the response of fluids to inhomogeneities directly, but are limited by the accuracy of the underlying interatomic potentials. Here, we use neural network potentials (NNPs) trained to ab initio simulations to accurately predict the inhomogeneous responses of two distinct fluids: liquid water and molten NaCl. Although NNPs can be readily trained to model complex bulk systems across a range of state points, we show that to appropriately model a fluid's response at an interface, relevant inhomogeneous configurations must be included in the training data. In order to sufficiently sample appropriate configurations of such inhomogeneous fluids, we develop protocols based on molecular dynamics simulations in the presence of external potentials. We demonstrate that NNPs trained on inhomogeneous fluid configurations can more accurately predict several key properties of fluids—including the density response, surface tension and size-dependent cavitation free energies—for liquid water and molten NaCl, compared to both empirical interatomic potentials and NNPs that are not trained on such inhomogeneous configurations. This work therefore provides a first demonstration and framework to extract the response of inhomogeneous fluids from first principles for classical density-functional treatment of fluids free from empirical potentials.

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

A Deep Learning-Driven Sampling Technique to Explore the Phase Space of an RNA Stem-Loop

The folding and unfolding of RNA stem-loops are critical biological processes; however, their computational studies are often hampered by the ruggedness of their folding landscape, necessitating long simulation times at the atomistic scale. Here, we adapted DeepDriveMD (DDMD), an advanced deep learning-driven sampling technique originally developed for protein folding, to address the challenges of RNA stem-loop folding. Although tempering- and order parameter-based techniques are commonly used for similar rare-event problems, the computational costs or the need for a priori knowledge about the system often present a challenge in their effective use. DDMD overcomes these challenges by adaptively learning from an ensemble of running MD simulations using generic contact maps as the raw input. DeepDriveMD enables on-the-fly learning of a low-dimensional latent representation and guides the simulation toward the undersampled regions while optimizing the resources to explore the relevant parts of the phase space. We showed that DDMD estimates the free energy landscape of the RNA stem-loop reasonably well at room temperature. Our simulation framework runs at a constant temperature without external biasing potential, hence preserving the information on transition rates, with a computational cost much lower than that of the simulations performed with external biasing potentials. Here, we also introduced a reweighting strategy for obtaining unbiased free energy surfaces and presented a qualitative analysis of the latent space. This analysis showed that the latent space captures the relevant slow degrees of freedom for the RNA folding problem of interest. Finally, throughout the manuscript, we outlined how different parameters are selected and optimized to adapt DDMD for this system. We believe this compendium of decision-making processes will help new users adapt this technique for the rare-event sampling problems of their interest.

Gupta, Ayush↗
DOE OSTI2024

Estimation of Surface Sensible Heat Flux due to Precipitation over CONUS and Its Impact on Urban Extreme Precipitation Modeling

The surface sensible heat flux induced by precipitation (Q P ) is a consequence of the temperature difference between the surface and the rain droplets. Despite its seemingly negligible nature, Q P is frequently omitted from both meteorological and climatological models. Nevertheless, it is important to acknowledge the numerous occasions in which the instantaneous values of Q P can be significant, particularly during extreme precipitation events. This study undertakes a comprehensive assessment of Q P across the contiguous United States (CONUS) utilizing high-resolution reanalysis, observational data, and numerical modeling to examine the influence of Q P on precipitation and the surface energy budget. The findings indicate that the spatial distribution of Q P climatology is analogous to that of precipitation, with magnitudes ranging from 2 to 3 W m -2 predominantly over the Midwest and Southeast regions. A seasonal analysis of QP reveals that the highest values occurring during the June–August (JJA) period, averaging 3.18 W m -2 . Peak Q P values of approximately 4 W m -2 are observed during JJA over the Great Plains region. We hypothesize that the Q P during an extreme precipitation event would be nonnegligible and have a significant impact on the local weather. To test this conjecture, we perform high-resolution simulations with and without Q P during an extreme precipitation event over the Chicago Metropolitan Area (CMA). The results show that the Q P may be a dominant factor compared to other components of surface heat flux during the zenith of precipitation hours. Also, Q P has the potential to not only diminish precipitation but also alter and reconfigure the remaining surface energy budget components.

54 ENVIRONMENTAL SCIENCES↗