Radio Interferometry Basics U Klein 1 Fourier optics
Radio Interferometry Basics U. Klein 1. Fourier optics 2. Fourier theorems 3. Visibility function 4. Aperture synthesis 5. Coordinates 6. Correlation interferometer 7. Imaging 8. Cleaning GLOW Interferometry School, Bielefeld 2011
Fundamental understanding of the functionality of radio interferometry GLOW Interferometry School, Bielefeld 2011
1. Fourier optics Aperture with constant extent perpendicular to screen, with extent a in this y-direction, constant current Gy(x) in y-direction. Infinitesimal far-field element d. E resulting from an E-field component E(x) along infinitesimal element dx and integrating over y and x leads to This is Huygens’ Principle for the far-field: superposition of field components along the aperture Far-field: R x (… ties up to Jackson, Chapt. 14) GLOW Interferometry School, Bielefeld 2011
Normalising by wavelength Lumping constants and R dependence into leads to simple Fourier intergral (x = x / , direction cosine ξ = sin ) GLOW Interferometry School, Bielefeld 2011
slit or uniformly illuminated aperture: sinc function circular aperture: Bessel function antenna diagram: P( , ) = |E( , )|2 GLOW Interferometry School, Bielefeld 2011
2. Fourier theorems Fourier transform F(s) of function f(x) with inverse transform GLOW Interferometry School, Bielefeld 2011
Convolution: GLOW Interferometry School, Bielefeld 2011
Consider the two domains: x time t [s] s frequency application [Hz] acoustics Oboe: waves angle [rad] Oboe: spectrum spatial frequ. u [rad-1] radio astronomy LMC: radio brightness LMC: spatial frequ. spectrum GLOW Interferometry School, Bielefeld 2011
Let F f(x) = F(s), F g(x) = G(s) addition theorem similarity theorem example: an aperture has beam size /D GLOW Interferometry School, Bielefeld 2011
shift theorem example: electronic steering GLOW Interferometry School, Bielefeld 2011
convolution theorem example: beam smearing convolution with broader beam multiplication with narrower ‘filter‘ GLOW Interferometry School, Bielefeld 2011
autocorrelation theorem example: antenna diagramme GLOW Interferometry School, Bielefeld 2011
convolution and auto-/cross-correlation theorem: reception pattern of an interferometer GLOW Interferometry School, Bielefeld 2011
sampling theorem: how to properly convert continuous into discretised functions comb or „Sha“ function: multiplication f(x) · III(x) converts continuous into discretised function GLOW Interferometry School, Bielefeld 2011
A function whose Fourier transform is zero for |s| > sc is fully specified by values spaced at equal intervals not exceeding ½ sc-1 save for any harmonic term with zeros at the sampling points. ‘aliasing‘! GLOW Interferometry School, Bielefeld 2011
GLOW Interferometry School, Bielefeld 2011
3. Visibility function Correlated power received from source with brightness B where (D = D/ ) D is the baseline. Integration yields GLOW Interferometry School, Bielefeld 2011
With it is convenient to introduce complex visibility such that GLOW Interferometry School, Bielefeld 2011
Comparing with 1 st version of P we arrive at Since AN 1 (the single-antenna diagramme is very broad), and writing we obtain This is the basic relation for aperture synthesis: In order to retrieve the brightness B(ξ, η), we need to measure the visibility V(u, v) for as many baselines Dλ= (u, v) as possible and Fourier-transform it. GLOW Interferometry School, Bielefeld 2011
example: double source (schematically): apply convolution theorem source brightness = convolution of -function with two -functions visibility = multiplication of cosine with sinc-function Now observe with different baselines GLOW Interferometry School, Bielefeld 2011
astronomical telescope = spatial filter: single dish: low-pass filter interferometer: high-pass filter GLOW Interferometry School, Bielefeld 2011
4. Aperture synthesis need many D = (u, v) consider 4 4 apertures, producing n (n-1) 2 = 120 baselines, most of which are redundant, while we need only 24 of them to cover the (u, v) plane aim: preferably little redundancy! GLOW Interferometry School, Bielefeld 2011
Earth-rotation synthesis need good sampling of (u, v) plane with many (nonredundant) baselines earth-rotation synthesis helps a lot due to Sir Martin Ryle (Nobel Prize 1974 with A. Hewish) GLOW Interferometry School, Bielefeld 2011
east-west array: e. g. WSRT 0 1 144 m 2 144 m 3 4 9 A B C D 36 m 9 144 m 300 m GLOW Interferometry School, Bielefeld 2011 9 144 m 180 m
earth rotation produces a good transfer function: (u, v) tracks e. g. pure east-west interferometer FT antenna diagramme („dirty beam“) GLOW Interferometry School, Bielefeld 2011
for the kth ring, the grating rings have functional form where sinc 1/2 implies the “half-order derivative“ of the sinc function derívative theorem: and from this: (intrigued? Bracewell & Thompson, Ap. J. 188, 77, 1973) GLOW Interferometry School, Bielefeld 2011
actually need max. 12 hrs. to obtain full (u, v) coverage VLA: 8 hrs. (Y-shape) V(u, v) is Hermitean, since the brightness is a real quantity. … we get the other half for free! GLOW Interferometry School, Bielefeld 2011
5. Coordinates Complex visibility reads Define coordinate system (u, v, w), which are related to unit vectors and baseline: We may have to deal with non-zero w-component! GLOW Interferometry School, Bielefeld 2011
, , are direction cosines and the visibility becomes GLOW Interferometry School, Bielefeld 2011
Coordinate system ( , ) corresponds to projection of celestial sphere onto a tangential plane at the origin of the field centre. possible projections: • tan: optical astronomy • arc: Schmidt telescopes, single dishes • sin : aperture synthesis (more: Calabretta & Greisen, A&A 395, 1077, 2002) GLOW Interferometry School, Bielefeld 2011
direction cosines and their relation to other coordinate systems: spherical coordinates , equatorial coordinates , GLOW Interferometry School, Bielefeld 2011
Using polar coordinates , , the solid-angle element d reads The releation between the infinitesimal elements is given by where |J| is the determinant of the Jacobi matrix which leads to GLOW Interferometry School, Bielefeld 2011
Antenna coordinates and baseline components defines a double ellipse: GLOW Interferometry School, Bielefeld 2011
Coverage of (u, v) plane e. g. VLA e. g. LOFAR a bit weird: component LZ gives rise to non-coplanar baselines not simple Fourier integral anymore! GLOW Interferometry School, Bielefeld 2011
6. Correlation interferometer has power response with phase lag due to geometric time delay GLOW Interferometry School, Bielefeld 2011
heterodyne principle: down-convert from high (RF or HF) to intermediate frequency (IF) produces frequency spectrum including signal frequ. (USB) image frequ. (LSB) GLOW Interferometry School, Bielefeld 2011
Effect of finite bandwidth (1): Delay compensation i indispensible in order to avoid ‚fringe washing‘. We measure correlation of two voltages V 1, V 2: T By definition, this is also the Fourier transform of the bandpass |H( )| Calculate power over the band , assuming rectangular bandpass and (roughly) constant power as a function of frequency: GLOW Interferometry School, Bielefeld 2011
Hence finally: So the correlated signal is modulated with a sinc function! = 50 MHz D = 1 km tracking source over 1. 6‘ without delay compensation produces ~1% of loss needs g 1. 6 ns GLOW Interferometry School, Bielefeld 2011
the phases from telescopes 1 and 2 are hence the correlated power is where = g - i and ‘ ‘ refers to USB and LSB. Measure complex visibility with quadrature network: GLOW Interferometry School, Bielefeld 2011
Fringe rotation/stopping: earth‘s rotation modulates correlated power quasisinusoidally with natural fringe rate where e = dh dt = angular rotation frequency of earth, = declination. At the VLA, the fringe frequency may exceed 150 Hz. Not interested in fringe rate, but rather in change of V; hence compensate fringe rotation by modulating LO phase. In the IF section, we have So control LO phase such that vanishes. This now requires a complex correlator in order to measure the amplitude and phase separately! GLOW Interferometry School, Bielefeld 2011
Effect of finite bandwidth (2): Fourier integral is precise only for monochromatic radiation! At the centre frequency 0 we have the precise relation Away from the centre frequency we have Generalised similarity theorem in n dimensions: GLOW Interferometry School, Bielefeld 2011
application to 2 -D Fourier relation between visibility and brightness: so that (“~“ indicating the influence of the bandpass) Delay tracking is exact for 0 = 0 at = 0. For signals at frequeny arriving from direction ( , ) there will be a lag error so that the phase is imprecise by an amount GLOW Interferometry School, Bielefeld 2011
Consider point source at location 0, assume rectangular bandpass. Then Fourier kernel V(u) where we have assumed ( 0)2 1. Thence, with du = ( 0) du 0, we obtain the convolution The Fourier transform of the sinc function is a box function where This effect called bandwidth smearing (chromatic aberration) is irreversible! GLOW Interferometry School, Bielefeld 2011
Short version of this: radiation within range of frequencies produces u-v tracks of different size produces images of different scale if referred to same frequency when doing the Fourier transform. v u re t n ld ‘ 13 to fie ce radial smearing GLOW Interferometry School, Bielefeld 2011
7. Imaging In practice: convert Fourier integral into a sum and apply an FFT. This requires a retabulation of the measured visibilities onto a 2 N 2 M grid. In what follows, we use the integral notation for the sake of better readability of the equations. The true visibilities are strongly modified owing to: 1. sampling by transfer function: S(u, v) 2. regridding onto regular grid for FFT with grid function: G(u, v) 3. weighting to shape the synthesized beam and control sensitivity: R(u, v) 4. FFT GLOW Interferometry School, Bielefeld 2011
transfer function: S(u, v) gives rise to ‘dirty beam‘ weighting with W(u, v) Plateau de Bure Interferometer W = R T D R: Aeff, Tsys, , T: taper (shapes the beam) D: density of visibilities brightness (‘dirty image‘) now is GLOW Interferometry School, Bielefeld 2011
now prepare for FFT by retabulation S(u, v); needs two steps: (i) convolve with appropriate function C(u, v) (ii) multiply with grid function G(u, v) (fakir‘s bed of nails) this is a significant modification of the data use convolution theorem to see modification in the image domain: GLOW Interferometry School, Bielefeld 2011
can ‘undo‘ gridding by dividing VD by the inverse transform of G: , yielding the ‘grid-corrected‘ dirty image: This is the true brightness distribution, convolved with the dirty beam GLOW Interferometry School, Bielefeld 2011
8. Cleaning Goal: deconvolution in order to get rid of sidelobes GLOW Interferometry School, Bielefeld 2011
basic clean algorithm (Högbom 1974): (i) First compute the dirty map BD( , ) via FT of VD(u, v) and the dirty beam FD( , ) via the FT of S(u, v) · W(u, v). (ii) Find position of maximum brightness and subtract ·FD( , ) · Bmax; “loop gain” controls speed of the cleaning process (0 < < 1); store Bmax at the corresponding position in a work array (creates “fakir’s bed of nails”). (iii) Go to (ii) until a user-defined level is reached (e. g. number of iterations or rms noise). (iv) Convolve the so obtained “point source model” (“fakir’s bed of nails”) with the “clean beam”, or synthesized beam, e. g. a Gaussian of width ~ / Dmax. (v) Add this convolved model to the residual map to yield the clean image BC( , ). There are variants with major (check in the u-v domain) and minor cycles (Clark 1980, Cotton-Schwab 1984) GLOW Interferometry School, Bielefeld 2011
GLOW Interferometry School, Bielefeld 2011
Zero-spacing problem Missing short spacings: interferometer measures zero total flux! hence since which cannot be measured “negative bowl” GLOW Interferometry School, Bielefeld 2011
need to add in single-dish measurements (in the u, v domain!) GLOW Interferometry School, Bielefeld 2011
Don‘t forget that the brightness distribution is multiplied by the enveloping ‘primary beam‘, i. e. the beam of the individual antennae. Need to correct for this degradation by dividing the clean image by the primary beam: Naturally gives rise to increasing noise at radial distances larger than the primary HPBW. GLOW Interferometry School, Bielefeld 2011
3 C 120 GLOW Interferometry School, Bielefeld 2011
Th a nk sf or yo ur att en ti on ! GLOW Interferometry School, Bielefeld 2011
- Slides: 56