No Arabic abstract
LBVs are massive evolved stars that suffer sporadic and violent mass-loss events. They have been proposed as the progenitors of some core-collapse SNe, but this idea is still debated due to the lack of direct evidence. Since SNRs can carry in their morphology the fingerprints of the progenitor stars as well as of the inhomogeneous CSM sculpted by the progenitors, the study of SNRs from LBVs could help to place core-collapse SNe in context with the evolution of massive stars. We investigate the physical, chemical and morphological properties of the remnants of SNe originating from LBVs, in order to search for signatures, revealing the nature of the progenitors, in the ejecta distribution and morphology of the remnants. As a template of LBVs, we considered the actual LBV candidate Gal 026.47+0.02. We selected a grid of models, which describe the evolution of a massive star with properties consistent with those of Gal 026.47+0.02 and its final fate as core-collapse SN. We developed a 3D HD model that follows the post-explosion evolution of the ejecta from the breakout of the shock wave at the stellar surface to the interaction of the SNR with a CSM characterized by two dense nested toroidal shells, parametrized in agreement with multi-wavelength observations of Gal 026.47+0.02. Our models show a strong interaction of the blast wave with the CSM which determines an important slowdown of the expansion of the ejecta in the equatorial plane where the two shells lay, determining a high degree of asymmetry in the remnant. After 10000 years of evolution the ejecta show an elongated shape forming a broad jet-like structure caused by the interaction with the shells and oriented along the axis of the toroidal shells.
We summarize our current understanding of gravitational wave emission from core-collapse supernovae. We review the established results from multi-dimensional simulations and, wherever possible, provide back-of-the-envelope calculations to highlight the underlying physical principles. The gravitational waves are predominantly emitted by protoneutron star oscillations. In slowly rotating cases, which represent the most common type of the supernovae, the oscillations are excited by multi-dimensional hydrodynamic instabilities, while in rare rapidly rotating cases, the protoneutron star is born with an oblate deformation due to the centrifugal force. The gravitational wave signal may be marginally visible with current detectors for a source within our galaxy, while future third-generation instruments will enable more robust and detailed observations. The rapidly rotating models that develop non-axisymmetric instabilities may be visible up to a megaparsec distance with the third-generation detectors. Finally, we discuss strategies for multi-messenger observations of supernovae.
Context. Luminous Blue Variables (LBVs) are thought to be in a transitory phase between O stars on the main-sequence and the Wolf-Rayet stage. Recent studies suggest that they might be formed through binary interaction. Only a few are known in binary systems but their multiplicity fraction is uncertain. Aims. This study aims at deriving the binary fraction among the Galactic (confirmed and candidate) LBV population. We combine multi-epoch spectroscopy and long-baseline interferometry. Methods. We use cross-correlation to measure their radial velocities. We identify spectroscopic binaries through significant RV variability (larger than 35 km/s). We investigate the observational biases to establish the intrinsic binary fraction. We use CANDID to detect interferometric companions, derive their parameters and positions. Results. We derive an observed spectroscopic binary fraction of 26 %. Considering period and mass ratio ranges from Porb=1 to 1000 days, and q = 0.1-1.0, and a representative set of orbital parameter distributions, we find a bias-corrected binary fraction of 62%. From interferometry, we detect 14 companions out of 18 objects, providing a binary fraction of 78% at projected separations between 1 and 120 mas. From the derived primary diameters, and the distances of these objects, we measure for the first time the exact radii of Galactic LBVs to be between 100 and 650 Rsun, making unlikely to have short-period systems. Conclusions. This analysis shows that the binary fraction among the Galactic LBV population is large. If they form through single-star evolution, their orbit must be initially large. If they form through binary channel that implies that either massive stars in short binary systems must undergo a phase of fully non-conservative mass transfer to be able to sufficiently widen the orbit or that LBVs form through merging in initially binary or triple systems.
Unraveling the mechanism for core-collapse supernova explosions is an outstanding computational challenge and the problem remains essentially unsolved despite more than four decades of effort. However, much progress in realistic modeling has occurred recently through the availability of multi-teraflop machines and the increasing sophistication of supernova codes. These improvements have led to some key insights which may clarify the picture in the not too distant future. Here we briefly review the current status of the three explosion mechanisms (acoustic, MHD, and neutrino heating) that are currently under active investigation, concentrating on the neutrino heating mechanism as the one most likely responsible for producing explosions from progenitors in the mass range ~10 to ~25 solar masses. We then briefly describe the CHIMERA code, a supernova code we have developed to simulate core-collapse supernovae in 1, 2, and 3 spatial dimensions. We finally describe the results of an ongoing suite of 2D simulations initiated from a 12, 15, 20, and 25 solar mass progenitor. These have all exhibited explosions and are currently in the expanding phase with the shock at between 5,000 and 10,000 km. We finally very briefly describe an ongoing simulation in 3 spatial dimensions initiated from the 15 solar mass progenitor.
The Luminous Blue Variable stars exhibit behavior ranging from light curve `microvariations on timescales of tens of days, to `outbursts accompanied by mass loss of up to 10e-03 solar masses per year, occurring decades apart, to `giant eruptions such as seen in Eta Carinae ejecting one or more solar masses and recurring on timescales of centuries. Here we review the work of the Los Alamos group since 1993 to investigate pulsations and instabilities in massive stars using linear pulsation models and non-linear hydrodynamic models. The models predict pulsational variability that may be associated with the microvariations. Using a nonlinear pulsation hydrodynamics code with a time-dependent convection treatment, we show that, in some circumstances, the Eddington limit is exceeded periodically in the pulsation driving region of the stellar envelope, accelerating the outer layers, and perhaps initiating mass loss or LBV outbursts. We discuss how pulsations and mass loss may be responsible for the location of the Humphreys-Davidson Limit in the H-R diagram. The `giant eruptions, however, must involve much deeper regions in the stellar core to cause such large amounts of mass to be ejected. We review and suggest some possible explanations, including mixing from gravity modes, secular instabilities, the epsilon mechanism, or the SASI instability as proposed for Type II supernovae. We outline future work and required stellar modeling capabilities to investigate these possibilities.
Nuclear shell burning in the final stages of the lives of massive stars is accompanied by strong turbulent convection. The resulting fluctuations aid supernova explosion by amplifying the non-radial flow in the post-shock region. In this work, we investigate the physical mechanism behind this amplification using a linear perturbation theory. We model the shock wave as a one-dimensional planar discontinuity and consider its interaction with vorticity and entropy perturbations in the upstream flow. We find that, as the perturbations cross the shock, their total turbulent kinetic energy is amplified by a factor of $sim!2$, while the average linear size of turbulent eddies decreases by about the same factor. These values are not sensitive to the parameters of the upstream turbulence and the nuclear dissociation efficiency at the shock. Finally, we discuss the implication of our results for the supernova explosion mechanism. We show that the upstream perturbations can decrease the critical neutrino luminosity for producing explosion by several percent.