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.

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

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.

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

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

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.

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.

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.

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.

The Position-Sensitive Detector (PSD) for photometrical and spectral observation on the 6-meter optical telescope of the Special Astrophysical Observatory (Russia) is described. The PSD consists of a position-sensitive tube, amplifiers of output signals, analog-to-digital converters (ADC) and a digital logic plate, which produces a signal for ADC start and an external strob pulse for reading information by registration system. If necessary, the thermoelectric cooler can be used. The position-sensitive tube has the following main elements: a photocathode, electrodes of inverting optics, a block of microchannel plates (MCP) and a position-sensitive collector of quadrant type. The main parameters of the PSD are the diameter of the sensitive surface is 25 mm, the spatial resolution is better than 100 (\mu)m in the centre and a little worse on the periphery; the dead time is near 0.5 (\mu)s; the detection quantum efficiency is defined by the photocathode and it is not less than 0.1, as a rule; dark current is about hundreds of cps, or less, when cooling. PSD spectral sensitivity depends on the type of photocathode and input window material. We use a multialkali photocathode and a fiber or UV-glass, which gives the short- wave cut of 360 nm or 250 nm, respectively.

Multicolor Panoramic Photometer-Polarimeter (MPPP) with a time resolution of 1 microsecond has been built based on a PSD and used at the 6-meter telescope in SAO (Russia). The device allows registration of the photon fluxes in four photometric bands simultaneously and finding values of 3 Stokes parameters. MPPP consists of Position-Sensitive Detector (PSD), acquisition MANIA-system, polarization unit and a set of dichroic filters. MPPP gives a possibility of detecting photons in 2 pupils with a size of 10 - 15 arc sec centered on the object and comparison star positions simultaneously. The first half of the object photon flux passes through the phase rotating plate and polarizer, and the second one through the polarizer alone. MPPP registers in each of the 4 filters four images of the object with different orientations of polarization plane and one image of a comparison star. It allows measuring instantaneous Stokes parameters. The main astrophysical problems to be solved with MPPP are as follows: investigation of optical pulsars; study of GRB phenomenon in the optical range; searching for single black holes; study of fast variability of X-ray binaries. As an illustration of MPPP use, the results of observations at the 6-meter telescope of Crab pulsar and soft gamma repeater are presented.

We consider the gravitational magnification of light for binary systems containing two compact objects: white dwarfs, a white dwarf and a neutron star or a white dwarf and a black hole. Light curves of the flares of the white dwarf caused by this effect were built in analytical approximations and by means of numerical calculations. We estimate the probability of the detection of these events in our Galaxy for different types of binaries and show that gravitational lensing provides a tool for detecting such systems. We propose to use the facilities of the Sloan Digital Sky Survey (SDSS) to search for these flares. It is possible to detect several dozens compact object pairs in such a programme over 5 years. This programme is apparently the best way to detect stellar mass black holes with open event horizons.

We propose using large Air Cerenkov Telescopes (ACT's) to search for optical, pulsed signals from extra-terrestrial intelligence. Such dishes collect tens of photons from a nanosecond-scale pulse of isotropic equivalent power of tens of solar luminosities at a distance of 100 pc. The field of view for giant ACT's can be on the order of ten square degrees, and they will be able to monitor 10 to 10$^2$ stars simultaneously for nanosecond pulses of about 6th mag or brighter. Using the Earth's diameter as a baseline, orbital motion of the planet could be detected by timing the pulse arrival times.

We report a phase-resolved photometric and morphological analysis of UBV data of the Crab pulsar obtained with the 2-d TRIFFID high speed optical photometer mounted on the Russian 6m telescope. By being able to accurately isolate the pulsar from the nebular background at an unprecedented temporal resolution (1 \mu s), the various light curve components were accurately fluxed via phase-resolved photometry. Within the $UBV$ range, our datasets are consistent with the existing trends reported elsewhere in the literature. In terms of flux and phase duration, both the peak Full Width Half Maxima and Half Width Half Maxima decrease as a function of photon energy. This is similarly the case for the flux associated with the bridge of emission. Power-law fits to the various light curve components are as follows; \alpha = 0.07 \pm 0.19 (peak 1), \alpha = -0.06 \pm 0.19 (peak 2) and \alpha = -0.44 \pm 0.19 (bridge) - the uncertainty here being dominated by the integrated CCD photometry used to independently reference the TRIFFID data. Temporally, the main peaks are coincident to \le 10 \mu s although an accurate phase lag with respect to the radio main peak is compromised by radio timing uncertainties. The plateau on the Crab's main peak was definitively determined to be \leq 55 \mu s in extent and may decrease as a function of photon energy. There is no evidence for non-stochastic activity over the light curves or within various phase regions, nor is there evidence of anything akin to the giant pulses noted in the radio. Finally, there is no evidence to support the existence of a reported 60 second modulation suggested to be as a consequence of free precession.

Based on observations of the Crab pulsar using the TRIFFID high speed imaging photometer in the UBV bands using the Special Astrophysical Observatory's 6m telescope in the Russian Caucasus, we report the detection of pronounced emission during the so-called `off' phase of emission. Following de-extinction, this unpulsed component of emission is shown to be consistent with a power law with an exponent of alpha = -0.60 +/- 0.37, the uncertainty being dominated by the error associated with the independent CCD photometry used to reference the TRIFFID data. This suggests a steeper power law form than that reported elsewhere in the literature for the total integrated spectrum, which is essentially flat with alpha ~ 0.1, although the difference in this case is only significant at the ~ 2 sigma level. Deeper reference integrated and TRIFFID phase-resolved photometry in these bands in conjunction with further observations in the UV and R region would constrain this fit further.

Observations of the 5 confirmed optical pulsars indicate that the peak emission scales according to the outer field strength. We show that this gives gives further confirmation that a simple phenomenological models such as Pacini and Salvati (1987) still have validity. Furthermore we show that the Crab pulsar exhibits unpulsed emission which further complicates any studies of the thermal emission from Isolated Neutron Stars.

Using the high speed 2-d TRIFFID photometer, we have obtained phase resolved photometry of the Crab pulsar in UBV that allows us to flux the unpulsed light curve component. Following de-extinction, weighted least-square fitting indicates a power-law exponent of alpha = -0.62 +/- 0.49. This is steeper than that reported for the peak components and its origin remains unclear with respect to contemporary magnetospheric theory.

Using time-resolved two-dimensional aperture photometry we have put upper limits on the pulsed emission from two proposed optical counterparts for PSR B1951+32. Our pulsed upper limits of m_{vpulsed}>23.3, m_{bpulsed}>24.4$, for the first candidate and m_{vpulsed}>23.6$, m_{bpulsed}>24.3 for the second, make it unlikely that either of these is, in fact, the pulsar. We discuss three further candidates, but also reject these on the basis of timing results. A search of a 5.5arcs x 5.5arcs area centred close to these stars failed to find any significant pulsations at the reported pulsar period.

We present optical data which shows that G", the optical counterpart of the gamma-ray pulsar Geminga, pulses in B with a period of 0.237 seconds. The similarity between the optical pulse shape and the gamma-ray light curve indicates that a large fraction of the optical emission is non-thermal in origin - contrary to recent suggestions based upon the total optical flux. The derived magnitude of the pulsed emission is m_B = 26.0 +/- 0.4. Whilst it is not possible to give an accurate figure for the pulsed fraction (due to variations in the sky background) we can give an upper limit of m_B ~ 27 for the unpulsed fraction.

The analysis of the light curves of the optical counterparts of GRB 970228 and GRB 970508 points out remarkable similarities. The spectral distribution obtained from the color indices shows that both the transients became bluer during the increasing stage, and redder after the maximum. A main difference concerns the behaviour of the optical fading which is well fitted by a single power law in the case of GRB 970508 but not in the case of GRB 970228.