Archival Hubble Space Telescope Observations Research Articles

HST SHEL: Enabling Comparative Exoplanetology with HST/STIS

The Hubble Space Telescope (HST) has been our most prolific tool to study exoplanet atmospheres. As the age of JWST begins, there are a wealth of HST archival data that are useful to strengthen our inferences from JWST. Notably, HST/Space Telescope Imaging Spectrograph (STIS), with its 0.3–1 μm wavelength coverage, extends past JWST’s 0.6 μm wavelength cutoff and holds an abundance of potential information: alkali (Na, K) and molecular (TiO, VO) species opacities, aerosol information, and the presence of stellar contamination. However, time-series observations with HST suffer from significant instrumental systematics and can be highly dependent on choices made during the transit fitting process. This makes comparing transmission spectra of planets with different data reduction methodologies challenging, as it is difficult to discern whether an observed trend is caused by differences in data reduction or underlying physical processes. Here we present the Sculpting Hubble’s Exoplanet Legacy (SHEL) program, which aims to build a consistent data reduction and light-curve analysis methodology and associated database of transmission spectra from archival HST observations. In this paper, we present the SHEL analysis framework for HST/STIS and its low-resolution spectroscopy modes, G430L and G750L. We apply our methodology to four notable hot Jupiters, WASP-39 b, WASP-121 b, WASP-69 b, and WASP-17 b, and use these examples to discuss nuances behind analysis with HST/STIS. Our results for WASP-39 b, WASP-121 b, and WASP-17 b are consistent with past publications, but our analysis of WASP-69 b differs and shows evidence of either a strong scattering slope or stellar contamination. The data reduction pipeline and tutorials are available on Github and Zenodo.

Recurring outbursts of the supernova impostor AT 2016blu in NGC 4559

ABSTRACT We present the first photometric analysis of the supernova (SN) impostor AT 2016blu in NGC 4559. This transient was discovered by the Lick Observatory Supernova Search in 2012 and has continued its outbursts since then. Optical and infrared photometry of AT 2016blu reveals at least 19 outbursts in 2012–2022. Similar photometry from 1999 to 2009 shows no outbursts, indicating that the star was relatively stable in the decade before discovery. Archival Hubble Space Telescope observations suggest that the progenitor had a minimum initial mass of $M \gtrsim 33$ M⊙ and a luminosity of $L \gtrsim 10^{5.7}$ L⊙. AT 2016blu’s outbursts show irregular variability with multiple closely spaced peaks having typical amplitudes of 1–2 mag and durations of 1–4 weeks. While individual outbursts have irregular light curves, concentrations of these peaks recur with a period of ∼113 ± 2 d. Based on this period, we predict upcoming outbursts in 2023 and 2024. AT 2016blu shares similarities with SN 2000ch in NGC 3432, where outbursts may arise from periastron encounters in an eccentric binary containing a luminous blue variable (LBV). We propose that AT 2016blu’s outbursts are also driven by interactions that intensify around periastron in an eccentric system. Intrinsic variability of the LBV-like primary star may cause different intensity and duration of binary interaction at each periastron passage. AT 2016blu also resembles the periastron encounters of η Carinae prior to its Great Eruption and the erratic pre-SN eruptions of SN 2009ip. This similarity and the onset of eruptions in the past decade hint that AT 2016blu may also be headed for a catastrophe, making it a target of great interest.

On the Contribution of the X-Ray Source to the Extended Nebular He ii Emission in IZW18

Abstract Nebular He ii emission implies the presence of energetic photons (E ≥ 54 eV). Despite the great deal of effort dedicated to understanding He ii ionization, its origin has remained mysterious, particularly in metal-deficient star-forming (SF) galaxies. Unfolding He ii-emitting, metal-poor starbursts at z ∼ 0 can yield insight into the powerful ionization processes occurring in the primordial universe. Here we present a new study on the effects that X-ray sources have on the He ii ionization in the extremely metal-poor galaxy IZw18 (Z ∼ 3% Z ⊙), whose X-ray emission is dominated by a single high-mass X-ray binary (HMXB). This study uses optical integral field spectroscopy, archival Hubble Space Telescope observations, and all of the X-ray data sets publicly available for IZw18. We investigate the time-variability of the IZw18 HMXB for the first time; its emission shows small variations on timescales from days to decades. The best-fit models for the HMXB X-ray spectra cannot reproduce the observed He ii ionization budget of IZw18, nor can recent photoionization models that combine the spectra of both very low metallicity massive stars and the emission from HMXB. We also find that the IZw18 HMXB and the He ii-emission peak are spatially displaced at a projected distance of ≃200 pc. These results reduce the relevance of X-ray photons as the dominant He ii ionizing mode in IZw18, which leaves uncertain what process is responsible for the bulk of its He ii ionization. This is in line with recent work discarding X-ray binaries as the main source responsible for He ii ionization in SF galaxies.

Constraining VLBI−optical offsets in high redshift galaxies using strong gravitational lensing

ABSTRACT We present a multiwavelength analysis of two highly magnified strong gravitationally lensed galaxies, CLASS B0712+472 and CLASS B1608+656, at redshifts 1.34 and 1.394, respectively, using new VLBI (very long baseline interferometry) and archival Hubble Space Telescope observations. We reconstruct the positions of the radio and optical emissions with their uncertainties using Monte Carlo sampling. We find that in CLASS B0712+472 the optical and radio emissions are co-spatial within 2 ± 5 mas (17 ± 42 pc at redshift of 1.34). But, in CLASS B1608+656, we reconstruct an optical–radio offset of 25 ± 16 mas (214 ± 137 pc at redshift of 1.394), among the smallest offsets measured for an AGN (active galactic nucleus) at such high redshift. The spectral features indicate that CLASS B1608+656 is a post-merger galaxy, which, in combination with the optical–VLBI offset reported here, makes CLASS B1608+656 a promising candidate for a high- z offset–AGN. Furthermore, the milliarcsecond angular resolution of the VLBI observations combined with the precise lens models allow us to spatially locate the radio emission at 0.05 mas precision (0.4 pc) in CLASS B0712+472, and 0.009 mas precision (0.08 pc) in CLASS B1608+656. The search for optical–radio offsets in high redshift galaxies will be eased by the upcoming synoptic all-sky surveys, including Extremely Large Telescope and Square Kilometre Array, which are expected to find ∼105 strongly lensed galaxies, opening an era of large strong lensing samples observed at high angular resolution.

Searching for the Donor Stars of ULX Pulsars

Abstract We report on our search for the optical counterparts of two ultraluminous X-ray pulsars with known orbital periods, M82 X-2 and NGC 5907 X-1, in new and archival Hubble Space Telescope observations, in an effort to characterize the donor stars in these systems. We detect five near-infrared sources consistent with the position of M82 X-2 that are too bright to be single stars. We also detect seven sources in the WFC3/UVIS F336W image whose photometry matches that of 10–15 M ⊙ stars turning off the main sequence. Such stars have densities consistent with the properties of the donor star of M82 X-2 as inferred from X-ray timing analysis, although it is also possible that the donor is a lower-mass star below our detection limit or that there is a significant contribution from the accretion disk to the optical emission. We detect three candidate counterparts to NGC 5907 X-1 in the near-infrared. All of these are too bright to be the donor star of the ultraluminous X-ray source (ULX), which, based on its orbital period, is a red giant. The high background at the location of NGC 5907 X-1 precludes us from detecting this expected donor star. The recently discovered NGC 5907 ULX-2 also falls within the field of view of the near-infrared imaging; we detect four sources in the error circle, with photometry that matches asymptotic giant branch stars. The star suggested to be the counterpart of NGC 5907 ULX-2 by Pintore et al. falls outside our 2σ error circle.

Inflows, Outflows, and a Giant Donor in the Remarkable Recurrent Nova M31N 2008-12a?—Hubble Space Telescope Photometry of the 2015 Eruption

Abstract The recurrent nova M31N 2008-12a experiences annual eruptions, contains a near-Chandrasekhar-mass white dwarf, and has the largest mass accretion rate in any nova system. In this paper, we present Hubble Space Telescope (HST) WFC3/UVIS photometry of the late decline of the 2015 eruption. We couple these new data with archival HST observations of the quiescent system and Keck spectroscopy of the 2014 eruption. The late-time photometry reveals a rapid decline to a minimum luminosity state, before a possible recovery/rebrightening in the run up to the next eruption. Comparison with accretion disk models supports the survival of the accretion disk during the eruptions, and uncovers a quiescent disk mass accretion rate of the order of 10 − 6 M ⊙ yr − 1 , which may rise beyond 10 − 5 M ⊙ yr − 1 during the super-soft source phase—both of which could be problematic for a number of well-established nova eruption models. Such large accretion rates, close to the Eddington limit, might be expected to be accompanied by additional mass loss from the disk through a wind and even through collimated outflows. The archival HST observations, combined with the disk modeling, provide the first constraints on the mass donor: L donor = 103 − 11 + 12 L ⊙ , R donor = 14.14 − 0.47 + 0.46 R ⊙ , and T eff , donor = 4890 ± 110 K, which may be consistent with an irradiated M31 red-clump star. Such a donor would require a system orbital period ≳ 5 days. Our updated analysis predicts that the M31N 2008-12a WD could reach the Chandrasekhar mass in < 20 kyr.

THE TYPE IIb SUPERNOVA 2013df AND ITS COOL SUPERGIANT PROGENITOR

We have obtained early-time photometry and spectroscopy of Supernova (SN) 2013df in NGC 4414. The SN is clearly of Type IIb, with notable similarities to SN 1993J. From its luminosity at secondary maximum light, it appears that less $^{56}$Ni ($\lesssim 0.06\ M_{\odot}$) was synthesized in the SN 2013df explosion than was the case for the SNe IIb 1993J, 2008ax, and 2011dh. Based on a comparison of the light curves, the SN 2013df progenitor must have been more extended in radius prior to explosion than the progenitor of SN 1993J. The total extinction for SN 2013df is estimated to be $A_V=0.30$ mag. The metallicity at the SN location is likely to be solar. We have conducted Hubble Space Telescope (HST) Target of Opportunity observations of the SN with the Wide Field Camera 3, and from a precise comparison of these new observations to archival HST observations of the host galaxy obtained 14 years prior to explosion, we have identified the progenitor of SN 2013df to be a yellow supergiant, somewhat hotter than a red supergiant progenitor for a normal Type II-Plateau SN. From its observed spectral energy distribution, assuming that the light is dominated by one star, the progenitor had effective temperature $T_{\rm eff} = 4250 \pm 100$ K and a bolometric luminosity $L_{\rm bol}=10^{4.94 \pm 0.06}\ L_{\odot}$. This leads to an effective radius $R_{\rm eff} = 545 \pm 65\ R_{\odot}$. The star likely had an initial mass in the range of 13 to 17 $M_{\odot}$; however, if it was a member of an interacting binary system, detailed modeling of the system is required to estimate this mass more accurately. The progenitor star of SN 2013df appears to have been relatively similar to the progenitor of SN 1993J.

Hα photometry of low-mass stars in 47 Tucanae: chromospheric activity and exotica

We have used archival Hubble Space Telescope observations obtained with the Advanced Camera for Surveys to study the Hα emission properties of main-sequence stars in the globular cluster 47 Tucanae. Using a combination of multiband observations in the F606W, F814W and F658N bands, we search for stars showing Hα excess emission. An accurate photometric measurement of their Hα equivalent width allows us to identify objects with large Hα emission, which we attribute to mass accretion rather than enhanced chromospheric activity. The spatial position of some of these stars is coincident with that of known X-ray sources, and their location in the colour–magnitude diagram allows us to classify them as active binaries or cataclysmic variables (CVs). We show that this method, commonly adopted to study accreting discs in young stellar objects, can be successfully used to identify and characterize candidate CVs.

THE NATURE OF STARBURSTS. III. THE SPATIAL DISTRIBUTION OF STAR FORMATION

We map the spatial distribution of recent star formation over a few x 100 Myr timescales in fifteen starburst dwarf galaxies using the location of young blue helium burning stars identified from optically resolved stellar populations in archival Hubble Space Telescope observations. By comparing the star formation histories from both the high surface brightness central regions and the diffuse outer regions, we measure the degree to which the star formation has been centrally concentrated during the galaxies' starbursts, using three different metrics for the spatial concentration. We find that the galaxies span a full range in spatial concentration, from highly centralized to broadly distributed star formation. Since most starbursts have historically been identified by relatively short timescale star formation tracers (e.g., Halpha emission), there could be a strong bias towards classifying only those galaxies with recent, centralized star formation as starbursts, while missing starbursts that are spatially distributed.

THE PANCHROMATIC HUBBLE ANDROMEDA TREASURY. II. TRACING THE INNER M31 HALO WITH BLUE HORIZONTAL BRANCH STARS

We attempt to constrain the shape of M31's inner stellar halo by tracing the surface density of blue horizontal branch (BHB) stars at galactocentric distances ranging from 2 kpc to 35 kpc. Our measurements make use of resolved stellar photometry from a section of the Panchromatic Hubble Andromeda Treasury (PHAT) survey, supplemented by several archival Hubble Space Telescope observations. We find that the ratio of BHB to red giant stars is relatively constant outside of 10 kpc, suggesting that the BHB is as reliable a tracer of the halo population as the red giant branch. In the inner halo, we do not expect BHB stars to be produced by the high metallicity bulge and disk, making BHB stars a good candidate to be a reliable tracer of the stellar halo to much smaller galactocentric distances. If we assume a power-law profile r^(-\alpha) for the 2-D projected surface density BHB distribution, we obtain a high-quality fit with a 2-D power-law index of \alpha=2.6^{+0.3}_{-0.2} outside of 3 kpc, which flattens to \alpha<1.2 inside of 3 kpc. This slope is consistent with previous measurements but is anchored to a radial baseline that extends much farther inward. Finally, assuming azimuthal symmetry and a constant mass-to-light ratio, the best-fitting profile yields a total halo stellar mass of 2.1^{+1.7}_{-0.4} x 10^9 M_sun. These properties are comparable with both simulations of stellar halo formation formed by satellite disruption alone, and with simulations that include some in situ formation of halo stars.

OBSERVATIONAL CONSTRAINTS ON RED AND BLUE HELIUM BURNING SEQUENCES

We derive the optical luminosity, colors, and ratios of the blue and red helium burning (HeB) stellar populations from archival Hubble Space Telescope observations of nineteen starburst dwarf galaxies and compare them with theoretical isochrones from Padova stellar evolution models across metallicities from Z=0.001 to 0.009. We find that the observational data and the theoretical isochrones for both blue and red HeB populations overlap in optical luminosities and colors and the observed and predicted blue to red HeB ratios agree for stars older than 50 Myr over the time bins studied. These findings confirm the usefulness of applying isochrones to interpret observations of HeB populations. However, there are significant differences, especially for the red HeB population. Specifically we find: (1) offsets in color between the observations and theoretical isochrones of order 0.15 mag (0.5 mag) for the blue (red) HeB populations brighter than M_V ~ -4 mag, which cannot be solely due to differential extinction; (2) blue HeB stars fainter than M_V ~ -3 mag are bluer than predicted; (3) the slope of the red HeB sequence is shallower than predicted by a factor of ~3; and (4) the models overpredict the ratio of the most luminous blue to red HeB stars corresponding to ages <50 Myr. Additionally, we find that for the more metal-rich galaxies in our sample (Z> 0.5 Zsolar) the red HeB stars overlap with the red giant branch stars in the color magnitude diagrams, thus reducing their usefulness as indicators of star formation for ages >100 Myr.

AHUBBLE SPACE TELESCOPESTUDY OF LYMAN LIMIT SYSTEMS: CENSUS AND EVOLUTION

We present a survey for optically thick Lyman limit absorbers at z<2.6 using archival Hubble Space Telescope observations with the Faint Object Spectrograph and Space Telescope Imaging Spectrograph. We identify 206 Lyman limit systems (LLSs) increasing the number of catalogued LLSs at z<2.6 by a factor of ~10. We compile a statistical sample of 50 tau_LLS > 2 LLSs drawn from 249 QSO sight lines that avoid known targeting biases. The incidence of such LLSs per unit redshift, l(z)=dn/dz, at these redshifts is well described by a single power law, l(z) = C1 (1+z)^gamma, with gamma=1.33 +/- 0.61 at z<2.6, or with gamma=1.83 +/- 0.21 over the redshift range 0.2 < z < 4.9. The incidence of LLSs per absorption distance, l(X), decreases by a factor of ~1.5 over the ~0.6 Gyr from z=4.9 to 3.5; l(X) evolves much more slowly at low redshifts, decreasing by a similar factor over the ~8 Gyr from z=2.6 to 0.25. We show that the column density distribution function, f(N(HI)), at low redshift is not well fitted by a single power law index (f(N(HI)) = C2 N(HI)^(-beta)) over the column density range 13 < log N(HI) < 22 or log N(HI) >17.2. While low and high redshift f(N(HI)) distributions are consistent for log N(HI)>19.0, there is some evidence that f(N(HI)) evolves with z for log N(HI) < 17.7, possibly due to the evolution of the UV background and galactic feedback. Assuming LLSs are associated with individual galaxies, we show that the physical cross section of the optically thick envelopes of galaxies decreased by a factor of ~9 from z~5 to 2 and has remained relatively constant since that time. We argue that a significant fraction of the observed population of LLSs arises in the circumgalactic gas of sub-L* galaxies.

Ionized gas and stars in the central kiloparsec of the type 2 Seyfert galaxy NGC 2110 - I. The data

In this paper, we present new results from an extensive set of Hubble Space Telescope (HST) and ground-based observations of the Seyfert galaxy NGC 2110. The HST data sets include Wide Field Planetary Camera 2 (WFPC2) observations as well as the first high-spatial resolution spectroscopy ([O i], [N ii], Hα and [S ii] lines) of this galaxy obtained using the Space Telescope Imaging Spectrograph (STIS). The ground-based data are three-dimensional (x, y, λ) spectrographic observations obtained using the integral field spectrograph OASIS on the Canada–France–Hawaii Telescope, complemented by near-infrared long-slit spectroscopy obtained using NIRSPEC on the Keck–II telescope. The OASIS observations cover regions containing both stellar absorption lines and major optical emission lines. The NIRSPEC observations cover the H and K bands. Combined with archival HST observations, the WFPC2 data provide us with a high-spatial resolution extinction map. The OASIS data allowed bidimensional mapping of the stellar and gaseous kinematics, as well as of the spectral properties of the ionized gas. These results are compared to those obtained in the near-infrared with NIRSPEC/Keck. Last, we used the STIS data to probe the ionized gas kinematics and properties in the inner 4 arcsec along PA = 156° at unprecedented spatial resolution. Our two-dimensional (2D) map of the stellar velocity field and the near-infrared stellar velocity profile are symmetric about the nucleus, confirming the results of previous long-slit observations. The asymmetry of the velocity field of the ionized gas is present at the same level for visible and near-infrared lines, indicating this is not a reddening effect. Multiple Gaussian fitting of the emission-line profile allowed the contributions of the broad and a narrow components to be disentangled. The intensity peak of the [O iii] narrow component is located north of the nucleus, indicating that the bulk of the narrow [O iii] emission comes from the jet-like structure (Mulchaey et al.) and not from the nucleus itself. We suggest that the northern arm is the anomalous one, contrary to what has been claimed earlier. Last, we also show that the elongated region of high gas velocity dispersion located close to the nucleus and discovered by Gonzalez Delgado et al. is intrinsic to the narrow component.

Star cluster formation and evolution in nearby starburst galaxies -- I. Systematic uncertainties

The large majority of extragalactic star cluster studies done to date have essentially used two or three-passband aperture photometry, combined with theoretical stellar population synthesis models, to obtain age, mass and extinction estimates, and sometimes also metallicities. The accuracy to which this can be done depends on the choice of (broad-band) passband combination and, crucially, also on the actual wavelengths and the wavelength range covered by the observations. Understanding the inherent systematic uncertainties (the main aim of this paper) is of the utmost importance for a well-balanced interpretation of the properties of extragalactic star cluster systems. We simultaneously obtain ages, metallicities and extinction values for ∼ 300 clusters in the nearby starburst galaxy NGC 3310, based on archival Hubble Space Telescope observations from the ultraviolet (UV) to the near-infrared (NIR). We show that, for ages 6 . log(age/yr) . 9, and if one can only obtain partial coverage of the spectral energy distribution (SED), an optical passband combination of at least four filters including both blue and red passbands results in the most representative age distribution, as compared to the better constrained ages obtained from the full UV–NIR SED coverage. We find that while blue-selected passband combinations lead to age distributions that are slightly biased towards younger ages due to the well-known age–metallicity degeneracy, red-dominated passband combinations should be avoided. NGC 3310 underwent a (possibly extended) global burst of cluster formation ∼ 3×10 7 yr ago. This coincides closely with the last tidal interaction or merger with a low-metallicity galaxy that likely induced the formation of the large fraction of clusters with (significantly) subsolar metallicities. The logarithmic slope of the Vband cluster luminosity function, for clusters in the range 17.7 . F606W . 20.2

Far-Ultraviolet Emission from Elliptical Galaxies at [CLC][ITAL]z[/ITAL][/CLC] = 0.33

We present far-ultraviolet (far-UV) images of the rich galaxy cluster ZwCl1358.1+6245, taken with the Space Telescope Imaging Spectrograph on board the Hubble Space Telescope (HST). When combined with archival HST observations, our data provide a measurement of the UV-to-optical flux ratio in 8 early-type galaxies at z=0.33. Because the UV flux originates in a population of evolved, hot, horizontal branch (HB) stars, this ratio is potentially one of the most sensitive tracers of age in old populations -- it is expected to fade rapidly with lookback time. We find that the UV emission in these galaxies, at a lookback time of 3.9 Gyr, is significantly weaker than it is in the current epoch, yet similar to that in galaxies at a lookback time of 5.6 Gyr. Taken at face value, these measurements imply different formation epochs for the massive ellipticals in these clusters, but an alternative explanation is a floor in the UV emission due to a dispersion in the parameters that govern HB morphology.

Surface brightness profiles and structural parameters for globular clusters in the Fornax and Sagittarius dwarf spheroidal galaxies

We present radial surface brightness profiles for all five globular clusters in the Fornax dwarf spheroidal galaxy, and for the four present members of the Sagittarius dwarf spheroidal galaxy. These profiles are derived from archival Hubble Space Telescope observations, and have been calculated using the same techniques with which we measured profiles in our previous studies of Large and Small Magellanic Cloud (LMC and SMC) clusters, apart from some small modifications. From the surface brightness profiles, we have determined structural parameters for each cluster, including core radii and luminosity and mass estimates. We also provide a brief summary of literature measurements of other parameters for these clusters, including their ages, metallicities and distances. Our core radius measurements are mostly in good agreement with those from previous lower resolution studies, although for several clusters our new values are significantly different. The profile for Fornax cluster 5 does not appear to be well fitted by a King-type model and we suggest that it is a post-core-collapse candidate. We examine the distribution of cluster core radii in each of the two dwarf galaxy systems and compare these with the distribution of core radii for old LMC clusters. The three distributions match within the limits of measurement errors and the small-sample sizes. We discuss the implications of this in the context of the radius—age trend we have previously highlighted for the Magellanic Cloud clusters.

A Search for the Optical Counterpart of the Luminous X-Ray Source in NGC 6652

We examine images of the field of X1832-330, the luminous (Lx ~ 10^36 erg/s) X-ray burst source near the center of the globular cluster NGC 6652, in order to identify the optical counterpart for further study. U and B ground-based images allow us to set a limit M_B > 3.5 for the counterpart at the time of those observations, provided that the color is (U-B)_0 ~ -1, similar to the sources known in other clusters. Archival Hubble Space Telescope observations survey most but not all of the 1 sigma X-ray error circle, and allow us to set limits M_B > 5.9 and M_B > 5.2 in the WF/PC and WFPC2 regions, respectively. In the WF/PC images we do weakly detect a faint object with UV-excess, but it is located 11.7'' from the ROSAT X-ray position. This considerable (2.3 sigma) discrepancy in position suggests that this candidate be treated with caution, but it remains the only reasonable one advanced thus far. We measure for this star m_439 = 20.2 +- 0.2, (m_336 - m_439) = -0.5 +- 0.2, and estimate M_B = 5.5, (U-B)_0 = -0.9, similar to other known optical counterparts. If this candidate is not the identification, our limits imply that the true counterpart, not yet identified, is probably the optically-faintest cluster source yet known, or alternatively that it did not show significant UV excess at the time of these observations. Finally, we assess the outlook for the identification of the remaining luminous globular cluster X-ray sources.