MARTIANS (MARs2020, TIANwen and So on) would see more potentially hazardous asteroids than Earthlings

Potentially Hazardous Asteroids (PHAs) are a special subset of Near-Earth Objects (NEOs) that can come close to the Earth and are large enough to cause significant damage in the event of an impact. Observations and researches of Earth-PHAs have been underway for decades. Here, we extend the concept of PHAs to Mars and study the feasibility of detecting Mars-PHAs in the near future. We focus on PHAs that truly undergo close approaches with a planet (dubbed CAPHAs) and aim to compare the actual quantities of Earth-CAPHAs and Mars-CAPHAs by conducting numerical simulations incorporating the Yarkovsky effect, based on observed data of the main asteroid belt. The estimated number of Earth-CAPHAs and Mars-CAPHAs are 4675 and 16910, respectively. The occurrence frequency of Mars-CAPHAs is about 52 per year, which is 2.6 times that of Earth-CAPHAs, indicating significant potential for future Mars-based observations. Furthermore, a few Mars-CAPHAs are predicted to be observable even from Earth around the time of next Mars opposition in 2025.


INTRODUCTION
Near-Earth objects (NEOs) are small bodies whose perihelion distances  ≤ 1.3 au and aphelion distances  ≥ 0.983 au, including near-Earth asteroids (NEAs) and near-Earth comet nuclei (NECs).Potentially hazardous asteroids (PHAs), a special subset of NEAs, are defined as having a minimum orbit intersection distance (MOID) to Earth orbit of less than 0.05 au and absolute magnitude  ≤ 22 mag (i.e., diameter  ≥ 140 m).PHAs are more likely to collide with Earth than other NEAs, and if such a collision were to occur, their size is large enough to cause destructive effects (for a review see Perna et al. 2013).However, PHAs are not necessarily highly 'dangerous' because orbital proximity does not imply proximity in actual distance (e.g. the case of Earth Trojan asteroids).Therefore, only those PHAs that truly undergo close approaches (hereafter called CAPHAs) to the Earth are the primary objects of human monitoring 1 , in which case the asteroid has a minimum distance  min < 0.05 au from Earth.
NEAs have been objects of human attention for decades, motivated by concerns for both scientific interests and our own safety.On one hand, such objects could have delivered water and organic-rich materials to the early Earth, which is crucial for our understanding about the origin of terrestrial life as well as the formation and early evolution of the Solar System (e.g.Marty 2012;Alexander et al. 2012;Altwegg et al. 2015).On the other hand, potential impact events can pose significant threats to the Earth.The most famous ★ E-mail: yufanz@smail.nju.edu.cn† E-mail: lizy@nju.edu.cn 1 https://cneos.jpl.nasa.gov/caone is undoubtedly the so-called 'K-T' event leading to the Cretaceous/Tertiary extinction that occurred about 65 million years ago (Kyte 1998;Bottke et al. 2007).The Tungus Explosion in 1908 (Ben-Menahem 1975) and the Chelyabinsk meteor falling in 2013 (Brown et al. 2013) are also well-known examples.
According to the Small-Body Database 2 (SBDB) of Jet Propulsion Laboratory (JPL), as of June 26, 2023, over 30000 NEAs have been observed.The dynamical lifetime of NEAs is typically between 10 6 and 10 8 years (Morbidelli & Gladman 1998), so asteroids must be continually transported to maintain the population we observe around Earth today (Mainzer et al. 2011b).Studies by Wetherill (1988), Rabinowitz (1997a,b), Bottke et al. (2002) suggest that NEAs are probably sourced from unstable regions in the main asteroid belt, such as the 3:1 mean motion resonance with Jupiter (Wisdom 1983;Moons 1996) and the  6 secular resonance (Morbidelli et al. 1994), which manifest as 'gaps' in the orbital distribution of asteroids.However, precisely because such regions are unstable, the number of asteroids in these locations is relatively low.Therefore, a mechanism that can push nearby asteroids into these resonance regions must be considered: the Yarkovsky effect.
The Yarkovsky effect (Burns et al. 1979;Bottke et al. 2006) was first discovered by the engineer Ivan Osipovich Yarkovsky in 1901.About a century later, it was shown to have an important influence on the motion of artificial satellites (Rubincam 1987(Rubincam , 1995) ) and small bodies (Farinella et al. 1998).The Yarkovsky effect is caused by the asymmetric thermal emission (with respect to the subsolar point) 2 https://ssd.jpl.nasa.gov/tools/sbdb_query.htmlfrom a rotating, atmosphereless body with finite thermal inertia, and it manifests itself as a recoil force.Radiation from the Sun heats up asteroids, with the side facing the Sun having the highest temperature.After a short period of time, the asteroid will re-radiate its thermal energy and generate a net radiation imbalance.As the asteroid is rotating, its hottest side may no longer be perfectly aligned with the Sun.Therefore, the recoil force has a transverse component, which will accelerate or decelerate the asteroid's orbital motion, thus change the semi-major axis of the asteroid over a reasonable long time scale (Peterson 1976;Farinella et al. 1998;Farinella & Vokrouhlicky 1999).
The above concepts are also applicable to asteroids near Mars.However, numerical simulations and observations of such objects are currently very limited (e.g.Migliorini et al. 1998;Michel et al. 2000;Alí-Lagoa & Delbo' 2017), and the focus has mainly been on Mars-crossing asteroids (MCAs), which are defined as asteroids that can cross the orbit of Mars.Here we focus on asteroids that truly close to Mars and propose the concept of Mars-CAPHAs.For Mars-CAPHAs, we set the distance criterion as  min < 0.036 au, instead of 0.05 au, to take into account the difference in the Hill radius between Earth and Mars (choosing 0.05 au for Mars is also acceptable and will result in more Mars-CAPHAs).
The study of Mars-CAPHAs will not be merely conceptual.Indeed, humanity possesses the capability to investigate them in situ.Since the 1960s, over forty Mars exploration missions have been conducted.With the advancement of technology, missions have progressed from early attempts at flybys, such as Mariner-4 (Huntress et al. 2003), to include orbiters, landers and rovers, such as 2001 Mars Odyssey (Saunders et al. 2004), Mars Express (Nielsen 2004) and Zhu Rong (Li et al. 2021).In the future, projects such as Mars 2020 from the USA (Perseverance has already arrived on Mars) (Farley et al. 2020), Tianwen-3 from China, ExoMars from Europe, among others, will utilize orbiters and rovers to carry out more in-depth explorations of Mars.Observing and imaging Mars-CAPHAs with Mars-based instruments, although challenging, may bring us new insights into the Martian environment.
Therefore, numerical simulations about Mars-CAPHAs are timely and necessary to assess their observational prospects.We conduct an N-body simulation incorporating the Yarkovsky effect for such a study.The methods are presented in section 2. We introduce the software and the implementation of the Yarkovsky effect.Information about the data used in this work and our simulation setup is also provided.In section 3, we estimate some observable properties of CAPHAs and address their implications.A brief summary is presented in section 4.

Mercury6 with the Yarkovsky force
The numerical simulation tool we use is the N-body software package Mercury6 (Chambers 1999), which can integrate any star-planet system like our Solar System, with an additional Yarkovsky force subroutine.The hybrid integrator is utilized, which is basically a symplectic integrator switching to Bulirsch-Stoer (an algorithm with adaptive timesteps) only when planetary close encounters occur.We incorporated the Yarkovsky force into the subroutine mfo_user using the same approach as Zhou et al. (2019).
The Yarkovsky effect results in a recoil force on asteroids, and the precise calculation of this force is highly complex and dependent on lots of physical parameters, many of which have not yet been accurately measured.In some studies (e.g.Vokrouhlicky 1998;Vokrouhlický 1999), the Yarkovsky effect is simply described as an equivalent drift rate of the semi-major axis d/d of an asteroid.Xu et al. (2020b) calculated that an asteroid located at  = 2.9 au with diameter  = 1 km has a drift rate of d/d = 0.256 au Gyr −1 , assuming rotation period  = 8 h, albedo  = 0.13, thermal conductivity  = 0.005 W m −1 K −1 , specific heat capacity  = 680 J kg −1 K −1 , surface density  s = 1.5 g cm −3 , bulk density  = 2.5 g cm −3 and spinning obliquity  = 0 deg (Vokrouhlický et al. 2006).Although Xu et al. (2020a) indicates that the relationship between the drift rate and semi-major axis may be complex, d/d can be considered inversely proportional to  2 within the range of the main belt (e.g.Nesvorný et al. 2023).Additionally, the drift rate is inversely proportional to the diameter .
According to Xu et al. (2020b), for an asteroid with diameter  = 1 km in the 7:3 mean motion resonance (2.958 au), we set its drift rate to d/d = 0.25 au Gyr −1 .Then following the same rules as Nesvorný et al. (2023), the formula for drift rate used in our work is (1)

Asteroid data sampling and grouping
Earth-CAPHAs and Mars-CAPHAs mainly originate from the main belt (e.g.Rabinowitz 1997a,b), which has been stable for a considerable period of time.During this time, the escape of asteroids from the main belt does not significantly change its distribution.Based on this assumption, we use observed data of main belt asteroids as the initial conditions for our simulations.From the SBDB of JPL, as of June 26, 2023, we downloaded the data containing a total of 1,223,263 asteroids.Fig. 1 shows the scatter plot of inclination  versus semi-major axis  of the main belt based on the 1,223,263 asteroids.Several strong mean motion resonances and the  6 secular resonance are marked.After removing those with large measurement errors ( a ≥ 0.0001 au), 1,144,461 asteroids remain.For a feasible computational cost, we randomly sampled 1% (11445 asteroids) as our 'sample data'.
We divide the 11445 asteroids into two categories named 'neargap' asteroids (709) and 'far-gap' asteroids (10736) (specific ranges of 'near-gap' and 'far-gap' regions can be found in column 3 of Ta-ble 1).Since different mean motion resonances have different widths, sizes of our six 'near-gap' regions should be set differently.We set the range of each 'near-gap' region to be significantly wider than its resonance width and to contain a sufficient number of asteroids.For different mean motion resonances, the lower its order is, the larger width we assign to it.It is worth mentioning that the 4:1 resonance is very close to the  6 resonance, so we set a relatively large width for it and extend the right edge to 2.200 au (Wetherill 1985) to embrace the  6 resonance (e.g.Scholl & Froeschle 1991).For other relatively strong resonances not shown in Fig. 1, such as the 5:1 resonance (1.780 au), even setting a particularly large width does not include more than 10 asteroids in our sample, so we do not consider them while studying 'near-gap' asteroids.

Simulation settings
The dynamical model used in our simulations consists of the Sun, eight planets from Mercury to Neptune and (massless) asteroids.The Earth-Moon barycentre is adopted instead of the separate Earth and Moon (e.g.Dvorak et al. 2012;Zhou et al. 2019).The initial data of eight planets at epoch JD 2457724.5 are obtained from the JPL HORIZONS system (Giorgini et al. 1996).Asteroids with different epochs from planets are integrated to the same epoch before the main integration begins.The initial timestep and accuracy parameter we set are 8 days and 1e-12, respectively.Our timestep is a bit longer compared to other work (Granvik et al. 2018), but here we focus on the statistical level results rather than high precision orbital integration of any specific body.In our simulations, an asteroid will no longer be tracked after being ejected from the Solar System (defined as being at a distance greater than 100 au from the Sun).
For 'near-gap' asteroids, we perform simulations with the Yarkovsky effect for 0.1 Gyr, which is sufficiently long for the effect to become evident.For 'far-gap' asteroids, we neglect the Yarkovsky effect (i.e.only with gravity) and run the simulations for 1 Myr only, a simplification in order to save computational effort.Granvik et al. (2017) showed that the Yarkovsky effect can deliver 100-metrediameter asteroids from most regions in the main belt to  = 1.3 au within tens of millions of years, but such a delivery predicts too many NEAs compared to observations.On the other hand, the combined Yarkovsky and YORP (Bottke et al. 2006) model yields results lower than the observed values, possibly due to the need for improvement in the canonical YORP model (Granvik et al. 2017).Therefore, currently there are some issues with the drift timescale of main belt asteroids caused by the Yarkovsky (and YORP) effect.Nevertheless, for 'near-gap' asteroids, the drift timescale does not matter much, as they are initially close to strong resonances, hence considering only the Yarkovsky effect can yield results consistent with observations (e.g.Xu et al. 2020a).Moreover, since we are mostly concerned with asteroids with sizes greater than 100-metre, which would take an even longer time to excite, 'far-gap' asteroids should contribute only a minor fraction of the resultant CAPHAs.Thus, our separate treatment for the 'near-gap' and 'far-gap' asteroids is justified.
When calculating the Yarkovsky drift rate of an asteroid, its semimajor axis  and diameter  are required.Here, for a 'near-gap' region,  of all asteroids within it are approximated by the value of the center of the region.However,  of most asteroids are still unknown.Additionally, observations of the main belt are incomplete: asteroids with small  (i.e., large absolute magnitude ) must have been missed so far, and such selection effect becomes more severe for more distant regions.Therefore, we generate multiple sets of cloned asteroids for each 'near-gap' region, assigning them different  values in the simulations, in order to prepare for subsequent correction for the selection effect. and  follow the empirical conversion formula: Here albedo  is set to 0.14 (Mainzer et al. 2011a), which is close to the value used in the calculation of the drift rate by Xu et al. (2020b).
The absolute magnitude range we consider is from 15 mag (the number of asteroids with  < 15 mag is small) to 22 mag (the upper limit of CAPHA's definition).We divide the range into 7 intervals, with each interval spanning 1 mag.For each interval, we approximate the entire range using the central magnitude value.
The specific settings for our simulations are shown in columns 1-3 of Table 1, where six 'near-gap' regions are designated by A-F and the 'far-gap' region is designated by G.For each 'near-gap' region, a total of 8 simulations are conducted.Simulation 0 serves as the control simulation without the Yarkovsky effect, while simulations 1-7 correspond to the  value ranging from 15.5 to 21.5.
Furthermore, in reality, depending on the direction of rotation of an asteroid, the Yarkovsky effect can either increase or decrease its semi-major axis.However, rotation directions of most asteroids are unknown.Therefore, for each 'near-gap' region, we set the direction of the Yarkovsky force such that the asteroid is drifted toward the center of the gap.Such a treatment artificially raises the number of CAPHAs by a factor of two, which is subsequently corrected in section 3. A related caveat is that this treatment prevents asteroids from drifting into more distant resonances in the opposite direction, but the fractional contribution of CAPHAs through this channel is expected to be negligible.

Quantities of CAPHAs
The direct outputs from our simulations are presented in columns 4-5 of Table 1, which are subject to further processing.First, for simulations 1-7, we calculate the increments Δ relative to simulation 0. Considering that our setting for the direction of the Yarkovsky force would result in the increments being doubling the actual values, we correct the results of simulations 1-7 to  * =  0 + 1 2 Δ, where  0 is the number of CAPHAs predicted by simulation 0.
Second, we introduce the active coefficient where  sim is the simulation duration and  act is the active lifetime.The active lifetime is the time difference between an asteroid's first and last occurrence as a CAPHA; for asteroids with just a single occurrence as a CAPHA, the active lifetime is set to be 5 years, which corresponds to the orbital period for  = 3 au (the center of our studied range).We introduce this coefficient because the quantities obtained from our simulations do not represent the current observable numbers.For example, if one CAPHA, whose active lifetime is 1 Myr, is generated in a simulation spanning 100 Myr, we consider its current survival probability to be 0.01, so its contribution to the total count is not 1 but rather 0.01.Average  values for different simulations are presented in columns 10-11 of Table 1.Third, we introduce the supplement coefficient () to correct the selection effect of the main belt.Considering the different completeness at different distances, we partition the main belt into five zones named , , ,  and , using 2.3 au, 2.7 au, 3.1 au and 3.5 au as dividing points.Our 'near-gap' simulations A, B, C&D, E, F respectively correspond to zone , , ,  and .Assuming that the actual Table 1.Simulation ID, Yarkovsky drift rate, region information (including the radial range, corresponding resonance and number of asteroids in our sample), number of CAPHAs ( ), difference caused by the Yarkovsky effect (Δ ), reasonable number ( * =  0 + 1 2 Δ ), average active coefficient ( ε ), supplement coefficient () and currently observable number ( obs ) of our simulations.Series A-F belong to 'near-gap' asteroids and correspond to six mean motion resonances (4:1, 3:1, 5:2, 7:3, 2:1, 3:2).Simulation 0 in each series corresponds to the control group, while simulations 1-7 correspond to the cases where the absolute magnitudes (diameters) are 15.5 mag (2.82 km), 16.5 mag (1.78 km), 17.5 mag (1.12 km), 18.5 mag (0.71 km), 19.5 mag (0.45 km), 20.5 mag (0.28 km) and 21.5 mag (0.18 km), respectively.Series G corresponds to 'far-gap' asteroids.Zero drift rate means no Yarkovsky effect.Except for series G with a duration of 0.001 Gyr, all other simulations have a duration of 0.1 Gyr.'EC' and 'MC' are abbreviations for Earth-CAPHA and Mars-CAPHA.The two  values for series G correspond to the values of Earth-CAPHAs and Mars-CAPHAs, respectively.The range (au) of series F (3:2 resonance) is 3.972 ± 0.050, whose result is not shown as it contributes insignificantly to the population of CAPHAs.cumulative magnitude distribution of asteroids follows an exponential function (e.g.Gladman et al. 2009), we conduct fittings using asteroids with 15 ≤  ≤  ′ in each of the five zones.The values of  ′ are chosen based on the observed cumulative magnitude distribution and are respectively set to 18.0, 18.0, 17.3, 17.0 and 16.2.We obtain:

ID
Furthermore, there is strong evidence showing that the sizefrequency distribution of the main belt flattens for  > 18 (e.g.Gladman et al. 2009), so we assume that the slopes of the five subequations are halved when  reaches 18.Therefore, for each zone, we can deduce the ratio of the number of asteroids in 7 magnitude intervals (corresponding to simulations 1-7), and we define  xk as the -th value among the 7 ratios of zone .For example, the ratio for zone  is 1. 00 : 2.73 : 7.46 : 7.68 : 12.70 : 20.98 : 34.68,where the first value ( iv1 , set to 1.00) corresponds to 15 mag-16 mag, and the last value ( iv7 ) corresponds to 21 mag-22 mag.However, due to the selection effect against faint objects, the observed distributions significantly differ from theoretical ones.Therefore, we introduce the supplement coefficient for zone , magnitude interval : x,1 and  x,total represent the number of observed asteroids with 15 ≤  < 16 and 15 ≤  < 22 in zone , respectively.For a 'near-gap' simulation (A-F), determining its associated zone and magnitude interval is sufficient to calculate .For the 'far-gap' simulation (G), it spans across multiple zones.Hence, starting from the outputs, we calculate the average value for zones belonged to each recorded CAPHA.Therefore,  of Earth-CAPHAs and Mars-CAPHAs for simulation G are different.Then if one simulation contributes  CAPHAs, we correct the quantity to .The  values of all simulations are shown in column 12 of Table 1.Lastly, we multiply the quantity  * by coefficients  and  and sum up results of all simulations to estimate the observable number of CAPHAs: We obtain a total count of 46.75 Earth-CAPHAs and 169.10 Mars-CAPHAs.Considering that our sample represents only 1% of the actual main belt, the predicted total numbers of observable Earth-CAPHAs and Mars-CAPHAs are 4675 and 16910, i.e., a ratio of approximately 1:3.6 (see Fig. 2).Then we estimate their expected occurrence frequencies (in units of yr −1 ).We first calculate the average frequency of each CAPHA for each simulation: Xk represents the number of CAPHAs generated in simulation Xk (e.g., A1),  denotes the repeating occurrence count of a CAPHA, and  act is its active lifetime.Next, for equation ( 6), we multiply each term in the summation by its corresponding f , yielding the total occurrence frequency: The resultant occurrence frequencies of Earth-CAPHAs and Mars-CAPHAs are 20 and 52 per year, respectively, making the latter approximately 2.6 times higher.It is noteworthy that the exact numbers and frequencies depend on the adopted intrinsic size distribution of the main belt asteroids given by equation ( 4), but the ratio of either quantity between Mars-CAPHAs and Earth-CAPHAs is insensitive to it.
Over 1300 independent (counted without repetition) Earth-CAPHAs have been observed and catalogued as of June 26, 2023, indicating that there are still at least 2.5 times more left to be discovered.The simulation predicted occurrence frequency of Earth-CAPHAs (20 per year) is in rough agreement with the observed frequency, which is about 27 per year averaged over 20 years, in view of the potential underestimation of the contribution from the 'far-gap' region (Section 2.3).Both in terms of number and frequency, Mars-CAPHAs significantly surpass Earth-CAPHAs, indicating tremendous observational prospects.
It is worth noting that a Mars-CAPHA may become an Earth-CAPHA later, and vise verse.For example, in simulation A2, a total of 72 Earth-CAPHAs and 125 Mars-CAPHAs were recorded, with 69 of them being in common.This raises an intriguing possibility of using CAPHAs as matter carriers between Earth and Mars.

Distances and velocities
We have further investigated the distance and velocity distributions of CAPHAs.Considering that the maximum distances  crit from the central planet to Earth-CAPHAs and Mars-CAPHAs are not the same (0.05 au for Earth and 0.036 au for Mars, see section 1), we focus on the relative distance, which is defined as the ratio of a CAPHA's close-approach distance  min and its corresponding planet's  crit .Fig. 2 shows that Earth-CAPHAs and Mars-CAPHAs exhibit a similar distribution in terms of relative distances, and the actual observed results for Earth-CAPHAs align well with this distribution.Certainly, in terms of absolute distance, the majority of Mars-CAPHAs are closer to the planet than Earth-CAPHAs, making them more favorable for observations.Moreover, the relative velocity (which remains nearly constant during a close approach of a given asteroid) of Mars-CAPHAs is systematically smaller compared to Earth-CAPHAs (see Fig. 2).It can be explained by the fact that, for the same orbit, when an asteroid is near Earth, it is closer to its perihelion, resulting in a higher velocity.This implies that the observing window is longer for Marsbased observations, or in other words, tracking CAPHAs from Mars is easier.The actually observed distribution of relative velocities for Earth-CAPHAs is also shown in Fig. 2, which matches our simulations well, demonstrating the validity of our simulations.

Active lifetime of different CAPHAs
The longer the active lifetime (defined in section 3.1) of a CAPHA, the higher the probability that it will coexist with the civilization on a particular planet, thereby increasing the chances of its observation.Our simulations reveal two main aspects of the active lifetime of CAPHAs.
First, the average active lifetime of Mars-CAPHAs (the longest reaching several tens of million years) is longer than that of Earth-CAPHAs (mostly ranges from several hundred thousand to several million years; see  values in Table 1), which can be explained by the difference in their eccentricities.Since Earth is farther from the main belt compared to Mars, an asteroid requires a higher eccentricity to become an Earth-CAPHA.A higher eccentricity implies a greater likelihood of close encounters with planets such as Jupiter (e.g.Granvik et al. 2018), which will kick asteroids out from the inner solar system, thus resulting in a shorter active lifetime.
Second, as the source region in the main belt becomes more distant, the average active lifetime of CAPHAs for a specific planet (Earth or Mars) decreases initially and then slightly increases, with a turning point occurring at around 2.96 au (see Table 1).Such a trend is determined by two factors: eccentricity and period.On one hand, the farther the source region is, the higher the eccentricity would be required for an asteroid to become a CAPHA, resulting in a shorter active lifetime.On the other hand, asteroids with larger semimajor axes have longer orbital periods, resulting in a lower frequency of close encounters with planets, and consequently, a longer active lifetime.For the nearer source regions, the first factor prevails, while as the distance increases, the second factor becomes more significant.

Evaluating the feasibility of Earth-based observations
In order to investigate the feasibility of Earth-based observations for Mars-CAPHAs, it is necessary to know their apparent magnitudes.Apparent magnitude  (, ) of an asteroid (or a planet) is determined by its absolute magnitude , distance  and phase angle , and can be expressed simply as:  (, ) =  +  (, ). (9) Due to the close proximity of Mars-CAPHAs and Mars, their  and  can be taken as the same.As a result, they follow the same apparent magnitude conversion rules.The next Mars opposition (corresponding to  = 0) will occur on January 16, 2025, so we focus on the time period around that date.Between December 10, 2024 and February 13, 2025, the apparent magnitude of Mars will be brighter than -0.52 mag, according to the ephemeris.As the absolute magnitude of Mars is -1.52 mag, the value of  (, ) during that period must be less than 1.Therefore, for Mars-CAPHAs appearing at that time, their apparent magnitude will not exceed their absolute magnitude plus one.Assuming the limiting magnitude of Earth-based observing facility is  lim = 20 mag, we need to consider Mars-CAPHAs with  ≤ 19 mag.Based on our simulations, the observable frequency of Mars-CAPHAs is at least 52 per year.According to equation (4), more than 22% of Mars-CAPHAs have an absolute magnitude brighter than 19 mag, which corresponds to over 11 occurrences per year.Therefore, during the aforementioned 65-day period, at least 2 Mars-CAPHAs can be observed from Earth.At that time, Earth-based telescopes, such as LSST (Ivezić et al. 2019) and WFST (Wang et al. 2023), may have the best chance to detect Mars-CAPHAs, which can provide guidance for the planning of Mars-based observations.

SUMMARY
We utilize the N-body software package Mercury6 with an additional Yarkovsky effect subroutine to investigate Earth-CAPHAs and Mars-CAPHAs.Based on the assumption that CAPHAs originate from the main belt, we randomly sample 1% of main belt asteroids as the initial conditions for our simulations.Considering special properties of gaps in the main belt, we divide our sample into 'far-gap' asteroids and 'near-gap' asteroids, and run them with different configurations.
The number and occurrence frequency of Mars-CAPHAs are a few times greater than that of Earth-CAPHAs, and the relative velocities of Mars-CAPHAs are favorable for Mars-based observations.Our results point to a new potential scientific objective for future Mars exploration missions (Mars 2020, Tianwen and so on).Observations of Mars-CAPHAs will deepen our understanding of the Martian environment, the interactions between asteroids and planets, and the evolutionary history of the inner Solar System.Lastly, a few Mars-CAPHAs are predicted to be observable even from Earth in 2025.

Figure 1 .
Figure 1.Scatter plot of inclination  versus semi-major axis  of the main belt (1,223,263 asteroids).The color reflects the number density of the scatter points.Several mean motion resonances and the  6 secular resonance are marked with reddish purple dashed lines.

Figure 2 .
Figure 2. Upper left (right) panel: Numbers (Frequencies) of Earth-CAPHAs and Mars-CAPHAs obtained from simulations and actual observations.Lower left panel: Distribution of relative distances between CAPHAs and the planet.The relative distance of a CAPHA in percentage is defined as the ratio between its  min and  crit (0.05 au for Earth-CAPHAs and 0.036 au for Mars-CAPHAs).Blue and reddish purple lines almost coincide with each other.The yellow histogram shows the actually observed distribution of Earth-CAPHAs.Lower right panel: Distribution of velocities of CAPHAs relative to the planet.The yellow histogram shows the actually observed distribution of Earth-CAPHAs.