X-ray Binaries Research Articles

No X-Rays or Radio from the Nearest Black Holes and Implications for Future Searches

Astrometry from the Gaia mission was recently used to discover the two nearest known stellar-mass black holes (BHs), Gaia BH1 and Gaia BH2. These objects are among the first stellar-mass BHs not discovered via X-rays or gravitational waves. Both systems contain ∼1 M ⊙ stars in wide orbits (a ≈ 1.4 au, 4.96 au) around ∼9 M ⊙ BHs, with both stars (solar-type main sequence star, red giant) well within their Roche lobes in Gaia BH1 and BH2, respectively. However, the BHs are still expected to accrete stellar winds, leading to potentially detectable X-ray or radio emission. Here, we report observations of both systems with the Chandra X-ray Observatory, the Very Large Array (for Gaia BH1) and MeerKAT (for Gaia BH2). We did not detect either system, leading to X-ray upper limits of L X < 9.4 × 1028 and L X < 4.0 × 1029 erg s−1 and radio upper limits of L r < 1.6 × 1025 and L r < 1.0 × 1026 erg s−1 for Gaia BH1 and BH2, respectively. For Gaia BH2, the non-detection implies that the accretion rate near the horizon is much lower than the Bondi rate, consistent with recent models for hot accretion flows. We discuss implications of these non-detections for broader BH searches, concluding that it is unlikely that isolated BHs will be detected via interstellar medium accretion in the near future. We also calculate evolutionary models for the binaries’ future evolution using Modules for Experiments in Stellar Astrophysics, and find that Gaia BH1 will be visible as a symbiotic BH X-ray binary for 5–50 Myr. Since no symbiotic BH X-ray binaries are known, this implies either that fewer than ∼104 Gaia BH1-like binaries exist in the Milky Way, or that they are common but have evaded detection.

An Optical Search for New Outbursting Low Mass X-Ray Binaries

Transient low-mass X-ray binaries (LMXBs) are discovered largely by X-ray and gamma-ray all-sky monitors. The X-ray outburst is also accompanied by an optical brightening, which empirically can precede the detection of X-rays. Newly sensitive optical synoptic surveys may offer a complementary pathway for discovery and potential for insight into the initial onset and propagation of the thermal instability that leads to the ionization of the accretion disk. We use the Zwicky Transient Facility (ZTF) alert stream to perform a comprehensive search at optical wavelengths for previously undiscovered outbursting LMXBs. Our pipeline first crossmatches the positions of the alerts to cataloged X-ray sources, and then analyzes the 30 day lightcurve of matched alerts by thresholding on differences with an 8 day exponentially weighted moving average. In addition to a 19 month long live search, we ran our pipeline over 4 yr of ZTF archival data, recovering four known LMXBs. We also independently detected an outburst of MAXI J1957+032 in the live search and found the first outburst of Swift J1943.4+0228, an unclassified X-ray transient, in 10 yr. Using Monte Carlo simulations of the Galactic LMXB population, we estimate that 29% of outbursting LMXBs are detectable by ZTF and that 4.4% of LMXBs would be present in the crossmatched X-ray catalogs, giving an estimated Galactic population of 3390−1930+3980 . We estimate that our current pipeline can detect 1.3% of all outbursting LMXBs, including those previously unknown, but that Rubin Observatory’s Legacy Survey of Space and Time will be able to detect 43% of outbursting LMXBs.

Analysis of Bright Source Hardness Ratios in the 4 yr Insight-HXMT Galactic Plane Scanning Survey Catalog

We conduct a statistical analysis of the hardness ratio (HR) for bright sources in the 4 yr Galactic Plane Scanning Survey catalog of Insight-HXMT. Depending on the stable (variable) flux F s (F v) or spectrum S s(S v) of each source, the bright sources are classified into three groups: F v&S v, F v&S s, and F s&S s. Our study of the HR characteristics in different types of sources reveals that accretion-powered neutron star (NS) low-mass X-ray binaries (LMXBs) exhibit softer energy spectra than NS high-mass X-ray binaries (HMXBs), but harder energy spectra than black hole binaries in most cases. This difference is probably due to their different magnetic field strengths. Additionally, F v&S v LMXBs tend to be harder than F v&S s LMXBs below 7 keV, while the opposite is true for HMXBs. Our results suggest that LMXBs may dominate unclassified sources, and NS binaries are likely to be the primary type of X-ray binaries with ambiguous compact stars. By comparing the HR of transient sources in their outburst and low-flux states, it is found that the averaged HR of four sources in the two states are roughly comparable within uncertainties. We also investigate the spatial properties of the three groups and find that F v&S v sources are mainly located in the longitude of −20° < l < 9°, F v&S s sources cross the Galactic Plane, and F s&S s sources are predominantly concentrated in 19° < l < 42°. In addition, analyzing the HR spatial distributions shows the absorption of soft X-rays (primarily below 2 keV) in the Galactic Plane.

AstroSat and NICER timing view of the Z−type neutron star X-ray binary GX 340 + 0

ABSTRACT The timing properties of the Z-type low-mass X-ray binaries provide insights into the emission components involved in producing the unique Z-shaped track in the hardness–intensity diagrams of these sources. In this work, we investigate the AstroSat and NICER observations of the GX 340 + 0 covering the complete ‘Z’-track from the horizontal branch (HB) to the extended flaring branch (EFB). For the first time, we present the Z-track as seen in soft X-rays using the AstroSat/Soft X-ray Telescope and NICER (the soft colour is defined as a ratio of 3–6 to 0.5–3 keV). The shape of the track is distinctly different in soft X-rays, strongly suggesting the presence of additional components active in soft X-rays. The detailed timing analysis revealed significant quasi-periodic oscillation throughout the HB and the normal branch (NB) using large area X-ray proportional counter and the first NICER detection of 33.1 ± 1.1 Hz HB oscillation (HBO) in 3–6 keV. The oscillations at the HB/NB vertex are observed to have higher frequencies (41–52 Hz) than the HBOs (16–31 Hz) and NB oscillations (6.2–8 Hz) but significantly lower rms (∼1.6 per cent). The HBO is also limited to the energy range 3–20 keV, indicating an association of HBO origin with the non-thermal component. It is also supported by earlier studies that found the strongest X-ray polarization during HB.

The First Polarimetric View on Quasiperiodic Oscillations in a Black Hole X-Ray Binary

We present the first polarimetric analysis of quasiperiodic oscillations (QPOs) in a black hole binary utilizing IXPE data. Our study focuses on Swift J1727.8–1613, which experienced a massive outburst that was observed by various telescopes across different wavelengths. The IXPE observation we studied was conducted during the hard-intermediate state. The polarization degree (PD) and polarization angle (PA) were measured at 4.28% ± 0.20% and 1.°9 ± 1.°4, respectively. Remarkably, significant QPO signals were detected during this observation, with a QPO frequency of approximately 1.34 Hz and a fractional rms amplitude of about 12.3%. Furthermore, we conducted a phase-resolved analysis of the QPO using the Hilbert–Huang transform technique. The photon index showed a strong modulation with respect to the QPO phase. In contrast, the PD and PA exhibit no modulations in relation to the QPO phase, which is inconsistent with the expectation of the Lense–Thirring precession of the inner flow. Further theoretical studies are needed to conform with the observational results.

Insight-HXMT observations of thermonuclear X-ray bursts in 4U 1636−53

ABSTRACT We conducted an analysis of 45 bursts observed from 4U 1636−53. To investigate the mechanism behind the light-curve profiles and the impact of thermonuclear X-ray bursts on the accretion environment in accreting neutron star low-mass X-ray binaries. This analysis employed both light-curve and time-resolved spectroscopy methodologies, with data collected by the Insight-Hard X-ray Modulation Telescope instrument. We found that 30 bursts exhibited similar light-curve profiles and were predominantly in the hard state, and two photospheric radius expansion (PRE) bursts were in the soft state. The light curves of most bursts did not follow a single exponential decay but displayed a dual-exponential behaviour. The initial exponent had a duration of approximately 6 s. We utilized both the standard method and the ‘fa’ method to fit the burst spectra. The majority of the ‘fa’ values exceeded 1, indicating an enhancement of the persistent emission during the burst. Under the two Comptonization components assumption, we suggest that the scattering of burst photons by the inner corona may mainly contribute to the persistent emission enhancement. We also observed an inverse correlation between the maximum fa and the persistent emission flux in the non-PRE burst. This anticorrelation suggests that when the accretion rate is lower, there is a greater enhancement of persistent emission during the burst peak. The prediction based on Poynting–Robertson drag (P–R drag) aligns with this observed anticorrelation.

The Common Envelope Evolution Outcome. II. Short-orbital-period Hot Subdwarf B Binaries Reveal a Clear Picture

Common envelope evolution (CEE) is vital for forming short-orbital-period compact binaries. It covers many objects, such as double compact merging binaries, Type Ia supernovae progenitors, binary pulsars, and X-ray binaries. Knowledge of the common envelope (CE) ejection efficiency still needs to be improved, though progress has been made recently. Short-orbital-period hot subdwarf B star (sdB) plus white dwarf (WD) binaries are the most straightforward samples with which to constrain CEE physics. We apply the known orbital period–WD mass relation to constrain the sdB progenitors of seven sdB+WD binaries with a known inclination angle. The average CE efficiency parameter is 0.32. This is consistent with previous studies. However, the CE efficiency need not be constant, but a function of the initial mass ratio, based on well-constrained sdB progenitor mass and evolutionary stage. Our results can be used as physical inputs for binary population synthesis simulations of related objects. A similar method can also be applied to study other short-orbital-period WD binaries.

Monitoring the X-Ray Variability of Bright X-Ray Sources in M33

We present a new five-epoch Chandra X-ray Observatory monitoring survey of the nearby spiral galaxy M33 which probes X-ray variability with time sampling between two weeks and four months. We characterize the X-ray variability of 55 bright point sources outside of the nucleus, many of which are expected to be high-mass X-ray binaries (HMXBs). We detect eight new candidate transients not detected in previous X-ray catalogs of M33 and discuss their possible nature. The final catalog includes 26 known HMXB candidates identified in the literature. We extend the baseline of the X-ray light curves up to 21 yr by including archival X-ray observations of these sources. We compare the detection and nondetection epochs of the sources to suites of simulated source duty cycles and infer that most of our detected sources have duty cycles >30%. We find only four sources whose detection patterns are consistent with having duty cycles below 30%. This large fraction of sources with high duty cycles is unexpected for a population of HMXBs; thus more frequent X-ray monitoring will likely reveal many more low duty cycle HMXBs in M33.

First Observational Evidence for an Interconnected Evolution between Time Lag and QPO Frequency among AGNs

Quasiperiodic oscillations (QPOs) have been widely observed in black hole X-ray binaries (BHBs), which often exhibit significant X-ray variations. Extensive research has explored the long-term evolution of the properties of QPOs in BHBs. In contrast, such evolution in active galactic nuclei (AGNs) has remained largely unexplored due to limited observational data. By using the 10 new XMM-Newton observations for the narrow-line Seyfert 1 galaxy RE J1034+396 from publicly available data, we analyze the characteristics of its X-ray QPOs and examine their long-term evolution. The hard-band (1–4 keV) QPOs are found in all 10 observations and the frequency of these QPOs evolves ranging at (2.47–2.83) × 10−4 Hz. Furthermore, QPO signals in the soft (0.3–1 keV) and hard bands exhibit strong coherence, although, at times, the variations in the soft band lead those in the hard band (the hard-lag mode), while at other times, it is the reverse (the soft-lag mode). The observations presented here serendipitously captured two ongoing lag reversals between these two modes within about two weeks, which are first seen in RE J1034+396 and also among all AGNs. A transition in QPO frequency also takes place within a two-week timeframe, two weeks prior to its corresponding lag reversal, indicating a possible coherence between the transitions of QPO frequency and lag mode with delay. The diagram of time lag versus QPO frequency clearly evidences this interconnected evolution with hysteresis, which is, for the first time, observed among AGNs.

A Fermi-LAT Study of Globular Cluster Dynamical Evolution in the Milky Way: Millisecond Pulsars as the Probe

Using archival Fermi-LAT data with a time span of ∼12 yr, we study the population of Millisecond Pulsars (MSPs) in Globular Clusters (GlCs) and investigate their dependence on cluster dynamical evolution in the Milky Way. We show that the γ-ray luminosity (L γ ) and emissivity (i.e., ϵ γ = L γ /M, with M the cluster mass) are good indicators of the population and abundance of MSPs in GlCs, and they are highly dependent on the dynamical evolution history of the host clusters. Specifically speaking, the dynamically older GlCs with more compact structures are more likely to have larger L γ and ϵ γ , and these trends can be summarized as strong correlations with cluster stellar encounter rate Γ and the specific encounter rate (Λ = Γ/M), with L γ ∝ Γ0.70±0.11 and ϵ γ ∝ Λ0.73±0.13 for dynamically normal GlCs. However, as GlCs evolve into deep core collapse, these trends are found to be reversed, implying that strong encounters may have lead to the disruption of Low-Mass X-ray Binaries and ejection of MSPs from core-collapsed systems. Besides, the GlCs are found to exhibit larger ϵ γ with increasing stellar mass function slope (ϵ γ ∝ 10(0.57±0.1)α ), decreasing tidal radius ( ϵγ∝Rt−1.0±0.22 ) and distances from the Galactic Center (GC, ϵγ∝Rgc−1.13±0.21 ). These correlations indicate that, as GlCs losing kinetic energy and spiral in toward the GC, tidal stripping and mass segregation have a preference in leading to the loss of normal stars from GlCs, while MSPs are more likely to concentrate to cluster center and be deposited into the GC. Moreover, we gauge ϵ γ of GlCs is ∼10–1000 times larger than the Galactic bulge, the latter is thought to reside thousands of unresolved MSPs and may be responsible for the GC γ-ray excess, which supports that GlCs are generous contributors to the population of MSPs in the GC.

A possible third body in the X-ray system GRS 1747−312 and models with higher order multiplicity

ABSTRACT GRS 1747−312 is a bright Low-Mass X-ray Binary in the globular cluster Terzan 6, located at a distance of 9.5 kpc from the Earth. It exhibits regular outbursts approximately every 4.5 months, during which periodic eclipses are known to occur. These eclipses have only been observed in the outburst phase, and are not clearly seen when the source is quiescent. Recent Chandra observations of the source were performed in 2019 June and 2021 April, June, and August. Two of these observations captured the source during its outburst, and showed clear flux decreases at the expected time of eclipse. The other two observations occurred when the source was quiescent. We present the discovery of a dip that occurred during the quiescent state. The dip is of longer duration and its time of occurrence does not fit the ephemeris of the shorter eclipses. We study the physical characteristics of the dip and determine that it has all the properties of an eclipse by an object with a well defined surface. We find that there are several possibilities for the nature of the object causing the 5.3 ks eclipse. First, GRS 1747−312 may be an X-ray triple, with an LMXB orbited by an outer third object, which could be an M-dwarf, brown dwarf, or planet. Secondly, there could be two LMXBs in close proximity to each other, potentially bound together. Whatever the true nature of the eclipser, its presence suggests that the GRS 1747−312 system is exotic.

Reflecting on naked singularities: iron line fitting as a probe of the cosmic censorship conjecture

ABSTRACT We demonstrate that the X-ray iron line fitting technique can be leveraged as a powerful probe of the cosmic censorship conjecture. We do this by extending existing emission line models to arbitrary spin parameters of the Kerr metric, no longer restricted to black hole metrics with |a•| &amp;lt; 1. We show that the emission lines from naked singularity metrics (|a•| &amp;gt; 1) show significant differences to their black hole counterparts, even for those metrics with identical locations of the innermost stable circular orbit, i.e. emission line fitting does not suffer from the degeneracy which affects continuum fitting approaches. These differences are entirely attributable to the disappearance of the event horizon for |a•| &amp;gt; 1. We highlight some novel emission line features of naked singularity metrics, such as ‘inverted’ emission lines (with sharp red wings and extended blue wings) and ‘triple lines’. The lack of detection of any of these novel features provides support of the cosmic censorship conjecture. We publicly release xspec packages skline and skconv which can now be used to probe the cosmic censorship conjecture in Galactic X-ray binaries and active galactic nuclei. The inclusion of super-extremal space–times can be alternatively posed as a way of stress testing conventional models of accretion.

The High Energy X-ray Probe (HEX-P): probing accretion onto stellar mass black holes

Accretion is a universal astrophysical process that plays a key role in cosmic history, from the epoch of reionization to galaxy and stellar formation and evolution. Accreting stellar-mass black holes in X-ray binaries are one of the best laboratories to study the accretion process and probe strong gravity—and most importantly, to measure the angular momentum, or spin, of black holes, and its role as a powering mechanism for relativistic astrophysical phenomena. Comprehensive characterization of the disk-corona system of accreting black holes, and their co-evolution, is fundamental to measurements of black hole spin. Here, we use simulated data to demonstrate how key unanswered questions in the study of accreting stellar-mass black holes will be addressed by the High Energy X-ray Probe (HEX-P). HEX-P is a probe-class mission concept that will combine high spatial resolution X-ray imaging and broad spectral coverage (0.2–80 keV) with a sensitivity superior to current facilities (including XMM-Newton and NuSTAR) to enable revolutionary new insights into a variety of important astrophysical problems. We illustrate the capability of HEX-P to: 1) measure the evolving structures of black hole binary accretion flows down to low (≲ 0.1%) Eddington-scaled luminosities via detailed X-ray reflection spectroscopy; 2) provide unprecedented spectral observations of the coronal plasma, probing its elusive geometry and energetics; 3) perform detailed broadband studies of stellar mass black holes in nearby galaxies, thus expanding the repertoire of sources we can use to study accretion physics and determine the fundamental nature of black holes; and 4) act as a complementary observatory to a range of future ground and space-based astronomical observatories, thus providing key spectral measurements of the multi-component emission from the inner accretion flows of black hole X-ray binaries.

Determining cosmological growth parameter for stellar-mass black holes

ABSTRACT It has recently been suggested that black holes (BHs) may grow with time, so that their mass is proportional to the cosmological scale factor to the power n, with suggested values n ≈ 3 for supermassive BHs in elliptical galaxies. Here, we test these predictions with stellar-mass BHs in X-ray binaries using their masses and ages. We perform two sets of tests to assess the compatible values of n. First, we assume that no compact object grows over the Tolman–Oppenheimer–Volkoff limit which marks the borderline between neutron stars and BHs. We show that half of BHs would be born with a mass below this limit if n = 3 applies. The possibility that all BHs were born above the limit is rejected at $4\, \sigma$ if n = 3 applies. In the second test, we assume that masses of BHs at their formation stay the same over cosmic history. We compare the mass distribution of the youngest BHs, which could have not grown yet, to their older counterparts. Distributions are compatible for $n = -0.9^{+1.3}_{-4.6}$, with n = 3 excluded formally with 87 per cent confidence. This result may be biased, because massive BHs tend to have a massive companion. Correcting for this bias yields n ≈ 0. We can therefore conclude that while our results are not a clear rejection of BH scaling with n = 3, we show that n = 0 is much more consistent with the data.

The high energy X-ray probe (HEX-P): probing the physics of the X-ray corona in active galactic nuclei

The hard X-ray emission in active galactic nuclei (AGN) and black hole X-ray binaries is thought to be produced by a hot cloud of electrons referred to as the corona. This emission, commonly described by a power law with a high-energy cutoff, is suggestive of Comptonization by thermal electrons. While several hypotheses have been proposed to explain the origin, geometry, and composition of the corona, we still lack a clear understanding of this fundamental component. NuSTAR has been playing a key role improving our knowledge of X-ray coronæ thanks to its unprecedented sensitivity above 10 keV. However, these constraints are limited to bright, nearby sources. The High Energy X-ray Probe (HEX-P) is a probe-class mission concept combining high spatial resolution X-ray imaging and broad spectral coverage (0.2–80 keV) with a sensitivity superior to current facilities. In this paper, we highlight the major role that HEX-P will play in further advancing our insights of X-ray coronæ notably in AGN. We demonstrate how HEX-P will measure key properties and track the temporal evolution of coronæ in unobscured AGN. This will allow us to determine their electron distribution and test the dominant emission mechanisms. Furthermore, we show how HEX-P will accurately estimate the coronal properties of obscured AGN in the local Universe, helping address fundamental questions about AGN unification. In addition, HEX-P will characterize coronæ in a large sample of luminous quasars at cosmological redshifts for the first time and track the evolution of coronæ in transient systems in real time. We also demonstrate how HEX-P will enable estimating the coronal geometry using spectral-timing techniques. HEX-P will thus be essential to understand the evolution and growth of black holes over a broad range of mass, distance, and luminosity, and will help uncover the black holes’ role in shaping the Universe.

Insights on the Optical and Infrared Nature of MAXI J0709-159: Implications for High-Mass X-ray Binaries

In our previous study (Bhattacharyya et al., 2022), HD 54786, the optical counterpart of the MAXI J0709-159 system, was identified to be an evolved star, departing from the main sequence, based on comparisons with non-X-ray binary systems. In this paper, using color-magnitude diagram (CMD) analysis for High-Mass X-ray Binaries (HMXBs) and statistical t-tests, we found evidence supporting HD 54786’s potential membership in both Be/X-ray binaries (BeXRBs) and supergiant X-ray binaries (SgXBs) populations of HMXBs. Hence, our study points towards dual optical characteristics of HD 54786, as an X-ray binary star and also belonging to a distinct evolutionary phase from BeXRB towards SgXB. Our further analysis suggests that MAXI J0709-159, associated with HD 54786, exhibits low-level activity during the current epoch and possesses a limited amount of circumstellar material. Although similarities with the previously studied BeXRB system LSI +61° 235 (Coe et al., 1994) are noted, continued monitoring and data collection are essential to fully comprehend the complexities of MAXI J0709-159 and its evolutionary trajectory within the realm of HMXBs.

From Feast to Famine: A Systematic Study of Accretion onto Oblique Pulsars with 3D GRMHD Simulations

Disk-fed accretion onto neutron stars can power a wide range of astrophysical sources ranging from X-ray binaries, to accretion-powered millisecond pulsars, ultraluminous X-ray sources, and gamma-ray bursts. A crucial parameter controlling the gas–magnetosphere interaction is the strength of the stellar dipole. In addition, coherent X-ray pulsations in many neutron star systems indicate that the star's dipole moment is oblique relative to its rotation axis. Therefore, it is critical to systematically explore the 2D parameter space of the star's magnetic field strength and obliquity, which is what this work does, for the first time, in the framework of 3D general-relativistic magnetohydrodynamics. If the accretion disk carries its own vertical magnetic field, this introduces an additional factor: the relative polarity of the disk and stellar magnetic fields. We find that depending on the strength of the stellar dipole and the star–disk relative polarity, the neutron star's jet power can either increase or decrease with increasing obliquity. For weak dipole strength (equivalently, high accretion rate), the parallel polarity results in a positive correlation between jet power and obliquity, whereas the antiparallel orientation displays the opposite trend. For stronger dipoles, the relative-polarity effect disappears, and jet power always decreases with increasing obliquity. The influence of the relative polarity gradually disappears as obliquity increases. Highly oblique pulsars tend to have an increased magnetospheric radius, a lower mass accretion rate, and enter the propeller regime at lower magnetic moments than aligned stars.

Statistical study of a large and cleaned sample of ultraluminous and hyperluminous X-ray sources

Context. Ultraluminous and hyperluminous X-ray (ULX and HLX) sources could constitute interesting laboratories to further improve our understanding of the supermassive black hole growth through super-Eddington accretion episodes and successive mergers of lighter holes. ULXs are thought to be powered by super-Eddington accretion onto stellar-mass compact objects, while HLXs are of an unknown nature, but they could be good candidates for accreting intermediate mass black holes (IMBHs). However, a significant portion of the sample of ULX and HLX candidates derived from catalogue searches are in fact background active galactic nuclei (AGN). Aims. Here we build samples of ULXs and HLXs from the three largest X-ray catalogues available, compiled from XMM-Newton, Swift-XRT, and Chandra detections, and the GLADE catalogue containing 1.7 million galaxies at D &lt; 1000 Mpc. We aim to characterise the frequency, environment, hardness, and variability of ULXs and HLXs to better assess their differences and understand their populations. Methods. After a thorough classification of these X-ray sources, we were able to remove ∼42% of sources with a signal-to-noise ratio (S/N) &gt; 3 which were shown to be contaminants, to obtain the cleanest sample of ULXs and HLXs to date. From a sample of 1342 ULXs and 191 HLXs detected with a S/N &gt; 3σ, we study the occupation fraction, hardness, variability, radial distribution, and preferred environment of the sources. We built their Malmquist-corrected X-ray luminosity functions (XLFs) and compared them with previous studies. Thanks to the unprecedented size of the sample, we were able to statistically compare ULXs and HLXs and assess the differences in their nature. The interpretation of HLXs as IMBHs is investigated. Results. A significant break is seen in the XLF at ∼1040 erg s−1. With our ULX sample, containing ≲2% of contaminants, we are able to confirm that ULXs are located preferentially in spiral galaxies and galaxies with higher star formation rates. While X-ray binaries (XRBs), ULXs, and most HLXs share common hardness and variability distributions, a fraction of HLXs appear significantly softer. Unlike ULXs, HLXs seem to reside equally in spiral as well as lenticular and elliptical galaxies. We note that 35% of the HLX candidates have an optical counterpart, and we estimate the mass of 120 of them to be in the range of 2 × 103 − 105 M⊙. Most of the HLX population is found to be consistent with an accreting massive black hole in a dwarf galaxy satellite of the primary host. This diverse nature needs to be confirmed with deeper optical and infrared observations, as well as upcoming X-ray facilities.

Three-dimensional GRMHD Simulations of Neutron Star Jets

Neutron stars and black holes in X-ray binaries are observed to host strong collimated jets in the hard spectral state. Numerical simulations can act as a valuable tool in understanding the mechanisms behind jet formation and its properties. Although there have been significant efforts in understanding black hole jets from general relativistic magnetohydrodynamic (GRMHD) simulations in recent years, neutron star jets still remain poorly explored. We present the results from three-dimensional GRMHD simulations of accreting neutron stars with oblique magnetospheres for the very first time. The jets in our simulations are produced due to the anchored magnetic field of the rotating star in analogy with the Blandford–Znajek process. We find that for accreting stars, the star–disk magnetic field interaction plays a significant role, and as a result, the jet power becomes directly proportional to , where Φjet is the open flux in the jet. The jet power decreases with increasing stellar magnetic inclination, and finally, for an orthogonal magnetosphere, it reduces by a factor of ≃2.95 compared to the aligned case. We also find that in the strong propeller regime, with a highly oblique magnetosphere, the disk-induced collimation of the open stellar flux preserves parts of the striped wind, resulting in a striped jet.

LONG TERM OPTICAL VARIABILITY OF TWO HMXB

We discuss the optical light curves of two Be X-ray Binaries, IGR J06074+2205 and SAX J2103.5+4545, recovered from the ATLAS, ZTF and ASAS-SN databases. Both sources show long term optical variability of 620 and 420 days, respectively, with color redder when brighter. We suggest that this is due to the precession of the circumstellar disk. Another possibility is the propagation of a density wave in the disk. Only these two sources show a large amplitude, periodic optical variability from a sample of 16 well-studied high-mass X-ray Binaries (HMXB) with a Main Sequence primary star.