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.

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 used the TAROT network of telescopes to search for the electromagnetic counterparts of GW150914, GW170104 and GW170814, which were reported to originate from binary black hole merger events by the LIGO and Virgo collaborations. Our goal is to constrain the emission from a binary black hole coalescence at visible wavelengths. We developed a simple and effective algorithm to detect new sources by matching the image data with the Gaia catalog data release 1. Machine learning was used and an algorithm was designed to locate unknown sources in a large field of view image. The angular distance between objects in the image and in the catalog was used to find new sources; we then process the candidates to validate them as possible new unknown celestial objects. Though several possible candidates were detected in the three gravitational wave source error boxes studied, none of them were confirmed as a viable counterpart. The algorithm was effective for the identification of unknown candidates in a very large field and provided candidates for GW150914, GW170104 and GW170814. The entire 90\% GW170814 error box was surveyed extensively within 0.6 days after the GW emission resulting in an absolute limiting R magnitude of -23.8. This strong limit excludes to a great extent a possible emission of a gamma-ray burst with an optical counterpart associated with GW170814.

The discovery of the first electromagnetic counterpart to a gravitational wave signal has generated follow-up observations by over 50 facilities world-wide, ushering in the new era of multi-messenger astronomy. In this paper, we present follow-up observations of the gravitational wave event GW170817 and its electromagnetic counterpart SSS17a/DLT17ck (IAU label AT2017gfo) by 14 Australian telescopes and partner observatories as part of Australian-based and Australian-led research programs. We report early- to late-time multi-wavelength observations, including optical imaging and spectroscopy, mid-infrared imaging, radio imaging, and searches for fast radio bursts. Our optical spectra reveal that the transient source afterglow cooled from approximately 6400K to 2100K over a 7-day period and produced no significant optical emission lines. The spectral profiles, cooling rate, and photometric light curves are consistent with the expected outburst and subsequent processes of a binary neutron star merger. Star formation in the host galaxy probably ceased at least a Gyr ago, although there is evidence for a galaxy merger. Binary pulsars with short (100 Myr) decay times are therefore unlikely progenitors, but pulsars like PSR B1534+12 with its 2.7 Gyr coalescence time could produce such a merger. The displacement (about 2.2 kpc) of the binary star system from the centre of the main galaxy is not unusual for stars in the host galaxy or stars originating in the merging galaxy, and therefore any constraints on the kick velocity imparted to the progenitor are poor.