Search NASASearch

SEARCH · Search NASA

Results for “high-throughput calculations”

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.

740 records · Page 3

Reduced fast-ion transport calculations of m = n = 1 fishbone-like instabilities in MAST-U

Fast-ion transport associated with an m = n = 1 fishbone-like burst in MAST-U discharge 47128 is investigated using a reduced guiding-center-based transport model (ORBIT-Kick) constrained by multi-diagnostic measurements. The two-dimensional beam-emission spectroscopy system provides measurements of the core poloidal mode structure and fluctuation amplitude, while EFIT++ reconstructions constrained by the motional Stark effect diagnostic indicate a flat q-profile with q 0 > 1⁠, indicating the absence of a resonant q = 1 surface and supporting a pressure-driven infernal-mode interpretation. Analytic m = n = 1 displacement profiles consistent with the measured core mode structure and equilibrium constraints are used as the mode structure inputs to ORBIT-Kick. The calculations show that the dominant resonances occur between the mode and co-passing fast ions, producing redistribution localized near the magnetic axis. Synthetic neutron camera signals from TRANSP-Kick recover up to 90% of the experimentally observed neutron deficit at the time of peak mode amplitude, indicating that the measured m = n = 1 mode is a dominant contributor to core fast-ion transport. However, the synthetic neutron signals recover rapidly, whereas the measured neutron emission continues to decrease after the peak amplitude. In conclusion, the remaining discrepancy may arise from contributions not included in the present single-harmonic model, including higher-m and higher-n harmonics, multi-harmonic interactions, and additional transport mechanisms, motivating future diagnostic development and modeling efforts to resolve and incorporate these additional contributions.

Wong, Henry H. [University of California, Los Ange

Equation of state for Hf, Ta, W, Re, Os, Ir, Pt, and Au to multi-terapascal pressures from density-functional theory

We present the zero-temperature equation of state (pressure dependence of compression) and phase stability predictions for the 5d-transition metals obtained from all-electron density-functional theory (DFT) calculations. The results compare favorably with experiments but extend beyond current experimental capabilities to 10 TPa. Our study reveals phase changes that are explained from the calculated electronic structure. The cubic face-centered and body-centered structures (fcc and bcc), together with two-, three-, and four-layered hexagonal structures, play major roles under compression. The results’ dependence on the electron exchange and correlation in the DFT approach is investigated, and it is shown that the impact of the choice, while significant at lower pressures, diminishes in the terapascal regime. We further illustrate that the normal parabolic trends in atomic volume and bulk modulus with atomic number, due to the occupation of bonding and anti-bonding 5d states, break down at TPa pressures, suggesting drastically different chemical bonding at these extreme conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND

Machine learning approach for vibronically renormalized electronic band structures

Here, we present a machine learning (ML) method for efficient computation of vibrational thermal expectation values of physical properties from first principles. Our approach is based on the nonperturbative frozen phonon formulation in which stochastic Monte Carlo algorithm is employed to sample configurations of nuclei in a supercell at finite temperatures based on a first-principles phonon model. A deep-learning neural network is trained to accurately predict physical properties associated with sampled phonon configurations, thus bypassing the time-consuming ab initio calculations. To incorporate the point-group symmetry of the electronic system into the ML model, group-theoretical methods are used to develop a symmetry-invariant descriptor for phonon configurations in the supercell. We apply our ML approach to compute the temperature dependent electronic energy gap of silicon based on density functional theory (DFT). We show that, with less than a hundred DFT calculations for training the neural network model, an order of magnitude larger number of sampling can be achieved for the computation of the vibrational thermal expectation values. Our work highlights the promising potential of ML techniques for finite temperature first-principles electronic structure methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND

Engineering 2D Square Lattice Hubbard Models in 90° Twisted GeX/SnX (X =S, Se) Moiré Superlattices

Because of the large-period superlattices emerging in moiré two-dimensional (2D) materials, electronic states in such systems exhibit low energy flat bands that can be used to simulate strongly correlated physics in a highly tunable setup. While many investigations have thus far focused on moiré flat bands and emergent correlated electron physics in triangular, honeycomb, and quasi-one-dimensional lattices, tunable moiré realizations of square lattices subject to strong correlations remain elusive. Here, in this work, we propose a feasible scheme to construct moiré square lattice systems by twisting two or more layers of 2D materials in a rectangular lattice by 90°. We demonstrate the concept with twisted GeX/SnX (X =S, Se) moiré superlattices and calculate their electronic structures from first principles. We show that the lowest conduction flat band in these systems can be described by a square lattice Hubbard model with parameters which can be controlled by varying the choice of host materials, number of layers, and external electric fields. In particular, twisted double bilayer GeSe realizes a square lattice Hubbard model with strong frustration due to the next-nearest-neighbor hopping that could host unconventional superconductivity, in close analogy to the Hubbard model for copper-oxygen planes of cuprate high-temperature superconductors. The presented scheme uses 90° twisted 2D materials with rectangular unit cells as a promising platform for realizing the physical phenomena of square lattice Hubbard models, establishing a new route for studying its rich phase diagram of magnetism, charge order, and unconventional superconductivity in a highly tunable setting.

2-dimensional systems

Integration of Online Cross-Section Generation Capability with Depletion and Transient Solvers in Griffin

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE)-based reactor multiphysics analysis application jointly developed by Argonne and Idaho National Laboratories under the DOENE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In FY25, an online crosssection generation capability based on the Self-Shielding Application Programming Interface (SSAPI) was demonstrated for TRISO-fueled reactor problems under steady-state conditions. This fiscal year, that capability was extended to support depletion and transient multiphysics calculations, enabling high-fidelity analyses that generate self-shielded cross sections on the fly from the actual evolving composition and temperature states rather than from pre-tabulated libraries. For depletion, a two-way coupling was established in which SSAPI computes compact-averaged self-shielded cross sections that the depletion solver then uses to advance the Bateman equations, with the updated compositions returned to SSAPI at each step; the depletion module was refactored to support both library-based and SSAPI-based cross sections, and additional logic was added to track daughter isotopes and to exclude minor isotopes for efficiency. For transient analysis, the SSAPI multigroup library was extended with the kinetics data required for time-dependent calculations, the Improved Quasi-Static (IQS) scheme was coupled with SSAPI, and several supporting capabilities were implemented, including a self-shielding treatment that lets control rods and drums move within a self-shielded model, which had previously been impossible and had ruled out rod- and drum-movement transients with on-the-fly cross sections altogether, a new mixing scheme for delayed-neutron precursor decay constants, a checkpoint-based restart workflow, and performance improvements such as pointwise cross-section interpolation and the bypassing of unnecessary Dancoff factor calculations. The implemented capabilities were verified against Serpent Monte Carlo solutions. For depletion, a prismatic pin-cell problem based on a Next Generation Nuclear Plant (NGNP) Very High Temperature Reactor benchmark showed excellent agreement, with eigenvalue differences within 200 pcm over the entire burnup range (up to 140 MWD/kgU) and fission-product and actinide inventories agreeing to within 0.8% and 2.5%, respectively; a heat-pipe microreactor assembly problem with a much higher fuel loading confirmed the same behavior and quantified the bias introduced when the multigroup equivalence effect is neglected. For transient analysis, a pin-cell problem with a step reactivity insertion and temperature feedback reproduced the analytically expected asymptotic power and showed close agreement between the direct and IQS solutions, and a two-dimensional microreactor core problem with control-drum rotation exercised the new moving-drum self-shielding treatment and demonstrated successful coupling of the online crosssection generation with both the direct and IQS transient methods. The capability was further exercised on a full-core pebble-bed problem, in which Griffin was coupled with the System Analysis Module (SAM) to simulate load-following operation of the gPBR with the Doppler feedback resolved at the TRISO fuel kernel temperature. These developments in Griffin provide a convenient, high-fidelity approach to cross-section generation for advanced thermal reactors with geometrically complex and highly heterogeneous configurations, including TRISO-fueled prismatic and pebble-bed systems, and support steady-state, depletion, and transient multiphysics calculations. They also enable self-shielded cross sections to be evaluated directly at the actual coupled state of the system, thereby establishing a foundation for high-fidelity, fully coupled multiphysics analysis of advanced reactors

Park, H.

Identifying Suitable Front Contacts for High‐Efficiency Cd(Se,Te) Solar Cells on Space‐Qualified Cover Glass

Deployment of photovoltaics in space requires devices that combine high-efficiency, low areal mass, and resilience to harsh environments. Historically, high-efficiency multijunction III–V materials have dominated space power systems; however, their high cost and limited manufacturing throughput motivate the exploration of scalable alternatives. While CdTe-based thin-film photovoltaics offer an attractive option, their performance on non-conventional substrates can suffer from front contact instability under higher-temperature processing. Here, the role of front contact chemistry in limiting cell performance is investigated using CdTe-based devices fabricated on 150 μm thick Ceria-doped space-qualified 0214 Corning glass. A matrix of four transparent conducting oxides (TCOs: CTO, AZO, ITO, IZO) combined with two n-type emitters (MZO, IGO) reveals chemical stability at the front interface—rather than absorber composition alone—governs recombination losses, voltage deficits, and device reproducibility. Chemically stable front contact combinations suppress elemental diffusion and interfacial degradation, resulting in significantly improved carrier lifetimes and junction quality. These insights are validated through record-certified Cd(Se,Te) cell efficiencies of 18.4% under AM1.5G and 16.2% under AM0 illumination on ultra-thin glass. Beyond CdTe, this work provides a general framework for the rational selection of TCO/emitter interfaces in superstrate thin-film photovoltaics, including emerging technologies like metal halide perovskites, while enabling high-efficiency, lightweight photovoltaics for space applications.

14 SOLAR ENERGY

High Pressure X-ray Diffraction and Equation of State of Hydrazine

Synchrotron X-ray diffraction has been used to investigate the structure and equation of state (EOS) of hydrazine (N 2 H 4 ) up to 54.3 GPa at 298 K. The diffraction patterns could be fit to a monoclinic unit-cell structure and put strong constraints on previously reported phase transitions documented by vibrational spectroscopy over this pressure range. Pressure–volume ( P–V ) data were fit using a Vinet EOS, yielding parameters: V 0 = 45.2 Å 3 /molecule (fixed), K 0 = 11.8(7) GPa, and K 0 ′ = 6.5(2). Previously measured high-pressure vibrational frequency shifts were used to estimate the vibrational free energy and model P–V–T isotherms from 0 to 1200 K. The results of the P–V–T isotherms are compared to existing shock Hugoniot data on hydrazine and 298 K isotherms for assemblages of possible decomposition products. This comparison suggests dissociation at high density under shock loading. Good correspondence was found between the static lattice EOS as calculated by the model and the previously reported EOS as calculated by density functional theory. Finally, these results resolve existing uncertainties about the EOS and crystal symmetry of hydrazine at high pressure and provide valuable baseline information on this important energetic material.

diffraction

Performance of Heterostructural TaC/AlGaN Schottky Diodes Based on First Principles Electronic Structure Properties

Advances in ultra-wide bandgap materials, such as high Al-content AlxGa1-xN (AlGaN), are essential for next generation power electronics, but the requirement for lattice matched substrates is currently a significant obstacle. Recently, conductive TaC has emerged as a promising virtual substrate for AlGaN heteroepitaxy, with wurtzite (0001) AlxGa1-xN lattice-matched to rocksalt (111) TaC at x ~ 0.5. Thus, understanding and controlling the electronic properties of the TaC/AlGaN interface is key for developing technological applications based on TaC/AlGaN devices. Using density functional theory and electronic structure calculations, we here investigate TaC/Al0.5Ga0.5N interfaces, where we include explicit alloy models in the slab calculations. We predict the Schottky barrier height and the electric field discontinuity resulting from interface charges. Considering all possible combinations of (Ta, C) substrate termination, (Al/Ga, N) nucleation, and (Al/Ga, N) polarity, we construct a chemical potential phase diagram to identify the stable interfaces that can be accessed through variation of the synthesis conditions. The predicted interface electronic properties are implemented in device performance simulations to demonstrate a practical design for a strain-free, high-efficiency TaC/AlGaN Schottky diode with a low barrier height and without interface charges, underscoring the potential of TaC as a substrate for ultra-wide bandgap devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC

Mu2e Production Target Thermal Test Project

The Mu2e experiment requires a production target that is capable of operating under extreme thermal conditions caused by an 8 GeV proton beam. This project’s objective supports the development of the Mu2e Pro- duction Target by testing the Stickman model’s thermal behavior. Angel Flores Luviano has assisted Jonathan Williams in progressing this project by contributing to the development of a Radiative Cooling Test Fixture (RCTF) that will be used to evaluate the thermal behavior of the new production target model. Engineering calculations were performed to an- alyze thermal performance and pressure drop within the cold well and wa- ter cooling circuit. These calculations also determined the optimal sizing for key components of the water system. An engineering note was made to document these calculations. CAD models of the water circuit piping, thermocouple mounting bars, and radiator chimney were developed, and prototype test rig components were fabricated using 3D printing. Future work will focus on continued RCTF development which includes control system integration, heater hardware design, and interfaces that can be scaled up to increase thermal capacity.

Flores Luviano, Angel Augusto [Unlisted, US]

Mechanistic origin of solvent-dependent thermal stability in sodiated Sn anodes for sodium-ion batteries

Understanding the thermal stability of high-energy density alloy anodes is critical for the safe deployment of sodium-ion batteries (SIBs). Here, accelerating rate calorimetry (ARC), post-mortem characterizations, and density functional theory (DFT) calculations are combined to understand the thermal reactivity of fully sodiated Sn, Sn-hard carbon (HC) blends, and HC anodes in carbonate- and ether-based electrolytes. ARC measurements show that propylene carbonate (PC) causes earlier self-heating rate (SHR) onset and higher reactivity than tetraethylene glycol dimethyl ether (TEGDME), indicating inferior thermal stability. Sodiated Sn exhibits better thermal stability than sodiated HC, while Sn-HC blends show intermediate behavior that improves with increasing Sn content. Post-ARC analyses reveal desodiation of Na15Sn4 to metallic Sn with particle coalescence, whereas Sn-HC blends and HC retain finer morphologies. PC promotes Sn oxidation to SnO, while TEGDME suppresses oxide formation; NaPF6-containing electrolytes additionally form NaF. DFT calculations show that PC adsorption lowers Na extraction energy and enhances interfacial electronic interactions, facilitating Na release and reductive decomposition. These results establish a direct correlation between solvent-dependent reaction pathways and thermal stability in SIB alloy anodes.

Accelerating rate calorimetry

Collective excitations and low-energy ionization signatures of relativistic particles in silicon detectors

Abstract Solid-state detectors with a low energy threshold have several applications, including searches of non-relativistic halo dark-matter particles with sub-GeV masses. When searching for relativistic, beyond-the-Standard-Model particles with enhanced cross sections for small energy transfers, a small detector with a low energy threshold may have better sensitivity than a larger detector with a higher energy threshold. In this paper, we calculate the low-energy ionization spectrum from high-velocity particles scattering in a dielectric material. We consider the full material response including the excitation of bulk plasmons. We generalize the energy-loss function to relativistic kinematics, and benchmark existing tools used for halo dark-matter scattering against electron energy-loss spectroscopy data. Compared to calculations commonly used in the literature, such as the Photo-Absorption-Ionization model or the free-electron model, including collective effects shifts the recoil ionization spectrum towards higher energies, typically peaking around 4–6 electron-hole pairs. We apply our results to the three benchmark examples: millicharged particles produced in a beam, neutrinos with a magnetic dipole moment produced in a reactor, and upscattered dark-matter particles. Our results show that the proper inclusion of collective effects typically enhances a detector’s sensitivity to these particles, since detector backgrounds, such as dark counts, peak at lower energies.

Physics

Electronegativity-Guided Site Differentiation in High-Entropy Alloy for pH-universal Hydrogen Evolution Reactions

Enhancing the intrinsic activity of transition metal catalysts for the hydrogen evolution reaction (HER) remains a critical challenge in sustainable energy conversion. Herein, we report an electronegativity-guided site differentiation strategy in a single-phase CoNiCuMoW high-entropy alloy (HEA) via electrodeposition by incorporating high-electronegativity 4d/5d orbital transition metals (Mo, W) into the face-centered cubic (fcc) matrix (CoNiCu). The as-synthesized HEA demonstrates exceptional HER performance in all pH conditions, delivering an outstanding overpotential of 65 mV (alkaline), 28 mV (acidic), and 155 mV (neutral) at a current density of 100 mA cm−2, showing performance comparable to commercial Pt/C and has excellent long-term stability at high current density (1 A cm−2, 1000 h). X-Ray absorption spectroscopy (XAS) and density functional theory (DFT) calculations reveal that the incorporation of Mo/W simultaneously alters the local coordination environment and induces element-dependent charge redistribution, accompanied by a system-level d-band center downshift, thereby optimizing the hydrogen binding strength across multimetallic sites. Meanwhile, oxophilic Mo/W sites lower the water dissociation energy barrier. These synergistic effects collectively enable efficient and durable pH-universal HER performance.

Wu, Yutong

Identifying Strain Stacking Boundaries between Multiphase Domains in Atomically Thin Two-Dimensional Magnets

Stacking engineering of van der Waals materials is an important strategy to control the materials’ properties, such as electronic correlations, ferroelectricity, and layer-dependent two-dimensional magnetism. A timely testbed for the study of the latter is atomically thin chromium trihalides (CrX 3 , X = Cl, Br, I). Notably, by understanding the sliding mechanism between different stacking sequences, control of the stacking arrangement, and thus magnetic properties in CrX 3 , can be achieved. Such insight, however, is currently lacking. Here, in this study, advanced electron microscopy methods are used to identify multiple stacking sequences corresponding to different bulk phases in atomically thin CrX 3 (X = Cl and Br) down to bilayer thickness and with lateral domain sizes as small as tens of nanometers. Indications of nanometer scale transitions and interactions at the stacking boundaries are found, including a universally preferred sliding direction that is consistent with density functional theory calculations and the strain fields at lateral heterostructure boundaries. This study demonstrates the necessity to consider local stacking structures when interpreting averaged magnetic properties measured with macroscale probes. Additionally, the preferred sliding direction insight from this study provides a strategy to control stacking sequence in atomically thin CrX 3 samples during the exfoliation and sample fabrication process.

DFT calculations

Quasicrystal stability and nucleation kinetics from density functional theory

The aperiodic order of quasicrystals bridges the amorphous and crystalline regime, so it has remained unclear whether quasicrystals are metastable or stable phases of matter. Density functional theory is often used to evaluate thermodynamic stability, but quasicrystals are long-range aperiodic and their energies cannot be calculated using conventional ab initio methods. Here, in this work, we perform first-principles calculations on quasicrystal nanoparticles of increasing size, from which we can directly extrapolate their bulk and surface energies. Using this technique, we determine with high confidence that the icosahedral quasicrystals ScZn 7.33 and YbCd 5.7 are ground-state phases, thus revealing that translational symmetry is not a necessary condition for the zero-temperature stability of inorganic solids. Although we found the ScZn 7.33 quasicrystal to be thermodynamically stable, we show on a mixed thermodynamic and kinetic phase diagram that its solidification from the melt is limited by nucleation, which illustrates why even stable materials may be kinetically challenging to grow. Our techniques broadly open the door to first-principles investigations into the structure–bonding–stability relationships of aperiodic materials.

density functional theory

Three-dimensional dispersion in the type-II Dirac semimetals PtTe 2 and PdTe 2 revealed through circular dichroism in angle-resolved photoemission spectroscopy

PtTe 2 and PdTe 2 are among the first transition metal dichalcogenides that were predicted to host type-II Dirac fermions, exotic particles prohibited in free space. These materials are layered and air stable, which makes them top candidates for technological applications that take advantage of their anisotropic magnetotransport properties. Here, in this work, we provide a detailed characterization of the electronic structure of PtTe 2 and PdTe 2 using angle-resolved photoemission spectroscopy (ARPES) and density functional theory calculations, offering an alternative interpretation for one of the Dirac-like dispersions in these materials. Through the use of circularly polarized light, we report a different behavior of such dispersion in PdTe 2 compared to PtTe 2 , that we relate to a symmetry analysis of the dipole matrix element. Such analysis reveals a link between the observed circular dichroism and the different momentum-dependent terms in the dispersion of these two compounds, despite their close similarity in crystal structure. Additionally, our data show a clear difference in the circular dichroic signal for the type-II Dirac cones characteristic of these materials, compared to their topologically protected surface states. Our paper provides a useful reference for the ARPES characterization of other transition metal dichalcogenides with topological properties and illustrates the use of circular dichroism as a guide to identify the topological character and attributes of two otherwise equivalent band dispersions.

angle-resolved photoemission spectroscopy

Systematic determination of a material’s magnetic ground state from first principles

Abstract We present a self-consistent method based on first-principles calculations to determine the magnetic ground state of materials, regardless of their dimensionality. Our methodology is founded on satisfying the stability conditions derived from the linear spin wave theory (LSWT) by optimizing the magnetic structure iteratively. We demonstrate the effectiveness of our method by successfully predicting the experimental magnetic structures of NiO, FePS 3 , FeP, MnF 2 , FeCl 2 , and CuO. In each case, we compared our results with available experimental data and existing theoretical calculations reported in the literature. Finally, we discuss the validity of the method and the possible extensions.

Chemistry

Stabilization Of The CN 3 5− Anion In Recoverable High‐pressure Ln 3 O 2 (CN 3 ) (Ln=La, Eu, Gd, Tb, Ho, Yb) Oxoguanidinates

Abstract A series of isostructural Ln 3 O 2 (CN 3 ) (Ln=La, Eu, Gd, Tb, Ho, Yb) oxoguanidinates was synthesized under high‐pressure (25–54 GPa) high‐temperature (2000–3000 K) conditions in laser‐heated diamond anvil cells. The crystal structure of this novel class of compounds was determined via synchrotron single‐crystal X‐ray diffraction (SCXRD) as well as corroborated by X‐ray absorption near edge structure (XANES) measurements and density functional theory (DFT) calculations. The Ln 3 O 2 (CN 3 ) solids are composed of the hitherto unknown CN 3 5− guanidinate anion—deprotonated guanidine. Changes in unit cell volumes and compressibility of Ln 3 O 2 (CN 3 ) (Ln=La, Eu, Gd, Tb, Ho, Yb) compounds are found to be dictated by the lanthanide contraction phenomenon. Decompression experiments show that Ln 3 O 2 (CN 3 ) compounds are recoverable to ambient conditions. The stabilization of the CN 3 5− guanidinate anion at ambient conditions provides new opportunities in inorganic and organic synthetic chemistry.

Chemistry

Modeling a generic TRISO-fueled heat pipe microreactor using SCALE: Depletion, transportation criticality, and shielding

This paper demonstrates the applicability of the SCALE code system to tristructural-isotropic (TRISO)-fueled heat pipe microreactors through depletion, transportation criticality, and shielding analyses of a generic reference design. The study conducted supports US Nuclear Regulatory Commission code readiness efforts for advanced non–light-water reactor technologies and is intended as a code capability demonstration rather than as an optimization of a specific microreactor design. The modeled reactor employs high-assay low-enriched uranium (HALEU) uranium oxycarbide (UCO) TRISO fuel and beryllium oxide (BeO) reflectors and operates at 7.5 MWth with a nominal lifetime of about 3 effective full power years. Representative cases for fresh and irradiated cores were selected to exercise SCALE methods relevant to reactor operation and post-irradiation transport. The discharged-core decay heat is approximately 6% of operating power immediately after shutdown. Transportation criticality calculations show that internal water ingress is the dominant reactivity effect, with fully flooded fresh core and irradiated core configurations remain above the subcriticality criterion, even with the available control mechanisms. Shielding calculations for a simplified transportation package indicate that normal-condition dose rates are governed mainly by shielding thickness and cooling time, whereas the breached hypothetical accident case is governed primarily by cooling time. Overall, the study shows that SCALE supports depletion, transportation criticality, and shielding evaluations efficiently for TRISO-fueled heat pipe microreactors within a single code system.

Criticality