Joint multi-messenger observations with gravitational waves and electromagnetic data offer new insights into the astrophysical studies of compact objects. The third Advanced LIGO and Advanced Virgo observing run began on April 1, 2019; during the eleven months of observation, there have been 14 compact binary systems candidates for which at least one component is potentially a neutron star. Although intensive follow-up campaigns involving tens of ground and space-based observatories searched for counterparts, no electromagnetic counterpart has been detected. Following on a previous study of the first six months of the campaign, we present in this paper the next five months of the campaign from October 2019 to March 2020. We highlight two neutron star - black hole candidates (S191205ah, S200105ae), two binary neutron star candidates (S191213g and S200213t) and a binary merger with a possible neutron star and a "MassGap" component, S200115j. Assuming that the gravitational-wave candidates are of astrophysical origin and their location was covered by optical telescopes, we derive possible constraints on the matter ejected during the events based on the non-detection of counterparts. We find that the follow-up observations during the second half of the third observing run did not meet the necessary sensitivity to constrain the source properties of the potential gravitational-wave candidate. Consequently, we suggest that different strategies have to be used to allow a better usage of the available telescope time. We examine different choices for follow-up surveys to optimize sky localization coverage vs.\ observational depth to understand the likelihood of counterpart detection.

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.

Short gamma-ray bursts result from mergers of two neutron stars or from collapsars, but probably at a smaller rate. In the first case, a kilonova occurs while in the second case a Type Ic supernova is expected. } {Even if future observations of kilonovae in association with gravitational wave events provide better data, detecting a kilonova during an afterglow follow-up would remain useful for exploring the diversity of the kilonova phenomenon. As supernovae produce a weaker gravitational signal, afterglow follow-up will be the only possible method to find one. In this work, we identify the conditions of the burst energy, external density, kilonova mass, supernova luminosity, that are necessary for the detection of a kilonova or supernova in the follow-up of short GRB afterglows.} {We have used a simple kilonova model to obtain the peak luminosities and times as a function of mass, expansion velocity and ejected matter opacity. Afterglow light curves are computed for a uniform medium and a stellar wind, in the kilonova and supernova cases, respectively.} {We represent, using diagrams of the burst kinetic energy vs. density of the external medium, the domains where the kilonova or supernova at maximum is brighter than the afterglow. %\LEt{Please check I have retained your intended meaning -> It's OK

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.

We use the correlation relation between the energy emitted by the GRBs in their prompt phases and the X-ray afterglow fluxes, in an effort to constrain cosmological parameters and aiming to construct a Hubble diagram at high redshifts, i.e. beyond those found with Type Ia supernovae. We use a sample of 126 \textit{Swift} GRBs, that we have selected among more than 800 long bursts observed until April 2015. The selection is based on a few observational constraints: GRB flux higher than 0.4 $photons/cm^2/s$ in the band 15-150 keV, spectrum fitted with simple power law, redshift accurately known and given, and X-ray afterglow observed and flux measured. The statistical method of maximum likelihood is then used to determine the best cosmological parameters ($\Omega_M$, $\Omega_{\Lambda}$) that give the best correlation for two relations: a) the Amati relation (between intrinsic spectral peak energy $E_{p,i}$ and the equivalent isotropic energy), b) the Dainotti relation, namely between the X-ray afterglow luminosity $L_X$ and the break time $T_a$, which is observed in the X-ray flux FX. Although the number of GRBs with high redshifts is rather small, and despite the notable dispersion found in the data, the results we have obtained are quite encouraging and promising. The results obtained using the Amati relation are close to those obtained using the Type Ia supernovae, and they appear to indicate a universe dominated by dark energy. However, those obtained with the correlation between the break time and the X-ray afterglow luminosity is consistent with the findings of the WMAP study of the cosmic microwave background radiation, and they seem to indicate a de Sitter-Einstein universe dominated by matter.

Studies of \textit{BATSE} bursts \citep{kouveliotou:93} have resulted in the widespread adoption of a two-group categorization: long bursts (those with durations $\geq 2$ seconds) and short bursts (those with durations $\leq 2$ seconds). This categorization, one must recall, used the observed $T_{90}$ time durations for bursts (during which 90\% of a burst's fluence is measured). In this work, we have explored two ideas: 1) a statistical search for a possible third, intermediate category of bursts (between the "short" and the "long" ones) among 2041 \textit{BATSE} GRBs and 757 \textit{Swift/BAT} ones, 2) a study of bursts' intrinsic durations, where durations in the bursts' reference frames (instead of the observed durations) are considered, for this, 248 \textit{Swift/BAT} bursts that have redshift measurements were statistically analyzed for the same categorization goal. We first use a Monte Carlo method to determine the proper binning of each GRB, considering that bursts come with different uncertainties on their durations. Then, using the method of minimization of chi-square $\chi^2$, we search for the best fit of the normalized frequency distributions $\frac{1}{N_0}\frac{dN}{d\ln{T}}$ of durations, this allows us to compare fits with two groups ("short" and "long") with fits with three groups ("short", "long", and "intermediate"). Our results indicate that the distributions of observed durations are better fitted by three groups than two groups for \textit{Swift/BAT} data, interestingly, the "intermediate" group appears rather clearly for both observed and intrinsic durations. For BATSE data, the statistical test does not prefer three groups over two. We discuss the results, their possible underlying causes, and reasonable interpretations.

We use a sample of \textit{Swift} gamma-ray bursts (GRBs) to analyze the Amati and Yonetoku correlations. The first relation is between $E_{p,i}$, the intrinsic peak energy of the prompt GRB emission, and $E_{iso}$, the equivalent isotropic energy. The second relation is between $E_{p,i}$ and $L_{iso}$, the isotropic peak luminosity. We select a sample of 71 \textit{Swift} GRBs that have a measured redshift and whose observed $E^{obs}_p$ is within the interval of energy 15-150 keV with a relative uncertainty of less than 70\%. We seek to find correlation relations for long-duration GRBs (LGRBs) with a peak photon flux $P_{ph}\geq 2.6~ \mathrm{ph/cm^{2}/s}$. Uncertainties (error bars) on the values of the calculated energy flux \textit{P}, the energy $E_{iso}$, and the peak isotropic luminosity $L_{iso}$ are estimated using a Monte Carlo approach. We find 27 \textit{Swift} LGRBs that satisfy all our constraints. Results of our analyses of the sample of 71 GRBs and the selected subsample (27 GRBs) are in good agreement with published results. The plots of the two relations for all bursts show a large dispersion around the best straight lines in the sample of 71 LGRBs but not so much in the subsample of 27 GRBs.

The first gamma-ray line originating from outside the solar system that was ever detected is the 511 keV emission from positron annihilation in the Galaxy. Despite 30 years of intense theoretical and observational investigation, the main sources of positrons have not been identified up to now. Observations in the 1990's with OSSE/CGRO showed that the emission is strongly concentrated towards the Galactic bulge. In the 2000's, the SPI instrument aboard ESA's INTEGRAL gamma-ray observatory allowed scientists to measure that emission across the entire Galaxy, revealing that the bulge/disk luminosity ratio is larger than observed in any other wavelength. This mapping prompted a number of novel explanations, including rather "exotic ones (e.g. dark matter annihilation). However, conventional astrophysical sources, like type Ia supernovae, microquasars or X-ray binaries, are still plausible candidates for a large fraction of the observed total 511 keV emission of the bulge. A closer study of the subject reveals new layers of complexity, since positrons may propagate far away from their production sites, making it difficult to infer the underlying source distribution from the observed map of 511 keV emission. However, contrary to the rather well understood propagation of high energy (>GeV) particles of Galactic cosmic rays, understanding the propagation of low energy (~MeV) positrons in the turbulent, magnetized interstellar medium, still remains a formidable challenge. We review the spectral and imaging properties of the observed 511 keV emission and we critically discuss candidate positron sources and models of positron propagation in the Galaxy.

We examine the annihilation of positrons on polycyclic aromatic hydrocarbon (PAH) molecules in interstellar medium conditions. We estimate the annihilation rates of positrons on PAHs by a semi-empirical approach. We show that PAHs can play a significant role in the overall galactic positron annihilation picture and use the annihilation rates and INTEGRAL galactic emission measurements to constrain the amount of PAHs present in the ISM. We find an upper limit of 4.6 x 10^-7 for the PAH abundance.

We consider the production of positrons in microquasars, i.e. X-ray binary systems that exhibit jets frequently, but not continuously. We estimate the production rate of positrons in microquasars, both by simple energy considerations and in the framework of various proposed models. We then evaluate the collective emissivity of the annihilation radiation produced by Galactic microquasars and we find that it might constitute a substantial contribution to the annihilation flux measured by INTEGRAL/SPI. We also discuss the possible spatial distribution of Galactic microquasars, on the basis of the (scarce) available data and the resulting morphology of the flux received on Earth. Finally, we consider nearby 'misaligned' microquasars, with jets occasionally hitting the atmosphere of the companion star; these would represent interesting point sources, for which we determine the annihilation flux and the corresponding light curve, as well as the line's spectral profile. We discuss the possibility of detection of such point sources by future instruments.

We present a measurement of the sky distribution of positronium (Ps) annihilation continuum emission obtained with the SPI spectrometer on board ESA's INTEGRAL observatory. The only sky region from which significant Ps continuum emission is detected is the Galactic bulge. The Ps continuum emission is circularly symmetric about the Galactic centre, with an extension of about 8 deg FWHM. Within measurement uncertainties, the sky distribution of the Ps continuum emission is consistent with that found by us for the 511 keV electron-positron annihilation line using SPI. Assuming that 511 keV line and Ps continuum emission follow the same spatial distribution, we derive a Ps fraction of 0.92 +/- 0.09. These results strengthen our conclusions regarding the origin of positrons in our Galaxy based on observations of the 511 keV line. In particular, they suggest that the main source of Galactic positrons is associated with an old stellar population, such as Type Ia supernovae, classical novae, or low-mass X-ray binaries. Light dark matter is a possible alternative source of positrons.

We present a spectral analysis of the e+e- annihilation emission from the Galactic Centre region based on the first year of measurements made with the spectrometer SPI of the INTEGRAL mission. We have found that the annihilation spectrum can be modelled by the sum of a narrow and a broad 511 keV line plus an ortho-Ps continuum. The broad line is detected with a flux of (0.35+/-0.11)e-3 s-1 cm-2. The measured width of 5.4+/-1.2 keV FWHM is in agreement with the expected broadening of 511 keV photons emitted in the annihilation of Ps that are formed by the charge exchange process of slowing down positrons with H atoms. The flux of the narrow line is (0.72+/-0.12)e-3 s-1 cm-2 and its width is 1.3+/-0.4 keV FWHM. The measured ortho-Ps continuum flux yields a fraction of Ps of (96.7+/-2.2)%. To derive in what phase of the interstellar medium positrons annihilate, we have fitted annihilation models calculated for each phase to the data. We have found that 49(+2,-23)% of the annihilation emission comes from the warm neutral phase and 51(+3,-2)% from the warm ionized phase. While we may not exclude that less than 23% of the emission might come from cold gas, we have constrained the fraction of annihilation emission from molecular clouds and hot gas to be less than 8% and 0.5%, respectively. We have compared our knowledge of the interstellar medium in the bulge and the propagation of positrons with our results and found that they are in good agreement if the sources are diffusively distributed and if the initial kinetic energy of positrons is lower than a few MeV. Despite its large filling factor, the lack of annihilation emission from the hot gas is due to its low density, which allows positrons to escape this phase.

We present a map of 511 keV electron-positron annihilation emission, based on data accumulated with the SPI spectrometer aboard ESA's INTEGRAL gamma-ray observatory, that covers approximately 95% of the celestial sphere. 511 keV line emission is significantly detected towards the galactic bulge region and, at a very low level, from the galactic disk. The bulge emission is highly symmetric and is centred on the galactic centre with an extension of 8 deg. The emission is equally well described by models that represent the stellar bulge or halo populations. The disk morphology is only weakly constrained by the present data, being compatible with both the distribution of young and old stellar populations. The 511 keV line flux from the bulge and disk components is 1.05e-3 ph cm-2 s-1 and 0.7e-3 ph cm-2 s-1, respectively, corresponding to a bulge-to-disk flux ratio in the range 1-3. Assuming a positronium fraction of 0.93 this translates into annihilation rates of 1.5e43 s-1 and 3e42 s-1, respectively. The ratio of the bulge luminosity to that of the disk is in the range 3-9. We find no evidence for a point-like source in addition to the diffuse emission, down to a typical flux limit of 1e-4 ph cm-2 s-1. We also find no evidence for the positive latitude enhancement that has been reported from OSSE measurements; our 3 sigma upper flux limit for this feature is 1.5e-4 ph cm-2 s-1. The disk emission can be attributed to the beta+ decay of the radioactive species 26Al and 44Ti. The bulge emission arises from a different source which has only a weak or no disk component. We suggest that Type Ia supernovae and/or low-mass X-ray binaries are the prime candidates for the source of the galactic bulge positrons. Light dark matter annihilation could also explain the observed 511 keV bulge emission characteristics.

We reexamine in detail the various processes undergone by positrons in the ISM from their birth to their annihilation using the most recent results of positron interaction cross sections with H, H2 and He. The positrons' lives are divided into two phases: the 'in-flight' phase and the thermal phase. The first phase is treated with a Monte Carlo simulation that allows us to determine the fraction of positrons that form positronium and annihilate as well as the characteristics of the annihilation emission as a function of the medium conditions. The second phase is treated with a binary reaction rate approach, with cross sections adopted from experimental measurement or theoretical calculations. An extensive search and update of the knowledge of positron processes was thus undertaken. New reaction rates and line widths have been obtained. We investigate the treatment of the complicated interactions between positrons and interstellar dust grains. New reaction rates and widths of the line resulting from the annihilation inside and outside of the grain have been obtained. The final results of our calculations showed that dust is only important in the hot phase of the ISM, where it dominates all other processes. Combining the new calculations, we have constructed annihilation spectra for each phase of the ISM, considering various grain contents, as well as an overall combined spectrum for the ISM as a whole.

The status of the analysis of the electron/positron annihilation radiation performed by INTEGRAL/SPI in the 4th GALACTIC QUADRANT is reported. We use data from the first two Galactic Center Deep Exposures (GCDE) and from the Galactic Plane Scans (GPS). The analysis presented here is focused on the spectroscopic aspects of the electron/positron annihilation radiation. Background substraction and model fitting methods are described, and the parameters of the 511 keV line (flux, energy, and width) are deduced.

During its first year in orbit the INTEGRAL observatory performed deep exposures of the Galactic Center region and scanning observations of the Galactic plane. We report on the status of our analysis of the positron annihilation radiation from the 4th Galactic quadrant with the spectrometer SPI, focusing on the sky distribution of the 511 keV line emission. The analysis methods are described; current constraints and limits on the Galactic bulge emission and the bulge-to-disk ratio are presented.

We consider the production of 2.22 MeV radiation resulting from the capture of neutrons in the atmosphere of the secondary in an X-ray binary system, where the neutrons are produced in the accretion disk around the compact primary star and radiated in all directions. We have considered several accretion disk models (ADAF, ADIOS, SLE, Uniform-Temperature) and a varity of parameters (accretion rate, mass of the compact object, mass, temperature and composition of the secondary star, distance between the two objects, etc.). The neutron rates are calculated by a network of nuclear reactions in the accretion disk, and this is handled by a reaction-rate formulation taking into account the structure equations given by each accretion model. The processes undergone by the neutrons in the atmosphere of the companion star are studied in great detail, including thermalization, elastic and inelastic scatterings, absorption, escape from the surface, decay, and capture by protons. The radiative transfer of the 2.22 MeV photons is treated separately, taking into consideration the composition and density of the star's atmosphere. The final flux of the 2.22 MeV radiation that can be detected from earth is calculated taking into account the distance to the source, the direction of observation with respect to the binary system frame, and the rotation of the source, as this can lead to an observable periodicity in the flux. We produce phasograms of the 2.22 MeV intensity as well as spectra of the line, where rotational Doppler shift effects can lead to changes in the spectra that are measurable by INTEGRAL's spectrometer (SPI).

We examine the processes which could lead to the observed enhancement of Li and possibly other light elements (Be, B) in the companions of a number of X-ray novae. We conclude that one of the most promising mechanisms is the spallation of CNO elements on the surface of the companion induced by the neutron flux produced in the hot accretion flow onto the compact object. Direct production of the observed Li and its deposition onto the dwarf companion seem less likely, mainly because of the possibility of its destruction in the production region itself and difficulties in its deposition associated with the configuration of the companion's magnetic field. We discuss other potential observables of the above scenario.