Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CHARGE CONSERVATION”

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

Integrable higher-spin deformations of sigma models from auxiliary fields

We construct a new infinite family of integrable deformations of the principal chiral model (PCM) parametrized by an interaction function of several variables, which extends the formalism of [C. Ferko and L. Smith, An infinite family of integrable sigma models using auxiliary fields, .] and includes deformations of the PCM by functions of both the stress tensor and higher-spin conserved currents. We show in detail that every model in this class admits a Lax representation for its equations of motion, and that the Poisson bracket of the Lax connection takes the Maillet form, establishing the existence of an infinite set of Poisson-commuting conserved charges. We argue that the non-Abelian T-dual of any model in this family is classically integrable, and that T-duality “commutes” with a general deformation in this class, in a sense which we make precise. Finally, we demonstrate that these higher-spin auxiliary field deformations can be extended to accommodate the addition of a Wess-Zumino term, and we exhibit the Lax connection in this case. Published by the American Physical Society 2025

Bielli, Daniele (ORCID:0009000640034768)↗
DOE OSTI2017

Constraining the hadronic spectrum through QCD thermodynamics on the lattice

Fluctuations of conserved charges allow us to study the chemical composition of hadronic matter. A comparison between lattice simulations and the hadron resonance gas (HRG) model suggested the existence of missing strange resonances. To clarify this issue we calculate the partial pressures of mesons and baryons with different strangeness quantum numbers using lattice simulations in the confined phase of QCD. In order to make this calculation feasible, we perform simulations at imaginary strangeness chemical potentials. We systematically study the effect of different hadronic spectra on thermodynamic observables in the HRG model and compare to lattice QCD results. We show that, for each hadronic sector, the well-established states are not enough in order to have agreement with the lattice results. Additional states, either listed in the Particle Data Group booklet (PDG) but not well established, or predicted by the quark model (QM), are necessary in order to reproduce the lattice data. For mesons, it appears that the PDG and the quark model do not list enough strange mesons, or that, in this sector, interactions beyond those included in the HRG model are needed to reproduce the lattice QCD results.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Large-Scale Description of Interacting One-Dimensional Bose Gases: Generalized Hydrodynamics Supersedes Conventional Hydrodynamics

The theory of generalized hydrodynamics (GHD) was recently developed as a new tool for the study of inhomogeneous time evolution in many-body interacting systems with infinitely many conserved charges. In this Letter, we show that it supersedes the widely used conventional hydrodynamics (CHD) of one-dimensional Bose gases. We illustrate this by studying “nonlinear sound waves” emanating from initial density accumulations in the Lieb-Liniger model. Here, we show that, at zero temperature and in the absence of shocks, GHD reduces to CHD, thus for the first time justifying its use from purely hydrodynamic principles. We show that sharp profiles, which appear in finite times in CHD, immediately dissolve into a higher hierarchy of reductions of GHD, with no sustained shock. CHD thereon fails to capture the correct hydrodynamics. We establish the correct hydrodynamic equations, which are finite-dimensional reductions of GHD characterized by multiple, disjoint Fermi seas. We further verify that at nonzero temperature, CHD fails at all nonzero times. Finally, we numerically confirm the emergence of hydrodynamics at zero temperature by comparing its predictions with a full quantum simulation performed using the NRG-TSA-abacus algorithm. The analysis is performed in the full interaction range, and is not restricted to either weak- or strong-repulsion regimes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Raman Response in the Nematic Phase of FeSe

Raman experiments on bulk FeSe revealed that the low-frequency part of the B 1g Raman response R B1 g (Ω), which probes nematic fluctuations, rapidly decreases below the nematic transition at T n ~ 85 K. Such behavior is expected when a gap opens up and at a first glance is inconsistent with the fact that FeSe remains a metal below Tn. We argue that the drop of R B1 g (Ω) can be ascribed to the fact that the nematic order drastically changes the orbital content of low-energy excitations near hole and electron pockets, making them nearly mono-orbital. In this situation, the B 1g Raman response gets reduced by the same vertex corrections that enforce charge conservation in the symmetric Raman channel. The reduction holds at low frequencies and gives rise to gaplike behavior of R B1 g (Ω). Here, we also show that the enhancement of the B 1g Raman response near T n is consistent with the sign change of the nematic order parameter between hole and electron pockets.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Pairwise Multiparticle States and the Monopole Unitarity Puzzle

We suggest a novel resolution for a decades old mystery-what happens when a positron scatters off a minimal grand-unification-theory monopole in an s wave, a puzzle first discussed by Callan in 1983. Using the language of on shell amplitudes and pairwise helicity we suggest that the final state contains two up quarks and a down quark in an entangled "pairwise" multiparticle state-the only particle final state that satisfies angular momentum and gauge charge conservation. The cross section for this process is as large as in the original Rubakov-Callan effect, only suppressed by the QCD scale. The final state we find cannot be seen in Callan's truncated 2D theory, since our new pairwise state appears only in more than two dimensions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Theory of Nonlinear Diffusion with a Physical Gapped Mode

In a system with one conserved charge the charge diffusion is modified by nonlinear self-interactions within an effective field theory (EFT) of diffusive fluctuations. We include the slowest ultraviolet (UV) mode, constructing a UV-regulated EFT. The relaxation time of this UV mode is protected from renormalization, as supported by experimental data in a bad metal system. Furthermore, the retarded density-density Green’s function acquires four branch points, eventually increasing the range of applicability. We discuss the fate of long time tails as well as implications for the quark gluon plasma. Published by the American Physical Society 2024

Physics↗
DOE OSTI2024

Infinite Family of Integrable Sigma Models Using Auxiliary Fields

We introduce a class of 2D sigma models which are parametrized by a function of one variable. In addition to the physical field g , these models include an auxiliary field v α which mediates interactions in a prescribed way. We prove that every theory in this family is classically integrable, in that it possesses an infinite set of conserved charges in involution, which can be constructed from a Lax representation for the equations of motion. This class includes the principal chiral model (PCM) and all deformations of the PCM by functions of the energy-momentum tensor. Published by the American Physical Society 2024

Physics↗
DOE OSTI2020

Effective spin-mixing conductance of topological-insulator/ferromagnet and heavy-metal/ferromagnet spin-orbit-coupled interfaces: A first-principles Floquet-nonequilibrium Green function approach

The spin-mixing conductance (SMC) is a key quantity determining efficiency of spin transport across interfaces. Thus, knowledge of its precise value is required for accurate measurement of parameters quantifying numerous effects in spintronics, such as spin-orbit torque, spin Hall magnetoresistance, spin Hall effect, and spin pumping. However, the standard expression for SMC, provided by the scattering theory in terms of the reflection probability amplitudes, is inapplicable when strong spin-orbit coupling (SOC) is present directly at the interface. This is the precisely the case of topological-insulator/ferromagnet and heavy-metal/ferromagnet interfaces of great contemporary interest. We introduce an approach where first-principles Hamiltonian of these interfaces, obtained from noncollinear density functional theory (ncDFT) calculations, is combined with charge-conserving Floquet-nonequilibrium-Green-function formalism to compute directly the pumped spin current $I^{S_z}_L$ into semi-infinite left lead of two-terminal heterostructures Cu/X/Co/Cu or Y/Co/Cu—where X = Bi 2 Se 3 and Y = Pt or W—due to microwave-driven steadily precessing magnetization of the Co layer. This allows us to extract an effective SMC as a prefactor in $I^{S_z}_L$ versus precession cone angle θ dependence, as long as it remains the same, $I^{S_z}$ ∝ sin 2 θ , as in the case where SOC is absent. By comparing calculations where SOC in switched off versus switched on in ncDFT calculations, we find that SOC consistently reduces the pumped spin current and therefore the effective SMC.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Emergent Dimer-Model Topological Order and Quasiparticle Excitations in Liquid Crystals: Combinatorial Vortex Lattices

Liquid crystals have proven to provide a versatile experimental and theoretical platform for studying topological objects such as vortices, skyrmions, and hopfions. In parallel, in hard condensed matter physics, the concept of topological phases and topological order has been introduced in the context of spin liquids to investigate emergent phenomena like quantum Hall effects and high-temperature superconductivity. Here, we bridge these two seemingly disparate perspectives on topology in physics. Combining experiments and simulations, we show how topological defects in liquid crystals can be used as versatile building blocks to create complex, highly degenerate topological phases, which we refer to as “combinatorial vortex lattices” (CVLs). CVLs exhibit extensive residual entropy and support locally stable quasiparticle excitations in the form of charge-conserving topological monopoles, which can act as mobile information carriers and be linked via Dirac strings. CVLs can be rewritten and reconfigured on demand, endowed with various symmetries, and modified through laser-induced topological surgery—an essential capability for information storage and retrieval. We demonstrate experimentally the realization, stability, and precise optical manipulation of CVLs, thus opening new avenues for understanding and technologically exploiting higher-hierarchy topology in liquid crystals and other ordered media.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Geometric origin of the energy-momentum tensor improvement terms

In a flat background, the canonical energy momentum tensor of Lorentz and conformally invariant matter field theories can be improved to a symmetric and traceless tensor that gives the same conserved charges. We argue that the geometric origin of this improvement process is unveiled when the matter theory is coupled to metric-affine gravity. In particular, we show that the Belinfante-Rosenfeld improvement terms correspond to the matter theory’s hypermomentum. The improvement terms in conformally invariant matter theories are also related to the hypermomentum; however, a general proof would require an extended investigation. We demonstrate our results through various examples, such as the free massless scalar, the Maxwell field, Abelian p-forms, the Dirac field, and a nonunitary massless scalar field. Possible applications of our method for theories that break Lorentz or special conformal invariance are briefly discussed.

classical solutions in field theory↗
DOE OSTI2015

Structures of the human Pals1 PDZ domain with and without ligand suggest gated access of Crb to the PDZ peptide-binding groove

This study characterizes the interaction between the carboxy-terminal (ERLI) motif of the essential polarity protein Crb and the Pals1/Stardust PDZ-domain protein. Structures of human Pals1 PDZ with and without a Crb peptide are described, explaining the highly conserved nature of the ERLI motif and revealing a sterically blocked peptide-binding groove in the absence of ligand. Many components of epithelial polarity protein complexes possess PDZ domains that are required for protein interaction and recruitment to the apical plasma membrane. Apical localization of the Crumbs (Crb) transmembrane protein requires a PDZ-mediated interaction with Pals1 (protein-associated with Lin7, Stardust, MPP5), a member of the p55 family of membrane-associated guanylate kinases (MAGUKs). This study describes the molecular interaction between the Crb carboxy-terminal motif (ERLI), which is required for Drosophila cell polarity, and the Pals1 PDZ domain using crystallography and fluorescence polarization. Only the last four Crb residues contribute to Pals1 PDZ-domain binding affinity, with specificity contributed by conserved charged interactions. Comparison of the Crb-bound Pals1 PDZ structure with an apo Pals1 structure reveals a key Phe side chain that gates access to the PDZ peptide-binding groove. Removal of this side chain enhances the binding affinity by more than fivefold, suggesting that access of Crb to Pals1 may be regulated by intradomain contacts or by protein–protein interaction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Exploring the QCD phase diagram through correlations and fluctuations

The exploration of the Quantum Chromodynamics (QCD) phase diagram is a central goal of relativistic heavy-ion collision experiments. This review focuses on the role of fluctuations and correlations as sensitive probes of the phase structure. We discuss theoretical advancements and experimental methodologies employed to map the QCD phase diagram, highlighting constraints derived from both lattice QCD calculations and existing experimental data. Key observables such as cumulants and factorial cumulants of conserved charges (e.g., net-proton, net-charge) are explored as promising signatures of phase transitions and the QCD critical point. We discuss how these quantities are measured experimentally and compared with theoretical predictions, addressing challenges and best practices for meaningful comparisons. Special attention is given to predictions and current experimental results at high baryon density, including recent findings from the STAR collaboration at RHIC. Finally, we identify open issues and future directions for fluctuation and correlation studies at lower collision energies, relevant for future measurements, for example by the CBM experiment.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Dynamics of hot QCD matter 2024 — Bulk properties

The second Hot QCD Matter 2024 conference at IIT Mandi focused on various ongoing topics in high-energy heavy-ion collisions, encompassing theoretical and experimental perspectives. This proceedings volume includes 19 contributions that collectively explore diverse aspects of the bulk properties of hot QCD matter. The topics encompass the dynamics of electromagnetic fields, transport properties, hadronic matter, spin hydrodynamics, and the role of conserved charges in high-energy environments. These studies significantly enhance our understanding of the complex dynamics of hot QCD matter, the quark–gluon plasma (QGP) formed in high-energy nuclear collisions. Advances in theoretical frameworks, including hydrodynamics, spin dynamics and fluctuation studies, aim to improve theoretical calculations and refine our knowledge of the thermodynamic properties of strongly interacting matter. Experimental efforts, such as those conducted by the ALICE and STAR collaborations, play a vital role in validating these theoretical predictions and deepening our insight into the QCD phase diagram, collectivity in small systems, and the early-stage behavior of strongly interacting matter. In conclusion, combining theoretical models with experimental observations offers a comprehensive understanding of the extreme conditions encountered in relativistic heavy-ion and proton-proton collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

ED-LPD, MATLAB standalone

This computer software is a MATLAB standalone code for implementation of the Lumped Particle Domain model. The model solves charge transfer kinetics, lithium diffusion dynamics, charge conservation and energy conservation occurring across the porous electrode pair unit cell sandwich to predict electrical and thermal response of a battery. Physico-chemical states of batteries are also evaluated by the model.

Kim, Gi-Heon↗
DOE OSTI2023

Hydrodynamics of higher-rank gauge theories

We extend recent work on hydrodynamics with global multipolar symmetries — known as “fracton hydrodynamics” — to systems in which the multipolar symmetries are gauged. We refer to the latter as “fracton magnetohydrodynamics”, in analogy to conventional magnetohydrodynamics (MHD), which governs systems with gauged charge conservation. We show that fracton MHD arises naturally from higher-rank Maxwell’s equations and in systems with one-form symmetries obeying certain constraints; while we focus on “minimal” higher-rank generalizations of MHD that realize diffusion, our methods may also be used to identify other, more exotic hydrodynamic theories (e.g., with magnetic subdiffusion). In contrast to semi-microscopic derivations of MHD, our approach elucidates the origin of the hydrodynamic modes by identifying the corresponding higher-form symmetries. Being rooted in symmetries, the hydrodynamic modes may persist even when the semi-microscopic equations no longer provide an accurate description of the system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Subsystem Rényi entropy of thermal ensembles for SYK-like models

The Sachdev-Ye-Kitaev model is an N N -modes fermionic model with infinite range random interactions. In this work, we study the thermal Rényi entropy for a subsystem of the SYK model using the path-integral formalism in the large- N N limit. The results are consistent with exact diagonalization and can be well approximated by thermal entropy with an effective temperature when subsystem size M\leq N/2 M ≤ N / 2 . We also consider generalizations of the SYK model with quadratic random hopping term or U(1) U ( 1 ) charge conservation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Development and Benchmarking of a Hybrid PIC Code For Dense Plasmas and Fast Ignition (Final Report)

Radiation processes play an important role in the study of both fast ignition and other inertial confinement schemes, such as plasma jet driven magneto-inertial fusion, both in their effect on energy balance, and in generating diagnostic signals. In the latter case, warm and hot dense matter may be produced by the convergence of a plasma shell formed by the merging of an assembly of high Mach number plasma jets. This innovative approach has the potential advantage of creating matter of high energy densities in voluminous amount compared with high power lasers or particle beams. An important application of this technology is as a plasma liner for the flux compression of magnetized plasma to create ultra-high magnetic fields and burning plasmas. HyperV Technologies Corp. has been developing plasma jet accelerator technology in both coaxial and linear railgun geometries to produce plasma jets of sufficient mass, density, and velocity to create such imploding plasma liners. An enabling tool for the development of this technology is the ability to model the plasma dynamics, not only in the accelerators themselves, but also in the resulting magnetized target plasma and within the merging/interacting plasma jets during transport to the target. Welch pioneered numerical modeling of such plasmas (including for fast ignition) using the LSP simulation code. Lsp is an electromagnetic, parallelized, plasma simulation code under development since 1995. It has a number of innovative features making it uniquely suitable for modeling high energy density plasmas including a hybrid fluid model for electrons that allows electrons in dense plasmas to be modeled with a kinetic or fluid treatment as appropriate. In addition to in-house use at Voss Scientific, several groups carrying out research in Fast Ignition (LLNL, SNL, UCSD, AWE (UK), and Imperial College (UK)) also use LSP. A collaborative team consisting of HyperV Technologies Corp., Voss Scientific LLC, FAR-TECH, Inc., Prism Computational Sciences, Inc. and Advanced Energy Systems Inc. joined efforts to develop new physics and numerical models for LSP in several key areas to enhance the ability of LSP to model high energy density plasmas (HEDP). This final report details those efforts. Areas addressed in this research effort include: adding radiation transport to LSP, first in 2D and then fully 3D, extending the EMHD model to 3D, implementing more advanced radiation and electrode plasma boundary conditions, and installing more efficient implicit numerical algorithms to speed complex 2-D and 3-D computations. The new capabilities allow modeling of the dominant processes in high energy density plasmas, and further assist the development and optimization of plasma jet accelerators, with particular attention to MHD instabilities and plasma/wall interaction (based on physical models for ion drag friction and ablation/erosion of the electrodes). In the first funding cycle we implemented a solver for the radiation diffusion equation. To solve this equation in 2-D, we used finite-differencing and applied the parallelized sparse-matrix solvers in the PETSc library (Argonne National Laboratory) to the resulting system of equations. A database of the necessary coefficients for materials of interest was assembled using the PROPACEOS and ATBASE codes from Prism. The model was benchmarked against Prism's 1-D radiation hydrodynamics code HELIOS, and against experimental data obtained from HyperV's separately funded plasma jet accelerator development program. Work in the second funding cycle focused on extending the radiation diffusion model to full 3-D, continued development of the EMHD model, optimizing the direct-implicit model to speed up calculations, add in multiply ionized atoms, and improved the way boundary conditions are handled in LSP. These new LSP capabilities were then used, along with analytic calculations and Mach2 runs, to investigate plasma jet merging, plasma detachment and transport, restrike and advanced jet accelerator design. In addition, a strong linkage to diagnostic measurements was made by modeling plasma jet experiments on PLX to support benchmarking of the code. A large number of upgrades and improvements advancing hybrid PIC algorithms were implemented in LSP during the second funding cycle. These include development of fully 3D radiation transport algorithms, new boundary conditions for plasma-electrode interactions, and a charge conserving equation of state that permits multiply ionized high-Z ions. The final funding cycle focused on 1) mitigating the effects of a slow-growing grid instability which is most pronounced in plasma jet frame expansion problems using the two-fluid Eulerian remap algorithm, 2) extension of the Eulerian Smoothing Algorithm to allow EOS/Radiation modeling, 3) simulations of collisionless shocks formed by jet merging, 4) simulations of merging jets using high-Z gases, 5) generation of PROPACEOS EOS/Opacity databases, 6) simulations of plasma jet transport experiments, 7) simulations of plasma jet penetration through transverse magnetic fields, and 8) GPU PIC code development. The tools developed during this project are applicable not only to the study of plasma jets, but also to a wide variety of HEDP plasmas of interest to DOE, including plasmas created in short-pulse laser experiments performed to study fast ignition concepts for inertial confinement fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

UCLA Final Technical Report for the "Community Petascale Project for Accelerator Science and Simulation”.

The UCLA Plasma Simulation Group is a major partner of the “Community Petascale Project for Accelerator Science and Simulation”. This is the final technical report. We include an overall summary, a list of publications, progress for the most recent year, and individual progress reports for each year. We have made tremendous progress during the three years. SciDAC funds have contributed to the development of a large number of skeleton codes that illustrate how to write PIC codes with a hierarchy of parallelism. These codes cover 2D and 3D as well as electrostatic solvers (which are used in beam dynamics codes and quasi-static codes) and electromagnetic solvers (which are used in plasma based accelerator codes). We also used these ideas to develop a GPU enabled version of OSIRIS. SciDAC funds were also contributed to the development of strategies to eliminate the Numerical Cerenkov Instability (NCI) which is an issue when carrying laser wakefield accelerator (LWFA) simulations in a boosted frame and when quantifying the emittance and energy spread of self-injected electron beams. This work included the development of a new code called UPIC-EMMA which is an FFT based electromagnetic PIC code and to new hybrid algorithms in OSIRIS. A new hybrid (PIC in r-z and gridless in φ) algorithm was implemented into OSIRIS. In this algorithm the fields and current are expanded into azimuthal harmonics and the complex amplitude for each harmonic is calculated separately. The contributions from each harmonic are summed and then used to push the particles. This algorithm permits modeling plasma based acceleration with some 3D effects but with the computational load of an 2D r-z PIC code. We developed a rigorously charge conserving current deposit for this algorithm. Very recently, we made progress in combining the speed up from the quasi-3D algorithm with that from the Lorentz boosted frame. SciDAC funds also contributed to the improvement and speed up of the quasi-static PIC code QuickPIC. We have also used our suite of PIC codes to make scientific discovery. Highlights include supporting FACET experiments which achieved the milestones of showing high beam loading and energy transfer efficiency from a drive electron beam to a witness electron beam and the discovery of a self-loading regime a for high gradient acceleration of a positron beam. Both of these experimental milestones were published in Nature together with supporting QuickPIC simulation results. Simulation results from QuickPIC were used on the cover of Nature in one case. We are also making progress on using highly resolved QuickPIC simulations to show that ion motion may not lead to catastrophic emittance growth for tightly focused electron bunches loaded into nonlinear wakefields. This could mean that fully self-consistent beam loading scenarios are possible. This work remains in progress. OSIRIS simulations were used to discover how 200 MeV electron rings are formed in LWFA experiments, on how to generate electrons that have a series of bunches on nanometer scale, and how to transport electron beams from (into) plasma sections into (from) conventional beam optic sections.

43 PARTICLE ACCELERATORS↗