The CIAO (Chandra Interactive Analysis of Observations) software package was first released in 1999 following the launch of the Chandra X-ray Observatory and is used by astronomers across the world to analyze Chandra data as well as data from other telescopes. From the earliest design discussions, CIAO was planned as a generalpurpose scientific data analysis system optimized for X-ray astronomy, and consists mainly of command line tools (allowing easy pipelining and scripting) with a parameter-based interface layered on a flexible data manipulation I/O library. The same code is used for the standard Chandra archive pipeline, allowing users to recalibrate their data in a consistent way.We will discuss the lessons learned from the first six years of the software's evolution. Our initial approach to documentation evolved to concentrate on recipe-based "threads" which have proved very successful. A multidimensional abstract approach to data analysis has allowed new capabilities to be added while retaining existing interfaces. A key requirement for our community was interoperability with other data analysis systems, leading us to adopt standard file formats and an architecture which was as robust as possible to the input of foreign data files, as well as re-using a number of external libraries. We support users who are comfortable with coding themselves via a flexible user scripting paradigm, while the availability of tightly constrained pipeline programs are of benefit to less computationally-advanced users. As with other analysis systems, we have found that infrastructure maintenance and re-engineering is a necessary and significant ongoing effort and needs to be planned in to any long-lived astronomy software.

We used the Spitzer Space Telescope's Infrared Spectrograph to map nearly the entire extent of Cassiopeia A between 5-40 µm. Using infrared and Chandra X-ray Doppler velocity measurements, along with the locations of optical ejecta beyond the forward shock, we constructed a 3-D model of the remnant. The structure of Cas A can be characterized into a spherical component, a tilted thick disk, and multiple ejecta jets/pistons and optical fast-moving knots all populating the thick disk plane. The Bright Ring in Cas A identifies the intersection between the thick plane/pistons and a roughly spherical reverse shock. The ejecta pistons indicate a radial velocity gradient in the explosion. Some ejecta pistons are bipolar with oppositely-directed flows about the expansion center while some ejecta pistons show no such symmetry. Some ejecta pistons appear to maintain the integrity of the nuclear burning layers while others appear to have punched through the outer layers. The ejecta pistons indicate a radial velocity gradient in the explosion. In 3-D, the Fe jet in the southeast occupies a "hole" in the Si-group emission and does not represent "overturning", as previously thought. Although interaction with the circumstellar medium affects the detailed appearance of the remnant and may affect the visibility of the southeast Fe jet, the bulk of the symmetries and asymmetries in Cas A are intrinsic to the explosion.

We present the 2-60 keV spectrum of the supernova remnant Cassiopeia A measured using the Proportional Counter Array and the High Energy X-ray Timing Experiment on the Rossi X-ray Timing Explorer satellite. In addition to the previously reported strong emission-line features produced by thermal plasmas, the broad-band spectrum has a high-energy "tail" that extends to energies at least as high as 120 keV. This tail may be described by a broken power law that has photon indices of Γ 1 = 1.8 +0.5 −0.6 and Γ 2 = 3.04 +0.15 −0.13 and a break energy of E b = 15.9 +0.3 −0.4 keV. We argue that the high-energy component, which dominates the spectrum above about 10 keV, is produced by synchrotron radiation from electrons that have energies up to at least 40 TeV. This conclusion supports the hypothesis that Galactic cosmic rays are accelerated predominantly in supernova remnants.

We present the results of a joint spectral analysis of RXTE PCA, ASCA SIS, and ROSAT PSPC data of the supernova remnant SN 1006. This work represents the first attempt to model both the thermal and nonthermal X-ray emission over the entire X-ray energy band. The thermal flux is described by a nonequilibrium ionization model with an electron temperature kT e = 0.6 keV, an ionization timescale n 0 t = 9 × 10 9 cm −3 s, and a relative elemental abundance of silicon that is 10-18 times larger than the solar abundance. The nonthermal X-ray spectrum is described by a broken power law model with low-and high-energy photon indices Γ 1 = 2.1 and Γ 2 = 3.0, respectively. Since the nonthermal X-ray spectrum steepens with increasing energy, the results of the present analysis corroborate previous claims that the nonthermal X-ray emission is produced by synchrotron radiation. We argue that the magnetic field strength is significantly larger than previous estimates of about 10 µG and arbitrarily use a value of 40 µG to estimate the parameters of the cosmic-ray electron, proton, and helium spectra of the remnant. The results for the ratio of the number densities of protons and electrons (R = 160 at 1 GeV), the total energy in cosmic rays (E cr = 1 × 10 50 ergs), and the spectral index of the electrons at 1 GeV (Γ e = 2.14 ± 0.12) are consistent with the hypothesis that Galactic cosmic rays are accelerated predominantly in the shocks of supernova remnants. Yet, the remnant may or may not accelerate nuclei to energies as high as the energy of the "knee," depending on the reason why the maximum energy of the electrons is only 10 TeV.