Source author record

Soumya D. Mohanty

Soumya D. Mohanty appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.IM gr-qc physics.data-an astro-ph astro-ph.GA astro-ph.HE Computation Machine Learning Methodology Neural and Evolutionary Computing physics.ed-ph

Catalog footprint

What is connected

10works

11topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Iterative time-domain method for resolving multiple gravitational wave sources in Pulsar Timing Array data

The sensitivity of ongoing searches for gravitational wave (GW) sources in the ultra-low frequency regime ($10^{-9}$ Hz to $10^{-7}$ Hz) using Pulsar Timing Arrays (PTAs) will continue to increase in the future as more well-timed pulsars are added to the arrays. It is expected that next-generation radio telescopes, namely, the Five-hundred-meter Aperture Spherical radio Telescope (FAST) and the Square Kilometer Array (SKA), will grow the number of well-timed pulsars to $O(10^3)$. The higher sensitivity will result in greater distance reach for GW sources, uncovering multiple resolvable GW sources in addition to an unresolved population. Data analysis techniques are, therefore, required that can search for and resolve multiple signals present simultaneously in PTA data. The multisource resolution problem in PTA data analysis poses a unique set of challenges such as non-uniformly sampled data, a large number of so-called pulsar phase parameters that arise from the inaccurately measured distances to the pulsars, and poor separation of signals in the Fourier domain due to a small number of cycles in the observed waveforms.

preprint2021arXiv

Extending the frequency reach of pulsar timing array based gravitational wave search without high cadence observations

Gravitational wave (GW) searches using pulsar timing arrays (PTAs) are assumed to be limited by the typical average observational cadence of $1/(2~{\rm weeks})$ for a single pulsar to GW frequencies $\lesssim 4\times 10^{-7}$ Hz. We show that this assumption is incorrect and that a PTA can detect signals with much higher frequencies, which are preserved in the data due to aliasing, by exploiting asynchronous observations from multiple pulsars. This allows an observation strategy that is scalable to future large-scale PTAs containing $O(10^3)$ pulsars, enabled by the Five-hundred meter Aperture Spherical Telescope and the Square Kilometer Array, without requiring a higher per-pulsar observation cadence. We show that higher frequency GW observations, reaching up to $4\times 10^{-4}$ Hz with an SKA-era PTA, have significant astrophysical implications, such as (i) a three orders of magnitude better constraint than current high-cadence observations on GW strain in the $[10,400]$ $μ{\rm Hz}$ band, and (ii) sensitive tests of the no-hair theorem in the mass range of supermassive black hole binaries using their inspiral, merger, and ringdown signals.

preprint2020arXiv

Adaptive spline fitting with particle swarm optimization

In fitting data with a spline, finding the optimal placement of knots can significantly improve the quality of the fit. However, the challenging high-dimensional and non-convex optimization problem associated with completely free knot placement has been a major roadblock in using this approach. We present a method that uses particle swarm optimization (PSO) combined with model selection to address this challenge. The problem of overfitting due to knot clustering that accompanies free knot placement is mitigated in this method by explicit regularization, resulting in a significantly improved performance on highly noisy data. The principal design choices available in the method are delineated and a statistically rigorous study of their effect on performance is carried out using simulated data and a wide variety of benchmark functions. Our results demonstrate that PSO-based free knot placement leads to a viable and flexible adaptive spline fitting approach that allows the fitting of both smooth and non-smooth functions.

preprint2020arXiv

Towards a real-time fully-coherent all-sky search for gravitational waves from compact binary coalescences using particle swarm optimization

While a fully-coherent all-sky search is known to be optimal for detecting gravitational wave signals from compact binary coalescences, its high computational cost has limited current searches to less sensitive coincidence-based schemes. Following up on previous work that has demonstrated the effectiveness of Particle Swarm Optimization in reducing the computational cost of this search, we present an implementation that achieves near real-time computational speed. This is achieved by combining the search efficiency of PSO with a significantly revised and optimized numerical implementation of the underlying mathematical formalism along with additional multi-threaded parallelization layers in a distributed computing framework. For a network of four second-generation detectors with $60$~min data from each, the runtime of the implementation presented here ranges between $\approx 1.4$ to $\approx 0.5$ times the data duration for network signal-to-noise ratios (SNRs) of $\gtrsim 10$ and $\gtrsim 12$, respectively. The reduced runtimes are obtained with small to negligible losses in detection sensitivity: for a false alarm rate of $\simeq 1$~event per year in Gaussian stationary noise, the loss in detection probability is $\leq 5\%$ and $\leq 2\%$ for SNRs of $10$ and $12$, respectively. Using the fast implementation, we are able to quantify frequentist errors in parameter estimation for signals in the double neutron star mass range using a large number of simulated data realizations. A clear dependence of parameter estimation errors and detection sensitivity on the condition number of the network antenna pattern matrix is revealed. Combined with previous work, this paper securely establishes the effectiveness of PSO-based fully-coherent all-sky search across the entire binary inspiral mass range that is relevant to ground-based detectors.

preprint2016arXiv

Coherent network analysis for continuous gravitational wave signals in a pulsar timing array: Pulsar phases as extrinsic parameters

Supermassive black hole binaries are one of the primary targets for gravitational wave searches using pulsar timing arrays. Gravitational wave signals from such systems are well represented by parametrized models, allowing the standard Generalized Likelihood Ratio Test (GLRT) to be used for their detection and estimation. However, there is a dichotomy in how the GLRT can be implemented for pulsar timing arrays: there are two possible ways in which one can split the set of signal parameters for semi-analytical and numerical extremization. The straightforward extension of the method used for continuous signals in ground-based gravitational wave searches, where the so-called pulsar phase parameters are maximized numerically, was addressed in an earlier paper (Wang et al. 2014). In this paper, we report the first study of the performance of the second approach where the pulsar phases are maximized semi-analytically. This approach is scalable since the number of parameters left over for numerical optimization does not depend on the size of the pulsar timing array. Our results show that, for the same array size (9 pulsars), the new method performs somewhat worse in parameter estimation, but not in detection, than the previous method where the pulsar phases were maximized numerically. The origin of the performance discrepancy is likely to be in the ill-posedness that is intrinsic to any network analysis method. However, scalability of the new method allows the ill-posedness to be mitigated by simply adding more pulsars to the array. This is shown explicitly by taking a larger array of pulsars.

preprint2016arXiv

Detection and localization of continuous gravitational waves with pulsar timing arrays: the role of pulsar terms

A pulsar timing array is a Galactic-scale detector of nanohertz gravitational waves (GWs). Its target signals contain two components: the `Earth term' and the `pulsar term' corresponding to GWs incident on the Earth and pulsar respectively. In this work we present a Frequentist method for the detection and localization of continuous waves that takes into account the pulsar term and is significantly faster than existing methods. We investigate the role of pulsar terms by comparing a full-signal search with an Earth-term-only search for non-evolving black hole binaries. By applying the method to synthetic data sets, we find that (i) a full-signal search can slightly improve the detection probability (by about five percent); (ii) sky localization is biased if only Earth terms are searched for and the inclusion of pulsar terms is critical to remove such a bias; (iii) in the case of strong detections (with signal-to-noise ratio $\gtrsim$ 30), it may be possible to improve pulsar distance estimation through GW measurements.

preprint2014arXiv

A coherent method for the detection and estimation of continuous gravitational wave signals using a pulsar timing array

The use of a high precision pulsar timing array is a promising approach to detecting gravitational waves in the very low frequency regime ($10^{-6} -10^{-9}$ Hz) that is complementary to the ground-based efforts (e.g., LIGO, Virgo) at high frequencies ($\sim 10 -10^3$ Hz) and space-based ones (e.g., LISA) at low frequencies ($10^{-4} -10^{-1}$ Hz). One of the target sources for pulsar timing arrays are individual supermassive black hole binaries that are expected to form in galactic mergers. In this paper, a likelihood based method for detection and estimation is presented for a monochromatic continuous gravitational wave signal emitted by such a source. The so-called pulsar terms in the signal that arise due to the breakdown of the long-wavelength approximation are explicitly taken into account in this method. In addition, the method accounts for equality and inequality constraints involved in the semi-analytical maximization of the likelihood over a subset of the parameters. The remaining parameters are maximized over numerically using Particle Swarm Optimization. Thus, the method presented here solves the monochromatic continuous wave detection and estimation problem without invoking some of the approximations that have been used in earlier studies.

preprint2011arXiv

Teaching introductory undergraduate Physics using commercial video games

Commercial video games are increasingly using sophisticated physics simulations to create a more immersive experience for players. This also makes them a powerful tool for engaging students in learning physics. We provide some examples to show how commercial off-the-shelf games can be used to teach specific topics in introductory undergraduate physics. The examples are selected from a course taught predominantly through the medium of commercial video games.

preprint2010arXiv

Particle Swarm Optimization and gravitational wave data analysis: Performance on a binary inspiral testbed

The detection and estimation of gravitational wave (GW) signals belonging to a parameterized family of waveforms requires, in general, the numerical maximization of a data-dependent function of the signal parameters. Due to noise in the data, the function to be maximized is often highly multi-modal with numerous local maxima. Searching for the global maximum then becomes computationally expensive, which in turn can limit the scientific scope of the search. Stochastic optimization is one possible approach to reducing computational costs in such applications. We report results from a first investigation of the Particle Swarm Optimization (PSO) method in this context. The method is applied to a testbed motivated by the problem of detection and estimation of a binary inspiral signal. Our results show that PSO works well in the presence of high multi-modality, making it a viable candidate method for further applications in GW data analysis.

preprint1999arXiv

Detecting an association between Gamma Ray and Gravitational Wave Bursts

If $γ$-ray bursts (GRBs) are accompanied by gravitational wave bursts (GWBs) the correlated output of two gravitational wave detectors evaluated in the moments just prior to a GRB will differ from that evaluated at times not associated with a GRB. We can test for this difference independently of any model of the GWB signal waveform. If we invoke a model for the GRB source population and GWB radiation spectral density we can find a confidence interval or upper limit on the root-mean-square GWB signal amplitude in the detector waveband. To illustrate we adopt a simple, physically motivated model and estimate that initial LIGO detector observations coincident with 1000 GRBs could lead us to exclude, with 95% confidence, associated GWBs with $h_{RMS} \gtrsim 1.7 \times 10^{-22}$. This result does not require the detector noise be Gaussian or that any inter-detector correlated noise be measured or measurable; it does not require advanced or a priori knowledge of the source waveform; and the limits obtained on the wave-strength improve with the number of observed GRBs.