1 / 40

Basics of mm interferometry

Basics of mm interferometry. Sébastien Muller Nordic ARC Onsala Space Observatory, Sweden. Turku Summer School – June 2009. Interests of mm radioastronomy. -> Cold Universe Giant Molecular Clouds -> COLD and DENSE phase Site of the STAR FORMATION -> Continuum emission of cold dust

Download Presentation

Basics of mm interferometry

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Basics of mm interferometry Sébastien Muller Nordic ARC Onsala Space Observatory, Sweden Turku Summer School – June 2009

  2. Interests of mm radioastronomy -> Cold Universe Giant Molecular Clouds -> COLD and DENSE phase Site of the STAR FORMATION -> Continuum emission of cold dust -> Molecular transitions - Diagnostics of the gas properties (temperature, density) - Kinematics (outflows, rotation)

  3. Where the atmosphere is relatively transparent Interests of CO Molecular gas  H2 But H2 symmetric ->electric dipolar momentum is 0 Most abundant molecule after H2 is CO [CO/H2] ~ 10-4 First rotational transitions of CO in the mm CO(1-0) @115 GHz CO(2-1) @230 GHz CO(3-2) @345 GHz E J=1,2,3 = 6, 17, 33 K Easily excited CO is difficult to destroy high ionization potential (14eV) and dissociation energy (11 eV)

  4. Handy formulae - HI line emission: N(HI) (cm-2) = 1.82 1018TBdv (K km/s) - Molecular line emission: N(H2) (cm-2) = X 1020TCOdv (K km/s) X = 0.5-3 Or use optically thin lines (13CO, C18O) - Visual extinction: N(HI)+2N (H2) (cm-2) = 2 1021AV (mag)

  5. Needs of angular resolution Resolution  /D (theory of diffraction)

  6. Would need very large single-dish antennas BUT - Surface accuracy(few 10s of microns !) -> technically difficult and expensive ! - Small field of view (1 pixel) - Pointing accuracy (fraction of the beam) Let’s fill in a large collecting area with small antennas And combine the signal they receive -> Interferometry (Aperture synthesis)

  7. Mm antennas need Good surface accuracy D APEX 12m <20 microns IRAM-30m 30m 55 microns (GBT 100m 300 microns) PdBI 15m <50 microns SMA 6m <20 microns ALMA 12m <25 microns

  8. Holography measurement

  9. source uv plane geometrical time delay antenna baseline Baseline, uv plane and spatial frequency • uv positions are the projection of the baseline vectors • Bij as seen from the source. • The distances (u2 + v2) are refered to as spatial frequencies • Interferometers can access the spatial frequencies • ONLY between Bmin and Bmax, the shortest and longest • projected baselines respectively.

  10. Interferometers measure VISIBILITIES V V(u,v) =  P(x,y) I(x,y) exp –i2(ux+vy) dxdy = FT{ P I } But astronomers want the SKY BRIGHTNESS DISTRIBUTION of the source : I(x,y) P(x,y) is the PRIMARY BEAM of the antennas - P has a finite support, so the field of view is limited - distorded source informations - P is in principle known ie. antenna characteristic

  11. I(x,y) P(x,y) =  V(u,v) exp i2(ux+vy) dudv Well, looks easy … BUT ! Interferometers have an irregular and limited uv sampling : - high spatial frequency (limit the resolution) - low spatial frequency (problem with wide field imaging) Incomplete sampling, non respect of the Nyquist’s criterion = LOSS of informations ! The direct deconvolution is not possible Need to use some smart algorithms (e.g. CLEAN)

  12. I x x0 Phase Amplitude S Slope = -2x0 u u Let’s take an easy example: 1D P = 1 I(x) = Dirac function: S(x-x0) S = constant V(u) = FT(I) = Sexp(-i2ux0) -> this is a complex value

  13. Illustration : dirty beam, dirty image and deconvolved (clean) image resulting in some interferometric observations of a source model

  14. Atmosphere « The atmosphere is the worst part of an astronomical instrument » - emits thermally, thus add noise - absorbs incoming radiation - is turbulent ! (seeing) Changes in refractive index introduce phase delay Phase noise -> DECORRELATION(more on long baselines) exp(-2/2) - Main enemy is water vapor (Scale height ~2 km)

  15. O2 H2O

  16. Calibration Vobs = G Vtrue + N Vobs = observed visibilities Vtrue = true visibilies = FT(sky) G = (complex) gains usually can be decomposed into antenna-based terms: G = Gij= Gi x Gj* N = noise After calibration: Vcorr = G’ –1 Vobs

  17. Calibration - Frequency-dependent response of the system Bandpass calibration -> Bright continuum source - Time-dependent response of the system Gain (phase and amplitude) -> Nearby quasars - Absolute flux scale calibration -> Flux calibrator

  18. Bandpass calibration

  19. Phase calibration

  20. Amplitude calibration

  21. From SMA Observer Center Tools http://sma1.sma.hawaii.edu/

  22. From SMA Observer Center Tools http://sma1.sma.hawaii.edu/

  23. From SMA Observer Center Tools http://sma1.sma.hawaii.edu/

  24. Quasars usually variable ! -> need reliable flux calibrator From SMA Observer Center Tools http://sma1.sma.hawaii.edu/

  25. Preparing a proposal 0) Search in Archives SMA: http://www.cfa.harvard.edu/cgi-bin/sma/smaarch.pl PdBI: http://vizier.cfa.harvard.edu/viz-bin/VizieR?-source=B/iram ALMA … 1) Science justifications  -> Model(s) to interpret the data 2) Technical feasibility: - Array configuration(s) (angular resolution, goals) - Sensitivity use Time Estimator ! Point source sensitivity Brightness sensitivity (extended sources)

  26. Array configuration Compact Detection Mapping of extended regions Intermediate Mapping Extended High angular resolution mapping Astrometry Very-extended Size measurements Astrometry

  27. PdBI 1 Jy = 10-26 W m-2 Hz-1

  28. For extended source: Take into account the synthesized beam -> brightness sensitivity T (K) = 2ln2c2/k2 x Flux density/majmin Use Time Estimator !

  29. Short spacings V(u,v) =  P(x,y) I(x,y) exp –i2(ux+vy) dxdy V(0,0) =  P(x,y) I(x,y) dxdy (Forget P), this is the total flux of the source And it is NOT measured by an interferometer ! -> Problem for extended sources !!! -> Try to fill in the short spacings

  30. Courtesy J. Pety

  31. Courtesy J. Pety

  32. Advantages of interferometers - High angular resolution - Large collecting area - Flatter baselines - Astrometry - Can filter out extended emission - Large field of view with independent pixels - Flexible angular resolution (different configuration)

  33. Disadvantages of interferometers - Require stable atmosphere - High altitude and ~flat site (usually difficult to access) - Lots of receivers to do - Complex correlator - Can filter out extended emission - Need time and different configuration to fill in the uv-plane

  34. Mm interferometry: summary - Essential to study the Cold Universe (SF) - Astrophysics: temperature, density, kinematics … - Robust technique High angular resolution High spectral/velocity resolution

  35. Let’s define - Sampling function S(u,v) = 1 at (u,v) points where visibilities are measured = 0 elsewhere - Weighting function W(u,v) = weights of the visibilities (arbitrary) We get : Iobs(x,y) =  V(u,v) S(u,v) W(u,v) exp i2(ux+vy) dudv

  36. Due to the Fourier Transform properties : FT { A B } = FT { A } ** FT { B } Can be rewritten as : where Iobs(x,y) =  V(u,v) S(u,v) W(u,v) exp i2(ux+vy) dudv Iobs(x,y) = P(x,y) I(x,y) ** D(x,y) D(x,y) =  S(u,v) W(u,v) exp i2(ux+vy) dudv = FT { S W }

  37. Iobs(x,y) = P(x,y) I(x,y) ** D(x,y) If Isou = (x,y) = Point source then Iobs(x,y) = D(x,y) That is : D is the image of a point source as seen by the interferometer. ~ Point Spread Function

  38. D(x,y) = FT { S W } D(x,y) is called DIRTY BEAM This dirty beam depends on : - the uv sampling (uv coverage) S - the weighting function W Note that :  D(x,y) dxdy = 0 because S(0,0) = 0 And that : D(0,0) > 0 because SW > 0 The dirty beam presents a positive peak at the center, surrounded by a complex pattern of positive and negative sidelobes, which depends on the uv coverage and the weighting function.

  39. Iobs(x,y) = P(x,y) I(x,y) ** D(x,y) Iobs(x,y) is called DIRTY IMAGE We want Iobs(x,y) I(x,y) This includes the two key issues for imaging : - Fourier Transform (to obtain Iobs from V and S) - Deconvolution (to obtain I from Iobs)

More Related