Search

The GPU-based High-order adaptive OpticS Testbench (GHOST) at the European Southern Observatory (ESO) is a new 2-stage extreme adaptive optics (XAO) testbench at ESO. The GHOST is designed to investigate and evaluate new control methods (machine learning, predictive control) for XAO which will be required for instruments such as the Planetary Camera and Spectrograph of ESOs Extremely Large Telescope. The first stage corrections are performed in simulation, with the residual wavefront error at each iteration saved. The residual wavefront errors from the first stage are then injected into the GHOST using a spatial light modulator. The second stage correction is made with a Boston Michromachines Corporation 492 actuator deformable mirror and a pyramid wavefront sensor. The flexibility of the bench also opens it up to other applications, one such application is investigating the flip-flop modulation method for the pyramid wavefront sensor., Comment: spie astronomical telescopes + instrumentation 2022

We propose a novel control approach that combines offline supervised learning to address the challenges posed by non-linear phase reconstruction using unmodulated pyramid wavefront sensors (P-WFS) and online reinforcement learning for predictive control. The control approach uses a high-order P-WFS to drive a tip-tilt stage and a high-dimensional mirror concurrently. Simulation results demonstrate that our method outperforms traditional control techniques, showing significant improvements in performance under challenging conditions such as faint stars and poor seeing, and exhibits robustness against variations in atmospheric conditions.

Differentiating between an exoplanet signal and residual speckle noise is a key challenge in high-contrast imaging. Speckles are due to a combination of fast, slow and static wavefront aberrations introduced by atmospheric turbulence and instrument optics. While wavefront control techniques developed over the last decade have shown promise in minimizing fast atmospheric residuals, slow and static aberrations such as non-common path aberrations (NCPAs) remain a key limiting factor for exoplanet detection. NCPA are not seen by the wavefront sensor (WFS) of the adaptive optics (AO) loop, hence the difficulty in correcting them. We propose to improve the identification and rejection of those aberrations. The algorithm DrWHO, performs frequent compensation of static and quasi-static aberrations to boost image contrast. By changing the WFS reference at every iteration of the algorithm, DrWHO changes the AO point of convergence to lead it towards a compensation of the static and slow aberrations. References are calculated using an iterative lucky-imaging approach, where each iteration updates the WFS reference, ultimately favoring high-quality focal plane images. We validate this concept through numerical simulations and on-sky testing on the SCExAO instrument at the 8.2-m Subaru telescope. Simulations show a rapid convergence towards the correction of 82% of the NCPAs. On-sky tests are performed over a 10-minute run in the visible (750 nm). We introduce a flux concentration (FC) metric to quantify the point spread function (PSF) quality and measure a 15.7% improvement. The DrWHO algorithm is a robust focal-plane wavefront sensing calibration method that has been successfully demonstrated on sky. It does not rely on a model nor requires wavefront sensor calibration or linearity. It is compatible with different wavefront control methods, and can be further optimized for speed and efficiency., Comment: Accepted in A&A. 12 pages, 15 figures. Abridged abstract

The central region of NGC 1068 is one of the closest and most studied Active Galactic Nucleus. It is known to be of type 2, meaning that its accretion disk is obscured by a large amount of dust and gas. The main properties of the obscuring structure are still to be determined. We model the inner edge of this structure, where the hot dust responsible for the near-infrared emission reaches its sublimation temperature. At first, we use simple geometrical models in image reconstructions to determine the main 2D geometrical properties of the source. In a second step, we try to reproduce the observables with K-band images produced by 3D radiative transfer simulations of a heated dusty disk. We explore various parameters to find an optimal solution and a model consistent with all the observables. The three methods indicate an elongated structure with 4 x 6 mas dimensions and its major axis along the NW/SE direction. All three of them suggest that the object looks like an elongated ring rather than an elongated thin disk, with its northeast edge less luminous than the South-West one. The best 3D model is a thick disk with an inner radius r=0.21 +/- 0.03 pc and a half-opening angle $\alpha_{1/2}=21\pm8\degree$ observed with an inclination $i=44_{-6}^{10}$\,$\degree$ and $PA=150_{-13}^{8}$\,$\degree$. A high density of dust $n=5_{-2.5}^{+5}\ M_{\odot}.pc^{-3}$ is required to explain the contrast between the two edges by self-absorption from the closer one. The overall structure is obscured by a large foreground obscuration $A_V \sim 75$. The hot dust is not responsible for the obscuration of the central engine. The geometry and the orientation of the structure are different from those of the previously observed maser and molecular disks. A single disk is unable to account for these differences, and we favor a description of the source where multiple rings are entangled around the central mass., Comment: Accepted for publication in Astronomy and Astrophysics

We propose the Correlation-Locking Optimization SchEme (CLOSE), a real-time adaptive filtering technique for adaptive optics (AO) systems controlled with integrators. CLOSE leverages the temporal autocorrelation of modal signals in the controller telemetry and drives the gains of the integral command law in a closed servo-loop. This supervisory loop is configured using only a few scalar parameters, and automatically controls the modal gains to closely match transfer functions achieving minimum variance control. This optimization is proven to work throughout the range of noise and seeing conditions relevant to the AO system. This technique has been designed while preparing the high-order AO systems for extremely large telescopes, in particular for tackling the optical gain (OG) phenomenon -- a sensitivity reduction induced by on-sky residuals -- which is a prominent issue with pyramid wavefront sensors (PWFS). CLOSE follows upon the linear modal compensation approach to OG, previously demonstrated to substantially improve AO correction with high order PWFS systems. Operating on modal gains through multiplicative increments, CLOSE naturally compensates for the recurring issue of unaccounted sensitivity factors throughout the AO loop. We present end-to-end simulations of the MICADO instrument single-conjugate AO to demonstrate the performances and capabilities of CLOSE. We demonstrate that a single configuration shall provide an efficient and versatile optimization of the modal integrator while accounting for OG compensation, and while providing significant robustness to transient effects impacting the PWFS sensitivity., Comment: 19 pages, 14 figures. To be published in Astronomy & Astrophysics

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., Comment: 12 pages - Accepted for publication

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., Comment: 20 pages. Revised version after referee's comments. Recommended for publication in A&A, after revision (this version)

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., Comment: 12 pages, 10 figures, regular paper accepted in Astronomy & Astrophysics (24 April 2018)

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.

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., Comment: Proceeding from conference SF2A 2016 - PNHE workshop

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., Comment: 13 pages, 7 figure. To appear in Adaptive Optics Systems V, eds. Marchetti E., Close L.M., V\'eran, J.-P., Proc. SPIE vol. 9909 id 54; 2016

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., Comment: Significantly expanded and updated version of previous ELT-MOS White Paper, so there is some textual overlap with arXiv:1303.0029

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., Comment: in press A&A (2014)

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., Comment: 17 pages, 7 figures, Accepted for publication in MNRAS

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., Comment: accepted for publication in ApJ

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., Comment: 19 pages, 14 figures, 2 tables, accepted for publication in Astrophysical Journal

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., Comment: 8 pages, 8 figures, 1 table, accepted for publication in Astronomy and Astrophysics

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., Comment: 12 pages & 13 figures. To be published in the proceedings of the SPIE conference Advances in Adaptive Optics - Astronomical Telescopes & Instrumentation, 24-31 May 2006, Orlando

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*., Comment: accepted in A&A

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., Comment: Submitted : 15 March 2003 ; accepted : 15 May 2003