Source author record

Varun Bhalerao

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

28works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

28 published item(s)

preprint2026arXiv

LUMOS : Linear programming Utility for Multi-messenger Optical Scheduling

The detection of gravitational-wave events by LIGO-Virgo-KAGRA has opened new avenues for multi-messenger astrophysics; however, electromagnetic counterparts remain elusive due to large localization uncertainties. Wide-field optical surveys like the Zwicky Transient Facility (ZTF) play a crucial role in follow-up, but efficient scheduling is essential. In this work, we present LUMOS, a Mixed Integer Linear Programming (MILP) approach that selects fields via a maximum coverage problem and schedules observations to maximize cumulative probability while respecting observability constraints. Using 1199 GW events from O4, we compare the LUMOS scheduler with gwemopt, showing an 84.7 percent higher mean cumulative probability and better performance in nearly all cases. While designed for ZTF, LUMOS's framework parallels the M4OPT toolkit for space missions, highlighting the broader applicability of MILP-based scheduling to both ground- and space-based follow-up.

preprint2023arXiv

High hard X-ray polarization in Cygnus X-1 confined to the intermediate hard state: evidence for a variable jet component

Cygnus X-1, the well-known accreting black hole system, exhibits several observational features hinting at an intricate interplay between the accretion disk, its atmosphere known as the corona and the putative relativistic jet. It has been extensively studied using all available observational methods, including using the newly available technique of sensitive X-ray polarimetry. X-ray polarization characteristics are distinct for coronal and jet emissions. The low X-ray polarization measured below $\sim$100 keV is understood as arising from the corona. In contrast, the high polarization measurements reported above $\sim$400 keV required a separate jet-dominated spectral component, which spectroscopy does not demonstrate conclusively. Here we report precise polarization measurements in the 100-380 keV region made during three different sub-classes of spectral states of the source using the CZTI instrument onboard {\em AstroSat}. A high polarization (23$\pm$4 \%) is found mainly in the Intermediate Hard State of the source, and the energy-resolved measurements smoothly connect the coronal and the jet regimes. When high polarization is observed, the simultaneous spectral data hints at a separate power law component above 100 keV. We examine the possible sources of this energy-dependent high polarization in Cygnus X-1.

preprint2022arXiv

Experimental verification of off-axis polarimetry with Cadmium Zinc Telluride detectors of AstroSat-CZT Imager

The Cadmium Zinc Telluride Imager (CZTI) onboard AstroSat consists of an array of a large number of pixellated CZT detectors capable of measuring the polarization of incident hard X-rays. The polarization measurement capability of CZTI for on-axis sources was experimentally confirmed before the launch. CZTI has yielded tantalizing results on the X-ray polarization of the Crab nebula and pulsar in the energy range of 100 - 380 keV. CZTI has also contributed to the measurement of prompt emission polarization for several Gamma-Ray Bursts (GRBs). However, polarization measurements of off-axis sources like GRBs are challenging. It is vital to experimentally calibrate the CZTI sensitivity to off-axis sources to enhance the credence of the measurements. In this context, we report the verification of the off-axis polarimetric capability of pixellated CZT detectors through the controlled experiments carried out with a CZT detector similar to that used in CZTI and extensive Geant4 simulations of the experimental set-up. Our current results show that the CZT detectors can be used to measure the polarization of bright GRBs up to off-axis angles of ~60 degrees. However, at incidence angles between 45-60 degrees, there might be some systematic effects which needs to be taken into account while interpreting the measured polarisation fraction.

preprint2022arXiv

Extending the energy range of AstroSat-CZTI up to 380 keV with Compton Spectroscopy

The CZTI (Cadmium Zinc Telluride Imager) onboard AstroSat is a high energy coded mask imager and spectrometer in the energy range of 20 - 100 keV. Above 100 keV, the dominance of Compton scattering cross-section in CZTI results in a significant number of 2-pixel Compton events and these have been successfully utilized for polarization analysis of Crab pulsar and nebula (and transients like Gamma-ray bursts) in 100 - 380 keV. These 2-pixel Compton events can also be used to extend the spectroscopic energy range of CZTI up to 380 keV for bright sources. However, unlike the spectroscopy in primary energy range, where simultaneous background measurement is available from masked pixels, Compton spectroscopy requires blank sky observation for background measurement. Background subtraction, in this case, is non-trivial because of the presence of both short-term and long-term temporal variations in the data, which depend on multiple factors like earth rotation and the effect of South Atlantic Anomaly (SAA) regions etc. We have developed a methodology of background selection and subtraction that takes into account for these effects. Here, we describe these background selection and subtraction techniques and validate them using spectroscopy of Crab in the extended energy range of 30 - 380 keV region, and compare the obtained spectral parameters with the INTEGRAL results. This new capability allows for the extension of the energy range of AstroSat spectroscopy and will also enable the simultaneous spectro-polarimetric study of other bright sources like Cygnus X-1.

preprint2022arXiv

GROWTH on S190426c II: GROWTH-India Telescope search for an optical counterpart with a custom image reduction and candidate vetting pipeline

S190426c / GW190426_152155 was the first probable neutron star - black hole merger candidate detected by the LIGO-Virgo Collaboration. We undertook a tiled search for optical counterparts of this event using the 0.7m GROWTH-India Telescope. Over a period of two weeks, we obtained multiple observations over a 22.1 deg^2 area, with a 17.5% probability of containing the source location. Initial efforts included obtaining photometry of sources reported by various groups, and a visual search for sources in all galaxies contained in the region. Subsequently, we have developed an image subtraction and candidate vetting pipeline with ~ 94% efficiency for transient detection. Processing the data with this pipeline, we find several transients, but none that are compatible with kilonova models. We present the details of our observations, working of our pipeline, results from the search, and our interpretations of the non-detections that will work as a pathfinder during the O4 run of LVK.

preprint2022arXiv

Hard X-ray polarization catalog for a 5-year sample of Gamma-Ray Bursts using AstroSat CZT-Imager

Cadmium Zinc Telluride Imager (CZTI) aboard AstroSat has been regularly detecting Gamma-Ray Bursts (GRBs) since its launch in 2015. Its sensitivity to polarization measurements at energies above 100 keV allows CZTI to attempt spectro-polarimetric studies of GRBs. Here, we present the first catalog of GRB polarization measurements made by CZTI during its first five years of operation. This presents the time integrated polarization measurements of the prompt emission of 20 GRBs in the energy range 100-600 keV. The sample includes the bright GRBs which were detected within an angle range of 0-60 degree and 120-180 degree where the instrument has useful polarization sensitivity and is less prone to systematics. We implement a few new modifications in the analysis to enhance polarimetric sensitivity of the instrument. Majority of the GRBs in the sample are found to possess less / null polarization across the total bursts' duration in contrast to a small fraction of five GRBs exhibiting high polarization. The low polarization across the bursts can be speculated to be either due to the burst being intrinsically weakly polarized or due to varying polarization angle within the burst even when it is highly polarized. In comparison to POLAR measurements, CZTI has detected a larger number of cases with high polarization. This may be a consequence of the higher energy window of CZTI observations which results in the sampling of smaller duration of burst emissions in contrast to POLAR, thereby, probing emissions of less temporal variations of polarization properties.

preprint2022arXiv

In search of short gamma-ray burst optical counterpart with the Zwicky Transient Facility

The Fermi Gamma-ray Burst Monitor (GBM) triggers on-board in response to $\sim$ 40 short gamma-ray bursts (SGRBs) per year; however, their large localization regions have made the search for optical counterparts a challenging endeavour. We have developed and executed an extensive program with the wide field of view of the Zwicky Transient Facility (ZTF) camera, mounted on the Palomar 48 inch Oschin telescope (P48), to perform target-of-opportunity (ToO) observations on 10 Fermi-GBM SGRBs during 2018 and 2020-2021. Bridging the large sky areas with small field of view optical telescopes in order to track the evolution of potential candidates, we look for the elusive SGRB afterglows and kilonovae (KNe) associated with these high-energy events. No counterpart has yet been found, even though more than 10 ground based telescopes, part of the Global Relay of Observatories Watching Transients Happen (GROWTH) network, have taken part in these efforts. The candidate selection procedure and the follow-up strategy have shown that ZTF is an efficient instrument for searching for poorly localized SGRBs, retrieving a reasonable number of candidates to follow-up and showing promising capabilities as the community approaches the multi-messenger era. Based on the median limiting magnitude of ZTF, our searches would have been able to retrieve a GW170817-like event up to $\sim$ 200 Mpc and SGRB afterglows to z = 0.16 or 0.4, depending on the assumed underlying energy model. Future ToOs will expand the horizon to z = 0.2 and 0.7 respectively.

preprint2022arXiv

India's first robotic eye for time domain astrophysics: the GROWTH-India telescope

We present the design and performance of the GROWTH-India telescope, a 0.7 m robotic telescope dedicated to time-domain astronomy. The telescope is equipped with a 4k back-illuminated camera giving a 0.82-degree field of view and sensitivity of m_g ~20.5 in 5-min exposures. Custom software handles observatory operations: attaining high on-sky observing efficiencies (>~ 80%) and allowing rapid response to targets of opportunity. The data processing pipelines are capable of performing PSF photometry as well as image subtraction for transient searches. We also present an overview of the GROWTH-India telescope's contributions to the studies of Gamma-ray Bursts, the electromagnetic counterparts to gravitational wave sources, supernovae, novae and solar system objects.

preprint2022arXiv

The long-active afterglow of GRB 210204A: Detection of the most delayed flares in a Gamma-Ray Burst

We present results from extensive broadband follow-up of GRB 210204A over the period of thirty days. We detect optical flares in the afterglow at 7.6 x 10^5 s and 1.1 x 10^6 s after the burst: the most delayed flaring ever detected in a GRB afterglow. At the source redshift of 0.876, the rest-frame delay is 5.8 x 10^5 s (6.71 d). We investigate possible causes for this flaring and conclude that the most likely cause is a refreshed shock in the jet. The prompt emission of the GRB is within the range of typical long bursts: it shows three disjoint emission episodes, which all follow the typical GRB correlations. This suggests that GRB 210204A might not have any special properties that caused late-time flaring, and the lack of such detections for other afterglows might be resulting from the paucity of late-time observations. Systematic late-time follow-up of a larger sample of GRBs can shed more light on such afterglow behaviour. Further analysis of the GRB 210204A shows that the late time bump in the light curve is highly unlikely due to underlying SNe at redshift (z) = 0.876 and is more likely due to the late time flaring activity. The cause of this variability is not clearly quantifiable due to the lack of multi-band data at late time constraints by the bad weather conditions. The flare of GRB 210204A is the latest flare detected to date.

preprint2021arXiv

Discovery and confirmation of the shortest gamma ray burst from a collapsar

Gamma-ray bursts (GRBs) are among the brightest and most energetic events in the universe. The duration and hardness distribution of GRBs has two clusters, now understood to reflect (at least) two different progenitors. Short-hard GRBs (SGRBs; T90 <2 s) arise from compact binary mergers, while long-soft GRBs (LGRBs; T90 >2 s) have been attributed to the collapse of peculiar massive stars (collapsars). The discovery of SN 1998bw/GRB 980425 marked the first association of a LGRB with a collapsar and AT 2017gfo/GRB 170817A/GW170817 marked the first association of a SGRB with a binary neutron star merger, producing also gravitational wave (GW). Here, we present the discovery of ZTF20abwysqy (AT2020scz), a fast-fading optical transient in the Fermi Satellite and the InterPlanetary Network (IPN) localization regions of GRB 200826A; X-ray and radio emission further confirm that this is the afterglow. Follow-up imaging (at rest-frame 16.5 days) reveals excess emission above the afterglow that cannot be explained as an underlying kilonova (KN), but is consistent with being the supernova (SN). Despite the GRB duration being short (rest-frame T90 of 0.65 s), our panchromatic follow-up data confirms a collapsar origin. GRB 200826A is the shortest LGRB found with an associated collapsar; it appears to sit on the brink between a successful and a failed collapsar. Our discovery is consistent with the hypothesis that most collapsars fail to produce ultra-relativistic jets.

preprint2021arXiv

Exploring Sub-MeV Sensitivity of AstroSat-CZTI for ON-axis Bright Sources

The Cadmium Zinc Telluride Imager (CZTI) onboard AstroSat is designed for hard X-ray imaging and spectroscopy in the energy range of 20 - 100 keV. The CZT detectors are of 5 mm thickness and hence have good efficiency for Compton interactions beyond 100 keV. The polarisation analysis using CZTI relies on such Compton events and have been verified experimentally. The same Compton events can also be used to extend the spectroscopy up to 380 keV. Further, it has been observed that about 20% pixels of the CZTI detector plane have low gain, and they are excluded from the primary spectroscopy. If these pixels are included, then the spectroscopic capability of CZTI can be extended up to 500 keV and further up to 700 keV with a better gain calibration in the future. Here we explore the possibility of using the Compton events as well as the low gain pixels to extend the spectroscopic energy range of CZTI for ON-axis bright X-ray sources. We demonstrate this technique using Crab observations and explore its sensitivity.

preprint2021arXiv

Initial Characterization of Active Transitioning Centaur, P/2019 LD2 (ATLAS), using Hubble, Spitzer, ZTF, Keck, APO and GROWTH Visible & Infrared Imaging and Spectroscopy

We present visible and mid-infrared imagery and photometry of temporary Jovian co-orbital comet P/2019 LD$_2$ taken with HST/WFC3, Spitzer/IRAC, the GROWTH telescope network, visible spectroscopy from Keck/LRIS and archival ZTF observations taken between 2019 April and 2020 August. Our observations indicate that the nucleus of LD$_2$ has a radius between 0.2-1.8 km assuming a 0.08 albedo and a coma dominated by $\sim$100$μ$ m-scale dust ejected at $\sim$1 m/s speeds with a $\sim$1'' jet pointing in the SW direction. LD$_2$ experienced a total dust mass loss of $\sim$10$^8$ kg at a loss rate of $\sim$6 kg/s with Af$ρ$/cross-section varying between $\sim$85 cm/125 km$^2$ and $\sim$200 cm/310 km$^2$ from 2019 April 9 to 2019 Nov 8. If the increase in Af$ρ$/cross-section remained constant, it implies LD$_2$'s activity began $\sim$2018 November when within 4.8 au of the Sun, implying the onset of H$_2$O sublimation. We measure CO/CO$_2$ gas production of $\lesssim$10$^{27}$ mol/s /$\lesssim$10$^{26}$ mol/s from our 4.5 $μ$m Spitzer observations, $g$-$r$ = 0.59$\pm$0.03, $r$-$i$ = 0.18$\pm$0.05, $i$-$z$ = 0.01$\pm$0.07 from GROWTH observations, H$_2$O gas production of $\lesssim$80 kg/s scaling from our estimated $C_2$ production of $Q_{C_2}\lesssim$7.5$\times10^{24}$ mol/s from Keck/LRIS spectroscopy. We determine that the long-term orbit of LD$_2$ is similar to Jupiter family comets having close encounters with Jupiter within $\sim$0.5 Hill radius in the last $\sim$3 y, within 0.8 Hill radius in $\sim$9 y. Additionally, 78.8$\%$ of our orbital clones are ejected from the Solar System within $1 \times 10^{6}$ years having a dynamical half-life of 3.4 $\times 10^5$ years.

preprint2021arXiv

SN 2020sck: deflagration in a carbon-oxygen white dwarf

We present optical UBVRI photometry and low-to-medium resolution spectroscopic observations of type Iax SN 2020sck spanning -5.5 d to +67 d from maximum light in the B-band. From the photometric analysis we find $Δm_{\rm{B}}$(15) = 2.03$\pm$0.05 mag and $M_{\rm{B}}$=-17.81$\pm$0.22 mag. Radiation diffusion model fit to the quasi-bolometric light curve indicates 0.13$\pm$0.02 $M_{\odot}$ of $^{56}$Ni and 0.34 $M_{\rm \odot}$ of ejecta are synthesized in the explosion. Comparing the observed quasi-bolometric light curve with angle-averaged bolometric light curve of three-dimensional pure deflagration explosion of $M_{\rm{ch}}$ carbon-oxygen white dwarf, we find agreement with a model in which 0.16 $M_{\odot}$ of $^{56}$Ni and 0.37 $M_{\odot}$ of ejecta is formed. By comparing the +1.4 day spectrum of SN 2020sck with synthetic spectrum generated using SYN++, we find absorption features due to C II, C III and O I. These are unburned materials in the explosion and indicate a C-O white dwarf. One dimensional radiative transfer modeling of the spectra with TARDIS shows higher density in the ejecta near the photosphere and a steep decrease in the outer layers with an ejecta composition dominated mostly by C, O, Si, Fe, and Ni. The star formation rate of the host galaxy computed from the luminosity of the H$α$ ($λ$6563) line is 0.09 $M_{\odot}$ yr$^{-1}$ indicating a relatively young stellar environment.

preprint2020arXiv

Characterization of Temporarily-Captured Minimoon 2020 CD$_3$ by Keck Time-resolved Spectrophotometry

We present time-resolved visible spectrophotometry of minimoon 2020 CD$_3$, the second asteroid known to become temporarily captured by the Earth-Moon system's gravitational field. The spectrophotometry was taken with Keck I/LRIS between wavelengths 434 nm and 912 nm in $B$, $g$, $V$, $R$, $I$ and RG850 filters as it was leaving the Earth-Moon system on 2020 March 23 UTC. The spectrophotometry of 2020 CD$_3$ most closely resembles the spectra of V-type asteroids and some Lunar rock samples with a reddish slope of ~18$\%$/100 nm between 434 nm and 761 nm corresponding to colors of $g$-$r$ = 0.62$\pm$0.08, $r$-$i$ = 0.21 $\pm$ 0.06 and an absorption band at ~900 nm corresponding to $i$-$z$ = -0.54$\pm$0.10. Combining our measured 31.9$\pm$0.1 absolute magnitude with an albedo of 0.35 typical for V-type asteroids, we determine 2020 CD$_3$'s diameter to be ~0.9$\pm$0.1 m making it the first minimoon and one of the smallest asteroids to be spectrally studied. We use our time-series photometry to detect periodic lightcurve variations with a $<$10$^{-4}$ false alarm probability corresponding to a lightcurve period of ~573 s and a lightcurve amplitude of ~1 mag implying 2020 CD$_3$ possesses a $b/a$ axial ratio of ~2.5. In addition, we extend the observational arc of 2020 CD$_3$ to 37 days between 2020 February 15 UTC and 2020 March 23 UTC. From the improved orbital solution for 2020 CD$_3$, we estimate its likely duration of its capture to be ~2 y, and we measure the non-gravitation perturbation on its orbit due to radiation pressure with an area-to-mass ratio of 6.9$\pm$2.4$\times$10$^{-4}$ m$^2$/kg implying a density of 2.3$\pm$0.8 g/cm$^3$, broadly compatible with the densities of other meter-scale asteroids and Lunar rock. We searched for pre-discovery detections of 2020 CD$_3$ in the ZTF archive as far back as 2018 October, but were unable to locate any positive detections.

preprint2020arXiv

Characterization of the Nucleus, Morphology and Activity of Interstellar Comet 2I/Borisov by Optical and Near-Infrared GROWTH, Apache Point, IRTF, ZTF and Keck Observations

We present visible and near-infrared photometric and spectroscopic observations of interstellar object 2I/Borisov taken from 2019 September 10 to 2019 November 29 using the GROWTH, the APO ARC 3.5 m and the NASA/IRTF 3.0 m combined with post and pre-discovery observations of 2I obtained by ZTF from 2019 March 17 to 2019 May 5. Comparison with imaging of distant Solar System comets shows an object very similar to mildly active Solar System comets with an out-gassing rate of $\sim$10$^{27}$ mol/sec. The photometry, taken in filters spanning the visible and NIR range shows a gradual brightening trend of $\sim0.03$ mags/day since 2019 September 10 UTC for a reddish object becoming neutral in the NIR. The lightcurve from recent and pre-discovery data reveals a brightness trend suggesting the recent onset of significant H$_2$O sublimation with the comet being active with super volatiles such as CO at heliocentric distances $>$6 au consistent with its extended morphology. Using the advanced capability to significantly reduce the scattered light from the coma enabled by high-resolution NIR images from Keck adaptive optics taken on 2019 October 04, we estimate a diameter of 2I's nucleus of $\lesssim$1.4 km. We use the size estimates of 1I/'Oumuamua and 2I/Borisov to roughly estimate the slope of the ISO size-distribution resulting in a slope of $\sim$3.4$\pm$1.2, similar to Solar System comets and bodies produced from collisional equilibrium.

preprint2020arXiv

Kilonova Luminosity Function Constraints based on Zwicky Transient Facility Searches for 13 Neutron Star Mergers

We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgo's third observing run. We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration. The GW triggers had a median localization of 4480 deg^2, median distance of 267 Mpc and false alarm rates ranging from 1.5 to 1e-25 per yr. The ZTF coverage had a median enclosed probability of 39%, median depth of 20.8mag, and median response time of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UVOIR photometric points, 64 OIR spectra, and 3 radio. We find no promising kilonova (radioactivity-powered counterpart) and we convert the upper limits to constrain the underlying kilonova luminosity function. Assuming that all kilonovae are at least as luminous as GW170817 at discovery (-16.1mag), we calculate our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than -16.6mag (extrapolated peak magnitude of GW170817) and fade at 1 mag/day (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations, the joint probability of zero detections, assuming all kilonovae are brighter than -16.6mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, <57% (<89%) of putative kilonovae could be brighter than -16.6mag assuming flat (fading) evolution, at 90% confidence. If we further account for the online terrestrial probability for each GW trigger, we find that <68% of putative kilonovae could be brighter than -16.6mag. Comparing to model grids, we find that some kilonovae must have Mej < 0.03 Msun or Xlan>1e-4 or phi>30deg to be consistent with our limits. (Abridged)

preprint2019arXiv

GROWTH on S190425z: Searching thousands of square degrees to identify an optical or infrared counterpart to a binary neutron star merger with the Zwicky Transient Facility and Palomar Gattini IR

The third observing run by LVC has brought the discovery of many compact binary coalescences. Following the detection of the first binary neutron star merger in this run (LIGO/Virgo S190425z), we performed a dedicated follow-up campaign with the Zwicky Transient Facility (ZTF) and Palomar Gattini-IR telescopes. The initial skymap of this single-detector gravitational wave (GW) trigger spanned most of the sky observable from Palomar Observatory. Covering 8000 deg$^2$ of the initial skymap over the next two nights, corresponding to 46\% integrated probability, ZTF system achieved a depth of $\approx$\,21 $m_\textrm{AB}$ in $g$- and $r$-bands. Palomar Gattini-IR covered 2200 square degrees in $J$-band to a depth of 15.5\,mag, including 32\% integrated probability based on the initial sky map. The revised skymap issued the following day reduced these numbers to 21\% for the Zwicky Transient Facility and 19\% for Palomar Gattini-IR. We narrowed 338,646 ZTF transient "alerts" over the first two nights of observations to 15 candidate counterparts. Two candidates, ZTF19aarykkb and ZTF19aarzaod, were particularly compelling given that their location, distance, and age were consistent with the GW event, and their early optical lightcurves were photometrically consistent with that of kilonovae. These two candidates were spectroscopically classified as young core-collapse supernovae. The remaining candidates were ruled-out as supernovae. Palomar Gattini-IR did not identify any viable candidates with multiple detections only after merger time. We demonstrate that even with single-detector GW events localized to thousands of square degrees, systematic kilonova discovery is feasible.

preprint2019arXiv

Prompt emission polarimetry of Gamma Ray Bursts with ASTROSAT CZT-Imager

X-ray and Gamma-ray polarization measurements of the prompt emission of Gamma-ray bursts (GRBs) are believed to be extremely important for testing various models of GRBs. So far, the available measurements of hard X-ray polarization of GRB prompt emission have not significantly constrained the GRB models, particularly because of the difficulty of measuring polarization in these bands. The CZT Imager (CZTI) onboard {\em AstroSat} is primarily an X-ray spectroscopic instrument that also works as a wide angle GRB monitor due to the transparency of its support structure above 100 keV. It also has experimentally verified polarization measurement capability in the 100 $-$ 300 keV energy range and thus provides a unique opportunity to attempt spectro-polarimetric studies of GRBs. Here we present the polarization data for the brightest 11 GRBs detected by CZTI during its first year of operation. Among these, 5 GRBs show polarization signatures with $\gtrapprox$3$σ$, and 1 GRB shows $\>$2$σ$ detection significance. We place upper limits for the remaining 5 GRBs. We provide details of the various tests performed to validate our polarization measurements. While it is difficult yet to discriminate between various emission models with the current sample alone, the large number of polarization measurements CZTI expects to gather in its minimum lifetime of five years should help to significantly improve our understanding of the prompt emission.

preprint2015arXiv

Broadband X-ray Properties of the Gamma-ray Binary 1FGL J1018.6-5856

We report on NuSTAR, XMM-Newton and Swift observations of the gamma-ray binary 1FGL J1018.6-5856. We measure the orbital period to be 16.544+/-0.008 days using Swift data spanning 1900 days. The orbital period is different from the 2011 gamma-ray measurement which was used in the previous X-ray study of An et al. (2013) using ~400 days of Swift data, but is consistent with a new gamma-ray solution reported in 2014. The light curve folded on the new period is qualitatively similar to that reported previously, having a spike at phase 0 and broad sinusoidal modulation. The X-ray flux enhancement at phase 0 occurs more regularly in time than was previously suggested. A spiky structure at this phase seems to be a persistent feature, although there is some variability. Furthermore, we find that the source flux clearly correlates with the spectral hardness throughout all orbital phases, and that the broadband X-ray spectra measured with NuSTAR, XMM-Newton, and Swift are well fit with an unbroken power-law model. This spectrum suggests that the system may not be accretion-powered.

preprint2015arXiv

Calibration of the NuSTAR High Energy Focusing X-ray Telescope

We present the calibration of the \textit{Nuclear Spectroscopic Telescope Array} (\nustar) X-ray satellite. We used the Crab as the primary effective area calibrator and constructed a piece-wise linear spline function to modify the vignetting response. The achieved residuals for all off-axis angles and energies, compared to the assumed spectrum, are typically better than $\pm 2$\% up to 40\,keV and 5--10\,\% above due to limited counting statistics. An empirical adjustment to the theoretical 2D point spread function (PSF) was found using several strong point sources, and no increase of the PSF half power diameter (HPD) has been observed since the beginning of the mission. We report on the detector gain calibration, good to 60\,eV for all grades, and discuss the timing capabilities of the observatory, which has an absolute timing of $\pm$ 3\,ms. Finally we present cross-calibration results from two campaigns between all the major concurrent X-ray observatories (\textit{Chandra}, \textit{Swift}, \textit{Suzaku} and \textit{XMM-Newton}), conducted in 2012 and 2013 on the sources 3C\,273 and PKS\,2155-304, and show that the differences in measured flux is within $\sim$10\% for all instruments with respect to \nustar.

preprint2015arXiv

NuSTAR Observations of X-Ray Binaries

As of 2014 August, the Nuclear Spectroscopic Telescope Array (NuSTAR) had observed ~30 X-ray binaries either as part of the planned program, as targets of opportunity, or for instrument calibration. The main science goals for the observations include probing the inner part of the accretion disk and constraining black hole spins via reflection components, providing the first observations of hard X-ray emission from quiescent Low Mass X-ray Binaries (LMXBs), measuring cyclotron lines from accreting pulsars, and studying type I X-ray bursts from neutron stars. Here, we describe the science objectives in more depth and give an overview of the NuSTAR observations that have been carried out to achieve the objectives. These include observation of four "IGR" High Mass X-ray Binaries (HMXBs) discovered by INTEGRAL. We also summarize the results that have been obtained and their implications. Among the IGR HMXBs, we focus on the discovery of a cyclotron line in the spectrum of IGR J17544-2619.

preprint2015arXiv

The Needle in the 100 deg2 Haystack: Uncovering Afterglows of Fermi GRBs with the Palomar Transient Factory

The Fermi Gamma-ray Space Telescope has greatly expanded the number and energy window of observations of gamma-ray bursts (GRBs). However, the coarse localizations of tens to a hundred square degrees provided by the Fermi GRB Monitor instrument have posed a formidable obstacle to locating the bursts' host galaxies, measuring their redshifts, and tracking their panchromatic afterglows. We have built a target-of-opportunity mode for the intermediate Palomar Transient Factory in order to perform targeted searches for Fermi afterglows. Here, we present the results of one year of this program: 8 afterglow discoveries out of 35 searches. Two of the bursts with detected afterglows (GRBs 130702A and 140606B) were at low redshift (z=0.145 and 0.384 respectively) and had spectroscopically confirmed broad-line Type Ic supernovae. We present our broadband follow-up including spectroscopy as well as X-ray, UV, optical, millimeter, and radio observations. We study possible selection effects in the context of the total Fermi and Swift GRB samples. We identify one new outlier on the Amati relation. We find that two bursts are consistent with a mildly relativistic shock breaking out from the progenitor star, rather than the ultra-relativistic internal shock mechanism that powers standard cosmological bursts. Finally, in the context of the Zwicky Transient Facility, we discuss how we will continue to expand this effort to find optical counterparts of binary neutron star mergers that may soon be detected by Advanced LIGO and Virgo.

preprint2014arXiv

NuSTAR Detection Of A Cyclotron Line In The Supergiant Fast X-ray Transient IGR J17544-2619

We present NuSTAR spectral and timing studies of the Supergiant Fast X-ray Transient (SFXT) IGR J17544-2619. The spectrum is well-described by a ~1 keV blackbody and a hard continuum component, as expected from an accreting X-ray pulsar. We detect a cyclotron line at 17 keV, confirming that the compact object in IGR J17544-2619 is indeed a neutron star. This is the first measurement of the magnetic field in a SFXT. The inferred magnetic field strength, B = (1.45 +/- 0.03) * 10^12 G * (1+z) is typical of neutron stars in X-ray binaries, and rules out a magnetar nature for the compact object. We do not find any significant pulsations in the source on time scales of 1-2000 s.

preprint2014arXiv

NuSTAR results and future plans for magnetar and rotation-powered pulsar observations

The Nuclear Spectroscopic Telescope Array (NuSTAR) is the first focusing hard X-ray mission in orbit and operates in the 3-79 keV range. NuSTAR's sensitivity is roughly two orders of magnitude better than previous missions in this energy band thanks to its superb angular resolution. Since its launch in 2012 June, NuSTAR has performed excellently and observed many interesting sources including four magnetars, two rotation-powered pulsars and the cataclysmic variable AE Aquarii. NuSTAR also discovered 3.76-s pulsations from the transient source SGR J1745-29 recently found by Swift very close to the Galactic Center, clearly identifying the source as a transient magnetar. For magnetar 1E 1841-045, we show that the spectrum is well fit by an absorbed blackbody plus broken power-law model with a hard power-law photon index of ~1.3. This is consistent with previous results by INTEGRAL and RXTE. We also find an interesting double-peaked pulse profile in the 25-35 keV band. For AE Aquarii, we show that the spectrum can be described by a multi-temperature thermal model or a thermal plus non-thermal model; a multi-temperature thermal model without a non-thermal component cannot be ruled out. Furthermore, we do not see a spiky pulse profile in the hard X-ray band, as previously reported based on Suzaku observations. For other magnetars and rotation-powered pulsars observed with NuSTAR, data analysis results will be soon available.

preprint2014arXiv

X-ray Spectral Components in the Afterglow of GRB 130925A

We have identified spectral features in the late-time X-ray afterglow of the unusually long, slow-decaying GRB 130925A using NuSTAR, Swift-XRT, and Chandra. A spectral component in addition to an absorbed power-law is required at $>4σ$ significance, and its spectral shape varies between two observation epochs at $2\times10^5$ and $10^6$ seconds after the burst. Several models can fit this additional component, each with very different physical implications. A broad, resolved Gaussian absorption feature of several keV width improves the fit, but it is poorly constrained in the second epoch. An additive black body or second power-law component provide better fits. Both are challenging to interpret: the blackbody radius is near the scale of a compact remnant ($10^8$ cm), while the second powerlaw component requires an unobserved high-energy cutoff in order to be consistent with the non-detection by Fermi-LAT.

preprint2013arXiv

The Nuclear Spectroscopic Telescope Array (NuSTAR) Mission

The Nuclear Spectroscopic Telescope Array (NuSTAR) mission, launched on 13 June 2012, is the first focusing high-energy X-ray telescope in orbit. NuSTAR operates in the band from 3 -- 79 keV, extending the sensitivity of focusing far beyond the ~10 keV high-energy cutoff achieved by all previous X-ray satellites. The inherently low-background associated with concentrating the X-ray light enables NuSTAR to probe the hard X-ray sky with a more than one-hundred-fold improvement in sensitivity over the collimated or coded-mask instruments that have operated in this bandpass. Using its unprecedented combination of sensitivity, spatial and spectral resolution, NuSTAR will pursue five primary scientific objectives, and will also undertake a broad program of targeted observations. The observatory consists of two co-aligned grazing-incidence X-ray telescopes pointed at celestial targets by a three-axis stabilized spacecraft. Deployed into a 600 km, near-circular, 6degree inclination orbit, the Observatory has now completed commissioning, and is performing consistent with pre-launch expectations. NuSTAR is now executing its primary science mission, and with an expected orbit lifetime of ten years, we anticipate proposing a guest investigator program, to begin in Fall 2014.

preprint2012arXiv

Constraints on the Compact Object Mass in the Eclipsing HMXB XMMU J013236.7+303228 in M33

We present optical spectroscopic measurements of the eclipsing High Mass X-ray Binary XMMU J013236.7+303228 in M33. Based on spectra taken at multiple epochs of the 1.73d binary orbital period we determine physical as well as orbital parameters for the donor star. We find the donor to be a B1.5IV sub-giant with effective temperature T=22,000-23,000 K. From the luminosity, temperature and known distance to M33 we derive a radius of R = 8.9 \pm 0.5 R_sun. From the radial--velocity measurements, we determine a velocity semi-amplitude of K_opt = 63 \pm 12 km/sec. Using the physical properties of the B-star determined from the optical spectrum, we estimate the star's mass to be M_opt = 11 \pm 1 M_sun. Based on the X-ray spectrum, the compact companion is likely a neutron star, although no pulsations have yet been detected. Using the spectroscopically derived B-star mass we find the neutron star companion mass to be M_X = 2.0 \pm 0.4 M_sun, consistent with the neutron star mass in the HMXB Vela X-1, but heavier than the canonical value of 1.4 M_sun found for many millisecond pulsars. We attempt to use as an additional constraint that the B star radius inferred from temperature, flux, and distance, should equate the Roche radius, since the system accretes by Roche lobe overflow. This leads to substantially larger masses, but from trying to apply the technique to known systems, we find that the masses are consistently overestimated. Attempting to account for that in our uncertainties, we derive M_X = 2.2^{+0.8}_{-0.6} M_sun and M_opt =13 \pm 4 M_sun. We conclude that precise constraints require detailed modeling of the shape of the Roche surface.