Context.
Adaptive optics (AO) systems greatly increase the resolution of large telescopes, but produce complex point spread function (PSF) shapes, varying in time and across the field of view. The ...PSF must be accurately known since it provides crucial information about optical systems for design, characterization, diagnostics, and image post-processing.
Aims.
We develop here a model of the AO long-exposure PSF, adapted to various seeing conditions and any AO system. This model is made to match accurately both the core of the PSF and its turbulent halo.
Methods.
The PSF model we develop is based on a parsimonious parameterization of the phase power spectral density, with only five parameters to describe circularly symmetric PSFs and seven parameters for asymmetrical ones. Moreover, one of the parameters is the Fried parameter
r
0
of the turbulence’s strength. This physical parameter is an asset in the PSF model since it can be correlated with external measurements of the
r
0
, such as phase slopes from the AO real time computer (RTC) or site seeing monitoring.
Results.
We fit our model against end-to-end simulated PSFs using the OOMAO tool, and against on-sky PSFs from the SPHERE/ZIMPOL imager and the MUSE integral field spectrometer working in AO narrow-field mode. Our model matches the shape of the AO PSF both in the core and the halo, with a relative error smaller than 1% for simulated and experimental data. We also show that we retrieve the
r
0
parameter with sub-centimeter precision on simulated data. For ZIMPOL data, we show a correlation of 97% between our
r
0
estimation and the RTC estimation. Finally, MUSE allows us to test the spectral dependency of the fitted
r
0
parameter. It follows the theoretical
λ
6/5
evolution with a standard deviation of 0.3 cm. Evolution of other PSF parameters, such as residual phase variance or aliasing, is also discussed.
ABSTRACT
High angular resolution disc-resolved images of (7) Iris collected by VLT/SPHERE instrument are allowed for the detailed shape modelling of this large asteroid revealing its surface ...features. If (7) Iris did not suffer any events catastrophic enough to disrupt the body (which is very likely) by studying its topography, we might get insights into the early Solar system’s collisional history. When it comes to internal structure and composition, thoroughly assessing the volume and density uncertainties is necessary. In this work, we propose a method of uncertainty calculation of asteroid shape models based on light curve and adaptive optics (AO) images. We apply this method on four models of (7) Iris produced from independent Shaping Asteroids using Genetic Evolution and All-Data Asteroid Modelling inversion techniques and multiresolution photoclinometry by deformation. Obtained diameter uncertainties stem from both the observations from which the models were scaled and the models themselves. We show that despite the availability of high-resolution AO images, the volume and density of (7) Iris have substantial error bars that were underestimated in the previous studies.
Context
. Access to knowledge of the point spread function (PSF) of adaptive optics(AO)-assisted observations is still a major limitation when processing AO data. This limitation is particularly ...important when image analysis requires the use of deconvolution methods. As the PSF is a complex and time-varying function, reference PSFs acquired on calibration stars before or after the scientific observation can be too different from the actual PSF of the observation to be used for deconvolution, and lead to artefacts in the final image.
Aims
. We improved the existing PSF-estimation method based on the so-called marginal approach by enhancing the object prior in order to make it more robust and suitable for observations of resolved extended objects.
Methods
. Our process is based on a two-step blind deconvolution approach from the literature. The first step consists of PSF estimation from the science image. For this, we made use of an analytical PSF model, whose parameters are estimated based on a marginal algorithm. This PSF was then used for deconvolution. In this study, we first investigated the requirements in terms of PSF parameter knowledge to obtain an accurate and yet resilient deconvolution process using simulations. We show that current marginal algorithms do not provide the required level of accuracy, especially in the presence of small objects. Therefore, we modified the marginal algorithm by providing a new model for object description, leading to an improved estimation of the required PSF parameters.
Results
. Our method fulfills the deconvolution requirement with realistic system configurations and different classes of Solar System objects in simulations. Finally, we validate our method by performing blind deconvolution with SPHERE/ZIMPOL observations of the Kleopatra asteroid.
ABSTRACT
One of the major limitations of using adaptive optics (AO) to correct image post-processing is the lack of knowledge about the system’s point spread function (PSF). The PSF is not always ...available as direct imaging on isolated point-like objects, such as stars. The use of AO telemetry to predict the PSF also suffers from serious limitations and requires complex and yet not fully operational algorithms. A very attractive solution is to estimate the PSF directly from the scientific images themselves, using blind or myopic post-processing approaches. We demonstrate that such approaches suffer from severe limitations when a joint restitution of object and PSF parameters is performed. As an alternative, here we propose a marginalized PSF identification that overcomes this limitation. In this case, the PSF is used for image post-processing. Here we focus on deconvolution, a post-processing technique to restore the object, given the image and the PSF. We show that the PSF estimated by marginalization provides good-quality deconvolution. The full process of marginalized PSF estimation and deconvolution constitutes a successful blind deconvolution technique. It is tested on simulated data to measure its performance. It is also tested on experimental AO images of the asteroid 4-Vesta taken by the Spectro-Polarimetric High-contrast Exoplanet Research (SPHERE)/Zurich Imaging Polarimeter (Zimpol) on the Very Large Telescope to demonstrate application to on-sky data.
Context.
Cybele asteroids constitute an appealing reservoir of primitive material genetically linked to the outer Solar System, and the physical properties (size and shape) of the largest members can ...be readily accessed by large (8m class) telescopes.
Aims.
We took advantage of the bright apparition of the most iconic member of the Cybele population, (65) Cybele, in July and August 2021 to acquire high-angular-resolution images and optical light curves of the asteroid with which we aim to analyse its shape and bulk properties.
Methods.
Eight series of images were acquired with VLT/SPHERE+ZIMPOL, seven of which were combined with optical light curves to reconstruct the shape of the asteroid using the
ADAM
,
MPCD
, and
SAGE
algorithms. The origin of the shape was investigated by means of N-body simulations.
Results.
Cybele has a volume-equivalent diameter of 263±3 km and a bulk density of 1.55 ± 0.19 g cm
−3
. Notably, its shape and rotation state are closely compatible with those of a Maclaurin equilibrium figure. The lack of a collisional family associated with Cybele and the higher bulk density of that body with respect to other large P-type asteroids suggest that it never experienced any large disruptive impact followed by rapid re-accumulation. This would imply that its present-day shape represents the original one. However, numerical integration of the long-term dynamical evolution of a hypothetical family of Cybele shows that it is dispersed by gravitational perturbations and chaotic diffusion over gigayears of evolution.
Conclusions.
The very close match between Cybele and an equilibrium figure opens up the possibility that
D
≥ 260 km (
M
≥ 1.5 × 10
19
kg) small bodies from the outer Solar System all formed at equilibrium. However, we cannot currently rule out an old impact as the origin of the equilibrium shape of Cybele. Cybele itself is found to be dynamically unstable, implying that it was ‘recently’ (<1 Gyr ago) placed on its current orbit either through slow diffusion from a relatively stable orbit in the Cybele region or, less likely, from an unstable, Jupiter-family-comet orbit in the planet-crossing region.
Context. Over the past decades, several interplanetary missions have studied small bodies in situ, leading to major advances in our understanding of their geological and geophysical properties. These ...missions, however, have had a limited number of targets. Among them, the NASA Dawn mission has characterised in detail the topography and albedo variegation across the surface of asteroid (4) Vesta down to a spatial resolution of ~20 m pixel−1 scale. Aims. Here our aim was to determine how much topographic and albedo information can be retrieved from the ground with VLT/SPHERE in the case of Vesta, having a former space mission (Dawn) providing us with the ground truth that can be used as a benchmark. Methods. We observed Vesta with VLT/SPHERE/ZIMPOL as part of our ESO large programme (ID 199.C-0074) at six different epochs, and deconvolved the collected images with a parametric point spread function (PSF). We then compared our images with synthetic views of Vesta generated from the 3D shape model of the Dawn mission, on which we projected Vesta’s albedo information. Results. We show that the deconvolution of the VLT/SPHERE images with a parametric PSF allows the retrieval of the main topographic and albedo features present across the surface of Vesta down to a spatial resolution of ~20–30 km. Contour extraction shows an accuracy of ~1 pixel (3.6 mas). The present study provides the very first quantitative estimate of the accuracy of ground-based adaptive-optics imaging observations of asteroid surfaces. Conclusions. In the case of Vesta, the upcoming generation of 30–40 m telescopes (ELT, TMT, GMT) should in principle be able to resolve all of the main features present across its surface, including the troughs and the north–south crater dichotomy, provided that they operate at the diffraction limit.
Context. Asteroid (16) Psyche is the target of the NASA Psyche mission. It is considered one of the few main-belt bodies that could be an exposed proto-planetary metallic core and that would thus be ...related to iron meteorites. Such an association is however challenged by both its near- and mid-infrared spectral properties and the reported estimates of its density. Aims. Here, we aim to refine the density of (16) Psyche to set further constraints on its bulk composition and determine its potential meteoritic analog. Methods. We observed (16) Psyche with ESO VLT/SPHERE/ZIMPOL as part of our large program (ID 199.C-0074). We used the high angular resolution of these observations to refine Psyche’s three-dimensional (3D) shape model and subsequently its density when combined with the most recent mass estimates. In addition, we searched for potential companions around the asteroid. Results. We derived a bulk density of 3.99 ± 0.26 g cm−3 for Psyche. While such density is incompatible at the 3-sigma level with any iron meteorites (∼7.8 g cm−3), it appears fully consistent with that of stony-iron meteorites such as mesosiderites (density ∼4.25 g cm−3). In addition, we found no satellite in our images and set an upper limit on the diameter of any non-detected satellite of 1460 ± 200 m at 150 km from Psyche (0.2% × RHill, the Hill radius) and 800 ± 200 m at 2000 km (3% × RHill). Conclusions. Considering that the visible and near-infrared spectral properties of mesosiderites are similar to those of Psyche, there is merit to a long-published initial hypothesis that Psyche could be a plausible candidate parent body for mesosiderites.
Context.
The recent estimates of the 3D shape of the M/Xe-type triple asteroid system (216) Kleopatra indicated a density of ~5 g cm
−3
, which is by far the highest for a small Solar System body. ...Such a high density implies a high metal content as well as a low porosity which is not easy to reconcile with its peculiar “dumbbell” shape.
Aims.
Given the unprecedented angular resolution of the VLT/SPHERE/ZIMPOL camera, here, we aim to constrain the mass (via the characterization of the orbits of the moons) and the shape of (216) Kleopatra with high accuracy, hence its density.
Methods.
We combined our new VLT/SPHERE observations of (216) Kleopatra recorded during two apparitions in 2017 and 2018 with archival data from the W. M. Keck Observatory, as well as lightcurve, occultation, and delay-Doppler images, to derive a model of its 3D shape using two different algorithms (ADAM, MPCD). Furthermore, an
N
-body dynamical model allowed us to retrieve the orbital elements of the two moons as explained in the accompanying paper.
Results.
The shape of (216) Kleopatra is very close to an equilibrium dumbbell figure with two lobes and a thick neck. Its volume equivalent diameter (118.75 ± 1.40) km and mass (2.97 ± 0.32) × 10
18
kg (i.e., 56% lower than previously reported) imply a bulk density of (3.38 ± 0.50) g cm
−3
. Such a low density for a supposedly metal-rich body indicates a substantial porosity within the primary. This porous structure along with its near equilibrium shape is compatible with a formation scenario including a giant impact followed by reaccumulation. (216) Kleopatra’s current rotation period and dumbbell shape imply that it is in a critically rotating state. The low effective gravity along the equator of the body, together with the equatorial orbits of the moons and possibly rubble-pile structure, opens the possibility that the moons formed via mass shedding.
Conclusions.
(216) Kleopatra is a puzzling multiple system due to the unique characteristics of the primary. This system certainly deserves particular attention in the future, with the Extremely Large Telescopes and possibly a dedicated space mission, to decipher its entire formation history.
Context. Asteroid (7) Iris is an ideal target for disk-resolved imaging owing to its brightness (V ~ 7–8) and large angular size of 0.33′′ during its apparitions. Iris is believed to belong to the ...category of large unfragmented asteroids that avoided internal differentiation, implying that its current shape and topography may record the first few 100 Myr of the solar system’s collisional evolution. Aims. We recovered information about the shape and surface topography of Iris from disk-resolved VLT/SPHERE/ZIMPOL images acquired in the frame of our ESO large program. Methods. We used the All-Data Asteroid Modeling (ADAM) shape reconstruction algorithm to model the 3D shape of Iris, using optical disk-integrated data and disk-resolved images from SPHERE and earlier AO systems as inputs. We analyzed the SPHERE images and our model to infer the asteroid’s global shape and the morphology of its main craters. Results. We present the 3D shape, volume-equivalent diameter Deq = 214 ± 5 km, and bulk density ρ = 2.7 ± 0.3 g cm−3 of Iris. Its shape appears to be consistent with that of an oblate spheroid with a large equatorial excavation. We identified eight putative surface features 20–40 km in diameter detected at several epochs, which we interpret as impact craters, and several additional crater candidates. Craters on Iris have depth-to-diameter ratios that are similar to those of analogous 10 km craters on Vesta. Conclusions. The bulk density of Iris is consistent with that of its meteoritic analog based on spectroscopic observations, namely LL ordinary chondrites. Considering the absence of a collisional family related to Iris and the number of large craters on its surface, we suggest that its equatorial depression may be the remnant of an ancient (at least 3 Gyr) impact. Iris’s shape further opens the possibility that large planetesimals formed as almost perfect oblate spheroids. Finally, we attribute the difference in crater morphology between Iris and Vesta to their different surface gravities, and the absence of a substantial impact-induced regolith on Iris.