Search NASA⌕ Search

SEARCH · Search NASA

Results for “fixed points”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Godiva IV Simulated Radiation Field Characterization and Variance Reduction

Godiva IV is a system comprised of highly enriched uranium alloyed with molybdenum in the form of fuel plate rings. The reactor, along with its predecessors, was designed with the unique ability to satisfy interests in the super-prompt-critical reactor operation space. Originally, the reactor was part of the Los Alamos Critical Experiments Facility (LACEF) at Technical Area-18 (TA-18). The radiation field around Godiva at this facility was well characterized and understood. As a fast neutron system, the neutron spectrum in and around Godiva was close to a Watt Fission spectrum. The Kiva where Godiva IV was located at LACEF was made of thin, sheet metal walls which did not contribute significantly to the neutron spectrum. Following the transition of LACEF to the National Critical Experiments and Research Center (NCERC) in Nevada, Godiva-IV was moved from TA-18 to the Device Assembly Facility (DAF) at the Nevada National Security Site (NNSS). Part of this move brought renewed interest in radiation field characterization. The new facility introduced significant changes to the environment surrounding Godiva, and preliminary foil irradiation results suggested that the room contribution to the neutron spectrum was significant. Unlike at TA-18, a large thermal neutron signature was added to the fast spectrum from Godiva due to significant room return. A primary goal due to the additional complexity that the room return adds to the Godiva IV radiation emission spectrum was the development of an efficient Monte Carlo N-Particle (MCNP) calculation capable of characterizing the neutron spectrum anywhere in the room around Godiva. A campaign of activation foil irradiations and analysis were completed to support the validation of the MCNP model. The modeling of these foils in MCNP can be easily done with a standard volumetric neutron flux tally. However, given the multitude of locations and reaction rates to be modeled, further steps must be taken to increase the efficiency of these calculations in MCNP. During this study, a benchmark model currently under development for Godiva IV was used. A qualitative assessment of the thermal neutron contributors was performed using spatial neutron distribution plots. Additional detail was added to the model based on the qualitative results showing the thermal spectrum’s large sensitivity to hydrogenous material. Neutron energy spectra was evaluated at discrete locations in the room around Godiva to quantify the relative contribution of various components. It was discovered that the concrete walls are the largest contributor to the thermal signature, with minor contributions from plastic components surrounding Godiva. Following these results, two different variance reduction techniques were implemented to improve the problem efficiency in these calculations. In the first approach, an F5 point detector tally was implemented in the standard Godiva IV criticality problem. The second approach involved a weight-window generator implementation with an F5 point detector tally in a fixed source problem. The weight window implementation reduced the runtime from 42739.55 minutes to 1803.34 minutes (computer time), compared to the F5 KCODE implementation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Correlations between the Neutron Star Mass–Radius Relation and the Equation of State of Dense Matter

We develop an analytic method of inverting the Tolman–Oppenheimer–Volkoff relations to high accuracy. In principle, a specified energy density–pressure relation gives a unique mass–radius (M–R) relation and vice versa. Our method is developed from the strong correlations that are shown to exist between the neutron star mass–radius curve and the equation of state (EOS) or pressure–energy density relation. Selecting points that have masses equal to fixed fractions of the maximum mass, we find a semi-universal power-law relation between the central energy densities, pressures, sound speeds, chemical potentials, and number densities of those stars, with the maximum mass and the radii of one or more fractional maximum mass points. Rms fitting accuracies, for EOSs without large first-order phase transitions, are typically 0.5% for all quantities at all mass points. The method also works well, although less accurately, in reconstructing the EOS of hybrid stars with first-order phase transitions. These results permit, in effect, an analytic method of inverting an arbitrary M–R curve to yield its underlying EOS. We discuss applications of this inversion technique to the inference of the dense matter EOS from measurements of neutron star masses and radii as a possible alternative to traditional Bayesian approaches.

Bayesian statistics↗
DOE OSTI2022

Detailed analysis of excited-state systematics in a lattice QCD calculation of 𝑔 𝐴

Excited state contamination remains one of the most challenging sources of systematic uncertainty to control in lattice QCD calculations of nucleon matrix elements and form factors: early time separations are contaminated by excited states and late times suffer from an exponentially bad signal-to-noise problem. High-statistics calculations at large time separations ≳ 1 fm are commonly used to combat these issues. In this work, focusing on g A , we explore the alternative strategy of utilizing a large number of relatively low-statistics calculations at short to medium time separations (0.2–1 fm), combined with a multistate analysis. On an ensemble with a pion mass of approximately 310 MeV and a lattice spacing of approximately 0.09 fm, we find this provides a more robust and economical method of quantifying and controlling the excited state systematic uncertainty. A quantitative separation of various types of excited states enables the identification of the transition matrix elements as the dominant contamination. The excited state contamination of the Feynman-Hellmann correlation function is found to reduce to the 1% level at approximately 1 fm while, for the more standard three-point functions, this does not occur until after 2 fm. Critical to our findings is the use of a global minimization, rather than fixing the spectrum from the two-point functions and using them as input to the three-point analysis. We find that the ground state parameters determined in such a global analysis are stable against variations in the excited state model, the number of excited states, and the truncation of early-time or late-time numerical data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Hydrogen Detection Strategies to Support H2@SCALE - The NREL Sensor Laboratory

Hydrogen represents a major pathway to decarbonize and stabilize the national and international energy industry and select manufacturing markets. To facilitate the development of hydrogen markets, the US Department of Energy initiated H2@Scale to bring together stakeholders to advance affordable hydrogen production, transport, storage, and utilization to increase revenue opportunities in multiple energy sectors. One major impediment to hydrogen implementation is cost. To expedite the use of hydrogen in energy and other markets, the United States announced in 2021 the Hydrogen Shot, which seeks to reduce the cost of clean hydrogen by 80% to $1 per 1 kilogram in 1 decade ("1 1 1"). As the cost of hydrogen drops, new applications will emerge that will require unique configurations of existing equipment and infrastructure, and eventually lead to advances in the generation and utilization of hydrogen. As the hydrogen economy expands, sensors and detection methods will need to adapt to changing infrastructure demands to address the primary targets of health & safety, emissions monitoring, and process control. The NREL Sensor Laboratory is playing a pivotal role in advancing the use of hydrogen sensors and detection methodologies in each of these categories to support DOE's mission for safe and efficient utilization in emerging markets. Health & safety monitors are required to ensure that operators and facilities can react to unintended hydrogen releases, either as GH2, LH2, or as a constituent of blends (e.g., natural gas or ammonia). Current detection methodologies focus on safety applications to detect near its lower flammable limit (4 vol %), and typically include point sensors in applications such as fixed or mobile detectors (e.g., personal gas monitors). Methodologies amenable for area detection include acoustic, emerging optical imaging methods, and flame detectors. Comparable detection strategies can be utilized for emissions monitoring and quantization, however few methods can simultaneously cover both low (emissions) and high (health & safety) levels. Deployment of emission level detectors will be required to 1) reduce product loss through small but potentially significant leaks from an environmental or cost perspective, 2) reduce downtime of high demand systems by early identification of eminent system failures (leaks through pump or compressor seals indicative of impending failure), and 3) address potential emission monitoring requirements that may be set by regulating bodies. The first two points should be adopted by industry to reduce the cost-of-goods-sold. The third main category for hydrogen detection relates to process control and may be advantageous for many existing applications. Two main applications are emerging. For example, the purity requirements for hydrogen that is dispensed from refueling systems for hydrogen fuel cell electric vehicles (FCEV) is rigorously regulated by the Standard SAE J2719, which prescribes maximum allowable levels of multiple impurities in the hydrogen fuel and must be verified by a regulatory body. Hydrogen contaminant detectors (HCD) integrated to the fueling station can assure this compliance. HCDs must be able operate in 100% H2 backgrounds and be able to distinguish between multiple contaminants at low ppm to low ppb levels. Secondly, as a strategy to decarbonize the natural gas grid, there are proposals to blend hydrogen with natural gas. This blending will affect transport applications (pipeline infrastructure), stationary combustion systems (turbines), and consumer and commercial appliances. In the short-term, hydrogen levels up to 20% are proposed. Variations in the hydrogen level can have dramatic impact on the combustion process and on the potential response of safety sensors. These mixtures may be regulated so that the concentration at a delivery point must be monitored with high precision. However, routine maintenance may introduce background gases such as ambient air (with water) or maintenance gases (introduced with welding processes or adhesive outgassing.) Therefore, the detection methodology must be robust enough to recover or respond to various contaminants. Several reviews can be found in literature addressing sensing and detection technologies, including their limitations and applications. However, for most applications, limitations can be alleviated by combining various detection techniques either through system integration or implementation of machine learning methods (artificial intelligence). In this presentation, we will discuss several applications, highlight their current approach for hydrogen detection, and suggest detection strategies to supplement their limitations.

ENERGY STORAGE,HYDROGEN↗
DOE OSTI2024

Pressure-stabilized fixed-stress iterative solutions of compositional poromechanics

We consider the numerical behavior of the fixed-stress splitting method for coupled poromechanics as undrained regimes are approached. We explain that pressure stability is related to the splitting error of the scheme, not the fact that the discrete saddle point matrix never appears in the fixed-stress approach. This observation reconciles previous results regarding the pressure stability of the splitting method. Using examples of compositional poromechanics with application to geological CO sequestration, we see that solutions obtained using the fixed-stress scheme with a low order finite element-finite volume discretization which is not inherently inf-sup stable can exhibit the same pressure oscillations obtained with the corresponding fully implicit scheme. Moreover, pressure jump stabilization can effectively remove these spurious oscillations in the fixed-stress setting, while also improving the efficiency of the scheme in terms of the number of iterations required at every time step to reach convergence.

42 ENGINEERING↗
DOE OSTI2021

Reconfigurable Two‐Dimensional DNA Lattices: Static and Dynamic Angle Control

Abstract Branched DNA motifs serve as the basic construction elements for all synthetic DNA nanostructures. However, precise control of branching orientation remains a key challenge to further heighten the overall structural order. In this study, we use two strategies to control the branching orientation. The first one is based on immobile Holliday junctions which employ specific nucleotide sequences at the branch points which dictate their orientation. The second strategy is to use angle‐enforcing struts to fix the branching orientation with flexible spacers at the branch points. We have also demonstrated that the branching orientation control can be achieved dynamically, either by canonical Watson–Crick base pairing or non‐canonical nucleobase interactions (e.g., i‐motif and G‐quadruplex). With precise angle control and feedback from the chemical environment, these results will enable novel DNA nanomechanical sensing devices, and precisely‐ordered three‐dimensional architectures.

Wang, Wen↗
DOE OSTI2023

Arrhenian to non-Arrhenian crossover in glass melt viscosity

The activation energy for glass melt viscosity is a function of temperature, i.e., non-Arrhenian, within a middle range of viscosity values (typically, between 10 2 and 10 10 Pa s) and constant, i.e., Arrhenian, outside this range. At the low-viscosity end, the Arrhenian to non-Arrhenian transition (termed the crossover) falls within the glass processing range: glass melting and fining occur within the Arrhenian range while glass forming occurs within the non-Arrhenian range. By the Adam-Gibbs equation, the configuration entropy is nearly constant at high temperatures (low viscosities), where the glass melt turns into a simple liquid, and decreases with decreasing temperature (increasing viscosity) as the glass is increasingly polymerized. In the present study, the power-law function and hyperbolic tangent function are used to express the configuration entropy as a function of temperature. The viscosity crossover is defined as the point at which onset of structural grouping is causing the configuration entropy to deviate from its high-temperature value. In conclusion, it appears that the limited range of the configuration entropy of the fully depolymerized melt imposes a restriction on the crossover viscosity; for well-defined glass families, such as float glasses or nuclear waste glasses, the crossover point can be defined using a fixed, i.e., composition-independent, viscosity value.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Reconfigurable Two-Dimensional DNA Lattices: Static and Dynamic Angle Control

Branched DNA motifs serve as the basic construction elements for all synthetic DNA nanostructures. However, precise control of branching orientation remains a key challenge to further heighten the overall structural order. Here in this study, we use two strategies to control the branching orientation. The first one is based on immobile Holliday junctions which employ specific nucleotide sequences at the branch points which dictate their orientation. The second strategy is to use angle-enforcing struts to fix the branching orientation with flexible spacers at the branch points. We have also demonstrated that the branching orientation control can be achieved dynamically, either by canonical Watson–Crick base pairing or non-canonical nucleobase interactions (e.g., i-motif and G-quadruplex). With precise angle control and feedback from the chemical environment, these results will enable novel DNA nanomechanical sensing devices, and precisely-ordered three-dimensional architectures.

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

Linear point and sound horizon as purely geometric standard rulers

The baryon acoustic oscillations feature (BAO) imprinted in the clustering correlation function is known to furnish us cosmic distance determinations that are independent of the cosmological-background model and the primordial perturbation parameters. These measurements can be accomplished rigorously by means of the purely geometric BAO methods. To date two different purely geometric BAO approaches have been proposed. The first exploits the linear-point standard ruler. The second, called correlation-function model-fitting, exploits the sound-horizon standard ruler. A key difference between them is that, when estimated from clustering data, the linear point makes use of a cosmological-model-independent procedure to extract the ratio of the ruler to the cosmic distance, while the correlation-function model-fitting relies on a phenomenological cosmological model for the correlation function. Nevertheless the two rulers need to be precisely defined independently of any specific observable (e.g., the BAO). We define the linear point and sound horizon and we fully characterize and compare the two rulers’ cosmological-parameter dependence. We find that they are both geometrical (i.e., independent of the primordial cosmological parameters) within the required accuracy, and that they have the same parameter dependence for a wide range of parameter values. We estimate the rulers’ best-fit values and errors, given the cosmological constraints obtained by the Planck Satellite team from their measurements of the cosmic microwave background temperature and polarization anisotropies. We do this for three different cosmological models encompassed by the purely geometric BAO methods. In each case we find that the relative errors of the two rulers coincide and they are insensitive to the assumed cosmological model. Interestingly both the linear point and the sound horizon shift by 0.5σ when we do not fix the spatial geometry to be flat in Λ CDM. Finally, this points toward a sensitivity of the rulers to different cosmological models when they are estimated from the cosmic microwave background.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A precision test of averaging in AdS/CFT

We reconsider the role of wormholes in the AdS/CFT correspondence. We focus on Euclidean wormholes that connect two asymptotically AdS or hyperbolic regions with $\mathbb{S}$ 1 × $\mathbb{S}$ d-1 boundary. There is no solution to Einstein’s equations of this sort, as the wormholes possess a modulus that runs to infinity. To find on-shell wormholes we must stabilize this modulus, which we can do by fixing the total energy on the two boundaries. Such a wormhole gives the saddle point approximation to a non-standard problem in quantum gravity, where we fix two asymptotic boundaries and constrain the common energy. Crucially the dual quantity does not factorize even when the bulk is dual to a single CFT, on account of the fixed energy constraint. From this quantity we extract a smeared version of the microcanonical spectral form factor. For a chaotic theory this quantity is self-averaging, i.e. well-approximated by averaging over energy windows, or over coupling constants. We go on to give a precision test involving the microcanonical spectral form factor where the two replicas have slightly different coupling constants. In chaotic theories this form factor is known to smoothly decay at a rate universally predicted in terms of one replica physics, provided that there is an average either over a window or over couplings. We compute the expected decay rate for holographic theories, and the form factor from a wormhole, and the two exactly agree for a wide range of two-derivative effective field theories in AdS. This gives a precision test of averaging in AdS/CFT. Our results interpret a number of confusing facts about wormholes and factorization in AdS and suggest that we should regard gravitational effective field theory as a mesoscopic description, analogous to semiclassical mesoscopic descriptions of quantum chaotic systems.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Counteracting Systemic Bias in the Lab, Field, and Classroom

Grappling with systemic discrimination and bias in geosciences can be overwhelming to the point that one may feel powerless to fix them. Despite the sweeping nature of this challenge, faculty, principal investigators, and other scientists with leadership roles have unparalleled power to mitigate harm in environments they oversee. Here, we identify ways that scientists in these roles can immediately address bias in three common spaces—the lab, field, and classroom. We highlight key actions that can be taken to improve the quality of life of marginalized students and other trainees quickly, while important but comparatively slow institutional changes proceed.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2026

Rapidity-dependent spin decomposition of the nucleon

We revisit the two-dimensional Fourier transform of generalized parton distributions (GPDs) at nonzero skewness. At 𝜂 = 0 it reduces to the standard impact-parameter density, while at 𝜂 ≠ 0 it is an off-forward amplitude that we interpret as a genuine parton–nucleon correlation. Its overall strength (the transverse-plane integral of the density) is fixed by the GPD at the kinematic point 𝑡 =−𝑐 𝜂 =−4⁢𝜂 2 ⁢𝑚$^{2}_{𝑁}$/(1 − 𝜂 2 ) and decreases monotonically with the rapidity gap Δ⁢𝑦 = ln⁡[(1+𝜂)/(1−𝜂)] = 2 artanh⁡(𝜂). This rapidity dependence implies rapidity-modified Ji identities that connect helicity, orbital, and total angular momenta of the correlation in closed form. To quantify these effects, we construct leading-twist quark and gluon GPDs in a string-based conformal framework: conformal moments are parametrized by linear open- and closed-string Regge trajectories with slopes constrained by parton distribution functions (PDFs), hadron/glueball spectroscopy, and form-factor data, and GPDs are reconstructed over the full (𝑥,𝜂,𝑡) domain by Mellin-Barnes inversion with next-to-leading order evolution. We find qualitative agreement (and fair quantitative agreement within quoted uncertainties) for several moments and selected nonsinglet 𝑥-space channels at 𝜇 = 2 GeV when compared with lattice QCD, while we also identify channels with visible tension and discuss likely sources (PDF priors and 𝑡-slope systematics).

Gauge-gravity dualities↗
DOE OSTI2022

Quantifying Errors in Effective Cluster Interactions of Lattice Gas Cluster Expansions

The promise of lattice gas (LG) cluster expansions (CEs) is that they can describe a given system property to any level of accuracy since the orthogonal “cluster basis functions” span the complete space of available configurations. Such an approach can be constructed to an arbitrarily large surface of a finite number of distinct adsorption sites. Unfortunately, this is only true for the case of an ideal, fixed lattice decorated with components at precise lattice points (the lattice “sites”) with no distortions or relaxations subsequently allowed. Since most systems, and surfaces specifically, do not conform to such an ideal set of constraints, errors in LG CEs must be expected or CE convergence severely hampered. Beyond this, numerical errors in the provided data can complicate the proper construction of a truly predictive and/or physically significant CE. In this work, we show here how reliance on typical statistical tools like confidence intervals cannot be expected to provide an accurate representation of the uncertainty of effective cluster interactions (ECIs) in the CE due to the nature of the target ab initio data and the nature of CEs themselves. We develop a method for estimating these errors that does not rely on statistical assumptions about the model or data. We then use these ECI errors to quantify fundamental consequences on the uncertainty of ECIs in CEs built from O/Fe(100) data whose surface and adsorbates have been allowed to relax in a typical manner and from O/Fe(100) data whose surface and adsorbates are fixed in ideal lattice positions. We also quantify the effect of using a different density functional theory exchange–correlation functional, using these ECI errors to assess the significance in any deviations. In both cases, our method is shown to have remarkable utility in the quantification of errors in the ECIs of CEs. While we stick to the lattice gas convention in this work, the method is equally applicable to the Ising convention or, in principle, any linear model of sufficient complexity.

08 HYDROGEN↗
DOE OSTI2025

A Particle-in-Cell Method for Plasmas with a Generalized Momentum Formulation, Part I: Model Formulation

Here, this paper formulates a new particle-in-cell method for the Vlasov–Maxwell system. Under the Lorenz gauge condition, Maxwell’s equations for the electromagnetic fields can be written as a collection of scalar and vector wave equations. The use of potentials for the fields motivates the adoption of a Hamiltonian formulation for particles that employs the generalized (conjugate) momentum. A notable advantage offered by the Hamiltonian formulation is the elimination of time derivatives in the Lorenz gauge formulation that are required by the standard Newton–Lorentz treatment of the particles. This allows the fields to retain the full time-accuracy guaranteed by the field solver. The resulting updates for particles require only knowledge of the fields and their spatial derivatives. An analytical method for constructing these spatial derivatives is presented that exploits the underlying integral solution used in the field solver for the wave equations. Moreover, these derivatives are demonstrated to converge at the same rate as the fields in both time and space. The Method of Lines Transpose field solver we consider in this work is globally first-order accurate in time and high-order accurate in space (e.g., fourth- and fifth-order) and belongs to a larger class of methods which are unconditionally stable, can address geometry, and leverage $\mathcal {O}(N)$ fast summation methods for efficiency. We demonstrate the method on several well-established benchmark problems on bounded domains, including a plasma sheath as well as a relativistic particle beam. The efficacy of the proposed formulation is established by comparing with a second-order accurate finite-difference time-domain method that employs a leapfrog time advance for particles and a charge conserving map suitable for bounded domains. The new method shows mesh-independent numerical heating properties even in cases where the plasma Debye length is smaller than the grid spacing. This is an important feature of the new method for problems defined on bounded domains, because it permits the use of coarser grids in space in the representation of the fields. Such a capability has significant implications for the simulation of plasmas in bounded domains with complex geometry, where the ratio between the largest and smallest cells can vary significantly. The use of high-order spatial approximations in the new method also means that fewer grid points are required in order to achieve a fixed accuracy. Our results also suggest that the new method can be used with fewer simulation particles per cell compared to the benchmark explicit method, which permits further computational savings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Wormholes and black hole microstates in AdS/CFT

It has long been known that the coarse-grained approximation to the black hole density of states can be computed using classical Euclidean gravity. In this work we argue for another entry in the dictionary between Euclidean gravity and black hole physics, namely that Euclidean wormholes describe a coarse-grained approximation to the energy level statistics of black hole microstates. To do so we use the method of constrained instantons to obtain an integral representation of wormhole amplitudes in Einstein gravity and in full-fledged AdS/CFT. These amplitudes are non-perturbative corrections to the two-boundary problem in AdS quantum gravity. The full amplitude is likely UV sensitive, dominated by small wormholes, but we show it admits an integral transformation with a macroscopic, weakly curved saddle-point approximation. The saddle is the “double cone” geometry of Saad, Shenker, and Stanford, with fixed moduli. In the boundary description this saddle appears to dominate a smeared version of the connected two-point function of the black hole density of states, and suggests level repulsion in the microstate spectrum. Using these methods we further study Euclidean wormholes in pure Einstein gravity and in IIB supergravity on Euclidean AdS 5 × S 5 . We address the perturbative stability of these backgrounds and study brane nucleation instabilities in 10d supergravity. In particular, brane nucleation instabilities of the Euclidean wormholes are lifted by the analytic continuation required to obtain the Lorentzian spectral form factor from gravity. Our results indicate a factorization paradox in AdS/CFT.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Measurements of Proton High-Order Cumulants in s N N = 3 GeV Au + Au Collisions and Implications for the QCD Critical Point

We report cumulants of the proton multiplicity distribution from dedicated fixed-target Au + Au collisions at $\sqrt{s_{NN}}$ = 3.0 GeV , measured by the STAR experiment in the kinematic acceptance of rapidity ( y ) and transverse momentum ( p T ) within - 0.5 < y < 0 and 0.4 < pT < 2.0 GeV / c . In the most central 0%–5% collisions, a proton cumulant ratio is measured to be C 4 / C 2 = - 0.85 ± 0.09 ( stat ) ± 0.82 ( syst ) , which is 2 σ below the Poisson baseline with respect to both the statistical and systematic uncertainties. The hadronic transport UrQMD model reproduces our C 4 / C 2 in the measured acceptance. Compared to higher energy results and the transport model calculations, the suppression in C 4 / C 2 is consistent with fluctuations driven by baryon number conservation and indicates an energy regime dominated by hadronic interactions. These data imply that the QCD critical region, if created in heavy-ion collisions, could only exist at energies higher than 3 GeV.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Covariance matrices for variance-suppressed simulations

ABSTRACT Cosmological N-body simulations provide numerical predictions of the structure of the Universe against which to compare data from ongoing and future surveys, but the growing volume of the Universe mapped by surveys requires correspondingly lower statistical uncertainties in simulations, usually achieved by increasing simulation sizes at the expense of computational power. It was recently proposed to reduce simulation variance without incurring additional computational costs by adopting fixed-amplitude initial conditions. This method has been demonstrated not to introduce bias in various statistics, including the two-point statistics of galaxy samples typically used for extracting cosmological parameters from galaxy redshift survey data, but requires us to revisit current methods for estimating covariance matrices of clustering statistics for simulations. In this work, we find that it is not trivial to construct covariance matrices analytically for fixed-amplitude simulations, but we demonstrate that ezmock (Effective Zel’dovich approximation mock catalogue), the most efficient method for constructing mock catalogues with accurate two- and three-point statistics, provides reasonable covariance matrix estimates for such simulations. We further examine how the variance suppression obtained by amplitude-fixing depends on three-point clustering, small-scale clustering, and galaxy bias, and propose intuitive explanations for the effects we observe based on the ezmock bias model.

79 ASTRONOMY AND ASTROPHYSICS↗