Marco FENUCCI

Dr. Marco FENUCCI

orcid

Research Focus & Metrics

Academic Career Length

9Years (2018-Present)

(Calculated once Scilit publication data is synced: current year − year of first publication + 1)

Representative Publications

View All Publications
20 August 2025
Verified
Open access
DSDarryl Z. Seligman
Darryl Z. Seligman
MMMarco Micheli
Marco Micheli
DFDavide Farnocchia
Davide Farnocchia
LDLarry Denneau
Larry Denneau
JNJohn W. Noonan
John W. Noonan
QZQicheng Zhang
Qicheng Zhang
We report initial observations aimed at the characterization of a third interstellar object. This object, 3I/ATLAS or C/2025 N1 (ATLAS), was discovered on 2025 July 1 UT and has an orbital eccentricity of e ∼ 6.1, perihelion of q ∼ 1.36 au, inclination of ∼175°, and hyperbolic velocity of V∞ ∼ 58 km s−1. We report deep stacked images obtained using the Canada–France–Hawaii Telescope and the Very Large Telescope that resolve a compact coma. Using images obtained from several smaller ground-based telescopes, we find minimal light-curve variation for the object over a ∼4 day time span. The visible/near-infrared spectral slope of the object is 17.1% ± 0.2%/100 nm, comparable to other interstellar objects and primitive solar system small bodies (comets and D-type asteroids). Moreover, 3I/ATLAS will be observable through early 2025 September, then unobservable by Earth-based observatories near perihelion due to low solar elongation. It will be observable again from the ground in late 2025 November. Although this limitation unfortunately prohibits detailed Earth-based observations at perihelion when the activity of 3I/ATLAS is likely to peak, spacecraft at Mars could be used to make valuable observations at this time.

4 November 2021 in The Astronomical Journal
Verified
Journal Article
Closed access
MFMarco FENUCCI
Marco FENUCCI
BNBojan Novaković
Bojan Novaković
The Near-Earth asteroid (469219) Kamo’oalewa (aka 2016 HO3) is an Earth coorbital and a potential space mission target. Its short-term dynamics are characterized by a periodic switching between quasisatellite and horseshoe configurations. Due to its small diameter of only about 36 m, the Yarkovsky effect may play a significant role in the long-term dynamics. In this work, we addressed this issue by studying the changes in the long-term motion of Kamo’oalewa caused by the Yarkovsky effect. We used an estimation of the magnitude of the Yarkovsky effect assuming different surface compositions and introduced the semimajor axis drift by propagating orbits of test particles representing the clones of the nominal orbit. Our simulations showed that the Yarkovsky effect may cause Kamo’oalewa to exit from the Earth coorbital region a bit faster when compared to a purely gravitational model. Nevertheless, it still could remain an Earth companion for at least 0.5 My in the future. Our results imply that Kamo’oalewa is the most stable Earth’s coorbital object known so far, not only from a short-term perspective but also on long timescales.

9 March 2021 in Astronomy & Astrophysics
Verified
Journal Article
Closed access
PDF file
MFMarco FENUCCI
Marco FENUCCI
BNBojan Novaković
Bojan Novaković
DVDavid Vokrouhlický
David Vokrouhlický
RWRobert J. Weryk
Robert J. Weryk
Context.Asteroids with a diameter of up to a few dozen meters may spin very fast and complete an entire rotation within a few minutes. These small and fast-rotating bodies are thought to be monolithic objects because the gravitational force due to their small size is not strong enough to counteract the strong centripetal force caused by the fast rotation. This argument means that the rubble-pile structure is not feasible for these objects. Additionally, it is not clear whether the fast spin prevents dust and small particles (regolith) from being kept on their surface.Aims.We develop a model for constraining the thermal conductivity of the surface of the small, fast-rotating near-Earth asteroids. This model may suggest whether regolith is likely present on these objects.Methods.Our approach is based on the comparison of the measured Yarkovsky drift and a predicted value using a theoretical model that depends on the orbital, physical and thermal parameters of the object. The necessary parameters are either deduced from statistical distribution derived for near-Earth asteroids population or determined from observations with associated uncertainty. With this information, we performed Monte Carlo simulations and produced a probability density distribution for the thermal conductivity.Results.Applying our model to the superfast rotator asteroid (499998) 2011 PT, we find that the measured Yarkovsky drift can only be achieved when the thermal conductivityKof the surface is low. The resulting probability density function for the conductivity is bimodal, with two most likely values being around 0.0001 and 0.005 W m−1K−1. Based on this, we find that the probability thatKis lower than 0.1 W m−1K−1is at least 95%. This low thermal conductivity might indicate that the surface of 2011 PT is covered with a thermal insulating layer, composed of a regolith-like material similar to lunar dust.

3 July 2020 in Astronomy & Astrophysics
Verified
Journal Article
Open access
PDF file
Giacomo LariGiacomo Lari
Giacomo Lari
MSMelaine Saillenfest
Melaine Saillenfest
MFMarco FENUCCI
Marco FENUCCI
Context. The Galilean satellites have very complex orbital dynamics due to the mean-motion resonances and the tidal forces acting in the system. The strong dissipation in the couple Jupiter–Io is spread to all the moons involved in the so-called Laplace resonance (Io, Europa, and Ganymede), leading to a migration of their orbits. Aims. We aim to characterize the future behavior of the Galilean satellites over the Solar System lifetime and to quantify the stability of the Laplace resonance. Tidal dissipation permits the satellites to exit from the current resonances or be captured into new ones, causing large variation in the moons’ orbital elements. In particular, we want to investigate the possible capture of Callisto into resonance. Methods. We performed hundreds of propagations using an improved version of a recent semi-analytical model. As Ganymede moves outwards, it approaches the 2:1 resonance with Callisto, inducing a temporary chaotic motion in the system. For this reason, we draw a statistical picture of the outcome of the resonant encounter. Results. The system can settle into two distinct outcomes: (A) a chain of three 2:1 two-body resonances (Io–Europa, Europa–Ganymede, and Ganymede–Callisto), or (B) a resonant chain involving the 2:1 two-body resonance Io–Europa plus at least one pure 4:2:1 three-body resonance, most frequently between Europa, Ganymede, and Callisto. In case A (56% of the simulations), the Laplace resonance is always preserved and the eccentricities remain confined to small values below 0.01. In case B (44% of the simulations), the Laplace resonance is generally disrupted and the eccentricities of Ganymede and Callisto can increase up to about 0.1, making this configuration unstable and driving the system into new resonances. In all cases, Callisto starts to migrate outward, pushed by the resonant action of the other moons. Conclusions. From our results, the capture of Callisto into resonance appears to be extremely likely (100% of our simulations). The exact timing of its entrance into resonance depends on the precise rate of energy dissipation in the system. Assuming the most recent estimate of the dissipation between Io and Jupiter, the resonant encounter happens at about 1.5 Gyr from now. Therefore, the stability of the Laplace resonance as we know it today is guaranteed at least up to about 1.5 Gyr.

1 January 2022 in Serbian Astronomical Journal
Verified
Journal Article
Open access
PDF file
MFMarco FENUCCI
Marco FENUCCI
BNBojan Novakovic
Bojan Novakovic
For studies of the long-term evolution of small Solar System objects, it is fundamental to add the Yarkovsky and Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effects in the dynamical model. Still, implementations of these effects in publicly available N-body codes is either lacking, or the effects are implemented using significantly simplified models. In this paper, we present an implementation of the coupled Yarkovsky/YORP effects in the mercury and orbfit N-body codes. Along with these two effects, we also included the effects of non-destructive collisions and rotationally induced breakups to model the asteroid spin state properly. Given the stochastic nature of many incorporated effects, the software is suitable for statistical dynamical studies. Here we primarily explained the scientific aspect of the implementation, while technical details will be made freely available along with the source codes.

11 July 2023 in Astronomy & Astrophysics
Verified
Journal Article
Open access
PDF file
MFMarco FENUCCI
Marco FENUCCI
BNBojan Novaković
Bojan Novaković
DMDušan Marčeta
Dušan Marčeta
Context. Asteroids smaller than about 100 m in diameter are observed to rotate very fast, with periods often much shorter than the critical spin limit of 2.2 h. Some of these super-fast rotators can also achieve a very large semimajor axis drift induced by the Yarkovsky effect, which, in turn, is determined by internal and surface physical properties. Aims. We consider here a small super-fast-rotating near-Earth asteroid, designated as 2016 GE1. This object rotates in just about 34 s, and a large Yarkovsky effect has been determined from astrometry. By using these results, we aim to constrain the thermal inertia of the surface of this extreme object. Methods. We used a recently developed statistical method to determine the thermal properties of near-Earth asteroids. The method is based on the comparison between the observed and the modeled Yarkovsky effect, and the thermal conductivity (inertia) is determined via a Monte Carlo approach. Parameters of the Yarkovsky effect model are fixed if their uncertainty is negligible, modeled with a Gaussian distribution of the errors if they are measured, or deduced from general properties of the population of near-Earth asteroids when they are unknown. Results. Using a well-established orbit determination procedure, we determined the Yarkovsky effect on 2016 GE1 and confirm a significant semimajor axis drift rate. Using a statistical method, we show that this semimajor axis drift rate can only be explained by low thermal inertia values below 100 J m−2 K−1 s−1/2. We benchmarked our statistical method using the well-characterized asteroid Bennu and find that only knowing the semimajor axis drift rate and the rotation period is generally insufficient for determining the thermal inertia. However, when the statistical method is applied to super-fast rotators, we find that the measured Yarkovsky effect can be achieved only for very low values of thermal inertia: namely, 90% of the probability density function of the model outcomes is contained at values smaller than 100 J m−2 K−1 s−1/2. Conclusions. We propose two possible interpretations for the extremely low thermal inertia of 2016 GE1: a high porosity or a cracked surface, or a thin layer of fine regolith on the surface. Though both possibilities seem somewhat unexpected, this opens up the possibility of a subclass of low-inertia, super-fast-rotating asteroids.