SEARCH · Search NASA

Results for “statistical models”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,207 records · Page 67

All sources NASA NTRS DOE OSTI

NASA NTRS1992

A stochastic damage model for the rupture prediction of a multi-phase solid. I - Parametric studies. II - Statistical approach

A stochastic damage model for predicting the rupture of a brittle multiphase material is developed, based on the microcrack-macrocrack interaction. The model, which incorporates uncertainties in locations, orientations, and numbers of microcracks, characterizes damage by microcracking and fracture by macrocracking. A parametric study is carried out to investigate the change of the stress intensity at the macrocrack tip by the configuration of microcracks. The inherent statistical distribution of the fracture toughness arising from the intrinsic random nature of microcracks is explored using a statistical approach. For this purpose, a computer simulation model is introduced, which incorporates a statistical characterization of geometrical parameters of a random microcrack array.

Lua, Yuan J.↗

NASA NTRS2006

Performance of Two Cloud-Radiation Parameterization Schemes in the Finite Volume General Circulation Model for Anomalously Wet May and June 2003 Over the Continental United States and Amazonia

An objective assessment of the impact of a new cloud scheme, called Microphysics of Clouds with Relaxed Arakawa-Schubert Scheme (McRAS) (together with its radiation modules), on the finite volume general circulation model (fvGCM) was made with a set of ensemble forecasts that invoke performance evaluation over both weather and climate timescales. The performance of McRAS (and its radiation modules) was compared with that of the National Center for Atmospheric Research Community Climate Model (NCAR CCM3) cloud scheme (with its NCAR physics radiation). We specifically chose the boreal summer months of May and June 2003, which were characterized by an anomalously wet eastern half of the continental United States as well as northern regions of Amazonia. The evaluation employed an ensemble of 70 daily 10-day forecasts covering the 61 days of the study period. Each forecast was started from the analyzed initial state of the atmosphere and spun-up soil moisture from the first-day forecasts with the model. Monthly statistics of these forecasts with up to 10-day lead time provided a robust estimate of the behavior of the simulated monthly rainfall anomalies. Patterns of simulated versus observed rainfall, 500-hPa heights, and top-of-the-atmosphere net radiation were recast into regional anomaly correlations. The correlations were compared among the simulations with each of the schemes. The results show that fvGCM with McRAS and its radiation package performed discernibly better than the original fvGCM with CCM3 cloud physics plus its radiation package. The McRAS cloud scheme also showed a reasonably positive response to the observed sea surface temperature on mean monthly rainfall fields at different time leads. This analysis represents a method for helpful systematic evaluation prior to selection of a new scheme in a global model.

Sud, Y. C.↗

DOE OSTI2023

Full Counting Statistics of Charge in Chaotic Many-Body Quantum Systems

We investigate the full counting statistics of charge transport in U(1)-symmetric random unitary circuits. We consider an initial mixed state prepared with a chemical potential imbalance between the left and right halves of the system and study the fluctuations of the charge transferred across the central bond in typical circuits. Using an effective replica statistical mechanics model and a mapping onto an emergent classical stochastic process valid at large on-site Hilbert space dimension, we show that charge transfer fluctuations approach those of the symmetric exclusion process at long times, with subleading t –1/2 quantum corrections. Here, we discuss our results in the context of fluctuating hydrodynamics and macroscopic fluctuation theory of classical nonequilibrium systems and check our predictions against direct matrix-product state calculations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS2004

The Contribution of Dynamic Interannual Variability to Ozone Trends

At middle latitudes the total column ozone and the lower stratospheric ozone mixing ratio exhibit natural variability. Models and statistical analyses of observations such as SAGE ozone profiles and TOMS column measurements show that seasonal cycle, solar cycle, and interannual dynamical variability and dynamical phenomena such as the quasi- biennial oscillation all contribute to ozone variability. These must be accounted for when deriving ozone trends. Systematic or random changes in the atmospheric circulation may also contribute to ozone trends. It is presently unclear how much of the ozone trend derived from observations is due to changes in the chemical composition of the stratosphere and how much is due to changes in the atmospheric circulation. We are attempting to resolve this issue by comparing a twenty-five year simulation of ozone with fixed source gas boundary conditions with an identical simulation with time dependent source gas boundary conditions. Both simulations are driven with output from a general circulation model that produces realistic interannual variability in dynamical forcing. The model trend in ozone due to changes in composition is determined from the difference in these simulations. We compare these trends with trends determined from observations and model output using the same analysis techniques. Initial results emphasize the complications to attribution of observed ozone trends to dynamical and photochemical effects that are due to interrelationships between trends in transport, temperature, and photochemical effects. It may not be possible to describe the ozone trend as a superposition of dynamical and photochemical contributions.

Douglass, Anne↗

DOE OSTI2021

Cosmic shear cosmology beyond two-point statistics: a combined peak count and correlation function analysis of DES-Y1

ABSTRACT We constrain cosmological parameters from a joint cosmic shear analysis of peak-counts and the two-point shear correlation functions, as measured from the Dark Energy Survey (DES-Y1). We find the structure growth parameter $S_8\equiv \sigma _8\sqrt{\Omega _{\rm m}/0.3} = 0.766^{+0.033}_{-0.038}$ which, at 4.8 per cent precision, provides one of the tightest constraints on S8 from the DES-Y1 weak lensing data. In our simulation-based method we determine the expected DES-Y1 peak-count signal for a range of cosmologies sampled in four w cold dark matter parameters (Ωm, σ8, h, w0). We also determine the joint covariance matrix with over 1000 realizations at our fiducial cosmology. With mock DES-Y1 data we calibrate the impact of photometric redshift and shear calibration uncertainty on the peak-count, marginalizing over these uncertainties in our cosmological analysis. Using dedicated training samples we show that our measurements are unaffected by mass resolution limits in the simulation, and that our constraints are robust against uncertainty in the effect of baryon feedback. Accurate modelling for the impact of intrinsic alignments on the tomographic peak-count remains a challenge, currently limiting our exploitation of cross-correlated peak counts between high and low redshift bins. We demonstrate that once calibrated, a fully tomographic joint peak-count and correlation functions analysis has the potential to reach a 3 per cent precision on S8 for DES-Y1. Our methodology can be adopted to model any statistic that is sensitive to the non-Gaussian information encoded in the shear field. In order to accelerate the development of these beyond-two-point cosmic shear studies, our simulations are made available to the community upon request.

Harnois-Déraps, Joachim↗

NASA NTRS2011

On the Causes and Dynamics of the Early Twentieth Century North American Pluvial

The early twentieth century North American pluvial (1905-1917) was one of the most extreme wet periods of the last five hundred years and directly led to overly generous water allotments in the water-limited American West. Here we examine the causes and dynamics of the pluvial event using a combination of observation-based data sets and general circulation model (GCM) experiments. The character of the moisture surpluses during the pluvial differed by region, alternately driven by increased precipitation (the Southwest), low evaporation from cool temperatures (the Central Plains), or a combination of the two (the Pacific Northwest). Cool temperature anomalies covered much of the west and persisted through most months, part of a globally extensive period of cooler land and sea surface temperatures (SST). Circulation during boreal winter favored increased moisture import and precipitation in the southwest, while other regions and seasons were characterized by near normal or reduced precipitation. Anomalies in the mean circulation, precipitation, and SST fields are partially consistent with the relatively weak El Nino forcing during the pluvial, and also reflect the impact of positive departures in the Arctic Oscillation that occurred in ten of the thirteen pluvial winters. Differences between the reanalysis dataset, an independent statistical drought model, and GCM simulations highlight some of the remaining uncertainties in understanding the full extent of SST forcing of North American hydroclimatic variability.

Cook, Benjamin I.↗

NASA NTRS1995

Comment on 'The Global Resurfacing of Venus' by R. G. Strom, G.G. Schaber, and D.D. Dawson

The distribution of impact craters on Venus has been the subject of a great deal of analysis since the return of Magellan data. Phillips el al. (1992) performed Monte Carlo two-dimensional (2-D) modeling of the areal distribution of craters, and the results of that exercise allowed a restricted, but still quite large, range of possible planetary resurfacing histories, including the possibility that the crater, were emplaced on a geologically inactive planet. However, the nonrandom distribution of embayed and deformed craters (Phillips el al., 1992), the hypsometric distribution of craters (Herrick and Phillips, 1994), the varied degradation states of craters (Izenberg et al., 1994), their nonrandom distribution with different geologic terrain types (Namiki and Solomon, 1994; Price et al, 1994), and three-dimensional resurfacing modeling (Bullock el al., 1993) all seem to argue against that particular possibility. In contrast, Strom el al. (1994) have collected a refined and more comprehensive data set of impact features, and they input these data into more sophisticated 2-D Monte Carlo modeling and statistical analyses of the areal distribution of craters, the hypsometric distribution of craters, and the number of embayed craters. They concluded that 'Venus experienced a global resurfacing event about 300 m.y. ago followed by a dramatic reduction of volcanism and tectonism. This global resurfacing event ended abruptly (less than 10 m.y.). The present crater population has accumulated since then and remains largely intact . . . only about 4%-6% of the planet has been volcanically resurfaced since the global event . . .' If these conclusions are well-founded, this work certainly represents a significant advancement in restricting tile number of plausible resurfacing histories for the planet. If Strom et al. (1994) are correct, it would also mean that all of the other aforementioned works are in error to various degrees, or at least represent overzealous interpretation of the data. However, we have identified apparent flaws in the observations, modeling, and interpretations presented by Strom el al. (1994) that lead us to question whether their conclusions are warranted. We limit our comments to three areas of their analysis: (1) observations pertaining to the number and area of disrupted and pristine craters and crater-related features, (2) modeling of the areal and elevation distribution of craters, and (3) interpretations of resurfacing models.

Herrick, Robert R.↗

DOE OSTI2026

Model Parameter Development for Complex Materials: Species-Specific Diffusion Barriers in 316 Stainless Steel from Systematic DFT Calculations

Vacancy-mediated diffusion barriers in 316 stainless steel have been systematically calculated using density functional theory to provide essential parameters for mesoscale microstructure evolution models. A statistical sampling approach employing 210 nudged elastic band calculations across multiple special quasi-random structures captures the effects of local chemical environments in this concentrated alloy. The computational methodology addresses challenges specific to chemically disordered systems, including proper magnetic treatment throughout multi-step calculations and validation against experimental structural properties. The calculated activation barriers reveal clear species-dependent diffusion behavior with the hierarchy Ni >> Fe ˜ Cr >> Mo. Nickel exhibits the highest barriers (0.74–1.31 eV, mean 1.045 eV), confirming its role as the slowest-diffusing major component. Iron and chromium show similar moderate barriers averaging 0.587 eV and 0.522 eV, respectively. Remarkably, molybdenum demonstrates exceptionally low barriers (0.12–0.28 eV, mean 0.194 eV), suggesting much higher mobility than previously recognized and potentially significant implications for precipitation kinetics and microstructure evolution. The barrier ranges remain consistent across different 316 SS compositions, supporting parameter transferability for modeling applications. The overall mean barrier of 0.64 eV provides a practical approximation for phase field simulations, while species-specific values enable detailed treatments of diffusion-controlled processes. This systematic approach establishes a validated framework for generating diffusion parameters in other concentrated alloys where experimental data are limited, while providing the first systematic set of species-specific barriers for predictive modeling of 316 stainless steel microstructure evolution.

36 MATERIALS SCIENCE↗

NASA NTRS1989

Periodicity in extinction and the problem of catastrophism in the history of life

The hypothesis that extinction events have recurred periodically over the last quarter billion years is greatly strengthened by new data on the stratigraphic ranges of marine animal genera. In the interval from the Permian to Recent, these data encompass some 13,000 generic extinctions, providing a more sensitive indicator of species-level extinctions than previously used familial data. Extinction time series computed from the generic data display nine strong peaks that are nearly uniformly spaced at 26 Ma intervals over the last 270 Ma. Most of these peaks correspond to extinction events recognized in more detailed, if limited, biostratigraphic studies. These new data weaken or negate most arguments against periodicity, which have involved criticisms of the taxonomic data base, sampling intervals, chronometric time scales, and statistical methods used in previous analyses. The criticisms are reviewed in some detail and various new calculations and simulations, including one assessing the effects of paraphyletic taxa, are presented. Although the new data strengthen the case for periodicity, they offer little new insight into the deriving mechanism behind the pattern. However, they do suggest that many of the periodic events may not have been catastrophic, occurring instead over several stratigraphic stages or substages.

NASA Discipline Number 52-40↗

DOE OSTI2025

A Computational Workflow of Elucidating Viral Impact on Mediating Microbial Response to In-situ Experimental Warming: Bridging microbial modeling to carbon and mineral modeling

Viruses are abundant in soils and shape microbial communities in ways that can potentially influence ecosystem processes, yet their contributions to carbon cycling and mineral transformations remain poorly understood. Here we present a multi-phase framework that links virus-host interactions to soil biogeochemistry by combining ecological simulations, genome- and community-scale metabolic modeling, and statistical and machine-learning analyses. We first calibrated microbial abundance profiles under explicit infection scenarios to capture how viral pressure alters community structure, then explored alternative interaction strategies, including kill-the-winner, piggyback-the-winner, and mixed lytic-lysogenic modes, through forward simulations. These ecological shifts were translated into metabolic consequences using exchange fluxes summarized into biologically meaningful categories, while integrated statistical and machine-learning screens elevated subtle but consistent signals. Application of this framework revealed that viral infections shift the balance between organic and inorganic fluxes, redirecting metabolism from diffuse organic transformations toward inorganic pools such as protons and CO 2 , directly linking viral regulation to respiration and soil carbon balance. The roll-up analysis also isolated perturbations in critical mineral ions, including magnesium, manganese, zinc, and copper, which serve as essential enzymatic cofactors. In piggyback-the-winner scenarios, uptake of these ions was strongly suppressed. Contrasting viral strategies produced distinct community structures and metabolic outcomes, from broad suppression under kill-the-winner dynamics to dramatic redistributions under high-lytic and high-gain lysogenic regimes that collapsed vulnerable microbial populations while promoting opportunists. Together, these results provide a tractable path to trace viral perturbations from host abundance shifts to metabolic flux adjustments and ecosystem-scale processes, offering a practical way to include viruses in earth system models.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2000

An architecture for a brain-image database

The widespread availability of methods for noninvasive assessment of brain structure has enabled researchers to investigate neuroimaging correlates of normal aging, cerebrovascular disease, and other processes; we designate such studies as image-based clinical trials (IBCTs). We propose an architecture for a brain-image database, which integrates image processing and statistical operators, and thus supports the implementation and analysis of IBCTs. The implementation of this architecture is described and results from the analysis of image and clinical data from two IBCTs are presented. We expect that systems such as this will play a central role in the management and analysis of complex research data sets.

NASA Discipline Neuroscience↗

NASA NTRS1994

Radiative interactions in multi-dimensional chemically reacting flows using Monte Carlo simulations

The Monte Carlo method (MCM) is applied to analyze radiative heat transfer in nongray gases. The nongray model employed is based on the statistical narrow band model with an exponential-tailed inverse intensity distribution. The amount and transfer of the emitted radiative energy in a finite volume element within a medium are considered in an exact manner. The spectral correlation between transmittances of two different segments of the same path in a medium makes the statistical relationship different from the conventional relationship, which only provides the non-correlated results for nongray methods is discussed. Validation of the Monte Carlo formulations is conducted by comparing results of this method of other solutions. In order to further establish the validity of the MCM, a relatively simple problem of radiative interactions in laminar parallel plate flows is considered. One-dimensional correlated Monte Carlo formulations are applied to investigate radiative heat transfer. The nongray Monte Carlo solutions are also obtained for the same problem and they also essentially match the available analytical solutions. the exact correlated and non-correlated Monte Carlo formulations are very complicated for multi-dimensional systems. However, by introducing the assumption of an infinitesimal volume element, the approximate correlated and non-correlated formulations are obtained which are much simpler than the exact formulations. Consideration of different problems and comparison of different solutions reveal that the approximate and exact correlated solutions agree very well, and so do the approximate and exact non-correlated solutions. However, the two non-correlated solutions have no physical meaning because they significantly differ from the correlated solutions. An accurate prediction of radiative heat transfer in any nongray and multi-dimensional system is possible by using the approximate correlated formulations. Radiative interactions are investigated in chemically reacting compressible flows of premixed hydrogen and air in an expanding nozzle. The governing equations are based on the fully elliptic Navier-Stokes equations. Chemical reaction mechanisms were described by a finite rate chemistry model. The correlated Monte Carlo method developed earlier was employed to simulate multi-dimensional radiative heat transfer. Results obtained demonstrate that radiative effects on the flowfield are minimal but radiative effects on the wall heat transfer are significant. Extensive parametric studies are conducted to investigate the effects of equivalence ratio, wall temperature, inlet flow temperature, and nozzle size on the radiative and conductive wall fluxes.

Liu, Jiwen↗

DOE OSTI2022

Numerical and experimental characterization of twin transmission across grain boundaries along the forward and lateral directions

Pervasive deformation twinning and transmission events across grain boundaries (GBs) affect the strength and failure of hexagonal close-packed (HCP) magnesium. A three-dimensional twin can transmit along the twinning shear direction, (forward), and along the direction perpendicular to both the twinning plane normal and the shear direction, (lateral). For the first time, phase-field calculations and electron backscatter diffraction (EBSD)-based statistical analysis are combined to investigate the effect of the twinned grain boundary characteristics on twin transmission (TT) along the forward and lateral directions. This combined analysis reveals that TT propensity decreases with increasing misorientation angle across the GB for both forward and lateral directions. Also, the TT is more favorable along the lateral than along the forward direction. Twin transmission seems harder across GBs with a misorientation axis closer to the twin -direction than the other directions ( and ). Further, the PF calculations reveal that, at the onset of a transmission process, the crystallography tends to be preserved in the case of lateral transmission, whereas, in the forward transmission case, the transmitted twin punches straight through the GBs and its morphology prevails. Additionally, the EBSD analysis finds that pure forward and lateral transmissions do not occur often, yet reveals a preference for lateral propagation consistent with PF simulations. Further, the local twin transmission configurations observed in the actual material do not correspond to pure tilt or twist GBs, which are most commonly considered as model cases.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Extracting model-independent nuclear level densities away from stability

The nuclear level density (NLD) is a fundamental measure of the complex structure of atomic nuclei at relatively high energies. Here, in this study, we present the first model-independent measurement of the absolute partial NLD for a short-lived nucleus. For this purpose we adapt the recently introduced “shape method” for β-decay experiments, providing the shape of the γ-ray strength function for exotic nuclei. In this work, we show that combining the shape method with the β-Oslo technique allows for the extraction of the NLD of the populated states without the need for theoretical input. This development opens the way for the extraction of experimental NLDs far from stability with major implications in astrophysical and other applications. We benchmark our approach using data for the stable 76 Ge nucleus, finding excellent agreement with previous experimental results. In addition, we present new experimental data and determine the absolute partial level density for the short-lived 88 Kr nucleus. Our results suggest a fivefold increase in the NLD for the case of 88 Kr, compared to the recommended values from semimicroscopic Hartree-Fock Bogoliubov calculations recommended by the RIPL3 nuclear data library. However, our results are in good agreement with other semimicroscopic level density models. We demonstrate the impact of our method on the 87 Kr(n, γ) neutron capture rate and show that our experimental uncertainties for NLDs fulfill the requirements needed for astrophysical calculations predicting r-process abundances.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Spin cutoff factor and level density for Ni 59 from an analysis of compound nuclear reactions

The spin cutoff parameter for 59 Ni has been studied from different types of experimental data, including neutron angular distributions from the 56 Fe(α,n) 59 Ni reaction, spin of discrete levels, the level density from the proton evaporation spectrum of the 54 Fe( 6 Li,p) 59 Ni reaction and neutron resonance spacing. Experimental data points were compared with calculations using models widely used in literature. It was found that the available empirical models overestimate data points in the energy region below the neutron separation energy, while microscopic calculations which take into account pairing correlations within Hartree-Fock +Bardeen-Cooper-Schrieffer approach are consistent with data. As a result, it confirmed earlier findings that pairing correlations play an important role and need to be taken into account when the spin cutoff parameter is calculated below the neutron separation energy.

(a,n) reaction↗

NASA NTRS2020

Using Two-Stream Theory to Capture Fluctuations of Satellite-Perceived TOA SW Radiances Reflected from Clouds Over Ocean

Shortwave (SW) fluxes estimated from broadband radiometry rely on empirically gathered and hemispherically resolved fields of outgoing top-of-atmosphere (TOA) radiances. This study aims to provide more accurate and precise fields of TOA SW radiances reflected from clouds over ocean by introducing a novel semiphysical model predicting radiances per narrow sun-observer geometry. This model was statistically trained using CERES-measured radiances paired with MODIS-retrieved cloud parameters as well as reanalysis-based geophysical parameters. By using radiative transfer approximations as a framework to ingest the above parameters, the new approach incorporates cloud-top effective radius and above-cloud water vapor in addition to traditionally used cloud optical depth, cloud fraction, cloud phase, and surface wind speed. A two-stream cloud albedo – serving to statistically incorporate cloud optical thickness and cloud-top effective radius – and Cox–Munk ocean reflectance were used to describe an albedo over each CERES footprint. Effective-radius-dependent asymmetry parameters were obtained empirically and separately for each viewing-illumination geometry. A simple equation of radiative transfer, with this albedo and attenuating above-cloud water vapor as inputs, was used in its log-linear form to allow for statistical optimization. We identified the two-stream functional form that minimized radiance residuals calculated against CERES observations and outperformed the state-of-the-art approach for most observer geometries outside the sun-glint and solar zenith angles between 20 and 70∘, reducing the median SD of radiance residuals per solar geometry by up to 13.2 % for liquid clouds, 1.9 % for ice clouds, and 35.8 % for footprints containing both cloud phases. Geometries affected by sun glint (constituting between 10 % and 1 % of the discretized upward hemisphere for solar zenith angles of 20 and 70∘, respectively), however, often showed weaker performance when handled with the new approach and had increased residuals by as much as 60 % compared to the state-of-the-art approach. Overall, uncertainties were reduced for liquid-phase and mixed-phase footprints by 5.76 % and 10.81 %, respectively, while uncertainties for ice-phase footprints increased by 0.34 %. Tested for a variety of scenes, we further demonstrated the plausibility of scene-wise predicted radiance fields. This new approach may prove useful when employed in angular distribution models and may result in improved flux estimates, in particular dealing with clouds characterized by small or large droplet/crystal sizes.

Two- stream theory↗

DOE OSTI2023

Component-Level Analysis of Heating and Cooling Loads in the U.S. Residential Building Stock

The residential building sector accounts for a substantial portion of total energy consumption in the United States and offers a significant opportunity for energy reduction and decarbonization through improvements in energy efficiency. Heating and air conditioning are the primary contributors to residential energy usage and electricity system peak demand. However, due to the diversity of the housing stock and the complexity of factors affecting heating and cooling demand, identifying the relative contributions to heating and cooling loads poses challenges. To address this, we applied the ResStock analysis tool to simulate 550,000 building energy models, providing statistical representation of residential buildings in the contiguous United States. We introduced outputs that quantified the heating and cooling influence of different components of a home, such as air leakage, envelope components (ceilings, walls, windows, foundations), internal heat gains from people, lighting, plug loads, and duct losses and gains. Leveraging the granularity of ResStock, we present a dataset to enable deeper understanding of the contributors to heating and cooling loads as a function of housing characteristics such as location, envelope efficiency, and building type. This work aims to support prioritization of research and development and informed decision-making for residential building decarbonization.

building simulation↗

DOE OSTI2024

Finite density QCD equation of state: Critical point and lattice-based 𝑇′ expansion

Here, we present a novel construction of the QCD equation of state (EoS) at finite baryon density. Our work combines a recently proposed resummation scheme for lattice QCD results with the universal critical behavior at the QCD critical point. This allows us to obtain a family of equations of state in the range 0 ≤ 𝜇 𝐵 ≤ 700 MeV and 25 MeV ≤ 𝑇 ≤ 800 MeV, which match lattice QCD results near 𝜇 𝐵 =0 while featuring a critical point in the 3D Ising model universality class. The position of the critical point can be chosen within the range accessible to beam-energy scan heavy-ion collision experiments. The strength of the singularity and the shape of the critical region are parametrized using a standard parameter set. We impose stability and causality constraints and discuss the available ranges of critical point parameter choices, finding that they extend beyond earlier parametric QCD EoS proposals. We present thermodynamic observables, including baryon density, pressure, entropy density, energy density, baryon susceptibility and speed of sound, that cover a wide range in the QCD phase diagram relevant for experimental exploration.

Astronomy & Astrophysics↗