Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computer Codes”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Advanced capabilities for materials modelling with Quantum ESPRESSO

Quantum ESPRESSO is an integrated suite of open-source computer codes for quantum simulations of materials using state-of-the-art electronic-structure techniques, based on density-functional theory, density-functional perturbation theory, and many-body perturbation theory, within the plane-wave pseudopotential and projector-augmented-wave approaches. Quantum ESPRESSO owes its popularity to the wide variety of properties and processes it allows to simulate, to its performance on an increasingly broad array of hardware architectures, and to a community of researchers that rely on its capabilities as a core open-source development platform to implement their ideas. In this paper we describe recent extensions and improvements, covering new methodologies and property calculators, improved parallelization, code modularization, and extended interoperability both within the distribution and with external software.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Correlation matrix renormalization theory in multi-band lattice systems

An appropriate treatment of electronic correlation effects plays an important role in accurate descriptions of physical and chemical properties of real materials. The recently proposed correlation matrix renormalization theory with sum rule correction (CMR) for studying correlated-electron materials has shown good performance in molecular systems and a periodic hydrogen chain in comparison with various quantum chemistry and quantum Monte Carlo calculations. Additionally, this work gives a detailed formulation and computational code implementation of CMR in multi-band periodic lattice systems. This lattice CMR ab initio theory is highly efficient, has no material specific adjustable parameters, and has no double counting issues faced by the hybrid approaches like LDA + U, DFT + DMFT and DFT + GA type theories. Benchmark studies on materials with s and p orbitals in this study show that CMR in its current implementation consistently performs well for these systems as the electron correlation increases from the bonding region to the bond breaking region.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Calculation of RF sheath properties from surface wave-fields: a post-processing method

In ion cyclotron range of frequency (ICRF) experiments in fusion research devices, radio frequency (RF) sheaths form where plasma, strong RF wave fields and material surfaces coexist. These RF sheaths affect plasma material interactions such as sputtering and localized power deposition, as well as the global RF wave fields themselves. RF sheaths may be modeled by employing a sheath boundary condition (BC) in place of the more customary conducting wall BC; however, there are still many ICRF computer codes that do not implement the sheath BC. In this paper we present a method for post-processing results obtained with the conducting wall BC. The post-processing method produces results that are equivalent to those that would have been obtained with the RF sheath BC, under certain assumptions. We found the post-processing method is also useful for verification of sheath BC implementations and as a guide to interpretation and understanding of the role of RF sheaths and their interactions with the waves that drive them.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

The Thermonuclear Runaway and the Classical Nova Outburst

Nova explosions occur on the white dwarf component of a cataclysmic variable binary stellar system that is accreting matter lost by its companion. When sufficient material has been accreted by the white dwarf, a thermonuclear runaway occurs and ejects material in what is observed as a classical nova explosion. Here, we describe both the recent advances in our understanding of the progress of the outburst and outline some of the puzzles that are still outstanding. We report on the effects of improving both the nuclear reaction rate library and including a modern nuclear reaction network in our one-dimensional, fully implicit, hydrodynamic computer code. In addition, there has been progress in observational studies of supernovae Ia with implications about the progenitors, and we discuss that in this review.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Practical overview of the ICRP Publication 145 phantoms for use with MCNP

Abstract This paper provides a practical overview of the mesh-type reference phantoms (MRCPs) published by the ICRP by providing supplementary guidance on the proper implementation of using the MRCPs in computational codes, such as Monte Carlo N-Particle. Various validation scenarios for photon and neutron exposures were also examined where published values were reproduced using the MRCPs for both external and internal exposures. Generally, the MRCPs performed well in the validation scenarios with the values reproduced being within 10% of the published values.

Hetrick, Lucas (ORCID:0000000178978007)↗
DOE OSTI2025

Competing ionization and dissociation: Extension of the energy-dependent frame transformation to the gerade symmetry of H 2

This article solves two major tasks that frequently arise in the theory of electron collisions with a target molecular cation. First, it extends the energy-dependent frame transformation (EDFT) treatment, which is needed to map fixed-nuclei electron-molecule scattering matrices into an energy-dependent laboratory-frame scattering matrix with vibrational channel indices. The EDFT mapping can now be carried out even when the target molecule possesses multiple low-energy potential curves, significantly transcending previous applications. Second, it implements a method to extract the rest of the full laboratory-frame scattering matrix, i.e., the columns and rows describing input and/or output dissociation channels. The treatment is benchmarked in this article against the essentially exact solution of a refined two-dimensional model of the singlet gerade Σ symmetry of H 2 . Our tests demonstrate that the theory accurately maps fixed-nuclei scattering information, of the type provided by existing electron-molecule computer codes, into a laboratory-frame scattering matrix that includes both ionization and dissociation. Furthermore, this treatment can provide a general framework applicable to a broad class of electron collision processes involving diatomic target ions, suitable for an accurate description of challenging processes such as dissociative recombination.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Neutron-neutron scattering length from the 6 He(p,pα)nn reaction

We propose a novel method to measure the neutron-neutron scattering length using the 6 He(p,pα)nn reaction in inverse kinematics at high energies. The method is based on the final-state interaction (FSI) between the neutrons after the sudden knockout of the α particle. We show that the details of the neutron-neutron relative-energy distribution allow for a precise extraction of the s-wave scattering length. Here, we present the state of the art in regard to the theory of this distribution. The distribution is calculated in two steps. First, we calculate the ground-state wave function of 6 He as a αnn three-body system. For this purpose we use Halo effective field theory, which also provides uncertainty estimates for the results. We compare our results at this stage to model calculations done with the computer code face. In a second step we determine the effects of the nn FSI using the nn t-matrix. We compare these FSI results to approximate FSI approaches based on standard FSI enhancement factors. While the final distribution is sensitive to the nn scattering length, it depends only weakly on the effective range. Throughout we emphasize the impact of theoretical uncertainties on the neutron-neutron relative-energy distribution, and discuss the extent to which those uncertainties limit the extraction of the neutron-neutron scattering length from the reaction 6 He(p,pα)nn.

6 ≤ A ≤ 19↗
DOE OSTI2015

Predicting the production of neutron-rich heavy nuclei in multinucleon transfer reactions using a semi-classical model including evaporation and fission competition, GRAZING-F

Multinucleon transfer reactions have recently attracted attention as a possible path to the synthesis of new neutron-rich heavy nuclei. Here, we study transfer reactions involving massive nuclei with the intention of understanding if the semi-classical model GRAZING coupled to an evaporation and fission competition model can satisfactorily reproduce experimental data on transfer reactions in which fission plays a role. We have taken the computer code grazing and have added fission competition to it (grazing-f) using our current understanding of Γ n /Γ f , fission barriers, and level densities. The code grazing-f seems to satisfactorily reproduce experimental data for +1p, +2p, and +3p transfers but has limitations in reproducing measurements of larger above-target and below-target transfers. Nonetheless, we use grazing-f to estimate production rates of neutron-rich N = 126 nuclei, actinides, and transactinides. The grazing code, with appropriate modifications to account for fission decay as well as neutron emission by excited primary fragments, does not predict large cross sections for multinucleon transfer reactions leading to neutron-rich transactinide nuclei but predicts opportunities to produce new neutron-rich actinide isotopes.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2024

Supersymmetry with scalar sequestering

Supersymmetric models with a strongly interacting superconformal hidden sector (HS) may drive soft supersymmetry (SUSY) breaking scalar masses, bilinear soft term B μ and Higgs combinations m H u , d 2 + μ 2 to small values at some intermediate scale, leading to unique sparticle mass spectra along with possibly diminished fine-tuning in spite of a large superpotential μ parameter. We set up a computer code to calculate such spectra, which are then susceptible to a variety of constraints: (1) possible charge-or-color breaking (CCB) minima in the scalar potential, (2) unbounded from below (UFB) scalar potential, (3) improper electroweak symmetry breaking, (4) a charged or sneutrino lightest SUSY particle (LSP), (5) generating m h ∼ 125 GeV , (6) consistency with LHC sparticle mass limits, and (7) naturalness. We find this bevy of constraints leaves little or no viable parameter space for the case where hidden sector dynamics dominates minimal supersymmetric standard model (MSSM) running, even for the case of nonuniversal gaugino masses. For the case with moderate HS running with comparable MSSM running, and with universal gaugino masses, then the fine-tuning is ameliorated, but nonetheless remains high. Viable spectra with moderate HS running and with low fine-tuning and large μ can be found for nonuniversal gaugino masses. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2021

First-principles equation of state database for warm dense matter computation

We put together a first-principles equation of state (FPEOS) database for matter at extreme conditions by combining results from path integral Monte Carlo and density functional molecular dynamics simulations of the elements H, He, B, C, N, O, Ne, Na, Mg, Al, and Si as well as the compounds LiF , B 4 C , BN , CH 4 , CH 2 , C 2 H 3 , CH , C 2 H , MgO , and MgSiO 3 . For all these materials, we provide the pressure and internal energy over a density-temperature range from ~ 0.5 to 50 g cm - 3 and from ~ 10 4 to 10 9 K, which are based on ~ 5000 different first-principles simulations. We compute isobars, adiabats, and shock Hugoniot curves in the regime of L - and K -shell ionization. Invoking the linear mixing approximation, we study the properties of mixtures at high density and temperature. Furthermore, we derive the Hugoniot curves for water and alumina as well as for carbon-oxygen, helium-neon, and CH-silicon mixtures. We predict the maximal shock compression ratios of H 2 O , H 2 O 2 , Al 2 O 3 , CO , and CO 2 to be 4.61, 4.64, 4.64, 4.89, and 4.83, respectively. Finally we use the FPEOS database to determine the points of maximum shock compression for all available binary mixtures. We identify mixtures that reach higher shock compression ratios than their end members. We discuss trends common to all mixtures in pressure-temperature and particle-shock velocity spaces. In the Supplemental Material, we provide all FPEOS tables as well as computer codes for interpolation, Hugoniot calculations, and plots of various thermodynamic functions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Interaction models and configurational entropies of binary MoTa and the MoNbTaW high entropy alloy

We introduce a simplified method to model the interatomic interactions of high entropy alloys based on a lookup table of cluster energies. Furthermore, these interactions are employed in replica exchange Monte Carlo simulations with histogram analysis to obtain thermodynamic properties across a broad temperature range. Kikuchi's cluster variation method entropy formalism is applied to directly calculate entropy from statistics on short- and long-range chemical order, and we discuss the convergence of the entropy as clusters of differing size are included. A high temperature series expansion aids in our understanding of the convergence. Computer codes implementing these methods, and supporting data, are freely available on the internet.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Evaluating Recombination Mechanisms in RbF Treated Cu(In${}_\mathrm{x}$Ga$_\mathrm{1-x}$)Se$_{2}$ Solar Cells

Rubidium fluoride (RbF) postdeposition treatment (PDT) has been shown to improve the performance of Cu(In x Ga 1-x )Se 2 (CIGS) photovoltaic devices. Here, in this study, temperature-dependent current voltage (JVT) and time-resolved photoluminescence (TRPL) experiments were combined with modeling using the solar cell capacitance simulator (SCAPS) computer code to investigate the effect of the RbF PDT. Two devices, one as-deposited and one with RbF PDT, were deposited by a three stage coevaporation process. JVT measurements suggest the dominant recombination mechanism may be tunneling-enhanced recombination via bandtail states, but that defect states in the bandgap can also be important. RbF PDT is shown to decrease the characteristic energy of the bandtails. TRPL data show an increase in the minority carrier lifetime after RbF PDT, leading to an improved open-circuit voltage. SCAPS modeling indicates that the dominant recombination mechanism is dependent on the specific defect makeup of a device, suggesting that small changes in processing conditions can impact device behavior. This explains the observation that, for some devices, defect states in the gap dominate while others, as is the case here, appear to be dominated by bandtails.

14 SOLAR ENERGY↗
DOE OSTI2024

Co-design for Particle Applications at Exascale

Co-design across the Exascale Computing Project (ECP) has been critical for both enabling science applications and bringing disparate communities together. Developing and porting applications to the various high-performance computing (HPC) architectures on pre-exascale and exascale computers has been quite challenging due to the diversity of hardware features and software stacks. The Co-design Center for Particle Applications (CoPA) has developed and enhanced the Cabana and PROGRESS/BML libraries to facilitate the creation of new particle applications, make existing particle applications exascale capable, and allow teams to explore new capabilities. Particle methods from atomistic, mesoscale, continuum, through cosmological scales have been built with Cabana, along with new possibilities for application coupling. Similarly, the PROGRESS/BML library has enabled quantum particle applications with linear algebra solvers to use advanced hardware. Across these CoPA-developed libraries, the co-design abstraction layer combines performance portability with math library support to facilitate separation of concerns and directly support science runs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Rapid and semi-analytical design and simulation of a toroidal magnet made with YBCO and MgB 2 superconductors

Recent advances in second-generation (YBCO) high-temperature superconducting wire could potentially enable the design of super high performance energy storage devices that combine the high energy density of chemical storage with the high power of superconducting magnetic storage. However, the high aspect ratio and the considerable filament size of these wires require the concomitant development of dedicated optimization methods that account for the critical current density in type-II superconductors. In this study, we report on the novel application and results of a CPU-efficient semianalytical computer code based on the Radia 3-D magnetostatics software package. Our algorithm is used to simulate and optimize the energy density of a superconducting magnetic energy storage device model, based on design constraints, such as overall size and number of coils. The rapid performance of the code is pivoted on analytical calculations of the magnetic field based on an efficient implementation of the Biot-Savart law for a large variety of 3-D “base” geometries in the Radia package. The significantly reduced CPU time and simple data input in conjunction with the consideration of realistic input variables, such as material-specific, temperature, and magnetic-field-dependent critical current densities, have enabled the Radia -based algorithm to outperform finite-element approaches in CPU time at the same accuracy levels. Comparative simulations of MgB 2 and YBCO-based devices are performed at 4.2 K, in order to ascertain the realistic efficiency of the design configurations.

2G↗
DOE OSTI2025

Serpentine Magnet Designs for the Interaction Region of the Electron-Ion Collider (EIC)

The Electron-Ion Collider (EIC), hosted by Brookhaven National Laboratory, is designed to deliver a peak luminosity of 1 × 10 34 cm −2 sec −1 . The interaction region (IR) of the EIC imposes several constraints in terms of field quality, aperture, and spatial layout, which necessitates the development of several unique superconducting serpentine direct wind magnets. These magnets are constructed using either a single strand or a small-diameter 6-around-1 NbTi cable, presenting unique challenges for design and optimization. This paper introduces a new computational code specifically developed to streamline and integrate the design process for these magnets, enabling faster design iterations while addressing their complex requirements. Here, in this paper, we first introduce the code, which builds on established electromagnetic fundamentals. The code incorporates tools for optimizing winding patterns and for correcting magnetic multipoles; additionally, it interfaces with established magnet design software. We also present the design of several serpentine magnets for the EIC IR, demonstrating the code’s capability to deliver precise and efficient solutions. These designs highlight the code’s ability to accelerate the development cycle, ensuring the serpentine magnets meet the demanding specifications of the EIC project.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Same Data, Different Audiences: Using Personas to Scope a Supercomputing Job Queue Visualization

Domain-specific visualizations sometimes focus on narrow, albeit important, tasks for one group of users. This focus limits the utility of a visualization to other groups working with the same data. While tasks elicited from other groups can present a design pitfall if not disambiguated, they also present a design opportunity—namely, the development of visualizations that support multiple groups. This development choice presents a trade-off of broadening the scope but limiting support for the more narrow tasks of any one group, which in some cases can enhance the overall utility of the visualization. We investigate this scenario through a design study where we develop Guidepost, a notebook-embedded visualization of data that helps scientists assess compute wait times, machine learning researchers understand prediction accuracy, and system maintainers analyze usage trends. We adapt the use of personas for visualization design from existing literature in the HCI and design domains, applying them to categorize tasks based on their uniqueness across stakeholder personas. Under this model, tasks shared between all groups should be supported by interactive visualizations and tasks unique to each group can be deferred to scripting with notebook-embedded visualization design. We evaluate our visualization through real-world case studies and a task-focused evaluation with nine participants. We observe that together, Guidepost's visual encodings, interactions, and export capabilities support the tasks of our differing personas.

97 MATHEMATICS AND COMPUTING↗