Source author record

T. Henning

T. Henning appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

170works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

170 published item(s)

preprint2025arXiv

Multi-band infrared imaging reveals dusty spiral arcs around the binary B[e] star 3 Puppis

3 Puppis is the brightest known B[e] star. Recent work classifies this A-type object as a supergiant, yet the impact of its binarity on the circumstellar environment (CE) remains hard to characterize. To resolve its dusty region at 5-10 mas, we obtained mid-IR interferometric observations with VLTI/MATISSE over 3-12 μm. Because the (u,v) coverage supports imaging, we introduce a statistical interferometric-imaging workflow based on MiRA to generate averaged images: this systematic approach enables the selection of an optimal set of reconstructions, improving the robustness and fidelity of the recovered features. We also use SPARCO, an independent tool well suited to bright central objects embedded in fainter extended emission. Images from both tools in the L, M, and N bands agree and reveal an asymmetric, elongated feature ~17 mas (~10 au at 631 pc) southeast of the star with ~20% density contrast. A second northwest asymmetry and a skewed inner rim are detected. Simple geometric modelling, guided by the MATISSE images, constrains the morphology, location, and flux of the CE and its asymmetries. The images are consistent with earlier VLTI measurements but expose a more complex CE with large-scale clumps in the southeast and northwest parts of the disc. Hydrodynamic modelling indicates that tidal spiral-wake perturbations from the central binary, dynamically excited at Lindblad resonances in the circumbinary disc, best explain the radial extent and curvature of the elongated structures seen in all bands.

preprint2022arXiv

A scaled-up planetary system around a supernova progenitor

Virtually all known exoplanets reside around stars with $M<2.3~M_\odot$; to clarify if the dearth of planets around more massive stars is real, we launched the direct-imaging B-star Exoplanet Abundance STudy (BEAST) survey targeting B stars ($M>2.4~M_\odot$) in the young (5-20 Myr) Scorpius-Centaurus association (Sco-Cen). Here we present the case of a massive ($M \sim 9~M_\odot$) BEAST target, $μ^2$ Sco. Based on kinematic information, we found that $μ^2$ Sco is a member of a small group which we label Eastern Lower Scorpius, refining in turn the precision on stellar parameters. Around this star we identified a robustly detected substellar companion ($14.4\pm 0.8 M_J$) at a projected separation of $290\pm 10$ au, and a probable second object ($18.5\pm 1.5 M_J$) at $21\pm 1$ au. The planet-to-star mass ratios of these objects are similar to that of Jupiter to the Sun, and their irradiation is similar to those of Jupiter and Mercury, respectively. The two companions of $μ^2$ Sco are naturally added to the giant planet b Cen b recently discovered by BEAST; although slightly more massive than the deuterium burning limit, their properties resemble those of giant planets around less massive stars and they are better reproduced by a formation under a planet-like, rather than a star-like scenario. Irrespective of the (needed) confirmation of the inner companion, $μ^2$ Sco is the first star that would end its life as a supernova that hosts such a system. The tentative high frequency of BEAST discoveries shows that giant planets or small-mass brown dwarfs can form around B stars. When putting this finding in the context of core accretion and gravitational instability, we conclude that the current modeling of both mechanisms is not able to produce this kind of companion. BEAST will pave the way for the first time to an extension of these models to intermediate and massive stars. (abridged)

preprint2022arXiv

Discovery and mass measurement of the hot, transiting, Earth-sized planet GJ 3929 b

We report the discovery of GJ 3929 b, a hot Earth-sized planet orbiting the nearby M3.5 V dwarf star, GJ 3929 (G 180--18, TOI-2013). Joint modelling of photometric observations from TESS sectors 24 and 25 together with 73 spectroscopic observations from CARMENES and follow-up transit observations from SAINT-EX, LCOGT, and OSN yields a planet radius of $R_b = 1.150 +/- 0.040$ R$_{earth}$, a mass of $M_b = 1.21 +/- 0.42$ M$_{earth}$, and an orbital period of $P_b = 2.6162745 +/- 0.0000030$ d. The resulting density of $ρ_b= 4.4 +/- 1.6$ g/cm$^{-3}$ is compatible with the Earth's mean density of about 5.5 g/cm$^{-3}$. Due to the apparent brightness of the host star (J=8.7 mag) and its small size, GJ 3929 b is a promising target for atmospheric characterisation with the JWST. Additionally, the radial velocity data show evidence for another planet candidate with $P_{[c]} = 14.303 +/- 0.035$ d, which is likely unrelated to the stellar rotation period, $P_{rot} = 122+/-13$ d, which we determined from archival HATNet and ASAS-SN photometry combined with newly obtained TJO data.

preprint2022arXiv

Host galaxy magnitude of OJ 287 from its colours at minimum light

OJ 287 is a BL Lacertae type quasar in which the active galactic nucleus (AGN) outshines the host galaxy by an order of magnitude. The only exception to this may be at minimum light when the AGN activity is so low that the host galaxy may make quite a considerable contribution to the photometric intensity of the source. Such a dip or a fade in the intensity of OJ 287 occurred in November 2017, when its brightness was about 1.75 magnitudes lower than the recent mean level. We compare the observations of this fade with similar fades in OJ 287 observed earlier in 1989, 1999, and 2010. It appears that there is a relatively strong reddening of the B$-$V colours of OJ 287 when its V-band brightness drops below magnitude 17. Similar changes are also seen V$-$R, V$-$I, and R$-$I colours during these deep fades. These data support the conclusion that the total magnitude of the host galaxy is $V=18.0 \pm 0.3$, corresponding to $M_{K}=-26.5 \pm 0.3$ in the K-band. This is in agreement with the results, obtained using the integrated surface brightness method, from recent surface photometry of the host. These results should encourage us to use the colour separation method also in other host galaxies with strongly variable AGN nuclei. In the case of OJ 287, both the host galaxy and its central black hole are among the biggest known, and its position in the black hole mass-galaxy mass diagram lies close to the mean correlation.

preprint2022arXiv

In-depth direct imaging and spectroscopic characterization of the young Solar System analog HD 95086

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

preprint2022arXiv

The CARMENES search for exoplanets around M dwarfs: Two Saturn-mass planets orbiting active stars

The CARMENES radial-velocity survey is currently searching for planets in a sample of 387 M dwarfs. Here we report on two Saturn-mass planets orbiting TYC 2187-512-1 ($M_\star = 0.50 M_\odot$) and TZ Ari ($M_\star = 0.15 M_\odot$), respectively. We obtained supplementary photometric time series, which we use along with spectroscopic information to determine the rotation periods of the two stars. In both cases, the radial velocities also show strong modulations at the respective rotation period. We thus modeled the radial velocities as a Keplerian orbit plus a Gaussian process representing the stellar variability. TYC 2187-512-1 is found to harbor a planet with a minimum mass of 0.33 $M_{\rm Jup}$ in a near-circular 692-day orbit. The companion of TZ Ari has a minimum mass of 0.21 $M_{\rm Jup}$, orbital period of 771 d, and orbital eccentricity of 0.46. We provide an overview of all known giant planets in the CARMENES sample, from which we infer an occurrence rate of giant planets orbiting M dwarfs with periods up to 2 years in the range between 2% to 6%. TZ Ari b is only the second giant planet discovered orbiting a host with mass less than 0.3 $M_\odot$. These objects occupy an extreme location in the planet mass versus host mass plane. It is difficult to explain their formation in core-accretion scenarios, so they may possibly have been formed through a disk fragmentation process.

preprint2022arXiv

The SEDIGISM survey: Molecular cloud morphology. I. Classification and star formation

We present one of the very first extensive classifications of a large sample of molecular clouds based on their morphology. This is achieved using a recently published catalogue of 10663 clouds obtained from the first data release of the SEDIGISM survey. The clouds are classified into four different morphologies by visual inspection and using an automated algorithm -- J plots. The visual inspection also serves as a test for the J plots algorithm, as this is the first time it has been used on molecular gas. Generally, it has been found that the structure of molecular clouds is highly filamentary and our observations indeed verify that most of our molecular clouds are elongated structures. Based on our visual classification of the 10663 SEDIGISM clouds, 15% are ring-like, 57% are elongated, 15% are concentrated and 10% are clumpy clouds. The remaining clouds do not belong to any of these morphology classes and are termed unclassified. We compare the SEDIGISM molecular clouds with structures identified through other surveys, i.e. ATLASGAL elongated structures and the bubbles from Milky Way Project (MWP). We find that many of the ATLASGAL and MWP structures are velocity coherent. ATLASGAL elongated structures overlap with ~ 21% of the SEDIGISM elongated structures (elongated and clumpy clouds) and MWP bubbles overlap with ~ 25% of the SEDIGISM ring-like clouds. We also analyse the star-formation associated with different cloud morphologies using two different techniques. The first technique examines star formation efficiency (SFE) and the dense gas fraction (DGF), based on SEDIGISM clouds and ATLASGAL clumps data. The second technique uses the high-mass star formation (HMSF) threshold for molecular clouds. The results indicate that clouds with ring-like and clumpy morphologies show a higher degree of star formation.

preprint2022arXiv

The SEDIGISM survey: Molecular cloud morphology. II. Integrated source properties

The Structure, Excitation, and Dynamics of the Inner Galactic InterStellar Medium (SEDIGISM) survey has produced high (spatial and spectral) resolution $^{13}$CO (2-1) maps of the Milky Way. It has allowed us to investigate the molecular interstellar medium in the inner Galaxy at an unprecedented level of detail and characterise it into molecular clouds. In a previous paper, we have classified the SEDIGISM clouds into four morphologies. However, how the properties of the clouds vary for these four morphologies is not well understood. Here, we use the morphological classification of SEDIGISM clouds to find connections between the cloud morphologies, their integrated properties, and their location on scaling relation diagrams. We observe that ring-like clouds show the most peculiar properties, having, on average, higher masses, sizes, aspect ratios and velocity dispersions compared to other morphologies. We speculate that this is related to the physical mechanisms that regulate their formation and evolution, for example, turbulence from stellar feedback can often results in the creation of bubble-like structures. We also see a trend of morphology with virial parameter whereby ring-like, elongated, clumpy and concentrated clouds have virial parameters in a decreasing order. Our findings provide a foundation for a better understanding of the molecular cloud behaviour based on their measurable properties.

preprint2021arXiv

An extended scattered light disk around AT Pyx -- Possible planet formation in a cometary globule

To understand how the multitude of planetary systems that have been discovered come to be, we need to study systems at different evolutionary stages, with different central stars but also in different environments. The most challenging environment for planet formation may be the harsh UV radiation field of nearby massive stars which quickly erodes disks by external photo-evaporation. We have observed the AT Pyx system, located in the head of a cometary globule in the Gum Nebula, to search for signs of ongoing planet formation. We used the extreme adaptive optics imager VLT/SPHERE to observe AT Pyx in polarized light as well as total intensity in the J, H and K-band. Additionally we employed VLT/NACO to observe the system in the L-band. We resolve the disk around AT Pyx in scattered light across multiple wavelengths. We find an extended (>126 au) disk, with an intermediate inclination between 35 deg and 42 deg. The disk shows complex sub-structure and we identify 2 and possibly 3 spiral-like features. Depending on the precise geometry of the disk (which we can not unambiguously infer from our data) the disk may be eccentric with an eccentricity of ~0.16 or partially self-shadowed. The spiral features and possible eccentricity are both consistent with signatures of an embedded gas giant planet equal in mass to Jupiter. Our own observations can rule out brown dwarf companions embedded in the resolved disk, but are not sensitive enough to detect gas giants. AT Pyx is the first disk in a cometray globule in the Gum Nebula which is spatially resolved. By comparison with disks in the Orion Nebula Cluster we note that the extension of the disk may be exceptional for this environment if the external UV radiation field is comparable to other cometary globules in the region. The signposts of ongoing planet formation are intriguing and need to be followed up with higher sensitivity.

preprint2021arXiv

Constraining the Nature of the PDS 70 Protoplanets with VLTI/GRAVITY

We present K-band interferometric observations of the PDS 70 protoplanets along with their host star using VLTI/GRAVITY. We obtained K-band spectra and 100 $μ$as precision astrometry of both PDS 70 b and c in two epochs, as well as spatially resolving the hot inner disk around the star. Rejecting unstable orbits, we found a nonzero eccentricity for PDS 70 b of $0.17 \pm 0.06$, a near-circular orbit for PDS 70 c, and an orbital configuration that is consistent with the planets migrating into a 2:1 mean motion resonance. Enforcing dynamical stability, we obtained a 95% upper limit on the mass of PDS 70 b of 10 $M_\textrm{Jup}$, while the mass of PDS 70 c was unconstrained. The GRAVITY K-band spectra rules out pure blackbody models for the photospheres of both planets. Instead, the models with the most support from the data are planetary atmospheres that are dusty, but the nature of the dust is unclear. Any circumplanetary dust around these planets is not well constrained by the planets' 1-5 $μ$m spectral energy distributions (SEDs) and requires longer wavelength data to probe with SED analysis. However with VLTI/GRAVITY, we made the first observations of a circumplanetary environment with sub-au spatial resolution, placing an upper limit of 0.3~au on the size of a bright disk around PDS 70 b.

preprint2021arXiv

Deep Images of the Galactic Center with GRAVITY

Stellar orbits at the Galactic Center provide a very clean probe of the gravitational potential of the supermassive black hole. They can be studied with unique precision, beyond the confusion limit of a single telescope, with the near-infrared interferometer GRAVITY. Imaging is essential to search the field for faint, unknown stars on short orbits which potentially could constrain the black hole spin. Furthermore, it provides the starting point for astrometric fitting to derive highly accurate stellar positions. Here, we present $\mathrm{G^R}$, a new imaging tool specifically designed for Galactic Center observations with GRAVITY. The algorithm is based on a Bayesian interpretation of the imaging problem, formulated in the framework of information field theory and building upon existing works in radio-interferometric imaging. Its application to GRAVITY observations from 2021 yields the deepest images to date of the Galactic Center on scales of a few milliarcseconds. The images reveal the complicated source structure within the central $100\,\mathrm{mas}$ around Sgr A*, where we detected the stars S29 and S55 and confirm S62 on its trajectory, slowly approaching Sgr A*. Furthermore, we were able to detect S38, S42, S60, and S63 in a series of exposures for which we offset the fiber from Sgr A*. We provide an update on the orbits of all aforementioned stars. In addition to these known sources, the images also reveal a faint star moving to the west at a high angular velocity. We cannot find any coincidence with any known source and, thus, we refer to the new star as S300. From the flux ratio with S29, we estimate its K-band magnitude as $m_\mathrm{K}\left(\mathrm{S300}\right)\simeq 19.0 - 19.3$. Images obtained with CLEAN confirm the detection.

preprint2021arXiv

Detection of faint stars near SgrA* with GRAVITY

The spin of the supermassive black hole that resides at the Galactic Centre can in principle be measured by accurate measurements of the orbits of stars that are much closer to SgrA* than S2, the orbit of which recently provided the measurement of the gravitational redshift and the Schwarzschild precession. The GRAVITY near-infrared interferometric instrument combining the four 8m telescopes of the VLT provides a spatial resolution of 2-4 mas, breaking the confusion barrier for adaptive-optics-assisted imaging with a single 8-10m telescope. We used GRAVITY to observe SgrA* over a period of six months in 2019 and employed interferometric reconstruction methods developed in radio astronomy to search for faint objects near SgrA*. This revealed a slowly moving star of magnitude 18.9 in K band within 30mas of SgrA*. The position and proper motion of the star are consistent with the previously known star S62, which is at a substantially larger physical distance, but in projection passes close to SgrA*. Observations in August and September 2019 easily detected S29, with K-magnitude of 16.6, at approximately 130 mas from SgrA*. The planned upgrades of GRAVITY, and further improvements in the calibration, hold the promise of finding stars fainter than magnitude 19 at K.

preprint2021arXiv

First light of the VLT planet finder SPHERE. III. New spectrophotometry and astrometry of the HR8799 exoplanetary system

The planetary system discovered around the young A-type HR8799 provides a unique laboratory to: a) test planet formation theories, b) probe the diversity of system architectures at these separations, and c) perform comparative (exo)planetology. We present and exploit new near-infrared images and integral-field spectra of the four gas giants surrounding HR8799 obtained with SPHERE, the new planet finder instrument at the Very Large Telescope, during the commissioning and science verification phase of the instrument (July-December 2014). With these new data, we contribute to completing the spectral energy distribution of these bodies in the 1.0-2.5 $μ$m range. We also provide new astrometric data, in particular for planet e, to further constrain the orbits. We used the infrared dual-band imager and spectrograph (IRDIS) subsystem to obtain pupil-stabilized, dual-band $H2H3$ (1.593 $μ$m, 1.667 $μ$m), $K1K2$ (2.110 $μ$m, 2.251 $μ$m), and broadband $J$ (1.245 $μ$m) images of the four planets. IRDIS was operated in parallel with the integral field spectrograph (IFS) of SPHERE to collect low-resolution ($R\sim30$), near-infrared (0.94-1.64 $μ$m) spectra of the two innermost planets HR8799d and e. The data were reduced with dedicated algorithms, such as the Karhunen-Loève image projection (KLIP), to reveal the planets. We used the so-called negative planets injection technique to extract their photometry, spectra, and measure their positions. We illustrate the astrometric performance of SPHERE through sample orbital fits compatible with SPHERE and literature data.

preprint2021arXiv

Improved GRAVITY astrometric accuracy from modeling of optical aberrations

The GRAVITY instrument on the ESO VLTI pioneers the field of high-precision near-infrared interferometry by providing astrometry at the $10 - 100\,μ$as level. Measurements at such high precision crucially depend on the control of systematic effects. Here, we investigate how aberrations introduced by small optical imperfections along the path from the telescope to the detector affect the astrometry. We develop an analytical model that describes the impact of such aberrations on the measurement of complex visibilities. Our formalism accounts for pupil-plane and focal-plane aberrations, as well as for the interplay between static and turbulent aberrations, and successfully reproduces calibration measurements of a binary star. The Galactic Center observations with GRAVITY in 2017 and 2018, when both Sgr A* and the star S2 were targeted in a single fiber pointing, are affected by these aberrations at a level of less than 0.5 mas. Removal of these effects brings the measurement in harmony with the dual beam observations of 2019 and 2020, which are not affected by these aberrations. This also resolves the small systematic discrepancies between the derived distance $R_0$ to the Galactic Center reported previously.

preprint2021arXiv

Probing inner and outer disk misalignments in transition disks

For several transition disks (TDs), dark regions interpreted as shadows have been observed in scattered light imaging and are hypothesized to originate from misalignments between distinct disk regions. We aim to investigate the presence of misalignments in TDs. We study the inner disk geometries of 20 well-known transition disks with VLTI/GRAVITY observations and use complementary $^{12}$CO and $^{13}$CO molecular line data from ALMA to derive the orientation of the outer disk regions. We fit simple models to the GRAVITY data to derive the inner disks inclination and position angles. The outer disk geometries were derived from Keplerian fits to the ALMA velocity maps and compared to the inner disk constraints. We also predicted the locations of shadows for significantly misaligned systems. Our analysis reveals six disks to exhibit significant misalignments between their inner and outer disks. The predicted shadow positions agree well with the scattered light images of HD100453 and HD142527, and we find supporting evidence for a shadow in the disk around CQ Tau. In the other three targets for which we infer significantly misaligned disks, V1247 Ori, V1366 Ori, and RY Lup, we do not see any evident sign of shadows in the scattered light images. The scattered light shadows observed in DoAr44, HD135344B, and HD139614 are consistent with our observations, yet the underlying morphology is likely too complex to be described by our models and the accuracy achieved by our observations. Whereas we can derive precise constraints on the potential shadow positions for well-resolved inner disks around HAeBe stars, the statistical uncertainties for the marginally resolved inner disks around the TTS of our sample make it difficult to extract conclusive constraints for the presence of shadows in these systems.

preprint2021arXiv

Signs of late infall and possible planet formation around DR Tau using VLT/SPHERE and LBTI/LMIRCam

Context. Protoplanetary disks around young stars often contain substructures like rings, gaps, and spirals that could be caused by interactions between the disk and forming planets. Aims. We aim to study the young (1-3 Myr) star DR Tau in the near-infrared and characterize its disk, which was previously resolved through sub-millimeter interferometry with ALMA, and to search for possible sub-stellar companions embedded into it. Methods. We observed DR Tau with VLT/SPHERE both in polarized light (H broad band) and total intensity (in Y, J, H, and K spectral bands). We also performed L' band observations with LBTI/LMIRCam on the Large Binocular Telescope (LBT). Results. We found two previously undetected spirals extending north-east and south of the star, respectively. We further detected an arc-like structure north of the star. Finally a bright, compact and elongated structure was detected at separation of 303 +/- 10 mas and position angle 21.2 +/- 3.7 degrees, just at the root of the north-east spiral arm. Since this feature is visible both in polarized light and in total intensity and has a flat spectrum it is likely caused by stellar light scattered by dust. Conclusions. The two spiral arms are at different separation from the star, have very different pitch angles, and are separated by an apparent discontinuity, suggesting they might have a different origin. The very open southern spiral arm might be caused by infalling material from late encounters with cloudlets into the formation environment of the star itself. The compact feature could be caused by interaction with a planet in formation still embedded in its dust envelope and it could be responsible for launching the north-east spiral. We estimate a mass of the putative embedded object of the order of few M_Jup .

preprint2021arXiv

The ExoGRAVITY project: using single mode interferometry to characterize exoplanets

Combining adaptive optics and interferometric observations results in a considerable contrast gain compared to single-telescope, extreme AO systems. Taking advantage of this, the ExoGRAVITY project is a survey of known young giant exoplanets located in the range of 0.1'' to 2'' from their stars. The observations provide astrometric data of unprecedented accuracy, being crucial for refining the orbital parameters of planets and illuminating their dynamical histories. Furthermore, GRAVITY will measure non-Keplerian perturbations due to planet-planet interactions in multi-planet systems and measure dynamical masses. Over time, repetitive observations of the exoplanets at medium resolution ($R=500$) will provide a catalogue of K-band spectra of unprecedented quality, for a number of exoplanets. The K-band has the unique properties that it contains many molecular signatures (CO, H$_2$O, CH$_4$, CO$_2$). This allows constraining precisely surface gravity, metallicity, and temperature, if used in conjunction with self-consistent models like Exo-REM. Further, we will use the parameter-retrieval algorithm petitRADTRANS to constrain the C/O ratio of the planets. Ultimately, we plan to produce the first C/O survey of exoplanets, kick-starting the difficult process of linking planetary formation with measured atomic abundances.

preprint2021arXiv

The SEDIGISM survey: a search for molecular outflows

Context. The formation processes of massive stars are still unclear but a picture is emerging involving accretion disks and molecular outflows in what appears to be a scaled-up version of low-mass star formation. A census of outflow activity towards massive star-forming clumps in various evolutionary stages has the potential to shed light on massive star formation (MSF). Aims. We conducted an outflow survey towards ATLASGAL clumps using SEDIGISM data and aimed to obtain a large sample of clumps exhibiting outflows in different evolutionary stages. Methods. We identify the high-velocity wings of the 13CO lines toward ATLASGAL clumps by (1) extracting the simultaneously observed 13CO and C18O spectra from SEDIGISM, and (2) subtracting Gaussian fits to the scaled C18O from the 13CO, line after considering opacity broadening. Results. We have detected high-velocity gas towards 1192 clumps out of a total sample of 2052, giving an overall detection rate of 58%. Outflow activity has been detected in the earliest quiescent clumps (i.e., 70$μ$m weak), to the most evolved HII region stages i.e., 8$μ$m bright with MSF tracers. The detection rate increases as a function of evolution (quiescent=51%, protostellar=47%, YSO=57%, UCHII regions=76%). Conclusion. Our sample is the largest outflow sample identified so far. The high-detection rate from this large sample is consistent with previous results and supports that outflows are a ubiquitous feature of MSF. The lower detection rate in early evolutionary stages may be due to that outflows in the early stages are weak and difficult to detect. We obtain a statistically significant sample of outflow clumps for every evolutionary stage, especially for outflow clumps in the 70$μ$m dark stage. This suggests that the absence of 70$μ$m emission is not a robust indicator of starless/pre-stellar cores.

preprint2021arXiv

The SEDIGISM survey: The influence of spiral arms on the molecular gas distribution of the inner Milky Way

The morphology of the Milky Way is still a matter of debate. In order to shed light on uncertainties surrounding the structure of the Galaxy, in this paper, we study the imprint of spiral arms on the distribution and properties of its molecular gas. To do so, we take full advantage of the SEDIGISM survey that observed a large area of the inner Galaxy in the $^{13}$CO(2-1) line at an angular resolution of 28". We analyse the influences of the spiral arms by considering the features of the molecular gas emission as a whole across the longitude-velocity map built from the full survey. Additionally, we examine the properties of the molecular clouds in the spiral arms compared to the properties of their counterparts in the inter-arm regions. Through flux and luminosity probability distribution functions, we find that the molecular gas emission associated with the spiral arms does not differ significantly from the emission between the arms. On average, spiral arms show masses per unit length of $\sim10^5-10^6$ M$_{\odot} $kpc$^{-1}$. This is similar to values inferred from data sets in which emission distributions were segmented into molecular clouds. By examining the cloud distribution across the Galactic plane, we infer that the molecular mass in the spiral arms is a factor of 1.5 higher than that of the inter-arm medium, similar to what is found for other spiral galaxies in the local Universe. We observe that only the distributions of cloud mass surface densities and aspect ratio in the spiral arms show significant differences compared to those of the inter-arm medium; other observed differences appear instead to be driven by a distance bias. By comparing our results with simulations and observations of nearby galaxies, we conclude that the measured quantities would classify the Milky Way as a flocculent spiral galaxy, rather than as a grand-design one.

preprint2020arXiv

A triple star in disarray -- Multi-epoch observations of T Tauri with VLT-SPHERE and LBT-LUCI

T Tauri remains an enigmatic triple star for which neither the evolutionary state of the stars themselves, nor the geometry of the complex outflow system is completely understood. Eight-meter class telescopes equipped with state-of-the-art adaptive optics provide the spatial resolution necessary to trace tangential motion of features over a timescale of a few years, and they help to associate them with the different outflows. We used J-, H-, and K-band high-contrast coronagraphic imaging with VLT-SPHERE recorded between 2016 and 2018 to map reflection nebulosities and obtain high precision near-infrared (NIR) photometry of the triple star. We also present molecular hydrogen emission maps of the 1-0 S(1) line at 2.122 micron obtained with LBT-LUCI during its commissioning period at the end of 2016. The data reveal a number of new features in the system, some of which are seen in reflected light and some are seen in H2 emission; furthermore, they can all be associated with the main outflows. The tangential motion of the features provides compelling evidence that T Tauri Sb drives the southeast-northwest outflow. T Tauri Sb has recently faded probably because of increased extinction as it passes through the southern circumbinary disk. While T Tauri Sb is approaching periastron, T Tauri Sa instead has brightened and is detected in all our J-band imagery for the first time.

preprint2020arXiv

Cloud formation in the atomic and molecular phase: HI self absorption (HISA) towards a Giant Molecular Filament

Molecular clouds form from the atomic phase of the interstellar medium. However, characterizing the transition between the atomic and the molecular interstellar medium (ISM) is a difficult observational task. Here we address cloud formation processes by combining HSIA with molecular line data. One scenario proposed by numerical simulations is that the column density probability density functions (N-PDF) evolves from a log-normal shape at early times to a power-law-like shape at later times. In this paper, we study the cold atomic component of the giant molecular filament GMF38a (d=3.4 kpc, length$\sim230$ pc). We identify an extended HISA feature, which is partly correlated with the 13CO emission. The peak velocities of the HISA and 13CO observations agree well on the eastern side of the filament, whereas a velocity offset of approximately 4 km s$^{-1}$ is found on the western side. The sonic Mach number we derive from the linewidth measurements shows that a large fraction of the HISA, which is ascribed to the cold neutral medium (CNM), is at subsonic and transonic velocities. The column density of the CNM is on the order of 10$^{20}$ to 10$^{21}$ cm$^{-2}$. The column density of molecular hydrogen is an order of magnitude higher. The N-PDFs from HISA (CNM), HI emission (WNM+CNM), and 13CO (molecular component) are well described by log-normal functions, which is in agreement with turbulent motions being the main driver of cloud dynamics. The N-PDF of the molecular component also shows a power law in the high column-density region, indicating self-gravity. We suggest that we are witnessing two different evolutionary stages within the filament. The eastern subregion seems to be forming a molecular cloud out of the atomic gas, whereas the western subregion already shows high column density peaks, active star formation and evidence of related feedback processes.

preprint2020arXiv

Detection of the Schwarzschild precession in the orbit of the star S2 near the Galactic centre massive black hole

The star S2 orbiting the compact radio source Sgr A* is a precision probe of the gravitational field around the closest massive black hole (candidate). Over the last 2.7 decades we have monitored the star's radial velocity and motion on the sky, mainly with the SINFONI and NACO adaptive optics (AO) instruments on the ESO VLT, and since 2017, with the four-telescope interferometric beam combiner instrument GRAVITY. In this paper we report the first detection of the General Relativity (GR) Schwarzschild Precession (SP) in S2's orbit. Owing to its highly elliptical orbit (e = 0.88), S2's SP is mainly a kink between the pre-and post-pericentre directions of motion ~ +- 1 year around pericentre passage, relative to the corresponding Kepler orbit. The superb 2017-2019 astrometry of GRAVITY defines the pericentre passage and outgoing direction. The incoming direction is anchored by 118 NACO-AO measurements of S2's position in the infrared reference frame, with an additional 75 direct measurements of the S2-Sgr A* separation during bright states ('flares') of Sgr A*. Our 14-parameter model fits for the distance, central mass, the position and motion of the reference frame of the AO astrometry relative to the mass, the six parameters of the orbit, as well as a dimensionless parameter f_SP for the SP (f_SP = 0 for Newton and 1 for GR). From data up to the end of 2019 we robustly detect the SP of S2, del phi = 12' per orbital period. From posterior fitting and MCMC Bayesian analysis with different weighting schemes and bootstrapping we find f_SP = 1.10 +- 0.19. The S2 data are fully consistent with GR. Any extended mass inside S2's orbit cannot exceed ~ 0.1% of the central mass. Any compact third mass inside the central arcsecond must be less than about 1000 M_sun.

preprint2020arXiv

Dynamically important magnetic fields near the event horizon of Sgr A*

We study the time-variable linear polarisation of Sgr A* during a bright NIR flare observed with the GRAVITY instrument on July 28, 2018. Motivated by the time evolution of both the observed astrometric and polarimetric signatures, we interpret the data in terms of the polarised emission of a compact region ('hotspot') orbiting a black hole in a fixed, background magnetic field geometry. We calculated a grid of general relativistic ray-tracing models, created mock observations by simulating the instrumental response, and compared predicted polarimetric quantities directly to the measurements. We take into account an improved instrument calibration that now includes the instrument's response as a function of time, and we explore a variety of idealised magnetic field configurations. We find that the linear polarisation angle rotates during the flare, which is consistent with previous results. The hotspot model can explain the observed evolution of the linear polarisation. In order to match the astrometric period of this flare, the near horizon magnetic field is required to have a significant poloidal component, which is associated with strong and dynamically important fields. The observed linear polarisation fraction of $\simeq 30\%$ is smaller than the one predicted by our model ($\simeq 50\%$). The emission is likely beam depolarised, indicating that the flaring emission region resolves the magnetic field structure close to the black hole.

preprint2020arXiv

Gap, shadows, spirals, streamers: SPHERE observations of binary-disk interactions in GG Tau A

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

preprint2020arXiv

Large-scale changes of the cloud coverage in the $ε$ Indi Ba,Bb system

We present the results of 14 nights of \textit{I}-band photometric monitoring of the nearby brown dwarf binary, $ε$ Indi Ba,Bb. Observations were acquired over 2 months, and total close to 42 hours of coverage at a typically high cadence of 1.4 minutes. At a separation of just $0.7''$, we do not resolve the individual components, and so effectively treat the binary as if it were a single object. However, $ε$ Indi Ba (spectral type T1) is the brightest known T-type brown dwarf, and is expected to dominate the photometric signal. We typically find no strong variability associated with the target during each individual night of observing, but see significant changes in mean brightness - by as much as $0.10$ magnitudes - over the 2 months of the campaign. This strong variation is apparent on a timescale of at least 2 days. We detect no clear periodic signature, which suggests we may be observing the T1 brown dwarf almost pole-on, and the days-long variability in mean brightness is caused by changes in the large-scale structure of the cloud coverage. Dynamic clouds will very likely produce lightning, and complementary high cadence \textit{V}-band and H\textit{$α$} images were acquired to search for the emission signatures associated with stochastic "strikes." We report no positive detections for the target in either of these passbands.

preprint2020arXiv

Orbital and spectral characterization of the benchmark T-type brown dwarf HD 19467B

Context. Detecting and characterizing substellar companions for which the luminosity, mass, and age can be determined independently is of utter importance to test and calibrate the evolutionary models due to uncertainties in their formation mechanisms. HD 19467 is a bright and nearby star hosting a cool brown dwarf companion detected with RV and imaging, making it a valuable object for such studies. Aims. We aim to further characterize the orbital, spectral, and physical properties of the HD 19467 system. Methods. We present new high-contrast imaging data with the SPHERE and NaCo instruments. We also analyze archival data from HARPS, NaCo, HIRES, UVES, and ASAS. We also use proper motion data of the star from Hipparcos and Gaia. Results. We refine the properties of the host star and derive an age of 8.0$^{+2.0}_{-1.0}$ Gyr based on isochrones, gyrochronology, and chemical and kinematic arguments. This estimate is slightly younger than previous estimates of ~9-11 Gyr. No orbital curvature is seen in the current imaging, RV, and astrometric data. From a joint fit of the data, we refine the orbital parameters for HD 19467B: period 398$^{+95}_{-93}$ yr, inclination 129.8$^{+8.1}_{-5.1}$ deg, eccentricity 0.56$\pm$0.09, longitude of the ascending node 134.8$\pm$4.5 deg, and argument of the periastron 64.2$^{+5.5}_{-6.3}$ deg. We assess a dynamical mass of 74$^{+12}_{-9}$ MJ. The fit with atmospheric models of the spectrophotometric data of HD 19467B indicates an atmosphere without clouds or with very thin clouds, an effective temperature of 1042$^{+77}_{-71}$ K, and a large surface gravity of 5.34$^{+0.08}_{-0.09}$ dex. The comparison to model predictions of the bolometric luminosity and dynamical mass of HD 19467B, assuming our system age estimate, indicates a better agreement with the Burrows et al. models; whereas the other evolutionary models used tend to underestimate its cooling rate.

preprint2020arXiv

Precise mass and radius of a transiting super-Earth planet orbiting the M dwarf TOI-1235: a planet in the radius gap?

We report the confirmation of a transiting planet around the bright, inactive M0.5 V star TOI-1235 (TYC 4384-1735-1, V = 11.5 mag), whose transit signal was detected in the photometric time series of Sectors 14, 20, and 21 of the TESS space mission. We confirm the planetary nature of the transit signal, which has a period of 3.44 d, by using precise radial velocity measurements with CARMENES and HARPS-N spectrographs. A comparison of the properties derived for TOI-1235 b's with theoretical models reveals that the planet has a rocky composition, with a bulk density slightly higher than Earth's. In particular, we measure a mass of M_p = 5.9+/-0.6 M_Earth and a radius of R_p = 1.69+/-0.08 R_Earth, which together result in a density of rho_p = 6.7+1.3-1.1 g/cm3. When compared with other well-characterized exoplanetary systems, the particular combination of planetary radius and mass puts our discovery in the radius gap, a transition region between rocky planets and planets with significant atmospheric envelopes, with few known members. While the exact location of the radius gap for M dwarfs is still a matter of debate, our results constrain it to be located at around 1.7 R_Earth or larger at the insolation levels received by TOI-1235 b (~60 S_Earth), which makes it an extremely interesting object for further studies of planet formation and atmospheric evolution.

preprint2020arXiv

Searching for the near infrared counterpart of Proxima c using multi-epoch high contrast SPHERE data at VLT

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

preprint2020arXiv

SPHERE+: Imaging young Jupiters down to the snowline

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.

preprint2020arXiv

Spirals inside the millimeter cavity of transition disk SR 21

Hydrodynamical simulations of planet-disk interactions suggest that planets may be responsible for a number of the sub-structures frequently observed in disks in both scattered light and dust thermal emission. Despite the ubiquity of these features, direct evidence of planets embedded in disks and of the specific interaction features like spiral arms within planetary gaps still remain rare. In this study we discuss recent observational results in the context of hydrodynamical simulations in order to infer the properties of a putative embedded planet in the cavity of a transition disk. We imaged the transition disk SR 21 in H-band in scattered light with SPHERE/IRDIS and in thermal dust emission with ALMA band 3 (3mm) observations at a spatial resolution of 0.1". We combine these datasets with existing band 9 (430um) and band 7 (870um) ALMA continuum data. The Band 3 continuum data reveals a large cavity and a bright ring peaking at 53 au strongly suggestive of dust trapping.The ring shows a pronounced azimuthal asymmetry, with a bright region in the north-west that we interpret as a dust over-density. A similarly-asymmetric ring is revealed at the same location in polarized scattered light, in addition to a set of bright spirals inside the mm cavity and a fainter spiral bridging the gap to the outer ring. These features are consistent with a number of previous hydrodynamical models of planet-disk interactions, and suggest the presence of a ~1 MJup planet at 44 au and PA=11°. This makes SR21 the first disk showing spiral arms inside the mm cavity, as well as one for which the location of a putative planet can be precisely inferred. With the location of a possible planet being well-constrained by observations, it is an ideal candidate for follow-up observations to search for direct evidence of a planetary companion still embedded in its disk.

preprint2020arXiv

The CARMENES search for exoplanets around M dwarfs. A super-Earth planet orbiting HD 79211 (GJ 338 B)

We report on radial velocity time series for two M0.0V stars, GJ338B and GJ338A, using the CARMENES spectrograph, complemented by ground-telescope photometry from Las Cumbres and Sierra Nevada observatories. We aim to explore the presence of small planets in tight orbits using the spectroscopic radial velocity technique. We obtained 159 and 70 radial velocity measurements of GJ338B and A, respectively, with the CARMENES visible channel. We also compiled additional relative radial velocity measurements from the literature and a collection of astrometric data that cover 200 a of observations to solve for the binary orbit. We found dynamical masses of 0.64$\pm$0.07M$_\odot$ for GJ338B and 0.69$\pm$0.07M$_\odot$ for GJ338A. The CARMENES radial velocity periodograms show significant peaks at 16.61$\pm$0.04 d (GJ338B) and 16.3$^{+3.5}_{-1.3}$ d (GJ338A), which have counterparts at the same frequencies in CARMENES activity indicators and photometric light curves. We attribute these to stellar rotation. GJ338B shows two additional, significant signals at 8.27$\pm$0.01 and 24.45$\pm$0.02 d, with no obvious counterparts in the stellar activity indices. The former is likely the first harmonic of the star's rotation, while we ascribe the latter to the existence of a super-Earth planet with a minimum mass of 10.27$^{+1.47}_{-1.38}$$M_{\oplus}$ orbiting GJ338B. GJ338B b lies inside the inner boundary of the habitable zone around its parent star. It is one of the least massive planets ever found around any member of stellar binaries. The masses, spectral types, brightnesses, and even the rotational periods are very similar for both stars, which are likely coeval and formed from the same molecular cloud, yet they differ in the architecture of their planetary systems.

preprint2020arXiv

The CARMENES search for exoplanets around M dwarfs. Characterization of the nearby ultra-compact multiplanetary system YZ Ceti

The nearby ultra-compact multiplanetary system YZ Ceti consists of at least three planets. The orbital period of each planet is the subject of discussion in the literature due to strong aliasing in the radial velocity data. The stellar activity of this M dwarf also hampers significantly the derivation of the planetary parameters. With an additional 229 radial velocity measurements obtained since the discovery publication, we reanalyze the YZ Ceti system and resolve the alias issues. We use model comparison in the framework of Bayesian statistics and periodogram simulations based on a method by Dawson and Fabrycky to resolve the aliases. We discuss additional signals in the RV data, and derive the planetary parameters by simultaneously modeling the stellar activity with a Gaussian process regression model. To constrain the planetary parameters further we apply a stability analysis on our ensemble of Keplerian fits. We resolve the aliases: the three planets orbit the star with periods of $2.02$ d, $3.06$ d, and $4.66$ d. We also investigate an effect of the stellar rotational signal on the derivation of the planetary parameters, in particular the eccentricity of the innermost planet. Using photometry we determine the stellar rotational period to be close to $68$ d. From the absence of a transit event with TESS, we derive an upper limit of the inclination of $i_\mathrm{max} = 87.43$ deg. YZ Ceti is a prime example of a system where strong aliasing hindered the determination of the orbital periods of exoplanets. Additionally, stellar activity influences the derivation of planetary parameters and modeling them correctly is important for the reliable estimation of the orbital parameters in this specific compact system. Stability considerations then allow additional constraints to be placed on the planetary parameters.

preprint2020arXiv

The CARMENES search for exoplanets around M dwarfs. Measuring precise radial velocities in the near infrared: the example of the super-Earth CD Cet b

The high-resolution, dual channel, visible and near-infrared spectrograph CARMENES offers exciting opportunities for stellar and exoplanetary research on M dwarfs. In this work we address the challenge of reaching the highest radial velocity precision possible with a complex, actively cooled, cryogenic instrument, such as the near-infrared channel. We describe the performance of the instrument and the work flow used to derive precise Doppler measurements from the spectra. The capability of both CARMENES channels to detect small exoplanets is demonstrated with the example of the nearby M5.0 V star CD Cet (GJ 1057), around which we announce a super-Earth ($4.0\pm0.4\,M_\oplus$) companion on a 2.29 d orbit.

preprint2020arXiv

The CARMENES search for exoplanets around M dwarfs. Radial velocities and activity indicators from cross-correlation functions with weighted binary masks

For years, the standard procedure to measure radial velocities (RVs) of spectral observations consisted in cross-correlating the spectra with a binary mask, that is, a simple stellar template that contains information on the position and strength of stellar absorption lines. The cross-correlation function (CCF) profiles also provide several indicators of stellar activity. We present a methodology to first build weighted binary masks and, second, to compute the CCF of spectral observations with these masks from which we derive radial velocities and activity indicators. These methods are implemented in a python code that is publicly available. To build the masks, we selected a large number of sharp absorption lines based on the profile of the minima present in high signal-to-noise ratio (S/N) spectrum templates built from observations of reference stars. We computed the CCFs of observed spectra and derived RVs and the following three standard activity indicators: full-width-at-half-maximum as well as contrast and bisector inverse slope.We applied our methodology to CARMENES high-resolution spectra and obtain RV and activity indicator time series of more than 300 M dwarf stars observed for the main CARMENES survey. Compared with the standard CARMENES template matching pipeline, in general we obtain more precise RVs in the cases where the template used in the standard pipeline did not have enough S/N. We also show the behaviour of the three activity indicators for the active star YZ CMi and estimate the absolute RV of the M dwarfs analysed using the CCF RVs.

preprint2020arXiv

The circumstellar environment of EX Lup: the SPHERE and SINFONI views

EX Lup is a well-studied T Tauri star that represents the prototype of young eruptive stars EXors. In this paper we analyze new adaptive optics imaging and spectroscopic observations of EX Lup and its circumstellar environment in near-infrared in its quiescent phase. We observed EX Lup with the high contrast imager SPHERE/IRDIS in the dual-beam polarimetric imaging mode to resolve the circumstellar environment in near-infrared scattered light. We complemented these data with earlier SINFONI spectroscopy. We resolve for the first time in scattered light a compact feature around EX Lup azimuthally extending from 280deg to 360deg, and radially extending from 0.3arcsec to 0.55arcsec in the plane of the disk. We explore two different scenarios for the detected emission. One accounting for the emission as coming from the brightened walls of the cavity excavated by the outflow whose presence was suggested by ALMA observations in the J=3-2 line of 12CO. The other accounts for the emission as coming from an inclined disk. We detect for the first time a more extended circumstellar disk in scattered light, which shows that a region between 10 and 30 au is depleted of mum-size grains. We compare the J-, H- and K-band spectra obtained with SINFONI in quiescence with the spectra taken during the outburst, showing that all the emission lines were due to the episodic accretion event. Conclusions. Based on the morphology analysis we favour the scenario in which the scattered light is coming from a circumstellar disk rather than the outflow around EX Lup. We analyze the origin of the observed feature either as coming from a continuous circumstellar disk with a cavity, or from the illuminated wall of the outer disk or from a shadowed disk. Moreover, we discuss what is the origin of the mum-size grains depleted region, exploring the possibility that a sub-stellar companion may be the cause of it.

preprint2020arXiv

The Flux Distribution of Sgr A*

The Galactic Center black hole Sagittarius A* is a variable NIR source that exhibits bright flux excursions called flares. The low-flux density turnover of the flux distribution is below the sensitivity of current single-aperture telescopes. We use the unprecedented resolution of the GRAVITY instrument at the VLTI. Our light curves are unconfused, overcoming the confusion limit of previous photometric studies. We analyze the light curves using standard statistical methods and obtain the flux distribution. We find that the flux distribution of SgrA* turns over at a median flux density of (1.1\pm0.3)mJy. We measure the percentiles of the flux distribution and use them to constrain the NIR K-band SED. Furthermore, we find that the flux distribution is intrinsically right-skewed to higher flux density in log space. Flux densities below 0.1mJy are hardly ever observed. In consequence, a single powerlaw or lognormal distribution does not suffice to describe the observed flux distribution in its entirety. However, if one takes into account a power law component at high flux densities, a lognormal distribution can describe the lower end of the observed flux distribution. We confirm the RMS-flux relation for Sgr~A* and find it to be linear for all flux densities in our observation. We conclude that Sgr~A* has two states: the bulk of the emission is generated in a lognormal process with a well-defined median flux density and this quiescent emission is supplemented by sporadic flares that create the observed power law extension of the flux distribution.

preprint2020arXiv

The GRAVITY Young Stellar Object survey IV. The CO overtone emission in 51 Oph at sub-au scales

51 Oph is a Herbig Ae/Be star that exhibits strong near-infrared CO ro-vibrational emission at 2.3 micron, most likely originating in the innermost regions of a circumstellar disc. We aim to obtain the physical and geometrical properties of the system by spatially resolving the circumstellar environment of the inner gaseous disc. We used the second-generation VLTI/GRAVITY to spatially resolve the continuum and the CO overtone emission. We obtained data over 12 baselines with the auxiliary telescopes and derive visibilities, and the differential and closure phases as a function of wavelength. We used a simple LTE ring model of the CO emission to reproduce the spectrum and CO line displacements. Our interferometric data show that the star is marginally resolved at our spatial resolution, with a radius of 10.58+-2.65 Rsun.The K-band continuum emission from the disc is inclined by 63+-1 deg, with a position angle of 116+-1 deg, and 4+-0.8 mas (0.5+-0.1 au) across. The visibilities increase within the CO line emission, indicating that the CO is emitted within the dust-sublimation radius.By modelling the CO bandhead spectrum, we derive that the CO is emitted from a hot (T=1900-2800 K) and dense (NCO=(0.9-9)x10^21 cm^-2) gas. The analysis of the CO line displacement with respect to the continuum allows us to infer that the CO is emitted from a region 0.10+-0.02 au across, well within the dust-sublimation radius. The inclination and position angle of the CO line emitting region is consistent with that of the dusty disc. Our spatially resolved interferometric observations confirm the CO ro-vibrational emission within the dust-free region of the inner disc. Conventional disc models exclude the presence of CO in the dust-depleted regions of Herbig AeBe stars. Ad hoc models of the innermost disc regions, that can compute the properties of the dust-free inner disc, are therefore required.

preprint2020arXiv

The GRAVITY young stellar object survey. II. First spatially resolved observations of the CO bandhead emission in a high-mass YSO

The inner regions of the discs of high-mass young stellar objects (HMYSOs) are still poorly known due to the small angular scales and the high visual extinction involved. We deploy near-infrared (NIR) spectro-interferometry to probe the inner gaseous disc in HMYSOs and investigate the origin and physical characteristics of the CO bandhead emission (2.3-2.4 $μ$m). We present the first GRAVITY/VLTI observations at high spectral (R=4000) and spatial (mas) resolution of the CO overtone transitions in NGC 2024 IRS2. The continuum emission is resolved in all baselines and is slightly asymmetric, displaying small closure phases ($\leq$8$^{\circ}$). Our best ellipsoid model provides a disc inclination of 34$^{\circ}$$\pm$1$^{\circ}$, a disc major axis position angle of 166$^{\circ}$$\pm$1$^{\circ}$, and a disc diameter of 3.99$\pm$0.09 mas (or 1.69$\pm$0.04 au, at a distance of 423 pc). The small closure phase signals in the continuum are modelled with a skewed rim, originating from a pure inclination effect. For the first time, our observations spatially and spectrally resolve the first four CO bandheads. Changes in visibility, as well as differential and closure phases across the bandheads are detected. Both the size and geometry of the CO-emitting region are determined by fitting a bidimensional Gaussian to the continuum-compensated CO bandhead visibilities. The CO-emitting region has a diameter of 2.74$\pm^{0.08}_{0.07}$ mas (1.16$\pm$0.03 au), and is located in the inner gaseous disc, well within the dusty rim, with inclination and $PA$ matching the dusty disc geometry, which indicates that both dusty and gaseous discs are coplanar. Physical and dynamical gas conditions are inferred by modelling the CO spectrum. Finally, we derive a direct measurement of the stellar mass of $M_*\sim$14.7$^{+2}_{-3.6}$ M$_{\odot}$ by combining our interferometric and spectral modelling results.

preprint2020arXiv

The search for disks or planetary objects around directly imaged companions: A candidate around DH Tau B

In recent decades, thousands of substellar companions have been discovered with both indirect and direct methods of detection. In this paper, we focus our attention on substellar companions detected with the direct imaging technique, with the primary goal of investigating their close surroundings and looking for additional companions and satellites, as well as disks and rings. Any such discovery would shed light on many unresolved questions, particularly with regard to their possible formation mechanisms. To reveal bound features of directly imaged companions we need to suppress the contribution from the source itself. Therefore, we developed a method based on the negative fake companion (NEGFC) technique that first estimates the position in the field of view (FoV) and the flux of the imaged companion, then subtracts a rescaled model point spread function (PSF) from the imaged companion. Next it performs techniques, such as angular differential imaging (ADI), to further remove quasi-static patterns of the star. We applied the method to the sample of substellar objects observed with SPHERE during the SHINE GTO survey. Among the 27 planets and brown dwarfs we analyzed, we detected a possible point source close to DH Tau B. This candidate companion was detected in four different SPHERE observations, with an estimated mass of $\sim 1$ M\textsubscript{Jup}, and a mass ratio with respect to the brown dwarf of $1/10$. This binary system, if confirmed, would be the first of its kind, opening up interesting questions for the formation mechanism, evolution, and frequency of such pairs. In order to address the latter, the residuals and contrasts reached for 25 companions in the sample of substellar objects observed with SPHERE were derived. If the DH Tau Bb companion is real, the binary fraction obtained is $\sim 7\%$, which is in good agreement with the results obtained for field brown dwarfs.

preprint2020arXiv

Three Short Period Jupiters from TESS

We report the confirmation and mass determination of three hot Jupiters discovered by the Transiting Exoplanet Survey Satellite (TESS) mission: HIP 65Ab (TOI-129, TIC-201248411) is an ultra-short-period Jupiter orbiting a bright (V=11.1 mag) K4-dwarf every 0.98 days. It is a massive 3.213 +/- 0.078 Mjup planet in a grazing transit configuration with an impact parameter of b = 1.17 +0.10/-0.08. As a result the radius is poorly constrained, 2.03 +0.61/-0.49 Rjup. The planet's distance to its host star is less than twice the separation at which it would be destroyed by Roche lobe overflow. It is expected to spiral into HIP 65A on a timescale ranging from 80 Myr to a few gigayears, assuming a reduced tidal dissipation quality factor of Qs' = 10^7 - 10^9. We performed a full phase-curve analysis of the TESS data and detected both illumination- and ellipsoidal variations as well as Doppler boosting. HIP 65A is part of a binary stellar system, with HIP 65B separated by 269 AU (3.95 arcsec on sky). TOI-157b (TIC 140691463) is a typical hot Jupiter with a mass of 1.18 +/- 0.13 Mjup and a radius of 1.29 +/- 0.02 Rjup. It has a period of 2.08 days, which corresponds to a separation of just 0.03 AU. This makes TOI-157 an interesting system, as the host star is an evolved G9 sub-giant star (V=12.7). TOI-169b (TIC 183120439) is a bloated Jupiter orbiting a V=12.4 G-type star. It has a mass of 0.79 +/- 0.06 Mjup and a radius of 1.09 +0.08/-0.05 Rjup. Despite having the longest orbital period (P = 2.26 days) of the three planets, TOI-169b receives the most irradiation and is situated on the edge of the Neptune desert. All three host stars are metal rich with [Fe/H] ranging from 0.18 - 0.24.

preprint2020arXiv

Ubiquitous velocity fluctuations throughout the molecular interstellar medium

The density structure of the interstellar medium (ISM) determines where stars form and release energy, momentum, and heavy elements, driving galaxy evolution. Density variations are seeded and amplified by gas motion, but the exact nature of this motion is unknown across spatial scale and galactic environment. Although dense star-forming gas likely emerges from a combination of instabilities, convergent flows, and turbulence, establishing the precise origin is challenging because it requires quantifying gas motion over many orders of magnitude in spatial scale. Here we measure the motion of molecular gas in the Milky Way and in nearby galaxy NGC 4321, assembling observations that span an unprecedented spatial dynamic range ($10^{-1}{-}10^3$ pc). We detect ubiquitous velocity fluctuations across all spatial scales and galactic environments. Statistical analysis of these fluctuations indicates how star-forming gas is assembled. We discover oscillatory gas flows with wavelengths ranging from $0.3{-}400$ pc. These flows are coupled to regularly-spaced density enhancements that likely form via gravitational instabilities. We also identify stochastic and scale-free velocity and density fluctuations, consistent with the structure generated in turbulent flows. Our results demonstrate that ISM structure cannot be considered in isolation. Instead, its formation and evolution is controlled by nested, interdependent flows of matter covering many orders of magnitude in spatial scale.

preprint2019arXiv

Demographics of disks around young very low-mass stars and brown dwarfs in Lupus

We present new 890 $μm$ continuum ALMA observations of 5 brown dwarfs (BDs) with infrared excess in Lupus I and III -- which, in combination with 4 BDs previously observed, allowed us to study the mm properties of the full known BD disk population of one star-forming region. Emission is detected in 5 out of the 9 BD disks. Dust disk mass, brightness profiles and characteristic sizes of the BD population are inferred from continuum flux and modeling of the observations. Only one source is marginally resolved, allowing for the determination of its disk characteristic size. We conduct a demographic comparison between the properties of disks around BDs and stars in Lupus. Due to the small sample size, we cannot confirm or disprove if the disk mass over stellar mass ratio drops for BDs, as suggested for Ophiuchus. Nevertheless, we find that all detected BD disks have an estimated dust mass between 0.2 and 3.2 $M_{\bigoplus}$; these results suggest that the measured solid masses in BD disks can not explain the observed exoplanet population, analogous to earlier findings on disks around more massive stars. Combined with the low estimated accretion rates, and assuming that the mm-continuum emission is a reliable proxy for the total disk mass, we derive ratios of $\dot{M}_{\mathrm{acc}} / M_{\mathrm{disk}}$ significantly lower than in disks around more massive stars. If confirmed with more accurate measurements of disk gas masses, this result could imply a qualitatively different relationship between disk masses and inward gas transport in BD disks.

preprint2019arXiv

Feedback in W49A diagnosed with radio recombination lines and models

We present images of radio recombination lines (RRLs) at wavelengths around 17 cm from the star-forming region W49A to determine the kinematics of ionized gas in the THOR survey (The HI/OH/Recombination line survey of the inner Milky Way) at an angular resolution of 16.8"x13.8". The distribution of ionized gas appears to be affected by feedback processes from the star clusters in W49A. The velocity structure of the RRLs shows a complex behavior with respect to the molecular gas. We find a shell-like distribution of ionized gas as traced by RRL emission surrounding the central cluster of OB stars in W49A. We describe the evolution of the shell with the recent feedback model code WARPFIELD that includes the important physical processes and has previously been applied to the 30 Doradus region in the Large Magellanic Cloud. The cloud structure and dynamics of W49A are in agreement with a feedback-driven shell that is re-collapsing. The shell may have triggered star formation in other parts of W49A. We suggest that W49A is a potential candidate for star formation regulated by feedback-driven and re-collapsing shells.

preprint2019arXiv

Peering into the formation history of beta Pictoris b with VLTI/GRAVITY long baseline interferometry

Our objective is to estimate the C/O ratio in the atmosphere of beta Pictoris b and obtain an estimate of the dynamical mass of the planet, as well as to refine its orbital parameters using high-precision astrometry. We used the GRAVITY instrument with the four 8.2 m telescopes of the Very Large Telescope Interferometer to obtain K-band spectro-interferometric data on $β$ Pic b. We extracted a medium resolution (R=500) K-band spectrum of the planet and a high-precision astrometric position. We estimated the planetary C/O ratio using two different approaches (forward modeling and free retrieval) from two different codes (ExoREM and petitRADTRANS, respectively). Finally, we used a simplified model of two formation scenarios (gravitational collapse and core-accretion) to determine which can best explain the measured C/O ratio. Our new astrometry disfavors a circular orbit for $β$ Pic b ($e=0.15^{+0.05}_{-0.04}$). Combined with previous results and with Hipparcos/GAIA measurements, this astrometry points to a planet mass of $M = 12.7\pm{}2.2\,M_\mathrm{Jup}$. This value is compatible with the mass derived with the free-retrieval code petitRADTRANS using spectral data only. The forward modeling and free-retrieval approches yield very similar results regarding the atmosphere of beta Pic b. In particular, the C/O ratios derived with the two codes are identical ($0.43\pm{}0.05$ vs $0.43^{+0.04}_{-0.03}$). We argue that if the stellar C/O in $β$ Pic is Solar, then this combination of a very high mass and a low C/O ratio for the planet suggests a formation through core-accretion, with strong planetesimal enrichment.

preprint2019arXiv

Shadowing and multiple rings in the protoplanetary disk of HD 139614

Shadows in scattered light images of protoplanetary disks are a common feature and support the presence of warps or misalignments between disk regions. These warps are possibly due to an inclined (sub-)stellar companion embedded in the disk. We study the morphology of the protoplanetary disk around the Herbig Ae star HD 139614 based on the first scattered light observations of this disk, which we model with the radiative transfer code MCMax3D. We obtained J- and H-band observations in polarized scattered light with VLT/SPHERE that show strong azimuthal asymmetries. In the outer disk, beyond ~30 au, a broad shadow spans a range of ~240° in position angle, in the East. A bright ring at ~16 au also shows an azimuthally asymmetric brightness, with the faintest side roughly coincidental with the brightest region of the outer disk. Additionally, two arcs are detected at ~34 au and ~50 au. We created a simple 4-zone approximation to a warped disk model of HD 139614 in order to qualitatively reproduce these features. The location and misalignment of the disk components were constrained from the shape and location of the shadows they cast. We find that the shadow on the outer disk covers a range of position angle too wide to be explained by a single inner misaligned component. Our model requires a minimum of two separate misaligned zones -- or a continuously warped region -- to cast this broad shadow on the outer disk. A small misalignment of ~4° between adjacent components can reproduce most of the observed shadow features. Multiple misaligned disk zones, potentially mimicing a warp, can explain the observed broad shadows in the HD 139614 disk. A planetary mass companion in the disk, located on an inclined orbit, could be responsible for such a feature and for the dust depleted gap responsible for a dip in the SED.

preprint2019arXiv

The HI/OH/Recombination line survey of the inner Milky Way (THOR): data release 2 and HI overview

With the $Karl~G.~Jansky$ Very Large Array (VLA) in C-configuration, we observed a large portion of the first Galactic quadrant ($l=14.0-67.4^\circ $ and $\lvert b \rvert \leq 1.25^\circ $) achieving an angular resolution of $\leq 40^{\prime\prime}$. At $L$ Band, the WIDAR correlator at the VLA was set to cover the 21~cm HI line, four OH transitions, a series of H$nα$ radio recombination lines (RRLs; $n=151$ to 186), and eight 128~MHz wide continuum spectral windows (SPWs) simultaneously. The HI emission shows clear filamentary substructures at negative velocities with low velocity crowding. The emission at positive velocities is more smeared-out likely due to higher spatial and velocity crowding of structures at the positive velocities. Comparing to the spiral arm model of the Milky Way, the atomic gas follows the Sagittarius and Perseus Arm well but with significant material in the inter-arm regions. With the C-configuration-only HI+continuum data, we produced a HI optical depth map of the THOR areal coverage from 228 absorption spectra with the nearest-neighbor method. With this $τ$ map, we corrected the HI emission for optical depth and the derived column density is 38% higher than the column density with optically thin assumption. The total HI mass with optical depth correction in the survey region is 4.7$\times10^8~M_\odot$, 31% more than the mass derived assuming the emission is optically thin. If we apply this 31% correction to the whole Milky Way, the total atomic gas mass would be 9.4-10.5$\times 10^9~M_\odot$. Comparing the HI with existing CO data, we find a significant increase in the atomic-to-molecular gas ratio from the spiral arms to the inter-arm regions.

preprint2019arXiv

The polarimetric imaging mode of VLT/SPHERE/IRDIS I: Description, data reduction and observing strategy

Context. Polarimetric imaging is one of the most effective techniques for high-contrast imaging and characterization of protoplanetary disks, and has the potential to be instrumental in characterizing exoplanets. VLT/SPHERE contains the InfraRed Dual-band Imager and Spectrograph (IRDIS) with a dual-beam polarimetric imaging (DPI) mode, which offers the capability to obtain linear polarization images at high contrast and resolution. Aims. We aim to provide an overview of IRDIS/DPI and study its optical design to improve observing strategies and data reduction. Methods. For H-band observations of TW Hya, we compare two data reduction methods that correct for instrumental polarization effects in different ways: a minimization of the noise image, and a polarimetric-model-based correction method that we present in Paper II of this study. Results. We use observations of TW Hya to illustrate the data reduction. In the images of the protoplanetary disk around this star we detect variability in the polarized intensity and angle of linear polarization with pointing-dependent instrument configuration. We explain these variations as instrumental polarization effects and correct for these effects using our model-based correction method. Conclusions. IRDIS/DPI has proven to be a very successful and productive high-contrast polarimetric imaging system. However, the instrument performance depends on the specific instrument configuration. We suggest adjustments to future observing strategies to optimize polarimetric efficiency in field tracking mode by avoiding unfavourable derotator angles. We recommend reducing on-sky data with the pipeline called IRDAP that includes the model-based correction method (described in Paper II) to optimally account for the remaining telescope and instrumental polarization effects and to retrieve the true polarization state of the incident light.

preprint2019arXiv

The polarimetric imaging mode of VLT/SPHERE/IRDIS II: Characterization and correction of instrumental polarization effects

Context. Circumstellar disks and self-luminous giant exoplanets or companion brown dwarfs can be characterized through direct-imaging polarimetry at near-infrared wavelengths. SPHERE/IRDIS at the Very Large Telescope has the capabilities to perform such measurements, but uncalibrated instrumental polarization effects limit the attainable polarimetric accuracy. Aims. We aim to characterize and correct the instrumental polarization effects of the complete optical system, i.e. the telescope and SPHERE/IRDIS. Methods. We create a detailed Mueller matrix model in the broadband filters Y-, J-, H- and Ks, and calibrate it using measurements with SPHERE's internal light source and observations of two unpolarized stars. We develop a data-reduction method that uses the model to correct for the instrumental polarization effects, and apply it to observations of the circumstellar disk of T Cha. Results. The instrumental polarization is almost exclusively produced by the telescope and SPHERE's first mirror and varies with telescope altitude angle. The crosstalk primarily originates from the image derotator (K-mirror). At some orientations, the derotator causes severe loss of signal (>90% loss in H- and Ks-band) and strongly offsets the angle of linear polarization. With our correction method we reach in all filters a total polarimetric accuracy of <0.1% in the degree of linear polarization and an accuracy of a few degrees in angle of linear polarization. Conclusions. The correction method enables us to accurately measure the polarized intensity and angle of linear polarization of circumstellar disks, and is a vital tool for detecting unresolved (inner) disks and measuring the polarization of substellar companions. We have incorporated the correction method in a highly-automatic end-to-end data-reduction pipeline called IRDAP which is publicly available at https://irdap.readthedocs.io.

preprint2019arXiv

The young stellar content of the giant HII regions M8, G333.6-0.2, and NGC6357 with VLT/KMOS

Context: The identification and characterisation of populations of young massive stars in (giant) HII regions provides important constraints on i) the formation process of massive stars and their early feedback on the environment, and ii) the initial conditions for population synthesis models predicting the evolution of ensembles of stars. Aims: We identify and characterise the stellar populations of the following young giant HII regions: M8, G333.6-0.2, and NGC6357. Methods: We acquired H- and K-band spectra of around 200 stars using The K-band KMOS on the ESO Very Large Telescope. The targets for M8 and NGC6357 were selected from the MYStIX project, which combines X-ray observations with near-infrared and mid-infrared data. For G333.6-0.2, the sample selection is based on the near-infrared colours combined with X-ray data. We introduce an automatic spectral classification method in order to obtain temperatures and luminosities for the observed stars. We analyse the stellar populations using their photometric, astrometric, and spectroscopic properties and compared the position of the stars in the Hertzprung-Russell diagram with stellar evolution models to constrain their ages and mass ranges. Results: We confirm the presence of candidate ionising sources in the three regions and report new ones, including the first spectroscopically identified O stars in G333.6-0.2. In M8 and NGC6357, two populations are identified: i) OB main-sequence stars ($M > 5~\rm{M_{\odot}}$) and ii) pre-main sequence stars ($M\approx0.5-5~\rm{M_{\odot}}$). The ages of the clusters are $\sim$1-3~Myr, $< 3$~Myr, and $\sim$0.5-3~Myr for M8, G333.6-0.2, and NGC6357, respectively. We show that MYStIX selected targets have $>$ 90\% probability of being members of the HII region, whereas a selection based on near infrared (NIR) colours leads to a membership probability of only $\sim$70\%.

preprint2016arXiv

A $Herschel-SPIRE$ Survey of the Mon R2 Giant Molecular Cloud: Analysis of the Gas Column Density Probability Density Function

We present a far-IR survey of the entire Mon R2 GMC with $Herschel-SPIRE$ cross-calibrated with $Planck-HFI$ data. We fit the SEDs of each pixel with a greybody function and an optimal beta value of 1.8. We find that mid-range column densities obtained from far-IR dust emission and near-IR extinction are consistent. For the entire GMC, we find that the column density histogram, or N-PDF, is lognormal below $\sim$10$^{21}$ cm$^{-2}$. Above this value, the distribution takes a power law form with an index of -2.16. We analyze the gas geometry, N-PDF shape, and YSO content of a selection of subregions in the cloud. We find no regions with pure lognormal N-PDFs. The regions with a combination of lognormal and one power law N-PDF have a YSO cluster and a corresponding centrally concentrated gas clump. The regions with a combination of lognormal and two power law N-PDF have significant numbers of typically younger YSOs but no prominent YSO cluster. These regions are composed of an aggregate of closely spaced gas filaments with no concentrated dense gas clump. We find that for our fixed scale regions, the YSO count roughly correlates with the N-PDF power law index. The correlation appears steeper for single power law regions relative to two power law regions with a high column density cut-off, as a greater dense gas mass fraction is achieved in the former. A stronger correlation is found between embedded YSO count and the dense gas mass among our regions.

preprint2016arXiv

An M-dwarf star in the transition disk of Herbig HD 142527; Physical parameters and orbital elements

HD 142527A is one of the most studied Herbig Ae/Be stars with a transitional disk, as it has the largest imaged gap in any protoplanetary disk: the gas is cleared from 30 to 90 AU. The HD142527 system is also unique in that it has a stellar companion with a small mass compared to the mass of the primary star. This factor of ~20 in mass ratio between the two objects makes this binary system different from any other YSO. The HD142527 system could therefore provide a valuable test bed. This low-mass stellar object may be responsible for both the gap and dust trapping observed by ALMA at longer distances. We observed this system with the NACO and GPI instruments using the aperture masking technique. Aperture masking is ideal for providing high dynamic range even at very small angular separations. We present the spectral energy distribution for HD142527A and B. Brightness of the companion is now known from the R band up to the M' band. We also followed the orbital motion of HD 142527B over a period of more than two years. The SED of the companion is compatible with a T=3000+/-100K object in addition to a 1700K blackbody environment (likely a circus-secondary disk). From evolution models, we find that it is compatible with an object of mass 0.13+/-0.03Msun, radius 0.90+/-0.15Rsun, and age $1.0^{+1.0}_{-0.75}$Myr. This age is significantly younger than the age previously estimated for HD142527A. Computations to constrain the orbital parameters found a semi major axis of $140^{+120}_{-70}$mas, an eccentricity of 0.5+/-0.2, an inclination of 125+/-15 degrees, and a position angle of the right ascending node of -5+/-40 degrees. Inclination and position angle of the ascending node are in agreement with an orbit coplanar with the inner disk, not coplanar with the outer disk. Despite its high eccentricity, it is unlikely that HD142527B is responsible for truncating the inner edge of the outer disk.

preprint2016arXiv

An overview of the mid-infrared spectro-interferometer MATISSE: science, concept, and current status

MATISSE is the second-generation mid-infrared spectrograph and imager for the Very Large Telescope Interferometer (VLTI) at Paranal. This new interferometric instrument will allow significant advances by opening new avenues in various fundamental research fields: studying the planet-forming region of disks around young stellar objects, understanding the surface structures and mass loss phenomena affecting evolved stars, and probing the environments of black holes in active galactic nuclei. As a first breakthrough, MATISSE will enlarge the spectral domain of current optical interferometers by offering the L and M bands in addition to the N band. This will open a wide wavelength domain, ranging from 2.8 to 13 um, exploring angular scales as small as 3 mas (L band) / 10 mas (N band). As a second breakthrough, MATISSE will allow mid-infrared imaging - closure-phase aperture-synthesis imaging - with up to four Unit Telescopes (UT) or Auxiliary Telescopes (AT) of the VLTI. Moreover, MATISSE will offer a spectral resolution range from R ~ 30 to R ~ 5000. Here, we present one of the main science objectives, the study of protoplanetary disks, that has driven the instrument design and motivated several VLTI upgrades (GRA4MAT and NAOMI). We introduce the physical concept of MATISSE including a description of the signal on the detectors and an evaluation of the expected performances. We also discuss the current status of the MATISSE instrument, which is entering its testing phase, and the foreseen schedule for the next two years that will lead to the first light at Paranal.

preprint2016arXiv

Continuum sources from the THOR survey between 1 and 2 GHz

We carried out a large program with the Karl G. Jansky Very Large Array (VLA): "THOR: The HI, OH, Recombination line survey of the Milky Way". We observed a significant portion of the Galactic plane in the first quadrant of the Milky Way in the 21cm HI line, 4 OH transitions, 19 radio recombination lines, and continuum from 1 to 2 GHz. In this paper we present a catalog of the continuum sources in the first half of the survey (l=14.0-37.9deg and l=47.1-51.2deg, |b|<1.1deg) at a spatial resolution of 10-25", with a spatially varying noise level of ~0.3-1 mJy/beam. The catalog contains ~4400 sources. Around 1200 of these are spatially resolved, and ~1000 are possible artifacts, given their low signal-to-noise ratios. Since the spatial distribution of the unresolved objects is evenly distributed and not confined to the Galactic plane, most of them are extragalactic. Thanks to the broad bandwidth of the observations from 1 to 2 GHz, we are able to determine a reliable spectral index for ~1800 sources. The spectral index distribution reveals a double-peaked profile with maxima at spectral indices of alpha = -1 and alpha = 0 , corresponding to steep declining and flat spectra, respectively. This allows us to distinguish between thermal and non-thermal emission, which can be used to determine the nature of each source. We examine the spectral index of ~300 known HII regions, for which we find thermal emission with spectral indices around alpha = 0. In contrast, supernova remnants (SNR) show non-thermal emission with alpha = -0.5 and extragalactic objects generally have a steeper spectral index of alpha = -1. Using the spectral index information of the THOR survey, we investigate potential SNR candidates. We classify the radiation of four SNR candidates as non-thermal, and for the first time, we provide strong evidence for the SNR origin of these candidates.

preprint2016arXiv

Hard X-ray irradiation of cosmic silicate analogs: structural evolution and astrophysical implications

Protoplanetary disks, interstellar clouds, and active galactic nuclei, contain X-ray dominated regions. X-rays interact with the dust and gas present in such environments. While a few laboratory X-ray irradiation experiments have been performed on ices, X-ray irradiation experiments on bare cosmic dust analogs have been scarce up to now. Our goal is to study the effects of hard X-rays on cosmic dust analogs via in-situ X-ray diffraction. By using a hard X-ray synchrotron nanobeam, we seek to simulate cumulative X-ray exposure on dust grains during their lifetime in these astrophysical environments, and provide an upper limit on the effect of hard X-rays on dust grain structure. We prepared enstatite nanograins, analogs to cosmic silicates, via the melting-quenching technique. These amorphous grains were then annealed to obtain polycrystalline grains. These were characterized via scanning electron microscopy and high-resolution transmission electron microscopy before irradiation. Powder samples were prepared in X-ray transparent substrates and were irradiated with hard X-rays nanobeams (29.4 keV) provided by beamline ID16B of the European Synchrotron Radiation Facility. X-ray diffraction images were recorded in transmission mode. We detected the amorphization of polycrystalline silicates embedded in an organic matrix after an accumulated X-ray exposure of 6.4 x 10$^{27}$ eV cm$^{-2}$. Pure crystalline silicate grains (without resin) did not exhibit amorphization. None of the amorphous silicate samples (pure and embedded in resin) underwent crystallization. We analyzed the evolution of the polycrystalline sample embedded in an organic matrix as a function of X-ray exposure. Loss of diffraction peak intensity, peak broadening, and the disappearance of discrete spots and arcs, revealed the amorphization of the resin embedded (originally polycrystalline) silicate sample.

preprint2016arXiv

HATS-11b and HATS-12b: Two transiting Hot Jupiters orbiting sub-solar metallicity stars selected for the K2 Campaign 7

We report the discovery of two transiting extrasolar planets from the HATSouth survey. HATS-11, a V=14.1 G0-star shows a periodic 12.9 mmag dip in its light curve every 3.6192 days and a radial velocity variation consistent with a Keplerian orbit. HATS-11 has a mass of 1.000 $\pm$ 0.060 M$_{\odot}$, a radius of 1.444 $\pm$ 0.057 M$_{\odot}$ and an effective temperature of 6060 $\pm$ 150 K, while its companion is a 0.85 $\pm$ 0.12 M$_J$, 1.510 $\pm$ 0.078 R$_J$ planet in a circular orbit. HATS-12 shows a periodic 5.1 mmag flux decrease every 3.1428 days and Keplerian RV variations around a V=12.8 F-star. HATS-12 has a mass of 1.489 $\pm$ 0.071 M$_{\odot}$, a radius of 2.21 $\pm$ 0.21 R$_{\odot}$, and an effective temperature of 6408 $\pm$ 75 K. For HATS-12, our measurements indicate that this is a 2.38 $\pm$ 0.11 M$_J$, 1.35 $\pm$ 0.17 R$_J$ planet in a circular orbit. Both host stars show sub-solar metallicity of -0.390 $\pm$ 0.060 dex and -0.100 $\pm$ 0.040 dex, respectively and are (slightly) evolved stars. In fact, HATS-11 is amongst the most metal-poor and, HATS-12 is amongst the most evolved stars hosting a hot Jupiter planet. Importantly, HATS-11 and HATS-12 have been observed in long cadence by Kepler as part of K2 campaign 7 (EPIC216414930 and EPIC218131080 respectively).

preprint2016arXiv

HATS-19b, HATS-20b, HATS-21b: Three Transiting Hot-Saturns Discovered by the HATSouth Survey

We report the discovery by the HATSouth exoplanet survey of three hot-Saturn transiting exoplanets: HATS-19b, HATS-20b, and HATS-21b. The planet host HATS-19 is a slightly evolved V = 13.0 G0 star with [Fe/H] = 0.240, a mass of 1.303 Msun, and a radius of 1.75 Rsun. HATS-19b is in an eccentric orbit (e = 0.30) around this star with an orbital period of 4.5697 days and has a mass of 0.427 Mjup and a highly inflated radius of 1.66 Rjup. The planet HATS-20b has a Saturn-like mass and radius of 0.273 Mjup and 0.776 Rjup respectively. It orbits the V = 13.8 G9V star HATS-20 (Ms = 0.910 Msun; Rs = 0.892 Rsun) with a period of 3.7993 days. Finally, HATS-21 is a V = 12.2 G4V star with [Fe/H] = 0.300, a mass of 1.080 Msun, and a radius of 1.021 Rsun. Its accompanying planet HATS-21b has a 3.5544-day orbital period, a mass of 0.332 Mjup, and a moderately inflated radius of 1.123 Rjup. With the addition of these three very different planets to the growing sample of hot-Saturns, we re-examine the relations between the observed giant planet radii, stellar irradiation, and host metallicity. We find a significant positive correlation between planet equilibrium temperature and radius, and a weak negative correlation between host metallicity and radius. To assess the relative influence of various physical parameters on observed planet radii, we train and fit models using Random Forest regression. We find that for hot-Saturns (0.1 < Mp < 0.5 Mjup), the planetary mass and equilibrium temperature play dominant roles in determining radii. For hot-Jupiters (0.5 < Mp < 2.0 Mjup), the most important parameter is equilibrium temperature alone. Finally, for irradiated higher-mass planets (Mp > 2.0 Mjup), we find that equilibrium temperature dominates in influence, with smaller contributions from planet mass and host metallicity.

preprint2016arXiv

HATS-25b through HATS-30b: A Half-dozen New Inflated Transiting Hot Jupiters from the HATSouth Survey

We report six new inflated hot Jupiters (HATS-25b through HATS-30b) discovered using the HATSouth global network of automated telescopes. The planets orbit stars with $V$ magnitudes in the range $\sim 12-14$ and have masses in the largely populated $0.5M_J-0.7M_J$ region of parameter space but span a wide variety of radii, from $1.17R_J$ to $1.75 R_J$. HATS-25b, HATS-28b, HATS-29b and HATS-30b are typical inflated hot Jupiters ($R_p = 1.17-1.26R_J$) orbiting G-type stars in short period ($P=3.2-4.6$ days) orbits. However, HATS-26b ($R_p = 1.75R_J$, $P = 3.3024$ days) and HATS-27b ($R_p=1.50R_J$, $P=4.6370$ days) stand out as highly inflated planets orbiting slightly evolved F stars just after and in the turn-off points, respectively, which are among the least dense hot Jupiters, with densities of $0.153$ g cm$^{-3}$ and $0.180$ g cm$^{-3}$, respectively. All the presented exoplanets but HATS-27b are good targets for future atmospheric characterization studies, while HATS-27b is a prime target for Rossiter-McLaughlin monitoring in order to determine its spin-orbit alignment given the brightness ($V = 12.8$) and stellar rotational velocity ($v \sin i \approx 9.3$ km/s) of the host star. These discoveries significantly increase the number of inflated hot Jupiters known, contributing to our understanding of the mechanism(s) responsible for hot Jupiter inflation.

preprint2016arXiv

HATS-31b Through HATS-35b: Five Transiting Hot Jupiters Discovered by the HATSouth Survey

We report the discovery of five new transiting hot Jupiter planets discovered by the HATSouth survey: HATS-31b through HATS-35b. These planets orbit moderately bright stars with V magnitudes within the range 11.9-14.4mag while the planets span a range of masses 0.88-1.22MJ, and have somewhat inflated radii between 1.23-1.64RJ.These planets can be classified as typical hot Jupiters, with HATS-31b and HATS-35b being moderately inflated gas giant planets with radii of $1.64 \pm 0.22$ RJ and 1.464+0.069-0.044RJ, respectively, that can be used to constrain inflation mechanisms. All five systems present a higher Bayesian evidence for a fixed circular orbit model than for an eccentric orbit. The orbital periods range from $1.8209993 \pm 0.0000016$ day for HATS-35b) to $3.377960 \pm 0.000012$ day for HATS-31b. Additionally, HATS-35b orbits a relatively young F star with an age of $2.13 \pm 0.51$ Gyr. We discuss the analysis to derive the properties of these systems and compare them in the context of the sample of well characterized transiting hot Jupiters known to date.

preprint2016arXiv

New transit observations for HAT-P-30 b, HAT-P-37 b, TrES-5 b, WASP-28 b, WASP-36 b, and WASP-39 b

We present new transit light curves for planets in six extrasolar planetary systems. They were acquired with 0.4-2.2 m telescopes located in west Asia, Europe, and South America. When combined with literature data, they allowed us to redetermine system parameters in a homogeneous way. Our results for individual systems are in agreement with values reported in previous studies. We refined transit ephemerides and reduced uncertainties of orbital periods by a factor between 2 and 7. No sign of any variations in transit times was detected for the planets studied.

preprint2016arXiv

Physical parameters of late M-type members of Chamaleon I and TW Hydrae Association: Dust settling, age dispersion and activity

Although mid-to-late type M dwarfs are the most common stars in our stellar neighborhood, our knowledge of these objects is still limited. Open questions include the evolution of their angular momentum, internal structures, dust settling in their atmospheres, age dispersion within populations. In addition, at young ages, late-type Ms have masses below the hydrogen burning limit and therefore are key objects in the debate on the brown dwarf mechanism of formation. In this work we determine and study in detail the physical parameters of two samples of young, late M-type sources belonging to either the Chamaeleon I Dark Cloud or the TW Hydrae Association and compare them with the results obtained in the literature for other young clusters and also for older, field, dwarfs. We used multi-wavelength photometry to construct and analyze SEDs to determine general properties of the photosphere and disk presence. We also used low resolution optical and near-infrared spectroscopy to study activity, accretion, gravity and effective temperature sensitive indicators. We propose a VO-based spectral index that is both temperature and age sensitive. We derived physical parameters using independent techniques confirming the already common feature/problem of the age/luminosity spread. In particular, we highlight two brown dwarfs showing very similar temperatures but clearly different surface gravity (explained invoking extreme early accretion). We also show how, despite large improvement in the dust treatment in theoretical models, there is still room for further progress in the simultaneous reproduction of the optical and near-infrared features of these cold young objects.

preprint2016arXiv

Polarimetry and flux distribution in the debris disk around HD 32297

We present high-contrast angular differential imaging (ADI) observations of the debris disk around HD 32297 in H-band, as well as the first polarimetric images for this system in polarized differential imaging (PDI) mode with Subaru/HICIAO. In ADI, we detect the nearly edge-on disk at >5sigma levels from ~0.45 arcsec to ~1.7 arcsec (50-192 AU) from the star and recover the spine deviation from the midplane already found in previous works. We also find for the first time imaging and surface brightness (SB) indications for the presence of a gapped structure on both sides of the disk at distances of ~0.75 arcsec (NE side) and ~0.65 arcsec (SW side). Global forward-modelling work delivers a best-fit model disk and well-fitting parameter intervals that essentially match previous results, with high-forward scattering grains and a ring located at 110 AU. However, this single ring model cannot account for the gapped structure seen in our SB profiles. We create simple double ring models and achieve a satisfactory fit with two rings located at 60 and 95 AU, respectively, low-forward scattering grains and very sharp inner slopes. In polarized light we retrieve the disk extending from ~0.25-1.6 arcsec, although the central region is quite noisy and high S/N are only found in the range ~0.75-1.2 arcsec. The disk is polarized in the azimuthal direction, as expected, and the departure from the midplane is also clearly observed. Evidence for a gapped scenario is not found in the PDI data. We obtain a linear polarization degree of the grains that increases from ~10% at 0.55 arcsec to ~25% at 1.6 arcsec. The maximum is found at scattering angles of ~90degrees, either from the main components of the disk or from dust grains blown out to larger radii.

preprint2016arXiv

Resolving the chemical substructure of Orion-KL

The Kleinmann-Low nebula in Orion (Orion-KL) is the nearest example of a high-mass star-forming environment. For the first time, we complemented 1.3 mm Submillimeter Array (SMA) interferometric line survey with IRAM 30 m single-dish observations of the Orion-KL region. Covering a 4 GHz bandwidth in total, this survey contains over 160 emission lines from 20 species (25 isotopologues), including 11 complex organic molecules (COMs). At a spatial resolution of 1200 AU, the continuum substructures are resolved. Extracting the spectra from individual substructures and providing the intensity-integrated distribution map for each species, we studied the small-scale chemical variations in this region. Our main results are: (1) We identify lines from the low-abundance COMs CH3COCH3 and CH3CH2OH, as well as tentatively detect CH3CHO and long carbon-chains C6H and HC7N. (2) We find that while most COMs are segregated by type, peaking either towards the hot core (e.g., N-bearing species) or the compact ridge (e.g., O-bearing species like HCOOCH3 and CH3OCH3), while the distributions of others do not follow this segregated structure (e.g., CH3CH2OH, CH3OH, CH3COCH3). (3) We find a second velocity component of HNCO, SO2, 34SO2, and SO lines, which may be associated with a strong shock event in the low-velocity outflow. (4) Temperatures and molecular abundances show large gradients between central condensations and the outflow regions, illustrating a transition between hot molecular core and shock-chemistry dominated regimes. Our observations of spatially resolved chemical variations in Orion-KL provide the nearest reference source for hot molecular core and outflow chemistry, which will be an important example for interpreting the chemistry of more distant HMSFRs.

preprint2016arXiv

Resolving the planet-hosting inner regions of the LkCa 15 disk

LkCa 15 hosts a pre-transitional disk as well as at least one accreting protoplanet orbiting in its gap. Previous disk observations have focused mainly on the outer disk, which is cleared inward of ~50 au. The planet candidates, on the other hand, reside at orbital radii around 15 au, where disk observations have been unreliable until recently. Here we present new J-band imaging polarimetry of LkCa 15 with SPHERE IRDIS, yielding the most accurate and detailed scattered-light images of the disk to date down to the planet-hosting inner regions. We find what appear to be persistent asymmetric structures in the scattering material at the location of the planet candidates, which could be responsible at least for parts of the signals measured with sparse-aperture masking. These images further allow us to trace the gap edge in scattered light at all position angles and search the inner and outer disks for morphological substructure. The outer disk appears smooth with slight azimuthal variations in polarized surface brightness, which may be due to shadowing from the inner disk or a two-peaked polarized phase function. We find that the near-side gap edge revealed by polarimetry matches the sharp crescent seen in previous ADI imaging very well. Finally, the ratio of polarized disk to stellar flux is more than six times larger in J-band than in the RI bands.

preprint2016arXiv

SEEDS direct imaging of the RV-detected companion to V450 Andromedae, and characterization of the system

We report the direct imaging detection of a low-mass companion to a young, moderately active star V450 And, that was previously identified with the radial velocity method. The companion was found in high-contrast images obtained with the Subaru Telescope equipped with the HiCIAO camera and AO188 adaptive optics system. From the public ELODIE and SOPHIE archives we extracted available high-resolution spectra and radial velocity (RV) measurements, along with RVs from the Lick planet search program. We combined our multi-epoch astrometry with these archival, partially unpublished RVs, and found that the companion is a low-mass star, not a brown dwarf, as previously suggested. We found the best-fitting dynamical masses to be $m_1=1.141_{-0.091}^{+0.037}$ and $m_2=0.279^{+0.023}_{-0.020}$ M$_\odot$. We also performed spectral analysis of the SOPHIE spectra with the iSpec code. The Hipparcos time-series photometry shows a periodicity of $P=5.743$ d, which is also seen in SOPHIE spectra as an RV modulation of the star A. We interpret it as being caused by spots on the stellar surface, and the star to be rotating with the given period. From the rotation and level of activity, we found that the system is $380^{+220}_{-100}$ Myr old, consistent with an isochrone analysis ($220^{+2120}_{-90}$ Myr). This work may serve as a test case for future studies of low-mass stars, brown dwarfs and exoplanets by combination of RV and direct imaging data.

preprint2016arXiv

Shadows and spirals in the protoplanetary disk HD 100453

Understanding the diversity of planets requires to study the morphology and the physical conditions in the protoplanetary disks in which they form. We observed and spatially resolved the disk around the ~10 Myr old protoplanetary disk HD 100453 in polarized scattered light with SPHERE/VLT at optical and near-infrared wavelengths, reaching an angular resolution of ~0.02", and an inner working angle of ~0.09". We detect polarized scattered light up to ~0.42" (~48 au) and detect a cavity, a rim with azimuthal brightness variations at an inclination of 38 degrees, two shadows and two symmetric spiral arms. The spiral arms originate near the location of the shadows, close to the semi major axis. We detect a faint spiral-like feature in the SW that can be interpreted as the scattering surface of the bottom side of the disk, if the disk is tidally truncated by the M-dwarf companion currently seen at a projected distance of ~119 au. We construct a radiative transfer model that accounts for the main characteristics of the features with an inner and outer disk misaligned by ~72 degrees. The azimuthal brightness variations along the rim are well reproduced with the scattering phase function of the model. While spirals can be triggered by the tidal interaction with the companion, the close proximity of the spirals to the shadows suggests that the shadows could also play a role. The change in stellar illumination along the rim, induces an azimuthal variation of the scale height that can contribute to the brightness variations. Dark regions in polarized images of transition disks are now detected in a handful of disks and often interpreted as shadows due to a misaligned inner disk. The origin of such a misalignment in HD 100453, and of the spirals, is unclear, and might be due to a yet-undetected massive companion inside the cavity, and on an inclined orbit.

preprint2016arXiv

Spiral Structure and Differential Dust Size Distribution in the LkHa 330 Disk

Dust trapping accelerates the coagulation of dust particles, and thus it represents an initial step toward the formation of planetesimals. We report $H$-band (1.6 um) linear polarimetric observations and 0.87 mm interferometric continuum observations toward a transitional disk around LkHa 330. As results, a pair of spiral arms were detected in the $H$-band emission and an asymmetric (potentially arm-like) structure was detected in the 0.87 mm continuum emission. We discuss the origin of the spiral arm and the asymmetric structure, and suggest that a massive unseen planet is the most plausible explanation. The possibility of dust trapping and grain growth causing the asymmetric structure was also investigated through the opacity index (beta) by plotting the observed SED slope between 0.87 mm from our SMA observation and 1.3 mm from literature. The results imply that grains are indistinguishable from ISM-like dust in the east side ($beta = 2.0 pm 0.5$), but much smaller in the west side $beta = 0.7^{+0.5}_{-0.4}$, indicating differential dust size distribution between the two sides of the disk. Combining the results of near-infrared and submillimeter observations, we conjecture that the spiral arms exist at the upper surface and an asymmetric structure resides in the disk interior. Future observations at centimeter wavelengths and differential polarization imaging in other bands (Y to K) with extreme AO imagers are required to understand how large dust grains form and to further explore the dust distribution in the disk.

preprint2016arXiv

The evolution of far-infrared CO emission from protostars

We investigate the evolution of far-IR CO emission from protostars observed with Herschel/PACS for 50 sources from the combined sample of HOPS and DIGIT Herschel key programs. From the uniformly sampled spectral energy distributions, we computed $L_{\rm{bol}}$, $T_{\rm{bol}}$ and $L_{\rm {bol}}/L_{\rm {smm}}$ for these sources to search for correlations between far-IR CO emission and protostellar properties. We find a strong and tight correlation between far-IR CO luminosity ($L^{\rm fir}_{\rm CO}$) and the bolometric luminosity ($L_{\rm{bol}}$) of the protostars with $L^{\rm fir}_{\rm CO}$ $\propto$ $L_{\rm{bol}}^{0.7}$. We, however, do not find a strong correlation between $L^{\rm fir}_{\rm CO}$ and protostellar evolutionary indicators, $T_{\rm{bol}}$ and $L_{\rm {bol}}/L_{\rm {smm}}$. FIR CO emission from protostars traces the currently shocked gas by jets/outflows, and $L^{\rm fir}_{\rm CO}$ is proportional to the instantaneous mass loss rate, $\dot{M}_{\rm{out}}$. The correlation between $L^{\rm fir}_{\rm CO}$ and $L_{\rm{bol}}$ is indicative of instantaneous $\dot{M}_{\rm{out}}$ tracking instantaneous $\dot{M}_{\rm{acc}}$. The lack of correlation between $L^{\rm fir}_{\rm CO}$ and evolutionary indicators $T_{\rm{bol}}$ and $L_{\rm {bol}}/L_{\rm {smm}}$ suggests that $\dot{M}_{\rm{out}}$ and, therefore, $\dot{M}_{\rm{acc}}$ do not show any clear evolutionary trend. These results are consistent with mass accretion/ejection in protostars being episodic. Taken together with the previous finding that the time-averaged mass ejection/accretion rate declines during the protostellar phase (e.g., Bontemps et al. 1996), our results suggest that the instantaneous accretion/ejection rate of protostars is highly time variable and episodic, but the amplitude and/or frequency of this variability decreases with time such that the time averaged accretion/ejection rate declines with system age.

preprint2016arXiv

The shadow of the Flying Saucer: A very low temperature for large dust grains

Dust determines the temperature structure of protoplanetary disks. However, dust temperature determinations almost invariably rely on a complex modeling of the Spectral Energy Distribution. We attempt a direct determination of the temperature of large grains emitting at mm wavelengths.} We observe the edge-on dust disk of the Flying Saucer, which appears in silhouette against the CO J=2-1 emission from a background molecular cloud in $ρ$ Oph. The combination of velocity gradients due to the Keplerian rotation of the disk and intensity variations in the CO background as a function of velocity allows us to directly measure the %absorbing dust temperature. The dust opacity can then be derived from the emitted continuum radiation. The dust disk absorbs the radiation from the CO clouds at several velocities. We derive very low dust temperatures, 5 to 7 K at radii around 100 au, which is much lower than most model predictions. The dust optical depth is $> 0.2$ at 230 GHz, and the scale height at 100 au is at least 8 au (best fit 13 au). However, the dust disk is very flat (flaring index -0.35), which is indicative of dust settling in the outer parts.

preprint2015arXiv

A 0.24+0.18 Msun double-lined eclipsing binary from the HATSouth survey

We report the discovery and characterisation of a new M-dwarf binary, with component masses and radii of M1 = 0.244 -0.003/+0.003 Msun, R1 = 0.261 -0.009/+0.006 Rsun, M2 = 0.179 -0.001/+0.002 Msun, R2 = 0.218 -0.011/+0.007 Rsun, and orbital period of ~4.1 days. The M-dwarf binary HATS551-027 (LP 837-20) was identified as an eclipsing binary by the HATSouth survey, and characterised by a series of high precision photometric observations of the eclipse events, and spectroscopic determinations of the atmospheric parameters and radial velocity orbits. HATS551-027 is one of few systems with both stellar components lying in the fully-convective regime of very low mass stars, and can serve as a test for stellar interior models. The radius of HATS551-027A is consistent with models to 1 sigma, whilst HATS551-027B is inflated by 9% at 2 sigma significance. We measure the effective temperatures for the two stellar components to be Teff,1 = 3190 +/- 100 K and Teff,2 = 2990+/-110 K, both are slightly cooler than theoretical models predict, but consistent with other M-dwarfs of similar masses that have previously been studied. We also measure significant Halpha emission from both components of the binary system, and discuss this in the context of the correlation between stellar activity and the discrepancies between the observed and model temperatures.

preprint2015arXiv

A multi-wavelength analysis for interferometric (sub-)mm observations of protoplanetary disks: radial constraints on the dust properties and the disk structure

Theoretical models of grain growth predict dust properties to change as a function of protoplanetary disk radius, mass, age and other physical conditions. We lay down the methodology for a multi-wavelength analysis of (sub-)mm and cm continuum interferometric observations to constrain self-consistently the disk structure and the radial variation of the dust properties. The computational architecture is massively parallel and highly modular. The analysis is based on the simultaneous fit in the uv-plane of observations at several wavelengths with a model for the disk thermal emission and for the dust opacity. The observed flux density at the different wavelengths is fitted by posing constraints on the disk structure and on the radial variation of the grain size distribution. We apply the analysis to observations of three protoplanetary disks (AS 209, FT Tau, DR Tau) for which a combination of spatially resolved observations in the range ~0.88mm to ~10mm is available (from SMA, CARMA, and VLA), finding evidence of a decreasing maximum dust grain size (a_max) with radius. We derive large a_max values up to 1 cm in the inner disk between 15 and 30 AU and smaller grains with a_max~1 mm in the outer disk (R > 80AU). In this paper we develop a multi-wavelength analysis that will allow this missing quantity to be constrained for statistically relevant samples of disks and to investigate possible correlations with disk or stellar parameters.

preprint2015arXiv

A narrow, edge-on disk resolved around HD 106906 with SPHERE

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.

preprint2015arXiv

Asymmetric features in the protoplanetary disk MWC758

The study of dynamical processes in protoplanetary disks is essential to understand planet formation. In this context, transition disks are prime targets because they are at an advanced stage of disk clearing and may harbor direct signatures of disk evolution. In this paper, we aim to derive new constraints on the structure of the transition disk MWC 758, to detect non-axisymmetric features and understand their origin. We obtained infrared polarized intensity observations of the protoplanetary disk MWC 758 with SPHERE/VLT at 1.04 microns to resolve scattered light at a smaller inner working angle (0.093") and a higher angular resolution (0.027") than previously achieved. We observe polarized scattered light within 0.53" (148 au) down to the inner working angle (26 au) and detect distinct non-axisymmetric features but no fully depleted cavity. The two small-scale spiral features that were previously detected with HiCIAO are resolved more clearly, and new features are identified, including two that are located at previously inaccessible radii close to the star. We present a model based on the spiral density wave theory with two planetary companions in circular orbits. The best model requires a high disk aspect ratio (H/r~0.20 at the planet locations) to account for the large pitch angles which implies a very warm disk. Our observations reveal the complex morphology of the disk MWC758. To understand the origin of the detected features, the combination of high-resolution observations in the submillimeter with ALMA and detailed modeling is needed.

preprint2015arXiv

Discovery of a Disk Gap Candidate at 20 AU in TW Hydrae

We present a new Subaru/HiCIAO high-contrast H-band polarized intensity (PI) image of a nearby transitional disk associated with TW Hydrae. The scattered light from the disk was detected from 0.2" to 1.5" (11 - 81 AU) and the PI image shows a clear axisymmetric depression in polarized intensity at ~ 0.4" (~ 20 AU) from the central star, similar to the ~ 80 AU gap previously reported from HST images. Azimuthal polarized intensity profile also shows the disk beyond 0.2" is almost axisymmetric. We discuss two possible scenarios explaining the origin of the polarized intensity depression: 1) a gap structure may exist at ~ 20 AU from the central star because of shallow slope seen in the polarized intensity profile, and 2) grain growth may be occurring in the inner region of the disk. Multi-band observations at NIR and millimeter/sub-millimeter wavelengths play a complementary role in investigating dust opacity and may help reveal the origin of the gap more precisely.

preprint2015arXiv

First light of the VLT planet finder SPHERE. I. Detection and characterization of the sub-stellar companion GJ 758 B

GJ758 B is a brown dwarf companion to a nearby (15.76 pc) solar-type, metal-rich (M/H = +0.2 dex) main-sequence star (G9V) that was discovered with Subaru/HiCIAO in 2009. From previous studies, it has drawn attention as being the coldest (~600K) companion ever directly imaged around a neighboring star. We present new high-contrast data obtained during the commissioning of the SPHERE instrument at the VLT. The data was obtained in Y-, J-, H-, and Ks-bands with the dual-band imaging (DBI) mode of IRDIS, providing a broad coverage of the full near-infrared (near-IR) range at higher contrast and better spectral sampling than previously reported. In this new set of high-quality data, we report the re-detection of the companion, as well as the first detection of a new candidate closer-in to the star. We use the new 8 photometric points for an extended comparison of GJ758 B with empirical objects and 4 families of atmospheric models. From comparison to empirical object, we estimate a T8 spectral type, but none of the comparison object can accurately represent the observed near-IR fluxes of GJ758 B. From comparison to atmospheric models, we attribute a Teff = 600K $\pm$ 100K, but we find that no atmospheric model can adequately fit all the fluxes of GJ758 B. The photometry of the new candidate companion is broadly consistent with L-type objects, but a second epoch with improved photometry is necessary to clarify its status. The new astrometry of GJ758 B shows a significant proper motion since the last epoch. We use this result to improve the determination of the orbital characteristics using two fitting approaches, Least-Square Monte Carlo and Markov Chain Monte Carlo. Finally, we analyze the sensitivity of our data to additional closer-in companions and reject the possibility of other massive brown dwarf companions down to 4-5 AU. [abridged]

preprint2015arXiv

First light of the VLT planet finder SPHERE. II. The physical properties and the architecture of the young systems PZ Tel and HD 1160 revisited

[Abridged] Context. The young systems PZ Tel and HD 1160, hosting known low-mass companions, were observed during the commissioning of the new planet finder SPHERE with several imaging and spectroscopic modes. Aims. We aim to refine the physical properties and architecture of both systems. Methods. We use SPHERE commissioning data and REM observations, as well as literature and unpublished data from VLT/SINFONI, VLT/NaCo, Gemini/NICI, and Keck/NIRC2. Results. We derive new photometry and confirm the nearly daily photometric variability of PZ Tel A. Using literature data spanning 38 yr, we show that the star also exhibits a long-term variability trend. The 0.63-3.8 mic SED of PZ Tel B allows us to revise its properties: spectral type M7+/-1, Teff=2700+/-100 K, log(g)<4.5 dex, log(L/L_Sun)=-2.51+/-0.10 dex, and mass 38-72 MJ. The 1-3.8 mic SED of HD 1160 B suggests a massive brown dwarf or a low-mass star with spectral type M5.5-7.0, Teff=3000+/-100 K, [M/H]=-0.5-0.0 dex, log(L/L_Sun)=-2.81+/-0.10 dex, and mass 39-168 MJ. We confirm the deceleration and high eccentricity (e>0.66) of PZ Tel B. For e<0.9, the inclination, longitude of the ascending node, and time of periastron passage are well constrained. The system is seen close to an edge-on geometry. We reject other brown dwarf candidates outside 0.25" for both systems, and massive giant planets (>4 MJ) outside 0.5" for the PZ Tel system. We also show that K1-K2 color can be used with YJH low-resolution spectra to identify young L-type companions, provided high photometric accuracy (<0.05 mag) is achieved. Conclusions. SPHERE opens new horizons in the study of young brown dwarfs and giant exoplanets thanks to high-contrast imaging capabilities at optical and near-infrared wavelengths, as well as high signal-to-noise spectroscopy in the near-infrared from low (R~30-50) to medium resolutions (R~350).

preprint2015arXiv

First light of the VLT planet finder SPHERE. IV. Physical and chemical properties of the planets around HR8799

The system of four planets around HR8799 offers a unique opportunity to probe the physics and chemistry at play in the atmospheres of self-luminous young (~30 Myr) planets. We recently obtained new photometry of the four planets and low-resolution (R~30) spectra of HR8799 d and e with the SPHERE instrument (paper III). In this paper (paper IV), we compare the available spectra and photometry of the planets to known objects and atmospheric models (BT-SETTL14, Cloud-AE60, Exo-REM) to characterize the atmospheric properties of the planets. We find that HR8799d and e properties are well reproduced by those of L6-L8 dusty dwarfs discovered in the field, among which some are candidate members of young nearby associations. No known object reproduces well the properties of planets b and c. Nevertheless, we find that the spectra and WISE photometry of peculiar and/or young early-T dwarfs reddened by submicron grains made of corundum, iron, enstatite, or forsterite successfully reproduce the SED of these two planets. Our analysis confirms that only the Exo-REM models with thick clouds fit (within 2σ) the whole set of spectrophotometric datapoints available for HR8799 d and e for Teff = 1200 K, log g in the range 3.0-4.5, and M/H=+0.5. The models still fail to reproduce the SED of HR8799c and b. The determination of the metallicity, log g, and cloud thickness are degenerate. We conclude that an enhanced content in dust and decreased CIA of H2 is certainly responsible for the deviation of the properties of the planet with respect to field dwarfs. The analysis suggests in addition that HR8799c and b have later spectral types than the two other planets, and therefore could both have lower masses.

preprint2015arXiv

HATS-15 b and HATS-16 b: Two massive planets transiting old G dwarf stars

We report the discovery of HATS-15 b and HATS-16 b, two massive transiting extrasolar planets orbiting evolved ($\sim 10$ Gyr) main-sequence stars. The planet HATS-15 b, which is hosted by a G9V star ($V=14.8$ mag), is a hot Jupiter with mass of $2.17\pm0.15\, M_{\mathrm{J}}$ and radius of $1.105\pm0.0.040\, R_{\mathrm{J}}$, and completes its orbit in nearly 1.7 days. HATS-16 b is a very massive hot Jupiter with mass of $3.27\pm0.19\, M_{\mathrm{J}}$ and radius of $1.30\pm0.15\, R_{\mathrm{J}}$; it orbits around its G3 V parent star ($V=13.8$ mag) in $\sim2.7$ days. HATS-16 is slightly active and shows a periodic photometric modulation, implying a rotational period of 12 days which is unexpectedly short given its isochronal age. This fast rotation might be the result of the tidal interaction between the star and its planet.

preprint2015arXiv

HATS-17b: A Transiting Compact Warm Jupiter in a 16.3 Days Circular Orbit

We report the discovery of HATS-17b, the first transiting warm Jupiter of the HATSouth network. HATS-17b transits its bright (V=12.4) G-type (M$_{\star}$=1.131 $\pm$ 0.030 M$_{\odot}$, R$_{\star}$=1.091$^{+0.070}_{-0.046}$ R$_{\star}$) metal-rich ([Fe/H]=+0.3 dex) host star in a circular orbit with a period of P=16.2546 days. HATS-17b has a very compact radius of 0.777 $\pm$ 0.056 R$_J$ given its Jupiter-like mass of 1.338 $\pm$ 0.065 M$_J$. Up to 50% of the mass of HATS-17b may be composed of heavy elements in order to explain its high density with current models of planetary structure. HATS-17b is the longest period transiting planet discovered to date by a ground-based photometric survey, and is one of the brightest transiting warm Jupiter systems known. The brightness of HATS-17b will allow detailed follow-up observations to characterize the orbital geometry of the system and the atmosphere of the planet.

preprint2015arXiv

HATS-7b: A Hot Super Neptune Transiting a Quiet K Dwarf Star

IW ../submit_V2/abstract.txt ( Row 1 Col 1 6:48 Ctrl-K H for help We report the discovery by the HATSouth network of HATS-7b, a transiting Super-Neptune with a mass of 0.120+/-0.012MJ, a radius of 0.563+/-(0.046,0.034)RJ, and an orbital period of 3.1853days. The host star is a moderately bright (V=13.340+/-0.010mag, K_S=10.976+/-0.026mag) K dwarf star with a mass of 0.849+/-0.027Msun , a radius of 0.815+/-(0.049,-0.035)Rsun, and a metallicity of [Fe/H]=+0.250+/-0.080. The star is photometrically quiet to within the precision of the HATSouth measurements and has low RV jitter. HATS-7b is the second smallest radius planet discovered by a wide-field ground-based transit survey, and one of only a handful of Neptune-size planets with mass and radius determined to 10% precision. Theoretical modeling of HATS-7b yields a hydrogen-helium fraction of 18+/-4% (rock-iron core and H2-He envelope), or 9+/-4% (ice core and H2-He envelope), i.e.it has a composition broadly similar to that of Uranus and Neptune, and very different from that of Saturn, which has 75% of its mass in H2-He. Based on a sample of transiting exoplanets with accurately (<20%) determined parameters, we establish approximate power-law relations for the envelopes of the mass-density distribution of exoplanets. HATS-7b, which, together with the recently discovered HATS-8b, is one of the first two transiting super-Neptunes discovered in the Southern sky, is a prime target for additional follow-up observations with Southern hemisphere facilities to characterize the atmospheres of Super-Neptunes (which we define as objects with mass greater than that of Neptune, and smaller than halfway between that of Neptune and Saturn, i.e. 0.054 MJ<Mp<0.18 MJ).

preprint2015arXiv

HATS-8b: A Low-Density Transiting Super-Neptune

HATS-8b is a low density transiting super-Neptune discovered as part of the HATSouth project. The planet orbits its solar-like G dwarf host (V=14.03 $\pm$ 0.10 and T$_{eff}$ =5679 $\pm$ 50 K) with a period of 3.5839 d. HATS-8b is the third lowest mass transiting exoplanet to be discovered from a wide-field ground based search, and with a mass of 0.138 $\pm$ 0.019 M$_J$ it is approximately half-way between the masses of Neptune and Saturn. However HATS-8b has a radius of 0.873 (+0.123,-0.075) R$_J$, resulting in a bulk density of just 0.259 $\pm$ 0.091 g.cm$^{-3}$. The metallicity of the host star is super-Solar ([Fe/H]=0.210 $\pm$ 0.080), arguing against the idea that low density exoplanets form from metal-poor environments. The low density and large radius of HATS-8b results in an atmospheric scale height of almost 1000 km, and in addition to this there is an excellent reference star of near equal magnitude at just 19 arcsecond separation on the sky. These factors make HATS-8b an exciting target for future atmospheric characterization studies, particularly for long-slit transmission spectroscopy.

preprint2015arXiv

HATS-9b and HATS-10b: Two Compact Hot Jupiters in Field 7 of the K2 Mission

We report the discovery of two transiting extrasolar planets by the HATSouth survey. HATS-9b orbits an old (10.8 $\pm$ 1.5 Gyr) V=13.3 G dwarf star, with a period P = 1.9153 d. The host star has a mass of 1.03 M$_{\odot}$, radius of 1.503 R$_\odot$ and effective temperature 5366 $\pm$ 70 K. The planetary companion has a mass of 0.837 M$_J$, and radius of 1.065 R$_J$ yielding a mean density of 0.85 g cm$^{-3}$ . HATS-10b orbits a V=13.1 G dwarf star, with a period P = 3.3128 d. The host star has a mass of 1.1 M$_\odot$, radius of 1.11 R$_\odot$ and effective temperature 5880 $\pm$ 120 K. The planetary companion has a mass of 0.53 M$_J$, and radius of 0.97 R$_J$ yielding a mean density of 0.7 g cm$^{-3}$ . Both planets are compact in comparison with planets receiving similar irradiation from their host stars, and lie in the nominal coordinates of Field 7 of K2 but only HATS-9b falls on working silicon. Future characterisation of HATS-9b with the exquisite photometric precision of the Kepler telescope may provide measurements of its reflected light signature.

preprint2015arXiv

Hunting for brown dwarf binaries and testing atmospheric models with X-Shooter

The determination of the brown dwarf binary fraction may contribute to the understanding of the substellar formation mechanisms. Unresolved brown dwarf binaries may be revealed through their peculiar spectra or the discrepancy between optical and near-infrared spectral type classification. We obtained medium-resolution spectra of 22 brown dwarfs with these characteristics using the X-Shooter spectrograph at the VLT. We aimed to identify brown dwarf binary candidates, and to test if the BT-Settl 2014 atmospheric models reproduce their observed spectra. To find binaries spanning the L-T boundary, we used spectral indices and compared the spectra of the selected candidates to single spectra and synthetic binary spectra. We used synthetic binary spectra with components of same spectral type to determine as well the sensitivity of the method to this class of binaries. We identified three candidates to be combination of L plus T brown dwarfs. We are not able to identify binaries with components of similar spectral type. In our sample, we measured minimum binary fraction of $9.1^{+9.9}_{-3.0}$. From the best fit of the BT-Settl models 2014 to the observed spectra, we derived the atmospheric parameters for the single objects. The BT-Settl models were able to reproduce the majority of the SEDs from our objects, and the variation of the equivalent width of the RbI (794.8 nm) and CsI (852.0 nm) lines with the spectral type. Nonetheless, these models did not reproduce the evolution of the equivalent widths of the NaI (818.3 nm and 819.5 nm) and KI (1253 nm) lines with the spectral type.

preprint2015arXiv

Performance of the VLT Planet Finder SPHERE II. Data analysis and Results for IFS in laboratory

We present the performance of the Integral Field Spectrograph (IFS) of SPHERE, the high-contrast imager for the ESO VLT telescope designed to perform imaging and spectroscopy of extrasolar planets, obtained from tests performed at the Institute de Planétologie et d'Astrophysique de Grenoble facility during the integration phase of the instrument.} {The tests were performed using the instrument software purposely prepared for SPHERE. The output data were reduced applying the SPHERE data reduction and handling software, adding an improved spectral deconvolution procedure. To this aim, we prepared an alternative procedure for the spectral subtraction exploiting the principal components analysis algorithm. Moreover, a simulated angular differential imaging procedure was also implemented to estimate how the instrument performed once this procedure was applied at telescope. The capability of the IFS to faithfully retrieve the spectra of the detected faint companions was also considered.} {We found that the application of the updated version of the spectral deconvolution procedure alone, when the algorithm throughput is properly taken into account, gives us a $5σ$ limiting contrast of the order of 5$\times$$10^{-6}$ or slightly better. The further application of the angular differential imaging procedure on these data should allow us to improve the contrast by one order of magnitude down to around 7$\times$$10^{-7}$ at a separation of 0.3 arcsec. The application of a principal components analysis procedure that simultaneously uses spectral and angular data gives comparable results. Finally, we found that the reproducibility of the spectra of the detected faint companions is greatly improved when angular differential imaging is applied in addition to the spectral deconvolution.

preprint2015arXiv

Relationship between the column density distribution and evolutionary class of molecular clouds as viewed by ATLASGAL

We present the first study of the relationship between the column density distribution of molecular clouds within nearby Galactic spiral arms and their evolutionary status as measured from their stellar content. We analyze a sample of 195 molecular clouds located at distances below 5.5 kpc, identified from the ATLASGAL 870 micron data. We define three evolutionary classes within this sample: starless clumps, star-forming clouds with associated young stellar objects, and clouds associated with HII regions. We find that the N(H2) probability density functions (N-PDFs) of these three classes of objects are clearly different: the N-PDFs of starless clumps are narrowest and close to log-normal in shape, while star-forming clouds and HII regions exhibit a power-law shape over a wide range of column densities and log-normal-like components only at low column densities. We use the N-PDFs to estimate the evolutionary time-scales of the three classes of objects based on a simple analytic model from literature. Finally, we show that the integral of the N-PDFs, the dense gas mass fraction, depends on the total mass of the regions as measured by ATLASGAL: more massive clouds contain greater relative amounts of dense gas across all evolutionary classes.

preprint2015arXiv

Scattered light images of spiral arms in marginally gravitationally unstable discs with an embedded planet

Scattered light images of transition discs in the near-infrared often show non-axisymmetric structures in the form of wide-open spiral arms in addition to their characteristic low-opacity inner gap region. We study self-gravitating discs and investigate the influence of gravitational instability on the shape and contrast of spiral arms induced by planet-disc interactions. Two-dimensional non-isothermal hydrodynamical simulations including viscous heating and a cooling prescription are combined with three-dimensional dust continuum radiative transfer models for direct comparison to observations. We find that the resulting contrast between the spirals and the surrounding disc in scattered light is by far higher for pressure scale height variations, i.e. thermal perturbations, than for pure surface density variations. Self-gravity effects suppress any vortex modes and tend to reduce the opening angle of planet-induced spirals, making them more tightly wound. If the disc is only marginally gravitationally stable with a Toomre parameter around unity, an embedded massive planet (planet-to-star mass ratio of $10^{-2}$) can trigger gravitational instability in the outer disc. The spirals created by this instability and the density waves launched by the planet can overlap resulting in large-scale, more open spiral arms in the outer disc. The contrast of these spirals is well above the detection limit of current telescopes.

preprint2015arXiv

The GAPS Programme with HARPS-N@TNG VI: The Curious Case of TrES-4b

We revisit the TrES-4 system parameters based on high-precision HARPS-N radial-velocity measurements and new photometric light curves. A combined spectroscopic and photometric analysis allows us to determine a spectroscopic orbit with an amplitude $K=51\pm3$ m s$^{-1}$. The derived mass of TrES-4b is found to be $M_{\rm p} = 0.49\pm0.04 \rm M_{Jup}$, significantly lower than previously reported. Combined with the large radius ($R_{\rm p} = 1.84_{-0.09}^{+0.08} \rm R_{Jup}$) inferred from our analysis, TrES-4b becomes the second-lowest density transiting hot Jupiter known. We discuss several scenarios to explain the puzzling discrepancy in the mass of TrES-4b in the context of the exotic class of highly inflated transiting giant planets.

preprint2015arXiv

The LEECH Exoplanet Imaging Survey. Further constraints on the planet architecture of the HR 8799 system

Context. Astrometric monitoring of directly-imaged exoplanets allows the study of their orbital parameters and system architectures. Because most directly-imaged planets have long orbital periods (>20 AU), accurate astrometry is challenging when based on data acquired on timescales of a few years and usually with different instruments. The LMIRCam camera on the LBT is being used for the LEECH survey to search for and characterize young and adolescent exoplanets in L' band, including their system architectures. Aims. We first aim to provide a good astrometric calibration of LMIRCam. Then, we derive new astrometry, test the predictions of the orbital model of 8:4:2:1 mean motion resonance proposed by Goździewski & Migaszewski, and perform new orbital fitting of the HR 8799 bcde planets. We also present deep limits on a putative fifth planet interior to the known planets. Methods. We use observations of HR 8799 and the Theta1 Ori C field obtained during the same run in October 2013. Results. We first characterize the distortion of LMIRCam. We determine a platescale and a true north orientation for the images of 10.707 +/- 0.012 mas/pix and -0.430 +/- 0.076 deg, respectively. The errors on the platescale and true north orientation translate into astrometric accuracies at a separation of 1 of 1.1 mas and 1.3 mas, respectively. The measurements for all planets are usually in agreement within 3 sigma with the ephemeris predicted by Goździewski & Migaszewski. The orbital fitting based on the new astrometric measurements favors an architecture for the planetary system based on 8:4:2:1 mean motion resonance. The detection limits allow us to exclude a fifth planet slightly brighter/more massive than HR 8799 b at the location of the 2:1 resonance with HR 8799 e (~9.5 AU) and about twice as bright as HR 8799 cde at the location of the 3:1 resonance with HR 8799 e (~7.5 AU).

preprint2015arXiv

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

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.

preprint2015arXiv

THOR - The HI, OH, Recombination Line Survey of the Milky Way - The pilot study: HI observations of the giant molecular cloud W43

To study the atomic, molecular and ionized emission of Giant Molecular Clouds (GMCs), we have initiated a Large Program with the VLA: 'THOR - The HI, OH, Recombination Line survey of the Milky Way'. We map the 21cm HI line, 4 OH lines, 19 H_alpha recombination lines and the continuum from 1 to 2 GHz of a significant fraction of the Milky Way (l=15-67deg, |b|<1deg) at ~20" resolution. In this paper, we focus on the HI emission from the W43 star-formation complex. Classically, the HI 21cm line is treated as optically thin with column densities calculated under this assumption. This might give reasonable results for regions of low-mass star-formation, however, it is not sufficient to describe GMCs. We analyzed strong continuum sources to measure the optical depth, and thus correct the HI 21cm emission for optical depth effects and weak diffuse continuum emission. Hence, we are able to measure the HI mass of W43 more accurately and our analysis reveals a lower limit of M=6.6x10^6 M_sun, which is a factor of 2.4 larger than the mass estimated with the assumption of optically thin emission. The HI column densities are as high as N(HI)~150 M_sun/pc^2 ~ 1.9x10^22 cm^-2, which is an order of magnitude higher than for low mass star formation regions. This result challenges theoretical models that predict a threshold for the HI column density of ~10 M_sun/pc^2, at which the formation of molecular hydrogen should set in. By assuming an elliptical layered structure for W43, we estimate the particle density profiles. While at the cloud edge atomic and molecular hydrogen are well mixed, the center of the cloud is strongly dominated by molecular hydrogen. We do not identify a sharp transition between hydrogen in atomic and molecular form. Our results are an important characterization of the atomic to molecular hydrogen transition in an extreme environment and challenges current theoretical models.

preprint2014arXiv

A Statistical Analysis of SEEDS and Other High-Contrast Exoplanet Surveys: Massive Planets or Low-Mass Brown Dwarfs?

We conduct a statistical analysis of a combined sample of direct imaging data, totalling nearly 250 stars. The stars cover a wide range of ages and spectral types, and include five detections ($κ$ And b, two $\sim$60 M$_{\rm J}$ brown dwarf companions in the Pleiades, PZ Tel B, and CD$-$35 2722B). For some analyses we add a currently unpublished set of SEEDS observations, including the detections GJ 504b and GJ 758B. We conduct a uniform, Bayesian analysis of all stellar ages using both membership in a kinematic moving group and activity/rotation age indicators. We then present a new statistical method for computing the likelihood of a substellar distribution function. By performing most of the integrals analytically, we achieve an enormous speedup over brute-force Monte Carlo. We use this method to place upper limits on the maximum semimajor axis of the distribution function derived from radial-velocity planets, finding model-dependent values of $\sim$30--100 AU. Finally, we model the entire substellar sample, from massive brown dwarfs to a theoretically motivated cutoff at $\sim$5 M$_{\rm Jup}$, with a single power law distribution. We find that $p(M, a) \propto M^{-0.65\pm0.60} a^{-0.85\pm0.39}$ (1$σ$ errors) provides an adequate fit to our data, with 1.0--3.1\% (68\% confidence) of stars hosting 5--70 $M_{\rm Jup}$ companions between 10 and 100 AU. This suggests that many of the directly imaged exoplanets known, including most (if not all) of the low-mass companions in our sample, formed by fragmentation in a cloud or disk, and represent the low-mass tail of the brown dwarfs.

preprint2014arXiv

Exoplanet Science with the European Extremely Large Telescope. The Case for Visible and Near-IR Spectroscopy at High Resolution

Exoplanet science is booming. In 20 years our knowledge has expanded considerably, from the first discovery of a Hot Jupiter, to the detection of a large population of Neptunes and super-Earths, to the first steps toward the characterization of exoplanet atmospheres. Between today and 2025, the field will evolve at an even faster pace with the advent of several space-based transit search missions, ground-based spectrographs, high-contrast imaging facilities, and the James Webb Space Telescope. Especially the ESA M-class PLATO mission will be a game changer in the field. From 2024 onwards, PLATO will find transiting terrestrial planets orbiting within the habitable zones of nearby, bright stars. These objects will require the power of Extremely Large Telescopes (ELTs) to be characterized further. The technique of ground-based high-resolution spectroscopy is establishing itself as a crucial pathway to measure chemical composition, atmospheric structure and atmospheric circulation in transiting exoplanets. A high-resolution spectrograph covering the visible and near-IR domains, mounted on the European ELT, will be able to detect molecules such as water vapour, carbon dioxide and oxygen in the atmospheres of habitable planets under favourable circumstances. E-ELT HiRES is the perfect ground-based match to the PLATO space mission and represents a unique opportunity for Europe to lead the world into the era of exploration of exoplanets with habitable conditions. HiRES will also be extremely complementary to other E-ELT planned instruments specialising in different kinds of planets, such as METIS and EPICS.

preprint2014arXiv

Grain opacity and the bulk composition of extrasolar planets. I. Results from scaling the ISM opacity

The opacity due to grains in the envelope of a protoplanet regulates the accretion rate of gas during formation, thus the final bulk composition of planets with primordial H/He is a function of it. Observationally, for exoplanets with known mass and radius it is possible to estimate the bulk composition via internal structure models. We first determine the reduction factor of the ISM grain opacity f_opa that leads to gas accretion rates consistent with grain evolution models. We then compare the bulk composition of synthetic low-mass and giant planets at different f_opa with observations. For f_opa=1 (full ISM opacity) the synthetic low-mass planets have too small radii, i.e., too low envelope masses compared to observations. At f_opa=0.003, the value calibrated with the grain evolution models, synthetic and actual planets occupy similar mass-radius loci. The mean enrichment of giant planets relative to the host star as a function of planet mass M can be approximated as Z_p/Z_star = beta*(M/M_Jup)^alpha. We find alpha=-0.7 independent of f_opa in synthetic populations in agreement with the observational result (-0.71+-0.10). The absolute enrichment level decreases from beta=8.5 at f_opa=1 to 3.5 at f_opa=0. At f_opa=0.003 one finds beta=7.2 which is similar to the observational result (6.3+-1.0). We thus find observational hints that the opacity in protoplanetary atmospheres is much smaller than in the ISM even if the specific value of the grain opacity cannot be constrained here. The result for the enrichment of giant planets helps to distinguish core accretion and gravitational instability. In the simplest picture of core accretion where first a critical core forms and afterwards only gas is added, alpha=-1. If a core accretes all planetesimals inside the feeding zone, alpha=-2/3. The observational result lies between these values, pointing to core accretion as the formation mechanism.

preprint2014arXiv

HATS-4b: A Dense Hot-Jupiter Transiting a Super Metal-Rich G Star

We report the discovery by the HATSouth survey of HATS-4b, an extrasolar planet transiting a V=13.46 mag G star. HATS-4b has a period of P = 2.5167 d, mass of Mp = 1.32 Mj, radius of Rp = 1.02 Rj and density of rho_p = 1.55 +- 0.16 g/cm^3 ~ 1.24 rhoj. The host star has a mass of 1.00 Msun, a radius of 0.92 Rsun and a very high metallicity [Fe/H]= 0.43 +- 0.08. HATS-4b is among the densest known planets with masses between 1-2 Mj and is thus likely to have a significant content of heavy elements of the order of 75 Mearth. In this paper we present the data reduction, radial velocity measurement and stellar classification techniques adopted by the HATSouth survey for the CORALIE spectrograph. We also detail a technique to estimate simultaneously vsini and macroturbulence using high resolution spectra.

preprint2014arXiv

HATS-5b: A Transiting hot-Saturn from the HATSouth Survey

We report the discovery of HATS-5b, a transiting hot-Saturn orbiting a G type star, by the HAT-South survey. HATS-5b has a mass of Mp=0.24 Mj, radius of Rp=0.91 Rj, and transits its host star with a period of P=4.7634d. The radius of HATS-5b is consistent with both theoretical and empirical models. The host star has a V band magnitude of 12.6, mass of 0.94 Msun, and radius of 0.87 Rsun. The relatively high scale height of HATS-5b, and the bright, photometrically quiet host star, make this planet a favourable target for future transmission spectroscopy follow-up observations. We reexamine the correlations in radius, equilibrium temperature, and metallicity of the close-in gas-giants, and find hot Jupiter-mass planets to exhibit the strongest dependence between radius and equilibrium temperature. We find no significant dependence in radius and metallicity for the close-in gas-giant population.

preprint2014arXiv

HATS-6b: A Warm Saturn Transiting an Early M Dwarf Star, and a Set of Empirical Relations for Characterizing K and M Dwarf Planet Hosts

We report the discovery by the HATSouth survey of HATS-6b, an extrasolar planet transiting a V=15.2 mag, i=13.7 mag M1V star with a mass of 0.57 Msun and a radius of 0.57 Rsun. HATS-6b has a period of P = 3.3253 d, mass of Mp=0.32 Mjup, radius of Rp=1.00 Rjup, and zero-albedo equilibrium temperature of Teq=712.8+-5.1 K. HATS-6 is one of the lowest mass stars known to host a close-in gas giant planet, and its transits are among the deepest of any known transiting planet system. We discuss the follow-up opportunities afforded by this system, noting that despite the faintness of the host star, it is expected to have the highest K-band S/N transmission spectrum among known gas giant planets with Teq < 750 K. In order to characterize the star we present a new set of empirical relations between the density, radius, mass, bolometric magnitude, and V, J, H and K-band bolometric corrections for main sequence stars with M < 0.80 Msun, or spectral types later than K5. These relations are calibrated using eclipsing binary components as well as members of resolved binary systems. We account for intrinsic scatter in the relations in a self-consistent manner. We show that from the transit-based stellar density alone it is possible to measure the mass and radius of a ~0.6 Msun star to ~7% and ~2% precision, respectively. Incorporating additional information, such as the V-K color, or an absolute magnitude, allows the precision to be improved by up to a factor of two.

preprint2014arXiv

Impacts of planet migration models on planetary populations. Effects of saturation, cooling and stellar irradiation

Context: Several recent studies have found that planet migration in adiabatic discs differs significantly from migration in isothermal discs. Depending on the thermodynamic conditions, i.e., the effectiveness of radiative cooling, and the radial surface density profile, planets migrate inward or outward. Clearly, this will influence the semimajor axis - mass distribution of planets as predicted by population synthesis simulations. Aims: Our goal is to study the global effects of radiative cooling, viscous torque desaturation and gap opening as well as stellar irradiation on the tidal migration of a synthetic planet population. Methods: We combine results from several analytical studies and 3D hydrodynamic simulations in a new semi-analytical migration model for the application in our planet population synthesis calculations. Results: We find a good agreement of our model with torques obtained in a 3D radiative hydrodynamic simulations. We find three convergence zones in a typical disc, towards which planets migrate from the in- and outside, affecting strongly the migration behavior of low-mass planets. Interestingly, this leads to slow type II like migration behavior for low-mass planets captured in those zones even without an ad hoc migration rate reduction factor or a yet to be defined halting mechanism. This means that the new prescription of migration including non-isothermal effects makes the preciously widely used artificial migration rate reduction factor obsolete. Conclusions: Outward migration in parts of a disc makes some planets survive long enough to become massive. The convergence zones lead to a potentially observable accumulations of low-mass planets at certain semimajor axes. Our results indicate that further studies of the mass where the corotation torque saturates will be needed since its value has a major impact on the properties of planet populations.

preprint2014arXiv

Improving the astrometric performance of VLTI-PRIMA

In the summer of 2011, the first on-sky astrometric commissioning of PRIMA-Astrometry delivered a performance of 3 m'' for a 10 '' separation on bright objects, orders of magnitude away from its exoplanet requirement of 50 μ'' ~ 20 μ'' on objects as faint as 11 mag ~ 13 mag in K band. This contribution focuses on upgrades and characterizations carried out since then. The astrometric metrology was extended from the Coudé focus of the Auxillary Telescopes to their secondary mirror, in order to reduce the baseline instabilities and improve the astrometric performance. While carrying out this extension, it was realized that the polarization retardance of the star separator derotator had a major impact on both the astrometric metrology and the fringe sensors. A local compensation of this retardance and the operation on a symmetric baseline allowed a new astrometric commissioning. In October 2013, an improved astrometric performance of 160 μ'' was demonstrated, still short of the requirements. Instabilities in the astrometric baseline still appear to be the dominating factor. In preparation to a review held in January 2014, a plan was developed to further improve the astrometric and faint target performance of PRIMA Astrometry. On the astrometric aspect, it involved the extension of the internal longitudinal metrology to primary space, the design and implementation of an external baseline metrology, and the development of an astrometric internal fringes mode. On the faint target aspect, investigations of the performance of the fringe sensor units and the development of an AO system (NAOMI) were in the plan. Following this review, ESO decided to take a proposal to the April 2014 STC that PRIMA be cancelled, and that ESO resources be concentrated on ensuring that Gravity and Matisse are a success. This proposal was recommended by the STC in May 2014, and endorsed by ESO.

preprint2014arXiv

Millimetre spectral indices of transition disks and their relation to the cavity radius

Transition disks are protoplanetary disks with inner depleted dust cavities and excellent candidates to investigate the dust evolution under the existence of a pressure bump. A pressure bump at the outer edge of the cavity allows dust grains from the outer regions to stop their rapid inward migration towards the star and efficiently grow to millimetre sizes. Dynamical interactions with planet(s) have been one of the most exciting theories to explain the clearing of the inner disk. We look for evidence of the presence of millimetre dust particles in transition disks by measuring their spectral index with new and available photometric data. We investigate the influence of the size of the dust depleted cavity on the disk integrated millimetre spectral index. We present the 3mm photometric observations carried out with PdBI of four transition disks: LkHa330, UXTauA, LRLL31, and LRLL67. We use available values of their fluxes at 345GHz to calculate their spectral index, as well as the spectral index for a sample of twenty transition disks. We compare the observations with two kind of models. In the first set of models, we consider coagulation and fragmentation of dust in a disk in which a cavity is formed by a massive planet located at different positions. The second set of models assumes disks with truncated inner parts at different radius and with power-law dust size distributions, where the maximum size of grains is calculated considering turbulence as the source of destructive collisions. We show that the integrated spectral index is higher for transition disks than for regular protoplanetary disks. For transition disks, the probability that the measured spectral index is positively correlated with the cavity radius is 95%. High angular resolution imaging of transition disks is needed to distinguish between the dust trapping scenario and the truncated disk case.

preprint2014arXiv

New constraints on the formation and settling of dust in the atmospheres of young M and L dwarfs

We obtained medium-resolution near-infrared spectra of seven young M9.5-L3 dwarfs classified in the optical. We aim to confirm the low surface gravity of the objects in the NIR. We also test whether atmospheric models correctly represent the formation and the settling of dust clouds in the atmosphere of young late-M and L dwarfs. We used ISAAC at VLT to obtain the spectra of the targets. We compared them to those of mature and young BD, and young late-type companions to nearby stars with known ages, in order to identify and study gravity-sensitive features. We computed spectral indices weakly sensitive to the surface gravity to derive near-infrared spectral types. Finally, we found the best fit between each spectrum and synthetic spectra from the BT-Settl 2010 and 2013 models. Using the best fit, we derived the atmospheric parameters of the objects and identify which spectral characteristics the models do not reproduce. We confirmed that our objects are young BD and we found NIR spectral types in agreement with the ones determined at optical wavelengths. The spectrum of the L2-gamma dwarf 2MASSJ2322-6151 reproduces well the spectrum of the planetary mass companion 1RXS J1609-2105b. BT-Settl models fit the spectra and the 1-5 $μ$m SED of the L0-L3 dwarfs for temperatures between 1600-2000 K. But the models fail to reproduce the shape of the H band, and the NIR slope of some of our targets. This fact, and the best fit solutions found with super-solar metallicity are indicative of a lack of dust, in particular at high altitude, in the cloud models. The modeling of the vertical mixing and of the grain growth will be revised in the next version of the BT-Settl models. These revisions may suppress the remaining non-reproducibilities.

preprint2014arXiv

On the structure of the transition disk around TW Hya

For over a decade, the structure of the inner cavity in the transition disk of TW Hydrae has been a subject of debate. Modeling the disk with data obtained at different wavelengths has led to a variety of proposed disk structures. Rather than being inconsistent, the individual models might point to the different faces of physical processes going on in disks, such as dust growth and planet formation. Our aim is to investigate the structure of the transition disk again and to find to what extent we can reconcile apparent model differences. A large set of high-angular-resolution data was collected from near-infrared to centimeter wavelengths. We investigated the existing disk models and established a new self-consistent radiative-transfer model. A genetic fitting algorithm was used to automatize the parameter fitting. Simple disk models with a vertical inner rim and a radially homogeneous dust composition from small to large grains cannot reproduce the combined data set. Two modifications are applied to this simple disk model: (1) the inner rim is smoothed by exponentially decreasing the surface density in the inner ~3 AU, and (2) the largest grains (>100 um) are concentrated towards the inner disk region. Both properties can be linked to fundamental processes that determine the evolution of protoplanetary disks: the shaping by a possible companion and the different regimes of dust-grain growth, respectively. The full interferometric data set from near-infrared to centimeter wavelengths requires a revision of existing models for the TW Hya disk. We present a new model that incorporates the characteristic structures of previous models but deviates in two key aspects: it does not have a sharp edge at 4 AU, and the surface density of large grains differs from that of smaller grains. This is the first successful radiative-transfer-based model for a full set of interferometric data.

preprint2014arXiv

Performance of the VLT Planet Finder SPHERE I. Photometry and astrometry precision with IRDIS and IFS in laboratory

The new planet finder for the Very Large Telescope (VLT), the Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE), just had its first light in Paranal. A dedicated instrument for the direct detection of planets, SPHERE, is composed of a polametric camera in visible light, the Zurich IMager POLarimeter (ZIMPOL), and two near-infrared sub-systems: the Infra-Red Dual-beam Imager and Spectrograph (IRDIS), a multi-purpose camera for imaging, polarimetry, and long-slit spectroscopy, and the integral field spectrograph (IFS), an integral field spectrograph. We present the results obtained from the analysis of data taken during the laboratory integration and validation phase, after the injection of synthetic planets. Since no continuous field rotation could be performed in the laboratory, this analysis presents results obtained using reduction techniques that do not use the angular differential imaging (ADI) technique. To perform the simulations, we used the instrumental point spread function (PSF) and model spectra of L and T-type objects scaled in contrast with respect to the host star. We evaluated the expected error in astrometry and photometry as a function of the signal to noise of companions, after spectral differential imaging (SDI) reduction for IRDIS and spectral deconvolution (SD) or principal component analysis (PCA) data reductions for IFS. We deduced from our analysis, for example, that $β$Picb, a 12~Myr old planet of $\sim$10~\MJ and semi-major axis of 9--10 AU, would be detected with IRDIS with a photometric error of 0.16~mag and with a relative astrometric position error of 1.1~mas. With IFS, we could retrieve a spectrum with error bars of about 0.15~mag on each channel and astrometric relative position error of 0.6~mas. For a fainter object such as HR8799d, a 13~\MJ planet at a distance of 27~AU, IRDIS could obtain a relative astrometric error of 3~mas.

preprint2014arXiv

The architecture of the LkCa 15 transitional disk revealed by high-contrast imaging

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

preprint2014arXiv

The Moving Group Targets of the SEEDS High-Contrast Imaging Survey of Exoplanets and Disks: Results and Observations from the First Three Years

We present results from the first three years of observations of moving group targets in the SEEDS high-contrast imaging survey of exoplanets and disks using the Subaru telescope. We achieve typical contrasts of ~10^5 at 1" and ~10^6 beyond 2" around 63 proposed members of nearby kinematic moving groups. We review each of the kinematic associations to which our targets belong, concluding that five, βPictoris (~20 Myr), AB Doradus (~100 Myr), Columba (~30 Myr), Tucana-Horogium (~30 Myr), TW Hydrae (~10 Myr), are sufficiently well-defined to constrain the ages of individual targets. Somewhat less than half of our targets are high-probability members of one of these moving groups. For all of our targets, we combine proposed moving group membership with other age indicators where available, including Ca II HK emission, X-ray activity, and rotation period, to produce a posterior probability distribution of age. SEEDS observations discovered a substellar companion to one of our targets, κAnd, a late B star. We do not detect any other substellar companions, but do find seven new close binary systems, of which one still needs to be confirmed. A detailed analysis of the statistics of this sample, and of the companion mass constraints given our age probability distributions and exoplanet cooling models, will be presented in a forthcoming paper.

preprint2014arXiv

The VLT/NaCo large program to probe the occurrence of exoplanets and brown dwarfs at wide orbits: II- Survey description, results and performances

In anticipation of the VLT/SPHERE planet imager guaranteed time programs, we have conducted a preparatory survey of 86 stars between 2009 and 2013 in order to identify new faint comoving companions to ultimately carry out a comprehensive analysis of the occurence of giant planets and brown dwarf companions at wide (10-2000 AU) orbits around young, solar-type stars. We used NaCo at VLT to explore the occurrence rate of giant planets and brown dwarfs between typically 0.1 and 8''. Diffraction-limited observations in H-band combined with angular differential imaging enabled us to reach primary star-companion brightness ratios as small as 10-6 at 1.5''. 12 systems were resolved as new binaries, including the discovery of a new white dwarf companion to the star HD8049. Around 34 stars, at least one companion candidate was detected in the observed field of view. More than 400 faint sources were detected, 90% of them in 4 crowded fields. With the exception of HD8049B, we did not identify any new comoving companions. The survey also led to spatially resolved images of the thin debris disk around HD\,61005 that have been published earlier. Finally, considering the survey detection limits, we derive a preliminary upper limit on the frequency of giant planets for semi-major axes of [10,2000] AU: typically less than 15% between 100 and 500 AU, and less than 10% between 50 and 500 AU for exoplanets more massive than 5 MJup and 10 MJup respectively, considering a uniform input distribution and with a confidence level of 95%.

preprint2014arXiv

The VLT/NaCo Large program to probe the occurrence of exoplanets and brown dwarfs in wide orbits: I- Sample definition and characterization

Young, nearby stars are ideal targets to search for planets using the direct imaging technique. The determination of stellar parameters is crucial for the interpretation of imaging survey results particularly since the luminosity of substellar objects has a strong dependence on system age. We have conducted a large program with NaCo at the VLT in order to search for planets and brown dwarfs in wide orbits around 86 stars. A large fraction of the targets observed with NaCo were poorly investigated in the literature. We performed a study to characterize the fundamental properties (age, distance, mass) of the stars in our sample. To improve target age determinations, we compiled and analyzed a complete set of age diagnostics. We measured spectroscopic parameters and age diagnostics using dedicated observations acquired with FEROS and CORALIE spectrographs at La Silla Observatory. We also made extensive use of archival spectroscopic data and results available in the literature. Additionally, we exploited photometric time-series, available in ASAS and Super-WASP archives, to derive rotation period for a large fraction of our program stars. We provided updated characterization of all the targets observed in the VLT NaCo Large program, a survey designed to probe the occurrence of exoplanets and brown dwarfs in wide orbits. The median distance and age of our program stars are 64 pc and 100 Myr, respectively. Nearly all the stars have masses between 0.70 and 1.50sun, with a median value of 1.01 Msun. The typical metallicity is close to solar, with a dispersion that is smaller than that of samples usually observed in radial velocity surveys. Several stars are confirmed or proposed here to be members of nearby young moving groups. Eight spectroscopic binaries are identified.

preprint2013arXiv

Characterization of the gaseous companion κ Andromedae b: New Keck and LBTI high-contrast observations

We previously reported the direct detection of a low mass companion at a projected separation of 55+-2 AU around the B9 type star κ Andromedae. The properties of the system (mass ratio, separation) make it a benchmark for the understanding of the formation and evolution of gas giant planets and brown dwarfs on wide-orbits. We present new angular differential imaging (ADI) images of the Kappa Andromedae system at 2.146 (Ks), 3.776 (L'), 4.052 (NB 4.05) and 4.78 μm (M') obtained with Keck/NIRC2 and LBTI/LMIRCam, as well as more accurate near-infrared photometry of the star with the MIMIR instrument. We derive a more accurate J = 15.86 +- 0.21, H = 14.95 +- 0.13, Ks = 14.32 +- 0.09 mag for κ And b. We redetect the companion in all our high contrast observations. We confirm previous contrasts obtained at Ks and L' band. We derive NB 4.05 = 13.0 +- 0.2 and M' = 13.3 +- 0.3 mag and estimate Log10(L/Lsun) = -3.76 +- 0.06. We build the 1-5 microns spectral energy distribution of the companion and compare it to seven PHOENIX-based atmospheric models in order to derive Teff = 1900+100-200 K. Models do not set constrains on the surface gravity. ``Hot-start" evolutionary models predict masses of 14+25-2 MJup based on the luminosity and temperature estimates, and considering a conservative age range for the system (30+120-10 Myr). ``warm-start" evolutionary tracks constrain the mass to M >= 11 MJup. Therefore, the mass of κ Andromedae b mostly falls in the brown-dwarf regime, due to remaining uncertainties in age and mass-luminosity models. According to the formation models, disk instability in a primordial disk could account for the position and a wide range of plausible masses of κ And b.

preprint2013arXiv

Direct Imaging Detection of Methane in the Atmosphere of GJ 504 b

Most exoplanets detected by direct imaging so far have been characterized by relatively hot (> ~1000 K) and cloudy atmospheres. A surprising feature in some of their atmospheres has been a distinct lack of methane, possibly implying non-equilibrium chemistry. Recently, we reported the discovery of a planetary companion to the Sun-like star GJ 504 using Subaru/HiCIAO within the SEEDS survey. The planet is substantially colder (<600 K) than previously imaged planets, and has indications of fewer clouds, which implies that it represents a new class of planetary atmospheres with expected similarities to late T-type brown dwarfs in the same temperature range. If so, one might also expect the presence of significant methane absorption, which is characteristic of such objects. Here, we report the detection of deep methane absorption in the atmosphere of GJ 504 b, using the Spectral Differential Imaging mode of HiCIAO to distinguish the absorption feature around 1.6 um. We also report updated JHK photometry based on new Ks-band data and a re-analysis of the existing data. The results support the notion that GJ 504 b has atmospheric properties distinct from other imaged exoplanets, and will become a useful reference object for future planets in the same temperature range.

preprint2013arXiv

Direct Imaging of a Cold Jovian Exoplanet in Orbit around the Sun-like Star GJ 504

Several exoplanets have recently been imaged at wide separations of >10 AU from their parent stars. These span a limited range of ages (<50 Myr) and atmospheric properties, with temperatures of 800--1800 K and very red colors (J - H > 0.5 mag), implying thick cloud covers. Furthermore, substantial model uncertainties exist at these young ages due to the unknown initial conditions at formation, which can lead to an order of magnitude of uncertainty in the modeled planet mass. Here, we report the direct imaging discovery of a Jovian exoplanet around the Sun-like star GJ 504, detected as part of the SEEDS survey. The system is older than all other known directly-imaged planets; as a result, its estimated mass remains in the planetary regime independent of uncertainties related to choices of initial conditions in the exoplanet modeling. Using the most common exoplanet cooling model, and given the system age of 160 [+350, -60] Myr, GJ 504 b has an estimated mass of 4 [+4.5, -1.0] Jupiter masses, among the lowest of directly imaged planets. Its projected separation of 43.5 AU exceeds the typical outer boundary of ~30 AU predicted for the core accretion mechanism. GJ 504 b is also significantly cooler (510 [+30, -20] K) and has a bluer color (J-H = -0.23 mag) than previously imaged exoplanets, suggesting a largely cloud-free atmosphere accessible to spectroscopic characterization. Thus, it has the potential of providing novel insights into the origins of giant planets, as well as their atmospheric properties.

preprint2013arXiv

HATS-2b: A transiting extrasolar planet orbiting a K-type star showing starspot activity

We report the discovery of HATS-2b, the second transiting extrasolar planet detected by the HATSouth survey. HATS-2b is moving on a circular orbit around a V=13.6 mag, K-type dwarf star (GSC 6665-00236), at a separation of 0.0230 \pm 0.0003 AU and with a period of 1.3541 days. The planetary parameters have been robustly determined using a simultaneous fit of the HATSouth, MPG/ESO~2.2\,m/GROND, Faulkes Telescope South/Spectral transit photometry and MPG/ESO~2.2\,m/FEROS, Euler~1.2\,m/CORALIE, AAT~3.9\,m/CYCLOPS radial-velocity measurements. HATS-2b has a mass of 1.37 \pm 0.16 M_J, a radius of 1.14 \pm 0.03 R_J and an equilibrium temperature of 1567 \pm 30 K. The host star has a mass of 0.88 \pm 0.04 M_Sun, radius of 0.89 \pm 0.02 R_Sun and shows starspot activity. We characterized the stellar activity by analysing two photometric follow-up transit light curves taken with the GROND instrument, both obtained simultaneously in four optical bands (covering the wavelength range of 3860-9520 Å). The two light curves contain anomalies compatible with starspots on the photosphere of the parent star along the same transit chord.

preprint2013arXiv

HATS-3b: An inflated hot Jupiter transiting an F-type star

We report the discovery by the HATSouth survey of HATS-3b, a transiting extrasolar planet orbiting a V=12.4 F-dwarf star. HATS-3b has a period of P = 3.5479d, mass of Mp = 1.07MJ, and radius of Rp = 1.38RJ. Given the radius of the planet, the brightness of the host star, and the stellar rotational velocity (vsini = 9.0km/s), this system will make an interesting target for future observations to measure the Rossiter-McLaughlin effect and determine its spin-orbit alignment. We detail the low/medium-resolution reconnaissance spectroscopy that we are now using to deal with large numbers of transiting planet candidates produced by the HATSouth survey. We show that this important step in discovering planets produces logg and Teff parameters at a precision suitable for efficient candidate vetting, as well as efficiently identifying stellar mass eclipsing binaries with radial velocity semi-amplitudes as low as 1 km/s.

preprint2013arXiv

High-fidelity view of the structure and fragmentation of the high-mass, filamentary IRDC G11.11-0.12

Star formation in molecular clouds is intimately linked to their internal mass distribution. We present an unprecedentedly detailed analysis of the column density structure of a high-mass, filamentary molecular cloud, namely IRDC G11.11-0.12 (G11). We use two novel column density mapping techniques: high-resolution (FWHM=2", or ~0.035 pc) dust extinction mapping in near- and mid-infrared, and dust emission mapping with the Herschel satellite. These two completely independent techniques yield a strikingly good agreement, highlighting their complementarity and robustness. We first analyze the dense gas mass fraction and linear mass density of G11. We show that G11 has a top heavy mass distribution and has a linear mass density (M_l ~ 600 Msun pc^{-1}) that greatly exceeds the critical value of a self-gravitating, non-turbulent cylinder. These properties make G11 analogous to the Orion A cloud, despite its low star-forming activity. This suggests that the amount of dense gas in molecular clouds is more closely connected to environmental parameters or global processes than to the star-forming efficiency of the cloud. We then examine hierarchical fragmentation in G11 over a wide range of size-scales and densities. We show that at scales 0.5 pc > l > 8 pc, the fragmentation of G11 is in agreement with that of a self-gravitating cylinder. At scales smaller than l < 0.5 pc, the results agree better with spherical Jeans' fragmentation. One possible explanation for the change in fragmentation characteristics is the size-scale-dependent collapse time-scale that results from the finite size of real molecular clouds: at scales l < 0.5 pc, fragmentation becomes sufficiently rapid to be unaffected by global instabilities.

preprint2013arXiv

Imaging Discovery of the Debris Disk Around HIP 79977

We present Subaru/HiCIAO H-band high-contrast images of the debris disk around HIP 79977, whose pres- ence was recently inferred from an infrared excess. Our images resolve the disk for the first time, allowing characterization of its shape, size, and dust grain properties. We use angular differential imaging (ADI) to reveal the disk geometry in unpolarized light out to a radius of ~2", as well as polarized differential imaging (PDI) to measure the degree of scattering polarization out to ~1.5". In order to strike a favorable balance between suppression of the stellar halo and conservation of disk flux, we explore the application of principal component analysis (PCA) to both ADI and reference star subtraction. This allows accurate forward modeling of the effects of data reduction on simulated disk images, and thus direct comparison with the imaged disk. The resulting best-fit values and well-fitting intervals for the model parameters are a surface brightness power-law slope of S_out = -3.2 [-3.6,-2.9], an inclination of i = 84° [81°,86°], a high Henyey-Greenstein forward-scattering parameter of g = 0.45 [0.35, 0.60], and a non-significant disk-star offset of u = 3.0 [-1.5, 7.5] AU = 24 [-13, 61] mas along the line of nodes. Furthermore, the tangential linear polarization along the disk rises from ~10% at 0.5" to ~45% at 1.5". These measurements paint a consistent picture of a disk of dust grains produced by collisional cascades and blown out to larger radii by stellar radiation pressure.

preprint2013arXiv

Kepler-91b: a planet at the end of its life. Planet and giant host star properties via light-curve variations

The evolution of planetary systems is intimately linked to the evolution of their host star. Our understanding of the whole planetary evolution process is based on the large planet diversity observed so far. To date, only few tens of planets have been discovered orbiting stars ascending the Red Giant Branch. Although several theories have been proposed, the question of how planets die remains open due to the small number statistics. In this work we study the giant star Kepler-91 (KOI-2133) in order to determine the nature of a transiting companion. This system was detected by the Kepler Space Telescope. However, its planetary confirmation is needed. We confirm the planetary nature of the object transiting the star Kepler-91 by deriving a mass of $ M_p=0.88^{+0.17}_{-0.33} ~M_{\rm Jup}$ and a planetary radius of $R_p=1.384^{+0.011}_{-0.054} ~R_{\rm Jup}$. Asteroseismic analysis produces a stellar radius of $R_{\star}=6.30\pm 0.16 ~R_{\odot}$ and a mass of $M_{\star}=1.31\pm 0.10 ~ M_{\odot} $. We find that its eccentric orbit ($e=0.066^{+0.013}_{-0.017}$) is just $1.32^{+0.07}_{-0.22} ~ R_{\star}$ away from the stellar atmosphere at the pericenter. Kepler-91b could be the previous stage of the planet engulfment, recently detected for BD+48 740. Our estimations show that Kepler-91b will be swallowed by its host star in less than 55 Myr. Among the confirmed planets around giant stars, this is the planetary-mass body closest to its host star. At pericenter passage, the star subtends an angle of $48^{\circ}$, covering around 10% of the sky as seen from the planet. The planetary atmosphere seems to be inflated probably due to the high stellar irradiation.

preprint2013arXiv

Low-velocity shocks traced by extended SiO emission along the W43 ridges: witnessing the formation of young massive clusters

The formation of high-mass stars is tightly linked to that of their parental clouds. We here focus on the high-density parts of W43, a molecular cloud undergoing an efficient event of formation. The cloud structure is studied with a column density image derived from Herschel continuum maps obtained at 70, 160, 250, 350, and 500 micron. We identify two high-column density filamentary clouds, quoted as the W43-MM1 and W43-MM2 ridges, which both account for 1.5x10^4 Msun gas mass above 10^23 cm-2 and within areas of 5 and 14pc^2, respectively. We used the N_2H^+ 1--0 line to confirm that the W43-MM1 and W43-MM2 ridges are structures coherent in velocity and gravitationally bound, despite their large velocity dispersion and ~5 kms line widths. The most intriguing result of the W43 large program is the bright wide-spread SiO 2--1 emission: 1--11 K kms$ stretching an area of ~28 pc^2. Concentrated toward the W43-MM1 and W43-MM2 ridges and their immediate surroundings, it leads to a total luminosity of L_SiO 2-1 ~4 10^4 K kms kpc^2pc^2. We measured a steep relation between the luminosity and velocity extent of the SiO~2--1 lines and propose to use it to distinguish the low-velocity shocks observed here from the more classical high-velocity ones associated with outflows of high-mass young stellar objects. We used state-of-the-art shock models to demonstrate that low-velocity (<10 kms^-1) shocks with a small amount (10%) of Si atoms initially in gas phase or in grain mantles can explain the observed SiO column density in W43. The spatial and velocity overlaps between the ridges high-density gas (n_H2>10^4-10^5 cm^-3) and the shocked SiO gas suggests that ridges could be forming via colliding flows driven by gravity and accompanied by low-velocity shocks. This mechanism may be the initial conditions for the formation of young massive clusters in these ridges.

preprint2013arXiv

New Techniques for High-Contrast Imaging with ADI: the ACORNS-ADI SEEDS Data Reduction Pipeline

We describe Algorithms for Calibration, Optimized Registration, and Nulling the Star in Angular Differential Imaging (ACORNS-ADI), a new, parallelized software package to reduce high-contrast imaging data, and its application to data from the SEEDS survey. We implement several new algorithms, including a method to register saturated images, a trimmed mean for combining an image sequence that reduces noise by up to ~20%, and a robust and computationally fast method to compute the sensitivity of a high-contrast observation everywhere on the field-of-view without introducing artificial sources. We also include a description of image processing steps to remove electronic artifacts specific to Hawaii2-RG detectors like the one used for SEEDS, and a detailed analysis of the Locally Optimized Combination of Images (LOCI) algorithm commonly used to reduce high-contrast imaging data. ACORNS-ADI is written in python. It is efficient and open-source, and includes several optional features which may improve performance on data from other instruments. ACORNS-ADI requires minimal modification to reduce data from instruments other than HiCIAO. It is freely available for download at www.github.com/t-brandt/acorns-adi under a BSD license.

preprint2013arXiv

Parallax measurements of cool brown dwarfs

Accurate parallax measurements allow us to determine physical properties of brown dwarfs, and help us to constrain evolutionary and atmospheric models, break the age-mass degeneracy and reveal unresolved binaries. We measured absolute trigonometric parallaxes and proper motions of 6 cool brown dwarfs using background galaxies to establish an absolute reference frame. We derive the absolute J-mag. The six T brown dwarfs in our sample have spectral types between T2.5 and T7.5 and magnitudes in J between 13.9 and 18.0, with photometric distances below 25 pc. The observations were taken in the J-band with the Omega-2000 camera on the 3.5 m telescope at Calar Alto, during a time period of 27 months, between March 2011 and June 2013. The number of epochs varied between 11 and 12 depending on the object. The reduction of the astrometric measurements was carried out with respect to the field stars. The relative parallax and proper motions were transformed into absolute measurements using the background galaxies in our fields. We obtained absolute parallaxes for our six brown dwarfs with a precision between 3 and 6 mas. We compared our results in a color-magnitude diagram with other brown dwarfs with determined parallax and with the BT-Settl 2012 atmospheric models. For four of the six targets we found a good agreement in luminosity with objects of similar spectral types. We obtained an improved accuracy in the parallaxes and proper motions in comparison to previous works. The object 2MASS J11061197+2754225 is more than 1 mag overluminous in all bands pointing to binarity or high order multiplicity.

preprint2013arXiv

The Mass-Radius Relationship for Very Low Mass Stars: Four New Discoveries from the HATSouth Survey

We report the discovery of four transiting F-M binary systems with companions between 0.1-0.2 Msun in mass by the HATSouth survey. These systems have been characterised via a global analysis of the HATSouth discovery data, combined with high-resolution radial velocities and accurate transit photometry observations. We determined the masses and radii of the component stars using a combination of two methods: isochrone fitting of spectroscopic primary star parameters, and equating spectroscopic primary star rotation velocity with spin-orbit synchronisation. These new very low mass companions are HATS550-016B (0.110 -0.006/+0.005 Msun, 0.147 -0.004/+0.003 Rsun), HATS551-019B (0.17 -0.01/+0.01 Msun, 0.18 -0.01/+0.01 Rsun), HATS551-021B (0.132 -0.005/+0.014 Msun, 0.154 -0.008/+0.006 Rsun), HATS553-001B (0.20 -0.02/+0.01 Msun, 0.22 -0.01/+0.01 Rsun). We examine our sample in the context of the radius anomaly for fully-convective low mass stars. Combining our sample with the 13 other well-studied very low mass stars, we find a tentative 5% systematic deviation between the measured radii and theoretical isochrone models.

preprint2013arXiv

Towards a complete stellar mass function of the Hyades. I. Pan-STARRS1 optical observations of the low-mass stellar content

The Hyades cluster is an ideal target to study the dynamical evolution of a star cluster over the entire mass range due to its intermediate age and proximity to the Sun. We wanted to extend the Hyades mass function towards lower masses down to 0.1 Msol and to use the full three-dimensional spatial information to characterize the dynamical evolution of the cluster. We performed a kinematic and photometric selection using the PPMXL and Pan-STARRS1 sky surveys, to search for cluster members up to 30 pc from the cluster centre. We determined our detection efficiency and field star contamination rate to derive the cluster luminosity and mass functions down to masses of 0.1 Msol. The thorough astrometric and photometric constraints minimized the contamination. A minimum spanning tree algorithm was used to quantify the mass segregation. We discovered 43 new Hyades member candidates with velocity perpendicular to the Hyades motion up to 2 km/s. They have mass estimates between 0.43 and 0.09 Msol, for a total mass of 10 Msol. This doubles the number of Hyades candidates with masses smaller than 0.15 Msol. We provide an additional list of 11 possible candidates with velocity perpendicular to the Hyades motion up to 4 km/s. The cluster is significantly mass segregated. The extension of the mass function towards lower masses provided an even clearer signature than estimated in the past. We also identified as likely Hyades member an L0 dwarf previously assumed to be a field dwarf. Finally we question the membership of a number of previously published candidates, including a L2.5-type dwarf.

preprint2012arXiv

ATLASGAL - Compact source catalogue: 330 < l < 21 degrees

Context. The APEX Telescope Large Area Survey of the Galaxy (ATLASGAL) is the first systematic survey of the inner Galactic plane in the sub-millimetre. The observations were carried out with the Large APEX Bolometer Camera (LABOCA), an array of 295 bolometers observing at 870 microns (345 GHz). Aims. Here we present a first version of the compact source catalogue extracted from this survey. This catalogue provides an unbiased database of dusty clumps in the inner Galaxy. Methods. The construction of this catalogue was made using the source extraction routine SExtractor. We have cross-associated the obtained sources with the IRAS and MSX catalogues, in order to constrain their nature. Results. We have detected 6639 compact sources in the range from 330 < l < 21 degrees and |b| < 1.5 degrees. The catalogue has a 99% completeness for sources with a peak flux above 6 sigma, which corresponds to a flux density of ~0.4 Jy/beam. The parameters extracted for sources with peak fluxes below the 6 sigma completeness threshold should be used with caution. Tests on simulated data find the uncertainty in the flux measurement to be ~12%, however, in more complex regions the flux values can be overestimated by a factor of 2 due to the additional background emission. Using a search radius of 30" we found that 40% of ATLASGAL compact sources are associated with an IRAS or MSX point source, but, ~50% are found to be associated with MSX 21 microns fluxes above the local background level, which is probably a lower limit to the actual number of sources associated with star formation. Conclusions. Although infrared emission is found towards the majority of the clumps detected, this catalogue is still likely to include a significant number of clumps that are devoid of star formation activity and therefore excellent candidates for objects in the coldest, earliest stages of (high-mass) star formation.

preprint2012arXiv

Characterization of exoplanets from their formation I: Models of combined planet formation and evolution

A first characterization of many exoplanets has recently been achieved by the observational determination of their radius. For some planets, a measurement of the luminosity has also been possible, with many more directly imaged planets expected in the future. The statistical characterization of exoplanets through their mass-radius and mass-luminosity diagram is thus becoming possible. This is for planet formation and evolution theory of similar importance as the mass-distance diagram. Our aim in this and a companion paper is to extend our formation model into a coupled formation and evolution model. We want to calculate in a self-consistent way all basic characteristics (M,a,R,L) of a planet and use the model for population synthesis calculations. Here we show how we solve the structure equations describing the gaseous envelope not only during the early formation phase, but also during gas runaway accretion, and during the evolutionary phase at constant mass on Gyr timescales. We then study the in situ formation and evolution of Jupiter, the mass-radius relationship of giants, the influence of the core mass on the radius and the luminosity both in the "hot start" and the "cold start" scenario. We put special emphasis on the comparison with other models. We find that our results agree very well with those of more complex models, despite a number of simplifications. The upgraded model yields the most important characteristics of a planet from its beginning as a seed embryo to a Gyr old planet. This is the case for all planets in a synthetic planetary population. Therefore, we can now use self-consistently the statistical constraints coming from all major observational techniques. This is important in a time where different techniques yield constraints on very diverse sub-populations of planets, and where its is challenging to put all these constraints together in one coherent picture.

preprint2012arXiv

Characterization of exoplanets from their formation II: The planetary mass-radius relationship

The research of exoplanets has entered an era in which we characterize extrasolar planets. This has become possible with measurements of radii and luminosities. Meanwhile, radial velocity surveys discover also very low-mass planets. Uniting all this observational data into one coherent picture to better understand planet formation is an important, but difficult undertaking. Our approach is to develop a model which can make testable predictions for all these observational methods. We continue to describe how we have extended our formation model into a self-consistently coupled formation and evolution model. We show how we calculate the internal structure of the solid core and radiogenic heating. We also improve the protoplanetary disk model. Finally, we conduct population synthesis calculations. We present how the planetary mass-radius relationship of planets with primordial H/He envelopes forms and evolves in time. The basic shape of the M-R relation can be understood from the core accretion model. Low-mass planets cannot bind massive envelopes, while super-critical cores necessarily trigger runway gas accretion, leading to "forbidden" zones in the M-R plane. For a given mass, there is a considerable diversity of radii. We compare the synthetic M-R relation with the observed one, finding good agreement for a>0.1 AU. The synthetic radius distribution is characterized by a strong increase towards small R, and a second, lower local maximum at ~1 Jovian radius. The increase towards small radii reflects the increase of the mass function towards low M. The second local maximum is due to the fact that radii are nearly independent of mass for giant planets. A comparison of the synthetic radius distribution with Kepler data shows agreement for R>2 Earth radii, but divergence for smaller radii. We predict that in the next few years, Kepler should find the second, local maximum at ~1 Jovian radius.

preprint2012arXiv

Chemistry in Disks VIII: the CS molecule as an analytic tracer of turbulence in disks

Turbulence is thought to be a key driver of the evolution of protoplanetary disks, regulating the mass accretion process, the transport of angular momentum, and the growth of dust particles. We intend to determine the magnitude of the turbulent motions in the outer parts (> 100 AU) of the disk surrounding DM Tau. Turbulent motions can be constrained by measuring the nonthermal broadening of line emission from heavy molecules. We used the IRAM Plateau de Bure interferometer to study emission from the CS molecule in the disk of DM Tau. High spatial (1.4 x 1 ") and spectral resolution (0.126 km/s) CS J=3-2 images provide constraints on the molecule distribution and velocity structure of the disk. A low sensitivity CS J=5-4 image was used in conjunction to evaluate the excitation conditions. We analyzed the data in terms of two parametric disk models, and compared the results with detailed time-dependent chemical simulations. The CS data confirm the relatively low temperature suggested by observations of other simple molecules. The intrinsic linewidth derived from the CS J=3-2 data is much larger than expected from pure thermal broadening. The magnitude of the derived nonthermal component depends only weakly on assumptions about the location of the CS molecules with respect to the disk plane. Our results indicate turbulence with a Mach number around 0.4 - 0.5 in the molecular layer. Geometrical constraints suggest that this layer is located near one scale height, in reasonable agreement with chemical model predictions.

preprint2012arXiv

CID: Chemistry In Disks VII. First detection of HC3N in protoplanetary disks

Molecular line emission from protoplanetary disks is a powerful tool to constrain their physical and chemical structure. Nevertheless, only a few molecules have been detected in disks so far. We take advantage of the enhanced capabilities of the IRAM 30m telescope by using the new broad band correlator (FTS) to search for so far undetected molecules in the protoplanetary disks surrounding the TTauri stars DM Tau, GO Tau, LkCa 15 and the Herbig Ae star MWC 480. We report the first detection of HC3N at 5 sigma in the GO Tau and MWC 480 disks with the IRAM 30-m, and in the LkCa 15 disk (5 sigma), using the IRAM array, with derived column densities of the order of 10^{12}cm^{-2}. We also obtain stringent upper limits on CCS (N < 1.5 x 10^{12} cm^{-3}). We discuss the observational results by comparing them to column densities derived from existing chemical disk models (computed using the chemical code Nautilus) and based on previous nitrogen and sulfur-bearing molecule observations. The observed column densities of HC3N are typically two orders of magnitude lower than the existing predictions and appear to be lower in the presence of strong UV flux, suggesting that the molecular chemistry is sensitive to the UV penetration through the disk. The CCS upper limits reinforce our model with low elemental abundance of sulfur derived from other sulfur-bearing molecules (CS, H2S and SO).

preprint2012arXiv

Direct Imaging Discovery of a `Super-Jupiter' Around the late B-Type Star Kappa And

We present the direct imaging discovery of an extrasolar planet, or possible low-mass brown dwarf, at a projected separation of 55 +/- 2 AU (1.058 +/- 0.007 arcsec) from the B9-type star Kappa And. The planet was detected with Subaru/HiCIAO during the SEEDS survey, and confirmed as a bound companion via common proper motion measurements. Observed near-infrared magnitudes of J = 16.3 +/- 0.3, H = 15.2 +/- 0.2, Ks = 14.6 +/- 0.4, and L' = 13.12 +/- 0.09 indicate a temperature of ~1700 K. The galactic kinematics of the host star are consistent with membership in the Columba association, implying a corresponding age of 30 +20 -10 Myr. The system age, combined with the companion photometry, points to a model-dependent companion mass ~12.8 MJup. The host star's estimated mass of 2.4-2.5 Msun places it among the most massive stars ever known to harbor an extrasolar planet or low-mass brown dwarf. While the mass of the companion is close to the deuterium burning limit, its mass ratio, orbital separation, and likely planet-like formation scenario imply that it may be best defined as a `Super-Jupiter' with properties similar to other recently discovered companions to massive stars.

preprint2012arXiv

Discovery of Small-Scale Spiral Structures in the Disk of SAO 206462 (HD 135344B): Implications for the Physical State of the Disk from Spiral Density Wave Theory

We present high-resolution, H-band, imaging observations, collected with Subaru/HiCIAO, of the scattered light from the transitional disk around SAO 206462 (HD 135344B). Although previous sub-mm imagery suggested the existence of the dust-depleted cavity at r~46AU, our observations reveal the presence of scattered light components as close as 0.2" (~28AU) from the star. Moreover, we have discovered two small-scale spiral structures lying within 0.5" (~70AU). We present models for the spiral structures using the spiral density wave theory, and derive a disk aspect ratio of h~0.1, which is consistent with previous sub-mm observations. This model can potentially give estimates of the temperature and rotation profiles of the disk based on dynamical processes, independently from sub-mm observations. It also predicts the evolution of the spiral structures, which can be observable on timescales of 10-20 years, providing conclusive tests of the model. While we cannot uniquely identify the origin of these spirals, planets embedded in the disk may be capable of exciting the observed morphology. Assuming that this is the case, we can make predictions on the locations and, possibly, the masses of the unseen planets. Such planets may be detected by future multi-wavelengths observations.

preprint2012arXiv

Extrasolar planet population synthesis IV. Correlations with disk metallicity, mass and lifetime

Context. This is the fourth paper in a series showing the results of planet population synthesis calculations. Aims. Our goal in this paper is to systematically study the effects of important disk properties, namely disk metallicity, mass and lifetime on fundamental planetary properties. Methods. For a large number of protoplanetary disks we calculate a population of planets with our core accretion formation model including planet migration and disk evolution. Results. We find a large number of correlations: Regarding the planetary initial mass function, metallicity, disk mass and disk lifetime have different roles: For high [Fe/H], giant planets are more frequent. For high disk masses, giant planets are more massive. For long disk lifetimes, giant planets are both more frequent and massive. At low metallicities, very massive giant planets cannot form, but otherwise giant planet mass and metallicity are uncorrelated. In contrast, planet masses and disk gas masses are correlated. The sweet spot for giant planet formation is at 5 AU. In- and outside this distance, higher planetesimals surface densities are necessary. Low metallicities can be compensated by high disk masses, and vice versa, but not ad infinitum. At low metallicities, giant planets only form outside the ice line, while at high metallicities, giant planet formation occurs throughout the disk. The extent of migration increases with disk mass and lifetime and usually decreases with metallicity. No clear correlation of metallicity and the semimajor axis of giant planets exists because in low [Fe/H] disks, planets start further out, but migrate more, whereas for high [Fe/H] they start further in, but migrate less. Close-in low mass planets have a lower mean metallicity than Hot Jupiters. Conclusions. The properties of protoplanetary disks are decisive for the properties of planets, and leave many imprints.

preprint2012arXiv

HATS-1b: The First Transiting Planet Discovered by the HATSouth Survey

We report the discovery of HATS-1b, a transiting extrasolar planet orbiting the moderately bright V=12.05 G dwarf star GSC 6652-00186, and the first planet discovered by HATSouth, a global network of autonomous wide-field telescopes. HATS-1b has a period P~3.4465 d, mass Mp~1.86MJ, and radius Rp~1.30RJ. The host star has a mass of 0.99Msun, and radius of 1.04Rsun. The discovery light curve of HATS-1b has near continuous coverage over several multi-day periods, demonstrating the power of using a global network of telescopes to discover transiting planets.

preprint2012arXiv

HATSouth: a global network of fully automated identical wide-field telescopes

HATSouth is the world's first network of automated and homogeneous telescopes that is capable of year-round 24-hour monitoring of positions over an entire hemisphere of the sky. The primary scientific goal of the network is to discover and characterize a large number of transiting extrasolar planets, reaching out to long periods and down to small planetary radii. HATSouth achieves this by monitoring extended areas on the sky, deriving high precision light curves for a large number of stars, searching for the signature of planetary transits, and confirming planetary candidates with larger telescopes. HATSouth employs 6 telescope units spread over 3 locations with large longitude separation in the southern hemisphere (Las Campanas Observatory, Chile; HESS site, Namibia; Siding Spring Observatory, Australia). Each of the HATSouth units holds four 0.18m diameter f/2.8 focal ratio telescope tubes on a common mount producing an 8.2x8.2 arcdeg field, imaged using four 4Kx4K CCD cameras and Sloan r filters, to give a pixel scale of 3.7 arcsec/pixel. The HATSouth network is capable of continuously monitoring 128 square arc-degrees. We present the technical details of the network, summarize operations, and present weather statistics for the 3 sites. On average each of the 6 HATSouth units has conducted observations on ~500 nights over a 2-year time period, yielding a total of more than 1million science frames at 4 minute integration time, and observing ~10.65 hours per day on average. We describe the scheme of our data transfer and reduction from raw pixel images to trend-filtered light curves and transiting planet candidates. Photometric precision reaches ~6 mmag at 4-minute cadence for the brightest non-saturated stars at r~10.5. We present detailed transit recovery simulations to determine the expected yield of transiting planets from HATSouth. (abridged)

preprint2012arXiv

High-Contrast NIR Polarization Imaging of MWC480

One of the key predictions of modeling from the IR excess of Herbig Ae stars is that for protoplanetary disks, where significant grain growth and settling has occurred, the dust disk has flattened to the point that it can be partially or largely shadowed by the innermost material at or near the dust sublimation radius. When the self-shadowing has already started, the outer disk is expected to be detected in scattered light only in the exceptional cases that the scale height of the dust disk at the sublimation radius is smaller than usual. High-contrast imaging combined with the IR spectral energy distribution allow us to measure the degree of flattening of the disk, as well as to determine the properties of the outer disk. We present polarimetric differential imaging in $H$ band obtained with Subaru/HiCIAO of one such system, MWC 480. The HiCIAO data were obtained at a historic minimum of the NIR excess. The disk is detected in scattered light from 0\farcs2-1\farcs0 (27.4-137AU). Together with the marginal detection of the disk from 1998 February 24 by HST/NICMOS, our data constrain the opening half angle for the disk to lie between 1.3$\leqθ\leq 2.2^\circ$. When compared with similar measures in CO for the gas disk from the literature, the dust disk subtends only $\sim$30% of the gas disk scale height (H/R$\sim$0.03). Such a dust disk is a factor of 5-7 flatter than transitional disks, which have structural signatures that giant planets have formed.

preprint2012arXiv

High-Resolution Near-Infrared Polarimetry of a Circumstellar Disk around UX Tau A

We present H-band polarimetric imagery of UX Tau A taken with HiCIAO/AO188 on the Subaru Telescope. UX Tau A has been classified as a pre-transitional disk object, with a gap structure separating its inner and outer disks. Our imagery taken with the 0.15 (21 AU) radius coronagraphic mask has revealed a strongly polarized circumstellar disk surrounding UX Tau A which extends to 120 AU, at a spatial resolution of 0.1 (14 AU). It is inclined by 46 \pm 2 degree as the west side is nearest. Although SED modeling and sub-millimeter imagery suggested the presence of a gap in the disk, with the inner edge of the outer disk estimated to be located at 25 - 30 AU, we detect no evidence of a gap at the limit of our inner working angle (23 AU) at the near-infrared wavelength. We attribute the observed strong polarization (up to 66 %) to light scattering by dust grains in the disk. However, neither polarization models of the circumstellar disk based on Rayleigh scattering nor Mie scattering approximations were consistent with the observed azimuthal profile of the polarization degrees of the disk. Instead, a geometric optics model of the disk with nonspherical grains with the radii of 30 micron meter is consistent with the observed profile. We suggest that the dust grains have experienced frequent collisional coagulations and have grown in the circumstellar disk of UX Tau A.

preprint2012arXiv

Narrow-angle astrometry with PRIMA

The Extrasolar Planet Search with PRIMA project (ESPRI) aims at characterising and detecting extrasolar planets by measuring the host star's reflex motion using the narrow-angle astrometry capability of the PRIMA facility at the Very Large Telescope Interferometer. A first functional demonstration of the astrometric mode was achieved in early 2011. This marked the start of the astrometric commissioning phase with the purpose of characterising the instrument's performance, which ultimately has to be sufficient for exoplanet detection. We show results obtained from the observation of bright visual binary stars, which serve as test objects to determine the instrument's astrometric precision, its accuracy, and the plate scale. Finally, we report on the current status of the ESPRI project, in view of starting its scientific programme.

preprint2012arXiv

Subaru Imaging of Asymmetric Features in a Transitional Disk in Upper Scorpius

We report high-resolution (0.07 arcsec) near-infrared polarized intensity images of the circumstellar disk around the star 2MASS J16042165-2130284 obtained with HiCIAO mounted on the Subaru 8.2 m telescope. We present our $H$-band data, which clearly exhibits a resolved, face-on disk with a large inner hole for the first time at infrared wavelengths. We detect the centrosymmetric polarization pattern in the circumstellar material as has been observed in other disks. Elliptical fitting gives the semimajor axis, semiminor axis, and position angle (P.A.) of the disk as 63 AU, 62 AU, and -14 $^{\circ}$, respectively. The disk is asymmetric, with one dip located at P.A.s of $\sim85^{\circ}$. Our observed disk size agrees well with a previous study of dust and CO emission at submillimeter wavelength with Submillimeter Array. Hence, the near-infrared light is interpreted as scattered light reflected from the inner edge of the disk. Our observations also detect an elongated arc (50 AU) extending over the disk inner hole. It emanates at the inner edge of the western side of the disk, extending inward first, then curving to the northeast. We discuss the possibility that the inner hole, the dip, and the arc that we have observed may be related to the existence of unseen bodies within the disk.

preprint2012arXiv

The missing cavities in the SEEDS polarized scattered light images of transitional protoplanetary disks: a generic disk model

Transitional circumstellar disks around young stellar objects have a distinctive infrared deficit around 10 microns in their Spectral Energy Distributions (SED), recently measured by the Spitzer Infrared Spectrograph (IRS), suggesting dust depletion in the inner regions. These disks have been confirmed to have giant central cavities by imaging of the submillimeter (sub-mm) continuum emission using the Submillimeter Array (SMA). However, the polarized near-infrared scattered light images for most objects in a systematic IRS/SMA cross sample, obtained by HiCIAO on the Subaru telescope, show no evidence for the cavity, in clear contrast with SMA and Spitzer observations. Radiative transfer modeling indicates that many of these scattered light images are consistent with a smooth spatial distribution for micron-sized grains, with little discontinuity in the surface density of the micron-sized grains at the cavity edge. Here we present a generic disk model that can simultaneously account for the general features in IRS, SMA, and Subaru observations. Particularly, the scattered light images for this model are computed, which agree with the general trend seen in Subaru data. Decoupling between the spatial distributions of the micron-sized dust and mm-sized dust inside the cavity is suggested by the model, which, if confirmed, necessitates a mechanism, such as dust filtration, for differentiating the small and big dust in the cavity clearing process. Our model also suggests an inwardly increasing gas-to-dust-ratio in the inner disk, and different spatial distributions for the small dust inside and outside the cavity, echoing the predictions in grain coagulation and growth models.

preprint2012arXiv

The Structure of Pre-transitional Protoplanetary Disks I: Radiative Transfer Modeling of the Disk+Cavity in the PDS 70 system

Through detailed radiative transfer modeling, we present a disk+cavity model to simultaneously explain both the SED and Subaru H-band polarized light imaging for the pre-transitional protoplanetary disk PDS 70. Particularly, we are able to match not only the radial dependence, but also the absolute scale, of the surface brightness of the scattered light. Our disk model has a cavity 65 AU in radius, which is heavily depleted of sub-micron-sized dust grains, and a small residual inner disk which produces a weak but still optically thick NIR excess in the SED. To explain the contrast of the cavity edge in the Subaru image, a factor of ~1000 depletion for the sub-micron-sized dust inside the cavity is required. The total dust mass of the disk may be on the order of 1e-4 M_sun, only weakly constrained due to the lack of long wavelength observations and the uncertainties in the dust model. The scale height of the sub-micron-sized dust is ~6 AU at the cavity edge, and the cavity wall is optically thick in the vertical direction at H-band. PDS 70 is not a member of the class of (pre-)transitional disks identified by Dong et al. (2012), whose members only show evidence of the cavity in the millimeter-sized dust but not the sub-micron-sized dust in resolved images. The two classes of (pre-)transitional disks may form through different mechanisms, or they may just be at different evolution stages in the disk clearing process.

preprint2012arXiv

Triggered/sequential star formation? A multi-phase ISM study around the prominent IRDC G18.93-0.03

G18.93-0.03 is a prominent dust complex within an 0.8deg long filament, with the molecular clump G18.93/m being IR dark from near IR wavelength up to 160mu. Spitzer composite images show an IR bubble spatially associated with G18.93. We use GRS 13CO and IRAM 30m H13CO+ data to disentangle the spatial structure of the region. From ATLASGAL submm data we calculate the gas mass, while we use the H13CO+ line width to estimate its virial mass. Using HERSCHEL data we produce temperature maps from fitting the SED. With the MAGPIS 20cm and SuperCOSMOS Halpha data we trace the ionized gas, and the VGPS HI survey provides information on the atomic hydrogen gas. We show that the bubble is spatially associated with G18.93, located at a kinematic near distance of 3.6kpc. With 280Msun, the most massive clump within G18.93 is G18.93/m. The virial analysis shows that it may be gravitationally bound and has neither Spitzer young stellar objects nor mid-IR point sources within. Fitting the SED reveals a temperature distribution that decreases towards its center, but heating from the ionizing source puts it above the general ISM temperature. We find that the bubble is filled by HII gas, ionized by an O8.5 star. Between the ionizing source and the IR dark clump G18.93/m we find a layered structure, from ionized to atomic to molecular hydrogen, revealing a PDR. Furthermore, we identify an additional velocity component within the bubble's 8mu emission rim at the edge of the infrared dark cloud and speculate that it might be shock induced by the expanding HII region. While the elevated temperature allows for the build-up of larger fragments, and the shock induced velocity component may lead to additional turbulent support, we do not find conclusive evidence that the massive clump G18.93/m is prone to collapse because of the expanding HII region.

preprint2011arXiv

A benchmark for multi-conjugated AO: VLT-MAD observations of the young massive cluster Trumpler 14

MAD is the first multi-conjugated adaptive optics system at the VLT. We present H and Ks observations of the young massive cluster Trumpler 14 revealing the power of MCAO systems by providing a homogeneous Strehl ratio over a large field of view. Mean Strehl ratios of 6.0 and 5.9 per cent with maximum Strehl ratios of 9.8 and 12.6 per cent in H and Ks, respectively, show significant improvement of the spatial PSF stability compared to single-conjugated adaptive optics systems. Photometry of our observations cover a dynamic range of ~10 mag including 2-3 times more sources than comparable seeing-limited observations. The colour-magnitude diagram reveals that the massive cluster originated in a recent starburst-like event 1+/-0.5 Myr ago. We tentatively detect hints for an older population of 3 Myr suggesting that low intensity star formation might have been going on in the HII region for a few Myr. We derive the luminosity function and mass function between 0.1 M_sun and 3.2 M_sun and identify a change of the power law slope of the mass function at m_c~0.53(+0.12/-0.10) M_sun. The MF appears shallow with power law slopes of Gamma1=-0.50+/-0.11 above m_c and Gamma2=0.63+/-0.32 below m_c.

preprint2011arXiv

A Planetary Companion around a Metal-Poor Star with Extragalactic Origin

We report the detection of a planetary companion around HIP 13044, a metal-poor star on the red Horizontal Branch. The detection is based on radial velocity observations with FEROS, a high-resolution spectrograph at the 2.2-m MPG/ESO telescope, located at ESO La Silla observatory in Chile. The periodic radial velocity variation of P = 16.2 days can be distinguished from the periods of the stellar activity indicators. We computed a minimum planetary mass of 1.25 MJup and an orbital semi-major axis of 0.116 AU for the planet. This discovery is unique in three aspects: First, it is the first planet detection around a star with a metallicity much lower than few percent of the solar value; second, the planet host star resides in a stellar evolutionary stage that is still unexplored in the exoplanet surveys; third, the star HIP 13044 belongs to one of the most significant stellar halo streams in the solar neighborhood, implying an extragalactic origin of the planetary system HIP 13044 in a disrupted former satellite of the Milky Way.

preprint2011arXiv

Application of recent results on the orbital migration of low mass planets: convergence zones

Previous models of the combined growth and migration of protoplanets needed large ad hoc reduction factors for the type I migration rate as found in the isothermal approximation. In order to eliminate these factors, a simple semi-analytical model is presented that incorporates recent results on the migration of low mass planets in non-isothermal disks. It allows for outward migration. The model is used to conduct planetary populations synthesis calculations. Two points with zero torque are found in the disks. Planets migrate both in- and outward towards these convergence zones. They could be important for accelerating planetary growth by concentrating matter in one point. We also find that the updated type I migration models allow the formation of both close-in low mass planets, but also of giant planets at large semimajor axes. The problem of too rapid migration is significantly mitigated.

preprint2011arXiv

Direct Imaging of Fine Structures in Giant Planet Forming Regions of the Protoplanetary Disk around AB Aurigae

We report high-resolution 1.6 $\micron$ polarized intensity ($PI$) images of the circumstellar disk around the Herbig Ae star AB Aur at a radial distance of 22 AU ($0."15$) up to 554 AU (3.$"$85), which have been obtained by the high-contrast instrument HiCIAO with the dual-beam polarimetry. We revealed complicated and asymmetrical structures in the inner part ($\lesssim$140 AU) of the disk, while confirming the previously reported outer ($r$ $\gtrsim$200 AU) spiral structure. We have imaged a double ring structure at $\sim$40 and $\sim$100 AU and a ring-like gap between the two. We found a significant discrepancy of inclination angles between two rings, which may indicate that the disk of AB Aur is warped. Furthermore, we found seven dips (the typical size is $\sim$45 AU or less) within two rings as well as three prominent $PI$ peaks at $\sim$40 AU. The observed structures, including a bumpy double ring, a ring-like gap, and a warped disk in the innermost regions, provide essential information for understanding the formation mechanism of recently detected wide-orbit ($r$ $>$20 AU) planets.

preprint2011arXiv

EChO - Exoplanet Characterisation Observatory

A dedicated mission to investigate exoplanetary atmospheres represents a major milestone in our quest to understand our place in the universe by placing our Solar System in context and by addressing the suitability of planets for the presence of life. EChO -the Exoplanet Characterisation Observatory- is a mission concept specifically geared for this purpose. EChO will provide simultaneous, multi-wavelength spectroscopic observations on a stable platform that will allow very long exposures. EChO will build on observations by Hubble, Spitzer and groundbased telescopes, which discovered the first molecules and atoms in exoplanetary atmospheres. EChO will simultaneously observe a broad enough spectral region -from the visible to the mid-IR- to constrain from one single spectrum the temperature structure of the atmosphere and the abundances of the major molecular species. The spectral range and resolution are tailored to separate bands belonging to up to 30 molecules to retrieve the composition and temperature structure of planetary atmospheres. The target list for EChO includes planets ranging from Jupiter-sized with equilibrium temperatures Teq up to 2000 K, to those of a few Earth masses, with Teq ~300 K. We have baselined a dispersive spectrograph design covering continuously the 0.4-16 micron spectral range in 6 channels (1 in the VIS, 5 in the IR), which allows the spectral resolution to be adapted from several tens to several hundreds, depending on the target brightness. The instrument will be mounted behind a 1.5 m class telescope, passively cooled to 50 K, with the instrument structure and optics passively cooled to ~45 K. EChO will be placed in a grand halo orbit around L2. We have also undertaken a first-order cost and development plan analysis and find that EChO is easily compatible with the ESA M-class mission framework.

preprint2011arXiv

Ex Lupi from Quiescence to Outburst: Exploring the LTE Approach in Modelling Blended H2O and OH Mid-Infrared Emission

We present a comparison of archival Spitzer spectra of the strongly variable T Tauri EX Lupi, observed before and during its 2008 outburst. We analyze the mid-infrared emission from gas-phase molecules thought to originate in a circumstellar disk. In quiescence the emission shows a forest of H2O lines, highly excited OH lines, and the Q branches of the organics C2H2, HCN, and CO2, similar to the emission observed toward several T Tauri systems. The outburst emission shows instead remarkable changes: H2O and OH line fluxes increase, new OH, H2, and HI transitions are detected, and organics are no longer seen. We adopt a simple model of a single-temperature slab of gas in local thermal equilibrium, a common approach for molecular analyses of Spitzer spectra, and derive the excitation temperature, column density, and emitting area of H2O and OH. We show how model results strongly depend on the selection of emission lines fitted, and that spectrally-resolved observations are essential for a correct interpretation of the molecular emission from disks, particularly in the case of water. Using H2O lines that can be approximated as thermalized to a single temperature, our results are consistent with a column density decrease in outburst while the emitting area of warm gas increases. A rotation diagram analysis suggests that the OH emission can be explained with two temperature components, which remarkably increase in column density in outburst. The relative change of H2O and OH emission suggests a key role for UV radiation in the disk surface chemistry.

preprint2011arXiv

Images of the Extended Outer Regions of the Debris Ring Around HR 4796 A

We present high-contrast images of HR 4796 A taken with Subaru/HiCIAO in H-band, resolving the debris disk in scattered light. The application of specialized angular differential imaging methods (ADI) allows us to trace the inner edge of the disk with high precision, and reveals a pair of "streamers" extending radially outwards from the ansae. Using a simple disk model with a power-law surface brightness profile, we demonstrate that the observed streamers can be understood as part of the smoothly tapered outer boundary of the debris disk, which is most visible at the ansae. Our observations are consistent with the expected result of a narrow planetesimal ring being ground up in a collisional cascade, yielding dust with a wide range of grain sizes. Radiation forces leave large grains in the ring and push smaller grains onto elliptical, or even hyperbolic trajectories. We measure and characterize the disk's surface brightness profile, and confirm the previously suspected offset of the disk's center from the star's position along the ring's major axis. Furthermore, we present first evidence for an offset along the minor axis. Such offsets are commonly viewed as signposts for the presence of unseen planets within a disk's cavity. Our images also offer new constraints on the presence of companions down to the planetary mass regime (~9 Jupiter masses at 0.5", ~3 Jupiter masses at 1").

preprint2011arXiv

Mass reservoirs surrounding massive infrared dark clouds: A view by near-infrared dust extinction

Context: Infrared Dark Clouds (IRDCs) harbor progenitors of high-mass stars. Little is known of the parental molecular clouds of the IRDCs. Aims: We demonstrate the feasibility of the near-infrared (NIR) dust extinction mapping in tracing the parental molecular clouds of IRDCs at the distances of D = 2.5 - 8 kpc. Methods: We derive NIR extinction maps for 10 prominent IRDC complexes using a color-excess mapping technique and NIR data from the UKIDSS/Galactic Plane Survey. We compare the resulting maps to the 13CO emission line data, to the 8 μm dust opacity data, and to the millimeter dust emission data. We derive distances for the clouds by comparing the observed NIR source densities to the Besancon stellar distribution model and compare them to the kinematic distance estimates. Results: The NIR extinction maps provide a view to the IRDC complexes over the dynamical range of Av = 2 - 40 mag, in spatial resolution of 30". The NIR extinction data correlate well with the 13CO data and probe a similar gas component, but also extend to higher column densities. The NIR data reveal a wealth of extended structures surrounding the dense gas traced by the 8 μm shadowing features and sub-mm dust emission, showing that the clouds contain typically > 10 times more mass than traced by those tracers. The IRDC complexes of our sample contain relatively high amount of high-column density material, and their cumulative column density distributions resemble active nearby star-forming clouds like Orion rather than less active clouds like California. Conclusions: NIR dust extinction data provide a new powerful tool to probe the mass distribution of the parental molecular clouds of IRDCs up to the distances of D = 8 kpc. This encourages for deeper NIR observations of IRDCs, because the sensitivity and resolution of the data can be directly enhanced with dedicated observations.

preprint2011arXiv

Near-Infrared Multi-Band Photometry of the Substellar Companion GJ 758 B

GJ 758 B is a cold (~600K) companion to a Sun-like star at 29 AU projected separation, which was recently detected with high-contrast imaging. Here we present photometry of the companion in seven photometric bands from Subaru/HiCIAO, Gemini/NIRI and Keck/NIRC2, providing a rich sampling of the spectral energy distribution in the 1-5 micron wavelength range. A clear detection at 1.58 micron combined with an upper limit at 1.69 micron shows methane absorption in the atmosphere of the companion. The mass of the companion remains uncertain, but an updated age estimate indicates that the most likely mass range is ~30-40 Mjup. In addition, we present an updated astrometric analysis that imposes tighter constraints on GJ 758 B's orbit and identifies the proposed second candidate companion, "GJ 758 C", as a background star.

preprint2011arXiv

Probing the evolution of molecular cloud structure II: From chaos to confinement

We present an analysis of the large-scale molecular cloud structure and of the stability of clumpy structures in nearby molecular clouds. In our recent work, we identified a structural transition in molecular clouds by studying the probability distributions of gas column densities in them. In this paper, we further examine the nature of this transition. The transition takes place at the visual extinction of A_V^tail = 2-4 mag, or equivalently, at Σ^tail = 40-80 Ms pc^{-2}. The clumps identified above this limit have wide ranges of masses and sizes, but a remarkably constant mean volume density of n = 10^3 cm^{-3}. This is 5-10 times larger than the density of the medium surrounding the clumps. By examining the stability of the clumps, we show that they are gravitationally unbound entities, and that the external pressure from the parental molecular cloud is a significant source of confining pressure for them. Then, the structural transition at A_V^tail may be linked to a transition between this population and the surrounding medium. The star formation rates in the clouds correlate strongly with the total mass in the clumps, i.e, with the mass above A_V^tail, dropping abruptly below that threshold. These results imply that the formation of pressure confined clumps introduces a prerequisite for star formation. Furthermore, they give a physically motivated explanation for the recently reported relation between the star formation rates and the amount of dense material in molecular clouds. Likewise, they give rise to a natural threshold for star formation at A_V^tail.

preprint2011arXiv

VLT/NACO astrometry of the HR8799 planetary system. L'-band observations of the three outer planets

HR8799 is so far the only directly imaged multiple exoplanet system. The orbital configuration would, if better known, provide valuable insight into the formation and dynamical evolution of wide-orbit planetary systems. We present L'-band observations of the HR8799 system obtained with NACO at VLT, adding to the astrometric monitoring of the planets HR8799b, c and d. We investigate how well the two simple cases of (i) a circular orbit and (ii) a face-on orbit fit the astrometric data for HR8799d over a total time baseline of ~2 years. The results indicate that the orbit of HR8799d is inclined with respect to our line of sight, and suggest that the orbit is slightly eccentric or non-coplanar with the outer planets and debris disk.

preprint2011arXiv

W43: the closest molecular complex of the Galactic Bar?

We used a large database extracted from Galaxy-wide surveys of H {\scriptsize I}, $^{13}$CO 1-0, $8 \micron$ and $870 \micron$ continuum to trace diffuse atomic gas, low- to medium-density molecular gas, high-density molecular gas, and star formation activity which we complemented by dedicated $^{12}$CO 2--1, 3--2 observations of the W43 region. From the detailed 3D (space-space-velocity) analysis of the molecular and atomic cloud tracers through the region and despite its wide velocity range (\emph{FWHM}$\sim22.3 \kms$ around $\vlsr\sim 95.9 \kms$), we identified W43 as a large (equivalent diameter $\sim 140$ pc) and coherent complex of molecular clouds which is surrounded by an atomic gas envelope (equivalent diameter $\sim 290$ pc). We measured the total mass of this newly-identified molecular complex ($M_{\tiny total}\sim 7.1 \times 10^6 \msun$), the mass contained in dense $870 \micron$ clumps ($<5$ pc dense cloud structures, $M_{\tiny clumps}\sim 8.4 \times 10^5 \msun$) and conclude that W43 is particularly massive and concentrated. The distance we assume for the W43 complex is 6 kpc from the Sun, which may place it at the meeting point of the Scutum-Centaurus (or Scutum-Crux) Galactic arm and the Bar, a dynamically complex region where high-velocity streams could easily collide. The star formation rate of W43 is suggested not to be steady but it is increasing from $\, \sim 0.01 \msun\, {yr}^{-1}$ (measured from its $8\,μ$m luminosity) to $\sim 0.1 \msun\, {yr}^{-1}$ (measured from its molecular content). From the global properties of W43, we claim it is an extreme molecular complex in the Milky Way and it could even be forming starburst clusters in the near future.

preprint2010arXiv

A new benchmark T8-9 brown dwarf and a couple of new mid-T dwarfs from the UKIDSS DR5+ LAS

Benchmark brown dwarfs are those objects for which fiducial constraints are available, including effective temperature, parallax, age, metallicity. We searched for new cool brown dwarfs in 186 sq.deg. of the new area covered by the data release DR5+ of the UKIDSS Large Area Survey. Follow-up optical and near-infrared broad-band photometry, and methane imaging of four promising candidates, revealed three objects with distinct methane absorption, typical of mid- to late-T dwarfs, and one possibly T4 dwarf. The latest-type object, classified as T8-9, shares its large proper motion with Ross 458 (BD+13o2618), an active M0.5 binary which is 102" away, forming a hierarchical low-mass star+brown dwarf system. Ross 458C has an absolute J-band magnitude of 16.4, and seems overluminous, particularly in the K band, compared to similar field brown dwarfs. We estimate the age of the system to be less than 1 Gyr, and its mass to be as low as 14 Jupiter masses for the age of 1 Gyr. At 11.4 pc, this new late T benchmark dwarf is a promising target to constrain the evolutionary and atmospheric models of very low-mass brown dwarfs. We present proper motion measurements for our targets and for 13 known brown dwarfs. Two brown dwarfs have velocities typical of the thick disk and may be old brown dwarfs.

preprint2010arXiv

CHESS, Chemical Herschel surveys of star forming regions:Peering into the protostellar shock L1157-B1

The outflow driven by the low-mass class 0 protostar L1157 is the prototype of the so-called chemically active outflows. The bright bowshock B1 in the southern outflow lobe is a privileged testbed of magneto-hydrodynamical (MHD) shock models, for which dynamical and chemical processes are strongly interdependent. We present the first results of the unbiased spectral survey of the L1157-B1 bowshock, obtained in the framework of the key program "Chemical Herschel Surveys of Star Forming Regions" (CHESS). The main aim is to trace the warm and chemically enriched gas and to infer the excitation conditions in the shock region. The CO 5-4 and H2O lines have been detected at high-spectral resolution in the unbiased spectral survey of the HIFI-Band 1b spectral window (555-636 GHz), presented by Codella et al. in this volume. Complementary ground-based observations in the submm window help establish the origin of the emission detected in the main-beam of HIFI, and the physical conditions in the shock.}{Both lines exhibit broad wings, which extend to velocities much higher than reported up to now. We find that the molecular emission arises from two regions with distinct physical conditions: an extended, warm (100K), dense (3e5 cm-3) component at low-velocity, which dominates the water line flux in Band~1; a secondary component in a small region of B1 (a few arcsec) associated with high-velocity, hot (> 400 K) gas of moderate density ((1.0-3.0)e4 cm-3), which appears to dominate the flux of the water line at 179mu observed with PACS. The water abundance is enhanced by two orders of magnitude between the low- and the high-velocity component, from 8e-7 up to 8e-5. The properties of the high-velocity component agree well with the predictions of steady-state C-shock models.

preprint2010arXiv

Dust-temperature of an isolated star-forming cloud: Herschel observations of the Bok globule CB244

We present Herschel observations of the isolated, low-mass star-forming Bok globule CB244. It contains two cold sources, a low-mass Class 0 protostar and a starless core, which is likely to be prestellar in nature, separated by 90 arcsec (~ 18000 AU). The Herschel data sample the peak of the Planck spectrum for these sources, and are therefore ideal for dust-temperature and column density modeling. With these data and a near-IR extinction map, the MIPS 70 micron mosaic, the SCUBA 850 micron map, and the IRAM 1.3 mm map, we model the dust-temperature and column density of CB244 and present the first measured dust-temperature map of an entire star-forming molecular cloud. We find that the column-averaged dust-temperature near the protostar is ~ 17.7 K, while for the starless core it is ~ 10.6K, and that the effect of external heating causes the cloud dust-temperature to rise to ~ 17 K where the hydrogen column density drops below 10^21 cm^-2. The total hydrogen mass of CB244 (assuming a distance of 200 pc) is 15 +/- 5 M_sun. The mass of the protostellar core is 1.6 +/- 0.1 M_sun and the mass of the starless core is 5 +/- 2 M_sun, indicating that ~ 45% of the mass in the globule is participating in the star-formation process.

preprint2010arXiv

First detection of ND in the solar-mass protostar IRAS16293-2422

In the past decade, much progress has been made in characterising the processes leading to the enhanced deuterium fractionation observed in the ISM and in particular in the cold, dense parts of star forming regions such as protostellar envelopes. Very high molecular D/H ratios have been found for saturated molecules and ions. However, little is known about the deuterium fractionation in radicals, even though simple radicals often represent an intermediate stage in the formation of more complex, saturated molecules. The imidogen radical NH is such an intermediate species for the ammonia synthesis in the gas phase. Herschel/HIFI represents a unique opportunity to study the deuteration and formation mechanisms of such species, which are not observable from the ground. We searched here for the deuterated radical ND in order to determine the deuterium fractionation of imidogen and constrain the deuteration mechanism of this species. We observed the solar-mass Class 0 protostar IRAS16293-2422 with the heterodyne instrument HIFI as part of the Herschel key programme CHESS (Chemical HErschel Surveys of Star forming regions). The deuterated form of the imidogen radical ND was detected and securely identified with 2 hyperfine component groups of its fundamental transition in absorption against the continuum background emitted from the nascent protostar. The 3 groups of hyperfine components of its hydrogenated counterpart NH were also detected in absorption. We derive a very high deuterium fractionation with an [ND]/[NH] ratio of between 30 and 100%. The deuterium fractionation of imidogen is of the same order of magnitude as that in other molecules, which suggests that an efficient deuterium fractionation mechanism is at play. We discuss two possible formation pathways for ND, by means of either the reaction of N+ with HD, or deuteron/proton exchange with NH.

preprint2010arXiv

Herschel/HIFI observations of spectrally resolved methylidyne signatures toward the high-mass star-forming core NGC6334I

In contrast to extensively studied dense star-forming cores, little is known about diffuse gas surrounding star-forming regions. We study molecular gas in the high-mass star-forming region NGC6334I, which contains diffuse, quiescent components that are inconspicuous in widely used molecular tracers such as CO. We present Herschel/HIFI observations of CH toward NGC6334I observed as part of the CHESS key program. HIFI resolves the hyperfine components of its J=3/2-1/2 transition, observed in both emission and absorption. The CH emission appears close to the systemic velocity of NGC6334I, while its measured linewidth of 3 km/s is smaller than previously observed in dense gas tracers such as NH3 and SiO. The CH abundance in the hot core is 7 10^-11, two to three orders of magnitude lower than in diffuse clouds. While other studies find distinct outflows in, e.g., CO and H2O toward NGC6334I, we do not detect outflow signatures in CH. To explain the absorption signatures, at least two absorbing components are needed at -3.0 and +6.5 km/s with N(CH)=7 10^13 and 3 10^13 cm^-2. Two additional absorbing clouds are found at +8.0 and 0.0 km/s, both with N(CH)=2 10^13 cm^-2. Turbulent linewidths for the four absorption components vary between 1.5 and 5.0 km/s in FWHM. We constrain physical properties of our CH clouds by matching our CH absorbers with other absorption signatures. In the hot core, molecules such as H2O and CO trace gas that is heated and dynamically influenced by outflow activity, whereas CH traces more quiescent material. The four CH absorbers have column densities and turbulent properties consistent with diffuse clouds: two are located near NGC6334, and two are unrelated foreground clouds. Local density and dynamical effects influence the chemical composition of physical components of NGC6334, causing some components to be seen in CH but not in other tracers, and vice versa.

preprint2010arXiv

Herschel/PACS Imaging of Protostars in the HH 1-2 Outflow Complex

We present 70 and 160 micron Herschel science demonstration images of a field in the Orion A molecular cloud that contains the prototypical Herbig-Haro objects HH 1 and 2, obtained with the Photodetector Array Camera and Spectrometer (PACS). These observations demonstrate Herschel's unprecedented ability to study the rich population of protostars in the Orion molecular clouds at the wavelengths where they emit most of their luminosity. The four protostars previously identified by Spitzer 3.6-40 micron imaging and spectroscopy are detected in the 70 micron band, and three are clearly detected at 160 microns. We measure photometry of the protostars in the PACS bands and assemble their spectral energy distributions (SEDs) from 1 to 870 microns with these data, Spitzer spectra and photometry, 2MASS data, and APEX sub-mm data. The SEDs are fit to models generated with radiative transfer codes. From these fits we can constrain the fundamental properties of the protostars. We find luminosities in the range 12-84 L_sun and envelope densities spanning over two orders of magnitude. This implies that the four protostars have a wide range of envelope infall rates and evolutionary states: two have dense, infalling envelopes, while the other two have only residual envelopes. We also show the highly irregular and filamentary structure of the cold dust and gas surrounding the protostars as traced at 160 microns.

preprint2010arXiv

Hier ist wahrhaftig ein Loch im Himmel - The NGC 1999 dark globule is not a globule

The NGC 1999 reflection nebula features a dark patch with a size of ~10,000 AU, which has been interpreted as a small, dense foreground globule and possible site of imminent star formation. We present Herschel PACS far-infrared 70 and 160mum maps, which reveal a flux deficit at the location of the globule. We estimate the globule mass needed to produce such an absorption feature to be a few tenths to a few Msun. Inspired by this Herschel observation, we obtained APEX LABOCA and SABOCA submillimeter continuum maps, and Magellan PANIC near-infrared images of the region. We do not detect a submillimer source at the location of the Herschel flux decrement; furthermore our observations place an upper limit on the mass of the globule of ~2.4x10^-2 Msun. Indeed, the submillimeter maps appear to show a flux depression as well. Furthermore, the near-infrared images detect faint background stars that are less affected by extinction inside the dark patch than in its surroundings. We suggest that the dark patch is in fact a hole or cavity in the material producing the NGC 1999 reflection nebula, excavated by protostellar jets from the V 380 Ori multiple system.

preprint2010arXiv

Imaging of a Transitional Disk Gap in Reflected Light: Indications of Planet Formation Around the Young Solar Analog LkCa 15

We present H- and Ks-band imaging data resolving the gap in the transitional disk around LkCa 15, revealing the surrounding nebulosity. We detect sharp elliptical contours delimiting the nebulosity on the inside as well as the outside, consistent with the shape, size, ellipticity, and orientation of starlight reflected from the far-side disk wall, whereas the near-side wall is shielded from view by the disk's optically thick bulk. We note that forward-scattering of starlight on the near-side disk surface could provide an alternate interpretation of the nebulosity. In either case, this discovery provides confirmation of the disk geometry that has been proposed to explain the spectral energy distributions (SED) of such systems, comprising an optically thick outer disk with an inner truncation radius of ~46 AU enclosing a largely evacuated gap. Our data show an offset of the nebulosity contours along the major axis, likely corresponding to a physical pericenter offset of the disk gap. This reinforces the leading theory that dynamical clearing by at least one orbiting body is the cause of the gap. Based on evolutionary models, our high-contrast imagery imposes an upper limit of 21 Jupiter masses on companions at separations outside of 0.1" and of 13 Jupiter masses outside of 0.2". Thus, we find that a planetary system around LkCa 15 is the most likely explanation for the disk architecture.

preprint2010arXiv

Lucky Imaging survey for southern M dwarf binaries

While M dwarfs are the most abundant stars in the Milky Way, there is still large uncertainty about their basic physical properties (mass, luminosity, radius, etc.) as well as their formation environment. Precise knowledge of multiplicity characteristics and how they change in this transitional mass region, between Sun-like stars on the one side and very low mass stars and brown dwarfs on the other, provide constraints on low mass star and brown dwarf formation. In the largest M dwarf binary survey to date, we search for companions to active, and thus preferentially young, M dwarfs in the solar neighbourhood. We study their binary/multiple properties, such as the multiplicity frequency and distributions of mass ratio and separation, and identify short period visual binaries, for which orbital parameters and hence dynamical mass estimates can be derived in the near future. The observations are carried out in the SDSS i' and z' band using the Lucky Imaging camera AstraLux Sur at the ESO 3.5 m New Technology Telescope. In the first part of the survey, we observed 124 M dwarfs of integrated spectral types M0-M6 and identified 34 new and 17 previously known companions to 44 stars. We derived relative astrometry and component photometry for these systems. More than half of the binaries have separations smaller than 1 arcsec and would have been missed in a simply seeing-limited survey. Correcting our sample for selection effects yields a multiplicity fraction of 32+/-6% for 108 M dwarfs within 52 pc and with angular separations of 0.1-6.0 arcsec, corresponding to projected separation 3-180 AU at median distance 30 pc. Compared to early-type M dwarfs (M>0.3M_Sun), later type (and hence lower mass) M dwarf binaries appear to have closer separations, and more similar masses.

preprint2010arXiv

Ortho-to-para ratio of interstellar heavy water

Despite the low elemental deuterium abundance in the Galaxy, enhanced molecular D/H ratios have been found in the environments of low-mass star forming regions, and in particular the Class 0 protostar IRAS 16293-2422. The CHESS (Chemical HErschel Surveys of Star forming regions) Key Program aims at studying the molecular complexity of the interstellar medium. The high sensitivity and spectral resolution of the HIFI instrument provide a unique opportunity to observe the fundamental 1,1,1 - 0,0,0 transition of the ortho-D2O molecule, inaccessible from the ground, and to determine the ortho-to-para D2O ratio. We have detected the fundamental transition of the ortho-D2O molecule at 607.35 GHz towards IRAS 16293-2422. The line is seen in absorption with a line opacity of 0.62 +/- 0.11 (1 sigma). From the previous ground-based observations of the fundamental 1,1,0 - 1,0,1 transition of para-D2O seen in absorption at 316.80 GHz we estimate a line opacity of 0.26 +/- 0.05 (1 sigma). We show that the observed absorption is caused by the cold gas in the envelope of the protostar. Using these new observations, we estimate for the first time the ortho to para D2O ratio to be lower than 2.6 at a 3 sigma level of uncertainty, to be compared with the thermal equilibrium value of 2:1.

preprint2010arXiv

Tau Boo b: Hunting for reflected starlight

Aims. We attempt to detect starlight reflected from the hot Jupiter orbiting the main-sequence star Tau Boo, in order to determine the albedo of the planetary atmosphere, the orbital inclination of the planetary system and the exact mass of the planetary companion. Methods. We analyze high-precision, high-resolution spectra, collected over two half nights using UVES at the VLT/UT2, by way of data synthesis. We interpret our data using two different atmospheric models for hot Jupiters. Results. Although a weak candidate signal appears near the most probable radial velocity amplitude, its statistical significance is insufficient for us to claim a detection. However, this feature agrees very well with a completely independently obtained result by another research group, which searched for reflected light from Tau Boo b. As a consequence of the non-detection of reflected light, we place upper limits to the planet-to-star flux ratio at the 99.9% significance level. For the most probable orbital inclination around i = 46 degrees, we can limit the relative reflected radiation to be less than 5.7 x 10^-5 for grey albedo. This implies a geometric albedo smaller than 0.40, assuming a planetary radius of 1.2 R_Jup .

preprint2010arXiv

Testing the theory of grain growth and fragmentation by millimeter observations of protoplanetary disks

Context. Observations at sub-millimeter and mm wavelengths will in the near future be able to resolve the radial dependence of the mm spectral slope in circumstellar disks with a resolution of around a few AU at the distance of the closest star-forming regions. Aims. We aim to constrain physical models of grain growth and fragmentation by a large sample of (sub-)mm observations of disks around pre-main sequence stars in the Taurus-Auriga and Ophiuchus star-forming regions. Methods. State-of-the-art coagulation/fragmentation and disk-structure codes are coupled to produce steady-state grain size distributions and to predict the spectral slopes at (sub-)mm wavelengths. Results. This work presents the first calculations predicting the mm spectral slope based on a physical model of grain growth. Our models can quite naturally reproduce the observed mm-slopes, but a simultaneous match to the observed range of flux levels can only be reached by a reduction of the dust mass by a factor of a few up to about 30 while keeping the gas mass of the disk the same. This dust reduction can either be due to radial drift at a reduced rate or during an earlier evolutionary time (otherwise the predicted fluxes would become too low) or due to efficient conversion of dust into larger, unseen bodies.

preprint2010arXiv

The distribution of water in the high-mass star-forming region NGC 6334I

We present observations of twelve rotational transitions of H2O-16, H2O-18, and H2O-17 toward the massive star-forming region NGC 6334 I, carried out with Herschel/HIFI as part of the guaranteed time key program Chemical HErschel Surveys of Star forming regions (CHESS). We analyze these observations to obtain insights into physical processes in this region. We identify three main gas components (hot core, cold foreground, and outflow) in NGC 6334 I and derive the physical conditions in these components. The hot core, identified by the emission in highly excited lines, shows a high excitation temperature of 200 K, whereas water in the foreground component is predominantly in the ortho- and para- ground states. The abundance of water varies between 4 10^-5 (outflow) and 10^-8 (cold foreground gas). This variation is most likely due to the freeze-out of water molecules onto dust grains. The H2O-18/H2O-17 abundance ratio is 3.2, which is consistent with the O-18/O-17 ratio determined from CO isotopologues. The ortho/para ratio in water appears to be relatively low 1.6(1) in the cold, quiescent gas, but close to the equilibrium value of three in the warmer outflow material (2.5(0.8)).

preprint2010arXiv

The methanol lines and hot core of OMC2-FIR4, an intermediate-mass protostar, with Herschel-HIFI

In contrast with numerous studies on the physical and chemical structure of low- and high-mass protostars, much less is known about their intermediate-mass counterparts, a class of objects that could help to elucidate the mechanisms of star formation on both ends of the mass range. We present the first results from a rich HIFI spectral dataset on an intermediate-mass protostar, OMC2-FIR4, obtained in the CHESS (Chemical HErschel SurveyS of star forming regions) key programme. The more than 100 methanol lines detected between 554 and 961 GHz cover a range in upper level energy of 40 to 540 K. Our physical interpretation focusses on the hot core, but likely the cold envelope and shocked regions also play a role in reality, because an analysis of the line profiles suggests the presence of multiple emission components. An upper limit of 10^-6 is placed on the methanol abundance in the hot core, using a population diagram, large-scale source model and other considerations. This value is consistent with abundances previously seen in low-mass hot cores. Furthermore, the highest energy lines at the highest frequencies display asymmetric profiles, which may arise from infall around the hot core.

preprint2010arXiv

The Photodetector Array Camera and Spectrometer (PACS) on the Herschel Space Observatory

The Photodetector Array Camera and Spectrometer (PACS) is one of the three science instruments on ESA's far infrared and submillimetre observatory. It employs two Ge:Ga photoconductor arrays (stressed and unstressed) with 16x25 pixels, each, and two filled silicon bolometer arrays with 16x32 and 32x64 pixels, respectively, to perform integral-field spectroscopy and imaging photometry in the 60-210μ m wavelength regime. In photometry mode, it simultaneously images two bands, 60-85μ m or 85-125μ\m and 125-210μ m, over a field of view of ~1.75'x3.5', with close to Nyquist beam sampling in each band. In spectroscopy mode, it images a field of 47"x47", resolved into 5x5 pixels, with an instantaneous spectral coverage of ~1500km/s and a spectral resolution of ~175km/s. We summarise the design of the instrument, describe observing modes, calibration, and data analysis methods, and present our current assessment of the in-orbit performance of the instrument based on the Performance Verification tests. PACS is fully operational, and the achieved performance is close to or better than the pre-launch predictions.

preprint2010arXiv

Worlds Beyond: A Strategy for the Detection and Characterization of Exoplanets

This is a scientific strategy for the detection and characterization of extrasolar planets; that is, planets orbiting other stars. As such, it maps out over a 15-year horizon the techniques and capabilities required to detect and measure the properties of planets as small as Earth around stars as large as our own Sun. It shows how the technology pieces and their development fit together to achieve the primary goal of the strategy: if planets like Earth exist around stars within some tens of light years of our own Solar System, those planets will be found and their basic properties characterized. Essential to this strategy is not only the search for and examination of individual planets, but also a knowledge of the arrangement, or architecture, of planetary systems around as large a number of stars as possible; this is the second goal of the strategy. The final goal of the strategy is the study of disks around stars, important both to understand the implications of the variety of exoplanet systems for planet formation, and to determine how many nearby stars have environments around them clean enough of debris that planets may be sought and, if found, characterized.

preprint2009arXiv

ATLASGAL - The APEX Telescope Large Area Survey of the Galaxy at 870 microns

(Abridged) Studying continuum emission from interstellar dust is essential to locating and characterizing the highest density regions in the interstellar medium. In particular, the early stages of massive star formation remain poorly understood. Our goal is to produce a large-scale, systematic database of massive pre- and proto-stellar clumps in the Galaxy, to understand how and under what conditions star formation takes place. A well characterized sample of star-forming sites will deliver an evolutionary sequence and a mass function of high-mass, star-forming clumps. This systematic survey at submm wavelengths also represents a preparatory work for Herschel and ALMA. The APEX telescope is ideally located to observe the inner Milky Way. The Large APEX Bolometer Camera (LABOCA) is a 295-element bolometer array observing at 870 microns, with a beam size of 19". Taking advantage of its large field of view (11.4') and excellent sensitivity, we started an unbiased survey of the Galactic Plane, with a noise level of 50-70 mJy/beam: the APEX Telescope Large Area Survey of the Galaxy (ATLASGAL). As a first step, we covered 95 sq. deg. These data reveal 6000 compact sources brighter than 0.25 Jy, as well as extended structures, many of them filamentary. About two thirds of the compact sources have no bright infrared counterpart, and some of them are likely to correspond to the precursors of (high-mass) proto-stars or proto-clusters. Other compact sources harbor hot cores, compact HII regions or young embedded clusters. Assuming a typical distance of 5 kpc, most sources are clumps smaller than 1 pc with masses from a few 10 to a few 100 M_sun. In this introductory paper, we show preliminary results from these ongoing observations, and discuss the perspectives of the survey.

preprint2009arXiv

Brown dwarfs and very low mass stars in the Praesepe open cluster: a dynamically unevolved mass function?

[Abridged] In this paper, we present the results of a photometric survey to identify low mass and brown dwarf members of the old open cluster Praesepe (age of 590[+150][-120]Myr and distance of 190[+6.0][-5.8]pc) and use this to infer its mass function which we compare with that of other clusters. We have performed an optical (Ic-band) and near-infrared (J and Ks-band) photometric survey of Praesepe with a spatial coverage of 3.1deg^2. With 5sigma detection limits of Ic=23.4 and J=20.0, our survey is sensitive to objects with masses from about 0.6 to 0.05Msol. The mass function of Praesepe rises from 0.6Msol down to 0.1Msol and then turns-over at ~0.1Msol. The rise observed is in agreement with the mass function derived by previous studies, including a survey based on proper motion and photometry. Comparing our mass function with that for another open cluster with a similar age, the Hyades (age ~ 600Myr), we see a significant difference. Possible reasons are that dynamical evaporation has not influenced the Hyades and Praesepe in the same way, or that the clusters did not have the same initial mass function, or that dynamical interactions have modified the evolution of one or both clusters. Although a difference in the binary fractions of the clusters could cause the observed (i.e. system) mass functions to differ, measurements in the literature give no evidence for a significant difference in the binary fractions of the two clusters. Of our cluster candidates, six have masses predicted to be equal to or below the stellar/substellar boundary at 0.072Msol.

preprint2009arXiv

Candidate free-floating super-Jupiters in the young sigma Orionis open cluster

Free-floating substellar candidates with estimated theoretical masses of as low as ~5 Jupiter masses have been found in the ~3 Myr old sigma Orionis open cluster. As the overlap with the planetary mass domain increases, the question of how these objects form becomes important. The determination of their number density and whether a mass cut-off limit exists is crucial to understanding their formation. We propose to search for objects of yet lower masses in the cluster and determine the shape of the mass function at low mass. Using new- and (re-analysed) published IZJHKs[3.6]-[8.0]-band data of an area of 840 arcmin2, we performed a search for LT-type cluster member candidates in the magnitude range J=19.5-21.5 mag, based on their expected magnitudes and colours. Besides recovering the T type object S Ori 70 and two other known objects, we find three new cluster member candidates, S Ori 72-74, with J=21 mag and within 12 arcmin of the cluster centre. They have theoretical masses of 4 (-2,+3) M_Jup and are among the least massive free-floating objects detected by direct imaging outside the Solar System. The photometry in archival Spitzer [3.6]-[5.8]-band images infers that S Ori 72 is an L/T transition candidate and S Ori 73 a T-type candidate, following the expected cluster sequence in the mid-infrared. Finally, the L-type candidate S Ori 74 with lower quality photometry is located at 11.8 arcsec (~4250 AU) of a stellar member of sigma Orionis and could be a companion. After contaminant correction in the area complete to J=21.1 mag, we estimate that there remain between zero and two cluster members in the mass interval 6-4 M_Jup. Our result suggests a possible turnover in the substellar mass spectrum below ~6 Jupiter masses, which could be investigated further by wider and deeper photometric surveys.

preprint2009arXiv

Probing the evolution of molecular cloud structure: From quiescence to birth

Aims: We derive the probability density functions (PDFs) of column density for a complete sample of prominent molecular cloud complexes closer than 200 pc. Methods: We derive near-infrared dust extinction maps for 23 molecular cloud complexes, using the "nicest" colour excess mapping technique and data from the 2MASS archive. The extinction maps are then used to examine the column density PDFs in the clouds. Results: The column density PDFs of most molecular clouds are well-fitted by log-normal functions at low column densities (0.5 mag < A_v < 3-5 mag). However, at higher column densities prominent, power-law-like wings are common. In particular, we identify a trend among the PDFs: active star-forming clouds always have prominent non-log-normal wings. In contrast, clouds without active star formation resemble log-normals over the whole observed column density range, or show only low excess of higher column densities. This trend is also reflected in the cumulative PDFs, showing that the fraction of high column density material is significantly larger in star-forming clouds. These observations are in agreement with an evolutionary trend where turbulent motions are the main cloud-shaping mechanism for quiescent clouds, but the density enhancements induced by them quickly become dominated by gravity (and other mechanisms) which is strongly reflected by the shape of the column density PDFs. The dominant role of the turbulence is restricted to the very early stages of molecular cloud evolution, comparable to the onset of active star formation in the clouds.

preprint2008arXiv

Binarity at the L/T brown dwarf transition

Current atmospheric models cannot reproduce some of the characteristics of the transition between the L dwarfs with cloudy atmospheres and the T dwarfs with dust-depleted photospheres. It has been proposed that a majority of the L/T transition brown dwarfs could actually be a combinaison of a cloudy L dwarf and a clear T dwarf. Indeed binarity seems to occur more frequently among L/T transition brown dwarfs. We aim to refine the statistical significance of the seemingly higher frequency of binaries. Co-eval binaries would also be interesting test-beds for evolutionary models. We obtained high-resolution imaging for six mid-L to late-T dwarfs, with photometric distances between 8 and 33pc, using the adaptive optics systems NACO at the VLT, and the Lick system, both with the laser guide star. We resolve none of our targets. Combining our data with published results, we obtain a frequency of resolved L/T transition brown dwarfs of (31+21-15)%, compared to (21+10-7)% and (14+14-7)% for mid-L and T dwarfs (90% of confidence level). These fractions do not significantly support, nor contradict, the hypothesis of a larger binary fraction in the L/T transition. None of our targets has companions with effective temperatures as low as 360-1000K at separations larger than 0.5".

preprint2008arXiv

Probing dust grain evolution in IM Lupi's circumstellar disc. Multi-wavelength observations and modelling of the dust disc

We present a panchromatic study, involving a multiple technique approach, of the circumstellar disc surrounding the T Tauri star IM Lupi (Sz 82). We have undertaken a comprehensive observational study of IM Lupi using photometry, spectroscopy, millimetre interferometry and multi-wavelength imaging. For the first time, the disc is resolved from optical and near-infrared wavelengths in scattered light, to the millimetre regime in thermal emission. Our data-set, in conjunction with existing photometric data, provides an extensive coverage of the spectral energy distribution, including a detailed spectrum of the silicate emission bands. We have performed a simultaneous modelling of the various observations, using the radiative transfer code MCFOST, and analysed a grid of models over a large fraction of the parameter space via Bayesian inference. We have constructed a model that can reproduce all of the observations of the disc. Our analysis illustrates the importance of combining a wide range of observations in order to fully constrain the disc model, with each observation providing a strong constraint only on some aspects of the disc structure and dust content. Quantitative evidence of dust evolution in the disc is obtained: grain growth up to millimetre-sized particles, vertical stratification of dust grains with micrometric grains close to the disc surface and larger grains which have settled towards the disc midplane, and possibly the formation of fluffy aggregates and/or ice mantles around grains.

preprint2000arXiv

Realizing 3D Spectral Imaging in the Far-Infrared: FIFI LS

We present a progress report on the design and construction of the Field-Imaging Far-Infrared Line Spectrometer (FIFI LS) for the SOFIA airborne observatory. The design of the instrument is driven by the goal of maximizing observing efficiency, especially for observations of faint, extragalactic objects. Thus, FIFI LS utilizes an integral field technique that uses slicer mirrors to optically re-arrange the two-dimensional field into a single slit for a long slit spectrometer. Effectively, a 5x5 pixel spatial field of view is imaged to a 25x1 pixel slit and dispersed to a 25x16 pixel, two-dimensional detector array, providing diffraction-limited spatial and spectral multiplexing. In this manner, the instrument employs two parallel, medium resolution (R~2000) grating spectrometers for simultaneous observations in two bands: a short wavelength band (42 to 110 micron) and a long wavelength band (110 to 210 micron). Overall, for each of the 25 spatial pixels, the instrument can cover a velocity range of ~1500 km/s around selected far-infrared lines with an estimated sensitivity of 2x10^-15 W Hz^1/2 per pixel. This arrangement provides good spectral coverage with high responsivity. ***This paper does not include Figures due to astro-ph size limitations. Please download entire file at http://fifi-ls.mpe-garching.mpg.de/fifils.ps.gz ****