Context. Only several doubly eclipsing quadruple stellar systems are known to date, and no dedicated effort to characterize population properties of these interesting objects has yet been made. Aims. Our first goal was to increase number of known doubly eclipsing systems such that the resulting dataset would allow us to study this category of objects via statistical means. In order to minimize biases, we used long-lasting, homogeneous, and well-documented photometric surveys. Second, a common problem of basically all known doubly eclipsing systems is the lack of proof that they constitute gravitationally bound quadruple system in the 2+2 architecture (as opposed to two unrelated binaries that are projected onto the same location in the sky by chance). When possible, we thus sought evidence for the relative motion of the two binaries. In that case, we tried to determine the relevant orbital periods and other parameters. Methods. We analysed photometric data for eclipsing binaries provided by the OGLE survey and we focused on the LMC fields. We found a large number of new doubly eclipsing systems (our discoveries are three times more numerous than the previously known cases in this dataset). In order to prove relative motion of the binaries about a common centre of mass, we made use of the fact that OGLE photometry covers several years. With a typical orbital period of days for the observed binaries, we sought eclipse time variations (ETVs) on the timescale comparable to a decade (this is the same method used for an archetype of the doubly eclipsing system, namely V994 Her). In the cases where we were able to detect the ETV period, the difference between the inner and outer periods in the quadruple system is large enough. This allows us to interpret ETVs primarily as the light-time effect, thus providing an interesting constraint on masses of the binaries.Results. In addition to significantly enlarging the database of known doubly eclipsing systems, we performed a thorough analysis of 72 cases. ETVs for 28 of them (39% of the studied cases) showed evidence of relative motion. Among these individual systems, we note OGLE BLG-ECL-145467, by far the most interesting case; it is bright (12.6 mag in I filter), consists of two detached binaries with periods of ≃ 3.3 d and ≃ 4.9 d (making it a candidate for a 3 : 2 resonant system) revolving about each other in only ≃ 1538 d. Distribution of the orbital period ratio P A /P B of binaries in 2+2 quadruples shows statistically significant excess at ≃ 1 and ≃ 1.5. The former is likely a natural statistical preference in weakly interacting systems with periods within the same range. The latter is thought to be evidence of a capture in the 3 : 2 mean motion resonance of the two binaries. This sets important constraints on evolutionary channels in these systems. Conclusions. The total number of doubly eclipsing systems increased to 146, more than 90% of which are at low declinations on the southern sky. This motivates us to use southern hemisphere facilities to further character...

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.

No abstract