Chemical evolution of protoplanetary disks - the effects of viscous accretion, turbulent mixing and disk winds

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 local thermodynamic equilibrium 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 that 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 sulfur-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 manner similar to that found when mixing is included.

Molecular line emission from a protoplanetary disk irradiated externally by a nearby massive star

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).

ALMA hints at the presence of two companions in the disk around HD 100546

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, for the first time, the gas and dust structure at (sub)mm wavelengths. The CO emission (at 345.795 GHz) originates from an extensive molecular disk (390 AU in radius) whereas the continuum emission is more compact (230 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 an azimuthally-symmetric disk, the continuum visibilities at long baselines (> 100 klambda) are reproduced by a compact ring with a width of 21 AU centered at 26 AU. An outer component is required to fit the short baselines: assuming a flat brightness distribution, the best-fit model is a ring with a width of 75 AU centered at 190 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 at 68 AU.

The Well-Aligned Orbit of Wasp-84b: Evidence for Disk Migration of a Hot Jupiter

We report the sky-projected orbital obliquity (spin–orbit angle) of WASP-84 b, a 0.69MJup planet in an 8.52 day orbit around a G9V/K0V star, to be λ = −0.3 ± 1.7°. We obtain a true obliquity of ψ = 17.3 ± 7.7° from a measurement of the inclination of the stellar spin axis with respect to the sky plane. Due to the young age and the weak tidal forcing of the system, we suggest that the orbit of WASP-84b is unlikely to have both realigned and circularized from the misaligned and/or eccentric orbit likely to have arisen from high-eccentricity migration. Therefore we conclude that the planet probably migrated via interaction with the protoplanetary disk. This would make it the first “hot Jupiter” (P d < 10 ) to have been shown to have migrated via this pathway. Further, we argue that the distribution of obliquities for planets orbiting cool stars (Teff < 6250 K) suggests that high-eccentricity migration is an important pathway for the formation of short-orbit, giant planets.

YSOVAR: Mid-infrared Variability Among YSOs in the Star Formation Region GGD12-15

We present an IR-monitoring survey with the Spitzer Space Telescope of the star-forming region GGD 12-15. More than 1000 objects were monitored, including about 350 objects within the central 5′, which is found to be especially dense in cluster members. The monitoring took place over 38 days and is part of the Young Stellar Object VARiability project. The region was also the subject of a contemporaneous 67 ks Chandra observation. The field includes 119 previously identified pre-main sequence star candidates. X-rays are detected from 164 objects, 90 of which are identified with cluster members. Overall, we find that about half the objects in the central 5′ are young stellar objects (YSOs) based on a combination of their spectral energy distribution, IR variability, and X-ray emission. Most of the stars with IR excess relative to a photosphere show large amplitude (>0.1 mag) mid-infrared (mid-IR) variability. There are 39 periodic sources, and all but one is found to be a cluster member. Almost half of the periodic sources do not show IR excesses. Overall, more than 85% of the Class I, flat spectrum, and Class II sources are found to vary. The amplitude of the variability is larger in more embedded YSOs. Most of the Class I/ II objects exhibit redder colors in a fainter state, which is compatible with time-variable extinction. A few become bluer when fainter, which can be explained with significant changes in the structure of the inner disk. A search for changes in the IR due to X-ray events is carried out, but the low number of flares prevented an analysis of the direct impact of X-ray flares on the IR light curves. However, we find that X-ray detected Class II sources have longer timescales for change in the MIR than a similar set of non-X-ray detected Class IIs.

First Detection of Gas-phase Methanol in a Protoplanetary Disk

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 the Atacama Large Millimeter/submillimeter Array, is spectrally resolved and detected across six velocity channels (>3σ), reaching a peak signal-to-noise of 5.5σ, with the kinematic pattern expected for TW Hya. Using an appropriate disk model, a fractional abundance of 3 x 10^-12 – 4 x 10^-11 (with respect to H₂) reproduces the stacked line profile and channel maps, with the favored 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 ≈30 au, reaching a column density ≈3–6 x 10¹² cm⁻². In the case of TW Hya, the larger (up to millimeter-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 into studying complex organic chemistry during planetary system formation.

A METHOD FOR THE STUDY OF ACCRETION DISK EMISSION IN CATACLYSMIC VARIABLES. I. THE MODEL

We have developed a spectrum synthesis method for modeling the ultraviolet (UV) emission from the accretion disk from cataclysmic variables (CVs). The disk is separated into concentric rings, with an internal structure from the Wade & Hubeny disk-atmosphere models. For each ring, a wind atmosphere is calculated in the comoving frame with a vertical velocity structure obtained from a solution of the Euler equation. Using simple assumptions, regarding rotation and the wind streamlines, these one-dimensional models are combined into a single 2.5-dimensional model for which we compute synthetic spectra. We find that the resulting line and continuum behavior as a function of the orbital inclination is consistent with the observations, and verify that the accretion rate affects the wind temperature, leading to corresponding trends in the intensity of UV lines. In general, we also find that the primary mass has a strong effect on the P Cygni absorption profiles, the synthetic emission line profiles are strongly sensitive to the wind temperature structure, and an increase in the mass-loss rate enhances the resonance line intensities. Synthetic spectra were compared with UV data for two high orbital inclination nova-like CVs-RW Tri and V347 Pup. We needed to include disk regions with arbitrary enhanced mass loss to reproduce reasonably well widths and line profiles. This fact and a lack of flux in some high ionization lines may be the signature of the presence of density-enhanced regions in the wind, or alternatively, may result from inadequacies in some of our simplifying assumptions.

Vortices in a mesoscopic superconducting disk of variable thickness

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Mean motion resonances and the stability of a circumbinary disk in a triple stellar system

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Mixing, demixing, and structure formation in a binary dipolar Bose-Einstein condensate

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

A study of low velocities regions near the Roche lobe during the gas giant planets formation

In this paper were studied regions close to the Roche lobe of a planet like Jupiter, in order to find regions with low velocities. We simulated a two dimensional and non-self-gravitating disk, where tidal and viscous torques are considered, using the hydrodynamic numerical integrator FARGO 2D. As stated earlier we are interested in find low velocities regions for in future works study the possibility of satellites formation in these regions.

Magnetic flux quantum in a Superconducting concave/convex disk

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

The star formation rate in the inner Milky Way Galaxy

The present star formation rate (SFR) in the inner Galaxy is puzzling for the chemical evolution models (CEM). No static CEM is able to reproduce the peak of the SFR in the 4 kpc ring. The main reason is probably a shortage of gas, which could be due to the dynamical effects produced by the galactic bar, not considered by these models. We developed a CEM that includes radial gas flows in order to mimic the effects of the galactic bar in the first 5 kpc of the galactic disk. In this model, the star formation (SF) is a two-step process: first, the diffuse gas forms molecular clouds. Then, stars form from cloud-cloud collisions or by the interaction between massive stars and the molecular gas. The former is called spontaneous and the latter induced SF. The mass in the different phases of each region changes by the processes associated with the stellar formation and death by: the SF due to spontaneous fragmentation of gas in the halo; formation of gas clouds in the disk from the diffuse gas; induced SF in the disk due to the interaction between massive stars and gas clouds; and finally, the restitution of the diffuse gas associated to these process of cloud and star formation. In the halo, the star formation rate for the diffuse gas follows a Schmidt law with a power n = 1.5. In the disk, the stars form in two steps: first, molecular clouds are formed from the diffuse gas also following a Schmidt law with n=1.5 and a proportionality factor. Including a specific pattern of radial gas flows, the CEM is able to reproduce with success the peak in the SFR at 4 kpc (fig. 1).

The Gemini planet-finding campaign: the frequency of giant planets aroud debris disk stars

We have completed a high-contrast direct imaging survey for giant planets around 57 debris disk stars as part of the Gemini NICI Planet-Finding Campaign. We achieved median H-band contrasts of 12.4 mag at 0.''5 and 14.1 mag at 1'' separation. Follow-up observations of the 66 candidates with projected separation <500 AU show that all of them are background objects. To establish statistical constraints on the underlying giant planet population based on our imaging data, we have developed a new Bayesian formalism that incorporates (1) non-detections, (2) single-epoch candidates, (3) astrometric and (4) photometric information, and (5) the possibility of multiple planets per star to constrain the planet population. Our formalism allows us to include in our analysis the previously known β Pictoris and the HR 8799 planets. Our results show at 95% confidence that <13% of debris disk stars have a ≥5 M Jup planet beyond 80 AU, and <21% of debris disk stars have a ≥3 M Jup planet outside of 40 AU, based on hot-start evolutionary models. We model the population of directly imaged planets as d 2 N/dMdavpropm α a β, where m is planet mass and a is orbital semi-major axis (with a maximum value of a max). We find that β < –0.8 and/or α > 1.7. Likewise, we find that β < –0.8 and/or a max < 200 AU. For the case where the planet frequency rises sharply with mass (α > 1.7), this occurs because all the planets detected to date have masses above 5 M Jup, but planets of lower mass could easily have been detected by our search. If we ignore the β Pic and HR 8799 planets (should they belong to a rare and distinct group), we find that <20% of debris disk stars have a ≥3 M Jup planet beyond 10 AU, and β < –0.8 and/or α < –1.5. Likewise, β < –0.8 and/or a max < 125 AU. Our Bayesian constraints are not strong enough to reveal any dependence of the planet frequency on stellar host mass. Studies of transition disks have suggested that about 20% of stars are undergoing planet formation; our non-detections at large separations show that planets with orbital separation >40 AU and planet masses >3 M Jup do not carve the central holes in these disks.

Molecular gas in the galaxy M83 : Its distribution, kinematics, and relation to star formation

The barred spiral galaxy M83 (NGC5236) has been observed in the 12CO J=1–0 and J=2–1 millimetre lines with the Swedish-ESO Submillimetre Telescope (SEST). The sizes of the CO maps are 100×100, and they cover the entire optical disk. The CO emission is strongly peaked toward the nucleus. The molecular spiral arms are clearly resolved and can be traced for about 360º. The total molecular gas mass is comparable to the total Hi mass, but H2 dominates in the optical disk. Iso-velocity maps show the signature of an inclined, rotating disk, but also the effects of streaming motions along the spiral arms. The dynamical mass is determined and compared to the gas mass. The pattern speed is determined from the residual velocity pattern, and the locations of various resonances are discussed. The molecular gas velocity dispersion is determined, and a trend of decreasing dispersion with increasing galactocentric radius is found. A total gas (H2+Hi+He) mass surface density map is presented, and compared to the critical density for star formation of an isothermal gaseous disk. The star formation rate (SFR) in the disk is estimated using data from various star formation tracers. The different SFR estimates agree well when corrections for extinctions, based on the total gas mass map, are made. The radial SFR distribution shows features that can be associated with kinematic resonances. We also find an increased star formation efficiency in the spiral arms. Different Schmidt laws are fitted to the data. The star formation properties of the nuclear region, based on high angular resolution HST data, are also discussed.