Showing 1–39 of 39 results for author: Tazaki, R

Support for fragile porous dust in a gravitationally self-regulated disk around IM Lup

Authors: Takahiro Ueda, Ryo Tazaki, Satoshi Okuzumi, Mario Flock, Prakruti Sudarshan

Abstract: Protoplanetary disks, the birthplace of planets, are expected to be gravitationally unstable in their early phase of evolution. IM Lup, a well-known T-Tauri star, is surrounded by a protoplanetary disk with spiral arms likely caused by gravitational instability. The IM Lup disk has been observed using various methods, but developing a unified explanatory model is challenging. Here we present a phy… ▽ More Protoplanetary disks, the birthplace of planets, are expected to be gravitationally unstable in their early phase of evolution. IM Lup, a well-known T-Tauri star, is surrounded by a protoplanetary disk with spiral arms likely caused by gravitational instability. The IM Lup disk has been observed using various methods, but developing a unified explanatory model is challenging. Here we present a physical model of the IM Lup disk that offers a comprehensive explanation for diverse observations spanning from near-infrared to millimeter wavelengths. Our findings underscore the importance of dust fragility in retaining the observed millimeter emission and reveal the preference for moderately porous dust to explain observed millimeter polarization. We also find that the inner disk region is likely heated by gas accretion, providing a natural explanation for bright millimeter emission within 20 au. The actively heated inner region in the model casts a 100-au-scale shadow, aligning seamlessly with the near-infrared scattered light observation. The presence of accretion heating also supports the fragile dust scenario in which accretion efficiently heat the disk midplane. Due to the fragility of dust, it is unlikely that a potential embedded planet at 100 au formed via pebble accretion in a smooth disk, pointing to local dust enhancement boosting pebble accretion or alternative pathways such as outward migration or gravitational fragmentation. △ Less

Submitted 11 June, 2024; originally announced June 2024.

Comments: Accepted for publication. Author version

The SPHERE view of the Chamaeleon I star-forming region

Authors: C. Ginski, A. Garufi, M. Benisty, R. Tazaki, C. Dominik, A. Ribas, N. Engler, T. Birnstiel, G. Chauvin, G. Columba, S. Facchini, A. Goncharov, J. Hagelberg, T. Henning, M. Hogerheijde, R. G. van Holstein, J. Huang, T. Muto, P. Pinilla, K. Kanagawa, S. Kim, N. Kurtovic, M. Langlois, C. Manara, J. Milli , et al. (10 additional authors not shown)

Abstract: We used VLT/SPHERE to observe 20 systems in the Cha I cloud in polarized scattered light in the near-infrared. We combined the scattered light observations with existing literature data on stellar properties and with archival ALMA continuum data to study trends with system age and dust mass. We also connected resolved near-infrared observations with the spectral energy distributions of the systems… ▽ More We used VLT/SPHERE to observe 20 systems in the Cha I cloud in polarized scattered light in the near-infrared. We combined the scattered light observations with existing literature data on stellar properties and with archival ALMA continuum data to study trends with system age and dust mass. We also connected resolved near-infrared observations with the spectral energy distributions of the systems. In 13 of the 20 systems included in this study we detected resolved scattered light signals from circumstellar dust. For the CR Cha, CT Cha, CV Cha, SY Cha, SZ Cha, and VZ Cha systems we present the first detailed descriptions of the disks in scattered light. The observations found typically smooth or faint disks, often with little substructure, with the notable exceptions of SZ Cha, which shows an extended multiple-ringed disk, and WW Cha, which shows interaction with the cloud environment. New high S/N K- band observations of the HD 97048 system in our survey reveal a significant brightness asymmetry that may point to disk misalignment and subsequent shadowing of outer disk regions, possibly related to the suggested planet candidate in the disk. We resolve for the first time the stellar binary in the CS Cha system. Multiple wavelength observations of the disk around CS Cha have revealed that the system contains small, compact dust grains that may be strongly settled, consistent with numerical studies of circumbinary disks. We find in our sample that there is a strong anti-correlation between the presence of a (close) stellar companion and the detection of circumstellar material with five of our seven nondetections located in binary systems. △ Less

Submitted 4 March, 2024; originally announced March 2024.

Comments: Accepted by A&A

Insight from laboratory measurements on dust in debris discs

Authors: Julien Milli, Olivier Poch, Jean-Baptiste Renard, Jean-Charles Augereau, Pierre Beck, Elodie Choquet, Jean-Michel Geffrin, Edith Hadamcik, Jérémie Lasue, François Ménard, Arthur Péronne, Clément Baruteau, Ryo Tazaki, Vanesa Tobon Valencia

Abstract: Extreme adaptive optics instruments have revealed exquisite details on debris discs, allowing to extract the optical properties of the dust particles such as the phase function, the degree of polarisation and the spectral reflectance. These are three powerful diagnostic tools to understand the physical properties of the dust : the size, shape and composition of the dust particles. This can inform… ▽ More Extreme adaptive optics instruments have revealed exquisite details on debris discs, allowing to extract the optical properties of the dust particles such as the phase function, the degree of polarisation and the spectral reflectance. These are three powerful diagnostic tools to understand the physical properties of the dust : the size, shape and composition of the dust particles. This can inform us on the population of parent bodies, also called planetesimals, which generate those particles through collisions. It is however very rare to be able to combine all those three observables for the same system, as this requires different high-contrast imaging techniques to suppress the starlight and reveal the faint scattered light emission from the dust. Due to its brightness, the ring detected around the A-type star HR 4796 is a notable exception, with both unpolarised and polarised images covering near-infrared wavelengths. Here, we show how measurements of dust particles in the laboratory can reproduce the observed near-infrared photo-polarimetric properties of the HR 4796 disc. Experimental characterisation of dust allows to bypass the current limitations of dust models to reproduce simultaneously the phase function, the degree of polarisation and the spectral reflectance. △ Less

Submitted 4 December, 2023; originally announced December 2023.

Comments: 4 pages, 4 figures, Proceedings of the Annual meeting of the French Society of Astronomy and Astrophysics

Journal ref: sf2a.conf, 275-278 (2023)

The polarisation properties of the HD 181327 debris ring. Evidence for sub-micron particles from scattered light observations

Authors: Julien Milli, Elodie Choquet, Ryo Tazaki, François Ménard, Jean-Charles Augereau, Johan Olofsson, Philippe Thébault, Olivier Poch, Anny-Chantal Levasseur-Regourd, Jérémie Lasue, Jean-Baptiste Renard, Edith Hadamcik, Clément Baruteau, Hans Martin Schmid, Natalia Engler, Rob G. van Holstein, Evgenij Zubko, Anne-Marie Lagrange, Sebastian Marino, Chirstophe Pinte, Carsten Dominik, Anthony Boccaletti, Maud Langlois, Alice Zurlo, Célia Desgrange , et al. (4 additional authors not shown)

Abstract: Polarisation is a powerful remote-sensing tool to study the nature of particles scattering the starlight. It is widely used to characterise interplanetary dust particles in the Solar System and increasingly employed to investigate extrasolar dust in debris discs' systems. We aim to measure the scattering properties of the dust from the debris ring around HD 181327 at near-infrared wavelengths. We… ▽ More Polarisation is a powerful remote-sensing tool to study the nature of particles scattering the starlight. It is widely used to characterise interplanetary dust particles in the Solar System and increasingly employed to investigate extrasolar dust in debris discs' systems. We aim to measure the scattering properties of the dust from the debris ring around HD 181327 at near-infrared wavelengths. We obtained high-contrast polarimetric images of HD 181327 in the H band with the SPHERE / IRDIS instrument on the Very Large Telescope (ESO). We complemented them with archival data from HST / NICMOS in the F110W filter reprocessed in the context of the Archival Legacy Investigations of Circumstellar Environments (ALICE) project. We developed a combined forward-modelling framework to simultaneously retrieve the scattering phase function in polarisation and intensity. We detected the debris disc around HD 181327 in polarised light and total intensity. We measured the scattering phase function and the degree of linear polarisation of the dust at 1.6 micron in the birth ring. The maximum polarisation is 23.6% +/- 2.6% and occurs between a scattering angle of 70 deg and 82 deg. We show that compact spherical particles made of a highly refractive and relatively absorbing material in a differential power-law size distribution of exponent $-3.5$ can simultaneously reproduce the polarimetric and total intensity scattering properties of the dust. This type of material cannot be obtained with a mixture of silicates, amorphous carbon, water ice, and porosity, and requires a more refracting component such as iron-bearing minerals. We reveal a striking analogy between the near-infrared polarisation of comets and that of HD 181327. The methodology developed here combining VLT/SPHERE and HST/NICMOS may be applicable in the future to combine the polarimetric capabilities of SPHERE with the sensitivity of JWST. △ Less

Submitted 4 December, 2023; originally announced December 2023.

Comments: 21 pages, 22 figures, accepted in Astronomy & Astrophysics

Aligned Grains and Scattered Light Found in Gaps of Planet-Forming Disk

Authors: Ian W. Stephens, Zhe-Yu Daniel Lin, Manuel Fernandez-Lopez, Zhi-Yun Li, Leslie W. Looney, Haifeng Yang, Rachel Harrison, Akimasa Kataoka, Carlos Carrasco-Gonzalez, Satoshi Okuzumi, Ryo Tazaki

Abstract: Polarized (sub)millimeter emission from dust grains in circumstellar disks was initially thought to be due to grains aligned with the magnetic field. However, higher resolution multi-wavelength observations along with improved models found that this polarization is dominated by self-scattering at shorter wavelengths (e.g., 870 $μ$m) and by grains aligned with something other than magnetic fields a… ▽ More Polarized (sub)millimeter emission from dust grains in circumstellar disks was initially thought to be due to grains aligned with the magnetic field. However, higher resolution multi-wavelength observations along with improved models found that this polarization is dominated by self-scattering at shorter wavelengths (e.g., 870 $μ$m) and by grains aligned with something other than magnetic fields at longer wavelengths (e.g., 3 mm). Nevertheless, the polarization signal is expected to depend on the underlying substructure, and observations hitherto have been unable to resolve polarization in multiple rings and gaps. HL Tau, a protoplanetary disk located 147.3 $\pm$ 0.5 pc away, is the brightest Class I or Class II disk at millimeter/submillimeter wavelengths. Here we show deep, high-resolution 870 $μ$m polarization observations of HL Tau, resolving polarization in both the rings and gaps. We find that the gaps have polarization angles with a significant azimuthal component and a higher polarization fraction than the rings. Our models show that the disk polarization is due to both scattering and emission from aligned effectively prolate grains. The intrinsic polarization of aligned dust grains is likely over 10%, which is much higher than what was expected in low resolution observations (~1%). Asymmetries and dust features are revealed in the polarization observations that are not seen in non-polarimetric observations. △ Less

Submitted 14 November, 2023; originally announced November 2023.

Comments: Published in Nature

JWST imaging of edge-on protoplanetary disks II. Appearance of edge-on disks with a tilted inner region: case study of IRAS04302+2247

Authors: M. Villenave, K. R. Stapelfeldt, G. Duchene, F. Menard, S. G. Wolff, M. D. Perrin, C. Pinte, R. Tazaki, D. Padgett

Abstract: We present JWST imaging from 2$μ$m to 21$μ$m of the edge-on protoplanetary disk around the embedded young star IRAS04302+2247. The structure of the source shows two reflection nebulae separated by a dark lane. The source extent is dominated by the extended filamentary envelope at $\sim$4.4$μ$m and shorter wavelengths, transitioning at 7$μ$m and longer wavelengths to more compact lobes of scattered… ▽ More We present JWST imaging from 2$μ$m to 21$μ$m of the edge-on protoplanetary disk around the embedded young star IRAS04302+2247. The structure of the source shows two reflection nebulae separated by a dark lane. The source extent is dominated by the extended filamentary envelope at $\sim$4.4$μ$m and shorter wavelengths, transitioning at 7$μ$m and longer wavelengths to more compact lobes of scattered light from the disk itself. The dark lane thickness does not vary significantly with wavelength, which we interpret as an indication for intermediate-sized ($\sim10μ$m) grains in the upper layers of the disk. Intriguingly, we find that the brightest nebula of IRAS40302 switches side between 12.8$μ$m and 21$μ$m. We explore the effect of a tilted inner region on the general appearance of edge-on disks. We find that radiative transfer models of a disk including a tilted inner region can reproduce an inversion in the brightest nebula. In addition, for specific orientations, the model also predicts strong lateral asymmetries, which can occur for more than half possible viewing azimuths. A large number of edge-on protoplanetary disks observed in scattered light show such lateral asymmetries (15/20), which suggests that a large fraction of protoplanetary disks might host a tilted inner region. Stellar spots may also induce lateral asymmetries, which are expected to vary over a significantly shorter timescale. Variability studies of edge-on disks would allow to test the dominant scenario for the origin of these asymmetries. △ Less

Submitted 13 November, 2023; originally announced November 2023.

Comments: Accepted for publication in ApJ

Protoplanetary disks in $K_s$-band total intensity and polarized light

Authors: Bin B. Ren, Myriam Benisty, Christian Ginski, Ryo Tazaki, Nicole L. Wallack, Julien Milli, Antonio Garufi, Jaehan Bae, Stefano Facchini, François Ménard, Paola Pinilla, C. Swastik, Richard Teague, Zahed Wahhaj

Abstract: Diverse protoplanetary disk morphology can result from planet-disk interaction, suggesting planetary presence. To date, most scattered light imaging campaigns have probed polarized light, which is only a fraction of the total light and not very sensitive to planets. To observe and characterize protoplanetary disk systems in the near-infrared in both polarized and total intensity light, we carried… ▽ More Diverse protoplanetary disk morphology can result from planet-disk interaction, suggesting planetary presence. To date, most scattered light imaging campaigns have probed polarized light, which is only a fraction of the total light and not very sensitive to planets. To observe and characterize protoplanetary disk systems in the near-infrared in both polarized and total intensity light, we carried out an unprecedented study of scattering properties of disks, as well as of any planetary companions. Using SPHERE with star-hopping at the Very Large Telescope, we observed 29 disk hosts and their reference stars in $K_s$-band polarized light. We extracted disks in total intensity by adopting the data imputation concept with sequential non-negative matrix factorization (DI-sNMF). We obtained high-quality disk images in total intensity for 15 systems and in polarized light for 23. For well-recovered disks in polarized light and total intensity, we parameterized the polarization fraction phase functions using scaled beta distribution: the peak of polarization fraction tentatively correlates with the peak scattering angle, which could be reproduced using certain compact dust, yet more detailed modeling studies are needed. We investigated the empirical DI-sNMF detectability of disks using logistic regression: total intensity detectability of disks primarily depends on host star brightness. For disks with SPHERE data in $Y$-/$J$-/$H$-band, we summarized their polarized color at ~90 deg scattering angle: most of disks are blue in polarized $J-K_s$ color, and they are relatively redder as stellar luminosity increases, indicating larger scatterers. High-quality disk imagery in both total intensity and polarized light thus allows for disk characterization in polarization fraction, and reduces the confusion between disk and planetary signals. △ Less

Submitted 12 October, 2023; originally announced October 2023.

Comments: 25 pages, 16 figures, 3 tables, A&A accepted. Data files in FITS format will be publicly available

Journal ref: A&A 680, A114 (2023)

JWST imaging of edge-on protoplanetary disks. I. Fully vertically mixed 10$μ$m grains in the outer regions of a 1000 au disk

Authors: G. Duchene, F. Menard, K. Stapelfeldt, M. Villenave, S. G. Wolff, M. D. Perrin, C. Pinte, R. Tazaki, D. L. Padgett

Abstract: Scattered light imaging of protoplanetary disks provides key insights on the geometry and dust properties in the disk surface. Here we present JWST 2--21\,$μ$m images of a 1000\,au-radius edge-on protoplanetary disk surrounding an 0.4\,$M_\odot$ young star in Taurus, 2MASS\,J04202144+2813491. These observations represent the longest wavelengths at which a protoplanetary disk is spatially resolved… ▽ More Scattered light imaging of protoplanetary disks provides key insights on the geometry and dust properties in the disk surface. Here we present JWST 2--21\,$μ$m images of a 1000\,au-radius edge-on protoplanetary disk surrounding an 0.4\,$M_\odot$ young star in Taurus, 2MASS\,J04202144+2813491. These observations represent the longest wavelengths at which a protoplanetary disk is spatially resolved in scattered light. We combine these observations with HST optical images and ALMA continuum and CO mapping. We find that the changes in the scattered light disk morphology are remarkably small across a factor of 30 in wavelength, indicating that dust in the disk surface layers is characterized by an almost gray opacity law. Using radiative transfer models, we conclude that grains up to $\gtrsim10\,μ$m in size are fully coupled to the gas in this system, whereas grains $\gtrsim100\,μ$m are strongly settled towards the midplane. Further analyses of these observations, and similar ones of other edge-on disks, will provide strong empirical constraints on disk dynamics and evolution and grain growth models. In addition, the 7.7 and 12.\,$μ$m JWST images reveal an X-shaped feature located above the warm molecular layer traced by CO line emission. The highest elevations at which this feature is detectable roughly match the maximal extent of the disk in visible wavelength scattered light as well as of an unusual kinematic signature in CO. We propose that these phenomena could be related to a disk wind entraining small dust grains. △ Less

Submitted 13 September, 2023; originally announced September 2023.

Comments: Accepted for publication in the Astronomical Journal

Revealing magnetic field structure at the surfaces of protoplanetary disks via near-infrared circular polarization

Authors: Ilse de Langen, Ryo Tazaki

Abstract: Context. Magnetic fields play a fundamental role in the dynamical evolution of protoplanetary disks, in particular via magnetically induced disk winds. The magnetic field structure at the disk surface is crucial for driving the disk winds; however, it is still poorly understood observationally. Aims. We explore a new method to probe the magnetic field structure at the disk surface using near-infra… ▽ More Context. Magnetic fields play a fundamental role in the dynamical evolution of protoplanetary disks, in particular via magnetically induced disk winds. The magnetic field structure at the disk surface is crucial for driving the disk winds; however, it is still poorly understood observationally. Aims. We explore a new method to probe the magnetic field structure at the disk surface using near-infrared (NIR) circular polarization. Near-infrared circular polarization arises when unpolarized stellar light is scattered by magnetically aligned grains at the disk surface. In this study, we aim to clarify to what extent the observed circular polarization pattern can be used to diagnose the magnetic field structure. Methods. We first calculated light scattering properties of aligned spheroids, and the results were then used to create expected observational images of the degree of circular polarization at a NIR wavelength. Results. Magnetically aligned grains can produce circular polarization, particularly when the field configuration deviates from a purely toroidal field. We find that disk azimuthal dependence of the degree of circular polarization tends to exhibit a double peaked profile when the field structure is favorable for driving disk winds by centrifugal force. We also find that even if the disk is spatially unresolved, a net circular polarization can possibly be nonzero. We also show that the amplitude of circular polarization is strongly dependent on grain composition and axis ratio. Conclusions. Our results suggest that circular polarization observations would be useful to study the magnetic field structure and dust properties at the disk surface. △ Less

Submitted 23 February, 2023; originally announced February 2023.

Comments: 6 pages, 6 figures; Published in A&A

Journal ref: A&A, Volume 670, February 2023, A168

Fractal aggregates of sub-micron-sized grains in the young planet-forming disk around IM Lup

Authors: Ryo Tazaki, Christian Ginski, Carsten Dominik

Abstract: Despite rapidly growing disk observations, it remains a mystery what primordial dust aggregates look like and what the physical and chemical properties of their constituent grains (monomers) are in young planet-forming disks. Confrontation of models with observations to answer this mystery has been a notorious task because we have to abandon a commonly used assumption, perfectly spherical grains,… ▽ More Despite rapidly growing disk observations, it remains a mystery what primordial dust aggregates look like and what the physical and chemical properties of their constituent grains (monomers) are in young planet-forming disks. Confrontation of models with observations to answer this mystery has been a notorious task because we have to abandon a commonly used assumption, perfectly spherical grains, and take into account particles with complex morphology. In this Letter, we present the first thorough comparison between near-infrared scattered light of the young planet-forming disk around IM Lup and the light-scattering properties of complex-shaped dust particles. The availability of scattering observations at multiple wavelengths and over a significant range of scattering angles allows for the first determination of the monomer size, fractal dimension, and size of dust aggregates in a planet-forming disk. We show that the observations are best explained by fractal aggregates with a fractal dimension of 1.5 and a characteristic radius larger than $\sim2~μ$m. We also determined the radius of the monomer to be $\sim200$ nm, and monomers much smaller than this size can be ruled out on the premise that the fractal dimension is less than 2. Furthermore, dust composition comprising amorphous carbon is found to be favorable to simultaneously account for the faint scattered light and the flared disk morphology. Our results support that planet formation begins with fractal coagulation of sub-micron-sized grains. All the optical properties of complex dust particles computed in this study are publicly available. △ Less

Submitted 9 February, 2023; v1 submitted 2 February, 2023; originally announced February 2023.

Comments: 19 pages, 11 figures; Accepted for publication in ApJ Letters. The optical properties of dust particles are available at https://github.com/rtazaki1205/AggScatVIR . Updated the caption of fig. 8 to improve clarity

Observed polarized scattered light phase functions of planet-forming disks

Authors: Christian Ginski, Ryo Tazaki, Carsten Dominik, Tomas Stolker

Abstract: Dust particles are the building blocks from which planetary bodies are made. A major goal of the studies of planet-forming disks is to constrain the properties of dust particles and aggregates in order to trace their origin, structure, and the associated growth and mixing processes in the disk. Observations of scattering and/or emission of dust in a location of the disk often lead to degenerate in… ▽ More Dust particles are the building blocks from which planetary bodies are made. A major goal of the studies of planet-forming disks is to constrain the properties of dust particles and aggregates in order to trace their origin, structure, and the associated growth and mixing processes in the disk. Observations of scattering and/or emission of dust in a location of the disk often lead to degenerate information about the kind of particles, such as size, porosity, or fractal dimension of aggregates. Progress can be made by deriving the full (polarizing) scattering phase function of such particles at multiple wavelengths. This has now become possible by careful extraction from scattered light images. Such an extraction requires knowledge about the shape of the scattering surface in the disk and we discuss how to obtain such knowledge as well as the associated uncertainties. We use a sample of disk images from observations with VLT/SPHERE to, for the first time, extract the phase functions of a whole sample of disks with broad phase angle coverage. We find that polarized phase functions come in two categories. Comparing the extracted functions with theoretical predictions from rigorous T-Matrix computations of aggregates, we show that one category can be linked back to fractal, porous aggregates, while the other is consistent with more compact, less porous aggregates. We speculate that the more compact particles become visible in disks where embedded planets trigger enhanced vertical mixing. △ Less

Submitted 11 January, 2023; originally announced January 2023.

Comments: Submitted to ApJ with parallel submission to arxiv for logistical reasons; We will update the arxiv submission with the final version once the manuscript was refereed and accepted by ApJ. Data files for all phase functions will be made available once refereeing is concluded. We welcome comments from the community

How large are the monomers of dust aggregates in planet-forming disks?: Insights from quantitative optical and near-infrared polarimetry

Authors: Ryo Tazaki, Carsten Dominik

Abstract: Context: The size of the constituent particles (monomers) of dust aggregates is one of the most uncertain parameters directly affecting collisional growth of aggregates in planet-forming disks. Despite its importance, the monomer size has not yet been meaningfully constrained by disk observations. Aims: We attempt to derive the monomer size from optical and near-infrared (IR) polarimetric observat… ▽ More Context: The size of the constituent particles (monomers) of dust aggregates is one of the most uncertain parameters directly affecting collisional growth of aggregates in planet-forming disks. Despite its importance, the monomer size has not yet been meaningfully constrained by disk observations. Aims: We attempt to derive the monomer size from optical and near-infrared (IR) polarimetric observations of planet-forming disks. Methods: We perform a comprehensive parameter survey on the degree of linear polarization of light scattered by dust aggregates, using an exact numerical method called the $T$-matrix method. We investigate the effect of the monomer size, aggregate size, porosity, and composition on the degree of polarization. The obtained results are then compared with observed polarization fractions of several planet-forming disks at optical and near-IR wavelengths. Results: It is shown that the degree of polarization of aggregates depends sensitively on the monomer size unless the monomer size parameter is smaller than one or two. Comparing the simulation results with the disk observations, we find that the monomer radius is no greater than $0.4~μ$m. The inferred monomer size is therefore similar to subunit sizes of the solar system dust aggregates and the maximum size of interstellar grains. Conclusions: Optical and near-IR quantitative polarimetry will provide observational grounds on the initial conditions for dust coagulation and thereby planetesimal formation in planet-forming disks. △ Less

Submitted 18 April, 2022; originally announced April 2022.

Comments: 17 pages, 12 figures, 1 table; Accepted for publication in A&A

The water-ice feature in near-infrared disk-scattered light around HD 142527: Micron-sized icy grains lifted up to the disk surface?

Authors: Ryo Tazaki, Koji Murakawa, Takayuki Muto, Mitsuhiko Honda, Akio K. Inoue

Abstract: We study the $3~μ$m scattering feature of water ice detected in the outer disk of HD 142527 by performing radiative transfer simulations. We show that an ice mass abundance at the outer disk surface of HD 142527 is much lower than estimated in a previous study. It is even lower than inferred from far-infrared ice observations, implying ice disruption at the disk surface. Next, we demonstrate that… ▽ More We study the $3~μ$m scattering feature of water ice detected in the outer disk of HD 142527 by performing radiative transfer simulations. We show that an ice mass abundance at the outer disk surface of HD 142527 is much lower than estimated in a previous study. It is even lower than inferred from far-infrared ice observations, implying ice disruption at the disk surface. Next, we demonstrate that a polarization fraction of disk-scattered light varies across the ice-band wavelengths depending on ice grain properties; hence, polarimetric spectra would be another tool for characterizing water-ice properties. Finally, we argue that the observed reddish disk-scattered light is due to grains with a few microns in size. To explain the presence of such grains at the disk surface, we need a mechanism that can efficiently oppose dust settling. If we assume turbulent mixing, our estimate requires $α\gtrsim2\times10^{-3}$, where $α$ is a non-dimensional parameter describing the vertical diffusion coefficient of grains. Future observations probing gas kinematics would be helpful to elucidate vertical grain dynamics in the outer disk of HD 142527. △ Less

Submitted 19 August, 2021; originally announced August 2021.

Comments: 21 pages, 14 figures, 1 table; Accepted for publication in ApJ

The formation of planetary systems with SPICA

Authors: I. Kamp, M. Honda, H. Nomura, M. Audard, D. Fedele, L. B. F. M. Waters, Y. Aikawa, A. Banzatti, J. E. Bowey, M. Bradford, C. Dominik, K. Furuya, E. Habart, D. Ishihara, D. Johnstone, G. Kennedy, M. Kim, Q. Kral, S. P. Lai, B. Larsson, M. McClure, A. Miotello, M. Momose, T. Nakagawa, D. Naylor , et al. (16 additional authors not shown)

Abstract: In this era of spatially resolved observations of planet forming disks with ALMA and large ground-based telescopes such as the VLT, Keck and Subaru, we still lack statistically relevant information on the quantity and composition of the material that is building the planets, such as the total disk gas mass, the ice content of dust, and the state of water in planetesimals. SPICA is an infrared spac… ▽ More In this era of spatially resolved observations of planet forming disks with ALMA and large ground-based telescopes such as the VLT, Keck and Subaru, we still lack statistically relevant information on the quantity and composition of the material that is building the planets, such as the total disk gas mass, the ice content of dust, and the state of water in planetesimals. SPICA is an infrared space mission concept developed jointly by JAXA and ESA to address these questions. The key unique capabilities of SPICA that enable this research are (1) the wide spectral coverage 10-220 micron, (2) the high line detection sensitivity of (1-2) 10-19 W m-2 with R~2000-5000 in the far-IR (SAFARI) and 10-20 W m-2 with R~29000 in the mid-IR (SMI, spectrally resolving line profiles), (3) the high far-IR continuum sensitivity of 0.45 mJy (SAFARI), and (4) the observing efficiency for point source surveys. This paper details how mid- to far-IR infrared spectra will be unique in measuring the gas masses and water/ice content of disks and how these quantities evolve during the planet forming period. These observations will clarify the crucial transition when disks exhaust their primordial gas and further planet formation requires secondary gas produced from planetesimals. The high spectral resolution mid-IR is also unique for determining the location of the snowline dividing the rocky and icy mass reservoirs within the disk and how the divide evolves during the build-up of planetary systems. Infrared spectroscopy (mid- to far-IR) of key solid state bands is crucial for assessing whether extensive radial mixing, which is part of our Solar System history, is a general process occurring in most planetary systems and whether extrasolar planetesimals are similar to our Solar System comets/asteroids. ... (abbreviated) △ Less

Submitted 25 June, 2021; originally announced June 2021.

Comments: accepted for publication in PASA

Analytic expressions for geometric cross-sections of fractal dust aggregates

Authors: Ryo Tazaki

Abstract: In protoplanetary discs and planetary atmospheres, dust grains coagulate to form fractal dust aggregates. The geometric cross-section of these aggregates is a crucial parameter characterizing aerodynamical friction, collision rates, and opacities. However, numerical measurements of the cross-section are often time-consuming as aggregates exhibit complex shapes. In this study, we derive a novel ana… ▽ More In protoplanetary discs and planetary atmospheres, dust grains coagulate to form fractal dust aggregates. The geometric cross-section of these aggregates is a crucial parameter characterizing aerodynamical friction, collision rates, and opacities. However, numerical measurements of the cross-section are often time-consuming as aggregates exhibit complex shapes. In this study, we derive a novel analytic expression for geometric cross-sections of fractal aggregates. If an aggregate consists of $N$ monomers of radius $R_0$, its geometric cross-section $G$ is expressed as \begin{equation} \frac{G}{NπR_0^2}=\frac{A}{1+(N-1)\tildeσ}, \nonumber \end{equation} where $\tildeσ$ is an overlapping efficiency, and $A$ is a numerical factor connecting the analytic expression to the small non-fractal cluster limit. The overlapping efficiency depends on the fractal dimension, fractal prefactor, and $N$ of the aggregate, and its analytic expression is derived as well. The analytic expressions successfully reproduce numerically measured cross-sections of aggregates. We also find that our expressions are compatible with the mean-field light scattering theory of aggregates in the geometrical optics limit. The analytic expressions greatly simplify an otherwise tedious calculation and will be useful in model calculations of fractal grain growth in protoplanetary discs and planetary atmospheres. △ Less

Submitted 14 April, 2021; originally announced April 2021.

Comments: 12 pages, 8 figures, 1 table; Accepted for publication in MNRAS

Scattering polarization of 3-$μ$m water-ice feature by large icy grains

Authors: Ryo Tazaki, Koji Murakawa, Takayuki Muto, Mitsuhiko Honda, Akio K. Inoue

Abstract: Water ice has a strong spectral feature at a wavelength of approximately $3~μ$m, which plays a vital role in our understanding of the icy universe. In this study, we investigate the scattering polarization of this water-ice feature. The linear polarization degree of light scattered by $μ$m-sized icy grains is known to be enhanced at the ice band; however, the dependence of this polarization enhanc… ▽ More Water ice has a strong spectral feature at a wavelength of approximately $3~μ$m, which plays a vital role in our understanding of the icy universe. In this study, we investigate the scattering polarization of this water-ice feature. The linear polarization degree of light scattered by $μ$m-sized icy grains is known to be enhanced at the ice band; however, the dependence of this polarization enhancement on various grain properties is unclear. We find that the enhanced polarization at the ice band is sensitive to the presence of $μ$m-sized grains as well as their ice abundance. We demonstrate that this enhancement is caused by the high absorbency of the water-ice feature, which attenuates internal scattering and renders the surface reflection dominant over internal scattering. Additionally, we compare our models with polarimetric observations of the low-mass protostar L1551 IRS 5. Our results show that scattering by a maximum grain radius of a few microns with a low water-ice abundance is consistent with observations. Thus, scattering polarization of the water-ice feature is a useful tool for characterizing ice properties in various astronomical environments. △ Less

Submitted 19 January, 2021; originally announced January 2021.

Comments: 19 pages, 17 figures, 1 table; Accepted for publication in ApJ

Haze Formation on Triton

Authors: Kazumasa Ohno, Xi Zhang, Ryo Tazaki, Satoshi Okuzumi

Abstract: The largest moon of Neptune, Triton, possess a cold and hazy atmosphere. Since the discovery of near-surface haze layer during the Voyager fly in 1989, the haze formation mechanism has not been investigated in detail. Here, we provide the first haze microphysical model on Triton. Our model solves the evolution of both size and porosity distributions of haze particles in a self-consistent manner. W… ▽ More The largest moon of Neptune, Triton, possess a cold and hazy atmosphere. Since the discovery of near-surface haze layer during the Voyager fly in 1989, the haze formation mechanism has not been investigated in detail. Here, we provide the first haze microphysical model on Triton. Our model solves the evolution of both size and porosity distributions of haze particles in a self-consistent manner. We simulated the formation of sphere and aggregate hazes with and without condensation of the C$_2$H$_4$ ice. The haze particles can grow into fractal aggregates with mass-equivalent sphere sizes of $\sim0.1$--$1~{\rm μm}$ and fractal dimension of $D_{\rm f} = 1.8$--$2.2$. The ice-free hazes cannot simultaneously explain both UV and visible observations of Voyager 2, while including the condensation of C$_2$H$_4$ ices provides two better solutions. For ice aggregates, the required total haze mass flux is $\sim2\times{10}^{-15}~{\rm g~{cm}^{-2}~s^{-1}}$. For the icy sphere scenario, the column integrated C$_2$H$_4$ production rate is $\sim8\times{10}^{-15}~{\rm g~{cm}^{-2}~s^{-1}}$, and the ice-free mass flux of $\sim6\times{10}^{-17}~{\rm g~{cm}^{-2}~s^{-1}}$. The UV occultation observations at short wavelength $<0.15~{\rm μm}$ may slightly favor the icy aggregates. Observations of the haze optical depth and the degree of forward scattering in UV and visible should be able to distinguish whether Triton's hazes are icy spheres or ice aggregates in future Triton missions. △ Less

Submitted 16 April, 2021; v1 submitted 22 December, 2020; originally announced December 2020.

Comments: 29 pages, 14 figures, accepted for publication in ApJ, several typos had been amended in the proof of correction

Distortion of Magnetic Fields in the Dense Core SL42 (CrA-E) in the Corona Australis Molecular Cloud Complex

Authors: Ryo Kandori, Motohide Tamura, Masao Saito, Kohji Tomisaka, Tomoaki Matsumoto, Ryo Tazaki, Tetsuya Nagata, Nobuhiko Kusakabe, Yasushi Nakajima, Jungmi Kwon, Takahiro Nagayama, Ken'ichi Tatematsu

Abstract: Detailed magnetic field structure of the dense core SL42 (CrA-E) in the Corona Australis molecular cloud complex was investigated based on near-infrared polarimetric observations of background stars to measure dichroically polarized light produced by magnetically aligned dust grains. The magnetic fields in and around SL42 were mapped using 206 stars and curved magnetic fields were identified. On t… ▽ More Detailed magnetic field structure of the dense core SL42 (CrA-E) in the Corona Australis molecular cloud complex was investigated based on near-infrared polarimetric observations of background stars to measure dichroically polarized light produced by magnetically aligned dust grains. The magnetic fields in and around SL42 were mapped using 206 stars and curved magnetic fields were identified. On the basis of simple hourglass (parabolic) magnetic field modeling, the magnetic axis of the core on the plane of sky was estimated to be $40^{\circ} \pm 3^{\circ}$. The plane-of-sky magnetic field strength of SL42 was found to be $22.4 \pm 13.9$ $μ$G. Taking into account the effects of thermal/turbulent pressure and the plane-of-sky magnetic field component, the critical mass of SL42 was obtained to be $M_{\rm cr} = 21.2 \pm 6.6$ M$_{\odot}$, which is close to the observed core mass of $M_{\rm core} \approx 20$ M$_{\odot}$. We thus conclude that SL42 is in a condition close to the critical state if the magnetic fields lie near the plane of the sky. Since there is a very low luminosity object (VeLLO) toward the center of SL42, it is unlikely this core is in a highly subcritical condition (i.e., magnetic inclination angle significantly deviated from the plane of sky). The core probably started to collapse from a nearly kinematically critical state. In addition to the hourglass magnetic field modeling, the Inoue \& Fukui (2013) mechanism may explain the origin of the curved magnetic fields in the SL42 region. △ Less

Submitted 28 July, 2020; originally announced July 2020.

Comments: 23 pages, 12 figures, accepted for publication in ApJ. arXiv admin note: text overlap with arXiv:1912.10029

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

Authors: Taichi Uyama, Thayne Currie, Valentin Christiaens, Jaehan Bae, Takayuki Muto, Sanemichi Z. Takahashi, Ryo Tazaki, Marie Ygouf, Jeremy N. Kasdin, Tyler Groff, Timothy D. Brandt, Jeffrey Chilcote, Masahiko Hayashi, Michael W. McElwain, Olivier Guyon, Julien Lozi, Nemanja Jovanovic, Frantz Martinache, Tomoyuki Kudo, Motohide Tamura, Eiji Akiyama, Charles A. Beichman, Carol A. Grady, Gillian R. Knapp, Jungmi Kwon , et al. (5 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 22 July, 2020; originally announced July 2020.

Comments: 17 pages, 14 figures, accepted for publication in ApJ

Distortion of Magnetic Fields in BHR 71

Authors: Ryo Kandori, Motohide Tamura, Masao Saito, Kohji Tomisaka, Tomoaki Matsumoto, Ryo Tazaki, Tetsuya Nagata, Nobuhiko Kusakabe, Yasushi Nakajima, Jungmi Kwon, Takahiro Nagayama, Ken'ichi Tatematsu

Abstract: The magnetic field structure of a star-forming Bok globule BHR 71 was determined based on near-infrared polarimetric observations of background stars. The magnetic field in BHR 71 was mapped from 25 stars. By using a simple 2D parabolic function, the plane-of-sky magnetic axis of the core was found to be $θ_{\rm mag} = 125^{\circ} \pm 11^{\circ}$. The plane-of-sky mean magnetic field strength of B… ▽ More The magnetic field structure of a star-forming Bok globule BHR 71 was determined based on near-infrared polarimetric observations of background stars. The magnetic field in BHR 71 was mapped from 25 stars. By using a simple 2D parabolic function, the plane-of-sky magnetic axis of the core was found to be $θ_{\rm mag} = 125^{\circ} \pm 11^{\circ}$. The plane-of-sky mean magnetic field strength of BHR 71 was found to be $B_{\rm pos} = 8.8 - 15.0$ $μ$G, indicating that the BHR 71 core is magnetically supercritical with $λ= 1.44 - 2.43$. Taking into account the effect of thermal/turbulent pressure and the plane-of-sky magnetic field component, the critical mass of BHR 71 was $M_{\rm cr} = 14.5-18.7$ M$_{\odot}$, which is consistent with the observed core mass of $M_{\rm core} \approx 14.7$ M$_{\odot}$ (Yang et al. 2017). We conclude that BHR 71 is in a condition close to a kinematically critical state, and the magnetic field direction lies close to the plane of sky. Since BHR 71 is a star-forming core, a significantly subcritical condition (i.e., the magnetic field direction deviating from the plane of sky) is unlikely, and collapsed from a condition close to a kinematically critical state. There are two possible scenarios to explain the curved magnetic fields of BHR 71, one is an hourglass-like field structure due to mass accumulation and the other is the Inoue \& Fukui (2013) mechanism, which proposes the interaction of the core with a shock wave to create curved magnetic fields wrapping around the core. △ Less

Submitted 27 February, 2020; originally announced February 2020.

Comments: Accepted for publication in ApJ

Dust Destruction by Charging: A Possible Origin of Grey Extinction Curves of Active Galactic Nuclei

Authors: Ryo Tazaki, Kohei Ichikawa, Mitsuru Kokubo

Abstract: Observed extinction curves of active galactic nuclei (AGNs) are significantly different from those observed in the Milky Way. The observations require preferential removal of small grains at the AGN environment; however, the physics for this remains unclear. In this paper, we propose that dust destruction by charging, or Coulomb explosion, may be responsible for AGN extinction curves. Harsh AGN ra… ▽ More Observed extinction curves of active galactic nuclei (AGNs) are significantly different from those observed in the Milky Way. The observations require preferential removal of small grains at the AGN environment; however, the physics for this remains unclear. In this paper, we propose that dust destruction by charging, or Coulomb explosion, may be responsible for AGN extinction curves. Harsh AGN radiation makes a dust grain highly charged through photoelectric emission, and grain fission via the Coulomb explosion occurs when the electrostatic tensile stress of a charge grain exceeds its tensile strength. We show that the Coulomb explosion can preferentially remove both small silicate and graphite grains and successfully reproduce both flat extinction curves and the absence of 2175Å~bump. △ Less

Submitted 19 February, 2020; originally announced February 2020.

Comments: 8 pages, 6 Figures; Accepted for publication in ApJ

Dust Destruction by Drift-Induced Sputtering in Active Galactic Nuclei

Authors: Ryo Tazaki, Kohei Ichikawa

Abstract: Recent mid-infrared high spatial resolution observations have revealed that active galactic nuclei (AGNs) may host a polar dust region with the size of several pc, and such dust may be carried by radiation from the central engine. The polar dust emission often exhibits very weak or absence of the silicate 10-$μ$m emission feature. A possible explanation is that the polar dust is dominated by micro… ▽ More Recent mid-infrared high spatial resolution observations have revealed that active galactic nuclei (AGNs) may host a polar dust region with the size of several pc, and such dust may be carried by radiation from the central engine. The polar dust emission often exhibits very weak or absence of the silicate 10-$μ$m emission feature. A possible explanation is that the polar dust is dominated by micron-sized large grains because these grains do not show the silicate feature, while it remains unclear how large grains are preferentially supplied to the polar region. We here propose a new scenario describing the prevalence of large grains at the polar region. We show that grains are accelerated to the hypersonic drift velocity by the radiation pressure from AGN, and the hypersonic drift results in dust destruction via kinetic sputtering. Sputtering destroys small grains faster than the large ones, and thus larger grains will be preferentially blown over longer distance. Although the hypersonic drift, or kinetic sputtering, tends to be suppressed for very small grains due to the Coulomb drag, they might also be disrupted by Coulomb explosion. Removal of small grains and/or survival of large grains may explain the lack of a silicate 10-$μ$m emission feature in polar dust emission. △ Less

Submitted 5 February, 2020; originally announced February 2020.

Comments: 17 pages, 14 Figures; Accepted for publication in ApJ

Distortion of Magnetic Fields in Barnard 335

Authors: Ryo Kandori, Masao Saito, Motohide Tamura, Kohji Tomisaka, Tomoaki Matsumoto, Ryo Tazaki, Tetsuya Nagata, Nobuhiko Kusakabe, Yasushi Nakajima, Jungmi Kwon, Takahiro Nagayama, Ken'ichi Tatematsu

Abstract: In this study, the detailed magnetic field structure of the dense protostellar core Barnard 335 (B335) was revealed based on near-infrared polarimetric observations of background stars to measure dichroically polarized light produced by magnetically aligned dust grains in the core. Magnetic fields pervading B335 were mapped using 24 stars after subtracting unrelated ambient polarization components… ▽ More In this study, the detailed magnetic field structure of the dense protostellar core Barnard 335 (B335) was revealed based on near-infrared polarimetric observations of background stars to measure dichroically polarized light produced by magnetically aligned dust grains in the core. Magnetic fields pervading B335 were mapped using 24 stars after subtracting unrelated ambient polarization components, for the first time revealing that they have an axisymmetrically distorted hourglass-shaped structure toward the protostellar core. On the basis of simple two- and three-dimensional magnetic field modeling, magnetic inclination angles in the plane-of-sky and line-of-sight directions were determined to be $90^{\circ} \pm 7^{\circ}$ and $50^{\circ} \pm 10^{\circ}$, respectively. The total magnetic field strength of B335 was determined to be $30.2 \pm 17.7$ $μ{\rm G}$. The critical mass of B335, evaluated using both magnetic and thermal/turbulent support against collapse, was determined to be $M_{\rm cr} = 3.37 \pm 0.94$ ${\rm M}_{\odot}$, which is identical to the observed core mass of $M_{\rm core}=3.67$ M$_{\odot}$. We thus concluded that B335 started its contraction from a condition near equilibrium. We found a linear relationship in the polarization versus extinction diagram, up to $A_V \sim 15$ mag toward the stars with the greatest obscuration, which verified that our observations and analysis provide an accurate depiction of the core. △ Less

Submitted 21 January, 2020; originally announced January 2020.

Comments: Accepted for publication in ApJ. arXiv admin note: substantial text overlap with arXiv:1910.09121, arXiv:1912.10029

Distortion of Magnetic Fields in the Dense Core CB81 (L1774, Pipe 42) in the Pipe Nebula

Authors: Ryo Kandori, Motohide Tamura, Masao Saito, Kohji Tomisaka, Tomoaki Matsumoto, Ryo Tazaki, Tetsuya Nagata, Nobuhiko Kusakabe, Yasushi Nakajima, Jungmi Kwon, Takahiro Nagayama, Ken'ichi Tatematsu

Abstract: The detailed magnetic field structure of the starless dense core CB81 (L1774, Pipe 42) in the Pipe Nebula was determined based on near-infrared polarimetric observations of background stars to measure dichroically polarized light produced by magnetically aligned dust grains in the core. The magnetic fields pervading CB81 were mapped using 147 stars and axisymmetrically distorted hourglass-like fie… ▽ More The detailed magnetic field structure of the starless dense core CB81 (L1774, Pipe 42) in the Pipe Nebula was determined based on near-infrared polarimetric observations of background stars to measure dichroically polarized light produced by magnetically aligned dust grains in the core. The magnetic fields pervading CB81 were mapped using 147 stars and axisymmetrically distorted hourglass-like fields were identified. On the basis of simple 2D and 3D magnetic field modeling, the magnetic inclination angles in the plane-of-sky and line-of-sight directions were determined to be $4^{\circ} \pm 8^{\circ}$ and $20^{\circ} \pm 20^{\circ}$, respectively. The total magnetic field strength of CB81 was found to be $7.2 \pm 2.3$ $μ{\rm G}$. Taking into account the effects of thermal/turbulent pressure and magnetic fields, the critical mass of CB81 was calculated to be $M_{\rm cr}=4.03 \pm 0.40$ M$_{\odot}$, which is close to the observed core mass of $M_{\rm core}=3.37 \pm 0.51$ M$_{\odot}$. We thus conclude that CB81 is in a condition close to the critical state. In addition, a spatial offset of $92''$ was found between the center of magnetic field geometry and the dust extinction distribution; this offset structure could not have been produced by self-gravity. The data also indicate a linear relationship between polarization and extinction up to $A_V \sim 30$ mag going toward the core center. This result confirms that near-infrared polarization can accurately trace the overall magnetic field structure of the core. △ Less

Submitted 19 December, 2019; originally announced December 2019.

Comments: Accepted for publication in ApJ. arXiv admin note: text overlap with arXiv:1910.09121

Distortion of Magnetic Fields in a Starless Core VI: Application of Flux Freezing Model and Core Formation of FeSt 1-457

Authors: Ryo Kandori, Kohji Tomisaka, Masao Saito, Motohide Tamura, Tomoaki Matsumoto, Ryo Tazaki, Tetsuya Nagata, Nobuhiko Kusakabe, Yasushi Nakajima, Jungmi Kwon, Takahiro Nagayama, Ken'ichi Tatematsu

Abstract: Observational data for the hourglass-like magnetic field toward the starless dense core FeSt 1-457 were compared with a flux freezing magnetic field model (Myers et al. 2018). Fitting of the observed plane-of-sky magnetic field using the flux freezing model gave a residual angle dispersion comparable with the results based on a simple three-dimensional parabolic model. The best-fit parameters for… ▽ More Observational data for the hourglass-like magnetic field toward the starless dense core FeSt 1-457 were compared with a flux freezing magnetic field model (Myers et al. 2018). Fitting of the observed plane-of-sky magnetic field using the flux freezing model gave a residual angle dispersion comparable with the results based on a simple three-dimensional parabolic model. The best-fit parameters for the flux freezing model were a line-of-sight magnetic inclination angle of $γ_{\rm mag} = 35^{\circ} \pm 15^{\circ}$ and a core center to ambient (background) density contrast of $ρ_{\rm c} / ρ_{\rm bkg} = 75$. The initial density for core formation ($ρ_0$) was estimated to be $ρ_{\rm c} / 75 = 4670$ cm$^{-3}$, which is about one order of magnitude higher than the expected density ($\sim 300$ cm$^{-3}$) for the inter-clump medium of the Pipe Nebula. FeSt 1-457 is likely to have been formed from the accumulation of relatively dense gas, and the relatively dense background column density of $A_V \simeq 5$ mag supports this scenario. The initial radius (core formation radius) $R_0$ and the initial magnetic field strength $B_0$ were obtained to be 0.15 pc ($1.64 R$) and $10.8-14.6$ $μ$G, respectively. We found that the initial density $ρ_0$ is consistent with the mean density of the nearly critical magnetized filament with magnetic field strength $B_0$ and radius $R_0$. The relatively dense initial condition for core formation can be naturally understood if the origin of the core is the fragmentation of magnetized filaments. △ Less

Submitted 5 December, 2019; originally announced December 2019.

Comments: Accepted for publication in ApJ

Distortion of Magnetic Fields in Barnard 68

Authors: Ryo Kandori, Motohide Tamura, Masao Saito, Kohji Tomisaka, Tomoaki Matsumoto, Nobuhiko Kusakabe, Jungmi Kwon, Takahiro Nagayama, Tetsuya Nagata, Ryo Tazaki, Ken'ichi Tatematsu

Abstract: The magnetic field structure, kinematical stability, and evolutionary status of the starless dense core Barnard 68 (B68) are revealed based on the near-infrared polarimetric observations of background stars, measuring the dichroically polarized light produced by aligned dust grains in the core. After subtracting unrelated ambient polarization components, the magnetic fields pervading B68 are mappe… ▽ More The magnetic field structure, kinematical stability, and evolutionary status of the starless dense core Barnard 68 (B68) are revealed based on the near-infrared polarimetric observations of background stars, measuring the dichroically polarized light produced by aligned dust grains in the core. After subtracting unrelated ambient polarization components, the magnetic fields pervading B68 are mapped using 38 stars and axisymmetrically distorted hourglass-like magnetic fields are obtained, although the evidence for the hourglass field is not very strong. On the basis of simple 2D and 3D magnetic field modeling, the magnetic inclination angles on the plane-of-sky and in the line-of-sight direction are determined to be $47^{\circ} \pm 5^{\circ}$ and $20^{\circ} \pm 10^{\circ}$, respectively. The total magnetic field strength of B68 is obtained to be $26.1 \pm 8.7$ $μ{\rm G}$. The critical mass of B68, evaluated using both magnetic and thermal/turbulent support, is $M_{\rm cr} = 2.30 \pm 0.20$ ${\rm M}_{\odot}$, which is consistent with the observed core mass of $M_{\rm core}=2.1$ M$_{\odot}$, suggesting nearly critical state. We found a relatively linear relationship between polarization and extinction up to $A_V \sim 30$ mag toward the stars with deepest obscuration. Further theoretical and observational studies are required to explain the dust alignment in cold and dense regions in the core. △ Less

Submitted 20 October, 2019; originally announced October 2019.

Comments: Accepted by PASJ

Clouds of Fluffy Aggregates: How They Form in Exoplanetary Atmospheres and Influence Transmission Spectra

Authors: Kazumasa Ohno, Satoshi Okuzumi, Ryo Tazaki

Abstract: Transmission spectrum surveys have suggested the ubiquity of high-altitude clouds in exoplanetary atmospheres. Theoretical studies have investigated the formation processes of the high-altitude clouds; however, cloud particles have been commonly approximated as compact spheres, which is not always true for solid mineral particles that likely constitute exoplanetary clouds. Here, we investigate how… ▽ More Transmission spectrum surveys have suggested the ubiquity of high-altitude clouds in exoplanetary atmospheres. Theoretical studies have investigated the formation processes of the high-altitude clouds; however, cloud particles have been commonly approximated as compact spheres, which is not always true for solid mineral particles that likely constitute exoplanetary clouds. Here, we investigate how the porosity of cloud particles evolve in exoplanetary atmospheres and influence the cloud vertical profiles. We first construct a porosity evolution model that takes into account the fractal aggregation and the compression of cloud particle aggregates. Using a cloud microphysical model coupled with the porosity model, we demonstrate that the particle internal density can significantly decrease during the cloud formation. As a result, fluffy-aggregate clouds ascend to altitude much higher than that for compact-sphere clouds assumed so far. We also examine how the fluffy-aggregate clouds affect transmission spectra. We find that the clouds largely obscure the molecular features and produce a spectral slope originated by the scattering properties of aggregates. Finally, we compare the synthetic spectra with the observations of GJ1214 b and find that its flat spectrum could be explained if the atmospheric metallicity is sufficiently high ($\ge100\times$ solar) and the monomer size is sufficiently small ($r_{\rm mon}<1~{\rm μm}$). The high-metallicity atmosphere may offer the clues to explore the gas accretion processes onto past GJ1214b. △ Less

Submitted 20 January, 2020; v1 submitted 6 August, 2019; originally announced August 2019.

Comments: 23 pages, 12 figures, accepted for publication in ApJ. v3: Figures 5-9 have been corrected

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

Authors: Ryo Tazaki, Hidekazu Tanaka, Akimasa Kataoka, Satoshi Okuzumi, Takayuki Muto

Abstract: 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 i… ▽ More 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. △ Less

Submitted 19 September, 2019; v1 submitted 29 June, 2019; originally announced July 2019.

Comments: 17 pages, 10 Figures, 1 Table, accepted for publication in ApJ

Nonsticky Ice at the Origin of the Uniformly Polarized Submillimeter Emission from the HL Tau Disk

Authors: Satoshi Okuzumi, Ryo Tazaki

Abstract: Recent (sub)millimeter polarimetric observations toward the young star HL Tau have successfully detected polarization emission from its circumstellar disk. The polarization pattern observed at 0.87 mm is uniform and parallel to the disk's minor axis, consistent with the self-scattering of thermal emission by dust particles whose maximum radius is $\approx 100~μm$. However, this maximum size is con… ▽ More Recent (sub)millimeter polarimetric observations toward the young star HL Tau have successfully detected polarization emission from its circumstellar disk. The polarization pattern observed at 0.87 mm is uniform and parallel to the disk's minor axis, consistent with the self-scattering of thermal emission by dust particles whose maximum radius is $\approx 100~μm$. However, this maximum size is considerably smaller than anticipated from dust evolution models that assume a high sticking efficiency for icy particles. Here, we propose that the unexpectedly small particle size can be explained if CO$_2$ ice covers the particles in the outer region of the HL Tau disk. CO$_2$ ice is one of the most major interstellar ices, and laboratory experiments show that it is poorly sticky. Based on dust evolution models accounting for CO$_2$ ice mantles as well as aggregate sintering and post-processing radiative transfer, we simulate the polarimetric observation of HL Tau at 0.87 mm. We find that the models with CO$_2$ ice mantles better match the observation. These models also predict that only particles lying between the H$_2$O and CO$_2$ snow lines can grow to millimeter to centimeter sizes, and that their rapid inward drift results in a local dust gap similar to the 10 au gap of the HL Tau disk. We also suggest that the millimeter spectral index for the outer part of the HL Tau disk is largely controlled by the optical thickness of this region and does not necessarily indicate dust growth to millimeter sizes. △ Less

Submitted 8 May, 2019; v1 submitted 8 April, 2019; originally announced April 2019.

Comments: 15 pages, 6 figures, accepted for publication in ApJ

Polarization in Disks

Authors: Ian Stephens, Zhi-Yun Li, Haifeng Yang, Akimasa Kataoka, Leslie W. Looney, Charles L. H. Hull, Manuel Fernández-López, Sarah I. Sadavoy, Woojin Kwon, Satoshi Ohashi, Ryo Tazaki, Dan Li, Thiem Hoang, Gesa H. -M. Bertrang, Carlos Carrasco-González, William R. F. Dent, Satoko Takahashi, Francesca Bacciotti, Felipe O. Alves, Josep M. Girart, Qizhou Zhang, Ramprasad Rao, Adriana Pohl, Marco Padovani, Daniele Galli , et al. (2 additional authors not shown)

Abstract: Polarized dust emission outside of disks reveal the magnetic field morphology of molecular clouds. Within disks, however, polarized dust emission can arise from very different mechanisms (e.g., self-scattering), and each of them are useful for constraining physical properties in the disk. For example, these mechanisms allow us to constrain the disk grain size distributions and grain/disk geometrie… ▽ More Polarized dust emission outside of disks reveal the magnetic field morphology of molecular clouds. Within disks, however, polarized dust emission can arise from very different mechanisms (e.g., self-scattering), and each of them are useful for constraining physical properties in the disk. For example, these mechanisms allow us to constrain the disk grain size distributions and grain/disk geometries, independent from current methods of measuring these parameters. To accurately model these features and disentangle the various polarization mechanisms, multiwavelength observations at very high resolution and sensitivity are required. With significant upgrades to current interferometric facilities, we can understand how grains evolve in disks during the planet formation process. △ Less

Submitted 13 March, 2019; originally announced March 2019.

Comments: Science white paper submitted to the Astro2020 Decadal Survey

Millimeter-wave polarization due to grain alignment by the gas flow in protoplanetary disks

Authors: Akimasa Kataoka, Satoshi Okuzumi, Ryo Tazaki

Abstract: Dust grains emit intrinsic polarized emission if they are elongated and aligned in the same direction. The direction of the grain alignment is determined by external forces, such as magnetic fields, radiation, and gas flow against the dust grains. In this letter, we apply the concept of the grain alignment by gas flow, which is called mechanical alignment, to the situation of a protoplanetary disk… ▽ More Dust grains emit intrinsic polarized emission if they are elongated and aligned in the same direction. The direction of the grain alignment is determined by external forces, such as magnetic fields, radiation, and gas flow against the dust grains. In this letter, we apply the concept of the grain alignment by gas flow, which is called mechanical alignment, to the situation of a protoplanetary disk. We assume that grains have a certain helicity, which results in the alignment with the minor axis parallel to the grain velocity against the ambient disk gas and discuss the morphology of polarization vectors in a protoplanetary disk. We find that the direction of the polarization vectors depends on the Stokes number, which denotes how well grains are coupled to the gas. If the Stokes number is less than unity, orientation of polarization is in the azimuthal direction since the dust velocity against the gas is in the radial direction. If the Stokes number is as large as unity, the polarization vectors show a leading spiral pattern since the radial and azimuthal components of the gas velocity against the dust grains are comparable. This suggests that if the observed polarization vectors show a leading spiral pattern, it would indicate that Stokes number of dust grains is around unity, which is presumably radially drifting. △ Less

Submitted 8 March, 2019; originally announced March 2019.

Comments: 7 pages, 5 figures, accepted for publication in The Astrophysical Journal Letters

Effect of dust size and structure on scattered light images of protoplanetary discs

Authors: Ryo Tazaki, Hidekazu Tanaka, Takayuki Muto, Akimasa Kataoka, Satoshi Okuzumi

Abstract: We study scattered light properties of protoplanetary discs at near-infrared wavelengths for various dust size and structure by performing radiative transfer simulations. We show that different dust structures might be probed by measuring disk polarisation fraction as long as the dust radius is larger than the wavelength. When the radius is larger than the wavelength, disc scattered light will be… ▽ More We study scattered light properties of protoplanetary discs at near-infrared wavelengths for various dust size and structure by performing radiative transfer simulations. We show that different dust structures might be probed by measuring disk polarisation fraction as long as the dust radius is larger than the wavelength. When the radius is larger than the wavelength, disc scattered light will be highly polarised for highly porous dust aggregates, whereas more compact dust structure tends to show low polarisation fraction. Next, roles of monomer radius and fractal dimension for scattered light colours are studied. We find that, outside the Rayleigh regime, as fractal dimension or monomer radius increases, colours of the effective albedo at near-infrared wavelengths vary from blue to red. Our results imply that discs showing grey or slightly blue colours and high polarisation fraction in near-infrared wavelengths might be explained by the presence of large porous aggregates containing sub-microns sized monomers. △ Less

Submitted 5 March, 2019; originally announced March 2019.

Comments: Accepted for publication in MNRAS, 18 pages, 19 figures

Hunting Axion Dark Matter with Protoplanetary Disk Polarimetry

Authors: Tomohiro Fujita, Ryo Tazaki, Kenji Toma

Abstract: We find that the polarimetric observations of protoplanetary disks are useful to search for ultra-light axion dark matter. Axion dark matter predicts the rotation of the linear polarization plane of propagating light, and protoplanetary disks are ideal targets to observe it. We show that a recent observation puts the tightest constraint on the axion-photon coupling constant for axion mass… ▽ More We find that the polarimetric observations of protoplanetary disks are useful to search for ultra-light axion dark matter. Axion dark matter predicts the rotation of the linear polarization plane of propagating light, and protoplanetary disks are ideal targets to observe it. We show that a recent observation puts the tightest constraint on the axion-photon coupling constant for axion mass $m\lesssim10^{-21}$eV. △ Less

Submitted 1 May, 2019; v1 submitted 8 November, 2018; originally announced November 2018.

Comments: 5 pages, 3 figures, Updated to match the version accepted by PRL

Journal ref: Phys. Rev. Lett. 122, 191101 (2019)

Distortion of Magnetic Fields in a Starless Core V: Near-infrared and Submillimeter Polarization in FeSt 1-457

Authors: Ryo Kandori, Tetsuya Nagata, Ryo Tazaki, Motohide Tamura, Masao Saito, Kohji Tomisaka, Tomoaki Matsumoto, Nobuhiko Kusakabe, Yasushi Nakajima, Jungmi Kwon, Takahiro Nagayama, Ken'ichi Tatematsu

Abstract: The relationship between submillimeter (submm) dust emission polarization and near-infrared (NIR) $H$-band polarization produced by dust dichroic extinction was studied for the cold starless dense core FeSt 1-457. The distribution of polarization angles ($90^{\circ}$-rotated for submm) and degrees were found to be very different between at submm and NIR wavelengths. The mean polarization angles fo… ▽ More The relationship between submillimeter (submm) dust emission polarization and near-infrared (NIR) $H$-band polarization produced by dust dichroic extinction was studied for the cold starless dense core FeSt 1-457. The distribution of polarization angles ($90^{\circ}$-rotated for submm) and degrees were found to be very different between at submm and NIR wavelengths. The mean polarization angles for FeSt 1-457 at submm and NIR wavelengths are $132.1^{\circ} \pm 22.0^{\circ}$ and $2.7^{\circ} \pm 16.2^{\circ}$, respectively. The correlation between $P_H$ and $A_V$ was found to be linear from outermost regions to relatively dense line of sight of $A_V \approx 25$ mag, indicating that NIR polarization reflects overall polarization (magnetic field) structure of the core at least in this density range. The flat $P_H/A_V$ versus $A_V$ correlations were confirmed, and the polarization efficiency was found to be comparable to the observational upper limit (Jones 1989). On the other hand, as reported by Alves et al., submm polarization degrees show clear linearly decreasing trend against $A_V$ from $A_V \approx 20$ mag to the densest center ($A_V \approx 41$ mag), appearing as "polarization hole" structure. The power law index for the $P_{\rm submm}$ versus $A_V$ relationship was obtained to be $\approx -1$, indicating that the alignment for the submm sensitive dust is lost. These very different polarization distributions at submm and NIR wavelengths suggest that (1) there is different radiation environment at these wavelengths or (2) submm-sensitive dust is localized or the combination of them. △ Less

Submitted 14 October, 2018; originally announced October 2018.

Comments: Accepted to the Astrophysical Journal (ApJ)

Light Scattering by Fractal Dust Aggregates. II. Opacity and Asymmetry Parameter

Authors: Ryo Tazaki, Hidekazu Tanaka

Abstract: Optical properties of dust aggregates are important at various astrophysical environments. To find a reliable approximation method for optical properties of dust aggregates, we calculate the opacity and the asymmetry parameter of dust aggregates by using a rigorous numerical method, the T-Matrix Method (TMM), and then the results are compared to those obtained by approximate methods; the Rayleigh-… ▽ More Optical properties of dust aggregates are important at various astrophysical environments. To find a reliable approximation method for optical properties of dust aggregates, we calculate the opacity and the asymmetry parameter of dust aggregates by using a rigorous numerical method, the T-Matrix Method (TMM), and then the results are compared to those obtained by approximate methods; the Rayleigh-Gans-Debye (RGD) theory, the effective medium theory (EMT), and the distribution of hollow spheres method (DHS). First of all, we confirm that the RGD theory breaks down when multiple scattering is important. In addition, we find that both EMT and DHS fail to reproduce the optical properties of dust aggregates with fractal dimension of 2 when the incident wavelength is shorter than the aggregate radius. In order to solve these problems, we test the mean field theory (MFT), where multiple scattering can be taken into account. We show that the extinction opacity of dust aggregates can be well reproduced by MFT. However, it is also shown that MFT is not able to reproduce the scattering and absorption opacities when multiple scattering is important. We successfully resolve this weak point of MFT, by newly developing a modified mean field theory (MMF). Hence, we conclude that MMF can be a useful tool to investigate radiative transfer properties of various astrophysical environments. We also point out an enhancement of the absorption opacity of dust aggregates in the Rayleigh domain, which would be important to explain the large millimeter-wave opacity inferred from observations of protoplanetary disks. △ Less

Submitted 19 June, 2018; v1 submitted 10 March, 2018; originally announced March 2018.

Comments: version accepted in ApJ

Journal ref: The Astrophysical Journal, 860:79 (2018)

Cooling timescale of dust tori in dying active galactic nuclei

Authors: Kohei Ichikawa, Ryo Tazaki

Abstract: We estimate the dust torus cooling timescale once the active galactic nucleus (AGN) is quenched. In a clumpy torus system, once the incoming photons are suppressed, the cooling timescale of one clump from $T_{\rm dust}=1000$ K to several $10$ K is less than $10$ years, indicating that the dust torus cooling time is mainly governed by the light crossing time of the torus from the central engine. Af… ▽ More We estimate the dust torus cooling timescale once the active galactic nucleus (AGN) is quenched. In a clumpy torus system, once the incoming photons are suppressed, the cooling timescale of one clump from $T_{\rm dust}=1000$ K to several $10$ K is less than $10$ years, indicating that the dust torus cooling time is mainly governed by the light crossing time of the torus from the central engine. After considering the light crossing time of the torus, the AGN torus emission at $12~μ$m becomes over two orders of magnitude fainter within $100$ years after the quenching. We also propose that those "dying" AGN could be found using the AGN indicators with different physical scale $R$ such as $12~μ$m band luminosity tracing AGN torus ($R \sim 10$ pc) and the optical [OIII]$\lambda5007$ emission line narrow line regions ($R=10^{2-4}$ pc). △ Less

Submitted 9 June, 2017; originally announced June 2017.

Comments: 8 pages, 4 figures, accepted for publication in ApJ

Radiative grain alignment in protoplanetary disks: Implications for polarimetric observations

Authors: Ryo Tazaki, Alexandre Lazarian, Hideko Nomura

Abstract: We apply the theory of radiative torque (RAT) alignment for studying protoplanetary disks around a T-Tauri star and perform 3D radiative transfer calculations to provide the expected maps of polarized radiation to be compared with observations, such as with ALMA. We revisit the issue of grain alignment for large grains expected in the protoplanetary disks and find that mm-sized grains at midplane… ▽ More We apply the theory of radiative torque (RAT) alignment for studying protoplanetary disks around a T-Tauri star and perform 3D radiative transfer calculations to provide the expected maps of polarized radiation to be compared with observations, such as with ALMA. We revisit the issue of grain alignment for large grains expected in the protoplanetary disks and find that mm-sized grains at midplane do not align with the magnetic field as the Larmor precession timescale for such large grains becomes longer than the gaseous damping timescale. Hence, for these grains the RAT theory predicts that the alignment axis is determined by the grain precession with respect to the radiative flux. As a result, we expect that the polarization will be in the azimuthal direction for a face-on disk. It is also shown that if dust grains have superparamagnetic inclusions, magnetic field alignment is possible for (sub-)micron grains at the surface layer of disks, and this can be tested by mid-infrared polarimetric observations. △ Less

Submitted 9 March, 2017; v1 submitted 8 January, 2017; originally announced January 2017.

Comments: 17 pages, 12 figures, 1 table, accepted for publication in ApJ

Light scattering by fractal dust aggregates: I. Angular dependence of scattering

Authors: Ryo Tazaki, Hidekazu Tanaka, Satoshi Okuzumi, Akimasa Kataoka, Hideko Nomura

Abstract: In protoplanetary disks, micron-sized dust grains coagulate to form highly porous dust aggregates. Because the optical properties of these aggregates are not completely understood, it is important to investigate how porous dust aggregates scatter light. In this study, the light scattering properties of porous dust aggregates were calculated using a rigorous method, the T-matrix method, and the res… ▽ More In protoplanetary disks, micron-sized dust grains coagulate to form highly porous dust aggregates. Because the optical properties of these aggregates are not completely understood, it is important to investigate how porous dust aggregates scatter light. In this study, the light scattering properties of porous dust aggregates were calculated using a rigorous method, the T-matrix method, and the results were then compared with those obtained using the Rayleigh-Gans-Debye (RGD) theory and Mie theory with the effective medium approximation (EMT). The RGD theory is applicable to moderately large aggregates made of nearly transparent monomers. This study considered two types of porous dust aggregates, ballistic cluster-cluster agglomerates (BCCAs) and ballistic particle-cluster agglomerates (BPCAs). First, the angular dependence of the scattered intensity was shown to reflect the hierarchical structure of dust aggregates; the large-scale structure of the aggregates is responsible for the intensity at small scattering angles, and their small-scale structure determines the intensity at large scattering angles. Second, it was determined that the EMT underestimates the backward scattering intensity by multiple orders of magnitude, especially in BCCAs, because the EMT averages the structure within the size of the aggregates. It was concluded that the RGD theory is a very useful method for calculating the optical properties of BCCAs. △ Less

Submitted 24 March, 2016; originally announced March 2016.

Comments: 16 pages, 15 figures, 1 table, accepted for publication in ApJ

Outward Motion of Porous Dust Aggregates by Stellar Radiation Pressure in Protoplanetary Disks

Authors: Ryo Tazaki, Hideko Nomura

Abstract: We study the dust motion at the surface layer of protoplanetary disks. Dust grains in surface layer migrate outward due to angular momentum transport via gas-drag force induced by the stellar radiation pressure. In this study, we calculate mass flux of the outward motion of compact grains and porous dust aggregates by the radiation pressure. The radiation pressure force for porous dust aggregates… ▽ More We study the dust motion at the surface layer of protoplanetary disks. Dust grains in surface layer migrate outward due to angular momentum transport via gas-drag force induced by the stellar radiation pressure. In this study, we calculate mass flux of the outward motion of compact grains and porous dust aggregates by the radiation pressure. The radiation pressure force for porous dust aggregates is calculated using the T-Matrix Method for the Clusters of Spheres. First, we confirm that porous dust aggregates are forced by strong radiation pressure even if they grow to be larger aggregates in contrast to homogeneous and spherical compact grains to which efficiency of radiation pressure becomes lower when their sizes increase. In addition, we find that the outward mass flux of porous dust aggregates with monomer size of 0.1 $μ$m is larger than that of compact grains by an order of magnitude at the disk radius of 1 AU, when their sizes are several microns. This implies that large compact grains like calcium-aluminum rich inclusions (CAIs) are hardly transported to outer region by stellar radiation pressure, whereas porous dust aggregates like chondritic-porous interplanetary dust particles (CP-IDPs) are efficiently transported to comet formation region. Crystalline silicates are possibly transported in porous dust aggregates by stellar radiation pressure from inner hot region to outer cold cometary region in the protosolar nebula. △ Less

Submitted 18 November, 2014; originally announced November 2014.

Comments: 9 pages, 5 figures, 1 table, accepted for publication in ApJ