A large fraction of stars is found to be part of binary or higher-order multiple systems. The ubiquity of planets found around single stars raises the question if and how planets in binary systems may form. Protoplanetary disks are the birthplaces of planets, and their characterization is crucial in order to understand the planet formation process. Our aim is to characterize the morphology of the GG Tau A disk, one of the largest and most massive circumbinary disks, and trace evidence for binary-disk interactions. We obtained observations in polarized scattered light of GG Tau A using the SPHERE/IRDIS instrument in the H-band filter. We analyze the observed disk morphology and substructures. We run 2D hydrodynamical models simulating the evolution of the circumbinary ring over the lifetime of the disk. The disk, as well as the cavity and the inner region are highly structured with several shadowed regions, spiral structures, and streamer-like filaments, some of them detected for the first time. The streamer-like filaments appear to connect the outer ring with the northern arc. Their azimuthal spacing suggests that they may be generated by periodic perturbations by the binary, tearing off material from the inner edge of the outer disk once during each orbit. By comparing observations to hydrodynamical simulations we find that the main features, in particular the gap size, as well as the spiral and streamer filaments, can be qualitatively explained by the gravitational interactions of a binary with semi-major axis of $\sim$35 au on an orbit coplanar with the circumbinary ring.

We calculated hydrogen recombination line luminosities (H-$\alpha$, Paschen-$\beta$ and Brackett-$\gamma$) from three dimensional thermo-hydrodynamical simulations of forming planets from 1 to 10 Jupiter-masses. We explored various opacities to estimate the line emissions with extinction in each cases. When realistic opacities are considered, only lines from planets $\ge$10 Jupiter-mass can be detected with current instrumentation, highlighting that from most planets one cannot expect detectable emission. This might explain the very low detection rate of H-$\alpha$ from forming planets from observations. While the line emission comes from both the forming planet and its circumplanetary disk, we found that only the disk component could be detected due to extinction. We examined the line variability as well, and found that it is higher for higher mass planets. Furthermore, we determine for the first time, the parametric relationship between the mass of the planet and the luminosity of the hydrogen recombination lines, as well as the equation between the accretion luminosity and hydrogen recombination line luminosities. These relationships for planets appear to be steeper than the analogous relation for young stars, suggesting that the accretion process around planets proceeds differently than around stars.

We investigated the formation and evolution of satellite systems in a cold, extended circumplanetary disc around a 10 $M_{\rm{Jupiter}}$ gas giant which was generated by gravitational instability. We use a population synthesis approach, in which we seed the disc with satellite embryos and let them migrate, accrete mass, collide, with subsequent generations being created until the disc is dissipated. In each run we choose randomly the dust-to-gas ratio, dispersion- and refilling time-scales, the number of embryos and their starting locations in a plausible range. The disc structure is the result of a 3D global disc SPH simulation \citep{Szulagyi16b}. We also investigate the effect of the planet's semi-major axis on the resulting satellite systems, taking 50 AU as the nominal case. In the nominal case we find that most satellites are close in mass to the Galilean ones, with a maximum at around 3 $M_{\rm{Earth}}$ and form typically on time-scales comparable to the dispersion time-scale. We also find that about 10 $M_{\rm{Earth}}$ worth of satellites migrate into the planet, polluting it with metals. The influence of a different planet semi-major axis is mainly felt through the change in disc size. We find that for the discs closer to the star, the satellites are lighter, formation time-scales are longer, and more moons are lost into the planet, as a smaller disc means easier migration for the massive satellites. The probability of detecting satellites like the ones in this work is very low ($\leq$ 3%) even with upcoming powerful telescopes like E-ELT.

Chondrules are crystallised droplets of silicates formed by rapid heating to high temperatures (> 1800 K) of solid precursors followed by hours or days of cooling. Dating of chondrules is consistent with the formation timescale of Jupiter in the core-accretion model (1-4 Myrs). Here we investigate if the shocks generated by a massive planet could generate flash heating episodes necessary to form chondrules using high resolution 2D simulations with the multi-fluid code ROSSBI. We use different radiative cooling prescriptions, i.e. different cooling rates and models, and vary planet mass, orbit and disk models. For long disk cooling rates (> 1000 orbits) and a massive protoplanet (> 0.75 jupiter masses), we obtain hot enough temperatures for chondrule formation, while using more realistic thermodynamics is not successful in the Minimum Mass Solar Nebula (MMSN) model. However, sufficient flash heating is achieved by a Jupiter mass planet in a 5 times more massive disk, which is a conceivable scenario for the young solar nebula and exoplanetary systems. Interestingly, a gap-forming massive planet triggers vortices, which can trap dust, i.e. chondrule precursors, and generates a high pressure environment that is consistent with cosmochemical evidence from chondrules. A massive gas giant can thus, in principle, both stimulate the concentration of chondrule precursors in the vicinity of the shocking regions, and flash-heat them via the shocks.

Debris disks are the natural by-products of the planet formation process. Scattered or polarized light observations are mostly sensitive to small dust grains that are released from the grinding down of bigger planetesimals. High angular resolution observations at optical wavelengths can provide key constraints on the radial and azimuthal distribution of the small dust grains. These constraints can help us better understand where most of the dust grains are released upon collisions. We present SPHERE/ZIMPOL observations of the debris disk around HR 4796 A, and model the radial profiles along several azimuthal angles of the disk with a code that accounts for the effect of stellar radiation pressure. This enables us to derive an appropriate description for the radial and azimuthal distribution of the small dust grains. Even though we only model the radial profiles along (or close to) the semi-major axis of the disk, our best-fit model is not only in good agreement with our observations but also with previously published datasets (from near-IR to sub-mm wavelengths). We find that the reference radius is located at $76.4\pm0.4$ au, and the disk has an eccentricity of $0.076_{-0.010}^{+0.016}$, with the pericenter located on the front side of the disk (north of the star). We find that small dust grains must be preferentially released near the pericenter to explain the observed brightness asymmetry. Even though parent bodies spend more time near the apocenter, the brightness asymmetry implies that collisions happen more frequently near the pericenter of the disk. Our model can successfully reproduce the shape of the outer edge of the disk, without having to invoke an outer planet shepherding the debris disk. With a simple treatment of the effect of the radiation pressure, we conclude that the parent planetesimals are located in a narrow ring of about $3.6$ au in width.

There are growing amount of very high-resolution polarized scattered light images of circumstellar disks. Naturally, the question arises whether the circumplanetary disk forming around nascent planets can be detected with the same technique. Here we created scattered light mock observations at 1.2 and 1.6 microns for instruments like SPHERE and GPI, for various planetary masses and disk inclinations. We found that the detection of a circumplanetary disk is significantly favored if the planet is massive ($\geq 5 \rm{M_{Jup}}$) and the system is nearly face-on ($\leq 30^\circ$). Its detection is hindered by the neighboring circumstellar disk that also provides a strong polarized flux. However, the comparison between the $PI$ and the $Q_\phi$ maps, as well as the contrasts between the J and H bands are viable tools to pinpoint the presence of the circumplanetary disk within the circumstellar disk, as the two disks are behaving differently on those images.

Detection of forming planets means detection of the circumplanetary disk (CPD) in reality, since the planet is still surrounded by a disk at this evolutionary stage. Yet, no comprehensive CPD modeling was done in near-infrared wavelengths, where high contrast imaging is a powerful tool to detect these objects. We combined 3D radiative hydrodynamic simulations of various embedded planets with RADMC-3D radiative transfer post-processing that includes scattering of photons on dust particles. We made synthetic images for VLT NaCo/ERIS in the Ks, L', M' bands as well as examined the spectral energy distributions (SEDs) of disks between 1 $\mu m$ and 10 cm. We found that the observed magnitudes from the planet's vicinity will mostly depend on the CPD parameters, not on the planet's. The CPD is 20-100x brighter than the embedded planet in near-IR. We also show how the CPD parameters, e.g. the dust-to-gas ratio will affect the resulting CPD magnitudes. According to the SEDs, the best contrast ratio between the CPD and circumstellar disks is in sub-mm/radio wavelengths and between 8-33 microns in case if the planet opened a resolvable, deep gap ($\ge 5 \rm{M_{Jup}}$), while the contrast is particularly poor in the near-IR. Hence, to detect the forming planet and its CPD, the best chance today is targeting the sub-mm/radio wavelengths and the 10-micron silicate feature vicinity. In order to estimate the forming planet's mass from the observed brightness, it is necessary to run system specific disk modeling.

Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been allowing for increasingly stringent tests of the theoretical models. Modeling efforts are crucial to support the interpretation of direct imaging analyses, not just for potential detections but also to put meaningful upper limits on mass accretion rates and other physical quantities in current and future large-scale surveys. This white paper addresses the questions of what efforts on the computational side are required in the next decade to advance our theoretical understanding, explain the observational data, and guide new observations. We identified the nature of accretion, ab initio planet formation, early evolution, and circumplanetary disks as major fields of interest in computational planet formation. We recommend that modelers relax the approximations of alpha-viscosity and isothermal equations of state, on the grounds that these models use flawed assumptions, even if they give good visual qualitative agreement with observations. We similarly recommend that population synthesis move away from 1D hydrodynamics. The computational resources to reach these goals should be developed during the next decade, through improvements in algorithms and the hardware for hybrid CPU/GPU clusters. Coupled with high angular resolution and great line sensitivity in ground based interferometers, ELTs and JWST, these advances in computational efforts should allow for large strides in the field in the next decade.

Aims: We want to detect and quantify observables related to accretion processes occurring locally in circumstellar disks, which could be attributed to young forming planets. We focus on objects known to host protoplanet candidates and/or disk structures thought to be the result of interactions with planets. Methods: We analyzed observations of 6 young stars (age $3.5-10$ Myr) and their surrounding environments with the SPHERE/ZIMPOL instrument on the VLT in the H$\alpha$ filter (656 nm) and a nearby continuum filter (644.9 nm). Results: We re-detect the known accreting M-star companion HD142527 B with the highest published signal to noise to date in both H$\alpha$ and the continuum. We derive new astrometry ($r = 62.8^{+2.1}_{-2.7}$ mas and $\text{PA} = (98.7\,\pm1.8)^\circ$) and photometry ($\Delta$N_Ha=$6.3^{+0.2}_{-0.3}$ mag, $\Delta$B_Ha=$6.7\pm0.2$ mag and $\Delta$Cnt_Ha=$7.3^{+0.3}_{-0.2}$ mag) for the companion in agreement with previous studies, and estimate its mass accretion rate ($\dot{M}\approx1-2\,\times10^{-10}\,M_\odot\text{ yr}^{-1}$). A faint point-like source around HD135344 B (SAO206462) is also investigated, but a second deeper observation is required to reveal its nature. No other companions are detected. In the framework of our assumptions we estimate detection limits at the locations of companion candidates around HD100546, HD169142 and MWC758 and calculate that processes involving H$\alpha$ fluxes larger than $\sim8\times10^{-14}-10^{-15}\,\text{erg/s/cm}^2$ ($\dot{M}>10^{-10}-10^{-12}\,M_\odot\text{ yr}^{-1}$) can be excluded. Furthermore, flux upper limits of $\sim10^{-14}-10^{-15}\,\text{erg/s/cm}^2$ ($\dot{M}<10^{-11}-10^{-12}\,M_\odot \text{ yr}^{-1}$) are estimated within the gaps identified in the disks surrounding HD135344B and TW Hya.

We present long baseline Atacama Large Millimeter/submillimeter Array (ALMA) observations of the 870$\,\mu$m dust continuum emission and CO (3-2) from the protoplanetary disk around the Herbig Ae/Be star HD 100546, which is one of the few systems claimed to have two young embedded planets. These observations achieve a resolution of 4 au (3.8 mas), an rms noise of 66$\mu$Jy/beam, and reveal an asymmetric ring between $\sim$20-40 au with largely optically thin dust continuum emission. This ring is well fit by two concentric and overlapping Gaussian rings of different widths and a Vortex. In addition, an unresolved component is detected at a position consistent with the central star, which may trace the central inner disk ($<$2au in radius). We report a lack of compact continuum emission at the positions of both claimed protoplanets. We use this result to constrain the circumplanetary disk (CPD) mass and size of 1.44M$_{\rm Earth}$ and 0.44au in the optically thin and thick regime, respectively, for the case of the previously directly imaged protoplanet candidate at $\sim$55 au (HD100546 b). We compare these empirical CPD constraints to previous numerical simulations. This suggests that HD100546 b is inconsistent with several planet accretion models, while gas-starved models are also still compatible. We estimate the planetary mass as 1.65 M$_J$ by using the relation between planet, circumstellar, and circumplanetary masses derived from numerical simulations. Finally, the CO integrated intensity map shows a possible spiral arm feature that would match the spiral features identified in Near-Infrared scattered light polarized emission, which suggests a real spiral feature in the disk surface that needs to be confirmed with further observations.

Satellites of giant planets thought to form in gaseous circumplanetary disks (CPDs) during the late planet-formation phase, but it was unknown so far whether smaller mass planets, such as the ice giants could form such disks, thus moons there. We combined radiative hydrodynamical simulations with satellite population synthesis to investigate the question in the case of Uranus and Neptune. For both ice giants we found that a gaseous CPD is created at the end of their formation. The population synthesis confirmed that Uranian-like, icy, prograde satellite-system could form in these CPDs within a couple of $10^5$ years. This means that Neptune could have a Uranian-like moon-system originally that was wiped away by the capture of Triton. Furthermore, the current moons of Uranus can be reproduced by our model without the need for planet-planet impact to create a debris disk for the moons to grow. These results highlight that even ice giants -- that among the most common mass-category of exoplanets -- can also form satellites, opening a way to a potentially much larger population of exomoons than previously thought.

It is believed that satellites of giant planets form in circumplanetary disks. Many of the previous contributions assumed that their formation process proceeds similarly to rocky planet formation, via accretion of the satellite seeds, called satellitesimals. However, the satellitesimal formation itself poses a nontrivial problem as the dust evolution in the circumplanetary disk is heavily impacted by fast radial drift and thus dust growth to satellitesimals is hindered. To address this problem, we connected state-of-the-art hydrodynamical simulations of a circumplanetary disk around a Jupiter-mass planet with dust growth and drift model in a post-processing step. We found that there is an efficient pathway to satellitesimal formation if there is a dust trap forming within the disk. Thanks to the natural existence of an outward gas flow region in the hydrodynamical simulation, a significant dust trap arises at the radial distance of 85~R$_{\rm J}$ from the planet, where the dust-to-gas ratio becomes high enough to trigger streaming instability. The streaming instability leads to the efficient formation of the satellite seeds. Because of the constant infall of material from the circumstellar disk and the very short timescale of dust evolution, the circumplanetary disk acts as a satellitesimal factory, constantly processing the infalling dust to pebbles that gather in the dust trap and undergo the streaming instability.

The satellites of Jupiter are thought to form in a circumplanetary disc. Here we address their formation and orbital evolution with a population synthesis approach, by varying the dust-to-gas ratio, the disc dispersal timescale and the dust refilling timescale. The circumplanetary disc initial conditions (density and temperature) are directly drawn from the results of 3D radiative hydrodynamical simulations. The disc evolution is taken into account within the population synthesis. The satellitesimals were assumed to grow via streaming instability. We find that the moons form fast, often within $10^4$ years, due to the short orbital timescales in the circumplanetary disc. They form in sequence, and many are lost into the planet due to fast type I migration, polluting Jupiter's envelope with typically 15 Earth-masses of metals. The last generation of moons can form very late in the evolution of the giant planet, when the disc has already lost more than the 99% of its mass. The late circumplanetary disc is cold enough to sustain water ice, hence not surprisingly the 85% of the moon population has icy composition. The distribution of the satellite-masses is peaking slightly above Galilean masses, up until a few Earth-masses, in a regime which is observable with the current instrumentation around Jupiter-analog exoplanets orbiting sufficiently close to their host stars. We also find that systems with Galilean-like masses occur in 20% of the cases and they are more likely when discs have long dispersion timescales and high dust-to-gas ratios.

Debris disks are usually detected through the infrared excess over the photospheric level of their host star. The most favorable stars for disk detection are those with spectral types between A and K, while the statistics for debris disks detected around low-mass M-type stars is very low, either because they are rare or because they are more difficult to detect. Terrestrial planets, on the other hand, may be common around M-type stars. Here, we report on the discovery of an extended (likely) debris disk around the M-dwarf GSC 07396-00759. The star is a wide companion of the close accreting binary V4046 Sgr. The system probably is a member of the $\beta$ Pictoris Moving Group. We resolve the disk in scattered light, exploiting high-contrast, high-resolution imagery with the two near-infrared subsystems of the VLT/SPHERE instrument, operating in the YJ bands and the H2H3 doublet. The disk is clearly detected up to 1.5" ($\sim110$ au) from the star and appears as a ring, with an inclination $i\sim83$ degree, and a peak density position at $\sim 70$ au. The spatial extension of the disk suggests that the dust dynamics is affected by a strong stellar wind, showing similarities with the AU Mic system that has also been resolved with SPHERE. The images show faint asymmetric structures at the widest separation in the northwest side. We also set an upper limit for the presence of giant planets to $2 M_J$. Finally, we note that the 2 resolved disks around M-type stars of 30 such stars observed with SPHERE are viewed close to edge-on, suggesting that a significant population of debris disks around M dwarfs could remain undetected because of an unfavorable orientation.

The SHINE program is a large high-contrast near-infrared survey of 600 young, nearby stars. It is aimed at searching for and characterizing new planetary systems using VLT/SPHERE's unprecedented high-contrast and high-angular resolution imaging capabilities. It also intends at placing statistical constraints on the occurrence and orbital properties of the giant planet population at large orbits as a function of the stellar host mass and age to test planet formation theories. We use the IRDIS dual-band imager and the IFS integral field spectrograph of SPHERE to acquire high-constrast coronagraphic differential near-infrared images and spectra of the young A2 star HIP65426. It is a member of the ~17 Myr old Lower Centaurus-Crux association. At a separation of 830 mas (92 au projected) from the star, we detect a faint red companion. Multi-epoch observations confirm that it shares common proper motion with HIP65426. Spectro-photometric measurements extracted with IFS and IRDIS between 0.95 and 2.2um indicate a warm, dusty atmosphere characteristic of young low surface-gravity L5-L7 dwarfs. Hot-start evolutionary models predict a luminosity consistent with a 6-12 MJup, Teff=1300-1600 K and R=1.5 RJup giant planet. Finally, the comparison with Exo-REM and PHOENIX BT-Settl synthetic atmosphere models gives consistent effective temperatures but with slightly higher surface gravity solutions of log(g)=4.0-5.0 with smaller radii (1.0-1.3 RJup). Given its physical and spectral properties, HIP65426b occupies a rather unique placement in terms of age, mass and spectral-type among the currently known imaged planets. It represents a particularly interesting case to study the presence of clouds as a function of particle size, composition, and location in the atmosphere, to search for signatures of non-equilibrium chemistry, and finally to test the theory of planet formation and evolution.

HD$\,$169142 is an excellent target to investigate signs of planet-disk interaction due to the previous evidence of gap structures. We performed J-band (~1.2{\mu}m) polarized intensity imaging of HD169142 with VLT/SPHERE. We observe polarized scattered light down to 0.16" (~19 au) and find an inner gap with a significantly reduced scattered light flux. We confirm the previously detected double ring structure peaking at 0.18" (~21 au) and 0.56" (~66 au), and marginally detect a faint third gap at 0.70"-0.73" (~82-85 au). We explore dust evolution models in a disk perturbed by two giant planets, as well as models with a parameterized dust size distribution. The dust evolution model is able to reproduce the ring locations and gap widths in polarized intensity, but fails to reproduce their depths. It, however, gives a good match with the ALMA dust continuum image at 1.3 mm. Models with a parameterized dust size distribution better reproduce the gap depth in scattered light, suggesting that dust filtration at the outer edges of the gaps is less effective. The pile-up of millimeter grains in a dust trap and the continuous distribution of small grains throughout the gap likely require a more efficient dust fragmentation and dust diffusion in the dust trap. Alternatively, turbulence or charging effects might lead to a reservoir of small grains at the surface layer that is not affected by the dust growth and fragmentation cycle dominating the dense disk midplane. The exploration of models shows that extracting planet properties such as mass from observed gap profiles is highly degenerate.

We present mock observations of forming planets with ALMA. The possible detections of circumplanetary disks (CPDs) were investigated around planets of Saturn, 1, 3, 5, and 10 Jupiter-masses that are placed at 5.2 AU from their star. The radiative, three dimensional hydrodynamic simulations were then post-processed with RADMC3D and the ALMA Observation Simulator. We found that even though the CPDs are too small to be resolved, they are hot due to the accreting planet in the optically thick limit, therefore the best chance to detect them with continuum observations in this case is at the shortest ALMA wavelengths, such as Band 9 (440 microns). Similar fluxes were found in the case of Saturn and Jupiter-mass planets, as for the 10 $\mathrm{M_{Jup}}$ gas-giant, due to temperature weighted optical depth effects: when no deep gap is carved, the planet region is blanketed by the optically thick circumstellar disk leading to a less efficient cooling there. A test was made for a 52 AU orbital separation, showed that optically thin CPDs are also detectable in band 7 but they need longer integration times ($>$5hrs). Comparing the gap profiles of the same simulation at various ALMA bands and the hydro simulation confirmed that they change significantly, first because the gap is wider at longer wavelengths due to decreasing optical depth; second, the beam convolution makes the gap shallower and at least 25% narrower. Therefore, caution has to be made when estimating planet masses based on ALMA continuum observations of gaps.

We have observed the protoplanetary disk of the well-known young Herbig star HD 142527 using ZIMPOL Polarimetric Differential Imaging with the VBB (Very Broad Band, ~600-900nm) filter. We obtained two datasets in May 2015 and March 2016. Our data allow us to explore dust scattering around the star down to a radius of ~0.025" (~4au). The well-known outer disk is clearly detected, at higher resolution than before, and shows previously unknown sub-structures, including spirals going inwards into the cavity. Close to the star, dust scattering is detected at high signal-to-noise ratio, but it is unclear whether the signal represents the inner disk, which has been linked to the two prominent local minima in the scattering of the outer disk, interpreted as shadows. An interpretation of an inclined inner disk combined with a dust halo is compatible with both our and previous observations, but other arrangements of the dust cannot be ruled out. Dust scattering is also present within the large gap between ~30 and ~140au. The comparison of the two datasets suggests rapid evolution of the inner regions of the disk, potentially driven by the interaction with the close-in M-dwarf companion, around which no polarimetric signal is detected.

Circumplanetary disks regulate the late accretion to the giant planet and serve as the birthplace for satellites. Understanding their characteristics via simulations also helps to prepare for their observations. Here we study disks around 1, 3, 5, 10 $\mathrm{M_{Jup}}$ planets with three dimensional, global radiative hydrodynamic simulations with sub-planet peak resolution, and various planetary temperatures. We found that as the 1 $\mathrm{M_{Jup}}$ planet radiates away its formation heat, the circumplanetary envelope transitions to a disk between $T_p = 6000$ K and 4000 K. In the case of 3-10 $\mathrm{M_{Jup}}$ planets a disk always forms. The temperature profile of the circumplanetary disks is very steep, the inner 1/6th is over the silicate condensation temperature and the entire disk is above water freezing point, making satellite formation impossible in this early stage ($<$1 Myr). Satellites might form much later and first in the outer parts of the disk migrating inwards later on. Our disk masses are $1, 7, 20, 40 \times 10^{-3}\mathrm{M_{Jup}}$ for the 1, 3, 5, 10 $\mathrm{M_{Jup}}$ gas giants respectively, and we provide an empirical formula to estimate the subdisk masses based on the planet- and circumstellar disk mass. Our finding is that the cooler the planet, the lower the temperature of the subdisk, the higher the vertical influx velocities, and the planetary gap is both deeper and wider. We also show that the gaps in 2D and 3D are different. The subdisk eccentricity increases with planetary mass and violently interacts with the circumstellar disk, making satellite-formation less likely, if $\mathrm{M_p} \gtrsim 5 \mathrm{M_{Jup}}$.

The transition disk around the T Tauri star T Cha possesses a large gap, making it a prime target for high-resolution imaging in the context of planet formation. We aim to find signs of disk evolutionary processes by studying the disk geometry and the dust grain properties at its surface, and to search for companion candidates. We analyze a set of VLT/SPHERE data at near-infrared and optical wavelengths. We performed polarimetric imaging of T Cha with IRDIS (1.6 $\mu$m) and ZIMPOL (0.5-0.9 $\mu$m), and obtained intensity images from IRDIS dual-band imaging with simultaneous spectro-imaging with IFS (0.9-1.3 $\mu$m). The disk around T Cha is detected in all observing modes and its outer disk is resolved in scattered light with unprecedented angular resolution and signal-to-noise. The images reveal a highly inclined disk with a noticeable east-west brightness asymmetry. The significant amount of non-azimuthal polarization signal in the $U_{\phi}$ images, with a $U_{\phi}$/$Q_{\phi}$ peak-to-peak value of 14%, is in accordance with theoretical studies on multiple scattering. Our optimal axisymmetric radiative transfer model considers two coplanar inner and outer disks, separated by a gap of 0.28" (~30au) in size. We derive a disk inclination of ~69 deg and PA of ~114 deg. In order to self-consistently reproduce the intensity and polarimetric images, the dust grains, responsible for the scattered light, need to be dominated by sizes of around ten microns. A point source is detected at an angular distance of 3.5" from the central star. It is, however, found not to be co-moving. We confirm that the dominant source of emission is forward scattered light from the near edge of the outer disk. Our point source analysis rules out the presence of a companion heavier than ~8.5 $M_{\mathrm{jup}}$ between 0.1" and 0.3". The detection limit decreases to ~2 $M_{\mathrm{jup}}$ for 0.3" to 4.0".

ASTEP (Antarctica Search for Transiting ExoPlanets) is a pilot project that aims at searching and characterizing transiting exoplanets from Dome C in Antarctica and to qualify this site for photometry in the visible. Two instruments were installed at Dome C and ran for six winters in total. The analysis of the collected data is nearly complete. We present the operation of the instruments, and the technical challenges, limitations, and possible solutions in light of the data quality. The instruments performed continuous observations during the winters. Human interventions are required mainly for regular inspection and ice dust removal. A defrosting system is efficient at preventing and removing ice on the mirrors. The PSF FWHM is 4.5 arcsec on average which is 2.5 times larger than the specification, and is highly variable; the causes are the poor ground-level seeing, the turbulent plumes generated by the heating system, and to a lower extent the imperfect optical alignment and focusing, and some astigmatism. We propose solutions for each of these aspects that would largely increase the PSF stability. The astrometric and guiding precisions are satisfactory and would deserve only minor improvements. Major issues are encountered with the camera shutter which did not close properly after two winters; we minimized this issue by heating the shutter and by developing specific image calibration algorithms. Finally, we summarize the site testing and science results obtained with ASTEP. Overall, the ASTEP experiment will serve as a basis to design and operate future optical and near-infrared telescopes in Antarctica.

The luminosity of young giant planets can inform about their formation and accretion history. The directly imaged planets detected so far are consistent with the "hot-start" scenario of high entropy and luminosity. If nebular gas passes through a shock front before being accreted into a protoplanet, the entropy can be substantially altered. To investigate this, we present high resolution, 3D radiative hydrodynamic simulations of accreting giant planets. The accreted gas is found to fall with supersonic speed in the gap from the circumstellar disk's upper layers onto the surface of the circumplanetary disk and polar region of the protoplanet. There it shocks, creating an extended hot supercritical shock surface. This shock front is optically thick, therefore, it can conceal the planet's intrinsic luminosity beneath. The gas in the vertical influx has high entropy which when passing through the shock front decreases significantly while the gas becomes part of the disk and protoplanet. This shows that circumplanetary disks play a key role in regulating a planet's thermodynamic state. Our simulations furthermore indicate that around the shock surface extended regions of atomic - sometimes ionized - hydrogen develop. Therefore circumplanetary disk shock surfaces could influence significantly the observational appearance of forming gas-giants.

Circumplanetary disks can be found around forming giant planets, regardless of whether core accretion or gravitational instability built the planet. We carried out state-of-the-art hydrodynamical simulations of the circumplanetary disks for both formation scenarios, using as similar initial conditions as possible to unveil possible intrinsic differences in the circumplanetary disk mass and temperature between the two formation mechanisms. We found that the circumplanetary disks mass linearly scales with the circumstellar disk mass. Therefore, in an equally massive protoplanetary disk, the circumplanetary disks formed in the disk instability model can be only a factor of eight more massive than their core-accretion counterparts. On the other hand, the bulk circumplanetary disk temperature differs by more than an order of magnitude between the two cases. The subdisks around planets formed by gravitational instability have a characteristic temperature below 100 K, while the core accretion circumplanetary disks are hot, with temperatures even greater than 1000 K when embedded in massive, optically thick protoplanetary disks. We explain how this difference can be understood as the natural result of the different formation mechanisms. We argue that the different temperatures should persist up to the point when a full-fledged gas giant forms via disk instability, hence our result provides a convenient criteria for observations to distinguish between the two main formation scenarios by measuring the bulk temperature in the planet vicinity.

We present three-dimensional simulations with nested meshes of the dynamics of the gas around a Jupiter mass planet with the JUPITER and FARGOCA codes. We implemented a radiative transfer module into the JUPITER code to account for realistic heating and cooling of the gas. We focus on the circumplanetary gas flow, determining its characteristics at very high resolution ($80\%$ of Jupiter's diameter). In our nominal simulation where the temperature evolves freely by the radiative module and reaches 13000 K at the planet, a circumplanetary envelope was formed filling the entire Roche-lobe. Because of our equation of state is simplified and probably overestimates the temperature, we also performed simulations with limited maximal temperatures in the planet region (1000 K, 1500 K, and 2000 K). In these fixed temperature cases circumplanetary disks (CPDs) were formed. This suggests that the capability to form a circumplanetary disk is not simply linked to the mass of the planet and its ability to open a gap. Instead, the gas temperature at the planet's location, which depends on its accretion history, plays also fundamental role. The CPDs in the simulations are hot and cooling very slowly, they have very steep temperature and density profiles, and are strongly sub-Keplerian. Moreover, the CPDs are fed by a strong vertical influx, which shocks on the CPD surfaces creating a hot and luminous shock-front. In contrast, the pressure supported circumplanetary envelope is characterized by internal convection and almost stalled rotation.

The identification of on-going planet formation requires the finest angular resolutions and deepest sensitivities in observations inspired by state-of-the-art numerical simulations. Hydrodynamic simulations of planet-disk interactions predict the formation of circumplanetary disks (CPDs) around accreting planetary cores. These CPDs have eluded unequivocal detection -their identification requires predictions in CPD tracers. In this work, we aim to assess the observability of embedded CPDs with ALMA as features imprinted in the gas kinematics. We use 3D Smooth Particle Hydrodynamic (SPH) simulations of CPDs around 1 and 5 M_Jup planets at large stellocentric radii, in locally isothermal and adiabatic disks. The simulations are then connected with 3D radiative transfer for predictions in CO isotopologues. Observability is assessed by corrupting with realistic long baseline phase noise extracted from the recent HL Tau ALMA data. We find that the presence of a CPD produces distinct signposts: 1) compact emission separated in velocity from the overall circumstellar disk's Keplerian pattern, 2) a strong impact on the velocity pattern when the Doppler shifted line emission sweeps across the CPD location, and 3) a local increase in the velocity dispersion. We test our predictions with a simulation tailored for HD 100546 -which has a reported protoplanet candidate. We find that the CPDs are detectable in all 3 signposts with ALMA Cycle 3 capabilities for both 1 and 5 M_Jup protoplanets, when embedded in an isothermal disk.

Planetary systems are born in the disks of gas, dust and rocky fragments that surround newly formed stars. Solid content assembles into ever-larger rocky fragments that eventually become planetary embryos. These then continue their growth by accreting leftover material in the disc. Concurrently, tidal effects in the disc cause a radial drift in the embryo orbits, a process known as migration. Fast inward migration is predicted by theory for embryos smaller than three to five Earth masses. With only inward migration, these embryos can only rarely become giant planets located at Earth's distance from the Sun and beyond, in contrast with observations. Here we report that asymmetries in the temperature rise associated with accreting infalling material produce a force (which gives rise to an effect that we call "heating torque") that counteracts inward migration. This provides a channel for the formation of giant planets and also explains the strong planet-metallicity correlation found between the incidence of giant planets and the heavy-element abundance of the host stars.

In the core-accretion model the nominal runaway gas-accretion phase brings most planets to multiple Jupiter masses. However, known giant planets are predominantly Jupiter-mass bodies. Obtaining longer timescales for gas accretion may require using realistic equations of states, or accounting for the dynamics of the circumplanetary disk (CPD) in low-viscosity regime, or both. Here we explore the second way using global, three-dimensional isothermal hydrodynamical simulations with 8 levels of nested grids around the planet. In our simulations the vertical inflow from the circumstellar disk (CSD) to the CPD determines the shape of the CPD and its accretion rate. Even without prescribed viscosity Jupiter's mass-doubling time is $\sim 10^4$ years, assuming the planet at 5.2 AU and a Minimum Mass Solar Nebula. However, we show that this high accretion rate is due to resolution-dependent numerical viscosity. Furthermore, we consider the scenario of a layered CSD, viscous only in its surface layer, and an inviscid CPD. We identify two planet-accretion mechanisms that are independent of the viscosity in the CPD: (i) the polar inflow -- defined as a part of the vertical inflow with a centrifugal radius smaller than 2 Jupiter-radii and (ii) the torque exerted by the star on the CPD. In the limit of zero effective viscosity, these two mechanisms would produce an accretion rate 40 times smaller than in the simulation.

In the last decade many kinematic groups of young stars (<100 Myr) were discovered in the solar neighbourhood. Since the most interesting period of planet formation overlaps with the age of these groups, their well dated members are attractive targets for exoplanet searches by direct imaging. We combined astrometric, photometric and X-ray data, and applied strict selection criteria to explore the stellar content of five nearby moving groups. We identified more than 100 potential new candidate members in the {\ss} Pic moving group, and in the Tucana-Horologium, Columba, Carina, and Argus associations. In order to further assess and confirm their membership status, we analysed radial velocity data and lithium equivalent widths extracted from high-resolution spectra of 54 candidate stars. We identified 35 new probable/possible young moving group members: 4 in the {\ss} Pic moving group, 11 in the Columba association, 16 in the Carina association, and 4 in the Argus association. We found serendipitiously a new AB Dor moving group member as well. For four Columba systems Hipparcos based parallaxes have already been available and as they are consistent with the predicted kinematic parallaxes, they can be considered as secure new members.

Recently, a new planet candidate was discovered on direct images around the young (10-17 Myr) A-type star HD95086. The strong infrared excess of the system indicates that, similarly to HR8799, {\ss} Pic, and Fomalhaut, the star harbors a circumstellar disk. Aiming to study the structure and gas content of the HD95086 disk, and to investigate its possible interaction with the newly discovered planet, here we present new optical, infrared and millimeter observations. We detected no CO emission, excluding the possibility of an evolved gaseous primordial disk. Simple blackbody modeling of the spectral energy distribution suggests the presence of two spatially separate dust belts at radial distances of 6 and 64 AU. Our resolved images obtained with the Herschel Space Observatory reveal a characteristic disk size of ~6.0x5.4 arcsec (540x490 AU) and disk inclination of ~25 degree. Assuming the same inclination for the planet candidate's orbit, its re-projected radial distance from the star is 62 AU, very close to the blackbody radius of the outer cold dust ring. The structure of the planetary system at HD95086 resembles the one around HR8799. Both systems harbor a warm inner dust belt and a broad colder outer disk and giant planet(s) between the two dusty regions. Modelling implies that the candidate planet can dynamically excite the motion of planetesimals even out to 270 AU via their secular perturbation if its orbital eccentricity is larger than about 0.4. Our analysis adds a new example to the three known systems where directly imaged planet(s) and debris disks co-exist.

The ASTEP (Antarctica Search for Transiting ExoPlanets) program was originally aimed at probing the quality of the Dome C, Antarctica for the discovery and characterization of exoplanets by photometry. In the first year of operation of the 40 cm ASTEP 400 telescope (austral winter 2010), we targeted the known transiting planet WASP-19b in order to try to detect its secondary transits in the visible. This is made possible by the excellent sub-millimagnitude precision of the binned data. The WASP-19 system was observed during 24 nights in May 2010. The photometric variability level due to starspots is about 1.8% (peak-to-peak), in line with the SuperWASP data from 2007 (1.4%) and larger than in 2008 (0.07%). We find a rotation period of WASP-19 of 10.7 +/- 0.5 days, in agreement with the SuperWASP determination of 10.5 +/- 0.2 days. Theoretical models show that this can only be explained if tidal dissipation in the star is weak, i.e. the tidal dissipation factor Q'star > 3.10^7. Separately, we find evidence for a secondary eclipse of depth 390 +/- 190 ppm with a 2.0 sigma significance, a phase consistent with a circular orbit and a 3% false positive probability. Given the wavelength range of the observations (420 to 950 nm), the secondary transit depth translates into a day side brightness temperature of 2690(-220/+150) K, in line with measurements in the z' and K bands. The day side emission observed in the visible could be due either to thermal emission of an extremely hot day side with very little redistribution of heat to the night side, or to direct reflection of stellar light with a maximum geometrical albedo Ag=0.27 +/- 0.13. We also report a low-frequency oscillation well in phase at the planet orbital period, but with a lower-limit amplitude that could not be attributed to the planet phase alone, and possibly contaminated with residual lightcurve trends.

The ASTEP project aims at detecting and characterizing transiting planets from Dome C, Antarctica, and qualifying this site for photometry in the visible. The first phase of the project, ASTEP South, is a fixed 10 cm diameter instrument pointing continuously towards the celestial South pole. Observations were made almost continuously during 4 winters, from 2008 to 2011. The point-to-point RMS of 1-day photometric lightcurves can be explained by a combination of expected statistical noises, dominated by the photon noise up to magnitude 14. This RMS is large, from 2.5 mmag at R=8 to 6% at R=14, because of the small size of ASTEP South and the short exposure time (30 s). Statistical noises should be considerably reduced using the large amount of collected data. A 9.9-day period eclipsing binary is detected, with a magnitude R=9.85. The 2-season lightcurve folded in phase and binned into 1000 points has a RMS of 1.09 mmag, for an expected photon noise of 0.29 mmag. The use of the 4 seasons of data with a better detrending algorithm should yield a sub-millimagnitude precision for this folded lightcurve. Radial velocity follow-up observations are conducted and reveal a F-M binary system. The detection of this 9.9-day period system with a small instrument such as ASTEP South and the precision of the folded lightcurve show the quality of Dome C for continuous photometric observations, and its potential for the detection of planets with orbital period longer than those usually detected from the ground.

Mid-infrared atomic and ionic line ratios measured in spectra of pre-main sequence stars are sensitive indicators of the hardness of the radiation field impinging on the disk surface. We present a low-resolution Spitzer IRS search for [Ar II] at 6.98 $\mu$m, [Ne II] at 12.81 $\mu$m, and [Ne III] 15.55 $\mu$m lines in 56 transitional disks. These objects, characterized by reduced near-infrared but strong far-infrared excess emission, are ideal targets to set constraints on the stellar radiation field onto the disk because their spectra are not contaminated by shock emission from jets/outflows or by molecular emission lines. After demonstrating that we can detect [Ne II] lines and recover their fluxes from the low-resolution spectra, here we report the first detections of [Ar II] lines towards protoplanetary disks. We did not detect [Ne III] emission in any of our sources. Our [Ne II]/[Ne III] line flux ratios combined with literature data suggest that a soft-EUV or X-ray spectrum produces these gas lines. Furthermore, the [Ar II]/[Ne II] line flux ratios point to a soft X-ray and/or soft-EUV stellar spectrum as the ionization source of the [Ar II] and [Ne II] emitting layer of the disk. If the soft X-ray component dominates over the EUV than we would expect larger photoevaporation rates hence a reduction of the time available to form planets.