AnisoCADO is a Python package for generating images of the point spread function (PSF) for the european extremely large telescope (ELT). The code allows the user to set many of the most important atmospheric and observational parameters that influence the shape and strehl ratio of the resulting PSF, including but not limited to: the atmospheric turbulence profile, the guide star position for a single conjugate adaptive optics (SCAO) solution, differential telescope pupil transmission, etc. Documentation can be found at https://anisocado.readthedocs.io/en/latest/

The Galactic Center black hole Sagittarius A* is a variable NIR source that exhibits bright flux excursions called flares. The low-flux density turnover of the flux distribution is below the sensitivity of current single-aperture telescopes. We use the unprecedented resolution of the GRAVITY instrument at the VLTI. Our light curves are unconfused, overcoming the confusion limit of previous photometric studies. We analyze the light curves using standard statistical methods and obtain the flux distribution. We find that the flux distribution of SgrA* turns over at a median flux density of (1.1\pm0.3)mJy. We measure the percentiles of the flux distribution and use them to constrain the NIR K-band SED. Furthermore, we find that the flux distribution is intrinsically right-skewed to higher flux density in log space. Flux densities below 0.1mJy are hardly ever observed. In consequence, a single powerlaw or lognormal distribution does not suffice to describe the observed flux distribution in its entirety. However, if one takes into account a power law component at high flux densities, a lognormal distribution can describe the lower end of the observed flux distribution. We confirm the RMS-flux relation for Sgr~A* and find it to be linear for all flux densities in our observation. We conclude that Sgr~A* has two states: the bulk of the emission is generated in a lognormal process with a well-defined median flux density and this quiescent emission is supplemented by sporadic flares that create the observed power law extension of the flux distribution.

The star S2 orbiting the compact radio source Sgr A* is a precision probe of the gravitational field around the closest massive black hole (candidate). Over the last 2.7 decades we have monitored the star's radial velocity and motion on the sky, mainly with the SINFONI and NACO adaptive optics (AO) instruments on the ESO VLT, and since 2017, with the four-telescope interferometric beam combiner instrument GRAVITY. In this paper we report the first detection of the General Relativity (GR) Schwarzschild Precession (SP) in S2's orbit. Owing to its highly elliptical orbit (e = 0.88), S2's SP is mainly a kink between the pre-and post-pericentre directions of motion ~ +- 1 year around pericentre passage, relative to the corresponding Kepler orbit. The superb 2017-2019 astrometry of GRAVITY defines the pericentre passage and outgoing direction. The incoming direction is anchored by 118 NACO-AO measurements of S2's position in the infrared reference frame, with an additional 75 direct measurements of the S2-Sgr A* separation during bright states ('flares') of Sgr A*. Our 14-parameter model fits for the distance, central mass, the position and motion of the reference frame of the AO astrometry relative to the mass, the six parameters of the orbit, as well as a dimensionless parameter f_SP for the SP (f_SP = 0 for Newton and 1 for GR). From data up to the end of 2019 we robustly detect the SP of S2, del phi = 12' per orbital period. From posterior fitting and MCMC Bayesian analysis with different weighting schemes and bootstrapping we find f_SP = 1.10 +- 0.19. The S2 data are fully consistent with GR. Any extended mass inside S2's orbit cannot exceed ~ 0.1% of the central mass. Any compact third mass inside the central arcsecond must be less than about 1000 M_sun.

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.

The inner regions of the discs of high-mass young stellar objects (HMYSOs) are still poorly known due to the small angular scales and the high visual extinction involved. We deploy near-infrared (NIR) spectro-interferometry to probe the inner gaseous disc in HMYSOs and investigate the origin and physical characteristics of the CO bandhead emission (2.3-2.4 $\mu$m). We present the first GRAVITY/VLTI observations at high spectral (R=4000) and spatial (mas) resolution of the CO overtone transitions in NGC 2024 IRS2. The continuum emission is resolved in all baselines and is slightly asymmetric, displaying small closure phases ($\leq$8$^{\circ}$). Our best ellipsoid model provides a disc inclination of 34$^{\circ}$$\pm$1$^{\circ}$, a disc major axis position angle of 166$^{\circ}$$\pm$1$^{\circ}$, and a disc diameter of 3.99$\pm$0.09 mas (or 1.69$\pm$0.04 au, at a distance of 423 pc). The small closure phase signals in the continuum are modelled with a skewed rim, originating from a pure inclination effect. For the first time, our observations spatially and spectrally resolve the first four CO bandheads. Changes in visibility, as well as differential and closure phases across the bandheads are detected. Both the size and geometry of the CO-emitting region are determined by fitting a bidimensional Gaussian to the continuum-compensated CO bandhead visibilities. The CO-emitting region has a diameter of 2.74$\pm^{0.08}_{0.07}$ mas (1.16$\pm$0.03 au), and is located in the inner gaseous disc, well within the dusty rim, with inclination and $PA$ matching the dusty disc geometry, which indicates that both dusty and gaseous discs are coplanar. Physical and dynamical gas conditions are inferred by modelling the CO spectrum. Finally, we derive a direct measurement of the stellar mass of $M_*\sim$14.7$^{+2}_{-3.6}$ M$_{\odot}$ by combining our interferometric and spectral modelling results.

Infrared observations of Sgr A* probe the region close to the event horizon of the black hole at the Galactic center. These observations can constrain the properties of low-luminosity accretion as well as that of the black hole itself. The GRAVITY instrument at the ESO VLTI has recently detected continuous circular relativistic motion during infrared flares which has been interpreted as orbital motion near the event horizon. Here we analyze the astrometric data from these flares, taking into account the effects of out-of-plane motion and orbital shear of material near the event horizon of the black hole. We have developed a new code to predict astrometric motion and flux variability from compact emission regions following particle orbits. Our code combines semi-analytic calculations of timelike geodesics that allow for out-of-plane or elliptical motions with ray tracing of photon trajectories to compute time-dependent images and light curves. We apply our code to the three flares observed with GRAVITY in 2018. We show that all flares are consistent with a hotspot orbiting at R$\sim$9 gravitational radii with an inclination of $i\sim140^\circ$. The emitting region must be compact and less than $\sim5$ gravitational radii in diameter. We place a further limit on the out-of-plane motion during the flare.

Our objective is to estimate the C/O ratio in the atmosphere of beta Pictoris b and obtain an estimate of the dynamical mass of the planet, as well as to refine its orbital parameters using high-precision astrometry. We used the GRAVITY instrument with the four 8.2 m telescopes of the Very Large Telescope Interferometer to obtain K-band spectro-interferometric data on $\beta$ Pic b. We extracted a medium resolution (R=500) K-band spectrum of the planet and a high-precision astrometric position. We estimated the planetary C/O ratio using two different approaches (forward modeling and free retrieval) from two different codes (ExoREM and petitRADTRANS, respectively). Finally, we used a simplified model of two formation scenarios (gravitational collapse and core-accretion) to determine which can best explain the measured C/O ratio. Our new astrometry disfavors a circular orbit for $\beta$ Pic b ($e=0.15^{+0.05}_{-0.04}$). Combined with previous results and with Hipparcos/GAIA measurements, this astrometry points to a planet mass of $M = 12.7\pm{}2.2\,M_\mathrm{Jup}$. This value is compatible with the mass derived with the free-retrieval code petitRADTRANS using spectral data only. The forward modeling and free-retrieval approches yield very similar results regarding the atmosphere of beta Pic b. In particular, the C/O ratios derived with the two codes are identical ($0.43\pm{}0.05$ vs $0.43^{+0.04}_{-0.03}$). We argue that if the stellar C/O in $\beta$ Pic is Solar, then this combination of a very high mass and a low C/O ratio for the planet suggests a formation through core-accretion, with strong planetesimal enrichment.

Mapping of the large-scale structure through cosmic time has numerous applications in the studies of cosmology and galaxy evolution. At $z > 2$, the structure can be traced by the neutral intergalactic medium (IGM) by way of observing the Ly$\alpha$, forest towards densely-sampled lines-of-sight of bright background sources, such as quasars and star forming galaxies. We investigate the scientific potential of MOSAIC, a planned multi-object spectrograph on the European Extremely Large Telescope (ELT), for the 3D mapping of the IGM at $z \gtrsim 3$. We simulate a survey of $3 \lesssim z \lesssim 4$ galaxies down to a limiting magnitude of $m_{r}\sim 25.5$ mag in an area of 1 degree$^2$ in the sky. Galaxies and their spectra (including the line-of-sight Ly$\alpha$ absorption) are taken from the lightcone extracted from the Horizon-AGN cosmological hydrodynamical simulation. The quality of the reconstruction of the original density field is studied for different spectral resolutions and signal-to-noise ratios of the spectra. We demonstrate that the minimum $S/N$ (per resolution element) of the faintest galaxies that such survey has to reach is $S/N = 4$. We show that a survey with such sensitivity enables a robust extraction of cosmic filaments and the detection of the theoretically-predicted galaxy stellar mass and star-formation rate gradients towards filaments. By simulating the realistic performance of MOSAIC we obtain $S/N(T_{\rm obs}, R, m_{r})$ scaling relations. We estimate that $\lesssim 35~(65)$ nights of observation time are required to carry out the survey with the instrument's high multiplex mode and with the spectral resolution of $R=1000~(2000)$. A survey with a MOSAIC-concept instrument on the ELT is found to enable the mapping of the IGM at $z > 3$ on Mpc scales, and as such will be complementary to and competitive with other planned IGM tomography surveys. [abridged]

The formation and the evolution of protoplanetary disks are important stages in the lifetime of stars. The processes of disk evolution and planet formation are intrinsically linked. We spatially resolve with GRAVITY/VLTI in the K-band the sub au-scale region of 27 stars to gain statistical understanding of their properties. We look for correlations with stellar parameters, such as luminosity, mass, temperature and age. Our sample also cover a range of various properties in terms of reprocessed flux, flared or flat morphology, and gaps. We developed semi-physical geometrical models to fit our interferometric data. Our best models correspond to smooth and wide rings, implying that wedge-shaped rims at the dust sublimation edge are favored, as found in the H-band. The closure phases are generally non-null with a median value of ~10 deg, indicating spatial asymmetries of the intensity distributions. Multi-size grain populations could explain the closure phase ranges below 20-25 deg but other scenarios should be invoked to explain the largest ones. Our measurements extend the Radius-Luminosity relation to ~1e4 Lsun and confirm the significant spread around the mean relation observed in the H-band. Gapped sources exhibit a large N-to-K band size ratio and large values of this ratio are only observed for the members of our sample that would be older than 1 Ma, less massive, and with lower luminosity. In the 2 Ms mass range, we observe a correlation in the increase of the relative age with the transition from group II to group I, and an increase of the N-to-K size ratio. However, the size of the current sample does not yet permit us to invoke a clear universal evolution mechanism across the HAeBe mass range. The measured locations of the K-band emission suggest that these disks might be structured by forming young planets, rather than by depletion due to EUV, FUV, and X-ray photo-evaporation.

Precise measurements of the S-stars orbiting SgrA* have set strong constraints on the nature of the compact object at the centre of the Milky Way. The presence of a black hole in that region is well established, but its neighboring environment is still an open debate. In that respect, the existence of dark matter in that central region may be detectable due to its strong signatures on the orbits of stars: the main effect is a Newtonian precession which will affect the overall pericentre shift of S2, the latter being a target measurement of the GRAVITY instrument. The exact nature of this dark matter (e.g., stellar dark remnants or diffuse dark matter) is unknown. This article assumes it to be an scalar field of toroidal distribution, associated with ultra-light dark matter particles, surrounding the Kerr black hole. Such a field is a form of "hair" expected in the context of superradiance, a mechanism that extracts rotational energy from the black hole. Orbital signatures for the S2 star are computed and shown to be detectable by GRAVITY. The scalar field can be constrained because the variation of orbital elements depends both on the relative mass of the scalar field to the black hole and on the field mass coupling parameter.

We present a 0.16% precise and 0.27% accurate determination of R0, the distance to the Galactic Center. Our measurement uses the star S2 on its 16-year orbit around the massive black hole Sgr A* that we followed astrometrically and spectroscopically for 27 years. Since 2017, we added near-infrared interferometry with the VLTI beam combiner GRAVITY, yielding a direct measurement of the separation vector between S2 and Sgr A* with an accuracy as good as 20 micro-arcsec in the best cases. S2 passed the pericenter of its highly eccentric orbit in May 2018, and we followed the passage with dense sampling throughout the year. Together with our spectroscopy, in the best cases with an error of 7 km/s, this yields a geometric distance estimate: R0 = 8178 +- 13(stat.) +- 22(sys.) pc. This work updates our previous publication in which we reported the first detection of the gravitational redshift in the S2 data. The redshift term is now detected with a significance level of 20 sigma with f_redshift = 1.04 +- 0.05.

To date, infrared interferometry at best achieved contrast ratios of a few times $10^{-4}$ on bright targets. GRAVITY, with its dual-field mode, is now capable of high contrast observations, enabling the direct observation of exoplanets. We demonstrate the technique on HR8799, a young planetary system composed of four known giant exoplanets. We used the GRAVITY fringe tracker to lock the fringes on the central star, and integrated off-axis on the HR8799e planet situated at 390 mas from the star. Data reduction included post-processing to remove the flux leaking from the central star and to extract the coherent flux of the planet. The inferred K band spectrum of the planet has a spectral resolution of 500. We also derive the astrometric position of the planet relative to the star with a precision on the order of 100$\,\mu$as. The GRAVITY astrometric measurement disfavors perfectly coplanar stable orbital solutions. A small adjustment of a few degrees to the orbital inclination of HR 8799 e can resolve the tension, implying that the orbits are close to, but not strictly coplanar. The spectrum, with a signal-to-noise ratio of $\approx 5$ per spectral channel, is compatible with a late-type L brown dwarf. Using Exo-REM synthetic spectra, we derive a temperature of $1150\pm50$\,K and a surface gravity of $10^{4.3\pm0.3}\,$cm/s$^{2}$. This corresponds to a radius of $1.17^{+0.13}_{-0.11}\,R_{\rm Jup}$ and a mass of $10^{+7}_{-4}\,M_{\rm Jup}$, which is an independent confirmation of mass estimates from evolutionary models. Our results demonstrate the power of interferometry for the direct detection and spectroscopic study of exoplanets at close angular separations from their stars.

During its orbit around the four million solar mass black hole Sagittarius A* the star S2 experiences significant changes in gravitational potential. We use this change of potential to test one part of the Einstein equivalence principle: the local position invariance (LPI). We study the dependency of different atomic transitions on the gravitational potential to give an upper limit on violations of the LPI. This is done by separately measuring the redshift from hydrogen and helium absorption lines in the stellar spectrum during its closest approach to the black hole. For this measurement we use radial velocity data from 2015 to 2018 and combine it with the gravitational potential at the position of S2, which is calculated from the precisely known orbit of S2 around the black hole. This results in a limit on a violation of the LPI of $|\beta_{He}-\beta_{H}| = (2.4 \pm 5.1) \cdot 10^{-2}$. The variation in potential that we probe with this measurement is six magnitudes larger than possible for measurements on Earth, and a factor ten larger than in experiments using white dwarfs. We are therefore testing the LPI in a regime where it has not been tested before.

Advanced adaptive optics (AO) instruments on ground-based telescopes require accurate knowledge of the atmospheric turbulence strength as a function of altitude. This information assists point spread function reconstruction, AO temporal control techniques and is required by wide-field AO systems to optimize the reconstruction of an observed wavefront. The variability of the atmosphere makes it important to have a measure of the optical turbulence profile in real time. This measurement can be performed by fitting an analytically generated covariance matrix to the cross-covariance of Shack-Hartmann wavefront sensor (SHWFS) centroids. In this study we explore the benefits of reducing cross-covariance data points to a covariance map region of interest (ROI). A technique for using the covariance map ROI to measure and compensate for SHWFS misalignments is also introduced. We compare the accuracy of covariance matrix and map ROI optical turbulence profiling using both simulated and on-sky data from CANARY, an AO demonstrator on the 4.2 m William Herschel telescope, La Palma. On-sky CANARY results are compared to contemporaneous profiles from Stereo-SCIDAR - a dedicated high-resolution optical turbulence profiler. It is shown that the covariance map ROI optimizes the accuracy of AO telemetry optical turbulence profiling. In addition, we show that the covariance map ROI reduces the fitting time for an extremely large telescope-scale system by a factor of 72. The software package we developed to collect all of the presented results is now open source.

The European Extremely Large Telescope will see first lights by the end of 2024. With a diameter of almost 40 meters, it will be the biggest optical telescope ever built from the ground. The ELT will open a brand new window in a sensitivity / spatial angular resolution parameter space. To take the full benefit of the scientific potential of this giant, all the instruments will be equipped with Adaptive Optics (AO), providing the sharpest images. This paper provides a quick overview of the AO capabilities of the future instruments to be deployed at the ELT, highlight some of the expected performance and describe a couple of technical challenges that are still to tackle for an optimal scientific use. This paper has been presented at the "Societe Francaise d'Astronomie et Astrophysique" symposium in Bordeaux 2018, it is then naturally biased toward the French contribution for the ELT.

We report the detection of continuous positional and polarization changes of the compact source SgrA* in high states ('flares') of its variable near- infrared emission with the near-infrared GRAVITY-Very Large Telescope Interferometer (VLTI) beam-combining instrument. In three prominent bright flares, the position centroids exhibit clockwise looped motion on the sky, on scales of typically 150 micro-arcseconds over a few tens of minutes, corresponding to about 30% the speed of light. At the same time, the flares exhibit continuous rotation of the polarization angle, with about the same 45(+/-15)-minute period as that of the centroid motions. Modelling with relativistic ray tracing shows that these findings are all consistent with a near face-on, circular orbit of a compact polarized 'hot spot' of infrared synchrotron emission at approximately six to ten times the gravitational radius of a black hole of 4 million solar masses. This corresponds to the region just outside the innermost, stable, prograde circular orbit (ISCO) of a Schwarzschild-Kerr black hole, or near the retrograde ISCO of a highly spun-up Kerr hole. The polarization signature is consistent with orbital motion in a strong poloidal magnetic field.

On 6 Sept. 2017, an exceptional coronal mass ejection departed from the Sun toward the Earth and Uranus, whose magnetospheres are sensitive to the solar wind. The resulting interplanetary shock triggered geomagnetic storm and intense aurora at Earth the next day and was predicted by MHD models to reach Uranus around 10-11 Nov. This event provided a unique opportunity to investigate the auroral response of the asymmetric Uranian magnetosphere in its intermediate equinox-to-solstice configuration. Coordinated multi-spectral observations were acquired with the Hubble Space Telescope (HST) in the far-UV (FUV), with the Very Large Telescope (VLT) and Gemini North in the near-IR (NIR) and with Chandra in the X-ray domain. In this study, we focus on the analysis of NIR images obtained between 9 and 17 Nov. 2017 which are compared to one FUV image acquired on 11 Nov. The latter reveals a bright southern auroral spot in the H2 bands, which we use as a reference to locate auroral precipitations. The NIR images were aimed at mapping H3+ emission from the Uranian ionosphere and at updating the results built from a couple of pioneer images taken 25 years ago. These new high resolution images reveal H3+ from the whole disc although brighter near the southern pole, but show no evidence of localized auroral emission.

Eta Car is one of the most intriguing luminous blue variables in the Galaxy. Observations and models at different wavelengths suggest a central binary with a 5.54 yr period residing in its core. 2D and 3D radiative transfer and hydrodynamic simulations predict a primary with a dense and slow stellar wind that interacts with the faster and lower density wind of the secondary. The wind-wind collision scenario suggests that the secondary's wind penetrates the primary's wind creating a low-density cavity in it, with dense walls where the two winds interact. We aim to trace the inner ~5-50 au structure of Eta Car's wind-wind interaction, as seen through BrG and, for the first time, through the He I 2s-2p line. We have used spectro-interferometric observations with GRAVITY at the VLTI. Our modeling of the continuum allows us to estimate its FWHM angular size close to 2 mas and an elongation ratio of 1.06 +/- 0.05 over a PA = 130 +/- 20 deg. Our CMFGEN modeling helped us to confirm that the role of the secondary should be taken into account to properly reproduce the observed BrG and He I lines. Chromatic images across BrG reveal a southeast arc-like feature, possibly associated to the hot post-shocked winds flowing along the cavity wall. The images of He I 2s-2p served to constrain the 20 mas structure of the line-emitting region. The observed morphology of He I suggests that the secondary is responsible for the ionized material that produces the line profile. Both the BrG and the He I 2s-2p maps are consistent with previous hydrodynamical models of the colliding wind scenario. Future dedicated simulations together with an extensive interferometric campaign are necessary to refine our constraints on the wind and stellar parameters of the binary, which finally will help us predict the evolutionary path of Eta Car.

MICADO will enable the ELT to perform diffraction limited near-infrared observations at first light. The instrument's capabilities focus on imaging (including astrometric and high contrast) as well as single object spectroscopy. This contribution looks at how requirements from the observing modes have driven the instrument design and functionality. Using examples from specific science cases, and making use of the data simulation tool, an outline is presented of what we can expect the instrument to achieve.

The highly elliptical, 16-year-period orbit of the star S2 around the massive black hole candidate Sgr A* is a sensitive probe of the gravitational field in the Galactic centre. Near pericentre at 120 AU, ~1400 Schwarzschild radii, the star has an orbital speed of ~7650 km/s, such that the first-order effects of Special and General Relativity have now become detectable with current capabilities. Over the past 26 years, we have monitored the radial velocity and motion on the sky of S2, mainly with the SINFONI and NACO adaptive optics instruments on the ESO Very Large Telescope, and since 2016 and leading up to the pericentre approach in May 2018, with the four-telescope interferometric beam-combiner instrument GRAVITY. From data up to and including pericentre, we robustly detect the combined gravitational redshift and relativistic transverse Doppler effect for S2 of z ~ 200 km/s / c with different statistical analysis methods. When parameterising the post-Newtonian contribution from these effects by a factor f, with f = 0 and f = 1 corresponding to the Newtonian and general relativistic limits, respectively, we find from posterior fitting with different weighting schemes f = 0.90 +/- 0.09 (stat) +\- 0.15 (sys). The S2 data are inconsistent with pure Newtonian dynamics.

The paper discusses an assessment study about the impact of the distortions on the astrometric observations with the Extremely Large Telescope originated from the optics positioning errors and telescope instabilities. Optical simulations combined with Monte Carlo approach reproducing typical inferred opto-mechanical and dynamical instabilities, show RMS distortions between $\sim$ 0.1-5 mas over 1 arcmin field of view. Over minutes timescales the plate scale variations from ELT-M2 caused by wind disturbances and gravity flexures and the field rotation from ELT-M4-M5 induce distortions and PSF jitter at the edge of 1 arcmin FoV (radius 35 arcsec) up to $\sim$ 5 mas comparable to the diffraction-limited PSF size $FWHM_H = 8.5$ mas. The RMS distortions inherent to the ELT design are confined to the 1$^{st}$-3$^{rd}$ order and reduce to an astrometric RMS residual post fit of $\sim$ 10-20 $\mu as$ for higher order terms. In this paper, we study which calibration effort has to be undertaken to reach an astrometric stability close to this level of higher order residuals. The amplitude and timescales of the assumed telescope tolerances indicate the need for frequent on-sky calibrations and MCAO stabilization of the plate scale to enable astrometric observations with ELT at the level of $\leq 50 \mu as$, which is one of the core science missions for the ELT / MICADO instrument.

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.

To investigate the inner regions of protoplanetary disks, we performed near-infrared interferometric observations of the classical TTauri binary system S CrA. We present the first VLTI-GRAVITY high spectral resolution ($R\sim$4000) observations of a classical TTauri binary, S CrA (composed of S CrA N and S CrA S and separated by $\sim$1.4"), combining the four 8-m telescopes in dual-field mode. Our observations in the near-infrared K-band continuum reveal a disk around each binary component, with similar half-flux radii of about 0.1 au at d$\sim$130 pc, inclinations ($i=$28$\pm$3$^o$\ and $i=$22$\pm$6$^o$), and position angles (PA=0$^o\pm$6$^o$ and PA=-2$^o\pm$12$^o$), suggesting that they formed from the fragmentation of a common disk. The S CrA N spectrum shows bright HeI and Br$\gamma$ line emission exhibiting inverse P-Cygni profiles, typically associated with infalling gas. The continuum-compensated Br$\gamma$ line visibilities of S CrA N show the presence of a compact Br$\gamma$ emitting region the radius of which is about $\sim$0.06 au, which is twice as big as the truncation radius. This component is mostly tracing a wind. Moreover, a slight radius change between the blue- and red-shifted Br$\gamma$ line components is marginally detected. The presence of an inverse P-Cygni profile in the HeI and Br$\gamma$ lines, along with the tentative detection of a slightly larger size of the blue-shifted Br$\gamma$ line component, hint at the simultaneous presence of a wind and magnetospheric accretion in S CrA N.

GRAVITY is a new instrument to coherently combine the light of the European Southern Observatory Very Large Telescope Interferometer to form a telescope with an equivalent 130 m diameter angular resolution and a collecting area of 200 m$^2$. The instrument comprises fiber fed integrated optics beam combination, high resolution spectroscopy, built-in beam analysis and control, near-infrared wavefront sensing, phase-tracking, dual beam operation and laser metrology [...]. This article gives an overview of GRAVITY and reports on the performance and the first astronomical observations during commissioning in 2015/16. We demonstrate phase tracking on stars as faint as m$_K$ ~ 10 mag, phase-referenced interferometry of objects fainter than m$_K$ ~ 15 mag with a limiting magnitude of m$_K$ ~ 17 mag, minute long coherent integrations, a visibility accuracy of better than 0.25 %, and spectro-differential phase and closure phase accuracy better than 0.5{\deg}, corresponding to a differential astrometric precision of better than 10 microarcseconds ({\mu}as). The dual-beam astrometry, measuring the phase difference of two objects with laser metrology, is still under commissioning. First observations show residuals as low as 50 {\mu}as when following objects over several months. We illustrate the instrument performance with the observations of archetypical objects for the different instrument modes. Examples include the Galactic Center supermassive black hole and its fast orbiting star S2 for phase referenced dual beam observations and infrared wavefront sensing, the High Mass X-Ray Binary BP Cru and the Active Galactic Nucleus of PDS 456 for few {\mu}as spectro-differential astrometry, the T Tauri star S CrA for a spectro-differential visibility analysis, {\xi} Tel and 24 Cap for high accuracy visibility observations, and {\eta} Car for interferometric imaging with GRAVITY.

We observe the HMXB BP Cru using interferometry in the near-infrared K band with VLTI/GRAVITY. Continuum visibilities are at most partially resolved, consistent with the predicted size of the hypergiant. Differential visibility amplitude ($\Delta |V| \sim 5\%$) and phase ($\Delta \phi \sim 2 \deg$) signatures are observed across the HeI $2.059 \mu$m and Br$\gamma$ lines, the latter seen strongly in emission, unusual for the donor star's spectral type. For a baseline $B \sim 100$m, the differential phase RMS $\sim 0.2 \deg$ corresponds to an astrometric precision of $\sim 2 \mu$as. A model-independent analysis in the marginally resolved limit of interferometry reveals asymmetric and extended emission with a strong wavelength dependence. We propose geometric models based on an extended and distorted wind and/or a high density gas stream, which has long been predicted to be present in this system. The observations show that optical interferometry is now able to resolve HMXBs at the spatial scale at which accretion takes place, and therefore probe the effects of the gravitational and radiation fields of the compact object on its environment.

We present the first optical observation at sub-milliarcsecond (mas) scale of the microquasar SS 433 obtained with the GRAVITY instrument on the VLT interferometer. The 3.5 hour exposure reveals a rich K-band spectrum dominated by hydrogen Br$\gamma\ $ and \ion{He}{i} lines, as well as (red-shifted) emission lines coming from the jets. The K-band continuum emitting region is dominated by a marginally resolved point source ($<$ 1 mas) embedded inside a diffuse background accounting for 10\% of the total flux. The jet line positions agree well with the ones expected from the jet kinematic model, an interpretation also supported by the consistent sign (i.e. negative/positive for the receding/approaching jet component) of the phase shifts observed in the lines. The significant visibility drop across the jet lines, together with the small and nearly identical phases for all baselines, point toward a jet that is offset by less than 0.5 mas from the continuum source and resolved in the direction of propagation, with a typical size of 2 mas. The jet position angle of $\sim$80$^{\circ}$ is consistent with the expected one at the observation date. Jet emission so close to the central binary system would suggest that line locking, if relevant to explain the amplitude and stability of the 0.26c jet velocity, operates on elements heavier than hydrogen. The Br$\gamma\ $ profile is broad and double peaked. It is better resolved than the continuum and the change of the phase signal sign across the line on all baselines suggests an East-West oriented geometry alike the jet direction and supporting a (polar) disk wind origin.

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) pathfinder for the future MOAO-assisted Integral-Field Units (IFU) proposed for Extremely Large Telescopes (ELT). The MOAO concept relies on tomographically reconstructing the turbulence using multiple measurements along different lines of sight. Tomography requires the knowledge of the statistical turbulence parameters, commonly recovered from the system telemetry using a dedicated profiling technique. For demonstration purposes with the MOAO pathfinder Canary , this identification is performed thanks to the Learn & Apply (L&A) algorithm, that consists in model- fitting the covariance matrix of WFS measurements dependent on relevant parameters: $C_n^2(h)$ profile, outer scale profile and system mis-registration. We explore an upgrade of this algorithm, the Learn 3 Steps (L3S) approach, that allows one to dissociate the identification of the altitude layers from the ground in order to mitigate the lack of convergence of the required empirical covariance matrices therefore reducing the required length of data time-series for reaching a given accuracy. For nominal observation conditions, the L3S can reach the same level of tomographic error in using five times less data frames than the L&A approach. The L3S technique has been applied over a large amount of Canary data to characterize the turbu- lence above the William Herschel Telescope (WHT). These data have been acquired the 13th, 15th, 16th, 17th and 18th September 2013 and we find 0.67"/8.9m/3.07m/s of total seeing/outer scale/wind-speed, with 0.552"/9.2m/2.89m/s below 1.5 km and 0.263"/10.3m/5.22m/s between 1.5 and 20 km. We have also de- termined the high altitude layers above 20 km, missed by the tomographic reconstruction on Canary , have a median seeing of 0.187" and have occurred 16% of observation time.

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.

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.

The spectral resolution of a dispersive astronomical spectrograph is limited by the trade-off between throughput and the width of the entrance slit. Photonic guided-wave transitions have been proposed as a route to bypass this trade-off, by enabling the efficient reformatting of incoherent seeing-limited light collected by the telescope into a linear array of single modes: a pseudo-slit which is highly multimode in one axis but diffraction-limited in the dispersion axis of the spectrograph. It is anticipated that the size of a single-object spectrograph fed with light in this manner would be essentially independent of the telescope aperture size. A further anticipated benefit is that such spectrographs would be free of `modal noise', a phenomenon that occurs in high-resolution multimode fibre-fed spectrographs due to the coherent nature of the telescope point-spread-function (PSF). We address these aspects by integrating a multicore fibre photonic lantern with an ultrafast laser inscribed three-dimensional waveguide interconnect to spatially reformat the modes within the PSF into a diffraction-limited pseudo-slit. Using the CANARY adaptive optics (AO) demonstrator on the William Herschel Telescope, and 1530 $\pm$ 80 nm stellar light, the device is found to exhibit a transmission of 47-53 % depending upon the mode of AO correction applied. We also show the advantage of using AO to couple light into such a device by sampling only the core of the CANARY PSF. This result underscores the possibility that a fully-optimised guided-wave device can be used with AO to provide efficient spectroscopy at high spectral resolution.

The spectral resolution of a dispersive spectrograph is dependent on the width of the entrance slit. This means that astronomical spectrographs trade-off throughput with spectral resolving power. Recently, optical guided-wave transitions known as photonic lanterns have been proposed to circumvent this trade-off, by enabling the efficient reformatting of multimode light into a pseudo-slit which is highly multimode in one axis, but diffraction-limited in the other. Here, we demonstrate the successful reformatting of a telescope point spread function into such a slit using a three-dimensional integrated optical waveguide device, which we name the photonic dicer. Using the CANARY adaptive optics demonstrator on the William Herschel Telescope, and light centred at 1530 nm with a 160 nm FWHM, the device shows a transmission of between 10 and 20\% depending upon the type of AO correction applied. Most of the loss is due to the overfilling of the input aperture in poor and moderate seeing. Taking this into account, the photonic device itself has a transmission of 57 $\pm$ 4\%.We show how a fully-optimised device can be used with AO to provide efficient spectroscopy at high spectral resolution.

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.

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

GRAVITY is a second generation instrument for the VLT Interferometer, designed to enhance the near-infrared astrometric and spectro-imaging capabilities of VLTI. Combining beams from four telescopes, GRAVITY will provide an astrometric precision of order 10 micro-arcseconds, imaging resolution of 4 milli-arcseconds, and low and medium resolution spectro-interferometry, pushing its performance far beyond current infrared interfero- metric capabilities. To maximise the performance of GRAVITY, adaptive optics correction will be implemented at each of the VLT Unit Telescopes to correct for the effects of atmospheric turbulence. To achieve this, the GRAVITY project includes a development programme for four new wavefront sensors (WFS) and NIR-optimized real time control system. These devices will enable closed-loop adaptive correction at the four Unit Telescopes in the range 1.4-2.4 {\mu}m. This is crucially important for an efficient adaptive optics implementation in regions where optically bright references sources are scarce, such as the Galactic Centre. We present here the design of the GRAVITY wavefront sensors and give an overview of the expected adaptive optics performance under typical observing conditions. Benefiting from newly developed SELEX/ESO SAPHIRA electron avalanche photodiode (eAPD) detectors providing fast readout with low noise in the near-infrared, the AO systems are expected to achieve residual wavefront errors of \leq400 nm at an operating frequency of 500 Hz.

The EAGLE and EVE Phase A studies for instruments for the European Extremely Large Telescope (E-ELT) originated from related top-level scientific questions, but employed different (yet complementary) methods to deliver the required observations. We re-examine the motivations for a multi-object spectrograph (MOS) on the E-ELT and present a unified set of requirements for a versatile instrument. Such a MOS would exploit the excellent spatial resolution in the near-infrared envisaged for EAGLE, combined with aspects of the spectral coverage and large multiplex of EVE. We briefly discuss the top-level systems which could satisfy these requirements in a single instrument at one of the Nasmyth foci of the E-ELT.

Astronomical adaptive optics systems with open-loop deformable mirror control have recently come on-line. In these systems, the deformable mirror surface is not included in the wavefront sensor paths, and so changes made to the deformable mirror are not fed back to the wavefront sensors. This gives rise to all sorts of linearity and control issues mainly centred on one question: Has the mirror taken the shape requested? Non-linearities in wavefront measurement and in the deformable mirror shape can lead to significant deviations in mirror shape from the requested shape. Here, wavefront sensor measurements made using a brightest pixel selection method are discussed along with the implications that this has for open-loop AO systems. Discussion includes elongated laser guide star spots and also computational efficiency.

GRAVITY is an adaptive optics assisted Beam Combiner for the second generation VLTI instrumentation. The instrument will provide high-precision narrow-angle astrometry and phase-referenced interferometric imaging in the astronomical K-band for faint objects. We describe the wide range of science that will be tackled with this instrument, highlighting the unique capabilities of the VLTI in combination with GRAVITY. The most prominent goal is to observe highly relativistic motions of matter close to the event horizon of Sgr A*, the massive black hole at center of the Milky Way. We present the preliminary design that fulfils the requirements that follow from the key science drivers: It includes an integrated optics, 4-telescope, dual feed beam combiner operated in a cryogenic vessel; near-infrared wavefrontsensing adaptive optics; fringe-tracking on secondary sources within the field of view of the VLTI and a novel metrology concept. Simulations show that 10 {\mu}as astrometry within few minutes is feasible for a source with a magnitude of mK = 15 like Sgr A*, given the availability of suitable phase reference sources (mK = 10). Using the same setup, imaging of mK = 18 stellar sources in the interferometric field of view is possible, assuming a full night of observations and the corresponding UV coverage of the VLTI.

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

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.

The very large telescope (VLT) interferometer (VLTI) in its current operating state is equipped with high-order adaptive optics (MACAO) working in the visible spectrum. A low-order near-infrared wavefront sensor (IRIS) is available to measure non-common path tilt aberrations downstream the high-order deformable mirror. For the next generation of VLTI instrumentation, in particular for the designated GRAVITY instrument, we have examined various designs of a four channel high-order near-infrared wavefront sensor. Particular objectives of our study were the specification of the near-infrared detector in combination with a standard wavefront sensing system. In this paper we present the preliminary design of a Shack-Hartmann wavefront sensor operating in the near-infrared wavelength range, which is capable of measuring the wavefronts of four telescopes simultaneously. We further present results of our design study, which aimed at providing a first instrumental concept for GRAVITY.

We present the second-generation VLTI instrument GRAVITY, which currently is in the preliminary design phase. GRAVITY is specifically designed to observe highly relativistic motions of matter close to the event horizon of Sgr A*, the massive black hole at center of the Milky Way. We have identified the key design features needed to achieve this goal and present the resulting instrument concept. It includes an integrated optics, 4-telescope, dual feed beam combiner operated in a cryogenic vessel; near infrared wavefront sensing adaptive optics; fringe tracking on secondary sources within the field of view of the VLTI and a novel metrology concept. Simulations show that the planned design matches the scientific needs; in particular that 10 microarcsecond astrometry is feasible for a source with a magnitude of K=15 like Sgr A*, given the availability of suitable phase reference sources.

Integral field spectrographs are major instruments to study the mechanisms involved in the formation and the evolution of early galaxies. When combined with multi-object spectroscopy, those spectrographs can behave as machines used to derive physical parameters of galaxies during their formation process. Up to now, there is only one available spectrograph with multiple integral field units, e.g. FLAMES/GIRAFFE on the VLT. However, current ground based instruments suffer from a degradation of their spatial resolution due to atmospheric turbulence. In this article we describe the performance of FALCON, an original concept of a new generation multi-object integral field spectrograph with adaptive optics for the ESO Very Large Telescope. The goal of FALCON is to combine high angular resolution (0.25 arcsec) and high spectral resolution (R > 5000) in J and H bands over a wide field of view (10x10 arcmin2) in the VLT Nasmyth focal plane. However, instead of correcting the whole field, FALCON will use multi-object adaptive optics (MOAO) to perform locally on each scientific target the adaptive optics correction. This requires then to use atmospheric tomography in order to use suitable natural guide stars for wavefront sensing. We will show that merging MOAO and atmospheric tomography allows us to determine the internal kinematics of distant galaxies up to z=2 with a sky coverage of 50%, even for objects observed near the galactic pole. The application of such a concept to Extremely Large Telescopes seems therefore to be a very promising way to study galaxy evolution from z = 1 to redshifts as high as z = 7.

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.

Context. The knowledge of the point-spread function compensated by adaptive optics is of prime importance in several image restoration techniques such as deconvolution and astrometric/photometric algorithms. Wavefront-related data from the adaptive optics real-time computer can be used to accurately estimate the point-spread function in adaptive optics observations. The only point-spread function reconstruction algorithm implemented on astronomical adaptive optics system makes use of particular functions, named $U\_{ij}$. These $U\_{ij}$ functions are derived from the mirror modes, and their number is proportional to the square number of these mirror modes. Aims. We present here two new algorithms for point-spread function reconstruction that aim at suppressing the use of these $U\_{ij}$ functions to avoid the storage of a large amount of data and to shorten the computation time of this PSF reconstruction. Methods. Both algorithms take advantage of the eigen decomposition of the residual parallel phase covariance matrix. In the first algorithm, the use of a basis in which the latter matrix is diagonal reduces the number of $U\_{ij}$ functions to the number of mirror modes. In the second algorithm, this eigen decomposition is used to compute phase screens that follow the same statistics as the residual parallel phase covariance matrix, and thus suppress the need for these $U\_{ij}$ functions. Results. Our algorithms dramatically reduce the number of $U\_{ij}$ functions to be computed for the point-spread function reconstruction. Adaptive optics simulations show the good accuracy of both algorithms to reconstruct the point-spread function.

In the last few years, new Adaptive Optics [AO] techniques have emerged to answer new astronomical challenges: Ground-Layer AO [GLAO] and Multi-Conjugate AO [MCAO] to access a wider Field of View [FoV], Multi-Object AO [MOAO] for the simultaneous observation of several faint galaxies, eXtreme AO [XAO] for the detection of faint companions. In this paper, we focus our study to one of these applications : high red-shift galaxy observations using MOAO techniques in the framework of Extremely Large Telescopes [ELTs]. We present the high-level specifications of a dedicated instrument. We choose to describe the scientific requirements with the following criteria : 40% of Ensquared Energy [EE] in H band (1.65um) and in an aperture size from 25 to 150 mas. Considering these specifications we investigate different AO solutions thanks to Fourier based simulations. Sky Coverage [SC] is computed for Natural and Laser Guide Stars [NGS, LGS] systems. We show that specifications are met for NGS-based systems at the cost of an extremely low SC. For the LGS approach, the option of low order correction with a faint NGS is discussed. We demonstrate that, this last solution allows the scientific requirements to be met together with a quasi full SC.

In this paper, we present a first comparison of different Adaptive Optics (AO) concepts to reach a given scientific specification for 3D spectroscopy on Extremely Large Telescope (ELT). We consider that a range of 30%-50% of Ensquarred Energy (EE) in H band (1.65um) and in an aperture size from 25 to 100mas is representative of the scientific requirements. From these preliminary choices, different kinds of AO concepts are investigated : Ground Layer Adaptive Optics (GLAO), Multi-Object AO (MOAO) and Laser Guide Stars AO (LGS). Using Fourier based simulations we study the performance of these AO systems depending on the telescope diameter.

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

We report VLT/NACO imaging observations of the stars GSC 08047-00232 and HIP 6856, probable members of the large Tucana-Horologium association. During our previous ADONIS/SHARPII deep imaging survey, a substellar candidate companion was discovered around each star. Based on VLT/NACO astrometric measurements, we find that GSC 08047-00232 and the faint candidate companion near to it share the same proper motion with a significance of 3.1sigma. On the contrary, the candidate companion to HIP 6856 is probably a background object with a significance of 4.3sigma. We also detect a new fainter and closer candidate companion to HIP 6856, but which is likely a background object too with a significance of 4sigma. Recent VLT/NACO spectroscopic measurements of GSC 08047-00232 finally confirm the substellar nature of this young brown dwarf with a derived spectral type M9.5. GSC 08047-00232 B, with an estimated mass of 25MJup and an effective temperature of 2100K, is presently the third substellar companion identified among young, nearby associations.

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.

FALCON is an original concept for a next generation spectrograph at ESO VLT or at future ELTs. It is a spectrograph including multiple small integral field units (IFUs) which can be deployed within a large field of view such as that of VLT/GIRAFFE. In FALCON, each IFU features an adaptive optics correction using off-axis natural reference stars in order to combine, in the 0.8-1.8 \mu m wavelength range, spatial and spectral resolutions (0.1-0.15 arcsec and R=10000+/-5000). These conditions are ideally suited for distant galaxy studies, which should be done within fields of view larger than the galaxy clustering scales (4-9 Mpc), i.e. foV > 100 arcmin2. Instead of compensating the whole field, the adaptive correction will be performed locally on each IFU. This implies to use small miniaturized devices both for adaptive optics correction and wavefront sensing. Applications to high latitude fields imply to use atmospheric tomography because the stars required for wavefront sensing will be in most of the cases far outside the isoplanatic patch.

We report on Galactic Center L'-band observations made during NAOS/CONICA Science Verification. Colors of the inner 2" stars reveal an infrared excess of S2, the closest star to the black hole, that could sign the first thermal infrared detection of Sgr A*. A multi-wavelength maximum likelihood analysis has allowed us to eliminate all but two of the candidates for gravitational lensing proposed by Alexander (2001). Our observations of the thin and intersecting structures of the Northern Arm could trace several shocks heating the neighbooring dust rather than a stream of matter in orbit around the central mass as previously thought.

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.

Near-infrared observations of line emission from excited H2 and in the continuum are reported in the direction of the Orion molecular cloud OMC1, using the European Southern Observatory Very Large Telescope UT4, equipped with the NAOS adaptive optics system and the CONICA infrared array camera. Spatial resolution has been achieved at close to the diffraction limit of the telescope (0.08" - 0.12") and images show a wealth of morphological detail. Structure is not fractal but shows two preferred scale sizes of 2.4" (1100 AU) and 1.2" (540 AU), where the larger scale may be associated with star formation.

We report the detection of a possible subsolar mass companion to the triple young system TY CrA using the NACO instrument at the VLT UT4 during its commissioning. Assuming for TY CrA a distance similar to that of the close binary system HD 176386, the photometric spectral type of this fourth stellar component candidate is consistent with an ~M4 star. We discuss the dynamical stability of this possible quadruple system as well as the possible location of dusty particles inside or outside the system.

We report on the first results obtained using adaptive optics measurements of the Galactic Centre done with NAOS/CONICA.

Measurements of stellar velocities and variable X-ray emission near the centre of the Milky Way have provided the strongest evidence so far that the dark mass concentrations seen in many galactic nuclei are likely supermassive black holes, but have not yet excluded several alternative configurations. Here we report ten years of high resolution astrometric imaging that allow us to trace two thirds of the orbit of the star currently closest to the compact radio source and massive black hole candidate SgrA*. In particular, we have observed both peri- and apocentre passages. Our observations show that the star is on a bound, highly elliptical Keplerian orbit around SgrA*, with an orbital period of 15.2 years and a peri-centre distance of only 17 light hours. The orbital elements require an enclosed point mass of 3.7+-1.5x10^6 solar masses. The data exclude with high confidence that the central dark mass consists of a cluster of astrophysical objects or massive, degenerate fermions, and strongly constrain the central density structure.

A large fraction of the present-day stellar mass was formed between z=0.5 and z~3 and our understanding of the formation mechanisms at work at these epochs requires both high spatial and high spectral resolution: one shall simultaneously} obtain images of objects with typical sizes as small as 1-2kpc(~0''.1), while achieving 20-50 km/s (R >= 5000) spectral resolution. The obvious instrumental solution to adopt in order to tackle the science goal is therefore a combination of multi-object 3D spectrograph with multi-conjugate adaptive optics in large fields. A partial, but still competitive correction shall be prefered, over a much wider field of view. This can be done by estimating the turbulent volume from sets of natural guide stars, by optimizing the correction to several and discrete small areas of few arcsec2 selected in a large field (Nasmyth field of 25 arcmin) and by correcting up to the 6th, and eventually, up to the 60th Zernike modes. Simulations on real extragalactic fields, show that for most sources (>80%), the recovered resolution could reach 0".15-0".25 in the J and H bands. Detection of point-like objects is improved by factors from 3 to >10, when compared with an instrument without adaptive correction. The proposed instrument concept, FALCON, is equiped with deployable mini-integral field units (IFUs), achieving spectral resolutions between R=5000 and 20000. Its multiplex capability, combined with high spatial and spectral resolution characteristics, is a natural ground based complement to the next generation of space telescopes.

We present diffraction-limited near-IR images in J, H and K of the nucleus of NGC 1068, obtained with the Adaptive Optics system {Pueo} at CFHT. The achieved resolution (0.12") reveals several components, particularly prominent on the [J-K] image: a) an unresolved, conspicuous core (size < 9 pc); b) an elongated structure at P.A. ~102 deg, beginning to show up at radius ~ 15 pc; c) a S-shaped structure with radial extent ~ 20 pc, including a bar-like central elongation at P.A. ~ 15 deg and two short spiral arms. The K core is at the location of the putative central engine (radio source S1) : the core is likely the emission from the hot inner walls of the dust/molecular torus. The extremely red colors of the central 0.2", [J-K]=7.0, [H-K]=3.8, lead to an extinction Av > 25. The elongated structure at P.A.~ 102 deg may trace the presence of cooler dust within and around the torus. This interpretation is supported by two facts : a) the elongated structure is perpendicular to the local radio jet originating at S1; b) its direction follows exactly that of the disk of ionized gas recently found with the VLBA. The S-shaped feature suggests an extremely compact barred spiral structure, that would be the innermost of a series of nested spiral structures, as predicted by simulations. This is supported by the inner stellar distribution - deduced from the J image - which clearly follows an exponential disk with a 19 pc scale-length, precisely that expected from the rotation of a bar.