SEARCH · Engineering Papers

Results for “THREE-DIMENSIONAL CALCULATIONS”

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 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Charge Transport in Solvated Donor–Acceptor Functionalized Peptoids: Molecular Dynamics and Rate Theory

Scalable solar-energy conversion requires photoactive materials that combine the efficiency of natural photosynthetic systems with the stability and processability needed for practical applications. Achieving reliable charge transport in soft, self-assembled organic materials remains challenging, as structural fluctuations and environmental effects strongly influence charge-transfer (CT) rates. Here, we present a broadly applicable computational framework for evaluating CT rates in the condensed phase, combining Fermi’s golden rule rate theory with inputs from all-atom molecular dynamics (MD) simulations and first-principles electronic-structure calculations. The approach does not rely on system-specific parametrization and is applicable to a wide range of soft and disordered materials. We demonstrate the applicability and usefulness of the framework on redox-active peptoids functionalized with iron–porphyrin (Fe–P) complexes, a bioinspired platform with programmable donor–acceptor units and tunable three-dimensional organization. The calculated CT rates exhibit strong sensitivity to molecular conformation, with variations spanning several orders of magnitude. This dependence is shown to arise from the pronounced variation in diabatic electronic coupling with the relative orientations and separations of the Fe–P complexes across the conformational ensemble. The framework provides a consistent route for connecting atomistic structure to CT kinetics in the condensed phase and enables analysis of structure–rate relations in organic semiconducting systems.

Charge transfer↗

DOE OSTI2026

Modeling of injected current stream-induced 3D perturbations in local helicity injection plasmas

Solenoid-free tokamak startup techniques are essential for spherical tokamaks and offer a pathway to cost reduction and design simplification in fusion energy systems. Local helicity injection (LHI) is one such approach, employing compact edge current sources to drive open field line current that initiates and sustains tokamak plasmas. The recently commissioned Pegasus-III Experiment provides a platform for advancing this and other solenoid-free startup methods. This study investigates the effect of LHI on magnetic topology in Pegasus-III plasmas. A helical filament model represents the injected current, and the linear plasma response to its three-dimensional field is calculated with M3D-C1. Poincaré mapping reveals substantial flux surface degradation in all modeled cases. The onset of overlapping magnetic structures and large-scale surface deformation begins near Ψ N ≈ 0.37, indicating a broad region of perturbed topology extending toward the edge. In rotating plasmas, both single-fluid and two-fluid models exhibit partial screening of the n = 1 perturbation, with two-fluid calculations showing stronger suppression near the edge. In contrast, the absence of rotation leads to strong resonant field amplification in the single-fluid case, while the two-fluid case with zero electron rotation mitigates this amplification and preserves edge screening. Magnetic probe measurements indicate that modeling the stream with spatial spreading—representing distributed current and/or oscillatory motion—better reproduces measured magnetic power profiles than a rigid filament model. The results underscore the role of rotation and two-fluid physics in screening stream perturbations and point to plasma flow measurements and refined stream models as key steps toward improving predictive fidelity.

Schaefer, Carolyn E. [Univ. of Wisconsin, Madison,↗

DOE OSTI2026

Modeling The Nucleosynthetic Imprint of Stellar Merger Phenomena - Final Technical Report

This DOE CAREER project developed a comprehensive, multi-physics framework for modeling stellar mergers and their observable consequences. The work successfully integrated analytical models, three-dimensional hydrodynamic simulations, stellar evolution calculations, nucleosynthesis, and radiation transport into a unified pipeline. Key scientific advances include demonstrating stellar mergers as a unifying explanation for systems such as Betelgeuse and R Coronae Borealis stars, quantifying merger-driven chemical signatures, and establishing the connection between merger physics and circumstellar environments. A major outcome of the project is the development of SuperLite, an open-source Monte Carlo radiation transport code that enables the generation of synthetic spectra for astrophysical transients. This work bridges dynamical, thermal, and radiative timescales and provides a powerful, DOE-relevant computational framework for interpreting observations of supernovae and related transient phenomena, while also contributing to workforce development through the training of postdoctoral researchers and graduate students.

Chatzopoulos, Emmanouil [Louisiana State Universit↗

DOE OSTI2026

Optical conductivity of layered topological semimetal TaNiTe 5

Here, we present an infrared spectroscopy study of the layered topological semimetal TaNiTe 5 , a material with a quasi-one-dimensional structure and strong in-plane anisotropy. Despite its structural features, infrared reflectivity and electronic transport measurements along the 𝑎 and 𝑐 crystallographic axes show metallic behavior without evidence of reduced dimensionality. Optical conductivity reveals an anisotropic but conventional metallic response with low scattering rates and a single sharp infrared-active phonon mode at 396cm −1 ⁢(49meV). Ab initio calculations closely match the experimental optical data and confirm a three-dimensional electronic structure. Our results demonstrate that TaNiTe 5 behaves as a three-dimensional anisotropic semimetal in its electronic and optical properties.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Quantum stress and torsion distributions in the deuteron

Stress distributions in the deuteron are related to form factors of the asymmetric energy-momentum tensor through three-dimensional Fourier transforms. There are eleven such form factors, which we calculate in an impulse approximation. We compare the obtained form factors to prior results for the six form factors that have been previously calculated. We then elaborate on the formalism for relating the form factors to internal distributions of mass, mass flux, momentum, stresses, and forces, and obtain results for all of these distributions. We obtain the principal stresses for the symmetric part of the stress tensor, and show that the antisymmetric part describes the reorientation of fermion spin by torsion stress when the nucleon moves between the S and D waves. Force distributions in the nucleons depend on the so-called nonconserved form factors through the Cauchy momentum equation, and are nonradial owing to the presence of tensor forces and spin-orbit coupling.

QCD phenomenology↗

DOE OSTI2026

Magnetization orientation-dependent Shubnikov-de Haas oscillations in ferromagnetic Weyl semimetal Co 3 Sn 2 S 2

We report a study of Shubnikov–de Haas oscillations in high-quality single crystals of ferromagnetic Weyl semimetal Co 3 Sn 2 S 2 . The Fermi surfaces resolved in our experiments are three-dimensional and reflect an underlying trigonal crystallographic symmetry. Combined with density functional calculations, we identify that multiple Fermi surfaces in the system—of both electron and hole nature—arise from the energy dispersion of the (spin-orbit gapped) mirror-protected nodal rings. We observe an evolution of the Fermi surfaces with in-plane magnetic fields, in contrast to field perpendicular to the kagome lattice planes, which has little effect. Viewed alongside the easy-axis anisotropy of the system, our observation reveals an evolution of the electronic structure of Co 3 Sn 2 S 2 —including the Weyl points—with the ferromagnetic moment orientation. Through the case study of Co 3 Sn 2 S 2 , our results provide concrete experimental evidence of an anisotropic interplay via spin-orbit coupling between the magnetic degrees of freedom and electronic band singularities, which has long been expected in semimetallic and metallic magnetic systems.

36 MATERIALS SCIENCE↗

DOE OSTI2026

ARC physics basis–magnetohydrodynamics

ARC is designed to produce 400 ⁢MW of net electricity and prove the commercial feasibility of a fusion power plant. In order to achieve this goal ARC has to operate with optimal core performance in a stationary scenario that minimises wear on the first wall and divertor. This requires avoiding or mitigating magnetohydrodynamic (MHD) instabilities which have the potential to not only degrade the plasma core but also lead to deleterious transient heat loads on plasma facing components. Therefore, this work aims at characterising the MHD stability of the high performance ARC scenario and inform the design of error field correction coils. Firstly, simulations of vertical displacement events show that an in-vessel coil is not needed and instead the poloidal shaping coils can be used to control vertical stability. These simulations also inform the demands on the corresponding coil power supplies. Stability analysis of the ideal kink mode with or without a conducting wall and kinetic effects suggests that the ARC baseline scenario operates deeply in the stable region. Using RDCON, tearing modes at the 𝑚/𝑛 =2/1 and 3/2 surfaces (with poloidal mode number 𝑚, and toroidal mode number 𝑛) are shown to be linearly stable, and including thermal transport effects in the rational surfaces lead to further stabilisation. However, other transient plasma instabilities can seed neoclassical tearing modes (NTMs). The marginally stable width of NTMs in ARC strongly depends on the internal inductance and can fall below 0.1% of the normalised poloidal flux. Furthermore, an empirical cross-machine model of the 𝑛 =1 error field leading to a disruption predicts a critical error field larger than SPARC but smaller than ITER. Three-dimensional coils can be designed with the Generalised Purturbed Equilbium Code based on a simple model that calculates the maximum correctable error field that is limited by the neoclassical toroidal viscosity torque. Broad scans of different coil geometries identify a set of 2 rows of off-midplane coils to be a suitable solution. It is also determined that such a set of three-dimensional coils is capable of correcting 𝑛 =2 error fields to some degree and creating strong enough 𝑛 =2 or 𝑛 =3 edge resonant perturbation fields for the suppression of edge-localised modes at reasonable coil currents. The final design of the first ARC will be further informed by results from SPARC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Direct Measurement of Diffusion Coefficients: Evidence for Diffusive Stochastic Heating in Collisionless Plasmas

Open questions in collisionless plasma dissipation can be addressed using space-based observations in different astrophysical environments, with implications for both astrophysical and laboratory plasma systems. We study a low-𝛽, highly imbalanced, sub-Alfvénic stream observed by Parker Solar Probe (PSP) to identify and distinguish between signatures of stochastic heating (SH) and resonant heating (RH) by parallel ion cyclotron waves (∥-ICWs). Prior work studying this stream [Trevor A. Bowen et al., Stochastic heating in the sub-Alfvénic solar wind, Phys. Rev. Lett. 135, 255201 (2025)] showed that the SH rate, accounting for intermittency, matched the amplitude of the local energy transfer (LET) rate, while the RH rate did not. This comparison relied on a number of assumptions regarding the nature of the diffusive process and the calculation of the LET rate. We introduce a novel technique of inverting the proton guiding center equation to empirically measure velocity-space diffusion coefficients using three-dimensional proton velocity distribution functions, from the ion electrostatic analyzer (the Solar Probe Analyzer for Ions) on PSP. Measured diffusion coefficients are used to determine phase-space heating rates, leading to a calculation of a fully kinetic heating rate independent of assumptions made in prior work. We show that scale-dependent analytic expressions for SH via noncoherent fluctuations match the empirical measurements from PSP data, provided that we account for intermittency in the heating calculation. In contrast, the derived heating rates for SH that accounts for the effects of the helicity barrier and heating rates for RH via ∥-ICWs do not peak in the same region of velocity space as the empirical measurements, nor do they reach the required magnitude. Our approach provides novel methodology to uniquely identify and constrain heating processes in collisionless plasmas and shows evidence of a Fokker-Planck-like diffusive process in the near-Sun solar wind.

Plasma kinetic theory↗

DOE OSTI2026

Observation of flow vector fluctuations in p–Pb collisions at $\sqrt{{s}_{NN}}=5.02$ TeV

Measurements of transverse momentum (pT) and pseudorapidity (η) dependent flow vector fluctuations in p–Pb collisions at $$\sqrt{{s}_{NN}}=5.02$$ TeV at the CERN Large Hadron Collider are presented. By studying long-range two-particle correlations with a template fit method, potential biases from non-flow effects such as jets and resonance decays are effectively suppressed. Significant pT- and η-dependent fluctuations of the second-harmonic flow vector are observed with more than 5σ confidence in p–Pb collisions, similar to the observations in Pb–Pb collisions. The influence of residual non-flow effects has been evaluated and cannot account for the observed fluctuations, thereby confirming the observation of flow vector fluctuations in small collision systems at the LHC. Comparisons to model calculations from 3DGlauber+MUSIC+UrQMD and the parton transport model from AMPT are also presented. The measurements provide constraints on the theoretical modelling of the three-dimensional initial geometry and its event-by-event fluctuations, offering critical insights into the origin of collective flow in small collision systems at the LHC.

Abdallah, D A H↗

DOE OSTI2026

Changing-look Active Galactic Nuclei Powered by Disk Tearing

Changing-look active galactic nuclei (CLAGN) feature order-of-magnitude variability in both the continuum and broad line luminosities on months-to-years-long timescales, and are currently unexplained. Simulations have demonstrated that rotating black holes (BHs) sometimes tear apart tilted accretion disks. These tearing events violently restructure the disk on timescales much shorter than a viscous timescale, hinting at a connection to CLAGN. Here, we show that disk tearing can power changing-look events. We report synthetic observations of an extremely high-resolution, three-dimensional, general-relativistic magnetohydrodynamic simulation of a geometrically thin, tilted accretion disk around a rapidly rotating, 10 8 M ⊙ BH. We perform ray-tracing calculations that follow the thermal disk emission to both a line-of-sight camera and to a distribution of cameras in a prescribed torus-like broad-line region (BLR). The continuum photoionizes the BLR, and we calculate the resulting spectrum. Both the continuum and line luminosities undergo order-of-magnitude swings on months-to-years-long timescales. We find shorter, week-long variability driven by the geometric precession of the inner disk and an intraday quasiperiodic oscillation driven by radial breathing of the inner disk. When the torn disk precesses, it causes asymmetric illumination of the BLR, driving time-evolving red-to-blue asymmetries of the broad emission lines that may be a smoking gun for disk tearing. We also make predictions for future photometric observations from ULTRASAT and Vera Rubin Observatory, given a Shakura–Sunyaev–like effective temperature profile, both of which may play an important role in detecting future changing-look events. Future work should consider nonthermal emission, especially from the time-variable corona.

accretion↗

DOE OSTI2026

Integrated machine learning-molecular dynamics framework for electrolyte property prediction

Electrochemical stability windows determine the operating range of battery electrolytes, yet accurate prediction remains challenging because stability emerges from statistical ensembles of local solvation environments rather than single ground-state molecular structures. Traditional density functional theory calculations on energy-minimized clusters cannot capture the thermal variations in local coordination environments and geometries that govern decomposition, while SMILES-based machine learning methods lack explicit representation of three-dimensional solvation structure and ion pairing. Here, we introduce a structure-aware machine learning framework that predicts frontier orbital energies (HOMO and LUMO) directly from molecular dynamics-sampled solvation configurations, achieving sub-0.6 eV accuracy at computational costs 3–4 orders of magnitude lower than first-principles methods. Across twelve representative battery electrolytes, we demonstrate that solvent-separated and contact ion pairs exhibit strong size- and local chemistry dependent electronic stability, with variations in coordination shifts of HOMO or LUMO level by 2–3 eV, and that extended solvation structure and partially desolvated environment further modulate stability by up to 3 eV. By encoding the statistical nature of electrochemical failure through ensemble sampling of explicit solvation geometries, our approach enables high-throughput screening and rational design of next-generation battery electrolytes with mechanistic understanding of structure–property relationships.

Energy - Storage↗

DOE OSTI2026

ML-Based Pebble Power Reconstruction for Pebble Bed Reactor Analysis

Pebble power reconstruction has been explored to complement the conventional homogenized modeling approach in pebble bed reactor (PBR) analysis, as detailed heterogeneous geometry calculations are computationally expensive. The random distribution of pebble fuels within the core challenges the application of conventional pin power reconstruction methods. To address this, we introduce a machine learning approach based on the transformer model, composed of encoder and decoder layers, to estimate the flux and power form functions for reconstructing individual pebble neutron fluxes and powers. The homogeneous neutron flux distribution within each spectral zone (SZ) is obtained from finite element solutions of global diffusion or transport calculations. Verification tests demonstrate that the trained transformer model accurately predicts power form functions over a range of conditions, including variations in pebble enrichment, location, type, SZ size, and burnup. In particular, verification using a three-dimensional PBR benchmark with burned pebbles shows good agreement in heterogeneous pebble power distributions between Griffin and Serpent. These results highlight the potential of applying conventional pin power reconstruction approaches to PBR cores with randomly distributed pebbles.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Short-range spin freezing in the double trillium lattice spin-liquid candidate KSrFe 2 ⁢(PO 4 ) 3 revealed via 31 P NMR

A comprehensive 31 P nuclear magnetic resonance (NMR) study, combined with thermodynamic measurements and first-principles band-structure calculations, has been conducted to explore the ground state of the 𝑆=5/2 double trillium lattice antiferromagnet KSrFe 2 ⁢(PO 4 ) 3 . Our experimental results indicate that the magnetic ground state is neither a conventional three-dimensional (3D) long-range order (LRO) nor a pure gapless spin-liquid state, as conjectured previously [Boya et al., APL Mater. 10, 101103 (2022)]. Specifically, the observation of a nearly field-independent NMR linewidth below 𝑇* = (3.5 ± 0.4)⁢ K, and a significant enhancement of spin-spin relaxation rate 1/𝑇 2 below 2⁢𝑇* (where 𝑇* is the characteristic temperature identified from the magnetic susceptibility), indicate a complex magnetic ground state where spin freezing coexists with persistent dynamics. Furthermore, we argue that the lack of magnetic LRO and the persistence of strong magnetic fluctuations in KSrFe 2 ⁢(PO 4 ) 3 are unlikely to originate from intersite K/Sr disorder, rather they arise due to intrinsic magnetic frustration. In conclusion, our findings position KSrFe 2 ⁢(PO 4 ) 3 into a broader family of geometrically frustrated magnets characterized by coexisting spin freezing and pronounced antiferromagnetic fluctuations, marking it as a promising platform for investigating exotic phenomena in 3D frustrated magnets.

Exotic phases of matter↗

DOE OSTI2026

Scalable Quantum Monte Carlo Method for Polariton Chemistry via Mixed Block Sparsity and Tensor Hypercontraction Method

We present a reduced-scaling auxiliary-field quantum Monte Carlo (AFQMC) framework designed for large molecular systems and ensembles, with or without coupling to optical cavities. Our approach leverages the natural block sparsity of the Cholesky decomposition (CD) of electron repulsion integrals in molecular ensembles and employs tensor hypercontraction (THC) to efficiently compress low-rank Cholesky blocks. By representing the Cholesky vectors in a mixed format, keeping high-rank blocks in block-sparse form and compressing low-rank blocks with THC, we reduce the scaling of exchange-energy evaluation from quartic to robust cubic in the number of molecular orbitals N, while lowering memory from cubic toward quadratic. Benchmark analyses on one-, two-, and three-dimensional molecular ensembles (up to ∼1,200 orbitals) show that (a) the number of nonzeros in Cholesky tensors grows linearly with system size across dimensions; (b) the average numerical rank increases sublinearly and does not saturate at these sizes; and (c) rank heterogeneity─some blocks nearly full rank and many low rank, naturally motivates the proposed mixed block sparsity and THC scheme for efficient calculation of exchange energy. In conclusion, we demonstrate that the mixed scheme yields cubic wall-time scaling with favorable prefactors and preserves AFQMC accuracy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Small Molecule Inhibitor-Modulated Al 2 O 3 Atomic Layer Deposition on Monolayer MoS 2 for Controlled Nucleation

The integration of ultrathin dielectrics on two-dimensional (2D) semiconductors is essential for advancing beyond-Si electronics. However, the intrinsic inertness of van der Waals 2D basal planes remains a primary bottleneck to achieving uniform dielectric nucleation and growth. Here, in this study, we introduce a small molecule inhibitor (SMI)-modulated thermal atomic layer deposition (ALD) strategy, exemplified by aluminum oxide (Al 2 O 3 ) ALD on monolayer molybdenum disulfide (1L MoS 2 ) with acetic acid (HAc) SMI. The ABC-type sequence comprises HAc inhibitor (A), trimethylaluminum (TMA) precursor (B), and deionized H 2 O coreactant (C). In situ quartz crystal microbalance (QCM) studies reveal robust HAc adsorption on Al 2 O 3 and suppression of subsequent oxide growth on HAc-passivated surfaces. When applied to 1L MoS 2 , this inhibitory pathway enables HAc to selectively passivate nascent Al 2 O 3 nuclei formed on the MoS 2 surface, limiting their three-dimensional (3D) island coarsening and redirecting precursor adsorption toward the uncovered basal plane. Consequently, nearly continuous ultrathin (∼1.5 nm) Al 2 O 3 films are achieved on 1L MoS 2 with markedly improved uniformity compared to standard Al 2 O 3 ALD using TMA and H 2 O, as validated by atomic force microscopy (AFM), cross-sectional scanning transmission electron microscopy (STEM), and energy-dispersive X-ray spectroscopy (EDS). Density functional theory (DFT) calculations further provide atomistic insight into HAc-modulated Al 2 O 3 nucleation, corroborating the energetic preference of HAc for Al 2 O 3 over MoS 2 and attenuated TMA adsorption on HAc-passivated surfaces. Spatially resolved Raman spectroscopy also confirms that the HAc-modulated process preserves the structural integrity of 1L MoS 2 , with only minimal strain and doping perturbations observed after dielectric deposition. This SMI-modulated approach offers a broadly applicable framework for controlling ALD nucleation across various inhibitors, ALD chemistries, and 2D materials, opening opportunities for reliable dielectric integration in next-generation nanoelectronics.

36 MATERIALS SCIENCE↗

DOE OSTI2026

4-Electron Oxygen Reduction Reaction (ORR) with Iron Phthalocyanine (FePc) Functionalized Nanowire Templated-3D Fuzzy Graphene (NT-3DFG)

Iron phthalocyanine (FePc) is a promising alternative to platinum-based catalysts for sustainable energy devices; however, the plane-symmetry of Fe-N 4 sites, random aggregation, and poor conductivity of FePc present major barriers for their application as oxygen reduction reaction (ORR) electrocatalysts. Here, we report the synergistic effects of FePc electrocatalysts supported by a nanowire-templated three-dimensional fuzzy graphene (FePc@NT-3DFG) substrate. The in situ functionalized oxygen groups (iFOGs) at the edge of NT-3DFG localize Fe active sites in FePc under alkaline ORR conditions. With a uniform FePc distribution through many single layers of graphene, the NT-3DFG substrates improve O 2 adsorption and catalytic activity while stabilizing the electrochemical activity during reactions. The FePc@NT-3DFG catalyst exhibits fast ORR kinetics with an extremely low Tafel slope of 28.3 ± 2.7 mV dec −1 , a higher half-wave potential of 0.911 ± 0.004 V (vs RHE), and notable long-term stability at 0.5 V (vs RHE) of 96.0 ± 0.4% retention after 30 h. Surface chemistry spectra validate electronic configuration modification of Fe at the iFOGs. Density functional theory calculations indicate that the extra layers of graphene improve oxygen adsorption. Moreover, additional exploration of other transition metal phthalocyanines supports the effects of iFOGs through the transition toward 4e − ORR. This work offers an expanded strategy for active site modification through edge-based graphene substrates for 4e − ORR.

X-ray absorption spectroscopy↗

DOE OSTI2026

Measurement of the longitudinal flow-plane decorrelation using multiplane cumulants in $\sqrt{s_{NN}}$ = 200 GeV Au+Au, Ru+Ru, and Zr+Zr collisions

Measurements of the variation of anisotropic flow-plane angles (𝛹 𝑛 ) with rapidity, commonly known as the flow-plane decorrelation, provide important insights into the initial conditions of the matter produced in heavy-ion collisions. Here, in this paper, using data collected by the STAR experiment, we report the first measurement of the four-plane correlator observable 𝑇 𝑛 ⁡{𝑏⁢𝑎;𝑑⁢𝑐}=⟨⟨sin⁡[𝑛⁢(𝛹$^𝑏_𝑛$ − 𝛹$^𝑎_𝑛$)]⁢sin⁡[𝑛⁢(𝛹$^𝑑_𝑛$ − 𝛹$^𝑐_𝑛$)]⟩⟩, where superscripts 𝑎, 𝑏, 𝑐, and 𝑑 denote sequential pseudorapidity (𝜂) regions with 𝑎 corresponding to the most backward region, 𝑏 and 𝑐 close to mid-rapidity with 𝜂 𝑏 < 0 and 𝜂 𝑐 > 0, and 𝑑 being the most forward. The measurement is performed for the elliptic and triangular flow (i.e.~𝑛 = 2 and 3) in Au+Au and isobar (Ru+Ru, Zr+Zr) collisions at $\sqrt{s_{NN}}$=200~GeV. The goal of calculating the correlation of the flow-plane angle variations from backward to mid-central, and from mid-central to forward regions, is to probe the systematic variation of flow angle over a wide 𝜂 range. In mid-central collisions (10−30% centrality), we find 𝑇 2 ⁡{𝑏⁢𝑎;𝑑⁢𝑐} =−0.004 ±0.001⁢(s⁢t⁢a⁢t) ±0.002⁢(s⁢y⁢s⁢t) independent of the collision system. Such a small value of 𝑇 2 favors a random-walk'' variation of the flow-plane angles, where the rapidity correlation length is smaller than the entire region under study. These measurements provide new information on the decorrelation patterns in the system and offer a quantitative estimate of possible systematic variations in anisotropic flow angles such astwist’’ between forward and backward regions. This opens new opportunities for understanding the three-dimensional structure and the time evolution of the quark-gluon plasma created in heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2026

Mesoscale atomic engineering in a crystal lattice

Controlling individual atoms using lasers, ion traps and scanning probe tips has transformed our understanding of matter and enabled breakthroughs in quantum science. Extending this control into three-dimensional (3D) solids and across mesoscopic scales, however, remains a foundational challenge. Electron irradiation in electron microscopes is known to induce atomic displacements, and atomic manipulation has been proposed and demonstrated. Yet repeated and deterministic control has remained elusive. Here, in this study, we demonstrate deterministic atomic engineering in a 3D crystal, creating ordered arrangements of more than 40,000 user-defined defects within minutes across a 150 nm × 100 nm × 13 nm volume. By steering individual Cr atoms in the magnetic semiconductor CrSBr into selected interstitial sites using an electron beam directed with sub-20-pm-scale accuracy, we create vacancy–interstitial complexes. The resulting impurity array forms a mesoscale crystal embedded within the host lattice, a new form of engineered artificial matter that remains stable at room temperature and outside the microscope. By tracking Cr atom displacements, we identify conditions under which the defect structures are predictable. Our calculations suggest that these defects form correlated impurity states with intra-defect optical transitions and inter-defect kinetic and Coulomb interactions. This establishes a generalizable platform for atomic defect engineering at mesoscopic, and potentially macroscopic, scales, opening opportunities for scalable quantum technologies, including deterministic colour-centre placement, quantum simulation of many-body lattice models and atomic-scale manufacturing.

74 ATOMIC AND MOLECULAR PHYSICS↗