This paper investigates the recent stellar merger transient M31LRN 2015 in the Andromeda galaxy. We analyze published optical photometry and spectroscopy along with a Hubble Space Telescope detection of the color and magnitude of the pre-outburst source. The transient outburst is consistent with dynamically driven ejecta at the onset of a common envelope episode, which eventually leads to the complete merger of a binary system. The light curve appears to contain two components: first $\sim10^{-2} M_\odot$ of fast ejecta driven by shocks at the onset of common envelope, and later, $\sim0.3 M_\odot$ of further ejecta as the secondary becomes more deeply engulfed within the primary. Just prior to merger, we find that the primary star is a $3-5.5 M_\odot$ sub-giant branch primary star with radius of $30-40R_\odot$. Its position in the color-magnitude diagram shows that it is growing in radius, consistent with a picture where it engulfs its companion. By matching the properties of the primary star to the transient outburst, we show that the optical transient lasts less than ten orbits of the original binary, which had a pre-merger period of $\sim 10$ days. We consider the possible orbital dynamics leading up to the merger, and argue that if the system merged due to the Darwin tidal instability it implies a lower mass main sequence companion of $0.1-0.6M_\odot$. This analysis represents a promising step toward a more detailed understanding of flows in common envelope episodes through direct observational constraints.

The most popular model for short gamma-ray bursts (sGRBs) involves the coalescence of binary neutron stars. Because the progenitor is actually hidden from view, we must consider under which circumstances such merging systems are capable of producing a successful sGRB. Soon after coalescence, winds are launched from the merger remnant. In this paper, we use realistic wind profiles derived from global merger simulations in order to investigate the interaction of sGRB jets with these winds using numerical simulations. We analyze the conditions for which these axisymmetric winds permit relativistic jets to breakout and produce a sGRB. We find that jets with luminosities comparable to those observed in sGRBs are only successful when their half-opening angles are below ~20{\deg}. This jet collimation mechanism leads to a simple physical interpretation of the luminosities and opening angles inferred for sGRBs. If wide, low luminosity jets are observed, they might be indicative of a different progenitor avenue such as the merger of a neutron star with a black hole. We also use the observed durations of sGRB to place constraints on the lifetime of the wind phase, which is determined by the time it takes the jet to breakout. In all cases we find that the derived limits argue against completely stable remnants for binary neutron star mergers that produce sGRBs.

The r-process nuclei are robustly synthesized in the material ejected during a neutron star binary merger (NSBM), as tidal torques transport angular momentum and energy through the outer Lagrange point in the form of a vast tidal tail. If NSBM are indeed solely responsible for the solar system r- process abundances, a galaxy like our own would require to host a few NSBM per million years, with each event ejecting, on average, about 5x10^{-2} M_sun of r-process material. Because the ejecta velocities in the tidal tail are significantly larger than in ordinary supernovae, NSBM deposit a comparable amount of energy into the interstellar medium (ISM). In contrast to extensive efforts studying spherical models for supernova remnant evolution, calculations quantifying the impact of NSBM ejecta in the ISM have been lacking. To better understand their evolution in a cosmological context, we perform a suite of three-dimensional hydrodynamic simulations with optically-thin radiative cooling of isolated NSBM ejecta expanding in environments with conditions adopted from Milky Way-like galaxy simulations. Although the remnant morphology is highly complex at early times, the subsequent radiative evolution that results from thermal instability in atomic gas is remarkably similar to that of a standard supernova blast wave. This implies that sub-resolution supernova feedback models can be used in galaxy-scale simulations that are unable to resolve the key evolutionary phases of NSBM blast waves. Among other quantities, we examine the radius, time, mass and kinetic energy content of the NSBM remnant at shell formation as well as the momentum injected to the ISM. We find that the shell formation epoch is attained when the swept-up mass is about 10^3 M_sun, at this point the mass fraction of r-process material is drastically enhanced up to two orders of magnitude in relation to a solar metallicity ISM.

Short gamma-ray bursts (sGRBs) are widely believed to result from the mergers of compact binaries. This model predicts an afterglow that bears the characteristic signatures of a constant, low density medium, including a smooth prompt-afterglow transition, and a simple temporal evolution. However, these expectations are in conflict with observations for a non-negligible fraction of sGRB afterglows. In particular, the onset of the afterglow phase for some of these events appears to be delayed and, in addition, a few of them exhibit late- time rapid fading in their lightcurves. We show that these peculiar observations can be explained independently of ongoing central engine activity if some sGRB progenitors are compact binaries hosting at least one pulsar. The Poynting flux emanating from the pulsar companion can excavate a bow-shock cavity surround- ing the binary. If this cavity is larger than the shock deceleration length scale in the undisturbed interstellar medium, then the onset of the afterglow will be delayed. Should the deceleration occur entirely within the swept-up thin shell, a rapid fade in the lightcurve will ensue. We identify two types of pulsar that can achieve the conditions necessary for altering the afterglow: low field, long lived pulsars, and high field pulsars. We find that a sizable fraction (~20-50%) of low field pulsars are likely to reside in neutron star binaries based on observations, while their high field counterparts are not. Hydrodynamical calculations motivated by this model are shown to be in good agreement with observations of sGRB afterglow lightcurves.

M17 JVLA 35 is a radio source detected in projection against the M17 HII region. In recent observations, its spectrum between 4.96 and 8.46 GHz was found to be positive and very steep, with $\alpha \geq 2.9 \pm 0.6$ ($S_\nu \propto \nu^\alpha$). Here we present Very Large Array observations made in the 18.5 to 36.5 GHz region that indicate a spectral turnover at $\sim$13 GHz and a negative spectral index ($\alpha \simeq -2.0$) at higher frequencies. The spectrum is consistent with that of an extragalactic High Frequency Peaker (HFP). However, M17 JVLA 35 has an angular size of $\sim0\rlap.{''}5$ at 8.46 GHz, while HFPs have extremely compact, milliarcsecond dimensions. We discuss other possible models for the spectrum of the source and do not find them feasible. Finally, we propose that M17 JVLA35 is indeed an HFP but that its angular size becomes broadened by plasma scattering as its radiation travels across M17. If our interpretation is correct, accurate measurements of the angular size of M17 JVLA35 across the centimeter range should reveal the expected $\nu^{-2}$ dependence.

The central engine of short gamma-ray bursts (sGRBs) is hidden from direct view, operating at a scale much smaller than that probed by the emitted radiation. Thus we must infer its origin not only with respect to the formation of the {\it trigger} - the actual astrophysical configuration that is capable of powering a sGRB - but also from the consequences that follow from the various evolutionary pathways that may be involved in producing it. Considering binary neutron star mergers we critically evaluate, analytically and through numerical simulations, whether the neutrino-driven wind produced by the newly formed hyper-massive neutron star can allow the collimated relativistic outflow that follows its collapse to actually produce a sGRB or not. Upon comparison with the observed sGRB duration distribution, we find that collapse cannot be significantly delayed (<=100 ms) before the outflow is choked, thus limiting the possibility that long-lived hyper-massive remnants can account for these events. In the case of successful breakthrough of the jet through the neutrino-driven wind, the energy stored in the cocoon could contribute to the precursor and extended emission observed in sGRBs.

The problem of explaining the X-ray emission properties of the massive, close binary WR 20a is discussed. Located near the cluster core of Westerlund 2, WR 20a is composed of two nearly identical Wolf- Rayet stars of 82 and 83 solar masses orbiting with a period of only 3.7 days. Although Chandra observations were taken during the secondary optical eclipse, the X-ray light curve shows no signs of a flux decrement. In fact, WR 20a appears slightly more X-ray luminous and softer during the optical eclipse, opposite to what has been observed in other binary systems. To aid in our interpretation of the data, we compare with the results of hydrodynamical simulations using the adaptive mesh refinement code Mezcal that includes radiative cooling and a radiative acceleration force term. It is shown that the X-ray emission can be successfully explained in models where the wind-wind collision interface in this system occurs while the outflowing material is still being accelerated. Consequently, WR 20a serves as a critical test-case for how radiatively-driven stellar winds initiate and interact. Our models not only procure a robust description of current Chandra data, which cover the orbital phases between 0.3 to 0.6, but provide detailed predictions over the entire orbit.

The classic HII region M17 is one of the best studied across the electromagnetic spectrum. We present sensitive, high angular resolution observations made with the Jansky Very Large Array (JVLA) at 4.96 and 8.46 GHz that reveal the presence of 38 compact radio sources, in addition to the well known hypercompact cometary HII region M17 UC1. For this last source we find that its spectral index of value $\sim$1 is due to a gradient in opacity across its face. Of the 38 compact radio sources detected, 19 have stellar counterparts detected in the infrared, optical, or X-rays. Finally, we discuss the nature of the radio emission from the massive binary system CEN 1a and 1b, concluding that both are most probably non-thermal emitters, although the first is strongly time variable and the second is steady.

Aims. The interaction of two isotropic stellar winds is studied in order to calculate the free-free emission from the wind collision region. The effects of the binary separation and the wind momentum ratio on the emission from the wind-wind interaction region are investigated. Methods. We developed a semi-analytical model for calculating the thermal emission from colliding stellar winds. Assuming radiative shocks for the compressed layer, which are expected in close binaries, we obtained the emission measure of the thin shell. Then, we computed the total optical depth along each line of sight to obtain the emission from the whole configuration. Results. Here, we present predictions of the free-free emission at radio frequencies from analytic, radiative shock models in colliding wind binaries. It is shown that the emission from the wind collision region mainly arises from the optically thick region of the compressed layer and scales as ~ D^{4/5}, where D is the binary separation. The predicted flux density from the wind collision region becomes more important as the frequency increases, showing higher spectral indices than the expected 0.6 value from the unshocked winds. We also investigate the emission from short-period WR+O systems calculated with our analytic formulation. In particular, we apply the model to the binary systems WR 98 and WR 113 and compare our results with the observations. Our theoretical results are in good agreement with the observed thermal spectra from these sources.

We present high angular resolution (FWHM_beam < 0.2") observations of the 23.1-GHz methanol (CH_3OH) transition toward the massive-star forming region NGC 7538 IRS 1. The two velocity components previously reported by Wilson et al. are resolved into distinct spatial features with brightness temperatures (T_B) greater than 10^4 K, proving their maser nature. Thus, NGC 7538 IRS 1 is the third region confirmed to show methanol maser emission at this frequency. The brighter 23.1-GHz spot coincides in position with a rare formaldehyde (H_2CO) maser, and marginally with a 22.2-GHz water (H_2O) maser, for which we report archival observations. The weaker CH_3OH spot coincides with an H_2O maser. The ratio of T_B for the 23.1-GHz masers to that of the well-known 12.2-GHz CH_3OH masers in this region roughly agrees with model predictions. However, the 23.1-GHz spots are offset in position from the CH_3OH masers at other frequencies. This is difficult to interpret in terms of models that assume that all the masers arise from the same clumps, but it may result from turbulent conditions within the gas or rapid variations in the background radiation field.

We present quasi-simultaneous, multi-frequency VLA observations at 4.8, 8.4, and 22.5 GHz, of a sample of 13 Wolf Rayet (WR) stars, aimed at disentangling the nature of their radio emission and the possible detection of a non-thermal behavior in close binary systems. We detected 12 stars from our sample, for which we derived spectral information and estimated their mass loss rates. From our data, we identified four thermal sources (WR 89, 113, 138, and 141), and three sources with a composite spectrum (similar contribution of thermal and non-thermal emission; WR 8, 98, and 156). On the other hand, from the comparison with previous observations, we confirm the non-thermal spectrum of one (WR 105), and also found evidence of a composite spectrum for WR 79a, 98a, 104, and 133. Finally, we discuss the possible scenarios to explain the nature of the emission for the observed objects.

The classical model for free-free emission from ionized stellar winds is based on the assumption of a stationary, isotropic and homogeneous wind. However, since there exist objects whose wind behaviour deviates from the standard model, during the last decade these assumptions have been questioned. In this work, we present results for 3 bright sources: P Cyg, Cyg OB2 No.12 and WR 147. These objects have been reported to possess winds that deviate from the basic assumptions. We have obtained flux densities, sizes, spectral indices and mass loss rates for each of the targets. These parameters allow us to analyze possible asymmetries, inhomogeneities and time variations in the flux densities. These features confirm the nonclassical behaviour of these winds.