Search NASA⌕ Search

SEARCH · Search NASA

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Cumulants of multiple conserved charges and global conservation laws

We analyze the behavior of cumulants of conserved charges in a subvolume of a thermal system with exact global conservation laws by extending a recently developed subensemble acceptance method (SAM) to multiple conserved charges. Explicit expressions for all diagonal and off-diagonal cumulants up to sixth order that relate them to the grand canonical susceptibilities are obtained. The derivation is presented for an arbitrary equation of state with an arbitrary number of different conserved charges. The global conservation effects cancel out in any ratio of two second order cumulants, in any ratio of two third order cumulants, as well as in a ratio of strongly intensive measures Σ and Δ involving any two conserved charges, making all these quantities particularly suitable for theory-to-experiment comparisons in heavy-ion collisions. We also show that the same cancellation occurs in correlators of a conserved charge, like the electric charge, with any non-conserved quantity such as net proton or net kaon number. The main results of the SAM are illustrated in the framework of the hadron resonance gas model. We also elucidate how net-proton and net-Λ fluctuations are affected by conservation of electric charge and strangeness in addition to baryon number.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Charge conservation and higher moments of charge fluctuations

Higher moments of distributions of net charge and baryon number in heavy-ion collisions have been proposed as signals of fundamental QCD phase transitions. In order to better understand background processes for these observables, models are presented which enable one to gauge the effects of local charge conservation, decays of resonances and clusters, Bose symmetrization, and volume fluctuations. Monte Carlo methods for generating samplings of particles consistent with local charge conservation are presented and are followed by a review of simple analytic models involving a single type of charge with a constant experimental efficiency. The main model consists of thermal emission superimposed onto a simple parametrization of collective flow, known as a blast wave, with emission being consistent with individual canonical ensembles. The spatial extent of local charge conservation is parameterized by the size and extent over which charge is conserved. Here, the sensitivity of third- and fourth-order moments, skewness and kurtosis, to these parameters, and to beam energy and baryon density is explored. Comparisons with STAR data show that a significant part of the observed non-Poissonian fluctuations in net-proton fluctuations are explained by charge and baryon-number conservation, but that measurements of the STAR collaboration for fluctuations of net electric charge significantly differ from expectations of the models presented here.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

On interpretation of fluctuations of conserved charges at high T

Abstract Fluctuations of conserved charges calculated on the lattice which can be measured experimentally, are well reproduced by a hadron resonanse gas model at temperatures below$$T_{ch} \sim 155$$ T ch ∼ 155 MeV and radically deviate from the hadron resonance gas predictions above the chiral restoration crossover. This behaviour is typically interpreted as an indication of deconfinement in the quark-gluon plasma regime. We present an argument that this interpretation may be too simple. The argument is based on the scaling of quantities with the number of colors: demonstration of deconfinement and QGP requires observable that is sensitive to$$\sim N_c^2$$ ∼ N c 2 gluons while the conserved charges are sensitive only to quarks and above$$T_{ch}$$ T ch scale as$$N_c^1$$ N c 1 . The latter scaling is consistent with the existence of an intermediate regime characterized by restored chiral symmetry and by approximate chiral spin symmetry which is a symmetry of confining interaction. In this regime the energy density, pressure and entropy density scale as$$N_c^1$$ N c 1 . In the large$$N_c$$ N c limit this regime might become a distinct phase separated from the hadron gas and from QGP by phase transitions. A natural observable that associates with deconfinement and is directly sensitive to deconfined$$N_c^2-1$$ N c 2 - 1 gluons is the Polyakov loop; in the$$N_c=3$$ N c = 3 world it remains very close to 0 at temperatures well above chiral crossover, reaches the value$$\sim 0.5$$ ∼ 0.5 around$$\sim 3T_{ch}$$ ∼ 3 T ch and the value close to 1 at temperatures$$\sim 1$$ ∼ 1 GeV.

Physics↗
DOE OSTI2022

Monte Carlo event generator for initial conditions of conserved charges in nuclear geometry

At top collider energies where baryon stopping is negligible, the initial state of heavy-ion collisions is overall charge neutral and predominantly composed of gluons. Nevertheless, there can also be significant local fluctuations of the baryon number, strangeness, and electric charge densities about zero, perturbatively corresponding to the production of quark/antiquark pairs. These previously ignored local charge fluctuations can permit the study of charge diffusion in the quark-gluon plasma (QGP), even at top collider energies. In this paper we present a new model denoted ICCING (initial conserved charges in nuclear geometry) which can reconstruct the initial conditions of conserved charges in the QGP by sampling a (g→q¯q) splitting probability over the initial energy density. We find that the new charge distributions generally differ from the bulk energy density; in particular, the strangeness distribution is significantly more eccentric than standard bulk observables and appears to be associated with the geometry of hot spots in the initial state. The new information provided by these conserved charges opens the door to studying a wealth of new charge- and flavor-dependent correlations in the initial state and ultimately the charge transport parameters of the QGP.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

New proxies for second-order cumulants of conserved charges in heavy-ion collisions within the EPOS4 framework

Proxies for cumulants of baryon number 𝐵, electric charge 𝑄, and strangeness 𝑆 are usually measured in heavy-ion collisions via moments of net-number distribution of given hadronic species. Since these cumulants of conserved charges are expected to be sensitive to the existence of a critical point in the phase diagram of nuclear matter, it is crucial to ensure that the proxies used as substitutes are as close to them as possible. Hence, we use the EPOS 4 framework to generate Au + Au collisions at several collision energies of the BNL Relativistic Heavy Ion Collider beam energy scan. We compute second-order net cumulants of 𝜋, 𝐾, and 𝑝, for which experimental data have been published as well as the corresponding conserved charge cumulants. We then compare them with proxies, defined in previous lattice QCD and hadron resonance gas model studies, which are shown to reproduce more accurately their associated conserved charge cumulants. We investigate the impact of hadronic rescatterings occurring in the late evolution of the system on these quantities, as well as the amount of signal actually originating from the bulk medium which endures a phase transition.

Physics↗
DOE OSTI2024

An implicit particle code with exact energy and charge conservation for studies of dense plasmas in axisymmetric geometries

A collisional particle code based on implicit energy- and charge-conserving methods in axisymmetric geometries is presented. A new particle pusher for axisymmetric systems is introduced that is compatible with exact energy and charge conservation and yields improved accuracy compared to other methods. How to appropriately treat all aspects of the algorithm near the r = 0 axis of symmetry is described in detail. Here, the axisymmetric model is verified by simulating the free expansion of a plasma sphere in 2D cylindrical and 1D spherical geometries. The algorithm's ability to study the dynamic compression of a dense plasma is illustrated by simulating the dynamic Z-pinch in 1D cylindrical geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

EZ: An efficient, charge conserving current deposition algorithm for electromagnetic particle-in-cell simulations

Here, we present EZ, a novel current deposition algorithm for particle-in-cell (PIC) simulations. EZ calculates the current density on the electromagnetic grid due to macro-particle motion within a time step by solving the continuity equation of electrodynamics. Being a charge conserving hybridization of Esirkepov's method and ZigZag, we refer to it as “EZ” as shorthand for “Esirkepov meets ZigZag”. Simulations of a warm, relativistic plasma with PIConGPU show that EZ achieves the same level of charge conservation as the commonly used method by Esirkepov, yet reaches higher performance for macro-particle assignment-functions up to third-order. In addition to a detailed description of the functioning of EZ, reasons for the expected and observed performance increase are given, and guidelines for its implementation aiming at highest performance on GPUs are provided.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Viscous Gubser flow with conserved charges to benchmark fluid simulations

We present semi-analytical solutions for the evolution of both the temperature and chemical potentials for viscous Gubser flow with conserved charges. Such a solution can be especially useful in testing numerical codes intended to simulate relativistic fluids with large chemical potentials. The freeze-out hypersurface profiles for constant energy density are calculated, along with the corresponding normal vectors, and presented as a new unit test for numerical codes. We also compare the influence of the equation of state on the semi-analytical solutions. We benchmark the newly developed smoothed particle hydrodynamics code ccake that includes both shear viscosity and three conserved charges. Here, the numerical solutions are in excellent agreement with the semi-analytical solution and are also able to accurately reproduce the hypersurface at freeze-out.

Hydrodynamic models↗
DOE OSTI2021

Time integrator agnostic charge conserving finite element PIC

Developing particle-in-cell (PIC) methods using finite element basis sets, and without auxiliary divergence cleaning methods, was a longstanding problem until recently. It was shown that if consistent spatial basis functions are used, one can indeed create a methodology that was charge conserving, albeit using a leapfrog time stepping method. While this is a significant advance, leapfrog schemes are only conditionally stable and time step sizes are closely tied to the underlying mesh. Ideally, to take full advantage of advances in finite element methods (FEMs), one needs a charge conserving PIC methodology that is agnostic to the time stepping method. This is the principal contribution of this paper. In what follows, we shall develop this methodology, prove that both charge and Gauss’ laws are discretely satisfied at every time step, provide the necessary details to implement this methodology for both the wave equation FEM and Maxwell solver FEM, and finally demonstrate its efficacy on a suite of test problems. The method will be demonstrated by single particle evolution, non-neutral beams with space-charge, and adiabatic expansion of a neutral plasma, where the Debye length has been resolved, and real mass ratios are used.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Cross-Correlators of Conserved Charges in QCD

We present cross-correlators of QCD conserved charges at µ B = 0 from lattice simulations and perform a Hadron Resonance Gas (HRG) model analysis to break down the hadronic contributions to these correlators. We construct a suitable hadronic proxy for the ratio –χ$^{BS}_{11}$ /χ$^{S}_{2}$ and discuss the dependence on the chemical potential and experimental cuts. Furthermore, we then perform a comparison to preliminary STAR results and comment on a possible direct comparison of lattice and experiment

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Pre-equilibrium evolution of conserved charges with initial conditions in the ICCING Monte Carlo event generator

Heavy-ion collisions can be well described through relativistic viscous hydrodynamics, but questions still remain when hydrodynamics is applicable because the initial state may begin very far from equilibrium. Thus, a pre-equilibrium evolution phase is used to bridge the gap between the initial state and hydrodynamics. KøMPøST is one such pre-equilibrium model that propagates the energy-momentum tensor by decomposing it into the background and fluctuations around that background, whose evolution is captured by Green's functions. We extend this formalism to include conserved charges and calculate the corresponding nonequilibrium Green's functions in the relaxation-time approximation. The ICCING algorithm initializes conserved charges in the initial state by sampling $g$ → $q$$\overline{q}$ splitting probabilities and is, thus, perfectly positioned to implement Green's functions for charge propagation. As a result, we show that this method alters the initial-state charge geometries and is applicable in central to mid-central collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

An implicit particle code with exact energy and charge conservation for electromagnetic studies of dense plasmas

A collisional particle code based on implicit energy- and charge-conserving methods is presented. A modified version of the particle-suppressed Jacobian-Free Newton-Krylov method that can enhance the solver efficiency is introduced. Mathematically, it is shown that this new approach can be viewed as a fixed-point iteration method for the particle positions. The model can exactly conserve global energy and local charge and can efficiently use time steps larger than the plasma period. In conclusion, the algorithm's ability to simulate dense plasmas accurately and efficiently is quantified by simulating the dynamic compression of a plasma slab via a magnetic piston in 1D planar geometry.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Constraining the hadronic spectrum and repulsive interactions in a hadron resonance gas via fluctuations of conserved charges

We simultaneously incorporate two common extensions of the hadron resonance gas model, namely the addition of extra, unconfirmed resonances to the particle list and the excluded volume repulsive interactions. We emphasize the complementary nature of these two extensions and identify combinations of conserved charge susceptibilities that allow to constrain them separately. In particular, ratios of second-order susceptibilities like $\chi_{11}^{BQ}/\chi_2^B$ and $\chi_{11}^{BS}/\chi_2^B$ are sensitive only to the baryon spectrum, while fourth-to-second order ratios like $\chi_4^B/\chi_2^B$, $\chi_{31}^{BS}/\chi_{11}^{BS}$, or $\chi_{31}^{BQ}/\chi_{11}^{BQ}$ are mainly determined by repulsive interactions. Analysis of the available lattice results suggests the presence of both the extra states in the baryon-strangeness sector and the repulsive baryonic interaction, with indications that hyperons have a smaller repulsive core than non-strange baryons. Furthermore, the modified hadron resonance gas model presented here significantly improves the description of lattice QCD susceptibilities at chemical freeze-out and can be used for the analysis of event-by-event fluctuations in heavy-ion collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

A pseudospectral implicit particle-in-cell method with exact energy and charge conservation

The standard particle-in-cell (PIC) method employs explicit finite-difference (FD) methods (e.g. the leap-frog scheme) for both spatial and temporal integrations. Here, we employ a pseudospectral method for solving the Poisson equation and a fully implicit time integration to achieve exact energy conservation. The advantage of a pseudospectral field solver is its spectral accuracy in solving field solutions. Earlier studies of implicit time integration of PIC FD equations can enforce exact energy exchange between field and particles, resulting in exact energy-conserving schemes. Here, we prove that the exact energy conservation property can be carried over to the pseudospectral scheme. Simultaneously, we provide a solution to ensure a pseudospectral charge continuity equation. We demonstrate the new scheme in a 2D electrostatic PIC code. In conclusion, theoretical results are confirmed via numerical examples.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Quasi-Helmholtz decomposition, Gauss' laws and charge conservation for finite element particle-in-cell

Development of particle-in-cell (PIC) methods using finite element based methods (FEMs) have been a topic of renewed interest; this has largely been driven by (a) the ability of finite element methods to better model geometry, (b) better understanding of function spaces that are necessary to represent all Maxwell quantities, and (c) more recently, the fundamental rubrics that should be obeyed in space and time so as to satisfy Gauss' laws and the equation of continuity. In that vein, methods have been developed recently that satisfy these equations and are agnostic to time stepping methods. While this development is indeed a significant advance, it should be noted that implicit FEM transient solvers support an underlying null space that corresponds to a gradient of a scalar potential ∇Φ(r) (or t∇Φ(r) in the case of wave equation solvers). While explicit schemes do not suffer from this drawback, they are only conditionally stable with time step sizes that are mesh dependent and very small. Furthermore, the null space produces spurious charge that can corrupt the desired physics of a PIC simulation. The way to overcome this bottleneck, and indeed, satisfy all four Maxwell's equation is to use a quasi-Helmholtz formulation on a tessellation. In the re-formulation presented, we strictly satisfy the equation of continuity and Gauss' laws for both the electric and magnetic flux densities. Results illustrating the efficacy of this scheme will be demonstrated by analyzing non-neutral beams with space-charge and the adiabatic expansion of a neutral plasma with realistic parameters (Debye length and real mass ratios).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗