Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “double shells”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Observation of Double Excitations in the Resonant Inelastic X-ray Scattering of Nitric Oxide

The nitrogen K-edge resonant inelastic X-ray scattering (RIXS) map of nitric oxide (NO) has been measured and simulated to provide a detailed analysis of the observed features. High-resolution experimental RIXS maps were collected using an in situ gas flow cell and a high-transmission soft X-ray spectrometer. Accurate descriptions of the ground, excited, and core-excited states are based upon restricted active space self-consistent-field calculations using second order multiconfigurational perturbation theory. The nitrogen K-edge RIXS map of NO shows a range of features that can be assigned to intermediate states arising from 1s → π* and 1s → Rydberg excitations; additional bands are attributed to doubly excited intermediate states comprising 1s → π* and π → π* excitations. These results provide a detailed picture of RIXS for an open-shell molecule and an extensive description of the core-excited electronic structure of NO, an important molecule in many chemical and biological processes.

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

Electronic structure analysis and DFT benchmarking of Rydberg-type alkali-metal-crown ether, -cryptand, and -adamanzane complexes

Density functional theory (DFT) and electron propagator theory (EPT) calculations were performed to study ground and excited electronic structures of alkali-metal (M) coordinated 9-crown-3, 24-crown-8, [2.1.1]cryptand, o-Me 2 -1.1.1, and 3 6 Adamanzane complexes. Each complex bears an expanded electron in the periphery and occupies diffuse 1p-, 1d-, 1f-type molecular orbitals (or superatomic 1P, 1D, 1F orbitals) in excited electronic states. The calculated superatomic shell model of the M(9-crown-3) 2 is 1S, 1P, 1D, 1F, 2S, 2P, 2D, 1G and it is held by all other complexes up to the studied 1F level. Due to the highly diffuse nature of the electron, the ionization energies of these complexes are significantly lower (1.6–2.0 eV) and hence these complexes belong to the superalkali category. Further, the ab initio EPT ionization energy and the excitation energies of the Li(9-crown-3) 2 were used to evaluate DFT errors associated with a series of exchange correlation functionals that span multiple rungs of Jacob's ladder (i.e., GGA, meta-GGA, global GGA hybrid, meta-GGA hybrid, range-separated hybrid, double-hybrid). Among these, the best performing functional is the range-separated hybrid CAM-B3LYP and the errors are within 6% of high-level ab initio EPT results. The accuracy of CAM-B3LYP is indeed transferable to similar complexes and hence the findings are expected to accelerate the progression of studies of Rydberg-type systems.

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

Tetragonal BaCoO 3 : A Co 4+ Ferromagnetic Mott Insulator

We report the stabilization of the metastable body-centered tetragonal (BCT) phase of BaCoO 3 (BCT-BaCoO 3 ) under high-pressure (15 GPa) and hightemperature (1200 °C) conditions using a mixture precursor. This double perovskite adopts the EuTiO 3 -type structure (space group I4/mcm), as confirmed by powder Xray diffraction and high-resolution STEM. X-ray photoelectron spectroscopy indicates a predominant Co 4+ oxidation state without detectable oxygen vacancies. Magnetization and heat capacity measurements reveal ferromagnetic ordering at T C ∼ 107 K, attributable to the BCT-BaCoO 3 phase. Above this temperature, the mixed-phase sample exhibits Curie−Weiss paramagnetism, a low-spin to high-spin crossover upon cooling, and a possible intermediate-spin state at elevated temperatures. Resistivity data indicates insulating behavior with weak magnetoresistance. DFT and DFT + DMFT calculations suggest that the insulating state originates from an orbitally selective transition sensitive to the nominal valence of the Co-d shell. The metastable BCT-BaCoO 3 phase cannot be retained in pure form at ambient pressure but can be stabilized by embedding it in a disordered mixture, offering a potential route to discover and preserve other high-pressure phases under ambient conditions.

Magnetic properties↗
DOE OSTI2022

β − decay of exotic P and S isotopes with neutron number near 28

β - decay of very neutron-rich isotopes of P and S, studied at the National Superconducting Cyclotron Laboratory (NSCL) using the Beta Counting Station (BCS) consisting of a Double Sided Strip Detector surrounded by clovers detectors for observing delayed γ transitions is reported here. β-decay half-lifes and delayed neutron emission probabilities were extracted for 42,43,44 P and 44,46 S by analyzing spatial and temporal correlations between implants and decay events in the Si detector with further coincidence with γ transitions. Detection of delayed γ rays allowed for the identification of negative parity 1p1h states in 42 S for the first time, also constraining the parent ( 42 P) spin/parity to 2 - or 3 - . For the most exotic isotope studied, 46 S, no strong γ transition was observed unlike lighter even-even S isotopes, thus implying the shift of Gamow Teller (GT) strength distribution to higher energies. In conclusion, comparison of experimental observations to detailed shell-model calculations using the SDPFSDG-MU interaction allowed us to infer the importance and role of First Forbidden (FF) β transitions as the neutron number approaches and then exceeds N = 28.

39 ≤ A ≤ 58↗
DOE OSTI2024

Low spin spectroscopy of neutron-rich Cl 43 , 44 , 45 via β - and β n decay

In this article, β - decay studies of neutron-rich isotopes 43,45 S performed at the National Superconducting Cyclotron Laboratory are reported. β -delayed $\gamma$ transitions were detected by an array of 16 clover detectors surrounding the Beta Counting Station, which consists of a 40×40 double-sided silicon strip detector followed by a single-sided silicon strip detector. β -decay half-lives were extracted for 43,45 S by correlating implants and decays in the pixelated implant detector and by considering further coincidences with $\gamma$ transitions in the daughter nuclei. Further, the level structures of 43,45 Cl are expanded by the addition of 20 and 8 new $\gamma$ transitions in 43 Cl and 45 Cl respectively, and core-excited negative-parity states were observed in both nuclei for the first time. For 45 S, a large fraction of the β -decay strength was observed feeding neutron-unbound states in 45 Cl, which, decaying by delayed neutrons, populated excited states in the β ⁢1⁢ n daughter, 44 Cl. Experimental observations were compared to detailed shell-model calculations using the SDPFSDG-MU interaction to highlight the role of the diminished N = 28 neutron shell gap and the near degeneracy of the proton s 1/2 and d 3/2 orbitals in the structure of the neutron-rich Cl isotopes. The current work also provides further support to a ground-state spin-parity assignment of 3/2 + in 45 Cl.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

A Theory of Ultrafast Charge Transfer Relaxation with Non-Innocent Solvent Molecules

In this work we revisit the photodynamics of tetracyanoethylene-hexamethylbenzene (TCNE- HMB), the molecular complex studied by Hochstrasser et al. [J. Chem. Phys. 100, 4797–4810, 1994] that has long challenged the applicability of Marcus’s theory of elec- tron transfer for predicting photochemical reactions. Using a novel black-box electronic structure algorithm (time-dependent density functional theory with one double, TD- DFT-1D) to efficiently run molecular dynamics that can treat charge recombination, we run ab initio surface hopping molecular dynamics and confirm that, for a polar solvent, charge recombination rates can be incredibly fast (indeed faster than the sol- vent relaxation time); for non-polar solvents, the rate is much slower. We demonstrate that, although Marcus theory cannot be directly applied, these nonequilibrium (and sometimes incredibly fast) photoexcited dynamics can be effectively explained within a two-state model without any evidence of a transition through a conical intersection. Most importantly, for this paradigmatic model system, we are able to identify two nuclear coordinates of interest (rather than the single coordinate predicted by Marcus or a full set of internal quantum modes studied by Bixon and Jortner): the solvent relaxation in the first shell (that strongly modulates the energies of the charge trans- fer state and differentiates time scales for relaxation) and a nuclear displacement in the TCNE-HMB complex arising from a handful of vibrations that induces non-Born Oppenheimer motion and eventually facilitates an abrupt electronic transition to the ground state. Altogether, these findings suggest a tractable generalization of Marcus theory for future simulations of photochemistry with non-innocent solvent environ- ments in the spirit of a Hamiltonian suggested by Stuchebrukhov (J. Chem. Phys. 107, 3821, 1997).

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

Magnetic ordering in terbium at high pressures and low temperatures

We have investigated magnetic ordering in the heavy lanthanide terbium (Tb) for pressures up to 70 GPa and temperatures down to 20 K by synchrotron x-ray diffraction at the Advanced Photon Source, Argonne National Laboratory. The x-ray diffraction studies were complemented by neutron diffraction studies at the Spallation Neutron Source, Oak Ridge National Laboratory. We focused on the higher-pressure phases of Tb, such as the double hexagonal close-packed (dhcp, above 16 GPa), distorted face-centered cubic (hR24, above 30 GPa), and orthorhombic (oF16, above 51 GPa) phases. We observed the phenomenon of spontaneous magnetostriction in the dhcp and hR24 phases, where magnetic ordering at low temperatures gives rise to subtle splitting in x-ray diffraction peaks. The neutron diffraction study of the dhcp phase at 20 GPa revealed magnetic peaks that can be assigned a propagation vector k = ($\frac{1}{2}$, 0, $\frac{1}{2}$), and a magnetic ordering temperature that is consistent with the magnetostriction effect observed with x-ray diffraction. The high-pressure and low-temperature x-ray diffraction study at 70 GPa in the collapsed oF16 phase did not show any evidence of magnetostriction, indicating a significant decrease in the magnetic moment of the 4f-shell due to band broadening at ultra-high compression.

36 MATERIALS SCIENCE↗
DOE OSTI2024

The effect of low temperature on poly(3-methyl- N -vinylcaprolactam)- b -poly( N -vinylpyrrolidone) diblock copolymer nanovesicles assembled from all-aqueous media

Nanosized polymeric vesicles (polymersomes) self-assembled from double hydrophilic copolymers of poly(3-methyl-N-vinylcaprolactam) n -b-poly(N-vinylpyrrolidone) m (PMVC n -b-PVPON m ) using all aqueous media are a promising platform for biomedical applications, because of their superior stability over liposomes in vivo and high loading capacity. Herein, we explored the temperature-sensitive behavior of PMVC 58 -b-PVPON 65 vesicles using transmission electron microscopy (TEM), dynamic light scattering (DLS), atomic force microscopy (AFM), and small-angle neutron scattering (SANS) in response to lowering the solution temperature from 37 to 25, 20, 14 and 4 °C. The copolymer vesicles with an average size of 350 nm at 37 °C were assembled from the diblock copolymer dissolved in aqueous solution at 4 °C. We show that while the polymersome's size gradually decreases upon the temperature decrease from 37 to 4 °C, the average shell thickness increases from 17 nm to 25 nm, respectively. SANS study revealed that the PMVC 58 -b-PVPON 65 vesicle undergoes a gradual structure evolution from a dense-shell vesicle at 37–25 °C to a highly-hydrated shell vesicle at 20–14 °C to molecular chain aggregates at 4 °C. From SANS contrast matching study, this vesicle behavior is found to be driven by the gradual rehydration of PMVC block at 37–14 °C. The shell hydration at 20–14 °C also correlated with the 4.4-fold decrease in the relative fluorescence intensity from vesicle-encapsulated fluorescent dye, indicating ~80% of the dye release within 12 hours after the vesicle exposure to 14 °C. No significant (<5%) dye release was observed for the vesicle solutions at 37–20 °C, indicating excellent cargo retention inside the vesicles. Our study provides new fundamental insights on temperature-sensitive polymer vesicles and demonstrates that the copolymer assembly into polymersomes can be achieved by decreasing a copolymer aqueous solution temperature below 14 °C followed by solution exposure to ≥20 °C. This type of all-aqueous assembly, instead of nanoprecipitation from organic solvents or solvent exchange, can be highly desirable for encapsulating a wide range of biological molecules, including proteins, peptides, and nucleic acids, into stable polymer vesicles without a need for organic solvents for dissolution of the copolymers that are amphiphilic at physiologically relevant temperatures of 20–37 °C.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Accurate prediction of core-level spectra of radicals at density functional theory cost via square gradient minimization and recoupling of mixed configurations

State-specific orbital optimized approaches are more accurate at predicting core-level spectra than traditional linear-response protocols, but their utility had been restricted due to the risk of "variational collapse" down to the ground state. We employ the recently developed square gradient minimization algorithm to reliably avoid variational collapse and study the effectiveness of orbital optimized density functional theory (DFT) at predicting second period element 1s core-level spectra of open-shell systems. Several density functionals (including SCAN, B3LYP, and ωB97X-D3) are found to predict excitation energies from the core to singly occupied levels with high accuracy (≤0.3 eV RMS error) against available experimental data. Higher excited states are, however, more challenging by virtue of being intrinsically multiconfigurational. We thus present a configuration interaction inspired route to self-consistently recouple single determinant mixed configurations obtained from DFT, in order to obtain approximate doublet states. This recoupling scheme is used to predict the C K-edge spectra of the allyl radical, the O K-edge spectra of CO + , and the N K-edge of NO 2 with high accuracy relative to experiment, indicating substantial promise in using this approach for the computation of core-level spectra for doublet species [vs more traditional time dependent DFT, equation of motion coupled cluster singles and doubles (EOM-CCSD), or using unrecoupled mixed configurations]. Finally, we also present general guidelines for computing core-excited states from orbital optimized DFT.

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

Determining g A / g V with High-Resolution Spectral Measurements Using a LiInSe 2 Bolometer

Neutrinoless double beta decay (0νββ) processes sample a wide range of intermediate forbidden nuclear transitions, which may be impacted by quenching of the axial vector coupling constant (g A /g V ), the uncertainty of which plays a pivotal role in determining the sensitivity reach of 0νββ experiments. Here we present measurements performed on a high-resolution LiInSe 2 bolometer in a “source = detector” configuration to measure the spectral shape of the fourfold forbidden β decay of 115 In . The value of g A /g V is determined by comparing the spectral shape of theoretical predictions to the experimental β spectrum taking into account various simulated background components as well as a variety of detector effects. We find evidence of quenching of g A /g V at >5 σ with a model-dependent quenching factor of 0.655 ± 0.002 as compared to the free-nucleon value for the interacting shell model. We also measured the 115 In half-life to be [5.18 ± 0.06(stat) $^{+0.005}_{- 0.015}$(sys)] x 10 14 yr within the interacting shell model framework. This Letter demonstrates the power of the bolometeric technique to perform precision nuclear physics single-β decay measurements, which along with improved nuclear modeling can help reduce the uncertainties in the calculation of several decay nuclear matrix elements including those used in 0νββ sensitivity calculations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Polarizabilities of neutral atoms and atomic ions with a noble gas electron configuration

Atomic polarizabilities play an important role in the development of force fields for molecular simulations, as well as for the development of qualitative concepts of atomic and molecular behavior. Coupled cluster theory at the coupled cluster singles doubles triples level with very large correlation-consistent basis sets with extended diffuse functions has been used to predict the polarizabilities of the atomic neutrals, mono-cations and mono-anions with a noble gas configuration. Additional corrections for scalar relativistic and spin–orbit effects were also included for the electron configurations of Kr, Xe, and Rn. The results are in excellent agreement with experiment or with other high level calculations where available. The current results for most of these species represent the best available values for the polarizabilities. The results show that the polarizability of H - is very difficult to calculate without extremely diffuse functions. The polarizability of H - is the largest value, 34.05 Å 3 , calculated for all species in the current study. The polarizabilities of the remaining halogen anions are also the best available values. The polarizabilities of the halogen anions (excluding F - ) and H - have a linear correlation with the electron affinity of the neutral atom. Spin– orbit effects, even for closed shell species, cannot be ignored for quantitative accuracy, and the inclusion of spin–orbit effects for Fr + , Rn, and At - increases the polarizability by 4%, 6%, and 15%, respectively.

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

The Landscape of Thermal Transients from Supernovae Interacting with a Circumstellar Medium

The interaction of supernova ejecta with a surrounding circumstellar medium (CSM) generates a strong shock, which can convert ejecta kinetic energy into observable radiation. Given the diversity of potential CSM structures (arising from diverse mass-loss processes such as late-stage stellar outbursts, binary interaction, and winds), the resulting transients can display a wide range of light-curve morphologies. We provide a framework for classifying the transients arising from interaction with a spherical CSM shell. The light curves are decomposed into five consecutive phases, starting from the onset of interaction and extending through shock breakout and subsequent shock cooling. The relative prominence of each phase in the light curve is determined by two dimensionless quantities representing the CSM-to-ejecta mass ratio η, and the breakout parameter ξ. These two parameters define four light-curve morphology classes, where each class is characterized by the location of the shock breakout and the degree of deceleration as the shock sweeps up the CSM. We compile analytic scaling relations connecting the luminosity and duration of each light-curve phase to the physical parameters. We then run a grid of radiation hydrodynamics simulations for a wide range of ejecta and CSM parameters to numerically explore the landscape of interaction light curves, and to calibrate and confirm the analytic scalings. We connect our theoretical framework to several case studies of observed transients, highlighting the relevance in explaining slow-rising and superluminous supernovae, fast blue optical transients, and double-peaked light curves.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The High Energy X-ray Probe (HEX-P): supernova remnants, pulsar wind nebulae, and nuclear astrophysics

HEX-P is a probe-class mission concept that will combine high spatial resolution X-ray imaging (< 10" full width at half maximum) and broad spectral coverage (0.2–80 keV) with an effective area far superior to current facilities (including XMM-Newton and NuSTAR) to enable revolutionary new insights into a variety of important astrophysical problems. HEX-P is ideally suited to address important problems in the physics and astrophysics of supernova remnants (SNRs) and pulsar wind nebulae (PWNe). For shell SNRs, HEX-P can greatly improve our understanding via more accurate spectral characterization and localization of non-thermal X-ray emission from both non-thermal-dominated SNRs and those containing both thermal and non-thermal components, and can discover previously unknown non-thermal components in SNRs. Multi-epoch HEX-P observations of several young SNRs (e.g., Cas A and Tycho) are expected to detect year-scale variabilities of X-ray filaments and knots, thus enabling us to determine fundamental parameters related to diffusive shock acceleration, such as local magnetic field strengths and maximum electron energies. For PWNe, HEX-P will provide spatially-resolved, broadband X-ray spectral data separately from their pulsar emission, allowing us to study how particle acceleration, cooling, and propagation operate in different evolution stages of PWNe. HEX-P is also poised to make unique and significant contributions to nuclear astrophysics of Galactic radioactive sources by improving detections of, or limits on, 44Ti in the youngest SNRs and by potentially discovering rare nuclear lines as evidence of double neutron star mergers. Throughout the paper, we present simulations of each class of objects, demonstrating the power of both the imaging and spectral capabilities of HEX-P to advance our knowledge of SNRs, PWNe, and nuclear astrophysics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

9 Li($\textit{d,p}$) reaction as a specific probe of 10 Li, the paradigm of parity-inverted nuclei around the $\textit{N}$ = 6 closed shell

In this work, we show, within the framework of renormalized nuclear field theory and of the induced reaction surrogate formalism, that the highly debated 10 Li structure, studied in a recent high statistics 9 Li ($\textit{d,p}$) 10 Li one-neutron transfer experiment, is consistent with, or better, requires, the presence of a virtual 1/2 + state of similar single-particle strength than that of the 1/2 – resonance at 0.45 ± 0.03 MeV. Based on continuum spectroscopy self-energy techniques, we find that the physical mechanism responsible for parity inversion in $^{10}_3$Li is the same as that at the basis of the similar phenomenon observed in $^{11}_4$Be and as that needed in 11 Li to have an important $\textit{s}$-wave ground-state component. In particular the strong dynamical coupling between the $s_{1/2}$ and the $d_{5/2}$ states, mediated by the quadrupole vibration of the core 9 Li. A phenomenon which also affects the strength distribution of the $d_{5/2}$ state, in particular, in the energy range of 3–4.5 MeV. Furthermore, this mechanism is also consistent with the (normal) sequence of the $1_{p_{1/2}}$ and $2_{s_{1/2}}$ levels in the $\textit{N}$ = 7 isotones $^{12}_5$B and $^{13}_6$C. The main aim of the present Rapid Communication is that of treating structure and reactions on equal footing and in a common language. In other words, the calculation of the 9 Li ($\textit{d,p}$) 10 Li reaction as a single conceptual step from individual single-particle motion and collective vibrations to absolute double differential cross sections of renormalized virtual and resonant final states, which can be directly compared with experiment.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Multipoint Observations of Dispersive Scale Alfvénic Field‐Line Resonances Associated With Substorm Auroral Beads

We present a case study of the field-aligned current (FAC) systems that transpire within the high-altitude auroral acceleration region of an “auroral bead” initiated double oval substorm observed on 23 February 2001 by the Cluster fleet. Conjunctive Cluster measurements and auroral images from IMAGE reveal that auroral bead current system formation and evolution is a multi-scale, injection-mediated process. The FACs at large scales vary on substorm evolution time scales (~minutes) in response to the injection and evolution of hotter denser magnetospheric plasma. Embedded within the large-scale FACs are intense short-scale (≲ few 10s of km) currents comprising dispersive scale Alfvén wave (DAW) fluctuations. The DAWs are a complex mixture of ingoing and reflected components that regularly interfere to form a broad spectrum of kinetic (dispersive) scale Alfvénic field-line resonances (KFLRs). The Alfvénic currents appear as a nested series of upward and downward FAC densities with amplitudes reaching a few 100 nA/m 2 . Energized field-aligned or counterstreaming electrons near keV energies and below are observed with parallel skews that vary in concert with variations in the DAW current sense. Positive correlations between DAW electric field energy densities and the energies of energized H + , He + , and O + outflow are observed, indicative of ion energization within the DAW fields. Due to their L-shell location (L~5.8–7.0) and associations with injections, the KFLRs are interpreted as the high-altitude auroral zone analog of KFLRs observed in the equatorial inner magnetosphere.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Compression of curium pyrrolidine-dithiocarbamate enhances covalency

Curium is unique in the actinide series because its half-filled 5f 7 shell has lower energy than other 5f n configurations, rendering it both redox-inactive and resistant to forming chemical bonds that engage the 5f shell. This is even more pronounced in gadolinium, curium’s lanthanide analogue, owing to the contraction of the 4f orbitals with respect to the 5f orbitals. However, at high pressures metallic curium undergoes a transition from localized to itinerant 5f electrons. This transition is accompanied by a crystal structure dictated by the magnetic interactions between curium atoms. Therefore, the question arises of whether the frontier metal orbitals in curium(III)–ligand interactions can also be modified by applying pressure, and thus be induced to form metal–ligand bonds with a degree of covalency. Here we report experimental and computational evidence for changes in the relative roles of the 5f/6d orbitals in curium–sulfur bonds in [Cm(pydtc) 4 ] – (pydtc, pyrrolidinedithiocarbamate) at high pressures (up to 11 gigapascals). We compare these results to the spectra of [Nd(pydtc) 4 ] – and of a Cm(III) mellitate that possesses only curium–oxygen bonds. Compared with the changes observed in the [Cm(pydtc) 4 ] – spectra, we observe smaller changes in the f–f transitions in the [Nd(pydtc) 4 ] – absorption spectrum and in the f–f emission spectrum of the Cm(III) mellitate upon pressurization, which are related to the smaller perturbation of the nature of their bonds. These results reveal that the metal orbital contributions to the curium–sulfur bonds are considerably enhanced at high pressures and that the 5f orbital involvement doubles between 0 and 11 gigapascal. Furthermore, our work implies that covalency in actinides is complex even when dealing with the same ion, but it could guide the selection of ligands to study the effect of pressure on actinide compounds.

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

matsim-agents v1.0

matsim-agents is a multi-agent AI framework for atomistic materials simulation and discovery. It orchestrates large language models (LLMs), machine-learned interatomic potentials (MLIPs), and DFT codes into a single agentic loop running on laptops and DOE leadership-class supercomputers. MULTI-AGENT ORCHESTRATION A LangGraph state machine with three nodes: a Planner that converts a natural-language research objective into structured tasks; an Executor that dispatches atomistic tools and loops until the queue is empty; and an Analyst that summarizes results into a human-readable report. State is checkpointed after every step and human-in-the-loop gates can be inserted at any edge. HYPOTHESIS-DRIVEN DISCOVERY CHAT An interactive REPL (matsim-agents chat) that couples LLM dialogue with atomistic simulation. Chemical formulas are automatically detected in conversation turns and trigger a full crystal-phase exploration: structure generation → relaxation → stability scoring → result injection back into the conversation, creating a closed hypothesis-refinement loop. CRYSTAL PHASE ENUMERATION Given a composition, the phase explorer enumerates prototypes by stoichiometry: elemental (fcc/bcc/hcp/sc/diamond), binary 1:1 (rocksalt/CsCl/zincblende/ wurtzite/fluorite/rutile), ternary 1:1:3 (cubic perovskite), ternary 1:2:4 (perovskite + spinel), quaternary 1:1:2:6 (Fm-3m double perovskite). 2-D prototypes (graphene, h-BN, MoS2 2H/1T) and multilayer stacking are also supported via --include-2d and --num-layers. SUPERCELL GENERATION AND SITE DECORATION Auto-tiling to a minimum atom count (--min-atoms), explicit NxNxN tiling (--supercell), symmetry-distinct site decorations (--n-orderings), and isotropic lattice-scale sweeps (--lattice-scales) for volume bracketing. MLFF RELAXATION AND STABILITY SCORING HydraGNN (multi-headed GNN) drives structure relaxation via ASE with FIRE, BFGS, or BFGSLineSearch. Stability output: delta-E/atom ranking across phases and a max-residual-force dynamical-stability proxy. Other MLIPs (MACE, NequIP, Orb) can be plugged in through the same interface. DFT BACKENDS Quantum ESPRESSO pw.x and VASP 6.6 are first-class labellers. Both have validated GPU builds and SLURM/PBS launchers for three DOE platforms: Frontier (AMD MI250X, ROCm), Aurora (Intel PVC, oneAPI), Perlmutter (NVIDIA A100, CUDA). QE produces ~100 binaries (pw.x, ph.x, epw.x, ...). VASP supports scf, relax, vc-relax, and vc-relax-shape run types. ACTIVE-LEARNING LOOP matsim-agents al run CONFIG.yaml drives an iterative HydraGNN-DFT loop: MD generates candidates → ensemble/MC-dropout uncertainty selects the most informative → DFT labels them in parallel inside one allocation → dataset grows → HydraGNN retrains → repeat. DFT backend is a single YAML toggle (dft.backend: vasp | qe). LLM-generated seed structures are supported (no curated POSCAR library needed). Config uses ${VAR}, ${VAR:-default}, ${VAR:?msg} shell-style substitution for cross-user/cross-site portability. LLM BACKENDS Ollama (local, default), vLLM (HPC multi-GPU serving), OpenAI, Anthropic, HuggingFace Transformers+Accelerate. Selected at runtime via flag or env var with no code changes. HPC PORTABILITY Same Python entry points run on Frontier (ROCm 7.2), Aurora (oneAPI), and Perlmutter (CUDA 12). DFT and ML stacks are never co-loaded in the same shell; they couple through the scheduler and filesystem. Advanced multi-node launchers (serve, discovery-chat, single-relaxation, active-learning, QE warm-start) are provided for all three platforms. CODABENCH COMPETITION BUNDLE A self-contained benchmark: 159 atomistic test structures across 11 material classes, 5 tasks (formation energy, forces, ML relaxation, AI-DFT relaxation, phase stability ranking), public/private leaderboard split (30/70), and four ready-to-run baselines: MACE-MP-0, HydraGNN, UMA, AllScAIP.

Lupo Pasini, Massimiliano [Oak Ridge National Labo↗
DOE OSTI2024

3D modeling of the Staged Z-pinch with the FLASH code

The Staged Z-pinch (SZP) fusion concept is a magneto-inertial compression scheme developed by Magneto-Inertial Fusion Technologies, Inc. (MIFTI), in which small amounts of fusion fuel are brought to fusion-relevant conditions by passing multi-million amperes strong current through a cylindrical shell of high atomic number material. One- and two-dimensional modeling performed by MIFTI with the MACH2 code, suggests that net fusion energy gain can be achieved when currents in the 10 million amperes range compress a 50%-50% mixture of deuterium and tritium gas. In this project we enhance and use the FLASH code to execute high-fidelity, simulations of various SZP configurations, a collaboration between the Flash Center for Computational Science at the University of Rochester and the Laboratory for Laser Energetics, and MIFTI. FLASH is a high-performance computing, multi-physics, radiation magnetohydrodynamic (MHD) code with extended physics capabilities, which is developed at the Flash Center. The goal of this project is to assess the shell/fuel stability of the pinch to two- and three-dimensional MHD instabilities and to utilize FLASH’s extended physics capabilities to understand how extended-MHD effects impact implosion dynamics and plasma conditions at stagnation. The project is a natural extension of an ongoing collaborative effort between MIFTI and the Flash Center through the U.S. DOE Advanced Research Projects Agency-Energy (ARPA-E) BETHE program and can provide MIFTI with simulation capabilities that are currently beyond their reach with MACH2. Ultimately, MIFTI wants FLASH to become one of their simulation workhorses for reliable, high-fidelity SZP platform design. During the project, the Flash Center team developed a suite of in one- and two-dimensional FLASH simulations to model different variants of the SZP platform, performed code-to-code comparisons with MACH2, and used the FLASH code to design and model SZP experiments like the Double Eagle experiments that MIFTI recently performed. Also, the Flash Center team developed the capabilities of FLASH to be able to do three-dimensional simulations of pulsed-power experiments with the code for the first time. Two publications from this effort are currently under review and two more are in preparation. MIFTI has also committed to the use of FLASH for future SZP simulation efforts.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗