The Spectrum Power of Interstellar Plasma Inhomogeneities in the Direction of Eleven Radio Pulsars

We have analyzed two-dimensional correlation functions from the dynamic spectra of 11 pulsars using the archival data of the “Radioastron” project. The time-sections of these functions were approximated by exponential functions with a power \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha $$\end{document}. It is shown that this approximation describes the shape of the correlation function much better than the Gaussian. The temporal structure function \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D(\Delta t)$$\end{document} for small values of the delay \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta t$$\end{document}is a power law with an index \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha $$\end{document}. The spectrum power of spatial inhomogeneities of the interstellar plasma is related to the power of the structure function as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n = \alpha + 2$$\end{document}. We have determined the characteristic scintillation time and the power \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n$$\end{document} in the direction of 11 pulsars. In the direction of three pulsars (B0329+54, B0823+26, and B1929+10), the spectrum power of spatial inhomogeneities of the interstellar plasma turned out to be very close to the value for the Kolmogorov spectrum (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n = 3.67$$\end{document}). For other pulsars, it ranges from 3.18 to 3.86. It is shown that the measured scintillation parameters are significantly influenced by the duration of the observation session, expressed by its ratio to the characteristic scintillation time. If this parameter is less than 10, the parameter estimates may be biased: the values of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha $$\end{document} and the characteristic scintillation time \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{t}_{{{\text{scint}}}}}$$\end{document} may decrease.


INTRODUCTION
Inhomogeneities in the interstellar plasma scatter radio emission of space sources. Intense studies of scattering effects began with the discovery of pulsars, since pulsars are ultra-compact objects with coherent radiation and the influence of the source structure can be neglected. The scattering is manifested as image blurring, pulse broadening, intensity modulation over time and frequency (scintillation), and distortion of radio spectra. The parameters characterizing these phenomena are scattering angle , scattering time , characteristic scintillation time , and characteristic decorrelation band . Theoretical studies of scattering effects have established certain relationships between their parameters (see, e.g., [1][2][3]). They must be compared with the parameters measured in radio astronomical observations of pulsars.
In the analysis of scintillations, it three modes of averaging of observed values are considered: snapshot mode, averaging mode, and ensemble averaging mode. First, snapshot mode corresponds to the averaging time , i.e., the averaging time is much shorter than the characteristic scintillation time. Second, analysis in the averaging mode corresponds to the opposite case, . Here, as scintillation time, we mean the diffraction scintillations time, which in most cases ranges from a few seconds to tens of minutes at meter and decimeter wavelengths. Third, there is also refractive scintillation with a characteristic time of several weeks and months [4]. Analysis of the observed parameters over such time intervals is named ensemble averaging [5]. In this study, we will analyze scintillations in the normal averaging mode, . We will use the observational data obtained during the implementation of the scientific program of the ground-space interferometer "Radioastron." These data have been used already in other studies [6][7][8]. Dynamic spectra and two-dimensional correlation functions will be analyzed.
2. OBSERVATIONS Pulsar studies were an important part of the scientific program of the "Radioastron" ground-space 3 67 n The observations were carried out mainly at the frequency of 324 MHz and, in some cases, at the frequency of 1668 MHz. The receiver frequency bands were 16 and 32 MHz, respectively. These studies were mainly aimed at the properties and spatial distribution of interstellar plasma. The distances to effective scattering screens were determined by comparing the angular broadening with the characteristic time of the pulse scattering [6][7][8][9]. The pulsars selected for our study and their parameters are listed in Table 1. In this paper, we will use the dynamical spectra of pulsars defined earlier for other purposes. The details concerning construction of dynamic spectra may be found in the publications cited above. The individual spectrum for the given time is calibrated as (1) where and are the spectra obtained in time windows during the pulse of the pulsar and out the pulse, respectively. Each dynamic spectrum contains values, where is the number of frequency channels and is the number of spectra in the given observational sample. Usually, the time interval between successive spectra is equal to the period of the pulsar, but in some cases, averaging over several periods was performed in order to smooth out the intrinsic variations in the intensity of the pulsar from pulse to pulse. All these parameters are listed in Table 1. In our list, two pulsars (B1749-28 and B1933+16) were observed at the frequency 1668 MHz and the rest of objects-at the frequency 324 MHz.

THE ANALYSIS METHOD OF DYNAMIC SPECTRA
As it was already noted in the Introduction, for a justified assessment of the scattering parameters, we analyzed the dynamic spectra in the time interval . For this purpose, two-dimensional correlation functions of dynamic spectra were calculated over the entire observation interval. The two-dimensional correlation function was calculated as (2) where and are frequency and time delays, while , where Then, we used normalized correlation functions . Our task was to estimate the characteristic scintillation time . Following Cordes [5], we tried to determine the value of the scintillation time from the time-section of twodimensional correlation functions of dynamic spectra taken with zero frequency delay. Traditionally [5], this cross-section is approximated by a Gaussian and the half-width of this function at the level is taken as the scintillation time. Other researchers, who used correlation functions for data analysis [2,10,11], have pointed out deviations of these functions from the Gaussian form. We suggest to use a more universal function for approximation of the time-sec- , namely, an exponential function with an arbitrary power: (4) Here, -amplitude, -the parameter characterizing the width of , the half-width at the level is defined as . The parameters , , were determined by approximation of computed by function (4) for positive shifts with the step of 2 counts from made for exclusion of noise peak.

ANALYSIS RESULTS
Analysis results are presented in Table 2. The first column of the table gives the designation of the pulsar, the second one -the ratio of the time at which the dynamic spectrum was obtained ( ) and characteristic scintillation time ( ), and the third one shows the characteristic scintillation time ( ) determined by the relation . The next two columns of the Table compare the values of the root mean square deviations (RMS) for two options for approximation of the calculated ones: by the traditional Gaussian and by the exponent with an arbitrary power; the last column shows this power. The numbers in the parentheses give the value of the formal RMS error of determination of this parameter as a result of the least squares fit. These errors refer to the last significant digits of the measured values. In fact, these parameters change over time in a much wider range of values.
Cordes et al. [12] called attention to the so-called "relative error" in their estimates of the decorrelation band, which arises in the analysis of correlation functions. They estimated the relative error by the value , where is the total number of scintillation details (scintiles) in the spectrum and is the bandwidth of the receiver. This error, when is determined from the CCF, usually is 10-20%. We call attention to the fact that in determination of , along with random errors as the indicated relative error, a systematic error also appears, which is also associated with restrictions on the observation time and the receiver bandwidth. This error arises due to the cutting of individual scintiles, which is due to the limited duration of the observational scan in time and frequency ( ). We suggest to introduce a "cutting parameter" , which is related to the number of scintiles truncated by the edges of the observation field in time. As such a truncation parameter, we will use the ratio value of the duration of the observational scan and the value of the characteristic scintillation time ; the values of are given in column 2 of Table 2. To compare the quality of the approximations of calculated by simple Gaussians and by exponents with an arbitrary power, we show the residual deviations for these cases in Table 2. It is seen that these residual deviations are always noticeably smaller (from 1.5 to 8.8 times) for the approximations by exponents with an arbitrary power. Figure 1 shows, as an example of comparison of these methods of approximation, the data for two pulsars: It is clearly seen that approximation by exponents with an arbitrary power reflects the observed much better.

DISCUSSION
Interstellar scattering is closely related to the properties of the interstellar scattering medium (ISM), since the scintillation effects are caused by fluctuations of the electron density in the latter. An analysis of a large amount of observational data on pulsars has shown that pulsar scintillations are statistically consistent with the power-law spectrum of spatial inhomogeneities with an index (i.e., the Kolmogorov spectrum) in a very wide range of spatial scales [13,14]. The structure function of oscillations of the phase of pulsars radio emission is related to the spatial spectrum of the electron density, which is described by the NCF NCF NCF = .
3 67 n  , where is the spatial frequency and characterizes the degree of turbulence.
The analysis of the structure functions of some pulsars has shown that the spectrum for different local directions in the Galaxy may differ from the Kolmogorov spectrum [8,15]. The relationship between the index of the spectrum of inhomogeneity and the power in the time dependence of the structure function was discussed in several publications [3,16,17]. In the paper by Shishov et al. [3], it was shown that, for the power-law spectrum in the regime of saturated scintillations, the temporal structure function has the form for and . The structure function can be expressed in terms of the normalized correlation function as . For the case of approximating the correlation function by an exponent, we obtain the expression for small delays: . Thus, for small delays, the structure function is a power-law with an index and the spectrum power of spatial inhomogeneities of the interstellar plasma in this approximation will be . The powers are shown in the last column of Table 2. The numbers in parentheses correspond to the curve negotiation errors of the power function.
It turned out that in the direction of three pulsars (B0329+54, B0823+26, and B1929+10) out of 11 pulsars in our list, the spectrum exponent of spatial inhomogeneities of the interstellar plasma is very close to the value of the Kolmogorov spectrum ( ). The highest value was obtained in the direction to the pulsar B1133+16 and the smallest value ( )-in the direction to the pulsar B1933+16. It should be noted that the correct approximation of the correlation function (power-law) gives the correct value of the power . It can also be noted that the difference between the correlation function and the Gaussian shows that the internal scale of turbulence is smaller than the characteristic spatial scale of scintillations.
As indicated in the previous paragraph, a systematic error in the determination of the characteristic scintillation time t scint and in the determination of the indicator can be introduced by the effect of the limited observational time at which the dynamic spectra. We propose to characterize this effect by the parameter n = .
3 18 n n  Table 2. We will discuss this influence based on the data presented in Fig. 2 This pulsar has the greatest value of the measure of dispersion (DM = 158 kpc/cm 3 ) and is located at a considerable distance (3.7 kpc). We also note that the RMS for the Gaussian in relation to the RMS for the power function increases with decreasing parameter R.
The duration of the observation session also affects the determination of the characteristic scintillation time t scint . Figure 3 shows a comparison of the NCF time cross sections for the full (7000 s) and truncated (100 s) observational scans: it can be seen that with a decrease in R, the cross section shape changes and t scint decreases. The smaller the R, the greater the effect is. Apparently, the reason for the dependence of the estimates of the exponent and the characteristic scintillation time on the relative sample length is the insufficient number of diffraction spots in the analyzed dynamic spectrum to obtain a statistically correct estimate of the parameters. With a small T obs /t scint ratio, some of the diffraction structures are cut off both in time and in frequency if f dif is only several times smaller than the receiver bandwidth.
6. CONCLUSIONS 1. We have carried out the analysis of normalized two-dimensional correlation functions of dynamic spectra for 11 pulsars using the archival data of the "Radioastron" project. The time-sections of these functions were approximated by exponential functions with a power α. It was shown that these functions describe the shape of NCF much better than the Gaussian functions. For all sources, the power α and characteristic scintillation time t scint were determined.
2. Out of 11 pulsars in our list, in the direction of three pulsars (B0329+54, B0823+26, and B1929+10), the spectrum power of spatial inhomogeneities of the interstellar plasma turned out to be very close to the value for the Kolmogorov spectrum ( ). For other pulsars, it ranges from 3.18 (PSR B1933+16) to 3.86 (PSR B1133+16). The average value of this power NCF = .
3 67 n for the entire list, with the exception of the pulsars B2016+28 and B1933+16, is 3.62.
3. It was shown that the duration of the observation session, expressed in the units of the characteristic scintillation time (parameter ), has a noticeable effect on the measured parameters. If this parameter is less than 10, biased estimates of the parameters can be obtained: a decrease in the values of the power and the characteristic scintillation time . Out of our list, the pulsar B2016+28 definitely belongs to this category.

ACKNOWLEDGMENTS
The project "Radioastron" was carried out by the Astro-Space Center of the P.N. Lebedev Physics Institute of the Russian Academy of Sciences and S.A. Lavochkin Research and Production Association under a contract with the State Corporation ROSCOSMOS together with many scientific and technical organizations in Russia and other countries.
In this publication, we used the results of pulsars observations carried out for the "Radioastron" program with the following radio telescopes: 100-m radio telescope of the Green Bank Observatory, which is an institution of the National Science Foundation (NSF) and is operated in accordance with a cooperation agreement with the Association of Universities (Associated Universities, Inc.); the Arecibo Observatory 300-m radio telescope which was operated by SRI International under a collaborative agreement with the National Science Foundation (AST-1100968) and in alliance with Ana G. Mendes University (Metropolitana University and the University Space Research Association); as well as the 64-m Parkes Radio Telescope is a part of the Australian Telescope (ATNF), which is funded by the Australian government as a national facility operated by CSIRO.
We acknowledge E.N. Fadeev (ASC of the Physical Institute) for providing dynamic spectra from the personal archive and useful advice during preparation of the paper.

OPEN ACCESS
This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.