Context. Atmospheric properties of exoplanets can be constrained with transit spectroscopy. At low spectral resolution, this technique is limited by the presence of clouds. The signature of atomic ...sodium (Na i), known to be present above the clouds, is a powerful probe of the upper atmosphere, where it can be best detected and characterized at high spectral resolution. Aims. Our goal is to obtain a high-resolution transit spectrum of HD 189733b in the region around the resonance doublet of Na i at 589 nm, to characterize the absorption signature that was previously detected from space at low resolution. Methods. We analyzed archival transit data of HD 189733b obtained with the HARPS spectrograph (ℛ = 115 000) at the ESO 3.6-m telescope. We performed differential spectroscopy to retrieve the transit spectrum and light curve of the planet, implementing corrections for telluric contamination and planetary orbital motion. We compared our results to synthetic transit spectra calculated from isothermal models of the planetary atmosphere. Results. We spectrally resolve the Na i D doublet and measure line contrasts of 0.64 ± 0.07% (D2) and 0.40 ± 0.07% (D1) and FWHMs of 0.52 ± 0.08 Å. This corresponds to a detection at the 10σ level of excess of absorption of 0.32 ± 0.03% in a passband of 2 × 0.75 Å centered on each line. We derive temperatures of 2600 ± 600 K and 3270 ± 330 K at altitudes of 9800 ± 2800 and 12 700 ± 2600 km in the Na i D1 and D2 line cores, respectively. We measure a temperature gradient of ~0.2 K km-1 in the region where the sodium absorption dominates the haze absorption from a comparison with theoretical models. We also detect a blueshift of 0.16 ± 0.04 Å (4σ) in the line positions. This blueshift may be the result of winds blowing at 8 ± 2 km s-1 in the upper layers of the atmosphere. Conclusions. We demonstrate the relevance of studying exoplanet atmospheres with high-resolution spectrographs mounted on 4-m-class telescopes. Our results pave the way for an in-depth characterization of physical conditions in the atmospheres of many exoplanetary systems with future spectrographs such as ESPRESSO on the VLT or HiReS and METIS on the E-ELT.
Context. Ground-based telescopes equipped with state-of-the-art spectrographs are able to obtain high-resolution transmission and emission spectra of exoplanets that probe the structure and ...composition of their atmospheres. Various atomic and molecular species, such as Na, CO, H2O have been already detected in a number of hot Jupiters. Molecular species have been observed only in the near-infrared while atomic species have been observed in the visible. In particular, the detection and abundance determination of water vapor bring important constraints to the planet formation process. Aims. We aim to search for water vapor in the atmosphere of the exoplanet HD 189733b using a high-resolution transmission spectrum in the visible obtained with HARPS. Methods. We used the atmospheric transmission code Molecfit to correct for telluric absorption features. Then we computed the high-resolution transmission spectrum of the planet using three transit datasets. We finally searched for water vapor absorption in the water band around 6500 Å using a cross-correlation technique that combines the signal of 600–900 individual lines. Results. Telluric features are corrected to the noise level. We place a 5-σ upper limit of 100 ppm on the strength of the 6500 Å water vapor band. The 1-σ precision of 20 ppm on the transmission spectrum demonstrates that space-like sensitivity can be achieved from the ground, even for a molecule that is a strong telluric absorber. Conclusions. This approach opens new possibilites for the detection of various atomic and molecular species with future instruments such as ESPRESSO at the VLT. Extrapolating from our results, we show that only one transit with ESPRESSO would be sufficient to detect water vapor on HD 189733b-like hot Jupiter with a cloud-free atmosphere. Upcoming near-IR spectrographs will be even more efficient and sensitive to a wider range of molecular species. Moreover, the detection of the same molecular species in different bands (e.g., visible and IR) is key to constrain the structure and composition of the atmosphere, such as the presence of Rayleigh scattering or aerosols (cloud and/or hazes).
Context.
Stellar activity and instrumental signals are the main limitations to the detection of Earth-like planets using the radial-velocity (RV) technique. Recent studies show that the key to ...mitigating those perturbing effects might reside in analysing the spectra themselves, rather than the RV time series and a few activity proxies.
Aims.
The goal of this paper is to demonstrate that we can reach further improvement in RV precision by performing a principal component analysis (PCA) decomposition of the shell time series, with the shell as the projection of a spectrum onto the space-normalised flux versus flux gradient.
Methods.
By performing a PCA decomposition of shell time series, it is possible to obtain a basis of first-order spectral variations that are not related to Keplerian motion. The time coefficients associated with this basis can then be used to correct for non-Dopplerian signatures in RVs.
Results.
We applied this new method on the YARARA post-processed spectra time series of HD 10700 (
τ
Ceti) and HD 128621 (
α
Cen B). On HD 10700, we demonstrate, thanks to planetary signal injections, that this new approach can successfully disentangle real Dopplerian signals from instrumental systematics. The application of this new methodology on HD 128621 shows that the strong stellar activity signal seen at the stellar rotational period and one-year aliases becomes insignificant in a periodogram analysis. The RV root mean square on the 5-yr data is reduced from 2.44 m s
−1
down to 1.73 m s
−1
. This new approach allows us to strongly mitigate stellar activity, however, noise injections tests indicate that rather high signal-to-noise ratio (S/N > 250) is required to correct for the observed activity signal on HD 128621.
Aims.
Even the most precise radial-velocity instruments gather high-resolution spectra that present systematic errors that a data reduction pipeline cannot identify and correct for efficiently by ...simply analysing a set of calibrations and a single science frame. In this paper we aim at improving the radial-velocity precision of HARPS measurements by ‘cleaning’ individual extracted spectra using the wealth of information contained in spectral time series.
Methods.
We developed YARARA, a post-processing pipeline designed to clean high-resolution spectra of instrumental systematics and atmospheric contamination. Spectra are corrected for: tellurics, interference patterns, detector stitching, ghosts, and fibre B contaminations, as well as more advanced spectral line-by-line corrections. YARARA uses principal component analysis on spectral time series with prior information to disentangle contaminations from real Doppler shifts. We applied YARARA to three systems, HD 10700, HD 215152, and HD 10180, and compared our results to the standard HARPS data reduction software and the SERVAL post-processing pipeline.
Results.
We ran YARARA on the radial-velocity dataset of three stars intensively observed with HARPS: HD 10700, HD 215152, and HD 10180. For HD 10700, we show that YARARA enables us to obtain radial-velocity measurements that present an rms smaller than 1 m s
−1
over the 13 years of the HARPS observations, which is 20% and 10% better than the HARPS data reduction software and the SERVAL post-processing pipeline, respectively. We also injected simulated planets into the data of HD 10700 and demonstrated that YARARA does not alter pure Doppler-shifted signals. For HD 215152, we demonstrated that the 1-year signal visible in the periodogram becomes marginal after processing with YARARA and that the signals of the known planets become more significant. Finally, for HD 10180, the six known exoplanets are well recovered, although different orbital parameters and planetary masses are provided by the new reduced spectra.
Conclusions.
The post-processing correction of spectra using spectral time series allows the radial-velocity precision of HARPS data to be significantly improved and demonstrates that for the extremely quiet star HD 10700 a radial-velocity rms better than 1 m s
−1
can be reached over the 13 years of HARPS observations. Since the processing proposed in this paper does not absorb planetary signals, its application to intensively followed systems is promising and will certainly result in advances in the detections of the lightest exoplanets.
Aims.We present a new list of thorium and argon emission lines in the visible obtained by analyzing high-resolution ($R=110\,000$) spectra of a ThAr hollow cathode lamp. The aim of this new line list ...is to allow significant improvements in the quality of wavelength calibration for medium- to high-resolution astronomical spectrographs. Methods.We use a series of ThAr lamp exposures obtained with the HARPS instrument (High Accuracy Radial-velocity Planet Searcher) to detect previously unknown lines, perform a systematic search for blended lines and correct individual wavelengths by determining the systematic offset of each line relative to the average wavelength solution. Results.We give updated wavelengths for more than 8400 lines over the spectral range 3785–6915 Å. The typical internal uncertainty on the line positions is estimated to be ~10 m s-1 (3.3 parts in 108 or 0.18 mÅ), which is a factor of 2–10 better than the widely used Los Alamos Atlas of the Thorium Spectrum (Palmer & Engleman 1983). The absolute accuracy of the global wavelength scale is the same as in the Los Alamos Atlas. Using this new line list on HARPS ThAr spectra, we are able to obtain a global wavelength calibration which is precise at the 20 cm s-1 level (6.7 parts in 1010 or 0.0037 mÅ). Conclusions.Several research fields in astronomy requiring high-precision wavelength calibration in the visible (e.g. radial velocity planet searches, variability of fundamental constants) should benefit from using the new line list.
Since the discovery of the transiting super-Earth CoRoT-7b, several investigations have yielded different results for the number and masses of planets present in the system, mainly owing to the ...star's high level of activity. We re-observed CoRoT-7 in 2012 January with both HARPS and CoRoT, so that we now have the benefit of simultaneous radial-velocity and photometric data. This allows us to use the off-transit variations in the star's light curve to estimate the radial-velocity variations induced by the suppression of convective blueshift and the flux blocked by starspots. To account for activity-related effects in the radial velocities which do not have a photometric signature, we also include an additional activity term in the radial-velocity model, which we treat as a Gaussian process with the same covariance properties (and hence the same frequency structure) as the light curve. Our model was incorporated into a Monte Carlo Markov Chain in order to make a precise determination of the orbits of CoRoT-7b and CoRoT-7c. We measure the masses of planets b and c to be 4.73 ± 0.95 and 13.56 ± 1.08 M⊕, respectively. The density of CoRoT-7b is (6.61 ± 1.72)(R
p
/1.58 R⊕)−3 g cm−3, which is compatible with a rocky composition. We search for evidence of an additional planet d, identified by previous authors with a period close to 9 d. We are not able to confirm the existence of a planet with this orbital period, which is close to the second harmonic of the stellar rotation at ∼7.9 d. Using Bayesian model selection, we find that a model with two planets plus activity-induced variations is most favoured.
KELT-9 b, the hottest known exoplanet, with
T
eq
~ 4400 K, is the archetype of a new planet class known as ultra-hot Jupiters. These exoplanets are presumed to have an atmosphere dominated by neutral ...and ionized atomic species. In particular, H
α
and H
β
Balmer lines have been detected in the KELT-9 b upper atmosphere, suggesting that hydrogen is filling the planetary Roche lobe and escaping from the planet. In this work, we detected
δ
Scuti-type stellar pulsation (with a period
P
puls
= 7.54 ± 0.12 h) and studied the Rossiter-McLaughlin effect (finding a spin-orbit angle
λ
= −85.01° ± 0.23°) prior to focussing on the Balmer lines (H
α
to H
ζ
) in the optical transmission spectrum of KELT-9 b. Our HARPS-N data show significant absorption for H
α
to H
δ
. The precise line shapes of the H
α
, H
β
, and H
γ
absorptions allow us to put constraints on the thermospheric temperature. Moreover, the mass loss rate, and the excited hydrogen population of KELT-9 b are also constrained, thanks to a retrieval analysis performed with a new atmospheric model. We retrieved a thermospheric temperature of
T
= 13 200
−720
+800
K and a mass loss rate of
Ṁ
= 10
12.8±0.3
g s
−1
when the atmosphere was assumed to be in hydrodynamical expansion and in local thermodynamic equilibrium (LTE). Since the thermospheres of hot Jupiters are not expected to be in LTE, we explored atmospheric structures with non-Boltzmann equilibrium for the population of the excited hydrogen. We do not find strong statistical evidence in favor of a departure from LTE. However, our non-LTE scenario suggests that a departure from the Boltzmann equilibrium may not be sufficient to explain the retrieved low number densities of the excited hydrogen. In non-LTE, Saha equilibrium departure via photo-ionization, is also likely to be necessary to explain the data.
The combination of high-contrast imaging and high-dispersion spectroscopy, which has successfully been use to detect the atmosphere of a giant planet, is one of the most promising potential probes of ...the atmosphere of Earth-size worlds. The forthcoming generation of extremely large telescopes (ELTs) may obtain sufficient contrast with this technique to detect O2 in the atmosphere of those worlds that orbit low-mass M dwarfs. This is strong motivation to carry out a census of planets around cool stars for which habitable zones can be resolved by ELTs, i.e. for M dwarfs within ~5 parsec. Our HARPS survey has been a major contributor to that sample of nearby planets. Here we report on our radial velocity observations of Ross 128 (Proxima Virginis, GJ447, HIP 57548), an M4 dwarf just 3.4 parsec away from our Sun. This source hosts an exo-Earth with a projected mass m sini = 1.35 M⊕ and an orbital period of 9.9 days. Ross 128 b receives less than 1.5 times as much flux as Earth from the Sun and its equilibrium ranges in temperature between 269 K for an Earth-like albedo and 213 K for a Venus-like albedo. Recent studies place it close to the inner edge of the conventional habitable zone. An 80-day long light curve from K2 campaign C01 demonstrates that Ross 128 b does not transit. Together with the All Sky Automated Survey (ASAS) photometry and spectroscopic activity indices, the K2 photometry shows that Ross 128 rotates slowly and has weak magnetic activity. In a habitability context, this makes survival of its atmosphere against erosion more likely. Ross 128 b is the second closest known exo-Earth, after Proxima Centauri b (1.3 parsec), and the closest temperate planet known around a quiet star. The 15 mas planet-star angular separation at maximum elongation will be resolved by ELTs (>3λ∕D) in the optical bands of O2.
When a planet transits its host star, it blocks regions of the stellar surface from view; this causes a distortion of the spectral lines and a change in the line-of-sight (LOS) velocities, known as ...the Rossiter-McLaughlin (RM) effect. Since the LOS velocities depend, in part, on the stellar rotation, the RM waveform is sensitive to the star-planet alignment (which provides information on the system’s dynamical history). We present a new RM modelling technique that directly measures the spatially-resolved stellar spectrum behind the planet. This is done by scaling the continuum flux of the (HARPS) spectra by the transit light curve, and then subtracting the in- from the out-of-transit spectra to isolate the starlight behind the planet. This technique does not assume any shape for the intrinsic local profiles. In it, we also allow for differential stellar rotation and centre-to-limb variations in the convective blueshift. We apply this technique to HD 189733 and compare to 3D magnetohydrodynamic (MHD) simulations. We reject rigid body rotation with high confidence (>99% probability), which allows us to determine the occulted stellar latitudes and measure the stellar inclination. In turn, we determine both the sky-projected (λ ≈ −0.4 ± 0.2°) and true 3D obliquity (ψ ≈ 7+12-4°). We also find good agreement with the MHD simulations, with no significant centre-to-limb variations detectable in the local profiles. Hence, this technique provides a new powerful tool that can probe stellar photospheres, differential rotation, determine 3D obliquities, and remove sky-projection biases in planet migration theories. This technique can be implemented with existing instrumentation, but will become even more powerful with the next generation of high-precision radial velocity spectrographs.