Aims. We investigate the fueling and the feedback of nuclear activity in the nearby (D = 14 Mpc) Seyfert 2 barred galaxy NGC 1068 by studying the distribution and kinematics of molecular gas in the ...torus and its connections to the host galaxy disk. Methods. We used the Atacama Large Millimeter Array (ALMA ) to image the emission of a set of molecular gas tracers in the circumnuclear disk (CND) and the torus of the galaxy using the CO(2–1), CO(3–2), and HCO+(4–3) lines and their underlying continuum emission with high spatial resolutions (0.03″ − 0.09″ ≃ 2 − 6 pc). These transitions, which span a wide range of physical conditions of molecular gas (n(H2)⊂103 − 107 cm−3), are instrumental in revealing the density radial stratification and the complex kinematics of the gas in the torus and its surroundings. Results. The ALMA images resolve the CND as an asymmetric ringed disk of D ≃ 400 pc in size and ≃1.4 × 108 M⊙ in mass. The CND shows a marked deficit of molecular gas in its central ≃130 pc region. The inner edge of the ring is associated with the presence of edge-brightened arcs of NIR polarized emission, which are identified with the current working surface of the ionized wind of the active galactic nucleus (AGN). ALMA proves the existence of an elongated molecular disk/torus in NGC 1068 of Mtorusgas ≃ 3 × 105 M⊙ M torus gas ≃ 3 × 10 5 M ⊙ $ M_{\mathrm{torus}}^{\mathrm{gas}}\simeq3\times10^{5}\,M_{{\odot}} $ , which extends over a large range of spatial scales D ≃ 10 − 30 pc around the central engine. The new observations evidence the density radial stratification of the torus: the HCO+(4–3) torus, with a full size DHCO+(4 − 3) = 11 ± 0.6 pc, is a factor of between two and three smaller than its CO(2–1) and CO(3–2) counterparts, which have full sizes of DCO(3 − 2) = 26 ± 0.6 pc and DCO(2 − 1) = 28 ± 0.6 pc, respectively. This result brings into light the many faces of the molecular torus. The torus is connected to the CND through a network of molecular gas streamers detected inside the CND ring. The kinematics of molecular gas show strong departures from circular motions in the torus, the gas streamers, and the CND ring. These velocity field distortions are interconnected and are part of a 3D outflow that reflects the effects of AGN feedback on the kinematics of molecular gas across a wide range of spatial scales around the central engine. In particular, we estimate through modeling that a significant fraction of the gas inside the torus ( ≃ 0.4 − 0.6 × Mtorusgas ≃ 0.4 − 0.6 × M torus gas $ {\simeq}0.4{-}0.6 \times M_{\mathrm{torus}}^{\mathrm{gas}} $ ) and a comparable amount of mass along the gas streamers are outflowing. However, the bulk of the mass, momentum, and energy of the molecular outflow of NGC 1068 is contained at larger radii in the CND region, where the AGN wind and the radio jet are currently pushing the gas assembled at the Inner Lindblad Resonance (ILR) ring of the nuclear stellar bar. Conclusions. In our favored scenario a wide-angle AGN wind launched from the accretion disk of NGC1068 is currently impacting a sizable fraction of the gas inside the torus. However, a large gas reservoir (≃1.2 − 1.8 × 105 M⊙), which lies close to the equatorial plane of the torus, remains unaffected by the feedback of the AGN wind and can therefore continue fueling the AGN for at least ≃1 − 4 Myr. Nevertheless, AGN fueling currently seems thwarted on intermediate scales (15 pc ≤r ≤ 50 pc).
We present a novel formulation of closed-loop adaptive optics (AO) control as a multi-agent reinforcement learning (MARL) problem in which the controller is able to learn a non-linear policy and does ...not need a priori information on the dynamics of the atmosphere. We identify the different challenges of applying a reinforcement learning (RL) method to AO and, to solve them, propose the combination of model-free MARL for control with an autoencoder neural network to mitigate the effect of noise. Moreover, we extend current existing methods of error budget analysis to include a RL controller. The experimental results for an 8m telescope equipped with a 40x40 Shack-Hartmann system show a significant increase in performance over the integrator baseline and comparable performance to a model-based predictive approach, a linear quadratic Gaussian controller with perfect knowledge of atmospheric conditions. Finally, the error budget analysis provides evidence that the RL controller is partially compensating for bandwidth error and is helping to mitigate the propagation of aliasing.
The pyramid wavefront sensor (PWFS) is the currently preferred design for high-sensitivity adaptive optics (AO) systems for extremely large telescopes (ELTs). Yet, nonlinearities of the signal ...retrieved from the PWFS pose a significant problem for achieving the full correction potential using this sensor, a problem that will only worsen with the increasing dimension of telescopes. This paper investigates the so-called optical gain (OG) phenomenon, a sensitivity reduction and an overall modification of the sensor response induced by the residual wavefront itself, with considerable effects in standard observation conditions for ELT-sized AO systems. Through extensive numerical analysis, this work proposes a formalism to measure and minimize the first-order nonlinearity error caused by optical gain variation, which uses a modal compensation technique of the calibrated reconstructor; this enables a notable increase in performance in faint guide stars or important seeing scenarios, for example from 16 to 30% H-band Strehl ratio for a sixteenth magnitude star in r0 = 13 cm turbulence. Beyond the performance demonstrated by this compensation, a complete algorithm for realistic operation conditions is designed, which from dithering a few deformable mirror modes retrieves the optimal gains and updates the command matrix accordingly. The performance of this self-updating technique – which successfully allows automatic OG compensation regardless of the turbulent conditions, and its minimal interference with the scientific instrument are demonstrated through extensive end-to-end numerical simulations, all at the scale of an ELT instrument single-conjugate AO system.
Here, we show that the approximately 10-million-year-old β Pictoris system hosts a massive giant planet, β Pictoris b, located 8 to 15 astronomical units from the star. This result confirms that gas ...giant planets form rapidly within disks and validates the use of disk structures as fingerprints of embedded planets. Among the few planets already imaged, β Pictoris b is the closest to its parent star. Its short period could allow for recording of the full orbit within 17 years.
The broadening of atomic emission lines by high-velocity motion of gas near accreting supermassive black holes is an observational hallmark of quasars
. Observations of broad emission lines could ...potentially constrain the mechanism for transporting gas inwards through accretion disks or outwards through winds
. The size of regions for which broad emission lines are observed (broad-line regions) has been estimated by measuring the delay in light travel time between the variable brightness of the accretion disk continuum and the emission lines
-a method known as reverberation mapping. In some models the emission lines arise from a continuous outflow
, whereas in others they arise from orbiting gas clouds
. Directly imaging such regions has not hitherto been possible because of their small angular size (less than 10
arcseconds
). Here we report a spatial offset (with a spatial resolution of 10
arcseconds, or about 0.03 parsecs for a distance of 550 million parsecs) between the red and blue photo-centres of the broad Paschen-α line of the quasar 3C 273 perpendicular to the direction of its radio jet. This spatial offset corresponds to a gradient in the velocity of the gas and thus implies that the gas is orbiting the central supermassive black hole. The data are well fitted by a broad-line-region model of a thick disk of gravitationally bound material orbiting a black hole of 3 × 10
solar masses. We infer a disk radius of 150 light days; a radius of 100-400 light days was found previously using reverberation mapping
. The rotation axis of the disk aligns in inclination and position angle with the radio jet. Our results support the methods that are often used to estimate the masses of accreting supermassive black holes and to study their evolution over cosmic time.
Aims. 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. Methods. 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. Results. On a SPHERE-like system, the comparison between a PSF computed from the error breakdown with a PSF obtained from classical end-to-end simulation shows that the statistics convergence limits converge very well, with a sub-percent difference in terms of Strehl ratio and ensquared energy at 5λ/D$5\frac{\lambda}{D}$5λD separation. 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.
We present new near-infrared VLTI/GRAVITY interferometric spectra that spatially resolve the broad Br
γ
emission line in the nucleus of the active galaxy IRAS 09149−6206. We use these data to measure ...the size of the broad line region (BLR) and estimate the mass of the central black hole. Using an improved phase calibration method that reduces the differential phase uncertainty to 0.05° per baseline across the spectrum, we detect a differential phase signal that reaches a maximum of ∼0.5° between the line and continuum. This represents an offset of ∼120
μ
as (0.14 pc) between the BLR and the centroid of the hot dust distribution traced by the 2.3
μ
m continuum. The offset is well within the dust sublimation region, which matches the measured ∼0.6 mas (0.7 pc) diameter of the continuum. A clear velocity gradient, almost perpendicular to the offset, is traced by the reconstructed photocentres of the spectral channels of the Br
γ
line. We infer the radius of the BLR to be ∼65
μ
as (0.075 pc), which is consistent with the radius–luminosity relation of nearby active galactic nuclei derived based on the time lag of the H
β
line from reverberation mapping campaigns. Our dynamical modelling indicates the black hole mass is ∼1 × 10
8
M
⊙
, which is a little below, but consistent with, the standard
M
BH
–
σ
*
relation.
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 was designed while preparing the high-order AO systems for extremely large telescopes, in particular to tackle the optical gain (OG) phenomenon. This is a sensitivity reduction induced by on-sky residuals and is a prominent issue with pyramid wavefront sensors (PWFS). CLOSE follows upon the linear modal compensation approach to OG that was 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 provides 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.
Context: Since the discovery of its dusty disk in 1984, β 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. Aims: Our goal was to investigate the close environment of β Pic, searching for planetary companion(s). Methods: Deep adaptive-optics L'-band images of β Pic were recorded using the NaCo instrument at the Very Large Telescope. Results: A faint point-like signal is detected at a projected distance of ≃8 AU from the star, within the northeastern extension of the dust disk. Various tests were made to rule out possible instrumental or atmospheric artefacts at a good confidence level. The probability of a foreground or background contaminant is extremely low, based in addition on the analysis of previous deep HST images. Its L'=11.2 apparent magnitude would indicate a typical temperature of ~1500 K and a mass of ~8 M_Jup. If confirmed, it could explain the main morphological and dynamical peculiarities of the β 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 β Pic disk suggest a formation process by core accretion or disk instabilities rather than binary-like formation processes. Based on observations collected at the European Southern Observatory, Chile, ESO (runs 072.C-0624(B) and 60.A-9026(A)) and on observations made with the Space Telescope Imaging Spectrograph onboard the NASA/ESA Hubble Space Telescope.
We present near-infrared interferometric data on the Seyfert 2 galaxy NGC 1068, obtained with the GRAVITY instrument on the European Southern Observatory Very Large Telescope Interferometer. The ...extensive baseline coverage from 5 to 60
M
λ
allowed us to reconstruct a continuum image of the nucleus with an unrivaled 0.2 pc resolution in the
K
-band. We find a thin ring-like structure of emission with a radius
r
= 0.24 ± 0.03 pc, inclination
i
= 70 ± 5°, position angle PA = −50 ± 4°, and
h
/
r
< 0.14, which we associate with the dust sublimation region. The observed morphology is inconsistent with the expected signatures of a geometrically and optically thick torus. Instead, the infrared emission shows a striking resemblance to the 22 GHz maser disc, which suggests they share a common region of origin. The near-infrared spectral energy distribution indicates a bolometric luminosity of (0.4–4.7) × 10
45
erg s
−1
, behind a large
A
K
≈ 5.5 (
A
V
≈ 90) screen of extinction that also appears to contribute significantly to obscuring the broad line region.