No Arabic abstract
The low water content of the terrestrial planets in the solar system suggests that the protoplanets formed within the water snow line. Accurate prediction of the snow line location moving with time provides a clue to constrain the formation process of the planets. In this paper, we investigate the migration of the snow line in protoplanetary disks whose accretion is controlled by laminar magnetic fields, which have been proposed by various nonideal magnetohydrodynamic (MHD) simulations. We propose an empirical model of the disk temperature based on our nonideal MHD simulations, which show that the accretion heating is significantly less efficient than in turbulent disks, and calculate the snow line location over time. We find that the snow line in the magnetically accreting laminar disks moves inside the current Earths orbit within 1 Myr after star formation, whereas the time for the conventional turbulent disk is much longer than 1 Myr. This result suggests that either the rocky protoplanets formed in such an early phase of the disk evolution, or the protoplanets moved outward to the current orbits after they formed close to the protosun.
[Abridged] Planet formation is expected to be enhanced around snowlines in protoplanetary disks, in particular around the water snowline. However, the close proximity of the water snowline to the host star and water in the Earths atmosphere makes a direct detection of the water snowline in protoplanetary disks challenging. Following earlier work on protostellar envelopes, the aim of this research is to investigate the validity of HCO$^+$ and H$^{13}$CO$^+$, as tracers of the water snowline in protoplanetary disks, as HCO$^+$ is destroyed by gas-phase water. Two small chemical networks are used to predict the HCO$^+$ abundance in a typical Herbig Ae disk. Subsequently, the corresponding emission profiles are modelled for H$^{13}$CO$^+$ and HCO$^+$ $J=2-1$, which provides the best balance between brightness and optical depth effects of the continuum emission. The HCO$^+$ abundance jumps by two orders of magnitude just outside the water snowline at 4.5 AU. We find that the emission of H$^{13}$CO$^+$ and HCO$^+$ is ring-shaped due to three effects: destruction of HCO$^+$ by gas-phase water, continuum optical depth, and molecular excitation effects. The presence of gas-phase water causes an additional drop of only $sim$13% and 24% in the center of the disk, for H$^{13}$CO$^+$ and HCO$^+$, respectively. For the much more luminous outbursting source V883Ori, our models predict that the effect of dust and excitation are not limiting if the snowline is located outside $sim$40 AU. Our analysis of ALMA observations of HCO$^+$ $J=3-2$ is consistent with the water snowline located around 100 AU. The HCO$^+$ abundance drops steeply around the water snowline, but dust and excitation can conceal the drop in HCO$^+$ emission due to the water snowline. Therefore, locating the water snowline with HCO$^+$ in Herbig disks is very difficult, but it is possible for outbursting sources like V883Ori.
We investigate the roles of magnetically driven disk wind (MDW) and thermally driven photoevaporative wind (PEW) in the long-time evolution of protoplanetary disks. We start simulations from the early phase in which the disk mass is $0.118,{mathrm{M}_{odot}}$ around a $1,{mathrm{M}_{odot}}$ star and track the evolution until the disk is completely dispersed. We incorporate the mass loss by PEW and the mass loss and magnetic braking (wind torque) by MDW, in addition to the viscous accretion, viscous heating, and stellar irradiation. We find that MDW and PEW respectively have different roles: magnetically driven wind ejects materials from an inner disk in the early phase, whereas photoevaporation has a dominant role in the late phase in the outer ($gtrsim1,$au) disk. The disk lifetime, which depends on the combination of MDW, PEW, and viscous accretion, shows a large variation of $sim1$-$20,$Myr; the gas is dispersed mainly by the MDW and the PEW in the cases with a low viscosity and the lifetime is sensitive to the mass-loss rate and torque of the MDW, whereas the lifetime is insensitive to these parameters when the viscosity is high. Even in disks with very weak turbulence, the cooperation of MDW and PEW enables the disk dispersal within a few Myr.
Mid-infrared molecular line emission detected with the Spitzer Space Telescope is often interpreted using slab models. However, we need to understand the mid-infrared line emission in 2D disk models, such that we gain information about from where the lines are being emitted and under which conditions, such that we gain information about number densities, temperatures, and optical depths in both the radial and vertical directions. In this paper, we introduce a series of 2D thermochemical models of a prototypical T Tauri protoplanetary disk, in order to examine how sensitive the line-emitting regions are to changes in the UV and X-ray fluxes, the disk flaring angle, dust settling, and the dust-to-gas ratio. These all affect the heating of the inner disk, and thus can affect the mid-infrared spectral lines. Using the ProDiMo and FLiTs codes, we produce a series of 2D thermochemical disk models. We find that there is often a significant difference between the gas and dust temperatures in the line emitting regions, and we illustrate that the size of the line emitting regions is relatively robust against changes in the stellar and disk parameters (namely, the UV and X-ray fluxes, the flaring angle, and dust settling). These results demonstrate the potential for localized variations in the line-emitting region to greatly affect the resulting spectra and line fluxes, and the necessity of allowing for such variations in our models.
It is usually thought that viscous torque works to align a circumbinary disk with the binarys orbital plane. However, recent numerical simulations suggest that the disk may evolve to a configuration perpendicular to the binary orbit (polar alignment) if the binary is eccentric and the initial disk-binary inclination is sufficiently large. We carry out a theoretical study on the long-term evolution of inclined disks around eccentric binaries, calculating the disk warp profile and dissipative torque acting on the disk. For disks with aspect ratio $H/r$ larger than the viscosity parameter $alpha$, bending wave propagation effectively makes the disk precess as a quasi-rigid body, while viscosity acts on the disk warp and twist to drive secular evolution of the disk-binary inclination. We derive a simple analytic criterion (in terms of the binary eccentricity and initial disk orientation) for the disk to evolve toward polar alignment with the eccentric binary. When the disk has a non-negligible angular momentum compared to the binary, the final polar alignment inclination angle is reduced from $90^circ$. For typical protoplanetary disk parameters, the timescale of the inclination evolution is shorter than the disk lifetime, suggesting that highly-inclined disks and planets may exist orbiting eccentric binaries.
We discovered a new growth mode of dust grains to km-sized bodies in protoplanetary disks that evolve by viscous accretion and magnetically driven disk winds (MDWs). We solved an approximate coagulation equation of dust grains with time-evolving disks that consist of both gas and solid components by a one-dimensional model. With the grain growth, all solid particles initially drift inward toward the central star by the gas drag force. However, the radial profile of gas pressure, $P$, is modified by the MDW that disperses the gas in an inside-out manner. Consequently, a local concentration of solid particles is created by the converging radial flux of drifting dust grains at the location with the convex upward profile of $P$. When the dimensionless stopping time, ${rm St}$, there exceeds unity, the solid particles spontaneously reach the growth dominated state because of the positive feedback between the suppressed radial drift and the enhanced accumulation of dust particles that drift from the outer part. Once the solid particles are in the drift limited state, the above-mentioned condition of ${rm St} gtrsim 1$ for the dust growth is equivalent with begin{equation} Sigma_{rm d}/Sigma_{rm g}gtrsim eta, onumber end{equation} where $Sigma_{rm d}/Sigma_{rm g}$ is the dust-to-gas surface-density ratio and $eta$ is dimensionless radial pressure-gradient force. As a consequence of the successful growth of dust grains, a ring-like structure containing planetesimal-sized bodies is formed at the inner part of the protoplanetary disks. Such a ring-shaped concentration of planetesimals is expected to play a vital role in the subsequent planet formation.