the development of high-resolution imaging in radio astronomy jim moran harvard-smithsonian center...
Post on 15-Jan-2016
218 views
TRANSCRIPT
![Page 1: The Development of High-Resolution Imaging in Radio Astronomy Jim Moran Harvard-Smithsonian Center for Astrophysics 7th IRAM Interferometry School, Grenoble,](https://reader036.vdocument.in/reader036/viewer/2022062518/56649d3e5503460f94a17549/html5/thumbnails/1.jpg)
The Development of High-Resolution Imaging
in Radio AstronomyJim Moran
Harvard-Smithsonian Center for Astrophysics
7th IRAM Interferometry School, Grenoble, October 4–7, 2010
It is an honor to give this lecture in the city where Joseph Fourier did the work that is so fundamental to our craft.
![Page 2: The Development of High-Resolution Imaging in Radio Astronomy Jim Moran Harvard-Smithsonian Center for Astrophysics 7th IRAM Interferometry School, Grenoble,](https://reader036.vdocument.in/reader036/viewer/2022062518/56649d3e5503460f94a17549/html5/thumbnails/2.jpg)
Outline of Talk
I. Origins of Interferometry
II. Fundamental Theorem of Interferometry
(Van Cittert-Zernike Theorem)
III. Limits to Resolution (uv plane coverage)
IV. Quest for High Resolution in the 1950s
V. Key Ideas in Image Calibration and Restoration
VI. Back to Basics – Imaging Sgr A* in 2010 and beyond
![Page 3: The Development of High-Resolution Imaging in Radio Astronomy Jim Moran Harvard-Smithsonian Center for Astrophysics 7th IRAM Interferometry School, Grenoble,](https://reader036.vdocument.in/reader036/viewer/2022062518/56649d3e5503460f94a17549/html5/thumbnails/3.jpg)
I. Origins of Interferometry
A. Young’s Two-Slit Experiment
B. Michelson’s Stellar Interferometer
C. Basic Radio Implementation
D. Ryle’s Correlator
E. Sea Cliff Interferometer
F. Earth Rotation Synthesis
Albert Michelson (1852–1931)
Martin Ryle (1918–1984)
John Bolton (1922–1993)
Martin Ryle (1918–1984)
Thomas Young (1773–1829)
![Page 4: The Development of High-Resolution Imaging in Radio Astronomy Jim Moran Harvard-Smithsonian Center for Astrophysics 7th IRAM Interferometry School, Grenoble,](https://reader036.vdocument.in/reader036/viewer/2022062518/56649d3e5503460f94a17549/html5/thumbnails/4.jpg)
Young’s Two-Slit Experiment (1805)
1. Move source shift pattern (phase)2. Change aperture hole spacing change period of fringes3. Enlarge source plane hole reduce visibility
FringesBright
Dark
I
X
X
Imax
Imin
I = (E1 + E2 )2 = E1
2 + E22 + 2(E1E2 )
interference term
Imax– Imin
Imax+ IminV =
t 1bandwidth effect
Source plane
Aperture plane(screen with 2 small holes)
Pupil or Image plane
quasi-monochromatic light source
![Page 5: The Development of High-Resolution Imaging in Radio Astronomy Jim Moran Harvard-Smithsonian Center for Astrophysics 7th IRAM Interferometry School, Grenoble,](https://reader036.vdocument.in/reader036/viewer/2022062518/56649d3e5503460f94a17549/html5/thumbnails/5.jpg)
Michelson-Pease Stellar Interferometer (1890-1920)
Image plane fringe pattern.Solid line: unresolved starDotted line: resolved star
Two outrigger mirrors on the Mount Wilson 100 inch telescope
Paths for on axis ray and slightly offset ray
![Page 6: The Development of High-Resolution Imaging in Radio Astronomy Jim Moran Harvard-Smithsonian Center for Astrophysics 7th IRAM Interferometry School, Grenoble,](https://reader036.vdocument.in/reader036/viewer/2022062518/56649d3e5503460f94a17549/html5/thumbnails/6.jpg)
Simple Radio Interferometer
R = I cos
= d s = cos 2
2 d
= e sin2 d
ddt
= sindd
2 d
= 2 = 2 d sin
d sin =
projected baseline
d sin
d
s
R“fringes”
![Page 7: The Development of High-Resolution Imaging in Radio Astronomy Jim Moran Harvard-Smithsonian Center for Astrophysics 7th IRAM Interferometry School, Grenoble,](https://reader036.vdocument.in/reader036/viewer/2022062518/56649d3e5503460f94a17549/html5/thumbnails/7.jpg)
Cygnus ACassiopeia A
Simple Adding Interferometer (Ryle, 1952)
R = (a + b)2 = a 2 +b
2 +2ab
![Page 8: The Development of High-Resolution Imaging in Radio Astronomy Jim Moran Harvard-Smithsonian Center for Astrophysics 7th IRAM Interferometry School, Grenoble,](https://reader036.vdocument.in/reader036/viewer/2022062518/56649d3e5503460f94a17549/html5/thumbnails/8.jpg)
Phase Switching Interferometer (Ryle, 1952)
R = (a + b)2 – (a – b)2 = 4ab
![Page 9: The Development of High-Resolution Imaging in Radio Astronomy Jim Moran Harvard-Smithsonian Center for Astrophysics 7th IRAM Interferometry School, Grenoble,](https://reader036.vdocument.in/reader036/viewer/2022062518/56649d3e5503460f94a17549/html5/thumbnails/9.jpg)
Sea Cliff Interferometer (Bolton and Stanley, 1948)
Response to Cygnus A at 100 MHz (Nature, 161, 313, 1948)
Time in hours
![Page 10: The Development of High-Resolution Imaging in Radio Astronomy Jim Moran Harvard-Smithsonian Center for Astrophysics 7th IRAM Interferometry School, Grenoble,](https://reader036.vdocument.in/reader036/viewer/2022062518/56649d3e5503460f94a17549/html5/thumbnails/10.jpg)
II. Fundamental Theorem
A. Van Cittert-Zernike Theorem
B. Projection Slice Theorem
C. Some Fourier Transforms
![Page 11: The Development of High-Resolution Imaging in Radio Astronomy Jim Moran Harvard-Smithsonian Center for Astrophysics 7th IRAM Interferometry School, Grenoble,](https://reader036.vdocument.in/reader036/viewer/2022062518/56649d3e5503460f94a17549/html5/thumbnails/11.jpg)
Van Cittert–Zernike Theorem (1934)
E1
E2
z r2
r1
x ′
y ′
( = x ′/z)
( = y ′/z)
x=u
y = v
source distribution I (, )
Ground Plane
Sky Plane
E1 = e i r
2
r1
1
E1 E2 = e i (r – r ) 2
r1 r2
21* 2
E2 = e i r
2
r2
2 Huygen’s principle
2 = I(,) r1r2 z 2
Integrate over source
V(u,v) = I(,) e i2 ( u + v) d d
Assumptions1. Incoherent source2. Far field z >d
2 / ; d = 104 km, = 1 mm, z > 3 pc !3. Small field of view4. Narrow bandwidth field =
max
( ) dmax
![Page 12: The Development of High-Resolution Imaging in Radio Astronomy Jim Moran Harvard-Smithsonian Center for Astrophysics 7th IRAM Interferometry School, Grenoble,](https://reader036.vdocument.in/reader036/viewer/2022062518/56649d3e5503460f94a17549/html5/thumbnails/12.jpg)
F (u,v) = f(x,y) e–i 2
(ux + vy) dx dy
F(u,0) = [ f(x,y) dy ] e–i 2
ux dx
F(u,0) fs (x)
Works for any arbitrary angle
Projection-Slice Theorem (Bracewell, 1956)
“strip” integral
Strip integrals, also called back projections, are the common link between radio interferometryand medical tomography.
![Page 13: The Development of High-Resolution Imaging in Radio Astronomy Jim Moran Harvard-Smithsonian Center for Astrophysics 7th IRAM Interferometry School, Grenoble,](https://reader036.vdocument.in/reader036/viewer/2022062518/56649d3e5503460f94a17549/html5/thumbnails/13.jpg)
Moran, PhD thesis, 1968
Visibility (Fringe) Amplitude Functions for Various Source Models
![Page 14: The Development of High-Resolution Imaging in Radio Astronomy Jim Moran Harvard-Smithsonian Center for Astrophysics 7th IRAM Interferometry School, Grenoble,](https://reader036.vdocument.in/reader036/viewer/2022062518/56649d3e5503460f94a17549/html5/thumbnails/14.jpg)
III. Limits to Resolution (uv plane coverage)
A. Lunar Occultation
B. uv Plane Coverage of a Single Aperture
![Page 15: The Development of High-Resolution Imaging in Radio Astronomy Jim Moran Harvard-Smithsonian Center for Astrophysics 7th IRAM Interferometry School, Grenoble,](https://reader036.vdocument.in/reader036/viewer/2022062518/56649d3e5503460f94a17549/html5/thumbnails/15.jpg)
H H
uF.T.
Lunar Occultation
Geometric optics one-dimension integrationof source intensity
MacMahon (1909)
Moon as a knife edge
Im( ) = I( ) H( )
Im(u) = I(u) H(u) where H(u) = 1/u
![Page 16: The Development of High-Resolution Imaging in Radio Astronomy Jim Moran Harvard-Smithsonian Center for Astrophysics 7th IRAM Interferometry School, Grenoble,](https://reader036.vdocument.in/reader036/viewer/2022062518/56649d3e5503460f94a17549/html5/thumbnails/16.jpg)
Criticized by Eddington (1909)
Same amplitude as response in geometric optics, but scrambled phase
F = 2R
(wiggles) (R = earth-moon distance)
H = e –iF u sign(u)2 1
F = 5 mas @ 0.5 wavelength
F = 2´ @ 10 m wavelength´
H
u
H
![Page 17: The Development of High-Resolution Imaging in Radio Astronomy Jim Moran Harvard-Smithsonian Center for Astrophysics 7th IRAM Interferometry School, Grenoble,](https://reader036.vdocument.in/reader036/viewer/2022062518/56649d3e5503460f94a17549/html5/thumbnails/17.jpg)
Occultation of Beta Capricorni withMt. Wilson 100 Inch Telescope and
Fast Photoelectric Detector
Whitford, Ap.J., 89, 472, 1939
![Page 18: The Development of High-Resolution Imaging in Radio Astronomy Jim Moran Harvard-Smithsonian Center for Astrophysics 7th IRAM Interferometry School, Grenoble,](https://reader036.vdocument.in/reader036/viewer/2022062518/56649d3e5503460f94a17549/html5/thumbnails/18.jpg)
Radio Occultation Curves (Hazard et al., 1963)
3C273
![Page 19: The Development of High-Resolution Imaging in Radio Astronomy Jim Moran Harvard-Smithsonian Center for Astrophysics 7th IRAM Interferometry School, Grenoble,](https://reader036.vdocument.in/reader036/viewer/2022062518/56649d3e5503460f94a17549/html5/thumbnails/19.jpg)
Single Aperture
Restore high spatial frequencies up to u = D/
no super resolution
Single Pixel Airy Pattern
Chinese Hat Function
= 1.2 /
DD
F.T.
u
vD
D
![Page 20: The Development of High-Resolution Imaging in Radio Astronomy Jim Moran Harvard-Smithsonian Center for Astrophysics 7th IRAM Interferometry School, Grenoble,](https://reader036.vdocument.in/reader036/viewer/2022062518/56649d3e5503460f94a17549/html5/thumbnails/20.jpg)
IV. Quest for High Resolution in the 1950s
A. Hanbury Brown’s Three Ideas
B. The Cygnus A Story
![Page 21: The Development of High-Resolution Imaging in Radio Astronomy Jim Moran Harvard-Smithsonian Center for Astrophysics 7th IRAM Interferometry School, Grenoble,](https://reader036.vdocument.in/reader036/viewer/2022062518/56649d3e5503460f94a17549/html5/thumbnails/21.jpg)
Hanbury Brown’s Three Ideasfor High Angular Resolution
1. Let the Earth Move (250 km/s, but beware the radiometer formula!)
2. Reflection off Moon (resolution too high)
3. Intensity Interferometer (inspired the field of quantum optics)
In about 1950, when sources were called “radio stars,” Hanbury Brown had several ideas of how to dramatically increase angular resolution to resolve them.
![Page 22: The Development of High-Resolution Imaging in Radio Astronomy Jim Moran Harvard-Smithsonian Center for Astrophysics 7th IRAM Interferometry School, Grenoble,](https://reader036.vdocument.in/reader036/viewer/2022062518/56649d3e5503460f94a17549/html5/thumbnails/22.jpg)
Intensity Interferometry
Fourth-order moment theorem for Gaussian processes
Normal Interferometer
Intensity Interferometer = v12v2
2
= v1v2 V
constant square of visibility
v1v2v3v4 = v1v2 v3v4 + v1v3 v2v4 + v1v4 v2v3
= v12v2
2 = v12 v2
2 + v1v2 2
= P1P2 + V 2
v2
v1
![Page 23: The Development of High-Resolution Imaging in Radio Astronomy Jim Moran Harvard-Smithsonian Center for Astrophysics 7th IRAM Interferometry School, Grenoble,](https://reader036.vdocument.in/reader036/viewer/2022062518/56649d3e5503460f94a17549/html5/thumbnails/23.jpg)
Observations of Cygnus A with Jodrell Bank Intensity Interferometer
Square of Visibilityat 125 MHz
Jennison and Das Gupta, 1952, see Sullivan 2010
![Page 24: The Development of High-Resolution Imaging in Radio Astronomy Jim Moran Harvard-Smithsonian Center for Astrophysics 7th IRAM Interferometry School, Grenoble,](https://reader036.vdocument.in/reader036/viewer/2022062518/56649d3e5503460f94a17549/html5/thumbnails/24.jpg)
Cygnus A with Cambridge 1-mile Telescope at 1.4 GHz
Ryle, Elsmore, and Neville, Nature, 205, 1259, 1965
3 telescopes
20 arcsec resolution
![Page 25: The Development of High-Resolution Imaging in Radio Astronomy Jim Moran Harvard-Smithsonian Center for Astrophysics 7th IRAM Interferometry School, Grenoble,](https://reader036.vdocument.in/reader036/viewer/2022062518/56649d3e5503460f94a17549/html5/thumbnails/25.jpg)
Cygnus A with Cambridge 5 km Interferometer at 5 GHz
Hargrave and Ryle, MNRAS, 166, 305, 1974
16 element E-W Array, 3 arcsec resolution
![Page 26: The Development of High-Resolution Imaging in Radio Astronomy Jim Moran Harvard-Smithsonian Center for Astrophysics 7th IRAM Interferometry School, Grenoble,](https://reader036.vdocument.in/reader036/viewer/2022062518/56649d3e5503460f94a17549/html5/thumbnails/26.jpg)
V. Key Ideas in Image Calibration and Restoration
A. CLEAN
B. Phase and Amplitude Closure
C. Self Calibration
D. Mosaicking
E. The Cygnus A Story Continued
Jan Högbom (1930–)
Roger Jennison (1922–2006)Alan Rogers (1942–)
Ron Eker (~1944–)Arnold Rots (1946–)
several
![Page 27: The Development of High-Resolution Imaging in Radio Astronomy Jim Moran Harvard-Smithsonian Center for Astrophysics 7th IRAM Interferometry School, Grenoble,](https://reader036.vdocument.in/reader036/viewer/2022062518/56649d3e5503460f94a17549/html5/thumbnails/27.jpg)
First Illustration of Clean Algorithm on 3C224.1 at 2.7 GHz with Green Bank Interferometer
J. Högbom, Astron. Astrophys. Suppl., 15, 417, 1974
Zero, 1, 2 and 6iterations
![Page 28: The Development of High-Resolution Imaging in Radio Astronomy Jim Moran Harvard-Smithsonian Center for Astrophysics 7th IRAM Interferometry School, Grenoble,](https://reader036.vdocument.in/reader036/viewer/2022062518/56649d3e5503460f94a17549/html5/thumbnails/28.jpg)
Closure Phase“Necessity is the Mother of Invention”
Observe a Point Source
Arbitrary Source Distribution
d12 + d23 + d31 = 0
C = m + m + m = v + v + v + noise12 23 12 2331 31
N stations baselines, (N – 1)(N – 2) closure conditionsN(N – 1)2
12
12 = d12 s + 1 – 2 2
C = 12 + 23 + 31 = [d12 + d23 + d31] s2
0
fraction of phases f = 1 –
2N N = 27, f 0.9
R. Jennison, 1952 (thesis): MNRAS, 118, 276, 1956; A. Rogers et al., Ap.J., 193, 293, 1974
to so
urce
1
3
3
2
2 d12
“cloud” with phase shift 1 = 2 t
m = v + (i – j) + ijijij
noise
![Page 29: The Development of High-Resolution Imaging in Radio Astronomy Jim Moran Harvard-Smithsonian Center for Astrophysics 7th IRAM Interferometry School, Grenoble,](https://reader036.vdocument.in/reader036/viewer/2022062518/56649d3e5503460f94a17549/html5/thumbnails/29.jpg)
Closure AmplitudeN ≥ 4
Unknown voltage gain factors for each antenna gi (i = 1–4)
![Page 30: The Development of High-Resolution Imaging in Radio Astronomy Jim Moran Harvard-Smithsonian Center for Astrophysics 7th IRAM Interferometry School, Grenoble,](https://reader036.vdocument.in/reader036/viewer/2022062518/56649d3e5503460f94a17549/html5/thumbnails/30.jpg)
A Half Century of Improvements in Imaging of Cygnus A
Kellermann and Moran, Ann. Rev. Astron., 39, 457, 2001
![Page 31: The Development of High-Resolution Imaging in Radio Astronomy Jim Moran Harvard-Smithsonian Center for Astrophysics 7th IRAM Interferometry School, Grenoble,](https://reader036.vdocument.in/reader036/viewer/2022062518/56649d3e5503460f94a17549/html5/thumbnails/31.jpg)
VI. Back to Basics
Imaging Sgr A* in 2010 and beyond
![Page 32: The Development of High-Resolution Imaging in Radio Astronomy Jim Moran Harvard-Smithsonian Center for Astrophysics 7th IRAM Interferometry School, Grenoble,](https://reader036.vdocument.in/reader036/viewer/2022062518/56649d3e5503460f94a17549/html5/thumbnails/32.jpg)
230 GHz Observations of SgrA*
VLBI program led by large consortium led by Shep Doeleman, MIT/Haystack
4030 km
908 km
4630 km
![Page 33: The Development of High-Resolution Imaging in Radio Astronomy Jim Moran Harvard-Smithsonian Center for Astrophysics 7th IRAM Interferometry School, Grenoble,](https://reader036.vdocument.in/reader036/viewer/2022062518/56649d3e5503460f94a17549/html5/thumbnails/33.jpg)
Days 96 and 97 (2010)
Visibility Amplitude on SgrA* at 230 GHz, March 2010
Doeleman et al., private communication
Model fits: (solid) Gaussian, 37 uas FWHM; (dotted) Annular ring, 105/48 as diameter – both with 25 as of interstellar scattering
![Page 34: The Development of High-Resolution Imaging in Radio Astronomy Jim Moran Harvard-Smithsonian Center for Astrophysics 7th IRAM Interferometry School, Grenoble,](https://reader036.vdocument.in/reader036/viewer/2022062518/56649d3e5503460f94a17549/html5/thumbnails/34.jpg)
New (sub)mm VLBI Sites
Phase 1: 7 Telescopes (+ IRAM, PdB, LMT, Chile)
Phase 2: 10 Telescopes (+ Spole, SEST, Haystack)
Phase 3: 13 Telescopes (+ NZ, Africa)
![Page 35: The Development of High-Resolution Imaging in Radio Astronomy Jim Moran Harvard-Smithsonian Center for Astrophysics 7th IRAM Interferometry School, Grenoble,](https://reader036.vdocument.in/reader036/viewer/2022062518/56649d3e5503460f94a17549/html5/thumbnails/35.jpg)
Progression to an Image
Doeleman et al., “The Event Horizon Telescope,” Astro2010: The Astronomy and Astrophysics Decadal Survey, Science White Papers, no. 68
GR Model 7 Stations 13 Stations