The identification of electromagnetic emission from gravitational-wave sources typically requires multiple follow-up observations due to the limited fields-of-view of follow-up observatories compared to the poorly localized direction of gravitational waves. Gravitational-wave localization regions are typically covered with multiple telescope pointings using a "honeycomb" structure, which is optimal only on an infinite, flat surface. Here we present a machine-learning algorithm which uses genetic algorithms along with Broyden-Fletcher-Goldfarb-Shanno (BFGS) optimization to find an optimal configuration of tiles to cover the gravitational-wave sky localization area on a spherical surface.

This paper presents EinsteinPy (version 0.3), a community-developed Python package for gravitational and relativistic astrophysics. Python is a free, easy to use a high-level programming language which has seen a huge expansion in the number of its users and developers in recent years. Specifically, a lot of recent studies show that the use of Python in Astrophysics and general physics has increased exponentially. We aim to provide a very high level of abstraction, an easy to use interface and pleasing user experience. EinsteinPy is developed keeping in mind the state of a theoretical gravitational physicist with little or no background in computer programming and trying to work in the field of numerical relativity or trying to use simulations in their research. Currently, EinsteinPy supports simulation of time-like and null geodesics and calculates trajectories in different background geometries some of which are Schwarzschild, Kerr, and KerrNewmann along with coordinate inter-conversion pipeline. It has a partially developed pipeline for plotting and visualization with dependencies on libraries like Plotly, matplotlib, etc. One of the unique features of EinsteinPy is a sufficiently developed symbolic tensor manipulation utilities which are a great tool in itself for teaching yourself tensor algebra which for many beginner students can be overwhelmingly tricky. EinsteinPy also provides few utility functions for hypersurface embedding of Schwarzschild spacetime which further will be extended to model gravitational lensing simulation.

The recent discovery of TeV emission from gamma-ray bursts (GRBs) by the MAGIC and H.E.S.S. Cherenkov telescopes confirmed that emission from these transients can extend to very high energies. The TeV energy domain reaches the most sensitive band of the Cherenkov Telescope Array (CTA). This newly anticipated, improved sensitivity will enhance the prospects of gravitational-wave follow-up observations by CTA to probe particle acceleration and high-energy emission from binary black hole and neutron star mergers, and stellar core-collapse events. Here we discuss the implications of TeV emission on the most promising strategies of choice for the gravitational-wave follow-up effort for CTA and Cherenkov telescopes more broadly. We find that TeV emission (i) may allow more than an hour of delay between the gravitational-wave event and the start of CTA observations; (ii) enables the use of CTA's small size telescopes that have the largest fields of view. We characterize the number of pointings needed to find a counterpart. (iii) We compute the annual follow-up time requirements and find that prioritization will be needed. (iv) Even a few telescopes could detect sufficiently nearby counterparts, raising the possibility of adding a handful of small-size or medium-size telescopes to the network at diverse geographic locations taking into account the positions of CTA and the LIGO-Virgo-KAGRA network. (v) The continued operation of VERITAS/H.E.S.S./MAGIC would be a useful compliment to CTA's follow-up capabilities by increasing the sky area that can be rapidly covered, especially for directions above and 'below' the United States in which the present network of gravitational-wave detectors is more sensitive.

Kilonovae are cosmic optical flashes produced in the aftermath of the merger of two neutron stars. While the typical radiant flux of a kilonova can be as high as $10^{34}$\,W, they typically occur at cosmological distances, requiring meter-class or larger telescopes for their observation. Here we explore how a kilonova would look like from Earth if it occurred in the Solar System's backyard, 1000 light years from Earth. This is a small distance on cosmological scales, with only one neutron-star merger expected to occur this close every 100,000,000 years. While humans will likely only see kilonovae at much greater distances, showing how such a nearby event would look on camera can help visualize these events, and demonstrate their unique spectral evolution.

The recent multi-messenger discovery of binary neutron star (BNS) merger GW170817 showed that $\gamma$-ray emission in short GRBs is wider than the central energetic narrow cone, and weakly expands out to tens of degrees. Here we explore some of the observational consequences of this structured emission, taking the reconstructed angular emission profile of gamma-ray burst GRB170817A to be typical. We calculate the expected fraction of gravitational-wave observations from BNS mergers that will have an observed GRB counterpart to be $\sim30\%$, implying that joint gravitational-wave and GRB observations will be common. Further, we find that $\sim10\%$ of observed short GRBs occur within 200 Mpc. Finally, we estimate a BNS merger rate of $\sim500$ Gpc$^{-3}$yr$^{-1}$.