-
Primordial magnetic relics and their signatures
Authors:
Arka Banerjee,
Lalit Singh Bhandari,
Ashwat Jain,
Arun M. Thalapillil
Abstract:
Primordial black holes bearing magnetic charges may bypass the constraints imposed by Hawking radiation, thereby enabling reasonable present-day populations, even for masses below $10^{15}\,\text{g}$ -- a range previously considered improbable. They could, therefore, conceivably contribute to a component of dark matter. We investigate novel Faraday rotation signatures exhibited by primordial magne…
▽ More
Primordial black holes bearing magnetic charges may bypass the constraints imposed by Hawking radiation, thereby enabling reasonable present-day populations, even for masses below $10^{15}\,\text{g}$ -- a range previously considered improbable. They could, therefore, conceivably contribute to a component of dark matter. We investigate novel Faraday rotation signatures exhibited by primordial magnetic black holes while also establishing new Parker-type bounds on their populations. For the latter, we bound the dark matter fraction from intergalactic magnetic fields in cosmic voids $\left(f_{\text{DM}} \lesssim 10^{-8}\right)$ and cosmic web filaments $\left(f_{\text{ DM}} \lesssim 10^{-4}\right)$, notably eclipsing previous bounds. Exploring Faraday rotation effects, we discern a pronounced rotation of the polarization angle and the rotation measure values for extremal primordial magnetic black holes with masses $M^{\text{ ex.}}_{\text{ BH}}\gtrsim 10^{-6}~ \text{M}_\odot$. This makes them potentially detectable in current observations. A comparative investigation finds that the effects are notably greater than for a neutron star, like a Magnetar, with a similar magnetic field at the surface. Moreover, the polarization angle maps for primordial magnetic black holes exhibit unique features, notably absent in other astrophysical magnetic configurations. In this context, we also introduce a simple integral measure, offering a quantitative measure for their discrimination in many scenarios. These traits potentially suggest a robust avenue for their observational detection and differentiation.
△ Less
Submitted 12 June, 2024;
originally announced June 2024.
-
Small-scale cosmic ray anisotropy observed by the GRAPES-3 experiment at TeV energies
Authors:
M. Chakraborty,
S. Ahmad,
A. Chandra,
S. R. Dugad,
U. D. Goswami,
S. K. Gupta,
B. Hariharan,
Y. Hayashi,
P. Jagadeesan,
A. Jain,
P. Jain,
S. Kawakami,
T. Koi,
H. Kojima,
S. Mahapatra,
P. K. Mohanty,
R. Moharana,
Y. Muraki,
P. K. Nayak,
T. Nonaka,
T. Nakamura,
A. Oshima,
B. P. Pant,
D. Pattanaik,
S. Paul
, et al. (13 additional authors not shown)
Abstract:
GRAPES-3 is a mid-altitude (2200 m) and near equatorial ($11.4^{\circ}$ North) air shower array, overlapping in its field of view for cosmic ray observations with experiments that are located in Northern and Southern hemispheres. We analyze a sample of $3.7\times10^9$ cosmic ray events collected by the GRAPES-3 experiment between 1 January 2013 and 31 December 2016 with a median energy of…
▽ More
GRAPES-3 is a mid-altitude (2200 m) and near equatorial ($11.4^{\circ}$ North) air shower array, overlapping in its field of view for cosmic ray observations with experiments that are located in Northern and Southern hemispheres. We analyze a sample of $3.7\times10^9$ cosmic ray events collected by the GRAPES-3 experiment between 1 January 2013 and 31 December 2016 with a median energy of $\sim16$ TeV for study of small-scale ($<60^{\circ}$) angular scale anisotropies. We observed two structures labeled as A and B, deviate from the expected isotropic distribution of cosmic rays in a statistically significant manner. Structure `A' spans $50^{\circ}$ to $80^{\circ}$ in the right ascension and $-15^{\circ}$ to $30^{\circ}$ in the declination coordinate. The relative excess observed in the structure A is at the level of $(6.5\pm1.3)\times10^{-4}$ with a statistical significance of 6.8 standard deviations. Structure `B' is observed in the right ascension range of $110^{\circ}$ to $140^{\circ}$. The relative excess observed in this region is at the level of $(4.9\pm1.4)\times10^{-4}$ with a statistical significance of 4.7 standard deviations. These structures are consistent with those reported by Milagro, ARGO-YBJ, and HAWC. These observations could provide a better understanding of the cosmic ray sources, propagation and the magnetic structures in our Galaxy.
△ Less
Submitted 23 October, 2023;
originally announced October 2023.
-
AOTF based spectro-polarimeter for observing Earth as an Exoplanet
Authors:
Bhavesh Jaiswal,
Swapnil Singh,
Anand Jain,
K Sankarasubramanian,
Anuj Nandi
Abstract:
Earth is the only known habitable planet and it serves as a testbed to benchmark the observations of temperate and more Earth-like exoplanets. It is required to observe the disc-integrated signatures of Earth for a large range of phase angles, resembling the observations of an exoplanet. In this work, an AOTF (Acousto-Optic Tunable Filter) based experiment is designed to observe the spectro-polari…
▽ More
Earth is the only known habitable planet and it serves as a testbed to benchmark the observations of temperate and more Earth-like exoplanets. It is required to observe the disc-integrated signatures of Earth for a large range of phase angles, resembling the observations of an exoplanet. In this work, an AOTF (Acousto-Optic Tunable Filter) based experiment is designed to observe the spectro-polarimetric signatures of Earth. The results of spectroscopic and polarimetric laboratory calibration are presented here along with a brief overview of a possible instrument configuration. Based on the results of the spectro-polarimetric calibration, simulations are carried out to optimize the instrument design for the expected signal levels for various observing conditions. The usefulness of an AOTF based spectro-polarimeter is established from this study and it is found that, in the present configuration, the instrument can achieve a polarimetric accuracy of $<0.3$\% for linear polarization for an integration time of 100 ms or larger. The design configuration of the instrument and the planning of conducting such observations from Lunar orbit are discussed.
△ Less
Submitted 21 February, 2023;
originally announced February 2023.
-
The Astropy Project: Sustaining and Growing a Community-oriented Open-source Project and the Latest Major Release (v5.0) of the Core Package
Authors:
The Astropy Collaboration,
Adrian M. Price-Whelan,
Pey Lian Lim,
Nicholas Earl,
Nathaniel Starkman,
Larry Bradley,
David L. Shupe,
Aarya A. Patil,
Lia Corrales,
C. E. Brasseur,
Maximilian Nöthe,
Axel Donath,
Erik Tollerud,
Brett M. Morris,
Adam Ginsburg,
Eero Vaher,
Benjamin A. Weaver,
James Tocknell,
William Jamieson,
Marten H. van Kerkwijk,
Thomas P. Robitaille,
Bruce Merry,
Matteo Bachetti,
H. Moritz Günther,
Thomas L. Aldcroft
, et al. (111 additional authors not shown)
Abstract:
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 $\texttt{astropy}$, which serves as the foundation for more specialized projects and packages. In this article, we summarize key features in the core package as…
▽ More
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 $\texttt{astropy}$, which serves as the foundation for more specialized projects and packages. In this article, we summarize key features in the core package as of the recent major release, version 5.0, and provide major updates for the Project. We then discuss supporting a broader ecosystem of interoperable packages, including connections with several astronomical observatories and missions. We also revisit the future outlook of the Astropy Project and the current status of Learn Astropy. We conclude by raising and discussing the current and future challenges facing the Project.
△ Less
Submitted 28 June, 2022;
originally announced June 2022.
-
New Strong Gravitational Lenses from the DESI Legacy Imaging Surveys Data Release 9
Authors:
C. Storfer,
X. Huang,
A. Gu,
W. Sheu,
S. Banka,
A. Dey,
A. Jain,
J. Kwon,
D. Lang,
V. Lee,
A. Meisner,
J. Moustakas,
A. D. Myers,
S. Tabares-Tarquinio,
E. F. Schlafly,
D. J. Schlegel
Abstract:
We have conducted a search for strong gravitational lensing systems in the Dark Energy Spectroscopic Instrument (DESI) Legacy Imaging Surveys Data Release 9. This is the third paper in a series (following Huang et al. 2020; Huang et al. 2021, Paper I & II, respectively). These surveys together cover $\sim$ 19,000 deg$^2$ visible from the northern hemisphere, reaching a z-band AB magnitude of…
▽ More
We have conducted a search for strong gravitational lensing systems in the Dark Energy Spectroscopic Instrument (DESI) Legacy Imaging Surveys Data Release 9. This is the third paper in a series (following Huang et al. 2020; Huang et al. 2021, Paper I & II, respectively). These surveys together cover $\sim$ 19,000 deg$^2$ visible from the northern hemisphere, reaching a z-band AB magnitude of $\sim$ 22.5. We use a deep residual neural network, trained on a compilation of known lensing systems and candidates as well as non-lenses in the same footprint. After applying our trained neural networks to the survey data, we visually inspect and rank images with probabilities above a threshold. We have found 1895 lens candidates. Out of these, 1512 are identified for the first time. Combining the discoveries from this work, Paper I (335) and II (1210), the total number of strong lens candidates from the Legacy Surveys that we have discovered is 3057.
△ Less
Submitted 26 April, 2023; v1 submitted 6 June, 2022;
originally announced June 2022.
-
EinsteinPy: A Community Python Package for General Relativity
Authors:
Shreyas Bapat,
Ritwik Saha,
Bhavya Bhatt,
Shilpi Jain,
Akshita Jain,
Sofía Ortín Vela,
Priyanshu Khandelwal,
Jyotirmaya Shivottam,
Jialin Ma,
Gim Seng Ng,
Pratyush Kerhalkar,
Hrishikesh Sudam Sarode,
Rishi Sharma,
Manvi Gupta,
Divya Gupta,
Tushar Tyagi,
Tanmay Rustagi,
Varun Singh,
Saurabh Bansal,
Naman Tayal,
Abhijeet Manhas,
Raphael Reyna,
Gaurav Kumar,
Govind Dixit,
Ratin Kumar
, et al. (24 additional authors not shown)
Abstract:
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 increas…
▽ More
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.
△ Less
Submitted 22 May, 2020;
originally announced May 2020.
-
Energy sensitivity of the GRAPES-3 EAS array for primary cosmic ray protons
Authors:
B. Hariharan,
S. Ahmad,
M. Chakraborty,
A. Chandra,
S. R. Dugad,
S. K. Gupta,
Y. Hayashi,
H. Kojima,
S. S. R. Inbanathan,
P. Jagadeesan,
A. Jain,
P. Jain,
V. B. Jhansi,
S. Kawakami,
P. K. Mohanty,
S. D. Morris,
P. K. Nayak,
A. Oshima,
D. Pattanaik,
P. S. Rakshe,
K. Ramesh,
B. S. Rao,
L. V. Reddy,
S. Shibata,
F. Varsi
, et al. (1 additional authors not shown)
Abstract:
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 showe…
▽ More
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.
△ Less
Submitted 3 May, 2020;
originally announced May 2020.
-
The angular resolution of GRAPES-3 EAS array after correction for the shower front curvature
Authors:
V. B. Jhansi,
S. Ahmad,
M. Chakraborty,
S. R. Dugad,
S. K. Gupta,
B. Hariharan,
Y. Hayashi,
P. Jagadeesan,
A. Jain,
P. Jain,
S. Kawakami,
H. Kojima,
S. Mahapatra,
P. K. Mohanty,
S. D. Morris,
P. K. Nayak,
A. Oshima,
D. Pattanaik,
P. S. Rakshe,
K. Ramesh,
B. S. Rao,
L. V. Reddy,
S. Shibata,
F. Varsi,
M. Zuberi
Abstract:
The angular resolution of an extensive air shower (EAS) array plays a critical role in determining its sensitivity for the detection of point $γ$-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 $γ$-rays in the TeV-PeV energy range. It comprises of a dense array of 400 plastic scin…
▽ More
The angular resolution of an extensive air shower (EAS) array plays a critical role in determining its sensitivity for the detection of point $γ$-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 $γ$-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.
△ Less
Submitted 12 November, 2019;
originally announced November 2019.
-
Modeling of rigidity dependent CORSIKA simulations for GRAPES-3
Authors:
B. Hariharan,
S. R. Dugad,
S. K. Gupta,
Y. Hayashi,
S. S. R. Inbanathan,
P. Jagadeesan,
A. Jain,
S. Kawakami,
P. K. Mohanty,
B. S. Rao
Abstract:
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 o…
▽ More
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.
△ Less
Submitted 16 August, 2019;
originally announced August 2019.
-
One-form superfluids and magnetohydrodynamics
Authors:
Jay Armas,
Akash Jain
Abstract:
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…
▽ More
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.
△ Less
Submitted 8 January, 2020; v1 submitted 12 November, 2018;
originally announced November 2018.
-
Magnetohydrodynamics as superfluidity
Authors:
Jay Armas,
Akash Jain
Abstract:
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 potenti…
▽ More
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.
△ Less
Submitted 22 March, 2019; v1 submitted 6 August, 2018;
originally announced August 2018.
-
Was the cosmic ray burst detected by the GRAPES-3 on 22 June 2015 caused by transient weakening of geomagnetic field or by an interplanetary anisotropy?
Authors:
P. K. Mohanty,
K. P. Arunbabu,
T. Aziz,
S. R. Dugad,
S. K. Gupta,
B. Hariharan,
P. Jagadeesan,
A. Jain,
S. D. Morris,
P. K. Nayak,
P. S. Rakshe,
K. Ramesh,
B. S. Rao,
M. Zuberi,
Y. Hayashi,
S. Kawakami,
P. Subramanian,
S. Raha,
S. Ahmad,
A. Oshima,
S. Shibata,
H. Kojima
Abstract:
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$σ$ 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 (stor…
▽ More
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$σ$ 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.
△ Less
Submitted 28 March, 2018;
originally announced March 2018.
-
Variational Inference as an alternative to MCMC for parameter estimation and model selection
Authors:
Geetakrishnasai Gunapati,
Anirudh Jain,
P. K. Srijith,
Shantanu Desai
Abstract:
Most applications of Bayesian Inference for parameter estimation and model selection in astrophysics involve the use of Monte Carlo techniques such as Markov Chain Monte Carlo (MCMC) and nested sampling. However, these techniques are time consuming and their convergence to the posterior could be difficult to determine. In this work, we advocate Variational inference as an alternative to solve the…
▽ More
Most applications of Bayesian Inference for parameter estimation and model selection in astrophysics involve the use of Monte Carlo techniques such as Markov Chain Monte Carlo (MCMC) and nested sampling. However, these techniques are time consuming and their convergence to the posterior could be difficult to determine. In this work, we advocate Variational inference as an alternative to solve the above problems, and demonstrate its usefulness for parameter estimation and model selection in Astrophysics. 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 characterize the difference. It takes advantage of fast optimization techniques, which make it ideal to deal with large datasets and makes it trivial to parallelize on a multicore platform. We also derive a new approximate evidence estimation based on variational posterior, and importance sampling technique called posterior weighted importance sampling for the calculation of evidence (PWISE), which is useful to perform Bayesian model selection. As a proof of principle, we apply variational inference to five different problems in astrophysics, where Monte Carlo techniques were previously used. These include assessment of significance of annual modulation in the COSINE-100 dark matter experiment, 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 and nested sampling techniques for most of these problems while providing competitive results. All our analysis codes have been made publicly available.
△ Less
Submitted 23 December, 2021; v1 submitted 17 March, 2018;
originally announced March 2018.
-
The Astropy Project: Building an inclusive, open-science project and status of the v2.0 core package
Authors:
The Astropy Collaboration,
A. M. Price-Whelan,
B. M. Sipőcz,
H. M. Günther,
P. L. Lim,
S. M. Crawford,
S. Conseil,
D. L. Shupe,
M. W. Craig,
N. Dencheva,
A. Ginsburg,
J. T. VanderPlas,
L. D. Bradley,
D. Pérez-Suárez,
M. de Val-Borro,
T. L. Aldcroft,
K. L. Cruz,
T. P. Robitaille,
E. J. Tollerud,
C. Ardelean,
T. Babej,
M. Bachetti,
A. V. Bakanov,
S. P. Bamford,
G. Barentsen
, et al. (112 additional authors not shown)
Abstract:
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 p…
▽ More
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.
△ Less
Submitted 16 January, 2018; v1 submitted 8 January, 2018;
originally announced January 2018.
-
Binary corrected X-ray light curve of Cygnus X-3 : implications for the timing properties of the compact binary system
Authors:
M. Choudhury,
A. R. Rao,
S. V. Vadawale,
A. K. Jain,
N. S. Singh
Abstract:
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 f…
▽ More
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.
△ Less
Submitted 22 March, 2004;
originally announced March 2004.
-
Correlated radio:X-ray emission in the hard states of Galactic microquasars
Authors:
M. Choudhury,
A. R. Rao,
S. V. Vadawale,
A. K. Jain
Abstract:
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 o…
▽ More
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.
△ Less
Submitted 23 June, 2003; v1 submitted 22 April, 2003;
originally announced April 2003.
-
Disk-Jet Connection in Cygnus X-3
Authors:
M. Choudhury,
A. R. Rao,
S. V. Vadawale,
C. H. Ishwara-Chandra,
A. K. Jain
Abstract:
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 an…
▽ More
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.
△ Less
Submitted 17 January, 2002; v1 submitted 21 November, 2001;
originally announced November 2001.
-
Neutrino-Induced Giant Air Showers in Large Extra Dimension Models
Authors:
Ambar Jain,
Pankaj Jain,
Douglas W. McKay,
John P. Ralston
Abstract:
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 sho…
▽ More
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.
△ Less
Submitted 7 February, 2001; v1 submitted 24 November, 2000;
originally announced November 2000.
