Search

The rotational distribution of asteroids as a function of their size is used {as a diagnostic of} their physical properties and evolution. Recent photometric surveys from the Gaia mission, allowing observation of asteroids with long spin periods (for example $\geq 24$h), found an excessive group of slow rotators and a gap separating them from faster rotators, which is unexplained by current theories. Here we developed an asteroid rotational evolution model capable of reproducing the observed distribution. {We suggest that this distribution is regulated by the competition between collisions and internal friction dampening of "tumblers" -asteroids with unstable rotation vectors, and that the slow rotator group is mainly populated by tumblers.} {We constrain the product of the rigidity and quality factor, which relates to the body's viscosity, to $\mu Q \sim 4 \times 10^9~$Pa. This number, two orders of magnitude smaller than the one assumed for monolithic boulders,} implies that {rubble pile} asteroids could have a porous structure or a thick regolith layer, and undergo stronger tidal effects., Comment: Published in Nature Astronomy

The binary Yarkovsky effect on the secondary asteroid (BYS) was recently discovered to influence binary asteroid systems by pushing the secondary asteroid toward a synchronous orbit on a short timescale. However, the binary Yarkovsky effect on the primary (BYP) remains less understood, partly due to non-linear effects from partial eclipses, but could have significant implications for singly synchronous binaries. In this work, we studied the BYP effect by numerical methods and estimated its induced orbital drifting rates for real binary asteroids. We find an empirical modified solution to estimate the effective BYP: the traditional BYP formula multiplied by $(r_s / r_p)^(\alpha -1 )$. We confirm that the BYP pushes the primary towards a synchronous orbit where its spin equals the mean motion. The parameter $\alpha$ is insensitive to the ratio of the spin rate to the mean motion and decreases slightly with increasing thermal inertia. For small binary systems with a typical thermal inertia of 200 tiu, $\alpha$ is approximately 1.7. The BYP is found to affect the mutual orbit of singly synchronous binaries with a timescale typically an order of magnitude longer than that of the BYS. Drift rates induced by the BYP for known small binary asteroids (primary radius < 1 km) range from -0.001 to -1 cm $yr^{-1}$. A comparative analysis with observed orbital drift rates shows agreement for pre-impact Didymos and 1996 FG$_3$ but discrepancies for 2001 SL$_9$ and 1999 KW$_4$, suggesting complex dynamics in these systems involving the BYP, the binary Yarkovsky-O'Keefe-Radzievskii-Paddack (BYORP) effect, and tides. The BYP is changing the mutual orbits of most discovered binary asteroids. We suggest that the BYP should be considered along with BYORP and tidal effects when studying binary systems' long-term dynamics., Comment: 7 pages, 5 figures. Published in A&A Letters

Context. About 25% -- 50% of white dwarfs are found to be contaminated by heavy elements, which are believed to originate from external sources such as planetary materials. Elemental abundances suggest that most of the pollutants are rocky objects and only a small fraction of white dwarfs bear traces of volatile accretion. Aims. In order to account for the scarcity of volatile pollution, we investigate the role of the white dwarfs' magnetospheres in shielding the volatile content of icy objects. Methods. We estimated the volatile sublimation of inward-drifting exocomets. We assume the orbits of the exocomets are circularized by the Alfven wing drag that is effective for long-period comets. Results. Volatile material can sublimate outside the corotation radius and be shielded by the magnetic field. {The two conditions for this volatile-shielded mechanism are that the magnetosphere radius must be larger than the corotation radius and that the volatiles are depleted outside the corotation radius, which requires a sufficiently slow orbital circularization process.} We applied our model to nine white dwarfs with known rotational periods, magnetic fields, and atmosphere compositions. Our volatile-shielded model may explain the excess of volatile elements such as C and S in the disk relative to the white dwarf atmosphere in WD2326+049 (G29-38). Nevertheless, given the sensitivity of our model to the circularization process and material properties of icy objects, there remains considerable uncertainty in our results. Conclusions. Our work suggests a possible explanation for the scarcity of volatile-accretion signatures among white dwarfs. We also identify a correlation between the magnetic field strength, the spin period, and the composition of pollutants in white dwarf atmospheres., Comment: 13 pages, 11 figures. Accepted for publication in A&A

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 small-scale roughness, boulders, and craters. Aims: The aim of this paper is to develop a computationally efficient semi-analytical model for the crater-induced YORP (CYORP) effect that can be used to investigate the functional dependence of this effect. Methods. This study obtains the temperature field in a crater over a rotational period in the form of a Fourier series, accounting for the effects of self-sheltering, self-radiation, and self-scattering. 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 depth-to-diameter 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. For a z-axis 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 small-scale 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. The semi-analytical 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., Comment: 13 pages, 11 figures. Published on A&A in February 2024

We explore the Yarkovsky effect on small binary asteroids. While significant attention has been given to the binary YORP effect, the Yarkovsky effect is often overlooked. We develop an analytical model for the binary Yarkovsky effect, considering both the Yarkovsky-Schach and planetary Yarkovsky components, and verify it against thermophysical numerical simulations. We find that the Yarkovsky force could change the mutual orbit when the asteroid's spin period is unequal to the orbital period. Our analysis predicts new evolutionary paths for binaries. For a prograde asynchronous secondary, the Yarkovsky force will migrate the satellite towards the location of the synchronous orbit on ~100 kyr timescales, which could be faster than other synchronization processes such as YORP and tides. For retrograde secondaries, the Yarkovsky force always migrates the secondary outwards, which could produce asteroid pairs with opposite spin poles. Satellites spinning faster than the Roche limit orbit period (e.g. from ~4h to ~10h) will migrate inwards until they disrupt, reshape, or form a contact binary. We also predict a short-lived equilibrium state for asynchronous secondaries where the Yarkovsky force is balanced by tides. We provide calculations of the Yarkovsky-induced drift rate for known asynchronous binaries. If the NASA DART impact broke Dimorphos from synchronous rotation, we predict that Dimorphos's orbit will shrink by \dot a ~ 7 cm/yr, which can be measured by the Hera mission. We also speculate that the Yarkovsky force may have synchronized the Dinkinesh-Selam system after a possible merger of Selam's two lobes., Comment: 14 pages, 3 figures. Accepted for publication in the Astrophysical Journal Letters

Multi-wavelength scattered light imaging of debris disks may inform dust properties including typical size and mineral composition. Existing studies have investigated a small set of individual systems across a variety of imaging instruments and filters, calling for uniform comparison studies to systematically investigate dust properties. We obtain the surface brightness of dust particles in debris disks by post-processing coronagraphic imaging observations, and compare the multi-wavelength reflectance of dust. For a sample of resolved debris disks, we perform a systematic analysis on the reflectance properties of their birth rings. We reduced the visible and near-infrared images of 23 debris disk systems hosted by A through M stars using two coronagraphs onboard the Hubble Space Telescope: the STIS instrument observations centering at 0.58 $\mu$m, and the NICMOS instrument at 1.12 $\mu$m or 1.60 $\mu$m. For proper recovery of debris disks, we used classical reference differential imaging for STIS, and adopted non-negative matrix factorization with forward modeling for NICMOS. By dividing disk signals by stellar signals to take into account of intrinsic stellar color effects, we systematically obtained and compared the reflectance of debris birth rings at ~90 deg scattering angle. Debris birth rings typically exhibit a blue color at ~90 deg scattering angle. As the stellar luminosity increases, the color tends to be more neutral. A likely L-shaped color-albedo distribution indicates a clustering of scatterer properties. The observed color trend correlates with the expected blow-out size of dust particles. The color-albedo clustering likely suggests different populations of dust in these systems. More detailed radiative transfer models with realistic dust morphology will contribute to explaining the observed color and color-albedo distribution of debris systems., Comment: 17 pages, 8 figures, 3 tables. A&A accepted

The Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effect plays an important role in the rotational properties and evolution of asteroids. While the YORP effect induced by the macroscopic shape of the asteroid and by the presence of surface boulders has been well studied, no investigation has been performed yet regarding how craters with given properties influence this effect. We introduce and estimate the crater-induced YORP effect (CYORP), which arises from the concave structure of the crater, to investigate the magnitude of the resulting torques as a function of varying properties of the crater and the asteroid by a semi-analytical method. The CYORP torque can be comparable to the normal YORP torque when the size of the crater is about one-tenth of the size of the asteroid, or equivalently when the crater/roughness covers one-tenth of the asteroid's surface. Although the torque decreases with the crater size, the combined contribution of all small craters can become non-negligible due to their large number when the commonly used power-law crater size distribution is considered. The CYORP torque of small concave structures, usually considered as surface roughness, is essential to the accurate calculation of the complete YORP torque. Under the CYORP effect that is produced by collisions, asteroids go through a random walk in spin rate and obliquity, with a YORP reset timescale typically of 0.4 Myr. This has strong implications for the rotational evolution and orbital evolution of asteroids. Craters and roughness on asteroid surfaces can influence the YORP torques and therefore the rotational properties and evolution of asteroids. We suggest that the CYORP effect should be considered in the future investigation of the YORP effect on asteroids., Comment: 13 pages, 7 figures. Accepted by A&A

We scrutinize the light sail scenario of the first interstellar object (ISO) 1I/2017 U1 (`Oumuamua) by making comparisons between physical models and observational constraints. These analyses can be generalized for future surveys of `Oumuamua-like objects. The light sail goes through a drift in the interstellar space due to the magnetic field and gas atoms, which poses challenges to the navigation system. When the light sail enters the inner solar system, the sideways radiation pressure leads to a considerable non-radial displacement. The immensely high dimensional ratio and the tumbling motion could cause a light curve with an extremely large amplitude and could even make the light sail invisible from time to time. These observational features allow us to examine the light sail scenario of interstellar objects. Our results show that the drift of the freely rotating light sail in the interstellar medium is $\sim 100\,$au even if the travel distance is only 1 pc. The probability is < 1.5\% for the expected brightness modulation of the light sail to match with `Oumuamua's observed variation amplitude ($\sim$ 2.5 -- 3). In addition, the probability is 0.4% for the tumbling light-sail to be visible (brighter than V=27) in all 55 observations spread over two months after discovery. Radiation pressure could cause a larger displacement that is normal to the orbital plane for a lightsail than that for `Oumuamua. Also, the ratio of anti-solar to sideways acceleration of `Oumuamua deviates from that of the light sail by ~ 1.5{\sigma}. We suggest that `Oumuamua is unlikely a light sail. The dynamics of an intruding light sail, if exist, has distinct observational signatures, which can be quantitatively identified and analyzed with our methods in future surveys., Comment: 17 pages, 13 figures Comments are welcome

The first interstellar object 'Oumuamua is discovered in 2017. When 'Oumuamua travels in interstellar space, it keeps colliding with interstellar medium (ISM). Given a sufficiently long interaction time, its rotation state may change significantly because of the angular momentum transfer with interstellar medium. Using generated Gaussian random spheres with the dimension ratios 6:1:1 and 5:5:1, this paper explores the ISM torque curve and proposes that ISM collision may account for 'Oumuamua's tumbling with the simple constant-torque analytical method. The statistic results show that the asymptotic obliquities distribute mostly at 0 and 180 degree and most cases spin down at the asymptotic obliquity, indicating the ISM collision effect is similar to the YORP effect with zero heat conductivity assumed. Given a long time of deceleration of the spin rate, an initial major-axis rotation may evolve into tumbling motion under ISM torque. Using a constant-torque analytical model, the timescales of evolving into tumbling for the sample of 200 shapes are found to range from several Gyrs to ens of Gyrs, highly dependent on the chosen shape. The mean value is about 8.5 Gyrs for prolate shapes and 7.3 Gyrs for oblate shapes. Rotation of asteroids in the Oort cloud might be also dominated by the ISM collision effect since the YORP effect is quite weak at such a long distance from the Sun. Although this paper assumes an ideal mirror reflection and a constant relative velocity of 'Oumuamua, the results still show the importance of the ISM collision effect., Comment: 11 pages, 5 figures, accepted by ApJ. Comments are welcome! Email: u3006384@connect.hku.hk