Source author record

Takayuki Muto

Takayuki Muto 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

45works
7topics
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

45 published item(s)

preprint2026arXiv

Bridging the gap: consistent modeling of protoplanetary disk heating and gap formation by planet-induced spiral shocks

A giant planet embedded in a protoplanetary disk excites spiral density waves, which steepen into shocks as they propagate away from the planet. These shocks lead to secular disk heating and gap opening, both of which can have important implications for the evolution of solids near the planet. To date, these two effects have largely been modeled independently. In this study, we present a self-consistent model that unifies these processes by linking shock heating and angular momentum deposition through the entropy jumps across the spiral shocks. We show that this model accurately reproduces the temperature and surface density profiles around the planet's orbit, as obtained from two-dimensional hydrodynamic simulations with standard $α$ viscosity and $β$ thermal relaxation prescriptions. Furthermore, by incorporating an empirically derived scaling law for the radial distribution of the entropy jump, we construct a fully analytic model that self-consistently predicts the temperature and surface density structures of disks hosting a giant planet. This work represents a first step toward understanding how a giant planet forming in the inner disk region influences the distribution and composition of second-generation planets and planetesimals in its vicinity.

preprint2024arXiv

Shadowing in the protoplanetary disk of ZZ Tau IRS with HST

An inner component misaligned from an outer component in a protoplanetary disk can result in the former casting shadows on the latter. We present a new instance of shadowing on the outer disk around a very low mass star, ZZ~Tau~IRS. Through the analysis of near-infrared (NIR) archival data at $λ=1.6$~$μ$m acquired with the Wide Field Camera 3 on the Hubble Space Telescope, we identified brightness asymmetries in the top and bottom halves of the highly inclined outer disk, separated by a dark lane. The brighter sides in the top and bottom halves are on the opposite sides, which we attributed to shadows cast by a misaligned inner disk. Radiative transfer modeling of the system with a misaligned angle of 15~deg between the inner and outer disks well reproduced the observations. Additionally, we found an elevated brightness temperature of $^{12}$CO~(3-2) at $r\sim30$~au on the brighter side in NIR wavelengths in the top half disk, which can be explained by the shadowing effect too. While the origin of the misaligned inner disk remains unclear, future monitoring observations to search for temporal variations in brightness asymmetries will likely provide useful clues.

preprint2022arXiv

A spatially-resolved large cavity of the J0337 protoplanetary disk in Perseus

We present Keck/NIRC2 $K_{\rm p}L_{\rm p}$ high-contrast imaging observations of a J0337 protoplanetary disk. The data discover the spatially-resolved large cavity, which is the second report among protoplanetary disks in the Perseus star forming region after the LkH$α$~330 system. Our data and forward modeling using RADMC-3D suggests $\sim80$~au for the cavity radius. There is discrepancy between J0337's SED and the modeled SED at $\sim10\micron$ and this suggests an unseen inner disk. We also searched for companions around J0337 but did not detect any companion candidates at separations between $0\farcs1$ and $2\farcs5$. The $L_{\rm p}$-band detection limit corresponds to $\sim20 M_{\rm Jup}$ at 60~au, $\sim9-10 M_{\rm Jup}$ at 90~au, and $\sim3 M_{\rm Jup}$ at $>120$~au. Compared with other young systems with large cavities such as PDS~70 and RX~J1604, multiple Jovian planets, a single eccentric Jovian planet, or a massive brown-dwarf at an inner separation could exist within the cavity.

preprint2022arXiv

ALMA High-resolution Multiband Analysis for the Protoplanetary Disk around TW Hya

We present a high-resolution (2.5 au) multiband analysis of the protoplanetary disk around TW Hya using ALMA long baseline data at Bands 3, 4, 6, and 7. We aim to reconstruct a high-sensitivity millimeter continuum image and revisit the spectral index distribution. The imaging is performed by combining new ALMA data at Bands 4 and 6 with available archive data. Two methods are employed to reconstruct the images; multi-frequency synthesis (MFS) and the fiducial image-oriented method, where each band is imaged separately and the frequency dependence is fitted pixel by pixel. We find that the MFS imaging with the second order of Taylor expansion can reproduce the frequency dependence of the continuum emission between Bands 3 and 7 in a manner consistent with previous studies and is a reasonable method to reconstruct the spectral index map. The image-oriented method provides a spectral index map consistent with the MFS imaging, but with a two times lower resolution. Mock observations of an intensity model were conducted to validate the images from the two methods. We find that the MFS imaging provides a high-resolution spectral index distribution with an uncertainty of $<10$~\%. Using the submillimeter spectrum reproduced from our MFS images, we directly calculated the optical depth, power-law index of the dust opacity coefficient ($β$), and dust temperature. The derived parameters are consistent with previous works, and the enhancement of $β$ within the intensity gaps is also confirmed, supporting a deficit of millimeter-sized grains within the gaps.

preprint2022arXiv

Dippers from TESS Full-frame Images. II. Spectroscopic Characterization of Four Young Dippers

Photometric monitoring by the Transiting Exoplanet Survey Satellite (TESS) has discovered not only periodic signals by transiting exoplanets but also episodic or quasi-periodic dimming around young stellar objects. The mechanisms of the dimming of these objects, so-called "dippers", are thought to be related to the property of the accretion or the structure of protoplanetary disks especially in regions close to the host star. Recently, we have created the catalog of dippers from the one year of TESS Full Frame Image (FFI) data. In this paper, we report spectral features of four newly found dippers in that catalog and show that they potentially shed light on the dimming mechanisms. We found that all of the targets exhibit the H alpha emission line, which is an indicator of an ongoing accretion. Based on its line profiles and/or their variability, we characterized the properties of the disks of each source, which can support the dimming mechanisms due to a dusty disk wind or an accretion warp. Also, we found an interesting dipper (TIC 317873721), "close-in binary dipper," showing the complex variability of the line profile and the large radial velocity variation. Since the dimming intervals are similar to the orbital period of the binary, we suggest that the dips are caused by dust in the warp of accretion from a circumbinary disk onto stars. Such a close-in ($<$ 0.1 au) binary dipper is rarely reported so far, further investigation will reveal the new aspect of the disk evolution and planetary formation.

preprint2022arXiv

Gap Opening and Inner Disk Structure in the Strongly Accreting Transition Disk of DM Tau

Large inner dust gaps in transition disks are frequently posited as evidence of giant planets sculpting gas and dust in the disk, or the opening of a gap by photoevaporative winds. Although the former hypothesis is strongly supported by the observations of planets and deep depletions in gas within the gap some disks, many T Tauri stars hosting transition disks accrete at rates typical for an undepleted disk, raising the question of how gap opening occurs in these objects. We thus present an analysis of the structure of the transition disk around the T Tauri star DM Tau, which is strongly accreting ($\sim 10^{-8.3}~\mathrm{M}_\odot~ \mathrm{yr}^{-1}$) and turbulent ($α=0.078 \pm 0.02$). Using the DALI thermochemical code, we fit disk models to simultaneously reproduce the accretion rate, high level of turbulence, the gas traced by ALMA band 6 observations of $^{12}$CO, $^{13}$CO, and C$^{18}$O J=2--1 lines, and the observed dust emission from the mm continuum and spectral energy distribution. We find a shallow depletion in gas surface density of $\sim 10$ relative to the outer disk and a gas rich inner disk is consistent with the observations. The planet mass of $<1$ M$_\mathrm{Jup}$ implied by the gap depth is in tension with predictions for dust trapping in a highly viscous disk, which requires a more massive planet of of $\sim10$M$_\mathrm{Jup}$. Photoevaporative models including a dead zone can qualitatively reproduce some features of the DM Tau disk, but still struggle to explain the high accretion rates and the observed mm continuum flux.

preprint2022arXiv

Images of Embedded Jovian Planet Formation At A Wide Separation Around AB Aurigae

Direct images of protoplanets embedded in disks around infant stars provide the key to understanding the formation of gas giant planets like Jupiter. Using the Subaru Telescope and Hubble Space Telescope, we find evidence for a jovian protoplanet around AB Aurigae orbiting at a wide projected separation (93 au), likely responsible for multiple planet-induced features in the disk. Its emission is reproducible as reprocessed radiation from an embedded protoplanet. We also identify two structures located at 430-580 au that are candidate sites of planet formation. These data reveal planet formation in the embedded phase and a protoplanet discovery at wide, > 50 au separations characteristic of most imaged exoplanets. With at least one clump-like protoplanet and multiple spiral arms, the AB Aur system may also provide the evidence for a long-considered alternative to the canonical model for Jupiter's formation: disk (gravitational) instability.

preprint2020arXiv

A new and simple prescription for planet orbital migration and eccentricity damping by planet-disc interactions based on dynamical friction

During planet formation gravitational interaction between a planetary embryo and the protoplanetary gas disc causes orbital migration of the planetary embryo, which plays an important role in shaping the final planetary system. While migration sometimes occurs in the supersonic regime, wherein the relative velocity between the planetary embryo and the gas is higher than the sound speed, migration prescriptions proposed thus far describing the planet-disc interaction force and the timescales of orbital change in the supersonic regime are inconsistent with one another. Here we discuss the details of existing prescriptions in the literature and derive a new simple and intuitive formulation for planet-disc interactions based on dynamical friction that can be applied in both supersonic and subsonic cases. While the existing prescriptions assume particular disc models, ours include the explicit dependence on the disc parameters; hence it can be applied to discs with any radial surface density and temperature dependence (except for the local variations with radial scales less than the disc scale height). Our prescription will reduce the uncertainty originating from different literature formulations of planet migration and will be an important tool to study planet accretion processes, especially when studying the formation of close-in low-mass planets that are commonly found in exoplanetary systems.

preprint2020arXiv

ALMA Observations of the Inner Cavity in the Protoplanetary Disk around Sz 84

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of a protoplanetary disk around the T Tauri star Sz~84 and analyses of the structures of the inner cavity in the central region of the dust disk. Sz~84's spectral energy distribution (SED) has been known to exhibit negligible infrared excess at $λ\lesssim$10~$μ$m due to the disk's cavity structure. Analyses of the observed visibilities of dust continuum at 1.3~mm and the SED indicate that the size of the cavity in the disk of large (millimeter size) dust grains is 8~au in radius and that in the disk of small (sub-micron size) dust grains is 60~au in radius. Furthermore, from the SED analyses, we estimate that the upper limit mass of small dust grains at $r<$60~au is less than $\sim$10$^{-3}$~$M_{\rm earth}$, which is $\lesssim$0.01~\% of the total (small~$+$~large) dust mass at $r<$60~au. These results suggest that large dust grains are dominant at $r<$60~au, implying that dust grains efficiently grow with less efficient fragmentation in this region, potentially due to weak turbulence and/or stickier dust grains. The balance of grain growth and dust fragmentation is an important factor for determining the size of large dust grains in protoplanetary disks, and thus Sz~84 could serve as a good testbed for investigations of grain growth in such disks.

preprint2020arXiv

GW Ori: Interactions Between a Triple-star System and its Circumtriple Disk in Action

GW Ori is a hierarchical triple system which has a rare circumtriple disk. We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of 1.3 mm dust continuum and 12CO J=2-1 molecular gas emission of the disk. For the first time, we identify three dust rings in the disk at ~46, 188, and 338 AU, with estimated dust mass of ~70-250 Earth masses, respectively. To our knowledge, the outer ring in GW Ori is the largest dust ring ever found in protoplanetary disks. We use visibility modelling of dust continuum to show that the disk has misaligned parts and the innermost dust ring is eccentric. The disk misalignment is also suggested by the CO kinematics modelling. We interpret these substructures as evidence of ongoing dynamical interactions between the triple stars and the circumtriple disk.

preprint2020arXiv

Model of a gap formed by a planet with fast inward migration

A planet is formed within a protoplanetary disk. Recent observations have revealed substructures such as gaps and rings, which may indicate forming planets within the disk. Due to disk--planet interaction, the planet migrates within the disk, which can affect a shape of the planet-induced gap. In this paper, we investigate effects of fast inward migration of the planet on the gap shape, by carrying out hydrodynamic simulations. We found that when the migration timescale is shorter than the timescale of the gap-opening, the orbital radius is shifted inward as compared to the radial location of the gap. We also found a scaling relation between the radial shift of the locations of the planet and the gap as a function of the ratio of the timescale of the migration and gap-opening. Our scaling relation also enables us to constrain the gas surface density and the viscosity when the gap and the planet are observed. Moreover, we also found the scaling relation between the location of the secondary gap and the aspect ratio. By combining the radial shift and the secondary gap, we may constrain the physical condition of the planet formation and how the planet evolves in the protoplanetary disk, from the observational morphology.

preprint2020arXiv

Near-Infrared Imaging of a Spiral in the CQ Tau Disk

We present $L^\prime$-band Keck/NIRC2 imaging and $H$-band Subaru/AO188+HiCIAO polarimetric observations of CQ Tau disk with a new spiral arm. Apart from the spiral feature our observations could not detect any companion candidates. We traced the spiral feature from the $r^2$-scaled HiCIAO polarimetric intensity image and the fitted result is used for forward modeling to reproduce the ADI-reduced NIRC2 image. We estimated the original surface brightness after throughput correction in $L^\prime$-band to be $\sim126$ mJy/arcsec$^2$ at most. We suggest that the grain temperature of the spiral may be heated up to $\sim$200 K in order to explain both of the $H$- and $L^{\prime}$-bands results. The $H$-band emission at the location of the spiral originates from the scattering from the disk surface while both scattering and thermal emission may contribute to the $L^{\prime}$-band emission. If the central star is only the light source of scattered light, the spiral emission at $L^\prime$-band should be thermal emission. If an inner disk also acts as the light source, the scattered light and the thermal emission may equally contribute to the $L^\prime$-band spiral structure.

preprint2020arXiv

SCExAO/CHARIS High-Contrast Imaging of Spirals and Darkening Features in the HD 34700 A Protoplanetary Disk

We present Subaru/SCExAO+CHARIS broadband ($JHK$-band) integral field spectroscopy of HD 34700 A. CHARIS data recover HD 34700 A's disk ring and confirm multiple spirals discovered in Monnier et al. (2019). We set limits on substellar companions of $\sim12\ M_{\rm Jup}$ at $0\farcs3$ (in the ring gap) and $\sim5\ M_{\rm Jup}$ at $0\farcs75$ (outside the ring). The data reveal darkening effects on the ring and spiral, although we do not identify the origin of each feature such as shadows or physical features related to the outer spirals. Geometric albedoes converted from the surface brightness suggests a higher scale height and/or prominently abundant sub-micron dust at position angle between $\sim45^\circ$ and $90^\circ$. Spiral fitting resulted in very large pitch angles ($\sim30-50^\circ$) and a stellar flyby of HD 34700 B or infall from a possible envelope is perhaps a reasonable scenario to explain the large pitch angles.

preprint2020arXiv

Super-resolution Imaging of the Protoplanetary Disk HD 142527 Using Sparse Modeling

With an emphasis on improving the fidelity even in super-resolution regimes, new imaging techniques have been intensively developed over the last several years, which may provide substantial improvements to the interferometric observation of protoplanetary disks. In this study, sparse modeling (SpM) is applied for the first time to observational data sets taken by the Atacama Large Millimeter/submillimeter Array (ALMA). The two data sets used in this study were taken independently using different array configurations at Band 7 (330 GHz), targeting the protoplanetary disk around HD 142527; one in the shorter-baseline array configuration (~ 430 m), and the other in the longer-baseline array configuration (~ 1570 m). The image resolutions reconstructed from the two data sets are different by a factor of ~ 3. We confirm that the previously known disk structures appear on the images produced by both SpM and CLEAN at the standard beam size. The image reconstructed from the shorter-baseline data using the SpM matches that obtained with the longer-baseline data using CLEAN, achieving a super-resolution image from which a structure finer than the beam size can be reproduced. Our results demonstrate that on-going intensive development in the SpM imaging technique is beneficial to imaging with ALMA.

preprint2019arXiv

The detection of dust gap-ring structure in the outer region of the CR Cha protoplanetary disk

We observe the dust continuum at 225 GHz and CO isotopologue (12CO, 13CO, and C18O) J=2-1 emission lines toward the CR Cha protoplanetary disk using the Atacama Large Millimeter/Submillimeter Array (ALMA). The dust continuum image shows a dust gap-ring structure in the outer region of the dust disk. A faint dust ring is also detected around 120 au beyond the dust gap. The CO isotopologue lines indicate that the gas disk is more extended than the dust disk. The peak brightness temperature of the 13CO line shows a small bump around 130 au while 12CO and C18O lines do not show. We investigate two possible mechanisms for reproducing the observed dust gap-ring structure and a gas temperature bump. First, the observed gap structure can be opened by a Jupiter mass planet using the relation between the planet mass and the gap depth and width. Meanwhile, the radiative transfer calculations based on the observed dust surface density profile show that the observed dust ring could be formed by dust accumulation at the gas temperature bump, that is, the gas pressure bump produced beyond the outer edge of the dust disk.

preprint2019arXiv

Unveiling Dust Aggregate Structure in Protoplanetary Disks by Millimeter-wave Scattering Polarization

Dust coagulation in a protoplanetary disk is the first step of planetesimal formation. However, a pathway from dust aggregates to planetesimals remains unclear. Both numerical simulations and laboratory experiments have suggested the importance of dust structure in planetesimal formation, but it is not well constrained by observations. We study how dust structure and porosity alters polarimetric images at millimeter wavelength by performing 3D radiative transfer simulations. Aggregates with different porosity and fractal dimension are considered. As a result, we find that dust aggregates with lower porosity and/or higher fractal dimension are favorable to explain observed millimeter-wave scattering polarization of disks. Aggregates with extremely high porosity fail to explain the observations. In addition, we also show that particles with moderate porosity show weak wavelength dependence of scattering polarization, indicating that multi-wavelength polarimetry is useful to constrain dust porosity. Finally, we discuss implications for dust evolution and planetesimal formation in disks.

preprint2016arXiv

ALMA Observations of a Gap and a Ring in the Protoplanetary Disk around TW Hya

We report the first detection of a gap and a ring in 336 GHz dust continuum emission from the protoplanetary disk around TW Hya, using the Atacama Large Millimeter/Submillimeter Array (ALMA). The gap and ring are located at around 25 and 41 au from the central star, respectively, and are associated with the CO snow line at ~30 au. The gap has a radial width of less than 15 au and a mass deficit of more than 23%, taking into account that the observations are limited to an angular resolution of ~15 au. In addition, the 13CO and C18O J = 3 - 2 lines show a decrement in CO line emission throughout the disk, down to ~10 au, indicating a freeze-out of gas-phase CO onto grain surfaces and possible subsequent surface reactions to form larger molecules. The observed gap could be caused by gravitational interaction between the disk gas and a planet with a mass less than super-Neptune (2M_{Neptune}), or could be the result of the destruction of large dust aggregates due to the sintering of CO ice.

preprint2016arXiv

Extreme Asymmetry in the Disk of V1247 Ori

We present the first near-infrared scattered-light detection of the transitional disk around V1247 Ori, which was obtained using high-resolution polarimetric differential imaging observations with Subaru/HiCIAO. Our imaging in the H band reveals the disk morphology at separations of ~0.14"-0.86" (54-330 au) from the central star. The polarized intensity (PI) image shows a remarkable arc-like structure toward the southeast of the star, whereas the fainter northwest region does not exhibit any notable features. The shape of the arm is consistent with an arc of 0.28" $\pm$ 0.09" in radius (108 au from the star), although the possibility of a spiral arm with a small pitch angle cannot be excluded. V1247 Ori features an exceptionally large azimuthal contrast in scattered, polarized light; the radial peak of the southeastern arc is about three times brighter than the northwestern disk measured at the same distance from the star. Combined with the previous indication of an inhomogeneous density distribution in the gap at $\lesssim$46 au, the notable asymmetry in the outer disk suggests the presence of unseen companions and/or planet-forming processes ongoing in the arc.

preprint2016arXiv

Grain size constraints on HL Tau with polarization signature

The millimeter-wave polarization of the protoplanetary disk around HL Tau has been interpreted as the emission from elongated dust grains aligned with the magnetic field in the disk. However, the self-scattering of thermal dust emission may also explain the observed millimeter-wave polarization. In this paper, we report a modeling of the millimeter-wave polarization of the HL Tau disk with the self-polarization. Dust grains are assumed to be spherical and to have a power-law size distribution. We change the maximum grain size with a fixed dust composition in a fixed disk model to find the grain size to reproduce the observed signature. We find that the direction of the polarization vectors and the polarization degree can be explained with the self-scattering. Moreover, the polarization degree can be explained only if the maximum grain size is $\sim 150 {\rm~μm}$. The obtained grain size from the polarization is different from a size of millimeter or larger that has been previously expected from the spectral index of the dust opacity coefficient if the emission is optically thin. We discuss that porous dust aggregates may solve the inconsistency of the maximum grain size between the two constraints.

preprint2016arXiv

Investigating the Relation between CO (3-2) and Far Infrared Luminosities for Nearby Merging Galaxies Using ASTE

We present the new single dish CO (3-2) emission data obtained toward 19 early stage and 7 late stage nearby merging galaxies using the Atacama Submillimeter Telescope Experiment (ASTE). Combining with the single dish and interferometric data of galaxies observed in previous studies, we investigate the relation between the CO (3-2) luminosity (L'CO(3-2)) and the far Infrared luminosity (LFIR) in a sample of 29 early stage and 31 late stage merging galaxies, and 28 nearby isolated spiral galaxies. We find that normal isolated spiral galaxies and merging galaxies have different slopes (alpha) in the log L'CO(3-2) - log LFIR plane (alpha ~ 0.79 for spirals and ~ 1.12 for mergers). The large slope (alpha > 1) for merging galaxies can be interpreted as an evidence for increasing Star Formation Efficiency (SFE=LFIR/L'CO(3-2)) as a function of LFIR. Comparing our results with sub-kpc scale local star formation and global star-burst activity in the high-z Universe, we find deviations from the linear relationship in the log L'CO(3-2) - log LFIR plane for the late stage mergers and high-z star forming galaxies. Finally, we find that the average SFE gradually increases from isolated galaxies, merging galaxies, and to high-z submillimeter galaxies / quasi-stellar objects (SMGs/QSOs). By comparing our findings with the results from numerical simulations, we suggest; (1) inefficient star-bursts triggered by disk-wide dense clumps occur in the early stage of interaction and (2) efficient star-bursts triggered by central concentration of gas occur in the final stage. A systematic high spatial resolution survey of diffuse and dense gas tracers is a key to confirm this scenario.

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.

preprint2016arXiv

Submillimeter polarization observation of the protoplanetary disk around HD 142527

We present the polarization observations toward the circumstellar disk around HD 142527 by using Atacama Large Millimeter/submillimeter Array (ALMA) at the frequency of 343 GHz. The beam size is $0.51 " \times 0.44 "$, which corresponds to the spatial resolution of $\sim$ 71 $\times$ 62 AU. The polarized intensity displays a ring-like structure with a peak located on the east side with a polarization fraction of $P= 3.26 \pm 0.02$ %, which is different from the peak of the continuum emission from the northeast region. The polarized intensity is significantly weaker at the peak of the continuum where $P= 0.220 \pm 0.010$ %. The polarization vectors are in the radial direction in the main ring of the polarized intensity, while there are two regions outside at the northwest and northeast areas where the vectors are in the azimuthal direction. If the polarization vectors represent the magnetic field morphology, the polarization vectors indicate the toroidal magnetic field configuration on the main ring and the poloidal fields outside. On the other hand, the flip of the polarization vectors is predicted by the self-scattering of thermal dust emission due to the change of the direction of thermal radiation flux. Therefore, we conclude that self-scattering of thermal dust emission plays a major role in producing polarization at millimeter wavelengths in this protoplanetary disk. Also, this puts a constraint on the maximum grain size to be approximately 150 ${\rm μm}$ if we assume compact spherical dust grains.

preprint2016arXiv

The Matryoshka Disk: Keck/NIRC2 Discovery of a Solar System-Scale, Radially Segregated Residual Protoplanetary Disk Around HD 141569A

Using Keck/NIRC2 $L^\prime$ (3.78 $μm$) data, we report the direct imaging discovery of a scattered light-resolved, solar system-scale residual protoplanetary disk around the young A-type star HD 141569A, interior to and concentric with the two ring-like structures at wider separations. The disk is resolved down to $\sim$ 0\farcs{}25 and appears as an arc-like rim with attached hook-like features. It is located at an angular separation intermediate between that of warm CO gas identified from spatially-resolved mid-infrared spectroscopy and diffuse dust emission recently discovered with the \textit{Hubble Space Telescope}. The inner disk has a radius of $\sim$ 39 AU, a position angle consistent with north-up, an inclination of $i$ $\sim$ 56$^{o}$, and has a center offset from the star. Forward-modeling of the disk favors a thick torus-like emission sharply truncated at separations beyond the torus' photocenter and heavily depleted at smaller separations. In particular, the best-fit density power law for the dust suggests that the inner disk dust and gas (as probed by CO) are radially segregated, a feature consistent with the dust trapping mechanism inferred from observations of "canonical" transitional disks. However, the inner disk component may instead be explained by radiation pressure-induced migration in optically-thin conditions, in contrast to the two stellar companion/planet-influenced ring-like structures at wider separations. HD 141569A's circumstellar environment --- with three nested, gapped, concentric dust populations --- is an excellent laboratory for understanding the relationship between planet formation and the evolution of both dust grains and disk architecture.

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

Millimeter-wave polarization of protoplanetary disks due to dust scattering

We present a new method to constrain the grain size in protoplanetary disks with polarization observations at millimeter wavelengths. If dust grains are grown to the size comparable to the wavelengths, the dust grains are expected to have a large scattering opacity and thus the continuum emission is expected to be polarized due to self-scattering. We perform 3D radiative transfer calculations to estimate the polarization degree for the protoplanetary disks having radial Gaussian-like dust surface density distributions, which have been recently discovered. The maximum grain size is set to be $100 {\rm~μm}$ and the observing wavelength to be 870 ${\rm μm}$. We find that the polarization degree is as high as 2.5 % with a subarcsec spatial resolution, which is likely to be detected with near-future ALMA observations. The emission is polarized due to scattering of anisotropic continuum emission. The map of the polarization degree shows a double peaked distribution and the polarization vectors are in the radial direction in the inner ring and in the azimuthal direction in the outer ring. We also find the wavelength dependence of the polarization degree: the polarization degree is the highest if dust grains have a maximum size of $a_{\rm max}\simλ/2π$, where $λ$ is the observing wavelength. Hence, multi-wave and spatially resolved polarization observations toward protoplanetary disks enable us to put a constraint on the grain size. The constraint on the grain size from polarization observations is independent of or may be even stronger than that from the opacity index.

preprint2015arXiv

Resolving the HD 100546 Protoplanetary System with the Gemini Planet Imager: Evidence for Multiple Forming, Accreting Planets

We report Gemini Planet Imager H band high-contrast imaging/integral field spectroscopy and polarimetry of the HD 100546, a 10 $Myr$-old early-type star recently confirmed to host a thermal infrared bright (super)jovian protoplanet at wide separation, HD 100546 b. We resolve the inner disk cavity in polarized light, recover the thermal-infrared (IR) bright arm, and identify one additional spiral arm. We easily recover HD 100546 b and show that much of its emission originates an unresolved, point source. HD 100546 b likely has extremely red infrared colors compared to field brown dwarfs, qualitatively similar to young cloudy superjovian planets, however, these colors may instead indicate that HD 100546 b is still accreting material from a circumplanetary disk. Additionally, we identify a second point source-like peak at $r_{proj}$ $\sim$ 14 AU, located just interior to or at inner disk wall consistent with being a 10--20 $M_{J}$ candidate second protoplanet-- "HD 100546 c" -- and lying within a weakly polarized region of the disk but along an extension of the thermal IR bright spiral arm. Alternatively, it is equally plausible that this feature is a weakly polarized but locally bright region of the inner disk wall. Astrometric monitoring of this feature over the next 2 years and emission line measurements could confirm its status as a protoplanet, rotating disk hot spot that is possibly a signpost of a protoplanet, or a stationary emission source from within the disk.

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

Recovery of the Candidate Protoplanet HD 100546 b with Gemini/NICI and Detection of Additional (Planet-Induced?) Disk Structure at Small Separations

We report the first independent, second-epoch (re-)detection of a directly-imaged protoplanet candidate. Using $L^\prime$ high-contrast imaging of HD 100546 taken with the Near-Infrared Coronagraph and Imager (NICI) on Gemini South, we recover `HD 100546 b' with a position and brightness consistent with the original VLT/NaCo detection from Quanz et al, although data obtained after 2013 will be required to decisively demonstrate common proper motion. HD 100546 b may be spatially resolved, up to $\approx$ 12-13 AU in diameter, and is embedded in a finger of thermal IR bright, polarized emission extending inwards to at least 0.3". Standard hot-start models imply a mass of $\approx$ 15 $M_{J}$. But if HD 100546 b is newly formed or made visible by a circumplanetary disk, both of which are plausible, its mass is significantly lower (e.g. 1--7 $M_{J}$). Additionally, we discover a thermal IR-bright disk feature, possibly a spiral density wave, at roughly the same angular separation as HD 100546 b but 90 degrees away. Our interpretation of this feature as a spiral arm is not decisive, but modeling analyses using spiral density wave theory implies a wave launching point exterior to $\approx$ 0.45" embedded within the visible disk structure: plausibly evidence for a second, hitherto unseen wide-separation planet. With one confirmed protoplanet candidate and evidence for 1--2 others, HD 100546 is an important evolutionary precursor to intermediate-mass stars with multiple super-jovian planets at moderate/wide separations like HR 8799.

preprint2014arXiv

SEEDS Adaptive Optics Imaging of the Asymmetric Transition Disk Oph IRS 48 in Scattered Light

We present the first resolved near infrared imagery of the transition disk Oph IRS 48 (WLY 2-48), which was recently observed with ALMA to have a strongly asymmetric sub-millimeter flux distribution. H-band polarized intensity images show a $\sim$60AU radius scattered light cavity with two pronounced arcs of emission, one from Northeast to Southeast and one smaller, fainter and more distant arc in the Northwest. K-band scattered light imagery reveals a similar morphology, but with a clear third arc along the Southwestern rim of the disk cavity. This arc meets the Northwestern arc at nearly a right angle, revealing the presence of a spiral arm or local surface brightness deficit in the disk, and explaining the East-West brightness asymmetry in the H-band data. We also present 0.8-5.4$μ$m IRTF SpeX spectra of this object, which allow us to constrain the spectral class to A0$\pm$1 and measure a low mass accretion rate of 10$^{-8.5}$M$_{\odot}$/yr, both consistent with previous estimates. We investigate a variety of reddening laws in order to fit the mutliwavelength SED of Oph IRS 48 and find a best fit consistent with a younger, higher luminosity star than previous estimates.

preprint2014arXiv

Surface Geometry of Protoplanetary Disks Inferred From Near-Infrared Imaging Polarimetry

We present a new method of analysis for determining the surface geometry of five protoplanetary disks observed with near-infrared imaging polarimetry using Subaru-HiCIAO. Using as inputs the observed distribution of polarized intensity (PI), disk inclination, assumed properties for dust scattering, and other reasonable approximations, we calculate a differential equation to derive the surface geometry. This equation is numerically integrated along the distance from the star at a given position angle. We show that, using these approximations, the local maxima in the PI distribution of spiral arms (SAO 206462, MWC 758) and rings (2MASS J16042165-2130284, PDS 70) are associated with local concave-up structures on the disk surface. We also show that the observed presence of an inner gap in scattered light still allows the possibility of a disk surface that is parallel to the light path from the star, or a disk that is shadowed by structures in the inner radii. Our analysis for rings does not show the presence of a vertical inner wall as often assumed in studies of disks with an inner gap. Finally, we summarize the implications of spiral and ring structures as potential signatures of ongoing planet formation.

preprint2014arXiv

Variability of Disk Emission in Pre-Main Sequence and Related Stars. III. Exploring Structural Changes in the Pre-transitional Disk in HD 169142

We present near-IR and far-UV observations of the pre-transitional (gapped) disk in HD 169142 using NASA's Infrared Telescope Facility and Hubble Space Telescope. The combination of our data along with existing data sets into the broadband spectral energy distribution reveals variability of up to 45% between ~1.5-10 μm over a maximum timescale of 10 years. All observations known to us separate into two distinct states corresponding to a high near-IR state in the pre-2000 epoch and a low state in the post-2000 epoch, indicating activity within the <1 AU region of the disk. Through analysis of the Pa β and Br γ lines in our data we derive a mass accretion rate in May 2013 of (1.5 - 2.7) x 10^-9 Msun/yr. We present a theoretical modeling analysis of the disk in HD 169142 using Monte-Carlo radiative transfer simulation software to explore the conditions and perhaps signs of planetary formation in our collection of 24 years of observations. We find that shifting the outer edge (r = 0.3 AU) of the inner disk by 0.05 AU toward the star (in simulation of accretion and/or sculpting by forming planets) successfully reproduces the shift in NIR flux. We establish that the ~40-70 AU dark ring imaged in the NIR by Quanz et al. (2013) and Momose et al. (2013) and at 7 mm by Osorio et al. (2014) may be reproduced with a 30% scaled density profile throughout the region, strengthening the link to this structure being dynamically cleared by one or more planetary mass bodies.

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

The Structure of Pre-transitional Protoplanetary Disks I: Radiative Transfer Modeling of the Disk+Cavity in the PDS 70 system

Through detailed radiative transfer modeling, we present a disk+cavity model to simultaneously explain both the SED and Subaru H-band polarized light imaging for the pre-transitional protoplanetary disk PDS 70. Particularly, we are able to match not only the radial dependence, but also the absolute scale, of the surface brightness of the scattered light. Our disk model has a cavity 65 AU in radius, which is heavily depleted of sub-micron-sized dust grains, and a small residual inner disk which produces a weak but still optically thick NIR excess in the SED. To explain the contrast of the cavity edge in the Subaru image, a factor of ~1000 depletion for the sub-micron-sized dust inside the cavity is required. The total dust mass of the disk may be on the order of 1e-4 M_sun, only weakly constrained due to the lack of long wavelength observations and the uncertainties in the dust model. The scale height of the sub-micron-sized dust is ~6 AU at the cavity edge, and the cavity wall is optically thick in the vertical direction at H-band. PDS 70 is not a member of the class of (pre-)transitional disks identified by Dong et al. (2012), whose members only show evidence of the cavity in the millimeter-sized dust but not the sub-micron-sized dust in resolved images. The two classes of (pre-)transitional disks may form through different mechanisms, or they may just be at different evolution stages in the disk clearing process.

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.

preprint2011arXiv

Resonant Trapping of Stars by Merging Massive Black Hole Binaries

A massive black hole binary might resonantly trap a star (e.g. a white dwarf) and the gas released by its tidal disruption might emit electromagnetic wave signals around the coalescence of the binary. With post-Newtonian equations of motion including gravitational radiation reaction, we numerically studied resonant trappings by black hole binaries with mass ratio 1/100. It is found that 2:1 (and simultaneously 4:2) mean motion resonances of the binaries would be strong and could, in principle, draw small third objects deep into relativistic regimes (e.g. ~10 Schwarzschild radii). The inclinations of the trapped objects could increase significantly and, in some cases, retrograde orbits could be realized eventually.

preprint2011arXiv

The Structure of a Self-Gravitating Protoplanetary Disk and its Implications to Direct Imaging Observations

We consider the effects of self-gravity on the hydrostatic balance in the vertical direction of a gaseous disk and discuss the possible signature of the self-gravity that may be captured by the direct imaging observations of protoplanetary disks in future. In this paper, we consider a vertically isothermal disk in order to isolate the effects of self-gravity. The specific disk model we consider in this paper is the one with a radial surface density gap, at which the Toomre's $Q$-parameter of the disk varies rapidly in the radial direction. We calculate the vertical structure of the disk including the effects of self-gravity. We then calculate the scattered light and the dust thermal emission. We find that if the disk is massive enough and the effects of self-gravity come into play, a weak bump-like structure at the gap edge appears in the near-infrared (NIR) scattered light, while no such bump-like structure is seen in the sub-mm dust continuum image. The appearance of the bump is caused by the variation of the height of the surface in the NIR wavelength. If such bump-like feature is detected in future direct imaging observations, with the combination of sub-mm observations, it will bring us useful information about the physical states of the disk.

preprint2010arXiv

Protoplanetary Disk Winds by Magnetorotational Instability : Formation of an Inner Hole and a Crucial Assist for Planet Formation

By constructing a global model based on 3D local magnetohydrodynamical (MHD) simulations, we show that the disk wind driven by magnetorotational instability (MRI) plays a significant role in the dispersal of the gas component of proto-planetary disks. Because the mass loss time scale by the MRI-driven disk winds is proportional to the local Keplerian rotation period, a gas disk dynamically evaporates from the inner region with possibly creating a gradually expanding inner hole, while a sizable amount of the gas remains in the outer region. The disk wind is highly time-dependent with quasi-periodicity of several times Keplerian rotation period at each radius, which will be observed as time-variability of protostar-protoplanetary disk systems. These features persistently hold even if a dead zone exists because the disk winds are driven from the surface regions where ionizing cosmic rays and high energy photons can penetrate. Moreover, the predicted inside-out clearing significantly suppresses the infall of boulders to a central star and the Type I migration of proto-planets which are favorable for the formation and survival of planets.

preprint2010arXiv

Relativistic Astrophysics with Resonant Multiple Inspirals

We show that a massive black hole binary might resonantly trap a small third body (e.g. a neutron star) down to a stage when the binary becomes relativistic due to its orbital decay by gravitational radiation. The final fate of the third body would be quite interesting for relativistic astrophysics. For example, the parent binary could expel the third body with a velocity more than 10% of the speed of light. We also discuss the implications of this three-body system for direct gravitational wave observation.

preprint2010arXiv

Two-Dimensional Study of the Propagation of Planetary Wake and the Indication to Gap Opening in an Inviscid Protoplanetary Disk

We analyze the physical processes of gap formation in an inviscid protoplanetary disk with an embedded protoplanet using two-dimensional local shearing-sheet model. Spiral density wave launched by the planet shocks and the angular momentum carried by the wave is transferred to the background flow. The exchange of the angular momentum can affect the mass flux in the vicinity of the planet to form an underdense region, or gap, around the planetary orbit. We first perform weakly non-linear analyses to show that the specific vorticity formed by shock dissipation of density wave can be a source of mass flux in the vicinity of the planet, and that the gap can be opened even for low-mass planets unless the migration of the planet is substantial. We then perform high resolution numerical simulations to check analytic consideration. By comparing the gap opening timescale and type I migration timescale, we propose a criterion for the formation of underdense region around the planetary orbit that is qualitatively different from previous studies. The minimum mass required for the planet to form a dip is twice as small as previous studies if we incorporate the standard values of type I migration timescale, but it can be much smaller if there is a location in the disk where type I migration is halted.

preprint2009arXiv

Orbital Evolution of a Particle Interacting with a Single Planet in a Protoplanetary Disk

We investigate the motion of a particle around a low mass planet embedded in a non-turbulent gaseous disk. We take into account the effect of the gas structure that is modified by the gravitational interaction between the planet. We derive an analytic formula that describes the change of the semi-major axis of the particle due to the encounter with the planet using local approximation in distant encounter regime. Our final formula includes the effects of steady, axisymmetric radial gas flow, the global gas pressure gradient in the disk, planet gravity, and the structure of the gas flow modified by the planet's gravity. We compare the analytic results with numerical calculations, and indicate that our formula well describes the secular evolution of the dust particles' semi-major axes well, especially for small particles with large drag coefficient. We discuss the conditions for dust gap opening around a low mass planet and radial distribution of dust particles. Our formula may provide a useful tool for calculating radial distribution of particles in a disk around the planet.