Search NASA⌕ Search

SEARCH · Search NASA

Results for “computing methodologies”

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 73 records · Page 4

Continuation of probability density functions using a generalized Lyapunov approach

Techniques from numerical bifurcation theory are very useful to study transitions between steady fluid flow patterns and the instabilities involved. Here, we provide computational methodology to use parameter continuation in determining probability density functions of systems of stochastic partial differential equations near fixed points, under a small noise approximation. Key innovation is the efficient solution of a generalized Lyapunov equation using an iterative method involving low-rank approximations. We apply and illustrate the capabilities of the method using a problem in physical oceanography, i.e. the occurrence of multiple steady states of the Atlantic Ocean circulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

An approach to combine neutron and ion irradiation data to accelerate material qualification for nuclear reactors

Next-generation nuclear power plants are generally characterized by higher operating temperatures, increased neutron fluences and energies, and distinct corrosive coolant environments versus the existing light water reactor fleet. Whether using existing materials in new environments, newly developed materials tailored for these environments, or new manufacturing methods, the traditional decades-long approach for materials qualification does not facilitate rapid deployment. Ion irradiation has demonstrated success in reproducing material microstructure and select property evolution resulting from neutron irradiation with three to four orders of magnitude reduction in time and cost, making it an ideal candidate for accelerated irradiation testing. Because microstructure has a large impact on bulk material properties, limited neutron irradiation data at lower damage levels can in principle be combined with accelerated ion testing results and modeling and simulation to form an accurate prediction of microstructure evolution and select properties under different neutron irradiation conditions and at higher damage levels. The objective of this work is to present a conceptual framework of specific steps to fulfill several technical challenges associated with qualifying materials for performance in radiation environments on an accelerated time frame. A brief review of the regulatory landscape for materials in nuclear environments is presented, followed by additional overviews to understand the current state of the art for correlation of materials properties across radiation environments using experimental and computational methodologies. Finally, the roles of academia, national laboratories, and industry in the advancement of this accelerated materials qualification framework are discussed as a path forward, with possible case studies presented.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Coupled neutronics, thermochemistry, corrosion modeling and sensitivity analyses for isotopic evolution in molten salt reactors

This study presents a computational methodology for analyzing isotopic evolution and associated uncertainties in molten salt reactors (MSRs), focusing on both fluoride- and chloride-based fuel salts. The primary goal is to enhance the understanding of isotopic behavior in MSRs and provide data to support future experimental efforts. The methodology integrates transport-coupled depletion calculations using OpenMC, equilibrium thermodynamics modeling with Thermochimica, and a corrosion model. Sensitivity analyses are performed to evaluate the impact of power density, air ingress, and humidity content on isotopic evolution in MSR concepts. This study examines representative F- and Cl-based MSR designs, highlighting the dominant influence of power density on isotopic composition, which significantly affects isotope production and depletion rates, accounting for approximately 76% of the observed variance in element concentration. Air ingress and humidity content also affect the redox potential, solubility of heavier elements, and corrosion rates, thereby altering the expected isotopic evolution in the reactor. On average, air ingress accounts for around 17% of the variance in element concentrations, while humidity explains the remaining 7%. These variances differ significantly from element to element, depending on the element’s role in depletion, redox potential evolution, and galvanic corrosion. The findings indicate that power density, air ingress, and humidity content are all critical factors for optimizing reactor design and operational strategies. Furthermore, the study provides expected ranges for key impurities in the fuel salt, which are crucial for guiding future experimental studies and refining MSR designs. Finally, this study demonstrates the importance of modeling depletion coupled with the evolution of redox potential and chemical interactions in MSR fuel salts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Site Disorder Drives Cyanide Dynamics and Fast Ion Transport in Li 6 PS 5 CN

Halide argyrodite solid-state electrolytes of the general formula Li 6 PS 5 X exhibit complex static and dynamic disorder that plays a crucial role in ion transport processes. Here, we unravel the rich interplay between site disorder and dynamics in the plastic crystal argyrodite Li 6 PS 5 CN and the impact on ion diffusion processes through a suite of experimental and computational methodologies, including temperature-dependent synchrotron powder X-ray diffraction, AC electrochemical impedance spectroscopy, 7 Li solid-state NMR, and machine learning-assisted molecular dynamics simulations. Sulfide and (pseudo)halide site disorder between the two anion sublattices unilaterally improves long-range lithium diffusion irrespective of the (pseudo)halide identity, which demonstrates the importance of site disorder in dictating bulk ionic conductivity in the argyrodite family. Furthermore, we find that anion site disorder modulates the presence and time scales of cyanide rotational dynamics. Ordered configurations of anions enable fast, quasi-free rotations of cyanides that occur on time scales of 10 11 Hz at T = 300 K. In contrast, we find that cyanide dynamics are slow or frozen in Li 6 PS 5 CN when site disorder between the cyanide and sulfide sublattices is present at T = 300 K. We rationalize the observed differences in cyanide dynamics in the context of elastic dipole interactions between neighboring cyanide anions and local strain induced by the configurations of site disorder that may impact the energetic landscape for cyanide rotational dynamics. Through this study, we find that anion disorder plays a decisive role in dictating the extent and time scales of both lithium ion and cyanide dynamics in Li 6 PS 5 CN.

36 MATERIALS SCIENCE↗

Importance of Electron Correlation on the Geometry and Electronic Structure of [2Fe–2S] Systems: A Benchmark Study of the [Fe 2 S 2 (SCH 3 ) 4 ] 2–,3–,4– , [Fe 2 S 2 (SCys) 4 ] 2– , [Fe 2 S 2 (S- p -tol) 4 ] 2– , and [Fe 2 S 2 (S- o -xyl) 4 ] 2– Complexes

Iron–sulfur clusters are crucial for biological electron transport and catalysis. Obtaining accurate geometries, energetics, manifolds of their excited electronic states, and reduction energies is important to understand their role in these processes. Here, using a [2Fe–2S] model complex with Fe II and Fe III oxidation states, which leads to different charges, i.e., [Fe 2 S 2 (SMe) 4 ] 2–,3–,4– , we benchmarked a variety of computational methodologies ranging from density functional theory (DFT) to post-Hartree–Fock methods, including complete active space self-consistent field (CASSCF), multireference configuration interaction, the second-order N-electron valence state perturbation theory (NEVPT2), and the linearized integrand approximation of adiabatic connection (AC0) approaches. Additionally, we studied three experimentally well-characterized complexes, [Fe 2 S 2 (SCys) 4 ] 2– , [Fe 2 S 2 (S-o-tol) 4 ] 2– , and [Fe 2 S 2 (S-o-xyl) 4 ] 2– , via DFT methods. We conclude that the dynamic electron correlation is important for accurately predicting the geometry of these complexes. Broken symmetry (BS) DFT correctly predicts experimental geometries of low-spin multiplicity, while CASSCF does not. However, BS-DFT significantly overestimates the difference between the low- and high-spin electronic states for a given oxidation state. At the same time, CASSCF underestimates it but provides relative energies closer to the reference NEVPT2 results. Finally, AC0 provides energetics of NEVPT2 quality with the additional advantage of being able to use large CASSCF sizes. NEVPT2 gives the best estimates of the Fe III /Fe III → Fe II /Fe III (4.27 eV) and Fe II /F III → Fe II /F II (7.72 eV) reduction energies. The results provide insight into the electronic structure of these complexes and assist in the understanding of their physical properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Water–Hydrocarbon Interactions in Anionic Pyrene Monohydrate

Interactions between water and polycyclic aromatic hydrocarbons are essential in many aspects of chemistry, from interstellar and atmospheric processes to interfacial hydrophobicity and wetting phenomena. Despite their growing importance, the intermolecular potentials of the water-hydrocarbon interactions are underdeveloped compared to water-water potentials, and there are similarly few experimental probes that are sensitive to the details of the water-hydrocarbon potential. We present a combined experimental and computational study of anionic pyrene monohydrate, one of the simplest water/hydrocarbon clusters. The action spectrum in the OH region of the mass-selected cluster ion provides a rigorous benchmark for intermolecular potentials and computational methodologies. We identify missing intermolecular interactions and shortcomings in conventional dynamics calculations by comparing experimental data to density functional theory and classical molecular dynamics calculations. Kinetic trapping is prevalent, even for one water molecule and one pyrene molecule, leading to slow equilibration in conventional molecular dynamics calculations, even on nanosecond timescales and at low temperatures (50 K). At constant energy, temperature fluctuations for the pair of molecules are substantial. Immersing the system in a bath of soft spheres and employing parallel tempering alleviates kinetic trapping and dampens temperature fluctuations, bringing the system closer to the thermodynamic limit. With such augmented sampling, a simple, flexible water model reproduces the linewidth and the asymmetric broadening of the symmetric OH stretching mode, which we assign to spectral diffusion. In the OH stretching region, dynamics calculations predict a more intense antisymmetric peak than experiments observe but do not predict the bimodal split symmetric peak that the experiments show. Furthermore, our work suggests that electronic polarization, missing in the empirical force field, is responsible for the first discrepancy and that quantum nuclear effects, captured neither in density functional theory nor in classical dynamics, may be responsible for the second.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Minimization of Cathode|Solid-Electrolyte Interfacial Delamination through the Application of Interphase Layers

Next-generation lithium-ion batteries are expected to use solid electrolytes (SEs) to enable higher energy density and extreme fast-charge capabilities. One major mode of degradation at the cathode|SE interface is delamination between the cathode active materials and SEs, which leads to performance decay. Experimental observations indicate that implementation of interphase layers can minimize the cathode|SE delamination induced capacity fade. A multiscale computational methodology is developed here to investigate the applicability of boron substituted lithium carbonate (Li 2+x B x C 1–x O 3 , x = 0.5, or LBCO) to minimize the delamination at the cathode|SE interface. Atomistic simulations indicate that the fracture energies at both the cathode|LBCO and LBCO|SE interfaces are higher than those at the cathode|SE interface, which reduces the extent of delamination. Mesoscale simulations indicate that, apart from increasing the fracture energy, decreasing the evolution of strain energy by lowering the elastic modulus of the interphase layer can also minimize the extent of delamination at the cathode|SE interface. However, the adoption of an interphase layer with high ionic conductivity is necessary to minimize the ohmic losses during operation at higher current densities. This study provides guidance on selecting interphase layers with specific properties and thicknesses to minimize both interfacial delamination and impedance growth.

LBCO↗

Unveiling Spatial and Temporal Dynamics of Plasmon-Enhanced Localized Fields in Metallic Nanoframes through Ultrafast Electron Microscopy

Plasmonic nanomaterials, particularly noble metal nanoframes (NFs), are important for applications such as catalysis, biosensing, and energy harvesting due to their ability to enhance localized electric fields and atomic efficiency via localized surface plasmon resonance (LSPR). Yet the fundamental structure-function relationships and plasmonic dynamics of the NFS are difficult to study experimentally and thus far rely predominately on computational methodologies, limiting their utilization. This study leverages the capabilities of ultrafast electron microscopy (UEM), specifically photon-induced near-field electron microscopy (PINEM), to probe the light-matter interactions within plasmonic NF structures. Here, the effects of shape, size, and plasmonic coupling of Pt@Au core-shell NFs on spatial and temporal characteristics of plasmon-enhanced localized electric fields are explored. Importantly, time-resolved PINEM analysis reveals that the plasmonic fields around hexagonal NF prisms exhibit a spatially dependent excitation and decay rate, indicating a nuanced interplay between the spatial geometry of the NF and the temporal evolution of the localized electric field. These results and observations uncover nanophotonic energy transfer dynamics in NFs and highlight their potential for applications in biosensing and photocatalysis.

Plasmonics↗

Electronic Structural Analysis of Copper(II)–TEMPO/ABNO Complexes Provides Evidence for Copper(I)–Oxoammonium Character

Copper/aminoxyl species are proposed as key intermediates in aerobic alcohol oxidation. Several possible electronic structural descriptions of these species are possible, and here we probe this issue by examining four crystallographically characterized Cu/aminoxyl halide complexes by Cu K-edge, Cu L 2,3 - edge, and Cl K-edge X-ray absorption spectroscopy. The mixing coefficients between Cu, aminoxyl, and halide orbitals are determined via these techniques with support from density functional theory. The emergent electronic structure picture reveals that Cu coordination confers appreciable oxoammonium character to the aminoxyl ligand. The computational methodology is extended to one of the putative intermediates invoked in catalytic Cu/aminoxyl-driven alcohol oxidation reactions, with similar findings. On the whole, the results have important implications for the mechanism of alcohol oxidation and the underlying basis for cooperativity in this co- catalyst system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

2DMatPedia, an open computational database of two-dimensional materials from top-down and bottom-up approaches

Two-dimensional (2D) materials have been a hot research topic in the last decade, due to novel fundamental physics in the reduced dimension and appealing applications. Systematic discovery of functional 2D materials has been the focus of many studies. Here, we present a large dataset of 2D materials, with more than 6,000 monolayer structures, obtained from both top-down and bottom-up discovery procedures. First, we screened all bulk materials in the database of Materials Project for layered structures by a topology-based algorithm and theoretically exfoliated them into monolayers. Then, we generated new 2D materials by chemical substitution of elements in known 2D materials by others from the same group in the periodic table. The structural, electronic and energetic properties of these 2D materials are consistently calculated, to provide a starting point for further material screening, data mining, data analysis and artificial intelligence applications. We present the details of computational methodology, data record and technical validation of our publicly available data (http://www.2dmatpedia.org/).

36 MATERIALS SCIENCE↗

High-throughput computation and evaluation of raman spectra

Abstract Raman spectroscopy is used ubiquitously in the characterization of condensed materials, spanning from biomaterials, minerals to polymers, as it provides a unique fingerprint of local bonding and environment. In this work, we design and demonstrate a robust, automatic computational workflow for Raman spectra that employs first-principle calculations based on density functional perturbation theory. A set of computational results are compared to Raman spectra obtained from established experimental databases to estimate the accuracy of the calculated properties across chemical systems and structures. Details regarding the computational methodology and technical validation are presented along with the format of our publicly available data record.

47 OTHER INSTRUMENTATION↗

High-throughput screening of inorganic compounds for the discovery of novel dielectric and optical materials

Dielectrics are an important class of materials that are ubiquitous in modern electronic applications. Even though their properties are important for the performance of devices, the number of compounds with known dielectric constant is on the order of a few hundred. Here, we use Density Functional Perturbation Theory as a way to screen for the dielectric constant and refractive index of materials in a fast and computationally efficient way. Our results constitute the largest dielectric tensors database to date, containing 1,056 compounds. Details regarding the computational methodology and technical validation are presented along with the format of our publicly available data. In addition, we integrate our dataset with the Materials Project allowing users easy access to material properties. Finally, we explain how our dataset and calculation methodology can be used in the search for novel dielectric compounds.

36 MATERIALS SCIENCE↗

Amino acid interacting network in the receptor-binding domain of SARS-CoV-2 spike protein

The relation between amino acid (AA) sequence and biologically active conformation controls the process of polypeptide chains folding into three-dimensional (3d) protein structures. The recent achievements in the resolution achieved in cryo-electron microscopy coupled with improvements in computational methodologies have accelerated the analysis of structures and properties of proteins. However, the detailed interaction between AAs has not been fully elucidated. Herein, we present a de novo method to evaluate inter-amino acid interactions based on the concept of accurately evaluating the amino acid bond pairs (AABP). The results obtained enabled the identification of complex 3d long-range interconnected AA interacting network in proteins. The method is applied to the receptor binding domain (RBD) of the SARS-CoV-2 spike protein. We show that although nearest-neighbor AAs in the primary sequence have large AABP, other nonlocal AAs make substantial contribution to AABP with significant participation of both covalent and hydrogen bonding. Detailed analysis of AABP in RBD reveals the pivotal role they play in sequence conservation with profound implications on residue mutations and for therapeutic drug design. This approach could be easily applied to many other proteins of biomedical interest in life sciences.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Detailed simulations of the first deuterium-tritium-filled double shell implosions on the National Ignition Facility

The first indirectly driven, liquid DT-filled double shell inertial confinement fusion (ICF) implosions have recently been successfully performed on the National Ignition Facility (NIF). Double shells are a class of alternative designs that use a low-Z outer shell to compress a foam cushion that accelerates a high-Z inner shell to efficiently compress a liquid DT core. Double shells are challenging to fabricate, field, and model. Important engineering features enabling double shell fabrication include a fill-tube penetrating all shells and a carefully designed and very narrow (few μm) step-joint in the ablator. Due to the higher density materials involved, high Atwood number instabilities are also important at many material interfaces. In this paper, numerical simulations of double shell implosions using the Los Alamos National Laboratory multi-physics radiation-hydrodynamics code xRAGE will be discussed. An extensive effort has been under way for several years to develop the code capabilities for ICF simulations in a common modeling framework to allow ease of simulation setup and standardization of the computational methodology. This paper will present a wide range of simulation results capturing, quantifying, and comparing the impact of all these degradation mechanisms on implosion performance. Brief comparisons with recent experimental results and suggestions for future improvements will also be discussed. Our results suggest that capsule surface roughness and the step-joint gap have the largest impact on implosion performance. Initial experimental data may suggest that the sensitivity to the step-joint gap could provide the dominant explanation for DT-filled double shell experiments that have been fielded on NIF thus far.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An Artificial Retina Processor for Track Reconstruction at the LHC Crossing Rate

The goal of the INFN-RETINA R&D project is to develop and implement a computational methodology that allows to reconstruct events with a large number (> 100) of charged-particle tracks in pixel and silicon strip detectors at 40 MHz, thus matching the requirements for processing LHC events at the full bunch-crossing frequency. Our approach relies on a parallel pattern-recognition algorithm, dubbed artificial retina, inspired by the early stages of image processing by the brain. In order to demonstrate that a track-processing system based on this algorithm is feasible, we built a sizable prototype of a tracking processor tuned to 3 000 patterns, based on already existing readout boards equipped with Altera Stratix III FPGAs. The detailed geometry and charged-particle activity of a large tracking detector currently in operation are used to assess its performances. We report on the test results with such a prototype.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Ionization potential of radium monofluoride

The ionization potential (IP) of radium monofluoride (RaF) was measured to be 4.969(2)[10] eV, revealing a relativistic enhancement in the series of alkaline earth monofluorides. The results are in agreement with a relativistic coupled-cluster prediction of 4.981(7) eV, incorporating up to quantum electrodynamics corrections. Using the same computational methodology, an improved calculation for the dissociation energy (𝐷 0 ) of 5.54(5) eV is presented. This confirms that RaF joins the group of diatomic molecules for which 𝐷 0 > IP, paving the way for precision control and interrogation of its Rydberg states.

Wilkins, S. G. [Massachusetts Inst. of Technology ↗

Magnetic borophenes from an evolutionary search

A computational methodology based on ab initio evolutionary algorithms and spin-polarized density functional theory was developed to predict two-dimensional magnetic materials. Its application to a model system borophene reveals an unexpected rich magnetism and polymorphism. A metastable borophene with nonzero thickness is an antiferromagnetic semiconductor from first-principles calculations, and can be further tuned into a half-metal by finite electron doping. In this borophene, the buckling and coupling among three atomic layers are not only responsible for magnetism, but also result in an out-of-plane negative Poisson's ratio under uniaxial tension, making it the first elemental material possessing auxetic and magnetic properties simultaneously.

2-dimensional systems↗