Source author record

Taku Takeuchi

Taku Takeuchi 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

17works
4topics
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

17 published item(s)

preprint2016arXiv

Mass constraint for a planet in a protoplanetary disk from the gap width

A giant planet creates a gap in a protoplanetary disk, which might explain the observed gaps in protoplanetary disks. The width and depth of the gaps depend on the planet mass and disk properties. We have performed two--dimensional hydrodynamic simulations for various planet masses, disk aspect ratios and viscosities, to obtain an empirical formula for the gap width. The gap width is proportional to the square root of the planet mass, -3/4 power of the disk aspect ratio and -1/4 power of the viscosity. This empirical formula enables us to estimate the mass of a planet embedded in the disk from the width of an observed gap. We have applied the empirical formula for the gap width to the disk around HL~Tau, assuming that each gap observed by ALMA observations is produced by planets, and discussed the planet masses within the gaps. The estimate of planet masses from the gap widths is less affected by the observational resolution and dust filtration than that from the gap depth.

preprint2016arXiv

Parametric Study of the Rossby Wave Instability in a Two-Dimensional Barotropic Disk

Protoplanetary disks with non-axisymmetric structures have been observed. The Rossby wave instability (RWI) is considered as one of the origins of the non-axisymmetric structures. We perform linear stability analyses of the RWI in barotropic flow using four representative types of the background flow on a wide parameter space. We find that the co-rotation radius is located at the background vortensity minimum with large concavity if the system is marginally stable to the RWI, and this allows us to check the stability against the RWI easily. We newly derive the necessary and sufficient condition for the onset of the RWI in semi-analytic form. We discuss the applicability of the new condition in realistic systems and the physical nature of the RWI.

preprint2015arXiv

Detailed structure of the outer disk around HD 169142 with polarized light in H-band

Coronagraphic imagery of the circumstellar disk around HD 169142 in H-band polarized intensity (PI) with Subaru/HiCIAO is presented. The emission scattered by dust particles at the disk surface in 0.2" <= r <= 1.2", or 29 <= r <= 174 AU, is successfully detected. The azimuthally-averaged radial profile of the PI shows a double power-law distribution, in which the PIs in r=29-52 AU and r=81.2-145 AU respectively show r^{-3}-dependence. These two power-law regions are connected smoothly with a transition zone (TZ), exhibiting an apparent gap in r=40-70 AU. The PI in the inner power-law region shows a deep minimum whose location seems to coincide with the point source at λ= 7 mm. This can be regarded as another sign of a protoplanet in TZ. The observed radial profile of the PI is reproduced by a minimally flaring disk with an irregular surface density distribution or with an irregular temperature distribution or with the combination of both. The depletion factor of surface density in the inner power-law region (r< 50 AU) is derived to be <= 0.16 from a simple model calculation. The obtained PI image also shows small scale asymmetries in the outer power-law region. Possible origins for these asymmetries include corrugation of the scattering surface in the outer region, and shadowing effect by a puffed up structure in the inner power-law region.

preprint2015arXiv

Formation of a disc gap induced by a planet: Effect of the deviation from Keplerian disc rotation

The gap formation induced by a giant planet is important in the evolution of the planet and the protoplanetary disc. We examine the gap formation by a planet with a new formulation of one-dimensional viscous discs which takes into account the deviation from Keplerian disc rotation due to the steep gradient of the surface density. This formulation enables us to naturally include the Rayleigh stable condition for the disc rotation. It is found that the derivation from Keplerian disc rotation promotes the radial angular momentum transfer and makes the gap shallower than in the Keplerian case. For deep gaps, this shallowing effect becomes significant due to the Rayleigh condition. In our model, we also take into account the propagation of the density waves excited by the planet, which widens the range of the angular momentum deposition to the disc. The effect of the wave propagation makes the gap wider and shallower than the case with instantaneous wave damping. With these shallowing effects, our one-dimensional gap model is consistent with the recent hydrodynamic simulations.

preprint2015arXiv

Mass Estimates of a Giant Planet in a Protoplanetary Disk from the Gap Structures

A giant planet embedded in a protoplanetary disk forms a gap. An analytic relationship among the gap depth, planet mass $M_{p}$, disk aspect ratio $h_p$, and viscosity $α$ has been found recently, and the gap depth can be written in terms of a single parameter $K= (M_{p}/M_{\ast})^2 h_p^{-5} α^{-1}$. We discuss how observed gap features can be used to constrain the disk and/or planet parameters based on the analytic formula for the gap depth. The constraint on the disk aspect ratio is critical in determining the planet mass so the combination of the observations of the temperature and the image can provide a constraint on the planet mass. We apply the formula for the gap depth to observations of HL~Tau and HD~169142. In the case of HL~Tau, we propose that a planet with $\gtrsim 0.3$ is responsible for the observed gap at $30$~AU from the central star based on the estimate that the gap depth is $\lesssim 1/3$. In the case of HD~169142, the planet mass that causes the gap structure recently found by VLA is $\gtrsim 0.4 M_J$. We also argue that the spiral structure, if observed, can be used to estimate the lower limit of the disk aspect ratio and the planet mass.

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

Radial Transport of Large-Scale Magnetic Fields in Accretion Disks. I. Steady Solutions and an Upper Limit on the Vertical Field Strength

Large-scale magnetic fields are key ingredients of magnetically driven disk accretion. We study how large-scale poloidal fields evolve in accretion disks, with the primary aim of quantifying the viability of magnetic accretion mechanisms in protoplanetary disks. We employ a kinematic mean-field model for poloidal field transport and focus on steady states where inward advection of a field balances with outward diffusion due to effective resistivities. We analytically derive the steady-state radial distribution of poloidal fields in highly conducting accretion disks. The analytic solution reveals an upper limit on the strength of large-scale vertical fields attainable in steady states. Any excess poloidal field will be diffused away within a finite time, and we demonstrate this with time-dependent numerical calculations of the mean-field equations. We apply this upper limit to large-scale vertical fields threading protoplanetary disks. We find that the maximum attainable strength is about 0.1 G at 1 AU, and about 1 mG at 10 AU from the central star. When combined with recent magnetic accretion models, the maximum field strength translates into the maximum steady-state accretion rate of $\sim 10^{-7} M_\odot {\rm yr}^{-1}$, in agreement with observations. We also find that the maximum field strength is ~ 1 kG at the surface of the central star provided that the disk extends down to the stellar surface. This implies that any excess stellar poloidal field of strength >~ kG can be transported to the surrounding disk. This might in part resolve the magnetic flux problem in star formation.

preprint2014arXiv

Radial Transport of Large-Scale Magnetic Fields in Accretion Disks. II. Relaxation to Steady States

We study the time evolution of a large-scale magnetic flux threading an accretion disk. Induction equation of the mean poloidal field is solved under the standard viscous disk model. Magnetic flux evolution is controlled by the two timescales: One is the timescale of the inward advection of the magnetic flux, tau_{adv}. This is induced by the dragging of the flux by the accreting gas. The other is the outward diffusion timescale of the magnetic flux tau_{dif}. We consider diffusion due to the Ohmic resistivity. These timescales can be significantly different from the disk viscous timescale tau_{disk}. The behaviors of the magnetic flux evolution is quite different depending on the magnitude relationship of the timescales tau_{adv}, τ_{dif}, and tau_{disk}. The most interesting phenomena occurs when tau_{adv} << tau_{dif}, tau_{disk}. In such a case, the magnetic flux distribution approaches a quasi-steady profile much faster than the viscous evolution of the gas disk, and also the magnetic flux has been tightly bundled to the inner part of the disk. In the inner part, although the poloidal magnetic field becomes much stronger than the interstellar magnetic field, the field strength is limited to the maximum value that is analytically given by our previous work (Okuzumi et al. 2014, ApJ, 785, 127). We also find a condition for that the initial large magnetic flux, which is a fossil of the magnetic field dragging during the early phase of star formation, survives for a duration in which significant gas disk evolution proceeds.

preprint2014arXiv

Rotational Instability in the Outer Region of Protoplanetary Disks

We analytically calculate the marginally stable surface density profile for rotational instability of protoplanetary disks. The derived profile can be utilized for considering the region in a rotating disk where radial pressure gradient force is comparable to the gravitational force, such as an inner edge, steep gaps or bumps and an outer region of the disk. In this paper we especially focus on the rotational instability in the outer region of disks. We find an protoplanetary disk with a surface density profile of similarity solution becomes rotationally unstable at a certain radius, depending on its temperature profile and a mass of the central star. If the temperature is relatively low and the mass of the central star is high, disks have rotationally stable similarity profiles. Otherwise, deviation from the similarity profiles of surface density could be observable, using facilities with high sensitivity, such as ALMA.

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.

preprint2012arXiv

Minimum Dust Abundances for Planetesimal Formation via Secular Gravitational Instabilities

We estimate minimum dust abundances required for secular gravitational instability (SGI) to operate at the midplane dust layer of protoplanetary disks. For SGI to be a viable process, the growth time of the instability T_grow must be shorter than the radial drift time of the dust T_drift. The growth time depends on the turbulent diffusion parameter alpha, because the modes with short wavelengths are stabilized by turbulent diffusion. Assuming that turbulence is excited via the Kelvin-Helmholtz or streaming instabilities in the dust layer, and that its strength is controlled by the energy supply rate from dust accretion, we estimate the diffusion parameter and the growth time of the instability. The condition T_grow < T_drift requires that the dust abundance must be greater than a critical abundance Z_min, which is a function of the Toomre parameter Q_g and aspect ratio h_g / r of the gas disk. For a wide range of parameter space, the required dust abundance is less than 0.1. A slight increase in dust abundance opens a possible route for the dust to directly collapse to planetesimals.

preprint2011arXiv

Electrostatic Barrier against Dust Growth in Protoplanetary Disks. I. Classifying the Evolution of Size Distribution

Collisional growth of submicron-sized dust grains into macroscopic aggregates is the first step of planet formation in protoplanetary disks. These grains are expected to carry nonzero negative charges in the weakly ionized disks, but its effect on their collisional growth has not been fully understood so far. In this paper, we investigate how the charging affects the evolution of the dust size distribution properly taking into account the charging mechanism in a weakly ionized gas as well as porosity evolution through low-energy collisions. To clarify the role of the size distribution, we divide our analysis into two steps. First, we analyze the collisional growth of charged aggregates assuming a monodisperse (i.e., narrow) size distribution. We show that the monodisperse growth stalls due to the electrostatic repulsion when a certain condition is met, as is already expected in the previous work. Second, we numerically simulate dust coagulation using Smoluchowski's method to see how the outcome changes when the size distribution is allowed to freely evolve. We find that, under certain conditions, the dust undergoes bimodal growth where only a limited number of aggregates continue to grow carrying the major part of the dust mass in the system. This occurs because remaining small aggregates efficiently sweep up free electrons to prevent the larger aggregates from being strongly charged. We obtain a set of simple criteria that allows us to predict how the size distribution evolves for a given condition. In Paper II (arXiv:1009.3101), we apply these criteria to dust growth in protoplanetary disks.

preprint2011arXiv

Electrostatic Barrier against Dust Growth in Protoplanetary Disks. II. Measuring the Size of the "Frozen" Zone

Coagulation of submicron-sized dust grains into porous aggregates is the initial step of dust evolution in protoplanetary disks. Recently, it has been pointed out that negative charging of dust in the weakly ionized disks could significantly slow down the coagulation process. In this paper, we apply the growth criteria obtained in Paper I to finding out a location ("frozen" zone) where the charging stalls dust growth at the fractal growth stage. For low-turbulence disks, we find that the frozen zone can cover the major part of the disks at a few to 100 AU from the central star. The maximum mass of the aggregates is approximately 10^{-7} g at 1 AU and as small as a few monomer masses at 100 AU. Strong turbulence can significantly reduce the size of the frozen zone, but such turbulence will cause the fragmentation of macroscopic aggregates at later stages. We examine a possibility that complete freezeout of dust evolution in low-turbulence disks could be prevented by global transport of dust in the disks. Our simple estimation shows that global dust transport can lead to the supply of macroscopic aggregates and the removal of frozen aggregates on a timescale of 10^6 yr. This overturns the usual understanding that tiny dust particles get depleted on much shorter timescales unless collisional fragmentation is effective. The frozen zone together with global dust transport might explain "slow" (\sim 10^6 yr) dust evolution suggested by infrared observation of T Tauri stars and by radioactive dating of chondrites.

preprint2011arXiv

Induced Turbulence and the Density Structure of the Dust Layer in a Protoplanetary Disk

We study the turbulence induced in the dust layer of a protoplanetary disk based on the energetics of dust accretion due to gas drag. We estimate turbulence strength from the energy supplied by dust accretion, using the radial drift velocity of the dust particles in a laminar disk. Our estimate of the turbulence strength agrees with previous analytical and numerical research on the turbulence induced by Kelvin-Helmholtz and/or streaming instabilities for particles whose stopping time is less than the Keplerian time. For such small particles, the strongest turbulence is expected to occur when the dust-to-gas ratio of the disk is ~C_eff^(1/2) (h_g / r) ~ 10^(-2), where C_eff ~ 0.2 represents the energy supply efficiency to turbulence and h_g / r ~ 5 x 10^(-2) is the aspect ratio of the gas disk. The maximum viscosity parameter is alpha_max ~ C_eff T_s (h_g / r)^2 ~ 10^(-4) T_s, where T_s (<1) is the non-dimensional stopping time of the dust particles. Modification in the dust-to-gas ratio from the standard value, 10^(-2), by any process, results in weaker turbulence and a thinner dust layer, and consequently may accelerate the growth process of the dust particles.

preprint2011arXiv

On the Interaction between a Protoplanetary Disk and a Planet in an Eccentric Orbit: Application of Dynamical Friction

We present a new analytic approach to the disk-planet interaction that is especially useful for planets with eccentricity larger than the disk aspect ratio. We make use of the dynamical friction formula to calculate the force exerted on the planet by the disk, and the force is averaged over the period of the planet. The resulting migration and eccentricity damping timescale agrees very well with the previous works in which the planet eccentricity is moderately larger than the disk aspect ratio. The advantage of this approach is that it is possible to apply this formulation to arbitrary large eccentricity. We have found that the timescale of the orbital evolution depends largely on the adopted disk model in the case of highly eccentric planets. We discuss the possible implication of our results to the theory of planet formation.

preprint2003arXiv

The 1.2 Millimeter Image of the beta Pictoris Disk

We present millimeter imaging observations in the 1200 micron continuum of the disk around beta Pictoris. With the 25 arcsec beam, the beta Pic disk is unresolved perpendicularly to the disk plane (< 10 arcsec), but slightly resolved in the northeast-southwest direction (26 arcsec). Peak emission is observed at the stellar position. A secondary maximum is found 1000 AU along the disk plane in the southwest, which does not positionally coincide with a similar feature reported earlier at 850 micron. Arguments are presented which could be seen in support of the reality of these features. The observed submm/mm emission is consistent with thermal emission from dust grains, which are significantly larger than those generally found in the interstellar medium, including mm-size particles, and thus more reminiscent of the dust observed in protostellar disks. Modelling the observed scattered light in the visible and the emission in the submm/mm provides evidence for the particles dominating the scattering in the visible/NIR and those primarily responsible for the thermal emission at longer wavelengths belonging to different populations.

preprint2002arXiv

Radial Flow of Dust Particles in Accretion Disks

We study the radial migration of dust particles in accreting protostellar disks analogous to the primordial solar nebula. This study takes account of the two dimensional (radial and normal) structure of the disk gas, including the effects of the variation in the gas velocity as a function of distance from the midplane. It is shown that the dust component of disks accretes slower than the gas component. At high altitude from the disk midplane, the gas rotates faster than particles because of the inward pressure gradient force, and its drag force causes particles to move outward in the radial direction. Viscous torque induces the gas within a scale height from the disk midplane to flow outward, carrying small (size < 100 micron at 10 AU) particles with it. Only particles at intermediate altitude or with sufficiently large sizes (> 1 mm at 10 AU) move inward. When the particles' radial velocities are averaged over the entire vertical direction, particles have a net inward flux. At large distances from the central star, particles migrate inward with a velocity much faster than the gas accretion velocity. However, their inward velocity is reduced below that of the gas in the inner regions of the disk. The rate of velocity decrease is a function of the particles' size. While larger particles retain fast accretion velocity until they approach closer to the star, 10 micron particles have slower velocity than the gas in the most part of the disk (r < 100 AU). This differential migration of particles causes the size fractionation. Dust disks composed mostly of small particles (size < 10 micron) accrete slower than gas disks, resulting in the increase in the dust-gas ratio during the gas accretion phase.