Search NASA⌕ Search

SEARCH · Search NASA

Results for “symplectic”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Minimizing phase-space energies

A primary technical challenge for harnessing fusion energy is to control and extract energy from a nonthermal distribution of charged particles. The fact that phase space evolves by symplectomorphisms fundamentally limits how a distribution may be manipulated. While the constraint of phase-space volume preservation is well understood, other constraints remain to be fully appreciated. To better understand these constraints, we study the problem of extracting energy from a distribution of particles using area-preserving and symplectic linear maps. When a quadratic potential is imposed, we find that the maximal extractable energy can be computed as trace minimization problems. We solve these problems and show that the extractable energy under linear symplectomorphisms may be much smaller than the extractable energy under special linear maps. As a result, the method introduced in the present study enables an energy-based proof of the linear Gromov nonsqueezing theorem.

First-principles calculations in plasma physics↗
DOE OSTI2025

Achievable Rates for Concatenated Square Gottesman-Kitaev-Preskill Codes

The Gottesman-Kitaev-Preskill (GKP) codes are known to achieve optimal rates under displacement noise and pure-loss channels, which establishes theoretical foundations for its optimality. However, such optimal rates are only known to be achieved at a discrete set of noise strengths with the current self-dual symplectic lattice construction. In this work, we develop a new coding strategy using concatenated continuous variable-discrete variable encodings to go beyond past results and establish GKP’s optimal rate over all noise strengths. In particular, for displacement noise, the rate is obtained through a constructive approach by concatenating GKP codes with a quantum polar code and analog decoding. For a pure-loss channel, we prove the existence of capacity-achieving GKP codes through a random coding approach. These results highlight the capability of concatenation-based GKP codes and provides new methods for constructing good GKP lattices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Exact spectral gaps of random one-dimensional quantum circuits

The spectral gap of local random quantum circuits is a fundamental property that determines how close the moments of the circuit's unitaries match those of a Haar random distribution. When studying spectral gaps, it is common to bound these quantities using tools from statistical mechanics or via quantum information-based inequalities. Here, by focusing on the second moment of one-dimensional unitary circuits where nearest-neighboring gates act on sets of qudits (with open and closed boundary conditions), we show that one can exactly compute the associated spectral gaps. Indeed, having access to their functional form allows us to prove several important results, such as the fact that the spectral gap for closed boundary condition is exactly the square of the gap for open boundaries, as well as improve on previously known bounds for approximate design convergence. Finally, we verify our theoretical results by numerically computing the spectral gap for systems of up to 70 qubits, as well as comparing them to gaps of random orthogonal and symplectic circuits.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Time-dependent variational principle for mixed matrix product states in the thermodynamic limit

Here, we describe a time evolution algorithm for quantum spin chains whose Hamiltonians are composed of an infinite uniform left and right bulk part, and an arbitrary finite region in between. The left and right bulk parts are allowed to be different from each other. The algorithm is based on the time-dependent variational principle (TDVP) of matrix product states. It is inversion-free and very simple to adapt from an existing TDVP code for finite systems. The importance of working in the projective Hilbert space is highlighted. We study the quantum Ising model as a benchmark and an illustrative example. The spread of information after a local quench is studied in both the ballistic and the diffusive case. We also offer a derivation of TDVP directly from symplectic geometry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Geometric integration of classical spin dynamics via a mean-field Schrödinger equation

The Landau-Lifshitz equation describes the time evolution of magnetic dipoles and can be derived by taking the classical limit of a quantum mechanical spin Hamiltonian. To take this limit, one constrains the many-body quantum state to a tensor product of coherent states, thereby neglecting entanglement between sites. Expectation values of the quantum spin operators produce the usual classical spin dipoles. One may also consider expectation values of polynomials of the spin operators, leading to quadrupole and higher-order spin moments, which satisfy a dynamical equation of motion that generalizes the Landau-Lifshitz dynamics [Zhang and Batista, Phys. Rev. B 104, 104409 (2021)]. Here we reformulate the dynamics of these N 2 –1 generalized spin components as a mean-field Schrödinger equation on the N-dimensional coherent state. Furthermore, this viewpoint suggests efficient integration methods that respect the local symplectic structure of the classical spin dynamics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Quantum Liouville theorem based on Haar measure

Liouville theorem (L theorem) reveals robust incompressibility of the distribution function in phase space, given arbitrary potentials. However, its quantum generalization, Wigner flow, is compressible, i.e., L theorem is only conditionally true (e.g., for perfect Harmonic potential). Here, we develop quantum L theorem (rigorous incompressibility) for arbitrary potentials (interacting or not) in Hamiltonians. Haar measure, instead of symplectic measure dp$\bigwedge$dq used in Wigner’s scheme, plays a central role. The argument is based on general measure theory, independent of specific spaces or coordinates. Comparison of classical and quantum is made: for instance, here we address why Haar measure and metric preservation do not work in the classical case. Applications of the theorems in statistics, topological phase transition, ergodic theory, etc., are discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Weak localization as a probe of intervalley coherence in graphene multilayers

Spontaneous intervalley coherence is suspected in several different graphene multilayer systems, but is difficult to confirm because of a paucity of convenient experimental signatures. Here we suggest that magneto-conductance features associated with quantum corrections to Drude conductivity can serve as a smoking gun for intervalley coherence that does not break time-reversal symmetry. In this class of ordered multilayer quantum transport corrections can produce weak localization or weak antilocalization, depending on whether the valley order belongs to the orthogonal or symplectic symmetry class. Our analysis motivates low-temperature weak-field magnetoresistance measurements in graphene multilayers in which time-reversal invariant intervalley coherent order is conjectured.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Clustering and α-capture reaction rate from ab initio 20 Ne

We introduce a new framework for studying clustering and for calculating α partial widths using ab initio wave functions. We demonstrate the formalism for 20 Ne, by calculating the overlap between the 16 O + α cluster configuration and states in 20 Ne computed in the ab initio symmetry-adapted no-core shell model. We present spectroscopic amplitudes and spectroscopic factors, and compare those to no-core symplectic shell-model results in larger model spaces, to gain insight into the underlying physics that drives α clustering. Specifically, we report on the α partial width of the lowest 1 – resonance in 20 Ne, which is found to be in good agreement with experiment.We also present first no-core shell-model estimates for asymptotic normalization coefficients for the ground state, as well as for the first excited 4 + state in 20 Ne that lies in a close proximity to the α + 16 O threshold. This outcome highlights the importance of correlations for developing cluster structures and for describing α widths. The widths can then be used to calculate α-capture reaction rates for narrow resonances of interest to astrophysics. Finally, we explore the reaction rate for the α-capture reaction 16 O(α, γ ) 20 Ne at astrophysically relevant temperatures and determine its impact on simulated x-ray burst abundances.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Discrete gravity with local Lorentz invariance

A novel structure-preserving algorithm for general relativity in vacuum is derived from a lattice gauge theoretic discretization of the tetradic Palatini action. Here, the resulting model of discrete gravity is demonstrated to preserve local Lorentz invariance and symplectic structure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

U(1) fields from qubits: An approach via D-theory algebra

A new quantum link microstructure was proposed for the lattice quantum chromodynamics (QCD) Hamiltonian, replacing the Wilson gauge links with a bilinear of fermionic qubits, later generalized to D-theory. This formalism provides a general framework for building lattice field theory algorithms for quantum computing. We focus mostly on the simplest case of a quantum rotor for a single compact U(1) field. We also make some progress for non-Abelian setups, making it clear that the ideas developed in the U(1) case extend to other groups. These in turn are building blocks for 1 + 0 -dimensional ( 1 + 0 -D) matrix models, 1 + 1 -D sigma models and non-Abelian gauge theories in 2 + 1 and 3 + 1 dimensions. By introducing multiple flavors for the U(1) field, where the flavor symmetry is gauged, we can efficiently approach the infinite-dimensional Hilbert space of the quantum O(2) rotor with increasing flavors. The emphasis of the method is on preserving the symplectic algebra exchanging fermionic qubits by sigma matrices (or hard bosons) and developing a formal strategy capable of generalization to a SU ( 3 ) field for lattice QCD and other non-Abelian 1 + 1 -D sigma models or 3 + 1 -D gauge theories. For U(1), we discuss briefly the qubit algorithms for the study of the discrete 1 + 1 -D sine-Gordon equation. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2025

Gromov ground state in phase space engineering for fusion energy

Phase space engineering by rf waves plays important roles in both thermal D-T fusion and nonthermal advanced fuel fusion, but not all phase space manipulation is allowed; certain fundamental limits exist. In addition to Liouville's theorem, which requires the manipulation to be volume preserving, Gromov's nonsqueezing theorem imposes another constraint. Here, the Gardner ground state is defined as the ground state accessible by smooth volume-preserving maps. However, the extra Gromov constraint should produce a higher-energy ground state. An example of a Gardner ground state forbidden by Gromov's nonsqueezing theorem is given. The challenge question is “What is the Gromov ground state, i.e., the lowest energy state accessible by smooth symplectic maps?” This is a difficult problem. As a simplification, we conjecture that the linear Gromov ground state problem is solvable.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Emergent Sp(3,$\mathbb{R}$) Dynamical Symmetry in the Nuclear Many-Body System from an Ab Initio Description

Ab initio nuclear theory provides not only a microscopic framework for quantitative description of the nuclear many-body system, but also a foundation for deeper understanding of emergent collective correlations. In this work, a symplectic Sp(3, $\mathbb{R}$) $\supset$ U(3) dynamical symmetry is identified in ab initio predictions, from a no-core configuration interaction approach, and found to provide a qualitative understanding of the spectrum of 7 Be. Low-lying states form an Elliott SU(3) spectrum, while an Sp(3, $\mathbb{R}$) excitation gives rise to an excited rotational band with strong quadrupole connections to the ground state band.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Gate-Based Quantum Simulation of Gaussian Bosonic Circuits on Exponentially Many Modes

We introduce a framework for simulating, on an ( n + 1 )-qubit quantum computer, the action of a Gaussian bosonic (GB) circuit on a state over 2 n modes. Specifically, we encode the initial bosonic state’s expectation values over quadrature operators (and their covariance matrix) as an input qubit state. This is then evolved by a quantum circuit that effectively implements the symplectic propagators induced by the GB gates. We find families of GB circuits and initial states leading to efficient quantum simulations. For this purpose, we introduce a dictionary that maps between GB and qubit gates such that particle- (non-particle-) preserving GB gates lead to real- (imaginary-) time evolutions at the qubit level. For the special case of particle-preserving circuits, we present a bounded-error-quantum-polynomial time (BQP)-complete GB decision problem, indicating that GB evolutions of Gaussian states on exponentially many modes are as powerful as universal quantum computers. We also perform numerical simulations of an interferometer on ∼ 8 × 10 9 modes, illustrating the power of our framework. Published by the American Physical Society 2025

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Theory of quantum circuits with Abelian symmetries

Quantum circuits with gates (local unitaries) respecting a global symmetry have broad applications in quantum information science and related fields, such as condensed-matter theory and quantum thermodynamics. However, despite their widespread use, fundamental properties of such circuits are not well understood. Recently, it was found that generic unitaries respecting a global symmetry cannot be realized, even approximately, using gates that respect the same symmetry. This observation raises important open questions: What unitary transformations can be realized with k -local gates that respect a global symmetry? In other words, in the presence of a global symmetry, how does the locality of interactions constrain the possible time evolution of a composite system? In this work, we address these questions for the case of Abelian (commutative) symmetries and develop constructive methods for synthesizing circuits with such symmetries. Remarkably, as a corollary, we find that, while the locality of interactions still imposes additional constraints on realizable unitaries, certain restrictions observed in the case of non-Abelian symmetries do not apply to circuits with Abelian symmetries. For instance, in circuits with a general non-Abelian symmetry such as SU ( d ) , the unitary realized in a subspace with one irreducible representation (charge) of the symmetry dictates the realized unitaries in multiple other sectors with inequivalent representations of the symmetry. Furthermore, in certain sectors, rather than all unitaries respecting the symmetry, the realizable unitaries are the symplectic or orthogonal subgroups of this group. We prove that none of these restrictions appears in the case of Abelian symmetries. This result suggests that global non-Abelian symmetries may affect the thermalization of quantum systems in ways not possible under Abelian symmetries. Published by the American Physical Society 2024

Marvian, Iman↗
DOE OSTI2026

New Time Integrators and Capabilities in SUNDIALS Versions 6.2.0-7.4.0

SUNDIALS is a well-established numerical library that provides robust and efficient time integrators and nonlinear solvers. This article overviews several significant improvements and new features added over the last 3 years to support scientific simulations run on high-performance computing systems. Notably, three new classes of one-step methods have been implemented: low storage Runge–Kutta, symplectic partitioned Runge–Kutta, and operator splitting. In addition, we describe new timestep adaptivity support for multirate methods, adjoint sensitivity analysis capabilities for explicit Runge–Kutta methods, additional options for Anderson acceleration in nonlinear solvers, and improved error handling and logging.

Computer science↗
DOE OSTI2021

Nuclear Dynamics and Reactions in the Ab Initio Symmetry-Adapted Framework

We review the ab initio symmetry-adapted (SA) framework for determining the structure of stable and unstable nuclei, along with related electroweak, decay, and reaction processes. This framework utilizes the dominant symmetry of nuclear dynamics, the shape-related symplectic [Formula: see text] symmetry, which has been shown to emerge from first principles and to expose dominant degrees of freedom that are collective in nature, even in the lightest species or seemingly spherical states. This feature is illustrated for a broad scope of nuclei ranging from helium to titanium isotopes, enabled by recent developments of the ab initio SA no-core shell model expanded to the continuum through the use of the SA basis and that of the resonating group method. The review focuses on energies, electromagnetic transitions, quadrupole and magnetic moments, radii, form factors, and response function moments for ground-state rotational bands and giant resonances. The method also determines the structure of reaction fragments that is used to calculate decay widths and α-capture reactions for simulated X-ray burst abundance patterns, as well as nucleon–nucleus interactions for cross sections and other reaction observables.

Physics↗
DOE OSTI2022

BIhNNs

The code enables to perform Bayesian inference in an efficient manner through the use of Hamiltonian Neural Networks (HNNs), Deep Neural Networks (DNNs), Neural ODEs, and Symplectic Neural Networks (SympNets) used with state-of-the-art sampling schemes like Hamiltonian Monte Carlo (HMC) and the No-U-Turn-Sampler (NUTS).

Dhulipala, Som↗
DOE OSTI2022

6-D ELEMENT-BY-ELEMENT PARTICLE TRACKING WITH CRAB CAVITY PHASE NOISE AND WEAK-STRONG BEAM-BEAM INTERACTION FOR THE HADRON STORAGE RING OF THE ELECTRON-ION COLLIDER

The Electron Ion Collider (EIC) presently under construction at Brookhaven National Laboratory will collide polarized high energy electron beams with hadron beams with luminosities up to 10^34 cm?^2 s^?1 in center mass energy range of 20-140 GeV. Crab cavities are used to compensate the geometric luminosity due to a large crossing angle 25 mrad in the EIC. It was found that the phase noise in crab cavities will generate a significant emittance growth for the hadron beams. The phase noise tolerance from the analytical calculation is very small for the Hadron Storage Ring (HSR) of the EIC. In this article, we present our preliminary numerical simulation results with a 6-d symplectic particle tracking to determine the proton emittance growth rate with crab cavity phase noise and beam-beam interaction for the HSR.

Luo, Y.↗