GRANDMA is a network of 25 telescopes of different sizes, including both photometric and spectroscopic facilities. The network aims to coordinate follow-up observations of gravitational-wave candidate alerts, especially those with large localisation uncertainties, to reduce the delay between the initial detection and the optical confirmation. In this paper, we detail GRANDMA's observational performance during Advanced LIGO/Advanced Virgo Observing Run 3 (O3), focusing on the second part of O3; this includes summary statistics pertaining to coverage and possible astrophysical origin of the candidates. To do so, we quantify our observation efficiency in terms of delay between gravitational-wave candidate trigger time, observations, and the total coverage. Using an optimised and robust coordination system, GRANDMA followed-up about 90 % of the gravitational-wave candidate alerts, i.e. 49 out of 56 candidates. This led to coverage of over 9000 deg2 during O3. The delay between the gravitational-wave candidate trigger and the first observation was below 1.5 hour for 50 % of the alerts. We did not detect any electromagnetic counterparts to the gravitational-wave candidates during O3, likely due to the very large localisation areas (on average thousands of degrees squares) and relatively large distance of the candidates (above 200 Mpc for 60 % of BNS candidates). We derive constraints on potential kilonova properties for two potential binary neutron star coalescences (GW190425 and S200213t), assuming that the events' locations were imaged.

Context: GR290 (M33 V0532=Romano's star) is a luminous M33 object undergoing photometric variability typical for luminous blue variable (LBV) stars. It lies inside Wolf-Rayet region in the Hertzsprung-Russell diagram and possesses a WN8 type spectrum at the light minima. Analysis of Gran Telescopio Canarias (GTC) spectra obtained in 2016 led to the conclusion that it is surrounded by an unresolved HII region formed mostly of ejected material from the central star, and disclosed the presence of a second, more extended asymmetrical emission region. Aims: The aim of this paper is to further explore the structure of the nearby environment of GR290. Methods: Long-slit spectra of GR290 were obtained with three slit orientations in the visual and red spectral regions. The emission-line distribution for each slit was analyzed. Results: We confirm the presence of an asymmetric HII region that extends ~50 pc to the south; ~30pc to the north and southeast; ~20 pc to the east and northwest and ~10pc to the west. We also present the first spectrum to be acquired of a star belonging to the neighboring OB88 association, J013501.87+304157.3, which we classify as a B-type supergiant with a possible binary companion.

Advanced knowledge of the detailed atmospheric properties of both the future sites of the Cherenkov Telescope Array is essential in preparation of the arrival of the first scientific data. Meteorological variables are studied using a dedicated characterization station installed at the southern site in Chile and a wealth of data from existing observatories around the northern site on the La Palma island. Campaigns using radiosondes launched on balloons are foreseen to complement these data in the near future. Cloudiness during the night has been continuously monitored at both sites for several years using All-sky Cameras which assess the presence of clouds based on detection of stars. The integrated aerosol optical depth over the southern site has been measured using a Sun/Moon Photometer since 2016 and the small robotic FRAM telescope since 2017; identical instruments have been deployed at the northern site in autumn 2018. Also in October 2018, the ARCADE Raman lidar (RL) has started to take measurements on routine basis at the northern site, providing data on the vertical profile of the aerosol optical properties (i.e., extinction and scattering) and of the water vapour mixing ratio. We present the data currently available from these instruments from both sites with emphasis on characteristics important for the (future) operation of Imaging Atmospheric Cherenkov Telescopes.

The atmospheric monitoring devices for the planned calibration system of the Cherenkov Telescope Array (CTA) are undergoing intensive development, prototyping and testing. The All-Sky Cameras, the Sun/Moon Photometers and the FRAM telescopes have been gradually deployed at the future CTA sites with the primary goal of site characterization, simultaneously allowing the assessment of their operational reliability in realistic environmental conditions. All three devices have shown the ability to work smoothly in both the extreme dryness and the large temperature variations of the southern site as well as in the occasional adverse weather during winter months at the northern site. The target availability of 95\% of time has not yet been reached mostly due to minor hardware failures that have proven difficult to fix due to the remoteness of the installation in the absence of the future CTA infrastructure. The experience gathered during the prototype operations will contribute to the improved reliability of the final instruments. The Raman LIDARs, described in separate proceedings of this conference, and the infrared Ceilometer, ready for testing in Prague, will complement the set of atmospheric calibration devices in near future. The final operational procedures for the atmospheric calibration of the CTA during its operation are being finalized foreseeing the use of the All-sky Cameras and the Ceilometer for the monitoring of clouds over the whole sky and the LIDARs and FRAMs for precision measurements of the atmospheric transmission as a function of altitude and position within the field-of-view of the CTA array.

Scientific CMOS image sensors are a modern alternative for a typical CCD detectors, as they offer both low read-out noise, large sensitive area, and high frame rates. All these makes them promising devices for a modern wide-field sky surveys. However, the peculiarities of CMOS technology have to be properly taken into account when analyzing the data. In order to characterize these, we performed an extensive laboratory testing of Andor Marana sCMOS camera. Here we report its results, especially on the temporal stability and linearity, and compare it to the previous versions of Andor sCMOS cameras. We also present the results of an on-sky testing of this sensor connected to a wide-field lens, and discuss its applications for an astronomical sky surveys.

To assess the potential of a small telescopes for a high temporal resolution astrophysics, we observed the field of FRB~121102 repeating source of fast radio bursts on a 50-cm D50 telescope of Ond\v{r}ejov observatory, equipped with a fast frame rate EMCCD detector. In a three nights of observations we did not detect any optical flares from the source, which allows us to place an upper limit of 10 mJy for a brightness of possible fast optical events on a time scale of 10 ms. We also characterize the apparent brightness stability of a field stars on the same time scale in order to investigate the potential of such telescopes for detecting faint optical variability on a sub-second time scales.

AG Draconis is a strongly interacting binary system which manifests characteristic symbiotic activity of alternating quiescent and active stages. The latter ones consist of the series of individual outbursts repeating at about a one-year interval. After seven years of flat quiescence following the 2006-2008 major outbursts, in the late spring of 2015, the symbiotic system AG Dra started to become brighter again toward what appeared to be a new minor outburst. The current outburst activity of AG Dra was confirmed by the following three outbursts in April 2016, May 2017 and April 2018. The photometric and spectroscopic observations suggest that all these outbursts are of the hot type. Such behaviour is considerably peculiar in almost 130-year history of observing of this object, because the major outbursts at the beginning of active stages are typically cool ones. In the present work, the current peculiar activity of the symbiotic binary AG Dra is described in detail.

Results of comparison of Gaia DR2 parallaxes with data derived from a combined analysis of 2MASS (Two Micron All-Sky Survey), SDSS (Sloan Digital Sky Survey), GALEX (Galaxy Evolution Explorer), and UKIDSS (UKIRT Infrared Deep Sky Survey) surveys in four selected high-latitude $|b|>48^{\circ}$ sky areas are presented. It is shown that multicolor photometric data from large modern surveys can be used for parameterization of stars closer than 4400 pc and brighter than $g_{SDSS} = 19.^m6$, including estimation of parallax and interstellar extinction value. However, the stellar luminosity class should be properly determined.

A FRAM (F/(Ph)otometric Robotic Atmospheric Monitor) telescope is a system of a robotic mount, a large-format CCD camera and a fast telephoto lens that can be used for atmospheric monitoring at any site when information about the atmospheric transparency is required with high spatial or temporal resolution and where continuous use of laser-based methods for this purpose would interfere with other observations. The original FRAM has been operated at the Pierre Auger Observatory in Argentina for more than a decade, while three more FRAMs are foreseen to be used by the Cherenkov Telescope Array (CTA). The CTA FRAMs are being deployed ahead of time to characterize the properties of the sites prior to the operation of the CTA telescopes; one FRAM has been running on the planned future CTA site in Chile for a year while two others are expected to become operational before the end of 2018. We report on the hardware and current status of operation and/or deployment of all the FRAM instruments in question as well as on some of the preliminary results of integral aerosol measurements by the FRAMs in Argentina and Chile

The idea of using stellar photometry for atmospheric monitoring for optical experiments in high-energy astrophysics is seemingly straightforward, but reaching high precision of the order of 0.01 in the determination of the vertical aerosol optical depth (VAOD) has proven difficult. Wide-field photometry over a large span of altitudes allows a fast determination of VAOD independently of the absolute calibration of the system, while providing this calibration as a useful by-product. Using several years of data taken by the FRAM (F/(Ph)otometric Robotic Atmospheric Monitor) telescope at the Pierre Auger Observatory in Argentina and about a year of data taken by a similar instrument deployed at the planned future Southern site of the Cherenkov Telescope Array in Chile, we have developed methods to improve the precision of this measurement technique towards and possibly beyond the 0.01 mark. Detailed laboratory measurements of the response of the whole system to both the spectrum and intensity of incoming light have proven indispensable in this analysis as the usual assumption of linearity of the CCD detectors is not valid anymore for the conditions of the observations.

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.

During our optical monitoring of UV Ceti, iconic late-type flaring star, with high temporal resolution using the Russian 6-m telescope in 2008 we detected a giant flare with the amplitude of about 3 magnitudes in $U$-band. Near flare maximum more than a dozen of spike bursts have been discovered with triangular shapes and durations from 0.6 to 1.2 s and maximal luminosities in the range $(1.5-8)\times10^{27}$ erg~s$^{-1}$. For the half of these events the linear polarization exceeds 35% with significance better than $5\sigma$. We argue that these events are synchrotron emission of electron streams with the energies of several hundred MeV moving in the magnetic field of about 1.4 kG. Emission from such ultrarelativistic (with energies far exceeding 10 MeV) particles is being routinely observed in solar flares, but has never been detected from UV Ceti type stars. This is the first ever detection of linearly polarized optical light from the UV Ceti-type stars which indicates that at least some fraction of the flaring events on these stars are powered by a non-thermal synchrotron emission mechanism.

To explain the nature of the anomalously high reddening (Av~10 mag) towards one of the most luminous stars in the Galaxy - Cyg OB2 #12 (B5 Ia-0), also known as MT304, we carried out spectro-photometric observations of 24 stars located in its vicinity. We included five new B-stars to the members of Cygnus OB2, and for five more photometrically selected stars we spectroscopically confirmed their membership. We constructed the map of interstellar extinction within 2.5 arcmin radius and found that interstellar extinction increases towards MT304. The increase of reddening suggests that the reddening excess may be caused by the circumstellar shell ejected by the star during its evolution. We also report the detection of a second component of MT304, and discovery of an even fainter third component, based on data of speckle interferometric observations taken with the Russian 6-m telescope.

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.

To explain the nature of the high reddening (Av~10 mag) towards one of the most luminous stars in the Galaxy - Cyg OB2 #12 (B5 Ia-0), also known as MT304, we carried out spectro-photometric observations of 24 stars located in its vicinity. We included five new B-stars among the members of Cygnus OB2, and for five more photometrically selected stars we spectroscopically confirmed their membership. We constructed the map of interstellar extinction within 2.5 arcmin radius and found that interstellar extinction increases towards MT304. According to our results the most reddened OB-stars in the association after MT304 are J203240.35+411420.1 and J203239.90+411436.2, located about 15 arcsec away from it. Interstellar extinction towards these stars is about 9 mag. The increase of reddening towards MT304 suggests that the reddening excess may be caused by the circumstellar shell ejected by the star during its evolution. This shell absorbs 1 mag, but its chemical composition and temperature are unclear. We also report the detection of a second component of MT304, and discovery of an even fainter third component, based on data of speckle interferometric observations taken with the Russian 6-m telescope.

We present the statistical analysis of the properties of gamma-ray bursts with measured host galaxy redshifts and peaked optical light curves in proper frames of reference. The optical transients are classified by comparing the time lag of the optical peak relative to the GRB trigger with the duration of the gamma-ray emission itself. The results of the correlation analysis of all possible pairs of energy, spectral, and temporal characteristics of both gamma-ray and optical emissions are given. We specify the pairs of the parameters with correlation coefficients greater than 50 % at significance levels better than 1 %. The following empirical relations, obtained for the first time, are specifically discussed: a correlation between the peak optical afterglow $R$ band luminosity and redshift $L_{R} \propto (z+1)^{5.39 \pm 0.74}$ and a correlation between the peak luminosity of the prompt optical emissions and the time of the peak $L_{R} \propto T_{\rm peak}^{-3.85 \pm 1.22}$. We also analyze the similarity of the relationships between the peak optical luminosity and the isotropic equivalent of the total energy of gamma-ray bursts for afterglows ($L_{R} \propto E_{\rm iso}^{1.06 \pm 0.22}$) and for prompt optical emissions ($L_{R} \propto E_{\rm iso}^{1.59 \pm 0.21}$).

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.

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.

Prompt {\gamma}-ray emissions from gamma-ray bursts (GRBs) exhibit a vast range of extremely complex temporal structures with a typical variability time-scale significantly short - as fast as milliseconds. This work aims to investigate the apparent randomness of the GRB time profiles making extensive use of nonlinear techniques combining the advanced spectral method of the Singular Spectrum Analysis (SSA) with the classical tools provided by the Chaos Theory. Despite their morphological complexity, we detect evidence of a non stochastic short-term variability during the overall burst duration - seemingly consistent with a chaotic behavior. The phase space portrait of such variability shows the existence of a well-defined strange attractor underlying the erratic prompt emission structures. This scenario can shed new light on the ultra-relativistic processes believed to take place in GRB explosions and usually associated with the birth of a fast-spinning magnetar or accretion of matter onto a newly formed black hole.

We analyze the randomness of the sky distribution of cosmic gamma-ray bursts. These events are associated with massive galaxies, spiral or elliptical, and therefore their positions should trace the large-scale structure, which, in turn, could show up in the sky distribution of fluctuations of the cosmic microwave background (CMB). We test this hypothesis by mosaic correlation mapping of the distributions of CMB peaks and burst positions, find the distribution of these two signals to be correlated, and interpret this correlation as a possible systematic effect.

We present the observed/intrinsic optical parameters and the variability analysis of the Naked-Eye Burst, GRB 080319B, observed by the TORTORA wide-field optical monitoring system. The event is extreme not only in observed properties but also intrinsically: it is the most luminous event ever recorded at optical wavelengths. The temporal properties suggest short-lived periodic activities of the internal engine. This is the fastest optically variable source detected at cosmological distances.

We imaged the position of the Naked-Eye Burst, GRB080319B, before, during and after its gamma-ray activity with sub-second temporal resolution and discovered the fast variability of its prompt optical emission. Its characteristics and similarity with properties of gamma emission temporal structure suggest that it reflects the behaviour of internal engine -- supposedly, a hyperaccreting solar-mass black hole formed in the collapse of a massive stellar core.

Highly supercritical accretion discs are probable sources of dense optically thick axisymmetric winds. We introduce a new approach based on diffusion approximation radiative transfer in a funnel geometry and obtain an analytical solution for the energy density distribution inside the wind assuming that all the mass, momentum and energy are injected well inside the spherization radius. This allows to derive the spectrum of emergent emission for various inclination angles. We show that self-irradiation effects play an important role altering the temperature of the outcoming radiation by about 20% and the apparent X-ray luminosity by a factor of 2-3. The model has been successfully applied to two ULXs. The basic properties of the high ionization HII-regions found around some ULXs are also easily reproduced in our assumptions.

Long duration gamma-ray bursts (GRBs) release copious amounts of energy across the entire electromagnetic spectrum, and so provide a window into the process of black hole formation from the collapse of a massive star. Over the last forty years, our understanding of the GRB phenomenon has progressed dramatically; nevertheless, fortuitous circumstances occasionally arise that provide access to a regime not yet probed. GRB 080319B presented such an opportunity, with extraordinarily bright prompt optical emission that peaked at a visual magnitude of 5.3, making it briefly visible with the naked eye. It was captured in exquisite detail by wide-field telescopes, imaging the burst location from before the time of the explosion. The combination of these unique optical data with simultaneous gamma-ray observations provides powerful diagnostics of the detailed physics of this explosion within seconds of its formation. Here we show that the prompt optical and gamma-ray emissions from this event likely arise from different spectral components within the same physical region located at a large distance from the source, implying an extremely relativistic outflow. The chromatic behaviour of the broadband afterglow is consistent with viewing the GRB down the very narrow inner core of a two-component jet that is expanding into a wind-like environment consistent with the massive star origin of long GRBs. These circumstances can explain the extreme properties of this GRB.

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 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.

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.

SS433 is the only known persistent supercritical accretor, it may be very important for understanding ultraluminous X-ray sources (ULXs) located in external galaxies. We describe main properties of the SS433 supercritical accretion disk and jets. Basing on observational data of SS433 and published 2D simulations of supercritical accretion disks we estimate parameters of the funnel in the disk/wind of SS 433. We argue that the UV radiation of the SS433 disk (~50000 K, ~10^{40}erg/s) is roughly isotropic, but X-ray radiation (~10^7 K, ~10^{40}erg/s) of the funnel is midly anisotropic. A face-on SS433 object has to be ultraluminous in X-rays (10^{40-41}erg/s). Typical time-scales of the funnel flux variability are estimated. Shallow and very broad (0.1-0.3c) and blue-shifted absorption lines are expected in the funnel X-ray spectrum.

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.

A new type of cosmic-ray experiment is under construction in the Pyh\"asalmi mine in the underground laboratory of the University of Oulu, Finland. It aims to study the composition of cosmic rays at and above the knee region. The experiment, called EMMA, will cover approximately 150 square-metres of detector area. The array is capable of measuring the multiplicity and the lateral distribution of underground muons, and the arrival direction of the air shower. The full-size detector is expected to run by the end of 2007.

We suggest that the ultraluminous X-ray sources located in external galaxies (ULXs) are supercritical accretion disks like that in SS433, observed close to the disk axis. We estimate parameters of the SS433 funnel, where the relativistic jets are formed. Emergent X-ray spectrum in the proposed model of the multicolor funnel (MCF) is calculated. The spectrum can be compared with those of ULXs. We predict a complex absorption-line spectrum with broad and shallow K$\alpha$/Kc blends of the most abundant heavy elements and particular temporal variability. Another critical idea comes from observations of nebulae around the ULXs. We present results of 3D-spectroscopy of nebulae of two ULXs located in Holmberg II and NGC6946. In both cases the nebula is found to be powered by the central black hole. The nebulae are compared with SS433 nebula (W50).

Magnetic field behaviour in a spherically-symmetric accretion flow for parameters typical of single black holes in the Galaxy is discussed. It is shown that in the majority of Galaxy volume, accretion onto single stellar-mass black holes will be spherical and have a low accretion rate ($10^{-6} - 10^{-9}$ of the Eddington rate). 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. In a framework of quasi-diffusion acceleration, the "energy-radius" electron distribution is computed and the function describing the shape of synchrotron radiation spectrum is constructed. It is shown that nonthermal electron emission leads to 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 accretion flow. For typical interstellar medium parameters, a black hole at 100 pc distance will be a 16-25$^{\rm m}$ optical source coinciding with the highly variable bright X-ray counterpart, while the variable component of optical emission will be about 18-27$^{\rm m}$. The typical time scale of the variability is $10^{-4}$ sec, with relative flare amplitudes of 0.2-6% in various spectral bands. Possible applications of these results to the problem of search for single black holes are discussed.

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

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.