Search NASA⌕ Search

SEARCH · Search NASA

Results for “Constrained parameter spaces”

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 415 records · Page 23

FIGS: spectral fitting constraints on the star formation history of massive galaxies since the cosmic noon

We constrain the stellar population properties of a sample of 52 massive galaxies – with stellar mass log (M(s)/M(⊙)) ≳ 10.5 – over the redshift range 0.5 < z < 2 by use of observer-frame optical and near-infrared slitless spectra from Hubble Space Telescope’s ACS and WFC3 grisms. The deep exposures (∼100 ks) allow us to target individual spectra of massive galaxies to F160W = 22.5 AB. Our spectral fitting approach uses a set of six base models adapted to the redshift and spectral resolution of each observation, and fits the weights of the base models, including potential dust attenuation, via a Markov Chain Monte Carlo method. Our sample comprises a mixed distribution of quiescent (19) and star-forming galaxies (33). We quantify the width of the age distribution (Δt) that is found to dominate the variance of the retrieved parameters according to principal component analysis. The population parameters follow the expected trend towards older ages with increasing mass, and Δt appears to weakly anticorrelate with stellar mass, suggesting a more efficient star formation at the massive end. As expected, the redshift dependence of the relative stellar age (measured in units of the age of the Universe at the source) in the quiescent sample rejects the hypothesis of a single burst (aka monolithic collapse). Radial colour gradients within each galaxy are also explored, finding a wider scatter in the star-forming subsample, but no conclusive trend with respect to the population parameters.

Ignacio Ferreras↗

High Power Laser Diode Arrays for 2-Micron Solid State Coherent Lidars Applications

Laser diode arrays are critical components of any diode-pumped solid state laser systems, constraining their performance and reliability. Laser diode arrays (LDAs) are used as the pump source for energizing the solid state lasing media to generate an intense coherent laser beam with a high spatial and spectral quality. The solid state laser design and the characteristics of its lasing materials define the operating wavelength, pulse duration, and power of the laser diodes. The pump requirements for high pulse energy 2-micron solid state lasers are substantially different from those of more widely used 1-micron lasers and in many aspects more challenging [1]. Furthermore, the reliability and lifetime demanded by many coherent lidar applications, such as global wind profiling from space and long-range clear air turbulence detection from aircraft, are beyond the capability of currently available LDAs. In addition to the need for more reliable LDAs with longer lifetime, further improvement in the operational parameters of high power quasi-cw LDAs, such as electrical efficiency, brightness, and duty cycle, are also necessary for developing cost-effective 2-micron coherent lidar systems for applications that impose stringent size, heat dissipation, and power constraints. Global wind sounding from space is one of such applications, which is the main driver for this work as part of NASA s Laser Risk Reduction Program. This paper discusses the current state of the 792 nm LDA technology and the technology areas being pursued toward improving their performance. The design and development of a unique characterization facility for addressing the specific issues associated with the LDAs for pumping 2-micron coherent lidar transmitters and identifying areas of technological improvement will be described. Finally, the results of measurements to date on various standard laser diode packages, as well as custom-designed packages with potentially longer lifetime, will be reported.

Amzajerdian, Farzin↗

Modeling the Spin Motor Current of the International Space Station's Control Moment Gyroscopes

The International Space Station (ISS) attitude control is provided by two means: The Russian Segment uses thrusters and the U.S. Segment uses double-gimbaled control moment gyroscopes (CMG). CMGs are used as momentum exchange devices, providing non propulsive attitude control for the vehicle. The CMGs are very important for the ISS program because, first, they save propellant - which needs to be transferred to the Station in special cargo vehicles - and, second, they provide the microgravity environment on the Station - which is necessary for scientific experiments planned for the ISS mission. Since 2002, when one of the CMG on the ISS failed, all CMGs are closely monitored. High gimbal rates, vibration spikes, unusual variations of spin motor current and bearing temperatures are of great concern, since these parameters are the CMG health indicators. The telemetry analysis of these and some other CMG parameters is used to determine constrains and make changes to the CMGs operation on board. These CMG limitations, in turn, may limit the ISS attitude control capabilities and may be critical to ISS operation. Therefore, it is important to know whether the CMG parameter is nominal or out of family, and why. The goal of this project is to analyze an important CMG parameter - spin motor current. Some operational decisions are made now based on the spin motor current signatures. The spin motor current depends on gimbal rates, ISS rates, and spin bearing friction. The spin bearing friction in turn depends on the bearing temperatures, wheel rates, normal load - which is a function of gimbal and wheel rates - lubrication, etc. The first task of this project is to create a spin motor current mathematical model based on CMG dynamics model and the current knowledge on bearing friction in microgravity.

Pereira, Miguel A.↗

Lunar Neutral Exposphere Properties from Pickup Ion Analysis

Composition and structure of neutral constituents in the lunar exosphere can be determined through measurements of phase space distributions of pickup ions borne from the exosphere [1]. An essential point made in an early study [ 1 ] and inferred by recent pickup ion measurements [2, 3] is that much lower neutral exosphere densities can be derived from ion mass spectrometer measurements of pickup ions than can be determined by conventional neutral mass spectrometers or remote sensing instruments. One approach for deriving properties of neutral exospheric source gasses is to first compare observed ion spectra with pickup ion model phase space distributions. Neutral exosphere properties are then inferred by adjusting exosphere model parameters to obtain the best fit between the resulting model pickup ion distributions and the observed ion spectra. Adopting this path, we obtain ion distributions from a new general pickup ion model, an extension of a simpler analytic description obtained from the Vlasov equation with an ion source [4]. In turn, the ion source is formed from a three-dimensional exospheric density distribution, which can range from the classical Chamberlain type distribution to one with variable exobase temperatures and nonthermal constituents as well as those empirically derived. The initial stage of this approach uses the Moon's known neutral He and Na exospheres to deriv e He+ and Na+ pickup ion exospheres, including their phase space distributions, densities and fluxes. The neutral exospheres used are those based on existing models and remote sensing studies. As mentioned, future ion measurements can be used to constrain the pickup ion model and subsequently improve the neutral exosphere descriptions. The pickup ion model is also used to estimate the exosphere sources of recently observed pickup ions on KAGUYA [3]. Future missions carrying ion spectrometers (e.g., ARTEMIS) will be able to study the lunar neutral exosphere with great sensitivity, yielding the necessary ion velocity spectra needed to further analysis of parent neutral exosphere properties.

Hartle, R. E.↗

Potential Multi-Component Structure of the Debris Disk Around HIP 17439 Revealed by Herschel DUNES

Context. The dust observed in debris disks is produced through collisions of larger bodies left over from the planet/planetesimal formation process. Spatially resolving these disks permits to constrain their architecture and thus that of the underlying planetary/planetesimal system. Aims. Our Herschel open time key program DUNES aims at detecting and characterizing debris disks around nearby, sun-like stars. In addition to the statistical analysis of the data, the detailed study of single objects through spatially resolving the disk and detailed modeling of the data is a main goal of the project. Methods. We obtained the first observations spatially resolving the debris disk around the sun-like star HIP 17439 (HD 23484) using the instruments PACS and SPIRE on board the Herschel Space Observatory. Simultaneous multi-wavelength modeling of these data together with ancillary data from the literature is presented. Results. A standard single component disk model fails to reproduce the major axis radial profiles at 70 μm, 100 μm, and 160 μm simultaneously. Moreover, the best-fit parameters derived from such a model suggest a very broad disk extending from few au up to few hundreds of au from the star with a nearly constant surface density which seems physically unlikely. However, the constraints from both the data and our limited theoretical investigation are not strong enough to completely rule out this model. An alternative, more plausible, and better fitting model of the system consists of two rings of dust at approx. 30 au and 90 au, respectively, while the constraints on the parameters of this model are weak due to its complexity and intrinsic degeneracies. Conclusions. The disk is probably composed of at least two components with different spatial locations (but not necessarily detached), while a single, broad disk is possible, but less likely. The two spatially well-separated rings of dust in our best-fit model suggest the presence of at least one high mass planet or several low-mass planets clearing the region between the two rings from planetesimals and dust.

sun-like stars↗

Lightning Charge Retrievals: Dimensional Reduction, LDAR Constraints, and a First Comparison with LIS Satellite Data

A "dimensional reduction" (DR) method is introduced for analyzing lightning field changes (DELTAEs) whereby the number of unknowns in a discrete two-charge model is reduced from the standard eight (x, y, z, Q, x', y', z', Q') to just four (x, y, z, Q). The four unknowns (x, y, z, Q) are found by performing a numerical minimization of a chi-square function. At each step of the minimization, an Overdetermined Fixed Matrix (OFM) method is used to immediately retrieve the best "residual source" (x', y', z', Q'), given the values of (x, y, z, Q). In this way, all 8 parameters (x, y, z, Q, x', y', z', Q') are found, yet a numerical search of only 4 parameters (x, y, z, Q) is required. The DR method has been used to analyze lightning-caused DeltaEs derived from multiple ground-based electric field measurements at the NASA Kennedy Space Center (KSC) and USAF Eastern Range (ER). The accuracy of the DR method has been assessed by comparing retrievals with data provided by the Lightning Detection And Ranging (LDAR) system at the KSC-ER, and from least squares error estimation theory, and the method is shown to be a useful "stand-alone" charge retrieval tool. Since more than one charge distribution describes a finite set of DELTAEs (i.e., solutions are non-unique), and since there can exist appreciable differences in the physical characteristics of these solutions, not all DR solutions are physically acceptable. Hence, an alternative and more accurate method of analysis is introduced that uses LDAR data to constrain the geometry of the charge solutions, thereby removing physically unacceptable retrievals. The charge solutions derived from this method are shown to compare well with independent satellite- and ground-based observations of lightning in several Florida storms.

Koshak, W. J.↗

Constraining the phase shift of relativistic species in DESI BAOs

In the early Universe, neutrinos decouple quickly from the primordial plasma and propagate without further interactions. The impact of free-streaming neutrinos is to create a temporal shift in the gravitational potential that impacts the acoustic waves known as baryon acoustic oscillations (BAOs), resulting in a non-linear spatial shift in the Fourier-space BAO signal. In this work, we make use of and extend upon an existing methodology to measure the phase shift amplitude $\beta _{\phi }$ and apply it to the Dark Energy Spectroscopic Instrument (DESI) Data Release 1 (DR1) BAOs with an anisotropic BAO fitting pipeline. We validate the fitting methodology by testing the pipeline with two publicly available fitting codes applied to highly precise cubic box simulations and realistic simulations representative of the DESI DR1 data. We find further study towards the methods used in fitting the BAO signal will be necessary to ensure accurate constraints on $\beta _{\phi }$ in future DESI data releases. Using DESI DR1, we present individual measurements of the anisotropic BAO distortion parameters and the $\beta _{\phi }$ for the different tracers, and additionally a combined fit to $\beta _{\phi }$ resulting in $\beta _{\phi } = 2.7 \pm 1.7$. After including a prior on the distortion parameters from constraints using Planck we find $\beta _{\phi } = 2.7^{+0.60}_{-0.67}$ suggesting $\beta _{\phi } > 0$ at 4.3$\sigma$ significance. This result may hint at a phase shift that is not purely sourced from the standard model expectation for $N_{\rm {eff}}$ or could be a upwards statistical fluctuation in the measured $\beta _{\phi }$; this result relaxes in models with additional freedom beyond Lambda-cold dark matter.

79 ASTRONOMY AND ASTROPHYSICS↗

A physics-constrained neural network for multiphase flows

The present study develops a physics-constrained neural network (PCNN) to predict sequential patterns and motions of multiphase flows (MPFs), which includes strong interactions among various fluid phases. To predict the order parameters, which locate individual phases in the future time, a neural network (NN) is applied to quickly infer the dynamics of the phases by encoding observations. The multiphase consistent and conservative boundedness mapping algorithm (MCBOM) is next implemented to correct the predicted order parameters. This enforces the predicted order parameters to strictly satisfy the mass conservation, the summation of the volume fractions of the phases to be unity, the consistency of reduction, and the boundedness of the order parameters. Then, the density of the fluid mixture is updated from the corrected order parameters. Finally, the velocity in the future time is predicted by another NN with the same network structure, but the conservation of momentum is included in the loss function to shrink the parameter space. The proposed PCNN for MPFs sequentially performs (NN)-(MCBOM)-(NN), which avoids nonphysical behaviors of the order parameters, accelerates the convergence, and requires fewer data to make predictions. Numerical experiments demonstrate that the proposed PCNN is capable of predicting MPFs effectively.

Mechanics↗

Lightning Charge Retrievals: Dimensional Reduction, LDAR Constraints, and a First Comparison with LIS Satellite Data

A "dimensional reduction" ("DR") method is introduced for analyzing lightning field changes ((Delta)Es) whereby the number of unknowns in a discrete two-charge model is reduced from the standard eight (x, y, z, Q, x', y', z', Q') to just four (x, y, z, Q). The four unknowns (x, y, z, Q) are found by performing a numerical minimization of a chi-square function. At each step of the minimization, an overdetermined fixed matrix (OFM) method is used to immediately retrieve the best "residual source" (x', y', z', Q'), given the values of (x, y, z, Q). In this way, all eight parameters Ix. y, z, Q, x', y', z', Q') are found, yet a numerical search of only four parameters (x, y, z, Q) is required. The DR method has been used to analyze lightning-caused (Delta)Es derived from multiple ground-based electric field measurements at the NASA Kennedy Space Center (KSC) and U.S. Air Force Eastern Range (ER). The accuracy of the DR method has been assessed by comparing retrievals with data provided by the lightning detection and ranging (LDAR) system at the KSC-ER. and from least squares error estimation theory, and the method is shown to be a useful "stand alone" charge retrieval tool. Since more than one charge distribution describes a finite set of (Delta)Es (i.e., solutions are nonunique), and since there can be appreciable differences in the physical characteristics of these solutions, not all DR solutions are physically acceptable. Hence. an alternative and more accurate method of analysis is introduced that uses LDAR data to constrain the geometry of the charge solutions. thereby removing physically unacceptable retrievals. The charge solutions derived from this method are shown to compare well with independent satellite- and ground-based observations of lightning m several Florida storms.

Koshak, W. J.↗

A physics-constrained deep learning treatment of runaway electron dynamics

An adjoint formulation leveraging a physics-informed neural network (PINN) is employed to advance the density moment of a runaway electron (RE) distribution forward in time. A distinguishing feature of this approach is that once the adjoint problem is solved, its solution can be used to project the RE density forward in time for an arbitrary initial momentum space distribution of REs. Furthermore, by employing a PINN, a parametric solution to the adjoint problem can be learned. Thus, once trained, this adjoint-deep learning framework is able to efficiently project the RE density forward in time across various plasma conditions while still including a fully kinetic description of RE dynamics. As an example application, the temporal evolution of the density of primary electrons is studied, with particular emphasis on evaluating the decay of a RE population when below threshold. Predictions from the adjoint-deep learning framework are found to be in good agreement with a traditional relativistic electron Fokker–Planck solver, for several distinct initial conditions, and across an array of physics parameters. Once trained, the PINN thus provides a means of generating RE density time histories with exceptionally low online execution time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Cosmic topology. Part I. Limits on orientable Euclidean manifolds from circle searches

The Einstein field equations of general relativity constrain the local curvature at every point in spacetime, but say nothing about the global topology of the Universe. Cosmic microwave background anisotropies have proven to be the most powerful probe of non-trivial topology since, within ΛCDM, these anisotropies have well-characterized statistical properties, the signal is principally from a thin spherical shell centered on the observer (the last scattering surface), and space-based observations nearly cover the full sky. The most generic signature of cosmic topology in the microwave background is pairs of circles with matching temperature and polarization patterns. No such circle pairs have been seen above noise in the WMAP or Planck temperature data, implying that the shortest non-contractible loop around the Universe through our location is longer than 98.5% of the comoving diameter of the last scattering surface. We translate this generic constraint into limits on the parameters that characterize manifolds with each of the nine possible non-trivial orientable Euclidean topologies, and provide a code which computes these constraints. In all but the simplest cases, the shortest non-contractible loop in the space can avoid us, and be shorter than the diameter of the last scattering surface by a factor ranging from 2 to at least 6. Further, this result implies that a broader range of manifolds is observationally allowed than widely appreciated. Probing these manifolds will require more subtle statistical signatures than matched circles, such as off-diagonal correlations of harmonic coefficients.

79 ASTRONOMY AND ASTROPHYSICS↗

Constraints on Cosmological Parameters From The Sunyaev-Zel'dovich Effect and Thermal Bremsstrahlung

We discuss the possibility of constraining cosmological parameters using the Sunyaev-Zel'dovich (SZ) effect and thermal bremsstrahlung caused by intra-cluster gas in clusters of galaxies. The new generation of X-ray satellites and ground based interferometers dedicated to SZ observations will enable one to reduce uncertainties in these measurements, and thus make this method potentially quite promising in the near future. The importance of this method is that, unlike most other methods, it is based on physical principles, no 'standard candles' or 'rulers' needed. We estimate the accuracy achievable in the determination of the matter density, OMEGA(sub m), the cosmological constant, OMEGA(sub LAMBDA), and the Hubble constant, h, using the redshift dependence of the angular diameter distance derived from observations in the near future. We demonstrate that constraints from the angular diameter distance are orthogonal to those from Cosmic Microwave Background (CMB) fluctuations in the parameter space defined by OMEGA(sub m), OMEGA(sub LAMBDA), and h. Assuming a statistical error of five percent in the angular diameter distance for each cluster in a sample of five hundred clusters, we show that the redshift dependence of the angular diameter distance combined with constraints from CMB fluctuations can put stringent constraints on OMEGA(sub m) (+/- 0.03), OMEGA(sub LAMBDA) (+/- 0.03) and h (+/- 0.03, 3(sigma) errors). We also show that, with as few as 50 clusters between redshifts 0.01 and 1.5 with an assumed 10% statistical error in the angular diameter distance determination, one can distinguish between models with zero cosmological constant and spatially flat models with a cosmological constant with high confidence level (independently from the supernova results) and put meaningful constraints on OMEGA(sub m) (+/- 0.01), OMEGA(sub LAMBDA) (+/- 0.01) and h (+/- 0.01, 3(sigma) errors). With the expected advances in observational technology, we will be limited by systematic errors. We discuss the sources of systematic errors and how they can be reduced.

Molnar, S. M.↗

Constraining time dependent dark matter signals from the Sun

Dark matter (DM) particles captured by the Sun can produce high-energy electrons outside the Sun through annihilating into meta-stable mediators. The corresponding cosmic-ray electron signals observed by the space-based experiments will be time dependent due to the orbital motion of the space-based detectors. The shape of this time dependence is predictable given the orbital information of the detectors. Since the high-energy CR electron (with energy E > 100 GeV) fluxes are expected to be constant in time, non-observation of such time variation can be used to place upper limits on the DM annihilation cross section. Here we analyze the time dependence of dark matter cosmic-ray signals in three space-based experiments: AMS-02, DAMPE and CALET. Under the assumption that no time dependent signal is observed, we derive the 95% C.L. exclusion limits on the signal strength from the current data. We map our limits onto the parameter space of the dark photon model and find that the constraints are comparable with that derived from the supernova SN1987A.

79 ASTRONOMY AND ASTROPHYSICS↗

An Expert System-Driven Method for Parametric Trajectory Optimization During Conceptual Design

During the early phases of engineering design, the costs committed are high, costs incurred are low, and the design freedom is high. It is well documented that decisions made in these early design phases drive the entire design's life cycle cost. In a traditional paradigm, key design decisions are made when little is known about the design. As the design matures, design changes become more difficult in both cost and schedule to enact. The current capability-based paradigm, which has emerged because of the constrained economic environment, calls for the infusion of knowledge usually acquired during later design phases into earlier design phases, i.e. bringing knowledge acquired during preliminary and detailed design into pre-conceptual and conceptual design. An area of critical importance to launch vehicle design is the optimization of its ascent trajectory, as the optimal trajectory will be able to take full advantage of the launch vehicle's capability to deliver a maximum amount of payload into orbit. Hence, the optimal ascent trajectory plays an important role in the vehicle's affordability posture yet little of the information required to successfully optimize a trajectory is known early in the design phase. Thus, the current paradigm of optimizing ascent trajectories involves generating point solutions for every change in a vehicle's design parameters. This is often a very tedious, manual, and time-consuming task for the analysts. Moreover, the trajectory design space is highly non-linear and multi-modal due to the interaction of various constraints. When these obstacles are coupled with the Program to Optimize Simulated Trajectories (POST), an industry standard program to optimize ascent trajectories that is difficult to use, expert trajectory analysts are required to effectively optimize a vehicle's ascent trajectory. Over the course of this paper, the authors discuss a methodology developed at NASA Marshall's Advanced Concepts Office to address these issues. The methodology is two-fold: first, capture the heuristics developed by human analysts over their many years of experience; and secondly, leverage the power of modern computing to evaluate multiple trajectories simultaneously and therefore enable the exploration of the trajectory's design space early during the pre- conceptual and conceptual phases of design. This methodology is coupled with design of experiments in order to train surrogate models, which enables trajectory design space visualization and parametric optimal ascent trajectory information to be available when early design decisions are being made.

Dees, Patrick D.↗

Gravitational Wave Sources as Timing References for LISA Data

In the megahertz gravitational-wave band, galactic ultra-compact binaries (UCBs) are continuous sources emitting at near-constant frequency. The signals from many of these galactic binaries will be sufficiently strong to be detectable by the Laser Interferometer Space Antenna (LISA) after approximately Omicron (1 week) of observing. In addition to their astrophysical value, these UCBs can be used to monitor the data quality of the observatory. This paper demonstrates the capabilities of galactic UCBs to be used as calibration sources for LISA by demanding signal coherence between adjacent week-long data segments separated by a gap in time of a priori unknown duration. A parameter for the gap duration is added to the UCB waveform model and used in a Markov-chain Monte Carlo algorithm simultaneously fitting for the astrophysical source parameters. Results from measurements of several UCBs are combined to produce a joint posterior on the gap duration. The measurement accuracy's dependence on how much is known about the UCBs through prior observing, and seasonal variations due to the LISA orbital motion, is quantified. The duration of data gaps in a two-week segment of data can be constrained to within approximately 0.2 s using Omicron (10) UCBs after one month of observing. The timing accuracy from UCBs improves to approximately or less than o.1 s after 1 year of mission operations. These results are robust to within a factor of approximately 2 when taking into account seasonal variations.

Gravitation Wave Sources↗

Constraining Cosmology with Simulation-based inference and Optical Galaxy Cluster Abundance

We test the robustness of simulation-based inference (SBI) in the context of cosmological parameter estimation from galaxy cluster counts and masses in simulated optical datasets. We construct ``simulations'' using analytical models for the galaxy cluster halo mass function (HMF) and for the observed richness (number of observed member galaxies) to train and test the SBI method. We compare the SBI parameter posterior samples to those from an MCMC analysis that uses the same analytical models to construct predictions of the observed data vector. The two methods exhibit comparable performance, with reliable constraints derived for the primary cosmological parameters, ($\Omega_m$ and $\sigma_8$), and richness-mass relation parameters. We also perform out-of-domain tests with observables constructed from galaxy cluster-sized halos in the Quijote simulations. Again, the SBI and MCMC results have comparable posteriors, with similar uncertainties and biases. Unsurprisingly, upon evaluating the SBI method on thousands of simulated data vectors that span the parameter space, SBI exhibits worsened posterior calibration metrics in the out-of-domain application. We note that such calibration tests with MCMC is less computationally feasible and highlight the potential use of SBI to stress-test limitations of analytical models, such as in the use for constructing models for inference with MCMC.

79 ASTRONOMY AND ASTROPHYSICS↗

Type Ia Supernova Growth-rate Measurement with LSST Simulations: Intrinsic Scatter Systematics

Measurement of the growth rate of structures (fσ 8 ) with Type Ia supernovae (SNe Ia) will improve our understanding of the nature of dark energy and enable tests of general relativity. In this paper, we generate simulations of the 10 yr SN Ia data set of the Rubin-LSST survey, including a correlated velocity field from an N-body simulation and realistic models of SNe Ia properties and their correlations with host-galaxy properties. We find, similar to SN Ia analyses that constrain the dark energy equation-of-state parameters w 0 w a , that constraints on fσ 8 can be biased depending on the intrinsic scatter of SNe Ia. While for the majority of intrinsic scatter models we recover fσ 8 with a precision of ∼13%–14%, for the most realistic dust-based model, we find that the presence of non-Gaussianities in Hubble diagram residuals leads to a bias on fσ 8 of ∼ −20%. When trying to correct for the dust-based intrinsic scatter, we find that the propagation of the uncertainty on the model parameters does not significantly increase the error on fσ 8 . We also find that while the main component of the error budget of fσ 8 is the statistical uncertainty (>75% of the total error budget), the systematic error budget is dominated by the uncertainty on the damping parameter, σ u , that gives an empirical description of the effect of redshift space distortions on the velocity power spectrum. Our results motivate a search for new methods to correct for the non-Gaussian distribution of the Hubble diagram residuals, as well as an improved modeling of the damping parameter.

Carreres, Bastien [Duke Univ., Durham, NC (United ↗

An analytic and complete equation of state for condensed phase materials

Analytic equations of state (EOS) are intended to reproduce theoretical and experimental data in a single phase portion of the thermodynamic space. We devise a complete and thermodynamically consistent model with four distinct features: (1) a reference isotherm that remains thermodynamically stable, (2) a flexible specific heat model based on a fourth-order rational polynomial, (3) a Grüneisen parameter that depends on specific volume and temperature, and (4) pressure and internal energy functions that can be inverted analytically in temperature. The model aims to improve the accuracy of existing equations of state while remaining computationally efficient. To demonstrate its features, we include calibrations for single-crystal pentaerythritol tetranitrate (PETN), liquid nitromethane (NM), and hexagonal close-packed beryllium (Be) metal. The parameter optimization uses the specific heat capacity, Grüneisen parameter, and static compression curves obtained from density functional theory for the crystalline solids and molecular dynamics simulations for liquid NM. We also present a velocity autocorrelation function that yields accurate phonon densities of states for the EOS calibration from the molecular dynamics trajectories. Each of the three calibrations is constrained to enforce the ambient state from experimental measurements and validated against experimental Hugoniot data from multiple sources. We also include one-dimensional hydrodynamic simulations of the isentropic compression experiments for beryllium conducted at the Z facility.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗