The simultaneous detection of electromagnetic and gravitational waves from the coalescence of two neutron stars (GW170817 and GRB170817A) has ushered in a new era of "multi-messenger" astronomy, with electromagnetic detections spanning from gamma to radio. This great opportunity for new scientific investigations raises the issue of how the available multi-messenger tools can best be integrated to constitute a powerful method to study the transient universe in particular. To facilitate the classification of possible optical counterparts to gravitational-wave events, it is important to optimize the scheduling of observations and the filtering of transients, both key elements of the follow-up process. In this work, we describe the existing workflow whereby telescope networks such as GRANDMA and GROWTH are currently scheduled; we then present modifications we have developed for the scheduling process specifically, so as to face the relevant challenges that have appeared during the latest observing run of Advanced LIGO and Advanced Virgo. We address issues with scheduling more than one epoch for multiple fields within a skymap, especially for large and disjointed localizations. This is done in two ways: by optimizing the maximum number of fields that can be scheduled, and by splitting up the lobes within the skymap by right ascension to be scheduled individually. In addition, we implement the ability to take previously observed fields into consideration when rescheduling. We show the improvements that these modifications produce in making the search for optical counterparts more efficient, and we point to areas needing further improvement.

The ever-increasing sensitivity of the network of gravitational-wave detectors has resulted in the accelerated rate of detections from compact binary coalescence systems in the third observing run of Advanced LIGO and Advanced Virgo. Not only has the event rate increased, but also the distances to which phenomena can be detected, leading to a rise in the required sky volume coverage to search for counterparts. Additionally, the improvement of the detectors has resulted in the discovery of more compact binary mergers involving neutron stars, revitalizing dedicated follow-up campaigns. While significant effort has been made by the community to optimize single telescope observations, using both synoptic and galaxy-targeting methods, less effort has been paid to coordinated observations in a network. This is becoming crucial, as the advent of gravitational-wave astronomy has garnered interest around the globe, resulting in abundant networks of telescopes available to search for counterparts. In this paper, we extend some of the techniques developed for single telescopes to a telescope network. We describe simple modifications to these algorithms and demonstrate them on existing network examples. These algorithms are implemented in the open-source software \texttt{gwemopt}, used by some follow-up teams, for ease of use by the broader community.

Ever since the insightful analysis of the durations of GRBs by [Kouveliotou:93], GRBs have most often been classified into two populations: "short bursts" (shorter than 2.0 seconds) and "long bursts" (longer than 2.0 seconds). However, recent works have suggested the existence of an intermediate population in the bursts observed by the Swift satellite. Moreover, some researchers have questioned the universality of the 2.0-second dividing line between short and long bursts: some bursts may be short but actually result from collapsars, the physical mechanism behind normally long bursts, and some long ones may originate from mergers, the usual progenitors of short GRBs. In this work, we focus on GRBs detected by the Fermi satellite and study the distribution of their durations measured in the observer's reference frame and, for those with known redshifts, in the bursts' reference frames. However, there are relatively few bursts with measured redshifts, and this makes an accurate study difficult. To overcome this problem, we follow [Zhang:2018] and determine a "pseudo-redshift" from the correlation relation between the luminosity $L_p$ and the energy $E_p$, both of which are calculated at the peak of the flux. Interestingly, we find that the uncertainties in the quantities observed and used in the determination of pseudo-redshifts, do affect the precision of the individual results significantly, but they keep the distribution of pseudo-redshifts very similar to that of the actual ones and thus allow us to use pseudo-redshifts for our statistical study. We use the reduced chi-square and the maximization of the log-likelihood to statistically analyze the distribution of Fermi GRB durations. Both methods show that the distribution of the observed and the intrinsic bursts durations are better represented by two groups, rather than three.