Showing 1–50 of 61 results for author: Barnard, R

Speech foundation models in healthcare: Effect of layer selection on pathological speech feature prediction

Authors: Daniela A. Wiepert, Rene L. Utianski, Joseph R. Duffy, John L. Stricker, Leland R. Barnard, David T. Jones, Hugo Botha

Abstract: Accurately extracting clinical information from speech is critical to the diagnosis and treatment of many neurological conditions. As such, there is interest in leveraging AI for automatic, objective assessments of clinical speech to facilitate diagnosis and treatment of speech disorders. We explore transfer learning using foundation models, focusing on the impact of layer selection for the downst… ▽ More Accurately extracting clinical information from speech is critical to the diagnosis and treatment of many neurological conditions. As such, there is interest in leveraging AI for automatic, objective assessments of clinical speech to facilitate diagnosis and treatment of speech disorders. We explore transfer learning using foundation models, focusing on the impact of layer selection for the downstream task of predicting pathological speech features. We find that selecting an optimal layer can greatly improve performance (~15.8% increase in balanced accuracy per feature as compared to worst layer, ~13.6% increase as compared to final layer), though the best layer varies by predicted feature and does not always generalize well to unseen data. A learned weighted sum offers comparable performance to the average best layer in-distribution (only ~1.2% lower) and had strong generalization for out-of-distribution data (only 1.5% lower than the average best layer). △ Less

Submitted 21 June, 2024; v1 submitted 2 February, 2024; originally announced February 2024.

Comments: Accepted to INTERSPEECH 2024

Detecting Speech Abnormalities with a Perceiver-based Sequence Classifier that Leverages a Universal Speech Model

Authors: Hagen Soltau, Izhak Shafran, Alex Ottenwess, Joseph R. JR Duffy, Rene L. Utianski, Leland R. Barnard, John L. Stricker, Daniela Wiepert, David T. Jones, Hugo Botha

Abstract: We propose a Perceiver-based sequence classifier to detect abnormalities in speech reflective of several neurological disorders. We combine this classifier with a Universal Speech Model (USM) that is trained (unsupervised) on 12 million hours of diverse audio recordings. Our model compresses long sequences into a small set of class-specific latent representations and a factorized projection is use… ▽ More We propose a Perceiver-based sequence classifier to detect abnormalities in speech reflective of several neurological disorders. We combine this classifier with a Universal Speech Model (USM) that is trained (unsupervised) on 12 million hours of diverse audio recordings. Our model compresses long sequences into a small set of class-specific latent representations and a factorized projection is used to predict different attributes of the disordered input speech. The benefit of our approach is that it allows us to model different regions of the input for different classes and is at the same time data efficient. We evaluated the proposed model extensively on a curated corpus from the Mayo Clinic. Our model outperforms standard transformer (80.9%) and perceiver (81.8%) models and achieves an average accuracy of 83.1%. With limited task-specific data, we find that pretraining is important and surprisingly pretraining with the unrelated automatic speech recognition (ASR) task is also beneficial. Encodings from the middle layers provide a mix of both acoustic and phonetic information and achieve best prediction results compared to just using the final layer encodings (83.1% vs. 79.6%). The results are promising and with further refinements may help clinicians detect speech abnormalities without needing access to highly specialized speech-language pathologists. △ Less

Submitted 16 October, 2023; originally announced October 2023.

Journal ref: Proc. ASRU, 2023

Modeling Information Flow with a Multi-Stage Queuing Mode

Authors: Richard C. Barnard, Mohammad Daneshvar, Cory Hauck, Ilya Timofeyev

Abstract: In this paper, we introduce a nonlinear stochastic model to describe the propagation of information inside a computer processor. In this model, a computational task is divided into stages, and information can flow from one stage to another. The model is formulated as a spatially-extended, continuous-time Markov chain where space represents different stages. This model is equivalent to a spatially-… ▽ More In this paper, we introduce a nonlinear stochastic model to describe the propagation of information inside a computer processor. In this model, a computational task is divided into stages, and information can flow from one stage to another. The model is formulated as a spatially-extended, continuous-time Markov chain where space represents different stages. This model is equivalent to a spatially-extended version of the M/M/s queue. The main modeling feature is the throttling function which describes the processor slowdown when the amount of information falls below a certain threshold. We derive the stationary distribution for this stochastic model and develop a closure for a deterministic ODE system that approximates the evolution of the mean and variance of the stochastic model. We demonstrate the validity of the closure with numerical simulations. △ Less

Submitted 4 August, 2023; originally announced August 2023.

MSC Class: 60J20; 60E05

Local sensing of absolute refractive index during protein-binding using microlasers with spectral encoding

Authors: Soraya Caixeiro, Casper Kunstmann-Olsen, Marcel Schubert, Joseph Hill, Isla R. M. Barnard, Matthew D. Simmons, Steven Johnson, Malte C. Gather

Abstract: Multiplexed, specific and sensitive detection of antigens is critical for the rapid and accurate diagnosis of disease and the informed development of personalized treatment plans. Here, we show that polymer microsphere lasers can be used as photonic sensors to monitor and quantify direct surface binding of biomolecules via changes in the refractive index. The unique spectral signature of each indi… ▽ More Multiplexed, specific and sensitive detection of antigens is critical for the rapid and accurate diagnosis of disease and the informed development of personalized treatment plans. Here, we show that polymer microsphere lasers can be used as photonic sensors to monitor and quantify direct surface binding of biomolecules via changes in the refractive index. The unique spectral signature of each individual laser can be used to find their size and effective refractive index which adds a new encoding dimension when compared to conventional fluorescent beads. We utilize antibody-functionalized microlasers to selectively detect protein binding. Different stages of the multilayer surface modification can be resolved, and protein binding is demonstrated for two different proteins, IgG and CRP. Moreover, by continuously monitoring single lasers, we demonstrate the possibility of real-time monitoring of binding dynamics between antigens in solution phase and the immobilized antibodies. For multiplexed detection, the microlasers are employed in a flow cytometer configuration, with fast spectral detection and identification of microlasers with and without antigen binding. We envision that by combining microlasers with well-established surface modification chemistries and flow geometries, the multiplexing ability of microbead immunoassays can be strongly increased while also opening avenues for single cell profiling within heterogenous cell populations. △ Less

Submitted 19 October, 2022; originally announced October 2022.

Comments: 18 pages 5 figures and SI

A mathematical model of asynchronous data flow in parallel computers

Authors: Richard C Barnard, Kai Huang, Cory Hauck

Abstract: We present a simplified model of data flow on processors in a high performance computing framework involving computations necessitating inter-processor communications. From this ordinary differential model, we take its asymptotic limit, resulting in a model which treats the computer as a continuum of processors and data flow as an Eulerian fluid governed by a conservation law. We derive a Hamilton… ▽ More We present a simplified model of data flow on processors in a high performance computing framework involving computations necessitating inter-processor communications. From this ordinary differential model, we take its asymptotic limit, resulting in a model which treats the computer as a continuum of processors and data flow as an Eulerian fluid governed by a conservation law. We derive a Hamilton-Jacobi equation associated with this conservation law for which the existence and uniqueness of solutions can be proven. We then present the results of numerical experiments for both discrete and continuum models; these show a qualitative agreement between the two and the effect of variations in the computing environment's processing capabilities on the progress of the modeled computation. △ Less

Submitted 23 October, 2019; v1 submitted 8 October, 2019; originally announced October 2019.

MSC Class: 35L65; 93A30; 70H20; 41A60

Three Convolution Inequalities on the Real Line with Connections to Additive Combinatorics

Authors: Richard C. Barnard, Stefan Steinerberger

Abstract: We discuss three convolution inequalities that are connected to additive combinatorics. Cloninger and the second author showed that for nonnegative $f \in L^1(-1/4, 1/4)$, $$ \max_{-1/2 \leq t \leq 1/2} \int_{\mathbb{R}}{f(t-x) f(x) dx} \geq 1.28 \left( \int_{-1/4}^{1/4}{f(x) dx}\right)^2$$ which is related to $g-$Sidon sets (1.28 cannot be replaced by 1.52). We prove a dual statement, related to… ▽ More We discuss three convolution inequalities that are connected to additive combinatorics. Cloninger and the second author showed that for nonnegative $f \in L^1(-1/4, 1/4)$, $$ \max_{-1/2 \leq t \leq 1/2} \int_{\mathbb{R}}{f(t-x) f(x) dx} \geq 1.28 \left( \int_{-1/4}^{1/4}{f(x) dx}\right)^2$$ which is related to $g-$Sidon sets (1.28 cannot be replaced by 1.52). We prove a dual statement, related to difference bases, and show that for $f \in L^1(\mathbb{R})$, $$ \min_{0 \leq t \leq 1}\int_{\mathbb{R}}{f(x) f(x+t) dx} \leq 0.42 \|f\|_{L^1}^2,$$ where the constant 1/2 is trivial, 0.42 cannot be replaced by 0.37. This suggests a natural conjecture about the asymptotic structure of $g-$difference bases. Finally, we show for all functions $f \in L^1(\mathbb{R}) \cap L^2(\mathbb{R})$, $$ \int_{-\frac{1}{2}}^{\frac{1}{2}}{ \int_{\mathbb{R}}{f(x) f(x+t) dx}dt} \leq 0.91 \|f\|_{L^1}\|f\|_{L^2}$$ △ Less

Submitted 19 February, 2020; v1 submitted 20 March, 2019; originally announced March 2019.

Journal ref: Journal of Number Theory 207, p. 42-55 (2020)

Epidemic threshold in pairwise models for clustered networks: closures and fast correlations

Authors: Rosanna C. Barnard, Luc Berthouze, Péter L. Simon, István Z. Kiss

Abstract: The epidemic threshold is probably the most studied quantity in the modelling of epidemics on networks. For a large class of networks and dynamics the epidemic threshold is well studied and understood. However, it is less so for clustered networks where theoretical results are mostly limited to idealised networks. In this paper we focus on a class of models known as pairwise models where, to our k… ▽ More The epidemic threshold is probably the most studied quantity in the modelling of epidemics on networks. For a large class of networks and dynamics the epidemic threshold is well studied and understood. However, it is less so for clustered networks where theoretical results are mostly limited to idealised networks. In this paper we focus on a class of models known as pairwise models where, to our knowledge, no analytical result for the epidemic threshold exists. We show that by exploiting the presence of fast variables and using some standard techniques from perturbation theory we are able to obtain the epidemic threshold analytically. We validate this new threshold by comparing it to the numerical solution of the full system. The agreement is found to be excellent over a wide range of values of the clustering coefficient, transmission rate and average degree of the network. Interestingly, we find that the analytical form of $R_0$ depends on the choice of closure, highlighting the importance of model choice when dealing with real-world epidemics. Nevertheless, we ex- pect that our method will extend to other systems in which fast variables are present. △ Less

Submitted 15 June, 2018; originally announced June 2018.

A Bilevel Approach for Parameter Learning in Inverse Problems

Authors: Gernot Holler, Karl Kunisch, Richard C. Barnard

Abstract: A learning approach to selecting regularization parameters in multi-penalty Tikhonov regularization is investigated. It leads to a bilevel optimization problem, where the lower level problem is a Tikhonov regularized problem parameterized in the regularization parameters. Conditions which ensure the existence of solutions to the bilevel optimization problem of interest are derived, and these condi… ▽ More A learning approach to selecting regularization parameters in multi-penalty Tikhonov regularization is investigated. It leads to a bilevel optimization problem, where the lower level problem is a Tikhonov regularized problem parameterized in the regularization parameters. Conditions which ensure the existence of solutions to the bilevel optimization problem of interest are derived, and these conditions are verified for two relevant examples. Difficulties arising from the possible lack of convexity of the lower level problems are discussed. Optimality conditions are given provided that a reasonable constraint qualification holds. Finally, results from numerical experiments used to test the developed theory are presented. △ Less

Submitted 5 February, 2018; originally announced February 2018.

Comments: 26 pages, 16 figures

MSC Class: 49K20; 35J47; 49N45; 93C20; 90C30; 65J20; 90C31; 90C26; 90C46

Edge-based compartmental modelling of an SIR epidemic on a dual-layer static-dynamic multiplex network with tunable clustering

Authors: Rosanna C Barnard, Istvan Z Kiss, Luc Berthouze, Joel C Miller

Abstract: The duration, type and structure of connections between individuals in real-world populations play a crucial role in how diseases invade and spread. Here, we incorporate the aforementioned heterogeneities into a model by considering a dual-layer static-dynamic multiplex network. The static network layer affords tunable clustering and describes an individual's permanent community structure. The dyn… ▽ More The duration, type and structure of connections between individuals in real-world populations play a crucial role in how diseases invade and spread. Here, we incorporate the aforementioned heterogeneities into a model by considering a dual-layer static-dynamic multiplex network. The static network layer affords tunable clustering and describes an individual's permanent community structure. The dynamic network layer describes the transient connections an individual makes with members of the wider population by imposing constant edge rewiring. We follow the edge-based compartmental modelling approach to derive equations describing the evolution of a susceptible - infected - recovered (SIR) epidemic spreading through this multiplex network of individuals. We derive the basic reproduction number, measuring the expected number of new infectious cases caused by a single infectious individual in an otherwise susceptible population. We validate model equations by showing convergence to pre-existing edge-based compartmental model equations in limiting cases and by comparison with stochastically simulated epidemics. We explore the effects of altering model parameters and multiplex network attributes on resultant epidemic dynamics. We validate the basic reproduction number by plotting its value against associated final epidemic sizes measured from simulation and predicted by model equations for a number of setups. Further, we explore the effect of varying individual model parameters on the basic reproduction number. We conclude with a discussion of the significance and interpretation of the model and its relation to existing research literature. We highlight intrinsic limitations and potential extensions of the present model and outline future research considerations, both experimental and theoretical. △ Less

Submitted 4 April, 2018; v1 submitted 4 January, 2018; originally announced January 2018.

Comments: 29 pages, 12 Figures. Submitted to Bulletin of Mathematical Biology on 29/03/18. Version 2 has minor changes to the content: a reduced and edited abstract, the author summary is removed, acknowledgements were edited, references to the matrix G were bolded, some Figure titles and captions have been edited, and a reference to providing Matlab code at the end of the Methods section has been removed

MSC Class: 92C42; 90C35

Scientific Data Interpolation with Low Dimensional Manifold Model

Authors: Wei Zhu, Bao Wang, Richard Barnard, Cory D. Hauck, Frank Jenko, Stanley Osher

Abstract: We propose to apply a low dimensional manifold model to scientific data interpolation from regular and irregular samplings with a significant amount of missing information. The low dimensionality of the patch manifold for general scientific data sets has been used as a regularizer in a variational formulation. The problem is solved via alternating minimization with respect to the manifold and the… ▽ More We propose to apply a low dimensional manifold model to scientific data interpolation from regular and irregular samplings with a significant amount of missing information. The low dimensionality of the patch manifold for general scientific data sets has been used as a regularizer in a variational formulation. The problem is solved via alternating minimization with respect to the manifold and the data set, and the Laplace-Beltrami operator in the Euler-Lagrange equation is discretized using the weighted graph Laplacian. Various scientific data sets from different fields of study are used to illustrate the performance of the proposed algorithm on data compression and interpolation from both regular and irregular samplings. △ Less

Submitted 22 June, 2017; originally announced June 2017.

Accelerating optimization-based computed tomography via sparse matrix approximations

Authors: Richard C. Barnard, Rick Archibald

Abstract: Variational formulations of reconstruction in computed tomography have the notable drawback of requiring repeated evaluations of both the forward Radon transform and either its adjoint or an approximate inverse transform which are relatively expensive. We look at two methods for reducing the effect of this resulting computational bottleneck via approximating the transform evaluation with sparse ma… ▽ More Variational formulations of reconstruction in computed tomography have the notable drawback of requiring repeated evaluations of both the forward Radon transform and either its adjoint or an approximate inverse transform which are relatively expensive. We look at two methods for reducing the effect of this resulting computational bottleneck via approximating the transform evaluation with sparse matrix multiplications. The first method is applicable for general iterative optimization algorithms. The second is applicable in error-forgetting algorithms such as split Bregman. We demonstrate these approximations significantly reduce the needed computational time needed for the iterative algorithms needed to solve the reconstruction problem while still providing good reconstructions. △ Less

Submitted 21 May, 2017; originally announced May 2017.

MSC Class: 44A12; 49N45

A convex penalty for switching control of partial differential equations

Authors: Christian Clason, Armin Rund, Karl Kunisch, Richard C. Barnard

Abstract: A convex penalty for promoting switching controls for partial differential equations is introduced; such controls consist of an arbitrary number of components of which at most one should be simultaneously active. Using a Moreau-Yosida approximation, a family of approximating problems is obtained that is amenable to solution by a semismooth Newton method. The efficiency of this approach and the str… ▽ More A convex penalty for promoting switching controls for partial differential equations is introduced; such controls consist of an arbitrary number of components of which at most one should be simultaneously active. Using a Moreau-Yosida approximation, a family of approximating problems is obtained that is amenable to solution by a semismooth Newton method. The efficiency of this approach and the structure of the obtained controls are demonstrated by numerical examples. △ Less

Submitted 24 February, 2017; originally announced February 2017.

Journal ref: Systems & Control Letters 89 (2016), 66-73

L1 penalization of volumetric dose objectives in optimal control of PDEs

Authors: Richard C. Barnard, Christian Clason

Abstract: This work is concerned with a class of optimal control problems governed by a partial differential equation that are motivated by an application in radiotherapy treatment planning, where the primary design objective is to minimize the volume where a functional of the state violates a prescribed level, but prescribing these levels in the form of pointwise state constraints can lead to infeasible pr… ▽ More This work is concerned with a class of optimal control problems governed by a partial differential equation that are motivated by an application in radiotherapy treatment planning, where the primary design objective is to minimize the volume where a functional of the state violates a prescribed level, but prescribing these levels in the form of pointwise state constraints can lead to infeasible problems. We therefore propose an alternative approach based on $L^1$ penalization of the violation. We establish well-posedness of the corresponding optimal control problem, derive first-order optimality conditions, and present a semismooth Newton method for the efficient numerical solution of these problems. The performance of this method for a model problem is illustrated and contrasted with the alternative approach based on (regularized) state constraints. △ Less

Submitted 11 February, 2017; v1 submitted 6 July, 2016; originally announced July 2016.

Comments: Computational Optimization and Applications, Online First (2017)

Journal ref: Computational Optimization and Applications (67) 2017, 401-419

On the Spin of the Black Hole in IC 10 X-1

Authors: James F. Steiner, Dominic J. Walton, Javier A. Garcia, Jeffrey E. McClintock, Silas G. T. Laycock, Matthew J. Middleton, Robin Barnard, Kristin K. Madsen

Abstract: The compact X-ray source in the eclipsing X-ray binary IC 10 X-1 has reigned for years as ostensibly the most massive stellar-mass black hole, with a mass estimated to be about twice that of its closest rival. However, striking results presented recently by Laycock et al. reveal that the mass estimate, based on emission-line velocities, is unreliable and that the mass of the X-ray source is essent… ▽ More The compact X-ray source in the eclipsing X-ray binary IC 10 X-1 has reigned for years as ostensibly the most massive stellar-mass black hole, with a mass estimated to be about twice that of its closest rival. However, striking results presented recently by Laycock et al. reveal that the mass estimate, based on emission-line velocities, is unreliable and that the mass of the X-ray source is essentially unconstrained. Using Chandra and NuSTAR data, we rule against a neutron-star model and conclude that IC 10 X-1 contains a black hole. The eclipse duration of IC 10 X-1 is shorter and its depth shallower at higher energies, an effect consistent with the X-ray emission being obscured during eclipse by a Compton-thick core of a dense wind. The spectrum is strongly disk-dominated, which allows us to constrain the spin of the black hole via X-ray continuum fitting. Three other wind-fed black-hole systems are known; the masses and spins of their black holes are high: M ~ 10-15 Msun and a*>0.8. If the mass of IC 10 X-1's black hole is comparable, then its spin is likewise high. △ Less

Submitted 10 December, 2015; originally announced December 2015.

Comments: 27 manuscript pages, 7 figures, 2 tables, submitted to ApJ

Chandra, XMM-Newton and HST observations of a transient in M31 with a possible asymmetric, precessing disk

Authors: R. Barnard, M. R. Garcia, S. S. Murray

Abstract: The transient X-ray source CXOM31 004205.77+411330.43 exhibited several outbursts during our long-term monitoring campaign of ~monthly observations of the M31 center with Chandra. However, the decay profile appears to be unlike those observed from Galactic transients. We followed up the 2011 outburst with two ~B band HST/ACS observations, one in outburst and the other in quiescence, and used diffe… ▽ More The transient X-ray source CXOM31 004205.77+411330.43 exhibited several outbursts during our long-term monitoring campaign of ~monthly observations of the M31 center with Chandra. However, the decay profile appears to be unlike those observed from Galactic transients. We followed up the 2011 outburst with two ~B band HST/ACS observations, one in outburst and the other in quiescence, and used difference imaging to search for a counterpart; this would be dominated by re-processed X-ray emission from the disk. We found a counterpart with B = 28.21+/-0.16. An XMM-Newton observation from a previous outburst yielded a spectrum that is well described by an absorbed power law with absorption equivalent to ~2.6 E+21 H atom cm^(-2) and photon index ~1.8; the highest quality Chandra spectrum yielded only $\sim$130 counts, and best fits consistent with the XMM-Newton results. We calculated an absolute V magnitude of +1.9 during outburst for a typical disk spectrum. An empirical relation between the ratio of X-ray to optical flux and orbital period suggests a period < ~4 hr for a black hole accretor. Such a short period is expected to yield an asymmetric, precessing disk, and we propose that the observed decay lightcurve is due to modulation of the mass transfer rate due to the disk precession. △ Less

Submitted 11 March, 2015; originally announced March 2015.

Comments: Accepted for publication in MNRAS, 7 pages, 4 figures. arXiv admin note: text overlap with arXiv:1305.1302

XMM-Newton and Chandra observations of the M31 globular cluster black hole candidate XB135: a heavyweight contender cut down to size

Authors: R. Barnard, F. Primini, M. R. Garcia, U. C. Kolb, S. S. Murray

Abstract: CXOM31 J004252.030+413107.87 is one of the brightest X-ray sources within the D_25 region of M31, and associated with a globular cluster (GC) known as B135; we therefore call this X-ray source XB135. XB135 is a low mass X-ray binary (LMXB) that apparently exhibited hard state characteristics at 0.3--10 keV luminosities 4--6 E+38 erg/s, and the hard state is only observed below ~10% Eddington. If t… ▽ More CXOM31 J004252.030+413107.87 is one of the brightest X-ray sources within the D_25 region of M31, and associated with a globular cluster (GC) known as B135; we therefore call this X-ray source XB135. XB135 is a low mass X-ray binary (LMXB) that apparently exhibited hard state characteristics at 0.3--10 keV luminosities 4--6 E+38 erg/s, and the hard state is only observed below ~10% Eddington. If true, the accretor would be a high mass black hole (BH) (> ~50 M_Sun); such a BH may be formed from direct collapse of a metal-poor, high mass star, and the very low metalicity of B135 (0.015 Z_Sun) makes such a scenario plausible. We have obtained new XMM-Newton and Chandra HRC observations to shed light on the nature of this object. We find from the HRC observation that XB135 is a single point source located close to the center of B135. The new XMM-Newton spectrum is consistent with a rapidly spinning ~10--20 M_Sun BH in the steep power law or thermal dominant state, but inconsistent with the hard state that we previously assumed. We cannot formally reject three component emission models that have been associated with high luminosity neutron star LMXBs (known as Z-sources); however, we prefer a BH accretor. We note that deeper observation of XB135 could discriminate against a neutron star accretor. △ Less

Submitted 26 January, 2015; originally announced January 2015.

Comments: Accepted for publication in ApJ; 8 pages, 6 figures

Swift reveals a ~5.7 day super-orbital period in the M31 globular cluster X-ray binary XB158

Authors: R. Barnard, M. R. Garcia, S. S. Murrray

Abstract: The M31 globular cluster X-ray binary XB158 (a.k.a. Bo 158) exhibits intensity dips on a 2.78 hr period in some observations, but not others. The short period suggests a low mass ratio, and an asymmetric, precessing disk due to additional tidal torques from the donor star since the disk crosses the 3:1 resonance. Previous theoretical 3D smoothed particle hydrodynamical modeling suggested a super-o… ▽ More The M31 globular cluster X-ray binary XB158 (a.k.a. Bo 158) exhibits intensity dips on a 2.78 hr period in some observations, but not others. The short period suggests a low mass ratio, and an asymmetric, precessing disk due to additional tidal torques from the donor star since the disk crosses the 3:1 resonance. Previous theoretical 3D smoothed particle hydrodynamical modeling suggested a super-orbital disk precession period 29$\pm$1 times the orbital period, i.e. $\sim$81$\pm$3 hr. We conducted a Swift monitoring campaign of 30 observations over ~1 month in order to search for evidence of such a super-orbital period. Fitting the 0.3--10 keV Swift XRT luminosity lightcurve with a sinusoid yielded a period of 5.65+/-0.05 days, and a >5$σ$ improvement in $χ^2$ over the best fit constant intensity model. A Lomb-Scargle periodogram revealed that periods 5.4--5.8 days were detected at a >3$σ$ level, with a peak at 5.6 days. We consider this strong evidence for a 5.65 day super-orbital period, ~70\% longer than the predicted period. The 0.3--10 keV luminosity varied by a factor ~5, consistent with variations seen in long-term monitoring from Chandra. We conclude that other X-ray binaries exhibiting similar long-term behaviour are likely to also be X-ray binaries with low mass ratios and super-orbital periods. △ Less

Submitted 8 January, 2015; originally announced January 2015.

Comments: Accepted for publication in ApJ; 4 pages, 3 figures

Sensitivity analysis for dose deposition in radiotherapy via a Fokker-Planck model

Authors: Richard C. Barnard, Martin Frank, Kai Krycki

Abstract: In this paper we study the sensitivities of electron dose calculations with respect to the stopping power and the transport coefficients. We focus on the application to radiotherapy simulations. We use a Fokker-Planck approximation to the Boltzmann transport equation. Equations for the sensitivities are derived by the adjoint method. The Fokker-Planck equation and its adjoint are solved numericall… ▽ More In this paper we study the sensitivities of electron dose calculations with respect to the stopping power and the transport coefficients. We focus on the application to radiotherapy simulations. We use a Fokker-Planck approximation to the Boltzmann transport equation. Equations for the sensitivities are derived by the adjoint method. The Fokker-Planck equation and its adjoint are solved numerically in slab geometry using the spherical harmonics expansion ($P_N$) and an HLL finite volume method. Our method is verified by comparison to finite difference approximations of the sensitivities. Finally, we present numerical results of the sensitivities for the normalized average dose deposition depth with respect to the stopping power and transport coefficients, demonstrating the increasing relative sensitivities as beam energy decreases. △ Less

Submitted 21 October, 2014; originally announced October 2014.

MSC Class: 65M08 92C50 78A70

Energy-dependent evolution in IC10 X-1: hard evidence for an extended corona, and implications

Authors: R. Barnard, J. F. Steiner, A. F. Prestwich, I. R. Stevens, J. S. Clark, U. C. Kolb

Abstract: We have analyzed a ~130 ks XMM-Newton observation of the dynamically confirmed black hole + Wolf-Rayet (BH+WR) X-ray binary (XB) IC10 X-1, covering ~1 orbital cycle. This system experiences periodic intensity dips every ~35 hours. We find that energy-independent evolution is rejected at a >5 sigma level. The spectral and timing evolution of IC10 X-1 are best explained by a compact disk blackbody a… ▽ More We have analyzed a ~130 ks XMM-Newton observation of the dynamically confirmed black hole + Wolf-Rayet (BH+WR) X-ray binary (XB) IC10 X-1, covering ~1 orbital cycle. This system experiences periodic intensity dips every ~35 hours. We find that energy-independent evolution is rejected at a >5 sigma level. The spectral and timing evolution of IC10 X-1 are best explained by a compact disk blackbody and an extended Comptonized component, where the thermal component is completely absorbed and the Comptonized component is partially covered during the dip. We consider three possibilities for the absorber: cold material in the outer accretion disk, as is well documented for Galactic neutron star (NS) XBs at high inclination; a stream of stellar wind that is enhanced by traveling through the L1 point; and a spherical wind. We estimated the corona radius (r_ADC) for IC10 X-1 from the dip ingress to be ~1 E+6 km, assuming absorption from the outer disk, and found it to be consistent with the relation between r_m ADC and 1--30 keV luminosity observed in Galactic NS XBs that spans 2 orders of magnitude. For the other two scenarios, the corona would be larger. Prior BH mass (M_BH) estimates range over 23--38 M_Sun, depending on the inclination and WR mass. For disk absorption, the inclination, i, is likely to be ~60--80 degrees, with M_BH ~24--41 M_Sun. Alternatively, the L1-enhanced wind requires i ~80 degrees, suggesting ~24--33 M_Sun. For a spherical absorber, i ~ 40 degrees, and M_BH ~50--65 M_Sun. △ Less

Submitted 21 July, 2014; originally announced July 2014.

Comments: Accepted for publication in ApJ; 10 pages, 4 figures

50 M31 black hole candidates identified by Chandra and XMM-Newton

Authors: R. Barnard, M. R. Garcia, F. Primini, S. S. Murray

Abstract: Over the last ~5 years we have identified ~35 black hole candidates (BHCs) in M31 from their X-ray spectra. Our BHCs exhibited 0.3--10 keV spectra consistent with the X-ray binary (XB) hard state, at luminosities that are above the upper limit for neutron star (NS) XBs. When our BHC spectra were modeled with a disk blackbody + blackbody model, for comparison with bright NS XBs, we found that the B… ▽ More Over the last ~5 years we have identified ~35 black hole candidates (BHCs) in M31 from their X-ray spectra. Our BHCs exhibited 0.3--10 keV spectra consistent with the X-ray binary (XB) hard state, at luminosities that are above the upper limit for neutron star (NS) XBs. When our BHC spectra were modeled with a disk blackbody + blackbody model, for comparison with bright NS XBs, we found that the BHCs inhabited a different parameter space to the NS XBs. However, BH XBs may also exhibit a thermally dominated (TD) state that has never been seen in NS XBs; this TD state is most often observed in X-ray transients. We examined the ~50 X-ray transients in our Chandra survey of M31, and found 13 with spectra suitable for analysis. We also examined 2 BHCs outside the field of view of our survey, in the globular clusters B045 and B375. We have 42 strong BHCs, and 8 plausible BHCs that may benefit from further observation. Of our 15 BHCs in globular clusters, 12 differ from NS spectra by >5 sigma. Due to improvements in our analysis, we have upgraded 10 previously identified plausible BHCs to strong BHCs. The mean maximum duty cycle of the 33 X-ray transients within 6' of M31* was 0.13; we estimate that >40% of the XBs in this region contain BH accretors. Remarkably, we estimate that BHCs contribute >90% of those XBs >1 E+38 erg/s. △ Less

Submitted 23 June, 2014; originally announced June 2014.

Comments: Accepted for publication in ApJ; 8 pages, 1 figure

Discovery of a new black hole transient within 100" (400 pc) of M31*

Authors: R. Barnard, M. R. Garcia, F. Primini, S. S. Murray

Abstract: We identified a new X-ray transient CXOM31 004252.457+411631.17 (T13) in M31 during a 2013 June Chandra observation. This system is particularly exciting because it is located within 100" of M31*; it is thought that this region of the bulge is sufficiently dense to form X-ray binaries dynamicaly, but only systems with black hole accretors and/or short periods are expected to survive. A follow-up X… ▽ More We identified a new X-ray transient CXOM31 004252.457+411631.17 (T13) in M31 during a 2013 June Chandra observation. This system is particularly exciting because it is located within 100" of M31*; it is thought that this region of the bulge is sufficiently dense to form X-ray binaries dynamicaly, but only systems with black hole accretors and/or short periods are expected to survive. A follow-up XMM-Newton observation yielded a soft spectrum, well described by a 0.39+/-0.02 keV disk blackbody; applying this model to the Chandra observation yielded an observed 0.3--10 keV luminosity peak of 6.2+/-0.6 E+37 erg/s (4.7 E+36 erg/s in the 2.0--10 keV band). Observing with HST/ACS did not reveal an optical counterpart, but allowed us to place an upper limit of B >26.9, corresponding to an absolute V band magnitude >2.0. From the 2--10 keV luminosity and absolute V magnitude, we estimate an orbital period $<$5 h from an empirical relation. Fitting a disk blackbody + blackbody model allows us to reject a neutron star accretor at a 14 sigma level. △ Less

Submitted 4 December, 2013; originally announced December 2013.

Comments: Accepted for publication in ApJ. 4 pages, 4 figures. arXiv admin note: substantial text overlap with arXiv:1305.1302

Around 200 new X-ray binary IDs from 13 years of Chandra observations of the M31 center

Authors: R. Barnard, M. R. Garcia, F. Primini, Z. Li, F. Baganoff, S. S. Murray

Abstract: We have created 0.3--10 keV, 13 year, unabsorbed luminosity lightcurves for 528 X-ray sources in the central 20' of M31. We have 174 Chandra observations spaced at ~1 month intervals thanks to our transient monitoring program, deeper observations of the M31 nucleus, and some public data from other surveys. We created 0.5--4.5 keV structure functions (SFs) for each source, for comparison with the e… ▽ More We have created 0.3--10 keV, 13 year, unabsorbed luminosity lightcurves for 528 X-ray sources in the central 20' of M31. We have 174 Chandra observations spaced at ~1 month intervals thanks to our transient monitoring program, deeper observations of the M31 nucleus, and some public data from other surveys. We created 0.5--4.5 keV structure functions (SFs) for each source, for comparison with the ensemble structure function of AGN. We find 220 X-ray sources with luminosities > ~1E+35 erg/s that have SFs with significantly more variability than the ensemble AGN SF, and are likely X-ray binaries (XBs). A further 30 X-ray sources were identified as XBs using other methods. We therefore have 250 probable XBs in total, including ~200 new identifications. This result represents great progress over the ~50 XBs and ~40 XB candidates previously identified out of the ~2000 X-ray sources within the D_25 region of M31; it also demonstrates the power of SF analysis for identifying XBs in external galaxies. We also identify a new transient black hole candidate, associated with the M31 globular cluster B128. △ Less

Submitted 7 November, 2013; originally announced November 2013.

Comments: Accepted for publication in ApJ; 20 pages, 7 figures; Tables 1 and 2 continue after the references (8 pages)

The second ULX transient in M31: Chandra, HST and XMM observations, and evidence for an extended corona

Authors: R. Barnard, M. Garcia, S. S. Murray

Abstract: XMM J004243.6+412519 is a transient X-ray source in M31, first discovered 2012 January 15. Different approaches to fitting the brightest follow-up observation gave luminosities 1.3--2.5E+39 erg/s, making it the second ultraluminous X-ray source (ULX) in M31, with a probable black hole accretor. These different models represent different scenarios for the corona: optically thick and compact, or opt… ▽ More XMM J004243.6+412519 is a transient X-ray source in M31, first discovered 2012 January 15. Different approaches to fitting the brightest follow-up observation gave luminosities 1.3--2.5E+39 erg/s, making it the second ultraluminous X-ray source (ULX) in M31, with a probable black hole accretor. These different models represent different scenarios for the corona: optically thick and compact, or optically thin and extended. We obtained Chandra ACIS and {\em HST} ACS observations of this object as part of our transient monitoring program, and also observed it serendipitously in a 120 ks XMM-Newton observation. We identify an optical counterpart at J2000 position 00:42:43.70 +41:25:18.54; its F435W (~ B band) magnitude was 25.97+/-0.03 in the 2012 March 7 observation, and >28.4 at the 4 sigma level during the 2012 September 7 observation, indicating a low mass donor. We created two alternative lightcurves, using the different corona scenarios, finding linear decay for the compact corona and exponential decay for the extended corona; linear decay implies a disk that is >5 magnitudes brighter than we observed. We therefore favor the extended corona scenario, but caution that there is no statistical preference for this model in the X-ray spectra alone. Using two empirical relations between the X-ray to optical ratio and the orbital period, we estimate a period of ~9--30 hr; this period is consistent with that of the first ULX in M31 (18 +5 -6 hr). △ Less

Submitted 6 May, 2013; originally announced May 2013.

Comments: Accepted for publication in ApJ; 9 pages, 4 figures

Chandra identification of 26 new black hole candidates in the central region of M31

Authors: R. Barnard, M. R. Garcia, S. S. Murray

Abstract: We have previously identified 10 M31 black hole candidates (BHCs) in M31, from their X-ray properties alone. They exhibit "hard state" emission spectra that are seen at luminosities < ~10% Eddington in X-ray binaries (XBs) containing a neutron star (NS) or black hole (BH), at luminosities that significantly exceed the NS threshold. Nine of these are associated with globular clusters (GCs); hence,… ▽ More We have previously identified 10 M31 black hole candidates (BHCs) in M31, from their X-ray properties alone. They exhibit "hard state" emission spectra that are seen at luminosities < ~10% Eddington in X-ray binaries (XBs) containing a neutron star (NS) or black hole (BH), at luminosities that significantly exceed the NS threshold. Nine of these are associated with globular clusters (GCs); hence, these are most likely low mass X-ray binaries (LMXBs); eight are included in this survey. We have recently discovered that analysis of the long term 0.5--4.5 keV variability of XBs via structure functions allows us to separate XBs from AGN, even though the emission spectra are often similar; this has enabled us to search for BHCs outside of GCs. We have identified 26 new BHCs (12 strong, 14 plausible) within 20' of the M31 nucleus (M31*), using 152 Chandra observations spaced over ~13 years; some of our classifications were enhanced with XMM-Newton observations. Of these, 7 appear within 100" of M31*; this supports the theory suggesting that this region experiences enhanced XB production via dynamical processes similar to those seen in GCs. We have found a parameter space where our black hole candidates are separated from Galactic neutron star binaries: we show that modelling a simulated hard state spectrum with a disk blackbody + blackbody model yields parameters that lie outside the space occupied by neutron star binaries that are modeled this way. The probability that our BHCs all lie within the NS parameter space is ~3E-29. △ Less

Submitted 29 April, 2013; originally announced April 2013.

Comments: Accepted for publication in ApJ, 22 pages, 17 figures, 4 pages. Figures 3 and 5 have multiple pages (6 figures each), and are continued after the references

Flow Invariance on Stratified Domains

Authors: Richard Barnard, Peter Wolenski

Abstract: This paper studies conditions for invariance of dynamical systems on stratified do- mains as originally introduced by Bressan and Hong. We establish Hamiltonian conditions for both weak and strong invariance of trajectories on systems with non-Lipschitz data. This is done via the identification of a new multifunction, the essential velocity multifunction. Properties of this multifunction are inves… ▽ More This paper studies conditions for invariance of dynamical systems on stratified do- mains as originally introduced by Bressan and Hong. We establish Hamiltonian conditions for both weak and strong invariance of trajectories on systems with non-Lipschitz data. This is done via the identification of a new multifunction, the essential velocity multifunction. Properties of this multifunction are investigated and used to establish the relevant invariance criteria. △ Less

Submitted 23 August, 2012; originally announced August 2012.

MSC Class: 34H05; 49J52

State Constrained Optimization with Partial Differential Equations via Generalized Gradients

Authors: Richard Barnard

Abstract: We consider optimization problems constrained by partial differential equations (PDEs) with additional constraints placed on the solution of the PDEs. We develop a general and versatile framework using infinite-valued penalization functions and Clarke subgradients and apply this to problems with box constraints as well as more general constraints arising in applications, such as constraints on the… ▽ More We consider optimization problems constrained by partial differential equations (PDEs) with additional constraints placed on the solution of the PDEs. We develop a general and versatile framework using infinite-valued penalization functions and Clarke subgradients and apply this to problems with box constraints as well as more general constraints arising in applications, such as constraints on the average value of the state in subdomains. The framework also allows for problems with discontinuous data in the constraints. We present numerical results of this algorithm for the elliptic case and compare with other state-constrained algorithms. △ Less

Submitted 27 February, 2013; v1 submitted 23 July, 2012; originally announced July 2012.

Comments: 19 pages, 3 figures

MSC Class: 49K20; 49J52; 49J53

A period distribution of X-ray binaries observed in the central region of M31 with Chandra and HST

Authors: R. Barnard, J. L. Galache, M. Garcia, N. Nooraee, P. J. Callanan, A. Zezas, S. S. Murray

Abstract: Almost all Galactic black hole binaries with low mass donor stars are transient X-ray sources; we expect most of the X-ray transients observed in external galaxies to be black hole binaries also. Obtaining period estimates for extra-galactic transients is challenging, but the resulting period distribution is an important tool for modeling the evolution history of the host galaxy. We have obtained… ▽ More Almost all Galactic black hole binaries with low mass donor stars are transient X-ray sources; we expect most of the X-ray transients observed in external galaxies to be black hole binaries also. Obtaining period estimates for extra-galactic transients is challenging, but the resulting period distribution is an important tool for modeling the evolution history of the host galaxy. We have obtained periods, or upper limits, for 12 transients in M31, using an updated relation between the optical and X-ray luminosities. We have monitored the central region of M31 with Chandra for the last ~12 years, and followed up promising transients with HST; 4σB magnitude limits for optical counterparts are ~26--29, depending on crowding. We obtain period estimates for each transient for both neutron star and black hole accretors. Periods range from <0.4 to 490+/-90 hours (<0.97 to <175 hrs if all are BH systems). These M31 transients appear to be somewhat skewed towards shorter periods than the Milky Way (MW) transients; indeed, comparing the M31 and MW transients with survival analysis techniques used to account for some data with only upper limits yield probabilities of ~0.02--0.08 that the two populations are drawn from the same distribution. We also checked for a correlation between orbital period and distance from the nucleus, finding a 12% probability of no correlation. Further observations of M31 transients will strengthen these results. △ Less

Submitted 25 July, 2012; v1 submitted 14 March, 2012; originally announced March 2012.

Comments: Accepted for publication in ApJ, 20 pages, 3 tables, 6 figures

12 years of X-ray variability in M31 globular clusters, including 8 black hole candidates, as seen by Chandra

Authors: R. Barnard, M. Garcia, S. S. Murray

Abstract: We examined 134 Chandra observations of the population of X-ray sources associated with globular clusters (GCs) in the central region of M31. These are expected to be X-ray binary systems (XBs), consisting of a neutron star or black hole accreting material from a close companion. We created long-term lightcurves for these sources, correcting for background, interstellar absorption and instrumental… ▽ More We examined 134 Chandra observations of the population of X-ray sources associated with globular clusters (GCs) in the central region of M31. These are expected to be X-ray binary systems (XBs), consisting of a neutron star or black hole accreting material from a close companion. We created long-term lightcurves for these sources, correcting for background, interstellar absorption and instrumental effects. We tested for variability by examining the goodness of fit for the best fit constant intensity. We also created structure functions (SFs) for every object in our sample, the first time this technique has been applied to XBs. We found significant variability in 28 out of 34 GCs and GC candidates; the other 6 sources had 0.3--10 keV luminosities fainter than ~2E+36 erg/s, limiting our ability to detect similar variability. The SFs of XBs with 0.3--10 keV luminosities ~2--50 E+36 generally showed considerably more variability than the published ensemble SF of AGN. Our brightest XBs were mostly consistent with the AGN SF; however, their 2--10 keV fluxes could be matched by <1 AGN per square degree. These encouraging results suggest that examining the long term lightcurves of other X-ray sources in the field may provide an important distinction between X-ray binaries and background galaxies, as the X-ray emission spectra from these two classes of X-ray sources are similar. Additionally, we identify 3 new black hole candidates (BHCs) using additional XMM-Newton data, bringing the total number of M31 GC BHCs to 9, with 8 covered in this survey. △ Less

Submitted 25 July, 2012; v1 submitted 12 March, 2012; originally announced March 2012.

Comments: Accepted for publication in ApJ; 30 pages, 16 figures, 3 tables

Reinstating the M31 X-ray system RX J0042.3+4115 as a black hole X-ray binary, and compelling evidence for an extended corona

Authors: R. Barnard, M. R. Garcia, S. S. Murray

Abstract: The M31 X-ray source RX J0042.3+4115 was originally identified as a black hole binary because it displayed characteristic low state variability at conspicuously high luminosities; unfortunately, this variability was later found to be artificial. However, analysis of 84 Chandra ACIS observations, an HST ACS/WFC observation, and a 60 ks XMM-Newton observation has supplied new evidence that RX J0042.… ▽ More The M31 X-ray source RX J0042.3+4115 was originally identified as a black hole binary because it displayed characteristic low state variability at conspicuously high luminosities; unfortunately, this variability was later found to be artificial. However, analysis of 84 Chandra ACIS observations, an HST ACS/WFC observation, and a 60 ks XMM-Newton observation has supplied new evidence that RX J0042.3+4115 is indeed a black hole binary. The brightest optical star within 3 sigma of the position of RX J0042.3+4115 had a F435W (~B) magnitude of 25.4+/-0.2; M_B > -0.4, hence we find a low mass donor likely. RX J0042.3+4115 was persistently bright over ~12 years. Spectral fits revealed characteristic black hole binary states: a low/hard state at 2.08+/-0.08 E+38 erg/s, and a steep power law state at 2.41+/-0.05 E+38 erg/s (0.3--10 keV). The high luminosity low state suggests a ~20 M_Sun primary; this is high, but within the range of known stellar black hole masses. The inner disk temperature during the steep power law state is 2.24$\pm$0.15 keV, high but strikingly similar to that of GRS 1915+105, the only known Galactic black hole binary with a low mass donor t o be persistently bright. Therefore RX J0042.3+4115 may be an analog for GRS 1915+105; however, other mechanisms may account for its behavior. We find compelling evidence for an extended corona during the steep power law state, because compact corona models where the seed photons for Comptonization are tied to the inner disc temperature are rejected. △ Less

Submitted 17 November, 2011; originally announced November 2011.

Comments: Accepted for publication in ApJ, 8 pages, 5 figures

The deep XMM-Newton Survey of M 31

Authors: H. Stiele, W. Pietsch, F. Haberl, D. Hatzidimitriou, R. Barnard, B. F. Williams, A. K. H. Kong, U. Kolb

Abstract: Aims. M 31 has been completely imaged for the first time obtaining a luminosity lower limit ~10^35 erg s-1 in the 0.2-4.5 keV band. Our XMM-Newton EPIC survey combines archival observations with observations taken between June 2006 and February 2008. The main goal of the paper is a study of the X-ray source population of M 31. Methods. An X-ray catalogue of 1897 sources was created; 914 were detec… ▽ More Aims. M 31 has been completely imaged for the first time obtaining a luminosity lower limit ~10^35 erg s-1 in the 0.2-4.5 keV band. Our XMM-Newton EPIC survey combines archival observations with observations taken between June 2006 and February 2008. The main goal of the paper is a study of the X-ray source population of M 31. Methods. An X-ray catalogue of 1897 sources was created; 914 were detected for the first time. Source classification and identification were based on X-ray hardness ratios, spatial extent of the sources, and by cross correlating with catalogues in the X-ray, optical, infrared and radio wavelengths. We also analysed the long-term variability of the X-ray sources. Results. Fourteen of the 30 supersoft source candidates represent supersoft emission of optical novae. Many of the 25 supernova remnants (SNRs) and 31 SNR candidates lie within the 10 kpc dust ring and other star forming regions in M 31; this connection between SNRs and star forming regions implies that most of the remnants originate in type II supernovae. The brightest sources in X-rays in M 31 belong to the class of X-ray binaries (XRBs). Ten low mass XRBs (LMXBs) and 26 LMXB candidates were identified based on their temporal variability. In addition 36 LMXBs and 17 LMXB candidates were identified due to correlations with globular clusters and globular cluster candidates. From optical and X-ray colour-colour diagrams, possible high mass XRB (HMXB) candidates were selected. Two of these candidates have an X-ray spectrum as is expected for an HMXB containing a neutron star primary. Conclusions. While our survey has greatly improved our understanding of the X-ray source populations in M 31, at this point 65% of the sources can still only be classified as "hard" sources. Deeper observations in X-ray and at other wavelengths would help classify these sources. ABRIDGED △ Less

Submitted 23 June, 2011; originally announced June 2011.

Comments: 51 pages, 36 figures, submitted to A&A

Optimal radiotherapy treatment planning using minimum entropy models

Authors: Richard Barnard, Martin Frank, Michael Herty

Abstract: We study the problem of finding an optimal radiotherapy treatment plan. A time-dependent Boltzmann particle transport model is used to model the interaction between radiative particles with tissue. This model allows for the modeling of inhomogeneities in the body and allows for anisotropic sources modeling distributed radiation---as in brachytherapy---and external beam sources---as in teletherapy.… ▽ More We study the problem of finding an optimal radiotherapy treatment plan. A time-dependent Boltzmann particle transport model is used to model the interaction between radiative particles with tissue. This model allows for the modeling of inhomogeneities in the body and allows for anisotropic sources modeling distributed radiation---as in brachytherapy---and external beam sources---as in teletherapy. We study two optimization problems: minimizing the deviation from a spatially-dependent prescribed dose through a quadratic tracking functional; and minimizing the survival of tumor cells through the use of the linear-quadratic model of radiobiological cell response. For each problem, we derive the optimality systems. In order to solve the state and adjoint equations, we use the minimum entropy approximation; the advantages of this method are discussed. Numerical results are then presented. △ Less

Submitted 26 May, 2011; originally announced May 2011.

Comments: 16 pages,3 figures

MSC Class: 92C50; 49K20; 35Q20

From X-ray dips to eclipse: Witnessing disk reformation in the recurrent nova USco

Authors: J. -U. Ness, B. E. Schaefer, A. Dobrotka, A. Sadowski, J. J. Drake, R. Barnard, A. Talavera, R. Gonzalez-Riestra, K. L. Page, M. Hernanz, G. Sala, S. Starrfield

Abstract: The 10th recorded outburst of the recurrent eclipsing nova USco was observed simultaneously in X-ray, UV, and optical by XMM-Newton on days 22.9 and 34.9 after outburst. Two full passages of the companion in front of the nova ejecta were observed, witnessing the reformation of the accretion disk. On day 22.9, we observed smooth eclipses in UV and optical but deep dips in the X-ray light curve whic… ▽ More The 10th recorded outburst of the recurrent eclipsing nova USco was observed simultaneously in X-ray, UV, and optical by XMM-Newton on days 22.9 and 34.9 after outburst. Two full passages of the companion in front of the nova ejecta were observed, witnessing the reformation of the accretion disk. On day 22.9, we observed smooth eclipses in UV and optical but deep dips in the X-ray light curve which disappeared by day 34.9, then yielding clean eclipses in all bands. X-ray dips can be caused by clumpy absorbing material that intersects the line of sight while moving along highly elliptical trajectories. Cold material from the companion could explain the absence of dips in UV and optical light. The disappearance of X-ray dips before day 34.9 implies significant progress in the formation of the disk. The X-ray spectra contain photospheric continuum emission plus strong emission lines, but no clear absorption lines. Both continuum and emission lines in the X-ray spectra indicate a temperature increase from day 22.9 to day 34.9. We find clear evidence in the spectra and light curves for Thompson scattering of the photospheric emission from the white dwarf. Photospheric absorption lines can be smeared out during scattering in a plasma of fast electrons. We also find spectral signatures of resonant line scattering that lead to the observation of the strong emission lines. Their dominance could be a general phenomenon in high-inclination systems such as Cal87. △ Less

Submitted 6 October, 2011; v1 submitted 13 May, 2011; originally announced May 2011.

Comments: Submitted to ApJ. 16 pages, 16 figures

Four new black hole candidates identified in M31 globular clusters with Chandra and XMM-Newton

Authors: R. Barnard, M. Garcia, Z. Li, F. Primini, S. S. Murray

Abstract: We have identified four new black hole candidates in M31 globular clusters using 123 Chandra, and 4 XMM-Newton observations of the M31 central region. The X-ray source associated with Bo 163 (XB163) is a recurrent transient, with the highest luminosity ~1.4E+38 erg/s, considerably brighter than any outbursts from the neutron star transients Aql X-1 or 4U 1608-452; the outburst apparently started ~… ▽ More We have identified four new black hole candidates in M31 globular clusters using 123 Chandra, and 4 XMM-Newton observations of the M31 central region. The X-ray source associated with Bo 163 (XB163) is a recurrent transient, with the highest luminosity ~1.4E+38 erg/s, considerably brighter than any outbursts from the neutron star transients Aql X-1 or 4U 1608-452; the outburst apparently started ~45 days earlier than the observed peak, hence the luminosity could have been considerably higher. We identified XB082, XB153 and XB185 as BHCs by observing "low state" emission spectra at luminosities that exceed the threshold for neutron star binaries. The probability that these are neutron star systems with anisotropic emission beamed toward us is <4E-4, and their variability suggests emission from a single source. We therefore conclude that these systems likely contain black holes rather than neutron stars. We have now identified 4 persistently bright BHCs in the region; the probability that these are all background AGN is <1E-20. According to theory, the donors could be tidally captured main sequence stars, or white dwarves in ultra-compact binaries. We find that GCs that are particularly massive (XB082) or metal rich (XB144) can host bright X-ray sources in addition to those that are both (XB163). Our method may reveal BHCs in other bright X-ray sources. △ Less

Submitted 5 April, 2011; originally announced April 2011.

Comments: Accepted for publication in ApJ. 17 pages, 5 figures

Observations of the recurrent M31 transient XMMU~J004215.8+411924 with Swift, Chandra, HST and Einstein

Authors: R. Barnard, M. Garcia, S. Murray, N. Nooraee, W. Pietsch

Abstract: The transient X-ray source XMMU J004215.8+411924 within M31 was found to be in outburst again in the 2010 May 27 Chandra observation. We present results from our four Chandra and seven Swift observations that covered this outburst. X-ray transient behaviour is generally caused by one of two things: mass accretion from a high mass companion during some restricted phase range in the orbital cycle, o… ▽ More The transient X-ray source XMMU J004215.8+411924 within M31 was found to be in outburst again in the 2010 May 27 Chandra observation. We present results from our four Chandra and seven Swift observations that covered this outburst. X-ray transient behaviour is generally caused by one of two things: mass accretion from a high mass companion during some restricted phase range in the orbital cycle, or disc instability in a low mass system. We aim to exploit Einstein, HST, Chandra and Swift observations to determine the nature of XMMU J004215.8+411924. We model the 2010 May spectrum, and use the results to convert from intensity to counts in the fainter Chandra observations, as well as the Swift observations; these data are used to create a lightcurve. We also estimate the flux in the 1979 January 13 Einstein observation. Additionally, we search for an optical counterpart in HST data. Our best X-ray positions from the 2006 and 2010 outbursts are 0.3" apart, and 1.6" from the Einstein source; these outbursts are likely to come from the same star system. We see no evidence for an optical counterpart with m_B < ~25.5; this new limit is 3.5 magnitudes fainter than the existing one. Furthermore, we see no V band counterpart with m_V < ~26. The local absorption is ~7 times higher than the Galactic line-of-sight, and provides ~2 magnitudes of extinction in the V band. Hence M_V > ~ -0.5. Fits to the X-ray emission spectrum suggest a black hole primary. We find that XMMU J004215.8+411924 is most likely to be a transient LMXB, rather than a HMXB as originaly proposed. The nature of the primary is unclear, although we argue that a black hole is likely. △ Less

Submitted 16 November, 2010; originally announced November 2010.

Comments: Accepted for publication in A&A, 5 pages, 3 figures

NGC 300 X-1 is a Wolf-Rayet/Black-Hole binary

Authors: P. A. Crowther, R. Barnard, S. Carpano, J. S. Clark, V. S. Dhillon, A. M. T. Pollock

Abstract: We present VLT/FORS2 time-series spectroscopy of the Wolf-Rayet star #41 in the Sculptor group galaxy NGC 300. We confirm a physical association with NGC 300 X-1, since radial velocity variations of the HeII 4686 line indicate an orbital period of 32.3 +/- 0.2 hr which agrees at the 2 sigma level with the X-ray period from Carpano et al. We measure a radial velocity semi-amplitude of 267 +/- 8 k… ▽ More We present VLT/FORS2 time-series spectroscopy of the Wolf-Rayet star #41 in the Sculptor group galaxy NGC 300. We confirm a physical association with NGC 300 X-1, since radial velocity variations of the HeII 4686 line indicate an orbital period of 32.3 +/- 0.2 hr which agrees at the 2 sigma level with the X-ray period from Carpano et al. We measure a radial velocity semi-amplitude of 267 +/- 8 km/s, from which a mass function of 2.6 +/- 0.3 Msun is obtained. A revised spectroscopic mass for the WN-type companion of 26+7-5 Msun yields a black hole mass of 20 +/- 4 Msun for a preferred inclination of 60-75 deg. If the WR star provides half of the measured visual continuum flux, a reduced WR (black hole) mass of 15 +4 -2.5 Msun (14.5 +3 -2.5 Msun) would be inferred. As such, #41/NGC 300 X-1 represents only the second extragalactic Wolf-Rayet plus black-hole binary system, after IC 10 X-1. In addition, the compact object responsible for NGC 300 X-1 is the second highest stellar-mass black hole known to date, exceeded only by IC 10 X-1. △ Less

Submitted 26 January, 2010; originally announced January 2010.

Comments: 5 pages, 6 figures, MNRAS Letters in press. Images and animations available from http://www.eso.org/public/news/eso1004/ at 11am GMT on 27 Jan 10

In-depth studies of the NGC 253 ULXs with XMM-Newton: remarkable variability in ULX1, and evidence for extended coronae

Authors: R. Barnard

Abstract: We examined the variability of three ultra-luminous X-ray sources (ULXs) in the 2003, 110 ks XMM-Newton observation of NGC253. Remarkably, we discovered ULX1 to be three times more variable than ULX2 in the 0.3--10 keV band, even though ULX2 is brighter. Indeed, ULX1 exhibits a power density spectrum that is consistent with the canonical high state or very high/steep power law state, but not the… ▽ More We examined the variability of three ultra-luminous X-ray sources (ULXs) in the 2003, 110 ks XMM-Newton observation of NGC253. Remarkably, we discovered ULX1 to be three times more variable than ULX2 in the 0.3--10 keV band, even though ULX2 is brighter. Indeed, ULX1 exhibits a power density spectrum that is consistent with the canonical high state or very high/steep power law state, but not the canonical low state. The 0.3--10 keV emission of ULX1 is predominantly non-thermal, and may be related to the very high state. We also fitted the ULX spectra with disc blackbody, slim disc and convolution Comptonization (SIMPL x DISKBB) models. The brightest ULX spectra are usually described by a two emission components (disc blackbody + Comptonized component); however, the SIMPL model results in a single emission component, and may help determine whether the well known soft excess is a feature of ULX spectra or an artifact of the two-component model. The SIMPL models were rejected for ULX3 (and also for the black hole + Wolf-Rayet binary IC10 X-1); hence, we infer that the observed soft-excesses are genuine features of ULX emission spectra. We use an extended corona scenario to explain the soft excess seen in all the highest quality ULX spectra, and provide a mechanism for stellar mass black holes to exhibit super-Eddington luminosities while remaining locally sub-Eddington. △ Less

Submitted 6 January, 2010; originally announced January 2010.

Comments: Accepted for publication in MNRAS. 7 pages, 5 figures

A second black hole candidate in a M31 globular cluster is identified with XMM-Newton

Authors: R. Barnard, U. Kolb

Abstract: We use arguments developed in previous work to identify a second black hole candidate associated with a M31 globular cluster, Bo 144, on the basis of X-ray spectral and timing properties. The 2002 XMM-Newton observation of the associated X-ray source (hereafter XBo 144) revealed behaviour that is common to all low-mass X-ray binaries (LMXBs) in the low-hard state. Studies have shown that neutron… ▽ More We use arguments developed in previous work to identify a second black hole candidate associated with a M31 globular cluster, Bo 144, on the basis of X-ray spectral and timing properties. The 2002 XMM-Newton observation of the associated X-ray source (hereafter XBo 144) revealed behaviour that is common to all low-mass X-ray binaries (LMXBs) in the low-hard state. Studies have shown that neutron star LMXBs exhibit this behaviour at 0.01-1000 keV luminosities <=10% of the Eddington limit (L_Edd). However, the unabsorbed 0.3-10 keV XBo 144 luminosity was ~0.30 L_Edd for a 1.4 M_sun neutron star, and the expected 0.01-1000 keV luminosity is 3-7 times higher. We therefore identify XBo 144 as a black hole candidate. Furthermore, it is the second black hole candidate to be consistent with formation via tidal capture of a mean sequence donor in a GC; such systems were previously though non-existent, because the donor was thought to be disrupted during the capture process. △ Less

Submitted 20 May, 2009; originally announced May 2009.

Comments: Accepted for publication in MRAS letters. Four pages, three figures

NGC300 X-1 and IC10 X-1: a new breed of black hole binary?

Authors: R. Barnard, J. S. Clark, U. C. Kolb

Abstract: [ABRIDGED] IC10 X-1 has recently been confirmed as a black hole (BH) + Wolf-Rayet (WR) X-ray binary, and NGC300 X-1 is thought to be. IC10 X-1 and NGC300 X-1 have similar X-ray properties, with luminosities ~10^38 erg/s, and orbital periods ~30 hr. We investigate similarities between these two, as well as differences between them and the known Galactic BH binary systems. We have examined XMM-New… ▽ More [ABRIDGED] IC10 X-1 has recently been confirmed as a black hole (BH) + Wolf-Rayet (WR) X-ray binary, and NGC300 X-1 is thought to be. IC10 X-1 and NGC300 X-1 have similar X-ray properties, with luminosities ~10^38 erg/s, and orbital periods ~30 hr. We investigate similarities between these two, as well as differences between them and the known Galactic BH binary systems. We have examined XMM-Newton observations of NGC300 X-1 and IC10 X-1. We extracted lightcurves and spectra; power density spectra (PDS) were constructed from the lightcurves, and the X-ray emission spectra were modeled. Each source exhibits PDS that are characteristic of disc-accreting X-ray binaries (XBs) in the high state. In this state, Galactic XBs with known BH primaries have soft, thermal emission; however the emission spectra of our targets are predominantly non-thermal. Furthermore, the Observation 1 spectrum of NGC300 X-1 is strikingly similar to that of IC10 X-1. The remarkable similarity between the behaviour of NGC300 X-1 in Observation 1 and that of IC10 X-1 lends strong evidence for NGC300 X-1 being a (BH+WR) binary. The unusual spectra of NGC300 X-1 and IC10 X-1 may be due to these systems existing in a persistently high state, whereas all known BH LMXBs are transient. BH XBs in a persistent high state could retain their corona, and hence exhibit a large non-thermal component. LMC X-1 is a BH XB that has only been observed in the high state, and its spectrum is remarkably similar to those of our targets. We therefore classify NGC300 X-1, IC10 X-1 and perhaps LMC X-1 as a new breed of BH XB, defined by their persistently high accretion rates and consequent stable disc configuration and corona. This scenario may also explain the lack of ultraluminous X-ray sources in the canonical soft state. △ Less

Submitted 3 July, 2008; originally announced July 2008.

Comments: Approved for publication in A&A. 8 pages, 5 figures

New XMM-Newton analysis of three bright X-ray sources in M31 globular clusters, including a new black hole candidate

Authors: R. Barnard, H. Stiele, D. Hatzidimitriou, A. K. H. Kong, B. F. Williams, W. Pietsch, U. C. Kolb, F. Haberl, G. Sala

Abstract: We present detailed analysis of three globular cluster X-ray sources in the XMM-Newton extended survey of M31. The X-ray counterpart to the M31 globular cluster Bo 45 (XBo 45) was observed with XMM-Newton on 2006 December 26. Its combined pn+MOS 0.3--10 keV lightcurve exhibited a r.m.s variability of ~10%, and its 0.3--7.0 keV emission spectrum was well described by an absorbed power law with ph… ▽ More We present detailed analysis of three globular cluster X-ray sources in the XMM-Newton extended survey of M31. The X-ray counterpart to the M31 globular cluster Bo 45 (XBo 45) was observed with XMM-Newton on 2006 December 26. Its combined pn+MOS 0.3--10 keV lightcurve exhibited a r.m.s variability of ~10%, and its 0.3--7.0 keV emission spectrum was well described by an absorbed power law with photon index 1.44$\pm$0.12. Its variability and emission is characteristic of low mass X-ray binaries (LMXBs) in the low-hard state, whether the accretor is a neutron star or black hole. Such behaviour is typically observed at luminosities $\la$10% Eddington. However, XBo 45 exhibited this behaviour at an unabsorbed, 0.3--10 keV luminosity of 2.5$\pm0.2\times 10^{38}$ erg s$^{-1}$, or{~140%} Eddington for a 1.4 $M_{\odot}$ neutron star accreting hydrogen. Hence, we identify XBo 45 as a new candidate black hole LMXB. XBo 45 appears to have been consistently bright for ~30 years, consistent with theoretical prediction for a globular cluster black hole binary formed via tidal capture. Bo 375 was observed in the 2007, January 2 XMM-Newton observation, and has a two-component spectrum that is typical for a bright neutron star LMXB. Bo 135 was observed in the same field as Bo 45, and could contain either a black hole or neutron star. △ Less

Submitted 21 August, 2008; v1 submitted 1 July, 2008; originally announced July 2008.

Comments: Accepted by ApJ, 16 pages, 5 figures. This version includes the final changes made at the request of the referee

An X-ray spectral survey of the disc of M31 with XMM-Newton

Authors: L. Shaw Greening, R. Barnard, U. Kolb, C. Tonkin, J. P. Osborne

Abstract: We present the results of a complete spectral survey of the X-ray point sources detected in five XMM-Newton observations along the major axis of M31 but avoiding the central bulge, aimed at establishing the population characteristics of X-ray sources in this galaxy. We obtained background subtracted spectra and lightcurves for each of the 335 X-ray point sources detected across the five observat… ▽ More We present the results of a complete spectral survey of the X-ray point sources detected in five XMM-Newton observations along the major axis of M31 but avoiding the central bulge, aimed at establishing the population characteristics of X-ray sources in this galaxy. We obtained background subtracted spectra and lightcurves for each of the 335 X-ray point sources detected across the five observations from 2002. We also correlate our source list with those of earlier X-ray surveys and radio, optical and infra-red catalogues. Sources with more than 50 source counts are individually spectrally fit in order to create the most accurate luminosity functions of M31 to date. Based on the spectral fitting of these sources with a power law model, we observe a broad range of best fit photon index. From this distribution of best fit index, we identify 16 strong high mass X-ray binary system candidates in M31. We show the first cumulative luminosity functions created using the best fit spectral model to each source with more than 50 source counts in the disc of M31. The cumulative luminosity functions show a prominent, statistically significant flattening in the X-ray luminosity LX interval 37.0 \lesssim log LX erg s-1 \lesssim 37.5. Such a feature may also be present in the X-ray populations of several other galaxies, but at a much lower statistical significance. We investigate the number of AGN present in our source list and find that above LX ~1.4x1036 erg s-1 the observed population is statistically dominated by the point source population of M31. △ Less

Submitted 13 August, 2008; v1 submitted 31 May, 2008; originally announced June 2008.

Comments: accepted by A&A. 24 pages, 6 figures, 7 tables

A multi-coloured survey of NGC 253 with XMM-Newton: testing the methods used for creating luminosity functions from low-count data

Authors: R. Barnard, L. Shaw Greening, U. Kolb

Abstract: NGC 253 is a local, star-bursting spiral galaxy with strong X-ray emission from hot gas, as well as many point sources. We have conducted a spectral survey of the X-ray population of NGC 253 using a deep XMM-Newton observation.NGC 253 only accounts for ~20% of the XMM-Newton EPIC field of view, allowing us to identify ~100 X-ray sources that are unlikely to be associated with NGC\thinspace 253.… ▽ More NGC 253 is a local, star-bursting spiral galaxy with strong X-ray emission from hot gas, as well as many point sources. We have conducted a spectral survey of the X-ray population of NGC 253 using a deep XMM-Newton observation.NGC 253 only accounts for ~20% of the XMM-Newton EPIC field of view, allowing us to identify ~100 X-ray sources that are unlikely to be associated with NGC\thinspace 253. Hence we were able to make a direct estimate of contamination from e.g. foreground stars and background galaxies. X-ray luminosity functions (XLFs) of galaxy populations are often used to characterise their properties. There are several methods for estimating the luminosities of X-ray sources with few photons. We have obtained spectral fits for the brightest 140 sources in the 2003 XMM-Newton observation of NGC 253, and compare the best fit luminosities of those 69 non-nuclear sources associated with NGC 253 with luminosities derived using other methods. We find the luminosities obtained from these various methods to vary systematically by a factor of up to three for the same data; this is largely due to differences in absorption. We therefore conclude that assuming Galactic absorption is probably unwise; rather, one should measure the absorption for the population. A remarkable correlation has been reported between the XLFs of galaxies and their star formation rates. However, the XLFs used in that study were obtained using several different methods. If the sample galaxies were revisited and a single method were applied, then this correlation may become stronger still. △ Less

Submitted 9 May, 2008; originally announced May 2008.

Comments: Accepted for publication in the Monthly Notices of the Royal Astronomical Society (MNRAS). 17 pages, 7 figures

GRB 070201: A possible Soft Gamma Ray Repeater in M31

Authors: E. O. Ofek, M. Muno, R. Quimby, S. R. Kulkarni, H. Stiele, W. Pietsch, E. Nakar, A. Gal-Yam, A. Rau, P. B. Cameron, S. B. Cenko, M. M. Kasliwal, D. B. Fox, P. Chandra, A. K. H. Kong, R. Barnard

Abstract: The gamma-ray burst (GRB) 070201 was a bright short-duration hard-spectrum GRB detected by the Inter-Planetary Network (IPN). Its error quadrilateral, which has an area of 0.124 sq. deg, intersects some prominent spiral arms of the nearby M31 (Andromeda) galaxy. Given the properties of this GRB, along with the fact that LIGO data argues against a compact binary merger origin in M31, this GRB is… ▽ More The gamma-ray burst (GRB) 070201 was a bright short-duration hard-spectrum GRB detected by the Inter-Planetary Network (IPN). Its error quadrilateral, which has an area of 0.124 sq. deg, intersects some prominent spiral arms of the nearby M31 (Andromeda) galaxy. Given the properties of this GRB, along with the fact that LIGO data argues against a compact binary merger origin in M31, this GRB is an excellent candidate for an extragalactic Soft Gamma-ray Repeater (SGR) giant flare, with energy of 1.4x10^45 erg. Analysis of ROTSE-IIIb visible light observations of M31, taken 10.6 hours after the burst and covering 42% of the GRB error region, did not reveal any optical transient down to a limiting magnitude of 17.1. We inspected archival and proprietary XMM-Newton X-ray observations of the intersection of the GRB error quadrilateral and M31, obtained about four weeks prior to the outburst, in order to look for periodic variable X-ray sources. No SGR or Anomalous X-ray Pulsar (AXP) candidates (periods in range 1 to 20 s) were detected. We discuss the possibility of detecting extragalactic SGRs/AXPs by identifying their periodic X-ray light curves. Our simulations suggest that the probability of detecting the periodic X-ray signal of one of the known Galactic SGRs/AXPs, if placed in M31, is about 10% (50%), using 50 ks (2 Ms) XMM-Newton exposures. △ Less

Submitted 20 December, 2007; originally announced December 2007.

Comments: 7 pages, submitted to ApJ (Fig. 2 resolution reduced)

Journal ref: Astrophys.J.681:1464-1469,2008

A Multi-coloured survey of NGC 253 with XMM-Newton

Authors: R. Barnard, L. Shaw Greening, U. Kolb

Abstract: There is a large body of work that has used the excellent Chandra observations of nearby galaxies with neglible low mass X-ray binary (LMXB) populations. This has culminated in a ``Universal'' X-ray luminosity function (XLF) for high mass X-ray binaries (HMXBs). However, a number of methods have been used to convert from source intensities to luminosities when creating these XLFs. We have taken… ▽ More There is a large body of work that has used the excellent Chandra observations of nearby galaxies with neglible low mass X-ray binary (LMXB) populations. This has culminated in a ``Universal'' X-ray luminosity function (XLF) for high mass X-ray binaries (HMXBs). However, a number of methods have been used to convert from source intensities to luminosities when creating these XLFs. We have taken advantage of the XMM-Newton observations of the nearby starbursting spiral galaxy NGC 253 to test some of these methods. We find the luminosities derived from these various methods to vary by a factor of $\sim$3. We also find the most influential factor in the conversion from intensity to luminosity to be the absorption. We therefore conclude that a more consistent approach is required for determining the true Universal XLF for HMXBs. Ideally, this would involve individual spectral fitting of each X-ray source. Certainly, the line-of-sight absorption should be determined from the observations rather than assuming Galactic absorption. We find the best approach for obtaining an XLF from low-count data to be the splitting of the X-ray sources into two or more intensity intervals, and obtaining a conversion from intensity to flux for each group from spectral modelling of the summed spectrum of that group. △ Less

Submitted 25 October, 2007; originally announced October 2007.

Comments: To appear in the proceedings of "X-rays from Nearby Galaxies", 4-7 Septembeer 2007, 4 pages

On the artificial nature of aperiodic variability in XMM-Newton observations of M31 X-ray sources and the ultraluminous X-ray source NGC 4559 ULX-7

Authors: R. Barnard, S. Trudolyubov, U. C. Kolb, C. A. Haswell, J. P. Osborne, W. C. Priedhorsky

Abstract: Context: Power density spectra (PDS) that are characteristic of low mass X-ray binaries (LMXBs) have been previously reported for M31 X-ray sources, observed by XMM-Newton. However, we have recently discovered that these PDS result from the improper addition/subtraction of non-simultaneous lightcurves. Aims: To understand the properties and origins of the artefact. Methods: We re-analysed our XM… ▽ More Context: Power density spectra (PDS) that are characteristic of low mass X-ray binaries (LMXBs) have been previously reported for M31 X-ray sources, observed by XMM-Newton. However, we have recently discovered that these PDS result from the improper addition/subtraction of non-simultaneous lightcurves. Aims: To understand the properties and origins of the artefact. Methods: We re-analysed our XMM-Newton observations of M31 with non-simultaneous and simultaneous lightcurves, then combined simulated lightcurves at various intensities with various offsets and found that the artefact is more dependent on the offset than the intensity. Results: The lightcurves produced by the XMM-Newton Science Analysis Software (SAS) are non-synchronised by default. This affects not only the combination of lightcurves from the three EPIC detectors (MOS1, MOS2 and pn), but also background subtraction in the same CCD. It is therefore imperative that all SAS-generated lightcurves are synchronised by time filtering, even if the whole observation is to be used. We also find that the reported timing behaviour for NGC 4559 ULX-7 was also contaminated by the artefact; there is no significant variability in the correctly-combined lightcurves of NGC 4559 ULX-7. Hence, the classification of this source as an intermediate-mass black hole is no longer justified. Conclusions: While previous timing results from M31 have been proven wrong, and also the broken power law PDS in NGC 4559 ULX-7, XMM-Newton was able to detect aperiodic variability in just 3 ks of observations of NGC 5408 ULX1. Hence XMM-Newton remains a viable tool for analysing variability in extra-galactic X-ray sources. △ Less

Submitted 6 March, 2007; originally announced March 2007.

Comments: Accepted for publication in A&A; 6 pages, 4 figures

Artificial variability in XMM-Newton observations of X-ray sources: M31 as a case study

Authors: R. Barnard, S. Trudolyubov, C. A. Haswell, U. C. Kolb, J. P. Osborne, W. H. Priedhorsky

Abstract: Power density spectra (PDS) that are characteristic of low mass X-ray binaries (LMXBs) have been previously reported for M31 X-ray sources observed by XMM-Newton. However, we have recently discovered that these PDS are false positives resulting from the improper manipulation of non-simultaneous lightcurves. The lightcurves produced by the XMM-Newton Science Analysis Software (SAS) are non-synchr… ▽ More Power density spectra (PDS) that are characteristic of low mass X-ray binaries (LMXBs) have been previously reported for M31 X-ray sources observed by XMM-Newton. However, we have recently discovered that these PDS are false positives resulting from the improper manipulation of non-simultaneous lightcurves. The lightcurves produced by the XMM-Newton Science Analysis Software (SAS) are non-synchronised by default. This affects not only the combination of lightcurves from the three EPIC detectors (MOS1, MOS2 and pn), but also background subtraction in the same CCD. It is therefore imperative that all SAS-generated lightcurves are synchronised by time filtering, even if the whole observation is to be used. We combined simulated lightcurves at various intensities with various offsets and found that the artefact is more dependent on the offset than the intensity. While previous timing results from M31 have been proven wrong, and also the broken power law PDS in NGC 4559 ULX-7, XMM-Newton was able to detect aperiodic variability in just 3 ks of observations of NGC 5408 ULX1. Hence XMM-Newton remains a viable tool for analysing variability in extra-galactic X-ray sources. △ Less

Submitted 2 October, 2006; originally announced October 2006.

Comments: 6 pages, 5 figures; to appear in proceedings for "The multicoloured landscape of compact objects and their explosive progenitors", Cefalu', Sicily, 2006 June 11-24, AIP, in press

Discovery and modelling of disc precession in the M31 X-ray binary Bo 158?

Authors: R. Barnard, S. B. Foulkes, C. A. Haswell, U. Kolb, J. P. Osborne, J. R. Murray

Abstract: The low mass X-ray binary (LMXB) associated with the M31 globular cluster Bo 158 is known to exhibit intensity dips on a ~2.78 hr period. This is due to obscuration of the X-ray source on the orbital period by material on the outer edge of the accretion disc. However, the depth of dipping varied from <10% to \~83% in three archival XMM-Newton observations of Bo 158. Previous work suggested that… ▽ More The low mass X-ray binary (LMXB) associated with the M31 globular cluster Bo 158 is known to exhibit intensity dips on a ~2.78 hr period. This is due to obscuration of the X-ray source on the orbital period by material on the outer edge of the accretion disc. However, the depth of dipping varied from <10% to \~83% in three archival XMM-Newton observations of Bo 158. Previous work suggested that the dip depth was anticorrelated with the X-ray luminosity. However, we present results from three new XMM-Newton observations that suggest that the evolution of dipping is instead due to precession of the accretion disc. Such precession is expected in neutron star LMXBs with mass ratios <0.3 (i.e. with orbital periods <4 hr), such as the Galactic dipping LMXB 4U 1916-053. We simulated the accretion disc of Bo 158 using cutting-edge 3D smoothed particle hydrodynamics (SPH), and using the observed parameters. Our results show disc variability on two time-scales. The disc precesses in a prograde direction on a period of 81+/-3 hr. Also, a radiatively-driven disc warp is present in the inner disc, which undergoes retrograde precesson on a \~31 hr period. From the system geometry, we conclude that the dipping evolution is driven by the disc precession. Hence we predict that the dipping behaviour repeats on a ~81 hr cycle. △ Less

Submitted 2 November, 2005; originally announced November 2005.

Comments: 6 pages, 4 figures, to appear in proceedings "The X-ray Universe 2005", San Lorenzo de El Escoriale (Madrid, Spain), 26-30 September 2005

Cumulative luminosity functions of the X-ray point source population in M31

Authors: L. Shaw Greening, C. Tonkin, R. Barnard, U. Kolb, J. P. Osborne

Abstract: We present preliminary results from a detailed analysis of the X-ray point sources in the XMM-Newton survey of M31. These sources are expected to be mostly X-ray binaries. We have so far studied 225 of the 535 sources found by automated source detection. Only sources which were present in all three EPIC images were considered. X-ray binaries are identified by their energy spectrum and power dens… ▽ More We present preliminary results from a detailed analysis of the X-ray point sources in the XMM-Newton survey of M31. These sources are expected to be mostly X-ray binaries. We have so far studied 225 of the 535 sources found by automated source detection. Only sources which were present in all three EPIC images were considered. X-ray binaries are identified by their energy spectrum and power density spectrum. Unlike in other surveys we have obtained source luminosities from freely fit emission models. We present uncorrected luminosity functions of the sources analysed so far. △ Less

Submitted 20 September, 2005; originally announced September 2005.

Comments: 2 pages, 1 figure. To appear in proceedings of IAUS230

XMM-Newton reveals ~100 new LMXBs in M31 from variability studies

Authors: R. Barnard, L. Shaw Greening, C. Tonkin, U. Kolb, J. P. Osborne

Abstract: We have conducted a survey of X-ray sources in XMM-Newton observations of M31, examining their power density spectra (PDS) and spectral energy distributions (SEDs). Our automated source detection yielded 535 good X-ray sources; to date, we have studied 225 of them. In particular, we examined the PDS because low mass X-ray binaries (LMXBs) exhibit two distinctive types of PDS. At low accretion ra… ▽ More We have conducted a survey of X-ray sources in XMM-Newton observations of M31, examining their power density spectra (PDS) and spectral energy distributions (SEDs). Our automated source detection yielded 535 good X-ray sources; to date, we have studied 225 of them. In particular, we examined the PDS because low mass X-ray binaries (LMXBs) exhibit two distinctive types of PDS. At low accretion rates, the PDS is characterised by a broken power law, with the spectral index changing from ~0 to ~1 at some frequency in the range \~0.01--1 Hz; we refer to such PDS as Type A. At higher accretion rates, the PDS is described by a simple power law; we call these PDS Type B. Of the 225 sources studied to date, 75 exhibit Type A variability, and are almost certainly LMXBs, while 6 show Type B but not Type A, and are likely LMXBs. Of these 81 candidate LMXBs, 71 are newly identified in this survey; furthermore, they are mostly found near the centre of M31. Furthermore, most of the X-ray population in the disc are associated with the spiral arms, making them likely high mass X-ray binaries (HMXBs). In general these HMXBs do not exhibit Type A variability, while many central X-ray sources (LMXBs) in the same luminosity range do. Hence the PDS may distinguish between LMXBs and HMXBs in this luminosity range. △ Less

Submitted 19 September, 2005; originally announced September 2005.

Comments: 4 pages, 2 figures. To appear in proceedings of IAUS230: "Populations of High Energy Sources in Galaxies", 14-19 August 2005, Dublin, Eds E.J.A. Meurs and G. Fabbiano

The X-ray Properties of M101 ULX-1 = CXOKM101 J140332.74+542102

Authors: K. Mukai, M. Still, R. H. D. Corbet, K. D. Kuntz, R. Barnard

Abstract: We report our analysis of X-ray data on M101 ULX-1, concentrating on high state Chandra and XMM-Newton observations. We find that the high state of M101 ULX-1 may have a preferred recurrence timescale. If so, the underlying clock may have periods around 160 or 190 days, or possibly around 45 days. Its short-term variations resemble those of X-ray binaries at high accretion rate. If this analogy… ▽ More We report our analysis of X-ray data on M101 ULX-1, concentrating on high state Chandra and XMM-Newton observations. We find that the high state of M101 ULX-1 may have a preferred recurrence timescale. If so, the underlying clock may have periods around 160 or 190 days, or possibly around 45 days. Its short-term variations resemble those of X-ray binaries at high accretion rate. If this analogy is correct, we infer that the accretor is a 20-40 Msun object. This is consistent with our spectral analysis of the high state spectra of M101 ULX-1, from which we find no evidence for an extreme (> 10^40 ergs/s) luminosity. We present our interpretation in the framework of a high mass X-ray binary system consisting of a B supergiant mass donor and a large stellar-mass black hole. △ Less

Submitted 12 August, 2005; originally announced August 2005.

Comments: 23 pages, 7 figures, accepted for publication in the Astrophysical Journal

Journal ref: Astrophys.J.634:1085-1092,2005

Timing the bright X-ray population of the core of M31 with XMM-Newton

Authors: R. Barnard, U. Kolb, J. P. Osborne

Abstract: All 63 discrete X-ray sources brighter than L_{\rm X} = 5x10^{36} erg/s in any of four XMM-Newton observations of the core of M31 were surveyed for time variability over time-scales of seconds to thousands of seconds, and for intensity variations between observations. This population is likely to be dominated by low mass X-ray binaries (LMXBs). Analysis of the shapes of power density spectra (… ▽ More All 63 discrete X-ray sources brighter than L_{\rm X} = 5x10^{36} erg/s in any of four XMM-Newton observations of the core of M31 were surveyed for time variability over time-scales of seconds to thousands of seconds, and for intensity variations between observations. This population is likely to be dominated by low mass X-ray binaries (LMXBs). Analysis of the shapes of power density spectra (PDS) of these sources allows us to determine whether they are accreting at a high or low rate; in the observed frequency range a broken power law (Type A) PDS indicates a low accretion rate while a simple power law (Type B) indicates a high accretion rate. We obtained the 0.3--10 keV luminosities of the sources by modelling their spectral energy distributions. The luminosity function for Type A PDS appears to consist of two populations, which we tentatively classify as neutron star and black hole LMXBs. We find that 46 sources are likely X-ray binaries, 13 with possible black hole primaries. Remarkably, 5 out of the 13 black hole candidates appear persistently bright; there are no persistently bright black hole LMXBs known in our Galaxy. The lightcurves of 55 X-ray sources had a probability of variability >99% in at least one observation. Also, 57 of the sources show a luminosity variation > ~5$σ$ between observations; 7 of these are classed as transients, since they are absent in at least one observation, and vary in luminosity by a factor of > ~10. Scaling this sample with the known Galactic LMXBs, we find that the M31 sample has ~50% of the dippers, Z-sources and transients. △ Less

Submitted 17 February, 2006; v1 submitted 12 August, 2005; originally announced August 2005.

Comments: 28 pages, 9 figures, accepted for publication in A&A. This version contains major revisions over the previous version, with significantly improved interpretation and analysis