The dispersal phase of planet-forming disks via winds driven by irradiation from the central star and/or magnetic fields in the disk itself is likely to play an important role in the formation and evolution of planetary systems. Current theoretical models lack predictive power to adequately constrain observations. We present PRIZMO, a code for evolving thermochemistry in protoplanetary disks capable of being coupled with hydrodynamical and multi-frequency radiative transfer codes. We describe the main features of the code, including gas and surface chemistry, photochemistry, microphysics, and the main cooling and heating processes. The results of a suite of benchmarks, which include photon-dominated regions, slabs illuminated by radiation spectra that include X-ray, and well-established cooling functions evaluated at different temperatures show good agreement both in terms of chemical and thermal structures. The development of this code is an important step to perform quantitative spectroscopy of disk winds, and ultimately the calculation of line profiles, which is urgently needed to shed light on the nature of observed disk winds.

ALMA surveys of nearby star-forming regions have shown that the dust mass in the disk is correlated with the stellar mass, but with a large scatter. This scatter could indicate either different evolutionary paths of disks or different initial conditions within a single cluster. We present ALMA Cycle 3 follow-up observations for 14 Class II disks that were low S/N detections or non-detections in our Cycle 2 survey of the $\sim 2$ Myr-old Chamaeleon I star-forming region. With 5 times better sensitivity, we detect millimeter dust continuum emission from six more sources and increase the detection rate to 94\% (51/54) for Chamaeleon I disks around stars earlier than M3. The stellar-disk mass scaling relation reported in \citet{pascucci2016} is confirmed with these updated measurements. Faint outliers in the $F_{mm}$--$M_*$ plane include three non-detections (CHXR71, CHXR30A, and T54) with dust mass upper limits of 0.2 M$_\oplus$ and three very faint disks (CHXR20, ISO91, and T51) with dust masses $\sim 0.5$ M$_\oplus$. By investigating the SED morphology, accretion property and stellar multiplicity, we suggest for the three millimeter non-detections that tidal interaction by a close companion ($<$100 AU) and internal photoevaporation may play a role in hastening the overall disk evolution. The presence of a disk around only the secondary star in a binary system may explain the observed stellar SEDs and low disk masses for some systems.

We present numerical simulations of molecular clouds (MCs) with self-consistent CO gas-phase and isotope chemistry in various environments. The simulations are post-processed with a line radiative transfer code to obtain $^{12}$CO and $^{13}$CO emission maps for the $J=1\rightarrow0$ rotational transition. The emission maps are analysed with commonly used observational methods, i.e. the $^{13}$CO column density measurement, the virial mass estimate and the so-called $X_{\textrm{CO}}$ (also CO-to-H$_2$) conversion factor, and then the inferred quantities (i.e. mass and column density) are compared to the physical values. We generally find that most methods examined here recover the CO-emitting H$_{2}$ gas mass of MCs within a factor of two uncertainty if the metallicity is not too low. The exception is the $^{13}$CO column density method. It is affected by chemical and optical depth issues, and it measures both the true H$_{2}$ column density distribution and the molecular mass poorly. The virial mass estimate seems to work the best in the considered metallicity and radiation field strength range, even when the overall virial parameter of the cloud is above the equilibrium value. This is explained by a systematically lower virial parameter (i.e. closer to equilibrium) in the CO-emitting regions; in CO emission, clouds might seem (sub-)virial, even when, in fact, they are expanding or being dispersed. A single CO-to-H$_{2}$ conversion factor appears to be a robust choice over relatively wide ranges of cloud conditions, unless the metallicity is low. The methods which try to take the metallicity dependence of the conversion factor into account tend to systematically overestimate the true cloud masses.

EX Lupi-type objects (EXors) form a sub-class of T Tauri stars, defined by sudden sporadic flare-ups of 1-5 magnitudes at optical wavelengths. These eruptions are attributed to enhanced mass accretion from the circumstellar disk to the star, and may constitute important events in shaping the structure of the inner disk and the forming planetary system. Although disk properties must play a fundamental role in driving the outbursts, they are surprisingly poorly known. In order to characterize the dust and gas components of EXor disks, here we report on observations of the $^{12}$CO J=3-2 and 4-3 lines, and the $^{13}$CO 3-2 line in EX Lup, the prototype of the EXor class. We reproduce the observed line fluxes and profiles with a line radiative transfer model, and compare the obtained parameters with corresponding ones of other T Tauri disks.

The $^{13}$CO molecule is often used as a column density tracer in regions where the $^{12}$CO emission saturates. The $^{13}$CO column density is then related to that of $^{12}$CO by a uniform isotopic ratio. A similar approximation is frequently used when deriving $^{13}$CO emission maps from numerical simulations of molecular clouds. To test this assumption we calculate the $^{12}$CO/$^{13}$CO ratio self-consistently, taking the isotope selective photodissociation and the chemical fractionation of CO into account. We model the coupled chemical, thermal and dynamical evolution and the emergent $^{13}$CO emission of isolated, starless molecular clouds in various environments. Selective photodissociation has a minimal effect on the ratio, while the chemical fractionation causes a factor of 2-3 decrease at intermediate cloud depths. The variation correlates with both the $^{12}$CO and the $^{13}$CO column densities. Neglecting the depth dependence results in $\leq$60 per cent error in $^{12}$CO column densities derived from $^{13}$CO. The same assumption causes $\leq$50 per cent disparity in the $^{13}$CO emission derived from simulated clouds. We show that the discrepancies can be corrected by a fitting formula. The formula is consistent with millimetre-wavelength isotopic ratio measurements of dense molecular clouds, but underestimates the ratios from the ultraviolet absorption of diffuse regions.

Circumstellar disks are expected to be the birthplaces of planets. The potential for forming one or more planets of various masses is essentially driven by the initial mass of the disks. We present and analyze Herschel/PACS observations of disk-bearing M-type stars that belong to the young ~2 Myr old Chamaleon-I star forming region. We used the radiative transfer code RADMC to successfully model the SED of 17 M-type stars detected at PACS wavelengths. We first discuss the relatively low detection rates of M5 and later spectral type stars with respect to the PACS sensitivity, and argue their disks masses, or flaring indices, are likely to be low. For M0 to M3 stars, we find a relatively broad range of disk masses, scale heights, and flaring indices. Via a parametrization of dust stratification, we can reproduce the peak fluxes of the 10 $\mu$m emission feature observed with Spitzer/IRS, and find that disks around M-type stars may display signs of dust sedimentation. The Herschel/PACS observations of low-mass stars in Cha-I provide new constraints on their disk properties, overall suggesting that disk parameters for early M-type stars are comparable to those for more massive stars (e.g., comparable scale height and flaring angles). However, regions of the disks emitting at about 100 $\mu$m may still be in the optically thick regime, preventing direct determination of disk masses. Thus the modeled disk masses should be considered as lower limits. Still, we are able to extend the wavelength coverage of SED models and start characterizing effects such as dust sedimentation, an effort leading the way towards ALMA observations of these low-mass stars.

(Abridged) Since the discovery of the first directly-imaged, planetary-mass object, 2MASS 1207 b, several works have sought to explain a disparity between its observed temperature and luminosity. Given its known age, distance, and spectral type, 2MASS 1207 b is under-luminous by a factor of ~10 (~2.5 mags) when compared to standard models of brown-dwarf/giant-planet evolution. In this paper, we study three possible sources of 2MASS 1207 b's under-luminosity. First, we investigate Mohanty et al. (2007)'s hypothesis that a near edge-on disk might be responsible for 2MASS 1207 b's under-luminosity. We conclude that the hypothesis is unlikely due to the lack of variability seen in multi-epoch photometry and unnecessary due to the increasing sample of under-luminous brown-dwarfs/giant-exoplanets that cannot be explained by an edge-on disk. Next, we test the analogous possibility that a spherical shell of dust, could explain 2MASS 1207 b's under-luminosity. Models containing enough dust to create ~2.5 mags of extinction, placed at reasonable radii, are ruled out by our new Gemini/T-ReCS 8.7 micron photometric upper-limit for 2MASS 1207 b. Finally, we investigate the possibility that 2MASS 1207 b is intrinsically cooler than the commonly used AMES-DUSTY fits to its spectrum, and thus it is not, in fact, under-luminous. New, thick cloud model grids by Madhusudhan et al. (2011) fit 2MASS 1207 b's 1-10 micron SED well, but they do not quite fit its near-infrared spectrum. However, we suggest that with some "tuning", they might be capable of simultaneously reproducing 2MASS 1207 b's spectral shape and luminosity. In this case, the whole class of young, under-luminous brown-dwarfs/giant-exoplanets might be explained by atmospheres that are able to suspend thick, dusty clouds in their photospheres at lower temperatures than field brown-dwarfs.

Previous studies suggest that the planet-forming disks around very-low-mass stars/brown dwarfs may be flatter than those around more massive stars, in contrast to model predictions of larger scale heights for gas-disks around lower-mass stars. We conducted a statistically robust study to determine whether there is evidence for stellar-mass-dependent disk structure in planet-forming disks. We find a statistically significant difference in the Spitzer/IRAC color distributions of disks around very-low-mass and low-mass stars all belonging to the same star-forming region, the Chamaeleon I star-forming region. We show that self consistently calculated disk models cannot fit the median spectral energy distributions (SEDs) of the two groups. These SEDs can be only explained by flatter disk models, consistent with the effect of dust settling in disks. We find that relative to the disk structure predicted for flared disks the required reduction in disk scale height is anti-correlated with the stellar mass, i.e. disks around lower-mass stars are flatter. Our results show that the initial and boundary conditions of planet formation are stellar-mass-dependent, an important finding that must be considered in planet formation models.