SEARCH · Engineering Papers

Results for “density”

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 343 records · Page 19

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Method for printing low-density polymer structures

A method of forming a low-density three-dimensional article is provided. The method includes printing a low-density composition on a substrate to form at least one layer comprising the low-density composition. The low-density composition includes a (P) polymer component and (M) a microsphere component in a ratio by volume (P):(M). The method also includes selectively controlling a density of the low-density composition during printing to give the at least one layer on the substrate. Selectively controlling the density of the low-density composition includes varying the ratio (P):(M) during printing. The method further includes repeating the printing and selectively controlling the density of the low-density composition to form additional layer(s), thereby forming the low-density three-dimensional article. A low-density three-dimensional article prepared in accordance with the method is also provided.

Kunc, Vlastimil↗

DOE OSTI2026

Propagation method and planting density influence canopy developmental transition and biomass productivity in Miscanthus × giganteus

Understanding how establishment practices influence the mechanisms underlying Miscanthus × giganteus (miscanthus) productivity and canopy development is critical for optimizing management. Data was collected during the juvenile (2011–2013) and mature (2024) phases of a long-term field experiment established in Urbana, Illinois, to evaluate the effects of propagation method (plug propagation [PP] and rhizome propagation [RP]), planting density (1.0, 0.75, and 0.25 plants m⁻²), and nitrogen application (0 and 67 kg N ha⁻¹) on end-of-season biomass yield, tiller mass, tiller density, and tiller height. Linear regression models identified the dominant predictors of yield across stand ages and management regimes. Planting density, nitrogen (N) application, and propagation method significantly influenced early yield and canopy development. During the juvenile phase, biomass yield was driven by tiller density due to canopy expansion; in the mature phase, yield became driven by tiller mass. The PP plots produced higher tiller density than the RP plots, resulting in faster canopy closure and higher juvenile-phase yields. Rhizome-propagated (RP) plots produced lower tiller density, but individual tillers were 3.3–6.4 g tiller −1 heavier than PP tillers. After the canopy reached equilibrium, the PP and RP yields were similar because greater RP tiller mass compensated for its lower tiller density. Higher planting density resulted in greater yield and tiller density during the second year (2012), but this effect was absent from the third year (2013) onward. In the juvenile phase, N fertilization enhanced yield by 1.6–3.4 Mg ha −1 . Initiating fertilization in 2013 on unfertilized plots produced biomass similar to that in fertilized plots, suggesting yield recovery in the mature phase. These findings revealed that establishment strategies, including propagation method and planting density, influence juvenile miscanthus canopy development and productivity, transitioning from tiller-density- to mass-dominated yields, but not mature phase productivity.

09 BIOMASS FUELS↗

DOE OSTI2022

Impact of the electron density and temperature gradient on drift-wave turbulence in the Large Plasma Device

In this paper we present an experimental study of edge turbulence in the Large Plasma Device at UCLA. We utilize a scan of discharge power and prefill pressure (neutral density) to show experimentally that turbulent density fluctuations decrease with decreasing density gradient, as predicted for resistive drift-wave turbulence (RDWT). As expected for RDWT, we observe that the cross-phase between the density and potential fluctuations is close to 0. Moreover, the addition of an electron temperature gradient leads to a reduction in the amplitude of the density fluctuations, as expected for RDWT. However, counter to theoretical expectations, we find that the potential fluctuations do not follow the same trends as the density fluctuations for changes either in density gradients or the addition of a temperature gradient. The disconnect between the density and potential fluctuations is connected to changes in the parallel flows as a result of differences in the prefill pressure, i.e. neutral density. Further analysis of the density and potential fluctuation spectra show that the electron temperature gradient reduces the low frequency fluctuations up to $10 \,{\rm kHz}$ and the introduction of a temperature gradient leads to an unexpected ${\sim }{\rm \pi}$ shift of the density–potential cross-phase at ${\sim }10\,{\rm kHz}$ , while maintaining the typical resistive drift-wave cross-phase at lower frequencies. These experiments partly confirm existing knowledge on resistive drift-wave turbulence, but also introduce new observations that indicate a need for dedicated nonlinear three-dimensional turbulence simulations that include neutrals.

Physics↗

DOE OSTI2025

Habitat-specific seasonal densities of sympatric raccoons and opossums in the southeastern United States

Raccoon and opossum densities have implications for rabies management, but estimates of seasonal densities of both species are lacking for rural nonagricultural habitats of the southeastern United States, a core portion of their geographic range. Consequently, it remains unclear whether the densities of 1 species limit the other, which is possible considering their substantial niche overlap. We carried out a mark–recapture study of raccoons and opossums in 4 rural nonagricultural habitats (bottomland hardwood forest, riparian forest, upland pine forest, and isolated wetlands) in South Carolina, United States (2020 to 2022), and combined this with previous data from the same habitats (2017 to 2019) to estimate habitat-specific spring and fall densities. Raccoon densities ranged from 5.17 ± 0.96 animals/km 2 (bottomland hardwood fall) to 1.63 ± 0.83 animals/km 2 (upland pine spring) and were on average 19% higher in fall compared to spring. Opossum densities ranged from 10.35 ± 1.98 animals/km 2 (bottomland hardwood fall) to 1.11 ± 1.55 animals/km 2 (upland pine spring) with divergent seasonal patterns among habitats. These low densities across all habitats compared to other studies are likely the result of low resource availability, consistent with other habitats that have minimal anthropogenic influence. We observed a positive association between raccoon and opossum densities across trapping grids, suggesting that raccoons do not suppress opossum densities, but that densities of both species increase with increasing resource availability. Furthermore, our results can be used to inform oral rabies vaccination efforts such as refining bait densities and timing of bait distribution in these habitats.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2021

Nuclear charge densities in spherical and deformed nuclei: Toward precise calculations of charge radii

Precise measurements of atomic transitions affected by electron-nucleus hyperfine interactions offer sensitivity to explore basic properties of the atomic nucleus and study fundamental symmetries, including the search for new physics beyond the standard model of particle physics. Specifically, such measurements, augmented by atomic and nuclear calculations, will enable extraction of the higher-order radial moments of the charge-density distribution in spherical and deformed nuclei. The new data impose higher precision requirements on a theoretical description. The nuclear charge density is composed of the proton point distribution folded with the nucleonic charge distributions. The latter induce subtle relativistic corrections due to the coupling of nucleon magnetic moments with the nuclear spin-orbit density. Additional corrections come from the effect of center-of-mass projection. We assess the precision of nuclear charge density calculations by studying the behavior of relativistic and center-of-mass motion corrections to the second and fourth charge radial moments. Special attention has been paid to the magnetic spin-orbit density associated with the local variations of the spin-orbit current. The calculations for semimagic and open-shell nuclei are performed in the framework of self-consistent mean-field theory using quantified energy density functionals and density-dependent pairing forces. We used the general expression for the spin-orbit form factor that is valid for spherical and deformed nuclei. We studied the impact of various correction terms on the charge radii, fourth radial moments, diffraction radii, and surface thickness of spherical and deformed nuclei. The spin-orbit corrections to charge radial moments and surface thickness show strong shell fluctuations which can make an appreciable effect when aiming at high-precision predictions of isotopic shifts. The inclusion of relativistic and center-of-mass corrections impacts the quality of energy density functionals optimized to charge radii data. To establish reliable constraints on the existence of new forces from isotope shift measurements, precise calculations of nuclear charge densities of deformed nuclei are needed. The proper inclusion of the spin-orbit charge density and other correction terms is essential when aiming at extraction of subtle effects which become particularly visible in isotopic trends. It is also important when developing high-quality nuclear energy density functionals optimized using heterogeneous datasets involving absolute charge radii, differential charge radii, and charge form factor properties deduced from electron-scattering data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS2010

Examining Periodic Solar-Wind Density Structures Observed in the SECCHI Heliospheric Imagers

We present an analysis of small-scale, periodic, solar-wind density enhancements (length scales as small as approximately equals 1000 Mm) observed in images from the Heliospheric Imager (HI) aboard STEREO-A. We discuss their possible relationship to periodic fluctuations of the proton density that have been identified at 1 AU using in-situ plasma measurements. Specifically, Viall, Kepko, and Spence examined 11 years of in-situ solar-wind density measurements at 1 AU and demonstrated that not only turbulent structures, but also nonturbulent, periodic density structures exist in the solar wind with scale sizes of hundreds to one thousand Mm. In a subsequent paper, Viall, Spence, and Kasper analyzed the alpha-to-proton solar-wind abundance ratio measured during one such event of periodic density structures, demonstrating that the plasma behavior was highly suggestive that either temporally or spatially varying coronal source plasma created those density structures. Large periodic density structures observed at 1 AU, which were generated in the corona, can be observable in coronal and heliospheric white-light images if they possess sufficiently high density contrast. Indeed, we identify such periodic density structures as they enter the HI field of view and follow them as they advect with the solar wind through the images. The smaller, periodic density structures that we identify in the images are comparable in size to the larger structures analyzed in-situ at 1 AU, yielding further evidence that periodic density enhancements are a consequence of coronal activity as the solar wind is formed.

Viall, Nicholeen M.↗

NASA NTRS2017

Meteoroid Bulk Density and Ceplecha Types

Determination of asteroid bulk density is an important aspect of NEO characterization, yet difficult to measure. As a fraction of meteoroids originate from asteroids (including some NEOs), a study of meteoroid bulk densities can potentially provide useful insights into the densities of NEOs and PHOs in lieu of mutual perturbations, satellite, or expensive spacecraft missions. NASA's Meteoroid Environment Office characterizes the meteoroid environment for the purpose of spacecraft risk and operations. To accurately determine the risk, a distribution of meteoroid bulk densities are needed. This is not trivial to determine. If the particle survives to the ground the bulk density can be directly measured, however only the most dense particles land on the Earth. The next best approach is to model the meteor's ablation, which is not straightforward. Clear deceleration is necessary to do this and there are discrepancies in results between models. One approach to a distribution of bulk density is to use a measured proxy for the densities, then calibrate the proxy with known densities from meteorite falls, ablation modelling, and other sources. An obvious proxy choice is the Ceplecha type, K(sub B), thought to indicate the strength of a meteoroid. KB is frequented cited as a good proxy for meteoroid densities, but we find it is poorly correlated with density. However, a distinct split by dynamical type was seen with Jovian Tisserand parameter, T(sub J), with meteoroids from Halley Type comets (T(sub J less than 2 ) exhibiting much lower densities than those originating from Jupiter and asteroids (T(sub J greater than 2).

Blaauw, R. C.↗

DOE OSTI2018

Understanding Clinical Mammographic Breast Density Assessment: a Deep Learning Perspective

Mammographic breast density has been established as an independent risk marker for developing breast cancer. Breast density assessment is a routine clinical need in breast cancer screening and current standard is using the Breast Imaging and Reporting Data System (BI-RADS) criteria including four qualitative categories (i.e., fatty, scattered density, heterogeneously dense, or extremely dense). In each mammogram examination, a breast is typically imaged with two different views, i.e., the mediolateral oblique (MLO) view and cranial caudal (CC) view. The BI-RADS-based breast density assessment is a qualitative process made by visual observation of both the MLO and CC views by radiologists, where there is a notable inter- and intra-reader variability. In order to maintain consistency and accuracy in BI-RADS-based breast density assessment, gaining understanding on radiologists’ reading behaviors will be educational. In this study, we proposed to leverage the newly emerged deep learning approach to investigate how the MLO and CC view images of a mammogram examination may have been clinically used by radiologists in coming up with a BI-RADS density category. We implemented a convolutional neural network (CNN)-based deep learning model, aimed at distinguishing the breast density categories using a large (15,415 images) set of real-world clinical mammogram images. Our results showed that the classification of density categories (in terms of area under the receiver operating characteristic curve) using MLO view images is significantly higher than that using the CC view. This indicates that most likely it is the MLO view that the radiologists have predominately used to determine the breast density BI-RADS categories. Our study holds a potential to further interpret radiologists’ reading characteristics, enhance personalized clinical training to radiologists, and ultimately reduce reader variations in breast density assessment.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2025

A Semi-Empirical Density Law for Ternary, Homogeneous PuCl 3 /HCl/H 2 O Solutions

Nuclear material operations pose unique hazards that are not encountered in other chemical, energy, or manufacturing industries. One of these hazards is the potential for a nuclear criticality accident when handling fissile isotopes such as 235 U and 239 Pu. These hazards are particularly high when fissile material is dissolved in solution as the neutron behaviors of the system can change rapidly with the physical and chemical changes accessible in solution. Current estimates of solution density used for criticality safety are outdated and hinder fissile material handling. Developing new estimates for these safety calculations requires experimental characterization and the derivation of empirical density models. We have derived a density law describing PuCl 3 /HCl/H 2 O solutions from experimental data characterizing solution density. Density data was treated using a Pitzer-derived eight-parameter equation, defining density as a function of analyte concentrations, temperature, and interactions between these variables. The model is predictive across the concentration and temperature ranges from which it was derived. The potential effects of varying oxidation states of plutonium, which are easily accessible in aqueous media, on the bulk solution density of the ternary system were also investigated. The resulting Pitzer-derived density law was applied to a nuclear criticality safety model, and the impact of the experimental characterization of solution density relative to previous estimates was demonstrated to be significant and suggest that the current approach to estimating density in nuclear criticality safety calculations may lead to overly conservative controls.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Quantitative radiography for determining density fluctuations in HED experiments

We have developed a method to extract density fluctuation measurements from x-ray radiographs of high-energy density (HED) instability growth and turbulence experiments. We use this information to calculate density fluctuation statistics for constraining the performance of turbulent mix models in HED systems. The density calculation combines image filtering, removal of systemic effects such as backlighter variation, calculation of transmission across multiple materials, and use of tracer materials to generate an approximate single-material density field. From the density map, we calculate both average density and a variance-like moment b (density-specific-volume covariance), which we compare to our models. We infer both quantities from a single image, which is significantly more information than the historic single scalar mix width measurements. We also develop a method of analyzing simulation outputs that incorporate both the density fluctuation metric from a turbulence model and the bulk material maps from the hydrodynamic code. This analysis helps address the question of how to initialize the simulations for best comparison to data from systems with large separations of scale in the mixing perturbation initial condition. We find that our data analysis method yields 1D average density and b curves with similar morphology and amplitudes as those from preliminary simulation comparisons.

47 OTHER INSTRUMENTATION↗

DOE OSTI2017

Quantitative Analysis of Ultralow-Density Materials Using Laboratory-Based Quasi-Monochromatic Radiography

The measurement of the density of materials, especially ultralow-density foams, is difficult in that the measurement must be precise and localizable. The density of the material is often governed by its cellular (i.e., porous) structure, and many techniques exist to create that structure. Often, the cellular structure can vary from one location within the material to another, and when at low densities (i.e., densities lower than ~500 mg/cm 3 ), it can vary due to shrinkage during syneresis, collapse under the weight of gravity, or gas/water vapor uptake. Quantifying this variation is important for a variety of applications, especially when used in plasma physics targets. Knowing the density and its variation across the sample is critical for experimental results to be accurately predicted by physics calculations and for modeling the results of the physics targets. The use of quasi-monochromatic radiography provides a means to image the two-dimensional (2-D) distribution of density variation within silica aerogel materials and to quantitatively measure that variation from sample to sample and lot to lot. For this work, two batches of silica aerogels with targeted densities of ~20 mg/cm 3 were created, one batch at Lawrence Livermore National Laboratory, and the other batch at Los Alamos National Laboratory. Outlined here is a quasi-monochromatic radiography system using various X-ray sources coupled to a doubly curved crystal optic and X-ray charge-coupled device camera to image and characterize these materials. It was found that measuring the density both gravimetrically and using quasi-monochromatic radiography were statistically identical, although the two batches were found to be slightly higher than their targeted density due to shrinkage. The radiography system also provided 2-D information as to the aerogel quality, i.e., presence of voids, chipped material, or inclusions.

36 MATERIALS SCIENCE↗

DOE OSTI2018

Convergence of Probability Densities Using Approximate Models for Forward and Inverse Problems in Uncertainty Quantification

Here, we analyze the convergence of probability density functions utilizing approximate models for both forward and inverse problems. We consider the standard forward uncertainty quantification problem where an assumed probability density on parameters is propagated through the approximate model to produce a probability density, often called a push-forward probability density, on a set of quantities of interest (QoI). The inverse problem considered in this paper seeks to update an initial probability density assumed on model input parameters such that the subsequent push-forward of this updated density through the parameter-to-QoI map matches a given probability density on the QoI. We prove that the densities obtained from solving the forward and inverse problems, using approximate models, converge to the true densities as the approximate models converge to the true models. Numerical results are presented to demonstrate convergence rates of densities for sparse grid approximations of parameter-to-QoI maps and standard spatial and temporal discretizations of PDEs and ODEs.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2013

Evidence of Temporal Variation of Titan Atmospheric Density in 2005-2013

One major science objective of the Cassini mission is an investigation of Titan's atmosphere constituent abundances. Titan's atmospheric density is of interest not only to planetary scientists but also to mission design and mission control engineers. Knowledge of the dependency of Titan's atmospheric density with altitude is important because any unexpectedly high atmospheric density has the potential to tumble the spacecraft during a flyby. During low-altitude Titan flyby, thrusters are fired to counter the torque imparted on the spacecraft due to the Titan atmosphere. The denser the Titan's atmosphere is, the higher are the duty cycles of the thruster firings. Therefore thruster firing telemetry data could be used to estimate the atmospheric torque imparted on the spacecraft. Since the atmospheric torque imparted on the spacecraft is related to the Titan's atmospheric density, atmospheric densities are estimated accordingly. In 2005-2013, forty-three low-altitude Titan flybys were executed. The closest approach altitudes of these Titan flybys ranged from 878 to 1,074.8 km. Our density results are also compared with those reported by other investigation teams: Voyager-1 (in November 1980) and the Huygens Atmospheric Structure Instrument, HASI (in January 2005). From our results, we observe a temporal variation of the Titan atmospheric density in 2005-2013. The observed temporal variation is significant and it isn't due to the estimation uncertainty (5.8%, 1 sigma) of the density estimation methodology. Factors that contributed to this temporal variation have been conjectured but are largely unknown. The observed temporal variation will require synergetic analysis with measurements made by other Cassini science instruments and future years of laboratory and modeling efforts to solve. The estimated atmospheric density results are given in this paper help scientists to better understand and model the density structure of the Titan atmosphere.

Atmosphere↗

DOE OSTI2023

Summary Analysis of Different Offshore Wind Capacity Density Drivers in Proposed U.S. Projects and Impacts on Progress Towards State and Federal Deployment Targets

Understanding the density of offshore wind development - or the "capacity density" - is important for state and federal offshore wind planning efforts and to track progress towards policy goals. NREL has traditionally assumed a capacity density of 3 megawatts per square kilometer (MW/km2) when estimating the U.S. project pipeline (Musial et al. 2022), but capacity densities for existing European offshore wind farms ranged from 2-19 MW/km2 (Borrmann et al. 2018). We track the expected capacity densities of the proposed U.S. projects in the East Coast and make an analysis by state and by developer to identify what are the capacity density norms in the industry and whether the state in which the development efforts are taking place has some impact on capacity density. We also analyze the effects of the area occupied by station keeping systems on capacity density for lease areas dedicated to floating projects and provide a list of main capacity density drivers. This summary analysis sets a baseline to study how different capacity density drivers could impact the progress towards state and federal deployment targets. Given the currently delineated Renewable Energy Areas in the Outer Continental Shelf (OCS) and the basis from the summary analysis, we estimate the expected offshore wind deployment over time federally and by state and compare it to the state and federal deployment targets. The proposed analysis helps identify gaps and needs to meet state and federal targets, which provides sufficient information to propose strategic supply chain and policy planning with the objective of meeting deployment goals.

capacity density↗

DOE OSTI2019

A semi-Lagrangian direct-interaction closure of the spectra of isotropic variable-density turbulence

A study of variable-density homogeneous stationary isotropic turbulence based on the sparse direct-interaction perturbation (SDIP) and supporting direct numerical simulations (DNS) is presented. The non-solenoidal flow considered here is an example of turbulent mixing of gases with different densities. The spectral statistics of this type of flow are substantially more difficult to understand theoretically than those of the similar solenoidal flows. In the approach described here, the nonlinearly coupled velocity and scalar (which determine the density of the fluid) equations are expanded in terms of a normalised density ratio parameter. A new set of coupled integro-differential SDIP equations are derived and then solved numerically for the first-order correction to the incompressible equations in the variable-density expansion parameter. By adopting a regular expansion approach, one obtains leading-order corrections that are universal and therefore interesting in their own right. The predictions are then compared with DNS of forced variable-density flow with different density contrasts. It is found that the velocity spectrum owing to variable density is indistinguishable from that of constant-density turbulence, as it is supported by a wealth of indirect experimental evidence, but the scalar spectra show significant deviations, and even loss of monotonicity, as a function of the type and strength of the large-scale source of the mixing. Furthermore, the analysis helps clarify what may be the proper approach to interpret the power spectrum of variable-density turbulence.

Mechanics↗

DOE OSTI2023

Structure of iso-density sets in supersonic isothermal turbulence

Context: The gas density structure of the cold molecular phase of the interstellar medium is the main controller of star formation. Aims: A theoretical framework is proposed to describe the structural content of the density field in isothermal supersonic turbulence. Methods: It makes use of correlation and structure functions of the phase indicator field defined for different iso-density values. The relations between these two-point statistics and the geometrical features of iso-density sets such as the volume fraction, the surface density, the curvature, and fractal characteristics are provided. An exact scale-by-scale budget equation is further derived revealing the role of the turbulent cascade and dilation on the structural evolution of the density field. Although applicable to many flow situations, this tool is here first invoked for characterising supersonic isothermal turbulence, using data from the currently best-resolved numerical simulation. Results: We show that iso-density sets are surface fractals rather than mass fractals, with dimensions that markedly differ between dilute, neutral, and dense regions. The surface–size relation is established for different iso-density values. We further find that the turbulent cascade of iso-density sets is directed from large towards smaller scales, in agreement with the classical picture that turbulence acts to concentrate more surface into smaller volumes. Intriguingly, there is no range of scales that complies with a constant transfer rate in the cascade, challenging our fundamental understanding of interstellar turbulence. Finally, we recast the virial theorem in a new formulation drawing an explicit relation between the aforementioned geometrical measures and the dynamics of iso-density sets.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

The effect of the isotope on the H-mode density limit

In order to understand the mechanisms for the H-mode density limit in machines with fully metallic walls, systematic investigations of H-mode density limit plasmas in experiments with deuterium and hydrogen external gas fuelling have been performed on JET-ILW. The observed H-mode density limit on JET in D- as well as in H-plasmas demonstrates similar operation phases: the stable H-mode phase, degrading H-mode, breakdown of the H-mode with energy confinement deterioration accompanied by a dithering cycling phase, followed by the L-mode phase. The density limit is not related to an inward collapse of the hot core plasma due to an overcooling of the plasma periphery by radiation. Indeed, independently of the isotopic effect, the total radiated power stay almost constant during the H-mode phase until the H-L back transition. It was observed in D- and H-plasmas that neither detachment, nor the X-point MARFE itself do trigger the H-L transition and that they thus do not present a limit on the plasma density. It is the plasma confinement, most likely determined by edge parameters, which is ultimately responsible for the H-mode DL. By comparing similar discharges but fuelled with either deuterium or hydrogen, we have found that the H-mode density limit exhibits a dependence on the isotope mass: the density limit is up to 35% lower in hydrogen compared to similar deuterium plasma conditions (the obtained density limit is in agreement with the Greenwald limit for D-plasma). In addition, the density limit is nearly independent of the applied power both in deuterium or hydrogen fuelling conditions. The measured Greenwald fractions are consistent with the predictions from a theoretical model based on an MHD instability theory in the near-SOL. The JET operational domains are significantly broadened when increasing the plasma effective mass (e.g. tritium or deuterium–tritium operation), i.e. the L to H power threshold is reduced whereas the density limit for the L-mode back transition is increased.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Direct astrophysical tests of chiral effective field theory at supranuclear densities

Recent observations of neutron stars with gravitational waves and X-ray timing provide unprecedented access to the equation of state (EoS) of cold dense matter at densities difficult to realize in terrestrial experiments. At the same time, predictions for the EoS equipped with reliable uncertainty estimates from chiral effective field theory ( x EFT) allow us to bound our theoretical ignorance. In this work, we analyze astrophysical data using a nonparametric representation of the neutron-star EoS conditioned on x EFT to directly constrain the underlying physical properties of the compact objects without introducing modeling systematics. We discuss how the data alone constrain the EoS at high densities when we condition on x EFT at low densities. Here, we also demonstrate how to exploit astrophysical data to directly test the predictions of x EFT for the EoS up to twice nuclear saturation density, in order to estimate the density at which these predictions might break down. We nd that the existence of massive pulsars, gravitational waves from GW170817, and NICER observations of PSR J0030+0451 favor x EFT predictions for the EoS up to nuclear saturation density over a more agnostic analysis by as much as a factor of 7 for the quantum Monte Carlo (QMC) calculations used in this work. While x EFT predictions using QMC are fully consistent with gravitational-wave data up to twice nuclear saturation density, NICER observations suggest that the EoS stiffens relative to these predictions at or slightly above nuclear saturation density. Additionally, for these QMC calculations, we marginalize over the uncertainty in the density at which x EFT begins to break down, constraining the radius of a 1.4M $\bigodot$ neutron star to R 1.4 = $11.40^{+1.38}_{–1.04}$ ($12.54^{+0.71}_{–0.63}$) km and the pressure at twice nuclear saturation density to p(2n sat ) = $14.2^{+18.1}_{–8.4}$ ($28.7^{+15.3}_{–15.0}$) MeV=fm 3 with massive pulsar and gravitational-wave (and NICER) data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗