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.

The detection of gas in debris disks raises the question of whether this gas is a remnant from the primordial protoplanetary phase, or released by the collision of secondary bodies. In this paper we analyze ALMA observations at 1-1.5" resolution of three debris disks where the $^{12}$CO(2-1) rotational line was detected: HD131835, HD138813, and HD156623. We apply the iterative Lucy-Richardson deconvolution technique to the problem of circumstellar disks to derive disk geometries and surface brightness distributions of the gas. The derived disk parameters are used as input for thermochemical models to test both primordial and cometary scenarios for the origin of the gas. We favor a secondary origin for the gas in these disks and find that the CO gas masses ($\sim 3\times10^{-3}$ M$_{\oplus}$) require production rates ($\sim 5\times 10^{-7}$ M$_{\oplus}$~yr$^{-1}$) similar to those estimated for the bona-fide gas rich debris disk $\beta$ Pic.

We present a survey of optical [OI] emission at 6300 \AA\ toward 65 TTauri stars at the spectral resolution of $\sim7$ km/s. Past work identified a highly blueshifted velocity component (HVC) tracing microjets, and a less blueshifted low-velocity-component (LVC) attributed to winds. We focus here on the LVC kinematics to investigate links between winds, jets, accretion, and disk dispersal. We track the behavior of four types of LVC components: a broad and narrow component ("BC" and "NC") in LVCs that are decomposed into two Gaussians, which typically have an HVC, and the single-Gaussian LVC profiles separated into those that have an HVC ("SCJ") and those that do not ("SC"). LVC centroid velocities and line widths correlate with HVC equivalent width and accretion luminosity, suggesting that LVC/winds and HVC/jets are kinematically linked and connected to accretion. The deprojected HVC velocity correlates with accretion luminosity, showing that faster jets come with higher accretion. BC and NC kinematics correlate and their blueshifts are maximum at $\sim 35^{\circ}$, suggesting a conical wind geometry with this semi-opening angle. Only SCs include $n_{13-31}$ up to $\sim3$ and their properties correlate with this infrared index, showing that [OI] emission recedes to larger radii as the inner dust is depleted, tracing less dense/hot gas and a decrease in wind velocity. All together, these findings support a scenario where optically thick, accreting inner disks launch radially-extended MHD disk winds that feed jets, and where inner disk winds recede to larger radii and jets disappear in concert with dust depletion.

Magnetohydrodynamic (MHD) and photoevaporative winds are thought to play an important role in the evolution and dispersal of planet-forming disks. We report the first high-resolution ($\Delta v\sim$6\kms) analysis of [S II] $\lambda$4068, [O I] $\lambda$5577, and [O I] $\lambda$6300 lines from a sample of 48 T Tauri stars. Following Simon et al. (2016), we decompose them into three kinematic components: a high-velocity component (HVC) associated with jets, and a low-velocity narrow (LVC-NC) and broad (LVC-BC) components. We confirm previous findings that many LVCs are blueshifted by more than 1.5 kms$^{-1}$ thus most likely trace a slow disk wind. We further show that the profiles of individual components are similar in the three lines. We find that most LVC-BC and NC line ratios are explained by thermally excited gas with temperatures between 5,000$-$10,000 K and electron densities $\sim10^{7}-10^{8}$ cm$^{-3}$. The HVC ratios are better reproduced by shock models with a pre-shock H number density of $\sim10^{6}-10^{7}$ cm$^{-3}$. Using these physical properties, we estimate $\dot{M}_{\rm wind}/\dot{M}_{\rm acc}$ for the LVC and $\dot{M}_{\rm jet}/\dot{M}_{\rm acc}$ for the HVC. In agreement with previous work, the mass carried out in jets is modest compared to the accretion rate. With the likely assumption that the NC wind height is larger than the BC, the LVC-BC $\dot{M}_{\rm wind}/\dot{M}_{\rm acc}$ is found to be higher than the LVC-NC. These results suggest that most of the mass loss occurs close to the central star, within a few au, through an MHD driven wind. Depending on the wind height, MHD winds might play a major role in the evolution of the disk mass.

Observations of multiple rotational transitions from a single molecule allow for unparalleled constraints on the physical conditions of the emitting region. We present an analysis of CS in TW~Hya using the $J=7-6$, $5-4$ and $3-2$ transitions imaged at $\sim 0.5^{\prime\prime}$ spatial resolution, resulting in a temperature and column density profile of the CS emission region extending out to 230~au, far beyond previous measurements. In addition, the 15~kHz resolution of the observations and the ability to directly estimate the temperature of the CS emitting gas, allow for one of the most sensitive searches for turbulent broadening in a disk to date. Limits of $v_{\rm turb} \lesssim 0.1 c_s$ can be placed across the entire radius of the disk. We are able to place strict limits of the local H$_2$ density due to the collisional excitations of the observed transitions. From these we find that a minimum disk mass of $3 \times 10^{-4}~M_{\rm sun}$ is required to be consistent with the CS excitation conditions and can uniquely constrain the gas surface density profile in the outer disk.

We present a suite of ALMA interferometric molecular line and continuum images that elucidate, on linear size scales of $\sim$30--40 AU, the chemical structure of the nearby, evolved, protoplanetary disk orbiting the close binary system V4046 Sgr. The observations were undertaken in the 1.1--1.4 mm wavelength range (ALMA Bands 6 and 7) with antenna configurations involving maximum baselines of several hundred meters, yielding subarcsecond-resolution images in more than a dozen molecular species and isotopologues. Isotopologues of CO and HCN display centrally peaked morphologies of integrated emission line intensity, whereas the line emission from complex nitrile group molecules (HC$_3$N, CH$_3$CN), deuterated molecules (DCN, DCO$^+$), hydrocarbons (as traced by C$_2$H), and potential CO ice line tracers (N$_2$H$^+$, and H$_2$CO) appears as a sequence of sharp and diffuse rings. The dimensions and morphologies of HC$_3$N and CH$_3$CN emission are suggestive of photodesorption of organic ices from the surfaces of dust grains, while the sequence of increasing radius of peak intensity represented by DCN (smallest), DCO$^+$, N$_2$H$^+$, and H$_2$CO (largest) is qualitatively consistent with the expected decline of midplane gas temperature with increasing disk radius. Empirical modeling indicates that the sharp-edged C$_2$H emission ring lies at relatively deep disk layers, leaving open the question of the origin of C$_2$H abundance enhancements in evolved disks. This study of the "molecular anatomy" of V4046 Sgr should serve as motivation for additional subarcsecond ALMA molecular line imaging surveys of nearby, evolved protoplanetary disks aimed at addressing major uncertainties in protoplanetary disk physical and chemical structure and molecular production pathways.

Determining the gas density and temperature structures of protoplanetary disks is a fundamental task to constrain planet formation theories. This is a challenging procedure and most determinations are based on model-dependent assumptions. We attempt a direct determination of the radial and vertical temperature structure of the Flying Saucer disk, thanks to its favorable inclination of 90 degrees. We present a method based on the tomographic study of an edge-on disk. Using ALMA, we observe at 0.5$"$ resolution the Flying Saucer in CO J=2-1 and CS J=5-4. This edge-on disk appears in silhouette against the CO J=2-1 emission from background molecular clouds in $\rho$ Oph. The combination of velocity gradients due to the Keplerian rotation of the disk and intensity variations in the CO background as a function of velocity provide a direct measure of the gas temperature as a function of radius and height above the disk mid-plane. The overall thermal structure is consistent with model predictions, with a cold ($< 15-12 $~K), CO-depleted mid-plane, and a warmer disk atmosphere. However, we find evidence for CO gas along the mid-plane beyond a radius of about 200\,au, coincident with a change of grain properties. Such a behavior is expected in case of efficient rise of UV penetration re-heating the disk and thus allowing CO thermal desorption or favoring direct CO photo-desorption. CO is also detected up to 3-4 scale heights while CS is confined around 1 scale height above the mid-plane. The limits of the method due to finite spatial and spectral resolutions are also discussed. This method appears to be very promising to determine the gas structure of planet-forming disks, provided that the molecular data have an angular resolution which is high enough, of the order of $0.3 - 0.1"$ at the distance of the nearest star forming regions.

Recent years have seen growing interest in the streaming instability as a candidate mechanism to produce planetesimals. However, these investigations have been limited to small-scale simulations. We now present the results of a global protoplanetary disk evolution model that incorporates planetesimal formation by the streaming instability, along with viscous accretion, photoevaporation by EUV, FUV, and X-ray photons, dust evolution, the water ice line, and stratified turbulence. Our simulations produce massive (60-130 $M_\oplus$) planetesimal belts beyond 100 au and up to $\sim 20 M_\oplus$ of planetesimals in the middle regions (3-100 au). Our most comprehensive model forms 8 $M_\oplus$ of planetesimals inside 3 au, where they can give rise to terrestrial planets. The planetesimal mass formed in the inner disk depends critically on the timing of the formation of an inner cavity in the disk by high-energy photons. Our results show that the combination of photoevaporation and the streaming instability are efficient at converting the solid component of protoplanetary disks into planetesimals. Our model, however, does not form enough early planetesimals in the inner and middle regions of the disk to give rise to giant planets and super-Earths with gaseous envelopes. Additional processes such as particle pileups and mass loss driven by MHD winds may be needed to drive the formation of early planetesimal generations in the planet forming regions of protoplanetary disks.

We present ALMA Cycle~2 observations at $0.5^{\prime\prime}$ resolution of TW Hya of CS $J=5-4$ emission. The radial profile of the integrated line emission displays oscillatory features outwards of $1.5^{\prime\prime}$ ($\approx 90$ au). A dip-like feature at $1.6^{\prime\prime}$ is coincident in location, depth and width with features observed in dust scattered light at near-infrared wavelengths. Using a thermochemical model indicative of TW Hya, gas-grain chemical modelling and non-LTE radiative transfer, we demonstrate that such a feature can be reproduced with a surface density depression, consistent with the modelling performed for scattered light observations of TW Hya. We further demonstrate that a gap in the dust distribution and dust opacity only cannot reproduce the observed CS feature. The outer enhancement at $3.1^{\prime\prime}$ is identified as a region of intensified desorption due to enhanced penetration of the interstellar FUV radiation at the exponential edge of the disk surface density, which intensifies the photochemical processing of gas and ices.

Using Keck/HIRES spectra {\Delta}v ~ 7 km/s, we analyze forbidden lines of [O I] 6300 {\AA}, [O I] 5577 {\AA} and [S II] 6731 {\AA} from 33 T Tauri stars covering a range of disk evolutionary stages. After removing a high velocity component (HVC) associated with microjets, we study the properties of the low velocity component (LVC). The LVC can be attributed to slow disk winds that could be magnetically (MHD) or thermally (photoevaporative) driven. Both of these winds play an important role in the evolution and dispersal of protoplanetary material. LVC emission is seen in all 30 stars with detected [O I] but only in 2 out of eight with detected [S II] , so our analysis is largely based on the properties of the [O I] LVC. The LVC itself is resolved into broad (BC) and narrow (NC) kinematic components. Both components are found over a wide range of accretion rates and their luminosity is correlated with the accretion luminosity, but the NC is proportionately stronger than the BC in transition disks. The FWHM of both the BC and NC correlates with disk inclination, consistent with Keplerian broadening from radii of 0.05 to 0.5 AU and 0.5 to 5 AU, respectively. The velocity centroids of the BC suggest formation in an MHD disk wind, with the largest blueshifts found in sources with closer to face-on orientations. The velocity centroids of the NC however, show no dependence on disk inclination. The origin of this component is less clear and the evidence for photoevaporation is not conclusive.

We present a CO(2-1) and 1240 um continuum survey of 23 debris disks with spectral types B9-G1, observed at an angular resolution of 0.5-1 arcsec with the Atacama Large Millimeter/Submillimeter Array (ALMA). The sample was selected for large infrared excess and age ~10 Myr, to characterize the prevalence of molecular gas emission in young debris disks. We identify three CO-rich debris disks, plus two additional tentative (3-sigma) CO detections. Twenty disks were detected in the continuum at the >3-sigma level. For the 12 disks in the sample that are spatially resolved by our observations, we perform an independent analysis of the interferometric continuum visibilities to constrain the basic dust disk geometry, as well as a simultaneous analysis of the visibilities and broad-band spectral energy distribution to constrain the characteristic grain size and disk mass. The gas-rich debris disks exhibit preferentially larger outer radii in their dust disks, and a higher prevalence of characteristic grain sizes smaller than the blowout size. The gas-rich disks do not exhibit preferentially larger dust masses, contrary to expectations for a scenario in which a higher cometary destruction rate would be expected to result in a larger mass of both CO and dust. The three debris disks in our sample with strong CO detections are all around A stars: the conditions in disks around intermediate-mass stars appear to be the most conducive to the survival or formation of CO.

We obtain high spatial and spectral resolution images of the CO J=2-1, CN N=2-1 and CS J=5-4 emission with ALMA in Cycle~2. The radial distribution of the turbulent broadening is derived with three approaches: two `direct' and one modelling. The first requires a single transition and derives \Tex{} directly from the line profile, yielding a \vturb{}. The second assumes two different molecules are co-spatial thus their relative linewidths allow for a calculation of \Tkin{} and \vturb{}. Finally we fit a parametric disk model where physical properties of the disk are described by power laws, to compare our `direct' methods with previous values. The two direct methods were limited to the outer $r > 40$~au disk due to beam smear. The direct method found \vturb{} ranging from $\approx$~\vel{130} at 40~au, dropping to $\approx$~\vel{50} in the outer disk, qualitatively recovered with the parametric model fitting. This corresponds to roughly $0.2 - 0.4~c_s$. CN was found to exhibit strong non-LTE effects outside $r \approx 140$~au, so \vturb{} was limited to within this radius. The assumption that CN and CS are co-spatial is consistent with observed linewidths only within $r \lesssim 100$~au, within which \vturb{} was found to drop from \vel{100} ($\approx~0.4~c_s$) to nothing at 100~au. The parametric model yielded a near constant \vel{50} for CS ($0.2 - 0.4~c_s$). We demonstrate that absolute flux calibration is and will be the limiting factor in all studies of turbulence using a single molecule. The magnitude of the dispersion is comparable with or below that predicted by the magneto-rotational instability theory. A more precise comparison would require to reach an absolute calibration precision of order 3\%, or to find a suitable combination of light and heavy molecules which are co-located in the disk.

To investigate the potential connection between the intense X-ray emission from young, low-mass stars and the lifetimes of their circumstellar, planet-forming disks, we have compiled the X-ray luminosities ($L_X$) of M stars in the $\sim$8 Myr-old TW Hya Association (TWA) for which X-ray data are presently available. Our investigation includes analysis of archival Chandra data for the TWA binary systems TWA 8, 9, and 13. Although our study suffers from poor statistics for stars later than M3, we find a trend of decreasing $L_X/L_{bol}$ with decreasing $T_{eff}$ for TWA M stars wherein the earliest-type (M0--M2) stars cluster near $\log{(L_X/L_{bol})} \approx -3.0$ and then $\log{(L_X/L_{bol})}$ decreases, and its distribution broadens, for types M4 and later. The fraction of TWA stars that display evidence for residual primordial disk material also sharply increases in this same (mid-M) spectral type regime. This apparent anticorrelation between the relative X-ray luminosities of low-mass TWA stars and the longevities of their circumstellar disks suggests that primordial disks orbiting early-type M stars in the TWA have dispersed rapidly as a consequence of their persistent large X-ray fluxes. Conversely, the disks orbiting the very lowest-mass pre-MS stars and pre-MS brown dwarfs in the Association may have survived because their X-ray luminosities and, hence, disk photoevaporation rates are very low to begin with, and then further decline relatively early in their pre-MS evolution.

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.

Protoplanetary disks dissipate rapidly after the central star forms, on time-scales comparable to those inferred for planet formation. In order to allow the formation of planets, disks must survive the dispersive effects of UV and X-ray photoevaporation for at least a few Myr. Viscous accretion depletes significant amounts of the mass in gas and solids, while photoevaporative flows driven by internal and external irradiation remove most of the gas. A reasonably large fraction of the mass in solids and some gas get incorporated into planets. Here, we review our current understanding of disk evolution and dispersal, and discuss how these might affect planet formation. We also discuss existing observational constraints on dispersal mechanisms and future directions.

We present a detailed analysis of narrow of NaI and KI absorption resonance lines toward nearly 40 T Tauri stars in Taurus with the goal of clarifying their origin. The NaI 5889.95 angstrom line is detected toward all but one source, while the weaker KI 7698.96 angstrom line in about two thirds of the sample. The similarity in their peak centroids and the significant positive correlation between their equivalent widths demonstrate that these transitions trace the same atomic gas. The absorption lines are present towards both disk and diskless young stellar objects, which excludes cold gas within the circumstellar disk as the absorbing material. A comparison of NaI and CO detections and peak centroids demonstrates that the atomic and molecular gas are not co-located, the atomic gas is more extended than the molecular gas. The width of the atomic lines corroborates this finding and points to atomic gas about an order of magnitude warmer than the molecular gas. The distribution of NaI radial velocities shows a clear spatial gradient along the length of the Taurus molecular cloud filaments. This suggests that absorption is associated with the Taurus molecular cloud. Assuming the gradient is due to cloud rotation, the rotation of the atomic gas is consistent with differential galactic rotation while the rotation of the molecular gas, although with the same rotation axis, is retrograde. Our analysis shows that narrow NaI and KI absorption resonance lines are useful tracers of the atomic envelope of molecular clouds. In line with recent findings from giant molecular clouds, our results demonstrate that the velocity fields of the atomic and molecular gas are misaligned. The angular momentum of a molecular cloud is not simply inherited from the rotating Galactic disk from which it formed but may be redistributed by cloud-cloud interactions.

We have used the Submillimeter Array to image, at ~1" resolution, C2H(3-2) emission from the molecule-rich circumstellar disks orbiting the nearby, classical T Tauri star systems TW Hya and V4046 Sgr. The SMA imaging reveals that the C2H emission exhibits a ring-like morphology within each disk, the inner hole radius of the C2H ring within the V4046 Sgr disk (~70 AU) is somewhat larger than than of its counterpart within the TW Hya disk (~45 AU). We suggest that, in each case, the C2H emission likely traces irradiation of the tenuous surface layers of the outer disks by high-energy photons from the central stars.

We have used the Submillimeter Array to image, at ~1.5" resolution, C2H (3-2) emission from the circumstellar disk orbiting the nearby (D = 54 pc), ~8 Myr-old, ~0.8 Msun classical T Tauri star TW Hya. The SMA imaging reveals that the C2H emission exhibits a ring-like morphology. Based on a model in which the C2H column density follows a truncated radial power-law distribution, we find that the inner edge of the ring lies at ~45 AU, and that the ring extends to at least ~120 AU. Comparison with previous (single-dish) observations of C2H (4-3) emission indicates that the C2H molecules are subthermally excited and, hence, that the emission arises from the relatively warm, tenuous upper atmosphere of the disk. We propose that the C2H emission most likely traces particularly efficient photo-destruction of small grains and/or photodesorption and photodissociation of hydrocarbons derived from grain ice mantles in the surface layers of the outer disk. The presence of a C2H ring in the TW Hya disk hence likely serves as a marker of dust grain processing and radial and vertical grain size segregation within the disk.

Protoplanetary disks are dispersed by viscous evolution and photoevaporation in a few million years; in the interim small, sub-micron sized dust grains must grow and form planets. The time-varying abundance of small grains in an evolving disk directly affects gas heating by far-ultraviolet photons, while dust evolution affects photoevaporation by changing the disk opacity and resulting penetration of FUV photons in the disk. Photoevaporative flows, in turn, selectively carry small dust grains leaving the larger particles---which decouple from the gas---behind in the disk. We study these effects by investigating the evolution of a disk subject to viscosity, photoevaporation by EUV, FUV and X-rays, dust evolution, and radial drift using a 1-D multi-fluid approach (gas + different dust grain sizes) to solve for the evolving surface density distributions. The 1-D evolution is augmented by 1+1D models constructed at each epoch to obtain the instantaneous disk structure and determine photoevaporation rates. The implementation of a dust coagulation/fragmentation model results in a marginal decrease in disk lifetimes when compared to models with no dust evolution; the disk lifetime is thus found to be relatively insensitive to the evolving dust opacity. We find that photoevaporation can cause significant reductions in the gas/dust mass ratio in the planet-forming regions of the disk as it evolves, and may result in a corresponding increase in heavy element abundances relative to hydrogen. We discuss implications for theories of planetesimal formation and giant planet formation, including the formation of gas-poor giants. After gas disk dispersal, $\sim 3\times 10^{-4}$ \ms\ of mass in solids typically remain, comparable to the solids inventory of our solar system.

In this paper we investigate the origin of the mid-infrared (IR) hydrogen recombination lines for a sample of 114 disks in different evolutionary stages (full, transitional and debris disks) collected from the {\it Spitzer} archive. We focus on the two brighter {H~{\sc i}} lines observed in the {\it Spitzer} spectra, the {H~{\sc i}}(7-6) at 12.37$\mu$m and the {H~{\sc i}}(9-7) at 11.32$\mu$m. We detect the {H~{\sc i}}(7-6) line in 46 objects, and the {H~{\sc i}}(9-7) in 11. We compare these lines with the other most common gas line detected in {\it Spitzer} spectra, the {[Ne~{\sc iii}]} at 12.81$\mu$m. We argue that it is unlikely that the {H~{\sc i}} emission originates from the photoevaporating upper surface layers of the disk, as has been found for the {[Ne~{\sc iii}]} lines toward low-accreting stars. Using the {H~{\sc i}}(9-7)/{H~{\sc i}}(7-6) line ratios we find these gas lines are likely probing gas with hydrogen column densities of 10$^{10}$-10$^{11}$~cm$^{-3}$. The subsample of objects surrounded by full and transitional disks show a positive correlation between the accretion luminosity and the {H~{\sc i}} line luminosity. These two results suggest that the observed mid-IR {H~{\sc i}} lines trace gas accreting onto the star in the same way as other hydrogen recombination lines at shorter wavelengths. A pure chromospheric origin of these lines can be excluded for the vast majority of full and transitional disks.We report for the first time the detection of the {H~{\sc i}}(7-6) line in eight young (< 20~Myr) debris disks. A pure chromospheric origin cannot be ruled out in these objects. If the {H~{\sc i}}(7-6) line traces accretion in these older systems, as in the case of full and transitional disks, the strength of the emission implies accretion rates lower than 10$^{-10}$M$_{\odot}$/yr. We discuss some advantages of extending accretion indicators to longer wavelengths.

The amount of high-energy stellar radiation reaching the surface of protoplanetary disks is essential to determine their chemistry and physical evolution. Here, we use millimetric and centimetric radio data to constrain the EUV luminosity impinging on 14 disks around young (~2-10Myr) sun-like stars. For each object we identify the long-wavelength emission in excess to the dust thermal emission, attribute that to free-free disk emission, and thereby compute an upper limit to the EUV reaching the disk. We find upper limits lower than 10$^{42}$ photons/s for all sources without jets and lower than $5 \times 10^{40}$ photons/s for the three older sources in our sample. These latter values are low for EUV-driven photoevaporation alone to clear out protoplanetary material in the timescale inferred by observations. In addition, our EUV upper limits are too low to reproduce the [NeII] 12.81 micron luminosities from three disks with slow [NeII]-detected winds. This indicates that the [NeII] line in these sources primarily traces a mostly neutral wind where Ne is ionized by 1 keV X-ray photons, implying higher photoevaporative mass loss rates than those predicted by EUV-driven models alone. In summary, our results suggest that high-energy stellar photons other than EUV may dominate the dispersal of protoplanetary disks around sun-like stars.

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.

The formation time, masses, and location of planets are strongly impacted by the physical mechanisms that disperse protoplanetary disks and the timescale over which protoplanetary material is cleared out. Accretion of matter onto the central star, protostellar winds/jets, magnetic disk winds, and photoevaporative winds operate concurrently. Hence, disentangling their relative contribution to disk dispersal requires identifying diagnostics that trace different star-disk environments. Here, we analyze the low velocity component (LVC) of the Oxygen optical forbidden lines, which is found to be blueshifted by a few km/s with respect to the stellar velocity. We find that the [OI] LVC profiles are different from those of [NeII] at 12.81\mu m and CO at 4.7\mu m lines pointing to different origins for these gas lines. We report a correlation between the luminosity of the [OI] LVC and the accretion luminosity Lacc. We do not find any correlation with the X-ray luminosity, while we find that the higher is the stellar FUV luminosity, the higher is the luminosity of the [OI] LVC. In addition, we show that the [OI]6300\AA/5577\AA\ ratio is low (ranging between 1 and 8). These findings favor an origin of the [OI] LVC in a region where OH is photodissociated by stellar FUV photons and argue against thermal emission from an X-ray-heated layer. Detailed modeling of two spectra with the highest S/N and resolution shows that there are two components within the LVC: a broad, centrally peaked component that can be attributed to gas arising in a warm disk surface in Keplerian rotation (with FWHM between ~40 and ~60 km/s), and a narrow component (with FWHM ~10 km/s and small blueshifts of ~2 km/s) that may arise in a cool (<1,000 K) molecular wind.

From the masses of planets orbiting our Sun, and relative elemental abundances, it is estimated that at birth our Solar System required a minimum disk mass of ~0.01 solar masses within ~100 AU of the star. The main constituent, gaseous molecular hydrogen, does not emit from the disk mass reservoir, so the most common measure of the disk mass is dust thermal emission and lines of gaseous carbon monoxide. Carbon monoxide emission generally probes the disk surface, while the conversion from dust emission to gas mass requires knowledge of the grain properties and gas-to-dust mass ratio, which likely differ from their interstellar values. Thus, mass estimates vary by orders of magnitude, as exemplified by the relatively old (3--10 Myr) star TW Hya, with estimates ranging from 0.0005 to 0.06 solar masses. Here we report the detection the fundamental rotational transition of hydrogen deuteride, HD, toward TW Hya. HD is a good tracer of disk gas because it follows the distribution of molecular hydrogen and its emission is sensitive to the total mass. The HD detection, combined with existing observations and detailed models, implies a disk mass >0.05 solar masses, enough to form a planetary system like our own.

We present Herschel PACS spectroscopy of the [OI] 63 micron gas-line for three circumstellar disk systems showing signs of significant disk evolution and/or planet formation: HR 8799, HD 377 and RX J1852.3-3700. [OI] is undetected toward HR 8799 and HD 377 with 3 sigma upper limits of 6.8 x 10^-18 W m^-2 and 9.9 x 10^-18 W m^-2 respectively. We find an [OI] detection for RX J1852.3-3700 at 12.3 +- 1.8 x 10^-18 W m^-2. We use thermo-chemical disk models to model the gas emission, using constraints on the [OI] 63 micron, and ancillary data to derive gas mass upper limits and constrain gas-to-dust ratios. For HD 377 and HR 8799, we find 3 sigma upper limits on the gas mass of 0.1-20 Mearth. For RX J1852.3-3700, we find two distinct disk scenarios that could explain the detection of [OI] 63 micron and CO(2-1) upper limits reported from the literature: (i) a large disk with gas co-located with the dust (16-500 AU), resulting in a large tenuous disk with ~16 Mearth of gas, or (ii) an optically thick gas disk, truncated at ~70 AU, with a gas mass of 150 Mearth. We discuss the implications of these results for the formation and evolution of planets in these three systems.

Recent infrared observations have demonstrated that photoevaporation driven by high-energy photons from the central star contributes to the dispersal of protoplanetary disks. Here, we show that photoevaporative winds should produce a detectable free-free continuum emission given the range of stellar ionizing photons and X-ray luminosities inferred for young sun-like stars. We point out that VLA observations of the nearby disk around TWHya might have already detected this emission at centimeter wavelengths and calculate the wind electron density and mass flow rate. We also estimate the intensities of H radio recombination lines tracing the wind and discuss which ones could be detected with current instrumentation. The detection and profiles of these recombination lines would unambiguously prove our inference of free-free emission from photoevaporating disks like TWHya. In addition, radio/millimeter data can help constraining wind parameters such as temperature and electron density that are fundamental in measuring mass flow rates.

Photoevaporation driven by the central star is expected to be a ubiquitous and important mechanism to disperse the circumstellar dust and gas from which planets form. Here, we present a detailed study of the circumstellar disk surrounding the nearby star TW Hya and provide observational constraints to its photoevaporative wind. Our new high-resolution (R ~ 30,000) mid-infrared spectroscopy in the [Ne II] 12.81 {\mu}m line confirms that this gas diagnostic traces the unbound wind component within 10AU from the star. From the blueshift and asymmetry in the line profile, we estimate that most (>80%) of the [Ne II] emission arises from disk radii where the midplane is optically thick to the redshifted outflowing gas, meaning beyond the 1 or 4AU dust rim inferred from other observations. We re-analyze high-resolution (R ~ 48, 000) archival optical spectra searching for additional transitions that may trace the photoevaporative flow. Unlike the [Ne II] line, optical forbidden lines from OI, SII, and MgI are centered at the stellar velocity and have symmetric profiles. The only way these lines could trace the photoevaporative flow is if they arise from a disk region physically distinct from that traced by the [Ne II] line, specifically from within the optically thin dust gap. However, the small (~10 km/s) FWHM of these lines suggest that most of the emitting gas traced at optical wavelengths is bound to the system rather than unbound. We discuss the implications of our results for a planet-induced versus a photoevaporation-induced gap.

We compare line emission calculated from theoretical disk models with optical to sub-millimeter wavelength observational data of the gas disk surrounding TW Hya and infer the spatial distribution of mass in the gas disk. The model disk that best matches observations has a gas mass ranging from $\sim10^{-4}-10^{-5}$\ms\ for $0.06{\rm AU} <r<3.5$AU and $\sim 0.06$\ms\ for $ 3.5 {\rm AU} <r<200$AU. We find that the inner dust hole ($r<3.5$AU) in the disk must be depleted of gas by $\sim 1-2$ orders of magnitude compared to the extrapolated surface density distribution of the outer disk. Grain growth alone is therefore not a viable explanation for the dust hole. CO vibrational emission arises within $r\sim 0.5$AU from thermal excitation of gas. [OI] 6300\AA\ and 5577\AA\ forbidden lines and OH mid-infrared emission are mainly due to prompt emission following UV photodissociation of OH and water at $r\lesssim0.1$AU and at $r\sim 4$AU. [NeII] emission is consistent with an origin in X-ray heated neutral gas at $r\lesssim 10$AU, and may not require the presence of a significant EUV ($h\nu>13.6$eV) flux from TW Hya. H$_2$ pure rotational line emission comes primarily from $r\sim 1-30$AU. [OI]63$\mu$m, HCO$^+$ and CO pure rotational lines all arise from the outer disk at $r\sim30-120$AU. We discuss planet formation and photoevaporation as causes for the decrease in surface density of gas and dust inside 4 AU. If a planet is present, our results suggest a planet mass $\sim 4-7$M$_J$ situated at $\sim 3$AU. Using our photoevaporation models and the best surface density profile match to observations, we estimate a current photoevaporative mass loss rate of $4\times10^{-9}$\ms\ yr$^{-1}$ and a remaining disk lifetime of $\sim 5$ million years.

We report sensitive Spitzer IRS spectroscopy in the 10-20 micron region of TW Hya, a nearby T Tauri star. The unusual spectral energy distribution of the source, that of a transition object, indicates that the circumstellar disk in the system has experienced significant evolution, possibly as a result of planet formation. The spectrum we measure is strikingly different from that of other classical T Tauri stars reported in the literature, displaying no strong emission features of water, acetylene, or HCN. The difference indicates that the inner planet formation region (within 5 AU) of the gaseous disk has evolved physically and/or chemically away from the classical T Tauri norm. Nevertheless, TW Hya does show a rich spectrum of emission features of atoms (HI, [NeII], and [NeIII]) and molecules (H2, OH, CO2, HCO+, and possibly CH3), some of which are also detected in classical T Tauri spectra. The properties of the neon emission are consistent with an origin for the emission in a disk irradiated by X-rays (with a possible role for additional irradiation by stellar EUV). The OH emission we detect, which also likely originates in the disk, is hot, arising from energy levels up to 23,000 K above ground, and may be produced by the UV photodissociation of water. The HI emission is surprisingly strong, with relative strengths that are consistent with case B recombination. While the absence of strong molecular emission in the 10-20 micron region may indicate that the inner region of the gaseous disk has been partly cleared by an orbiting giant planet, chemical and/or excitation effects may be responsible instead. We discuss these issues and how our results bear on our understanding of the evolutionary state of the TW Hya disk.

We present the time evolution of viscously accreting circumstellar disks as they are irradiated by ultraviolet and X-ray photons from a low-mass central star. Our model is a hybrid of a 1D time-dependent viscous disk model coupled to a 1+1D disk vertical structure model used for calculating the disk structure and photoevaporation rates. We find that disks of initial mass 0.1M_o around 1M_o stars survive for 4x10^6 years, assuming a viscosity parameter $\alpha=0.01$, a time-dependent FUV luminosity $L_{FUV}~10^{-2}-10^{-3}$ L_o and with X-ray and EUV luminosities $L_X \sim L_{EUV} ~ 10^{-3}$L_o. We find that FUV/X-ray-induced photoevaporation and viscous accretion are both important in depleting disk mass. Photoevaporation rates are most significant at ~ 1-10 AU and at >~ 30 AU. Viscosity spreads the disk which causes mass loss by accretion onto the central star and feeds mass loss by photoevaporation in the outer disk. We find that FUV photons can create gaps in the inner, planet-forming regions of the disk (~ 1-10 AU) at relatively early epochs in disk evolution while disk masses are still substantial. EUV and X-ray photons are also capable of driving gaps, but EUV can only do so at late, low accretion-rate epochs after the disk mass has already declined substantially. Disks around stars with predominantly soft X-ray fields experience enhanced photoevaporative mass loss. We follow disk evolution around stars of different masses, and find that disk survival time is relatively independent of mass for stars with M <~ 3M_o; for M >~ 3M_o the disks are short-lived(~10^5 years).

Extreme ultraviolet (EUV, 13.6 eV $< h\nu \lta 100$ eV) and X-rays in the 0.1-2 keV band can heat the surfaces of disks around young, low mass stars to thousands of degrees and ionize species with ionization potentials greater than 13.6 eV. Shocks generated by protostellar winds can also heat and ionize the same species close to the star/disk system. These processes produce diagnostic lines (e.g., [NeII] 12.8 $\mu$m and [OI] 6300 \AA) that we model as functions of key parameters such as EUV luminosity and spectral shape, X-ray luminosity and spectral shape, and wind mass loss rate and shock speed. Comparing our models with observations, we conclude that either internal shocks in the winds or X-rays incident on the disk surfaces often produce the observed [NeII] line, although there are cases where EUV may dominate. Shocks created by the oblique interaction of winds with disks are unlikely [NeII] sources because these shocks are too weak to ionize Ne. Even if [NeII] is mainly produced by X-rays or internal wind shocks, the neon observations typically place upper limits of $\lta 10^{42}$ s$^{-1}$ on the EUV photon luminosity of these young low mass stars. The observed [OI] 6300 \AA line has both a low velocity component (LVC) and a high velocity component. The latter likely arises in internal wind shocks. For the former we find that X-rays likely produce more [OI] luminosity than either the EUV layer, the transition layer between the EUV and X-ray layer, or the shear layer where the protostellar wind shocks and entrains disk material in a radial flow across the surface of the disk. Our soft X-ray models produce [OI] LVCs with luminosities up to $10^{-4}$ L$_\odot$, but may not be able to explain the most luminous LVCs.

We calculate the rate of photoevaporation of a circumstellar disk by energetic radiation (FUV, 6eV $<h\nu<$13.6eV; EUV, 13.6eV $<h\nu<$0.1keV; and Xrays, $h\nu>0.1$keV) from its central star. We focus on the effects of FUV and X-ray photons since EUV photoevaporation has been treated previously, and consider central star masses in the range $0.3-7 {\rm M}_{\odot}$. Contrary to the EUV photoevaporation scenario, which creates a gap at about $r_g\sim 7\ (M_*/1{\rm M}_{\odot})$ AU and then erodes the outer disk from inside out, we find that FUV photoevaporation predominantly removes less bound gas from the outer disk. Heating by FUV photons can cause significant erosion of the outer disk where most of the mass is typically located. X-rays indirectly increase the mass loss rates (by a factor $\sim 2$) by ionizing the gas, thereby reducing the positive charge on grains and PAHs and enhancing FUV-induced grain photoelectric heating. FUV and X-ray photons may create a gap in the disk at $\sim 10$ AU under favourable circumstances. Photoevaporation timescales for M$_* \sim 1{\rm M}_{\odot}$ stars are estimated to be $\sim 10^6$ years, after the onset of disk irradiation by FUV and X-rays. Disk lifetimes do not vary much for stellar masses in the range $0.3-3$M$_{\odot}$. More massive stars ($\gtrsim 7 {\rm M}_{\odot}$) lose their disks rapidly (in $\sim 10^5$ years) due to their high EUV and FUV fields. Disk lifetimes are shorter for shallow surface density distributions and when the dust opacity in the disk is reduced by processes such as grain growth or settling. The latter suggests that the photoevaporation process may accelerate as the dust disk evolves.

We present self-consistent models of gas in optically-thick dusty disks and calculate its thermal, density and chemical structure. The models focus on an accurate treatment of the upper layers where line emission originates, and at radii $\gtrsim 0.7$ AU. We present results of disks around $\sim 1{\rm M}_{\odot}$ stars where we have varied dust properties, X-ray luminosities and UV luminosities. We separately treat gas and dust thermal balance, and calculate line luminosities at infrared and sub-millimeter wavelengths from all transitions originating in the predominantly neutral gas that lies below the ionized surface of the disk. We find that the [ArII] 7$\mu$m, [NeII] 12.8$\mu$m, [FeI] 24$\mu$m, [SI] 25$\mu$m, [FeII] 26$\mu$m, [SiII] 35 $\mu$m, [OI] 63$\mu$m and pure rotational lines of H$_2$, H$_2$O and CO can be quite strong and are good indicators of the presence and distribution of gas in disks. We apply our models to the disk around the nearby young star, TW Hya, and find good agreement between our model calculations and observations. We also predict strong emission lines from the TW Hya disk that are likely to be detected by future facilities. A comparison of CO observations with our models suggests that the gas disk around TW Hya may be truncated to $\sim 120 $ AU, compared to its dust disk of 174 AU. We speculate that photoevaporation due to the strong stellar FUV field from TW Hya is responsible for the gas disk truncation.

We report the detection of [Ne II] emission at 12.81 micron in four out of the six optically thick dust disks observed as part of the FEPS Spitzer Legacy program. In addition, we detect a H I(7-6) emission line at 12.37 micron from the source RXJ1852.3-3700. Detections of [Ne II] lines are favored by low mid-infrared excess emission. Both stellar X-rays and extreme UV (EUV) photons can sufficiently ionize the disk surface to reproduce the observed line fluxes, suggesting that emission from Ne+ originates in the hot disk atmosphere. On the other hand, the H I(7-6) line is not associated with the gas in the disk surface and magnetospheric accretion flows can account only for at most ~30% of the observed flux. We conclude that accretion shock regions and/or the stellar corona could contribute to most of the H I(7-6)emission. Finally, we discuss the observations necessary to identify whether stellar X-rays or EUV photons are the dominant ionization mechanism for Ne atoms. Because the observed [Ne II] emission probes very small amounts of gas in the disk surface (~10^{-6} Jupiter masses) we suggest using this gas line to determine the presence or absence of gas in more evolved circumstellar disks.

We provide an overview of the Spitzer Legacy Program ``Formation and Evolution of Planetary Systems" (FEPS) which was proposed in 2000, begun in 2001, and executed aboard the Spitzer Space Telescope between 2003 and 2006. This program exploits the sensitivity of Spitzer to carry out mid-infrared spectrophotometric observations of solar-type stars. With a sample of ~ 328 stars ranging in age from ~ 3 Myr to ~ 3 Gyr, we trace the evolution of circumstellar gas and dust from primordial planet-building stages in young circumstellar disks through to older collisionally generated debris disks. When completed, our program will help define the time scales over which terrestrial and gas giant planets are built, constrain the frequency of planetesimal collisions as a function of time, and establish the diversity of mature planetary architectures. In addition to the observational program, we have coordinated a concomitant theoretical effort aimed at understanding the dynamics of circumstellar dust with and without the effects of embedded planets, dust spectral energy distributions, and atomic and molecular gas line emission. Together with the observations, these efforts will provide astronomical context for understanding whether our Solar System - and its habitable planet - is a common or a rare circumstance. Additional information about the FEPS project can be found on the team website: feps.as.arizona.edu

We have carried out a sensitive search for gas emission lines at infrared and millimeter wavelengths for a sample of 15 young sun-like stars selected from our dust disk survey with the Spitzer Space Telescope. We have used mid-infrared lines to trace the warm (300-100 K) gas in the inner disk and millimeter transitions of 12CO to probe the cold (~20 K) outer disk. We report no gas line detections from our sample. Line flux upper limits are first converted to warm and cold gas mass limits using simple approximations allowing a direct comparison with values from the literature. We also present results from more sophisticated models following Gorti and Hollenbach (2004) which confirm and extend our simple analysis. These models show that the SI line at 25.23 micron can set constraining limits on the gas surface density at the disk inner radius and traces disk regions up to a few AU. We find that none of the 15 systems have more than 0.04 MJ of gas within a few AU from the disk inner radius for disk radii from 1 AU up to ~40 AU. These gas mass upper limits even in the 8 systems younger than ~30 Myr suggest that most of the gas is dispersed early. The gas mass upper limits in the 10-40 AU region, that is mainly traced by our CO data, are <2 Mearth. If these systems are analogs of the Solar System, either they have already formed Uranus- and Neptune-like planets or they will not form them beyond 100 Myr. Finally, the gas surface density upper limits at 1 AU are smaller than 0.01% of the minimum mass solar nebula for most of the sources. If terrestrial planets form frequently and their orbits are circularized by gas, then circularization occurs early.

We report infrared spectroscopic observations of HD 105, a nearby ($\sim 40$ pc) and relatively young ($\sim 30$ Myr) G0 star with excess infrared continuum emission, which has been modeled as arising from an optically thin circumstellar dust disk with an inner hole of size $\gtrsim 13$ AU. We have used the high spectral resolution mode of the Infrared Spectrometer (IRS) on the Spitzer Space Telescope to search for gas emission lines from the disk. The observations reported here provide upper limits to the fluxes of H$_2$ S(0) 28$\mu$m, H$_2$ S(1) 17$\mu$m, H$_2$ S(2) 12 $\mu$m, [FeII] 26$\mu$m, [SiII] 35$\mu$m, and [SI] 25$\mu$m infrared emission lines. The H$_2$ line upper limits directly place constraints on the mass of warm molecular gas in the disk: $M({\rm H_2})< 4.6$, 3.8$\times 10^{-2}$, and $3.0\times 10^{-3}$ M$_J$ at $T= 50$, 100, and 200 K, respectively. We also compare the line flux upper limits to predictions from detailed thermal/chemical models of various gas distributions in the disk. These comparisons indicate that if the gas distribution has an inner hole with radius $r_{i,gas}$, the surface density at that inner radius is limited to values ranging from $\lesssim 3$ gm cm$^{-2}$ at $r_{i,gas}=0.5$ AU to 0.1 gm cm$^{-2}$ at $r_{i,gas}= 5-20$ AU. These values are considerably below the value for a minimum mass solar nebula, and suggest that less than 1 M$_J$ of gas (at any temperature) exists in the 1-40 AU planet-forming region. Therefore, it is unlikely that there is sufficient gas for gas giant planet formation to occur in HD 105 at this time.

We present 3-160 micron photometry obtained with the IRAC and MIPS instruments for the first five targets from the Spitzer Legacy Science Program "Formation and Evolution of Planetary Systems" and 4-35 micron spectro-photometry obtained with the IRS for two sources. We discuss in detail our observations of the debris disks surrounding HD 105 (G0V, 30 +- 10 Myr) and HD 150706 (G3V, ~ 700 +- 300 Myr). For HD 105, possible interpretations include large bodies clearing the dust inside of 45 AU or a reservoir of gas capable of sculpting the dust distribution. The disk surrounding HD 150706 also exhibits evidence of a large inner hole in its dust distribution. Of the four survey targets without previously detected IR excess, spanning ages 30 Myr to 3 Gyr, the new detection of excess in just one system of intermediate age suggests a variety of initial conditions or divergent evolutionary paths for debris disk systems orbiting solar-type stars.

We model gas and dust emission from regions 0.3-20 AU from a central low mass star in intermediate-aged(~10Myr) disks with optically thin dust. The models treat thermal balance and chemistry self-consistently and calculate the vertical density temperature structure of the gas in a disk. The gas and dust temperatures are calculated separately. We focus on IR gas emission lines and find that the [SI]25.23 micron line is the strongest emission line for a wide range of disk and stellar parameters, followed by emission from [SiII]34.8 um, [FeII]26 um and [OI]63um. [FeI]24um and rotational lines of OH and H2O are strong when gas masses are high (>0.1 M_J). Emission from the rotational lines of H2 is difficult to detect, unless disk gas masses are substantial (>0.1 M_J).

When stars form in small groups (N = 100 - 500 members), their circumstellar disks are exposed to little EUV radiation but a great deal of FUV radiation from massive stars in the group. This paper calculates mass loss rates for circumstellar disks exposed to external FUV radiation. Previous work treated large disks and/or intense radiation fields in which the disk radius exceeds the critical radius (supercritical disks) where the sound speed in the FUV heated layer exceeds the escape speed. This paper shows that significant mass loss still takes place for subcritical systems. Some of the gas extends beyond the disk edge (above the disk surface) to larger distances where the temperature is higher, the escape speed is lower, and an outflow develops. The evaporation rate is a sensitive function of the stellar mass and disk radius, which determine the escape speed, and the external FUV flux, which determines the temperature structure of the flow. Disks around red dwarfs are readily evaporated and shrink to disk radii of 15 AU on short time scales (10 Myr) when exposed to moderate FUV fields with $G_0$ = 3000. Although disks around solar type stars are more durable, these disks shrink to 15 AU in 10 Myr for intense FUV radiation fields with $G_0$ = 30,000; such fields exist in the central 0.7 pc of a cluster with N = 4000 stars. If our solar system formed in the presence of such strong FUV radiation fields, this mechanism could explain why Neptune and Uranus in our solar system are gas poor, whereas Jupiter and Saturn are gas rich. This mechanism for photoevaporation can also limit the production of Kuiper belt objects and can suppress giant planet formation in sufficiently large clusters, such as the Hyades, especially for disks associated with low mass stars.

We present the results of an investigation of the effects of Far Ultraviolet (FUV) radiation from hot early type OB stars on clumps in star-forming molecular clouds. Clumps in Photodissociation regions (PDRs) undergo external heating which, if rapid, creates strong photoevaporative mass flows off the clump surfaces, and drives shocks into the clumps, compressing them to high densities. The clumps lose mass on relatively short timescales. The evolution of an individual clump is found to be sensitive to its initial colunm density, the temperature of the heated surface and the ratio of the ``turn-on time'' $t_{FUV}$ of the heating flux on a clump to its initial sound crossing-time $t_{c}$. In this paper, we use spherical 1-D numerical hydrodynamic models as well as approximate analytical models to study the evolution of turbulence-generated and pressure-confined clumps in PDRs. Turbulent clumps evolve so that their column densities are equal to a critical value determined by the local FUV field, and typically have short photoevaporation timescales, $\sim 10^{4-5}$ years for a 1 M$_{\odot}$ clump in a typical star-forming region. Clumps that are confined by an interclump medium may either get completely photoevaporated, or may preserve a shielded core with a warm, dissociated, protective shell that absorbs the incident FUV flux. We compare our results with observations of some well-studied PDRs: the Orion Bar, M17SW, NGC 2023 and the Rosette Nebula. The data are consistent with both interpretations of clump origin, with a slight indication for favouring the turbulent model for clumps over pressure-confined clumps.

We will utilize the sensitivity of SIRTF through the Legacy Science Program to carry out spectrophotometric observations of solar-type stars aimed at (1) defining the timescales over which terrestrial and gas giant planets are built, from measurements diagnostic of dust/gas masses and radial distributions; and (2) establishing the diversity of planetary architectures and the frequency of planetesimal collisions as a function of time through observations of circumstellar debris disks. Together, these observations will provide an astronomical context for understanding whether our solar system - and its habitable planet - is a common or a rare circumstance. Achieving our science goals requires measuring precise spectral energy distributions for a statistically robust sample capable of revealing evolutionary trends and the diversity of system outcomes. Our targets have been selected from two carefully assembled databases of solar-like stars: (1) a sample located within 50 pc of the Sun spanning an age range from 100-3000 Myr for which a rich set of ancillary measurements (e.g. metallicity, stellar activity, kinematics) are available; and (2) a selection located between 15 and 180 pc and spanning ages from 3 to 100 Myr. For stars at these distances SIRTF is capable of detecting stellar photospheres with SNR >30 at lambda < 24 microns for our entire sample, as well as achieving SNR >5 at the photospheric limit for over 50% of our sample at lambda=70 microns. Thus we will provide a complete census of stars with excess emission down to the level produced by the dust in our present-day solar system. More information concerning our program can be found at: http://gould.as.arizona.edu/feps

The evolution of a protostellar binary system is investigated while it is embedded in its parent molecular cloud core and acted upon by gas drag due to dynamical friction. Approximate analytical results are obtained for the energy and angular momentum evolution of the orbit in the limiting cases, where the velocity is much smaller than, and much larger than the velocity dispersion of the gas. The general case is solved numerically. Dissipation causes a decay of the orbit to smaller separations and orbital eccentricity increases with time. As a binary probably spends about a few $10^6$ years embedded in the cloud, significant evolution of the orbit is expected before the dispersal of the cloud. Binary populations have been statistically generated and evolved for comparison with observations of MS and PMS binaries. Dynamical friction changes the initial frequency distribution with period of the binaries, and a maximum is obtained in the distribution as seen observationally. Decay of the orbit due to gas drag causes an evolution of long-period orbits to shorter periods and thus enables the formation of close binaries.