We obtained thorough photometric observations of two binary near-Earth asteroids (66391) Moshup = 1999 KW4 and (88710) 2001 SL9 from June 2000 to June 2019 and derived physical and dynamical properties of the binary systems. We found that the data for 1999 KW4 are inconsistent with a constant orbital period and we obtained unique solution with a quadratic drift of the mean anomaly of the satellite of -0.65 +/- 0.16 deg/yr2 (all quoted uncertainties are 3sigma). This means that the semimajor axis of the mutual orbit of the components of this binary system increases in time with a mean rate of 1.2 +/- 0.3 cm/yr. The data for 2001 SL9 are also inconsistent with a constant orbital period and we obtained two solutions for the quadratic drift of the mean anomaly: 2.8 +/- 0.2 and 5.2 +/- 0.2 deg/yr2, implying that the semimajor axis of the mutual orbit of the components decreases in time with a mean rate of -2.8 +/- 0.2 or -5.1 +/- 0.2 cm/yr for the two solutions, respectively. The expanding orbit of 1999 KW4 may be explained by mutual tides interplaying with binary YORP (BYORP) effect (McMahon and Scheeres, 2010). However, a modeling of the BYORP drift using radar-derived shapes of the binary components predicted a much higher value of the orbital drift than the observed one. It suggests that either the radar-derived shape model of the secondary is inadequate for computing the BYORP effect, or the present theory of BYORP overestimates it. It is possible that the BYORP coefficient has instead an opposite sign than predicted; in that case, the system may be moving into an equilibrium between the BYORP and the tides. In the case of 2001 SL9, the BYORP effect is the only known physical mechanism that can cause the inward drift of its mutual orbit.

We present the Global Rapid Advanced Network Devoted to the Multi-messenger Addicts (GRANDMA). The network consists of 21 telescopes with both photometric and spectroscopic facilities. They are connected together thanks to a dedicated infrastructure. The network aims at coordinating the observations of large sky position estimates of transient events to enhance their follow-up and reduce the delay between the initial detection and the optical confirmation. The GRANDMA program mainly focuses on follow-up of gravitational-wave alerts to find and characterise the electromagnetic counterpart during the third observational campaign of the Advanced LIGO and Advanced Virgo detectors. But it allows for any follow-up of transient alerts involving neutrinos or gamma-ray bursts, even with poor spatial localisation. We present the different facilities, tools, and methods we developed for this network, and show its efficiency using observations of LIGO/Virgo S190425z, a binary neutron star merger candidate. We furthermore report on all GRANDMA follow-up observations performed during the first six months of the LIGO-Virgo observational campaign, and we derive constraints on the kilonova properties assuming that the events' locations were imaged by our telescopes.

We studied the rotational properties of the dwarf planet Makemake. The photometric observations were carried out at different telescopes between 2006 and 2017. Most of the measurements were acquired in BVRI broad-band filters of a standard Johnson-Cousins photometric system. We found that Makemake rotates more slowly than was previously reported. A possible lightcurve asymmetry suggests a double-peaked period of P = 22.8266$\pm$0.0001~h. A small peak-to-peak lightcurve amplitude in R-filter A = 0.032$\pm$0.005 mag implies an almost spherical shape or near pole-on orientation. We also measured BVRI colours and the R-filter phase-angle slope and revised the absolute magnitudes. The absolute magnitude of Makemake has remained unchanged since its discovery in 2005. No direct evidence of a newly discovered satellite was found in our photometric data; however, we discuss the possible existence of another larger satellite.

We investigate the prompt emission and the afterglow properties of short duration gamma-ray burst (sGRB) 130603B and another eight sGRB events during 2012-2015, observed by several multi-wavelength facilities including the GTC 10.4m telescope. Prompt emission high energy data of the events were obtained by INTEGRAL/SPI/ACS, Swift/BAT and Fermi/GBM satellites. The prompt emission data by INTEGRAL in the energy range of 0.1-10 MeV for sGRB 130603B, sGRB 140606A, sGRB 140930B, sGRB 141212A and sGRB 151228A do not show any signature of the extended emission or precursor activity and their spectral and temporal properties are similar to those seen in case of other short bursts. For sGRB130603B, our new afterglow photometric data constraints the pre jet-break temporal decay due to denser temporal coverage. For sGRB 130603B, the afterglow light curve, containing both our new as well as previously published photometric data is broadly consistent with the ISM afterglow model. Modeling of the host galaxies of sGRB 130603B and sGRB 141212A using the LePHARE software supports a scenario in which the environment of the burst is undergoing moderate star formation activity. From the inclusion of our late-time data for 8 other sGRBs we are able to; place tight constraints on the non-detection of the afterglow, host galaxy or any underlying kilonova emission. Our late-time afterglow observations of the sGRB 170817A/GW170817 are also discussed and compared with the sub-set of sGRBs.

We studied 93 asteroid pairs. We estimated times elapsed since separation of pair members that are between 7*10^3 and a few 10^6 yr. We derived the rotation periods for all the primaries and a sample of secondaries. We derived the absolute magnitude differences of the asteroid pairs that provide their mass ratios. We refined their WISE geometric albedos and estimated their taxonomic classifications. For 17 pairs, we determined their pole positions. In 2 pairs where we obtained the spin poles for both components, we saw the same sense of rotation for both components and constrained the angles between their original spin vectors at the time of their separation. We found that the primaries of 13 pairs are actually binary or triple systems, i.e., they have one or two bound secondaries (satellites). As by-product, we found 3 new young asteroid clusters (each of them consisting of three known asteroids on highly similar orbits). We compared the obtained asteroid pair data with theoretical predictions and discussed their implications. We found that 86 of the 93 studied pairs follow the trend of primary rotation period vs mass ratio that was found by Pravec et al. (2010). Of the 7 outliers, 3 appear insignificant (may be due to our uncertain or incomplete knowledge), but 4 are high mass ratio pairs that were unpredicted by the theory of asteroid pair formation by rotational fission. We discuss a (remotely) possible way that they could be created by rotational fission of flattened parent bodies followed by re-shaping of the formed components. The 13 pairs with binary primaries are particularly interesting systems that place important constraints on formation and evolution of asteroid pairs. We present two hypotheses for their formation: The pairs having both bound and unbound secondaries could be `failed asteroid clusters', or they could be formed by a cascade primary spin fission process.

The rotation states of small asteroids are affected by a net torque arising from an anisotropic sunlight reflection and thermal radiation from the asteroids' surfaces. On long timescales, this so-called YORP effect can change asteroid spin directions and their rotation periods. We analyzed lightcurves of four selected near-Earth asteroids with the aim of detecting secular changes in their rotation rates that are caused by YORP. We use the lightcurve inversion method to model the observed lightcurves and include the change in the rotation rate $\mathrm{d} \omega / \mathrm{d} t$ as a free parameter of optimization. We collected more than 70 new lightcurves. For asteroids Toro and Cacus, we used thermal infrared data from the WISE spacecraft and estimated their size and thermal inertia. We also used the currently available optical and radar astrometry of Toro, Ra-Shalom, and Cacus to infer the Yarkovsky effect. We detected a YORP acceleration of $\mathrm{d}\omega / \mathrm{d} t = (1.9 \pm 0.3) \times 10^{-8}\,\mathrm{rad}\,\mathrm{d}^{-2}$ for asteroid Cacus. For Toro, we have a tentative ($2\sigma$) detection of YORP from a significant improvement of the lightcurve fit for a nonzero value of $\mathrm{d}\omega / \mathrm{d} t = 3.0 \times 10^{-9}\,\mathrm{rad}\,\mathrm{d}^{-2}$. For asteroid Eger, we confirmed the previously published YORP detection with more data and updated the YORP value to $(1.1 \pm 0.5) \times 10^{-8}\,\mathrm{rad}\,\mathrm{d}^{-2}$. We also updated the shape model of asteroid Ra-Shalom and put an upper limit for the change of the rotation rate to $|\mathrm{d}\omega / \mathrm{d} t| \lesssim 1.5 \times 10^{-8}\,\mathrm{rad}\,\mathrm{d}^{-2}$. Ra-Shalom has a greater than $3\sigma$ Yarkovsky detection with a theoretical value consistent with observations assuming its size and/or density is slightly larger than the nominally expected values.

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.

We present optical observations of SN 2013dx, related to the Fermi burst GRB 130702A occurred at a redshift z = 0.145. It is the second-best sampled GRB-SN after SN~1998bw: the observational light curves contain more than 280 data points in uBgrRiz filters until 88 day after the burst, and the data were collected from our observational collaboration (Maidanak Observatory, Abastumani Observatory, Crimean Astrophysical Observatory, Mondy Observatory, National Observatory of Turkey, Observatorio del Roque de los Muchachos) and from the literature. We model numerically the multicolour light curves using the one-dimensional radiation hydrodynamical code STELLA, previously widely implemented for the modelling of typical non-GRB SNe. The best-fitted model has the following parameters: pre-supernova star mass M = 25 M_Sun, mass of a compact remnant M_CR = 6 M_Sun, total energy of the outburst E_oburst = 3.5 x 10^(52) erg, pre-supernova star radius R = 100 R_Sun, M_56Ni = 0.2 M_Sun which is totally mixed through the ejecta; M_O = 16.6 M_Sun, M_Si = 1.2 M_Sun, and M_Fe = 1.2 M_Sun, and the radiative efficiency of the SN is 0.1 per cent.

How black holes accrete surrounding matter is a fundamental, yet unsolved question in astrophysics. It is generally believed that matter is absorbed into black holes via accretion disks, the state of which depends primarily on the mass-accretion rate. When this rate approaches the critical rate (the Eddington limit), thermal instability is supposed to occur in the inner disc, causing repetitive patterns of large-amplitude X-ray variability (oscillations) on timescales of minutes to hours. In fact, such oscillations have been observed only in sources with a high mass accretion rate, such as GRS 1915+105. These large-amplitude, relatively slow timescale, phenomena are thought to have physical origins distinct from X-ray or optical variations with small amplitudes and fast ($\lesssim$10 sec) timescales often observed in other black hole binaries (e.g., XTE J1118+480 and GX 339-4). Here we report an extensive multi-colour optical photometric data set of V404 Cygni, an X-ray transient source containing a black hole of nine solar masses (and a conpanion star) at a distance of 2.4 kiloparsecs. Our data show that optical oscillations on timescales of 100 seconds to 2.5 hours can occur at mass-accretion rates more than ten times lower than previously thought. This suggests that the accretion rate is not the critical parameter for inducing inner-disc instabilities. Instead, we propose that a long orbital period is a key condition for these large-amplitude oscillations, because the outer part of the large disc in binaries with long orbital periods will have surface densities too low to maintain sustained mass accretion to the inner part of the disc. The lack of sustained accretion -- not the actual rate -- would then be the critical factor causing large-amplitude oscillations in long-period systems.

Astrometric observations performed by the Gaia Follow-Up Network for Solar System Objects (Gaia-FUN-SSO) play a key role in ensuring that moving objects first detected by ESA's Gaia mission remain recoverable after their discovery. An observation campaign on the potentially hazardous asteroid (99 942) Apophis was conducted during the asteroid's latest period of visibility, from 12/21/2012 to 5/2/2013, to test the coordination and evaluate the overall performance of the Gaia-FUN-SSO . The 2732 high quality astrometric observations acquired during the Gaia-FUN-SSO campaign were reduced with the Platform for Reduction of Astronomical Images Automatically (PRAIA), using the USNO CCD Astrograph Catalogue 4 (UCAC4) as a reference. The astrometric reduction process and the precision of the newly obtained measurements are discussed. We compare the residuals of astrometric observations that we obtained using this reduction process to data sets that were individually reduced by observers and accepted by the Minor Planet Center. We obtained 2103 previously unpublished astrometric positions and provide these to the scientific community. Using these data we show that our reduction of this astrometric campaign with a reliable stellar catalog substantially improves the quality of the astrometric results. We present evidence that the new data will help to reduce the orbit uncertainty of Apophis during its close approach in 2029. We show that uncertainties due to geolocations of observing stations, as well as rounding of astrometric data can introduce an unnecessary degradation in the quality of the resulting astrometric positions. Finally, we discuss the impact of our campaign reduction on the recovery process of newly discovered asteroids.

Using our photometric observations taken between 1996 and 2013 and other published data, we derived properties of the binary near-Earth asteroid (175706) 1996 FG3 including new measurements constraining evolution of the mutual orbit with potential consequences for the entire binary asteroid population. We also refined previously determined values of parameters of both components, making 1996 FG3 one of the most well understood binary asteroid systems. We determined the orbital vector with a substantially greater accuracy than before and we also placed constraints on a stability of the orbit. Specifically, the ecliptic longitude and latitude of the orbital pole are 266{\deg} and -83{\deg}, respectively, with the mean radius of the uncertainty area of 4{\deg}, and the orbital period is 16.1508 +/- 0.0002 h (all quoted uncertainties correspond to 3sigma). We looked for a quadratic drift of the mean anomaly of the satellite and obtained a value of 0.04 +/- 0.20 deg/yr^2, i.e., consistent with zero. The drift is substantially lower than predicted by the pure binary YORP (BYORP) theory of McMahon and Scheeres (McMahon, J., Scheeres, D. [2010]. Icarus 209, 494-509) and it is consistent with the theory of an equilibrium between BYORP and tidal torques for synchronous binary asteroids as proposed by Jacobson and Scheeres (Jacobson, S.A., Scheeres, D. [2011]. ApJ Letters, 736, L19). Based on the assumption of equilibrium, we derived a ratio of the quality factor and tidal Love number of Q/k = 2.4 x 10^5 uncertain by a factor of five. We also derived a product of the rigidity and quality factor of mu Q = 1.3 x 10^7 Pa using the theory that assumes an elastic response of the asteroid material to the tidal forces. This very low value indicates that the primary of 1996 FG3 is a 'rubble pile', and it also calls for a re-thinking of the tidal energy dissipation in close asteroid binary systems.

Events such as GRB130606A at z=5.91, offer an exciting new window into pre-galactic metal enrichment in these very high redshift host galaxies. We study the environment and host galaxy of GRB 130606A, a high-z event, in the context of a high redshift population of GRBs. We have obtained multiwavelength observations from radio to gamma-ray, concentrating particularly on the X-ray evolution as well as the optical photometric and spectroscopic data analysis. With an initial Lorentz bulk factor in the range \Gamma_0 ~ 65-220, the X-ray afterglow evolution can be explained by a time-dependent photoionization of the local circumburst medium, within a compact and dense environment. The host galaxy is a sub-DLA (log N (HI) = 19.85+/-0.15), with a metallicity content in the range from ~1/7 to ~1/60 of solar. Highly ionized species (N V and Si IV) are also detected. This is the second highest redshift burst with a measured GRB-DLA metallicity and only the third GRB absorber with sub-DLA HI column density. GRB ' lighthouses' therefore offer enormous potential as backlighting sources to probe the ionization and metal enrichment state of the IGM at very high redshifts for the chemical signature of the first generation of stars.

The spin state of small asteroids can change on a long timescale by the Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effect, the net torque that arises from anisotropically scattered sunlight and proper thermal radiation from an irregularly-shaped asteroid. The secular change in the rotation period caused by the YORP effect can be detected by analysis of asteroid photometric lightcurves. We analyzed photometric lightcurves of near-Earth asteroids (1865) Cerberus, (2100) Ra-Shalom, and (3103) Eger with the aim to detect possible deviations from the constant rotation caused by the YORP effect. We carried out new photometric observations of the three asteroids, combined the new lightcurves with archived data, and used the lightcurve inversion method to model the asteroid shape, pole direction, and rotation rate. The YORP effect was modeled as a linear change in the rotation rate in time d\omega /dt. Values of d\omega/ dt derived from observations were compared with the values predicted by theory. We derived physical models for all three asteroids. We had to model Eger as a nonconvex body because the convex model failed to fit the lightcurves observed at high phase angles. We probably detected the acceleration of the rotation rate of Eger d\omega / dt = (1.4 +/- 0.6) x 10^{-8} rad/d (3\sigma error), which corresponds to a decrease in the rotation period by 4.2 ms/yr. The photometry of Cerberus and Ra-Shalom was consistent with a constant-period model, and no secular change in the spin rate was detected. We could only constrain maximum values of |d\omega / dt| < 8 x 10^{-9} rad/d for Cerberus, and |d\omega / dt| < 3 x 10^{-8} rad/d for Ra-Shalom.

A new method for the determination of open cluster membership based on a cumulative effect is proposed. In the field of a plate the relative x and y coordinate positions of each star with respect to all the other stars are added. The procedure is carried out for two epochs t_1 and t_2 separately, then one sum is subtracted from another. For a field star the differences in its relative coordinate positions of two epochs will be accumulated. For a cluster star, on the contrary, the changes in relative positions of cluster members at t_1 and t_2 will be very small. On the histogram of sums the cluster stars will gather to the left of the diagram, while the field stars will form a tail to the right. The procedure allows us to efficiently discriminate one group from another. The greater the distance between t_1 and t_2 and the more cluster stars present, the greater is the effect. The accumulation method does not require reference stars, determination of centroids and modelling the distribution of field stars, necessary in traditional methods. By the proposed method 240 open clusters have been processed, including stars up to m<13. The membership probabilities have been calculated and compared to those obtained by the most commonly used Vasilevskis-Sanders method. The similarity of the results acquired the two different approaches is satisfactory for the majority of clusters.