(Abridged) The quiet M2.5 star GJ 436 hosts a warm Neptune that displays an extended atmosphere that dwarfs its own host star. Predictions of atmospheric escape in such planets state that H atoms escape from the upper atmosphere in a collisional regime and that the flow can drag heavier atoms to the upper atmosphere. It is unclear, however, what astrophysical mechanisms drive the process. Our objective is to leverage the extensive coverage of HST/COS observations of the far-ultraviolet (FUV) spectrum of GJ 436 to search for signals of metallic ions in the upper atmosphere of GJ 436 b. We analyzed flux time-series of species present in the FUV spectrum of GJ 436, as well as the Lyman-$\alpha$ line. GJ 436 displays FUV flaring events with a rate of $\sim$10 d$^{-1}$. There is evidence for a possibly long-lived active region or longitude that modulates the FUV metallic lines of the star with amplitudes up to 20%. Despite the strong geocoronal contamination in the COS spectra, we detected in-transit excess absorption signals of $\sim$50% and $\sim$30% in the blue and red wings, respectively, of the Lyman-$\alpha$ line. We rule out a wide range of excess absorption levels in the metallic lines of the star during the transit. The large atmospheric loss of GJ 436 b observed in Lyman-$\alpha$ transmission spectra is stable over the timescale of a few years, and the red wing signal supports the presence of a variable hydrogen absorption source besides the stable exosphere. The previously claimed in-transit absorption in the Si III line is likely an artifact resulting from the stellar magnetic cycle. The non-detection of metallic ions in absorption could indicate that the escape is not hydrodynamic or that the atmospheric mixing is not efficient in dragging metals high enough for sublimation to produce a detectable escape rate of ions to the exosphere.

We present HST near-ultraviolet (NUV) transits of the hot Jupiter WASP-121b, acquired as part of the PanCET program. Time series spectra during two transit events were used to measure the transmission spectra between 2280 and 3070 Angstroms at a resolution of 30,000. Using HST data from 61 STIS visits, we show that data from HST's Pointing Control System can be used to decorrelate the instrument systematic errors (Jitter Decorrelation), which we used to fit the WASP-121b light curves. The NUV spectrum show very strong absorption features, with the NUV white light curve found to be larger than the average optical and near-infrared value at 6-$\sigma$ confidence. We identify and spectrally resolve absorption from the Mg ii doublet in the planetary exosphere at a 5.9-$\sigma$ confidence level. The Mg ii doublet is observed to reach altitudes of $R_{pl}/R_{star}=0.284\pm0.037$ for the 2796 Angstrom line and $0.242\pm0.0431$ in the 2804 Angstrom line, which exceeds the Roche lobe size as viewed in transit geometry ($R_{\rm eqRL}/R_{star}$ = 0.158). We also detect and resolve strong features of the Fe ii UV1 and UV2 multiplets, and observe the lines reaching altitudes of $R_{pl}/R_{star}\approx0.3$. At these high altitudes, the atmospheric Mg ii and Fe ii gas is not gravitationally bound to the planet, and these ionized species may be hydrodynamically escaping or could be magnetically confined. Refractory Mg and Fe atoms at high altitudes also indicates that these species are not trapped into condensate clouds at depth, which places constraints on the deep interior temperature.

We present an atmospheric transmission spectrum for the ultra-hot Jupiter WASP-121b, measured using the Space Telescope Imaging Spectrograph (STIS) onboard the Hubble Space Telescope (HST). Across the 0.47-1 micron wavelength range, the data imply an atmospheric opacity comparable to - and in some spectroscopic channels exceeding - that previously measured at near-infrared wavelengths (1.15-1.65 micron). Wavelength-dependent variations in the opacity rule out a gray cloud deck at a confidence level of 3.8-sigma and may instead be explained by VO spectral bands. We find a cloud-free model assuming chemical equilibrium for a temperature of 1500K and metal enrichment of 10-30x solar matches these data well. Using a free-chemistry retrieval analysis, we estimate a VO abundance of -6.6(-0.3,+0.2) dex. We find no evidence for TiO and place a 3-sigma upper limit of -7.9 dex on its abundance, suggesting TiO may have condensed from the gas phase at the day-night limb. The opacity rises steeply at the shortest wavelengths, increasing by approximately five pressure scale heights from 0.47 to 0.3 micron in wavelength. If this feature is caused by Rayleigh scattering due to uniformly-distributed aerosols, it would imply an unphysically high temperature of 6810+/-1530K. One alternative explanation for the short-wavelength rise is absorption due to SH (mercapto radical), which has been predicted as an important product of non-equilibrium chemistry in hot Jupiter atmospheres. Irrespective of the identity of the NUV absorber, it likely captures a significant amount of incident stellar radiation at low pressures, thus playing a significant role in the overall energy budget, thermal structure, and circulation of the atmosphere.

We present a thermal emission spectrum of the bloated hot Jupiter HAT-P-32Ab from a single eclipse observation made in spatial scan mode with the Wide Field Camera 3 (WFC3) aboard the Hubble Space Telescope (HST). The spectrum covers the wavelength regime from 1.123 to 1.644 microns which is binned into 14 eclipse depths measured to an averaged precision of 104 parts-per million. The spectrum is unaffected by a dilution from the close M-dwarf companion HAT-P-32B, which was fully resolved. We complemented our spectrum with literature results and performed a comparative forward and retrieval analysis with the 1D radiative-convective ATMO model. Assuming solar abundance of the planet atmosphere, we find that the measured spectrum can best be explained by the spectrum of a blackbody isothermal atmosphere with Tp = 1995 +/- 17K, but can equally-well be described by a spectrum with modest thermal inversion. The retrieved spectrum suggests emission from VO at the WFC3 wavelengths and no evidence of the 1.4 micron water feature. The emission models with temperature profiles decreasing with height are rejected at a high confidence. An isothermal or inverted spectrum can imply a clear atmosphere with an absorber, a dusty cloud deck or a combination of both. We find that the planet can have continuum of values for the albedo and recirculation, ranging from high albedo and poor recirculation to low albedo and efficient recirculation. Optical spectroscopy of the planet's day-side or thermal emission phase curves can potentially resolve the current albedo with recirculation degeneracy.

Infrared radiation emitted from a planet contains information about the chemical composition and vertical temperature profile of its atmosphere. If upper layers are cooler than lower layers, molecular gases will produce absorption features in the planetary thermal spectrum. Conversely, if there is a stratosphere - where temperature increases with altitude - these molecular features will be observed in emission. It has been suggested that stratospheres could form in highly irradiated exoplanets, but the extent to which this occurs is unresolved both theoretically and observationally. A previous claim for the presence of a stratosphere remains open to question, owing to the challenges posed by the highly variable host star and the low spectral resolution of the measurements. Here we report a near-infrared thermal spectrum for the ultrahot gas giant WASP-121b, which has an equilibrium temperature of approximately 2,500 kelvin. Water is resolved in emission, providing a detection of an exoplanet stratosphere at 5-sigma confidence. These observations imply that a substantial fraction of incident stellar radiation is retained at high altitudes in the atmosphere, possibly by absorbing chemical species such as gaseous vanadium oxide and titanium oxide.

We present results from the first observations of the Hubble Space Telescope (HST) Panchromatic Comparative Exoplanet Treasury (PanCET) program for WASP-101b, a highly inflated hot Jupiter and one of the community targets proposed for the James Webb Space Telescope (JWST) Early Release Science (ERS) program. From a single HST Wide Field Camera 3 (WFC3) observation, we find that the near-infrared transmission spectrum of WASP-101b contains no significant H$_2$O absorption features and we rule out a clear atmosphere at 13{\sigma}. Therefore, WASP-101b is not an optimum target for a JWST ERS program aimed at observing strong molecular transmission features. We compare WASP-101b to the well studied and nearly identical hot Jupiter WASP-31b. These twin planets show similar temperature-pressure profiles and atmospheric features in the near-infrared. We suggest exoplanets in the same parameter space as WASP-101b and WASP-31b will also exhibit cloudy transmission spectral features. For future HST exoplanet studies, our analysis also suggests that a lower count limit needs to be exceeded per pixel on the detector in order to avoid unwanted instrumental systematics.

The Voyagers (V) 1 and 2 Ultraviolet Spectrometers (UVS) data harvest covers outer planets encounters, heliosphere sky-background measurements, and stellar spectrophotometry. Because their operation period overlaps with many ultraviolet missions, the V1 and V2 UVS calibration with other spectrometers are invaluable. Here we revisit the UVS calibration to assess the intriguing 243 % (V1) and 156 % (V2) sensitivity enhancements recently proposed. Using the Saturn Lyman-$\alpha$ airglow, observed in-situ by both Voyagers, and remotely by IUE, we match the Voyager values to IUE, taking into account the shape of the Saturn Lyman-$\alpha$ line observed with the Goddard High Resolution Spectrograph onboard the Hubble Space Telescope. For all known ranges of the interplanetary hydrogen density, we show that the V1 and V2 UVS sensitivities cannot be enhanced by the amounts thus far proposed. The same diagnostic holds for distinct channels covering the diffuse HeI 58.4 nm emission. Our prescription is to keep the original calibration of the Voyager UVS with a maximum uncertainty of 30 %, making both instruments some of the most stable EUV/FUV spectrographs of the history of space exploration. In that frame, we reassess the Lyman-$\alpha$ emission excess detected by Voyager UVS deep in the heliosphere, to show its consistency with the heliospheric but not the galactic origin. Our finding confirms results obtained nearly two decades ago--namely, the UVS discovery of the heliosphere distortion and the corresponding local interstellar magnetic field's obliquity ($\sim40^{\circ}$ from upwind) in the solar system neighborhood-- without requiring any revision of the Voyager UVS calibration.

By varying the profiles of stellar extreme ultraviolet (EUV) spectral energy distribution (SED), we tested the influences of stellar EUV SEDs on the physical and chemical properties of the escaping atmosphere. We apply our model to study four exoplanets, HD\,189733b, HD\,209458b, GJ \,436b, and Kepler-11b. We found that the total mass loss rates of an exoplanet, which are determined mainly by the integrated fluxes, are moderately affected by the profiles of the EUV SED, but the composition and species distributions in the atmosphere can be dramatically modified by the different profiles of the EUV SED. For exoplanets with a high hydrodynamic escape parameter ($\lambda$), the amount of atomic hydrogen produced by photoionization at different altitudes can vary by one to two orders of magnitude with the variation of stellar EUV SEDs. The effect of photoionization of H is prominent when the EUV SED is dominated by the low-energy spectral region (400-900${\AA}$), which pushes the transition of H/H$^{+}$ to low altitudes. On the contrary, the transitions of H/H$^{+}$ move to higher altitudes when most photons concentrate in the high-energy spectral region (50-400${\AA}$). For exoplanets with a low $\lambda$, the lower temperatures of the atmosphere make many chemical reactions so important that photoionization alone can no longer determine the composition of the escaping atmosphere. For HD 189733b, it is possible to explain the time variability of \lya\ between 2010 and 2011 by a change in the EUV SED of the host K star, yet invoking only thermal H\,I in the atmosphere.

The discovery of OI atoms and CII ions in the upper atmosphere of HD 209458b, made with the Hubble Space Telescope Imaging Spectrograph (STIS) using the G140L grating, showed that these heavy species fill an area comparable to the planet's Roche lobe. The derived ~10% transit absorption depths require super-thermal processes and/or supersolar abundances. From subsequent Cosmic Origins Spectrograph (COS) observations, CII absorption was reported with tentative velocity signatures, and absorption by SiIII ions was also claimed in disagreement with a negative STIS G140L detection. Here, we revisit the COS dataset showing a severe limitation in the published results from having contrasted the in-transit spectrum against a stellar spectrum averaged from separate observations, at planetary phases 0.27, 0.72, and 0.49. We find variable stellar SiIII and CII emissions that were significantly depressed not only during transit but also at phase 0.27 compared to phases 0.72 and 0.49. Their respective off-transit 7.5 and 3.1% flux variations are large compared to their reported 8.2+/-1.4% and 7.8+/-1.3% transit absorptions. Significant variations also appear in the stellar line shapes, questioning reported velocity signatures. We furthermore present archive STIS G140M transit data consistent with no SiIII absorption, with a negative result of 1.7+/-18.7 including ~15% variability. Silicon may still be present at lower ionization states, in parallel with the recent detection of extended magnesium, as MgI atoms. In this frame, the firm detection of OI and CII implying solar or supersolar abundances contradicts the recent inference of potential x20-125 subsolar metallicity for HD 209458b.

The sun and giant planets are generally thought to have the same helium abundance as that in the solar nebula from which they were formed 4.6 billion years ago. In contrast, the interstellar medium reflects current galactic conditions. The departure of current abundances from the primordial and protosolar values may help trace the processes that drive the nucleosynthesis evolution of the galaxy and planetary interior formation and evolution. The Galileo probe measured the He abundance in situ the atmosphere of Jupiter, showing that He is only slightly depleted compared to the solar value. For Saturn, contradictory estimates from past Voyager observations make its He abundance very uncertain. Here, we use He 58.4 nm dayglow measured from the outer planets by the Voyager ultraviolet spectrometers to derive the He abundance in the atmosphere of Jupiter and Saturn. We also use the solar He 58.4 nm line measured by the Solar Heliospheric Observatory to derive the He abundance inside the focusing cone. Finally, we compare He abundances derived here with primordial and protosolar values, stressing the unique opportunity offered by inner heliosphere observations and future Voyager in situ local interstellar medium measurements to derive the He abundance in the very interstellar cloud in which we reside.

In the solar system, moons largely exceed planets in number. The Kepler database has been shown to be sensitive to exomoon detection down to the mass of Mars, but the first search has been unsuccessful. Here, we use a particles-in-cell code to predict the transit of the plasma torus produced by a satellite. Despite the small size of a moon, the spatial extent of its plasma torus can be large enough to produce substantial transit absorptions. The model is used for the interpretation of Hubble Space Telescope early ingress absorptions apparently observed during WASP-12b and HD 189733b UV transits for which no consistent explanation exists. For HD 189733b an exomoon transiting $\sim 16$ $R_p$ ahead of the planet and loading $\sim 10^{29}$ C II ions/s into space is required to explain the tentative early ingress absorption observed for C II. For WASP-12b, a moon transiting $\sim 6$ $R_p$ ahead from the planet and ejecting $\sim 10^{28}$ Mg II ions per second is required to explain the NUV early ingress absorption feature. Interestingly, both HD 189733b and WASP-12b predicted satellites are outside the Hill sphere of their planets, an indication that the moons, if present, were not formed in situ but probably captured later. Finally, our simulations show a strong electromagnetic coupling between the polar regions of planets and the orbital position of the moons, an expected outcome of the unipolar induction DC circuit model. Future observations should test our predictions with a potential opportunity to unambiguously detect the first exomoon plasma torus.

We report a new diagnostic between two different states of the local interstellar medium (LISM) near our solar system using a sensitivity study constrained by several distinct and complementary observations of the LISM, solar wind, and inner heliosphere. Assuming the Interstellar Boundary Explorer (IBEX) He flow parameters for the LISM, we obtain a strength of $\sim 2.7\pm0.2$ $\mu$G and a direction pointing away from galactic coordinates $(28, 52)\pm 3^\circ$ for the interstellar magnetic field as resulting from fitting Voyager 1 & 2 in situ plasma measurements and IBEX energetic neutral atoms ribbon}. When using Ulysses parameters for the LISM He flow, we recently reported the same direction but a strength of $2.2\pm0.1$ $\mu$G. First, we notice that with Ulysses He flow, our solution is in the expected hydrogen deflection plane (HDP). In contrast, for the IBEX He flow, the solution is $\sim 20^{\circ}$ away from the corresponding HDP plane. Second, the long-term monitoring of the interplanetary H I flow speed shows a value of $\sim 26$ km/s measured at upwind from the Doppler-shift in the strong Lyman-$\alpha $ sky background emission line. All elements of diagnostics seem therefore to support Ulysses He flow parameters for the interstellar state. In that frame, we argue that reliable discrimination between superfast, subfast, or superslow states of the interstellar flow should be based on most existing in situ and remote observations used together with global modelling of the heliosphere. For commonly accepted LISM ionization rates, we show that a fast interstellar bow-shock should be standing-off upstream of the heliopause.

Detecting heavy atoms in the inflated atmospheres of giant exoplanets that orbit close to their parent stars is a key factor for understanding their bulk composition, and the processes that drive their expansion and interaction with the impinging stellar wind. Here, we use archive data obtained with the Cosmic Origins Spectrograph onboard the Hubble Space Telescope to report an absorption of ~6.4+/-1.8% by neutral oxygen during the HD 189733b transit. Scaling published HI results from a simple hydrodynamic model of HD 189733b, a vertical OI column density of ~8x10^15/cm2 produces only a 3.5% attenuation, implying that non-thermal line broadening or super-solar abundances are required. We also report evidence of short-time variability in the measured stellar flux, a variability that we analyze and compare to solar flaring activity. In that frame, we find that non-statistical uncertainties in the measured fluxes are not negligible, which calls for caution when reporting transit absorptions. Despite these uncertainties, we also show a possible detection for both a transit and early-ingress absorption in the ion CII 133.5nm lines. If confirmed, this would be the second exoplanet for which an early-ingress absorption is reported. Assuming the HD 189733b magnetosphere to be at the origin of the early absorption, we use the Parker model for the stellar wind and a particle-in-cell code for the magnetosphere to show that its orientation should be deflected ~10-30{\deg} from the planet-star line, while its nose's position should be at least ~16.7 Rp upstream of the exoplanet in order to fit the CII transit light curve. The derived stand-off distance is consistent with a surface magnetic field strength of ~5.3 Gauss for the exoplanet, and a supersonic stellar wind impinging at ~250 km/s, with a temperature of 1.2x10^5K and a density ~6.3x10^6/cm^3 at the planetary orbit, yet the fit is not unique.

Since its formation 4.6 billion years ago, our solar system has most likely crossed numerous magnetized interstellar clouds and bubbles of different sizes and contents on its path through the Milky Way. Having a reference model for how the heliosphere and interstellar winds interact is critical for understanding our current Galactic environment, and it requires untangling the roles of two major actors: the time-variable solar wind and the local interstellar magnetic field. Numerical simulations predict a distortion of the heliosphere caused by both solar wind anisotropy and interstellar magnetic field orientation. However, model comparison to deep space probes' measurements led to contradictory reports by Voyager 1 and Voyager 2 of both several crossings of the solar wind's termination shock and of the strength of the local interstellar field, with values ranging from 1.8 to 5.7 {\mu}G. Here, we show that Voyager 1 & 2 plasma, fields, and Lyman-{\alpha} sky background measurements, as well as space observations of high-energy particles of heliospheric origin, may all be explained by a rather weak interstellar field 2.2 +/- 0.1 {\mu}G pointing from Galactic coordinates (l,b) \sim (28, 52)+/- 3{\deg}. For the 2000 epoch Ulysses-based helium parameters assumed thus far, the interstellar bow shock must exist. By contrast, using the 2010 epoch IBEX-based He parameters and a stronger magnetic field leads to a plasma configuration that is not consistent with the Voyagers TS crossings. For the newly proposed interstellar He parameters, more simulations are required before one may determine whether the interstellar bow shock truly does disappear under those assumptions.

We present a follow up study of the sensitivity of the Earth's magnetosphere to solar wind activity using a particles-in-cell model [Baraka and Ben Jaffel, 2007], but here during northward IMF. The formation of the magnetospheric cavity and its elongation is obtained with the classical structure of a magnetosphere with parallel lobes. An impulsive disturbance is then applied to the system by changing the bulk velocity of the solar wind to simulate a decrease in the solar wind dynamic pressure followed by its recovery. In response to the imposed disturbance, a gap [abrupt depression] in the incoming solar wind plasma appears moving toward the Earth. The gap's size is a ~15 RE and is comparable to the sizes previously obtained for both Bz<0 and Bz =0. During the initial phase of the disturbance, the dayside magnetopause (MP) expands slower than the previous cases of IMF orientations as a result of the depression. The size of the MP expands nonlinearly due to strengthening of its outer boundary by the northward IMF. Also, during the initial 100 {\Delta}t, the MP shrank down from 13.3 RE to ~9.2 RE before it started expanding; a phenomenon that was also observed for southern IMF conditions but not during the no IMF case. As soon as they felt the solar wind depression, cusps widened at high altitude while dragged in an upright position. For the field's topology, the reconnection between magnetospheric and magnetosheath fields is clearly observed in both northward and southward cusps areas. Also, the tail region in the northward IMF condition is more confined, in contrast to the fishtail-shape obtained in the southward IMF case. An X-point is formed in the tail at ~110 RE compared to ~103 RE and ~80 RE for Bz =0 and Bz <0 respectively. Our findings are consistent with existing reports from many space observatories for which predictions are proposed to test furthermore our simulation technique.

Stellar irradiation and particles forcing strongly affect the immediate environment of extrasolar giant planets orbiting near their parent stars. Here, we use far-ultraviolet emission spectra from HD209458 in the wavelength range (1180-1710)A to bring new insight to the composition and energetic processes in play in the gas nebula around the transiting planetary companion. In that frame, we consider up-to-date atmospheric models of the giant exoplanet where we implement non-thermal line broadening to simulate the impact on the transit absorption of superthermal atoms (HI, OI, and CII) populating the upper layers of the nebula. Our sensitivity study shows that for all existing models, a significant line broadening is required for OI and probably for CII lines in order to fit the observed transit absorptions. In that frame, we show that OI and CII are preferentially heated compared to the background gas with effective temperatures as large as T_{OI}/T_B~10 for OI and T_{CII}/T_B~5 for CII. By contrast, the situation is much less clear for HI because several models could fit the Lyman-a observations including either thermal HI in an atmosphere that has a dayside vertical column [HI]~1.05x10^{21} cm^{-2}, or a less extended thermal atmosphere but with hot HI atoms populating the upper layers of the nebula. If the energetic HI atoms are either of stellar origin or populations lost from the planet and energized in the outer layers of the nebula, our finding is that most models should converge toward one hot population that has an HI vertical column in the range [HI]_{hot}(2-4)x10^{13} cm^{-2} and an effective temperature in the range T_{HI}(1-1.3)x10^6 K, but with a bulk velocity that should be rather slow.

All far ultraviolet observations of HD209458 tend to support a scenario in which the inflated hydrogen atmosphere of its planetary companion strongly absorbs the stellar \lya flux during transit. However, it was not clear how the transit absorption depends on the selected wavelength range in the stellar line profile, nor how the atomic hydrogen cloud was distributed spatially around HD209458b. Here we report a sensitivity study of observed time and spectral variations of the stellar flux. In particular, the sensitivity of the absorption depth during transit to the assumed spectral range in the stellar line profile is shown to be very weak, leading to a transit depth in the range $(8.4-8.9)%\pm 2.0%$ for all possible wavelength ranges, and thereby confirming our initially-reported absorption rate. Taking the ratio of the line profile during transit to the unperturbed line profile, we also show that the spectral signature of the absorption by the exoplanetary hydrogen nebula is symmetric and typical of a Lorentzian, optically thick medium. Our results question the adequacy of models that require a huge absorption and/or a strong asymmetry between the blue and red side of the absorption line during transit as no such features could be detected in the HST FUV absorption profile. Finally, we show that standard atmospheric models of HD209458b provide a good fit to the observed absorption profile during transit. Other hybrid models assuming a standard model with a thin layer of superthermal hydrogen on top remain possible.

Many extrasolar planets orbit closely to their parent star. Their existence raises the fundamental problem of loss and gain in their mass. For exoplanet HD209458b, reports on an unusually extended hydrogen corona and a hot layer in the lower atmosphere seem to support the scenario of atmospheric inflation by the strong stellar irradiation. However, difficulties in reconciling evaporation models with observations call for a reassessment of the problem. Here, we use HST archive data to report a new absorption rate of ~8.9% +/- 2.1% by atomic hydrogen during the HD209458b transit, and show that no sign of evaporation could be detected for the exoplanet. We also report evidence of time variability in the HD209458 Lyman-a flux, a variability that was not accounted for in previous studies, which corrupted their diagnostics. Mass loss rates thus far proposed in the literature in the range 5x(10^{10}-10^{11} g s^{-1}) must induce a spectral signature in the Lyman-a line profile of HD209458 that cannot be found in the present analysis. Either an unknown compensation effect is hiding the expected spectral feature or else the mass loss rate of neutrals from HD209458 is modest.

We report a quantitative study of the escape of Lyman-alpha photons from an inhomogeneous optically thick medium that mimics the structure of damped Lyman-alpha absorbers. Modeling the optically thick disk with 3 components (massive stars and HII regions, dust, and neutral hydrogen), we study the resulting emission line profile that may arise near the extended damped absorption profile.