Source author record

Takashi Tsukagoshi

Takashi Tsukagoshi 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

34works
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

34 published item(s)

preprint2022arXiv

A new method for direct measurement of isotopologue ratios in protoplanetary disks: a case study of the $^{12}$CO/$^{13}$CO ratio in the TW Hya disk

Planetary systems are thought to be born in protoplanetary disks. Isotope ratios are a powerful tool for investigating the material origin and evolution from molecular clouds to planetary systems via protoplanetary disks. However, it is challenging to measure the isotope (isotopologue) ratios, especially in protoplanetary disks, because the emission lines of major species are saturated. We developed a new method to overcome these challenges by using optically thin line wings induced by thermal broadening. As a first application of the method, we analyzed two carbon monoxide isotopologue lines, $^{12}$CO $3-2$ and $^{13}$CO $3-2$, from archival observations of a protoplanetary disk around TW Hya with the Atacama Large Millimeter/sub-millimeter Array. The $^{12}$CO/$^{13}$CO ratio was estimated to be ${ 20\pm5}$ at disk radii of ${ 70-110}$ au, which is significantly smaller than the value observed in the local interstellar medium, $\sim69$. It implies that an isotope exchange reaction occurs in a low-temperature environment with $\rm C/O>1$ . In contrast, it is suggested that $^{12}$CO/$^{13}$CO is higher than $\sim{ 84}$ in the outer disk ($r > { 130}$ au), which can be explained by the difference in the binding energy of the isotopologues on dust grains and the CO gas depletion processes. Our results imply that the gas-phase $^{12}$CO/$^{13}$CO can vary by a factor of ${ > 4}$ even inside a protoplanetary disk, and therefore, can be used to trace material evolution in disks.

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

CO Line Emission Surfaces and Vertical Structure in Mid-Inclination Protoplanetary Disks

High spatial resolution CO observations of mid-inclination (30-75°) protoplanetary disks offer an opportunity to study the vertical distribution of CO emission and temperature. The asymmetry of line emission relative to the disk major axis allows for a direct mapping of the emission height above the midplane, and for optically-thick, spatially-resolved emission in LTE, the intensity is a measure of the local gas temperature. Our analysis of ALMA archival data yields CO emission surfaces, dynamically-constrained stellar host masses, and disk atmosphere gas temperatures for the disks around: HD 142666, MY Lup, V4046 Sgr, HD 100546, GW Lup, WaOph 6, DoAr 25, Sz 91, CI Tau, and DM Tau. These sources span a wide range in stellar masses (0.50-2.10 M$_{\odot}$), ages (${\sim}$0.3-23 Myr), and CO gas radial emission extents (${\approx}$200-1000 au). This sample nearly triples the number of disks with mapped emission surfaces and confirms the wide diversity in line emitting heights ($z/r\approx0.1$ to ${\gtrsim}0.5$) hinted at in previous studies. We compute radial and vertical CO gas temperature distributions for each disk. A few disks show local temperature dips or enhancements, some of which correspond to dust substructures or the proposed locations of embedded planets. Several emission surfaces also show vertical substructures, which all align with rings and gaps in the millimeter dust. Combining our sample with literature sources, we find that CO line emitting heights weakly decline with stellar mass and gas temperature, which, despite large scatter, is consistent with simple scaling relations. We also observe a correlation between CO emission height and disk size, which is due to the flared structure of disks. Overall, CO emission surfaces trace ${\approx}2$-$5\times$ gas pressure scale heights (H$_{\rm{g}}$) and could potentially be calibrated as empirical tracers of H$_{\rm{g}}$.

preprint2022arXiv

Detection of HC18O+ in a protoplanetary disk: exploring oxygen isotope fractionation of CO

The oxygen isotope fractionation scenario, which has been developed to explain the oxygen isotope anomaly in the solar system materials, predicts that CO gas is depleted in 18O in protoplanetary disks, where segregation between solids and gas inside disks had already occurred. Based on ALMA observations, we report the first detection of HC18O+(4-3) in a Class II protoplanetary disk (TW Hya). This detection allows us to explore the oxygen isotope fractionation of CO in the TW Hya disk from optically thin HCO+ isotopologues as a proxy of optically thicker CO isotopologues. Using the H13CO+(4-3) data previously obtained with SMA, we find that the H13CO+/HC18O+ ratio in the central <100 au regions of the disk is 10.3 +- 3.2. We construct a chemical model of the TW Hya disk with carbon and oxygen isotope fractionation chemistry, and estimate the conversion factor from H13CO+/HC18O+ to 13CO/C18O. With the conversion factor (= 0.8), the 13CO/C18O ratio is estimated to be 8.3 +- 2.6, which is consistent with the elemental abundance ratio in the local ISM (8.1 +- 0.8) within error margin. Then there is no clear evidence of 18O depletion in CO gas of the disk, although we could not draw any robust conclusion due to large uncertainties. In conclusion, optically thin lines of HCO+ isotopologues are useful tracers of CO isotopic ratios, which are hardly constrained directly from optically thick lines of CO isotopologues. Future higher sensitivity observations of H13CO+ and HC18O+ would be able to allow us to better constrain the oxygen fractionation in the disk.

preprint2022arXiv

Massive compact dust disk with a gap around CW Tau revealed by ALMA multi-band observations

Compact protoplanetary disks with a radius of $\lesssim$ 50 au are common around young low-mass stars. We report high resolution ALMA dust continuum observations toward a compact disk around CW Tau at Band 4 ($λ=2.2$ mm), 6 (1.3 mm), 7 (0.89 mm) and 8 (0.75 mm). The SED shows the spectral slope of $2.0\pm0.24$ between 0.75 and 1.3 mm, while it is $3.7\pm0.29$ between 2.17 and 3.56 mm. The steep slope between 2.17 and 3.56 mm is consistent with that of optically thin emission from small grains ($\lesssim$ 350 ${\rm μm}$). We perform parametric fitting of the ALMA data to characterize the dust disk. Interestingly, if the dust-to-gas mass ratio is 0.01, the Toomre's Q parameter reaches $\sim$ 1-3, suggesting that the CW Tau disk might be marginally gravitationally unstable. The total dust mass is estimated as $\sim250M_{\oplus}$ for the maximum dust size of 140 ${\rm μm}$ that is inferred from the previous Band 7 polarimetric observation and at least $80M_{\oplus}$ even for larger grain sizes. This result shows that the CW Tau disk is quite massive in spite of its smallness. Furthermore, we clearly identify a gap structure located at $\sim20$ au, which might be induced by a giant planet. In spite of these interesting characteristics, the CW Tau disk has normal disk luminosity, size and spectral index at ALMA Band 6, which could be a clue to the mass budget problem in Class II disks.

preprint2022arXiv

Molecules with ALMA at Planet-forming Scales (MAPS) III: Characteristics of Radial Chemical Substructures

The Molecules with ALMA at Planet-forming Scales (MAPS) Large Program provides a detailed, high resolution (${\sim}$10-20 au) view of molecular line emission in five protoplanetary disks at spatial scales relevant for planet formation. Here, we present a systematic analysis of chemical substructures in 18 molecular lines toward the MAPS sources: IM Lup, GM Aur, AS 209, HD 163296, and MWC 480. We identify more than 200 chemical substructures, which are found at nearly all radii where line emission is detected. A wide diversity of radial morphologies - including rings, gaps, and plateaus - is observed both within each disk and across the MAPS sample. This diversity in line emission profiles is also present in the innermost 50 au. Overall, this suggests that planets form in varied chemical environments both across disks and at different radii within the same disk. Interior to 150 au, the majority of chemical substructures across the MAPS disks are spatially coincident with substructures in the millimeter continuum, indicative of physical and chemical links between the disk midplane and warm, elevated molecular emission layers. Some chemical substructures in the inner disk and most chemical substructures exterior to 150 au cannot be directly linked to dust substructure, however, which indicates that there are also other causes of chemical substructures, such as snowlines, gradients in UV photon fluxes, ionization, and radially-varying elemental ratios. This implies that chemical substructures could be developed into powerful probes of different disk characteristics, in addition to influencing the environments within which planets assemble. This paper is part of the MAPS special issue of the Astrophysical Journal Supplement.

preprint2022arXiv

Unveiling the outer dust disc of TW Hya with deep ALMA observations

The radial extent of millimetre dust in protoplanetary discs is often far smaller than that of their gas, mostly due to processes such as dust growth and radial drift. However, it has been suggested that current millimetre continuum observations of discs do not trace their full extent due to limited sensitivity. In this Letter, we present deep (19 $μ$Jy beam$^{-1}$) moderate resolution (0.37") ALMA observations at 1 mm of the nearest protoplanetary disc, TW Hya. Using the visibility analysis tool `frank', we reveal a structured millimetre intensity distribution out to 100 au, well beyond previous estimates of 60-70 au. Our analysis suggests the presence of a new millimetre continuum gap at 82 au, coincident with similar features seen in optical/near-infrared scattered light and millimetre molecular line observations. Examination of the fit residuals confirms the presence of the previously reported au-scale continuum excess at 52 au (P.A. = 242.5 degrees). Our results demonstrate the utility of combining deep, moderate resolution observations with super-resolution analysis techniques to probe the faintest regions of protoplanetary discs.

preprint2020arXiv

$^{14}$N/$^{15}$N isotopic ratio in CH$_3$CN of Titan's atmosphere measured with ALMA

Each of the nitriles present in the atmosphere of Titan can be expected to exhibit different \ce{^{14}N/^{15}N} values depending on their production processes, primarily because of the various \ce{N2} dissociation processes induced by different sources such as ultraviolet radiation, magnetospheric electrons, and galactic cosmic rays. For \ce{CH3CN}, one photochemical model predicted a \ce{^{14}N/^{15}N} value as 120--130 in the lower stratosphere. This is much higher than that for \ce{HCN} and \ce{HC3N}, $\sim$67--94. By analyzing archival data obtained by the Atacama Large Millimeter/submillimeter Array (ALMA), we successfully detected submillimeter rotational transitions of \ce{CH3C^{15}N} ($J$ = 19--18) locate at the 338 GHz band in Titan's atmospheric spectra. By comparing those observations with the simultaneously observed \ce{CH3CN} ($J$ = 19--18) lines at the 349 GHz band, which probe from 160 to $\sim$400 km altitude, we then derived \ce{^{14}N/^{15}N} in \ce{CH3CN} as 125$^{+145}_{-44}$. Although the range of the derived value shows insufficient accuracy due to data quality limitations, the best-fit value suggests that \ce{^{14}N/^{15}N} for \ce{CH3CN} is higher than values that have been previously observed and theoretically predicted for \ce{HCN} and \ce{HC3N}. This may be explained by the different \ce{N2} dissociation sources according to the altitudes, as suggested by a recent photochemical model.

preprint2020arXiv

DESHIMA on ASTE: On-sky Responsivity Calibration of the Integrated Superconducting Spectrometer

We are developing an ultra-wideband spectroscopic instrument, DESHIMA (DEep Spectroscopic HIgh-redshift MApper), based on the technologies of an on-chip filter-bank and Microwave Kinetic Inductance Detector (MKID) to investigate dusty star-burst galaxies in the distant universe at millimeter and submillimeter wavelength. An on-site experiment of DESHIMA was performed using the ASTE 10-m telescope. We established a responsivity model that converts frequency responses of the MKIDs to line-of-sight brightness temperature. We estimated two parameters of the responsivity model using a set of skydip data taken under various precipitable water vapor (PWV, 0.4-3.0 mm) conditions for each MKID. The line-of-sight brightness temperature of sky is estimated using an atmospheric transmission model and the PWVs. As a result, we obtain an average temperature calibration uncertainty of $1σ=4$%, which is smaller than other photometric biases. In addition, the average forward efficiency of 0.88 in our responsivity model is consistent with the value expected from the geometrical support structure of the telescope. We also estimate line-of-sight PWVs of each skydip observation using the frequency response of MKIDs, and confirm the consistency with PWVs reported by the Atacama Large Millimeter/submillimeter Array.

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

Scattering-induced intensity reduction: large mass content with small grains in the inner region of the TW Hya disk

Dust continuum observation is one of the best methods to constrain the properties of protoplanetary disks. Recent theoretical studies have suggested that the dust scattering at the millimeter wavelength potentially reduces the observed intensity, which results in an underestimate in the dust mass. We investigate whether the dust scattering indeed reduces the observed continuum intensity by comparing the ALMA archival data of the TW Hya disk at Band 3, 4, 6, 7 and 9 to models obtained by radiative transfer simulations. We find that the model with scattering by 300 ${\rm μm}$-sized grains well reproduces the observed SED of the central part of the TW Hya disk while the model without scattering is also consistent within the errors of the absolute fluxes. To explain the intensity at Band 3, the dust surface density needs to be $\sim$ 10 ${\rm g\,cm^{-2}}$ at 10 au in the model with scattering, which is 26 times more massive than previously predicted. The model without scattering needs 2.3 times higher dust mass than the model with scattering because it needs lower temperature. At Band 7, scattering reduces the intensity by $\sim$ 35% which makes the disk looks optically thin even though it is optically thick. Our study suggests the TW Hya disk is still capable of forming cores of giant planets at where the current solar system planets exist.

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.

preprint2020arXiv

TW Hya: an old protoplanetary disc revived by its planet

Dark rings with bright rims are the indirect signposts of planets embedded in protoplanetary discs. In a recent first, an azimuthally elongated AU-scale blob, possibly a planet, was resolved with ALMA in TW Hya. The blob is at the edge of a cliff-like rollover in the dust disc rather than inside a dark ring. Here we build time-dependent models of TW Hya disc. We find that the classical paradigm cannot account for the morphology of the disc and the blob. We propose that ALMA-discovered blob hides a Neptune mass planet losing gas and dust. We show that radial drift of mm-sized dust particles naturally explains why the blob is located on the edge of the dust disc. Dust particles leaving the planet perform a characteristic U-turn relative to it, producing an azimuthally elongated blob-like emission feature. This scenario also explains why a 10 Myr old disc is so bright in dust continuum. Two scenarios for the dust-losing planet are presented. In the first, a dusty pre-runaway gas envelope of about 40 Earth mass Core Accretion planet is disrupted, e.g., as a result of a catastrophic encounter. In the second, a massive dusty pre-collapse gas giant planet formed by Gravitational Instability is disrupted by the energy released in its massive core. Future modelling may discriminate between these scenarios and allow us to study planet formation in an entirely new way -- by analysing the flows of dust and gas recently belonging to planets, informing us about the structure of pre-disruption planetary envelopes.

preprint2019arXiv

Nobeyama 45-m mapping observations toward nearby molecular clouds, Orion A, Aquila Rift, and M17: Project overview

We carried out mapping observations toward three nearby molecular clouds, Orion A, Aquila Rift, and M17, using a new 100 GHz receiver, FOREST, on the Nobeyama 45-m telescope. In the present paper, we describe the details of the data obtained such as intensity calibration, data sensitivity, angular resolution, and velocity resolution. Each target contains at least one high-mass star-forming region. The target molecular lines were $^{12}$CO ($J = 1 - 0$), $^{13}$CO ($J = 1 - 0$), C$^{18}$O ($J = 1 - 0$), N$_2$H$^+$ ($J=1-0$), and CCS ($J_N=8_7-7_6$), with which we covered the density range of 10$^2$ cm$^{-3}$ to 10$^6$ cm$^{-3}$ with an angular resolution of $\sim 20\arcsec$ and a velocity resolution of $\sim$ 0.1 km s$^{-1}$. Assuming the representative distances of 414 pc, 436 pc, and 2.1 kpc, the maps of Orion A, Aquila Rift, and M17 cover most of the densest parts with areas of about 7 pc $\times$ 15 pc, 7 pc $\times$ 7 pc, and 36 pc $\times$ 18 pc, respectively. On the basis of the $^{13}$CO column density distribution, the total molecular masses are derived to be $3.86 \times 10^4 M_\odot$, $2.67 \times 10^4 M_\odot$, and $8.1\times 10^5 M_\odot$ for Orion A, Aquila Rift, and M17, respectively. For all the clouds, the H$_2$ column density exceeds the theoretical threshold for high-mass star formation of $\gtrsim$ 1 g cm$^{-2}$, only toward the regions which contain current high-mass star-forming sites. For other areas, further mass accretion or dynamical compression would be necessary for future high-mass star formation. This is consistent with the current star formation activity. Using the $^{12}$CO data, we demonstrate that our data have enough capability to identify molecular outflows, and for Aquila Rift, we identify 4 new outflow candidates. The scientific results will be discussed in details in separate papers.

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.

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

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

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

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 intrinsic abundance ratio and X-factor of CO isotopologues in L1551 shielded from FUV photodissociation

L1551 is chosen because it is relatively isolated in the Taurus molecular cloud shielded from FUV photons, providing an ideal environment for studying the target properties. Our observations cover ~40'x40' with resolution ~30", which are the maps with highest spatial dynamical range to date. We derive the X(13CO)/X(C18O) value on the sub-parsec scales in the range of ~3-27 with a mean value of 8.0+-2.8. Comparing to the visual extinction map derived from the Herschel observations, we found that the abundance ratio reaches its maximum at low Av (i.e., Av ~ 1-4mag), and decreases to the typical solar system value of 5.5 inside L1551 MC. The high X(13CO)/X(C18O) value at the boundary of the cloud is most likely due to the selective FUV photodissociation of C18O. This is in contrast with Orion-A where its internal OB stars keep the abundance ratio at a high level greater than ~10. In addition, we explore the variation of the X-factor, because it is an uncertain but widely used quantity in extragalactic studies. We found that X-factor is proportional to N(H2)^1.0 which is consistent with previous simulations. Excluding the high density region, the average X-factor is similar to the Milky Way average value.

preprint2015arXiv

First Detection of [C I] $^3$P$_1$-$^3$P$_0$ Emission from a Protoplanetary Disk

We performed single point [C I] $^3$P$_1$-$^3$P$_0$ and CO J=4-3 observations toward three T Tauri stars, DM Tau, LkCa 15, and TW Hya, using the Atacama Large Millimeter/submillimeter Array (ALMA) Band 8 qualification model receiver installed on the Atacama Submillimeter Telescope Experiment (ASTE). Two protostars in the Taurus L1551 region, L1551 IRS 5 and HL Tau, were also observed. We successfully detected [C I] emission from the protoplanetary disk around DM Tau as well as the protostellar targets. The spectral profile of the [C I] emission from the protoplanetary disk is marginally single-peaked, suggesting that atomic carbon (C) extends toward the outermost disk. The detected [C I] emission is optically thin and the column densities of C are estimated to be <~10$^{16}$ cm$^{-2}$ and ~10$^{17}$ cm$^{-2}$ for the T Tauri star targets and the protostars, respectively. We found a clear difference in the total mass ratio of C to dust, $M$(C)/$M$(dust), between the T Tauri stars and protostellar targets; the $M$(C)/$M$(dust) ratio of the T Tauri stars is one order of magnitude smaller than that of the protostars. The decrease of the estimated $M$(C)/$M$(dust) ratios for the disk sources is consistent with a theoretical prediction that the atomic C can survive only in the near surface layer of the disk and C$^+$/C/CO transition occurs deeper into the disk midplane.

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

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

Dense Molecular Clumps associated with the LMC Supergiant Shells LMC 4 \& LMC 5

We investigate the effects of Supergiant Shells (SGSs) and their interaction on dense molecular clumps by observing the Large Magellanic Cloud (LMC) star forming regions N48 and N49, which are located between two SGSs, LMC 4 and LMC 5. $^{12}$CO ($J$=3-2, 1-0) and $^{13}$CO ($J$=1-0) observations with the ASTE and Mopra telescopes have been carried out towards these regions. A clumpy distribution of dense molecular clumps is revealed with 7 pc spatial resolution. Large velocity gradient analysis shows that the molecular hydrogen densities ($n({\rm H}_2)$) of the clumps are distributed from low to high density ($10^3$-$10^5$ cm$^{-3}$) and their kinetic temperatures ($T_{\rm kin}$) are typically high (greater than $50$ K). These clumps seem to be in the early stages of star formation, as also indicated from the distribution of H$α$, young stellar object candidates, and IR emission. We found that the N48 region is located in the high column density HI envelope at the interface of the two SGSs and the star formation is relatively evolved, whereas the N49 region is associated with LMC 5 alone and the star formation is quiet. The clumps in the N48 region typically show high $n({\rm H}_2)$ and $T_{\rm kin}$, which are as dense and warm as the clumps in LMC massive cluster-forming areas (30 Dor, N159). These results suggest that the large-scale structure of the SGSs, especially the interaction of two SGSs, works efficiently on the formation of dense molecular clumps and stars.

preprint2014arXiv

High abundance ratio of $^{13}$CO to C$^{18}$O toward photon-dominated regions in the Orion-A giant molecular cloud

Aims. We derive physical properties such as the optical depths and the column densities of $^{13}$CO and C$^{18}$O to investigate the relationship between the far ultraviolet (FUV) radiation and the abundance ratios between $^{13}$CO and C$^{18}$O. Method. We have carried out wide-field (0.4 deg$^2$) observations with an angular resolution of 25.8 arcsec ($\sim$ 0.05 pc) in $^{13}$CO ($J$=1--0) and C$^{18}$O ($J$=1--0) toward the Orion-A giant molecular cloud using the Nobeyama 45 m telescope in the on-the-fly mode. Results. Overall distributions and velocity structures of the $^{13}$CO and C$^{18}$O emissions are similar to those of the $^{12}$CO ($J$=1--0) emission. The optical depths of the $^{13}$CO and C18O emission lines are estimated to be 0.05 $<$ $τ_{\rm ^{13}CO}$ $<$ 1.54 and 0.01 $<$ $τ_{\rm C^{18}O}$ $<$ 0.18, respectively. The column densities of the $^{13}$CO and C$^{18}$O emission lines are estimated to be 0.2 $\times$ 10$^{16}$ $<$ $N_{\rm ^{13}CO}$ $<$ 3.7 $\times$ 10$^{17}$ cm$^{-2}$ and 0.4 $\times$ 10$^{15}$ $<$ $N_{\rm C^{18}O}$ $<$ 3.5 $\times$ 10$^{16}$ cm$^{-2}$, respectively. The abundance ratios between $^{13}$CO and C$^{18}$O, $X_{\rm ^{13}CO}$/$X_{\rm C^{18}O}$, are found to be 5.7 - 33.0. The mean value of $X_{\rm ^{13}CO}$/$X_{\rm C^{18}O}$ in the nearly edge-on photon-dominated regions is found to be 16.47 $\pm$ 0.10, which is a third larger than that the solar system value of 5.5. The mean value of $X_{\rm ^{13}CO}$/$X_{\rm C^{18}O}$ in the other regions is found to be 12.29 $\pm$ 0.02. The difference of the abundance ratio is most likely due to the selective FUV photodissociation of C$^{18}$O.

preprint2014arXiv

High-Resolution Submillimeter and Near-Infrared Studies of the Transition Disk around Sz 91

To reveal the structures of a transition disk around a young stellar object in Lupus, Sz 91, we have performed aperture synthesis 345 GHz continuum and CO(3--2) observations with the Submillimeter Array ($\sim1\arcsec$--3$\arcsec$ resolution), and high-resolution imaging of polarized intensity at the $K_s$-band by using the HiCIAO instrument on the Subaru Telescope ($0\farcs25$ resolution). Our observations successfully resolved the inner and outer radii of the dust disk to be 65 AU and 170 AU, respectively, which indicates that Sz 91 is a transition disk source with one of the largest known inner holes. The model fitting analysis of the spectral energy distribution reveals an H$_2$ mass of $2.4\times10^{-3}$ $M_\sun$ in the cold ($T<$30 K) outer part at $65<r<170$ AU by assuming a canonical gas-to-dust mass ratio of 100, although a small amount ($>3\times10^{-9}$ $M_\sun$) of hot ($T\sim$180 K) dust possibly remains inside the inner hole of the disk. The structure of the hot component could be interpreted as either an unresolved self-luminous companion body (not directly detected in our observations) or a narrow ring inside the inner hole. Significant CO(3--2) emission with a velocity gradient along the major axis of the dust disk is concentrated on the Sz 91 position, suggesting a rotating gas disk with a radius of 420 AU. The Sz 91 disk is possibly a rare disk in an evolutionary stage immediately after the formation of protoplanets because of the large inner hole and the lower disk mass than other transition disks studied thus far.

preprint2013arXiv

Extensive [CI] Mapping toward the Orion-A Molecular Cloud

We have carried out wide-field (0.17 degree^2) and high-angular resolution (21.3 arcsec ~ 0.04 pc) observations in [CI] line toward the Orion-A giant molecular cloud with the Atacama Submillimeter Telescope Experiment (ASTE) 10 m telescope in the On-The-Fly (OTF) mode. Overall features of the [CI] emission are similar to those of the CO (1--0) emission in Shimajiri et al. (2011); the total intensity ratio of the [CI] to CO emission ranges from 0.05 to 0.2. The optical depth of the [CI] emission is found to be 0.1 -- 0.75, suggesting optically thin emission. The column density of the [CI] emission is estimated to be (1.0 -- 19) x 10^17 cm^-2. These results are consistent with the results of the previous [CI] observations with a low-angular resolution of 2.2 arcmin (e.g. Ikeda et al. 1999). In the nearly edge-on PDRs and their candidates of the Orion Bar, DLSF, M 43 Shell, and Region D, the distributions of the [CI] emission coincide with those of the CO emission, inconsistent with the prediction by the plane-parallel PDR model (Hollenbach & Tielens 1999). In addition, the [CI] distribution in the Orion A cloud is found to be more similar to those of the ^{13}CO (1--0), C^{18}O (1--0), and H^{13}CO^+ (1--0) lines than that of the CO (1--0) line, suggesting that the [CI] emission is not limited to the cloud surface, but is tracing the dense, inner parts of the cloud.

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.

preprint2011arXiv

Evidence For Cloud-Cloud Collision and Parsec-Scale Stellar Feedback Within the L1641-N Region

We present high spatial resolution $^{12}$CO ($J=1-0$) images taken by the Nobeyama 45m telescope toward a $48' \times 48'$ area including the L1641-N cluster. The effective spatial resolution of the maps is $21"$, corresponding to 0.04 pc at a distance of 400 pc. A recent 1.1 mm dust continuum map reveals that the dense gas is concentrated in several thin filaments. We find that a few dust filaments are located at the parts where $^{12}$CO ($J=1-0$) emission drops sharply. Furthermore, the filaments have two-components with different velocities. The velocity difference between the two-components is about 3 km s$^{-1}$, corresponding to a Mach number of 10, significantly larger than the local turbulent velocity in the cloud. These facts imply that the collision of the two components (hereafter, the cloud-cloud collision) possibly contributed to the formation of these filaments. Since the two components appear to overlap toward the filaments on the plane of the sky, the collision may have occurred almost along the line of sight. Star formation in the L1641-N cluster was probably triggered by such a collision. We also find several parsec-scale CO shells whose centers are close to either the L1641-N cluster or V 380 Ori cluster. We propose that these shells were created by multiple winds and/or outflows from cluster YSOs, i.e., "protocluster winds." One exceptional dust filament located at the western cloud edge lies along a shell; it is presumably a part of the expanding shell. Both the cloud-cloud collision and protocluster winds are likely to influence the cloud structure and kinematics in this region.

preprint2011arXiv

Molecular Outflows From the Protocluster, Serpens South

We present the results of CO ($J=3-2$) and HCO$^+$ ($J=4-3$) mapping observations toward a nearby embedded cluster, Serpens South, using the ASTE 10 m telescope. Our CO ($J=3-2$) map reveals that many outflows are crowded in the dense cluster-forming clump that can be recognized as a HCO$^+$ clump with a size of $\sim$ 0.2 pc and mass of $\sim$ 80 M$_\odot$. The clump contains several subfragments with sizes of $\sim$ 0.05 pc. By comparing the CO ($J=3-2$) map with the 1.1 mm dust continuum image taken by AzTEC on ASTE, we find that the spatial extents of the outflow lobes are sometimes anti-correlated with the distribution of the dense gas and some of the outflow lobes apparently collide with the dense gas. The total outflow mass, momentum, and energy are estimated at 0.6 $M_\odot$, 8 $M_\odot$ km s$^{-1}$, and 64 $M_\odot$ km$^2$ s$^{-2}$, respectively. The energy injection rate due to the outflows is comparable to the turbulence dissipation rate in the clump, implying that the protostellar outflows can maintain the supersonic turbulence in this region. The total outflow energy seems only about 10 percent the clump gravitational energy. We conclude that the current outflow activity is not enough to destroy the whole cluster-forming clump, and therefore star formation is likely to continue for several or many local dynamical times.

preprint2010arXiv

Detection of Strong Millimeter Emission from the Circumstellar Dust Disk Around V1094 Sco: Cold and Massive Disk around a T Tauri Star in a Quiescent Accretion Phase?

We present the discovery of a cold massive dust disk around the T Tauri star V1094 Sco in the Lupus molecular cloud from the 1.1 millimeter continuum observations with AzTEC on ASTE. A compact ($r\lesssim$320 AU) continuum emission coincides with the stellar position having a flux density of 272 mJy which is largest among T Tauri stars in Lupus. We also present the detection of molecular gas associated with the star in the five-point observations in $^{12}$CO J=3--2 and $^{13}$CO J=3--2. Since our $^{12}$CO and $^{13}$CO observations did not show any signature of a large-scale outflow or a massive envelope, the compact dust emission is likely to come from a disk around the star. The observed SED of V1094 Sco shows no distinct turnover from near infrared to millimeter wavelengths, which can be well described by a flattened disk for the dust component, and no clear dip feature around 10 $\micron$ suggestive of absence of an inner hole in the disk. We fit a simple power-law disk model to the observed SED. The estimated disk mass ranges from 0.03 to $\gtrsim$0.12 $M_\sun$, which is one or two orders of magnitude larger than the median disk mass of T Tauri stars in Taurus.

preprint2010arXiv

The Molecular Outflows in the rho Ophiuchi Main Cloud: Implications For Turbulence Generation

We present the results of CO (J=3-2) and CO (J=1-0) mapping observations toward the active cluster forming clump, L1688, in the rho Ophiuchi molecular cloud. From the CO (J=3-2) and CO (J=1-0) data cubes, we identify five outflows, whose driving sources are VLA 1623, EL 32, LFAM 26, EL 29, and IRS 44. Among the identified outflows, the most luminous outflow is the one from the prototypical Class 0 source, VLA 1623. We also discover that the EL 32 outflow located in the Oph B2 region has very extended blueshifted and redshifted lobes with wide opening angles. This outflow is most massive and have the largest momentum among the identified outflows in the CO (J=1-0) map. We estimate the total energy injection rate due to the molecular outflows identified by the present and previous studies to be about 0.2 L_solar, larger than or at least comparable to the turbulence dissipation rate [~(0.03 - 0.1) L_solar]. Therefore, we conclude that the protostellar outflows are likely to play a significant role in replenishing the supersonic turbulence in this clump.