AU Mic b is a Neptune-sized planet on an 8.47-day orbit around the nearest pre-main sequence (~20 Myr) star to the Sun, the bright (
V
= 8.81) M dwarf AU Mic. The planet was preliminary detected in ...Doppler radial velocity time series and recently confirmed to be transiting with data from the TESS mission. AU Mic b is likely to be cooling and contracting and might be accompanied by a second, more massive planet, in an outer orbit. Here, we present the observations of the transit of AU Mic b using ESPRESSO on the Very Large Telescope. We obtained a high-resolution time series of spectra to measure the Rossiter-McLaughlin effect, to constrain the spin-orbit alignment of the star and planet, and to simultaneously attempt to retrieve the planet’s atmospheric transmission spectrum. These observations allowed us to study, for the first time, the early phases of the dynamical evolution of young systems. We applied different methodologies to derive the spin-orbit angle of AU Mic b, and all of them retrieve values consistent with the planet being aligned with the rotation plane of the star. We determined a conservative spin-orbit angle
λ
value of −2.96
−10.30
+10.44
degrees, indicative that the formation and migration of the planets of the AU Mic system occurred within the disc. Unfortunately, and despite the large signal-to-noise ratio of our measurements, the degree of stellar activity prevented us from detecting any features from the planetary atmosphere. In fact, our results suggest that transmission spectroscopy for recently formed planets around active young stars is going to remain very challenging, if at all possible, for the near future.
Ultra-hot Jupiters are emerging as a new class of exoplanets. Studying their chemical compositions and temperature structures will improve our understanding of their mass loss rate as well as their ...formation and evolution. We present the detection of ionized calcium in the two hottest giant exoplanets – KELT-9b and WASP-33b. By using transit datasets from CARMENES and HARPS-N observations, we achieved high-confidence-level detections of Ca II using the cross-correlation method. We further obtain the transmission spectra around the individual lines of the Ca II H&K doublet and the near-infrared triplet, and measure their line profiles. The Ca II H&K lines have an average line depth of 2.02 ± 0.17% (effective radius of 1.56 Rp) for WASP-33b and an average line depth of 0.78 ± 0.04% (effective radius of 1.47 Rp) for KELT-9b, which indicates that the absorptions are from very high upper-atmosphere layers close to the planetary Roche lobes. The observed Ca II lines are significantly deeper than the predicted values from the hydrostatic models. Such a discrepancy is probably a result of hydrodynamic outflow that transports a significant amount of Ca II into the upper atmosphere. The prominent Ca II detection with the lack of significant Ca I detection implies that calcium is mostly ionized in the upper atmospheres of the two planets.
Stellar activity can induce signals in the radial velocities of stars, complicating the detection of orbiting low-mass planets. We present a method to determine the number of planetary signals ...present in radial-velocity datasets of active stars, using only radial-velocity observations. Instead of considering separate fits with different number of planets, we use a birth-death Markov chain Monte Carlo algorithm to infer the posterior distribution for the number of planets in a single run. In a natural way, the marginal distributions for the orbital parameters of all planets are also inferred. This method is applied to HARPS data of CoRoT-7. We confidently recover the orbits of both CoRoT-7b and CoRoT-7c although the data show evidence for the presence of additional signals.
Stellar-activity features such as spots can complicate the determination of planetary parameters through spectroscopic and photometric observations. The overlap of a transiting planet and a stellar ...spot, for instance, can produce anomalies in the transit light-curves that may lead to an inaccurate estimation of the transit duration, depth, and timing. These inaccuracies can for instance affect the precise derivation of the planet radius. We present the results of a quantitative study on the effects of stellar spots on high-precision transit light-curves. We show that spot anomalies can lead to an estimate of a planet radius that is 4% smaller than the real value. Likewise, the transit duration may be estimated about 4%, longer or shorter. Depending on the size and distribution of spots, anomalies can also produce transit-timing variations (TTVs) with significant amplitudes. For instance, TTVs with signal amplitudes of 200 s can be produced when the spot is completely dark and has the size of the largest Sun spot. Our study also indicates that the smallest size of a stellar spot that still has detectable effects on a high-precision transit light-curve is around 0.03 time the stellar radius for typical Kepler Telescope precision. We also show that the strategy of including more free parameters (such as transit depth and duration) in the fitting procedure to measure the transit time of each individual transit will not produce accurate results for active stars.
Ultra-hot Jupiters (UHJs) are highly irradiated giant exoplanets with extremely high day-side temperatures, which lead to thermal dissociation of most molecular species. It is expected that the ...neutral hydrogen atom is one of the main species in the upper atmospheres of UHJs. Neutral hydrogen has been detected in several UHJs by observing their Balmer line absorption. In this work, we report four transit observations of the UHJ 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 ± 0.05%, which corresponds to an effective radius of 1.31 ± 0.01
R
p
. 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 local thermodynamic equilibrium. We retrieved a thermosphere temperature 12 200
−1000
+1300
K and a mass-loss rate
Ṁ
= 1011.8
−0.5
+0.6
g s
−1
. The retrieved high mass-loss rate is compatible with the “Balmer-driven” atmospheric escape scenario, in which the stellar Balmer continua radiation in the near-ultraviolet is substantially absorbed by excited hydrogen atoms in the planetary thermosphere.
This paper presents a fast computational technique based on the wavelet collocation method for the numerical solution of an optimal control problem governed by elliptic variational inequalities of ...obstacle type. In this problem, the solution divides the domain into contact and noncontact sets. The boundary between the contact and noncontact sets is a free boundary, which is not known a priori and the solution is not smooth on it. Accordingly, a very fine grid is needed in order to obtain a solution with a reasonable accuracy. In this paper, our aim is to propose an adaptive scheme in order to generate an appropriate and economic irregular dyadic mesh for finding the optimal control and state functions. The irregular mesh will be generated such that its density around the free boundary is higher than in other places and high-resolution computations are focused on these zones. To this aim, we use an adaptive wavelet collocation method and take advantage of the fast wavelet transform of compact-supported interpolating wavelets to develop a multi-level algorithm, which generates an adaptive computational grid. Using this adaptive grid takes less CPU time than using a full regular mesh. At each step of the algorithm, the active set method is used for solving the optimality system of the obstacle problem on the adapted mesh. Finally, the numerical examples are presented to show the validity and efficiency of the technique.
ABSTRACT Observations of the Rossiter-McLaughlin (RM) effect provide information on star-planet alignments, which can inform planetary migration and evolution theories. Here, we go beyond the ...classical RM modeling and explore the impact of a convective blueshift that varies across the stellar disk and non-Gaussian stellar photospheric profiles. We simulated an aligned hot Jupiter with a four-day orbit about a Sun-like star and injected center-to-limb velocity (and profile shape) variations based on radiative 3D magnetohydrodynamic simulations of solar surface convection. The residuals between our modeling and classical RM modeling were dependent on the intrinsic profile width and v sin i; the amplitude of the residuals increased with increasing v sin i and with decreasing intrinsic profile width. For slowly rotating stars the center-to-limb convective variation dominated the residuals (with amplitudes of 10 s of cm s−1 to ∼1 m s−1); however, for faster rotating stars the dominant residual signature was due a non-Gaussian intrinsic profile (with amplitudes from 0.5 to 9 m s−1). When the impact factor was 0, neglecting to account for the convective center-to-limb variation led to an uncertainty in the obliquity of ∼10°-20°, even though the true v sin i was known. Additionally, neglecting to properly model an asymmetric intrinsic profile had a greater impact for more rapidly rotating stars (e.g., v sin i = 6 km s−1) and caused systematic errors on the order of ∼20° in the measured obliquities. Hence, neglecting the impact of stellar surface convection may bias star-planet alignment measurements and consequently theories on planetary migration and evolution.
Aims. One of the main scientific drivers for ESPRESSO, Échelle SPectrograph, is the detection and characterization of Earth-class exoplanets. With this goal in mind, the ESPRESSO guaranteed time ...observations (GTO) Catalog identifies the best target stars for a blind search for the radial velocity (RV) signals caused by Earth-class exoplanets. Methods. Using the most complete stellar catalogs available, we screened for the most suitable G, K, and M dwarf stars for the detection of Earth-class exoplanets with ESPRESSO. For most of the stars, we then gathered high-resolution spectra from new observations or from archival data. We used these spectra to spectroscopically investigate the existence of any stellar binaries, both bound or background stars. We derived the activity level using chromospheric activity indexes using log (R′HK) $(R' _{\textrm{HK}})$ (RHK′) , as well as the projected rotational velocity v sin i. For the cases where planet companions are already known, we also looked at the possibility that additional planets may exist in the host’s habitable zone using dynamical arguments. Results. We estimated the spectroscopic contamination level, v sin i, activity, stellar parameters and chemical abundances for 249 of the most promising targets. Using these data, we selected 45 stars that match our criteria for detectability of a planet like Earth. The stars presented and discussed in this paper constitute the ESPRESSO GTO catalog for the RV blind search for Earth-class planets. They can also be used for any other work requiring a detailed spectroscopic characterization of stars in the solar neighborhood.
The Rossiter-McLaughlin (hereafter RM) effect is a key tool for measuring the projected spin-orbit angle between stellar spin axes and orbits of transiting planets. However, the measured radial ...velocity (RV) anomalies produced by this effect are not intrinsic and depend on both instrumental resolution and data reduction routines. Using inappropriate formulas to model the RM effect introduces biases, at least in the projected velocity Vsini⋆ compared to the spectroscopic value. Currently, only the iodine cell technique has been modeled, which corresponds to observations done by, e.g., the HIRES spectrograph of the Keck telescope. In this paper, we provide a simple expression of the RM effect specially designed to model observations done by the Gaussian fit of a cross-correlation function (CCF) as in the routines performed by the HARPS team. We derived a new analytical formulation of the RV anomaly associated to the iodine cell technique. For both formulas, we modeled the subplanet mean velocity vp and dispersion βp accurately taking the rotational broadening on the subplanet profile into account. We compare our formulas adapted to the CCF technique with simulated data generated with the numerical software SOAP-T and find good agreement up to Vsini⋆ ≲ 20 km s-1. In contrast, the analytical models simulating the two different observation techniques can disagree by about 10σ in Vsini⋆ for large spin-orbit misalignments. It is thus important to apply the adapted model when fitting data.