Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solutions center”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Module-Level Solutions For Degradation by Ionization Damage

After years of improving module efficiency while targeting mean degradation rates of 0.5 to 0.6 %/y for crystalline silicon technology, there is much evidence that the degradation rates are now increasing significantly. Contributing factors include:Radiation damage (UV-Light Induced Degradation): Trina solar: -4.5%/y degradation rate in Singapore, DOE National Laboratory Regional Test Centers showed degradation of -1% < r < -2%/y in crystalline silicon modules, Jinko Solar: -4% < r < -7% efficiency loss from 540 MJ/m-2 of UV-A light, ISFH: 15% relative power loss during of 1.8 GJ/m-2 UV exposure, attributed to H+. Electrical bias from positive system voltage (e.g., +1000 V, +1500 V) can drive ions and metallization through the encapsulation, NREL: Ion transport can affect the cell passivation, resulting in power loss of 5% to 40% in p-PERC+ (bifacial), NREL: Damage at cell rear with up to 17% power at cell fronts in n-PERT modules. Delamination can also occur due to cell surface reactions driven by light and bias.

module efficiency↗
DOE OSTI2026

Studies of Localization Effects on Transverse Kinematic Imbalance in the GENIE Generator

The forthcoming Deep Underground Neutrino Experiment (DUNE) requires precise modeling of neutrino--nucleus interactions to achieve its neutrino-oscillation measurement goals. GENIE, the Monte Carlo event generator used by DUNE and other Fermilab-based experiments as a central value model for neutrino-nucleus cross sections, exhibits known discrepancies in comparison to MicroBooNE cross section data, especially when considering neutrino scattering events' transverse kinematic imbalance (TKI). For quasielastic-like measurements, the data show a larger transverse missing momentum tail while maintaining a peak similar to that predicted by the baseline GENIE model. Previous attempts at a solution include variations of final-state interaction (FSI) strengths, incorrectly decreasing both the peak and tail, thus leaving the discrepancy unresolved. This work investigates whether inconsistencies between local and global treatments of intranuclear physics contribute to the observed mismodeling. In this work, GENIE FSI routines have been modified to directly utilize intranuclear particle positions, thereby introducing a newly localized momentum treatment in correlation with the nominal density in both the hA and hN intranuclear cascade models for their 2018 and 2025 variants. Meson-exchange-current (MEC) localization was also tested for the hN 2018 model. Comparisons of the intranuclear nucleons' ("scattering center") momentum and its radial dependence validate the localization implementation. However, while such consistent localization produces modest changes in the proton momentum spectrum at low-to-medium momenta, only small changes are observed in the TKI distributions. These changes are insufficient to account for the discrepancy with MicroBooNE data. Although a localized treatment improves the internal consistency of the GENIE model, the origin of the TKI discrepancy remains unresolved by such a solution.

Bulla, Braden [Ctr. Coll., Danville; Fermilab]↗
DOE OSTI2027

Conceptual design of highly-constrained splitters for the FFA@CEBAF energy upgrade study

The Continuous Electron Beam Accelerator Facility (CEBAF) at Jefferson Lab is investigating a significant energy upgrade utilizing Fixed-Field Alternating-gradient (FFA) recirculating arcs. This upgrade requires the design of complex horizontal beam splitters to manage up to six concurrent beam passes. This paper presents the conceptual design of these splitters, which are subject to severe physical constraints imposed by the existing accelerator tunnel and multifaceted beam dynamics requirements for matching into the permanent-magnet FFA arcs. The design methodology, centered on multi-pass simulations in the Bmad toolkit, is detailed from the initial geometric layout through the advanced optics matching. Key results include a robust geometric arrangement that fits within the spatial boundaries and the development of multiple, flexible optics matching solutions. Furthermore, the design integrates a viable scheme for extracting high-energy beams for the experimental halls, a critical operational requirement. This work establishes a comprehensive and viable conceptual design, forming a baseline for future engineering and performance optimization studies.

Bodenstein, R.M. [Thomas Jefferson National Accele↗
DOE OSTI2020

Application and Certification of Comparative Vacuum Monitoring Sensors for Structural Health Monitoring of 737 Wing Box Fittings

Multi-site fatigue damage, hidden cracks in hard-to-reach locations, disbonded joints, erosion, impact, and corrosion are among the major flaws encountered in today's extensive fleet of aging aircraft and space vehicles. The use of in-situ sensors for real-time health monitoring of aircraft structures are a viable option to overcome inspection impediments stemming from accessibility limitations, complex geometries, and the location and depth of hidden damage. Reliable, structural health monitoring systems can automatically process data, assess structural condition, and signal the need for human intervention. Prevention of unexpected flaw growth and structural failure can be improved if on-board health monitoring systems are used to continuously assess structural integrity. Such systems are able to detect incipient damage before catastrophic failures occurs. Condition-based maintenance practices could be substituted for the current time-based maintenance approach. Other advantages of on-board distributed sensor systems are that they can eliminate costly, and potentially damaging, disassembly, improve sensitivity by producing optimum placement of sensors and decrease maintenance costs by eliminating more time- consuming manual inspections. This report presents a Sandia Labs-aviation industry effort to move SHM into routine use for aircraft maintenance. This program addressed formal SHM technology validation and certification issues so that the full spectrum of concerns, including design, deployment, performance and certification were appropriately considered. The Airworthiness Assurance NDI Validation Center (AANC) at Sandia Labs, in conjunction with Boeing, Delta Air Lines, Structural Monitoring Systems Ltd., Anodyne Electronics Manufacturing Corp. and the Federal Aviation Administration (FAA) carried out a certification program to formally introduce Comparative Vacuum Monitoring (CVM) as a structural health monitoring solution to a specific aircraft wing box application. Validation tasks were designed to address the SHM equipment, the health monitoring task, the resolution required, the sensor interrogation procedures, the conditions under which the monitoring will occur, the potential inspector population, adoption of CVM into an airline maintenance program and the document revisions necessary to allow for routine use of CVM as an alternate means of performing periodic structural inspects. To carry out the validation process, knowledge of aircraft maintenance practices was coupled with an unbiased, independent evaluation. Sandia Labs designed, implemented, and analyzed the results from a focused and statistically-relevant experimental effort to quantify the reliability of the CVM system applied to the Boeing 737 Wing Box fitting application. All factors that affect SHM sensitivity were included in this program: flaw size, shape, orientation and location relative to the sensors, as well as operational and environmental variables. Statistical methods were applied to performance data to derive Probability of Detection (POD) values for CVM sensors in a manner that agrees with current nondestructive inspection (NDI) validation requirements and also is acceptable to both the aviation industry and regulatory bodies. This report presents the use of several different statistical methods, some of them adapted from NDI performance assessments and some proposed to address the unique nature of damage detection via SHM systems, and discusses how they can converge to produce a confident quantification of SHM performance An important element in developing SHM validation processes is a clear understanding of the regulatory measures needed to adopt SHM solutions along with the knowledge of the structural and maintenance characteristics that may impact the operational performance of an SHM system. This report describes the major elements of an SHM validation approach and differentiates the SHM elements from those found in NDI validation. The activities conducted in this program demonstrated the feasibility of routine SHM usage in general and CVM in particular for the application selected. They also helped establish an optimum OEM-airline-regulator process and determined how to safely adopt SHM solutions. This formal SHM validation will allow aircraft manufacturers and airlines to confidently make informed decisions about the proper utilization of CVM technology. It will also streamline the regulatory actions and formal certification measures needed to assure the safe application of SHM solutions.

42 ENGINEERING↗
DOE OSTI2022

Science and Solutions for a Planet Under Pressure: Synthetic Research, Analytics, and Community Engagement on Five Priority Theme Areas

Support from the Department of Energy proved instrumental in supporting the work of Global Council for Science and the Environment (GCSE) on synthetic research, analysis, and community engagement, oriented around the themes and discussions from the GCSE Drawdown 2021 Conference — Research to Action: Science and Solutions for a Planet Under Pressure. These efforts helped inform and enhance the societal relevance of the knowledge that has been generated through the conference and, at the same time, advance the scientific insights acquired through synthesis and scientific translation efforts. The work centered on three main objectives: 1) documentation of findings from the completed GCSE 2021, culminating in the publication of a synthesis report; 2) development and vetting of integrated use cases that illustrate synthesis methods and enable scientific outreach, communications, and translational discussions among universities and the broader research community; and 3) strengthening the scientific basis for energy education research by sharing findings in university-educator fora.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Computationally-Guided ODS Refractory HEAs via Additive Manufacturing

This document is the final report for ARPA-E ULTIMATE project DE-AR0001422. To achieve higher efficiency turbine operation, the RTX Technology Research Center (RTRC) explored the use of additive manufacturing (AM) to produce test specimens of novel high entropy alloys (HEA) enhanced with oxide dispersion strengthening (ODS). HEAs, with multiple principal elements, offer a vast alloy design space rich with unique properties especially enhanced solid solution strengthening. ODS alloys historically are some of the best performing materials for high temperature applications. Leveraging AM enables near net shape part production and incorporation of ODS into the alloy. Density, melting temperature, and elevated temperature strength properties of several novel HEA compositions were found to meet or exceed program targets; however, obtaining sufficient room temperature ductility in these alloys proved challenging. In order to obtain an alloy with sufficient ductility, the team explored more conventional refractory alloy compositions, focusing on a niobium based alloy. This inherently ductile alloy also satisfied program metrics for cost, density, and fracture toughness. Results showed the oxide dispersions increased high temperature strength, and improved creep resistance compared to the conventional alloy, though the program creep metric was not met. Further tuning of the alloy chemistry, and optimization of the oxide dispersion, is likely required to improve creep resistance.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Understanding the microstructural stability in a y’-strengthened Ni-Fe-Cr-Al-Ti alloy

Ni-Fe-Cr-Al-Ti alloys, with Ni levels near 50 wt.%, have the potential to develop a microstructure consisting of a face-centered cubic ? matrix with the homogeneous precipitation of fine ordered ?’ precipitates similar to traditional Ni-based superalloys with significantly greater Ni content. Scanning electron microscopy (SEM), transmission electron microscopy (TEM), atom probe tomography (APT), and CALPHAD -based thermodynamic modeling were employed to understand the phase stabilities and microstructural evolution in an age-hardenable Ni-27Fe-18Cr-1Co-1.6Al-3.75Ti-1.2Mo-0.03C (wt%) alloy. The primary heat-treatment of solution annealing at 1121°C for 4h and age-hardening treatment at 760 °C for 16h resulted in a microstructure consisting of fine ?’ precipitates in an austenitic matrix along with grain boundary carbides, consistent with thermodynamic calculations. Long-term aging at 900 °C for 250h resulted in the coarsening of ?’ along with a change in the morphology of the precipitates from spherical to a more cuboidal shape. In addition, ? phase formation was observed concomitant with the partial dissolution of the ?’ phase. The ability of computational thermodynamic models to predict microstructural characteristics is discussed.

Gwalani, Bharat↗
DOE OSTI2025

Surface engineering of Pt nanocatalysts with transition metal oleates for selective catalysis: a case study on the hydrogenation of α,β-unsaturated aldehydes

Selective and active catalysts enable effective use of feedstocks, reduced energy consumption and waste generation. Tuning the electronic structure of heterogeneous metal nanocatalysts via their surface modifications is a promising strategy to design highly selective and active catalysts for the synthesis of harder to make and more cost-efficient products. Here, we introduce transition metal oleates as a new class of ligands to engineer the catalytically active and very selective surface in organic solvents. Using citral hydrogenation and 5 nm Pt NPs as a model reaction and model catalytic system, respectively, we show that surface engineering of Pt nanocatalysts with metal oleates allows synthesis of desired partially hydrogenated product (geraniol) with ∼90% conversion with selectivity over 93%. We demonstrate that the selective synthesis of the unsaturated alcohols catalyzed by Pt NPs modified by adsorption of the transition metal salts cannot be explained by the widely accepted mechanism of preferred coordination of C$=$O groups by Lewis acids (e.g. partially oxidized transition surface metals). Our results indicate that C$=$O groups prefer to bind to negatively charged surfaces. We propose the explanation on how the adsorption of transition metal oleates can result in the increased electron density at the surface of Pt nanoparticles. Our study not only provides reliable solutions to selective hydrogenation but opens a new possibility of using metal oleates for the electronic ligand effect.

Kwon, Soon Gu [Argonne National Laboratory (ANL), ↗
DOE OSTI2024

Regulating the Third Metal to Design and Engineer Multilayered NiFeM (M: Co, Mn, and Cu) Nanofoam Anode Catalysts for Anion-Exchange Membrane Water Electrolyzers

In this study, alkaline anion-exchange membrane water electrolyzers (AEMWEs) for green hydrogen production have received intensive attention due to their feasibility of using earth-abundant platinum group metal (PGM)-free catalysts. Herein, the third metal is incorporated into NiFe-based catalysts to regulate their electronic structures and morphologies, aiming to achieve sufficient oxygen evolution reaction (OER) activity and performance in AEMWEs. The ternary NiFeM (M: Cu, Co, or Mn) catalysts are featured with multiple layered structures and nanofoam network morphologies, consisting of highly OER-active amorphous Ni-rich oxide shells and electrically conductive metallic alloy cores. The physical and electronic perturbations to the NiFe induced by a third element lead to a fine-tuning of the redox ability of the metal sites at the reaction centers, which breaks the scaling relationship between OH* and O* intermediates at the reaction centers. Thus, the unique structural configuration and electronic regulation simultaneously benefit catalytic activity and performance improvements. These NiFeM nanofoam catalysts demonstrated promising anode performance in actual AEMWEs, comparable to the IrO 2 reference, especially at high current densities. Notably, using various electrolytes (e.g., KOH solution or pure water) for AEMWEs exhibited a different performance trend among studied NiFeM catalysts, likely due to dynamic changes of catalysts under various OER environments.

25 ENERGY STORAGE↗
DOE OSTI2025

Conceptualizing Coproduced Climate Research as Care: Practical Lessons Learned With Women Farmland‐Owners in the Central Midwest United States

ABSTRACT This article's team of interdisciplinary researchers and conservation educator‐practitioners learned with, and from, a group of women farmland‐owners regarding how to conceptualize coproduced climate research by putting “care” at the center—care for the soil, for relationships, for data. We outline the creation and evolution of a storytelling‐based conservation program that allowed our diverse group to discover how the language of care could integrate climate analysis, conservation, and relationship‐building to foster tangible solutions. As a result of the project, the women landowners took actions that supported social‐environmental resilience—from planting cover crops to fostering watershed/neighborhood relationships. Our diverse group of women landowners and researchers had very different experiences with conservation and often very different views on climate change itself, but, through storytelling and the language of care, we not only coproduced knowledge but also created relationships and action. This article outlines specific practices for how to inflect a coproduced process for climate resilience with practices that promote care and yield action projects.

Shenk, Linda [Department of English Iowa State Uni↗
DOE OSTI2024

M IV /Co –I (M = Zr, Hf) Bis(phosphinoamide) Complexes with η 6 - and η 4 -Arenes

Despite their relevance to catalysis, low-/subvalent cobalt complexes are difficult to synthesize and isolate. Consequently, very few “cobaltate” complexes are known and there is a lack of architectural diversity in this field. Lewis acidic d 0 group IV metals have been demonstrated to stabilize Co –I centers via metal–metal bonds. Herein, we report the synthesis of bis(phosphinoamide) heterobimetallic M IV /Co –I arene complexes (M = Zr, Hf). The driving force to maintain the aromaticity of the arenes dictates the hapticity of the Co-bound arene ligands and influences the binding arrangement of the phosphinoamide ligands. Benzene and toluene were observed to bind η 6 to the Co –I center, forcing dissociation of one of the phosphinoamide ligands, whereas η 4 -coordination of anthracene allows both phosphinoamide ligands to remain bound to the Co center. The identity of both the arene and the group IV metal ion starkly influence the lability of the arene. For instance, the Co-bound benzene ligand in the Zr IV /Co –I benzene complex rapidly exchanges with C 6 D 6 in solution, whereas toluene/C 6 D 6 exchange is much slower and appreciable C 6 H 6 /C 6 D 6 exchange is not observed for the Hf analogue. The Zr IV /Co –I benzene complex loses benzene upon repeated exposure to vacuum to form an arene-free tetrametallic dimer.

aromatic compounds↗
DOE OSTI2025

Structure and ultrafast dynamics of tri-nuclear Ag-/Tl–Pt 2 POP 4 complexes in solution

The energetics and dynamics of ion assembly in solution has broad influence in nanomaterials and inorganic synthesis. To investigate the fundamental processes involved, we present a time-resolved x-ray solution scattering (TR-XSS) study of the trinuclear silver and thallium complexes of the diplatinum ion PtPOP [Pt 2 (H 2 P 2 O 5 )$_4^{4−}$] in aqueous solution. These complexes, their structural properties, and their electronic structure are not well understood and afford a unique opportunity to study the metal–metal bond formation that influences molecular and material assembly in solution. We present model-independent analysis of the observed dynamics as well as an analysis incorporating time-resolved structural refinements of key bond lengths with $<$100 fs time resolution. We find that upon photoexcitation, the Pt atoms contract ∼0.25 Å toward the center of both the Ag- and the Tl-PtPOP complexes, as previously observed for the PtPOP anion. For the AgPtPOP system, an ultrafast Ag-Pt bond expansion of ∼0.2 Å is observed, whereas in contrast, the TlPtPOP system exhibits a Tl-Pt bond contraction of ∼0.3 Å upon photoexcitation. For both complexes, the change in electronic state leads to coherent (“wave-packet”) oscillations along the metal–Pt coordinates. Based on these structural dynamics, we propose an electronic structure model that describes the metal–metal bonding behavior in both the ground and excited state for both complexes.

Lenzen, Philipp [Technical Univ. of Denmark, Lyngb↗
DOE OSTI2021

Gyrokinetic and extended-MHD simulations of a flow shear stabilized Z-pinch experiment

Axisymmetric (m = 0) gyrokinetic and extended-MHD simulations of a sheared-flow Z-pinch plasma are performed with the high-order finite volume code COGENT. The present gyrokinetic model solves the long-wavelength limit of the gyrokinetic equation for both ion and electron species coupled to the electrostatic gyro-Poisson equation for the electrostatic potential. The extended-MHD model is electromagnetic and includes the effects of the gyro-viscous pressure tensor, diamagnetic electron and ion heat fluxes, and generalized Ohm's law. Additionally, a prominent feature of this work is that the radial profiles for the plasma density and temperature are taken from the fusion Z-pinch experiment (FuZE), and the magnetic field profile is obtained as a solution of the MHD force balance equation. Such an approach allows to address realistic plasma parameters and provide insights into the current and planned experiments. In particular, it is demonstrated that the radial profiles play an important role in stabilization, as the embedded guiding center (E×B ) drift has a strong radial shear, which can contribute to the Z-pinch stabilization even in the absence of the fluid flow shear. The results of simulations for the FuZE plasma parameters show a decrease in the linear growth rate with an increase in the flow shear; however, full stabilization in the linear regime is not observed even for large (comparable to the Alfvén velocity) radial variations of the axial flow. Nonlinear stability properties of the FuZE plasmas are also studied, and it is found that profile broadening can have a pronounced stabilizing effect in the nonlinear regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Sierra/SolidMechanics Verification Report

This document presents the Sierra/SolidMechanics (Sierra/SM) verification plan. This plan centers on the tests in the Sierra/SM verification test suite, a subset of which are documented in the Sierra/SM Verification Tests Manual. Most of these tests are run nightly with the Sierra/SM code suite, and the results of the tests are checked against analytic solutions. For each of the tests presented in the Verification Tests Manual, the test setup, a description of the analytic solution, and comparison of the Sierra/SM code results to the analytic solution is provided. Mesh convergence is also checked on a nightly basis for several of these tests. This verification plan discusses these various types of tests and what they mean for Sierra/SM verification. Many other activities also contribute to Sierra/SM quality. These address code and solution quality and range from low-level unit tests, run nightly, up to full-fidelity acceptance tests, used to verify release stability. This acceptance test suite checks that new versions of Sierra/SM continue to yield the same answers for high-resolution analyst problems. Further code quality measures include an extensive suite of intermediate-size regression tests and automated nightly code quality checks. While these additional activities do not fall under a strict definition of verification, they greatly add to the quality, stability, and reliability of Sierra/SM, and are discussed here as well.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Kauhale 'O Kipapa: Kahikinui Community Cultural, Educational, and Resilience Center Preliminary Design Technical Assistance Summary Report

The Kahikinui Homestead Community in Maui faces energy resilience challenges due to its remote, off-grid location. This report summarizes preliminary technical assistance to design a microgrid for six modular homes using solar PV, wind, and battery storage. Findings show hybrid PV and battery systems, supplemented by small wind turbines, offer cost-effective, resilient solutions. This work supports enhanced energy independence and disaster preparedness for the community.

Quiroz, Jimmy Edward [Sandia National Laboratories↗
DOE OSTI2021

Interplays of electron and nuclear motions along CO dissociation trajectory in myoglobin revealed by ultrafast X-rays and quantum dynamics calculations

Significance In this multidisciplinary study, we investigated local and global structural responses upon carbon monoxide (CO) dissociation from the heme iron(II) in myoglobin in solution using ultrafast X-ray absorption and scattering. CO dissociation was optically triggered and pulsed X-rays probed structural and electronic changes as a function of time, revealing an iron(II) spin state transition and nuclear motions propagating from the active site throughout the entire protein. Quantum mechanical computational results were used to assign the iron K-edge absorption features, and to map out electronic and nuclear trajectories during Fe–CO bond elongation and heme macrocycle doming. The study has important implications for heme–protein function and the interaction of small molecule substrates with the metal active centers in enzymes.

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

Holistic energy analysis method for thermal management architectures of data centers

Modern high-performance computing (HPC) data centers (DCs), particularly those supporting energy-intensive artificial intelligence (AI) workloads, face escalating thermal management challenges that degrade performance through thermal throttling and drive up cooling power consumption and operational costs. To address this challenge, many have developed a wide variety of thermal management solutions (single-phase, two-phase, direct, indirect, hybrid, and more) which attempt to cool HPC DCs effectively while attempting to minimize overall system power consumption. However, the analysis of these solutions and methods to effectively compare one with another is lacking. Overall power usage effectiveness (PUE) and total-power usage effectiveness (TUE) provide a metric to quantify power consumption but fail to identify components in the system which require further optimization. To address this, we propose a holistic analytical framework – the waterfall diagram (WFD) – which leverages a waterfall chart methodology, offering a comprehensive visualization of both the thermal management system loop and heat flow pathways from individual server components to the outdoor ambient. Use of the WFD enables graphical estimations of power efficiency and cooling performance across each component of a DC cooling system and complements Sankey-style energy flow visualizations by additionally resolving stage-wise temperature changes and incremental TUE contributions. The framework is used in conjunction with simulation-based approaches, to conduct a detailed pressure drop and flow distribution analysis aimed at identifying the optimal coolant distribution architecture for a single-phase direct-to-chip water-cooled DC, which serves as the baseline for subsequent WFD analysis. Among the evaluated architectures, the 3 U modular coolant distribution architecture is found to demonstrate the best performance, considering minimal pressure drop and uniform flow distribution. In addition, TUE is calculated for each cooling loop component based on its associated pressure drop and corresponding pumping power, which are integrated into the WFD. This correlation between TUE and local temperature offers immediate insight into the power efficiency and thermal performance contributions of individual components, facilitating further development and optimization. Examples of WFD applications are presented under varying thermal loads and ambient conditions, demonstrating reasonable cooling strategies. Notably, the 3 U modular architecture maintains a consistent chip case temperature of 85°C, achieving a TUE of 1.016 at ambient temperature of 47°C, and a TUE of 1.026 at ambient temperature of 52°C. The WFD methodology provides an efficient, holistic, and streamlined framework for DC thermal management architecture assessment and enables design optimization which is important for addressing the thermal-fluidic energy challenges of current and next-generation DCs.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Condensed Phase Water Molecular Multipole Moments from Deep Neural Network Models Trained on Ab Initio Simulation Data

Ionic solvation phenomena in liquids involve intense interactions in the inner solvation shell. For interactions beyond the first shell, the ion–solvent interaction energies result from the sum of many smaller-magnitude contributions that can still include polarization effects. Deep neural network (DNN) methods have recently found wide application in developing efficient molecular models that maintain near-quantum accuracy. In this work, we extend the DeePMD-kit code to produce accurate molecular multipole moments in the bulk and near interfaces. The new method is validated by comparing the DNN moments with those generated by ab initio simulations. The moments are used to compute the electrostatic potential at the center of a molecular-sized hydrophobic cavity in water. The results show that the fields produced by the DNN models are in quantitative agreement with the AIMD-derived values. These efficient methods will open the door to more accurate solvation models for large solutes such as proteins.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗