We report on the substructure of visibility functions in the delay domain of PSRs B0329+54, B0823+26, B0834+06, B1933+16 and B0833-45 (Vela) observed with earth-earth and RadioAstron space-earth two-element interferometers at frequencies of 324 MHz and 1668 MHz. All visibility functions display unresolved spikes distributed over a range of delays. They are due to band-limited scintillation noise and related to the scattering time. The envelopes for each but the Vela pulsar are well fit by a single Lorentzian which we interpret as being indicative of isotropic scattering on the plane of the sky due to a thin scattering screen between the pulsar and us. In contrast, the envelope for the Vela pulsar needs to be mostly fit by at least two Lorentzians, a narrow and a broad one at the same zero delay. We interpret this characteristic as indicative of anisotropic scattering due to more complex structure of scattering screens in the supernova remnant. The possibility of describing the delay visibility functions by Lorentzians is likely a general property of pulsars and offers a new way of describing scattering parameters of the intervening interstellar medium. Furthermore, for all our pulsars, the unresolved spikes in visibility functions of similar projected baselines were well correlated indicating that the telescopes are located in the same diffraction spot. The correlation vanished for visibilities from largely different baselines, when some radio telescopes are not in the same spot.

We show that narrow filaments or sheets of over- or under-dense plasma, or "noodles," with fluctuations of scattering phase of less than a radian, can form the scintillation arcs seen for many pulsars. The required local fluctuations of electron density are indefinitely small. We assume a cosine profile for the electron column and find the scattered field by analytic Kirchhoff integration. For a large electron column, corresponding to large amplitude of phase variation, the stationary-phase approximation is accurate; we call this regime "ray optics". For smaller-amplitude phase variation, the stationary-phase approximation is inaccurate or inapplicable; we call this regime "wave optics". We show that scattering is most efficient when the width of the strip equals that of one pair of Fresnel zones, and in the wave-optics regime. We show that the resolution of present observations is about 100 Fresnel zones on the scattering screen. Incoherent superposition of strips within a resolution element tends to increase the scattered field. We find that observations match a single noodle per resolution element with phase of up to 12 radians; or many noodles per resolution element with arbitrarily small phase variation each, for net phase of less than a radian. Observations suggest a minimum radius for noodles of about 650 km, comparable to the ion inertial scale or the ion cyclotron radius in the scattering plasma.

I show that narrow, parallel strips of phase-changing material, or "noodles," generically produce parabolic structures in the delay-rate domain. Such structures are observed as "scintillation arcs" for many pulsars. The model assumes the strips have widths of a few Fresnel zones or less, and are much longer than they are wide. I use the Kirchhoff integral to find the scattered field. Along the strips, integration leads to a stationary-phase point where the strip is closest to the line of sight. Across the strip, the integral leads to a 1D Fourier transform. In the limit of narrow bandwidth and short integration time, the integral reproduces the observed scintillation arcs and secondary arclets. The set of scattered paths follows the pulsar as it moves. Cohorts of noodles parallel to different axes produce multiple arcs, as often observed. A single strip canted with respect to the rest produces features off the main arc. I present calculations for unrestricted frequency ranges and integration times; behavior of the arcs matches that observed, and can blur the arcs. Physically, the noodles may correspond to filaments or sheets of over- or under-dense plasma, with a normal perpendicular to the line of sight. The noodles may lie along parallel magnetic field lines that carry density fluctuations, perhaps in reconnection sheets. If so, observations of scintillation arcs would allow visualization of magnetic fields in reconnection regions.

Radio images of the Galactic Center supermassive black hole, Sagittarius A* (Sgr A*), are dominated by interstellar scattering. Previous studies of Sgr A* have adopted an anisotropic Gaussian model for both the intrinsic source and the scattering, and they have extrapolated the scattering using a purely $\lambda^2$ scaling to estimate intrinsic properties. However, physically motivated source and scattering models break all three of these assumptions. They also predict that refractive scattering effects will be significant, which have been ignored in standard model fitting procedures. We analyze radio observations of Sgr A* using a physically motivated scattering model, and we develop a prescription to incorporate refractive scattering uncertainties when model fitting. We show that an anisotropic Gaussian scattering kernel is an excellent approximation for Sgr A* at wavelengths longer than 1cm, with an angular size of $(1.380 \pm 0.013) \lambda_{\rm cm}^2\,{\rm mas}$ along the major axis, $(0.703 \pm 0.013) \lambda_{\rm cm}^2\,{\rm mas}$ along the minor axis, and a position angle of $81.9^\circ \pm 0.2^\circ$. We estimate that the turbulent dissipation scale is at least $600\,{\rm km}$, with tentative support for $r_{\rm in} = 800 \pm 200\,{\rm km}$, suggesting that the ion Larmor radius defines the dissipation scale. We find that the power-law index for density fluctuations in the scattering material is $\beta < 3.47$, shallower than expected for a Kolmogorov spectrum ($\beta=11/3$), and we estimate $\beta = 3.38^{+0.08}_{-0.04}$ in the case of $r_{\rm in} = 800\,{\rm km}$. We find that the intrinsic structure of Sgr A* is nearly isotropic over wavelengths from 1.3mm to 1.3cm, with a size that is roughly proportional to wavelength. We discuss implications for models of Sgr A*, for theories of interstellar turbulence, and for imaging Sgr A* with the Event Horizon Telescope.

We carried out observations of pulsar PSR B1919+21 at 324 MHz to study the distribution of interstellar plasma in the direction of this pulsar. We used the RadioAstron (RA) space radiotelescope together with two ground telescopes: Westerbork (WB) and Green Bank (GB). The maximum baseline projection for the space-ground interferometer was about 60000 km. We show that interstellar scintillation of this pulsar consists of two components: diffractive scintillations from inhomogeneities in a layer of turbulent plasma at a distance $z_{1} = 440$ pc from the observer or homogeneously distributed scattering material to pulsar; and weak scintillations from a screen located near the observer at $z_{2} = 0.14 \pm 0.05$ pc. Furthermore, in the direction to the pulsar we detected a prism that deflects radiation, leading to a shift of observed source position. We show that the influence of the ionosphere can be ignored for the space-ground baseline. Analysis of the spatial coherence function for the space-ground baseline (RA-GB) yielded the scattering angle in the observer plane: $\theta_{scat}$ = 0.7 mas. An analysis of the time-frequency correlation function for weak scintillations yielded the angle of refraction in the direction to the pulsar: $\theta_{ref, 0}$ = 110 ms and the distance to the prism $z_{prism} \le 2$ pc.

We have resolved the scatter-broadened image of PSR B0329+54 and detected substructure within it. These results are not influenced by any extended structure of a source but instead are directly attributed to the interstellar medium. We obtained these results at 324 MHz with the ground-space interferometer RadioAstron which included the space radio telescope (SRT), ground-based Westerbork Synthesis Radio Telescope and 64-m Kalyazin Radio Telescope on baseline projections up to 330,000 km in 2013 November 22 and 2014 January 1 to 2. At short 15,000 to 35,000 km ground-space baseline projections the visibility amplitude decreases with baseline length providing a direct measurement of the size of the scattering disk of 4.8$\pm$0.8 mas. At longer baselines no visibility detections from the scattering disk would be expected. However, significant detections were obtained with visibility amplitudes of 3 to 5% of the maximum scattered around a mean and approximately constant up to 330,000 km. These visibilities reflect substructure from scattering in the interstellar medium and offer a new probe of ionized interstellar material. The size of the diffraction spot near Earth is 17,000$\pm$3,000 km. With the assumption of turbulent irregularities in the plasma of the interstellar medium, we estimate that the effective scattering screen is located 0.6$\pm$0.1 of the distance from Earth toward the pulsar.

We discovered fine-scale structure within the scattering disk of PSR B0329+54 in observations with the RadioAstron ground-space radio interferometer. Here, we describe this phenomenon, characterize it with averages and correlation functions, and interpret it as the result of decorrelation of the impulse-response function of interstellar scattering between the widely-separated antennas. This instrument included the 10-m Space Radio Telescope, the 110-m Green Bank Telescope, the 14x25-m Westerbork Synthesis Radio Telescope, and the 64-m Kalyazin Radio Telescope. The observations were performed at 324 MHz, on baselines of up to 235,000 km in November 2012 and January 2014. In the delay domain, on long baselines the interferometric visibility consists of many discrete spikes within a limited range of delays. On short baselines it consists of a sharp spike surrounded by lower spikes. The average envelope of correlations of the visibility function show two exponential scales, with characteristic delays of $\tau_1=4.1\pm 0.3\ \mu{\rm s}$ and $\tau_2=23\pm 3\ \mu{\rm s}$, indicating the presence of two scales of scattering in the interstellar medium. These two scales are present in the pulse-broadening function. The longer scale contains 0.38 times the scattered power of the shorter one. We suggest that the longer tail arises from highly-scattered paths, possibly from anisotropic scattering or from substructure at large angles.

Earth-space interferometry with RadioAstron provides the highest direct angular resolution ever achieved in astronomy at any wavelength. RadioAstron detections of the classic quasar 3C273 on interferometric baselines up to 171,000 km suggest brightness temperatures exceeding expected limits from the "inverse-Compton catastrophe" by two orders of magnitude. We show that at 18 cm, these estimates most probably arise from refractive substructure introduced by scattering in the interstellar medium. We use the scattering properties to estimate an intrinsic brightness temperature of 7*10^12 K, which is consistent with expected theoretical limits, but which is ~15 times lower than estimates that neglect substructure. At 6.2 cm, the substructure influences the measured values appreciably but gives an estimated brightness temperature that is comparable to models that do not account for the substructure. At 1.35 cm, the substructure does not affect the extremely high inferred brightness temperatures, in excess of 10^13 K. We also demonstrate that for a source having a Gaussian surface brightness profile, a single long-baseline estimate of refractive substructure determines an absolute minimum brightness temperature, if the scattering properties along a given line of sight are known, and that this minimum accurately approximates the apparent brightness temperature over a wide range of total flux densities.

At radio wavelengths, scattering in the interstellar medium distorts the appearance of astronomical sources. Averaged over a scattering ensemble, the result is a blurred image of the source. However, Narayan & Goodman (1989) and Goodman & Narayan (1989) showed that for an incomplete average, scattering introduces refractive substructure in the image of a point source that is both persistent and wideband. We show that this substructure is quenched but not smoothed by an extended source. As a result, when the scatter-broadening is comparable to or exceeds the unscattered source size, the scattering can introduce spurious compact features into images. In addition, we derive efficient strategies to numerically compute realistic scattered images, and we present characteristic examples from simulations. Our results show that refractive substructure is an important consideration for ongoing missions at the highest angular resolutions, and we discuss specific implications for RadioAstron and the Event Horizon Telescope.

Almost 50 years after radio pulsars were discovered in 1967, our understanding of these objects remains incomplete. On the one hand, within a few years it became clear that neutron star rotation gives rise to the extremely stable sequence of radio pulses, that the kinetic energy of rotation provides the reservoir of energy, and that electromagnetic fields are the braking mechanism. On the other hand, no consensus regarding the mechanism of coherent radio emission or the conversion of electromagnetic energy to particle energy yet exists. In this review, we report on three aspects of pulsar structure that have seen recent progress: the self-consistent theory of the magnetosphere of an oblique magnetic rotator; the location, geometry, and optics of radio emission; and evolution of the angle between spin and magnetic axes. These allow us to take the next step in understanding the physical nature of the pulsar activity.

The SKA will discover tens of thousands of pulsars and provide unprecedented data quality on these, as well as the currently known population, due to its unrivalled sensitivity. Here, we outline the state of the art of our understanding of magnetospheric radio emission from pulsars and how we will use the SKA to solve the open problems in pulsar magnetospheric physics.

The image of the emission surrounding the black hole in the center of the Milky Way is predicted to exhibit the imprint of general relativistic (GR) effects, including the existence of a shadow feature and a photon ring of diameter ~50 microarcseconds. Structure on these scales can be resolved by millimeter-wavelength very long baseline interferometry (VLBI). However, strong-field GR features of interest will be blurred at lambda >= 1.3 mm due to scattering by interstellar electrons. The scattering properties are well understood over most of the relevant range of baseline lengths, suggesting that the scattering may be (mostly) invertible. We simulate observations of a model image of Sgr A* and demonstrate that the effects of scattering can indeed be mitigated by correcting the visibilities before reconstructing the image. This technique is also applicable to Sgr A* at longer wavelengths.

We have detected substructure within the smooth scattering disk of the celebrated Galactic Center radio source Sagittarius A* (SgrA*). We observed this structure at 1.3 cm wavelength with the Very Long Baseline Array together with the Green Bank Telescope, on baselines of up to 3000 km, long enough to completely resolve the average scattering disk. Such structure is predicted theoretically, as a consequence of refraction by large-scale plasma fluctuations in the interstellar medium. Along with the much-studied $\theta_\mathrm{d}\propto \lambda^2$ scaling of angular broadening $\theta_\mathrm{d}$ with observing wavelength $\lambda$, our observations indicate that the spectrum of interstellar turbulence is shallow, with an inner scale larger than 300 km. The substructure is consistent with an intrinsic size of about 1 mas at 1.3 cm wavelength, as inferred from deconvolution of the average scattering. Further observations of the substructure can set stronger constraints on the properties of scattering material and on the intrinsic size of SgrA*. These constraints will guide understanding of effects of scatter-broadening and emission physics of the black hole, in images with the Event Horizon Telescope at millimeter wavelengths.

RadioAstron space-ground VLBI observations of the pulsar B0950+08, conducted with the 10-m space radio telescope in conjunction with the Arecibo 300-m telescope and Westerbork Synthesis Radio Telescope at a frequency of 324 MHz, were analyzed in order to investigate plasma inhomogeneities in the direction of this nearby pulsar. The observations were conducted at a spacecraft distance of 330,000 km, resulting in a projected baseline of 220,000 km, providing the greatest angular resolution ever achieved at meter wavelengths. Our analysis is based on fundamental behavior of structure and coherence functions. We find that the pulsar shows scintillation on two frequency scales, both much less than the observing frequency; but modulation is less than 100%. We infer that the scattering is weak, but a refracting wedge disperses the scintillation pattern. The refraction angle of this "cosmic prism" is measured as theta_0=1.1 - 4.4 mas, with the refraction direction being approximately perpendicular to the observer velocity. We show that the observed parameters of scintillation effects indicate that two plasma layers lie along the line of sight to the pulsar, at distances of 4.4 - 16.4 pc and 26 - 170 pc, and traveling in different directions relative to the line of sight. Spectra of turbulence for the two layers are found to follow a power law with the indices gamma_1 = gamma_2 = 3.00 +/- 0.08, significantly different from the index expected for a Kolmogorov spectrum of turbulence, gamma=11/3.

We observe an extraordinary link in the Crab pulsar between the enhancement of an optical pulse and the timing of the corresponding giant radio pulse. At optical through infrared wavelengths, our observations use the high time resolution of ARCONS, a unique superconducting energy-resolving photon-counting array at the Palomar 200-inch telescope. At radio wavelengths, we observe with the Robert C. Byrd Green Bank Telescope and the GUPPI backend. We see an $11.3\pm2.5\%$ increase in peak optical flux for pulses that have an accompanying giant radio pulse arriving near the peak of the optical main pulse, in contrast to a $3.2\pm0.5\%$ increase when an accompanying giant radio pulse arrives soon after the optical peak. We also observe that the peak of the optical main pulse is $2.8\pm0.8\%$ enhanced when there is a giant radio pulse accompanying the optical interpulse. We observe no statistically significant spectral differences between optical pulses accompanied by and not accompanied by giant radio pulses. Our results extend previous observations of optical-radio correlation to the time and spectral domains. Our refined temporal correlation suggests that optical and radio emission are indeed causally linked, and the lack of spectral differences suggests that the same mechanism is responsible for all optical emission.

We derive the distribution of interferometric visibility for a source exhibiting strong diffractive scintillation, with particular attention to spectral resolution at or near the Nyquist limit. We also account for arbitrary temporal averaging, intrinsic variability within the averaging time, and the possibility of spatially-extended source emission. We demonstrate that the interplay between scintillation and self-noise induces several remarkable features, such as a broad "skirt" in the visibility distribution. Our results facilitate the interpretation of interferometric observations of pulsars at meter and decimeter wavelengths.

Using a maximum-likelihood criterion, we derive optimal correlation strategies for signals with and without digitization. We assume that the signals are drawn from zero-mean Gaussian distributions, as is expected in radio-astronomical applications, and we present correlation estimators both with and without a priori knowledge of the signal variances. We demonstrate that traditional estimators of correlation, which rely on averaging products, exhibit large and paradoxical noise when the correlation is strong. However, we also show that these estimators are fully optimal in the limit of vanishing correlation. We calculate the bias and noise in each of these estimators and discuss their suitability for implementation in modern digital correlators.

We compare the noise in interferometric measurements of the Vela pulsar from ground- and space-based antennas with theoretical predictions. The noise depends on both the flux density and the interferometric phase of the source. Because the Vela pulsar is bright and scintillating, these comparisons extend into both the low and high signal-to-noise regimes. Furthermore, our diversity of baselines explores the full range of variation in interferometric phase. We find excellent agreement between theoretical expectations and our estimates of noise among samples within the characteristic scintillation scales. Namely, the noise is drawn from an elliptical Gaussian distribution in the complex plane, centered on the signal. The major axis, aligned with the signal phase, varies quadratically with the signal, while the minor axis, at quadrature, varies with the same linear coefficients. For weak signal, the noise approaches a circular Gaussian distribution. Both the variance and covariance of the noise are also affected by artifacts of digitization and correlation. In particular, we show that gating introduces correlations between nearby spectral channels.

We present measurements of the linear diameter of the emission region of the Vela pulsar at observing wavelength lambda=18 cm. We infer the diameter as a function of pulse phase from the distribution of visibility on the Mopra-Tidbinbilla baseline. As we demonstrate, in the presence of strong scintillation, finite size of the emission region produces a characteristic W-shaped signature in the projection of the visibility distribution onto the real axis. This modification involves heightened probability density near the mean amplitude, decreased probability to either side, and a return to the zero-size distribution beyond. We observe this signature with high statistical significance, as compared with the best-fitting zero-size model, in many regions of pulse phase. We find that the equivalent full width at half maximum of the pulsar's emission region decreases from more than 400 km early in the pulse to near zero at the peak of the pulse, and then increases again to approximately 800 km near the trailing edge. We discuss possible systematic effects, and compare our work with previous results.

Using a novel technique, we achieve ~100 picoarcsecond resolution and set an upper bound of less than 4 km for the characteristic size of the Vela pulsar's emission region. Specifically, we analyze flux-density statistics of the Vela pulsar at 760 MHz. Because the pulsar exhibits strong diffractive scintillation, these statistics convey information about the spatial extent of the radio emission region. We measure both a characteristic size of the emission region and the emission sizes for individual pulses. Our results imply that the radio emission altitude for the Vela pulsar at this frequency is less than 340 km.

We derive the distribution of flux density of a compact source exhibiting strong diffractive scintillation. Our treatment accounts for arbitrary spectral averaging, spatially-extended source emission, and the possibility of intrinsic variability within the averaging time, as is typical for pulsars. We also derive the modulation index and present a technique for estimating the self-noise of the distribution, which can be used to identify amplitude variations on timescales shorter than the spectral accumulation time. Our results enable a for direct comparison with ultra-high resolution observations of pulsars, particularly single-pulse studies with Nyquist-limited resolution, and can be used to identify the spatial emission structure of individual pulses at a small fraction of the diffractive scale.

We compare signal and noise for observations of the scintillating pulsar B0834+06, using very-long baseline interferometry and a single-dish spectrometer. Comparisons between instruments and with models suggest that amplitude variations of the pulsar strongly affect the amount and distribution of self-noise. We show that noise follows a quadratic polynomial with flux density, in spectral observations. Constant coefficients, indicative of background noise, agree well with expectation; whereas second-order coefficients, indicative of self-noise, are about 3 times values expected for a pulsar with constant on-pulse flux density. We show that variations in flux density during the 10-sec integration account for the discrepancy. In the secondary spectrum, about 97% of spectral power lies within the pulsar's typical scintillation bandwidth and timescale; an extended scintillation arc contains about 3%. For a pulsar with constant on-pulse flux density, noise in the dynamic spectrum will appear as a uniformly-distributed background in the secondary spectrum. We find that this uniform noise background contains 95% of noise in the dynamic spectrum for interferometric observations; but only 35% of noise in the dynamic spectrum for single-dish observations. Receiver and sky dominate noise for our interferometric observations, whereas self-noise dominates for single-dish. We suggest that intermittent emission by the pulsar, on timescales < 300 microseconds, concentrates self-noise near the origin in the secondary spectrum, by correlating noise over the dynamic spectrum. We suggest that intermittency sets fundamental limits on pulsar astrometry or timing. Accounting of noise may provide means for detection of intermittent sources, when effects of propagation are unknown or impractical to invert.

We show that intermittency of noiselike emission, after propagation through a scattering medium, affects the distribution of noise in the observed correlation function. Intermittency also affects correlation of noise among channels of the spectrum, but leaves the average spectrum, average correlation function, and distribution of noise among channels of the spectrum unchanged. Pulsars are examples of such sources: intermittent and affected by interstellar propagation. We assume that the source emits Gaussian white noise, modulated by a time-envelope. Propagation convolves the resulting time series with an impulse-response function that represents effects of dispersion, scattering, and absorption. We assume that this propagation kernel is shorter than the time for an observer to accumulate a single spectrum. We show that rapidly-varying intermittent emission tends to concentrate noise near the central lag of the correlation function. We derive mathematical expressions for this effect, in terms of the time envelope and the propagation kernel. We present examples, discuss effects of background noise, and compare our results with observations.

We report the results of simultaneous observations of the Vela pulsar in X-rays and radio from the RXTE satellite and the Mount Pleasant Radio Observatory in Tasmania. We sought correlations between the Vela's X-ray and radio flux densities and radio arrival times on a pulse by pulse basis. We found significantly higher flux density in Vela's main X-ray peak during radio pulses that arrived early. This excess flux shifts to the 'trough' following the 2nd X-ray peak during radio pulses that arrive later. We suggest that the mechanism producing the radio pulses is intimately connected to the mechanism producing X-rays. Current models using resonant absorption in the outer magnetosphere as a cause of the radio emission, and less directly of the X-ray emission, are explored as a possible explanation for the correlation.

We report the results of simultaneous observations of the Vela pulsar in X-rays and radio from the RXTE satellite and the Mount Pleasant Radio Observatory in Tasmania. We sought correlations between the Vela's X-ray emission and radio arrival times on a pulse by pulse basis. At a confidence level of 99.8% we have found significantly higher flux density in Vela's main X-ray peak during radio pulses that arrived early. This excess flux shifts to the 'trough' following the 2nd X-ray peak during radio pulses that arrive later. Our results suggest that the mechanism producing the radio pulses is intimately connected to the mechanism producing X-rays. Current models using resonant absorption of radio emission in the outer magnetosphere as a cause of the X-ray emission are explored as a possible explanation for the correlation.

Topics of the Colloquium: a) Interplanetary scintillation b) Interstellar scintillation c) Modeling and physical origin of the interplanetary and the interstellar plasma turbulence d) Scintillation as a tool for investigation of radio sources e) Seeing through interplanetary and interstellar turbulent media Ppt-presentations are available on the Web-site: http://www.prao.ru/conf/Colloquium/main.html

We sought to determine the scale of scintillation in the interstellar plasma of PSR J0437-4715. We used the Very Long Baseline Array to obtain scintillation amplitude and phase data, from dynamic spectra at 327 MHz. We observe two scales of scintillation of pulsar PSR J0437-4715, differing by more than an order of magnitude in scintillation bandwidth. The wider-bandwidth scale of scintillation that we observe indicates less scattering for this pulsar than for other nearby pulsars, except for PSR B0950+08.

We studied the turbulence spectrum of the local interstellar plasma in the direction of PSR J0437-4715, on the basis of our observations and those reported earlier by others. We combine these data to form a structure function for the variations of phase along the line of sight to the pulsar. For observations that did not report them, we infer modulation indices from a theoretical model. We find that all of the observations fit a power-law spectrum of turbulence with index n=3.46+/-0.20. We suggest that differences among reported values for scintillation bandwidth and timescale for this pulsar arise from differences in observing parameters. We suggest that refractive effects dominate for this line of sight, with refraction angle about twice the diffraction angle at 330 MHz observing frequency. We suggest that the scattering of this pulsar lies in a layer of enhanced turbulence, about 10 pc from the Sun. We propose that the flux variations of the extragalactic source PKS 0405-385 arise in the same scattering layer.

I calculate the noise in the measured correlation functions and spectra of digitized, noiselike signals. In the spectral domain, the signals are drawn from a Gaussian distribution with variance that depends on frequency. Nearly all astrophysical signals have noiselike statistics of this type, many with important spectral variations. The noise is digitized and correlated to form a measured correlation function, which is then Fourier transformed to produce a measured spectrum. When averaged over many samples, the elements of the correlation function and of the spectrum, follow Gaussian distributions. For each element, the mean of that distribution is the deterministic part of the measurement. The standard deviation of the Gaussian is the noise. Here I calculate that noise, as a function of the parameters of digitization. The noise of the correlation function is related to the underlying spectrum, by constants that depend on the digitization parameters. Noise affects variances of elements of the correlation function and covariances between them. In the spectral domain, noise also produces variances and covariances. I show that noise is correlated between spectral channels, for digitized spectra, and calculate the correlation. These statistics of noise are important for understanding of signals sampled with very high signal-to-noise ratio, or signals with rapidly-changing levels such as pulsars.

I calculate the statistics of correlation of two digitized noiselike signals, which are drawn from complex Gaussian distributions, sampled, quantized, correlated, and averaged. Averaged over many such samples, the correlation r approaches a Gaussian distribution. The mean and variance of r fully characterize the distribution of r. The mean corresponds to the reproducible part of the measurement, and the variance corresponds to the random part, or noise. I investigate the case of nonnegligible covariance rho between the signals. Noise in the correlation can increase or decrease, depending on quantizer parameters, when rho increases. This contrasts with the correlation of continuously valued or unquantized signals, for which the noise in phase with rho increases with increasing rho, and noise out of phase decreases. Indeed, for some quantizer parameters, I find that the correlation of quantized signals provides a more accurate estimate of rho than would correlation without quantization. I present analytic results in exact form and as polynomial expansions, and compare these mathematical results with results of computer simulations.

Radio waves propagating from distant pulsars in the interstellar medium (ISM), are refracted by electron density inhomogeneities, so that the intensity of observed pulses fluctuates with time. The theory relating the observed pulse time-shapes to the electron-density correlation function has developed for 30 years, however, two puzzles have remained. First, observational scaling of pulse broadening with the pulsar distance is anomalously strong; it is consistent with the standard model only when non-uniform statistics of electron fluctuations along the line of sight are assumed. Second, the observed pulse shapes are consistent with the standard model only when the scattering material is concentrated in a narrow slab between the pulsar and the Earth. We propose that both paradoxes are resolved at once if one assumes stationary and uniform, but non-Gaussian statistics of the electron-density distribution. Such statistics must be of Levy type, and the propagating ray should exhibit a Levy flight. We propose that a natural realization of such statistics may be provided by the interstellar medium with random electron-density discontinuities. We develop a theory of wave propagation in such a non-Gaussian random medium, and demonstrate its good agreement with observations. The qualitative introduction of the approach and the resolution of the anomalous-scaling paradox was presented earlier in [PRL 91, 131101 (2003); ApJ 584, 791 (2003)].

Observations of radio signals from distant pulsars provide a valuable tool for investigation of interstellar turbulence. The time-shapes of the signals are the result of pulse broadening by the fluctuating electron density in the interstellar medium. While the scaling of the shapes with the signal frequency is well understood, the observed anomalous scaling with respect to the pulsar distance has remained a puzzle for more than 30 years. We propose a new model for interstellar electron density fluctuations, which explains the observed scaling relations. We suggest that these fluctuations obey Levy statistics rather than Gaussian statistics, as assumed in previous treatments of interstellar scintillations.

We report on our analysis of a 274 ks observation of the Vela pulsar with the Rossi X-Ray Timing Explorer (RXTE). The double-peaked, pulsed emission at 2 - 30 keV, which we had previously detected during a 93 ks observation, is confirmed with much improved statistics. There is now clear evidence, both in the spectrum and the light curve, that the emission in the RXTE band is a blend of two separate non-thermal components. The spectrum of the harder component connects smoothly with the OSSE, COMPTEL and EGRET spectrum and the peaks in the light curve are in phase coincidence with those of the high-energy light curve. The spectrum of the softer component is consistent with an extrapolation to the pulsed optical flux, and the second RXTE pulse is in phase coincidence with the second optical peak. In addition, we see a peak in the 2-8 keV RXTE pulse profile at the radio phase.

Temporal broadening of pulsar signals results from electron density fluctuations in the interstellar medium that cause the radiation to travel along paths of different lengths. The Gaussian theory of fluctuations predicts that the pulse temporal broadening should scale with the wavelength as lambda^4, and with the dispersion measure (corresponding to distance to the pulsar) as DM^2. For large dispersion measure, DM > 20 pc/cm^3, the observed scaling is lambda^4 DM^4, contradicting the conventional theory. Although the problem has existed for 30 years, there has been no resolution to this paradox. We suggest that scintillations for distant pulsars are caused by non-Gaussian, spatially intermittent density fluctuations with a power-like probability distribution. This probability distribution does not have a second moment in a large range of density fluctuations, and therefore the previously applied conventional Fokker-Planck theory does not hold. Instead, we propose to apply the theory of Levy distributions (so-called Levy flights). Using the scaling analysis (confirmed by numerical simulations of ray propagation) we show that the observed scaling is recovered for large DM, if the density differences, delta N, have Levy distribution decaying as |delta N|^{-5/3}.

I suggest that radio-wave scattering by the interstellar plasma, in combination with subsonic gradients in the Doppler velocity of interstellar HI, is responsible for the observed small-scale variation in HI absorption spectra of pulsars. Velocity gradients on the order of 0.05 to 0.3 km/s across 1 AU can produce the observed variations. I suggest observational tests to distinguish between this model and the traditional picture of small-scale opacity variations from cloudlets.

In the relativistic plasma surrounding a pulsar, a subluminal ordinary-mode electromagnetic wave will propagate along a magnetic field line. After some distance, it can break free of the field line and escape the magnetosphere to reach an observer. We describe a simple model of pulsar radio emission based on this scenario and find that applying this model to the case of the Vela pulsar reproduces qualitative characteristics of the observed Vela pulse profile.

We present spectroscopic results of 146 water maser outbursts in W49 N, obtained with the Metsahovi radio telescope at 22 GHz. Combining this data with Gwinn's VLBI results, we were able to fix the free parameters in the shock model of Hollenbach and McKee and the maser model of Elitzur, Hollenbach, and McKee. This enabled a straightforward determination of some 20 shock and maser parameters, among others, the kinetic temperature, postshock density, water density, water abundance, preshock and postshock magnetic field strengths as well as the inclination of the mean field with respect to the line of sight. A step-by-step presentation of our diagnostic method is given and the relation between observations and model parameters is discussed. One uniquely powerful outburst feature, referred to as the "big flare feature", showed also the narrowest linewidth. Observations indicate that the velocity of this feature lies in the plane of the sky, whereas preshock and postshock magnetic fields are directed nearly along the line of sight. Consequently, Alfvenic wave fluctuations along the line of sight, and linewidth, are minimal, and a very high aspect ratio is achieved. Furthermore, the big flare feature stands out through its low space velocity, higher temperature (480 K), and larger preshock magnetic field strength (8.2 mG).

We describe measurements of the size of the Vela pulsar via scintillation, using both fits to the distribution of intensity and measurements of the modulation index. We briefly discuss systematic effects other than source size that can affect the distribution, including gain variations, self-noise, scintillation shot noise, and correlator saturation. Modulation index, a single number, can be biased by all of these, whereas the distribution of intensity is affected in different ways by different effects, providing means of distinguishing among them. Self-noise and gain variations are likely more important at long observing wavelengths, and correlator saturation and scintillation shot noise at short wavelengths. We find a size of about 500 km at decimeter wavelengths. Interestingly, this agrees with measurements of modulation index by Roberts & Ables at the same wavelength. Their results (and more recently that reported by Macquart et al.) suggest that size decreases with increasing wavelength. Although consistent with the observations, this conclusion is perhaps surprising from the standpoint of the traditional radius-to-frequency mapping. However, these measurements are of size rather than height; and of course systematic effects may play a role.

We report on our analysis of a 300 ks observation of the Vela pulsar with the Rossi X-Ray Timing Explorer (RXTE). The double-peaked, pulsed emission at 2 - 30 keV, which we had previously detected during a 93 ks observation, is confirmed with much improved statistics. There is now clear evidence, both in the spectrum and the light curve, that the emission in the RXTE band is a blend of two separate components. The spectrum of the harder component connects smoothly with the OSSE, COMPTEL and EGRET spectrum and the peaks in the light curve are in phase coincidence with those of the high-energy light curve. The spectrum of the softer component is consistent with an extrapolation to the pulsed optical flux, and the second RXTE pulse is in phase coincidence with the second optical peak. In addition, we see a peak in the 2-8 keV RXTE pulse profile at the radio phase.

We present measurements of the size of the Vela pulsar in 3 gates across the pulse, from observations of the distribution of intensity. We calculate the effects on this distribution of noise in the observing system, and measure and remove it using observations of a strong continuum source. We also calculate and remove the expected effects of averaging in time and frequency. We find that effects of variations in pulsar flux density and instrumental gain, self-noise, and one-bit digitization are undetectably small. Effects of normalization of the correlation are detectable, but do not affect the fitted size. The size of the pulsar declines from 440 +/- 90 km (FWHM of best-fitting Gaussian distribution) to less than 200 km across the pulse. We discuss implications of this size for theories of pulsar emission.

A consistent approach for an exhaustive solution of the problem of propagation of light rays in the field of gravitational waves emitted by a localized source of gravitational radiation is developed in the first post-Minkowskian and quadrupole approximation of General Relativity. We demonstrate that the equations of light propagation in the retarded gravitational field of an arbitrary localized source emitting quadrupolar gravitational waves can be integrated exactly. The influence of the gravitational field on the light propagation is examined not only in the wave zone but also in cases when light passes through the intermediate and near zones of the source. Explicit analytic expressions for light deflection and integrated time delay (Shapiro effect) are obtained accounting for all possible retardation effects and arbitrary relative locations of the source of gravitational waves, that of light rays, and the observer. It is shown that the ADM and harmonic gauge conditions can both be satisfied simultaneously outside the source of gravitational waves. Their use drastically simplifies the integration of light propagation equations and those for the motion of light source and observer in the field of the source of gravitational waves, leading to the unique interpretation of observable effects. The two limiting cases of small and large values of impact parameter are elaborated in more detail. Explicit expressions for Shapiro effect and deflection angle are obtained in terms of the transverse-traceless part of the space-space components of the metric tensor. We also discuss the relevance of the developed formalism for interpretation of radio interferometric and timing observations, as well as for data processing algorithms for future gravitational wave detectors.

The BN/KL region in the Orion molecular cloud is an archetype in the study of the formation of stars much more massive than the Sun. This region contains luminous young stars and protostars, but it is difficult to study because of overlying dust and gas. Our basic expectations are shaped to some extent by the present theoretical picture of star formation, the cornerstone of which is that protostars acrete gas from rotating equatorial disks, and shed angular momentum by ejecting gas in bipolar outflows. The main source of the outflow in the BN/KL region may be an object known as radio source I, which is commonly believed to be surrounded by a rotating disk of molecular material. Here we report high-resolution observations of silicon monoxide (SiO) and water maser emission from the gas surrounding source I; we show that within 60 AU (about the size of the Solar System), the region is dominated by a conical bipolar outflow, rather than the expected disk. A slower outflow, close to the equatorial plane of the protostellar system, extends to radii of 1,000 AU.

We present a parsec-scale image of the OH maser in the nucleus of the active galaxy IIIZw35, made using the Very Long Baseline Array at a wavelength of 18 cm. We detected two distinct components, with a projected separation of 50 pc (for D=110 Mpc) and a separation in Doppler velocity of 70 km/s, which contain 50% of the total maser flux. Velocity gradients within these components could indicate rotation of clouds with binding mass densities of ~7000 solar masses per cubic parsec, or total masses of more than 500,000 solar masses. Emission in the 1665-MHz OH line is roughly coincident in position with that in the 1667-MHz line, although the lines peak at different Doppler velocities. We detected no 18 cm continuum emission; our upper limit implies a peak apparent optical depth greater than 3.4, assuming the maser is an unsaturated amplifier of continuum radiation.

We report observational upper limits on the mass-energy of the cosmological gravitational-wave background, from limits on proper motions of quasars. Gravitational waves with periods longer than the time span of observations produce a simple pattern of apparent proper motions over the sky, composed primarily of second-order transverse vector spherical harmonics. A fit of such harmonics to measured motions yields a 95%-confidence limit on the mass-energy of gravitational waves with frequencies <2e-9 Hz, of <0.11/h*h times the closure density of the universe.

We conducted a very long baseline interferometric observation of 5 nearby pulsars and did not resolve the scattering disks of any of these sources. Using our upper limits on the angular diameters of these scattering disks and published values of the broadening times and proper motion velocities, we constrain the possible distributions of scattering material. The material responsible for scattering these sources is neither uniformly distributed nor concentrated at the surface of the Local Bubble. We argue that these pulsars themselves influence their environments to produce this scattering material.

We have made the first VLBI synthesis images of the H2O maser emission associated with the central engine of the Seyfert galaxy NGC 1068. Emission extends about +/-300 km/s from the systemic velocity. Images with submilliarcsecond angular resolution show that the red-shifted emission lies along an arc to the northwest of the systemic emission. (The blue-shifted emission has not yet been imaged with VLBI.) Based on the maser velocities and the relative orientation of the known radio jet, we propose that the maser emission arises on the surface of a nearly edge-on torus, where physical conditions are conducive to maser action. The visible part of the torus is axially thick, with comparable height and radius. The velocity field indicates sub-Keplerian differential rotation around a central mass of about 1e7 Msun that lies within a cylindrical radius of about 0.65 pc. The estimated luminosity of the central engine is about 0.5 of the Eddington limit. There is no detectable compact radio continuum emission near the proposed center of the torus (T_B< 5e6 K on size scales of about 0.1 pc), so that the observed flat-spectrum core cannot be direct self-absorbed synchrotron radiation.

We use interstellar scattering of the Vela pulsar to determine the size of its emission region. From interferometric phase variations on short baselines, we find that radio-wave scattering broadens the source by 3.4+/-0.3 milliarcseconds along the major axis at position angle 81+/-3 degrees. The ratio of minor axis to major axis is 0.51+/-0.03. Comparison of angular and temporal broadening indicates that the scattering material lies in the Vela-X supernova remnant surrounding the pulsar. From the modulation of the pulsar's scintillation on very short baselines, we infer a size of 500 km for the pulsar's emission region. We suggest that radio-wave refraction within the pulsar's magnetosphere may plausibly explain this size.

Gravitational waves affect the observed direction of light from distant sources. At telescopes, this change in direction appears as periodic variations in the apparent positions of these sources on the sky; that is, as proper motion. A wave of a given phase, traveling in a given direction, produces a characteristic pattern of proper motions over the sky. Comparison of observed proper motions with this pattern serves to test for the presence of gravitational waves. A stochastic background of waves induces apparent proper motions with specific statistical properties, and so, may also be sought. In this paper we consider the effects of a cosmological background of gravitational radiation on astrometric observations. We derive an equation for the time delay measured by two antennae observing the same source in an Einstein-de Sitter spacetime containing gravitational radiation. We also show how to obtain similar expressions for curved Friedmann-Robertson-Walker spacetimes.