Asymmetric fundamental band CO lines as a sign of an embedded giant planet

We investigate the formation of double-peaked asymmetric line profiles of CO in the fundamental band spectra emitted by young (1-5Myr) protoplanetary disks hosted by a 0.5-2 Solar mass star. Distortions of the line profiles can be caused by the gravitational perturbation of an embedded giant planet with q=4.7 10^-3 stellar-to-planet mass ratio. Locally isothermal, 2D hydrodynamic simulations show that the disk becomes globally eccentric inside the planetary orbit with stationary ~0.2-0.25 average eccentricity after ~2000 orbital periods. For orbital distances 1-10 AU, the disk eccentricity is peaked inside the region where the fundamental band of CO is thermal excitated. Hence, these lines become a sensitive indicators of the embedded planet via their asymmetries (both in flux and wavelength). We find that the line shape distortions (e.g. distance, central dip, asymmetry and positions of peaks) of a given transition depend on the excitation energy (i.e. on the rotational quantum number J). The magnitude of line asymmetry is increasing/decreasing with J if the planet orbits inside/outside the CO excitation zone (R_CO<=3, 5 and 7 AU for a 0.5,1 and 2 Solar mass star, respectively), thus one can constrain the orbital distance of a giant planet by determining the slope of peak asymmetry-J profile. We conclude that the presented spectroscopic phenomenon can be used to test the predictions of planet formation theories by pushing the age limits for detecting the youngest planetary systems.

Trapping of giant-planet cores - I. Vortex aided trapping at the outer dead zone edge

In this paper the migration of a 10 Earth-mass planetary core is investigated at the outer boundary of the dead zone of a protoplanetary disc by means of 2D hydrodynamic simulations done with the graphics processor unit version of the FARGO code. In the dead zone, the effective viscosity is greatly reduced due to the disc self-shielding against stellar UV radiation, X-rays from the stellar magnetosphere and interstellar cosmic rays. As a consequence, mass accumulation occurs near the outer dead zone edge, which is assumed to trap planetary cores enhancing the efficiency of the core-accretion scenario to form giant planets. Contrary to the perfect trapping of planetary cores in 1D models, our 2D numerical simulations show that the trapping effect is greatly dependent on the width of the region where viscosity reduction is taking place. Planet trapping happens exclusively if the viscosity reduction is sharp enough to allow the development of large-scale vortices due to the Rossby wave instability. The trapping is only temporarily, and its duration is inversely proportional to the width of the viscosity transition. However, if the Rossby wave instability is not excited, a ring-like axisymmetric density jump forms, which cannot trap the 10 Earth-mass planetary cores. We revealed that the stellar torque exerted on the planet plays an important role in the migration history as the barycentre of the system significantly shifts away from the star due to highly non-axisymmetric density distribution of the disc. Our results still support the idea of planet formation at density/pressure maximum, since the migration of cores is considerably slowed down enabling them further growth and runaway gas accretion in the vicinity of an overdense region.

Possible planet-forming regions on submillimetre images

Submillimetre images of transition discs are expected to reflect the distribution of the optically thin dust. Former observation of three transition discs LkHa330, SR21N, and HD1353444B at submillimetre wavelengths revealed images which cannot be modelled by a simple axisymmetric disc. We show that a large-scale anticyclonic vortex that develops where the viscosity has a large gradient (e.g., at the edge of the disc dead zone), might be accountable for these large-scale asymmetries. We modelled the long-term evolution of vortices being triggered by the Rossby wave instability. We found that a horseshoe-shaped (azimuthal wavenumber m=1) large-scale vortex forms by coalescing of smaller vortices within 5x10^4 yr, and can survive on the disc life-time (~5x10^6 yr), depending on the magnitude of global viscosity and the thickness of the viscosity gradient. The two-dimensional grid-based global disc simulations with local isothermal approximation and compressible-gas model have been done by the GPU version of hydrodynamic code FARGO (GFARGO). To calculate the dust continuum image at submillimetre wavelengths, we combined our hydrodynamical results with a 3D radiative transfer code. By the striking similarities of the calculated and observed submillimetre images, we suggest that the three transition discs can be modelled by a disc possessing a large-scale vortex formed near the disc dead zone edge. Since the larger dust grains (larger than mm in size) are collected in these vortices, the non-axisymmetric submillimetre images of the above transition discs might be interpreted as active planet and planetesimal forming regions situated far (> 50 AU) from the central stars.

Spectral signatures of disk eccentricity in young binary systems: I. Circumprimary case

Star formation occurs via fragmentation of molecular clouds, which means that the majority of stars born are a members of binaries. There is growing evidence that planets might form in circumprimary disks of medium-separation binaries. The tidal forces caused by the secondary generally act to distort the originally circular disk to an eccentric one. To infer the disk eccentricity from high-res NIR spectroscopy, we calculate the fundamental band emission lines of the CO molecule emerging from the atmosphere of the disk. We model circumprimary disk evolution under the gravitational perturbation of the orbiting secondary using a 2D grid-based hydrodynamical code, assuming alpha-type viscosity. The hydrodynamical results are combined with our spectral code based on the double-layer disk model to calculate the CO molecular line profiles. We find that the orbital velocity distribution of the gas parcels differs significantly from the circular Keplerian fashion, thus the line profiles are asymmetric in shape. The magnitude of asymmetry is insensitive to the binary mass ratio, the magnitude of viscosity, and the disk mass. In contrast, the disk eccentricity, thus the level of the line profile asymmetry, is influenced significantly by the binary eccentricity and the disk geometrical thickness. We demonstrate that the disk eccentricity profile in the planet-forming region can be determined by fitting the high-resolution CO line profile asymmetry using a simple 2D spectral model that accounts for the velocity distortions caused by the disk eccentricity. Thus, with our novel approach the disk eccentricity can be inferred with high-resolution near-IR spectroscopy prior to the era of high angular resolution optical or radio direct-imaging. By determining the disk eccentricity in medium-separation young binaries, we might be able to constrain the planet formation theories.

Fundamental Vibrational Transition of CO during the Outburst of EX Lupi in 2008

We report monitoring observations of the T Tauri star EX Lupi during its outburst in 2008 in the CO fundamental band at 4.6-5.0 um. The observations were carried out at the VLT and the Subaru Telescope at six epochs from April to August 2008, covering the plateau of the outburst and the fading phase to a quiescent state. The line flux of CO emission declines with the visual brightness of the star and the continuum flux at 5 um, but composed of two subcomponents that decay with different rates. The narrow line emission (50 km s-1 in FWHM) is near the systemic velocity of EX Lupi. These emission lines appear exclusively in v=1-0. The line widths translate to a characteristic orbiting radius of 0.4 AU. The broad line component (FWZI ~ 150 km s-1) is highly excited upto v<=6. The line flux of the component decreases faster than the narrow line emission. Simple modeling of the line profiles implies that the broad-line emitting gas is orbiting around the star at 0.04-0.4 AU. The excitation state, the decay speed of the line flux, and the line profile, indicate that the broad-line emission component is physically distinct from the narrow-line emission component, and more tightly related to the outburst event.