Classical Wolf-Rayet (WR) stars are massive, hydrogen depleted, post main-sequence stars that exhibit emission-line dominated spectra. For a given metallicity Z, stars exceeding a certain initial mass M_single(Z) can reach the WR phase through intrinsic mass-loss (single-star channel). Stars of lower masses can reach the WR phase via binary mass transfer (binary channel). It is commonly assumed that the binary channel dominates the formation of WR stars in environments with low Z such as the SMC and LMC. However, their reported WR binary fractions of 30-40% are comparable to that of the Galaxy. Here, we explain this apparent contradiction by considering the minimum initial mass M_spec(Z) needed for the stripped product to appear as a WR star. We calibrate M_spec(Z) using the lowest-luminosity WR stars in the Clouds and the Galaxy. A range of M_single(Z) values are explored using various evolution codes. We estimate the additional contribution of the binary channel by considering the interval [M_spec(Z), M_single(Z)], which characterises the initial-mass range in which binaries can form additional WR stars. Results: The WR-phenomenon ceases below luminosities of logL = 4.9, 5.25, and 5.6 [Lsun] in the Galaxy, LMC, and SMC, which translates to He-star masses of 7.5, 11, 17 Msun and initial masses of M_spec = 18, 23, 37 Msun. Stripped stars with lower initial masses in the respective galaxies would tend to not appear as WR stars. M_single lies in the range 20-30, 30-60, and > 40 Msun for the Galaxy, LMC, and SMC. We find that that the additional contribution of the binary channel is a non-trivial function of Z that cannot be conclusively claimed to be monotonically increasing with decreasing Z. Hence, one should not a-priori expect that binary interactions become increasingly important in forming WR stars at low Z, or that the WR binary fraction grows with decreasing Z.

We simulate the evolution of binary systems with a massive primary star of 15Mo where we introduce an enhanced mass loss due to jets that the secondary star might launch, and find that in many cases the enhanced mass loss brings the binary system to experience the grazing envelope evolution (GEE) and form a progenitor of Type IIb supernova (SN IIb). The jets, the Roche lobe overflow (RLOF), and a final stellar wind remove most of the hydrogen-rich envelope, leaving a blue-compact SN IIb progenitor. In many cases without this jet-driven mass loss the system enters a common envelope evolution (CEE) and does not form a SN IIb progenitor. We use the stellar evolutionary code MESA binary and mimic the jet-driven mass loss with a simple prescription and some free parameters. Our results show that the jet-driven mass loss, that some systems have during the GEE, increases the parameter space for stellar binary systems to form SN IIb progenitors. We estimate that the binary evolution channel with GEE contributes about a quarter of all SNe IIb, about equal to the contribution of each of the other three channels, binary evolution without a GEE, fatal CEE (where the secondary star merges with the core of the giant primary star), and the single star channel.

From stellar evolution simulations (using MESA) we conclude that the fatal common envelope evolution (CEE) channel for the formation of Type IIb core collapse supernova (SN IIb) progenitors can indeed account for some SNe IIb. In the fatal CEE channel for SNe IIb a low mass main sequence secondary star inspirals inside the giant envelope of the massive primary star and removes most of the giant envelope before it merges with the giant core. The key ingredient of the scenario studied here is that the tidally destroyed secondary star forms a new giant envelope. The mass-loss process in a wind during the evolution from the merger process until core collapse, i.e., until the explosion, leaves little hydrogen mass at explosion as inferred from observations of SNe IIb. In the case of a massive primary star with a zero age main sequence mass of M(ZAMS)=16Mo that during its giant phase swallows a main sequence star of mass M2=0.5Mo, we find at explosion a hydrogen mass of M(H)=0.02-0.09Mo, depending on the rotation we assume. We find similar values for M(ZAMS)=12Mo.

We find that applying a theoretical wind mass-loss rate from Monte Carlo radiative transfer models for hydrogen-deficient stars results in significantly more leftover hydrogen following stable mass transfer through Roche-lobe overflow than when we use an extrapolation of an empirical fit for Galactic Wolf-Rayet stars, for which a negligible amount of hydrogen remains in a large set of binary stellar evolution computations. These findings have implications for modelling progenitors of Type Ib and Type IIb supernovae. Most importantly, our study stresses the sensitivity of the stellar evolution models to the assumed mass-loss rates and the need to develop a better theoretical understanding of stellar winds.

We propose a common-envelope jets supernova (CEJSN) scenario for the fast-rising blue optical transient AT2018cow. In a CEJSN a neutron star (NS) spirals-in inside the extended envelope of a massive giant star and enters the core. The NS accretes mass from the core through an accretion disc and launches jets. These jets explode the core and the envelope. In the specific polar CEJSN scenario that we propose here the jets clear the polar regions of the giant star before the NS enters the core. The jets that the NS launches after it enters the core expand almost freely along the polar directions that contain a small amount of mass. This, we suggest, explains the fast rise to maximum and the fast ejecta observed at early times of the enigmatic transient AT2018cow. The slower later time ejecta is the more massive equatorial outflow. We roughly estimate the accretion phase onto the NS during the explosion phase to last for a time of about 1000 seconds, during which the average mass accretion rate is about 0.0001 Mo/sec. We outline the possible diversity of CEJSNe by listing five other scenarios in addition to the polar CEJSN scenario.

We propose a new type of repeating transient outburst initiated by a neutron star (NS) entering the envelope of an evolved massive star, accreting envelope material and subsequently launching jets which interact with their surroundings. This interaction is the result of either a rapid expansion of the primary star due to an instability in its core near the end of its nuclear evolution, or due to a dynamical process which rapidly brings the NS into the primary star. The ejecta can reach velocities of 10,000 km/s despite not being a supernova, and might explain such velocities in the 2011 outburst of the luminous blue variable progenitor of SN 2009ip. The typical transient duration and kinetic energy are weeks to months, and up to approximately 1 foe, respectively. The interaction of a NS with a giant envelope might be a phase in the evolution of the progenitors of most NS-NS binary systems that later undergo a merger event. If the NS spirals in all the way to the core of the primary star and brings about its complete disruption we term this a `common envelope jets supernova' (CEJSN), which is a possible explanation for the peculiar supernova iPTF14hls. For a limited interaction of the NS with the envelope we get a less luminous transient, which we term a CEJSN impostor.

We use the binary module of the MESA code to study the evolution of an evolved binary system where we assume that a main sequence companion removes the outskirts of the envelope of an asymptotic giant branch (AGB) star by launching jets, and explore the characteristics of this grazing envelope evolution (GEE). We base our assumption that jets launched by the secondary star remove a substantial fraction of the outskirts of the envelope of an AGB star on earlier hydrodynamical simulations. We find that in many cases that we study, but not in all cases, the binary system experiences the GEE rather than entering the common envelope phase, under our assumptions of jet-driven mass removal. To prevent the common envelope phase, we assume the secondary star may accrete a large amount of mass in a short time while avoiding rapid inflation, the feasibility of which requires further study. Because of our simplifying assumptions we cannot yet present the parameter space for the GEE. Although the incorporation of the GEE into population synthesis numerical codes requires further studies of the GEE, we conclude that analyses of population synthesis studies of evolved binary stars should include the GEE.

We evolve models of rotating massive stars up to the stage of iron core collapse using the MESA code and find a shell with a mixed composition of primarily helium and oxygen in some cases. In the parameter space of initial masses of 13-40 Mo and initial rotation velocities of 0-450 km/s that we investigate, we find a mixed helium-oxygen shell with a significant mass only for a small fraction of the models. While the shell is formed due to rotational mixing, the pre-collapse rotation rate is not very high, as required for an energetic collapse-induced thermonuclear explosion. Our results suggest that the collapse-induced thermonuclear explosion mechanism that was revisited recently can account for at most a small fraction of core-collapse supernovae. The presence of such a mixed helium-oxygen shell still might have some implications for core-collapse supernovae, such as some nucleosynthesis processes when jets are present, or might result in peculiar sub-luminous core-collapse supernovae.

We propose a common envelope jets supernova scenario for the enigmatic supernova iPTF14hls where a neutron star that spirals-in inside the envelope of a massive giant star accretes mass and launches jets that power the ejection of the circumstellar shell and a few weeks later the explosion itself. To account for the kinetic energy of the circumstellar gas and the explosion, the neutron star should accrete a mass of approximately 0.3 solar masses. The tens of solar masses of circumstellar gas that accounts for some absorption lines is ejected while the neutron star orbits for about one to several weeks inside the envelope of the giant star. In the last hours of the interaction the neutron star merges with the core, accretes mass, and launches jets that eject the core and the inner envelope to form the explosion itself and the medium where the supernova photosphere resides. The remaining neutron star accretes fallback gas and further powers the supernova. We attribute the 1954 pre-explosion outburst to an eccentric orbit and temporary mass accretion by the neutron star at periastron passage prior to the onset of the common envelope phase.

Non-axisymmetric features are found in the core collapse of a rapidly rotating massive star, which might have important implications for magnetic field amplification and production of a bipolar outflow that can explode the star, as well as for r-process nucleosynthesis and natal kicks. The collapse of an evolved rapidly rotating massive star is followed in three-dimensional hydrodynamic simulations using the FLASH code with neutrino leakage. A rotating proto-neutron star (PNS) forms with a non-zero linear velocity. This can contribute to the natal kick of the remnant compact object. The PNS is surrounded by a turbulent medium, where high shearing is likely to amplify magnetic fields, which in turn can drive a bipolar outflow. Neutron-rich material in the PNS vicinity might induce strong r-process nucleosynthesis. The rapidly rotating PNS possesses a rotational energy of E>10foe. Magnetar formation proceeding in a similar fashion will be able to deposit a portion of this energy later on in the supernova ejecta through a spin down mechanism. These processes can be important for rare supernovae generated by rapidly rotating progenitors, even though a complete explosion is not simulated in the present study.

We analyze recent magnetar light-curve modeling of 38 hydrogen-poor superluminous supernovae (SLSNe), and find that the energies of the explosions themselves, that take place before the magnetar energy is released, are more than what the neutrino-driven explosion mechanism can supply for about half of the systems. These SLSNe must have been exploded by a different process than the delayed neutrino mechanism, most likely the jet feedback mechanism (JFM). The conclusion for magnetar modeling of SLSNe is that jets launched at magnetar birth cannot be ignored, not at the explosion itself and not later when mass fall-back might occur. More generally, the present analysis strengthens the call for a paradigm shift from neutrino-driven to jet-driven explosion models of all core collapse supernovae.

We evolve stellar models to study the rotational profiles of the pre-explosion cores of single massive stars that are progenitors of core collapse supernovae (CCSNe), and find large rotational shear above the iron core that might play an important role in the jet feedback explosion mechanism by amplifying magnetic fields before and after collapse. Initial masses of 15 Mo and 30 Mo and various values of the initial rotation velocity are considered, as well as a reduced mass-loss rate along the evolution and the effect of core-envelope coupling through magnetic fields. We find that the rotation profiles just before core collapse differ between models, but share the following properties. (1) There are narrow zones of very large rotational shear adjacent to convective zones. (2) The rotation rate of the inner core is slower than required to form a Keplerian accretion disk. (3) The outer part of the core and the envelope have non-negligible specific angular momentum compared to the last stable orbit around a black hole (BH). Our results suggest the feasibility of magnetic field amplification which might aid a jet-driven explosion leaving behind a neutron star. Alternatively, if the inner core fails in exploding the star, an accretion disk from the outer parts of the core might form and lead to a jet-driven CCSN which leaves behind a BH.

We propose a speculative scenario where dynamo amplification of magnetic fields in the core convective shells of massive stars, tens of years to hours before they explode, leads to envelope expansion and enhanced mass loss rate, resulting in pre-explosion outbursts (PEOs). The convective luminosity in the burning shells of carbon, neon, oxygen, and then silicon, are very high. Based on the behavior of active main sequence stars we speculate that the convective shells can trigger magnetic activity with a power of about 0.001 times the convective luminosity. Magnetic flux tubes might buoy outward, and deposit their energy in the outer parts of the envelope. This in turn might lead to the expansion of the envelope and to an enhanced mass loss rate. If a close binary companion is present, mass transfer might take place and lead to an energetic outburst. The magnetic activity requires minimum core rotation and that the stochastic magnetic activity be on its high phase. Only in rare cases these conditions are met, accounting for that only the minority of core collapse supernovae (CCSNe) experience PEO. Such a pre-explosion magnetic activity might have implications for the explosion mechanism itself.

We find significant fluctuations of angular momentum within the convective helium shell of a pre-collapse massive star - a core-collapse supernova progenitor - which may facilitate the formation of accretion disks and jets that can explode the star. The convective flow in our model of an evolved M_ZAMS=15Msun star, computed with the sub-sonic hydrodynamic solver MAESTRO, contains entire shells with net angular momentum in different directions. This phenomenon may have important implications for the late evolutionary stages of massive stars, and for the dynamics of core-collapse.

We suggest that the energetic radiation from core-collapse super-energetic supernovae (SESNe) is due to a long lasting accretion process onto the newly born neutron star (NS), resulting from an inefficient operation of the jet-feedback mechanism. The jets that are launched by the accreting NS or black hole (BH) maintain their axis due to a rapidly rotating pre-collapse core, and do not manage to eject core material from near the equatorial plane. The jets are able to eject material from the core along the polar directions, and reduce the gravity near the equatorial plane. The equatorial gas expands, and part of it falls back over a timescale of minutes to days to prolong the jets-launching episode. According to the model for SESNe proposed in the present paper, the principal parameter that distinguishes between the different cases of CCSN explosions, such as between normal CCSNe and SESNe, is the efficiency of the jet-feedback mechanism. This efficiency in turn depends on the pre-collapse core mass, envelope mass, core convection, and most of all on the angular momentum profile in the core. One prediction of the inefficient jet-feedback mechanism for SESNe is the formation of a slow equatorial outflow in the explosion. Typical velocity and mass of this outflow are estimated to be approximately 1000 km/s and greater than about 1 solar mass, respectively, though quantitative values will have to be checked in future hydrodynamic simulations.

We show that the standing accretion shock instability (SASI) that has been used to ease the shock revival in core collapse supernovae (CCSNe) neutrino-driven explosion models, might play a much more decisive role in supplying the stochastic angular momentum required to trigger an explosion with jittering jets. We find that if the kinetic energy associated with the transverse (non radial) motion of the SASI is larger than about ten percent of the energy associated with the energy of the accreted gas, then the stochastic angular momentum can reach about five percent of the Keplerian specific angular momentum around the newly born neutron star. Such an accretion flow leaves an open conical region along the poles with an average opening angle of about 5 degrees. The outflow from the open polar regions powers an explosion according to the jittering-jets model.

We show that turbulence in core collapse supernovae (CCSNe) which has been shown recently to ease shock revival might also lead to the formation of intermittent thick accretion disks, or accretion belts, around the newly born neutron star (NS). The accretion morphology is such that two low density funnels are formed along the polar directions. The disks then are likely to launch jets with a varying axis direction, i.e., jittering-jets, through the two opposite funnels. The energy contribution of jets in this jittering jets mechanism might result in an explosion energy of E>10^51erg, even without reviving the stalled shock. We strengthen the jittering jets mechanism as a possible explosion mechanism of CCSNe.

We find that convective regions of collapsing massive stellar cores possess sufficient stochastic angular momentum to form intermittent accretion disks around the newly born neutron star (NS) or black hole (BH), as required by the jittering-jets model for core-collapse supernova (CCSN) explosions. To reach this conclusion we derive an approximate expression for stochastic specific angular momentum in convection layers of stars, and using the mixing-length theory apply it to four stellar models at core-collapse epoch. In all models, evolved using the stellar evolution code MESA, the convective helium layer has sufficient angular momentum to form an accretion disk. The mass available for disk formation around the NS or BH is 0.1-1.2Mo; stochastic accretion of this mass can form intermittent accretion disks that launch jets powerful enough to explode the star according to the jittering-jets model. Our results imply that even if no explosion occurs after accretion of the inner ~2-5Mo of the core onto the NS or BH (the mass depends on the stellar model), accretion of outer layers of the core will eventually lead to an energetic supernova explosion.

By simulating jet-inflated bubbles in cooling flows with the PLUTO hydrodynamic code we show that mixing of high entropy shocked jet's material with the intra-cluster medium (ICM) is the major heating process perpendicular to the jets' axis. Heating by the forward shock is not significant. The mixing is very efficient in heating the ICM in all directions, to distances of ~10kpc and more. Although the jets are active for a time period of only 20 Myr, the mixing and heating near the equatorial plane, as well as along the symmetry axis, continues to counter radiative cooling for times of >10^8 yr after the jets have ceased to exist. We discuss some possible implications of the results. (i) The vigorous mixing is expected to entangle magnetic field lines, hence to suppress any global heat conduction in the ICM near the center. (ii) The vigorous mixing forms multi-phase ICM in the inner cluster regions, where the coolest parcels of gas will eventually cool first, flow inward, and feed the active galactic nucleus to set the next jet-activity episode. This further supports the cold feedback mechanism. (iii) In cases where the medium outside the region of r~10kpc is not as dense as in groups and clusters of galaxies, like during the process of galaxy formation, the forward shock and the high pressure of the shocked jets' material might expel gas from the system.

We study the similarities of jet-medium interactions in several quite different astrophysical systems using 2D and 3D hydrodynamical numerical simulations, and find many similarities. The systems include cooling flow (CF) clusters of galaxies, core-collapse supernovae (CCSNe), planetary nebulae (PNe), and common envelope (CE) evolution. The similarities include hot bubbles inflated by jets in a bipolar structure, vortices on the sides of the jets, vortices inside the inflated bubbles, fragmentation of bubbles to two and more bubbles, and buoyancy of bubbles. The activity in many cases is regulated by a negative feedback mechanism. Namely, higher accretion rate leads to stronger jet activity that in turn suppresses the accretion process. After the jets power decreases the accretion resumes, and the cycle restarts. In the case of CF in galaxies and clusters of galaxies we also study the accretion process, which is most likely by cold clumps, i.e., the cold feedback mechanism. In CF clusters we find that heating of the intra-cluster medium (ICM) is done by mixing hot shocked jet gas with the ICM, and not by shocks. Our results strengthen the jet feedback mechanism (JFM) as a common process in many astrophysical objects.