Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “core crystallization”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Integral X-ray diffuse scattering for studying irradiation-induced dislocation loops in single crystals

Numerical diffuse scattering cross-section calculations are used to establish a rigorous basis for determining the concentration and size distribution of dis­location loops in irradiated single crystals from integral X-ray diffuse scattering (XRDS) measurements. Differential XRDS intensities for prismatic {111} type dislocation loops are numerically calculated as a function of loop radius R and wavevectors q relative to Bragg reflections in tungsten. The results show the well known 1/q 2 Huang scattering form at small q that transitions to a ∼1/q 4 dependence associated with the Stokes–Wilson approximation for q ≳ 1/R. More importantly, they show further that the 1/q 4 falloff is not the asymptotic large-q form of the diffuse scattering for small loops (R < 200 Å) as has often been assumed. Rather, for loop sizes as small as R ≃ 5 Å with strong curvature, the calculations show definitively that the scattering transitions to a robust 1/q 5 falloff at larger q that arises due to the local strains near the dislocation core defining the circumference of the dislocation loops. The presence of this 1/q 5 asymptotic form for both small and large loops is experimentally confirmed using an integral XRDS measurement around the 110 reflection on self-ion-irradiated tungsten combined with numerically calculated integral XRDS cross-sections. Accordingly, the historical two-region theoretical treatment of the cross-sections for integral XRDS is extended to a three-region model that has direct sensitivity to the (first-moment) dislocation line lengths of dislocation loops. In conclusion, these developments enable the use of both numerical and analytically modeled cross-sections to make accurate integral XRDS determinations of dislocation loop sizes and concentrations using modest-intensity laboratory X-ray sources.

42 ENGINEERING↗
NASA NTRS2014

Modeling the Deterioration of Engine and Low Pressure Compressor Performance During a Roll Back Event Due to Ice Accretion

The main focus of this study is to apply a computational tool for the flow analysis of the engine that has been tested with ice crystal ingestion in the Propulsion Systems Laboratory (PSL) of NASA Glenn Research Center. A data point was selected for analysis during which the engine experienced a full roll back event due to the ice accretion on the blades and flow path of the low pressure compressor. The computational tool consists of the Numerical Propulsion System Simulation (NPSS) engine system thermodynamic cycle code, and an Euler-based compressor flow analysis code, that has an ice particle melt estimation code with the capability of determining the rate of sublimation, melting, and evaporation through the compressor blade rows. Decreasing the performance characteristics of the low pressure compressor (LPC) within the NPSS cycle analysis resulted in matching the overall engine performance parameters measured during testing at data points in short time intervals through the progression of the roll back event. Detailed analysis of the fan-core and LPC with the compressor flow analysis code simulated the effects of ice accretion by increasing the aerodynamic blockage and pressure losses through the low pressure compressor until achieving a match with the NPSS cycle analysis results, at each scan. With the additional blockages and losses in the LPC, the compressor flow analysis code results were able to numerically reproduce the performance that was determined by the NPSS cycle analysis, which was in agreement with the PSL engine test data. The compressor flow analysis indicated that the blockage due to ice accretion in the LPC exit guide vane stators caused the exit guide vane (EGV) to be nearly choked, significantly reducing the air flow rate into the core. This caused the LPC to eventually be in stall due to increasing levels of diffusion in the rotors and high incidence angles in the inlet guide vane (IGV) and EGV stators. The flow analysis indicating compressor stall is substantiated by the video images of the IGV taken during the PSL test, which showed water on the surface of the IGV flowing upstream out of the engine, indicating flow reversal, which is characteristic of a stalled compressor.

Veres, Joseph P.↗
NASA NTRS2014

Modeling the Deterioration of Engine and Low Pressure Compressor Performance During a Roll Back Event Due to Ice Accretion

The main focus of this study is to apply a computational tool for the flow analysis of the engine that has been tested with ice crystal ingestion in the Propulsion Systems Laboratory (PSL) of NASA Glenn Research Center. A data point was selected for analysis during which the engine experienced a full roll back event due to the ice accretion on the blades and flow path of the low pressure compressor. The computational tool consists of the Numerical Propulsion System Simulation (NPSS) engine system thermodynamic cycle code, and an Euler-based compressor flow analysis code, that has an ice particle melt estimation code with the capability of determining the rate of sublimation, melting, and evaporation through the compressor blade rows. Decreasing the performance characteristics of the low pressure compressor (LPC) within the NPSS cycle analysis resulted in matching the overall engine performance parameters measured during testing at data points in short time intervals through the progression of the roll back event. Detailed analysis of the fan-core and LPC with the compressor flow analysis code simulated the effects of ice accretion by increasing the aerodynamic blockage and pressure losses through the low pressure compressor until achieving a match with the NPSS cycle analysis results, at each scan. With the additional blockages and losses in the LPC, the compressor flow analysis code results were able to numerically reproduce the performance that was determined by the NPSS cycle analysis, which was in agreement with the PSL engine test data. The compressor flow analysis indicated that the blockage due to ice accretion in the LPC exit guide vane stators caused the exit guide vane (EGV) to be nearly choked, significantly reducing the air flow rate into the core. This caused the LPC to eventually be in stall due to increasing levels of diffusion in the rotors and high incidence angles in the inlet guide vane (IGV) and EGV stators. The flow analysis indicating compressor stall is substantiated by the video images of the IGV taken during the PSL test, which showed water on the surface of the IGV flowing upstream out of the engine, indicating flow reversal, which is characteristic of a stalled compressor.

Jorgenson, Philip C. E.↗
NASA NTRS1993

High-resolution X-ray spectroscopy of the supernova remnant N132D

A joint nonequilibrium ionization analysis of spectral data from the Einstein Observatory of the SNR N132D in the LMC is presented on the basis of data from the Focal Plane Crystal Spectrometer (FPCS) and the Solid State Spectrometer (SSS), and lower spectral resolution data from the IPC and the Monitor Proportional Counter (MPC). The FPCS detected individual emission lines of O VII, O VIII, Ne IX, Ne X, Fe XVII, and possibly Fe XX. Measured line widths for the oxygen lines suggest Doppler broadening that is roughly consistent with optically measured expansion velocities of 2250 km/s. At the SSS/IPC temperature, FPCS flux ratios constrain the O/Fe abundance to be at least 1.9 times its solar value and the O/Ne abundance to be 0.2-1.0 times its solar value. Models for remnants with progenitor masses of 20 and 25 solar masses are completely consistent with the data, while remnants with progenitor masses of 13 and 15 solar masses can be made consistent if the progenitors are required to eject a large fraction of their iron cores.

Hwang, Una↗
DOE OSTI2026

Elongation Factor Tu Acts as a Chaperone to Activate an Antibacterial RNase Toxin

Many Gram-negative bacterial species use contact-dependent growth inhibition (CDI) systems to deliver toxic proteins into neighboring competitors. CDI + strains deploy CdiA effector proteins, which translocate their C-terminal toxin (CT) domains into target bacteria through a receptor-mediated delivery pathway. To protect against auto-intoxication, CDI + bacteria also produce CdiI immunity proteins that neutralize CT toxin activity. Here, we present the crystal structure of the CT·CdiI O32:H37 complex from Escherichia coli O32:H37. CT O32:H37 adopts the same fold as the tRNase domain of colicin D, and the nucleases share similar catalytic centers. However, unlike colicin D, which cleaves the anticodon loops of tRNA Arg isoacceptors, CT O32:H37 exhibits nonspecific RNase activity. Notably, we find that endogenous elongation factor Tu (EF-Tu) co-purifies with the over-produced CT·CdiI O32:H37 complex. Although EF-Tu does not bind stably to CT O32:H37 in the absence of CdiI O32:H37 , the translation factor is required for toxic RNase activity in vitro. AlphaFold 3 modeling and site-directed mutagenesis indicate that CT O32:H37 interacts with the N-terminal GTPase domain of EF-Tu. EF-Tu appears to stabilize residue Trp52 within the hydrophobic core of the toxin, which in turn supports the RNase active site through an unusual hydrogen-bonding interaction with the catalytic His67 residue. Furthermore, EF-Tu is hijacked as an essential co-factor to organize the toxin's catalytic center.

Nhan, Dinh Quan [University of California, Santa B↗
DOE OSTI2021

Experimental investigation on the coolability of nuclear reactor debris beds using seawater

During the Fukushima nuclear reactor accident, seawater was injected into the reactor core to cool the decay heat from the heat generating porous debris bed. However, the impact of dissolved salts in water on the coolability of a debris bed is not well understood. This paper is one of the first works to investigate the cooling mechanisms for debris bed using seawater. An experimental system was built where the packed bed test section was volumetrically heated via a high frequency induction heater. The temperature along the axis of the test section was measured using a high resolution temperature sensor based on Optical Frequency Domain Reflectometry. Experimental results show that the dryout heat flux increases with an increase in concentration, with pure water having a dryout volumetric heat flux of 1.8 MW/m and salt water of concentration 7% having a dryout volumetric heat flux of 2.8 MW/m. However, this enhanced coolability using salt water is short term, as the crystallization fouling mechanism becomes significant over time and leads to gradual temperature excursion in the bed. The fouling mechanism eventually leads to a cross section of the debris bed to plug with the deposited salts, and this leads to a pressure build up in the lower regions of the debris bed.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS1993

Modeling study on the cleavage step of the self-splicing reaction in group I introns

A three-dimensional model of the Tetrahymena thermophila group I intron is used to further explore the catalytic mechanism of the transphosphorylation reaction of the cleavage step. Based on the coordinates of the catalytic core model proposed by Michel and Westhof (Michel, F., Westhof, E. J. Mol. Biol. 216, 585-610 (1990)), we first converted their ligation step model into a model of the cleavage step by the substitution of several bases and the removal of helix P9. Next, an attempt to place a trigonal bipyramidal transition state model in the active site revealed that this modified model for the cleavage step could not accommodate the transition state due to insufficient space. A lowering of P1 helix relative to surrounding helices provided the additional space required. Simultaneously, it provided a better starting geometry to model the molecular contacts proposed by Pyle et al. (Pyle, A. M., Murphy, F. L., Cech, T. R. Nature 358, 123-128. (1992)), based on mutational studies involving the J8/7 segment. Two hydrated Mg2+ complexes were placed in the active site of the ribozyme model, using the crystal structure of the functionally similar Klenow fragment (Beese, L.S., Steitz, T.A. EMBO J. 10, 25-33 (1991)) as a guide. The presence of two metal ions in the active site of the intron differs from previous models, which incorporate one metal ion in the catalytic site to fulfill the postulated roles of Mg2+ in catalysis. The reaction profile is simulated based on a trigonal bipyramidal transition state, and the role of the hydrated Mg2+ complexes in catalysis is further explored using molecular orbital calculations.

NASA Discipline Exobiology↗
DOE OSTI2020

A Chip-Scale Oscillation-Mode Optomechanical Inertial Sensor Near the Thermodynamical Limits

Modern navigation systems integrate the global positioning system (GPS) with an inertial navigation system (INS), which complement each other for correct attitude and velocity determination. Here, the core of the INS integrates accelerometers and gyroscopes used to measure forces and angular rate in the vehicular inertial reference frame. With the help of gyroscopes and by integrating the acceleration to compute velocity and distance, precision and compact accelerometers with sufficient accuracy can provide small–error location determination. Solid–state implementations, through coherent readout, can provide a platform for high performance acceleration detection. In contrast to prior accelerometers using piezoelectric or capacitive readout techniques, optical readout provides narrow–linewidth high–sensitivity laser detection along with low–noise resonant optomechanical transduction near the thermodynamical limits. Here an optomechanical inertial sensor with an 8.2 µg Hz –1/2 velocity random walk (VRW) at an acquisition rate of 100 Hz and 50.9 µg bias instability is demonstrated, suitable for applications, such as, inertial navigation, inclination sensing, platform stabilization, and/or wearable device motion detection. Driven into optomechanical sustained–oscillation, the slot photonic crystal cavity provides radio–frequency readout of the optically–driven transduction with an enhanced 625 µg Hz –1 sensitivity. Measuring the optomechanically–stiffened oscillation shift, instead of the optical transmission shift, provides a 220× VRW enhancement over pre–oscillation mode detection.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Topological photonic crystal fiber

Photonic crystal fibers (PCFs) provide a versatile platform for various applications, thanks to the flexibility with which light guiding can be customized by modifying the fiber geometry. We realize a PCF with guided modes produced by photonic band structure topology rather than conventional mode-trapping mechanisms. The design, which is compatible with the stack-and-draw fabrication process, consists of a cross-sectional photonic topological crystalline insulator with a disclination. A bulk-defect correspondence produces degenerate topological modes, lying below the cladding light line. We use various theoretical methods to confirm their topological origins, including a spectral localizer that makes minimal assumptions about the band structure. Our experiments on the fabricated fiber show it transmitting visible to near-infrared light with low losses of 10 to 20 decibels per kilometer, which do not increase substantially when the fiber is bent. A comparable solid-core PCF of conventional design exhibits substantially higher bending losses. Optical fibers based on topological modes hold promise for improved performance and versatile functionalities.

Zhu, Bofeng [Nanyang Technological Univ. (Singapor↗
DOE OSTI2022

Experimental X-ray Charge-Density Studies–A Suitable Probe for Superconductivity? A Case Study on MgB 2

Case studies of 1T-TiSe 2 and YBa 2 Cu 3 O 7-δ have demonstrated that X-ray diffraction (XRD) studies can be used to trace even subtle structural phase transitions which are inherently connected with the onset of superconductivity in these benchmark systems. However, the utility of XRD in the investigation of superconductors like MgB 2 lacking an additional symmetry-breaking structural phase transition is not immediately evident. Nevertheless, high-resolution powder XRD experiments on MgB 2 in combination with maximum entropy method analyses hinted at differences between the electron density distributions at room temperature and 15 K, that is, below the T c of approx. 39 K. The high-resolution single-crystal XRD experiments in combination with multipolar refinements presented here can reproduce these results but show that the observed temperature-dependent density changes are almost entirely due to a decrease of atomic displacement parameters as a natural consequence of a reduced thermal vibration amplitude with decreasing temperature. Our investigations also shed new light on the presence or absence of magnesium vacancies in MgB 2 samples–a defect type claimed to control the superconducting properties of the compound. Here, we propose that previous reports on the tendency of MgB 2 to form non-stoichiometric Mg 1–x B 2 phases (1 – x ~ 0.95) during high-temperature (HT) synthesis might result from the interpretation of XRD data of insufficient resolution and/or usage of inflexible refinement models. Indeed, advanced refinements based on an Extended Hansen–Coppens multipolar model and high-resolution X-ray data, which consider explicitly the contraction of core and valence shells of the magnesium cations, do not provide any significant evidence for the formation of non-stoichiometric Mg 1–x B 2 phases during HT synthesis.

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

Nickel-rich, volatile depleted iron meteorites: Relationships and formation processes

Ungrouped iron meteorites Tishomingo, Willow Grove, and Chinga, and group IVB iron meteorites, are Ni-rich. Similarities include enrichments of 10–100 × CI for some refractory siderophile elements, and equivalent depletions in more volatile siderophile elements. Superimposed on the overall enrichment/depletion trend, certain siderophile elements (P, W, Fe, Mo) are depleted relative to elements of similar volatility. All three ungrouped irons derive from parent bodies formed in the early Solar System. Willow Grove and Chinga are characterized by cosmic ray exposure corrected 182 W/ 184 W consistent with metal-silicate segregation on their parent bodies within 1–3 Myr of Solar System formation, within the age range determined for segregation of magmatic iron meteorite parent bodies, including group IVB irons. Tishomingo is characterized by a younger model age 4–5 Myr subsequent to Solar System formation, reflecting either late stage melting resulting from 26 Al decay, or an impact resetting. The discovery of stishovite in Tishomingo, indicating exposure to a minimum shock pressure of 8–9 GPa, is consistent with the latter. Stishovite in Tishomingo and chromite included in troilite-daubréelite in IVB irons allows oxygen isotopic composition comparison between these meteorites. Different mass independent oxygen isotopic compositions of IVB irons and Tishomingo indicate genetically distinct parent bodies. By contrast, mass independent Mo isotopic compositions overlap within analytical uncertainties, indicating a similar, carbonaceous chondrite (CC) type genetic heritage. Molybdenum and 183 W isotopic data for Chinga and Willow Grove indicate derivation from CC type parent bodies. Willow Grove shares Mo and 183 W isotopic compositions with the Ni-rich South Byron Trio (SBT) grouplet and the Milton pallasite. These Ni-rich meteorites likely formed in the same general nebular environment as other CC planetesimals, likely the outer Solar System. Highly siderophile element (HSE) abundances of Willow Grove and Tishomingo are similar to some IVB meteorites, consistent with formation by moderate degrees of fractional crystallization from initial metallic melts with low S and P, and modestly fractionated HSE. The comparable HSE abundances of Tishomingo and Willow Grove to some IVB irons, yet substantially higher Ni concentrations, indicate formation on parent bodies with lower bulk HSE abundances or HSE concentration in proportionally smaller volumes of metal. HSE abundances in Chinga are considerably lower than in IVB irons, highly fractionated, and processes responsible for these remain elusive. For IVB irons and these ungrouped irons, high temperature condensation likely dominated the enrichment and depletion of the refractory and volatile siderophile elements, respectively. Parent body degassing may have also played a role. Relative depletion of volatile siderophile elements is not, however, a universal feature of high-Ni meteorites. The SBT and Milton pallasite are Ni-rich, but less depleted in the more volatile siderophile elements. Nickel enrichment was likely driven by oxidation of Fe metal during parent body accretion or core segregation. Oxidation of the Tishomingo and Willow Grove parent bodies may have occurred at ~ IW + 1, indicated by relative Mo and W depletions due to metal/water reaction during differentiation. Late-stage reduction, indicated by the presence of Cr-bearing sulfides in Tishomingo and IVB irons, may have resulted from exhaustion of the oxidant.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Prediction and Control of Atomic Ordering in Electrodeposited Binary Alloy Films: Direct Synthesis of L10 Magnetic Phases

Permanent magnet materials are essential components in renewable energy technology, underlying, electric vehicles, hybrids, wind turbines and more. Today commercial permanent magnets based on Nd-Fe-B exhibit an energy product of ~ 400 kJ/m 3 , capable to run wind turbines efficiently. Unfortunately, the supplies of rare-earth elements are centralized, and their availability may see volatility over time. This state of affairs stimulates the search for the search for rare-earth-free permanent magnetic materials. The L1 0 crystal structure underlines an important class of chemical ordered alloys that exhibit uniaxial magnetocrystalline anisotropy. This material has been considered an excellent candidate for rare-earth-free permanent magnets. In this program, a series of systematic experiments were performed and analysed in order to understand the effect of each. In order to synthesize L 10-Fe-Ni, one has to circumvent the sluggish kinetic and the limited driving force resulted from the low order-disorder temperature (~320°C) of the system. A strategy was found to trigger a martensitic transformation from the FCC phase to one of the non-cubic phases by high-strain methods. The body-centered tetragonal Fe-Ni obtained by this strategy was considered as a precursor of L10 Fe-Ni. The formation of BCT Fe-Ni via a displacive process was also implemented at the nanoscale by using tailored AuNi@FeNi core@shell.

36 MATERIALS SCIENCE↗
NASA NTRS1994

Hydrogen storage in Earth's mantle and core

Two different approaches to explaining how hydrogen might be stored in the mantle are illustrated by a number of papers published over the past 25-30 years, but there has been little attempt to provide objective comparisons of the two. One approach invokes the presence in the mantle of dense hydrous magnesium silicates (DHMS) stable at elevated pressures and temperatures. The other involves nominally anhydrous minerals (NAM) that contain hydrogen as a minor constituent on the ppm level. Experimental studies on DHMS indicate these phases may be stable to pressures and temperatures as high at 16 GPa and 1200 C. This temperature is lower than that indicated by a mantle geotherm at 16 GPa, but may be reasonable for a subducting slab. It is possible that other DHMS could be stable to even higher pressures, but little is known about maximum temperature limits. For NAM, small amounts of hydrogen (up to several hundred ppm) have been detected in olivine, orthopyroxene, clinopyroxene, and garnet recovered from xenoliths in kimberlites, eclogites, and alkali basalts; it has been demonstrated that synthetic wadsleyite and perovskite can accommodate significant amounts of hydrogen. A number of problems are associated with each possibility. For NAM originating in the mantle, one would like to assume that the hydrogen measured in samples recovered on Earth's surface was incorporated when the phase-crystallized at high temperatures and pressures, but it could have been introduced during transport to the surface. Major problems for the DHMS proponents are that none of these phases have been found as minerals and little is yet known about their stabilities in systems containing other cations such as Fe, Al, and Ca.

Prewitt, Charles T.↗
DOE OSTI2023

Polyhexatic and polycrystalline states of skyrmion lattices

We report Monte Carlo studies of lattices of up to 10 5 skyrmions treated as particles with negative core energy and repulsive interaction obtained from a microscopic spin model. Temperature dependence of translational and orientational correlations has been investigated for different experimental protocols and initial conditions. Cooling the skyrmion liquid from a fully disordered high-temperature state results in the formation of a skyrmion polycrystal. A perfect skyrmion lattice prepared at T = 0, on raising temperature undergoes a first-order melting transition into a polyhexatic state that consists of large orientationally ordered domains of fluctuating shape. On the further increasing temperature, these domains decrease in size, leading to a fully disordered liquid of skyrmions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Reducing waste of the hydride-dehydride process for U-6 wt% Nb spherical powders through lower impact and targeted milling

The breakdown of solid metal into powder during the hydride-dehydride process is commercially important for the formation of titanium and other metal powders. Typically, the milling of the brittle hydride powder occurs in a ball mill, where milling media impacts powder particles to break them down. The milling media can impact particles that are larger than desired, as desired, or smaller than desired, indiscriminately making all particles smaller. In this work, we investigate how to minimize waste powder production during milling using two different milling methods, planetary ball milling and milling in a sieve shaker (sieve-milling). Both processes yielded similar amounts of 20–75 μm diameter powder (the target size range); however, sieve-milling generated a significantly smaller amount of undersized waste powder. The powders were characterized by X-ray diffraction, SEM/STEM, and magnetic susceptibility. Several differences between ball milling and sieve-milling processes are discussed. We then conclude that the decreased yield of undersized powder in sieve-milling was due to a combination of lower impact energy in sieve-milling, unreacted metallic cores in the hydride flakes, and the ability to mill target particle sizes during sieve-milling. While these results are from the milling of brittle hydride powder, similar methods may be applicable to other brittle powders, including ceramics or salts.

Chemistry↗
NASA NTRS2019

A Coordinated Microstructural and Isotopic Study ofa Wark-Lovering Rim on a Vigarano CAI

We carried out a coordinated mineralogical and isotopic study of a Wark-Lovering (WL) rim on a Ca,Al-rich inclusion (CAI) from the reduced CV3 chondrite Vigarano. The outermost edge of the CAI mantle is mineralogically and texturally distinct compared to the underlying mantle that is composed of coarse, zoned melilite (Åk~10-60) grains. The mantle edge contains fine-grained gehlenite with hibonite and rare grossite that likely formed by rapid crystallization from a melt enriched in Ca and Al. These gehlenite and hibonite layers are surrounded by successive layers of spinel, zoned melilite (Åk~0-10), zoned diopside that grades outwards from Al,Ti-rich to Al,Ti-poor, and forsteritic olivine intergrown with diopside. These layered textures are indicative of sequential condensation of spinel, melilite, diopside, and forsterite onto hibonite. Anorthite occurs as a discontinuous layer that corrodes adjacent melilite and Al2 diopside, and appears to have replaced them, probably even later than the forsterite layer formation. Based on these observations, we conclude that the WL rim formation was initiated by flash melting and extensive evaporation of the original inclusion edge, followed by subsequent gas-solid reactions under highly dynamic conditions. All the WL rim minerals are 16O-rich (Δ17O = ~−23‰), indicating their formation in an 16O-rich nebular reservoir. Our Al-Mg measurements of hibonite, spinel, and diopside from the WL rim, as well as spinel and Al,Ti-diopside in the core, define a single, well-correlated isochron with an inferred initial 26Al/27Al ratio of (4.94 ± 0.12) × 10−5. This indicates that the WL rim formed shortly after the host CAI. In contrast, the lack of 26Mg excesses in the WL rim anorthite suggest its later formation or later isotopic disturbance in the solar nebula, after 26Al had decayed.

Jangmi Han↗
DOE OSTI2022

A deep potential model with long-range electrostatic interactions

Machine learning models for the potential energy of multi-atomic systems, such as the deep potential (DP) model, make molecular simulations with the accuracy of quantum mechanical density functional theory possible at a cost only moderately higher than that of empirical force fields. However, the majority of these models lack explicit long-range interactions and fail to describe properties that derive from the Coulombic tail of the forces. To overcome this limitation, we extend the DP model by approximating the long-range electrostatic interaction between ions (nuclei + core electrons) and valence electrons with that of distributions of spherical Gaussian charges located at ionic and electronic sites. The latter are rigorously defined in terms of the centers of the maximally localized Wannier distributions, whose dependence on the local atomic environment is modeled accurately by a deep neural network. In the DP long-range (DPLR) model, the electrostatic energy of the Gaussian charge system is added to short-range interactions that are represented as in the standard DP model. The resulting potential energy surface is smooth and possesses analytical forces and virial. Missing effects in the standard DP scheme are recovered, improving on accuracy and predictive power. By including long-range electrostatics, DPLR correctly extrapolates to large systems the potential energy surface learned from quantum mechanical calculations on smaller systems. We illustrate the approach with three examples: the potential energy profile of the water dimer, the free energy of interaction of a water molecule with a liquid water slab, and the phonon dispersion curves of the NaCl crystal.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗