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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

1,272 records · Page 3

All sources NASA NTRS DOE OSTI

NASA NTRS2026

Diffusion Quantum Monte Carlo Calculation of the Austenite and Martensite Phases of NiTi

NiTi is a promising material for smart and active technologies due to its exhibition of the shape memory effect, superelasticity, and biocompatibility. The shape memory effect is tied to the reversible transition between the austenite and martensite phases. A major research direction is to alloy NiTi with Zr, Hf, Pd, Pt, etc., in order to tune the martensitic transition temperature (MTT). Modeling the MTT from first principles is challenging because the lattice dynamics is complicated by anharmonicity and various low-energy structures. Using density functional theory, the energy difference between the austenite and martensite phases of NiTi varies by up to 100 meV/atom depending on the choice of density functional, which is of the same order of the energy difference itself. Consequently, free energy calculations with different functionals can result in estimates of the MTT that vary by several hundred K. Using diffusion quantum Monte Carlo, we calculated the energy difference between the B2 and B19' structures of NiTi to be 70.9 +- 2.5 meV/atom.

Kevin K Ly↗

DOE OSTI2026

Comparative Assessment of Battery Carbon Footprint Calculation Frameworks to Support U.S. Battery Manufacturing

A “battery passport” is a digital record that provides comprehensive information about an individual battery across its life cycle. This concept has been introduced by Battery Regulation (EU) 2023/1542 [1], which applies to batteries for electric vehicles (EVs), industrial batteries with a capacities greater than 2 kWh, and batteries for light means of transport (LMT) greater than 2kWh sold in the European Union (EU) market – impacting manufacturers and exporters across multiple jurisdictions globally. Among its reporting mandates, a key feature of battery passports is the requirement for a carbon footprint (CF) calculation methodology supported by enhanced data granularity to enable traceable and verifiable CF results. While no other jurisdictions have yet adopted formal battery passport requirements like the EU’s, some are developing CF calculation methods to comply with the EU Battery Regulation or are creating CF-related regulations that could evolve in a similar direction, reflecting the growing importance of battery CF guidelines for manufacturers seeking to remain competitive in global markets.

Zhang, Jingyi [Argonne National Laboratory (ANL), ↗

NASA NTRS1992

Lightweight bipolar storage battery

An apparatus [10] is disclosed for a lightweight bipolar battery of the end-plate cell stack design. Current flow through a bipolar cell stack [12] is collected by a pair of copper end-plates [16a,16b] and transferred edgewise out of the battery by a pair of lightweight, low resistance copper terminals [28a,28b]. The copper terminals parallel the surface of a corresponding copper end-plate [16a,16b] to maximize battery throughput. The bipolar cell stack [12], copper end-plates [16a,16b] and copper terminals [28a,28b] are rigidly sandwiched between a pair of nonconductive rigid end-plates [20] having a lightweight fiber honeycomb core which eliminates distortion of individual plates within the bipolar cell stack due to internal pressures. Insulating foam [30] is injected into the fiber honeycomb core to reduce heat transfer into and out of the bipolar cell stack and to maintain uniform cell performance. A sealed battery enclosure [ 22] exposes a pair of terminal ends [26a,26b] for connection with an external circuit.

Rowlette, John J.↗

DOE OSTI2024

Calculating a minimum overlap period for successful intercalibration of soil moisture sensors

Abstract Long‐term in situ soil moisture monitoring inevitably requires sensors to be replaced. Ensuing discontinuities in the data record can be mitigated by intercalibration, however it is unclear how long the existing sensor needs to remain alongside the newly installed before there is enough overlapping data to generate a robust intercalibration. We used 154 pairs of established and newly installed sensors within the Marena, Oklahoma, In Situ Sensor Testbed to determine if there is a minimum overlap time that should be considered when planning upcoming replacements. Hourly observations of the existing sensor were linearly calibrated to those of the newly installed sensor with coefficients determined from overlap periods incremented by 30 days until a reference period of 2 years was reached. The resulting bias, root‐mean‐square error, and correlation coefficient for sensor pairs indicate that a minimum of 6 to 9 months of overlapping data are required to generate a successful intercalibration. Extending that to a full year before decommissioning the old sensor results in a stable intercalibration with higher confidence.

Agriculture↗

DOE OSTI2025

CASM Monte Carlo: Calculations of the thermodynamic and kinetic properties of complex multicomponent crystals

Monte Carlo techniques play a central role in statistical mechanics approaches that connect macroscopic thermodynamic and kinetic properties to the electronic structure of a material. This paper describes the implementation of Monte Carlo techniques for the study of multicomponent crystalline materials within the Clusters Approach to Statistical Mechanics (CASM) software suite, and demonstrates their use in model systems to calculate free energies and kinetic coefficients, study phase transitions, and construct phase diagrams from first principles. Many crystal structures are complex, with multiple sublattices occupied by differing sets of chemical species, along with the presence of vacancies or interstitial species. This imposes constraints on concentration variables, the form of thermodynamic potentials, and the values of kinetic transport coefficients. The framework used by CASM to formulate thermodynamic potentials and kinetic transport coefficients accounting for arbitrarily complex crystal structures is presented and demonstrated with examples of increasing complexity. Additionally, an overview of the capabilities of the CASM software specific to Monte Carlo methods is given, and a new CASM software package is introduced, casm-flow, which helps automate the setup, submission, management, and analysis of Monte Carlo simulations.

Cluster expansion↗

DOE OSTI2026

The CP-PAW Code Package for First-Principles Calculations from a User’s Perspective

CP-PAW is a combined electronic structure and ab initio molecular dynamics code to perform mixed quantum and classical simulations of atomistic condensed phase systems, such as solids, liquids, and molecular systems. As the name suggests, the CP-PAW code unifies the all-electron projector augmented-wave (PAW) method with the Car–Parrinello (CP) approach to determine not only the electronic and nuclear ground states of condensed matter but also to study their properties and dynamics. In addition to briefly outlining the underlying theory, the focus will be on the unique aspects of CP-PAW and how to correctly employ them as a user. How to install CP-PAW using the new build system will also be briefly mentioned.

Blöchl, Peter E [Institute for Theoretical Physic↗

DOE OSTI2026

Reducing experimental uncertainty in the calculation of cavity swelling in a transmission electron microscope through crystallographically aided void volume tracking

Accurate measurement of cavity swelling in a transmission electron microscope is essential to define material performance under irradiation, and the conventionally applied spherical assumption for the calculation of cavity volumes in irradiated materials can result in errors ranging between a 25% underestimation and 72% overestimation of volume purely based on the assumption of shape. This assumption is undeniably expedient for calculation but does not sufficiently account for the 3D nature of polyhedral cavities and their shape projection in the transmission electron microscope, and therefore presents too large of an associated uncertainty in swelling determination for faceted cavities. This uncertainty has been defined for many common cavity shapes in FCC and BCC crystal systems, and has been tabulated across the cubic fundamental region. A revised methodology for crystallographically aided void volume tracking, or CAVV-T, is presented and demonstrated on a specimen of neutron irradiated Ni. In-depth discussion on the application of this technique is provided along with resources to allow for conversion between the spherical assumption and this revised method. This work seeks to increase experimental confidence in the characterization and quantification of critical aspects of irradiation damage in materials by applying a crystallographically-resolved approach for cavity swelling calculation.

Haag IV, James V. [Pacific Northwest National Labo↗

DOE OSTI2026

Electronic structure of the kagome compound CaTi 3⁢ Bi 4 using high-field torque magnetometry and density functional theory

Here, we report systematic torque magnetometry measurements to investigate the electronic properties of the newly discovered kagome compound CaTi 3 ⁢Bi 4 . Electrical transport, magnetic susceptibility, and thermal measurements reveal no evidence of a magnetic ground state in this material. Torque data obtained in magnetic fields up to 41.5 T exhibit clear de Haas–van Alphen (dHvA) oscillations, with nine distinct frequencies ranging from 13 to 6164 T. Angular-dependent dHvA measurements show that, with the exception of the lowest frequency (13 T), all observed frequencies nearly follow a 1/cos ⁡𝜃 dependence, where 𝜃 is the angle between the crystallographic 𝑐 axis and the magnetic field direction. This behavior is characteristic of quasi-two-dimensional Fermi-surface sheets with nearly circular cross sections. To further elucidate the electronic structure, we performed density functional theory (DFT) calculations of the band structure and Fermi surface. The calculated bands reveal the presence of multiple Dirac points (DP), flat bands (FB), and van Hove singularities (VHS) near the Fermi level. The resulting Fermi surface consists of several quasi-two-dimensional cylindrical sheets, consistent with the experimentally observed 1/cos⁡ 𝜃 dependence. Notably, the theoretical dHvA frequencies, derived from extremal Fermi-surface cross-sectional areas, agree well with the experimental values and reproduce their angular dependence. Remarkably, all experimentally observed frequencies are captured by the DFT predictions. Pressure-dependent calculations up to 10 GPa show that the electronic features—DP, FB, and VHS—evolve systematically with pressure. In particular, the VHS shifts closer to the Fermi level, demonstrating that pressure acts as an effective tuning parameter in this material. These combined experimental and theoretical results provide a comprehensive understanding of the electronic structure of CaTi 3 ⁢Bi 4 and demonstrate how pressure can be used to tune its key electronic features, offering valuable guidance for exploring related kagome materials.

Shtefiienko, Kyryl [West Texas A & M University, C↗

DOE OSTI2026

Extreme confinement of hydrogen gas within fullerenelike nanoporous carbon

Nanoporous carbons and carbon nanostructures can store hydrogen at cryogenic temperatures but lack the volumetric and gravimetric capacity to be industrially significant. Recent inelastic neutron scattering experiments suggest a highly dense phase of hydrogen at temperatures well above the melting point of solid hydrogen. However, it remains unclear how pore geometry and intermolecular interactions enable these dense phases to exist, with dispersion (van der Waals) or electrostatic/induction suggested to be the key effects in slit and curved pores but their relative contributions have yet to be quantified. In this paper, we perform benchmark electronic structure calculations allowing the interactions between planar and curved aromatic molecules with hydrogen to be accurately determined. Dispersion was found to dominate over electrostatic and inductive effects with some many-body charge transfer (Dobson type-A) effects needed to capture the most highly curved structures. Density functional methods that include type-A many-body effects were found to accurately describe the intermolecular interactions at a fraction of the cost of coupled-cluster simulations and these approaches were used to calculate the energies inside large carbon bowl and slit pores. The interaction energies inside the bowl pores were found to depend on the orientation of the hydrogen molecule. This rotational barrier, modeled as a quantum hindered rotor, could reproduce the peak splitting observed in inelastic neutron scattering experiments, with weak splitting arising from bowl-like fullerene pores and strong splitting from highly confining nanotubelike pores. Increasing the fraction of such curved pores in nanoporous carbons may therefore offer a pathway to enhance their hydrogen-storage capacity. Moreover, the preferential adsorption of ortho hydrogen on nanotubelike pores could enable the storage of high-density hydrogen without the need to remove heat produced during the ortho-para hydrogen conversion.

36 MATERIALS SCIENCE↗

NASA NTRS1977

Calculations of the Moon's Thermal History at Different Concentrations of Radioactive Elements, Taking Into Account Differentiation on Melting

Calculations of the thermal history of the moon were carried out by solving the thermal conductivity equation for the case in which the heat sources are the long-lived radioactive elements Th, U, and K 40 . The concentrations of these elements were adjusted to give four variations of the heat flow: 1.35 x 10 -8 and 0.91 x 10 -6 cal cm -2 s -1 (I and I**), 0.61 x 10 -6 (variant II, the terrestrial mixture of Lyubimova and Starkova), and 0.236 x 10 -6 (variant III, the chondrite model of Urey and MacDonald). In the same calculations, we considered layering of the differentiated material with transport to the surface of the radioactive elements after the temperature of the layer rose to 200 K above the melting temperature, which is considered in five variants that differ in the amount of transported radioactive elements: 100 percent ( n = 1), 80 percent n = 0.8), 40 percent, 20 percent ( n = 0.6, 0.4, 0.2). During fusion the heat capacity and heat conductivity were changed. We considered two variants of an initially cold (273 K) and hot (900 K) Moon. Calculations show that the interior of the Moon was heated to melting during the first 0.7 to 2.3 x 10 9 years. The maximum fusion involved practically the entire Moon to a distance from 15 to 45 km beneath the surface, and started 3.5 to 4.0 x 10 9 years ago (I, I**), or 2.5 to 3.0 x 10 9 years ago (II, III) and continued for 1 to 2 x 10 9 years. Today the Moon is cooling. The current thickness of the solid crust is from 150 to 200 km and the heat flow exceeds the stationary value 1.5 fold. Apparently the most realistic variant is II (terrestrial mixture) for an initially hot Moon, and gives, regardless of the moderate concentration of radioactive elements, a heat flow of 0.9 to 0.95 x 10 -6 cal cm -2 s -1 , which agrees with the radioastronomical measurements of Troitsky and Krotikov and with the Apollo 15 data.

O I Ornatskaya↗

NASA NTRSDate not supplied

Laser Beam Welding Benchmark Experiments Performed in Reduced Gravity and Vacuum

Laser beam welding (LBW) is affected by the extreme temperatures, reduced pressure, and reduced gravity present in space environments. Gravity and pressure especially influence its melt pool and solidification dynamics. A compact, modular vacuum chamber adaptable to flight platforms from parabolic to orbital currently hosts an experiment to investigate the combined influence of reduced gravity and pressure on LBW. A swappable cartridge contains a rotating platen on which customizable workpieces can be welded under vacuum, greatly increasing experimental throughput. Instrumentation includes weld and thermal cameras observing the process, thermocouples placed on workpieces, accelerometers, and vacuum sensors. Experimental data gathered during the welding process will be combined with post-flight nondestructive evaluation, metallography, and mechanical testing to provide validation datasets for computational modeling. Phase I of this effort involves a parabolic flight campaign in low gravity while an anticipated Phase II would proceed to in-space demonstration to access extended duration microgravity.

in-space welding↗

NASA NTRS2026

Assessment of Near-Angle Scatter on Exo-Earth Coronagraphy

Near-Angle Scatter (NAS) of the host star’s light into the coronagraph dark hole may limit the ability of a potential Habitable Worlds Observatory (HWO) to detect and characterize an Earth-like planet around a Sun-like star. This paper summarizes a 5-year investigation of the impact of NAS on high-contrast coronagraphy. Science requirements were flowed down to a draft flux ratio noise ratio (FRN) error budget with allocations for NAS. The modeled sources of NAS include correlated (low and mid-spatial) surface structure, uncorrelated surface microroughness and coating structure, particulate contamination, micrometeoroid impacts and polarization leakage. The paper develops specifications for these sources of scattered light that meet their FRN error budget allocations, distinguishing between scatter that can and cannot be modulated. The development process utilizes an analytical expression that predicts scatter throughput into the dark hole based on BRDF/BSDF. Using this relationship, Photon Engineering performed a FRED® straylight analysis of the HWO EAC-1 (Exploratory Analytical Case #1) coronagraph to parametrically study predicted scattered light contrast in the dark hole. The model revealed that primary mirror microroughness or contamination is not a significant source of scattered light; the most critical surfaces are the smallest and those located closest to the coronagraph focal plane mask. The assessment found that uncorrelated microstructure, contamination, and micrometeoroid impacts require more investigation and potential significant sources of unmodulated scattered light, and highlights that slow thermal wavefront drift is an order of magnitude more impactful than fast mechanical jitter.

coronagraphy↗

DOE OSTI2025

Carbon in GaN as a nonradiative recombination center

Trap-assisted nonradiative recombination has been shown to limit the efficiency of optoelectronic devices. While substitutional carbon (⁠C N ) has been suggested to be a nonradiative recombination center in GaN devices, a complete recombination cycle including the two charge-state transition levels has not been previously described. In this work, we investigate the trap-assisted recombination process due to C N in GaN, including multiphonon emission, radiative recombination, trap-assisted Auger–Meitner (TAAM) recombination as well as thermal emission of holes. Our study shows the key role of TAAM processes at the high carrier densities relevant for devices. We also reveal the carrier-density regimes where thermal emission and radiative recombination are expected to play an observable role. Our results highlight that carbon concentrations exceeding ~10 17 cm −3 can have a noticeable impact on device efficiency, not just in GaN active layers but also in InGaN and AlGaN. Furthermore, our comprehensive formalism not only offers detailed results for carbon but also provides a general framework for assessing the multiple processes that participate in trap-assisted recombination in semiconductors.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Distinguishing Thermal Fluctuations from Polaron Formation in Halide Perovskites

Recent angle-resolved photoelectron spectroscopy (ARPES) measurements of the hole effective mass in CsPbBr$_3$ revealed an enhancement of $\sim$50 % compared to the bare mass computed from first principles for CsPbBr$_3$ at $T = 0 K$. This large enhancement was interpreted as evidence of polaron formation. Employing accurate finite-temperature first-principles calculations, we show that the calculated hole effective mass of CsPbBr$_3$ at $T = 300 K$ can explain experimental results without invoking polarons. Thermal fluctuations are particularly strong in halide perovskites compared to conventional semiconductors such as Si and GaAs, and cannot be ignored when comparing with experiment. We not only resolve the debate on polaron formation in halide perovskites, but also demonstrate the general importance of including thermal fluctuations in first-principles calculations for strongly anharmonic materials.

36 MATERIALS SCIENCE↗

DOE OSTI2026

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↗

DOE OSTI2026

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↗

DOE OSTI2026

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↗

NASA NTRS1977

Riddles About the Origin and Thermal History of the Moon

Analyses of lunar rocks have confirmed that the interior of the Moon has been subjected to magmatic differentiation, as has been determined by calculation of the thermal history of the Moon. However, it appears that differentiation occurred so early that it could not have been related to heating by long-lived radioactive elements, but rather indicates a high initial temperature of the Moon. The reasons for the high initial temperature remain unknown. The presence of a central, partially melted area revealed by seismic observations forces us to reject convective models of the Moon, while the dimensions of this area, together with the measured heat fluxes from the interior, serve as "boundary conditions" for any new calculations of the thermal history for a two-layered differentiated model of the Moon operative over the past 3.5 billion years. If we start with the early differentiation of the Moon and the high content of long-lived radioactive elements in it, explanation of its current thermal properties represents no difficulty. But these prerequisites themselves remain unexplained. Due to the continuing controversy concerning the true nature of the process of magmatic differentiation of the Moon, estimation of the composition of the entire Moon on the basis of analyses of surface rock remains unreliable, and use of these estimates to clarify the origin of the Moon is premature. Existing hypotheses concerning the origin of the Moon encounter dynamic difficulties, as well as difficulties with explanation of the high initial temperature, to say nothing of the difficulty involved in explaining the proposed differences in the composition of the Earth and Moon.

B Yu Levin↗