Mass loss processes are a key uncertainty in the evolution of massive stars. They determine the amount of mass and angular momentum retained by the star, thus influencing its evolution and ...presupernova structure. Because of the high complexity of the physical processes driving mass loss, stellar evolution calculations must employ parametric algorithms, and usually only include wind mass loss. We carried out an extensive parameter study of wind mass loss and its effects on massive star evolution using the open-source stellar evolution code MESA. We provide a systematic comparison of wind mass loss algorithms for solar-metallicity, nonrotating, single stars in the initial mass range of 15 M⊙ to 35 M⊙. We consider combinations drawn from two hot phase (i.e., roughly the main sequence) algorithms, three cool phase (i.e., post-main-sequence) algorithms, and two Wolf-Rayet mass loss algorithms. We discuss separately the effects of mass loss in each of these phases. In addition, we consider linear wind efficiency scale factors of 1, 0.33, and 0.1 to account for suggested reductions in mass loss rates due to wind inhomogeneities. We find that the initial to final mass mapping for each zero-age main-sequence (ZAMS) mass has a ~ 50% uncertainty if all algorithm combinations and wind efficiencies are considered. The ad-hoc efficiency scale factor dominates this uncertainty. While the final total mass and internal structure of our models vary tremendously with mass loss treatment, final luminosity and effective temperature are much less sensitive for stars with ZAMS mass ≲ 30 M⊙. This indicates that uncertainty in wind mass loss does not negatively affect estimates of the ZAMS mass of most single-star supernova progenitors from pre-explosion observations. Our results furthermore show that the internal structure of presupernova stars is sensitive to variations in both main sequence and post main-sequence mass loss. The compactness parameter ξ ∝ ℳ /R(ℳ) has been identified as a proxy for the “explodability” of a given presupernova model. We find that ξ varies by as much as 30% for models of the same ZAMS mass evolved with different wind efficiencies and mass loss algorithm combinations. This suggests that the details of the mass loss treatment might bias the outcome of detailed core-collapse supernova calculations and the predictions for neutron star and black hole formation.
We present a state-of-the-art scenario for newly born magnetars as strong sources of gravitational waves (GWs) in the early days after formation. We address several aspects of the astrophysics of ...rapidly rotating, ultra-magnetized neutron stars (NSs), including early cooling before transition to superfluidity, the effects of the magnetic field on the equilibrium shape of NSs, the internal dynamical state of a fully degenerate, oblique rotator and the strength of the electromagnetic torque on the newly born NS. We show that our scenario is consistent with recent studies of supernova remnant surrounding Anomalous X-ray Pulsars (AXPs) and Soft Gamma-Ray Repeaters (SGRs) in the Galaxy that constrains the electromagnetic energy input from the central NS to be ≤ 1051 erg. We further show that if this condition is met, then the GW signal from such sources is potentially detectable with the forthcoming generation of GW detectors up to Virgo cluster distances where an event rate ∼1 yr−1 can be estimated. Finally, we point out that the decay of an internal magnetic field in the 1016 G range couples strongly with the NS cooling at very early stages, thus significantly slowing down both processes: the field can remain this strong for at least 103 yr, during which the core temperature stays higher than several times 108 K.
ABSTRACT We report the Fermi Large Area Telescope (LAT) detections of high-energy (>100 MeV) γ-ray emission from two recent optically bright classical novae, V1369 Centauri 2013 and V5668 Sagittarii ...2015. At early times, Fermi target-of-opportunity observations prompted by their optical discoveries provided enhanced LAT exposure that enabled the detections of γ-ray onsets beginning ∼2 days after their first optical peaks. Significant γ-ray emission was found extending to 39-55 days after their initial LAT detections, with systematically fainter and longer-duration emission compared to previous γ-ray-detected classical novae. These novae were distinguished by multiple bright optical peaks that encompassed the time spans of the observed γ-rays. The γ-ray light curves and spectra of the two novae are presented along with representative hadronic and leptonic models, and comparisons with other novae detected by the LAT are discussed.
Context. The binary star system RS Ophiuchi is a recurrent nova, with outbursts occurring about every 22 years. It consists of a red giant star (RG) and a wind accreting white dwarf close to the ...Chandrasekhar limit. This system is considered a prime candidate for evolving into an SNIa. For its most recent outbursts in 1985 and 2006, exquisite multiwavelength observational data are available. Aims. Deeper physical insight is needed regarding the inter-outburst accretion phase and the dynamical effects of the subsequent nova explosion in order to improve the interpretation of the observed data and to shed light on whether the system is an SNIa progenitor. Methods. We present a 3D hydrodynamic simulation of the quiescent accretion with the subsequent explosive phase. Results. The computed circumstellar mass distribution in the quiescent phase is highly structured with a mass enhancement in the orbital plane of about a factor of 2 as compared to the poleward directions. The simulated nova remnant evolves aspherically, propagating faster toward the poles. The shock velocities derived from the simulations agree with those derived from observations. For $v_{\rm RG} = 20$ km s-1 and for nearly isothermal flows, we find that 10% of the mass lost by the RG is transfered to the WD. For an RG mass loss of $10^{-7}~M_{\odot}$ yr-1, the orbit of the system decays by 3% per million years. With the derived mass transfer rate, multi-cycle nova models provide a qualitatively correct recurrence time, amplitude, and fastness of the nova. Conclusions. Our 3D hydrodynamic simulations provide, along with the observations and nova models, the third ingredient for a deeper understanding of the recurrent novae of the RS Oph type. In combination with recent multi-cycle nova models, our results suggest that the WD in RS Oph will increase in mass. Several speculative outcomes then seem plausible. The WD may reach the Chandrasekhar limit and explode as an SN Ia. Alternatively, the mass loss of the RG could result in a smaller Roch volume, a common envelope phase, and a narrow WD + WD system. Angular momentum loss due to gravitational wave emission could trigger the merger of the two WDs and – perhaps – an SN Ia via the double degenerate scenario.
ABSTRACT
On 2021 August 8, the recurrent nova RS Ophiuchi (RS Oph) erupted again, after an interval of 15.5 yr. Regular monitoring by the Neil Gehrels Swift Observatory began promptly, on August 9.9 ...(0.37 d after the optical peak), and continued until the source passed behind the Sun at the start of November, 86 d later. Observations then restarted on day 197, once RS Oph emerged from the Sun constraint. This makes RS Oph the first Galactic recurrent nova to have been monitored by Swift throughout two eruptions. Here we investigate the extensive X-ray data sets between 2006 and 2021, as well as the more limited data collected by the European X-ray Observatory Satellite (EXOSAT) in 1985. The hard X-rays arising from shock interactions between the nova ejecta and red giant wind are similar following the last two eruptions. In contrast, the early supersoft source (SSS) in 2021 was both less variable and significantly fainter than in 2006. However, 0.3–1 keV light curves from 2021 reveal a 35 s quasi-periodic oscillation consistent in frequency with the 2006 data. The Swift X-ray spectra from 2021 are featureless, with the soft emission typically being well parametrized by a simple blackbody, while the 2006 spectra showed much stronger evidence for superimposed ionized absorption edges. Considering the data after day 60 following each eruption, during the supersoft phase the 2021 spectra are hotter, with smaller effective radii and lower wind absorption, leading to an apparently reduced bolometric luminosity. We explore possible explanations for the gross differences in observed SSS behaviour between the 2006 and 2021 outbursts.
Context. We study the convection zones in the outer envelope of hot massive stars which are caused by opacity peaks associated with iron and helium ionization. Aims. We determine the occurrence and ...properties of these convection zones as function of the stellar parameters. We then confront our results with observations of OB stars. Methods. A stellar evolution code is used to compute a grid of massive star models at different metallicities. In these models, the mixing length theory is used to characterize the envelope convection zones. Results. We find the iron convection zone (FeCZ) to be more prominent for lower surface gravity, higher luminosity and higher initial metallicity. It is absent for luminosities below about $10^{3.2}\,{L}_\odot$, $10^{3.9}\,{L}_\odot$, and $10^{4.2}\,{L}_\odot$ for the Galaxy, LMC and SMC, respectively. We map the strength of the FeCZ on the Hertzsprung-Russell diagram for three metallicities, and compare this with the occurrence of observational phenomena in O stars: microturbulence, non-radial pulsations, wind clumping, and line profile variability. Conclusions. The confirmation of all three trends for the FeCZ as function of stellar parameters by empirical microturbulent velocities argues for a physical connection between sub-photospheric convective motions and small scale stochastic velocities in the photosphere of O- and B-type stars. We further suggest that clumping in the inner parts of the winds of OB stars could be caused by the same mechanism, and that magnetic fields produced in the FeCZ could appear at the surface of OB stars as diagnosed by discrete absorption components in ultraviolet absorption lines.
Context. Classical novae are the product of thermonuclear runaway-initiated explosions occurring on accreting white dwarfs. Aims. V339 Del (Nova Delphinus 2013) was one of the brightest classical ...novae of the last hundred years. Spectroscopy and photometry are available from γ-rays through infrared at stages that have frequently not been observed well. The complete data set is intended to provide a benchmark for comparison with modeling and for understanding more sparsely monitored historical classical and recurrent novae. This paper is the first in the series of reports on the development of the nova. We report here on the early stages of the outburst, through the X-ray active stage. Methods. A time sequence of optical, flux calibrated high resolution spectra was obtained with the Nordic Optical Telescope (NOT) using FIES simultaneously, or contemporaneously, with the Space Telescope Imaging Spectrograph (STIS) aboard the Hubble Space Telescope during the early stages of the outburst. These were supplemented with Mercator/HERMES optical spectra. High resolution IUE ultraviolet spectra of OS And 1986, taken during the Fe curtain phase, served as a template for the distance determination. We used standard plasma diagnostics (e.g., O III and N II line ratios, and the Hβ line flux) to constrain electron densities and temperatures of the ejecta. Using Monte Carlo modeling of the ejecta, we derived the structure, filling factor, and mass from comparisons of the optical and ultraviolet line profiles. Results. We derive an extinction of E(B − V) = 0.23 ± 0.05 from the spectral energy distribution, the interstellar absorption, and H I emission lines. The distance, about 4−4.5 kpc, is in agreement with the inferred distance from near infrared interferometry. The maximum velocity was about 2500 km s-1, measured from the UV resonance and optical profiles. The ejecta showed considerable fine structure in all transitions, much of which persisted as emission knots. The line profiles were modeled using a bipolar conical structure for the ejecta within a relatively restricted range of parameters. For V339 Del, we find that an inclination to the line of sight of about 35°−55°, an opening angle of 60°−80°, and an inner radius ΔR/R(t) ≈ 0.3 based on vrad,max matches the permitted and intercombination lines. The filling factor is f ≈ 0.1, and the derived range in the ejecta mass is (2−3) × 10-5M⊙.
Abstract
We report the Fermi LAT
γ
-ray detection of the 2021 outburst of the symbiotic recurrent nova RS Ophiuchi. In this system, unlike classical novae from cataclysmic binaries, the ejecta from ...the white dwarf form shocks when interacting with the dense circumstellar wind environment of the red giant companion. We find the LAT spectra from 50 MeV to ∼20–23 GeV, the highest-energy photons detected in some subintervals, are consistent with
π
0
-decay emission from shocks in the ejecta as proposed by Tatischeff & Hernanz for its previous 2006 outburst. The LAT light curve displayed a fast rise to its peak >0.1 GeV flux of ≃6 × 10
−6
ph cm
−2
s
−1
beginning on day 0.745 after its optically constrained eruption epoch of 2021 August 8.50. The peak lasted for ∼1 day and exhibited a power-law decline up to the final LAT detection on day 45. We analyze the data on shorter timescales at early times and found evidence of an approximate doubling of emission over ∼200 minutes at day 2.2, possibly indicating a localized shock-acceleration event. Comparing the data collected by the American Association of Variable Star Observers, we measured a constant ratio of ∼ 2.8 × 10
−3
between the
γ
-ray and optical luminosities except for a ∼5×smaller ratio within the first day of the eruption likely indicating attenuation of
γ
rays by ejecta material and lower high-energy proton fluxes at the earliest stages of the shock development. The hard X-ray emission due to bremsstrahlung from shock-heated gas traced by the Swift-XRT 2–10 keV light curve peaked at day ∼6, later than at GeV and optical energies. Using X-ray derived temperatures to constrain the velocity profile, we find the hadronic model reproduces the observed >0.1 GeV light curve.
We present 5-28 m SOFIA FORECAST spectroscopy complemented by panchromatic X-ray through infrared observations of the CO nova V5668 Sgr documenting the formation and destruction of dust during ∼500 ...days following outburst. Dust condensation commenced by 82 days after outburst at a temperature of ∼1090 K. The condensation temperature indicates that the condensate was amorphous carbon. There was a gradual decrease of the grain size and dust mass during the recovery phase. Absolute parameter values given here are for an assumed distance of 1.2 kpc. We conclude that the maximum mass of dust produced was 1.2 × 10−7 M if the dust was amorphous carbon. The average grain radius grew to a maximum of ∼2.9 m at a temperature of ∼720 K around day 113 when the shell visual optical depth was τv ∼ 5.4. Maximum grain growth was followed by a period of grain destruction. X-rays were detected with Swift from day 95 to beyond day 500. The Swift X-ray count rate due to the hot white dwarf peaked around day 220, when its spectrum was that of a kT = 35 eV blackbody. The temperature, together with the supersoft X-ray turn-on and turn-off times, suggests a white dwarf mass of ∼1.1 M . We show that the X-ray fluence was sufficient to destroy the dust. Our data show that the post-dust event X-ray brightening is not due to dust destruction, which certainly occurred, as the dust is optically thin to X-rays.
Aims.This paper develops a method for obtaining the star formation histories of a mixed, resolved population through the use of color-magnitude diagrams (CMDs). The method provides insight into the ...local star formation rate, analyzing the observations of the Hipparcos satellite through a comparison with synthetic CMDs computed for different histories with an updated stellar evolution library. Methods.Parallax and photometric uncertainties are included explicitly and corrected using the Bayesian Richardson-Lucy algorithm. We first describe our verification studies using artificial data sets. From this sensitivity study, the critical factors determining the success of a recovery for a known star formation rate are a partial knowledge of the IMF and the age-metallicity relation, and sample contamination by clusters and moving groups (special populations whose histories are different than that of the whole sample). Unresolved binaries are less important impediments. We highlight how these limit the method. Results.For the real field sample, complete to $M_V < 3.5$, we find that the solar neighborhood star formation rate has a characteristic timescale for variation of about 6 Gyr, with a maximum activity close to 3 Gyr ago. The similarity of this finding with column integrated star formation rates may indicate a global origin, possibly a collision with a satellite galaxy. We also discuss applications of this technique to general photometric surveys of other complex systems (e.g. Local Group dwarf galaxies) where the distances are well known.