# On the similarity of rotational motion of dust particles in the inner atmosphere of comets

Stavro L. Ivanovski, V. V. Zakharov, F. Moreno, N. Y. Bykov, O. Muñoz, M.

Fulle, A. Rotundi, V. Della Corte, A. V. Rodionov

► To cite this version:

Stavro L. Ivanovski, V. V. Zakharov, F. Moreno, N. Y. Bykov, O. Muñoz, et al. On the similarity of rotational motion of dust particles in the inner atmosphere of comets. Monthly Notices of the Royal Astronomical Society, 2024, 528, pp.5723-5729. 10.1093/mnras/stad2665. insu-04853464

HAL Id: insu-04853464

https://insu.hal.science/insu-04853464v1

Submitted on 23 Dec 2024

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>

Distributed under a Creative Commons Attribution 4.0 International License
# On the similarity of rotational motion of dust particles in the inner atmosphere of comets

Stavro L. Ivanovski,<fnref n="1" />★ V. V. Zakharov,<fnref n="2" /> F. Moreno,<fnref n="3" /> N. Y. Bykov,<fnref n="4" /> O. Muñoz,<fnref n="3" /> M. Fulle,<fnref n="1" /> A. Rotundi,<fnref n="5" />,<fnref n="6" /> V. Della Corte<fnref n="5" />,<fnref n="7" /> and A. V. Rodionov<fnref n="8" />

Accepted 2023 September 1. Received 2023 September 1; in original form 2023 June 3

## ABSTRACT

In situ measurements of individual dust grain parameters in the immediate vicinity of a cometary nucleus are being carried by the Rosetta spacecraft at comet 67P/Churyumov–Gerasimenko. For interpretation of these observational data, a model of dust grain motion as realistic as possible is requested. In particular, the results of Stardust mission and analysis of samples of interplanetary dust have shown that these particles are highly non-spherical. In many cases precise simulations of non-spherical grain’s dynamics is either impossible or computationally too expensive. In such situation it is proposed to use available experimental or numerical data obtained for other conditions and rescale them considering similarity of the physical process. In the present paper we focus on the derivation of scaling laws of rotational motion applicable for any shape of particles. We use a set of universal, dimensionless parameters characterizing the dust motion in the inner cometary coma. The scaling relations for translational and rotational motion of dust grains in a cometary environment are proposed. The scaled values are compared with numerically computed ones in our previous works.

**Key words:** methods: numerical – comets: general.

## 1 INTRODUCTION

The Rosetta probe accumulated plenty of novel scientific data from the closest vicinity of the comet 67P/Churyumov–Gerasimenko (67P). In October 2014, at 3.0 au heliocentric distance and in bound orbit around the comet of about 10 km, the Optical, Spectroscopic, and Infrared Remote Imaging System camera (OSIRIS; Keller et al. (2007)) onboard Rosetta gained images showing a coma that was completely resolved in single particles of subpixel apparent sizes. The vectorial sum of the dust proper motion and spacecraft motion produced long rectilinear tracks in the images. Many of these tracks showed regular brightness variations, probably due to the rotation of particles during the exposure.

The first attempt to interpret the dust tracks of 67P with brightness periodicity by a model of dust dynamics assuming simple non-spherical grains (ellipsoids) was made in Fulle et al. (2015). In order to constrain the rotational motion it was necessary to perform a large number of numerical simulations for a vast number of physical and dynamical parameters of the dust particles and the comet, for

example, comet gas production rate, particle mass, size, orientation, initial velocity, etc.

A detailed analysis of the dynamics of non-spherical dust particles in the vicinity of a cometary nucleus was performed in Ivanovski (2017) (spheroidal shape particles) and Moreno et al. (2022) (irregularly shaped particles). For some detected 67P particles with effective radii less than 10 µm, the rotation frequencies can exceed more than 500 Hz, which in return leads to a very small time-step and makes the numerical simulation not feasible.

In many cases the precise simulation of non-spherical grain’s dynamics is either impossible or computationally too expensive. In such situation, it is proposed to use available experimental or numerical data for other conditions and rescale them considering the similarity of the physical process.

In the present paper, we focus on the derivation of scaling laws of rotational motion applicable for any shape of particles. We use a set of universal, dimensionless parameters characterizing the dust motion in the inner cometary coma. This approach was proposed in Zakharov et al. (2018) and successfully applied in Zakharov et al. (2021) to the dynamics of spherical grains. Here, we extend this approach on the translational and rotational motion of non-spherical particles. Our approach allows to reduce the number of parameters

<footnote n="1">¹INAF – Osservatorio Astronomico di Trieste, Via Tiepolo 11, I-34143 Trieste, Italy</footnote>

<footnote n="2">²LESIA, Observatoire de Paris, Universite PSL, Sorbonne Universite, Universite de Paris, 5 place Jules Janssen, F-92195 Meudon, France</footnote>

<footnote n="3">³Instituto de Astrofísica de Andalucía, CSIC, Glorieta de la Astronomía s/n, E-18008 Granada, Spain</footnote>

<footnote n="4">⁴Peter the Great St. Petersburg Polytechnic University, Polytechnicheskaya street, 29, 195251 St. Petersburg, Russia</footnote>

<footnote n="5">⁵Università degli Studi di Napoli Parthenope, Dip. di Scienze e Tecnologie, CDN IC4, I-80143 Naples, Italy</footnote>

<footnote n="6">⁶INAF – Istituto di Astrofisica e Planetologia Spaziali, Area Ricerca Tor Vergata, Via Fosso del Cavaliere 100, I-00133 Rome, Italy</footnote>

<footnote n="7">⁷INAF – Osservatorio Astronomico di Capodimonte, Salita Moariello, 16, I-80131 Napoli, Italy</footnote>

<footnote n="8">⁸RFNC – All-Russian Scientific Research Institute of Experimental Physics, prospect Mira 37, 607188 Sarov, Russia</footnote>

<footnote>* E-mail: stavro.ivanovski@inaf.it</footnote>
for analysis, to reveal similarities of the dust flows and to rescale the available numerical solutions.

## 2 THE MODEL

We assume that the dusty-gas coma is formed by the gas sublimating from the nucleus (from the surface and/or from the interior) and solid particles (mineral or/and icy) released from the nucleus with zero initial velocity and entrained by the gas flow. It is assumed that the dusty-gas flow is coupled in one way only – the gas drags the dust (i.e. the presence of dust in the coma does not affect the gas motion), and that the dust particles do not collide with each other. The dust particles are assumed to be isothermal, internally homogeneous with invariable size and mass.

In our analysis of the dust dynamics we use two frames (i) a fixed cometocentric coordinate system XYZ with basis $i_1 i_2 i_3$; (ii) a moving frame attached to the grain with an origin in its centre of inertia with axes along the principal moments of inertia xyz with basis $e_1 e_2 e_3$. Orientation of the grain in XYZ is given by three Euler's angles: the angle of precession $\psi$, the angle of nutation $\theta$ and the angle of self-rotation $\varphi$ (see Fig. 1).

The motion of the centre of mass of a dust grain is given by

$$m_d \frac{d^2 \mathbf{r}}{dt^2} = m_d \frac{d\mathbf{v}_d}{dt} = \mathbf{R} \equiv \mathbf{F}_N + \mathbf{F}_a, \quad (1)$$

where $\mathbf{r}$ and $\mathbf{v}_d$ are the radius vector and the velocity of the centre of mass in the fixed cometocentric frame, $m_d$ is the grain's mass, and $\mathbf{R}$ is a resultant vector of external forces, that is, the sum of the gravity $\mathbf{F}_N$ and the aerodynamic force $\mathbf{F}_a$.

From the law of variation of the angular momentum, we can get in the moving frame the Euler dynamic equations

$$\mathbf{I} \frac{d\omega}{dt} + \omega \times \mathbf{I} \omega = \mathbf{M}, \quad (2)$$

and the Euler kinematic equation

$$\omega = \frac{d\psi}{dt} \mathbf{i}_3 + \frac{d\theta}{dt} \mathbf{i}'_1 + \frac{d\varphi}{dt} \mathbf{e}_3, \quad (3)$$

where $\mathbf{M} = \mathbf{M}_N + \mathbf{M}_a$ is a resultant torque of external forces, that is, the sum of the gravitational $\mathbf{M}_N$ and aerodynamic $\mathbf{M}_a$ torques, $\omega$ is a vector of instantaneous angular velocity, $\mathbf{I}$ is a tensor of inertia, and $\mathbf{i}_3, \mathbf{i}'_1, \mathbf{e}_3$ are the unit vectors (see Fig. 1).

The equations (1), (2), and (3) form a close system of differential equations that describes motion of a grain.

<figure><img src="image_2.png" /><figcaption>Figure 1. Euler's angles.</figcaption></figure>

The gravitational and aerodynamic force $\vec{F}_N$ and $\vec{F}_a$ are given by

$$\mathbf{F}_N = -G \frac{m_N m_d}{r^3} \mathbf{r} \quad (4)$$

$$\mathbf{F}_a = -\int (p\mathbf{n} + \tau [(\mathbf{v}_r \times \mathbf{n}) \times \mathbf{n}] / |(\mathbf{v}_r \times \mathbf{n}) \times \mathbf{n}|) dS, \quad (5)$$

and the gravitational and aerodynamic torque $\mathbf{M}_N$ and $\mathbf{M}_a$ are given by

$$\mathbf{M}_N = -G m_N \rho_d \int_{V_d} \frac{l \times \mathbf{r}}{r^3} d\Gamma \quad (6)$$

$$\mathbf{M}_a = -\int l \times (p\mathbf{n} + \tau [(\mathbf{v}_r \times \mathbf{n}) \times \mathbf{n}] / |(\mathbf{v}_r \times \mathbf{n}) \times \mathbf{n}|) dS, \quad (7)$$

where $p, \tau$ are the gas pressure and the shear stress of the elementary surface element with area $dS$, $\mathbf{n}$ is the unit vector of outward normal to the element, $\mathbf{v}_r$ is the gas–grain relative velocity, $\rho_d$ is the specific density of the dust grain, $d\Gamma$ is the elementary volume of the grain, $V_d$ is the volume of the grain, and $l$ is the radius vector from grain's centre of mass to $d\Gamma$ or $dS$. In case of a cometary atmosphere, the mean-free path of the molecules $\lambda$ is much larger than the grain size $a_d$ (i.e. $\lambda \gg a_d$). Therefore, the pressure $p$ and the shear stress $\tau$ can be computed as in a free molecular flow:

$$p/p_g = \left[ \begin{aligned} & s' \cos \beta / \sqrt{\pi} + \frac{1}{2} \sqrt{\frac{T_d}{T_g}} \\ & + \left[ \left( \frac{1}{2} + s'^2 \cos^2 \beta \right) + \frac{1}{2} \sqrt{\pi} s' \cos \beta \sqrt{\frac{T_d}{T_g}} \right] \\ & \times \left[ 1 + \erf(s' \cos \beta) \right] \end{aligned} \right] \exp(-s'^2 \cos^2 \beta) \quad (8)$$

$$\tau/p_g = -\frac{s' \sin \beta}{\sqrt{\pi}} \times \left[ \exp(-s'^2 \cos^2 \beta) + \sqrt{\pi} s' \cos \beta \{1 + \erf(s' \cos \beta)\} \right]. \quad (9)$$

Here it is assumed that molecules are reflected diffusely, $p_g = \rho v_r^2 / (2s'^2)$, $s' = v_r' \sqrt{m_g / (2k_B T_g)}$ is gas–dust speed ratio, $v_r' = v_g - v_d$ being the gas–dust relative velocity ($v_d$ is the velocity of the surface element accounting rotation of the grain), $\beta$ is an angle between $v_r'$ and inward normal of the element. The pressure $p$ is co-directional to the inward normal of the element, the shear stress $\tau$ is co-directional with the projection of $v_r'$ on the element's plane.

For given shapes of grains the integrals in equations (5) and (7) can be pre-computed for different orientations of the grain with respect to the flow and various values of $T_g/T_g$ and $s'$. Then the approximate expressions for the aerodynamic force and torque are

$$\mathbf{F}_a = \frac{1}{2} (\mathbf{v}_g - \mathbf{v}_d)^2 \rho_d \sigma_d C_D \frac{\mathbf{v}_g - \mathbf{v}_d}{|\mathbf{v}_g - \mathbf{v}_d|} \quad (10)$$

$$\mathbf{M}_a = \frac{1}{2} \rho_g (\mathbf{v}_g - \mathbf{v}_d)^2 \sigma_d a_d \mathbf{C}_M, \quad (11)$$

where $C_D$ and $C_M$ are the averaged coefficients of aerodynamic drag and torque, $a_d$ and $\sigma_d$ are the characteristic size and cross-section of the particle.

This paper neglects the solar pressure and the solar tidal forces that affect significantly the dust motion out of the acceleration zone over long times (weeks and longer). In addition, for particles equal and larger than mm-sized these effects are important also just outside the dust acceleration zone. The rigorous computation in the heliocentric
reference frame of the dust motion out of cometary Hill's sphere can be found in Fulle et al. (1995).

Following Zakharov et al. (2018) and Zakharov et al. (2021) we introduce the following dimensionless variables: $\tilde{v}_g = v_g/v_g^{\text{max}}$, $\tilde{v}_d = v_d/v_g^{\text{max}}$, $\tilde{r} = r/R_N$, $\tilde{t} = t/\Delta t$, $\tilde{\rho}_g = \rho_g/\rho_s$. Here $v_g^{\text{max}} = \sqrt{\gamma \frac{\gamma+1}{\gamma-1} \frac{k_b}{m_g} T_s}$ is the theoretical maximum velocity of gas expansion, $\rho_s$ and $T_s$ are the gas density and temperature on the sonic surface (i.e. where the gas velocity is equal the local sound velocity $v_* = \sqrt{\gamma k_B T_g/m_g}$) at the subsolar point, $\gamma$ is the specific heat ratio, $R_N$ is the characteristic linear scale (e.g. the radius of the nucleus), and $\Delta t = R_N/v_g^{\text{max}}$. In addition, we introduce $\tilde{g}$ via normalizing $g$ on $Gm_N/R_N$ (therefore, in the case of a spherical nucleus $\tilde{g} = -\tilde{r}/|\tilde{r}^3|)$, $\tilde{I} = I/(m_d a_d^2)$ and $\tilde{\omega} = \omega/(v_g^{\text{max}}/a_d)$.

Using the dimensionless variables, the equations of translational motion of dust grain can be rewritten in dimensionless form

$$ \frac{d\tilde{v}_d}{d\tilde{t}} = \tilde{\rho}_g(\tilde{v}_g - \tilde{v}_d)^2 \frac{\tilde{v}_g - \tilde{v}_d}{|\tilde{v}_g - \tilde{v}_d|} C_D \text{ Iv} + \text{Fu} \tilde{g}. \quad (12) $$

$$ \frac{d\tilde{r}}{d\tilde{t}} = \tilde{v}_d. \quad (13) $$

The equation of the rotational motion in the dimensionless form is

$$ \tilde{I} \frac{d\tilde{\omega}}{d\tilde{t}} + \tilde{\omega} \times \tilde{I} \tilde{\omega} \frac{R_N}{a_d} = \tilde{\rho}_g (\tilde{v}_g - \tilde{v}_d)^2 C_M \text{ Iv} + \text{Fu} \int \frac{\tilde{I} \times \tilde{r}}{\tilde{r}^3} d\tilde{\Gamma}. \quad (14) $$

These equations contain three dimensionless parameters Iv, Fu, and $R_N/a_d$. The parameter Iv characterizes the efficiency of entrainment of the particle within the gas flow (i.e. the ability of a dust particle to adjust to the gas velocity)

$$ \text{Iv} = \frac{1}{2} \frac{\rho_s \sigma_d R_N}{m_d} = \frac{3}{8} \frac{\rho_s R_N}{a_d \rho_d} = \frac{3 Q_g m_g}{32 R_N a_d \rho_d \pi \sqrt{T_s \gamma k_B / m_g}}. \quad (15) $$

Here, $Q_g$ is the total gas production rate [$s^{-1}$] of a spherical nucleus of radius $R_N$ with uniform gas production over the surface equal to the gas production at subsolar point. The parameter Fu characterizes the efficiency of gravitational attraction

$$ \text{Fu} = \frac{Gm_N}{R_N} \frac{1}{(v_g^{\text{max}})^2}. \quad (16) $$

The dimensionless parameters Iv, Fu for the first time were introduced in Zakharov et al. (2018). In order to define Iv, Fu it is necessary to know: $m_g$, $\gamma$, $\rho_s$, $T_s$, $R_N$, $m_N$, $\sigma_d$, and $m_d$ (or $a_d$ and $\rho_d$).

### 3 RESULTS

In a given gas flow field, if two dust grains have same Iv and Fu their translational motion is similar. If, in addition, they have same $R_N/a_d$ then their rotational motion is similar as well. We remind that the influence of solar pressure is assumed negligible.

In the following section, we derive scaling laws for the terminal dust velocity $\tilde{v}_d^t$ and rotational frequency $\tilde{\omega}^t$, that is, the velocity and rotational frequency that a dust grain has after decoupling with the gas flow. As was derived from observational, the dust acceleration limited within six nuclear radii for a broad range of particle sizes (Gerig et al. 2018). The numerical studies (e.g. Zakharov et al. 2018; Ivanovski 2017) also show that gas and dust decoupling occurs within first ten radii of the nucleus.

At distances $r/R_N > 10$ the time derivatives $\frac{d\tilde{\omega}}{d\tilde{t}} \to 0$ and $\frac{d\tilde{v}_d}{d\tilde{t}} \to 0$, the gas velocity is close to $v_g^{\text{max}}$ therefore $\tilde{v}_g \approx 1$. As was shown in Zakharov et al. (2021) for the range of conditions Iv < 0.1 and Fu <

<figure><img src="image_3.png" /><figcaption>**Figure 2.** Model particles used in the dynamical calculations in Ivanovski (2017). Prolate (right) and oblate (left).</figcaption></figure>

0.01 · Iv, the ratio of terminal velocities of two dust particles (labelled 0 and 1) in a given gas flow field is

$$ \frac{\tilde{v}_{d1}^t}{\tilde{v}_{d0}^t} = \sqrt{\frac{\text{Iv}_1}{\text{Iv}_0}}. \quad (17) $$

Therefore, the scaling law for the dust particle terminal velocity is

$$ v_{d1}^t = v_{d0}^t \sqrt{\frac{Q_g m_g}{R_{N0} a_{d0} \rho_{d0} v_{*0}}} \frac{R_{N0} a_{d0} \rho_{d0} v_{*0}}{Q_g m_g}, \quad (18) $$

where $\tilde{v}_{d1}^t = v_{d1}^t / v_g^{\text{max}}$.

The scaling law for the dust particle rotational frequency follows from the equation (14) under assumption of negligible influence of gravitational torque (Iv $\gg$ Fu)

$$ \frac{\tilde{\omega}_0^t}{\tilde{\omega}_1^t} = \frac{1 - \tilde{v}_{d1}^t}{1 - \tilde{v}_{d0}^t} \sqrt{\frac{\text{Iv}_1 a_{d1} R_{N0}}{\text{Iv}_0 a_{d0} R_{N1}}} \quad (19) $$

$$ \tilde{\omega}_1^t = \tilde{\omega}_0^t \frac{1 - \tilde{v}_{d1}^t}{1 - \tilde{v}_{d0}^t} \sqrt{\frac{\text{Iv}_1 a_{d1} R_{N0}}{\text{Iv}_0 a_{d0} R_{N1}}} = \tilde{\omega}_0^t \frac{1 - \tilde{v}_{d1}^t}{1 - \tilde{v}_{d0}^t} \sqrt{\frac{Q_1 m_{g1} v_{*0} \rho_{d0} R_{N0}^2}{Q_0 m_{g0} v_{*1} \rho_{d1} R_{N1}^2}}, \quad (20) $$

where $\tilde{\omega}_1^t = \omega_0^t a_{d1} / v_g^{\text{max}}$. The value of $\tilde{v}_d^t$ is taken from equation (18) or, if $\tilde{v}_d^t \ll 1$, the equation (20) becomes

$$ \tilde{\omega}_1^t = \tilde{\omega}_0^t \sqrt{\frac{Q_1 m_{g1} v_{*0} \rho_{d0} R_{N0}^2}{Q_0 m_{g0} v_{*1} \rho_{d1} R_{N1}^2}}. \quad (21) $$

In order to verify these scaling relations we use the numerically computed data from Ivanovski (2017) and Moreno et al. (2022). Ellipsoids of revolution with different aspect ratios were considered in Ivanovski (2017) (see Fig. 2). Irregularly shaped particles were considered in Moreno et al. (2022) (see Fig. 3). The gas coma was approximated by a spherically symmetric expanding gas (pure $H_2O$) flow. It was assumed that the nucleus has radius $R_N = 2$ [km], mass $m_N = 10^{13}$ [kg] and the surface temperature $T_s = 200$ [K].

At first, we test the scaling on the results of numerical simulations given in Ivanovski (2017). As the referential case for the scaling we use numerical solution for $a_d = 10^{-4}$ [m], $Q_g = 10^{26}$ [s$^{-1}$], $\rho_d = 100$ [kg m$^{-3}$], which gives $v_d^t = 4.2$ [m s$^{-1}$] and $\omega^t = 0.09$ [s$^{-1}$] for the prolate spheroid (case a03g1d1pr), and $v_d^t = 5.0$ [m s$^{-1}$] and $\omega^t = 0.1$ [s$^{-1}$] for the oblate spheroid (case a03g1d1ob). This case corresponds to Iv = 1.280 · 10$^{-5}$ and Fu = 3.860 · 10$^{-7}$ (i.e. Iv/Fu = 33.16). Tables 1 and 2 show the comparison of scaled (via eqs. (17) and (20)) and numerically computed values from Ivanovski
<figure><img src="image_4.png" /><figcaption>**Figure 3.** Model particles used in the dynamical calculations in Moreno et al. (2022). On the top row of panels, BEGONIA shape model particle, having a flattened shape, on the middle row GLASENAPPIA shape model particle (elongated), and on the bottom row, AURORA shape model particle displaying a rounded shape. Approximate outer dimensions of the particles relative are (1.00 × 0.72 × 0.36), (1.00 × 0.43 × 0.36), and (1.00 × 1.09 × 1.14), respectively.</figcaption></figure>

<table><caption>**Table 1.** Comparison of scaled (via eqs. (17) and (20) with numerically computed values from Ivanovski (2017). The referential case is $a_03g1d1pr$ (prolate spheroid) with $a_d = 10^{-4}$ [m], $Q_g = 10^{26}$ [s$^{-1}$], $\rho_d = 100$ [kg m$^{-3}$], $v'_d = 4.2$ [m s$^{-1}$], $\omega' = 0.09$ [s$^{-1}$]. The case's labels correspond to the Tables B.4–B.6 in Ivanovski (2017).</caption><thead><tr><th rowspan="2">Case</th><th rowspan="2">a<sub>d</sub><br>[m]</th><th rowspan="2">Q<sub>g</sub><br>[s<sup>-1</sup>]</th><th rowspan="2">ρ<sub>d</sub><br>[kg m<sup>-3</sup>]</th><th rowspan="2">Iv</th><th rowspan="2">Iv/Fu</th><th colspan="2">v'<sub>d</sub>[m s<sup>-1</sup>]</th><th colspan="2">ω'[s<sup>-1</sup>]</th><th rowspan="2">ε<sub>v</sub><br>per cent</th><th rowspan="2">ε<sub>ω</sub><br>per cent</th></tr><tr><th>Scaled</th><th>Computed</th><th>Scaled</th><th>Computed</th></tr></thead><tbody><tr><td>a03g1d2pr</td><td>10<sup>-4</sup></td><td>10<sup>26</sup></td><td>1000.0</td><td>1.28·10<sup>-6</sup></td><td>3.315</td><td>1.328</td><td>1.0</td><td>0.0285</td><td>0.02</td><td>32.80</td><td>42.50</td></tr><tr><td>a03g2d1pr</td><td>10<sup>-4</sup></td><td>10<sup>27</sup></td><td>100.0</td><td>1.28·10<sup>-4</sup></td><td>3.315·10<sup>2</sup></td><td>13.282</td><td>13.7</td><td>0.2818</td><td>0.23</td><td>3.05</td><td>22.52</td></tr><tr><td>a03g2d2pr</td><td>10<sup>-4</sup></td><td>10<sup>27</sup></td><td>1000.0</td><td>1.28·10<sup>-5</sup></td><td>3.315·10<sup>1</sup></td><td>4.200</td><td>4.2</td><td>0.0900</td><td>0.09</td><td>0.0</td><td>0.0</td></tr><tr><td>a03g3d1pr</td><td>10<sup>-4</sup></td><td>10<sup>28</sup></td><td>100.0</td><td>1.28·10<sup>-3</sup></td><td>3.315·10<sup>3</sup></td><td>42.000</td><td>43.0</td><td>0.8632</td><td>0.7</td><td>2.33</td><td>23.31</td></tr><tr><td>a03g3d2pr</td><td>10<sup>-4</sup></td><td>10<sup>28</sup></td><td>1000.0</td><td>1.28·10<sup>-4</sup></td><td>3.315·10<sup>2</sup></td><td>13.282</td><td>13.7</td><td>0.2818</td><td>0.23</td><td>3.05</td><td>22.52</td></tr><tr><td>a03g4d1pr</td><td>10<sup>-4</sup></td><td>10<sup>29</sup></td><td>100.0</td><td>1.28·10<sup>-2</sup></td><td>3.315·10<sup>4</sup></td><td>132.816</td><td>127.5</td><td>2.4504</td><td>2.0</td><td>4.17</td><td>22.52</td></tr><tr><td>a03g4d2pr</td><td>10<sup>-4</sup></td><td>10<sup>29</sup></td><td>1000.0</td><td>1.28·10<sup>-3</sup></td><td>3.315·10<sup>3</sup></td><td>42.000</td><td>43.0</td><td>0.8632</td><td>0.7</td><td>2.33</td><td>23.31</td></tr><tr><td>b03g1d1pr</td><td>10<sup>-3</sup></td><td>10<sup>26</sup></td><td>100.0</td><td>1.28·10<sup>-6</sup></td><td>3.315</td><td>1.328</td><td>1.0</td><td>0.0090</td><td>0.01</td><td>32.80</td><td>10.00</td></tr><tr><td>b03g2d1pr</td><td>10<sup>-3</sup></td><td>10<sup>27</sup></td><td>100.0</td><td>1.28·10<sup>-5</sup></td><td>3.315·10<sup>1</sup></td><td>4.200</td><td>4.2</td><td>0.0285</td><td>0.025</td><td>0.0</td><td>14.00</td></tr><tr><td>b03g2d2pr</td><td>10<sup>-3</sup></td><td>10<sup>27</sup></td><td>1000.0</td><td>1.28·10<sup>-6</sup></td><td>3.315</td><td>1.328</td><td>1.0</td><td>0.0090</td><td>0.008</td><td>32.80</td><td>12.50</td></tr><tr><td>b03g3d1pr</td><td>10<sup>-3</sup></td><td>10<sup>28</sup></td><td>100.0</td><td>1.28·10<sup>-4</sup></td><td>3.315·10<sup>2</sup></td><td>13.282</td><td>13.5</td><td>0.0891</td><td>0.08</td><td>1.62</td><td>11.37</td></tr><tr><td>b03g3d2pr</td><td>10<sup>-3</sup></td><td>10<sup>28</sup></td><td>1000.0</td><td>1.28·10<sup>-5</sup></td><td>3.315·10<sup>1</sup></td><td>4.200</td><td>4.2</td><td>0.0285</td><td>0.025</td><td>0.0</td><td>14.00</td></tr><tr><td>b03g4d1pr</td><td>10<sup>-3</sup></td><td>10<sup>29</sup></td><td>100.0</td><td>1.28·10<sup>-3</sup></td><td>3.315·10<sup>3</sup></td><td>42.000</td><td>43.0</td><td>0.2730</td><td>0.2</td><td>2.33</td><td>36.50</td></tr><tr><td>b03g4d2pr</td><td>10<sup>-3</sup></td><td>10<sup>29</sup></td><td>1000.0</td><td>1.28·10<sup>-4</sup></td><td>3.315·10<sup>2</sup></td><td>13.282</td><td>14.0</td><td>0.0891</td><td>0.08</td><td>5.13</td><td>11.37</td></tr><tr><td>c03g2d1pr</td><td>10<sup>-2</sup></td><td>10<sup>27</sup></td><td>100.0</td><td>1.28·10<sup>-6</sup></td><td>3.315</td><td>1.328</td><td>1.0</td><td>0.0029</td><td>0.002</td><td>32.80</td><td>45.00</td></tr><tr><td>c03g3d1pr</td><td>10<sup>-2</sup></td><td>10<sup>28</sup></td><td>100.0</td><td>1.28·10<sup>-5</sup></td><td>3.315·10<sup>1</sup></td><td>4.200</td><td>4.2</td><td>0.0090</td><td>0.008</td><td>0.0</td><td>12.50</td></tr><tr><td>c03g3d2pr</td><td>10<sup>-2</sup></td><td>10<sup>28</sup></td><td>1000.0</td><td>1.28·10<sup>-6</sup></td><td>3.315</td><td>1.328</td><td>1.0</td><td>0.0029</td><td>0.003</td><td>32.80</td><td>3.33</td></tr><tr><td>c03g4d1pr</td><td>10<sup>-2</sup></td><td>10<sup>29</sup></td><td>100.0</td><td>1.28·10<sup>-4</sup></td><td>3.315·10<sup>2</sup></td><td>13.282</td><td>13.7</td><td>0.0282</td><td>0.03</td><td>3.05</td><td>6.00</td></tr><tr><td>c03g4d2pr</td><td>10<sup>-2</sup></td><td>10<sup>29</sup></td><td>1000.0</td><td>1.28·10<sup>-5</sup></td><td>3.315·10<sup>1</sup></td><td>4.200</td><td>4.2</td><td>0.0090</td><td>0.008</td><td>0.0</td><td>12.50</td></tr></tbody></table>
<table><caption>**Table 2.** Comparison of scaled (via eqs. (17) and (20)) with numerically computed values from Ivanovski (2017). The referential case is *a03g1d1ob* (oblate spheroid) with $a_d = 10^{-4}$ [m], $Q_g = 10^{26}$ [s$^{-1}$], $\rho_d = 100$ [kg m$^{-3}$], $v_d^t = 5.0$ [m s$^{-1}$], $\omega^t = 0.1$ [s$^{-1}$]. The case's labels correspond to the Tables B.8–B.10 in Ivanovski (2017).</caption><thead><tr><th rowspan="2">Case</th><th rowspan="2">a<sub>d</sub><br>[m]</th><th rowspan="2">Q<sub>g</sub><br>[s<sup>-1</sup>]</th><th rowspan="2">ρ<sub>d</sub><br>[kg m<sup>-3</sup>]</th><th rowspan="2">Iv</th><th rowspan="2">Iv/Fu</th><th colspan="2">v<sub>d</sub><sup>t</sup>[m s<sup>-1</sup>]</th><th colspan="2">ω<sup>t</sup>[s<sup>-1</sup>]</th><th rowspan="2">ε<sub>v</sub><br>per cent</th><th rowspan="2">ε<sub>ω</sub><br>per cent</th></tr><tr><th>Scaled</th><th>Computed</th><th>Scaled</th><th>Computed</th></tr></thead><tbody><tr><td>a03g1d2ob</td><td>10<sup>-4</sup></td><td>10<sup>26</sup></td><td>1000.0</td><td>1.28·10<sup>-6</sup></td><td>3.315</td><td>1.581</td><td>1.3</td><td>0.0317</td><td>0.03</td><td>21.6</td><td>5.80</td></tr><tr><td>a03g2d1ob</td><td>10<sup>-4</sup></td><td>10<sup>27</sup></td><td>100.0</td><td>1.28·10<sup>-4</sup></td><td>3.315·10<sup>2</sup></td><td>15.811</td><td>14.0</td><td>0.3125</td><td>0.30</td><td>12.9</td><td>4.18</td></tr><tr><td>a03g2d2ob</td><td>10<sup>-4</sup></td><td>10<sup>27</sup></td><td>1000.0</td><td>1.28·10<sup>-5</sup></td><td>3.315·10<sup>1</sup></td><td>5.000</td><td>5.0</td><td>0.1000</td><td>0.10</td><td>0.00</td><td>0.00</td></tr><tr><td>a03g3d1ob</td><td>10<sup>-4</sup></td><td>10<sup>28</sup></td><td>100.0</td><td>1.28·10<sup>-3</sup></td><td>3.315·10<sup>3</sup></td><td>50.000</td><td>48.0</td><td>0.9513</td><td>1.00</td><td>4.17</td><td>4.87</td></tr><tr><td>a03g3d2ob</td><td>10<sup>-4</sup></td><td>10<sup>28</sup></td><td>1000.0</td><td>1.28·10<sup>-4</sup></td><td>3.315·10<sup>2</sup></td><td>15.811</td><td>15.0</td><td>0.3125</td><td>0.10</td><td>5.41</td><td>213.0</td></tr><tr><td>a03g4d1ob</td><td>10<sup>-4</sup></td><td>10<sup>29</sup></td><td>100.0</td><td>1.28·10<sup>-2</sup></td><td>3.315·10<sup>4</sup></td><td>158.114</td><td>140.0</td><td>2.6385</td><td>2.20</td><td>12.9</td><td>19.9</td></tr><tr><td>a03g4d2ob</td><td>10<sup>-4</sup></td><td>10<sup>29</sup></td><td>1000.0</td><td>1.28·10<sup>-3</sup></td><td>3.315·10<sup>3</sup></td><td>50.000</td><td>48.0</td><td>0.9513</td><td>1.00</td><td>4.17</td><td>4.87</td></tr><tr><td>b03g1d1ob</td><td>10<sup>-3</sup></td><td>10<sup>26</sup></td><td>100.0</td><td>1.28·10<sup>-6</sup></td><td>3.315</td><td>1.581</td><td>1.3</td><td>0.0100</td><td>0.01</td><td>21.6</td><td>0.370</td></tr><tr><td>b03g2d1ob</td><td>10<sup>-3</sup></td><td>10<sup>27</sup></td><td>100.0</td><td>1.28·10<sup>-5</sup></td><td>3.315·10<sup>1</sup></td><td>5.000</td><td>4.8</td><td>0.0316</td><td>0.03</td><td>4.17</td><td>5.41</td></tr><tr><td>b03g2d2ob</td><td>10<sup>-3</sup></td><td>10<sup>27</sup></td><td>1000.0</td><td>1.28·10<sup>-6</sup></td><td>3.315</td><td>1.581</td><td>1.3</td><td>0.0100</td><td>0.01</td><td>21.6</td><td>0.370</td></tr><tr><td>b03g3d1ob</td><td>10<sup>-3</sup></td><td>10<sup>28</sup></td><td>100.0</td><td>1.28·10<sup>-4</sup></td><td>3.315·10<sup>2</sup></td><td>15.811</td><td>15.0</td><td>0.0988</td><td>0.10</td><td>5.41</td><td>1.17</td></tr><tr><td>b03g3d2ob</td><td>10<sup>-3</sup></td><td>10<sup>28</sup></td><td>1000.0</td><td>1.28·10<sup>-5</sup></td><td>3.315·10<sup>1</sup></td><td>5.000</td><td>4.8</td><td>0.0316</td><td>0.03</td><td>4.17</td><td>5.41</td></tr><tr><td>b03g4d1ob</td><td>10<sup>-3</sup></td><td>10<sup>29</sup></td><td>100.0</td><td>1.28·10<sup>-3</sup></td><td>3.315·10<sup>3</sup></td><td>50.000</td><td>48.0</td><td>0.3008</td><td>0.30</td><td>4.17</td><td>0.278</td></tr><tr><td>b03g4d2ob</td><td>10<sup>-3</sup></td><td>10<sup>29</sup></td><td>1000.0</td><td>1.28·10<sup>-4</sup></td><td>3.315·10<sup>2</sup></td><td>15.811</td><td>15.0</td><td>0.0988</td><td>0.10</td><td>5.41</td><td>1.17</td></tr><tr><td>c03g2d1ob</td><td>10<sup>-2</sup></td><td>10<sup>27</sup></td><td>100.0</td><td>1.28·10<sup>-6</sup></td><td>3.315</td><td>1.581</td><td>1.2</td><td>0.0032</td><td>0.003</td><td>31.8</td><td>5.80</td></tr><tr><td>c03g3d1ob</td><td>10<sup>-2</sup></td><td>10<sup>28</sup></td><td>100.0</td><td>1.28·10<sup>-5</sup></td><td>3.315·10<sup>1</sup></td><td>5.000</td><td>4.6</td><td>0.0100</td><td>0.010</td><td>8.70</td><td>0.00</td></tr><tr><td>c03g3d2ob</td><td>10<sup>-2</sup></td><td>10<sup>28</sup></td><td>1000.0</td><td>1.28·10<sup>-6</sup></td><td>3.315</td><td>1.581</td><td>1.2</td><td>0.0032</td><td>0.003</td><td>31.8</td><td>5.80</td></tr><tr><td>c03g4d1ob</td><td>10<sup>-2</sup></td><td>10<sup>29</sup></td><td>100.0</td><td>1.28·10<sup>-4</sup></td><td>3.315·10<sup>2</sup></td><td>15.811</td><td>15.5</td><td>0.0313</td><td>0.030</td><td>2.01</td><td>4.18</td></tr><tr><td>c03g4d2ob</td><td>10<sup>-2</sup></td><td>10<sup>29</sup></td><td>1000.0</td><td>1.28·10<sup>-5</sup></td><td>3.315·10<sup>1</sup></td><td>5.000</td><td>5.0</td><td>0.0100</td><td>0.008</td><td>0.00</td><td>25.0</td></tr></tbody></table>

<figure><img src="image_5.png" /><figcaption>**Figure 4.** Left: scaling precision of the velocity and the rotation for the prolate cases in Table 1 and right: scaling precision of the velocity and the rotation for the oblate cases in Table 2.</figcaption></figure>

(2017) for different values of $a_d$, $Q_g$, and $\rho_d$. Using the two groups of results in Tables 1 and 2, we summarize the corresponding scaling precision in Fig. 4, left and 4, right, respectively. We note that in the cases of prolate spheroids (Table 1 and Fig. 4, right) high-gas production rate and small particles lead to better precision in scaling for the velocity than for the rotation. Nevertheless, the rotation precision is within maximum 40 per cent. In the cases of oblate spheroid (Table 2 and Fig. 4, right) an increase in the gas production rate produces better precision in scaling for the velocity than for the rotation. All this trend is expected since higher gas production rate is more efficient for particles with large cross-sections than

the ones of smaller sizes, provided equal particle masses for all considered particles. Furthermore, in case of large particles we have an opposite trend, less precision in velocity scaling, and better one for rotation. Only 2 cases of 38 considered show precision worse than 35 per cent and one worse than 100 per cent for the rotation. However, for the remaining cases, we note that the scaling achieves precision of more than 60–70 per cent for both prolate and oblate spheroids and more than 95 per cent for the velocity of prolate spheroids.

The precision of the scaling is estimated by the difference of numerically computed and scaled values with respect to numerically
<table><caption>**Table 3.** Comparison of scaled (via eqs. (17) and (20)) with numerically computed values from Moreno et al. (2022). The referential case is $a_d = 10^{-3}$ [m], $Q_g = 3 \cdot 10^{28}$ [s$^{-1}$], $\varrho_d = 800$ [kg m$^{-3}$] which gives for $v'_d$ [m s$^{-1}$] and $\omega^t$ [s$^{-1}$] the following values: 10.0 and 1.6 for Begonia; 10.0 and 2.7 for Glasenappia; 9.34 and 0.9 for Aurora.</caption><thead><tr><th rowspan="2">Case</th><th rowspan="2">a<sub>d</sub><br>[m]</th><th rowspan="2">Q<sub>g</sub><br>[s<sup>-1</sup>]</th><th rowspan="2">Iv</th><th rowspan="2">Iv/Fu</th><th colspan="2">v'<sub>d</sub>[m s<sup>-1</sup>]</th><th colspan="2">ω<sup>t</sup>[s<sup>-1</sup>]</th><th rowspan="2">ε<sub>v</sub><br>per cent</th><th rowspan="2">ε<sub>ω</sub><br>per cent</th></tr><tr><th>Scaled</th><th>Computed</th><th>Scaled</th><th>Computed</th></tr></thead><tbody><tr><td>Begonia</td><td>10<sup>-3</sup></td><td>1.10<sup>27</sup></td><td>1.6·10<sup>-6</sup></td><td>4.144</td><td>1.83</td><td>1.6</td><td>0.29</td><td>0.3</td><td>14.37</td><td>3.33</td></tr><tr><td>Glasenappia</td><td>10<sup>-3</sup></td><td>1.10<sup>27</sup></td><td>1.6·10<sup>-6</sup></td><td>4.144</td><td>1.83</td><td>1.6</td><td>0.50</td><td>0.5</td><td>14.37</td><td>0.0</td></tr><tr><td>Aurora</td><td>10<sup>-3</sup></td><td>1.10<sup>27</sup></td><td>1.6·10<sup>-6</sup></td><td>4.144</td><td>1.71</td><td>1.45</td><td>0.17</td><td>0.15</td><td>17.93</td><td>13.33</td></tr><tr><td>Begonia</td><td>10<sup>-2</sup></td><td>3.10<sup>28</sup></td><td>4.8·10<sup>-6</sup></td><td>12.43</td><td>3.16</td><td>3.06</td><td>0.16</td><td>0.09</td><td>3.27</td><td>77.78</td></tr><tr><td>Glasenappia</td><td>10<sup>-2</sup></td><td>3.10<sup>28</sup></td><td>4.8·10<sup>-6</sup></td><td>12.43</td><td>3.16</td><td>3.10</td><td>0.27</td><td>0.20</td><td>1.94</td><td>35.00</td></tr><tr><td>Aurora</td><td>10<sup>-2</sup></td><td>3.10<sup>28</sup></td><td>4.8·10<sup>-6</sup></td><td>12.43</td><td>2.95</td><td>2.85</td><td>0.09</td><td>0.04</td><td>3.51</td><td>125.00</td></tr><tr><td>Begonia</td><td>10<sup>-4</sup></td><td>3.10<sup>28</sup></td><td>4.8·10<sup>-4</sup></td><td>1.243·10<sup>3</sup></td><td>31.62</td><td>31.50</td><td>15.62</td><td>97.5</td><td>0.38</td><td>83.98</td></tr><tr><td>Glasenappia</td><td>10<sup>-5</sup></td><td>3.10<sup>28</sup></td><td>4.8·10<sup>-3</sup></td><td>1.243·10<sup>4</sup></td><td>100.00</td><td>98.20</td><td>243.57</td><td>780.00</td><td>1.83</td><td>68.77</td></tr><tr><td>Begonia</td><td>10<sup>-6</sup></td><td>3.10<sup>28</sup></td><td>4.8·10<sup>-2</sup></td><td>1.243·10<sup>5</sup></td><td>316.23</td><td>–</td><td>1067.13</td><td>–</td><td>–</td><td>–</td></tr><tr><td>Glasenappia</td><td>10<sup>-6</sup></td><td>3.10<sup>28</sup></td><td>4.8·10<sup>-2</sup></td><td>1.243·10<sup>5</sup></td><td>316.23</td><td>–</td><td>1800.78</td><td>–</td><td>–</td><td>–</td></tr><tr><td>Aurora</td><td>10<sup>-6</sup></td><td>3.10<sup>28</sup></td><td>4.8·10<sup>-2</sup></td><td>1.243·10<sup>5</sup></td><td>295.36</td><td>–</td><td>620.24</td><td>–</td><td>–</td><td>–</td></tr></tbody></table>

computed value

$$ \epsilon_v = \frac{|v_{\text{computed}} - v_{\text{scaled}}|}{v_{\text{computed}}} \qquad (22) $$

$$ \epsilon_{\omega} = \frac{|\omega_{\text{computed}} - \omega_{\text{scaled}}|}{\omega_{\text{computed}}}. \qquad (23) $$

For both prolate and oblate grains the scaled values of $v'_d$ and $\omega^t$ are in good agreement with numerically computed ones if Iv/Fu > 10. The large difference in scaled and numerically computed frequency in case a03g3d2ob is probably due to inaccuracy of numerical solution.

And further, we test the scaling on the results of numerical simulations given in Moreno et al. (2022). As the referential case for the scaling we use numerical solution for $a_d = 10^{-3}$ [m], $Q_g = 3 \cdot 10^{28}$ [s$^{-1}$], $\varrho_d = 800$ [kg m$^{-3}$], which gives $v'_d = 10.0$ [m s$^{-1}$] and $\omega^t = 1.6$ [s$^{-1}$] for the Begonia shape, and $v'_d = 10.0$ [m s$^{-1}$] and $\omega^t = 2.7$ [s$^{-1}$] for the Glasenappia shape, and $v'_d = 9.34$ [m s$^{-1}$] and $\omega^t = 0.9$ [s$^{-1}$] for the Aurora shape. This case corresponds to Iv = 4.799·10$^{-5}$ and Fu = 3.860·10$^{-7}$ (i.e. Iv/Fu = 124.3). Table 3 shows the comparison of scaled (via eqs. (17) and (20)) and numerically computed values from Moreno et al. (2022) for different values of $a_d$ and $Q_g$.

Table 3 shows that precision of scaling of $v'_d$ increases with increase of Iv/Fu.

In order to estimate the precision of the scaling with respect to the shape irregularity we have calculated the dispersion of the moments of inertia of the considered shapes. The dispersion was calculated as a standard deviation of the three moments of inertia of a given shape divided on their mean value. The purpose of this estimation is that the dispersion calculated in this way is dimensionless and can be easily used for comparison in all available cases. The general trend is: the higher the dispersion is, the better the scaling is. The only deviation from this observation is the case of a prolate spheroid and this can be owing to the fact that in the cases of Tables 1 and 2 the rotational motion is happening only in one plane as discussed in Ivanovski (2017). In cases of Moreno et al. (2022), the dispersion for the three shapes is the following: 31 per cent for Begonia, 47 per cent for Glasenappia and 5 per cent for Aurora. As can be checked in Table 3 the scaling works better for more irregular shapes like Glasenappia and Begonia ones. The shapes with a smaller dispersion can easily and often change their axis of rotation that results in more complex rotational motion, that is, seems more difficult to be precisely obtained by scaling.

In addition, Table 3 shows scaled $v'_d$ and $\omega^t$ for micron size grains for which the numerical solution is not given in Moreno et al. (2022) due to excessive computational demands. It should be noted that too fast-rotating grains may violate the condition (10) from Ivanovski (2017) (it constrains applicability of free molecular approximation for evaluation of pressure $p$ and shear stress $\tau$).

## 4 CONCLUSIONS

In the present paper, we focus on derivation of scaling laws of rotational motion applicable for any shape of particles. We use a set of universal, dimensionless parameters characterizing the dust motion in the inner cometary coma. Based on the dimensionless description of the translational and rotational motion of dust particles we derived scaling relations for the terminal velocity and rotational frequency of non-spherical grains. The comparison of scaled values with available numerically computed data shows that the proposed scaling relations provide estimations of $v'_d$ and $\omega^t$ with sufficient precision. Therefore, having one numerically computed ‘reference case’ one can get order of magnitude estimation for other cases via simple scaling. This allows avoiding long-time and huge electric power consumption numerical simulations.

## ACKNOWLEDGEMENTS

This research was supported by the Italian Space Agency (ASI) within the ASI-INAF agreements I/032/05/0, I/024/12/0, 2020-4-HH.0, and N. 2023-14-HH.0. The work of NYB is funded by the Ministry of Science and Higher Education of the Russian Federation as part of the World-class Research Centre programme: Advanced Digital Technologies (contract no. 075-15-2022-311 dated 2022 April 20). This work was supported by the International Space Science Institute (ISSI) through the ISSI International Team ‘Characterization of cometary activity of 67P/Churyumov-Gerasimenko comet’.

## DATA AVAILABILITY

This work uses simulated data, generated as described diligently in the text. If you want to present additional material which would interrupt the flow of the main paper, it can be placed in an Appendix which appears after the list of references.
REFERENCES

Fulle M., Colangeli L., Mennella V., Rotundi A., Bussoletti E., 1995, A&A,
304, 622

Fulle M., Ivanovski S. L., Bertini I. et al., 2015, A&A, 583, 8

Gerig S.-B. et al., 2018, *Icarus*, **311**, 1

Ivanovski S. L., Zakharov V. V., Della Corte V., Crifo J.-F., Rotundi A., Fulle
M., 2017, *Icarus*, **282**, 333

Keller H. U., Barbieri C., Lamy P. L. et al., 2007, *Space Science Reviews*,
128, 433

Moreno F. et al., 2022, MNRAS, 510, 5142

Zakharov V. V., Ivanovski S. L., Crifo J.-F., Della Corte V., Rotundi A., Fulle
M., 2018, *Icarus*, **312**, 121

Zakharov V. V., Rotundi A., Della Corte V., Fulle M., Ivanovski S. L., Rodionov A. V., Bykov N. Y., 2021, *Icarus*, **364**, 114476

This paper has been typeset from a *TEX/LATEX* file prepared by the author.