Search NASA⌕ Search

SEARCH · Search NASA

Results for “lie algebra”

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.

70 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Coupling and recoupling coefficients for Wigner’s U(4) supermultiplet symmetry

A novel procedure for evaluating Wigner coupling coefficients and Racah recoupling coefficients for U(4) in two group–subgroup chains is presented. The canonical U(4) > U(3) > U(2) > U(1) coupling and recoupling coefficients are applicable to any system that possesses U(4) symmetry, while the physical U(4) coupling coefficients are more specific to nuclear structure studies that utilize Wigner’s supermultiplet symmetry concept. The procedure that is proposed sidesteps the use of binomial coefficients and alternating sum series and consequently enables fast and accurate computation of any and all U(4)-underpinned features. The inner multiplicity of a (S, T) pair within a single U(4) > SU S (2) Ⓧ SU T (2) irreducible representation is obtained from the dimension of the null space of the SU(2) raising generators, while the resolution for the outer multiplicity follows from the work of Alex et al. on U(N) . It is anticipated that a C++ library will ultimately be available for determining generic coupling and recoupling coefficients associated with both the canonical and the physical group–subgroup chains of U(4).

Cross-Coupling Reaction↗
DOE OSTI2022

Nucleon-pair truncation of the shell model for medium-heavy nuclei

Background: Establishing computationally tractable models of atomic nuclei is a long-time goal of nuclear structure physics. A flexible framework which easily includes excited states and many-body correlations is the configuration-interaction shell model (SM), but the exponential growth of the basis means one needs an efficient truncation scheme, ideally one that includes both deformation and pairing correlations. Purpose: Here, we propose an efficient truncation scheme of the SM: starting from a pair condensate variationally defined by Hartree-Fock single-particle states and the particle-number conserved Bardeen-Cooper-Schrieffer (NBCS) approximation, we carry out projection of states with good angular momentum. Methods: After generating Hartree-Fock single-particle states with Kramers degeneracy in a SM space, we optimize the pair amplitudes in the NBCS by minimizing the energy, and then use linear algebra projection (LAP) of states with good angular momentum. Both NBCS and LAP are computationally fast. Results: Our calculations yield good agreement with full configuration-interaction SM calculations for low-lying states of transitional and rotational nuclei with axially symmetric and triaxial deformation in medium- and heavy-mass regions: 44,46,48 Ti, 48,50 Cr, 52 Fe, 60,62,64 Zn, 66,68 Ge, 68 Se, and 108,110 Xe. We predict low-lying states of 112–114 Ba and 116–120 Ce, nuclei difficult to reach by large-scale SM calculations. Conclusions: Both pair correlation and the configuration mixing between different intrinsic states play a key role in reproducing collectivity and shape coexistence, demonstrating the utility of this truncation scheme of the SM to study transitional and deformed nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Structure of Odd-A Ag Isotopes Studied via Algebraic Approaches

The structure of the odd-A silver isotopes Ag 103–115 is discussed within the frame of the interacting boson–fermion model (IBFM). An overview of their key properties is presented, with a particular attention paid to the “ J -1 anomaly”, represented by an abnormal ordering of the lowest 7/2 + and 9/2 + states. By examining previously published data and newly performed calculations, it is demonstrated that the experimentally known level schemes and electromagnetic properties of Ag 103–115 can be reproduced well within IBFM-1 by using a consistent set of model parameters. The contribution of different single-particle orbitals to the structure of the lowest-lying excited nuclear states in Ag 103–115 is discussed. Given that the J-1 anomaly brings down the 7/2 + level from the j −3 multiplet to energies, which can be thermally populated in hot stellar environments, the importance of low-lying excited states in odd-A silver isotopes for astrophysical processes is outlined.

IBFM-1↗
DOE OSTI2024

Open-closed string duality, branes, and topological recursion

We consider matrix models exhibiting open-closed string duality in two-dimensional string theories with various amounts of supersymmetry. In particular, a relationship between matrix models in the β = 2 Wigner-Dyson class and models in the (1 + 2Γ, 2) Altland-Zirnbauer class relates the perturbative solutions of the two systems’ string equations. Point-like operator insertions in the closed string theory are mapped to the topological expansion of the free energy in the open string theory. We compute correlation functions of macroscopic loop operators and FZZT branes in a general topological gravity background. The relationship between the topological recursion of moduli space volumes and branes is discussed by analyzing the Virasoro conditions in the matrix models.

2D Gravity↗
DOE OSTI2020

η-pairing in Hubbard models: From spectrum generating algebras to quantum many-body scars

In this work, we revisit the η-pairing states in Hubbard models and explore their connections to quantum many-body scars to discover a universal scars mechanism. η-pairing occurs due to an algebraic structure known as a spectrum generating algebra (SGA), giving rise to equally spaced towers of eigenstates in the spectrum. We generalize the original η-pairing construction and show that several Hubbard-like models on arbitrary graphs exhibit SGAs, including ones with disorder and spin-orbit coupling. We further define a restricted spectrum generating algebra (RSGA) and give examples of perturbations to the Hubbard-like models that preserve an equally spaced tower of the original model as eigenstates. The states of the surviving tower exhibit a subthermal entanglement entropy, and we analytically obtain parameter regimes for which they lie in the bulk of the spectrum, showing that they are exact quantum many-body scars. The RSGA framework also explains the equally spaced towers of eigenstates in several well-known models of quantum scars, including the Affleck-Kennedy-Lieb-Tasaki model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Design Considerations for GPU-based Mixed Integer Programming on Parallel Computing Platforms

Mixed Integer Programming (MIP) is a powerful abstraction in combinatorial optimization that finds real-life application across many significant sectors. The recent proliferation of graphical processing unit (GPU)-based accelerated computing architectures in large-scale parallel computing or supercomputing presents new opportunities as well as challenges in the advancement of MIP solver technology to effectively use the new accelerated computing platforms and scale to large parallel systems. Here, we recount the conventional processor-based strategies and focus on configurations where the most promising intersection lies between parallel MIP solver approaches and the specific strengths of accelerated parallel platforms. We note that the best potential lies in solving problems whose individual matrix sizes (of the linear program relaxation) fit entirely within one accelerator's memory and whose branch-and-bound (or branch-and-cut) trees cannot be fully contained within a small number of computational nodes. Additionally, we identify ideal features of computational linear algebra support on GPU accelerators that would help advance this direction of scalable parallel solution of MIP problems on GPU-based accelerated computing architectures.

Perumalla, Kalyan↗
DOE OSTI2022

Landau discriminants

Scattering amplitudes in quantum field theories have intricate analytic properties as functions of the energies and momenta of the scattered particles. In perturbation theory, their singularities are governed by a set of nonlinear polynomial equations, known as Landau equations, for each individual Feynman diagram. The singularity locus of the associated Feynman integral is made precise with the notion of the Landau discriminant, which characterizes when the Landau equations admit a solution. In order to compute this discriminant, we present approaches from classical elimination theory, as well as a numerical algorithm based on homotopy continuation. These methods allow us to compute Landau discriminants of various Feynman diagrams up to 3 loops, which were previously out of reach. For instance, the Landau discriminant of the envelope diagram is a reducible surface of degree 45 in the three-dimensional space of kinematic invariants. We investigate geometric properties of the Landau discriminant, such as irreducibility, dimension and degree. In particular, we find simple examples in which the Landau discriminant has codimension greater than one. Furthermore, we describe a numerical procedure for determining which parts of the Landau discriminant lie in the physical regions. In order to study degenerate limits of Landau equations and bounds on the degree of the Landau discriminant, we introduce Landau polytopes and study their facet structure. Finally, we provide an efficient numerical algorithm for the computation of the number of master integrals based on the connection to algebraic statistics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

From large to small $$ \mathcal{N} $$ = (4, 4) superconformal surface defects in holographic 6d SCFTs

Abstract Two-dimensional (2d)$$ \mathcal{N} $$ N = (4, 4) Lie superalgebras can be either “small” or “large”, meaning their R-symmetry is either$$ \mathfrak{so} $$ so (4) or$$ \mathfrak{so} $$ so (4) ⊕$$ \mathfrak{so} $$ so (4), respectively. Both cases admit a superconformal extension and fit into the one-parameter family$$ \mathfrak{d} $$ d (2,1;γ) ⊕$$ \mathfrak{d} $$ d (2,1;γ), with parameterγ∈ (−∞,∞). The large algebra corresponds to generic values ofγ, while the small case corresponds to a degeneration limit withγ→ −∞. In 11d supergravity, we study known solutions with superisometry algebra$$ \mathfrak{d} $$ d (2,1;γ) ⊕$$ \mathfrak{d} $$ d (2,1;γ) that are asymptotically locally AdS 7 ×𝕊 4 . These solutions are holographically dual to the 6d maximally superconformal field theory with 2d superconformal defects invariant under$$ \mathfrak{d} $$ d (2,1;γ) ⊕$$ \mathfrak{d} $$ d (2,1;γ). We show that a limit of these solutions, in whichγ→ −∞, reproduces another known class of solutions, holographically dual tosmall$$ \mathcal{N} $$ N = (4, 4) superconformal defects. We then use this limit to generate new small$$ \mathcal{N} $$ N = (4, 4) solutions with finite Ricci scalar, in contrast to the known small$$ \mathcal{N} $$ N = (4, 4) solutions. We then use holography to compute the entanglement entropy of a spherical region centered on these small$$ \mathcal{N} $$ N = (4, 4) defects, which provides a linear combination of defect Weyl anomaly coefficients that characterizes the number of defect-localized degrees of freedom. We also comment on the generalization of our results to include$$ \mathcal{N} $$ N = (0,4) surface defects through orbifolding.

Physics↗
DOE OSTI2023

Generalized dualities and supergroups

Abstract Using a recently developed formulation of double field theory in superspace, the graviton,B-field, gravitini, dilatini, and Ramond-Ramond bispinor are encoded in a single generalized supervielbein. Duality transformations are encoded as orthosymplectic transformations, extending the bosonic O(D,D) duality group, and these act on all constituents of the supervielbein in an easily computable way. We first review conventional non-abelian T-duality in the Green-Schwarz superstring and describe the dual geometries in the language of double superspace. Since dualities are related to super-Killing vectors, this includes as special cases both abelian and non-abelian fermionic T-duality. We then extend this approach to include Poisson-Lie T-duality and its generalizations, including the generalized coset construction recently discussed in [arXiv:1912.11036]. As an application, we construct the supergeometries associated with the integrableλandηdeformations of the AdS 5 × S 5 superstring. The deformation parametersλandηare identified with the possible one-parameter embeddings of the supergravity frame within the doubled supergeometry. In this framework, the Ramond-Ramond bispinors are directly computable purely from the algebraic data of the supergroup.

Physics↗
DOE OSTI2020

Massive celestial fermions

In an effort to further the study of amplitudes in the celestial CFT (CCFT), we construct conformal primary wavefunctions for massive fermions. Upon explicitly calculating the wavefunctions for Dirac fermions, we deduce the corresponding transformation of momentum space amplitudes to celestial amplitudes. The shadow wavefunctions are shown to have opposite spin and conformal dimension 2 - Δ. The Dirac conformal primary wave- functions are delta function normalizable with respect to the Dirac inner product provided they lie on the principal series with conformal dimension Δ = 1 + iλ for λ ∈ $\mathbb{R}$. It is shown that there are two choices of a complete basis: single spin J=$\frac{1}{2}$ or J=-$\frac{1}{2}$ and λ ∈ $\mathbb{R}$ or multiple spin J=±$\frac{1}{2}$ and λ $\mathbb{R}$ +∪0 . The massless limit of the Dirac conformal primary wavefunctions is shown to agree with previous literature. The momentum generators on the celestial sphere are derived and, along with the Lorentz generators, form a representation of the Poincaré algebra. Finally, we show that the massive spin-1 conformal primary wavefunctions can be constructed from the Dirac conformal primary wavefunctions using the standard Clebsch-Gordan coefficients. We use this procedure to write the massive spin-$\frac{3}{2}$, Rarita-Schwinger, conformal primary wavefunctions. This provides a prescription for constructing all massive fermionic and bosonic conformal primary wavefunctions starting from spin-$\frac{1}{2}$.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Scalable learning of potentials to predict time-dependent Hartree–Fock dynamics

We propose a framework to learn the time-dependent Hartree–Fock (TDHF) inter-electronic potential of a molecule from its electron density dynamics. Although the entire TDHF Hamiltonian, including the inter-electronic potential, can be computed from first principles, we use this problem as a testbed to develop strategies that can be applied to learn a priori unknown terms that arise in other methods/approaches to quantum dynamics, e.g., emerging problems such as learning exchange–correlation potentials for time-dependent density functional theory. We develop, train, and test three models of the TDHF inter-electronic potential, each parameterized by a four-index tensor of size up to 60 × 60 × 60 × 60. Two of the models preserve Hermitian symmetry, while one model preserves an eight-fold permutation symmetry that implies Hermitian symmetry. Across seven different molecular systems, we find that accounting for the deeper eight-fold symmetry leads to the best-performing model across three metrics: training efficiency, test set predictive power, and direct comparison of true and learned inter-electronic potentials. All three models, when trained on ensembles of field-free trajectories, generate accurate electron dynamics predictions even in a field-on regime that lies outside the training set. To enable our models to scale to large molecular systems, we derive expressions for Jacobian-vector products that enable iterative, matrix-free training.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Finite-size effects in periodic coupled cluster calculations

Here, we provide the first rigorous study of the finite-size error in the simplest and representative coupled cluster theory, namely the coupled cluster doubles (CCD) theory, for gapped periodic systems. Given exact Hartree-Fock orbitals and their corresponding orbital energies, we demonstrate that the correlation energy obtained from the approximate CCD method, after a finite number of fixed-point iterations over the amplitude equation, exhibits a finite-size error scaling as $\mathcal{O}(N^{-\frac{1}{3}}_k)$. Here $N_k$ is the number of discretization points in the Brillouin zone and characterizes the system size. Under additional assumptions ensuring the convergence of the fixed-point iterations, we demonstrate that the CCD correlation energy also exhibits a finite-size error scaling as $\mathcal{O}(N^{-\frac{1}{3}}_k)$. Our analysis shows that the dominant error lies in the coupled cluster amplitude calculation, and the convergence of the finite-size error in energy calculation can be boosted to $\mathcal{O}(N^{-1}_k)$ with accurate amplitudes. This also provides the first proof of the scaling of the finite-size error in the third order Møller-Plesset perturbation theory (MP3) for periodic systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A Flexible Power Method for Solving Infinite Dimensional Tensor Eigenvalue Problems

We propose a flexible power method for computing the leftmost, i.e., algebraically smallest, eigenvalue of an infinite dimensional tensor eigenvalue problem, $H x = \lambda x$, where the infinite dimensional symmetric matrix $H$ exhibits a translational invariant structure. We assume the smallest eigenvalue of $H$ is simple and apply a power iteration of $e^{-H}$ with the eigenvector represented in a compact way as a translational invariant infinite Tensor Ring (iTR). Hence, the infinite dimensional eigenvector can be represented by a finite number of iTR cores of finite rank. In order to implement this power iteration, we use a small parameter $t$ so that the infinite matrix-vector operation $e^{-Ht}x$ can efficiently be approximated by the Lie product formula, also known as Suzuki--Trotter splitting, and we employ a low rank approximation through a truncated singular value decomposition on the iTR cores in order to keep the cost of subsequent power iterations bounded. We also use an efficient way for computing the iTR Rayleigh quotient and introduce a finite size iTR residual which is used to monitor the convergence of the Rayleigh quotient and to modify the timestep $t$. In this paper, we discuss 2 different implementations of the flexible power algorithm and illustrate the automatic timestep adaption approach for several numerical examples.

Beeumen, Roel Van↗
DOE OSTI2025

Preserving nonlinear constraints in variational flow filtering data assimilation

Data assimilation aims to estimate the states of a dynamical system by optimally combining sparse and noisy observations of the physical system with uncertain forecasts produced by a computational model. The states of many dynamical systems of interest obey nonlinear physical constraints, and the corresponding dynamics is confined to a certain sub-manifold of the state space. Standard data assimilation techniques applied to such systems yield posterior states lying outside the manifold, violating the physical constraints. This work focuses on particle flow filters which use stochastic differential equations to evolve state samples from a prior distribution to samples from an observation-informed posterior distribution. The variational Fokker-Planck (VFP)—a generic particle flow filtering framework—is extended to incorporate non-linear, equality state constraints in the analysis. To this end, two algorithmic approaches that modify the VFP stochastic differential equation are discussed: (i) VFPSTAB, to inexactly preserve constraints with the addition of a stabilizing drift term, and (ii) VFPDAE, to exactly preserve constraints by treating the VFP dynamics as a stochastic differential-algebraic equation (SDAE). Additionally, an implicit-explicit time integrator is developed to evolve the VFPDAE dynamics. The strength of the proposed approach for constraint preservation in data assimilation is demonstrated on three test problems: the double pendulum, Korteweg-de-Vries, and the incompressible Navier-Stokes equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Constructing the infrared conformal generators on the fuzzy sphere

We investigate the conformal algebra on the fuzzy sphere, and in particular the generators of translations and special conformal transformations which are emergent symmetries in the infinite IR but are broken along the RG flow. We show how to extract these generators using the energy momentum tensor, which is complicated by the fact that one does not have a priori access to the energy momentum tensor of the CFT limit but rather must construct it numerically. We discuss and quantitatively analyze the main sources of corrections to the conformal generators due to the breaking of scale-invariance at finite energy, and develop efficient methods for removing these corrections. The resulting generators have matrix elements that match CFT predictions with accuracy varying from sub-percent level for the lowest-lying states up to several percent accuracy for states with dimension \sim 5 ∼ 5 with N=16 N = 16 fermions. We show that the generators can be used to accurately identify primary operators vs descendant operators in energy ranges where the spectrum is too dense to do the identification solely based on the approximate integer spacing within conformal multiplets.

Fardelli, Giulia (ORCID:0000000217998124)↗
DOE OSTI2021

Scattering amplitudes and simple canonical forms for simple polytopes

We provide an efficient recursive formula to compute the canonical forms of arbitrary d-dimensional simple polytopes, which are convex polytopes such that every vertex lies precisely on d facets. For illustration purposes, we explicitly derive recursive formulae for the canonical forms of Stokes polytopes, which play a similar role for a theory with quartic interaction as the Associahedron does in planar bi-adjoint φ 3 theory. As a by-product, our formula also suggests a new way to obtain the full planar amplitude in φ 4 theory by taking suitable limits of the canonical forms of constituent Stokes polytopes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗