BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Satoshi Okuzumi

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

ResearcherUnclaimed source record
astro-ph.EPastro-ph.SRastro-ph.GAphysics.plasm-ph

Catalog footprint

What is connected

45works
4topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

45 published item(s)

preprint2026arXiv

A Dichotomy of the Mass-Metallicity Relation of Exoplanetary Atmospheres Demarcated by their Birthplace

Atmospheric observations by JWST raise a growing evidence that atmospheric metallicity exhibits an anti-correlation with masses of giant exoplanets. While such a trend was anticipated by planetesimal-based planet formation models, it remains unclear what kind of atmospheric metallicity trends emerge from pebble-based planet formation. Moreover, while recent studies of solar system Jupiter suggest that uppermost observable atmosphere may not represent the bulk envelope composition, it remains uncertain how the envelope inhomogeneity influences the atmospheric metallicity trend. In this study, we develop disk evolution and planet formation models to investigate the possible atmospheric metallicity trends of giant exoplanets formed via pebble accretion and how they depend on the metallicity inhomogeneity within the envelope. We find that pebble-based planet formation produces two distinct mass-metallicity relations depending on planetary birthplace. Planets formed beyond the H2O snowline exhibit a mass-metallicity anti-correlation similar to that predicted by planetesimal-based models if their atmospheres are fully convective. This anti-correlation disappears if the convective mixing is inefficient. In contrast, planets formed inside the H2O snowline show a shallower mass-metallicity anti-correlation, regardless of the efficiency of atmospheric mixing. Many gas giants observed by JWST observations lie around the mass-metallicity relation predicted for formation at close-in orbits, although some planets with sub-stellar atmospheric metallicity appear to require unmixed envelopes and formation beyond the H2O snowline. We also examine the relationship between bulk and atmospheric metallicity and find a clear correlation that closely follows atmospheric metallicity that is comparable to bulk metallicity.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Structured Distributions of Gas and Solids in Protoplanetary Disks

Recent spatially-resolved observations of protoplanetary disks revealed a plethora of substructures, including concentric rings and gaps, inner cavities, misalignments, spiral arms, and azimuthal asymmetries. This is the major breakthrough in studies of protoplanetary disks since Protostars and Planets VI and is reshaping the field of planet formation. However, while the capability of imaging substructures in protoplanetary disks has been steadily improving, the origin of many substructures are still largely debated. The structured distributions of gas and solids in protoplanetary disks likely reflect the outcome of physical processes at work, including the formation of planets. Yet, the diverse properties among the observed protoplanetary disk population, for example, the number and radial location of rings and gaps in the dust distribution, suggest that the controlling process may differ between disks and/or the outcome may be sensitive to stellar or disk properties. In this review, we (1) summarize the existing observations of protoplanetary disk substructures collected from the literature; (2) provide a comprehensive theoretical review of various processes proposed to explain observed protoplanetary disk substructures; (3) compare current theoretical predictions with existing observations and highlight future research directions to distinguish between different origins; and (4) discuss implications of state-of-the-art protoplanetary disk observations to protoplanetary disk and planet formation theory.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A global two-layer radiative transfer model for axisymmetric, shadowed protoplanetary disks

Understanding the thermal structure of protoplanetary disks is crucial for modeling planet formation and interpreting disk observations. We present a new two-layer radiative transfer model for computing the thermal structure of axisymmetric irradiated disks. Unlike the standard two-layer model, our model accounts for the radial as well as vertical transfer of the starlight reprocessed at the disk surface. The model thus allows us to compute the temperature below "shadowed" surfaces receiving no direct starlight. Thanks to the assumed axisymmetry, the reprocessed starlight flux is given in one-dimensional integral form that can be computed at a low cost. Furthermore, our model evolves the midplane temperature using a time-dependent energy equation and can therefore treat thermal instabilities. We apply our global two-layer model to disks with a planetary induced gap and confirm that the model reproduces the disks' temperature profiles obtained from more computationally expensive Monte Carlo radiative transfer calculations to an accuracy of less than 20%. We also apply the model to study the long-term behavior of the thermal wave instability in irradiated disks. Being simple and computationally efficient, the global two-layer model will be suitable for studying the interplay between disks' thermal evolution and dust evolution.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Color Dependence of the Transit Detectability for Young Active M-dwarfs

We investigate the planetary transit detectability in the presence of stellar rotational activity from light curves for young M-dwarfs and estimate improvements of the detection at near-infrared (NIR) wavelengths. Making maps of the transit signal detection efficiency over the orbital period and planetary radius with light curves of members of four clusters, Hyades, Praesepe, Pleiades, and Upper Scorpius observed by the K2 mission, we evaluate the detectability for the rotation period and modulation semi-amplitude. We find that the detection efficiency remarkably decreases to about 20% for rapidly rotators with P_{rot} <= 1 d and the lack of planets in Pleiades is likely due to the high fraction of rapidly rotating M-dwarfs. We also evaluate the improvements of the planet detection with NIR photometry via tests using mock light curves assuming that the signal amplitude of stellar rotation decreases at NIR wavelengths. Our results suggest that NIR photometric monitoring would double relative detection efficiency for transiting planetary candidates with P_{rot} <= 1 d and find planets around M-dwarfs with approximately 100 Myr missing in the past transit surveys from the space.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Impact of dust size distribution including large dust grains on magnetic resistivity: an analytical approach

This paper investigates the impact of dust size distribution on magnetic resistivity. In particular, we focus on its impact when the maximum dust size significantly increases from sub-micron. The first half of the paper describes our calculation method for magnetic resistivity based on the model of \citet{1987ApJ...320..803D} and shows that the method reproduces the results of a more realistic chemical reaction network calculations reasonably well. Then, we describe the results of the resistivity calculations for dust distributions with large maximum dust grains. Our results show that resistivity tends to decrease with dust growth, which is particularly true when the dust size power exponent $q$ is $q=2.5$. On the other hand, the decrease is less pronounced when the dust size power exponent $q$ is $q=3.5$, i.e., when the small dust is also responsible for the dust cross-section. Our results suggest that detailed dust coagulation and fragmentation processes play a vital role in the magnetic resistivities in protostar formation.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Planetesimal formation around the snow line. II. Dust or pebbles?

Around the snow line, icy pebbles and silicate dust may locally pile-up and form icy and rocky planetesimals via streaming instability and/or gravitational instability. We perform 1D diffusion-advection simulations that include the back-reaction to radial drift and diffusion of icy pebbles and silicate dust, ice sublimation, release of silicate dust, and their recycling through recondensation and sticking onto pebbles outside the snow line. We use a realistic description of the scale height of silicate dust obtained from Ida et al. and that of pebbles including the effects of a Kelvin-Helmholtz instability. We study the dependence of solid pile-up on distinct effective viscous parameters for turbulent diffusions in the radial and vertical directions ($α_{\rm Dr}$ and $α_{\rm Dz}$) and for the gas accretion to the star ($α_{\rm acc}$) as well as that on the pebble-to-gas mass flux ($F_{\rm p/g}$). We derive the sublimation width of drifting icy pebbles which is a critical parameter to characterize the pile-up of silicate dust and pebbles around the snow line. We identify a parameter space (in the $F_{\rm p/g}-α_{\rm acc}-α_{\rm Dz}(=α_{\rm Dr})$ space) where pebbles no longer drift inward to reach the snow line due to the back-reaction that slows down radial velocity of pebbles. We show that the pile-up of solids around the snow line occurs in a broader range of parameters for $α_{\rm acc}=10^{-3}$ than for $α_{\rm acc}=10^{-2}$. Above a critical $F_{\rm p/g}$ value, the runaway pile-up of silicate dust inside the snow line is favored for $α_{\rm Dr}/α_{\rm acc} \ll 1$, while that of pebbles outside the snow line is favored for $α_{\rm Dr}/α_{\rm acc} \sim 1$. Our results imply that a distinct evolutionary path could produce a diversity of outcomes in terms of planetesimal formation around the snow line.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Steady-state accretion in magnetized protoplanetary disks

[abridged] We present a 1+1D global magnetically-driven disk accretion model that captures the essence of the MRI-driven accretion, without resorting to 3D global non-ideal MHD simulations. The gas dynamics is assumed to be solely controlled by the MRI and hydrodynamic instabilities. For given stellar and disk parameters, the Shakura-Sunyaev viscosity parameter $α$ is determined self-consistently under the framework of viscously-driven accretion from detailed considerations of the MRI with non-ideal MHD effects (Ohmic resistivity and ambipolar diffusion), accounting for disk heating by stellar irradiation, non-thermal sources of ionization, and dust effects on the ionization chemistry. Additionally, the magnetic field strength is constrained and adopted to maximize the MRI activity. We demonstrate the use of our framework by investigating steady-state MRI-driven accretion in a fiducial protoplanetary disk model around a solar-type star. We find that the equilibrium solution displays no pressure maximum at the dead zone outer edge, except if a sufficient amount of dust particles have accumulated there before the disk reaches a steady-state accretion regime. Furthermore, the steady-state accretion solution describes a disk that displays a spatially extended long-lived inner disk gas reservoir (the dead zone) accreting a few $10^{-9}\, M_{\odot}.\rm{yr}^{-1}$. By conducting a detailed parameter study, we find that the extend to which the MRI can drive efficient accretion is primarily determined by the total disk gas mass, the representative grain size, the vertically-integrated dust-to-gas mass ratio, and the stellar X-ray luminosity. A self-consistent time-dependent coupling between gas, dust, stellar evolution models and our general framework on million-year timescales is required to fully understand the formation of dead zones and their potential to trap dust particles.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Planetesimal formation around the snow line: I. Monte Carlo simulations of silicate dust pile-up in a turbulent disk

Context: The formation of rocky planetesimals is a long-standing problem in planet formation theory. One of the possibilities is that it results from gravitational instability as a result of pile-up of small silicate dust particles released from sublimating icy pebbles that pass the snow line. Aims: We want to understand and quantify the role of the water snow line for the formation of rock-rich and ice-rich planetesimals. In this paper, we focus on the formation of rock-rich planetesimals. A companion paper examines the combined formation of both rock-rich and ice-rich planetesimals. Methods: We develop a new Monte Carlo code to calculate the radial evolution of silicate particles in a turbulent accretion disk, accounting for the back-reaction (i.e., inertia) of the particles on their radial drift velocity and diffusion. Results depend in particular on the particle injection width (determined from the radial sublimation width of icy pebbles), the pebble scale height and the pebble mass flux through the disk. The scale height evolution of the silicate particles, which is the most important factor for the runaway pile-up, is automatically calculated in this Lagrange method. Results: From the numerical results, we derive semi-analytical relations for the scale height of the silicate dust particles and the particles-to-gas density ratio at the midplane, as functions of a pebble-to-gas mass flux ratio and the $α$ parameters for disk gas accretion and vertical/radial diffusion. We find that the runaway pile-up of the silicate particles (formation of rocky planetesimals) occurs if the pebble-to-gas mass flux ratio is $> [(α_{Dz}/α_{acc})/3 \times 10^{-2}]^{1/2}$ where $α_{Dz}$ and $α_{acc}$ are the $α$ parameters for vertical turbulent diffusion and disk gas accretion.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The Galilean Satellites Formed Slowly from Pebbles

It is generally accepted that the four major (Galilean) satellites formed out of the gas disk that accompanied Jupiter's formation. However, understanding the specifics of the formation process is challenging as both small particles (pebbles) as well as the satellites are subject to fast migration processes. Here, we hypothesize a new scenario for the origin of the Galilean system, based on the capture of several planetesimal seeds and subsequent slow accretion of pebbles. To halt migration, we invoke an inner disk truncation radius, and other parameters are tuned for the model to match physical, dynamical, compositional, and structural constraints. In our scenario it is natural that Ganymede's mass is determined by pebble isolation. Our slow-pebble-accretion scenario then reproduces the following characteristics: (1) the mass of all the Galilean satellites; (2) the orbits of Io, Europa, and Ganymede captured in mutual 2:1 mean motion resonances; (3) the ice mass fractions of all the Galilean satellites; (4) the unique ice-rock partially differentiated Callisto and the complete differentiation of the other satellites. Our scenario is unique to simultaneously reproduce these disparate properties.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the water delivery to terrestrial embryos by ice pebble accretion

Standard accretion disk models suggest that the snow line in the solar nebula migrated interior to the Earth's orbit in a late stage of nebula evolution. In this late stage, a significant amount of ice could have been delivered to 1 AU from outer regions in the form of mm to dm-sized pebbles. This raises the question why the present Earth is so depleted of water (with the ocean mass being as small as 0.023% of the Earth mass). Here we quantify the amount of icy pebbles accreted by terrestrial embryos after the migration of the snow line assuming that no mechanism halts the pebble flow in outer disk regions. We use a simplified version of the coagulation equation to calculate the formation and radial inward drift of icy pebbles in a protoplanetary disk. The pebble accretion cross section of an embryo is calculated using analytic expressions presented by recent studies. We find that the final mass and water content of terrestrial embryos strongly depends on the radial extent of the gas disk, the strength of disk turbulence, and the time at which the snow lines arrives at 1 AU. The disk's radial extent sets the lifetime of the pebble flow, while turbulence determines the density of pebbles at the midplane where the embryos reside. We find that the final water content of the embryos falls below 0.023 wt% only if the disk is compact (< 100 AU), turbulence is strong at 1 AU, and the snow line arrives at 1 AU later than 2-4 Myr after disk formation. If the solar nebula extended to 300 AU, initially rocky embryos would have evolved into icy planets of 1-10 Earth masses unless the snow-line migration was slow. If the proto-Earth contained water of ~ 1 wt% as might be suggested by the density deficit of the Earth's outer core, the formation of the proto-Earth was possible with weaker turbulence and with earlier (> 0.5-2 Myr) snow-line migration.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Development of a Method for the Observation of Lightning in Protoplanetary Disks Using Ion Lines

In this paper, we propose observational methods for detecting lightning in protoplanetary disks. We do so by calculating the critical electric field strength in the lightning matrix gas (LMG), the parts of the disk where the electric field is strong enough to cause lightning. That electric field accelerates multiple positive ion species to characteristic terminal velocities. In this paper, we present three distinct discharge models, with corresponding critical electric fields. We simulate the position-velocity diagrams and the integrated emission maps for the models. We calculate the measure of sensitivity values for detection of the models, and for distinguishing between the models. At the distance of TW-Hya (54pc), LMG that occupies $2π$ in azimuth and $25 \mathrm{au}<r<50 \mathrm{au}$ is $1200σ$- to $4000σ$-detectable. The lower limits of the radii of $5σ$-detectable LMG clumps are between 1.6 au and 5.3 au, depending on the models.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Effects of Ohmic and ambipolar diffusion on the formation and evolution of the first cores, protostars and circumstellar discs

We investigate the formation and evolution of a first core, protostar, and circumstellar disc with a three-dimensional non-ideal (including both Ohmic and ambipolar diffusion) radiation magnetohydrodynamics simulation. We found that the magnetic flux is largely removed by magnetic diffusion in the first core phase and that the plasma $β$ of the centre of the first core becomes large, $β>10^4$. Thus, proper treatment of first core phase is crucial in investigating the formation of protostar and disc. On the other hand, in an ideal simulation, $β\sim 10$ at the centre of the first core. The simulations with magnetic diffusion show that the circumstellar disc forms at almost the same time of protostar formation even with a relatively strong initial magnetic field (the value for the initial mass-to-flux ratio of the cloud core relative to the critical value is $μ=4$). The disc has a radius of $r \sim 1$ AU at the protostar formation epoch. We confirm that the disc is rotationally supported. We also show that the disc is massive ($Q\sim 1$) and that gravitational instability may play an important role in the subsequent disc evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Electron Heating in the Magnetorotational Instability: Implications for Turbulence Strength in Outer Regions of Protoplanetary Disks

The magnetorotational instability (MRI) drives vigorous turbulence in a region of protoplanetary disks where the ionization fraction is sufficiently high. It has recently been shown that the electric field induced by the MRI can heat up electrons and thereby affect the ionization balance in the gas. In particular, in a disk where abundant dust grains are present, the electron heating causes a reduction of the electron abundance, thereby preventing further growth of the MRI. By using the nonlinear Ohm's law that takes into account electron heating, we investigate where in protoplanetary disks this negative feedback between the MRI and ionization chemistry becomes important. We find that the "e-heating zone," the region where the electron heating limits the saturation of the MRI, extends out up to 80 AU in the minimum-mass solar nebula with abundant submicron-sized grains. This region is considerably larger than the conventional dead zone whose radial extent is $\sim20$ AU in the same disk model. Scaling arguments show that the MRI turbulence in the e-heating zone should have a significantly lower saturation level. Submicron-sized grains in the e-heating zone are so negatively charged that their collisional growth is unlikely to occur. Our present model neglects ambipolar and Hall diffusion, but our estimate shows that ambipolar diffusion would also affect the MRI in the e-heating zone.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On shocks driven by high-mass planets in radiatively inefficient disks. II. Three-dimensional global disk simulations

Recent high-resolution, near-infrared images of protoplanetary disks have shown that these disks often present spiral features. Spiral arms are among the structures predicted decades ago by numerical simulations of disk-planet interaction and thus it is tempting to suspect that planetary perturbers are responsible for the observed signatures. However, such interpretation is not free of problems. The spirals are found to have large pitch angles, and in at least one case (HD 100546) the spiral feature appears effectively unpolarized, which implies thermal emission of the order of 1000K (465$\pm$40K at closer inspection). We have recently shown in two-dimensional models that shock dissipation in the supersonic wake of high-mass planets can lead to significant heating if the disk is sufficiently adiabatic. In this paper we extend this analysis to three dimensions in thermodynamically evolving disks. We use the Pencil Code in spherical coordinates for our models, with a prescription for thermal cooling based on the optical depth of the local vertical gas column. We use a 5$M_J$ planet, and show that shocks in the region around the planet where the Lindblad resonances occur heat the gas to substantially higher temperatures than the ambient disk gas at that radius. The gas is accelerated vertically away from the midplane by the shocks to form shock bores, and the gas falling back toward the midplane breaks up into a turbulent surf near the Lindblad resonances. This turbulence, although localized, has high $α$ values, reaching 0.05 in the inner Lindblad resonance, and 0.1 in the outer one. We also find evidence that the disk regions heated up by the planetary shocks eventually becomes superadiabatic, generating convection far from the planet's orbit.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Radiation Magnetohydrodynamic Simulations of Protostellar Collapse: Non-Ideal Magnetohydrodynamic Effects and Early Formation of Circumstellar Disks

The transport of angular momentum by magnetic fields is a crucial physical process in formation and evolution of stars and disks. Because the ionization degree in star forming clouds is extremely low, non-ideal magnetohydrodynamic (MHD) effects such as ambipolar diffusion and Ohmic dissipation work strongly during protostellar collapse. These effects have significant impacts in the early phase of star formation as they redistribute magnetic flux and suppress angular momentum transport by magnetic fields. We perform three-dimensional nested-grid radiation magnetohydrodynamic (RMHD) simulations including Ohmic dissipation and ambipolar diffusion. Without these effects, magnetic fields transport angular momentum so efficiently that no rotationally supported disk is formed even after the second collapse. Ohmic dissipation works only in a relatively high density region within the first core and suppresses angular momentum transport, enabling formation of a very small rotationally supported disk after the second collapse. With both Ohmic dissipation and ambipolar diffusion, these effects work effectively in almost the entire region within the first core and significant magnetic flux loss occurs. As a result, a rotationally supported disk is formed even before a protostellar core forms. The size of the disk is still small, about 5 AU at the end of the first core phase, but this disk will grow later as gas accretion continues. Thus the non-ideal MHD effects can resolve the so-called magnetic braking catastrophe while maintaining the disk size small in the early phase, which is implied from recent interferometric observations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the Viability of the Magnetorotational Instability in Circumplanetary Disks

We examine whether the magnetorotational instability (MRI) can serve as a mechanism of angular momentum transport in circumplanetary disks. For the MRI to operate the ionization degree must be sufficiently high and the magnetic pressure must be sufficiently lower than the gas pressure. We calculate the spatial distribution of the ionization degree and search for the MRI-active region where the two criteria are met. We find that there can be thin active layers at the disk surface depending on the model parameters, however, we find hardly any region which can sustain well-developed MRI turbulence; when the magnetic field is enhanced by MRI turbulence at the disk surface layer, a magnetically dominated atmosphere encroaches on a lower altitude and a region of well-developed MRI turbulence becomes smaller. We conclude that if there are no angular momentum transfer mechanisms other than MRI in gravitationally stable circumplanetary disks, gas is likely to pile up until disks become gravitationally unstable, and massive disks may survive for a long time.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Opacity of fluffy dust aggregates

Context. Dust grains coagulate to form dust aggregates in protoplanetary disks. Their porosity can be extremely high in the disks. Although disk emission may come from fluffy dust aggregates, the emission has been modeled with compact grains. Aims. We aim to reveal the mass opacity of fluffy aggregates from infrared to millimeter wavelengths with the filling factor ranging from 1 down to $10^{-4}$. Methods. We use Mie calculations with an effective medium theory. The monomers are assumed to be 0.1 ${\rm μm}$ sized grains, which is much shorter than the wavelengths that we focus on. Results. We find that the absorption mass opacity of fluffy aggregates are characterized by the product $a\times f$, where $a$ is the dust radius and $f$ is the filling factor, except for the interference structure. The scattering mass opacity is also characterized by $af$ at short wavelengths while it is higher in more fluffy aggregates at long wavelengths. We also derive the analytic formula of the mass opacity and find that it reproduces the Mie calculations. We also calculate the expected difference of the emission between compact and fluffy aggregates in protoplanetary disks with a simple dust growth and drift model. We find that compact grains and fluffy aggregates can be distinguished by the radial distribution of the opacity index $β$. The previous observation of the radial distribution of $β$ is consistent with the fluffy case, but more observations are required to distinguish between fluffy or compact. In addition, we find that the scattered light would be another way to distinguish between compact grains and fluffy aggregates.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Nonlinear Ohm's Law: Plasma Heating by Strong Electric Fields and its Effects on the Ionization Balance in Protoplanetary Disks

The ionization state of the gas plays a key role in the MHD of protoplanetary disks. However, the ionization state can depend on the gas dynamics, because electric fields induced by MHD turbulence can heat up plasmas and thereby affect the ionization balance. To study this nonlinear feedback, we construct an ionization model that includes plasma heating by electric fields and impact ionization by heated electrons, as well as charging of dust grains. We show that when plasma sticking onto grains is the dominant recombination process, the electron abundance in the gas decreases with increasing electric field strength. This is a natural consequence of electron-grain collisions whose frequency increases with electron's random velocity. The decreasing electron abundance may lead to a self-regulation of MHD turbulence. In some cases, not only the electron abundance but also the electric current decreases with increasing field strength in a certain field range. The resulting N-shaped current--field relation violates the fundamental assumption of the non-relativistic MHD that the electric field is uniquely determined by the current density. At even higher field strengths, impact ionization causes an abrupt increase of the electric current as expected by previous studies. We find that this discharge current is multi-valued (i.e., the current--field relation is S-shaped) under some circumstances, and that the intermediate branch is unstable. The N/S-shaped current--field relations may yield hysteresis in the evolution of MHD turbulence in some parts of protoplanetary disks.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Fluffy dust forms icy planetesimals by static compression

Context: In planetesimal formation theory, several barriers have been proposed, which are bouncing, fragmentation, and radial drift problems. To understand the structure evolution of dust aggregates is a key in the planetesimal formation. Dust grains become fluffy by coagulation in protoplanetary disks. However, once they become fluffy, they are not sufficiently compressed by collisional compression to form compact planetesimals. Aims: We aim to reveal the pathway of the dust structure evolution from dust grains to compact planetesimals. Methods: Using the compressive strength formula, we analytically investigate how fluffy dust aggregates are compressed by static compression due to ram pressure of the disk gas and self gravity of the aggregates in protoplanetary disks. Results: We reveal the pathway of the porosity evolution from dust grains via fluffy aggregates to form planetesimals, circumventing the barriers in planetesimal formation. The aggregates are compressed by the disk gas to the density of 10^{-3} g/cm^3 in coagulation, which is more compact than the case with collisional compression. Then, they are compressed more by self gravity to 10^{-1} g/cm^3 when the radius is 10 km. Although the gas compression decelerate the growth, they grow enough rapidly to avoid the radial drift barrier when the orbital radius is < 6 AU in a typical disk. Conclusions: We propose fluffy dust growth scenario from grains to planetesimals. It enables the icy planetesimal formation in a wide range beyond the snowline in protoplanetary disks. This result proposes a concrete initial condition of planetesimals for the later stages of the planet formation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Interdependence of Electric Discharge and The Magnetorotational Instability in Protoplanetary Disks

We study how the magnetorotational instability (MRI) in protoplanetary disks is affected by the electric discharge caused by the electric field in the resistive MHD. We have performed three-dimensional shearing box simulations with various values of plasma beta and electrical breakdown models. We find the self-sustainment of the MRI in spite of the high resistivity. The instability gives rise to the large electric field that causes the electrical breakdown, and the breakdown maintains the high ionization degree required for the instability. The condition for this self-sustained MRI is set by the balance between the energy supply from the shearing motion and the energy consumed by the Ohmic dissipation. We apply the condition to various disk models and study where the active, self-sustained, and dead zones of MRI are. In the fiducial minimum-mass solar nebula (MMSN) model, the newly-found sustained zone occupies only the limited volume of the disk. In the late-phase gas-depleted disk models, however, the sustained zone occupies larger volume of the disk.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Static compression of porous dust aggregates

Context: In protoplanetary disks, dust grains coagulate with each other and grow to form aggregates. As these aggregates grow by coagulation, their filling factor ϕdecreases down to ϕ<< 1. However, comets, the remnants of these early planetesimals, have ϕ~ 0.1. Thus, static compression of porous dust aggregates is important in planetesimal formation. However, the static compression strength has been investigated only for relatively high density aggregates (ϕ> 0.1). Aims: We investigate and find the compression strength of highly porous aggregates (ϕ<< 1). Methods: We perform three dimensional N-body simulations of aggregate compression with a particle-particle interaction model. We introduce a new method of static compression: the periodic boundary condition is adopted and the boundaries move with low speed to get closer. The dust aggregate is compressed uniformly and isotropically by themselves over the periodic boundaries. Results: We empirically derive a formula of the compression strength of highly porous aggregates (ϕ<< 1). We check the validity of the compression strength formula for wide ranges of numerical parameters, such as the size of initial aggregates, the boundary speed, the normal damping force, and material. We also compare our results to the previous studies of static compression in the relatively high density region (ϕ> 0.1) and confirm that our results consistently connect to those in the high density region. The compression strength formula is also derived analytically.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Fate of Planetesimals in Turbulent Disks with Dead Zones. I. The Turbulent Stirring Recipe

Turbulence in protoplanetary disks affects planet formation in many ways. While small dust particles are mainly affected by the aerodynamical coupling with turbulent gas velocity fields, planetesimals and larger bodies are more affected by gravitational interaction with gas density fluctuations. For the latter process, a number of numerical simulations have been performed in recent years, but a fully parameter-independent understanding has not been yet established. In this study, we present simple scaling relations for the planetesimal stirring rate in turbulence driven by magnetorotational instability (MRI), taking into account the stabilization of MRI due to Ohmic resistivity. We begin with order-of-magnitude estimates of the turbulence-induced gravitational force acting on solid bodies and associated diffusion coefficients for their orbital elements. We then test the predicted scaling relations using the results of recent Ohmic-resistive MHD simulations by Gressel et al. We find that these relations successfully explain the simulation results if we properly fix order-of-unity uncertainties within the estimates. We also update the saturation predictor for the density fluctuation amplitude in MRI-driven turbulence originally proposed by Okuzumi & Hirose. Combination of the scaling relations and saturation predictor allows to know how the turbulent stirring rate of planetesimals depends on disk parameters such as the gas column density, distance from the central star, vertical resistivity distribution, and net vertical magnetic flux. In Paper II, we apply our recipe to planetesimal accretion to discuss its viability in turbulent disks.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The fate of planetesimals in turbulent disks with dead zones. II. Limits on the viability of runaway accretion

A critical phase in the standard model for planet formation is the runaway growth phase. During runaway growth bodies in the 0.1--100 km size range (planetesimals) quickly produce a number of much larger seeds. The runaway growth phase is essential for planet formation as the emergent planetary embryos can accrete the leftover planetesimals at large gravitational focusing factors. However, torques resulting from turbulence-induced density fluctuations may violate the criterion for the onset of runaway growth, which is that the magnitude of the planetesimals' random (eccentric) motions are less than their escape velocity. This condition represents a more stringent constraint than the condition that planetesimals survive their mutual collisions. To investigate the effects of MRI turbulence on the viability of the runaway growth scenario, we apply our semi-analytical recipes of Paper I, which we augment by a coagulation/fragmentation model for the dust component. We find that the surface area-equivalent abundance of 0.1 micron particles is reduced by factors 10^2--10^3, which tends to render the dust irrelevant to the turbulence. We express the turbulent activity in the midplane regions in terms of a size s_run above which planetesimals will experience runaway growth. We find that s_run is mainly determined by the strength of the vertical net field that threads the disks and the disk radius. At disk radii beyond 5 AU, s_run becomes larger than ~100 km and the collision times among these bodies longer than the duration of the nebula phase. Our findings imply that the classical, planetesimal-dominated, model for planet formation is not viable in the outer regions of a turbulent disk.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Geometrical Cross Sections of Dust Aggregates and a Compression Model for Aggregate Collisions

Geometrical cross sections of dust aggregates determine their coupling with disk gas, which governs their motions in protoplanetary disks. Collisional outcomes also depend on geometrical cross sections of initial aggregates. In the previous paper, we performed three-dimensional N-body simulations of sequential collisions of aggregates composed of a number of sub-micron-sized icy particles and examined radii of gyration (and bulk densities) of the obtained aggregates. We showed that collisional compression of aggregates is not efficient and that aggregates remain fluffy. In the present study, we examine geometrical cross sections of the aggregates. Their cross sections decreases due to the compression as well as their gyration radii. It is found that a relation between the cross section and the gyration radius proposed by Okuzumi et al. is valid for the compressed aggregates. We also refine the compression model proposed in our previous paper. The refined model enables us to calculate the evolution of both gyration radii and cross sections of growing aggregates and reproduces well our numerical results of sequential aggregate collisions. The refined model can describe non-equal-mass collisions as well as equal-mass case. Although we do not take into account oblique collisions in the present study, oblique collisions would further hinder compression of aggregates.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Planetesimal Formation in Magnetorotationally Dead Zones: Critical Dependence on the Net Vertical Magnetic Flux

Turbulence driven by magnetorotational instability (MRI) affects planetesimal formation by inducing diffusion and collisional fragmentation of dust particles. We examine conditions preferred for planetesimal formation in MRI-inactive "dead zones" using an analytic dead-zone model based on our recent resistive MHD simulations. We argue that successful planetesimal formation requires not only a sufficiently large dead zone (which can be produced by tiny dust grains) but also a sufficiently small net vertical magnetic flux (NVF). Although often ignored, the latter condition is indeed important since the NVF strength determines the saturation level of turbulence in MRI-active layers. We show that direct collisional formation of icy planetesimal across the fragmentation barrier is possible when the NVF strength is lower than 10 mG (for the minimum-mass solar nebula model). Formation of rocky planetesimals via the secular gravitational instability is also possible within a similar range of the NVF strength. Our results indicate that the fate of planet formation largely depends on how the NVF is radially transported in the initial disk formation and subsequent disk accretion processes.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Radiation Magnetohydrodynamic Simulations of Protostellar Collapse: Protostellar Core Formation

We report the first three-dimensional radiation magnetohydrodynamic (RMHD) simulations of protostellar collapse with and without Ohmic dissipation.We take into account many physical processes required to study star formation processes, including a realistic equation of state. We follow the evolution from molecular cloud cores until protostellar cores are formed with sufficiently high resolutions without introducing a sink particle. The physical processes involved in the simulations and adopted numerical methods are described in detail.We can calculate only about one year after the formation of the protostellar cores with our direct three-dimensional RMHD simulations because of the extremely short timescale in the deep interior of the formed protostellar cores, but successfully describe the early phase of star formation processes. The thermal evolution and the structure of the first and second (protostellar) cores are consistent with previous one-dimensional simulations using full radiation transfer, but differ considerably from preceding multi-dimensional studies with the barotropic approximation. The protostellar cores evolve virtually spherically symmetric in the ideal MHD models because of efficient angular momentum transport by magnetic fields, but Ohmic dissipation enables the formation of the circumstellar disks in the vicinity of the protostellar cores as in previous MHD studies with the barotropic approximation. The formed disks are still small (less than 0.35 AU) because we simulate only the earliest evolution. We also confirm that two different types of outflows are naturally launched by magnetic fields from the first cores and protostellar cores in the resistive MHD models.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Fast and Accurate Calculation Scheme for Ionization Degrees in Protoplanetary and Circumplanetary Disks with Charged Dust Grains

We develop a fast and accurate calculation method for ionization degrees in protoplanetary and circumplanetary disks including dust grains. We apply our method to calculate the ionization degree of circumplanetary disks. It is important to understand the structure and evolution of protoplanetary/circumplanetary disks since they are thought to be the sites of planet/satellite formation. The turbulence that causes gas accretion is supposed to be driven by magnetorotational instability (MRI) that occurs only when the ionization degree is high enough for magnetic field to be coupled to gas. We calculate the ionization degrees in circumplanetary disks to estimate the sizes of MRI-inactive regions. We properly include the effect of dust grains because they efficiently capture charged particles and make ionization degree lower. Inclusion of dust grains complicates the reaction equations and requires expensive computation. In order to accelerate the calculation of ionization reactions, we develop a semianalytic method based on the charge distribution model proposed previously. This method enables us to study the ionization state of disks for a wide range of model parameters. For a previous model of circum-Jovian disk, we find that an MRI-inactive region covers almost all regions even without dust grains. This suggests that the gas accretion rates in circumplanetary disks are much smaller than previously thought.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Modeling Magnetorotational Turbulence in Protoplanetary Disks with Dead Zones

Turbulence driven by magnetorotational instability (MRI) crucially affects the evolution of solid bodies in protoplanetary disks. On the other hand, small dust particles stabilize MRI by capturing ionized gas particles needed for the coupling of the gas and magnetic fields. To provide an empirical basis for modeling the coevolution of dust and MRI, we perform three-dimensional, ohmic-resistive MHD simulations of a vertically stratified shearing box with an MRI-inactive "dead zone" of various sizes and with a net vertical magnetic flux of various strengths. We find that the vertical structure of turbulence is well characterized by the vertical magnetic flux and three critical heights derived from the linear analysis of MRI in a stratified disk. In particular, the turbulent structure depends on the resistivity profile only through the critical heights and is insensitive to the details of the resistivity profile. We discover scaling relations between the amplitudes of various turbulent quantities (velocity dispersion, density fluctuation, vertical diffusion coefficient, and outflow mass flux) and vertically integrated accretion stresses. We also obtain empirical formulae for the integrated accretion stresses as a function of the vertical magnetic flux and the critical heights. These empirical relations allow to predict the vertical turbulent structure of a protoplanetary disk for a given strength of the magnetic flux and a given resistivity profile.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Electric charging of dust aggregates and its effect on dust coagulation in protoplanetary disks

Mutual sticking of dust aggregates is the first step toward planetesimal formation in protoplanetary disks. In spite that the electric charging of dust particles is well recognized in some contexts, it has been largely ignored in the current modeling of dust coagulation. In this study, we present a general analysis of the dust charge state in protoplanetary disks, and then demonstrate how the electric charging could dramatically change the currently accepted scenario of dust coagulation. First, we describe a new semianalytical method to calculate the dust charge state and gas ionization state self-consistently. This method is far more efficient than previous numerical methods, and provides a general and clear description of the charge state of gas-dust mixture. Second, we apply this analysis to early evolutionary stages where the dust has been thought to grow into fractal ($D \sim 2$) aggregates with a quasi-monodisperse (i.e., narrow) size distribution. We find that, for a wide range of model parameters, the fractal growth is strongly inhibited by the electric repulsion between colliding aggregates and eventually "freezes out" on its way to the subsequent growth stage involving collisional compression. Strong disk turbulence would help the aggregates to overcome this growth barrier, but then it would cause catastrophic collisional fragmentation in later growth stages. These facts suggest that the combination of electric repulsion and collisional fragmentation would impose a serious limitation on dust growth in protoplanetary disks. We propose a possible scenario of dust evolution after the freeze-out. Finally, we point out that the fractal growth of dust aggregates tends to maintain a low ionization degree and, as a result, a large magnetorotationally stable region in the disk.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Numerical Modeling of the Coagulation and Porosity Evolution of Dust Aggregates

Porosity evolution of dust aggregates is crucial in understanding dust evolution in protoplanetary disks. In this study, we present useful tools to study the coagulation and porosity evolution of dust aggregates. First, we present a new numerical method for simulating dust coagulation and porosity evolution as an extension of the conventional Smoluchowski equation. This method follows the evolution of the mean porosity for each aggregate mass simultaneously with the evolution of the mass distribution function. This method reproduces the results of previous Monte Carlo simulations with much less computational expense. Second, we propose a new collision model for porous dust aggregates on the basis of our N-body experiments on aggregate collisions. We first obtain empirical data on porosity changes between the classical limits of ballistic cluster-cluster and particle-cluster aggregation. Using the data, we construct a recipe for the porosity change due to general hit-and-stick collisions as well as formulae for the aerodynamical and collisional cross sections. Simple coagulation simulations using the extended Smoluchowski method show that our collision model explains the fractal dimensions of porous aggregates observed in a full N-body simulation and a laboratory experiment. Besides, we discover that aggregates at the high-mass end of the distribution can have a considerably small aerodynamical cross section per unit mass compared with aggregates of lower masses. We point out an important implication of this discovery for dust growth in protoplanetary disks.

0 citations0 reviewsNo code linkedOpen access