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.

Low energy ground-based cosmic ray air shower experiments generally have energy threshold in the range of a few tens to a few hundreds of TeV. The shower observables are measured indirectly with an array of detectors. The atmospheric absorption of low energy secondaries limits their detection frequencies at the Earth's surface. However, due to selection effects, a tiny fraction of low energy showers, which are produced in the lower atmosphere can reach the observational level. But, due to less information of shower observables, the reconstruction of these showers are arduous. Hence, it is believed that direct measurements by experiments aboard on satellites and balloon flights are more reliable at low energies. Despite having very small efficiency ($\sim$0.1%) at low energies, the large acceptance ($\sim$5 m$^2$sr) of GRAPES-3 experiment allows observing primary cosmic rays down below to $\sim$1 TeV and opens up the possibility to measure primary energy spectrum spanning from a few TeV to beyond cosmic ray knee (up to 10$^{16}$ eV), covering five orders of magnitude. The GRAPES-3 energy threshold for primary protons through Monte Carlo simulations are calculated, which gives reasonably good agreement with data. Furthermore, the total efficiencies and acceptance are also calculated for protons primaries. The ability of GRAPES-3 experiment to cover such a broader energy range may provide a unique handle to bridge the energy spectrum between direct measurements at low energies and indirect measurements at ultra-high energies.

We use the framework of generalised global symmetries to study various hydrodynamic regimes of hot electromagnetism. We formulate the hydrodynamic theories with an unbroken or a spontaneously broken U(1) one-form symmetry. The latter of these describes a one-form superfluid, which is characterised by a vector Goldstone mode and a two-form superfluid velocity. Two special limits of this theory have been studied in detail: the string fluid limit where the U(1) one-form symmetry is partly restored, and the electric limit in which the symmetry is completely broken. The transport properties of these theories are investigated in depth by studying the constraints arising from the second law of thermodynamics and Onsager's relations at first order in derivatives. We also construct a hydrostatic effective action for the Goldstone modes in these theories and use it to characterise the space of all equilibrium configurations. To make explicit contact with hot electromagnetism, the traditional treatment of magnetohydrodynamics, where the electromagnetic photon is incorporated as dynamical degrees of freedom, is extended to include parity-violating contributions. We argue that the chemical potential and electric fields are not independently dynamical in magnetohydrodynamics, and illustrate how to eliminate these within the hydrodynamic derivative expansion using Maxwell's equations. Additionally, a new hydrodynamic theory of non-conducting, but polarised, plasmas is formulated, focusing primarily on the magnetically dominated sector. Finally, it is shown that the different limits of one-form superfluids formulated in terms of generalised global symmetries are exactly equivalent to magnetohydrodynamics and the hydrodynamics of non-conducting plasmas at the non-linear level.

The angular resolution of an extensive air shower (EAS) array plays a critical role in determining its sensitivity for the detection of point $\gamma$-ray sources in the multi-TeV energy range. The GRAPES-3 an EAS array located at Ooty in India (11.4$^{\circ}$N, 76.7$^{\circ}$E, 2200 m altitude) is designed to study $\gamma$-rays in the TeV-PeV energy range. It comprises of a dense array of 400 plastic scintillators deployed over an area of 25000 m$^2$ and a large area (560 m$^2$) muon telescope. A new statistical method allowed real time determination of the propagation delay of each detector in the GRAPES-3 array. The shape of shower front is known to be curved and here the details of a new method developed for accurate measurement of the shower front curvature is presented. These two developments have led to a sizable improvement in the angular resolution of GRAPES-3 array. It is shown that the curvature depends on the size and age of an EAS. By employing two different techniques, namely, the odd-even and the left-right methods, independent estimates of the angular resolution are obtained. The odd-even method estimates the best achievable resolution of the array. For obtaining the angular resolution, the left-right method is used after implementing the size and age dependent curvature corrections. A comparison of the angular resolution as a function of EAS energy by these two methods shows them be virtually indistinguishable. The angular resolution of GRAPES-3 array is 47$^{\prime}$ for energies E$>$5 TeV and improves to 17$^{\prime}$ at E$>$100 TeV and finally approaching 10$^{\prime}$ at E$>$500 TeV.

The GRAPES-3 muon telescope located in Ooty, India records 4x10^9 muons daily. These muons are produced by interaction of primary cosmic rays (PCRs) in the atmosphere. The high statistics of muons enables GRAPES-3 to make precise measurement of various sun-induced phenomenon including coronal mass ejections (CME), Forbush decreases, geomagnetic storms (GMS) and atmosphere acceleration during the overhead passage of thunderclouds. However, the understanding and interpretation of observed data requires Monte Carlo (MC) simulation of PCRs and subsequent development of showers in the atmosphere. CORSIKA is a standard MC simulation code widely used for this purpose. However, these simulations are time consuming as large number of interactions and decays need to be taken into account at various stages of shower development from top of the atmosphere down to ground level. Therefore, computing resources become an important consideration particularly when billion of PCRs need to be simulated to match the high statistical accuracy of the data. During the GRAPES-3 simulations, it was observed that over 60% of simulated events don't really reach the Earth's atmosphere. The geomagnetic field (GMF) creates a threshold to PCRs called cutoff rigidity Rc, a direction dependent parameter below which PCRs can't reach the Earth's atmosphere. However, in CORSIKA there is no provision to set a direction dependent threshold. We have devised an efficient method that has taken into account of this Rc dependence. A reduction by a factor ~3 in simulation time and ~2 in output data size was achieved for GRAPES-3 simulations. This has been incorporated in CORSIKA version v75600 onwards. Detailed implementation of this along the potential benefits are discussed in this work.

We show that relativistic magnetohydrodynamics (MHD) can be recast as a novel theory of superfluidity. This new theory formulates MHD just in terms of conservation equations, including dissipative effects, by introducing appropriate variables such as a magnetic scalar potential, and providing necessary and sufficient conditions to obtain equilibrium configurations. We show that this scalar potential can be interpreted as a Goldstone mode originating from the spontaneous breaking of a one-form symmetry, and present the most generic constitutive relations at one derivative order for a parity-preserving plasma in this new superfluid formulation.

Most applications of Bayesian Inference for parameter estimation and model selection in astrophysics involve the use of Markov Chain Monte Carlo (MCMC) techniques. In this work, we introduce Variational Inference as an alternative to solve these problems, and compare how the results hold up to MCMC methods. Variational Inference converts the inference problem into an optimization problem by approximating the posterior from a known family of distributions and using Kullback-Leibler divergence to measure closeness. Variational Inference takes advantage of fast optimization techniques, which make it ideal to deal with large datasets and also makes it trivial to parallelize. As a proof of principle, we apply Variational Inference to four different problems in astrophysics, where MCMC techniques were previously used. These include measuring exoplanet orbital parameters from radial velocity data, tests of periodicities in measurements of Newton's constant $G$, assessing the significance of a turnover in the spectral lag data of GRB 160625B , and estimating the mass of a galaxy cluster using weak gravitational lensing. We find that Variational Inference is much faster than MCMC for these problems.

The GRAPES-3 muon telescope in Ooty, India had claimed detection of a 2 hour (h) high-energy ($\sim$20 GeV) burst of galactic cosmic-rays (GCRs) through a $>$50$\sigma$ surge in GeV muons, was caused by reconnection of the interplanetary magnetic field (IMF) in the magnetosphere that led to transient weakening of Earth's magnetic shield. This burst had occurred during a G4-class geomagnetic storm (storm) with a delay of $\frac{1}{2}$h relative to the coronal mass ejection (CME) of 22 June 2015 (Mohanty et al., 2016). However, recently a group interpreted the occurrence of the same burst in a subset of 31 neutron monitors (NMs) to have been the result of an anisotropy in interplanetary space (Evenson et al., 2017) in contrast to the claim in (Mohanty et al., 2016). A new analysis of the GRAPES-3 data with a fine 10.6$^{\circ}$ angular segmentation shows the speculation of interplanetary anisotropy to be incorrect, and offers a possible explanation of the NM observations. The observed 28 minutes (min) delay of the burst relative to the CME can be explained by the movement of the reconnection front from the bow shock to the surface of Earth at an average speed of 35 km/s, much lower than the CME speed of 700 km/s. This measurement may provide a more accurate estimate of the start of the storm.

The Astropy project supports and fosters the development of open-source and openly-developed Python packages that provide commonly-needed functionality to the astronomical community. A key element of the Astropy project is the core package Astropy, which serves as the foundation for more specialized projects and packages. In this article, we provide an overview of the organization of the Astropy project and summarize key features in the core package as of the recent major release, version 2.0. We then describe the project infrastructure designed to facilitate and support development for a broader ecosystem of inter-operable packages. We conclude with a future outlook of planned new features and directions for the broader Astropy project.

We study the binary corrected X-ray light curve of Cygnus X-3 using the RXTE-ASM data and RXTE-PCA data. We find that the X-ray template derived from the EXOSAT observations adequately explains the binary variations. Using the recently determined quadratic binary ephemeris we correct the RXTE-PCA lightcurve and obtain the power spectrum which shows distinct features of shifting towards the low frequency regime vis-a-vis the Galactic X-ray binaries. The power spectrum doesn't vary from the hard to the soft state. We also obtain the binary corrected RXTE-ASM monitoring lightcurve and examine the previously obtained correlation between soft X-ray (2-12 keV from RXTE-ASM), hard X-ray (20-100 keV from CGRO-BATSE) and radio (2.2 GHz from GBI) after the binary correction. We find that the correlation time scale is less than a day.

We present results of our study of correlated radio and X-ray emission in two black hole candidates and Galactic microquasars GRS 1915+105 and Cygnus X-1 in their steady long term hard states, along with Cygnus X-3 (using data obtained from RXTE-ASM, CGRO-BATSE and GBI). We detect a pivotal behavior in the X-ray spectrum of GRS 1915+105, correlated to the radio emission. Similar to the results obtained for Cygnus X-3, the flux of X-rays softer than the pivoting point correlates with the radio emission, while the corresponding harder X-ray flux anti-correlates with both the radio and the softer X-ray emission, in this state. We examine all the previously reported correlations of X-ray properties with the radio emission in Galactic microquasars and argue that these are consistent with a general picture where a spectral pivoting is a common feature in these sources with the shape of the spectrum determining the flux of radio emission, during the hard states. We also detect a general monotonic increase in the radio emission of these sources with the soft X-ray emission spanning about 5 orders of magnitude. We qualitatively explain these findings with a Two Component Advective Flow model where the location of a boundary layer between the thin disk and the Comptonizing region determines the spectral shape and also the amount of outflow.

We present the results of a detailed correlation study between the soft X-ray, hard X-ray, and radio emission (obtained from RXTE ASM, BATSE, and GBI observations, respectively) of the bright radio emitting Galactic X-ray binary Cygnus X-3. We detect a very strong positive correlation between the soft X-ray and radio emission during the low-hard and minor flaring periods of the source, and an anti-correlation between the soft and hard X-ray emissions. We present statistical arguments to suggest that the anti-correlation between the radio and hard X-ray emission, reported earlier, is primarily due to their correlation and anti-correlation, respectively, with the soft X-ray emission. We make a wide band X-ray spectral study using the pointed RXTE observations and detect a pivotal behaviour in the X-ray spectrum. We argue that this X-ray spectral pivoting is responsible for the anti-correlation between the soft and hard X-ray emissions. The strong correlation between the soft X-ray and radio emission suggests a close link between the accreting mechanism, plasma cloud surrounding the compact object and the radio emission.

In models based on large extra dimensions where massive spin 2 exchange can dominate at high energies, the neutrino-proton cross section can rise to typical hadronic values at energies above 10^20 eV. The neutrino then becomes a candidate for the primary that initiates the highest energy cosmic ray showers. We investigate characteristics of neutrino-induced showers compared to proton-induced showers. The comparison includes study of starting depth, profile with depth, lateral particle distribution at ground and muon lateral distribution at ground level. We find that for cross sections above 20 mb there are regions of parameter space where the two types of showers are essentially indistinguishable. We conclude that the neutrino candidate hypothesis cannot be ruled out on the basis of shower characteristics.