We present new observations of the active star-formation region NGC 1333 in the Perseus molecular cloud complex from the James Clerk Maxwell Telescope (JCMT) B-Fields In Star-forming Region Observations (BISTRO) survey with the POL-2 instrument. The BISTRO data cover the entire NGC 1333 complex (~ 1.5 pc x 2 pc) at 0.02 pc resolution and spatially resolve the polarized emission from individual filamentary structures for the first time. The inferred magnetic field structure is complex as a whole, with each individual filament aligned at different position angles relative to the local field orientation. We combine the BISTRO data with low- and high- resolution data derived from Planck and interferometers to study the multi-scale magnetic field structure in this region. The magnetic field morphology drastically changes below a scale of ~ 1 pc and remains continuous from the scales of filaments (~ 0.1 pc) to that of protostellar envelopes (~ 0.005 pc or ~ 1000 au). Finally, we construct simple models in which we assume that the magnetic field is always perpendicular to the long axis of the filaments. We demonstrate that the observed variation of the relative orientation between the filament axes and the magnetic field angles are well reproduced by this model, taking into account the projection effects of the magnetic field and filaments relative to the plane of the sky. These projection effects may explain the apparent complexity of the magnetic field structure observed at the resolution of BISTRO data towards the filament network.

We study the morpho-kinematics of the circumbinary envelope of Mira Ceti between $\sim$100 and $\sim$350 au from the stars using ALMA observations of the SiO ($\nu$=0, $J$=5-4) and CO ($\nu$=0, $J$=3-2) emissions with the aim of presenting an accurate and reliable picture of what cannot be ignored when modelling the dynamics at stake. A critical study of the uncertainties attached to imaging is presented. The line emissions are shown to be composed of a few separated fragments. They are described in detail and plausible interpretations of their genesis are discussed. Evidence for a focusing effect of the Mira A wind by Mira B over the past century is presented; it accounts for only a small fraction of the overall observed emission but its accumulation over several orbital periods may have produced an enhancement of CO emission in the orbital plane of Mira B. We identify a South-western outflow and give arguments for the anti-correlation observed between CO and SiO emissions being the result of a recent mass ejection accompanied by a shock wave. We discuss the failure of simple scenarios that have been proposed earlier to explain some of the observed features and comment on the apparent lack of continuity between the present observations and those obtained in the close environment of the stars. Evidence is obtained for the presence of large Doppler velocity components near the line of sight aiming to the star, possibly revealing the presence of important turbulence at $\sim$5 to 10 au away from Mira A.

We study the morpho-kinematics in the nascent wind of AGB star R Doradus in the light of high Doppler velocity wings observed in the spectral lines of several species. We probe distances from the star between ~10 and ~100 au using ALMA observations of the emission of five different molecular lines. High Doppler velocity enhancements of the line emission are observed in the vicinity of the line of sight crossing the star, reminiscent of those recently interpreted as gas streams in the nascent wind of a similar AGB star, EP Aqr. They are present in both blue-shifted and red-shifted hemispheres but are not exactly back-to-back. They are accelerated at a typical rate of 0.7 km s$^{-1}$ au$^{-1}$ up to some 20 km s$^{-1}$. Important differences are observed between the emissions of different molecules. We exclude an effect of improper continuum subtraction. However, in contrast to EP Aqr, the line of sight plays no particular role in the R Dor morpho-kinematics, shedding doubt on the validity of a gas stream interpretation.We discuss possible interpretations in terms of stellar pulsations or of rotation of the gas in the environment of the star. We conclude that, in the state of current knowledge, no fully convincing picture of the physics governing the production of such high velocities, typically twice as large as the terminal velocity, can be reliably drawn. New high resolution analyses of observations of the nascent wind of oxygen-rich AGB stars are needed to clarify the issue.

We used new ALMA $^{13}$CO and C$^{18}$O(3-2) observations obtained at high angular resolution ($\sim$0.2") together with previous CO(3-2) and (6-5) ALMA data and continuum maps at 1.3 and 0.8 mm in order to determine the gas properties (temperature, density, and kinematics) in the cavity and to a lesser extent in the outer disk of GG Tau A, the prototype of a young triple T Tauri star that is surrounded by a massive and extended Keplerian outer disk. By deprojecting, we studied the radial and azimuthal gas distribution and its kinematics. We also applied a new method to improve the deconvolution of the CO data and in particular better quantify the emission from gas inside the cavity. We perform local and nonlocal thermodynamic equilibrium studies in order to determine the excitation conditions and relevant physical parameters inside the ring and in the central cavity. Residual emission after removing a smooth-disk model indicates unresolved structures at our angular resolution, probably in the form of irregular rings or spirals. The outer disk is cold, with a temperature $<20$ K beyond 250 au that drops quickly (r$^{-1}$). The kinematics of the gas inside the cavity reveals infall motions at about 10% of the Keplerian speed. We derive the amount of gas in the cavity, and find that the brightest clumps, which contain about 10% of this mass, have kinetic temperatures 40$-$80 K, CO column densities of a few 10$^{17}$ cm$^{-2}$, and H$_2$ densities around 10$^7$ cm$^{-3}$. Although the gas in the cavity is only a small fraction of the disk mass, the mass accretion rate throughout the cavity is comparable to or higher than the stellar accretion rate. It is accordingly sufficient to sustain the circumstellar disks on a long timescale.

We analyse MESSENGER reflectance measurements covering the northern polar region of Mercury, the least studied region of the northern mercurian hemisphere. We use observations from the Mercury Dual Imaging System Wide-Angle Camera (MDIS/WAC) and the Mercury Atmospheric and Surface Composition Spectrometer (MASCS/VIRS) to study the spectral dependence of the surface reflectance. The results obtained from the observations made by both instruments are remarkably consistent. We find that a second degree polynomial description of the measured reflectance spectra gives very good fits to the data and that the information that they carry can best be characterized by two parameters, the mean reflectance and the mean relative spectral slope, averaged over the explored range of wavelengths. The properties of the four main types of terrains known to form Mercury's regolith in the northern region, smooth plains (SP), heavily cratered terrain (HCT), fresh ejecta/materials and red pitted ground (RPG) are examined in terms of these two parameters. The results are compared, and found consistent with those obtained by earlier studies in spite of difficulties met in obtaining accurate reflectance measurements under the large incidence angle condition characteristic of polar regions. These results will help with the preparation of the BepiColombo mission and with supporting its observational strategy.

We analyse ALMA observations of the SO($J_K=6_5-5_4$) emission of the circumstellar envelope of oxygen-rich AGB star R Dor, probing distances between 20 and 100 au from the star where the nascent wind is building up. We give evidence for the slow wind to host, in addition to a previously observed rotating disc, a radial outflow covering very large solid angles and displaying strong inhomogeneity both in direction and radially: the former takes the form of multiple cores and the latter displays a radial dependence suggesting an episode of enhanced mass loss having occurred a century or so ago.

HD 163296 is one of the few protoplanetary discs displaying rings in the dust component. The present work uses ALMA observations of the 0.9 mm continuum emission having significantly better spatial resolution (~8 au) than previously available, providing new insight on the morphology of the dust disc and its double ring structure. The disc is shown to be thin and its position angle and inclination with respect to the sky plane are accurately measured as are the locations and shapes that characterize the observed ring/gap structure. Significant modulation of the intensity of the outer ring emission have been revealed and discussed. In addition, earlier ALMA observations of the emission of three molecular lines, CO(2-1), C18O(2-1), and DCO+(3-2), having a resolution of ~70 au, are used to demonstrate the Keplerian motion of the gas, found consistent with a central mass of 2.3 solar masses. An upper limit of ~9% of the rotation velocity is placed on the in-fall velocity. The beam size is shown to give the dominant contribution to the line widths, accounting for both their absolute values and their dependence on the distance to the central star.

Using ALMA observations of $^{12}$CO(2-1), $^{28}$SiO(5-4) and $^{32}$SO$_2$(16$_{6,10}$-17$_{5,13}$) emissions of the circumstellar envelope of AGB star EP Aqr, we describe the morpho-kinematics governing the nascent wind. Main results are: 1) Two narrow polar structures, referred to as jets, launched from less than 25 au away from the star, build up between $\sim$ 20 au and $\sim$ 100 au to a velocity of $\sim$ 20 \kms. They fade away at larger distances and are barely visible in CO data. 2) SO$_2$, SiO and CO emissions explore radial ranges reaching respectively $\sim$30 au, 250 au and 1000 au from the star, preventing the jets to be detected in SO$_2$ data. 3) Close to the star photosphere, rotation (undetected in SiO and CO data) and isotropic radial expansion combine with probable turbulence to produce a broad SO$_2$ line profile ($\sim$ 7.5 \kms\ FWHM). 4) A same axis serves as axis of rotation close to the star, as jet axis and as axi-symmetry axis at large distances. 5) A radial wind builds up at distances up to $\sim$ 300 au from the star, with larger velocity near polar than equatorial latitudes. 6) A sharp depletion of SiO and CO emissions, starting near the star, rapidly broadens to cover the whole blue-western quadrant, introducing important asymmetry in the CO and particularly SiO observations. 7) The $^{12}$C/$^{13}$C abundance ratio is measured as 9$\pm$2. 8) Plausible interpretations are discussed, in particular assuming the presence of a companion.

ALMA observations of CO(1-0) and CO(2-1) emissions of the circumstellar envelope of EP Aqr, an oxygen-rich AGB star, are reported. A thorough analysis of their properties is presented using an original method based on the separation of the data-cube into a low velocity component associated with an equatorial outflow and a faster component associated with a bipolar outflow. A number of important and new results are obtained concerning the distribution in space of the effective emissivity, the temperature, the density and the flux of matter. A mass loss rate of (1.6$\pm$0.4)10$^{-7}$ solar masses per year is measured. The main parameters defining the morphology and kinematics of the envelope are evaluated and uncertainties inherent to de-projection are critically discussed. Detailed properties of the equatorial region of the envelope are presented including a measurement of the line width and a precise description of the observed inhomogeneity of both morphology and kinematics. In particular, in addition to the presence of a previously observed spiral enhancement of the morphology at very small Doppler velocities, a similarly significant but uncorrelated circular enhancement of the expansion velocity is revealed, both close to the limit of sensitivity. The results of the analysis place significant constraints on the parameters of models proposing descriptions of the mass loss mechanism, but cannot choose among them with confidence.

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.

Studying molecular species in protoplanetary disks is very useful to characterize the properties of these objects, which are the site of planet formation. We attempt to constrain the chemistry of S-bearing molecules in the cold parts of circumstellar disk of GG Tau A. We searched for H$_2$S, CS, SO, and SO$_2$ in the dense disk around GG Tau A with the NOrthem Extended Millimeter Array (NOEMA) interferometer. We detected H$_2$S emission from the dense and cold ring orbiting around GG Tau A. This is the first detection of H$_2$S in a protoplanetary disk. We also detected HCO$^+$, H$^{13}$CO$^+$, and DCO$^+$ in the disk. Upper limits for other molecules, CCS, SO$_2$, SO, HC$_3$N, and $c$-C$_3$H$_2$ are also obtained. The observed DCO$^+$/HCO$^+$ ratio is similar to those in other disks. The observed column densities, derived using our radiative transfer code DiskFit, are then compared with those from our chemical code Nautilus. The column densities are in reasonable agreement for DCO$^{+}$, CS, CCS, and SO$_2$. For H$_2$S and SO, our predicted vertical integrated column densities are more than a factor of 10 higher than the measured values. Our results reinforce the hypothesis that only a strong sulfur depletion may explain the low observed H$_2$S column density in the disk. The H$_2$S detection in GG Tau A is most likely linked to the much larger mass of this disk compared to that in other T Tauri systems.

The problem of de-projection of radio line observations of axi-symmetric expanding circumstellar envelopes is studied with the aim of easing their analysis in terms of physics models. The arguments developed rest on the remark that, in principle, when the wind velocity distribution is known, the effective emissivity can be calculated at any point in space. The paper provides a detailed study of how much this is true in practice. The wind velocity distribution assumed to be axi-symmetric and in expansion, is described by four parameters: the angles defining the orientation of the symmetry axis, an overall velocity scale and a parameter measuring the elongation (prolateness) of the distribution. Tools are developed that allow for measuring, or at least constraining, each of the four parameters. The use of effective emissivity as relevant quantity, rather than temperature and density being considered separately, implies important assumptions and simplifications meaning that the approach being considered here is only a preliminary to, and by no means a replacement for, a physics analysis accounting for radiative transfer and hydrodynamics arguments. While most considerations are developed using simulated observations as examples, two case studies (EP Aqr, observed with ALMA, and RS Cnc, with NOEMA) are presented that illustrate their usefulness in practical cases.

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.