SEARCH · Engineering Papers

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

SourceResearch connectionFrom yearThrough year

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

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Automated label flows for excited states of correlation functions in lattice gauge theory

Extracting excited states from lattice gauge theory correlation functions can be achieved through $χ^2$ minimization fits or algebraic approaches such as the variational method and Prony's method. Performing any kind of error analysis often leads to overlapping confidence regions of model parameters, even when the spectrum is not particularly dense. To correctly estimate errors, one must beware of mislabeling the states. In this work, we provide an algorithm that we call automated label flows which consistently and systematically identifies a deterministic labeling of states. This is a black-box approach in the sense that it gives a sensible set of labels without user guidance. Further, as an example, we pair one black box method with another, analyzing a lattice correlation function from real data using automated label flows in the context of Prony's method.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

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 OSTI2025

String membrane nets from higher-form gauging: An alternate route to 𝑝-string condensation

We present a new perspective on the $p$-string condensation procedure for constructing 3+1D fracton phases by implementing this process via the gauging of higher-form symmetries. Specifically, we show that gauging a 1-form symmetry in 3+1D that is generated by Abelian anyons in isotropic stacks of 2+1D topological orders naturally results in a 3+1D $p$-string condensed phase, providing a controlled non-perturbative construction that realizes fracton orders. This approach clarifies the symmetry principles underlying $p$-string condensation and generalizes the familiar connection between anyon condensation and one-form gauging in two spatial dimensions. We demonstrate this correspondence explicitly in both field theories and lattice models: in field theory, we derive the foliated field theory description of the $\mathbb{Z}_N$ X-Cube model by gauging a higher-form symmetry in stacks of 2+1D $\mathbb{Z}_N$ gauge theories; on the lattice, we show how gauging a diagonal 1-form symmetry in isotropic stacks of $G$-graded string-net models leads to string-membrane-nets hosting restricted mobility excitations. This perspective naturally generalizes to spatial dimensions $d \geq 2$ and provides a step towards building an algebraic theory of $p$-string condensation.

Anyons↗

NASA NTRS2008

Application of Fast Multipole Methods to the NASA Fast Scattering Code

The NASA Fast Scattering Code (FSC) is a versatile noise prediction program designed to conduct aeroacoustic noise reduction studies. The equivalent source method is used to solve an exterior Helmholtz boundary value problem with an impedance type boundary condition. The solution process in FSC v2.0 requires direct manipulation of a large, dense system of linear equations, limiting the applicability of the code to small scales and/or moderate excitation frequencies. Recent advances in the use of Fast Multipole Methods (FMM) for solving scattering problems, coupled with sparse linear algebra techniques, suggest that a substantial reduction in computer resource utilization over conventional solution approaches can be obtained. Implementation of the single level FMM (SLFMM) and a variant of the Conjugate Gradient Method (CGM) into the FSC is discussed in this paper. The culmination of this effort, FSC v3.0, was used to generate solutions for three configurations of interest. Benchmarking against previously obtained simulations indicate that a twenty-fold reduction in computational memory and up to a four-fold reduction in computer time have been achieved on a single processor.

Dunn, Mark H.↗

DOE OSTI2020

Machine learning Frenkel Hamiltonian parameters to accelerate simulations of exciton dynamics

In this manuscript, we develop multiple machine learning (ML) models to accelerate a scheme for parameterizing site-based models of exciton dynamics from all-atom configurations of condensed phase sexithiophene systems. This scheme encodes the details of a system’s specific molecular morphology in the correlated distributions of model parameters through the analysis of many single-molecule excited-state electronic-structure calculations. These calculations yield excitation energies for each molecule in the system and the network of pair-wise intermolecular electronic couplings. In this work, we demonstrate that the excitation energies can be accurately predicted using a kernel ridge regression (KRR) model with Coulomb matrix featurization. We present two ML models for predicting intermolecular couplings. The first one utilizes a deep neural network and bi-molecular featurization to predict the coupling directly, which we find to perform poorly. The second one utilizes a KRR model to predict unimolecular transition densities, which can subsequently be analyzed to compute the coupling. We find that the latter approach performs excellently, indicating that an effective, generalizable strategy for predicting simple bimolecular properties is through the indirect application of ML to predict higher-order unimolecular properties. Such an approach necessitates a much smaller feature space and can incorporate the insight of well-established molecular physics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Quantum Dynamics at Scale: Ultrafast Control of Emergent Functional Materials

Confluence of extreme-scale quantum dynamics simulations (i.e. quantum@scale) and cutting-edge x-ray free-electron laser experiments are revolutionizing materials science. An archetypal example is the exciting concept of using picosecond light pulses to control emergent material properties on demand in atomically-thin layered materials. This paper describes efforts to scale our quantum molecular dynamics engine toward the United States' first exaflop/s computer, under an Aurora Early Science Program project named "Metascalable layered material genome". Key algorithmic and computing techniques incorporated are: (1) globally-scalable and locally-fast solvers within a linear-scaling divide-conquer-recombine algorithmic framework; (2) algebraic `BLASification' of computational kernels; and (3) data alignment and loop restructuring, along with register and cache blocking, for enhanced vectorization and efficient memory access. The resulting weak-scaling parallel efficiency was 0.93 on 131,072 Intel Xeon Phi cores for a 56.6 million atom (or 169 million valence-electron) system, whereas the various code transformations achieved 5-fold speedup. The optimized simulation engine allowed us for the first time to establish a significant effect of substrate on the dynamics of layered material upon electronic excitation.

Algebraic BLASification↗

NASA NTRS1988

Sensitivity of control-augmented structure obtained by a system decomposition method

The verification of a method for computing sensitivity derivatives of a coupled system is presented. The method deals with a system whose analysis can be partitioned into subsets that correspond to disciplines and/or physical subsystems that exchange input-output data with each other. The method uses the partial sensitivity derivatives of the output with respect to input obtained for each subset separately to assemble a set of linear, simultaneous, algebraic equations that are solved for the derivatives of the coupled system response. This sensitivity analysis is verified using an example of a cantilever beam augmented with an active control system to limit the beam's dynamic displacements under an excitation force. The verification shows good agreement of the method with reference data obtained by a finite difference technique involving entire system analysis. The usefulness of a system sensitivity method in optimization applications by employing a piecewise-linear approach to the same numerical example is demonstrated. The method's principal merits are its intrinsically superior accuracy in comparison with the finite difference technique, and its compatibility with the traditional division of work in complex engineering tasks among specialty groups.

Sobieszczanskisobieski, Jaroslaw↗

NASA NTRS1991

Sensitivity of control-augmented structure obtained by a system decomposition method

Sobieszczanski-Sobieski, Jaroslaw↗

NASA NTRS1990

On the receptivity and non-parallel stability of travelling disturbances in rotating disk flow

The generation and evolution of small amplitude wavelength traveling disturbances in rotating disk flow is discussed. The steady rotational speed of the disk is perturbed so as to introduce high frequency oscillations in the flow field. Secondly, surface imperfections are introduced on the disk such as roughness elements. The interaction of these two disturbances will generate the instability waves whose evolution is governed by parabolic partial differential equations that are solved numerically. For the class of disturbances considered (wavelength on the order of Reynolds number), it is found that eigensolutions exist which decay or grow algebraically in the radial direction. However, these solutions grow only for frequencies larger than 4.58 times the steady rotational speed of the disk. The computed receptivity coefficient shows that there is an optimum size of roughness for which these modes are excited the most. The width of these roughness elements in the radial direction is about .1 r(sub 0) where r(sub 0) is the radial location of the roughness. It is also found that the receptivity coefficient is larger for a negative spanwise wavenumber than for a positive one. Typical wave angles found for these disturbances are about -26 degrees.

Balakumar, P.↗

NASA NTRS1990

Nonlinear spatial evolution of externally excited instability waves in free shear layers

The effects of critical-layer nonlinearity on spatially growing instability waves on shear layers between parallel streams are considered. In the two-dimensional incompressible case, the flow in the critical layer is governed by a nonequilibrium, nonlinear vorticity equation. The initial exponential growth of the instability wave is converted into algebraic growth during the streamwise 'aging' of the critical layer into a quasi-equilibrium state. This leads to a next stage of evolution where the instability-wave growth is affected by both mean-flow-divergence and nonlinear effects and is eventually converted to decay. A uniformly valid composite formula for the instability wave amplitude, accounting for both nonparallel and nonlinear effects, is given and compared to experimental results. The inviscid solutions to these equations always end in a singularity at finite downstream distance.

Hultgren, L. S.↗

DOE OSTI2020

Controllable vortex shedding from dissipative exchange flows in ferromagnetic channels

Ferromagnetic channels subject to spin injection at one extremum sustain long-range coherent textures that carry spin currents known as dissipative exchange flows (DEFs). In the weak injection regime, spin currents carried by DEFs decay algebraically and extend through the length of the channel, a regime known as spin superfluidity. Similar to fluids, these structures are prone to phase slips that manifest as vortex-antivortex pairs. In this work, we numerically study vortex shedding from DEFs excited in a magnetic nanowire with a physical obstacle. Using micromagnetic simulations, we find regimes of laminar flow and vortex shedding as a function of obstacle position tunable by the spin injection sign and magnitude. Vortex-antivortex pairs translate forward (VF regime) or backward (VB regime) with respect to the detector's extremum, resulting in well-defined spectral features. Qualitatively similar results are obtained when temperature, anisotropy, and weak nonlocal dipole fields are included in the simulations. These results provide clear features associated with DEFs that may be detected experimentally in devices with nominally identical boundary conditions. Furthermore, our results suggest that obstacles can be considered as DEF control gates, opening an avenue to manipulate DEFs via physical defects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Almost strong zero modes at finite temperature

Interacting fermionic chains exhibit extended regions of topological degeneracy of their ground states as a result of the presence of Majorana or parafermionic zero modes localized at the edges. In the opposite limit of infinite temperature, the corresponding nonintegrable spin chains, obtained via generalized Jordan-Wigner mapping, are known to host so-called almost strong zero modes, which are long-lived with respect to any bulk excitations. Here we study the fairly unexplored territory that bridges these two extreme cases of zero and infinite temperature. We blend two established techniques for states, the Lanczos series expansion and a tensor network ansatz, uplifting them to the level of operator algebra. This allows us to efficiently simulate large system sizes for arbitrarily long timescales and to extract the temperature-dependent decay rates. We observe that for the Kitaev-Hubbard model, the decay rate of the edge mode depends exponentially on the inverse temperature 𝛽, and on an effective energy scale Δ eff that is greater than the thermodynamic gap of the system Δ.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Beam Dynamics of the Muon $g\textrm{-}2$ Experiment

The Muon $g\textrm{-}2$ Experiment (E989) at Fermilab aims to measure the muon anomalous magnetic moment $a_{\mu}$ with unprecedented precision, potentially uncovering physics beyond the Standard Model of particle physics. The result based on Runs 1-3, released in 2023, achieved a precision of 0.20 ppm. The experiment circulates muons in a storage ring, measuring $a_{\mu}$ from decay positron time and energy measurements collected with calorimeters. To achieve the required accuracy, it is crucial to measure and control the magnetic field in the ring with high precision. Beam dynamics corrections are necessary for muons not orbiting exactly in the midplane, for their oscillations, and for electric field effects. Highly accurate beam dynamics simulations are instrumental for quantifying and validating the beam dynamics corrections, ultimately improving the precision of the $a_{\mu}$ measurement and facilitating the achievement of the ambitious $70\:\mathrm{ppb}$ systematic uncertainty goal. The measured field data was incorporated into models for simulations using three codes: \texttt{gm2ringsim} (an internal Geant4-based code), \textit{COSY INFINITY}, and \textit{BMAD}. The advantages of \texttt{gm2ringsim} include using CAD-based geometry and modelling the detector effects. \textit{COSY INFINITY} is a highly accurate and efficient code that uses high-order differential-algebraic transfer maps, precise fringe field calculations, and advanced symplectification methods. Symplectification is important for maintaining the physical correctness of the muon beam behaviour with high precision over the storage time, ensuring conservation of phase space volume and preventing artificial damping or excitation of particle motion. The experiment completed its final Run 6 in July 2023, collecting 21 times more data than the previous BNL experiment. Analyses of data from Runs 4-6 are ongoing, with results planned for release in 2025, potentially resolving the current tension between experiment and theory.

43 PARTICLE ACCELERATORS↗

NASA NTRS1991

On the receptivity and nonparallel stability of traveling disturbances in rotating-disk flow

The generation and evolution of small amplitude long wavelength traveling disturbances in rotating-disk flow are discussed. The steady rotational speed of the disk is perturbed so as to introduce high-frequency oscillations in the flow field. Secondly, surface imperfections on the disk such as roughness elements are introduced. The interaction of these two disturbances will generate the instability waves whose evolution is governed by parabolic partial differential equations which are solved numerically. It is found that, for the class of disturbances considered here (wavelength on the order of Reynolds number), eigensolutions exist which decay or grow algebraically in the radial direction. However, these solutions grow only for frequencies larger than 4.58 times the steady rotational speed of the disk. The computed receptivity coefficient shows that there is an optimum size of roughness for which these modes are preferentially excited.

Balakumar, P.↗

NASA NTRS1991

Calculations of steady and transient channel flows with a time-accurate L-U factorization scheme

Calculations of steady and unsteady, transonic, turbulent channel flows with a time accurate, lower-upper (L-U) factorization scheme are presented. The L-U factorization scheme is formally second-order accurate in time and space, and it is an extension of the steady state flow solver (RPLUS) used extensively to solve compressible flows. A time discretization method and the implementation of a consistent boundary condition specific to the L-U factorization scheme are also presented. The turbulence is described by the Baldwin-Lomax algebraic turbulence model. The present L-U scheme yields stable numerical results with the use of much smaller artificial dissipations than those used in the previous steady flow solver for steady and unsteady channel flows. The capability to solve time dependent flows is shown by solving very weakly excited and strongly excited, forced oscillatory, channel flows.

Kim, S.-W.↗

NASA NTRS1991

Post-glacial relaxation of a viscously stratified compressible mantle

The postglacial relaxation phenomenon of a viscously compressible nonself-gravitating spherical-shell model is investigated. Analytical solutions to an exponentially depth-dependent viscosity are developed with various types of compressible density models such as those with an exponential dependence of the radius and an algebraic-root dependence of the radius. The solutions of a viscously compressible multilayered model with constant thermodynamic properties and viscosity are developed by the propagator matrix method. The results show that inferences of deep mantle viscosity from postglacial rebound would be hampered by mantle compressibility for long-wavelength harmonics because of the smaller excitation of compressible eigenfunctions in the lower mantle. The relaxation times and velocity fields are more sensitive to higher viscosity contrast for the exponentially varying viscosity than for models with discrete jumps in the viscosity structure.

Wu, J.↗

DOE OSTI2025

Closed-form expressions for unitaries of spin-adapted fermionic operators

One of the open challenges in quantum computing simulations of problems of chemical interest is the proper enforcement of spin symmetry. Efficient quantum circuits implementing unitaries generated by spin-adapted operators remain elusive, while naïve Trotterization schemes break spin symmetry. Here, in this work, we analyze the mathematical structure of spin-adapted operators and derive closed-form expressions for unitaries generated by singlet spin-adapted generalised single and double excitations. These results represent significant progress toward the economical enforcement of spin symmetry in quantum simulations.

fermionic algebra↗

DOE OSTI2020

SU(1,1) Echoes for Breathers in Quantum Gases

Though the celebrated spin echoes have been widely used to reverse quantum dynamics, they are not applicable to systems whose constituents are beyond the control of the su(2) algebra. In this work, we design echoes to reverse quantum dynamics of breathers in three-dimensional unitary fermions and two-dimensional bosons and fermions with contact interactions, which are governed by an underlying su(1, 1) algebra. Geometrically, SU(1, 1) echoes produce closed trajectories on a single or multiple Poincaré disks and thus could recover any initial states without changing the sign of the Hamiltonian. In particular, the initial shape of a breather determines the superposition of trajectories on multiple Poincaré disks and whether the revival time has period multiplication. Our work provides physicists with a recipe to tailor collective excitations of interacting many-body systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗