No Arabic abstract
Alignment of non-spherical grains with magnetic fields is an important problem as it lays the foundation of probing magnetic fields with polarized dust thermal emissions. In this paper, we investigate the feasibility of magnetic alignment in protoplanetary disks (PPDs). We use an alignment condition that Larmor precession should be fast compared with the damping timescale. We first show that the Larmor precession timescale is some three orders of magnitude longer than the damping time for millimeter-sized grains under conditions typical of PPDs, making the magnetic alignment unlikely. The precession time can be shortened by superparamagnetic inclusions (SPIs), but the reduction factor strongly depends on the size of the SPI clusters, which we find is limited by the so-called N{e}els relaxation process. In particular, the size limit of SPIs is set by the so-called anisotropic energy constant of the SPI material, which describes the energy barrier needed to change the direction of the magnetic moment of an SPI. For the most common iron-bearing materials, we find maximum SPI sizes corresponding to a reduction factor of the Larmor precession timescale of order $10^3$. We also find that reaching this maximum reduction factor requires fine-tuning on the SPI sizes. Lastly, we illustrate the effects of the SPI size limits on magnetic alignment of dust grains with a simple disk model, and we conclude that it is unlikely for relatively large grains of order 100 $mu$m or more to be aligned with magnetic fields even with SPIs.
The mechanisms causing millimeter-wave polarization in protoplanetary disks are under debate. To disentangle the polarization mechanisms, we observe the protoplanetary disk around HL Tau at 3.1 mm with the Atacama Large Millimeter/submillimeter Array (ALMA), which had polarization detected with CARMA at 1.3 mm. We successfully detect the ring-like azimuthal polarized emission at 3.1 mm. This indicates that dust grains are aligned with the major axis being in the azimuthal direction, which is consistent with the theory of radiative alignment of elongated dust grains, where the major axis of dust grains is perpendicular to the radiation flux. Furthermore, the morphology of the polarization vectors at 3.1 mm is completely different from those at 1.3 mm. We interpret that the polarization at 3.1 mm to be dominated by the grain alignment with the radiative flux producing azimuthal polarization vectors, while the self-scattering dominates at 1.3 mm and produces the polarization vectors parallel to the minor axis of the disk. By modeling the total polarization fraction with a single grain population model, the maximum grain size is constrained to be $100{rm~mu m}$, which is smaller than the previous predictions based on the spectral index between ALMA at 3 mm and VLA at 7 mm.
High-energy irradiation of the circumstellar material might impact the structure and the composition of a protoplanetary disk and hence the process of planet formation. In this paper, we present a study on the possible influence of the stellar irradiation, indicated by X-ray emission, on the crystalline structure of the circumstellar dust. The dust crystallinity is measured for 42 class II T Tauri stars in the Taurus star-forming region using a decomposition fit of the 10 micron silicate feature, measured with the Spitzer IRS instrument. Since the sample includes objects with disks of various evolutionary stages, we further confine the target selection, using the age of the objects as a selection parameter. We correlate the X-ray luminosity and the X-ray hardness of the central object with the crystalline mass fraction of the circumstellar dust and find a significant anti-correlation for 20 objects within an age range of approx. 1 to 4.5 Myr. We postulate that X-rays represent the stellar activity and consequently the energetic ions of the stellar winds which interact with the circumstellar disk. We show that the fluxes around 1 AU and ion energies of the present solar wind are sufficient to amorphize the upper layer of dust grains very efficiently, leading to an observable reduction of the crystalline mass fraction of the circumstellar, sub-micron sized dust. This effect could also erase other relations between crystallinity and disk/star parameters such as age or spectral type.
We study the effects of grain surface reactions on the chemistry of protoplanetary disks where gas, ice surface layers and icy mantles of dust grains are considered as three distinct phases. Gas phase and grain surface chemistry is found to be mainly driven by photo-reactions and dust temperature gradients. The icy disk interior has three distinct chemical regions: (i) the inner midplane with low FUV fluxes and warm dust ($gtrsim 15$K) that lead to the formation of complex organic molecules, (ii) the outer midplane with higher FUV from the ISM and cold dust where hydrogenation reactions dominate and, (iii) a molecular layer above the midplane but below the water condensation front where photodissociation of ices affects gas phase compositions. Some common radicals, e.g., CN and C$_2$H, exhibit a two-layered vertical structure and are abundant near the CO photodissociation front and near the water condensation front. The 3-phase approximation in general leads to lower vertical column densities than 2-phase models for many gas-phase molecules due to reduced desorption, e.g., H$_2$O, CO$_2$, HCN and HCOOH decrease by $sim$ two orders of magnitude. Finally, we find that many observed gas phase species originate near the water condensation front; photo-processes determine their column densities which do not vary significantly with key disk properties such as mass and dust/gas ratio.
ALMA observations of protoplanetary disks confirm earlier indications that there is a clear difference between the dust and gas radial extents. The origin of this difference is still debated, with both radial drift of the dust and optical depth effects suggested in the literature. In this work, the feedback of realistic dust particle distributions onto the gas chemistry and molecular emissivity is investigated, with a particular focus on CO isotopologues. The radial dust grain size distribution is determined using dust evolution models that include growth, fragmentation and radial drift. A new version of the code DALI is used to take into account how dust surface area and density influence the disk thermal structure, molecular abundances and excitation. The difference of dust and gas radial sizes is largely due to differences in the optical depth of CO lines and millimeter continuum, without the need to invoke radial drift. The effect of radial drift is primarily visible in the sharp outer edge of the continuum intensity profile. The gas outer radius probed by $^{12}$CO emission can easily differ by a factor of $sim 2$ between the models for a turbulent $alpha$ ranging between typical values. Grain growth and settling concur in thermally decoupling the gas and dust components, due to the low collision rate with large grains. As a result, the gas can be much colder than the dust at intermediate heights, reducing the CO excitation and emission, especially for low turbulence values. Also, due to disk mid-plane shadowing, a second CO thermal desorption (rather than photodesorption) front can occur in the warmer outer mid-plane disk. The models are compared to ALMA observations of HD 163296 as a test case. In order to reproduce the observed CO snowline of the system, a binding energy for CO typical of ice mixtures needs to be used rather than the lower pure CO value.
Grain surface chemistry is key to the composition of protoplanetary disks around young stars. The temperature of grains depends on their size. We evaluate the impact of this temperature dependence on the disk chemistry. We model a moderately massive disk with 16 different grain sizes. We use POLARIS to calculate the dust grain temperatures and the local UV flux. We model the chemistry using the 3-phase astrochemical code NAUTILUS. Photoprocesses are handled using frequency-dependent cross-sections, and a new method to account for self and mutual shielding. The multi-grain model outputs are compared to those of single-grain size models (0.1 $mu$m), with two different assumptions for their equivalent temperature. We find that the Langmuir-Hinshelwood (LH) mechanism at equilibrium temperature is not efficient to form H$_2$ at 3-4 scale heights ($H$), and adopt a parametric fit to a stochastic method to model H$_2$ formation instead. We find the molecular layer composition (1-3 $H$) to depend on the amount of remaining H atoms. Differences in molecular surface densities between single and multi-grain models are mostly due to what occurs above 1.5 $H$. At 100 au, models with colder grains produce H$_2$O and CH$_4$ ices in the midplane, and warmer ones produce more CO$_2$ ices, both allowing efficient depletion of C and O as soon as CO sticks on grain surfaces. Complex organic molecules (COMs) production is enhanced by the presence of warmer grains in the multi-grain models. Using a single grain model mimicking grain growth and dust settling fails to reproduce the complexity of gas-grain chemistry. Chemical models with a single grain size are sensitive to the adopted grain temperature, and cannot account for all expected effects. A spatial spread of the snowlines is expected to result from the ranges in grain temperature. The amplitude of the effects will depend on the dust disk mass.