Search

Principal regular satellites of gas giants are thought to be formed by the accumulation of solid materials in circumplanetary disks (CPDs). While there has been significant progress in the study of satellite formation in CPDs, details of the supply of satellite building blocks to CPDs remain unclear. We performed orbital integration of solid particles in the protoplanetary disk (PPD) approaching a planet, considering the gas drag force using the results of three-dimensional hydrodynamical simulations of a local region around the planet. We investigated planetary-mass dependence of the capture positions and capture rates of dust particles accreting onto the CPD. We also examined the degree of dust retention in accreting gas onto the CPD, which is important for determining the ratio of dust-to-gas inflow rates, a key parameter in satellite formation. We found that the degree of dust retention increases with increasing planetary mass for a given dust scale height in the PPD. In the case of a small planet ($M_{\rm p}=0.2M_{\rm Jup}$), most particles with insufficient initial altitudes in the PPD are isolated from the gas in the accreting region. On the other hand, in the case of a massive planet ($M_{\rm p}=1M_{\rm Jup}$), dust particles can be coupled to the vertically accreting gas, even when the dust scale height is about $10-30$\% of the gas scale height. The results of this study can be used for models of dust delivery and satellite formation in the CPDs of gas giants of various masses, including exoplanets., Comment: 27 pages, 10 figures, accepted to publication in Astrophysical Journal

Gas accretion onto the circumplanetary disks and the source region of accreting gas are important to reveal dust accretion that leads to satellite formation around giant planets. We performed local three-dimensional high-resolution hydrodynamic simulations of isothermal and inviscid gas flow around a planet to investigate planetary-mass dependence of gas accretion band width and gas accretion rate onto circumplanetary disks. We examined cases with various planetary masses corresponding to M_p=0.05-1M_{Jup} at 5.2 au, where M_{Jup} is the current Jovian mass. We found that the radial width of the gas accretion band is proportional to M_p^{1/6} for the low-mass regime with M_p < 0.2 M_{Jup} while it is proportional to M_p for the high-mass regime with M_p > 0.2M_{Jup}. We found that the ratio of the mass accretion rate onto the circumplanetary disk to that into the Hill sphere is about 0.4 regardless of planetary mass for the cases we examined. Combining our results with the gap model obtained from global hydrodynamic simulations, we derive semi-analytical formulae of mass accretion rate onto circumplanetary disks. We found that the mass dependence of our three-dimensional accretion rates is the same as the previously-obtained two-dimensional case, although the qualitative behavior of accretion flow onto the CPD is quite different between the two cases., Comment: accepted for publication in Astrophysical Journal, 24 pages, 13 figures

Furrows are a concentric system of tectonic troughs, and are the oldest recognizable surface feature on Ganymede. We analyzed the distribution of furrows utilizing Voyager and Galileo images and found that furrows over Ganymede's surface are part of a global concentric circular structure. If this multi-ring structure is impact origin, this is the largest impact structure identified so far in the solar system. Deviations of the shapes of the furrows from the concentricity are small everywhere, which implies that the relative location of the blocks of the dark terrains over the entire surface of Ganymede has not changed appreciably even during formation of the bright terrains. The estimate of the impactor size is difficult, but an 150km-radius impactor is consistent with the observed properties of furrows. The furrow-forming impact should have significant effects on the satellite's geological and internal evolution, which are expected to be confirmed by future explorations of Jupiter's icy moons, such as the JUICE (Jupiter Icy moon Explorer) or Europa Clipper mission., Comment: 35 pages 12 figures

Spacecraft and ground-based observations show that the main rings of Saturn lack particles larger than 10 m. Tidal or collisional destruction of satellites/comets have been proposed as the origin of the main rings; however, Saturn's tide alone cannot grind km-sized fragments into submeter-sized particles because of the high mechanical strength of water ice and rock. The question arises as to why such large particles are not left in the current ring. It is known that thermal stress induced by diurnal and seasonal temperature variations can cause weathering and fragmentation of boulders and contribute to dust and regolith production on the Moon and terrestrial planets, and then such thermal stress can break particles larger than a critical radius while cannot smaller than the critical radius. In this study, we examined the role of thermal stress acting on Saturn's ring particles. We found that thermal stress can grind porous ring particles larger than 10-20 m, which explains the lack of particles larger than 10 m in Saturn's ring. Also, fragmentation by thermal stress can be adoptable for the Epsilon rings of Uranus. Furthermore, thermal stress caused by diurnal or seasonal temperature variation acting on boulders on surfaces of icy satellites and asteroids may play an important role in the evolution of their sizes. Our calculations explain the lack of boulders on icy satellites, except in the geologically active provinces such as the tiger stripes of Enceladus, where boulders are supplied by recent geological activity. We predict that future observations can find numerous boulders around Europa's geologically active cracks., Comment: 35 pages, 8 figures

The asteroid (3200) Phaethon is known to be the parent body of the Geminids, although meteor showers are commonly associated with the activity of periodic comets. What is most peculiar to the asteroid is its comet-like activity in the ejection of micrometer-sized dust particles at every perihelion passage, while the activity of the asteroid has never been identified outside the near-perihelion zone at $0.14~\mathrm{au}$ from the Sun. From the theoretical point of view, we argue that the activity of the asteroid is well explained by the electrostatic lofting of micrometer-sized dust particles with the aid of mobile alkali ions at high temperatures. The mass-loss rates of micrometer-sized particles from the asteroid in our model is entirely consistent with the values inferred from visible observations of Phaethon's dust tail. For millimeter-sized particles, we predict three orders of magnitudes higher mass-loss rates, which could also account for the total mass of the Geminid meteoroid stream by the electrostatic lofting mechanism., Comment: 7 pages. 2 figures, to appear in Icarus Letters

We present an analysis of survey observations of the trailing L5 Jupiter Trojan swarm using the wide-field Hyper Suprime-Cam CCD camera on the 8.2 m Subaru Telescope. We detected 189 L5 Trojans from our survey that covered about 15 deg^2 of sky with a detection limit of m_r = 24.1 mag, and selected an unbiased sample consisting of 87 objects with absolute magnitude 14 < H_r < 17 corresponding to diameter 2 km < D < 10 km for analysis of size distribution. We fit their differential magnitude distribution to a single-slope power-law with an index \alpha = 0.37 +- 0.01, which corresponds to a cumulative size distribution with an index of b = 1.85 +- 0.05. Combining our results with data for known asteroids, we obtained the size distribution of L5 Jupiter Trojans over the entire size range for 9 < H_V < 17, and found that the size distributions of the L4 and L5 swarms agree well with each other for a wide range of sizes. This is consistent with the scenario that asteroids in the two swarms originated from the same primordial population. Based on the above results, the ratio of the total number of asteroids with D > 2 km in the two swarms is estimated to be N_L4/N_L5=1.40 +- 0.15, and the total number of L_5 Jupiter Trojans with D > 1 km is estimated to be 1.1 x 10^5 by extrapolating the obtained distribution.

We performed a wide-field survey observation of small asteroids using the Hyper Suprime-Cam installed on the 8.2 m Subaru Telescope. We detected more than 3,000 main-belt asteroids with a detection limit of 24.2 mag in the r-band, which were classified into two groups (bluish C-like and reddish S-like) by the g-r color of each asteroid and obtained size distributions of each group. We found that the shapes of size distributions of asteroids with the C-like and S-like colors agree with each other in the size range of 0.4-5 km in diameter. Assuming the asteroid population in this size range is under collision equilibrium, our results indicate that compositional difference hardly affects the size dependence of impact strength, at least for the size range between several hundred meters and several kilometers. This size range corresponds to the size range of ``spin-barrier'', an upper limit observed in the rotation rate distribution. Our results are consistent with the view that most asteroids in this size range have a rubble-pile structure., Comment: 23 pages, 11 figures, accepted for publication in The Astronomical Journal

Rotation periods of 53 small (diameters $2 < D < 40$ km) Jupiter Trojans (JTs) were derived using the high-cadence light curves obtained by the FOSSIL phase I survey, a Subaru/Hyper Suprime-Cam intensive program. These are the first reported periods measured for JTs with $D < 10$ km. We found a lower limit of the rotation period near 4 hr, instead of the previously published result of 5 hr (Ryan et al. 2017; Szabo et al. 2017, 2020) found for larger JTs. Assuming a rubble-pile structure for JTs, a bulk density of 0.9 gcm$^{-3}$ is required to withstand this spin rate limit, consistent with the value $0.8-1.0$ gcm$^{-3}$ (Marchis et al. 2006; Mueller et al. 2010; Buie et al. 2015; Berthier et al. 2020) derived from the binary JT system, (617) Patroclus-Menoetius system., Comment: Accepted by PSJ on July 9th, 2021. 15 pages, 7 figures, and 3 tables

Centaurs have orbits between Jupiter and Neptune and are thought to originate from the trans-Neptunian region. Observations of surface properties of Centaurs and comparison with those of trans-Neptunian objects (TNOs) would provide constraints on their origin and evolution. We analyzed imaging data of nine known Centaurs observed by the Hyper Suprime-Cam (HSC) installed on the Subaru Telescope with the g and i band filters. Using the data available in the public HSC data archive as well as those obtained by the HSC Subaru Strategic Program (HSC-SSP) by the end of June, 2017, we obtained the g-i colors of the nine Centaurs. We compared them with those of known TNOs in the HSC-SSP data obtained by Terai et al. (2018). We found that the color distribution of the nine Centaurs is similar to that of those TNOs with high orbital inclinations, but distinct from those TNOs with low orbital inclinations. We also examined correlations between the colors of these Centaurs and their orbital elements and absolute magnitude. The Centaurs' colors show a moderate positive correlation with semi-major axis, while no significant correlations between the color and other orbital elements or absolute magnitude were found for these Centaurs. On the other hand, recent studies on Centaurs with larger samples show interesting correlations between their color and absolute magnitude or orbital inclination. We discuss how our data fit in these previous studies, and also discuss implications of these results for their origin and evolution., Comment: 12 pages, 4 figures. Accepted for publication in Publications of the Astronomical Society of Japan

We present a visible multi-band photometry of trans-Neptunian objects (TNOs) observed by the Subaru Telescope in the framework of Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) from March in 2014 to September in 2016. We measured the five broad-band (g, r, i, z, and Y) colors over the wavelength range from 0.4 um to 1.0 um for 30 known TNOs using the HSC-SSP survey data covering ~500 deg2 of sky within +/-30 deg of ecliptic latitude. This dataset allows us to characterize the dynamical classes based on visible reflectance spectra as well as to examine the relationship between colors and the other parameters such as orbital elements. Our results show that the hot classical and scattered populations share similar color distributions, while the cold classical population has a reflective decrease toward shorter wavelength below the i band. Based on the obtained color properties, we found that the TNO sample examined in the present work can be separated into two groups by inclination (I), the low-I population consisting of cold classical objects and high-I population consisting of hot classical and scattered objects. The whole sample exhibits an anti-correlation between colors and inclination, but no significant correlation between colors and semi-major axis, perihelion distance, eccentricity, or absolute magnitude. The color-inclination correlation does not seem to be continuous over the entire inclination range. Rather, it is seen only in the high-I population. We found that the low- and high-I populations are distinguishable in the g-i vs. eccentricity plot, but four high-I objects show g-i colors similar to those of the low-I population. If we exclude these four objects, the high-I objects show a positive correlation between g-i and eccentricity and a negative correlation between g-i and inclination with high significance levels., Comment: 29 pages, 11 figure, 5 tables, submitted to HSC special issue in PASJ

Regular satellites of giant planets are formed by accretion of solid bodies in circumplanetary disks. Planetesimals that are moving on heliocentric orbits and are sufficiently large to be decoupled from the flow of the protoplanetary gas disk can be captured by gas drag from the circumplanetary disk. In the present work, we examine the distribution of captured planetesimals in circumplanetary disks using orbital integrations. We find that the number of captured planetesimals reaches an equilibrium state as a balance between continuous capture and orbital decay into the planet. The number of planetesimals captured into retrograde orbits is much smaller than those on prograde orbits, because the former ones experience strong headwind and spiral into the planet rapidly. We find that the surface number density of planetesimals at the current radial location of regular satellites can be significantly enhanced by gas drag capture, depending on the velocity dispersions of planetesimals and the width of the gap in the protoplanetary disk. Using a simple model, we also examine the ratio of the surface densities of dust and captured planetesimals in the circumplanetary disk, and find that solid material at the current location of regular satellites can be dominated by captured planetesimals when the velocity dispersion of planetesimals is rather small and a wide gap is not formed in the protoplanetary disk. In this case, captured planetesimals in such a region can grow by mutual collision before spiraling into the planet, and would contribute to the growth of regular satellites., Comment: 13 pages, 15 figures, accepted for publication in ApJ

The origin of rings around giant planets remains elusive. Saturn's rings are massive and made of 90-95% of water ice. In contrast, the much less massive rings of Uranus and Neptune are dark and likely to have higher rock fraction. Here we investigate, for the first time, the tidal disruption of a passing object, including the subsequent formation of planetary rings. First, we perform SPH simulations of the tidal destruction of big differentiated objects ($M_{\rm body}=10^{21-23}$) that experience close encounters with Saturn or Uranus. We find that about $0.1-10$% of the mass of the passing body is gravitationally captured around the planet. However, these fragments are initially big chunks and have highly eccentric orbits around the planet. Then, we perform N-body simulations including the planet's oblateness, starting with data obtained from the SPH simulations. Our N-body simulations show that the chunks are tidally destroyed during their next several orbits. Their individual orbits then start to precess incoherently around the planet's equator, which enhances their encounter velocities on longer-term evolution, resulting in more destructive impacts. These collisions would damp their eccentricities resulting in a progressive collapse of the debris cloud into a thin equatorial and low-eccentricity ring. These high energy impacts are expected to be catastrophic enough to produce small particles. Our numerical results also show that the mass of formed rings is large enough to explain current rings including inner regular satellites around Saturn and Uranus. In the case of Uranus, a body can go deeper inside the planet's Roche limit resulting in a more efficient capture of rocky material compared to Saturn's case in which mostly ice is captured. Thus, our results can naturally explain the compositional difference between the rings of Saturn, Uranus and Neptune., Comment: 33 pages, 21 figures; accepted in Icarus

Centaurs are minor planets orbiting between Jupiter and Neptune that have or had crossing orbits with one or more giant planets. Recent observations and reinterpretation of previous observations have revealed the existence of ring systems around 10199 Chariklo and 2060 Chiron. However, the origin of the ring systems around such a minor planet is still an open question. Here, we propose that the tidal disruption of a differentiated object that experiences a close encounter with a giant planet could naturally form diverse ring-satellite systems around the Centaurs. During the close encounter, the icy mantle of the passing object is preferentially ripped off by the planet's tidal force and the debris is distributed mostly within the Roche limit of the largest remnant body. Assuming the existence of $20-50$wt% silicate core below the icy mantle, a disk of particles is formed when the objects pass within $0.4-0.8$ of the planet's Roche limit with the relative velocity at infinity $3-6$km s$^{-1}$ and 8h initial spin period of the body. The resultant ring mass is $0.1-10$% of the central object's mass. Such particle disks are expected to spread radially, and materials spreading beyond the Roche limit would accrete into satellite(s). Our numerical results suggest that ring formation would be a natural outcome of such extreme close encounters and Centaurs can naturally have such ring systems because they cross the orbits of the giant planets., Comment: 8 pages, 3 figures, published in ApJL

When gas giant protoplanets grow sufficiently massive, circumplanetary disks would form. While solid bodies captured by the circumplanetary disks likely contribute to the growth of the planets and regular satellites around them, some of captured bodies would remain in planet-centered orbits after the dispersal of the disk. We examine capture and subsequent orbital evolution of planetesimals in waning circumplanetary gas disks using three-body orbital integration. We find that some of captured planetesimals can survive in the circumplanetary disk for a long period of time under such weak gas drag. Captured planetesimals have semi-major axes smaller than about one third of the planet's Hill radius. Distributions of their eccentricities and inclinations after disk dispersal depend on the strength of gas drag and the timescale of disk dispersal, and initially strong gas drag and quick disk dispersal facilitates capture and survival of planetesimals. However, in such a case, final orbital eccentricities and inclinations of captured bodies remain rather large. Although our results suggest that some of the present irregular satellites of gas giant planets with small semi-major axes would have been captured by gas drag, other mechanisms are required to fully explain their current orbital characteristics., Comment: 46 pages, 18 figures, accepted for publication in ApJ

Sufficiently massive growing giant planets have circumplanetary disks, and the capture of solid bodies by the disks would likely influence the growth of the planets and formation of satellite systems around them. In addition to dust particles that are supplied to the disk with inflowing gas, recent studies suggest the importance of capture of planetesimals whose motion is decoupled from the gas, but orbital evolution of captured bodies in circumplanetary gas disks has not been studied in detail. In the present work, using three-body orbital integration and analytic calculations, we examine orbital characteristics and subsequent dynamical evolution of planetesimals captured by gas drag from circumplanetary gas disks. We find that the semi-major axes of the planet-centered orbits of planetesimals at the time of permanent capture are smaller than about one third of the planet's Hill radius in most cases. Typically, captured bodies rapidly spiral into the planet, and the rate of the orbital decay is faster for the retrograde orbits due to the strong headwind from the circumplanetary gas. When a planetesimal captured into a retrograde orbit suffers from sufficiently strong gas drag before spiraling into the planet, its orbit turns to the prograde direction at a radial location that can be explained using the Stokes number. We also find that those captured into certain types of orbits can survive for a long period of time even under gas drag both in the prograde and retrograde cases, which may be important for the origin of irregular satellites of giant planets., Comment: 37 pages, 11 figures, accepted for publication in AJ

Circumplanetary particle disks would be created in the late stage of planetary formation either by impacts of planetary bodies or disruption of satellites or passing bodies, and satellites can be formed by accretion of disk particles spreading across the Roche limit. Previous N-body simulation of lunar accretion focused on the formation of single-satellite systems from disks with large disk-to-planet mass ratios, while recent models of the formation of multiple-satellite systems from disks with smaller mass ratios do not take account of gravitational interaction between formed satellites. In the present work, we investigate satellite accretion from particle disks with various masses, using N-body simulation. In the case of accretion from somewhat less massive disks than the case of lunar accretion, formed satellites are not massive enough to clear out the disk, but can become massive enough to gravitationally shepherd the disk outer edge and start outward migration due to gravitational interaction with the disk. When the radial location of the 2:1 mean motion resonance of the satellite reaches outside the Roche limit, the second satellite can be formed near the disk outer edge, and then the two satellites continue outward migration while being locked in the resonance. Co-orbital satellites are found to be occasionally formed on the orbit of the first satellite. Our simulations also show that stochastic nature involved in gravitational interaction and collision between aggregates in the tidal environment can lead to diversity in the final mass and orbital architecture, which would be expected in satellite systems of exoplanets., Comment: 50 pages, 17 figures, 1 table, accepted for publication in ApJ

The thermal inertia values of Saturn's main rings (the A, B, and C rings and the Cassini division) are derived by applying our thermal model to azimuthally scanned spectra taken by the Cassini Composite Infrared Spectrometer (CIRS). Model fits show the thermal inertia of ring particles to be 16, 13, 20, and 11 Jm$^{-2}$K$^{-1}$s$^{-1/2}$ for the A, B, and C rings, and the Cassini division, respectively. However, there are systematic deviations between modeled and observed temperatures in Saturn's shadow depending on solar phase angle, and these deviations indicate that the apparent thermal inertia increases with solar phase angle. This dependence is likely to be explained if large slowly spinning particles have lower thermal inertia values than those for small fast spinning particles because the thermal emission of slow rotators is relatively stronger than that of fast rotators at low phase and vise versa. Additional parameter fits, which assume that slow and fast rotators have different thermal inertia values, show the derived thermal inertia values of slow (fast) rotators to be 8 (77), 8 (27), 9 (34), 5 (55) Jm$^{-2}$K$^{-1}$s$^{-1/2}$ for the A, B, and C rings, and the Cassini division, respectively. The values for fast rotators are still much smaller than those for solid ice with no porosity. Thus, fast rotators are likely to have surface regolith layers, but these may not be as fluffy as those for slow rotators, probably because the capability of holding regolith particles is limited for fast rotators due to the strong centrifugal force on surfaces of fast rotators. Other additional parameter fits, in which radii of fast rotators are varied, indicate that particles less than $\sim$ 1 cm should not occupy more than a half of the cross section for the A, B, and C rings., Comment: 47 pages, 7 tables, 13 figures

We investigate gas accretion flow onto a circumplanetary disk from a protoplanetary disk in detail by using high-resolution three-dimensional nested-grid hydrodynamic simulations, in order to provide a basis of formation processes of satellites around giant planets. Based on detailed analyses of gas accretion flow, we find that most of gas accretion onto circumplanetary disks occurs nearly vertically toward the disk surface from high altitude, which generates a shock surface at several scale heights of the circumplanetary disk. The gas that has passed through the shock surface moves inward because its specific angular momentum is smaller than that of the local Keplerian rotation, while gas near the midplane in the protoplanetary disk cannot accrete to the circumplanetary disk. Gas near the midplane within the planet's Hill sphere spirals outward and escapes from the Hill sphere through the two Lagrangian points L$_1$ and L$_2$. We also analyze fluxes of accreting mass and angular momentum in detail and find that the distributions of the fluxes onto the disk surface are well described by power-law functions and that a large fraction of gas accretion occurs at the outer region of the disk, i.e., at about 0.1 times the Hill radius. The nature of power-law functions indicates that, other than the outer edge, there is no specific radius where gas accretion is concentrated. These source functions of mass and angular momentum in the circumplanetary disk would provide us with useful constraints on the structure and evolution of the circumplanetary disk, which is important for satellite formation., Comment: 22pages, 17 figures, accepted for publication in ApJ

When protoplanets growing by accretion of planetesimals have atmospheres, small planetesimals approaching the protoplanets lose their energy by gas drag from the atmospheres, which leads them to be captured within the Hill sphere of the protoplanets. As a result, growth rates of the protoplanets are enhanced. In order to study the effect of an atmosphere on planetary growth rates, we performed numerical integration of orbits of planetesimals for a wide range of orbital elements and obtained the effective accretion rates of planetesimals onto planets that have atmospheres. Numerical results are obtained as a function of planetesimals' eccentricity, inclination, planet's radius, and non-dimensional gas-drag parameters which can be expressed by several physical quantities such as the radius of planetesimals and the mass of the protoplanet. Assuming that the radial distribution of the gas density near the surface can be approximated by a power-law, we performed analytic calculation for the loss of planetesimals' kinetic energy due to gas drag, and confirmed agreement with numerical results. We confirmed that the above approximation of the power-law density distribution is reasonable for accretion rate of protoplanets with one to ten Earth-masses, unless the size of planetesimals is too small. We also calculated the accretion rates of planetesimals averaged over a Rayleigh distribution of eccentricities and inclinations, and derived a semi-analytical formula of accretion rates, which reproduces the numerical results very well. Using the obtained expression of the accretion rate, we examined the growth of protoplanets in nebular gas. We found that the effect of atmospheric gas drag can enhance the growth rate significantly, depending on the size of planetesimals., Comment: 41 pages, 14 figures, accepted for publication in Icarus