EX Lup is a well-studied T Tauri star that represents the prototype of young eruptive stars EXors. In this paper we analyze new adaptive optics imaging and spectroscopic observations of EX Lup and its circumstellar environment in near-infrared in its quiescent phase. We observed EX Lup with the high contrast imager SPHERE/IRDIS in the dual-beam polarimetric imaging mode to resolve the circumstellar environment in near-infrared scattered light. We complemented these data with earlier SINFONI spectroscopy. We resolve for the first time in scattered light a compact feature around EX Lup azimuthally extending from 280deg to 360deg, and radially extending from 0.3arcsec to 0.55arcsec in the plane of the disk. We explore two different scenarios for the detected emission. One accounting for the emission as coming from the brightened walls of the cavity excavated by the outflow whose presence was suggested by ALMA observations in the J=3-2 line of 12CO. The other accounts for the emission as coming from an inclined disk. We detect for the first time a more extended circumstellar disk in scattered light, which shows that a region between 10 and 30 au is depleted of mum-size grains. We compare the J-, H- and K-band spectra obtained with SINFONI in quiescence with the spectra taken during the outburst, showing that all the emission lines were due to the episodic accretion event. Conclusions. Based on the morphology analysis we favour the scenario in which the scattered light is coming from a circumstellar disk rather than the outflow around EX Lup. We analyze the origin of the observed feature either as coming from a continuous circumstellar disk with a cavity, or from the illuminated wall of the outer disk or from a shadowed disk. Moreover, we discuss what is the origin of the mum-size grains depleted region, exploring the possibility that a sub-stellar companion may be the cause of it.

Shadows in scattered light images of protoplanetary disks are a common feature and support the presence of warps or misalignments between disk regions. These warps are possibly due to an inclined (sub-)stellar companion embedded in the disk. We study the morphology of the protoplanetary disk around the Herbig Ae star HD 139614 based on the first scattered light observations of this disk, which we model with the radiative transfer code MCMax3D. We obtained J- and H-band observations in polarized scattered light with VLT/SPHERE that show strong azimuthal asymmetries. In the outer disk, beyond ~30 au, a broad shadow spans a range of ~240{\deg} in position angle, in the East. A bright ring at ~16 au also shows an azimuthally asymmetric brightness, with the faintest side roughly coincidental with the brightest region of the outer disk. Additionally, two arcs are detected at ~34 au and ~50 au. We created a simple 4-zone approximation to a warped disk model of HD 139614 in order to qualitatively reproduce these features. The location and misalignment of the disk components were constrained from the shape and location of the shadows they cast. We find that the shadow on the outer disk covers a range of position angle too wide to be explained by a single inner misaligned component. Our model requires a minimum of two separate misaligned zones -- or a continuously warped region -- to cast this broad shadow on the outer disk. A small misalignment of ~4{\deg} between adjacent components can reproduce most of the observed shadow features. Multiple misaligned disk zones, potentially mimicing a warp, can explain the observed broad shadows in the HD 139614 disk. A planetary mass companion in the disk, located on an inclined orbit, could be responsible for such a feature and for the dust depleted gap responsible for a dip in the SED.

Observations of circumstellar environments to look for the direct signal of exoplanets and the scattered light from disks has significant instrumental implications. In the past 15 years, major developments in adaptive optics, coronagraphy, optical manufacturing, wavefront sensing and data processing, together with a consistent global system analysis have enabled a new generation of high-contrast imagers and spectrographs on large ground-based telescopes with much better performance. One of the most productive is the Spectro-Polarimetic High contrast imager for Exoplanets REsearch (SPHERE) designed and built for the ESO Very Large Telescope (VLT) in Chile. SPHERE includes an extreme adaptive optics system, a highly stable common path interface, several types of coronagraphs and three science instruments. Two of them, the Integral Field Spectrograph (IFS) and the Infra-Red Dual-band Imager and Spectrograph (IRDIS), are designed to efficiently cover the near-infrared (NIR) range in a single observation for efficient young planet search. The third one, ZIMPOL, is designed for visible (VIR) polarimetric observation to look for the reflected light of exoplanets and the light scattered by debris disks. This suite of three science instruments enables to study circumstellar environments at unprecedented angular resolution both in the visible and the near-infrared. In this work, we present the complete instrument and its on-sky performance after 4 years of operations at the VLT.

Context. Circumstellar disks and self-luminous giant exoplanets or companion brown dwarfs can be characterized through direct-imaging polarimetry at near-infrared wavelengths. SPHERE/IRDIS at the Very Large Telescope has the capabilities to perform such measurements, but uncalibrated instrumental polarization effects limit the attainable polarimetric accuracy. Aims. We aim to characterize and correct the instrumental polarization effects of the complete optical system, i.e. the telescope and SPHERE/IRDIS. Methods. We create a detailed Mueller matrix model in the broadband filters Y-, J-, H- and Ks, and calibrate it using measurements with SPHERE's internal light source and observations of two unpolarized stars. We develop a data-reduction method that uses the model to correct for the instrumental polarization effects, and apply it to observations of the circumstellar disk of T Cha. Results. The instrumental polarization is almost exclusively produced by the telescope and SPHERE's first mirror and varies with telescope altitude angle. The crosstalk primarily originates from the image derotator (K-mirror). At some orientations, the derotator causes severe loss of signal (>90% loss in H- and Ks-band) and strongly offsets the angle of linear polarization. With our correction method we reach in all filters a total polarimetric accuracy of <0.1% in the degree of linear polarization and an accuracy of a few degrees in angle of linear polarization. Conclusions. The correction method enables us to accurately measure the polarized intensity and angle of linear polarization of circumstellar disks, and is a vital tool for detecting unresolved (inner) disks and measuring the polarization of substellar companions. We have incorporated the correction method in a highly-automatic end-to-end data-reduction pipeline called IRDAP which is publicly available at https://irdap.readthedocs.io.

Planetary systems hold the imprint of the formation and of the evolution of planets especially at young ages, and in particular at the stage when the gas has dissipated leaving mostly secondary dust grains. The dynamical perturbation of planets in the dust distribution can be revealed with high-contrast imaging in a variety of structures. SPHERE, the high-contrast imaging device installed at the VLT, was designed to search for young giant planets in long period, but is also able to resolve fine details of planetary systems at the scale of astronomical units in the scattered-light regime. As a young and nearby star, NZ Lup was observed in the course of the SPHERE survey. A debris disk had been formerly identified with HST/NICMOS. We observed this system in the near-infrared with the camera in narrow and broad band filters and with the integral field spectrograph. High contrasts are achieved by the mean of pupil tracking combined with angular differential imaging algorithms. The high angular resolution provided by SPHERE allows us to reveal a new feature in the disk which is interpreted as a superimposition of two belts of planetesimals located at stellocentric distances of $\sim$85 and $\sim$115\,au, and with a mutual inclination of about 5$\degb$. Despite the very high inclination of the disk with respect to the line of sight, we conclude that the presence of a gap, that is, a void in the dust distribution between the belts, is likely. We discuss the implication of the existence of two belts and their relative inclination with respect to the presence of planets.

With an orbital distance comparable to that of Saturn in the solar system, \bpic b is the closest (semi-major axis $\simeq$\,9\,au) exoplanet that has been imaged to orbit a star. Thus it offers unique opportunities for detailed studies of its orbital, physical, and atmospheric properties, and of disk-planet interactions. With the exception of the discovery observations in 2003 with NaCo at the Very Large Telescope (VLT), all following astrometric measurements relative to \bpic have been obtained in the southwestern part of the orbit, which severely limits the determination of the planet's orbital parameters. We aimed at further constraining \bpic b orbital properties using more data, and, in particular, data taken in the northeastern part of the orbit. We used SPHERE at the VLT to precisely monitor the orbital motion of beta \bpic b since first light of the instrument in 2014. We were able to monitor the planet until November 2016, when its angular separation became too small (125 mas, i.e., 1.6\,au) and prevented further detection. We redetected \bpic b on the northeast side of the disk at a separation of 139\,mas and a PA of 30$^{\circ}$ in September 2018. The planetary orbit is now well constrained. With a semi-major axis (sma) of $a = 9.0 \pm 0.5$ au (1 $\sigma $), it definitely excludes previously reported possible long orbital periods, and excludes \bpic b as the origin of photometric variations that took place in 1981. We also refine the eccentricity and inclination of the planet. From an instrumental point of view, these data demonstrate that it is possible to detect, if they exist, young massive Jupiters that orbit at less than 2 au from a star that is 20 pc away.

Context. A low-mass brown dwarf has been recently imaged around HR 2562 (HD 50571), a star hosting a debris disk resolved in the far infrared. Interestingly, the companion location is compatible with an orbit coplanar with the disk and interior to the debris belt. This feature makes the system a valuable laboratory to analyze the formation of substellar companions in a circumstellar disk and potential disk-companion dynamical interactions. Aims. We aim to further characterize the orbital motion of HR 2562 B and its interactions with the host star debris disk. Methods. We performed a monitoring of the system over ~10 months in 2016 and 2017 with the VLT/SPHERE exoplanet imager. Results. We confirm that the companion is comoving with the star and detect for the first time an orbital motion at high significance, with a current orbital motion projected in the plane of the sky of 25 mas (~0.85 au) per year. No orbital curvature is seen in the measurements. An orbital fit of the SPHERE and literature astrometry of the companion without priors on the orbital plane clearly indicates that its orbit is (quasi-)coplanar with the disk. To further constrain the other orbital parameters, we used empirical laws for a companion chaotic zone validated by N-body simulations to test the orbital solutions that are compatible with the estimated disk cavity size. Non-zero eccentricities (>0.15) are allowed for orbital periods shorter than 100 yr, while only moderate eccentricities up to ~0.3 for orbital periods longer than 200 yr are compatible with the disk observations. A comparison of synthetic Herschel images to the real data does not allow us to constrain the upper eccentricity of the companion.

Young circumstellar disks are of prime interest to understand the physical and chemical conditions under which planet formation takes place. Only very few detections of planet candidates within these disks exist, and most of them are currently suspected to be disk features. In this context, the transition disk around the young star PDS 70 is of particular interest, due to its large gap identified in previous observations, indicative of ongoing planet formation. We aim to search for the presence of planets and search for disk structures indicative for disk-planet interactions and other evolutionary processes. We analyse new and archival near-infrared (NIR) images of the transition disk PDS 70 obtained with the VLT/SPHERE, VLT/NaCo and Gemini/NICI instruments in polarimetric differential imaging (PDI) and angular differential imaging (ADI) modes. We detect a point source within the gap of the disk at about 195 mas (about 22 au) projected separation. The detection is confirmed at five different epochs, in three filter bands and using different instruments. The astrometry results in an object of bound nature, with high significance. The comparison of the measured magnitudes and colours to evolutionary tracks suggests that the detection is a companion of planetary mass. We confirm the detection of a large gap of about 54 au in size within the disk in our scattered light images, and detect a signal from an inner disk component. We find that its spatial extent is very likely smaller than about 17 au in radius. The images of the outer disk show evidence of a complex azimuthal brightness distribution which may in part be explained by Rayleigh scattering from very small grains. Future observations of this system at different wavelengths and continuing astrometry will allow us to test theoretical predictions regarding planet-disk interactions, planetary atmospheres and evolutionary models.

The G-type star GJ504A is known to host a 3 to 35 MJup companion whose temperature, mass, and projected separation all contribute to make it a test case for the planet formation theories and for atmospheric models of giant planets and light brown dwarfs. We collected data from the CHARA interferometer, SOPHIE spectrograph, and VLT/SPHERE high contrast imager to revisit the properties of the system. We measure a radius of 1.35+/- 0.04Rsun for GJ504A which yields isochronal ages of 21+/-2Myr or 4.0+/-1.8Gyr for the system and line-of-sight stellar rotation axis inclination of $162.4_{-4.3}^{+3.8}$ degrees or $18.6_{-3.8}^{+4.3}$ degrees. We re-detect the companion in the Y2, Y3, J3, H2, and K1 dual band SPHERE images. The complete 1-4 $\mu$m SED shape of GJ504b is best reproduced by T8-T9.5 objects with intermediate ages ($\leq1.5$Gyr), and/or unusual dusty atmospheres and/or super-solar metallicities. All six atmospheric models used yield $\mathrm{T_{eff}=550 \pm 50}$K for GJ504b and point toward a low surface gravity (3.5-4.0 dex). The accuracy on the metallicity value is limited by model-to-model systematics. It is not degenerate with the C/O ratio. We derive $\mathrm{log\:L/L_{\odot}=-6.15\pm0.15}$ dex for the companion compatible with masses of $\mathrm{M=1.3^{+0.6}_{-0.3}M_{Jup}}$ and $\mathrm{M=23^{+10}_{-9} M_{Jup}}$ for the young and old age ranges, respectively. The semi-major axis (sma) is above 27.8 au and the eccentricity lower than 0.55. The posterior on GJ~504b's orbital inclination suggests a misalignment with GJ~504A rotation axis. We combine the radial velocity and multi-epoch imaging data to exclude additional objects (90\% prob.) more massive than 2.5 and 30 $\mathrm{M_{Jup}}$ with sma in the range 0.01-80 au for the young and old system ages, respectively. The companion is in the envelope of the population of planets synthetized with our core-accretion model.

Evolutionary models are widely used to infer the mass of stars, brown dwarfs, and giant planets. Their predictions are thought to be less reliable at young ages ($<$ 200 Myr) and in the low-mass regime ($\mathrm{<1~M_{\odot}}$). GJ 2060 AB and TWA 22 AB are two rare astrometric M-dwarf binaries respectively members of the AB Doradus and Beta Pictoris moving groups. As their dynamical mass can be measured within a few years, they can be used to calibrate the evolutionary tracks and set new constraints on the age of young moving groups. We find a total mass of $\mathrm{0.18\pm 0.02~M_\odot}$ for TWA 22. That mass is in good agreement with model predictions at the age of the Beta Pic moving group. We obtain a total mass of $\mathrm{1.09 \pm 0.10~M_{\odot}}$ for GJ 2060. We estimate a spectral type of M$1\pm0.5$, $\mathrm{L/L_{\odot}=-1.20\pm0.05}$ dex, and $\mathrm{T_{eff}=3700\pm100}$ K for GJ 2060 A. The B component is a M$3\pm0.5$ dwarf with $\mathrm{L/L_{\odot}=-1.63\pm0.05}$ dex and $\mathrm{T_{eff}=3400\pm100}$ K. The dynamical mass of GJ 2060 AB is inconsistent with the most recent models predictions (BCAH15, PARSEC) for an ABDor age in the range 50-150 Myr. It is 10 to 20\% (1-2 sigma, depending on the assumed age) above the models predictions, corresponding to an underestimation of $0.10$ to $0.20~\mathrm{M_\odot}$. Coevality suggests a young age for the system ($\sim$ 50 Myr) according to most evolutionary models. TWA 22 validates the predictions of recent evolutionary tracks at $\sim$20 Myr. On the other hand, we evidence a 1-2 sigma mismatch between the predicted and observed mass of GJ 2060 AB. This slight departure may indicate that one of the star hosts a tight companion. Alternatively, this would confirm the models tendency to underestimate the mass of young low-mass stars.

The consortium of the Spectro-Polarimetric High-contrast Exoplanet REsearch installed at the Very Large Telescope (SPHERE/VLT) has been operating its guaranteed observation time (260 nights over five years) since February 2015. The main part of this time (200 nights) is dedicated to the detection and characterization of young and giant exoplanets on wide orbits. The large amount of data must be uniformly processed so that accurate and homogeneous measurements of photometry and astrometry can be obtained for any source in the field. To complement the European Southern Observatory pipeline, the SPHERE consortium developed a dedicated piece of software to process the data. First, the software corrects for instrumental artifacts. Then, it uses the speckle calibration tool (SpeCal) to minimize the stellar light halo that prevents us from detecting faint sources like exoplanets or circumstellar disks. SpeCal is meant to extract the astrometry and photometry of detected point-like sources (exoplanets, brown dwarfs, or background sources). SpeCal was intensively tested to ensure the consistency of all reduced images (cADI, Loci, TLoci, PCA, and others) for any SPHERE observing strategy (ADI, SDI, ASDI as well as the accuracy of the astrometry and photometry of detected point-like sources. SpeCal is robust, user friendly, and efficient at detecting and characterizing point-like sources in high contrast images. It is used to process all SPHERE data systematically, and its outputs have been used for most of the SPHERE consortium papers to date. SpeCal is also a useful framework to compare different algorithms using various sets of data (different observing modes and conditions). Finally, our tests show that the extracted astrometry and photometry are accurate and not biased.

We report the discovery of a bright, brown dwarf companion to the star HIP 64892, imaged with VLT/SPHERE during the SHINE exoplanet survey. The host is a B9.5V member of the Lower-Centaurus-Crux subgroup of the Scorpius Centaurus OB association. The measured angular separation of the companion ($1.2705\pm0.0023$") corresponds to a projected distance of $159\pm12$ AU. We observed the target with the dual-band imaging and long-slit spectroscopy modes of the IRDIS imager to obtain its SED and astrometry. In addition, we reprocessed archival NACO L-band data, from which we also recover the companion. Its SED is consistent with a young (<30 Myr), low surface gravity object with a spectral type of M9$_{\gamma}\pm1$. From comparison with the BT-Settl atmospheric models we estimate an effective temperature of $T_{\textrm{eff}}=2600 \pm 100$ K, and comparison of the companion photometry to the COND evolutionary models yields a mass of $\sim29-37$ M$_{\text{J}}$ at the estimated age of $16^{+15}_{-7}$ Myr for the system. HIP 64892 is a rare example of an extreme-mass ratio system ($q\sim0.01$) and will be useful for testing models relating to the formation and evolution of such low-mass objects.

The SHINE program is a large high-contrast near-infrared survey of 600 young, nearby stars. It is aimed at searching for and characterizing new planetary systems using VLT/SPHERE's unprecedented high-contrast and high-angular resolution imaging capabilities. It also intends at placing statistical constraints on the occurrence and orbital properties of the giant planet population at large orbits as a function of the stellar host mass and age to test planet formation theories. We use the IRDIS dual-band imager and the IFS integral field spectrograph of SPHERE to acquire high-constrast coronagraphic differential near-infrared images and spectra of the young A2 star HIP65426. It is a member of the ~17 Myr old Lower Centaurus-Crux association. At a separation of 830 mas (92 au projected) from the star, we detect a faint red companion. Multi-epoch observations confirm that it shares common proper motion with HIP65426. Spectro-photometric measurements extracted with IFS and IRDIS between 0.95 and 2.2um indicate a warm, dusty atmosphere characteristic of young low surface-gravity L5-L7 dwarfs. Hot-start evolutionary models predict a luminosity consistent with a 6-12 MJup, Teff=1300-1600 K and R=1.5 RJup giant planet. Finally, the comparison with Exo-REM and PHOENIX BT-Settl synthetic atmosphere models gives consistent effective temperatures but with slightly higher surface gravity solutions of log(g)=4.0-5.0 with smaller radii (1.0-1.3 RJup). Given its physical and spectral properties, HIP65426b occupies a rather unique placement in terms of age, mass and spectral-type among the currently known imaged planets. It represents a particularly interesting case to study the presence of clouds as a function of particle size, composition, and location in the atmosphere, to search for signatures of non-equilibrium chemistry, and finally to test the theory of planet formation and evolution.

HD$\,$169142 is an excellent target to investigate signs of planet-disk interaction due to the previous evidence of gap structures. We performed J-band (~1.2{\mu}m) polarized intensity imaging of HD169142 with VLT/SPHERE. We observe polarized scattered light down to 0.16" (~19 au) and find an inner gap with a significantly reduced scattered light flux. We confirm the previously detected double ring structure peaking at 0.18" (~21 au) and 0.56" (~66 au), and marginally detect a faint third gap at 0.70"-0.73" (~82-85 au). We explore dust evolution models in a disk perturbed by two giant planets, as well as models with a parameterized dust size distribution. The dust evolution model is able to reproduce the ring locations and gap widths in polarized intensity, but fails to reproduce their depths. It, however, gives a good match with the ALMA dust continuum image at 1.3 mm. Models with a parameterized dust size distribution better reproduce the gap depth in scattered light, suggesting that dust filtration at the outer edges of the gaps is less effective. The pile-up of millimeter grains in a dust trap and the continuous distribution of small grains throughout the gap likely require a more efficient dust fragmentation and dust diffusion in the dust trap. Alternatively, turbulence or charging effects might lead to a reservoir of small grains at the surface layer that is not affected by the dust growth and fragmentation cycle dominating the dense disk midplane. The exploration of models shows that extracting planet properties such as mass from observed gap profiles is highly degenerate.

Transition disks correspond to a short stage between the young protoplanetary phase and older debris phase. Along this evolutionary sequence, the gas component disappears leaving room for a dust-dominated environment where already-formed planets signpost their gravitational perturbations. We endeavor to study the very inner region of the well-known and complex debris, but still gas-rich disk, around HD 141569A using the exquisite high-contrast capability of SPHERE at the VLT. Recent near-infrared (IR) images suggest a relatively depleted cavity within ~200 au, while former mid-IR data indicate the presence of dust at separations shorter than ~100 au. We obtained multi-wavelength images in the near-IR in J, H2, H3 and Ks bands with the IRDIS camera and a 0.95-1.35 micrometers spectral data cube with the IFS. Data were acquired in pupil-tracking mode, thus allowing for angular differential imaging. We discovered several new structures inside 1", of which the most prominent is a bright ring with sharp edges (semi-major axis: 0.4") featuring a strong north-south brightness asymmetry. Other faint structures are also detected from 0.4" to 1" in the form of concentric ringlets and at least one spiral arm. Finally, the VISIR data at 8.6 micrometers suggests the presence of an additional dust population closer in. Besides, we do not detect companions more massive than 1-3 mass of Jupiter. The performance of SPHERE allows us to resolve the extended dust component, which was previously detected at thermal and visible wavelengths, into very complex patterns with strong asymmetries ; the nature of these asymmetries remains to be understood. Scenarios involving shepherding by planets or dust-gas interactions will have to be tested against these observations.

A measurement of the cosmic-ray spectrum for energies exceeding $4{\times}10^{18}$ eV is presented, which is based on the analysis of showers with zenith angles greater than $60^{\circ}$ detected with the Pierre Auger Observatory between 1 January 2004 and 31 December 2013. The measured spectrum confirms a flux suppression at the highest energies. Above $5.3{\times}10^{18}$ eV, the "ankle", the flux can be described by a power law $E^{-\gamma}$ with index $\gamma=2.70 \pm 0.02 \,\text{(stat)} \pm 0.1\,\text{(sys)}$ followed by a smooth suppression region. For the energy ($E_\text{s}$) at which the spectral flux has fallen to one-half of its extrapolated value in the absence of suppression, we find $E_\text{s}=(5.12\pm0.25\,\text{(stat)}^{+1.0}_{-1.2}\,\text{(sys)}){\times}10^{19}$ eV.

Neutrinos in the cosmic ray flux with energies near 1 EeV and above are detectable with the Surface Detector array of the Pierre Auger Observatory. We report here on searches through Auger data from 1 January 2004 until 20 June 2013. No neutrino candidates were found, yielding a limit to the diffuse flux of ultra-high energy neutrinos that challenges the Waxman-Bahcall bound predictions. Neutrino identification is attempted using the broad time-structure of the signals expected in the SD stations, and is efficiently done for neutrinos of all flavors interacting in the atmosphere at large zenith angles, as well as for "Earth-skimming" neutrino interactions in the case of tau neutrinos. In this paper the searches for downward-going neutrinos in the zenith angle bins $60^\circ-75^\circ$ and $75^\circ-90^\circ$ as well as for upward-going neutrinos, are combined to give a single limit. The $90\%$ C.L. single-flavor limit to the diffuse flux of ultra-high energy neutrinos with an $E^{-2}$ spectrum in the energy range $1.0 \times 10^{17}$ eV - $2.5 \times 10^{19}$ eV is $E_\nu^2 dN_\nu/dE_\nu < 6.4 \times 10^{-9}~ {\rm GeV~ cm^{-2}~ s^{-1}~ sr^{-1}}$.

The planetary system discovered around the young A-type HR8799 provides a unique laboratory to: a) test planet formation theories, b) probe the diversity of system architectures at these separations, and c) perform comparative (exo)planetology. We present and exploit new near-infrared images and integral-field spectra of the four gas giants surrounding HR8799 obtained with SPHERE, the new planet finder instrument at the Very Large Telescope, during the commissioning and science verification phase of the instrument (July-December 2014). With these new data, we contribute to completing the spectral energy distribution of these bodies in the 1.0-2.5 $\mu$m range. We also provide new astrometric data, in particular for planet e, to further constrain the orbits. We used the infrared dual-band imager and spectrograph (IRDIS) subsystem to obtain pupil-stabilized, dual-band $H2H3$ (1.593 $\mu$m, 1.667 $\mu$m), $K1K2$ (2.110 $\mu$m, 2.251 $\mu$m), and broadband $J$ (1.245 $\mu$m) images of the four planets. IRDIS was operated in parallel with the integral field spectrograph (IFS) of SPHERE to collect low-resolution ($R\sim30$), near-infrared (0.94-1.64 $\mu$m) spectra of the two innermost planets HR8799d and e. The data were reduced with dedicated algorithms, such as the Karhunen-Lo\`eve image projection (KLIP), to reveal the planets. We used the so-called negative planets injection technique to extract their photometry, spectra, and measure their positions. We illustrate the astrometric performance of SPHERE through sample orbital fits compatible with SPHERE and literature data.

The system of four planets around HR8799 offers a unique opportunity to probe the physics and chemistry at play in the atmospheres of self-luminous young (~30 Myr) planets. We recently obtained new photometry of the four planets and low-resolution (R~30) spectra of HR8799 d and e with the SPHERE instrument (paper III). In this paper (paper IV), we compare the available spectra and photometry of the planets to known objects and atmospheric models (BT-SETTL14, Cloud-AE60, Exo-REM) to characterize the atmospheric properties of the planets. We find that HR8799d and e properties are well reproduced by those of L6-L8 dusty dwarfs discovered in the field, among which some are candidate members of young nearby associations. No known object reproduces well the properties of planets b and c. Nevertheless, we find that the spectra and WISE photometry of peculiar and/or young early-T dwarfs reddened by submicron grains made of corundum, iron, enstatite, or forsterite successfully reproduce the SED of these two planets. Our analysis confirms that only the Exo-REM models with thick clouds fit (within 2{\sigma}) the whole set of spectrophotometric datapoints available for HR8799 d and e for Teff = 1200 K, log g in the range 3.0-4.5, and M/H=+0.5. The models still fail to reproduce the SED of HR8799c and b. The determination of the metallicity, log g, and cloud thickness are degenerate. We conclude that an enhanced content in dust and decreased CIA of H2 is certainly responsible for the deviation of the properties of the planet with respect to field dwarfs. The analysis suggests in addition that HR8799c and b have later spectral types than the two other planets, and therefore could both have lower masses.

The ExTrA facility, located at La Silla observatory, will consist of a near-infrared multi-object spectrograph fed by three 60-cm telescopes. ExTrA will add the spectroscopic resolution to the traditional differential photometry method. This shall enable the fine correction of color-dependent systematics that would otherwise hinder ground-based observations. With both this novel method and an infrared-enabled efficiency, ExTrA aims to find transiting telluric planets orbiting in the habitable zone of bright nearby M dwarfs. It shall have the versatility to do so by running its own independent survey and also by concurrently following-up on the space candidates unveiled by K2 and TESS. The exoplanets detected by ExTrA will be amenable to atmospheric characterisation with VLTs, JWST, and ELTs and could give our first peek into an exo-life laboratory.

Spherical harmonic moments are well-suited for capturing anisotropy at any scale in the flux of cosmic rays. An unambiguous measurement of the full set of spherical harmonic coefficients requires full-sky coverage. This can be achieved by combining data from observatories located in both the northern and southern hemispheres. To this end, a joint analysis using data recorded at the Telescope Array and the Pierre Auger Observatory above $10^{19}$ eV is presented in this work. The resulting multipolar expansion of the flux of cosmic rays allows us to perform a series of anisotropy searches, and in particular to report on the angular power spectrum of cosmic rays above $10^{19}$ eV. No significant deviation from isotropic expectations is found throughout the analyses performed. Upper limits on the amplitudes of the dipole and quadrupole moments are derived as a function of the direction in the sky, varying between 7% and 13% for the dipole and between 7% and 10% for a symmetric quadrupole.

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.

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

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 visitor instrument PIONIER provides VLTI with improved imaging capabilities and sensitivity. The instrument started routinely delivering scientific data in November 2010, that is less than 12 months after being approved by the ESO Science and Technical Committee. We recall the challenges that had to be tackled to design, built and commission PIONIER. We summarize the typical performances and some astrophysical results obtained so far. We conclude this paper by summarizing lessons learned.

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.