SEARCH · Search NASA

Results for “Computer simulation”

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.

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Tunable Nanostructures from Inverse Surfactants

Hierarchical materials in the natural world are often made through the self-assembly of amphiphilic molecules. Achieving similar structural complexity in synthetic materials requires understanding how various molecular parameters affect assembly behavior. In recent years, inverse surfactants─molecules with hydrophobic head groups and hydrophilic macromolecular tails─have been shown to self-assemble into supramolecular assemblies in aqueous solutions that show promise for a number of applications, including drug delivery. Here, we build an understanding of the morphological phase diagram of inverse surfactants using insights from scattering experiments, computer simulations, and statistical mechanics. The scattering and simulation results reveal that changing the headgroup size is an important molecular knob in controlling morphological transitions. The molecular size ratio of the hydrophobic group to the hydrophilic group emerges as a crucial dimensionless quantity in our theory and plays a determining role in setting the micelle structure and the transition from mesoscale to macroscale aggregates. Our minimal theory is able to qualitatively explain the key features of the morphological phase diagram, including the prevalence of fiber-like structures in comparison to spherical and planar micelles. Together, these findings provide a more complete picture of the molecular dependencies of assemblies of inverse surfactants, which we hope may aid in the de novo design of supramolecular structures.

Christakopoulos, Panagiotis [ORNL] (ORCID:00090004↗

DOE OSTI2026

Electrolyte-Dependent, “Microscopically Irreversible” H-Atom Transfer Kinetics of Ce-Based Metal–Organic Framework, Ce-MOF-808

Redox reactions at the interface of metal oxides and protic electrolytes almost always involve protons and electrons in equal amounts. Given the stoichiometry, these proton-coupled electron transfer (PCET) reactions are thermochemically equivalent to net H-atom transfer (HAT) reactions. The correlation between the chemical nature of solid catalysts and HAT kinetics has been employed for decades as the design principle for energy-relevant reactions (e.g., reactions of 2H + /H 2 ). More recently, chemists have experimentally determined that a change in liquid electrolytes that alters the microenvironment at the redox-active sites has an equally profound impact on electrocatalysis involving PCET/HAT. Yet, precise correlations between the chemical nature of electrolytes and the PCET kinetics are, to date, rare in the literature. Herein, we report our findings using the Ce-based metal−organic framework, Ce-MOF-808, as a model system. Each Ce 6 (μ 3 −O) 4 (μ 3 − OH) 4 (OH) 6 (H 2 O) 6 node of this MOF undergoes a 1H + /1e − redox reaction. Using chronoamperometry and the Cottrell analysis, we have determined that the PCET hopping kinetics within the pores of Ce-MOF-808 can change by orders of magnitude by altering the buffer species and the proton activity of the electrolyte. Furthermore, in all buffers, reductive reactions were ∼3−10 times faster in kinetics than the reverse oxidative reaction with the same electrochemical driving force, suggesting that the system, at first glance, violates the principle of microscopic reversibility. Isothermal titration calorimetry (ITC) and computational simulations corroborated that the buffer-node binding thermodynamics are quite distinct, depending on the chemical nature of the buffer and the oxidation state of the node. Together, these results suggest that the substrate and the product during the oxidative vs reductive reaction of Ce-MOF-808 are chemically different species, which explains the apparent ‘microscopic irreversibility.’ Thus, the rational modulation of electrolytes can dramatically enhance PCET kinetics, even though the solid electrodes remain identical. Implications of these findings are contrasted with the electrochemical/electrocatalytic behavior of other redox-active MOFs, heterogeneous catalysts, and enzymatic systems at the solid−liquid interface.

Ce-based MOF↗

DOE OSTI2026

Effect of Lithium Doping on MgO Hydroxylation and Carbonation

Recovery of magnesium from brines can potentially be used to source MgO (periclase) as a CO 2 sorbent or for Mg-based cements. However, it is not clear how common impurities in brines, such as lithium, affect the resulting MgO reactivity. Here, to test the effect of lithium incorporation on MgO reactivity for hydration and carbonation, we combined computational simulations with experiments. Experimentally altered (Mg,Li)O with a low dopant concentration (0.012 ± 0.002% w/w Li) was characterized using synchrotron-based X-ray scattering and high-resolution electron microscopy to measure reaction layer formation on (Mg,Li)O. Single-crystal X-ray diffraction analysis of (Mg,Li)O demonstrates that the incorporation of lithium leads to the formation of oxygen vacancies. The presence of vacancies is likely causing faster hydroxylation rates as predicted by ab initio molecular dynamics simulations. However, the faster hydroxylation rates likely lead to faster passivation of the surface because we observe thinner reaction layers on (Mg,Li)O samples both over short time periods (30 days) and over long time periods (28 years). After 28 years, the reaction layer on the (Mg,Li)O sample was less than one-third of the thickness of that of the pure MgO sample. In addition, over 30 days, reaction layers on (Mg,Li)O samples primarily formed at steps rather than on terraces, in contrast to our previous observations on MgO. Based on our results, naturally occurring impurities in MgO modify its reactivity even at very low concentrations and need to be considered for accurate reaction rate prediction for application of MgO as a CO 2 sorbent or in cements.

MgO↗

DOE OSTI2024

Electrochemically Determined and Structurally Justified Thermochemistry of H atom Transfer on Ti-Oxo Nodes of the Colloidal Metal–Organic Framework Ti-MIL-125

Titanium dioxide (TiO 2 ) has long been employed as (photo)electrodes for reactions relevant to energy storage and renewable energy synthesis. Proton-coupled electron transfer (PCET) reactions with equimolar amounts of protons and electrons at the TiO 2 surface or within the bulk structure lie at the center of these reactions. Because a proton and an electron are thermochemically equivalent to an H-atom, these reactions are essentially H-atom transfer reactions. Thermodynamics of H-atom transfer has a complex dependence on the synthetic protocol and chemical history of the electrode, the reaction medium, and many others; together, these complications preclude the understanding of the H-atom transfer thermochemistry with atomic-level structural knowledge. Herein, we report our success in employing open-circuit potential (E OCP ) measurements to quantitatively determine the H-atom transfer thermochemistry at structurally well-defined Ti-oxo clusters within a colloidally stabilized metal–organic framework (MOFs), Ti-MIL-125. The free energy to transfer H-atom, Ti 3+ O–H bond dissociation free energy (BDFE), was measured to be 68(2) kcal mol -1 . To the best of our understanding, this is the first report on using E OCP measurements to quantify thermochemistry on any MOFs. The proton topology, the structural change upon the redox reaction, and BDFE values were further quantitatively corroborated using computational simulations. Furthermore, comparisons of the E OCP -derived BDFEs of Ti-MIL-125 to similar parameters in the literature suggest that E OCP should be the preferred method for quantitatively accurate BDFE calculations. Here, the reported success in employing E OCP for nanosized Ti- MIL-125 should lay the ground for thermochemical measurements of other colloidal systems, which are otherwise challenging. Implications of these measurements on Ti-MIL-125 as an H-atom acceptor in chemical reactions and comparisons with other MOFs/metal oxides are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Interactions of Polar and Nonpolar Groups of Alcohols in Zeolite Pores

Understanding the quantitative interactions among zeolite pore walls, Bro̷nsted acid sites, and molecules with both polar and nonpolar regions is essential for scoping out the potential of zeolites as sorbents and catalysts. Purely siliceous zeolites (MFI and Beta in the present study) are hydrophobic, whereas those containing aluminum are considered hydrophilic, preferentially adsorbing organic molecules even in aqueous environments. To characterize these interactions, we use primary alcohols of increasing molecular weight, quantifying their specific interactions in the confined pore space of the alkyl (CH x ) and OH groups. Three types of interactions were identified: (i) alkyl CH x groups interacting with the zeolite pore walls (approximately 10 kJ mol −1 per carbon), (ii) alcohol OH groups interacting with the pore walls (30−35 kJ mol −1 ), and (iii) alcohol OH groups interacting with Bro̷nsted acid sites (37 kJ mol −1 ). All three interactions were well mirrored by computational simulations. The contribution of the alkyl CH x groups was inferred from the incremental increase in sorption enthalpy with increasing molecular weight; the interaction strength of the OH groups was determined by extrapolating the global adsorption enthalpy of the alcohols to a hypothetical OH group without an alkyl group. This value was identical to the adsorption enthalpy of water. The experiments demonstrated that only water has an adsorption enthalpy on zeolite pore walls lower than its condensation enthalpy (30−35 kJ mol −1 vs 45 kJ mol −1 ), limiting the concentration of water that can be adsorbed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Size-Dependent Optical Band Gaps in Metal–Organic Framework Nanoparticles

Decades of research into size-dependent semiconductor optical gaps have focused on quantum confinement as the dominant mechanism. Emerging reports indicate that lattice strain─intentional or incidental─can impart optical shifts similar or greater in magnitude. Here, we report evidence of optical absorption and photoluminescence spectra of M(1,2,3-triazolate)2 (M = Mg, Cr, Mn, Fe, Co, Cu, Zn, or Cd) nanoparticles that blueshift from bulk values with decreasing particle sizes in a manner that defies explanation by conventional quantum confinement. Here, the phenomenon persists for particle sizes as large as 200 nm, whereas quantum confinement generally ceases beyond 20–30 nm diameters and follows a weaker dependence on the particle radius. Computational simulations and crystallographic analysis suggest that this behavior arises from size-dependent changes to metal–linker bonding that manifest in strain values comparable to literature reports of strain-induced optical shifts in other classes of materials. This behavior appears beyond this family of materials in other notable examples of metal–organic frameworks (MOFs), including the well-studied Cu3(trimesate)2 (CuBTC), where smaller sizes correlate with blueshifted optical gaps. Taken together, these results represent one of the few examples of size-dependent strain in crystalline materials and reinforce the emerging view that MOFs become softer materials when isolated as nanoparticles.

Electrical conductivity↗

DOE OSTI2024

Electrochemical formation of bis(fluorosulfonyl)imide-derived solid-electrolyte interphase at Li-metal potential

Lithium bis(fluorosulfonyl)imide-based liquid electrolytes are promising for realizing high Coulombic efficiency and long cycle life for next-generation Li-metal batteries. However, the role of anions in the formation of solid-electrolyte interphase remains unclear. Herein, we combine electrochemical analyses and X-ray photoelectron spectroscopy measurements both with and without sample washing, together with computational simulations, to propose the reaction pathways of electrolyte decomposition and correlate the interphase component solubility with the efficacy of passivation. Additionally, we discovered that not all the products derived from interphase-forming reactions were incorporated into the resulting passivation layer, with a notable portion present in the liquid electrolyte. Further, we found that the high-performance electrolytes can afford a sufficiently passivating interphase with minimized electrolyte decomposition, by incorporating more anion-decomposition products. Overall, this work presents a systematic approach of coupling electrochemical and surface analyses to paint a comprehensive picture of solid-electrolyte interphase formation, while identifying key attributes of high-performance electrolytes for guiding future designs.

Yu, Weilai [Stanford University, CA (United States↗

DOE OSTI2026

Lithium metal-mediated electrochemical reduction of per- and poly-fluoroalkyl substances

Per- and poly-fluoroalkyl substances (PFAS) have substantial environmental and health hazards. Unfortunately, current degradation routes require high temperatures or corrosive conditions and/or lead to incomplete defluorination and the generation of shorter alkyl chains. Inspired by the lithium-metal battery literature, here we develop an electrochemical degradation process that leverages reactive metals and highly reducing environments. Here, we show that electrodeposited lithium metal can enable 95% degradation and 94% defluorination of perfluorooctanoic acid to LiF without forming any shorter C 2 –C 6 PFAS as end products. Using computational simulations, we find that electron transfer from lithium to perfluorooctanoic acid leads to rapid C–F bond cleavage, fluoride formation and carbon chain fragmentation. We expand the scope to other PFAS compounds and demonstrate substantial degrees of degradation on over 22 different PFAS, plus complete mineralization to inorganic fluorides. Finally, we use the mineralized F − as a fluorine source for the synthesis of fluorinated non-PFAS compounds to complete a circular fluorine loop from waste to valuable product.

Sarkar, Bidushi [Univ. of Chicago, IL (United Stat↗

DOE OSTI2025

Space-time crystals from particle-like topological solitons

Time crystals are unexpected states of matter that spontaneously break time-translation symmetry either in a discrete or continuous manner. However, spatially mesoscale space-time crystals that break both space and time symmetries have not been reported. Here we report a continuous space-time crystal in a nematic liquid crystal driven by ambient-power, constant-intensity unstructured light. Our numerically constructed four-dimensional configurations exhibit good agreement with these experimental findings. Although meeting the established criteria to identify time-crystalline order, both experiments and computer simulations reveal a space-time crystallization phase formed by particle-like topological solitons. The robustness against temporal perturbations and spatiotemporal dislocations shows the stability and rigidity of the studied space-time crystals, which relates to their locally topological nature and many-body interactions between emergent spontaneously twisted, particle-like solitonic building blocks. Their potential technological utility includes optical devices, photonic space-time crystal generators, telecommunications and anti-counterfeiting designs, among others.

Liquid crystals↗

DOE OSTI2025

Broadband coherent XUV light from e − / e + microbunching in an intense laser pulse

Attosecond pulses of coherent extreme ultraviolet (XUV) light are instrumental for investigating subatomic dynamics and can be produced using a free-electron laser (FEL). It has been suggested that an optical FEL, which employs a laser pulse in place of a conventional undulator, could enable a dramatically more compact implementation of such a light source. Yet, the high electron density and subsequent high emittance implied by an optical FEL makes this concept challenging to realize with an electron beam. There has been impressive progress in recent years producing collimated dense and relativistic beams of electrons and positrons in the laboratory. As we demonstrate here, the inherent stability of a quasi-neutral electron-positron beam mitigates Coulomb expansion, and renders it a promising alternative source of coherent light. Specifically, we show via computer simulations that broadband coherent light in the XUV domain, which takes the form of 8-as pulses at 92-as intervals, can be generated by microbunching of relativistic electrons and positrons in a laser pulse. This process occurs over a sub-millimeter length scale, enabling the development of light sources which are orders-of-magnitude more compact than existing sources, with potential applications in physics, chemistry, biology, and industry.

47 OTHER INSTRUMENTATION↗

DOE OSTI2025

Elucidating the geometric and electronic structure of a fully sulfided analog of an Anderson polyoxomolybdate cluster

The catalytic activity of transition metal sulfide (TMS) clusters in small molecule activation, redox transformations, and charge transfer has inspired the design of novel TMS-based materials for energy-related catalysis and chemical applications. Polyoxometalates (POMs), known for their structural diversity, can in principle be transformed into TMS clusters; however, fully sulfided analogs are rarely isolated, likely due to the strong tendency of uncapped TMS clusters to agglomerate. Here, we report the geometric and electronic structure of a capping ligand-free fully sulfided analog of heptamolybdate Anderson POM [Mo VI 7 O 24 ] 6− , synthesized through the sulfidation of a nanoconfined POM secured within a porous Zr-metal organic framework (NU-1000). A combined computational and experimental analysis indicates that the sulfided counterpart of the Anderson POM is geometrically and electronically more sophisticated than the parent POM. Comparison of experimental pair distribution function (PDF) data with computational simulations confirms that, unlike the oxygen-only [Mo VI 7 O 24 ] 6− cluster, the [Mo IV 7 (μ 3 -S) 6 (μ 2 -SH) 6 (S 2 ) 6 ] 2− polythiometalate (PTM) exhibits diverse sulfur anions (S 2− , HS − , S 2 2− ). DFT calculations indicate that H 2 S acts as a reducing agent, and together with terminal disulfide (S 2 2− ) ligands in the PTM structure, facilitates the complete reduction of all seven Mo VI centers in the parent POM to Mo IV . These findings are supported by X-ray photoelectron spectroscopy (XPS), which confirms exclusive Mo IV , and elemental analysis, which shows quantitative sulfur incorporation. Difference envelope density (DED) mapping further reveals that the PTM clusters are spatially confined within the MOF pores, preventing agglomeration and preserving molecular integrity.

Rabbani, S. M. Gulam [The Ohio State University, C↗

DOE OSTI2024

Preface for the Sherwood Fusion Theory 2023 special collection

Pursuing controlled fusion is a monumental endeavor that promises to unlock a transformational source of clean energy. Here, the International Sherwood Fusion Theory Conference is a crucial platform for researchers to discuss the physics governing plasmas in fusion-relevant configurations. From exploring the fundamental properties of plasma to cutting-edge computer simulation techniques, this conference fosters innovation and collaboration, propelling us closer to realizing the dream of a fusion-powered future.

Beidler, Matthew T. [Oak Ridge National Laboratory↗

DOE OSTI2025

Laser wavelength dependence of particle acceleration mechanisms in high intensity laser–solid density plasma interactions

We investigate the generation of relativistic electrons and the subsequent ion acceleration due to target-normal sheath acceleration when ultra-intense (⁠ I > 10 18 W/cm 2 ⁠) short pulse (⁠ τ L < 10ps⁠) lasers are incident onto solid density targets as laser wavelength is varied. Scaling laws for the hot electron temperature, T hot ⁠, and the maximum ion energy, E max ⁠, are recast as a function of laser wavelength. These predictions are compared to results from particle-in-cell computer simulations in a variety of geometries, including cases where realistic plasma density profiles as determined by a radiation hydrodynamics code are used. It is found that the wavelength dependence observed in simulation is less pronounced than what is predicted from the well-established scaling laws. An assessment of how switching to longer laser wavelengths, specifically 2 μm Tm:YLF technology, would impact current high energy density science applications and diagnostics is made.

Electromagnetism↗

DOE OSTI2026

Algorithm to extract direction in 2D discrete distributions and a continuous Frobenius norm

In this study, we present a novel algorithm for determining directionality in 2D distributions of discrete data. We compare a reference dataset with a known direction to a measured dataset with an unknown direction by the Frobenius norm of the difference (FND) to find the unknown direction. To generalize this concept, we develop a continuous Frobenius norm of the difference (CFND) as a continuous analog of the FND and derive its analytical expression. By relating fitted and normalized 2D Gaussian distributions, we show that the CFND approximates the FND, and we validate this relationship with computer simulations. We find that a first-order approximation of the CFND between two similar Gaussian distributions takes the form of an absolute sine function, offering a simple analytical form with potential applications in specialized areas such as segmented inverse beta decay neutrino detectors, astronomy, machine learning, and more. Our methodology consists of modeling a 2D Gaussian distribution, binning the data into a histogram, and encoding it as a square matrix. Rotating this matrix around its geometric center and comparing it to a measured dataset using the FND gives us rotational data that we fit with an absolute sine function. The location of the minimum of this fit is the angle closest to the true angle of the direction in the measured dataset. We present the derivation and discuss initial applications of the CFND in our novel algorithm, demonstrating its success in approximating directionality in 2D distributions.

Physics↗

DOE OSTI2026

Novel approach to general curvilinear coordinates for plasma fluid applications

In general geometry, plasma fluid equations include nonlinear geometric sources associated with fictitious forces, which pose significant challenges to computer simulations. We reformulate the plasma fluid hierarchy to rigorously preserve geometry and conservation properties critical to numerical simulations, while concealing the geometric sources. In their discrete form, the reformulated models conserve mass, angular momentum, and energy naturally, by simple analogy with the continuum equations. These conservation properties have minimal requirements in discrete space, namely, the anti-symmetry of the first derivative and the orthogonality of the scalar and cross products. By decoupling magnetic geometry, coordinate systems, and numerical discretization, this enables maximum flexibility while preserving physics fidelity. As a testbed, we apply the novel representation to the resistive magnetohydrodynamic system, which involves a complete set of curvilinear operations. We verify the correctness of the approach using steady state liquid metal flows and the classic Orszag–Tang vortex.

Halpern, Federico D. [General Atomics, San Diego, ↗

DOE OSTI2024

Active doping controls the mode of failure in dense colloidal gels

Mechanical properties of disordered materials are governed by their underlying free energy landscape. In contrast to external fields, embedding a small fraction of active particles within a disordered material generates nonequilibrium internal fields, which can help to circumvent kinetic barriers and modulate the free energy landscape. In this work, we investigate through computer simulations how the activity of active particles alters the mechanical response of deeply annealed polydisperse colloidal gels. We show that the “swim force” generated by the embedded active particles is responsible for determining the mode of mechanical failure, i.e., brittle vs. ductile. We find, and theoretically justify, that at a critical swim force the mechanical properties of the gel decrease abruptly, signaling a change in the mode of mechanical failure. The weakening of the elastic modulus above the critical swim force results from the change in gel porosity and distribution of attractive forces among gel particles, while below the critical swim force, the ductility enhancement is caused by an increase of gel structural disorder. Above the critical swim force, the gel develops a pronounced heterogeneous structure characterized by multiple pore spaces, and the mechanical response is controlled by dynamical heterogeneities. We contrast these results with those of a simulated monodisperse gel that exhibits a nonmonotonic trend of ductility modulation with increasing swim force, revealing a complex interplay between the gel energy landscape and embedded activity.

Zhou, Tingtao (ORCID:000000021766719X)↗

DOE OSTI2025

Unraveling the structure of molten salts: a review of high-energy x-ray diffraction (HEXRD) studies

High energy x-ray diffraction (HEXRD) has emerged as a pivotal technique for examining the structural properties of molten salts, which play vital roles in various high-temperature applications, such as energy storage, nuclear reactors, and metallurgy. This review provides a comprehensive analysis of recent progress in HEXRD studies of various molten salt systems, including molten halide salts, their mixtures, and salts with molecular anions, focusing on the insights into their atomic-scale arrangements. Additionally, we discuss methodological advancements that have improved the resolution and accuracy of HEXRD measurements, as well as recent developments in computational simulations, such as machine learning interatomic potentials, for interpreting results. This review aims to serve as a valuable resource for researchers in this field, offering a detailed overview of the progress of HEXRD studies of molten salts and suggesting pathways for further exploration.

Blockmon, Avery L. [Argonne National Laboratory (A↗

DOE OSTI2026

Roadmap for warm dense matter physics

This roadmap presents the state-of-the-art, current challenges and near future developments anticipated in the thriving field of warm dense matter (WDM) physics. Originating from strongly coupled plasma physics, high pressure physics and high energy density science, the WDM physics community has recently taken a giant leap forward. This is due to spectacular developments in laser technology, diagnostic capabilities, and computer simulation techniques. Only in the last decade has it become possible to perform accurate enough simulations & experiments to truly verify theoretical results as well as to reliably design experiments based on predictions. Consequently, this roadmap discusses recent developments of and contemporary challenges for theoretical methods and experimental techniques needed to describe, create and diagnose WDM. A large part of this roadmap is dedicated to specific WDM systems and applications in astrophysics, inertial confinement fusion and novel material synthesis.

dense astrophysical objects↗