all radio signals from cosmic sources are distorted by the plasma turbulence in the interstellar medium ( ism ) .
understanding of this turbulence is therefore essential for the proper interpretation of astronomical radio observations .
the properties and characteristics of this turbulence can best be studied by observing point - like radio sources , where the results are not influenced by the extended structure of the source , but instead are directly attributable to the effect of the ism itself .
pulsars are such sources . dispersion and scattering affect radio emission from pulsars . whereas dispersion in the plasma column introduces delays in arrival time that depend upon frequency and results in smearing of the pulse , scattering by density inhomogeneities causes angular broadening , pulse broadening , intensity modulation or scintillation , and distortion of radio spectra in the form of diffraction patterns .
the scattering effects have already been studied extensively theoretically ( see , e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) and observationally with ground vlbi of sgra@xmath3 @xcite and pulsars ( see , e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) , as well as with ground - space vlbi of psr  b0329 + 54  ( _ halca _ , * ? ? ? * ) and the quasar 3c273 ( _ radioastron _ , * ? ? ?
whereas the vsop pulsar observations were done at a relatively high frequency of 1.7  ghz and with baselines of @xmath425,000  km and less , ground - space vlbi with _ radioastron
_  allows observations at one - fifth the frequency , where propagation effects are expected to be much stronger , and with baselines @xmath5 times longer @xcite .
such observations can resolve the scatter - broadened image of a pulsar and reveal new information about the scattering medium @xcite .    in this paper , we study the scattered image of the pulsar b0329 + 54 with _
radioastron_. we demonstrate that the pulsar is detected on baselines that fully resolve the scattering disk .
the interferometric visibility on these long baselines takes the form of random phase and amplitude variations that vary randomly with observing frequency and time . in the fourier - conjugate domain of delay and fringe rate ,
the visibility forms a localized , extended region around the origin , composed of many random spikes .
we characterize the shape of this region using averages and correlation functions .
we argue theoretically that its extent in delay is given by the average envelope of the impulse - response function of interstellar scattering , sometimes called the pulse - broadening function .
we find that the observed distribution is well - fit by a model that is derived from an impulse - response function that has two different exponential scales .
we discuss possible origins of the two scales .
llllr 2012 nov 26 through 29 & 1  hr / day & gb & rcp@xmath6lcp & 570  s + 2014 jan 1 and 2 & 12 hr & wb , kl & rcp & 1170  s [ tab : diary ]
our fundamental observable is the interferometric visibility @xmath7 . in the domain of frequency @xmath8 , this is the product of electric fields at two antennas @xmath9 and @xmath10 : @xmath11 this representation of the visibility is known as the cross spectrum , or cross - power spectrum .
because electric fields at the antennas are complex and different , @xmath12 is complex .
usually visibility is averaged over multiple accumulations of the spectrum , to reduce noise from background and the noiselike electric field of the source .
the second argument @xmath13 allows for the possibility that the visibility changes in time , as it does for a scintillating source , over times longer than the time to accumulate a single spectrum .
such a spectrum that changes in time is known as a `` dynamic spectrum '' @xcite .
the correlator used to analyze our data , as discussed in sections  [ sec : obs ] and  [ sec : red ] , calculates @xmath14 @xcite .
hereafter , we omit the baseline subscript indicating baseline @xmath15 in this paper , except in sections of the appendix where the baseline is important .     in various domains , and functions derived from it .
the fundamental observable is the visibility in the domain of frequency @xmath8 and time @xmath13 , @xmath16 ; this is known as the cross - power spectrum , or cross spectrum .
an inverse fourier transform of @xmath8 to delay @xmath17 leads to the visibility @xmath18 ; this is the cross - correlation function of electric fields in the time domain ( see equation  [ eq : vtau_def ] ) . a forward fourier transform of @xmath13 to fringe rate @xmath19 leads to @xmath20 .
a forward transform of @xmath17 back to @xmath8 produces @xmath21 , and an inverse fourier transform of @xmath19 to @xmath13 returns to @xmath16 the square modulus of @xmath18 is @xmath22 .
the cross - correlation function in @xmath17 of @xmath23 for right- and @xmath24 for left - circular polarization is @xmath25 .
we denote the fourier transform by @xmath26 , and quantities in the domain of frequency @xmath8 by the accent @xmath27 .
[ fig : ftcircle],scaledwidth=48.0% ]    under the assumptions that the source is pointlike , and that we can ignore background and source noise , the impulse - response function of interstellar scattering @xmath28 determines the visibility of the source .
a single delta - function impulse of electric field at the source is received as a function @xmath29 of time @xmath30 at the observer . here , @xmath30 is fourier - conjugate to @xmath8 and varies at the nyquist rate .
the visibility is the product of fourier transforms of @xmath28 at the two antennas : @xmath31 where @xmath32 is the fourier transform of @xmath29 .
we denote the typical duration of @xmath29 as @xmath33 , the broadening time for a sharp pulse . within this time span
, @xmath29 has a complicated amplitude and phase .
the function @xmath29 changes over longer times , as the line of sight shifts with motions of source , observer , and medium .
this change takes place on a timescale @xmath34 , and over a spatial scale @xmath35 .
the shorter and longer timescales @xmath33 and @xmath36 lead to our use of dual time variables : @xmath30 , of up to a few times @xmath37 and fourier - conjugate to @xmath8 ; and @xmath13 , of a fraction of @xmath36 or more and fourier - conjugate to @xmath19 .
this duality is commonly expressed via the `` dynamic spectrum '' ( see section  [ sec : appimpluseresponsefunction ] ) . if the scattering material remains nearly at rest while the line of sight travels through it at velocity @xmath38 , then one spatial dimension in the observer plane maps into time , and @xmath39 the averaged square modulus of @xmath28 is the pulse - broadening function @xmath40 . here , the subscripted angular brackets @xmath41 indicate an average over realizations of the scattering .
this function is the average observed intensity for a single sharp pulse emitted at the source .
an average over time is usually assumed to approximate the desired average over an ensemble of statistically - identical realizations of scattering .
we derive a number of representations of the visibility and quantities derived from it , and show that these provide straightforward means to extract the impulse - response function .
these functions are summarized in figure  [ fig : ftcircle ] , and discussed briefly here , and in detail in section  [ sec : impulse - response_visibility ] of the appendix . in particular , visibility in the domain of delay @xmath17 and time @xmath13
is @xmath18 .
this is the correlation function of electric field at the two antennas @xmath9 and @xmath10 ( equation  [ eq : vtau_def ] ) , and is the inverse fourier transform of @xmath42 from @xmath8 to @xmath17 .
we are also concerned with the square modulus of @xmath18 ( see  section  [ sec : acf_visibility ] ) : @xmath43 we calculate @xmath44 for right- and left - circular polarizations separately , and then correlate them in delay @xmath17 to form @xmath45 , the cross - correlation between polarizations : @xmath46 here , @xmath45 is the correlation of a single measurement of @xmath23 and @xmath24 , and @xmath47 is the number of samples in @xmath23 and @xmath24 .    when averaged over many realizations of the scattering material , @xmath48
is related to the statistics of the pulse - broadening function @xmath49 .
most commonly , the average over many realizations of scattering material is approximated by averaging over a time much longer than @xmath34 ; for this reason we omit the time argument for @xmath50 .
equivalently , evaluation of @xmath51 at the fringe rate @xmath52 of the maximum magnitude of @xmath45 yields the same time average . for this theoretical discussion , @xmath53
; for practical observations , instrumental factors can offset the fringe rate from zero , so that @xmath52 provides the most reliable time average .
for a baseline that extends much further than the scale of scattering @xmath35 ( see equation  [ eq : k_result ] ) : @xmath54 here , we introduce the symbol @xmath55 to indicate convolution , and denote the time - reverse of @xmath49 as @xmath56 .
our analysis method differs somewhat from @xcite , who used structure functions of intensity , visibility , and visibility squared to study scattering of pulsar b0950 + 08 on an extremely long baseline to _
radioastron_. the two methods are closely related theoretically .
structure functions are particularly valuable when the characteristic bandwidth approaches the instrumental bandwidth , and can be extended to cases where the signal - to - noise ratio is low , as they discuss .
lcc 2012 nov 26 & 60 & 60 + 2012 nov 27 & 90 & 60 + 2012 nov 28 & 175 & 60 + 2012 nov 29 & 235 & 60 + 2014 jan 1 & 20 & 60 + 2014 jan 2 & 70 & 100 + 2014 jan 2 & 90 & 120 [ tab : baselines ]
the observations were made in two sessions : the first for one hour each on the four successive days november 26 to 29 , 2012 , and the second for a total of 12 hours on the two days january 1 and 2 , 2014 .
the first session used the 10-m _ radioastron _
space radio telescope ( ra ) together with the 110-m robert c.  byrd green bank telescope ( gb ) .
the second session used the ra together with the @xmath0-m westerbork synthesis radio telescope ( wb ) , and the 64-m kalyazin radio telescopes ( kl ) .
both right ( rcp ) and left circular polarizations ( lcp ) were recorded in november 2012 , and only one polarization channel ( rcp ) was recorded in january 2014 . because of an ra peculiarity at 324  mhz , the 316332  mhz observing band was recorded as a single upper sideband , with one - bit digitization at the ra and with two - bit digitization at the gb , wb , and kl .
science data from the ra were transmitted in real time to the telemetry station in pushchino @xcite and then recorded with the _ radioastron _
data recorder ( rdr ) .
this type of recorder was also used at the kl , while the mk5b recording system was used at the gb and wb .
table  [ tab : diary ] summarizes the observations .
the data were transferred via internet to the astro space center ( asc ) in moscow and then processed with the asc correlator with gating and dedispersion applied @xcite . to determine the phase of the gate in the pulsar period , the average pulse profile was computed for every station by integrating the autocorrelation spectra obtained from the asc correlator .
the autocorrelation spectra @xmath57 are the square modulus of electric field at a single antenna .
in november 2012 the projected baselines to the space radio telescope were about 60 , 90 , 175 , and 235 thousand kilometers for the four consecutive days , respectively .
data were recorded in 570-second scans , with 30-second gaps between scans . in january 2014
the projected baselines were about 20 , 70 , and 90 thousand kilometers during the 12-hour session .
data were recorded in 1170-second scans .
the ra operated only during three sets of scans of 60 , 100 and 120 min each , with large gaps in between caused by thermal constraints on the spacecraft . the auto - level ( agc ) , phase cal , and noise diode were turned off during our observations to avoid interference with pulses from the pulsar
. table  [ tab : baselines ] gives parameters of the earth - space baselines observed .
all of the recorded data were correlated with the asc correlator using 4096 channels for the november 2012 session and 2048 channels for the january 2014 session , with gating and dedispersion activated .
the on - pulse window was centered on the main component of the average profile , with a width of 5  ms in the november 2012 session and 8  ms in the january 2014 session .
these compare with a 7-ms pulse width at 50% of the peak flux density @xcite .
the off - pulse window was offset from the main pulse by half a period and had the same width as the on - pulse window .
the correlator output was always sampled synchronously with the pulsar period of 0.714  s ( single pulse mode ) .
we used ephemerides computed with the program tempo for the earth center @xcite .
the results of the correlation were tabulated as cross power spectra , @xmath16 , written in standard fits format .
rrrrrcc nov 2012 & @xmath58 & @xmath59 & @xmath60 & @xmath61 & @xmath62 & @xmath63 + jan 2014 & @xmath64 & @xmath65 & @xmath66 & @xmath67 & @xmath68 &  [ tab : scatt_params ]      using autocorrelation spectra at gb , kl , and wb , we measured the scintillation time @xmath36 and bandwidth @xmath69 .
the results are given in table  [ tab : scatt_params ] .
our analysis using interferometric data , for which the noise baseline is absent and the spectral resolution was higher , is more accurate for the constants @xmath70 and @xmath71 as discussed below , so we quote those values in table  [ tab : scatt_params ] .
the asc correlator calculates the cross - power spectrum , @xmath16 , as discussed in sections  [ sec : theory_intro ] and  [ sec : visibility ] .
the resolution of the resulting cross - power spectra is 3.906  khz for the 2012 observations and 7.812  khz for the 2014 observations . because the scintillation bandwidth was comparable to the channel bandwidth for the 2014 observations , as shown in table  [ tab : scatt_params ] , and because the single recorded polarization at that epoch prevented us from correlating polarizations to form @xmath45 , as discussed in section  [ sec : two_scales ] , we focus our analysis and interpretation on the 2012 observations .
we investigated the scattering of the pulsar from the visibility in the delay - fringe - rate domain , @xmath20 .
we studied the statistics of visibility @xmath72 as a function of delay , fringe rate , and baseline length .
if there were no scattering material between the pulsar and the observer , we would expect for @xmath73 one spike at zero delay and fringe rate with magnitude that remains constant as a function of baseline length , and with width equal to the inverse of the observed bandwidth in delay , and the inverse of the scan length in fringe rate . scattering material in between changes this picture .
first we expect the spike at zero delay and fringe rate to decrease in magnitude with increasing baseline length , perhaps to the point where it would become invisible .
second , we expect additional spikes to appear around the spike at zero delay and fringe rate .
the distribution of these spikes give us invaluable information about the statistics of the scattering material .
as we discuss in this section , we fitted models to the distribution of visibility , as measured by the correlation function @xmath45 , and thus derived scintillation parameters that describe the impulse - response function for propagation along the line of sight from the pulsar .
we also computed the maximum visibility as a function of projected baseline length , as we discuss in detail in a separate paper ( paper ii : popov et al .
, in preparation ) .    for strong single pulses the visibility in the cross spectra , @xmath16 , had signal - to - noise ratios sufficiently large for a useful analysis
. however , we decided to analyze the data from the time series of multiple pulses .
fourier transform of the cross spectrum , @xmath16 , to the delay / fringe - rate domain yields @xmath74 and concentrates the signal into a central region , and thus provides a high signal - to - noise ratio .
the sampling rate of individual cross spectra in the time series was the pulse period of 0.714  s , as noted in section  [ sec : corr ] .
the time span of cross spectra used to form @xmath72 varied , ranging from 71.4  s to 570  s , depending upon the application .      in figure
[ fig : frrt_delay ] we display the magnitude of the visibility in the delay / fringe - rate domain , @xmath73 , for a 500-s time span .
the data were obtained on 29 november 2012 in the rcp channel for a projected @xmath75 gb - ra baseline .
the cross spectra , @xmath16 , from which we obtained @xmath73 were sampled with 4096 spectral channels across the 16-mhz band , at the pulsar period of 0.714 s ; consequently , the resolution was 0.03125  @xmath76s in delay , and 2  mhz in fringe rate .
as figure  [ fig : frrt_delay ] shows , no dominant central spike is visible at zero delay and fringe rate , as would be expected for an unresolved source .
our long baseline interferometer completely resolves the scattering disk .
instead we see a distribution of spikes around zero delay and fringe rate that is concentrated in a relatively limited region of the delay - fringe rate domain .
the locations of the various spikes appear to be random .
because the scattering disk is completely resolved on our long baseline , we conclude that the spikes are a consequence of random reinforcement or cancellation of paths to the different locations of the two telescopes , and hence interferometer phase .    in figure
[ fig : frrt_delay ] , the distribution of the magnitude of visibility is relatively broad along the delay axis and relatively narrow along the fringe rate axis .
the extent is limited in delay to about the inverse of the scintillation bandwidth , @xmath77 ; and in fringe rate to about the inverse of the diffractive timescale @xmath78 . within this region
, the visibility shows many narrow , discrete spikes .
if statistics of the random phase and amplitude of scintillation are gaussian , and the phases of the fourier transform randomize the different sums that comprise the visibility in the delay - fringe rate domain , then the square modulus of @xmath72 should be drawn from an exponential distribution , multiplied by the envelope defined by the deterministic part of the impulse - response function , as discussed in the appendix .    along the delay axis
, @xmath79 takes the general form suggested by figures  [ fig : frrt_delay ] and [ fig : delay_fine_structure ] : a narrow spike surrounded by a broad distribution .
we found that the central spike takes the form of a sinc function in both delay and fringe rate coordinates , as expected for uniform visibility across a square passband @xcite .
the widths are somewhat larger than values expected from observing bandwidth of 16  mhz and time span of 71.4  s , of @xmath80  ns and @xmath81  mhz respectively , probably because of the non - uniformity of receiver bandpasses and pulse - to - pulse intensity variations , respectively .
the broader part of the distribution takes an exponential form along the fringe - rate axis in this case ; more generally , the form can be complicated , particularly over times longer than 600  s. traveling ionospheric disturbances may affect the time behavior of our 92-cm observations ; in particular , they may be responsible for the 20 to 25  mhz width of the narrow component in fringe rate , as noted in table  [ tab : scatt_params ] .
we do not analyze the broader distribution in fringe rate further in this paper ; we will discuss this distribution , and the influence of traveling ionospheric disturbances , in a separate publication ( paper iii , popov et al . in preparation ) .
because of the relatively small optical path length of the ionosphere , even at @xmath82  cm , they can not affect the cross spectrum @xcite .
the distribution of the magnitude of the visibility in delay / fringe - rate domain changes with baseline length .
figure  [ fig : delay_fine_structure ] displays cross - sections through the maximum of the distribution of magnitude for a range of baseline lengths , as a function of delay . the maxima lie near zero fringe rate , as expected . under the plausible and usual assumption that the correct fringe rate lies at the fringe rate , @xmath83 , where the distribution peaks ,
the cross - section represents the visibility averaged over the time span of the sample : @xmath84 the top panel of figure  [ fig : delay_fine_structure ] shows this cross - section through figure  [ fig : frrt_delay ] .
the next lower panel shows the cross - section for the slightly shorter kl - ra baseline .
the three lower plots give the equivalent cross - sections for 10 times and 100 times shorter projected baselines .
these three short - baseline cross - sections are qualitatively different from the long baseline cross sections : the visibility has a central spike resulting from the component of the cross - spectrum that has a constant phase over frequency , as well as the broad distribution from the component that has a varying phase over frequency .
the central spike is strongest for the shortest baseline and weaker for the next longer baselines , as expected based on the results of sections  [ sec : visibility ] and  [ sec : acf_visibility ] . at very long baselines the central spike
is absent even after averaging the visibility over the whole observing period , and only the broad component is present .
as expected from figure  [ fig : frrt_delay ] , in the delay / time domain the broad component appears as spikes distributed over a range of about @xmath85s in delay .
these spikes keep their position in delay for the scintillation time of about 100 to 115  s , as listed in table  [ tab : scatt_params ] .
the character of the broad component changes with baseline length as well : mean and mean square visibility are the same for short and long baselines ; but excursions to small and large visibilities are more common for a long baseline ( * ? ? ? * eq .
12 ) .
averages of the visibility , and averages of the correlation function of visibility , extract the parameters of the broad and narrow components of visibility .
such averages approximate the statistical averages discussed in sections  [ sec : theory_intro ] and [ sec : impulse - response_visibility ] .
they seek to reduce noise from the observing system and emission of the source , as well as variations from the finite number of scintillations sampled , while preserving the statistics of scintillation .
the averages and correlation functions allow the inference of parameters of the impulse - response function of propagation from the statistics of visibility .
the mean square modulus of visibility , @xmath86 , provides useful and simple characterization of visibility . to approximate the average over realizations of scattering @xmath41 , we average over many samples in time @xmath13 and over bins in delay @xmath17 .
we realize the average over time by evaluating @xmath20 at the fringe rate of maximum amplitude @xmath52 , as discussed in section  [ sec : theory_intro ] .
we also average over 16 lags in delay @xmath17 .
the resulting average shows a broad component surrounding the origin ; on shorter baselines , it shows a spike at the origin .
the broad distribution samples the properties of the fine structure seen in figures  [ fig : frrt_delay ] and [ fig : delay_fine_structure ] , and the spike to those seen on the shorter baselines in figure [ fig : delay_fine_structure ] . we argue in section  [ sec : impulse - response_visibility ] that the spike in @xmath87 is related to the average visibility , and the broad component to the impulse - response function .    figure  [ fig : cross - section_fr - dl ] shows an example of the broad component of @xmath87 .
this is estimated as @xmath88 , by selecting the peak fringe rate @xmath52 to average in time for each of 6 scans , averaging the results for the scans , and averaging over 16 lags of delay to smooth the data .
these averaging procedures serve to approximate the average over an ensemble of realizations of scattering .
background noise adds complex , zero - mean noise to @xmath72 , with uniform variance at all lags ; this adds a constant offset to the average @xmath89 .     along the delay axis , at the fringe rate @xmath52 where the magnitude of visibility peaks , close to zero mhz .
the visibilities for the gb - ra baseline on 2012 nov 28 at 21:40 ut are shown as open circles .
the visibilities were computed by an inverse fourier transform of the spectra , @xmath42 , over 71.4 s time spans , and then by averaging over 6 observing scans , each 570  s long .
they were then further averaged in delay , over 16 points or @xmath90  @xmath76s , to smooth fluctuations .
the dashed horizontal line shows the offset contributed by background noise .
the solid gray line shows the reconstructed form given by equation [ eq : ctau_2exp_form ] , offset by the noise level , with parameters taken from the fit shown in figure  [ fig : ktwoscaleexamples ] .
the light dashed curve shows only the narrow component of the two - exponential model .
units of visibility are correlator units.,scaledwidth=47.0% ]      using equation  [ eq : krldef ] , we estimated @xmath92 , the averaged cross - correlation function between the square modulus of right - circular polarized ( rcp ) and of left - circular polarized ( lcp ) of visibility in the delay domain .
( note that @xmath93 is not the correlation function of the average @xmath94 , but rather the average of the correlation function @xmath95 . ) because the background noise in the two circular polarizations is uncorrelated , they do not contribute an offset to @xmath50 .
this allows us to follow the effects of the impulse - response function to much lower levels than for @xmath96 .
the correlation function @xmath50 is thus less subject to effects of noise , and is more sensitive to the broad component of the distribution , than @xmath96 .    to compute an estimate of @xmath92
, we calculated the squared sum of real and imaginary components of @xmath97 , the inverse fourier transform of the cross - power spectrum .
we formed these for each strong pulse , and normalized them by the autocorrelation functions at each antenna . from these we formed the un - averaged correlation function @xmath98
we then averaged @xmath45 over 570-sec scans to form @xmath99 . averaging in the time domain approximated an average over realizations of the impulse - response function for the scattering medium .
each 570-sec scan included 100 to 250 strong pulses , yielding one averaged sample of @xmath100 for each scan .
we obtained 22 measurements in total , with 6 samples of @xmath101 for november 26 , 28 , 29 observing sessions .
we obtained only 4 such samples for november 27 because of no significant detections of @xmath7 for two scans on that date .
examination of the averaged cross correlation function , @xmath102 , revealed a spike at the origin and two exponential scales for the broad component , a large one and a small one .
figure  [ fig : ktwoscaleexamples ] shows an example .    the spike at the origin arises from the fine structure of scintillation in the broad component of visibility , as seen in figures  [ fig : frrt_delay ] and  [ fig : delay_fine_structure ] .
this structure is identical in right- and left - circular polarizations , so its correlation leads to the spike .
the two exponential scales are apparent as the slopes of the steeper and narrower parts of the distribution .
we see these two scales even for single pulses , which are strong enough to show the two - scale structure .
we did not observe these scales without doubt in spectra from single - dish observations , because the resulting correlation functions are more subject to noise , gain fluctuations , and interference .
the scales are both present for @xmath96 , but the longer scale is seen more clearly in @xmath50 ( as comparison of equations  [ eq : ctau_2exp_form ] through [ eq : k_2exponential_model_params ] shows ) .
we formalized the two exponential scales seen for @xmath103 with a model fit .
the model assumed a pulse - broadening function with two exponential scales . under this assumption
, a short pulse appears at the observer with average shape : @xmath104 the pulse rises rapidly , and falls as the sum of the two exponentials .
the assumed form for @xmath49 leads to predictions for the forms of @xmath105 and @xmath106 , as discussed in section  [ sec : appendix2exponentials ] .
for @xmath105 , we expect a cusp at the origin , and two exponentials with scales @xmath107 and @xmath108 and different weights on either side . for @xmath106 , correlation smooths the cusp at the origin , producing a smooth peak , with the same exponential scales appearing to either side .
figure  [ fig : ktwoscaleexamples ] shows the best - fitting model of this form for the data shown there .
this model has parameters : @xmath109 the model reproduces the two scales , and the smooth peak , well .
the model also predicts the magnitude of the spike accurately , with zero average visibility @xmath110 .    the model shown in figure  [ fig : cross - section_fr - dl ] shows the model for @xmath44 , reconstructed using equation  [ eq : ctau_2exp_form ] with parameters from the fit to figure  [ fig : ktwoscaleexamples ] .
the two scales appear in the model , although the offset from noise contributes at large delay @xmath17 .
as the figure shows , a single exponential does not fit the model well : the narrow component is satisfactory at small @xmath17 , but falls well under the data at larger @xmath17 .
a high - winged function such as a lorentzian can fit @xmath44 well , but the rounded peak leads to a very wide peak for @xmath91 that can not match the data , and the inversion to a @xmath111 that remains finite , and is zero for @xmath112 as causality demands , is problematic .    the best - fitting scales and the magnitudes of the two contributions varied from scan to scan , but in a manner that was consistent with our finite sample of the scintillation pattern , and the inhomogeneous averaging of pulses with different intensities .
we show a histogram of the results of our fits to 570-sec intervals in figure  [ fig : scaledistributions ] . on 26 to 29 november 2012
, the shorter scale averaged to @xmath113 , and the longer scale to @xmath114 .
the scales had a relative power of @xmath115 .
on a long baseline that fully resolves the scattering disk , as figure  [ fig : delay_fine_structure ] shows , we observe multiple sharp spikes in the visibility @xmath20 is a consequence of the variation of the amplitude and phase of visibility .
( see also paper ii , popov et al . in preparation ) .
the characteristic region of that variation , @xmath116 , reflects the product of the inverses of the scintillation bandwidth @xmath117 and the scintillation timescale @xmath118 .
these quantities are the width in time of the impulse - response function , and the time for the impulse - response function to change as the line of sight to the observer moves through the scattering material .
detailed examination of the correlation function of visibility @xmath119 reveals the presence of two characteristic , exponential scales .
both scales are visible in the single - pulse correlation functions of right and left circular polarization , as well as in the correlation function @xmath120 averaged over 570  s shown in figure  [ fig : ktwoscaleexamples ] . for an assumed screen distance of half the pulsar distance of @xmath121  kpc @xcite ,
the two scales correspond to diffractive scales of : @xmath122 the diffractive scale is the lateral distance at the screen where phases decorrelate by a radian @xcite .
the refractive scale gives the scale of the scattering disk : @xmath123 in contrast , @xcite measured angular broadening for psr  b0329 + 54  of @xmath124  mas at @xmath125  mhz , where @xmath126 is the full width of the scattered image at half the maxium intensity .
this corresponds to a refractive scale of @xmath127 .
this upper limit is somewhat smaller than the values obtained from our observations , even if one takes into account the facts that the larger scale contains only 0.38 of the power of the shorter one , and that the scattering material may be somewhat closer to the pulsar than to the observer .
@xcite studied the scattering properties of psr  b0329 + 54  in detail , using single - antenna observations at 102  mhz , 610  mhz , 5  ghz , and 10.6  ghz to form structure functions of the scintillation in time and frequency on a wide range of scales .
they concluded that the scattering material has a power - law spatial spectrum with index @xmath128 , marginally consistent with the value of @xmath129 expected for a kolmogorov spectrum , with an outer scale of @xmath130 . using vlbi , @xcite observed psr  b0329 + 54  at 2.3  ghz and set limits on the separation on the emission regions corresponding to different components of the pulse profile . @xcite observed psr  b0329 + 54  at @xmath131  ghz with ground - space baselines to halca , and found that the source varied strongly with time .
they ascribed this variation to scintillation , with the scintillation bandwidth comparable to the observing bandwidth at their observing frequency .
self - calibration with timespans less than the scintillation time returned a pointlike image , as expected .
@xcite analyzed these data , including ground - ground baselines .
they studied both single - antenna autocorrelation functions @xmath132 and cross - correlation functions @xmath133 .
they detected two timescales for the scintillation pattern , of 20  min and 1  min .
they found that the properties of scattering could not be explained by a single , thin screen , and further that velocities indicated relative motions within the scattering medium .
@xcite had previously observed two coexisting scales of scattering for psr  b0329 + 54  .
they found scintillation bandwidths of @xmath134  hz and @xmath135  hz , measured as the @xmath136 point of the correlation function of intensity at an observing frequency of 102  mhz , using the large cophase array of puschino observatory .
the ratio of these scales , @xmath137 is larger than the ratio of @xmath138 that we observe . scaled to our observing frequency of 324  mhz , using the @xmath139 scaling appropriate for a kolmogorov spectrum , and converting from @xmath140 to @xmath17 using the uncertainty relation @xmath141 , we find that these values correspond to 1.3 and 8  @xmath76s , respectively , about a factor of 3 smaller that the scales we observe .
of course , interpolation over a factor of 3 in observing frequency and the different observing techniques may introduce biases , and scattering parameters likely vary over the years between the two measurements .
two scales of scattering have also been observed for other pulsars @xcite .
two scales of scattering may be a consequence of a variety of factors .
non - gaussian statistics of scattering can produce multiple scales , although this usually appears as a continuum of scales rather than two different individual scales , as in a power - law distribution or a levy flight @xcite . a kolmogorov model for scattering in a thin screen
does not fit as well as our model based upon a two - exponential impulse - response function , or even as one based upon one exponential . a model with two discrete scales
appears to fit our data better .
one explanation is anisotropic scattering .
this can produce two scales , corresponding to the major and minor axes of the scattering disk , as discussed in section  [ sec : anisotropic ] of the appendix .
the ratio of the scales of @xmath142 corresponds to the parameter @xmath143 , and an axial ratio of @xmath144 . in a simple model for anisotropic scattering in a thin screen
, we expect the ratio of power in the scales to be approximately @xmath145 , as shown in section  [ sec : anisotropic ] of the appendix .
this compares well with our observed ratio of @xmath146
. however , our observations for ground - space baselines at a variety of orientations do not show anisotropy .
a variety of models , involving material with varying anisotropy distributed along the line of sight , and strong anisotropy that slips between our long baselines , might match our data .
a second explanation is the complicated structure observed within dynamic spectra : most commonly observed as `` scintillation arcs '' @xcite .
recently , it has been suggested that this structure arises from interference among subimages , resulting from refraction by interstellar reconnection sheets @xcite .
this complicated structure produces time and frequency variations on a wide range of scales .
of course , we are considering very long baselines , where the scintillation - arc patterns should be completely uncorrelated between antennas .
this may lead to blurring , resulting in a 2-scale correlation function without particularly strong structure corresponding to the discrete arcs seen on shorter baselines @xcite .
we do not see any direct evidence of scintillation arcs , as such .
the magnitude of the visibility shows a featureless decline with increase of either of the 2 dimensions @xmath147 and @xmath148 .
the gb autocorrelation functions do not show scintillation arcs either , for our observations .
we made vlbi observations of psr  b0329 + 54  with _ radioastron _  at 324 mhz on projected baselines of up to 235,000  km .
our goal was to investigate scattering by the interstellar medium .
these properties affect radio observations of all celestial sources .
while the results of such observations are in general influenced by the convolution of source structure with the scattering processes , pulsars are virtually point - like sources and signatures in the observational results can be directly related to the scattering properties of the interstellar medium .    on long baselines , in the domain of delay @xmath17 and fringe - rate @xmath19 ,
the correlation function of visibility @xmath72 is a collection of narrow spikes , located within a region defined by the inverses of the scintillation bandwidth @xmath117 and the scintillation timescale @xmath118 . for shorter baselines ,
a sharp spike at the center of this region represents the average visibility ; on long baselines where the average visibility drops to near zero , this spike is absent .    the mean square visibility , @xmath149 , is well fit with a smooth model , indicating that the visibility spikes are the result of random interference of many scattered rays . to form a quantity less subject to effects of noise , we convolve the mean - square left- and right - circular polarized visibility to form @xmath150 .
the average correlation function @xmath50 shows two exponentials with different characteristic timescales .
the forms of @xmath151 and @xmath152 are well fit with a simple model , that assumes that the average pulse - broadening function @xmath49 is the sum of two exponentials with different timescales .
on 2012 nov 26 to 29 , the shorter timescale was @xmath153 , and the longer timescale was @xmath154 , with the longer - scale exponential containing approximately 0.38 times the power of the shorter - scale exponential .
this double exponential may arise from anisotropic scattering ; or from scattered radiation at large angle , perhaps corresponding to the subimages seen in single - dish and shorter - baseline observations .
further investigation of the properties of the image of the scattered pulsar on long and short baselines , using these data , will help to clarify the origin of the two scales .
the radioastron project is led by the astro space center of the lebedev physical institute of the russian academy of sciences and the lavochkin scientific and production association under a contract with the russian federal space agency , in collaboration with partner organizations in russia and other countries .
the national radio astronomy observatory is a facility of the national science foundation operated under cooperative agreement by associated universities , inc .
this study was supported by russian foundation for basic research grant 13 - 02 - 00460 and basic research program p-7 of the presidium of the russian academy of sciences .
c.r.g .  acknowledges support of the us national science foundation ( ast-1008865 ) .
n.b .  was supported by nserc .
we thank the referee for constructive comments which helped to improve the manuscript .    _
facilities : _ .
under general assumptions , refraction and scattering convolve the electric field of a source with an impulse - response function @xmath28 . this function varies with position in the observer plane , decorrelating over some lateral scale ; and with time , as the line of sight to the source moves with respect to the scattering material , and as the scattering material evolves .
the task of this section is to relate the impulse - response function to the statistics of visibility , as given by the functions @xmath44 and @xmath91 introduced in section  [ sec : theory_intro ] above .
the visibility @xmath155 is the conjugated product of electric fields at two antennas ( equation  [ eq : vdef ] ) .
we usually omit the subscripts indicating baseline @xmath15 on @xmath7 , unless they are important for the immediate argument .
we denote the fourier transform from the time or delay domain ( @xmath13 or @xmath17 ) to the frequency or fringe rate domain ( @xmath8 or @xmath19 ) by @xmath156 , and its inverse by @xmath157 .
we accent symbols with tilde `` @xmath27 '' to denote quantities that depend on observing frequency @xmath8 , and the same symbols without accent for the fourier - conjugate domain of delay @xmath17 or time @xmath30 .
we assume that the variables describing time and frequency @xmath158 , are discrete .
they range from @xmath159 to @xmath160 , where @xmath47 is the number of samples in the time or frequency span . for @xmath17 and
@xmath30 one sample is the inverse of the nyquist rate , and they can span the time to accumulate a single realization of the spectrum ; for @xmath13 and its fourier conjugate @xmath19 one sample is the averaging time for one spectrum , and they can span one observation .
our convention for normalization of the fourier transform is that a function @xmath161 normalized to unit area in the delay domain has value unity at zero frequency : @xmath162 .
conversely , if @xmath163 is normalized to unit area in the frequency domain , @xmath164 .
this is the `` @xmath165 '' convention of wolfram mathematica @xcite . with this convention
, parseval s theorem takes the form : @xmath166      as noted above , the observed electric field of a pulsar @xmath167 is the convolution of the electric field emitted at the source with a kernel @xmath28 that depends on scattering : @xmath168 where we introduce the symbol @xmath55 for convolution .
the kernel @xmath28 is the impulse - response function ; in other words , if the pulsar emits a sharp spike , then the observed electric field of the pulse is simply a copy of @xmath28 . because of this convolution
, @xmath28 is also known as the propagation kernel ; it is also known as the green s function , and the s matrix ( * ? ? ? * and references therein ) .
both @xmath169 and @xmath28 vary at the nyquist rate : the inverse of the total observed bandwidth .
usually , we assume that the intrinsic electric field of the source is white noise at the nyquist rate : it is drawn from a gaussian distribution in the complex plane at each instant @xcite .
the impulse - response function extends over a time span of a few times @xmath33 , representing the time over which a sharp pulse at the source would be received .
it is zero outside this relatively narrow time window .
if the statistics of the scattering material are stationary , the characteristic shape and scales of @xmath28 will remain fixed , while details of amplitude and phase vary on the timescale @xmath34 ( equation [ eq : t_sc ] ) .
an average of the squared electric field over many impulses emitted by the source over times longer than @xmath34 will reveal the characteristic form .
one simple model form for @xmath28 that includes deterministic and random parts is the product of a non - varying envelope @xmath170 , and a random function @xmath171 , that varies rapidly with @xmath30 during the course of each pulse : @xmath172 both @xmath171 and @xmath173 span a few times @xmath33 . over that time span
, @xmath171 varies wildly , and randomly ; however , it exhibits nearly the same form for the next pulse . over the longer timescale @xmath174 , the form of the random function @xmath171 changes slowly . for typical observations of a pulsar , such as those described in this paper
, @xmath33 is some fraction of the width of one pulse , or a few microseconds ; whereas @xmath34 is many pulsar periods , or many seconds .
such situations , where a convolution may have a slowly - varying kernel , are commonly treated as `` dynamic spectra '' ( see * ? ? ?
the intensity received by an observer for an electric - field impulse at the source , averaged over many such impulses with different realizations of the scattering material , is the square modulus of the deterministic part of @xmath28 , which we call @xmath49 : @xmath175 here , the subscripted angular brackets @xmath176 indicate a statistical average over realizations of the scattering medium , for example as approximated by an average over pulses spanning a time greater than @xmath34 .
often , @xmath49 is called the pulse - broadening function .
so that propagation kernel leaves the intensity of the source unchanged , when averaged over time , we set : @xmath177    for strong scattering , as is observed for most pulsars at most wavelengths , we expect that many different paths , with random amplitude and phase , will contribute to the received pulse at each instant @xmath30 .
therefore , we expect that the random part @xmath171 will have the statistics of a random walk at each instant : the observed electric field will be drawn from a circularly gaussian distribution in the complex plane , with zero mean . on the other hand , the deterministic part @xmath173 sets the standard deviation of @xmath28 ; it reflects how many paths , and with what strength , contribute at each delay .
this model for scintillation is closely related to the amplitude - modulated - noise ( amn ) model for pulsar emission @xcite . in this model ,
the electric field emitted by the pulsar is the product of noise , drawn from a zero - mean gaussian distribution in the complex plane , with a more - slowly varying envelope that determines the standard deviation of the noise at each instant .
we suppose that the random part of the propagation kernel is completely uncorrelated in time , at the nyquist rate , within its span of a few @xmath33 .
then , at a location `` @xmath9 , '' @xmath178 on the other hand , @xmath179 is nearly the same for each emitted pulse ; it changes only over the longer timescale @xmath34 .
the question of how this slower variation of @xmath171 with time depends upon baseline length is much more complicated , and we discuss it briefly below .
however , if the lateral separation of the two stations @xmath9 and @xmath10 is much greater than the scale of the scattering pattern , then the random parts of @xmath28 for the two stations , @xmath179 and @xmath180 , are completely uncorrelated .
as the previous discussion shows , the impulse - response function involves three timescales : @xmath171 changes at the nyquist rate ; @xmath173 varies over the typical span of the impulse - response function @xmath33 ; and the time for the random variations of @xmath171 to change is @xmath34 .
the dynamic cross - power spectrum provides a useful description for these different variations @xcite .
a single sample of the cross - power spectrum , when averaged over time less than @xmath34 , has the characteristic scale @xmath181 , resulting from the finite span of the impulse - response function and the uncertainty principle .
the time variation of the cross - power spectrum over times @xmath182 captures the changes of @xmath171 .
visibility @xmath183 in the domain of frequency @xmath8 and time @xmath13 is the product of the fourier transforms of electric fields at stations a and b ( see equation  [ eq : vdef ] ) : @xmath184 we suppose that each sample of the cross - spectrum is averaged over many realizations of the source electric field @xmath169 , over a time short compared with @xmath34 .
this reduces noise from the source and backgrounds .
one may represent the visibility in four domains , linked by fourier transforms of frequency @xmath8 to delay @xmath17 , and time @xmath13 to fringe rate @xmath19 , as figure  [ fig : ftcircle ] illustrates . in this paper
we are particularly concerned with @xmath185 , visibility in the domain of delay @xmath17 and time @xmath13 .
this is the fourier transform of @xmath186 .
the convolution theorem for fourier transforms shows that @xmath97 is the cross - correlation function of electric fields in the time domain : @xmath187 here , @xmath30 indexes individual samples of electric field , over a short interval near the index time of the measurement of the cross - power spectrum , @xmath13 .
visibility in the delay - rate domain , conjugate to the frequency - time domain , takes the form : @xmath188 \right ] \label{eq : vtauftvnu}\end{aligned}\ ] ] searches for interference fringes are often conducted in this domain : because absolute calibration of delay and fringe rate are usually impossible for very - long baseline interferometry , the peak of @xmath73 can be used to determine them @xcite .
visibility depends on the separation @xmath189 of stations @xmath9 and @xmath10 , as well as on delay and rate , or time and frequency . from equations  [ eq : g_conv ] and  [ eq : vtau_def ] , and the assumption that the electric field of the source @xmath169 is a stationary random variable without correlation in time , we find that visibility in the delay domain is the cross - correlation function of @xmath28 at the two stations : @xmath190 this leads to the expected form of @xmath191 : a spike at @xmath192 , with average magnitude equal to the average correlation @xmath193 ; and a broad component of width @xmath33 , with random amplitude , and phase variations that increase with baseline length , corresponding to the random character of @xmath171 and its decorrelation with increasing baseline . equations  [ eq : normalize_gra ] and  [ eq : normalize_gd ] show that @xmath194 for @xmath195 , and @xmath196 for uncorrelated @xmath179 and @xmath180 .
figures  [ fig : frrt_delay ] and [ fig : delay_fine_structure ] show examples .    on intermediate baselines ,
the time structure of the correlation of @xmath171 is more complicated , in a way that depends on the geometrical distribution of the paths that contribute to @xmath171 . for scattering material concentrated in a thin screen , for example
, the shortest - length paths result in small @xmath30 in the impulse - response function , and also tend to appear at small angles at the observer .
thus , at small delays correlation is high even for rather long baselines ; whereas at long delays correlation is poor even for shorter baselines .
thus , correlation between antennas should decrease at later times @xmath30 within @xmath171 .
this correlation is imprecise for scattering material distributed along the line of sight , where many deflections along the line of sight lead to a large time lag @xmath30 , but little or no angular deflection at the observer s interferometer .
moreover , in the frequency domain , dynamic single - dish spectra can show slants and complicated patterns @xcite , suggesting complicated correlations of time and delay in the observer plane .
analysis of the visibility is thus easiest on very short baselines and very long ones .
equations  [ eq : vtauftvnu ] and  [ eq : vabtau_as_g ] provide the relation of the visibility in the frequency - time domain to the impulse - response function : @xmath197 averaging the visibility over many scintillations yields the average visibility , @xmath193 : @xmath198      the average of @xmath191 over many realizations of scintillations leaves the delta - function at the origin that corresponds to the average visibility : @xmath199 here , @xmath200 is the kronecker delta - function , with value 1 if @xmath192 and 0 otherwise .
the secondary spectrum @xmath201 may be defined as the square modulus of @xmath72 . for zero baselines
@xmath16 is real , so that @xmath202 , and our and their expressions are identical .
their expression includes phase information in an elegant way for their short baseline , where departures from zero phase are small .
for observations on long baselines , their expression is impractical because identification of the origin of @xmath203 is not possible , as figure  [ fig : frrt_delay ] shows .
consequently , the pair @xmath203 and @xmath204 can not be combined reliably .
they also use the accent @xmath27 for visibility in the delay domain , although not for @xmath44 ; whereas we use the accent for quantities in the frequency domain of @xmath8 . ] : @xmath205 this function provides information similar to @xmath73 , shown in figure  [ fig : frrt_delay ] , but is easier to deal with statistically . because @xmath206 is the square modulus of the complex visibility @xmath20 , background noise adds noise to @xmath44 , except for the central lag , where it may contribute a constant offset .
similarly , self - noise will add noise with an envelope that follows the average form of @xmath44 .
this behavior is in contrast to that of correlation functions of single - dish quantities such as the intensity , where noise can contribute to the mean correlation function .
our long baselines fully resolve the scattered image .
the phases of scintillation elements in @xmath16 appear to be random , and the phases of @xmath201 show no discernible patterns .
an inverse fourier transform from fringe rate @xmath19 to time @xmath13 leads to @xmath22 : @xmath207 .\end{aligned}\ ] ] evaluated at the fringe rate @xmath52 of its peak magnitude , the secondary spectrum @xmath208 is a time average , that approximates an average over realizations of the scintillation pattern of @xmath22 : @xmath209 the autocorrelation function of @xmath44 is @xmath91 : @xmath210 note that @xmath22 appears without averaging in this expression .
conveniently , the correlation @xmath211 between the secondary spectra @xmath44 in the right and left circular polarizations eliminates some effects of noise and interference , as noted in section  [ sec : two_scales ] .
the behaviors of @xmath44 and @xmath91 are simplest to describe on very short baselines , where @xmath212 for all @xmath213 , and @xmath194 ; and for very long baselines , where the random parts of @xmath214 and @xmath215 are completely uncorrelated , so that @xmath110 .
if the correlation of the random parts of the propagation kernels @xmath193 is constant , then the ensemble - average values of these correlation functions are : @xmath216 again , @xmath200 is the kronecker delta - function , with value 1 if @xmath192 and 0 otherwise ; and similarly for @xmath217 .
the time - reversed pulse - broadening function is @xmath218 , given by @xmath219 .
the constant @xmath220 is the mean square of @xmath49 : @xmath221 .
we consider a situation where the impulse - response function is the sum of two exponentials with different time constants @xmath107 and @xmath108 , with a rapid rise from @xmath222 at @xmath223 .
as we discuss below , this may result in a variety of circumstances .
we parametrize the impulse - response function : @xmath224 in the text , we also make use of the inverse scales @xmath225 and @xmath226 ; these can provide better physical insight .
the autocorrelation of @xmath49 provides the form for @xmath227 , the mean square visibility in the delay - time domain , as given by equation  [ eq : c_result ] . for this impulse - response function , under the assumption that the baseline is
so long that @xmath110 , this takes the form : @xmath228 the autocorrelation function @xmath91 of the mean square visibility for this impulse - response function then takes the form : @xmath229 where again @xmath230 is the kronecker delta function , and : @xmath231 thus , the two exponential scales @xmath232 appear again , in @xmath91 , as do the weights @xmath233 , @xmath234 .
an observer sees an anisotropic distribution of radiation from a screen at distance @xmath235 from the observer .
the probability of receiving radiation from the screen at position @xmath236 is : @xmath237 we suppose without loss of generality that @xmath238 .
if the source is at infinite distance beyond the screen , the delay @xmath17 along this path ( ignoring any contribution from the screen ) is : @xmath239 where @xmath240 is the speed of light .
if the source is at distance @xmath241 beyond the screen , then @xmath235 is replaced by @xmath242 in this and subsequent equations .
we convert the distribution of angles in equation  [ eq : anis_angle_dist ] to coordinates @xmath243 , where @xmath244 .
the resulting distribution of @xmath245 is : @xmath246 where @xmath247 parametrizes the anisotropy .
we integrate over @xmath248 to find the pulse - broadening function : @xmath249 thus , @xmath250 here , @xmath251 is the regular modified bessel function of order 0 , and @xmath252 .
note that for @xmath253 , this distribution becomes the familiar exponential distribution , with scale @xmath71 , as expected for @xmath254 .
@xcite present a similar expression .
in the case of @xmath255 , at small @xmath17 the distribution has the limit : @xmath256 at large values , @xmath257 so that @xmath258 at a particular scale @xmath17 , the logarithmic derivative of @xmath259 is @xmath260 , although the coefficient depends weakly on @xmath17 .
thus , @xmath259 exhibits two exponential scales : @xmath71 at large @xmath17 , and @xmath261 at small @xmath17 .
the relative strength of the scales is about @xmath262 near the larger scale @xmath71 , with the larger scale being the weaker .