sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Reggiani, Maddalena M.

Normalized to: Reggiani, M.

31 article(s) in total. 252 co-authors, from 1 to 12 common article(s). Median position in authors list is 11,0.

[1]  [pdf] - 2118321
Is HR 6819 a triple system containing a black hole? -- An alternative explanation
Comments: Submitted to A&A, 13 pages (16 figures and 2 tables); 4 pages supplementary material (4 figures and 4 tables). Comments are welcome
Submitted: 2020-06-18
HR 6819 was recently proposed to be a triple system consisting of an inner B-type giant + black hole binary with an orbital period of 40d and an outer Be tertiary. This interpretation is mainly based on two inferences: that the emission attributed to the outer Be star is stationary, and that the inner star, which is used as mass calibrator for the black hole, is a B-type giant. We re-investigate the properties of HR 6819 by spectral disentangling and an atmosphere analysis of the disentangled spectra to search for a possibly simpler alternative explanation for HR 6819. Disentangling implies that the Be component is not a static tertiary, but rather a component of the binary in the 40-d orbit. The inferred radial velocity amplitudes imply an extreme mass ratio of M_2/M_1 = 15 +/- 3. We infer spectroscopic masses of 0.4$^{+0.3}_{-0.1}$ Msun and 6$^{+5}_{-3}$ Msun for the primary and secondary, which agree well with the dynamical masses for an inclination of i = 32 deg. This indicates that the primary might be a stripped star rather than a B-type giant. Evolutionary modelling suggests that a possible progenitor system would be a tight (P_i ~ 2d) B+B binary system that experienced conservative mass transfer. While the observed nitrogen enrichment of the primary conforms with the predictions of the evolutionary models, we find no indications for the predicted He enrichment. We suggest that HR 6819 is a binary system consisting of a stripped B-type primary and a rapidly-rotating Be star that formed from a previous mass-transfer event. In the framework of this interpretation, HR 6819 does not contain a black hole. Interferometry can distinguish between these two scenarios by providing an independent measurement of the separation between the visible components.
[2]  [pdf] - 2107385
The Carina High-Contrast Imaging Project for massive Stars (CHIPS) I. Methodology and proof of concept on QZ Car ($\equiv$ HD93206)
Comments: Submitted and accepted by A&A. 17 pages, 10 figures
Submitted: 2020-06-02
Massive stars like company. However, low-mass companions have remained extremely difficult to detect at angular separations ($\rho$) smaller than 1" (~1000-3000 au considering typical distance to nearby massive stars) given the large brightness contrast between the companion and the central star. Constraints on the low-mass end of the companions mass-function for massive stars are however needed, for example to help distinguishing between various scenarios for the formation of massive stars. To obtain statistically significant constraint on the presence of low-mass companions beyond the typical detection limit of current surveys ($\Delta \mathrm{mag} \lesssim 5$ at $\rho \lesssim 1$"), we initiated a survey of O and Wolf-Rayet stars in the Carina region using the SPHERE coronagraphic instrument on the VLT. In this first paper, we aim to introduce the survey, to present the methodology and to demonstrate the capability of SPHERE for massive stars using the multiple system QZ~Car. High-contrast imaging techniques, such as angular- and spectral-differential imaging techniques as well as PSF-fitting, were applied to detect and measure the relative flux of companions in each spectral channel of the instrument. We detected 19 sources around the QZ~Car system with detection limits of 9~mag at $\rho > 200$~mas for IFS and as faint as 13~mag at $\rho$ > 1".8 for IRDIS (corresponding to sub-solar masses for potential companions). All but two are reported here for the first time. Based on this proof of concept, we showed that VLT/SPHERE allows us to reach the sub-solar mass regime of the companion mass function. This paves the way for this type of observation with a large sample of massive stars to provide novel constraints on the multiplicity of massive stars in a region of the parameter space that has remained inaccessible so far.
[3]  [pdf] - 2084386
The hidden companion in LB-1 unveiled by spectral disentangling
Comments: Submitted to A&A. 5 pages (4 figures, 1 table); 2 pages of supporting material (3 figures and 1 table)
Submitted: 2020-04-27
The intriguing binary LS V + 22 25 (LB-1) drew much attention after being claimed to be a 79-day binary comprising a B-type star orbiting a ~70Ms black hole - the most massive stellar black hole reported to date. Subsequent studies showed that evidence for such a large companion mass is lacking. Recent analyses implied that the primary star is a stripped, He-rich star with peculiar sub-solar abundances of heavy elements such as Mg and Fe. However, the nature of the secondary, which was proposed to be a black hole, a neutron star, or a main sequence star, remains unknown. Relying on 26 newly acquired spectroscopic observations secured with the HERMES and FEROS spectrographs, we perform an orbital analysis and spectral disentangling of LB-1 to elucidate the nature of the system. We find that LB-1 contains two non-degenerate stars. The "hidden" secondary is a rapidly rotating (vsini ~ 300 km/s) B3 V star with a decretion disk - a Be star. The properties of the primary match those predicted for stripped stars: it is He- and N-rich and portrays significant Balmer-line emission, presumably stemming from its wind. Adopting a light contribution in the optical of 55% for the stripped primary, the abundances of heavy elements are found to be solar. With the derived value of K2 = 11.2 +- 1.0 km/s and by calibrating the mass of the B3 Ve secondary to M2 = 7 +- 2Ms, we derive an orbital mass for the stripped primary of M1 = 1.5 +- 0.4Ms. The orbital inclination of 39 +- 4deg implies a near-critical rotation for the Be secondary Veq ~ 470 km/s. Hence, LB-1 does not contain a compact object. Instead, it is a rare Be binary system consisting of a stripped donor star and a Be mass accretor rotating at near its critical velocity. This system is a clear example that binary interactions play a decisive role in the production of rapid stellar rotators and Be stars.
[4]  [pdf] - 2057764
Near-Infrared Imaging of a Spiral in the CQ Tau Disk
Comments: accepted for publication in AJ
Submitted: 2019-10-16, last modified: 2020-01-24
We present $L^\prime$-band Keck/NIRC2 imaging and $H$-band Subaru/AO188+HiCIAO polarimetric observations of CQ Tau disk with a new spiral arm. Apart from the spiral feature our observations could not detect any companion candidates. We traced the spiral feature from the $r^2$-scaled HiCIAO polarimetric intensity image and the fitted result is used for forward modeling to reproduce the ADI-reduced NIRC2 image. We estimated the original surface brightness after throughput correction in $L^\prime$-band to be $\sim126$ mJy/arcsec$^2$ at most. We suggest that the grain temperature of the spiral may be heated up to $\sim$200 K in order to explain both of the $H$- and $L^{\prime}$-bands results. The $H$-band emission at the location of the spiral originates from the scattering from the disk surface while both scattering and thermal emission may contribute to the $L^{\prime}$-band emission. If the central star is only the light source of scattered light, the spiral emission at $L^\prime$-band should be thermal emission. If an inner disk also acts as the light source, the scattered light and the thermal emission may equally contribute to the $L^\prime$-band spiral structure.
[5]  [pdf] - 2014974
Discovery of a directly imaged disk in scattered light around the Sco-Cen member Wray 15-788
Comments: Accepted for publication in A&A (16 pages, 11 figures)
Submitted: 2019-02-20, last modified: 2019-12-16
As part of our on-going survey we have carried out high-contrast imaging with VLT/SPHERE/IRDIS to obtain polarized and total intensity images of the young ($11^{+16}_{-7}$Myr old) K3IV star Wray 15-788 within the Lower Centaurus Crux subgroup of Sco-Cen. For the total intensity images, we remove the stellar halo by an approach based on reference star differential imaging in combination with principal component analysis. Both total intensity and polarimetric data resolve a disk around Wray 15-788. Modeling of the stellar spectral energy distribution suggests that this is a protoplanetary disk at a transition stage. We detect a bright, outer ring at a projected separation of $\sim$370mas ($\approx$56au), hints for inner substructures at $\sim$170mas ($\approx$28au) and a gap in between. Only within a position angle range of $60^\circ<\varphi<240^\circ$, we are confident at 5$\sigma$ level to detect actual scattered light flux from the outer ring of the disk; the remaining part is indistinguishable from background noise. For the detected part of the outer ring we determine a disk inclination of $i$=21$^\circ\pm$6$^\circ$ and a position angle of $\varphi$=76$^\circ\pm$16$^\circ$. Furthermore, we find that Wray 15-788 is part of a binary system with the A2V star HD 98363 at a separation of $\sim$50'' ($\approx$6900au). The detection of only half of the outer ring might be due to shadowing by a misaligned inner disk. A potential substellar companion can cause the misalignment of the inner structures and can be responsible for clearing the detected gap from scattering material. We can not, however, rule out the possibility of a non-detection due to our limited signal to noise ratio, combined with brightness azimuthal asymmetry. From our data, we can exclude companions more massive than 10$M_\mathrm{jup}$ within the gap at a separation of $\sim$230mas ($\approx$35au).
[6]  [pdf] - 2086447
No signature of the orbital motion of a putative 70 solar mass black hole in LB-1
Comments: 10 pages, 8 figures
Submitted: 2019-12-09
Liu et al. (2019) recently reported the detection of a 68 [+11/-13] solar mass (Msun) black hole (BH) paired with an 8.2 [+0.9/-1.2] Msun B-type sub-giant star in the 78.9-day spectroscopic binary system LB-1. Such a black hole is over twice as massive as any other known stellar-mass black hole with non-compact companions2 and its mass approaches those that result from BH-BH coalescences that are detected by gravitational wave interferometers. Its presence in a solar-like metallicity environment challenges conventional theories of massive binary evolution, stellar winds and core-collapse supernovae, so that more exotic scenarios seem to be needed to explain the existence and properties of LB-1. Here, we show that the observational diagnostics used to derive the BH mass results from the orbital motion of the B-type star, not that of the BH. As a consequence, no evidence for a massive BH remains in the data, therefore solving the existing tension with formation models of such a massive BH at solar metallicity and with theories of massive star evolution in general.
[7]  [pdf] - 2030688
The Young Suns Exoplanet Survey: Detection of a wide orbit planetary mass companion to a solar-type Sco-Cen member
Comments: Accepted for publication in MNRAS (15 pages, 9 figures)
Submitted: 2019-12-09
The Young Suns Exoplanet Survey (YSES) consists of a homogeneous sample of 70 young, solar-mass stars located in the Lower Centaurus-Crux subgroup of the Scorpius-Centaurus association with an average age of $15\pm3\,$Myr. We report the detection of a co-moving companion around the K3IV star TYC 8998-760-1 (2MASSJ13251211-6456207) that is located at a distance of $94.6\pm0.3\,$pc using SPHERE/IRDIS on the VLT. Spectroscopic observations with VLT/X-SHOOTER constrain the mass of the star to $1.00\pm0.02\,M_{\odot}$ and an age of $16.7\pm1.4\,$Myr. The companion TYC 8998-760-1 b is detected at a projected separation of 1.71'', which implies a projected physical separation of $162\,$au. Photometric measurements ranging from $Y$ to $M$ band provide a mass estimate of $14\pm3\,M_\mathrm{jup}$ by comparison to BT-Settl and AMES-dusty isochrones, corresponding to a mass ratio of $q=0.013\pm0.003$ with respect to the primary. We rule out additional companions to TYC 8998-760-1 that are more massive than $12\,M_\mathrm{jup}$ and farther than $12\,$au away from the host. Future polarimetric and spectroscopic observations of this system with ground and space based observatories will facilitate testing of formation and evolution scenarios shaping the architecture of the circumstellar environment around this 'young Sun'.
[8]  [pdf] - 1863978
Exploring the RCrA environment with SPHERE: Discovery of a new stellar companion
Comments: 11 pages, 13 figures, accepted by A&A
Submitted: 2019-02-07
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.
[9]  [pdf] - 1833262
Reference star differential imaging of close-in companions and circumstellar disks with the NIRC2 vortex coronagraph at W.M. Keck Observatory
Comments: Accepted for publication in The Astronomical Journal
Submitted: 2019-01-13, last modified: 2019-01-14
Reference star differential imaging (RDI) is a powerful strategy for high contrast imaging. Using example observations taken with the vortex coronagraph mode of Keck/NIRC2 in $L^\prime$ band, we demonstrate that RDI provides improved sensitivity to point sources at small angular separations compared to angular differential imaging (ADI). Applying RDI to images of the low-mass stellar companions HIP 79124 C (192 mas separation, $\Delta L^\prime$=4.01) and HIP 78233 B (141 mas separation, $\Delta L^\prime$=4.78), the latter a first imaging detection, increases the significance of their detections by up to a factor of 5 with respect to ADI. We compare methods for reference frames selection and find that pre-selection of frames improves detection significance of point sources by up to a factor of 3. In addition, we use observations of the circumstellar disks around MWC 758 and 2MASS J16042165-2130284 to show that RDI allows for accurate mapping of scattered light distributions without self-subtraction artifacts.
[10]  [pdf] - 1822965
High-Resolution ALMA Observations of HD100546: Asymmetric Circumstellar Ring, and Circumplanetary Disk Upper Limits
Comments: 14 pages, accepted by ApJ
Submitted: 2018-11-26
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.
[11]  [pdf] - 1811143
Deep exploration of $\epsilon$ Eridani with Keck Ms-band vortex coronagraphy and radial velocities: mass and orbital parameters of the giant exoplanet
Comments: Accepted to AJ
Submitted: 2018-10-08, last modified: 2018-10-29
We present the most sensitive direct imaging and radial velocity (RV) exploration of $\epsilon$ Eridani to date. $\epsilon$ Eridani is an adolescent planetary system, reminiscent of the early Solar system. It is surrounded by a prominent and complex debris disk which is likely stirred by one or several gas giant exoplanets. The discovery of the RV signature of a giant exoplanet was announced 15 years ago, but has met with scrutiny due to possible confusion with stellar noise. We confirm the planet with a new compilation and analysis of precise RV data spanning 30 years, and combine it with upper limits from our direct imaging search, the most sensitive ever performed. The deep images were taken in the Ms band (4.7$\mu$m) with the vortex coronagraph recently installed in W.M. Keck Observatory's infrared camera NIRC2, which opens a sensitive window for planet searches around nearby adolescent systems. The RV data and direct imaging upper limit maps were combined in an innovative joint Bayesian analysis, providing new constraints on the mass and orbital parameters of the elusive planet. $\epsilon$ Eridani b has a mass of $0.78^{+0.38}_{-0.12}$ $M_{Jup}$ and is orbiting $\epsilon$ Eridani at about $3.48\pm 0.02$ AU with a period of $7.37 \pm 0.07$ years. The eccentricity of $\epsilon$ Eridani b's orbit is $0.07^{+0.06}_{-0.05}$, an order of magnitude smaller than early estimates and consistent with a circular orbit. We discuss our findings from the standpoint of planet-disk interactions and prospects for future detection and characterization with the James Webb Space Telescope.
[12]  [pdf] - 1719138
Hi-5: a potential high-contrast thermal near-infrared imager for the VLTI
Comments: 15 pages, 3 figures, SPIE proceedings, based on arXiv:1801.04148
Submitted: 2018-07-19
Hi-5 is a high-contrast (or high dynamic range) infrared imager project for the VLTI. Its main goal is to characterize young extra-solar planetary systems and exozodiacal dust around southern main-sequence stars. In this paper, we present an update of the project and key technology pathways to improve the contrast achieved by the VLTI. In particular, we discuss the possibility to use integrated optics, proven in the near-infrared, in the thermal near-infrared (L and M bands, 3-5~$\mu$m) and advanced fringe tracking strategies. We also address the strong exoplanet science case (young exoplanets, planet formation, and exozodiacal disks) offered by this wavelength regime as well as other possible science cases such as stellar physics (fundamental parameters and multiplicity) and extragalactic astrophysics (active galactic nuclei and fundamental constants). Synergies and scientific preparation for other potential future instruments such as the Planet Formation Imager are also briefly discussed.
[13]  [pdf] - 1703380
High-contrast imaging of HD 163296 with the Keck/NIRC2 L'-band vortex coronograph
Comments: 9 pages, 6 figures. Accepted by MNRAS
Submitted: 2018-06-21
We present observations of the nearby (D$\sim$100\,pc) Herbig star HD~163296 taken with the vortex coronograph at Keck/NIRC2 in the L' band (3.7~$\mu$m), to search for planetary mass companions in the ringed disc surrounding this pre-main sequence star. The images reveal an arc-like region of scattered light from the disc surface layers that is likely associated with the first bright ring detected with ALMA in the $\lambda$=1.3mm dust continuum at $\sim$65~au. We also detect a point-like source at $\sim$0\farcs5 projected separation in the North-East direction, close to the inner edge of the second gap in the millimetre images. Comparing the point source photometry with the atmospheric emission models of non-accreting giant planets, we obtain a mass of 6--7~M$_J$ for a putative protoplanet, assuming a system age of 5~Myr. Based on the contrast at a 95\% level of completeness calculated on the emission-free regions of our images, we set upper limits for the masses of giant planets of 8--15~M$_J$, 4.5--6.5~M$_J$ and 2.5-4.0~M$_J$ at the locations of the first, second and third gap in the millimetre dust continuum, respectively. Further deep, high resolution thermal IR imaging of the HD~163296 system are warranted, to confirm the presence and nature of the point source and to better understand the structure of the dust disc.
[14]  [pdf] - 1709313
The path towards high-contrast imaging with the VLTI: the Hi-5 project
Comments: 24 pages, 2 figures, accepted for publication in Experimental Astronomy
Submitted: 2018-01-12, last modified: 2018-05-25
The development of high-contrast capabilities has long been recognized as one of the top priorities for the VLTI. As of today, the VLTI routinely achieves contrasts of a few 10$^{-3}$ in the near-infrared with PIONIER (H band) and GRAVITY (K band). Nulling interferometers in the northern hemisphere and non-redundant aperture masking experiments have, however, demonstrated that contrasts of at least a few 10$^{-4}$ are within reach using specific beam combination and data acquisition techniques. In this paper, we explore the possibility to reach similar or higher contrasts on the VLTI. After reviewing the state-of-the-art in high-contrast infrared interferometry, we discuss key features that made the success of other high-contrast interferometric instruments (e.g., integrated optics, nulling, closure phase, and statistical data reduction) and address possible avenues to improve the contrast of the VLTI by at least one order of magnitude. In particular, we discuss the possibility to use integrated optics, proven in the near-infrared, in the thermal near-infrared (L and M bands, 3-5 $\mu$m), a sweet spot to image and characterize young extra-solar planetary systems. Finally, we address the science cases of a high-contrast VLTI imaging instrument and focus particularly on exoplanet science (young exoplanets, planet formation, and exozodiacal disks), stellar physics (fundamental parameters and multiplicity), and extragalactic astrophysics (active galactic nuclei and fundamental constants). Synergies and scientific preparation for other potential future instruments such as the Planet Formation Imager are also briefly discussed.
[15]  [pdf] - 1659573
Discovery of a point-like source and a third spiral arm in the transition disk around the Herbig Ae star MWC 758
Comments: 10 pages, 9 figures, Accepted for publication in Astronomy & Astrophysics
Submitted: 2017-10-31, last modified: 2017-12-01
Transition disks offer the extraordinary opportunity to look for newly born planets and investigate the early stages of planet formation. In this context we observed the Herbig A5 star MWC 758 with the L band vector vortex coronagraph installed in the near-infrared camera and spectrograph NIRC2 at the Keck II telescope, with the aim of unveiling the nature of the spiral structure by constraining the presence of planetary companions in the system. Our high-contrast imaging observations show a bright (delta L=7.0+/-0.3 mag) point-like emission, south of MWC 758 at a deprojected separation of about 20 au (r=0.111+/- 0. 004 arcsec) from the central star. We also recover the two spiral arms (south-east and north-west), already imaged by previous studies in polarized light, and discover a third one to the south-west of the star. No additional companions were detected in the system down to 5 Jupiter masses beyond 0.6 arcsec from the star. We propose that the bright L band emission could be caused by the presence of an embedded and accreting protoplanet, although the possibility of it being an asymmetric disk feature cannot be excluded. The spiral structure is probably not related to the protoplanet candidate, unless on an inclined and eccentric orbit, and it could be due to one (or more) yet undetected planetary companions at the edge of or outside the spiral pattern. Future observations and additional simulations will be needed to shed light on the true nature of the point-like source and its link with the spiral arms.
[16]  [pdf] - 1682455
M Dwarf Exoplanet Surface Density Distribution: A Log-Normal Fit from 0.07-400 AU
Comments: 4 pages with 2 figures. Accepted for publication in Astronomy and Astrophysics Letters
Submitted: 2017-07-17
We fit a log-normal function to the M dwarf orbital surface density distribution of gas giant planets, over the mass range 1-10 times that of Jupiter, from 0.07-400 AU. We use a Markov Chain Monte Carlo approach to explore the likelihoods of various parameter values consistent with point estimates of the data given our assumed functional form. This fit is consistent with radial velocity, microlensing, and direct imaging observations, is well-motivated from theoretical and phenomenological viewpoints, and makes predictions of future surveys. We present probability distributions for each parameter as well as a Maximum Likelihood Estimate solution. We suggest this function makes more physical sense than other widely used functions, and explore the implications of our results on the design of future exoplanet surveys.
[17]  [pdf] - 1581953
The VLT/NaCo large program to probe the occurrence of exoplanets and brown dwarfs at wide orbits. IV. Gravitational instability rarely forms wide, giant planets
Comments: 20 pages, 11 figures. Accepted for publication in A&A
Submitted: 2017-03-15
Understanding the formation and evolution of giant planets ($\ge$1 $M_{Jup}$) at wide orbital separation ($\ge$5 AU) is one of the goals of direct imaging. Over the past 15 years, many surveys have placed strong constraints on the occurrence rate of wide-orbit giants, mostly based on non-detections, but very few have tried to make a direct link with planet formation theories. In the present work, we combine the results of our previously published VLT/NaCo large program with the results of 12 past imaging surveys to constitute a statistical sample of 199 FGK stars within 100 pc, including 3 stars with sub-stellar companions. Using Monte Carlo simulations and assuming linear flat distributions for the mass and semi-major axis of planets, we estimate the sub-stellar companion frequency to be within 0.75-5.7% at the 68% confidence level (CL) within 20-300 AU and 0.5-75 $M_{Jup}$, which is compatible with previously published results. We also compare our results with the predictions of state-of-the-art population synthesis models based on the gravitational instability (GI) formation scenario by Forgan & Rice (2013), with and without scattering. We estimate that in both the scattered and non-scattered populations, we would be able to detect more than 30% of companions in the 1-75 $M_{Jup}$ range (95% CL). With the 3 sub-stellar detections in our sample, we estimate the fraction of stars that host a planetary system formed by GI to be within 1.0-8.6% (95% CL). We also conclude that even though GI is not common, it predicts a mass distribution of wide-orbit massive companions that is much closer to what is observed than what the core accretion scenario predicts. Finally, we associate the present paper with the release of the Direct Imaging Virtual Archive (DIVA,, a public database that aims at gathering the results of past, present, and future direct imaging surveys.
[18]  [pdf] - 1533197
Characterization of the inner disk around HD 141569 A from Keck/NIRC2 L-band vortex coronagraphy
Comments: 12 pages, 6 figures
Submitted: 2016-12-09
HD 141569 A is a pre-main sequence B9.5 Ve star surrounded by a prominent and complex circumstellar disk, likely still in a transition stage from protoplanetary to debris disk phase. Here, we present a new image of the third inner disk component of HD 141569 A made in the L' band (3.8 micron) during the commissioning of the vector vortex coronagraph recently installed in the near-infrared imager and spectrograph NIRC2 behind the W.M. Keck Observatory Keck II adaptive optics system. We used reference point spread function subtraction, which reveals the innermost disk component from the inner working distance of $\simeq 23$ AU and up to $\simeq 70$ AU. The spatial scale of our detection roughly corresponds to the optical and near-infrared scattered light, thermal Q, N and 8.6 micron PAH emission reported earlier. We also see an outward progression in dust location from the L'-band to the H-band (VLT/SPHERE image) to the visible (HST/STIS image), likely indicative of dust blowout. The warm disk component is nested deep inside the two outer belts imaged by HST NICMOS in 1999 (respectively at 406 and 245 AU). We fit our new L'-band image and spectral energy distribution of HD 141569 A with the radiative transfer code MCFOST. Our best-fit models favor pure olivine grains, and are consistent with the composition of the outer belts. While our image shows a putative very-faint point-like clump or source embedded in the inner disk, we did not detect any true companion within the gap between the inner disk and the first outer ring, at a sensitivity of a few Jupiter masses.
[19]  [pdf] - 1533199
The W. M. Keck Observatory infrared vortex coronagraph and a first image of HIP79124 B
Comments: 21 pages, 5 figures
Submitted: 2016-12-09
An optical vortex coronagraph has been implemented within the NIRC2 camera on the Keck II telescope and used to carry out on-sky tests and observations. The development of this new L'-band observational mode is described, and an initial demonstration of the new capability is presented: a resolved image of the low-mass companion to HIP79124, which had previously been detected by means of interferometry. With HIP79124 B at a projected separation of 186.5 mas, both the small inner working angle of the vortex coronagraph and the related imaging improvements were crucial in imaging this close companion directly. Due to higher Strehl ratios and more relaxed contrasts in L' band versus H band, this new coronagraphic capability will enable high-contrast small-angle observations of nearby young exoplanets and disks on a par with those of shorter-wavelength extreme adaptive optics coronagraphs.
[20]  [pdf] - 1531052
Three years of harvest with the vector vortex coronagraph in the thermal infrared
Comments: To appear in SPIE proceedings vol. 9908
Submitted: 2016-07-18
For several years, we have been developing vortex phase masks based on sub-wavelength gratings, known as Annular Groove Phase Masks. Etched onto diamond substrates, these AGPMs are currently designed to be used in the thermal infrared (ranging from 3 to 13 {\mu}m). Our AGPMs were first installed on VLT/NACO and VLT/VISIR in 2012, followed by LBT/LMIRCam in 2013 and Keck/NIRC2 in 2015. In this paper, we review the development, commissioning, on-sky performance, and early scientific results of these new coronagraphic modes and report on the lessons learned. We conclude with perspectives for future developments and applications.
[21]  [pdf] - 1354831
The VLT/NaCo large program to probe the occurrence of exoplanets and brown dwarfs at wide orbits. III. The frequency of brown dwarfs and giant planets as companions to solar-type stars
Comments: 9 pages, 3 figures. Accepted for publication by Astronomy & Astrophysics
Submitted: 2015-10-28
In recent years there have been many attempts to characterize the occurrence of stellar, BD and planetary-mass companions to solar-type stars, with the aim of constraining formation mechanisms. From RV observations a dearth of companions with masses between 10-40 MJup has been noticed at close separations, suggesting the possibility of a distinct formation mechanism for objects above and below this range. We present a model for the substellar companion mass function (CMF). It consists of the superposition of the planet and BD companion mass distributions, assuming that we can extrapolate the RV measured companion mass function for planets to larger separations and the stellar companion mass-ratio distribution over all separations into the BD mass regime. By using both the results of the VLT/NaCo large program and the complementary archive datasets that probe the occurrence of planets and BDs on wide orbits around solar-type stars, we place some constraints on the planet and BD distributions. We developed a MC simulation tool to predict the outcome of a given survey, depending on the shape of the orbital parameter distributions. Comparing the predictions with the results of the observations, we calculate how likely different models are and which can be ruled out. Current observations are consistent with the proposed model for the CMF, as long as a sufficiently small outer truncation radius is introduced for the planet separation distribution. The results of the direct imaging surveys searching for substellar companions around Sun-like stars are consistent with a combined substellar mass spectrum of planets and BDs. This mass distribution has a minimum between 10 and 50 MJup, in agreement with RV measurements. The dearth of objects in this mass range would naturally arise from the shape of the mass distribution, without the introduction of any distinct formation mechanism for BDs.
[22]  [pdf] - 1277149
Searching for gas giant planets on Solar System scales - A NACO/APP L'-band survey of A- and F-type Main Sequence stars
Comments: Accepted for publication in MNRAS, 8 pages, 4 figures
Submitted: 2015-08-03
We report the results of a direct imaging survey of A- and F-type main sequence stars searching for giant planets. A/F stars are often the targets of surveys, as they are thought to have more massive giant planets relative to solar-type stars. However, most imaging is only sensitive to orbital separations $>$30 AU, where it has been demonstrated that giant planets are rare. In this survey, we take advantage of the high-contrast capabilities of the Apodizing Phase Plate coronagraph on NACO at the Very Large Telescope. Combined with optimized principal component analysis post-processing, we are sensitive to planetary-mass companions (2 to 12 $M_{\rm Jup}$) at Solar System scales ($\leq$30 AU). We obtained data on 13 stars in L'-band and detected one new companion as part of this survey: an M$6.0\pm0.5$ dwarf companion around HD 984. We re-detect low-mass companions around HD 12894 and HD 20385, both reported shortly after the completion of this survey. We use Monte Carlo simulations to determine new constraints on the low-mass ($<$80 $M_{\rm Jup}$) companion frequency, as a function of mass and separation. Assuming solar-type planet mass and separation distributions, normalized to the planet frequency appropriate for A-stars, and the observed companion mass-ratio distribution for stellar companions extrapolated to planetary masses, we derive a truncation radius for the planetary mass companion surface density of $<$135 AU at 95% confidence.
[23]  [pdf] - 904108
Long-lasting dust rings in gas-rich disks: sculpting by single and multiple planets
Comments: Submitted to ApJ. Comments very welcome
Submitted: 2014-11-19
We propose a mechanism by which dust rings in protoplanetary disks can form and be long-lasting compared to gas rings. This involves the existence of a pressure maximum which traps dust either in between two gap-opening planets or at the outermost gap edge of a single or multiple planet system, combined with the decoupling of large dust particles from the gas. We perform 2D gas hydrodynamical simulations of disks with one and two giant planets which may open deep or partial gaps. A gas ring forms in between two planets such that the surface mass density is higher than on either side of it. This ring is a region of pressure maximum where we expect large grains, which are marginally coupled to the gas and would otherwise be subject to radial drift, to collect. Such a pressure maximum also occurs at the outermost gap edge in a disk with one or more planets. We infer the dust evolution in these regions as the gas disk evolves, to understand the longer term behavior of the resulting dust rings. Over time the gas surface density in the ring(s) decreases, which may cause the larger trapped particles to decouple. Consequently, these particles are expected to stay in ring structure(s) longer than the gas. For a Minimum Mass Solar Nebula model, we expect that millimeter and centimeter-sized grains in the outer O(10) au would most likely undergo this trapping and decoupling process.
[24]  [pdf] - 1216103
Discovery of a Companion Candidate in the HD169142 Transition Disk and the Possibility of Multiple Planet Formation
Comments: Accepted to ApJL, see also Biller et al. 2014
Submitted: 2014-08-04
We present L' and J-band high-contrast observations of HD169142, obtained with the VLT/NACO AGPM vector vortex coronagraph and the Gemini Planet Imager, respectively. A source located at 0".156+/-0".032 north of the host star (PA=7.4+/-11.3 degrees) appears in the final reduced L' image. At the distance of the star (~145 pc), this angular separation corresponds to a physical separation of 22.7+/-4.7 AU, locating the source within the recently resolved inner cavity of the transition disk. The source has a brightness of L'=12.2+/-0.5 mag, whereas it is not detected in the J band (J>13.8 mag). If its L' brightness arose solely from the photosphere of a companion and given the J-L' color constraints, it would correspond to a 28-32 MJupiter object at the age of the star, according to the COND models. Ongoing accretion activity of the star suggests, however, that gas is left in the inner disk cavity from which the companion could also be accreting. In this case the object could be lower in mass and its luminosity enhanced by the accretion process and by a circumplanetary disk. A lower mass object is more consistent with the observed cavity width. Finally, the observations enable us to place an upper limit on the L'-band flux of a second companion candidate orbiting in the disk annular gap at ~50 AU, as suggested by millimeter observations. If the second companion is also confirmed, HD169142 might be forming a planetary system, with at least two companions opening gaps and possibly interacting with each other.
[25]  [pdf] - 1215885
Imaging the Inner and Outer Gaps of the Pre-Transitional Disk of HD 169142 at 7 mm
Comments: Accepted by ApJ Letters, 16 pages, 3 figures, ApJ Letters 2014
Submitted: 2014-07-24
We present Very Large Array observations at 7 mm that trace the thermal emission of large dust grains in the HD 169142 protoplanetary disk. Our images show a ring of enhanced emission of radius ~25-30 AU, whose inner region is devoid of detectable 7 mm emission. We interpret this ring as tracing the rim of an inner cavity or gap, possibly created by a planet or a substellar companion. The ring appears asymmetric, with the western part significantly brighter than the eastern one. This azimuthal asymmetry is reminiscent of the lopsided structures that are expected to be produced as a consequence of trapping of large dust grains. Our observations also reveal an outer annular gap at radii from ~40 to ~70 AU. Unlike other sources, the radii of the inner cavity, the ring, and the outer gap observed in the 7 mm images, which trace preferentially the distribution of large (mm/cm sized) dust grains, coincide with those obtained from a previous near-infrared polarimetric image, which traces scattered light from small (micron- sized) dust grains. We model the broad-band spectral energy distribution and the 7 mm images to constrain the disk physical structure. From this modeling we infer the presence of a small (radius ~0.6 AU) residual disk inside the central cavity, indicating that the HD 169142 disk is a pre-transitional disk. The distribution of dust in three annuli with gaps in between them suggests that the disk in HD 169142 is being disrupted by at least two planets or substellar objects.
[26]  [pdf] - 1165918
The Universality of the Companion Mass Ratio Distribution
Comments: 3 pages, 1 figure. Accepted for publication in A&A (research note)
Submitted: 2013-04-11, last modified: 2013-04-23
We present new results regarding the companion mass-ratio distribution (CMRD) of stars, as a follow-up of our previous work. We used a maximum-likelihood-estimation method to re-derive the field CMRD power law avoiding dependence on the arbitrary binning. We also considered two new surveys of multiples in the field for solar-type stars and M dwarfs to test the universality of the CMRD. We found no significant differences in the CMRD for M dwarfs and solar-type stars compared with previous results over the common mass ratio and separation range. The new best-fit power law of the CMRD in the field, combining two previous sets of data, is $dN/dq \propto q^{\beta}$, with $\beta=0.25\pm0.29$.
[27]  [pdf] - 1166045
The Hubble Space Telescope Treasury Program on the Orion Nebula Cluster
Comments: Accepted for publication on the Astrophysical Journal Supplement Series, March 27, 2013
Submitted: 2013-04-17
The Hubble Space Telescope (HST) Treasury Program on the Orion Nebula Cluster has used 104 orbits of HST time to image the Great Orion Nebula region with the Advanced Camera for Surveys (ACS), the Wide-Field/Planetary Camera 2 (WFPC2) and the Near Infrared Camera and Multi Object Spectrograph (NICMOS) instruments in 11 filters ranging from the U-band to the H-band equivalent of HST. The program has been intended to perform the definitive study of the stellar component of the ONC at visible wavelengths, addressing key questions like the cluster IMF, age spread, mass accretion, binarity and cirumstellar disk evolution. The scanning pattern allowed to cover a contiguous field of approximately 600 square arcminutes with both ACS and WFPC2, with a typical exposure time of approximately 11 minutes per ACS filter, corresponding to a point source depth AB(F435W) = 25.8 and AB(F775W)=25.2 with 0.2 magnitudes of photometric error. We describe the observations, data reduction and data products, including images, source catalogs and tools for quick look preview. In particular, we provide ACS photometry for 3399 stars, most of them detected at multiple epochs, WFPC2 photometry for 1643 stars, 1021 of them detected in the U-band, and NICMOS JH photometry for 2116 stars. We summarize the early science results that have been presented in a number of papers. The final set of images and the photometric catalogs are publicly available through the archive as High Level Science Products at the STScI Multimission Archive hosted by the Space Telescope Science Institute.
[28]  [pdf] - 1165883
The binary companion mass ratio distribution: an imprint of the star formation process?
Comments: 8 pages, 4 figures, accepted for publication in MNRAS
Submitted: 2013-04-10
We explore the effects of dynamical evolution in dense clusters on the companion mass ratio distribution (CMRD) of binary stars. Binary systems are destroyed by interactions with other stars in the cluster, lowering the total binary fraction and significantly altering the initial semi-major axis distribution. However, the shape of the CMRD is unaffected by dynamics; an equal number of systems with high mass ratios are destroyed compared to systems with low mass ratios. We might expect a weak dependence of the survivability of a binary on its mass ratio because its binding energy is proportional to both the primary and secondary mass components of the system. However, binaries are broken up by interactions in which the perturbing star has a significantly higher energy (by a factor of >10, depending on the particular binary properties) than the binding energy of the binary, or through multiple interactions in the cluster. We therefore suggest that the shape of the observed binary CMRD is an outcome of the star formation process, and should be measured in preference to the distributions of orbital parameters, such as the semi-major axis distribution.
[29]  [pdf] - 1117353
Direct imaging constraints on planet populations detected by microlensing
Comments: Accepted for publication in A&A as Research Note, 3 pages
Submitted: 2012-03-16
Results from gravitational microlensing suggested the existence of a large population of free-floating planetary mass objects. The main conclusion from this work was partly based on constraints from a direct imaging survey. This survey determined upper limits for the frequency of stars that harbor giant exoplanets at large orbital separations. Aims. We want to verify to what extent upper limits from direct imaging do indeed constrain the microlensing results. We examine the current derivation of the upper limits used in the microlensing study and re-analyze the data from the corresponding imaging survey. We focus on the mass and semi-major axis ranges that are most relevant in context of the microlensing results. We also consider new results from a recent M-dwarf imaging survey as these objects are typically the host stars for planets detected by microlensing. We find that the upper limits currently applied in context of the microlensing results are probably underestimated. This means that a larger fraction of stars than assumed may harbor gas giant planets at larger orbital separations. Also, the way the upper limit is currently used to estimate the fraction of free-floating objects is not strictly correct. If the planetary surface density of giant planets around M-dwarfs is described as df_Planet ~ a^beta da, we find that beta ~ 0.5 - 0.6 is consistent with results from different observational studies probing semi-major axes between ~0.03 - 30 AU. Having a higher upper limit on the fraction of stars that may have gas giant planets at orbital separations probed by the microlensing data implies that more of the planets detected in the microlensing study are potentially bound to stars rather than free-floating. The current observational data are consistent with a rising planetary surface density for giant exoplanets around M-dwarfs out to ~30 AU.
[30]  [pdf] - 1083208
Quantitative Evidence for an Intrinsic Age Spread in the Orion Nebula Cluster
Comments: 12 pages, 12 figures, Accepted for publication in Astronomy and Astrophysics
Submitted: 2011-08-04
Aims. We present a study of the distribution of stellar ages in the Orion Nebula Cluster (ONC) based on accurate HST photometry taken from the HST Treasury Program observations of the ONC utilizing the most recent estimate of the cluster's distance (Menten et al. 2007). We investigate the presence of an intrinsic age spread in the region and a possible trend of age with the spatial distribution. Methods. We estimate the extinction and accretion luminosity towards each source by performing synthetic photometry on an empirical calibration of atmospheric models (Da Rio et al. 2010) using the package Chorizos (Maiz-Apellaniz 2004). The position of the sources in the HR-diagram is compared with different theoretical isochrones to estimate the mean cluster age and age dispersion. Through Monte Carlo simulations we quantify the amount of intrinsic age spread in the region, taking into account uncertainties on the distance, spectral type, extinction, unresolved binaries, accretion and photometric variability. Results. According to Siess et al. (2000) evolutionary models the mean age of the Cluster is 2.2 Myr with a scatter of few Myrs. With Monte Carlo simulations we find that the observed age spread is inconsistent with a coeval stellar population, but is in agreement with a star formation activity between 1.5 and 3.5 Myrs. We also observe light evidence for a trend of ages with spatial distribution.
[31]  [pdf] - 1077312
Binary Formation Mechanisms: Constraints from the Companion Mass Ratio Distribution
Comments: 9 pages, 9 figures, Accepted for publication in The Astrophysical Journal
Submitted: 2011-06-15
We present a statistical comparison of the mass ratio distribution of companions, as observed in different multiplicity surveys, to the most recent estimate of the single object mass function (Bochanski et al. 2010). The main goal of our analysis is to test whether or not the observed companion mass ratio distribution (CMRD) as a function of primary star mass and star formation environment is consistent with having been drawn from the field star IMF. We consider samples of companions for M dwarfs, solar type and intermediate mass stars, both in the field as well as clusters or associations, and compare them with populations of binaries generated by random pairing from the assumed IMF for a fixed primary mass. With regard to the field we can reject the hypothesis that the CMRD was drawn from the IMF for different primary mass ranges: the observed CMRDs show a larger number of equal-mass systems than predicted by the IMF. This is in agreement with fragmentation theories of binary formation. For the open clusters {\alpha} Persei and the Pleiades we also reject the IMF random- pairing hypothesis. Concerning young star-forming regions, currently we can rule out a connection between the CMRD and the field IMF in Taurus but not in Chamaeleon I. Larger and different samples are needed to better constrain the result as a function of the environment. We also consider other companion mass functions (CMF) and we compare them with observations. Moreover the CMRD both in the field and clusters or associations appears to be independent of separation in the range covered by the observations. Combining therefore the CMRDs of M and G primaries in the field and intermediate mass primary binaries in Sco OB2 for mass ratios, q = M2/M1, from 0.2 to 1, we find that the best chi-square fit follows a power law dN/dq \propto q^{\beta}, with {\beta} = -0.50 \pm 0.29, consistent with previous results.