Search

We present 3 years of high-contrast imaging of the PDS 70 b and c accreting protoplanets with the new extreme AO system MagAO-X as part of the MaxProtoPlanetS survey of H$\alpha$ protoplanets. In 2023 and 2024 our sharp (25-27 mas FWHM); well AO corrected (20-26% Strehl), deep (2-3.6hr) images detect compact (r~30 mas; r~3 au) circumplanetary disks (CPDs) surrounding both protoplanets. Starlight scattering off the dusty outer edges of these CPDs is the likely source of the bright compact continuum light detected within ~30 mas of both planets in our simultaneously obtained continuum 668 nm filter images. After subtraction of contaminating continuum and PSF residuals with pyKLIP ADI and SDI we obtained high-contrast ASDI H$\alpha$ images of both planets in 2022, 2023 and 2024. We find the H$\alpha$ line flux of planet b fell by (8.1$\pm$1.6)x10$^{-16}$ ergs/s/cm$^2$ a factor of 4.6 drop in flux from 2022 to 2023. In March 2024, planet b continued to be faint with just a slight 1.6x rise to an H$\alpha$ line flux of (3.64$\pm$0.87)x10$^{-16}$ ergs/s/cm$^2$. For c we measure a significant increase of (2.74$\pm$0.51)x10$^{-16}$ ergs/s/cm$^2$ from 2023 to 2024 which is a factor of 2.3x increase. So both protoplanets have recently experienced significant H$\alpha$ variability with ~1 yr sampling. In 2024, planet c is brighter than b: as c is brightening and b generally fading. We also tentatively detect one new point source "CC3" inside the inner disk (~49 mas; at PA~295 deg; 2024) with orbital motion roughly consistent with a ~5.6 au orbit., Comment: 50 pages, 20 figures, published in January 2025 Astronomical Journal. NOTE Figure D1 has correct planetary mass accretion scale compared to published AJ version

We carried out very deep VLT/SPHERE imaging polarimetry of the nearby system Eps Eri based on 38.5 hours of integration time with a 600 - 900 nm broadband filter to search for polarized scattered light from a planet or from circumstellar dust using AO, coronagraphy, high precision differential polarimetry, and angular differential imaging. We have improved several data reduction and post-processing techniques and also developed new ones to further increase the sensitivity of SPHERE/ZIMPOL. The data provide unprecedented contrast limits, but no significant detection of a point source or an extended signal from circumstellar dust. For each observing epoch, we obtained a point source contrast for the polarized intensity between $2\cdot 10^{-8}$ and $4\cdot 10^{-8}$ at the expected separation of the planet Eps Eri b of 1'' near quadrature phase. The polarimetric contrast limits are about six to 50 times better than the intensity limits because polarimetric imaging is much more efficient in speckle suppression. Combining the entire 14-month data set to the search for a planet moving on a Keplerian orbit with the K-Stacker software further improves the contrast limits by a factor of about two, to about $8 \cdot 10^{-9}$ at 1''. This would allow the detection of a planet with a radius of about 2.5 Jupiter radii. The surface brightness contrast limits achieved for the polarized intensity from an extended scattering region are about 15 mag arcsec$^{-2}$ at 1'', or up to 3 mag arcsec$^{-2}$ deeper than previous limits. For Eps Eri, these limits exclude the presence of a narrow dust ring and they constrain the dust properties. This study shows that the polarimetric contrast limits for reflecting planets with SPHERE/ZIMPOL can be improved to a level $<10^{-8}$ simply by collecting more data over many nights and using the K-Stacker software., Comment: 26 pages, 20 figures, 5 tables. Accepted for publication in A&A

The sample of planet-forming disks observed by high-contrast imaging campaigns over the last decade is mature enough to enable the demographical analysis of individual star-forming regions. We present the full census of Taurus sources with VLT/SPHERE polarimetric images available. The whole sample sums up to 43 targets (of which 31 have not been previously published) corresponding to one-fifth of the Class II population in Taurus and about half of such objects that are observable. A large fraction of the sample is apparently made up of isolated faint disks (equally divided between small and large self-shadowed disks). Ambient signal is visible in about one-third of the sample. This probes the interaction with the environment and with companions or the outflow activity of the system. The central portion of the Taurus region almost exclusively hosts faint disks, while the periphery also hosts bright disks interacting with their surroundings. The few bright disks are found around apparently older stars. The overall picture is that the Taurus region is in an early evolutionary stage of planet formation. Yet, some objects are discussed individually, as in an intermediate or exceptional stage of the disk evolution. This census provides a first benchmark for the comparison of the disk populations in different star forming regions., Comment: Accepted by A&A

(Pre-)transitional disks show gaps and cavities that can be related with on-going planet formation. According to theory, young embedded planets can accrete material from the circumplanetary and circumstellar disks, so that they could be detected in accretion tracers, like the H$_{\alpha}$ emission line. In this work, we present spectral angular differential imaging AO-assisted observations of five (pre-)transitional disks obtained with SPHERE/ZIMPOL at the Very Large Telescope (VLT). They were obtained in the H$_{\alpha}$ line and the adjacent continuum. We have combined spectral and angular differential imaging techniques to increase the contrast in the innermost regions close to the star, and search for the signature of young accreting protoplanets. As a result, the reduced images do not show any clear H$_{\alpha}$ point source around any of the targets. We report faint H$_{\alpha}$ emissions around TW Hya and HD163296: while the former is most probably an artifact related with a spike, the nature of the latter remains unclear. The spectral and angular differential images yield contrasts of 6--8 magnitudes at separations of $\sim$ 100 mas from the central stars, except in the case of LkCa15, with values of $\sim$3 mag. We have estimated upper limits to the accretion luminosity of potential protoplanets, obtaining that planetary models provide an average value of $L_{\rm acc} \sim 10^{-4}$ $L_{\odot}$ at 200 mas, which is $\sim$2 orders of magnitude higher than the $L_{\rm acc}$ estimated from the extrapolation of the $L_{H_{\alpha}}$ - $L_{acc}$ stellar relationship. We explain the lack of protoplanet detections as a combination of different factors, like e.g. episodic accretion, extinction from the circumstellar and circumplanetray disks, and/or a majority of low-mass, low-accreting planets., Comment: 17 pages, Accepted for publication in A&A

Context. HD 95086 is a young nearby Solar System analog hosting a giant exoplanet orbiting at 57 au from the star between an inner and outer debris belt. The existence of additional planets has been suggested as the mechanism that maintains the broad cavity between the two belts. Aims. We present a dedicated monitoring of HD 95086 with the VLT/SPHERE instrument to refine the orbital and atmospheric properties of HD 95086 b, and to search for additional planets in this system. Methods. SPHERE observations, spread over ten epochs from 2015 to 2019 and including five new datasets, were used. Combined with archival observations, from VLT/NaCo (2012-2013) and Gemini/GPI (2013-2016), the extended set of astrometric measurements allowed us to refine the orbital properties of HD 95086 b. We also investigated the spectral properties and the presence of a circumplanetary disk around HD 95086 b by using the special fitting tool exploring the diversity of several atmospheric models. In addition, we improved our detection limits in order to search for a putative planet c via the K-Stacker algorithm. Results. We extracted for the first time the JH low-resolution spectrum of HD 95086 b by stacking the six best epochs, and confirm its very red spectral energy distribution. Combined with additional datasets from GPI and NaCo, our analysis indicates that this very red color can be explained by the presence of a circumplanetary disk around planet b, with a range of high-temperature solutions (1400-1600 K) and significant extinction (Av > 10 mag), or by a super-solar metallicity atmosphere with lower temperatures (800-1300 K), and small to medium amount of extinction (Av < 10 mag). We do not find any robust candidates for planet c, but give updated constraints on its potential mass and location., Comment: 29 pages, 20 figures, A&A, accepted

[Abridged] Most disks observed at high angular resolution show substructures. Knowledge about the gas surface density and temperature is essential to understand these. The aim of this work is to constrain the gas temperature and surface density in two transition disks: LkCa15 and HD 169142. We use new ALMA observations of the $^{13}$CO $J=6-5$ transition together with archival $J=2-1$ data of $^{12}$CO, $^{13}$CO and C$^{18}$O to observationally constrain the gas temperature and surface density. Furthermore, we use the thermochemical code DALI to model the temperature and density structure of a typical transition disk. The $6-5/2-1$ line ratio in LkCa15 constrains the gas temperature in the emitting layers inside the dust cavity to be up to 65 K, warmer than in the outer disk at 20-30 K. For the HD 169142, the peak brightness temperature constrains the gas in the dust cavity of HD 169142 to be 170 K, whereas that in the outer disk is only 100 K. Models also show that a more luminous central star, a lower abundance of PAHs and the absence of a dusty inner disk increase the temperature of the emitting layers and hence the line ratio in the gas cavity. The gas column density in the LkCa15 dust cavity drops by a factor >2 compared to the outer disk, with an additional drop of an order of magnitude inside the gas cavity at 10 AU. In the case of HD 169142, the gas column density drops by a factor of 200$-$500 inside the gas cavity, which could be due to a massive companion of several M$_{\mathrm{J}}$. The broad dust-depleted gas region from 10-68 AU for LkCa15 may imply several lower mass planets. This work demonstrates that knowledge of the gas temperature is important to determine the gas surface density and thus whether planets, and if so what kind of planets, are the most likely carving the dust cavities., Comment: Accepted for publication in Astronomy and astrophysics

We carried out 3D dust+gas radiative hydrodynamic simulations of forming planets. We investigated a parameter grid of Neptune-, Saturn-, Jupiter-, and 5 Jupiter-mass planets at 5.2, 30, 50 AU distance from their star. We found that the meridional circulation \citep{Szulagyi14,FC16} drives a strong vertical flow for the dust as well, hence the dust is not settled in the midplane, even for mm-sized grains. The meridional circulation will deliver dust and gas vertically onto the circumplanetary region, efficiently bridging over the gap. The Hill-sphere accretion rates for the dust are $\sim10^{-8}$ to $10^{-10}$ $\rm{M_{Jup}/yr}$, increasing with planet-mass. For the gas component, the gain is $10^{-6}$ to $10^{-8}$ $\rm{M_{Jup}/yr}$. The difference between the dust and gas accretion rates is smaller with decreasing planetary mass. In the vicinity of the planet, the mm-grains can get trapped easier than the gas, which means the circumplanetary disk might be enriched with solids in comparison to the circumstellar disk. We calculated the local dust-to-gas ratio (DTG) everywhere in the circumstellar disk and identified the altitude above the midplane where the DTG is 1, 0.1, 0.01, 0.001. The larger the planetary mass, the higher the mm-sized dust is delivered and a larger fraction of the dust disk is lifted by the planet. The stirring of mm-dust is negligible for Neptune-mass planets or below, but significant above Saturn-mass., Comment: ApJ accepted

The detection of a wide range of substructures such as rings, cavities and spirals has become a common outcome of high spatial resolution imaging of protoplanetary disks, both in the near-infrared scattered light and in the thermal millimetre continuum emission. The most frequent interpretation of their origin is the presence of planetary-mass companions perturbing the gas and dust distribution in the disk (perturbers), but so far the only bona-fide detection has been the two giant planets around PDS 70. Here, we collect a sample of 15 protoplanetary disks showing substructures in SPHERE scattered light images and present a homogeneous derivation of planet detection limits in these systems. We also estimate the mass of these perturbers through a Hill radius prescription and a comparison to ALMA data. Assuming that one single planet carves each substructure in scattered light, we find that more massive perturbers are needed to create gaps within cavities than rings, and that we might be close to a detection in the cavities of RX J1604, RX J1615, Sz Cha, HD 135344B and HD 34282. We reach typical mass limits in these cavities of 3-10 Mjup. For planets in the gaps between rings, we find that the detection limits of SPHERE are about an order of magnitude away in mass, and that the gaps of PDS 66 and HD 97048 seem to be the most promising structures for planet searches. The proposed presence of massive planets causing spiral features in HD 135344B and HD 36112 are also within SPHERE's reach assuming hot-start models.These results suggest that current detection limits are able to detect hot-start planets in cavities, under the assumption that they are formed by a single perturber located at the centre of the cavity. More realistic planet mass constraints would help to clarify whether this is actually the case, which might point to perturbers not being the only way of creating substructures., Comment: 26 pages, 33 figures, accepted for publication in A&A

One of the long-term goals of exoplanet science is the atmospheric characterization of dozens of small exoplanets in order to understand their diversity and search for habitable worlds and potential biosignatures. Achieving this goal requires a space mission of sufficient scale. We seek to quantify the exoplanet detection performance of a space-based mid-infrared nulling interferometer that measures the thermal emission of exoplanets. For this, we have developed an instrument simulator that considers all major astrophysical noise sources and coupled it with Monte Carlo simulations of a synthetic exoplanet population around main-sequence stars within 20 pc. This allows us to quantify the number (and types) of exoplanets that our mission concept could detect over a certain time period. Two different scenarios to distribute the observing time among the stellar targets are discussed and different apertures sizes and wavelength ranges are considered. Within a 2.5-year initial search phase, an interferometer consisting of four 2 m apertures with a total instrument throughput of 5% covering a wavelength range between 4 and 18.5 $\mu$m could detect up to ~550 exoplanets with radii between 0.5 and 6 R$_\oplus$ with an integrated SNR$\ge$7. At least ~160 of the detected exoplanets have radii $\le$1.5 R$_\oplus$. Depending on the observing scenario, ~25-45 rocky exoplanets (objects with radii between 0.5 and 1.5 $_{\oplus}$) orbiting within the empirical habitable zone (eHZ) of their host stars are among the detections. With an aperture size of 3.5 m, the total number of detections can increase to up to ~770, including ~60-80 rocky, eHZ planets. With 1 m aperture size, the maximum detection yield is ~315 exoplanets, including $\le$20 rocky, eHZ planets. In terms of predicted detection yield, such a mission can compete with large single-aperture reflected light missions. (abridged), Comment: Accepted for publication by A&A - some typos corrected and affiliations updated; 14 pages main text (incl. 14 figures); first paper in the LIFE paper series; papers II (arXiv:2203.00471) and III (arXiv:2112.02054) are also available

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., Comment: Accepted for publication in A&A

Proxima Centauri is known to host an earth-like planet in its habitable zone; very recently a second candidate planet was proposed based on radial velocities. At quadrature, the expected projected separation of this new candidate is larger than 1 arcsec, making it a potentially interesting target for direct imaging. While difficult, identification of the optical counterpart of this planet would allow detailed characterization of the closest planetary system. We searched for a counterpart in SPHERE images acquired during four years through the SHINE survey. In order to account for the large orbital motion of the planet, we used a method that assumes the circular orbit obtained from radial velocities and exploits the sequence of observations acquired close to quadrature in the orbit. We checked this with a more general approach that considers keplerian motion, K-stacker. We did not obtain a clear detection. The best candidate has S/N=6.1 in the combined image. A statistical test suggests that the probability that this detection is due to random fluctuation of noise is < 1% but this result depends on the assumption that distribution of noise is uniform over the image. The position of this candidate and the orientation of its orbital plane fit well with observations in the ALMA 12m array image. However, the astrometric signal expected from the orbit of the candidate we detected is 3-sigma away from the astrometric motion of Proxima as measured from early Gaia data. This, together with the unexpectedly high flux associated with our direct imaging detection, means we cannot confirm that our candidate is indeed Proxima c. On the other hand, if confirmed, this would be the first observation in imaging of a planet discovered from radial velocities and the second one (after Fomalhaut b) of reflecting circumplanetary material. Further confirmation observations should be done as soon as possible., Comment: 14 pages, 5 figures. Accepted for publication on Astronomy and Astrophysics

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 assuming boundary layer accretion. 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., Comment: ApJ accepted. 13 pages, 5 figures, 4 tables. This version: typo corrections

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 formed by gravitational instability at 50\,AU from its star. The disc parameters were from a 3D global SPH simulation. We used a population synthesis approach, where we placed satellite embryos in this disc, and let them accrete mass, migrate, collide until the gaseous disc is dissipated. In each run we changed the initial dust-to-gas ratio, dispersion- and refilling time-scales within reasonable limits, as well as the number of embryos and their starting locations. We found that most satellites have mass similar to the Galilean ones, but very few can reach a maximum of 3 $M_{\rm{Earth}}$ due to the massive circumplanetary disc. Large moons are often form as far as 0.5 $R_{\rm{disc}}$. The migration rate of satellites are fast, hence during the disc lifetime, an average of 10 $M_{\rm{Earth}}$ worth of moons will be engulfed by the planet, increasing greatly its metallicity. We also investigated the effect of the planet's semi-major axis on the resulting satellite systems by re-scaling our model. This test revealed that for the discs closer to the star, the formed moons are lighter, and a larger amount of satellites are lost into the planet due to the even faster migration. Finally, we checked the probability of detecting satellites like our population, which resulted in a low number of $\leq$ 3\% even with upcoming powerful telescopes like E-ELT., Comment: 16 pages, 24 figures, accepted for publication by MNRAS

RefPlanets is a guaranteed time observation (GTO) programme that uses the Zurich IMaging POLarimeter (ZIMPOL) of SPHERE/VLT for a blind search for exoplanets in wavelengths from 600-900 nm. The goals of this study are the characterization of the unprecedented high polarimetic contrast and polarimetric precision capabilities of ZIMPOL for bright targets, the search for polarized reflected light around some of the closest bright stars to the Sun and potentially the direct detection of an evolved cold exoplanet for the first time. For our observations of Alpha Cen A and B, Sirius A, Altair, Eps Eri and Tau Ceti we used the polarimetric differential imaging (PDI) mode of ZIMPOL which removes the speckle noise down to the photon noise limit for angular separations >0.6". We describe some of the instrumental effects that dominate the noise for smaller separations and explain how to remove these additional noise effects in post-processing. We then combine PDI with angular differential imaging (ADI) as a final layer of post-processing to further improve the contrast limits of our data at these separations. For good observing conditions we achieve polarimetric contrast limits of 15.0-16.3 mag at the effective inner working angle of about 0.13", 16.3-18.3 mag at 0.5" and 18.8-20.4 mag at 1.5". The contrast limits closer in (<0.6") depend significantly on the observing conditions, while in the photon noise dominated regime (>0.6"), the limits mainly depend on the brightness of the star and the total integration time. We compare our results with contrast limits from other surveys and review the exoplanet detection limits obtained with different detection methods. For all our targets we achieve unprecedented contrast limits. Despite the high polarimetric contrasts we are not able to find any additional companions or extended polarized light sources in the data that has been taken so far., Comment: 23 pages, 17 figures, 2 table, Accepted for publication in A&A

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., Comment: 14 pages, 10 figures, accepted by A&A

There are growing amount of very high-resolution polarized scattered light images of circumstellar disks. Nascent giant planets planets are surrounded by their own circumplanetary disks which may scatter and polarize both the planetary and stellar light. Here we investigate whether we could detect circumplanetary disks with the same technique and what can we learn from such detections. Here we created scattered light mock observations at 1.245 microns (J band) for instruments like SPHERE and GPI, for various planetary masses (0.3, 1.0, 5.0, 10.0 $\rm{M_{Jup}}$), disk inclinations (90, 60, 30, 0 degrees) and planet position angles (0, 45, 90 degrees). We found that the detection of a circumplanetary disk at 50AU from the star is significantly favored if the planet is massive ($\geq 5 \rm{M_{Jup}}$) and the system is nearly face-on ($\leq 30^\circ$). In these cases the accretion shock front on the surface of the circumplanetary disks are strong and bright enough to help the visibility of this subdisk. 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 is a viable tool to pinpoint the presence of the circumplanetary disk within the circumstellar disk, as the two disks are behaving differently on those images., Comment: Accepted at MNRAS

Thanks to recent high resolution ALMA observations, there is an accumulating evidence for presence of giant planets with masses from $\sim 0.01$ Jupiter mass to a few Jupiter mass with separations up to $ 100$~AU in the annular structures observed in young protoplanetary discs. We point out that these observations set unique "live" constraints on the process of gas accretion onto sub-Jovian planets that were not previously available. Accordingly, we use a population synthesis approach in a new way: we build time-resolved models and compare the properties of the synthetic planets with the ALMA data at the same age. Applying the widely used gas accretion formulae leads to a deficit of sub-Jovian planets and an over-abundance of a few Jupiter mass planets compared to observations. We find that gas accretion rate onto planets needs to be suppressed by about an order of magnitude to match the observed planet mass function. This slower gas giant growth predicts that the planet mass should correlate positively with the age of the protoplanetary disc, albeit with a large scatter. This effect is not clearly present in the ALMA data but may be confirmed in the near future with more observations., Comment: Accepted to MNRAS Letters

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., Comment: Version after the 2nd referee report

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., Comment: Accepted for publication in A&A (20 pages, 18 figures)

For a comprehensive understanding of planetary formation and evolution, we need to investigate the environment in which planets form: circumstellar disks. Here we present high-contrast imaging observations of V4046 Sagittarii, a 20-Myr-old close binary known to host a circumbinary disk. We have discovered the presence of rotating shadows in the disk, caused by mutual occultations of the central binary. Shadow-like features are often observed in disks\cite{garufi,marino15}, but those found thus far have not been due to eclipsing phenomena. We have used the phase difference due to light travel time to measure the flaring of the disk and the geometrical distance of the system. We calculate a distance that is in very good agreement with the value obtained from the Gaia mission's Data Release 2 (DR2), and flaring angles of $\alpha = 6.2 \pm 0.6 $ deg and $\alpha = 8.5 \pm 1.0 $ deg for the inner and outer disk rings, respectively. Our technique opens up a path to explore other binary systems, providing an independent estimate of distance and the flaring angle, a crucial parameter for disk modelling., Comment: Authors' version of the article published in Nature Astronomy. Final version available at https://www.nature.com/natastron/

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., Comment: 14 pages, accepted by ApJ

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., Comment: In press at ApJL

With an orbital distance comparable to that of Saturn in the solar system, \bpic b is the closest (semi-major axis $\simeq$\,9\,au) exoplanet that has been imaged to orbit a star. Thus it offers unique opportunities for detailed studies of its orbital, physical, and atmospheric properties, and of disk-planet interactions. With the exception of the discovery observations in 2003 with NaCo at the Very Large Telescope (VLT), all following astrometric measurements relative to \bpic have been obtained in the southwestern part of the orbit, which severely limits the determination of the planet's orbital parameters. We aimed at further constraining \bpic b orbital properties using more data, and, in particular, data taken in the northeastern part of the orbit. We used SPHERE at the VLT to precisely monitor the orbital motion of beta \bpic b since first light of the instrument in 2014. We were able to monitor the planet until November 2016, when its angular separation became too small (125 mas, i.e., 1.6\,au) and prevented further detection. We redetected \bpic b on the northeast side of the disk at a separation of 139\,mas and a PA of 30$^{\circ}$ in September 2018. The planetary orbit is now well constrained. With a semi-major axis (sma) of $a = 9.0 \pm 0.5$ au (1 $\sigma $), it definitely excludes previously reported possible long orbital periods, and excludes \bpic b as the origin of photometric variations that took place in 1981. We also refine the eccentricity and inclination of the planet. From an instrumental point of view, these data demonstrate that it is possible to detect, if they exist, young massive Jupiters that orbit at less than 2 au from a star that is 20 pc away., Comment: accepted by A&A

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., Comment: 10 pages, 10 figures, in press by A&A

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., Comment: 15 pages, 17 figures. Accepted by MNRAS, please check the final published version

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., Comment: 21 pages, 14 figures. Accepted for publication in ApJ

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., Comment: Accepted for publication at MNRAS. Typos are corrected since previous version. 11 pages, 5 tables, 4 figures

We present observations of the known edge-on debris disk around HIP 79977 (HD 146897, F star in Upper Sco, 123 pc), taken with the ZIMPOL differential polarimeter of the SPHERE instrument in the Very Broad Band filter ($\lambda_c=735$ nm, $\Delta\lambda=290$ nm) with a spatial resolution of about 25 mas. We measure the polarization flux along and perpendicular to the disk spine of the highly inclined disk for projected separations between 0.2" (25 AU) and 1.6" (200 AU) and investigate the diagnostic potential of such data with model simulations. The polarized flux contrast ratio for the disk is $F_{pol}/F_\ast= (5.5 \pm 0.9) 10^{-4}$. The surface brightness reaches a maximum of 16.2 mag arcsec$^{-2}$ at a separation of $0.2"-0.5"$ along the disk spine with a maximum surface brightness contrast of 7.64 mag arcsec$^{-2}$. The polarized flux has a minimum near the star $<0.2"$ because no or only little polarization is produced by forward or backward scattering in the disk section lying in front of or behind the star. The data are modeled as a circular dust belt with an inclination $i=85(\pm 1.5)^\circ$ and a radius between $r_0$ = 60 AU and 90 AU. The radial density dependence is described by $(r/r_0)^{\alpha}$ with a steep power law index $\alpha=5$ inside $r_0$ and a more shallow index $\alpha=-2.5$ outside $r_0$. The scattering asymmetry factor lies between $g$ = 0.2 and 0.6 adopting a scattering angle-dependence for the fractional polarization as for Rayleigh scattering. Our data are qualitatively very similar to the case of AU Mic and they confirm that edge-on debris disks have a polarization minimum at a position near the star and a maximum near the projected separation of the main debris belt. The comparison of the polarized flux contrast ratio $F_{pol}/F_{\ast}$ with the fractional infrared excess provides strong constraints on the scattering albedo of the dust., Comment: Corrected typos, changed affiliation

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., Comment: 10 pages, 8 figures, 4 tables, accepted for publication in A&A, file updated

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., Comment: 11 pages, 7 figures, accepted for publication in AJ

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}}$., Comment: Accepted for publication at ApJ

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"., Comment: 17 pages, 14 figures, accepted for publication in A&A. Revised version after language editing

R Aqr is a symbiotic binary system consisting of a mira variable, a hot companion with a spectacular jet outflow, and an extended emission line nebula. We have used R Aqr as test target for the visual camera subsystem ZIMPOL, which is part of the new extreme adaptive optics (AO) instrument SPHERE at the Very Large Telescope (VLT). We compare our observations with data from the Hubble Space Telescope (HST) and illustrate the complementarity of the two instruments. We determine from the Halpha emission the position, size, geometric structure, and line fluxes of the jet source and the clouds in the innermost region (<2") of R Aqr and determine Halpha emissivities mean density, mass, recombination time scale, and other cloud parameters. Our data resolve for the first time the R Aqr binary and we measure for the jet source a relative position 46+/-1 mas West of the mira. The central jet source is the strongest Halpha component. North east and south west from the central source there are many clouds with very diverse structures. We see in the SW a string of bright clouds arranged in a zig-zag pattern and, further out, more extended bubbles. In the N and NE we see a bright, very elongated filamentary structure and faint perpendicular "wisps" further out. Some jet clouds are also detected in the ZIMPOL [OI] and He I filters, as well as in the HST line filters for Halpha, [OIII], [NII], and [OI]. We determine jet cloud parameters and find a very well defined anti-correlation between cloud density and distance to the central binary. Future Halpha observations will provide the orientation of the orbital plane of the binary and allow detailed hydrodynamical investigations of this jet outflow and its interaction with the wind of the red giant companion., Comment: 24 pages, 14 figures (accepted for publication in Astronomy and Astrophysics)

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., Comment: 10 pages, 6 figures, SPIE conference proceedings, Ground-based and Airborne Instrumentation for Astronomy VI, 2016

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., Comment: 12 pages, 9 figures, 1 table, accepted for publication at MNRAS

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., Comment: 5 pages, 3 figures, 1 table, accepted for publication at MNRAS Letters

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., Comment: Accepted for publication at MNRAS

For the very first time we present 3D simulations of planets embedded in stellar irradiated discs. It is well known that thermal effects could reverse the direction of planetary migration from inwards to outwards, potentially saving planets in the inner, optically thick parts of the protoplanetary disc. When considering stellar irradiation in addition to viscous friction as a source of heating, the outer disc changes from a shadowed to a flared structure. Using a suited analytical formula it has been shown that in the flared part of the disc the migration is inwards; planets can migrate outwards only in shadowed regions of the disc, { because the radial gradient of entropy is stronger there}. In order to confirm this result numerically, we have computed the total torque acting on planets held on fixed orbits embedded in stellar irradiated 3D discs using the hydrodynamical code FARGOCA. We find qualitatively good agreement between the total torque obtained with numerical simulations and the one predicted by the analytical formula. For large masses (>20 Earth masses) we find quantitative agreement, and we obtain outwards migration regions for planets up to 60 Earth masses in the early stages of accretional discs. We find nevertheless that the agreement with the analytic formula is quite fortuitous because the formula underestimates the size of the horseshoe region; this error is compensated by imperfect estimates of other terms, most likely the cooling rate and saturation., Comment: 14 pages , 10 figures

We examine the gas circulation near a gap opened by a giant planet in a protoplanetary disk. We show with high resolution 3D simulations that the gas flows into the gap at high altitude over the mid-plane, at a rate dependent on viscosity. We explain this observation with a simple conceptual model. From this model we derive an estimate of the amount of gas flowing into a gap opened by a planet with Hill radius comparable to the scale-height of a layered disk (i. e. a disk with viscous upper layer and inviscid midplane). Our estimate agrees with modern MRI simulations(Gressel et al., 2013). We conclude that gap opening in a layered disk can not slow down significantly the runaway gas accretion of Saturn to Jupiter-mass planets., Comment: in press as a Note in Icarus

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., Comment: 11 pages, 10 figures. Accepted for publication in Astrophysical Journal

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., Comment: Accepted for publication in ApJ Letters, 6 pages, 2 figures, 3 tables (emulateapj style)

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., Comment: Accepted for publication in MNRAS, 14 pages, 3 figures, 3 tables

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., Comment: Accepted for publication in Astronomy and Astrophysics, 13 pages, 13 figures

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., Comment: 5 pages, 3 figures, IAUS 288 conference proceedings

Due to the strong effect of systematics/trends in variable star observations, we employ the Trend Filtering Algorithm (TFA) on a subset of the MACHO database and search for variable stars. TFA has been applied successfully in planetary transit searches, where weak, short-lasting periodic dimmings are sought in the presence of noise and various systematics (due to, e.g., imperfect flat fielding, crowding, etc). These latter effects introduce colored noise in the photometric time series that can lead to a complete miss of the signal. By using a large number of available photometric time series of a given field, TFA utilizes the fact that the same types of systematics appear in several/many time series of the same field. As a result, we fit each target time series by a (least-square-sense) optimum linear combination of templates and frequency-analyze the residuals. Once a signal is found, we reconstruct the signal by employing the full model, including the signal, systematics and noise. We apply TFA on the brightest ~5300 objects from subsets of each of the MACHO Large Magellanic Cloud fields #1 and #79. We find that the Fourier frequency analysis performed on the original data detect some 60% of the objects as trend-dominated. This figure decreases essentially to zero after using TFA. Altogether, We detect 387 variables in the two fields, 183 of which would have remained undetected without using TFA. Where possible, we give preliminary classification of the variables found., Comment: 12 pages, 15 figures, 3 tables with online material; to appear in Astronomy and Astrophysics