Alan M. Watson

THE PHOTOMETRIC PERFORMANCE AND CALIBRATION OF WFPC2

We discuss the photometric performance and calibration of the Wide Field Planetary Camera 2 (WFPC2) on the Hubble Space Telesopce (HST). The stability and accuracy of WFPC2 photometric measurements is discussed, with particular attention given to charge transfer efficiency (CTE) effects, contamination effects in the ultraviolet (UV), and flat field accuracy and normalization. Observational data are presented from both WFPC2 observations and ground observations using a system similar to that flown. WFPC2 photometric systems are defined both for the ground and flight systems. Transformations between these systems and the Landolt UBVRI system are presented. These transformations are sensitive to details in the spectra being transformed, and these sensitivities are quantified and discussed. On-orbit observations are used to revise the prelaunch estimates of response curves to best match synthetic photometry results with observations, and the accuracy of the resulting synthetic photometry is discussed. Synthetic photometry is used to determine zeropoints and transformations for all of the fight filters, and also to derive interstellar extinction values for the WFPC2 system. Using stellar interior and atmosphere models, isochrones in the WFPC2 system are calculated and compared with several observations.

THE PERFORMANCE AND CALIBRATION OF WFPC2 ON THE HUBBLE SPACE TELESCOPE

The WFPC2 was installed in the Hubble Space Telescope (HST) in 1993 December. Since then, the instrument has been providing high-quality images. A significant amount of calibration data has been collected to aid in the understanding of the on-orbit performance of the instrument. Generally, the behavior of the camera is similar to its performance during the system-level thermal vacuum test at JPL in 1993 May. Surprises were a significant charge-transfer-efficiency (CTE) problem and a significant growth rate in hot pixels at the original operating temperature of the CCDs (-76 °C). The operating temperature of the WFPC2 CCDs was changed to -88 °C on 1994 April 23, and significant improvements in CTE and hot pixels are seen at this temperature. In this paper we describe the on-orbit performance of the WFPC2. We discuss the optical and thermal history, the instrument throughput and stability, the PSF, the effects of undersampling on photometry, the properties of cosmic rays observed on-orbit, and the geometric distortion in the camera. We present the best techniques for the reduction of WFPC2 data, and describe the construction of calibration products including superbiases, superdarks, and fiat fields.

Far-Ultraviolet Imaging of Jupiter's Aurora and the Io “Footprint”

Far-ultraviolet images of Jupiter from the Hubble Space Telescope Wide Field Planetary Camera 2 reveal polar auroral emissions at 300 kilometer resolution and three times higher sensitivity than previously achieved. Persistent features include a main oval containing most of the emission and magnetically connected to the middle magnetosphere, diffuse and variable emissions poleward of the main oval, and discrete emission from Ios magnetic footprint equatorward of the oval. The auroral emissions are variable, exhibit magnetic conjugacy, and are visible above the planet limb. All emissions approximately co-rotate with Jupiter except the Io “footprint,” which is fixed along Ios magnetic flux tube.

Benchmark problems for continuum radiative transfer. High optical depths, anisotropic scattering, and polarisation

Solving the continuum radiative transfer equation in high opacity media requires sophisticated numerical tools. In order to test the reliability of such tools, we present a benchmark of radiative transfer codes in a 2D disc configuration. We test the accuracy of seven independently developed radiative transfer codes by comparing the temperature structures, spectral energy distributions, scattered light images, and linear polarisation maps that each model predicts for a variety of disc opacities and viewing angles. The test cases have been chosen to be numerically challenging, with midplane optical depths up 10^6, a sharp density transition at the inner edge and complex scattering matrices. We also review recent progress in the implementation of the Monte Carlo method that allow an efficient solution to these kinds of problems and discuss the advantages and limitations of Monte Carlo codes compared to those of discrete ordinate codes. For each of the test cases, the predicted results from the radiative transfer codes are within good agreement. The results indicate that these codes can be confidently used to interpret present and future observations of protoplanetary discs.

Hubble Space Telescope Observations of the SN 1987A Triple Ring Nebula

We have observed SN 1987A with the optically corrected WFPC2 on the Hubble Space Telescope both in emission lines and in the UV and optical continuum. The previously observed outer nebular structure is shown to be part of two closed unresolved loops. These loops were flash-ionized by the supernova itself. They are not caused by limb brightening of an hourglass shell produced by the interaction of the winds from the progenitor. The inner ring is seen to be extended and may be connected to the new outer rings by sheets of material. However, beyond the outer rings, emission is not seen, implying a very low density (n 1000. This density contrast of at least 100 is difficult to reconcile with the conventional picture of the progenitor evolution. Two models for the rings are presented, but each is deficient in important respects. A proper understanding of this system will require new physical insight.

The Luminosity Function and Initial Mass Function in the Galactic Bulge

We present deep photometry obtained with the Hubble Space Telescope in a field in Baades window in the Galactic bulge. We derive a luminosity function down to I ~ 24.3, or V ~ 27.5, corresponding to M ~ 0.3 M_☉. The luminosity function from the turnoff down to this level appears remarkably similar to that observed in the solar neighborhood. We derive a mass function using both an empirical local mass-luminosity relation and a mass-luminosity relation from recent stellar model calculations, allowing for the presence of binaries and photometric errors. The mass function has a power-law form with dN/dM ∝ M^(-2.2) for M ≳ 0.7 M_☉. However, we find strong evidence for a break in the mass function slope around 0.5–0.7 M_☉, with a significantly shallower slope at lower masses. The value of the slope for the low masses depends on the assumed binary fraction and the accuracy of our completeness correction. This mass function should directly reflect the initial mass function.

The X-ray counterpart to the gravitational-wave event GW170817

A long-standing paradigm in astrophysics is that collisions—or mergers—of two neutron stars form highly relativistic and collimated outflows (jets) that power γ-ray bursts of short (less than two seconds) duration. The observational support for this model, however, is only indirect. A hitherto outstanding prediction is that gravitational-wave events from such mergers should be associated with γ-ray bursts, and that a majority of these bursts should be seen off-axis, that is, they should point away from Earth. Here we report the discovery observations of the X-ray counterpart associated with the gravitational-wave event GW170817. Although the electromagnetic counterpart at optical and infrared frequencies is dominated by the radioactive glow (known as a ‘kilonova’) from freshly synthesized rapid neutron capture (r-process) material in the merger ejecta, observations at X-ray and, later, radio frequencies are consistent with a short γ-ray burst viewed off-axis. Our detection of X-ray emission at a location coincident with the kilonova transient provides the missing observational link between short γ-ray bursts and gravitational waves from neutron-star mergers, and gives independent confirmation of the collimated nature of the γ-ray-burst emission.

Detection of the Tip of Red Giant Branc in NGC 5128

We present a color-magnitude diagram of more than 10,000 stars in the halo of the galaxy NGC 5128 (Centaurus A), based on WFPC2 images through V-band and I-band filters. The position of the red giant branch (RGB) stars is compared with the loci of the red giant branch in six well-studied globular clusters and in the dwarf elliptical galaxy NGC 185; the tip of the RGB is signalled by an observed turn-up in the luminosity function at I≃24.1 ± 0.1 mag; this yields a distance modulus (m - M)_0 = 27.8 ± 0.2 for NGC 5128 (i.e., a distance of 3.6 ± 0.2 Mpc), in agreement with previous determinations based on the planetary nebulae luminosity function and on the surface brightness fluctuations technique. The presence of an intermediate-age stellar population (~5 Gyr) is suggested by the luminosity function of the asymptotic giant branch stars, extending up to I= 22.6 mag (for V - I > 2) and M_(bot) ~ -5 mag; however, the number of these stars constrains the intermediate-age stellar population in the halo of NGC 5128 to be less than ~ 10% of the total. The color distribution at constant I magnitude, albeit affected by the completeness level of our sample, strongly suggests a mean value of [Fe/H] > -0.9 dex, possibly similar to the value found in M31 and higher than that observed in NGC 185. Like the M31 halo, the halo of NGC 5128 exhibits a broad range of levels of chemical enrichment.

Stellar populations in three outer fields of the Large Magellanic Cloud

We present Hubble Space Telescope photometry for three fields in the outer disk of the Large Magellanic Cloud (LMC) extending approximately 4 mag below the faintest main-sequence turnoff. We cannot detect any strongly significant differences in the stellar populations of the three fields based on the morphologies of the color-magnitude diagrams, the luminosity functions, and the relative numbers of stars in different evolutionary stages. Our observations therefore suggest similar star formation histories in these regions, although some variations are certainly allowed. The fields are located in two regions of the LMC: one is in the northeast and two are located in the northwest. Under the assumption of a common star formation history, we combine the three fields with ground-based data at the same location as one of the fields to improve statistics for the brightest stars. We compare this stellar population with those predicted from several simple star formation histories suggested in the literature, using a combination of the R-method of Bertelli et al. (1992) and comparisons with the observed luminosity function. The only model we consider that is not rejected by the observations is one in which the star formation rate is roughly constant for most of the LMCs history and then increases by a factor of 3 about 2 Gyr ago. Such a model has roughly equal numbers of stars older and younger than 4 Gyr, and thus is not dominated by young stars. This star formation history, combined with a closed-box chemical evolution model, is consistent with observations that the metallicity of the LMC has doubled in the past 2 Gyr.

Stellar populations in the Large Magellanic Cloud: Evidence for a significant number of older stars or a steeper IMF?

We present deep photometry obtained with the HST in an outer LMC field. A well-defined main sequence is seen down to V > 26. We derive a luminosity function from the data and use it to constrain the IMF and the star formation history. We derive limits on the IMF slope, ɑ(with dN/∝dM^ɑ), from stars on the main sequence which are fainter than the oldest turnoff. For most choices of star formation history and metallicity, we derive slopes which are consistent the Salpeter (ɑ= -2.35) or local solar neighborhood IMF, although the preferred values are steeper. We can rule out IMF slopes shallower than -1.6 and steeper than -3.1 for the mass range 0.6≾M≾1.1 M_☉. Assuming a Salpeter IMF over the entire observed mass range, we derive star formation histories from the entire luminosity function, which covers the mass range 0.6≾M≾1.1 3 M_☉. We find that the luminosity function is inconsistent with the scenario in which the bulk of the field stars in the LMC are younger than 4 Gyr. Instead, we find that there must be a comparable number of stars older and younger than 4 Gyr. Our best model has a star formation rate which is roughly constant for 10 Gyr then increases by about a factor of three for the past 2 Gyr. Such a model is also roughly consistent with the distribution of stars in the color-magnitude diagram. Similar model parameters are derived if we adopt the Kroupa, Tout, and Gilmore solar neighborhood LMF instead of a Salpeter slope. Alternatively, we can fit the luminosity function with a predominantly young population if we use a steeper single power law IMF slope with ɑ ~-2.75 over the entire range of observed masses.