Recent high-contrast imaging surveys, looking for planets in young, nearby systems showed evidence of a small number of giant planets at relatively large separation beyond typically 20 au where those surveys are the most sensitive. Access to smaller physical separations between 5 and 20 au is the next step for future planet imagers on 10 m telescopes and ELTs in order to bridge the gap with indirect techniques (radial velocity, transit, astrometry with Gaia). In that context, we recently proposed a new algorithm, Keplerian-Stacker, combining multiple observations acquired at different epochs and taking into account the orbital motion of a potential planet present in the images to boost the ultimate detection limit. We showed that this algorithm is able to find planets in time series of simulated images of SPHERE even when a planet remains undetected at one epoch. Here, we validate the K-Stacker algorithm performances on real SPHERE datasets, to demonstrate its resilience to instrumental speckles and the gain offered in terms of true detection. This will motivate future dedicated multi-epoch observation campaigns in high-contrast imaging to search for planets in emitted and reflected light. Results. We show that K-Stacker achieves high success rate when the SNR of the planet in the stacked image reaches 7. The improvement of the SNR ratio goes as the square root of the total exposure time. During the blind test and the redetection of HD 95086 b, and betaPic b, we highlight the ability of K-Stacker to find orbital solutions consistent with the ones derived by the state of the art MCMC orbital fitting techniques, confirming that in addition to the detection gain, K-Stacker offers the opportunity to characterize the most probable orbital solutions of the exoplanets recovered at low signal to noise.

SPHERE (Beuzit et al,. 2019) has now been in operation at the VLT for more than 5 years, demonstrating a high level of performance. SPHERE has produced outstanding results using a variety of operating modes, primarily in the field of direct imaging of exoplanetary systems, focusing on exoplanets as point sources and circumstellar disks as extended objects. The achievements obtained thus far with SPHERE (~200 refereed publications) in different areas (exoplanets, disks, solar system, stellar physics...) have motivated a large consortium to propose an even more ambitious set of science cases, and its corresponding technical implementation in the form of an upgrade. The SPHERE+ project capitalizes on the expertise and lessons learned from SPHERE to push high contrast imaging performance to its limits on the VLT 8m-telescope. The scientific program of SPHERE+ described in this document will open a new and compelling scientific window for the upcoming decade in strong synergy with ground-based facilities (VLT/I, ELT, ALMA, and SKA) and space missions (Gaia, JWST, PLATO and WFIRST). While SPHERE has sampled the outer parts of planetary systems beyond a few tens of AU, SPHERE+ will dig into the inner regions around stars to reveal and characterize by mean of spectroscopy the giant planet population down to the snow line. Building on SPHERE's scientific heritage and resounding success, SPHERE+ will be a dedicated survey instrument which will strengthen the leadership of ESO and the European community in the very competitive field of direct imaging of exoplanetary systems. With enhanced capabilities, it will enable an even broader diversity of science cases including the study of the solar system, the birth and death of stars and the exploration of the inner regions of active galactic nuclei.

The occurrence rate of long-period giant planets around young stars is highly uncertain since it is not only governed by the protoplanetary disc structure and planet formation process, but also reflects dynamical re-structuring processes after planet formation as well as possible capture of planets not formed in-situ. Direct imaging is currently the only feasible method to detect such wide-orbit planets and constrain their occurrence rate. We carry out a large L'-band high-contrast direct imaging survey for giant planets around young stars with protoplanetary or debris discs using the NACO instrument at the ESO Very Large Telescope on Cerro Paranal in Chile. We use very deep angular differential imaging observations with typically >60 deg field rotation, and employ a vector vortex coronagraph where feasible to achieve the best possible point source sensitivity down to an inner working angle of about 100mas. This paper introduces our NACO Imaging Survey for Planets around Young stars ("NACO-ISPY"), its goals and strategy, the target list, and data reduction scheme, and presents preliminary results from the first 2.5 survey years. We achieve a mean 5 sigma L' contrast of 6.4mag at 150mas and a background limit of 16.5mag at >1.5". Our detection probability is >50\% for companions with 8\,M$_{\rm Jup}$\ at semi-major axes 80-200au. It thus compares well to the detection space of other state-of-the-art high-contrast imaging surveys. We have contributed to the characterisation of two new planets originally discovered by VLT/SPHERE, but we have not yet independently discovered new planets around any of our target stars. We report the discovery of close-in low-mass stellar companions around four young stars and show L'-band scattered light images of the discs around eleven stars, six of which have never been imaged at L'-band before.

Context. PDS 70 is a young (5.4 Myr), nearby (~113 pc) star hosting a known transition disk with a large gap. Recent observations with SPHERE and NACO in the near-infrared (NIR) allowed us to detect a planetary mass companion, PDS70b, within the disk cavity. Moreover, observations in H_alpha with MagAO and MUSE revealed emission associated to PDS70b and to another new companion candidate, PDS70c, at a larger separation from the star. Aims. Our aim is to confirm the discovery of the second planet PDS70c using SPHERE at VLT, to further characterize its physical properties, and search for additional point sources in this young planetary system. Methods. We re-analyzed archival SPHERE NIR observations and obtained new data in Y, J, H and K spectral bands for a total of four different epochs. The data were reduced using the data reduction and handling pipeline and the SPHERE data center. We then applied custom routines (e.g. ANDROMEDA and PACO) to subtract the starlight. Results. We re-detect both PDS 70 b and c and confirm that PDS70c is gravitationally bound to the star. We estimate this second planet to be less massive than 5 M Jup and with a T_eff around 900 K. Also, it has a low gravity with log g between 3.0 and 3.5 dex. In addition, a third object has been identified at short separation (~0.12") from the star and gravitationally bound to the star. Its spectrum is however very blue, so that we are probably seeing stellar light reflected by dust and our analysis seems to demonstrate that it is a feature of the inner disk. We, however, cannot completely exclude the possibility that it is a planetary mass object enshrouded by a dust envelope. In this latter case, its mass should be of the order of few tens of M_Earth. Moreover, we propose a possible structure for the planetary system based on our data that, however, cannot be stable on a long timescale.

Context. The interaction between low-mass companions and the debris discs they reside in is still not fully understood. A debris disc can evolve due to self-stirring, a process in which planetesimals can excite their neighbours to the point of destructive collisions. In addition, the presence of a companion could further stir the disc (companion-stirring). Additional information is necessary to understand this fundamental step in the formation and evolution of a planetary system, and at the moment of writing only a handful of systems are known where a companion and a debris disc have both been detected and studied at the same time. Aims. Our primary goal is to augment the sample of these systems and to understand the relative importance between self-stirring and companion-stirring. Methods. In the course of the VLT/NaCo Imaging Survey for Planets around Young stars (ISPY), we observed HD 193571, an A0 debris disc hosting star at a distance of 68 pc with an age between 60 and 170 Myr. We obtained two sets of observations in L' band and a third epoch in H band using the GPI instrument at Gemini-South. Results. A companion was detected in all three epochs at a projected separation of 11 au (0.17 arcsec), and co-motion was confirmed through proper motion analysis. Given the inferred disc size of 120 au, the companion appears to reside within the gap between the host star and the disc. Comparison between the L' and H band magnitude and evolutionary tracks suggests a mass of 0.31 - 0.39 solar masses . Conclusions. We discovered a previously unknown M-dwarf companion around HD 193571, making it the third low-mass stellar object discovered within a debris disc. A comparison to self- and companion-stirring models suggests that the companion is likely responsible for the stirring of the disc.

PDS 70b is the most robustly detected young planet imaged in the gap of a transition disk so far, found at a separation of ~195 mas (~22 au) from its host star and at a position angle of about 155 deg. This system is therefore a unique laboratory to characterize the properties of young planetary systems at the stage of their formation. We aim to trace direct and indirect imprints of PDS 70b on the gas and dust emission of the circumstellar disk in order to study the properties of this ~5 Myr young planetary system. We obtained ALMA band 7 observations of PDS 70 in dust continuum and $^{12}CO$ (3-2) and combined them with archival data resulting in an unprecedented angular resolution of about 70 mas (~8 au). We derive an upper limit on circumplanetary material at the location of PDS 70b of ~0.01 $M_{\oplus}$ and find a highly structured circumstellar disk in both dust and gas. The outer dust ring peaks at ~0.65" (74 au) and reveals a possible second unresolved peak at ~0.53" (60 au). The integrated intensity of CO also shows evidence of a depletion of emission at ~0.2" (23 au) with a width of ~0.1" (11 au). The gas kinematics show evidence of a deviation from Keplerian rotation inside $\lesssim$0.8" (91 au). This implies a pressure gradient that can account for the location of the dust ring well beyond the location of PDS 70b. Further in, we detect an inner disk which appears to be connected to the outer disk by a possible bridge feature in the North-West region in both gas and dust. We compare the observations to hydrodynamical simulations including a planet with different masses covering the estimated mass range previously derived from NIR photometry (~5-9 M$_{Jup}$). We find that even a planet with a mass of 10 M$_{Jup}$ may not be sufficient to explain the extent of the wide gap and an additional low-mass companion may be needed to account for the observed disk morphology.

Context. The 51 Eridani system harbors a complex architecture with its primary star forming a hierarchical system with the binary GJ 3305AB at a projected separation of 2000 au, a giant planet orbiting the primary star at 13 au, and a low-mass debris disk around the primary star with possibly a cold component and a warm component inferred from the spectral energy distribution. Aims. We aim to better constrain the orbital parameters of the known giant planet. Methods. We monitored the system over three years from 2015 to 2018 with the VLT/SPHERE exoplanet imaging instrument. Results. We measure an orbital motion for the planet of ~130 mas with a slightly decreasing separation (~10 mas) and find a hint of curvature. This potential curvature is further supported at 3$\sigma$ significance when including literature GPI astrometry corrected for calibration systematics. Fits of the SPHERE and GPI data using three complementary approaches provide broadly similar results. The data suggest an orbital period of 32$^{+17}_{-9}$ yr (i.e. 12$^{+4}_{-2}$ au in semi-major axis), an inclination of 133$^{+14}_{-7}$ deg, an eccentricity of 0.45$^{+0.10}_{-0.15}$, and an argument of periastron passage of 87$^{+34}_{-30}$ deg [mod 180 deg]. The time at periastron passage and the longitude of node exhibit bimodal distributions because we do not detect yet if the planet is accelerating or decelerating along its orbit. Given the inclinations of the planet's orbit and of the stellar rotation axis (134-144 deg), we infer alignment or misalignment within 18 deg for the star-planet spin-orbit. Further astrometric monitoring in the next 3-4 years is required to confirm at a higher significance the curvature in the planet's motion, determine if the planet is accelerating or decelerating on its orbit, and further constrain its orbital parameters and the star-planet spin-orbit.

Aims. R Coronae Australis (R CrA) is the brightest star of the Coronet nebula of the Corona Australis (CrA) star forming region. It has very red colors, probably due to dust absorption and it is strongly variable. High contrast instruments allow for an unprecedented direct exploration of the immediate circumstellar environment of this star. Methods. We observed R CrA with the near-IR channels (IFS and IRDIS) of SPHERE at VLT. In this paper, we used four different epochs, three of them from open time observations while one is from the SPHERE guaranteed time. The data were reduced using the DRH pipeline and the SPHERE Data Center. On the reduced data we implemented custom IDL routines with the aim to subtract the speckle halo.We have also obtained pupil-tracking H-band (1.45-1.85 micron) observations with the VLT/SINFONI near-infrared medium-resolution (R~3000) spectrograph. Results. A companion was found at a separation of 0.156" from the star in the first epoch and increasing to 0.18400 in the final one. Furthermore, several extended structures were found around the star, the most noteworthy of which is a very bright jet-like structure North-East from the star. The astrometric measurements of the companion in the four epochs confirm that it is gravitationally bound to the star. The SPHERE photometry and the SINFONI spectrum, once corrected for extinction, point toward an early M spectral type object with a mass between 0.3 and 0.55 M?. The astrometric analyis provides constraints on the orbit paramenters: e~0.4, semi-major axis at 27-28 au, inclination of ~ 70{\deg} and a period larger than 30 years. We were also able to put constraints of few MJup on the mass of possible other companions down to separations of few tens of au.

HIP 65426 b is a recently discovered exoplanet imaged during the course of the SPHERE-SHINE survey. Here we present new $L'$ and $M'$ observations of the planet from the NACO instrument at the VLT from the NACO-ISPY survey, as well as a new $Y-H$ spectrum and $K$-band photometry from SPHERE-SHINE. Using these data, we confirm the nature of the companion as a warm, dusty planet with a mid-L spectral type. From comparison of its SED with the BT-Settl atmospheric models, we derive a best-fit effective temperature of $T_{\text{eff}}=1618\pm7$ K, surface gravity $\log g=3.78^{+0.04}_{-0.03}$ and radius $R=1.17\pm0.04$ $R_{\text{J}}$ (statistical uncertainties only). Using the DUSTY and COND isochrones we estimate a mass of $8\pm1$ $M_{\text{J}}$. Combining the astrometric measurements from our new datasets and from the literature, we show the first indications of orbital motion of the companion (2.6$\sigma$ significance) and derive preliminary orbital constraints. We find a highly inclined orbit ($i=107^{+13}_{-10}$ deg) with an orbital period of $800^{+1200}_{-400}$ yr. We also report SPHERE sparse aperture masking observations that investigate the possibility that HIP 65426 b was scattered onto its current orbit by an additional companion at a smaller orbital separation. From this data we rule out the presence of brown dwarf companions with masses greater than 16 $M_{\text{J}}$ at separations larger than 3 AU, significantly narrowing the parameter space for such a companion.

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.

The nearby Herbig Be star HD100546 is known to be a laboratory for the study of protoplanets and their relation with the circumstellar disk that is carved by at least 2 gaps. We observed the HD100546 environment with high contrast imaging exploiting several different observing modes of SPHERE, including datasets with/without coronagraphs, dual band imaging, integral field spectroscopy and polarimetry. The picture emerging from these different data sets is complex. Flux-conservative algorithms images clearly show the disk up to 200au. More aggressive algorithms reveal several rings and warped arms overlapping the main disk. The bright parts of this ring lie at considerable height over the disk mid-plane at about 30au. Our images demonstrate that the brightest wings close to the star in the near side of the disk are a unique structure, corresponding to the outer edge of the intermediate disk at ~40au. Modeling of the scattered light from the disk with a geometrical algorithm reveals that a moderately thin structure can well reproduce the light distribution in the flux-conservative images. We suggest that the gap between 44 and 113 au span between the 1:2 and 3:2 resonance orbits of a massive body located at ~70au that might coincide with the candidate planet HD100546b detected with previous thermal IR observations. In this picture, the two wings can be the near side of a ring formed by disk material brought out of the disk at the 1:2 resonance with the same massive object. While we find no clear evidence confirming detection of the planet candidate HD100546c in our data, we find a diffuse emission close to the expected position of HD100546b. This source can be described as an extremely reddened substellar object surrounded by a dust cloud or its circumplanetary disk. Its astrometry is broadly consistent with a circular orbital motion on the disk plane.

Context. A low-mass brown dwarf has been recently imaged around HR 2562 (HD 50571), a star hosting a debris disk resolved in the far infrared. Interestingly, the companion location is compatible with an orbit coplanar with the disk and interior to the debris belt. This feature makes the system a valuable laboratory to analyze the formation of substellar companions in a circumstellar disk and potential disk-companion dynamical interactions. Aims. We aim to further characterize the orbital motion of HR 2562 B and its interactions with the host star debris disk. Methods. We performed a monitoring of the system over ~10 months in 2016 and 2017 with the VLT/SPHERE exoplanet imager. Results. We confirm that the companion is comoving with the star and detect for the first time an orbital motion at high significance, with a current orbital motion projected in the plane of the sky of 25 mas (~0.85 au) per year. No orbital curvature is seen in the measurements. An orbital fit of the SPHERE and literature astrometry of the companion without priors on the orbital plane clearly indicates that its orbit is (quasi-)coplanar with the disk. To further constrain the other orbital parameters, we used empirical laws for a companion chaotic zone validated by N-body simulations to test the orbital solutions that are compatible with the estimated disk cavity size. Non-zero eccentricities (>0.15) are allowed for orbital periods shorter than 100 yr, while only moderate eccentricities up to ~0.3 for orbital periods longer than 200 yr are compatible with the disk observations. A comparison of synthetic Herschel images to the real data does not allow us to constrain the upper eccentricity of the companion.

Young circumstellar disks are of prime interest to understand the physical and chemical conditions under which planet formation takes place. Only very few detections of planet candidates within these disks exist, and most of them are currently suspected to be disk features. In this context, the transition disk around the young star PDS 70 is of particular interest, due to its large gap identified in previous observations, indicative of ongoing planet formation. We aim to search for the presence of planets and search for disk structures indicative for disk-planet interactions and other evolutionary processes. We analyse new and archival near-infrared (NIR) images of the transition disk PDS 70 obtained with the VLT/SPHERE, VLT/NaCo and Gemini/NICI instruments in polarimetric differential imaging (PDI) and angular differential imaging (ADI) modes. We detect a point source within the gap of the disk at about 195 mas (about 22 au) projected separation. The detection is confirmed at five different epochs, in three filter bands and using different instruments. The astrometry results in an object of bound nature, with high significance. The comparison of the measured magnitudes and colours to evolutionary tracks suggests that the detection is a companion of planetary mass. We confirm the detection of a large gap of about 54 au in size within the disk in our scattered light images, and detect a signal from an inner disk component. We find that its spatial extent is very likely smaller than about 17 au in radius. The images of the outer disk show evidence of a complex azimuthal brightness distribution which may in part be explained by Rayleigh scattering from very small grains. Future observations of this system at different wavelengths and continuing astrometry will allow us to test theoretical predictions regarding planet-disk interactions, planetary atmospheres and evolutionary models.

The G-type star GJ504A is known to host a 3 to 35 MJup companion whose temperature, mass, and projected separation all contribute to make it a test case for the planet formation theories and for atmospheric models of giant planets and light brown dwarfs. We collected data from the CHARA interferometer, SOPHIE spectrograph, and VLT/SPHERE high contrast imager to revisit the properties of the system. We measure a radius of 1.35+/- 0.04Rsun for GJ504A which yields isochronal ages of 21+/-2Myr or 4.0+/-1.8Gyr for the system and line-of-sight stellar rotation axis inclination of $162.4_{-4.3}^{+3.8}$ degrees or $18.6_{-3.8}^{+4.3}$ degrees. We re-detect the companion in the Y2, Y3, J3, H2, and K1 dual band SPHERE images. The complete 1-4 $\mu$m SED shape of GJ504b is best reproduced by T8-T9.5 objects with intermediate ages ($\leq1.5$Gyr), and/or unusual dusty atmospheres and/or super-solar metallicities. All six atmospheric models used yield $\mathrm{T_{eff}=550 \pm 50}$K for GJ504b and point toward a low surface gravity (3.5-4.0 dex). The accuracy on the metallicity value is limited by model-to-model systematics. It is not degenerate with the C/O ratio. We derive $\mathrm{log\:L/L_{\odot}=-6.15\pm0.15}$ dex for the companion compatible with masses of $\mathrm{M=1.3^{+0.6}_{-0.3}M_{Jup}}$ and $\mathrm{M=23^{+10}_{-9} M_{Jup}}$ for the young and old age ranges, respectively. The semi-major axis (sma) is above 27.8 au and the eccentricity lower than 0.55. The posterior on GJ~504b's orbital inclination suggests a misalignment with GJ~504A rotation axis. We combine the radial velocity and multi-epoch imaging data to exclude additional objects (90\% prob.) more massive than 2.5 and 30 $\mathrm{M_{Jup}}$ with sma in the range 0.01-80 au for the young and old system ages, respectively. The companion is in the envelope of the population of planets synthetized with our core-accretion model.

Aims: We aim to characterize the orbital and atmospheric properties of PDS 70 b, which was first identified on May 2015 in the course of the SHINE survey with SPHERE, the extreme adaptive-optics instrument at the VLT. Methods: We obtained new deep SPHERE/IRDIS imaging and SPHERE/IFS spectroscopic observations of PDS 70 b. The astrometric baseline now covers 6 years which allows us to perform an orbital analysis. For the first time, we present spectrophotometry of the young planet which covers almost the entire near-infrared range (0.96 to 3.8 micrometer). We use different atmospheric models covering a large parameter space in temperature, log(g), chemical composition, and cloud properties to characterize the properties of the atmosphere of PDS 70 b. Results: PDS 70 b is most likely orbiting the star on a circular and disk coplanar orbit at ~22 au inside the gap of the disk. We find a range of models that can describe the spectrophotometric data reasonably well in the temperature range between 1000-1600 K and log(g) no larger than 3.5 dex. The planet radius covers a relatively large range between 1.4 and 3.7 R_jupiter with the larger radii being higher than expected from planet evolution models for the age of the planet of 5.4 Myr. Conclusions: This study provides a comprehensive dataset on the orbital motion of PDS 70 b, indicating a circular orbit and a motion coplanar with the disk. The first detailed spectral energy distribution of PDS 70 b indicates a temperature typical for young giant planets. The detailed atmospheric analysis indicates that a circumplanetary disk may contribute to the total planet flux.

Evolutionary models are widely used to infer the mass of stars, brown dwarfs, and giant planets. Their predictions are thought to be less reliable at young ages ($<$ 200 Myr) and in the low-mass regime ($\mathrm{<1~M_{\odot}}$). GJ 2060 AB and TWA 22 AB are two rare astrometric M-dwarf binaries respectively members of the AB Doradus and Beta Pictoris moving groups. As their dynamical mass can be measured within a few years, they can be used to calibrate the evolutionary tracks and set new constraints on the age of young moving groups. We find a total mass of $\mathrm{0.18\pm 0.02~M_\odot}$ for TWA 22. That mass is in good agreement with model predictions at the age of the Beta Pic moving group. We obtain a total mass of $\mathrm{1.09 \pm 0.10~M_{\odot}}$ for GJ 2060. We estimate a spectral type of M$1\pm0.5$, $\mathrm{L/L_{\odot}=-1.20\pm0.05}$ dex, and $\mathrm{T_{eff}=3700\pm100}$ K for GJ 2060 A. The B component is a M$3\pm0.5$ dwarf with $\mathrm{L/L_{\odot}=-1.63\pm0.05}$ dex and $\mathrm{T_{eff}=3400\pm100}$ K. The dynamical mass of GJ 2060 AB is inconsistent with the most recent models predictions (BCAH15, PARSEC) for an ABDor age in the range 50-150 Myr. It is 10 to 20\% (1-2 sigma, depending on the assumed age) above the models predictions, corresponding to an underestimation of $0.10$ to $0.20~\mathrm{M_\odot}$. Coevality suggests a young age for the system ($\sim$ 50 Myr) according to most evolutionary models. TWA 22 validates the predictions of recent evolutionary tracks at $\sim$20 Myr. On the other hand, we evidence a 1-2 sigma mismatch between the predicted and observed mass of GJ 2060 AB. This slight departure may indicate that one of the star hosts a tight companion. Alternatively, this would confirm the models tendency to underestimate the mass of young low-mass stars.

The study of high contrast imaged brown dwarfs and exoplanets depends strongly on evolutionary models. To estimate the mass of a directly imaged substellar object, its extracted photometry or spectrum is used and adjusted with model spectra together with the estimated age of the system. These models still need to be properly tested and constrained. HD 4747B is a brown dwarf close to the H burning mass limit, orbiting a nearby, solar-type star and has been observed with the radial velocity method over almost two decades now. Its companion was also recently detected by direct imaging, allowing a complete study of this particular object. We aim to fully characterize HD 4747B by combining a well constrained dynamical mass and a study of its observed spectral features in order to test evolutionary models for substellar objects and characterize its atmosphere. We combine the radial velocity measurements of HIRES and CORALIE taken over two decades and high contrast imaging of several epochs from NACO, NIRC2 and SPHERE to obtain a dynamical mass. From the SPHERE data we obtain a low resolution spectrum of the companion from Y to H band, as well as two narrow band-width photometric measurements in the K band. A study of the primary star allows in addition to constrain the age of the system as well as its distance. Thanks to the new SPHERE epoch and NACO archival data combined with previous imaging data and high precision radial velocity measurements, we have been able to derive a well constrained orbit. We derive a dynamical mass of mB=70.0$\pm$1.6 MJup which is higher than a previous study, but in better agreement with the models. By comparing the object with known brown dwarfs spectra, we derive a spectral type of L9 and an effective temperature of 1350$\pm$50 K. With a retrieval analysis we constrain the oxygen and carbon abundances and compare them with the ones from the HR 8799 planets.

The nearby and young M star AU Mic is surrounded by a debris disk in which we previously identified a series of large-scale arch-like structures that have never been seen before in any other debris disk and that move outward at high velocities. We initiated a monitoring program with the following objectives: 1) track the location of the structures and better constrain their projected speeds, 2) search for new features emerging closer in, and ultimately 3) understand the mechanism responsible for the motion and production of the disk features. AU Mic was observed at 11 different epochs between August 2014 and October 2017 with the IR camera and spectrograph of SPHERE. These high-contrast imaging data were processed with a variety of angular, spectral, and polarimetric differential imaging techniques to reveal the faintest structures in the disk. We measured the projected separations of the features in a systematic way for all epochs. We also applied the very same measurements to older observations from the Hubble Space Telescope (HST) with the visible cameras STIS and ACS. The main outcomes of this work are 1) the recovery of the five southeastern broad arch-like structures we identified in our first study, and confirmation of their fast motion (projected speed in the range 4-12 km/s); 2) the confirmation that the very first structures observed in 2004 with ACS are indeed connected to those observed later with STIS and now SPHERE; 3) the discovery of two new very compact structures at the northwest side of the disk (at 0.40" and 0.55" in May 2015) that move to the southeast at low speed; and 4) the identification of a new arch-like structure that might be emerging at the southeast side at about 0.4" from the star (as of May 2016). Abridged.

HD95086 (A8V, 17Myr) hosts a rare planetary system for which a multi-belt debris disk and a giant planet of 4-5MJup have been directly imaged. Our study aims to characterize the physical and orbital properties of HD95086b, search for additional planets at short and wide orbits and image the cold outer debris belt in scattered light. We used HARPS at the ESO 3.6m telescope to monitor the radial velocity of HD95086 over 2 years and investigate the existence of giant planets at less than 3au orbital distance. With the IRDIS dual-band imager and the IFS integral field spectrograph of SPHERE at VLT, we imaged the faint circumstellar environment beyond 10au at six epochs between 2015 and 2017. We do not detect additional giant planets around HD95086. We identified the nature (bound companion or background contaminant) of all point-like sources detected in the IRDIS field of view. None of them correspond to the ones recently discovered near the edge of the cold outer belt by ALMA. HD95086b is resolved for the first time in J-band with IFS. Its near-infrared spectral energy distribution is well fitted by a few dusty and/or young L7-L9 dwarf spectral templates. The extremely red 1-4um spectral distribution is typical of low-gravity objects at the L/T spectral type transition. The planet's orbital motion is resolved between January 2015 and May 2017. Together with past NaCo measurements properly re-calibrated, our orbital fitting solutions favor a retrograde low to moderate-eccentricity orbit e=0.2 (0.0 to 0.5), with a semi-major axis 52au corresponding to orbital periods of 288$ yrs and an inclination that peaks at i = 141deg, which is compatible with a planet-disk coplanar configuration. Finally, we report the detection in polarimetric differential imaging of the cold outer debris belt between 100 and 300au, consistent in radial extent with recent ALMA 1.3mm resolved observations.

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.

We present observations of the Herbig Ae star HD169142 with VLT/SPHERE instruments InfraRed Dual-band Imager and Spectrograph (IRDIS) ($K1K2$ and $H2H3$ bands) and the Integral Field Spectrograph (IFS) ($Y$, $J$ and $H$ bands). We detect several bright blobs at $\sim$180 mas separation from the star, and a faint arc-like structure in the IFS data. Our reference differential imaging (RDI) data analysis also finds a bright ring at the same separation. We show, using a simulation based on polarized light data, that these blobs are actually part of the ring at 180 mas. These results demonstrate that the earlier detections of blobs in the $H$ and $K_S$ bands at these separations in Biller et al. as potential planet/substellar companions are actually tracing a bright ring with a Keplerian motion. Moreover, we detect in the images an additional bright structure at $\sim$93 mas separation and position angle of 355$^{\circ}$, at a location very close to previous detections. It appears point-like in the $YJ$ and $K$ bands but is more extended in the $H$ band. We also marginally detect an inner ring in the RDI data at $\sim$100 mas. Follow-up observations are necessary to confirm the detection and the nature of this source and structure.

A large number of systems harboring a debris disk show evidence for a double belt architecture. One hypothesis for explaining the gap between the belts is the presence of one or more planets dynamically carving it. This work aims to investigate this scenario in systems harboring two components debris disks. All the targets in the sample were observed with the SPHERE instrument which performs high-contrast direct imaging. Positions of the inner and outer belts were estimated by SED fitting of the infrared excesses or, when available, from resolved images of the disk. Very few planets have been observed so far in debris disks gaps and we intended to test if such non-detections depend on the observational limits of the present instruments. This aim is achieved by deriving theoretical predictions of masses, eccentricities and semi-major axes of planets able to open the observed gaps and comparing such parameters with detection limits obtained with SPHERE. The relation between the gap and the planet is due to the chaotic zone around the orbit of the planet. The radial extent of this zone depends on the mass ratio between the planet and the star, on the semi-major axis and on the eccentricity of the planet and it can be estimated analytically. We apply the formalism to the case of one planet on a circular or eccentric orbit. We then consider multi-planetary systems: 2 and 3 equal-mass planets on circular orbits and 2 equal-mass planets on eccentric orbits in a packed configuration. We then compare each couple of values (M,a), derived from the dynamical analysis of single and multiple planetary models, with the detection limits obtained with SPHERE. Our results show that the apparent lack of planets in gaps between double belts could be explained by the presence of a system of two or more planets possibly of low mass and on an eccentric orbits whose sizes are below the present detection limits.

In this paper we present the results of the SPHERE observation of the HD 284149 system, aimed at a more detailed characterisation of both the primary and its brown dwarf companion. We observed HD 284149 in the near-infrared with SPHERE, using the imaging mode (IRDIS+IFS) and the long-slit spectroscopy mode (IRDIS-LSS). The data were reduced using the dedicated SPHERE pipeline, and algorithms such as PCA and TLOCI were applied to reduce the speckle pattern. The IFS images revealed a previously unknown low-mass (~0.16$M_{\odot}$) stellar companion (HD 294149 B) at ~0.1$^{\prime\prime}$, compatible with previously observed radial velocity differences, as well as proper motion differences between Gaia and Tycho-2 measurements. The known brown dwarf companion (HD 284149 b) is clearly visible in the IRDIS images. This allowed us to refine both its photometry and astrometry. The analysis of the medium resolution IRDIS long slit spectra also allowed a refinement of temperature and spectral type estimates. A full reassessment of the age and distance of the system was also performed, leading to more precise values of both mass and semi-major axis. As a result of this study, HD 284149 ABb therefore becomes the latest addition to the (short) list of brown dwarfs on wide circumbinary orbits, providing new evidence to support recent claims that object in such configuration occur with a similar frequency to wide companions to single stars.

Context. The SAO 206462 (HD 135344B) disk is one of the few known transitional disks showing asymmetric features in scattered light and thermal emission. Near-infrared scattered-light images revealed two bright outer spiral arms and an inner cavity depleted in dust. Giant protoplanets have been proposed to account for the disk morphology. Aims. We aim to search for giant planets responsible for the disk features and, in the case of non-detection, to constrain recent planet predictions using the data detection limits. Methods. We obtained new high-contrast and high-resolution total intensity images of the target spanning the Y to the K bands (0.95-2.3 mic) using the VLT/SPHERE near-infrared camera and integral field spectrometer. Results. The spiral arms and the outer cavity edge are revealed at high resolutions and sensitivities without the need for image post-processing techniques, which introduce photometric biases. We do not detect any close-in companions. For the derivation of the detection limits on putative giant planets embedded in the disk, we show that the knowledge of the disk aspect ratio and viscosity is critical for the estimation of the attenuation of a planet signal by the protoplanetary dust because of the gaps that these putative planets may open. Given assumptions on these parameters, the mass limits can vary from ~2-5 to ~4-7 Jupiter masses at separations beyond the disk spiral arms. The SPHERE detection limits are more stringent than those derived from archival NaCo/L' data and provide new constraints on a few recent predictions of massive planets (4-15 MJ) based on the spiral density wave theory. The SPHERE and ALMA data do not favor the hypotheses on massive giant planets in the outer disk (beyond 0.6). There could still be low-mass planets in the outer disk and/or planets inside the cavity.

51 Eridani b is an exoplanet around a young (20 Myr) nearby (29.4 pc) F0-type star, recently discovered by direct imaging. Being only 0.5" away from its host star it is well suited for spectroscopic analysis using integral field spectrographs. We aim to refine the atmospheric properties of this and to further constrain the architecture of the system by searching for additional companions. Using the SPHERE instrument at the VLT we extend the spectral coverage of the planet to the complete Y- to H-band range and provide photometry in the K12-bands (2.11, 2.25 micron). The object is compared to other cool and peculiar dwarfs. Furthermore, the posterior probability distributions of cloudy and clear atmospheric models are explored using MCMC. We verified our methods by determining atmospheric parameters for the two benchmark brown dwarfs Gl 570D and HD 3651B. For probing the innermost region for additional companions, archival VLT-NACO (L') SAM data is used. We present the first spectrophotometric measurements in the Y- and K-bands for the planet and revise its J-band flux to values 40% fainter than previous measurements. Cloudy models with uniform cloud coverage provide a good match to the data. We derive the temperature, radius, surface gravity, metallicity and cloud sedimentation parameter f_sed. We find that the atmosphere is highly super-solar (Fe/H~1.0) with an extended, thick cloud cover of small particles. The model radius and surface gravity suggest planetary masses of about 9 M_jup. The evolutionary model only provides a lower mass limit of >2 M_jup (for pure hot-start). The cold-start model cannot explain the planet's luminosity. The SPHERE and NACO/SAM detection limits probe the 51 Eri system at Solar System scales and exclude brown-dwarf companions more massive than 20 M_jup beyond separations of ~2.5 au and giant planets more massive than 2 M_jup beyond 9 au.

We find a new substellar companion to the Pleiades member star, Pleiades HII 3441, using the Subaru telescope with adaptive optics. The discovery is made as part of the high-contrast imaging survey to search for planetary-mass and substellar companions in the Pleiades and young moving groups. The companion has a projected separation of 0".49 +/- 0".02 (66 +/- 2 AU) and a mass of 68 +/- 5 M_J based on three observations in the J-, H-, and K_S-band. The spectral type is estimated to be M7 (~2700 K), and thus no methane absorption is detected in the H band. Our Pleiades observations result in the detection of two substellar companions including one previously reported among 20 observed Pleiades stars, and indicate that the fraction of substellar companions in the Pleiades is about 10.0 +26.1/-8.8 %. This is consistent with multiplicity studies of both the Pleiades stars and other open clusters.

Debris disks offer valuable insights into the latest stages of circumstellar disk evolution, and can possibly help us to trace the outcomes of planetary formation processes. In the age range 10 to 100\,Myr, most of the gas is expected to have been removed from the system, giant planets (if any) must have already been formed, and the formation of terrestrial planets may be on-going. Pluto-sized planetesimals, and their debris released in a collisional cascade, are under their mutual gravitational influence, which may result into non-axisymmetric structures in the debris disk. High angular resolution observations are required to investigate these effects and constrain the dynamical evolution of debris disks. Furthermore, multi-wavelength observations can provide information about the dust dynamics by probing different grain sizes. Here we present new VLT/SPHERE and ALMA observations of the debris disk around the 40\,Myr-old solar-type star HD\,61005. We resolve the disk at unprecedented resolution both in the near-infrared (in scattered and polarized light) and at millimeter wavelengths. Thanks to the new observations, we propose a solution for both the radial and azimuthal distribution of the dust grains in the debris disk. We find that the disk has a moderate eccentricity ($e \sim 0.1$) and that the dust density is two times larger at the pericenter compared to the apocenter. With no giant planets detected in our observations, we investigate alternative explanations besides planet-disk interactions to interpret the inferred disk morphology. We postulate that the morphology of the disk could be the consequence of a massive collision between $\sim$\,1000\,km-sized bodies at $\sim$\,61\,au. If this interpretation holds, it would put stringent constraints on the formation of massive planetesimals at large distances from the star.

HD~106906AB is so far the only young binary system around which a planet has been imaged and a debris disk evidenced thanks to a strong IR excess. As such, it represents a unique opportunity to study the dynamics of young planetary systems. We aim at further investigating the close (tens of au scales) environment of the HD~106906AB system. We used the extreme AO fed, high contrast imager SPHERE recently installed on the VLT to observe HD~106906. Both the IRDIS imager and the Integral Field Spectrometer were used. We discovered a very inclined, ring-like disk at a distance of 65~au from the star. The disk shows a strong brightness asymmetry with respect to its semi-major axis. It shows a smooth outer edge, compatible with ejection of small grains by the stellar radiation pressure. We show furthermore that the planet's projected position is significantly above the disk's PA. Given the determined disk inclination, it is not excluded though that the planet could still orbit within the disk plane if at a large separation (2000--3000 au). We identified several additional point sources in the SPHERE/IRDIS field-of-view, that appear to be background objects. We compare this system with other debris disks sharing similarities, and we briefly discuss the present results in the framework of dynamical evolution.

We compared the number of faint stars detected in deep survey fields with the current stellar distribution model of the Galaxy and found that the detected number in the H band is significantly smaller than the predicted number. This indicates that M-dwarfs, the major component, are fewer in the halo and the thick disk. We used archived data of several surveys in both the north and south field of GOODS (Great Observatories Origins Deep Survey), MODS in GOODS-N, and ERS and CANDELS in GOODS-S. The number density of M-dwarfs in the halo has to be 20+/-13% relative to that in the solar vicinity, in order for the detected number of stars fainter than 20.5 mag in the H band to match with the predicted value from the model. In the thick disk, the number density of M-dwarfs must be reduced (52+/-13%) or the scale height must be decreased (~600 pc). Alternatively, overall fractions of the halo and thick disks can be significantly reduced to achieve the same effect, because our sample mainly consists of faint M-dwarfs. Our results imply that the M-dwarf population in regions distant from the Galactic plane is significantly smaller than previously thought. We then discussed the implications this has on the suitability of the model predictions for the prediction of non-companion faint stars in direct imaging extrasolar planet surveys by using the best-fit number densities.

In this work, we want to find out if the IMF can be determined from colour images, integrated colours, or mass-to-light ratios, especially at high redshift, where galaxies cannot be resolved into individual stars, which would enable us to investigate dependencies of the IMF on cosmological epoch. We use chemo-dynamical models to investigate the influence of the Initial Mass Function (IMF) on the evolution of a Milky Way-type disk galaxy, in particular of its colours. We find that the effect of the IMF on the internal gas absorption is larger than its effect on the light from the stellar content. However, the two effects work in the opposite sense: An IMF with more high mass stars leads to brighter and bluer star-light, but also to more interstellar dust and thus to more absorption, causing a kind of "IMF degeneracy". The most likely wavelength region in which to detect IMF effects is the infrared (i. e., JHK). We also provide photometric absorption and inclination corrections in the SDSS ugriz and the HST WFPC2 and NICMOS systems.

We investigate the origin of the clumpy structures observed at high redshift, like the chain galaxies. We use a three dimensional chemodynamical simulation describing the dynamics of stars and a two-phase interstellar medium, as well as feedback processes from the stars. For high efficiency of energy dissipation in the cold cloud medium, the initially gaseous disk fragments and develops several massive clumps of gas and stars. We follow the evolution of the individual clumps and determine their masses, metallicities and velocities. A few dynamical times after fragmentation of the disk, the clumps merge to build a massive bulge. Calculating HST- and UBVRIJHKLM-colors, including absorption by interstellar dust, we determine the morphologies and colors of this model in HST images. Several peculiar morphological structures seen in the HDF can be well-explained by a fragmented galactic disk model, including chain galaxies and objects consisting of several nearby knots.

We investigate the evolution of star-forming gas-rich disks, using a 3D chemodynamical model including a dark halo, stars, and a two-phase interstellar medium with feedback processes from the stars. We show that galaxy evolution proceeds along very different routes depending on whether it is the gas disk or the stellar disk which first becomes unstable, as measured by the respective Q-parameters. This in turn depends on the uncertain efficiency of energy dissipation of the cold cloud component from which stars form. When the cold gas cools efficiently and drives the instability, the galactic disk fragments and forms a number of massive clumps of stars and gas. The clumps spiral to the center of the galaxy in a few dynamical times and merge there to form a central bulge component in a strong starburst. When the kinetic energy of the cold clouds is dissipated at a lower rate, stars form from the gas in a more quiescent mode, and an instability only sets in at later times, when the surface density of the stellar disk has grown sufficiently high. The system then forms a stellar bar, which channels gas into the center, evolves, and forms a bulge whose stars are the result of a more extended star formation history. We investigate the stability of the gas-stellar disks in both regimes, as well as the star formation rates and element enrichment. We study the morphology of the evolving disks, calculating spatially resolved colours from the distribution of stars in age and metallicity, including dust absorption. We then discuss morphological observations such as clumpy structures and chain galaxies at high redshift as possible signatures of fragmenting, gas-rich disks. Finally, we investigate abundance ratio distributions as a means to distinguish the different scenarios for bulge formation.

The astrophysical nature of r-process sites is a long standing mystery and many probable sources have been suggested in the past, among them lower-mass core-collapse supernovae (in the range 8 - 10 Msol), higher-mass core-collapse supernovae (with masses > 20 Msol) and neutron star mergers. In this work, we present a detailed inhomogeneous chemical evolution study that considers for the first time neutron star mergers as major r-process sources, and compare this scenario to the ones in which core-collapse supernovae act as dominant r-process sites. We conclude that, due to the lack of reliable iron and r-process yields as function of progenitor mass, it is not possible at present to distinguish between the lower-mass and higher-mass supernovae scenarios within the framework of inhomogeneous chemical evolution. However, neutron-star mergers seem to be ruled out as the dominant r-process source, since their low rates of occurrence would lead to r-process enrichment that is not consistent with observations at very low metallicities. Additionally, the considerable injection of r-process material by a single neutron-star merger leads to a scatter in [r-process/Fe] ratios at later times which is much too large compared to observations.

We present a dynamical model for the formation and evolution of a massive disk galaxy, within a growing dark halo whose mass evolves according to cosmological simulations of structure formation. The galactic evolution is simulated with a new 3D chemo-dynamical code, including dark matter, stars and a multi-phase ISM. The simulations start at redshift z=4.85 with a small dark halo in a LCDM universe and we follow the evolution until the present epoch. The energy release by massive stars and SNe prevents a rapid collapse of the baryonic matter and delays the maximum star formation until z=1. The galaxy forms radially from inside-out and vertically from halo to disk. The first galactic component that forms is the halo, followed by the bulge, the disk-halo transition region, and the disk. At z=1, a bar begins to form which later turns into a triaxial bulge. There is a pronounced deficiency of low-metallicity disk stars due to pre-enrichment of the disk ISM with metal-rich gas from the bulge and inner disk (G-dwarf problem). The mean rotation and the distribution of orbital eccentricities for all stars as a function of metallicity are not very different from those observed in the solar neighbourhood, showing that homogeneous collapse models are oversimplified. The approach presented here provides a detailed description of the formation and evolution of an isolated disk galaxy in a LCDM universe, yielding new information about the kinematical and chemical history of the stars and the ISM, but also about the evolution of the luminosity, the colours and the morphology of disk galaxies.

Inhomogeneous chemical evolution models of galaxies which try to reproduce the scatter seen in element-to-iron ratios of metal-poor halo stars are heavily dependent on theoretical nucleosynthesis yields of core-collapse supernovae. Hence inhomogeneous models present themselves as a test for stellar nucleosyn- thesis calculations. Applying an inhomogeneous chemical evolution model to our Galaxy reveals a number of shortcomings of existing theoretical nucleosynthesis yields. One problem is the predicted scatter in [O/Fe] and [Mg/Fe] which is too large compared to the one observed in metal-poor halo stars. This can be either due to the O or Mg yields or due to the Fe yields (or both). However, O and Mg are alpha-elements that are produced mainly during hydrostatic burning and thus are not affected by the theoretical uncertainties afflicting the collapse and explosion of a massive star. Stellar iron yields, on the other hand, depend heavily on the choice of the mass-cut between ejecta and proto neutron star and are therefore very uncertain. We present Fe yield distributions as function of progenitor mass that are consistent with the abundance distribution of metal- poor halo stars and are in agreement with observed Ni yields of SNe II with known progenitor masses. The iron yields of lower-mass SNe II (in the range 10-20 Msol) are well constrained by those observations. Present observations, however, do not allow to determine a unique solution for higher-mass SNe. Nevertheless, the main dependence of the stellar iron yield as function of progenitor mass may be derived and can be used as constraint for future supernova/hypernova models. The results are of importance for the earliest stages of galaxy formation when the ISM is dominated by chemical inhomogeneities and the instantaneous mixing approximation is not valid.

We calculate synthetic UBVRIJHKLM images, integrated spectra and colours for the disk galaxy formation models of Samland & Gerhard (2002), from redshift z=4 to z=0. Two models are considered, an accretion model based on LambdaCDM structure formation simulations, and a classical collapse model in a dark matter halo. Both models provide the star formation history and dynamics of the baryonic component within a three-dimensional chemo-dynamical description. To convert to spectra and colours, we use the latest, metallicity-calibrated spectral library of Westera et al. (2002), including internal absorption. As a first application, we compare the derived colours with Hubble Deep Field North bulge colours and find good agreement. With our model, we disentangle metallicity effects and absorption effects on the integrated colours, and find that absorption effects are dominant for redshift z < 1.5. Furthermore, we confirm the quality of m_K as a mass tracer, and find indications for a correlation between (J-K) and metallicity gradients.

We introduce a stochastic halo formation model to compute the early chemical enrichment of the interstellar medium (ISM) of the halo. Contrary to 1-zone chemical evolution models, we are able to resolve local inhomogeneities in the ISM caused by single core-collapse supernovae. These inhomogeneities lead to different element abundance patterns in very metal-poor stars, which can be seen as scatter in the abundances of halo stars with metallicities [Fe/H] < 2.0. The early chemical evolution of the halo proceeds in different enrichment phases: At [Fe/H] < -3.0, the halo ISM is unmixed and dominated by local inhomogeneities caused by single core-collapse supernova (SN) events. For metallicities [Fe/H] >-2.0 the halo ISM is well mixed, showing an element abundance pattern integrated over the initial mass function. In the range -3.0 < [Fe/H] < -2.0 a continuous transition from the unmixed to the well mixed ISM occurs. For some elements (Si, Ca, Eu), the scatter in the element-to-iron ratio [El/Fe] seen in the observations of metal-poor halo stars can be reproduced. Stellar yields of other elements predict a scatter which, compared to the observations, is too large (O, Mg) or too small (Ni). Cr and Mn show a decreasing trend for lower metallicities, which can not be explained by metallicity independent yields, provided that the mixing of the ejecta with the interstellar medium does not depend on progenitor mass. This demonstrates the need for revised, self-consistent core-collapse SN yields.