Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “STARS”

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 757 records · Page 42

LONG-ORBITAL-PERIOD PREPOLARS CONTAINING EARLY K-TYPE DONOR STARS. BOTTLENECK ACCRETION MECHANISM IN ACTION

We studied two objects identified as cataclysmic variables (CVs) with periods exceeding the natural boundary for Roche-lobe-filling zero-age main sequence (ZAMS) secondary stars. We present observational results for V1082 Sgr with a 20.82 hr orbital period, an object that shows a low luminosity state when its flux is totally dominated by a chromospherically active K star with no signs of ongoing accretion. Frequent accretion shutoffs, together with characteristics of emission lines in a high state, indicate that this binary system is probably detached, and the accretion of matter on the magnetic white dwarf takes place through stellar wind from the active donor star via coupled magnetic fields. Its observational characteristics are surprisingly similar to V479 And, a 14.5 hr binary system. They both have early K-type stars as donor stars. We argue that, similar to the shorter-period prepolars containing M dwarfs, these are detached binaries with strong magnetic components. Their magnetic fields are coupled, allowing enhanced stellar wind from the K star to be captured and channeled through the bottleneck connecting the two stars onto the white dwarf’s magnetic pole, mimicking a magnetic CV. Hence, they become interactive binaries before they reach contact. This will help to explain an unexpected lack of systems possessing white dwarfs with strong magnetic fields among detached white+red dwarf systems.

79 ASTRONOMY AND ASTROPHYSICS↗

Evolution Models of Helium White Dwarf–Main-sequence Star Merger Remnants

It is predicted that orbital decay by gravitational-wave radiation and tidal interaction will cause some close binary stars to merge within a Hubble time. The merger of a helium-core white dwarf with a main-sequence (MS) star can produce a red giant branch star that has a low-mass hydrogen envelope when helium is ignited and thus become a hot subdwarf. Because detailed calculations have not been made, we compute post-merger models with a stellar evolution code. We find the evolutionary paths available to merger remnants and find the pre-merger conditions that lead to the formation of hot subdwarfs. We find that some such mergers result in the formation of stars with intermediate helium-rich surfaces. These stars later develop helium-poor surfaces owing to diffusion. Combining our results with a model population and comparing to observed stars, we find that some observed intermediate helium-rich hot subdwarfs can be explained as the remnants of the mergers of helium-core white dwarfs with low-mass MS stars.

79 ASTRONOMY AND ASTROPHYSICS↗

On the Origin of Sub-subgiant Stars. I. Demographics

Sub-subgiants are stars that are observed to be redder than normal main-sequence stars and fainter than normal subgiant (and giant) stars in an optical color–magnitude diagram (CMD). The red straggler stars, which lie redward of the red giant branch, may be related and are often grouped together with the sub-subgiants in the literature. These stars defy our standard theory of single-star evolution and are important tests for binary evolution and stellar collision models. In total, we identify 65 sub-subgiants (SSG) and red stragglers (RS) in 16 open and globular star clusters from the literature; 50 of these, including 43 sub-subgiants, pass our strict membership selection criteria (though the remaining sources may also be cluster members). In addition to their unique location on the CMD, we find that at least 58% (25/43) of sub-subgiants in this sample are X-ray sources with typical 0.5–2.5 keV luminosities of order 10{sup 30}–10{sup 31} erg s{sup −1}. Their X-ray luminosities and optical–to–X-ray flux ratios are similar to those of RS CVn active binaries. At least 65% (28/43) of the sub-subgiants in our sample are variables, 21 of which are known to be radial-velocity binaries. Typical variability periods are ≲15 days. At least 33% (14/43) of the sub-subgiants are H α emitters. These observational demographics provide strong evidence that binarity is important for sub-subgiant formation. Finally, we find that the number of sub-subgiants per unit mass increases toward lower-mass clusters, such that the open clusters in our sample have the highest specific frequencies of sub-subgiants.

79 ASTRONOMY AND ASTROPHYSICS↗

Star formation rate and stellar mass calibrations based on infrared photometry and their dependence on stellar population age and extinction

The stellar mass (M $\star$ ) and the star formation rate (SFR) are among the most important features that characterize galaxies. Measuring these fundamental properties accurately is critical for understanding the present state of galaxies, their history, and future evolution. Infrared (IR) photometry is widely used to measure the M $\star$ and SFR of galaxies because the near-IR traces the continuum emission of the majority of their stellar populations (SPs), and the mid/far-IR traces the dust emission powered by star-forming activity. This work explores the dependence of the IR emission of galaxies on their extinction, and the age of their SPs. It aims to provide accurate and precise IR-photometry SFR and M $\star$ calibrations that account for SP age and extinction while providing quantification of their scatter. We used the CIGALE spectral energy distribution (SED) fitting code to create model SEDs of galaxies with a wide range of star formation histories, dust content, and interstellar medium properties. We fit the relations between M $\star$ and SFR with IR and optical photometry of the model-galaxy SEDs with the Markov chain Monte Carlo (MCMC) method. As an independent confirmation of the MCMC fitting method, we performed a machine-learning random forest (RF) analysis on the same data set. The RF model yields similar results to the MCMC fits, thus validating the latter. This work provides calibrations for the SFR using a combination of the WISE bands 1 and 3, or the JWST NIR-F200W and MIRI-F2100W. It also provides mass-to-light ratio calibrations based on the WISE band-1, the JWST NIR-F200W, and the optical u - r or g - r colors. These calibrations account for the biases attributed to the SP age, while they are given in the form of extinction-dependent and extinction-independent relations. The proposed calibrations show robust estimations while minimizing the scatter and biases throughout a wide range of SFRs and stellar masses. The SFR calibration offers better results, especially in dust-free or passive galaxies where the contributions of old SPs or biases from the lack of dust are significant. Similarly, the M $\star$ calibration yields significantly better results for dusty and high-SFR galaxies where dust emission can otherwise bias the estimations.

79 ASTRONOMY AND ASTROPHYSICS↗

UCAC4 nearby star survey: A search for our stellar neighbors

We use data from the U.S. Naval Observatory fourth CCD Astrograph Catalog (UCAC4) in combination with photometry from the AAVSO Photometric All-Sky Survey and Two Micron All-Sky Survey to identify stars within 25 pc of the Sun. A sample of nearby stars with accurate trigonometric parallaxes from the Research Consortium On Nearby Stars is used to generate a set of 16 new photometric color–M{sub K{sub s}} relations that provide distance estimates with uncertainties of 15%. This work expands the available suites of well-calibrated photometric distance relations that can be used to identify nearby stellar systems. The distance relations are used with quality cuts to extract an initial sample of stars from the UCAC4 estimated to be within 25 pc. Color, proper motion, and existing literature sources are then used to obtain a clean sample of red dwarfs, while limiting the amount of contamination from background giants, resulting in a sample of 1761 candidate nearby stars within 25 pc. Of these, 339 are new discoveries with no previously known published parallax or distance estimate, primarily with proper motions less than 0.2 arcsec yr{sup −1}. Five stars are estimated to be within 10 pc, with the nearest, TYC 3980 1081 1, with V=10.50 estimated to be at 5.9 pc. That several hundred new stars have been revealed so close to the Sun illustrates once again that there is considerable work yet to be done to map the solar neighborhood and that additional nearby stars are likely still to be discovered.

79 ASTRONOMY AND ASTROPHYSICS↗

ROTATIONAL CORRECTIONS TO NEUTRON-STAR RADIUS MEASUREMENTS FROM THERMAL SPECTRA

We calculate the rotational broadening in the observed thermal spectra of neutron stars spinning at moderate rates in the Hartle-Thorne approximation. These calculations accurately account for the effects of the second-order Doppler boosts as well as for the oblate shapes and the quadrupole moments of the neutron stars. We find that fitting the spectra and inferring the bolometric fluxes under the assumption that a star is not rotating causes an underestimate of the inferred fluxes and, thus, radii. The correction depends on the stellar spin, mass, radius, and the observer's inclination. For a 10 km, 1.4 M {sub ☉} neutron star spinning at 600 Hz, the rotational correction to the flux is ∼1%-4%, while for a 15 km neutron star with the same spin period, the correction ranges from 2% for pole-on sources to 12% for edge-on sources. We calculate the inclination-averaged corrections to inferred radii as a function of the neutron-star radius and mass and provide an empirical formula for the corrections. For realistic neutron-star parameters (1.4 M {sub ☉}, 12 km, 600 Hz), the stellar radius is on the order of 4% larger than the radius inferred under the assumption that the star is not spinning.

79 ASTRONOMY AND ASTROPHYSICS↗

Star-forming galaxies as the origin of the ICECUBE PeV neutrinos

Star-forming galaxies, due to their high star formation rates, and hence large number of supernova remnants (SNRs) therein, are huge reservoirs of cosmic rays (CRs). These CRs collide with gases in galaxies and produce high-energy neutrinos through proton–proton collisions. In this paper, we calculate the neutrino production efficiency in star-forming galaxies by considering realistic galaxy properties, such as the gas density and galactic wind in star-forming galaxies. To calculate the accumulated neutrino flux, we use the infrared luminosity function of star-forming galaxies recently obtained by the Herschel PEP/HerMES survey. The intensity of CRs producing PeV neutrinos in star-forming galaxies is normalized with the observed CR flux at EeV (1 EeV = 10{sup 18} eV), assuming that SNR or hypernova remnants in star-forming galaxies can accelerate protons to EeV energies. Our calculations show that the accumulated neutrino emission produced by CRs in star-forming galaxies can account for the flux and spectrum of the sub-PeV/PeV neutrinos under reasonable assumptions on the CR confinement time in these galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

The fastest unbound stars in the universe

The discovery of hypervelocity stars (HVSs) leaving our galaxy with speeds of nearly 10{sup 3} km s{sup −1} has provided strong evidence of the existence of a massive compact object at the galaxy’s center. HVSs ejected via the disruption of stellar binaries can occasionally yield a star with v{sub ∞} ≲ 10{sup 4} km s{sup −1}; here we show that this mechanism can be extended to massive black hole (MBH) mergers, where the secondary star is replaced by a MBH with mass M{sub 2}≳10{sup 5}M{sub ⊙}. We find that stars that are originally bound to the secondary MBH are frequently ejected with v{sub ∞}>10{sup 4} km s{sup −1}, and occasionally with velocities ∼10{sup 5} km s{sup −1} (one third the speed of light). For this reason we refer to stars ejected from these systems as “semi-relativistic” hypervelocity stars (SHSs). Bound to no galaxy, the velocities of these stars are so great that they can cross a significant fraction of the observable universe in the time since their ejection (several Gpc). We demonstrate that if a significant fraction of MBH mergers undergoes a phase in which their orbital eccentricity is ≳0.5 and their periapse distance is tens of the primary’s Schwarzschild radius, the space density of fast-moving (v{sub ∞}>10{sup 4} km s{sup −1}) SHSs may be as large as 10{sup 3} Mpc{sup −3}. Hundreds of SHSs will be giant stars that can be detected by future all-sky infrared surveys such as WFIRST or Euclid and proper motion surveys such as LSST, with spectroscopic follow-up being possible with the James Webb Space Telescope.

79 ASTRONOMY AND ASTROPHYSICS↗

Linking galaxies to dark matter halos at z ∼ 1: Dependence of galaxy clustering on stellar mass and specific star formation rate

We study the dependence of angular two-point correlation functions on stellar mass (M{sub *}) and specific star formation rate (sSFR) of M{sub ∗}>10{sup 10}M{sub ⊙} galaxies at z∼1. The data from the UK Infrared Telescope Infrared Deep Sky Survey Deep eXtragalactic Survey and Canada–France–Hawaii Telescope Legacy Survey cover 8.2 deg{sup 2} sample scales larger than 100 h{sup −1} Mpc at z∼1, allowing us to investigate the correlation between clustering, M{sub *}, and star formation through halo modeling. Based on halo occupation distributions (HODs) of M{sub *} threshold samples, we derive HODs for M{sub *} binned galaxies, and then calculate the M{sub ∗}/M{sub halo} ratio. The ratio for central galaxies shows a peak at M{sub halo}∼10{sup 12}h{sup −1}M{sub ⊙}, and satellites predominantly contribute to the total stellar mass in cluster environments with M{sub ∗}/M{sub halo} values of 0.01–0.02. Using star-forming galaxies split by sSFR, we find that main sequence galaxies (log sSFR/yr{sup −1}∼−9) are mainly central galaxies in ∼10{sup 12.5}h{sup −1}M{sub ⊙} halos with the lowest clustering amplitude, while lower sSFR galaxies consist of a mixture of both central and satellite galaxies where those with the lowest M{sub *} are predominantly satellites influenced by their environment. Considering the lowest M{sub halo} samples in each M{sub *} bin, massive central galaxies reside in more massive halos with lower sSFRs than low mass ones, indicating star-forming central galaxies evolve from a low M{sub *}–high sSFR to a high M{sub *}–low sSFR regime. We also find that the most rapidly star-forming galaxies (log sSFR/yr{sup −1}>−8.5) are in more massive halos than main sequence ones, possibly implying galaxy mergers in dense environments are driving the active star formation. These results support the conclusion that the majority of star-forming galaxies follow secular evolution through the sustained but decreasing formation of stars.

79 ASTRONOMY AND ASTROPHYSICS↗

An upper bound on neutron star masses from models of short gamma-ray bursts

The discovery of two neutron stars with gravitational masses ≈2 M{sub ⊙} has placed a strong lower limit on the maximum mass of nonrotating neutron stars, and with it a strong constraint on the properties of cold matter beyond nuclear density. Current upper mass limits are much looser. Here, we note that if most short gamma-ray bursts are produced by the coalescence of two neutron stars, and if the merger remnant collapses quickly, then the upper mass limit is constrained tightly. If the rotation of the merger remnant is limited only by mass-shedding (which seems probable based on numerical studies), then the maximum gravitational mass of a nonrotating neutron star is ≈2−2.2 M{sub ⊙} if the masses of neutron stars that coalesce to produce gamma-ray bursts are in the range seen in Galactic double neutron star systems. These limits would be increased by ∼4% in the probably unrealistic case that the remnants rotate at ∼30% below mass-shedding, and by ∼15% in the extreme case that the remnants do not rotate at all. Future coincident detection of short gamma-ray bursts with gravitational waves will strengthen these arguments because they will produce tight bounds on the masses of the components for individual events. If these limits are accurate, then a reasonable fraction of double neutron star mergers might not produce gamma-ray bursts. In that case, or in the case that many short bursts are produced instead by the mergers of neutron stars with black holes, the implied rate of gravitational wave detections will be increased.

79 ASTRONOMY AND ASTROPHYSICS↗

TECHNIQUES FOR HIGH-CONTRAST IMAGING IN MULTI-STAR SYSTEMS. I. SUPER-NYQUIST WAVEFRONT CONTROL

Direct imaging of extra-solar planets is now a reality with the deployment and commissioning of the first generation of specialized ground-based instruments (GPI, SPHERE, P1640, and SCExAO). These systems allow of planets 10{sup 7} times fainter than their host star. For space-based missions (EXCEDE, EXO-C, EXO-S, WFIRST), various teams have demonstrated laboratory contrasts reaching 10{sup −10} within a few diffraction limits from the star. However, all of these current and future systems are designed to detect faint planets around a single host star, while most non-M-dwarf stars such as Alpha Centauri belong to multi-star systems. Direct imaging around binaries/multiple systems at a level of contrast allowing detection of Earth-like planets is challenging because the region of interest is contaminated by the host star's companion in addition to the host itself. Generally, the light leakage is caused by both diffraction and aberrations in the system. Moreover, the region of interest usually falls outside the correcting zone of the deformable mirror (DM) with respect to the companion. Until now, it has been thought that removing the light of a companion star is too challenging, leading to the exclusion of many binary systems from target lists of direct imaging coronographic missions. In this paper, we will show new techniques for high-contrast imaging of planets around multi-star systems and detail the Super-Nyquist Wavefront Control (SNWC) method, which allows wavefront errors to be controlled beyond the nominal control region of the DM. Our simulations have demonstrated that, with SNWC, raw contrasts of at least 5 × 10{sup −9} in a 10% bandwidth are possible.

79 ASTRONOMY AND ASTROPHYSICS↗

Extracting radial velocities of A- AND B-type stars from ECHELLE spectrograph calibration spectra

We present a technique to extract radial velocity (RV) measurements from echelle spectrograph observations of rapidly rotating stars (V sin i≳50 km s{sup −1}). This type of measurement is difficult because the line widths of such stars are often comparable to the width of a single echelle order. To compensate for the scarcity of lines and Doppler information content, we have developed a process that forward-models the observations, fitting the RV shift of the star for all echelle orders simultaneously with the echelle blaze function. We use our technique to extract RV measurements from a sample of rapidly rotating A- and B-type stars used as calibrator stars observed by the California Planet Survey observations. We measure absolute RVs with a precision ranging from 0.5–2.0 km s{sup −1} per epoch for more than 100 A- and B-type stars. In our sample of 10 well-sampled stars with RV scatter in excess of their measurement uncertainties, three of these are single-lined binaries with long observational baselines. From this subsample, we present detections of two previously unknown spectroscopic binaries and one known astrometric system. Our technique will be useful in measuring or placing upper limits on the masses of sub-stellar companions discovered by wide-field transit surveys, and conducting future spectroscopic binarity surveys and Galactic space-motion studies of massive and/or young, rapidly rotating stars.

79 ASTRONOMY AND ASTROPHYSICS↗

Gravity-mode period spacings as a seismic diagnostic for a sample oF γ DORADUS stars from Kepler space photometry and high-resolution ground-based spectroscopy

Gamma Doradus stars (hereafter γ Dor stars) are gravity-mode pulsators of spectral type A or F. Such modes probe the deep stellar interior, offering a detailed fingerprint of their structure. Four-year high-precision space-based Kepler photometry of γ Dor stars has become available, allowing us to study these stars with unprecedented detail. We selected, analyzed, and characterized a sample of 67 γ Dor stars for which we have Kepler observations available. For all the targets in the sample we assembled high-resolution spectroscopy to confirm their F-type nature. We found fourteen binaries, among which are four single-lined binaries, five double-lined binaries, two triple systems, and three binaries with no detected radial velocity variations. We estimated the orbital parameters whenever possible. For the single stars and the single-lined binaries, fundamental parameter values were determined from spectroscopy. We searched for period spacing patterns in the photometric data and identified this diagnostic for 50 of the stars in the sample, 46 of which are single stars or single-lined binaries. We found a strong correlation between the spectroscopic vsini and the period spacing values, confirming the influence of rotation on γ Dor-type pulsations as predicted by theory. We also found relations between the dominant g-mode frequency, the longest pulsation period detected in series of prograde modes, vsini, and log T{sub eff}.

79 ASTRONOMY AND ASTROPHYSICS↗

THE HABITABLE ZONES OF PRE-MAIN-SEQUENCE STARS

We calculate the pre-main-sequence habitable zone (HZ) for stars of spectral classes F-M. The spatial distribution of liquid water and its change during the pre-main-sequence phase of protoplanetary systems is important for understanding how planets become habitable. Such worlds are interesting targets for future missions because the coolest stars could provide habitable conditions for up to 2.5 billion years post-accretion. Moreover, for a given star type, planetary systems are more easily resolved because of higher pre-main-sequence stellar luminosities, resulting in larger planet-star separation for cool stars than is the case for the traditional main-sequence (MS) HZ. We use one-dimensional radiative-convective climate and stellar evolutionary models to calculate pre-main-sequence HZ distances for F1-M8 stellar types. We also show that accreting planets that are later located in the traditional MS HZ orbiting stars cooler than a K5 (including the full range of M stars) receive stellar fluxes that exceed the runaway greenhouse threshold, and thus may lose substantial amounts of water initially delivered to them. We predict that M-star planets need to initially accrete more water than Earth did, or, alternatively, have additional water delivered later during the long pre-MS phase to remain habitable. Our findings are also consistent with recent claims that Venus lost its water during accretion.

79 ASTRONOMY AND ASTROPHYSICS↗

Variable stars in Local Group galaxies – V. The fast and early evolution of the low-mass Eridanus II dSph galaxy

ABSTRACT We present a detailed study of the variable star population of Eridanus II (Eri II), an ultra-faint dwarf galaxy that lies close to the Milky Way virial radius. We analyse multi-epoch g, r, i ground-based data from Goodman and the Dark Energy Camera, plus F475W, F606W, F814W space data from the Advanced Camera for Surveys. We report the detection of 67 RR Lyrae (RRL) stars and 2 Anomalous Cepheids, most of them new discoveries. With the RRL stars, we measure the distance modulus of Eri II, μ0 = 22.84 ± 0.05 mag (D⊙ = 370 ± 9 kpc), and derive a metallicity spread of 0.3 dex (0.2 dex intrinsic). The colour distribution of the horizontal branch (HB) and the period distribution of the RRL stars can be nicely reproduced by a combination of two stellar models of [Fe/H] = (−2.62, −2.14). The overall low metallicity is consistent with the red giant branch bump location, 0.65 mag brighter than the HB. These results are in agreement with previous spectroscopic studies. The more metal-rich RRL and the RRab stars have greater central concentration than the more metal-poor RRL and the RRc stars that are mainly located outside ∼1 rh. This is similar to what is found in larger dwarf galaxies such as Sculptor, and in agreement with an outside-in galaxy formation scenario. This is remarkable in such a faint dwarf galaxy with an apparently single and extremely short (<1 Gyr) star formation burst. Finally, we have derived new and independent structural parameters for Eri II and its star cluster using our new data that are in very good agreement with previous estimates.

79 ASTRONOMY AND ASTROPHYSICS↗

Spectroscopic survey of faint planetary-nebula nuclei – II. The subdwarf O central star of Fr 2-30

Fr 2-30 = PN? G126.8−15.5 is a faint emission nebula, hosting a 14th-mag central star that we identify here for the first time. Deep Hα and [O III ] images reveal a roughly elliptical nebula with dimensions of at least 22 arcmin × 14 arcmin, fading into a surrounding network of even fainter emission. Optical spectrograms of the central star show it to have a subdwarf O spectral type, with a Gaia parallax distance of 890 pc. A model-atmosphere analysis gives parameters of T eff = 60 000 K, log g = 6.0, and a low helium content of n He /n H = 0.0017. The location of the central star in the log g–T eff plane is inconsistent with a post-asymptotic-giant-branch evolutionary status. Two alternatives are that it is a helium-burning post-extreme-horizontal-branch object, or a hydrogen-burning post-red-giant-branch star. In either case, the evolutionary ages are so long that a detectable planetary nebula (PN) should not be present. We find evidence for a variable radial velocity (RV), suggesting that the star is a close binary. However, there are no photometric variations, and the spectral-energy distribution rules out a companion earlier than M2 V. The RVs of the star and surrounding nebula are discordant, and the nebula lacks typical PN morphology. We suggest that Fr 2-30 is a ‘PN mimic’ – the result of a chance encounter between the hot sdO star and an interstellar cloud. However, we note the puzzling fact that there are several nuclei of genuine PNe that are known to be in evolutionary states similar to that of the Fr 2-30 central star.

79 ASTRONOMY AND ASTROPHYSICS↗

Neutron stars with exceptionally light QCD axions

We present a comprehensive study of axion condensed neutron stars that arise in models of an exceptionally light axion that couples to quantum chromodynamics (QCD). These axions solve the strong-charge-parity (𝐶⁢𝑃) problem, but have a mass-squared lighter than that due to QCD by a factor of 𝜖 < 1. Inside dense matter, the axion potential is altered, and much of the matter in neutron stars resides in the axion condensed phase where the strong-𝐶⁢𝑃 parameter 𝜃 = 𝜋 and 𝐶⁢𝑃 remains a good symmetry. In these regions, masses and interactions of nucleons are modified, in turn changing the equation of state (EOS), structure, and phenomenology of the neutron stars. We take the first steps toward the study of the EOS of neutron star matter at 𝜃 = 𝜋 within chiral effective field theory and use relativistic mean field theory to deduce the resulting changes to nuclear matter and the neutron star low-density EOS. We derive constraints on the exceptionally light axion parameter space based on observations of the thermal relaxation of accreting neutron stars, isolated neutron star cooling, and pulsar glitches, excluding the region up to 5 ×10 −7 ≲ 𝜖 ≲ 0.2 for 𝑚 𝑎 ≳ 2 × 10 −9 eV. We comment on potential changes to the neutron star mass-radius relationship, and discuss the possibility of novel, nuclear-density compact objects with 𝜃 = 𝜋 that are stabilized not by gravity but by the axion potential.

axions↗

Unified nonparametric equation-of-state inference from the neutron-star crust to perturbative-QCD densities

Perturbative quantum chromodynamics (pQCD), while valid only at densities exceeding those found in the cores of neutron stars, could provide constraints on the dense-matter equation of state (EOS). Here, in this work, we examine the impact of pQCD information on the inference of the EOS using a nonparametric framework based on Gaussian processes (GPs). We examine the application of pQCD constraints through a ``pQCD likelihood,'' and verify the findings of previous works; namely, a softening of the EOS at the central densities of the most massive neutron stars and a reduction in the maximum neutron-star mass. Although the pQCD likelihood can be easily integrated into existing EOS inference frameworks, this approach requires an arbitrary selection of the density at which the constraints are applied. The EOS behavior is also treated differently on either side of the chosen density. To mitigate these issues, we extend the EOS model to higher densities, thereby constructing a ``unified'' description of the EOS from the neutron-star crust to densities relevant for pQCD. In this approach the pQCD constraints effectively become part of the prior. Since the EOS is unconstrained by any calculation or data between the densities applicable to neutron stars and pQCD, we argue for maximum modeling flexibility in that regime. We compare the unified EOS with the traditional pQCD likelihood, and although we confirm the EOS softening, we do not see a reduction in the maximum neutron-star mass or any impact on macroscopic observables. Though residual model dependence cannot be ruled out, we find that pQCD suggests the speed of sound in the densest neutron-star cores has already started decreasing toward the asymptotic limit; we find that the speed of sound squared at the center of the most massive neutron star has an upper bound of $\sim 0.5$ at the 90% level.

equations of state of nuclear matter↗