scispace - formally typeset
Search or ask a question
Journal ArticleDOI

HST hot-Jupiter transmission spectral survey: haze in the atmosphere of WASP-6b

TL;DR: Nikolov et al. as discussed by the authors reported the optical to near-infrared transmission spectroscopy of the hot-Jupiter WASP-6b, measured with the Space Telescope Imaging Spectrograph and Spitzer's InfraRed Array Camera.
Abstract: Author(s): Nikolov, N; Sing, DK; Burrows, AS; Fortney, JJ; Henry, GW; Pont, F; Ballester, GE; Aigrain, S; Wilson, PA; Huitson, CM; Gibson, NP; Desert, JM; Lecavelier des Etangs, A; Showman, AP; Vidal-Madjar, A; Wakeford, HR; Zahnle, K | Abstract: We report Hubble Space Telescope optical to near-infrared transmission spectroscopy of the hot-Jupiter WASP-6b, measured with the Space Telescope Imaging Spectrograph and Spitzer's InfraRed Array Camera. The resulting spectrum covers the range 0.29-4.5 μm. We find evidence for modest stellar activity of WASP-6 and take it into account in the transmission spectrum. The overall main characteristic of the spectrum is an increasing radius as a function of decreasing wavelength corresponding to a change of Δ (Rp/R*) = 0.0071 from 0.33 to 4.5 μm. The spectrum suggests an effective extinction cross-section with a power law of index consistent with Rayleigh scattering, with temperatures of 973 ± 144K at the planetary terminator. We compare the transmission spectrum with hot-Jupiter atmospheric models including condensate-free and aerosol-dominated models incorporating Mie theory. While none of the clear-atmosphere models is found to be in good agreement with the data, we find that the complete spectrum can be described by models that include significant opacity from aerosols including Fe-poor Mg2 SiO4, MgSiO3, KCl and Na2S dust condensates.WASP- 6b is the second planet after HD 189733b which has equilibrium temperatures near ~1200K and shows prominent atmospheric scattering in the optical.

Summary (5 min read)

1.1 The WASP-6b system

  • An analysis with new radial velocity data from Husnoo et al. (2012) brought evidence for non-significant eccentricity.
  • Typical STIS G430L and G750L spectra (blue and red continuous lines, respectively).
  • The atmosphere of WASP-6b has recently been probed in the optical regime from the ground by Jordán et al. (2013) .
  • The authors found evidence for an atmospheric haze characterized by a decreasing apparent planetary size with wavelength and no evidence for the pressure-broadened alkali features.

2.1 HST STIS spectroscopy

  • The three HST visits were scheduled such that the second and third spacecraft orbits occurred between the second and third contacts of the planetary transit in order to provide good sampling of the planetary radius while three orbits secured the stellar flux level before and after the transit.
  • The raw STIS data was reduced (bias-, dark-and flat-corrected) using the latest version of the CALSTIS 1 pipeline and the relevant up-to-date calibration frames.
  • In addition, the authors also corrected all pixels identified by CALSTIS as 'bad' with the same procedure, which together with the cosmic ray identified pixels resulted in a total of ∼11 per cent interpolated pixels.
  • The STIS spectra were then used to extract both white-light spectrophotometric time series and custom wavelength bands after integrating the appropriate flux from each bandpass.

2.3 Stellar variability monitoring

  • Stellar activity can complicate the interpretation of exoplanet transmission spectra, especially when multi-instrument multi-epoch data sets are combined (Pont et al. 2013) .
  • Each nightly observation of WASP-6 consists of 4-10 consecutive exposures that include several comparison stars in the same field of view.
  • The standard deviations of groups 2, 4 and 5 significantly exceed the measurement precision and thus indicate low-level activity due to star spots.
  • The authors take this period to be their best determination of the star's rotation period.

3 ANALYSIS

  • The authors adopt similar analysis methods for the whole HST survey that are similar to the approach detailed in Sing et al.
  • The authors fit each HST and Spitzer transit light curve in flux employing a two-component model that consists of a transit model multiplied by a systematics model.
  • The systematics models are different for the HST and Spitzer data and are detailed below.
  • The authors initially fixed the orbital period of the planet to its literature value (Gillon et al. 2009 ) before being updated to the final value reported in Section 3.2.

3.0.1 STIS white light-curve modelling

  • The three STIS white-light curves were fitted simultaneously with common inclination, semimajor axis and a planet-to-star radius contrast.
  • The final systematics model for the red data include a fourth-order polynomial of HST orbital phase, a linear time term, and ω with BIC of 77.
  • Finally, the authors find that the fitted transit parameters are very well constrained (including the transit depth) and are consistent regardless of the exact systematics model (see Fig. 5 ).
  • The final results for the uncertainties of the fitted parameters were taken from MPFIT after the authors rescaled the errors per data point based on the standard deviation of the residuals.
  • In these fits, the authors report the fitted parameters with uncertainties assuming pure photon noise and no rescaling.

3.1 System parameters and transit ephemeris

  • The authors combined the central transit times from their HST data set with the transit times reported by Dragomir et al. (2011) and Jordán et al. (2013) to derive an updated transit ephemeris of WASP-6 b.
  • In addition, the authors obtained the archival light curve with full transit coverage reported in Gillon et al. (2009) 5 and performed a fit to determine its central time.
  • The authors excluded the two Spitzer IRAC light curves as these originate from incomplete transits and are not useful for determination of accurate transit ephemeris.
  • The three measured transit radii in a bin between 5250-5400 Å taken between the two STIS G430L visits and the STIS G750L visit.
  • In addition, the authors refined the system parameters including the orbital inclination, i and normalized semimajor axis, (a/R * ).

3.2 An evaluation of the stellar activity from the HST data

  • The overlapping region between the two STIS gratings offers the opportunity to measure the photometric variability of the host star over the time interval of 40 d covered by the three HST visits, as significant stellar activity from stellar spots would change the measured transit radii.
  • 1 from Nikolov et al. 2014 ) and measured the quantity R p /R * from a fit, keeping the system parameters, ephemeris and limb darkening to their best-fitting white light values and theoretical coefficients, respectively.
  • The authors result suggests that the three HST visits occurred at similar stellar flux levels implying no significant offset between the blue and red STIS sections of the transmission spectrum.
  • Finally, the very similar radius values for the two blue visits obtained six days apart imply the lack of variability related with the planetary atmosphere itself.

3.3 Fits to transmission spectra

  • Light curves from various spectral bins were extracted from the STIS/G430L and G750L data to construct a broad-band transmission spectrum and to search for expected strong absorbers (see Figs 8 and 9 ).
  • The activity of the parent star was taken into account by inflating the error bars of the measured R p /R * of each spectral bin assuming stellar variability of 0.3 per cent.
  • This number also takes into account the contribution of unocculted spots.
  • In the second approach, the authors removed the common-mode systematics from each spectral bin before fitting for residual trends with a parametrized model but with fewer parameters.
  • It was found that this method reduces the amplitude of the breathing systematics.

4 D I S C U S S I O N 4.1 A search for narrow-band spectral signatures

  • The authors performed a thorough search for narrow-band absorption features including the expected signatures of Na, K, Hα and Hβ.
  • The authors find only tentative evidence for sodium and potassium line cores with the largest significance levels of 1.2σ and 2.7σ in 50 and 170 Å bins, respectively.
  • The authors also searched but did not find evidence for Na and K pressure-broadened line wings.

4.2 Comparison to existing cloud-free theoretical models

  • The authors averaged the model within the transmission spectrum wavelength bins and fitted these theoretical values to the data with a single free parameter that controls their vertical position.
  • The models from Burrows et al. (2010) and Hothe authors & Burrows (2012) assume 1D dayside temperature-pressure (T-P) profile with stellar irradiation, in radiative, chemical, and hydrostatic equilibrium.
  • Chemical mixing ratios and corresponding opacities assume solar metallicity and local thermodynamical chemical equilibrium accounting for condensation with no ionization, using the opacity data base from Sharp & Burrows (2007) and the equilibrium chemical abundances from Burrows & Sharp (1999) and Burrows et al. (2001) .
  • Models were also generated both with and without the inclusion of TiO and VO opacities.
  • Because each of the aforementioned models started at wavelength λ = 3500 Å (due to a restriction of the employed opacity data base), the authors extrapolated the models in the range 2700-3500 Å to enable self-content comparisons to the bluest STIS measurements along with all models.

4.3 Interpreting the transmission spectrum with aerosols

  • The authors followed Lecavelier Des Etangs et al. (2008a) and performed a fit with a pure scattering model to the transmission spectrum with excluded sodium and potassium measurements as they can bias the result towards smaller slopes.
  • These measurements are in modest disagreement (at the 2.2σ confidence level) with the prediction of Rayleigh scattering from the ground-based result of Jordán et al. (2013) who found αT = −10670 ± 3015 K.
  • Note that their result significantly improves the error bar of the product αT due to the much broader wavelength coverage.

4.3.1 Rayleigh scattering

  • No secondary eclipse observations have been reported in the literature, which prevent an empiric comparison to the planet dayside brightness temperature.
  • This value is consistent with their results at the ∼1σ confidence level.

4.3.2 Mie scattering dust

  • Good fits to the data are possible assuming a significant opacity from aerosols, i.e. colloids of fine solid particles or liquid droplets in a gas.
  • Compared to the enstatite and forsterite, the alkali sulphides and chlorides are expected to be considerably less abundant in solarcomposition atmospheres.
  • A variety of different atmospheric models were tested assuming aerosols are the main spectral feature throughout the transmission spectrum.
  • In all fits, the authors excluded the radius measurements corresponding to the sodium and potassium spectral bins, as these are not predicted by the Mie scattering theoretical models.
  • The authors also estimated condensate grain sizes assuming a planetary temperature of 1145 K and fitting for a and the baseline radius.

4.3.3 High-altitude haze and lower altitude clouds

  • The authors also compared the transmission spectrum of WASP-6b with a suite of Fortney et al. (2010) atmospheric models with either an artificially added Rayleigh scattering component simulating a scattering haze, or an added 'cloud deck' produced by a grey flat line at specific altitude.
  • The authors performed comparison to these models as they can help better understand their observational results, as clouds and hazes can mask or completely block the expected atomic and molecular features in a transmission spectrum depending on the altitude distribution and particle size.
  • In the Solar system planets, the scaleheight of haze is often smaller than the scaleheight of the atmosphere itself, because of sedimentation and the structure of the vertical mixing.
  • That would result in a smaller Rayleigh slope for the haze signature in the transmission spectrum.
  • It would also add a gas-to-grain scaleheight ratio factor in equation 1 of Lecavelier Des Etangs et al. (2008a) , thus severing the relation between the measured slope and the temperature.

4.4 The atmosphere of WASP-6b

  • All of the cloud-free models have been found to produce unsatisfactory fits indicating low level of molecular absorption in their data.
  • Rayleigh scattering can potentially be produced by molecular hydrogen, a potential scenario explored for HD 209458b (Lecavelier Des Etangs et al. 2008b ).
  • Molecular hydrogen is disfavoured for WASP-6b by the lack of pressure-broadened line wings in its transmission spectrum and the low altitudes of any molecular features as apparent in Fig. 11 .
  • The first five best-fitting Mie scattering aerosol models are practically indistinguishable in quality between each other providing good representation of the observed transmission spectrum (Table 5 ).
  • The authors find that the particle sizes of all models are sub-micron, regardless the temperature.

4.4.1 Comparison to a ground-based result

  • Recently, Jordán et al. (2013) reported WASP-6b transmission spectrum with the Magellan IMACS instrument.
  • When fitting the light curves, the authors fixed the system parameters to their best-fitting values.
  • The authors find no evidence for the observed discrepancy to be due to the different system parameters (i.e. a/R * and i) nor the limbdarkening laws adopted.
  • 6 in Jordán et al. (2013) where white light curves are compared after subtracting the best-fitting systematics models, it can be seen that the final Magellan transit depth varies greatly between the different systematics models used.
  • Fig. 13 shows that the three STIS transit light curves obtained at three individual epochs are all in excellent agreement with each other, illustrating the measured HST transit depth is highly repeatable.

4.4.2 Comparison to HD 189733b

  • Complemented by the transmission spectra of HD 209458b and HD 189733b, the first results from the large HST programme cou- pled by Spitzer IRAC transit radii start building up an observational data base of hot-Jupiter atmospheres.
  • The similarities and differences in the spectra and system parameters of the studied objects could potentially indicate distinct categories of gas giant atmospheres.
  • Instead, absorption from the narrow cores of the Na I and K I doublets are evidenced for HD 189733b in the medium resolution STIS data and the revised ACS data from Huitson et al. (2012) and Pont et al. (2013) , respectively.
  • In addition, HD 189733b's transmission spectrum exhibits a steeper scattering slope and relatively higher altitude for the hazes on WASP-6b (Fig. 15 ).
  • A higher surface gravity means that a given pressure corresponds to less mass.

5 C O N C L U S I O N

  • With a broad-coverage optical transmission spectrum measured from HST and Spitzer broad-band transit spectrophotometry, WASP-6b joins the small but highly valuable family of hot-Jupiter exoplanets with atmospheric constraints.
  • The data point on the potassium spectral line is ∼2.7σ higher than the spectrum at that wavelength region.
  • The authors find no evidence of pressure-broadened wings which is in agreement with the haze model as clouds and hazes can reside at high altitudes which can significantly mask or even completely block the signatures of collisionally broadened wings.
  • Given the rather low temperature of WASP-6b, aerosol species are expected to originate from high-temperature condensates similar to these in the prototype HD 189733b.
  • Given their condensation temperatures and the results from their analysis, MgSiO 4 , KCl, MgSiO 3 and Na 2 S are candidate condensates.

Did you find this useful? Give us your feedback

Figures (17)

Content maybe subject to copyright    Report

HST hot-Jupiter transmission spectral survey: haze in the atmosphere
of WASP-6b
Nikolov, N., Sing, D. K., Burrows, A. S., Fortney, J. J., Henry, G. W., Pont, F., Ballester, G. E., Aigrain, S.,
Wilson, P. A., Huitson, C. M., Gibson, N. P., Désert, J-M., Lecavelier Des Etangs, A., Showman, A. P., Vidal-
Madjar, A., Wakeford, H. R., & Zahnle, K. (2015). HST hot-Jupiter transmission spectral survey: haze in the
atmosphere of WASP-6b.
Monthly Notices of the Royal Astronomical Society
,
447
(1), 463-478.
https://doi.org/10.1093/mnras/stu2433
Published in:
Monthly Notices of the Royal Astronomical Society
Document Version:
Publisher's PDF, also known as Version of record
Queen's University Belfast - Research Portal:
Link to publication record in Queen's University Belfast Research Portal
Publisher rights
© 2014 The Authors.
This work is made available online in accordance with the publisher’s policies. Please refer to any applicable terms of use of the publisher.
General rights
Copyright for the publications made accessible via the Queen's University Belfast Research Portal is retained by the author(s) and / or other
copyright owners and it is a condition of accessing these publications that users recognise and abide by the legal requirements associated
with these rights.
Take down policy
The Research Portal is Queen's institutional repository that provides access to Queen's research output. Every effort has been made to
ensure that content in the Research Portal does not infringe any person's rights, or applicable UK laws. If you discover content in the
Research Portal that you believe breaches copyright or violates any law, please contact openaccess@qub.ac.uk.
Download date:10. Aug. 2022

MNRAS 447, 463–478 (2015) doi:10.1093/mnras/stu2433
HST hot-Jupiter transmission spectral survey: haze in the
atmosphere of WASP-6b
N. Nikolov,
1
D. K. Sing,
1
A. S. Burrows,
2
J. J. Fortney,
3
G. W. Henry,
4
F. Pont,
1
G. E. Ballester,
5
S. Aigrain,
6
P. A. Wilson,
1,7
C. M. Huitson,
8
N. P. Gibson,
9
J.-M. D
´
esert,
8
A. Lecavelier des Etangs,
7
A. P. Showman,
5
A. Vidal-Madjar,
7
H. R. Wakeford
1
and K. Zahnle
10
1
Astrophysics Group, School of Physics, University of Exeter, Stocker Road, Exeter EX4 4QL, UK
2
Department of Astrophysical Sciences, Peyton Hall, Princeton University, Princeton, NJ 08544, USA
3
Department of Astronomy and Astrophysics, University of California, Santa Cruz, CA 95064, USA
4
Tennessee State University, 3500 John A. Merritt Blvd, PO Box 9501, Nashville, TN 37209, USA
5
Lunar and Planetary Laboratory, University of Arizona, Tucson, AZ 85721, USA
6
Department of Physics, University of Oxford, Denys Wilkinson Building, Keble Road, Oxford OX1 3RH, UK
7
CNRS, Institut dAstrophysique de Paris, UMR 7095, 98bis boulevard Arago, F-75014 Paris, France
8
CASA, Department of Astrophysical and Planetary Sciences, University of Colorado, 389-UCB, Boulder, CO 80309, USA
9
European Southern Observatory, Karl-Schwarzschild-Str. 2, D-85748 Garching bei M
¨
unchen, Germany
10
NASA Ames Research Center, Moffett Field, CA 94035, USA
Accepted 2014 November 17. Received 2014 November 17; in original form 2014 July 25
ABSTRACT
We report Hubble Space Telescope optical to near-infrared transmission spectroscopy of
the hot-Jupiter WASP-6b, measured with the Space Telescope Imaging Spectrograph and
Spitzer’s InfraRed Array Camera. The resulting spectrum covers the range 0.29–4.5 µm.
We find evidence for modest stellar activity of WASP-6 and take it into account in the
transmission spectrum. The overall main characteristic of the spectrum is an increasing radius
as a function of decreasing wavelength corresponding to a change of (R
p
/R
) = 0.0071
from 0.33 to 4.5 µm. The spectrum suggests an effective extinction cross-section with a power
law of index consistent with Rayleigh scattering, with temperatures of 973 ± 144 K at the
planetary terminator. We compare the transmission spectrum with hot-Jupiter atmospheric
models including condensate-free and aerosol-dominated models incorporating Mie theory.
While none of the clear-atmosphere models is found to be in good agreement with the data, we
find that the complete spectrum can be described by models that include significant opacity
from aerosols including Fe-poor Mg
2
SiO
4
, MgSiO
3
, KCl and Na
2
S dust condensates. WASP-
6b is the second planet after HD 189733b which has equilibrium temperatures near 1200 K
and shows prominent atmospheric scattering in the optical.
Key words: methods: observational techniques: spectroscopic planets and satellites: at-
mospheres stars: activity planets and satellites: individual: WASP-6b.
1 INTRODUCTION
Multiwavelength transit observations of hot Jupiters provide a
unique window to the chemistry and structure of the atmospheres
of these distant alien worlds. During planetary transits, a small
fraction of the stellar light is transmitted through the planetary at-
mosphere and signatures of atmospheric constituents are imprinted
on the stellar spectrum (Seager & Sasselov 2000). Transmission
E-mail: nikolov.nkn@gmail.com
spectroscopy, where the planet radius is measured as a function
of wavelength, has revolutionized our understanding of extrasolar
gas-giant atmospheres. Numerous studies from both space- and
ground-based observations have led to the characterization and
detection of atomic and molecular features as well as hazes and
clouds in the atmospheres of several hot Jupiters revealing a huge
diversity (Charbonneau et al. 2002; Grillmair et al. 2008; Pont
et al. 2008, 2013; Redfield et al. 2008; Snellen et al. 2008, 2010;
Sing et al. 2011, 2012, 2013; Brogi et al. 2012; Bean et al. 2013;
Birkby et al. 2013; Deming et al. 2013; Gibson et al. 2013a,b; Huit-
son et al. 2013; Wakeford et al. 2013; Nikolov et al. 2014; Stevenson
et al. 2014).
C
2014 The Authors
Published by Oxford University Press on behalf of the Royal Astronomical Society
Downloaded from https://academic.oup.com/mnras/article-abstract/447/1/463/989512
by Queen's University of Belfast user
on 08 February 2018

464 N. Nikolov et al.
Atmospheric hazes have currently been detected in hot-Jupiter
exoplanets over a wide range of atmospheric temperature regimes
including HD 189733b and WASP-12b representing the cool and
hot extremes for hot Jupiters at equilibrium temperatures around
1200 and 3000 K, respectively. Hazes are considered to originate
from dust condensation (forming aerosols in the atmospheres of
exoplanets) or a result from photochemistry (Marley et al. 2013).
Transmission spectroscopy studies may constrain the most prob-
able condensates responsible for these two planets, with iron-free
silicates (e.g. MgSiO
3
) and corundum (Al
2
O
3
) being two prime can-
didates for the above temperature regimes (Lecavelier Des Etangs
et al. 2008a;Singetal.2013). Cloud-free hot-Jupiter atmospheric
models predict sodium and potassium to be the dominant absorbing
features in optical transmission spectra (Seager & Sasselov 2000;
Fortney et al. 2005).
In this paper, we present new results for WASP-6b from a large
Hubble Space Telescope (HST) hot-Jupiter transmission spectral
survey comprising eight transiting planets. The ultimate aim of the
project is to explore the variety of hot-Jupiter atmospheres, i.e.
clear/hazy/cloudy, delve into the presence/lack of TiO/alkali and
other molecular features, probe the diversities between possible
sub-classes and perform comparative exoplanetology. Initial results
from four exoplanets have been presented so far in Huitson et al.
(2013) for WASP-19b, Wakeford et al. (2013)andNikolovetal.
(2013) for HAT-P-1b and Sing et al. (2013) for WASP-12b and
Sing et al. (2015) for WASP-31b. Already a wide diversity among
exoplanet atmospheres is observed. In this paper, we report new
HST optical transit observations with the Space Telescope Imaging
Spectrograph (STIS) instrument and combine them with Spitzer
InfraRed Array Camera (IRAC) broad-band transit photometry to
calculate a near-UV to near-infrared transmission spectrum, capable
of detecting atmospheric constituents. We describe the observations
in Section 2, present the light-curve analysis in Section 3, discuss
the results in Section 4 and conclude in Section 5.
1.1 The WASP-6b system
Discovered by the Wide Angle Search for Planets, WASP-6 b is
an inflated sub-Jupiter-mass transiting extrasolar planet orbiting a
moderately bright V = 11.9 solar-type, mildly metal-poor star (with
T
eff
= 5375 ± 65 K, log g = 4.61 ± 0.07 and [Fe/H]=−0.20 ± 0.09)
located in the southern part of the constellation Aquarius (Gillon
et al. 2009; Doyle et al. 2013). The planet is moving on an orbit with
a period of P 3.6 d and semimajor axis a 0.042 au. The sky-
projected angle between the stellar spin and the planetary orbital
axis has been determined through observations of the Rossiter–
McLaughlin effect by Gillon et al. (2009) indicating a good align-
ment (β = 11
+14
18
deg) that consequently favours a planet migration
scenario via the spin–orbit preserving tidal interactions with a pro-
toplanetary disc. Dragomir et al. (2011) refined the WASP-6 system
parameters and orbital ephemeris from a single ground-based transit
observation finding good agreement between their results and the
discovery paper. Although Gillon et al. (2009) claimed evidence for
non-zero orbital eccentricity, an analysis with new radial velocity
data from Husnoo et al. (2012) brought evidence for non-significant
eccentricity. Finally, Doyle et al. (2013) refined the spectroscopic
parameters of WASP-6b’s host star.
A cloud-free hot-Jupiter theoretical transmission spectrum pre-
dicts strong optical absorbers, dominated by Na
I and K I absorption
lines and H
2
Rayleigh scattering for a planetary system with the
physical properties of WASP-6b, including an effective planetary
temperature (assuming zero albedo and f = 1/4 heat redistribution)
Figure 1. Typical STIS G430L and G750L spectra (blue and red continuous
lines, respectively). A fringe-corrected spectrum (offset by 5 × 10
3
counts)
is portrayed above the uncorrected red spectrum for comparison.
of 1194 K and surface gravity g 8ms
2
(Fortney et al. 2008;
Burrows et al. 2010). The atmosphere of WASP-6b has recently
been probed in the optical regime from the ground by Jord
´
an et al.
(2013). The authors found evidence for an atmospheric haze char-
acterized by a decreasing apparent planetary size with wavelength
and no evidence for the pressure-broadened alkali features.
2 OBSERVATIONS AND CALIBRATIONS
2.1 HST STIS spectroscopy
We acquired low-resolution (R = λ/λ = 500–1040) HST STIS
spectra (Proposal ID GO-12473, P.I., D. Sing) during three transits
of WASP-6 b on
UT 2012 June 10 (visit 9) and 16 (visit 10) with
grating G430L (2.7 Å pixel
1
) and 2012 July 23 (visit 21) with
grating G750L (4.9 Å pixel
1
). When combined, the blue and
red STIS data provide complete wavelength coverage from 2900 to
10 300 Å with a small overlap region between them from 5240 to
5700 Å (Fig. 1). Each visit consisted of five 96 min orbits, during
which data collection was truncated with 45 min gaps due to Earth
occultations. Incorporating wide 52 × 2 arcsec slit to minimize slit
light losses and an integration time of 278 s, a total of 43 spectra
were obtained during each visit. Data acquisition overheads were
minimized by reading-out a reduced portion of the CCD with a size
of 1024 × 128. This observing strategy has proven to provide high
signal-to-noise ratio (S/N) spectra that are photometrically accurate
near the Poisson limit during a transit event (Brown et al. 2001; Sing
et al. 2011, 2013; Huitson et al. 2012). The three HST visits were
scheduled such that the second and third spacecraft orbits occurred
between the second and third contacts of the planetary transit in
order to provide good sampling of the planetary radius while three
orbits secured the stellar flux level before and after the transit.
The data reduction and analysis process is somewhat uniform
for the complete large HST programme and follows the general
methodology detailed in Huitson et al. (2013), Sing et al. (2013)
and Nikolov et al. (2014). The raw STIS data was reduced (bias-,
dark- and flat-corrected) using the latest version of the CALSTIS
1
pipeline and the relevant up-to-date calibration frames. Correction
of the significant fringing effect was performed for the G750L data
using contemporaneous fringe flat frame obtained at the end of the
1
CALSTIS comprises software tools developed for the calibration of STIS
data (Katsanis & McGrath 1998) inside the
IRAF (Image Reduction and
Analysis Facility; http://iraf.noao.edu/) environment.
MNRAS 447, 463–478 (2015)
Downloaded from https://academic.oup.com/mnras/article-abstract/447/1/463/989512
by Queen's University of Belfast user
on 08 February 2018

HST STIS transit spectroscopy of WASP-6b 465
observing sequence and the procedure detailed in (Goudfrooij &
Christensen 1998, Fig. 1).
Due to the relatively long STIS integration time (i.e. 287 s), the
data contain large number of cosmic ray events, which were iden-
tified and removed following the procedures described in Nikolov
et al. (2013). It was found that the total number of pixels affected
by cosmic ray events comprise 4 per cent of the total number of
pixels of each 2D spectrum. In addition, we also corrected all pixels
identified by CALSTIS as ‘bad’ with the same procedure, which
together with the cosmic ray identified pixels resulted in a total
of 11 per cent interpolated pixels.
Spectral extractions were performed in
IRAF employing the APALL
procedure using the calibrated .flt science files after fringe and
cosmic ray correction. We performed spectral extraction with aper-
ture sizes in the range 6.0–18.0 pixels with a step of 0.2. The best
aperture for each grating was selected based on the resulting low-
est light-curve residual scatter after fitting the white light curves
(see Section 3 for details). We found that aperture sizes 12.0, 12.2
and 8.6 pixels satisfy this criterium for visits 9, 10 and 21, respec-
tively. We then placed the extracted spectra to a common Doppler
corrected rest frame through cross-correlation to prevent sub-pixel
wavelength shifts in the dispersion direction. A wavelength solution
was obtained from the x1d files from CALSTIS. The STIS spectra
were then used to extract both white-light spectrophotometric time
series and custom wavelength bands after integrating the appropri-
ate flux from each bandpass.
The raw STIS light curves exhibit instrumental systematics s imi-
lar to those described by Gilliland, Goudfrooij & Kimble (1999)and
Brown et al. (2001). In summary, the major source of the systemat-
ics is related with the orbital motion of the telescope. In particular,
the HST focus is known to experience quite noticeable variations
on the spacecraft orbital time-scale, which are attributed to ther-
mal contraction/expansion (often referred to as the ‘breathing ef-
fect’) as the spacecraft warms up during its orbital day and cools
down during orbital night (Hasan & Bely 1993, 1994; Suchkov &
Hershey 1998). We take into account the systematics associated
with the telescope temperature variations in the transit light-curve
fits by fitting a baseline function depending on various parameters
(Section 3.0.1).
2.2 Spitzer IRAC data
Photometric data were collected for WASP-6 during two transits
of its planet on
UT 2013 January 14 and 21 with the Spitzer space
telescope (Werner et al. 2004) employing, respectively, the 4.5 and
3.6 µm channels of the IRAC (Fazio et al. 2004)aspartofpro-
gramme 90092 (P.I. J.-M. Desert). Both observations started shortly
before ingress with effective integration times of 1.92 s per image
and were terminated 2.2 h after egress resulting in 8 320 images
(see Fig. 2), which were calibrated by the Spitzer pipeline (version
S19.1.0) and are available in the form of Basic Calibrated Data
(.bcd)files.
After organizing the data, we converted the images from flux in
mega-Jansky per steradian (MJy sr
1
) to photon counts (i.e. elec-
trons) using the information provided in the
FITS headers. In partic-
ular, we multiplied the images by the gain and individual exposure
time (
FITS header key words SAMPTIME and GAIN) and divided
by the flux conversion factor (FLUXCONV). Timing of each im-
age was computed using the UTC-based Barycentric Julian Date
(BJD
UTC
) from the FITS header key word BMJD_OBS, transform-
ing these time stamps into BJD based on the BJD terrestrial time
(TT) using the following conversion BJD
TT
BJD
UTC
+ 66.184 s
Figure 2. Warm Spitzer IRAC 3.6 and 4.5µm photometry (left and right, respectively). Top panels: raw flux and the best-fitting transit and systematics model;
middle panels: detrended light curves and the best-fitting transit models and binned by 8 min; lower panels: light-curve residuals.
MNRAS 447, 463–478 (2015)
Downloaded from https://academic.oup.com/mnras/article-abstract/447/1/463/989512
by Queen's University of Belfast user
on 08 February 2018

466 N. Nikolov et al.
(Eastman, Siverd & Gaudi 2010). This conversion is preferable
as leap seconds are occasionally added to the BJD
UTC
standard
(Knutson et al. 2012; Todorov et al. 2013).
We performed outlier filtering for hot (energetic) or lower pixels
in the data by following each pixel through time. This task was
performed in two steps, first flagging all pixels with intensity 8σ or
more away from the median value computed from the five preced-
ing and five following images. The values of these flagged pixels
were replaced with the local median value. In the second pass, we
flagged and replaced outliers above the 4σ level, following the same
procedure. The total fraction of corrected pixels was 0.26 per cent
for the 3.6 µm and 0.06 per cent for the 4.5 µm channel.
We estimated and subtracted the background flux from each im-
age of the time series. To do this, we performed an iterative 3σ
outlier clipping for each image to remove the pixels with values
associated with the stellar point spread function (PSF), background
stars or hot pixels, created a histogram from the remaining pixels
and fitted a Gaussian to determine the sky background.
We measured the position of the star on the detector in each im-
age incorporating the flux-weighted centroiding method
2
using the
background-subtracted pixels from each image for a circular region
with radius 3 pixels centred on the approximate position of the star.
While we could perform PSF centroiding using alternative methods
such as fitting a two-dimensional Gaussian function to the stellar
image, previous experiences with warm Spitzer photometry showed
that the flux-weighted centroiding method is either equivalent or su-
perior (Beerer et al. 2011;Knutsonetal.2012; Lewis et al. 2013).
The variation of the x and y positions of the PSF on the detectors
were measured to be 0.20 and 0.21 for the 3.6 µ m channel and 0.89
and 0.8 pixels for the 4.5 µm channel.
We extract photometric measurements from our data following
two methods. First, aperture photometry was performed with
IDL
routine APER using circular apertures ranging in radius from 1.5 to
3.5 pixels in increments of 0.1. We filtered the resulting light curves
for 5σ outliers with a width of 20 data points.
We also performed photometry with time-variable aperture with
a size scaled by the value of a quantity known as the noise pixel
parameter (Mighell 2005; Knutson et al. 2012), which is described
in section 2.2.2 of the IRAC instrument handbook and has been
used in previous exoplanet studies to improve the results of warm
Spitzer photometry (see e.g. Charbonneau et al. 2008; Knutson
et al. 2012; Todorov et al. 2013; O’Rourke et al. 2014). The noise
parameter depends on the full width at half-maximum of the stellar
PSF squared and is defined as
˜
β =
(
i
I
i
)
2
i
I
2
i
, (1)
where I
i
is the intensity detected by the ith pixel. We use each
image to measure the noise pixel parameter applying an aperture
with radius of 4 pixels, ensuring that each pixel is considered should
the border of the aperture cross that pixel. We extracted source fluxes
from both channels with the aperture radii following the relation
r =
˜
βa
0
+ a
1
, (2)
where a
0
and a
1
were varied in the ranges 0.8 to 1.2 and 0.4 to
+0.4 with step 0.1.
The best results from both photometric methods were identified
by examining both the residual root mean square (rms) after fitting
2
As implemented in the IDL box centroider.pro, provided on the Spitzer
homepage: http://irsa.ipac.caltech.edu.
the light curves from each channel as well as the white and red
noise components measured with the wavelet technique detailed in
Carter & Winn (2009). The second photometric method resulted in
the lowest white and random red noise correlated with data points
co-added in time or spectral sampling for the 3.6 µmdatawith
a
0
= 0.9 and a
1
=−0.1. For the 4.5 µm data, the first method gave
better results with aperture radius 2.4 pixels and sky annulus defined
between radii 2.40 and 6.48 pixels.
Finally, we also performed a visual inspection of the Spitzer
IRAC images for obvious stellar companions by aligning and stack-
ing 3000 of the available images in both channels, finding no ev-
idence for bright stellar sources within the 38 × 38 arcsec
2
field of
view.
2.3 Stellar variability monitoring
Stellar activity can complicate the interpretation of exoplanet trans-
mission spectra, especially when multi-instrument multi-epoch data
sets are combined (Pont et al. 2013). As the star rotates, active re-
gions on its surface enter into and out of view, causing the measured
flux to exhibit a quasi-periodic variability that introduces variation
of the measured transit depth (as R
p
/R
). This becomes particularly
important when transit observations made over several months or
years are combined to construct an exoplanet transmission spec-
trum, as in our HST/Spitzer study.
We obtained high-resolution spectra of WASP-6 from the pub-
licly available HARPS data to search for evidence of chromospheric
activity. All spectra in the wavelength region of the Ca
II H&K
lines show evidence for emission, implying a moderate stellar ac-
tivity compared to the same region of HARPS spectra obtained
for HD189733. An examination of the HST and Spitzer white and
spectral light curves show no evidence for spot crossings.
We obtained nightly photometry of WASP-6 to monitor and char-
acterize the stellar activity over the past three observing seasons with
the Tennessee State University Celestron 14-inch (C14) automated
imaging telescope (AIT) located at Fairborn Observatory in Arizona
(Henry 1999). The AIT uses an SBIG STL-1001E CCD camera
and exposes through a Cousins R filter. Each nightly observation
of WASP-6 consists of 4–10 consecutive exposures that include
several comparison stars in the same field of view. The individual
nightly frames were co-added and reduced to differential magni-
tudes (i.e. WASP-6 minus the mean brightness of nine constant
comparison stars). Each nightly observation has been corrected for
bias, flat-fielding, pier-side offset and for differential atmospheric
extinction.
A total of 258 nightly observations (excluding a few isolated
transit observations) were acquired in five groups during three ob-
serving seasons between 2011 September and 2014 January. The
five groups are plotted in Fig. 3, where groups 2 through 5 have
been normalized to have the same mean magnitude as the first
group. The standard deviations of the five groups with respect to
their corresponding means are 0.0032, 0.0065, 0.0034, 0.0078 and
0.0053 mag, respectively. The precision of a single measurement in
good photometric conditions with C14 is typically 0.002–0.003 mag
(see, e.g. Sing et al. 2013). The standard deviations of groups 2, 4
and 5 significantly exceed the measurement precision and thus in-
dicate low-level activity due to star spots. A periodogram analysis
has been performed of each data set with trial frequencies ranging
between 0.005 and 0.95 c d
1
, corresponding to a period range be-
tween 1 and 200 d. Period analysis of group 5 gives the clearest
indication of rotational modulation with a period of 23.6 ± 0.5 d
and a peak-to-peak amplitude of 0.01 mag. We take this period to
MNRAS 447, 463–478 (2015)
Downloaded from https://academic.oup.com/mnras/article-abstract/447/1/463/989512
by Queen's University of Belfast user
on 08 February 2018

Citations
More filters
Journal ArticleDOI
07 Jan 2016-Nature
TL;DR: The difference between the planetary radius measured at optical and infrared wavelengths is an effective metric for distinguishing different atmosphere types, so that strong water absorption lines are seen in clear-atmosphere planets and the weakest features are associated with clouds and hazes.
Abstract: Thousands of transiting exoplanets have been discovered, but spectral analysis of their atmospheres has so far been dominated by a small number of exoplanets and data spanning relatively narrow wavelength ranges (such as 1.1-1.7 micrometres). Recent studies show that some hot-Jupiter exoplanets have much weaker water absorption features in their near-infrared spectra than predicted. The low amplitude of water signatures could be explained by very low water abundances, which may be a sign that water was depleted in the protoplanetary disk at the planet's formation location, but it is unclear whether this level of depletion can actually occur. Alternatively, these weak signals could be the result of obscuration by clouds or hazes, as found in some optical spectra. Here we report results from a comparative study of ten hot Jupiters covering the wavelength range 0.3-5 micrometres, which allows us to resolve both the optical scattering and infrared molecular absorption spectroscopically. Our results reveal a diverse group of hot Jupiters that exhibit a continuum from clear to cloudy atmospheres. We find that the difference between the planetary radius measured at optical and infrared wavelengths is an effective metric for distinguishing different atmosphere types. The difference correlates with the spectral strength of water, so that strong water absorption lines are seen in clear-atmosphere planets and the weakest features are associated with clouds and hazes. This result strongly suggests that primordial water depletion during formation is unlikely and that clouds and hazes are the cause of weaker spectral signatures.

955 citations

Journal ArticleDOI
TL;DR: In this article, the authors present a consistent optimal estimation retrieval analysis of 10 hot Jupiter exoplanets, each with transmission spectral data spanning the visible to near-infrared wavelength range.
Abstract: We present a consistent optimal estimation retrieval analysis of 10 hot Jupiter exoplanets, each with transmission spectral data spanning the visible to near-infrared wavelength range. Using the NEMESIS radiative transfer and retrieval tool, we calculate a range of possible atmospheric states for WASP-6b, WASP-12b, WASP-17b, WASP-19b, WASP-31b, WASP-39b, HD 189733b, HD 209458b, HAT-P-1b, and HAT-P-12b. We find that the spectra of all 10 planets are consistent with the presence of some atmospheric aerosol; WASP-6b, WASP-12b, WASP-17b, WASP-19b, HD 189733b, and HAT-P-12b are all fit best by Rayleigh scattering aerosols, whereas WASP-31b, WASP-39b and HD 209458b are better represented by a gray cloud model. HAT-P-1b has solutions that fall into both categories. WASP-6b, HAT-P-12b, HD 189733b, and WASP-12b must have aerosol extending to low atmospheric pressures (below 0.1 mbar). In general, planets with equilibrium temperatures between 1300 and 1700 K are best represented by deeper, gray cloud layers, whereas cooler or hotter planets are better fit using high Rayleigh scattering aerosol. We find little evidence for the presence of molecular absorbers other than H2O. Retrieval methods can provide a consistent picture across a range of hot Jupiter atmospheres with existing data, and will be a powerful tool for the interpretation of James Webb Space Telescope observations.

235 citations


Cites background from "HST hot-Jupiter transmission spectr..."

  • ...…details of observations, see Pont et al. 2013; Line et al. 2013a; Huitson et al. 2013; Mandell et al. 2013; Sing et al. 2013, 2015; Wakeford et al. 2013; Nikolov et al. 2014, 2015; McCullough et al. 2014) represent a variety of atmospheres, interpreted as a continuum of clear to cloudy conditions....

    [...]

Journal ArticleDOI
12 May 2017-Science
TL;DR: HAT-P-26b’s atmosphere is primordial and obtained its gaseous envelope late in its disk lifetime, with little contamination from metal-rich planetesimals, according to a detailed atmospheric study of the transiting Neptune-mass exoplanet.
Abstract: A correlation between giant-planet mass and atmospheric heavy elemental abundance was first noted in the past century from observations of planets in our own Solar System and has served as a cornerstone of planet-formation theory. Using data from the Hubble and Spitzer Space Telescopes from 0.5 to 5 micrometers, we conducted a detailed atmospheric study of the transiting Neptune-mass exoplanet HAT-P-26b. We detected prominent H 2 O absorption bands with a maximum base-to-peak amplitude of 525 parts per million in the transmission spectrum. Using the water abundance as a proxy for metallicity, we measured HAT-P-26b’s atmospheric heavy element content ( 4.8 − 4.0 + 21.5 times solar). This likely indicates that HAT-P-26b’s atmosphere is primordial and obtained its gaseous envelope late in its disk lifetime, with little contamination from metal-rich planetesimals.

195 citations

Journal ArticleDOI
TL;DR: In this paper, the authors used ground-based photometric light curves from the NASA/ESA Hubble Space Telescope to obtain the ground-truth light curves of the WASP-121 discovery team.
Abstract: The authors would like to thank the referee for their prompt and thoughtful review. This work is based on observations made with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-26555. These observations are associated with program GO-14468. The authors are grateful to the WASP-121 discovery team for generously providing the ground-based photometric light curves. The research leading to these results received funding from the European Research Council under the European Union Seventh Framework Program (FP7/2007-2013) ERC grant agreement no. 336792. H.R.W. acknowledges support by an appointment to the NASA Postdoctoral Program at Goddard Space Flight Center, administered by ORAU and USRA through a contract with NASA. N.P.G. gratefully acknowledges support from the Royal Society in the form of a University Research Fellowship.

193 citations


Cites background from "HST hot-Jupiter transmission spectr..."

  • ...Transmission spectra at optical wavelengths show evidence for Rayleigh scattering by molecular hydrogen and high-altitude aerosols, and absorption by alkali metals (e.g. Charbonneau et al. 2002; Pont et al. 2008; Sing et al. 2013, 2015; Nikolov et al. 2014, 2015)....

    [...]

Journal ArticleDOI
TL;DR: The theory of atmospheric chemistry in all classes of exoplanets discovered to date is reviewed, from highly irradiated gas giants, ice giants, and super-Earths, to directly imaged giant planets at large orbital separations.
Abstract: Characterizing the atmospheres of extrasolar planets is the new frontier in exoplanetary science. The last two decades of exoplanet discoveries have revealed that exoplanets are very common and extremely diverse in their orbital and bulk properties. We now enter a new era as we begin to investigate the chemical diversity of exoplanets, their atmospheric and interior processes, and their formation conditions. Recent developments in the field have led to unprecedented advancements in our understanding of atmospheric chemistry of exoplanets and the implications for their formation conditions. We review these developments in the present work. We review in detail the theory of atmospheric chemistry in all classes of exoplanets discovered to date, from highly irradiated gas giants, ice giants, and super-Earths, to directly imaged giant planets at large orbital separations. We then review the observational detections of chemical species in exoplanetary atmospheres of these various types using different methods, including transit spectroscopy, Doppler spectroscopy, and direct imaging. In addition to chemical detections, we discuss the advances in determining chemical abundances in these atmospheres and how such abundances are being used to constrain exoplanetary formation conditions and migration mechanisms. Finally, we review recent theoretical work on the atmospheres of habitable exoplanets, followed by a discussion of future outlook of the field.

182 citations


Cites result from "HST hot-Jupiter transmission spectr..."

  • ...Similar inferences have been made for several other hot Jupiters in recent years (Sing et al. 2013, 2015; Nikolov et al. 2015)....

    [...]

References
More filters
Journal ArticleDOI
TL;DR: In this paper, the problem of selecting one of a number of models of different dimensions is treated by finding its Bayes solution, and evaluating the leading terms of its asymptotic expansion.
Abstract: The problem of selecting one of a number of models of different dimensions is treated by finding its Bayes solution, and evaluating the leading terms of its asymptotic expansion. These terms are a valid large-sample criterion beyond the Bayesian context, since they do not depend on the a priori distribution.

38,681 citations

01 Jan 2005
TL;DR: In this paper, the problem of selecting one of a number of models of different dimensions is treated by finding its Bayes solution, and evaluating the leading terms of its asymptotic expansion.
Abstract: The problem of selecting one of a number of models of different dimensions is treated by finding its Bayes solution, and evaluating the leading terms of its asymptotic expansion. These terms are a valid large-sample criterion beyond the Bayesian context, since they do not depend on the a priori distribution.

36,760 citations


"HST hot-Jupiter transmission spectr..." refers methods in this paper

  • ...Such models were found statistically unjustified based on the Bayesian Information Criterion (BIC; Schwarz 1978)....

    [...]

Journal ArticleDOI
TL;DR: In this article, solar photospheric and meteoritic CI chondrite abundance determinations for all elements are summarized and the best currently available photosphere abundances are selected, including the meteoritic and solar abundances of a few elements (e.g., noble gases, beryllium, boron, phosphorous, sulfur).
Abstract: Solar photospheric and meteoritic CI chondrite abundance determinations for all elements are summarized and the best currently available photospheric abundances are selected. The meteoritic and solar abundances of a few elements (e.g., noble gases, beryllium, boron, phosphorous, sulfur) are discussed in detail. The photospheric abundances give mass fractions of hydrogen (X ¼ 0:7491), helium (Y ¼ 0:2377), and heavy elements (Z ¼ 0:0133), leading to Z=X ¼ 0:0177, which is lower than the widely used Z=X ¼ 0:0245 from previous compilations. Recent results from standard solar models considering helium and heavy-element settling imply that photospheric abundances and mass fractions are not equal to protosolar abundances (representative of solar system abundances). Protosolar elemental and isotopic abundances are derived from photospheric abundances by considering settling effects. Derived protosolar mass fractions are X0 ¼ 0:7110, Y0 ¼ 0:2741, and Z0 ¼ 0:0149. The solar system and photospheric abundance tables are used to compute self-consistent sets of condensation temperatures for all elements. Subject headings: astrochemistry — meteors, meteoroids — solar system: formation — Sun: abundances — Sun: photosphere

4,305 citations


"HST hot-Jupiter transmission spectr..." refers background in this paper

  • ...…adopted here due to the lack of observations on the planet dayside spectrum. c© 2013 RAS, MNRAS 000, 1–?? of WASP-6b the prime condensate candidate materials include MgSiO3 (enstatite), Mg2SiO4 (forsterite), sulphides and chlorides of the alkali metals Na and K (Lodders 2003; Morley et al. 2012)....

    [...]

Journal ArticleDOI
TL;DR: The Infrared Array Camera (IRAC) is one of the three focal plane instruments in the Spitzer Space Telescope as discussed by the authors, which is a four-channel camera that obtains simultaneous broadband images at 3.6, 4.5, 5.8, and 8.0 microns.
Abstract: The Infrared Array Camera (IRAC) is one of three focal plane instruments in the Spitzer Space Telescope. IRAC is a four-channel camera that obtains simultaneous broad-band images at 3.6, 4.5, 5.8, and 8.0 microns. Two nearly adjacent 5.2x5.2 arcmin fields of view in the focal plane are viewed by the four channels in pairs (3.6 and 5.8 microns; 4.5 and 8 microns). All four detector arrays in the camera are 256x256 pixels in size, with the two shorter wavelength channels using InSb and the two longer wavelength channels using Si:As IBC detectors. IRAC is a powerful survey instrument because of its high sensitivity, large field of view, and four-color imaging. This paper summarizes the in-flight scientific, technical, and operational performance of IRAC.

3,696 citations


"HST hot-Jupiter transmission spectr..." refers methods in this paper

  • ...…data was collected for WASP-6 during two transits of its planet on UT 2013 January 14 and 21 with the Spitzer space telescope (Werner et al. 2004) employing respectively the 4.5 and 3.6µm channels of the Infrared Array Camera (IRAC, Fazio et al. 2004) as part of program 90092 (P.I. J.-M. Desert)....

    [...]

Journal ArticleDOI
TL;DR: The Infrared Array Camera (IRAC) is one of three focal plane instruments on the Spitzer Space Telescope as mentioned in this paper, which is a four-channel camera that obtains simultaneous broadband images at 3.6, 4.5, 5.8, and 8.0 m.
Abstract: The Infrared Array Camera (IRAC) is one of three focal plane instruments on the Spitzer Space Telescope. IRAC is a four-channel camera that obtains simultaneous broadband images at 3.6, 4.5, 5.8, and 8.0 � m. Two nearly adjacent 5A2 ; 5A2 fields of view in the focal plane are viewed by the four channels in pairs (3.6 and 5.8 � m; 4.5 and 8 � m). All four detector arrays in the camera are 256 ; 256 pixels in size, with the two shorter wavelength channels using InSb and the two longer wavelength channels using Si:As IBC detectors. IRAC is a powerful survey instrument because of its high sensitivity, large field of view, and four-color imaging. This paper summarizes the in-flight scientific, technical, and operational performance of IRAC.

3,567 citations