-
Low-frequency wideband timing of InPTA pulsars observed with the uGMRT
Authors:
K Nobleson,
Nikita Agarwal,
Raghav Girgaonkar,
Arul Pandian,
Bhal Chandra Joshi,
M A Krishnakumar,
Abhimanyu Susobhanan,
Shantanu Desai,
T Prabu,
Adarsh Bathula,
Timothy T Pennucci,
Sarmistha Banik,
Manjari Bagchi,
Neelam Dhanda Batra,
Arpita Choudhary,
Subhajit Dandapat,
Lankeswar Dey,
Yashwant Gupta,
Shinnosuke Hisano,
Ryo Kato,
Divyansh Kharbanda,
Tomonosuke Kikunaga,
Neel Kolhe,
Yogesh Maan,
Piyush Marmat
, et al. (7 additional authors not shown)
Abstract:
High-precision measurements of the pulsar dispersion measure (DM) are possible using telescopes with low-frequency wideband receivers. We present an initial study of the application of the wideband timing technique, which can simultaneously measure the pulsar times of arrival (ToAs) and DMs, for a set of five pulsars observed with the upgraded Giant Metrewave Radio Telescope (uGMRT) as part of the…
▽ More
High-precision measurements of the pulsar dispersion measure (DM) are possible using telescopes with low-frequency wideband receivers. We present an initial study of the application of the wideband timing technique, which can simultaneously measure the pulsar times of arrival (ToAs) and DMs, for a set of five pulsars observed with the upgraded Giant Metrewave Radio Telescope (uGMRT) as part of the Indian Pulsar Timing Array (InPTA) campaign. We have used the observations with the 300-500 MHz band of the uGMRT for this purpose. We obtain high precision in DM measurements with precisions of the order 10^{-6}cm^{-3}pc. The ToAs obtained have sub-μs precision and the root-mean-square of the post-fit ToA residuals are in the sub-μs range. We find that the uncertainties in the DMs and ToAs obtained with this wideband technique, applied to low-frequency data, are consistent with the results obtained with traditional pulsar timing techniques and comparable to high-frequency results from other PTAs. This work opens up an interesting possibility of using low-frequency wideband observations for precision pulsar timing and gravitational wave detection with similar precision as high-frequency observations used conventionally.
△ Less
Submitted 23 February, 2022; v1 submitted 13 December, 2021;
originally announced December 2021.
-
High Precision Measurements of Interstellar Dispersion Measure with the upgraded GMRT
Authors:
M. A. Krishnakumar,
P. K. Manoharan,
Bhal Chandra Joshi,
Raghav Girgaonkar,
Shantanu Desai,
Manjari Bagchi,
K. Nobleson,
Lankeswar Dey,
Abhimanyu Susobhanan,
Sai Chaitanya Susarla,
Mayuresh P. Surnis,
Yogesh Maan,
A. Gopakumar,
Avishek Basu,
Neelam Dhanda Batra,
Arpita Choudhary,
Kishalay De,
Yashwant Gupta,
Arun Kumar Naidu,
Dhruv Pathak,
Jaikhomba Singha,
T. Prabu
Abstract:
Pulsar radio emission undergoes dispersion due to the presence of free electrons in the interstellar medium (ISM). The dispersive delay in the arrival time of pulsar signal changes over time due to the varying ISM electron column density along the line of sight. Correcting for this delay accurately is crucial for the detection of nanohertz gravitational waves using Pulsar Timing Arrays. In this wo…
▽ More
Pulsar radio emission undergoes dispersion due to the presence of free electrons in the interstellar medium (ISM). The dispersive delay in the arrival time of pulsar signal changes over time due to the varying ISM electron column density along the line of sight. Correcting for this delay accurately is crucial for the detection of nanohertz gravitational waves using Pulsar Timing Arrays. In this work, we present in-band and inter-band DM estimates of four pulsars observed with uGMRT over the timescale of a year using two different template alignment methods. The DMs obtained using both these methods show only subtle differences for PSR 1713+0747 and J1909$-$3744. A considerable offset is seen in the DM of PSR J1939+2134 and J2145$-$0750 between the two methods. This could be due to the presence of scattering in the former and profile evolution in the latter. We find that both methods are useful but could have a systematic offset between the DMs obtained. Irrespective of the template alignment methods followed, the precision on the DMs obtained is about $10^{-3}$ pc cm$^{-3}$ using only BAND3 and $10^{-4}$ pc cm$^{-3}$ after combining data from BAND3 and BAND5 of the uGMRT. In a particular result, we have detected a DM excess of about $5\times10^{-3}$ pc cm$^{-3}$ on 24 February 2019 for PSR J2145$-$0750. This excess appears to be due to the interaction region created by fast solar wind from a coronal hole and a coronal mass ejection (CME) observed from the Sun on that epoch. A detailed analysis of this interesting event is presented.
△ Less
Submitted 1 July, 2021; v1 submitted 13 January, 2021;
originally announced January 2021.
-
Building Halo Merger Trees from the Q Continuum Simulation
Authors:
Esteban Rangel,
Nicholas Frontiere,
Salman Habib,
Katrin Heitmann,
Wei-keng Liao,
Ankit Agrawal,
Alok Choudhary
Abstract:
Cosmological N-body simulations rank among the most computationally intensive efforts today. A key challenge is the analysis of structure, substructure, and the merger history for many billions of compact particle clusters, called halos. Effectively representing the merging history of halos is essential for many galaxy formation models used to generate synthetic sky catalogs, an important applicat…
▽ More
Cosmological N-body simulations rank among the most computationally intensive efforts today. A key challenge is the analysis of structure, substructure, and the merger history for many billions of compact particle clusters, called halos. Effectively representing the merging history of halos is essential for many galaxy formation models used to generate synthetic sky catalogs, an important application of modern cosmological simulations. Generating realistic mock catalogs requires computing the halo formation history from simulations with large volumes and billions of halos over many time steps, taking hundreds of terabytes of analysis data. We present fast parallel algorithms for producing halo merger trees and tracking halo substructure from a single-level, density-based clustering algorithm. Merger trees are created from analyzing the halo-particle membership function in adjacent snapshots, and substructure is identified by tracking the "cores" of merging halos -- sets of particles near the halo center. Core tracking is performed after creating merger trees and uses the relationships found during tree construction to associate substructures with hosts. The algorithms are implemented with MPI and evaluated on a Cray XK7 supercomputer using up to 16,384 processes on data from HACC, a modern cosmological simulation framework. We present results for creating merger trees from 101 analysis snapshots taken from the Q Continuum, a large volume, high mass resolution, cosmological simulation evolving half a trillion particles.
△ Less
Submitted 19 August, 2020;
originally announced August 2020.
-
pinta: The uGMRT Data Processing Pipeline for the Indian Pulsar Timing Array
Authors:
Abhimanyu Susobhanan,
Yogesh Maan,
Bhal Chandra Joshi,
T. Prabu,
Shantanu Desai,
K. Nobleson,
Sai Chaitanya Susarla,
Raghav Girgaonkar,
Lankeswar Dey,
Neelam Dhanda Batra,
Yashwant Gupta,
A. Gopakumar,
Manjari Bagchi,
Avishek Basu,
Suryarao Bethapudi,
Arpita Choudhary,
Kishalay De,
M. A. Krishnakumar,
P. K. Manoharan,
Arun Kumar Naidu,
Dhruv Pathak,
Jaikhomba Singha,
Mayuresh P. Surnis
Abstract:
We introduce pinta, a pipeline for reducing the upgraded Giant Metre-wave Radio Telescope (uGMRT) raw pulsar timing data, developed for the Indian Pulsar Timing Array experiment. We provide a detailed description of the workflow and usage of pinta, as well as its computational performance and RFI mitigation characteristics. We also discuss a novel and independent determination of the relative time…
▽ More
We introduce pinta, a pipeline for reducing the upgraded Giant Metre-wave Radio Telescope (uGMRT) raw pulsar timing data, developed for the Indian Pulsar Timing Array experiment. We provide a detailed description of the workflow and usage of pinta, as well as its computational performance and RFI mitigation characteristics. We also discuss a novel and independent determination of the relative time offsets between the different back-end modes of uGMRT and the interpretation of the uGMRT observation frequency settings, and their agreement with results obtained from engineering tests. Further, we demonstrate the capability of pinta to generate data products which can produce high-precision TOAs using PSR J1909-3744 as an example. These results are crucial for performing precision pulsar timing with the uGMRT.
△ Less
Submitted 16 April, 2021; v1 submitted 6 July, 2020;
originally announced July 2020.
-
GiRaFFE: An Open-Source General Relativistic Force-Free Electrodynamics Code
Authors:
Zachariah B. Etienne,
Mew-Bing Wan,
Maria C. Babiuc,
Sean T. McWilliams,
Ashok Choudhary
Abstract:
We present GiRaFFE, the first open-source general relativistic force-free electrodynamics (GRFFE) code for dynamical, numerical-relativity generated spacetimes. GiRaFFE adopts the strategy pioneered by McKinney and modified by Paschalidis and Shapiro to convert a GR magnetohydrodynamic (GRMHD) code into a GRFFE code. In short, GiRaFFE exists as a modification of IllinoisGRMHD, a user-friendly, ope…
▽ More
We present GiRaFFE, the first open-source general relativistic force-free electrodynamics (GRFFE) code for dynamical, numerical-relativity generated spacetimes. GiRaFFE adopts the strategy pioneered by McKinney and modified by Paschalidis and Shapiro to convert a GR magnetohydrodynamic (GRMHD) code into a GRFFE code. In short, GiRaFFE exists as a modification of IllinoisGRMHD, a user-friendly, open-source, dynamical-spacetime GRMHD code. Both GiRaFFE and IllinoisGRMHD leverage the Einstein Toolkit's highly-scalable infrastructure to make possible large-scale simulations of magnetized plasmas in strong, dynamical spacetimes on adaptive-mesh refinement (AMR) grids. We demonstrate that GiRaFFE passes a large suite of both flat and curved-spacetime code tests passed by a number of other state-of-the-art GRFFE codes, and is thus ready for production-scale simulations of GRFFE phenomena of key interest to relativistic astrophysics.
△ Less
Submitted 1 October, 2017; v1 submitted 30 March, 2017;
originally announced April 2017.
-
Hubble imaging of V1331 Cygni: Proper motion study of its circumstellar structures
Authors:
Arpita Choudhary,
Bringfried Stecklum,
Hendrik Linz
Abstract:
The young star V1331 Cyg received previous attention because it is surrounded by an optical, arc-like reflection nebula. V1331 Cyg is commonly considered to be a candidate for an object that has undergone an FU-Ori (FUOR) the outbreak in the past. This in turn could lead to a time-varying appearance of the dusty arcs that may be revealed by multi-epoch imaging. In particular, a radial colour analy…
▽ More
The young star V1331 Cyg received previous attention because it is surrounded by an optical, arc-like reflection nebula. V1331 Cyg is commonly considered to be a candidate for an object that has undergone an FU-Ori (FUOR) the outbreak in the past. This in turn could lead to a time-varying appearance of the dusty arcs that may be revealed by multi-epoch imaging. In particular, a radial colour analysis of the dust arcs can then be attempted to check whether radial grain size distribution was modified by a previous FUOR wind. Second-epoch imaging of V1331 Cyg was obtained by us in 2009 using the Hubble Space Telescope (HST). By comparing this to archival HST data from 2000, we studied the time evolution of the circumstellar nebulae. After a point spread function subtraction using model point spread functions, we used customised routines to perform a proper motion analysis. The nebula expansion was first derived by deconvolving and correlating the two-epoch radial brightness profiles. Additional data from other facilities TLS, UKIDSS, SPITZER, and HERSCHEL were also incorporated to improve our understanding of the star in terms of environment, viewing angle, bipolar outflow length, and the FUOR phenomenon. The derived radial colour profiles do not indicate a spatial separation of the dust grain sizes. The HERSCHEL 160 micron images show for the time thermal emission from dust probably residing in the outer arc. By viewing V1331 Cyg almost pole-on, the length of the bipolar outflow exceeds previous estimates by far. The outer arc expansion timescale is consistent with the implantation time of the CO torus, which supports the hypothesis of an outburst that occurred a few thousand years ago. The azimuthal colour variation of the outer arc is probably due to changes of the scattering angle, imposed by a tilt or helical geometry of the dust configuration.
△ Less
Submitted 14 April, 2016; v1 submitted 13 April, 2016;
originally announced April 2016.
-
A New Hybrid Technique for Modeling Dense Star Clusters
Authors:
Carl L. Rodriguez,
Bharath Pattabiraman,
Sourav Chatterjee,
Alok Choudhary,
Wei-keng Liao,
Meagan Morscher,
Frederic A. Rasio
Abstract:
The "gravitational million-body problem," to model the dynamical evolution of a self-gravitating, collisional N-body system with ~10^6 particles over many relaxation times, remains a major challenge in computational astrophysics. Unfortunately, current techniques to model such systems suffer from severe limitations. A direct N-body simulation with more than 10^5 particles can require months or eve…
▽ More
The "gravitational million-body problem," to model the dynamical evolution of a self-gravitating, collisional N-body system with ~10^6 particles over many relaxation times, remains a major challenge in computational astrophysics. Unfortunately, current techniques to model such systems suffer from severe limitations. A direct N-body simulation with more than 10^5 particles can require months or even years to complete, while an orbit-sampling Monte Carlo approach cannot adequately model the dynamics in a dense cluster core, particularly in the presence of many black holes. We have developed a new technique combining the precision of a direct N-body integration with the speed of a Monte Carlo approach. Our Rapid And Precisely Integrated Dynamics code, the RAPID code, statistically models interactions between neighboring stars and stellar binaries while integrating directly the orbits of stars or black holes in the cluster core. This allows us to accurately simulate the dynamics of the black holes in a realistic globular cluster environment without the burdensome N^2 scaling of a full N-body integration. We compare RAPID models of idealized globular clusters to identical models from the direct N-body and Monte Carlo methods. Our tests show that RAPID can reproduce the half-mass radii, core radii, black hole ejection rates, and binary properties of the direct N-body models far more accurately than a standard Monte Carlo integration while remaining significantly faster than a full N-body integration. With this technique, it will be possible to create more realistic models of Milky Way globular clusters with sufficient rapidity to explore the full parameter space of dense stellar clusters.
△ Less
Submitted 13 November, 2018; v1 submitted 2 November, 2015;
originally announced November 2015.
-
A Parallel Monte Carlo Code for Simulating Collisional N-body Systems
Authors:
Bharath Pattabiraman,
Stefan Umbreit,
Wei-Keng Liao,
Alok Choudhary,
Vassiliki Kalogera,
Gokhan Memik,
Frederic A. Rasio
Abstract:
We present a new parallel code for computing the dynamical evolution of collisional N-body systems with up to N~10^7 particles. Our code is based on the the Henon Monte Carlo method for solving the Fokker-Planck equation, and makes assumptions of spherical symmetry and dynamical equilibrium. The principal algorithmic developments involve optimizing data structures, and the introduction of a parall…
▽ More
We present a new parallel code for computing the dynamical evolution of collisional N-body systems with up to N~10^7 particles. Our code is based on the the Henon Monte Carlo method for solving the Fokker-Planck equation, and makes assumptions of spherical symmetry and dynamical equilibrium. The principal algorithmic developments involve optimizing data structures, and the introduction of a parallel random number generation scheme, as well as a parallel sorting algorithm, required to find nearest neighbors for interactions and to compute the gravitational potential. The new algorithms we introduce along with our choice of decomposition scheme minimize communication costs and ensure optimal distribution of data and workload among the processing units. The implementation uses the Message Passing Interface (MPI) library for communication, which makes it portable to many different supercomputing architectures. We validate the code by calculating the evolution of clusters with initial Plummer distribution functions up to core collapse with the number of stars, N, spanning three orders of magnitude, from 10^5 to 10^7. We find that our results are in good agreement with self-similar core-collapse solutions, and the core collapse times generally agree with expectations from the literature. Also, we observe good total energy conservation, within less than 0.04% throughout all simulations. We analyze the performance of the code, and demonstrate near-linear scaling of the runtime with the number of processors up to 64 processors for N=10^5, 128 for N=10^6 and 256 for N=10^7. The runtime reaches a saturation with the addition of more processors beyond these limits which is a characteristic of the parallel sorting algorithm. The resulting maximum speedups we achieve are approximately 60x, 100x, and 220x, respectively.
△ Less
Submitted 15 November, 2012; v1 submitted 25 June, 2012;
originally announced June 2012.
-
The role of axisymmetric flow configuration in the estimation of the analogue surface gravity and related Hawking like temperature
Authors:
Neven Bilic,
Arpita Choudhary,
Tapas K Das,
Sankhasubhra Nag
Abstract:
For axially symmetric flow of dissipationless inhomogeneous fluid onto a non rotating astrophysical black hole under the influence of a generalized pseudo-Schwarzschild gravitational potential, we investigate the influence of the background flow configuration on determining the salient features of the corresponding acoustic geometry. The acoustic horizon for the aforementioned flow structure has b…
▽ More
For axially symmetric flow of dissipationless inhomogeneous fluid onto a non rotating astrophysical black hole under the influence of a generalized pseudo-Schwarzschild gravitational potential, we investigate the influence of the background flow configuration on determining the salient features of the corresponding acoustic geometry. The acoustic horizon for the aforementioned flow structure has been located and the corresponding acoustic surface gravity $κ$ as well as the associated analogue Hawking temperature $T_{\rm AH}$ has been calculated {\it analytically}. The dependence of $κ$ on the flow geometry as well as on the nature of the back ground black hole space time (manifested through the nature of the pseudo-Schwarzschild potential used) has been discussed. Dependence of the value of $κ$ on various initial boundary conditions governing the dynamic and the thermodynamic properties of the background fluid flow has also been studied.
△ Less
Submitted 24 May, 2012;
originally announced May 2012.