Robotic observation pipeline for small bodies in the solar system based on open-source software and commercially available telescope hardware

The observation of small bodies in the Space Environment is an ongoing important task in astronomy. While nowadays new objects are mostly detected in larger sky surveys, several follow-up observations are usually needed for each object to improve the accuracy of orbit determination. In particular objects orbiting close to Earth, so called Near-Earth Objects (NEOs) are of special concern as a small but not negligible fraction of them can have a non-zero impact probability with Earth. Additionally, the observation of manmade space debris and tracking of satellites falls in the same class measurements. Telescopes for these follow-up observations are mainly in a aperture class between 1 m down to approximately 25 cm. These telescopes are often hosted by amateur observatories or dedicated companies like 6ROADS specialized on this type of observation. With upcoming new NEO search campaigns by very wide field of view telescopes, like the Vera C. Rubin Observatory, NASA’s NEO surveyor space mission and ESA’s Flyeye telescopes, the number of NEO discoveries will increase dramatically. This will require an increasing number of useful telescopes for follow-up observations at different geographical locations. While well-equipped amateur astronomers often host instruments which might be capable of creating useful measurements, both observation planning and scheduling, and also analysis are still a major challenge for many observers. In this work we present a fully robotic planning, scheduling and observation pipeline that extends the widely used open-source cross-platform software KStars/Ekos for Instrument Neutral Distributed Interface (INDI) devices. The method consists of algorithms which automatically select NEO candidates with priority according to ESA’s Near-Earth Object Coordination Centre (NEOCC). It then analyses detectable objects (based on limiting magnitudes, geographical position, and time) with preliminary ephemeris from the Minor Planet Center (MPC). Optimal observing slots during the night are calculated and scheduled. Immediately before the measurement the accurate position of the minor body is recalculated and finally the images are taken. Besides the detailed description of all components, we will show a complete robotic hard- and software solution based on our methods.TS-R acknowledges funding from the NEO-MAPP project (H2020-EU-2-1-6/870377). This work was (partially) funded by the Spanish MICIN/AEI/10.13039/501100011033 and by “ERDF A way of making Europe” by the “European Union” through grant RTI2018-095076-B-C21, and the Institute of Cosmos Sciences University of Barcelona (ICCUB, Unidad de Excelencia “María de Maeztu”) through grant CEX2019-000918-M

Asteroids seen by JWST-MIRI: Radiometric size, distance, and orbit constraints

Infrared measurements of asteroids are crucial for the determination of physical and thermal properties of individual objects, and for understanding the small-body populations in the solar system as a whole. However, standard radiometric methods can only be applied if the orbit of an object is known, hence its position at the time of the observation. With JWST-MIRI observations the situation will change and many unknown, often very small, solar system objects will be detected. Later orbit determinations are difficult due to the faintness of the objects and the lack of dedicated follow-up concepts. We present MIRI observations of the outer-belt asteroid (10920) 1998 BC1 and an unknown object, detected in all nine MIRI bands in close apparent proximity to (10920). We developed a new method called STM-ORBIT to interpret the multi-band measurements without knowing the object’s true location. The power of the new technique is that it determines the most-likely heliocentric and observer-centric distance and phase angle ranges, allowing us to make a radiometric size estimate. The application to the MIRI fluxes of (10920) was used to validate the method. It leads to a confirmation of the known radiometric size-albedo solution, and puts constraints on the asteroid’s location and orbit in agreement with its true orbit. To back up the validation of the method, we obtained additional ground-based light curve observations of (10920), combined with Gaia data, which indicate a very elongated object (a/b ≥ 1.5), with a spin-pole at (λ, β)ecl = (178°, +81°), with an estimated error of about 20°, and a rotation period of 4.861191 ± 0.000015 h. A thermophysical study of all available JWST-MIRI and WISE measurements leads to a size of 14.5–16.5 km (diameter of an equal-volume sphere), a geometric albedo pV between 0.05 and 0.10, and a thermal inertia in the range 9–35 (best value 15) J m−2 s−0.5 K−1. For the newly discovered MIRI object, the STM-ORBIT method revealed a size of 100–230 m. The new asteroid must be on a low-inclination orbit (0.7° < i < 2.0°) and it was located in the inner main-belt region during JWST observations. A beaming parameter η larger than 1.0 would push the size even below 100 meters, a main-belt regime that has escaped IR detections so far. This kind of MIRI observations can therefore contribute to formation and evolution studies via classical size-frequency studies, which are currently limited to objects larger than about one kilometer in size. We estimate that MIRI frames with pointings close to the ecliptic and short integration times of only a few seconds will always include a few asteroids; most of them will be unknown objects.TSR acknowledges funding from the NEO-MAPP project (H2020-EU-2-1-6/870377). This work was (partially) funded by the Spanish MICIN/AEI/10.13039/501100011033 and by “ERDF A way of making Europe” by the European Union through grant RTI2018-095076-B-C21, and the Institute of Cosmos Sciences University of Barcelona (ICCUB, Unidad de Excelencia ‘María de Maeztu’) through grant CEX2019-000918-M. PPB acknowledges funding through the Spanish Government retraining plan ‘María Zambrano 2021-2023’ at the University of Alicante (ZAMBRANO22-04)

Mining archival data from wide-field astronomical surveys in search of near-Earth objects

Context. Increasing our knowledge of the orbits and compositions of near-earth objects (NEOs) is important for a better understanding of the evolution of the Solar System and life. The detection of serendipitous NEO appearances among the millions of archived exposures from large astronomical imaging surveys can provide a contribution which is complementary to NEO surveys. Aims. Using the ASTROWISE information system, this work aims to assess the detectability rate, the achieved recovery rate and the quality of astrometry when data mining the European Southern Observatory (ESO) archive for the OmegaCAM wide-field imager at the VLT Survey Telescope (VST). Methods. We developed an automatic pipeline that searches for NEO appearances inside the ASTROWISE environment. Throughout the recovery process the pipeline uses several public web tools (SSOIS, NEODyS, JPL Horizons) to identify possible images that overlap with the positions of NEOs, and acquires information on the NEOs’ predicted position and other properties (e.g. magnitude, rate, and direction of motion) at the time of observations. Considering these properties, the pipeline narrows down the search to potentially detectable NEOs, searches for streak-like objects across the images, and finds a matching streak for the NEOs. Results. We recovered 196 appearances of NEOs from a set of 968 appearances predicted to be recoverable. It includes appearances for three NEOs that were on the impact risk list at that point. These appearances occurred well before their discovery. The subsequent risk assessment using the extracted astrometry removes these NEOs from the risk list. More generally, we estimate a detectability rate of ~0.05 per NEO at a signal-to-noise ratio higher than 3 for NEOs in the OmegaCAM archive. Our automatic recovery rates are 40% and 20% for NEOs on the risk list and the full list, respectively. The achieved astrometric and photometric accuracy is on average 0.12″ and 0.1 mag. Conclusions. These results show the high potential of the archival imaging data of the ground-based wide-field surveys as useful instruments for the search, (p)recovery, and characterization of NEOs. Highly automated approaches, as possible using ASTROWISE, make this undertaking feasible.This work was executed as part of ESA contract no. 4000134667/21/D/MRP (CARMEN) with their Planetary Defence Office. The Big Data Layer of the Target Field Lab project “Mining Big Data” was used. The Target Field Lab is supported by the Northern Netherlands Alliance (SNN) and is financially supported by the European Regional Development Fund. The data science software system ASTROWISE runs on powerful databases and computing clusters at the Donald Smits Center of the University of Groningen and is supported, among other parties, by NOVA (the Dutch Research School for Astronomy). TSR acknowledges funding from the NEO-MAPP project (H2020-EU-2-1-6/870377). This work was (partially) supported by the Spanish MICIN/AEI/10.13039/501100011033 and by “ERDF A way of making Europe” by the “European Union” through grant PID2021-122842OB-C21, and the Institute of Cosmos Sciences University of Barcelona (ICCUB, Unidad de Excelencia “María de Maeztu”) through grant CEX2019-000918-M

Unusual polarimetric properties for interstellar comet 2I/Borisov

So far, only two interstellar objects have been observed within our Solar System. While the first one, 1I/‘Oumuamua, had asteroidal characteristics, the second one, 2I/Borisov, showed clear evidence of cometary activity. We performed polarimetric observations of comet 2I/Borisov using the European Southern Observatory Very Large Telescope to derive the physical characteristics of its coma dust particles. Here we show that the polarization of 2I/Borisov is higher than what is typically measured for Solar System comets. This feature distinguishes 2I/Borisov from dynamically evolved objects such as Jupiter-family and all short- and long-period comets in our Solar System. The only object with similar polarimetric properties as 2I/Borisov is comet C/1995 O1 (Hale-Bopp), an object that is believed to have approached the Sun only once before its apparition in 1997. Unlike Hale-Bopp and many other comets, though, comet 2I/Borisov shows a polarimetrically homogeneous coma, suggesting that it is an even more pristine object.This work is based on observations made with ESO Telescopes at the La Silla Paranal Observatory thanks to Director Discretionary Time under programme ID 2104.C-5003. The work of T.S.-R. was carried out through grant APOSTD/2019/046 by Generalitat Valenciana (Spain). This work was supported by the MINECO (Spanish Ministry of Economy) through grant RTI2018-095076-B-C21 (MINECO/FEDER, UE)

A Deep and Wide Twilight Survey for Asteroids Interior to Earth and Venus

We are conducting a survey using twilight time on the Dark Energy Camera with the Blanco 4 m telescope in Chile to look for objects interior to Earth's and Venus' orbits. To date we have discovered two rare Atira/Apohele asteroids, 2021 LJ4 and 2021 PH27, which have orbits completely interior to Earth's orbit. We also discovered one new Apollo-type Near Earth Object (NEO) that crosses Earth's orbit, 2022 AP7. Two of the discoveries have diameters ≳1 km. 2022 AP7 is likely the largest Potentially Hazardous Asteroid (PHA) discovered in about eight years. To date we have covered 624 square degrees of sky near to and interior to the orbit of Venus. The average images go to 21.3 mag in the r band, with the best images near 22nd mag. Our new discovery 2021 PH27 has the smallest semimajor axis known for an asteroid, 0.4617 au, and the largest general relativistic effects (53 arcsec/century) known for any body in the solar system. The survey has detected ∼15% of all known Atira NEOs. We put strong constraints on any stable population of Venus co-orbital resonance objects existing, as well as the Atira and Vatira asteroid classes. These interior asteroid populations are important to complete the census of asteroids near Earth, including some of the most likely Earth impactors that cannot easily be discovered in other surveys. Comparing the actual population of asteroids found interior to Earth and Venus with those predicted to exist by extrapolating from the known population exterior to Earth is important to better understand the origin, composition, and structure of the NEO population.D.T. was supported by NASA grant 80NSSC21K0807. P.P. was supported by NASA ISFM EIMM award, the NASA Cooperative Agreement 80GSFC21M0002 and NASA solar system Workings award 80NSSC21K0153. T.S.R. acknowledges funding from the NEO-MAPP project (H2020-EU-2-1-6/870377). C.M. was partially supported by the international Gemini Observatory, a program of NSF's NOIRLab, which is managed by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation, on behalf of the Gemini partnership of Argentina, Brazil, Canada, Chile, the Republic of Korea, and the United States of America. This work was (partially) funded by the Spanish MICIN/AEI/10.13039/501100011033 and by "ERDF A way of making Europe" by the "European Union" through grant RTI2018-095076-B-C21, and the Institute of Cosmos Sciences University of Barcelona (ICCUB, Unidad de Excelencia María de Maeztu') through grant CEX2019-000918-M

The equilibrium shape of (65) Cybele: primordial or relic of a large impact?

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 × 1019 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.This work has been supported by the Czech Science Foundation through grants 20-08218S (J. Hanuš) and 21-11058S (M. Brož), as well as by the National Science Foundation under Grant No. 1743015 (F. Marchis). T. Santana-Ros acknowledges funding from the NEO-MAPP project (H2020-EU-2-1-6/870377). In addition, this work was partially funded by the Spanish MICIN/AEI/10.13039/501100011033 and by “ERDF A way of making Europe” by the “European Union” through grant RTI2018-095076-B-C21, and the Institute of Cosmos Sciences University of Barcelona (ICCUB, Unidad de Excelencia ‘María de Maeztu’) through grant CEX2019-000918-M. This research has made use of the Asteroid Families Portal maintained at the Department of Astronomy, University of Belgrade. TRAPPIST is a project funded by the Belgian Fonds (National) de la Recherche Scientifique (F.R.S.-FNRS) under grant PDR T.0120.21. TRAPPIST-North is a project funded by the University of Liège, in collaboration with the Cadi Ayyad University of Marrakech (Morocco). E. Jehin is F.R.S.-FNRS Senior Research Associate

A large topographic feature on the surface of the trans-Neptunian object (307261) 2002 MS4 measured from stellar occultations

Rommel, F. L., et al.[Context] The physical characterization of trans-Neptunian objects is essential for improving our understanding of the formation and evolution of our Solar System. Stellar occultation is a ground-based technique that can be successfully used to determine some of the TNOs’ fundamental physical properties with high precision, such as size and shape.[Aims] This work is aimed at constraining the size, shape, and geometric albedo of the dwarf planet candidate (307261) 2002 MS4 through the analysis of nine stellar occultation events. Using multichord detection, we also study the object’s topography by analyzing the obtained limb and residuals between the observed chords and the best-fit ellipse.[Metods] We predicted and organized the observational campaigns of nine stellar occultations by 2002 MS4 between 2019 and 2022, resulting in two single-chord events, four double-chord detections, and three events with between 3 and 61 positive chords. We derived the occultation light curves using differential aperture photometry, from which the star ingress and egress instants were calculated. Using 13 selected chords from the 8 August 2020 event, we determined the global elliptical limb of 2002 MS4. The best-fit ellipse, combined with the object’s rotational information from the literature, sets constraints on the object’s size, shape, and albedo. Additionally, we developed a new method to characterize the topography features on the object’s limb.[Results] The global limb has a semi-major axis of a′ = 412 ± 10 km, a semi-minor axis of b′= 385 ± 17 km, and the position angle of the minor axis is 121° ± 16°. From this instantaneous limb, we obtained 2002 MS4’s geometric albedo of pV = 0.1 ± 0.025, using HV = 3.63 ± 0.05 mag and a projected area-equivalent diameter of 796 ± 24 km. Significant deviations from the fitted ellipse in the northernmost limb were detected from multiple sites, highlighting three distinct topographic features: one 11 km depth depression, followed by a 25−5+4 km height elevation next to a crater-like depression, with an extension of 322 ± 39 km and 45.1 ± 1.5 km deep.[Conclusions] Our results indicate the presence of an object that is ≈138 km smaller in diameter than that derived from thermal data, possibly indicating the presence of a thus-far unknown satellite. However, within the error bars, the geometric albedo in the V-band is in agreement with the results published in the literature, even with the radiometric-derived albedo. This stellar occultation has allowed for the first multichord measurement of a large topography in a TNO.This study was financed in part by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior – Brazil (CAPES) --Finance Code 001, the National Institute of Science and Technology of the e-Universe project (INCT do e-Universo, CNPq grant 465376/2014-2), the Spanish MICIN/AEI/10.13039/501100011033, the Institute of Cosmos Sciences University of Barcelona (ICCUB, Unidad de Excelencia ’María de Maeztu’) through grant CEX2019-000918-M, the “ERDF A way of making Europe” by the “European Union” through grant PID2021-122842OB-C21, and within the “Lucky Star” umbrella that agglomerates the efforts of the Paris, Granada, and Rio teams, which the European Research Council funds under the European Community’s H2020 (ERC Grant Agreement No. 669416). The following authors acknowledge the respective CNPq grants: F. L. R. 103096/2023-0; F.B.-R. 314772/2020-0; R.V.-M. 307368/2021-1; J.I.B.C. 308150/2016-3 and 305917/2019-6; M.A. 427700/2018-3, 310683/2017-3, 473002/2013-2; G.M. 128580/2020-8; B.E.M. 150612/2020-6; and O.C.W. 305210/2018-1. The following authors acknowledge the respective grants: B.E.M. thanks the CAPES/Cofecub-394/2016-05; G.M. thanks the CAPES grant 88887.705245/2022-00; G.B-R. acknowledges CAPES – FAPERJ/PAPDRJ grant E26/203.173/2016 and the scholarship granted in the scope of the Program CAPES-PrInt, process number 88887.310463/2018-00, Mobility number 88887.571156/2020-00; M.A. acknowledges FAPERJ grant E-26/111.488/2013; A.R.G.Jr acknowledges FAPESP grant 2018/11239-8; O.C.W. and R.S. acknowledge FAPESP grant 2016/24561-0; K.B. acknowledges the scholarship funded by F.R.S.-FNRS grant T.0109.20 and by the Francqui Foundation; D.N. acknowledges the support from the French Centre National d’Etudes Spatiales (CNES); D. S. thank to Fulbright Visiting Scholar (2022-2023) at the University of California, Berkeley; A.P. and R.S. thank to the National Research, Development and Innovation Office (NKFIH, Hungary) grants K-138962 and K-125015. Partial funding for the computational infrastructure and database servers is received from the grant KEP-7/2018 of the Hungarian Academy of Sciences; R.D., J.L.O., P.S.-S., N.M., R.H., A.S.L., and J.M. T.-R. acknowledge the MCIN/AEI/10.13039/501100011033 under the grant respective grants: CEX2021-001131-S, PID2019-109467GB-I00, and PID2021-128062NB-I00; T.S.R. acknowledges funding from the NEO-MAPP project (H2020-EU-2-1-6/870377); K.H. was supported by the project R.V.O.: 67985815; A.K. thanks to the IRAP, Midi-Pyrenees Observatory, CNRS, University of Toulouse, France; J.M.O. acknowledges the Portuguese Foundation for Science and Technology (FCT) and the European Social Fund (ESF) through the Ph.D. grant SFRH/BD/131700/2017; J.d.W. and MIT acknowledge the Heising-Simons Foundation, Dr. and Mrs. Colin Masson, and Dr. Peter A. Gilman for Artemis, the first telescope of the SPECULOOS network situated in Tenerife, Spain. The ULiege’s contribution to SPECULOOS has received funding from the ERC under the European Union’s Seventh Framework Programme (FP/2007–2013) (grant number 336480/SPECULOOS); J.L. acknowledges the ACIISII, Consejería de Economia, Conocimiento y Empleo del Gobierno de Canarias, and the European Regional Development Fund (ERDF) under the grant ProID2021010134, also the Agencia Estatal de Investigación del Ministerio de Ciencia e Innovacion (AEI-MCINN) under the grant PID2020-120464GB-100; D.T., R.K., M.H., and T.P. were supported by the Slovak Grant Agency for Science grants number VEGA 2/0059/22, and VEGA 2/0031/22; M. P. was supported by a grant from the Romanian National Authority for Scientific Research — UEFISCDI, project number PN-III-P1-1.1-TE-2019-1504; P.B., M.M., and M.D.G. thank the support of the Italian Amateur Astronomers Union (UAI); C.F. acknowledges the support from ANID BASAL project FB210003, and grant MAS IC120009. This work has made use of data from the European Space Agency (ESA) mission Gaia (https://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www.cosmos.esa.int/web/gaia/dpac/consortium). Funding for the DPAC has been provided by national institutions, in particular, the institutions participating in the Gaia Multilateral Agreement. The Joan Oró Telescope (TJO) of the Montsec Observatory (OdM) is owned by the Catalan Government and operated by the Institute for Space Studies of Catalonia (IEEC). TCH telescope is financed by the Island Council of Ibiza. İST60 and IST40 are the observational facilities of the Istanbul University Observatory, funded by the Scientific Research Projects Coordination Unit of Istanbul University with project numbers BAP-3685 and FBG-2017-23943 and the Presidency of Strategy and Budget of the Republic of Turkey with the project 2016K12137. TRAPPIST is a project funded by the Belgian FNRS grant PDR T.0120.21 and the ARC grant for Concerted Research Actions, financed by the Wallonia-Brussels Federation. E.J. is an FNRS Senior Research Associate. TRAPPIST-North is funded by the University of Liège and performed in collaboration with the Cadi Ayyad University of Marrakesh. This work made use of observations obtained at the 1.6 m telescope on the Pico dos Dias Observatory of the National Laboratory of Astrophysics (LNA/Brazil), at the Copernicus and Schmidt telescopes (Asiago, Italy) of the INAF-Astronomical Observatory of Padova, and at the Southern Astrophysical Research (SOAR) telescope, which is a joint project of the Ministério da Ciência, Tecnologia, e Inovação (MCTI) da República Federativa do Brasil, the U.S. National Optical Astronomy Observatory (NOAO), the University of North Carolina at Chapel Hill (UNC), and Michigan State University (MSU).With funding from the Spanish government through the "Severo Ochoa Centre of Excellence" accreditation (CEX2019-000918-M).With funding from the Spanish government through the "Severo Ochoa Centre of Excellence" accreditation (CEX2021-001131-S).Peer reviewe

Asteroid (16) Psyche’s primordial shape: A possible Jacobi ellipsoid

Context. Asteroid (16) Psyche is the largest M-type asteroid in the main belt and the target of the NASA Psyche mission. It is also the only asteroid of this size (D >  200 km) known to be metal rich. Although various hypotheses have been proposed to explain the rather unique physical properties of this asteroid, a perfect understanding of its formation and bulk composition is still missing. Aims. We aim to refine the shape and bulk density of (16) Psyche and to perform a thorough analysis of its shape to better constrain possible formation scenarios and the structure of its interior. Methods. We obtained disk-resolved VLT/SPHERE/ZIMPOL images acquired within our ESO large program (ID 199.C-0074), which complement similar data obtained in 2018. Both data sets offer a complete coverage of Psyche’s surface. These images were used to reconstruct the three-dimensional (3D) shape of Psyche with two independent shape modeling algorithms (MPCD and ADAM). A shape analysis was subsequently performed, including a comparison with equilibrium figures and the identification of mass deficit regions. Results. Our 3D shape along with existing mass estimates imply a density of 4.20  ±  0.60 g cm−3, which is so far the highest for a solar system object following the four telluric planets. Furthermore, the shape of Psyche presents small deviations from an ellipsoid, that is, prominently three large depressions along its equator. The flatness and density of Psyche are compatible with a formation at hydrostatic equilibrium as a Jacobi ellipsoid with a shorter rotation period of ∼3h. Later impacts may have slowed down Psyche’s rotation, which is currently ∼4.2 h, while also creating the imaged depressions. Conclusions. Our results open the possibility that Psyche acquired its primordial shape either after a giant impact while its interior was already frozen or while its interior was still molten owing to the decay of the short-lived radionuclide 26Al.Based on observations collected at the European Organisation for Astronomical Research in the Southern Hemisphere under ESO programme 199.C-0074 (principal investigator: P. Vernazza). P. Vernazza, A. Drouard, M. Ferrais and B. Carry were supported by CNRS/INSU/PNP. J.H. and J.D. were supported by grant 18-09470S of the Czech Science Foundation and by the Charles University Research Programme no. UNCE/SCI/023. E.J. is F.R.S.-FNRS Senior Research Associate. The work of TSR was carried out through grant APOSTD/2019/046 by Generalitat Valenciana (Spain). This work was supported by the MINECO (Spanish Ministry of Economy) through grant RTI2018-095076-B-C21 (MINECO/FEDER, UE)

Binary asteroid (31) Euphrosyne: ice-rich and nearly spherical

Aims. Asteroid (31) Euphrosyne is one of the biggest objects in the asteroid main belt and it is also the largest member of its namesake family. The Euphrosyne family occupies a highly inclined region in the outer main belt and contains a remarkably large number of members, which is interpreted as an outcome of a disruptive cratering event. Methods. The goals of this adaptive-optics imaging study are threefold: to characterize the shape of Euphrosyne, to constrain its density, and to search for the large craters that may be associated with the family formation event. Results. We obtained disk-resolved images of Euphrosyne using SPHERE/ZIMPOL at the ESO 8.2 m VLT as part of our large program (ID: 199.C-0074, PI: Vernazza). We reconstructed its 3D shape via the ADAM shape modeling algorithm based on the SPHERE images and the available light curves of this asteroid. We analyzed the dynamics of the satellite with the Genoid meta-heuristic algorithm. Finally, we studied the shape of Euphrosyne using hydrostatic equilibrium models. Conclusions. Our SPHERE observations show that Euphrosyne has a nearly spherical shape with the sphericity index of 0.9888 and its surface lacks large impact craters. Euphrosyne’s diameter is 268 ± 6 km, making it one of the top ten largest main belt asteroids. We detected a satellite of Euphrosyne – S/2019 (31) 1 – that is about 4 km across, on a circular orbit. The mass determined from the orbit of the satellite together with the volume computed from the shape model imply a density of 1665 ± 242 kg m−3, suggesting that Euphrosyne probably contains a large fraction of water ice in its interior. We find that the spherical shape of Euphrosyne is a result of the reaccumulation process following the impact, as in the case of (10) Hygiea. However, our shape analysis reveals that, contrary to Hygiea, the axis ratios of Euphrosyne significantly differ from those suggested by fluid hydrostatic equilibrium following reaccumulation.This work has been supported by the Czech Science Foundation through grant 18-09470S (J. Hanuš, O. Chrenko, P. Ševeček) and by the Charles University Research program No. UNCE/SCI/023. M.B. was supported by the Czech Science Foundation grant 18-04514J. Computational resources were supplied by the Ministry of Education, Youth and Sports of the Czech Republic under the projects CESNET (LM2015042) and IT4Innovations National Supercomputing Centre (LM2015070). P. Vernazza, A. Drouard, M. Ferrais and B. Carry were supported by CNRS/INSU/PNP. M.M. was supported by the National Aeronautics and Space Administration under grant No. 80NSSC18K0849 issued through the Planetary Astronomy Program. The work of TSR was carried out through grant APOSTD/2019/046 by Generalitat Valenciana (Spain). This work was supported by the MINECO (Spanish Ministry of Economy) through grant RTI2018-095076-B-C21 (MINECO/FEDER, UE). The research leading to these results has received funding from the ARC grant for Concerted Research Actions, financed by the Wallonia-Brussels Federation. TRAPPIST is a project funded by the Belgian Fonds (National) de la Recherche Scientifique (F.R.S.-FNRS) under grant FRFC 2.5.594.09.F. TRAPPIST-North is a project funded by the Université de Liège, and performed in collaboration with Cadi Ayyad University of Marrakesh. E. Jehin is a FNRS Senior Research Associate