Extensive air showers, originating from ultra-high energy cosmic rays, have been successfully measured through the use of arrays of water-Cherenkov detectors (WCDs). Sophisticated analyses exploiting WCD data have made it possible to demonstrate that shower simulations, based on different hadronic-interaction models, cannot reproduce the observed number of muons at the ground. The accurate knowledge of the WCD response to muons is paramount in establishing the exact level of this discrepancy. In this work, we report on a study of the response of a WCD of the Pierre Auger Observatory to atmospheric muons performed with a hodoscope made of resistive plate chambers (RPCs), enabling us to select and reconstruct nearly 600 thousand single muon trajectories with zenith angles ranging from 0$^\circ$ to 55$^\circ$. Comparison of distributions of key observables between the hodoscope data and the predictions of dedicated simulations allows us to demonstrate the accuracy of the latter at a level of 2%. As the WCD calibration is based on its response to atmospheric muons, the hodoscope data are also exploited to show the long-term stability of the procedure.

One of the major puzzles in the study of stellar evolution is the formation process of bipolar and multi-polar planetary nebulae. There is growing consensus that collimated jets create cavities with dense walls in the slowly-expanding (10--20 ~km~s$^{-1}$) envelope ejected in previous evolutionary phases, leading to the observed morphologies. However, the launching of the jet and the way it interacts with the circumstellar material to create such asymmetric morphologies have remained poorly known. Here we present for the first time CO emission from the asymptotic giant branch star W43A that traces the whole stream of a jet, from the vicinity of its driving stellar system out to the regions where it shapes the circumstellar envelope. We found that the jet has a launch velocity of 175~km~s$^{-1}$ and decelerates to a velocity of 130~km~s$^{-1}$ as it interacts with circumstellar material. The continuum emission reveals a bipolar shell with a compact bright dot in the centre that pinpoints the location of the driving source of the jet. The kinematical ages of the jet and the bipolar shell are equal, $\tau$$\sim$60~years, indicating that they were created simultaneously, probably by a common underlying mechanism, and in an extremely short time. These results provide key initial conditions for the theoretical models that aim to explain the formation of bipolar morphologies in the circumstellar envelopes of low and intermediate mass stars.

We report Atacama-Large-millimeter-Array (ALMA) Band 6 full-polarization observations of the lensed blazar PKS1830-211 during its record-breaking radio and gamma-ray flare in the spring of 2019. The observations were taken close to the peak of the gamma activity and show a clear difference in polarization state between the two time-delayed images. The leading image has a fractional polarization about three times lower than the trailing image, implying that significant depolarization occurred during the flare. In addition, we observe clear intra-hour variability of the polarization properties between the two lensed images, with a quasi-linear increase of the differential electric-vector position angle at a rate of about two degrees per hour, associated with changes of the relative fractional polarization of ~10%. This variability, combined with the lower polarization close to the peak of gamma activity, is in agreement with models of magnetic turbulence to explain polarization variability in blazar jets. Finally, the comparison of results from the full and differential polarization analysis confirms that the differential polarization technique (Marti-Vidal et al. 2016) can provide useful information on the polarization state of sources like gravitationally lensed radio-loud quasars.

We search for signals of magnetically-induced effects in the arrival directions of ultra-high-energy cosmic rays detected at the Pierre Auger Observatory. We apply two different methods. One is a search for sets of events that show a correlation between their arrival direction and the inverse of their energy, which would be expected if they come from the same point-like source, they have the same electric charge and their deflection is relatively small and coherent. We refer to these sets of events as "multiplets". The second method, called "thrust", is a principal axis analysis aimed to detect the elongated patterns in a region of interest. We study the sensitivity of both methods using a benchmark simulation and we apply them to data in two different searches. The first search is done assuming as source candidates a list of nearby active galactic nuclei and starburst galaxies. The second is an all-sky blind search. We report the results and we find no statistically significant features. We discuss the compatibility of these results with the indications on the mass composition inferred from data of the Pierre Auger Observatory.

We present measurements of the large-scale cosmic-ray anisotropies in right ascension, using data collected by the surface detector array of the Pierre Auger Observatory over more than 14 years. We determine the equatorial dipole component, $\vec{d}_\perp$, through a Fourier analysis in right ascension that includes weights for each event so as to account for the main detector-induced systematic effects. For the energies at which the trigger efficiency of the array is small, the ``East-West'' method is employed. Besides using the data from the array with detectors separated by 1500 m, we also include data from the smaller but denser sub-array of detectors with 750 m separation, which allows us to extend the analysis down to $\sim 0.03$ EeV. The most significant equatorial dipole amplitude obtained is that in the cumulative bin above 8~EeV, $d_\perp=6.0^{+1.0}_{-0.9}$%, which is inconsistent with isotropy at the 6$\sigma$ level. In the bins below 8 EeV, we obtain 99% CL upper-bounds on $d_\perp$ at the level of 1 to 3 percent. At energies below 1 EeV, even though the amplitudes are not significant, the phases determined in most of the bins are not far from the right ascension of the Galactic center, at $\alpha_{\rm GC}=-94^\circ$, suggesting a predominantly Galactic origin for anisotropies at these energies. The reconstructed dipole phases in the energy bins above 4 EeV point instead to right ascensions that are almost opposite to the Galactic center one, indicative of an extragalactic cosmic ray origin.

One of the key questions in the study of relativistic jets is how magnetic reconnection occurs and whether it can effectively accelerate electrons in the jet. We performed 3D particle-in-cell (PIC) simulations of a relativistic electron-proton jet of relatively large radius that carries a helical magnetic field. We focussed our investigation on the interaction between the jet and the ambient plasma and explore how the helical magnetic field affects the excitation of kinetic instabilities such as the Weibel instability (WI), the kinetic Kelvin-Helmholtz instability (kKHI), and the mushroom instability (MI). In our simulations these kinetic instabilities are indeed excited, and particles are accelerated. At the linear stage we observe recollimation shocks near the center of the jet. As the electron-proton jet evolves into the deep nonlinear stage, the helical magnetic field becomes untangled due to reconnection-like phenomena, and electrons are repeatedly accelerated as they encounter magnetic-reconnection events in the turbulent magnetic field.

We present an initial study of one of the first ALMA Band 3 observations of the Sun with the aim to characterise the diagnostic potential of brightness temperatures measured with ALMA on the Sun. The observation covers 48min at a cadence of 2s targeting a Quiet Sun region at disk-centre. Corresponding time series of brightness temperature maps are constructed with the first version of the Solar ALMA Pipeline (SoAP) and compared to simultaneous SDO observations. The angular resolution of the observations is set by the synthesized beam (1.4x2.1as). The ALMA maps exhibit network patches, internetwork regions and also elongated thin features that are connected to large-scale magnetic loops as confirmed by a comparison with SDO maps. The ALMA Band 3 maps correlate best with the SDO/AIA 171, 131 and 304 channels in that they exhibit network features and, although very weak in the ALMA maps, imprints of large-scale loops. A group of compact magnetic loops is very clearly visible in ALMA Band 3. The brightness temperatures in the loop tops reach values of about 8000-9000K and in extreme moments up to 10 000K. ALMA Band 3 interferometric observations from early observing cycles already reveal temperature differences in the solar chromosphere. The weak imprint of magnetic loops and the correlation with the 171, 131, and 304 SDO channels suggests though that the radiation mapped in ALMA Band 3 might have contributions from a larger range of atmospheric heights than previously assumed but the exact formation height of Band 3 needs to be investigated in more detail. The absolute brightness temperature scale as set by Total Power measurements remains less certain and must be improved in the future. Despite these complications and the limited angular resolution, ALMA Band 3 observations have large potential for quantitative studies of the small-scale structure and dynamics of the solar chromosphere.

PG 1553+113 is the first blazar showing an approximately two-year quasi-periodic pattern in its gamma-ray light curve. Such quasi-periodicity might have a geometrical origin, possibly related to the precessing nature of the jet, or could be intrinsic to the source and related to pulsational accretion flow instabilities. By means of a ~2yr very long baseline array (VLBA) monitoring at 15, 24, and 43 GHz we investigate the source pc-scale properties during an entire cycle of gamma-ray activity in the period 2015-2017. In contrast to the well-defined periodicity in the gamma-ray emission, at radio frequencies no clear periodic pattern can be recognized. The jet position angle, constrained by means of the total intensity ridge line, varies across the different observing epochs in the range 40-60 deg. We also investigate the time evolution of the source polarization properties, including the rotation measure. The brightness temperature is found to decrease as the frequency increases with an intrinsic value of ~1.5 x 10^10 K and the estimated Doppler factor is ~1.4.

We present high spatial resolution observations of ground-state OH masers achieved with the Australia Telescope Compact Array (ATCA). These observations targeted 253 pointing centres containing OH maser candidates at all four ground-state OH transitions identified in the Southern Parkes Large-Area Survey in Hydroxyl (SPLASH) across 96 square degrees of the Southern Galactic plane (332\degree$<l<$334\degree and $-$2\degree$<b<+$2\degree, 344\degree$<l<$355\degree and $-$2\degree$<b<+$2\degree, 358\degree$<l<$4\degree and $+$2\degree$<b<+$6\degree, 5\degree$<l<$10\degree and $-$2\degree$<b<+$2\degree). We detect maser emission towards 236 fields and suggest that 7 out of 17 non-detections are due to the slightly lower sensitivity of the ATCA observations, combined with some temporal variability. The superior resolution provided by the ATCA data has allowed us to identify 362 OH maser sites in the 236 target fields. Almost half (160 of 362) of these masers have been detected for the first time. Comparison between these 362 maser sites with information presented in the literature allowed us to categorize 238 sites as evolved star sites (66\%), 63 as star formation (17\%), eight as supernova remnants and 53 unknown maser sites (15\%). We present analysis of the OH masers across the full SPLASH survey range (176 square degrees) and find that the detection rate of 1.7 GHz radio continuum sources (18\%) is lower than that previously found at 8.2 and 9.2 GHz (38\%). We also find that the velocity separations of evolved star sites with symmetric 1612 MHz maser profiles are generally smaller than those with asymmetric profiles.

We present a multiwavelength study of the flat-spectrum radio quasar CTA 102 during 2013-2017. We use radio-to-optical data obtained by the Whole Earth Blazar Telescope, 15 GHz data from the Owens Valley Radio Observatory, 91 and 103 GHz data from the Atacama Large Millimeter Array, near-infrared data from the Rapid Eye Monitor telescope, as well as data from the Swift (optical-UV and X-rays) and Fermi ($\gamma$ rays) satellites to study flux and spectral variability and the correlation between flux changes at different wavelengths. Unprecedented $\gamma$-ray flaring activity was observed during 2016 November-2017 February, with four major outbursts. A peak flux of (2158 $\pm$ 63)$\times$10$^{-8}$ ph cm$^{-2}$ s$^{-1}$, corresponding to a luminosity of (2.2 $\pm$ 0.1)$\times$10$^{50}$ erg s$^{-1}$, was reached on 2016 December 28. These four $\gamma$-ray outbursts have corresponding events in the near-infrared, optical, and UV bands, with the peaks observed at the same time. A general agreement between X-ray and $\gamma$-ray activity is found. The $\gamma$-ray flux variations show a general, strong correlation with the optical ones with no time lag between the two bands and a comparable variability amplitude. This $\gamma$-ray/optical relationship is in agreement with the geometrical model that has successfully explained the low-energy flux and spectral behaviour, suggesting that the long-term flux variations are mainly due to changes in the Doppler factor produced by variations of the viewing angle of the emitting regions. The difference in behaviour between radio and higher energy emission would be ascribed to different viewing angles of the jet regions producing their emission.

With the Surface Detector array (SD) of the Pierre Auger Observatory we can detect neutrinos with energy between $10^{17}\,$eV and $10^{20}\,$eV from point-like sources across the sky, from close to the Southern Celestial Pole up to $60^\circ$ in declination, with peak sensitivities at declinations around $\sim -53^\circ$ and $\sim+55^\circ$, and an unmatched sensitivity for arrival directions in the Northern hemisphere. A search has been performed for highly-inclined air showers induced by neutrinos of all flavours with no candidate events found in data taken between 1 Jan 2004 and 31 Aug 2018. Upper limits on the neutrino flux from point-like steady sources have been derived as a function of source declination. An unrivaled sensitivity is achieved in searches for transient sources with emission lasting over an hour or less, if they occur within the field of view corresponding to the zenith angle range between $60^\circ$ and $~95^\circ$ where the SD of the Pierre Auger Observatory is most sensitive to neutrinos.

The determination of the primary energy of extensive air showers using the fluorescence detection technique requires an estimation of the energy carried away by particles that do not deposit all their energy in the atmosphere. This estimation is typically made using Monte Carlo simulations and thus depends on the assumed primary particle mass and on model predictions for neutrino and muon production. In this work we present a new method to obtain the invisible energy from events detected by the Pierre Auger Observatory. The method uses measurements of the muon number at ground level, and it allows us to reduce significantly the systematic uncertainties related to the mass composition and the high energy hadronic interaction models, and consequently to improve the estimation of the energy scale of the Observatory.

Neutrinos with energies above $10^{17}$ eV are detectable with the Surface Detector Array of the Pierre Auger Observatory. The identification is efficiently performed for neutrinos of all flavors interacting in the atmosphere at large zenith angles, as well as for Earth-skimming $\tau$ neutrinos with nearly tangential trajectories relative to the earth. No neutrino candidates were found in $\sim\,14.7$ years of data taken up to 31 August 2018. This leads to restrictive upper bounds on their flux. The $90\%$ C.L. single-flavor limit to the diffuse flux of ultra-high-energy neutrinos with an $E_\nu^{-2}$ spectrum in the energy range $1.0 \times 10^{17}~{\rm eV} - 2.5 \times 10^{19}~{\rm eV}$ is $E^2 {\rm d}N_\nu/{\rm d}E_\nu < 4.4 \times 10^{-9}~{\rm GeV~cm^{-2}~s^{-1}~sr^{-1}}$, placing strong constraints on several models of neutrino production at EeV energies and on the properties of the sources of ultra-high-energy cosmic rays.

Contributions of the Pierre Auger Collaboration to the 36th International Cosmic Ray Conference (ICRC 2019), 24 July - 1 August 2019, Madison, Wisconsin, USA.

The Event Horizon Telescope (EHT) recently produced the first horizon-scale image of a supermassive black hole. Expanding the array to include a 3-meter space telescope operating at >200 GHz enables mass measurements of many black holes, movies of black hole accretion flows, and new tests of general relativity that are impossible from the ground.

We present a detailed multi-wavelength characterization of the multi-ring disk of HD 169142. We report new ALMA observations at 3 mm and analyze them together with archival 0.89 and 1.3 mm data. Our observations resolve three out of the four rings in the disk previously seen in high-resolution ALMA data. A simple parametric model is used to estimate the radial profile of the dust optical depth, temperature, density, and particle size distribution. We find that the multiple ring features of the disk are produced by annular accumulations of large particles, probably associated with gas pressure bumps. Our model indicates that the maximum dust grain size in the rings is $\sim1$ cm, with slightly flatter power-law size distributions than the ISM-like size distribution ($p\sim3.5$) found in the gaps. In particular, the inner ring ($\sim26$ au) is associated with a strong and narrow buildup of dust particles that could harbor the necessary conditions to trigger the streaming instability. According to our analysis, the snowlines of the most important volatiles do not coincide with the observed substructures. We explore different ring formation mechanisms and find that planet-disk interactions are the most likely scenario to explain the main features of HD 169142. Overall, our multi-wavelength analysis provides some of the first unambiguous evidence of the presence of radial dust traps in the rings of HD 169142. A similar analysis in a larger sample of disks could provide key insights on the impact that disk substructures have on the dust evolution and planet formation processes.

The profile of the longitudinal development of showers produced by ultra-high energy cosmic rays carries information related to the interaction properties of the primary particles with atmospheric nuclei. In this work, we present the first measurement of the average shower profile in traversed atmospheric depth at the Pierre Auger Observatory. The shapes of profiles are well reproduced by the Gaisser-Hillas parametrization within the range studied, for E > 10^{17.8} eV. A detailed analysis of the systematic uncertainties is performed using 10 years of data and a full detector simulation. The average shape is quantified using two variables related to the width and asymmetry of the profile, and the results are compared with predictions of hadronic interaction models for different primary particles.

An overview of the multi-messenger capabilities of the Pierre Auger Observatory is presented. The techniques and performance of searching for Ultra-High Energy neutrinos, photons and neutrons are described. Some of the most relevant results are reviewed, such as stringent upper bounds that were placed to a flux of diffuse cosmogenic neutrinos and photons, bounds placed on neutrinos emitted from compact binary mergers that were detected by LIGO and Virgo during their first and second observing runs, as well as searches for high energy photons and neutrons from the Galactic center that constrain the properties of the putative Galactic PeVatron, observed by the H.E.S.S.\ collaboration. The observation of directional correlations between ultra-high energy cosmic rays and either high energy astrophysical neutrinos or specific source populations, weighted by their electromagnetic radiation, are also discussed. They constitute additional multi-messenger approaches aimed at identifying the sources of high energy cosmic rays.

Water fountains are evolved stars showing early stages of collimated mass loss during transition from the asymptotic giant branch, providing valuable insight into the formation of asymmetric planetary nebulae. We report the results of multi-epoch VLBI observations, which determine the spatial and three-dimensional kinematic structure of H2O masers associated with the water fountain IRAS 18113-2503. The masers trace three pairs of high-velocity (~150-300 km/s) bipolar bow shocks on a scale of 0.18'' (~2000 au). The expansion velocities of the bow shocks exhibit an exponential decrease as a function of distance from the central star, which can be explained by an episodic, jet-driven outflow decelerating due to drag forces in a circumstellar envelope. Using our model, we estimate an initial ejection velocity ~840 km/s, a period for the ejections ~10 yr, with the youngest being ~12 yr old, and an average envelope density within the H2O maser region n(H2) ~ 10^6 cm^(-3). We hypothesize that IRAS 18113-2503 hosts a binary central star with a separation of ~10 au, revealing novel clues about the launching mechanisms of high-velocity collimated outflows in water fountains.

We present continuum and molecular line (CO, C$^{18}$O, HCO$^+$) observations carried out with the Atacama Large Millimeter/submillimeter Array toward the "water fountain" star IRAS 15103-5754, an object that could be the youngest PN known. We detect two continuum sources, separated by $0.39\pm 0.03$ arcsec. The emission from the brighter source seems to arise mainly from ionized gas, thus confirming the PN nature of the object. The molecular line emission is dominated by a circumstellar torus with a diameter of $\simeq 0.6$ arcsec (2000 au) and expanding at $\simeq 23$ km s$^{-1}$. We see at least two gas outflows. The highest-velocity outflow (deprojected velocities up to 250 km s$^{-1}$), traced by the CO lines, shows a biconical morphology, whose axis is misaligned $\simeq 14^\circ$ with respect to the symmetry axis of the torus, and with a different central velocity (by $\simeq 8$ km s$^{-1}$). An additional high-density outflow (traced by HCO$^+$) is oriented nearly perpendicular to the torus. We speculate that IRAS 15103-5754 was a triple stellar system that went through a common envelope phase, and one of the components was ejected in this process. A subsequent low-collimation wind from the remaining binary stripped out gas from the torus, creating the conical outflow. The high velocity of the outflow suggests that the momentum transfer from the wind was extremely efficient, or that we are witnessing a very energetic mass-loss event.

Our understanding of protoplanetary disks is rapidly departing from the classical view of a smooth, axisymmetric disk. This is in part thanks to the high angular resolution that (sub)mm observations can provide. Here we present the combined results of ALMA (0.9 mm) and VLA (7 mm) dust continuum observations toward the protoplanetary disk around the solar analogue GM Aur. Both images clearly resolve the $\sim$35 au inner cavity. The ALMA observations also reveal a fainter disk that extends up to $\sim250$ au. We model our observations using two approaches: an analytical fit to the observed deprojected visibilities, and a physical disk model that fits the SED as well as the VLA and ALMA observations. Despite not being evident in the deconvolved images, the VLA and ALMA visibilities can only be fitted with two bright rings of radii $\sim$40 and $\sim$80 au. Our physical model indicates that this morphology is the result of an accumulation or trapping of large dust grains, probably due to the presence of two pressure bumps in the disk. Even though alternative mechanisms cannot be discarded, the multiple rings suggest that forming planets may have cleared at least two gaps in the disk. Finally, our analysis suggests that the inner cavity might display different sizes at 0.9 mm and 7 mm. This discrepancy could be caused by the presence of free-free emission close to the star at 7 mm, or by a more compact accumulation of the large dust grains at the edge of the cavity.

We present the results of our power spectral density analysis for the BL Lac object OJ\,287, utilizing the {\it Fermi}-LAT survey at high-energy $\gamma$-rays, {\it Swift}-XRT in X-rays, several ground-based telescopes and the {\it Kepler} satellite in the optical, and radio telescopes at GHz frequencies. The light curves are modeled in terms of continuous-time auto-regressive moving average (CARMA) processes. Owing to the inclusion of the {\it Kepler} data, we were able to construct \emph{for the first time} the optical variability power spectrum of a blazar without any gaps across $\sim6$ dex in temporal frequencies. Our analysis reveals that the radio power spectra are of a colored-noise type on timescales ranging from tens of years down to months, with no evidence for breaks or other spectral features. The overall optical power spectrum is also consistent with a colored noise on the variability timescales ranging from 117 years down to hours, with no hints of any quasi-periodic oscillations. The X-ray power spectrum resembles the radio and optical power spectra on the analogous timescales ranging from tens of years down to months. Finally, the $\gamma$-ray power spectrum is noticeably different from the radio, optical, and X-ray power spectra of the source: we have detected a characteristic relaxation timescale in the {\it Fermi}-LAT data, corresponding to $\sim 150$\,days, such that on timescales longer than this, the power spectrum is consistent with uncorrelated (white) noise, while on shorter variability timescales there is correlated (colored) noise.

Here we present deep (16 mumJy), very high (40 mas) angular resolution 1.14 mm, polarimetric, Atacama Large Millimeter/submillimeter Array (ALMA) observations towards the massive protostar driving the HH 80-81 radio jet. The observations clearly resolve the disk oriented perpendicular to the radio jet, with a radius of ~0.171 arcsec (~291 au at 1.7 kpc distance). The continuum brightness temperature, the intensity profile, and the polarization properties clearly indicate that the disk is optically thick for a radius of R<170 au. The linear polarization of the dust emission is detected almost all along the disk and its properties suggest that dust polarization is produced mainly by self-scattering. However, the polarization pattern presents a clear differentiation between the inner (optically thick) part of the disk and the outer (optically thin) region of the disk, with a sharp transition that occurs at a radius of 0.1 arcsec (~170 au). The polarization characteristics of the inner disk suggest that dust settling has not occurred yet with a maximum dust grain size between 50 and 500 mum. The outer part of the disk has a clear azimuthal pattern but with a significantly higher polarization fraction compared to the inner disk. This pattern is broadly consistent with self-scattering of a radiation field that is beamed radially outward, as expected in the optically thin outer region, although contribution from non-spherical grains aligned with respect to the radiative flux cannot be excluded.

We present total and linearly polarized 3 mm Global mm-VLBI Array images of a sample of blazars and radio galaxies from the VLBA-BU-BLAZAR 7 mm monitoring program designed to probe the innermost regions of active galactic nuclei (AGN) jets and locate the sites of gamma-ray emission observed by the Fermi-LAT. The lower opacity at 3 mm and improved angular resolution, on the order of 50 microarcseconds, allow us to distinguish features in the jet not visible in the 7 mm VLBA data. We also compare two different methods used for the calibration of instrumental polarisation and we analyze the resulting images for some of the sources in the sample.

AMIGA (Auger Muons and Infill for the Ground Array) is an upgrade of the Pierre Auger Observatory designed to extend its energy range of detection and to directly measure the muon content of the cosmic ray primary particle showers. The array will be formed by an infill of surface water-Cherenkov detectors associated with buried scintillation counters employed for muon counting. Each counter is composed of three scintillation modules, with a 10 m$^2$ detection area per module. In this paper, a new generation of detectors, replacing the current multi-pixel photomultiplier tube (PMT) with silicon photo sensors (aka. SiPMs), is proposed. The selection of the new device and its front-end electronics is explained. A method to calibrate the counting system that ensures the performance of the detector is detailed. This method has the advantage of being able to be carried out in a remote place such as the one where the detectors are deployed. High efficiency results, i.e. 98 % efficiency for the highest tested overvoltage, combined with a low probability of accidental counting ($\sim$2 %), show a promising performance for this new system.

We analyse the circular polarisation data accumulated in the first 7 years of the POLAMI project introduced in an accompanying paper (Agudo et al.). In the 3mm wavelength band, we acquired more than 2600 observations, and all but one of our 37 sample sources were detected, most of them several times. For most sources, the observed distribution of the degree of circular polarisation is broader than that of unpolarised calibrators, indicating that weak (<0.5%) circular polarisation is present most of the time. Our detection rate and the maximum degree of polarisation found, 2.0%, are comparable to previous surveys, all made at much longer wavelengths. We argue that the process generating circular polarisation must not be strongly wavelength dependent, and we propose that the widespread presence of circular polarisation in our short wavelength sample dominated by blazars is mostly due to Faraday conversion of the linearly polarised synchrotron radiation in the helical magnetic field of the jet. Circular polarisation is variable, most notably on time scales comparable to or shorter than our median sampling interval <1 month. Longer time scales of about one year are occasionally detected, but severely limited by the weakness of the signal. At variance with some longer wavelength investigations we find that the sign of circular polarisation changes in most sources, while only 7 sources, including 3 already known, have a strong preference for one sign. The degrees of circular and linear polarisation do not show any systematic correlation. We do find however one particular event where the two polarisation degrees vary in synchronism during a time span of 0.9 years. The paper also describes a novel method for calibrating the sign of circular polarisation observations.

We report on the first results of the POLAMI program, a simultaneous 3.5 and 1.3mm full-Stokes-polarisation monitoring of a sample of 36 of the brightest active galactic nuclei in the northern sky with the IRAM 30m Telescope. Through a systematic statistical study of data taken from October 2006 (from December 2009 for the case of the 1.3mm observations) to August 2014, we characterise the variability of the total flux density and linear polarisation. We find that all sources in the sample are highly variable in total flux density at both 3.5 and 1.3mm, as well as in spectral index, that is found to be optically thin in general. The total flux-density variability at 1.3mm is found, in general, to be faster, and to have larger amplitude and flatter PSD slopes than 3.5mm. The polarisation degree is on average larger at 1.3mm than at 3.5mm, by a factor of 2.6. The variability of linear polarisation degree is faster and has higher fractional amplitude than for total flux density, with the typical time scales during prominent polarisation peaks being significantly faster at 1.3mm than at 3.5mm. The polarisation angle at both 3.5 and 1.3mm is highly variable. Most of the sources show one or two excursions of >180 deg. on time scales from a few weeks to about a year during the course of our observations. The 3.5 and 1.3mm polarisation angle evolution follow rather well each other, although the 1.3mm data show a clear preference to more prominent variability on the short time scales, i.e. weeks. The data are compatible with multi-zone models of conical jets involving smaller emission regions for the shortest-wavelength emitting sites. Such smaller emitting regions should also be more efficient in energising particle populations. The data also favours the integrated emission at 1.3mm to have better ordered magnetic fields than the one at 3.5mm.

We describe the POLAMI program for the monitoring of all four Stokes parameters of a sample of bright radio-loud active galactic nuclei with the IRAM 30m telescope at 3.5 and 1.3mm. The program started in October 2006 and accumulated, until August 2014, 2300 observations at 3.5mm, achieving a median time sampling interval of 22 days for the sample of 37 sources. This first paper explains the source selection, mostly blazars, the observing strategy and data calibration, and gives the details of the instrumental polarisation corrections. The sensitivity (1sigma) reached at 3.5mm is 0.5% (linear polarisation degree), 4.7 deg. (polarisation angle), and 0.23% (circular polarisation), while the corresponding values at 1.3mm are 1.7%, 9.9 deg., and 0.72%, respectively. The data quality is demonstrated by the time sequences of our calibrators Mars and Uranus. For the quasar 3C286, widely used as a linear polarisation calibrator, we give improved estimates of its linear polarisation, and show for the first time occasional detections of its weak circular polarisation, which suggests a small level of variability of the source at millimeter wavelengths.

We present interferometric, full-polarization observations of the four ground-state transitions of OH, toward five confirmed and one candidate OH-emitting planetary nebulae (OHPNe). OHPNe are believed to be very young PNe, and information on their magnetic fields (provided by their polarization) could be key to understand the early evolution of PNe. We detect significant circular and linear polarization in four and two objects, respectively. Possible Zeeman pairs are seen in JaSt 23 and IRAS 17393-2727, resulting in estimates of magnetic field strengths between 0.8 and 24 mG. We also report the new detection of OH emission at 1720 MHz toward Vy 2-2, making it the third known PN with this type of emission. We suggest that younger PNe have spectra dominated by narrow maser features and higher degrees of polarization. Shock-excited emission at 1720 MHz seems to be more common in PNe than in early evolutionary phases, and could be related to equatorial ejections during the early PN phase.

In this study we investigate jet interaction at a microscopic level in a cosmological environment, which responds to a key open question in the study of relativistic jets. Using small simulation systems during prior research, we initially studied the evolution of both electron-proton and electron-positron relativistic jets containing helical magnetic fields, by focusing on their interactions with an ambient plasma. Here, using larger jet radii, we have performed simulations of global jets containing helical magnetic fields in order to examine how helical magnetic fields affect kinetic instabilities such as the Weibel instability, the kinetic Kelvin-Helmholtz instability (kKHI) and the Mushroom instability (MI). We found that the evolution of global jets strongly depends on the size of the jet radius. For example, phase bunching of jet electrons, in particular in the electron-proton jet, is mixed with larger jet radius due to the more complicated structures of magnetic fields with excited kinetic instabilities. In our simulation study these kinetic instabilities lead to new types of instabilities in global jets. In the electron-proton jet simulation a modified recollimation occurs and jet electrons are strongly perturbed. In the electron-positron jet simulation mixed kinetic instabilities occur at early times followed by a turbulence-like structure. Simulations using much larger (and longer) systems are further required in order to thoroughly investigate the evolution of global jets containing helical magnetic fields.

Observations of fast TeV $\gamma$-ray flares from blazars reveal the extreme compactness of emitting regions in blazar jets. Combined with very-long-baseline radio interferometry measurements, they probe the structure and emission mechanism of the jet. We report on a fast TeV $\gamma$-ray flare from BL Lacertae observed by VERITAS, with a rise time of about 2.3 hours and a decay time of about 36 minutes. The peak flux at $>$200 GeV measured with the 4-minute binned light curve is $(4.2 \pm 0.6) \times 10^{-6} \;\text{photons} \;\text{m}^{-2}\, \text{s}^{-1}$, or $\sim$180% the Crab Nebula flux. Variability in GeV $\gamma$-ray, X-ray, and optical flux, as well as in optical and radio polarization was observed around the time of the TeV $\gamma$-ray flare. A possible superluminal knot was identified in the VLBA observations at 43 GHz. The flare constrains the size of the emitting region, and is consistent with several theoretical models with stationary shocks.

Water fountain stars (WFs) are evolved objects with water masers tracing high-velocity jets (up to several hundreds of km s$^{-1}$). They could represent one of the first manifestations of collimated mass-loss in evolved objects and thus, be a key to understanding the shaping mechanisms of planetary nebulae. Only 13 objects had been confirmed so far as WFs with interferometer observations. We present new observations with the Australia Telescope Compact Array and archival observations with the Very Large Array of four objects that are considered to be WF candidates, mainly based on single-dish observations. We confirm IRAS 17291-2147 and IRAS 18596+0315 (OH 37.1-0.8) as bona fide members of the WF class, with high-velocity water maser emission consistent with tracing bipolar jets. We argue that IRAS 15544-5332 has been wrongly considered as a WF in previous works, since we see no evidence in our data nor in the literature that this object harbours high-velocity water maser emission. In the case of IRAS 19067+0811, we did not detect any water maser emission, so its confirmation as a WF is still pending. With the result of this work, there are 15 objects that can be considered confirmed WFs. We speculate that there is no significant physical difference between WFs and obscured post-AGB stars in general. The absence of high-velocity water maser emission in some obscured post-AGB stars could be attributed to a variability or orientation effect.

We carried out multiwavelength (0.7-5 cm), multiepoch (1994-2015) Very Large Array (VLA) observations toward the region enclosing the bright far-IR sources FIR 3 (HOPS 370) and FIR 4 (HOPS 108) in OMC-2. We report the detection of 10 radio sources, seven of them identified as young stellar objects. We image a well-collimated radio jet with a thermal free-free core (VLA 11) associated with the Class I intermediate-mass protostar HOPS 370. The jet presents several knots (VLA 12N, 12C, 12S) of non-thermal radio emission (likely synchrotron from shock-accelerated relativistic electrons) at distances of ~7,500-12,500 au from the protostar, in a region where other shock tracers have been previously identified. These knots are moving away from the HOPS 370 protostar at ~ 100 km/s. The Class 0 protostar HOPS 108, which itself is detected as an independent, kinematically decoupled radio source, falls in the path of these non-thermal radio knots. These results favor the previously proposed scenario where the formation of HOPS 108 has been triggered by the impact of the HOPS 370 outflow with a dense clump. However, HOPS 108 presents a large proper motion velocity of ~ 30 km/s, similar to that of other runaway stars in Orion, whose origin would be puzzling within this scenario. Alternatively, an apparent proper motion could result because of changes in the position of the centroid of the source due to blending with nearby extended emission, variations in the source shape, and /or opacity effects.

We present a survey for water maser emission toward a sample of 44 low-luminosity young objects, comprising (proto-)brown dwarfs, first hydrostatic cores (FHCs), and other young stellar objects (YSOs) with bolometric luminosities lower than 0.4 L$_\odot$. Water maser emission is a good tracer of energetic processes, such as mass-loss and/or accretion, and is a useful tool to study this processes with very high angular resolution. This type of emission has been confirmed in objects with L$_{\rm bol}\ge 1$ L$_\odot$. Objects with lower luminosities also undergo mass-loss and accretion, and thus, are prospective sites of maser emission. Our sensitive single-dish observations provided a single detection when pointing toward the FHC L1448 IRS 2E. However, follow-up interferometric observations showed water maser emission associated with the nearby YSO L1448 IRS 2 { (a Class 0 protostar of L$_{\rm bol}\simeq 3.6-5.3$ L$_\odot$)}, and did not find any emission toward L1448 IRS 2E. The upper limits for water maser emission determined by our observations are one order of magnitude lower than expected from the correlation between water maser luminosities and bolometric luminosities found for YSOs. This suggests that this correlation does not hold at the lower end of the (sub)stellar mass spectrum. Possible reasons are that the slope of this correlation is steeper at L$_{\rm bol}\le 1$ L$_\odot$, or that there is an absolute luminosity threshold below which water maser emission cannot be produced. Alternatively, if the correlation still stands at low luminosity, the detection rates of masers would be significantly lower than the values obtained in higher-luminosity Class 0 protostars.

We report Very Large Array observations at 7 mm, 9 mm, and 3 cm toward the pre-transitional disk of the Herbig Ae star HD 169142. These observations have allowed us to study the mm emission of this disk with the highest angular resolution so far ($0\rlap."12\times0\rlap."09$, or 14 au$\times$11 au, at 7 mm). Our 7 and 9 mm images show a narrow ring of emission at a radius of $\sim25$ au tracing the outer edge of the inner gap. This ring presents an asymmetric morphology that could be produced by dynamical interactions between the disk and forming planets. Additionally, the azimuthally averaged radial intensity profiles of the 7 and 9 mm images confirm the presence of the previously reported gap at $\sim45$ au, and reveal a new gap at $\sim85$ au. We analyzed archival DCO$^+$(3-2) and C$^{18}$O(2-1) ALMA observations, showing that the CO snowline is located very close to this third outer gap. This suggests that growth and accumulation of large dust grains close to the CO snowline could be the mechanism responsible for this proposed outer gap. Finally, a compact source of emission is detected at 7 mm, 9 mm, and 3 cm toward the center of the disk. Its flux density and spectral index indicate that it is dominated by free-free emission from ionized gas, which could be associated with either the photoionization of the inner disk, an independent object, or an ionized jet.

On September 14, 2015 the Advanced LIGO detectors observed their first gravitational-wave (GW) transient GW150914. This was followed by a second GW event observed on December 26, 2015. Both events were inferred to have arisen from the merger of black holes in binary systems. Such a system may emit neutrinos if there are magnetic fields and disk debris remaining from the formation of the two black holes. With the surface detector array of the Pierre Auger Observatory we can search for neutrinos with energy above 100 PeV from point-like sources across the sky with equatorial declination from about -65 deg. to +60 deg., and in particular from a fraction of the 90% confidence-level (CL) inferred positions in the sky of GW150914 and GW151226. A targeted search for highly-inclined extensive air showers, produced either by interactions of downward-going neutrinos of all flavors in the atmosphere or by the decays of tau leptons originating from tau-neutrino interactions in the Earth's crust (Earth-skimming neutrinos), yielded no candidates in the Auger data collected within $\pm 500$ s around or 1 day after the coordinated universal time (UTC) of GW150914 and GW151226, as well as in the same search periods relative to the UTC time of the GW candidate event LVT151012. From the non-observation we constrain the amount of energy radiated in ultrahigh-energy neutrinos from such remarkable events.

We report a first measurement for ultra-high energy cosmic rays of the correlation between the depth of shower maximum and the signal in the water Cherenkov stations of air-showers registered simultaneously by the fluorescence and the surface detectors of the Pierre Auger Observatory. Such a correlation measurement is a unique feature of a hybrid air-shower observatory with sensitivity to both the electromagnetic and muonic components. It allows an accurate determination of the spread of primary masses in the cosmic-ray flux. Up till now, constraints on the spread of primary masses have been dominated by systematic uncertainties. The present correlation measurement is not affected by systematics in the measurement of the depth of shower maximum or the signal in the water Cherenkov stations. The analysis relies on general characteristics of air showers and is thus robust also with respect to uncertainties in hadronic event generators. The observed correlation in the energy range around the `ankle' at $\lg(E/{\rm eV})=18.5-19.0$ differs significantly from expectations for pure primary cosmic-ray compositions. A light composition made up of proton and helium only is equally inconsistent with observations. The data are explained well by a mixed composition including nuclei with mass $A > 4$. Scenarios such as the proton dip model, with almost pure compositions, are thus disfavoured as the sole explanation of the ultrahigh-energy cosmic-ray flux at Earth.

Ultrahigh energy cosmic ray air showers probe particle physics at energies beyond the reach of accelerators. Here we introduce a new method to test hadronic interaction models without relying on the absolute energy calibration, and apply it to events with primary energy 6-16 EeV (E_CM = 110-170 TeV), whose longitudinal development and lateral distribution were simultaneously measured by the Pierre Auger Observatory. The average hadronic shower is 1.33 +- 0.16 (1.61 +- 0.21) times larger than predicted using the leading LHC-tuned models EPOS-LHC (QGSJetII-04), with a corresponding excess of muons.

We present a search for ultra-relativistic magnetic monopoles with the Pierre Auger Observatory. Such particles, possibly a relic of phase transitions in the early universe, would deposit a large amount of energy along their path through the atmosphere, comparable to that of ultrahigh-energy cosmic rays (UHECRs). The air shower profile of a magnetic monopole can be effectively distinguished by the fluorescence detector from that of standard UHECRs. No candidate was found in the data collected between 2004 and 2012, with an expected background of less than 0.1 event from UHECRs. The corresponding 90% confidence level (C.L.) upper limits on the flux of ultra-relativistic magnetic monopoles range from $10^{-19}$ (cm$^{2}$ sr s)$^{-1}$ for a Lorentz factor $\gamma=10^9$ to $2.5 \times10^{-21}$ (cm$^{2}$ sr s)$^{-1}$ for $\gamma=10^{12}$. These results - the first obtained with a UHECR detector - improve previously published limits by up to an order of magnitude.

In the study of relativistic jets one of the key open questions is their interaction with the environment on the microscopic level. Here, we study the initial evolution of both electron$-$proton ($e^{-}-p^{+}$) and electron$-$positron ($e^{\pm}$) relativistic jets containing helical magnetic fields, focusing on their interaction with an ambient plasma. We have performed simulations of "global" jets containing helical magnetic fields in order to examine how helical magnetic fields affect kinetic instabilities such as the Weibel instability, the kinetic Kelvin-Helmholtz instability (kKHI) and the Mushroom instability (MI). In our initial simulation study these kinetic instabilities are suppressed and new types of instabilities can grow. In the $e^{-}-p^{+}$ jet simulation a recollimation-like instability occurs and jet electrons are strongly perturbed. In the $e^{\pm}$ jet simulation a recollimation-like instability occurs at early times followed by a kinetic instability and the general structure is similar to a simulation without helical magnetic field. Simulations using much larger systems are required in order to thoroughly follow the evolution of global jets containing helical magnetic fields.

Photoevaporation is probably the main agent for gas dispersal during the last stages of protoplanetary disk evolution. However, the overall mass loss rate in the photoevaporative wind and its driving mechanism are still not well understood. Here we report multi-configuration Very Large Array observations at 0.7, 3, and 5 cm towards the transitional disk of GM Aur. Our radio continuum observations allow us to image and spatially resolve, for the first time, the three main components at work in this stage of the disk evolution: the disk of dust, the ionized radio jet perpendicular to it, and the photoevaporative wind arising from the disk. The mass loss rate inferred from the flux density of the radio jet is consistent with the ratio between ejection and accretion rates found in younger objects, suggesting that transitional disks can power collimated ejections of material apparently following the same physical mechanisms as much younger protostars. Our results indicate that extreme-UV (EUV) radiation is the main ionizing mechanism of the photoevaporative wind traced by the free-free emission. The required low EUV photon luminosity of $\sim6\times10^{40}$ s$^{-1}$ would produce a photoevaporation rate of only $\dot{M}_{\rm w,EUV}\simeq1.3\times10^{-10}~M_{\odot}$ yr$^{-1}$. Therefore, other mechanisms are required to disperse the disk in the timescale imposed by observations.

We measure the energy emitted by extensive air showers in the form of radio emission in the frequency range from 30 to 80 MHz. Exploiting the accurate energy scale of the Pierre Auger Observatory, we obtain a radiation energy of 15.8 \pm 0.7 (stat) \pm 6.7 (sys) MeV for cosmic rays with an energy of 1 EeV arriving perpendicularly to a geomagnetic field of 0.24 G, scaling quadratically with the cosmic-ray energy. A comparison with predictions from state-of-the-art first-principle calculations shows agreement with our measurement. The radiation energy provides direct access to the calorimetric energy in the electromagnetic cascade of extensive air showers. Comparison with our result thus allows the direct calibration of any cosmic-ray radio detector against the well-established energy scale of the Pierre Auger Observatory.

The Auger Engineering Radio Array (AERA) is part of the Pierre Auger Observatory and is used to detect the radio emission of cosmic-ray air showers. These observations are compared to the data of the surface detector stations of the Observatory, which provide well-calibrated information on the cosmic-ray energies and arrival directions. The response of the radio stations in the 30 to 80 MHz regime has been thoroughly calibrated to enable the reconstruction of the incoming electric field. For the latter, the energy deposit per area is determined from the radio pulses at each observer position and is interpolated using a two-dimensional function that takes into account signal asymmetries due to interference between the geomagnetic and charge-excess emission components. The spatial integral over the signal distribution gives a direct measurement of the energy transferred from the primary cosmic ray into radio emission in the AERA frequency range. We measure 15.8 MeV of radiation energy for a 1 EeV air shower arriving perpendicularly to the geomagnetic field. This radiation energy -- corrected for geometrical effects -- is used as a cosmic-ray energy estimator. Performing an absolute energy calibration against the surface-detector information, we observe that this radio-energy estimator scales quadratically with the cosmic-ray energy as expected for coherent emission. We find an energy resolution of the radio reconstruction of 22% for the data set and 17% for a high-quality subset containing only events with at least five radio stations with signal.

We report the discovery of a dwarf protoplanetary disk around the star XZ Tau B that shows all the features of a classical transitional disk but on a much smaller scale. The disk has been imaged with the Atacama Large Millimeter/Submillimeter Array (ALMA), revealing that its dust emission has a quite small radius of ~ 3.4 au and presents a central cavity of ~ 1.3 au in radius that we attribute to clearing by a compact system of orbiting (proto)planets. Given the very small radii involved, evolution is expected to be much faster in this disk (observable changes in a few months) than in classical disks (observable changes requiring decades) and easy to monitor with observations in the near future. From our modeling we estimate that the mass of the disk is large enough to form a compact planetary system.

The azimuthal asymmetry in the risetime of signals in Auger surface detector stations is a source of information on shower development. The azimuthal asymmetry is due to a combination of the longitudinal evolution of the shower and geometrical effects related to the angles of incidence of the particles into the detectors. The magnitude of the effect depends upon the zenith angle and state of development of the shower and thus provides a novel observable, $(\sec \theta)_\mathrm{max}$, sensitive to the mass composition of cosmic rays above $3 \times 10^{18}$ eV. By comparing measurements with predictions from shower simulations, we find for both of our adopted models of hadronic physics (QGSJETII-04 and EPOS-LHC) an indication that the mean cosmic-ray mass increases slowly with energy, as has been inferred from other studies. However, the mass estimates are dependent on the shower model and on the range of distance from the shower core selected. Thus the method has uncovered further deficiencies in our understanding of shower modelling that must be resolved before the mass composition can be inferred from $(\sec \theta)_\mathrm{max}$.