Stellar Surface Inhomogeneities as a Potential Source of the Atmospheric Signal Detected in the K2-18b Transmission Spectrum

Transmission spectroscopy of transiting exoplanets is a proven technique that can yield information on the composition and structure of a planet’s atmosphere. However, transmission spectra may be compromised by inhomogeneities in the stellar photosphere. The sub-Neptune-sized habitable zone planet K2-18b has water absorption detected in its atmosphere using data from the Hubble Space Telescope (HST). Herein, we examine whether the reported planetary atmospheric signal seen from HST transmission spectroscopy of K2-18b could instead be induced by time-varying starspots. We built a time-variable spectral model of K2-18 that is designed to match the variability amplitude seen in K2 photometric data, and we used this model to simulate 1000 HST data sets that follow the K2-18b observation strategy. More than 1% of these provide a better fit to the data than the best-fitting exoplanet atmosphere model. After resampling our simulations to generate synthetic HST observations, we find that 40% of random draws would produce an atmospheric detection at a level at least as significant as that seen in the actual HST data of K2-18b. This work illustrates that the inferred detection of an atmosphere on K2-18b may alternatively be explained by stellar spectral contamination due to the inhomogeneous photosphere of K2-18. We do not rule out a detection of water in the planet’s atmosphere, but we provide a plausible alternative that should be considered and conclude that more observations are needed to fully rule out stellar contamination.

Short-term stability of particles in the WD J0914+1914 white dwarf planetary system

ABSTRACT Nearly all known white dwarf planetary systems contain detectable rocky debris in the stellar photosphere. A glaring exception is the young and still evolving white dwarf WD J0914+1914, which instead harbours a giant planet and a disc of pure gas. The stability boundaries of this disc and the future prospects for this white dwarf to be polluted with rocks depend upon the mass and orbit of the planet, which are only weakly constrained. Here, we combine an ensemble of plausible planet orbits and masses to determine where observers should currently expect to find the outer boundary of the gas disc. We do so by performing a sweep of the entire plausible phase space with short-term numerical integrations. We also demonstrate that particle-star collisional trajectories, which would lead to the (unseen) signature of rocky metal pollution, occupy only a small fraction of the phase space, mostly limited to particle eccentricities above 0.75. Our analysis reveals that a highly inflated planet on a near-circular orbit is the type of planet which is most consistent with the current observations.

Computational tools for the spectroscopic analysis of white dwarfs

ABSTRACT The spectroscopic features of white dwarfs are formed in the thin upper layer of their stellar photosphere. These features carry information about the white dwarf’s surface temperature, surface gravity, and chemical composition (hereafter ‘labels’). Existing methods to determine these labels rely on complex ab-initio theoretical models, which are not always publicly available. Here, we present two techniques to determine atmospheric labels from white dwarf spectra: a generative fitting pipeline that interpolates theoretical spectra with artificial neural networks and a random forest regression model using parameters derived from absorption line features. We test and compare our methods using a large catalogue of white dwarfs from the Sloan Digital Sky Survey (SDSS), achieving the same accuracy and negligible bias as compared to previous studies. We package our techniques into an open-source python module ‘wdtools’ that provides a computationally inexpensive way to determine stellar labels from white dwarf spectra observed from any facility. We will actively develop and update our tool as more theoretical models become publicly available. We discuss applications of our tool in its present form to identify interesting outlier white dwarf systems including those with magnetic fields, helium-rich atmospheres, and double-degenerate binaries.

The stellar photosphere–hydrogen ionization front interaction in classical pulsators: a theoretical explanation for observed period–colour relations

ABSTRACT Period–colour and amplitude–colour (PCAC) relations can be used to probe both the hydrodynamics of outer envelope structure and evolutionary status of Cepheids and RR Lyraes. In this work, we incorporate the PCAC relations for RR Lyraes, BL Her, W Vir, and classical Cepheids in a single unifying theory that involves the interaction of the hydrogen ionization front (HIF) and stellar photosphere and the theory of stellar evolution. PC relations for RR Lyraes and classical Cepheids using the Optical Gravitational Lensing Experiment (OGLE-IV) data are found to be consistent with this theory: RR Lyraes have shallow/sloped relations at minimum/maximum light, whilst long-period (P > 10 d) Cepheids exhibit sloped/flat PC relations at minimum/maximum light. The differences in the PC relations for Cepheids and RR Lyraes can be explained based on the relative location of the HIF and stellar photosphere which changes depending on their position on the Hertzsprung–Russell diagram. We also extend our analysis of PCAC relations for type II Cepheids in the Galactic bulge, Large and Small Magellanic Clouds using OGLE-IV data. We find that BL Her stars have sloped PC relations at maximum and minimum light similar to short-period (P < 10 d) classical Cepheids. W Vir stars exhibit sloped/flat PC relation at minimum/maximum light similar to long-period classical Cepheids. We also compute state-of-the-art 1D radiation hydrodynamic models of RR Lyraes, BL Her and classical Cepheids using the radial stellar pulsation code in mesa to further test these ideas theoretically and find that the models are generally consistent with this picture. We are thus able to explain PC relations at maximum and minimum light across a broad spectrum of variable star types.

Asymmetric shocks in χ Cygni observed with linear spectropolarimetry

Aims. We derive information about the dynamics of the stellar photosphere, including pulsation, from a coherent interpretation of the linear polarisation detected in the spectral lines of the Mira star χ Cyg. Methods. From spectropolarimetric observations of χ Cyg, we performed a careful analysis of the polarisation signals observed in atomic and molecular lines, both in absorption and emission, using radiative transfer in the context of polarisation produced through two mechanisms: intrinsic polarisation and continuum depolarisation. We also explain the observed line doubling phenomenon in terms of an expanding shell in spherical geometry, which allows us to pinpoint the coordinates over the stellar disc with enhanced polarisation. Results. We find that the polarised spectrum of χ Cyg is dominated by intrinsic polarisation and has a negligible continuum depolarisation. The observed polarised signals can only be explained by assuming that this polarisation is locally enhanced by velocity fields. During the pulsation, radial velocities are not homogeneous over the disc. We map these regions of enhanced velocities. Conclusions. We set an algorithm to distinguish the origin of this polarisation in any stellar spectra of linear polarisation and to find a way to increase the signal by coherently adding many lines with an appropriated weight. Applied to the Mira star χ Cyg, we reached the unexpected result that during the pulsation, velocities are radial but not homogeneous over the disc. The reason for these local velocity enhancements are probably related to the interplay between the atmospheric pulsation dynamics and the underlying stellar convection.

Resolving the extended stellar atmospheres of asymptotic giant branch stars at (sub)millimetre wavelengths

Context. The initial conditions for mass loss during the asymptotic giant branch (AGB) phase are set in their extended atmospheres, where, among others, convection and pulsation driven shocks determine the physical conditions. Aims. High resolution observations of AGB stars at (sub)millimetre wavelengths can now directly determine the morphology, activity, density, and temperature close to the stellar photosphere. Methods. We used Atacama Large Millimeter/submillimeter Array (ALMA) high angular resolution observations to resolve the extended atmospheres of four of the nearest AGB stars: W Hya, Mira A, R Dor, and R Leo. We interpreted the observations using a parameterised atmosphere model. Results. We resolve all four AGB stars and determine the brightness temperature structure between 1 and 2 stellar radii. For W Hya and R Dor we confirm the existence of hotspots with brightness temperatures > 3000 to 10 000 K. All four stars show deviations from spherical symmetry. We find variations on a timescale of days to weeks, and for R Leo we directly measure an outward motion of the millimetre wavelength surface with a velocity of at least 10.6 ± 1.4 km s−1. For all objects but W Hya we find that the temperature-radius and size-frequency relations require the existence of a (likely inhomogeneous) layer of enhanced opacity. Conclusions. The ALMA observations provide a unique probe of the structure of the extended AGB atmosphere. We find highly variable structures of hotspots and likely convective cells. In the future, these observations can be directly compared to multi-dimensional chromosphere and atmosphere models that determine the temperature, density, velocity, and ionisation structure between the stellar photosphere and the dust formation region. However, our results show that for the best interpretation, both very accurate flux calibration and near-simultaneous observations are essential.

Light variations due to the line-driven wind instability and wind blanketing in O stars

A small fraction of the radiative flux emitted by hot stars is absorbed by their winds and redistributed towards longer wavelengths. This effect, which leads also to the heating of the stellar photosphere, is termed wind blanketing. For stars with variable winds, the effect of wind blanketing may lead to the photometric variability. We have studied the consequences of line driven wind instability and wind blanketing for the light variability of O stars. We combined the results of wind hydrodynamic simulations and of global wind models to predict the light variability of hot stars due to the wind blanketing and instability. The wind instability causes stochastic light variability with amplitude of the order of tens of millimagnitudes and a typical timescale of the order of hours for spatially coherent wind structure. The amplitude is of the order of millimagnitudes when assuming that the wind consists of large number of independent concentric cones. The variability with such amplitude is observable using present space borne photometers. We show that the simulated light curve is similar to the light curves of O stars obtained using BRITE and CoRoT satellites.

Circumstellar chemistry of Si-C bearing molecules in the C-rich AGB star IRC+10216

Silicon carbide together with amorphous carbon are the main components of dust grains in the atmospheres of C-rich AGB stars. Small gaseous Si-C bearing molecules (such as SiC, SiCSi, and SiC2) are efficiently formed close to the stellar photosphere. They likely condense onto dust seeds owing to their highly refractory nature at the lower temperatures (i.e., below about 2500 K) in the dust growth zone which extends a few stellar radii from the photosphere. Beyond this region, the abundances of Si-C bearing molecules are expected to decrease until they are eventually reformed in the outer shells of the circumstellar envelope, owing to the interaction between the gas and the interstellar UV radiation field. Our goal is to understand the time-dependent chemical evolution of Si-C bond carriers probed by molecular spectral line emission in the circumstellar envelope of IRC+10216 at millimeter wavelengths.

Molecular dust precursors in envelopes of oxygen-rich AGB stars and red supergiants

Condensation of circumstellar dust begins with formation of molecular clusters close to the stellar photosphere. These clusters are predicted to act as condensation cores at lower temperatures and allow efficient dust formation farther away from the star. Recent observations of metal oxides, such as AlO, AlOH, TiO, and TiO2, whose emission can be traced at high angular resolutions with ALMA, have allowed first observational studies of the condensation process in oxygen-rich stars. We are now in the era when depletion of gas-phase species into dust can be observed directly. I review the most recent observations that allow us to identify gas species involved in the formation of inorganic dust of AGB stars and red supergiants. I also discuss challenges we face in interpreting the observations, especially those related to non-equilibrium gas excitation and the high complexity of stellar atmospheres in the dust-formation zone.

Spectroscopic Parameters and atmosphEric ChemIstriEs of Stars (SPECIES)

Context. The detection and subsequent characterisation of exoplanets are intimately linked to the characteristics of their host star. Therefore, it is necessary to study the star in detail in order to understand the formation history and characteristics of their companion(s). Aims. Our aims are to develop a community tool that allows the automated calculation of stellar parameters for a large number of stars, using high resolution echelle spectra and minimal photometric magnitudes, and introduce the first catalogue of these measurements in this work. Methods. We measured the equivalent widths of several iron lines and used them to solve the radiative transfer equation assuming local thermodynamic equilibrium in order to obtain the atmospheric parameters (Teff, [Fe/H], logg, and ξt). We then used these values to derive the abundance of 11 chemical elements in the stellar photosphere (Na, Mg, Al, Si, Ca, Ti, Cr, Mn, Ni, Cu, and Zn). Rotation and macroturbulent velocity were obtained using temperature calibrators and synthetic line profiles to match the observed spectra of five absorption lines. Finally, by interpolating in a grid of MIST isochrones, we were able to derive the mass, radius, and age for each star using a Bayesian approach. Results. SPECIES obtains bulk parameters that are in good agreement with measured values from different existing catalogues, including when different methods are used to derive them. We find discrepancies in the chemical abundances for some elements with respect to other works, which could be produced by differences in Teff, or in the line list or the atomic line data used to derive them. We also obtained analytic relations to describe the correlations between different parameters, and we implemented new methods to better handle these correlations, which provides a better description of the uncertainties associated with the measurements.