Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TRAJECTORIES”

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

The hidden structure of hydrodynamic transport in random fracture networks

We study the large-scale dynamics and prediction of hydrodynamic transport in random fracture networks. The flow and transport behaviour is characterized by first passage times and displacement statistics, which show heavy tails and anomalous dispersion with a strong dependence on the injection condition. The origin of these behaviours is investigated in terms of Lagrangian velocities sampled equidistantly along particle trajectories, unlike classical sampling strategies at a constant rate. The velocity series are analysed by their copula density, the joint distribution of the velocity unit scores, which reveals a simple, albeit hidden, correlation structure that can be described by a Gaussian copula. Based on this insight, we derive a Langevin equation for the evolution of equidistant particle speeds. In this framework, particle motion is quantified by a stochastic time-domain random walk, the joint density of particle position, and speed satisfies a Klein–Kramers equation. The upscaled theory quantifies particle motion in terms of the characteristic fracture length scale and the distribution of Eulerian flow velocities. That is, it is predictive in the sense that it does not require the a priori knowledge of transport attributes. The upscaled model captures non-Fickian transport features, and their dependence on the injection conditions in terms of the velocity point statistics and average fracture length. It shows that the first passage times and displacement moments are dominated by extremes occurring at the first step. The presented approach integrates the interaction of flow and structure into a predictive model for large-scale transport in random fracture networks.

42 ENGINEERING↗
DOE OSTI2023

Chaotic and integrable magnetic fields in one-dimensional hybrid Vlasov–Maxwell equilibria

The construction of kinetic equilibrium states is important for studying stability and wave propagation in collisionless plasmas. Thus, many studies over the past decades have been focused on calculating Vlasov–Maxwell equilibria using analytical and numerical methods. However, the problem of kinetic equilibrium of hybrid models is less studied, and self-consistent treatments often adopt restrictive assumptions ruling out cases with irregular and chaotic behaviour, although such behaviour is observed in spacecraft observations of space plasmas. In this paper, we develop a one-dimensional (1-D), quasineutral, hybrid Vlasov–Maxwell equilibrium model with kinetic ions and massless fluid electrons and derive associated solutions. The model allows for an electrostatic potential that is expressed in terms of the vector potential components through the quasineutrality condition. The equilibrium states are calculated upon solving an inhomogeneous Beltrami equation that determines the magnetic field, where the inhomogeneous term is the current density of the kinetic ions and the homogeneous term represents the electron current density. We show that the corresponding 1-D system is Hamiltonian, with position playing the role of time, and its trajectories have a regular, periodic behaviour for ion distribution functions that are symmetric in the two conserved particle canonical momenta. For asymmetric distribution functions, the system is nonintegrable, resulting in irregular and chaotic behaviour of the fields. The electron current density can modify the magnetic field phase space structure, inducing orbit trapping and the organization of orbits into large islands of stability. Thus, the electron contribution can be responsible for the emergence of localized electric field structures that induce ion trapping. We also provide a paradigm for the analytical construction of hybrid equilibria using a rotating two-dimensional harmonic oscillator Hamiltonian, enabling the calculation of analytic magnetic fields and the construction of the corresponding distribution functions in terms of Hermite polynomials.

Kaltsas, Dimitrios A. (ORCID:0000000300769015)↗
DOE OSTI2023

Fast-ion transport in quasisymmetric equilibria in the presence of a resonant Alfvénic perturbation

Significant progress has been made in designing magnetic fields that provide excellent confinement of the guiding-centre trajectories of alpha particles using quasisymmetry (QS). Given the reduction in this transport channel, we assess the impact of resonant Alfvén eigenmodes (AEs) on the guiding-centre motion. The AE amplitudes are chosen to be consistent with experimental measurements and large-scale simulations. We evaluate the drift resonance condition, phase-space island width and island overlap criterion for quasisymmetric configurations. Kinetic Poincaré plots elucidate features of the transport, including stiff transport above a critical perturbation amplitude. Our analysis highlights key departures from the AE-driven transport in tokamaks, such as the avoidance of phase-space island overlap in quasihelical configurations and the enhanced transport due to wide phase-space islands in low magnetic shear configurations. In configurations that are closer to QS, with QS deviations $\delta B/B_0 \lesssim 10^{-3}$ , the transport is primarily driven by the AE, while configurations that are further from QS, $\delta B/B_0 \sim 10^{-2}$ , experience significant transport due to the QS-breaking fields in addition to the AE.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Structure and dynamics of magneto-inertial, differentially rotating laboratory plasmas

We present a detailed characterization of the structure and evolution of differentially rotating plasmas driven on the MAGPIE pulsed-power generator (1.4 MA peak current, 240 ns rise time). The experiments were designed to simulate physics relevant to accretion discs and jets on laboratory scales. A cylindrical aluminium wire array Z pinch enclosed by return posts with an overall azimuthal off-set angle was driven to produce ablation plasma flows that propagate inwards in a slightly off-radial trajectory, injecting mass, angular momentum and confining ram pressure to a rotating plasma column on the axis. However, the plasma is free to expand axially, forming a collimated, differentially rotating axial jet that propagates at ${\approx }100\,{\rm km}\,{\rm s}^{-1}$ . The density profile of the jet corresponds to a dense shell surrounding a low-density core, which is consistent with the centrifugal barrier effect being sustained along the jet's propagation. We show analytically that, as the rotating plasma accretes mass, conservation of mass and momentum implies plasma radial growth scaling as $r \propto t^{1/3}$ . As the characteristic moment of inertia increases, the rotation velocity is predicted to decrease and settle on a characteristic value ${\approx }20\,{\rm km}\,{\rm s}^{-1}$ . We find that both predictions are in agreement with Thomson scattering and optical self-emission imaging measurements.

Valenzuela-Villaseca, V. (ORCID:0000000167721441)↗
DOE OSTI2021

Formation and surface melting of nanoparticle superlattices in a solution

The wisdom in the saying of “There are no two snowflakes alike” lies in the importance of history or kinetic pathways in the phase transitions of solids. Likewise, “artificial solids,” namely superlattices consisting of functional nanoparticles, have lattice size, surface morphology, crystallinity, symmetry, and structural reconfiguration (for example, transition into a disordered state) highly dependent on the kinetic pathways as the nanoparticles interact with each other in solution [1]. Great progresses have been made in understanding the formation pathways of superlattices using liquid-phase transmission electron microscopy (TEM) [2-4]. For example, by tracking single nanoparticle’s trajectories, especially aided by U-net neural network-based machine learning, previous studies mapped the fundamental nanoparticle interactions at nanometer resolution [5]. Nonclassical, two-step nucleation pathway has also been elucidated in the system of nanoprisms, by optimizing protocols such as loading nanoparticle suspensions over the supersaturation threshold and minimizing particle‒substrate interaction [2]. Surface morphologies or exposed facets of superlattices have been shown to follow the principles of Wulff construction rule, where the facet-dependent surface energy can be measured based on the capillary wave theory [4]. However, the reverse process of crystallization of superlattices, the conversion from crystalline to disordered state, has been much less explored. On one hand, the melting of nanoparticle superlattices can provide a preferred pathway to induce structural reorganization or shuffling of building blocks for them to transform into different types of crystal structures. On the other hand, understanding nanoscale superlattice melting and comparing such behaviors with the prevailing surface melting theories developed for atomic/molecular solids can provide a potent way to engineer phase transitions of supra- and hierarchical structures constructed from nanoscale entities (e.g., DNA-coated nanoparticles, proteins), for their applications in reprogrammable and switchable materials with multifunctional properties [6, 7]. The experimental challenges to observe melting of superlattices are twofold. Practically it is difficult to load the initial superlattice form, in an intact manner, into the highly confined liquid-phase TEM chamber for in-situ observation. Here, the triggering of melting also needs meticulous manipulation of nanoparticle concentration, interparticle interaction, and solution environment.

Kim, Ahyoung↗
DOE OSTI2021

Application of Dislocation Theory to Minimize Defects in Artificial Solids Built with Nanocrystal Building Blocks

Conspectus. Oriented atomic attachment of colloidal inorganic nanocrystals represents a powerful synthetic method for preparing complex inorganic superstructures. Examples include fusion of nanocrystals into dimer and superlattice structures. If the attachment were perfect throughout, then the resulting materials would have single crystal-like alignment of the individual nanocrystals' atomic lattices. While individual colloidal nanocrystals typically are free of many defects, there are a multitude of pathways that can generate defects upon nanocrystal attachment. These attachment generated defects are typically undesirable, and thus developing strategies to favor defect-free attachment or heal defective interfaces are essential. There may also be some cases where attachment-derived defects are desirable. In this Account, we summarize our current understanding of how these defects arise, in order to offer guidance to those who are designing nanocrystal derived solids. The small size of inorganic nanocrystals means short diffusion lengths to the surface, which favor the formation of nanocrystal building blocks with pristine atomic structures. Upon attachment, however, there are numerous pathways that can lead to atomic scale defects, and bulk crystal dislocation theory provides an invaluable guide to understanding these phenomena. As an example, an atomic step edge can be incorporated into the interface leading to an extra half-plane of atoms, known as an edge dislocation. These dislocations can be well described by the Burgers vector description of dislocations, which geometrically identifies planes in which a dislocation can move. Our in situ measurements have verified that bulk dislocation theory predictions for 1D defects hold true at few-nanometer length scales in PbTe and CdSe nanocrystal interfaces. Ultimately, the applicability of dislocation theory to nanocrystal attachment enables the predictive design of attachment to prevent or facilitate healing of defects upon nanocrystal attachment. We applied similar logic to understand formation of planar (2D) defects such as stacking faults upon nanocrystal attachment. Again concepts from bulk crystal defect crystallography can identify attachment pathways that can prevent or deterministically form planar defects upon nanocrystal attachment. The concepts we discuss work well for identifying favorable attachment geometries for nanocrystal pairs ; however it is currently unclear how to translate these ideas to near-simultaneous multiparticle attachment. Geometric frustration, which prevents nanocrystal rotation, and yet-to-be considered defect generation pathways unique to multiparticle attachment complicate defect-free superlattice attachment. New imaging methods now allow for the direct observation of local attachment trajectories and may enable improved understanding of such multiparticle phenomena. With further refinement, a unified framework for understanding and ultimately eliminating structural defects in fused nanocrystal superstructures may well be achievable in coming years.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

A Miniature Multilevel Structures for Lossless Ion Manipulations Ion Mobility Spectrometer with Wide Mobility Range Separation Capabilities

Here, ion mobility spectrometry employing structures for lossless ion manipulations (SLIM-IMS) is an attractive gas-phase separation technique due to its ability to achieve unprecedented effective ion path lengths (>1 km) and IMS resolving powers in a small footprint. The emergence of multilevel SLIM technology, where ions are transferred between vertically stacked SLIM electrode surfaces, has subsequently allowed for ultralong single pass path lengths (>40 m) to be achieved, enabling ultrahigh resolution IMS measurements to be performed over the entire mobility range in a single experiment. The implementation of multiple SLIM levels requires very little additional space for the SLIM system, and we developed a miniature SLIM module (miniSLIM) based on multilevel SLIM technology and report the performance here. The module is 11.1 cm x 6.7 cm x 1.4 cm (L x W x H) and consists of three SLIM levels totaling 1-meter path length. Ion trajectory simulations were used to optimize the SLIM board spacings and SLIM board thicknesses of the miniSLIM IMS system. Methods of efficiently transferring ions between SLIM levels were examined and a new approach using asymmetric traveling waves (TWs) was developed and implemented in the miniSLIM. We experimentally characterized the performance of the --meter multilevel miniSLIM IMS-MS relative to a drift tube IM--MS using an Agilent tuning mixture and tetraalkylammonium cations. The 1-meter miniSLIM achieved a resolving power of up to 131 (CCS/ΔCCS), ~1.5x higher than achievable with a 78 cm path length drift tube IMS, and successfully transmitted the entire ion mobility range in a single separation. We also demonstrated the miniSLIM’s performance as a standalone IMS system (i.e., without MS), showing baseline separation between Agilent tuning mixture cations with the full mobility range observed, and a standard peptide mixture with different charge states readily differentiated. Overall the miniSLIM provides a compact alternative to high performance IMS instruments possessing similar path lengths.

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

Local Electric Field Effects on Water Dissociation in Bipolar Membranes Studied Using Core–Shell Catalysts

The local electric field strength is thought to affect the rate of water dissociation (WD) in bipolar membranes (BPMs) at the catalyst–nanoparticle surfaces. Here, we study core–shell nanoparticles, where the core is metallic, semiconducting, or insulating, to understand this effect. The nanoparticle cores were coated with a WD catalyst layer (TiO2 or HfO2) via atomic layer deposition (ALD), and the morphology was imaged with transmission electron microscopy. Irrespective of the core material, these core–shell catalysts displayed comparable WD overpotentials at optimal mass loading, despite the hypothesized differences in the electric field strength across the catalyst particle suggested by continuum electrostatic simulations. Substantial atomic interdiffusion between the core and shell was ruled out by X-ray absorption spectroscopy, X-ray photoelectron spectroscopy, and diffuse reflectance optical measurements. However, the optimal mass loading of catalyst was roughly 1 order of magnitude higher for the conductive and high dielectric core materials than for the low dielectric insulating cores. These findings are consistent with the hypothesis that electric field screening within the core material focuses the electric field drop between particles such that larger film thicknesses can be tolerated. Collectively, these data support the idea that it is the local electric field at the molecular level that controls proton-transfer rates and that the metal core/dielectric-shell constructs introduced here modulate that field. Further materials and synthetic design may enable optimization of the electric field strength across the proton-transfer trajectory at the material surface.

Sarma, Prasad V↗
DOE OSTI2025

Structure–Activity Relationships in Ether-Functionalized Solid-State Metal–Organic Framework Electrolytes

The structure–property relationships of metal–organic framework (MOF)-based solid-state electrolytes are not well understood. Herein, a systematic investigation of 12 Zr(IV)-based UiO-66 MOFs with varying ether-chain functional groups was carried out to elucidate the critical microscopic interactions that facilitate improved solid-state electrolyte performance. Enhanced sampling molecular dynamics (MD) simulations were employed and revealed a three-tier ion hopping mechanism: linker–linker hopping, linker-counterion hopping, and counterion-counterion hopping. Detailed structural analysis of the MD trajectories revealed that the chemistry and morphology of the linker groups affect the relative stability and population distribution of the electrolyte components, such that crown-ether-based linker groups enhance the probability of extended, low-barrier ion percolation pathways. As a result, we were able to tune the ionic conductivities by rationally manipulating the counterion distributions, linker binding strengths, and the configurational entropy (multivariability of the linkers). The resulting performance of these MOF-based solid-state electrolytes was significantly enhanced, with a methoxy-functionalized framework (UiO-66-L1 100 ) achieving high ionic conductivities of 2.32 × 10 –4 S/cm and 2.07 × 10 –3 S/cm at 30 °C and 90 °C, respectively, an order of magnitude greater than other all-solid-state MOF electrolyte systems. The electrolyte stability was evaluated with LiIn|LPSCl|MOF:LiTFSI|LPSCl|LiIn symmetric cells, showing excellent Li plating/stripping processes for over 2 months.

36 MATERIALS SCIENCE↗
DOE OSTI2020

New Concepts in Electrolytes

Over the past decades, Li-ion battery (LIB) has turned into one of the most important advances in the history of technology due to its extensive and in-depth impact on our life. Its omnipresence in all electric vehicles, consumer electronics and electric grids relies on the precisely tuned electrochemical dynamics and interactions among the electrolytes and the diversified anode and cathode chemistries therein. Furthermore, with consumers' demand for battery performance ever increasing, more and more stringent requirements are being imposed upon the established equilibria among these LIB components, and it became clear that the state-of-the-art electrolyte systems could no longer sustain the desired technological trajectory. Driven by such gap, researchers started to explore more unconventional electrolyte systems. From superconcentrated solvent-in-salt electrolytes to solid-state electrolytes, the current research realm of novel electrolyte systems has grown to unprecedented levels. In this review, we will avoid discussions on current state-of-the-art electrolytes but instead focus exclusively on unconventional electrolyte systems that represent new concepts.

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

Octane Requirements of Lean Mixed-Mode Combustion in a Direct-Injection Spark-Ignition Engine

Here, this study investigates the octane requirements of a hybrid flame propagation and controlled autoignition mode referred to as mixed-mode combustion (MMC), which allows for strong control over combustion parameters via a spark-initiated deflagration phase. Due to the throughput limitations associated with both experiments and 3-D computational fluid dynamics calculations, a hybrid 0-D and 1-D modeling methodology was developed, supported by experimental validation data. This modeling approach relied on 1-D, two-zone engine simulations to predict bulk in-cylinder thermodynamic conditions over a range of engine speeds, compression ratios, intake pressures, trapped residual levels, fueling rates, and spark timings. Those predictions were then transferred to a 0-D chemical kinetic model, which was used to evaluate the autoignition behavior of fuels when subjected to temperature–pressure trajectories of interest. Finally, the predicted autoignition phasings were screened relative to the progress of the modeled deflagration-based combustion in order to determine if an operating condition was feasible or infeasible due to knock or stability limits. The combined modeling and experimental results reveal that MMC has an octane requirement similar to modern stoichiometric spark-ignition engines in that fuels with high research octane number (RON) and high octane sensitivity (S) enable higher loads. Experimental trends with varying RON and S were well predicted by the model for 1000 and 1400 rpm, confirming its utility in identifying the compatibility of a fuel’s autoignition behavior with an engine configuration and operating strategy. However, the model was not effective in predicting (nor designed to predict) operability limits due to cycle-to-cycle variations, which experimentally inhibited operation of some fuels at 2000 rpm. Putting the operable limits and efficiency from MMC in the context of a state-of-the-art engine, the MMC showed superior efficiencies over the range investigated, demonstrating the potential to further improve fuel economy.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Food–Energy–Water Crises in the United States and China: Commonalities and Asynchronous Experiences Support Integration of Global Efforts

Food, energy, and water (FEW) systems have been recognized as an issue of critical global importance. Understanding the mechanisms that govern the FEW nexus is essential to develop solutions and avoid humanitarian crises of displacement, famine, and disease. The U.S. and China are the world’s leading economies. Although the two nations are shaped by fundamentally different political and economic systems, they share FEW trajectories in several complementary ways. These realities place the U.S. and China in unique positions to engage in problem definition, dialogue, actions, and transdisciplinary convergence of research to achieve productive solutions addressing FEW challenges. By comparing the nexus and functions of the FEW systems in the two nations, this perspective aims to facilitate collaborative innovations that reduce disciplinary silos, mitigate FEW challenges, and enhance environmental sustainability. Finally, the review of experiences and challenges facing the U.S. and China also offers valuable insights for other nations seeking to achieve sustainable development goals.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Physicochemical and Molecular Insights into the Boundary Layer and Free Troposphere Aerosol Interactions over the Southern Great Plains

Ambient aerosols’ vertical profiles are critical for evaluating the role of aerosols in atmospheric chemistry and radiative transfer, but limited data on these profiles hinders our ability to fully assess their impact on the Earth's radiative balance. Here, in this study, we investigated the size-, time-, and altitude resolved composition of individual particles and bulk molecular composition of particle samples collected by an uncrewed aerial system–ArcticShark over the Southern Great Plains. Single particle microanalysis shows that, the free tropospheric (FT) samples are dominated (56-66%) by carbonaceous sulfate particles, while boundary layer (BL) samples are dominated (57-74%) by carbonaceous particles. Back trajectory simulations suggest that FT particles are likely influenced by long-range transport and have undergone aqueous-phase processing. Conversely, in-situ size distribution data shows evidence of particle growth in the upper BL and just below the FT. This observation may indicate vertical transport of particles from an elevated aerosol layer in the FT, possibly linked to a new particle formation event. This observation is further supported by high resolution molecular composition data, which reveals particle volatility increasing with increasing size, which aligns with the growth event. This study aids in fundamental understanding of the compositional and molecular specificity of vertically resolved organic aerosols to provide insights into particle size evolution for future atmospheric models.

ArcticShark↗
DOE OSTI2026

Three-Dimensional Settling Dynamics of Environmental Microplastics

Resolving the three-dimensional settling dynamics of microplastic (MP) particles is essential for developing comprehensive models of MP transport in rivers-both vertically within the water column and laterally across the channel. While previous research has largely examined one-dimensional vertical settling velocities, little is known about the lateral drifting, settling paths, and horizontal velocities of MPs. To address this, we investigated the full three-dimensional settling behavior of environmental MPs collected from rivers and ocean water, as well as from estuarine and ocean sediment. Geometric properties of 127 environmental MPs were quantified by a dynamic particle image analyzer, and their settling trajectories were recorded and reconstructed via a multicamera tracking algorithm. This enabled quantification of a particle's horizontal drift, tortuosity, amplitude and settling pattern, as well as vertical and horizontal velocities. Results showed that spherical MPs settled with minimal lateral displacement, whereas elongated particles, such as rod- and blade-shaped MPs, displayed pronounced lateral movements, reaching up to 65 times their equivalent diameter and averaging more than twice that of spheres. These dynamics suggest that elongated MPs may have a greater probability for wider lateral dispersion in rivers, increasing their likelihood for interactions with riverbanks and channel boundaries compared to more spherical shaped MPs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Coordination Sphere of Lanthanide Aqua Ions Resolved with Ab Initio Molecular Dynamics and X-ray Absorption Spectroscopy

To resolve the fleeting structures of lanthanide Ln 3+ aqua ions in solution, we (i) performed the first ab initio molecular dynamics (AIMD) simulations of the entire series of Ln 3+ aqua ions in explicit water solvent using pseudopotentials and basis sets recently optimized for lanthanides and (ii) measured the symmetry of the hydrating waters about Ln 3+ ions (Nd 3+ , Dy 3+ , Er 3+ , Lu 3+ ) for the first time with extended X-ray absorption fine structure (EXAFS). EXAFS spectra were measured experimentally and generated from AIMD trajectories to directly compare simulation, which concurrently considers the electronic structure and the atomic dynamics in solution, with experiment. In this study, we performed a comprehensive evaluation of EXAFS multiple-scattering analysis (up to 6.5 Å) to measure Ln–O distances and angular correlations (i.e., symmetry) and elucidate the molecular geometry of the first hydration shell. This evaluation, in combination with symmetry-dependent L 3 - and L 1 -edge spectral analysis, shows that the AIMD simulations remarkably reproduces the experimental EXAFS data. The error in the predicted Ln–O distances is less than 0.07 Å for the later lanthanides, while we observed excellent agreement with predicted distances within experimental uncertainty for the early lanthanides. Our analysis revealed a dynamic, symmetrically disordered first coordination shell, which does not conform to a single molecular geometry for most lanthanides. This work sheds critical light on the highly elusive coordination geometry of the Ln 3+ aqua ions.

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

Ensemble Simulations and Experimental Free Energy Distributions: Evaluation and Characterization of Isoxazole Amides as SMYD3 Inhibitors

Optimization of binding affinities for ligands to their target protein is a primary objective in rational drug discovery. Herein, we report on a collaborative study that evaluates various compounds designed to bind to the SET and MYND domain-containing protein 3 (SMYD3). SMYD3 is a histone methyltransferase and plays an important role in transcriptional regulation in cell proliferation, cell cycle, and human carcinogenesis. Experimental measurements using the scintillation proximity assay show that the distributions of binding free energies from a large number of independent measurements exhibit non-normal properties. We use ESMACS (enhanced sampling of molecular dynamics with approximation of continuum solvent) and TIES (thermodynamic integration with enhanced sampling) protocols to predict the binding free energies and to provide a detailed chemical insight into the nature of ligand-protein binding. Our results show that the 1-trajectory ESMACS protocol works well for the set of ligands studied here. Although one unexplained outlier exists, we obtain excellent statistical ranking across the set of compounds from the ESMACS protocol and good agreement between calculations and experiments for the relative binding free energies from the TIES protocol. ESMACS and TIES are again found to be powerful protocols for the accurate comparison of the binding free energies.

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

Ice Phase Classification Made Easy with Score-Based Denoising

Accurate identification of ice phases is essential for understanding various physicochemical phenomena. However, such classification for structures simulated with molecular dynamics is complicated by the complex symmetries of ice polymorphs and thermal fluctuations. For this purpose, both traditional order parameters and data-driven machine learning approaches have been employed, but they often rely on expert intuition, specific geometric information, or large training data sets. In this work, we present an unsupervised phase classification framework that combines a score-based denoiser model with a subsequent model-free classification method to accurately identify ice phases. Further, the denoiser model is trained on perturbed synthetic data of ideal reference structures, eliminating the need for large data sets and labeling efforts. The classification step utilizes the smooth overlap of atomic position (SOAP) descriptors as the atomic fingerprint, ensuring Euclidean symmetries and transferability to various structural systems. Our approach achieves a remarkable 100% accuracy in distinguishing ice phases of test trajectories using only seven ideal reference structures of ice phases as model inputs. This demonstrates the generalizability of the score-based denoiser model in facilitating phase identification for complex molecular systems. The proposed classification strategy can be broadly applied to investigate structural evolution and phase identification for a wide range of materials, offering new insights into the fundamental understanding of water and other complex systems.

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

Quantum-Based Molecular Dynamics Simulations Using Tensor Cores

Tensor cores, along with tensor processing units, represent a new form of hardware acceleration specifically designed for deep neural network calculations in artificial intelligence applications. Tensor cores provide extraordinary computational speed and energy efficiency but with the caveat that they were designed for tensor contractions (matrix–matrix multiplications) using only low-precision floating-point operations. Despite this perceived limitation, we demonstrate how tensor cores can be applied with high efficiency to the challenging and numerically sensitive problem of quantum-based Born–Oppenheimer molecular dynamics, which requires highly accurate electronic structure optimizations and conservative force evaluations. The interatomic forces are calculated on-the-fly from an electronic structure that is obtained from a generalized deep neural network, where the computational structure naturally takes advantage of the exceptional processing power of the tensor cores and allows for high performance in excess of 100 Tflops on a single Nvidia A100 GPU. Stable molecular dynamics trajectories are generated using the framework of extended Lagrangian Born–Oppenheimer molecular dynamics, which combines computational efficiency with long-term stability, even when using approximate charge relaxations and force evaluations that are limited in accuracy by the numerically noisy conditions caused by the low-precision tensor core floating-point operations. A canonical ensemble simulation scheme is also presented, where the additional numerical noise in the calculated forces is absorbed into a Langevin-like dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗