# Asteroid spin-states of a 4 Gyr collisional family,

D. Athanasopoulos, J. Hanuš, Chrysa Avdellidou, R. Bonamico, M. Delbo, M.

Conjat, A. Ferrero, K. Gazeas, J. Rivet, N. Sioulas, et al.

► To cite this version:

D. Athanasopoulos, J. Hanuš, Chrysa Avdellidou, R. Bonamico, M. Delbo, et al., Asteroid spin-states of a 4 Gyr collisional family, Astronomy and Astrophysics - A&A, 2022, 666, pp.A116. 10.1051/0004-6361/202243905 hal-03854451

HAL Id: hal-03854451

https://hal.science/hal-03854451

Submitted on 16 Nov 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
# Asteroid spin-states of a 4 Gyr collisional family ★,★★

D. Athanasopoulos<fnref n="1" />,★★★, J. Hanuš², C. Avdellidou³, R. Bonamico⁴, M. Delbo³, M. Conjat³, A. Ferrero⁵,
K. Gazeas<fnref n="1" />, J. P. Rivet³, N. Sioulas⁶, G. van Belle⁷, P. Antonini⁸, M. Audejean⁹, R. Behrend¹⁰, L. Bernasconi¹¹,
J. W. Brinsfield¹², S. Brouillard¹³, L. Brunetto¹⁴, M. Fauvaud¹⁵, S. Fauvaud¹⁵, R. González¹⁶, D. Higgins¹⁷,
T. W.-S. Holoien¹⁸, G. Kober¹⁹, R. A. Koff²⁰, A. Kryszczynska²¹, F. Livet²², A. Marciniak²¹, J. Oey²³,
O. Pejcha²⁴, J. J. Rives¹³, and R. Roy²⁵

*(Affiliations can be found after the references)*

Received 29 April 2022 / Accepted 30 June 2022

## ABSTRACT

**Context.** Families of asteroids generated by the collisional fragmentation of a common parent body have been identified using clustering methods of asteroids in their proper orbital element space. However, there is growing evidence that some of the real families are larger than the corresponding cluster of objects in orbital elements, and there are families that escaped identification by clustering methods. An alternative method has been developed in order to identify collisional families from the correlation between the asteroid fragment sizes and their proper semi-major axis distance from the family centre (V-shape). This method has been shown to be effective in the cases of the very diffuse families that formed billions of years ago.

**Aims.** Here we use multiple techniques for observing asteroids to provide corroborating evidence that one of the groups of asteroids identified as a family from the correlation between size and proper semi-major axis of asteroids are real fragments of a common parent body, and thus form a collisional family.

**Methods.** We obtained photometric observations of asteroids in order to construct their rotational light curves; we combine them with the literature light curves and sparse-in-time photometry; we input these data in the light curve inversion methods, which allow us to determine a convex approximation to the 3D shape of the asteroids and their orientation in space, from which we extract the latitude (or obliquity) of the spin pole in order to assess whether an object is prograde or retrograde. We included in the analysis spin pole solutions already published in the literature aiming to increase the statistical significance of our results. The ultimate goal is to assess whether we find an excess of retrograde asteroids on the inward side of the V-shape of a 4 Gyr asteroid family identified via the V-shape method. This excess of retrograde rotators is predicted by the theory of asteroid family evolution.

**Results.** We obtained the latitude of the spin poles for 55 asteroids claimed to belong to a 4 Gyr collisional family of the inner main belt that consists of low-albedo asteroids. After re-evaluating the albedo and spectroscopic information, we found that nine of these asteroids are interlopers in the 4 Gyr family. Of the 46 remaining asteroids, 31 are found to be retrograde and 15 prograde. We also found that these retrograde rotators have a very low probability (1.29%) of being due to random sampling from an underlying uniform distribution of spin poles.

**Conclusions.** Our results constitute corroborating evidence that the asteroids identified as members of a 4 Gyr collisional family have a common origin, thus strengthening their family membership.

**Key words.** minor planets, asteroids: general – astronomical databases: miscellaneous

## 1. Introduction

The study of asteroid families has been an active field of research since the discovery of the first groupings of asteroids in orbital element space (Hirayama 1918). As more asteroids were discovered, these initial groupings became more numerous, thus substantiating their significance. At the same time, more asteroid groupings (i.e. families) were discovered. Studies of the physical properties of asteroids highlighted that the families were also homogenous in colour, albedo, and spectral properties, in general (see Masiero et al. 2015, for a review). This corroborated the idea that these groups of asteroids in orbital element space were

fragments of a common parent body. This the reason why they are also called collisional asteroid families (see Nesvorný et al. 2015, for a review on the subject).

The identification of the collisional families has been done using classical methods, such as the hierarchical clustering method (HCM; Zappala et al. 1990, 1995). Surveys of identification of asteroid collisional families have found that about one-third of the known asteroid population is associated with over 120 collisional families (see e.g. the Minor Planet Physical Properties Catalogue, MP3C)<fnref n="1" />. However, it is well known that the large majority of asteroid family identification surveys are conservative in order to clearly identify the core of the family and keep a good separation between the nearby families as well (Nesvorný et al. 2015).

A very important question is how many of the asteroids that are included in the background population (i.e. that do not belong to families) of the main belt are instead collisional members that have not been associated with known families. There is evidence

<footnote>* Table B.5 is only available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/cat/J/A+A/666/A116</footnote>

<footnote>** This article is dedicated to the memory of Gianfranco Marcon, whose telescopes have contributed to the development of astronomy and to the observations collected for this study.</footnote>

<footnote>*** Corresponding authors: D. Athanasopoulos, email: dimathanaso@phys.uoa.gr</footnote>

<footnote n="1">¹ mp3c.oca.eu</footnote>
that this number is very large, implying (i) that collisional fam-
ilies are larger than is reported in our current catalogues (Brož
& Morbidelli 2013; Tsirvoulis et al. 2018; Dermott et al. 2021)
and (ii) that there are still undiscovered families. This problem
derives from the difficulty in linking the collisional members
that currently reside far away in orbital elements to the family
core (Milani et al. 2014). It has also been shown that there is an
unexpected lack of asteroid families that formed more than 2 Gyr
years ago from a parent body larger than about 100 km (Brož &
Morbidelli 2013). Moreover, Tsirvoulis et al. (2018) attempted
to identify families using an aggressive version of the HCM,
such that it was very generous in linking large number of aster-
oids to their respective families. These authors found evidence
of undetected family members in the investigated region of the
main belt from the study of the size frequency distribution of
the non-family members. These observations also corroborate
hypotheses (i) and (ii).

As a collisional family ages, a non-gravitational force known
as the Yarkovsky effect (see e.g. the review in Vokrouhlický et al.
2015) pushes asteroids away from the centre of the family with a
drift rate *da/dt* that is proportional to the inverse diameter (1/*D*),
where *a* is the asteroid orbital semi-major axis. Prograde rotat-
ing asteroids have *da/dt* > 0 and move to a larger semi-major
axis, while retrograde asteroids, with *da/dt* < 0, move to a smaller
semi-major axis. This creates correlations of points in the (1/*D*
vs *a*) plane called V-shapes, because they resemble the letter V,
whose slope (*K*) indicates the family age. As an asteroid fam-
ily spreads in the orbital semi-major axis, family members cross
orbital resonances with planets, which perturb their eccentric-
ity and inclination. Hence, old families are less compact than
younger ones in all three proper orbital element space. This
makes old families (whose members have large orbital element
spreading) more difficult to be identified by the HCM (e.g. Bolin
et al. 2017), and cause families to overlap to each other.

Based on the V-shape characteristic, a method was devel-
oped to discover old and dispersed asteroid families (Bolin et al.
2017; Delbo et al. 2017). This V-shape identification method has
already been used to discover five families of the inner main belt:
the Eulalia and New Polana families (Walsh et al. 2013); a ‘pri-
mordial’ family with a nominal age of 4.0<sup>+1.7</sup><sub>-1.1</sub> Gyr, but that could
be as old as the Solar System; the Athor family, 3.0<sup>+0.5</sup><sub>-0.4</sub> Gyr; and
the Zita family, 5.0<sup>+1.6</sup><sub>-1.3</sub> Gyr (Delbo et al. 2017, 2019).

Although the V-shape family identification method is indeed
a powerful tool, its efficiency decreases with increasing family
age (Deienno et al. 2021) because the most dispersed family
members can no longer be distinguished from the background,
which might also consist of old and dispersed overlapping fam-
ilies. Given the current limitations of the detection methods, we
are looking for independent ways to confirm family members
and thus the borders of the families identified on the basis of
their V-shapes. One of these independent methods is to check
the spin state of the family members. According to the theory
of the Yarkovsky effect, it is expected that most of the asteroids
on the inward side (the left side of the V) of the family are retro-
grade, and prograde on the outward side (the right side of the V).
This has been shown in notable HCM asteroid families (Hanuš
et al. 2013a, 2018).

One of the most effective techniques to identify the spin state
of asteroids is the inversion of their photometric light curves
(Kaasalainen & Torppa 2001; Kaasalainen et al. 2001). This
method requires the acquisition of large datasets of photomet-
ric measurements, where most of them are currently retrieved
from all-sky surveys. However, these data are sparse in time by

nature and often do not allow the unambiguous determination
of asteroid rotational periods due to their low photometric accu-
racy of typically 0.1 magnitude. Often the addition of classical
dense-in-time optical light curves to the sparse dataset leads to
the removal of such ambiguity.

In this work we study the spin states of asteroids that
belong to the innermost border of the inward side of a pri-
mordial family of the inner main belt that was reported by
Delbo et al. (2017). This primordial family has orbital elements
(*a*, *e*, *i*) roughly (2.26, 0.14, 5.75°) for the inward wing as only
this has been identified and is suspected to contain low-albedo
C-complex asteroids. For this, we ran our own observing pro-
gramme (called *Ancient Asteroids*) and obtained dense-in-time
light curves, which we combined with existing sparse-in-time
photometric data.

In Sect. 2 we present the optical datasets used in this work
and present our observing campaign. In Sects. 3 and 4 we
describe the light curve inversion method used for the analysis
of the optical data and for the determination of asteroidal physi-
cal properties. In Sect. 5 we present our results on the spin poles
of our dataset indicating also the potential asteroid interlopers of
the primordial family, while in Sect. 6 we discuss the results.

## 2. Datasets

In order to verify membership in the primordial family, we focus
our analysis on the asteroids that are located between the inward
border of the primordial family and the respective border of the
Polana family (see Fig. 1). The latter, which is also called New
Polana after its reassessment by Walsh et al. (2013), is another
large but younger family in the inner main belt of similar car-
bonaceous composition and low albedo. We use the V-shapes of
Walsh et al. (2013) and Delbo et al. (2017) to distinguish between
the two families. This group of objects should belong either to
the vast and extensive primordial family or to the unidentified
background population. In order to study the spin poles of the
primordial family we combined a great deal of data, includ-
ing asteroid light curves that we collected from the databases,
sparse photometric data obtained from different surveys, existing
complete or incomplete shape models, and our own photometric
observations.

2.1. Currently available asteroid models

The vast majority of asteroid shape models have been produced
by the light curve inversion method (Kaasalainen & Torppa
2001; Kaasalainen et al. 2001). The Database of Asteroid Models
from Inversion Techniques (DAMIT)<fnref n="2" /> contains ~6000 asteroid
models for ~3460 asteroids (March 2022) that are publicly avail-
able (Durech et al. 2010). Shape models for 15 asteroids from our
list are already included in DAMIT (see Table B.1). Furthermore,
we also considered the partial models published in Durech et al.
(2020). For these models only the sidereal rotation period, the
ecliptic latitude of the spin axis, and its range are reported. This
information is sufficient for our purposes, as it often allows us to
securely decide whether the asteroid is a prograde or a retrograde
rotator.

## 2.2. Archival photometric data

Photometric data of asteroids are scattered throughout various
public databases or were provided to us directly by the observers.

<footnote n="2">² https://astro.troja.mff.cuni.cz/projects/damit/</footnote>
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Primordial family members presented in proper semi-major axis vs. inverse diameter plane, along with the low-albedo asteroids located in the innermost region of the main belt. Panel A: the yellow diamonds indicate members with known spin pole from the literature. The plus signs, crossed, and squares are members for which there are dense light curves only from literature, from this work, and from both sources, respectively (see details in Table B.2). Panel B: left side of the V-shape of the primordial family. The prograde asteroids are shown in red and the retrograde asteroids in blue.</figcaption></figure>

During the past two decades a large internal database of optical light curves maintained at the Institute of Astronomy of Charles University has been routinely used for shape modelling. A large amount of data was obtained from the Asteroid light curve Data Exchange Format (ALCDEF) database<fnref n="3" /> (Stephens & Warner 2018; 147 dense photometric light curves that exist for 19 of our asteroid targets). Additional dense light curves were downloaded from the Asteroid Photometric Catalogue (APC; Lagerkvist & Magnusson 2011) and the Courbes de rotation d'astéroïdes et de comètes database (CdR<fnref n="4" />) or were provided directly by the observers (see Table B.2).

The sparse-in-time data that come from various sky-surveys were obtained from corresponding databases and archives connected to the publications. We used data from the US Naval Observatory in Flagstaff (USNO-Flagstaff, IAU code 689), the Catalina Sky Survey (CSS, IAU code 703; Larson et al. 2003), Gaia Data Release 2 (GaiaDR2; Gaia Collaboration 2018), the All-Sky Automated Survey for Supernovae (ASAS-SN; Shappee et al. 2014; Kochanek et al. 2017; Hanuš et al. 2021), the Asteroid Terrestrial-impact Last Alert System (ATLAS; Tonry et al. 2018; Ďurech et al. 2020), the Zwicky Transient Facility (ZTF, IAU code 141; Bellm et al. 2019), the Palomar Transient Facility Survey (PTF; Chang et al. 2015; Waszczak et al. 2015), and TESS (Ricker et al. 2015; Pál et al. 2020). Tables B.3 and B.4 summarise the typical number of measurements from these surveys that were available for our targets. In general, data from USNO-Flagstaff, GaiaDR2, ZTF, TESS, and PTF are rather limited for our targets, while hundreds of individual measurements are available from the CSS, ASAS-SN, and ATLAS surveys. The CSS, USNO-Flagstaff, and ZTF data were obtained through the AstDys-2 database<fnref n="5" />.

## 2.3. ASAS-SN *g*-band data

As discussed in Sect. 2.2, we use *V*-band sparse data from the ASAS-SN survey through the catalogue of Hanuš et al. (2021). However, since 2018 ASAS-SN has used the SLOAN *g* filter and significantly expanded to more telescope units and sites;

we utilised these data in our work as well. We accessed and processed the *g*-band data following the same procedure as in Hanuš et al. (2021). The time coverage is already more than three years. Moreover, the *g*-band limiting magnitude is larger by about one magnitude than that of the *V*-band. Therefore, the *g*-band dataset is often comparable in terms of the number of measurements to the *V*-band dataset for brighter objects, and outperforms the *V*-band dataset for fainter objects. Both filters are treated independently in the shape modelling.

## 2.4. Ancient Asteroids: An international observing campaign

In order to enlarge our input dataset used for the shape modelling, which would potentially lead to new and improved shape solutions, we performed additional ground-based photometric observations. An international observing campaign has been initiated in the framework of our international initiative called *Ancient Asteroids*<fnref n="6" /> in order to collect dense photometric data for asteroids that belong to the oldest asteroid families (Athanasopoulos et al. 2021). *Ancient Asteroids* establishes a network of astronomers, currently from four countries, who follow a common observing plan. In the following we present the observing facilities and their corresponding equipment that participated in this work and provided data for 35 asteroids in our dataset.

*Bonamico Star Adventure Astronomical Observatory.* The Bonamico Star Adventure (BSA) is an amateur observatory located in Savigliano, Italy. For this study BSA used a robotic 0.3 m (f/8) Ritchey-Chretien telescope and an open-filter SBIG ST-9 XME CCD detector.

*Lowell Observatory.* Lowell is located in Arizona, United States, and for this project it operated two robotic telescopes: the Titan Monitor Telescope (TiMo), a 20" PlaneWave CDK20 telescope equipped with a Moravian instruments G3-6300 CCD detector, and the 1 m PlaneWave (PW1) telescope equipped with a Finger Lakes Instruments ML-16803 CCD imager. TiMo is located on Lowell's main Mars Hill campus (IAU Code 690), whereas the PW1 is on Lowell's Anderson Mesa campus

<footnote n="3">³ https://alcdef.org/</footnote>

<footnote n="4">⁴ https://obswww.unige.ch/~behrend/page3cou.html</footnote>

<footnote n="5">⁵ https://newton.spacedys.com/astdys/</footnote>

<footnote n="6">⁶ http://users.uoa.gr/~kgaze/ancient_asteroids.html</footnote>
(IAU Code 688). All the observations were performed in the Sloan r' filter.

**Bigmuskie Observatory.** Bigmuskie is an amateur observatory located in Mombercelli-Asti, Italy. Bigmuskie utilises a 40 cm (f/8.25) Ritchey-Chrétien telescope. The observations were performed unfiltered by using a Moravian G3-1000 CCD detector.

**Observatoire de la Côte d'Azur.** The observations were performed at two stations, which belong to the Observatoire de la Côte d'Azur (OCA) in France. The first station is the C2PU facility in Calern, at an altitude of a 1300 m. C2PU operated the 1.04 m Cassegrain telescope (known as Omicron@C2PU) with an f/3.2 parabolic, prime focus and with a three-lens Wynne coma corrector using an unfiltered QHY600 CMOS camera (Bendjoya et al. 2012). The second station is located on Mont Gros, the historical site of the OCA, on the east-side hills of the town of Nice, France. Mont Gros station operated the 0.4 m (f/5) diameter telescope (called Schaumasse) equipped with a QSI 583ws CCD camera.

**University of Athens Observatory.** The University of Athens Observatory (UOAO) belongs to the National and Kapodistrian University of Athens in Greece, which utilises a robotic 0.4 m (f/8) Cassegrain telescope equipped with an SBIG ST-10 XME CCD detector (Gazeas 2016). All the observations were performed unfiltered.

**Helmos Observatory.** Helmos observatory is operated by the National Observatory of Athens and is located on Mount Helmos (Aroania) in Greece, at an altitude of 2340 m. It utilises a robotic 2.3 m (f/8) Ritchey-Chrétien telescope (called Aristarchos; Goudis et al. 2010). All the observations were performed unfiltered by using the Princeton Instruments VersArray 2048B LN CCD camera.

**NOAK Observatory.** The NOAK observatory is located in the city of Ioannina, Greece. It utilises a 0.25 m (f/4.7) robotic Newtonian telescope. All the observations were performed unfiltered by using an ATIK 460EXM CCD camera.

**BlueEye 600 Observatory.** The BlueEye 600 robotic observatory (BE600) is operated by the Astronomical Institute of the Charles University and is located in Ondřejov, Czech Republic. It utilises a 60 cm Ritchey-Chrétien telescope (Officina Stellare). All the observations were performed by Martin Lehky<fnref n="7" /> utilising the standard Johnson *R* filter and the E2V42–40 CCD camera (Durech et al. 2018b).

**Pic de Château-Renard Observatory.** The Pic de Château-Renard Observatory (ChR) is located at an altitude of 2936 m in Saint-Véran in the French Alps. The facility is operated by the Paris-Meudon Observatory and AstroQueyras, an amateur association. Observations were performed unfiltered by using a 0.5 m (f/8) Ritchey-Chrétien telescope equipped with a SBIG STX 16803 camera.

**Observatoire du Bois de Bardon.** The Observatoire du Bois de Bardon (OBdB) is an amateur observatory located in Taponnat, France. OBdB used a 0.28 m (f/3) Schmidt-Cassegrain telescope equipped with a SBIG ST-402 ME CCD camera. Observations were performed with an *r'* (Sloan) filter.

**Blue Mountains Observatory.** The Blue Mountains Observatory is located at an altitude of 900 m in Leura, Australia. The

photometric observations were done with a classical Celestron Schmidt Cassegrain telescope, 0.35 m in diameter operating at f/5. All images were taken unfiltered using a SBIG CCD camera ST8-XME at bin 1 × 1.

### 3. Photometric reduction

The photometric datasets used in this work include both dense photometric data from ground-based facilities (retrieved from the literature or from our observing campaign), as well as sparse data from several sky surveys and space missions, as described above. These two different datasets require different analysis techniques.

#### 3.1. Observations in the Ancient Asteroids programme

Our observations were performed mainly in clear filter in order to increase the signal-to-noise ratio in our light curves, while keeping the exposure time as short as possible and increasing the sampling frequency. The exposure time varied between 30 s and 240 s, depending on the brightness of the target, telescope aperture, and observing conditions.

The collected data from all the observatories were reduced, following the standard image processing procedure of calibration and aperture photometry (e.g. Massey 1997; Gallaway 2020). The calibration was performed for all light frames in three steps: bias subtraction, dark subtraction, and flat-field correction. Aperture photometry is a quite simple technique and most applicable to stellar fields that are relatively sparse. This procedure was compiled by utilising *AIP4Win* software (Berry & Burnell 2005) for the fields observed by the C2PU, UOAO, Helmos, and NOAK observatories; IRIS software<fnref n="8" /> for the images performed by the OBdB and ChR observatories; and *MPO Canopus* software (Warner 2015) for the remaining images (see Table B.5).

The differential photometry was performed either with five bright field stars or by estimating an artificial comparison star, following the methodology presented by Broeg et al. (2005). The resulting measurements were provided in differential magnitudes with a photometric accuracy of 0.02–0.1 mag. In the case of OBdB and ChR, the differential photometric data were performed by estimating an artificial comparison star following the methodology described by Fauvaud & Fauvaud (2013, 2014).

We used the sigma-clipping method (see Gallaway 2020) to remove the prominent outliers in our measurements, which were usually caused by cosmic rays or satellites passing through the field. In the case where an asteroid was passing near a field star within a range of an aperture size (typically of the order of 5–7 arcseconds) we trimmed the light curve and we kept only the ‘clear’ parts.

We used the Pogson equation ($m = -2.5 \log(F) + c$) to transform the differential magnitude (*m*) to relative flux (*F*). The relative flux values were normalised by defining the average flux of each light curve as one. For all the epochs of the observed light curves we performed the light travel time correction (from the asteroid to the observer) and computed the ecliptic Cartesian coordinates (*x*, *y*, *z*) of the Sun and of the Earth, with the asteroid as the reference point, in [au] via the Miriade service (Berthier et al. 2009). This format is required by the convex inversion (CI) method that we used for the shape modelling, as described in Sect. 4.

<footnote n="7">⁷ Deceased November 18, 2020.</footnote>

<footnote n="8">⁸ http://www.astrosurf.com/buil/iris-software.html</footnote>
## 3.2. Adopted data

The dense photometric data from databases such as the Asteroid Lightcurve Data Exchange Format database (ALCDEF) or APC are in magnitude values, so we converted them to relative fluxes by following the procedure described in Sect. 3.1.

In addition to the dense photometric data, we included sparse-in-time photometric measurements from various sources as they proved to be useful in constraining the asteroid models, despite their usually low photometric accuracy of ~0.1 mag (Durech et al. 2009, 2016; Hanuš et al. 2011, 2013b). In order to use the sparse data for the shape modelling by the CI method, we processed them following the procedure of Hanuš et al. (2011). For more details, we refer to the most recent description of the procedure applied to ASAS-SN data by Hanuš et al. (2021). All individual measurements within each sparse dataset (i.e. specific survey and photometric filter) are internally calibrated; therefore, we process each dataset separately. First, the sparse data are usually available in magnitude values, which we transform into fluxes utilising the Pogson equation and, for convenience, setting the zero magnitude to 15. We then apply the light travel time correction to each epoch. Next we normalise the fluxes to a referenced one astronomical unit distance of the asteroid to the Earth and the Sun. The final steps were sigma-clipping to reject the outliers and estimating the relative weights of each sparse dataset with respect to the dense data (see Hanuš et al. 2021).

## 4. Determination of the spin poles

We used the CI method developed by Kaasalainen & Torppa (2001); Kaasalainen et al. (2001). This gradient-based inversion technique is based on shape model parametrisation by a set of facets and their normal vectors and their optimisation such that they fit to the observed light curves. Assuming a convex shape representation of the asteroid shape, the inversion problem is unique. However, adding the rotation state (i.e. sidereal rotation period and spin axis orientation) as additional free parameters, we lose the uniqueness of the solution, and the parameter space becomes full of local minima. We have to search the parameter space on a grid of input parameters and find the local minimum that corresponds to the global minimum, and thus the correct set of searched parameters. The production of the model light curve, that is compared by the method to the observed light curves, is performed by using an empirical light-scattering model, which is a combination of single Lommel-Seeliger and multiple Lambert scattering models (Kaasalainen et al. 2002). To date, the CI has been used to derive asteroid models for more than 3460 asteroids that are stored in the DAMIT database.

We assume that the shape effects in the light curves are independent of the photometric filters used while covering the reflected-dominated spectral range. Therefore, we can treat all dense light curves as relative (i.e. normalised to unity). Although each sparse dataset is, in principle, internally calibrated in a different photometric system, we also use the sparse data as normalised to unity. The only caveat is that we assume that the phase function is the same in each photometric system. This is not fully correct; for example, Durech et al. (2020) found statistically significant differences in the ATLAS data taken in the c and o filters, in accordance with the phase reddening effect (Millis et al. 1976; Lumme & Bowell 1981). However, for our purposes it is sufficient to have a single phase function for all sparse datasets.

We combined all the available datasets and applied CI to them. We weighted individual light curves and sparse datasets

based on their accuracy (expected rms). The individual weights $w_i$ were normalised such that $\sum(w_i) = N$, where N is the number of dense light curves plus the number of sparse datasets. The process that we followed is described in Hanuš et al. (in prep.).

### 4.1. Rotation period

For asteroids with previously known rotation periods we searched for the best-fit model with the period parameter varying between boundaries defined by 5, 10, and 20% of the previously reported period in the LightCurve DataBase (LCDB; Warner et al. 2009), depending on the reliability flag provided for each period estimate 3 and 3-, 2+, 2. Rotation periods for asteroids with other reliability flags were considered unknown.

Fourier-based algorithms are unable to efficiently estimate the rotation period of asteroids in cases of only sparse photometric data availability or extreme slow rotators, where their rotation period far exceeds the night duration. In these cases we performed a dense scanning in the rotation period parameter space from 2 h, which has been observed as an approximate lower limit of asteroids (>150 m in size) by Pravec et al. (2002), up to 5 000 h, which was motivated by the recent discovery of superslow rotating asteroids by Erasmus et al. (2021).

It should be noted that the period search is actually a full shape and rotation state optimisation by the CI. However, this procedure was performed on a grid of pole orientations limited to only ten values evenly distributed on a sphere and rather rough shape model resolution. We only recorded the rotation period and the best-fit rms (and $\chi^2$) value within the grid of pole directions. We considered the best-fitting period searched for on a selected period interval as unique if its $\chi^2_{min}$ is the only solution below the threshold defined as

$$ \chi^2_{\text{tr}} = \left(1 + 0.5 \sqrt{\frac{2}{\nu}}\right) \chi^2_{\text{min}}, \qquad (1) $$

where $\nu$ corresponds to the number of degrees of freedom (number of observations minus the number of free parameters).

Interestingly, the difference between two local minima in the period parameter space $\Delta P$ is dependent only on the time span of the data T and the sidereal rotation period P itself

$$ \frac{\Delta P}{P} = \frac{1}{2} \frac{P}{T}, \qquad (2) $$

and thus, for instance, independent of the shape resolution. We sample the period with a step of $0.5\Delta P$ that we recompute every step. We note that ~1/20 of $\Delta P$ is a reliable uncertainty for the derived rotation period as it corresponds to a rotation phase offset of about 10°, which is a typical uncertainty on the pole direction.

### 4.2. Spin pole and shape

After the determination of the rotation period, we applied the CI method on a much denser grid of initial spin directions (ecliptic longitude $\lambda$ and latitude $\beta$) and a higher shape resolution<fnref n="9" />. The best-fitting shape and spin pole solution is considered unique if

<footnote n="9">⁹ For the search for the rotation period we fit $\gamma^2$ coefficients of the spherical harmonic expansion and start with ten different pole directions uniformly distributed on the sphere. The denser grid utilises $\gamma^2$ coefficients of the expansion and 48 input pole directions. The typical polyhedron representing the final shape solution contains about 1000 vertices and 2000 facets.</footnote>
it fulfils the condition of Eq. (1). Most of the cases presented two symmetrical solutions with respect to ecliptic longitude (λ±180°), the so-called pole ambiguity (Kaasalainen & Lamberg 2006).

If photometric data are rich enough in observing geometries and have reasonable accuracy relative to the amplitude of the asteroid light curve, the unique solution can be often derived. Sometimes, however, the photometric data are insufficient to derive a unique solution, but they allow us to derive a unique period and three or four pole solutions that have similar values of the ecliptic latitude β. These so-called partial models are still useful, especially for our study, as it is possible to decide whether the asteroid is a prograde or a retrograde rotator. Therefore, we made an effort to identify such cases here together with the unique solutions.

# 5. Results

## 5.1. New period estimates

We derived improved values of the sidereal rotation period for 29 asteroids, as presented in Tables B.3 and B.4 (see also Fig. A.1 for an example of typical periodograms). All these period values are consistent with the synodic periods reported in the LCDB database<fnref n="10" />. The revised asteroid models have almost the same period as the previous solutions. Moreover, we measured rotation periods for seven asteroids for the first time. These periods were derived using sparse photometric data. Of these asteroids one is a super-slow rotator with *P* = 3253.5 h and one a slow rotator with *P* = 152.62 h, which are (2776) Baikal and (8315) Bajin, respectively. All the others, namely (12722) Petrarca, (13066) 1991 PM13, (23495) 1991 UQ1, (49863) 1999 XK 104, and (70184) 1999 RU3, are rather fast rotators, with periods between 3 and 10 h.

## 5.2. Spin pole directions

By combining new and literature data, we successfully determined the shapes and spin states for 55 asteroids that belong to the nominal population of the primitive primordial family of the inner main belt (Delbo et al. 2017). This corresponds to 51% of the population in the sliver between the left-wing border of the Polana family and the primordial family (see Fig. 1). In particular, we calculated 25 new complete asteroid models, 20 revised and 4 new partial models (see Tables B.3 and B.4). Specifically, 34 asteroids have retrograde rotation, while 21 are prograde. The case of (2575) Bulgaria is presented as an example of model fits to the dense and sparse photometric data in Figs. A.2 and A.3, respectively. Figure A.4 shows an example of four shape models derived from our analysis.

## 5.3. Identification of interlopers

In order to proceed and study the spin states of the inner main belt primordial family, we need to eliminate any potential interlopers beforehand. These are defined as objects that may reside inside the V-shape of a family and/or that are grabbed by the HCM methods in the same cluster, but their physical properties, such as the geometric visible albedo *p*<sub>V</sub> and spectral class,

are totally different from the bulk composition of the family (Nesvorný et al. 2015). These asteroids that are clearly non-family members, may belong to the background population or to other nearby families.

In our case seven asteroids have reported spectral (visible or near-infrared) or spectrophotometric (e.g. SDSS, ECAS, and MOVIS) data that do not match the primordial family composition, which is composed primarily of dark albedo and featureless spectra. Specifically, asteroids (1806) Derice, (2171) Kiev, (2575) Bulgaria, and (6125) Singto are classified as S-complex asteroids; (2768) Gorky as an A-type; and (5524) Lecacheux and (15415) Rika as V-types in any of the main classification schemes of Tholen (Tholen 1989), Bus (Bus & Binzel 2002), and Bus&DeMeo (DeMeo et al. 2009). Of these, asteroids (1806), (2171), (2575), and (2768) had been already assigned to the nearby bright family of Flora that consists of S-complex asteroid members, which leaves no doubts that they are interlopers of the low-albedo primordial family of Delbo et al. (2017). Although asteroids (2536) Kozyrev and (2705) Wu have no spectral or spectrophotometric information, their *p*<sub>V</sub> values are moderate and beyond the 12% that is generally defined as the separation between the S and C spectroscopic complexes (Delbo et al. 2017). In the studied sample there are seven asteroids (220, 428, 917, 1244, 1544, 1700, 3633) whose *p*<sub>V</sub> is in agreement with the primordial family but are classified as members of the X-complex. Using only the visible part of the spectrum (or spectrophotometric data) is not always sufficient to distinguish the slope between the X- and C-complex asteroids, the latter being the main components of the primordial family. So, although these seven asteroids could indeed be interlopers, at this stage we cannot definitely exclude them from the family.

On the other hand, in the population that is studied in this work there are asteroids that had been previously assigned using the clustering methods to the other nearby families of Flora and Vesta (Nesvorný et al. 2015), such as the asteroids (428) Monachia, (4524) Barklajdetolli, (2839) Annette, and (3633) Mira. However, their low (<10%) *p*<sub>V</sub> values and their featureless spectrophotometric data are in contrast to this assignment, and therefore remain in the primordial family.

The above analysis indicated nine interlopers in the sample of 55 studied objects. From these 46 confirmed asteroid members of the primordial family, 31 asteroid models (67%) have retrograde rotation and 15 prograde, including the partial solutions. Excluding the seven X-complex asteroids the abundance of retrograde asteroids reaches 72%. The distribution of the sense of rotation (i.e. prograde or retrograde) within the family is presented in Fig. 1 (Panel B). Additionally, Fig. 2 illustrates the distribution of spin axis directions. Table B.6 presents all the objects studied in this work, indicating their prograde or retrograde spin along with their physical properties. Diameters (*D*) and geometric albedos (*p*<sub>V</sub>) were calculated as the weighted averages of all the available measurements in the literature and were retrieved from the Minor Planet Physical Properties Catalogue. All uncertainty-weighted averages use $1/\sigma^2$ as weights, where $\sigma$ is the error of each measurement.

# 6. Discussion

## 6.1. The cases of (2171) Kiev, (7132) Casuli and (2705) Wu

Most of the asteroids studied in our sample are found to be single objects with no evidence of any close companion body or satellite. However, two asteroids in our sample host a satellite, while another is probably a non-principal-axis (NPA) rotator.

<footnote n="10">¹⁰ The only exception is asteroid (1159) Granada with the reported LCDB period of 77.28 h. However, this period adopted from the CdR database is preliminary and inconsistent with the period reported, for example, by Waszczak et al. (2015), which is similar to the value we derived.</footnote>
<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Distribution of the complete spin pole solutions for the primordial family members. The main plot is a sinusoidal equal-area cartographic representation, where the vertical grey lines define the longitude (λ) and the horizontal curves define the latitude (β). The right histogram represents the latitude (β) of prograde (red) and retrograde (blue) rotators. The bottom histogram represents the longitude (λ) of prograde (red) and retrograde (blue) rotators. The specific values for each asteroid are shown in Tables B.3 and B.4.</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Histogram of spin rate for members of the primordial asteroid family. The prograde rotators are in red and the retrograde in blue. The light grey bars represent the primordial family members whose rotational period is known either from this study or the literature.</figcaption></figure>

The first two cases are (2171) Kiev and (7132) Casuli, which have been reported as binaries by Loera-González et al. (2020) and Franco et al. (2020a), respectively. Our dense photometric data confirm that (2171) Kiev and (7132) Casuli are indeed binary asteroids; however, the observed light curves are not enough to further constrain the orbits of secondaries. In each case the eclipsing part of the light curve was removed and the CI method was applied to the rotational light curves of the primary body. Thus, the results presented in Tables B.6 and B.3 are only for the primary bodies. Although the derived parameters do not deviate from the rest of the sample, care should be taken given that the companion body can alter the spin axis and rotational characteristics of the primary body.

The third case, (2705) Wu, is a slow retrograde rotator, as our analysis and that of Durech et al. (2020) have shown by using sparse photometric data. Previous dense-in-time observations have shown that it is possibly a NPA rotator (i.e. tumbler, Oey 2010). As has been noted, some deviations from the single periodicity are clearly seen, but not at a conclusive level, while more photometric data are needed to resolve the second period. In this study no further dense-in-time observations were obtained. The characteristic timescale of damping of the excited NPA rotation can be estimated as $\tau_d = P_{[h]}/(C^3 \cdot D_{[km]}^2)$ [Gyr] by Harris (1994), where $C = 17 \pm 2.5$ Gyr $\cdot$ km$^2$/h$^3$. For this asteroid, the damping timescale is estimated to be $12 \pm 5$ Gyr, which is greater than the age of our Solar System. This is statistically common for NPAs, which have a diameter larger than ~0.4 km (Pravec et al. 2005). So, if (2705) Wu is a tumbler, the spin solution of our study related to the sense of rotation can be trusted.

## 6.2. Distribution of the spin poles

All asteroid models with retrograde solutions, except the solutions for (933) Susi and (49863) 1999 XK104, have large ecliptic pole latitude values $|\beta| \ge 30^\circ$ with a large predominance towards the YORP end state values approaching $\beta \sim -90^\circ$. The latitude distribution for prograde rotators differs slightly from the retrograde rotators by having more values with $|\beta| \le 60^\circ$. This is

likely due to various resonances acting only on prograde rotators. Similar behavior is also observed in other asteroid families (Hanusš et al. 2013a).

The distribution of ecliptic pole longitudes is bi-modal for prograde rotators (with two peaks, at ~45° and ~225°) and irregular for retrograde asteroids (with two peaks, at ~75° and ~255°). Previous studies have estimated that the longitude distribution for main belt asteroids is uniform with no statistically significant features (Kryszczyńska et al. 2007; Hanuš et al. 2011). On the contrary, more recent studies present an anisotropic longitude distribution with two symmetrical maxima around ~50° and ~230° and minima around ~140° and ~320° (Bowell et al. 2014; Cibulková et al. 2016).

Retrograde and prograde asteroids also have different period distributions. As Fig. 3 shows, the majority of slow rotators are retrograde. Moreover, the periods of retrograde asteroids have a non-Maxwellian distribution with excesses at the fast and the slow rotations. This bimodal distribution is in agreement with that of main belt asteroids <40 km, as a result of the YORP effect (Pravec & Harris 2000; Pravec et al. 2002). The period distribution of prograde asteroids is irregular, with small peaks for slow, moderate, and fast rotations. Thus, prograde and retrograde seem to have different spin rate distributions. A simple model by Pravec et al. (2008) estimated that a uniform distribution for <40 km main belt asteroids could happen on large timescales. Moreover, the different spin rate distributions could signify a different YORP evolution (Pravec & Harris 2000; Pravec et al. 2002, 2008, and references therein).

## 6.3. Statistical predominance of the retrograde spin poles

It is possible to test whether the observed predominance of the retrograde spin poles could be due to chance, created by random sampling of an equal-probability population of prograde and retrograde asteroids. In particular, we test the probability of obtaining 15 or fewer prograde rotators from 46 observed asteroids, drawing from a population having an equal probability (0.5) of being retrograde or prograde:

$$p(\le 15, 46) = \sum_{j=1}^{15} \binom{46}{j} 0.5^j (1 - 0.5)^{46-j}. \quad (3)$$
An evaluation of Eq. (3) gives $p(\le 15, 46) = 1.29\%$. This shows that our observations rule out the null hypothesis (i.e. there is no statistical predominance of the retrograde spin poles) at 98.71% probability.

## 6.4. The YORP reorientation timescale for the primordial family members

Although a predominance of retrograde spinning asteroids is expected in the inward wing of the V-shape of a collisional family, it is also likely to find prograde rotators. During the family lifespan of a few billion years, there are a number of processes that can re-orient the spin vector on an asteroid, such as non-catastrophic collisions, activity, or topographic changes, for example due to mass wasting and/or movement (see e.g. Paolicchi & Knežević 2016; Bottke et al. 2015, for a discussion). Torques on asteroids due to the unbalanced emitted and reflected radiation (which cause the YORP effect Rubincam 2000) are very sensitive to topographic features, and so are the strength and sign of the YORP effect (which can be reversed; Statler 2009). Statler (2009) theoretically showed that YORP is sensitive to small surface features, and hence even a small impact or mass movement could alter the shape sufficiently and change the sense of pole-period evolution. Moreover, as the YORP effect drives some of the family members towards the critical threshold for fast rotation, topographic instability might easily occur, leading to new YORP coefficients, which can drive the spin-pole in the opposite direction. The older a family is, the higher the cumulative probability is that some family members may have undergone shape changes or received non-catastrophic impacts, hence the higher the probability of detecting prograde rotators in the wing of the V-shape with a predominance of retrograde rotators.

Given the above, we used the model of Vokrouhlický et al. (2006) to estimate the evolution of the spin vectors of the family member asteroids. The model of Vokrouhlický et al. (2006) takes into account the change in orbital semi-major axis due to the Yarkovsky effect and the changes in the rotation rate $\frac{d\omega}{dt}$ and spin axis obliquity<fnref n="11" /> $\omega\frac{d\epsilon}{dt}$ of the spin axis due to the YORP torques. In particular, the values of $\frac{d\omega}{dt}$ and $\omega\frac{d\epsilon}{dt}$ are multiplied by a constant named $C_{YORP}$, which Bottke et al. (2015) have proposed to be between 0.5 and 0.7. However, Vokrouhlický et al. (2006) found that values of $C_{YORP} \sim 1$ could produce model asteroid families that well represent the real ones. Hence, for simplicity, we assumed here that $C_{YORP} = 1$. The values of $\frac{d\omega}{dt}$ and $\omega\frac{d\epsilon}{dt}$ are functions of the orbital and physical parameters of asteroids, as described by Vokrouhlický et al. (2006), which we follow hereafter.

Following Vokrouhlický et al. (2006), we modelled the evolution of the 46 family members for which we have a pole solution. We assumed they all started with an initial retrograde spin direction, which we randomly assigned uniformly between 90° and 180° obliquity. We initialised model asteroids with the known diameters and current spin periods. We evolved the model with a time step of 10 Myr for 4 Gyr. At each time step the values of $a$, $\omega$, and $\epsilon$ were updated by summing their respective time derivatives multiplied by the step in time. At each time step we also evaluated the probability that an asteroid could suffer a non-catastrophic collision capable of changing its spin axis. This probability is given by the 10 Myr time step divided by the

re-orientation timescale, $\tau_{reor}$, which is estimated as

$$ \tau_{reor} = 0.845 \left(\frac{5}{P_{[h]}}\right)^{5/6} \left(\frac{D_{[km]}}{2}\right)^{4/3} \text{ [Gyr]}. \quad (4) $$

Then, for each asteroid (at each time step), we extracted a random number, uniformly distributed between 0 and 1, and we re-oriented the spin state of the asteroid when the random number was smaller than said probability. The spin axis reorientation is performed by picking a new random direction of the obliquity uniformly distributed between 0° and 180°.

After the model was completed, we counted how many of the initially retrograde asteroids became prograde rotators due to spin evolution. We ran the model 10 000 times and we found that the probability of having 15 or fewer prograde rotators of the initial 46 retrograde objects is 13.5%; on average, we found 11.4 prograde asteroids.

However, Delbo et al. (2017) showed that $D \ge 35$ km asteroids could be primordial objects that accreted as planetesimals from the dust of our protoplanetary disk. It is therefore possible that some of them are within the inward wing of the V-shape of the primordial family. Hence, we also considered the above model only for $D < 35$ asteroids, of which 14 are prograde rotators and 29 retrograde rotators. In this case we found a 20% probability of having 14 or fewer prograde rotators of the initial 43 retrograde objects is 20%.

We can conclude that it is possible to observe the current mix of retrograde and prograde rotators within the inward wing of the 4 Gyr collisional family.

## 7. Conclusions

We carried out a campaign of photometric observations of those asteroids that have been classified as members of one of the oldest collisional (primordial) families in the Solar System (Delbo et al. 2017). We constructed photometric time series, the light curves, for 49 asteroids. This corresponds to 46% of the members of the primordial family. We combined our light curves with those from the literature, and with sparse-in-time photometry in order to create multi-epoch photometric datasets to be used as inputs for the convex inversion method. We obtained 49 new and revised shape models and their spin vector solutions. We combined this with the literature spin vectors (for six objects).

We reassessed the albedo values for the observed asteroids and studied their literature spectra. This allowed us to find nine interlopers among the initial list of family members of Delbo et al. (2017). After removing these interlopers, we find that 31 and 15 out of the remaining 46 asteroids are retrograde and prograde, respectively. We show that this predominance of retrograde compared to prograde asteroids is very unlikely (1.29% probability) to be due to sampling a distribution of objects with equal probability of being prograde and retrograde. This corroborates the hypothesis that the statistical predominance of the retrograde spin poles is due to a physical process, as was claimed by Delbo et al. (2017), namely formation as collisional fragments of a common parent body, a subsequent dynamical evolution driven by the Yarkovsky effect.

*Acknowledgements.* M.D. and C.A. acknowledge support from ANR "ORIGINS" (ANR-18-CE31-0014). This work is based on data provided by the Minor Planet Physical Properties Catalogue (MP3C) of the Observatoire de la Côte d'Azur. The research of JH has been supported by the Czech Science Foundation through grant 20-08218S. The work of OP has been supported by INTER-EXCELLENCE grant LTAUSA18093 from the Ministry of Education, Youth, and Sports. Support for T.W.-S.H. was provided by NASA through the NASA Hubble Fellowship

<footnote n="11">¹¹ Obliquity ($\epsilon$) is defined as the angle between the equatorial and orbital planes of an asteroid. For small orbital inclinations, an obliquity of 0° is ~90°, and 180° is ~90°.</footnote>
grant HST-HF2-51458.001-A awarded by the Space Telescope Science Institute (STScI), which is operated by the Association of Universities for Research in Astronomy, Inc., for NASA, under contract NAS5-26555. We thank the Las Cumbres Observatory and their staff for its continuing support of the ASAS-SN project. ASAS-SN is supported by the Gordon and Betty Moore Foundation through grant GBMF5490 to the Ohio State University, and funded in part by the Alfred P. Sloan Foundation grant G-2021-14192 and NSF grant AST-1908570. Development of ASAS-SN has been supported by NSF grant AST-0908816, the Mt. Cuba Astronomical Foundation, the Center for Cosmology and Astrophysics Physics at the Ohio State University, the Chinese Academy of Sciences South America Center for Astronomy (CAS-SACA), the Villum Foundation, and George Skestos.

## References

Alvarez-Candal, A., Duffard, R., Lazzaro, D., & Michtchenko, T. 2006, A&A, 459, 969

Athanasopoulos, D., Bonamico, R., Van Belle, G., et al. 2021, Eur. Planet. Sci. Cong., 15, 335

Bellm, E. C., Kulkarni, S. R., Graham, M. J., et al. 2019, PASP, 131, 018002

Bendjoya, P., Abe, L., Rivet, J. P., et al. 2012, in SF2A-2012: Proceedings of the Annual meeting of the French Society of Astronomy and Astrophysics, eds. S. Boissier, P. de Laverny, N. Nardetto, R. Samadi, D. Valls-Gabaud, & H. Wozniak, 643

Berry, R., & Burnell, J. 2005, The Handbook of Astronomical Image Processing (USA: Willmann-Bell), 2

Berthier, J., Hestroffer, D., Carry, B., et al. 2009, Eur. Planet. Sci. Cong., 2009, 676

Bolin, B. T., Delbo, M., Morbidelli, A., & Walsh, K. J. 2017, Icarus, 282, 290

Bonamico, R., & van Belle, G. 2021, Minor Planet Bulletin, 48, 210

Bottke, W. F., Vokrouhlický, D., Walsh, K. J., et al. 2015, Icarus, 247, 191

Bowell, E., Oszkiewicz, D. A., Wasserman, L. H., et al. 2014, Meteorit. Planet. Sci., 49, 95

Broeg, C., Fernández, M., & Neuhäuser, R. 2005, Astron. Nachr., 326, 134

Brož, M., & Morbidelli, A. 2013, Icarus, 223, 844

Bus, S. J., & Binzel, R. P. 2002, Icarus, 158, 146

Carvano, J. M., Hasselmann, P. H., Lazzaro, D., & Mothé-Diniz, T. 2010, A&A, 510, A43

Chang, C.-K., Ip, W.-H., Lin, H.-W., et al. 2015, ApJS, 219, 27

Cibulková, H., Dürech, J., Vokrouhlický, D., Kaasalainen, M., & Oszkiewicz, D. A. 2016, A&A, 596, A57

Cooney, W., Benishek, V., Pravec, P., et al. 2017, Central Bureau Electronic Telegrams, 4401, 1

Deienno, R., Walsh, K. J., & Delbo, M. 2021, Icarus, 357, 11428

Delbo, M., Walsh, K., Bolin, B., Avdellidou, C., & Morbidelli, A. 2017, Science, 357, 1026

Delbo, M., Avdellidou, C., & Morbidelli, A. 2019, A&A, 624, A69

DeMeo, F. E., Binzel, R. P., Slivan, S. M., & Bus, S. J. 2009, Icarus, 202, 160

Dermott, S. F., Li, D., Christou, A. A., et al. 2021, MNRAS, 505, 1917

Dürech, J., Kaasalainen, M., Warner, B. D., et al. 2009, A&A, 493, 291

Durech, J., Sidorin, V., & Kaasalainen, M. 2010, A&A, 513, A46

Dürech, J., Hanuš, J., Oszkiewicz, D., & Vančo, R. 2016, A&A, 587, A48

Dürech, J., Hanuš, J., & Alí-Lagoa, V. 2018a, A&A, 617, A57

Dürech, J., Hanuš, J., Brož, M., et al. 2018b, Icarus, 304, 101

Dürech, J., Hanuš, J., & Vančo, R. 2019, A&A, 631, A2

Dürech, J., Tonry, J., Erasmus, N., et al. 2020, A&A, 643, A59

Erasmus, N., Navarro-Meza, S., McNeill, A., et al. 2020, ApJS, 247, 13

Erasmus, N., Kramer, D., McNeill, A., et al. 2021, MNRAS, 506, 3872

Fauvaud, S., & Fauvaud, M. 2013, Minor Planet Bull., 40, 224

Fauvaud, S., & Fauvaud, M. 2014, Minor Planet Bull., 41, 143

Ferrero, A. 2021, Minor Planet Bull., 48, 215

Franco, L., Marchini, A., Bonnoli, G., et al. 2020a, ATel, 13590, 1

Franco, L., Marchini, A., Saya, L.-F., et al. 2020b, Minor Planet Bull., 47, 242

Gaia Collaboration, (Spoto, F., et al.) 2018, A&A, 616, A13

Gallaway, M. 2020, in An Introduction to Observational Astrophysics (Berlin: Springer), 169

Gazeas, K. 2016, Rev. Mex. Astron. Astrofis. Conf. Ser., 48, 22

Goudis, C., Hantzios, P., Boumis, P., et al. 2010, ASP Conf. Ser., 424, 422

Hanuš, J., Dürech, J., Brož, M., et al. 2011, A&A, 530, A134

Hanuš, J., Brož, M., Dürech, J., et al. 2013a, A&A, 559, A134

Hanuš, J., Dürech, J., Brož, M., et al. 2013b, A&A, 551, A67

Hanuš, J., Dürech, J., Oszkiewicz, D. A., et al. 2016, A&A, 586, A108

Hanuš, J., Delbo, M., Alí-Lagoa, V., et al. 2018, Icarus, 299, 84

Hanuš, J., Pejcha, O., Shappee, B. J., et al. 2021, A&A, 654, A48

Harris, A. W. 1994, Icarus, 107, 209

Hirayama, K. 1918, AJ, 31, 185

Kaasalainen, M., & Lamberg, L. 2006, Inv. Prob., 22, 749

Kaasalainen, M., & Torppa, J. 2001, Icarus, 153, 24

Kaasalainen, M., Torppa, J., & Muinonen, K. 2001, Icarus, 153, 37

Kaasalainen, M., Mottola, S., & Fulchignoni, M. 2002, in Asteroids III (Tucson: University of Arizona press), 139

Kochanek, C. S., Shappee, B. J., Stanek, K. Z., et al. 2017, PASP, 129, 104502

Kryszczyńska, A., La Spina, A., Paolicchi, P., et al. 2007, Icarus, 192, 223

Lagerkvist, C. I., & Magnusson, P. 2011, NASA Planetary Data System, EAR

Larson, S., Beshore, E., Hill, R., et al. 2003, BAAS, 35, 982

Lazzaro, D., Angeli, C. A., Carvano, J. M., et al. 2004, Icarus, 172, 179

Loera-González, P., Olguín, L., Saucedo-Morales, J., & Nuñez-López, R. 2020, Minor Planet Bull., 47, 348

Lumme, K., & Bowell, E. 1981, AJ, 86, 1705

Mahlke, M., Carry, B., & Denneau, L. 2021, Icarus, 354, 114094

Masiero, J. R., Carruba, V., Mainzer, A., Bauer, J. M., & Nugent, C. 2015, ApJ, 809, 179

Massey, P. 1997, National Optical Astronomy Observatory (Tucson: University of Arizona Press)

Milani, A., Cellino, A., Knežević, Z., et al. 2014, Icarus, 239, 46

Millis, R. L., Bowell, E., & Thompson, D. T. 1976, Icarus, 28, 53

Morate, D., de León, J., De Prá, M., et al. 2019, A&A, 630, A141

Nesvorný, D., Brož, M., & Carruba, V. 2015, Asteroids IV, Identification and Dynamical Properties of Asteroid Families, eds. P. Michel, F. E. DeMeo, & W. F. Bottke (Tucson: University of Arizona Press), 297

Oey, J. 2010, Minor Planet Bull., 37, 53

Pál, A., Szakáts, R., Kiss, C., et al. 2020, ApJS, 247, 26

Paolicchi, P., & Knežević, Z. 2016, Icarus, 274, 314

Polakis, T. 2021, Minor Planet Bull., 48, 158

Popescu, M., Licandro, J., Carvano, J. M., et al. 2018, A&A, 617, A12

Pravec, P., & Harris, A. W. 2000, Icarus, 148, 12

Pravec, P., Harris, A. W., & Michalowski, T. 2002, in Asteroids III (Tucson: University of Arizona Press), 113

Pravec, P., Harris, A. W., Scheirich, P., et al. 2005, Icarus, 173, 108

Pravec, P., Harris, A. W., Vokrouhlický, D., et al. 2008, Icarus, 197, 497

Ricker, G. R., Winn, J. N., Vanderspek, R., et al. 2015, J. Astron. Teles. Instrum. Syst., 1, 014003

Rubincam, D. P. 2000, Icarus, 148, 2

Shappee, B. J., Prieto, J. L., Grupe, D., et al. 2014, ApJ, 788, 48

Statler, T. S. 2009, Icarus, 202, 502

Stephens, R. D. 2011, Minor Planet Bull., 38, 23

Stephens, R., & Warner, B. D. 2018, AAS/Div. Planet. Sci. Meeting Abs., 50, 417.03

Stephens, R. D., & Warner, B. D. 2019, Minor Planet Bull., 46, 449

Stephens, R. D., & Warner, B. D. 2020, Minor Planet Bull., 47, 224

Tholen, D. J. 1989, in Asteroids II, eds. R. P. Binzel, T. Gehrels, & M. S. Matthews (University of Arizona Press), 1139

Tonry, J. L., Denneau, L., Heinze, A. N., et al. 2018, PASP, 130, 064505

Tsirvoulis, G., Morbidelli, A., Delbo, M., & Tsiganis, K. 2018, Icarus, 304, 14

Vokrouhlický, D., Brož, M., Bottke, W. F., Nesvorný, D., & Morbidelli, A. 2006, Icarus, 182, 118

Vokrouhlický, D., Bottke, W. F., Chesley, S. R., Scheeres, D. J., & Statler, T. S. 2015, in Asteroids IV, eds. P. Michel, F. E. DeMeo, & W. F. Bottke (The University of Arizona Press), 509

Walsh, K. J., Delbó, M., Bottke, W. F., Vokrouhlický, D., & Lauretta, D. S. 2013, Icarus, 225, 283

Warner, B., D. 2015, MPO Canopus software

Warner, B. D., Harris, A. W., & Pravec, P. 2009, Icarus, 202, 134

Waszczak, A., Chang, C.-K., Ofek, E. O., et al. 2015, AJ, 150, 75

Zappala, V., Cellino, A., Farinella, P., & Knezevic, Z. 1990, AJ, 100, 2030

Zappala, V., Bendjoya, P., Cellino, A., Farinella, P., & Froeschlé, C. 1995, Icarus, 116, 291

Zellner, B., Tholen, D. J., & Tedesco, E. F. 1985, Icarus, 61, 355

<footnote n="1">¹ Section of Astrophysics, Astronomy and Mechanics, Department of Physics, National and Kapodistrian University of Athens, Zografos GR 15784, Athens, Greece</footnote>

<footnote n="2">² Charles University, Faculty of Mathematics and Physics, Institute of Astronomy, V Holešovičkách 2, CZ-18000, Prague, Czech Republic e-mail: josef.hanus@mff.cuni.cz</footnote>
³ Université Côte d'Azur, Observatoire de la Côte d'Azur, CNRS, Laboratoire Lagrange, Nice, France

⁴ BSA Osservatorio (K76), Strada Collarelle 53, 12038 Savigliano, Cuneo, Italy

⁵ Bigmuskie Observatory (B88), via Italo Aresca 12, 14047 Moberelli, Asti, Italy

⁶ NOAK Observatory (L02), Stavraki Ioannina, Greece

⁷ Lowell Observatory, 1400 West Mars Hill Road, Flagstaff, AZ 86001, USA

⁸ Observatoire des Hauts Patys, 84410 Bédoin, France

⁹ Observatoire de Chinon, Mairie de Chinon, 37500 Chinon, France

¹⁰ Geneva Observatory, CH-1290 Sauverny, Switzerland

¹¹ Observatoire des Engarouines, 1606 chemin de Rigoy, 84570 Maletort-du-Comtat, France

¹² Via Capote Observatory, Thousand Oaks, CA 91320, USA

¹³ AstroQueyras, 05530 Saint-Véran, France

¹⁴ Le Florian, Villa 4, 880 chemin de Ribac-Estagnol, 06600 Antibes, France

¹⁵ Observatoire du Bois de Bardon, 16110 Taponnat, France

¹⁶ Observatorio Uraniborg, 41400 Écija, Sevilla, Spain

¹⁷ Hunters Hill Observatory, 7 Mawalan Street, Ngunnawal ACT 2913, Australia

¹⁸ The Observatories of the Carnegie Institution for Science, 813 Santa Barbara St., Pasadena, CA 91101, USA; NHFP Einstein Fellow

¹⁹ Observatoire de Durtal, 49430 Durtal, France

²⁰ Antelope Hills Observatory, 980 Antelope DR W, Bennett, CO 80102 USA

²¹ Astronomical Observatory Institute, Faculty of Physics, A. Mickiewicz University, Słoneczna 36, 60-286 Poznań, Poland

²² Institut d'Astrophysique de Paris, 98 bis boulevard Arago, UMR 7095 CNRS et Sorbonne Universités, 75014 Paris, France

²³ Blue Mountains Observatory, Leura, Australia

²⁴ Charles University, Faculty of Mathematics and Physics, Institute of Theoretical Physics, V Holešovičkách 2, 18000 Prague, Czech Republic

²⁵ Observatoire de Blauvac, 293 chemin de St Guillaume, 84570 Blauvac, France
**Appendix A: Supplementary figures**

<figure><img src="image_5.png" /><figcaption>**Fig. A.1.** Typical periodograms of six different asteroids and their new shape model determinations. Each dot represents one trial run that samples all the local minima at a fixed rotation period (Eq. 2) within the searched interval. The vertical lines indicate the best-fit values. The horizontal line represents the $\chi^2$ threshold defined by Eq. 1.</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>**Fig. A.2.** Example of fits to the dense light curves for asteroid (2575) Bulgaria. All five light curves were obtained in the framework of this study (observed by Andrea Ferrero).</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>**Fig. A.3.** Example of fits to the sparse datasets for asteroid (2575) Bulgaria. The function $f(\alpha)$ is our fit of a semi-empirical phase function (Kaasalainen et al. 2002). The datasets (from top to bottom) are ASAS-SN V-band, ASAS-SN g-band, ATLAS c-band, and ATLAS o-band.</figcaption></figure>
<figure><img src="image_8.png" /><figcaption>**Fig. A.4.** Examples of several shape models. The first panel (top) is the equatorial view with a 90° rotation offset; the third panel is north pole-on view. **The shape models correspond to the first pole solution ($\lambda_1, \beta_1$) for each asteroid.** The asteroids (from top to bottom) are (1159) Granada, (1700) Zvezdara, (2575) Bulgaria, and (2792) Ponomarev.</figcaption></figure>
## Appendix B: Tables

<table><caption>**Table B.1.** Asteroids whose spin poles are known from the literature.</caption><thead><tr><th>Number</th><th>Asteroid Name/Designation</th><th>P (h)</th><th>λ<sub>1</sub> (deg)</th><th>β<sub>1</sub> (deg)</th><th>λ<sub>2</sub> (deg)</th><th>β<sub>2</sub> (deg)</th><th>Original model publication</th></tr></thead><tbody><tr><td>220</td><td>Stephania</td><td>18.209</td><td>26</td><td>-50</td><td>223</td><td>-62</td><td>Hanuš et al. (2013b)</td></tr><tr><td>249</td><td>Ilse</td><td>84.995</td><td>2</td><td>85</td><td>222</td><td>41</td><td>Hanuš et al. (2016)</td></tr><tr><td>282</td><td>Clorinde</td><td>49.36</td><td>353</td><td>-66</td><td>184</td><td>-47</td><td>Ďurech et al. (2020)</td></tr><tr><td>370</td><td>Modestia</td><td>22.5411</td><td>-</td><td>-50 ± 9</td><td>-</td><td>-</td><td>Ďurech et al. (2020)</td></tr><tr><td>428</td><td>Monachia</td><td>3.63360</td><td>-</td><td>51 ± 9</td><td>-</td><td>-</td><td>Ďurech et al. (2020)</td></tr><tr><td>933</td><td>Susi</td><td>4.6224</td><td>301</td><td>-10</td><td>125</td><td>-15</td><td>Ďurech et al. (2020)</td></tr><tr><td>1216</td><td>Askania</td><td>6.53713</td><td>-</td><td>44 ± 14</td><td>-</td><td>-</td><td>Ďurech et al. (2020)</td></tr><tr><td>1244</td><td>Deira</td><td>216.98</td><td>107</td><td>-56</td><td>314</td><td>-46</td><td>Hanuš et al. (2016)</td></tr><tr><td>1705</td><td>Tapio</td><td>25.544</td><td>265</td><td>-48</td><td>106</td><td>-57</td><td>Ďurech et al. (2018a)</td></tr><tr><td>2012</td><td>Guo Shou-Jing</td><td>228.33</td><td>-</td><td>-59 ± 18</td><td>-</td><td>-</td><td>Ďurech et al. (2020)</td></tr><tr><td>2536</td><td>Kozyrev</td><td>7.189</td><td>257</td><td>16</td><td>79</td><td>18</td><td>Ďurech et al. (2020)</td></tr><tr><td>2705</td><td>Wu</td><td>150.8</td><td>356</td><td>-81</td><td>138</td><td>-55</td><td>Ďurech et al. (2020)</td></tr><tr><td>2772</td><td>Dugan</td><td>235.72</td><td>-</td><td>-58 ± 20</td><td>-</td><td>-</td><td>Ďurech et al. (2020)</td></tr><tr><td>2839</td><td>Annette</td><td>10.4609</td><td>154</td><td>-36</td><td>341</td><td>-49</td><td>Hanuš et al. (2013b)</td></tr><tr><td>4231</td><td>Fireman</td><td>339.5</td><td>72</td><td>-43</td><td>258</td><td>-36</td><td>Ďurech et al. (2019)</td></tr><tr><td>4524</td><td>Barklajdetolli</td><td>965.9</td><td>-</td><td>49 ± 16</td><td>-</td><td>-</td><td>Ďurech et al. (2020)</td></tr><tr><td>5081</td><td>Sanguin</td><td>10.26460</td><td>-</td><td>-49 ± 5</td><td>-</td><td>-</td><td>Ďurech et al. (2020)</td></tr><tr><td>5524</td><td>Lecacheux</td><td>8.41706</td><td>-</td><td>-57 ± 9</td><td>-</td><td>-</td><td>Ďurech et al. (2020)</td></tr><tr><td>5924</td><td>Teruo</td><td>9.9918</td><td>340</td><td>-44</td><td>164</td><td>-34</td><td>Ďurech et al. (2019)</td></tr><tr><td>6125</td><td>Singto</td><td>10.2642</td><td>-</td><td>43 ± 17</td><td>-</td><td>-</td><td>Ďurech et al. (2020)</td></tr><tr><td>9723</td><td>Binyang</td><td>12.388</td><td>-</td><td>55 ± 2</td><td>-</td><td>-</td><td>Ďurech et al. (2018a)</td></tr><tr><td>20771</td><td>2000 QY150</td><td>8.3014</td><td>4</td><td>-47</td><td>172</td><td>-48</td><td>Ďurech et al. (2019)</td></tr><tr><td>28736</td><td>2000 GE133</td><td>4.6544</td><td>249</td><td>-52</td><td>134</td><td>-84</td><td>Hanuš et al. (2016)</td></tr><tr><td>30596</td><td>Amdeans</td><td>23.134</td><td>114</td><td>35</td><td>294</td><td>37</td><td>Ďurech et al. (2018a)</td></tr><tr><td>59072</td><td>1998 VV9</td><td>7.2982</td><td>41</td><td>40</td><td>223</td><td>30</td><td>Ďurech et al. (2018a)</td></tr></tbody></table>

**Notes.** *P* is the sidereal period, and *λ* and *β* the ecliptic longitude and latitude of the spin axis, respectively.
<table><caption>**Table B.2.** Photometric observations of targets in the frame of the *Ancient Asteroids* campaign.</caption><thead><tr><th>Number</th><th>Asteroid Name/Designation</th><th>P<sub>synodic</sub> (h)</th><th>Observing season</th><th>N<sub>LC</sub></th><th>N<sub>app</sub></th><th>Site*</th></tr></thead><tbody><tr><td>282</td><td>Clorinde</td><td>49.352(4)<sup>a</sup></td><td>Nov 2020 - Mar 2021</td><td>7+20</td><td>1</td><td>BSA+Lowell</td></tr><tr><td>370</td><td>Modestia</td><td>22.5299(1)<sup>c</sup></td><td>Aug 2021 - Oct 2021</td><td>6+1</td><td>1</td><td>NOAK+OBdB</td></tr><tr><td>428</td><td>Monachia</td><td>3.6343(5)<sup>b</sup></td><td>Mar 2018</td><td>1</td><td>1</td><td>BE600</td></tr><tr><td>853</td><td>Nansenia</td><td>7.931(2)<sup>d</sup></td><td>Dec 2020 - Mar 2021</td><td>15</td><td>1</td><td>Lowell</td></tr><tr><td>917</td><td>Lyka</td><td>7.8838(3)<sup>d</sup></td><td>Oct 2018 - Nov 2018</td><td>8</td><td>1</td><td>OCA(Mont Gros)</td></tr><tr><td>933</td><td>Susi</td><td>4.6222(4)<sup>f</sup></td><td>Feb 2018 - Jan 2021</td><td>3+5+1</td><td>2</td><td>BSA+Lowell+BE600</td></tr><tr><td>1159</td><td>Granada</td><td>77.28(5)<sup>d</sup></td><td>Oct 2021</td><td>1+2</td><td>1</td><td>OCA(C2PU)+UOAO</td></tr><tr><td>1700</td><td>Zvezdara</td><td>9.098(2)<sup>d</sup></td><td>Nov 2020 - Jan 2021</td><td>2+11</td><td>1</td><td>BSA+Lowell</td></tr><tr><td>1806</td><td>Derice</td><td>3.22443(1)<sup>g</sup></td><td>May 2021</td><td>3</td><td>1</td><td>BSA</td></tr><tr><td>1924</td><td>Horus</td><td>6.177(14)<sup>h</sup></td><td>Sep 2021 - Nov 2021</td><td>6</td><td>1</td><td>BO</td></tr><tr><td>2012</td><td>Guo Shou-Jing</td><td>–</td><td>Oct 2021 - Dec 2021</td><td>15</td><td>1</td><td>BO</td></tr><tr><td>2171</td><td>Kiev</td><td>3.1714(2)<sup>i</sup></td><td>Jan 2022</td><td>3</td><td>1</td><td>BO</td></tr><tr><td>2259</td><td>Sofievka</td><td>63.0918(5)<sup>b</sup></td><td>Dec 2020 - Jan 2021</td><td>9</td><td>1</td><td>BSA</td></tr><tr><td>2322</td><td>Kitt-Peak</td><td>8.460(6)<sup>j</sup></td><td>Nov 2020 - Jan 2021</td><td>4+9</td><td>1</td><td>BSA+Lowell</td></tr><tr><td>2575</td><td>Bulgaria</td><td>8.618(7)<sup>k</sup></td><td>Nov 2021 - Dec 2021</td><td>5</td><td>1</td><td>BO</td></tr><tr><td>2768</td><td>Gorky</td><td>4.5118(7)<sup>d</sup></td><td>Dec 2019 - May 2021</td><td>4+3+9</td><td>2</td><td>BSA+BO+UOAO</td></tr><tr><td>2772</td><td>Dugan</td><td>235.0(5)<sup>e</sup></td><td>Dec 2007</td><td>1</td><td>1</td><td>BMO</td></tr><tr><td>2773</td><td>Brooks</td><td>4.838(1)<sup>d</sup></td><td>Dec 2018 - Oct 2021</td><td>3+4</td><td>2</td><td>OCA(Mont Gros) + BSA</td></tr><tr><td>2776</td><td>Baikal</td><td>–</td><td>Jan 2007</td><td>1</td><td>1</td><td>BMO</td></tr><tr><td>2778</td><td>Tangshan</td><td>3.468(3)<sup>f</sup></td><td>Jan 2018 - Oct 2021</td><td>5+3+6</td><td>2</td><td>OCA(Mont Gros)+BSA+Lowell</td></tr><tr><td>2839</td><td>Annette</td><td>10.459(5)<sup>b</sup></td><td>Apr 2020</td><td>5</td><td>1</td><td>UOAO</td></tr><tr><td>3633</td><td>Mira</td><td>19.17(2)<sup>h</sup></td><td>Nov 2021 - Jan 2022</td><td>12</td><td>1</td><td>BO</td></tr><tr><td>3723</td><td>Voznesenskij</td><td>7.9640(85)<sup>h</sup></td><td>Oct 2021</td><td>2+1</td><td>1</td><td>OCA(C2PU)+Helmos</td></tr><tr><td>4231</td><td>Fireman</td><td>28.0(2)<sup>d</sup></td><td>Dec 2020 - Apr 2021</td><td>19</td><td>1</td><td>Lowell</td></tr><tr><td>4422</td><td>Jarre</td><td>7.013(1)<sup>m</sup></td><td>Apr 2021 - May 2021</td><td>5+3+9+3</td><td>1</td><td>BSA+Lowell+BO+UOAO</td></tr><tr><td>5081</td><td>Sanguin</td><td>10.2619(5)<sup>b</sup></td><td>May 2021 - Jun 2021</td><td>6</td><td>1</td><td>BSA</td></tr><tr><td>5333</td><td>Kanaya</td><td>3.69(5)<sup>d</sup></td><td>Jun 2021 - Jul 2021</td><td>3+1</td><td>1</td><td>OCA(C2PU)+ChR</td></tr><tr><td>6125</td><td>Singto</td><td>10.2642(1)<sup>b</sup></td><td>Jan 2022</td><td>3</td><td>1</td><td>BO</td></tr><tr><td>6647</td><td>Josse</td><td>5.9498(3)<sup>d</sup></td><td>Aug 2018</td><td>6</td><td>1</td><td>OCA(Mont Gros)</td></tr><tr><td>7132</td><td>Casulli</td><td>3.5238(2)<sup>n</sup></td><td>Oct 2021</td><td>3+1</td><td>1</td><td>OCA(C2PU)+Helmos</td></tr><tr><td>9086</td><td>1995 SA3</td><td>–</td><td>Oct 2019 - Nov 2019</td><td>22</td><td>1</td><td>UOAO</td></tr><tr><td>9972</td><td>Minoruoda</td><td>3.4221(2)<sup>o</sup></td><td>Oct 2021</td><td>1</td><td>1</td><td>Helmos</td></tr><tr><td>10542</td><td>Ruckers</td><td>–</td><td>Feb 2022 - Mar 2022</td><td>4</td><td>1</td><td>BO</td></tr><tr><td>15985</td><td>1998 WU20</td><td>18.2252(5)<sup>b</sup></td><td>Oct 2021</td><td>4+1</td><td>1</td><td>C2PU+Helmos</td></tr><tr><td>25343</td><td>1999 RA44</td><td>590.5(5)<sup>m</sup></td><td>Jan 2021 - Apr 2021</td><td>27</td><td>1</td><td>BO</td></tr></tbody></table>

**Notes.** The literature value of the (synodic) rotational period is given, along with the observing log, where $N_{LC}$ is the number of individual light curves obtained by each corresponding observing site given in the last column, while $N_{app}$ is the number of apparitions. The standard error for each value is expressed in brackets, in units of the last decimal digit quoted.

(*)BSA: Astronomical Observatory BSA, Lowell: Lowell Observatory, BO: Bigmuskie Observatory, OCA: Observatoire de la Côte d'Azur - Calern station (C2PU) and Mont Gros station, UOAO: University of Athens Observatory, Helmos: Helmos Observatory of National Observatory of Athens, NOAK: NOAK Observatory, BE600: BlueEye 600 Observatory, ChR: Pic de Château-Renard Observatory, OBdB: Observatoire du Bois de Bardon, BMO:Blue Mountains Observatory

(a)Bonamico & van Belle (2021), (b)Pál et al. (2020), (c)Stephens (2011), (d)http://obswww.unige.ch/~behrend/page_cou.html, (e)https://www.asu.cas.cz/~ppravec/, (f)https://web.archive.org/web/20081004205615/http://www.david-higgins.com/, (g)Stephens & Warner (2020), (h)Waszczak et al. (2015), (i)Loera-González et al. (2020), (j)Polakis (2021), (k)Erasmus et al. (2020), (l)Stephens & Warner (2019), (m)Ferrero (2021), (n)Franco et al. (2020b), (o)Cooney et al. (2017)
<table><caption>**Table B.3.** Rotation state properties and summary of the optical dataset for asteroids for which we derived a new shape solution or estimated the sense of rotation.</caption><thead><tr><th>Number</th><th>Name/Designation</th><th>λ<sub>1</sub><br>(deg)</th><th>β<sub>1</sub><br>(deg)</th><th>λ<sub>2</sub><br>(deg)</th><th>β<sub>2</sub><br>(deg)</th><th>P<br>(h)</th><th>N<sub>LC</sub></th><th>N<sub>app</sub></th><th>N<sub>689</sub></th><th>N<sub>703</sub></th><th>N<sub>GAI</sub></th><th>N<sub>ASA</sub></th><th>N<sub>ATL</sub></th><th>N<sub>I41</sub></th><th>N<sub>PTF</sub></th></tr></thead><tbody><tr><td colspan="16" style="text-align:center;">New asteroid models</td></tr><tr><td>917</td><td>Lyka</td><td>44</td><td>49</td><td>230</td><td>37</td><td>7.88167</td><td>13</td><td>1</td><td>116</td><td>205</td><td>9</td><td>599</td><td>582</td><td>26</td><td>0</td></tr><tr><td>1159</td><td>Granada</td><td>70</td><td>-54</td><td>142</td><td>-85</td><td>72.789</td><td>5</td><td>2</td><td>130</td><td>196</td><td>14</td><td>382</td><td>508</td><td>77</td><td>49</td></tr><tr><td>1544</td><td>Vinterhansenia</td><td>55</td><td>22</td><td>238</td><td>23</td><td>13.78270</td><td>0</td><td>0</td><td>104</td><td>215</td><td>9</td><td>389</td><td>341</td><td>0</td><td>0</td></tr><tr><td>1700</td><td>Zvezdara</td><td>66</td><td>62</td><td>248</td><td>48</td><td>9.11217</td><td>43</td><td>2</td><td>116</td><td>176</td><td>20</td><td>354</td><td>324</td><td>51</td><td>0</td></tr><tr><td>1806</td><td>Derice</td><td>42</td><td>43</td><td>216</td><td>38</td><td>3.223531</td><td>17</td><td>2</td><td>59</td><td>286</td><td>0</td><td>406</td><td>488</td><td>20</td><td>0</td></tr><tr><td>1924</td><td>Horus</td><td>131</td><td>-59</td><td>315</td><td>-54</td><td>6.17893</td><td>6</td><td>1</td><td>0</td><td>178</td><td>0</td><td>163</td><td>459</td><td>29</td><td>21</td></tr><tr><td>2171</td><td>Kiev</td><td>145</td><td>54</td><td></td><td></td><td>3.171587</td><td>7</td><td>2</td><td>0</td><td>277</td><td>12</td><td>339</td><td>451</td><td>0</td><td>0</td></tr><tr><td>2322</td><td>KittPeak</td><td>195</td><td>-78</td><td></td><td></td><td>8.46780</td><td>13</td><td>1</td><td>0</td><td>191</td><td>0</td><td>319</td><td>425</td><td>0</td><td>0</td></tr><tr><td>2575</td><td>Bulgaria</td><td>66</td><td>91</td><td>246</td><td>70</td><td>8.61720</td><td>5</td><td>1</td><td>20</td><td>205</td><td>0</td><td>309</td><td>445</td><td>0</td><td>28</td></tr><tr><td>2768</td><td>Gorky</td><td>54</td><td>34</td><td></td><td></td><td>4.50814</td><td>18</td><td>2</td><td>37</td><td>217</td><td>11</td><td>151</td><td>545</td><td>0</td><td>0</td></tr><tr><td>2773</td><td>Brooks</td><td>85</td><td>-86</td><td>241</td><td>-71</td><td>4.83726</td><td>7</td><td>2</td><td>0</td><td>162</td><td>14</td><td>110</td><td>328</td><td>0</td><td>0</td></tr><tr><td>2776</td><td>Baikal</td><td>81</td><td>-41</td><td>264</td><td>-27</td><td>3253.5</td><td>1</td><td>1</td><td>0</td><td>143</td><td>0</td><td>312</td><td>479</td><td>44</td><td>0</td></tr><tr><td>2792</td><td>Ponomarev</td><td>72</td><td>-58</td><td>230</td><td>-82</td><td>138.28</td><td>0</td><td>0</td><td>0</td><td>193</td><td>0</td><td>223</td><td>409</td><td>92</td><td>60</td></tr><tr><td>3633</td><td>Mira</td><td>170</td><td>-58</td><td>346</td><td>-50</td><td>19.1832</td><td>12</td><td>1</td><td>0</td><td>275</td><td>10</td><td>124</td><td>331</td><td>0</td><td>29</td></tr><tr><td>3684</td><td>Berry</td><td>96</td><td>-57</td><td>255</td><td>-49</td><td>11.9074</td><td>0</td><td>0</td><td>0</td><td>186</td><td>0</td><td>0</td><td>332</td><td>0</td><td>0</td></tr><tr><td>3723</td><td>Voznesenskij</td><td>75</td><td>-35</td><td>259</td><td>-34</td><td>7.96485</td><td>3</td><td>1</td><td>0</td><td>185</td><td>0</td><td>52</td><td>292</td><td>0</td><td>20</td></tr><tr><td>5333</td><td>Kanaya</td><td>157</td><td>-56</td><td></td><td></td><td>3.80247</td><td>9</td><td>2</td><td>0</td><td>172</td><td>25</td><td>210</td><td>533</td><td>0</td><td>0</td></tr><tr><td>6647</td><td>Josse</td><td>38</td><td>-68</td><td>222</td><td>-62</td><td>5.94972</td><td>6</td><td>1</td><td>0</td><td>165</td><td>0</td><td>24</td><td>271</td><td>0</td><td>0</td></tr><tr><td>8022</td><td>Scotcrossfield</td><td>186</td><td>30</td><td>359</td><td>40</td><td>4.00267</td><td>0</td><td>0</td><td>0</td><td>228</td><td>18</td><td>84</td><td>293</td><td>0</td><td>76</td></tr><tr><td>8315</td><td>Bajin</td><td>6</td><td>71</td><td></td><td></td><td>152.62</td><td>0</td><td>0</td><td>0</td><td>167</td><td>0</td><td>0</td><td>262</td><td>74</td><td>0</td></tr><tr><td>12722</td><td>Petrarca</td><td>143</td><td>24</td><td>346</td><td>41</td><td>3.538576</td><td>0</td><td>0</td><td>0</td><td>138</td><td>0</td><td>0</td><td>175</td><td>41</td><td>0</td></tr><tr><td>15415</td><td>Rika</td><td>117</td><td>-79</td><td></td><td></td><td>6.36228</td><td>0</td><td>0</td><td>0</td><td>164</td><td>0</td><td>0</td><td>275</td><td>0</td><td>0</td></tr><tr><td>23495</td><td>1991 UQ1</td><td>84</td><td>-70</td><td></td><td></td><td>10.50559</td><td>0</td><td>0</td><td>0</td><td>124</td><td>0</td><td>26</td><td>370</td><td>0</td><td>0</td></tr><tr><td>49863</td><td>1999 XK104</td><td>32</td><td>-18</td><td>212</td><td>-22</td><td>8.00795</td><td>0</td><td>0</td><td>0</td><td>72</td><td>0</td><td>0</td><td>94</td><td>0</td><td>0</td></tr><tr><th>Number</th><th>Name/Designation</th><th>λ<sub>1</sub><br>(deg)</th><th>β<sub>1</sub><br>(deg)</th><th>λ<sub>2</sub><br>(deg)</th><th>β<sub>2</sub><br>(deg)</th><th>P<br>(h)</th><th>N<sub>LC</sub></th><th>N<sub>app</sub></th><th>N<sub>689</sub></th><th>N<sub>703</sub></th><th>N<sub>GAI</sub></th><th>N<sub>ASA</sub></th><th>N<sub>ATL</sub></th><th>N<sub>I41</sub></th><th>N<sub>PTF</sub></th></tr><tr><td colspan="16" style="text-align:center;">Partial asteroid models</td></tr><tr><td>2778</td><td>Tangshan</td><td>-65</td><td>18</td><td></td><td></td><td>3.460693</td><td>19</td><td>4</td><td>0</td><td>169</td><td>0</td><td>130</td><td>263</td><td>0</td><td>0</td></tr><tr><td>11975</td><td>1995 FA1</td><td>65</td><td>9</td><td></td><td></td><td>6.21338</td><td>0</td><td>0</td><td>0</td><td>119</td><td>0</td><td>0</td><td>235</td><td>0</td><td>0</td></tr><tr><td>13066</td><td>1991 PM13</td><td>-57</td><td>19</td><td></td><td></td><td>3.93765</td><td>0</td><td>0</td><td>0</td><td>160</td><td>0</td><td>90</td><td>259</td><td>33</td><td>0</td></tr><tr><td>70184</td><td>1999 RU3</td><td>-59</td><td>14</td><td></td><td></td><td>5.27716</td><td>0</td><td>0</td><td>0</td><td>150</td><td>0</td><td>46</td><td>320</td><td>0</td><td>0</td></tr></tbody></table>

**Notes.** The new determinations are listed first, then the cases where we estimated only the ecliptic latitude of the spin axis (i.e. partial models). Listed are the physical properties: ecliptic longitudes and latitudes of the spin axis directions λ and β for one or two possible solutions; the sidereal rotation period *P*; the mean value of the ecliptic latitude *β*<sub>p</sub> and 1/2 of the range in latitude within the multiple pole solutions *δ*β (for partial models). The uncertainty of the pole direction is usually about 10° and of the period is of the order of the last decimal digit. Also given is information about the light curve dataset: number of dense light curves N<sub>LC</sub> with the number of covered apparitions N<sub>app</sub>; the number of measurements in each sparse dataset (N<sub>689</sub>: UNSO-Flagstaff; N<sub>703</sub>: CSS; N<sub>GAI</sub>: Gaia DR2; N<sub>ASA</sub>: ASAS-SN; N<sub>ATL</sub>: ATLAS; N<sub>I41</sub>: ZTF; N<sub>PTF</sub>: PTF).
<table><caption>Table B.4. Updated shape models.</caption><thead><tr><th>Number</th><th>Name/Designation</th><th>λ<sub>1</sub><br>(deg)</th><th>β<sub>1</sub><br>(deg)</th><th>λ<sub>2</sub><br>(deg)</th><th>β<sub>2</sub><br>(deg)</th><th>P<br>(h)</th><th>N<sub>LC</sub></th><th>N<sub>app</sub></th><th>N<sub>689</sub></th><th>N<sub>703</sub></th><th>N<sub>GAI</sub></th><th>N<sub>ASA</sub></th><th>N<sub>ATL</sub></th><th>N<sub>I41</sub></th><th>N<sub>PTF</sub></th></tr></thead><tbody><tr><td>220</td><td>Stephania</td><td>36</td><td>-51</td><td></td><td></td><td>18.2087</td><td>9</td><td>2</td><td>137</td><td>242</td><td>12</td><td>376</td><td>599</td><td>30</td><td>0</td></tr><tr><td>249</td><td>Ilse</td><td>85</td><td>79</td><td>266</td><td>69</td><td>84.996</td><td>29</td><td>3</td><td>135</td><td>254</td><td>9</td><td>499</td><td>529</td><td>0</td><td>0</td></tr><tr><td>282</td><td>Clorinde</td><td>189</td><td>-36</td><td>346</td><td>-59</td><td>49.3597</td><td>29</td><td>2</td><td>175</td><td>197</td><td>13</td><td>659</td><td>548</td><td>0</td><td>0</td></tr><tr><td>370</td><td>Modestia</td><td>268</td><td>-92</td><td></td><td></td><td>22.5403</td><td>21</td><td>3</td><td>176</td><td>202</td><td>11</td><td>564</td><td>472</td><td>0</td><td>0</td></tr><tr><td>428</td><td>Monachia</td><td>98</td><td>90</td><td>283</td><td>53</td><td>3.633613</td><td>28</td><td>16</td><td>137</td><td>260</td><td>25</td><td>519</td><td>396</td><td>0</td><td>0</td></tr><tr><td>933</td><td>Susi</td><td>300</td><td>-9</td><td></td><td></td><td>4.62241</td><td>18</td><td>3</td><td>114</td><td>149</td><td>0</td><td>385</td><td>439</td><td>33</td><td>0</td></tr><tr><td>1216</td><td>Askania</td><td>62</td><td>55</td><td>260</td><td>66</td><td>6.53708</td><td>0</td><td>0</td><td>55</td><td>138</td><td>0</td><td>191</td><td>370</td><td>0</td><td>0</td></tr><tr><td>2012</td><td>Guo Shou-Jing</td><td>46</td><td>-79</td><td>230</td><td>-76</td><td>228.29</td><td>15</td><td>1</td><td>0</td><td>200</td><td>12</td><td>83</td><td>422</td><td>0</td><td>0</td></tr><tr><td>2536</td><td>Kozyrev</td><td>80</td><td>22</td><td>256</td><td>21</td><td>7.18887</td><td>3</td><td>1</td><td>27</td><td>138</td><td>14</td><td>348</td><td>494</td><td>0</td><td>0</td></tr><tr><td>2705</td><td>Wu</td><td>86</td><td>-70</td><td>265</td><td>-82</td><td>150.78</td><td>0</td><td>0</td><td>0</td><td>174</td><td>0</td><td>191</td><td>429</td><td>55</td><td>0</td></tr><tr><td>2772</td><td>Dugan</td><td>157</td><td>-71</td><td></td><td></td><td>235.76</td><td>1</td><td>1</td><td>0</td><td>174</td><td>0</td><td>120</td><td>360</td><td>0</td><td>0</td></tr><tr><td>2839</td><td>Annette</td><td>338</td><td>-47</td><td></td><td></td><td>10.46100</td><td>13</td><td>17</td><td>0</td><td>242</td><td>13</td><td>319</td><td>557</td><td>0</td><td>0</td></tr><tr><td>4231</td><td>Fireman</td><td>250</td><td>-83</td><td></td><td></td><td>339.57</td><td>16</td><td>1</td><td>26</td><td>228</td><td>13</td><td>248</td><td>511</td><td>88</td><td>0</td></tr><tr><td>4524</td><td>Barklajdetolli</td><td>76</td><td>70</td><td>257</td><td>93</td><td>966.5</td><td>19</td><td>2</td><td>0</td><td>175</td><td>0</td><td>188</td><td>391</td><td>0</td><td>0</td></tr><tr><td>5081</td><td>Sanguin</td><td>242</td><td>-39</td><td></td><td></td><td>10.26457</td><td>11</td><td>20</td><td>41</td><td>250</td><td>0</td><td>248</td><td>627</td><td>0</td><td>0</td></tr><tr><td>5524</td><td>Lecacheux</td><td>254</td><td>-86</td><td></td><td></td><td>8.41707</td><td>5</td><td>1</td><td>0</td><td>196</td><td>9</td><td>89</td><td>468</td><td>31</td><td>0</td></tr><tr><td>5924</td><td>Teruo</td><td>184</td><td>-60</td><td>359</td><td>-71</td><td>9.99174</td><td>6</td><td>1</td><td>20</td><td>153</td><td>13</td><td>217</td><td>447</td><td>0</td><td>0</td></tr><tr><td>9723</td><td>Binyang</td><td>36</td><td>45</td><td>215</td><td>38</td><td>12.38806</td><td>2</td><td>1</td><td>0</td><td>127</td><td>0</td><td>0</td><td>357</td><td>0</td><td>26</td></tr><tr><td>20771</td><td>2000 QY150</td><td>159</td><td>-63</td><td></td><td></td><td>8.30140</td><td>0</td><td>0</td><td>0</td><td>197</td><td>14</td><td>26</td><td>386</td><td>0</td><td>0</td></tr><tr><td>28736</td><td>2000 GE133</td><td>241</td><td>-65</td><td></td><td></td><td>4.65441</td><td>3</td><td>1</td><td>0</td><td>237</td><td>0</td><td>53</td><td>493</td><td>0</td><td>0</td></tr></tbody></table>

Notes. Same as in Table B.3.
<table><caption>**Table B.6.** Physical properties of asteroids of the inner main belt primordial family presented in this study.</caption><thead><tr><th>Number</th><th>Asteroid Name/Designation</th><th>D (km)</th><th>σD (km)</th><th>p<sub>V</sub></th><th>σp<sub>V</sub></th><th>Spectral Class</th><th>Ref.</th></tr></thead><tbody><tr><td colspan="8" style="text-align:center; font-weight:bold;">Prograde rotators</td></tr><tr><td>249</td><td>Ilse</td><td>31.57</td><td>0.300</td><td>0.054</td><td>0.0014</td><td>Ch</td><td>Lazzaro et al. (2004)</td></tr><tr><td>428</td><td>Monachia</td><td>20.55</td><td>0.129</td><td>0.066</td><td>0.0055</td><td>X</td><td>Alvarez-Candal et al. (2006)</td></tr><tr><td>917</td><td>Lyka</td><td>35.61</td><td>0.135</td><td>0.043</td><td>0.0040</td><td>X</td><td>Lazzaro et al. (2004)</td></tr><tr><td>1216</td><td>Askania</td><td>10.35</td><td>0.086</td><td>0.086</td><td>0.0075</td><td>-</td><td></td></tr><tr><td>1544</td><td>Vinterhansenia</td><td>24.78</td><td>0.074</td><td>0.049</td><td>0.0036</td><td>X,D</td><td>Carvano et al. (2010); Alvarez-Candal et al. (2006)</td></tr><tr><td>1700</td><td>Zvezdara</td><td>20.54</td><td>0.181</td><td>0.039</td><td>0.0016</td><td>X</td><td>Zellner et al. (1985)</td></tr><tr><td>1806</td><td>Derice</td><td>10.67</td><td>0.060</td><td>0.219</td><td>0.0512</td><td>S1</td><td>Lazzaro et al. (2004)</td></tr><tr><td>2171</td><td>Kiev</td><td>8.30</td><td>0.055</td><td>0.101</td><td>0.0059</td><td>S</td><td>Avdellidou et al. (in prep.)</td></tr><tr><td>2536</td><td>Kozyrev</td><td>9.59</td><td>0.218</td><td>0.195</td><td>0.0227</td><td>-</td><td></td></tr><tr><td>2575</td><td>Bulgaria</td><td>7.92</td><td>0.063</td><td>0.270</td><td>0.0300</td><td>Sr,S</td><td>Bus & Binzel (2002); Popescu et al. (2018)</td></tr><tr><td>2768</td><td>Gorky</td><td>10.67</td><td>0.085</td><td>0.258</td><td>0.0367</td><td>A</td><td>Alvarez-Candal et al. (2006)</td></tr><tr><td>4024</td><td>Ronan</td><td>11.90</td><td>0.073</td><td>0.055</td><td>0.0032</td><td>-</td><td></td></tr><tr><td>4524</td><td>Barklajdetolli</td><td>12.62</td><td>0.146</td><td>0.087</td><td>0.0069</td><td>U</td><td>Popescu et al. (2018)</td></tr><tr><td>6125</td><td>Singto</td><td>6.30</td><td>0.479</td><td>0.109</td><td>0.0382</td><td>S</td><td>Popescu et al. (2018)</td></tr><tr><td>8022</td><td>Scottcrossfield</td><td>8.51</td><td>0.147</td><td>0.046</td><td>0.0067</td><td>C</td><td>Carvano et al. (2010)</td></tr><tr><td>8315</td><td>Bajin</td><td>7.44</td><td>0.665</td><td>0.039</td><td>0.0086</td><td>-</td><td></td></tr><tr><td>9723</td><td>Binyang</td><td>3.69</td><td>0.095</td><td>0.117</td><td>0.0320</td><td>-</td><td></td></tr><tr><td>11975</td><td>1995 FA1</td><td>4.360</td><td>1.440</td><td>0.060</td><td>0.0500</td><td>-</td><td></td></tr><tr><td>12722</td><td>Petrarca</td><td>4.52</td><td>0.746</td><td>0.095</td><td>0.0384</td><td>C,U</td><td>Carvano et al. (2010); Popescu et al. (2018)</td></tr><tr><td>30596</td><td>Amdeans</td><td>5.11</td><td>0.158</td><td>0.056</td><td>0.0074</td><td>-</td><td></td></tr><tr><td>59072</td><td>1998 VV9</td><td>4.59</td><td>0.070</td><td>0.038</td><td>0.0045</td><td>-</td><td></td></tr><tr><td colspan="8" style="text-align:center; font-weight:bold;">Retrograde rotators</td></tr><tr><td>220</td><td>Stephania</td><td>32.54</td><td>0.138</td><td>0.060</td><td>0.0032</td><td>X,Xk</td><td>Carvano et al. (2010); Lazzaro et al. (2004)</td></tr><tr><td>282</td><td>Clorinde</td><td>34.01</td><td>0.148</td><td>0.045</td><td>0.0019</td><td>B,C</td><td>Bus & Binzel (2002); Popescu et al. (2018)</td></tr><tr><td>370</td><td>Modestia</td><td>38.11</td><td>0.106</td><td>0.059</td><td>0.0046</td><td>-</td><td></td></tr><tr><td>933</td><td>Susi</td><td>23.42</td><td>0.245</td><td>0.042</td><td>0.0032</td><td>C</td><td>Carvano et al. (2010); Popescu et al. (2018)</td></tr><tr><td>1159</td><td>Granada</td><td>30.14</td><td>0.099</td><td>0.046</td><td>0.0019</td><td>-</td><td></td></tr><tr><td>1244</td><td>Deira</td><td>33.97</td><td>0.152</td><td>0.046</td><td>0.0029</td><td>X</td><td>Lazzaro et al. (2004)</td></tr><tr><td>1705</td><td>Tapio</td><td>11.81</td><td>0.061</td><td>0.093</td><td>0.0068</td><td>B,U</td><td>Bus & Binzel (2002); Popescu et al. (2018)</td></tr><tr><td>1924</td><td>Horus</td><td>12.90</td><td>0.130</td><td>0.070</td><td>0.0036</td><td>-</td><td></td></tr><tr><td>2012</td><td>Guo Shou-Jing</td><td>11.91</td><td>0.076</td><td>0.048</td><td>0.0016</td><td>C</td><td>Carvano et al. (2010); Popescu et al. (2018)</td></tr><tr><td>2322</td><td>Kitt-Peak</td><td>11.91</td><td>0.085</td><td>0.058</td><td>0.0060</td><td>-</td><td></td></tr><tr><td>2705</td><td>Wu</td><td>7.79</td><td>0.319</td><td>0.163</td><td>0.0330</td><td>-</td><td></td></tr><tr><td>2772</td><td>Dugan</td><td>9.58</td><td>0.135</td><td>0.057</td><td>0.0075</td><td>B</td><td>Bus & Binzel (2002)</td></tr><tr><td>2773</td><td>Brooks</td><td>13.37</td><td>0.100</td><td>0.042</td><td>0.0031</td><td>-</td><td></td></tr><tr><td>2776</td><td>Baikal</td><td>13.37</td><td>0.100</td><td>0.042</td><td>0.0031</td><td>-</td><td></td></tr><tr><td>2778</td><td>Tangshan</td><td>12.66</td><td>0.144</td><td>0.062</td><td>0.0080</td><td>Cb</td><td>Bus & Binzel (2002)</td></tr><tr><td>2792</td><td>Ponomarev</td><td>12.52</td><td>0.222</td><td>0.056</td><td>0.0112</td><td>-</td><td></td></tr><tr><td>2839</td><td>Annette</td><td>7.44</td><td>0.094</td><td>0.059</td><td>0.0076</td><td>-</td><td></td></tr><tr><td>3633</td><td>Mira</td><td>10.09</td><td>0.258</td><td>0.045</td><td>0.0058</td><td>X</td><td>Carvano et al. (2010)</td></tr><tr><td>3684</td><td>Berry</td><td>9.64</td><td>0.495</td><td>0.053</td><td>0.0096</td><td>C</td><td>Carvano et al. (2010); Bus & Binzel (2002)</td></tr><tr><td>3723</td><td>Voznesenskij</td><td>9.54</td><td>0.029</td><td>0.041</td><td>0.0015</td><td>C</td><td>Carvano et al. (2010)</td></tr><tr><td>4231</td><td>Fireman</td><td>17.96</td><td>0.557</td><td>0.022</td><td>0.0014</td><td>-</td><td></td></tr><tr><td>5081</td><td>Sanguin</td><td>17.19</td><td>0.263</td><td>0.056</td><td>0.0063</td><td>Ch</td><td>Bus & Binzel (2002)</td></tr><tr><td>5333</td><td>Kanaya</td><td>13.70</td><td>0.035</td><td>0.040</td><td>0.0011</td><td>Ch</td><td>Bus & Binzel (2002); Morate et al. (2019)</td></tr><tr><td>5524</td><td>Lecacheux</td><td>19.90</td><td>12.770</td><td>0.034</td><td>0.1020</td><td>V</td><td>Carvano et al. (2010); Popescu et al. (2018)</td></tr><tr><td>5924</td><td>Teruo</td><td>13.16</td><td>0.080</td><td>0.059</td><td>0.0015</td><td>-</td><td></td></tr><tr><td>6647</td><td>Josse</td><td>6.42</td><td>0.148</td><td>0.049</td><td>0.0069</td><td>C</td><td>Carvano et al. (2010)</td></tr><tr><td>13066</td><td>1991 PM13</td><td>8.60</td><td>0.069</td><td>0.038</td><td>0.0058</td><td>B</td><td>Popescu et al. (2018)</td></tr><tr><td>15415</td><td>Rika</td><td>2.83</td><td>0.488</td><td>0.605</td><td>0.1924</td><td>V</td><td>Popescu et al. (2018)</td></tr><tr><td>15998</td><td>1999 AG2</td><td>7.13</td><td>0.046</td><td>0.087</td><td>0.0062</td><td>-</td><td></td></tr><tr><td>20771</td><td>2000 QY150</td><td>9.08</td><td>0.115</td><td>0.047</td><td>0.0068</td><td>-</td><td></td></tr><tr><td>23495</td><td>1991 UQ1</td><td>7.94</td><td>0.088</td><td>0.058</td><td>0.0054</td><td>-</td><td></td></tr><tr><td>28736</td><td>2000 GE133</td><td>7.03</td><td>0.680</td><td>0.084</td><td>0.0186</td><td>C</td><td>Carvano et al. (2010)</td></tr><tr><td>49863</td><td>1999 XK104</td><td>3.93</td><td>0.189</td><td>0.054</td><td>0.0053</td><td>-</td><td></td></tr><tr><td>70184</td><td>1999 RU3</td><td>4.31</td><td>0.629</td><td>0.137</td><td>0.0557</td><td>-</td><td></td></tr></tbody></table>

Notes. *D* and *p<sub>V</sub>* are the diameter and the geometric visible albedo, respectively. These values are from the Minor Planet Physical Properties Catalogue and represent the uncertainty-weighted average values of each asteroid. Where the spectral class is known, it is given with its reference.