Lanthanide/Actinide Opacities for Neutron Star Merger Models
Explore the source record for details and available documents.
Hungerford, Aimee L.↗
SEARCH · Search NASA
Results for “STAR MODELS”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Modeling Neutron Star Oscillations in a Fixed General Relativistic Background Including Solid Crust Dynamics
Explore the source record for details and available documents.
Tsao, Bing-Jyun↗
The age of the directly imaged planet host star κ Andromedae determined from interferometric observations
κ Andromedae, an early-type star that hosts a directly imaged low-mass companion, is expected to be oblate due to its rapid rotational velocity (v sin i = ∼162 km s{sup −1}). We observed the star with the CHARA Array’s optical beam combiner, PAVO, measuring its size at multiple orientations and determining its oblateness. The interferometric measurements, combined with photometry and this v sin i value are used to constrain an oblate star model that yields the fundamental properties of the star and finds a rotation speed that is ∼85% of the critical rate and a low inclination of ∼30°. Three modeled properties (the average radius, bolometric luminosity, and equatorial velocity) are compared to MESA evolution models to determine an age and mass for the star. In doing so, we determine an age for the system of 47{sub −40}{sup +27} Myr. Based on this age and previous measurements of the companion’s temperature, the BHAC15 evolution models imply a mass for the companion of 22{sub −9}{sup +8} M {sub J}.
79 ASTRONOMY AND ASTROPHYSICS↗
Star-forming Filament Models
New models of star-forming filamentary clouds are presented in order to quantify their properties and to predict their evolution. These 2D axisymmetric models describe filaments that have no core, one low-mass core, and one cluster-forming core. They are based on Plummer-like cylinders and spheroids that are bounded by a constant-density surface of finite extent. In contrast to 1D Plummer-like models, they have specific values of length and mass, they approximate observed column density maps, and their distributions of column density ( N -pdfs) are pole-free. Each model can estimate the star-forming potential of a core-filament system by identifying the zone of gas dense enough to form low-mass stars and by counting the number of enclosed thermal Jeans masses. This analysis suggests that the Musca central filament may be near the start of its star-forming life, with enough dense gas to make its first ∼3 protostars, while the Coronet filament is near the midpoint of its star formation, with enough dense gas to add ∼8 protostars to its ∼20 known stars. In contrast, L43 appears to be near the end of its star-forming life, since it lacks enough dense gas to add any new protostars to the two young stellar objectsalready known.
79 ASTRONOMY AND ASTROPHYSICS↗
THE STELLAR OBLIQUITY, PLANET MASS, AND VERY LOW ALBEDO OF QATAR-2 FROM K2 PHOTOMETRY
The Qatar-2 transiting exoplanet system was recently observed in short-cadence mode by Kepler as part of K2 Campaign 6. We identify dozens of starspot-crossing events, when the planet eclipses a relatively dark region of the stellar photosphere. The observed patterns of these events demonstrate that the planet always transits over the same range of stellar latitudes and, therefore, that the stellar obliquity is less than about 10°. We support this conclusion with two different modeling approaches: one based on explicit identification and timing of the events and the other based on fitting the light curves with a spotted-star model. We refine the transit parameters and measure the stellar rotation period (18.5 ± 1.9 days), which corresponds to a “gyrochronological” age of 1.4 ± 0.3 Gyr. Coherent flux variations with the same period as the transits are well modeled as the combined effects of ellipsoidal light variations (15.4 ± 4.8 ppm) and Doppler boosting (14.6 ± 5.1 ppm). The magnitudes of these effects correspond to a planetary mass of 2.6±0.9 M{sub Jup} and 3.9±1.5 M{sub Jup}, respectively. Both of these independent mass estimates agree with the mass determined by the spectroscopic Doppler technique (2.487±0.086 M{sub Jup}). No occultations are detected, giving a 2σ upper limit of 0.06 on the planet’s visual geometric albedo. We find no evidence for orbital decay, although we are only able to place a weak lower bound on the relevant tidal quality factor: Q{sub ⋆}{sup ′}>1.5×10{sup 4} (95% confidence).
79 ASTRONOMY AND ASTROPHYSICS↗
Two-dimensional simulations of internal gravity waves in a $\mathrm{5M_⊙}$ zero-age-main-sequence model
Main-sequence intermediate-mass stars present a radiative envelope that supports internal gravity waves (IGWs). Excited at the boundary with the convective core, IGWs propagate towards the stellar surface and are suspected to impact physical processes such as rotation and chemical mixing. Using the fully compressible time-implicit code MUSIC, we study IGWs in two-dimensional simulations of a zero-age-main-sequence 5 solar mass star model up to 91 percent of the stellar radius with different luminosity and radiative diffusivity enhancements. Our results show that low-frequency waves excited by core convection are strongly impacted by radiative effects as they propagate. This impact depends on the radial profile of radiative diffusivity which increases by almost 5 orders of magnitude between the centre of the star and the top of the simulation domain. In the upper layers of the simulation domain, we observe an increase of the temperature. Our study suggests that this is due to heat added in these layers by IGWs damped by radiative diffusion. We show that non-linear effects linked to large amplitude IGWs may be relevant just above the convective core. Both these effects are intensified by the artificial enhancement of the luminosity and radiative diffusivity, with enhancement factors up to 10 4 times the realistic values. Our results also highlight that direct comparison between numerical simulations with enhanced luminosity and observations must be made with caution. Finally, our work suggests that thermal effects linked to the damping of IGWs could have a non-negligible impact on stellar structure.
79 ASTRONOMY AND ASTROPHYSICS↗
Assessing equation of state-independent relations for neutron stars with nonparametric models
Relations between neutron star properties that do not depend on the nuclear equation of state offer insights on neutron star physics and have practical applications in data analysis. Such relations are obtained by fitting to a range of phenomenological or nuclear physics equation of state models, each of which may have varying degrees of accuracy. In this study we revisit commonly used relations and reassess them with a very flexible set of phenomenological nonparametric equation of state models that are based on Gaussian processes. Our models correspond to two sets: equations of state which mimic hadronic models, and equations of state with rapidly changing behavior that resemble phase transitions. Here we quantify the accuracy of relations under both sets and discuss their applicability with respect to expected upcoming statistical uncertainties of astrophysical observations. We further propose a goodness-of-fit metric which provides an estimate for the systematic error introduced by using the relation to model a certain equation-of-state set. Overall, the nonparametric distribution is more poorly fit with existing relations, with the I–Love–Q relations retaining the highest degree of universality. Fits degrade for relations involving the tidal deformability, such as the binary-Love and compactness-Love relations, and when introducing phase transition phenomenology. For most relations, systematic errors are comparable to current statistical uncertainties under the nonparametric equation of state distributions.
79 ASTRONOMY AND ASTROPHYSICS↗
Simultaneous Constraints on the Mass and Radius of Aql X–1 from Quiescence and X-Ray Burst Observations
The measurement of neutron star mass and radius is one of the most direct ways to distinguish between various dense matter equations of state. The mass and radius of accreting neutron stars hosted in low-mass X-ray binaries can be constrained by several methods, including photospheric radius expansion from type I X-ray bursts and from quiescent spectra. In this paper, we apply for the first time these two methods simultaneously to constrain the mass and radius of Aql X–1. The quiescent spectra from Chandra and XMM-Newton , and photospheric radius expansion bursts from RXTE are used. The determination of the mass and radius of Aql X–1 is also used to verify the consistency between the two methods and to narrow down the uncertainties of the neutron star mass and radius. It is found that the distance to Aql X–1 should be in the range of 4.0–5.75 kpc, based on the overlapping confidence regions between photospheric radius expansion burst and quiescent spectra methods. In addition, we show that the mass and radius determined for the compact star in Aql X–1 are compatible with strange star equations of state and conventional neutron star models.
79 ASTRONOMY AND ASTROPHYSICS↗
Physical Properties of the Star-forming Regions in the Interacting Galaxies NGC 3395/NGC 3396
We report long-slit spectroscopy of the early major merger galaxies NGC 3395/NGC 3396. The spectra are consistent with those for star-forming galaxies, but there is some indication of LINER-like active galactic nucleus activity in the center of NGC 3396. The total star formation rate in the regions observed is 2.83 M {sub ⊙} yr{sup −1}, consistent with estimates for the entire galaxies. The highest abundances are in the centers of the galaxies, with the abundances decreasing with distance. There is a correlation between high abundance and high ionization parameter, both of which can be attributed to the presence of massive stars. Modeling with SB 99 indicates the star-forming regions are younger than 10 Myr. There are 1000–2000 WNL stars in the system, along with several thousand O stars, consistent with the ages of the star-forming regions. The highest electron densities are found in young regions with high star formation rates. The electron temperatures are higher than results for non-interacting galaxies, which is probably due to shock waves produced by the galaxy–galaxy interaction, the outflow of gas from massive stars, and/or collisions between gas clouds in the galaxies. There is star formation in the bridge of material between the galaxies. These regions are among the youngest in the system and have low abundances, suggesting the gas was pulled from the outer parts of the galaxies. X-ray point sources, probably high-mass X-ray binaries, are associated with several star-forming regions.
79 ASTRONOMY AND ASTROPHYSICS↗
Millihertz quasi-periodic oscillations in 4U 1636-536: putting possible constraints on the neutron star size
Based on previous studies of quasi-periodic oscillations (QPOs) in neutron star (NS) low-mass X-ray binaries, mHz QPOs are believed to be related to “marginally stable” burning on the NS surface. Our study of phase-resolved energy spectra of these oscillations in 4U 1636–53 shows that the oscillations are not caused by variations in the blackbody temperature of the NS, but reveals a correlation between the change of the count rate during the mHz QPO pulse and the spatial extent of a region emitting blackbody emission. The maximum size of the emission area, R{sub BB}{sup 2}=216.7{sub −86.4}{sup +93.2} km{sup 2}, provides direct evidence that the oscillations originate from a variable surface area constrained on the NS and are therefore not related to instabilities in the accretion disk. The obtained lower limit on the size of the NS (11.0 km) rules out equations of state that prefer small NS radii. Observations of mHz QPOs in NS LMXBs with NICER and eXTP will reduce the statistical uncertainty in the lower limit on the NS radius, which together with better estimates of the hardening factor and distance, will allow for improved discrimination between different equations of state and compact star models. Furthermore, future missions will allow us to measure the peak blackbody emission area for a single mHz QPO pulse, which will push the lower limit to larger radii.
79 ASTRONOMY AND ASTROPHYSICS↗
On the spin of the black hole in IC 10 X–1
The compact X-ray source in the eclipsing X-ray binary IC 10 X–1 has reigned for years as ostensibly the most massive stellar-mass black hole, with a mass estimated to be about twice that of its closest rival. However, striking results presented recently by Laycock et al. reveal that the mass estimate, based on emission-line velocities, is unreliable and that the mass of the X-ray source is essentially unconstrained. Using Chandra and NuSTAR data, we rule against a neutron-star model and conclude that IC 10 X–1 contains a black hole. The eclipse duration of IC 10 X–1 is shorter and its depth shallower at higher energies, an effect consistent with the X-ray emission being obscured during eclipse by a Compton-thick core of a dense wind. The spectrum is strongly disk-dominated, which allows us to constrain the spin of the black hole via X-ray continuum fitting. Three other wind-fed black hole systems are known; the masses and spins of their black holes are high: M∼10--15M{sub ⊙} and a{sub ∗}>0.8. If the mass of IC 10 X–1's black hole is comparable, then its spin is likewise high.
79 ASTRONOMY AND ASTROPHYSICS↗
Rapidly rotating neutron stars with a massive scalar field—structure and universal relations
We construct rapidly rotating neutron star models in scalar-tensor theories with a massive scalar field. The fact that the scalar field has nonzero mass leads to very interesting results since the allowed range of values of the coupling parameters is significantly broadened. Deviations from pure general relativity can be very large for values of the parameters that are in agreement with the observations. We found that the rapid rotation can magnify the differences several times compared to the static case. The universal relations between the normalized moment of inertia and quadrupole moment are also investigated both for the slowly and rapidly rotating cases. The results show that these relations are still EOS independent up to a large extend and the deviations from pure general relativity can be large. This places the massive scalar-tensor theories amongst the few alternative theories of gravity that can be tested via the universal I -Love- Q relations.
79 ASTRONOMY AND ASTROPHYSICS↗
New measurements of low-energy resonances in the Ne 22 ( p , γ ) Na 23 reaction
The Ne 22 ( p , γ ) Na 23 reaction is one of the most uncertain reactions in the NeNa cycle and plays a crucial role in the creation of Na 23 , the only stable Na isotope. Uncertainties in the low-energy rates of this and other reactions in the NeNa cycle lead to ambiguities in the nucleosynthesis predicted from models of thermally pulsing asymptotic giant branch (AGB) stars. This in turn complicates the interpretation of anomalous Na-O trends in globular cluster evolutionary scenarios. Prior studies of the Ne 22 ( p , γ ) Na 23 , Ne 22 ( He 3 , d ) Na 23 , and C 12 ( C 12 , p ) Na 23 reactions disagree on the strengths, spins, and parities of low-energy resonances in Na 23 and the direct-capture Ne 22 ( p , γ ) Na 23 reaction rate contains large uncertainties as well. Here, we present new measurements of resonances at E r c.m. = 417 , 178, and 151 keV and of the direct-capture process in the Ne 22 ( p , γ ) Na 23 reaction. The resulting total Ne 22 ( p , γ ) Na 23 rate is approximately a factor of 20 higher than the rate listed in a recent compilation at temperatures relevant to hot-bottom burning in AGB stars. Although our rate is close to that derived from a recent Ne 22 ( p , γ ) Na 23 measurement by Cavanna et al. in 2015, we find that this large rate increase results in only a modest 18% increase in the Na 23 abundance predicted from a 5 M ⊙ thermally pulsing AGB star model from Ventura and D'Antona (2005). The estimated astrophysical impact of this rate increase is in marked contrast to the factor of ~ 3 increase in Na 23 abundance predicted by Cavanna et al. and is attributed to the interplay between the Na 23 ( p , α ) Ne 20 and Ne 20 ( p , γ ) Na 21 reactions, both of which remain fairly uncertain at the relevant temperature range.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Application of a theory and simulation-based convective boundary mixing model for AGB star evolution and nucleosynthesis
The s-process nucleosynthesis in Asymptotic giant branch (AGB) stars depends on the modeling of convective boundaries. We present models and s-process simulations that adopt a treatment of convective boundaries based on the results of hydrodynamic simulations and on the theory of mixing due to gravity waves in the vicinity of convective boundaries. Hydrodynamics simulations suggest the presence of convective boundary mixing (CBM) at the bottom of the thermal pulse-driven convective zone. Similarly, convection-induced mixing processes are proposed for the mixing below the convective envelope during third dredge-up (TDU), where the {sup 13}C pocket for the s process in AGB stars forms. In this work, we apply a CBM model motivated by simulations and theory to models with initial mass M = 2 and M=3 M{sub ⊙}, and with initial metal content Z = 0.01 and Z = 0.02. As reported previously, the He-intershell abundances of {sup 12}C and {sup 16}O are increased by CBM at the bottom of the pulse-driven convection zone. This mixing is affecting the {sup 22}Ne(α, n){sup 25}Mg activation and the s-process efficiency in the {sup 13}C-pocket. In our model, CBM at the bottom of the convective envelope during the TDU represents gravity wave mixing. Furthermore, we take into account the fact that hydrodynamic simulations indicate a declining mixing efficiency that is already about a pressure scale height from the convective boundaries, compared to mixing-length theory. We obtain the formation of the {sup 13}C-pocket with a mass of ≈10{sup −4} M{sub ⊙}. The final s-process abundances are characterized by 0.36<[s/Fe]<0.78 and the heavy-to-light s-process ratio is −0.23<[hs/ls]<0.45. Finally, we compare our results with stellar observations, presolar grain measurements and previous work.
79 ASTRONOMY AND ASTROPHYSICS↗
Gravitational waves from binary neutron star mergers with a spectral equation of state
In numerical simulations of binary neutron star systems, the equation of state of the dense neutron star matter is an important factor in determining both the physical realism and the numerical accuracy of the simulations. Some equations of state used in simulations are C 2 or smoother in the pressure/density relationship function, such as a polytropic equation of state, but may not have the flexibility to model stars or remnants of different masses while keeping their radii within known astrophysical constraints. Other equations of state, such as tabular or piece-wise polytropic, may be flexible enough to model additional physics and multiple stars' masses and radii within known constraints, but are not as smooth, resulting in additional numerical error. We will study in this paper a recently developed family of equation of state, using a spectral expansion with sufficient free parameters to allow for a larger flexibility than current polytropic equations of state, and with sufficient smoothness to reduce numerical errors compared to tabulated or piece-wise polytropic equations of state. We perform simulations at three mass ratios with a common chirp mass, using two distinct spectral equations of state, and at multiple numerical resolutions. We evaluate the gravitational waves produced from these simulations, comparing the phase error between resolutions and equations of state, as well as with respect to analytical models. From our simulations we estimate that the phase difference at merger for binaries with a dimensionless weighted tidal deformability difference greater than Δ$\tilde{Λ}$=55 can be captured by the SpEC code for these equations of state.
79 ASTRONOMY AND ASTROPHYSICS↗
Three-dimensional Hydrodynamic Simulations of Convective Nuclear Burning in Massive Stars Near Iron Core Collapse
Nonspherical structure in massive stars at the point of iron core collapse can have a qualitative impact on the properties of the ensuing core-collapse supernova explosions and the multimessenger signals they produce. Strong perturbations can aid successful explosions by strengthening turbulence in the postshock region. Here we report on a set of 4π 3D hydrodynamic simulations of O- and Si-shell burning in massive star models of varied initial masses using MESA and the FLASH simulation framework. We evolve four separate 3D models for roughly the final 10 minutes prior to and including iron core collapse. We consider initial 1D MESA models with masses of 14, 20, and 25 M ⊙ to survey a range of O/Si-shell density and compositional configurations. Further, we characterize the convective shells in our 3D models and compare them to the corresponding 1D models. In general, we find that the angle-average convective speeds in our 3D simulations near collapse are three to four times larger than the convective speeds predicted by MESA at the same epoch for our chosen mixing length parameter of α MLT = 1.5. In three of our simulations, we observe significant power in the spherical harmonic decomposition of the radial velocity field at harmonic indices of ℓ = 1–3 near collapse. Our results suggest that large-scale modes are common in massive stars near collapse and should be considered a key aspect of presupernova progenitor models.
79 ASTRONOMY AND ASTROPHYSICS↗
A Neutron Star Binary Merger Model for GW170817/GRB 170817A/SSS17a
The merging neutron star gravitational-wave event GW170817 has been observed throughout the entire electromagnetic spectrum from radio waves to γ -rays. The resulting energetics, variability, and light curves are shown to be consistent with GW170817 originating from the merger of two neutron stars, in all likelihood followed by the prompt gravitational collapse of the massive remnant. The available γ -ray, X-ray, and radio data provide a clear probe for the nature of the relativistic ejecta and the non-thermal processes occurring within, while the ultraviolet, optical, and infrared emission are shown to probe material torn during the merger and subsequently heated by the decay of freshly synthesized r -process material. The simplest hypothesis, that the non-thermal emission is due to a low-luminosity short γ -ray burst (sGRB), seems to agree with the present data. While low-luminosity sGRBs might be common, we show here that the collective prompt and multi-wavelength observations are also consistent with a typical, powerful sGRB seen off-axis. Detailed follow-up observations are thus essential before we can place stringent constraints on the nature of the relativistic ejecta in GW170817.
79 ASTRONOMY AND ASTROPHYSICS↗
Fully compressible simulations of waves and core convection in main-sequence stars
Context: Recent, nonlinear simulations of wave generation and propagation in full-star models have been carried out in the anelastic approximation using spectral methods. Although it makes long time steps possible, this approach excludes the physics of sound waves completely and requires rather high artificial viscosity and thermal diffusivity for numerical stability. A direct comparison with observations is thus limited. Aims: In this study, we explore the capabilities of our compressible multidimensional Seven-League Hydro (SLH) code to simulate stellar oscillations. Methods: We compare some fundamental properties of internal gravity and pressure waves in 2D SLH simulations to linear wave theory using two test cases: (1) an interval gravity wave packet in the Boussinesq limit and (2) a realistic 3 M ⊙ stellar model with a convective core and a radiative envelope. Oscillation properties of the stellar model are also discussed in the context of observations. Results: Our tests show that specialized low-Mach techniques are necessary when simulating oscillations in stellar interiors. Basic properties of internal gravity and pressure waves in our simulations are in good agreement with linear wave theory. As compared to anelastic simulations of the same stellar model, we can follow internal gravity waves of much lower frequencies. The temporal frequency spectra of velocity and temperature are flat and compatible with the observed spectra of massive stars. Conclusion: The low-Mach compressible approach to hydrodynamical simulations of stellar oscillations is promising. Our simulations are less dissipative and require less luminosity boosting than comparable spectral simulations. The fully-compressible approach allows for the coupling of gravity and pressure waves in the outer convective envelopes of evolved stars to be studied in the future.
79 ASTRONOMY AND ASTROPHYSICS↗