Black Hole Low-mass X-ray Binaries Research Articles

An Alternative Channel to Black Hole Low-mass X-Ray Binaries: Dynamical Friction of Dark Matter?

Both the anomalous magnetic braking of Ap/Bp stars and the surrounding circumbinary disk models can account for the formation of black hole (BH) low-mass X-ray binaries (LMXBs), while the simulated effective temperatures of the donor stars are significantly higher than the observed values. Therefore, the formation of BH LMXBs is still not completely understood. In this work, we diagnose whether the dynamical friction between dark matter and the companion stars can drive BH binaries to evolve toward the observed BH LMXBs and alleviate the effective temperature problem. Assuming that there exists a density spike of dark matter around BH, the dynamical friction can produce an efficient angular momentum loss, driving BH binaries with an intermediate-mass companion star to evolve into BH LMXBs for a spike index higher than γ = 1.58. Our detailed stellar evolution models show that the calculated effective temperatures can match the observed value of most BH LMXBs for a spike index range of γ = 1.7–2.1. However, the simulated mass-transfer rates when γ = 2.0 and 2.1 are too high to be consistent with the observed properties showing that BH LMXBs appear as soft X-ray transients. Therefore, the dynamical friction of dark matter can only alleviate the effective temperature problem of those BH LMXBs with a relatively short orbital period.

On the Magnetic Braking Law in Black Hole Low-mass X-Ray Binaries

Magnetic braking (MB) plays an important role in the evolution of close low-mass X-ray binaries (LMXBs). It is also essential to the formation of ultracompact X-ray binaries (UCXBs). There have been lively investigations on the MB mechanism(s) in both single stars and close binaries including cataclysmic variables and neutron star (NS) LMXBs but with diverse conclusions. In this paper we explore the effect of MB on the black hole (BH) LMXB evolution. We combine binary population synthesis with detailed binary evolution to obtain the expected properties of the Galactic BH LMXB population. The simulated results are compared with the observational data including the BH mass, companion mass, companion temperature, orbital period, and mean accretion rate. Our results reveal that the MB laws with relatively low efficiency (i.e., RM12 and RVJ83) exhibit better agreement with observations, contrary to what was found for NS LMXBs. This raises the interesting question about whether MB really follows the same unified law in different types of binaries. We also predict that only a very small fraction (≲2.5%) of BH LMXBs can evolve to be UCXBs. This explains why no BH UCXB has been discovered so far.

A dependable distance estimator to black hole low-mass X-ray binaries

ABSTRACT Black hole low-mass X-ray binaries (BH-LMXBs) are excellent observational laboratories for studying many open questions in accretion physics. However, determining the physical properties of BH-LMXBs necessitates knowing their distances. With the increased discovery rate of BH-LMXBs, many canonical methods cannot produce accurate distance estimates at the desired pace. In this study, we develop a versatile statistical framework to obtain robust distance estimates soon after discovery. Our framework builds on previous methods where the soft spectral state and the soft-to-hard spectral state transitions, typically present in an outbursting BH-LMXB, are used to place constraints on mass and distance. We further develop the traditional framework by incorporating general relativistic corrections, accounting for spectral/physical parameter uncertainties, and employing assumptions grounded in current theoretical and observational knowledge. We tested our framework by analysing a sample of 50 BH-LMXB sources using X-ray spectral data from the Swift/XRT, MAXI/GSC, and RXTE/PCA missions. By modelling their spectra, we applied our framework to 26 sources from the 50. Comparison of our estimated distances to previous distance estimates indicates that our findings are dependable and in agreement with the accurate estimates obtained through parallax and H i absorption methods. Investigating the accuracy of our constraints, we have found that estimates obtained using both the soft and transition spectral information have a median uncertainty (1σ) of 20 per cent, while estimates obtained using only the soft spectral state spectrum have a median uncertainty (1σ) of around 50 per cent. Furthermore, we have found no instrument-specific biases.

WIMP constraints from black hole low-mass X-ray binaries

The abnormally fast orbital decay observed in the black hole (BH) Low-Mass X-ray binaries (BH-LMXB) A0620-00 and XTE J1118+480 can be explained by the dynamical friction between Dark Matter (DM) and the companion star orbiting around the low-mass BH (≃ a few M ⊙) of the system. In this case the value of the index γ sp of the DM spike surrounding the BH can be pinned down with an accuracy of ≃ a few percent, way better than that for much bigger systems such as the super massive BHs (SMBHs) in the Galactic Center or in M87. We have used data from XTE J1118+480 to put bounds on the WIMP annihilation cross section times velocity ⟨σv⟩, assuming that DM annihilation is driven by the χχ → bb̅ annihilation channel and that it proceeds in s-wave. The bounds are driven by the radio synchrotron signal produced by e ± final states propagating in the magnetic field in the vicinity of the BH. We find that for DM masses m χ up to the TeV scale XTE J1118+480 allows to constrain ⟨σv⟩ well below the standard value ⟨σv⟩thermal, corresponding to the observed DM relic density in the Universe for a thermal WIMP. On the other hand, for m χ ≳ 15 GeV, the bounds from the SMBHs in the GC or in M87 do not reach ⟨σv⟩thermal when the very large uncertainties on the corresponding spike indices are taken into account, in spite of potentially producing much larger DM densities compared to XTE J1118+480. Our bounds for XTE J1118+480 have a mild sensitivity on the effect of spatial diffusion (which implies at most a weakening of the bounds of a factor ≲ 6 at large m χ). However, diffusion is instrumental in enhancing the sensitivity of the results upon the intensity of the magnetic field. In particular, our bounds rest on the assumption that the magnetic field B reaches the equipartition value B eq. We find that a reduction factor of the magnetic field B eq/B larger than about 14 at low m χ, becoming progressively smaller at higher WIMP masses, would be sufficient to relax the XTE J1118+480 bound to the level of other existing bounds. Recent estimates, albeit not conclusive, may suggest values of B eq/B in BH-LMXB systems as large as 20. This implies that the intensity of the magnetic field in BH-LMXB systems represents the major uncertainty in using them as an alternative to heavier BHs to search for WIMPs.

The current state of disk wind observations in BHLMXBs through X-ray absorption lines in the iron band

The presence of blueshifted absorption lines in the X-ray spectra of black hole low-mass X-ray binaries is the telltale mark of massive outflows called winds. These signatures are found almost exclusively in soft states of high-inclined systems, hinting at equatorial ejections originating from the accretion disk and deeply intertwined with the evolution of the outburst patterns displayed by these systems. In the wake of the launch of the new generation of X-ray spectrometers, studies of wind signatures remain mostly restricted to single sources and outbursts, with some of the recent detections departing from the commonly expected behaviors. We thus give an update to the current state of iron band absorption line detections through the analysis of all publicly available XMM-Newton-pn and Chandra-HETG exposures of known black hole low-mass X-ray binary candidates. Our results agree with previous studies, as our wind detections are exclusively found in dipping, high-inclined sources and almost exclusively in bright (LX > 0.01LEdd) soft (HR < 0.8) states with blueshift values generally restricted to a few 100 km s−1. The line parameters indicate similar properties between objects and outbursts of single sources, and despite more than 20 yr of data, very few sources have the HID sampling necessary to properly study the evolution of the wind during a single outburst. We provide an online tool with details of the wind signatures and outburst evolution data for all sources in our sample.

Unveiling hidden variability components in accreting X-ray binaries using both the Fourier power and cross-spectra

ABSTRACT We present a novel method for measuring the lags of (weak) variability components in neutron-star and black-hole low-mass X-ray binaries (LMXBs). For this we assume that the power and cross-spectra of these sources consists of a number of components that are coherent in different energy bands, but are incoherent with one another. The technique is based on fitting simultaneously the power spectrum (PS) and the Real and Imaginary parts of the cross-spectrum (CS) with a combination of Lorentzian functions. We show that, because the PS of LMXBs is insensitive to signals with a large Imaginary part and a small Real part in the CS, this approach allows us to uncover new variability components that are only detected in the CS. We also demonstrate that, contrary to earlier claims, the frequency of the type-C quasi-periodic oscillation (QPO) in the black-hole binary GRS 1915+105 does not depend on energy. Rather, the apparent energy dependence of the QPO frequency can be explained by the presence of a separate QPO component with a slightly higher frequency than that of the QPO, whose rms amplitude increases faster with energy than the rms amplitude of the QPO. From all the above we conclude that, as in the case of the PS, the CS of black-hole and neutron-star binaries can be fitted by a combination of Lorentzian components. Our findings provide evidence that the frequency-dependent part of the transfer function of these systems can be described by a combination of responses, each of them acting over relatively well-defined time-scales. This conclusion challenges models that assume that the main contribution to the lags comes from a global, broadband, transfer function of the accreting system.

Simultaneous NICER and NuSTAR Observations of the Ultracompact X-Ray Binary 4U 0614+091

We present the first joint NuSTAR and NICER observations of the ultracompact X-ray binary 4U 0614+091. This source shows quasiperiodic flux variations on the timescale of ∼days. We use reflection modeling techniques to study various components of the accretion system as the flux varies. We find that the flux of the reflected emission and the thermal components representing the disk and the compact object trend closely with the overall flux. However, the flux of the power-law component representing the illuminating X-ray corona scales in the opposite direction, increasing as the total flux decreases. During the lowest flux observation, we see evidence of accretion disk truncation from roughly 6 gravitational radii to 11.5 gravitational radii. This is potentially analogous to the truncation seen in black hole low-mass X-ray binaries, which tends to occur during the low/hard state at sufficiently low Eddington ratios.

An empirical connection between line-emitting regions and X-rays heating the accretion disc in BH-LMXB MAXI J1820+070

ABSTRACT The recurring transient outbursts in low-mass X-ray binaries (LMXBs) provide ideal laboratories to study the accretion process. Unlike their supermassive relatives, LMXBs are far too small and distant to be imaged directly. Fortunately, phase-resolved spectroscopy can provide an alternative diagnostic to study their highly complex, time-dependent accretion discs. The primary spectral signature of LMXBs are strong, disc-formed emission lines detected at optical wavelengths. The shape, profile, and appearance/disappearance of these lines change throughout a binary orbit, and thus, can be used to trace how matter in these discs behaves and evolves over time. By combining a Swift multiwavelength monitoring campaign, phase-resolved spectroscopy from the Gran Telescopio Canarias (GTC) and Liverpool Telescope, and modern astrotomography techniques, we find a clear empirical connection between the line emitting regions and physical properties of the X-rays heating the disc in the black hole LMXB MAXI J1820+070 during its 2018 outburst. In this paper, we show how these empirical correlations can be used as an effective observational tool for understanding the geometry and structure of a LMXB accretion disc and present further evidence for an irradiation-driven warped accretion disc present in this system.

Chance Coincidences between Black Hole Low-mass X-Ray Binaries and Supernova Remnants

I argue that black hole low-mass X-ray binaries (BH LMXBs) are very unlikely to be physically associated with supernova remnants (SNRs). The timescales of BH LMXBs are so much longer than those of SNRs, that there is only a 0.2% chance of any BH LMXB being identified within its natal SNR. However, the probability of a BH LMXB being projected within an SNR is significant; I estimate that two BH LMXBs should be projected within SNRs from our perspective. I look more closely at the suggestion by Balakrishnan and collaborators of an association between the BH X-ray binary Swift J1728.9-3613 and the SNR G351.9-0.9, and show that this is most likely a chance coincidence.

NuSTAR View of the R − Γ Correlation in the Hard State of Black Hole Low-mass X-Ray Binaries

The power law and reflection emission have been observed in the X-ray spectra of both black hole X-ray binaries (BHXRBs) and active galactic nuclei (AGNs), indicating a common physical origin of the X-ray emission from these two types of sources. The relevant parameters describing the shape of both components and the potential correlation between these parameters can provide important clues on the geometric and physical properties of the disk and the corona in these sources. In this work, we present a positive correlation between the photon index Γ and the reflection strength R for the low-mass BHXRBs in the hard state by modeling NuSTAR data, which is qualitatively consistent with the previous studies. We compare our results with the predictions from different theoretical disk-corona models. We show that the R − Γ correlation found in this work seems to favor the moving corona model proposed by Beloborodov. Our results indicate that the coronal geometry varies significantly among BHXRBs. We further compare our results with that of AGNs. We find that the reflection strength R is smaller than unity in the hard state of BHXRBs, while it can be as large as ∼5 in AGNs, which implies that the variations of the disk-coronal geometry of AGNs are more vigorous than that of the BHXRBs in the hard state.

Close encounters of tight binary stars with stellar-mass black holes

ABSTRACT Strong dynamical interactions among stars and compact objects are expected in a variety of astrophysical settings, such as star clusters and the disks of active galactic nuclei. Via a suite of three-dimensional hydrodynamics simulations using the moving-mesh code arepo, we investigate the formation of transient phenomena and their properties in close encounters between an $2\, {\rm M}_{\odot }$ or $20\, {\rm M}_{\odot }$ equal-mass circular binary star and single $20\, {\rm M}_{\odot }$ black hole (BH). Stars can be disrupted by the BH during dynamical interactions, naturally producing electromagnetic transient phenomena. Encounters with impact parameters smaller than the semimajor axis of the initial binary frequently lead to a variety of transients whose electromagnetic signatures are qualitatively different from those of ordinary disruption events involving just two bodies. These include the simultaneous or successive disruptions of both stars and one full disruption of one star accompanied by successive partial disruptions of the other star. On the contrary, when the impact parameter is larger than the semimajor axis of the initial binary, the binary is either simply tidally perturbed or dissociated into bound and unbound single stars (‘micro-Hills’ mechanism). The dissociation of $20\, {\rm M}_{\odot }$ binaries can produce a runaway star and an active BH moving away from one another. Also, the binary dissociation can either produce an interacting binary with the BH, or a non-interacting, hard binary; both could be candidates of BH high- and low-mass X-ray binaries. Hence, our simulations especially confirm that strong encounters can lead to the formation of the (generally difficult to form) BH low-mass X-ray binaries.

Dual-corona Comptonization model for the type-b quasi-periodic oscillations in GX 339-4

ABSTRACT Characterizing the fast variability in black hole low-mass X-ray binaries (BHXBs) can help us to understand the geometrical and physical nature of the innermost regions of these sources. Particularly, type-B quasi-periodic oscillations (QPOs), observed in BHXBs during the soft-intermediate state (SIMS) of an outburst, are believed to be connected to the ejection of a relativistic jet. The X-ray spectrum of a source in the SIMS is characterized by a dominant soft blackbody-like component – associated with theVaccretion disc – and a hard component – associated with a Comptonizing region or corona. Strong type-B QPOs were observed by NICER and AstroSat in GX 339-4 during its 2021 outburst. We find that the fractional rms spectrum of the QPO remains constant at ∼1 per cent for energies below ∼1.8 keV, and then increases with increasing energy up to ∼17 per cent at 20–30 keV. We also find that the lag spectrum is ‘U-shaped’, decreasing from ∼1.2 rad at 0.7 keV to 0 rad at ∼3.5 keV, and increasing again at higher energies up to ∼0.6 rad at 20–30 keV. Using a recently developed time-dependent Comptonization model, we fit simultaneously the fractional rms and lag spectra of the QPO, and the time-averaged energy spectrum of GX 339-4 to constrain the physical parameters of the region responsible for the variability we observe. We suggest that the radiative properties of the type-B QPOs observed in GX 339-4 can be explained by two physically-connected comptonizing regions that interact with the accretion disc via a feedback loop of X-ray photons.

Low mass X-ray binary simulation data release

We present DABS (Database of Accreting Binary Simulations), an open-access database of modelled Low Mass X-ray Binaries (LMXBs). DABS has been created using evolutionary tracks of neutron star and black hole LMXBs, spanning a large set of initial conditions for the accretor mass, donor mass, and orbital period. The LMXBs are evolved with the Convection and Rotation Boosted Magnetic Braking prescription. The most important asset of this online database is the tool PEAS (Progenitor Extractor for Accreting Systems). This tool can be used to predict the progenitors of any user-entered LMXB system and view their properties before the start of mass transfer. This prediction can facilitate preliminary searches for the progenitors of observed LMXBs, which can help in streamlining further detailed analyses. The PEAS tool can also be used to constrain population synthesis techniques that specialize in supernova kicks in binaries and common envelope outcomes.

Optical and X-ray follow-up to a globular cluster ultraluminous X-ray source in NGC 4472

ABSTRACT NGC 4472 is home to five ultraluminous X-ray sources hosted by globular clusters. These sources have been suggested as good black hole candidates in extragalactic globular clusters—a highly sought after population that may provide observational information regarding the progenitors of merging black hole binaries. In this body of work, we present X-ray and optical follow-up to one of these sources, CXOUJ1229410+075744 (GCU1). We find no evidence of [OIII] optical emission in GCU1 which indicates a lack of significant evidence for super-Eddington outflows, unlike what is seen in a handful of ULXs in extragalactic GCs. X-ray monitoring from 2019 to 2021 shows no detected X-ray emission above a few × 1038 erg/s. Comparisons of the multiwavelength properties to disc-dominated, near Eddington Galactic black hole low-mass X-ray binaries (GRS 1915+105 and XTEJ1817-330) suggests that GCU1 may show similar behaviour to GRS 1915+105 in terms of X-ray variability and similar relationships between LX and kT, with GCU1 showing maximum X-ray luminosities on order of higher magnitude.

Radio monitoring of transient Be/X-ray binaries and the inflow–outflow coupling of strongly magnetized accreting neutron stars

ABSTRACT Strongly magnetized (B ≥ 1012 G) accreting neutron stars (NSs) are prime targets for studying the launching of jets by objects with a solid surface; while classical jet-launching models predict that such NSs cannot launch jets, recent observations and models argue otherwise. Transient Be/X-ray binaries (BeXRBs) are critical laboratories for probing this poorly explored parameter space for jet formation. Here, we present the coordinated monitoring campaigns of three BeXRBs across four outbursts: giant outbursts of SAX 2103.5+4545, 1A 0535+262, and GRO J1008–57, as well as a Type-I outburst of the latter. We obtain radio detections of 1A 0535+262 during ten out of twenty observations, while the other targets remained undetected at typical limits of 20–50 $\mu$Jy. The radio luminosity of 1A 0535+262 positively correlates with its evolving X-ray luminosity, and inhabits a region of the LX–LR plane continuing the correlation observed previously for the BeXRB Swift J0243.6+6124. We measure a BeXRB LX–LR coupling index of β = 0.86 ± 0.06 ($L_R \propto L_X^\beta$), similar to the indices measured in NS and black hole low-mass X-ray binaries. Strikingly, the coupling’s LR normalization is ∼275 and ∼6.2 × 103 times lower than in those two comparison samples, respectively. We conclude that jet emission likely dominates during the main peak of giant outbursts, but is only detectable for close-by or super-Eddington systems at current radio sensitivities. We discuss these results in the broader context of X-ray binary radio studies, concluding that our results suggest how supergiant X-ray binaries may host a currently unidentified additional radio emission mechanism.

Population synthesis of AX J1745.6−2901 X-ray nova-type binaries with rapidly decreasing orbital periods

ABSTRACT The neutron star low-mass X-ray binary (LMXB) AX J1745.6−2901 was detected with an anomalously fast decrease of its orbital period. The decreasing rate of the orbital period exceeds the contribution of all processes extracting angular momentum from the binary star in the standard model. Using the scenario machine code, we conducted a population synthesis study of X-ray novae with neutron stars to investigate a possible formation and evolution of such binaries. Such close LMXBs should experience a preceding common envelope stage, in which the magnetic fields of the low-mass main-sequence donor stars can be dramatically amplified. Our calculations show that the magnetic stellar wind of the optical companion can efficiently extract angular momentum from the binary systems, and produce the observed orbital-period derivatives of AX J1745.6−2901 and black hole LMXBs. The estimated values of the required magnetic field induction are the following: Bd ≈ 400 G (AX J1745.6−2901), Bd ≈ 1500 G (KV UMa), Bd ≈ 400 G (A0620−00) and Bd ≈ 1800 G (Nova Muscae). We successfully reproduced the current observational abundance of such anomalous neutron star X-ray novae, and computed the appropriate value of the parameter of magnetic braking λMSW (0.8−0.6 for Roche lobe filling stars and 0.4−0.15 for binaries with partial Roche lobe filling).

The evolving properties of the corona of GRS 1915+105: a spectral-timing perspective through variable-Comptonization modelling

ABSTRACT The inverse Compton process by which soft photons are up-scattered by hot electrons in a corona plays a fundamental role in shaping the X-ray spectra of black hole (BH) low-mass X-ray binaries (LMXBs), particularly in the hard and hard-intermediate states. In these states, the power-density spectra of these sources typically show Type-C low-frequency quasi-periodic oscillations (QPOs). Although several models have been proposed to explain the dynamical origin of their frequency, only a few of those models predict the spectral-timing radiative properties of the QPOs. Here, we study the physical and geometrical properties of the corona of the BH-LMXB GRS 1915+105 based on a large sample of observations available in the RXTE archive. We use a recently developed spectral-timing Comptonization model to fit simultaneously the energy-dependent fractional rms amplitude and phase-lag spectra of the Type-C QPO in 398 observations. For this, we include spectral information gathered from fitting a Comptonization model to the corresponding time-averaged spectra. We analyse the dependence of the physical and geometrical properties of the corona upon the QPO frequency and spectral state of the source, the latter characterized by the hardness ratio. We find consistent trends in the evolution of the corona size, temperature, and feedback (the fraction of the corona photons that impinge back on to the disc) that persist for roughly 15 yr. By correlating our observations with simultaneous radio-monitoring of the source at 15 GHz, we propose a scenario in which the disc–corona interactions connect with the launching mechanism of the radio jet in this source.

The NICER “Reverberation Machine”: A Systematic Study of Time Lags in Black Hole X-Ray Binaries

We perform the first systematic search of all NICER archival observations of black hole (and candidate) low-mass X-ray binaries for signatures of reverberation. Reverberation lags result from the light travel time difference between the direct coronal emission and the reflected disk component, and therefore their properties are a useful probe of the disk-corona geometry. We detect new signatures of reverberation lags in eight sources, increasing the total sample from three to 11, and study the evolution of reverberation lag properties as the sources evolve in outbursts. We find that in all of the nine sources with more than one reverberation lag detection, the reverberation lags become longer and dominate at lower Fourier frequencies during the hard-to-soft state transition. This result shows that the evolution in reverberation lags is a global property of the state transitions of black hole low-mass X-ray binaries, which is valuable in constraining models of such state transitions. The reverberation lag evolution suggests that the corona is the base of a jet that vertically expands and/or gets ejected during state transition. We also discover that in the hard state, the reverberation lags get shorter, just as the quasiperiodic oscillations (QPOs) move to higher frequencies, but then in the state transition, while the QPOs continue to higher frequencies, the lags get longer. We discuss the implications of the coronal geometry and physical models of QPOs in light of this new finding.

The Disk Veiling Effect of the Black Hole Low-mass X-Ray Binary A0620-00*

Abstract The optical light curves of quiescent black hole low-mass X-ray binaries often exhibit significant nonellipsoidal variabilities, showing the photospheric radiation of the companion star is veiled by other sources of optical emission. Assessing this “veiling” effect is critical to the black hole mass measurement. Here in this work, we carry out a strictly simultaneous spectroscopic and photometric campaign on the prototype of black hole low-mass X-ray binary A0620-00. We find that for each observation epoch, the extra optical flux beyond a pure ellipsoidal modulation is positively correlated with the fraction of veiling emission, indicating the accretion disk contributes most of the nonellipsoidal variations. Meanwhile, we also obtain a K2V spectral classification of the companion, as well as the measurements of the companion’s rotational velocity v sin i = 83.8 ± 1.9 km s−1 and the mass ratio between the companion and the black hole q = 0.063 ± 0.004.

The Observed Mass Distribution of Galactic Black Hole LMXBs Is Biased against Massive Black Holes

Abstract The discovery of gravitational-wave radiation from merging black holes (BHs) also uncovered BHs with masses in the range of ≈20–160 M ⊙. In contrast, the most massive Galactic stellar-mass BH currently known has a mass of ≈21 M ⊙. While low-mass X-ray binaries (LMXBs) will never independently evolve into a binary BH system, and binary evolution effects can play an important role in explaining the different BH masses found through studies of X-ray binaries and gravitational-wave events, (electromagnetic) selection effects may also play a role in this discrepancy. Assuming BH LMXBs originate in the Galactic plane, we show that the spatial distributions of the current samples of confirmed and candidate BH LMXBs are both biased to sources that lie at a large distance from the plane. Specifically, most of the confirmed and candidate BH LMXBs are found at a Galactic height larger than three times the scale height for massive star formation. In addition, the confirmed BH LMXBs are found at larger distances to the Galactic center than the candidate BH LMXBs. Interstellar absorption makes candidate BH LMXBs in the plane and bulge too faint for a dynamical mass measurement using current instrumentation. Given the observed and theoretical evidence for BH natal and/or Blaauw kicks, their relation with BH mass and binary orbital period, and the relation between outburst recurrence time and BH mass, the observational selection effects imply that the current sample of confirmed BH LMXBs is biased against the most massive BHs.