# Presentation: 182756

## High Resolution X-ray Spectra

- Randall K. Smith
- Chandra X-ray Center

- The Glory and the Grandeur

## Introduction

- It has been said that Chandra and XMM/Newton have finally brought high-resolution spectra to X-ray astronomy, but this is not strictly true.

- Chandra & XMM/Newton have made high-resolution spectra common, but the first high-resolution spectra to be published was:

- A 7 ksec exposure of Capella with the *Einstein* Solid State Spectrometer (SSS) (Holt *et al*. 1979)

## Introduction

- Of course, it could be argued that the Einstein SSS is not really high-resolution; it is only ~ CCD resolution. The first published grating spectrum was:

- A 42 ksec observation of **Capella** with the Einstein Objective Grating Spectrometer (OGS) by Mewe *et al*. (1982)

## Introduction

- Close on their heels with an even higher resolution spectrometer was:

- Vedder & Canizares (1983) with a 59 ksec observation of (yes!) **Capella** with the *Einstein* Focal Plane Crystal Spectrometer (FPCS) :

## Introduction

- After the *Einstein* mission, EXOSAT was launched.  Although it primary claim to fame is timing, it did manage to get a high resolution spectrum of:

- Capella, using the Transmission Grating Spectrometer (TGS) in a 85 ksec observation (Lemen *et al*. 1989)

## Introduction

- Coming into the more modern era (i.e., when I started in this business), there was a joint analysis of EUV/X-ray emission from:

- Capella, using a 21 ksec with the ASCA Silicon Imaging Spectrometer  (SIS) and a 120 ksec with the EUVE spectrometer  (Brickhouse *et al*. 2000)

- Capella

## Introduction

- Once *Chandra* was launched, two new gratings became available.  Here is spectrum of Capella observed for 95 ksec with the Low Energy Transmission Grating (LETG) using the High-Resolution Camera (HRC) as a detector. (Brinkman *et al*. 2000)

## Introduction

- And, of course, the High-Energy Transmission Grating (HETG) also observed Capella, with the Advanced CCD Imaging Spectrometer (ACIS) for 89 ksec (Canizares *et al.* 2000)

- Capella

## Introduction

- With the launch of XMM/Newton, the Reflection Grating Spectrometer (RGS) was finally able to achieve the one true goal of all spectroscopists: a 52 ksec observation of

- Audard *et al*. 2001

- CAPELLA!

- Capella.

## Grating Basics

*Chandra*; notice the position of the removable gratings:

- So, where do these gratings go? Here is a diagram of

## Grating Basics

- Here is a the HETG, pre-installation.  The many grating facets are all carefully positioned to catch the light from the 4 nested mirrors.

## Grating Basics

- The LETG, pre-installation.  Note the support structures, which will later show up as scattered light in deep LETG observations:

## Grating Basics

- Canizares et al. (2000)

- The two arms of the HEG and MEG gratings can be clearly seen here.  Notice the larger effective area of the MEG; the exploded region shows the effect of the HEG’s higher resolution.

## Grating Basics

- Canizares et al. (2000)

- X-ray gratings work identically to optical gratings; both follow the grating equation:

- sin  = m/p

- where  is the wavelength,  is the dispersion angle (measured from the zero-order image),  p is the spatial period of the grating itself, and m is the order (1st, 2nd, etc).

## Grating Basics

- High-resolution (grating) spectra on Chandra cover a huge range of wavelengths: from 1.2-170Å, over two orders of magnitude.  Note that although “wavelength” is the natural unit for grating data, many (older) X-ray astronomers use energy (keV) units anyway.

- LETG/HRC-S observation of NGC6624

- If ACIS is used, the CCD resolution can distinguish between the different orders; on the HRC, this must be modeled.

- The spatial and spectral elements are tightly coupled.  If the zero-order image is slightly displaced (*i.e.*, if the source is piled-up), the +/- order wavelengths will be offset from each other.  If so, it calls for reprocessing.

- sin  = m/p

## Grating Basics

- The light blue lines show the extraction region in the cross-dispersion direction.  Notice the enhanced background on the rightmost chips; this can be (largely) eliminated by using the energy resolution of the CCDs.  Also, one chip has failed; this also happened on RGS2, but fortunately at a different position.

- The star AB Dor observed with the XMM RGS1

## Grating Basics

- So, when you download grating data, what do you actually have?

| ACIS/HETG | 12 spectra; ±1,2,3 orders for both HEG and MEG |
| --- | --- |
| ACIS/LETG | 6 spectra; ±1,2,3 orders for the LEG |
| HRC/LETG | 2 spectra; ±∑ni orders |
| HRC/HETG | Not recommended |
| XMM/RGS | 4 spectra; +1,2 from each of two RGSs. |

- For Chandra data,  these are stored in one file with multiple spectra, a PHA Type II file; for XMM data, each spectrum is in a different PHA Type I file.

## Grating Basics

| Grating | 
(Å) | Resolution @ 1keV | Eff. Area @ 1keV (cm-2) |
| --- | --- | --- | --- |
| HEG 1st | 0.012 | 1000 | 10 |
| MEG 1st | 0.023 | 540 | 34 |
| LEG 1st | 0.05 | 250 | 25 |
| RGS 1st | ~0.06 | 200 | 150 |

- The characteristics of the different gratings in their prime arrangement:

## Grating Basics

- Consider a line at 1.0 keV (12.398Å), on top of a 0.5 keV blackbody

- So, what can we do with grating data that we can’t with CCDs?

## Grating Basics

- Now consider the same spectrum observed by each of the four gratings:

- If this data were fit assuming a pure Gaussian response, the scattering wings of the RGS would mean we’d measure only 1/2 the true power!

## Processing

- This talk really isn’t about actual data processing; if you find yourself needing to reduce *Chandra* and/or *XMM/Newton* grating data,  I recommend:

| RGS | http://heasarc.gsfc.nasa.gov/docs/xmm/abc/node8.html |
| --- | --- |
| HETG/ACIS | http://cxc.harvard.edu/ciao/threads/spectra_hetgacis.html |
| LETG/ACIS | http://cxc.harvard.edu/ciao/threads/spectra_letgacis.html |
| LETG/HRC | http://cxc.harvard.edu/ciao/threads/spectra_letghrc.html |

- The end product will be either multiple Type I or a Type II PHA file(s); basically, a histogram of counts as a function of wavelength, along with the instrumental response.

## Processing

- Identify the dispersed source(s).  In some cases, multiple sources may be present; both CIAO and SAS software can identify and reduce more than one source, but the user must select which source to extract.
- Calculate the dispersion distance for each event, and (if available) resolve the order by comparing the dispersion wavelength to the measured CCD energy.
- Extract the spectrum into a PHA file, suitable for use in Sherpa or XSPEC.
- Create the detector response for the source.  This depends on the position of the source on the detector, the aspect solution, and the operating mode of the telescope.  The variations are not generally large, so using another observation's response is OK for quick work.

- The XMM/Newton RGS is a reflection grating, while the Chandra LETG and HETG are transmission gratings.  However, the basic processing steps are similar:

- SAS does steps 1-3 with one command, **rgsproc**.  *Chandra* breaks
- steps 1-3 into subprograms: **tg_create_mask**, **tg_resolve_events**,
- and **tgextract**.

## Processing

- Know your Data: RGS event files

- unix% dmlist rgs_evt2.fits cols

## Processing

- Know your Data: HETG/ACIS event files

- unix% dmlist acis_evt2.fits cols

## Processing

- unix% dmlist hrc_evt2.fits cols

- Know your Data: LETG/HRC event files

## Processing

- For bright sources on ACIS-S, the background is likely negligible.  However, on the HRC-S or with the RGS it isn’t:

**Source/Background regions for the HRC/LETG**

- Proper background subtraction is still a topic of some debate!

## Processing

- Using a detector with even moderate energy resolution can substantially reduce the background, because one can cross-correlate the wavelength measured from the dispersed position with the energy measured in the CCD using the relation:

- XMM RGS1 spectra of the star AB Dor

- EkeV = 12.398 / Å

## Processing

- What if you want both high-resolution spectra AND timing?

- 1) You’d better have a bright source

- 2) With *Chandra*, the ACIS detectors can be run in continuous clocking mode while either of the gratings are in place.  The time resolution is about 3 ms.

- 3) With XMM/Newton, the RGS detectors can be run in CC mode, but as of this writing (May 2003) this mode is not available to users for technical reasons.

- 4) The HRC has good time resolution, but due to an electronics problem, the time tags for each event suffer from an off-by-one problem and actually refer to the previous event.  It is possible to run the HRC in a “timing” mode where this is mitigated, but not with the full grating array.

## Processing

- Unfortunately, not every X-ray source is as bright as Capella or NGC6624.  What can you do to enhance the S/N?
- In this case, you have 4 choices:

**Co-add plus/minus orders of the same grating**  *Can broaden lines if zero-order is offset.*
**Co-adding different gratings, such as HEG and MEG or RGS1 and RGS2**  *Complicates line shape function.*
**Co-add separate observations**  *Instrumental background can vary, plus the same issues of zero-order offsets.*
**Co-add separate observations and instruments**   *All of the above*

## Grating Analysis

- unix% sherpa
- sherpa> data acis_pha2.fits
- sherpa> paramprompt off
- sherpa> rsp[hm1]
- sherpa> rsp[hp1]
- sherpa> rsp[mm1]
- sherpa> rsp[mp1]
- sherpa> hm1.rmf = acisheg1D1999-07-22rmfN0004.fits
- sherpa> hm1.arf = acisf01318HEG_-1_garf.fits
- sherpa> hp1.rmf = acisheg1D1999-07-22rmfN0004.fits
- sherpa> hp1.arf = acisf01318HEG_1_garf.fits
- sherpa> mm1.rmf = acismeg1D1999-07-22rmfN0004.fits
- sherpa> mm1.arf = acisf01318MEG_-1_garf.fits
- sherpa> mp1.rmf = acismeg1D1999-07-22rmfN0004.fits
- sherpa> mp1.arf = acisf01318MEG_1_garf.fits
- sherpa> instrument 3 = hm1
- sherpa> instrument 4 = hp1
- sherpa> instrument 9 = mm1
- sherpa> instrument 10= mp1
- sherpa> ignore allsets all
- sherpa> notice allsets wave 14.9:15.4
- sherpa> source 3,4,9,10 = poly[b1] + delta1d[l1] + delta1d[l2] + delta1d[l3]
- sherpa> l1.pos = 15.014
- sherpa> l2.pos = 15.079
- sherpa> l3.pos = 15.2610
- sherpa> freeze l1.pos
- sherpa> freeze l2.pos
- sherpa> freeze l3.pos
- sherpa> fit
- sherpa> lp 4 fit 3 fit 4 fit 9 fit 10
- sherpa> import(``guide'')
- sherpa> mdl2latex
- \begin{tabular}{lllllll}
- ModelName & Line Model & Position & Flux & Flux Error & Fit Data & Label \\
-  &  & Angstrom & ph/cm$^2$/s & ph/cm$^2$/s &  \\
- l1 & delta1d & 15.014 & 0.00308923 & 6.7101e-05 & 3,4,9,10 &  \\
- l2 & delta1d & 15.079 & 0.000270431 & 2.81612e-05 & 3,4,9,10 &  \\
- l3 & delta1d & 15.261 & 0.00125857 & 4.79625e-05 & 3,4,9,10 &  \\

## Grating Analysis

- Results from Sherpa fits; notice the low background in the HEG data, and the higher resolution.

## Guide

- GUIDE is a collection of scripts which access the atomic database ATOMDB.  GUIDE provides a number of informational functions:

| identify | Print finding chart of wavelengths |
| --- | --- |
| strong | List strong lines at a given temperature |
| describe | Describe atomic parameters of a line |
| mdl2latex | Convert fit parameters into a latex table |
| ionbal | Output ionization balance fractions for a given ion |

- GUIDE works in Sherpa or Chips; initialize it with **import(“guide”)**

## WebGUIDE

## WebGUIDE

## WebGUIDE

## Grating Analysis

- When faced with the task of understanding a high-resolution spectrum, especially a line-dominated one, it is likely that simple equilibrium models will not work adequately, despite possibly giving low reduced  2 values--usually because the counts/bin is low and so the errors are overestimated, not because the fits are good.
- In this case, some new strategies are needed.  Here are some suggestions for starting points.
**First**: Determine if your plasma is dominated by photoionization or collisional ionization.  For example, the initial analysis of ASCA\ data of Cygnus X-3 used a collisional model, even though the emission is due to photoionization (see Liedahl & Paerels 1996)
- Make a list of all processes that could be affecting line emission. For example, in helium-like ions, the three dominant lines are the resonance, the forbidden, and the intercombination lines.  These can be excited by direct excitation, radiative recombination, dielectronic recombination of hydrogen-like ions, innershell ionization of lithium-like ions,  cascades from higher levels, or by photoexcitation or photoionization.  Once created, the lines can be absorbed or scattered by like ions or by different ions.  In many cases, simple
- physical arguments can be used to limit or exclude various processes, reducing the parameter space that must be searched.

## Grating Analysis

**Second**: Attempt to determine the true continuum level with confidence.  The continuum in a hot plasma is not necessarily dominated by bremsstrahlung:

- Weak emission lines will also blend in to make the continuum seem larger, which will lead to a systematic underestimate of line fluxes, and therefore elemental and ionic abundances.  After you have found an acceptable spectral model, search for regions in the model with no or few lines, and compare the model to the data in this region.  If the model continuum overestimates the continuum here, it is likely it overestimates it everywhere due to unresolved line emission

## AstroAtom

## Conclusions

-  Reprocessing grating data is no longer absolutely required, but has gotten far easier and provides a sense of confidence about the data.
-  Co-Adding and/or binning grating data should be avoided when possible.  Remember that, statistically, nothing is gained by it, although it may be much faster to fit it and easier to see the results.
-  A number of new facilities for atomic data analysis have been created for Sherpa and XSPEC However, remember to check the caveats on this data before trusting it totally!  For the ATOMDB, they are at _[http://asc.harvard.edu/atomdb/doc/caveats.html](http://asc.harvard.edu/atomdb/doc/caveats.html)_
-  Global fitting of generic equilibrium may be useful for guiding the analysis, but any project should begin with a physics-based approach, followed ideally by a line-based analysis and finally by checking regions which should well-understood (such as line-free areas or those dominated by a single line).