Magnetic fields of strongly magnetized stars can trap conducting matter due to frozen-in condition. In the force-free regime, the motion of the matter along the field lines may be considered in the "bead on a wire" approximation. Such a motion, if gravity and centrifugal forces are taken into account, has equilibrium points, some of which are stable. In most cases, stability is possible in about several per cent of the possible locations. Corresponding oscillation frequencies span the range from zero to $\sqrt{3}$ of the spin frequency. We suggest that this variability mode may be excited in some X-ray pulsars during the outbursts and create the peaked broad-band noise component near the break frequency in the power density spectrum, as well as produce some of the quasi-periodic oscillation features in this frequency range. Existence of this variability does not require any changes in mass accretion rate and involves only a small amount of matter infiltrating from the disc and magnetic flow due to interchange instabilities.

At present, there is no direct information about the evolution of inclination angle $\chi$ between magnetic and rotational axes in radio pulsars. As to theoretical models of pulsar evolution, they predict both the alignment, i.e. evolution of inclination angle $\chi$ to $0^{\circ}$, and its counter-alignment, i.e. evolution to $90^{\circ}$. In this paper, we demonstrate that the statistics of interpulse pulsars can give us the key test to solve the alignment/counter-alignment problem as the number of orthogonal interpulse pulsars ($\chi \approx 90^{\circ}$) drastically depends on the evolutionary trajectory.

Neutron stars are extremely strong cosmic magnets which fields are expected to decay with time. Here we report on the simple test of this process. Adopting a novel approach, we have estimated surface magnetic fields $B$ for 76 radiopulsars (the most numerous subclass of the known isolated neutron stars) which ages $t$ are known independently. Focusing on the accurate evaluation of the precision of both quantities, we determined a significant power-law trend $B(t) \propto t^{-\beta}$ with index $\beta = 0.19^{+0.05}_{-0.06}$ at 95\% C.L. The effects of the observational selection turn this value into the upper limit for the intrinsic field decay rate. If so, then neutron star crusts are close to the "impurity-free crystals", which results in a relatively slow magnetic fields decay.

The ejecta composition is an open question in gamma-ray bursts (GRB) physics. Some GRBs possess a quasi-thermal spectral component in the time-resolved spectral analysis, suggesting a hot fireball origin. Others show a featureless non-thermal spectrum known as the "Band" function, consistent with a synchrotron radiation origin and suggesting that the jet is Poynting-flux-dominated at the central engine and likely in the emission region as well. There are also bursts showing a sub-dominant thermal component and a dominant synchrotron component suggesting a likely hybrid jet composition. Here we report an extraordinarily bright GRB 160625B, simultaneously observed in gamma-rays and optical wavelengths, whose prompt emission consists of three isolated episodes separated by long quiescent intervals, with the durations of each "sub-burst" being $\sim$ 0.8 s, 35 s, and 212 s, respectively. Its high brightness (with isotropic peak luminosity L$_{\rm p, iso}\sim 4\times 10^{53}$ erg/s) allows us to conduct detailed time-resolved spectral analysis in each episode, from precursor to main burst and to extended emission. The spectral properties of the first two sub-bursts are distinctly different, allowing us to observe the transition from thermal to non-thermal radiation between well-separated emission episodes within a single GRB. Such a transition is a clear indication of the change of jet composition from a fireball to a Poynting-flux-dominated jet.

Numerical simulations of realistic non-vacuum magnetospheres of isolated neutron stars have shown that pulsar spin-down luminosities depend weakly on the magnetic obliquity. This result provides the opportunity to estimate the surface magnetic field for a given radiopulsar with quite a high accuracy. In the current work, we present a refinement of the classical `magneto-dipolar' formula for pulsar magnetic fields $B_{\rm md} = (3.2\times 10^{19}\mbox{ G})\sqrt{P\dot P}$, where $P$ is the neutron star spin period. The new, robust timing-based estimator is introduced as $\log B = \log B_{\rm md} + \Delta_{\rm B}(M, \alpha)$, where the correction $\Delta_{\rm B}$ depends on the equation of state (EOS) of dense matter, the individual pulsar obliquity $\alpha$ and the mass $M$. Adopting state-of-the-art statistics for $M$ and $\alpha$ we calculate the distributions of $\Delta_{\rm B}$ for a representative subset of 22 EOSs that do not contradict observations. It has been found that $\Delta_{\rm B}$ is distributed nearly normally, with the average in the range -0.5 to -0.25 dex and standard deviation $\sigma[\Delta_{\rm B}] \approx$ 0.06 to 0.09 dex, depending on the adopted EOS. The latter quantity represents a formal uncertainty of the corrected estimation of $\log B$ because $\Delta_{\rm B}$ is weakly correlated with $\log B_{\rm md}$. At the same time, if it is assumed that every considered EOS has the same chance of occurring in nature, then another, more generalized, estimator $B^* \approx 3B_{\rm md}/7$ can be introduced providing an unbiased value of the pulsar surface magnetic field with $\sim$30 per cent uncertainty with 68 per cent confidence. Finally, we discuss the possible impact of pulsar timing irregularities on the timing-based estimation of $B$ and review the astrophysical applications of the obtained results.

Here we report on the results of meteor observations with 9-channel Mini-MegaTORTORA (MMT-9) wide-field optical monitoring system with high temporal resolution. During first 1.5 years of operation more than 90 thousands of meteors have been detected, at a rate of 300-350 per night, with durations from 0.1 to 2.5 seconds and angular velocities up to 38 degrees per second. The faintest detected meteors has the peak brightness about 10 mag, while the majority - from 4 to 8 mag. Some of the meteors have been observed in BVR filters simultaneously. Color variations along the trail for them are determined. All parameters of detected meteors are published online. The database also includes the information on 10 thousands meteors detected by our previous FAVOR camera in 2006-2009 years.

We propose an improved method for the atmospheric extinction reduction within optical photometry. Our method is based on the simultaneous multicolor observations of photometric standards. Such data are now available within the modern wide-field sky surveys and contain a large amount of information about instant atmospheric conditions. So, it became possible to estimate the extinction parameters on the basis of a quite short observational dataset and, hence, to trace the rapid stars twinkling accurately. Having been developed for a new MiniMegaTORTORA observational system, the proposed method can be adopted for a wide range of modern observational programs.

Here we describe the first light of the novel 9-channel wide-field optical monitoring system with sub-second temporal resolution, Mini-MegaTORTORA, which is being tested now at Special Astrophysical Observatory on Russian Caucasus. The system is able to observe the sky simultaneously in either wide (~900 square degrees) or narrow (~100 square degrees) fields of view, either in clear light or with any combination of color (Johnson B, V or R) polarimetric filters installed, with exposure times ranging from 100 ms to 100 s. The primary goal of the system is the detection of rapid -- with sub-second characteristic time-scales -- optical transients, but it may be also used for studying the variability of the sky objects on longer time scales.

We have developed a method for analyzing the kinematic association of isolated relativistic objects - possible remnants of disrupted close binary systems. We investigate pairs of fairly young radio pulsars with known proper motions and estimated distances (dispersion measures) that are spaced no more than 2-3 kpc apart. Using a specified radial velocity distribution for these objects, we have constructed 100-300 thousand trajectories of their possible motion in the Galactic gravitational field on a time scale of several million years. The probabilities of their close encounters at epochs consistent with the age of the younger pulsar in the pair are analyzed. When these probabilities exceed considerably their reference values obtained by assuming a purely random encounter between the pulsars under consideration, we conclude that the objects may have been gravitationally bound in the past. As a result, we have detected six pulsar pairs (J0543+2329/J0528+2200, J1453-6413/J1430-6623, J2354+6155/J2321+6024, J1915+1009/J1909+1102, J1832-0827/J1836-1008, and J1917+1353/J1926+1648) that are companions in disrupted binary systems with a high probability. Estimates of their kinematic ages and velocities at binary disruption and at the present epoch are provided.

In this article we revise the problem of anomalous values of pulsars' braking indices n_{obs} and frequency second derivatives F2 arising in observations. The intrinsic evolutionary braking is buried deep under superimposed irregular processes, that prevent direct estimations of its parameters for the majority of pulsars. We re-analyze the distribution of "ordinary" radio pulsars on a F2-F1, F2-F0, F1-F0 and n_{obs}-tau_{ch} diagrams assuming their spin-down to be the superposition of a "true" monotonous term and a symmetric oscillatory term. We demonstrate that their effects may be clearly separated using simple ad hoc arguments. Using maximum likelihood estimator we derive the parameters of both components. We find characteristic timescales of such oscillations to be of the order of 1e3-1e4 years, while its amplitudes are large enough to modulate the observed spin-down rate up to 0.5-5 times and completely dominate the second frequency derivatives. On the other hand, pulsars' secular evolution is consistent with classical magnetodipolar model with braking index n ~ 3. So, observed pulsars' characteristic ages (and similar estimators that depend on the observed F1) are also affected by long term cyclic process and differ up to 0.5-5 times from their monotonous values. This fact naturally resolves the discrepancy of characteristic and independently estimated physical ages of several objects, as well as explains very large, up to 1e8 years, characteristic ages of some pulsars. We discuss the possible physical connection of long term oscillation with a complex neutron star rotation relative its magnetic axis due to influence of the near-field part of magnetodipolar torque.

The measurements of pulsar frequency second derivatives have shown that they are 1e2...1e6 times larger than expected for standard pulsar spin-down law. Moreover, the second derivatives as well as braking indices are even negative for about half the pulsars. We explain these paradoxical results on the basis of the statistical analysis of the rotational parameters f0, f1 and f2 of the subset of 295 pulsars taken mostly from the ATNF database. We have found a strong correlation between f2 and f1 for both f2 > 0 (correlation coefficient r ~ 0.9) and f2 < 0 (r ~ 0.85), as well as between f0 and f1 (r ~ 0.6...0.7). We interpret these dependencies as evolutionary ones due to f1 being nearly proportional to the pulsars' age. The derived statistical relations as well as "anomalous" values of f2 are well described by assuming the existence of long-time variations of the spin-down rate. The pulsar frequency evolution, therefore, consists of secular change of f0_{ev}(t), f1_{ev}(t) and f2_{ev}(t) according to the power law with n ~ 5, the irregularities, observed within the timespan as timing noise, and the non-monotonous variations on the timescale of several tens of years, which is larger than that of the timespan. It is possible that the nature of long-term variations is similar to that of short-term ones. The idea of non-constant secular pulsars' braking index n is also analysed.

General properties of accretion onto isolated stellar mass black holes in the Galaxy are discussed. An analysis of plasma internal energy growth during the infall is performed. Adiabatic heating of collisionless accretion flow due to magnetic adiabatic invariant conservation is 25% more efficient than in the standard non-magnetized gas case. It is shown that magnetic field line reconnections in discrete current sheets lead to significant nonthermal electron component formation, which leads to a formation of a hard (UV, X-ray, up to gamma), highly variable spectral component in addition to the standard synchrotron optical component first derived by Shvartsman generated by thermal electrons in the magnetic field of the accretion flow. Properties of accretion flow emission variability are discussed. Observation results of two single black hole candidates - gravitational lens MACHO-1999-BLG-22 and radio-loud x-ray source with featureless optical spectrum J1942+10 - in optical band with high temporal resolution are presented and interpreted in the framework of the proposed model.

The stability of the optical pulse of the Crab pulsar is analyzed based on the 1 $\mu$s resolution observations with the Russian 6-meter and William Hershel telescopes equipped with different photon-counting detectors. The search for the variations of the pulse shape along with its arrival time stability is performed. Upper limits on the possible short time scale free precession of the pulsar are placed. The evidence of pulse time of arrival (TOA) variations on 1.5-2 hours time scale is presented, along with evidence of small light curve (shape and separation of main and secondary peaks) changes between data sets, on time scale of years. Also, the fine structure of the main pulse is studied.

We present a correlation between two intrinsic parameters of GRB optical afterglows. These are the isotropic luminosity at the maximum of the light curve (Lpeak) and the time-integrated isotropic energy (Eiso) radiated after the observed maximum. We test the correlation between the logarithms of (Eiso) and (Lpeak) and finally we value the effect of the different samples of GRBs in according with the first optical observation reduced to proper time.

The measurements of pulsar frequency second derivatives have shown that they are $10^2-10^6$ times larger than expected for standard pulsar spin-down law, and are even negative for about half of pulsars. We explain these paradoxical results on the basis of the statistical analysis of the rotational parameters $\nu$, $\dot \nu$ and $\ddot \nu$ of the subset of 295 pulsars taken mostly from the ATNF database. We have found a strong correlation between $\ddot \nu$ and $\dot \nu$ for both $\ddot\nu > 0$ and $\ddot\nu < 0$, as well as between $\nu$ and $\dot\nu$. We interpret these dependencies as evolutionary ones due to $\dot\nu$ being nearly proportional to the pulsars' age. The derived statistical relations as well as "anomalous" values of $\ddot\nu$ are well described by assuming the long-time variations of the spin-down rate. The pulsar frequency evolution, therefore, consists of secular change of $\nu_{ev}(t)$, $\dot\nu_{ev}(t)$ and $\ddot\nu_{ev}(t)$ according to the power law with $n \approx 5$, the irregularities, observed within a timespan as a timing noise, and the variations on the timescale larger than that timespan -- several tens of years.

The region of SGR 1806-20 localization was observed during its gamma-ray activity in 2001. The observations have been performed on the 6-meter telescope of the Special Astrophysical Observatory, using the Panoramic Photometer-Polarimeter (PPP). The search for variability was performed on the $10^{-6}$ - 10 s time scale, and its results were compared to the properties of corresponding x-ray flares.

The fine structure and the variations of the optical pulse shape and phase of the Crab pulsar are studied on various time scales. The observations have been carried out on 4-m William Hershel and 6-m BTA telescopes with APD photon counter, photomultiplier based 4-channel photometer and PSD based panoramic spectrophotopolarimeter with 1$\mu$s time resolution in 1994, 1999, 2003 and 2005-2006 years. The upper limit on the pulsar precession on Dec 2, 1999 is placed in the 10 s - 2 hours time range. The evidence of a varying from set to set fine structure of the main pulse is found in the 1999 and 2003 years data. No such fine structure is detected in the integral pulse shape of 1994, 1999 and 2003 years. The drastic change of the pulse shape in the 2005-2006 years set is detected along with the pulse shape variability and quasi-periodic phase shifts.

The recent massive measurements of pulsar frequency second derivatives have shown that they are 100-1000 times larger than expected for standard pulsar slowdown low. Moreover, the second derivatives as well as braking indices are even negative for about half of pulsars. We explain these paradoxical results on the basis of the statistical analysis of the rotational parameters $\nu$, $\dot \nu$ and $\ddot \nu$ of the subset of 295 pulsars taken mostly from the ATNF database. We have found strong correlation of $\ddot \nu$ and $\dot \nu$ either for $\ddot\nu > 0$ (correlation coefficient $r\approx0.9$) and $\ddot\nu < 0$ ($r\approx0.85$), and of $\nu$ and $\dot\nu$ ($r\approx0.7$). We interpret these dependencies as evolutionary ones due to $\dot\nu$ being nearly proportional to characteristic age $\tau_{ch}$. The derived statistical relations as well as "anomalous" values of $\ddot\nu$ are well explained in the framework of the simple model of cyclic evolution of the rotational frequency of the pulsars. It combines the secular change of $\nu_{tr}(t)$, $\dot\nu_{tr}(t)$ and $\ddot\nu_{tr}(t)$ according to the power law with $n\approx5$ and harmonic oscillations of 100--1000 years period with an amplitude from $10^{-3}$ Hz for young pulsars to $10^{-10}$ Hz for elder ones. The physical nature of these cyclic variations of the rotational frequency may be similar to the well-known red timing noise, however, with much larger characteristic time scale.

The wide field optical camera with high temporal resolution for the continuous monitoring of the sky in order to catch the initial stages of GRBs is described.

An astronomical complex intended to detect optical transients (OTs) in a wide field and follow them up with high time resolution investigation is described.