Substructures in PPDs, whose ubiquity was unveiled by recent ALMA observations, are widely discussed regarding their possible origins. We carry out global full magnetohydrodynamic (MHD) simulations in axisymmetry, coupled with self-consistent ray-tracing radiative transfer, thermochemistry, and non-ideal MHD diffusivities. The abundance profiles of grains are also calculated based on the global dust evolution calculation, including sintering effects. We found that dust size plays a crucial role in the ring formation around the snow lines of protoplanetary disks (PPDs) through the accretion process. Disk ionization structures and thus tensorial conductivities depend on the size of grains.When grains are significantly larger than PAHs, the non-ideal MHD conductivities change dramatically across each snow line of major volatiles, leading to a sudden change of the accretion process across the snow lines and the subsequent formation of gaseous rings/gaps there. On the other hand,the variations of conductivities are a lot less with only PAH sized grains in disks and then these disks retain smoother radial density profiles across snow lines.

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 {\mu}m}$). The high-metallicity atmosphere may offer the clues to explore the gas accretion processes onto past GJ1214b.

Recent high resolution observations unveil ring structures in circumstellar disks. The origin of these rings has been widely investigated under various theoretical scenarios. In this work we perform global 3D non-ideal MHD simulations including effects from both Ohmic resistivity and ambipolar diffusion (AD) to model the HL Tau disk. The non-ideal MHD diffusion profiles are calculated based on the global dust evolution calculation including sintering effects. Disk ionization structure changes dramatically across the snow line due to the change of dust size distribution close to snow line of major volatiles. We find that accretion is mainly driven by disk wind. Gaps and rings can be quickly produced from different accretion rates across snow line. Furthermore, ambipolar diffusion (AD) leads to highly preferential accretion at midplane, followed by magnetic reconnection. This results a local zone of decretion that drains of mass in the field reconnection area, which leaves a gap and an adjacent ring just outside it. Overall, under the favorable condition, both snow lines and non-ideal MHD effects can lead to gaseous gaps and rings in protoplanetary disks.

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.

The Earth's ocean mass is only 2.3 x 10^{-4} of the whole planet mass. Even including water in the interior, it would be at most 10^{-3}-10^{-2}. Ancient Mars may have had a similar or slightly smaller water fraction. It is important to clarify the water delivery mechanism to rocky planets in habitable zones in exoplanetary systems, as well as that to the Earth and Mars. Here, we consider water delivery to planets by icy pebbles after the snowline inwardly passes the planetary orbits and derive the water mass fraction (f_{water}) of the final planet as a function of disk parameters and discuss the parameters that reproduce f_{water} comparable to that inferred for the Earth and ancient Mars. We calculate the growth of icy pebbles and their radial drift with a 1D model, and accretion of icy pebbles onto planets, by simultaneously solving the snowline migration and the disk dissipation, to evaluate f_{water} of the planets. We find that f_{water} is regulated by the total mass (M_{res}) of icy dust materials preserved in the outer disk regions at the timing (t = t_{snow}) of the snowline passage of the planetary orbit. Because M_{res} decays rapidly after the pebble formation front reaches the disk outer edge (at t = t_{pff}), f_{water} is sensitive to the ratio t_{snow}/t_{pff}, which is determined by the disk parameters. We find t_{snow}/t_{pff} < 10 or > 10 is important. Deriving an analytical formula for f_{water} that reproduces the numerical results, we find that f_{water} of a rocky planet near 1 au is ~ 10^{-4}-10^{-2}, in the disks with initial disk size ~ 30-50 au and the initial disk mass accretion rate ~ (10^{-8}-10^{-7}) M_sun/y. Because these disks may be median or slightly compact/massive disks, the water fraction of rocky planets in habitable zones may be often similar to that of the Earth, if the icy pebble accretion is responsible for the water delivery.

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.

We report the detection of an excess in dust continuum emission at 233~GHz (1.3~mm in wavelength) in the protoplanetary disk around TW~Hya revealed through high-sensitivity observations at $\sim$3~au resolution with the Atacama Large Millimeter/submillimeter Array (ALMA). The sensitivity of the 233~GHz image has been improved by a factor of 3 with regard to that of our previous cycle 3 observations. The overall structure is mostly axisymmetric, and there are apparent gaps at 25 and 41 au as previously reported. The most remarkable new finding is a few au-scale excess emission in the south-west part of the protoplanetary disk. The excess emission is located at 52 au from the disk center and is 1.5 times brighter than the surrounding protoplanetary disk at a significance of 12$\sigma$. We performed a visibility fitting to the extracted emission after subtracting the axisymmetric protoplanetary disk emission and found that the inferred size and the total flux density of the excess emission are 4.4$\times$1.0~au and 250~$\mu$Jy, respectively. The dust mass of the excess emission corresponds to 0.03~$M_\oplus$ if a dust temperature of 18~K is assumed. Since the excess emission can also be marginally identified in the Band 7 image at almost the same position, the feature is unlikely to be a background source. The excess emission can be explained by a dust clump accumulated in a small elongated vortex or a massive circumplanetary disk around a Neptune mass forming-planet.

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~\mu 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.

The magnetohydrodynamics (MHD) of protoplanetary disks are strongly subject to the non-ideal MHD effects arising from the low ionization fraction of the disk gas. A strong electric field induced by gas motions can heat ionized gas particles and can thereby affect the ionization balance in the disks. Our previous studies revealed that in dusty protoplanetary disks, the Ohmic conductivity decreases with increasing electric field strength until the electrical breakdown of the disk gas occurs. In this study, we extend our previous work to more general cases where both electric and magnetic fields affect the motion of plasma particles, allowing us to study the impacts of plasma heating on all non-ideal MHD effects: Ohmic, Hall, and ambipolar diffusion. We find that the upper limit on the electric current we previously derived applies even in the presence of magnetic fields. Although the Hall and ambipolar resistivities can either increase or decrease with electric field strength depending on the abundance of charged dust grains, the Ohmic resistivity always increases with electric field strength. An order-of-magnitude estimate suggests that a large-scale electric current generated by gas motions in the inner part of protoplanetary disks could exceed the upper limit. This implies that MHD motions of the inner disk, such as the motion driven by the Hall-shear instability, could either get suppressed or trigger electrical breakdown (lightning discharge). This may have important implications for gas accretion and chondrule formation in the inner part of protoplanetary disks.

The poor stickiness of silicate dust grains is a major obstacle to the formation of rocky planetesimals. In this study, we examine the possibility that silicate grains with an organic mantle, which we call Organic-Mantled Grains (OMGs), form planetesimals through direct coagulation. Organic mantles are commonly found in interplanetary dust particles, and laboratory experiments show that they are softer than silicates, in particular in warm environments. This, combined with the theory of particle adhesion, implies that OMGs are stickier than bare silicate grains. Because organic mantles can survive up to 400 K, silicate grains inside the water snow line in protoplanetary disks can in principle hold such mantles. We construct a simple grain adhesion model to estimate the threshold collision velocity below which aggregates of OMGs can grow. The model shows that aggregates of 0.1 \micron-sized OMGs can overcome the fragmentation barrier in protoplanetary disks if the mantles are as thick as those in interplanetary dust particles and if the temperature is above $\sim$ 200 K. We use this adhesion model to simulate the global evolution of OMG aggregates in the inner part of a protoplanetary disk, demonstrating that OMG aggregates can indeed grow into planetesimals under favorable conditions. Because organic matter is unstable at excessively high temperatures, rocky planetesimal formation by the direct sticking of OMGs is expected to occur in a disk annulus corresponding to the temperature range $ \sim$ 200--400 K. The organic-rich planetesimals may grow into carbon-poor rocky planetesimals by accreting a large amount of carbon-poor chondrules.

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.

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.

The gas temperature in protoplanetary disks (PPDs) is determined by a combination of irradiation heating and accretion heating, with the latter conventionally attributed to turbulent dissipation. However, recent studies have suggested that the inner disk (a few AU) is largely laminar, with accretion primarily driven by magnetized disk winds, as a result of nonideal magnetohydrodynamic (MHD) effects from weakly ionized gas, suggesting an alternative heating mechanism by Joule dissipation. We perform local stratified MHD simulations including all three nonideal MHD effects (ohmic, Hall, and ambipolar diffusion) and investigate the role of Joule heating and the resulting disk vertical temperature profiles. We find that in the inner disk, as ohmic and ambipolar diffusion strongly suppress electrical current around the midplane, Joule heating primarily occurs at several scale heights above the midplane, making the midplane temperature much lower than that with the conventional viscous heating model. Including the Hall effect, Joule heating is enhanced/reduced when the magnetic fields threading the disks are aligned/anti-aligned with the disk rotation, but it is overall ineffective. Our results further suggest that the midplane temperature in the inner PPDs is almost entirely determined by irradiation heating, unless viscous heating can trigger thermal ionization in the disk innermost region to self-sustain magnetorotational instability turbulence.

The Hall effect plays a significant role in star formation because it induces rotation in the infalling envelope, which in turn affects the formation and evolution of the circumstellar disk. The importance of the Hall effect varies with the Hall coefficient, and this coefficient is determined by the fractional abundances of charged species. These abundance values are primarily based on the size and quantity of dust grains as well as the cosmic ray intensity, which respectively absorb and create charged species. Thus, the Hall coefficient varies with both the properties of dust grains and the cosmic ray rate (or ionization source). In this study, we explore the dependence of the Hall coefficient on the grain size and cosmic ray ionization rate using a simplified chemical network model. Following this, using an analytic model, we estimate the typical size of a circumstellar disk induced solely by the Hall effect. The results show that the disk grows during the main accretion phase to a size of ${\sim}$ 3 - 100 au, with the actual size depending on the parameters. These findings suggest that the Hall effect greatly affects circumstellar disk formation, especially in the case that the dust grains have a typical size of ${\sim}$ 0.025 ${\mu}$m - 0.075 ${\mu}$m.

We perform simulations of the dust and gas disk evolution to investigate the observational features of a dust-pileup at the dead-zone inner edge. We show that the total mass of accumulated dust particles is sensitive to the turbulence strength in the dead zone, $\alpha_{\rm dead}$, because of the combined effect of turbulence-induced particle fragmentation (which suppresses particle radial drift) and turbulent diffusion. For a typical critical fragmentation velocity of silicate dust particles of $1~{\rm m~s^{-1}}$, the stress to pressure ratio $\alpha_{\rm dead}$ needs to be lower than $3 \times 10^{-4}$ for dust trapping to operate. The obtained dust distribution is postprocessed using the radiative transfer code RADMC-3D to simulate infrared scattered-light images of the inner part of protoplanetary disks with a dust pileup. We find that a dust pileup at the dead-zone inner edge, if present, casts a shadow extending out to $\sim 10~{\rm au}$. In the shadowed region the temperature significantly drops, which in some cases yields even multiple water snow lines. We also find that even without a dust pileup at the dead-zone inner edge, the disk surface can become thermally unstable, and the excited waves can naturally produce shadows and ring-like structures in observed images. This mechanism might account for the ring-like structures seen in the scattered-light images of some disks, such as the TW Hya disk.

When a planet forms a deep gap in a protoplanetary disk, dust grains cannot pass through the gap. As a consequence, the density of the dust grains can increase up to the same level of the density of the gas at the outer edge. The feedback on the gas from the drifting dust grains is not negligible, in such a dusty region. We carried out two-dimensional two-fluid (gas and dust) hydrodynamic simulations. We found that when the radial flow of the dust grains across the gap is halted, a broad ring of the dust grains can be formed because of the dust feedback and the diffusion of the dust grains. The minimum mass of the planet to form the broad dust ring is consistent with the pebble-isolation mass, in the parameter range of our simulations. The broad ring of the dust grains is good environment for the formation of the protoplanetary solid core. If the ring is formed in the disk around the sun-like star at $\sim 2\ \mbox{AU}$, a massive solid core ($\sim 50M_{\oplus}$) can be formed within the ring, which may be connected to the formation of Hot Jupiters holding a massive solid core such as HD 149026b. In the disk of the dwarf star, a number of Earth-sized planets can be formed within the dust ring around $\sim 0.5\ \mbox{AU}$, which potentially explain the planet system made of multiple Earth-sized planets around the dwarf star such as TRAPPIST-1.

The effect of misalignment between the magnetic field $\magB$ and the angular momentum $\Jang$ of molecular cloud cores on the angular momentum evolution during the gravitational collapse is investigated by ideal and non-ideal MHD simulations. For the non-ideal effect, we consider the ohmic and ambipolar diffusion. Previous studies that considered the misalignment reported qualitatively contradicting results. Magnetic braking was reported as being either strengthened or weakened by misalignment in different studies. We conducted simulations of cloud-core collapse by varying the stability parameter $\alpha$ (the ratio of the thermal to gravitational energy of the core) with and without including magnetic diffusion The non-ideal MHD simulations show the central angular momentum of the core with $\theta=0^\circ$ ($\Jang \parallel \magB$) being always greater than that with $\theta=90^\circ$ ($\Jang \perp \magB$), independently of $\alpha$, meaning that circumstellar disks form more easily form in a core with $\theta=0^\circ$. The ideal MHD simulations, in contrast, show the the central angular momentum of the core with $\theta=90^\circ$ being greater than with $\theta=0^\circ$ for small $\alpha$, and is smaller for large $\alpha$. Inspection of the angular momentum evolution of the fluid elements reveals three mechanisms contributing to the evolution of the angular momentum: (i) magnetic braking in the isothermal collapse phase, (ii) selective accretion of the rapidly (for $\theta=90^\circ$ ) or slowly (for $\theta=0^\circ$) rotating fluid elements to the central region, and (iii) magnetic braking in the first-core and the disk. The difference between the ideal and non-ideal simulations arises from the different efficiencies of (iii).

The ubiquity of clouds in the atmospheres of exoplanets, especially of super-Earths, is one of the outstanding issues for transmission spectra survey. The understanding about the formation process of clouds in super-Earths is necessary to interpret the observed spectra correctly. In this study, we investigate the vertical distributions of particle size and mass density of mineral clouds in super-Earths using a microphysical model that takes into account the vertical transport and growth of cloud particles in a self-consistent manner. We demonstrate that the vertical profiles of mineral clouds significantly vary with the concentration of cloud condensation nuclei and atmospheric metallicity. We find that the height of the cloud top increases with increasing metallicity as long as the metallicity is lower than a threshold. If the metallicity is larger than the threshold, the cloud-top height no longer increases appreciably with metallicity because coalescence yields larger particles of higher settling velocities. We apply our cloud model to GJ1214 b and GJ436 b for which recent transmission observations suggest the presence of high-altitude opaque clouds. For GJ436 b, we show that KCl particles can ascend high enough to explain the observation. For GJ1214 b, by contrast, the height of KCl clouds predicted from our model is too low to explain its flat transmission spectrum. Clouds made of highly porous KCl particles could explain the observations if the atmosphere is highly metal-rich, and hence the particle microstructure might be a key to interpret the flat spectrum of GJ1214 b.

The discovery of thousands of exoplanets over the last couple of decades has shown that the birth of planets is a very efficient process in nature. Theories invoke a multitude of mechanisms to describe the assembly of planets in the disks around pre-main-sequence stars, but observational constraints have been sparse on account of insufficient sensitivity and resolution. Understanding how planets form and interact with their parental disk is crucial also to illuminate the main characteristics of a large portion of the full population of planets that is inaccessible to current and near-future observations. This White Paper describes some of the main issues for our current understanding of the formation and evolution of planets, and the critical contribution expected in this field by the Next Generation Very Large Array.

Star formation in magnetically subcritical clouds is investigated using a three-dimensional non-ideal magneto-hydrodynamics simulation. Since rapid cloud collapse is suppressed until the magnetic flux is sufficiently removed from the initially magnetically subcritical cloud by ambipolar diffusion, it takes > 5-10t_ff to form a protostar, where t_ff is the freefall timescale of the initial cloud. The angular momentum of the star forming cloud is efficiently transferred to the interstellar medium before the rapid collapse begins, and the collapsing cloud has a very low angular momentum. Unlike the magnetically supercritical case, no large-scale low-velocity outflow appears in such a collapsing cloud due to the short lifetime of the first core. Following protostar formation, a very weak high-velocity jet, which has a small momentum and might disappear at a later time, is driven near the protostar, while the circumstellar disc does not grow during the early mass accretion phase. The results show that the star formation process in magnetically subcritical clouds is qualitatively different from that in magnetically supercritical clouds.

We perform three-dimensional radiation non-ideal magnetohydrodynamics simulations and investigate the impact of the Hall effect on the angular momentum evolution in the collapsing cloud cores in which the magnetic field $\magB$ and angular momentum $\angJ$ are misaligned with each other. We find that the Hall effect notably changes the magnetic torques in the pseudo-disk, and strengthens and weakens the magnetic braking in cores with an acute and obtuse relative angles between $\magB$ and $\angJ$, respectively. This suggests that the bimodal evolution of the disk size may occur in early disk evolutionary phase even if $\magB$ and $\angJ$ are randomly distributed. We show that a counter-rotating envelope form in the upper envelope of the pseudo-disk in cloud cores with obtuse relative angles. We also find that a counter-rotating region forms at the midplane of the pseudo-disk in cloud cores with acute relative angles. The former and latter types of counter-rotating envelopes may be associated with the YSOs with a large ($r\sim100$ AU) and small ($r\lesssim10$ AU) disks, respectively.

Magnetorotational instability (MRI) has a potential to generate the vigorous turbulence in protoplanetary disks, although its turbulence strength and accretion stress remains debatable because of the uncertainty of MRI with low ionization fraction. We focus on the heating of electrons by strong electric fields which amplifies nonideal magnetohydrodynamic effects. The heated electrons frequently collide with and stick to dust grains, which in turn decreases the ionization fraction and is expected to weaken the turbulent motion driven by MRI. In order to quantitatively investigate the nonlinear evolution of MRI including the electron heating, we perform magnetohydrodynamical simulation with the unstratified shearing box. We introduce a simple analytic resistivity model depending on the current density by mimicking resistivity given by the calculation of ionization. Our simulation confirms that the electron heating suppresses magnetic turbulence when the electron heating occurs with low current density. We find a clear correlation between magnetic stress and its current density, which means that the magnetic stress is proportional to the squared current density. When the turbulent motion is completely suppressed, laminar accretion flow is caused by ordered magnetic field. We give an analytical description of the laminar state by using a solution of linear perturbation equations with resistivity. We also propose a formula that successfully predicts the accretion stress in the presence of the electron heating.

Chondrules in primitive meteorites likely formed by recrystallisation of dust aggregates that were flash-heated to nearly complete melting. Chondrules may represent the building blocks of rocky planetesimals and protoplanets in the inner regions of protoplanetary discs, but the source of ubiquitous thermal processing of their dust aggregate precursors remains elusive. Here we demonstrate that escape of positrons released in the decay of the short-lived radionuclide $^{26}$Al leads to a large-scale charging of dense pebble structures, resulting in neutralisation by lightning discharge and flash-heating of dust and pebbles. This charging mechanism is similar to a nuclear battery where a radioactive source charges a capacitor. We show that the nuclear battery effect operates in circumplanetesimal pebble discs. The extremely high pebble densities in such discs are consistent with conditions during chondrule heating inferred from the high abundance of sodium within chondrules. The sedimented mid-plane layer of the protoplanetary disc may also be prone to charging by the emission of positrons, if the mass density of small dust there is at least an order of magnitude above the gas density. Our results imply that the decay energy of $^{26}$Al can be harvested to drive intense lightning activity in protoplanetary discs. The total energy stored in positron emission is comparable to the energy needed to melt all solids in the protoplanetary disc. The efficiency of transferring the positron energy to the electric field nevertheless depends on the relatively unknown distribution and scale-dependence of pebble density gradients in circumplanetesimal pebble discs and in the protoplanetary disc mid-plane layer.

Recently, an Atacama Large Millimeter/submillimeter Array (ALMA) observation of the water snow line in the protoplanetary disk around the FU Orionis star V883 Ori was reported. The radial variation of the spectral index at mm-wavelengths around the snow line was interpreted as being due to a pileup of particles interior to the snow line. However, radial transport of solids in the outer disk operates on timescales much longer than the typical timescale of an FU Ori outburst ($10^{1}$--$10^{2}$ yr). Consequently, a steady-state pileup is unlikely. We argue that it is only necessary to consider water evaporation and re-coagulation of silicates to explain the recent ALMA observation of V883 Ori because these processes are short enough to have had their impact since the outburst. Our model requires the inner disk to have already been optically thick before the outburst, and our results suggest that the carbon content of pebbles is low.

The icy satellites around Jupiter are considered to have formed in a circumplanetary disk. While previous models focused on the formation of satellites starting from satellitesimals, the question of how satellitesimals form from smaller dust particles has not been addressed so far. In this work, we study the possibility that satellitesimals form in situ in a circumplanetary disk. We calculate the radial distribution of the surface density and representative size of icy dust particles that grow by colliding with each other and drift toward the central planet in a steady circumplanetary disk with a continuous supply of gas and dust from the parent protoplanetary disk. The radial drift barrier is overcome if the ratio of the dust to gas accretion rates onto the circumplanetary disk, $\dot{M}_{\mathrm{d}}/\dot{M}_{\mathrm{g}}$, is high and the strength of turbulence, $\alpha$, is not too low. The collision velocity is lower than the critical velocity of fragmentation when $\alpha$ is low. Taken together, we find that the conditions for satellitesimal formation via dust coagulation are given by $\dot{M}_{\mathrm{d}}/\dot{M}_{\mathrm{g}}\ge1$ and $10^{-4}\le\alpha<10^{-2}$. The former condition is generally difficult to achieve, suggesting that the in-situ satellitesimal formation via particle sticking is viable only under an extreme condition. We also show that neither satellitesimal formation via the collisional growth of porous aggregates nor via streaming instability is viable as long as $\dot{M}_{\mathrm{d}}/\dot{M}_{\mathrm{g}}$ is low.

The recent rapid progress in observations of circumstellar disks and extrasolar planets has reinforced the importance of understanding an intimate coupling between star and planet formation. Under such a circumstance, it may be invaluable to attempt to specify when and how planet formation begins in star-forming regions and to identify what physical processes/quantities are the most significant to make a link between star and planet formation. To this end, we have recently developed a couple of projects. These include an observational project about dust growth in Class 0 YSOs and a theoretical modeling project of the HL Tauri disk. For the first project, we utilize the archive data of radio interferometric observations, and examine whether dust growth, a first step of planet formation, occurs in Class 0 YSOs. We find that while our observational results can be reproduced by the presence of large ($\sim$ mm) dust grains for some of YSOs under the single-component modified blackbody formalism, an interpretation of no dust growth would be possible when a more detailed model is used. For the second project, we consider an origin of the disk configuration around HL Tauri, focusing on magnetic fields. We find that magnetically induced disk winds may play an important role in the HL Tauri disk. The combination of these attempts may enable us to move towards a comprehensive understanding of how star and planet formation are intimately coupled with each other.

The mechanism of angular momentum transport in protoplanetary disks is fundamental to understand the distributions of gas and dust in the disks. The unprecedented, high spatial resolution ALMA observations taken toward HL Tau and subsequent radiative transfer modeling reveal that a high degree of dust settling is currently achieved at the outer part of the HL Tau disk. Previous observations however suggest a high disk accretion rate onto the central star. This configuration is not necessarily intuitive in the framework of the conventional viscous disk model, since efficient accretion generally requires a high level of turbulence, which can suppress dust settling considerably. We develop a simplified, semi-analytical disk model to examine under what condition these two properties can be realized in a single model. Recent, non-ideal MHD simulations are utilized to realistically model the angular momentum transport both radially via MHD turbulence and vertically via magnetically induced disk winds. We find that the HL Tau disk configuration can be reproduced well when disk winds are properly taken into account. While the resulting disk properties are likely consistent with other observational results, such an ideal situation can be established only if the plasma $\beta$ at the disk midplane is $\beta_0 \simeq 2 \times 10^4$ under the assumption of steady accretion. Equivalently, the vertical magnetic flux at 100 au is about 0.2 mG. More detailed modeling is needed to fully identify the origin of the disk accretion and quantitatively examine plausible mechanisms behind the observed gap structures in the HL Tau disk.

The total amount of dust (or "metallicity") and the dust distribution in protoplanetary disks are crucial for planet formation. Dust grains radially drift due to gas--dust friction, and the gas is affected by the feedback from dust grains. We investigate the effects of the feedback from dust grains on the viscous evolution of the gas, taking into account the vertical dust settling. The feedback from the grains pushes the gas outward. When the grains are small and the dust-to-gas mass ratio is much smaller than unity, the radial drift velocity is reduced by the feedback effect but the gas still drifts inward. When the grains are sufficiently large or piled-up, the feedback is so effective that forces the gas flows outward. Although the dust feedback is affected by dust settling, we found that the 2D approximation reasonably reproduces the vertical averaged flux of gas and dust. We also performed the 2D two-fluid hydrodynamic simulations to examine the effect of the feedback from the grains on the evolution of the gas disk. We show that when the feedback is effective, the gas flows outward and the gas density at the region within $\sim 10\ \mbox{AU}$ is significantly depleted. As a result, the dust-to-gas mass ratio at the inner radii may significantly excess unity, providing the environment where planetesimals are easily formed via, e.g., streaming instability. We also show that a simplified 1D model well reproduces the results of the 2D two-fluid simulations, which would be useful for future studies.

We investigate the dust structure of gravitationally unstable disks undergoing mass accretion from the envelope envisioning the application to Class 0/I young stellar objects (YSOs) We find that the dust disk quickly settles into a steady state and that, compared to a disk with interstellar medium (ISM) dust-to-gas mass ratio and micron-sized dust, the dust mass in the steady-state decreases by a factor of 1/2 to 1/3, and the dust thermal emission decreases by a factor of 1/3 to 1/5. The latter decrease is caused by dust depletion and opacity decrease owing to dust growth. Our results suggest that the masses of gravitationally unstable disks in the Class 0/I YSOs are underestimated by a factor of 1/3 to 1/5 when calculated from the dust thermal emission assuming an ISM dust-to-gas mass ratio and micron-sized dust opacity, and that a larger fraction of disks in Class 0/I YSOs is gravitationally unstable than was previously believed. We also investigate the orbital radius $r_{\rm P}$ within which planetesimals form via coagulation of porous dust aggregates and show that $r_{\rm P}$ becomes $\sim 20$ AU for a gravitationally unstable disk around a solar mass star. Because $r_{\rm P}$ increases as the gas surface density increases and a gravitationally unstable disk has a maximum gas surface density, $r_{\rm P}\sim 20$ AU is the theoretical maximum radius. We suggest that planetesimals formation in the Class 0/I phase is preferable to that in the Class II phase because large gas surface density is expected and large amount of dust is supplied by envelope-to-disk accretion.

We analytically derive the expressions for the structure of the inner region of protoplanetary disks based on the results from the recent hydrodynamical simulations. The inner part of a disk can be divided into four regions: dust-free region with gas temperature in the optically thin limit, optically thin dust halo, optically thick condensation front and the classical optically thick region in order from the inside. We derive the dust-to-gas mass ratio profile in the dust halo using the fact that partial dust condensation regulates the temperature to the dust evaporation temperature. Beyond the dust halo, there is an optically thick condensation front where all the available silicate gas condenses out. The curvature of the condensation surface is determined by the condition that the surface temperature must be nearly equal to the characteristic temperature $\sim 1200{\,\rm K}$. We derive the mid-plane temperature in the outer two regions using the two-layer approximation with the additional heating by the condensation front for the outermost region. As a result, the overall temperature profile is step-like with steep gradients at the borders between the outer three regions. The borders might act as planet traps where the inward migration of planets due to gravitational interaction with the gas disk stops. The temperature at the border between the two outermost regions coincides with the temperature needed to activate magnetorotational instability, suggesting that the inner edge of the dead zone must lie at this border. The radius of the dead-zone inner edge predicted from our solution is $\sim$ 2-3 times larger than that expected from the classical optically thick temperature.

A number of transiting exoplanets have featureless transmission spectra that might suggest the presence of clouds at high altitudes. A realistic cloud model is necessary to understand the atmospheric conditions under which such high-altitude clouds can form. In this study, we present a new cloud model that takes into account the microphysics of both condensation and coalescence. Our model provides the vertical profiles of the size and density of cloud and rain particles in an updraft for a given set of physical parameters, including the updraft velocity and the number density of cloud condensation nuclei (CCN). We test our model by comparing with observations of trade-wind cumuli on the Earth and ammonia ice clouds in Jupiter. For trade-wind cumuli, the model including both condensation and coalescence gives predictions that are consistent with observations, while the model including only condensation overestimates the mass density of cloud droplets by up to an order of magnitude. For Jovian ammonia clouds, the condensation-coalescence model simultaneously reproduces the effective particle radius, cloud optical thickness, and cloud geometric thickness inferred from Voyager observations if the updraft velocity and CCN number density are taken to be consistent with the results of moist convection simulations and Galileo probe measurements, respectively. These results suggest that the coalescence of condensate particles is important not only in terrestrial water clouds but also in Jovian ice clouds. Our model will be useful to understand how the dynamics, compositions, and nucleation processes in exoplanetary atmospheresaffects the vertical extent and optical thickness of exoplanetary clouds via cloud microphysics.

We report $\sim$3 au resolution imaging observations of the protoplanetary disk around TW Hya at 145 and 233 GHz with the Atacama Large Millimeter/Submillimeter Array. Our observations revealed two deep gaps ($\sim$25--50 \%) at 22 and 37~au and shallower gaps (a few \%) at 6, 28, and 44~au, as recently reported by Andrews et al. (2016). The central hole with a radius of $\sim3$~au was also marginally resolved. The most remarkable finding is that the spectral index $\alpha (R)$ between bands 4 and 6 peaks at the 22~au gap. The derived power-law index of the dust opacity $\beta (R)$ is $\sim1.7$ at the 22~au gap and decreases toward the disk center to $\sim0$. The most prominent gap at 22~au could be caused by the gravitational interaction between the disk and an unseen planet with a mass of $\lesssim$1.5 $M_\mathrm{Neptune}$ although other origins may be possible. The planet-induced gap is supported by the fact that $\beta (R)$ is enhanced at the 22~au gap, indicating a deficit of $\sim$mm-sized grains within the gap due to dust filtration by a planet.

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.

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.

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.

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 \mu 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%.

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.

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.

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 $\alpha$ 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.

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\pi$ in azimuth and $25 \mathrm{au}<r<50 \mathrm{au}$ is $1200\sigma$- to $4000\sigma$-detectable. The lower limits of the radii of $5\sigma$-detectable LMG clumps are between 1.6 au and 5.3 au, depending on the models.

The formation process of circumstellar disks is still controversial because of the interplay of complex physical processes that occurs during the gravitational collapse of prestellar cores. In this study, we investigate the effect of the Hall current term on the formation of the circumstellar disk using three- dimensional simulations. In our simulations, all non-ideal effects as well as the radiation transfer are considered. The size of the disk is significantly affected by a simple difference in the inherent properties of the prestellar core, namely whether the rotation vector and the magnetic field are parallel or anti- parallel. In the former case, only a very small disk (< 1 AU) is formed. On the other hand, in the latter case, a massive and large (> 20 AU) disk is formed in the early phase of protostar formation. Since the parallel and anti-parallel properties do not readily change, we expect that the parallel and anti-parallel properties are also important in the subsequent disk evolution and the difference between the two cases is maintained or enhanced. This result suggests that the disk size distribution of the Class 0 young stellar objects is bimodal. Thus, the disk evolution can be categorized into two cases and we may call the parallel and anti-parallel systems as Ortho-disk and Para-disk, respectively. We also show that the anti-rotating envelopes against the disk-rotation appear with a size of 200 AU. We predict that the anti-rotating envelope will be found in the future observations.

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 $\beta$ of the centre of the first core becomes large, $\beta>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, $\beta\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 $\mu=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.

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.

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.

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}, \tau_{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.

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 \mu 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 $\beta$. The previous observation of the radial distribution of $\beta$ 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.

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.

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.

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.

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.

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.

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.

Context: In protoplanetary disks, dust grains coagulate with each other and grow to form aggregates. As these aggregates grow by coagulation, their filling factor \phi decreases down to \phi << 1. However, comets, the remnants of these early planetesimals, have \phi ~ 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 (\phi > 0.1). Aims: We investigate and find the compression strength of highly porous aggregates (\phi << 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 (\phi << 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 (\phi > 0.1) and confirm that our results consistently connect to those in the high density region. The compression strength formula is also derived analytically.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

We investigate the nonlinear evolution of black hole ringdown in the framework of higher-order metric perturbation theory. By solving the initial-value problem of a simplified nonlinear field model analytically as well as numerically, we find that (i) second-order quasinormal modes (QNMs) are indeed excited at frequencies different from those of first-order QNMs, as predicted recently. We also find serendipitously that (ii) late-time evolution is dominated by a new type of power-law tail. This ``second-order power-law tail'' decays more slowly than any late-time tails known in the first-order (i.e., linear) perturbation theory, and is generated at the wavefront of the first-order perturbation by an essentially nonlinear mechanism. These nonlinear components should be particularly significant for binary black hole coalescences, and could open a new precision science in gravitational wave studies.

Quasinormal ringing of acoustic black holes in Laval nozzles is discussed. The equation for sounds in a transonic flow is written into a Schr\"{o}dinger-type equation with a potential barrier, and the quasinormal frequencies are calculated semianalytically. From the results of numerical simulations, it is shown that the quasinormal modes are actually excited when the transonic flow is formed or slightly perturbed, as well as in the real black hole case. In an actual experiment, however, the purely-outgoing boundary condition will not be satisfied at late times due to the wave reflection at the end of the apparatus, and a late-time ringing will be expressed as a superposition of "boxed" quasinormal modes. It is shown that the late-time ringing damps more slowly than the ordinary quasinormal ringing, while its central frequency is not greatly different from that of the ordinary one. Using this fact, an efficient way for experimentally detecting the quasinormal ringing of an acoustic black hole is discussed.