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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Continuum electromagnetic gyrokinetic simulations of turbulence in the tokamak scrape-off layer and laboratory devices

Here, we present algorithms and results from Gkeyll, a full-f continuum, electromagnetic gyrokinetic code, designed to study turbulence in the edge region of fusion devices. The edge is computationally very challenging, requiring robust algorithms that can handle large-amplitude fluctuations and stable interactions with plasma sheaths. We present an energy-conserving high-order discontinuous Galerkin scheme that solves gyrokinetic equations in Hamiltonian form. Efficiency is improved by a careful choice of basis functions and automatically generated computation kernels. Previous verification tests were performed in the straight-field-line large plasma device [Shi et al., J. Plasma Phys. 83, 905830304 (2017)] and the Texas Helimak, a simple magnetized torus [Bernard et al., Phys. Plasmas 26, 042301 (2019)], including the effect of end-plate biasing on turbulence. Results for the scrape-off layer for NSTX parameters with a model helical magnetic geometry with bad curvature have been obtained [Shi et al., Phys. Plasmas 26, 012307 (2019)]. In this paper, we present algorithms for the two formulations of electromagnetic gyrokinetics: the Hamiltonian and the symplectic. We describe each formulation and show results of benchmark tests. Although our scheme works for the Hamiltonian formulation, the presence of spurious numerical modes for high-β and large $k^2_⊥ρ^2_s$ regimes shows that the symplectic formulation is more robust. We then review our recent algorithm for the symplectic formulation [Mandell et al., J. Plasma Phys. 86, 905860109 (2020)], along with example application of this new capability. Maintaining positivity of the distribution function can be challenging, and we describe a new and novel exponential recovery based algorithm to address this.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Quantum area fluctuations from gravitational phase space

We study the gravitational phase space associated to a stretched horizon within a finite-sized causal diamond in (d + 2)-dimensional spacetimes. By imposing the Raychaudhuri equation, we obtain its constrained symplectic form using the covariant phase space formalism and derive the relevant quantum commutators by inverting the symplectic form and quantizing. Finally, we compute the area fluctuations of the causal diamond by taking a Carrollian limit of the stretched horizon in pure Minkowski spacetime, and derive the relationship $\left\langle {\left(\Delta A\right)}^2\right\rangle \ge \frac{2\pi G}{d}\left\langle A\right\rangle$, showing that the variance of the area fluctuations is proportional to the area itself.

classical theories of gravity↗
DOE OSTI2024

Dynamical edge modes and entanglement in Maxwell theory

Abstract Previous work on black hole partition functions and entanglement entropy suggests the existence of “edge” degrees of freedom living on the (stretched) horizon. We identify a local and “shrinkable” boundary condition on the stretched horizon that gives rise to such degrees of freedom. They can be interpreted as the Goldstone bosons of gauge transformations supported on the boundary, with the electric field component normal to the boundary as their symplectic conjugate. Applying the covariant phase space formalism for manifolds with boundary, we show that both the symplectic form and Hamiltonian exhibit a bulk-edge split. We then show that the thermal edge partition function is that of a codimension-two ghost compact scalar living on the horizon. In the context of a de Sitter static patch, this agrees with the edge partition functions found by Anninos et al. in arbitrary dimensions. It also yields a 4D entanglement entropy consistent with the conformal anomaly. Generalizing to Proca theory, we find that the prescription of Donnelly and Wall reproduces existing results for its edge partition function, while its classical phase space does not exhibit a bulk-edge split.

Physics↗
DOE OSTI2019

Chapter Six - Derivation, cross-validation, and comparison of analytic formulas for electrostatic deflector aberrations

Here, we derived first and second order analytic aberration formulas in the horizontal transverse plane for an electrostatic deflector specified by the reference orbit radius, the central angle spanning the deflector, and its inhomogeneity coefficients. The derivation was performed using an iterative order-by-order perturbation method in a Frenet–Serret beamline coordinate system. We produced a C program edabrt for calculation of the first and second order aberrations using the formulas that we derived. The electrostatic deflector aberration formulas from Wollnik (1965) disagree with COSY INFINITY and the aberration formulas that we derived here, and they also deviate from the first and second order symplecticity conditions. The reason for this discrepancy is that the aberration formulas from Wollnik (1965) consider only motion in the main field of the deflector and do not properly account for fringe field effects, in particular the necessarily occurring change of kinetic energy due to the change in potential. The code GIOS uses electrostatic aberration formulas from Wollnik (1965) and exhibits the same discrepancy. We performed a comparison of first and second order electrostatic deflector aberrations for (1) the analytic formulas derived in this work, (2) differential-algebraic (DA) numerical integration of the equations of motion using the code COSY INFINITY, (3) the aberration formulas from Wollnik (1965), (4) aberrations computed using the code GIOS, and (5) the aberration formulas from Wollnik (1965) adjusted to account for the occurring change in potential. An electrostatic spherical deflector and an electrostatic cylindrical deflector were used as test cases for this comparison. There is excellent agreement between methods (1), (2), and (5), and these three methods satisfy the first and second order symplecticity conditions.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Slow manifolds of classical Pauli particle enable structure-preserving geometric algorithms for guiding center dynamics

Since variational symplectic integrators for the guiding center was proposed, structure-preserving geometric algorithms have become an active research field in plasma physics. We found that the slow manifolds of the classical Pauli particle enable a family of structure-preserving geometric algorithms for guiding center dynamics with long-term stability and accuracy. This discovery overcomes the difficulty associated with the unstable parasitic modes for variational symplectic integrators when applied to the degenerate guiding center Lagrangian. It is a pleasant surprise that Pauli's Hamiltonian for electrons, which predated the Dirac equation and marks the beginning of particle physics, reappears in classical physics as an effective algorithm for solving an important plasma physics problem. This technique is applicable to other degenerate Lagrangians reduced from regular Lagrangians.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

The geometric theory of charge conservation in particle-in-cell simulations

In recent years, several gauge-symmetric particle-in-cell (PIC) methods have been developed whose simulations of particles and electromagnetic fields exactly conserve charge. While it is rightly observed that these methods’ gauge symmetry gives rise to their charge conservation, this causal relationship has generally been asserted via ad hoc derivations of the associated conservation laws. In this work, we develop a comprehensive theoretical grounding for charge conservation in gauge-symmetric Lagrangian and Hamiltonian PIC algorithms. For Lagrangian variational PIC methods, we apply Noether’s second theorem to demonstrate that gauge symmetry gives rise to a local charge conservation law as an off-shell identity. For Hamiltonian splitting methods, we show that the momentum map establishes their charge conservation laws. We define a new class of algorithms – gauge-compatible splitting methods – that exactly preserve the momentum map associated with a Hamiltonian system’s gauge symmetry – even after time discretization. This class of algorithms affords splitting schemes a decided advantage over alternative Hamiltonian integrators. We apply this general technique to design a novel, explicit, symplectic, gauge-compatible splitting PIC method, whose momentum map yields an exact local charge conservation law. Finally, our study clarifies the appropriate initial conditions for such schemes and examines their symplectic reduction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Rotational bands beyond the Elliott model

Rotational bands are commonplace in the spectra of atomic nuclei. Inspired by early descriptions of these bands by quadrupole deformations of a liquid drop, Elliott constructed discrete nucleon representations of SU(3) from fermionic creation and annihilation operators. Ever since, Elliott's model has been foundational to descriptions of rotation in nuclei. Later work, however, suggested the symplectic extension Sp(3, R) provides a more unified picture. We decompose no-core shell-model nuclear wave functions into symmetry-defined subspaces for several beryllium isotopes, as well as 20 Ne, using the quadratic Casimirs of both Elliott's SU(3) and Sp(3, R). The band structure, delineated by strong B(E2) values, has a more consistent description in Sp(3, R) rather than SU(3). In particular, we confirm previous work finding in some nuclides strongly connected upper and lower bands with the same underlying symplectic structure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

GFiRe—Gauge Field integrator for Reheating

We present a new numerical algorithm and code, GFiRe, for solving the nonlinear evolution of Abelian gauge fields coupled to complex scalar fields in homogeneous and isotropic spacetimes. We adopt a hybrid approach to solving the system: the spatial derivatives are discretized using standard Lattice Gauge Field Theory techniques, whereas the time evolution of the fields and scale factor is implemented with explicit, composite, symplectic integrators. An important property of our compound algorithm is that the discretized Gauss constraint is respected exactly, regardless of the order of the symplectic integrator. This remains true even when the background expansion is computed “self-consistently”; that is, when the expansion history is computed using spatial averaged components of the energy momentum tensor in the simulation volume. Hence, our code can also be used in cases where the fields dominate the energy density of the universe, for example, during reheating after inflation. We test the algorithm in scenarios of reheating where the inflaton is a complex scalar field with a potential ∝(2|Φ| 2 - v 2 ) 2 and is coupled to an Abelian gauge field. Tracing the evolution of the system through complex dynamics (including resonant excitation of fields, backreaction, formation of solitons, and changes in the equation of state) in a self-consistently expanding universe, we find the energy conservation violation (<10 -4 ) to be very stable and the Gauss constraint violation (<10 -6 ) to be dominated by differencing noise.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Explicit structure-preserving geometric particle-in-cell algorithm in curvilinear orthogonal coordinate systems and its applications to whole-device 6D kinetic simulations of tokamak physics

Explicit structure-preserving geometric particle-in-cell (PIC) algorithm in curvilinear orthogonal coordinate systems is developed. The work reported represents a further development of the structure-preserving geometric PIC algorithm achieving the goal of practical applications in magnetic fusion research. The algorithm is constructed by discretizing the field theory for the system of charged particles and electromagnetic field using Whitney forms, discrete exterior calculus, and explicit non-canonical symplectic integration. In addition to the truncated infinitely dimensional symplectic structure, the algorithm preserves exactly many important physical symmetries and conservation laws, such as local energy conservation, gauge symmetry and the corresponding local charge conservation. As a result, the algorithm possesses the long-term accuracy and fidelity required for first-principles-based simulations of the multiscale tokamak physics. The algorithm has been implemented in the SymPIC code, which is designed for high-efficiency massively-parallel PIC simulations in modern clusters. The code has been applied to carry out whole-device 6D kinetic simulation studies of tokamak physics. A self-consistent kinetic steady state for fusion plasma in the tokamak geometry is numerically found with a predominately diagonal and anisotropic pressure tensor. The state also admits a steady-state sub-sonic ion flow in the range of 10 km s -1 , agreeing with experimental observations and analytical calculations Kinetic ballooning instability in the self-consistent kinetic steady state is simulated. It is shown that high-n ballooning modes have larger growth rates than low-n global modes, and in the nonlinear phase the modes saturate approximately in 5 ion transit times at the 2% level by the E × B flow generated by the instability. These results are consistent with early and recent electromagnetic gyrokinetic simulations.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Nonlinear photoconductivities and quantum geometry of chiral multifold fermions

Here, chiral multifold fermions are quasiparticles that appear only in chiral crystals such as transition metal silicides in the cubic B20 structure (i.e., the CoSi family), and they may show exotic physical properties. Here we study the injection and shift photoconductivities and also the related geometrical quantities for several types of chiral multifold fermions, including spin-1/2 as well as pseudospin-1 and -3/2 fermions, dubbed as Kramers Weyl, triple-point, and Rarita-Schwinger-Weyl (RSW) fermions, respectively. We utilize the minimal symmorphic model to describe the triple-point fermions (TPF). We also consider the more realistic model Hamiltonian for the CoSi family including both linear and quadratic terms. We find that injection currents due to circularly polarized light are quantized as a result of the Chern numbers carried by the multifold fermions within the linear models. Surprisingly, we discover that in the TPF model, the linear shift conductivities, responsible for the shift current generation by linearly polarized light, are proportional to the pseudo spin-orbit coupling and independent of photon frequency. In contrast, for the RSW and Kramer Weyl fermions, the linear shift conductivity is linearly proportional to photon frequency. The numerical results agree with the power-counting analysis for quadratic Hamiltonians. The frequency independence of the linear shift conductivity could be attributed to the strong resonant symplectic Christoffel symbols of the flat bands. Moreover, the calculated symplectic Christoffel symbols show significant peaks at the nodes, suggesting that the shift currents are due to the strong geometrical response near the topological nodes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Topological quantum computation on a chiral Kondo chain

We describe the chiral Kondo chain model based on the symplectic Kondo effect and demonstrate that it has a quantum critical ground state populated by non-Abelian anyons. We show that the fusion channel of two arbitrary anyons can be detected by locally coupling the two anyons to an extra single channel of chiral current and measuring the corresponding conductance at a finite frequency. Based on such measurements, we propose that the chiral Kondo chain model with symplectic symmetry can be used for the implementation of measurement-only topological quantum computations, and it possesses several distinct features favorable for such applications. Furthermore, the sources and effects of errors in the proposed system are analyzed and possible material realizations are discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

A general framework for gravitational charges and holographic renormalization

We develop a general framework for constructing charges associated with diffeomorphisms in gravitational theories using covariant phase space techniques. This framework encompasses both localized charges associated with space–time subregions, as well as global conserved charges of the full space–time. Expressions for the charges include contributions from the boundary and corner terms in the subregion action, and are rendered unambiguous by appealing to the variational principle for the subregion, which selects a preferred form of the symplectic flux through the boundaries. The Poisson brackets of the charges on the subregion phase space are shown to reproduce the bracket of Barnich and Troessaert for open subsystems, thereby giving a novel derivation of this bracket from first principles. In the context of asymptotic boundaries, we show that the procedure of holographic renormalization can be always applied to obtain finite charges and fluxes once suitable counterterms have been found to ensure a finite action. This enables the study of larger asymptotic symmetry groups by loosening the boundary conditions imposed at infinity. We further present an algorithm for explicitly computing the counterterms that renormalize the action and symplectic potential, and, as an application of our framework, demonstrate that it reproduces known expressions for the charges of the generalized Bondi–Metzner–Sachs algebra.

Physics↗
DOE OSTI2024

Asymptotic structure of higher dimensional Yang-Mills theory

Using the covariant phase space formalism, we construct the phase space for non-Abelian gauge theories in (d+2)-dimensional Minkowski spacetime for any d ≥ 2, including the edge modes that symplectically pair to the low energy degrees of freedom of the gauge field. Despite the fact that the symplectic form in odd and even-dimensional spacetimes appear ostensibly different, we demonstrate that both cases can be treated in a unified manner by utilizing the shadow transform. Upon quantization, we recover the algebra of the vacuum sector of the Hilbert space and derive a Ward identity that implies the leading soft gluon theorem in (d+2)-dimensional spacetime.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Modeling Transverse Space Charge effects in IOTA with pyORBIT

The role and mitigation of space charge effects are important aspects of the beam physics research to be performed at the Integrable Optics Test Accelerator (IOTA) at Fermilab. The impact of nonlinear integrability (partial and complete) on space charge driven incoherent and coherent resonances will be a new feature of this accelerator and has the potential to influence the design of future low energy proton synchrotrons. In this report we model transverse space charge effects using the PIC code pyORBIT. First we benchmark the single particle tracking results against MADX with checks on symplecticity, tune footprints, and dynamic aperture in a partially integrable lattice realized with a special octupole string insert. Our space charge calculations begin with an examination of the 4D symplecticity. Short term tracking is done first with the initial bare lattice and then with a matching of the rms values with space charge. Next, we explore slow initialization of charge as a technique to establish steady state and reduce emittance growth and beam loss following injection into a mismatched lattice. We establish values of space charge simulation parameters so as to ensure numerical convergence. Finally, we compare the simulated space charge induced tune shifts and footprints against theory.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Optimized Lie–Trotter–Suzuki decompositions for two and three non-commuting terms

Lie–Trotter–Suzuki decompositions are an efficient way to approximate operator exponentials exp ( t H ) when H is a sum of n (non-commuting) terms which, individually, can be exponentiated easily. They are employed in time-evolution algorithms for tensor network states, digital quantum simulation protocols, path integral methods like quantum Monte Carlo, and splitting methods for symplectic integrators in classical Hamiltonian systems. Here, we provide optimized decompositions up to order t 6 . The leading error term is expanded in nested commutators (Hall bases) and we minimize the 1-norm of the coefficients. For n = 2 terms, several of the optima we find are close to those in McLachlan (1995). Generally, our results substantially improve over unoptimized decompositions by Forest, Ruth, Yoshida, and Suzuki. We explain why these decompositions are sufficient to efficiently simulate any one- or two-dimensional lattice model with finite-range interactions. This follows by solving a partitioning problem for the interaction graph.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Generalization of QCD3 symmetry-breaking and flavored quiver dualities

We extend the recently proposed symmetry breaking scenario of QCD 3 to the so-called “master” (2 + 1)d bosonization duality, which has bosonic and fermionic matter on both ends. Using anomaly arguments, a phase diagram emerges with several novel regions. We then construct 2+1 dimensional dualities for flavored quivers using node-by-node dualization. Such dualities are applicable to theories which live on domain walls in QCD 4 -like theories with dynamical quarks. We also derive dualities for quivers based on orthogonal and symplectic gauge groups. Lastly, we support the conjectured dualities using holographic constructions, even though several aspects of this holographic construction remain mostly qualitative.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Trialities of minimally supersymmetric 2d gauge theories

We study dynamics of two-dimensional N = (0, 1) supersymmetric gauge theories. In particular, we propose that there is an infrared triality between certain triples of theories with orthogonal and symplectic gauge groups. The proposal is supported by matching of anomalies and elliptic genera. This triality can be viewed as a (0, 1) counterpart of the (0, 2) triality proposed earlier by two of the authors and A. Gadde. We also describe the relation between global anomalies in gauge theoretic and sigma-model descriptions, filling in a gap in the present literature.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

The soft $ \mathcal{S} $-matrix in gravity

The gravitational $ \mathcal{S}$-matrix defined with an infrared (IR) cutoff factorizes into hard and soft factors. The soft factor is universal and contains all the IR and collinear divergences. Here we show, in a momentum space basis, that the intricate expression for the soft factor is fully reproduced by two boundary currents, which live on the celestial sphere. The first of these is the supertranslation current, which generates spacetime supertranslations. The second is its symplectic partner, the Goldstone current for spontaneously broken supertranslations. The current algebra has an off-diagonal level structure involving the gravitational cusp anomalous dimension and the logarithm of the IR cutoff. It is further shown that the gravitational memory effect is contained as an IR safe observable within the soft $ \mathcal{S}$-matrix.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗