Reducing LRIS longslit spectra in IRAF

The word IRAF scares me to date, though I am slowly getting used to using it. Here I am putting together my notes for reducing LRIS longslit spectra in IRAF. I will try to be general, but only to the extent that I expect my usage to vary. The document will often contain things which are pertinent to my current folders setup, gratings I use, etc. The scripts referred to in this document are available on request. I might eventually upload a tarball containing all the codes

Calcium-rich gap transients in the remote outskirts of galaxies

From the first two seasons of the Palomar Transient Factory, we identify three peculiar transients (PTF09dav, PTF10iuv, PTF11bij) with five distinguishing characteristics: peak luminosity in the gap between novae and supernovae (M_R ≈ - 15.5 to -16.5), rapid photometric evolution (t_(rise) ≈12-15 days), large photospheric velocities (≈6000 to 11000 km s^(-1)), early spectroscopic evolution into nebular phase (≈1 to 3 months) and peculiar nebular spectra dominated by Calcium. We also culled the extensive decade-long Lick Observatory Supernova Search database and identified an additional member of this group, SN 2007ke. Our choice of photometric and spectroscopic properties was motivated by SN 2005E (Perets et al. 2010). To our surprise, as in the case of SN 2005E, all four members of this group are also clearly offset from the bulk of their host galaxy. Given the well-sampled early and late-time light curves, we derive ejecta masses in the range of 0.4--0.7 M_⊙. Spectroscopically, we find that there may be a diversity in the photospheric phase, but the commonality is in the unusual nebular spectra. Our extensive follow-up observations rule out standard thermonuclear and standard core-collapse explosions for this class of "Calcium-rich gap" transients. If the progenitor is a white dwarf, we are likely seeing a detonation of the white dwarf core and perhaps, even shock-front interaction with a previously ejected nova shell. In the less likely scenario of a massive star progenitor, a very non-standard channel specific to a low-metallicity environment needs to be invoked (e.g., ejecta fallback leading to black hole formation). Detection (or lack thereof) of a faint underlying host (dwarf galaxy, cluster) will provide a crucial and decisive diagnostic to choose between these alternatives

Neutron Stars and NuSTAR: A Systematic Survey of Neutron Star Masses in High Mass X-ray Binaries & Characterization of CdZnTe Detectors for NuSTAR

My thesis centers around the study of neutron stars, especially those in massive binary systems. To this end, it has two distinct components: the observational study of neutron stars in massive binaries with a goal of measuring neutron star masses and participation in NuSTAR, the first imaging hard X-ray mission, one that is extremely well suited to the study of massive binaries and compact objects in our Galaxy. The Nuclear Spectroscopic Telescope Array (NuSTAR) is a NASA Small Explorer mission that will carry the first focusing high energy X-ray telescope to orbit. NuSTAR has an order-of-magnitude better angular resolution and has two orders of magnitude higher sensitivity than any currently orbiting hard X-ray telescope. I worked to develop, calibrate, and test CdZnTe detectors for NuSTAR. I describe the CdZnTe detectors in comprehensive detail here — from readout procedures to data analysis. Detailed calibration of detectors is necessary for analyzing astrophysical source data obtained by the NuSTAR. I discuss the design and implementation of an automated setup for calibrating flight detectors, followed by calibration procedures and results. Neutron stars are an excellent probe of fundamental physics. The maximum mass of a neutron star can put stringent constraints on the equation of state of matter at extreme pressures and densities. From an astrophysical perspective, there are several open questions in our understanding of neutron stars. What are the birth masses of neutron stars? How do they change in binary evolution? Are there multiple mechanisms for the formation of neutron stars? Measuring masses of neutron stars helps answer these questions. Neutron stars in high-mass X-ray binaries have masses close to their birth mass, providing an opportunity to disentangle the role of "nature" and "nurture" in the observed mass distributions. In 2006, masses had been measured for only six such objects, but this small sample showed the greatest diversity in masses among all classes of neutron star binaries. Intrigued by this diversity — which points to diverse birth masses — we undertook a systematic survey to measure the masses of neutron stars in nine high-mass X-ray binaries. In this thesis, I present results from this ongoing project. While neutron stars formed the primary focus of my work, I also explored other topics in compact objects. Appendix A describes the discovery and complete characterization of a 1RXS J173006.4+033813, a polar cataclysmic variable. Appendix B describes the discovery of a diamond planet orbiting a millisecond pulsar, and our search for its optical counterpart.</p

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$\sigma$, and 1 GRB shows $\>$2$\sigma$ 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.Comment: Accepted for Publication in ApJ ; a figure has been update

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.Comment: 11 pages emulateapj, 9 figures, ApJ accepte

The polar Catalysmic Variable 1RXS J173006.4+033813

We report the discovery of 1RXS J173006.4+033813, a polar cataclysmic variable with a period of 120.21 min. The white dwarf primary has a magnetic field of B = 42+6-5 MG, and the secondary is a M3 dwarf. The system shows highly symmetric double peaked photometric modulation in the active state as well as in quiescence. These arise from a combination of cyclotron beaming and ellipsoidal modulation. The projected orbital velocity of the secondary is K2 = 390+-4 km/s. We place an upper limit of 830+-65 pc on the distance.Comment: ApJ Accepted. 12 Pages, 13 Figures, 6 table

Tracing the Orphan Stream to 55 kpc with RR Lyrae Stars

We report positions, velocities and metallicities of 50 ab-type RR Lyrae (RRab) stars observed in the vicinity of the Orphan stellar stream. Using about 30 RRab stars classified as being likely members of the Orphan stream, we study the metallicity and the spatial extent of the stream. We find that RRab stars in the Orphan stream have a wide range of metallicities, from -1.5 dex to -2.7 dex. The average metallicity of the stream is -2.1 dex, identical to the value obtained by Newberg et al. (2010) using blue horizontal branch stars. We find that the most distant parts of the stream (40-50 kpc from the Sun) are about 0.3 dex more metal-poor than the closer parts (within ~30 kpc), suggesting a possible metallicity gradient along the stream's length. We have extended the previous studies and have mapped the stream up to 55 kpc from the Sun. Even after a careful search, we did not identify any more distant RRab stars that could plausibly be members of the Orphan stream. If confirmed with other tracers, this result would indicate a detection of the end of the leading arm of the stream. We have compared the distances of Orphan stream RRab stars with the best-fit orbits obtained by Newberg et al. (2010). We find that model 6 of Newberg et al. (2010) cannot explain the distances of the most remote Orphan stream RRab stars, and conclude that the best fit to distances of Orphan stream RRab stars and to the local circular velocity is provided by potentials where the total mass of the Galaxy within 60 kpc is M_{60}~2.7x10^{11} Msun, or about 60% of the mass found by previous studies. More extensive modelling that would consider non-spherical potentials and the possibility of misalignment between the stream and the orbit, is highly encouraged.Comment: Submitted to ApJ, 15 pages in emulateapj format, three tables in machine-readable format (download "Source" from "Other formats"

NuSTAR Discovery of a Cyclotron Line in the Be/X-ray Binary RX J0520.5-6932 During Outburst

We present spectral and timing analysis of NuSTAR observations of RX J0520.5$-$6932 in the 3-79 keV band collected during its outburst in January 2014. The target was observed on two epochs and we report the detection of a cyclotron resonant scattering feature with central energies of $E_\mathrm{CRSF} = 31.3_{-0.7}^{+0.8}$ keV and $31.5_{-0.6}^{+0.7}$ keV during the two observations, respectively, corresponding to a magnetic field of $B \approx 2 \times10^{12}$ G. The 3-79 keV luminosity of the system during the two epochs assuming a nominal distance of 50 kpc was $3.667\pm0.007\times 10^{38}\,\mathrm{erg\,s^{-1}}$ and $3.983\pm0.007\times10^{38}\,\mathrm{erg\,s^{-1}}$. Both values are much higher than the critical luminosity of $\approx1.5\times10^{37}\,\mathrm{erg\,s^{-1}}$ above which a radiation dominated shock front may be expected. This adds a new object to the sparse set of three systems that have a cyclotron line observed at luminosities in excess of $10^{38}\,\mathrm{erg\,s^{-1}}$. A broad ($\sigma\approx0.45$ keV) Fe emission line is observed in the spectrum at a central energy of $6.58_{-0.05}^{+0.05}$ keV in both epochs. The pulse profile of the pulsar was observed to be highly asymmetric with a sharply rising and slowly falling profile of the primary peak. We also observed minor variations in the cyclotron line energy and width as a function of the rotation phase.% As in observations of other cyclotron absorption line sources, there is a small ($\Delta\phi\lesssim0.1$) phase difference between the peak of the cyclotron energy variation and the peak of the flux variation.Comment: 9 pages, 4 figures, 6 tables. Accepted to Ap

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