Black hole hunting in the Andromeda Galaxy

We present a new technique for identifying stellar mass black holes in low mass X-ray binaries (LMXBs), and apply it to XMM-Newton observations of M31. We examine X-ray time series variability seeking power density spectra (PDS) typical of LMXBs accreting at a low accretion rate (which we refer to as Type A PDS); these are very similar for black hole and neutron star LMXBs. Galactic neutron star LMXBs exhibit Type A PDS at low luminosities (~10^36--10^37 erg/s) while black hole LMXBs can exhibit them at luminosities >10^38 erg/s. We propose that Type A PDS are confined to luminosities below a critical fraction of the Eddington limit, $l_c$ that is constant for all LMXBs; we have examined asample of black hole and neutron star LMXBs and find they are all consistent with $l_c$ = 0.10+/-0.04 in the 0.3--10 keV band. We present luminosity and PDS data from 167 observations of X-ray binaries in M31 that provide strong support for our hypothesis. Since the theoretical maximum mass for a neutron star is \~3.1 M_Sun, we therefore assert that any LMXB that exhibits a Type A PDS at a 0.3--10 keV luminosity greater than 4 x 10^37 erg/s is likely to contain a black hole primary. We have found eleven new black hole candidates in M31 using this method. We focus on XMM-Newton observations of RX J0042.4+4112, an X-ray source in M31 and find the mass of the primary to be 7+/-2 M_Sun, if our assumptions are correct. Furthermore, RX J0042.4+4112 is consistently bright in \~40 observations made over 23 years, and is likely to be a persistently bright LMXB; by contrast all known Galactic black hole LMXBs are transient. Hence our method may be used to find black holes in known, persistently bright Galactic LMXBs and also in LMXBs in other galaxies.Comment: 6 Pages, 6 figures. To appear in the conference proceedings of "Interacting Binaries: Accretion, Evolution and Outcomes" (Cefalu, July 4-10 2004

Exploring Interacting Quantum Many-Body Systems by Experimentally Creating Continuous Matrix Product States in Superconducting Circuits

Improving the understanding of strongly correlated quantum many body systems such as gases of interacting atoms or electrons is one of the most important challenges in modern condensed matter physics, materials research and chemistry. Enormous progress has been made in the past decades in developing both classical and quantum approaches to calculate, simulate and experimentally probe the properties of such systems. In this work we use a combination of classical and quantum methods to experimentally explore the properties of an interacting quantum gas by creating experimental realizations of continuous matrix product states - a class of states which has proven extremely powerful as a variational ansatz for numerical simulations. By systematically preparing and probing these states using a circuit quantum electrodynamics (cQED) system we experimentally determine a good approximation to the ground-state wave function of the Lieb-Liniger Hamiltonian, which describes an interacting Bose gas in one dimension. Since the simulated Hamiltonian is encoded in the measurement observable rather than the controlled quantum system, this approach has the potential to apply to exotic models involving multicomponent interacting fields. Our findings also hint at the possibility of experimentally exploring general properties of matrix product states and entanglement theory. The scheme presented here is applicable to a broad range of systems exploiting strong and tunable light-matter interactions.Comment: 11 pages, 9 figure

A search for thermal X-ray signatures in Gamma-Ray Bursts I: Swift bursts with optical supernovae

The X-ray spectra of Gamma-Ray Bursts can generally be described by an absorbed power law. The landmark discovery of thermal X-ray emission in addition to the power law in the unusual GRB 060218, followed by a similar discovery in GRB 100316D, showed that during the first thousand seconds after trigger the soft X-ray spectra can be complex. Both the origin and prevalence of such spectral components still evade understanding, particularly after the discovery of thermal X-ray emission in the classical GRB 090618. Possibly most importantly, these three objects are all associated with optical supernovae, begging the question of whether the thermal X-ray components could be a result of the GRB-SN connection, possibly in the shock breakout. We therefore performed a search for blackbody components in the early Swift X-ray spectra of 11 GRBs that have or may have associated optical supernovae, accurately recovering the thermal components reported in the literature for GRBs 060218, 090618 and 100316D. We present the discovery of a cooling blackbody in GRB 101219B/SN2010ma, and in four further GRB-SNe we find an improvement in the fit with a blackbody which we deem possible blackbody candidates due to case-specific caveats. All the possible new blackbody components we report lie at the high end of the luminosity and radius distribution. GRB 101219B appears to bridge the gap between the low-luminosity and the classical GRB-SNe with thermal emission, and following the blackbody evolution we derive an expansion velocity for this source of order 0.4c. We discuss potential origins for the thermal X-ray emission in our sample, including a cocoon model which we find can accommodate the more extreme physical parameters implied by many of our model fits.Comment: 16 pages, 6 figures, accepted for MNRA

Discovery of disc precession in the M31 dipping X-ray binary Bo 158

We present results from three XMM-Newton observations of the M31 low mass X-ray binary XMMU J004314.4+410726.3 (Bo 158), spaced over 3 days in 2004, July. Bo 158 was the first dipping LMXB to be discovered in M31. Periodic intensity dips were previously seen to occur on a 2.78-hr period, due to absorption in material that is raised out of the plane of the accretion disc. The report of these observations stated that the dip depth was anti-correlated with source intensity. However, our new observations do not favour a strict intensity dependance, but rather suggest that the dip variation is due to precession of the accretion disc. This is to be expected in LMXBs with a mass ratio <~ 0.3 (period <~ 4 hr), as the disc reaches the 3:1 resonance with the binary companion, causing elongation and precession of the disc. A smoothed particle hydrodynamics simulation of the disc in this system shows retrograde rotation of a disc warp on a period of ~11 P_orb, and prograde disc precession on a period of ~29 P_orb. This is consistent with the observed variation in the depth of the dips. We find that the dipping behaviour is most likely to be modified by the disc precession, hence we predict that the dipping behaviour repeats on a 81+/-3 hr cycle.Comment: 9 pages, 6 figures, accepted for publication by MNRAS, changed conten

Space tug propulsion system failure mode, effects and criticality analysis

For purposes of the study, the propulsion system was considered as consisting of the following: (1) main engine system, (2) auxiliary propulsion system, (3) pneumatic system, (4) hydrogen feed, fill, drain and vent system, (5) oxygen feed, fill, drain and vent system, and (6) helium reentry purge system. Each component was critically examined to identify possible failure modes and the subsequent effect on mission success. Each space tug mission consists of three phases: launch to separation from shuttle, separation to redocking, and redocking to landing. The analysis considered the results of failure of a component during each phase of the mission. After the failure modes of each component were tabulated, those components whose failure would result in possible or certain loss of mission or inability to return the Tug to ground were identified as critical components and a criticality number determined for each. The criticality number of a component denotes the number of mission failures in one million missions due to the loss of that component. A total of 68 components were identified as critical with criticality numbers ranging from 1 to 2990

The SSS phase of RS Ophiuchi observed with Chandra and XMM-Newton I.: Data and preliminary Modeling

The phase of Super-Soft-Source (SSS) emission of the sixth recorded outburst of the recurrent nova RS Oph was observed twice with Chandra and once with XMM-Newton. The observations were taken on days 39.7, 54.0, and 66.9 after outburst. We confirm a 35-sec period on day 54.0 and found that it originates from the SSS emission and not from the shock. We discus the bound-free absorption by neutral elements in the line of sight, resonance absorption lines plus self-absorbed emission line components, collisionally excited emission lines from the shock, He-like intersystem lines, and spectral changes during an episode of high-amplitude variability. We find a decrease of the oxygen K-shell absorption edge that can be explained by photoionization of oxygen. The absorption component has average velocities of -1286+-267 km/s on day 39.7 and of -771+-65 km/s on day 66.9. The wavelengths of the emission line components are consistent with their rest wavelengths as confirmed by measurements of non-self absorbed He-like intersystem lines. We have evidence that these lines originate from the shock rather than the outer layers of the outflow and may be photoexcited in addition to collisional excitations. We found collisionally excited emission lines that are fading at wavelengths shorter than 15A that originate from the radiatively cooling shock. On day 39.5 we find a systematic blue shift of -526+-114 km/s from these lines. We found anomalous He-like f/i ratios which indicates either high densities or significant UV radiation near the plasma where the emission lines are formed. During the phase of strong variability the spectral hardness light curve overlies the total light curve when shifted by 1000sec. This can be explained by photoionization of neutral oxygen in the line of sight if the densities of order 10^{10}-10^{11} cm^{-3}.Comment: 16 pages, 10 figures, 4 tables. Accepted by ApJ; v2: Co-author Woodward adde

Longitudinal LASSO: Jointly Learning Features and Temporal Contingency for Outcome Prediction

Longitudinal analysis is important in many disciplines, such as the study of behavioral transitions in social science. Only very recently, feature selection has drawn adequate attention in the context of longitudinal modeling. Standard techniques, such as generalized estimating equations, have been modified to select features by imposing sparsity-inducing regularizers. However, they do not explicitly model how a dependent variable relies on features measured at proximal time points. Recent graphical Granger modeling can select features in lagged time points but ignores the temporal correlations within an individual's repeated measurements. We propose an approach to automatically and simultaneously determine both the relevant features and the relevant temporal points that impact the current outcome of the dependent variable. Meanwhile, the proposed model takes into account the non-{\em i.i.d} nature of the data by estimating the within-individual correlations. This approach decomposes model parameters into a summation of two components and imposes separate block-wise LASSO penalties to each component when building a linear model in terms of the past $\tau$ measurements of features. One component is used to select features whereas the other is used to select temporal contingent points. An accelerated gradient descent algorithm is developed to efficiently solve the related optimization problem with detailed convergence analysis and asymptotic analysis. Computational results on both synthetic and real world problems demonstrate the superior performance of the proposed approach over existing techniques.Comment: Proceedings of the 21th ACM SIGKDD International Conference on Knowledge Discovery and Data Mining. ACM, 201

1SXPS: A deep Swift X-ray Telescope point source catalog with light curves and spectra

We present the 1SXPS (Swift-XRT Point Source) catalog of 151,524 X-ray point-sources detected by the Swift-XRT in 8 years of operation. The catalog covers 1905 square degrees distributed approximately uniformly on the sky. We analyze the data in two ways. First we consider all observations individually, for which we have a typical sensitivity of ~3e-13 erg/cm2/s (0.3--10 keV). Then we co-add all data covering the same location on the sky: these images have a typical sensitivity of ~9e-14 erg/cm2/s (0.3--10 keV). Our sky coverage is nearly 2.5 times that of 3XMM-DR4, although the catalog is a factor of ~1.5 less sensitive. The median position error is 5.5" (90% confidence), including systematics. Our source detection method improves on that used in previous XRT catalogs and we report >68,000 new X-ray sources. The goals and observing strategy of the Swift satellite allow us to probe source variability on multiple timescales, and we find ~30,000 variable objects in our catalog. For every source we give positions, fluxes, time series (in four energy bands and two hardness ratios), estimates of the spectral properties, spectra and spectral fits for the brightest sources, and variability probabilities in multiple energy bands and timescales.Comment: 27 pages, 19 figures; accepted for publication in ApJS. The accompanying website, http://www.swift.ac.uk/1SXPS is live; the Vizier entry should be available shortl

Early multi-wavelength emission from Gamma-ray Bursts: from Gamma-ray to X-ray

The study of the early high-energy emission from both long and short Gamma-ray bursts has been revolutionized by the Swift mission. The rapid response of Swift shows that the non-thermal X-ray emission transitions smoothly from the prompt phase into a decaying phase whatever the details of the light curve. The decay is often categorized by a steep-to-shallow transition suggesting that the prompt emission and the afterglow are two distinct emission components. In those GRBs with an initially steeply-decaying X-ray light curve we are probably seeing off-axis emission due to termination of intense central engine activity. This phase is usually followed, within the first hour, by a shallow decay, giving the appearance of a late emission hump. The late emission hump can last for up to a day, and hence, although faint, is energetically very significant. The energy emitted during the late emission hump is very likely due to the forward shock being constantly refreshed by either late central engine activity or less relativistic material emitted during the prompt phase. In other GRBs the early X-ray emission decays gradually following the prompt emission with no evidence for early temporal breaks, and in these bursts the emission may be dominated by classical afterglow emission from the external shock as the relativistic jet is slowed by interaction with the surrounding circum-burst medium. At least half of the GRBs observed by Swift also show erratic X-ray flaring behaviour, usually within the first few hours. The properties of the X-ray flares suggest that they are due to central engine activity. Overall, the observed wide variety of early high-energy phenomena pose a major challenge to GRB models.Comment: Accepted for publication in the New Journal of Physics focus issue on Gamma Ray Burst