We present multi-scale and multi-wavelength observations of the Galactic HII region S305, which is excited by massive O8.5V and O9.5V stars. Infrared images reveal an extended sphere-like shell (extension ~7.5 pc; at T_d = 17.5-27 K) enclosing the S305 HII region (size ~5.5 pc; age ~1.7 Myr). The extended structure observed in the Herschel temperature map indicates that the molecular environment of S305 is heated by the massive O-type stars. Regularly spaced molecular condensations and dust clumps are investigated toward the edges of the infrared shell, where the PAH and H_2 emission is also observed. The molecular line data show a signature of an expanding shell of molecular gas in S305. GMRT 610 and 1280 MHz continuum maps reveal overdensities of the ionized emission distributed around two O-type stars, which are surrounded by the horseshoe envelope (extension ~2.3 pc). A molecular gas deficient region/cavity is identified toward the center of the horseshoe envelope, which is well traced with PAH, H_2, molecular, and dust emission. The edges of the infrared shell are found to be located in the front of the horseshoe envelope. All these outcomes provide the observational evidence of the feedback of O-type stars in S305. Moreover, non-thermal radio emission is detected in S305 with an average spectral index alpha ~-0.45. The variations in alpha, ranging from -1.1 to 1.3, are explained due to soft synchrotron emission and either optically-thicker thermal emission at high frequencies or a suppression of the low-frequency emission by the Razin-Tsytovich effect.

We conducted an observational campaign towards one of the most massive and luminous colliding wind binaries in the Galaxy, HD~93129A, close to its periastron passage in 2018. During this time the source was predicted to be in its maximum of high-energy emission. Here we present our data analysis from the X-ray satellites \textit{Chandra} and \textit{NuSTAR} and the $\gamma$-ray satellite \textit{AGILE}. High-energy emission coincident with HD~93129A was detected in the X-ray band up to $\sim$18~keV, whereas in the $\gamma$-ray band only upper limits were obtained. We interpret the derived fluxes using a non-thermal radiative model for the wind-collision region. We establish a conservative upper limit for the fraction of the wind kinetic power that is converted into relativistic electron acceleration, $f_\mathrm{NT,e} < 0.02$. In addition, we set a lower limit for the magnetic field in the wind-collision region as $B_\mathrm{WCR} > 0.3$~G. We also argue a putative interpretation of the emission from which we estimate $f_\mathrm{NT,e} \approx 0.006$ and $B_\mathrm{WCR} \approx 0.5$~G. We conclude that multi-wavelength, dedicated observing campaigns during carefully selected epochs are a powerful tool for characterising the relativistic particle content and magnetic field intensity in colliding wind binaries.

The Argentine Institute of Radio astronomy (IAR) is equipped with two single-dish 30mts radio antennas capable of performing daily observations of pulsars and radio transients in the southern hemisphere at 1.4 GHz. We aim to introduce to the international community the upgrades performed and to show that IAR observatory has become suitable for investigations in numerous areas of pulsar radio astronomy, such as pulsar timing arrays, targeted searches of continuous gravitational waves sources, monitoring of magnetars and glitching pulsars, and studies of short time scale interstellar scintillation. We refurbished the two antennas at IAR to achieve high-quality timing observations. We gathered more than $1\,000$ hours of observations with both antennas to study the timing precision and sensitivity they can achieve. We introduce the new developments for both radio telescopes at IAR. We present observations of the millisecond pulsar J0437$-$4715 with timing precision better than 1~$\mu$s. We also present a follow-up of the reactivation of the magnetar XTE J1810--197 and the measurement and monitoring of the latest (Feb. 1st. 2019) glitch of the Vela pulsar (J0835--4510). We show that IAR is capable of performing pulsar monitoring in the 1.4 GHz radio band for long periods of time with a daily cadence. This opens the possibility of pursuing several goals in pulsar science, including coordinated multi-wavelength observations with other observatories. In particular, observations of the millisecond pulsar J0437$-$4715 will increase the gravitational wave sensitivity of the NANOGrav array in their current blind spot. We also show IAR's great potential for studying targets of opportunity and transient phenomena such as magnetars, glitches, and fast-radio-burst sources.

The supersonic winds produced by massive stars carry a large amount of kinetic power. In numerous scenarios such winds have been proven to produce shocks in which relativistic particles are accelerated emitting non-thermal radiation. Here, we report the first detection of non-thermal emission from a single stellar bubble, G2.4+1.4, associated with a WO star. We observed this source with the uGMRT in Band 4 ($550-850$ MHz) and Band 5 ($1050-1450$ MHz). We present intensity and spectral index maps for this source that are consistent with synchrotron emission (average spectral index, $\alpha = -0.83 \pm 0.10$). The fraction of the available kinetic wind power that is converted into cosmic ray acceleration is estimated to be of the order of a few per cent. This finding constitutes an observational breakthrough and gives new insight on the non-thermal physical processes taking place in the environments of isolated massive stars. In particular, our results show that non-runaway isolated massive stars are capable of accelerating relativistic particles and are therefore confirmed as sources of Galactic cosmic rays.

High-mass microquasars (HMMQs) are systems from which relativistic jets are launched. At the scales of several times the binary system size, the jets are expected to follow a helical path caused by the interaction with a strong stellar wind and orbital motion. Such a trajectory has its influence on the non-thermal emission of the jets, which also depends strongly on the observing angle due to Doppler boosting effects. We explore how the expected non-thermal emission of HMMQ jets at small scales is affected by the impact of the stellar wind and the orbital motion on the jet propagation. We studied the broadband non-thermal emission, from radio to gamma rays, produced in HMMQ jets up to a distance of several orbital separations, taking into account a realistic jet trajectory, different model parameters, and orbital modulation. The jet trajectory is computed by considering momentum transfer with the stellar wind. Electrons are injected at the position where a recollimation shock in the jets is expected due to the wind impact. Their distribution along the jet path is obtained assuming local acceleration at the recollimation shock, and cooling via adiabatic, synchrotron, and inverse Compton processes. The synchrotron and inverse Compton emission is calculated taking into account synchrotron self-absorption within the jet, free-free absorption with the stellar wind, and absorption by stellar photons via pair production. Asymmetric light curves are obtained owing to the helical trajectory of the jets. The presence of helical shaped jets could be inferred from asymmetries in the light curves, which become noticeable only for large jet Lorentz factors and low magnetic fields. Model parameters could be constrained if accurate phase-resolved light curves from GeV to TeV energies were available.

We present observations of the Cygnus OB2 region obtained with the Giant Metrewave Radio Telescope (GMRT) at the frequencies of 325 MHz and 610 MHz. In this contribution we focus on the study of proplyd-like objects (also known as free-floating Evaporating Gas Globules or frEGGs) that typically show an extended cometary morphology. We identify eight objects previously studied at other wavelengths and derive their physical properties by obtaining their optical depth at radio-wavelengths. Using their geometry and the photoionization rate needed to produce their radio-continuum emission, we find that these sources are possibly ionized by a contribution of the stars Cyg OB2 #9 and Cyg OB2 #22. Spectral index maps of the eight frEGGs were constructed, showing a flat spectrum in radio frequencies in general. We interpret these as produced by optically-thin ionized gas, although it is possible that a combination of thermal emission, not necessarily optically thin, produced by a diffuse gas component and the instrument response (which detects more diffuse emission at low frequencies) can artificially generate negative spectral indices. In particular, for the case of the Tadpole we suggest that the observed emission is not of non-thermal origin despite the presence of regions with negative spectral indices in our maps.

The massive binary system WR11 has been recently proposed as the counterpart of a Fermi source. If correct, it would be the second colliding wind binary detected in GeV gamma-rays. However, the reported flux measurements from 1.4 to 8.64GHz fail to establish the presence of non-thermal (synchrotron) emission from this source. Moreover, WR11 is not the only radio source within the Fermi detection box. Other possible counterparts have been identified in archival data, some of which present strong non-thermal radio emission. We conducted "-resolution observations towards WR11 at very low frequencies (150 to 1400~MHz) where the NT emission is expected to dominate, and present a catalog of more than 400 radio-emitters, among which a significant part is detected at more than one frequency, including limited spectral index information. A search for counterparts for this last group pointed at MOST0808-471, a source 2' away from WR11, as a promising candidate for high-energy emission, with resolved structure along 325 - 1390 MHz. For it, we reprocessed archive interferometric data up to 22.3 GHz and obtained a non-thermal radio spectral index of -0.97 +- 0.09. However, multiwavelength observations of this source are required to establish its nature and to assess whether it can produce (part of) the observed gamma-rays. WR11 spectrum follows a spectral index of 0.74 +- 0.03 from 150 to 230 GHz, consistent with thermal emission. We interpret that any putative synchrotron radiation from the colliding-wind region of this relatively short-period system is absorbed in the photospheres of the individual components. Notwithstanding, the new radio data allowed to derive a mass loss rate of 0.000025 Mo/yr, which, according to the latest models for gamma-ray emission in WR 11, would suffice to provide the required kinetic power to feed non-thermal radiation processes.

The innermost parts of powerful jets in active galactic nuclei are surrounded by dense, high-velocity clouds from the broad-line region, which may penetrate into the jet and lead to the formation of a strong shock. Such jet-cloud interactions are expected to have measurable effects on the $\gamma$-ray emission from blazars. We characterise the dynamics of a typical cloud-jet interaction scenario, and the evolution of its radiative output in the 0.1-30~GeV energy range, to assess to what extent these interactions can contribute to the $\gamma$-ray emission in blazars. We use semi-analytical descriptions of the jet-cloud dynamics, taking into account the expansion of the cloud inside the jet and its acceleration. Assuming that electrons are accelerated in the interaction and making use of the hydrodynamical information, we then compute the high-energy radiation from the cloud, including the absorption of $\gamma$-rays in the ambient photon field through pair creation. Jet-cloud interactions can lead to significant $\gamma$-ray fluxes in blazars with a broad-line region, in particular when the cloud expansion and acceleration inside the jet are taken into account. This is caused by 1) the increased shocked area in the jet, which leads to an increase in the energy budget for the non-thermal emission; 2) a more efficient inverse Compton cooling with the boosted photon field of the broad-line region; and 3) an increased observer luminosity due to Doppler boosting effects. For typical broad-line region parameters, either (i) jet-cloud interactions contribute significantly to the persistent $\gamma$-ray emission from blazars or (ii) the broad-line region is far from spherical or the fraction of energy deposited in non-thermal electrons is small.

Context: Runaway stars produce bowshocks that are usually observed at infrared (IR) wavelengths. Non-thermal radio emission has been detected so far only from the bowshock of BD+43{\deg}3654, whereas the detection of non-thermal radiation from these bowshocks at high energies remains elusive. Aims: We aim at characterising in detail the radio, X-ray, and gamma-ray emission from stellar bowshocks accounting for the structure of the region of interaction between the stellar wind and its environment. Methods: We develop a broadband-radiative, multi-zone model for stellar bowshocks that takes into account the spatial structure of the emitting region and the observational constraints. The model predicts the evolution and the emission of the relativistic particles accelerated and streaming together with the shocked flow. Results: We present broadband non-thermal spectral energy distributions for different scenarios, synthetic radio-cm synchrotron maps that reproduce the morphology of BD+43{\deg}3654, and updated predictions in X-ray and gamma-ray energy ranges. We also compare the results of the multi-zone model applied in this work with those of a refined one-zone model. Conclusions: A multi-zone model provides better constraints than a one-zone model on the relevant parameters, namely the magnetic field intensity and the amount of energy deposited in non-thermal particles. However, one-zone models can be improved by carefully characterising the intensity of the IR dust photon field and the escape rate of the plasma from the shocked region. Finally, comparing observed radio maps with those obtained from a multi-zone model enables constraints to be obtained on the direction of stellar motion with respect to the observer.

The stellar winds of the massive stars in high-mass microquasars are thought to be inhomogeneous. The interaction of these inhomogeneities, or clumps, with the jets of these objects may be a major factor in gamma-ray production. Our goal is to characterize a typical scenario of clump-jet interaction, and calculate the contribution of these interactions to the gamma-ray emission from these systems. We use axisymmetric, RHD simulations to model the emitting flow in a typical clump-jet interaction. Using the simulation results we perform a numerical calculation of the high-energy emission from one of these interactions. The radiative calculations are done for relativistic electrons locally accelerated at the jet shock, and the synchrotron and IC radiation spectra are computed for different stages of the shocked clump evolution. We also explore different parameter values, such as viewing angle and magnetic field strength. We generalize phenomenologically the results derived from one clump-jet interaction to multiple interactions under different wind models, estimating the clump-jet interaction rates, and the resulting luminosities in the GeV range. If particles are efficiently accelerated in clump-jet interactions, the apparent gamma-ray luminosity through IC scattering with the stellar photons can be significant even for rather strong magnetic fields and thus efficient synchrotron cooling. Moreover, despite the standing nature or slow motion of the jet shocks for most of the interaction stage, Doppler boosting in the postshock flow is relevant even for mildly relativistic jets. For clump-to-average wind density contrasts $\gtrsim 10$, clump-jet interactions could be bright enough to match the observed GeV luminosity in Cyg~X-1 and Cyg~X-3 when a jet is present in these sources, with required non-thermal-to-total available power fractions $\sim 0.01$ and $0.1$, respectively.

Several massive early-type binaries exhibit non-thermal emission which has been attributed to synchrotron radiation from particles accelerated by diffusive shock acceleration (DSA) in the wind-collision region (WCR). If the magnetic field in the strong shocks is ordered, its component parallel to the shock front should be enhanced, and the resultant synchrotron radiation would be polarized. However, such polarization has never been measured. We performed spatially-unresolved radio continuum observations of WR 146, a well known a WC6+O8 system, at 5 cm and 20 cm with the Karl G. Jansky Very Large Array. No linearly polarized radio emission was detected. The data constrain the fractional linear polarization to less than 0.6% between 1 to 8 GHz. This is compatible with a high level of turbulence and a dominant random component in the magnetic field. In this case the relativistic particles could be produced by turbulent magnetic reconnection (MR). In order for this scenario to satisfy the required non-thermal energy budget, we find through numerical modeling that the strength of the magnetic field in the WCR must be as high as ~150 mG. However, if the magnetic field is ordered and DSA is ongoing, then a combination of internal and external Faraday rotation could equally account for the depolarization of the emission. It is not clear whether it is possible to develop the high level of turbulence and strong magnetic fields required for efficient MR in a long-period binary such as WR 146. This scenario might also have trouble explaining the low-frequency cutoff in the spectrum. We therefore favor a scenario where particles are accelerated through DSA and the depolarization is produced by mechanisms other than a large ratio between random to regular magnetic fields.

The binary stellar system HD 93129A is one of the most massive known binaries in our Galaxy. This system presents non-thermal emission in the radio band, which can be used to infer its physical conditions and predict its emission in the high-energy band. We intend to constrain some of the unknown parameters of HD 93129A through modelling the non-thermal emitter, and also to analyse the detectability of this source in hard X-rays and $\gamma$-rays. We develop a broadband radiative model for the wind-collision region taking into account the evolution of the accelerated particles streaming along the shocked region, the emission by different radiative processes, and the attenuation of the emission propagating through the local matter and radiation fields. From the analysis of the radio emission, we find that the binary HD~93129A is more likely to have a low inclination and a high eccentricity. The minimum energy of the non-thermal electrons seems to be between $\sim 20 - 100$MeV, depending on the intensity of the magnetic field in the wind-collision region. The latter can be in the range $\sim 20 - 1500$ mG. Our model is able to reproduce the observed radio emission, and predicts that the non-thermal radiation from HD~93129A will increase in the near future. With instruments such as \textit{NuSTAR}, \textit{Fermi}, and CTA, it will be possible to constrain the relativistic particle content of the source, and other parameters such as the magnetic field strength in the wind collision zone, which in turn can be used to obtain upper-limits of the magnetic field on the surface of the very massive stars, thereby inferring whether magnetic field amplification is taking place in the particle acceleration region.

Context: Several binary systems hosting massive stars present gamma-ray emission. In most of these systems, despite detailed observational information being available, the nature and the structure of the emitter are still poorly known. Aims: We investigate the validity of the so-called one-zone approximation for the high-energy emitter in binary systems hosting a massive star. In particular, the case of LS 5039 is considered. Methods: Assuming a point-like emitter at rest, the presence of a nearby massive star, and the observed MeV and GeV fluxes as a reference, a non-thermal leptonic model is systematically applied for different locations, magnetic fields, and non-radiative losses. This allows us to identify both the emitter configurations that are most compatible with observations and inconsistencies between model predictions and the available data. Results: In the case of LS 5039, the best parameter combination is fast non-radiative cooling and a low magnetic field. However, discrepancies appear when comparing the model results at the MeV and GeV energy ranges with the observed fluxes. Predictions fail when the orbital motion is included in the analysis, because emitters and energy budgets that are too large are required. Values of X-ray and TeV fluxes that are too high are predicted in such a case, along half of the orbit. Conclusions: We show that the radiation in LS 5039 does not come from only one electron population, and the emitter is likely extended and inhomogeneous with a low magnetic field. We suggest that the emitter moves at relativistic velocities with Doppler boosting playing a significant role. Key words. Stars: high-mass, X-rays: binaries, gamma-rays: binaries