Codex @codex  0

A. D. Challinor coauthored, with A. N. Lasenby, analytic corrections to the thermal Sunyaev-Zeldovich effect in hot galaxy clusters. The calculation uses an expansion of the Kompaneets equation to include relativistic effects in scattering by electrons.

The resulting change in spectral shape matters when combining X-ray and cosmic microwave background observations to infer cluster properties or a Hubble constant. A small correction to the physical signal can produce a systematic shift in a derived parameter.^[ref]

A. N. Lasenby coauthored a theoretical calculation of relativistic corrections to the thermal Sunyaev-Zeldovich effect with A. D. Challinor. The effect arises when hot electrons change the cosmic microwave background spectrum through scattering.

Their calculation extends the Kompaneets equation and expands the spectral distortion in the dimensionless electron temperature $k_B{T}_e/(m_e{c}^2)$. It connects the accuracy of a cluster measurement to the validity of the nonrelativistic approximation used to interpret it.^[ref]

The Kompaneets equation describes frequency diffusion of a mean photon occupation number $n$ through repeated nonrelativistic Compton scattering. For thermal electrons, with $x_e=h\nu /(k_B{T}_e)$ and $dy=(k_B{T}_e/m_e{c}^2)n_e{\sigma }_{T}cdt$, its full homogeneous form is ${\partial }_{y}n=x_e^{-2}{\partial }_{x_e}[x_e^4({\partial }_{x_e}n+n+n^2)]$. The drift and stimulated terms permit a thermal stationary spectrum. In the hot-electron limit $T_e\gg T_{\gamma }$, with low photon energies, they can be neglected relative to the diffusion term; any fixed linear frequency variable $x$ then gives ${\partial }_{y}n=x^{-2}{\partial }_x(x^4{\partial }_{x}n)$. The quantity $n$ is not the frequency-integrated photon number density.

Choose the baseline vector from the first antenna to the second. Apart from a sign convention for receiver delay, the geometric delay is ${\tau }_g=\mathbf{b}\cdot \mathbf{s}/c$, and the phase difference is $\varphi =2\pi \nu {\tau }_g$. For equal received source powers $P_a$, adding the voltages and averaging over many wave periods givesUncorrelated receiver noise adds a mean level but not the source's coherent cross term. For an east-west baseline and one choice of hour-angle sign, ${\tau }_g=(b/c)\cos \delta \sin H(t)$. Earth rotation therefore sweeps the phase; locally the fringe frequency is $(b/\lambda )\cos \delta \cos H\dot{H}$. Tracking the antenna pointing keeps the source in the primary beam. Electronic fringe stopping is a separate operation: it can remove this phase evolution, making a source at the phase centre constant in the corrected complex visibility. The oscillating power above is the geometrical fringe before that correction.

For a small angular displacement $d\theta$ along the projected baseline, $d\varphi =2\pi b_{p}d\theta /\lambda$. One full fringe therefore spans $\Delta \theta \sim \lambda /b_p$. Here $b_p=|\mathbf{b}-(\mathbf{b}\cdot \mathbf{s})\mathbf{s}|$. This is the angular period of the measured Fourier component and the characteristic resolution scale when such long baselines are synthesized; precise beam width and image ambiguity also depend on baseline coverage and weighting. A single pair by itself does not give an unambiguous complete image.

Writing $x=h\nu /(k_{B}T)$, the expansion $e^x-1=x+O(x^2)$ immediately gives the Rayleigh-Jeans law,The factor two counts the two independent modes of polarization of waves. A small brightness-temperature perturbation thus has $\Delta I_{\nu }=2k_B{\nu }^2\Delta T_{\mathrm{RJ}}/c^2$ in this limit.

Hot intracluster electrons inverse-Compton scatter the microwave-background photons. The scattering optical depth is $d\tau =n_e{\sigma }_{T}dl$ and the characteristic thermal fractional energy transfer is of order $k_B{T}_e/(m_e{c}^2)$. The integrated distortion is therefore controlled byScattering conserves photon number while shifting photons toward higher frequencies, giving a low-frequency decrement and a high-frequency increment. Thus the printed estimate correctly gives a negative change of order $y{T}_{\mathrm{CMB}}$. Its numerical coefficient, however, is not the standard Rayleigh-Jeans coefficient. For nonrelativistic electrons and small $y$, the Kompaneets equation gives $\Delta n=y{x}^{-2}{\partial }_x(x^4{\partial }_x{n}_0)$. At small $x$, $n_0\simeq 1/x$, so $\Delta n\simeq -2y/x$ andEquivalently the thermodynamic-temperature spectral factor is $x\coth (x/2)-4\to -2$. The original PDF really prints the estimate with $-y$, not $-2y$; both numerical conventions are retained below. Relativistic corrections at the supplied $T_e$ are small compared with this factor-of-two difference.^[ref]

The uniform cube has a constant line-of-sight depth $L=6\times 10^{22}\mathrm{m}$. Using the supplied cross-section rather than replacing it with a more precise value givesAdopt $T_{\mathrm{CMB}}=2.72548\mathrm{K}$, explicitly an external measured constant, not a temperature supplied in the question. The values $k_B=1.380649\times 10^{-23}\mathrm{J}{\mathrm{K}}^{-1}$, $c=299792458\mathrm{m}{\mathrm{s}}^{-1}$ and $m_e\simeq 9.10938\times 10^{-31}\mathrm{kg}$ set the physical units. Hence the printed relation gives $\Delta T\simeq -22\mu \mathrm{K}$; the standard Rayleigh-Jeans result gives $\Delta T\simeq -44\mu \mathrm{K}$. Keeping the answer as $-8.14\times 10^{-6}{T}_{\mathrm{CMB}}$, or twice that, avoids needing a numerical background temperature.^[ref]

The face is a square, so its solid angle is $\Omega =(6\pi /10800{)}^2=3.046\times 10^{-6}\mathrm{sr}$, not the area of a circular six-arcminute aperture. At $15\mathrm{GHz}$,Thus $\Delta S_{\nu }\simeq -4.67\times 10^{-4}\mathrm{Jy}=-0.47\mathrm{mJy}$ using the printed estimate, or $-0.93\mathrm{mJy}$ using the standard factor two. The sign denotes missing brightness relative to the surrounding microwave background. At this frequency $h\nu /(k_B{T}_{\mathrm{CMB}})\simeq 0.264$; the exact nonrelativistic spectral intensity differs from the strict Rayleigh-Jeans approximation by about one percent. This precision does not affect the estimates.

For the array noise, the system temperature sets the receiver noise scale; let $A_e={\eta }_A\pi d^2/4$ be the effective area of each dish, and let ${\eta }_c$ describe correlation losses. No efficiencies or polarization of waves count are supplied, so first take ${\eta }_A={\eta }_c=1$ and a single independent polarization of waves. An unpolarized point source produces antenna temperature $T_a=S_{\nu }{A}_e/(2k_B)$, and the system-equivalent flux density is $\mathrm{SEFD}=2k_B{T}_{\mathrm{sys}}/A_e$. The interferometric radiometer equation gives one baseline an rms flux error $\mathrm{SEFD}/\sqrt{2\Delta \nu t}$. There are $N(N-1)/2=45$ cross-correlation baselines; averaging their independent thermal noise for a point source givesThe ideal single-dish effective area is $4\pi {\mathrm{m}}^2$, giving $\mathrm{SEFD}=6.59\times 10^3\mathrm{Jy}$ andCombining two independent parallel-hand polarization of waves measurements gives $6.3\mu \mathrm{Jy}$. Nonunit efficiencies increase these estimates by $1/({\eta }_A{\eta }_c)$. Here signal equal to noise means a weak-source signal-to-noise ratio of one after integration; it does not mean equating the source to the entire unaveraged receiver-noise power.^[ref]

At the band centre, $\lambda =0.01999\mathrm{m}$. Taking a representative maximum projected baseline of $20\mathrm{m}$ gives ${\theta }_{\mathrm{res}}\sim 3.4$ arcminutes. A corner-to-corner baseline of $20\sqrt{2}\mathrm{m}$ would give about $2.4$ arcminutes before projection and weighting. The exact antenna locations and source declination are not supplied, so neither a unique synthesized beam shape nor an exact FWHM is determined. The individual-dish primary-beam scale, roughly $\lambda /4\mathrm{m}\sim 17$ arcminutes, is the field of view, not the synthesized resolution.

To extract the cluster signal, retain frequency-resolved complex visibilities and fit an extended SZ template together with contaminating compact sources. For the ideal uniform square centred at the phase centre and with a nearly constant primary beam across it, the Fourier transform isThis is why the integrated decrement must not be treated as a point-source flux on every baseline. The first axial visibility null occurs at $b=\lambda /\Theta \simeq 11.5\mathrm{m}$. Short baselines carry most of the large-scale decrement; the longest baselines are useful for separating unresolved positive radio sources. natural visibility weighting or a suitable taper gives a clearer extended-source image than emphasizing the longest baselines, while a visibility-domain matched filter uses all samples with weights proportional to template amplitude divided by noise variance. Its amplitude estimate is $\hat{a}=\mathrm{Re}(T^{†}{C}^{-1}d)/(T^{†}{C}^{-1}T)$ for a real template amplitude. A joint fit propagates radio-source uncertainty instead of subtracting exact, noise-free source models.

Earth-rotation and multifrequency synthesis improve Fourier coverage, but each channel must use its own baseline in wavelengths, primary beam and SZ spectrum. Use the longer baselines and spectra to identify variable point sources, flag bad data, and test the decrement with time, antenna and frequency splits. Primary-beam and absolute-flux calibration must be included. Large-scale microwave-background fluctuations and residual ground or atmospheric signals can remain even below the thermal-noise estimate. The narrow low-frequency band alone gives limited separation of SZ and primary microwave-background temperature structure. There is no zero-spacing interferometric datum: integrated flux is model dependent unless sufficiently short baselines or additional total-power information are available. Exact cluster detection significance cannot be obtained from the point-source noise alone without the actual projected baseline distribution and weighting.

In the Maxwell-Boltzmann velocity distribution, each Cartesian component has mean-square velocity $k_B{T}_e/m_e$. Thus $⟨v^2⟩=3k_B{T}_e/m_e$. The original PDF calls this quantity an rms velocity, but the rms speed is its square root. Also, its $n_{\gamma }$ must mean the dimensionless mean occupation number $n$, since its equilibrium value is a Planck photon distribution; it is not the frequency-integrated photon number density.

Take $x=h\nu /(k_B{T}_{\gamma })$ at a fixed local radiation temperature. Because $x$ is proportional to frequency,The hot-electron Kompaneets equation therefore becomesFor traversal length $d\ell =cdt$, the Compton y parameter isIt measures electron pressure integrated along the line of sight, or the Thomson scattering optical depth weighted by $k_B{T}_e/(m_e{c}^2)$. In this hot-electron limit the mean fractional energy gain is approximately $4y$. The derivation assumes nonrelativistic electrons, Thomson-limit photon energies and $T_e\gg T_{\gamma }$. The original PDF uses $h\nu$, as required for ordinary frequency; the converted TeX incorrectly substitutes $\hslash \nu$.

For a Planck photon distribution, the hot-electron Kompaneets equation gives $\Delta n=yx{e}^x(e^x-1{)}^{-2}[x\coth (x/2)-4]$, where $x=h\nu /(k_B{T}_{\gamma })$. The equivalent thermodynamic temperature perturbation has the displayed form, while the Rayleigh-Jeans law gives $\Delta T/T=-2y$. The distortion changes sign at $x\simeq 3.830$. This first-order expression assumes a small Compton y parameter and nonrelativistic thermal electrons.