X-ray astronomy

X-ray astronomy infers the properties of celestial sources from their X-rays. It requires instruments above the strongly absorbing layers of Earth's terrestrial atmosphere. Riccardo Giacconi's 2002 Nobel lecture connects the early rocket detections to later grazing-incidence X-ray optics, which collect photons over a large entrance area onto a small detector.^[ref]

The measurement is an inverse problem: a count rate mixes incident photon arrival rate, detector response, and background. A point-like astronomical source produces a direction-dependent excess as the instrument scans; diffuse X-ray background contributes over many directions. Neither a count excess nor a narrow feature by itself determines the emission mechanism. The useful progression is from detecting an excess, to locating it, to measuring its spectrum and light curve.

In 1962, Riccardo Giacconi and colleagues flew three Geiger counters to search for lunar X-ray fluorescence. Different mica thicknesses tested penetration; lampblack blocked ultraviolet. Each counter sat in an anticoincidence detector using a scintillation counter to reject cosmic rays. There was no mechanical collimator. Rotation scanned the sky; optical sensors supplied orientation.^[ref]

Counts were sorted by azimuth; a faulty counter and tumbling descent data were excluded. A directional peak lay above a diffuse component. Window, altitude and angular responses tested X-rays against charged particles; the Moon's nondetection checked optical leakage. The authors argued for an origin outside the Solar system, while residual cosmic ray background was not directly measured and source size, position and spectrum remained assumption-dependent.^[ref]

A simple original forward problem makes the distinction precise. Let $\varphi$ be the scan azimuth, let $g(\varphi )$ be a dimensionless directional response with maximum one, let $F$ be a source's band-integrated photon arrival rate per unit area, let $A_0$ be its on-axis X-ray effective area, and let $b$ be a constant accepted background rate. The expected rate isFor a broad illustrative response $g(\varphi )=\max (0,\cos \varphi )$, the full width at half maximum is $120^{\circ }$: $\cos \varphi =1/2$ at $\varphi =\pm 60^{\circ }$. For a one-dimensional ideal collimator model $g(\varphi )=\max (0,1-|\varphi |/{\theta }_c)$, the full width at half maximum is ${\theta }_c$. Choosing ${\theta }_c=20^{\circ }$ narrows the response sixfold without focusing any additional photons. These are teaching response functions, not reconstructions of the rocket's window transmission or geometry. A narrower scan peak can reflect the instrument response even when the source is unresolved.

Scorpius X-1 is the celestial X-ray source associated with the 1962 discovery, subsequently located more precisely. The discovery scan did not supply a resolved image of this X-ray binary.^[ref]

Diffuse X-ray background is X-ray emission spread over the instrument's observed sky rather than assigned to a single resolved source. Diffuse appearance at finite angular resolution does not prove emission by a continuous medium: many unresolved point-like astronomical sources can also give a smooth signal. Instrumental background is a separate contribution, even though both enter measured counts.

X-ray effective area is the area that converts an incident photon arrival rate per unit area into an accepted count rate. For a fixed direction and a sufficiently narrow energy interval, writeHere $A_{\mathrm{geom}}$ is projected geometric area, $T$ is transmittance through windows or optics, $\eta$ is quantum efficiency, ${ϵ}_s$ is the probability that a detected source event passes event selection, and $t$ is live observing time. $N_s$ is an expected value, not a guarantee of an integer count. Multiplication follows by counting incident photons over area and time, then applying the successive conditional acceptance probabilities. Live time already excludes intervals when the instrument cannot record events; multiplying by a live-time fraction again would count the same loss twice.

For an invented narrow-band example, choose $A_{\mathrm{geom}}=40\mathrm{c}{\mathrm{m}}^2$, $T=0.8$, $\eta =0.75$, and ${ϵ}_s=5/6$. Then $A_{\mathrm{eff}}=20\mathrm{c}{\mathrm{m}}^2$. A source with $F=0.30\mathrm{photons}\mathrm{c}{\mathrm{m}}^{-2}{\mathrm{s}}^{-1}$ observed for $t=100\mathrm{s}$ gives $N_s=600$. $F$ here counts photons; converting it into observed radiative flux requires their energies. If all have energy $E_0$, the corresponding radiative flux is $E_{0}F$.

A broad spectrum instead requireswhere $F_E$ is the differential photon arrival rate per area and energy. A single band-averaged X-ray effective area depends on the assumed spectrum. This is why a count rate alone cannot give an assumption-free radiative flux.

Background subtraction must account for uncertainty in the background measurement itself. In an original idealized experiment, an on-source exposure contains source and background, while an independent off-source exposure contains background alone. Assume stationary rates and identical background acceptance per unit live time. With $\alpha =t_{\mathrm{on}}/t_{\mathrm{off}}$, let independent counts obey Poisson distributions:Then linearity of expectation and addition of variances for independent random variables giveThe result concerns accepted source counts $S$, after all efficiencies. For large counts the normal approximation gives the descriptive signal-to-noise ratio in photon countingThis plug-in ratio is not an exact null-hypothesis significance. Small counts require the discrete Poisson model; searching many sky directions also changes a detection claim's p-value.

Choose $t_{\mathrm{on}}=100\mathrm{s}$ and $t_{\mathrm{off}}=200\mathrm{s}$, with observed counts $N_{\mathrm{on}}=1000$ and $N_{\mathrm{off}}=800$. Thus $\alpha =1/2$, $\hat{S}=600$, and the estimated standard deviation is $\sqrt{1200}=34.64$ counts. The inferred $F=\hat{S}/(A_{\mathrm{eff}}{t}_{\mathrm{on}})=0.300\pm 0.0173\mathrm{photons}\mathrm{c}{\mathrm{m}}^{-2}{\mathrm{s}}^{-1}$ uses the preceding $A_{\mathrm{eff}}=20\mathrm{c}{\mathrm{m}}^2$ and includes counting uncertainty only. The descriptive ratio is $17.32$. Treating the scaled off-source background as exact would omit 200 counts squared from the estimated variance.

For a selection-design comparison, suppose before a veto the expected source contribution is $S_0=750$, the on-source background is $B_0=2000$, and $\alpha =1/2$. Let source acceptance be ${ϵ}_s=0.8$ and background acceptance be ${ϵ}_b=0.2$. Then $S=600$, $B=400$, and $S/\sqrt{S+(1+\alpha )B}=17.32$, whereas omitting the veto gives $750/\sqrt{750+3000}=12.25$. Rejection helps because the background is reduced enough to compensate for lost source events. A different veto can make matters worse; source efficiency and background rejection are separate probabilities. This comparison is an independent teaching example, not the efficiency calibration used above.

If both exposure times grow with fixed $\alpha$ and stationary rates, $S=s{t}_{\mathrm{on}}$ and $B=b{t}_{\mathrm{on}}$ imply $Z\simeq s\sqrt{t_{\mathrm{on}}}/\sqrt{s+(1+\alpha )b}$. Four times the live time doubles this ideal ratio. Misestimated background rates, uncertain X-ray effective area, and source variability can defeat that scaling; longer exposure does not remove systematic error.

Grazing-incidence X-ray optics redirect X-rays by shallow specular reflections along mirror surfaces. Riccardo Giacconi's lecture describes coaxial paraboloidal and hyperboloidal mirrors using two reflections, and the later Einstein observatory's imaging capability. Focusing separates collecting area from detector area: the collected source signal can remain large while the detector background is drawn from a smaller region.^[ref]

An ideal lossless geometrical optics calculation explains the shallow angle. Take vacuum refractive index one and a material with real refractive index $n=1-\delta$, where $0<\delta \ll 1$. Let $\theta$ be the grazing angle measured from the surface. Snell's law written for angles from the normal requires $\cos \theta =n\sin {\theta }_t$. A propagating transmitted ray exists only if $\cos \theta \le n$. Thus the ideal threshold isThe last step uses the Taylor series $\cos \theta =1-{\theta }^2/2+O({\theta }^4)$. For a chosen $\delta =2\times 10^{-5}$, ${\theta }_c\simeq 0.006325$ radians, or $0.3624^{\circ }$. Absorption makes the refractive index complex and softens this ideal threshold; surface roughness and alignment errors also reduce throughput or worsen the point spread function. These are reasons to calibrate X-ray effective area and angular resolution, rather than assume that reflection is perfect.

The benefit can be quantified without specifying a historical telescope. Suppose a source retains fraction $f$ of its counts inside an extraction region, while a uniform detector background contributes $qat$ counts over detector area $a$. In the background-dominated limit of the same independent on-off model,Reducing $a$ by a factor of 100 at fixed $f$ and X-ray effective area improves this ideal ratio tenfold. A collimator restricts admitted directions; focusing optics bring selected rays together. An image also separates nearby sources and identifies clean background regions, connecting optical design directly to a better constrained inverse problem.