One possible explanation of the cavity in debris discs is the gravitational perturbation of an embedded giant planet. Planetesimals passing close to a massive body are dynamically stirred resulting in a cleared region known as the chaotic zone. Theory of overlapping mean-motion resonances predicts the width of this cavity. To test whether this cavity is identical to the chaotic zone, we investigate the formation of cavities by means of collisionless N-body simulations assuming a 1.25-10 Jupiter mass planet with eccentricities of 0-0.9. Synthetic images at millimetre wavelengths are calculated to determine the cavity properties by fitting an ellipse to 14 percent contour level. Depending on the planetary eccentricity, e_pl, the elliptic cavity wall rotates as the planet orbits with the same (e_pl<0.2) or half (e_pl>0.2) period that of the planet. The cavity centre is offset from the star along the semi-major axis of the planet with a distance of d=0.1q^-0.17e_pl^0.5 in units of cavity size towards the planet's orbital apocentre, where q is the planet-to-star mass ratio. Pericentre (apocentre) glow develops for e_pl<0.05 (e_pl>0.1), while both are present for 0.05<=e_pl<=0.1. Empirical formulae are derived for the sizes of the cavities: da_cav=2.35q^0.36 and da_cav=7.87q^0.37e_pl^0.38 for e_pl<=0.05 and e_pl>0.05, respectively. The cavity eccentricity, e_cav, equals to that of the planet only for 0.3<=e_pl<=0.6. A new method based on ALMA observations for estimating the orbital parameters and mass of the planet carving the cavity is also given.

Recent sub-millimeter observations show non-axisymmetric brightness distributions with a horseshoe-like morphology for more than a dozen transition disks. The most accepted explanation for the observed asymmetries is the accumulation of dust in large-scale vortices. Protoplanetary disks vortices can form by the excitation of Rossby-wave instability in the vicinity of a steep pressure gradient, which can develop at the edges of a giant planet carved gap or at the edges of an accretionally inactive zone. We studied the formation and evolution of vortices formed in these two distinct scenarios by means of two-dimensional locally isothermal hydrodynamic simulations. We found that the vortex formed at the edge of a planetary gap is short-lived, unless the disk is nearly inviscid. In contrast, the vortex formed at the outer edge of a dead zone is long-lived. The vortex morphology can be significantly different in the two scenarios: the vortex radial and azimuthal extensions are ~1.5 and ~3.5 times larger for the dead zone edge compared to gap models. In some particular cases, the vortex aspect ratios can be similar in the two scenarios, however, the vortex azimuthal extensions can be used to distinguish the vortex formation mechanisms. We calculate predictions for vortex observability in the sub-millimeter continuum with ALMA. We found that the azimuthal and radial extent of brightness asymmetry correlates with vortex formation process, within the limitations of alpha-viscosity prescription.

Horseshoe-shaped brightness asymmetries of several transitional discs are thought to be caused by large-scale vortices. Anticyclonic vortices are efficiently collect dust particles, therefore they can play a major role in planet formation. Former studies suggest that the disc self-gravity weakens vortices formed at the edge of the gap opened by a massive planet in discs whose masses are in the range of 0.01<=M_disc/M_*<=0.1. Here we present an investigation on the long-term evolution of the large-scale vortices formed at the viscosity transition of the discs' dead zone outer edge by means of two-dimensional hydrodynamic simulations taking disc self-gravity into account. We perform a numerical study of low mass, 0.001<=M_disc/M_*<=0.01, discs, for which cases disc self-gravity was previously neglected. The large-scale vortices are found to be stretched due to disc self-gravity even for low-mass discs with M_disc/M_*>=0.005 where initially the Toomre Q-parameter was <=50 at the vortex distance. As a result of stretching, the vortex aspect ratio increases and a weaker azimuthal density contrast develops. The strength of the vortex stretching is proportional to the disc mass. The vortex stretching can be explained by a combined action of a non-vanishing gravitational torque caused by the vortex, and the Keplerian shear of the disc. Self-gravitating vortices are subject to significantly faster decay than non-self-gravitating ones. We found that vortices developed at sharp viscosity transitions of self-gravitating discs can be described by a GNG model as long as the disc viscosity is low, i.e. alpha_dz<=10^-5.

Grid-based hydrodynamics simulations of circumstellar disks are often performed in the curvilinear coordinate system, in which the center of the computational domain coincides with the motionless star. However, the center of mass may be shifted from the star due to the presence of any non-axisymmetric mass distribution. As a result, the system exerts a gravity force on the star, causing the star to move in response, which can affect the evolution of the circumstellar disk. We aim at studying the effects of stellar motion on the evolution of protostellar and protoplanetary disks. In protostellar disks, a non-axisymmetric distribution of matter in the form of spiral arms or massive clumps can form due to gravitational instability. Protoplanetary disks can also feature non-axisymmetric structures caused by a high-mass planet or a large-scale vortex. We use 2D grid-based hydrodynamic simulations to explore the effect of stellar motion. We adopt a non-inertial polar coordinate system centered on the star, in which the stellar motion is taken into account by calculating the indirect potential caused by the non-axisymmetric disk, a high-mass planet, or a large-scale vortex. We found that the stellar motion has a moderate effect on the evolution history in protostellar disks, reducing somewhat the disk size and mass, while having a profound effect on the collapsing envelope, changing its inner shape from an initially axisymmetric to a non-axisymmetric configuration. Protoplanetary disk simulations show that the stellar motion slightly reduces the width of the gap opened by a high-mass planet, decreases the planet migration rate, and strengthens the large-scale vortices formed at the viscosity transition. We conclude that the inclusion of the indirect potential is recommended in grid-based hydrodynamics simulations of circumstellar disks which use the curvilinear coordinate system.

The Planet Formation Imager (PFI) project aims to provide a strong scientific vision for ground-based optical astronomy beyond the upcoming generation of Extremely Large Telescopes. We make the case that a breakthrough in angular resolution imaging capabilities is required in order to unravel the processes involved in planet formation. PFI will be optimised to provide a complete census of the protoplanet population at all stellocentric radii and over the age range from 0.1 to about 100 Myr. Within this age period, planetary systems undergo dramatic changes and the final architecture of planetary systems is determined. Our goal is to study the planetary birth on the natural spatial scale where the material is assembled, which is the "Hill Sphere" of the forming planet, and to characterise the protoplanetary cores by measuring their masses and physical properties. Our science working group has investigated the observational characteristics of these young protoplanets as well as the migration mechanisms that might alter the system architecture. We simulated the imprints that the planets leave in the disk and study how PFI could revolutionise areas ranging from exoplanet to extragalactic science. In this contribution we outline the key science drivers of PFI and discuss the requirements that will guide the technology choices, the site selection, and potential science/technology tradeoffs.

Motivated by recent observational and numerical studies suggesting that collapsing protostellar cores may be replenished from the local environment, we explore the evolution of protostellar cores submerged in the external counter-rotating environment. These models predict the formation of counter-rotating disks with a deep gap in the gas surface density separating the inner disk (corotating with the star) and the outer counter-rotating disk. The properties of these gaps are compared to those of planet-bearing gaps that form in disks hosting giant planets. We employ numerical hydrodynamics simulations of collapsing cores that are replenished from the local counter-rotating environment, as well as numerical hydrodynamic simulations of isolated disks hosting giant planets, to derive the properties of the gaps that form in both cases. Our numerical simulations demonstrate that counter-rotating disks can form for a wide range of mass and angular momentum available in the local environment. The gap that separates both disks has a depletion factor smaller than 1%, can be located at a distance from ten to over a hundred AU from the star, and can propagate inward with velocity ranging from 1 AU/Myr to >100 AU/Myr. Unlike our previous conclusion, the gap can therefore be a long-lived phenomenon, comparable in some cases to the lifetime of the disk itself. For a proper choice of the planetary mass, the viscous \alpha-parameter and the disk mass, the planet-bearing gaps and the gaps in counter-rotating disks may show a remarkable similarity in the gas density profile and depletion factor, which may complicate their observational differentiation.

Context: Anticyclonic vortices are considered as a favourable places for trapping dust and forming planetary embryos. On the other hand, they are massive blobs that can interact gravitationally with the planets in the disc. Aims: We aim to study how a vortex interacts gravitationally with a planet which migrates toward it or a planet which is created inside the vortex. Methods: We performed hydrodynamical simulations of a viscous locally isothermal disc using GFARGO and FARGO-ADSG. We set a stationary Gaussian pressure bump in the disc in a way that RWI is triggered. After a large vortex is established, we implanted a low mass planet in the outer disc or inside the vortex and allowed it to migrate. We also examined the effect of vortex strength on the planet migration and checked the validity of the final result in the presence of self-gravity. Results: We noticed regardless of the planet's initial position, the planet is finally locked to the vortex or its migration is stopped in a farther orbital distance in case of a stronger vortex. For the model with the weaker vortex, we studied the effect of different parameters such as background viscosity, background surface density, mass of the planet and different planet positions. In these models, while the trapping time and locking angle of the planet vary for different parameters, the main result, which is the planet-vortex locking, remains valid. We discovered that even a planet with a mass less than 5 * 10^{-7} M_{\star} comes out from the vortex and is locked to it at the same orbital distance. For a stronger vortex, both in non-self-gravitated and self-gravitating models, the planet migration is stopped far away from the radial position of the vortex. This effect can make the vortices a suitable place for continual planet formation under the condition that they save their shape during the planetary growth.

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.

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.

EX Lup is the prototype of the EXor class of eruptive young stars. These objects show optical outbursts which are thought to be related to runaway accretion onto the star. In a previous study we observed in-situ crystal formation in the disk of EX Lup during its latest outburst in 2008, making the object an ideal laboratory to investigate circumstellar crystal formation and transport. This outburst was monitored by a campaign of ground-based and Spitzer Space Telescope observations. Here we modeled the spectral energy distribution of EX Lup in the outburst from optical to millimeter wavelengths with a 2D radiative transfer code. Our results showed that the shape of the SED at optical wavelengths was more consistent with a single temperature blackbody than a temperature distribution. We also found that this single temperature component emitted 80-100 % of the total accretion luminosity. We concluded that a thermal instability, the most widely accepted model of EXor outbursts, was likely not the triggering mechanism of the 2008 outburst of EX Lup. Our mid-infrared Spitzer spectra revealed that the strength of all crystalline bands between 8 and 30 um increased right after the end of the outburst. Six months later, however, the crystallinity in the 10 um silicate feature complex decreased. Our modeling of the mid-infrared spectral evolution of EXLup showed that, although vertical mixing should be stronger during the outburst than in the quiescent phase, fast radial transport of crystals (e.g., by stellar/disk wind) was required to reproduce the observed mid-infrared spectra.

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.

KOI-13.01, a planet-sized companion in an optical double star was announced as one of the 1235 Kepler planet candidates in February 2011. The transit curves show significant distortion that was stable over the ~130 days time-span of the data. Here we investigate the phenomenon via detailed analyses of the two components of the double star and a re-reduction of the Kepler data with pixel-level photometry. Our results indicate that KOI-13 is a common proper motion binary, with two rapidly rotating components (v sin i ~ 65--70 km/s). We identify the host star of KOI-13.01 and conclude that the transit curve asymmetry is consistent with a companion orbiting a rapidly rotating, possibly elongated star on an oblique orbit. After correcting the Kepler light curve to the second light of the optical companion star, we derive a radius of 2.2 R_J for the transiter, implying an irradiated late-type dwarf, probably a hot brown dwarf rather than a planet. KOI-13 is the first example for detecting orbital obliquity for a substellar companion without measuring the Rossiter-McLaughlin effect from spectroscopy.

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.

The number of known transiting exoplanets is rapidly increasing, which has recently inspired significant interest as to whether they can host a detectable moon. Although there has been no such example where the presence of a satellite was proven, several methods have already been investigated for such a detection in the future. All these methods utilize post-processing of the measured light curves, and the presence of the moon is decided by the distribution of a timing parameter. Here we propose a method for the detection of the moon directly in the raw transit light curves. When the moon is in transit, it puts its own fingerprint on the intensity variation. In realistic cases, this distortion is too little to be detected in the individual light curves, and must be amplified. Averaging the folded light curve of several transits helps decrease the scatter, but it is not the best approach because it also reduces the signal. The relative position of the moon varies from transit to transit, the moon's wing will appear in different positions on different sides of the planet's transit. Here we show that a careful analysis of the scatter curve of the folded light curves enhances the chance of detecting the exomoons directly.

Theories of planet formation predict the birth of giant planets in the inner, dense, and gas-rich regions of the circumstellar disks around young stars. These are the regions from which strong CO emission is expected. Observations have so far been unable to confirm the presence of planets caught in formation. We have developed a novel method to detect a giant planet still embedded in a circumstellar disk by the distortions of the CO molecular line profiles emerging from the protoplanetary disk's surface. The method is based on the fact that a giant planet significantly perturbs the gas velocity flow in addition to distorting the disk surface density. We have calculated the emerging molecular line profiles by combining hydrodynamical models with semianalytic radiative transfer calculations. Our results have shown that a giant Jupiter-like planet can be detected using contemporary or future high-resolution near-IR spectrographs such as VLT/CRIRES or ELT/METIS. We have also studied the effects of binarity on disk perturbations. The most interesting results have been found for eccentric circumprimary disks in mid-separation binaries, for which the disk eccentricity - detectable from the asymmetric line profiles - arises from the gravitational effects of the companion star. Our detailed simulations shed new light on how to constrain the disk kinematical state as well as its eccentricity profile. Recent findings by independent groups have shown that core-accretion is severely affected by disk eccentricity, hence detection of an eccentric protoplanetary disk in a young binary system would further constrain planet formation theories.

In this paper we intend to provide an indirect method to detect Jovian planets by studying near infrared emission spectra originating in the protoplanetary disks around T Tauri stars. Our idea is to investigate whether a massive planet could induce any observable effect on the spectral lines emerging in the disks atmosphere. As a tracer molecule we propose CO, which is excited in the ro-vibrational fundamental band in the disk atmosphere to a distance of ~2-3 AU (depending on the stellar mass) where terrestrial planets are thought to form. The synthetic molecular spectral line profiles were calculated by an own developed semi-analytical double layer disk model. 2D gas dynamics were incorporated in the calculation of synthetic spectral lines. We demonstrate that a massive planet embedded in a protoplanetary disk strongly influences the originally circular Keplerian gas dynamics. The perturbed motion of the gas can be detected by comparing the CO line profiles in emission, which is emerge from planet-bearing to those of planet-free disk models. The planet signal has two major characteristics: a permanent line profile asymmetry, and short timescale variability correlated with the orbital phase of the giant planet. We have found that the strength of the asymmetry depends on the physical parameters of the star-planet-disk system, such as the disk inclination angle, the planetary and stellar masses, the orbital distance, and the size of the disk inner cavity. The permanent line profile asymmetry is caused by a disk in an eccentric state in the gap opened by the giant planet. However, the variable component is a consequence of the local dynamical perturbation by the orbiting giant planet. We show that a forming giant planet, still embedded in the protoplanetary disk, can be detected using contemporary or future high-resolution near-IR spectrographs like VLT/CRIRES and ELT/METIS.

We studied the statistical properties of the luminosity distance and internal extinction data of type Ia supernovae in the lists published by Tonry et al. (2003) and Barris et al. (2004). After selecting the luminosity distance in an empty Universe as a reference level we divided the sample into low $z<0.25$ and high $z \ge 0.25$ parts. We further divided these subsamples by the median of the internal extinction. Performing sign tests using the standardized residuals between the estimated logarithmic luminosity distances and those of an empty universe, on the four subsamples separately, we recognized that the residuals were distributed symmetrically in the low redshift region, independently from the internal extinction. On the contrary, the low extinction part of the data of $z \ge 0.25$ clearly showed an excess of the points with respect to an empty Universe which was not the case in the high extinction region. This diversity pointed to an interrelation between the estimated luminosity distance and internal extinction. To characterize quantitatively this interrelation we introduced a hidden variable making use of the technics of factor analysis. After subtracting that part of the residual which was explained by the hiddenmaking use of the technics of factor analysis. After subtracting that part of the residual which was explained by the hidden variable we obtained luminosity distances which were already free from interrelation with internal extinction. Fitting the corrected luminosity distances with cosmological models we concluded that the SN Ia data alone did not exclude the possibility of the $\Lambda=0$ solution.