it is well known that radio waves propagating in the interstellar medium ( ism ) are scattered by the irregularities in the galactic electron density .
the scattering in turn gives rise to a number of observable phenomena . among others ,
these include angular broadening and intensity fluctuations ( both in time and frequency ) of compact radio sources . while a nuisance in many radio astronomical observations , these phenomena can be used to investigate the nature of the irregularities in the interstellar plasma density .
these density irregularities are in turn believed to follow the fluctuations in the interstellar kinetic and magnetic energies .
ideally one would like to invert observations of radio scintillation and scattering to determine the statistics of the plasma density . as noted by narayan ( 1992 )
this inverse problem is not well posed and one must rely on modeling methods .
a complete prediction of scintillation observables requires an _ a priori _ knowledge of both the form of the spatial power spectrum of the electron density fluctuations and the distribution of the scattering material along the line of sight . for a given profile of the distribution of the scattering material ,
one may compare observations and predictions and so constrain the functional form of the spectrum .
the power spectrum provides useful insight into the physics of the plasma irregularities .
hence a knowledge of the form of the density spectrum becomes central in both predicting scintillation phenomena and understanding the physics of the interstellar plasma . in this paper , we revisit the investigation of the form of the density spectrum .    a commonly used model for the density spectrum has been based on a power - law model with a large range between `` inner '' and `` outer '' scales ( e.g.  rickett 1977 ) : @xmath1 \ , \mbox{. } \label{eq : extpowerlaw } \ ] ] here @xmath2 is the magnitude of the three - dimensional wavenumber @xmath3 .
@xmath4 denotes the strength of fluctuations ( with a weak dependence on distance @xmath5 ) .
@xmath6 is the spectral exponent , and @xmath7 and @xmath8 are the inner and outer scales respectively . for @xmath9 , we obtain the _ simple _ power - law model : @xmath10 .
armstrong , rickett , & spangler ( 1995 ) have constructed an empirical density spectrum by combining radio scintillation observations in the local ism ( @xmath11 1  kpc ) with measurements of the differential faraday rotation angles and large - scale electron density gradients .
they have shown the power spectrum to be consistent with a simple kolmogorov power - law model ( @xmath12 ) over an astronomical 10 orders of magnitude in wavenumber scale ( @xmath13m@xmath14m@xmath15 ) .
the kolmogorov spectrum in density suggests a turbulent cascade in the magnetic and kinetic energies . this has lead to several theoretical investigations of the generation and maintenance of hydromagnetic turbulence in the ism ( e.g.  pouquet , 1978 ; higdon 1984 & 1986 ; biskamp , 1993 : sridhar & goldreich 1994 ; goldreich & sridhar 1995 and 1997 ) . however , the armstrong et al .  study combined observations from many lines of sight and the scatter among them leaves a substantial uncertainty in the exponent @xmath6 .
a list of symbols is given in table 3 .    in spite of the positive evidence for the simple kolmogorov spectrum , substantial observational inconsistencies remain .
for instance , long - term refractive intensity scintillations of some pulsars have modulation indices as much as a factor of 2 larger than predicted by the simple kolmogorov model ( cf.gupta , rickett , & coles 1993 ) .
other discrepancies are revealed in the diffractive dynamic spectra of pulsars . on some occasions , periodic fringes
are observed , which are not predicted by the simple kolmogorov model ( cf .  roberts & ables 1982 ; cordes & wolszczan 1986 ; rickett , lyne , & gupta 1997 ; gupta , bhat & rao , 1999 ) ; in addition , sloping bands in the dynamic spectra often persist longer than predicted by the model ( cf.gupta , rickett , & lyne 1994 ; bhat , rao & gupta , 1999b bhat , gupta & rao , 1999c ) .
further for some pulsars , the decorrelation bandwidth has larger amplitude variations than predicted for the kolmogorov spectrum ( bhat , gupta & rao , 1999c ) .
such anomalies suggest the presence of large refractive structures giving rise to the focusing and defocusing of the scattered ray bundles .
the interference of the ray bundles can also explain the occasional fringes observed in the dynamic spectra of some pulsars . given the relatively frequent occurrence of such events
, one can ask whether they can be considered as mere occasional anomalies or should be considered as a widespread phenomenon intrinsic to the spectrum on a grand scale .
the observational inconsistencies suggest the need for an enhancement in the power on the large `` refractive '' ( 10@xmath16  m to 10@xmath17  m ) spatial scales relative to the power on the small `` diffractive '' ( 10@xmath18  m to 10@xmath19  m ) scales .
there are several means by which this ratio may be enhanced .
one is to include the inner scale cut - off in the density spectrum , which reduces the power at small scales ( coles et al .
1987 ) ; these authors proposed inner scales of @xmath20 m. , though this does not correspond to any obvious physical scale . physically
, the inner scale corresponds to the scale at which the turbulent cascade dissipates and becomes a source of heating for the plasma ( spangler , 1991 ) .
the value of the inner scale is largely unknown . using different methods , spangler & gwinn ( 1990 ) , kaspi & stinebring ( 1992 ) , and gupta et al .
( 1993 ) have reported values for the inner scale ranging from @xmath21 to @xmath22 meters . in proposing the smaller values ,
spangler & gwinn ( 1990 ) argued that the inner scale is the larger of the ion inertial length , @xmath23 ( where @xmath24 is the alfvn speed , and @xmath25 is the ion cyclotron frequency ) , and the ion larmor radius , @xmath26 ( where @xmath27 is the ion thermal speed ) ; they obtained parameters for the warm ionized medium in reasonable agreement with observations . in a recent discussion , minter & spangler ( 1997 )
have suggested ion - neutral collisional damping and wave - packet steepening as possible dissipation mechanisms for the turbulence in the diffuse ionized gas , which would make the mean - free path for ion - neutral collisions a possible value for the inner scale .
however , this is thought to be larger than the maximum values proposed to explain the observations , making it a less convincing dissipation mechanism .
observationally , the kolmogorov model with a large inner scale predicts refractive modulation indices consistent with pulsar measurements ( gupta et al .
it has also been proposed to explain the occasional periodic fringes , with an inner scale on the order of the fresnel scale ( cordes , pidwerbetsky , & lovelace 1986 ; goodman et al .
. however , rickett et al .  (
1997 ) reported a fringe event for the pulsar b0834 + 06 that could not be explained as the effect of a large inner - scale spectrum .
the event requires similar conditions to those needed to explain the extreme scattering events ( fiedler et al .
1987 ) .
another way to enhance the ratio of the power between the refractive and diffractive scales in the spectrum is to steepen the spectrum  with spectral exponents @xmath28 ( blandford & narayan 1985 ; goodman & narayan 1985 ; romani , narayan , & blandford 1986 ) .
while power - law spectra with @xmath29 11/3 have a turbulence connotation , spectra with @xmath30 might involve some forms of turbulence , but are also consistent with a distribution of non - turbulent structures with a range of spatial scales .
such steep spectra with @xmath31 predict refractive modulation indices close to unity ( goodman & narayan 1985 ) , which is substantially larger than the range 30 to 40 % observed from the nearby pulsars . on this basis rickett & lyne ( 1990 ) and armstrong et al .
( 1995 ) have rejected spectra steeper than 4 for the interstellar plasma .
however , the special case of @xmath32 has been given little attention .
we can conceive of a power - law model with spectral exponent @xmath32 given by : @xmath33 hereafter , we refer to this model as the `` @xmath32 model . ''
blandford & narayan ( 1985 ) briefly discussed this special case without including a cut - off at low wavenumbers .
it is interesting to note that even though the kolmogorov spectral exponent @xmath12 is very close to 4 , the @xmath32 model has a very different physical implication
. its physical origin has been rarely discussed , and it has not been formally compared with observations .
physically , this spectrum suggests the random distribution in location and orientation of discrete discontinuous objects across the line of sight .
an `` outer '' scale , @xmath34 , is included to account for the typical size of such objects .
the `` inner '' scale here would correspond to the scale of the sharpness of a typical discontinuity .
we assume the `` inner '' scale to be smaller than the diffractive scale of scintillations , and hence has no significant effect on the scintillations .
the @xmath32 model could characterize stellar wind boundaries , supernova shock fronts , sharp boundaries of hii regions at the strmgren radius , or any plasma `` cloud '' with sufficiently sharp boundaries ( transition regions shorter than the diffractive scintillation scale ) which may cross the line of sight .
note that , though turbulence is not necessarily implied by the model , strong turbulence which has steepened to form shocks would also be described by the @xmath32 model . in these models a single discrete object crossing the line of sight could also explain the `` extreme scattering events '' observed in the flux density variations of extra - galactic sources ( e.g.  fiedler et al .
we note that the analysis of such events have been primarily based on geometrical optics involving a single `` cloud '' ( goodman et al .
1987 ; cordes et al .  1986 ; cordes & wolszczan 1986 ; roberts & ables 1982 ; ewing et al .  1970 ) .
our analysis of the @xmath32 model includes the `` wave optics '' effects when the line of sight passes through very many such clouds .
the ism is assumed to consist of a random assembly of discrete structures with abrupt density steps which may be independent of each other .
if the @xmath32 model were compatible with all of the scintillation observations , it would remove the implication of interstellar plasma turbulence with an inertial range spanning as much as 10 orders of magnitude in scale , which has become the canonical model for iss phenomena .
it would , however , still be consistent with turbulence that has steepened into shocks .
we start with a simple derivation of the @xmath32 model in section [ sec : deriv ] . in section
[ sec : structfn ] , we give equations for the wave structure function for the @xmath32 model and compare predictions in section [ sec : decbw ] , the variation of the diffractive decorrelation bandwidth with frequency is used to test the @xmath32 model against the simple kolmogorov model and the kolmogorov model with an inner - scale ( hereafter , the `` inner - scale '' model ) . in section
[ sec : mr ] , the observed variation of the refractive scintillation index with the normalized diffractive decorrelation bandwidth is compared with theoretical predictions for the simple kolmogorov , inner - scale , and @xmath32 models .
we use theoretical results from a previous paper ( lambert & rickett 1999 ; hereafter , lr99 ) , in which we developed the theory of diffractive scintillations in a medium modeled by these spectra , and we add details of the theory for refractive scintillation in appendix a. in appendix b , we discuss the use of the square of the second moment to approximate the intensity correlation function for diffractive scintillation with @xmath32 . in section [ sec : conclusion ] we give a discussion and our conclusions .
in this section , we give a simple derivation of the @xmath32 model .
consider first a random distribution in space of identical plasma structures ( blobs ) .
the resulting electron density may be described as the convolution of a three - dimensional poisson point process with the density profile of one individual blob .
this is the spatial analog of shot noise , in which each charge carrier contributes the same temporal profile of current arriving at randomly distributed times . in the wavenumber domain , the power spectrum of the electron density
is then given by the spectrum of the poisson point process multiplied by the squared magnitude of the fourier transform of one blob .
the spectrum of the poisson point process is a constant equal to the number density , @xmath35 , of the blobs in space ( papoulis 1991 ) .
thus the power spectrum of the medium follows the same shape as that of an individual blob .
a similar description applies to water droplets in a fog and many naturally occurring media , as discussed by ratcliffe ( 1956 ) .
if the blobs are asymmetrical and randomly oriented , one must also average the spectrum over the possible angles of orientation .
the basic feature of structures with discontinuous boundaries is that their power spectrum has an asymptotic behavior at large wavenumbers that varies as ( wavenumber)@xmath36 . this is the spatial analog of the idea that the temporal spectrum of any pulse shape with an abrupt rise or fall has a high - frequency asymptote as ( frequency)@xmath37 .
the prime example is a rectangular pulse , for which the spectrum is a sinc - squared function , whose high - frequency envelope follows this law .
the simplest spatial example is spherical blobs of radius @xmath38 with uniform density @xmath39 inside and zero outside ; clearly , an added uniform background density does not affect the power spectrum .
the three - dimensional spatial fourier transform of an _ isotropic function @xmath40 of radial distance @xmath41 is given by ( tatarskii 1961 ) : @xmath42 where @xmath2 is the magnitude of the wavenumber .
for any function @xmath40 that falls to zero beyond a radius @xmath38 , the integral will vary as @xmath43 in the high - wavenumber limit @xmath44 , making @xmath45 . for the spherical blobs of size @xmath38
we have : @xmath46 for the squared magnitude of the fourier transform , for very small and large values of @xmath47 : @xmath48 { \displaystyle   q^{-4 } \ ; a^2 f_0 ^ 2 /(8\pi^4 ) \ , \mbox { , }   } & q a > > 1 \ , \mbox{. } \\
\end{array } \right .
\label{eq : magsqspherefourier } \ ] ] the spectrum of the electron density is then @xmath49 .
we can approximate this by equation ( [ eq : beta4 ] ) , which exactly matches equation ( [ eq : magsqspherefourier ] ) at small and large values of @xmath50 , with @xmath51 and @xmath52 .
_    there are evident generalizations to make the model more like a real medium .
the spheres could have a probability distribution for their radii , which would weight the average of the square of equation ( [ eq : spherefourier ] ) .
the spheres could have smoothly varying density inside and an abrupt boundary at radius @xmath38 .
as with the pulse example , the high - wavenumber behavior would not change .
more realistic models would include anisotropic structures with random orientations and a distribution of scales .
for example , consider a circular cylinder as a model for an anisotropic blob , with a particular radius and thickness .
its fourier transform can be computed with respect to its major axes .
this can then be transformed to tilted axes and the result averaged over an isotropic distribution of angles of tilt . evidently the result will also be isotropic ; we have solved this approximately and find a spectrum with the same high - wavenumber _ asymptotic form as equation ( [ eq : beta4 ] ) , where the `` outer scale , '' @xmath53 , is the smaller of the thickness and radius .
this result could be averaged over a probability distribution of radii and thicknesses and a similar conclusion reached , with @xmath53 as a weighted average of the smallest dimension .
planar sheets at random angles are modeled as very flat cylinders and provide a crude representation of a distribution of shocks .
another simplified model would be ellipsoids of different sizes , eccentricities , and orientations , and again we expect a similar result . _    we conclude that the @xmath32 spectrum is an approximate description of the power spectrum for a medium with a random distribution of discontinuous structures , in the limit for wavenumbers greater than the reciprocal of the smallest dimension across the structures .
we note that for wavenumbers at or below about @xmath54 the functional form depends on the detailed shapes , though it will be much less steep than @xmath55 .
thus equation ( [ eq : beta4 ] ) will not be generally applicable at the lowest wavenumbers .
similarly , in practice there is a finite scale over which the density jump occurs ( e.g.  shock thickness ) .
this can be modeled by convolving the idealized discontinuous structures with a suitable narrow , say , gaussian function which thus provides a very high - wavenumber cut - off beyond which the spectrum falls rapidly to zero .
consequently , we consider a spectrum with an inverse fourth power law between inner and outer scales , though these do not have the same connotations as inner and outer scales in a turbulent cascade .
we do not discuss the physical origin of the supposed discontinuous plasma structures beyond the ideas formulated in the introduction , such as supernova shock fronts , stellar wind boundaries or sharp boundaries of hii regions at the strmgren radius .
there needs to be no physical coupling between the individual blob structures , and there is no turbulent cascade implied by the @xmath32 spectrum .
however , we note that shocks developing from the steepening of very strong turbulence could also be described by the model , and so if successful , the model would not rule out such strong turbulence , though in this case the @xmath32 range does not correspond to an inertial cascade as does the 11/3 spectrum .    in the following sections , we give the necessary scintillation theory needed to develop equations for the scintillation observables corresponding to the @xmath32 model , which we then compare with iss observations .
we concentrate on measurements which are sensitive to the form of the density spectrum and are relatively insensitive to the distribution of scattering material along the line of sight .
central to the description of the second moment and the coherence function for intensity is the phase structure function .
the longitudinal gradient of the phase structure function is related to the electron density spectrum through ( cf .
coles et al .   1987 ) : @xmath56 \kappa d\kappa \ , \mbox{. }
\label{eq : long_grad_phase_strucfunc } \ ] ] here @xmath57 is the speed of light , @xmath58 is the classical electron radius , and @xmath59 is the radio frequency . for a plane wave incident on a scattering medium , the line of sight integral of @xmath60 gives the structure function of the geometric optics phase , also called the _ wave _ structure function , @xmath61 .
for the spherical wave geometry applicable to a pulsar at distance @xmath62 , the separation @xmath63 at the observer s plane is projected back to the point of integration , giving @xmath64 .
the integrals are evaluated explicitly in equations ( 13 ) , ( 17 ) , and ( 21 ) of lr99 for a uniform scattering medium with the simple power - law , inner - scale , and @xmath32 spectral models , respectively . in the extended scattering geometry ,
the field coherence length at the observer , @xmath65 , is defined by : @xmath66 . for a screen of thickness @xmath67 at distance @xmath68 from the pulsar ,
the wave structure function at the observer becomes @xmath69 , which has an explicit dependence on the screen location . in lr99
, we found it useful to define a coherence scale @xmath70 for the screen by its phase structure function @xmath71 , which is independent of the screen location . note that in applying equation ( [ eq : long_grad_phase_strucfunc ] ) , we are assuming that the screen thickness , @xmath67 , is much greater than the outer scale @xmath53 and small compared to total distance @xmath62 .    here
, we repeat the equation for the phase structure function of the @xmath32 model in a thin layer ( screen ) . substituting the @xmath32 model for the density spectrum into equation ( [ eq : long_grad_phase_strucfunc ] ) , integrating , and multiplying by the screen thickness , @xmath67
, we obtain : @xmath72 \
, \mbox { , } \label{eq : struc4exact } \ ] ] where @xmath73 is the first order modified bessel function of the second kind , and @xmath74 is the scattering measure . equation ( [ eq : struc4exact ] ) may be approximated by the following logarithmic expression : @xmath75 \ , \mbox { , } \label{eq : struc4 } \ ] ] which matches the full bessel expression at both large and small values of the argument .
the detailed shape , near where the structure function flattens , is governed by the shape of the low - frequency turnover in our model spectrum , equation ( [ eq : beta4 ] ) ; however , as noted above , the turnover shape depends on details of the density profiles in the plasma blobs which are not constrained in our density model .
the logarithmic structure function for the @xmath32 model has far - reaching consequences .
as the outer scale becomes larger , the structure function approaches a square - law .
however , due to the presence of the logarithm , this occurs only slowly ; that is , only at the limit as @xmath63 goes to zero does the @xmath32 model structure function become exactly square - law . a square - law structure function for the medium
is a convenient mathematical model , but as noted in lr99 , in a real medium the structure function must eventually saturate at very large separations .
nevertheless , the square - law structure function has been widely used since it provides a valid approximation for scales smaller than the reciprocal of the wavenumber at which the spectrum is cut - off . the fact that the @xmath32 model structure function deviates from a pure square - law , even for very large values of the outer scale , makes this model interesting and an independent investigation worthwhile .    in figure
[ fig : slope ] , we have plotted the structure function of the @xmath32 model and its effective local logarithmic slope , @xmath76 , versus the transverse spatial lag , @xmath63 , for various values of the outer scale , @xmath53 ( see equation 20 of lr99 ) . to illustrate the shapes , we have chosen @xmath77  km for all curves .
the thin dashed lines correspond to the simple kolmogorov model with logarithmic slope of 5/3 , and the thin dotted lines correspond to the square - law structure function .
the plot shows how slowly the local slope approaches 2 ; even for @xmath78 , the slope is 1.95distinctly less than 2 .
the points in the upper panel correspond to measured local exponent of the wave structure function obtained from vlbi observations , as tabulated by spangler & gwinn ( 1990 ) . in principle , this method of comparison provides a good test for the form of the density spectrum .
these authors presented a similar plot corresponding to the inner - scale model , for which @xmath76 is close to 2.0 for @xmath63 less than the inner scale and close to 5/3 for @xmath63 greater than the inner scale , with a transition occurring over about a decade in @xmath63 .
they derived estimates of the inner scale ranging from 50,000 to 200,000  meters , for the highly scattered sources that they studied . whereas there is some uncertainty in some of the estimates of @xmath79
, there are several examples in which values greater than 5/3 are reliably observed ( e.g.  trotter et al . , 1998 ) .
as can be seen , there is also reasonable agreement with the @xmath32 model . assuming the model to be correct , a separate model - fitting to each observation yields an outer scale in the range @xmath80 to @xmath81  meters .
these values are much smaller than the parsec scales of supernova remnants and interstellar clouds .
however based on this comparison alone , we are unable to discriminate between the inner - scale and @xmath32 models since they both provide equally satisfactory agreement with the observations
. hence further independent tests are needed , and in the following sections , we present comparisons with two other scintillation observables which in the end argue against the @xmath32 model as a universal description of the ionized ism .
the diffractive decorrelation bandwidth , @xmath82 , is perhaps the easiest iss observable to measure .
it is the frequency - difference for a decorrelation to 50% of the correlation function for diffractive intensity scintillations .
these are usually recorded as a dynamic spectrum centered on a given radio frequency , @xmath59 , for a pulsar .
observations of @xmath82 for one pulsar at a wide range of radio frequencies have provided an important test of power - law models for the interstellar density spectrum on that particular line of sight .
when the diffractive scale is far from the inner and outer scales , one expects @xmath83 , which for the simple kolmogorov model is @xmath84 , in both screen and extended medium geometries . in this section
we re - examine the published measurements and compare them with theory for the three different spectral models .
since the scaling laws are all so steep , in figures [ fig : vela ] and [ fig : other ] , we have plotted the observations logarithmically as @xmath85 .
overplotted are lines giving the theoretical predictions for our three spectral models , with the simple kolmogorov model as a line of slope 0.4 . before giving our conclusions
, we must discuss the error bars and various corrections to the data and also the method for deriving the theoretical predictions .
decorrelation bandwidth data , collected together by cordes , weisberg , & boriakoff ( 1985 ; hereafter , cwb ) , for pulsars psr b0833@xmath8645 ( vela ) , psr b0329@xmath8754 , psr b1642@xmath8803 , psr b1749@xmath8628 , and psr b1933@xmath8716 are shown in figures [ fig : vela ] and [ fig : other ] .
we have also included recently measured points by johnston et al .
( 1998 ) for the vela pulsar .
these are @xmath89 and @xmath90 .
the data are shown as solid circles if derived from dynamic spectra and as solid stars if derived from pulse broadening measurements .
the latter become easier to estimate at low frequencies , where the resolution bandwidth required for the former becomes too narrow for an adequate signal - to - noise ratio .
the conversion of an estimate of a pulse broadening time to a decorrelation bandwidth relies on the uncertainty relation @xmath91 .
unfortunately , the `` constant '' @xmath92 takes on different values for different geometries and spectral models . for the @xmath32 model ,
it is also very weakly dependent on frequency . in table 1 of lr99
, we gave numerical values for @xmath92 for the three spectral models with spherical waves in both a screen and an extended medium geometry .
for a given geometry , the smallest @xmath92 is for the kolmogorov spectrum and the largest ( by about 50% ) is for the square - law structure function .
thus in figures [ fig : vela ] and [ fig : other ] we plot , as solid and open stars , the values converted using @xmath92 for the kolmogorov and square - law structure function models , respectively .
we display the predictions for extended medium and note that the screen values are only about 20% less for each spectrum .
thus in comparing data with the theory for the @xmath0 and inner - scale models , the data should lie between the extremes of the open and solid stars , depending on the outer and inner scales , respectively .
note that the value of @xmath92 in equation ( 7 ) of the taylor et al .
( 1993 ) catalog is more than 50% greater than values computed by lr99 , since it concerns the mean pulse delay rather than the 1/e decay time of the pulse .    the second correction to the data concerns the effects of refractive scintillation at the higher radio frequencies , which are closer to the transition to weak scintillation .
gupta et al .  ( 1994 ) gave a heuristic theory of the effect and expressions for the bias to the diffractive @xmath82 due to refractive shifts ; lr99 discussed the effect from the work of codona et al .
( 1986 ; hereafter , ccffh ) . in figures [ fig : vela ] and
[ fig : other ] , we have plotted as open circles @xmath82 corrected for this effect using equation ( d6 ) of gupta et al .
( 1994 ) , who also estimated the variability in estimates of @xmath82 due to the refractive modulation .
we used this to estimate an error bar on the open circles , which is typically larger than the error bars quoted by the original observers , plotted on the solid circles .    in summary ,
the open circles and their error bars provide the best direct estimates of @xmath82 and the indirect estimates lie between the open and solid stars .
the theory of diffractive scintillations has been discussed by many authors ; we will use the results of section 6 in lr99 for a point source ( spherical waves ) either in an extended scattering medium or with a screen at @xmath68 from the pulsar and @xmath93 from the observer ( with @xmath94 ) . following usual practice , we assume the strong scattering limit , in which the frequency decorrelation function for intensity is the squared magnitude of the diffractive component of the two - frequency second moment of the field .
the validity of this approximation is examined in appendix b , where we find that it introduces a small error for the @xmath32 model , which however is negligible compared to the errors in the observations .
the decorrelation bandwidth , @xmath82 , can be defined in terms of a normalized bandwidth , @xmath95 : @xmath96 here @xmath59 is the ( geometric ) mean of the two frequencies , and the parameter @xmath97 determines the strength of scattering ; @xmath57 is the speed of light and @xmath98 is the effective scattering distance . for a uniform scattering medium , @xmath99 , and for a screen geometry ,
@xmath101 is the field coherence scale where the appropriate phase structure function equals unity ; for the uniform scattering medium , it is defined at the observer , and for a screen , it is @xmath70 which is independent of the distances to the pulsar and screen , and leads to @xmath102 as the appropriate strength of scattering parameter .
the computations presented in lr99 include a tabulation of the normalized bandwidth @xmath95 for the various models under discussion , and associated plots of the intensity decorrelation function itself .
for the simple kolmogorov model , a constant value for @xmath95 is obtained regardless of the frequency or distance ( 0.773 and 0.654 for the extended medium and screen geometries , respectively ) .
thus the frequency dependence of @xmath82 reflects the frequency dependence of @xmath65 , giving @xmath83 .
cwb used this to estimate @xmath6 from observations of @xmath82 .
we now compare the observations with theoretical scaling laws for the two other spectral models .
there are , however , some complications .
the quantity @xmath95 is no longer exactly independent of frequency , and @xmath65 depends on both frequency and the outer scale , @xmath103 , or the inner scale , @xmath104 , that parameterize the other two models .
consider the details for the discontinuity spectrum . in figure
[ fig : vela ] , we plot @xmath105 versus @xmath106 as observed for the vela pulsar . in figure
[ fig : vela]a , we have plotted theoretical curves for @xmath32 model with various outer scales .
we use equation ( [ eq : delta_nu_d ] ) for the models and so need the variation of both @xmath65 and @xmath95 with @xmath59 .
for a screen , we determine @xmath70 from equation ( [ eq : struc4 ] ) for the @xmath32 model , and eliminate the scattering measure @xmath107 using the same equation for @xmath108 at a reference frequency of 1 ghz .
this gives a relation between @xmath109 and @xmath106 , with @xmath110 as a parameter .
for an extended scattering medium , we use equation ( 21 ) of lr99 to determine @xmath65 at the observer , and the same method as for the screen to eliminate @xmath107 . in figure [ fig : vd ] , we show the variation of @xmath95 with @xmath111 at a fixed frequency ( computed by lambert 1998 and described in lr99 ) .
this is combined with the @xmath65@xmath106 relation to obtain the relation of @xmath95 to @xmath106 . for the curves in figure [ fig : vela ] , @xmath108 is determined so that , for each value of @xmath53 , the model @xmath112 fits the measured values near 1ghz . for most of the models of interest @xmath113 , and then
the logarithm functions vary much more slowly than the @xmath114 term in equation ( [ eq : struc4 ] ) . in which case there is nearly a linear relation @xmath115 .
furthermore , figure [ fig : vd ] shows that as @xmath65 changes , the parameter @xmath95 changes very slowly .
consequently , for a very wide range of the outer scale , a good approximation is @xmath116 , which reflects the fact that the underlying structure function then approaches a square - law behavior .
consider the plot for the vela pulsar psr  b0833@xmath8645 ( figure [ fig : vela ] ) for which the data are most extensive .
although we have computed curves and data corrections for both the screen and extended scattering geometries , we only show the extended medium results in figure [ fig : vela ] because the two plots are so similar . in the comparison with the screen model ,
the star points are lowered about 0.08 vertical units ( 20% ) , and the @xmath32 model curves are very slightly steeper functions of frequency .
however , the scatter among the observations is greater than the differences in the models between screen and extended geometries . in the vela plot
, we see that the @xmath32 model agrees somewhat better with the observations than does the simple kolmogorov model .
this pulsar is known to lie in a highly scattered region , and hence the presence of discontinuities , as incorporated in the @xmath32 model , is perhaps reasonable .
however , the conclusion is not strong , since the data are also reasonably fitted by an inner scale in the range @xmath117  km to @xmath118  km , as can be found from figure [ fig : vela]b .    for the other pulsars ( figure [ fig : other ] ) , the data show a stronger frequency dependence than predicted by both the simple kolmogorov and @xmath32 models .
however , the inconsistencies between the measurements at different frequencies are even greater , and better observations are needed , before such comparisons could discriminate between the models . a recent series of measurements at 327 mhz ( bhat et al .
1999a , b , c ) documents the variability of @xmath82 , and a similar long sequence of such measurements at other frequencies is needed on the same pulsars , before reliable conclusions can be reached from the frequency scaling observations .
slow variations in the flux of pulsars are caused by refractive interstellar scintillation ( riss ) , and can be characterized by the rms deviation in flux density normalized by its mean , or scintillation index , @xmath119 .
riss is due to inhomogeneities in the interstellar electron density on scales much larger than those responsible for diffractive scintillation , as characterized by the decorrelation bandwidth .
there are now measurements of both phenomena on a substantial number of pulsars , from which we can constrain the density spectrum on each line of sight .
the relation of the refractive scintillation index , @xmath119 , to the normalized diffractive decorrelation bandwidth , @xmath120 , depends on the ratio of the power in the density spectrum at the large refractive scales ( @xmath121  m to @xmath122  m ) to the power at the smaller diffractive scales ( @xmath123  m to @xmath81  m ) .
we compare pulsar measurements gathered from the literature with theoretical predictions by plotting @xmath119 versus @xmath120 ( cf .
rickett & lyne 1990 and gupta et al .
related tests have been made by armstrong et al .
( 1995 ) , bhat et al .
( 1999a ) and smirnova , shishov , & stinebring ( 1998 ) .
@xmath119 measurements at 610 mhz are listed in table [ table : mr610 ] and include the results of long - term monitoring observations by smirnova et al .
( 1998 ) , and older measurements near 100  mhz are listed in table [ table : mr100 ] .
the tables also include the decorrelation bandwidth obtained from diffractive scintillation , typically observed at a different radio frequency , scaled to the observing frequency for @xmath119 .
we used the kolmogorov scaling law @xmath124 , and we note that the minor differences that result from changes in the scaling law for other spectral models are insignificant in the comparison .
the 100  mhz data consist of @xmath119 measurements at 73.8  mhz , 81.5  mhz , and 156  mhz .
the pulsars observed at 610 mhz are primarily located at distances @xmath125  kpc , whereas the pulsars observed at 100 mhz are nearer at distances @xmath126  kpc .
figures [ fig : mr4 ] and [ fig : mrin ] show plots of the theoretical and measured refractive scintillation index versus the normalized decorrelation bandwidth in separate panels for the two sets of measurements .
theoretical curves are also shown in the figures for different spectral models in both the screen and extended medium geometries .
the dashed line corresponds to the simple kolmogorov model ( @xmath127 ) , and the solid lines in figures [ fig : mr4 ] and [ fig : mrin ] correspond to the @xmath32 and inner - scale models , respectively .    in section [ sec : decbw ] we noted that variable refraction tends to decrease the measured decorrelation bandwidth and used an expression from gupta et al .
( 1994 ) to apply a nominal correction to measured values . in the table of data from near 100 mhz ,
we have made use of the recently published measurements of bhat et al .
( 1999a ) .
they monitored the apparent decorrelation bandwidth near 327 mhz for 20 nearby pulsars , many of which have @xmath119 measurements in table [ table : mr100 ] .
they found the bandwidth to vary by factors 3 - 5 and discussed the influence of varying refraction as a refractive bias . from their measurements , they derived a corrected decorrelation bandwidth , which we have scaled ( using the @xmath84 scaling law ) to the frequency at which @xmath119 was observed .
these gave values typically 25 times bigger than obtained from earlier measurements ( e.g.  cordes 1986 ) .
the pulsar b0809 + 74 was not observed by bhat et al .
( 1999a ) ; however , the recent weak scintillation observations of this pulsar by rickett et al .
( 1999 ) similarly suggests that earlier decorrelation bandwidth measurements overestimated the strength of scattering .
the effect of these changes is to shift the plots to the right in figures [ fig : mr4 ] and [ fig : mrin ] by about half a decade , compared to figure 5 of gupta et al .
there is only one pulsar ( b0329 + 54 ) in table [ table : mr610 ] common to the bhat et al .
( 1999a ) observations , and its decorrelation bandwidth has also been corrected for the refractive bias .
since the pulsars observed at 610 mhz were mostly more heavily scattered , the refractive bias correction is substantially smaller ( see gupta et al .
1994 ) and has been ignored .
the theory of refractive scintillations has been known since the 1970s and was applied to pulsar flux variations by rickett et al .
for example , prokhorov et al .
( 1975 ) described how the modulation index for intensity can exceed unity in strong scattering in their equations ( 4.53 ) _ et .
we use the notation of coles et al .
( 1987 ) and confine our discussion to spherical wave sources , propagating in a scattering plasma which is either concentrated in a screen or extended uniformly between source and observer .
the `` low - frequency '' approximation for the intensity covariance function is given by equation ( 10 ) of coles et al .
( 1987 ) , from which we obtain the normalized refractive variance , @xmath128 , by setting the spatial offset equal to zero : @xmath129 dy   \right\ } \times   \nonumber \\ & & \sin^2\left [ 0.5 \kappa^2 r_{{\rm f},l}^2 x(1-x)\right ] \ , \kappa d\kappa \ , dx \ , \mbox { , } \label{eq : crextmed } \ ] ] where : @xmath130   \\
r_{{\rm f},l } & = & \sqrt{lc/2 \pi \nu } \mbox{. }   \nonumber \label{eq : h } \ ] ] here the transverse wavenumber is @xmath131 , and @xmath132 , where @xmath5 is the distance along the line of sight measured from the source . for the screen geometry , both the density spectrum and the gradient in the phase structure function @xmath133
are concentrated in a thin layer of thickness @xmath67 at distance @xmath68 from the source and @xmath93 from the observer .
thus the @xmath134 and @xmath135 integrations in equation ( [ eq : crextmed ] ) give @xmath136 , and we find it useful to define an equivalent fresnel scale @xmath137 : @xmath138 the @xmath139 integration in equation ( [ eq : crextmed ] ) involves the product of the density spectrum , decreasing as a steep power of @xmath139 , times the high - pass fresnel filter @xmath140 up to @xmath141 , times the low - pass exponential term that cuts off wavenumbers above @xmath142 , where @xmath143 is the radius of the scattering disc . hence in strong scattering ( @xmath144 ) we only consider the @xmath145 part of the fresnel filter .
consequently , the integration depends on the ratio @xmath146 or @xmath147 . for an extended scattering medium the line - of - sight integration softens the exponential cut - off in the @xmath139 integration , but the basic relationships remain the same . as others have noted the level of refractive scintillation is greater in an extended medium than in a screen with the same observed diffractive scintillation .
_    the canonical spectral model in iss studies of pulsars has been the simple kolmogorov spectrum @xmath148 . by fitting this model to diffractive scintillation observations
, observers have estimated the scattering measure , @xmath149 , toward many pulsars . when divided by the pulsar distance
, this gives a line - of - sight average of @xmath4 , which is found to vary greatly from one direction to another and to increase dramatically for distant pulsars seen at low galactic latitudes ( e.g.  cwb ) .
taylor & cordes ( 1993 ) have developed a smoothed model for the galactic plasma density distribution , which includes enhanced scattering in spiral arms and toward the galactic center
. however , there are also large random variations on much finer spatial scales , which would produce a scatter in a plot of @xmath119 against distance of dispersion measure . for a simple power - law spectrum both @xmath119 and @xmath120 depend on @xmath107 through a single strength of scattering parameter .
thus for the kolmogorov spectrum the variation of @xmath119 with @xmath120 is independent of distance or frequency and is given by a single dashed line in each plot in figures [ fig : mr4 ] and [ fig : mrin ] .
if the scattering medium is uniform , the strength of scattering is @xmath150 . for the screen it becomes @xmath151 , and in that case
, the dashed theoretical line is independent of the location of the screen between the source and the observer .
this point is substantiated in appendix  a , where the theory is laid out in more detail . according to this screen model , along each line of sight
there is a scattering layer with a certain @xmath107 , which determines both @xmath119 and @xmath120 , but @xmath107 is not necessarily related to the pulsar distance .
the theoretical curve with @xmath107 as the variable is independent of where the layer is located along the line of sight .    for the @xmath32 and inner - scale models , the theoretical curves for @xmath119 versus @xmath120 depend on the extra parameter @xmath53 or @xmath104 , respectively .
details of the theory are given in appendix  a , where the relevant parameters are shown to be @xmath152 and @xmath153 .
since @xmath70 depends on frequency , @xmath107 , and distance , then @xmath119 also depends somewhat on frequency and distance . in order to fix the frequency dependence of the theoretical @xmath119 values ,
we have separated the measurements into two groups ( at 610  mhz and near 100  mhz ) . the distance dependence is dealt with by assuming that @xmath154 distance .
whereas this is clearly appropriate for the extended scattering medium , it is less clear for the screen model , since the screen model supposes that @xmath107 is not necessarily related to distance .
however , it is reasonable that , if on a long line of sight there is a single region that dominates the scattering , its @xmath107 value will statistically increase with line - of - sight distance .
indeed the experimentally derived scattering measure increases faster with distance than if the medium were uniform ( see figure 1 of cordes et al .
( 1991 ) .    in figure
[ fig : mr4 ] we show the theoretical curves for @xmath119 in the @xmath32 model , which are relatively flat where @xmath155 . as @xmath97 increases ( @xmath120 decreases ) , @xmath156 increases , and when @xmath157 the scintillations are suppressed . in a screen
this occurs when @xmath158
.    inspection of figure [ fig : mr4 ] shows that the measured @xmath119 values at 100  mhz are above the prediction of the simple kolmogorov model for the extended medium geometry .
these measurements can , however , be modeled with the @xmath32 spectrum with suitable specific choices of outer scale , @xmath53 .
the screen geometry is inconsistent with the measurements at 100  mhz for both the simple kolmogorov and @xmath32 models .
for the 610  mhz data , the measured @xmath119 values are mostly in agreement with the prediction of the simple kolmogorov model for the extended medium geometry ; however , they lie below the predictions of the @xmath32 model . for the screen geometry ,
the measured @xmath119 values are above the simple kolmogorov model curve . however , for the @xmath32 model , though less convincingly , an agreement with the measurements can be found by suitable choices of the outer scale , @xmath53 .
the most striking feature of the figure is the good agreement with the uniform extended kolmogorov model for most of the pulsars .
this result , relies most heavily on the excellent data from smirnova et al .
( 1998 ) , who also note this agreement .
we are persuaded by this figure that the @xmath32 model can not be viewed as a global alternative to the kolmogorov spectrum .
we also conclude that for the medium - distance pulsars measured at 610  mhz , @xmath119 does not agree with the screen geometry .
it appears that even if the scattering medium is not uniform on scales of a few kpc , there is not a single region that dominates the scattering .
computations of @xmath119 due to a patchy distribution along the line of sight could test what distribution would start to approximate the uniformly extended medium .
turning to a comparison with the alternative inner - scale spectrum in figure [ fig : mrin ] , the 100  mhz observations could be explained by very large values of @xmath104 ( @xmath159 m to @xmath121 m ) in a screen configuration or by more modest , but still large , values of @xmath104 ( @xmath81 m ) in an extended medium .
the latter was essentially the proposal made by coles et al .
( 1987 ) for the nearby lines of sight .
the 610 mhz points are relatively lower than those at 100 mhz .
for an inner scale spectrum in a screen , large values of @xmath104 ( @xmath81 m to @xmath159 m ) would be necessary ; and in an extended medium 60% of the points lie near the simple kolmogorov model , with the rest requiring inner scale values ( @xmath123 m to @xmath81 m ) .
the 610 mhz data in our analysis come from riss observations of 21 pulsars made by stinebring and are described in more detail by smirnova et al .
( 1998 ) , who also compare the results with various models for the spectrum and spatial distribution of the electron density .
they note that 4 pulsars which are seen through known hii regions or supernovae remnants show relatively elevated values for @xmath119 , and find 3 other pulsars with similar behavior .
these are the points which lie above the kolmogorov line in figure [ fig : mrin]b .
their interpretation is that for these objects the scattering is concentrated into regions either near the pulsar or near the earth , and that these regions are characterized by an inner scale near @xmath160 m .
their spectrum model is very similar to our inner - scale model , except that the cut - off is characterized by a steep power - law rather than by an exponential function .
we find a somewhat larger numerical value for the inner scale required to match those objects .
we also note that the theoretical @xmath119 values are consistently higher for an extended medium , and so require a more modest inner scale .
however , an important result of our analysis is an alternative explanation for the relatively high @xmath119 values seen for some pulsars .
we suggest that pulsar lines of sight can pass through discrete clouds with increased plasma density on large scales , that steepen the low - wavenumber spectrum as opposed to cutting off the high - wavenumbers . a discontinuity spectrum ( @xmath32 )
is one way that the spectrum can be steepened , but smoother structures on scales larger than @xmath121 m ( @xmath161 au ) would also boost the low - wavenumbers and could cause similar enhancements of refractive compared to diffractive scintillations .
such enhancements are likely to be associated with hii regions or supernova remnants .
a similar idea was also discussed by lestrade , rickett , & cognard ( 1998 ) in the context of extreme scattering events in pulsar timing measurements .
to take this idea further , a composite between the kolmogorov and the @xmath32 or steeper spectra should be investigated , particularly one in which the line of sight is not uniformly weighted .
smirnova et al .
( 1998 ) also note that the strength of scattering increases very much faster with distance and dispersion measure than if the medium were statistically uniform .
the distant pulsars in their sample are mostly observed at low galactic latitudes and so are subject to the enhanced density and turbulence described by the taylor & cordes ( 1993 ) model .
it appears that most distant pulsars follow the uniform extended kolmogorov model with inner scale smaller than @xmath162 m .
thus although these lines of sight are subject to enhanced scattering in the inner galactic plane , the spectrum effectively follows the kolmogorov law and the plasma is dispersed enough to approximate a uniform scattering medium .
there are presumably discrete refracting structures present along these lines of sight , but their contribution to the density spectrum is masked by the higher densities of the inner galaxy , which still follow an apparently turbulent spectral form .
llr@.lr@.lr@.lr@.lr@.l    rickett & lyne ( 1990 ) & b0531 + 21 & 2&0 & 610&0 & 2&96@xmath163 & 0&327 & 0&052    & b0531 + 21 & 2&0 & 610&0 & 2&96@xmath163 & 0&402 & 0&095    kaspi & stinebring ( 1992 ) & b0329 + 54 & 1&4 & 610&0 & 0&44 & 0&39 & 0&07    & b0833@xmath8645 & 0&55 & 610&0 & 1&29@xmath164 & 0&11 & 0&01    & b1749@xmath8628 & 1&2 & 610&0 & 5&13@xmath163 & 0&26 & 0&05    & b1911@xmath8604 & 2&29 & 610&0 & 1&91@xmath165 & 0&20 & 0&03    & b1933 + 16 & 7&8 & 610&0 & 1&82@xmath165 & 0&18 & 0&04    & b2111 + 46 & 5&22 & 610&0 & 2&69@xmath165 & 0&16 & 0&05    & b2217 + 47 & 2&31 & 610&0 & 0&26 & 0&21 & 0&02    stinebring et al .
( 1999 ) & b0136@xmath16657 & 2&9 & 610&0 & 5&7@xmath165 & 0&15 & 0&02    & b0329@xmath16654 & 1&4 & 610&0 & 0&350 & 0&37 & 0&02    & b0525@xmath16621 & 2&3 & 610&0 & 0&30 & 0&31 & 0&01    & b0531@xmath16621 & 2&0 & 610&0 & 4&0@xmath163 & 0&32 & 0&01    & b0736@xmath8640 & 2&1 & 610&0 & 9&0@xmath167 & 0&03 & 0&01    & b0740@xmath8628 & 1&9 & 610&0 & 3&5@xmath165 & 0&13 & 0&01    & b0818@xmath8613 & 2&5 & 610&0 & 0&061 & 0&23 & 0&01    & b0833@xmath8645 & 0&5 & 610&0 & 1&50@xmath164 & 0&24 & 0&01    & b0835@xmath8641 & 4&2 & 610&0 & 7&0@xmath164 & 0&21 & 0&02    & b1641@xmath8645 & 4&6 & 610&0 & 7&0@xmath168 & @xmath1690&1 & 0&05    & b1642@xmath8603 & 0&5 & 610&0 & 0&770 & 0&46 & 0&04    & b1749@xmath8628 & 1&5 & 610&0 & 6&0@xmath163 & 0&25 & 0&02    & b1818@xmath8604 & 1&6 & 610&0 & 4&0@xmath164 & @xmath1690&1 & 0&05    & b1859@xmath8603 & 8&1 & 610&0 & 5&0@xmath167 & @xmath1690&05 & 0&03    & b1911@xmath8604 & 3&2 & 610&0 & 8&0@xmath164 & 0&22 & 0&01    & b1933@xmath16616 & 3&5 & 610&0 & 2&0@xmath165 & 0&24 & 0&02    & b1946@xmath16635 & 7&9 & 610&0 & 3&0@xmath170 & @xmath1690&05 & 0&03    & b2111@xmath16646 & 5&0 & 610&0 & 3&0@xmath165 & 0&15 & 0&01    & b2217@xmath16647 & 2&5 & 610&0 & 0&20 & 0&30 & 0&01    llr@.lr@.lr@.lr@.lr@.l    cole et al .
( 1970 ) & b0809 + 74 & 0&31 & 81&5 & 1&38@xmath165 & 0&45 & 0&18    & b0834 + 06 & 0&72 & 81&5 & 1&3@xmath165 & 0&39 & 0&11    & b1919 + 21 & 0&66 & 81&5 & 1&2@xmath165 & 0&53 & 0&18    helfand et al .
( 1977 ) & b0329 + 54 & 1&4 & 156&0 & 5&7@xmath165 & 0&35 & 0&08    & b0823 + 26 & 0&38 & 156&0 & 1&0@xmath163 & 0&46 & 0&11    & b1133 + 16 & 0&27 & 156&0 & 2&9@xmath163 & 0&50 & 0&08    & b1508 + 55 & 1&93 & 156&0 & 7&0@xmath165 & 0&41 & 0&07    & b1919 + 21 & 0&66 & 156&0 & 1&4@xmath163 & 0&50 & 0&07    & b2217 + 47 & 2&31 & 156&0 & 6&45@xmath164 & 0&41 & 0&11    gupta et al .
( 1993 ) & b0329 + 54 & 1&4 & 73&8 & 3&1@xmath164 & @xmath1710&15 &    & b0809 + 74 & 0&31 & 73&8 & 8&90@xmath164 & 0&34 & 0&06    & b0834 + 06 & 0&72 & 73&8 & 8&6@xmath164 & 0&16 & 0&03    & b0950 + 08 & 0&12 & 73&8 & 0&17 & 0&45 & 0&05    & b1133 + 16 & 0&27 & 73&8 & 1&6@xmath165 & 0&18 & 0&02    & b1237 + 25 & 0&56 & 73&8 & 1&8@xmath165 & 0&30 & 0&06    & b1508 + 55 & 1&93 & 73&8 & 3&8@xmath164 & 0&28 & 0&09    & b1919 + 21 & 0&66 & 73&8 & 7&5@xmath164 & @xmath1710&21 &
in this paper , the theory of the @xmath32 model for the electron density spectrum was derived for discontinuous density structures and compared with pulsar observations .
a new feature of our analysis is the inclusion of an `` outer scale '' needed in any realistic model .
the model is characterized by an effective exponent @xmath76 of the structure function , which remains between 1.95 and 1.6 over a very wide range of @xmath53 values ( cf .
figure [ fig : slope ] ) .
this at first seems a promising explanation for the spread in the estimates of @xmath79 derived from vlbi observations of the angular broadening profile , as observed on heavily scattered lines of sight .    as discussed in section [ sec : decbw ] , from figures [ fig : vela ] and [ fig : other ] , we find that the @xmath32 model provides a somewhat better agreement with the measurements of the diffractive decorrelation bandwidth versus frequency for pulsar psr b0833@xmath8645 ( vela ) than does the simple kolmogorov model
. this might arise from refractive scattering effects caused in the supernova remnant associated with the vela pulsar .
four other pulsars with decorrelation bandwidths measured against frequency show an appreciably stronger frequency dependence than the predictions of both the simple kolmogorov and @xmath32 models
. however , there are substantial inconsistencies among the measurements and better observations are clearly needed , especially in view of the variability in @xmath112 documented by bhat at al .
( 1999c ) .
the predictions of the @xmath32 model for the variation of the refractive scintillation index with the diffractive decorrelation bandwidth are in partial agreement with the observations . as discussed in section [ sec : mr ] , the values of @xmath119 measured near 100  mhz are above the prediction of the simple kolmogorov model and had previously been explained as the effect of an inner scale , substantially larger than the values invoked by spangler & gwinn 1990 .
however , the measurements are also consistent with the @xmath32 model with suitable choices of the outer scale .
for the 610  mhz data , most of the measured @xmath119 values are in good agreement with the prediction of the simple kolmogorov model for the extended medium , and they lie below the curves of the @xmath32 model . however , there are a significant number of good - quality observations , which lie somewhat above the kolmogorov line .
we suggest an alternative to a large inner scale on those lines of sight , that they pass through regions of enhanced density , which cause enhanced refractive scattering ; these regions must have less small - scale substructure than in a turbulent medium and could include discontinuities .    based on the above considerations
, we reject the @xmath32 model as a universal spectral model for the interstellar electron density fluctuations .
the corollary is to strengthen the evidence for the kolmogorov density spectrum , which in turn suggests a turbulent process in the interstellar plasma .
however , the simple kolmogorov spectrum is not a universal model either , since it disagrees with several of the @xmath119 observations .
since the @xmath32 model provides reasonable agreement for many of these discrepant observations , we propose that enhancements in the large scale part of the spectrum ( which need not be described by discontinuities in density ) occur on these lines of sight . with such enhancements causing the increase in refractive scintillation , there is no need to invoke the relatively large inner scales proposed by coles et al .
( 1987 ) . as proposed by spangler & gwinn ( 1990 ) ,
a relatively small inner scale is then likely , controlled by the ion inertial length or larmor radius .    [ cols= " <
, < " , ]     it appears that different spectral models need to be considered for different lines of sight .
a widely distributed turbulent plasma with occasional large ionized structures that increase the effective average power density , @xmath172 , at low wavenumbers ( large scales : @xmath121  m to @xmath173  m ) is thus a model that needs further formal investigation .
very similar conclusions have been reached by lestrade et al .
( 1998 ) and by bhat et al .
( 1999b , c ) .
this model could also explain the occasional `` extreme scattering events '' and episodes of fringes in dynamic spectra when a line of sight passes through a particular discrete density enhancement .
new theoretical work is needed to quantify the expected statistics of these propagation events .
it is likely that numerical modeling will be necessary to model non - stationary scattering media .
* acknowledgements : this work was supported by the nsf under grant ast-9414144 . *
here we derive expressions for the refractive modulation index , @xmath119 , as a function of the normalized diffractive decorrelation bandwidth , @xmath120 , for the inner - scale and @xmath32 models .
the simple kolmogorov model is included as a special case where the inner scale goes to zero .
a general integral expression for the refractive scintillation index is given by equation ( [ eq : crextmed ] ) .
with the medium concentrated into a thin layer ( thickness @xmath67 at @xmath68 from the source ) and an assumption of isotropy in the density spectrum , the equation becomes : @xmath174 \nonumber \\ & & \times \sin^2\left(\frac{r_{\rm f , scr}^2 \kappa^2}{2}\right ) \kappa \ , d\kappa \
, \mbox{. } \label{eq : app : mrscreen } \ ] ] for both the inner - scale and @xmath32 models , the exponential function in equation ( [ eq : app : mrscreen ] ) may be approximated by the gaussian ( coles et al .
1987 ) : @xmath175 \simeq \exp\left [ -(r_{\rm f , scr } u_{\rm scr})^2 \kappa^2 \right ] \ , \mbox { , } \label{eq : app : gauss } \ ] ] where @xmath137 is the fresnel scale and @xmath102 is the strength of scattering , as defined in section [ sec : mr ] .
we now substitute models for the density spectrum @xmath176 , given by equations ( [ eq : extpowerlaw ] ) ( with @xmath177 ) and ( [ eq : beta4 ] ) and obtain relations between @xmath119 and the scattering measure @xmath74 . to these
we add the connections between @xmath107 and @xmath70 , for each spectral model , from the equations for the phase structure functions ( 17 ) , ( 11 ) , and ( 21 ) of lr99 .    for the inner - scale model ,
we obtain : @xmath178^{(2-\alpha)/\alpha } \ , \mbox { , } \label{eq : app : sminscale } \ ] ] where @xmath179^{\alpha/(\alpha-2 ) } \ , \mbox{. } \label{eq : app : mualpha } \ ] ] then substituting for @xmath107 we obtain an integral for @xmath128 , which can be solved analytically using standard techniques ( see e.g.  appendix 2 of rickett 1973 ) : @xmath180^{(2-\alpha)/\alpha } \times \nonumber \\ & & \left [ ( \cos\psi)^{-\alpha/2 } \cos(\psi\alpha/2 ) - 1 \right ]   \ , \times \nonumber \\   & & \left [ ( l_i/2s_{\rm 0,scr})^2 + u_{\rm scr}^4 \right]^{\alpha /2 } \ , \mbox { , } \label{eq : app : mrscreeninscale } \ ] ] where @xmath151 , and @xmath181 is given by : @xmath182 \ , \mbox{. } \label{eq : app : thetae } \ ] ] we see that @xmath119 depends on @xmath102 with @xmath153 as a parameter .
similarly the diffractive decorrelation bandwidth , @xmath120 , depends on @xmath102 through equation ( [ eq : delta_nu_d ] ) with @xmath95 as a parameter , which is a very slow function of @xmath153 as discussed in section [ sec : decbw ] .    for the simple power - law model we set @xmath183 , and in the strong scattering limit ( @xmath184 ) , equation ( [ eq : app : thetae ] ) gives @xmath185 , and we obtain the asymptotic expression :
@xmath186 putting this together for the simple kolmogorov model with @xmath127 , we have the simple relation : @xmath187 , plotted as the thick dashed lines in the screen plots in figures [ fig : mr4 ] and [ fig : mrin ] .    for the general inner - scale model , with @xmath104 non - zero , we show computed curves in figure [ fig : mrin ]
. we can recognize two regimes .
consider first the case of small @xmath104 . with sufficiently strong scattering , @xmath188 , @xmath189 .
however , there is a complication in portraying this behavior in a plot versus @xmath190 for fixed inner scale , because of the variation of @xmath137 with distance , which is not specified by the horizontal variable @xmath190 .
we deal with this by obtaining an approximate scaling of @xmath137 with @xmath190 . when @xmath191 , equation ( [ eq : app : sminscale ] ) relates @xmath192 , or @xmath193 .
we now argue that on average @xmath154 pulsar distance @xmath194 , and eliminating @xmath107 , we see that , @xmath195 . with these scalings and @xmath188 ,
we see : @xmath196 here @xmath197 is the fresnel scale for a `` reference '' pulsar , and the last version in equation ( [ eq : mr.asymptote1 ] ) comes from using the kolmogorov exponent @xmath198 , which gives the asymptotic slope of 1/8 for curves at the lower left of figure [ fig : mrin ] . now with @xmath199 and @xmath137 fixed , let @xmath70 increase ( i.e.  less scattering ) until @xmath200 .
the inner scale is no longer important , and we get the same relation as for the simple kolmogorov spectrum : @xmath201 ; visible where the curves steepen with increasing @xmath190 .    for larger inner scales , the curves in figure [ fig : mrin ]
show a pronounced peak .
these occur for inner scales greater than the fresnel scale and will be associated with focusing and caustics .
since our treatment only includes the first order term in the low wavenumber expansion , it is not reliable in the region of the peak .
goodman et al .
( 1987 ) have discussed caustics at length for the same inner - scale spectrum .
they note that when @xmath202 , the scintillation power spectrum starts to fill in at scales intermediate between the diffractive and refractive wavenumbers . their equation ( 2.5.7 )
gives an estimate of the variance in this extra term as : @xmath203 with fixed @xmath104 , @xmath204 increases much more steeply with decreasing @xmath102 than does @xmath128 .
thus as @xmath205 decreases , there are fewer and fewer independent phase perturbations across the scattering disc .
when @xmath206 , focusing represented by @xmath207 occurs .
we note that the theoretical @xmath119 versus @xmath208 plots in gupta et al .
( 1993 ) omit the focusing condition and are wrong for inner scales greater than the fresnel scale .
the effect of higher order terms in the low - wavenumber expansion has been studied by dashen & wang ( 1993 ) .
they obtain a more efficient expansion scheme that gives improved accuracy near the peak in scintillation index .
nevertheless , it seems that a reliable prediction for the behavior near the peak in scintillations requires numerical evaluation .
this becomes even more necessary in treating an extended scattering medium .
we also note that the drop in @xmath119 as scattering gets weaker past the peak in figure [ fig : mrin ] is real .
it represents the fact that when a square - law structure function applies for scales from @xmath70 up to the scattering disc size , there is insufficient phase curvature and the scintillations remain weak even though @xmath144 . in such circumstances
the `` scattering disc '' is a misnomer , since an observer would see only a single angle of arrival , that could wander over a region of scale @xmath205 .    turning to the screen analysis of the @xmath32 model
, we relate @xmath70 to @xmath107 using equation ( 21 ) from lr99 , where the phase structure function equals one .
this gives the following analog of equation ( [ eq : app : sminscale ] ) for the inner - scale model : @xmath209 \right\}^{-1 } \mbox{. } \label{eq : app : smbeta4 } \ ] ] substituting the @xmath32 model for the density spectrum in equation ( [ eq : app : mrscreen ] ) and letting @xmath210 , we obtain : @xmath211 \ ; \sin^2 ( r_{\rm f , scr}^2 \eta/2 ) \ ; d\eta \label{eq : app : mrscreenb4_1 } \ ] ] in order to evaluate this integral , we let @xmath212 represent the integrand , and we let @xmath213 be the same as @xmath212 but with @xmath214 set equal to zero
. we can then rewrite the integral in equation ( [ eq : app : mrscreenb4_1 ] ) as : @xmath215 \right\}^{-1 }   \times \nonumber \\ & & \left\ { \int_0^{\infty } g_1(\eta )   d\eta - \int_0^{\infty } [ g_1(\eta ) - g(\eta ) ] d\eta \right\ } \label{eq : app : mrscreenb4_2 } \ ] ] the first integral can be evaluated analytically ( see e.g.   & 1965 ) . in strong scattering , the exponential term cuts off the oscillations of the sine term , which can be approximated by its argument , and the @xmath216 becomes negligible for values of @xmath217 larger than @xmath218 . with these approximations
we can also do the second integral . putting these together ,
we obtain : @xmath219   \right\}^{-1 } \times \nonumber \\ & & \left\ { 4 u_{\rm scr}^2 \tan^{-1}(u_{\rm scr}^{-2 } ) - 2 u_{\rm scr}^4\ln ( 1 + u_{\rm scr}^{-4 } ) -   \right .
\nonumber \\   & & \left .
2 \zeta \ , [ ( 2 + \zeta ) \
, e^{\zeta } \ , { \rm e_1}(\zeta ) - 1 ] \right\ } \ , \mbox { , } \label{eq : app : mrscreenb4 } \ ] ] where @xmath220 is the exponential integral , @xmath221 . in the latter form @xmath222 is the refractive scale ( equal to scattering disk radius ) .
again , @xmath119 is related to @xmath120 through : @xmath223 .
the resulting curves are shown in the screen panels of figure [ fig : mr4 ] .
consider the asymptotic behavior for @xmath128 in equation ( [ eq : app : mrscreenb4 ] ) as a function of @xmath102 .
as the strength of scattering @xmath102 increases , @xmath120 decreases ( @xmath224 ) . in equation ( [ eq : app : mrscreenb4_1 ] ) ,
the exponential term cuts off the integral ( at @xmath225 ) before the oscillations of the sin@xmath226 fresnel filter , which then approximates @xmath227 .
if also @xmath228 , we can ignore @xmath218 in the denominator and the remaining @xmath229 cancels the @xmath230 from the fresnel filter , and the integral depends only on @xmath231 . in this approximation , @xmath128 is then simply proportional to the slowly varying logarithmic term , which explains the relatively flat part of the curves in figure [ fig : mr4 ] ; under these conditions , in equation ( [ eq : app : mrscreenb4 ] ) the first two terms in the curly brackets sum to 2 and the last term is negligible . with @xmath53
fixed , now let @xmath102 increase , making @xmath156 increase . eventually @xmath232 becomes greater than one when the scattering disk becomes greater than the outer scale . at this point
the exponential term cuts off the integral below @xmath218 , where the spectrum flattens .
as @xmath102 increases even further , the integral decreases steeply , causing the down - turn at very small @xmath120 .
we again note that our expressions rely on the first order of an expansion and will not be reliable near the peak in the scintillation index
. however , there is not the same focusing condition that applied for very large inner scales .      in order to obtain expressions for @xmath119 for the two spectrum models in the _
uniform extended medium geometry , we must complete the line - of - sight integrals in equation ( [ eq : crextmed ] ) in addition to following the steps used in the screen geometry .
for each distance @xmath233 in the line of sight , there is also an integration over variable @xmath135 in the exponential cut - off . if @xmath234 , this @xmath135-integration yields @xmath235 , where @xmath236 .
this again provides a low - pass cut - off at the reciprocal of the radius of the effective scattering disk , where @xmath237 .
we define the scattering strength by @xmath238 , with the field coherence scale @xmath239 , as in lr99 , defined where the spherical wave structure function equals unity , measured in the observing plane .
for the other spectrum models there is not such a simple relation for the @xmath135-integration , but there is still an effective cutoff given by a similar equation .
when the @xmath134-integration is completed , the effective fresnel scale is actually smaller than @xmath240 due to averaging over the @xmath241 . _    in analogy with the screen geometry ,
we make use of identities similar to those given by equations ( [ eq : app : sminscale ] ) and ( [ eq : app : smbeta4 ] ) . for the extended medium
, these identities are in turn derived from the _ wave _ structure functions for the inner - scale and @xmath32 models ( cf .
the identities obtained thus will be similar to those for the screen geometry , except that for the inner - scale model , ( 1 ) there will be a factor of 3 on the right side of equation ( [ eq : app : sminscale ] ) , and ( 2 ) the 1 in square brackets is replaced by @xmath242^{\alpha/(\alpha-2 ) } \approx 1.8 $ ] for the kolmogorov exponent . for the @xmath32 model
, the only change will be a factor of 3 on the right side of equation ( [ eq : app : smbeta4 ] ) .
we use all of the aforementioned identities for the inner - scale and @xmath32 models and compute the @xmath134 integral numerically , since it can not be carried out analytically .
the shapes of the curves bear a close relationship to the screen results , though the extended medium values generally lie above the associated screen values at the same @xmath190 .
the second moment of intensity is needed to describe the fluctuations of intensity . under strong scintillation conditions ,
separate forms can be used for refractive and diffractive fluctuations , since their spatial scales differ by several orders of magnitude . in lr99 , as elsewhere in the iss literature , the correlation of diffractive scintillations is approximated by the squared magnitude of the second moment of the field , leading to the simple result that the spatial scale of the diffractive scintillations is equal to the scale where the phase structure function equals unity ( @xmath65 ) .
however , goodman and narayan ( 1985 ) showed that for steep spectra ( @xmath243 ) this is no longer the case and the diffractive scale can be larger than @xmath65 . here
we examine this question for the @xmath32 spectrum . we give the details for a phase screen with plane wave source , which are readily generalized to a spherical wave source .
the two - frequency intensity cross - spectrum at wavenumber @xmath244 for a screen at distance @xmath62 is given by the fourier - like integral equation ( 17 ) of ccffh .
this depends on the combination of structure functions @xmath245 , which for a plasma screen can be written as : @xmath246 where @xmath247 and @xmath248 are the two radio wavenumbers , @xmath249 is their geometric mean , @xmath250 is their arithmetic mean , @xmath251 , @xmath252 , and @xmath253 is evaluated at @xmath249 ; @xmath254 is a spatial offset which is the variable of integration .
consider first the single frequency case @xmath255 ) . in the limit of very large @xmath139 ,
the first four structure functions saturate and sum to zero .
the last two are equal and @xmath256 , which gives the simple diffractive limit mentioned above .
this is the zero - order term of an expansion , which is obtained in terms of the sum of the first four terms as a small quantity .
the zero order result requires full saturation , which requires @xmath257 . in diffractive scintillation @xmath258
; hence , the condition becomes that the refractive scale @xmath259 .
our concern here is to consider what happens when the diffractive @xmath139 is not large enough for saturation of @xmath253 . for shallow density spectra ( @xmath260 ) ,
small argument approximations to the structure function follow an exponent @xmath261 , and the zero - order term gives a good approximation even when @xmath262 . however , for steep spectra goodman and narayan ( 1985 ) showed that the leading term in the structure function follows a square law , which exactly cancels in @xmath245 ; the result is that the high wavenumber limit depends on the next term in the structure function expansion , which yields a diffractive scale that is greater than the scale @xmath65 ( defined by the square law term ) .
now we consider the case for @xmath32 model , when the scattering disc @xmath156 is smaller than the outer scale @xmath53 . here
we can approximate equation ( [ eq : struc4 ] ) by : @xmath263 as for the steep spectra , the leading term in the structure function follows a square law , which cancels when substituted into equation ( [ eq : v4 ] ) .
@xmath245 can then be approximated for large @xmath139 by expanding in @xmath264 .
the result is : @xmath265 } { s_0 ^ 2 \ln(4/s_0 ^ 2\kappa_o^2 ) } \ , \mbox{. } \label{eq : v4approx}\ ] ] here @xmath97 is the strength of scattering defined in equation ( [ eq : delta_nu_d ] ) ; @xmath266 is the angle between vectors @xmath139 and @xmath267 .
we note that for a kolmogorov spectrum in the high - wavenumber limit , @xmath245 also includes terms in @xmath268 , which would be accounted for in the higher order terms of the expansion . a result similar to equation ( [ eq : v4approx ] )
is given by dashen & wang ( 1993 ) , though for a 1-dimensional phase screen .
in considering the spectrum of intensity fluctuations ( eq .  17 of ccffh ) , one can show that the dominant wavenumber is approximately @xmath269 with these substitutions in @xmath245 , which we then set @xmath270 , we solve for @xmath271 ; this gives an approximate equation for the diffractive spatial scale @xmath272 : @xmath273 under the condition assumed in this approximation , @xmath274 , we find @xmath275 . however , in practice the ratio @xmath276 never becomes large . with a large outer scale ,
say @xmath277  pc , and typical observing conditions @xmath278 m and @xmath279 , we find @xmath280 . thus the diffractive scale could be 70% greater than @xmath65 and would slowly approach @xmath65 for smaller outer scales .    turning to the two - frequency intensity correlation ( @xmath281 ) in the diffractive limit of large @xmath139
, the results of ccffh still apply .
namely , that the last two terms of equation ( [ eq : v4 ] ) largely control the decorrelation versus frequency .
they group the remaining terms into a filter that depends only on @xmath139 and a smaller term that becomes the basis of the expansion .
the filter term was discussed by lr99 and shown to be important only as the strength of scattering decreases .
it is the last two terms in @xmath245 that determine the zero - order result , so we looked at the effect of the higher order terms .
the quantity that we are ultimately concerned with is the cross - correlation of intensity at offset frequencies at the same observing point .
this comes from the integral of the cross - spectrum .
equations ( 31 ) through ( 34 ) of ccffh give the zero and first order terms of the cross spectrum in terms of the spectrum of refractive index fluctuations in the layer .
for the @xmath32 spectrum we reduced these to a sum of confluent hypergeometric functions , which can be explicitly computed . for a sample observing condition we found that the higher order terms for the cross - spectrum itself were significant compared to the zero - order term
; however when integrated to give intensity cross - correlation , they only had a minor effect on the decorrelation bandwidth itself ( @xmath282% increase ) .
the reason for this appears to be the dominant effect of the last two terms in the @xmath245 summation with unequal frequencies .
to summarize we find a modest ( logarithmic ) increase in the diffractive scale relative to the field coherence scale @xmath65 , but that this remains less than a factor of 1.7 for the likely iss parameters .
this is accompanied by a smaller increase in the decorrelation bandwidth relative to the calculations of lr99 , which relied on the normal zero - order expansion at high wavenumbers .
this small offset in the decorrelation bandwidth is negligible compared to the measurement errors for the observations under consideration .
we assume that the conclusions reached here for a screen would also apply for an extended scattering medium .
stinebring , d. r. , smirnova , t. v. , hovis , j. , kempner , j. c. , myers , e. b. , hankins , t. h. , kaspi , v. m. , & nice , d. j. 1996 , in `` pulsars : problems and progress '' , ed .
johnston , s. , walker , m. a. , & bailes , m. , proceedings of iau colloquium 160 , asp conference series , 105 , 455