Researcher profile

Shantanu Desai

Shantanu Desai contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
33works
0followers
10topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

33 published item(s)

preprint2025arXiv

Effects of Coronal Mass Ejection on PSR J1022+1001 and Possible Mode Change of PSR J2145-0750 in the InPTA DR2

The Indian Pulsar Timing Array (InPTA) has recently published its second data release (DR2), comprising the timing analysis of seven years of data on 27 millisecond pulsars (MSPs), observed simultaneously in the 300-500 MHz (band 3) and 1260-1460 MHz (band 5), using the upgraded Giant Metrewave Radio Telescope (uGMRT). The low-frequency data, particularly in band 3, is highly sensitive to propagation effects such as dispersion measure (DM) fluctuations, which can be imprints of some astrophysical phenomena (scientific outliers). Here, we analyze the two outliers of possible astrophysical origin coming from the band 3 DM time series of two pulsars: PSR J1022+1001, with an ecliptic latitude of -0.06 degree, and PSR J2145-0750, one of the brightest MSPs, with multi-component profile morphology. Our study reveals compelling evidence for a coronal mass ejection (CME) event traced in the data of PSR J1022+1001, and reports evidence for a potential mode-changing event in PSR J2145-0750. By contrasting these two cases, we show that DM fluctuations due to CME interacions and intrinsic mode-changing events produce distinct observational signatures, enabling a physically informed classification of scientific outliers in PTA datasets. Extending the analyses presented here to the full sample of InPTA-DR2 pulsars is expected to reveal additional CME events, and possible mode-changing events. Such detections will not only improve our understanding of solar and pulsar magnetospheric plasma interactions but will also enable more accurate modelling of DM variations, leading to improved pulsar timing solutions, which are crucial for high-precision Pulsar Timing Array (PTA) science.

preprint2025arXiv

Testing the cosmic distance duality relation using model-independent approach

In this work, we test the cosmic distance duality relation (CDDR) using the arbitrary redshift pivot Padé-(2,1) expansion methodology developed in arXiv:2509.16196. This approach allows us to constrain cosmography parameters and test CDDR at any redshift. Further, it does not rely on data reconstructions or extrapolations of the cosmography parameters to higher redshifts. We employ observational data from the Dark Energy Spectroscopic Instrument (DESI) Baryon Acoustic Oscillation dataset, cosmic chronometers (CC), and Type Ia supernovae from the Pantheon Plus (PP) and Dark Energy Survey Year 5 (DESY5) compilations. We find no significant deviations from the standard CDDR relation in the range $0\lesssim z \lesssim 1$ when considering DESI$+r_d$ dataset in combination with PP$+$CC and DESY5$+$CC datasets. However, on imposing a Gaussian prior on $M_B \in \mathcal{N}(-19.253, 0.027)$ (instead of treating it as a free parameter) in the dataset combination PP$+$CC, we find CDDR violation at a level of $(3-5)σ.$

preprint2024arXiv

The Gravitational Lensing Imprints of DES Y3 Superstructures on the CMB: A Matched Filtering Approach

$ $Low density cosmic voids gravitationally lens the cosmic microwave background (CMB), leaving a negative imprint on the CMB convergence $κ$. This effect provides insight into the distribution of matter within voids, and can also be used to study the growth of structure. We measure this lensing imprint by cross-correlating the Planck CMB lensing convergence map with voids identified in the Dark Energy Survey Year 3 data set, covering approximately 4,200 deg$^2$ of the sky. We use two distinct void-finding algorithms: a 2D void-finder which operates on the projected galaxy density field in thin redshift shells, and a new code, Voxel, which operates on the full 3D map of galaxy positions. We employ an optimal matched filtering method for cross-correlation, using the MICE N-body simulation both to establish the template for the matched filter and to calibrate detection significances. Using the DES Y3 photometric luminous red galaxy sample, we measure $A_κ$, the amplitude of the observed lensing signal relative to the simulation template, obtaining $A_κ= 1.03 \pm 0.22$ ($4.6σ$ significance) for Voxel and $A_κ= 1.02 \pm 0.17$ ($5.9σ$ significance) for 2D voids, both consistent with $Λ$CDM expectations. We additionally invert the 2D void-finding process to identify superclusters in the projected density field, for which we measure $A_κ= 0.87 \pm 0.15$ ($5.9σ$ significance). The leading source of noise in our measurements is Planck noise, implying that future data from the Atacama Cosmology Telescope (ACT), South Pole Telescope (SPT) and CMB-S4 will increase sensitivity and allow for more precise measurements.

preprint2022arXiv

A test of galaxy cluster fundamental plane for the X-COP sample

We test the galaxy cluster fundamental plane using the X-COP sample of 12 clusters. The fundamental plane is given by the relation $T_X \propto M_s^α r_s^β$, where $T_X$, $M_s$, and $r_s$ correspond to the gas temperature, NFW halo mass, and scale radius, respectively. We did this analysis using two different temperatures: the error-weighted temperature in $(50-500)h^{-1}$ kpc as well as the mass-weighted temperature in the same range. With both these temperatures, we find a very tight fundamental plane with dispersion of about 0.02 dex. The best-fit values for $α$ and $β$ are in-between those expected from virial equilibrium and self-similarity solution for secondary infall and collapse, with $α$ being closer to the virial expectation. Our best-fit values are also consistent with a recent re-analyses of the fundamental plane for the CLASH sample, after excluding the hottest clusters.

preprint2022arXiv

A test of the evolution of gas depletion factor in galaxy clusters using strong gravitational lensing systems

In this letter, we discuss a new method to probe the redshift evolution of the gas depletion factor, i.e. the ratio by which the gas mass fraction of galaxy clusters is depleted with respect to the universal mean of baryon fraction. The dataset we use for this purpose consists of 40 gas mass fraction measurements measured at $r_{2500}$ using Chandra X-ray observations, strong gravitational lensing sub-samples obtained from SLOAN Lens ACS + BOSS Emission-line Lens Survey (BELLS) + Strong Legacy Survey SL2S + SLACS. For our analysis, the validity of cosmic distance duality relation is assumed. We find a mildly decreasing trend for the gas depletion factor as a function of redshift at about 2.7$σ$. This is the first result in literature which does not find a constant gas depletion factor as a function of redshift using gas mass fraction measurements at $r_{2500}$.

preprint2022arXiv

A test of the standard dark matter density evolution law using galaxy clusters and cosmic chronometers

In this letter, we implement a test of the standard law for the dark matter density evolution. For this purpose, only a flat universe and the validity of the FRW metric are assumed. A deformed dark matter density evolution law is considered, given by $ρ_c(z) \propto (1+z)^{3+ε}$, and constraints on $ε$ are obtained by using galaxy cluster gas mass fractions, and cosmic chronometers measurements. We find that $ε=0$ within 2$σ$ c.l., in full agreement with other recent analyses.

preprint2022arXiv

Characterization of the GRB prompt fundamental plane using Fermi-GBM data

A recent work has found a tight global relation between the GRB prompt fluence, peak flux (based on the optimum time scale determined from the Bayesian Blocks-based analysis) and duration using data from the Fermi Gamma-Ray burst monitor, which they have dubbed as "Fundamental Plane". We quantitatively characterize the tightness of this Fundamental Plane relation by calculating the scatter in dex. We also check for a fundamental plane using the peak flux over time scales of 64 ms, 256 ms and 1024 ms. For our analysis, we incorporate the uncertainties in the above observables and carried out both a PCA as well as regression-based analysis. We find that the scatter in the fundamental plane is 0.16-0.17 dex.

preprint2022arXiv

Combined significance of spatial coincidence of high energy neutrinos from PSR B1509-58 by Super-Kamiokande and MACRO

In their searches for astrophysical point sources of high energy neutrinos, both the Super-Kamiokande and MACRO neutrino detectors saw the largest angular excess from the same source, viz. PSR B1509-58. We estimate the probability for the observed number of events by {\it both} Super-Kamiokande and MACRO to be a chance coincidence due to atmospheric neutrino background. We find that this probability is about 0.4%, corresponding to 2.6$σ$ significance. We also propose some additional tests to ascertain if this excess corresponds to an astrophysical signal or is only a background event.

preprint2022arXiv

Cosmic Ray Rejection with Attention Augmented Deep Learning

Cosmic Ray (CR) hits are the major contaminants in astronomical imaging and spectroscopic observations involving solid-state detectors. Correctly identifying and masking them is a crucial part of the image processing pipeline, since it may otherwise lead to spurious detections. For this purpose, we have developed and tested a novel Deep Learning based framework for the automatic detection of CR hits from astronomical imaging data from two different imagers: Dark Energy Camera (DECam) and Las Cumbres Observatory Global Telescope (LCOGT). We considered two baseline models namely deepCR and Cosmic-CoNN, which are the current state-of-the-art learning based algorithms that were trained using Hubble Space Telescope (HST) ACS/WFC and LCOGT Network images respectively. We have experimented with the idea of augmenting the baseline models using Attention Gates (AGs) to improve the CR detection performance. We have trained our models on DECam data and demonstrate a consistent marginal improvement by adding AGs in True Positive Rate (TPR) at 0.01% False Positive Rate (FPR) and Precision at 95% TPR over the aforementioned baseline models for the DECam dataset. We demonstrate that the proposed AG augmented models provide significant gain in TPR at 0.01% FPR when tested on previously unseen LCO test data having images from three distinct telescope classes. Furthermore, we demonstrate that the proposed baseline models with and without attention augmentation outperform state-of-the-art models such as Astro-SCRAPPY, Maximask (that is trained natively on DECam data) and pre-trained ground-based Cosmic-CoNN. This study demonstrates that the AG module augmentation enables us to get a better deepCR and Cosmic-CoNN models and to improve their generalization capability on unseen data.

preprint2022arXiv

Do Pulsar and Fast Radio Burst dispersion measures obey Benford's law?

We check if the first significant digit of the dispersion measure of pulsars and Fast Radio Bursts (using the CHIME catalog) is consistent with the Benford distribution. We find a large disagreement with Benford's law with $χ^2$ close to 80 for 8 degrees of freedom for both these aforementioned datasets. This corresponds to a discrepancy of about 7$σ$. Therefore, we conclude that the dispersion measures of pulsars and FRBs do not obey Benford's law.

preprint2022arXiv

Low-frequency wideband timing of InPTA pulsars observed with the uGMRT

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.

preprint2022arXiv

Search for a distance-dependent Baryonic Tully-Fisher Relation at low redshifts

A recent work (arXiv:2104.14481) has found a statistically significant transition in the Baryonic Tully-Fisher relation (BTFR) using low redshift data ($z<0.1$), with the transitions occurring at about 9 and 17 Mpc. Motivated by this finding, we carry out a variant of this analysis by fitting the data to an augmented BTFR, where both the exponent as well as normalization constant vary as a function of distance. We find that both the exponent and normalization constant show only a marginal variation with distance, and are consistent with a constant value, to within $2σ$. We also checked to see if there is a statistically significant difference between the BTFR results after bifurcating the dataset at distances of 9 and 17 Mpc. We find that almost all the sets of subsamples obey the BTFR with $χ^2$/dof close to 1 and the best-fit parameters consistent across the subsamples. Only the subsample with $D<17$ Mpc shows a marginal discrepancy (at $1.75σ$) with respect to the BTFR. Therefore, we do not find any evidence for statistically significant differences in the BTFR at distances of 9 and 17 Mpc.

preprint2022arXiv

Search for Cosmological time dilation from Gamma-Ray Bursts -- A 2021 status update

We carry out a search for signatures of cosmological time dilation in the light curves of Gamma Ray Bursts (GRBs), detected by the Neil Gehrels Swift Observatory. For this purpose, we calculate two different durations ($T_{50}$ and $T_{90}$) for a sample of 247 GRBs in the fixed rest frame energy interval of 140-350 keV, similar to Zhang et al. We then carry out a power law-based regression analysis between the durations and redshifts. This search is done using both the unbinned as well as the binned data, where both the weighted mean and the geometric mean was used. For each analysis, we also calculate the intrinsic scatter to determine the tightness of the relation. We find that weighted mean-based binned data for long GRBs and the geometric mean-based binned data is consistent with the cosmological time dilation signature, whereas the analyses using unbinned durations show a very large scatter. We also make our analysis codes and the procedure for obtaining the light curves and estimation of $T_{50}$/$T_{90}$ publicly available.

preprint2022arXiv

SOAR/Goodman Spectroscopic Assessment of Candidate Counterparts of the LIGO-Virgo Event GW190814

On 2019 August 14 at 21:10:39 UTC, the LIGO/Virgo Collaboration (LVC) detected a possible neutron star-black hole merger (NSBH), the first ever identified. An extensive search for an optical counterpart of this event, designated GW190814, was undertaken using the Dark Energy Camera (DECam) on the 4m Victor M. Blanco Telescope at the Cerro Tololo Inter-American Observatory. Target of Opportunity interrupts were issued on 8 separate nights to observe 11 candidates using the 4.1m Southern Astrophysical Research (SOAR) telescope&#39;s Goodman High Throughput Spectrograph in order to assess whether any of these transients was likely to be an optical counterpart of the possible NSBH merger. Here, we describe the process of observing with SOAR, the analysis of our spectra, our spectroscopic typing methodology, and our resultant conclusion that none of the candidates corresponded to the gravitational wave merger event but were all instead other transients. Finally, we describe the lessons learned from this effort. Application of these lessons will be critical for a successful community spectroscopic follow-up program for LVC observing run 4 (O4) and beyond.

preprint2022arXiv

Two Dimensional Clustering of Gamma-Ray Bursts using durations and hardness

Gamma-Ray Bursts (GRBs) have been traditionally divided into two categories: &#34;short&#34; and &#34;long&#34; with durations less than and greater than two seconds, respectively. However, there is a lot of literature (with conflicting results) regarding the existence of a third intermediate class. To investigate this issue, we carry out a two-dimensional classification using the GRB hardness and duration, and also incorporating the uncertainties in both the variables, by using an extension of Gaussian Mixture Model called Extreme Deconvolution (XDGMM). We carry out this analysis on datasets from two detectors, viz. BATSE and Fermi-GBM. We consider the duration and hardness features in log-scale for each of these datasets and determine the best-fit parameters using XDGMM. This is followed by information theoretic criterion-based tests (AIC and BIC) to determine the optimum number of classes. For BATSE, we find that both AIC and BIC show preference for two components with close to decisive and decisive significance, respectively. For Fermi-GBM, AIC shows preference for three components with decisive significance, whereas BIC does not find any significant difference between two and three components. Our analysis codes have been made publicly available.

preprint2021arXiv

Galaxy Morphology Classification using Neural Ordinary Differential Equations

We introduce a continuous depth version of the Residual Network (ResNet) called Neural ordinary differential equations (NODE) for the purpose of galaxy morphology classification. We carry out a classification of galaxy images from the Galaxy Zoo 2 dataset, consisting of five distinct classes, and obtained an accuracy between 91-95\%, depending on the image class. We train NODE with different numerical techniques such as adjoint and Adaptive Checkpoint Adjoint (ACA) and compare them against ResNet. While ResNet has certain drawbacks, such as time consuming architecture selection (e.g. the number of layers) and the requirement of a large dataset needed for training, NODE can overcome these limitations. Through our results, we show that that the accuracy of NODE is comparable to ResNet, and the number of parameters used is about one-third as compared to ResNet, thus leading to a smaller memory footprint, which would benefit next generation surveys.

preprint2021arXiv

The Gravity Collective: A Search for the Electromagnetic Counterpart to the Neutron Star-Black Hole Merger GW190814

We present optical follow-up imaging obtained with the Katzman Automatic Imaging Telescope, Las Cumbres Observatory Global Telescope Network, Nickel Telescope, Swope Telescope, and Thacher Telescope of the LIGO/Virgo gravitational wave (GW) signal from the neutron star-black hole (NSBH) merger GW190814. We searched the GW190814 localization region (19 deg$^{2}$ for the 90th percentile best localization), covering a total of 51 deg$^{2}$ and 94.6% of the two-dimensional localization region. Analyzing the properties of 189 transients that we consider as candidate counterparts to the NSBH merger, including their localizations, discovery times from merger, optical spectra, likely host-galaxy redshifts, and photometric evolution, we conclude that none of these objects are likely to be associated with GW190814. Based on this finding, we consider the likely optical properties of an electromagnetic counterpart to GW190814, including possible kilonovae and short gamma-ray burst afterglows. Using the joint limits from our follow-up imaging, we conclude that a counterpart with an $r$-band decline rate of 0.68 mag day$^{-1}$, similar to the kilonova AT 2017gfo, could peak at an absolute magnitude of at most $-17.8$ mag (50% confidence). Our data are not constraining for &#39;&#39;red&#39;&#39; kilonovae and rule out &#39;&#39;blue&#39;&#39; kilonovae with $M>0.5 M_{\odot}$ (30% confidence). We strongly rule out all known types of short gamma-ray burst afterglows with viewing angles $<$17$^{\circ}$ assuming an initial jet opening angle of $\sim$$5.2^{\circ}$ and explosion energies and circumburst densities similar to afterglows explored in the literature. Finally, we explore the possibility that GW190814 merged in the disk of an active galactic nucleus, of which we find four in the localization region, but we do not find any candidate counterparts among these sources.

preprint2021arXiv

Variational Inference as an alternative to MCMC for parameter estimation and model selection

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&#39;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.

preprint2020arXiv

An independent search for annual modulation and its significance in ANAIS-112 data

We perform an independent search for sinusoidal-based modulation in the recently released ANAIS-112 data, which could be induced by dark matter scatterings. We then evaluate this hypothesis against the null hypothesis that the data contains only background, using four different model comparison techniques. These include frequentist, Bayesian, and two information theory-based criteria (AIC and BIC). This analysis was done on both the residual data (by subtracting the exponential fit obtained from the ANAIS-112 Collaboration) as well as the total (non-background subtracted) data. We find that according to the Bayesian model comparison test, the null hypothesis of no modulation is decisively favored over a cosine-based annual modulation for the non-background subtracted dataset in 2-6 keV energy range. None of the other model comparison tests decisively favor any one hypothesis over another. This is the first application of Bayesian and information theory techniques to test the annual modulation hypothesis in ANAIS-112 data, extending our previous work on the DAMA/LIBRA and COSINE-100 data. Our analysis codes have also been made publicly available.

preprint2020arXiv

Constraints on variation of fine structure constant from joint SPT-SZ and XMM-Newton observations

We search for a variation of the electromagnetic fine structure constant ($α\equiv e^2/\hbar c$) using a sample of 58 SZ selected clusters in the redshift range ($0.2<z<1.5$) detected by the South Pole Telescope, along with X-ray measurements using the XMM-Newton observatory. We use the ratio of the integrated SZ Compto-ionization to its X-ray counterpart as our observable for this search. We first obtain a model-independent constraint on $α$ of about 0.7%, using the fact that the aforementioned ratio is constant as a function of redshift. We then look for a logarithmic dependence of $α$ as a function of redshift: $Δα/α= -γ\ln(1+z)$, as this is predicted by runaway dilaton models. We find that $γ$ = $-0.046 \pm 0.1$, which indicates that there is no logarithmic variation of $α$ as a function of redshift. We also search for a dipole variation of the fine structure constant using the same cluster sample. We do not find any evidence for such a spatial variation.

preprint2020arXiv

Dark Energy Survey Identification of A Low-Mass Active Galactic Nucleus at Redshift 0.823 from Optical Variability

We report the identification of a low-mass AGN, DES J0218$-$0430, in a redshift $z = 0.823$ galaxy in the Dark Energy Survey (DES) Supernova field. We select DES J0218$-$0430 as an AGN candidate by characterizing its long-term optical variability alone based on DES optical broad-band light curves spanning over 6 years. An archival optical spectrum from the fourth phase of the Sloan Digital Sky Survey shows both broad Mg II and broad H$β$ lines, confirming its nature as a broad-line AGN. Archival XMM-Newton X-ray observations suggest an intrinsic hard X-ray luminosity of $L_{\rm 2-12\,keV}\sim7.6\pm0.4\times10^{43}$ erg s$^{-1}$, which exceeds those of the most X-ray luminous starburst galaxies, in support of an AGN driving the optical variability. Based on the broad H$β$ from SDSS spectrum, we estimate a virial BH mass of $M_{\bullet}\approx10^{6.43}$-$10^{6.72}M_{\odot}$ (with the error denoting 1$σ$ statistical uncertainties only), consistent with the estimation from OzDES, making it the lowest mass AGN with redshift $>$ 0.4 detected in optical. We estimate the host galaxy stellar mass to be $M_{\ast}\sim10^{10.5\pm0.3}M_{\odot}$ based on modeling the multi-wavelength spectral energy distribution. DES J0218$-$0430 extends the $M_{\bullet}$-$M_{\ast}$ relation observed in luminous AGNs at $z\sim1$ to masses lower than being probed by previous work. Our work demonstrates the feasibility of using optical variability to identify low-mass AGNs at higher redshift in deeper synoptic surveys with direct implications for the upcoming Legacy Survey of Space and Time at Vera C. Rubin Observatory.

preprint2020arXiv

Dark Energy Survey Year 1 Results: Cosmological Constraints from Cluster Abundances and Weak Lensing

We perform a joint analysis of the counts and weak lensing signal of redMaPPer clusters selected from the Dark Energy Survey (DES) Year 1 dataset. Our analysis uses the same shear and source photometric redshifts estimates as were used in the DES combined probes analysis. Our analysis results in surprisingly low values for $S_8 =σ_8(Ω_{\rm m}/0.3)^{0.5}= 0.65\pm 0.04$, driven by a low matter density parameter, $Ω_{\rm m}=0.179^{+0.031}_{-0.038}$, with $σ_8-Ω_{\rm m}$ posteriors in $2.4σ$ tension with the DES Y1 3x2pt results, and in $5.6σ$ with the Planck CMB analysis. These results include the impact of post-unblinding changes to the analysis, which did not improve the level of consistency with other data sets compared to the results obtained at the unblinding. The fact that multiple cosmological probes (supernovae, baryon acoustic oscillations, cosmic shear, galaxy clustering and CMB anisotropies), and other galaxy cluster analyses all favor significantly higher matter densities suggests the presence of systematic errors in the data or an incomplete modeling of the relevant physics. Cross checks with X-ray and microwave data, as well as independent constraints on the observable--mass relation from SZ selected clusters, suggest that the discrepancy resides in our modeling of the weak lensing signal rather than the cluster abundance. Repeating our analysis using a higher richness threshold ($λ\ge 30$) significantly reduces the tension with other probes, and points to one or more richness-dependent effects not captured by our model.

preprint2020arXiv

Dust Reverberation Mapping in Distant Quasars from Optical and Mid-Infrared Imaging Surveys

The size of the dust torus in Active Galactic Nuclei (AGN) and their high-luminosity counterparts, quasars, can be inferred from the time delay between UV/optical accretion disk continuum variability and the response in the mid-infrared (MIR) torus emission. This dust reverberation mapping (RM) technique has been successfully applied to $\sim 70$ $z\lesssim 0.3$ AGN and quasars. Here we present first results of our dust RM program for distant quasars covered in the SDSS Stripe 82 region combining $\sim 20$-yr ground-based optical light curves with 10-yr MIR light curves from the WISE satellite. We measure a high-fidelity lag between W1-band (3.4 $μ$m) and $g$ band for 587 quasars over $0.3\lesssim z\lesssim 2$ ($\left<z\right>\sim 0.8$) and two orders of magnitude in quasar luminosity. They tightly follow (intrinsic scatter $\sim 0.17$ dex in lag) the IR lag-luminosity relation observed for $z<0.3$ AGN, revealing a remarkable size-luminosity relation for the dust torus over more than four decades in AGN luminosity, with little dependence on additional quasar properties such as Eddington ratio and variability amplitude. This study motivates further investigations in the utility of dust RM for cosmology, and strongly endorses a compelling science case for the combined 10-yr Vera C. Rubin Observatory Legacy Survey of Space and Time (optical) and 5-yr Nancy Grace Roman Space Telescope 2$μ$m light curves in a deep survey for low-redshift AGN dust RM with much lower luminosities and shorter, measurable IR lags. The compiled optical and MIR light curves for 7,384 quasars in our parent sample are made public with this work.

preprint2020arXiv

Galaxy cluster hydrostatic masses using Tolman-Oppenheimer-Volkoff equation

Motivated by previous studies in literature about the potential importance of relativistic corrections to galaxy cluster hydrostatic masses, we calculate the masses of 12 relaxed clusters (with Chandra X-ray data) using the Tolman-Oppenheimer-Volkov (TOV) equation of hydrostatic equilibrium and the ideal gas equation of state. Analytical formulae for gas density and temperature profiles for these clusters, previously derived by Vikhlinin et al (astro-ph/0507092) were used to obtain these masses. We compare the TOV-based masses with those obtained using the corresponding Newtonian equation of hydrostatic equilibrium. We find that the fractional relative difference between the two masses are negligible, corresponding to $\sim \mathcal{O}(10^{-5})$.

preprint2020arXiv

Generalized Lomb-Scargle analysis of $\rm{^{36}Cl}$ decay rate measurements at PTB and BNL

Recently, Pomme et al did an analysis of $\rm{^{36}Cl} $ radioactive decay data from measurements at the Physikalisch-Technische Bundesanstalt (PTB), in order to verify the claims by Sturrock and collaborators of an influence on beta-decay rates measured at Brookehaven National Lab (BNL) due to the rotation-induced modulation of the solar neutrino flux. Their analysis excluded any sinusoidal modulations in the frequency range from 0.2-20/year. We carry out an independent analysis of the same PTB and BNL data using the generalized Lomb-Scargle periodogram, to look for any statistically significant peaks in the range from 0 to 14 per year, and by evaluating the significance of every peak using multiple methods. Our results for the PTB data are in agreement with those by Pomme et al. For BNL data, we do find peaks at some of the same frequencies as Sturrock et al, but the significance is much lower. All our analysis codes and datasets have been made publicly available.

preprint2020arXiv

Model Comparison of $Λ$CDM vs $R_h=ct$ using Cosmic Chronometers

In 2012, Bilicki and Seikel (arXiv:1206.5130) showed that $H(z)$ data reconstructed using Gaussian Process Regression from cosmic chronometers and BAO, conclusively rules out the $R_h=ct$ model. These results were disputed by Melia and collaborators in two different works (arXiv:1304.1802 and arXiv:1802.02255), who showed using both an unbinned analysis and Gaussian Process reconstructed $H(z)$ data from chronometers, that $R_h=ct$ is favored over $Λ$CDM model. To resolve this imbroglio, we carry out model comparison of $Λ$CDM versus $R_h=ct$ by independently reproducing the above claims using latest chronometer data. We perform model selection between these two models using Bayesian model comparison. We find that no one model between $Λ$CDM and $R_h=ct$ is decisively favored when uniform priors on $Λ$CDM parameters are used. However, if we use priors centered around the Planck best-fit values, then $Λ$CDM is very strongly preferred over $R_h=ct$.

preprint2020arXiv

Model Comparison tests of modified gravity from the Eöt-Wash experiment

Perivolaropoulos et al (arXiv:1904.09462) (P19) have argued that the residual torque data in the Eöt-Wash experiment is consistent with an oscillating signal. This could either be a signature of non-local modified gravity theories or some other systematic error in the data. We independently assess the viability of such an oscillating signal in the same data using Bayesian and information theoretical criterion, to complement the frequentist analysis in P19. We fit this data to three different parametrizations (an offset Newtonian, Yukawa, and an oscillating model), and assess the significance of the oscillating model using AIC, BIC, WAIC, and Bayes factor. All these techniques provide decisive evidence for the oscillating model compared to the Newtonian model, provided the phase is fixed to the same value as in P19. If the phase is allowed to vary, then significance from BIC, WAIC, and Bayes based tests reduces to strong evidence, whereas only AIC still shows decisive evidence. Our analysis codes have been made publicly available.

preprint2020arXiv

Robust model comparison tests of DAMA/LIBRA annual modulation

We evaluate the statistical significance of the DAMA/LIBRA claims for annual modulation using three independent model comparison techniques, viz frequentist, information theory, and Bayesian analysis. We fit the data from the DAMA/LIBRA experiment to both cosine and a constant model, and carry out model comparison by choosing the constant model as the null hypothesis. For the frequentist test, we invoke Wilk&#39;s theorem and calculate the significance using $Δχ^2$ between the two models. For information theoretical tests, we calculate the difference in Akaike Information Criterion (AIC) and Bayesian Information criterion (BIC) between the two models. We also compare the two models in a Bayesian context by calculating the Bayes factor. We also search for higher harmonics in the DAMA/LIBRA data using generalized Lomb-Scargle periodogram. We finally test the sensitivity of these model comparison techniques in discriminating between pure noise and a cosine signal using synthetic data. This is the first proof of principles application of AIC, BIC as well as Bayes factor to the DAMA data. All our analysis codes along with the data used in this work have been made publicly available.

preprint2020arXiv

Scaling relations for dark matter core density and radius from Chandra X-ray cluster sample

A large number of studies have found that the dark matter surface density, given by the product of the dark matter core radius ($r_c$) and core density ($ρ_c$) is approximately constant for a wide range of galaxy systems. However, there has been only one systematic study of this {\it ansatz} for galaxy clusters by Chan (arXiv:1403.4352), who found that the surface density for clusters is not constant and $ρ_c \sim r_c^{-1.46}$. We carry out this test for an X-ray sample of 12 relaxed clusters from Chandra observations, implementing the same procedure as Chan, but also accounting for the gas and star mass. We find that $ρ_c \propto r_c^{-1.08 \pm 0.055}$, with an intrinsic scatter of about 18%. Therefore, the dark matter surface density for our cluster data shows deviations from a constant value at only about 1.4$σ$.

preprint2020arXiv

The Curious Case of PHL 293B: A Long-Lived Transient in a Metal-Poor Blue Compact Dwarf Galaxy

We report on small-amplitude optical variability and recent dissipation of the unusually persistent broad emission lines in the blue compact dwarf galaxy PHL 293B. The galaxy&#39;s unusual spectral features (P Cygni-like profiles with $\sim$800 km s$^{-1}$ blueshifted absorption lines) have resulted in conflicting interpretations of the nature of this source in the literature. However, analysis of new Gemini spectroscopy reveals the broad emission has begun to fade after being persistent for over a decade prior. Precise difference imaging light curves constructed with the Sloan Digital Sky Survey and the Dark Energy Survey reveal small-amplitude optical variability of $\sim$0.1 mag in the g band offset by $100\pm21$ pc from the brightest pixel of the host. The light curve is well-described by an active galactic nuclei (AGN)-like damped random walk process. However, we conclude that the origin of the optical variability and spectral features of PHL 293B is due to a long-lived stellar transient, likely a Type IIn supernova or non-terminal outburst, mimicking long-term AGN-like variability. This work highlights the challenges of discriminating between scenarios in such extreme environments, relevant to searches for AGNs in dwarf galaxies. This is the second long-lived transient discovered in a blue compact dwarf, after SDSS1133. Our result implies such long-lived stellar transients may be more common in metal-deficient galaxies. Systematic searches for low-level variability in dwarf galaxies will be possible with the upcoming Legacy Survey of Space and Time at Vera C. Rubin Observatory.

preprint2019arXiv

A meta-analysis of neutron lifetime measurements

We calculated the median as well as weighted mean central estimates for the neutron lifetime, from a subset of measurements compiled in the 2019 update of the Particle Data Group (PDG). We then reconstruct the error distributions for the residuals using three different central estimates and then check for consistency with a Gaussian distribution. We find that although the error distributions using the weighted mean as well as median estimate are consistent with a Gaussian distribution, the Student&#39;s $t$ and Cauchy distribution provide a better fit. This median statistic estimate of the neutron lifetime from these measurements is given by $881.5 \pm 0.47$ seconds. This can be used as an alternate estimate of the neutron lifetime. We also note that the discrepancy between beam and bottle-based measurements using median statistics of the neutron lifetime persists with a significance between 4-8$σ$, depending on which combination of measurements is used.

preprint2019arXiv

Recent bounds on graviton mass using galaxy clusters

Although galaxy clusters have proved to be wonderful laboratories for testing a plethora of modified gravity theories and other exotic alternatives to $Λ$CDM, until a few years ago, there was only one paper (from 1974), which obtained a limit on graviton mass of $\mathcal{O}(10^{-29})$ eV with clusters. To rectify this, in the last few years multiple works have obtained different bounds on graviton mass using single cluster data as well as stacking galaxy catalogs. We review these recent limits on graviton mass using galaxy clusters obtained using disparate methods.

preprint2019arXiv

Spectral Variability of a Sample of Extreme Variability Quasars and Implications for the MgII Broad-line Region

We present new Gemini/GMOS optical spectroscopy of 16 extreme variability quasars (EVQs) that dimmed by more than 1.5 mag in the $g$ band between the Sloan Digital Sky Survey (SDSS) and the Dark Energy Survey (DES) epochs (separated by a few years in the quasar rest frame). The quasar sample covers a redshift range of $0.5 < z < 2.1$. Nearly half of these EVQs brightened significantly (by more than 0.5 mag in the $g$ band) in a few years after reaching their previous faintest state, and some EVQs showed rapid (non-blazar) variations of greater than 1-2 mag on timescales of only months. Leveraging on the large dynamic range in continuum variability between the earlier SDSS and the new GMOS spectra, we explore the associated variations in the broad Mg II,$\lambda2798$ line, whose variability properties have not been well studied before. The broad Mg II flux varies in the same direction as the continuum flux, albeit with a smaller amplitude, which indicates at least some portion of Mg II is reverberating to continuum changes. However, the width (FWHM) of Mg II does not vary accordingly as continuum changes for most objects in the sample, in contrast to the case of the broad Balmer lines. Using the width of broad Mg II to estimate the black hole mass therefore introduces a luminosity-dependent bias.