Comets are thought to have information about the formation process of our solar system. Recently, detailed information about comet 67P/Churyumov-Gerasimenko has been found by a spacecraft mission Rosetta. It is remarkable that its tensile strength was estimated. In this paper, we measure and formulate the tensile strength of porous dust aggregates using numerical simulations, motivated by porous dust aggregation model of planetesimal formation. We perform three-dimensional numerical simulations using a monomer interaction model with periodic boundary condition. We stretch out a dust aggregate with a various initial volume filling factor between $10^{-2}$ and 0.5. We find that the tensile stress takes the maximum value at the time when the volume filling factor decreases to about a half of the initial value. The maximum stress is defined to be the tensile strength. We take an average of the results with 10 different initial shapes to smooth out the effects of initial shapes of aggregates. Finally, we numerically obtain the relation between the tensile strength and the initial volume filling factor of dust aggregates. We also use a simple semi-analytical model and successfully reproduce the numerical results, which enables us to apply to a wide parameter range and different materials. The obtained relation is consistent with previous experiments and numerical simulations about silicate dust aggregates. We estimate that the monomer radius of comet 67P has to be about 3.3--220 $\mathrm{\mu m}$ to reproduce its tensile strength using our model.

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.

Planetesimal formation is one of the most important unsolved problems in planet formation theory. In particular, rocky planetesimal formation is difficult because silicate dust grains are easily broken when they collide. Recently, it has been proposed that they can grow as porous aggregates when their monomer radius is smaller than $\sim$ 10 nm, which can also avoid the radial drift toward the central star. However, the stability of a layer composed of such porous silicate dust aggregates has not been investigated. Therefore, we investigate the gravitational instability of this dust layer. To evaluate the disk stability, we calculate Toomre's stability parameter $Q$, for which we need to evaluate the equilibrium random velocity of dust aggregates. We calculate the equilibrium random velocity considering gravitational scattering and collisions between dust aggregates, drag by mean flow of gas, stirring by gas turbulence, and gravitational scattering by gas density fluctuation due to turbulence. We derive the condition of the gravitational instability using the disk mass, dust-to-gas ratio, turbulent strength, orbital radius, and dust monomer radius. We find that, for the minimum mass solar nebula model at 1 au, the dust layer becomes gravitationally unstable when the turbulent strength $\alpha\lesssim10^{-5}$. If the dust-to-gas ratio is increased twice, the gravitational instability occurs for $\alpha\lesssim10^{-4}$. We also find that the dust layer is more unstable in disks with larger mass, higher dust-to-gas ratio, and weaker turbulent strength, at larger orbital radius, and with a larger monomer radius.

We validate a $R_p=2.32\pm 0.24R_\oplus$ planet on a close-in orbit ($P=2.260455\pm 0.000041$ days) around K2-28 (EPIC 206318379), a metal-rich M4-type dwarf in the Campaign 3 field of the K2 mission. Our follow-up observations included multi-band transit observations from the optical to the near infrared, low-resolution spectroscopy, and high-resolution adaptive-optics (AO) imaging. We perform a global fit to all the observed transits using a Gaussian process-based method and show that the transit depths in all passbands adopted for the ground-based transit follow-ups ($r'_2, z_\mathrm{s,2}, J, H, K_\mathrm{s}$) are within $\sim 2\sigma$ of the K2 value. Based on a model of the background stellar population and the absence of nearby sources in our AO imaging, we estimate the probability that a background eclipsing binary could cause a false positive to be $< 2\times 10^{-5}$. We also show that K2-28 cannot have a physically associated companion of stellar type later than M4, based on the measurement of almost identical transit depths in multiple passbands. There is a low probability for a M4 dwarf companion ($\approx 0.072_{-0.04}^{+0.02}$), but even if this were the case, the size of K2-28b falls within the planetary regime. K2-28b has the same radius (within $1\sigma$) and experiences a similar irradiation from its host star as the well-studied GJ~1214b. Given the relative brightness of K2-28 in the near infrared ($m_\mathrm{Kep}=14.85$ mag and $m_H=11.03$ mag) and relatively deep transit ($0.6-0.7\%$), a comparison between the atmospheric properties of these two planets with future observations would be especially interesting.