SPHERE (Beuzit et al,. 2019) has now been in operation at the VLT for more than 5 years, demonstrating a high level of performance. SPHERE has produced outstanding results using a variety of operating modes, primarily in the field of direct imaging of exoplanetary systems, focusing on exoplanets as point sources and circumstellar disks as extended objects. The achievements obtained thus far with SPHERE (~200 refereed publications) in different areas (exoplanets, disks, solar system, stellar physics...) have motivated a large consortium to propose an even more ambitious set of science cases, and its corresponding technical implementation in the form of an upgrade. The SPHERE+ project capitalizes on the expertise and lessons learned from SPHERE to push high contrast imaging performance to its limits on the VLT 8m-telescope. The scientific program of SPHERE+ described in this document will open a new and compelling scientific window for the upcoming decade in strong synergy with ground-based facilities (VLT/I, ELT, ALMA, and SKA) and space missions (Gaia, JWST, PLATO and WFIRST). While SPHERE has sampled the outer parts of planetary systems beyond a few tens of AU, SPHERE+ will dig into the inner regions around stars to reveal and characterize by mean of spectroscopy the giant planet population down to the snow line. Building on SPHERE's scientific heritage and resounding success, SPHERE+ will be a dedicated survey instrument which will strengthen the leadership of ESO and the European community in the very competitive field of direct imaging of exoplanetary systems. With enhanced capabilities, it will enable an even broader diversity of science cases including the study of the solar system, the birth and death of stars and the exploration of the inner regions of active galactic nuclei.

We present near-infrared interferometric data on the Seyfert 2 galaxy NGC 1068, obtained with the GRAVITY instrument on the European Southern Observatory Very Large Telescope Interferometer. The extensive baseline coverage from 5 to 60 M\lambda allowed us to reconstruct a continuum image of the nucleus with an unrivaled 0.2 pc resolution in the K-band. We find a thin ring-like structure of emission with a radius r = 0.24+/-0.03 pc, inclination i = 70+/-5 deg, position angle PA = -50+/-4 deg, and h/r < 0.14, which we associate with the dust sublimation region. The observed morphology is inconsistent with the expected signatures of a geometrically and optically thick torus. Instead, the infrared emission shows a striking resemblance to the 22 GHz maser disc, which suggests they share a common region of origin. The near-infrared spectral energy distribution indicates a bolometric luminosity of (0.4-4.7) x 10^45 erg/s, behind a large A_K ~ 5.5 (A_V ~ 90) screen of extinction that also appears to contribute significantly to obscuring the broad line region.

We investigate the fueling and the feedback of nuclear activity in the Seyfert 2 galaxy NGC1068, by studying the distribution and kinematics of molecular gas in the torus and its connections. We use ALMA to image the emission of a set of molecular gas tracers in the circumnuclear disk (CND) and the torus of the galaxy using the CO(2-1), CO(3-2) and HCO+(4-3) lines with spatial resolutions ~0.03"-0.09"(2-6pc). ALMA resolves the CND as an asymmetric ringed disk of D~400pc-size and mass of ~1.4x10^8 Msun. The inner edge of the ring is associated with edge-brightened arcs of NIR polarized emission identified with the working surface of the AGN ionized wind. ALMA proves the existence of a molecular torus of M_torus ~ 3x10^5Msun, which extends over a large range of spatial scales D=10-30pc around the central engine. The new observations evidence the density radial stratification of the torus: the HCO+(4-3) torus, with a full size D=11pc, is a factor of 2-3 smaller than its CO(2-1) and CO(3-2) counterparts, which have full-sizes D=26pc and D=28pc, respectively. The torus is connected to the CND through a network of gas streamers. The kinematics of molecular gas show strong departures from circular motions in the torus, the gas streamers, and the CND. These velocity distortions are interconnected and are part of a 3D outflow that reflects the effects of AGN feedback on the kinematics of molecular gas across a wide range of spatial scales. We conclude that a wide-angle AGN wind launched from the accretion disk is impacting a sizeable fraction of the gas inside the torus (~0.4-0.6 x M_torus). However, a large gas reservoir (~1.2-1.8 x 10^5Msun) close to the equatorial plane of the torus remains unaffected by the AGN wind and can continue fueling the AGN for ~1-4Myr. AGN fueling seems nevertheless thwarted on intermediate scales (15pc < r < 50pc).

We use VLTI/GRAVITY near-infrared interferometry measurements of 8 bright, Type 1 AGN to study the size and structure of hot dust heated by the central engine. We partially resolve each source, and report Gaussian FWHM sizes in the range 0.3-0.8 milliarcseconds. In all but one object, we find no evidence for significant elongation or asymmetry (closure phases < 1 deg). The effective physical radius increases with bolometric luminosity as found from past reverberation and interferometry measurements. The measured sizes for Seyfert galaxies are systematically larger than for the two quasars in our sample when measured relative to the previously reported R ~ L^1/2 relationship explained by emission at the sublimation radius. This could be evidence of evolving near-infrared emission region structure as a function of central luminosity.

Aims: We characterise the properties of stars, dust, and gas and their spatial distribution in the central region of the Seyfert 2 galaxy NGC1068. Method: Our study is based on near-infrared (YJH, 0.95-1.650 um, R=350) long-slit spectroscopy observations of the central region of NGC 1068 with a 0.4" spatial resolution. We decomposed the observed continuum emission into three components: hot dust, stars, and scattered light from the central engine. We measured their contributions at various distances from the nucleus. We also measured fluxes and Doppler shifts for the emission lines in our spectrum to probe the physical conditions of the narrow line region. Results: Dust and stars are the main sources of continuum emission, but scattered light from the central engine has also been detected in the very central region. Together, these three components reproduce the observed continuum well. The dust emission is compatible with a 830~K blackbody. It has only been detected in the very central region and is not spatially resolved. The stellar content is ubiquitous. It harbours a 250 pc cusp centred around the nucleus, over-imposed on a young stellar background. The spectrum of the cusp is consistent with a 120 Myr old single stellar population. Finally, the emission lines exhibit a significant Doppler shift that is consistent with a radial outflow from the nucleus in a biconical structure. The $\left[Fe\ II\right]$ behaviour strongly differs from other lines, indicating that it arises from a different structure.

Establishing precisely how stars and interstellar medium distribute within the central 100 pc area around an AGN, down to the pc scale, is key for understanding how the very final transfer of matter from kpc scale to the sub-parsec size of the accretion disc is achieved. Using AO-assisted (SPHERE-VLT) near-IR images in H and Ks and narrow-band of the Seyfert 2 galaxy NGC1068 we analyse the radial distribution of brightness in the central r < 100 pc area down to the pc scale. The median-averaged radial profiles are adjusted by a cusp (power-law) plus a central point-source. A simple radiative transfer model is used to interpret the data. We find that the fit of profiles beyond 10pc is done quite precisely at Ks by a cusp of exponent -2.0 plus a central point-source and by a cusp of exponent -1.2 at H. The difference between H and Ks can be explained by differential extinction, provided that the distribution of dust is itself cuspy as r^-1. However, the required stellar density follows a r^-4 cusp, much steeper than any other cusp theoretically predicted and the mass of the cluster is unreasonable, even introducing a segregation in the stellar population with an excess of giant stars inward. A much more acceptable solution is found with a K profile dominated by warm dust emission, while the H profile corresponds to a stellar cusp. NGC1068 is shown to satisfy a relationship between half-light radius, cusp luminosity and exponent which suggests that the cusp is the remnant of a recent starbust. We identify the central point-like source with the very hot dust at the internal wall of the putative torus and derive an intrinsic luminosity that requires an overall extinction AK ~ 8, a value consistent with predictions by several models.

The broadening of atomic emission lines by high-velocity motion of gas near accreting supermassive black holes is an observational hallmark of quasars. Observations of broad emission lines could potentially constrain the mechanism for transporting gas inwards through accretion disks or outwards through winds. The size of this broad-line region has been estimated by measuring the light travel time delay between the variable nuclear continuum and the emission lines - a method known as reverberation mapping. In some models the emission lines arise from a continuous outflow, whereas in others they are produced by orbiting gas clouds. Directly imaging such regions has not hitherto been possible because of their small angular sizes (< 0.1 milli-arcseconds). Here we report a spatial offset (with a spatial resolution of ten micro-arcseconds or about 0.03 parsecs for a distance of 550 million parsecs) between the red and blue photo-centres of the broad Paschen-{\alpha} line of the quasar 3C 273 perpendicular to the direction of its radio jet. This spatial offset corresponds to a gradient in the velocity of the gas and thus implies that the gas is orbiting the central supermassive black hole. The data are well fitted by a broad-line-region model of a thick disk of gravitationally bound material orbiting a black hole of 300 million solar masses. We infer a disk radius of 150 light days; a radius of 100-400 light days was found previously using reverberation mapping. The rotation axis of the disk aligns in inclination and position angle with the radio jet. Our results support the methods that are often used to estimate the masses of accreting supermassive black holes and to study their evolution over cosmic time.

Adaptive optics (AO) system performance is improved using post-processing techniques, such as point spread function (PSF) deconvolution. The PSF estimation involves characterization of the different wavefront (WF) error sources in the AO system. We propose a numerical error breakdown estimation tool that allows studying AO error source behavior such as their correlations. We also propose a new analytical model for anisoplanatism and bandwidth errors that were validated with the error breakdown estimation tool. This model is the first step for a complete AO residual error model that is expressed in deformable mirror space, leading to practical usage such as PSF reconstruction or turbulent parameters identification. We have developed in the computing platform for adaptive optics systems (COMPASS) code, which is an end-to-end simulation code using graphics processing units (GPU) acceleration, an estimation tool that provides a comprehensive error breakdown by the outputs of a single simulation run. We derive the various contributors from the end-to-end simulator at each iteration step: this method provides temporal buffers of each contributor. Then, we use this tool to validate a new model of anisoplanatism and bandwidth errors including their correlation. This model is based on a statistical approach that computes the error covariance matrices using structure functions. A correlation analysis shows significant correlations between some contributors, especially WF measurement deviation error and bandwidth error due to centroid gain, and the well-known correlation between bandwidth and anisoplanatism errors is also retrieved. The model we propose for the two latter errors shows an SR and EE difference of about one percent compared to the end-to-end simulation, even if some approximations exist.

Aims: In this paper, we aim to constrain the properties of dust structures in the central first parsecs of Active Galactic Nuclei (AGN). Our goal is to study the required optical depth and composition of different dusty and ionised structures. Methods: We developed a radiative transfer code, MontAGN, optimised for polarimetric observations in the infrared. With both this code and STOKES, designed to be relevant from the hard X-ray band to near-infrared wavelengths, we investigate the polarisation emerging from a characteristic model of the AGN environment. For that purpose, we compare predictions of our models with previous infrared observations of NGC 1068, and try to reproduce several key polarisation patterns revealed by polarisation mapping. Results: We constrain the required dust structures as well as their densities. More precisely, we find out that the electron density inside the ionisation cone is about $2.0 \times 10^9$ m$^{-3}$ . With structures constituted of spherical grains of constant density, we also highlight that the torus should be thicker than 20 in term of K band optical depth to block direct light from the centre. It should also have a stratification in density, with a lesser dense outer rim with an optical depth at 2.2 {\mu}m typically between 0.8 and 4 for observing the double scattering effect previously proposed. Conclusions: We bring constraints on the dust structures in the inner parsecs of an AGN model supposed to describe NGC 1068. When compared to observations, this leads to optical depth of at least 20 in Ks band for the torus of NGC 1068, corresponding to $\tau_V \approx 170$, which is within the range of current estimation based on observations. In the future, we will improve our study by including non uniform dust structures and aligned elongated grains to constrain other possible interpretations of the observations.

The performance of adaptive optics systems is partially dependant on the algorithms used within the real-time control system to compute wavefront slope measurements. We demonstrate use of a matched filter algorithm for the processing of elongated laser guide star (LGS) Shack-Hartmann images, using the CANARY adaptive optics instrument on the 4.2m William Herschel Telescope and the European Southern Observatory Wendelstein LGS Unit placed 40m away. This algorithm has been selected for use with the forthcoming Thirty Meter Telescope, but until now had not been demonstrated on-sky. From the results of a first observing run, we show that the use of matched filtering improves our adaptive optics system performance, with increases in on-sky H-band Strehl measured up to about a factor of 1.1 with respect to a conventional centre of gravity approach. We describe the algorithm used, and the methods that we implemented to enable on-sky demonstration.

CANARY is the multi-object adaptive optics (MOAO) on-sky pathfinder developed in the perspective of Multi-Object Spectrograph on Extremely Large Telescopes~(ELTs). In 2013, CANARY was operated on-sky at the William Herschel telescope~(WHT), using three off-axis natural guide stars~(NGS) and four off-axis Rayleigh laser guide stars~(LGS), in open-loop, with the on-axis compensated turbulence observed with a H-band imaging camera and a Truth wave-front sensor~(TS) for diagnostic purposes. Our purpose is to establish a reliable and accurate wave-front error breakdown for LGS MOAO. This will enable a comprehensive analysis of \cana on-sky results and provide tools for validating simulations of MOAO systems for ELTs. To evaluate the MOAO performance, we compared the CANARY on-sky results running in MOAO, in Single Conjugated Adaptive Optics~(SCAO) and in Ground Layer Adaptive Optics~(GLAO) modes, over a large set of data acquired in 2013. We provide a statistical study of the seeing. We also evaluated the wave-front error breakdown from both analytic computations, one based on a MOAO system modelling and the other on the measurements from the CANARY TS. We have focussed especially on the tomographic error and we detail its vertical error decomposition~(VED). We show that CANARY obtained 30.1\%, 21.4\% and 17.1\% H-band Strehl ratios in SCAO, MOAO and GLAO respectively, for median seeing conditions with 0.66" of total seeing including 0.59" at the ground. Moreover, we get 99\% of correlation over 4,500 samples, for any AO modes, between two analytic computations of residual phase variance. Based on these variances, we obtain a reasonable Strehl-ratio~(SR) estimation when compared to the measured IR image SR. We evaluate the gain in compensation for the altitude turbulence brought by MOAO when compared to GLAO.

In preparation of future Multi-Object Spectrographs (MOS) whose one of the major role is to provide an extensive statistical studies of high redshifted galaxies surveyed, the demonstrator Canary has been designed to tackle technical challenges related to open-loop Adaptive-Optics (AO) control with jointed Natural Guide Star (NGS) and Laser Guide Star (LGS) tomography. We have developed a Point Spread Function (PSF)-Reconstruction algorithm dedicated to MOAO systems using system telemetry to estimate the PSF potentially anywhere in the observed field, a prerequisite to post- process AO-corrected observations in Integral Field Spectroscopy (IFS). In this paper we show how to handle off-axis data to estimate the PSF using atmospheric tomography and compare it to a classical approach that uses on-axis residual phase from a truth sensor observing a natural bright source. We have reconstructed over 450 on-sky Canary PSFs and we get bias/1-${\sigma}$ standard-deviation (std) of 1.3/4.8 on the H-band Strehl ratio (SR) with 92.3% of correlation between reconstructed and sky SR. On the Full Width at Half Maximum (FWHM), we get respectively 2.94 mas, 19.9 mas and 88.3% for the bias, std and correlation. The reference method achieves 0.4/3.5/95% on the SR and 2.71 mas/14.9 mas/92.5% on the FWHM for the bias/std/correlation.

In this first research note of a series of two, we present a comparison between two Monte Carlo radiative transfer codes: MontAGN and STOKES. Both were developed in order to better understand the observed polarisation of Active Galactic Nuclei (AGN). Our final aim is to use these radiative transfer codes to simulate the polarisation maps of a prototypical type-2 radio-quiet AGN on a wide range of wavelengths, from the infrared band with MontAGN to the X-ray energies with STOKES. Doing so, we aim to analyse in depth the recent SPHERE/IRDIS polarimetric observations conducted on NGC 1068. In order to validate the codes and obtain preliminary results, we set for both codes a common and simple AGN model, and compared their polaro-imaging results.

In this second research note of a series of two, we present the first near-infrared results we obtained when modeling Active Galactic Nuclei (AGN). Our first proceedings showed the comparison between the MontAGN and STOKES Monte Carlo codes. Now we use our radiative transfer codes to simulate the polarization maps of a prototypical, NGC 1068-like, type-2 radio-quiet AGN. We produced high angular resolution infrared (1 micron) polarization images to be compared with recent observations in this wavelength range. Our preliminary results already show a good agreement between the models and observations but cannot account for the peculiar linear polarization angle of the torus such as observed. Gratadour et al. 2015 found a polarization position angle being perpendicular to the bipolar outflows axis. Further work is needed to improve the models by adding physical phenomena such as dichroism and clumpiness.

With the aim of paving the road for future accurate astrometry with MICADO at the European-ELT, we performed an astrometric study using two different but complementary approaches to investigate two critical components that contribute to the total astrometric accuracy. First, we tested the predicted improvement in the astrometric measurements with the use of an atmospheric dispersion corrector (ADC) by simulating realistic images of a crowded Galactic globular cluster. We found that the positional measurement accuracy should be improved by up to ~2 mas with the ADC, making this component fundamental for high-precision astrometry. Second, we analysed observations of a globular cluster taken with the only currently available Multi-Conjugate Adaptive Optics assisted camera, GeMS/GSAOI at Gemini South. Making use of previously measured proper motions of stars in the field of view, we were able to model the distortions affecting the stellar positions. We found that they can be as large as ~200 mas, and that our best model corrects them to an accuracy of ~1 mas. We conclude that future astrometric studies with MICADO requires both an ADC and an accurate modelling of distortions to the field of view, either through an a-priori calibration or an a-posteriori correction.

MICADO will equip the E-ELT with a first light capability for diffraction limited imaging at near-infrared wavelengths. The instrument's observing modes focus on various flavours of imaging, including astrometric, high contrast, and time resolved. There is also a single object spectroscopic mode optimised for wavelength coverage at moderately high resolution. This contribution provides an overview of the key functionality of the instrument, outlining the scientific rationale for its observing modes. The interface between MICADO and the adaptive optics system MAORY that feeds it is summarised. The design of the instrument is discussed, focussing on the optics and mechanisms inside the cryostat, together with a brief overview of the other key sub-systems.

Recent advances in adaptive optics (AO) have led to the implementation of wide field-of-view AO systems. A number of wide-field AO systems are also planned for the forthcoming Extremely Large Telescopes. Such systems have multiple wavefront sensors of different types, and usually multiple deformable mirrors (DMs). Here, we report on our experience integrating cameras and DMs with the real-time control systems of two wide-field AO systems. These are CANARY, which has been operating on-sky since 2010, and DRAGON, which is a laboratory adaptive optics real-time demonstrator instrument. We detail the issues and difficulties that arose, along with the solutions we developed. We also provide recommendations for consideration when developing future wide-field AO systems.

This White Paper presents the scientific motivations for a multi-object spectrograph (MOS) on the European Extremely Large Telescope (E-ELT). The MOS case draws on all fields of contemporary astronomy, from extra-solar planets, to the study of the halo of the Milky Way and its satellites, and from resolved stellar populations in nearby galaxies out to observations of the earliest 'first-light' structures in the partially-reionised Universe. The material presented here results from thorough discussions within the community over the past four years, building on the past competitive studies to agree a common strategy toward realising a MOS capability on the E-ELT. The cases have been distilled to a set of common requirements which will be used to define the MOSAIC instrument, entailing two observational modes ('high multiplex' and 'high definition'). When combined with the unprecedented sensitivity of the E-ELT, MOSAIC will be the world's leading MOS facility. In analysing the requirements we also identify a high-multiplex MOS for the longer-term plans for the E-ELT, with an even greater multiplex (>1000 targets) to enable studies of large-scale structures in the high-redshift Universe. Following the green light for the construction of the E-ELT the MOS community, structured through the MOSAIC consortium, is eager to realise a MOS on the E-ELT as soon as possible. We argue that several of the most compelling cases for ELT science, in highly competitive areas of modern astronomy, demand such a capability. For example, MOS observations in the early stages of E-ELT operations will be essential for follow-up of sources identified by the James Webb Space Telescope (JWST). In particular, multi-object adaptive optics and accurate sky subtraction with fibres have both recently been demonstrated on sky, making fast-track development of MOSAIC feasible.

MICADO is a near-IR camera for the Europea ELT, featuring an extended field (75" diameter) for imaging, and also spectrographic and high contrast imaging capabilities. It has been chosen by ESO as one of the two first-light instruments. Although it is ultimately aimed at being fed by the MCAO module called MAORY, MICADO will come with an internal SCAO system that will be complementary to it and will deliver a high performance on axis correction, suitable for coronagraphic and pupil masking applications. The basis of the pupil masking approach is to ensure the stability of the optical transfer function, even in the case of residual errors after AO correction (due to non common path errors and quasi-static aberrations). Preliminary designs of pupil masks are presented. Trade-offs and technical choices, especially regarding redundancy and pupil tracking, are explained.

Over the past 18 months we have revisited the science requirements for a multi-object spectrograph (MOS) for the European Extremely Large Telescope (E-ELT). These efforts span the full range of E-ELT science and include input from a broad cross-section of astronomers across the ESO partner countries. In this contribution we summarise the key cases relating to studies of high-redshift galaxies, galaxy evolution, and stellar populations, with a more expansive presentation of a new case relating to detection of exoplanets in stellar clusters. A general requirement is the need for two observational modes to best exploit the large (>40 sq. arcmin) patrol field of the E-ELT. The first mode ('high multiplex') requires integrated-light (or coarsely resolved) optical/near-IR spectroscopy of >100 objects simultaneously. The second ('high definition'), enabled by wide-field adaptive optics, requires spatially-resolved, near-IR of >10 objects/sub-fields. Within the context of the conceptual study for an ELT-MOS called MOSAIC, we summarise the top-level requirements from each case and introduce the next steps in the design process.

The first on-sky results obtained by CANARY, the Multi-Object Adaptive Optics (MOAO) demonstrator, are analysed. The data were recorded at the William Herschel Telescope, at the end of September 2010. We describe the command and calibrations algorithms used during the run and present the observing conditions. The processed data are MOAO-loop engaged or disengaged slope buffers, comprising the synchronised measurements of the four Natural Guide Stars (NGS) wavefront sensors running in parallel, and near Infra-Red (IR) images. We describe the method we use to establish the error budget of CANARY. We are able to evaluate the to- mographic and the open loop errors, having median values around 216 nm and 110 nm respectively. In addition, we identify an unexpected residual quasi-static field aberration term of mean value 110 nm. We present the detailed error budget analysed for three sets of data for three different asterisms. We compare the experimental budgets with the numerically simulated ones and demonstrate a good agreement. We find also a good agreement between the computed error budget from the slope buffers and the measured Strehl ratio on the IR images, ranging between 10% and 20% at 1 530 nm. These results make us confident in our ability to establish the error budget of future MOAO instruments.

We present recent results from the initial testing of an Artificial Neural Network (ANN) based tomographic reconstructor Complex Atmospheric Reconstructor based on Machine lEarNing (CARMEN) on Canary, an Adaptive Optics demonstrator operated on the 4.2m William Herschel Telescope, La Palma. The reconstructor was compared with contemporaneous data using the Learn and Apply (L&A) tomographic reconstructor. We find that the fully optimised L&A tomographic reconstructor outperforms CARMEN by approximately 5% in Strehl ratio or 15nm rms in wavefront error. We also present results for Canary in Ground Layer Adaptive Optics mode to show that the reconstructors are tomographic. The results are comparable and this small deficit is attributed to limitations in the training data used to build the ANN. Laboratory bench tests show that the ANN can out perform L&A under certain conditions, e.g. if the higher layer of a model two layer atmosphere was to change in altitude by ~300~m (equivalent to a shift of approximately one tenth of a subaperture).

The Gemini Multi-conjugate adaptive optics System (GeMS) at the Gemini South telescope in Cerro Pach{\'o}n is the first sodium-based multi-Laser Guide Star (LGS) adaptive optics system. It uses five LGSs and two deformable mirrors to measure and compensate for atmospheric distortions. The GeMS project started in 1999, and saw first light in 2011. It is now in regular operation, producing images close to the diffraction limit in the near infrared, with uniform quality over a field of view of two square arcminutes. The present paper (I) is the first one in a two-paper review of GeMS. It describes the system, explains why and how it was built, discusses the design choices and trade-offs, and presents the main issues encountered during the course of the project. Finally, we briefly present the results of the system first light.

We present spatially resolved photometric and spectroscopic observations of two wide brown dwarf binaries uncovered by the SIMP near-infrared proper motion survey. The first pair (SIMP J1619275+031350AB) has a separation of 0.691" (15.2 AU) and components T2.5+T4.0, at the cooler end of the ill-understood J-band brightening. The system is unusual in that the earlier-type primary is bluer in J-Ks than the later-type secondary, whereas the reverse is expected for binaries in the late-L to T dwarf range. This remarkable color reversal can possibly be explained by very different cloud properties between the two components. The second pair (SIMP J1501530-013506AB) consists of an L4.5+L5.5 (separation 0.96", 30-47 AU) with a surprisingly large flux ratio (Delta J =1.79 mag) considering the similar spectral types of its components. The large flux ratio could be explained if the primary is itself an equal-luminosity binary, which would make it one of the first known triple brown dwarf systems. Adaptive optics observations could not confirm this hypothesis, but it remains a likely one, which may be verified by high-resolution near-infrared spectroscopy. These two systems add to the handful of known brown dwarf binaries amenable to resolved spectroscopy without the aid of adaptive optics and constitute prime targets to test brown dwarf atmosphere models.

We present a detailed analysis of high resolution near-infrared imaging and spectroscopy of the potential star cluster IRS13E very close to the massive black hole in the Galactic Center. We detect 19 objects in IRS13E from Ks-band images, 15 of which are also detected reliably in H-band. We derive consistent proper motions for these objects from the two bands. Most objects share a similar westward proper motion. We characterize the objects using spectroscopy (1.45 to 2.45 micrometer) and (narrow-band) imaging from H- (1.66 mircrometer) to L'-band (3.80 micrometer). Nine of the objects detected in both Ks- and H-band are very red, and we find that they are all consistent with being warm dust clumps. The dust emission may be caused by the colliding winds of the two Wolf-Rayet stars in the cluster. Three of the six detected stars do not share the motion or spectral properties of the three bright stars. This leaves only the three bright, early-type stars as potential cluster members. It is unlikely that these stars are a chance configuration. Assuming the presence of an IMBH, a mass of about 14000 solar masses follows from the velocities and positions of these three stars. However, our acceleration limits make such an IMBH nearly as unlikely as a chance occurrence of such a star association. Furthermore, there is no variable X-ray source in IRS13E despite the high density of dust and gas. Therefore, we conclude that is unlikely that IRS13E hosts a black hole massive enough to bind the three stars.

Here we show that the ~10 Myr Beta Pictoris system hosts a massive giant planet, Beta Pictoris b, located 8 to 15 AU from the star. This result confirms that gas giant planets form rapidly within disks and validates the use of disk structures as fingerprints of embedded planets. Among the few planets already imaged, Beta Pictoris b is the closest to its parent star. Its short period could allow recording the full orbit within 17 years.

EAGLE is the multi-object, spatially-resolved, near-IR spectrograph instrument concept for the E-ELT, relying on a distributed Adaptive Optics, so-called Multi Object Adaptive Optics. This paper presents the results of a phase A study. Using 84x84 actuator deformable mirrors, the performed analysis demonstrates that 6 laser guide stars and up to 5 natural guide stars of magnitude R<17, picked-up in a 7.3' diameter patrol field of view, allow us to obtain an overall performance in terms of Ensquared Energy of 35% in a 75x75 mas^2 spaxel at H band, whatever the target direction in the centred 5' science field for median seeing conditions. The computed sky coverage at galactic latitudes |b|~60 is close to 90%.

With the development of high contrast imaging techniques and infrared detectors, vast efforts have been devoted during the past decade to detect and characterize lighter, cooler and closer companions to nearby stars, and ultimately image new planetary systems. Complementary to other observing techniques (radial velocity, transit, micro-lensing, pulsar-timing), this approach has opened a new astrophysical window to study the physical properties and the formation mechanisms of brown dwarfs and planets. I here will briefly present the observing challenge, the different observing techniques, strategies and samples of current exoplanet imaging searches that have been selected in the context of the LAOG-Planet Imaging Surveys. I will finally describe the most recent results that led to the discovery of giant planets probably formed like the ones of our solar system, offering exciting and attractive perspectives for the future generation of deep imaging instruments.

We recently reported on the detection of a possible planetary-mass companion to Beta Pictoris at a projected separation of 8 AU from the star, using data taken in November 2003 with NaCo, the adaptive-optics system installed on the Very Large Telescope UT4. Eventhough no second epoch detection was available, there are strong arguments to favor a gravitationally bound companion rather than a background object. If confirmed and located at a physical separation of 8 AU, this young, hot (~1500 K), massive Jovian companion (~8 Mjup) would be the closest planet to its star ever imaged, could be formed via core-accretion, and could explain the main morphological and dynamical properties of the dust disk. Our goal was to return to Beta Pic five years later to obtain a second-epoch observation of the companion or, in case of a non-detection, constrain its orbit. Deep adaptive-optics L'-band direct images of Beta Pic and Ks-band Four-Quadrant-Phase-Mask (4QPM) coronagraphic images were recorded with NaCo in January and February 2009. We also use 4QPM data taken in November 2004. No point-like signal with the brightness of the companion candidate (apparent magnitudes L'=11.2 or Ks ~ 12.5) is detected at projected distances down to 6.5 AU from the star in the 2009 data. As expected, the non-detection does not allow to rule out a background object; however, we show that it is consistent with the orbital motion of a bound companion that got closer to the star since first observed in 2003 and that is just emerging from behind the star at the present epoch. We place strong constraints on the possible orbits of the companion and discuss future observing prospects.

Since the discovery of its dusty disk in 1984, Beta Pictoris has become the prototype of young early-type planetary systems, and there are now various indications that a massive Jovian planet is orbiting the star at ~ 10 AU. However, no planets have been detected around this star so far. Our goal was to investigate the close environment of Beta Pic, searching for planetary companion(s). Deep adaptive-optics L'-band images of Beta Pic were recorded using the NaCo instrument at the Very Large Telescope. A faint point-like signal is detected at a projected distance of ~ 8 AU from the star, within the North-East side of the dust disk. Various tests were made to rule out with a good confidence level possible instrumental or atmospheric artifacts. The probability of a foreground or background contaminant is extremely low, based in addition on the analysis of previous deep Hubble Space Telescope images. The object L'=11.2 apparent magnitude would indicate a typical temperature of ~1500 K and a mass of ~ 8 Jovian masses. If confirmed, it could explain the main morphological and dynamical peculiarities of the Beta Pic system. The present detection is unique among A-stars by the proximity of the resolved planet to its parent star. Its closeness and location inside the Beta Pic disk suggest a formation process by core accretion or disk instabilities rather than a binary-like formation process.

Adaptive optics (AO) allows one to derive the point spread function (PSF) simultaneously to the science image, which is a major advantage in post-processing tasks such as astrometry/photometry or deconvolution. Based on the algorithm of \citet{veran97}, PSF reconstruction has been developed for four different AO systems so far: PUEO, ALFA, Lick-AO and Altair. A similar effort is undertaken for NAOS/VLT in a collaboration between the group PHASE (Onera and Observatoire de Paris/LESIA) and ESO. In this paper, we first introduce two new algorithms that prevent the use of the so-called "$U\_{ij}$ functions" to: (1) avoid the storage of a large amount of data (for both new algorithms), (2) shorten the PSF reconstruction computation time (for one of the two) and (3) provide an estimation of the PSF variability (for the other one). We then identify and explain issues in the exploitation of real-time Shack-Hartmann (SH) data for PSF reconstruction, emphasising the large impact of thresholding in the accuracy of the phase residual estimation. Finally, we present the data provided by the NAOS real-time computer (RTC) to reconstruct PSF ({\em (1)} the data presently available, {\em (2)} two NAOS software modifications that would provide new data to increase the accuracy of the PSF reconstruction and {\em (3)} the tests of these modifications) and the PSF reconstruction algorithms we are developing for NAOS on that basis.

We have collected in 2004 adaptive optics corrected L' images of the Galactic Center region with NAOS-CONICA at VLT. A strong variability was observed as well as a correlation between the photocenter positions and fluxes of the L'-band counterpart of Sgr A*. It is interpreted as the combined emission of a point-like flaring source at the position of Sgr A*/IR itself and an extended dust structure, 75 mas south west of Sgr A*/IR, which we name Sgr A*-f. We examine the different possible mechanisms to explain this dual Sgr A* L' emission and conclude it is likely a flaring emission associated to energetic events in the close environment of the black hole plus a quiescent emission resulting from the collision of Sgr A*-f by a jet from Sgr A*.

L'-band (lambda=3.8 microns) and M'-band (lambda=4.8 microns) observations of the Galactic Center region, performed in 2003 at VLT (ESO) with the adaptive optics imager NACO, have lead to the detection of an infrared counterpart of the radio source Sgr A* at both wavelengths. The measured fluxes confirm that the Sgr A* infrared spectrum is dominated by the synchrotron emission of nonthermal electrons. The infrared counterpart exhibits no significant short term variability but demonstrates flux variations on daily and yearly scales. The observed emission arises away from the position of the dynamical center of the S2 orbit and would then not originate from the closest regions of the black hole.

VLT J- to M\p-band adaptive optics observations of the Circinus Galaxy on parsec scales resolve a central bright Ks-band source with a FWHM size of 1.9 $\pm$ 0.6 pc. This source is only visible at wavelengths longward of 1.6 $\mu$m and coincides in position with the peak of the [Si VII]~2.48 $\mu$m coronal line emission. With respect to the peak of the central optical emission, the source is shifted by $\sim$ 0.15\arcsec (2.8 pc) to the south-east. Indeed, it defines the vertex of a fairly collimated beam which extends for $\sim$ 10 pc, and which is seen in both continuum light shortward of 1.6 $\mu$m and in H$\alpha$ line emission. The source also lies at the center of a $\sim$ 19 pc size [Si VII] ionization {\it bicone}. Identifying this source as the nucleus of Circinus, its size is compatible with a putative parsec-scale torus. Its spectral energy distribution, characterized by a prominent narrow peak, is compatible with a dust temperature of 300 K. Hotter dust within a 1 pc radius of the center is not detected. The AGN luminosity required to heat this dust is in the range of X-ray luminosities that have been measured toward the central source. This in turn supports the existence of highly obscuring material, with column densities of $10^{24}$ cm$^{-2}$, that must be located within 1 pc of the core.

We present K, L and M diffraction-limited images of NGC 1068 obtained with NAOS+CONICA at VLT/YEPUN over a 3.5" field around the central engine. Hot dust (Tcol = 550-650 K) is found in three different regions : (a) in the true nucleus, seen as a slightly NS elongated, core of extremely hot dust, "resolved" in K and L with respective diameters of ~5 pc and 8.5 pc ; (b) along the NS direction, as a "spiral arm" and a southern tongue ; (c) as a set of parallel elongated nodules ("wave-like") bracketting the jet. Several structures observed on radio maps, mid-IR or HST UV-visible maps are seen, so that a precise registration can be done from UV to 6 cm. These results do support the current interpretion that source (a) corresponds to emission from dust near sublimation temperature delimiting the walls of the cavity in the central obscuring torus. Structure (b) is thought to be a mixture of hot dust and active star forming regions along a micro spiral structure that could trace the tidal mechanism bringing matter to the central engine. Structure c)which was not known, exhibits too high a temperature for "classical'' grains ; it is most probably the signature of transiently heated very small dust grains (VSG) : "nano-diamonds", which are resistant and can form in strong UV field or in shocks, are very attractive candidates. The "waves'' can be condensations triggered by jet induced shocks, as predicted by recent models. First estimates, based on a simple VSG model and on a detailed radiative transfer model, do agree with those interpretations, both qualitatively and quantitatively.