Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “invariance”

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

3d-3d correspondence for mapping tori

One of the main challenges in 3d-3d correspondence is that no existent approach offers a complete description of 3d N = 2 SCFT T [ M 3 ] — or, rather, a “collection of SCFTs” as we refer to it in the paper — for all types of 3-manifolds that include, for example, a 3-torus, Brieskorn spheres, and hyperbolic surgeries on knots. The goal of this paper is to overcome this challenge by a more systematic study of 3d-3d correspondence that, first of all, does not rely heavily on any geometric structure on M 3 and, secondly, is not limited to a particular supersymmetric partition function of T [ M 3 ]. In particular, we propose to describe such “collection of SCFTs” in terms of 3d N = 2 gauge theories with “non-linear matter” fields valued in complex group manifolds. As a result, we are able to recover familiar 3-manifold invariants, such as Turaev torsion and WRT invariants, from twisted indices and half-indices of T [ M 3 ], and propose new tools to compute more recent q -series invariants Z * ( M 3 ) in the case of manifolds with b 1 > 0. Although we use genus-1 mapping tori as our “case study,” many results and techniques readily apply to more general 3-manifolds, as we illustrate throughout the paper.

Topological Field Theories↗
DOE OSTI2020

Interface conformal anomalies

We consider two d ≥ 2 conformal field theories (CFTs) glued together along a codimension one conformal interface. The conformal anomaly of such a system contains both bulk and interface contributions. In a curved-space setup, we compute the heat kernel coefficients and interface central charges in free theories. The results are consistent with the known boundary CFT data via the folding trick. In d = 4, two interface invariants generally allowed as anomalies turn out to have vanishing interface charges. These missing invariants are constructed from components with odd parity with respect to flipping the orientation of the defect. We conjecture that all invariants constructed from components with odd parity may have vanishing coefficient for symmetric interfaces, even in the case of interacting interface CFT.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Celestial amplitudes as AdS-Witten diagrams

Both celestial and momentum space amplitudes in four dimensions are beset by divergences resulting from spacetime translation and sometimes scale invariance. In this paper we consider a (linearized) marginal deformation of the celestial CFT for Yang-Mills theory which preserves 2D conformal invariance but breaks both spacetime translation and scale invariance and involves a chirally coupled massive scalar. The resulting MHV celestial amplitudes are completely finite (apart from the usual soft and collinear divergences and isolated poles in the sum of the weights) and take the canonical CFT form. Moreover, we show they can be simply rewritten in terms of AdS 3 -Witten contact diagrams which evaluate to the well-known D-functions, thereby forging a direct connection between flat and AdS holography.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Nearly Periodic Maps and Geometric Integration of Noncanonical Hamiltonian Systems

Abstract M. Kruskal showed that each continuous-time nearly periodic dynamical system admits a formal U (1)-symmetry, generated by the so-called roto-rate. When the nearly periodic system is also Hamiltonian, Noether’s theorem implies the existence of a corresponding adiabatic invariant. We develop a discrete-time analog of Kruskal’s theory. Nearly periodic maps are defined as parameter-dependent diffeomorphisms that limit to rotations along a U (1)-action. When the limiting rotation is non-resonant, these maps admit formal U (1)-symmetries to all orders in perturbation theory. For Hamiltonian nearly periodic maps on exact presymplectic manifolds, we prove that the formal U (1)-symmetry gives rise to a discrete-time adiabatic invariant using a discrete-time extension of Noether’s theorem. When the unperturbed U (1)-orbits are contractible, we also find a discrete-time adiabatic invariant for mappings that are merely presymplectic, rather than Hamiltonian. As an application of the theory, we use it to develop a novel technique for geometric integration of non-canonical Hamiltonian systems on exact symplectic manifolds.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Effect of Heave Plate Hydrodynamic Force Parameterization on a Two-Body Wave Energy Converter

Heave plates are one approach to generating the reaction force necessary to harvest energy from ocean waves. In a Morison equation description of the hydrodynamic force, the components of drag and added mass depend primarily on the heave plate oscillation. These terms may be parameterized in three ways: (1) as a single coefficient invariant across sea state, most accurate at the reference sea state, (2) coefficients dependent on the oscillation amplitude, but invariant in phase, that are most accurate for relatively small amplitude motions, and (3) coefficients dependent on both oscillation amplitude and phase, which are accurate for all oscillation amplitudes. We validate a MATLAB model for a two-body point absorber wave energy converter against field data and a dynamical model constructed in ProteusDS. We then use the MATLAB model to evaluate the effect of these parameterizations on estimates of heave plate motion, tension between the float and heave plate, and wave energy converter electrical power output. We find that power predictions using amplitude-dependent coefficients differ by up to 30% from models using invariant coefficients for regular waves ranging in height from 0.5 to 1.9 m. Amplitude- and phase-dependent coefficients, however, yield less than a 5% change when compared with coefficients dependent on amplitude only. This suggests that amplitude-dependent coefficients can be important for accurate wave energy converter modeling, but the added complexity of phase-dependent coefficients yields little further benefit. We show similar, though less pronounced, trends in maximum tether tension, but note that heave plate motion has only a weak dependence on coefficient fidelity. Finally, we emphasize the importance of using experimentally derived added mass over that calculated from boundary element methods, which can lead to substantial under-prediction of power output and peak tether tension.

dynamical model↗
DOE OSTI2022

Krylov subspace recycling for evolving structures

Krylov subspace recycling is a powerful tool when solving a long series of large, sparse linear systems that change only slowly over time. In PDE constrained shape optimization, these series appear naturally, as typically hundreds or thousands of optimization steps are needed with only small changes in the geometry. In this setting, however, applying Krylov subspace recycling can be a difficult task. As the geometry evolves, in general, so does the finite element mesh defined on or representing this geometry, including the numbers of nodes and elements and element connectivity. This is especially the case if re-meshing techniques are used. As a result, the number of algebraic degrees of freedom in the system changes, and in general the linear system matrices resulting from the finite element discretization change size from one optimization step to the next. Changes in the mesh connectivity also lead to structural changes in the matrices. In the case of re-meshing, even if the geometry changes only a little, the corresponding mesh might differ substantially from the previous one. Obviously, this prevents any straightforward mapping of the approximate invariant subspace of the linear system matrix (the focus of recycling in this work) from one optimization step to the next; similar problems arise for other selected subspaces. In this paper, we present an algorithm to map an approximate invariant subspace of the linear system matrix for the previous optimization step to an approximate invariant subspace of the linear system matrix for the current optimization step, for general meshes. This is achieved by exploiting the map from coefficient vectors to finite element functions on the mesh, combined with interpolation or approximation of functions on the finite element mesh. We demonstrate the effectiveness of our approach numerically with several proof of concept studies for a specific meshing technique.

42 ENGINEERING↗
DOE OSTI2025

Dynamics of McMillan mappings I. McMillan multipoles

In this article, we consider two dynamical systems: the McMillan sextupole and octupole integrable mappings, originally proposed by Edwin McMillan. Both represent the simplest symmetric McMillan maps, characterized by a single intrinsic parameter. While these systems find numerous applications across various domains of mathematics and physics, some of their dynamical properties remain unexplored. We aim to bridge this gap by providing a comprehensive description of all stable trajectories, including the parametrization of invariant curves, Poincaré rotation numbers, and canonical action–angle variables. In the second part, we establish connections between these maps and general chaotic maps in standard form. Our investigation reveals that the McMillan sextupole and octupole serve as first-order approximations of the dynamics around the fixed point, akin to the linear map and quadratic invariant (known as the Courant–Snyder invariant in accelerator physics), which represents zeroth-order approximations (referred to as linearization). Furthermore, we propose a novel formalism for nonlinear Twiss parameters, which accounts for the dependence of rotation number on amplitude. This stands in contrast to conventional betatron phase advance used in accelerator physics, which remains independent of amplitude. Notably, in the context of accelerator physics, this new formalism demonstrates its capability in predicting dynamical aperture around low-order resonances for flat beams, a critical aspect in beam injection/extraction scenarios.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

High precision tests of QCD without scale or scheme ambiguities: The 40th anniversary of the Brodsky–Lepage–Mackenzie method

A key issue in making precise predictions in QCD is the uncertainty in setting the renormalization scale μ r and thus determining the correct values of the QCD running coupling α s (μ r ) at each order in the perturbative expansion of a QCD observable. It has often been conventional to simply set the renormalization scale to the typical scale of the process Q and vary it in the range μ r $\in$ [Q/2, 2Q] in order to estimate the theoretical error. This is the practice of Conventional Scale Setting (CSS). The resulting CSS prediction will however depend on the theorist’s choice of renormalization scheme and the resulting pQCD series will diverge factorially. It will also disagree with renormalization scale setting used in QED and electroweak theory thus precluding grand unification. A solution to the renormalization scale-setting problem is offered by the Principle of Maximum Conformality (PMC), which provides a systematic way to eliminate the renormalization scale-and-scheme dependence in perturbative calculations. The PMC method has rigorous theoretical foundations, it satisfies Renormalization Group Invariance (RGI) and preserves all self-consistency conditions derived from the renormalization group. The PMC cancels the renormalon growth, reduces to the Gell-Mann–Low scheme in the N c → 0 Abelian limit and leads to scale- and scheme-invariant results. The PMC has now been successfully applied to many high-energy processes. In this article we summarize recent developments and results in solving the renormalization scale and scheme ambiguities in perturbative QCD. In particular, we present a recently developed method the PMC ∞ and its applications, comparing the results with CSS. The method preserves the property of renormalizable SU(N)/U(1) gauge theories defined as Intrinsic Conformality (iCF). This property underlies the scale invariance of physical observables and leads to a remarkably efficient method to solve the conventional renormalization scale ambiguity at every order in pQCD. This new method reflects the underlying conformal properties displayed by pQCD at NNLO, eliminates the scheme dependence of pQCD predictions and is consistent with the general properties of the PMC. A new method to identify conformal and β-terms, which can be applied either to numerical or to theoretical calculations is also shown. We present results for the thrust and C-parameter distributions in e + e - annihilation showing errors and comparison with the CSS. We also show results for a recent innovative comparison between the CSS and the PMC ∞ applied to the thrust distribution investigating both the QCD conformal window and the QED N c → 0 limit. In order to determine the thrust distribution along the entire renormalization group flow from the highest energies to zero energy, we consider the number of flavors near the upper boundary of the conformal window. In this flavor-number regime the theory develops a perturbative infrared interacting fixed point. These results show that PMC ∞ leads to higher precision and introduces new interesting features in the PMC. In fact, this method preserves with continuity the position of the peak, showing perfect agreement with the experimental data already at NNLO. We also show a detailed comparison of the PMC ∞ with the other PMC approaches: the multi-scale-setting approach (PMCm) and the single-scale-setting approach (PMCs) by comparing their predictions for three important fully integrated quantities R e+e- , R$_Τ$ and Γ (H→$b\bar{b}$) up to the four-loop accuracy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Modeling the Effects of Drift Orbit Bifurcation on Radiation Belt Electrons

Abstract Drift orbit bifurcation (DOB) has been suggested to play a major role in the loss and transport of radiation belt electrons since it violates the second adiabatic invariants of particles and makes the third invariant undefined. Results from our guiding‐center test particle simulations using the Tsyganenko‐1989c magnetic field model show that the DOB could affect a broad region of the outer radiation belt, which can penetrate inside the geosynchronous orbit at Kp ≥ 3, and its effects are more significant further away from the Earth, at a higher Kp and for higher electron energies. Specifically, the short‐term simulation results after one electron drift show both traditional and nontraditional DOB transport of electrons, with the nontraditional DOB, caused by a third minimum of the magnetic field strength near the equator, reported here for the first time. Moreover, our results show large ballistic jumps in the second invariant and radial distance for electrons at high equatorial pitch angles after one drift. In addition, long‐term DOB transport coefficients of electrons over many drifts are calculated based on our simulation results. We find that the pitch angle and radial diffusion coefficients of electrons due to DOB could be comparable to or even larger than those caused by electron interactions with chorus and Ultra‐Low‐Frequency waves, respectively. In sum, our results demonstrate that DOB could cause effective loss and transport of radiation belt electrons even in the absence of waves.

Astronomy & Astrophysics↗
DOE OSTI2021

Time evolution of ML-MCTDH wavefunctions. I. Gauge conditions, basis functions, and singularities

We derive a family of equations-of-motion (EOMs) for evolving multi-layer multiconfiguration time-dependent Hartree (ML-MCTDH) wavefunctions that, unlike the standard ML-MCTDH EOMs, never require the evaluation of the inverse of singular matrices. All members of this family of EOMs make use of alternative static gauge conditions than those used for standard ML-MCTDH. These alternative conditions result in an expansion of the wavefunction in terms of a set of potentially arbitrary orthonormal functions, rather than in terms of a set of non-orthonormal and potentially linearly dependent functions, as is the case for standard ML-MCTDH. We show that the EOMs used in the projector splitting integrator (PSI) and the invariant EOM approaches are two special cases of this family obtained from different choices for the dynamic gauge condition, with the invariant EOMs making use of a choice that introduces potentially unbounded operators into the EOMs. As a consequence, all arguments for the existence of parallelizable integration schemes for the invariant EOMs can also be applied to the PSI EOMs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Lightweight and effective tensor sensitivity for atomistic neural networks

Atomistic machine learning focuses on the creation of models that obey fundamental symmetries of atomistic configurations, such as permutation, translation, and rotation invariances. In many of these schemes, translation and rotation invariance are achieved by building on scalar invariants, e.g., distances between atom pairs. There is growing interest in molecular representations that work internally with higher rank rotational tensors, e.g., vector displacements between atoms, and tensor products thereof. Here, we present a framework for extending the Hierarchically Interacting Particle Neural Network (HIP-NN) with Tensor Sensitivity information (HIP-NN-TS) from each local atomic environment. Crucially, the method employs a weight tying strategy that allows direct incorporation of many-body information while adding very few model parameters. We show that HIP-NN-TS is more accurate than HIP-NN, with negligible increase in parameter count, for several datasets and network sizes. As the dataset becomes more complex, tensor sensitivities provide greater improvements to model accuracy. In particular, HIP-NN-TS achieves a record mean absolute error of 0.927 $\frac{\text{kcal}}{\text{mol}}$ for conformational energy variation on the challenging COMP6 benchmark, which includes a broad set of organic molecules. We also compare the computational performance of HIP-NN-TS to HIP-NN and other models in the literature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Optical N-plasmon: topological hydrodynamic excitations in graphene from repulsive Hall viscosity

Abstract Edge states occurring in Chern and quantum spin-Hall phases are signatures of the topological electronic band structure in two-dimensional (2D) materials. Recently, a new topological electromagnetic phase of graphene characterized by the optical N-invariant was proposed. Optical N-invariant arises from repulsive Hall viscosity in hydrodynamic many-body electron systems, distinct from the Chern and Z 2 invariants. In this paper, we introduce the topologically protected edge excitation—optical N-plasmon of interacting many-body electron systems in the topological optical N-phase. These optical N-plasmons are signatures of the topological plasmonic band structure in 2D materials. We demonstrate that optical N-plasmons exhibit unique dispersion relations, stability against various boundary conditions, and edge profiles when compared with the topologically trivial edge magneto plasmons. Based on the optical N-plasmon, we design an ultra sub-wavelength broadband topological hydrodynamic circulator, which is a chiral quantum radio-frequency circuit component crucial for information routing and interfacing quantum–classical computing systems. Furthermore, we reveal that optical N-plasmons can be effectively tuned by the neighboring dielectric environment without breaking the topological properties. Our work provides a smoking gun signature of topological electromagnetic phases occurring in 2D materials arising from repulsive Hall viscosity.

Physics↗
DOE OSTI2021

Asymptotic reversibility of thermal operations for interacting quantum spin systems via generalized quantum Stein’s lemma

Abstract For quantum spin systems in any spatial dimension with a local, translation-invariant Hamiltonian, we prove that asymptotic state convertibility from a quantum state to another one by a thermodynamically feasible class of quantum dynamics, called thermal operations, is completely characterized by the Kullback–Leibler (KL) divergence rate, if the state is translation-invariant and spatially ergodic. Our proof consists of two parts and is phrased in terms of a branch of the quantum information theory called the resource theory. First, we prove that any states, for which the min and max Rényi divergences collapse approximately to a single value, can be approximately reversibly converted into one another by thermal operations with the aid of a small source of quantum coherence. Second, we prove that these divergences collapse asymptotically to the KL divergence rate for any translation-invariant ergodic state. We show this via a generalization of the quantum Stein’s lemma for quantum hypothesis testing beyond independent and identically distributed situations. Our result implies that the KL divergence rate serves as a thermodynamic potential that provides a complete characterization of thermodynamic convertibility of ergodic states of quantum many-body systems in the thermodynamic limit, including out-of-equilibrium and fully quantum situations.

Physics↗
DOE OSTI2023

Obstructions to gapped phases from noninvertible symmetries

Quantum systems in 3+1 dimensions that are invariant under gauging a one-form symmetry enjoy novel noninvertible duality symmetries encoded by topological defects. These symmetries are renormalization group invariants which constrain dynamics. We show that such noninvertible symmetries often forbid a symmetry-preserving vacuum state with a gapped spectrum. In particular, we prove that a self-dual theory with $\mathbb{Z}^{(1)}_N$ one-form symmetry is gapless or spontaneously breaks the self-duality symmetry unless N = k 2 ⁢ℓ where –1 is a quadratic residue modulo . We also extend these results to noninvertible symmetries arising from invariance under more general gauging operations including, e.g., triality symmetries. Along the way, we discover how duality defects in symmetry-protected topological phases have a hidden time-reversal symmetry that organizes their basic properties. These noninvertible symmetries are realized in lattice gauge theories, which serve to illustrate our results.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Thermal transfer enhancement by hydrodynamic plasmons in electron bilayers

Here, we develop a theory of heat transfer induced by thermal charge fluctuations in two-dimensional electron double layers. We consider pristine systems comprised of identical layers and focus on the regime of sufficiently high temperatures and interlayer distances 𝑑, where the relevant charge fluctuations may be described using the hydrodynamic approach. In this limit heat transfer is dominated by the plasmon resonances. For systems with Galilean-invariant electron dispersion the interlayer thermal conductance 𝜘 is proportional to the kinematic viscosity of the electron liquid and decreases as 1/𝑑 4 . In the absence of Galilean invariance 𝜘∝𝜎/𝑑 3 , where 𝜎 is the intrinsic conductivity of the liquid. This strong enhancement can be traced to a drastically different broadening of plasmon resonances in systems with and without Galilean invariance.

2-dimensional systems↗
DOE OSTI2024

Characterization and classification of interacting ( 2 + 1 )-dimensional topological crystalline insulators with orientation-preserving wallpaper groups

While free fermion topological crystalline insulators have been largely classified, the analogous problem in the strongly interacting case has been only partially solved. In this work, we develop a characterization and classification of interacting, invertible fermionic topological phases in (2+1) dimensions with charge conservation, discrete magnetic translation and M-fold point group rotation symmetries, which form the group G f = U(1) f × Φ [Z 2 $\rtimes$Z M ] for M = 1,2,3,4, and 6. Φ is the magnetic flux per unit cell. We derive a topological response theory in terms of background crystalline gauge fields, which gives a complete classification of different phases and a physical characterization in terms of quantized response to symmetry defects. We then derive the same classification in terms of a set of real space invariants {$Θ^±_o$} that can be obtained from ground state expectation values of suitable partial rotation operators. We explicitly relate these real space invariants to the quantized coefficients in the topological response theory, and find the dependence of the invariants on the chiral central charge c – of the invertible phase. Finally, when Φ = 0 we derive an explicit map between the free and interacting classifications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Lorentz symmetry in ghost-free massive gravity

The role of Lorentz symmetry in ghost-free massive gravity is studied, emphasizing features emerging in approximately Minkowski spacetime. The static extrema and saddle points of the potential are determined and their Lorentz properties identified. Solutions preserving Lorentz invariance and ones breaking four of the six Lorentz generators are constructed. Locally, globally, and absolutely stable Lorentz-invariant extrema are found to exist for certain parameter ranges of the potential. Gravitational waves in the linearized theory are investigated. Deviations of the fiducial metric from the Minkowski metric are shown to lead to pentarefringence of the five wave polarizations, which can include superluminal modes and subluminal modes with negative energies in certain observer frames. The Newton limit of ghost-free massive gravity is explored. The propagator is constructed and used to obtain the gravitational potential energy between two point masses. Here, the result extends the Fierz-Pauli limit to include corrections generically breaking both rotation and boost invariance.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Finding evidence for inflation and the origin of galactic magnetic fields with CMB surveys

The origin of the μG magnetic fields observed in galaxies is unknown. One promising scenario is that magnetic fields generated during inflation, larger than 0.1 nG on Mpc scales, were adiabatically compressed to μG strengths in galaxies during structure formation. Thus, detecting a scale-invariant primordial magnetic field (PMF) above 0.1 nG on Mpc scales just after recombination would indicate an inflationary origin of galactic magnetic fields. Furthermore, this would also provide compelling evidence that inflation occurred since only an inflationary mechanism could generate such a strong, scale-invariant magnetic field on Mpc scales. In contrast, constraining the scale-invariant PMF strength to be below 0.1 nG would imply an inflationary scenario is not the primary origin, since such weak PMFs cannot be amplified enough via adiabatic compression to produce the strength of the galactic fields we observe today. We find that measurements of anisotropic birefringence by future CMB surveys will be able to improve the sensitivity to Mpc-scale inflationary PMFs by an order of magnitude, and, in particular, that CMB-HD would lower the upper bound to 0.072 nG at the 95% C.L., which is below the critical 0.1 nG threshold for ruling out a purely inflationary origin. If inflationary PMFs exist, we find that a CMB-HD survey would be able to detect them with about 3σ significance or higher, providing evidence for inflation itself.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗