3D modifications to the well-studied 2D flow topology around an embedded planet have the potential to resolve long-standing problems in planet formation theory. We present a detailed analysis of the 3D isothermal flow field around a 5 Earth-mass planet on a fixed circular orbit, simulated using our multi-GPU hydrodynamics code PEnGUIn. We find that, overall, the horseshoe region has a columnar structure extending vertically much beyond the Hill sphere of the planet. This columnar structure is only broken for some of the widest horseshoe streamlines, along which high altitude fluid descends rapidly into the planet's Bondi sphere, performs one horseshoe turn, and exits the Bondi sphere radially in the midplane. A portion of this flow exits the horseshoe region altogether, which we refer to as the "transient" horseshoe flow. The flow continues as it rolls up into a pair of up-down symmetric horizontal vortex lines shed into the wake of the planet. This flow, unique to 3D, affects both planet accretion and migration. It prevents the planet from sustaining a hydrostatic atmosphere due to its intrusion into the Bondi sphere, and leads to a significant corotation torque on the planet, unanticipated by 2D analysis. In the reported simulation, starting with a $\Sigma\sim r^{-3/2}$ radial surface density profile, this torque is positive and partially cancels with the negative differential Lindblad torque, resulting in a factor of 3 slower planet migration rate. Finally, we report 3D effects can be suppressed by a sufficiently large disk viscosity, leading to results similar to 2D.

An instability can potentially operate in highly irradiated disks where the disk sharply transitions from being radially transparent to opaque (the 'transition region'). Such conditions may exist at the inner edges of transitional disks around T Tauri stars and accretion disks around AGNs. We derive the criterion for this instability, which we term the 'irradiation instability', or IRI. We also present the linear growth rate as a function of beta, the ratio between radiation force and gravity, and c_s, the sound speed of the disk, obtained using two methods: a semi-analytic analysis of the linearized equations and a numerical simulation using the GPU-accelerated hydrodynamical code PEnGUIn. In particular, we find that IRI occurs at beta~0.1 if the transition region extends as wide as ~0.05r, and at higher beta values if it is wider. This threshold value applies to c_s ranging from 3% of the Keplerian orbital speed to 5%, and becomes higher if c_s is lower. Furthermore, in the nonlinear evolution of the instability, disks with a large beta and small c_s exhibit 'clumping', extreme local surface density enhancements that can reach over ten times the initial disk surface density.

We have carried out high contrast imaging of 70 young, nearby B and A stars to search for brown dwarf and planetary companions as part of the Gemini NICI Planet-Finding Campaign. Our survey represents the largest, deepest survey for planets around high-mass stars (~1.5-2.5 M_sun) conducted to date and includes the planet hosts beta Pic and Fomalhaut. We obtained follow-up astrometry of all candidate companions within 400 AU projected separation for stars in uncrowded fields and identified new low-mass companions to HD 1160 and HIP 79797. We have found that the previously known young brown dwarf companion to HIP 79797 is itself a tight (3 AU) binary, composed of brown dwarfs with masses 58 (+21, -20) M_Jup and 55 (+20, -19) M_Jup, making this system one of the rare substellar binaries in orbit around a star. Considering the contrast limits of our NICI data and the fact that we did not detect any planets, we use high-fidelity Monte Carlo simulations to show that fewer than 20% of 2 M_sun stars can have giant planets greater than 4 M_Jup between 59 and 460 AU at 95% confidence, and fewer than 10% of these stars can have a planet more massive than 10 M_Jup between 38 and 650 AU. Overall, we find that large-separation giant planets are not common around B and A stars: fewer than 10% of B and A stars can have an analog to the HR 8799 b (7 M_Jup, 68 AU) planet at 95% confidence. We also describe a new Bayesian technique for determining the ages of field B and A stars from photometry and theoretical isochrones. Our method produces more plausible ages for high-mass stars than previous age-dating techniques, which tend to underestimate stellar ages and their uncertainties.

We present the discovery of a wide (67 AU) substellar companion to the nearby (21 pc) young solar-metallicity M1 dwarf CD-35 2722, a member of the ~100 Myr AB Doradus association. Two epochs of astrometry from the NICI Planet-Finding Campaign confirm that CD-35 2722 B is physically associated with the primary star. Near-IR spectra indicate a spectral type of L4\pm1 with a moderately low surface gravity, making it one of the coolest young companions found to date. The absorption lines and near-IR continuum shape of CD-35 2722 B agree especially well the dusty field L4.5 dwarf 2MASS J22244381-0158521, while the near-IR colors and absolute magnitudes match those of the 5 Myr old L4 planetary-mass companion, 1RXS J160929.1-210524 b. Overall, CD-35 2722 B appears to be an intermediate-age benchmark for L-dwarfs, with a less peaked H-band continuum than the youngest objects and near-IR absorption lines comparable to field objects. We fit Ames-Dusty model atmospheres to the near-IR spectra and find T=1700-1900 K and log(g) =4.5\pm0.5. The spectra also show that the radial velocities of components A and B agree to within \pm10 km/s, further confirming their physical association. Using the age and bolometric luminosity of CD-35 2722 B, we derive a mass of 31\pm8 Mjup from the Lyon/Dusty evolutionary models. Altogether, young late-M to mid-L type companions appear to be over-luminous for their near-IR spectral type compared to field objects, in contrast to the under-luminosity of young late-L and early-T dwarfs.

Our team is carrying out a multi-year observing program to directly image and characterize young extrasolar planets using the Near-Infrared Coronagraphic Imager (NICI) on the Gemini-South 8.1-meter telescope. NICI is the first instrument on a large telescope designed from the outset for high-contrast imaging, comprising a high-performance curvature adaptive optics system with a simultaneous dual-channel coronagraphic imager. Combined with state-of-the-art observing methods and data processing, NICI typically achieves ~2 magnitudes better contrast compared to previous ground-based or space-based programs, at separations inside of ~2 arcsec. In preparation for the Campaign, we carried out efforts to identify previously unrecognized young stars, to rigorously construct our observing strategy, and to optimize the combination of angular and spectral differential imaging. The Planet-Finding Campaign is in its second year, with first-epoch imaging of 174 stars already obtained out of a total sample of 300 stars. We describe the Campaign's goals, design, implementation, performance, and preliminary results. The NICI Campaign represents the largest and most sensitive imaging survey to date for massive (~1 Mjup) planets around other stars. Upon completion, the Campaign will establish the best measurements to date on the properties of young gas-giant planets at ~5-10 AU separations. Finally, Campaign discoveries will be well-suited to long-term orbital monitoring and detailed spectrophotometric followup with next-generation planet-finding instruments.

We report the discovery of a tight substellar companion to the young solar analog PZ Tel, a member of the Beta Pictoris moving group observed with high contrast adaptive optics imaging as part of the Gemini NICI Planet-Finding Campaign. The companion was detected at a projected separation of 16.4 +/- 1.0 AU (0.33 +/- 0.01") in April 2009. Second-epoch observations in May 2010 demonstrate that the companion is physically associated and shows significant orbital motion. Monte Carlo modeling constrains the orbit of PZ Tel B to eccentricities > 0.6. The near-IR colors of PZ Tel B indicate a spectral type of M7+/-2 and thus this object will be a new benchmark companion for studies of ultracool, low-gravity photospheres. Adopting an age of 12 +8 -4 Myr for the system, we estimate a mass of 36 +/- 6 Mjup based on the Lyon/DUSTY evolutionary models. PZ Tel B is one of few young substellar companions directly imaged at orbital separations similar to those of giant planets in our own solar system. Additionally, the primary star PZ Tel A shows a 70 um emission excess, evidence for a significant quantity of circumstellar dust that has not been disrupted by the orbital motion of the companion.

We present 3-D models of dust distribution around beta Pictoris that produce the best fits to the Hubble Space Telescope Advanced Camera for Surveys' (HST/ACS) images obtained by Golimowski and co-workers. We allow for the presence of either one or two separate axisymmetric dust disks. The density models are analytical, radial two-power-laws joined smoothly at a cross-over radius with density exponentially decreasing away from the mid-plane of the disks. Two-disk models match the data best, yielding a reduced chi^2 of ~1.2. Our two-disk model reproduces many of the asymmetries reported in the literature and suggests that it is the secondary (tilted) disk which is largely responsible for them. Our model suggests that the secondary disk is not constrained to the inner regions of the system (extending out to at least 250 AU) and that it has a slightly larger total area of dust than the primary, as a result of slower fall-off of density with radius and height. This surprising result raises many questions about the origin and dynamics of such a pair of disks. The disks overlap, but can coexist owing to their low optical depths and therefore long mean collision times. We find that the two disks have dust replenishment times on the order of 10^4 yr at ~100 AU, hinting at the presence of planetesimals that are responsible for the production of 2nd generation dust. A plausible conjecture, which needs to be confirmed by physical modeling of the collisional dynamics of bodies in the disks, is that the two observed disks are derived from underlying planetesimal disks; such disks would be anchored by the gravitational influence of planets located at less than 70 AU from beta Pic that are themselves in slightly inclined orbits.

We present a numerical study of rapid, so called type III migration for Jupiter-sized planets embedded in a protoplanetary disc. We limit ourselves to the case of outward migration, and study in detail its evolution and physics, concentrating on the structure of the co-rotation and circumplanetary regions, and processes for stopping migration. We also consider the dependence of the migration behaviour on several key parameters. We perform this study using global, two-dimensional hydrodynamical simulations with adaptive mesh refinement. We find that the outward directed type III migration can be started if the initial conditions support $Z > 1$, that corresponds to initial value $M_\rmn{\Delta} \ga 1.5$. Unlike the inward directed migration, in the outward migration the migration rate increases due to the growing of the volume of the co-orbital region. We find the migration to be strongly dependent on the rate of the mass accumulation in the circumplanetary disc, leading to two possible regimes of migration, fast and slow. The structure of the co-orbital region and the stopping mechanism differ between these two regimes.

We investigate the fast (type III) migration regime of high-mass protoplanets orbiting in protoplanetary disks. This type of migration is dominated by corotational torques. We study the details of flow structure in the planet's vicinity, the dependence of migration rate on the adopted disc model, and the numerical convergence of models (independence of certain numerical parameters such as gravitational softening). We use two-dimensional hydrodynamical simulations with adaptive mesh refinement,based on the FLASH code with improved time-stepping scheme. We perform global disk simulations with sufficient resolution close to the planet, which is allowed to freely move throughout the grid. We employ a new type of equation of state in which the gas temperature depends on both the distance to the star and planet, and a simplified correction for self-gravity of the circumplanetary gas. We find that the migration rate in the type III migration regime depends strongly on the gas dynamics inside the Hill sphere (Roche lobe of the planet) which, in turn, is sensitive to the aspect ratio of the circumplanetary disc. Furthermore, corrections due to the gas self-gravity are necessary to reduce numerical artifacts that act against rapid planet migration. Reliable numerical studies of Type III migration thus require consideration of both the thermal andthe self-gravity corrections, as well as a sufficient spatial resolution and the calculation of disk-planet attraction both inside and outside the Hill sphere. With this proviso, we find Type III migration to be a robust mode of migration, astrophysically promising because of a speed much faster than in the previously studied modes of migration.

We present a numerical study of rapid, so called type III migration for Jupitersized planets embedded in a protoplanetary disc. We limit ourselves to the case of inward migration, and study in detail its evolution and physics, concentrating on the structure of the corotation and circumplanetary regions, and processes for stopping migration. We also consider the dependence of the migration behaviour on several key parameters. We perform this study using the results of global, two-dimensional hydrodynamical simulations with adaptive mesh refinement. The initial conditions are chosen to satisfy the condition for rapid inward migration. We find that type III migration can be divided into two regimes, fast and slow. The structure of the coorbital region, mass accumulation rate, and migration behaviour differ between these two regimes. All our simulations show a transition from the fast to the slow regime, ending type III migration well before reaching the star. The stopping radius is found to be larger for more massive planets and less massive discs. A sharp density drop is also found to be an efficient stopping mechanism. In the fast migration limit the migration rate and induced eccentricity are lower for less massive discs, but almost do not depend on planet mass. Eccentricity is damped on the migration time scale.

We put theoretical constraints on the presence and survival of icy grains in debris discs. Particular attention is paid to UV sputtering of water ice, which has so far not been studied in detail in this context. We present a photosputtering model based on available experimental and theoretical studies. We quantitatively estimate the erosion rate of icy and ice-silicate grains, under the influence of both sublimation and photosputtering, as a function of grain size, composition and distance from the star. The effect of erosion on the grain's location is investigated through numerical simulations coupling the grain size to its dynamical evolution. Our model predicts that photodesorption efficiently destroy ice in optically thin discs, even far beyond the sublimation snow line. For the reference case of beta Pictoris, we find that only > 5mm grains can keep their icy component for the age of the system in the 50-150AU region. When taking into account the collisional reprocessing of grains, we show that the water ice survival on grains improves (grains down to ~ 20 um might be partially icy). However, estimates of the amount of gas photosputtering would produce on such a hypothetical population of big icy grains lead to values for the OI column density that strongly exceed observational constraints for beta Pic, thus ruling out the presence of a significant amount of icy grains in this system. Erosion rates and icy grains survival timescales are also given for a set of 11 other debris disc systems. We show that, with the possible exception of M stars, photosputtering cannot be neglected in calculations of icy grain lifetimes.

Vortices in protoplanetary disks can capture solid particles and form planetary cores within shorter timescales than those involved in the standard core-accretion model. We investigate vortex generation in thin unmagnetized protoplanetary disks with an embedded giant planet with planet to star mass ratio $10^{-4}$ and $10^{-3}$. Two-dimensional hydrodynamical simulations of a protoplanetary disk with a planet are performed using two different numerical methods. The results of the non-linear simulations are compared with a time-resolved modal analysis of the azimuthally averaged surface density profiles using linear perturbation theory. Finite-difference methods implemented in polar coordinates generate vortices moving along the gap created by Neptune-mass to Jupiter-mass planets. The modal analysis shows that unstable modes are generated with growth rate of order $0.3 \Omega_K$ for azimuthal numbers m=4,5,6, where $\Omega_K$ is the local Keplerian frequency. Shock-capturing Cartesian-grid codes do not generate very much vorticity around a giant planet in a standard protoplanetary disk. Modal calculations confirm that the obtained radial profiles of density are less susceptible to the growth of linear modes on timescales of several hundreds of orbital periods. Navier-Stokes viscosity of the order $\nu=10^{-5}$ (in units of $a^2 \Omega_p$) is found to have a stabilizing effect and prevents the formation of vortices. This result holds at high resolution runs and using different types of boundary conditions. Giant protoplanets of Neptune-mass to Jupiter-mass can excite the Rossby wave instability and generate vortices in thin disks. The presence of vortices in protoplanetary disks has implications for planet formation, orbital migration, and angular momentum transport in disks.

We quantitatively investigate how collisional avalanches may developin debris discs as the result of the initial break-up of a planetesimal or comet-like object, triggering a collisional chain reaction due to outward escaping small dust grains. We use a specifically developed numerical code that follows both the spatial distribution of the dust grains and the evolution of their size-frequency distribution due to collisions. We investigate how strongly avalanche propagation depends on different parameters (e.g., amount of dust released in the initial break-up, collisional properties of dust grains and their distribution in the disc). Our simulations show that avalanches evolve on timescales of ~1000 years, propagating outwards following a spiral-like pattern, and that their amplitude exponentially depends on the number density of dust grains in the system. We estimate a probability for witnessing an avalanche event as a function of disc densities, for a gas-free case around an A-type star, and find that features created by avalanche propagation can lead to observable asymmetries for dusty systems with a beta Pictoris-like dust content or higher. Characteristic observable features include: (i) a brightness asymmetry of the two sides for a disc viewed edge-on, and (ii) a one-armed open spiral or a lumpy structure in the case of face-on orientation. A possible system in which avalanche-induced structures might have been observed is the edge-on seen debris disc around HD32297, which displays a strong luminosity difference between its two sides.

We perform numerical simulations of a disc-planet system using various grid-based and smoothed particle hydrodynamics (SPH) codes. The tests are run for a simple setup where Jupiter and Neptune mass planets on a circular orbit open a gap in a protoplanetary disc during a few hundred orbital periods. We compare the surface density contours, potential vorticity and smoothed radial profiles at several times. The disc mass and gravitational torque time evolution are analyzed with high temporal resolution. There is overall consistency between the codes. The density profiles agree within about 5% for the Eulerian simulations while the SPH results predict the correct shape of the gap although have less resolution in the low density regions and weaker planetary wakes. The disc masses after 200 orbital periods agree within 10%. The spread is larger in the tidal torques acting on the planet which agree within a factor 2 at the end of the simulation. In the Neptune case the dispersion in the torques is greater than for Jupiter, possibly owing to the contribution from the not completely cleared region close to the planet.

The discovery of close orbiting extrasolar giant planets led to extensive studies of disk planet interactions and the forms of migration that can result as a means of accounting for their location. Early work established the type I and type II migration regimes for low mass embedded planets and high mass gap forming planets respectively. While providing an attractive means of accounting for close orbiting planets intially formed at several AU, inward migration times for objects in the earth mass range were found to be disturbingly short, making the survival of giant planet cores an issue. Recent progress in this area has come from the application of modern numerical techniques which make use of up to date supercomputer resources. These have enabled higher resolution studies of the regions close to the planet and the initiation of studies of planets interacting with disks undergoing MHD turbulence. This work has led to indications of how the inward migration of low to intermediate mass planets could be slowed down or reversed. In addition, the possibility of a new very fast type III migration regime, that can be directed inwards or outwards, that is relevant to partial gap forming planets in massive disks has been investigated.

We examine the effect of a rapidly migrating protoplanet on a ring of planetesimals. The eccentricities of the planetesimals are usually increased by $\Delta e \in (0.01, 0.1)$, with the exact increase being proportional to the protoplanet's mass, and inversely proportional to its migration rate. The eccentricity distribution is also substantially changed from a Rayleigh distribution. We discuss the possible implications for further planet formation, and suggest that the rapid passage of a protoplanet may not prevent the planetesimal disc from forming further planets.

We present millimeter imaging observations in the 1200 micron continuum of the disk around beta Pictoris. With the 25 arcsec beam, the beta Pic disk is unresolved perpendicularly to the disk plane (< 10 arcsec), but slightly resolved in the northeast-southwest direction (26 arcsec). Peak emission is observed at the stellar position. A secondary maximum is found 1000 AU along the disk plane in the southwest, which does not positionally coincide with a similar feature reported earlier at 850 micron. Arguments are presented which could be seen in support of the reality of these features. The observed submm/mm emission is consistent with thermal emission from dust grains, which are significantly larger than those generally found in the interstellar medium, including mm-size particles, and thus more reminiscent of the dust observed in protostellar disks. Modelling the observed scattered light in the visible and the emission in the submm/mm provides evidence for the particles dominating the scattering in the visible/NIR and those primarily responsible for the thermal emission at longer wavelengths belonging to different populations.

Since a majority of young low-mass stars are members of multiple systems, the study of their stellar and disk configurations is crucial to our understanding of both star and planet formation processes. Here we present near-infrared adaptive optics observations of the young multiple star system VW Cha. The previously known 0.7 arcsec binary is clearly resolved already in our raw J and K band images. We report the discovery of a new, faint companion to the secondary, at an apparent separation of only 0.1 arcsec or 16 AU. Our high-resolution photometric observations also make it possible to measure the J-K colors of each of the three components individually. We detect an infrared excess in the primary, consistent with theoretical models of a circumprimary disk. Analytical and numerical calculations of orbital stability show that VW Cha may be a stable triple system. Using models for the age and total mass of the secondary pair, we estimate the orbital period to be 74 years. Thus, follow-up astrometric observations might yield direct dynamical masses within a few years, and constrain evolutionary models of low-mass stars. Our results demonstrate that adaptive optics imaging in conjunction with deconvolution techniques is a powerful tool for probing close multiple systems.

We analyze the dynamics of gas-dust coupling in the presence of stellar radiation pressure in circumstellar gas disks, which are in a transitional stage between the gas-dominated, optically thick, primordial nebulae, and the dust-dominated, optically thin Vega-type disks. Dust undergo radial migration, seeking a stable equilibrium orbit in corotation with gas. The migration of dust gives rise to radial fractionation of dust and creates a variety of possible observed disk morphologies, which we compute by considering the equilibrium between the dust production and the dust-dust collisions removing particles from their equilibrium orbits. Sand-sized and larger grains are distributed throughout most of the gas disk, with concentration near the gas pressure maximum in the inner disk. Smaller grains (typically in the range of 10 to 200 micron) concentrate in a prominent ring structure in the outer region of the gas disk (presumably at radius 100 AU), where gas density is rapidly declining with radius. The width and density, as well as density contrast of the dust ring with respect to the inner dust disk depend on the distribution of gas. Our results open the prospect for deducing the distribution of gas in circumstellar disks by observing their dust. We have qualitatively compared our models with two observed transitional disks around HR 4796A and HD 141569A. Dust migration can result in observation of a ring or a bimodal radial dust distribution, possibly very similar to the ones produced by gap-opening planet(s) embedded in the disk, or shepherding it from inside or outside. We conclude that a convincing planet detection via dust imaging should include specific non-axisymmetric structure following from the dynamical simulations of perturbed disks.

We analyze the nonlinear, two-dimensional response of a gaseous, viscous protoplanetary disk to the presence of a planet of one Jupiter mass (1 M_J) and greater that orbits a 1 solar mass star by using the ZEUS hydrodynamics code with high resolution near the planet's Roche lobe. The planet is assumed to be in a circular orbit about the central star and is not allowed to migrate. A gap is formed about the orbit of the planet, but there is a nonaxisymmetric flow through the gap and onto the planet. The gap partitions the disk into an inner (outer) disk that extends inside (outside) the planet's orbit. For a 1 M_J planet and typical disk parameters, the accretion through the gap onto the planet is highly efficient. For typical disk parameters, the mass doubling time scale is less than 10^5 years, considerably shorter than the disk lifetime. Following shocks near the L1 and L2 Lagrange points, disk material enters the Roche lobe in the form of two gas streams. Shocks occur within the Roche lobe as the gas streams collide, and shocks lead to rapid inflow towards the planet within much of planet's Roche lobe. Shocks also propagate in the inner and outer disks that orbit the star. For higher mass planets (of order 6 M_J), the flow rate onto the planet is considerably reduced, which suggests an upper mass limit to planets in the range of 10 M_J. This rate reduction is related to the fact that the gap width increases relative to the Roche (Hill sphere) radius with increasing planetary mass. The flow in the gap affects planetary migration. For the 1 M_J planet case, mass can penetrate from the outer disk to the inner disk, so that the inner disk is not depleted. The results suggest that most of the mass in gas giant planets is acquired by flows through gaps.