Source author record

Hiroshi Kobayashi

Hiroshi Kobayashi appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

31works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

31 published item(s)

preprint2022arXiv

Formation of dust clumps with sub-Jupiter mass and cold shadowed region in gravitationally unstable disk around Class 0/I protostar in L1527 IRS

We have investigated the protostellar disk around a Class 0/I protostar, L1527 IRS, using multi-wavelength observations of the dust continuum emission at $λ=0.87$, 2.1, 3.3, and 6.8 mm obtained by the Atacama Large Millimeter/submillimeter Array (ALMA) and the Jansky Very Large Array (VLA). Our observations achieved a spatial resolution of $3-13$ au and revealed an edge-on disk structure with a size of $\sim80-100$ au. The emission at 0.87 and 2.1 mm is found to be optically thick within a projected disk radius of $ r_{\rm proj}\lesssim50$ au. The emission at 3.3 and 6.8 mm shows that the power-law index of the dust opacity ($β$) is $β\sim1.7$ around $ r_{\rm proj}\sim 50$ au, suggesting that grain growth has not yet begun. The dust temperature ($T_{\rm dust}$) shows a steep decrease with $T_{\rm dust}\propto r_{\rm proj}^{-2}$ outside of the VLA clumps previously identified at $r_{\rm proj}\sim20$ au. Furthermore, the disk is gravitationally unstable at $r_{\rm proj}\sim20$ au, as indicated by a Toomre {\it Q} parameter value of $Q\lesssim1.0$. These results suggest that the VLA clumps are formed via gravitational instability, which creates a shadow on the outside of the substructure, resulting in the sudden drop in temperature. The derived dust masses for the VLA clumps are $\gtrsim0.1$ $M_{\rm J}$. Thus, we suggest that Class 0/I disks can be massive enough to be gravitationally unstable, which might be the origin of gas-giant planets in a 20 au radius. Furthermore, the protostellar disks can be cold due to shadowing.

preprint2022arXiv

No evidence of the significant grain growth but tentative discovery of disk substructure in a disk around the Class I Protostar L1489 IRS

For revealing the first step of the plant formation, it is important to understand how and when dust grains become larger in a disk around a protostar. To investigate the grain growth, we analyze dust continuum emission toward a disk around the Class I protostar, L1489 IRS at 0.9 and 1.3 mm wavelengths obtained by the Atacama Large Millimeter/submillimeter Array. The dust continuum emission extends to a disk radius ($r$) of $r\sim300$ au, and the spectral index ($α$) is derived to be $α\sim3.6$ at a radius of $r\sim100-300$ au, as similar to the interstellar dust. Therefore, the grain growth does not occur significantly in the outer disk ($r\sim100-300$ au). Furthermore, we tentatively identify a ring-like substructure at $r\sim90$ au even though the spatial resolution and sensitivity are not enough to determine this structure. If this is the real ring structure, the ring position and small dust in the disk outer part are consistent with the idea of the growth front. These results suggest that the L1489 protostellar disk may be the beginning of the planet formation.

preprint2020arXiv

Ring formation by coagulation of dust aggregates in early phase of disk evolution around a protostar

Ring structures are observed by (sub-)millimeter dust continuum emission in various circumstellar disks from early stages of Class 0 and I to late stage of Class II young stellar objects (YSOs). In this paper, we study one of the possible scenarios of such ring formation in early stage, which is coagulation of dust aggregates. The dust grains grow in an inside-out manner because the growth timescale is roughly proportional to the orbital period. The boundary of the dust evolution can be regarded as the growth front, where the growth time is comparable to the disk age. With radiative transfer calculations based on the dust coagulation model, we find that the growth front can be observed as a ring structure because dust surface density is sharply changed at this position. Furthermore, we confirm that the observed ring positions in the YSOs with an age of $\lesssim1$ Myr are consistent with the growth front. The growth front could be important to create the ring structure in particular for early stage of the disk evolution such as Class 0 and I sources.

preprint2019arXiv

High-Resolution Simulations of Catastrophic Disruptions: Resultant Shape Distributions

The members of asteroid families have various shapes. We investigate the origin of their shapes by high-resolution impact simulations for catastrophic disruptions using a Smoothed Particle Hydrodynamics code. Collisional remnants produced through our simulations of the catastrophic disruptions mainly have spherical or bilobed shapes. However, no flat remnants with the ratio of minor to major axis lengths c/a < 0.5 are formed. The results of the simulations provide various shapes of asteroids and explain most of the shapes in asteroid families that are supposed to be produced through catastrophic disruptions. However, the present simulations do not explain significantly flat asteroids. We suggest that these flat asteroids may be interlopers or formed through low-velocity collisions between member asteroids.

preprint2019arXiv

Importance of Giant Impact Ejecta for Orbits of Planets Formed during the Giant Impact Era

Terrestrial planets are believed to be formed via giant impacts of Mars-sized protoplanets. Planets formed via giant impacts have highly eccentric orbits. A swarm of planetesimals around the planets may lead to eccentricity damping for the planets via the equipartition of random energies (dynamical friction). However, dynamical friction increases eccentricities of planetesimals, resulting in high velocity collisions between planetesimals. The collisional cascade grinds planetesimals to dust until dust grains are blown out due to radiation pressure. Therefore, the total mass of planetesimals decreases due to collisional fragmentation, which weakens dynamical friction. We investigate the orbital evolution of protoplanets in a planetesimal disk, taking into account collisional fragmentation of planetesimals. For 100 km-sized or smaller planetesimals, dynamical friction is insignificant for eccentricity damping of planets because of collisional fragmentation. On the other hand, giant impacts eject collisional fragments. Although the total mass of giant impact ejecta is 0.1-0.3 Earth masses, the largest impact ejecta are ~ 1,000 km in size. We also investigate the orbital evolution of single planets with initial eccentricities 0.1 in a swarm of such giant impact ejecta. Although the total mass of giant impact ejecta decreases by a factor of 3 in 30 Myrs, eccentricities of planets are damped down to the Earth level (~0.01) due to interaction with giant impact ejecta. Therefore, giant impact ejecta play an important role for determination of terrestrial planet orbits.

preprint2016arXiv

Cohesion of Amorphous Silica Spheres: Toward a Better Understanding of the Coagulation Growth of Silicate Dust Aggregates

Adhesion forces between submicrometer-sized silicate grains play a crucial role in the formation of silicate dust agglomerates, rocky planetesimals, and terrestrial planets. The surface energy of silicate dust particles is the key to their adhesion and rolling forces in a theoretical model based on the contact mechanics. Here we revisit the cohesion of amorphous silica spheres by compiling available data on the surface energy for hydrophilic amorphous silica in various circumstances. It turned out that the surface energy for hydrophilic amorphous silica in a vacuum is a factor of 10 higher than previously assumed. Therefore, the previous theoretical models underestimated the critical velocity for the sticking of amorphous silica spheres, as well as the rolling friction forces between them. With the most plausible value of the surface energy for amorphous silica spheres, theoretical models based on the contact mechanics are in harmony with laboratory experiments. Consequently, we conclude that silicate grains with a radius of $0.1~μ$m could grow to planetesimals via coagulation in a protoplanetary disk. We argue that the coagulation growth of silicate grains in a molecular cloud is advanced either by organic mantles rather than icy mantles or, if there are no mantles, by nanometer-sized grain radius.

preprint2016arXiv

RZ Leonis Minoris Bridging between ER Ursae Majoris-Type Dwarf Nova and Novalike System

We observed RZ LMi, which is renowned for the extremely (~19d) short supercycle and is a member of a small, unusual class of cataclysmic variables called ER UMa-type dwarf novae, in 2013 and 2016. In 2016, the supercycles of this object substantially lengthened in comparison to the previous measurements to 35, 32, 60d for three consecutive superoutbursts. We consider that the object virtually experienced a transition to the novalike state (permanent superhumper). This observed behavior extremely well reproduced the prediction of the thermal-tidal instability model. We detected a precursor in the 2016 superoutburst and detected growing (stage A) superhumps with a mean period of 0.0602(1)d in 2016 and in 2013. Combined with the period of superhumps immediately after the superoutburst, the mass ratio is not as small as in WZ Sge-type dwarf novae, having orbital periods similar to RZ LMi. By using least absolute shrinkage and selection operator (Lasso) two-dimensional power spectra, we detected possible negative superhumps with a period of 0.05710(1)d. We estimated the orbital period of 0.05792d, which suggests a mass ratio of 0.105(5). This relatively large mass ratio is even above ordinary SU UMa-type dwarf novae, and it is also possible that the exceptionally high mass-transfer rate in RZ LMi may be a result of a stripped core evolved secondary which are evolving toward an AM CVn-type object.

preprint2016arXiv

Sintering-induced Dust Ring Formation in Protoplanetary Disks: Application to the HL Tau Disk

The latest observation of HL Tau by ALMA revealed spectacular concentric dust rings in its circumstellar disk. We attempt to explain the multiple ring structure as a consequence of aggregate sintering. Sintering is known to reduce the sticking efficiency of dust aggregates and occurs at temperatures slightly below the sublimation point of their constituent material. We here present a dust growth model incorporating sintering and use it to simulate global dust evolution due to sintering, coagulation, fragmentation, and radial inward drift in a modeled HL Tau disk. We show that aggregates consisting of multiple species of volatile ices experience sintering, collisionally disrupt, and pile up at multiple locations slightly outside the snow lines of the volatiles. At wavelengths of 0.87--1.3 mm, these sintering zones appear as bright, optically thick rings with a spectral slope of $\approx 2$, whereas the non-sintering zones as darker, optically thinner rings of a spectral slope of $\approx$ 2.3--2.5. The observational features of the sintering and non-sintering zones are consistent with those of the major bright and dark rings found in the HL Tau disk, respectively. Radial pileup and vertical settling occur simultaneously if disk turbulence is weak and if monomers constituting the aggregates are $\sim 1~{\rm μm}$ in radius. For the radial gas temperature profile of $T = 310(r/1~{\rm AU})^{-0.57}~{\rm K}$, our model perfectly reproduces the brightness temperatures of the optically thick bright rings, and reproduces their orbital distances to an accuracy of $\lesssim$ 30%.

preprint2015arXiv

Formation of terrestrial planets in disks evolving via disk winds and implications for the origin of the solar system's terrestrial planets

Recent three-dimensional magnetohydrodynamical simulations have identified a disk wind by which gas materials are lost from the surface of a protoplanetary disk, which can significantly alter the evolution of the inner disk and the formation of terrestrial planets. A simultaneous description of the realistic evolution of the gaseous and solid components in a disk may provide a clue for solving the problem of the mass concentration of the terrestrial planets in the solar system. We simulate the formation of terrestrial planets from planetary embryos in a disk that evolves via magnetorotational instability and a disk wind. The aim is to examine the effects of a disk wind on the orbital evolution and final configuration of planetary systems. We perform N-body simulations of sixty 0.1 Earth-mass embryos in an evolving disk. The evolution of the gas surface density of the disk is tracked by solving a one-dimensional diffusion equation with a sink term that accounts for the disk wind. We find that even in the case of a weak disk wind, the radial slope of the gas surface density of the inner disk becomes shallower, which slows or halts the type I migration of embryos. If the effect of the disk wind is strong, the disk profile is significantly altered (e.g., positive surface density gradient, inside-out evacuation), leading to outward migration of embryos inside ~ 1 AU. Disk winds play an essential role in terrestrial planet formation inside a few AU by changing the disk profile. In addition, embryos can undergo convergent migration to ~ 1 AU in certainly probable conditions. In such a case, the characteristic features of the solar system's terrestrial planets (e.g., mass concentration around 1 AU, late giant impact) may be reproduced.

preprint2015arXiv

From Planetesimals to Planets in Turbulent Protoplanetary Disks I. Onset of Runaway Growth

When planetesimals grow via collisions in a turbulent disk, stirring through density fluctuation caused by turbulence effectively increases the relative velocities between planetesimals, which suppresses the onset of runaway growth. We investigate the onset of runaway growth in a turbulent disk through simulations that calculate the mass and velocity evolution of planetesimals. When planetesimals are small, the average relative velocity between planetesimals, $v_{\rm r}$, is much greater than their surface escape velocity, $v_{\rm esc}$, so that runaway growth does not occur. As planetesimals become large via collisional growth, $v_{\rm r}$ approaches $v_{\rm esc}$. When $v_{\rm r} \approx 1.5 v_{\rm esc}$, runaway growth of the planetesimals occurs. During the oligarchic growth subsequent to runaway growth, a small number of planetary embryos produced via runaway growth become massive through collisions with planetesimals with radii of that at the onset of runaway growth, $r_{\rm p,run}$. We analytically derive $r_{\rm p,run}$ as a function of the turbulent strength. Growing $\sim 10\,M_\oplus$ embryos that are suitable to become the cores of Jupiter and Saturn requires $r_{\rm p,run} \sim 100$\,km, which is similar to the proposed fossil feature in the size distribution of main belt asteroids. In contrast, the formation of Mars as quickly as suggested from Hf-W isotope studies requires small planetesimals at the onset of runaway growth. Thus, the conditions required to form Mars, Jupiter, and Saturn and the size distribution of the main-belt asteroids indicate that the turbulence increased in amplitude relative to the sound speed with increasing distance from the young Sun.

preprint2015arXiv

Orbital evolution of planetesimals in gaseous disks

Planets are formed from collisional growth of small bodies in a protoplanetary disk. Bodies much larger than approximately $1$\,m are mainly controlled by the gravity of the host star and experience weak gas drag; their orbits are mainly expressed by orbital elements: semimajor axes $a$, eccentricities $e$, and inclinations $i$, which are modulated by gas drag. In a previous study, $\dot a$, $\dot e$, and $\dot i$ were analytically derived for $e \ll 1$ and $i \ll H/a$, where $H$ is the scale height of the disk. Their formulae are valid in the early stage of planet formation. However, once massive planets are formed, $e$ and $i$ increase greatly. Indeed, some small bodies in the solar system have very large $e$ and $i$. Therefore, in this paper, I analytically derive formulae for $\dot a$, $\dot e$, and $\dot i$ for $1-e^2 \ll 1$ and $i \ll H/a$ and for $i \gg H/a$. The formulae combined from these limited equations will represent the results of orbital integration unless $e \geq 1$ or $i > π- H/a$. Since the derived formulae are applicable for bodies not only in a protoplanetary disk but also in a circumplanetary disk, I discuss the possibility of the capture of satellites in a circumplanetary disk using the formulae.

preprint2015arXiv

Planetary and meteoritic Mg/Si and d30Si variations inherited from solar nebula chemistry

The bulk chemical compositions of planets are uncertain, even for major elements such as Mg and Si. This is due to the fact that the samples available for study all originate from relatively shallow depths. Comparison of the stable isotope compositions of planets and meteorites can help overcome this limitation. Specifically, the non-chondritic Si isotope composition of the Earth's mantle was interpreted to reflect the presence of Si in the core, which can also explain its low density relative to pure Fe-Ni alloy. However, we have found that angrite meteorites display a heavy Si isotope composition similar to the lunar and terrestrial mantles. Because core formation in the angrite parent-body (APB) occurred under oxidizing conditions at relatively low pressure and temperature, significant incorporation of Si in the core is ruled out as an explanation for this heavy Si isotope signature. Instead, we show that equilibrium isotopic fractionation between gaseous SiO and solid forsterite at 1370 K in the solar nebula could have produced the observed Si isotope variations. Nebular fractionation of forsterite should be accompanied by correlated variations between the Si isotopic composition and Mg/Si ratio following a slope of 1, which is observed in meteorites. Consideration of this nebular process leads to a revised Si concentration in the Earth's core of 3.6 (+6.0/-3.6) wt% and provides estimates of Mg/Si ratios of bulk planetary bodies.

preprint2015arXiv

Significant Gas-to-Dust Ratio Asymmetry and Variation in the Disk of HD 142527 and the Indication of Gas Depletion

We investigate the dust and gas distribution in the disk around HD 142527 based on ALMA observations of dust continuum, 13CO(3-2), and C18O(3-2) emission. The disk shows strong azimuthal asymmetry in the dust continuum emission, while gas emission is more symmetric. In this paper, we investigate how gas and dust are distributed in the dust-bright northern part of the disk and in the dust-faint southern part. We construct two axisymmetric disk models. One reproduces the radial profiles of the continuum and the velocity moments 0 and 1 of CO lines in the north and the other reproduces those in the south. We have found that the dust is concentrated in a narrow ring having ~50AU width (in FWHM; w_d=30AU in our parameter definition) located at ~170-200AU from the central star. The dust particles are strongly concentrated in the north. We have found that the dust surface density contrast between the north and south amounts to ~70. Compared to the dust, the gas distribution is more extended in the radial direction. We find that the gas component extends at least from ~100AU to ~250AU from the central star, and there should also be tenuous gas remaining inside and outside of these radii. The azimuthal asymmetry of gas distribution is much smaller than dust. The gas surface density differs only by a factor of ~3-10 between the north and south. Hence, gas-to-dust ratio strongly depends on the location of the disk: ~30 at the location of the peak of dust distribution in the south and ~3 at the location of the peak of dust distribution in the north. Despite large uncertainties, the overall gas-to-dust ratio is inferred to be ~10-30, indicating that the gas depletion may have already been under way.

preprint2014arXiv

Debris disc formation induced by planetary growth

Several hundred stars older than 10 million years have been observed to have infrared excesses. These observations are explained by dust grains formed by the collisional fragmentation of hidden planetesimals. Such dusty planetesimal discs are known as debris discs. In a dynamically cold planetesimal disc, collisional coagulation of planetesimals produces planetary embryos which then stir the surrounding leftover planetesimals. Thus, the collisional fragmentation of planetesimals that results from planet formation forms a debris disc. We aim to determine the properties of the underlying planetesimals in debris discs by numerically modelling the coagulation and fragmentation of planetesimal populations. The brightness and temporal evolution of debris discs depend on the radial distribution of planetesimal discs, the location of their inner and outer edges, their total mass, and the size of planetesimals in the disc. We find that a radially narrow planetesimal disc is most likely to result in a debris disc that can explain the trend of observed infrared excesses of debris discs around G-type stars, for which planet formation occurs only before 100 million years. Early debris disc formation is induced by planet formation, while the later evolution is explained by the collisional decay of leftover planetesimals around planets that have already formed. Planetesimal discs with underlying planetesimals of radii $\sim 100\,$km at $\approx 30$ AU most readily explain the Spitzer Space Telescope 24 and 70$ μ$m fluxes from debris discs around G-type stars.

preprint2014arXiv

N-body Simulations of Terrestrial Planet Formation under the Influence of a Hot Jupiter

We investigate the formation of multiple-planet systems in the presence of a hot Jupiter using extended N-body simulations that are performed simultaneously with semi-analytic calculations. Our primary aims are to describe the planet formation process starting from planetesimals using high-resolution simulations, and to examine the dependences of the architecture of planetary systems on input parameters (e.g., disk mass, disk viscosity). We observe that protoplanets that arise from oligarchic growth and undergo type I migration stop migrating when they join a chain of resonant planets outside the orbit of a hot Jupiter. The formation of a resonant chain is almost independent of our model parameters, and is thus a robust process. At the end of our simulations, several terrestrial planets remain at around 0.1 AU. The formed planets are not equal-mass; the largest planet constitutes more than 50 percent of the total mass in the close-in region, which is also less dependent on parameters. In the previous work of this paper (Ogihara et al. 2013), we have found a new physical mechanism of induced migration of the hot Jupiter, which is called a crowding-out. If the hot Jupiter opens up a wide gap in the disk (e.g., owing to low disk viscosity), crowding-out becomes less efficient and the hot Jupiter remains. We also discuss angular momentum transfer between the planets and disk.

preprint2014arXiv

Study of Negative and Positive Superhumps in ER Ursae Majoris

We carried out the photometric observations of the SU UMa-type dwarf nova ER UMa during 2011 and 2012, which showed the existence of persistent negative superhumps even during the superoutburst. We performed two-dimensional period analysis of its light curves by using a method called "least absolute shrinkage and selection operator" (Lasso) and "phase dispersion minimization" (PDM) analysis, and we found that the period of negative superhumps systematically changed between a superoutburst and the next superoutburst. The trend of the period change can beinterpreted as reflecting the change of the disk radius. This change of the disk radius is in good agreement with the predicted change of the disk radius by the thermal-tidal instability (TTI) model. The normal outbursts within a supercycle showed a general trend that the rising rate to maximum becomes slower as the next superoutburst approaches. The change can be interpreted as the consequence of the increased gas-stream flow onto the inner region of the disk as the result of the tilted disk. Some of the superoutbursts were found to be triggered by a precursor normal outburst when the positive superhumps appeared to develop. The positive and negative superhumps co-existed during the superoutburst. The positive superhumps were prominent only during four or five days after the supermaximum, while the signal of the negative superhumps became strong after the middle phase of the superoutburst plateau. A simple combination of the positive and negative superhumps was found to be insufficient in reproducing the complex profile variation. We were able to detect the developing phase of positive superhumps (stage A superhumps) for the first time in ER UMa-type dwarf novae. Using the period of stage A superhumps, we obtained a mass ratio of 0.100(15), which indicates that ER UMa is on the ordinary evolutional track of CVs.

preprint2013arXiv

Condition for Capture into First-order Mean Motion Resonances and Application to Constraints on Origin of Resonant Systems

We investigate the condition for capture into first-order mean motion resonances using numerical simulations with a wide range of various parameters. In particular, we focus on deriving the critical migration timescale for capture into the 2:1 resonance; additional numerical experiments for closely spaced resonances (e.g., 3:2) are also performed. We find that the critical migration timescale is determined by the planet-to-stellar mass ratio, and its dependence exhibits power-law behavior with index -4/3. This dependence is also supported by simple analytic arguments. We also find that the critical migration timescale for systems with equal-mass bodies is shorter than that in the restricted problem; for instance, for the 2:1 resonance between two equal-mass bodies, the critical timescale decreases by a factor of 10. In addition, using the obtained formula, the origin of observed systems that include first-order commensurabilities is constrained. Assuming that pairs of planets originally form well separated from each other and then undergo convergent migration and are captured in resonances, it is possible that a number of exoplanets experienced rapid orbital migration. For systems in closely spaced resonances, the differential migration timescale between the resonant pair can be constrained well; it is further suggested that several exoplanets underwent migration that can equal or even exceed the type I migration rate predicted by the linear theory. This implies that some of them may have formed in situ. Future observations and the use of our model will allow us to statistically determine the typical migration speed in a protoplanetary disk.

preprint2013arXiv

Crowding Out of Giants by Dwarfs: An Origin for the Lack of Companion Planets in Hot Jupiter Systems

We investigate formation of close-in terrestrial planets from planetary embryos under the influence of a hot Jupiter (HJ) using gravitational N-body simulations that include gravitational interactions between the gas disk and the terrestrial planet (e.g., type I migration). Our simulations show that several terrestrial planets efficiently form outside the orbit of the HJ, making a chain of planets, and all of them gravitationally interact directly or indirectly with the HJ through resonance, which leads to inward migration of the HJ. We call this mechanism of induced migration of the HJ as "crowding out." The HJ is eventually lost by collision with the central star, and only several terrestrial planets remain. We also find that the efficiency of the crowding-out effect depends on model parameters; for example, the heavier the disk is, the more efficient the crowding out is. When planet formation occurs in a massive disk, the HJ can be lost to the central star and is never observed. On the other hand, for a less massive disk, the HJ and terrestrial planets can coexist; however, the companion planets can be below the detection limit of current observations. In both cases, systems with the HJ and terrestrial planets have little chance for detection. Therefore, our model naturally explains the lack of companion planets in HJ systems regardless of the disk mass. In effect, our model provide a theoretical prediction for future observations; additional planets can be discovered just outside the HJ, and their masses should generally be small.

preprint2013arXiv

Evolution of dust grain size distribution by shattering in the interstellar medium: robustness and uncertainty

Shattering of dust grains in the interstellar medium is a viable mechanism of small grain production in galaxies. We examine the robustness or uncertainty in the theoretical predictions of shattering. We identify $P_1$ (the critical pressure above which the deformation destroys the original lattice structures) as the most important quantity in determining the timescale of small grain production, and confirm that the same $P_1/t$ ($t$ is the duration of shattering) gives the same grain size distribution [$n(a)$, where $a$ is the grain radius] after shattering within a factor of 3. The uncertainty in the fraction of shocked material that is eventually ejected as fragments causes uncertainties in $n(a)$ by a factor of 1.3 and 1.6 for silicate and carbonaceous dust, respectively. The size distribution of shattered fragments have minor effects as long as the power index of the fragment size distribution is less than ~ 3.5, since the slope of grain size distribution $n(a)$ continuously change by shattering and becomes consistent with $n(a)\propto a^{-3.5}$. The grain velocities as a function of grain radius can have an imprint in the grain size distribution especially for carbonaceous dust. We also show that the formulation of shattering can be simplified without losing sufficient precision.

preprint2013arXiv

Local Enhancement of Surface Density in the Protoplanetary Ring Surrounding HD 142527

We report ALMA observations of dust continuum, 13CO J=3--2, and C18O J=3--2 line emission toward a gapped protoplanetary disk around HD 142527. The outer horseshoe-shaped disk shows the strong azimuthal asymmetry in dust continuum with the contrast of about 30 at 336 GHz between the northern peak and the southwestern minimum. In addition, the maximum brightness temperature of 24 K at its northern area is exceptionally high at 160 AU from a star. To evaluate the surface density in this region, the grain temperature needs to be constrained and was estimated from the optically thick 13CO J=3--2 emission. The lower limit of the peak surface density was then calculated to be 28 g cm-2 by assuming a canonical gas-to-dust mass ratio of 100. This finding implies that the region is locally too massive to withstand self-gravity since Toomre's Q <~1--2, and thus, it may collapse into a gaseous protoplanet. Another possibility is that the gas mass is low enough to be gravitationally stable and only dust grains are accumulated. In this case, lower gas-to-dust ratio by at least 1 order of magnitude is required, implying possible formation of a rocky planetary core.

preprint2013arXiv

Small Planetesimals in a Massive Disk Formed Mars

Mars is likely to be a planetary embryo formed through collisions with planetesimals, which can explain its small mass and rapid formation timescale obtained from 182Hf-182$W chronometry. In the classical theory of planet formation, the final embryo mass is determined only by the solid surface density. However, embryos can stir surrounding planetesimals, leading to fragmentation through erosive (cratering) collisions. We find that radial drift of small fragments can drastically reduce the solid surface density. On the other hand, embryo growth is accelerated by fragment accretion. Since collisional fragmentation efficiency depends on the initial size of planetesimals, the final embryo mass and its growth time are determined by the initial planetesimal size and disk surface density. We have investigated the effect of these two parameters on the mass of Mars and the predicted radiogenic excess of 182W in the martian mantle. Two scenarios can explain the rapid formation of small Mars: (i) it formed by accretion of small planetesimals in a massive disk or (ii) it formed from large planetesimals but its growth was arrested by the inward then outward migration of Jupiter. Taking into account all constraints, we conclude that Mars is likely to have formed in a massive disk of about ~ 0.1 solar mass from planetesimals smaller than ~ 10 km in radius. Such small planetesimal size cannot explain core accretion of Jupiter, suggesting that there may have been a heliocentric gradient in planetesimal size in the solar nebula.

preprint2013arXiv

Survey of Period Variations of Superhumps in SU UMa-Type Dwarf Novae. IV: The Fourth Year (2011-2012)

Continuing the project described by Kato et al. (2009, arXiv:0905.1757), we studied 86 SU UMa-type dwarf novae. We confirmed the general trends such as the relation between period derivatives and orbital periods. There are some systems showing positive period derivatives despite the long orbital periods. We observed the 2011 outburst of the WZ Sge-type dwarf nova BW Scl, and recorded an O-C diagram similar to those of previously known WZ Sge-type dwarf novae. The WZ Sge-type dwarf nova OT J184228.1+483742 showed an unusual pattern of double outbursts composed of an outburst with early superhumps and one with ordinary superhumps. We propose an interpretation that a very small growth rate of the 3:1 resonance due to an extremely low mass-ratio led to a quenching of the superoutburst before the ordinary superhumps appeared. We studied ER UMa-type dwarf novae and found that V1159 Ori showed positive superhumps similar to ER UMa in the 1990s. The recently recognized ER UMa-type object BK Lyn dominantly showed negative superhumps, and its behavior was very similar to the present-day state of ER UMa. The pattern of period variations in AM CVn-type objects were very similar to short-period hydrogen-rich SU UMa-type dwarf novae, making them helium analogue of hydrogen-rich SU UMa-type dwarf novae. SBS 1108+574, a peculiar hydrogen-rich dwarf nova below the period minimum, showed a very similar pattern of period variations to those of short-period SU UMa-type dwarf novae. The mass-ratio derived from the detected orbital period suggests that this secondary is a somewhat evolved star whose hydrogen envelope was mostly stripped during the mass-exchange. CC Scl, MASTER OT J072948.66+593824.4 and OT J173516.9+154708 showed only low-amplitude superhumps with complex profiles. These superhumps are likely a combination of closely separated two periods.

preprint2013arXiv

Survey of Period Variations of Superhumps in SU UMa-Type Dwarf Novae. V: The Fifth Year (2012-2013)

Continuing the project described by Kato et al. (2009a, arXiv:0905.1757), we collected times of superhump maxima for SU UMa-type dwarf novae mainly observed during the 2012-2013 season. We found three objects (V444 Peg, CSS J203937 and MASTER J212624) having strongly positive period derivatives despite the long orbital period (Porb). By using the period of growing stage (stage A) superhumps, we obtained mass ratios for six objects. We characterized nine new WZ Sge-type dwarf novae. We made a pilot survey of the decline rate of slowly fading part of SU UMa-type and WZ Sge-type outbursts. The decline time scale was found to generally follow the expected Porb^(1/4) dependence and WZ Sge-type outbursts also generally follow this trend. There are some objects which show slower decline rates, and we consider these objects good candidates for period bouncers. We also studied unusual behavior in some objects, including BK Lyn which made a transition from an ER UMa-type state to the novalike (standstill) state in 2013 and unusually frequent occurrence of superoutbursts in NY Ser and CR Boo. We applied least absolute shrinkage and selection operator (Lasso) power spectral analysis, which has been proven to be very effective in analyzing the Kepler data, to ground-based photometry of BK Lyn and detected the dramatic disappearance of the signal of negative superhumps in 2013. We suggested that the mass-transfer rates did not vary strongly between the ER UMa-type state and novalike state in BK Lyn, and this transition was less likely caused by a systematic variation of the mass-transfer rate.

preprint2012arXiv

An improved model of the Edgeworth-Kuiper debris disk

(Abridged) We access the expected EKB dust disk properties by modeling. We treat the debiased population of the known transneptunian objects (TNOs) as parent bodies and generate the dust with our collisional code. The resulting dust distributions are modified to take into account the influence of gravitational scattering and resonance trapping by planets on migrating dust grains as well as the effect of sublimation. A difficulty is that the amount and distribution of dust are largely determined by sub-kilometer-sized bodies. These are directly unobservable, and their properties cannot be accessed by collisional modeling, because objects larger than 10...60m in the present-day EKB are not in a collisional equilibrium. To place additional constraints, we use in-situ measurements of the New Horizons spacecraft within 20AU. We show that the TNO population has to have a break in the size distribution at s<70km. However, even this still leaves us with several models that all correctly reproduce a nearly constant dust impact rates in the region of giant planet orbits and do not violate the constraints from the non-detection of the EKB dust thermal emission by the COBE spacecraft. The modeled EKB dust disks, which conform to the observational constraints, can either be transport-dominated or intermediate between the transport-dominated and collision-dominated regime. The in-plane optical depth of such disks is tau(r>10AU)~10^-6 and their fractional luminosity is f_d~10^-7. Planets and sublimation are found to have little effect on dust impact fluxes and dust thermal emission. The spectral energy distribution of an EKB analog, as would be seen from 10pc distance, peaks at wavelengths of 40...50\mum at F~0.5mJy, which is less than 1% of the photospheric flux at those wavelengths. Therefore, exact EKB analogs cannot be detected with present-day instruments such as Herschel/PACS.

preprint2012arXiv

Rapid Coagulation of Porous Dust Aggregates Outside the Snow Line: A Pathway to Successful Icy Planetesimal Formation

Rapid orbital drift of macroscopic dust particles is one of the major obstacles against planetesimal formation in protoplanetary disks. We reexamine this problem by considering porosity evolution of dust aggregates. We apply a porosity model based on recent N-body simulations of aggregate collisions, which allows us to study the porosity change upon collision for a wide range of impact energies. As a first step, we neglect collisional fragmentation and instead focus on dust evolution outside the snow line, where the fragmentation has been suggested to be less significant than inside the snow line because of a high sticking efficiency of icy particles. We show that dust particles can evolve into highly porous aggregates (with internal densities of much less than 0.1 g/cm^3) even if collisional compression is taken into account. We also show that the high porosity triggers significant acceleration in collisional growth. This acceleration is a natural consequence of particles' aerodynamical property at low Knudsen numbers, i.e., at particle radii larger than the mean free path of the gas molecules. Thanks to this rapid growth, the highly porous aggregates are found to overcome the radial drift barrier at orbital radii less than 10 AU (assuming the minimum-mass solar nebula model). This suggests that, if collisional fragmentation is truly insignificant, formation of icy planetesimals is possible via direct collisional growth of submicron-sized icy particles.

preprint2012arXiv

Rapid Formation of Saturn after Jupiter Completion

We have investigated Saturn's core formation at a radial pressure maximum in a protoplanetary disk, which is created by gap opening by Jupiter. A core formed via planetesimal accretion induces the fragmentation of surrounding planetesimals, which generally inhibits further growth of the core by removal of the resulting fragments due to radial drift caused by gas drag. However, the emergence of the pressure maximum halts the drift of the fragments, while their orbital eccentricities and inclinations are efficiently damped by gas drag. As a result, the core of Saturn rapidly grows via accretion of the fragments near the pressure maximum. We have found that in the minimum-mass solar nebula, kilometer sized planetesimals can produce a core exceeding 10 Earth masses within two million years. Since Jupiter may not have undergone significant type II inward migration, it is likely that Jupiter's formation was completed when the local disk mass has already decayed to a value comparable to or less than Jovian mass. The expected rapid growth of Saturn's core on a timescale comparable to or shorter than observationally inferred disk lifetime enables Saturn to acquire the current amount of envelope gas before the disk gas is completely depleted. The high heat energy release rate onto the core surface due to the rapid accretion of the fragments delays onset of runaway gas accretion until the core mass becomes somewhat larger than that of Jupiter, which is consistent with the estimate based on interior modeling. Therefore, the rapid formation of Saturn induced by gap opening of Jupiter can account for the formation of multiple gas giants (Jupiter and Saturn) without significant inward migration and larger core mass of Saturn than that of Jupiter.

preprint2012arXiv

Survey of Period Variations of Superhumps in SU UMa-Type Dwarf Novae. III: The Third Year (2010--2011)

Continuing the project described by Kato et al. (2009, PASJ 61, S395, arXiv:0905.1757), we collected times of superhump maxima for 51 SU UMa-type dwarf novae mainly observed during the 2010-2011 season. Although most of the new data for systems with short superhump periods basically confirmed the findings by Kato et al. (2009) and Kato et al. (2010, PASJ 62, 1525, arXiv:1009.5444), the long-period system GX Cas showed an exceptionally large positive period derivative. An analysis of public Kepler data of V344 Lyr and V1504 Cyg yielded less striking stage transitions. In V344 Lyr, there was prominent secondary component growing during the late stage of superoutbursts, and the component persisted at least for two more cycles of successive normal outbursts. We also investigated the superoutbursts of two conspicuous eclipsing objects: HT Cas and the WZ Sge-type object SDSS J080434.20+510349.2. Strong beat phenomena were detected in both objects, and late-stage superhumps in the latter object had an almost constant luminosity during the repeated rebrightenings. The WZ Sge-type object SDSS J133941.11+484727.5 showed a phase reversal around the rapid fading from the superoutburst. The object showed a prominent beat phenomenon even after the end of the superoutburst. A pilot study of superhump amplitudes indicated that the amplitudes of superhumps are strongly correlated with orbital periods, and the dependence on the inclination is weak in systems with inclinations smaller than 80 deg.

preprint2011arXiv

Planetary Core Formation with Collisional Fragmentation and Atmosphere to Form Gas Giant Planets

Massive planetary cores ($\sim 10$ Earth masses) trigger rapid gas accretion to form gas giant planets \rev{such as} Jupiter and Saturn. We investigate the core growth and the possibilities for cores to reach such a critical core mass. At the late stage, planetary cores grow through collisions with small planetesimals. Collisional fragmentation of planetesimals, which is induced by gravitational interaction with planetary cores, reduces the amount of planetesimals surrounding them, and thus the final core masses. Starting from small planetesimals that the fragmentation rapidly removes, less massive cores are formed. However, planetary cores acquire atmospheres that enlarge their collisional cross section before rapid gas accretion. Once planetary cores exceed about Mars mass, atmospheres significantly accelerate the growth of cores. We show that, taking into account the effects of fragmentation and atmosphere, initially large planetesimals enable formation of sufficiently massive cores. On the other hand, because the growth of cores is slow for large planetesimals, a massive disk is necessary for cores to grow enough within a disk lifetime. If the disk with 100\,km-sized initial planetesimals is 10 times as massive as the minimum mass solar nebula, planetary cores can exceed 10 Earth masses in the Jovian planet region ($>5\,$AU).

preprint2011arXiv

Understanding how planets become massive. I. Description and validation of a new toy model

The formation of giant planets requires accumulation of ~10 Earth mass in solids; but how do protoplanets acquire their mass? There are many, often competing processes that regulate the accretion rate of protoplanets. To assess their effects we present a new, publicly-available toy model. The rationale behind the toy model is that it encompasses as many physically-relevant processes as possible, but at the same time does not compromise its simplicity, speed, and physical insight. The toy model follows a modular structure, where key features -- e.g. planetesimal fragmentation, radial orbital decay, nebula turbulence -- can be switched on or off. Our model assumes three discrete components (fragments, planetesimals, and embryos) and is zero dimensional in space. We have tested the outcomes of the toy model against literature results and generally find satisfactory agreement. We include, for the first time, model features that capture the three-way interactions among small particles, gas, and protoplanets. Collisions among planetesimals will result in fragmentation, transferring a substantial amount of the solid mass to small particles, which couple strongly to the gas. Our results indicate that the efficiency of the accretion process then becomes very sensitive to the gas properties -- especially to the turbulent state and the magnitude of the disk headwind (the decrease of the orbital velocity of the gas with respect to Keplerian) -- as well as to the characteristic fragment size.

preprint2010arXiv

Planetary Growth with Collisional Fragmentation and Gas Drag

As planetary embryos grow, gravitational stirring of planetesimals by embryos strongly enhances random velocities of planetesimals and makes collisions between planetesimals destructive. The resulting fragments are ground down by successive collisions. Eventually the smallest fragments are removed by the inward drift due to gas drag. Therefore, the collisional disruption depletes the planetesimal disk and inhibits embryo growth. We provide analytical formulae for the final masses of planetary embryos, taking into account planetesimal depletion due to collisional disruption. Furthermore, we perform the statistical simulations for embryo growth (which excellently reproduce results of direct $N$-body simulations if disruption is neglected). These analytical formulae are consistent with the outcome of our statistical simulations. Our results indicate that the final embryo mass at several AU in the minimum-mass solar nebula can reach about $\sim 0.1$ Earth mass within $10^7$ years. This brings another difficulty in formation of gas giant planets, which requires cores with $\sim 10$ Earth masses for gas accretion. However, if the nebular disk is 10 times more massive than the minimum-mass solar nebula and the initial planetesimal size is larger than 100 km, as suggested by some models of planetesimal formation, the final embryo mass reaches about 10 Earth masses at 3-4 AU. The enhancement of embryos' collisional cross sections by their atmosphere could further increase their final mass to form gas giant planets at 5-10 AU in the solar system.

preprint2010arXiv

Survey of Period Variations of Superhumps in SU UMa-Type Dwarf Novae. II: The Second Year (2009-2010)

As an extension of the project in Kato et al. (2009, arXiv:0905.1757), we collected times of superhump maxima for 61 SU UMa-type dwarf novae mainly observed during the 2009-2010 season. The newly obtained data confirmed the basic findings reported in Kato et al. (2009): the presence of stages A-C, as well as the predominance of positive period derivatives during stage B in systems with superhump periods shorter than 0.07 d. There was a systematic difference in period derivatives for systems with superhump periods longer than 0.075 d between this study and Kato et al. (2009). We suggest that this difference is possibly caused by the relative lack of frequently outbursting SU UMa-type dwarf novae in this period regime in the present study. We recorded a strong beat phenomenon during the 2009 superoutburst of IY UMa. The close correlation between the beat period and superhump period suggests that the changing angular velocity of the apsidal motion of the elliptical disk is responsible for the variation of superhump periods. We also described three new WZ Sge-type objects with established early superhumps and one with likely early superhumps. We also suggest that two systems, VX For and EL UMa, are WZ Sge-type dwarf novae with multiple rebrightenings. The O-C variation in OT J213806.6+261957 suggests that the frequent absence of rebrightenings in very short-Porb objects can be a result of sustained superoutburst plateau at the epoch when usual SU UMa-type dwarf novae return to quiescence preceding a rebrightening. We also present a formulation for a variety of Bayesian extension to traditional period analyses.