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 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.

Stars are bad neighbors: they often disturb their surroundings. They sometimes travel very fast through the interstellar medium (ISM). They frequently undergo violent ejection events which leave an imprint on their neighborhood (jets, winds, supernovae). These supersonic flows generate shocks both in the ejected material (bow shock) and in the stellar environment (termination shock). The global purpose of this thesis is to link the properties of stars and the ISM by studying these shocks. We model them with the Paris-Durham planar shock code, which incorporates a wealth of micro-physics and chemical processes relevant to the magnetized ISM. In the first part, we model 3D magnetized axis-symmetric bow shock by a statistical distribution of 1D planar shocks computed with the Paris-Durham model. For the first time, we examine systematically the effect of the geometry, age, and various other parameters on the H2 excitation diagram and emission line profiles. Then, we will show that this 3D bow shock model unprecedentedly improves interpretations of the H2 observations inOrion BN-KL and BHR71, and show how spectral resolved H2 line emission profiles on the Herbig-Haro objects HH54 can be used to extract a wealth of dynamical information. In the second part, we model 1D steady stellar winds from AGB stars, which trigger the termination shock, by including relevant physical and chemical ingredients in the Paris-Durham code. For the first time, we examine the time-dependent conversion of atomic to molecular hydrogen along the wind trajectories of "hot" and "cold" AGB stars. We suggest that the low abundance of HI inferred from observations is due to hydrogen locked in its molecular form. Then, we try to reproduce HI 21-cm line profile in a "hot" AGB (called Y CVn) and a "cold" AGB (called CW Leo).

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.

When a fast moving star or a protostellar jet hits an interstellar cloud, the surrounding gas gets heated and illuminated: a bow shock is born which delineates the wake of the impact. In such a process, the new molecules that are formed and excited in the gas phase become accessible to observations. In this article, we revisit models of H2 emission in these bow shocks. We approximate the bow shock by a statistical distribution of planar shocks computed with a magnetized shock model. We improve on previous works by considering arbitrary bow shapes, a finite irradiation field, and by including the age effect of non-stationary C-type shocks on the excitation diagram and line profiles of H2. We also examine the dependence of the line profiles on the shock velocity and on the viewing angle: we suggest that spectrally resolved observations may greatly help to probe the dynamics inside the bow shock. For reasonable bow shapes, our analysis shows that low velocity shocks largely contribute to H2 excitation diagram. This can result in an observational bias towards low velocities when planar shocks are used to interpret H2 emission from an unresolved bow. We also report a large magnetization bias when the velocity of the planar model is set independently. Our 3D models reproduce excitation diagrams in BHR71 and Orion bow shocks better than previous 1D models. Our 3D model is also able to reproduce the shape and width of the broad H2 1-0S(1) line profile in an Orion bow shock.

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.

Imaging studies with the VLA have revealed HI emission associated with the extended circumstellar shells of red giants. We analyse the spectral map obtained on Y CVn, a J-type carbon star on the AGB. The HI line profiles can be interpreted with a model of a detached shell resulting from the interaction of a stellar outflow with the local interstellar medium. We reproduce the spectral map by introducing a distortion along a direction corresponding to the star's motion in space. We then use this fitting to simulate observations expected from the FAST radiotelescope, and discuss its potential for improving ourdescription of the outer regions of circumstellar shells.

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.

We report on observations obtained at IRAM on two semi-regular variable Asymptotic Giant Branch (AGB) stars, RS Cnc and EP Aqr, undergoing mass loss at an intermediate rate of ~ 10^-7 solar mass per year. Interferometric data obtained with the Plateau-de-Bure interferometer (NOEMA) have been combined with On-The-Fly maps obtained with the 30-m telescope in the CO(1-0) and (2-1) rotational lines. The spectral maps of spatially resolved sources reveal an axisymmetric morphology in which matter is flowing out at a low velocity (~ 2 km/s) in the equatorial planes, and at a larger velocity (~ 8 km/s) along the polar axes. There are indications that this kind of morpho-kinematics is relatively frequent among stars at the beginning of their evolution on the Thermally-Pulsing AGB, in particular among those that show composite CO line profiles, and that it might be caused by the presence of a companion. We discuss the progress that could be expected for our understanding of the mass loss mechanisms in this kind of sources by increasing the spatial resolution of the observations with ALMA or NOEMA.

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 present HI line profiles for various models of circumstellar shells around red giants. In the calculations we take into account the effect of the background at 21 cm, and show that in some circumstances it may have an important effect on the shape and intensity of the observed line profiles. We show that self-absorption should also be considered depending on the mass loss rate and the temperature reached by circumstellar gas. HI emission from circumstellar shells has been mostly reported from stars with mass loss rates around 10$^{-7}$ solar masses per year. We discuss the possible reasons for the non detection of many sources with larger mass loss rates that are hallmarks of the end of the AGB phase. Although radiative transfer effects may weaken the line emission, they cannot alone account for this effect. Therefore, it seems likely that molecular hydrogen, rather than atomic hydrogen, dominates the composition of matter expelled by stars at the end of their evolution on the Asymptotic Giant Branch. However sensitive HI observations can still yield important information on the kinematics and physical properties of the circumstellar material at large distances from central stars with heavy mass loss, despite the low abundance of atomic hydrogen.

We present a phenomenological study of CO(1-0) and CO(2-1) emission from the circumstellar envelope (CSE) of the Asymptotic Giant Branch (AGB) star RS\,Cnc. It reveals departures from central symmetry that turn out to be efficient tools for the exploration of some of the CSE properties. We use a wind model including a bipolar flow with a typical wind velocity of $\sim$8 km\,s$^{-1}$ decreasing to $\sim$2 km\,s$^{-1}$ near the equator to describe Doppler velocity spectral maps obtained by merging data collected at the IRAM Plateau de Bure Interferometer and Pico Veleta single dish radio telescope. Parameters describing the wind morphology and kinematics are obtained, together with the radial dependence of the gas temperature in the domain of the circumstellar envelope probed by the CO observations. Significant north-south central asymmetries are revealed by the analysis, which we quantify using a simple phenomenological description. The origin of such asymmetries is unclear.

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.

We present a detailed study of the circumstellar gas distribution and kinematics of the semi-regular variable star RS Cnc on spatial scales ranging from ~1" (~150 AU) to ~6' (~0.25 pc). New modeling of CO1-0 and CO2-1 imaging observations leads to a revised characterization of RS Cnc's previously identified axisymmetric molecular outflow. Rather than a simple disk-outflow picture, we find that a gradient in velocity as a function of latitude is needed to fit the spatially resolved spectra, and in our preferred model, the density and the velocity vary smoothly from the equatorial plane to the polar axis. In terms of density, the source appears quasi-spherical, whereas in terms of velocity the source is axi-symmetric with a low expansion velocity in the equatorial plane and faster outflows in the polar directions. The flux of matter is also larger in the polar directions than in the equatorial plane. An implication of our model is that the stellar wind is still accelerated at radii larger than a few hundred AU, well beyond the radius where the terminal velocity is thought to be reached in an asymptotic giant branch star. The HI data show the previously detected head-tail morphology, but also supply additional detail about the atomic gas distribution and kinematics. We confirm that the `head' seen in HI is elongated in a direction consistent with the polar axis of the molecular outflow, suggesting that we are tracing an extension of the molecular outflow well beyond the molecular dissociation radius (up to ~0.05 pc). The 6'-long HI `tail' is oriented at a PA of 305{\deg}, consistent with the space motion of the star. We measure a total mass of atomic hydrogen ~0.0055 solar mass and estimate a lower limit to the timescale for the formation of the tail to be ~6.4x10^4 years. (abridged)

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

We report on high angular resolution observations of the CO(7-6) line and millimeter continuum in the host galaxy of the gravitationally lensed (z~2.8) quasar RXJ0911.4+0551 using the Plateau de Bure Interferometer. Our CO observations resolve the molecular disk of the source. Using a lens model based on HST observations we fit source models to the observed visibilities. We estimate a molecular disk radius of 1$\pm$0.2 kpc and an inclination of 69$\pm$6\deg, the continuum is more compact and is only marginally resolved by our observations. The relatively low molecular gas mass, $Mgas=(2.3\pm 0.5)\times 10^{9}$ Msolar, and far infrared luminosity, $LFIR=(7.2\pm 1.5) \times 10^{11}$ Lsolar, of this quasar could be explained by its relatively low dynamical mass, $Mdyn=(3.9\pm 0.9)\times 10^9$ Msolar. It would be a scaled-down version the QSOs usually found at high-z. The FIR and CO luminosities lie on the correlation found for QSOs from low to high redshifts and the gas-to-dust ratio ($45\pm 17$) is similar to the one measured in the z=6.4 QSO, SDSS J1148+5251. Differential magnification affects the continuum-to-line luminosity ratio, the line profile and possibly the spectral energy distribution.

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.

Studies of the cosmic ray energy spectrum at the highest energies with the Pierre Auger Observatory

Technical reports on operations and monitoring of the Pierre Auger Observatory

Astrophysical observations of ultra-high-energy cosmic rays with the Pierre Auger Observatory