The Crab pulsar has striking radio emission properties, with the two dominant pulse components -- the main pulse and the interpulse -- consisting entirely of giant pulses. The emission is scattered in both the Crab nebula and the interstellar medium, causing multi-path propagation and thus scintillation. We study the scintillation of the Crab's giant pulses using phased Westerbork data at 1668\,MHz. We find that giant pulse spectra correlate at only $\sim 2 \%$, much lower than the $1/3$ correlation expected from a randomized signal imparted with the same impulse response function. In addition, we find that the main pulse and the interpulse appear to scintillate differently, with the scintillation of the interpulse correlating to the main pulse at a lower amplitude over a longer time and frequency range. These lines of evidence suggest that the giant pulse emission regions are extended, and that the main pulse and interpulse arise in physically distinct regions which are resolved by the scattering screen. Assuming the scattering takes place in the nebular filaments, the emission regions are of order a light cylinder radius, as projected on the sky. With further VLBI and multi-frequency data, it may be possible to measure the distance to the scattering screens, the size of giant pulse emission regions, and the physical separation between the pulse components.

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.

Characterizing stellar convection in multiple dimensions is a topic at the forefront of stellar astrophysics. Numerical simulations are an essential tool for this task. We present an extension of the existing numerical tool-kit A-MaZe that enables such simulations of stratified flows in a gravitational field. The finite-volume based, cell-centered, and time-explicit hydrodynamics solver of A-MaZe was extended such that the scheme is now well-balanced in both momentum and energy. The algorithm maintains an initially static balance between gravity and pressure to machine precision. Quasi-stationary convection in slab-geometry preserves gas energy (internal plus kinetic) on average despite strong local up- and down-drafts. By contrast, a more standard numerical scheme is demonstrated to result in substantial gains of energy within a short time on purely numerical grounds. The test is further used to point out the role of dimensionality, viscosity, and Rayleigh number for compressible convection. Applications to a young sun in 2D and 3D, covering a part of the inner radiative zone as well as the outer convective zone, demonstrate that the scheme meets its initial design goal. Comparison with results obtained for a physically identical setup with a time-implicit code show qualitative agreement.

The aim of our work was to study the spatial structure of inhomogeneities of interstellar plasma in the directions of five pulsars: B0823+26, B0834+06, B1237+25, B1929+10, and B2016+28. Observations of these pulsars were made with RadioAstron space-ground radio interferometer at 324 MHz. We measured the angular size of the scattering disks to be in range between 0.63 and 3.2 mas. We determined the position of scattering screens on the line of sight. Independent estimates of the distances to the screens were made from the curvature of parabolic arcs revealed in the secondary spectra of four pulsars. The model of uniform distribution of inhomogeneities on the line of sight is not suitable. According to the results, we came to the conclusion that scattering is mainly produced by compact plasma layers and the uniform model of inhomogeneties distribution on the line of sight in not applicable.

We study lithium depletion in low-mass and solar-like stars as a function of time, using a new diffusion coefficient describing extra-mixing taking place at the bottom of a convective envelope. This new form is motivated by multi-dimensional fully compressible, time implicit hydrodynamic simulations performed with the MUSIC code. Intermittent convective mixing at the convective boundary in a star can be modeled using extreme value theory, a statistical analysis frequently used for finance, meteorology, and environmental science. In this letter, we implement this statistical diffusion coefficient in a one-dimensional stellar evolution code, using parameters calibrated from multi-dimensional hydrodynamic simulations of a young low-mass star. We propose a new scenario that can explain observations of the surface abundance of lithium in the Sun and in clusters covering a wide range of ages, from $\sim$ 50 Myr to $\sim$ 4 Gyr. Because it relies on our physical model of convective penetration, this scenario has a limited number of assumptions. It can explain the observed trend between rotation and depletion, based on a single additional assumption, namely that rotation affects the mixing efficiency at the convective boundary. We suggest the existence of a threshold in stellar rotation rate above which rotation strongly prevents the vertical penetration of plumes and below which rotation has small effects. In addition to providing a possible explanation for the long standing problem of lithium depletion in pre-main sequence and main sequence stars, the strength of our scenario is that its basic assumptions can be tested by future hydrodynamic simulations.

We examine a penetration layer formed between a central radiative zone and a large convection zone in the deep interior of a young low-mass star. Using the Multidimensional Stellar Implicit Code (MUSIC) to simulate two-dimensional compressible stellar convection in a spherical geometry over long times, we produce statistics that characterize the extent and impact of convective penetration in this layer. We apply extreme value theory to the maximal extent of convective penetration at any time. We compare statistical results from simulations which treat non-local convection, throughout a large portion of the stellar radius, with simulations designed to treat local convection in a small region surrounding the penetration layer. For each of these situations, we compare simulations of different resolution, which have different velocity magnitudes. We also compare statistical results between simulations that radiate energy at a constant rate to those that allow energy to radiate from the stellar surface according to the local surface temperature. Based on the frequency and depth of penetrating convective structures, we observe two distinct layers that form between the convection zone and the stable radiative zone. We show that the probability density function of the maximal depth of convective penetration at any time corresponds closely in space with the radial position where internal waves are excited. We find that the maximal penetration depth can be modeled by a Weibull distribution with a small shape parameter. Using these results, and building on established scalings for diffusion enhanced by large-scale convective motions, we propose a new form for the diffusion coefficient that may be used for one-dimensional stellar evolution calculations in the large P\'eclet number regime. These results should contribute to the 321D link.

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 present the results of a numerical benchmark study for the MUlti-dimensional Stellar Implicit Code (MUSIC) based on widely applicable two- and three-dimensional compressible hydrodynamics problems relevant to stellar interiors. MUSIC is an implicit large eddy simulation code that uses implicit time integration, implemented as a Jacobian-free Newton Krylov method. A physics based preconditioning technique which can be adjusted to target varying physics is used to improve the performance of the solver. The problems used for this benchmark study include the Rayleigh-Taylor and Kelvin-Helmholtz instabilities, and the decay of the Taylor-Green vortex. Additionally we show a test of hydrostatic equilibrium, in a stellar environment which is dominated by radiative effects. In this setting the flexibility of the preconditioning technique is demonstrated. This work aims to bridge the gap between the hydrodynamic test problems typically used during development of numerical methods and the complex flows of stellar interiors. A series of multi-dimensional tests are performed and analysed. Each of these test cases is analysed with a simple, scalar diagnostic, with the aim of enabling direct code comparisons. As the tests performed do not have analytic solutions we verify MUSIC by comparing to established codes including ATHENA and the PENCIL code. MUSIC is able to both reproduce behaviour from established and widely-used codes as well as results expected from theoretical predictions. This benchmarking study concludes a series of papers describing the development of the MUSIC code and provides confidence in the future applications.

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.

Context: We study the impact of two-dimensional spherical shells on compressible convection. Realistic profiles for density and temperature from a one-dimensional stellar evolution code are used to produce a model of a large stellar convection zone representative of a young low-mass star. Methods: We perform hydrodynamic implicit large-eddy simulations of compressible convection using the MUltidimensional Stellar Implicit Code (MUSIC). Because MUSIC has been designed to use realistic stellar models produced from one-dimensional stellar evolution calculations, MUSIC simulations are capable of seamlessly modeling a whole star. Simulations in two-dimensional spherical shells that have different radial extents are performed over hundreds of convective turnover times, permitting the collection of well-converged statistics. Results: We evaluate basic statistics of the convective turnover time, the convective velocity, and the overshooting layer. These quantities are selected for their relevance to one-dimensional stellar evolution calculations, so that our results are focused toward the 321D link. The inclusion in the spherical shell of the boundary between the radiative and convection zones decreases the amplitude of convective velocities in the convection zone. The inclusion of near-surface layers in the spherical shell can increase the amplitude of convective velocities, although the radial structure of the velocity profile established by deep convection is unchanged. The impact from including the near-surface layers depends on the speed and structure of small-scale convection in the near-surface layers. Larger convective velocities in the convection zone result in a commensurate increase in the overshooting layer width and decrease in the convective turnover time. These results provide support for non-local aspects of convection.

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.

This work is the first attempt to describe the multi-dimensional structure of accreting young stars based on fully compressible time implicit multi-dimensional hydrodynamics simulations. One major motivation is to analyse the validity of accretion treatment used in previous 1D stellar evolution studies. We analyse the effect of accretion on the structure of a realistic stellar model of the young Sun. Our work is inspired by the numerical work of Kley \& Lin (1996, ApJ, 461, 933) devoted to the structure of the boundary layer in accretion disks. We analyse the redistribution of accreted material with a range of values of specific entropy relative to the bulk specific entropy of the material in the accreting object's convective envelope. A primary goal is to understand whether and how accreted energy deposited onto a stellar surface is redistributed in the interior. This study focusses on the high accretion rates characteristic of FU Ori systems. We find that the highest entropy cases produce a distinctive behaviour in the mass redistribution, rms velocities, and enthalpy flux in the convective envelope. This change in behaviour is characterised by the formation of a hot layer on the surface of the accreting object, which tends to suppress convection in the envelope. We analyse the long-term effect of such a hot buffer zone on the structure and evolution of the accreting object with 1D stellar evolution calculations. We study the relevance of the assumption of redistribution of accreted energy into the stellar interior used in the literature. One conclusion is that, for a given amount of accreted energy transferred to the accreting object, a treatment assuming accretion energy redistribution throughout the stellar interior could significantly overestimate the effects on the stellar structure, in particular, on the resulting expansion.

This work is a continuation of our efforts to develop an efficient implicit solver for multidimensional hydrodynamics for the purpose of studying important physical processes in stellar interiors, such as turbulent convection and overshooting. We present an implicit solver that results from the combination of a Jacobian-Free Newton-Krylov method and a preconditioning technique tailored to the inviscid, compressible equations of stellar hydrodynamics. We assess the accuracy and performance of the solver for both 2D and 3D problems for Mach numbers down to $10^{-6}$. Although our applications concern flows in stellar interiors, the method can be applied to general advection and/or diffusion-dominated flows. The method presented in this paper opens up new avenues in 3D modeling of realistic stellar interiors allowing the study of important problems in stellar structure and evolution.

Giant pulses from Crab pulsar were observed together with Radioastron space radiotelescope, Global EVN radio telescopes and Kvazar-KVO VLBI stations to study the scattering effects in the ISM. Five observing sessions were conducted at 18 cm (EVN codes: EG060A, EG060B, EG067B, EG075) and one at 92 cm (EVN code: GS033A). We have estimated distance to the scattering screen, angular size of scattering disk, scattering time. All estimations were done in an assumption of single thin scattering screen. Also, a significant change in the shape of cross-correlation functions for space-ground baselines was found (starting from 4 up to 12 Earth diameters).

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 have studied the statistical and polarization properties of giant pulses (GPs) emitted by the millisecond pulsar B1937+21, with high sensitivity and time resolution. The observations were made in June 2005 with the 100-m Robert C. Byrd Green Bank Telescope at S-band (2052-2116 MHz) using the Mk5A VLBI recording system, with formal time resolution of 16 ns. The total observing time was about 4.5 hours; the rate of detection of GPs was about 130 per hour at the average longitudes of the main pulse (MPGPs) and 60 per hour at the interpulse (IPGPs). While the average profile shows well-defined polarization behavior, with regular evolution of the linear polarization position angle (PA), GPs exhibit random properties, occasionally having high linear or circular polarization. Neither MPGPs nor IPGPs show a preferred PA. The cumulative probability distribution (CPD) of GP pulse energy was constructed down to the level where GPs merge with regular pulses and noise. For both MPGPs and IPGPs, the CPD follows a power law with a break, the power index changing from -2.4 at high energy to -1.6 for low energy. Pulse smearing due to scattering masks the intrinsic shape and duration of the detected GPs. The smearing time varied during the observing session within a range of a few hundred nanoseconds. The measured polarization and statistical properties of GPs impose strong constraints on physical models of GPs. Some of these properties support a model in which GPs are generated by the electric discharge caused by magnetic reconnection of field lines connecting the opposite magnetic poles of a neutron star.

We present a study of shape, spectra and polarization properties of giant pulses (GPs) from the Crab pulsar at the very high frequencies of 8.5 and 15.1 GHz. Studies at 15.1 GHz were performed for the first time. Observations were conducted with the 100-m radio telescope in Effelsberg in Oct-Nov 2007 at the frequencies of 8.5 and 15.1 GHz as part of an extensive campaign of multi-station multi-frequency observations of the Crab pulsar. A selection of the strongest pulses was recorded with a new data acquisition system, based on a fast digital oscilloscope, providing nanosecond time resolution in two polarizations in a bandwidth of about 500 MHz. We analyzed the pulse shapes, polarisation and dynamic spectra of GPs as well as the cross-correlations between their LHC and RHC signals. No events were detected outside main pulse and interpulse windows. GP properties were found to be very different for GPs emitted at longitudes of the main pulse and the interpulse. Cross-correlations of the LHC and RHC signals show regular patterns in the frequency domain for the main pulse, but these are missing for the interpulse GPs. We consider consequences of application of the rotating vector model to explain the apparent smooth variation in the position angle of linear polarization for main pulse GPs. We also introduce a new scenario of GP generation as a direct consequence of the polar cap discharge. We find further evidence for strong nano-shot discharges in the magnetosphere of the Crab pulsar. The repetitive frequency spectrum seen in GPs at the main pulse phase is interpreted as a diffraction pattern of regular structures in the emission region. The interpulse GPs however have a spectrum that resembles that of amplitude modulated noise. Propagation effects may be the cause of the differences.

Simultaneous multifrequency observations of the Crab pulsar giant pulses (GPs) were performed with the 64-m Kalyazin radio telescope at four frequencies 0.6, 1.4, 2.2 and 8.3 GHz using the K5 VLBI recording terminal. The K5 terminal provided continuous recording in 16 4-MHz wide frequency channels distributed over 4 frequency bands. Several thousands of GPs were detected during about 6 hours of observations in two successive days in July 2005. Radio spectra of single GPs were analysed at separate frequencies and over whole frequency range. These spectra manifest notable modulation over frequency ranges, $\Delta\nu$, both on large ($\Delta\nu/\nu\approx 0.5$) and small ($\Delta\nu/\nu\approx 0.01$) frequency scales. Cross-correlation analysis of GPs at 2.2 GHz showed that their pulse shapes can be interpreted as an ensemble of unresolved bursts grouped together at time scales of $\approx 1$ mcs being well-correlated over a 60-MHz band. The corresponding GP cross-correlation functions do not obey the predictions of the amplitude-modulated noise model of Rickett (1975), thus indicating that unresolved components represent a small number of elementary emitters.

Stable, accurate, divergence-free simulation of magnetized supersonic turbulence is a severe test of numerical MHD schemes and has been surprisingly difficult to achieve due to the range of flow conditions present. Here we present a new, higher order-accurate, low dissipation numerical method which requires no additional dissipation or local "fixes" for stable execution. We describe PPML, a local stencil variant of the popular PPM algorithm for solving the equations of compressible ideal magnetohydrodynamics. The principal difference between PPML and PPM is that cell interface states are evolved rather that reconstructed at every timestep, resulting in a compact stencil. Interface states are evolved using Riemann invariants containing all transverse derivative information. The conservation laws are updated in an unsplit fashion, making the scheme fully multidimensional. Divergence-free evolution of the magnetic field is maintained using the higher order-accurate constrained transport technique of Gardiner and Stone. The accuracy and stability of the scheme is documented against a bank of standard test problems drawn from the literature. The method is applied to numerical simulation of supersonic MHD turbulence, which is important for many problems in astrophysics, including star formation in dark molecular clouds. PPML accurately reproduces in three-dimensions a transition to turbulence in highly compressible isothermal gas in a molecular cloud model. The low dissipation and wide spectral bandwidth of this method make it an ideal candidate for direct turbulence simulations.

We present a review of observed parameters of giant radio pulses, based on the observations conducted by our group during recent years. The observations cover a broad frequency range of about 3 octaves, concentrating between 600 and 4850 MHz. Giant pulses of both the Crab pulsar and the millisecond pulsar B1937+21 were studied with the 70-m Tidbinbilla, the 100-m GBT, 64-m Kalyazin and Westerbork radio telescopes. We discuss pulse energy distribution, dependence of peak flux density from the pulse width, peculiarities of radio spectra, and polarization properties of giant radio pulses.

We have studied the statistics of giant pulses from the Crab pulsar for the first time with particular reference to their widths. We have analyzed data collected during 3.5 hours of observations conducted with the Westerbork Synthesis Radio Telescope operated in a tied-array mode at a frequency of 1200 MHz. The PuMa pulsar backend provided voltage recording of X and Y linear polarization states in two conjugate 10 MHz bands. We restricted the time resolution to 4 microseconds to match the scattering on the interstellar inhomogeneities. In total about 18000 giant pulses (GP) were detected in full intensity with a threshold level of 6 sigma. Cumulative probability distributions (CPD) of giant pulse energies were analyzed for groups of GPs with different effective widths in the range 4 to 65 microseconds. The CPDs were found to manifest notable differences for the different GP width groups. The slope of a power-law fit to the high-energy portion of the CPDs evolves from -1.7 to -3.2 when going from the shortest to the longest GPs. There are breaks in the CPD power-law fits indicating flattening at low energies with indices varying from -1.0 to -1.9 for the short and long GPs respectively. The GPs with a stronger peak flux density were found to be of shorter duration. We compare our results with previously published data and discuss the importance of these peculiarities in the statistical properties of GPs for the heoretical understanding of the emission mechanism responsible for GP generation.

VLBI observations of the Crab pulsar with the 64-m radio telescope at Kalyazin (Russia) and the 46-m radio telescope of the Algonquin Radio Observatory (Canada) at 2.2 GHz and single-dish observations of the millisecond pulsar B1937+21 with the GBT (USA) at 2.1 GHz were conducted to probe the interstellar medium and study the properties of giant pulses. The VLBI data were processed with a dedicated software correlator, which allowed us to obtain the visibility of single giant pulses. Two frequency scales of 50 and 450 kHz were found in the diffraction spectra of giant pulses from the Crab pulsar. The location of the scattering region was estimated to be close to the outer edge of the nebula. No correlation was found between the power spectra of giant pulses at left- and right-hand circular polarization. We explain this lack of correlation through the influence of the strong magnetic field on circularly polarized emission in the region close to the Crab pulsar. Combining the measurement of the decorrelation bandwidth with that of the scattering time of giant pulses for B1937+21, we found three frequency scales of 1.7, 3.8, and 16.5 MHz. The scattering time of giant pulses of B1937+21 at 2.1 GHz was found to be 40+-4 ns. We obtained an upper limit of the intrinsic width of giant pulses from B1937+21 of less than 8 ns. The frequency dependences of the scattering times for the Crab pulsar and PSR B1937+21 were found to be different. They are characterized by exponents of -3.5 and -4.2, respectively. We attribute the difference to the large influence of scattering in the Crab nebula.

The second fastest millisecond pulsar, B1937+21, is one of several pulsars known to emit giant pulses (GPs). GPs are characterized by their huge energy, power-law cumulative energy distribution, and particular longitudes of occurrence. All these characteristics are different from those of regular pulses. Here, we present a study of GPs from our observations of the pulsar B1937+21 with the GBT at 2.1 GHz in both left and right circular polarization with a time resolution of 8 ns. The Mark5 data acquisition system was used for the first time in single-dish observations with the GBT. This allowed us to obtain continuous and uniform recording for 7.5 hours with a data rate of 512 Mbps. As a result, more than 6000 GPs were found above a detection threshold of 200 Jy. We report on instantaneous spectra of GPs, as well as on a comparison with scintillation spectra of regular emission, on the distribution of GP energies, and on polarization properties of GPs.

The results of simultaneous multifrequency observations of giant radio pulses from the Crab pulsar, PSR B0531+21, at 23, 111, and 600 MHz are presented and analyzed. Giant pulses were detected at a frequency as low as 23 MHz for the first time. Of the 45 giant pulses detected at 23 MHz, 12 were identified with counterparts observed simultaneously at 600 MHz. Of the 128 giant pulses detected at 111 MHz, 21 were identified with counterparts observed simultaneously at 600 MHz. The spectral indices for the power-law frequency dependence of the giant-pulse energies are from -3.1 to -1.6. The mean spectral index is -2.7 +/- 0.1 and is the same for both frequency combinations (600-111 MHz and 600-23 MHz). The large scatter in the spectral indices of the individual pulses and the large number of unidentified giant pulses suggest that the spectra of the individual giant pulses do not actually follow a simple power law. The observed shapes of the giant pulses at all three frequencies are determined by scattering on interstellar plasma irregularities. The scatter broadening of the pulses and its frequency dependence were determined as tau_sc=20*(f/100)^(-3.5 +/- 0.1) ms, where the frequency f is in MHz.

The paper presents an analysis of dual-polarization observations of the Crab pulsar obtained on the 64-m Kalyazin radio telescope at 600 MHz with a time resolution of 250 ns. A lower limit for the intensities of giant pulses is estimated by assuming that the pulsar radio emission in the main pulse and interpulse consists entirely of giant radio pulses; this yields estimates of 100 Jy and 35 Jy for the peak flux densities of giant pulses arising in the main pulse and interpulse, respectively. This assumes that the normal radio emission of the pulse occurs in the precursor pulse. In this case, the longitudes of the giant radio pulses relative to the profile of the normal radio emission turn out to be the same for the Crab pulsar and the millisecond pulsar B1937+21, namely, the giant pulses arise at the trailing edge of the profile of the normal radio emission. Analysis of the distribution of the degree of circular polarization for the giant pulses suggests that they can consist of a random mixture of nanopulses with 100% circular polarization of either sign, with, on average, hundreds of such nanopulses within a single giant pulse.

A flare of maser radio emission in the OH-line 1665 MHz has been discovered in the star forming region W75 N in 2003, with the flux density of about 1000 Jy. At the time it was the strongest OH maser detected during the whole history of observations since the discovery of cosmic masers in 1965. The flare emission is linearly polarized with a degree of polarization near 100%. A weaker flare with a flux of 145 Jy was observed in this source in 2000 - 2001, which was probably a precursor of the powerful flare. Intensity of two other spectral features has decreased after beginning of the flare. Such variation of the intensity of maser condensation emission (increasing of one and decreasing of the other) can be explained by passing of the magneto hydrodynamic shock across regions of enhanced gas concentration.

Giant radio pulses of the millisecond pulsar B1937+21 were recorded with the S2 VLBI system at 1.65 GHz with NASA/JPL's 70-m radio telescope at Tidbinbilla, Australia. These pulses have been observed as strong as 65000 Jy with widths <= 15 ns, corresponding to a brightness temperature T_b >= 5 x 10^39 K, the highest observed in the universe. The vast majority of these pulses occur in a 5.8 mcs and 8.2 mcs window at the very trailing edges of the regular main pulse and interpulse profiles, respectively. Giant pulses occur in general with a single spike. Only in one case out of 309 was the structure clearly more complex. The cumulative distribution is fit by a power law with index -1.40 +/- 0.01 with a low-energy but no high-energy cutoff. We estimate that giant pulses occur frequently but are only rarely detected. When corrected for the directivity factor, 25 giant pulses are estimated to be generated in one neutron star revolution alone. The intensities of the giant pulses of the main pulses and interpulses are not correlated with each other nor with the intensities or energies of the main pulses and interpulses themselves. Their radiation energy density can exceed 300 times the plasma energy density at the surface of the neutron star and can even exceed the magnetic field energy density at that surface. We therefore do not think that the generation of giant pulses is linked to the plasma mechanisms in the magnetosphere. Instead we suggest that it is directly related to discharges in the polar cap region of the pulsar.

A numerical method presented by Imshennik et al. (2002) is used to solve the two dimensional axisymmetric hydrodynamic problem on the formation of a toroidal atmosphere during the collapse of an iron stellar core and outer stellar layers. An evolutionary model from Boyes et al. (1999) with a total mass of $25M_{\odot}$ is used as the initial data for the distribution of thermodynamic quantities in the outer shells of a high-mass star. We analyze in detail the results of three calculations in which the difference mesh and the location of the inner boundary of the computational region are varied. In the initial data, we roughly specify an angular velocity distribution that is actually justified by the final result - the formation of a hydrostatic equilibrium toroidal atmosphere with reasonable total mass, $M^{tot} = (0.117 \div 0.122)M_{\odot}$, and total angular momentum, $J^{tot} = (0.445 \div 0.472) x 10^{50} erg \cdot s$, for the two main calculations. We compare the numerical solution with our previous analytical solution in the form of toroidal atmospheres (Imshennik and Manukovskii 2000). This comparison indicates that they are identical if we take into account the more general and complex equation of state with a nonzero temperature and self-gravitation effects in the atmosphere. Our numerical calculations, first, prove the stability of toroidal atmospheres on characteristic hydrodynamic time scales and, second, show the possibility of sporadic fragmentation of these atmospheres even after a hydrodynamic equilibrium is established. The calculations were carried out under the assumption of equatorial symmetry of the problem and up to relatively long time scales $(\approx 10s)$.

We analyze the behavior of the outer envelope in a massive star during and after the collapse of its iron core into a protoneutron star (PNS) in terms of the equations of one-dimensional spherically symmetric ideal hydrodynamics. The profiles obtained in the studies of the evolution of massive stars up to the final stages of their existence, immediately before a supernova explosion (Boyes et al. 1999), are used as the initial data for the distribution of thermodynamic quantities in the envelope.We use a complex equation of state for matter with allowances made for arbitrary electron degeneracy and relativity, the appearance of electron-positron pairs, the presence of radiation, and the possibility of iron nuclei dissociating into free nucleons and helium nuclei. We performed calculations with the help of a numerical scheme based on Godunov's method. These calculations allowed us to ascertain whether the emersion of the outer envelope in a massive star is possible through the following two mechanisms: first, the decrease in the gravitational mass of the central PNS through neutrino-signal emission and, second, the effect of hot nucleon bubbles, which are most likely formed in the PNS corona, on the envelope emersion. We show that the second mechanism is highly efficient in the range of acceptable masses of the nucleon bubbles ($\leq 0.01M_\odot$) simulated in our hydrodynamic calculations in a rough, spherically symmetric approximation.

We report on a study of microstructure and its quasi-periodicities of three pulsars at 1.65 GHz with the S2 VLBI system at a resolution of 62.5 ns, by far the highest for any such statistical study yet. For PSR B1929+10 we found in the average cross-correlation function (CCF) broad microstructure with a characteristic timescale of 95+-10 mcs and confirmed microstructure with characteristic timescales between 100 and 450 mcs for PSRs B0950+08 and B1133+16. On a finer scale PSRs B0950+08, B1133+16 (component II) and B1929+10 show narrow microstructure with a characteristic timescale in the CCFs of ~10 mcs, the shortest found in the average CCF or autocorrelation function (ACF) for any pulsar, apart perhaps for the Crab pulsar. Histograms of microstructure widths are skewed heavily toward shorter timescales but display a sharp cutoff. The shortest micropulses have widths between 2 and 7 mcs. No nanopulses or unresolved pulse spikes were detected. Cross-power spectra of single pulses show a large range of complexity with single spectral features representing classic quasi-periodicities and broad and overlapping features with essentially no periodicities at all. Significant differences were found for the two components of PSR B1133+16 in every aspect of our statistical analysis of micropulses and their quasi-periodicities. Asymmetries in the magnetosphere and the hollow cone of emission above the polar cap of the neutron star may be responsible for these differences.

We report on the first space-VLBI observations of the OH34.26+0.15 maser in two main line OH transitions at 1665 and 1667 MHz. The observations involved the space radiotelescope on board the Japanese satellite HALCA and an array of ground radio telescopes. The map of the maser region and images of individual maser spots were produced with an angular resolution of 1 milliarcsec which is several times higher than the angular resolution available on the ground. The maser spots were only partly resolved and a lower limit to the brightness temperature 6x10^{12} K was obtained. The maser seems to be located in the direction of low interstellar scattering, an order of magnitude lower than the scattering of a nearby extragalactic source and pulsar.