We investigate the morphology of the colour-magnitude diagram (CMD) of the open cluster NGC 2509 in comparison with other Galactic open clusters of similar age using Gaia photometry. At $\sim900$ Myr Galactic open clusters in our sample all show an extended main sequence turn off (eMSTO) with the exception of NGC 2509, which presents an exceptionally narrow CMD. Our analysis of the Gaia data rules out differential extinction, stellar density, and binaries as a cause for the singular MSTO morphology in this cluster. We interpret this feature as a consequence of the stellar rotation distribution within the cluster and present the analysis with MIST stellar evolution models that include the effect of stellar rotation on which we based our conclusion. In particular, these models point to an unusually narrow range of stellar rotation rates ($\Omega/\Omega_{\rm{crit,ZAMS}} = [0.4, 0.6]$) within the cluster as the cause of this singular feature in the CMD of NGC 2509. Interestingly, models that do not include rotation are not as good at reproducing the morphology of the observed CMD in this cluster.

We present an analysis of the relatively low mass ($\sim2400$~M$_{\odot}$), $\sim800$~Myr, Galactic open cluster, NGC~2818, using Gaia DR2 results combined with VLT/FLAMES spectroscopy. Using Gaia DR2 proper motions and parallax measurements we are able to select a clean sample of cluster members. This cluster displays a clear extended main sequence turn-off (eMSTO), a phenomenon previously studied mainly in young and intermediate age massive clusters in the Magellanic clouds. The main sequence of NGC~2818 is extremely narrow, with a width of $\sim0.01$ magnitudes (G$_{\rm BP} - $ G$_{\rm RP}$), suggesting very low levels of differential extinction. Using VLT/FLAMES spectroscopy of 60 cluster members to measure the rotational velocity of the stars (Vsini) we find that stars on the red side of the eMSTO have high Vsini ($>160$~km/s) while stars on the blue side have low Vsini ($<160$~km/s), in agreement with model predictions. The cluster also follows the previously discovered trend between the age of the cluster and the extent of the eMSTO. We conclude that stellar rotation is the likely cause of the eMSTO phenomenon.

Potential signatures of proto-planets embedded in their natal protoplanetary disk are radial gaps or cavities in the continuum emission in the IR-mm wavelength range. ALMA observations are now probing spatially resolved rotational line emission of CO and other chemical species. These observations can provide complementary information on the mechanism carving the gaps in dust and additional constraints on the purported planet mass. We post-process 2D hydrodynamical simulations of planet-disk models, where the dust densities and grain size distributions are computed with a dust evolution code. The simulations explore different planet masses ($1\,M_{\rm J}\leq M_{\rm p}\leq15\,M_{\rm J}$) and turbulent parameters. The outputs are post-processed with the thermo-chemical code DALI, accounting for the radially and vertically varying dust properties as in Facchini et al. (2017). We obtain the gas and dust temperature structures, chemical abundances, and synthetic emission maps of both thermal continuum and CO rotational lines. This is the first study combining hydro simulations, dust evolution and chemistry to predict gas emission of disks hosting massive planets. All radial intensity profiles of the CO main isotopologues show a gap at the planet location. The ratio between the location of the gap as seen in CO and the peak in the mm continuum at the pressure maximum outside the orbit of the planet shows a clear dependence on planet mass. Due to the low dust density in the gaps, the dust and gas components can become thermally decoupled, with the gas being colder than the dust. The gaps seen in CO are due to a combination of gas temperature dropping at the location of the planet, and of the underlying surface density profile. In none of the models is a CO cavity observed, only CO gaps, indicating that one single massive planet is not able to explain the CO cavities observed in transition disks.

In this paper we describe an unmanned aerial system equipped with a thermal-infrared camera and software pipeline that we have developed to monitor animal populations for conservation purposes. Taking a multi-disciplinary approach to tackle this problem, we use freely available astronomical source detection software and the associated expertise of astronomers, to efficiently and reliably detect humans and animals in aerial thermal-infrared footage. Combining this astronomical detection software with existing machine learning algorithms into a single, automated, end-to-end pipeline, we test the software using aerial video footage taken in a controlled, field-like environment. We demonstrate that the pipeline works reliably and describe how it can be used to estimate the completeness of different observational datasets to objects of a given type as a function of height, observing conditions etc. -- a crucial step in converting video footage to scientifically useful information such as the spatial distribution and density of different animal species. Finally, having demonstrated the potential utility of the system, we describe the steps we are taking to adapt the system for work in the field, in particular systematic monitoring of endangered species at National Parks around the world.

[Abridged] Regions of low ionisation where the activity of the magneto-rotational instability is suppressed, the so-called dead zones, have been suggested to explain gaps and asymmetries of transition disks. We investigate the gas and dust evolution simultaneously assuming simplified prescriptions for a dead zone and a magnetohydrodynamic (MHD) wind acting on the disk. We explore whether the resulting gas and dust distribution can create signatures similar to those observed in transition disks. For the dust evolution, we included the transport, growth, and fragmentation of dust particles. To compare with observations, we produced synthetic images in scattered optical light and in thermal emission at mm wavelengths. In all models with a dead zone, a bump in the gas surface density is produced that is able to efficiently trap large particles ($\gtrsim 1$ mm) at the outer edge of the dead zone. The gas bump reaches an amplitude of a factor of $\sim5$, which can be enhanced by the presence of an MHD wind that removes mass from the inner disk. While our 1D simulations suggest that such a structure can be present only for $\sim$1 Myr, the structure may be maintained for a longer time when more realistic 2D/3D simulations are performed. In the synthetic images, gap-like low-emission regions are seen at scattered light and in thermal emission at mm wavelengths, as previously predicted in the case of planet-disk interaction. As a conclusion, main signatures of transition disks can be reproduced by assuming a dead zone in the disk, such as gap-like structure in scattered light and millimetre continuum emission, and a lower gas surface density within the dead zone. Previous studies showed that the Rossby wave instability can also develop at the edge of such dead zones, forming vortices and also creating asymmetries.

We investigate the effect that the turbulent mixing strength parameter $\alpha_{\rm{turb}}$ plays on near-infrared polarimetric and sub-millimetre interferometric imaging observations of transitional discs (TDs) with a gap carved by a planet. We generate synthetic observations of these objects with ALMA and VLT/SPHERE-ZIMPOL by combining hydrodynamical, dust evolution, radiative transfer and instrument models for values of $\alpha_{\rm{turb}}=[10^{-4}, 10^{-3}, 10^{-2}]$. We find that, through a combination of effects on the viscosity of the gas, the turbulent mixing and dust evolution processes, $\alpha_{\rm{turb}}$ strongly affects the morphology of the dust distribution that can be traced with these observations. We constrain the value of $\alpha_{\rm{turb}}$ to be within an order of magnitude of $10^{-3}$ in TD sources that show cavities in sub-mm continuum images while featuring continuous distribution of dust or smaller cavities in NIR-polarimetric images.

Context. Some of transition disks show asymmetric structures in thermal sub-millimetre emission and optical scattered light. These structures can be the result of planet(s) or companions embedded in the disk. Aims. We aim to detect and analyse the scattered light of the transition disk J160421.7-213028, identify disk structures, and compare the results with previous observations of this disk at other wavelengths. Methods. We obtained and analysed new polarised intensity observations of the transition disk J160421.7-213028 with VLT/SPHERE using the visible light instrument ZIMPOL at $R'$ band (0.626$\mu$m). We probe the disk gap down to a radius of confidence of 0.1'' (${\sim}15$ AU at 145 pc). We interpret the results in the context of dust evolution when planets interact with the parental disk. Results. We observe a gap from 0.1 to 0.3'' (${\sim}15$ to 40AU) and a bright annulus as previously detected by HiCIAO $H$ band observations at $1.65\mu$m. The radial width of the annulus is around $40$AU, and its centre is at ${\sim}61$AU from the central star. The peak of the reflected light at 0.626$\mu$m is located 20 AU inward of the cavity detected in the submillimetre. In addition, we detect a dip at a position angle of ${\sim}46.2 \pm 5.4^\circ$. A dip was also detected with HiCIAO, but located at ${\sim}85^\circ$. If the dip observed with HiCIAO is the same, this suggests an average dip rotation of ${\sim}12^\circ/$year, which is inconsistent with the local Keplerian angular velocity of $\sim$0.8$^\circ$/yr at $\sim$61AU. Conclusions. The spatial discrepancy in the radial emission in J160421.7-213028 at different wavelengths is consistent with dust filtration at the outer edge of a gap carved by a massive planet. The dip rotation can be interpreted by fast variability of the inner disk and/or the presence of a warp or circumplanetary material of a planet at ${\sim}9.6$ AU.

Understanding the formation of planetary nebulae remains elusive because in the preceding asymtotic giant branch (AGB) phase these stars are heavily enshrouded in an optically thick dusty envelope. To further understand the morphology of the circumstellar environments of AGB stars we observe the closest carbon-rich AGB star IRC+10216 in scattered light. When imaged in scattered light at optical wavelengths, IRC+10216 surprisingly shows a narrow equatorial density enhancement, in contrast to the large-scale spherical rings that have been imaged much further out. We use radiative transfer models to interpret this structure in terms of two models: firstly, an equatorial density enhancement, commonly observed in the more evolved post-AGB stars, and secondly, in terms of a dust rings model, where a local enhancement of mass-loss creates a spiral ring as the star rotates. We conclude that both models can be used to reproduce the dark lane in the scattered light images, which is caused by an equatorially density enhancement formed by dense dust rather than a bipolar outflow as previously thought. We are unable to place constraints on the formation of the equatorial density enhancement by a binary system.

Transition disks have dust-depleted inner regions and may represent an intermediate step of an on-going disk dispersal process, where planet formation is probably in progress. Recent millimetre observations of transition disks reveal radially and azimuthally asymmetric structures, where micron- and millimetre-sized dust particles may not spatially coexist. These properties can be the result of particle trapping and grain growth in pressure bumps originating from the disk interaction with a planetary companion. The multiple features observed in some transition disks, such as SR 21, suggest the presence of more than one planet. We study the gas and dust distributions of a disk hosting two massive planets as function of different disk and dust parameters. Observational signatures, such as the spectral energy distribution, sub-millimetre, and polarised images are simulated for the various parameters. We confirm that planets can lead to particle trapping, although for a disk with high viscosity ($\alpha_{\rm{turb}}=10^{-2}$), the planet should be more massive than $5 M_{\rm{Jup}}$ and dust fragmentation should occur with low efficiency ($v_{f}\sim30\rm{m s}^{-1}$). This will lead to a ring-like feature as observed in transition disks in the millimetre. When trapping occurs, we find that a smooth distribution of micron sized grains throughout the disk, sometimes observed in scattered light, can only happen if the combination of planet mass and turbulence is such that small grains are not fully filtered out. A high disk viscosity ($\alpha_{\rm{turb}}=10^{-2}$) ensures a replenishment of the cavity in micron-sized dust, while for lower viscosity ($\alpha_{\rm{turb}}=10^{-3}$), the planet mass is constrained to be less than $5 M_{\rm{Jup}}$. In these cases, the gas distribution is likely to show low-amplitude azimuthal asymmetries caused by disk eccentricity rather than by long-lived vortices.

Most stars form in a clustered environment. Therefore, it is important to assess how this environment influences the evolution of protoplanetary discs around young stars. In turn, this affects their ability to produce planets and ultimately life. We present here for the first time 3D SPH/N-body simulations that include both the hydrodynamical evolution of the discs around their natal stars, as well as the dynamics of the stars themselves. The discs are viscously evolving, accreting mass onto the central star and spreading. We find penetrating encounters to be very destructive for the discs as in previous studies, although the frequency of such encounters is low. We also find, however, that encounter influence the disc radii more strongly than other disc properties such as the disc mass. The disc sizes are set by the competition between viscous spreading and the disruptive effect of encounters. As discs spread, encounters become more and more important. In the regime of rapid spreading encounters simply truncate the discs, stripping the outer portions. In the opposite regime, we find that the effect of many distant encounters is able to limit the disc size. Finally, we predict from our simulations that disc sizes are limited by encounters at stellar densities exceeding $\sim 2-3 \times 10^3 \ \mathrm{pc}^{-2}$.

We present four new epochs of Ks-band images of the young pre-transitional disk around LkCa 15, and perform extensive forward modeling to derive the physical parameters of the disk. We find indications of strongly anisotropic scattering (Henyey-Greenstein g = 0.67 [-0.11,+0.18]) and a significantly tapered gap edge ('round wall'), but see no evidence that the inner disk, whose existence is predicted by the spectral energy distribution, shadows the outer regions of the disk visible in our images. We marginally confirm the existence of an offset between the disk center and the star along the line of nodes; however, the magnitude of this offset (x = 27 [-20,+19] mas) is notably lower than that found in our earlier H-band images (Thalmann et al. 2010). Intriguingly, we also find, at high significance, an offset of y = 69 [-25, +49] mas perpendicular to the line of nodes. If confirmed by future observations, this would imply a highly elliptical -- or otherwise asymmetric -- disk gap with an effective eccentricity of e = ~0.3. Such asymmetry would most likely be the result of dynamical sculpting by one or more unseen planets in the system. Finally, we find that the bright arc of scattered light we see in direct imaging observations originates from the near side of the disk, and appears brighter than the far side because of strong forward scattering.

The ~39-m European Extremely Large Telescope (E-ELT) will be the largest telescope ever built. This makes it particularly suitable for sensitive polarimetric observations, as polarimetry is a photon-starved technique. However, the telescope mirrors may severely limit the polarimetric accuracy of instruments on the Nasmyth platforms by creating instrumental polarisation and/or modifying the polarisation signal of the object. In this paper we characterise the polarisation effects of the two currently considered designs for the E-ELT Nasmyth ports as well as the effect of ageing of the mirrors. By means of the Mueller matrix formalism, we compute the response matrices of each mirror arrangement for a range of zenith angles and wavelengths. We then present two techniques to correct for these effects that require the addition of a modulating device at the polarisation-free intermediate focus that acts either as a switch or as a part of a two-stage modulator. We find that the values of instrumental polarisation, Stokes transmission reduction and cross- talk vary significantly with wavelength, and with pointing, for the lateral Nasmyth case, often exceeding the accuracy requirements for proposed polarimetric instruments. Realistic ageing effects of the mirrors after perfect calibration of these effects may cause polarimetric errors beyond the requirements. We show that the modulation approach with a polarimetric element located in the intermediate focus reduces the instrumental polarisation effects down to tolerable values, or even removes them altogether. The E-ELT will be suitable for sensitive and accurate polarimetry, provided frequent calibrations are carried out, or a dedicated polarimetric element is installed at the intermediate focus.

Transitional discs are a special type of protoplanetary discs where planet formation is thought to be taking place. These objects feature characteristic inner cavities and/or gaps of a few tens of AUs in the sub-millimitre images of the disc. This signature suggests a localised depletion of matter in the disc that could be caused by planet formation processes. However, recent observations have revealed differences in the structures imaged at different wavelengths in some of these discs. In this paper, we aim to explain these observational differences using self-consistent physical 2-D hydrodynamical and dust evolution models of such objects, assuming their morphology is indeed generated by the presence of a planet. We use these models to derive the distribution of gas and dust in a theoretical planet-hosting disc, for various planet masses and orbital separations. We then simulate observations of the emitted and scattered light from these models with VLT/SPHERE ZIMPOL, Subaru/HiCIAO, VLT/VISIR and ALMA. We do this by first computing the full resolution images of the models at different wavelengths, and then simulating the observations accounting for the characteristics of each particular instrument. The presence of the planet generates pressure bumps in the gas distribution of the disc whose characteristics strongly depend on the planet mass and position. These bumps cause large grains to accumulate while small grains are allowed into inner regions. This spatial differentiation of the grain sizes explains the differences in the observations since different wavelengths and observing techniques trace different parts of the dust size distribution. Based on this effect, we conclude that the combination of visible/near-infrared polarimetric and sub-mm images is the best strategy to constrain the properties of the unseen planet responsible for the disc structure.

Well over 700 exoplanets have been detected to date. Only a handful of these have been observed directly. Direct observation is extremely challenging due to the small separation and very large contrast involved. Imaging polarimetry offers a way to decrease the contrast between the unpolarized starlight and the light that has become linearly polarized after scattering by circumstellar material. This material can be the dust and debris found in circumstellar disks, but also the atmosphere or surface of an exoplanet. We present the design, calibration approach, polarimetric performance and sample observation results of the Extreme Polarimeter, an imaging polarimeter for the study of circumstellar environments in scattered light at visible wavelengths. The polarimeter uses the beam-exchange technique, in which the two orthogonal polarization states are imaged simultaneously and a polarization modulator swaps the polarization states of the two beams before the next image is taken. The instrument currently operates without the aid of Adaptive Optics. To reduce the effects of atmospheric seeing on the polarimetry, the images are taken at a frame rate of 35 fps, and large numbers of frames are combined to obtain the polarization images. Four successful observing runs have been performed using this instrument at the 4.2 m William Herschel Telescope on La Palma, targeting young stars with protoplanetary disks as well as evolved stars surrounded by dusty envelopes. In terms of fractional polarization, the instrument sensitivity is better than 10^-4. The contrast achieved between the central star and the circumstellar source is of the order 10^-6. We show that our calibration approach yields absolute polarization errors below 1%.

We present observational evidence of environmental effects on the formation and evolution of planetary systems. We combine catalogues of resolved protoplanetary discs (PPDs) and young stellar objects in the solar neighbourhood to analyse the PPD size distribution as a function of ambient stellar density. By running Kolmogorov-Smirnov tests between the PPD radii at different densities, we find empirical evidence, at the >97% confidence level, for a change in the PPD radius distribution at ambient stellar densities \Sigma > 10^{3.5} pc^{-2}. This coincides with a simple theoretical estimate for the truncation of PPDs or planetary systems by dynamical encounters. If this agreement is causal, the ongoing disruption of PPDs and planetary systems limits the possible existence of planets in the habitable zone, with shorter lifetimes at higher host stellar masses and ambient densities. Therefore, habitable planets are not likely to be present in long-lived stellar clusters, and may have been ejected altogether to form a population of unbound, free-floating planets. We conclude that, while highly suggestive, our results should be verified through other methods. Our simple model shows that truncations should lead to a measurable depletion of the PPD mass function that can be detected with ALMA observations of the densest nearby and young clusters.

Context. The Kepler object KIC 12557548 b is peculiar. It exhibits transit-like features every 15.7 hours that vary in depth between 0.2% and 1.2%. Rappaport et al. (2012) explain the observations in terms of a disintegrating, rocky planet that has a trailing cloud of dust created and constantly replenished by thermal surface erosion. The variability of the transit depth is then a consequence of changes in the cloud optical depth. Aims. We aim to validate the disintegrating-planet scenario by modeling the detailed shape of the observed light curve, and thereby constrain the cloud particle properties to better understand the nature of this intriguing object. Methods. We analysed the six publicly-available quarters of raw Kepler data, phase-folded the light curve and fitted it to a model for the trailing dust cloud. Constraints on the particle properties were investigated with a light-scattering code. Results. The light curve exhibits clear signatures of light scattering and absorption by dust, including a brightening in flux just before ingress correlated with the transit depth and explained by forward scattering, and an asymmetry in the transit light curve shape, which is easily reproduced by an exponentially decaying distribution of optically thin dust, with a typical grain size of 0.1 micron. Conclusions. Our quantitative analysis supports the hypothesis that the transit signal of KIC 12557548 b is due to a variable cloud of dust, most likely originating from a disintegrating object.

Although different approaches to model a polarimeter's accuracy have been described before, a complete error budgeting tool for polarimetric systems has not been yet developed. Based on the framework introduced by Keller & Snik, in 2009, we have developed the M&m's code as a first attempt to obtain a generic tool to model the performance and accuracy of a given polarimeter, including all the potential error contributions and their dependencies on physical parameters. The main goal of the code is to provide insight on the combined influence of many polarization errors on the accuracy of any polarimetric instrument. In this work we present the mathematics and physics based on which the code is developed as well as its general structure and operational scheme. Discussion of the advantages of the M&m's approach to error budgeting and polarimetric performance simulation is carried out and a brief outlook of further development of the code is also given.

We present recent polarimetric images of the highly variable star R CrB using ExPo and archival WFPC2 images from the HST. We observed R CrB during its current dramatic minimum where it decreased more than 9 mag due to the formation of an obscuring dust cloud. Since the dust cloud is only in the line-of-sight, it mimics a coronograph allowing the imaging of the star's circumstellar environment. Our polarimetric observations surprisingly show another scattering dust cloud at approximately 1.3" or 2000 AU from the star. We find that to obtain a decrease in the stellar light of 9 mag and with 30% of the light being reemitted at infrared wavelengths (from R CrB's SED) the grains in R CrB's circumstellar environment must have a very low albedo of approximately 0.07%. We show that the properties of the dust clouds formed around R CrB are best fitted using a combination of two distinct populations of grains size. The first are the extremely small 5 nm grains, formed in the low density continuous wind, and the second population of large grains (~0.14 {\mu}m) which are found in the ejected dust clouds. The observed scattering cloud, not only contains such large grains, but is exceptionally massive compared to the average cloud.