Researcher profile

Th. Henning

Th. Henning contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
58works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

58 published item(s)

preprint2026arXiv

Characterization of two new transiting sub-Neptunes and a terrestrial planet around M-dwarf hosts

We report the confirmation of three transiting exoplanets orbiting TOI-1243 (LSPM~J0902+7138), TOI-4529 (G~2--21), and TOI-5388 (Wolf~346) that were initially detected by TESS through ground-based photometry and radial velocity follow-up measurements with CARMENES. The planets present short orbital periods of $4.65$, $5.88$, and $2.59$ days, and they orbit early-M dwarfs (M2.0V, M1.5V, and M3.0V, respectively). We were able to precisely determine the radius of all three planets with a precision of $< 7\, \%$, the mass of TOI-1243 b with a precision of $19\, \%$, and upper mass limits for TOI-4529 b and TOI-5388 b. The radius of TOI-1243 b is $2.33\pm0.12\, {R_{\oplus}}$, its mass is $7.7 \pm 1.5\,{M_{\oplus}}$, and the mean density is $0.61 \pm 0.15 \, {ρ_\oplus}$. The radius of TOI-4529 b is $1.77 ^{+0.09}_{-0.08} \, {R_{\oplus}}$, the $3 σ$ upper mass limit is $4.9 \, {M_{\oplus}}$, and the $3 σ$ upper density limit is $0.88\, {ρ_\oplus}$. The third planet, TOI-5388 b, is Earth-sized with a radius of $0.99 ^{+0.07}_{-0.06} \, {R_{\oplus}}$, a $3 σ$ upper mass limit of $2.2 \, {M_{\oplus}}$, and a $3 σ$ upper density limit of $2.2\, {ρ_\oplus}$. While TOI-5388 b is most probably rocky, given its Earth-like radius, TOI-1243 b and TOI-4529 b are located in a highly degenerate region in the mass-radius space. TOI-4529 b appears to lean toward a water-world composition. TOI-1243 b has enough mass to host a significant H-He envelope, although a water-world and pure rocky compositions are also consistent with the data. Our analysis indicates that future atmospheric observations using JWST can aid in determining their real composition. The sample of small planets around M dwarfs is widely used to understand planet formation and composition theories, and our study adds three planets to this sample.

preprint2026arXiv

JWST/MIRI coronagraphic search for planets in systems with gapped exoKuiper belts and proper motion anomalies

Over the past decade, ALMA has uncovered a range of substructures within exoKuiper belts, pointing to a population of undetected planets. With JWST&#39;s sensitivity, we now have the opportunity to identify these planets thought to be responsible for the observed substructures in debris discs. We present Cycle 1 JWST/MIRI 11.4 μm coronagraphic observations of three exoKuiper belts that exhibit gaps in their radial structures: HD 92945, HD 107146, and HD 206893, to determine whether planets are responsible for carving these structures, as seen in our Solar System with the gas giants. We reduce the JWST/MIRI data using spaceKLIP, and introduce new routines to mitigate the Brighter-Fatter effect and persistence. We do not detect any planet candidates, and all detected objects in the field-of-view are consistent with background stars or galaxies. However, by combining JWST mass limits, archival observational constraints, and astrometric accelerations, we rule out a significant portion of planet parameter space, placing tight constraints on the planets possibly responsible for these gaps. To interpret these results, we explore multiple gap-carving scenarios in discs, either massless or with non-zero mass, including clearing by in-situ planet(s), as well as shaping by inner planets through mean-motion or secular apsidal resonances. Finally, we conclude that the planets causing the proper motion anomaly in these systems must reside within the inner 20 au.

preprint2026arXiv

The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) I: Motivation, sample, data reduction, and results overview

The outer regions of planetary systems host dusty debris discs analogous to the Kuiper belt (exoKuiper belts), which provide crucial constraints on planet formation and evolution processes. ALMA dust observations have revealed a great diversity, and that some belts contain CO gas, whose origin and implications are uncertain. Most of this progress, however, has been limited by low-resolution observations. We conducted the first ALMA large programme dedicated to debris discs: the ALMA survey to Resolve exoKuiper belt Substructures (ARKS). We selected the 24 most promising belts to constrain their detailed radial and vertical structure, and to characterise the gas content. We constrained the radial and vertical distribution of dust, as well as the presence of asymmetries. For a subset of six belts with CO gas, we constrained the gas distribution and kinematics. To interpret these observations, we used a wide range of dynamical models. The first ARKS results are presented as a series of ten papers. We discovered that up to 33% of our sample exhibits multiple dusty rings. For highly inclined belts, we found that non-Gaussian vertical distributions are common and are indicative of multiple dynamical populations. We also found that 10 of the 24 belts present asymmetries. We find that the CO gas is radially broader than the dust, but this could be an effect of optical depth. At least one system shows non-Keplerian kinematics due to strong pressure gradients, which may have triggered a vortex that trapped dust in an arc. Finally, we find evidence that the micron-sized grains may be affected by gas drag in gas rich systems. ARKS has revealed a great diversity of structures in exoKuiper belts that may arise when they are formed in protoplanetary discs or subsequently via interactions with planets and/or gas. We encourage the community to explore the reduced data and data products.

preprint2026arXiv

The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) V: Comparison between scattered light and thermal emission

Debris discs are analogues to our own Kuiper belt around main-sequence stars and are therefore referred to as exoKuiper belts. They have been resolved at high angular resolution at wavelengths spanning the optical to the submillimetre-millimetre regime. Short wavelengths probe the light scattered by such discs, which is dominated by micron-sized dust particles, while millimetre wavelengths probe the thermal emission of millimetre-sized particles. Determining differences in the dust distribution between millimetre- and micron-sized dust is fundamental to revealing the dynamical processes affecting the dust in debris discs. We aim to compare the scattered light from the discs of the ALMA survey to Resolve exoKuiper belt Substructures (ARKS) with the thermal emission probed by ALMA. We focus on the radial distribution of the dust. We used high-contrast scattered light observations obtained with VLT/SPHERE, GPI, and the HST to uniformly study the dust distribution in those systems and compare it to the dust distribution extracted from the ALMA observations carried out in the course of the ARKS project. We also set constraints on the presence of planets by using these high-contrast images combined with exoplanet evolutionary models. 15 of the 24 discs comprising the ARKS sample are detected in scattered light, with TYC9340-437-1 being imaged for the first time at near-infrared wavelengths. For 6 of those 15 discs, the dust surface density seen in scattered light peaks farther out compared to that observed with ALMA. These 6 discs except one are known to also host cold CO gas. Conversely, the systems without significant offsets are not known to host gas, except one. This observational study suggests that the presence of gas in debris discs may affect the small and large grains differently, pushing the small dust to greater distances where the gas is less abundant.

preprint2026arXiv

The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) VII: Optically thick gas with broad CO gaussian local line profiles in the HD 121617 disc

CO gas has been detected in $\sim$20 debris discs. We present ALMA observations of the CO-rich HD 121617 debris disc from the ARKS survey. Using high-resolution Band 7 observations of $^{12}CO \ J=3-2$, we analyse local CO line profiles to investigate optical depth, CO mass, and temperature. Spectra are aligned and stacked in concentric annuli to produce local line profiles. The resulting profiles are Gaussian-shaped and broadened by Keplerian shear. The line profiles are modelled using both a simplified toy model and a RADMC-3D model including projection effects and Keplerian shear. Fitting the RADMC-3D model to the $^{13}$CO data, we find that an optically thick model with a temperature of 38 K and a CO mass of $2 \times 10^{-3}$ M$_{\oplus}$ reproduces the observations. The model reproduces the enhanced emission at orbital azimuths of $\sim \pm45^{\circ}$ and $\pm135^{\circ}$, forming an X-shaped structure in the velocity-integrated intensity map, as well as the broader $^{12}$CO linewidth relative to $^{13}$CO. Scaling the model by the ISM abundance ratio ($\sim$77) also reproduces the $^{12}$CO data, though high optical depths and model assumptions limit mass constraints. We find that azimuthally averaged local line profiles appear Gaussian regardless of optical depth, cautioning against their use for distinguishing optically thin and thick emission. We constrain the mean molecular weight to $12.6_{-1.1}^{+1.3}$, dependent on model assumptions. Our $^{13}$CO results suggest that C$^{18}$O may also be optically thick in CO-rich debris discs and that the mean molecular weight is significantly higher than if H$_2$ were the dominant gas species, suggesting a non-primordial composition.

preprint2026arXiv

The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) VIII: A dust arc and non-Keplerian gas kinematics in HD 121617

ExoKuiper belts around young A-type stars often host CO gas, whose origin is still unclear. The ALMA survey to Resolve exoKuiper belt Substructures (ARKS) includes 6 of these gas-bearing belts, to characterise their dust and gas distributions and investigate the gas origin. As part of ARKS, we observed the gas-rich system HD121617 and discovered an arc of enhanced dust density. In this paper, we analyse in detail the dust and gas distributions and the gas kinematics of this system. We extracted radial and azimuthal profiles of the dust (in the millimetre and near-infrared) and gas emission ($^{12}$CO and $^{13}$CO) from reconstructed images. To constrain the morphology of the arc, we fitted an asymmetric model to the dust emission. To characterise the gas kinematics, we fitted a Keplerian model to the velocity map and extracted the azimuthal velocity profile by deprojecting the data. We find that the dust arc is narrow (1-5 au wide at a radius of 75 au), azimuthally extended, and asymmetric; the emission is more azimuthally compact in the direction of the system&#39;s rotation, and represents 13% of the total dust mass (0.2$M_\oplus$). The arc is much less pronounced or absent for small grains and gas. Finally, we find strong non-Keplerian azimuthal velocities at the inner and outer wings of the ring, as was expected due to strong pressure gradients. The dust arc resembles the asymmetries found in protoplanetary discs, often interpreted as the result of dust trapping in vortices. If the gas disc mass is high enough ($\gtrsim20M_\oplus$, requiring a primordial gas origin), both the radial confinement of the ring and the azimuthal arc may result from dust grains responding to gas drag. Alternatively, it could result from planet-disc interactions via mean motion resonances. Further studies should test these hypotheses and may provide a dynamical gas mass estimate in this CO-rich exoKuiper belt.

preprint2026arXiv

The ALMA survey to Resolve exoKuiper belt Substructures (ARKS): VI. Asymmetries and offsets

Asymmetries in debris discs provide unique clues to understand the evolution and architecture of planetary systems.** The aim of the ALMA survey to Resolve exoKuiper belt Substructures (ARKS) is to expand our understanding of radial and vertical dust structures, as well as gas distributions and kinematics, in debris discs.** Here, in ARKS~VI, we present a systematic analysis of the asymmetries and stellocentric offsets present in the ALMA continuum data for the ARKS survey. Our aims are to identify asymmetries in debris disc dust distributions, quantify debris disc asymmetry properties, and discuss the potential origins of debris disc asymmetries.** We utilised empirical methods to identify emission asymmetries** and the presence of offset emission between disc centres and the locations of the host stars, via an analysis of their calibration procedures and disc properties. We associated observational asymmetry types** and plausible physical classes** associated with each source. We show that there are ten systems, almost half of the ARKS sample, that host either a continuum emission asymmetry or offset emission. Three systems host offsets (HD15115, HD32297, and HD109573 (HR4796)), four host azimuthal asymmetries (HD9672 (49Ceti), HD92945, HD107146, and HD121617), two host an asymmetry in their major axis (HD10647 (q$^1$ Eri), and HD39060 ($β$ Pic)), and one hosts an asymmetry in their minor axis (HD61005). We attribute the offset asymmetries to non-zero eccentricities, and three of the azimuthal asymmetries to arcs. The presence of an asymmetry or offset in the ARKS sample appears to be correlated with the fractional luminosity of cold dust.** Conclusions: This study demonstrates that debris disc asymmetries in the ARKS sample are common, and plausibly so in the wider population of debris discs at (sub)-millimetre wavelengths.** ** = ABRIDGED FOR ARXIV: FULL ABSTRACT IN PAPER

preprint2026arXiv

The ALMA survey to Resolve exoKuiper belt Substructures (ARKS). X. Interpreting the peculiar dust rings around HD 131835

Dusty discs detected around main-sequence stars are thought to be signs of planetesimal belts in which the dust distribution is shaped by collisional and dynamical processes, including interactions with gas if present. The debris disc around the young A-type star HD 131835 is composed of two dust rings at ~65 au and ~100 au, a third unconstrained innermost component, and a gaseous component centred at ~65 au. New ALMA observations show that the inner of the two dust rings is brighter than the outer one, in contrast with previous observations in scattered light. We explore two scenarios that could explain these observations: the two dust rings might represent distinct planetesimal belts with different collisional properties, or only the inner ring might contain planetesimals while the outer ring consists entirely of dust that has migrated outwards due to gas drag. To explore the first scenario, we employed a state-of-the-art collisional evolution code. To test the second scenario, we used a simple dynamical model of dust grain evolution in an optically thin gaseous disc. Collisional models of two planetesimal belts cannot fully reproduce the observations by only varying their dynamical excitation, and matching the data through a different material strength requires an extreme difference in dust composition. The gas-driven scenario can reproduce the location of the outer ring and the brightness ratio of the two rings from scattered light observations, but the resulting outer ring is too faint overall in both scattered light and sub-millimetre emission. The dust rings in HD 131835 could be produced from two planetesimal belts, although how these belts would attain the required extremely different properties needs to be explained. The dust-gas interaction is a plausible alternative explanation and deserves further study using a more comprehensive model.

preprint2026arXiv

The ExoGRAVITY survey: A K-band spectral library of giant exoplanet and brown dwarf companions

Direct observations of exoplanet and brown dwarf companions with near-infrared interferometry, first enabled by the dual-field mode of VLTI/GRAVITY, provide unique measurements of the objects&#39; orbital motions and atmospheric compositions. Here, we compile a homogeneous library of all exoplanet and brown dwarf K-band spectra observed by GRAVITY thus far. We re-reduced all the available GRAVITY dual-field high-contrast data and, where companions are detected, extract their ~2.0-2.4 $μ$m K-band contrast spectra. We then derived stellar model atmospheres for all employed flux references, which we used to convert the companion contrast into companion flux spectra. Solely from the resulting GRAVITY spectra, we extracted spectral types, spectral indices, and bulk physical properties for all companions. Finally, and with the help of age constraints from the literature, we also derived isochronal masses for most companions using evolutionary models. The resulting library contains R ~ 500 GRAVITY spectra of 39 substellar companions from late M to late T spectral types, including the entire L-T transition. Throughout this transition, a shift from CO-dominated late M- and L-type dwarfs to CH4-dominated T-type dwarfs can be observed in the K-band. The GRAVITY spectra alone constrain the objects&#39; bolometric luminosity to typically within $\pm$0.15 dex. The derived isochronal masses agree with dynamical masses from the literature where available, except for HD 4113 c for which we confirm its previously reported potential underluminosity. Medium-resolution spectroscopy of substellar companions with GRAVITY provides insight into the carbon chemistry and the cloudiness of these objects&#39; atmospheres. It also constrains these objects&#39; bolometric luminosities which can yield measurements of their formation entropy if combined with dynamical masses, for instance from Gaia and GRAVITY astrometry.

preprint2025arXiv

JWST/MIRI observations of the young TWA 27 system: hydrocarbon disk chemistry, silicate clouds, evidence for a CPD

The Mid-Infrared Instrument (MIRI) on the James Webb Space Telescope (JWST) enables the characterisation of young self-luminous gas giants at previously inaccessible wavelengths, revealing physical processes in gas, dust, and clouds. We characterise the young planetary system TWA 27 (2M1207) in the mid-infrared (MIR), studying the atmosphere and disk spectra of the M9 brown dwarf TWA 27A and its L6 planetary-mass companion TWA 27b. We obtained data with the MIRI Medium Resolution Spectrometer (MRS) from 4.9 to 20 um, and MIRI Imaging in the F1000W and F1500W filters. We applied high-contrast imaging methods to extract the companion&#39;s spectral energy distribution up to 15 um at 0.78 arcsec separation and a contrast of 60. Combining these with published JWST/NIRSpec spectra, we analysed the 1-20 um range using self-consistent atmospheric grids and 0D slab models for molecular disk emission. The atmosphere of TWA 27A is well fitted by a BT-SETTL model with Teff 2780 K, log g 4.3, plus a 740 K blackbody for the inner disk rim. The disk shows at least 11 organic molecules, with no water or silicate dust emission detected. The atmosphere of TWA 27b is matched by a Teff 1400 K low-gravity model with extinction, best fit by the ExoREM grid. MIRI spectra and photometry for TWA 27b reveal a silicate cloud absorption feature between 8-10 um and significant (>5 sigma) infrared excess at 15 um consistent with circumplanetary disk emission. These MIRI observations provide new insights into TWA 27, revealing diverse features to study the formation and evolution of circumplanetary disks and young dusty atmospheres.

preprint2024arXiv

The elusive atmosphere of WASP-12 b / High-resolution transmission spectroscopy with CARMENES

To date, the hot Jupiter WASP-12 b has been the only planet with confirmed orbital decay. The late F-type host star has been hypothesized to be surrounded by a large structure of circumstellar material evaporated from the planet. We obtained two high-resolution spectral transit time series with CARMENES and extensively searched for absorption signals by the atomic species Na, H, Ca, and He using transmission spectroscopy, thereby covering the He I triplet with high resolution for the first time. We apply SYSREM for atomic line transmission spectroscopy, introduce the technique of signal protection to improve the results for individual absorption lines, and compare the outcomes to those of established methods. No transmission signals were detected and the most stringent upper limits as of yet were derived for the individual indicators. Nonetheless, we found variation in the stellar Halpha and He I lines, the origin of which remains uncertain but is unlikely to be activity. To constrain the enigmatic activity state of WASP-12, we analyzed XMM-Newton X-ray data and found the star to be moderately active at most. We deduced an upper limit for the X-ray luminosity and the irradiating X-ray and extreme ultraviolet (XUV) flux of WASP-12 b. Based on the XUV flux upper limit and the lack of the He I signal, our hydrodynamic models slightly favor a moderately irradiated planet with a thermospheric temperature of <= 12,000 K, and a conservative upper limit of <= 4e12 g/s on the mass-loss rate. Our study does not provide evidence for an extended planetary atmosphere or absorption by circumstellar material close to the planetary orbit.

preprint2022arXiv

A pathway to peptides in space through the condensation of atomic carbon

Organic molecules are widely present in the dense interstellar medium, and many have been synthesized in the laboratory on Earth under the conditions typical for an interstellar environment. Until now, however, only relatively small molecules of biological interest have been demonstrated to form experimentally under typical space conditions. Here we prove experimentally that the condensation of carbon atoms on the surface of cold solid particles (cosmic dust) leads to the formation of isomeric polyglycine monomers (aminoketene molecules). Following encounters between aminoketene molecules, they polymerize to produce peptides of different lengths. The chemistry involves three of the most abundant species (CO, C and NH$_3$) present in star-forming molecular clouds, and proceeds via a novel pathway that skips the stage of amino acid formation in protein synthesis. The process is efficient, even at low temperatures, without irradiation or the presence of water. The delivery of biopolymers formed by this chemistry to rocky planets in the habitable zone might be an important element in the origins of life.

preprint2022arXiv

A tentative detection of He I in the atmosphere of GJ 1214b

The $\ion{He}{I}$ 10833\,$Å$ triplet is a powerful tool for characterising the upper atmosphere of exoplanets and tracing possible mass loss. Here, we analysed one transit of GJ\,1214\,b observed with the CARMENES high-resolution spectrograph to study its atmosphere via transmission spectroscopy around the $\ion{He}{I}$ triplet. Although previous studies using lower resolution instruments have reported non-detections of $\ion{He}{I}$ in the atmosphere of GJ\,1214\,b, we report here the first potential detection. We reconcile the conflicting results arguing that previous transit observations did not present good opportunities for the detection of $\ion{He}{I}$, due to telluric H$_2$O absorption and OH emission contamination.We simulated those earlier observations, and show evidence that the planetary signal was contaminated. From our single non-telluric-contaminated transit, we determined an excess absorption of 2.10$^{+0.45}_{-0.50}$\,\% (4.6\,$σ$) with a full width at half maximum (FWHM) of 1.30$^{+0.30}_{-0.25}$\,Å. The detection of \ion{He}{I} is statistically significant at the 4.6\,$σ$ level, but repeatability of the detection could not be confirmed due to the availability of only one transit. By applying a hydrodynamical model and assuming an H/He composition of 98/2, we found that GJ\,1214\,b would undergo hydrodynamic escape in the photon-limited regime, losing its primary atmosphere with a mass-loss rate of (1.5--18)\,$\times$\,10$^{10}$\,g\,s$^{-1}$ and an outflow temperature in the range of 2900--4400\,K. Our $\ion{He}{I}$ excess absorption is the first tentative detection of a chemical species in the atmosphere of this benchmark sub-Neptune planet.

preprint2022arXiv

Detection of CO emission lines in the dayside atmospheres of WASP-33b and WASP-189b with GIANO

Ultra-hot Jupiters (UHJs) are expected to possess temperature inversion layers in their dayside atmospheres. Recent thermal emission observations have discovered several atomic and molecular species along with temperature inversions in UHJs. We observed the thermal emission spectra of two UHJs (WASP-33b and WASP-189b) with the GIANO-B high-resolution near-infrared spectrograph. Using the cross-correlation technique, we detected carbon monoxide (CO) in the dayside atmospheres of both planets. The detected CO lines are in emission, which agrees with previous discoveries of iron emission lines and temperature inversions in the two planets. This is the first detection of CO lines in emission with high-resolution spectroscopy. Further retrieval work combining the CO lines with other spectral features will enable a comprehensive understanding of the atmospheric properties such as temperature structures and C/O ratios. The detected CO and iron emission lines of WASP-189b have redshifted radial velocities of several km/s, which likely originate from a dayside to nightside wind in its atmosphere. Such a redshifted velocity has not been detected for the emission lines of WASP-33b, suggesting that the atmospheric circulation patterns of the two UHJs may be different.

preprint2022arXiv

Detection of iron emission lines and a temperature inversion on the dayside of the ultra-hot Jupiter KELT-20b

Ultra-hot Jupiters (UHJs) are gas giants with very high equilibrium temperatures. In recent years, multiple chemical species, including various atoms and ions, have been discovered in their atmospheres. Most of these observations have been performed with transmission spectroscopy, although UHJs are also ideal targets for emission spectroscopy due to their strong thermal radiation. We present high-resolution thermal emission spectroscopy of the transiting UHJ KELT-20b/MASCARA-2b. The observation was performed with the CARMENES spectrograph at orbital phases before and after the secondary eclipse. We detected atomic Fe using the cross-correlation technique. The detected Fe lines are in emission, which unambiguously indicates a temperature inversion on the dayside hemisphere. We furthermore retrieved the temperature structure with the detected Fe lines. The result shows that the atmosphere has a strong temperature inversion with a temperature of $4900\pm{700}$ K and a pressure of $10^{-4.8_{-1.1}^{+1.0}}$ bar at the upper layer of the inversion. A joint retrieval of the CARMENES data and the TESS secondary eclipse data returns a temperature of $2550_{-250}^{+150}$ K and a pressure of $10^{-1.5_{-0.6}^{+0.7}}$ bar at the lower layer of the temperature inversion. The detection of such a strong temperature inversion is consistent with theoretical simulations that predict an inversion layer on the dayside of UHJs. The joint retrieval of the CARMENES and TESS data demonstrates the power of combing high-resolution emission spectroscopy with secondary eclipse photometry in characterizing atmospheric temperature structures.

preprint2022arXiv

First Light for GRAVITY Wide: Large Separation Fringe Tracking for the Very Large Telescope Interferometer

GRAVITY+ is the upgrade of GRAVITY and the Very Large Telescope Interferometer (VLTI) with wide-separation fringe tracking, new adaptive optics, and laser guide stars on all four 8~m Unit Telescopes (UTs), for ever fainter, all-sky, high contrast, milliarcsecond interferometry. Here we present the design and first results of the first phase of GRAVITY+, called GRAVITY Wide. GRAVITY Wide combines the dual-beam capabilities of the VLTI and the GRAVITY instrument to increase the maximum separation between the science target and the reference star from 2 arcseconds with the 8 m UTs up to several 10 arcseconds, limited only by the Earth&#39;s turbulent atmosphere. This increases the sky-coverage of GRAVITY by two orders of magnitude, opening up milliarcsecond resolution observations of faint objects, and in particular the extragalactic sky. The first observations in 2019 - 2022 include first infrared interferometry of two redshift $z\sim2$ quasars, interferometric imaging on the binary system HD 105913A, and repeated observations of multiple star systems in the Orion Trapezium Cluster. We find the coherence loss between the science object and fringe-tracking reference star well described by the turbulence of the Earth&#39;s atmosphere. We confirm that the larger apertures of the UTs result in higher visibilities for a given separation due to larger overlap of the projected pupils on sky and give predictions for visibility loss as a function of separation to be used for future planning.

preprint2022arXiv

Formation of the simplest amide in molecular clouds: formamide (NH$_{2}$CHO) and its derivatives in H$_2$O-rich and CO-rich interstellar ice analogs upon VUV irradiation

The astronomical detection of formamide (NH$_2$CHO) toward various star-forming regions and in cometary material implies that the simplest amide might have an early origin in dark molecular clouds at low temperatures. Laboratory studies have proven the efficient NH$_2$CHO formation in interstellar CO:NH$_3$ ice analogs upon energetic processing. However, it is still under debate, whether the proposed radical-radical recombination reactions forming complex organic molecules remain valid in an abundant H$_2$O environment. The aim of this work was to investigate the formation of NH$_2$CHO in H$_2$O- and CO-rich ices under conditions prevailing in molecular clouds. Therefore, different ice mixtures composed of H$_2$O:CO:NH$_3$ (10:5:1), CO:NH$_3$ (4:1), and CO:NH$_3$ (0.6:1) were exposed to vacuum ultraviolet photons in an ultra-high vacuum chamber at 10 K. Fourier-transform infrared spectroscopy was utilized to monitor in situ the initial and newly formed species as a function of photon fluence. The infrared spectral identifications are complementarily secured by a temperature-programmed desorption experiment combined with a quadrupole mass spectrometer. The energetic processing of CO:NH$_3$ ice mixtures mainly leads to the NH$_2$CHO formation, along with its chemical derivatives such as isocyanic acid (HNCO) and cyanate ion (OCN$^-$). The formation kinetics of NH$_2$CHO shows an explicit dependency on ice ratios and compositions; the highest yield is found in H$_2$O-rich ice. The astronomical relevance of the resulting reaction network is discussed.

preprint2022arXiv

Magnetism, rotation, and nonthermal emission in cool stars -- Average magnetic field measurements in 292 M dwarfs

Stellar dynamos generate magnetic fields that are of fundamental importance to the variability and evolution of Sun-like and low-mass stars, and for the development of their planetary systems. We report measurements of surface-average magnetic fields in 292 M dwarfs from a comparison with radiative transfer calculations; for 260 of them, this is the first measurement of this kind. Our data were obtained from more than 15,000 high-resolution spectra taken during the CARMENES project. They reveal a relation between average field strength, <B>, and Rossby number, $Ro$, resembling the well-studied rotation-activity relation. Among the slowly rotating stars, we find that magnetic flux, $Φ_\textrm{B}$, is proportional to rotation period, $P$, and among the rapidly rotating stars that average surface fields do not grow significantly beyond the level set by the available kinetic energy. Furthermore, we find close relations between nonthermal coronal X-ray emission, chromospheric H$α$ and Ca H&K emission, and magnetic flux. Taken together, these relations demonstrate empirically that the rotation-activity relation can be traced back to a dependence of the magnetic dynamo on rotation. We advocate the picture that the magnetic dynamo generates magnetic flux on the stellar surface proportional to rotation rate with a saturation limit set by the available kinetic energy, and we provide relations for average field strengths and nonthermal emission that are independent of the choice of the convective turnover time. We also find that Ca H&K emission saturates at average field strengths of $\langle B \rangle \approx 800$ G while H$α$ and X-ray emission grow further with stronger fields in the more rapidly rotating stars. This is in conflict with the coronal stripping scenario predicting that in the most rapidly rotating stars coronal plasma would be cooled to chromospheric temperatures.

preprint2022arXiv

MATISSE, the VLTI mid-infrared imaging spectro-interferometer

Context:Optical interferometry is at a key development stage. ESO&#39;s VLTI has established a stable, robust infrastructure for long-baseline interferometry for general astronomical observers. The present second-generation instruments offer a wide wavelength coverage and improved performance. Their sensitivity and measurement accuracy lead to data and images of high reliability. Aims:We have developed MATISSE, the Multi AperTure mid-Infrared SpectroScopic Experiment, to access high resolution imaging in a wide spectral domain and explore topics such: stellar activity and mass loss; planet formation and evolution in the gas and dust disks around young stars; accretion processes around super massive black holes in AGN. Methods:The instrument is a spectro-interferometric imager covering three atmospheric bands (L,M,N) from 2.8 to 13.0 mu, combining four optical beams from the VLTI&#39;s telscopes. Its concept, related observing procedure, data reduction and calibration approach are the product of 30 years of instrumental research. The instrument utilizes a multi-axial beam combination that delivers spectrally dispersed fringes. The signal provides the following quantities at several spectral resolutions: photometric flux, coherent fluxes, visibilities, closure phases, wavelength differential visibilities and phases, and aperture-synthesis imaging. Results:We provide an overview of the physical principle of the instrument and its functionalities, the characteristics of the delivered signal, a description of the observing modes and of their performance limits. An ensemble of data and reconstructed images are illustrating the first acquired key observations. Conclusion:The instrument has been in operation at Cerro Paranal, ESO, Chile since 2018, and has been open for science use by the international community since April 2019. The first scientific results are being published now.

preprint2022arXiv

New binaries from the SHINE survey

We present the multiple stellar systems observed within the SpHere INfrared survey for Exoplanet (SHINE). SHINE searched for substellar companions to young stars using high contrast imaging. Although stars with known stellar companions within SPHERE field of view (<5.5 arcsec) were removed from the original target list, we detected additional stellar companions to 78 of the 463 SHINE targets observed so far. 27% of the systems have three or more components. Given the heterogeneity of the sample in terms of observing conditions and strategy, tailored routines were used for data reduction and analysis, some of which were specifically designed for these data sets. We then combined SPHERE data with literature and archival ones, TESS light curves and Gaia parallaxes and proper motions, to characterise these systems as completely as possible. Combining all data, we were able to constrain the orbits of 25 systems. We carefully assessed the completeness of our sample for the separation range 50-500 mas (period range a few years - a few tens of years), taking into account the initial selection biases and recovering part of the systems excluded from the original list due to their multiplicity. This allowed us to compare the binary frequency for our sample with previous studies and highlight some interesting trends in the mass ratio and period distribution. We also found that, for the few objects for which such estimate was possible, the values of the masses derived from dynamical arguments were in good agreement with the model predictions. Stellar and orbital spins appear fairly well aligned for the 12 stars having enough data, which favour a disk fragmentation origin. Our results highlight the importance of combining different techniques when tackling complex problems such as the formation of binaries and show how large samples can be useful for more than one purpose.

preprint2022arXiv

The disk of FU Orionis viewed with MATISSE/VLTI: first interferometric observations in $L$ and $M$ bands

The disk of FU Orionis is marginally resolved with MATISSE, suggesting that the region emitting in the thermal infrared is rather compact. An upper limit of $\sim1.3\pm0.1$ mas (in $L$) can be given for the diameter of the disk region probed in the $L$ band, corresponding to 0.5 au at the adopted Gaia EDR3 distance. This represents the hot, gaseous region of the accretion disk. The $N$-band data indicate that the dusty passive disk is silicate-rich. Only the innermost region of said dusty disk is found to emit strongly in the $N$ band, and it is resolved at an angular size of $\sim5$ mas, which translates to a diameter of about 2 au. The observations therefore place stringent constraints for the outer radius of the inner accretion disk. Dust radiative transfer simulations with RADMC-3D provide adequate fits to the spectral energy distribution from the optical to the submillimeter and to the interferometric observables when opting for an accretion rate $\dot{M}\sim 2\times 10^{-5}\, M_\odot$ yr$^{-1}$ and assuming $M_*=0.6\, M_\odot$. Most importantly, the hot inner accretion disk&#39;s outer radius can be fixed at 0.3 au. The outer radius of the dusty disk is placed at 100 au, based on constraints from scattered-light images in the literature. The dust mass contained in the disk is $2.4\times10^{-4}\, M_\odot$, and for a typical gas-to-dust ratio of 100, the total mass in the disk is approximately 0.02 $M_\odot$. We did not find any evidence for a nearby companion in the current interferometric data, and we tentatively explored the case of disk misalignment. For the latter, our modeling results suggest that the disk orientation is similar to that found in previous imaging studies by ALMA. Should there be an asymmetry in the very compact, inner accretion disk, this might be resolved at even smaller spatial scales ($\leq1$ mas).

preprint2022arXiv

The GAPS Programme at TNG: XXXVI. Measurement of the Rossiter-McLaughlin effect and revising the physical and orbital parameters of the HAT-P-15, HAT-P-17, HAT-P-21, HAT-P-26, HAT-P-29 eccentric planetary systems

Aim: We aim to refine the orbital and physical parameters and determine the sky-projected planet orbital obliquity of five eccentric transiting planetary systems: HAT-P-15, HAT-P-17, HAT-P-21, HAT-P-26, and HAT-P-29. Each of the systems hosts a hot Jupiter, except for HAT-P-26 which hosts a Neptune-mass planet. Methods: We observed transit events of these planets with the HARPS-N spectrograph, obtaining high-precision radial velocity measurements that allow us to measure the Rossiter-McLaughlin effect for each of the target systems. We used these new HARPS-N spectra and archival data, including those from Gaia, to better characterise the stellar atmospheric parameters. The photometric parameters for four of the hot Jupiters were recalculated using 17 new transit light curves, obtained with an array of medium-class telescopes, and data from the TESS space telescope. HATNet time-series photometric data were checked for the signatures of rotation periods of the target stars and their spin axis inclination. Results: From the analysis of the Rossiter-McLaughlin effect, we derived a sky-projected obliquity of 13, -26.3, -0.7, -26 degree for HAT-P-15b, HAT-P-17b, HAT-P-21b and HAT-P-29b, respectively. Due to the quality of the data, we were not able to well constrain the sky-projected obliquity for HAT-P-26b, although a prograde orbit is favoured. The stellar activity of HAT-P-21 indicates a rotation period of 15.88 days, which allowed us to determine its true misalignment angle (25 degree). Our new analysis of the physical parameters of the five exoplanetary systems returned values compatible with those existing in the literature. Using TESS and the available transit light curves, we reviewed the orbital ephemeris for the five systems and confirmed that the HAT-P-26 system shows transit timing variations, which may tentatively be attributed to the presence of a third body.

preprint2022arXiv

The HD 260655 system: Two rocky worlds transiting a bright M dwarf at 10 pc

We report the discovery of a multi-planetary system transiting the M0 V dwarf HD 260655 (GJ 239, TOI-4599). The system consists of at least two transiting planets, namely HD 260655 b, with a period of 2.77 d, a radius of R$_b$ = 1.240$\pm$0.023 R$_\oplus$, a mass of M$_b$ = 2.14$\pm$0.34 M$_\oplus$, and a bulk density of $ρ_b$ = 6.2$\pm$1.0 g cm$^{-3}$, and HD 260655 c, with a period of 5.71 d, a radius of R$_c$ = 1.533$^{+0.051}_{-0.046}$ R$_\oplus$, a mass of M$_c$ = 3.09$\pm$0.48 M$_\oplus$, and a bulk density of $ρ_c$ = 4.7$^{+0.9}_{-0.8}$ g cm$^{-3}$. The planets were detected in transit by the TESS mission and confirmed independently with archival and new precise radial velocities obtained with the HIRES and CARMENES instruments since 1998 and 2016, respectively. At a distance of 10 pc, HD 260655 becomes the fourth closest known multi-transiting planet system after HD 219134, LTT 1445 A, and AU Mic. Due to the apparent brightness of the host star (J = 6.7 mag), both planets are among the most suitable rocky worlds known today for atmospheric studies with the JWST, both in transmission and emission.

preprint2021arXiv

A nearby transiting rocky exoplanet that is suitable for atmospheric investigation

Spectroscopy of transiting exoplanets can be used to investigate their atmospheric properties and habitability. Combining radial velocity (RV) and transit data provides additional information on exoplanet physical properties. We detect a transiting rocky planet with an orbital period of 1.467 days around the nearby red dwarf star Gliese 486. The planet Gliese 486 b is 2.81 Earth masses and 1.31 Earth radii, with uncertainties of 5%, as determined from RV data and photometric light curves. The host star is at a distance of ~8.1 parsecs, has a J-band magnitude of ~7.2, and is observable from both hemispheres of Earth. On the basis of these properties and the planet&#39;s short orbital period and high equilibrium temperature, we show that this terrestrial planet is suitable for emission and transit spectroscopy.

preprint2021arXiv

A relation between the radial velocity dispersion of young clusters and their age: Evidence for hardening as the formation scenario of massive close binaries

The majority of massive stars ($>8$ $\rm{M_{\odot}}$) in OB associations are found in close binary systems. Nonetheless, the formation mechanism of these close massive binaries is not understood yet. Using literature data, we measured the radial-velocity dispersion ($σ_\mathrm{RV}$) as a proxy for the close binary fraction in ten OB associations in the Galaxy and the Large Magellanic Cloud, spanning an age range from 1 to 6 Myrs. We find a positive trend of this dispersion with the cluster&#39;s age, which is consistent with binary hardening. Assuming a universal binary fraction of $f_\mathrm{bin}$ = 0.7, we converted the $σ_\mathrm{RV}$ behavior to an evolution of the minimum orbital period $P_\mathrm{cutoff}$ from $\sim$9.5 years at 1 Myr to $\sim$1.4 days for the oldest clusters in our sample at $\sim$6 Myr. Our results suggest that binaries are formed at larger separations, and they harden in around 1 to 2 Myrs to produce the period distribution observed in few million year-old OB binaries. Such an inward migration may either be driven by an interaction with a remnant accretion disk or with other young stellar objects present in the system. Our findings constitute the first empirical evidence in favor of migration as a scenario for the formation of massive close binaries.

preprint2021arXiv

A SPHERE survey of self-shadowed planet-forming disks

To date, nearly two hundred planet-forming disks have been imaged with high resolution. Our propensity to study bright and extended objects is however biasing our view of the disk demography. In this work, we contribute to alleviate this bias by analyzing fifteen disks targeted with VLT/SPHERE that look faint in scattered light. Sources were selected based on a low far-IR excess from the spectral energy distribution. The comparison with the ALMA images available for a few sources shows that the scattered light surveyed by these datasets is only detected from a small portion of the disk extent. The mild anti-correlation between the disk brightness and the near-IR excess demonstrates that these disks are self-shadowed: the inner disk rim intercepts much starlight and leaves the outer disk in penumbra. Based on the uniform distribution of the disk brightness in scattered light across all spectral types, self-shadowing would act similarly for inner rims at a different distance from the star. We discuss how the illumination pattern of the outer disk may evolve with time. Some objects in the sample are proposed to be at an intermediate stage toward bright disks from the literature with either no shadow or with sign of azimuthally confined shadows.

preprint2021arXiv

Exoplanet Imaging Data Challenge: benchmarking the various image processing methods for exoplanet detection

The Exoplanet Imaging Data Challenge is a community-wide effort meant to offer a platform for a fair and common comparison of image processing methods designed for exoplanet direct detection. For this purpose, it gathers on a dedicated repository (Zenodo), data from several high-contrast ground-based instruments worldwide in which we injected synthetic planetary signals. The data challenge is hosted on the CodaLab competition platform, where participants can upload their results. The specifications of the data challenge are published on our website. The first phase, launched on the 1st of September 2019 and closed on the 1st of October 2020, consisted in detecting point sources in two types of common data-set in the field of high-contrast imaging: data taken in pupil-tracking mode at one wavelength (subchallenge 1, also referred to as ADI) and multispectral data taken in pupil-tracking mode (subchallenge 2, also referred to as ADI mSDI). In this paper, we describe the approach, organisational lessons-learnt and current limitations of the data challenge, as well as preliminary results of the participants submissions for this first phase. In the future, we plan to provide permanent access to the standard library of data sets and metrics, in order to guide the validation and support the publications of innovative image processing algorithms dedicated to high-contrast imaging of planetary systems.

preprint2021arXiv

H$α$ and He I absorption in HAT-P-32 b observed with CARMENES -- Detection of Roche lobe overflow and mass loss

We analyze two high-resolution spectral transit time series of the hot Jupiter HAT-P-32 b obtained with the CARMENES spectrograph. Our new XMM-Newton X-ray observations of the system show that the fast-rotating F-type host star exhibits a high X-ray luminosity of 2.3e29~erg/s (5-100 A), corresponding to a flux of 6.9e4 erg/cm**2/s at the planetary orbit, which results in an energy-limited escape estimate of about 1e13 g/s for the planetary mass-loss rate. The spectral time series show significant, time-dependent absorption in the Halpha and He I triplet lines with maximum depths of about 3.3% and 5.3%. The mid-transit absorption signals in the Halpha and He I lines are consistent with results from one-dimensional hydrodynamic modeling, which also yields mass-loss rates on the order of 1e13 g/s. We observe an early ingress of a redshifted component of the transmission signal, which extends into a redshifted absorption component, persisting until about the middle of the optical transit. While a super-rotating wind can explain redshifted ingress absorption, we find that an up-orbit stream, transporting planetary mass in the direction of the star, also provides a plausible explanation for the pre-transit signal. This makes HAT-P-32 a benchmark system for exploring atmospheric dynamics via transmission spectroscopy.

preprint2021arXiv

HATS-74Ab, HATS-75b, HATS-76b and HATS-77b: Four Transiting Giant Planets around K and M Dwarfs

The relative rarity of giant planets around low mass stars compared with solar-type stars is a key prediction from core accretion planet formation theory. In this paper we report on the discovery of four gas giant planets that transit low mass late K and early M dwarfs. The planets HATS-74Ab (TOI 737b), HATS-75b (TOI 552b), HATS-76b (TOI 555b), and HATS-77b (TOI 730b), were all discovered from the HATSouth photometric survey and followed-up using TESS and other photometric facilities. We use the new ESPRESSO facility at the VLT to confirm and systems and measure their masses. We find that that planets have masses of 1.46 +- 0.14 MJ , 0.491 +- 0.039 MJ , 2.629 +- 0.089 MJ and 1.374 +0.100-0.074 MJ , respectively, and radii of 1.032 +- 0.021 RJ , 0.884 +- 0.013 RJ , 1.079 +- 0.031 RJ , and 1.165 +- 0.021 RJ, respectively. The planets all orbit close to their host stars with orbital periods ranging from 1.7319 d to 3.0876 d. With further work we aim to test core accretion theory by using these and further discoveries to quantify the occurrence rate of giant planets around low mass host stars.

preprint2021arXiv

ISPY -- NaCo Imaging Survey for Planets around Young stars. CenteR: the impact of centering and frame selection

Abridged: Direct imaging has made significant progress over the past decade leading to the detection of several giant planets. Observing strategies and data rates vary from instrument to instrument and wavelength, and can result in tens of thousands of images to be combined. We here present a new approach, tailored for VLT/NaCo observations performed with the Annular Groove Phase Mask (AGPM) coronagraph at $L&#39;$ filter. Our pipeline aims at improving the post-processing of the observations on two fronts: identifying the location of the star behind the AGPM to better align the science frames and performing frame selection. Our method relies on finding the position of the AGPM in the sky frame observations, and correlate it with the circular aperture of the coronagraphic mask. This relationship allows us to retrieve the location of the AGPM in the science frames, in turn allowing us to estimate the position of the star. In the process we also gather additional information useful for our frame selection approach. We tested our pipeline on several targets, and find that we improve the S/N of companions around $β$ Pictoris and R CrA by $24\pm3$ \% and $117\pm11$ \% respectively, compared to other state-of-the-art reductions. The astrometry of the point sources is slightly different but remains compatible within $3σ$ compared to published values. Finally, we find that even for NaCo observations with tens of thousands of frames, frame selection yields just marginal improvement for point sources but may improve the final images for objects with extended emission such as disks. We proposed a novel approach to identify the location of the star behind a coronagraph even when it cannot easily be determined by other methods, leading to better S/N for nearby point sources, and led a thorough study on the importance of frame selection, concluding that the improvements are marginal in most case.

preprint2021arXiv

Metallicities in M dwarfs: Investigating different determination techniques

Deriving metallicities for solar-like stars follows well-established methods, but for cooler stars such as M dwarfs, the determination is much more complicated due to forests of molecular lines that are present. Several methods have been developed in recent years to determine accurate stellar parameters for these cool stars ($T_{\rm eff} \lesssim$ 4000 K). However, significant differences can be found at times when comparing metallicities for the same star derived using different methods. In this work, we determine the effective temperatures, surface gravities, and metallicities of 18 well-studied M dwarfs observed with the CARMENES high-resolution spectrograph following different approaches, including synthetic spectral fitting, analysis of pseudo-equivalent widths, and machine learning. We analyzed the discrepancies in the derived stellar parameters, including metallicity, in several analysis runs. Our goal is to minimize these discrepancies and find stellar parameters that are more consistent with the literature values. We attempted to achieve this consistency by standardizing the most commonly used components, such as wavelength ranges, synthetic model spectra, continuum normalization methods, and stellar parameters. We conclude that although such modifications work quite well for hotter main-sequence stars, they do not improve the consistency in stellar parameters for M dwarfs, leading to mean deviations of around 50-200 K in temperature and 0.1-0.3 dex in metallicity. In particular, M dwarfs are much more complex and a standardization of the aforementioned components cannot be considered as a straightforward recipe for bringing consistency to the derived parameters. Further in-depth investigations of the employed methods would be necessary in order to identify and correct for the discrepancies that remain.

preprint2021arXiv

Silicon in the dayside atmospheres of two ultra-hot Jupiters

Atmospheres of highly irradiated gas giant planets host a large variety of atomic and ionic species. Here we observe the thermal emission spectra of the two ultra-hot Jupiters WASP-33b and KELT-20b/MASCARA-2b in the near-infrared wavelength range with CARMENES. Via high-resolution Doppler spectroscopy, we searched for neutral silicon (Si) in their dayside atmospheres. We detect the Si spectral signature of both planets via cross-correlation with model spectra. Detection levels of 4.8$\mathrmσ$ and 5.4$\mathrmσ$, respectively, are observed when assuming a solar atmospheric composition. This is the first detection of Si in exoplanet atmospheres. The presence of Si is an important finding due to its fundamental role in cloud formation and, hence, for the planetary energy balance. Since the spectral lines are detected in emission, our results also confirm the presence of an inverted temperature profile in the dayside atmospheres of both planets.

preprint2021arXiv

The extended atmosphere and circumstellar environment of the cool evolved star VX Sagittarii as seen by MATISSE

Context. VX Sgr is a cool, evolved, and luminous red star whose stellar parameters are difficult to determine, which affects its classification. Aims. We aim to spatially resolve the photospheric extent as well as the circumstellar environment. Methods. We used interferometric observations obtained with the MATISSE instrument in the L (3 to 4 μm), M (4.5 to 5 μm), and N (8 to 13 μm) bands. We reconstructed monochromatic images using the MIRA software. We used 3D radiation-hydrodynamics (RHD) simulations carried out with CO5BOLD and a uniform disc model to estimate the apparent diameter and interpret the stellar surface structures. Moreover, we employed the radiative transfer codes Optim3D and Radmc3D to compute the spectral energy distribution for the L, M, and N bands, respectively. Results. MATISSE observations unveil, for the first time, the morphology of VX Sgr across the L, M, and N bands. The reconstructed images show a complex morphology with brighter areas whose characteristics depend on the wavelength probed. We measured the angular diameter as a function of the wavelength and showed that the photospheric extent in the L and M bands depends on the opacity through the atmosphere. In addition to this, we also concluded that the observed photospheric inhomogeneities can be interpreted as convection-related surface structures. The comparison in the N band yielded a qualitative agreement between the N band spectrum and simple dust radiative transfer simulations. However, it is not possible to firmly conclude on the interpretation of the current data because of the difficulty in constraing the model parameters using the limited accuracy of our absolute flux calibration. Conclusions. MATISSE observations and the derived reconstructed images unveil the appearance of the stellar surface and circumstellar environment across a very large spectral domain for the first time.

preprint2021arXiv

The ultra-hot-Jupiter KELT-16 b: Dynamical Evolution and Atmospheric Properties

We present broad-band photometry of 30 planetary transits of the ultra-hot Jupiter KELT-16b, using five medium-class telescopes. The transits were monitored through standard B, V, R, I filters and four were simultaneously observed from different places, for a total of 36 new light curves. We used these new photometric data and those from the TESS space telescope to review the main physical properties of the KELT-16 planetary system. Our results agree with previous measurements but are more precise. We estimated the mid-transit times for each of these transits and combined them with others from the literature to obtain 69 epochs, with a time baseline extending over more than four years, and searched for transit time variations. We found no evidence for a period change, suggesting a lower limit for orbital decay at 8 Myr, with a lower limit on the reduced tidal quality factor of $Q^{\prime}_{\star}>(1.9 \pm 0.8) \times 10^5$ with $95\%$ confidence. We built up an observational, low-resolution transmission spectrum of the planet, finding evidence of the presence of optical absorbers, although with a low significance. Using TESS data, we reconstructed the phase curve finding that KELT-16b has a phase offset of $25.25 \pm 14.03$ $^{\circ}$E, a day- and night-side brightness temperature of $3190 \pm 61$ K and $2668 \pm 56$ K, respectively. Finally, we compared the flux ratio of the planet over its star at the TESS and Spitzer wavelengths with theoretical emission spectra, finding evidence of a temperature inversion in the planet&#39;s atmosphere, the chemical composition of which is preferably oxygen-rich rather than carbon-rich.

preprint2021arXiv

VLTI-MATISSE L- and N-band aperture-synthesis imaging of the unclassified B[e] star FS Canis Majoris

Context: FS Canis Majoris (FS CMa, HD 45677) is an unclassified B[e] star surrounded by an inclined dust disk. The evolutionary stage of FS CMa is still debated. Perpendicular to the circumstellar disk, a bipolar outflow was detected. Infrared aperture-synthesis imaging provides us with a unique opportunity to study the disk structure. Aims: Our aim is to study the intensity distribution of the disk of FS CMa in the mid-infrared L and N bands. Methods: We performed aperture-synthesis imaging of FS CMa with the MATISSE instrument (Multi AperTure mid-Infrared SpectroScopic Experiment) in the low spectral resolution mode to obtain images in the L and N bands. We computed radiative transfer models that reproduce the L- and N-band intensity distributions of the resolved disks. Results: We present L- and N-band aperture-synthesis images of FS CMa reconstructed in the wavelength bands of 3.4-3.8 and 8.6-9.0 micrometer. In the L-band image, the inner rim region of an inclined circumstellar disk and the central object can be seen with a spatial resolution of 2.7 milliarcsec (mas). An inner disk cavity with an angular diameter of 6x12mas is resolved. The L-band disk consists of a bright northwestern (NW) disk region and a much fainter southeastern (SE) region. The images suggest that we are looking at the bright inner wall of the NW disk rim, which is on the far side of the disk. In the N band, only the bright NW disk region is seen. In addition to deriving the inclination and the inner disk radius, fitting the reconstructed brightness distributions via radiative transfer modeling allows one to constrain the innermost disk structure, in particular the shape of the inner disk rim.

preprint2020arXiv

A hot terrestrial planet orbiting the bright M dwarf L 168-9 unveiled by TESS

We report the detection of a transiting super-Earth-sized planet (R=1.39+-0.09 Rearth) in a 1.4-day orbit around L 168-9 (TOI-134),a bright M1V dwarf (V=11, K=7.1) located at 25.15+-0.02 pc. The host star was observed in the first sector of the Transiting Exoplanet Survey Satellite (TESS) mission and, for confirmation and planet mass measurement, was followed up with ground-based photometry, seeing-limited and high-resolution imaging, and precise radial velocity (PRV) observations using the HARPS and PFS spectrographs. Combining the TESS data and PRV observations, we find the mass of L168-9 b to be 4.60+-0.56 Mearth, and thus the bulk density to be 1.74+0.44-0.33 times larger than that of the Earth. The orbital eccentricity is smaller than 0.21 (95% confidence). This planet is a Level One Candidate for the TESS Mission&#39;s scientific objective - to measure the masses of 50 small planets - and is one of the most observationally accessible terrestrial planets for future atmospheric characterization.

preprint2020arXiv

Dense Gas in a Giant Molecular Filament

Recent surveys of the Galactic plane in the dust continuum and CO emission lines reveal that large ($\gtrsim 50$~pc) and massive ($\gtrsim 10^5$~$M_\odot$) filaments, know as giant molecular filaments (GMFs), may be linked to galactic dynamics and trace the mid-plane of the gravitational potential in the Milky Way. We have imaged one entire GMF located at $l\sim$52--54$^\circ$ longitude, GMF54 ($\sim$68~pc long), in the empirical dense gas tracers using the HCN(1--0), HNC(1--0), HCO$^+$(1--0) lines, and their $^{13}$C isotopologue transitions, as well as the N$_2$H$^+$(1--0) line. We study the dense gas distribution, the column density probability density functions (N-PDFs) and the line ratios within the GMF. The dense gas molecular transitions follow the extended structure of the filament with area filling factors between 0.06 and 0.28 with respect to $^{13}$CO(1--0). We constructed the N-PDFs of H$_2$ for each of the dense gas tracers based on their column densities and assumed uniform abundance. The N-PDFs of the dense gas tracers appear curved in log-log representation, and the HCO$^+$ N-PDF has the largest log-normal width and flattest power-law slope index. Studying the N-PDFs for sub-regions of GMF54, we found an evolutionary trend in the N-PDFs that high-mass star forming and Photon-Dominate Regions (PDRs) have flatter power-law indices. The integrated intensity ratios of the molecular lines in GMF54 are comparable to those in nearby galaxies. In particular, the N$_2$H$^+$/$^{13}$CO ratio, which traces the dense gas fraction, has similar values in GMF54 and all nearby galaxies except ULIRGs.

preprint2020arXiv

Detection of the hydrogen Balmer lines in the ultra-hot Jupiter WASP-33b

Ultra-hot Jupiters (UHJs) are highly irradiated giant exoplanets with extremely high day-side temperatures, which lead to thermal dissociation of most of the molecular species. It is expected that the neutral hydrogen atom is one of the main species in the upper atmospheres of ultra-hot Jupiters. Neutral hydrogen has been detected in several UHJs by observing its Balmer line absorption. Here, we report four transit observations of the ultra-hot Jupiter WASP-33b, performed with the CARMENES and HARPS-North spectrographs, and the detection of the H$α$, H$β$, and H$γ$ lines in the planetary transmission spectrum. The combined H$α$ transmission spectrum of the four transits has an absorption depth of 0.99$\pm$0.05 %, which corresponds to an effective radius of 1.31$\pm$0.01 Rp . The strong H$α$ absorption indicates that the line probes the high-altitude thermosphere. We further fitted the three Balmer lines using the PAWN model, assuming that the atmosphere is hydrodynamic and in LTE. We retrieved a thermosphere temperature $12200^{+1300}_{-1000}$ K and a mass-loss rate ${\rm \dot{M}}=10^{11.8^{+0.6}_{-0.5}}$ g/s. The retrieved large mass-loss rate is compatible with the &#34;Balmer-driven&#34; atmospheric escape scenario, in which the stellar Balmer continua radiation in the near-ultraviolet is substantially absorbed by the excited hydrogen atoms in the planetary thermosphere.

preprint2020arXiv

Dynamical cloud formation traced by atomic and molecular gas

Context: Atomic and molecular cloud formation is a dynamical process. However, kinematic signatures of these processes are still observationally poorly constrained. Methods: Targeting the cloud-scale environment of the prototypical infrared dark cloud G28.3, we employ spectral line imaging observations of the two atomic lines HI and [CI] as well as molecular lines observations in 13CO in the 1--0 and 3--2 transitions. The analysis comprises investigations of the kinematic properties of the different tracers, estimates of the mass flow rates, velocity structure functions, a Histogram of Oriented Gradients (HOG) study as well as comparisons to simulations. Results: The central IRDC is embedded in a more diffuse envelope of cold neutral medium (CNM) traced by HI self-absorption (HISA) and molecular gas. The spectral line data as well as the HOG and structure function analysis indicate a possible kinematic decoupling of the HI from the other gas compounds. Spectral analysis and position-velocity diagrams reveal two velocity components that converge at the position of the IRDC. Estimated mass flow rates appear rather constant from the cloud edge toward the center. The velocity structure function analysis is consistent with gas flows being dominated by the formation of hierarchical structures. Conclusions: The observations and analysis are consistent with a picture where the IRDC G28 is formed at the center of two converging gas flows. While the approximately constant mass flow rates are consistent with a self-similar, gravitationally driven collapse of the cloud, external compression by, e.g., spiral arm shocks or supernovae explosions cannot be excluded yet. Future investigations should aim at differentiating the origin of such converging gas flows.

preprint2020arXiv

Exocomets: A spectroscopic survey

While exoplanets are now routinely detected, the detection of small bodies in extrasolar systems remains challenging. Since the discovery of sporadic events interpreted as exocomets (Falling Evaporating Bodies) around $β$ Pic in the early 80s, only $\sim$20 stars have been reported to host exocomet-like events. We aim to expand the sample of known exocomet-host stars, as well as to monitor the hot-gas environment around stars with previously known exocometary activity. We have obtained high-resolution optical spectra of a heterogeneous sample of 117 main-sequence stars in the spectral type range from B8 to G8. The data have been collected in 14 observing campaigns expanding over 2 years from both hemispheres. We have analysed the Ca ii K&H and Na i D lines in order to search for non-photospheric absorptions originated in the circumstellar environment, and for variable events that could be caused by outgassing of exocomet-like bodies. We have detected non-photospheric absorptions towards 50% of the sample, attributing a circumstellar origin to half of the detections (i.e. 26% of the sample). Hot circumstellar gas is detected in the metallic lines inspected via narrow stable absorptions, and/or variable blue-/red-shifted absorption events. Such variable events were found in 18 stars in the Ca ii and/or Na i lines; 6 of them are reported in the context of this work for the first time. In some cases the variations we report in the Ca ii K line are similar to those observed in $β$ Pic. While we do not find a significant trend with the age or location of the stars, we do find that the probability of finding CS gas in stars with larger vsin i is higher. We also find a weak trend with the presence of near-infrared excess, and with anomalous ($λ$ Boo-like) abundances, but this would require confirmation by expanding the sample.

preprint2020arXiv

Formation of complex molecules in translucent clouds: Acetaldehyde, vinyl alcohol, ketene, and ethanol via nonenergetic processing of C2H2 ice

Complex organic molecules (COMs) have been identified toward high- and low-mass protostars as well as molecular clouds, suggesting that these interstellar species originate from the early stage(s) of starformation. The reaction pathways resulting in COMs described by the formula C$_2$H$_\text{n}$O are still under debate. In this work, we investigate the laboratory possible solid-state reactions that involve simple hydrocarbons and OH-radicals along with H$_2$O ice under translucent cloud conditions (1$\leq$A$_V$$\leq$5 and \textit{n}$_\text{H}$$\sim$10$^3$ cm$^{-3}$). We focus on the interactions of C$_2$H$_2$ with H-atoms and OH-radicals, which are produced along the H$_2$O formation sequence on grain surfaces at 10 K. Ultra-high vacuum (UHV) experiments were performed to study the surface chemistry observed during C$_2$H$_2$ + O$_2$ + H codeposition, where O$_2$ was used for the in-situ generation of OH-radicals. Reflection absorption infrared spectroscopy (RAIRS) was applied to in situ monitor the initial and newly formed species. After that, a temperature-programmed desorption experiment combined with a Quadrupole mass spectrometer (TPD-QMS) was used as a complementary analytical tool. The investigated 10 K surface chemistry of C$_2$H$_2$ with H-atoms and OH-radicals not only results in semi and fully saturated hydrocarbons, such as ethylene (C$_2$H$_4$) and ethane (C$_2$H$_6$), but it also leads to the formation of COMs, such as vinyl alcohol, acetaldehyde, ketene, ethanol, and possibly acetic acid. It is concluded that OH-radical addition reactions to C$_2$H$_2$, acting as a molecular backbone, followed by isomerization (i.e., keto-enol tautomerization) via an intermolecular pathway and successive hydrogenation provides a so far experimentally unreported solid-state route for the formation of these species without the need of energetic input.

preprint2020arXiv

HATS-47b, HATS-48Ab, HATS-49b and HATS-72b: Four Warm Giant Planets Transiting K Dwarfs

We report the discovery of four transiting giant planets around K dwarfs. The planets HATS-47b, HATS-48Ab, HATS-49b, and HATS-72b have masses of $0.369_{-0.021}^{+0.031}$ $M_{J}$, $0.243_{-0.030}^{+0.022}$ $M_{J}$, $0.353_{-0.027}^{+0.038}$ $M_{J}$ and $0.1254\pm0.0039$ $M_{J}$, respectively, and radii of $1.117\pm0.014$ $R_{J}$, $0.800\pm0.015$ $R_{J}$, $0.765\pm0.013$ $R_{J}$, and $0.7224\pm0.0032$ $R_{J}$, respectively. The planets orbit close to their host stars with orbital periods of $3.9228$ d, $3.1317$ d, $4.1480$ d and $7.3279$ d, respectively. The hosts are main sequence K dwarfs with masses of $0.674_{-0.012}^{+0.016}$ $M_{\odot}$, $0.7279\pm0.0066$ $M_{\odot}$, $0.7133\pm0.0075$ $M_{\odot}$, and $0.7311\pm0.0028$ $M_{\odot}$ and with $V$-band magnitudes of $V = 14.829\pm0.010$, $14.35\pm0.11$, $14.998\pm0.040$ and $12.469\pm0.010$. The Super-Neptune HATS-72b (a.k.a.\ WASP-191b and TOI 294.01) was independently identified as a transiting planet candidate by the HATSouth, WASP and TESS surveys, and we present a combined analysis of all of the data gathered by each of these projects (and their follow-up programs). An exceptionally precise mass is measured for HATS-72b thanks to high-precision radial velocity (RV) measurements obtained with VLT/ESPRESSO, FEROS, HARPS and Magellan/PFS. We also incorporate TESS observations of the warm Saturn-hosting systems HATS-47 (a.k.a. TOI 1073.01), HATS-48A and HATS-49. HATS-47 was independently identified as a candidate by the TESS team, while the other two systems were not previously identified from the TESS data. The RV orbital variations are measured for these systems using Magellan/PFS. HATS-48A has a resolved $5.\!\!^{\prime\prime}4$ neighbor in Gaia~DR2, which is a common-proper-motion binary star companion to HATS-48A with a mass of $0.22$ $M_{\odot}$ and a current projected physical separation of $\sim$1,400 au.

preprint2020arXiv

HD 117214 debris disk: scattered-light images and constraints on the presence of planets

We performed observations of the Sco-Cen F star HD 117214 aiming at a search for planetary companions and the characterization of the debris disk structure. HD 117214 was observed with the SPHERE subsystems IRDIS, IFS and ZIMPOL at optical and near-IR wavelengths using angular and polarimetric differential imaging techniques. This provided the first images of scattered light from the debris disk with a spatial resolution reaching 25 mas and an inner working angle $< 0.1&#39;&#39;$. With the observations with IRDIS and IFS we derive detection limits for substellar companions. The geometrical parameters of the detected disk are constrained by fitting 3D models for the scattering of an optically thin dust disk. Investigating the possible origin of the disk gap, we introduced putative planets therein and modeled the planet-disk and planet-planet dynamical interactions. The obtained planetary architectures are compared with the detection limit curves. The debris disk has an axisymmetric ring structure with a radius of $0.42(\pm 0.01)&#39;&#39;$ or $\sim45$ au and an inclination of $71(\pm 2.5)^\circ$ and exhibits a $0.4&#39;&#39;$ ($\sim40$ au) wide inner cavity. From the polarimetric data, we derive a polarized flux contrast for the disk of $(F_{\rm pol})_{\rm disk}/F_{\rm \ast}> (3.1 \pm 1.2)\cdot 10^{-4}$ in the RI band. The fractional scattered polarized flux of the disk is eight times smaller than the fractional infrared flux excess. This ratio is similar to the one obtained for the debris disk HIP 79977 indicating that dust radiation properties are not very different between these two disks. Inside the disk cavity we achieve the high sensitivity limits on planetary companions with a mass down to $\sim 4 M_{\rm J}$ at projected radial separations between $0.2&#39;&#39;$ and $0.4&#39;&#39;$. We can exclude the stellar companions at a radial separation larger than 75 mas from the star.

preprint2020arXiv

ISPY -- NACO Imaging Survey for Planets around Young stars: Survey description and results from the first 2.5 years of observations

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

preprint2020arXiv

Multi-wavelength modelling of the circumstellar environment of the massive proto-star AFGL 2591 VLA 3

We have studied the dust density, temperature and velocity distributions of the archetypal massive young stellar object (MYSO) AFGL 2591. Given its high luminosity ($L=2 \times 10^5$ L$_\odot$) and distance ($d=3.3$ kpc), AFGL 2591 has one of the highest $\sqrt{L}/d$ ratio, giving better resolved dust emission than any other MYSO. As such, this paper provides a template on how to use resolved multi-wavelength data and radiative transfer to obtain a well-constrained 2-D axi-symmetric analytic rotating infall model. We show for the first time that the resolved dust continuum emission from Herschel 70 $μ$m observations is extended along the outflow direction, whose origin is explained in part from warm dust in the outflow cavity walls. However, the model can only explain the kinematic features from CH$_3$CN observations with unrealistically low stellar masses ($<15$ M$_\odot$), indicating that additional physical processes may be playing a role in slowing down the envelope rotation. As part of our 3-step continuum and line fitting, we have identified model parameters that can be further constrained by specific observations. High-resolution mm visibilities were fitted to obtain the disc mass (6 M$_\odot$) and radius (2200 au). A combination of SED and near-IR observations were used to estimate the luminosity and envelope mass together with the outflow cavity inclination and opening angles.

preprint2020arXiv

Ongoing flyby in the young multiple system UX Tauri

We present observations of the young multiple system UX Tauri to look for circumstellar disks and for signs of dynamical interactions. We obtained SPHERE/IRDIS deep differential polarization images in the J and H bands. We also used ALMA archival CO data. Large extended spirals are well detected in scattered light coming out of the disk of UX Tau A. The southern spiral forms a bridge between UX Tau A and C. These spirals, including the bridge connecting the two stars, all have a CO (3-2) counterpart seen by ALMA. The disk of UX Tau C is detected in scattered light. It is much smaller than the disk of UX Tau A and has a major axis along a different position angle, suggesting a misalignment. We performed PHANTOM SPH hydrodynamical models to interpret the data. The scattered light spirals, CO emission spirals and velocity patterns of the rotating disks, and the compactness of the disk of UX Tau C all point to a scenario in which UX Tau A has been perturbed very recently (about 1000 years) by the close passage of UX Tau C.

preprint2020arXiv

Population synthesis of exocometary gas around A stars

The presence of CO gas around 10-50 Myr old A stars with debris discs has sparked debate on whether the gas is primordial or secondary. Since secondary gas released from planetesimals is poor in H$_2$, it was thought that CO would quickly photodissociate never reaching the high levels observed around the majority of A stars with bright debris discs. Kral et al. 2019 showed that neutral carbon produced by CO photodissociation can effectively shield CO and potentially explain the high CO masses around 9 A stars with bright debris discs. Here we present a new model that simulates the gas viscous evolution, accounting for carbon shielding and how the gas release rate decreases with time as the planetesimal disc loses mass. We find that the present gas mass in a system is highly dependant on its evolutionary path. Since gas is lost on long timescales, it can retain a memory of the initial disc mass. Moreover, we find that gas levels can be out of equilibrium and quickly evolving from a shielded onto an unshielded state. With this model, we build the first population synthesis of gas around A stars, which we use to constrain the disc viscosity. We find a good match with a high viscosity ($α\sim0.1$), indicating that gas is lost on timescales $\sim1-10$ Myr. Moreover, our model also shows that high CO masses are not expected around FGK stars since their planetesimal discs are born with lower masses, explaining why shielded discs are only found around A stars. Finally, we hypothesise that the observed carbon cavities could be due to radiation pressure or accreting planets.

preprint2020arXiv

Retrieving scattering clouds and disequilibrium chemistry in the atmosphere of HR 8799e

Clouds are ubiquitous in exoplanet atmospheres and represent a challenge for the model interpretation of their spectra. Complex cloud models are too numerically costly for generating a large number of spectra, while more efficient models may be too strongly simplified. We aim to constrain the atmospheric properties of the directly imaged planet HR 8799e with a free retrieval approach. We use our radiative transfer code petitRADTRANS for generating spectra, which we couple to the PyMultiNest tool. We added the effect of multiple scattering which is important for treating clouds. Two cloud model parameterizations are tested: the first incorporates the mixing and settling of condensates, the second simply parameterizes the functional form of the opacity. In mock retrievals, using an inadequate cloud model may result in atmospheres that are more isothermal and less cloudy than the input. Applying our framework on observations of HR 8799e made with the GPI, SPHERE and GRAVITY, we find a cloudy atmosphere governed by disequilibrium chemistry, confirming previous analyses. We retrieve that ${\rm C/O}=0.60_{-0.08}^{+0.07}$. Other models have not yet produced a well constrained C/O value for this planet. The retrieved C/O values of both cloud models are consistent, while leading to different atmospheric structures: cloudy, or more isothermal and less cloudy. Fitting the observations with the self-consistent Exo-REM model leads to comparable results, while not constraining C/O. With data from the most sensitive instruments, retrieval analyses of directly imaged planets are possible. The inferred C/O ratio of HR 8799e is independent of the cloud model and thus appears to be a robust. This C/O is consistent with stellar, which could indicate that the HR 8799e formed outside the CO$_2$ or CO iceline. As it is the innermost planet of the system, this constraint could apply to all HR 8799 planets.

preprint2020arXiv

Searching for a cometary belt around Trappist-1 with ALMA

Low mass stars might offer today the best opportunities to detect and characterise planetary systems, especially those harbouring close-in low mass temperate planets. Among those stars, TRAPPIST-1 is exceptional since it has seven Earth-sized planets, of which three could sustain liquid water on their surfaces. Here we present new and deep ALMA observations of TRAPPIST-1 to look for an exo-Kuiper belt which can provide clues about the formation and architecture of this system. Our observations at 0.88 mm did not detect dust emission, but can place an upper limit of 23$μ$Jy if the belt is smaller than 4 au, and 0.15 mJy if resolved and 100 au in radius. These limits correspond to low dust masses of $\sim10^{-5}-10^{-2}$ $M_\oplus$, which are expected after 8 Gyr of collisional evolution unless the system was born with a $>20$ $M_\oplus$ belt of 100 km-sized planetesimals beyond 40 au or suffered a dynamical instability. This $20$ $M_\oplus$ mass upper limit is comparable to the combined mass in TRAPPIST-1 planets, thus it is possible that most of the available solid mass in this system was used to form the known planets. A similar analysis of the ALMA data on Proxima Cen leads us to conclude that a belt born with a mass $\gtrsim1$ $M_\oplus$ in 100 km-sized planetesimals could explain its putative outer belt at 30 au. We recommend that future characterisations of debris discs around low mass stars should focus on nearby and young systems if possible.

preprint2020arXiv

The asymmetric inner disk of the Herbig Ae star HD 163296 in the eyes of VLTI/MATISSE: evidence for a vortex?

Context. The inner few au region of planet-forming disks is a complex environment. High angular resolution observations have a key role in understanding the disk structure and the dynamical processes at work. Aims. In this study we aim to characterize the mid-infrared brightness distribution of the inner disk of the young intermediate-mass star HD 163296, from VLTI/MATISSE observations. Methods. We use geometric models to fit the data. Our models include a smoothed ring, a flat disk with inner cavity, and a 2D Gaussian. The models can account for disk inclination and for azimuthal asymmetries as well. We also perform numerical hydro-dynamical simulations of the inner edge of the disk. Results. Our modeling reveals a significant brightness asymmetry in the L-band disk emission. The brightness maximum of the asymmetry is located at the NW part of the disk image, nearly at the position angle of the semimajor axis. The surface brightness ratio in the azimuthal variation is $3.5 \pm 0.2$. Comparing our result on the location of the asymmetry with other interferometric measurements, we confirm that the morphology of the $r<0.3$ au disk region is time-variable. We propose that this asymmetric structure, located in or near the inner rim of the dusty disk, orbits the star. For the physical origin of the asymmetry, we tested a hypothesis where a vortex is created by Rossby wave instability, and we find that a unique large scale vortex may be compatible with our data. The half-light radius of the L-band emitting region is $0.33\pm 0.01$ au, the inclination is ${52^\circ}^{+5^\circ}_{-7^\circ}$, and the position angle is $143^\circ \pm 3^\circ$. Our models predict that a non-negligible fraction of the L-band disk emission originates inside the dust sublimation radius for $μ$m-sized grains. Refractory grains or large ($\gtrsim 10\ μ$m-sized) grains could be the origin for this emission.

preprint2020arXiv

The CARMENES search for exoplanets around M dwarfs: Convective shift and starspot constraints from chromatic radial velocities

Context. Variability caused by stellar activity represents a challenge to the discovery and characterization of terrestrial exoplanets and complicates the interpretation of atmospheric planetary signals. Aims. We aim to use a detailed modeling tool to reproduce the effect of active regions on radial velocity measurements, which aids the identification of the key parameters that have an impact on the induced variability. Methods. We analyzed the effect of stellar activity on radial velocities as a function of wavelength by simulating the impact of the properties of spots, shifts induced by convective motions, and rotation. We focused our modeling effort on the active star YZ CMi (GJ 285), which was photometrically and spectroscopically monitored with CARMENES and the Telescopi Joan Oró. Results. We demonstrate that radial velocity curves at different wavelengths yield determinations of key properties of active regions, including spot filling factor, temperature contrast, and location, thus solving the degeneracy between them. Most notably, our model is also sensitive to convective motions. Results indicate a reduced convective shift for M dwarfs when compared to solar-type stars (in agreement with theoretical extrapolations) and points to a small global convective redshift instead of blueshift. Conclusions. Using a novel approach based on simultaneous chromatic radial velocities and light curves, we can set strong constraints on stellar activity, including an elusive parameter such as the net convective motion effect.

preprint2020arXiv

The CARMENES search for exoplanets around M dwarfs. Dynamical characterization of the multiple planet system GJ 1148 and prospects of habitable exomoons around GJ 1148 b

Context. GJ 1148 is an M-dwarf star hosting a planetary system composed of two Saturn-mass planets in eccentric orbits with periods of 41.38 and 532.02 days. Aims. We reanalyze the orbital configuration and dynamics of the GJ 1148 multi-planetary system based on new precise radial velocity (RV) measurements taken with CARMENES. Methods. We combined new and archival precise Doppler measurements from CARMENES with those available from HIRES for GJ 1148 and modeled these data with a self-consistent dynamical model. We studied the orbital dynamics of the system using the secular theory and direct N-body integrations. The prospects of potentially habitable moons around GJ 1148 b were examined. Results. The refined dynamical analyses show that the GJ 1148 system is long-term stable in a large phase-space of orbital parameters with an orbital configuration suggesting apsidal alignment, but not in any particular high-order mean-motion resonant commensurability. GJ 1148 b orbits inside the optimistic habitable zone (HZ). We find only a narrow stability region around the planet where exomoons can exist. However, in this stable region exomoons exhibit quick orbital decay due to tidal interaction with the planet. Conclusions. The GJ 1148 planetary system is a very rare M-dwarf planetary system consisting of a pair of gas giants, the inner of which resides in the HZ. We conclude that habitable exomoons around GJ 1148 b are very unlikely to exist.

preprint2020arXiv

The challenge of measuring the phase function of debris disks. Application to HR\,4796

Abridged: Debris disks are valuable systems to study dust properties. Because they are optically thin at all wavelengths, we have direct access to the properties of dust grains. One very promising technique to study them is to measure their phase function. Disks that are highly inclined are promising targets as a wider range of scattering angles can be probed. The phase function is usually either inferred by comparing the observations to synthetic disk models assuming a parametrized phase function, or estimating it from the surface brightness of the disk. We argue here that the latter approach can be biased due to projection effects leading to an increase in column density along the major axis of a non flat disk. We present a novel approach to account for those column density effects. The method remains model dependent, as one still requires a disk model to estimate the density variations as a function of the scattering angle. This method allows us however to estimate the shape of the phase function without having to invoke any parametrized form. We apply our method to SPHERE/ZIMPOL observations of HR\,4796 and highlight the differences with previous measurements. Our modelling results suggest that the disk is not vertically flat at optical wavelengths. We discuss some of the caveats of the approach, mostly that our method remains blind to real local increase of the dust density, and that it cannot yet be readily applied to angular differential imaging observations. Similarly to previous studies on HR\,4796, we still cannot reconcile the full picture using a given scattering theory to explain the shape of the phase function, a long lasting problem for debris disks. Nonetheless, we argue that similar effects as the ones highlighted in this study can also bias the determination of the phase function in total intensity.

preprint2019arXiv

RefPlanets: Search for reflected light from extra-solar planets with SPHERE/ZIMPOL

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

preprint2019arXiv

The highly inflated giant planet WASP-174b

The transiting exoplanetary system WASP-174 was reported to be composed by a main-sequence F star (V=11.8 mag) and a giant planet, WASP-174b (orbital period 4.23 days). However only an upper limit was placed on the planet mass (<1.3 Mj), and a highly uncertain planetary radius (0.7-1.7 Rj) was determined. We aim to better characterise both the star and the planet and precisely measure their orbital and physical parameters. In order to constrain the mass of the planet, we obtained new measurements of the radial velocity of the star and joined them with those from the discovery paper. Photometric data from the HATSouth survey and new multi-band, high-quality (precision reached up to 0.37~mmag) photometric follow-up observations of transit events were acquired and analysed for getting accurate photometric parameters. We fit the model to all the observations, including data from the TESS space telescope, in two different modes: incorporating the stellar isochrones into the fit, and using an empirical method to get the stellar parameters. The two modes resulted to be consistent with each other to within 2 sigma. We confirm the grazing nature of the WASP-174b transits with a confidence level greater than 5 sigma, which is also corroborated by simultaneously observing the transit through four optical bands and noting how the transit depth changes due to the limb-darkening effect. We estimate that ~76% of the disk of the planet actually eclipses the parent star at mid-transit of its transit events. We find that WASP-174b is a highly-inflated hot giant planet with a mass of 0.330 Mj and a radius of 1.435 Rj, and is therefore a good target for transmission-spectroscopy observations. With a density of 0.135 g/cm^3, it is amongst the lowest-density planets ever discovered with precisely measured mass and radius.

preprint2018arXiv

Resolving faint structures in the debris disk around TWA7

Debris disks are the intrinsic by-products of the star and planet formation processes. Most likely due to instrumental limitations and their natural faintness, little is known about debris disks around low-mass stars, especially when it comes to spatially resolved observations. We present new VLT/SPHERE IRDIS Dual-Polarization Imaging (DPI) observations in which we detect the dust ring around the M2 spectral type star TWA\,7. Combined with additional Angular Differential Imaging observations we aim at a fine characterization of the debris disk and setting constraints on the presence of low-mass planets. We model the SPHERE DPI observations and constrain the location of the small dust grains, as well as the spectral energy distribution of the debris disk, using the results inferred from the observations, and perform simple N-body simulations. We find that the dust density distribution peaks at 25 au, with a very shallow outer power-law slope, and that the disk has an inclination of 13 degrees with a position angle of 90 degrees East of North. We also report low signal-to-noise detections of an outer belt at a distance of ~52 au from the star, of a spiral arm in the Southern side of the star, and of a possible dusty clump at 3.9 au. These findings seem to persist over timescales of at least a year. Using the intensity images, we do not detect any planets in the close vicinity of the star, but the sensitivity reaches Jovian planet mass upper limits. We find that the SED is best reproduced with an inner disk at 7 au and another belt at 25 au. We report the detections of several unexpected features in the disk around TWA\,7. A yet undetected 100 M$_\oplus$ planet with a semi-major axis at 20-30 au could possibly explain the outer belt as well as the spiral arm. We conclude that stellar winds are unlikely to be responsible for the spiral arm.

preprint2018arXiv

The co-existence of hot and cold gas in debris discs

Debris discs have often been described as gas-poor discs as the gas-to-dust ratio is expected to be considerably lower than in primordial,protoplanetary discs. However, recent observations have confirmed the presence of a non-negligible amount of cold gas in the circumstellar (CS) debris discs around young main-sequence stars.This cold gas has been suggested to be related to the outgassing of planetesimals and cometary-like objects. The aim of the paper is to investigate the presence of hot gas in the surroundings of stars bearing cold-gas debris discs. High-resolution optical spectra of all currently known cold-gas-bearing debris-disc systems, with the exception of $β$ Pic and Fomalhaut, have been obtained from different observatories.We have analysed the Ca II H & K and the Na I D lines searching for non-photospheric absorptions of CS origin, usually attributed to cometary-like activity. Narrow, stable Ca II and/or Na I absorption features have been detected superimposed to the photospheric lines in 10 out of the 15 observed cold-gas-bearing debris disc.Features are found at the radial velocity of the stars, or slightly blue- or red-shifted, and/or at the velocity of the local interstellar medium (ISM). Some stars also present transient variable events or absorptions extended towards red wavelengths. These are the first detections of such Ca II features in 7 out of the 15 observed stars. In some of these stars, results suggest that the stable and variable absorptions arise from relatively hot gas located in the CS close-in environment. This hot gas is detected in at least ~80%, of edge-on cold-gas-bearing debris discs, while in only ~10% of the discs seen close to face-on. We interpret this as a geometrical effect, and suggest that the non-detection of hot gas absorptions is due to the disc inclination rather than to the absence of the hot-gas component.

preprint2017arXiv

HATS-36b and 24 other transiting/eclipsing systems from the HATSouth - K2 Campaign 7 program

We report on the result of a campaign to monitor 25 HATSouth candidates using the K2 space telescope during Campaign 7 of the K2 mission. We discover HATS-36b (EPIC 215969174b), a hot Jupiter with a mass of 2.79$\pm$0.40 M$_J$ and a radius of 1.263$\pm$0.045 R$_J$ which transits a solar-type G0V star (V=14.386) in a 4.1752d period. We also refine the properties of three previously discovered HATSouth transiting planets (HATS-9b, HATS-11b, and HATS-12b) and search the K2 data for TTVs and additional transiting planets in these systems. In addition we also report on a further three systems that remain as Jupiter-radius transiting exoplanet candidates. These candidates do not have determined masses, however pass all of our other vetting observations. Finally we report on the 18 candidates which we are now able to classify as eclipsing binary or blended eclipsing binary systems based on a combination of the HATSouth data, the K2 data, and follow-up ground-based photometry and spectroscopy. These range in periods from 0.7 days to 16.7 days, and down to 1.5 mmag in eclipse depths. Our results show the power of combining ground-based imaging and spectroscopy with higher precision space-based photometry, and serve as an illustration as to what will be possible when combining ground-based observations with TESS data.