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.

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.

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.

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$\sigma$ improvement in $\chi^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$\sigma$ 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.

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.

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.

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.

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.

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).

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.

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\sigma 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.

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.

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.

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.

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

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.

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.

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.

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.

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.

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.

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.

[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.

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.

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.

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.

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.

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.

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$\sigma$ 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.

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

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.

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.

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.

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.

The well-studied X-ray/Optical/Radio supernova remnant DDB 1-15 (CXOM31 J004327.8+411829; r3-63) in M31 has been investigated with archival XMM-Newton and Chandra observations. The timing data from XMM-Newton reveals a power density spectrum (PDS) characteristic of accreting compact objects in X-ray binaries (XRBs). The PDS shows features typical of Roche lobe overflow accretion, hinting that the XRB is low-mass. The Chandra observations resolve the SNR into a shell and show a variable count rate at the 94% confidence level in the northwest quadrant. Together, these XMM-Newton and Chandra data suggest that there is an XRB in the SNR r3-63 and that the XRB is located in the northwestern portion of the SNR. The currently-available X-ray and optical data show no evidence that the XRB is high-mass. If the XRB is low-mass, r3-63 would be the first SNR found to contain a low-mass X-ray binary.

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

Stochastic variability in two out of four XMM-Newton observations of XMMU J004303+4115 along with its power spectra and X-ray luminosities suggest a low-mass X-ray binary (LMXB) with a black hole primary. However, Chandra observations resolve the object into two point sources. We use data from 35 Chandra observations to analyse the contributions of each source, and attribute the variability to CXOM31 J004303.2+411528 (known as r2-3), which varies in intensity by a factor of ~100 between observations. We assume that the power density spectra of LMXBs are governed by the luminosity, and that the transition between types of power density spectra occurs at some critical luminosity in Eddington units, l_c, that applies to all LMXBs. We use results from these XMM-Newton observations and past results from the available literature to estimate this transition luminosity, and find that all results are consistent with l_c ~0.1 in the 0.3--10 keV band. CXOM31 J004303.2+411528 exhibits a low accretion rate power density spectrum at a 0.3--10 keV luminosity of (5.3+/-0.6)x10^{37} erg/s. Known stellar mass black holes have masses of 4--15 M_{\odot}; hence our observations of CXOM31 J004303.2+411528 are consistent with l_c ~0.1 if it has a black hole primary.

Four XMM-Newton observations of the core of M31, spaced at 6 month intervals, show that the brightest point X-ray source, RX J0042.6+4115, has a 0.4-10 keV luminosity of ~5 10^38 erg/s, and exhibits significant variability in intensity and X-ray spectrum over a time scale of ~100 s including hard flares; such behaviour is only observed in Z-sources and transient blackhole binaries in our Galaxy. The lightcurves, X-ray spectra and hardness-intensity data from the four XMM-Newton observations all strongly suggest that it is a Z-source, bringing the total number of known Z-sources to nine.

We present results of a detailed study of X-ray flaring in the Z-track source Sco X-1 in a highly super-Eddington state made using high quality Rossi-XTE data from the PCA and HEXTE instruments. The emission model successfully used to explain the dipping LMXB, and other classes of LMXB in recent years, was applied to study the physical evolution along the Z-track which remains a major problem. This model consists of blackbody emission from the neutron star plus Comptonized emission from an extended accretion disk corona. As found in earlier work, major changes take place in the neutron star blackbody emission with kT increasing in flaring, and the blackbody radius R_BB increasing substantially to a maximum value of 9.4 +/- 0.6 km, consistent with the radius of the neutron star, after which R_BB decreases. Thus this result is a measurement of neutron star radius. The behaviour of Sco X-1 in flaring is compared with our previous results for the strong flaring that takes place in the bright dipping, flaring LMXB X 1624-490. Remarkably, during movement along the Normal Branch towards the apex with the Flaring Branch, the luminosities of both spectral components decrease, suggesting the possibility that Mdot may decrease on the Normal Branch, contrary to the widely-held view that Mdot increases monotonically along the Z-track. During flaring, we detect for the first time an increase of the Comptonization cut-off energy which may suggest heating of the ADC plasma by the neutron star flare. The energy of a broad Gaussian line at 6.4 keV does not change, but the intensity of the line increases in flaring suggesting either an increase in ADC size in flaring or the effects of irradiation by the neutron star.

Four XMM-Newton X-ray observations of the central region of the Andromeda Galaxy (M31) have revealed an X-ray source that varies in luminosity over \~1--3, 10^{38} erg s^{-1} between observations and also displays significant variability over time-scales of a few hundred seconds. The power density spectra of lightcurves obtained in the 0.3--10 keV energy band from the three EPIC instruments on board XMM-Newton are typical of disc-accreting X-ray binaries at low accretion rates, observed in neutron star binaries only at much lower luminosities (~10^{36} erg s^{-1}). However X-ray binaries with massive black hole primaries have exhibited such power spectra for luminosities >10^{38} erg s^{-1}. We discuss alternative possibilities where RX J0042.3+4115 may be a background AGN or foreground object in the field of view, but conclude that it is located within M31 and hence use the observed power spectra and X-ray luminosities to identify the primary as a black hole candidate.

The XMM-Newton satellite has made four observations of the core of M31 as part of a survey that will span the entire galaxy. The majority of the X-ray emission is contributed by point sources, chiefly X-ray Binaries (XB). Exciting early results from a survey of time variability in three XMM exposures of the M31 core reveal XB exhibiting Sco X-1 like flaring, dips and a burst. The flaring source in M31 also exhibits variability in intensity and spectral hardness that is characteristic of a Z-source.

We report results from analysis of the observations of the dipping low mass X-ray binary XB 1323-619 made with BeppoSAX} and Rossi-XTE. The dust-scattered halo contributes significantly in this source, and the observation made with BeppoSAX on 1997 August was used to provide MECS radial intensity profiles at several energies. From these, the halo fractions were obtained and thus an optical depth to dust scattering of 1.8 +/- 0.4 derived. In the Rossi-XTE observation of April 25-28, 1997, seven X-rays dips were observed together with 7 bursts repeating approximately periodically. Non-dip and dip PCA spectra can be well-described by assuming the emission consists of point-like blackbody emission identified with the neutron star, plus Comptonized emission from an extended ADC. The blackbody temperature is 1.79 +/- 0.21 keV and the cut-off power law photon index 1.61 +/- 0.04. Spectral evolution in dipping is well described by progressive covering of the extended Comptonizing region by absorber plus more rapid removal of the point-like blackbody. The effects of dust scattering and of the X-ray pulsar 1SAX J1324.4-6200 also in the field of view are included in the fitting. We detect an iron line at ~6.4 keV and its probable origin in the ADC is discussed.

We present results of an investigation of the physical changes taking place in the emission regions of the LMXB X1624-490 during strong flaring in RXTE observations. Based on the detailed light curve, we propose that the flaring consists of a superposition of X-ray bursts. It is shown that major changes take place in the blackbody emission component, the temperature kT_BB increasing to ~2.2 keV in flaring. Remarkably, the blackbody area decreases by a factor of ~5 in flaring. During flare evolution, the blackbody luminosity remains approximately constant, constituting a previously unknown Eddington limiting effect which we propose is due to radiation pressure of the blackbody as kT_BB increases affecting the inner disk or accretion flow resulting in a decreased emitting area on the star. We argue that the large decrease in area cannot be explained in terms of modification of the blackbody spectrum by electron scattering in the atmosphere of the neutron star. The height of the emitting region on the non-flaring neutron star is shown to agree with the height of the inner radiatively-supported disk as found for sources in the ASCA survey of LMXB of Church & Balucinska-Church (2001). The decrease in height during flaring is discussed in terms of possible models, including radial accretion flow onto the stellar surface and the theory of accretion flow spreading on the neutron star surface of Inogamov & Sunyaev (1999). We demonstrate that the intensity of the broad iron line at 6.4 keV is strongly correlated with the luminosity of the blackbody emission from the neutron star, and discuss the probable origin of this line in the ADC. Finally, possible reasons for non-detection of a reflection component in this source, and LMXB in general, are discussed.

We present results of a comprehensive study of the distribution of absorption dips with orbital phase in Cygnus X-1. Firstly, the distribution was obtained using archival data from all major X-ray observatories and corrected for the selection effect that phase zero (superior conjunction of the black hole) has been preferentially observed. Dip occurrence was seen to vary strongly with orbital phase \phi, with a peak at \phi ~ 0.95, i.e. was not symmetric about phase zero. Secondly, the RXTE ASM has provided continuous coverage of the Low State of Cygnus X-1 since Sept. 1996, and we have selected dip data based on increases in hardness ratio. The distribution, with much increased numbers of dip events, confirms that the peak is at \phi ~ 0.95, and we report the discovery of a second peak at \phi ~ 0.6. We attribute this peak to absorption in an accretion stream from the companion star HDE 226868. We have estimated the ionization parameter at different positions showing that radiative acceleration of the wind is suppressed by photoionization in particular regions in the binary system. To obtain the variation of column density with phase, we make estimates of neutral wind density for the extreme cases that acceleration of the wind is totally suppressed, or not suppressed at all. An accurate description will lie between these extremes. In each case, a strong variation of column density with orbital phase resulted, similar to the variation of dip occurrence. This provides evidence that formation of the blobs in the wind which lead to absorption dips depends on the density of the neutral component in the wind, suggesting possible mechanisms for blob growth.