On February 17 2016, the IceCube real-time neutrino search identified, for the first time, three muon neutrino candidates arriving within 100 s of one another, consistent with coming from the same point in the sky. Such a triplet is expected once every 13.7 years as a random coincidence of background events. However, considering the lifetime of the follow-up program the probability of detecting at least one triplet from atmospheric background is 32%. Follow-up observatories were notified in order to search for an electromagnetic counterpart. Observations were obtained by Swift's X-ray telescope, by ASAS-SN, LCO and MASTER at optical wavelengths, and by VERITAS in the very-high-energy gamma-ray regime. Moreover, the Swift BAT serendipitously observed the location 100 s after the first neutrino was detected, and data from the Fermi LAT and HAWC observatory were analyzed. We present details of the neutrino triplet and the follow-up observations. No likely electromagnetic counterpart was detected, and we discuss the implications of these constraints on candidate neutrino sources such as gamma-ray bursts, core-collapse supernovae and active galactic nucleus flares. This study illustrates the potential of and challenges for future follow-up campaigns.

We report the discovery and multicolor (VRIW) photometry of a rare explosive star MASTER OT J004207.99+405501.1 - a luminous red nova - in the Andromeda galaxy M31N2015-01a. We use our original light curve acquired with identical MASTER Global Robotic Net telescopes in one photometric system: VRI during first 30 days and W (unfiltered) during 70 days. Also we added publishied multicolor photometry data to estimate the mass and energy of the ejected shell, and discuss the likely formation scenarios of outbursts of this type. We propose the interpretation of the explosion, that is consistent with the evolutionary scenario where star merger is a natural stage of the evolution of close-mass stars and may serve as an extra channel for the formation of nova outbursts.

The Advanced LIGO observatory recently reported (Abbott et al., 2016a) the first direct detection of gravitational waves predicted by Einstein (1916). The detection of this event was predicted in 1997 on the basis of the Scenario Machine population synthesis calculations (Lipunov et al., NewA, 2, 43). Now we discuss the parameters of binary black holes and event rates predicted by different scenarios of binary evolution. We give a simple explanation of the big difference between detected black hole masses and the mean black hole masses observed in of X-ray Nova systems. The proximity of the masses of the components of GW150914 is in good agreement with the observed initial mass ratio distribution in massive binary systems, as is used in Scenario Machine calculations for massive binaries.

On 2015 June 15 the Swift space observatory discovered that the Galactic black hole candidate V404 Cyg was undergoing another active X-ray phase, after 25 years of inactivity (Barthelmy et al. 2015). Twelve telescopes of the MASTER Global Robotic Net located at six sites across four continents were the first ground based observatories to start optical monitoring of the microquasar after its gamma-ray wakeup at 18h 34m 09s U.T. on 2015 June 15 (Lipunov et al. 2015). In this paper we report, for the first time, the discovery of variable optical linear polarization, changing by 4-6% over a timescale of approximately 1 h, on two different epochs. We can conclude that the additional variable polarization arisies from the relativistic jet generated by the black hole in V404Cyg. The polarization variability correlates with optical brightness changes, increasing when the flux decreases.

The Advanced LIGO observatory recently reported the first direct detection of gravitational waves predicted by Einstein (1916). We report on the first optical observations of the Gravitational Wave (GW) source GW150914 error region with the Global MASTER Robotic Net. We detected several optical transients, which proved to be unconnected with the GW event. Our result is consistent with the assumption that gravitational waves were produced by a binary black hole merger. The detection of the event confirmed the main prediction of the population synthesis performed with the "Scenario Machine" formulated in Lipunov1997b.

We report on the MASTER Global Robotic Net discovery of an eclipsing binary, MASTER OT J095310.04+335352.8, previously known as unremarkable star TYC 2505-672-1, which displays extreme orbital parameters. The orbital period P=69.1 yr is more than 2.5 times longer than that of epsilon-Aurigae, which is the previous record holder. The light curve is characterized by an extremely deep total eclipse with a depth of more than 4.5 mag, which is symmetrically shaped and has a total duration of 3.5 yrs. The eclipse is essentially gray. The spectra acquired with the Russian 6 m BTA telescope both at minimum and maximum light mainly correspond to an M0-1III--type red giant, but the spectra taken at the bottom of eclipse show small traces of a sufficiently hot source. The observed properties of this system can be better explained as the red giant eclipsed by a large cloud (the disk) of small particles surrounding the invisible secondary companion.

The non-linear dynamics of bending instability and vertical structure of a galactic stellar disc embedded into a spherical halo are studied with N-body numerical modelling. Development of the bending instability in stellar galactic disc is considered as the main factor that increases the disc thickness. Correlation between the disc vertical scale height and the halo-to-disc mass ratio is predicted from the simulations. The method of assessment of the spherical-to-disc mass ratio for edge-on spiral galaxies with small bulges is considered. Modelling of eight edge-on galaxies: NGC 891, NGC 4738, NGC 5170, UGC 6080, UGC 7321, UGC 8286, UGC 9422 and UGC 9556 is performed. Parameters of stellar discs, dark haloes and bulges are estimated. The lower limit of the dark-to-luminous mass ratio in our galaxies is of the order of 1 within the limits of their stellar discs. The dark haloes dominate by mass in the galaxies with very thin stellar discs (NGC 5170, UGC 7321 and UGC 8286).

We present results of numerical modeling made for the galactic stellar disk embedded in the spherical halo. The non-linear dynamics of bending instabilities developed in the disk is studied. The axisymmetrical bending mode is considered as a main factor increasing the disk thickness. The model dependence between the disk thickness and the halo to disk mass ratio $M_h/M_d$ is inferred. The modeling results enable us to estimate the dark halo mass for several edge-on galaxies. The obtained values of $M_h/M_d$ for the galaxies considered are of order 1 - 2.