Issue 
A&A
Volume 682, February 2024



Article Number  A130  
Number of page(s)  13  
Section  Planets and planetary systems  
DOI  https://doi.org/10.1051/00046361/202346970  
Published online  13 February 2024 
A semianalytical thermal model for craters with application to the craterinduced YORP effect
^{1}
Université Côte d’Azur, Observatoire de la Côte d’Azur, CNRS, Laboratoire Lagrange,
96 Bd de l’Observatoire,
06304
Nice,
France
email: wenhan.zhou@oca.eu
^{2}
University of Tokyo, Department of Systems Innovation, School of Engineering,
Tokyo,
Japan
Received:
23
May
2023
Accepted:
4
December
2023
Context. The YORP effect is the thermal torque generated by radiation from the surface of an asteroid. The effect is sensitive to surface topology, including smallscale roughness, boulders, and craters.
Aims. The aim of this paper is to develop a computationally efficient semianalytical model for the craterinduced YORP (CYORP) effect that can be used to investigate the functional dependence of this effect.
Methods. This study linearizes the thermal radiation term as a function of the temperature in the boundary condition of the heat conductivity, and obtains the temperature field in a crater over a rotational period in the form of a Fourier series, accounting for the effects of selfsheltering, selfradiation, and selfscattering. By comparison with a numerical model, we find that this semianalytical model for the CYORP effect works well for K > 0.1 Wm^{−1} K^{−1}. This semianalytical model is computationally threeordersofmagnitude more efficient than the numerical approach.
Results. We obtain the temperature field of a crater, accounting for the thermal inertia, crater shape, and crater location. We then find that the CYORP effect is negligible when the depthtodiameter ratio is smaller than 0.05. In this case, it is reasonable to assume a convex shape for YORP calculations. Varying the thermal conductivity yields a consistent value of approximately 0.01 for the spin component of the CYORP coefficient, while the obliquity component is inversely related to thermal inertia, declining from 0.004 in basalt to 0.001 in metal. The CYORP spin component peaks at an obliquity of 0°, 90°, or 180°, while the obliquity component peaks at an obliquity of around 45° or 135°. For a zaxis symmetric shape, the CYORP spin component vanishes, while the obliquity component persists. Our model confirms that the total YORP torque is damped by a few tens of percent by uniformly distributed smallscale surface roughness. Furthermore, for the first time, we calculate the change in the YORP torque at each impact on the surface of an asteroid explicitly and compute the resulting stochastic spin evolution more precisely.
Conclusions. This study shows that the CYORP effect due to smallscale surface roughness and impact craters is significant during the history of asteroids. The semianalytical method that we developed, which benefits from fast computation, offers new perspectives for future investigations of the YORP modeling of real asteroids and for the complete rotational and orbital evolution of asteroids accounting for collisions. Future research employing our CYORP model may explore the implications of spacevarying roughness distribution, roughness in binary systems, and the development of a comprehensive rotational evolution model for asteroid groups.
Key words: minor planets / asteroids: general
© The Authors 2024
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.
1 Introduction
The YarkovskyO’KeefeRadzievskiiPaddack (YORP) effect is a thermal torque that can alter the spin state of an asteroid over time (Rubincam 2000; Vokrouhlickỳ & Čapek 2002; Vokrouhlicky et al. 2015), and is caused by the asymmetric reemission of solar radiation by the irregular surface of the asteroid, resulting in a net torque that can spin up or spin down the asteroid’s rotation. The absorption of solar radiation makes no contributions to the YORP torque as it is averaged out over the spin and orbital periods for any asteroid shapes (Nesvornỳ & Vokrouhlickỳ 2008). So far, 11 asteroids showing timevarying rotational periods have been detected (Lowry et al. 2007; Taylor et al. 2007; Ďurech et al. 2022; Tian et al. 2022).
The YORP effect has important implications for an asteroid’s longterm rotational evolution. This effect can either spin down the asteroid to an extremely slow rotation, triggering a tumbling motion (Pravec et al. 2005), or spin up the asteroid to its rotation limit (e.g. rotational period of 2.2 h for rubble pile asteroids), leading to resurfacing (Sánchez & Scheeres 2020) and rotation disruption (Scheeres 2007; Fatka et al. 2020; Veras & Scheeres 2020). YORPinduced rotational disruption is supported by the observed asteroid pairs (Vokrouhlickỳ & Nesvornỳ 2008; Polishook 2014) and binary asteroids (Jacobson & Scheeres 2011; Delbo et al. 2011; Jacobson et al. 2013, 2016), including contactbinary asteroids (Rozek et al. 2019; Zegmott et al. 2021) and binary comets (Agarwal et al. 2020), which evolve under tidal effects and the binary YORP (BYORP) effect after the binary system is formed (Ćuk & Burns 2005; Steinberg & Sari 2011). Further potential observational evidence is the abnormal spin distribution and the obliquity distribution of nearEarth asteroids (Vokrouhlickỳ et al. 2003; Pravec et al. 2008; Rozitis & Green 2013b; Lupishko & Tielieusova 2014) and main belt asteroids (Lupishko et al. 2019), although a recent study points out that collisions might reproduce the observed distribution without the involvement of the YORP effect (Holsapple 2022). The YORP effect can influence the orbital evolution through the Yarkovsky effect, which is a thermal force that depends on the rotational state of the asteroid (Vokrouhlickỳ et al. 2000; Bottke et al. 2006). Therefore, understanding the YORP effect is important for correctly estimating the ages of asteroid families based on how much time is needed for the Yarkovsky effect to cause the observed orbital dispersion of their members from the original orbit following the disruption of their parent bodies (Vokrouhlickỳ et al. 2006; Ćuk et al. 2015; Carruba et al. 2014, 2015, 2016; Lowry et al. 2020; Marzari et al. 2020).
However, accurately calculating the YORP effect on a real asteroid remains a challenge, as it has been demonstrated to be highly sensitive to surface topology (Statler 2009; Breiter et al. 2009), such as uniform smallscale roughness (Rozitis & Green 2012), boulders (Golubov & Krugly 2012; Golubov et al. 2014, 2021; Ševeček et al. 2015; Golubov 2017; Golubov & Scheeres 2019; Golubov & Lipatova 2022), and craters (Zhou et al. 2022). Although the YORP torque caused by boulders and the tangential radiative force has been well studied, the YORP effect caused by concave structures has not yet been fully explored. More specifically, there are two kinds of concave structures on asteroid surfaces. The first one corresponds to craters, which result from impacts that an asteroid’s surface undergoes during its history and which can cover a large size range and be distributed in various ways. The second corresponds to surface roughness, which corresponds to uniformly distributed smallscale concave features that originate from the continuous effect of various processes that take place at the surface, such as thermal fatigue and space weathering. In our study, we consider both craters and surface roughness. While a pioneering study by Rozitis & Green (2012) used numerical simulations to investigate the effect of roughness, the high computational expense of such simulations prevents a comprehensive exploration of the functional dependence of this effect and its application to the global spin evolution of asteroid populations.
In addition to the precise calculation of the complete YORP effect, the longterm evolution of asteroids needs to account for stochastic collisions that affect this evolution, because each collision introduces a YORP torque due to the resulting crater. Bottke et al. (2015) performed a first study of the YORP effect accounting for craters caused by collisions, finding strong implications in the age estimate of asteroid families. However, their introduction of the concept of the stochastic YORP effect due to collisions assumed an arbitrary reset timescale for the YORP cycle. In reality, this reset timescale depends on the actual occurrence of each impact causing a crater on the asteroid’s surface and the resulting change in the YORP effect. In summary, collisions and the YORP and Yarkovsky effects are all coupled with each other in a way that is so far not well understood.
To account explicitly for the YORP torque caused by roughness and craters, Zhou et al. (2022) developed a semianalytical model that is computationally efficient for the craterinduced YORP (CYORP) torque. The CYORP torque is defined as the torque difference between the crater and the flat ground: (1)
In general, it takes the form of the following scaling rule with the radius of the crater R_{0} and of the asteroid radius R_{ast}: (2)
where Φ is the flux of solar radiation at the asteroid’s semimajor axis, and W is the CYORP coefficient, which in turn depends on the obliquity and irregularity of the asteroid, and the depthtodiameter ratio, location, thermal inertia, and albedo of the crater. As a first step, this model assumed a zero thermal conductivity.
Zhou et al. (2022) found that roughness or craters that cover 10% of the asteroid’s surface area could produce a CYORP torque comparable to the normal YORP (NYORP) torque, which arises from the macroscopic shape, with ignorance of the fine surface structure. Based on this number, which assumes that all craters have a depthtodiameter ratio of 0.16, the reset of the YORP torque by the CYORP torque could be as short as 0.4 Myr.
However, the effects of finite thermal conductivity, selfradiation, and selfscattering were not considered in Zhou et al. (2022). In the present paper, we propose a semianalytical model that accounts for the effects of selfsheltering, selfradiation, selfscattering, and nonzero thermal conductivity. This new model allows a more general exploration of the functional dependence of the CYORP effect. Moreover, as this semianalytical model is much more efficient than a purely numerical one, it can be used to study the combined influence of the YORP effect, collisions, and the Yarkovsky effect by incorporating the CYORP effect into the standard evolution model of asteroid families. We assume that the craters considered in this work are on a convex asteroid. It is possible to apply the model derived in this work to a moderately concave asteroid by approximating concave structures as craters, but this is beyond the scope of this paper.
In this paper, we describe our analytical model for the temperature field in a crater in Sect. 2. In Sect. 3, we introduce the numerical model that we developed to validate the analytical model. The main results are shown in Sects. 4 and 5. In Sect. 4, we analyse the effects of selfsheltering, selfradiation, and selfscattering of the crater, and in Sect. 5 we show the results of the calculation of the CYORP torque as a function of different variables. For the purpose of illustration, in Sect. 6 we provide an example of the analysis of CYORP considering a specific real asteroid shape and its surface roughness as well as the consequential rotational evolution. In Sect. 7, we summarise our main findings and draw conclusions.
2 Analytical model
2.1 Calculation of temperature distribution in a crater
2.1.1 Linearized analytical solution
The temperature u for the surface and the layer beneath is governed by (3)
with two boundary conditions, (4) (5)
where t is the time, z is the depth of the crater, ω is the angular velocity, K is the thermal conductivity, C is the specific heat capacity, ρ is the bulk density of the asteroid, e is the emissivity, and σ is the StefanBoltzmann constant. In the following, we use the spin angle β = ωt to replace time for simplicity. The effect of the seasonal wave is marginal and is ignored in this work.
While the onedimensional heat diffusion equation with the boundary condition of radiation has no complete analytical solution, it could be solved by linearizing the fourth power of the temperature, with the assumption that the temperature does not change significantly during a rotational period. For an arbitrary point in the hemispherical crater surface, the location of which is described by the polar angle θ and the azimuthal angle ϕ, the solution of the temperature is the real part of the expression (7)
with M = K/ρC. The function F_{n} is the nthorder of the Fourier series of the absorbed radiation flux F, which is expressed as (13)
We see that the only unknown variable is the absorbed radiation flux F. The absorbed radiation on a surface element contains three parts: solar radiation E(θ, ϕ, β), radiation from the crater itself H(θ, ϕ, β), and the scattering flux from the crater G(θ, ϕ, β), which leads us to (14)
where E_{n}(θ, ϕ, β) and H_{n}(θ, ϕ, β) denote the nthorder Fourier modes of E(θ, ϕ, t) and H(θ, ϕ, t), respectively. Here the albedo A is assumed to be 0.1 for these three flux components for the sake of simplicity, although the albedo at the thermalinfrared wavelengths is almost zero. Following Eqs. (8), (9), and (14), we obtain (15) (16)
Therefore, to solve the temperature u(θ, ϕ, β), we need to obtain the Fourier series of E(θ, ϕ, β), H(θ, ϕ, β), and G(θ, ϕ, β).
Fig. 1 Three coordinate systems in this paper: coordinate system oxyz for calculating the illuminated domain in the crater, PABC for calculating the effective radiative force of an arbitrary surface element, and OXYZ for averaging the YORP torque over the spin and orbital motion. 
2.1.2 Coordinate systems
We use three coordinate systems to calculate the radiation and the force received by the crater, as shown in Fig. 1. The coordinate system OXYZ is an inertial frame used to calculate the averaged YORP torque over the spin and orbital motion. Based on the axis OZ, we define the coordinate system oxyz fixed with the crater to calculate the instant solar flux. Finally, based on the axis oz, we define the coordinate system PABC to calculate the selfsheltering effect of the crater. The selfsheltering effect for a point in the crater refers to the nonworking moment of the photons reabsorbed by the shelter (i.e. the crater itself), which leads to the effective radiation force of the surface being tilted relative to the surface normal with a modified magnitude. This selfsheltering effect on the point P depends on the geometry of the surrounding shelter. We adopt a simple sphere with a radius of R_{1} to depict the crater shape.
In the coordinate system OXYZ, OZ points along the rotation axis and OXY is the equatorial plane. The axis OY is chosen such that the normal vector of the orbit plane lies in the plane OYZ. The axis OX follows the righthand rule. There are three crucial vectors determining the CYORP torque expressed in polar coordinates in the coordinate system OXYZ: the crater position vector r_{0} (denoted by α and β), the crater normal vector n_{0} (described by α, β, δ, and Δ), and the solar position s (described bye and η), as shown in the right panel of Fig. 1. As the polar angle is between 0° and 180° by definition, we limit 0° < α < 180° and 0° < α−Δ < 180° in our results.
The coordinate system oxyz with the origin o located at the sphere centre is fixed with the crater in order to simplify the calculation of the solar flux on the crater. The axis oz points in the opposite direction to the surface normal vector n_{0}, which is also the symmetry axis of the spherical crater. The direction of axis oy is along e_{OZ} × e_{oz} and ox follows the righthand rule. In this coordinate system, a crater with a depth of h can be defined as (17)
with cos θ_{0} = (R_{1} − h)/R_{1}.
The coordinate system PABC with the origin P at a chosen point on the crater is used to calculate the selfsheltering effect. The axis PC is along the direction of Po, PA points in the direction of e_{oz} × e_{PC}, and PB follows the righthand rule. In our code, the effective force felt by the point P is calculated first in the coordinate system PABC for simplicity and this is then transformed to the coordinate system oxyz by a rotation matrix.
2.1.3 Solar radiation E
In this section, we show how we derive the solar radiation received by an arbitrary point P in the crater, whose coordinates in oxyz system are r_{P} = (sin θ cos ϕ, sin θ sin ϕ, cos θ). The unit position vector of the Sun in the coordinate system OXYZ is described by (18)
where ϵ is the obliquity and η is the angle of orbital motion. The transform matrix between the coordinate systems oxyz and OXYZ is set to (19)
Therefore, the coordinates of the vector s in the coordinate system oxyz is (20)
On the other hand, we can use the angle λ and ϕ′ to represent s_{oxyz}: (21)
The absorbed radiation flux is (24)
Here, H is the Heaviside function defined by (29)
Using the substitution β′ = β − β_{1}, the absorbed radiation flux has the form (30)
Expanding Eq. (30) in Fourier series, we obtain (31)
The nth coefficient of the Fourier series of E(θ, ϕ) is (35)
2.1.4 Selfheating effect
Due to the concave structure of the crater, the surface element in the crateralso receives photons emitted by the crateritself, which is a process referred to as ‘selfheating’. In this section, we derive the selfradiation H and the selfscattering G as a function of the position in the crater.
Let us consider a surface element dS_{1} receiving the radiation from another surface element dS_{2}. In the reference frame oxyz, the positions of dS_{1} and dS_{2} are (37) (38)
The displacement from dS_{1} to dS_{2} is (39)
The incident angle ζ_{1} and the emission angle ζ_{2} are defined as (40) (41)
respectively. The radiation flux at the location of dS_{1} produced by the thermal radiation of dS_{2} is (42)
Substituting Eqs. (37)–(41) into Eq. (42), we obtain (43)
For an arbitrary point, the radiation flux caused by the whole crater is (44)
where 𝒵 is the crater surface, and is defined as (45)
Similarly, we obtain the selfscattering term: (46)
Therefore, both H(θ, ϕ) and G(θ, ϕ) can be expressed in terms of u(θ, ϕ) and E(θ, ϕ). As E(θ, ϕ) is derived in Sect. 2.1.3, the only unknown variable is the temperature u(θ, ϕ), the solution for which is discussed in the following section.
2.2 Solution for temperature
We obtained the Fourier series of the solar radiation flux (Sect. 2.1.3) and the radiation flux produced by the crater (Sect. 2.1.4), which allows us to return to Eq. (14) to solve the temperature distribution u(θ, ϕ). We note that the selfradiation term H contains the unknown temperature distribution that is to be solved.
The Fourier coefficients of the temperature of the crater follow (47)
Here, terms G and H are the scattering flux and selfradiation flux, respectively.
2.2.1 Solution of a general form
Equation (47) have a general form: (49)
Substituting Eq. (51) into Eq. (50), we obtain (52)
Rearranging Eq. (52), we find k is (53)
with which f(θ, ϕ) is solved out by Eq. (51).
2.2.2 Solution for temperature
In the case of u_{0}, (54) (55) (56)
and in the case of u_{n}, (57) (58) (59)
2.3 The CYORP torque
The CYORP torque is defined as the YORP torque difference between the crater and the flat ground: (60)
2.3.1 The YORP torque of the crater
The average radiative torque produced by the crater should be calculated in the inertial frame OXYZ: (61)
Here, f_{OXYZ} is the radiative force, and (62)
where n_{f} is the corrected force direction vector for each surface element. If there is no shelter around the surface element, n_{f} is equal to the surface unit normal vector n_{0}. However, the surface element in a crater has a sky sheltered by other elements, resulting in the reabsorption of the emitted photons along the direction of the shelter.
For the surface element dS(θ, ϕ), the radiative force is (63)
Without sheltering (e.g. for convex asteroids), ℋ is replaced by the hemisphere (i.e. θ ∈ (0,π/2), ϕ ∈ (0, 2π)). In this case, the force is reduced to 2Φn_{0}/3c.
2.3.2 The YORP torque of the flat portion of the surface
The absorbed radiation flux for a flat ground with the normal vector n_{0} is (64)
Here, β′ = β −β_{1}, with (65) (66) (67)
and β_{min} and β_{max} are the negative and positive values of arccos(−m_{2}/m_{1}), respectively.
3 Numerical model for examination
In the above analytical method, we adopted the assumption of a “small” temperature variation during a rotation period, which allows us to linearize the biquadrate of the temperature (i.e. ). This assumption is equivalent to a high value of the thermal parameter, which is defined as (68)
In the case of a low thermal parameter, the analytical model should be used with caution. To examine the appropriate range of the thermal inertia for which our analytical model is valid, we built a onedimensional thermophysical numerical model to perform crossvalidation.
3.1 Numerical model
We used a finite difference numerical method to solve Eq. (3) with the secondorder CrankNicholson scheme: (69)
Here, u_{i,j,k} represents the temperature at the depth of (j − 1)δz below the ith facet at the kth time step, where i, j, and k are integrals starting from 1 to i_{max}, j_{max}, and k_{max}. The coefficient c_{n} is (70)
The value of c_{n} should be smaller than 0.5 for the convergence of the iteration.
The surface temperature is determined by the first boundary condition (Eq. (4)): (71)
which can be solved by a NewtonianRaphson iterative method. The solar flux E_{i,k} on the ith facet at the kth time step is (72)
where Φ_{0} is 1364 W m^{−2} at the distance of 1 au to the Sun. Here, n_{i} is the normal vector of the ith facet and s_{k} is the unit solar position vector at the kth time step. Whether or not the facet is sheltered is judged according to the projections of other facets on the plane of the ith facet along the solar position vector s_{k}. The selfradiation term H_{i,k} is the sum of radiation from other facets: (73)
where r_{i,i′}, is the vector from the centre of the ith facet to the centre of i′th facet, and S_{i′} is the area of the ith facet. The scattering term S_{i,k} is given by (74)
The second boundary condition (Eq. (5)) translates into (75)
Combining Eqs. (69), (71), and (75), we can obtain the solution for the temperature in the next time step u_{i,j,k + 1}·,based on the temperature in the current time step u_{i,j,k} We adopted an initial temperature of u_{i,j,0} = 280 K. The maximum depth is set to be a few thermal skin depths and the total number of layers is set as j_{max} = 50. In order to make sure that the surface temperature is in an equilibrium state, we set the total time as k_{max}δt ~ 20 spin periods. We divided the considered crater into about 1000 facets and solve the temperature for each facet using the above method.
Fig. 2 Radiative force of the total crater averaged over a rotational period (8 h by default) obtained from the analytical method (solid lines) and the numerical method (dashed lines), as a function of the crater colatitude α for K = 2.65 Wm^{−1} K^{−1} (left panel) and K = 0.1 Wm^{−1} K^{−1} (right panel). The x, y, and 𝓏 components of the radiative force are shown in blue, red, and green, respectively. The distance of the crater here is 1 au from the Sun. 
3.2 Comparison with the analytical model
The thermal parameters of asteroids can vary widely depending on their composition and structure. For example, the thermal conductivity of a porous material is much lower than that of a dense metal. The thermal conductivity of stony asteroids, which are made mostly of silicates, can range from about 0.1 Wm K^{−1} to 1 WmK^{−1}, while the thermal conductivity of metallic asteroids is generally much higher, in the range of 20–50 WmK^{−1}. Asteroids that are composed of a mixture of rock and metal will have thermal conductivity values between those of pure rock and pure metal.
Here, we test three typical types of asteroid materials: regolith, solid basalt, and metal, whose properties are shown in Table 1. We calculate the radiative force of the total crater averaged over a rotational period (8 h by default), as a function of the colatitude of the crater. The craters in the test are placed at 1 au from the Sun. For simplicity, we set the obliquity to ϵ = 0. The results computed from the analytical method and the numerical method are shown in Fig. 2. We can see that the analytical result is consistent with the numerical result to a high degree, while a large discrepancy shows up when the thermal conductivity decreases to 0.1 Wm^{−1} K^{−1}. This coincides with our prerequisite of the application of our analytical method, which is that the temperature variation should be small. Therefore, our method is appropriate for solid basalt and metal materials.
Regarding regolith material, with a thermal conductivity of as low as 0.001 Wm^{−1} K^{−1}, we test the model described in Zhou et al. (2022), where zero thermal conductivity is assumed. The result is shown in Fig. 3, which indicates that this latter model works well for regolith materials. Therefore, for three materials representing asteroid surfaces, our two methods, namely the one in the present work and that described in Zhou et al. (2022), behave well in modelling the YORP effect.
Thermal parameter for three typical materials on asteroids, taken from Farinella et al. (1998).
4 Discussion on selfmodification effects
For a concave structure, there are three selfmodification effects: the selfsheltering, selfradiation, and selfscattering effect. The first one refers to the radiative force modification on the surface element due to the reabsorption of photons by the crater. The second and third ones refer to the temperature increase due to the emitted photons and scattered photons from the crater itself. Previous research shows that these selfmodification effects could be important for the YORP torque of the crater (Statler 2009; Rozitis & Green 2012, 2013a), but a quantitative description is still lacking. For example, it is not clear how deep the crater needs to be so that these selfmodification effects can no longer be ignored. This is crucial to the validity of the commonly made assumption that the asteroid can be considered as a moderately convex shape with the craters being overlooked when evaluating the YORP torque. An essential metric in this respect is the depthtodiameter ratio of the crater. We explore the radiative force (leading to the YORP torque) of a crater based on its depthtodiameter ratio under varying conditions  one for each selfmodification effect. Figure 4 shows our findings, situating the craters at the equator (α = π/2) with a spin axis of zerodegree obliquity.
Our analysis reveals that both selfradiation and selfscattering amplify the force, raising the temperature by as much as 50% when the crater mirrors a hemisphere. In contrast, selfsheltering diminishes the force by up to 30%. Notably, the impact of selfscattering remains negligible for typical asteroid surface albedos ranging between 0.1 and 0.2. When integrating all selfmodification effects, the radiative force increases by 15%.
Our result also shows that when the depthtodiameter ratio h/D_{0} < 0.05, the force increase is less than 1%. Therefore, for those shallow concave structures with h/D_{0} < 0.05, no selfmodification effects are needed in the YORP model for it to remain accurate. These can then be efficiently approximated as flat surfaces.
Fig. 3 Same as Fig. 2 but for K = 0 Wm^{−1} K^{−1} in the analytical model and K = 0.001 Wm^{−1} K^{−1} in the numerical model. 
Fig. 4 Ratio of the radiative force including selfmodification effects to that without selfmodification effects, as a function of crater depthtodiameter ratio and accounting for different selfmodification effects. 
5 Analysis of the CYORP torque
As shown in a previous work (Zhou et al. 2022), the CYORP torque depends on many factors, including the depthtodiameter ratio and location (or the colatitude α for example) of the crater as well as the obliquity and thermal parameter of the asteroid. In this section, we discuss the dependence of the CYORP torque on these factors. In the following, except for in Sect. 5.1, we assume the depthtodiameter ratio to be ~0.16.
5.1 Depthtodiameter ratio
Asteroid craters exhibit a range of distinct features in size and shape, with diameters ranging from a few centimetres to several kilometres for large asteroids. These craters generally display bowlshaped structures, containing central peaks and terraced walls when produced in the gravity regime. To simplify the modelling, a semisphere approximation is often used to represent the shape of craters.
According to the definition of CYORP torque, if the depthtodiameter ratio reaches zero, the CYORP torque is zero (Eq. (60)). Figure 5 shows the CYORP torques generated by craters with various depthtodiameter ratios. The parameters δ = ∆ = π/6, K = 2.65 Wm^{−1} K^{−1}, and α = π/2 are used. Higher depthtodiameter ratios correspond to larger spin components of the CYORP torque. A crater with h/D_{0} < 0.05 produces an insignificant CYORP torque, which may be disregarded. Furthermore, the depthtodiameter ratio also impacts the obliquity component, influencing both the torque magnitude and shape of the torque curve. For instance, when the depthtodiameter ratio is low, the asymptotic obliquity is 90°, while for higher depthtodiameter ratios, new asymptotic obliquities arise around 0° and 180°.
5.2 Crater latitude α and asteroid obliquity ϵ
Figure 6 displays the CYORP torque components as a function of the crater latitude and asteroid obliquity. The values of the parameters δ and ∆ are set to a representative value of π/4. The spin component W_{𝓏} exhibits symmetry about the axis 𝓏 = 90°, while the obliquity component W_{𝓏} is antisymmetric. The minimum and maximum values of the torque occur when the obliquity is 0^{°} or 90^{°}, with the absolute value of these extrema reaching up to 0.02. The coefficient of the obliquity component of the CYORP torque is considerably smaller, with a maximum value of 0.004.
For comparison, the typical value of the normal YORP spin coefficient is 0.005 for type I/II and <0.001 for type III/IV asteroids^{1}. The ratio of the CYORP torque to the normal YORP torque scales as (76)
Setting the ratio to 1, we find that the total area of concave structures needs to be as large as 1/4 and 1/20 of the asteroid surface area for type I/II and type III/IV asteroids, respectively.
Fig. 5 Spin component (left panel) and the obliquity component (right panel) of the CYORP coefficient as a function of the asteroid obliquity e, accounting for different depthtodiameter ratios, which are denoted by different colours. 
Fig. 6 Spin component (left panel) and the obliquity component (right panel) of the CYORP coefficient as a function of the asteroid obliquity e and the crater colatitude α(> ∆), for K = 2.65Wm^{−1} K^{−1}. Here, d/D_{0} = 0.16 and δ = ∆ = π/4. 
5.3 Thermal parameter
When the asteroid rotates quickly and has high heat conductivity, a higher value of the thermal parameter arises, resulting in a smaller variation in temperature. To explore the role of the thermal parameter, we employ the same parameters as in Sect. 5.2, but with K = 40 Wm^{−1} K^{−1} for metal materials, and the resulting CYORP torques are depicted in Fig. 7. The comparison with Fig. 6 reveals that the spin component remains relatively unchanged, while the obliquity component displays significant variation. This observation aligns with the prior assertion that the thermal parameter mainly impacts the obliquity component (Vokrouhlickỳ & Čapek 2002). In the regime of high thermal conductivity, the obliquity component diminishes as the thermal conductivity increases.
5.4 Irregularity δ and ∆
The angular parameters δ and ∆ are used to describe the irregularity of the asteroid, where δ = 0 and ∆ = 0 correspond to a perfect sphere. We explore the CYORP torque as a function of δ and ∆ with fixed asteroid obliquity and crater colatitude of π/4. The results are presented in Fig. 8. We can see that for the spin component, δ controls the torque magnitude while ∆ controls the torque direction.
Zhou et al. (2022) demonstrates that the CYORP torque vanishes for δ = 0. However, in the presence of finite thermal inertia, the obliquity component of the CYORP torque arises while the spin component remains negligible. Figure 9 illustrates the variation of the CYORP obliquity component with the asteroid obliquity and the crater colatitude when δ and ∆ are both zero. The CYORP torque has a tendency to lead the asteroid obliquity to 90° when the crater is near the poles, while it leads to an asymptotic obliquity of 0° or 180° when the crater is near the equator.
Fig. 8 Spin component (left panel) and the obliquity component (right panel) of the CYORP coefficient as a function of ∆(<α) and δ, for K = 2.65 Wm^{−1} K^{−1}. Here, d/D_{0} = 0.16 and a = ϵ = π/6. 
6 Application of the CYORP effect on a real asteroid
6.1 Roughness
The surface roughness of asteroids is produced by several processes, including micrometeorite impacts, thermal fatigue, ejecta, or outgassing. It was found that the YORP torque is extremely sensitive to the smallscale surface structures (Statler 2009; Breiter et al. 2009). The microscopic beaming effect of regolith grainsizescale roughness (<1 mm) was shown – using the Hapke reflectance and emissivity model (Breiter & Vokrouhlickỳ 2011) – to have a marginal influence on the YORP effect. The transverse heat conduction across thermal skin depth (~1 cm) causes an asymmetric thermal emission of the east and west sides of a boulder, giving rise to a systematic positive YORP torque (Golubov & Krugly 2012; Golubov & Lipatova 2022). The importance of the various selfmodification effects of a concave feature on the surface was considered gradually and numerical approaches were taken to study it (Statler 2009; Rozitis & Green 2012, 2013a). It was found that the concave feature of surface roughness could dampen the YORP torque by tens of percent. While the pioneering studies by Rozitis & Green (2012) shed light on the effects of roughness, the computational expense and difficulty in studying the functional dependence means that there are severe limitations to the numerical method.
In contrast, the analytical method that we introduce in the present study, and its computational efficiency, simplify the application of roughnessinduced YORP effects to real asteroids or binary systems. Given the objective of our semianalytical method to provide a rapid assessment of the impact of surface roughness, it is particularly wellsuited for models of asteroids with rough convex shapes derived from lightcurve observations. However, when dealing with highresolution shape models, especially those that possess a resolution of a few centimetres (the scale of the thermal skin), a 3D thermal model becomes essential for accurate calculations, owing to the presence of the tangential YORP (TYORP) effect, which accounts for the transverse heat conduction inside a boulder.
For illustrative purposes, to demonstrate the application of our method, we randomly selected the main belt asteroid (272) Antonia as an example. This choice is representative of the majority of asteroids lacking detailed information obtained through in situ observations. We used the shape model obtained from Hanus et al. (2013). To optimise the performance of our model, we assume a relatively high thermal conductivity of 1 Wm^{−1} K^{−1}. We uniformly distributed the roughness across the asteroid’s surface and investigated the total YORP torque (NYORP + CYORP). To do so, we input the normal vector and position vector of surface elements in the shape model to our CYORP model to obtain the CYORP torque coefficient of each surface element. We then used the area of each surface element to calculate the CYORP torque and sum up all CYORP torques generated by all surface elements. The depthtodiameter is assumed to be 0.5, following the assumption made by Rozitis & Green (2012). This has been compared with the sole NYORP torque. The result is depicted in Fig. 10. Our findings confirm that the roughnessinduced CYORP torque damps the normal YORP torque. Specifically, for the asteroid (272) Antonia, the spin component of the torque experiences a damping effect of approximately 35%, while the obliquity component is damped by approximately 64%.
Fig. 9 Obliquity component (right panel) of the CYORP coefficient as a function of the asteroid obliquity e and the crater colatitude α for K = 2.65 Wm^{−1} K^{−1}. Here, d/D_{0} = 0.16 and δ = ∆ = 0. 
6.2 Rotational evolution
The rotational dynamics of asteroids are primarily governed by two key processes: collisions and the YORP effect. The timescale for reorientation of the spin axis resulting from angular momentum transfer during a collision can be expressed as (Farinella et al. 1998) (77)
On the other hand, the typical timescale associated with the YORP effect is approximately given by: (78)
Clearly, the YORP timescale is shorter than τ_{imp,ire–ori}, although its specific value exhibits considerable variation across different asteroids. Consequently, it is widely accepted that the YORP effect primarily governs the rotational evolution of small objects, while collisions play a dominant role in larger objects. The classic static YORP model – which assumes an invariable YORP torque until the asteroid spins up to disruption or spins down to a quasistatic rotation – has been used to study the longterm rotational evolution of asteroids (Rubincam 2000; Pravec et al. 2005; Bottke et al. 2015). A more intricate model, referred to as the ‘stochastic YORP model’ (Bottke et al. 2015), takes into account the resetting of the YORP torque caused by collisions, which arises from the high sensitivity of the YORP effect to surface morphology. Although a suggested timescale of 1 Myr has been proposed for YORP reset (Bottke et al. 2015), a quantitative assessment of the torque changes resulting from craters is yet to be explored. The CYORP effect offers a powerful tool for investigating the stochastic YORP effect. While a comprehensive examination of the stochastic YORP effect is beyond the scope of this paper, we present an example of integrating the CYORP effect into the static YORP model.
In our simulation, a random YORP coefficient is assigned within the range of −0.005 to 0.005, with the coefficient’s sign matching that of the torque. The CYORP torque is introduced specifically when a collision event takes place. The timescale for the impact by an asteroid with the size R_{imp} is (79)
Here, P_{i} = 2.85 × 10^{−18} km^{−2} yr^{−1} is the intrinsic collision probability, C_{R} = 6 × 10^{5}, and b_{R} = 2.2 (Holsapple 2022). In the strength regime formulation, the crater produced by an impactor with the size of R_{imp} has a size of (81)
with Y = 100 Pa and υ_{imp} = 5.3 km s^{−1} .In the gravity regime, (82)
where 𝑔_{ast} is the surface gravity of the asteroid. In each time step (~10^{3} yr), we calculate the numbers of impact craters of different sizes, according to Eqs. (79) and (81). We then assign each crater with a random surface element of the polyhedron model of the asteroid Antonia, after which we can calculate the CYORP torque. Finally, we add the CYORP torque to the normal YORP torque directly calculated from the shape model in order to obtain the total YORP torque. The spin rate evolves following (83)
with I being the maximum moment of inertia and T_{z} being the torque component that is along a spin vector. The obliquity evolves according to (84)
where T_{ϵ} is the torque component that changes the obliquity.
There exist two possible end states in a YORP cycle: either the asteroid’s rotation slows down until it reaches a quasinonrotational state, or it accelerates to the spin threshold for shape change or disruption with a period of approximately 2.2 h.
Upon completing a YORP cycle, the asteroid’s rotational state is updated by assigning a new random rotational speed and obliquity. The impact of introducing CYORP torques can be observed in the evolution of a 10 km asteroid, as depicted in Fig. 11. Notably, significant differences arise when considering the inclusion of CYORP torques.
Nonetheless, the rotational evolution of asteroids currently lacks a standardised model. Some models propose that after spinning down to a nonrotational state, the asteroid’s spin rate is assigned a new random value within a specified range (Hanus et al. 2011; Bottke et al. 2015), while some assume it continues to spin up under the YORP effect (Pravec et al. 2008; Marzari et al. 2020). By selecting an initial spin rate for a new rotational state, Holsapple (2022) reproduces the spin evolution without the YORP effect. Hence, rather than attempting to address the entire complexity of the problem, our objective in this study is to present an illustrative example of the interaction between CYORP and the conventional YORP effect. Furthermore, we underscore the significance of the CYORP effect in the longterm rotational evolution of asteroids. A comprehensive investigation of the rotational evolution of asteroid groups is left for future research.
Fig. 10 YORP torque damped by the CYORP effect in the case of asteroid (272) Antonia. The spin component and the obliquity component are shown in the left and middle panels, respectively. The right panel shows the CYORP coefficient distribution over the asteroid’s surface. 
Fig. 11 Evolution of the spin rate (left panel) and the obliquity (right panel) of a 10 km synthetic asteroid. In the presence of the static YORP torque (blue line), the asteroid gradually decelerates until it reaches a quasinonrotational state. Subsequently, we impose a new rotational state on the asteroid by assigning random values of rotational speed and obliquity. Conversely, when incorporating the CYORP torque (red line), the asteroid follows a different path, exhibiting random fluctuations in its spin rate due to the occurrence of impacts, creating new craters that lead to changes in the CYORP torque. As a result, the asteroid experiences intermittent transitions between spin up and spin down. 
7 Summary and conclusions
The YORP effect is a thermal torque produced by radiation from the irregular surface of the asteroid. It has been demonstrated that this effect is highly sensitive to surface topology (Statler 2009; Breiter et al. 2009), including smallscale roughness (Rozitis & Green 2012), boulders (Golubov & Krugly 2012), and craters (Zhou et al. 2022). In this study, we developed a semianalytical model for calculating the temperature field of a crater, which accounts for the effects of selfsheltering, selfradiation, selfscattering, and nonzero thermal conductivity. Using this model, we investigated the craterinduced YORP (CYORP) effect in a computationally efficient manner (about three orders of magnitude faster than the numerical method), allowing for a comprehensive exploration of the functional dependence of the CYORP effect and its incorporation into the rotational and orbital evolution of asteroids. The main results and conclusions of this study can be summarised as follows.
Our semianalytical model for the CYORP effect is valid in the highthermalconductivity regime (K > 0.1 Wm^{−1} K^{−1}). This suggests that the model is suitable for application to materials such as solid basalt and metal, which are usually beneath the regolith on asteroid surfaces but may be exposed to sunlight due to the formation of deep craters.
The CYORP effect is significant when the crater depthtodiameter ratio is greater than 0.05. The selfmodification effects of a concave structure, including the selfsheltering effect, selfradiation effect, and selfscattering effect, are stronger with a higher depthtodiameter ratio. For concave structures with a depthtodiameter ratio of smaller than 0.05, the surface can be treated as a convex shape without introducing significant inaccuracies. The typical value of the CYORP coefficient for the spin component is 0.01, which is insensitive to the thermal parameter, while the obliquity component decreases from 0.004 for basalt to 0.001 for metal. For a zaxis symmetric shape (e.g. a spinning top shape), the spin component of the CYORP torque vanishes while the obliquity component survives, which implies that the spin acceleration of such symmetric shapes does not change significantly under the effect of crater formation.
Using our semianalytical method, we confirm that the YORP torque can be damped by the surface roughness, which was first discovered by Rozitis & Green (2012). The fast computation of our semianalytical model allows us to consider more flexible configurations of surface roughness, such as a spacevarying roughness distribution, roughness on components of binary asteroids, and so on.
The magnitude and directional change of the YORP coefficient at each impact are assessed for the first time using our CYORP model. While a complete investigation of the spin evolution of asteroids is left for future work, we show that rotational evolution can be severely affected by collisions.
Acknowledgements
We acknowledge support from the Université Côte d’Azur. We thank Yun Zhang, Xiaoran Yan, and Marco Delbo for the useful discussion. WenHan Zhou would like to acknowledge the funding support from the Chinese Scholarship Council (No. 202110320014). Patrick Michel acknowledges funding support from the French space agency CNES and from the European Union’s Horizon 2020 research and innovation program under grant agreement no. 870377 (project NEOMAPP).
References
 Agarwal, J., Kim, Y., Jewitt, D., et al. 2020, A&A, 643, A152 [EDP Sciences] [Google Scholar]
 Bottke, Jr, W. F., Vokrouhlickỳ, D., Rubincam, D. P., & Nesvornỳ, D. 2006, Annu. Rev. Earth Planet. Sci., 34, 157 [CrossRef] [Google Scholar]
 Bottke, Jr, W. F., Vokrouhlickỳ, D., Walsh, K. J., et al. 2015, Icarus, 247, 191 [CrossRef] [Google Scholar]
 Breiter, S., & Vokrouhlickỳ, D. 2011, MNRAS, 410, 2807 [NASA ADS] [CrossRef] [Google Scholar]
 Breiter, S., Bartczak, P., Czekaj, M., Oczujda, B., & Vokrouhlickỳ, D. 2009, A&A, 507, 1073 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Carruba, V., Domingos, R., Huaman, M., Santos, C. d., & Souami, D. 2014, MNRAS, 437, 2279 [NASA ADS] [CrossRef] [Google Scholar]
 Carruba, V., Nesvornỳ, D., Aljbaae, S., & Huaman, M. E. 2015, MNRAS, 451, 244 [NASA ADS] [CrossRef] [Google Scholar]
 Carruba, V., Nesvornỳ, D., & Vokrouhlickỳ, D. 2016, AJ, 151, 164 [NASA ADS] [CrossRef] [Google Scholar]
 Ćuk, M., & Burns, J. A. 2005, Icarus, 176, 418 [Google Scholar]
 Ćuk, M., Christou, A. A., & Hamilton, D. P. 2015, Icarus, 252, 339 [CrossRef] [Google Scholar]
 Delbo, M., Walsh, K., Mueller, M., Harris, A. W., & Howell, E. S. 2011, Icarus, 212, 138 [NASA ADS] [CrossRef] [Google Scholar]
 Ďurech, J., Vokrouhlickỳ, D., Pravec, P., et al. 2022, A&A, 657, A5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Farinella, P., Vokrouhlickỳ, D., & Hartmann, W. K. 1998, Icarus, 132, 378 [NASA ADS] [CrossRef] [Google Scholar]
 Fatka, P., Pravec, P., & Vokrouhlickỳ, D. 2020, Icarus, 338, 113554 [NASA ADS] [CrossRef] [Google Scholar]
 Golubov, O. 2017, AJ, 154, 238 [NASA ADS] [CrossRef] [Google Scholar]
 Golubov, O., & Krugly, Y. N. 2012, ApJ, 752, L11 [NASA ADS] [CrossRef] [Google Scholar]
 Golubov, O., & Lipatova, V. 2022, A&A, 666, A146 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Golubov, O., & Scheeres, D. J. 2019, ApJ, 157, 105 [CrossRef] [Google Scholar]
 Golubov, O., Scheeres, D., & Krugly, Y. N. 2014, ApJ, 794, 22 [NASA ADS] [CrossRef] [Google Scholar]
 Golubov, O., Unukovych, V., & Scheeres, D. J. 2021, AJ, 162, 8 [NASA ADS] [CrossRef] [Google Scholar]
 Hanuš, J., Ďurech, J., Brož, M., et al. 2011, A&A, 530, A134 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hanuš, J., Ďurech, J., Brož, M., et al. 2013, A&A, 551, A67 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Holsapple, K. A. 2022, Planet. Space Sci., 219, 105529 [NASA ADS] [CrossRef] [Google Scholar]
 Jacobson, S. A., & Scheeres, D. J. 2011, Icarus, 214, 161 [CrossRef] [Google Scholar]
 Jacobson, S. A., Scheeres, D. J., & McMahon, J. 2013, ApJ, 780, 60 [NASA ADS] [CrossRef] [Google Scholar]
 Jacobson, S. A., Marzari, F., Rossi, A., & Scheeres, D. J. 2016, Icarus, 277, 381 [NASA ADS] [CrossRef] [Google Scholar]
 Lowry, S. C., Fitzsimmons, A., Pravec, P., et al. 2007, Science, 316, 272 [NASA ADS] [CrossRef] [Google Scholar]
 Lowry, V. C., Vokrouhlickỳ, D., Nesvornỳ, D., & Campins, H. 2020, ApJ, 160, 127 [CrossRef] [Google Scholar]
 Lupishko, D., & Tielieusova, I. 2014, Meteor. Planet. Sci., 49, 80 [NASA ADS] [CrossRef] [Google Scholar]
 Lupishko, D., Mikhalchenko, O., & Chiorny, V. 2019, Solar Syst. Res., 53, 208 [NASA ADS] [CrossRef] [Google Scholar]
 Marzari, F., Rossi, A., Golubov, O., & Scheeres, D. J. 2020, AJ, 160, 128 [NASA ADS] [CrossRef] [Google Scholar]
 Nesvornỳ, D., & Vokrouhlickỳ, D. 2008, A&A, 480, 1 [CrossRef] [EDP Sciences] [Google Scholar]
 Polishook, D. 2014, Icarus, 241, 79 [NASA ADS] [CrossRef] [Google Scholar]
 Pravec, P., Harris, A. W., Scheirich, P., et al. 2005, Icarus, 173, 108 [NASA ADS] [CrossRef] [Google Scholar]
 Pravec, P., Harris, A., Vokrouhlickỳ, D., et al. 2008, Icarus, 197, 497 [NASA ADS] [CrossRef] [Google Scholar]
 Rozitis, B., & Green, S. F. 2012, MNRAS, 423, 367 [NASA ADS] [CrossRef] [Google Scholar]
 Rozitis, B., & Green, S. F. 2013a, MNRAS, 433, 603 [NASA ADS] [CrossRef] [Google Scholar]
 Rozitis, B., & Green, S. F. 2013b, MNRAS, 430, 1376 [NASA ADS] [CrossRef] [Google Scholar]
 Rożek, A., Lowry, S., Nolan, M., et al. 2019, A&A, 631, A149 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Rubincam, D. P. 2000, Icarus, 148, 2 [Google Scholar]
 Sánchez, P., & Scheeres, D. J. 2020, Icarus, 338, 113443 [Google Scholar]
 Scheeres, D. J. 2007, Icarus, 189, 370 [NASA ADS] [CrossRef] [Google Scholar]
 Ševeček, P., Brož, M., Capek, D., & Durech, J. 2015, MNRAS, 450, 2104 [CrossRef] [Google Scholar]
 Statler, T. S. 2009, Icarus, 202, 502 [NASA ADS] [CrossRef] [Google Scholar]
 Steinberg, E. & Sari, R. 2011, AJ, 141, 55 [NASA ADS] [CrossRef] [Google Scholar]
 Taylor, P. A., Margot, J.L., Vokrouhlicky, D., et al. 2007, Science, 316, 274 [NASA ADS] [CrossRef] [Google Scholar]
 Tian, J., Zhao, H.B., & Li, B. 2022, Res. Astron. Astrophys., 22, 125004 [CrossRef] [Google Scholar]
 Veras, D., & Scheeres, D. J. 2020, MNRAS, 492, 2437 [NASA ADS] [CrossRef] [Google Scholar]
 Vokrouhlickỳ, D., & Capek, D. 2002, Icarus, 159, 449 [CrossRef] [Google Scholar]
 Vokrouhlickỳ, D. & Nesvornỳ, D. 2008, AJ, 136, 280 [CrossRef] [Google Scholar]
 Vokrouhlickỳ, D., Milani, A., & Chesley, S. 2000, Icarus, 148, 118 [CrossRef] [Google Scholar]
 Vokrouhlickỳ, D., Nesvornỳ, D., & Bottke, W. F. 2003, Nature, 425, 147 [CrossRef] [Google Scholar]
 Vokrouhlickỳ, D., Brož, M., Bottke, W., Nesvornỳ, D., & Morbidelli, A. 2006, Icarus, 182, 118 [CrossRef] [Google Scholar]
 Vokrouhlicky, D., Bottke, W. F., Chesley, S. R., Scheeres, D. J., & Statler, T. S. 2015, Asteroids IV, eds. P. Michel, F. E. DeMeo, & W. F. Bottke (Tucson: University of Arizona Press), 509 [Google Scholar]
 Zegmott, T. J., Lowry, S., Rozek, A., et al. 2021, MNRAS, 507, 4914 [NASA ADS] [CrossRef] [Google Scholar]
 Zhou, W.H., Zhang, Y., Yan, X., & Michel, P. 2022, A&A, 668, A70 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
See the definition by Vokrouhlickỳ & Čapek (2002).
All Tables
Thermal parameter for three typical materials on asteroids, taken from Farinella et al. (1998).
All Figures
Fig. 1 Three coordinate systems in this paper: coordinate system oxyz for calculating the illuminated domain in the crater, PABC for calculating the effective radiative force of an arbitrary surface element, and OXYZ for averaging the YORP torque over the spin and orbital motion. 

In the text 
Fig. 2 Radiative force of the total crater averaged over a rotational period (8 h by default) obtained from the analytical method (solid lines) and the numerical method (dashed lines), as a function of the crater colatitude α for K = 2.65 Wm^{−1} K^{−1} (left panel) and K = 0.1 Wm^{−1} K^{−1} (right panel). The x, y, and 𝓏 components of the radiative force are shown in blue, red, and green, respectively. The distance of the crater here is 1 au from the Sun. 

In the text 
Fig. 3 Same as Fig. 2 but for K = 0 Wm^{−1} K^{−1} in the analytical model and K = 0.001 Wm^{−1} K^{−1} in the numerical model. 

In the text 
Fig. 4 Ratio of the radiative force including selfmodification effects to that without selfmodification effects, as a function of crater depthtodiameter ratio and accounting for different selfmodification effects. 

In the text 
Fig. 5 Spin component (left panel) and the obliquity component (right panel) of the CYORP coefficient as a function of the asteroid obliquity e, accounting for different depthtodiameter ratios, which are denoted by different colours. 

In the text 
Fig. 6 Spin component (left panel) and the obliquity component (right panel) of the CYORP coefficient as a function of the asteroid obliquity e and the crater colatitude α(> ∆), for K = 2.65Wm^{−1} K^{−1}. Here, d/D_{0} = 0.16 and δ = ∆ = π/4. 

In the text 
Fig. 7 Same as Fig. 6 but for K = 40 Wm^{−1} K^{−1}. 

In the text 
Fig. 8 Spin component (left panel) and the obliquity component (right panel) of the CYORP coefficient as a function of ∆(<α) and δ, for K = 2.65 Wm^{−1} K^{−1}. Here, d/D_{0} = 0.16 and a = ϵ = π/6. 

In the text 
Fig. 9 Obliquity component (right panel) of the CYORP coefficient as a function of the asteroid obliquity e and the crater colatitude α for K = 2.65 Wm^{−1} K^{−1}. Here, d/D_{0} = 0.16 and δ = ∆ = 0. 

In the text 
Fig. 10 YORP torque damped by the CYORP effect in the case of asteroid (272) Antonia. The spin component and the obliquity component are shown in the left and middle panels, respectively. The right panel shows the CYORP coefficient distribution over the asteroid’s surface. 

In the text 
Fig. 11 Evolution of the spin rate (left panel) and the obliquity (right panel) of a 10 km synthetic asteroid. In the presence of the static YORP torque (blue line), the asteroid gradually decelerates until it reaches a quasinonrotational state. Subsequently, we impose a new rotational state on the asteroid by assigning random values of rotational speed and obliquity. Conversely, when incorporating the CYORP torque (red line), the asteroid follows a different path, exhibiting random fluctuations in its spin rate due to the occurrence of impacts, creating new craters that lead to changes in the CYORP torque. As a result, the asteroid experiences intermittent transitions between spin up and spin down. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.