Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computer Modeling”

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

Generating forces in confinement via polymerization

Understanding how to produce forces using biomolecular building blocks is essential for the development of adaptive synthetic cells and living materials. Here we ask whether a dynamic polymer system can generate deformation forces in soft shells by pure self-assembly, motivated by the fact that biological polymer networks like the cytoskeleton can exert forces, move objects, and deform membranes by simply growing, even in the absence of molecular motors. We address this question by investigating polymer force generation by varying the release rate, the structure, and the interactions of self-assembling monomers. First, we develop a toy computational model of polymerization in a soft elastic shell that reveals the emergence of spontaneous bundling which enhances shell deformation. We then extend our model to account more explicitly for monomer binding dynamics. We find that the rate at which monomers are released into the interior of the shell is a crucial parameter for achieving deformation through polymer growth. Lastly, we demonstrate that the introduction of multivalent particles that can join polymers can either improve or impede polymer performance, depending on the amount and on the structure of the multivalent particles. Our results provide guidance for the experimental realization of polymer systems that can perform work at the nanoscale, for example through rationally designed self-assembling proteins or nucleic acids.

Osmanović, Dino [University of California at Los A↗
DOE OSTI2026

Accelerating catalytic advancements through the precision of high-throughput experiments & calculations

The growing demand for energy-efficient processes to support a sustainable future drives the need for research to rapidly explore chemical and material space through accelerated catalyst discovery initiatives. Recent breakthroughs in high-throughput experimental and computational methods are transforming the catalysis field, surpassing traditional approaches to manipulating variables in catalytic processes. Key advancements in innovation include the integration of machine learning for efficient catalyst screening, high-throughput experimentation, data-driven methodologies employing comprehensive databases, and in situ and in operando techniques for realistic observations. This progress has undoubtedly been intertwined with a collaborative framework across disciplines, reshaping catalyst discovery methods in both industry and academia. This Opinion article presents a multifaceted perspective from coauthors with expertise spanning various stages of the Technology Readiness Level spectrum, highlighting both opportunities and persistent challenges in integrating computational and experimental approaches in catalysis. These challenges span from obtaining high-quality experimental data, scaling simulations to industrially relevant materials and process conditions to navigating the complexity and predictive accuracy of computational models.

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

Identifying mechanistic differences between co-fed CO 2 hydrogenation and reactive CO 2 capture using Ru and Pd dual function materials

Dual function materials (DFMs) enable reactive carbon capture (RCkeC), an intensified approach to carbon dioxide capture and utilization for cost and energy input reductions. Yet, there is a fundamental lack of understanding of mechanisms around CO 2 adsorption and subsequent conversion on these materials, hindering further development. Herein, we investigated several supported alkaline metal oxides for their CO 2 adsorption characteristics to find that Na/Al 2 O 3 had the highest CO 2 adsorption capacity, accompanied by a variety of CO 2 adsorption geometries as identified by in situ DRIFTS and computational modeling. The addition of catalytic metals (Ru, Pd) increased the adsorption capacity of Na/Al 2 O 3 without altering binding modes. In the subsequent reactive desorption step, acetate and formate intermediates were observed. Notably, this mechanistic investigation identified that the formation of acetate species was unique to RCC on a DFM, as these species were not observed in co-fed hydrogenation over the DFM or RCC over a Na-free catalyst.

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

Illuminating the mechanistic impacts of an Fe-quaterpyridine functionalized crystalline poly(triazine imide) semiconductor for photocatalytic CO 2 reduction

The strategy of incorporating earth-abundant catalytic centers into light-absorbing architectures is desirable from the viewpoint of low cost, low toxicity, and versatility for activating small molecules to produce solar-based fuels. Herein, we show that an Fe-quaterpyridine molecular catalyst can be anchored to a light-absorbing, crystalline, carbon nitride (PTI), to yield a molecular-catalyst/material hybrid, Fe-qpy-PTI, capable of facilitating CO 2 reduction to CO selectively (up to ∼97–98%) in aqueous solution under low-intensity light irradiation. This hybrid material leverages the ability of the Fe-qpy catalyst to bind CO 2 upon a one-electron reduction, as achieved by transfer of excited electrons from the carbon-nitride semiconductor. At a low incident power density of only 50 mW cm −2 , the catalytic activity of the hybrid material was measured across a range of catalyst loadings from 0.1–3.8 wt%, yielding CO rates of up to 596 μmol g −1 h −1 for a 3.8 wt% loading during a 3 h experiment. Over the course of 8 h, the hybrid material attained a CO evolution rate of 608 μmol g −1 h −1 and 305 turnovers for a TOF of ∼38 h −1 and an apparent quantum yield of 2.6%. Higher light intensities provided an initial increase in activity but negatively impacted photocatalytic rates with time, with an AQY of 0.6% at 150 mW cm −2 and 0.4% at 250 mW cm −2 . Transient absorption spectroscopy results showed electron survival probabilities consistent with the trends in observed product rates. Computational modeling was also used to evaluate and understand the mechanistic pathway of the high product selectivity for CO versus H 2 . These results thus help unveil key factors for leveraging the mechanistic understanding of molecular catalysts for CO 2 reduction for pairing with light absorbing semiconductors and establishing optimal conditions to attain maximal rates in aqueous solution.

McGuigan, Scott [Baylor University, Waco, TX (Unit↗
DOE OSTI2025

Subsurface hydrogen, curvature, and strain: lessons from electro-reduction of benzaldehyde on nano-structured Pd catalysts

The unique ability of palladium (Pd) to absorb hydrogen and form a bulk hydride is vital for chemical transformations that involve hydrogenation reactions. Nano-structured Pd catalysts offer a promise of tuning these reaction rates by exploiting variations of reactant binding energies depending on the surface structure and morphological constraints that result in inhomogeneous strain. However, the interplay between the nano-structure of Pd and the ability of Pd to adsorb (and absorb) hydrogen as well as other reactive species needs to be better understood for a rational understanding of competitive chemical transformations at Pd surfaces. We consider the effects of the surface corrugation, strain, and subsurface Pd hydride on the reduction of benzaldehyde to benzyl alcohol in two qualitatively different samples – Pd nanoparticles and Pd gels formed by quasi-one-dimensional chains of these nanoparticles. Our electrochemical measurements and computational modelling suggest that surface concave sites, inherent to Pd gels, facilitate hydrogen transfer to the Pd subsurface region, thus weakening benzaldehyde binding to the surface. This effect is further modulated by the strain, depending on the local coordination environment on the corrugated surface. Furthermore, these findings demonstrate how structurally complex samples in the form of gels provide degrees of freedom for controlling the behavior of metal catalysts that are not available in isolated nanoparticles, which paves the way for new approaches in the design of catalytic materials and synthesis of metal hydrides.

Padavala, Sri Krishna Murthy [University of Minnes↗
DOE OSTI2024

Formation of a spherical plasma liner for plasma-jet-driven magneto-inertial fusion

Plasma-jet-driven magneto-inertial fusion is an alternative approach to controlled nuclear fusion, which aims to utilize a line-replaceable dense plasma liner as a repetitive spherical compression driver. In this experiment, first measurements of the formation of a spherical argon plasma liner formed from 36 discrete pulsed plasma jets are obtained on the Plasma Liner Experiment. Properties including liner uniformity and morphology, plasma density, temperature, and ram pressure are assessed as a function of time throughout the implosion process and indicate an apparent transition from initial kinetic inter-jet interpenetration to collisional regime near stagnation times, in accordance with theoretical expectation. A lack of primary shock structures between adjacent jets during flight implies that arbitrarily smooth liners may be formed by way of corresponding improvements in jet parameters and control. The measurements facilitate the benchmarking of computational models and understanding the scaling of plasma liners toward fusion-relevant energy density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

An updated technique to obtain explosive kinetics data on microsecond timescales

There are few techniques available for chemists to obtain time-to-explosion data with known temperature inputs at the early stages of the design and synthesis of new explosives. In the 1960s, a technique was developed to rapidly heat milligram-quantities of confined explosives to ~1000 K on microsecond timescales. Wenograd loaded explosives inside stainless steel hypodermic needles, connected them to a fireset and rapidly discharged a capacitor through the steel. He obtained the temperature by measuring the needle resistance in a Wheatstone bridge arrangement and the time to explosion from a needle rupture. However, owing to the narrow-gauge needles used in the original research, the experiment was only possible with melt-castable explosives; it was never replicated, and modern diagnostics are now available with advances beyond the 1960s. Here, we report the development of the High Explosives Initiation Time (HEIT) test, which utilizes a 250 J pulsed power system to heat the needles. This work extends the Wenograd approach by using optical diagnostics, computational modeling, and advanced techniques to measure needle resistance and needle rupture. Preliminary rate information for pentaerythritol tetranitrate (PETN) will be presented.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2025

Effect of particle shape on stratification in drying films of binary colloidal mixtures

The role of particle shape in evaporation-induced auto-stratification in polydisperse colloidal suspensions is explored with molecular dynamics simulations of mixtures of spheres and aspherical particles. A unified framework based on the competition between diffusion and diffusiophoresis is proposed to understand the effects of shape and size dispersity. In general, particles diffusing more slowly (e.g., larger particles) tend to accumulate more strongly at the evaporation front. However, larger particles have larger surface areas and therefore greater diffusiophoretic mobility. Hence, they are more likely to be driven away from the evaporation front via diffusiophoresis. For a rapidly dried bidisperse suspension containing small and large spheres, the competition leads to “small-on-top” stratification. Here, we employ a computational model in which the diffusion coefficient is inversely proportional to particle mass. For a mixture of spheres and aspherical particles with similar mass, the diffusion contrast is reduced, and the spheres are always enriched at the evaporation front as they have the smallest surface area for a given mass and, therefore, the lowest diffusiophoretic mobility. Furthermore, for a mixture of solid and hollow spheres that have the same outer radius and thus the same surface area, the diffusiophoretic contrast is suppressed, and the system is dominated by diffusion. Consequently, the solid spheres, which have a larger mass and diffuse more slowly, accumulate on top of the hollow spheres. Finally, for a mixture of thin disks and long rods that differ significantly in shape but have similar mass and surface area, both diffusion and diffusiophoresis contrasts are suppressed, and the mixture does not stratify.

Classical molecular dynamic simulations↗
DOE OSTI2025

Interpreting dynamic-compression experiments to uncover the time dependence of freezing: Application to gallium

Using pulsed-power magnetic field sources to compress gallium to gigapascal pressures on nanosecond timescales, we report here experiments on shockless dynamic compression of a liquid metal. Time-resolved velocimetry data reveal signatures of rapid freezing from a metastable liquid state, and we demonstrate that the kinetics of this nonequilibrium solidification can be accurately simulated with a computational modeling framework we have developed in previous studies, where classical nucleation theory is coupled with hydrodynamics. Notably, velocity traces in some of our experiments show evidence of a phase transition, while others do not, even though other types of evidence suggest that solidification may be occurring in all of them. We explain how predictions made by our models regarding the presence or absence of these phase-transition signatures motivated additional experiments that later confirmed the theoretical predictions. Our analysis shows that due to the rapid, quasi-isentropic nature of the loading path, our experiments were able to compress liquid gallium to metastable states that are undercooled below the equilibrium melt temperature by more than 300 K and exhibit pressures that approach five times the equilibrium melt pressure. The understanding gained in this study should form the basis for future dynamic-compression experiments aimed at interrogating melt curves at high pressures.

gallium↗
DOE OSTI2025

Design of a High-Assay Low-Enriched Uranium Tri-Structural Isotropic Critical Experiment for Advanced Reactor Validation

High-assay low-enriched uranium (HALEU) fuel is a key component of many small modular reactor designs. Critical experiments are an important way to understand the neutronic performance of systems by obtaining nuclear data validations through measurements. Data reduce uncertainty and risk by showing that systems respond as predicted to changes such as temperature, subsequently advancing the overall technology readiness level of the materials within. Numerous critical experiments have been performed at the National Criticality Experiments Research Center (NCERC) operated by Los Alamos National Laboratory at the Nevada National Security Site since it became operational in 2011. However, the first experiment with HALEU fuel did not occur until 2024. Through extensive engineering, the experiment described in this paper was successfully designed and executed for the Comet vertical lift assembly at NCERC to perform measurements with HALEU tri-structural isotropic fuel that will assist in validation of nuclear data and computational modeling of small modular reactors for years to come.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

The impact of optical measurement techniques on measured aerosol particle size distributions

Ambient aerosol particle size distributions measured by the Ultra-High Sensitivity Aerosol Spectrometer (UHSAS) at various sites around the world exhibit modes at optical diameters near 600 nm and 850 nm. These modes are not present in concurrent measurements with the Grimm 11-D Optical Particle Counter (OPC). Here, in this study, we argue that these modes result from the optical measurement technique itself, and we explain why they appear in measurements by the UHSAS but not in those by the Grimm OPC. We construct computer models of the UHSAS and Grimm (“digital UHSAS” and “digital Grimm”) and use these to investigate the size distribution that would result from measurements of artificial aerosol particle size distributions that do not contain modes. The appearance of modes for the structureless incoming size distributions sampled by the digital UHSAS is explained by the nonlinear behavior of partial scattering cross sections of uniform spherical particles as a function of their diameter. The absence of modes in the digital Grimm is explained by the coarser size resolution of that instrument. Detailed analysis of the relationship between optical and geometric diameters for uniform spherical particles reveals two important results. First, these diameters generally have different numerical values for the same particle, and second, the relationship is nonlinear; thus, widths of size bins in terms of optical diameter differ from those in geometric diameter. These results explain the modes observed in the ambient size distributions and highlight concerns with attempts to create a merged size distribution by combining measurements from different instruments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Wall heating by subcritical energetic electrons generated by the runaway electron avalanche source *

Abstract Subcritical energetic electrons (SEEs) produced by the runaway electron (RE) avalanche source at energies below the runaway threshold are found to be the primary contributor to surface heating of plasma-facing components (PFCs) during final loss events. This finding is supported by theoretical analysis, computational modeling with the Kinetic Orbit Runaway electrons Code (KORC), and qualitative agreement with DIII-D experimental observations. The avalanche source generates significantly more secondary electrons below the runaway threshold, which thermalize rapidly when well-confined. However, during a final loss event, the RE beam impacts the first wall, and SEEs are deconfined before they can thermalize. Additionally, because the energy deposition length decreases faster than energy, the deposited energy density, and thus the maximum PFC surface temperature change, is larger for SEEs than REs. KORC simulations employ an analytic first wall to model particle deconfinement onto a non-axisymmetric wall composed of individual tiles. PFC surface heating is calculated using a 1D model extended to include an energy-dependent deposition length scale. Simulations of DIII-D qualitatively agree with infrared (IR) imaging only when SEEs from the avalanche source are included. These results demonstrate that SEEs are the dominant contributor to PFC surface heating and indicate that the avalanche source plays a critical role in the PFC damage caused during final loss events. The prominence of SEEs also has important implications for interpreting IR imaging, one of the primary diagnostics for RE-wall interaction diagnosis, despite REs dominating the energy and current density. This result improves predictions of wall damage due to post-disruption REs to estimate material lifetime and design RE mitigation systems for ITER and future reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Mapping causal pathways with structural modes fingerprint for perovskite oxides

Abstract Causality is innate to the determination of the fundamental mechanism controlling any physical phenomena. However, combining causality within the standard practices of computational modelling to understand structure-functionality connections is extremely rare. This work proposes a fingerprint based on key structural modes for ABO 3 -type perovskite oxides and its derivatives, combined with causal models, for predicting Kohn–Sham energies. Our study of causal models captures the inherent coupling between structural modes such as rotation, tilt and antiferroelectric displacements, responsible for phase transition, polarization, magnetization and metal–insulator transition, exhibited by these materials. Although developed for modelling specific functionality, this method is universally applicable to derive other functionalities and even different material classes while tracking hidden causal mechanisms via structural distortions.

42 ENGINEERING↗
DOE OSTI2025

The impact of curation errors in the PDBBind Database on machine learning predictions of protein–protein binding affinity

The PDBBind database has been widely utilized for the computational prediction of protein–protein binding affinities. While the accuracy of the PDBBind-curated equilibrium dissociation constants (K D ) has been reported for the protein–ligand subset of the PDBBind database, the curation accuracy has not been reported for the protein–protein subset. Here, we present a detailed manual analysis for the subset of PDBBind records with PubMed Central Open Access primary publications and find that ~19% of these records had K D values that were not supported by their primary publications. The impact of these putative curation errors on the machine learning-based prediction of K D from experimental protein–protein 3D structures was evaluated and correcting the curation errors improved the Pearson correlation coefficient between measured and random forest-predicted log 10 (K D ) values by ~8 percentage points. This finding underscores the importance of dataset accuracy for computational modelling and highlights the need for more stringent curation processes when extracting information from the scientific literature.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Three-dimensional lattice modulations in the charge density wave system Lu 2 Ir 3 Si 5

Using total and resonant x-ray scattering coupled to large-scale computer modeling, we study the lattice modulations in the complex charge density wave (CDW) material Lu 2 ⁢Ir 3 ⁢Si 5 . Here, we find that it is a unique quantum system where periodic lattice modulations related to emergent CDW order occur in three orthogonal atomic planes of the crystal lattice, leading to the emergence of an unusual three-dimensional (3D) pattern of short and long Ir-Ir and Lu-Lu bonds. The 3D character of observed lattice modulations explains the largely isotropic character of the changes in the electronic properties occurring when the CDW order sets in, demonstrating the strong electron-lattice coupling in Lu 2 ⁢Ir 3 ⁢Si 5 . The result is supported by DFT calculations based on the experimental structure data. Altogether, our work provides strong evidence for the presence of a relationship between the dimensionality of emergent lattice distortions and that of concurrent changes in the electronic properties of CDW materials. The relationship may need to be accounted for when these materials are explored for practical applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Self-assembled reconfigurable pump architectures via magnetic colloidal swarms

Self-assembled swarms of interactive active units, which are adaptive and dynamically reconfigurable to accommodate different functionalities, represent a promising platform for the development of next-generation robotics. Here, we utilize the emergent collective behavior of active magnetic colloids confined in quasi-two-dimensional arrays of overlapping wells to demonstrate the self-organization of a colloidal swarm into a dynamic pump architecture capable of controlled transport of passive cargo particles. This dynamic architecture provides a global unidirectional looping flow pattern along the entire length of the system. We show that the flow direction of the dynamic swarm-based pump can be externally controlled by a phase shift of a driving magnetic field energizing the swarm. The experimental observations are supported by computational modeling based on phenomenological coarse-grained particle dynamics coupled to shallow-water Navier-Stokes hydrodynamics. In conclusion, our findings demonstrate how the emergent collective behavior of a swarm can be orchestrated into a desired functionality by exploiting the interplay between activity and confinement potentials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Single spin asymmetry in forward p A collisions from the Pomeron-odderon interference

Working in the hybrid framework of the high energy p A collisions we identify a new contribution to transverse single spin asymmetry (SSA). The phase necessary for the SSA is provided by the Pomeron-odderon interference in the dense nuclear target. The complete formula for the p A → h X polarized cross section also contains the transversity distribution for the polarized projectile as well as the real part of the twist-three fragmentation function. We numerically estimate the asymmetry A N and its nuclear dependence. Based on a model computation we find that A N can be a percent level in the forward and low- P h ⊥ region. For large nuclei we find significant suppression, with A N ∝ A − 7 / 6 parametrically. As a notable feature we find a node of A N as a function of the P h ⊥ around the values of the initial saturation scale that could be used to test this mechanism experimentally. Published by the American Physical Society 2025

Benić, Sanjin↗
DOE OSTI2025

Structural and electronic changes in L⁢i 2 ⁢Ru⁢O 3 induced by lithium intercalation

Despite extensive research on oxide battery cathodes that transcend classical cationic redox activity, the detailed interplay between structural transformations and electronic redox processes remains insufficiently understood. We report a detailed study of the sequential structural and electronic changes in Li 2 RuO 3 upon lithium intercalation, characterized by powder x-ray and neutron diffraction alongside Ru and O K-edge x-ray absorption spectroscopy (XAS), and guided by operando synchrotron x-ray diffraction. During delithiation, Li 2 RuO 3 evolves from a well-defined monoclinic state to a complex trigonal phase via multiple intermediate structures, marked by significant changes in Ru-O bond distances that closely track the transition from a classical cationic redox to an unconventional process centered at oxygen states. Armed with high-quality atomic structural descriptions, computational models of the O K-edge XAS closely reproduce the experimentally observed spectral shifts. Lastly, we relate observations of electrochemical hysteresis with concurrent changes in the pathways of structural and electronic transitions. In conclusion, our results not only clarify the mechanisms underpinning voltage hysteresis in a model for lattice oxygen redox but also underscore the importance of structural fidelity in modeling redox behavior when this type of complex reactivity is present.

Li, Haifeng [Univ. of Illinois, Chicago, IL (Unite↗