Low-mass Exoplanets Research Articles

Unexpectedly strong effect of supergranulation on the detectability of Earth twins orbiting Sun-like stars with radial velocities

Context. Magnetic activity and surface flows at different scales pertub radial velocity measurements. This affects the detectability of low-mass exoplanets. Aims. In these flows, the effect of supergranulation is not as well characterized as the other flows, and we wish to estimate its effect on the detection of Earth-like planets in the habitable zone of Sun-like stars. Methods. We produced time series of radial velocities due to oscillations, granulation, and supergranulation, and estimated the detection limit for a G2 star and a period of 300 days. We also studied in detail the behavior of the power when the signal of a 1 MEarth planet was superposed on the signal from the stellar flows. Results. We find that the detection rate does not reach 100% except for the supergranulation level we assume, which is still optimistic, and for an excellent sampling. Conclusions. We conclude that with current knowledge, it is a very challenging task to find Earth twins around Sun-like stars with our current capabilities.

Recent advances in experimental laboratory astrophysics for stellar astrophysics applications and future data needs

AbstractAccurate atomic data for line wavelengths, energy levels, line broadening such as hyperfine structure and isotope structure, and f-values, particularly for the line rich iron group elements, are needed for stellar astrophysics applications, and examples of recent measurements are given. These atomic data are essential for determination of elemental abundances in astronomical objects. With modern facilities, telescopes and spectrographs, access to underexplored regions (IR, UV, VUV), and improved stellar atmosphere models (3D, NLTE), and extremely large datasets, astronomers are tackling problems ranging from studying Galactic chemical evolution, to low mass stars and exoplanets. Such advances require improved accuracy and completeness of the atomic database for analyses of astrophysical spectra.

The metallicity–period–mass diagram of low-mass exoplanets

The number of exoplanet detections continues to grow following the development of better instruments and missions. Key steps for the understanding of these worlds comes from their characterization and its statistical studies. We explore the metallicity-period-mass diagram for known exoplanets by using an updated version of The Stellar parameters for stars With ExoplanETs CATalog (SWEET-Cat), a unique compilation of precise stellar parameters for planet-host stars provided for the exoplanet community. Here we focus on the planets with minimum mass below 30 M$_{\oplus}$ which seems to present a possible correlation in the metallicity-period-mass diagram where the mass of the planet increases with both metallicity and period. Our analysis suggests that the general observed correlation may be not fully explained by observational biases. Additional precise data will be fundamental to confirm or deny this possible correlation.

Dusty Outflows in Planetary Atmospheres: Understanding “Super-puffs” and Transmission Spectra of Sub-Neptunes

Abstract “Super-puffs” are planets with anomalously low mean densities ( ). With a low surface gravity, the extended atmosphere is susceptible to extreme hydrodynamic mass loss (“boil-off”) on a timescale that is much shorter than the system’s age. Even more puzzling, super-puffs are estimated to have a scale height of , yet recent observations revealed completely flat transmission spectra for Kepler 51b and 51d. We investigate a new scenario that explains both observations: non-static outflowing ( ) atmospheres that carry very small dust grains (∼10 Å in size, in mass fraction) to high altitudes ( ). Dust at high altitudes inflates the observed transit radius of the planet while flattening the transmission spectra. Previous static atmospheric models struggle to achieve cloud elevation and production of photochemical haze at such high altitudes. We propose to test this scenario by extending the wavelength coverage of transmission spectra. If true, dusty atmospheric outflows may affect many young ( ), low-mass ( ) exoplanets, thereby limiting our ability to study the atmospheric composition in transmission, and inflate the observed transit radius of a planet, hence obscuring the underlying mass–radius relationship.

Predictions of the WFIRST Microlensing Survey. I. Bound Planet Detection Rates

The Wide Field InfraRed Survey Telescope (WFIRST) is the next NASA astrophysics flagship mission, to follow the James Webb Space Telescope (JWST). The WFIRST mission was chosen as the top-priority large space mission of the 2010 astronomy and astrophysics decadal survey in order to achieve three primary goals: to study dark energy via a wide-field imaging survey, to study exoplanets via a microlensing survey, and to enable a guest observer program. Here we assess the ability of the several WFIRST designs to achieve the goal of the microlensing survey to discover a large sample of cold, low-mass exoplanets with semimajor axes beyond roughly one AU, which are largely impossible to detect with any other technique. We present the results of a suite of simulations that span the full range of the proposed WFIRST architectures, from the original design envisioned by the decadal survey, to the current design, which utilizes a 2.4-m telescope donated to NASA. By studying such a broad range of architectures, we are able to determine the impact of design trades on the expected yields of detected exoplanets. In estimating the yields we take particular care to ensure that our assumed Galactic model predicts microlensing event rates that match observations, consider the impact that inaccuracies in the Galactic model might have on the yields, and ensure that numerical errors in lightcurve computations do not bias the yields for the smallest mass exoplanets. For the nominal baseline WFIRST design and a fiducial planet mass function, we predict that a total of ${\sim}1400$ bound exoplanets with mass greater than ${\sim}0.1~M_{\oplus}$ should be detected, including ${\sim}200$ with mass ${\lesssim}3~M_{\oplus}$. WFIRST should have sensitivity to planets with mass down to ${\sim}0.02~M_{\oplus}$, or roughly the mass of Ganymede.

Measuring precise radial velocities and cross-correlation function line-profile variations using a Skew Normal density

Context.Stellar activity is one of the primary limitations to the detection of low-mass exoplanets using the radial-velocity (RV) technique. Stellar activity can be probed by measuring time-dependent variations in the shape of the cross-correlation function (CCF). It is therefore critical to measure with high-precision these shape variations to decorrelate the signal of an exoplanet from spurious RV signals caused by stellar activity.Aims.We propose to estimate the variations in shape of the CCF by fitting a Skew Normal (SN) density which, unlike the commonly employed Normal density, includes a Skewness parameter to capture the asymmetry of the CCF induced by stellar activity and the convective blueshift.Methods.We compared the performances of the proposed method to the commonly employed Normal density using both simulations and real observations with different levels of activity and signal-to-noise ratios.Results.When considering real observations, the correlation between the RV and the asymmetry of the CCF and between the RV and the width of the CCF are stronger when using the parameters estimated with the SN density rather than those obtained with the commonly employed Normal density. In particular, the strongest correlations have been obtained when using the mean of the SN as an estimate for the RV. This suggests that the CCF parameters estimated using a SN density are more sensitive to stellar activity, which can be helpful when estimating stellar rotational periods and when characterizing stellar activity signals. Using the proposed SN approach, the uncertainties estimated on the RV defined as the median of the SN are on average 10% smaller than the uncertainties calculated on the mean of the Normal. The uncertainties estimated on the asymmetry parameter of the SN are on average 15% smaller than the uncertainties measured on the Bisector Inverse Slope Span (BIS SPAN), which is the commonly used parameter to evaluate the asymmetry of the CCF. We also propose a new model to account for stellar activity when fitting a planetary signal to RV data. Based on simple simulations, we were able to demonstrate that this new model improves the planetary detection limits by 12% compared to the model commonly used to account for stellar activity.Conclusions.The SN density is a better model than the Normal density for characterizing the CCF since the correlations used to probe stellar activity are stronger and the uncertainties of the RV estimate and the asymmetry of the CCF are both smaller.

Data Reduction Pipeline of the TOU Optical Very High Resolution Spectrograph and Its sub-m s−1 Performance

TOU is an extremely high resolution optical spectrograph (R=$100,000$, 380-900~nm), which is designed to detect low mass exoplanets using the radial velocity technique.We describe an IDL-based radial velocity (RV) data reduction pipeline for the TOU spectrograph and its performance with stable stars. This pipeline uses a least-squares fitting algorithm to match observed stellar spectra to a high signal-to-noise template created for each star. By carefully controlling all of the error contributions to RV measurements in both the hardware and data pipeline, we have achieved $\sim$0.9\ms long-term RV precision with one of the most RV stable stars, Tau Ceti, similar to what has been achieved with HARPS. This paper presents steps and details in our data pipeline on how to reach the \sms RV precision and also all major error sources which contribute to the final RV measurement uncertainties. The lessons learned in this pipeline development can be applied to other environmentally controlled, very high resolution optical spectrographs to improve RV precision.

Aliasing in the Radial Velocities of YZ Ceti: An Ultra-short Period for YZ Ceti c?

Abstract Mid-to-late M stars are opportunistic targets for the study of low-mass exoplanets in transit because of the high planet-to-star radius ratios of their planets. Recent studies of such stars have shown that, like their early-M counterparts, they often host multi-resonant networks of small planets. Here, I reanalyze radial velocity measurements of YZ Ceti, an active M4 dwarf for which the HARPS exoplanet survey recently discovered three exoplanets on short-period (P = 4.66, 3.06, 1.97 days) orbits. My analysis finds that the orbital periods of the inner two planets cannot be uniquely determined using the published HARPS velocities. In particular, it appears likely that the 3.06 day period of YZ Ceti c is an alias and that its true period is 0.75 days. If so, the revised minimum mass of this planet is less than 0.6 Earth masses and its geometric transit probability increases to 10%. I encourage additional observations to determine the true periods of YZ Ceti b and c and suggest a search for transits at the 0.75 day period in TESS light curves.

Trapping Low-mass Planets at the Inner Edge of the Protostellar Disk

Abstract The formation of multiple close-in low-mass exoplanets is still a mystery. The challenge is to build a system wherein the outermost planet is beyond 0.2 au from the star. Here, we investigate how the prescription for type I planet migration affects the ability to trap multiple planets in a resonant chain near the inner edge of the protostellar disk. A sharp edge modeled as a hyperbolic tangent function coupled with supersonic corrections to the classical type I migration torques results in the innermost planets being pushed inside the cavity through resonant interaction with farther planets because migration is starward at slightly supersonic eccentricities. Planets below a few Earth masses are generally trapped in a resonant chain with the outermost planet near the disk edge, but long-term stability is not guaranteed. For more massive planets the migration is so fast that the eccentricity of the innermost resonant pair is excited to highly supersonic levels due to decreased damping on the innermost planet as it is pushed inside the cavity; collisions frequently occur, and the system consists of one or two intermediate-mass planets residing closer to the star than the disk’s inner edge. We found a neat pileup of resonant planets outside the disk edge only if the corotation torque does not rapidly diminish at high eccentricity. We call for detailed studies on planet migration near the disk’s inner edge, which is still uncertain, and for an improved understanding of eccentricity damping and disk torques in the supersonic regime.

Fast Coherent Differential Imaging on Ground-based Telescopes Using the Self-coherent Camera

Direct imaging and spectral characterization of exoplanets using extreme adaptive optics (ExAO) is a key science goal of future extremely large telescopes and space observatories. However, quasi-static wavefront errors will limit the sensitivity of this endeavor. Additional limitations for ground-based telescopes arise from residual AO-corrected atmospheric wavefront errors, generating millisecond-lifetime speckles that average into a halo over a long exposure. A solution to both of these problems is to use the science camera of an ExAO system as a wavefront sensor to perform a fast measurement and correction method to minimize these aberrations as soon as they are detected. We develop the framework for one such method based on the self-coherent camera (SCC) to be applied to ground-based telescopes, called Fast Atmospheric SCC Technique (FAST). We show that with the use of a specially designed coronagraph and coherent differential imaging algorithm, recording images every few milliseconds allows for a subtraction of atmospheric and static speckles while maintaining a close to unity algorithmic exoplanet throughput. Detailed simulations reach a contrast close to the photon noise limit after 30 seconds for a 1 % bandpass in H band on both 0$^\text{th}$ and 5$^\text{th}$ magnitude stars. For the 5th magnitude case, this is about 110 times better in raw contrast than what is currently achieved from ExAO instruments if we extrapolate for an hour of observing time, illustrating that sensitivity improvement from this method could play an essential role in the future detection and characterization of lower mass exoplanets.

A gap in the planetesimal disc around HD 107146 and asymmetric warm dust emission revealed by ALMA

While detecting low mass exoplanets at tens of au is beyond current instrumentation, debris discs provide a unique opportunity to study the outer regions of planetary systems. Here we report new ALMA observations of the 80-200 Myr old Solar analogue HD 107146 that reveal the radial structure of its exo-Kuiper belt at wavelengths of 1.1 and 0.86 mm. We find that the planetesimal disc is broad, extending from 40 to 140 au, and it is characterised by a circular gap extending from 60 to 100 au in which the continuum emission drops by about 50%. We also report the non-detection of the CO J=3-2 emission line, confirming that there is not enough gas to affect the dust distribution. To date, HD 107146 is the only gas-poor system showing multiple rings in the distribution of millimeter sized particles. These rings suggest a similar distribution of the planetesimals producing small dust grains that could be explained invoking the presence of one or more perturbing planets. Because the disk appears axisymmetric, such planets should be on circular orbits. By comparing N-body simulations with the observed visibilities we find that to explain the radial extent and depth of the gap, it would be required the presence of multiple low mass planets or a single planet that migrated through the disc. Interior to HD 107146's exo-Kuiper belt we find extended emission with a peak at ~20 au and consistent with the inner warm belt that was previously predicted based on 22$\mu$m excess as in many other systems. This warm belt is the first to be imaged, although unexpectedly suggesting that it is asymmetric. This could be due to a large belt eccentricity or due to clumpy structure produced by resonant trapping with an additional inner planet.

Compositional Imprints in Density–Distance–Time: A Rocky Composition for Close-in Low-mass Exoplanets from the Location of the Valley of Evaporation

Abstract We use an end-to-end model of planet formation, thermodynamic evolution, and atmospheric escape to investigate how the statistical imprints of evaporation depend on the bulk composition of planetary cores (rocky versus icy). We find that the population-wide imprints like the location of the “evaporation valley” in the distance–radius plane and the corresponding bimodal radius distribution clearly differ depending on the bulk composition of the cores. Comparison with the observed position of the valley suggests that close-in low-mass Kepler planets have a predominantly Earth-like rocky composition. Combined with the excess of period ratios outside of MMR, this suggests that low-mass Kepler planets formed inside of the water iceline but were still undergoing orbital migration. The core radius becomes visible for planets losing all primordial H/He. For planets in this “triangle of evaporation” in the distance–radius plane, the degeneracy in composition is reduced. In the observed planetary mass–mean density diagram, we identify a trend to more volatile-rich compositions with an increasing radius (R/R ⊕ ≲ 1.6 rocky; 1.6–3.0 ices, and/or H/He; ≳3: H/He). The mass–density diagram contains important information about formation and evolution. Its characteristic broken V-shape reveals the transitions from solid planets to low-mass core-dominated planets with H/He and finally to gas-dominated giants. Evaporation causes the density and orbital distance to be anticorrelated for low-mass planets in contrast to giants, where closer-in planets are less dense, likely due to inflation. The temporal evolution of the statistical properties reported here will be of interest for the PLATO 2.0 mission, which will observe the temporal dimension.

Rotational Period Measurement of Planet Host HD 147379

The star HD 147379 was recently announced as a host of a low-mass exoplanet by Reiners et al. (2017). That announcement included a discussion of the possible detection of a periodic signal in the radial velocity observations around 21.1 days, which was interpreted as rotational variability of the star. We have examined the light curve of this star as observed by the KELT transit survey, and we find a significant periodic signal in the light curve of 22.0064 days, supporting the interpretation of Reiners et al. (2017) of the rotational variation of the star.

Self-consistent atmosphere modeling with cloud formation for low-mass stars and exoplanets

Context: Low-mass stars and extrasolar planets have ultra-cool atmospheres where a rich chemistry occurs and clouds form. The increasing amount of spectroscopic observations for extrasolar planets requires self-consistent model atmosphere simulations to consistently include the formation processes that determine cloud formation and their feedback onto the atmosphere. Aims: Complement the MARCS model atmosphere suit with simulations applicable to low-mass stars and exoplanets in preparation of E-ELT, JWST, PLATO and other upcoming facilities. Methods: The MARCS code calculates stellar atmosphere models, providing self-consistent solutions of the radiative transfer and the atmospheric structure and chemistry. We combine MARCS with DRIFT, a kinetic model that describes cloud formation in ultra-cool atmospheres (seed formation, growth/ evaporation, gravitational settling, convective mixing, element depletion). Results: We present a small grid of self-consistently calculated atmosphere models for $T_ \text{eff} = 2000 - 3000$ K with solar initial abundances and $\log(g) = 4.5$. Cloud formation in stellar and sub-stellar atmospheres appears for $T_\text{eff} < 2700$ K and has a significant effect on the structure and the spectrum of the atmosphere for $T_\text{eff} < 2400$ K. We have compared the synthetic spectra of our models with observed spectra and found that they fit the spectra of mid to late type M-dwarfs and early type L-dwarfs well. We also test DRIFT-MARCS for an example exoplanet and demonstrate that our simulations reproduce the Spitzer observations for WASP-19b rather well for $T_{\rm eff}=2600$ K, $\log(g)=3.2$ and solar abundances. Our model points at an exoplanet with a deep cloud-free atmosphere with a substantial day-night energy transport and no temperature inversion.

Understanding stellar activity-induced radial velocity jitter using simultaneous K2 photometry and HARPS RV measurements

One of the best ways to improve our understanding of the stellar activity-induced signal in radial velocity (RV) measurements is through simultaneous high-precision photometric and RV observations. This is of prime importance to mitigate the RV signal induced by stellar activity and therefore unveil the presence of low-mass exoplanets. The K2 Campaign 7 and 8 field-of-views were located in the southern hemisphere, and provided a unique opportunity to gather unprecedented simultaneous high precision photometric observation with K2 and high-precision RV measurements with the HARPS spectrograph to study the relationship between photometric variability and RV jitter. We observed nine stars with different levels of activity; from quiet to very active. We probe the presence of any meaningful relation between measured RV jitter and the simultaneous photometric variation, and also other activity indicators (e.g. BIS, FWHM, $logR'_{HK}$, and F8), by evaluating the strength and significance of the correlation between RVs and each indicator. We found that for the case of very active stars, strong and significant correlations exist between almost all the observables and measured RVs; however, for lower activity levels the correlations become random. Except for the F8 which its strong correlation with RV jitter persists over a wide range of stellar activity level, and thus our result suggests that F8 might be a powerful proxy for activity induced RV jitter. Moreover, we examine the capability of two state-of-the-art modeling techniques, namely the FF' method and SOAP2.0, in accurately predicting the RV jitter amplitude using the simultaneous photometric observation. We found that for the very active stars both techniques can reasonably well predict the amplitude of the RV jitter, however, at lower activity levels the FF' method underpredicts the RV jitter amplitude.

Water in Extrasolar Planets and Implications for Habitability

Exoplanet detection missions have found thousands of planets or planet candidates outside of the Solar System—some of which are in the habitable zone, where liquid water is possible at the surface. We give an overview of the recent progress in observations of water-rich exoplanets, detection of water in the atmosphere of gas giants and less-massive targets, and modelling of the interior and evolution of water layers in exoplanets. We summarise the possible habitability of water-rich planets and discuss the potential of future missions and telescopes towards the detection of water in the atmosphere of low-mass exoplanets or on their surface.

H2 Fluorescence in M Dwarf Systems: A Stellar Origin

Abstract Observations of molecular hydrogen (H2) fluorescence are a potentially useful tool for measuring the H2 abundance in exoplanet atmospheres. This emission was previously observed in dwarfs with planetary systems. However, low signal-to-noise prevented a conclusive determination of its origin. Possible sources include exoplanetary atmospheres, circumstellar gas disks, and the stellar surface. We use observations from the “Measurements of the Ultraviolet Spectral Characteristics of Low-mass Exoplanet Host Stars” Treasury Survey to study H2 fluorescence in dwarfs. We detect fluorescence in Hubble Space Telescope spectra of 8/9 planet-hosting and 5/6 non-planet-hosting dwarfs. The detection statistics, velocity centroids, and line widths of the emission suggest a stellar origin. We calculate H2-to-stellar-ion flux ratios to compare flux levels between stars. For stars with planets, we find an average ratio of , using the fluxes of the brightest H2 feature and two stellar C iv lines. This is compared to for stars without planets, showing that the planet-hosting dwarfs do not have significant excess H2 emission. This claim is supported by the direct FUV imaging of GJ 832, where no fluorescence is observed at the expected star–planet separation. Additionally, the 3σ upper limit of 4.9 × 10−17 erg cm−2 s−1 from these observations is two orders of magnitude below the spectroscopically observed H2 flux. We constrain the location of the fluorescing H2 using 1D radiative transfer models, and find that it could reside in starspots or a ∼2500–3000 region in the lower chromosphere. The presence of this emission could complicate efforts to quantify the atmospheric abundance of H2 in exoplanets orbiting dwarfs.

Optical and Near-infrared Spectra of σ Orionis Isolated Planetary-mass Objects

Abstract We have obtained low-resolution optical (0.7–0.98 μm) and near-infrared (1.11–1.34 μm and 0.8–2.5 μm) spectra of 12 isolated planetary-mass candidates (J = 18.2–19.9 mag) of the 3 Myr σ Orionis star cluster with the aim of determining the spectroscopic properties of very young, substellar dwarfs and assembling a complete cluster mass function. We have classified our targets by visual comparison with high- and low-gravity standards and by measuring newly defined spectroscopic indices. We derived L0–L4.5 and M9–L2.5 using high- and low-gravity standards, respectively. Our targets reveal clear signposts of youth, thus corroborating their cluster membership and planetary masses (6–13 M Jup). These observations complete the σ Orionis mass function by spectroscopically confirming the planetary-mass domain to a confidence level of ∼75%. The comparison of our spectra with BT-Settl solar metallicity model atmospheres yields a temperature scale of 2350–1800 K and a low surface gravity of log g ≈ 4.0 [cm s−2], as would be expected for young planetary-mass objects. We discuss the properties of the cluster’s least-massive population as a function of spectral type. We have also obtained the first optical spectrum of S Ori 70, a T dwarf in the direction of σ Orionis. Our data provide reference optical and near-infrared spectra of very young L dwarfs and a mass function that may be used as templates for future studies of low-mass substellar objects and exoplanets. The extrapolation of the σ Orionis mass function to the solar neighborhood may indicate that isolated planetary-mass objects with temperatures of ∼200–300 K and masses in the interval 6–13 M Jup may be as numerous as very low-mass stars.

Stacked Bayesian general Lomb-Scargle periodogram: Identifying stellar activity signals

Context. Distinguishing between a signal induced by either stellar activity or a planet is currently the main challenge in radial velocity searches for low-mass exoplanets. Even when the presence of a transiting planet and hence its period are known, stellar activity can be the main barrier to measuring the correct amplitude of the radial velocity signal. Several tools are being used to help understand which signals come from stellar activity in the data. Aims. We aim to present a new tool that can be used for the purpose of identifying periodicities caused by stellar activity, and show how it can be used to track the signal-to-noise ratio (SNR) of the detection over time. The tool is based on the principle that stellar activity signals are variable and incoherent. Methods. We calculate the Bayesian general Lomb-Scargle periodogram for subsets of data and by adding one extra data point we track what happens to the presence and significance of periodicities in the data. Publicly available datasets from HARPS and HARPS-N were used for this purpose. Additionally, we analysed a synthetic dataset that we created with SOAP2.0 to simulate pure stellar activity and a mixture of stellar activity and a planet. Results. We find that this tool can easily be used to identify unstable and incoherent signals, such as those introduced by stellar activity. The SNR of the detection grows approximately as the square root of the number of data points, in the case of a stable signal. This can then be used to make decisions on whether it is useful to keep observing a specific object. The tool is relatively fast and easy to use, and thus lends itself perfectly to a quick analysis of the data.

A Low-mass Exoplanet Candidate Detected by K2 Transiting the Praesepe M Dwarf JS 183

Abstract We report the discovery of a repeating photometric signal from a low-mass member of the Praesepe open cluster that we interpret as a Neptune-sized transiting planet. The star is JS 183 (HSHJ 163, EPIC 211916756), with T eff = 3325 ± 100 K, M * = 0.44 ± 0.04 M ⊙, R * = 0.44 ± 0.03 R ⊙, and . The planet has an orbital period of 10.134588 days and a radius of R P = 0.32 ± 0.02 R J. Since the star is faint at V = 16.5 and J = 13.3, we are unable to obtain a measured radial velocity orbit, but we can constrain the companion mass to below about 1.7 M J, and thus well below the planetary boundary. JS 183b (since designated as K2-95b) is the second transiting planet found with K2 that resides in a several-hundred-megayear open cluster; both planets orbit mid-M dwarf stars and are approximately Neptune sized. With a well-determined stellar density from the planetary transit, and with an independently known metallicity from its cluster membership, JS 183 provides a particularly valuable test of stellar models at the fully convective boundary. We find that JS 183 is the lowest-density transit host known at the fully convective boundary, and that its very low density is consistent with current models of stars just above the fully convective boundary but in tension with the models just below the fully convective boundary.