The formation and early evolution of low mass young stellar objects (YSOs) are investigated using three-dimensional non-ideal magneto-hydrodynamics simulations. We investigate the evolution of YSOs up to ~ 10^4 yr after protostar formation, at which protostellar mass reaches ~ 0.1 M_\odot . We particularly focus on the impact of the dust model on the evolution. We found that a circumstellar disk is formed in all simulations regardless of the dust model. Disk size is approximately 10 AU at the protostar formation epoch, and it increases to several tens of AU at ~ 10^4 yr after protostar formation. Disk mass is comparable to central protostellar mass and gravitational instability develops. In the simulations with small dust size, the warp of the pseudodisk develops ~ 10^4 yr after protostar formation. The warp strengthens magnetic braking in the disk and decreases disk size. Ion-neutral drift can occur in the infalling envelope under the conditions that the typical dust size is a \gtrsim 0.2{\mu}m and the protostar (plus disk) mass is M \gtrsim 0.1 M_\odot. The outflow activity is anti-correlated to the dust size and the strong outflow appears with small dust grains.

The cometary materials are thought to be the reservoir of primitive materials in the solar system. The detection of glycine and CH$_3$NH$_2$ by the ROSINA mass spectrometer in the coma of 67P/Churyumov-Gerasimenko suggests that amino acids and their precursors may have been formed in an early evolutionary phase of the Solar System. We investigate if we can explain "Glycine/CH$_3$NH$_2$" ratio by assuming that glycine and CH$_3$NH$_2$ are inherited from the parentual molecular cloud, and the importance of the concentration process of non-volatile species. We use our astrochemical model representative of physio-chemical evolution during the low-mass (or Solar-type) star-formation to investigate if we can explain the observed "Glycine/CH$_3$NH$_2$" ratio. In addition, we develop a Simplified Cometary Concentration Model (SCCM) to discuss this concentration process with the hypothesis that glycine concentrates on the cometary surface during the numerous passages of perihelion after the formation of the Solar System. The predicted molecular ratio of "Glycine/CH$_3$NH$_2$" by our astrochemical model is less than 10$^{-2}$, while the ROSINA measurement has shown this value to be unity. On the other hand, SCCM results suggest that non-volatile species highly likely concentrate the cometary surface due to its high binding energy. If a small amount of non-volatile species is ejected via comet activity, "Glycine/CH$_3$NH$_2$" ratio would be highly enhanced.

GW Ori is a hierarchical triple system which has a rare circumtriple disk. We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of 1.3 mm dust continuum and 12CO J=2-1 molecular gas emission of the disk. For the first time, we identify three dust rings in the disk at ~46, 188, and 338 AU, with estimated dust mass of ~70-250 Earth masses, respectively. To our knowledge, the outer ring in GW Ori is the largest dust ring ever found in protoplanetary disks. We use visibility modelling of dust continuum to show that the disk has misaligned parts and the innermost dust ring is eccentric. The disk misalignment is also suggested by the CO kinematics modelling. We interpret these substructures as evidence of ongoing dynamical interactions between the triple stars and the circumtriple disk.

A planet is formed within a protoplanetary disk. Recent observations have revealed substructures such as gaps and rings, which may indicate forming planets within the disk. Due to disk--planet interaction, the planet migrates within the disk, which can affect a shape of the planet-induced gap. In this paper, we investigate effects of fast inward migration of the planet on the gap shape, by carrying out hydrodynamic simulations. We found that when the migration timescale is shorter than the timescale of the gap-opening, the orbital radius is shifted inward as compared to the radial location of the gap. We also found a scaling relation between the radial shift of the locations of the planet and the gap as a function of the ratio of the timescale of the migration and gap-opening. Our scaling relation also enables us to constrain the gas surface density and the viscosity when the gap and the planet are observed. Moreover, we also found the scaling relation between the location of the secondary gap and the aspect ratio. By combining the radial shift and the secondary gap, we may constrain the physical condition of the planet formation and how the planet evolves in the protoplanetary disk, from the observational morphology.

We observe the dust continuum at 225 GHz and CO isotopologue (12CO, 13CO, and C18O) J=2-1 emission lines toward the CR Cha protoplanetary disk using the Atacama Large Millimeter/Submillimeter Array (ALMA). The dust continuum image shows a dust gap-ring structure in the outer region of the dust disk. A faint dust ring is also detected around 120 au beyond the dust gap. The CO isotopologue lines indicate that the gas disk is more extended than the dust disk. The peak brightness temperature of the 13CO line shows a small bump around 130 au while 12CO and C18O lines do not show. We investigate two possible mechanisms for reproducing the observed dust gap-ring structure and a gas temperature bump. First, the observed gap structure can be opened by a Jupiter mass planet using the relation between the planet mass and the gap depth and width. Meanwhile, the radiative transfer calculations based on the observed dust surface density profile show that the observed dust ring could be formed by dust accumulation at the gas temperature bump, that is, the gas pressure bump produced beyond the outer edge of the dust disk.

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.

We present the observational result of a glycine precursor, methylamine (CH$_3$NH$_2$), together with methanol (CH$_3$OH) and methanimine (CH$_2$NH) towards high-mass star-forming regions, NGC6334I, G10.47+0.03, G31.41+0.3, and W51~e1/e2 using ALMA. The molecular abundances toward these sources were derived using the rotational diagram method and compared with our state-of-the-art chemical model. We found that the observed ratio of "CH$_3$NH$_2$/CH$_3$OH" is in between 0.11 and 2.2. We also found that the observed "CH$_3$NH$_2$/CH$_3$OH" ratio agrees well with our chemical model by considering the formation of CH$_3$NH$_2$ on the grain surface via hydrogenation process to HCN. This result clearly shows the importance of hydrogenation processes to form CH$_3$NH$_2$. NGC63343I MM3, where CH$_3$NH$_2$ was not detected in this study and showed "CH$_3$NH$_2$/CH$_3$OH" ratio of less than 0.02, is clearly distinguished from the other cores.

Measurements of the gas mass are necessary to determine the planet formation potential of protoplanetary disks. Observations of rare CO isotopologues are typically used to determine disk gas masses; however, if the line emission is optically thick this will result in an underestimated disk mass. With ALMA we have detected the rarest stable CO isotopologue, 13C17O, in a protoplanetary disk for the first time. We compare our observations with the existing detections of 12CO, 13CO, C18O and C17O in the HD163296 disk. Radiative transfer modelling using a previously benchmarked model, and assuming interstellar isotopic abundances, significantly underestimates the integrated intensity of the 13C17O J=3-2 line. Reconciliation between the observations and the model requires a global increase in CO gas mass by a factor of 3.5. This is a factor of 2-6 larger than previous gas mass estimates using C18O. We find that C18O emission is optically thick within the CO snow line, while the 13C17O emission is optically thin and is thus a robust tracer of the bulk disk CO gas mass.

We analyze $^{12}$CO($J$=2-1), $^{13}$CO($J$=2-1), C$^{18}$O ($J$=2-1), and 1.3 mm continuum maps of the $\rho$ Ophiuchi A photo-dissociation region (PDR) obtained with ALMA. Layered structures of the three CO isotopologues with an angular separation of 10 arcsec = 6.6$\times$10$^{-3}$ pc = 1400 au are clearly detected around the Be star, S1 (i.e., each front of emission shifts from the near to far side in order of $^{12}$CO, $^{13}$CO, and C$^{18}$O). We estimate the spatial variations of $X$($^{13}$CO)/$X$(C$^{18}$O) abundance ratios, and find that the abundance ratio is as high as 40 near the emission front, and decreases to the typical value in the solar system of 5.5 in a small angular scale of 4 arcsec = 2.6$\times$10$^{-3}$ pc = 560 au. We also find that the $I$($^{12}$CO(2-1))/$I$($^{13}$CO(2-1)) intensity ratio is very high ($>$21) in the flat-spectrum young stellar object, GY-51, located in the PDR. The enhancement of the ratios indicates that the UV radiation significantly affects the CO isotopologues via selective dissociation in the overall $\rho$ Ophiuchi A PDR, and that the $\rho$ Ophiuchi A PDR has a plane-parallel structure.

In this paper, we analyse the upper limit fluxes of sub-millimeter ortho-H$_{2}$$^{16}$O 321 GHz, para-H$_{2}$$^{18}$O 322 GHz, and HDO 335 GHz lines from the protoplanetary disk around the Herbig Ae star HD 163296, using the Atacama Large Millimeter/Submillimeter Array (ALMA). These water lines are considered to be the best candidate sub-millimeter lines to locate the position of the H$_{2}$O snowline, on the basis of our previous model calculations. We compare the upper limit fluxes with the values calculated by our models with dust emission included, and we constrain the line emitting region and the dust opacity from the observations. We conclude that, if the outer edge of the region with high water vapor abundance and if the position of the water snowline are beyond 8 au also, the mm dust opacity will have a value larger than 2.0 cm$^{2}$ g$^{-1}$. In addition, the position of the water snowline will be inside 20 au, if the mm dust opacity is 2.0 cm$^{2}$ g$^{-1}$. Future observations of the dust continuum emission at higher angular resolution and sub-millimeter water lines with longer observation time are required to clarify the detailed structures and the position of the H$_{2}$O snowline in the disk midplane.

Analyzing multiband observations of dust continuum emission is one of the useful tools to constrain dust properties which help us to understand the physical properties of the disks. We perform the synthetic ALMA multiband analysis to find the best ALMA band set for constraining the dust properties of the TW Hya protoplanetary disk. We find that the Band [10,6,3] set is the best set among the possible combinations of ALMA Band [3,4,5,6,7,8,9,10]. We also find two conditions for the good ALMA band sets providing narrow constraint ranges on dust properties; (1) Band 9 or 10 is included in the band set and (2) Enough frequency intervals between the bands. These are related with the conditions which give good constraints on dust properties: the combination of optically thick and thin bands are required, and large beta (beta is the power-law index of dust opacity), and low dust temperature are preferable. To examine our synthetic analysis results, we apply the multiband analysis to ALMA archival data of the TW Hya disk at Band 4, 6, 7, and 9. Band [9,6,4] set provides the dust properties close to the model profile, while Band [7,6,4] set gives the dust properties deviating from the model at all radii with too broad constraint range to specify the accurate values of dust temperature, optical depth, and \b{eta}. Since these features are expected by the synthetic multiband analysis, we confirm that the synthetic multiband analysis is well consistent with the results derived from real data.

We examined the excitations of planetesimals caused by the resonances of a giant planet in a protoplanetary gas disk. The highly excited planetesimals generate bow shocks, the mechanism of which results in chondrule formation, crystallization of silicate dust, and evaporation of icy planetesimals. The planetesimals beyond 2:1 resonance migrate owing to the gas drag and obtain the maximum eccentricity around 3:1 resonance, which is located at approximately half the planetary distance. The eccentricity depends on the parameters of the planetesimals and the Jovian planet, such as size and location, and gas density of the disk. The maximum relative velocity of a 100-km-sized planetesimal with respect to the gas disk reaches up to ~12 km/s in the case of Jupiter owing to secular resonance, which occurs because of the disk's gravity. We find that if a Jovian mass planet is located within 10 au, the planetesimals larger than 100 km gain sufficient velocity to cause the melting of chondrule precursors and crystallization of the silicate. The maximum velocity is higher for large planetesimals and eccentric planets. Planetesimals are trapped temporarily in the resonances and continue to have high speed over >1 Myr after the formation of a Jovian planet. This duration fits into the timescale of chondrule formation suggested by the isotopic data. The evaporation of icy planetesimals occurs when a Jovian planet is located within 15 au. This mechanism can be a new indicator of planet formation in exosystems if some molecules ejected from icy planetesimals are detected.

The bulk composition of the Earth is dramatically carbon poor compared to that of the interstellar medium, and this phenomenon extends to the asteroid belt. To interpret this carbon deficit problem, the carbonaceous component in grains must have been converted into the gas-phase in the inner regions of protoplanetary disks prior to planetary formation. We examine the effect of carbon grain destruction on the chemical structure of disks by calculating the molecular abundances and distributions using a comprehensive chemical reaction network. When carbon grains are destroyed and the elemental abundance of the gas becomes carbon-rich, the abundances of carbon-bearing molecules, such as HCN and carbon-chain molecules, increase dramatically near the midplane, while oxygen-bearing molecules, such as H$_2$O and CO$_2$, are depleted. We compare the results of these model calculations with the solid carbon fraction in the solar system. Although we find a carbon depletion gradient, there are some quantitative discrepancies: the model shows a higher value at the position of asteroid belt and a lower value at the location of the Earth. In addition, using the obtained molecular abundances distributions, coupled with line radiative transfer calculations, we make predictions for ALMA to potentially observe the effect of carbon grain destruction in nearby protoplanetary disks. The results indicate that HCN, H$^{13}$CN and c-C$_3$H$_2$ may be good tracers.

We perform a suite of radiation hydrodynamics simulations of photoevaporating disks with varying the metallicity in a wide range of $10^{-3} \, Z_\odot \leq Z \leq 10^{0.5} \, Z_\odot $. We follow the disk evolution for over $\sim 5000$ years by solving hydrodynamics, radiative transfer, and non-equilibrium chemistry. Our chemistry model is updated from the first paper of this series by adding X-ray ionization and heating. We study the metallicity dependence of the disk photoevaporation rate and examine the importance of X-ray radiation. In the fiducial case with solar metallicity, including the X-ray effects does not significantly increase the photoevaporation rate when compared to the case with ultra-violet (UV) radiation only. At sub-solar metallicities in the range of $Z \gtrsim 10^{-1.5} \, Z_\odot $, the photoevaporation rate increases as metallicity decreases owing to the reduced opacity of the disk medium. The result is consistent with the observational trend that disk lifetimes are shorter in low metallicity environments. Contrastingly, the photoevaporation rate decreases at even lower metallicities of $Z \lesssim 10^{-1.5} \, Z_\odot $, because dust-gas collisional cooling remains efficient compared to far UV photoelectric heating whose efficiency depends on metallicity. The net cooling in the interior of the disk suppresses the photoevaporation. However, adding X-ray radiation significantly increases the photoevaporation rate, especially at $Z \sim 10^{-2}\, Z_\odot$. Although the X-ray radiation itself does not drive strong photoevaporative flows, X-rays penetrate deep into the neutral region in the disk, increase the ionization degree there, and reduce positive charges of grains. Consequently, the effect of photoelectric heating by far UV radiation is strengthened by the X-rays and enhances the disk photoevaporation.

In this paper, we extend the results presented in our former papers (Notsu et al. 2016, 2017) on using ortho-$\mathrm{H_2}$$^{16}\mathrm{O}$ line profiles to constrain the location of the $\mathrm{H_2O}$ snowline in T Tauri and Herbig Ae disks, to include sub-millimeter para-$\mathrm{H_2}$$^{16}\mathrm{O}$ and ortho- and para-$\mathrm{H_2}$$^{18}\mathrm{O}$ lines. Since the number densities of the ortho- and para-H$_{2}$$^{18}$O molecules are about 560 times smaller than their $^{16}$O analogues, they trace deeper into the disk than the ortho-H$_{2}$$^{16}$O lines (down to $z=0$, i.e., the midplane). Thus these H$_{2}$$^{18}$O lines are potentially better probes of the position of the H$_{2}$O snowline at the disk midplane, depending on the dust optical depth. The values of the Einstein $A$ coefficients of sub-millimeter candidate water lines tend to be lower (typically $<$$10^{-4}$ s$^{-1}$) than infrared candidate water lines (Notsu et al. 2017). Thus in the sub-millimeter candidate water line cases, the local intensity from the outer optically thin region in the disk is around $10^{4}$ times smaller than that in the infrared candidate water line cases. Therefore, in the sub-millimeter lines, especially H$_{2}$$^{18}$O and para-H$_{2}$$^{16}$O lines with relatively lower upper state energies ($\sim$ a few 100K) can also locate the position of the $\mathrm{H_2O}$ snowline. We also investigate the possibility of future observations with ALMA to identify the position of the water snowline. There are several candidate water lines that trace the hot water vapor inside the $\mathrm{H_2O}$ snowline in ALMA Bands $5-10$.

Protoplanetary disks are thought to have lifetimes of several million years in the solar neighborhood, but recent observations suggest that the disk lifetimes are shorter in a low metallicity environment. We perform a suite of radiation hydrodynamics simulations of photoevaporation of protoplanetary disks to study the disk structure and its long-term evolution of $\sim 10000$ years, and the metallicity dependence of mass-loss rate. Our simulations follow hydrodynamics, extreme and far ultra-violet radiative transfer, and non-equilibrium chemistry in a self-consistent manner. Dust grain temperatures are also calculated consistently by solving the radiative transfer of the stellar irradiation and grain (re-)emission. We vary the disk gas metallicity over a wide range of $10^{-4}~ Z_\odot \leq Z \leq 10 ~Z_\odot$. The photoevaporation rate is lower with higher metallicity in the range of $10^{-1} \,Z_\odot \lesssim Z \lesssim 10 \,Z_\odot$, because dust shielding effectively prevents far-ultra violet (FUV) photons from penetrating into and heating the dense regions of the disk. The photoevaporation rate sharply declines at even lower metallicities in $10^{-2} \,Z_\odot \lesssim Z \lesssim 10^{-1}\,Z_\odot$, because FUV photoelectric heating becomes less effective than dust-gas collisional cooling. The temperature in the neutral region decreases, and photoevaporative flows are excited only in an outer region of the disk. At $10^{-4}\,Z_\odot \leq Z \lesssim 10^{-2}\,Z_\odot$, HI photoionization heating acts as a dominant gas heating process and drives photoevaporative flows with roughly a constant rate. The typical disk lifetime is shorter at $Z=0.3~Z_\odot$ than at $Z = Z_\odot$, being consistent with recent observations of the extreme outer galaxy.

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

Observationally measuring the location of the H$_{2}$O snowline is crucial for understanding the planetesimal and planet formation processes, and the origin of water on Earth. In disks around Herbig Ae stars ($T_{\mathrm{*}}\sim$ 10,000K, $M_{\mathrm{*}}\gtrsim$ 2.5$M_{\bigodot}$), the position of the H$_{2}$O snowline is further from the central star compared with that around cooler, and less massive T Tauri stars. Thus, the H$_{2}$O emission line fluxes from the region within the H$_{2}$O snowline are expected to be stronger. In this paper, we calculate the chemical composition of a Herbig Ae disk using chemical kinetics. Next, we calculate the H$_{2}$O emission line profiles, and investigate the properties of candidate water lines across a wide range of wavelengths (from mid-infrared to sub-millimeter) that can locate the position of the H$_{2}$O snowline. Those line identified have small Einstein $A$ coefficients ($\sim 10^{-6} -10^{-3}$ s$^{-1}$) and relatively high upper state energies ($\sim$ 1000K). The total fluxes tend to increase with decreasing wavelengths. We investigate the possibility of future observations (e.g., ALMA, SPICA/SMI-HRS) to locate the position of the H$_{2}$O snowline. Since the fluxes of those identified lines from Herbig Ae disks are stronger than those from T Tauri disks, the possibility of a successful detection is expected to increase for a Herbig Ae disk.

Transitional disks show a lack of excess emission at infrared wavelengths due to a large dust cavity, that is often corroborated by spatially resolved observations at ~ mm wavelengths. We present the first spatially resolved ~ mm-wavelength images of the disk around the Herbig Ae/Be star, HD 97048. Scattered light images show that the disk extends to ~640 au. The ALMA data reveal a circular-symmetric dusty disk extending to ~350 au, and a molecular disk traced in CO J=3-2 emission, extending to ~750 au. The CO emission arises from a flared layer with an opening angle ~ 30 deg - 40 deg. HD 97048 is another source for which the large (~ mm-sized) dust grains are more centrally concentrated than the small (~ {\mu}m-sized) grains and molecular gas, likely due to radial drift. The images and visibility data modelling suggests a decrement in continuum emission within ~50 au, consistent with the cavity size determined from mid-infrared imaging (34 +/- 4 au). The extracted continuum intensity profiles show ring-like structures with peaks at ~50, 150, and 300 au, with associated gaps at ~100 and 250 au. This structure should be confirmed in higher-resolution images (FWHM ~ 10 - 20 au). These data confirm the classification of HD 97048 as a transitional disk that also possesses multiple ring-like structures in the dust continuum emission. Additional data are required at multiple and well-separated frequencies to fully characterise the disk structure, and thereby constrain the mechanism(s) responsible for sculpting the HD 97048 disk.

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.

The first detection of gas-phase methanol in a protoplanetary disk (TW Hya) is presented. In addition to being one of the largest molecules detected in disks to date, methanol is also the first disk organic molecule with an unambiguous ice chemistry origin. The stacked methanol emission, as observed with ALMA, is spectrally resolved and detected across six velocity channels ($>3 \sigma$), reaching a peak signal-to-noise of $5.5\sigma$, with the kinematic pattern expected for TW~Hya. Using an appropriate disk model, a fractional abundance of $3\times 10^{-12} - 4 \times 10^{-11}$ (with respect to H$_2$) reproduces the stacked line profile and channel maps, with the favoured abundance dependent upon the assumed vertical location (midplane versus molecular layer). The peak emission is offset from the source position suggesting that the methanol emission has a ring-like morphology: the analysis here suggests it peaks at $\approx 30$~AU reaching a column density $\approx 3-6\times10^{12}$~cm$^{-2}$. In the case of TW Hya, the larger (up to mm-sized) grains, residing in the inner 50~AU, may thus host the bulk of the disk ice reservoir. The successful detection of cold gas-phase methanol in a protoplanetary disk implies that the products of ice chemistry can be explored in disks, opening a window to studying complex organic chemistry during planetary system formation.

Inside the H$_{2}$O snowline of protoplanetary disks, water evaporates from the dust-grain surface into the gas phase, whereas it is frozen out on to the dust in the cold region beyond the snowline. H$_{2}$O ice enhances the solid material in the cold outer part of a disk, which promotes the formation of gas-giant planet cores. We can regard the H$_{2}$O snowline as the surface that divides the regions between rocky and gaseous giant planet formation. Thus observationally measuring the location of the H$_{2}$O snowline is crucial for understanding the planetesimal and planet formation processes, and the origin of water on Earth. In this paper, we find candidate water lines to locate the H$_{2}$O snowline through future high-dispersion spectroscopic observations. First, we calculate the chemical composition of the disk and investigate the abundance distributions of H$_{2}$O gas and ice, and the position of the H$_{2}$O snowline. We confirm that the abundance of H$_{2}$O gas is high not only in the hot midplane region inside the H$_{2}$O snowline but also in the hot surface layer of the outer disk. Second, we calculate the H$_{2}$O line profiles and identify those H$_{2}$O lines which are promising for locating the H$_{2}$O snowline: the identified lines are those which have small Einstein $A$ coefficients and high upper state energies. The wavelengths of the candidate H$_{2}$O lines range from mid-infrared to sub-millimeter, and they overlap with the regions accessible to ALMA and future mid-infrared high dispersion spectrographs (e.g., TMT/MICHI, SPICA).

Protoplanetary disks with non-axisymmetric structures have been observed. The Rossby wave instability (RWI) is considered as one of the origins of the non-axisymmetric structures. We perform linear stability analyses of the RWI in barotropic flow using four representative types of the background flow on a wide parameter space. We find that the co-rotation radius is located at the background vortensity minimum with large concavity if the system is marginally stable to the RWI, and this allows us to check the stability against the RWI easily. We newly derive the necessary and sufficient condition for the onset of the RWI in semi-analytic form. We discuss the applicability of the new condition in realistic systems and the physical nature of the RWI.

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.

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.

We investigate the dust and gas distribution in the disk around HD 142527 based on ALMA observations of dust continuum, 13CO(3-2), and C18O(3-2) emission. The disk shows strong azimuthal asymmetry in the dust continuum emission, while gas emission is more symmetric. In this paper, we investigate how gas and dust are distributed in the dust-bright northern part of the disk and in the dust-faint southern part. We construct two axisymmetric disk models. One reproduces the radial profiles of the continuum and the velocity moments 0 and 1 of CO lines in the north and the other reproduces those in the south. We have found that the dust is concentrated in a narrow ring having ~50AU width (in FWHM; w_d=30AU in our parameter definition) located at ~170-200AU from the central star. The dust particles are strongly concentrated in the north. We have found that the dust surface density contrast between the north and south amounts to ~70. Compared to the dust, the gas distribution is more extended in the radial direction. We find that the gas component extends at least from ~100AU to ~250AU from the central star, and there should also be tenuous gas remaining inside and outside of these radii. The azimuthal asymmetry of gas distribution is much smaller than dust. The gas surface density differs only by a factor of ~3-10 between the north and south. Hence, gas-to-dust ratio strongly depends on the location of the disk: ~30 at the location of the peak of dust distribution in the south and ~3 at the location of the peak of dust distribution in the north. Despite large uncertainties, the overall gas-to-dust ratio is inferred to be ~10-30, indicating that the gas depletion may have already been under way.

(Abridged) Near- to mid-IR observations of protoplanetary disks show that the inner regions (<10AU) are rich in small organic volatiles (e.g., C2H2 and HCN). Trends in the data suggest that disks around cooler stars (~3000K) are potentially more carbon- and molecule-rich than their hotter counterparts. Our aims are to explore the composition of the planet-forming region of disks around stars from M dwarf to Herbig Ae and compare with the observed trends. Models of the disk physical structure are coupled with a gas-grain chemical network to map the abundances in the planet-forming zone. N2 self shielding, X-ray-induced chemistry, and initial abundances, are investigated. The composition in the 'observable' atmosphere is compared with that in the midplane where the planet-building reservoir resides. M dwarf disk atmospheres are relatively more molecule rich than those for T Tauri or Herbig Ae disks. The weak far-UV flux helps retain this complexity which is enhanced by X-ray-induced ion-molecule chemistry. N2 self shielding has only a small effect and does not explain the higher C2H2/HCN ratios observed towards cooler stars. The models underproduce the OH/H2O column density ratios constrained in Herbig Ae disks, despite reproducing the absolute value for OH: H2O self shielding only increases this discrepency. The disk midplane content is sensitive to the initial main elemental reservoirs. The gas in the inner disk is generally more carbon rich than the midplane ices and is most significant for disks around cooler stars. The atmospheric C/O ratio appears larger than it actually is when calculated using observable tracers only because gas-phase O2 is predicted to be a significant oxygen reservoir. The models suggest that the gas in the inner regions of disks around cooler stars is more carbon rich; however, calculations of the molecular emission are necessary to confirm the observed trends.

We present ALMA Cycle 1 observations of the central kpc region of the luminous type-1 Seyfert galaxy NGC 7469 with unprecedented high resolution (0.5$"$ $\times$ 0.4$"$ = 165 pc $\times$ 132 pc) at submillimeter wavelengths. Utilizing the wide-bandwidth of ALMA, we simultaneously obtained HCN(4-3), HCO$^+$(4-3), CS(7-6), and partially CO(3-2) line maps, as well as the 860 $\mu$m continuum. The region consists of the central $\sim$ 1$"$ component and the surrounding starburst ring with a radius of $\sim$ 1.5$"$-2.5$"$. Several structures connect these components. Except for CO(3-2), these dense gas tracers are significantly concentrated towards the central $\sim$ 1$"$, suggesting their suitability to probe the nuclear regions of galaxies. Their spatial distribution resembles well those of centimeter and mid-infrared continuum emissions, but it is anti-correlated with the optical one, indicating the existence of dust obscured star formation. The integrated intensity ratios of HCN(4-3)/HCO$^+$(4-3) and HCN(4-3)/CS(7-6) are higher at the AGN position than at the starburst ring, which is consistent to our previous findings (submm-HCN enhancement). However, the HCN(4-3)/HCO$^+$(4-3) ratio at the AGN position of NGC 7469 (1.11$\pm$0.06) is almost half of the corresponding value of the low-luminosity type-1 Seyfert galaxy NGC 1097 (2.0$\pm$0.2), despite the more than two orders of magnitude higher X-ray luminosity of NGC 7469. But the ratio is comparable to that of the close vicinity of the AGN of NGC 1068 ($\sim$ 1.5). Based on these results, we speculate that some other heating mechanisms than X-ray (e.g., mechanical heating due to AGN jet) can contribute significantly for shaping the chemical composition in NGC 1097.

We investigate the chemistry of ion molecules in protoplanetary disks, motivated by the detection of N$_2$H$^+$ ring around TW Hya. While the ring inner radius coincides with the CO snow line, it is not apparent why N$_2$H$^+$ is abundant outside the CO snow line in spite of the similar sublimation temperatures of CO and N$_2$. Using the full gas-grain network model, we reproduced the N$_2$H$^+$ ring in a disk model with millimeter grains. The chemical conversion of CO and N$_2$ to less volatile species (sink effect hereinafter) is found to affect the N$_2$H$^+$ distribution. Since the efficiency of the sink depends on various parameters such as activation barriers of grain surface reactions, which are not well constrained, we also constructed the no-sink model; the total (gas and ice) CO and N$_2$ abundances are set constant, and their gaseous abundances are given by the balance between adsorption and desorption. Abundances of molecular ions in the no-sink model are calculated by analytical formulas, which are derived by analyzing the full-network model. The N$_2$H$^+$ ring is reproduced by the no-sink model, as well. The 2D (R-Z) distribution of N$_2$H$^+$, however, is different among the full-network model and no-sink model. The column density of N$_2$H$^+$ in the no-sink model depends sensitively on the desorption rate of CO and N$_2$, and the flux of cosmic ray. We also found that N$_2$H$^+$ abundance can peak at the temperature slightly below the CO sublimation, even if the desorption energies of CO and N$_2$ are the same.

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

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.

HD 100546 is a well-studied Herbig Be star-disk system that likely hosts a close-in companion with compelling observational evidence for an embedded protoplanet at 68 AU. We present ALMA observations of the HD 100546 disk which resolve the gas and dust structure at (sub)mm wavelengths. The CO emission (at 345.795 GHz) originates from an extensive molecular disk (390+/-20 AU in radius) whereas the continuum emission is more compact (230+/-20 AU in radius) suggesting radial drift of the mm-sized grains. The CO emission is similar in extent to scattered light images indicating well-mixed gas and um-sized grains in the disk atmosphere. Assuming azimuthal symmetry, a single-component power-law model cannot reproduce the continuum visibilities. The visibilities and images are better reproduced by a double-component model: a compact ring with a width of 21 AU centered at 26 AU and an outer ring with a width of 75+/-3 AU centered at 190+/-3 AU. The influence of a companion and protoplanet on the dust evolution is investigated. The companion at 10 AU facilitates the accumulation of mm-sized grains within a compact ring, ~ 20 - 30 AU, by ~ 10 Myr. The injection of a protoplanet at 1 Myr hastens the ring formation (~ 1.2 Myr) and also triggers the development of an outer ring (~ 100 - 200 AU). These observations provide additional evidence for the presence of a close-in companion and hint at dynamical clearing by a protoplanet in the outer disk.

We analytically calculate the marginally stable surface density profile for rotational instability of protoplanetary disks. The derived profile can be utilized for considering the region in a rotating disk where radial pressure gradient force is comparable to the gravitational force, such as an inner edge, steep gaps or bumps and an outer region of the disk. In this paper we especially focus on the rotational instability in the outer region of disks. We find an protoplanetary disk with a surface density profile of similarity solution becomes rotationally unstable at a certain radius, depending on its temperature profile and a mass of the central star. If the temperature is relatively low and the mass of the central star is high, disks have rotationally stable similarity profiles. Otherwise, deviation from the similarity profiles of surface density could be observable, using facilities with high sensitivity, such as ALMA.

(Abridged) The birth environment of the Sun will have influenced the conditions in the pre-solar nebula, including the attainable chemical complexity, important for prebiotic chemistry. The formation and distribution of complex organic molecules (COMs) in a disk around a T Tauri star is investigated for two scenarios: (i) an isolated disk, and (ii) a disk irradiated externally by a nearby massive star. The chemistry is calculated along the accretion flow from the outer disk inwards using a comprehensive network. Two simulations are performed, one beginning with complex ices and one with simple ices only. For the isolated disk, COMs are transported without major alteration into the inner disk where they thermally desorb into the gas reaching an abundance representative of the initial assumed ice abundance. For simple ices, COMs efficiently form on grain surfaces under the conditions in the outer disk. Gas-phase COMs are released into the molecular layer via photodesorption. For the irradiated disk, complex ices are also transported inwards; however, they undergo thermal processing caused by the warmer conditions in the irradiated disk which tends to reduce their abundance along the accretion flow. For simple ices, grain-surface chemistry cannot synthesise COMs in the outer disk because the necessary grain-surface radicals, which tend to be particularly volatile, are not sufficiently abundant on the grain surfaces. Gas-phase COMs are formed in the inner region of the irradiated disk via gas-phase chemistry induced by the desorption of strongly bound molecules such as methanol; hence, the abundances are not representative of the initial molecular abundances injected into the outer disk. These results suggest that the composition of comets formed in isolated disks may differ from those formed in externally irradiated disks with the latter composed of more simple ices.

(Abridged) Protoplanetary disks are vital objects in star and planet formation, possessing all the material which may form a planetary system orbiting the new star. We investigate the synthesis of complex organic molecules (COMs) in disks to constrain the achievable chemical complexity and predict species and transitions which may be observable with ALMA. We have coupled a 2D model of a protoplanetary disk around a T Tauri star with a gas-grain chemical network including COMs. We compare compare synthesised line intensities and calculated column densities with observations and determine those COMs which may be observable in future. COMs are efficiently formed in the disk midplane via grain-surface chemical reactions, reaching peak grain-surface fractional abundances 1e-6 - 1e-4 that of the H nuclei number density. COMs formed on grain surfaces are returned to the gas phase via non-thermal desorption; however, gas-phase species reach lower fractional abundances than their grain-surface equivalents, 1e-12 - 1e-7. Including the irradiation of grain mantle material helps build further complexity in the ice through the replenishment of grain-surface radicals which take part in further grain-surface reactions. There is reasonable agreement with several line transitions of H2CO observed towards several T Tauri star-disk systems. The synthesised line intensities for CH3OH are consistent with upper limits determined towards all sources. Our models suggest CH3OH should be readily observable in nearby protoplanetary disks with ALMA; however, detection of more complex species may prove challenging. Our grain-surface abundances are consistent with those derived from cometary comae observations providing additional evidence for the hypothesis that comets (and other planetesimals) formed via the coagulation of icy grains in the Sun's natal disk.

Protoplanetary disks composed of dust and gas are ubiquitous around young stars and are commonly recognized as nurseries of planetary systems. Their lifetime, appearance, and structure are determined by an interplay between stellar radiation, gravity, thermal pressure, magnetic field, gas viscosity, turbulence, and rotation. Molecules and dust serve as major heating and cooling agents in disks. Dust grains dominate the disk opacities, reprocess most of the stellar radiation, and shield molecules from ionizing UV/X-ray photons. Disks also dynamically evolve by building up planetary systems which drastically change their gas and dust density structures. Over the past decade significant progress has been achieved in our understanding of disk chemical composition thanks to the upgrade or advent of new millimeter/Infrared facilities (SMA, PdBI, CARMA, Herschel, e-VLA, ALMA). Some major breakthroughs in our comprehension of the disk physics and chemistry have been done since PPV. This review will present and discuss the impact of such improvements on our understanding of the disk physical structure and chemical composition.

We investigate water and deuterated water chemistry in turbulent protoplanetary disks. Chemical rate equations are solved with the diffusion term, mimicking turbulent mixing in vertical direction. Water near the midplane is transported to the disk atmosphere by turbulence and destroyed by photoreactions to produce atomic oxygen, while the atomic oxygen is transported to the midplane and reforms water and/or other molecules. We find that this cycle significantly decreases column densities of water ice at r < 30 AU, where dust temperatures are too high to reform water ice effectively. The radial extent of such region depends on the desorption energy of atomic hydrogen. Our model indicates that water ice could be deficient even outside the sublimation radius. Outside this radius, the cycle decreases the D/H ratio of water ice from 2x10^-2, which is set by the collapsing core model, to 10^-4-10^-2 in 10^6 yr, without significantly decreasing the water ice column density. The resultant D/H ratios depend on the strength of mixing and the radial distance from the central star. Our finding suggests that the D/H ratio of cometary water (10^-3-10^-4) could be established (i.e. cometary water could be formed) in the solar nebula, even if the D/H ratio of water ice delivered to the disk was very high (10^-2).

We report ALMA observations of dust continuum, 13CO J=3--2, and C18O J=3--2 line emission toward a gapped protoplanetary disk around HD 142527. The outer horseshoe-shaped disk shows the strong azimuthal asymmetry in dust continuum with the contrast of about 30 at 336 GHz between the northern peak and the southwestern minimum. In addition, the maximum brightness temperature of 24 K at its northern area is exceptionally high at 160 AU from a star. To evaluate the surface density in this region, the grain temperature needs to be constrained and was estimated from the optically thick 13CO J=3--2 emission. The lower limit of the peak surface density was then calculated to be 28 g cm-2 by assuming a canonical gas-to-dust mass ratio of 100. This finding implies that the region is locally too massive to withstand self-gravity since Toomre's Q <~1--2, and thus, it may collapse into a gaseous protoplanet. Another possibility is that the gas mass is low enough to be gravitationally stable and only dust grains are accumulated. In this case, lower gas-to-dust ratio by at least 1 order of magnitude is required, implying possible formation of a rocky planetary core.

GASPS is a far-infrared line and continuum survey of protoplanetary and young debris disks using PACS on the Herschel Space Observatory. The survey includes [OI] at 63 microns, as well as 70, 100 and 160um continuum, with the brightest objects also studied in [OI]145um, [CII]157um, H2O and CO. Targets included T Tauri stars and debris disks in 7 nearby young associations, and a sample of isolated Herbig AeBe stars. The aim was to study the global gas and dust content in a wide disk sample, systemically comparing the results with models. In this paper we review the main aims, target selection and observing strategy. We show initial results, including line identifications, sources detected, and a first statistical study. [OI]63um was the brightest line in most objects, by a factor of ~10. Detection rates were 49%, including 100% of HAeBe stars and 43% of T Tauri stars. Comparison with published dust masses show a dust threshold for [OI]63um detection of ~1e-5 M_solar. Normalising to 140pc distance, 32% with mass 1e-6 - 1e-5 M_solar, and a small number with lower mass were also detected. This is consistent with moderate UV excess and disk flaring. In most cases, continuum and line emission is spatially and spectrally unresolved, suggesting disk emission. ~10 objects were resolved, likely from outflows. Detection rates in [OI]145um, [CII]157um and CO J=18-17 were 20-40%, but [CII] was not correlated with disk mass, suggesting it arises instead from a compact envelope. [OI] detection rates in T Tauri associations of ages 0.3-4Myr were ~50%. ~2 stars were detectable in associations of 5-20Myr, with no detections in associations of age >20Myr. Comparing with the total number of young stars, and assuming a ISM-like gas/dust ratio, this indicates that ~18% of stars retain a gas-rich disk of total mass >1M_Jupiter for 1-4Myr, 1-7% keep such disks for 5-10Myr, and none remain beyond 10-20Myr.

Star formation often occurs within or nearby stellar clusters. Irradiation by nearby massive stars can photoevaporate protoplanetary disks around young stars (so-called proplyds) which raises questions regarding the ability of planet formation to take place in these environments. We investigate the two-dimensional physical and chemical structure of a protoplanetary disk surrounding a low-mass (T Tauri) star which is irradiated by a nearby massive O-type star to determine the survivability and observability of molecules in proplyds. Compared with an isolated star-disk system, the gas temperature ranges from a factor of a few (in the disk midplane) to around two orders of magnitude (in the disk surface) higher in the irradiated disk. Although the UV flux in the outer disk, in particular, is several orders of magnitude higher, the surface density of the disk is sufficient for effective shielding of the disk midplane so that the disk remains predominantly molecular in nature. We also find that non-volatile molecules, such as HCN and H2O, are able to freeze out onto dust grains in the disk midplane so that the formation of icy planetesimals, e.g., comets, may also be possible in proplyds. We have calculated the molecular line emission from the disk assuming LTE and determined that multiple transitions of atomic carbon, CO (and isotopologues, 13CO and C18O), HCO+, CN, and HCN may be observable with ALMA, allowing characterization of the gas column density, temperature, and optical depth in proplyds at the distance of Orion (~400 pc).

We investigate the impact of photochemistry and X-ray ionization on the molecular composition of, and ionization fraction in, a protoplanetary disk surrounding a typical T Tauri star. We use a sophisticated physical model, which includes a robust treatment of the radiative transfer of UV and X-ray radiation, and calculate the time-dependent chemical structure using a comprehensive chemical network. In previous work, we approximated the photochemistry and X-ray ionization, here, we recalculate the photoreaction rates using the explicit UV wavelength spectrum and wavelength-dependent reaction cross sections. We recalculate the X-ray ionization rate using our explicit elemental composition and X-ray energy spectrum. We find photochemistry has a larger influence on the molecular composition than X-ray ionization. Observable molecules sensitive to the photorates include OH, HCO+, N2H+, H2O, CO2 and CH3OH. The only molecule significantly affected by the X-ray ionization is N2H+ indicating it is safe to adopt existing approximations of the X-ray ionization rate in typical T Tauri star-disk systems. The recalculation of the photorates increases the abundances of neutral molecules in the outer disk, highlighting the importance of taking into account the shape of the UV spectrum in protoplanetary disks. A recalculation of the photoreaction rates also affects the gas-phase chemistry due to the adjustment of the H/H2 and C+/C ratios. The disk ionization fraction is not significantly affected by the methods adopted to calculate the photochemistry and X-ray ionization. We determine there is a probable 'dead zone' where accretion is suppressed, present in a layer, Z/R <~ 0.1 - 0.2, in the disk midplane, within R \approx 200 AU.

We calculate the chemical evolution of protoplanetary disks considering radial viscous accretion, vertical turbulent mixing and vertical disk winds. We study the effects on the disk chemical structure when different models for the formation of molecular hydrogen on dust grains are adopted. Our gas-phase chemistry is extracted from the UMIST Database for Astrochemistry (Rate06) to which we have added detailed gas-grain interactions. We use our chemical model results to generate synthetic near- and mid-infrared LTE line emission spectra and compare these with recent Spitzer observations. Our results show that if H2 formation on warm grains is taken into consideration, the H2O and OH abundances in the disk surface increase significantly. We find the radial accretion flow strongly influences the molecular abundances, with those in the cold midplane layers particularly affected. On the other hand, we show that diffusive turbulent mixing affects the disk chemistry in the warm molecular layers, influencing the line emission from the disk and subsequently improving agreement with observations. We find that NH3, CH3OH, C2H2 and sulphur-containing species are greatly enhanced by the inclusion of turbulent mixing. We demonstrate that disk winds potentially affect the disk chemistry and the resulting molecular line emission in a similar manner to that found when mixing is included.

We have developed a high resolution combined physical and chemical model of a protoplanetary disk surrounding a typical T Tauri star. Our aims were to use our model to calculate the chemical structure of disks on small scales (sub-milli-arcsecond in the inner disk for objects at the distance of Taurus, ~ 140 pc) to investigate the various chemical processes thought to be important in disks and to determine potential molecular tracers of each process. Our gas-phase network was extracted from the UMIST Database for Astrochemistry to which we added gas-grain interactions including freeze out and thermal and non-thermal desorption (cosmic-ray induced desorption, photodesorption and X-ray desorption) and a grain-surface network. We find that cosmic-ray induced desorption has the least effect on our disk chemical structure while photodesorption has a significant effect, enhancing the abundances of most gas-phase molecules throughout the disk and affecting the abundances and distribution of HCN, CN and CS, in particular. In the outer disk, we also see enhancements in the abundances of H2O and CO2. X-ray desorption is a potentially powerful mechanism in disks, acting to homogenise the fractional abundances of gas-phase species across the depth and increasing the column densities of most molecules although there remain significant uncertainties in the rates adopted for this process. The addition of grain-surface chemistry enhances the fractional abundances of several small complex organic molecules including CH3OH, HCOOCH3 and CH3OCH3 to potentially observable values (i.e. a fractional abundance of >~ 1.0E-11).

Planets are formed in disks around young stars. With an age of ~10 Myr, TW Hya is one of the nearest T Tauri stars that is still surrounded by a relatively massive disk. In addition a large number of molecules has been found in the TW Hya disk, making TW Hya the perfect test case in a large survey of disks with Herschel-PACS to directly study their gaseous component. We aim to constrain the gas and dust mass of the circumstellar disk around TW Hya. We observed the fine-structure lines of [OI] and [CII] as part of the Open-time large program GASPS. We complement this with continuum data and ground-based 12CO 3-2 and 13CO 3-2 observations. We simultaneously model the continuum and the line fluxes with the 3D Monte-Carlo code MCFOST and the thermo-chemical code ProDiMo to derive the gas and dust masses. We detect the [OI] line at 63 micron. The other lines that were observed, [OI] at 145 micron and [CII] at 157 micron, are not detected. No extended emission has been found. Preliminary modeling of the photometric and line data assuming [12CO]/[13CO]=69 suggests a dust mass for grains with radius < 1 mm of ~1.9 times 10^-4 Msun (total solid mass of 3 times 10^-3 Msun) and a gas mass of (0.5--5) times 10^-3 Msun. The gas-to-dust mass may be lower than the standard interstellar value of 100.

Context - Circumstellar discs are ubiquitous around young stars, but rapidly dissipate their gas and dust on timescales of a few Myr. The Herschel space observatory allows for the study of the warm disc atmosphere, using far-infrared spectroscopy to measure gas content and excitation conditions, and far-IR photometry to constrain the dust distribution. Aims - We aim to detect and characterize the gas content of circumstellar discs in four targets as part of the Herschel science demonstration phase. Methods - We carried out sensitive medium resolution spectroscopy and high sensitivity photometry at lambda ~60-190 micron using the Photodetector Array Camera and Spectrometer instrument on the Herschel space observatory. Results - We detect [OI] 63 micron emission from the young stars HD 169142, TW Hydrae, and RECX 15, but not HD 181327. No other lines, including [CII] 158 and [OI] 145, are significantly detected. All four stars are detected in photometry at 70 and 160 micron. Extensive models are presented in associated papers.

The Herschel GASPS Key Program is a survey of the gas phase of protoplanetary discs, targeting 240 objects which cover a large range of ages, spectral types, and disc properties. To interpret this large quantity of data and initiate self-consistent analyses of the gas and dust properties of protoplanetary discs, we have combined the capabilities of the radiative transfer code MCFOST with the gas thermal balance and chemistry code ProDiMo to compute a grid of 300 000 disc models (DENT). We present a comparison of the first Herschel/GASPS line and continuum data with the predictions from the DENT grid of models. Our objective is to test some of the main trends already identified in the DENT grid, as well as to define better empirical diagnostics to estimate the total gas mass of protoplanetary discs. Photospheric UV radiation appears to be the dominant gas-heating mechanism for Herbig stars, whereas UV excess and/or X-rays emission dominates for T Tauri stars. The DENT grid reveals the complexity in the analysis of far-IR lines and the difficulty to invert these observations into physical quantities. The combination of Herschel line observations with continuum data and/or with rotational lines in the (sub-)millimetre regime, in particular CO lines, is required for a detailed characterisation of the physical and chemical properties of circumstellar discs.

In an effort to simultaneously study the gas and dust components of the disc surrounding the young Herbig Ae star HD 169142, we present far-IR observations obtained with the PACS instrument onboard the Herschel Space Observatory. This work is part of the Open Time Key Project GASPS, which is aimed at studying the evolution of protoplanetary discs. To constrain the gas properties in the outer disc, we observed the star at several key gas-lines, including [OI] 63.2 and 145.5 micron, [CII] 157.7 micron, CO 72.8 and 90.2 micron, and o-H2O 78.7 and 179.5 micron. We only detect the [OI] 63.2 micron line in our spectra, and derive upper limits for the other lines. We complement our data set with PACS photometry and 12/13CO data obtained with the Submillimeter Array. Furthermore, we derive accurate stellar parameters from optical spectra and UV to mm photometry. We model the dust continuum with the 3D radiative transfer code MCFOST and use this model as an input to analyse the gas lines with the thermo-chemical code ProDiMo. Our dataset is consistent with a simple model in which the gas and dust are well-mixed in a disc with a continuous structure between 20 and 200 AU, but this is not a unique solution. Our modelling effort allows us to constrain the gas-to-dust mass ratio as well as the relative abundance of the PAHs in the disc by simultaneously fitting the lines of several species that originate in different regions. Our results are inconsistent with a gas-poor disc with a large UV excess; a gas mass of 5.0 +/- 2.0 times 10^(-3) Msun is still present in this disc, in agreement with earlier CO observations.

Despite a surface dominated by carbon-based life, the bulk composition of the Earth is dramatically carbon poor when compared to the material available at formation. Bulk carbon deficiency extends into the asteroid belt representing a fossil record of the conditions under which planets are born. The initial steps of planet formation involve the growth of primitive sub-micron silicate and carbon grains in the Solar Nebula. We present a solution wherein primordial carbon grains are preferentially destroyed by oxygen atoms ignited by heating due to stellar accretion at radii < 5 AU. This solution can account for the bulk carbon deficiency in the Earth and meteorites, the compositional gradient within the asteroid belt, and for growing evidence for similar carbon deficiency in rocks surrounding other stars.

We investigate the damping of the baryon acoustic oscillations in the matter power spectrum due to the quasinonlinear clustering and redshift-space distortions by confronting the models with the observations of the Sloan Digital Sky Survey luminous red galaxy sample. The chi-squared test suggests that the observed power spectrum is better matched by models with the damping of the baryon acoustic oscillations rather than the ones without the damping.

We have studied the dependence of the profiles of molecular abundances and line emission on the accretion flow in the hot ($\ga 100$K) inner region of protoplanetary disks. The gas-phase reactions initiated by evaporation of the ice mantle on dust grains are calculated along the accretion flow. We focus on methanol, a molecule that is formed predominantly through the evaporation of warm ice mantles, to show how the abundance profile and line emission depend on the accretion flow. Our results show that some evaporated molecules keep high abundances only when the accretion velocity is large enough, and that methanol could be useful as a diagnostic of the accretion flow by means of ALMA observations at the disk radius of $\la 10$AU.

We investigate the damping of the baryon acoustic oscillations (BAO) signature in the matter power spectrum due to the quasi-nonlinear clustering of density perturbations. On the basis of the third order perturbation theory, we construct a fitting formula of the damping in an analytic way. This demonstrates that the damping is closely related with the growth factor and the amplitude of the matter power spectrum. Then, we investigate the feasibility of constraining the growth factor through a measurement of the damping of the BAO signature. An extension of our formula including higher order corrections of density perturbations is also discussed.

An acoustic oscillation of the primeval photon-baryon fluid around the decoupling time imprints a characteristic scale in the galaxy distribution today, known as the baryon acoustic oscillation (BAO) scale. Several on-going and/or future galaxy surveys aim at detecting and precisely determining the BAO scale so as to trace the expansion history of the universe. We consider nonlinear and redshift-space distortion effects on the shifts of the BAO scale in $k$-space using perturbation theory. The resulting shifts are indeed sensitive to different choices of the definition of the BAO scale, which needs to be kept in mind in the data analysis. We present a toy model to explain the physical behavior of the shifts. We find that the BAO scale defined as in Percival et al. (2007) indeed shows very small shifts ($\lesssim$ 1%) relative to the prediction in {\it linear theory} in real space. The shifts can be predicted accurately for scales where the perturbation theory is reliable.

Detailed models for the density and temperature profiles of gas and dust in protoplanetary disks are constructed by taking into account X-ray and ultraviolet (UV) irradiation from a central T Tauri star, as well as dust size growth and settling toward the disk midplane. The spatial and size distributions of dust grains in the disks are numerically computed by solving the coagulation equation for settling dust particles. The level populations and line emission of molecular hydrogen are calculated using the derived physical structure of the disks. X-ray irradiation is the dominant heating source of the gas in the inner disk region and in the surface layer, while the far UV heating dominates otherwise. If the central star has strong X-ray and weak UV radiation, the H2 level populations are controlled by X-ray pumping, and the X-ray induced transition lines could be observable. If the UV irradiation is strong, the level populations are controlled by thermal collisions or UV pumping, depending on the properties of the dust grains in the disks. As the dust particles evolve in the disks, the gas temperature at the disk surface drops because the grain photoelectric heating becomes less efficient, while the UV radiation fields become stronger due to the decrease of grain opacity. This makes the H2 level populations change from local thermodynamic equilibrium (LTE) to non-LTE distributions, which results in changes to the line ratios of H2 emission. Our results suggest that dust evolution in protoplanetary disks could be observable through the H2 line ratios. The emission lines are strong from disks irradiated by strong UV and X-rays and possessing small dust grains; such disks will be good targets in which to observe H2 emission.

We calculate the physical structure of protoplanetary disks by evaluating the gas density and temperature self-consistently and solving separately for the dust temperature. The effect of grain growth is taken into account by assuming a power-law size distribution and varying the maximum radius of grains a_max. In our fiducial model with a_max=10um, the gas is warmer than the dust in the surface layer of the disk, while the gas and dust have the same temperature in deeper layers. In the models with larger a_max, the gas temperature in the surface layer is lower than in the fiducial model because of reduced photo-electric heating rates from small grains, while the deeper penetration of stellar radiation warms the gas at intermediate height. A detailed chemical reaction network is solved at outer radii (r \ge 50 AU). Vertical distributions of some molecular species at different radii are similar, when plotted as a function of hydrogen column density Sigma_H from the disk surface. Consequently, molecular column densities do not much depend on disk radius. In the models with larger a_max, the lower temperature in the surface layer makes the geometrical thickness of the disk smaller, and the gaseous molecules are confined to smaller heights. However, if we plot the vertical distributions of molecules as a function of Sigma_H, they do not significantly depend on a_max. The dependence of the molecular column densities on a_max is not significant, either. Notable exceptions are HCO+, H3+ and H2D+, which have smaller column densities in the models with larger a_max.

We have studied dust evolution in a quiescent or turbulent protoplanetary disk by numerically solving coagulation equation for settling dust particles, using the minimum mass solar nebular model. As a result, if we assume an ideally quiescent disk, the dust particles settle toward the disk midplane to form a gravitationally unstable layer within 2x10^3 - 4x10^4 yr at 1 - 30 AU, which is in good agreement with an analytic calculation by Nakagawa, Sekiya, & Hayashi (1986) although they did not take into account the particle size distribution explicitly. In an opposite extreme case of a globally turbulent disk, on the other hand, the dust particles fluctuate owing to turbulent motion of the gas and most particles become large enough to move inward very rapidly within 70 - 3x10^4 yr at 1 - 30 AU, depending on the strength of turbulence. Our result suggests that global turbulent motion should cease for the planetesimal formation in protoplanetary disks.

We have modeled self-consistently the density and temperature profiles of gas and dust in protoplanetary disks, taking into account irradiation from a central star. Making use of this physical structure, we have calculated the level populations of molecular hydrogen and the line emission from the disks. As a result, we can reproduce the observed strong line spectra of molecular hydrogen from protoplanetary disks, both in the ultraviolet (UV) and the near-infrared, but only if the central star has a strong UV excess radiation.

We have made self-consistent models of the density and temperature profiles of the gas and dust surrounding embedded luminous objects using a detailed radiative transfer model together with observations of the spectral energy distribution of hot molecular cores. Using these profiles we have investigated the hot core chemistry which results when grain mantles are evaporated, taking into account the different binding energies of the mantle molecules, as well a model in which we assume that all molecules are embedded in water ice and have a common binding energy. We find that most of the resulting column densities are consistent with those observed toward the hot core G34.3+0.15 at a time around 10$^4$ years after central luminous star formation. We have also investigated the dependence of the chemical structure on the density profile which suggests an observational possibility of constraining density profiles from determination of the source sizes of line emission from desorbed molecules.

We investigate the structure and the stabilities of a protoplanetary disk, which is heated by viscous process in itself and by its central star. The disk is set to rotate with the Keplerian velocity and has the surface density distribution of the minimum mass solar nebula. We assume the vertical hydrostatic equilibrium and the radiative equilibrium at each point, and solve the two-dimensional radiative transfer equation by means of the Short Characteristic method in the spherical coordinate in order to determine the disk structure. Our calculation shows that at the outer region of the disk with a distance from the central star of x>1AU the radiative heating from the inner disk dominates the viscous heating even near the midplane. It is because of the high temperature distribution in the optically thin surface layer and the relatively high disk height as a consequence of the irradiation from the inner hot region of the disk. In addition, we examine the convective and the magnetorotational instabilities of the disk. As a result, the whole disk is convectively stable since the dusty region is not heated by the viscous dissipation from the midplane but by the radial radiative heating. On the other hand, almost all the disk is magnetorotationally unstable except for the region near the equatorial plane of from 2AU to 10AU. Finally we discuss the growth and the size distribution of dust particles in the disk, which suggests that there exist cm-sized particles in the surface layer, namely, in the exposed region of the disk.

We consider the propagation of hydromagnetic waves generated by a compact turbulent source in low ionized plasmas, applying the Lighthill theory. We assume the plasma to be isothermal, and adopt a uniform, stationary medium thread by ordered magnetic fields as an initial condition. Then, the distinct properties of the hydromagnetic waves originating from a source oscillating with a fixed frequency are studied in the linear regime. As is well known, in low ionized plasmas, the generated waves dissipate due to ion-neutral damping. In this paper, the dependence of the dissipation rate on the frequency of the oscillating source is investigated. The larger the frequency becomes, the more substantial is the wave dissipation. Implications of our results on the energy source in molecular clouds are also discussed. Interestingly, since the outflow lobes associated with young stellar objects act as compact turbulent sources, hydromagnetic waves are generated by them. From our order-estimations, about 70% of the energy of the outflow itself propagates as waves or turbulences, while the remaining 30% dissipates and heats the neutrals via ion-neutral damping. Then, we confirm that the outflows are significant energy sources in molecular clouds in the context of the Lighthill theory.