Search NASA⌕ Search

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.

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

At least 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

SPARC-X: Quantum simulations at extreme scale - reactive dynamics from first principles

We have developed the massively parallel electronic structure code SPARC-X: a computational framework for performing Kohn-Sham Density Functional Theory (DFT) calculations that can scale linearly with the number of atoms in the system, while being able to leverage petascale and emerging exascale parallel computers to study chemical phenomena at unprecedented length and time scales. SPARC-X exploits a recent breakthrough in electronic structure methodologies: systematically improvable, strictly local, orthonormal, discontinuous real-space bases that efficiently and systematically capture the local chemistry of the system. With further adaptation using new machine-learning techniques and the use of the massively parallel Spectral Quadrature (SQ) electronic structure method, the algorithmic complexity and prefactor associated with DFT calculations involving semilocal as well as hybrid functionals are dramatically reduced. Using petascale computational resources, SPARC-X enables quantum mechanical simulations at length and time scales previously accessible only by empirical approaches, e.g., 1,000,000 atoms for a few picoseconds using semilocal functionals or 1,000 atoms for a few picoseconds using hybrid functionals. Using exascale resources, the sizes and times targeted are two orders of magnitude larger. Such a capability has applications in a wide variety of chemical sciences, including reactive interfaces where large length- and/or long time-scales are needed and traditional force fields fail. This is particularly important in dynamic catalysis, where bond breaking and formation must be understood in detail. We developed, tested, and employed the SPARC-X framework to understand the photocatalytic properties of TiO 2 nanoparticles, revealing finite size effects that cannot be captured with standard model systems or functionals. This integrated development and application strategy ensures that SPARC-X remains a robust, efficient, and scalable software package for quantum simulations on current petascale and emerging exascale computing resources.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Machine Learning Predictions of Transition Probabilities in Atomic Spectra

Forward modeling of optical spectra with absolute radiometric intensities requires knowledge of the individual transition probabilities for every transition in the spectrum. In many cases, these transition probabilities, or Einstein A-coefficients, quickly become practically impossible to obtain through either theoretical or experimental methods. Complicated electronic orbitals with higher order effects will reduce the accuracy of theoretical models. Experimental measurements can be prohibitively expensive and are rarely comprehensive due to physical constraints and sheer volume of required measurements. Due to these limitations, spectral predictions for many element transitions are not attainable. In this work, we investigate the efficacy of using machine learning models, specifically fully connected neural networks (FCNN), to predict Einstein A-coefficients using data from the NIST Atomic Spectra Database. For simple elements where closed form quantum calculations are possible, the data-driven modeling workflow performs well but can still have lower precision than theoretical calculations. For more complicated nuclei, deep learning emerged more comparable to theoretical predictions, such as Hartree–Fock. Unlike experiment or theory, the deep learning approach scales favorably with the number of transitions in a spectrum, especially if the transition probabilities are distributed across a wide range of values. It is also capable of being trained on both theoretical and experimental values simultaneously. In addition, the model performance improves when training on multiple elements prior to testing. The scalability of the machine learning approach makes it a potentially promising technique for estimating transition probabilities in previously inaccessible regions of the spectral and thermal domains on a significantly reduced timeline.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Search for baryon junctions in e+A collisions at the electron ion collider

Constituent quarks in a nucleon are the essential elements in the standard “quark model” associated with the electric charge, spin, mass, and baryon number of a nucleon. Quantum chromodynamics (QCD) describes nucleon as a composite object containing current quarks (valence quarks and sea (anti-)quarks) and gluons. These subatomic elements and their interactions are known to contribute in complex ways to the overall nucleon spin and mass. In the early development of QCD theory in the 1970s, an alternative hypothesis postulated that the baryon number might manifest itself through a non-perturbative configuration of gluon fields forming a Y-shaped topology known as the gluon junction. In this work, we propose to test such hypothesis by measuring (i) the Regge intercept of the net-baryon distributions for e+(p)Au collisions, (ii) baryon and charge transport in the isobaric ratio between e+Ru and e+Zr collisions, and (iii) target flavor dependence of proton and antiproton yields at large rapidity, transported from the hydrogen and deuterium targets in e+p(d) collisions. Our study indicates that these measurements at the EIC can help determine what carries the baryon number.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Amplitude/operator basis in chiral perturbation theory

We establish a systematic construction of the on-shell amplitude/operator basis for Chiral Perturbation Theory (ChPT) in D = 4 spacetime dimensions and with an arbitrary number of flavors N f . For kinematic factors, we employ spinor-helicity variables to construct the soft blocks, which are local amplitudes satisfying the Adler’s zero condition, as well as to take into account the reduction in the kinematic basis due to the Gram determinant, which arises at O(p 10 ) when the number of multiplicity N in an amplitude becomes large: N > D. For flavor factors, we include group-theoretic relations at small N f , N f ≤ N, which decreases the flavor basis. The result is obtained by adapting the Young tensor method of constructing the operator basis for generic effective field theories to the case of non-linearly realized symmetries. Working in the massless quark limit, we present purely mesonic operators for both even- and odd-parity at O(p 6 ) and O(p 8 ) for N = 6 and arbitrary N f , and establish a direct correspondence between the amplitude basis and the operator basis. Furthermore, the redundancy due to the Gram determinant is studied at O(p 10 ) for N = 6, 8 and 10.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Star Block Copolymers at Homopolymer Interfaces: Conformation and Compatibilization

Star block copolymers (s-BCPs), comprised of multiple linear diblock copolymers joined at a central point, are shown to segregate to the interface between two immiscible homopolymers that are identical to the blocks of the s-BCPs. The s-BCPs undergo a configurational transition at the interface, with different blocks of copolymers being embedded in their respective homopolymers, thereby bridging the interface and promoting adhesion. A series of 4-arm s-BCPs were synthesized with hydrogenated or deuterated polystyrene (PS/dPS) as the core block and poly(2-vinylpyridine) (P2VP) as the corona block, which was directly placed at the interface between the two homopolymers. Neutron reflectivity (NR) was used to determine the concentration profiles of the PS homopolymer, s-BCP core blocks, and P2VP total segments under equilibrium. The investigation varies the molecular weight (MW) and the total number of s-BCPs at the interface. Self-consistent-field theory (SCFT) was also employed to calculate the concentration profiles of the components at the interface, which were in excellent agreement with experimental results. The NR showed that the interfacial width between the homopolymers increased with the increasing number of s-BCPs at the interface up to a saturation limit. Beyond this limit, additional s-BCPs were released into the corona-miscible phase as unimolecular micelles. For a comparable interlayer thickness of s-BCPs at the interface, lower MW s-BCPs generated a broader interface. SCFT analysis suggested that, at the same packing density, the arms of the low MW s-BCPs align more parallel to the interface, while the arms of high MW s-BCPs adopt a more normal orientation, like their linear BCP counterparts. Furthermore, it was also observed that the core blocks, constrained by the junction points, were oriented more parallel and closer to the interface than the corona blocks. The phase behavior of the polymer blends revealed that s-BCP additives can efficiently reduce the domain size, with the low MW yielding smaller domain sizes due to the greater reduction in the interfacial energy and the high MW arresting phase separation due to their higher binding energy and a jamming of the interfacial assemblies. Asymmetric double cantilever beam (ADCB) tests demonstrated that s-BCPs promoted adhesion more efficiently than their linear BCP counterparts due to stronger binding energy per molecule, suggesting a more efficient compatibilizer for polymer upcycling. In conclusion, the results from these studies provide fundamental insights into the assembly of s-BCPs at homopolymer interfaces, the reduction of domain size, and promotion of adhesion, providing a strategy for the use of s-BCPs as stealth surfactants and universal compatibilizers.

interface↗
DOE OSTI2021

Moiré commensurability and the quantum anomalous Hall effect in twisted bilayer graphene on hexagonal boron nitride

The quantum anomalous Hall (QAH) effect is sometimes observed in twisted bilayer graphene (tBG) when it is nearly aligned with an encapsulating hexagonal boron nitride (hBN) layer. We propose that the appearance or absence of the QAH effect in individual devices could be related to commensurability between the graphene/graphene and graphene/hBN moiré patterns. We identify a series of points in the (θ GG , θ GBN ) twist-angle space at which the two moiré patterns are commensurate, allowing moiré band theory to be applied, and we show that the band Chern numbers are in this case sensitive to a rigid in-plane hBN displacement. Given this property, we argue that the QAH effect is likely only when (i) the (θ GG , θ GBN ) twist-angle-pair is close enough to a commensurate point that the two moiré patterns yield a supermoiré pattern with a sufficiently long length scale, and (ii) the supermoiré has a percolating topologically nontrivial QAH phase. For twist angles far from commensurability, the hBN layer acts as a source of disorder that can destroy the QAH effect. Our proposal can explain a number of current experimental observations. Finally, further experimental studies that can test this proposal more directly are suggested.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Automatic enhancement in 6D supergravity and F-theory models

We observe that in many F-theory models, tuning a specific gauge group G and matter content M under certain circumstances leads to an automatic enhancement to a larger gauge group G ' ⊂ G and matter content M ' ⊂ M . We propose that this is true for any theory G, M whenever there exists a containing theory G ' , M ' that cannot be Higgsed down to G, M . We give a number of examples including non-Higgsable gauge factors, nonabelian gauge factors, abelian gauge factors, and exotic matter. In each of these cases, tuning an F-theory model with the desired features produces either an enhancement or an inconsistency, often when the associated anomaly coefficient becomes too large. This principle applies to a variety of models in the apparent 6D supergravity swampland, including some of the simplest cases with U(1) and SU( N ) gauge groups and generic matter, as well as infinite families of U(1) models with higher charges presented in the prior literature, potentially ruling out all these apparent swampland theories.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Computation of Auger Electron Spectra in Organic Molecules with Multiconfiguration Pair-Density Functional Theory

Efficient and accurate computation of molecular Auger electron spectra for larger systems is limited by the rapid increase in the number of doubly ionized final states as the system size grows. Here, in this work, we benchmark the application of multiconfiguration pair-density functional theory with a restricted active space (RAS) reference wave function for computing the carbon K-edge decay spectra of 20 organic molecules. Decay rates are computed within the one-center approximation. We evaluate the performance of different basis sets and on-top functionals and find that multiconfiguration pair-density functional theory achieves accuracy comparable to RAS followed by second-order perturbation theory, but at significantly lower computational cost.

Fouda, Adam E. A. [Argonne National Laboratory (AN↗
DOE OSTI2022

Holographic Kolmogorov-Sinai entropy and the quantum Lyapunov spectrum

In classical chaotic systems the entropy, averaged over initial phase space distributions, follows a universal behavior. While approaching thermal equilibrium it passes through a stage where it grows linearly, while the growth rate, the Kolmogorov-Sinai entropy (rate), is given by the sum over all positive Lyapunov exponents. A natural question is whether a similar relation is valid for quantum systems. We argue that the Maldacena-Shenker-Stanford bound on quantum Lyapunov exponents implies that the upper bound on the growth rate of the entropy, averaged over states in Hilbert space that evolve towards a thermal state with temperature T, should be given by πT times the thermal state’s von Neumann entropy. Strongly coupled, large N theories with black hole duals should saturate the bound. To test this we study a large number of isotropization processes of random, spatially homogeneous, far from equilibrium initial states in large N, $\mathcal{N}$ = 4 Super Yang Mills theory at strong coupling and compute the ensemble averaged growth rate of the dual black hole’s apparent horizon area. We find both an analogous behavior as in classical chaotic systems and numerical evidence that the conjectured bound on averaged entropy growth is saturated granted that the Lyapunov exponents are degenerate and given by λ i = ±2πT. This fits to the behavior of classical systems with plus/minus symmetric Lyapunov spectra, a symmetry which implies the validity of Liouville’s theorem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Rotational excitations in rare-earth nuclei: A comparative study within three cranking models with different mean fields and treatments of pairing correlations

High-spin rotational bands in rare-earth Er (Z = 68), Tm (Z = 69), and Yb (Z = 70) isotopes are investigated by three different nuclear models. These are (i) the cranked relativistic Hartree-Bogoliubov (CRHB) approach with approximate particle number projection by means of the Lipkin-Nogami (LN) method, (ii) the cranking covariant density functional theory (CDFT) with pairing correlations treated by a shell-model-like approach (SLAP) or the so-called particle-number conserving (PNC) method, and (iii) cranked shell model (CSM) based on the Nilsson potential with pairing correlations treated by the PNC method. A detailed comparison between these three models in the description of the ground-state rotational bands of even-even Er and Yb isotopes is performed. The similarities and differences between these models in the description of the moments of inertia, the features of band crossings, and equilibrium deformations and pairing energies of the even-even nuclei under study are discussed. These quantities are considered as a function of rotational frequency and proton and neutron numbers. The changes in the properties of the first band crossings with increasing neutron number in this mass region are investigated. On average, a comparable accuracy of the description of available experimental data is achieved in these models. However, the differences between model predictions become larger above the first band crossings. Because of the time-consuming nature of numerical calculations in the CDFT-based models, a systematic study of the rotational properties of both ground-state and excited-state bands in odd-mass Tm nuclei is carried out only by the PNC-SCM. With few exceptions, the rotational properties of experimental 1-quasiparticle and 3-quasiparticle bands in 165,167,169,171 Tm are reproduced reasonably well. As a result, the appearance of backbendings or upbendings in these nuclei is well understood from the analysis of the variations of the occupation probabilities of the single-particle states and their contributions to total angular momentum alignment with rotational frequency.

150 ≤ A ≤ 189↗
DOE OSTI2021

Mass enhancement in 3 d and s - p perovskites from symmetry breaking

In some d-electron oxides, the measured effective mass m* exptl has long been known to be significantly larger than the model effective mass m* model deduced from mean-field band theory, i.e., m* exptl =βm* model , where β>1 is the “mass-enhancement” or “mass-renormalization” factor. Previous applications of density functional theory (DFT), based on a symmetry-restricted structure with the smallest number of possible magnetic, orbital, and structural degrees of freedom, missed such mass enhancement. This fact has been taken as evidence of strong electronic correlation, often described via the symmetry-restricted dynamic mean-field approach of the many-body theory, being the exclusive enabling physics. This paper uses instead a static density functional approach that does not restrict positional or spin degrees of freedom (symmetry-broken structures). This approach analyzes the contributions of different symmetry-broken modalities to mass enhancement for a few nominally highly correlated d-electron perovskites as well as the nominally uncorrelated, closed-shell s-p bonding perovskites. It shows that the energy-lowering symmetry-broken spin effects (e.g., nonzero local moment in the paramagnetic phase) and structural effects (e.g., atomic displacement) as described in mean-field DFT already manifest mass enhancement for both electrons and holes in a range of d-electron perovskites SrVO 3 , SrTiO 3 , BaTiO 3 , and LaMnO 3 , as well as p-electron perovskites CsPbI 3 and SrBiO 3 , including both metals (SrVO 3 ) and insulators (the rest). This is revealed only when enlarged unit cells of the same parent global symmetry, which are large enough to allow for symmetry-breaking distortions and concomitant variations in spin order, are explored for their ability to lower the total energy. Positional symmetry breaking that leads to mass enhancement includes octahedral rotation in halide perovskites such as CsPbI 3 , Jahn-Teller-like Q${+}\atop{2}$ distortion in LaMnO 3 , and bond disproportionation in SrBiO 3 , while magnetic symmetry breaking resulting in mass enhancement includes the formation of a distribution of local moments in SrVO3 that averages to zero in the paramagnetic phase. Not all symmetry breaking leads to significant mass enhancement, e.g., the rather small octahedral rotations in the nearly perfectly cubic SrTiO 3 cause negligible mass enhancement, as do the paraelectric displacements in BaTiO 3 . In principle, by ergodicity, the two descriptions, i.e., the symmetry-restricted dynamic approach with a single time-fluctuating magnetic moment and the symmetry-broken mean-field approach with a static distribution of spatially fluctuated local moments, are not mutually exclusive but are a choice of representation and consequently, a choice of computational efficiency. In approximate implementations, the symmetry-broken mean-field approach appears to remove much of what was strong correlation in dynamically correlated symmetry-restricted solutions, leaving smaller (“weak”) residual correlation with respect to the exact solution.

36 MATERIALS SCIENCE↗
DOE OSTI2021

General Quantum Field Theory of Flavor Mixing and Oscillations

We review the canonical transformation in quantum physics known as the Bogoliubov transformation and present its application to the general theory of quantum field mixing and oscillations with an arbitrary number of mixed particles with either boson or fermion statistics. The mixing relations for quantum states are derived directly from the definition of mixing for quantum fields and the unitary inequivalence of the Fock space of energy and flavor eigenstates is shown by a straightforward algebraic method. The time dynamics of the interacting fields is then explicitly solved and the flavor oscillation formulas are derived in a unified general formulation with emphasis on antiparticle content and effect introduced by nontrivial flavor vacuum.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Field moment expansion method for interacting bosonic systems

In this work, we introduce a numerical method and python package, https://github.com/andillio/CHiMES, that simulates quantum systems initially well approximated by mean field theory using a second order extension of the classical field approach. We call this the field moment expansion method. In this way, we can accurately approximate the evolution of first and second field moments beyond where the mean field theory breaks down. This allows us to estimate the quantum breaktime of a classical approximation without any calculations external to the theory. We investigate the accuracy of the field moment expansion using a number of well studied quantum test problems. Interacting Bosonic systems similar to scalar field dark matter are chosen as test problems. We find that successful application of this method depends on two conditions: the quantum system must initially be well described by the classical theory, and that the growth of the higher order moments be hierarchical.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Band-limited photodetection of temporal coherence

The quantum theory of optical coherence plays a ubiquitous role in identifying optical emitters. An unequivocal identification, however, presumes that the photon number statistics is resolved from timing uncertainties. We demonstrate from first principle that the observed nth-order temporal coherence is a n-fold convolution of the instrument responses and the expected coherence. The consequence is detrimental in which the photon number statistics is masked from the unresolved coherence signatures. The experimental investigations are thus far consistent with the theory developed. We envision the present theory will mitigate the false identification of optical emitters and enlarge the coherence deconvolution to an arbitrary order.

42 ENGINEERING↗
DOE OSTI2021

Scalable ab initio fragmentation methods based on a truncated expansion of the non-orthogonal molecular orbital model

An alternative formulation of the non-orthogonal molecular orbital model of electronic structure theory is developed based on the expansion of the inverse molecular orbital overlap matrix. From this model, a hierarchy of ab initio fragment-based quantum chemistry methods, referred to as the nth-order expanded non-orthogonal molecular orbital methods, are developed using a minimal number of approximations, each of which is frequently employed in intermolecular interaction theory. These novel methods are compared to existing fragment-based quantum chemistry methods, and the implications of those significant differences, where they exist, between the methods developed herein and those already existing methods are examined in detail. Computational complexities and theoretical scaling are also analyzed and discussed. Future extensions for the hierarchy of methods, to account for additional intrafragment and interfragment interactions, are outlined.

Westheimer, Bryce M. (ORCID:0000000267645714)↗
DOE OSTI2022

Orbital design of flat bands in non-line-graph lattices via line-graph wave functions

Line-graph (LG) lattices are known for having flat bands (FBs) from the destructive interference of Bloch wave functions encoded in only lattice symmetry. In this work, we develop a generic atomic/molecular orbital design principle for FBs in non-LG lattices. Based on linear combination of atomic orbital theory, we demonstrate that the underlying wave-function symmetry of FBs in a LG lattice can be transformed into the atomic/molecular orbital symmetry in a non-LG lattice. We illustrate such orbital-designed topological FBs in three 2D non-LG, square, trigonal, and hexagonal lattices, where the designed orbitals faithfully reproduce the corresponding lattice symmetries of checkerboard, kagome, and diatomic-kagome lattices, respectively. Interestingly, systematic design of FBs with a high Chern number is also achieved based on the same principle. Fundamentally our theory enriches the FB physics; practically, it significantly expands the scope of FB materials, since most materials have multiple atomic/molecular orbitals at each lattice site, rather than a single s orbital mandated in graph theory and generic lattice models.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Navigating collinear superspace

We introduce a new set of effective field theory rules for constructing Lagrangians with \( \mathcal{N} \) = 1 supersymmetry in collinear superspace. In the standard superspace treatment, superfields are functions of the coordinates \( \left({x}^{\mu },{\theta}^{\alpha },{\theta}^{\dagger \overset{\cdot }{\alpha }}\right) \) , and supersymmetry preservation is manifest at the Lagrangian level in part due to the inclusion of auxiliary F - and D -term components. By contrast, collinear superspace depends on a smaller set of coordinates ( x μ , η , η † ), where η is a complex Grassmann number without a spinor index. This provides a formulation of supersymmetric theories that depends exclusively on propagating degrees of freedom, at the expense of obscuring Lorentz invariance and introducing inverse momentum scales. After establishing the general framework, we construct collinear superspace Lagrangians for free chiral matter and non-Abelian gauge fields. For the latter construction, an important ingredient is a superfield representation that is simultaneously chiral, anti-chiral, and real; this novel object encodes residual gauge transformations on the light cone. Additionally, we discuss a fundamental obstruction to constructing inter- acting theories with chiral matter; overcoming these issues is the subject of our companion paper, where we introduce a larger set of superfields to realize the full range of interactions compatible with \( \mathcal{N} \) = 1. Along the way, we provide a novel framing of reparametrization invariance using a spinor decomposition, which provides insight into this important light-cone symmetry.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Prediction of Solute Segregation at Metal/Oxide Interfaces Using Machine Learning Approaches

The atomic structure and chemistry at metal/oxide interfaces play a crucial role in determining their properties. However, studying semi-coherent metal/oxide interfaces that include misfit dislocations through density functional theory (DFT) is often computationally expensive due to the large number of atoms involved, ranging from hundreds to thousands. In this study, we explore solute segregation behavior at the Fe/Y 2 O 3 interface—an important model interface for cladding applications in nuclear fission reactors—by combining DFT calculations with a machine learning (ML) approach. ML models are trained using DFT-calculated segregation energies (𝐸 𝑆𝑒𝑔 ) to identify the key chemical and geometric factors influencing solute segregation at metal/oxide interfaces, revealing the competition between these features in determining 𝐸 𝑆𝑒𝑔 . Moreover, the segregation behavior at a specific Fe/Y 2 O 3 interface is predicted with high accuracy using ML models trained on data from this interface. Furthermore, it is found that the ML models could also predict solute segregation at a different Fe/Y 2 O 3 interface with a new orientation relationship (OR), at a computational cost of less than 1/45 of that required for similar DFT calculations.

36 - MATERIALS SCIENCE↗