Atacama Large Millimetre/sub-millimetre Array (ALMA) observations of the CO(1-0) and CO(2-1) emissions from the circumstellar envelope of the Asymptotic Giant Branch (AGB) star EP Aqr have been made with four times better spatial resolution than previously available. They are analysed with emphasis on the de-projection in space of the effective emissivity and flux of matter using as input a prescribed configuration of the velocity field, assumed to be radial. The data are found to display an intrinsic axi-symmetry with respect to an axis making a small angle with the line of sight. A broad range of wind configurations, from prolate (bipolar) to oblate (equatorial) has been studied and found to be accompanied by significant equatorial emission. Qualitatively, the effective emissivity is enhanced near the equator to produce the central narrow component observed in the Doppler velocity spectra and its dependence on star latitude generally follows that of the wind velocity with the exception of an omni-present depression near the poles. In particular, large equatorial expansion velocities produce a flared disc or a ring of effective emissivity and mass loss. The effect on the determination of the orientation of the star axis of radial velocity gradients and possibly competing rotation and expansion in the equatorial disc is discussed. In general, the flux of matter is found to reach a broad maximum at distances of the order of 500 au from the star. Arguments are given that may be used to prefer one wind velocity distribution to another. As a result of the improved quality of the data, a deeper understanding of the constraints imposed on morphology and kinematics has been obtained.

Observations by the Atacama Large Millimeter/sub-millimeter Array of the dust continuum and $^{13}$CO(3-2) millimeter emissions of the triple stellar system GG Tau A are analysed, giving evidence for a rotating gas disc and a concentric and coplanar dust ring. The present work complements an earlier analysis (Tang et al. 2016) by exploring detailed properties of the gas disc. A 95% confidence level upper limit of 0.24 arcsec (34 au) is placed on the disc scale height at a distance of 1 arcsec (140 au) from the central stars. Evidence for Keplerian rotation of the gas disc is presented, the rotation velocity reaching ~3.1 kms$^{-1}$ at 1 arcsec from the central stars, and a 99% confidence level upper limit of 9% is placed on a possible in-fall velocity relative contribution. Variations of the intensity across the disc area are studied in detail and confirm the presence of a hot spot in the south-eastern quadrant. However several other significant intensity variations, in particular a depression in the northern direction, are also revealed. Variations of the intensity are found to be positively correlated to variations of the line width. Possible contributions to the measured line width are reviewed, suggesting an increase of the disc temperature and opacity with decreasing distance from the stars.

We present high resolution ALMA observations of the CO(3-2) and 350 GHz continuum emissions of the debris disc of 49 Ceti, known to be particularly rich in molecular gas in spite of its age. The main new results are: i) both CO and dust discs share a same position angle and a same inclination but the gas disc is more homogeneous, more central and thinner than the dust disc; ii) evidence is obtained for a significant deficit of observed CO(3-2) emission at Doppler velocities differing from the star systemic velocity by less than 1 \kms; iii) gas velocities are accurately measured and found Keplerian over a broad range of disc radii; iv) the observed CO(3-2) line width is dominated by Keplerian shear and upper limits are obtained to the intrinsic line width. Simple phenomenological models of both CO(3-2) and \mbox{350 GHz} continuum emissions are presented, requiring the use of only very few parameters. The results are discussed in the frame of currently favoured models.

Observations by the Atacama Large Millimetre/sub-millimetre Array of the 358 GHz continuum emission of the gravitationally lensed quasar host RX J0911.4+0551 have been analysed. They complement earlier Plateau de Bure Interferometer observations of the CO(7-6) emission. The good knowledge of the lensing potential obtained from Hubble Space Telescope observations of the quasar makes a joint analysis of the three emissions possible. It gives evidence for the quasar source to be concentric with the continuum source within 0.31 kpc and with the CO(7-6) source within 1.10 kpc. It also provides a measurement of the size of the continuum source, 0.76 $\pm$ 0.04 kpc FWHM, making RX J0911.4+0551 one of the few high redshift galaxies for which the dust and gas components are resolved with dimensions being measured. Both are found to be very compact, the former being smaller than the latter by a factor of $\sim$3.4$\pm$0.4. Moreover, new measurements of the CO ladder $-$ CO(10-9) and CO(11-10) $-$ are presented that confirm the extreme narrowness of the CO line width (107$\pm$20 km s$^{-1}$ on average). Their mere detection implies higher temperature and/or density than for typical quasar hosts at this redshift and suggests a possible contribution of the central AGN to gas and dust heating. The results are interpreted in terms of current understanding of galaxy evolution at the peak of star formation. They suggest that RX J0911.4+0551 is a young galaxy in an early stage of its evolution, having experienced no recent major mergers, star formation being concentrated in its centre.

Using ALMA observations of the C$^{18}$O(2-1) line emission of the gas envelope of protostar L1527, we have reconstructed its morphology and kinematics under the assumption of axisymmetry about the west-east axis. The main original contribution to our understanding of the formation process of L1527 is the presentation of a simple 3D parameterisation based solely on regions that are not dominated by absorption. In the explored range ($\sim$0.7 to 5 arcsec from the star) the model reproduces observations better than earlier attempts. The main results include: a measurement of the rotation velocity that confirms its evolution to Keplerian toward short distances; a measurement of the mean in-fall velocity, 0.43$\pm$0.10 kms$^{-1}$, lower than free fall velocity, with no evidence for the significant $r$-dependence suggested by an earlier analysis; a measurement of the central mass, 0.23$\pm$0.06 M$_{\odot}$ within a distance of 1.5 arcsec from the star, in agreement with earlier estimates obtained from a different range of distances; evidence for a strong disc plane depression of the in-falling flux resulting in an $X$ shaped flow possibly caused by the freeze-out of CO molecules on dust grains; a measurement of the accretion rate, 3.5$\pm$1.0 10$^{-7}$ M$_{\odot}$yr$^{-1}$ at a distance of 1 arcsec (140 au) from the star; evidence for a 10$^\circ$ tilt of the symmetry plane of the envelope about the line of sight, cancelling below $\sim$3 arcsec from the star, but matching infrared observations and being also apparent on the sky map of the mean Doppler velocity.

The unprecedented quality of the observations available from the Atacama Large Millimetre/sub-millimetre Array (ALMA) calls for analysis methods making the best of them. Reconstructing in space the morphology and kinematics of radio sources is an underdetermined problem that requires imposing additional constraints for its solution. The hypothesis of rotational invariance about a well-defined star axis, which is a good approximation to the description of the gas envelopes of many evolved stars and protostars, is particularly efficient in this role. In the first part of the article, a systematic use of simulated observations allows for identifying the main problems and for constructing quantities aimed at solving them. In particular the evaluation of the orientation of the star axis in space and the differentiation between expansion along the star axis and rotation about it are given special attention. The use of polar rather than Cartesian sky coordinates is shown to better match the morphology and kinematics of actual stars. The radial dependence of the gas density and temperature and the possible presence of velocity gradients are briefly considered. In the second part, the results obtained in the first part are applied to a few stars taken as examples with the aim of evaluating their usefulness when applied to concrete cases. A third part takes stock of what precedes and formulates some guidelines for modelling the radio emission of axisymmetric radio sources, limited however to the mathematics and geometry of the problem, physics considerations being ignored.

Observations of 12CO(3-2) emission of the circumbinary envelope of Mira Ceti, made by ALMA are analysed. The observed Doppler velocity distribution is made of three components: a blue-shifted south-eastern arc, which can be described as a ring in slow radial expansion, ~1.7 km/s, making an angle of ~50 deg with the plane of the sky and born some 2000 years ago; a few arcs, probably born at the same epoch as the blue-shifted arc, all sharing Doppler velocities red-shifted by approximately 3 +/- 2 km/s with respect to the main star; the third, central region dominated by the circumbinary envelope, displaying two outflows in the south-western and north-eastern hemispheres. At short distances from the star, up to ~1.5", these hemispheres display very different morphologies: the south-western outflow covers a broad solid angle, expands radially at a rate between 5 and 10 km/s and is slightly red shifted; the north-eastern outflow consists of two arms, both blue-shifted, bracketing a broad dark region where emission is suppressed. At distances between ~1.5" and ~2.5" the asymmetry between the two hemispheres is significantly smaller and detached arcs, particularly spectacular in the north-eastern hemisphere are present. Close to the stars, we observe a mass of gas surrounding Mira B, with a size of a few tens of AU, and having Doppler velocities with respect to Mira B reaching +/-1.5 km/s, which we interpret as gas flowing from Mira A toward Mira B.

Observations of the $^{12}$CO(3-2) emission of the circumstellar envelope (CSE) of the variable star $\pi^1$ Gru using the compact array (ACA) of the ALMA observatory have been recently made accessible to the public. An analysis of the morphology and kinematics of the CSE is presented with a result very similar to that obtained earlier for $^{12}$CO(2-1) emission by Chiu et al. (2006) using the Sub-Millimeter Array. A quantitative comparison is made using their flared disk model. A new model is presented that provides a significantly better description of the data, using radial winds and smooth evolutions of the radio emission and wind velocity from the stellar equator to the poles.

We analyse ALMA observations of the 12CO(3-2) emission of the circumstellar envelope (CSE) of the Mira variable binary star W Aql. These provide, for the first time, spatially resolved Doppler velocity spectra of the CSE up to angular distances to the central star of ~ 5" (meaning some 2000 AU). The exploratory nature of the observations (only five minutes in each of two different configurations) does not allow for a detailed modelling of the properties of the CSE but provides important qualitative information on its morphology and kinematics. Emission is found to be enhanced along an axis moving from east/west to north-east/south-west when the angular distance from the central star projected on the plane of the sky increases from zero to four arcseconds. In parallel, the Doppler velocity distribution displays asymmetry along an axis moving from east/west to north-west/south-east. The results are discussed in the context of earlier observations, in particular of the dust morphology.

Studies of the CO and HI radio emission of some evolved stars are presented using data collected by the IRAM Plateau de Bure interferometer and Pico Veleta telescope, the Nan\c{c}ay Radio Telescope and the JVLA and ALMA arrays. Approximate axial symmetry of the physical and kinematic properties of the circumstellar envelope (CSE) are observed in CO emission, in particular, from RS Cnc, EP Aqr and the Red Rectangle. A common feature is the presence of a bipolar outflow causing an enhanced wind velocity in the polar directions. HI emission extends to larger radial distances than probed by CO emission and displays features related to the interaction between the stellar outflow and interstellar matter. With its unprecedented sensitivity, FAST will open a new window on such studies. Its potential in this domain is briefly illustrated.

The CO(1-0) and (2-1) emission of the circumstellar envelope of the AGB star EP Aqr has been observed using the IRAM PdBI and the IRAM 30-m telescope. The line profiles reveal the presence of two distinct components centered on the star velocity, a broad component extending up to ~10 km/s and a narrow component indicating an expansion velocity of ~2 km/s. An early analysis of these data was performed under the assumption of isotropic winds. The present study revisits this interpretation by assuming instead a bipolar outflow nearly aligned with the line of sight. A satisfactory description of the observed flux densities is obtained with a radial expansion velocity increasing from ~2 km/s at the equator to ~10 km/s near the poles. The angular aperture of the bipolar outflow is ~45 deg with respect to the star axis, which makes an angle of ~13 deg with the line of sight. A detailed study of the CO(1-0) to CO(2-1) flux ratio reveals a significant dependence of the temperature on the star latitude, smaller and steeper at the poles than at the equator at large distances from the star. Under the hypothesis of radial expansion and of rotation invariance about the star axis, the effective density has been evaluated in space as a function of star coordinates. Evidence is found for an enhancement of the effective density in the northern hemisphere of the star at angular distances in excess of ~3" and covering the whole longitudinal range. The peak velocity of the narrow component is observed to vary slightly with position on the sky, a variation consistent with the model and understood as the effect of the inclination of the star axis with respect to the line of sight. While the phenomenological model presented here reproduces well the general features of the observations, significant differences are also revealed, which would require a better spatial resolution to be properly described.

Observations of an unprecedented quality made by ALMA on the Red Rectangle of CO(3-2) and CO(6-5) emissions are analysed jointly with the aim of obtaining as simple as possible a description of the gas morphology and kinematics. Evidence is found for polar conical outflows and for a broad equatorial torus in rotation and expansion. Simple models of both are proposed. Comparing CO(6-5) and CO(3-2) emissions provides evidence for a strong temperature enhancement over the polar outflows. Continuum emission (dust) is seen to be enhanced in the equatorial region. Observed asymmetries are briefly discussed.

We describe the method devised to reconstruct inclined cosmic-ray air showers with zenith angles greater than $60^\circ$ detected with the surface array of the Pierre Auger Observatory. The measured signals at the ground level are fitted to muon density distributions predicted with atmospheric cascade models to obtain the relative shower size as an overall normalization parameter. The method is evaluated using simulated showers to test its performance. The energy of the cosmic rays is calibrated using a sub-sample of events reconstructed with both the fluorescence and surface array techniques. The reconstruction method described here provides the basis of complementary analyses including an independent measurement of the energy spectrum of ultra-high energy cosmic rays using very inclined events collected by the Pierre Auger Observatory.

A flux of neutrons from an astrophysical source in the Galaxy can be detected in the Pierre Auger Observatory as an excess of cosmic-ray air showers arriving from the direction of the source. To avoid the statistical penalty for making many trials, classes of objects are tested in combinations as nine "target sets", in addition to the search for a neutron flux from the Galactic Center or from the Galactic Plane. Within a target set, each candidate source is weighted in proportion to its electromagnetic flux, its exposure to the Auger Observatory, and its flux attenuation factor due to neutron decay. These searches do not find evidence for a neutron flux from any class of candidate sources. Tabulated results give the combined p-value for each class, with and without the weights, and also the flux upper limit for the most significant candidate source within each class. These limits on fluxes of neutrons significantly constrain models of EeV proton emission from non-transient discrete sources in the Galaxy.

Measurements of air showers made using the hybrid technique developed with the fluorescence and surface detectors of the Pierre Auger Observatory allow a sensitive search for point sources of EeV photons anywhere in the exposed sky. A multivariate analysis reduces the background of hadronic cosmic rays. The search is sensitive to a declination band from -85{\deg} to +20{\deg}, in an energy range from 10^17.3 eV to 10^18.5 eV. No photon point source has been detected. An upper limit on the photon flux has been derived for every direction. The mean value of the energy flux limit that results from this, assuming a photon spectral index of -2, is 0.06 eV cm^-2 s^-1, and no celestial direction exceeds 0.25 eV cm^-2 s^-1. These upper limits constrain scenarios in which EeV cosmic ray protons are emitted by non-transient sources in the Galaxy.

A multipath mechanism similar to that used in Australia sixty years ago by the Sea-cliff Interferometer is shown to generate correlations between the periods of oscillations observed by two distant radio telescopes pointed to the Sun. The oscillations are the result of interferences between the direct wave detected in the main antenna lobe and its reflection on ground detected in a side lobe. A model is made of such oscillations in the case of two observatories located at equal longitudes and opposite tropical latitudes, respectively in Ha Noi (Viet Nam) and Learmonth (Australia), where similar radio telescopes are operated at 1.4 GHz. Simple specular reflection from ground is found to give a good description of the observed oscillations and to explain correlations that had been previously observed and for which no satisfactory interpretation, instrumental or other, had been found.

The Pierre Auger Observatory is the world's largest cosmic ray observatory. Our current exposure reaches nearly 40,000 km$^2$ str and provides us with an unprecedented quality data set. The performance and stability of the detectors and their enhancements are described. Data analyses have led to a number of major breakthroughs. Among these we discuss the energy spectrum and the searches for large-scale anisotropies. We present analyses of our X$_{max}$ data and show how it can be interpreted in terms of mass composition. We also describe some new analyses that extract mass sensitive parameters from the 100% duty cycle SD data. A coherent interpretation of all these recent results opens new directions. The consequences regarding the cosmic ray composition and the properties of UHECR sources are briefly discussed.

Joint contributions of the Pierre Auger and Telescope Array Collaborations to the 33rd International Cosmic Ray Conference, Rio de Janeiro, Brazil, July 2013: cross-calibration of the fluorescence telescopes, large scale anisotropies and mass composition.

Contributions of the Pierre Auger Collaboration to the 33rd International Cosmic Ray Conference, Rio de Janeiro, Brazil, July 2013

After a brief reminder of the mechanism of gravitational lensing of extended sources, the particular case of the host galaxy of QSO RXJ0911, a high redshift (z~2.8) quadruply imaged quasar, is explored. The non linearity of the problem, together with the proximity of the source to a cusp of the lens inner caustic, have important consequences on the dependence of the image appearance on the size and shape of the source. Their expected main features and their interpretation in terms of source extension and shape are investigated in a spirit of simplicity and in preparation for the analysis of high sensitivity and spatial resolution images that will soon be within reach with the completion of the Atacama Large Millimeter/submillimeter Array (ALMA). In particular, the information on source size carried by relative image brightness is discussed. Extension of the results to other types of quadruply imaged quasars is briefly considered.

The Pierre Auger Observatory is exploring the potential of the radio detection technique to study extensive air showers induced by ultra-high energy cosmic rays. The Auger Engineering Radio Array (AERA) addresses both technological and scientific aspects of the radio technique. A first phase of AERA has been operating since September 2010 with detector stations observing radio signals at frequencies between 30 and 80 MHz. In this paper we present comparative studies to identify and optimize the antenna design for the final configuration of AERA consisting of 160 individual radio detector stations. The transient nature of the air shower signal requires a detailed description of the antenna sensor. As the ultra-wideband reception of pulses is not widely discussed in antenna literature, we review the relevant antenna characteristics and enhance theoretical considerations towards the impulse response of antennas including polarization effects and multiple signal reflections. On the basis of the vector effective length we study the transient response characteristics of three candidate antennas in the time domain. Observing the variation of the continuous galactic background intensity we rank the antennas with respect to the noise level added to the galactic signal.

The Pierre Auger Observatory is a facility built to detect air showers produced by cosmic rays above 10^17 eV. During clear nights with a low illuminated moon fraction, the UV fluorescence light produced by air showers is recorded by optical telescopes at the Observatory. To correct the observations for variations in atmospheric conditions, atmospheric monitoring is performed at regular intervals ranging from several minutes (for cloud identification) to several hours (for aerosol conditions) to several days (for vertical profiles of temperature, pressure, and humidity). In 2009, the monitoring program was upgraded to allow for additional targeted measurements of atmospheric conditions shortly after the detection of air showers of special interest, e.g., showers produced by very high-energy cosmic rays or showers with atypical longitudinal profiles. The former events are of particular importance for the determination of the energy scale of the Observatory, and the latter are characteristic of unusual air shower physics or exotic primary particle types. The purpose of targeted (or "rapid") monitoring is to improve the resolution of the atmospheric measurements for such events. In this paper, we report on the implementation of the rapid monitoring program and its current status. The rapid monitoring data have been analyzed and applied to the reconstruction of air showers of high interest, and indicate that the air fluorescence measurements affected by clouds and aerosols are effectively corrected using measurements from the regular atmospheric monitoring program. We find that the rapid monitoring program has potential for supporting dedicated physics analyses beyond the standard event reconstruction.

The development of extensive air showers at extreme energies is studied using a simulation model much simpler and cruder, but also more transparent and flexible, than existing sophisticated codes. Evidence for its satisfactory performance is presented. As an illustration, shower elongation rates are evaluated in the $10^{18}$ to $10^{20}$ eV region and compared with recently published data. Lateral distribution functions of both muons and electrons/photons are also briefly discussed. Reliable results are obtained in the comparison between proton-induced and iron-induced showers.

The Surface Detector of the Pierre Auger Observatory is sensitive to neutrinos of all flavours above 0.1 EeV. These interact through charged and neutral currents in the atmosphere giving rise to extensive air showers. When interacting deeply in the atmosphere at nearly horizontal incidence, neutrinos can be distinguished from regular hadronic cosmic rays by the broad time structure of their shower signals in the water-Cherenkov detectors. In this paper we present for the first time an analysis based on down-going neutrinos. We describe the search procedure, the possible sources of background, the method to compute the exposure and the associated systematic uncertainties. No candidate neutrinos have been found in data collected from 1 January 2004 to 31 May 2010. Assuming an E^-2 differential energy spectrum the limit on the single flavour neutrino is (E^2 * dN/dE) < 1.74x10^-7 GeV cm^-2 s^-1 sr^-1 at 90% C.L. in the energy range 1x10^17 eV < E < 1x10^20 eV.

Atmospheric conditions at the site of a cosmic ray observatory must be known for reconstructing observed extensive air showers. The Global Data Assimilation System (GDAS) is a global atmospheric model predicated on meteorological measurements and numerical weather predictions. GDAS provides altitude-dependent profiles of the main state variables of the atmosphere like temperature, pressure, and humidity. The original data and their application to the air shower reconstruction of the Pierre Auger Observatory are described. By comparisons with radiosonde and weather station measurements obtained on-site in Malarg\"ue and averaged monthly models, the utility of the GDAS data is shown.

We present a comprehensive study of the influence of the geomagnetic field on the energy estimation of extensive air showers with a zenith angle smaller than $60^\circ$, detected at the Pierre Auger Observatory. The geomagnetic field induces an azimuthal modulation of the estimated energy of cosmic rays up to the ~2% level at large zenith angles. We present a method to account for this modulation of the reconstructed energy. We analyse the effect of the modulation on large scale anisotropy searches in the arrival direction distributions of cosmic rays. At a given energy, the geomagnetic effect is shown to induce a pseudo-dipolar pattern at the percent level in the declination distribution that needs to be accounted for.

The surface detector array of the Pierre Auger Observatory consists of 1600 water-Cherenkov detectors, for the study of extensive air showers (EAS) generated by ultra-high-energy cosmic rays. We describe the trigger hierarchy, from the identification of candidate showers at the level of a single detector, amongst a large background (mainly random single cosmic ray muons), up to the selection of real events and the rejection of random coincidences. Such trigger makes the surface detector array fully efficient for the detection of EAS with energy above $3\times 10^{18}$ eV, for all zenith angles between 0$^\circ$ and 60$^\circ$, independently of the position of the impact point and of the mass of the primary particle. In these range of energies and angles, the exposure of the surface array can be determined purely on the basis of the geometrical acceptance.

In this paper we introduce the concept of Lateral Trigger Probability (LTP) function, i.e., the probability for an extensive air shower (EAS) to trigger an individual detector of a ground based array as a function of distance to the shower axis, taking into account energy, mass and direction of the primary cosmic ray. We apply this concept to the surface array of the Pierre Auger Observatory consisting of a 1.5 km spaced grid of about 1600 water Cherenkov stations. Using Monte Carlo simulations of ultra-high energy showers the LTP functions are derived for energies in the range between 10^{17} and 10^{19} eV and zenith angles up to 65 degs. A parametrization combining a step function with an exponential is found to reproduce them very well in the considered range of energies and zenith angles. The LTP functions can also be obtained from data using events simultaneously observed by the fluorescence and the surface detector of the Pierre Auger Observatory (hybrid events). We validate the Monte-Carlo results showing how LTP functions from data are in good agreement with simulations.

We present the results of an analysis of data recorded at the Pierre Auger Observatory in which we search for groups of directionally-aligned events (or `multiplets') which exhibit a correlation between arrival direction and the inverse of the energy. These signatures are expected from sets of events coming from the same source after having been deflected by intervening coherent magnetic fields. The observation of several events from the same source would open the possibility to accurately reconstruct the position of the source and also measure the integral of the component of the magnetic field orthogonal to the trajectory of the cosmic rays. We describe the largest multiplets found and compute the probability that they appeared by chance from an isotropic distribution. We find no statistically significant evidence for the presence of multiplets arising from magnetic deflections in the present data.