Source author record

Craig O. Heinke

Craig O. Heinke 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

45works
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

45 published item(s)

preprint2022arXiv

1RXH J082623.6-505741: a new long-period cataclysmic variable with an evolved donor and a low mass transfer rate

We report the discovery of 1RXH J082623.6-505741, a 10.4 hr orbital period compact binary. Modeling extensive optical photometry and spectroscopy reveals a $\sim 0.4 M_{\odot}$ K-type secondary transferring mass through a low-state accretion disk to a non-magnetic $\sim 0.8 M_{\odot}$ white dwarf. The secondary is overluminous for its mass and dominates the optical spectra at all epochs, and must be evolved to fill its Roche Lobe at this orbital period. The X-ray luminosity $L_X \sim 1$-$2 \times 10^{32}$ erg s$^{-1}$ derived from both new XMM-Newton and archival observations, although high compared to most CVs, still only requires a modest accretion rate onto the white dwarf of $\dot{M} \sim 3 \times 10^{-11}$ to $3 \times 10^{-10} M_{\odot}$ yr$^{-1}$, lower than expected for a cataclysmic variable with an evolved secondary. No dwarf nova outbursts have yet been observed from the system, consistent with the low derived mass transfer rate. Several other cataclysmic variables with similar orbital periods also show unexpectedly low mass transfer rates, even though selection effects disfavor the discovery of binaries with these properties. This suggests the abundance and evolutionary state of long-period, low mass transfer rate cataclysmic variables is worthy of additional attention.

preprint2022arXiv

A Census of X-ray Millisecond Pulsars in Globular Clusters

We present a comprehensive census of X-ray millisecond pulsars (MSPs) in 29 Galactic globular clusters (GCs), including 68 MSPs with confirmed X-ray luminosities and 107 MSPs with X-ray upper limits. We compile previous X-ray studies of GC MSPs, and add new analyses of six MSPs (PSRs J1326$-$4728A, J1326$-$4728B, J1518$+$0204C, J1717$+$4308A, J1737$-$0314A, and J1807$-$2459A) discovered in five GCs. Their X-ray spectra are well described by a single blackbody model, a single power-law model, or a combination of them, with X-ray luminosities ranging from 1.9$\times$10$^{30}$ erg s$^{-1}$ to 8.3$\times$10$^{31}$ erg s$^{-1}$. We find that most detected X-ray MSPs have luminosities between $\sim 10^{30}$ erg s$^{-1}$ to $3 \times 10^{31}$ erg s$^{-1}$. Redback pulsars are a relatively bright MSP population with X-ray luminosities of $\sim2\times10^{31}$--$3\times10^{32}$ erg s$^{-1}$. Black widows show a bi-modal distribution in X-ray luminosities, with eclipsing black widows between $\sim 7\times10^{30}$ and $2\times10^{31}$ erg s$^{-1}$, while the two confirmed non-eclipsing black widows are much fainter, with $L_X$ of $1.5-3\times10^{30}$ erg s$^{-1}$, suggesting an intrinsic difference in the populations. We estimate the total number of MSPs in 36 GCs by considering the correlation between the number of MSPs and stellar encounter rate in GCs, and suggest that between 600--1500 MSPs exist in these 36 GCs. Finally, we estimate the number of X-ray detectable MSPs in the Galactic bulge, finding that 1--86 MSPs with $L_X > 10^{33}$ erg s$^{-1}$, and 20--900 MSPs with $L_X > 10^{32}$ erg s$^{-1}$, should be detectable there.

preprint2022arXiv

Investigating the low-flux states in six Intermediate Polars

We present optical photometry of six intermediate polars that exhibit transitions to a low-flux state. For four of these systems, DW Cnc, V515 And, V1223 Sgr and RX J2133.7+5107, we are able to perform timing analysis in and out of the low states. We find that, for DW Cnc and V515 And, the dominant periodicities in the light curves change as the flux decreases, indicating a change in the sources' accretion properties as they transition to the low state. For V1223 Sgr we find that the variability is almost completely quenched at the lowest flux, but do not find evidence for a changing accretion geometry. For RX J2133.7+5107, the temporal properties do not change in the low state, but we do see a period of enhanced accretion that is coincident with increased variability on the beat frequency, which we do not associate with a change in the accretion mechanisms in the system.

preprint2022arXiv

Millisecond Pulsar Kicks Cause Difficulties in Explaining the Galactic Center Gamma-Ray Excess

The unexplained excess gamma-ray emission from the Milky Way's Galactic Center has puzzled astronomers for nearly a decade. Two theories strive to explain the origin of this excess: self-annihilating dark matter particles or an unresolved population of radio millisecond pulsars. We examine the plausibility of a pulsar origin for the GeV excess using N-body simulations. We simulated millisecond pulsars in a realistic dynamical environment: (i) pulsars were born from the known stellar mass components of our Galaxy; (ii) pulsars were given natal velocity kicks as empirically observed from two different studies (or, for comparison, without kicks); (iii) pulsars were evolved in a Galactic gravitational potential consistent with observations. Multiple populations of pulsars (with different velocity kicks) were simulated over 1 Gyr. With final spatial distributions of pulsars, we constructed synthetic gamma-ray surface brightness profiles. From comparisons with published Fermi-LAT surface brightness profiles, our pulsar simulations cannot reproduce the concentrated emission in the central degrees of the Bulge, though models without natal velocity kicks approach the data. We considered additive combinations of our (primordial MSP) simulations with models where pulsars are deposited from destroyed globular clusters in the Bulge, and a simple model for pulsars produced in the nuclear star cluster. We can reasonably reproduce the measured central gamma-ray surface brightness distribution of Horiuchi and collaborators using several combinations of these models, but we cannot reproduce the measured distribution of Di Mauro with any combination of models. Our fits provide constraints on potential pathways to explain the gamma-ray excess using MSPs.

preprint2022arXiv

The MAVERIC survey: A catalogue of radio sources in southern globular clusters from the Australia Telescope Compact Array

Radio continuum observations offer a new window on compact objects in globular clusters compared to typical X-ray or optical studies. As part of the MAVERIC survey, we have used the Australia Telescope Compact Array to carry out a deep (median central noise level of approximately 4 microJy per beam) radio continuum survey of 26 southern globular clusters at central frequencies of 5.5 and 9.0 GHz. This paper presents a catalogue of 1285 radio continuum sources in the fields of these 26 clusters. Considering the surface density of background sources, we find significant evidence for a population of radio sources in seven of the 26 clusters, and also identify at least 11 previously known compact objects (6 pulsars and 5 X-ray binaries). While the overall density of radio continuum sources with 7.25-GHz flux densities greater than about 20 microJy in typical globular clusters is relatively low, the survey has already led to the discovery of several exciting compact binaries, including a candidate ultracompact black hole X-ray binary in 47 Tuc. Many of the unclassified radio sources near the centres of the clusters are likely to be true cluster sources, and multi-wavelength follow-up will be necessary to classify these objects and better understand the demographics of accreting compact binaries in globular clusters.

preprint2022arXiv

The MAVERIC Survey: The first radio and X-ray limits on the detached black holes in NGC 3201

The Galactic globular cluster NGC 3201 is the first Galactic globular cluster observed to host dynamically-confirmed stellar-mass black holes, containing two confirmed and one candidate black hole. This result indicates that globular clusters can retain black holes, which has important implications for globular cluster evolution. NGC 3201 has been observed as part of the MAVERIC survey of Galactic globular clusters. We use these data to confirm that there is no radio or X-ray detection of the three black holes, and present the first radio and X-ray limits on these sources. These limits indicate that any accretion present is at an extremely low rate and may be extremely inefficient. In particular, for the system ACS ID #21859, by assuming the system is tidally locked and any accretion is through the capture of the companion's winds, we constrain the radiative efficiency of any accretion to $\leq1.5\times10^{-5}$. We also combine the radio and X-ray source catalogues from the MAVERIC survey with the existing MUSE spectroscopic surveys and the HUGS catalogue of NGC 3201 to provide a catalogue of 42 multiwavelength sources in this cluster. We identify a new red straggler source with X-ray emission, and investigate the multiwavelength properties of the sub-subgiant population in the cluster.

preprint2021arXiv

Chandra and HST Studies of Six Millisecond Pulsars in the Globular Cluster M13

We analyse 55 ks of Chandra X-ray observations of the Galactic globular cluster M13. Using the latest radio timing positions of six known millisecond pulsars (MSPs) in M13 from Wang et al. (2020), we detect confident X-ray counterparts to five of the six MSPs at X-ray luminosities of $L_X$(0.3-8 keV)$\sim 3 \times 10^{30} - 10^{31}~{\rm erg~s^{-1}}$, including the newly discovered PSR J1641+3627F. There are limited X-ray counts at the position of PSR J1641+3627A, for which we obtain an upper limit $L_X<1.3 \times 10^{30}~{\rm erg~s^{-1}}$. We analyse X-ray spectra of all six MSPs, which are well-described by either a single blackbody or a single power-law model. We also incorporate optical/UV imaging observations from the Hubble Space Telescope (HST) and find optical counterparts to PSR J1641+3627D and J1641+3627F. Our colour-magnitude diagrams indicate the latter contains a white dwarf, consistent with the properties suggested by radio timing observations. The counterpart to J1641+3627D is only visible in the V band; however, we argue that the companion to J1641+3627D is also a white dwarf, since we see a blackbody-like X-ray spectrum, while MSPs with nondegenerate companions generally show non-thermal X-rays from shocks between the pulsar and companion winds. Our work increases the sample of known X-ray and optical counterparts of MSPs in globular clusters.

preprint2021arXiv

Dynamical Modeling of CXOGBS J175553.2-281633: A 10 Hour Long Orbital Period Cataclysmic Variable

We present modeling of the long-term optical light curve and radial velocity curve of the binary stellar system CXOGBS J175553.2-281633, first detected in X-rays in the \textit{Chandra} Galactic Bulge Survey. We analyzed 7 years of optical I-band photometry from OGLE and found long-term variations from year to year. These long-term variations can most likely be explained with by either variations in the luminosity of the accretion disk or a spotted secondary star. The phased light curve has a sinusoidal shape, which we interpret as being due to ellipsoidal modulations. We improve the orbital period to be $P = 10.34488 \pm 0.00006$ h with a time of inferior conjunction of the secondary star $T_0 = {\rm HJD\ } 2455260.8204 \pm 0.0008$. Moreover, we collected 37 spectra over 6 non-consecutive nights. The spectra show evidence for an evolved K7 secondary donor star, from which we obtain a semi-amplitude for the radial velocity curve of $K_2 = 161 \pm 6 $ km s$^{-1}$. Using the light curve synthesis code {\tt XRbinary}, we derive the most likely orbital inclination for the binary of $i = 63.0\pm0.7$ deg, a primary mass of $M_1 = 0.83 \pm 0.06$ M$_\odot$, consistent with a white dwarf accretor, and a secondary donor mass of $M_2 = 0.65 \pm 0.07$ M$_\odot$, consistent with the spectral classification. Therefore, we identify the source as a long orbital period cataclysmic variable star.

preprint2020arXiv

Spectral analysis of the quiescent low-mass X-ray binary in the globular cluster M30

We present a recent Chandra observation of the quiescent low-mass X-ray binary containing a neutron star, located in the globular cluster M30. We fit the thermal emission from the neutron star to extract its mass and radius. We find no evidence of flux variability between the two observations taken in 2001 and 2017, nor between individual 2017 observations, so we analyse them together to increase the signal to noise. We perform simultaneous spectral fits using standard light-element composition atmosphere models (hydrogen or helium), including absorption by the interstellar medium, correction for pile-up of X-ray photons on the detector, and a power-law for count excesses at high photon energy. Using a Markov-chain Monte Carlo approach, we extract mass and radius credible intervals for both chemical compositions of the atmosphere: $R_{\textrm{NS}}=7.94^{+0.76}_{-1.21}$ km and $M_{\textrm{NS}}<1.19$ M$_{\odot}$ assuming pure hydrogen, and $R_{\textrm{NS}}=10.50^{+2.88}_{-2.03}$ km and $M_{\textrm{NS}}<1.78$ M$_{\odot}$ for helium, where the uncertainties represent the 90% credible regions. For H, the small radius is difficult to reconcile with most current nuclear physics models (especially for nucleonic equations of state) and with other measurements of neutron star radii, with recent preferred values generally in the 11-14 km range. Whereas for He, the measured radius is consistent with this range. We discuss possible sources of systematic uncertainty that may result in an underestimation of the radius, identifying the presence of surface temperature inhomogeneities as the most relevant bias. According to this, we conclude that either the atmosphere is composed of He, or it is a H atmosphere with a significant contribution of hot spots to the observed radiation.

preprint2020arXiv

The ELM Survey. VIII. 98 Double White Dwarf Binaries

We present the final sample of 98 detached double white dwarf (WD) binaries found in the Extremely Low Mass (ELM) Survey, a spectroscopic survey targeting <0.3 Msun He-core WDs completed in the Sloan Digital Sky Survey footprint. Over the course of the survey we observed ancillary low mass WD candidates like GD278, which we show is a P=0.19 d double WD binary, as well as candidates that turn out to be field blue straggler/subdwarf A-type stars with luminosities too large to be WDs given their Gaia parallaxes. Here, we define a clean sample of ELM WDs that is complete within our target selection and magnitude range 15<g_0<20 mag. The measurements are consistent with 100% of ELM WDs being 0.0089 < P < 1.5 d double WD binaries, 35% of which belong to the Galactic halo. We infer these are mostly He+CO WD binaries given the measurement constraints. The merger rate of the observed He+CO WD binaries exceeds the formation rate of stable mass transfer AM CVn binaries by a factor of 25, and so the majority of He+CO WD binaries must experience unstable mass transfer and merge. The shortest-period systems like J0651+2844 are signature {\it LISA} verification binaries that can be studied with gravitational waves and light.

preprint2020arXiv

The Flare-Dominated Accretion Mode of a Radio-Bright Candidate Transitional Millisecond Pulsar

We report new simultaneous X-ray and radio continuum observations of 3FGL J0427.9-6704, a candidate member of the enigmatic class of transitional millisecond pulsars. These XMM-Newton and Australia Telescope Compact Array observations of this nearly edge-on, eclipsing low-mass X-ray binary were taken in the sub-luminous disk state at an X-ray luminosity of ~10^33(d/2.3kpc)^2 erg/s. Unlike the few well-studied transitional millisecond pulsars, which spend most of their disk state in a characteristic high or low accretion mode with occasional flares, 3FGL J0427.9-6704 stayed in the flare mode for the entire X-ray observation of ~20 hours, with the brightest flares reaching ~2*10^34 erg/s. The source continuously exhibited flaring activity on time-scales of ~10-100 sec in both the X-ray and optical/UV. No measurable time delay between the X-ray and optical/UV flares is observed, but the optical/UV flares last longer, and the relative amplitudes of the X-ray and optical/UV flares show a large scatter. The X-ray spectrum can be well-fit with a partially-absorbed power-law (Gamma~1.4-1.5), perhaps due to the edge-on viewing angle. Modestly variable radio continuum emission is present at all epochs, and is not eclipsed by the secondary, consistent with the presence of a steady radio outflow or jet. The simultaneous radio/X-ray luminosity ratio of 3FGL J0427.9-6704 is higher than any known transitional millisecond pulsars and comparable to that of stellar-mass black holes of the same X-ray luminosity, providing additional evidence that some neutron stars can be as radio-loud as black holes.

preprint2020arXiv

The MAVERIC survey: A hidden pulsar and a black hole candidate in ATCA radio imaging of the globular cluster NGC 6397

Using a 16.2 hr radio observation by the Australia Telescope Compact Array (ATCA) and archival Chandra data, we found $>5σ$ radio counterparts to 4 known and 3 new X-ray sources within the half-light radius ($r_\mathrm{h}$) of the Galactic globular cluster NGC 6397. The previously suggested millisecond pulsar (MSP) candidate, U18, is a steep-spectrum ($S_ν\propto ν^α$; $α=-2.0^{+0.4}_{-0.5}$) radio source with a 5.5 GHz flux density of $54.7\pm 4.3~\mathrm{μJy}$. We argue that U18 is most likely a "hidden" MSP that is continuously hidden by plasma shocked at the collision between the winds from the pulsar and companion star. The nondetection of radio pulsations so far is probably the result of enhanced scattering in this shocked wind. On the other hand, we observed 5.5 GHz flux of the known MSP PSR J1740-5340 (U12) to decrease by a factor of $>2.8$ during epochs of 1.4 GHz eclipse, indicating that the radio flux is absorbed in its shocked wind. If U18 is indeed a pulsar whose pulsations are scattered, we note the contrast with U12's flux decrease in eclipse, which argues for two different eclipse mechanisms at the same radio frequency. In addition to U12 and U18, we also found radio associations for 5 other Chandra X-ray sources, four of which are likely background galaxies. The last, U97, which shows strong H$α$ variability, is mysterious; it may be either a quiescent black hole low-mass X-ray binary, or something more unusual.

preprint2019arXiv

LRP2020: The cosmic origin and evolution of the elements

The origin of many elements of the periodic table remains an unsolved problem. While many nucleosynthetic channels are broadly understood, significant uncertainties remain regarding certain groups of elements such as the intermediate and rapid neutron-capture processes, the p-process, or the origin of odd-Z elements in the most metal-poor stars. Canada has a long tradition of leadership in nuclear astrophysics, dating back to the work of Alastair Cameron in the 1950s. Recent faculty hires have further boosted activity in the field, including transient observation and theory, survey science on galactic nucleosynthesis, and nuclear experiments. This white paper contains a brief overview of recent activity in the community, highlighting strengths in each sub-field, and provides recommendations to improve interdisciplinary collaboration. Sustaining Canadian leadership in the next decade will require, on the observational side, access to transient and non-transient surveys like LSST, SKA, or MSE, support for target-of-opportunity observing in current and future Canadian telescopes, and participation in next-generation X-ray telescopes such as ATHENA. State-of-the-art theoretical predictions will require an ambitious succession plan for the Niagara supercomputer to support large parallel jobs. We propose a funding instrument for postdoctoral training that reflects the interdisciplinary nature of nuclear astrophysics research, and the creation of a national collaborative funding program that allows for joint projects and workshop organization.

preprint2016arXiv

A new gamma-ray loud, eclipsing low-mass X-ray binary

We report the discovery of an eclipsing low-mass X-ray binary at the center of the 3FGL error ellipse of the unassociated Fermi/Large Area Telescope gamma-ray source 3FGL J0427.9-6704. Photometry from OGLE and the SMARTS 1.3-m telescope and spectroscopy from the SOAR telescope have allowed us to classify the system as an eclipsing low-mass X-ray binary (P = 8.8 hr) with a main sequence donor and a neutron star accretor. Broad double-peaked H and He emission lines suggest the ongoing presence of an accretion disk. Remarkably, the system shows shows separate sets of absorption lines associated with the accretion disk and the secondary, and we use their radial velocities to find evidence for a massive (~ 1.8-1.9 M_sun) neutron star primary. In addition to a total X-ray eclipse of duration ~ 2200 s observed with NuSTAR, the X-ray light curve also shows properties similar to those observed among known transitional millisecond pulsars: short-term variability, a hard power-law spectrum (photon index ~ 1.7), and a comparable 0.5-10 keV luminosity (~ 2.4 x 10^33 erg/s). We find tentative evidence for a partial (~ 60%) gamma-ray eclipse at the same phase as the X-ray eclipse, suggesting the gamma-ray emission may not be confined to the immediate region of the compact object. The favorable inclination of this binary is promising for future efforts to determine the origin of gamma-rays among accreting neutron stars.

preprint2016arXiv

Massive Double White Dwarfs and the AM CVn Birthrate

We present Chandra and Swift X-ray observations of four extremely low-mass (ELM) white dwarfs with massive companions. We place stringent limits on X-ray emission from all four systems, indicating that neutron star companions are extremely unlikely and that the companions are almost certainly white dwarfs. Given the observed orbital periods and radial velocity amplitudes, the total masses of these binaries are greater than 1.02 to 1.39 Msun. The extreme mass ratios between the two components make it unlikely that these binary white dwarfs will merge and explode as Type Ia or underluminous supernovae. Instead, they will likely go through stable mass transfer through an accretion disk and turn into interacting AM CVn. Along with three previously known systems, we identify two of our targets, J0811 and J2132, as systems that will definitely undergo stable mass transfer. In addition, we use the binary white dwarf sample from the ELM Survey to constrain the inspiral rate of systems with extreme mass ratios. This rate, 0.00017/year, is consistent with the AM CVn space density estimated from the Sloan Digital Sky Survey. Hence, stable mass transfer double white dwarf progenitors can account for the entire AM CVn population in the Galaxy.

preprint2016arXiv

Neutron Star Mass-Radius Constraints of the Quiescent Low-mass X-ray Binaries X7 and X5 in the Globular Cluster 47 Tuc

We present Chandra ACIS-S sub-array observations of the quiescent neutron star low-mass X-ray binaries X7 and X5 in the globular cluster 47 Tuc. The large reduction in photon pile-up compared to previous deep exposures enables a substantial improvement in the spectroscopic determination of the neutron star radius and mass of these neutron stars. Modeling the thermal emission from the neutron star surface with a non-magnetized hydrogen atmosphere and accounting for numerous sources of uncertainties, we obtain for the neutron star in X7 a radius of $R=11.1^{+0.8}_{-0.7}$ km for an assumed stellar mass of $M=1.4$ M$_{\odot}$ (68% C.L.). We argue, based on astrophysical grounds, that the presence of a He atmosphere is unlikely for this source. Due to eclipses and variable absorption, the quiescent low-mass X-ray binary X5 provides less stringent constraints, leading to a radius of $R=9.6^{+0.9}_{-1.1}$ km, assuming a hydrogen atmosphere and a mass of $M=1.4$ M$_{\odot}$. When combined with all other existing spectroscopic radius measurements, these measurements strongly favor radii in the 9.9-11.2 km range for a ~1.5 M$_{\odot}$ neutron star and point to a dense matter equation of state that is somewhat softer than the nucleonic ones that are consistent with laboratory experiments at low densities.

preprint2015arXiv

Limits on thermal variations in a dozen quiescent neutron stars over a decade

In quiescent low-mass X-ray binaries (qLMXBs) containing neutron stars, the origin of the thermal X-ray component may be either release of heat from the core of the neutron star, or continuing low-level accretion. In general, heat from the core should be stable on timescales $<10^4$ years, while continuing accretion may produce variations on a range of timescales. While some quiescent neutron stars (e.g. Cen X-4, Aql X-1) have shown variations in their thermal components on a range of timescales, several others, particularly those in globular clusters with no detectable nonthermal hard X-rays (fit with a powerlaw), have shown no measurable variations. Here, we constrain the spectral variations of 12 low mass X-ray binaries in 3 globular clusters over $\sim10$ years. We find no evidence of variations in 10 cases, with limits on temperature variations below 11% for the 7 qLMXBs without powerlaw components, and limits on variations below 20% for 3 other qLMXBs that do show non-thermal emission. However, in 2 qLMXBs showing powerlaw components in their spectra (NGC 6440 CX 1 & Terzan 5 CX 12) we find marginal evidence for a 10% decline in temperature, suggesting the presence of continuing low-level accretion. This work adds to the evidence that the thermal X-ray component in quiescent neutron stars without powerlaw components can be explained by heat deposited in the core during outbursts. Finally, we also investigate the correlation between hydrogen column density (N$_H$) and optical extinction (A$_V$) using our sample and current models of interstellar X-ray absorption, finding $N_H ({\rm cm}^{-2}) = (2.81\pm0.13)\times10^{21} A_V$.

preprint2015arXiv

Tests of the nuclear equation of state and superfluid and superconducting gaps using the Cassiopeia A neutron star

The observed rapid cooling of the Cassiopeia A neutron star can be interpreted as being caused by neutron and proton transitions from normal to superfluid and superconducting states in the stellar core. Here we present two new Chandra ACIS-S Graded observations of this neutron star and measurements of the neutron star mass M and radius R found from consistent fitting of both the X-ray spectra and cooling behavior. This comparison is only possible for individual nuclear equations of state. We test phenomenological superfluid and superconducting gap models which mimic many of the known theoretical models against the cooling behavior. Our best-fit solution to the Cassiopeia A data is one in which the (M,R) = (1.44 Msun,12.6 km) neutron star is built with the BSk21 equation of state, strong proton superconductor and moderate neutron triplet superfluid gap models, and a pure iron envelope or a thin carbon layer on top of an iron envelope, although there are still large observational and theoretical uncertainties.

preprint2014arXiv

On the coexistence of stellar-mass and intermediate-mass black holes in globular clusters

In this paper, we address the question: What is the probability of stellar-mass black hole (BH) binaries co-existing in a globular cluster with an intermediate-mass black hole (IMBH)? Our results suggest that the detection of one or more BH binaries can strongly constrain the presence of an IMBH in most Galactic globular clusters. More specifically, the detection of one or more BH binaries could strongly indicate against the presence of an IMBH more massive than $\gtrsim 10^3$ M$_{\rm \odot}$ in roughly 80\% of the clusters in our sample. To illustrate this, we use a combination of N-body simulations and analytic methods to weigh the rate of formation of BH binaries against their ejection and/or disruption rate via strong gravitational interactions with the central (most) massive BH. The eventual fate of a sub-population of stellar-mass BHs (with or without binary companions) is for all BHs to be ejected from the cluster by the central IMBH, leaving only the most massive stellar-mass BH behind to form a close binary with the IMBH. During each phase of evolution, we discuss the rate of inspiral of the central BH-BH pair as a function of both the properties of the binary and its host cluster.

preprint2014arXiv

Revealing a new symbiotic X-ray binary with Gemini NIFS

We use K-band spectroscopy of the counterpart to the rapidly variable X-ray transient XMMU J174445.5-295044 to identify it as a new symbiotic X-ray binary. XMMU J174445.5-295044 has shown a hard X-ray spectrum (we verify its association with an Integral/IBIS 18-40 keV detection in 2013 using a short Swift/XRT observation), high and varying N$_H$, and rapid flares on timescales down to minutes, suggesting wind accretion onto a compact star. We observed its near-infrared counterpart using the Near-infrared Integral Field Spectrograph (NIFS) at Gemini-North, and classify the companion as ~ M2 III. We infer a distance of $3.1^{+1.8}_{-1.1}$ kpc (conservative 1-sigma errors), and therefore calculate that the observed X-ray luminosity (2-10 keV) has reached to at least 4$\times10^{34}$ erg/s. We therefore conclude that the source is a symbiotic X-ray binary containing a neutron star (or, less likely, black hole) accreting from the wind of a giant.

preprint2014arXiv

The Light Curve and Internal Magnetic Field of the Mode-Switching Pulsar PSR B0943+10

A number of radio pulsars exhibit intriguing mode-switching behavior. Recent observations of PSR B0943+10 revealed correlated radio and X-ray mode switches, providing a new avenue for understanding this class of objects. The large X-ray pulse fraction observed during the radio quiet phase (Q mode) was previously interpreted as a result of changing obscuration of X-rays by dense magnetosphere plasma. We show that the large X-ray pulse fraction can be explained by including the beaming effect of a magnetic atmosphere, while remaining consistent with the dipole field geometry constrained by radio observations. We also explore a more extreme magnetic field configuration, where a magnetic dipole displaced from the center of the star produces two magnetic polar caps of different sizes and magnetic field strengths. These models are currently consistent with data in radio and X-rays and can be tested or constrained by future X-ray observations.

preprint2014arXiv

The state of globular clusters at birth II: primordial binaries

(abridged) In this paper, we constrain the properties of primordial binary populations in Galactic globular clusters using the MOCCA Monte Carlo code for cluster evolution. Our results are compared to the observations of Milone et al. (2012) using the photometric binary populations as proxies for the true underlying distributions, in order to test the hypothesis that the data are consistent with an universal initial binary fraction near unity and the binary orbital parameter distributions of Kroupa (1995). With the exception of a few possible outliers, we find that the data are to first-order consistent with the universality hypothesis. Specifically, the present-day binary fractions inside the half-mass radius r$_{\rm h}$ can be reproduced assuming either high initial binary fractions near unity with a dominant soft binary component as in the Kroupa distribution combined with high initial densities (10$^4$-10$^6$ M$_{\odot}$ pc$^{-3}$), or low initial binary fractions ($\sim$ 5-10%) with a dominant hard binary component combined with moderate initial densities near their present-day values (10$^2$-10$^3$ M$_{\odot}$ pc$^{-3}$). This apparent degeneracy can be broken using the binary fractions outside r$_{\rm h}$- only high initial binary fractions with a significant soft component combined with high initial densities can contribute to reproducing the observed anti-correlation between the binary fractions outside r$_{\rm h}$ and the total cluster mass. We further illustrate using the simulated present-day binary orbital parameter distributions and the technique introduced in Leigh et al. (2012) that the relative fractions of hard and soft binaries can be used to further constrain the initial cluster density and mass-density relation. Our results favour an initial mass-density relation of the form r$_{\rm h} \propto$ M$_{\rm clus}^α$ with $α<$ 1/3.

preprint2013arXiv

A Deep Chandra X-ray Limit on the Putative IMBH in Omega Centauri

We report a sensitive X-ray search for the proposed intermediate mass black hole (IMBH) in the massive Galactic cluster, Omega Centauri (NGC 5139). Combining Chandra X-ray Observatory data from Cycles 1 and 13, we obtain a deep (~291 ks) exposure of the central regions of the cluster. We find no evidence for an X-ray point source near any of the cluster's proposed dynamical centers, and place an upper limit on the X-ray flux from a central source of f_X(0.5-7.0 keV) <= 5.0x10^-16 erg cm^-2 s^-1, after correcting for absorption. This corresponds to an unabsorbed X-ray luminosity of L_X(0.5-7.0 keV) <= 1.6x10^30 erg s^-1, for a cluster distance of 5.2 kpc, Galactic column density N_H = 1.2x10^21 cm^-2, and powerlaw spectrum with Gamma = 2.3. If a ~10^4 M_sun IMBH resides in the cluster's core, as suggested by some stellar dynamical studies, its Eddington luminosity would be L_Edd ~10^42 erg s^-1. The new X-ray limit would then establish an Eddington ratio of L_X/L_Edd <~ 10^-12, a factor of ~10 lower than even the quiescent state of our Galaxy's notoriously inefficient supermassive black hole Sgr A*, and imply accretion efficiencies as low as eta <~ 10^-6 - 10^-8. This study leaves open three possibilities: either Omega Cen does not harbor an IMBH or, if an IMBH does exist, it must experience very little or very inefficient accretion.

preprint2013arXiv

Constraints on physics of neutron stars from X-ray observations

I summarize some constraints on the physics of neutron stars arising from X-ray observations of the surfaces of neutron stars, focusing on using models of low-magnetic-field neutron star atmospheres to interpret their X-ray spectra. I discuss observations of spectral lines, pulsation profiles, X-ray bursts, radius measurements of transiently accreting neutron stars in quiescence, crust and core cooling measurements of transiently accreting neutron stars, and cooling of young neutron stars. These observations have constrained the neutron star mass and radius (and thus the internal composition, and dense matter equation of state), the superfluidity and neutrino emissivity properties of the core, and the composition and superfluid state of the crust.

preprint2013arXiv

Discovery of the Third Transient X-ray Binary in the Galactic Globular Cluster Terzan 5

We report and study the outburst of a new transient X-ray binary (XRB) in Terzan 5, the third detected in this globular cluster, Swift J174805.3-244637 or Terzan 5 X-3. We find clear spectral hardening in Swift/XRT data during the outburst rise to the hard state, thanks to our early coverage (starting at L_X ~ 4x10^{34} ergs/s) of the outburst. This hardening appears to be due to the decline in relative strength of a soft thermal component from the surface of the neutron star (NS) during the rise. We identify a {Type I X-ray burst} in Swift/XRT data with a long (16 s) decay time, indicative of {hydrogen burning on the surface of the} NS. We use Swift/BAT, Maxi/GSC, Chandra/ACIS, and Swift/XRT data to study the spectral changes during the outburst, identifying a clear hard-to-soft state transition. We use a Chandra/ACIS observation during outburst to identify the transient's position. Seven archival Chandra/ACIS observations show evidence for variations in Terzan 5 X-3's non-thermal component, but not the thermal component, during quiescence. The inferred long-term time-averaged mass accretion rate, from the quiescent thermal luminosity, suggests that if this outburst is typical and only slow cooling processes are active in the neutron star core, such outbursts should recur every ~10 years.

preprint2013arXiv

Found: The Progenitors of AM CVn and Supernovae .Ia

We present optical and X-ray observations of two tidally distorted, extremely low-mass white dwarfs (WDs) with massive companions. There is no evidence of neutron stars in our Chandra and XMM observations of these objects. SDSS J075141.18$-$014120.9 (J0751) is an eclipsing double WD binary containing a 0.19 Msol WD with a 0.97 Msol companion in a 1.9 h orbit. J0751 becomes the fifth eclipsing double WD system currently known. SDSS J174140.49+652638.7 (J1741) is another binary containing a 0.17 Msol WD with an unseen M > 1.11 Msol WD companion in a 1.5 h orbit. With a mass ratio of ~0.1, J1741 will have stable mass transfer through an accretion disk and turn into an interacting AM Canum Venaticorum (AM CVn) system in the next ~160 Myr. With a mass ratio of 0.2, J0751 is likely to follow a similar evolutionary path. These are the first known AM CVn progenitor binary systems and they provide important constraints on the initial conditions for AM CVn. Theoretical studies suggest that both J0751 and J1741 may create thermonuclear supernovae in ~10^8 yr, either .Ia or Ia. Such explosions can account for ~1% of the Type Ia supernova rate.

preprint2013arXiv

Mass/Radius Constraints on the Quiescent Neutron Star in M13 Using Hydrogen and Helium Atmospheres

The mass and radius of the neutron star (NS) in low-mass X-ray binaries can be obtained by fitting the X-ray spectrum of the NS in quiescence, and the mass and radius constrains the properties of dense matter in NS cores. A critical ingredient for spectral fits is the composition of the NS atmosphere: hydrogen atmospheres are assumed in most prior work, but helium atmospheres are possible if the donor star is a helium white dwarf. Here we perform spectral fits to XMM, Chandra, and ROSAT data of a quiescent NS in the globular cluster M13. This NS has the smallest inferred radius from previous spectral fitting. Assuming an atmosphere composed of hydrogen, we find a significantly larger radius, more consistent with those from other quiescent NSs. With a helium atmosphere (an equally acceptable fit), we find even larger values for the radius.

preprint2013arXiv

Optical counterparts of the nearest ultraluminous X-ray sources

We present a photometric survey of the optical counterparts of ultraluminous X-ray sources (ULXs) observed with the Hubble Space Telescope in nearby (< 5 Mpc) galaxies. Of the 33 ULXs with Hubble & Chandra data, 9 have no visible counterpart, placing limits on their M_V of ~ -4 to -9, enabling us to rule out O-type companions in 4 cases. The refined positions of two ULXs place them in the nucleus of their host galaxy. They are removed from our sample. Of the 22 remaining ULXs, 13 have one possible optical counterpart, while multiple are visible within the error regions of other ULXs. By calculating the number of chance coincidences, we estimate that 13 +/- 5 are the true counterparts. We attempt to constrain the nature of the companions by fitting the SED and M_V to obtain candidate spectral types. We can rule out O-type companions in 20 cases, while we find that one ULX (NGC 253 ULX2) excludes all OB-type companions. Fitting with X-ray irradiated models provides constraints on the donor star mass and radius. For 7 ULXs, we are able to impose inclination-dependent upper and/or lower limits on the black holes mass, if the extinction to the assumed companion star is not larger than the Galactic column. These are NGC 55 ULX1, NGC 253 ULX1, NGC 253 ULX2, NGC 253 XMM6, Ho IX X-1, IC342 X-1 & NGC 5204 X-1. This suggests that 10 ULXs do not have O companions, while none of the 18 fitted rule out B-type companions.

preprint2013arXiv

Stellar Encounter Rate in Galactic Globular Clusters

The high stellar densities in the cores of globular clusters cause significant stellar interactions. These stellar interactions can produce close binary mass-transferring systems involving compact objects and their progeny, such as X-ray binaries and radio millisecond pulsars. Comparing the numbers of these systems and interaction rates in different clusters drives our understanding of how cluster parameters affect the production of close binaries. In this paper we estimate stellar encounter rates (Gamma) for 124 Galactic globular clusters based on observational data as opposed to the methods previously employed, which assumed "King-model" profiles for all clusters. By deprojecting cluster surface brightness profiles to estimate luminosity density profiles, we treat "King-model" and "core-collapsed" clusters in the same way. In addition, we use Monte-Carlo simulations to investigate the effects of uncertainties in various observational parameters (distance, reddening, surface brightness) on Gamma, producing the first catalog of GC stellar encounter rates with estimated errors. Comparing our results with published observations of likely products of stellar interactions (numbers of X-ray binaries, numbers of radio millisecond pulsars, and gamma-ray luminosity) we find both clear correlations and some differences with published results.

preprint2013arXiv

The hottest superfluid and superconductor in the Universe: Discovery and nuclear physics implications

We present recent work on using astronomical observations of neutron stars to reveal unique insights into nuclear matter that cannot be obtained from laboratories on Earth. First, we discuss our measurement of the rapid cooling of the youngest neutron star in the Galaxy; this provides the first direct evidence for superfluidity and superconductivity in the supra-nuclear core of neutron stars. We show that observations of thermonuclear X-ray bursts on neutron stars can be used to constrain properties of neutron superfluidity and neutrino emission. We describe the implications of rapid neutron star rotation rates on aspects of nuclear and superfluid physics. Finally, we show that entrainment coupling between the neutron superfluid and the nuclear lattice leads to a less mobile crust superfluid; this result puts into question the conventional picture of pulsar glitches as being solely due to the crust superfluid and suggests that the core superfluid also participates.

preprint2013arXiv

The Runaway Binary LP 400-22 is Leaving the Galaxy

We present optical spectroscopy, astrometry, radio, and X-ray observations of the runaway binary LP 400-22. We refine the orbital parameters of the system based on our new radial velocity observations. Our parallax data indicate that LP 400-22 is significantly more distant (3 sigma lower limit of 840 pc) than initially predicted. LP 400-22 has a tangential velocity in excess of 830 km/s; it is unbound to the Galaxy. Our radio and X-ray observations fail to detect a recycled millisecond pulsar companion, indicating that LP 400-22 is a double white dwarf system. This essentially rules out a supernova runaway ejection mechanism. Based on its orbit, a Galactic center origin is also unlikely. However, its orbit intersects the locations of several globular clusters; dynamical interactions between LP 400-22 and other binary stars or a central black hole in a dense cluster could explain the origin of this unusual binary.

preprint2013arXiv

The X-ray Properties of Five Galactic Supernova Remnants Detected by the Spitzer GLIMPSE Survey

(Abbreviated) We present a study of the X-ray properties of five Galactic supernova remnants (SNRs) -- Kes 17 (G304.6$+$0.1), G311.5$-$0.3, G346.6$-$0.2, CTB 37A (G348.5$+$0.1) and G348.5$-$0.0 -- that were detected in the infrared by Reach et al. (2006) in an analysis of data from the Galactic Legacy Infrared Mid-Plane Survey Extraordinaire (GLIMPSE) that was conducted by the Spitzer Space Telescope. We present and analyze archival ASCA observations of Kes 17, G311.5$-$0.3 and G346.6$-$0.2, archival XMM-Newton observations of Kes 17, CTB 37A and G348.5$-$0.0 and an archival Chandra observation of CTB 37A. All of the SNRs are clearly detected in the X-ray possibly except for G348.5$-$0.0. Our study reveals that the four detected SNRs all feature center-filled X-ray morphologies and that the observed emission from these sources is thermal in all cases. We argue that these SNRs should be classified as mixed-morphology SNRs (MM SNRs): our study strengthens the correlation between MM SNRs and SNRs interacting with molecular clouds and suggests that the origin of mixed-morphology SNRs may be due to the interactions between these SNRs and adjacent clouds. Our ASCA analysis of G311.5$-$0.3 reveals for the first time X-ray emission from this SNR: the X-ray emission is center-filled within the radio and infrared shells and thermal in nature ($kT$ $\sim$ 0.98 keV), thus motivating its classification as an MM SNR. We find considerable spectral variations in the properties associated with the plasmas of the other X-ray-detected SNRs, such as a possible overabundance of magnesium in the plasma of Kes 17. Finally, we also estimate such properties as electron density $n$$_e$, radiative age $t$$_{rad}$ and swept-up mass $M$$_X$ for each of the four X-ray-detected SNRs.

preprint2012arXiv

A 2.15 Hour Orbital Period for the Low Mass X-Ray Binary XB 1832-330 in the Globular Cluster NGC 6652

We present a candidate orbital period for the low mass X-ray binary XB 1832-330 in the globular cluster NGC 6652 using a 6.5 hour Gemini South observation of the optical counterpart of the system. Light curves in g' and r' for two LMXBs in the cluster, sources A and B in previous literature, were extracted and analyzed for periodicity using the ISIS image subtraction package. A clear sinusoidal modulation is evident in both of A's curves, of amplitude ~0.11 magnitudes in g' and ~0.065 magnitudes in r', while B's curves exhibit rapid flickering, of amplitude ~1 magnitude in g' and ~0.5 magnitudes in r'. A Lomb-Scargle test revealed a 2.15 hour periodic variation in the magnitude of A with a false alarm probability less than 10^-11, and no significant periodicity in the light curve for B. Though it is possible saturated stars in the vicinity of our sources partially contaminated our signal, the identification of A's binary period is nonetheless robust.

preprint2012arXiv

The ELM Survey. IV. 24 White Dwarf Merger Systems

We present new radial velocity and X-ray observations of extremely low-mass (ELM, 0.2 Msol) white dwarf candidates in the Sloan Digital Sky Survey (SDSS) Data Release 7 area. We identify seven new binary systems with 1-18 h orbital periods. Five of the systems will merge due to gravitational wave radiation within 10 Gyr, bringing the total number of merger systems found in the ELM Survey to 24. The ELM Survey has now quintupled the known merger white dwarf population. It has also discovered the eight shortest period detached binary white dwarf systems currently known. We discuss the characteristics of the merger and non-merger systems observed in the ELM Survey, including their future evolution. About half of the systems have extreme mass ratios. These are the progenitors of the AM Canum Venaticorum systems and supernovae .Ia. The remaining targets will lead to the formation of extreme helium stars, subdwarfs, or massive white dwarfs. We identify three targets that are excellent gravitational wave sources. These should be detected by the Laser Interferometer Space Antenna (LISA)-like missions within the first year of operation. The remaining targets are important indicators of what the Galactic foreground may look like for gravitational wave observatories.

preprint2012arXiv

X-ray Binaries in the Ultrahigh Encounter Rate Globular Cluster NGC 6388

We report the results of a joint \chandra-\hst study of the X-ray binary population in the massive, high-density globular cluster NGC 6388. NGC 6388 has one of the highest predicted X-ray binary production rate of any Galactic cluster. We detected a large population of 61 \chandra sources within the half-mass radius with L$_X > 5 \times 10^{30}$ \ergs. From the X-ray colors, luminosities, (lack of) variability, and spectral fitting, we identify five as likely quiescent low-mass X-ray binaries. Due to the extremely crowded nature of the core of NGC 6388, finding optical identifications to \chandra sources is challenging. We have identified four blue, optically variable counterparts to spectrally hard X-ray sources, evidence that these are bright cataclysmic variables (CVs). One showed variability of 2 magnitudes in V, indicative of a dwarf nova eruption. One other likely CV is identified by its X-ray spectrum (partial covering with high $N_H$) and strong variability, making five likely CVs identified in this cluster. The relatively bright optical magnitudes of these sources put them in the same class as CV1 in M15 and the brightest CVs in 47 Tuc.

preprint2012arXiv

XMM-Newton Finds That SAX J1750.8-2900 May Harbor the Hottest, Most Luminous Known Neutron Star

We have performed the first sensitive X-ray observation of the low-mass X-ray binary SAX J1750.8-2900 in quiescence with XMM-Newton. The spectrum was fit to both a classical black body model, and a non-magnetized, pure hydrogen neutron star atmosphere model. A power law component was added to these models, but we found that it was not required by the fits. The distance to SAX J1750.8-2900 is known to be D = 6.79 kpc from a previous analysis of photospheric radius expansion bursts. This distance implies a bolometric luminosity (as given by the NS atmosphere model) of (1.05 +/- 0.12) x 10^34 (D/6.79 kpc)^2 erg s^-1, which is the highest known luminosity for a NS LMXB in quiescence. One simple explanation for this surprising result could be that the crust and core of the NS were not in thermal equilibrium during the observation. We argue that this was likely not the case, and that the core temperature of the NS in SAX J1750.8-2900 is unusually high.

preprint2011arXiv

Chandra X-ray Observations of 12 Millisecond Pulsars in the Globular Cluster M28

We present a Chandra X-ray Observatory investigation of the millisecond pulsars (MSPs) in the globular cluster M28 (NGC 6626). In what is one of the deepest X-ray observations of a globular cluster, we firmly detect seven and possibly detect two of the twelve known M28 pulsars. With the exception of PSRs B1821-24 and J1824-2452H, the detected pulsars have relatively soft spectra, with X-ray luminosities 10^30-31 ergs s^-1 (0.3-8 keV),similar to most "recycled" pulsars in 47 Tucanae and the field of the Galaxy, implying thermal emission from the pulsar magnetic polar caps. We present the most detailed X-ray spectrum to date of the energetic PSR B1821-24. It is well described by a purely non-thermal spectrum with spectral photon index 1.23 and luminosity 1.4x10^33Theta(D/5.5 kpc)^2 ergs s^-1 (0.3-8 keV), where Theta is the fraction of the sky covered by the X-ray emission beam(s). We find no evidence for the previously reported line emission feature around 3.3 keV, most likely as a consequence of improvements in instrument calibration. The X-ray spectrum and pulse profile of PSR B1821--24 suggest that the bulk of unpulsed emission from this pulsar is not of thermal origin, and is likely due to low-level non-thermal magnetospheric radiation, an unresolved pulsar wind nebula, and/or small-angle scattering of the pulsed X-rays by interstellar dust grains. The peculiar binary PSR J1824-2452H shows a relatively hard X-ray spectrum and possible variability at the binary period, indicative of an intrabinary shock formed by interaction between the relativistic pulsar wind and matter from its non-degenerate companion star.

preprint2011arXiv

Cooling neutron star in the Cassiopeia~A supernova remnant: Evidence for superfluidity in the core

According to recent results of Ho & Heinke (2009) and Heinke & Ho (2010), the Cassiopeia A supernova remnant contains a young neutron star which has carbon atmosphere and shows noticeable decline of the effective surface temperature. We report a new (November 2010) Chandra observation which confirms the previously reported decline rate. The decline is naturally explained if neutrons have recently become superfluid (in triplet-state) in the NS core, producing a splash of neutrino emission due to Cooper pair formation (CPF) process that currently accelerates the cooling. This scenario puts stringent constraints on poorly known properties of NS cores: on density dependence of the temperature $T_\mathrm{cn}(ρ)$ for the onset of neutron superfluidity [$T_\mathrm{cn}(ρ)$ should have a wide peak with maximum $\approx (7-9)\times 10^8$ K], on the reduction factor $q$ of CPF process by collective effects in superfluid matter ($q > 0.4$), and on the intensity of neutrino emission before the onset of neutron superfluidity (30--100 times weaker than the standard modified Urca process). This is serious evidence for nucleon superfluidity in NS cores that comes from observations of cooling NSs.

preprint2011arXiv

Identification of Faint Chandra X-ray Sources in the Core-Collapsed Globular Cluster NGC 6397: Evidence for a Bimodal Cataclysmic Variable Population

We have searched for optical identifications for 79 Chandra X-ray sources that lie within the half-mass radius of the nearby, core-collapsed globular cluster NGC 6397, using deep Hubble Space Telescope Advanced Camera for Surveys Wide Field Channel imaging in H-alpha, R, and B. Photometry of these images allows us to classify candidate counterparts based on color-magnitude diagram location. In addition to recovering nine previously detected cataclysmic variables (CVs), we have identified six additional faint CV candidates, a total of 42 active binaries (ABs), two millisecond pulsars (MSPs), one candidate active galactic nucleus, and one candidate interacting galaxy pair. Of the 79 sources, 69 have a plausible optical counterpart. The 15 likely and possible CVs in NGC 6397 mostly fall into two groups: a brighter group of six for which the optical emission is dominated by contributions from the secondary and accretion disk, and a fainter group of seven for which the white dwarf dominates the optical emission. There are two possible transitional objects that lie between these groups. The faintest CVs likely lie near the minimum of the CV period distribution, where an accumulation is expected. The spatial distribution of the brighter CVs is much more centrally concentrated than those of the fainter CVs and the active binaries. This may represent the result of an evolutionary process in which CVs are produced by dynamical interactions, such as exchange reactions, near the cluster center and are scattered to larger orbital radii, over their lifetimes, as they age and become fainter.

preprint2011arXiv

X-ray Sources in Galactic Globular Clusters

I review recent work on X-ray sources in Galactic globular clusters, identified with low-mass X-ray binaries (LMXBs), cataclysmic variables (CVs), millisecond pulsars (MSPs) and coronally active binaries by Chandra. Faint transient LMXBs have been identified in several clusters, challenging our understanding of accretion disk instabilities. Spectral fitting of X-rays from quiescent LMXBs offers the potential to constrain the interior structure of neutron stars. The numbers of quiescent LMXBs scale with the dynamical interaction rates of their host clusters, indicating their dynamical formation. Large numbers of CVs have been discovered, including a very faint population in NGC 6397 that may be at or beyond the CV period minimum. Most CVs in dense clusters seem to be formed in dynamical interactions, but there is evidence that some are primordial binaries. Radio millisecond pulsars show thermal X-rays from their polar caps, and often nonthermal X-rays, either from magnetospheric emission, or from a shock between the pulsar wind and material still flowing from the companion. Chromospherically active binaries comprise the largest number of X-ray sources in globular clusters, and their numbers generally scale with cluster mass, but their numbers seem to be reduced in all globular clusters compared to other old stellar populations.

preprint2010arXiv

Cooling rates of neutron stars and the young neutron star in the Cassiopeia A supernova remnant

We explore the thermal state of the neutron star in the Cassiopeia A supernova remnant using the recent result of Ho & Heinke (Nature, 462, 71 (2009)) that the thermal radiation of this star is well-described by a carbon atmosphere model and the emission comes from the entire stellar surface. Starting from neutron star cooling theory, we formulate a robust method to extract neutrino cooling rates of thermally relaxed stars at the neutrino cooling stage from observations of thermal surface radiation. We show how to compare these rates with the rates of standard candles -- stars with non-superfluid nucleon cores cooling slowly via the modified Urca process. We find that the internal temperature of standard candles is a well-defined function of the stellar compactness parameter $x=r_g/R$, irrespective of the equation of state of neutron star matter ($R$ and $r_g$ are circumferential and gravitational radii, respectively). We demonstrate that the data on the Cassiopeia A neutron star can be explained in terms of three parameters: $f_\ell$, the neutrino cooling efficiency with respect to the standard candle; the compactness $x$; and the amount of light elements in the heat blanketing envelope. For an ordinary (iron) heat blanketing envelope or a low-mass ($\lesssim 10^{-13}\,M_\odot$) carbon envelope, we find the efficiency $f_\ell \sim 1$ (standard cooling) for $x \lesssim 0.5$ and $f_\ell \sim 0.02$ (slower cooling) for a maximum compactness $x\approx 0.7$. A heat blanket containing the maximum mass ($\sim 10^{-8}\,M_\odot$) of light elements increases $f_\ell$ by a factor of 50. We also examine the (unlikely) possibility that the star is still thermally non-relaxed.

preprint2010arXiv

Direct Observation of the Cooling of the Cassiopeia A Neutron Star

The cooling rate of young neutron stars gives direct insight into their internal makeup. Although the temperatures of several young neutron stars have been measured, until now a young neutron star has never been observed to decrease in temperature over time. We fit 9 years of archival Chandra ACIS spectra of the likely neutron star in the ~330 years old Cassiopeia A supernova remnant with our non-magnetic carbon atmosphere model. Our fits show a relative decline in the surface temperature by 4% (5.4 sigma, from 2.12+-0.01*10^6 K in 2000 to 2.04+-0.01*10^6 K in 2009) and observed flux (by 21%). Using a simple model for neutron star cooling, we show that this temperature decline could indicate that the neutron star became isothermal sometime between 1965 and 1980, and constrains some combinations of neutrino emission mechanisms and envelope compositions. However, the neutron star is likely to have become isothermal soon after formation, in which case the temperature history suggests episodes of additional heating or more rapid cooling. Observations over the next few years will allow us to test possible explanations for the temperature evolution.

preprint2009arXiv

A Chandra X-ray Observatory Study of PSR J1740--5340 and Candidate Millisecond Pulsars in the Globular Cluster NGC 6397

We present a deep Chandra X-ray Observatory study of the peculiar binary radio millisecond pulsar PSR J1740--5340 and candidate millisecond pulsars (MSPs) in the globular cluster NGC 6397. The X-rays from PSR J1740--5340 appear to be non-thermal and exhibit variability at the binary period. These properties suggest the presence of a relativistic intrabinary shock formed due to interaction of a relativistic rotation-powered pulsar wind and outflow from the unusual "red-straggler/sub-subgiant" companion. We find the X-ray source U18 to show similar X-ray and optical properties to those of PSR J1740--5340, making it a strong MSP candidate. It exhibits variability on timescales from hours to years, also consistent with an intrabinary shock origin of its X-ray emission. The unprecedented depth of the X-ray data allows us to conduct a complete census of MSPs in NGC 6397. Based on the properties of the present sample of X-ray--detected MSPs in the Galaxy we find that NGC 6397 probably hosts no more than 6 MSPs.

preprint2009arXiv

A Neutron Star with a Carbon Atmosphere in the Cassiopeia A Supernova Remnant

The surface of hot neutron stars is covered by a thin atmosphere. If there is accretion after neutron star formation, the atmosphere could be composed of light elements (H or He); if no accretion takes place or if thermonuclear reactions occur after accretion, heavy elements (for example, Fe) are expected. Despite detailed searches, observations have been unable to confirm the atmospheric composition of isolated neutron stars. Here we report an analysis of archival observations of the compact X-ray source in the centre of the Cassiopeia A supernova remnant. We show that a carbon atmosphere neutron star (with low magnetic field) produces a good fit to the spectrum. Our emission model, in contrast with others, implies an emission size consistent with theoretical predictions for the radius of neutron stars. This result suggests that there is nuclear burning in the surface layers and also identifies the compact source as a very young (~330-year-old) neutron star.

preprint2008arXiv

The Discovery of a Very Faint X-ray Transient in the Globular Cluster M15

We have identified an X-ray transient (hereafter M15 X-3) in the globular cluster M15 from an archival Chandra grating observation. M15 X-3 appears at an X-ray luminosity of 6*10^{33} ergs/s with a spectrum consistent with an absorbed power law of photon index 1.51+-0.14. The object is identifiable in archival Chandra HRC-I observations with an X-ray luminosity of 2-6*10^{31} ergs/s and apparently soft colors, suggesting a neutron star low-mass X-ray binary in quiescence. We also observe it in outburst in a 2007 Chandra HRC-I observation, and in archival 1994-1995 ROSAT HRI observations. We identify a likely optical/UV counterpart with a (possibly transient) UV excess from archival HST data, which suggests a main sequence companion. We argue that M15 X-3's behavior is similar to that of the very faint X-ray transients which have been observed in the Galactic Center. We discuss several explanations for its very low X-ray luminosity, with the assumption that we have detected its companion. M15 X-3's uniquely low extinction and well-determined distance make it an excellent target for future studies.