Search NASA⌕ Search

SEARCH · Search NASA

Results for “Matrix Theory”

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

Quantum Monte Carlo Approaches to Na Intercalation on Bilayer Graphene

We have performed Quantum Monte Carlo (QMC) simulations on Na-intercalated bilayer graphene to study the evolution of electronic and optical properties upon Na intercalation into hard carbon layers. The objective was to model the optimal configuration of Na intercalation into a hard carbon matrix containing graphene regions. Our study showed that Na intercalation can be energetically stabilized at large interlayer distances (over 6 Å) in both AA- and AB-stacked bilayer graphene. In the QMC results, we found a significant band gap opening at the equilibrium interlayer distance of Na-intercalated bilayer graphene, while corresponding density functional theory (DFT) results showed no gap. This difference between DFT and QMC results indicates that the gap opening induced by Na intercalation into a hard carbon is underestimated within the DFT framework. In addition, a zigzag configuration of Na atoms was found to be energetically stable at interlayer distances up to 10 Å, leading us to predict the existence of a local minimum of Na intercalation at large interlayer distance. These computation and modeling results can provide guidance on how to synthesize and optimize hard carbon with bilayer graphene regions that permit a zigzag intercalation configuration that will maximize and stabilize sodium hosting.

Binding energy↗

Systematic Study of the Self-Renormalized Nucleon Gluon PDF in Large-Momentum Effective Theory

We present a systematic study of the nucleon gluon parton distribution function (PDF) using the self-renormalized large-momentum effective theory (LaMET) approach in lattice QCD. This work extends previous gluon-PDF extractions by performing a detailed analysis of key systematic effects, including gauge-link smearing, lattice spacing, pion mass, and nucleon boost momentum. The self-renormalization framework mitigates ultraviolet divergences associated with Wilson-line self-energy and renormalon contributions by combining lattice matrix elements with perturbative short-distance information, thereby preserving the correct infrared structure. Calculations are performed on $N_f=2+1+1$ HISQ ensembles generated by the MILC Collaboration at three lattice spacings and two pion masses, with boosted nucleon states reaching momenta up to 2.2~GeV. We determine renormalization factors from zero-momentum matrix elements and apply hybrid renormalization to suppress discretization artifacts. After extrapolating large-separation behavior and performing Fourier transforms, we reconstruct quasi-PDFs and match them to lightcone PDFs using next-to-leading order Wilson coefficients. Our results demonstrate that smearing and lattice-spacing effects are under control, and pion-mass and lattice-spacing dependence is mild relative to the current $O(10^6)$ statistics; however, momentum dependence remains a significant source of uncertainty. Future work including even larger boost momenta will be essential to reduce systematics in lattice determinations of the gluon PDF and to advance toward precision QCD phenomenology at the LHC and the future Electron-Ion Collider.

FOS: Physical sciences↗

Praseodymium in the formal +5 oxidation state

Praseodymium in the +5 oxidation state is a long-sought connection between lanthanide, early-transition and actinide metal redox chemistries. Unique among the lanthanide series, evidence for molecular pentavalent praseodymium species has been observed in the gas phase and noble gas matrix isolation conditions. Here we report the low-temperature synthesis and characterization of a molecular praseodymium complex in the formal +5 oxidation state, [Pr 5+ (NP t Bu 3 ) 4 ][X − ] (where t Bu = tert -butyl and X − = tetrakis(pentafluorophenyl)borate or hexafluorophosphate). Single-crystal X-ray diffraction, solution-state spectroscopic, solution magnetometric, density functional theory and multireference wavefunction-based methods indicate a highly multiconfigurational singlet ground state. Finally, an inverted ligand field drives this unique electronic structure, which establishes a critical link in understanding the bonding of high-valent metal complexes across the periodic table.

chemical bonding↗

Separable physics-informed DeepONet: Breaking the curse of dimensionality in physics-informed machine learning

The deep operator network (DeepONet) has shown remarkable potential in solving partial differential equations (PDEs) by mapping between infinite-dimensional function spaces using labeled datasets. However, in scenarios lacking labeled data, the physics-informed DeepONet (PI-DeepONet) approach, which utilizes the residual loss of the governing PDE to optimize the network parameters, faces significant computational challenges, particularly due to the curse of dimensionality. This limitation has hindered its application to high-dimensional problems, making even standard 3D spatial with 1D temporal problems computationally prohibitive. Additionally, the computational requirement increases exponentially with the discretization density of the domain. Here, to address these challenges and enhance scalability for high-dimensional PDEs, we introduce the Separable physics-informed DeepONet (Sep-PI-DeepONet). This framework employs a factorization technique, utilizing sub-networks for individual one-dimensional coordinates, thereby reducing the number of forward passes and the size of the Jacobian matrix required for gradient computations. By incorporating forward-mode automatic differentiation (AD), we further optimize computational efficiency, achieving linear scaling of computational cost with discretization density and dimensionality, making our approach highly suitable for high-dimensional PDEs. We demonstrate the effectiveness of Sep-PI-DeepONet through three benchmark PDE models: the viscous Burgers’ equation, Biot’s consolidation theory, and a parameterized heat equation. Our framework maintains accuracy comparable to the conventional PI-DeepONet while reducing training time by two orders of magnitude. Notably, for the heat equation solved as a 4D problem, the conventional PI-DeepONet was computationally infeasible (estimated 289.35 h), while the Sep-PI-DeepONet completed training in just 2.5 h. These results underscore the potential of Sep-PI-DeepONet in efficiently solving complex, high-dimensional PDEs, marking a significant advancement in physics-informed machine learning.

Neural operator↗

Comparison between explicit and implicit discretization strategies for a dissipative thermal environment

We investigate strategies for simulating open quantum systems coupled to dissipative baths by comparing explicit wave function-based discretization [via multi-layer multi-configuration time-dependent Hartree (ML-MCTDH)] and the implicit density matrix-based master equation method [via tree tensor network hierarchical equations of motion (TTN-HEOM)]. For dissipative baths characterized by exponentially decaying bath correlation functions, the implicit discretization approach of HEOM—rooted in bath correlation function decompositions—proves significantly more efficient than explicit discretization of the bath into discrete harmonic modes. Explicit methods, like ML-MCTDH, require extensive mode discretization to approximate continuum baths, leading to computational bottlenecks. Case studies for two-level systems and a Fenna–Matthews–Olson complex model highlight TTN-HEOM’s superiority in capturing dissipative dynamics with relaxations with a minimal number of auxiliary modes, while the explicit methods are as exact as the HEOM in pure dephasing regimes. This comparison is enabled by the TENSO package, which has both ML-MCTDH and TTN-HEOM implemented using the same computational structure and propagation strategy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Affinity and Treatment of LDR Organics in Low Activity Waste by Cementitious Materials – 26109

The physiochemical interactions between organic contaminants, minerals, and polymer additives in grout/cementitious materials are being identified. This research supports a proposed Resource Conservation and Recovery Act (RCRA) treatment variance relevant to the Land Disposal Restriction (LDR) organics present in Hanford tank waste. It provides the technical basis to address the aspects documented by the Environmental Protection Agency in a draft memo necessary to establish that organics can be immobilized/stabilized in a cementitious matrix. The organic species may interact with cementitious minerals, including slag, fly ash, cement, and other components/dopants like activated carbon. A list of 132 reasonably expected LDR organics associated with Hanford tank waste (RPP-RPT-63493, Rev 1a) was screened by functional groups, octanol-water partitioning coefficients, and detection frequency in Hanford tank waste samples. A subset of 15 compounds spanning the identified properties was tested using sorption and leachate tests. These compounds were subjected to traditional batch sorption tests and leaching tests on/from Cast Stone cementitious material with and without activated carbon (a potential organic adsorption additive). The interactions between organic compounds were observed in traditional batch sorption and modified TCLP leachate tests for surrogate Cast Stone material with and without activated carbon addition. Specifically, phthalic acid (negatively charged) and 4-chloroaniline (polar) sorbed strongly to unmodified Cast Stone, while a group of three phenolic compounds (2,4,6-Trichlorophenol, Pentachlorophenol, o-Cresol) sorbed to and were retained by Cast Stone with a 1% activated carbon amendment. Mixtures of these compounds demonstrated visible retention of the organics within the matrix with a quantifiable affinity of the organics to the cement. This novel testing supports the notion that organics within Hanford tank waste would be retained onto the cement phase with the addition of simple materials easily added to the formulation as part of the treatment process. To further provide additional support to this claim, models employing Flory-Huggins theory in conjunction with macromolecular and organic molecule solubility parameters and species activity coefficients are being utilized to better understand these interactions. Within this context, Van Krevelen (VK) and Hansen organic and polymer solubility parameters were used to inform the Flory-Huggins interaction parameter, allowing for a quantifiable method of evaluating polymer-solvent compatibility and, by extension, interactions with proposed binding additives (e.g., activated carbon). These models and measurement should test the hypothesis for our systems understanding and applicability of Flory-Huggins theory, and whether the prediction of organics in the polymer phase versus the aqueous phase can be confidently made using this strategy. The estimation of the organic molecular activity coefficient in these high salt-saturated aqueous systems is another factor that will be assessed. This project aims to predict and understand complicated organic interactions in these heterogenous systems, guiding optimal treatment and stabilization methodologies.

Jolin, Will [Savannah River National Laboratory (S↗

Multiscale Modeling of Nanoparticle Precipitation in Oxide Dispersion-Strengthened Steels Produced by Laser Powder Bed Fusion

Laser Powder Bed Fusion (LPBF) enables the efficient production of near-net-shape oxide dispersion-strengthened (ODS) alloys, which possess superior mechanical properties due to oxide nanoparticles (e.g., yttrium oxide, Y-O, and yttrium-titanium oxide, Y-Ti-O) embedded in the alloy matrix. To better understand the precipitation mechanisms of the oxide nanoparticles and predict their size distribution under LPBF conditions, we developed an innovative physics-based multiscale modeling strategy that incorporates multiple computational approaches. These include a finite volume method model (Flow3D) to analyze the temperature field and cooling rate of the melt pool during the LPBF process, a density functional theory model to calculate the binding energy of Y-O particles and the temperature-dependent diffusivities of Y and O in molten 316L stainless steel (SS), and a cluster dynamics model to evaluate the kinetic evolution and size distribution of Y-O nanoparticles in as-fabricated 316L SS ODS alloys. The model-predicted particle sizes exhibit good agreement with experimental measurements across various LPBF process parameters, i.e., laser power (110–220 W) and scanning speed (150–900 mm/s), demonstrating the reliability and predictive power of the modeling approach. The multiscale approach can be used to guide the future design of experimental process parameters to control oxide nanoparticle characteristics in LPBF-manufactured ODS alloys. Additionally, our approach introduces a novel strategy for understanding and modeling the thermodynamics and kinetics of precipitation in high-temperature systems, particularly molten alloys.

Wang, Zhengming (ORCID:0000000241627112)↗

A flavor of SO(10) unification with a spinor Higgs

We investigate Higgs Parity unification — a realization of SO(10) grand unification based on the Higgs Parity mechanism in which the Standard Model (SM) Higgs resides in a spinor representation. The theory has an intermediate left-right symmetric stage where the SU(2)R symmetry breaking scale is fixed by the vanishing of the SM Higgs quartic coupling. The strong CP problem is solved by parity. Gauge coupling unification successfully predicts αs(MZ) to within 1%. The spinor Higgs naturally leads to a seesaw origin for SM flavor observables. We identify a novel mechanism where large mixing of third generation fermions with additional heavy vector-like fermions accounts for the anarchical nature of the PMNS matrix and the lack of hierarchy in the neutrino mass spectrum, relative to the up-quarks. A fit to quark and lepton masses and mixings, with a minimal parameter set, predicts 1) A testable relation between the top quark mass and αs(MZ) which is about (1 – 2)σ from current best fit values, 2) The order of magnitude of the baryon asymmetry of the universe, via leptogenesis from second-generation right-handed neutrino decays. 3) The proton decay and the neutron EDM are likely observable in next generation experiments, and 4) A normal ordered neutrino mass spectrum where 0νββ decay and the mass of the lightest neutrino are out of reach of next generation experiments.

Baryo-and Leptogenesis↗

Polarized 𝐽/𝜓 production in semi-inclusive DIS at large 𝑄 2 : Comparing quark fragmentation and photon-gluon fusion

We compare the relative importance of different mechanisms for polarized 𝐽/𝜓 production in semi-inclusive deep inelastic scattering processes at large 𝑄 2 . The transverse momentum dependent (TMD) factorization framework and nonrelativistic quantum chromodynamics are used to study the leading contributions from light quark fragmentation to polarized 𝐽/𝜓, and compared to direct production via photon-gluon fusion, which can proceed through color-singlet as well as color-octet mechanisms. We identify kinematic regimes where light quark fragmentation dominates, allowing for the extraction of the 3 𝑆$^{[8]}_{1}$ matrix element, as well as regimes where photon gluon fusion dominates, suggesting that the gluon TMD parton distribution function can be probed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Ab Initio Uncertainty Quantification of Neutrinoless Double-Beta Decay in 76 Ge

The observation of neutrinoless double-beta (0⁢𝜈⁢𝛽⁢𝛽) decay would offer proof of lepton number violation, demonstrating that neutrinos are Majorana particles, while also helping us understand why there is more matter than antimatter in the Universe. If the decay is driven by the exchange of the three known light neutrinos, a discovery would, in addition, link the observed decay rate to the neutrino mass scale through a theoretical quantity known as the nuclear matrix element (NME). Accurate values of the NMEs for all nuclei considered for use in 0⁢𝜈⁢𝛽⁢𝛽 experiments are therefore crucial for designing and interpreting those experiments. Here, we report the first comprehensive ab initio uncertainty quantification of the 0⁢𝜈⁢𝛽⁢𝛽-decay NME, in the key nucleus 76 Ge. Here, our method employs nuclear strong and weak interactions derived within chiral effective field theory and recently developed many-body emulators. Our result, with a conservative treatment of uncertainty, is an NME of 2.60$^{+1.28}_{−1.36}$, which, together with the best-existing half-life sensitivity and phase-space factor, sets an upper limit for effective neutrino mass of 187$^{+205}_{−62}$ meV. The result is important for designing next generation germanium detectors aiming to cover the entire inverted hierarchy region of neutrino masses.

Ab initio calculations↗

Improved Kelbg Potentials for Z > 1 and Application to Carbon Plasmas

In this work, we present a general form for the electron‐ion diffractive potential derived from the quantum pair density matrix and fit to the improved Kelbg potential for atomic numbers up to $Z = 54$. We apply classical molecular dynamics using the improved Kelbg potential for carbon with various forms of the Pauli potential to compute internal energies and pressures for hot, dense plasma conditions. Our results are compared to an equation of state model based on path integral Monte Carlo and density functional theory simulations to examine the extent to which the improved Kelbg potential reproduces the internal energy and pressure of carbon plasmas. The regions of validity for carbon agree generally with those derived previously for hydrogen once pressure ionization effects are incorporated. Based on our carbon results and previously published hydrogen studies, we discuss the general applicability and limitations of these potentials for equation of state studies in warm dense matter and high energy density plasmas.

general physics↗

DXRD : a user-friendly suite of two- and multiple-beam dynamical X-ray diffraction programs

The DXRD program suite consisting of a series of dynamical theory programs is introduced for computing dynamical X-ray diffraction from single crystals. Its interactive graphical user interfaces (GUIs) allow general users to make complicated calculations with minimal effort. It can calculate plane-wave Darwin curves of single crystals (or multiple crystals) for both the Bragg and Laue cases, including grazing-incidence diffraction and backward diffraction (with Bragg angles approaching 90°). It is also capable of simulating rocking curves for divergent incident X-ray beams with finite bandwidths. A unique feature of DXRD is that it provides a convenient GUI-based multiple-beam diffraction program that can accurately compute arbitrary N-beam diffraction of any geometry using a universal 4N × 4N matrix method. DXRD also provides a mapping program for plotting all the multiple-beam diffraction lines (monochromator glitches) in the azimuth–energy coordinate system. All these functions make DXRD a convenient and powerful software tool for designing crystal-based synchrotron/X-ray optics (monochromators, analyzers, polarizers, phase plates etc.) and for crystal characterization, X-ray spectroscopy and X-ray diffraction teaching.

Bragg reflection↗

Measurement of differential t-channel single top (anti)quark production cross-sections at 13 TeV with the ATLAS detector

The production of single top quarks and top antiquarks via the t-channel exchange of a virtual W boson is measured in proton-proton collisions at a centre-of-mass energy of 13 TeV at the Large Hadron Collider. The full Run 2 data sample recorded with the ATLAS detector in the years 2015–2018 is used, corresponding to an integrated luminosity of 140 fb−1. The absolute and normalised production cross-sections are measured differentially as a function of the transverse momentum and absolute rapidity of the top quark and top antiquark. In addition, the ratio of top quark to top antiquark production cross-sections is measured. The measured distributions are compared with next-to-leading-order quantum chromodynamics predictions obtained with different combinations of matrix-element generators, parton-shower programs and proton parton distribution functions, as well as to next-to-next-to-leading-order calculations. Overall, good agreement is observed between the measurements and the theoretical predictions. For most measured distributions, the sensitivity to differences between the predictions is limited by the systematic uncertainties in the measurement. The measured differential distributions are also interpreted in an effective field theory approach to constrain the Wilson-Coefficient CQq3,1$$ {C}_{Qq}^{3,1} $$ associated with a four-quark operator. The interpretation accounts for the effect of the selection efficiency, which is altered significantly by non-zero contributions from CQq3,1$$ {C}_{Qq}^{3,1} $$.

Aad, G↗

Accurate models of the added mass force of a uniform random distribution of spherical particles or bubbles

The added mass force resulting from the acceleration of a body in a fluid is of fundamental and practical interest in dispersed multiphase flows. Euler–Lagrange (EL) and Euler–Euler (EE) simulations require closure terms for the added mass force in order to accurately couple the conserved variables between phases. Presently, a more thorough understanding of the added mass force in a multi-particle system is developed based on potential flow resulting in a resistance matrix formulation analogous to Stokesian dynamics. This formulation is then used to generate a dataset of added mass resistance matrices for large systems of randomly generated particles. This methodology is used to create a volume fraction corrected binary model for predicting the added mass force in large systems as well as generate statistics of the added mass force in such systems. This work provides clarification to the theory of the added mass force for particle clouds, and modelling options that may be implemented in existing EL and EE codes.

42 ENGINEERING↗

Bayesian mixture model approach to quantifying the empirical nuclear saturation point

The equation of state (EOS) in the limit of infinite symmetric nuclear matter exhibits an equilibrium density, $n_0 \approx 0.16 \, \mathrm{fm}^{-3}$, at which the pressure vanishes and the energy per particle attains its minimum, $E_0 \approx -16 \, \mathrm{MeV}$. Although not directly measurable, the nuclear saturation point $(n_0,E_0)$ can be extrapolated by density functional theory (DFT), providing tight constraints for microscopic interactions derived from chiral effective field theory (EFT). However, when considering several DFT predictions for $(n_0,E_0)$ from Skyrme and Relativistic Mean Field (RMF) models together, a discrepancy between these model classes emerges at high confidence levels that each model prediction's uncertainty cannot explain. How can we leverage these DFT constraints to rigorously benchmark nuclear saturation properties of chiral interactions? To address this question, we present a Bayesian mixture model that combines multiple DFT predictions for $(n_0,E_0)$ using an efficient conjugate prior approach. The inferred posterior distribution for the saturation point's mean and covariance matrix follows a Normal-inverse-Wishart class, resulting in posterior predictives in the form of correlated, bivariate $t$-distributions. The DFT uncertainty reports are then used to mix these posteriors using an ordinary Monte Carlo approach. At the 95\% credibility level, we estimate $n_0 \approx 0.157 \pm 0.010 \, \mathrm{fm}^{-3}$ and $E_0 \approx -15.97 \pm 0.40 \, \mathrm{MeV}$ for the marginal (univariate) $t$-distributions. Combined with chiral EFT calculations of the pure neutron matter EOS, we obtain bivariate normal distributions for the nuclear symmetry energy and its slope parameter evaluated at $n_0$: $S_v \approx 32.0 \pm 1.1 \, \mathrm{MeV}$ and $L\approx 52.6\pm 8.1 \, \mathrm{MeV}$ (95\%), respectively. Furthermore, our Bayesian framework is publicly available, so practitioners can readily use and extend our results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Rapid Initial-State Preparation for the Quantum Simulation of Strongly Correlated Molecules

Studies on quantum algorithms for ground-state energy estimation often assume perfect ground-state preparation; however, in reality the initial state will have imperfect overlap with the true ground state. Here, we address that problem in two ways: by faster preparation of matrix-product-state (MPS) approximations and by more efficient filtering of the prepared state to find the ground-state energy. We show how to achieve unitary synthesis with a Toffoli complexity about 7 × lower than that in prior work and use that to derive a more efficient MPS-preparation method. For filtering, we present two different approaches: sampling and binary search. For both, we use the theory of window functions to avoid large phase errors and minimize the complexity. We find that the binary-search approach provides better scaling with the overlap at the cost of a larger constant factor, such that it will be preferred for overlaps less than about 0.003 . Finally, we estimate the total resources to perform ground-state energy estimation of Fe - S cluster systems, including the Fe Mo cofactor by estimating the overlap of different MPS initial states with potential ground states of the Fe Mo cofactor using an extrapolation procedure. With a modest MPS bond dimension of 4000 , our procedure produces an estimate of approximately 0.9 overlap squared with a candidate ground state of the Fe Mo cofactor, producing a total resource estimate of 7.3 × 10 10 Toffoli gates; neglecting the search over candidates and assuming the accuracy of the extrapolation, this validates prior estimates that have used perfect ground-state overlap. This presents an example of a practical path to prepare states of high overlap in a challenging-to-compute chemical system. Published by the American Physical Society 2025

Berry, Dominic W. (ORCID:0000000334461449)↗

2025 EIC-France Workshop: Physics Highlights and Perspectives

This document presents a synthesis of the theory contributions and discussions from the 2nd EIC-France Workshop, held at IJCLab (Orsay) on 1-3 December 2025. The workshop brought together members of the French hadron-physics community to review recent theoretical developments relevant to the future Electron-Ion Collider (EIC) and to coordinate national efforts in preparation for its early physics program. The report first summarizes the collider's initial running conditions and luminosity performance, as outlined in the EIC Early Science Matrix. It then provides concise overviews of the theoretical presentations on inclusive, semi-inclusive, exclusive, heavy-flavor, and small-x physics. Based on these discussions, two measurements emerged as especially well suited for early EIC operation and strongly aligned with areas of established French expertise: inclusive diffraction and inclusive quarkonium production. These channels offer clean signatures, robust theoretical interpretability, and direct sensitivity to fundamental QCD phenomena such as gluon saturation, heavy-quark dynamics, and the small-x structure of hadrons and nuclei. In addition, the workshop identified longer-term physics opportunities that will benefit from the full capabilities of the EIC after its ramp-up phase. These include accessing the three-dimensional structure of the pion through the Sullivan process and a broader program of exclusive three-body final states, both of which represent high-impact avenues for exploring hadronic structure and non-perturbative QCD. Together, the elements summarized in this report provide a coherent overview of the strategic priorities and scientific ambitions shaping the French community's contribution to the EIC physics program.

FOS: Physical sciences↗

Surrogate models for linear response

Linear response theory is a well-established method in physics and chemistry for exploring excitations of many-body systems. In particular, the quasiparticle random-phase approximation (QRPA) provides a powerful microscopic framework by building excitations on top of the mean-field vacuum; however, its high computational cost limits model calibration and uncertainty quantification studies. Here, we present two complementary QRPA surrogate models and apply them to study response functions of finite nuclei. One is a reduced-order model that exploits the underlying QRPA structure, while the other utilizes the recently developed parametric matrix model algorithm to construct a map between the system’s Hamiltonian and observables. Our benchmark applications, the calculation of the electric dipole polarizability of 180 Yb and the 𝛽-decay half-life of 80 Ni, show that both emulators can achieve 0.1%–1% accuracy while offering a 6–7 orders of magnitude speedup compared to state-of-the-art QRPA solvers. These results demonstrate that the developed QRPA emulators are well positioned to enable Bayesian calibration and large-scale studies of computationally expensive physics models describing the properties of many-body systems.

Beta decay↗