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

We investigate the patterns observed in ejecta curtain induced by hypervelocity impact (2-6 km/s) with a variety of the size and shape of target particles. We characterize the patterns by an angle, defined as the ratio of the characteristic length of the pattern obtained by Fourier transformation to the distance from the impact point. This angle is found to be almost the same as that obtained by the reanalysis of the patterns in the previous study at lower impact velocities (Kadono et al., 2015, Icarus 250, 215-221), which are consistent with lunar crater-ray systems. Assuming that the pattern is formed by mutual collision of particles with fluctuation velocity in excavation flow, we evaluate an angle at which the pattern growth stops and show that this angle is the same in the order of magnitude as the ratio of the fluctuation velocity and the radial velocity. This relation is confirmed in the results of experiments and numerical simulations. Finally, we discuss the dependence of the patterns on impact conditions. The experiments show no dependence of the angle on impact velocity. This indicates that the ratio between the fluctuation and radial velocity components in excavation flow does not depend on impact velocity. Moreover, the independences on particle size and particle shape suggest that the angle characterizing the structure of the patterns does not depend on cohesive force. Since cohesive forces should be related with elastic properties of particles, the structure does not depend on elastic properties, though inelastic collisions are important for the persistence and contrast of the patterns., Comment: 22 pages and 7 figures, accepted for publication in Icarus

We performed impact experiments to observe patterns in an ejecta curtain with targets consisting of small sand particles and large inclusions comparable to or smaller than the size of the projectiles. The spatial intensity distributions in the ejecta at early stages of crater formation depend on the size of the inclusions. Our numerical simulations of radially spreading particles with different sizes support this result. Based on the results, we proposed a procedure for evaluating the subsurface structures of celestial bodies from the images of ejecta curtains obtained from space-impact experiments., Comment: 16 pages and 3 figures, accepted for publication in The Astrophysical Journal Letters

In most of the previous studies related to collisional disruption of planetesimals in the gravity regime, Smoothed Particle Hydrodynamics (SPH) simulations have been used. On the other hand, impact simulations using grid-based hydrodynamic code have not been sufficiently performed. In the present study, we execute impact simulations in the gravity regime using the shock-physics code iSALE, which is a grid-based Eulerian hydrocode. We examine the dependence of the critical specific impact energy Q_RD* on impact conditions for a wide range of specific impact energy (Q_R) from disruptive collisions to erosive collisions, and compare our results with previous studies. We find collision outcomes of the iSALE simulation agree well with those of the SPH simulation. Detailed analysis mainly gives three results. (1) The value of Q_RD* depends on numerical resolution, and is close to convergence with increasing numerical resolution. The difference in converged value of Q_RD* between the iSALE code and the SPH code is within 30%. (2) Ejected mass normalized by total mass (M_ej/M_tot) generally depends on various impact conditions. However, when Q_R is normalized by Q_RD* that is calculated for each impact simulation, M_ej/M_tot can be scaled by Q_R/Q_RD*, and is independent of numerical resolution, impact velocity and target size. (3) This similarity law for Q_R/Q_RD* is confirmed for a wide range of specific impact energy. We also derive a semi-analytic formula for Q_RD* based on the similarity law and the crater scaling law. We find that the semi-analytic formula for the case with a non-porous object is consistent with numerical results., Comment: Accepted for publication in Icarus

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

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