Search NASA⌕ Search

SEARCH · Search NASA

Results for “radius equation”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Constraining the Dense Matter Equation of State with Joint Analysis of NICER and LIGO/Virgo Measurements

The Neutron Star Interior Composition Explorer collaboration recently published a joint estimate of the mass and the radius of PSR J0030+0451, derived via X-ray pulse-profile modeling. Raaijmakers et al. explored the implications of this measurement for the dense matter equation of state (EOS) using two parameterizations of the high-density EOS: a piecewise-polytropic model, and a model based on the speed of sound in neutron stars (NSs). In this work we obtain further constraints on the EOS following this approach, but we also include information about the tidal deformability of NSs from the gravitational wave signal of the compact binary merger GW170817. Further, we compare the constraints on the EOS to those set by the recent measurement of a 2.14 Me pulsar, included as a likelihood function approximated by a Gaussian, and find a small increase in information gain. To show the flexibility of our method, we also explore the possibility that GW170817 was a NS–black hole merger, which yields weaker constraints on the EOS.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Confronting new NICER mass-radius measurements with phase transition in dense matter and twin compact stars

The (re)analysis of data on the X-ray emitting pulsars PSR J0030+0451 and J0740+6620, as well as new results on PSR J0437-4715 and J1231-1411, are confronted with the predictions of the equation of state (EoS) models allowing for strong first-order phase transition for the mass-radius (M-R) diagram. Here, we use models that are based on a covariant density functional (CDF) EoS for nucleonic matter at low densities and a quark matter EoS, parameterized by the speed of sound, at higher densities. To account for the variations in the ellipses for PSR J0030+0451 obtained from different analyses, we examined three scenarios to assess their consistency with our models, focusing particularly on the potential formation of twin stars. We found that in two scenarios, where the ellipses for PSR J0030+0451 and J0437-4715 with masses close to the canonical mass ∼ 1.4 M ⊙ are significantly separated, our models allow for the presence of twin stars as a natural explanation for potential differences in the radii of these stars.

X-ray binaries↗
DOE OSTI2023

The attractor flow for AdS 5 black holes in $\mathcal{N}$ = 2 gauged supergravity

We study the flow equations for BPS black holes in $\mathcal{N}$ = 2 five-dimensional gauged supergravity coupled to any number of vector multiplets via FI couplings. We develop the Noether-Wald procedure in this context and exhibit the conserved charges as explicit integrals of motion, in the sense that they can be computed at any radius on the rotating spacetime. The boundary conditions needed to solve the first order differential equations are discussed in great detail. We extremize the entropy function that controls the near horizon geometry and give explicit formulae for all geometric variables at their supersymmetric extrema. We have also considered a complexification of the near-horizon variables that elucidates some features of the theory from the near-horizon perspective.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Neural simulation-based inference of the neutron star equation of state directly from telescope spectra

Neutron stars provide a unique opportunity to study strongly interacting matter under extreme density conditions. The intricacies of matter inside neutron stars and their equation of state are not directly visible, but determine bulk properties, such as mass and radius, which affect the star's thermal X-ray emissions. However, the telescope spectra of these emissions are also affected by the stellar distance, hydrogen column, and effective surface temperature, which are not always well-constrained. Uncertainties on these nuisance parameters must be accounted for when making a robust estimation of the equation of state. In this study, we develop a novel methodology that, for the first time, can infer the full posterior distribution of both the equation of state and nuisance parameters directly from telescope observations. This method relies on the use of neural likelihood estimation, in which normalizing flows use samples of simulated telescope data to learn the likelihood of the neutron star spectra as a function of these parameters, coupled with Hamiltonian Monte Carlo methods to efficiently sample from the corresponding posterior distribution. Our approach surpasses the accuracy of previous methods, improves the interpretability of the results by providing access to the full posterior distribution, and naturally scales to a growing number of neutron star observations expected in the coming years.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Generalized kinetic equation for tokamak plasma equilibrium distribution function

A generalised kinetic equation for the equilibrium distribution function in a finite beta, arbitrary tokamak plasma is derived. The equation is correct to second order in ρ/L (ρ is the particle Larmor radius and L is the system size). Resolving finite Larmor radius length scales with no restriction on the ratio of poloidal to total equilibrium magnetic field, B $\vartheta$ /B, it generalises the drift kinetic theory of [Hazeltine Phys. Plasmas 15 (1973) 77] to the limit of B $\vartheta$ /B ~ 1 (e.g. to ensure validity for spherical tokamaks). Two cases are considered. The first provides the equilibrium distribution function, consistent with the generalised gyrokinetic formalism of [Dudkovskaia et al. Plasma Phys. Control. Fusion 65 (2023) 045010], derived specifically to capture neoclassical equilibrium currents in gyrokinetic stability analyses in strong gradient regions. The second assumes short length scales in the direction perpendicular to the magnetic field, which can occur as a result of small coherent magnetic structures in the plasma, such as neoclassical tearing mode magnetic islands close to threshold. This then extends the drift island equations of [Dudkovskaia et al. Nucl. Fusion 63 (2023) 016020] for the plasma response to magnetic islands to a spherical tokamak plasma configuration. In conclusion, resolving ρ ~ ρ $\vartheta$ (or B $\vartheta$ ~ B), where ρ $\vartheta$ is the particle poloidal Larmor radius, is also expected to influence calculations of the magnetic island propagation frequency and the associated contributions to the island onset conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

PSR J0030+0451 Mass and Radius from NICER Data and Implications for the Properties of Neutron Star Matter

Neutron stars are not only of astrophysical interest, but are also of great interest to nuclear physicists because their attributes can be used to determine the properties of the dense matter in their cores. One of the most informative approaches for determining the equation of state (EoS) of this dense matter is to measure both a star's equatorial circumferential radius R e and its gravitational mass M. Here we report estimates of the mass and radius of the isolated 205.53 Hz millisecond pulsar PSR J0030+0451 obtained using a Bayesian inference approach to analyze its energy-dependent thermal X-ray waveform, which was observed using the Neutron Star Interior Composition Explorer (NICER). This approach is thought to be less subject to systematic errors than other approaches for estimating neutron star radii. We explored a variety of emission patterns on the stellar surface. In this work, our best-fit model has three oval, uniform-temperature emitting spots and provides an excellent description of the pulse waveform observed using NICER. The radius and mass estimates given by this model are ${R}_{e}={13.02}_{-1.06}^{+1.24}$ km and $M={1.44}_{-0.14}^{+0.15}\,{M}_{\odot }$ (68%). The independent analysis reported in the companion paper by Riley et al. explores different emitting spot models, but finds spot shapes and locations and estimates of R e and M that are consistent with those found in this work. We show that our measurements of R e and M for PSR J0030+0451 improve the astrophysical constraints on the EoS of cold, catalyzed matter above nuclear saturation density.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The antesonic condition for the explosion of core-collapse supernovae – II. Rotation and turbulence

In the problem of steady free fall on to a standing shockwave around a central mass, the ‘antesonic’ condition limits the regime of stable accretion to $c_T^2/v_\mathrm{esc}^2\le 3/16$, where cT is the isothermal sound speed in the subsonic post-shock flow, and vesc is the escape velocity at the shock radius. Above this limit, it is impossible to satisfy both the Euler equation and the shock jump conditions, and the system transitions to a wind. This physics explains the existence of a critical neutrino luminosity in steady-state models of accretion in the context of core-collapse supernovae. Here, we extend the antesonic condition to flows with rotation and turbulence using a simple one-dimensional formalism. Both effects decrease the critical post-shock sound speed required for explosion. While quite rapid rotation is required for a significant change to the critical condition, we show that the level of turbulence typically achieved in supernova simulations can greatly impact the critical value of $c_T^2/v_\mathrm{esc}^2$. A core angular velocity corresponding to a millisecond rotation period after contraction of the proto-neutron star results in only a ~5 percent reduction of the critical curve. In contrast, near-sonic turbulence with specific turbulent kinetic energy $K/c_T^2=0.5-1$, leads to a decrease in the critical value of $c_T^2/v_{\rm esc}^2$ by ~20 to 40 percent. This analysis provides a framework for understanding the role of post-shock turbulence in instigating explosions in models that would otherwise fail and helps explain why multidimensional simulations explode more easily than their one-dimensional counterparts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Bayesian Inference of Hybrid Star Properties from Future High-precision Measurements of Their Radii

Future high-precision X-ray and gravitational-wave observations of neutron stars (NSs) are expected to constrain NS radii with uncertainties as small as σ ≃ 0.1 km. Such unprecedented precision offers a unique opportunity to extract new information about the nature and equation of state (EOS) of supradense matter in NS cores. Using mock radius data with uncertainties ranging from σ = 1.0 to 0.1 km, together with a flexible meta-model NS EOS that allows for a first-order hadron–quark phase transition, we perform a Bayesian statistical analysis to assess the impact of radius measurements on EOS constraints. We find that high-precision radius measurements, particularly for massive NSs, significantly tighten constraints on the hadron–quark transition density ρ t , the quark matter mass fraction in NS cores, and several parameters characterizing the EOS of supranuclear hadronic matter, although the degree of improvement depends on the assumed prior range of ρ t . In contrast, even with the highest precision considered, NS radii—including those of massive stars—remain largely insensitive to the stiffness of quark matter, independent of the measurement accuracy or the prior range adopted for ρ t .

Neutron stars↗
DOE OSTI2021

The Radius of PSR J0740+6620 from NICER and XMM-Newton Data

PSR J0740+6620 has a gravitational mass of 2.08 ± 0.07 M ⊙ , which is the highest reliably determined mass of any neutron star. As a result, a measurement of its radius will provide unique insight into the properties of neutron star core matter at high densities. Here we report a radius measurement based on fits of rotating hot spot patterns to Neutron Star Interior Composition Explorer (NICER) and X-ray Multi-Mirror (XMM-Newton) X-ray observations. We find that the equatorial circumferential radius of PSR J0740+6620 is 13.7$^{+2.6}_{–1.5}$ km (68%). We apply our measurement, combined with the previous NICER mass and radius measurement of PSR J0030+0451, the masses of two other ~2 M ⊙ pulsars, and the tidal deformability constraints from two gravitational wave events, to three different frameworks for equation-of-state modeling, and find consistent results at ~1.5-5 times nuclear saturation density. For a given framework, when all measurements are included, the radius of a 1.4 M ⊙ neutron star is known to ±4% (68% credibility) and the radius of a 2.08 M ⊙ neutron star is known to ±5%. The full radius range that spans the ±1σ credible intervals of all the radius estimates in the three frameworks is 12.45 ± 0.65 km for a 1.4 M ⊙ neutron star and 12.35 ± 0.75 km for a 2.08 M ⊙ neutron star.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Beam Envelope Variation due to Space-Charge Neutralization

Our non-invasive, time-resolved diamagnetic loop measurements of beam size during the electron current-pulse flattop at the exit of our 20-MeV linear induction accelerator revealed that the beam size varied by about 13% during the current flattop. This was an unexpected result, since both current and energy were constant during the interval of radius variation. One possible cause, poor vacuum, was investigated using a time resolved envelope equation. It was found that sufficiently high residual-gas pressure in the downstream transport region could result in the observed variation.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Verification and validation of detonation-shock-dynamics relations for explosives described by general equation of state and chemical reaction models

Detonation shock dynamics is a powerful method to model the behaviour of High Explosives (HE). However in order to use this method, the underlying relationship between the local radius of curvature and the detonation speed must be known. Previous work has developed methods to calculate this effect using simple, single-step Arrhenius and polytropic gas, models for the chemical reaction and the equation of state, respectively. In recent years, more complex models for both reaction rates and equations of state have been developed which show better agreement with experimental data than these simple models, especially when considering condensed phase explosives.. This work presents the governing equations for solving these problems in a way that is generalised to use arbitrary equations of state as well as reaction models which may have more than a single step and multiple product species. This implementation is verified against exact solutions, demonstrating that the equations were implemented properly. The verified algorithm is then validated against experimental data and high fidelity simulations, showing that it is able to make accurate predictions in a regime where the underlying assumptions of the governing equations are valid. Importantly, this approach has many applications: from creating equivalent detonation shock dynamics models for existing reactive burn calibrations for HE; to developing new functional forms and calibrations of reactive burn models for condensed phase high explosives.

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

Relativistic hybrid stars in light of the NICER PSR J0740 + 6620 radius measurement

We explore the implications of the recent radius determination of PSR J0740 + 6620 by the NICER experiment, combined with the neutron skin measurement by the PREX-II experiment and the associated inference of the slope of symmetry energy, for the structure of hybrid stars with a strong first-order phase transition from nucleonic to quark matter. Furthermore, we combine a covariant density-functional nucleonic equation of state (EOS) with a constant-speed-of-sound EOS for quark matter. We show that the radius and tidal deformability ranges obtained from GW170817 can be reconciled with the implication of the PREX-II experiment if there is a phase transition to quark matter in the low-mass compact star. In the high-mass segment, the EOS needs to be stiff to comply with the large-radius inference for PSR J0740 + 6620 and J0030 + 0451 with masses M ≃ 2 M ⊙ and M≃1.4 M⊙. We show that twin stars are not excluded, but the mass and radius ranges (with M ≥ M ⊙ ) are restricted to narrow domains ΔM twin ≲ 0.05 M ⊙ and ΔR twin ~ 1.0 km. We also show that the existence of twin configurations is compatible with the light companion in the GW190814 event being a hybrid star in the case of values of the sound-speed square s = 0.6 and s = 1/3.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Bayesian inference of multi-messenger astrophysical data: Joint and coherent inference of gravitational waves and kilonovae

Multi-messenger observations of binary neutron star mergers can provide information on the neutron star’s equation of state (EOS) above the nuclear saturation density by directly constraining the mass-radius diagram. We present a Bayesian framework for joint and coherent analyses of multi-messenger binary neutron star signals. As a first application, we analyze the gravitational-wave GW170817 and the kilonova (kN) AT2017gfo data. These results are then combined with the most recent X-ray pulsar analyses of PSR J0030+0451 and PSR J0740+6620 to obtain new EOS constraints.We extend the bajes infrastructure with a joint likelihood for multiple datasets, support for various semi-analytical kN models, and numerical-relativity (NR)-informed relations for the mass ejecta, as well as a technique to include and marginalize over modeling uncertainties. The analysis of GW170817 used the TEOBResumS effective-one-body waveform template to model the gravitational-wave signal. The analysis of AT2017gfo used a baseline multicomponent spherically symmetric model for the kN light curves. Various constraints on the mass-radius diagram and neutron star properties were then obtained by resampling over a set of ten million parameterized EOSs, which was built under minimal assumptions (general relativity and causality).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

A Quantum-mechanical Treatment of Electron Broadening in Strong Magnetic Fields

Abstract Spectral line-shape fitting is an extremely useful tool in determining the gravity of white dwarf stars. This method is so far limited to nonmagnetic white dwarfs largely because the theory of line broadening in high magnetic fields is not as complete as in the nonmagnetic case. Current Stark+Zeeman models treat plasma particles classically and ignore the motion of the nucleus. We develop the formalism for a quantum-mechanical treatment of the perturbing electrons and include the nuclear motion as part of the broadening and explore their relative importance. The conditions we explore are those found in white dwarf and neutron star atmospheres. We find that, contrary to previous studies, the quantized perturbing electrons create more broadening than perturbers on a straight-path trajectory. Additionally, the quantization of the plasma electrons gives rise to resonances away from the line center. The nuclear motion creates an additional electric field, which also leads to an increase in line broadening; however, this effect in neutron star atmospheres is not as large as previously estimated. This suggests that neutron star spectral lines are sensitive to density and that their mass and radius can be obtained from spectral line fitting, which would help constrain the neutron star equation of state.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Nuclear and hybrid equations of state in light of the low-mass compact star in HESS J1731-347

Here, we sample over 500 relativistic mean-field theories constrained by chiral effective-field theory and properties of isospin-symmetric nuclear matter and test them against known stellar structure constraints. This includes a recent mass and radius measurement of a compact object in supernova remnant HESS J1731-347, with an unusually low mass of M = 0.77$^{+0.20}_{–0.17}$M ⊙ and a compact radius of R = 10.4$^{+0.86}_{–0.78}$ km. We show that none of the sampled nuclear models meet all constraints at the 68% credibility level, but that hybrid equations of state with a quark-matter inner core and nuclear outer core easily can. This indicates a tension between astrophysical constraints and low-energy nuclear theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Non-local model for surface tension in fluid-fluid simulations

We propose a non-local model for surface tension obtained in the form of an integral of a molecular-force-like function with support 3.5e added to the Navier-Stokes momentum conservation equation. We demonstrate analytically and numerically that with the non-local model interfaces with a radius of curvature larger than the support length behave macroscopically and microscopically, otherwise. For static droplets, the pressure difference Pe,in - Pe,out satisfies the Young-Laplace law for droplet radius greater than 3.5e and other- wise deviates from the Young-Laplace law. The latter indicates that the surface tension in the proposed model decreases with decreasing radius of curvature, which agrees with molecular dynamics and experimental studies of nanodroplets. Using the non-local model we perform numerical simulations of droplets under dynamic conditions, including a rising droplet, a droplet in shear flow, and two colliding droplets in shear flow, and compare results with a standard Navier-Stokes model subject to the Young-Laplace boundary condition at the fluid-fluid interface implemented via the Conservative Level Set (CLS) method. We find good agreement with existing numerical methods and analytical results for a rising macroscopic droplet and a droplet in a shear flow. For colliding droplets in shear flow, the non- local model converges (with respect to the grid size) to the correct behavior, including sliding, coalescence, and merging and breaking of two droplets depending on the capillary number. In contrast, we find that the results of the CLS model are highly grid-size dependent.

Howard, Amanda A.↗
DOE OSTI2022

Structure and density of silicon carbide to 1.5 TPa and implications for extrasolar planets

There has been considerable recent interest in the high-pressure behavior of silicon carbide, a potential major constituent of carbon-rich exoplanets. In this work, the atomic-level structure of SiC was determined through in situ X-ray diffraction under laser-driven ramp compression up to 1.5 TPa; stresses more than seven times greater than previous static and shock data. Here we show that the B1-type structure persists over this stress range and we have constrained its equation of state (EOS). Using this data we have determined the first experimentally based mass-radius curves for a hypothetical pure SiC planet. Interior structure models are constructed for planets consisting of a SiC-rich mantle and iron-rich core. Carbide planets are found to be ~10% less dense than corresponding terrestrial planets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

KELT-9 b’s Asymmetric TESS Transit Caused by Rapid Stellar Rotation and Spin–Orbit Misalignment

KELT-9 b is an ultra-hot Jupiter transiting a rapidly rotating, oblate early-A-type star in a polar orbit. We model the effect of rapid stellar rotation on KELT-9 b’s transit light curve using photometry from the Transiting Exoplanet Survey Satellite to constrain the planet’s true spin–orbit angle and to explore how KELT-9 b may be influenced by stellar gravity darkening. We constrain the host star’s equatorial radius to be 1.089 ± 0.017 times as large as its polar radius and its local surface brightness to vary by ∼38% between its hot poles and cooler equator. We model the stellar oblateness and surface brightness gradient and find that it causes the transit light curve to lack the usual symmetry around the time of minimum light. We take advantage of the light-curve asymmetry to constrain KELT-9 b’s true spin–orbit angle (87{sup ∘}{sub −11{sup ∘+10{sup {sup ∘}}}}), agreeing with Gaudi et al. that KELT-9 b is in a nearly polar orbit. We also apply a gravity-darkening correction to the spectral energy distribution model from Gaudi et al. and find that accounting for rapid rotation gives a better fit to available spectroscopy and yields a more reliable estimate for the star’s polar effective temperature.

79 ASTRONOMY AND ASTROPHYSICS↗