The Rosetta mission to comet 67P/Churyumov-Gerasimenko has provided new data to better understand what comets are made of. The weak tensile strength of the cometary surface materials suggests that the comet is a hierarchical dust aggregate formed through gravitational collapse of a bound clump of small dust aggregates so-called ``pebbles'' in the gaseous solar nebula. Since pebbles are the building blocks of comets, which are the survivors of planetesimals in the solar nebula, estimating the size of pebbles using a combination of thermal observations and numerical calculations is of great importance to understand the planet formation in the outer solar system. In this study, we calculated the thermal inertias and thermal skin depths of the hierarchical aggregates of pebbles, for both diurnal and orbital variations of the temperature. We found that the thermal inertias of the comet 67P/Churyumov-Gerasimenko are consistent with the hierarchical aggregate of cm- to dm-sized pebbles. Our findings indicate that the icy planetesimals may have formed via accretion of cm- to dm-sized pebbles in the solar nebula.

Understanding the physical properties of dust aggregates is of great importance in planetary science. In this study, we revisited the sticking property of submillimetre-sized aggregates. We revealed that the "effective surface energy" model used in previous studies underestimates the critical pulling force needed to separate two sticking aggregates. We also derived a new and simple model of the critical pulling force based on the canonical theory of two contacting spheres. Our findings indicate that we do not need to consider the "effective surface energy" of dust aggregates when discussing the physical properties of loose agglomerates of submillimetre-sized aggregates.

We herein report a theoretical study of the geometrical structure of porous dust aggregates formed via ballistic cluster-cluster aggregation (BCCA). We calculated the gyration radius $R_{\rm gyr}$ and the graph-based geodesic radius $R_{\rm geo}$ as a function of the number of constituent particles $N$. We found that $R_{\rm gyr} / r_{0} \sim N^{0.531 \pm 0.011}$ and $R_{\rm geo} / r_{0} \sim N^{0.710 \pm 0.013}$, where $r_{0}$ is the radius of constituent particles. Furthermore, we defined two constants that characterize the geometrical structure of fractal aggregates: $D_{\rm f}$ and $\alpha$. The definition of $D_{\rm f}$ and $\alpha$ are $N \sim {( R_{\rm gyr} / r_{0} )}^{D_{\rm f}}$ and ${R_{\rm geo}} / {r_{0}} \sim {\left( {R_{\rm gyr}} / {r_{0}} \right)}^{\alpha}$, respectively. Our study revealed that $D_{\rm f} \simeq 1.88$ and $\alpha \simeq 1.34$ for the clusters of the BCCA. In addition, we also studied the filling factor dependence of thermal conductivity of statically compressed fractal aggregates. From this study, we reveal that the thermal conductivity of statically compressed aggregates $k$ is given by $k \sim 2 k_{\rm mat} {( r_{\rm c} / r_{0} )} \phi^{(1 + \alpha) / (3 - D_{\rm f})}$, where $k_{\rm mat}$ is the material thermal conductivity, $r_{\rm c}$ is the contact radius of constituent particles, and $\phi$ is the filling factor of dust aggregates.

Jupiter has four large regular satellites called the Galilean satellites: Io, Europa, Ganymede, and Callisto. The inner three of the Galilean satellites orbit in a 4:2:1 mean motion resonance; therefore their orbital configuration may originate from the stopping of the migration of Io near the bump in the surface density distribution and following resonant trapping of Europa and Ganymede. The formation mechanism of the bump near the orbit of the innermost satellite, Io, is not yet understood, however. Here, we show that photophoresis in the circumjovian disk could be the cause of the bump, using analytic calculations of steady-state accretion disks. We propose that photophoresis in the circumjovian disk could stop the inward migration of dust particles near the orbit of Io. The resulting dust depleted inner region would have a higher ionization fraction, and thus admit increased magnetorotational-instability-driven accretion stress than the outer region. The increase of the accretion stress at the photophoretic dust barrier would form a bump in the surface density distribution, halting the migration of Io.

Recent studies have revealed that all large (over 1000 km in diameter) trans-Neptunian objects (TNOs) form satellite systems. Although the largest Plutonian satellite, Charon, is thought to be an intact fragment of an impactor directly formed via a giant impact, whether giant impacts can explain the variations in secondary-to-primary mass ratios and spin/orbital periods among all large TNOs remains to be determined. Here we systematically perform hydrodynamic simulations to investigate satellite formation via giant impacts. We find that the simulated secondary-to-primary mass ratio varies over a wide range, which overlaps with observed mass ratios. We also reveal that the satellite systems' current distribution of spin/orbital periods and small eccentricity can be explained only when their spins and orbits tidally evolve: initially as fluid-like bodies, but finally as rigid bodies. These results suggest that all satellites of large TNOs were formed via giant impacts in the early stage of solar system formation, before the outward migration of Neptune, and that they were fully or partially molten during the giant impact era.

Shock-wave heating within the solar nebula is one of the leading candidates for the source of chondrule-forming events. Here, we examine the possibility of compound chondrule formation via optically thin shock waves. Several features of compound chondrules indicate that compound chondrules are formed via the collisions of supercooled precursors. We evaluate whether compound chondrules can be formed via the collision of supercooled chondrule precursors in the framework of the shock-wave heating model by using semi-analytical methods and discuss whether most of the crystallized chondrules can avoid destruction upon collision in the post-shock region. We find that chondrule precursors immediately turn into supercooled droplets when the shock waves are optically thin and they can maintain supercooling until the condensation of evaporated fine dust grains. Owing to the large viscosity of supercooled melts, supercooled chondrule precursors can survive high-speed collisions on the order of $1\ {\rm km}\ {\rm s}^{-1}$ when the temperature is below $\sim 1400\ {\rm K}$. From the perspective of the survivability of crystallized chondrules, shock waves with a spatial scale of $\sim 10^{4}\ {\rm km}$ may be potent candidates for the chondrule formation mechanism. Based on our results from one-dimensional calculations, a fraction of compound chondrules can be reproduced when the chondrule-to-gas mass ratio in the pre-shock region is $\sim 2 \times 10^{-3}$, which is approximately half of the solar metallicity.

Understanding the heat transfer mechanism within dust aggregates is of great importance for many subjects in planetary science. We calculated the coordination number and the thermal conductivity through the solid network of compressed dust aggregates. We found a simple relationship between the coordination number and the filling factor and revealed that the thermal conductivity through the solid network of aggregates is represented by a power-law function of the filling factor and the coordination number.

$\mathit{Context.}$ The thermal conductivity of highly porous dust aggregates is a key parameter for many subjects in planetary science; however, it is not yet fully understood. $\mathit{Aims.}$ In this study, we investigate the thermal conductivity of fluffy dust aggregates with filling factors of less than $10^{-1}$. $\mathit{Methods.}$ We determine the temperature structure and heat flux of the porous dust aggregates calculated by $N$-body simulations of static compression in the periodic boundary condition. $\mathit{Results.}$ We derive an empirical formula for the thermal conductivity through the solid network $k_{\rm sol}$ as a function of the filling factor of dust aggregates $\phi$. The results reveal that $k_{\rm sol}$ is approximately proportional to ${\phi}^{2}$, and the thermal conductivity through the solid network is significantly lower than previously assumed. In light of these findings, we must reconsider the thermal histories of small planetary bodies.

Chondritic meteorites primarily contain millimeter-sized spherical objects, the so-called chondrules; however, the co-accretion process of chondrules and matrix grains is not yet understood. In this study, we investigate the ejection process of chondrules via collisions of fluffy aggregates composed of chondrules and matrices. We reveal that fluffy aggregates cannot grow into planetesimals without losing chondrules if we assume that the chondrite parent bodies are formed via direct aggregation of similar-sized aggregates. Therefore, an examination of other growth pathways is necessary to explain the formation of rocky planetesimals in our solar system.

Several pieces of evidence suggest that silicate grains in primitive meteorites are not interstellar grains but condensates formed in the early solar system. Moreover, the size distribution of matrix grains in chondrites implies that these condensates might be formed as nanometer-sized grains. Therefore, we propose a novel scenario for rocky planetesimal formation in which nanometer-sized silicate grains are produced by evaporation and recondensation events in early solar nebula, and rocky planetesimals are formed via aggregation of these nanograins. We reveal that silicate nanograins can grow into rocky planetesimals via direct aggregation without catastrophic fragmentation and serious radial drift, and our results provide a suitable condition for protoplanet formation in our solar system.

We present a novel model showing that compound chondrules are formed by collisions of supercooled droplets. This model reproduces two prominent observed features of compound chondrules: the nonporphyritic texture and the size ratio between two components.