SEARCH · Search NASA

Results for “Number 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 361 records · Page 20

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Fermi surfaces of the topological semimetal CaSn 3 probed through de Haas van Alphen oscillations

Abstract In the search of topological superconductors, nailing down the Fermiology of the normal state is as crucial a prerequisite as unraveling the superconducting pairing symmetry. In particular, the number of time-reversal-invariant momenta (TRIM) in the Brillouin zone enclosed by Fermi surfaces is closely linked to the topological class of time-reversal-invariant systems, and can experimentally be investigated. We report here a detailed study of de Haas van Alphen quantum oscillations in single crystals of the topological semimetal CaSn 3 with torque magnetometry in high magnetic fields up to 35 T. In conjunction with density functional theory based calculations, the observed quantum oscillations frequencies indicate that the Fermi surfaces of CaSn 3 enclose an odd number of TRIM, satisfying one of the proposed criteria to realize topological superconductivity. Nonzero Berry phases extracted from the magnetic oscillations also support the nontrivial topological nature of CaSn 3 .

Physics↗

DOE OSTI2026

Coupled Lindblad Pseudomode Theory for Simulating Open Quantum Systems

Coupled Lindblad pseudomode theory is a promising approach for simulating non-Markovian quantum dynamics on both classical and quantum platforms, with dynamics that can be realized as a quantum channel. We provide theoretical evidence that the number of coupled pseudomodes only needs to scale as polylog⁡(𝑇/𝜖) in the simulation time 𝑇 and precision 𝜖. Inspired by the realization problem in control theory, we also develop a robust numerical algorithm for constructing the coupled modes that avoid the nonconvex optimization required by existing approaches. We demonstrate the effectiveness of our method by computing population dynamics and absorption spectra for the spin-boson model. Furthermore, this Letter provides a significant theoretical and computational improvement to the coupled Lindblad framework, which impacts a broad range of applications from classical simulations of quantum impurity problems to quantum simulations on near-term quantum platforms.

Anderson impurity model↗

DOE OSTI2020

Density-Matrix Based Extended Lagrangian Born–Oppenheimer Molecular Dynamics

Extended Lagrangian Born–Oppenheimer molecular dynamics [ Phys. Rev. Lett. 2008, 100, 123004] is presented for Hartree–Fock theory, where the extended electronic degrees of freedom are represented by a density matrix, including fractional occupation numbers at elevated electronic temperatures. In contrast to regular direct Born–Oppenheimer molecular dynamics simulations, no iterative self-consistent field optimization is required prior to the force evaluations. To sample regions of the potential energy landscape where the gap is small or vanishing, which leads to particular convergence problems in regular direct Born–Oppenheimer molecular dynamics simulations, an adaptive integration scheme for the extended electronic degrees of freedom is presented. The integration scheme is based on a tunable, low-rank approximation of a fourth-order kernel, K, that determines the metric tensor, T ≡ K T K, used in the extended harmonic oscillator of the Lagrangian that generates the dynamics of the electronic degrees of freedom. Here, the formulation and algorithms provide a general guide to implement extended Lagrangian Born–Oppenheimer molecular dynamics for quantum chemistry, density functional theory, and semiempirical methods using a density matrix formalism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Empirical investigation of nuclear correlation function distributions in lattice QCD

Two-point correlation functions of systems with baryon number 𝐵 ∈ {1,2,3,4} are investigated using lattice quantum chromodynamics (QCD). In particular, the empirical distributions of importance-sampling Monte-Carlo samples of these correlation functions are examined as a function of the spacetime separation between the two points and the baryon number. While the exact forms of these distributions are not known for QCD, recent work has determined asymptotic expressions for analogous correlation function distributions in simpler theories such as scalar field theory and the disordered phase of the 𝑂⁡(𝑁) model. The theoretical 𝑂⁡(𝑁) model distributions are found to provide an accurate description of the empirical QCD distributions at zero momentum over a wide range of temporal separations for each baryon number when assessed with a range of different statistical tests. In particular, the behavior of the baryon number 𝐵 QCD correlation function at large temporal separation is well reproduced by the 𝑂⁡(𝑁 ∼ 2/𝐵) model distribution.

Lattice field theory↗

DOE OSTI2025

Graph reinforcement learning for exploring model spaces beyond the standard model

We present a methodology for performing scans of beyond the standard model (BSM) parameter spaces with reinforcement learning. We identify a novel procedure using graph neural networks that is capable of exploring spaces of models without the user specifying a fixed particle content, allowing broad classes of BSM models to be explored—in theory, the technique is applicable to nearly any model space with a prespecified gauge group. We provide a generic procedure by which a suitable graph grammar can be developed for any BSM model that features user-specified symmetry groups and a finite number of different possible particle species, the use of which is applicable to a variety of machine learning tasks over the actions of BSM theories beyond our particular reinforcement learning use case. As a proof of concept, we construct the graph grammar for theories with vectorlike leptons that may or may not be charged under a dark U ( 1 ) group, inspired by portal matter extensions of the sub-GeV vector portal/kinetic mixing simplified dark matter models. We then use this graph grammar to create a reinforcement learning environment tasked with creating models with these vectorlike leptons that are consistent with a list of a variety of precision observables. The reinforcement learning agent succeeds in developing models that can address the observed muon anomalous magnetic moment discrepancy while remaining consistent with flavor violation and electroweak precision observables, including both constructions that have previously been studied as well as new models that have not, to our knowledge, previously been identified. By inspecting the resulting ensembles of models that the agent produces and experimenting with different configurations for our reinforcement learning environment and graph grammar, we also infer various lessons about the development of these environments that can be transferable to reinforcement learning scans of more complicated model spaces and comment on future directions for the development of this technique into a more mature tool. Published by the American Physical Society 2025

Wojcik, George N.↗

DOE OSTI2020

Spectroscopic Probe Molecule Selection Using Quantum Theory, First-Principles Calculations, and Machine Learning

Probe molecule vibrational spectra have a long history of being used to characterize materials including metals, oxides, metal-organic frameworks, and even human proteins. Furthermore, recent advances in machine learning have enabled computationally generated spectra to aid in detailed characterization of complex surfaces with probe molecules. Despite widespread use of probe molecules, the science of probe molecule selection is underdeveloped. Here, we develop physical concepts, including orbital interaction energy and the energy overlap integral, to explain and predict the ability of probe molecules to discriminate structural descriptors. We resolve the crystal orbital overlap population (COOP) to specific molecular orbitals and quantify their bonding character, which directly influences vibrational frequencies. Using only a single adsorbate calculation from density function theory (DFT), we compute the interaction energy of individual adsorbate molecular orbitals with adsorption site atomic orbitals across many different sites. Combining the molecular orbital resolved COOP and changes in orbital interaction energy enables probe molecule selection for improved discrimination of various sites. We demonstrate these concepts by comparing the predicted effectiveness of carbon monoxide (CO), nitric oxide (NO), and ethylene (C 2 H 4 ) to probe Pt adsorption sites. Finally, using a previously developed machine learning framework, we show that models trained on hundreds of thousand C 2 H 4 spectra, computed from DFT, which regress surface binding-type and generalized coordination number (GCN), outperform those trained using CO and NO spectra. Lastly, a python package, pDOS_overlap, for implementing the electron density based analysis on any combination of adsorbates and materials, is also made available.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2023

Space-Dependent Calculation of the Multiplicity Moments for Shells With the Inclusion of Scattering

In recent work, we extended the methodology of multiplicity counting in nuclear safeguards by elaborating the one-speed stochastic transport theory of the calculation of the so-called multiplicity moments, i.e., the factorial moments of the number of neutrons emitted from a fissile item, following a source event from an internal neutron source [spontaneous fission and (α; n) reactions]. Calculations were made for solid spheres and cylinders, with the source being homogeneously distributed within the item. Recent measurements of the Rocky Flats Shells during the Measurement of Uranium Subcritical and Critical (MUSIC) campaign conducted by Los Alamos National Laboratory and assisted by the University of Michigan inspired us to extend the model to spherical shell geometry with a point source in the middle of the central cavity. Comparison of the calculated results with the experimental ones indicated that accounting for fission as the only neutron reaction (the standard procedure in the point model, adapted also in our work so far) was not sufficient for reaching good agreement with measurements. The model was therefore extended to include elastic scattering into the one speed formalism, whereas the effect of inelastic scattering was accounted for in an empirical way. After these extensions, good agreement was found between the calculated and the measured values. The paper describes the extension of the theory and provides concrete quantitative results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Revisiting the connection of baryon number, lepton number, and operator dimension

The effects of heavy new particles beyond the Standard Model can be conveniently captured through higher-dimensional effective operators. As noted long ago by Weinberg, the amount of baryon and lepton number an operator can carry is intricately connected to its mass dimension. We derive an improved inequality for this connection and compare it to explicit operator constructions up to mass dimension 25. For the effective field theory of Standard Model plus right-handed neutrinos, our relationship is even an equality up to high mass dimension.

Baryon number violation↗

DOE OSTI2023

General quantum algorithms for Hamiltonian simulation with applications to a non-Abelian lattice gauge theory

With a focus on universal quantum computing for quantum simulation, and through the example of lattice gauge theories, we introduce rather general quantum algorithms that can efficiently simulate certain classes of interactions consisting of correlated changes in multiple (bosonic and fermionic) quantum numbers with non-trivial functional coefficients. In particular, we analyze diagonalization of Hamiltonian terms using a singular-value decomposition technique, and discuss how the achieved diagonal unitaries in the digitized time-evolution operator can be implemented. The lattice gauge theory studied is the SU(2) gauge theory in 1+1 dimensions coupled to one flavor of staggered fermions, for which a complete quantum-resource analysis within different computational models is presented. The algorithms are shown to be applicable to higher-dimensional theories as well as to other Abelian and non-Abelian gauge theories. The example chosen further demonstrates the importance of adopting efficient theoretical formulations: it is shown that an explicitly gauge-invariant formulation using loop, string, and hadron degrees of freedom simplifies the algorithms and lowers the cost compared with the standard formulations based on angular-momentum as well as the Schwinger-boson degrees of freedom. The loop-string-hadron formulation further retains the non-Abelian gauge symmetry despite the inexactness of the digitized simulation, without the need for costly controlled operations. Such theoretical and algorithmic considerations are likely to be essential in quantumly simulating other complex theories of relevance to nature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Predictive Model of Charge Mobilities in Organic Semiconductor Small Molecules with Force-Matched Potentials

Charge mobility of crystalline organic semiconductors (OSC) is limited by local dynamic disorder. Recently, the charge mobility for several high mobility OSCs, including TIPS-pentacene, were accurately predicted from a density functional theory (DFT) simulation constrained by the crystal structure and the inelastic neutron scattering spectrum, which provide direct measures of the structure and the dynamic disorder in the length scale and energy range of interest. However, the computational expense required for calculating all of the atomic and molecular forces is prohibitive. In this work, we demonstrate the use of density functional tight binding (DFTB), a semiempirical quantum mechanical method that is 2 to 3 orders of magnitude more efficient than DFT. We show that force matching a many-body interaction potential to DFT derived forces yields highly accurate DFTB models capable of reproducing the low-frequency intricacies of experimental inelastic neutron scattering (INS) spectra and accurately predicting charge mobility. We subsequently predicted charge mobilities from our DFTB model of a number of previously unstudied structural analogues to TIPS-pentacene using dynamic disorder from DFTB and transient localization theory. The approach we establish here could provide a truly rapid simulation pathway for accurate materials properties prediction, in our vision applied to new OSCs with tailored properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Simulating a numerical UV completion of quartic Galileons

The Galileon theory is a prototypical effective field theory that incorporates the Vainshtein screening mechanism—a feature that arises in some extensions of general relativity, such as massive gravity. The Vainshtein effect requires that the theory contain higher order derivative interactions, which results in Galileons, and theories like them, failing to be technically well posed. While this is not a fundamental issue when the theory is correctly treated as an effective field theory, it nevertheless poses significant practical problems when numerically simulating this model. These problems can be tamed using a number of different approaches: introducing an active low-pass filter and/or constructing a UV completion at the level of the equations of motion, which controls the high momentum modes. These methods have been tested on cubic Galileon interactions, and have been shown to reproduce the correct low-energy behavior. Here we show how the numerical UV-completion method can be applied to quartic Galileon interactions, and present the first simulations of the quartic Galileon model using this technique. We demonstrate that our approach can probe physics in the regime of the effective field theory in which the quartic term dominates, while successfully reproducing the known results for cubic interactions. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

DOE OSTI2020

Probing Sulfur Chemical and Electronic Structure with Experimental Observation and Quantitative Theoretical Prediction of Kα and Valence-to-Core Kβ X-ray Emission Spectroscopy

An extensive experimental and theoretical study of the Kα and Kβ high-resolution X-ray emission spectroscopy (XES) of sulfur-bearing systems is presented here. This study encompasses a wide range of organic and inorganic compounds, including numerous experimental spectra from both prior published work and new measurements. Employing a linear-response time-dependent density functional theory (LR-TDDFT) approach, strong quantitative agreement is found in the calculation of energy shifts of the core-to-core Kα as well as the full range of spectral features in the valence-to-core Kβ spectrum. The ability to accurately calculate the sulfur Kα energy shift supports the use of sulfur Kα XES as a bulk-sensitive tool for assessing sulfur speciation. The fine structure of the sulfur Kβ spectrum, in conjunction with the theoretical results, is shown to be sensitive to the local electronic structure including effects of symmetry, ligand type and number, and, in the case of organosulfur compounds, to the nature of the bonded organic moiety. This agreement between theory and experiment, augmented by the potential for high-access XES measurements with the latest generation of laboratory-based spectrometers, demonstrates the possibility of broad analytical use of XES for sulfur and nearby third-row elements. The effective solution of the forward problem, i.e., successful prediction of detailed spectra from known molecular structure, also suggests future use of supervised machine learning approaches to experimental inference, as has seen recent interest for interpretation of X-ray absorption near-edge structure (XANES).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Circles and triangles, the NLSM and Tr(Φ 3 )

A surprising connection has recently been made between the amplitudes for Tr(Φ 3 ) theory and the non-linear sigma model (NLSM). A simple shift of kinematic variables naturally suggested by the associahedron/stringy representation of Tr(Φ 3 ) theory yields pion amplitudes at all loops. In this note we provide an elementary motivation and proof for this link going in the opposite direction, starting from the non-linear sigma model and discovering its formulation as a sum over triangulations of surfaces with simple numerator factors. This uses an ancient connection between “circles” and “triangles”, interpreting the equation y = $\sqrt{1 – x^2}$ both as parametrizing points a circle as well as generating the number of triangulations of polygons. A further simplification of the numerator factors exposes them as arising from the kinematically shifted Tr(Φ3) theory, and gives rise to novel tropical representations of NLSM amplitudes. The connection to Tr(Φ 3 ) theory defines a natural notion of “surface-soft limit” intrinsic to curves on surfaces. Remarkably, with this definition, the soft limit of pion amplitudes vanishes directly at the level of the integrand, via obvious pairwise cancellations. We also give simple, explicit expressions for the multi-soft factors for tree and loop-level integrands in the limit as any number of pions are taken “surface-soft”.

Scattering Amplitudes↗

DOE OSTI2022

Fast Particle-Wave Interactions and Alfvén Eigenmodes in JET Tokamak Plasmas

This document serves as the closeout report of DOE Grant Award No. DE-FG02-99ER54563 with project period 1 April 2015 through 31 March 2020. The project comprised the international collaboration between MIT and EU scientists on the JET facility to improve our understanding of the physics of energetic particle-wave interactions by measuring the damping rates of stable Alfvén Eigenmodes (AEs) and unstable energetic particle driven modes. More specifically, this project involved the continued participation of MIT, the Culham Centre for Fusion Energy, the Swiss Plasma Center, and theorists from various European laboratories and from UC Irvine. Past contributions from the University of São Paulo are gratefully acknowledged. The Alfvén Eigenmode Active Diagnostic (AEAD) on JET discharges was successfully upgraded in this period and has extracted physically useful information, in primarily deuterium plasma discharges. It is expected that such experiments would be continued during the DT campaign in CY 2021 to assess the damping rates of similar modes in the presence of alpha particles, a product of fusing burning plasma. Only JET would carry out such experiments in the near term in the world. It is important to note that this collaboration was continued under DOE Grant Award No. DE-SC0014264 from 1 April 2020 through the present. In this grant period, the AEAD was upgraded with a set of individual amplifiers for each of a set of six antennas in two toroidally opposite locations. These new amplifiers allowed targeted selection of the antennas’ toroidal spectrum for AEs of interest with toroidal mode number | n | ≤ 20. The resonant detection and measurement of the damping rates of AEs was obtained from magnetic probes that could be compared with theory and simulations. This is a key topic of investigation for ITER and all other next-step fusion experiments, where such modes will interact with energetic particles produced by fusion reactions (alpha particles), Neutral Beam Injection (NBI) and Ion Cyclotron Resonance Heating (ICRH). The majority of the observations during this period were of Toroidal Alfvén Eigenmodes (TAEs) as these are most commonly observed on JET; a number of measurements were made during dedicated TAE experiments. However, a new set of lower frequency band filters were procured with the goal of studying Geodesic Acoustic Modes (GAMs), Beta (Acoustic) AEs (BAE/BAAEs), and Reverse Shear AEs (RSAEs), also predicted by theory. Commissioning and optimizations were successfully completed in this period, and the diagnostic has been in operation during the more recent JET campaigns. Initial results have been obtained in dedicated TAE experiments and successfully compared to drift-kinetic theory. A wide range of theoretical studies have been undertaken in support of the upcoming JET campaigns experiments. To supplement the ongoing use of ideal MHD codes, such as MISHKA during the studies of TAEs, gyrokinetic simulations of low frequency AEs were performed in collaboration with UC Irvine. The Gyrokinetic Toroidal Code (GTC) was used to determine the structure, frequencies and stability of AEs in JET plasmas. Thus, a solid scientific foundation has been laid for future DT campaigns in JET in the 2021 operation period when the damping rates of relevant modes in the presence of alpha particles could be assessed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

All loop scattering for all multiplicity

We study the recently introduced curve integral formalism that defines a new family of formulas for the scattering amplitudes of the colored scalar trϕ 3 theory. We find that the curve integral manifests a very surprising fact about these amplitudes: the dependence on the number of particles, n, and the loop order, L, is effectively decoupled in these formulas. We derive the curve integrals at tree-level for all n. We then show that, for higher loop-order, it suffices to study the curve integrals for L-loop tadpole-like amplitudes, which have just one particle per color trace-factor. By combining these tadpole-like formulas with the tree-level results, we find formulas for the all n amplitudes at L loops. We illustrate this result by giving explicit curve integrals for all the amplitudes in the theory, including the non-planar amplitudes, through to two loops, for all n.

1/N Expansion↗

DOE OSTI2023

Flavor-changing light bosons with accidental longevity

We consider a model with a complex scalar field that couples to (e, μ) or (μ, τ) within the “longevity” window: [|m l 1 – m l 2 |, m l 1 + m l 2 ] in which l 1 and l 2 are the two different charged leptons. Within such a mass window, even a relatively large coupling (e.g. of the size commensurate with the current accuracy/discrepancy in the muon g – 2 experiment) leads to long lifetimes and macroscopic propagation distance between production and decay points. We propose to exploit several existing neutrino experiments and one future experiment to probe the parameter space of this model. For the μ – e sector, we exploit the muonium decay branching ratio and the production and decay sequence at the LSND experiment, excluding the parametric region suggested by g μ – 2 anomaly. For the τ – μ sector, we analyze three main production mechanisms of scalars at beam dump experiments: the Drell-Yan process, the heavy meson decay, and the muon scattering. We explore the constraints from the past CHARM and NuTeV experiments, and evaluate sensitivity for the proposed beam dump experiment, SHiP. The latter can thoroughly probe the parameter space relevant for the g μ – 2 anomaly.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Longitudinal short-distance constraints for the hadronic light-by-light contribution to (g - 2) μ with large-N c Regge models

While the low-energy part of the hadronic light-by-light (HLbL) tensor can be constrained from data using dispersion relations, for a full evaluation of its contribution to the anomalous magnetic moment of the muon (g – 2) μ also mixed- and high-energy regions need to be estimated. Both can be addressed within the operator product expansion (OPE), either for configurations where all photon virtualities become large or one of them remains finite. Imposing such short-distance constraints (SDCs) on the HLbL tensor is thus a major aspect of a model-independent approach towards HLbL scattering. Here, we focus on longitudinal SDCs, which concern the amplitudes containing the pseudoscalar-pole contributions from π o , η, η'. Since these conditions cannot be fulfilled by a finite number of pseudoscalar poles, we consider a tower of excited pseudoscalars, constraining their masses and transition form factors from Regge theory, the OPE, and phenomenology. Implementing a matching of the resulting expressions for the HLbL tensor onto the perturbative QCD quark loop, we are able to further constrain our calculation and significantly reduce its model dependence. We find that especially for the π0 the corresponding increase of the HLbL contribution is much smaller than previous prescriptions in the literature would imply. Overall, we estimate that longitudinal SDCs increase the HLbL contribution by Δa$^{LSDC}_{μ}$ = 13(6) × 10 -11 . This number does not include the contribution from the charm quark, for which we find a$^{c–quark}_{μ}$ = 3(1) × 10 –11 .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Null states and time evolution in a toy model of black hole dynamics

Abstract Spacetime wormholes can provide non-perturbative contributions to the gravitational path integral that make the actual number of statese S in a gravitational system much smaller than the number of states$$ {e}^{S_{\textrm{p}}} $$ e S p predicted by perturbative semiclassical effective field theory. The effects on the physics of the system are naturally profound in contexts in which the perturbative description actively involvesN=O(e S ) of the possible$$ {e}^{S_{\textrm{p}}} $$ e S p perturbative states; e.g., in late stages of black hole evaporation. Such contexts are typically associated with the existence of non-trivial quantum extremal surfaces. However, by forcing a simple topological gravity model to evolve in time, we find that such effects can also have large impact forN≪e S (in which case no quantum extremal surfaces can arise). In particular, even for smallN, the insertion of generic operators into the path integral can cause the non-perturbative time evolution to differ dramatically from perturbative expectations. On the other hand, this discrepancy is small for the special case where the inserted operators are non-trivial only in a subspace of dimensionD≪e S . We thus study this latter case in detail. We also discuss potential implications for more realistic gravitational systems.

Physics↗