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.

Quasar microlensing offers a unique opportunity to resolve tiny sources in distant active galactic nuclei and study compact object populations in lensing galaxies. We therefore searched for microlensing-induced variability of the gravitationally lensed quasar QSO 2237+0305 (Einstein Cross) using 4374 optical frames taken with the 2.0 m Liverpool Telescope and the 1.5 m Maidanak Telescope. These $gVrRI$ frames over the 2006$-$2019 period were homogeneously processed to generate accurate long-term multi-band light curves of the four quasar images A-D. Through difference light curves, we found strong microlensing signatures. We then focused on the analytical modelling of two putative caustic-crossing events in image C, finding compelling evidence that this image experienced a double caustic crossing. Additionally, our overall results indicate that a standard accretion disc accounts reasonably well for the brightness profile of UV continuum emission sources and for the growth in source radius when the emission wavelength increases: $R_{\lambda} \propto \lambda^{\alpha}$, $\alpha$ = 1.33 $\pm$ 0.09. However, we caution that numerical microlensing simulations are required before firm conclusions can be reached on the UV emission scenario because the $VRI$-band monitoring during the first caustic crossing and one of our two $\alpha$ indicators lead to a few good solutions with $\alpha \approx$ 1.

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

UX Orionis stars (UXors) are Herbig Ae/Be or T Tauri stars exhibiting sporadic occultation of stellar light by circumstellar dust. GM\,Cephei is such a UXor in the young ($\sim4$~Myr) open cluster Trumpler\,37, showing prominent infrared excess, emission-line spectra, and flare activity. Our photometric monitoring (2008--2018) detects (1)~an $\sim$3.43~day period, likely arising from rotational modulation by surface starspots, (2)~sporadic brightening on time scales of days due to accretion, (3)~irregular minor flux drops due to circumstellar dust extinction, and (4)~major flux drops, each lasting for a couple of months with a recurrence time, though not exactly periodic, of about two years. The star experiences normal reddening by large grains, i.e., redder when dimmer, but exhibits an unusual "blueing" phenomenon in that the star turns blue near brightness minima. The maximum extinction during relatively short (lasting $\leq 50$~days) events, is proportional to the duration, a consequence of varying clump sizes. For longer events, the extinction is independent of duration, suggestive of a transverse string distribution of clumps. Polarization monitoring indicates an optical polarization varying $\sim3\%$--8$\%$, with the level anticorrelated with the slow brightness change. Temporal variation of the unpolarized and polarized light sets constraints on the size and orbital distance of the circumstellar clumps in the interplay with the young star and scattering envelope. These transiting clumps are edge-on manifestations of the ring- or spiral-like structures found recently in young stars with imaging in infrared of scattered light, or in submillimeter of thermalized dust emission.

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 near-Earth asteroid 3200 Phaethon (1983 TB) is an attractive object not only from a scientific viewpoint but also because of JAXA's DESTINY+ target. The rotational lightcurve and spin properties were investigated based on the data obtained in the ground-based observation campaign of Phaethon. We aim to refine the lightcurves and shape model of Phaethon using all available lightcurve datasets obtained via optical observation, as well as our time-series observation data from the 2017 apparition. Using eight 1-2-m telescopes and an optical imager, we acquired the optical lightcurves and derived the spin parameters of Phaethon. We applied the lightcurve inversion method and SAGE algorithm to deduce the convex and non-convex shape model and pole orientations. We analysed the optical lightcurve of Phaethon and derived a synodic and a sidereal rotational period of 3.6039 h, with an axis ratio of a/b = 1.07. The ecliptic longitude (lambda) and latitude (beta) of the pole orientation were determined as (308, -52) and (322, -40) via two independent methods. A non-convex model from the SAGE method, which exhibits a concavity feature, is also presented.

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.

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.

Gamma-ray bursts (GRBs) are the most luminous explosions in the universe, yet the nature and physical properties of their energy sources are far from understood. Very important clues, however, can be inferred by studying the afterglows of these events. We present optical and X-ray observations of GRB 130831A obtained by Swift, Chandra, Skynet, RATIR, Maidanak, ISON, NOT, LT and GTC. This burst shows a steep drop in the X-ray light-curve at $\simeq 10^5$ s after the trigger, with a power-law decay index of $\alpha \sim 6$. Such a rare behaviour cannot be explained by the standard forward shock (FS) model and indicates that the emission, up to the fast decay at $10^5$ s, must be of "internal origin", produced by a dissipation process within an ultrarelativistic outflow. We propose that the source of such an outflow, which must produce the X-ray flux for $\simeq 1$ day in the cosmological rest frame, is a newly born magnetar or black hole. After the drop, the faint X-ray afterglow continues with a much shallower decay. The optical emission, on the other hand, shows no break across the X-ray steep decrease, and the late-time decays of both the X-ray and optical are consistent. Using both the X-ray and optical data, we show that the emission after $\simeq 10^5$ s can be explained well by the FS model. We model our data to derive the kinetic energy of the ejecta and thus measure the efficiency of the central engine of a GRB with emission of internal origin visible for a long time. Furthermore, we break down the energy budget of this GRB into the prompt emission, the late internal dissipation, the kinetic energy of the relativistic ejecta, and compare it with the energy of the associated supernova, SN 2013fu.

We present observations of the dark Gamma-Ray Burst GRB 051008 provided by Swift/BAT, Swift/XRT, Konus-WIND, INTEGRAL/SPI-ACS in the high-energy domain and the Shajn, Swift/UVOT, Tautenburg, NOT, Gemini and Keck I telescopes in the optical and near-infrared bands. The burst was detected only in gamma- and X-rays and neither a prompt optical nor a radio afterglow were detected down to deep limits. We identified the host galaxy of the burst, which is a typical Lyman-break Galaxy (LBG) with R-magnitude of 24.06 +/- 0.10. A redshift of the galaxy of z = 2.77 (-0.20,+0.15) is measured photometrically due to the presence of a clear, strong Lyman-break feature. The host galaxy is a small starburst galaxy with moderate intrinsic extinction (A_V = 0.3 mag) and has a SFR of ~ 60 M_Sun / yr typical for LBGs. It is one of the few cases where a GRB host has been found to be a classical Lyman-break galaxy. Using the redshift we estimate the isotropic-equivalent radiated energy of the burst to be E_iso = (1.15 +/- 0.20) x 10^54 erg. We also provide evidence in favour of the hypothesis that the darkness of GRB 051008 is due to local absorption resulting from a dense circumburst medium.

We present optical and near-infrared (NIR) photometry for three gamma-ray burst supernovae (GRB-SNe): GRB 120729A, GRB 130215A / SN 2013ez and GRB 130831A / SN 2013fu. In the case of GRB 130215A / SN 2013ez, we also present optical spectroscopy at t-t0=16.1 d, which covers rest-frame 3000-6250 Angstroms. Based on Fe II (5169) and Si (II) (6355), our spectrum indicates an unusually low expansion velocity of 4000-6350 km/s, the lowest ever measured for a GRB-SN. Additionally, we determined the brightness and shape of each accompanying SN relative to a template supernova (SN 1998bw), which were used to estimate the amount of nickel produced via nucleosynthesis during each explosion. We find that our derived nickel masses are typical of other GRB-SNe, and greater than those of SNe Ibc that are not associated with GRBs. For GRB 130831A / SN 2013fu, we use our well-sampled R-band light curve (LC) to estimate the amount of ejecta mass and the kinetic energy of the SN, finding that these too are similar to other GRB-SNe. For GRB 130215A, we take advantage of contemporaneous optical/NIR observations to construct an optical/NIR bolometric LC of the afterglow. We fit the bolometric LC with the millisecond magnetar model of Zhang & Meszaros (2001), which considers dipole radiation as a source of energy injection to the forward shock powering the optical/NIR afterglow. Using this model we derive an initial spin period of P=12 ms and a magnetic field of B=1.1 x 10^15 G, which are commensurate with those found for proposed magnetar central engines of other long-duration GRBs.

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.

We study brightness variations in the double lensed quasar UM673 (Q0142-100) with the aim of measuring the time delay between its two images. In the paper we combine our previously published observational data of UM673 obtained during the 2003 - 2005 seasons at the Maidanak Observatory with archival and recently observed Maidanak and CTIO UM673 data. We analyze the V, R and I-band light curves of the A and B images of UM673, which cover ten observational seasons from August 2001 to November 2010. We also analyze the time evolution of the difference in magnitudes between images A and B of UM673 over more than ten years. We find that the quasar exhibits both short-term (with amplitude of \sim 0.1 mag in the R band) and high-amplitude (\sim 0.3 mag) long-term variability on timescales of about several months and several years, respectively. These brightness variations are used to constrain the time delay between the images of UM673. From cross-correlation analysis of the A and B quasar light curves and error analysis we measure the mean time delay and its error of 89 \pm11 days. Given the input time delay of 88 days, the most probable value of the delay that can be recovered from light curves with the same statistical properties as the observed R-band light curves of UM673 is 95{+5/-16}{+14/-29} days (68 and 95 % confidence intervals). Analysis of the V - I color variations and V, R and I-band magnitude differences of the quasar images does not show clear evidence of the microlensing variations between 1998 and 2010.

Main belt asteroids (6070) Rheinland and (54827) 2001NQ8 belong to a small population of couples of bodies which reside on very similar heliocentric orbits. Vokrouhlicky & Nesvorny (2008, AJ 136, 280) promoted a term "asteroid pairs", pointing out their common origin within the past tens to hundreds of ky. Previous attempts to reconstruct the initial configuration of Rheinland and 2001NQ8 at the time of their separation have led to the prediction that Rheinland's rotation should be retrograde. Here we report extensive photometric observations of this asteroid and use the lightcurve inversion technique to directly determine its rotation state and shape. We confirm the retrograde sense of rotation of Rheinland, with obliquity value constrained to be >= 140 deg. The ecliptic longitude of the pole position is not well constrained as yet. The asymmetric behavior of Rheinland's lightcurve reflects a sharp, near-planar edge in our convex shape representation of this asteroid. Our calibrated observations in the red filter also allow us to determine $H_R = 13.68\pm 0.05$ and $G = 0.31\pm 0.05$ values of the H-G system. With the characteristic color index $V-R = 0.49\pm 0.05$ for the S-type asteroids, we thus obtain $H = 14.17\pm 0.07$ for the absolute magnitude of (6070) Rheinland. This a significantly larger value than previously obtained from analysis of the astrometric survey observations. We next use the obliquity constraint for Rheinland to eliminate some degree of uncertainty in the past propagation of its orbit. This is because the sign of the past secular change of its semimajor axis due to the Yarkovsky effect is now constrained. Determination of the rotation state of the secondary component, asteroid (54827) 2001NQ8, is the key element in further constraining the age of the pair and its formation process.

We report the results of our multicolor observations of PG 1115+080 with the 1.5-m telescope of the Maidanak Observatory (Uzbekistan, Central Asia) in 2001-2006. Monitoring data in filter R spanning the 2004, 2005 and 2006 seasons (76 data points) demonstrate distinct brightness variations of the source quasar with the total amplitude of almost 0.4 mag. Our R light curves have shown image C leading B by 16.4d and image (A1+A2) by 12d that is inconsistent with the previous estimates obtained by Schechter et al. in 1997 - 24.7d between B and C and 9.4d between (A1+A2) and C. The new values of time delays in PG 1115+080 must result in larger values for the Hubble constant, thus reducing difference between its estimates taken from the gravitational lenses and with other methods. Also, we analyzed variability of the A2/A1 flux ratio, as well as color changes in the archetypal "fold" lens PG 1115+080. We found the A1/A2 flux ratio to grow during 2001-2006 and to be larger at longer wavelengths. In particular, the A2/A1 flux ratio reached 0.85 in filter I in 2006. We also present evidence that both the A1 and A2 images might have undergone microlensing during 2001-2006, with the descending phase for A1 and initial phase for A2. We find that the A2/A1 flux ratio anomaly in PG 1115 can be well explained both by microlensing and by finite distance of the source quasar from the caustic fold.

We report new estimates of the time delays in the quadruple gravitationally lensed quasar PG1115+080, obtained from the monitoring data in filter R with the 1.5-m telescope at the Maidanak Mountain (Uzbekistan, Central Asia) in 2004-2006. The time delays are 16.4 days between images C and B, and 12 days between C and A1+A2, with image C being leading for both pairs. The only known estimates of the time delays in PG1115 are those based on observations by Schechter et al. (1997) -- 23.7 and 9.4 days between images C and B, C and A1+A2, respectively, as calculated by Schechter et al., and 25 and 13.3 days as revised by Barkana (1997) for the same image components with the use of another method. The new values of time delays in PG 1115+080 may be expected to provide larger estimates of the Hubble constant thus decreasing a diversity between the H_0 estimates taken from gravitationally lensed quasars and with other methods.

We present and analyse new R-band frames of the gravitationally lensed double quasar FBQ 0951+2635. These images were obtained with the 1.5m AZT-22 Telescope at Maidanak (Uzbekistan) in the 2001-2006 period. Previous results in the R band (1999-2001 period) and the new data allow us to discuss the dominant kind of microlensing variability in FBQ 0951+2635. The time evolution of the flux ratio A/B does not favour the continuous production of short-timescale (months) flares in the faintest quasar component B (crossing the central region of the lensing galaxy). Instead of a rapid variability scenario, the observations are consistent with the existence of a long-timescale fluctuation. The flux ratio shows a bump in the 2003-2004 period and a quasi-flat trend in more recent epochs. Apart from the global behaviour of A/B, we study the intra-year variability over the first semester of 2004, which is reasonably well sampled. Short-timescale microlensing is not detected in that period. Additional data in the i band (from new i-band images taken in 2007 with the 2m Liverpool Robotic Telescope at La Palma, Canary Islands) also indicate the absence of short-timescale events in 2007.

The paper has been withdrawn because double checking and comparison with other data sets after the original submission showed that a broken R-band filter at the Maidanak telescope had affected our quasar monitoring observations in the years 2004 and 2005. They had led to partially spurious measurements, hence our original analysis and conclusions are not reliable.

The time delays between the components of a lensed quasar are basic tools to analyze the expansion of the Universe and the structure of the main lens galaxy halo. In this paper, we focus on the variability and time delay of the double system SBS 0909+532A,B as well as the time behaviour of the field stars. We use VR optical observations of SBS 0909+532A,B and the field stars in 2003. The frames were taken at Calar Alto, Maidanak and Wise observatories, and the VR light curves of the field stars and quasar components are derived from aperture and point-spread function fitting methods. We measure the R-band time delay of the system from the chi-square and dispersion techniques and 1000 synthetic light curves based on the observed records. One nearby field star (SBS 0909+532c) is found to be variable, and the other two nearby field stars are non-variable sources. With respect to the quasar components, the R-band records seem more reliable and are more densely populated than the V-band ones. The observed R-band fluctuations permit a pre-conditioned measurement of the time delay. From the chi-square minimization, if we assume that the quasar emission is observed first in B and afterwards in A (in agreement with basic observations of the system and the corresponding predictions), we obtain a delay of - 45 (+ 1)/(- 11) days (95% confidence interval). The dispersion technique leads to a similar delay range. A by-product of the analysis is the determination of a totally corrected flux ratio in the R band (corrected by the time delay and the contamination due to the galaxy light). Our 95% measurement of this ratio (0.575 +/- 0.014 mag) is in excellent agreement with previous results from contaminated fluxes at the same time of observation.

We present new brightness monitoring observations of the 4 components of gravitationally lensed system Q2237+0305, which show detection of an intrinsic quasar brightness fluctuation at a time of subdued microlensing activity, between June 27 and October 12, 2003. These data were used to determine the time delays between the arrivals of the four images. The measured delays are -6, 35, and 2 hours for images B, C and D relative to A, respectively, so they confirm that the long history of brightness monitoring has produced significant detection of microlensing. However the error bars associated with the delays, of order 2 days, are too large to discriminate between competing macro-imaging models. Moreover, our simulations show that for the amplitude of this intrinsic fluctuation and for photometric errors intrinsic to optical monitoring from our 1.5-m telescope or from the OGLE monitoring, a daily sampled brightness record cannot produce reliable lags for model discrimination. We use our simulations to devise a strategy for future delay determination with optical data. Nevertheless, we regard these first estimates to be significant, since they are the first direct measurements of time delays made for this system from ground-based observations in the visual wavelengths. Our results provide the most convincing confirmation of the gravitational-lens nature of Q2237+0305, and give observational justification to the extensive literature which attributes the quasar's previously observed brightness fluctuations to microlensing.

We present the results of a monitoring campaign of the double quasar SBS1520+530 at Maidanak observatory from April 2003 to August 2004. We obtained light curves in V and R filters that show small-amplitude \Delta m~0.1 mag intrinsic variations of the quasar on time scales of about 100 days. The data set is consistent with the previously determined time delay of \Delta t=(130+-3) days by Burud et al. (2002). We find that the time delay corrected magnitude difference between the quasar images is now larger by (0.14+-0.03) mag than during the observations by Burud et al. (2002). This confirms the presence of gravitational microlensing variations in this system.

We have observed the gravitational lens system Q2237+0305 from the Maidanak Observatory over the period from August 2002 to November 2003. Here we report the results of our observations. We implemented a two-stage technique that has been developed specifically for the purpose of gravitational lens image reconstruction. The technique is based on the Tikhonov regularization approach and allows one to obtain astrometric and photometric characteristics of the gravitational lens system. Light curves with 78 data points for the four quasar components are obtained. Slow brightness variations over the observational period are found in all components. Images A, C, D have a tendency to decrease in brightness. Image B does not vary more than 0.05mag. The observations did not reveal evidence for large variations in brightness of the components due to microlensing effects. To provide an overall picture of the photometry behaviour, our data are combined with the Maidanak observations published for 1995 -- 2000.

Photometry of the Q2237+0305gravitational lens in VRI spectral bands with the 1.5-m telescope of the high-altitude Maidanak observatory in 1995-2000 is presented. Monitoring of Q2237+0305 in July-October 2000, made at nearly daily basis, did not reveal rapid (night-to-night and intranight) variations of brightness of the components during this time period. Rather slow changes of magnitudes of the components were observed, such as 0.08 mag fading of B and C components and 0.05 mag brightening of D in R band during July 23 - October 7, 2000. By good luck three nights in 1999 were almost at the time of the strong brightness peak of image C, and approximately in the middle of the ascending slope of the image A brightness peak. The C component was the most blue one in the system in 1998 and 1999, having changed its (V-I) color from 0.56 mag to 0.12 mag since August 1997, while its brightness increased almost 1.2 mag during this time period. The A component behaved similarly between August 1998 and August 2000, having become 0.47 mag brighter in R, and at the same time, 0.15 mag bluer. A correlation between the color variations and variations of magnitudes of the components is demonstrated to be significant and reaches 0.75, with a regression line slope of 0.33. A color (V-I) vrs color (V-R) plot shows the components settled in a cluster, stretched along a line with a slope of 1.31. Both slopes are noticeably smaller than those expected if a standard galactic interstellar reddening law were responsible for the differences between the colors of images and their variations over time. We attribute the brightness and color changes to microlensing of the quasar's structure, which we conclude is more compact at shorter wavelengths, as predicted by most quasar models featuring an energizing central source.