1 / 28

Microstructure in the Ionized ISM from radio scattering observations.

Microstructure in the Ionized ISM from radio scattering observations. Barney Rickett UC San Diego O’Dell Symposium Lake Geneva WI April 2007. Radio Probing of the ionized ISM.

velma
Download Presentation

Microstructure in the Ionized ISM from radio scattering observations.

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. Microstructure in the Ionized ISM from radio scattering observations. Barney Rickett UC San Diego O’Dell Symposium Lake Geneva WI April 2007

  2. Radio Probing of the ionized ISM At radio frequencies the ionized ISM is dispersive - in that its refractive index varies strongly with frequency µ = 1 - fp2/f2. Since fp2 ne we have a probe of the electron density in the ISM The simplest observable is from measuring the dispersion in pulse arrival times from pulsars at different frequencies. This gives the column density of ne - called the dispersion measure: DM = ∫0L ne dl accurately known over 1000 lines of sight toward each pulsar. Other observables include: Rotation Measure RM = ∫0L neB.dl from Faraday Rotation Emission Measure EM = ∫0L ne2 dl from Ha emission and radio free-free absorption Scattering Measure SM = ∫0L Cn2 dl from radio scintillation and scattering - which is the topic of my talk. Cn2 gives the strength of variations in ne

  3. P ( ) k 3ne k-(3.67) k k outer inner log(k) The radio scintillation of pulsars and AGNs probe the fine structure in the interstellar electron density ne(s,z) versus transverse position s and distance z from the Earth. From observations one can estimate the structure function of density versus a transverse separation s Dne(s) = ∫0L<[ne(s,z) - ne(s+s,z)]2 > dz (for a pulsar at distance L) This is related to the Power spectrum of density versus transverse wavenumber k Dne(s) = ∫0L∫ P3ne(k,kz=0) {1-eik.s}d2kdz Electron density and its spectrum For scales far from the inner and outer scales in an isotropic Kolmogorov spectrum P3ne(k) = Cn2k-(2+a) ; kinner > k > kouter{kinner = 1/linner ; kouter =1/louter} Dne(s) sawitha=5/3 linner < s < louter Dne2 = ∫P3ne(k) d3k Cn2 louter(2/3) ~ <ne>2? Tiny for ISS is 100 km The scales probed by ISS are 100km to 30 AU

  4. What can we learn from ISS? Interstellar Scintillation (ISS) can probe the form of the spectrum => inner & outer scales, anisotropy for one line of sight And explore its rms amplitude versus Galactic coords => Cn2 or more accurately Scattering Measure SM = ∫Cn2 dz versus latitude, longitude and scattered path length Scattering Measure should be closely related to the Emission Measure The first order description is an isotropic Kolmogorov spectrum that is pervasive throughout the Galaxy - with some volume filling factor. Pulsar dispersion measures combined with independent distance estimates suggest that the local volume average <ne> ~ 0.03 cm-3 Estimates of the volume filling factor are 0.1-0.25 and so imply local electron densities in the ionized regions 0.1-0.3cm-3

  5. Note the relatively minor difference between the Kolmogorov spectrum k-11/3 and k-4, which corresponds to a random superposition of abrupt density steps such as due to shocks, discontinuities etc. k-4 k-11/3 Local Density Spectrum(Armstrong Rickett & Spangler, ApJ 1995) 1 pc 1 AU 100km

  6. ensity Spectrum(Armstrong Rickett & Spangler, ApJ 1995)

  7. Ramachandran et al 2006: Slope = 1.66 Variation of Dispersion Measure psr B1937+21 The solid line gives the best fit line for the data in the time interval of 5 days to 2000 days. The derived values of the the power law index a =1.66 ±0.04 Remarkably close to Kolmogorov value 5/3 When extrapolated down to the diffractive scale ~0.0001 AU, this is within a factor 2 of independent ISS observations. Giving evidence of a Kolmogorov spectrum over range 10-4 to 102 AU (with a possible inner scale at 10-3 AU)

  8. Waves from a point source are scattered by density inhomogeneities and appear to arrive from a range of angles This can be seen as a scattered brightness function of whose width is the diffractive scattering angle qd, which is a very small angle. There is an extra travel time for each scattered path. Hence an ideal impulse is broadened into a one-sided pulse with decay like an exponential of width td td = (Lscatt-L)/c = (z1q12 + z2qd2 )/2c = Lqd2 z2/(2z1c) which is called the diffractive pulse broadening time Angular Broadening and Temporal Broadening z1 q1 qd is related to the diffractive scale sd which is the transverse separation for an observer over which there is an rms difference of 1 radian in propagation phase. qd = l/(2πsd) <<1radian z2 qd Formally sd is defined via the structure function of the propagation phase f(r)at transverse coordinate r , which is: Df(r) = < [f(r+r) - f(r)]2 >  Dne(s) For isotropic Kolmogorov model Df(r) => (r/sd)awith a=5/3

  9. Pulse Broadening If the receiver bandwidth B Scattering time td If td > 1/B the pulse is broadened by scattering, its width tdl4.4 If td < 1/B the pulse does not appear to be broadened, but its amplitude becomes variable. We call this scintillation, which can also be studied as a probe of the ISM. Then the bandwidth over which the scintillations are correlated provides an another way to estimate td .

  10. Bhat et al. ApJ, 2004 (scaled to 1 GHz) Pulse broadening vs DM Observations of many pulsars show that td ( scaled to 1 GHz) depends strongly on DM. For a Kolmogorov spectrum over a path length L we predict: td(msec) ~ Lkpc SM1.2nGHz-4.4 if the inhomogeneities are uniformly distributed, SM  DM  L And we predict tdDM2.2 The observed td increases much more steeply with DM. The large DM pulsars are at low Galactic latitudes and mostly toward the inner Galactic disk. Since DM measures the mean ne, we conclude that the ratioDne/ne increases toward the inner Galactic disk by a factor of more than 20

  11. So there is an enormous increase in “plasma turbulence” toward the inner Galaxy. Dne/ne can be greater than one if the filling factor is less than one - that is the turbulence is sparsely distributed (or “intermittent”) DM dependence 2 Cordes et al 1991 describe this as electron “clouds” which fill a fraction f of the volume and are turbulent internally (with a Kolmogorov spectrum). They show that the turbulence parameter F ~ (Dne2/<ne>2 ) /(f louter2/3) increases by a factor 550 in the inner Galactic plane. louter is the outer scale which could be no larger than the size of the turbulent regions. The 550 factor increase in F implies a 23 times increase in the fractional fluctuations of density or a 550 times decrease in filling factor or a 104 decrease in outer scale… or more reasonably some combination of these. In recent papers Boldyrev and Gwinn characterize this higher turbulence by proposing non-Gaussian statistics for the interstellar electron density - specifically a Levy flight distribution, in which its probability distribution has a power law tail ne-b. They find b ~ 0.7, which is a rather extreme distribution for which both the mean and variance diverge. Boldyrev and Konigl have further shown that such a distribution can arise from lines of sight passing through randomly distributed thin shells of enhanced electron density (~100 cm-2) and/or enhanced turbulence.

  12. Inner Scale Estimates Spangler & Gwinn ApJ 1990: Angular broadening measurements of strongly scattered extra-galactic sources. They measured the precise shape of the visibility function Found inner scales ~ 100 km Suggested the ion inertial scale as the scale where turbulence is dissipated (cutting off the density spectrum) lioninertial = Alfven speed/(ion Larmor frequency) = (ne cm-3 )-0.5 230 km => ne ~ 5 cm-3 If ISS occurs where ne ~ 0.2 cm-3 we expect inner scale ~ 500 km The shape of the far-out tail of scatter-broadened pulses provides another diagnostic. For observations of pulsar J1644-45 at Parkes I found linner ≥ ~100 km in agreement with Spangler and Gwinn.

  13. Anisotropy • There has been increasing evidence in recent years that the scattering plasma often shows evidence for anisotropy. • Scattered images can appear elongated Axial ratio A~1.2-2 • Rapid ISS of quasars (IHV) appears to be “oscillatory” B0405-385 => A>4 • Quasar J1819+38 has annual changes in its ISS timescale with a 6-mo and 12-mo periods that require anisotropy A~6 (maybe source influence) • Scintillation arcs are prominent A ~2-5 (not modeled quantitatively) • Correlated ISS of the two pulsars in J0737-3039 requires anisotropy A>4

  14. Intra-Hour Variables ISS is normally observed in pulsars, but it can also be seen as a low level variation of 1-10% in some compact extra-galactic AGN jet sources. These are called intra-day variables since their flux density varies on time-scales of a day. There are about 5 AGNs which vary on times of I hr or less at 5 GHz and are called intra-hour variables. It has been established that these are due to scintillation in the local ISM - within 1-30 pc of the Earth.

  15. ISS of PKS B0405-385 observed with ATCA ( Rickett, Kedziora-Chudczer & Jauncey ApJ 2002) 4.8 GHz Source Diameter / Screen Distance trade-off 8.6 GHz: modulation index 0.08 < mc < 0.37 Tb constant and time scale 0.31hr <tc< 0.51hr mc constant Conclude that scattering region is between 3-30 pc from Earth. This suggests a relationship with local interstellar clouds. tc constant

  16. ISS Time scale for quasar B1257-326 from Bignall et al.Transverse velocities for local clouds from Linsky and Redfield Aura and Gem clouds match the ISS data Linsky has proposed that collisions between clouds may drive the plasma turbulence responsible for the ISS. Their relative velocities are ~10 km/s. (see Linsky & Redfield poster)

  17. PSR B1133+16 at Arecibo(Stinebring et al.) Scintillation Arcs Pulsar intensity varies with time and frequency due to interstellar Scintillation (ISS). When this “image” is Fourier analyzed the arcs are revealed when the Fourier amplitude is plotted on a log scale. The arcs follow a parabolic curve that relates delay and Doppler shift of the scattered waves. These imply discrete scattering “screens” (here 3 screens) along the line of sight - occupying a small fraction of the path. Primary Dynamic Spectrum “Secondary Spectrum” (S2) with three scintillation arcs

  18. The Puzzle of the “Arc-lets” Hill, Stinebring et al. (2005) showed this example of the arcs observed for pulsar B0834+06. In addition to the main forward arc (following the dotted curve) there are “reverse arclets”. Those labelled a-d are particularly striking. They followed these over 25 days and found that they moved in the secondary spectrum plots, and that the movement was due to the known pulsar proper motion and was consistent with scattering from isolated structures that were stationary in the ISM and survived for at least 20 days.

  19. The Puzzle of the “Arc-lets” 2 Predicted for plasma refraction The right plot shows how the ft values vary with observing frequency. Open circles at 334 MHz and filled circles at 321 MHz. Remarkably this shows that the spatial location of the scatterers is independent of frequency. They DO NOT show the expected shift due to the dispersive nature of plasma refraction. The left plot shows the angular position of the structures (in mas) responsible for each reverse arclet, mapped from the Doppler frequency ft . The lines have the slope expected for the known pulsar proper motion.

  20. The Puzzle of the “Arc-lets” 3 My first thought on seeing these arclets was that they are due to a multipath condition in which four extra ray-paths through the irregular plasma exist at angular offsets further from the unscattered path than the angular width of the scattering disk. Such a multi-path could exist if there are large scale gradients and curvature across the scattering disk, with an amplitude higher than expected from a Kolmogorov medium ( in which the rms phase gradient on a scale s decreases as s-1/6). BUT this idea is entirely incompatible with Dan’s result that the reverse arclets come from a fixed position in space - one that does not scale with frequency as do the stationary phase points that would govern multiple ray-paths in a plasma. Instead there must be isolated “tiny” structures that scatter or refract the waves through angles of the order of 10 mas at 330 MHz. But they must subtend an angle several times smaller than 10 mas or their signatures would overlap in the secondary spectrum.

  21. PSR B0834+06 - VLBI GB-Arecibo at 314 MHz Phase Amplitude Differential Delay Differential Doppler frequency We detect many “reverse arclets” one at a delay of 1 msec, which requires scattering ~40 mas from the pulsar direction. The observations confirm our basic arc interpretation since the fringe phase changes sign with the Doppler frequency.

  22. 0 Dec (mas) -20 -40 0 20 RA (mas) Interferometry of the ArcsW.F. Brisken , J.-P. Macquart, A.T. Deller, C. West, B.J. Rickett, W.A. Coles & S.J. Tingay We can use the interferometer phase to determine the discrete scattering positions that constitute a scattered image. The result is a remarkably elongated image (axial ratio > 10:1) AND an offset feature separated by ~5 times the half-width (7 mas) of the main feature. What is the corresponding spatial structure of the electron density ? (On this scale the electrons still follow the protons). Two possibilities: The image maps the location of extremely fine structure in the electron density, which scatters the radiation by angles larger than their angular extent (otherwise they could not be detected). => Filaments parallel to image axis OR Fine filamentary structure transverse to the image axis, causing the elongated image, but which has strong modulations that cause the offset feature. This is a preliminary estimate of the scattered image.

  23. Cartoons of Scattering Geometry The structures must be narrower than 0.5 AU and have internal fine structure down to scale of 2x106 m What confines it to 0.5 AU? B-field? But this would favor scattering transverse to the elongated image B- field transverse to elongated image, creating fine density structure which scatters perpendicular to B-Field (causing the elongated image). Strong modulation in turbulence transverse to B-field.

  24. Cygnus and Cirrus Nebula NGC 6992

  25. The Puzzle of the “Arc-lets” 4 20 mas at 300pc => 6 AU from the pulsar line of sight; but the ionized “cloud” is 10 times smaller -- dimension a ≤ 0.5 AU say. As for ESEs there are two possibilities pseudo-lens or scatterers: Consider a spherical lens of radius a and electron density nea ( see Hill et al. 2005) It refracts by an angle qr ~ (l/2p) f Roughly f ~ f/a ~ rel nea which is independent of scale a Substituting 10 mas for qr gives nea ~ 2pqr re-1l-2 ~ 200 cm-3 ! If cloud is elongated along the propagation path by an axial ratio R => nea ~ 200/R cm-3 Even if R~10 nea ~ 20 cm-3 is uncomfortably high for a structure only ~0.5 AU Now consider an elliptical scatterer of width a and length Ra with mean electron density nea Let it have a well developed Kolmogorov turbulence interior to a such that dnea ~ nea with an outer scale ~ a. Its scattering angle will be qsc ~ l2.2 (re nea)1.2 a0.2 R0.6 Hence nea ~ 2pqsc re-1l-2 {l/(2pqsca)}1/6 R-.5 With 20 mas for qsc and a = 0.5 AU nea ~ 40 /R0.5 cm-3 Even for R=10 nea ~ 12 cm-3 still uncomfortably high for a structure only ~0.5 AU.

  26. H-alpha map Composite Ha map compiled by Finkbeiner from WHAM, VTSS & SHASSA.

  27. Simulations of Turbulence Kritsuk, Norman et al, 2002 Box is 5 pc 2563 grid points (400AU steps) Starts at 106 K => cools and fragments into regions near 300 and 2x104 K The color coding is log particle density: Dense blobs at the intersections of the filaments, >60 cm-3 , are light blue; Stable cold phase, 6-60 cm-3 , is blue; Unstable density regime, 1.2-6 cm-3 is yellow to brown; Low-density gas, including the stable warm phase <1.2 cm-3, is a transparent red

  28. Summary Kolmogorov spectrum for the interstellar electron density is only a first approximation for the various ISS phenomena. This suggests Plasma Turbulence but does not require it. The medium is veryclumpy => denser and more turbulent regions becoming more common toward the inner Galactic Plane Inner scales consistent with the ion inertial scale are in the range 100-500km Isolated regions may often be anisotropic with axial ratios A>2. Presumably this implies that the magnetic field controls the plasma (ie b < 1) Reverse arclets imply discreteanisotropic structures on sub-AU scales, that have very high electron densities 10-200 cm-3. Suggested geometry is magnetically-controlled filaments of plasma, where the scattering is predominantly perpendicular to B-field and is highly modulated in that direction

More Related