Compact Objects Research Articles

New charged anisotropic solution in f(Q)-gravity and effect of non-metricity and electric charge parameters on constraining maximum mass of self-gravitating objects

In the present article, A new class of singularity-free charged anisotropic stars is derived in f(Q)-gravity regime. To solve the field equations, we assume a particular form of anisotropy along with an electric field and obtain a new exact solution in f(Q)-gravity. The explicit mathematical expression for the model parameters is derived by the smooth joining of the obtained solutions with the exterior Reissner–Nordstrom de-Sitter solution across the bounding surface of a compact star along with the requirement that the radial pressure vanishes at the boundary. We have modeled four self-gravitating pulsar objects such as LMC X-4, PSR J1903+327, PSR J1614-2230, and GW190814 in our current study and predict the radii of these objects that fall between 8 and 10 km. Furthermore, the physical validity of the solution is performed for self-gravitating object PSR J1614-2230 with mass 1.97±0.04M⊙\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$1.97\\pm 0.04~M_{\\odot }$$\\end{document} with radius 10 km. The solution successfully fulfills all the physical requirements along with the stability and hydrostatic equilibrium conditions for a well-behaved model. The non-metricity f(Q)-parameter χ1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\chi _{_1}$$\\end{document} and electric charge parameter η\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\eta $$\\end{document} play an important role in the maximum mass of the objects. The maximum mass increases when χ1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\chi _{_1}$$\\end{document} and η\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\eta $$\\end{document} increase but a non-collapsing stable object can be obtained when χ1≤0.0205\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\chi _{_1}\\le 0.0205$$\\end{document} and η≤0.0006\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\eta \\le 0.0006$$\\end{document}.

From higher-spin gauge interactions to Compton amplitudes for root-Kerr

We develop massive higher-spin theory as a framework for describing dynamics of rotating compact objects, such as Kerr black holes. In this paper, we explore gauge interactions up to quartic order and corresponding Compton amplitudes of higher-spin massive objects coupled to electromagnetism and Yang-Mills theory. Their classical counterparts are known as root-Kerr gauge-theory solutions, whose amplitudes are closely related to those of Kerr black holes. We use three distinct approaches: (i) massive higher-spin gauge symmetry to introduce cubic interactions for all spins and the quartic interactions up to spin 3, which is implemented both off shell and via Ward identities; (ii) a chiral higher-spin approach to construct quartic Lagrangians with correct degrees of freedom to all spins; (iii) on-shell functional patterns before and after taking the classical limit to constrain the Compton amplitudes. As final results, we arrive at simple local formulae for the candidate root-Kerr Compton amplitudes both in the quantum regime and classical limit, to all orders in spin. This is a precursor to the gravitational Kerr case, which is presented in a follow-up paper.

Signatures of low-mass black hole–neutron star mergers

ABSTRACT The recent observation of the GW230529 event indicates that black hole–neutron star binaries can contain low-mass black holes. Since lower mass systems are more favourable for tidal disruption, such events are promising candidates for multimessenger observations. In this study, we employ five finite-temperature, composition-dependent matter equations of state and present results from ten 3D general relativistic hydrodynamic simulations for the mass ratios $q = 2.6$ and 5. Two of these simulations target the chirp mass and effective spin parameter of the GW230529 event, while the remaining eight contain slightly higher mass black holes, including both spinning ($a_{\mathrm{ BH}} = 0.7$) and non-spinning ($a_{\mathrm{ BH}} = 0$) models. We discuss the impact of the equation of state, spin, and mass ratio on black hole–neutron star mergers by examining both gravitational-wave and ejected matter properties. For the low-mass ratio model we do not see fast-moving ejecta for the softest equation of state model, but the stiffer model produces on the order of $10^{-6}\,\mathrm{ M}_\odot$ of fast-moving ejecta, expected to contribute to an electromagnetic counterpart. Notably, the high-mass ratio model produces nearly the same amount of total dynamical ejecta, but yields 52 times more fast-moving ejecta than the low-mass ratio system. In addition, we observe that the black-hole spin tends to decrease the amount of fast-moving ejecta while increasing significantly the total ejected mass. Finally, we note that the disc mass tends to increase as the neutron star compactness decreases.

Model of a charged compact star: A brief study in f(T) gravity

Using modified [Formula: see text] gravity and the diagonal tetrad, we propose a new kind of analytical solution to describe anisotropic charged stellar structures in this study. We obtain precise solutions by applying the linear form of [Formula: see text] as [Formula: see text] in the framework of Tolman–Kuchowicz physically sound metric potentials. The stellar structures of the compact star EXO 1785-248 are then investigated using the model by incorporating a linear equation of state relating the radial pressure [Formula: see text] and the matter density [Formula: see text]. We computed the closed-form expressions for the model parameters and illustrated their characteristics. The comprehensive graphical analysis demonstrates the scientific plausibility, causality, and dynamical stability of the model describing charged anisotropic stars. Finally, using the M-R curve, we estimate the numerical values of mass for different values of [Formula: see text]. The observational data of the neutron star candidates 4U 1820-30, PSR J1614-2230, and PSR J2215 + 5135 are covered by the maximum allowable masses for different values of [Formula: see text] additionally, when [Formula: see text], the charged compact star’s maximum permissible mass is calculated to be [Formula: see text], which could indicate the secondary component of the GW190814 event discovered by LIGO/VIRGO experiments.

Merger Precursor: Year-long Transients Preceding Mergers of Low-mass Stripped Stars with Compact Objects

Binary mass transfer can occur at high rates due to rapid expansion of the donor’s envelope. In the case where mass transfer is unstable, the binary can rapidly shrink its orbit and lead to a merger. In this work we consider the appearance of the system preceding merger, specifically for the case of a low-mass ( ≈2.5- 3M⊙) helium star with a neutron star (NS) companion. Modeling the mass transfer history as well as the wind launched by super-Eddington accretion onto the NS, we find that such systems can power slowly rising transients with timescales as long as years, and luminosities of ∼1040- 1041 erg s −1 from optical to UV. The final explosion following the merger (or core-collapse of the helium star in some cases) leads to an interaction-powered transient with properties resembling Type Ibn supernovae (SNe), possibly with a bright early peak powered by shock cooling emission for merger-powered explosions. We apply our model to the Type Ibn SN 2023fyq, that displayed a long-term precursor activity from years before the terminal explosion.

Bose–Einstein Condensation dark matter models generated by gravitational decoupling

A study on massive compact objects influenced by Bose–Einstein condensation (BEC) dark matter (DM) with a gravitational decoupling approach is presented within the context of f(Q) gravity, assuming a linear form of f(Q). In order to obtain physically valid solutions to the basic field equations of the decoupled systems, we take the Tolman IV potential ansatz for the seed system and the BEC-DM density profile for the new source into consideration. After that, we apply the boundary conditions on the decoupled system with Schwarzschild-de Sitter spacetime as the exterior metric. Eventually, we obtain a complete non-singular solution to the decoupled system, consisting of a strange stellar configuration with a BEC-DM density profile. All the physical properties, such as the effective energy density, pressure along the radial as well as the tangential direction, and anisotropy in the system, are rigorously investigated. Analyzing the stability conditions of the decoupled stellar configuration, we study mass–radius (M−R) relations with observational constraints related to the supermassive compact stars, such as PSR J0740+6620, PSR J2215+5135, GW190814, and GW200210, with observed masses larger than or equal to 2 M⊙. Notably, we examine the influence of DM on constraining the M−R measurements of the observed stars, which satisfy the MIT bag model equation of state as well as the condition of mimicking, i.e., ρθ=ρDM, in the background of f(Q) gravity. Interestingly, any possible conversion of massive compact stars to smaller black holes can be restrained, as the presence of a DM halo in the strange stellar system reduces the maximum mass effectively, as indicated by the M−R curves. This result regarding the maximum mass of the compact star can be associated with the astrophysical data concerning the mass gap of the events GW190814 and GW200210. Specifically, the present model predicts the radius for PSR J0740+6620 as 13.69−0.10+0.09 km, which nearly matches the NICER and XMM-Newton data.

Axial gravitational perturbations of slowly rotating compact objects in general relativity and beyond

Axial gravitational perturbations of slowly rotating compact objects in general relativity and beyond

Accreting Schwarzschild-like compact object: Plasma-photon interaction and stability

Accreting Schwarzschild-like compact object: Plasma-photon interaction and stability

Anisotropic Durgapal-Fuloria compact stars in f(R) gravity

This study presented a new exact solution for anisotropic compact stellar objects within the framework of f(R)=R+αR2 gravity. In this context, the Durgapal-Fuloria metric potential has been employed to solve the field equation derived for f(R) theory. Furthermore, we have derived the generalized Darmois-Israel junction condition necessary for seamlessly connecting the interior region to the Schwarzschild exterior metric across the boundary hypersurface of the star in the context of f(R) gravity, and the interior solution is matched with the Schwarzschild exterior metric over the bounding surface of a compact star. These junction conditions stipulate that the pressure must not be zero at the boundary and should be proportional to the non-linear terms of f(R) gravity, a crucial aspect often overlooked by many researchers when investigating compact stellar models. Additionally, we derived the values of these parameters by using observational data of various compact stars (CSs), namely Her X-1, SAX J1808.4-3658, SMC X-1, LMC X-4, Cen X-3, 4U 1820-30, PSR J1903+327, 4U 1608-52, Vela X-1, and PSR J1416-2230. This approach enables us to investigate the comprehensive analysis of solutions numerically and graphically. We conducted various physical tests, including gradient of energy density and pressures, anisotropy, stability, equilibrium conditions, energy-density constraints, mass function, compactness, redshift, and adiabatic index, to assess the feasibility of our models. Our findings demonstrate the consistent behavior of our models provides a satisfactory physical situation as far as the observational results are confirmed.

Fast inspirals and the treatment of orbital resonances

Abstract Extreme mass ratio inspirals (EMRIs), where a compact object orbits a massive black hole, are a key source of gravitational waves for the future Laser Interferometer Space Antenna (LISA). Due to their small mass ratio, ε ~ 10-4 -- 10-7 , the binary evolves slowly and EMRI signals will be in-band for years. Additionally, astrophysical EMRIs are expected to have complex dynamics featuring both spin-precession and eccentricity. A standard approach to modelling these inspirals is via the method of osculating geodesics (OG) which we employ along with a toy model for the gravitational self-force. Using this method requires resolving tens of thousands radial and polar orbital librations over the long duration of the signal which makes the inspiral trajectory expensive to compute.In this work we accelerate these calculations by employing Near-Identity (averaging) Transformations. 
However, this averaging technique breaks down at orbital resonances where the radial and polar frequencies are an integer ratio of each other. Thus, we switch to a partial averaging transformation in the vicinity of the resonance where the dynamics are characterised by the slow evolution of the so-called ``resonant phase". Additionally, we develop an optimal switching criterion to minimise the computation time while maximising accuracy. We find the error in the waveform phase is improved from O(ε-1/2) in the fully averaged scheme to O(ε4/7) in the switching scheme. At the same time, this scheme improves the scaling of the computation time from being inversely proportional to ε using OG, to a very weak scaling with ε. This results in a speed-up of at least two orders of magnitude for LISA EMRIs with room for further optimisation.

Physical characteristics of anisotropic solutions in f(Q,T) gravity under the vanishing complexity, embedding class one, conformally flat and conformally Killing conditions

The foundation of this paper is to present a comparative study of the properties of anisotropic solutions for compact stars in the framework of modified f(Q,T) gravity for the first time under the schemes of vanishing complexity, embedding class one, conformally flat and conformally Killing, where Q is the nonmetricity scalar and T is the trace of the energy-momentum tensor. The model f(Q,T)=ϕQ+αQn+1−γ(1−e−Q/γ)+βT is chosen in order to compare the results of the various forms f(Q,T) and f(Q) gravity theories. The f(Q,T) and f(Q) gravities in this model are then further simplified to depend on parameters ϕ,α,γ&β. Then, we discussed several approaches to determining the spherically symmetric spacetime components. Schwarzschild geometry represents the exterior spacetime, while a Tolman-Kuchowiz spacetime is used to examine the spherically symmetric spacetime inside the interior geometry. Many properties of compact stars are investigated, like energy density, energy conditions, pressure profiles, sound speeds, gradients, adiabatic index, TOV equation, mass function, compactness, and redshift function are all carefully examined in order to complete the analysis. We have compiled results for both stable and unstable scenarios, derived from gravity theories categorized into various cases, with solutions considered under different frameworks such as Tolman-Kuchowiz spacetime, embedding class 1 spacetime, the vanishing complexity condition, the conformally flat condition, and the conformal Killing vector condition.

Long-term evolutionary links between the isolated neutron star populations

ABSTRACT We have investigated the evolutionary connections of the isolated neutron star (NS) populations including radio pulsars (RPs), anomalous X-ray pulsars (AXPs), soft gamma repeaters (SGRs), dim isolated NSs (XDINs), ‘high-magnetic field’ RPs (‘HBRPs’), central compact objects (CCOs), rotating radio transients (RRATs), and long-period pulsars (LPPs) in the fallback disc model. The model can reproduce these NS families as a natural outcome of different initial conditions (initial period, disc mass, and dipole moment, μ) with a continuous μ distribution in the $\sim 10^{27} - 5 \times 10^{30}$ G cm$^3$ range. Results of our simulations can be summarized as follows: (1) A fraction of ‘HBRPs’ with relatively high μ evolve into the persistent AXP/SGR properties, and subsequently become LPPs. (2) Persistent AXP/SGRs do not have evolutionary links with CCOs, XDINs, and RRATs. (3) For a wide range of μ, most RRATs evolve passing through RP or ‘HBRP’ properties during their early evolutionary phases. (4) A fraction of RRATs which have the highest estimated birth rate seem to be the progenitors of XDINs. (5) LPPs, whose existence was predicted by the fallback disc model, are the sources evolving in the late stage of evolution before the discs become inactive. These results provide concrete support to the ideas proposing evolutionary connections between the NS families to account for the ‘birth rate problem’, the discrepancy between the cumulative birth rate estimated for these systems and the core-collapse supernova rate.

CXOU J005245.0-722844: Discovery of a be star / white dwarf binary system in the SMC via a very fast, super-eddington X-ray outburst event

Abstract CXOU J005245.0-722844 is an X-ray source in the Small Magellanic Cloud (SMC) that has long been known as a Be/X-ray binary (BeXRB) star, containing an OBe main sequence star and a compact object. In this paper, we report on a new very fast X-ray outburst from CXOU J005245.0-722844. X-ray observations taken by Swift constrain the duration of the outburst to less than 16 days and find that the source reached super-Eddington X-ray luminosities during the initial phases of the eruption. The XRT spectrum of CXOU J005245.0-722844 during this outburst reveals a super-soft X-ray source, best fit by an absorbed thermal blackbody model. Optical and Ultraviolet follow-up observations from the Optical Gravitational Lensing Experiment (OGLE), Asteroid Terrestrial-impact Last Alert System (ATLAS), and Swift identify a brief ∼0.5 magnitude optical burst coincident with the X-ray outburst that lasted for less than 7 days. Optical photometry additionally identifies the orbital period of the system to be 17.55 days and identifies a shortening of the period to 17.14 days in the years leading up to the outburst. Optical spectroscopy from the Southern African Large Telescope (SALT) confirms that the optical companion is an early-type OBe star. We conclude from our observations that the compact object in this system is a white dwarf (WD), making this the seventh candidate Be/WD X-ray binary. The X-ray outburst is found to be the result of a very-fast, ultra-luminous nova similar to the outburst of MAXI J0158-744.

Imaging fermionic dark matter cores at the centre of galaxies

ABSTRACT Current images of the supermassive black hole (SMBH) candidates at the centre of our Galaxy and M87 have opened an unprecedented era for studying strong gravity and the nature of relativistic sources. Very-long-baseline interferometry data show images consistent with a central SMBH within General Relativity (GR). However, it is essential to consider whether other well-motivated dark compact objects within GR could produce similar images. Recent studies have shown that dark matter (DM) haloes modelled as self-gravitating systems of neutral fermions can harbour very dense fermionic cores at their centres, which can mimic the space–time features of a black hole (BH). Such dense, horizonless DM cores can satisfy the observational constraints: they can be supermassive and compact and lack a hard surface. We investigate whether such cores can produce similar observational signatures to those of BHs when illuminated by an accretion disc. We compute images and spectra of the fermion cores with a general-relativistic ray tracing technique, assuming the radiation originates from standard $\alpha$ discs, which are self-consistently solved within the current DM framework. Our simulated images possess a central brightness depression surrounded by a ring-like feature, resembling what is expected in the BH scenario. For Milky Way-like haloes, the central brightness depressions have diameters down to ${\sim} 35\, \mu \text{as}$ as measured from a distance of approximately $8\,$ kpc. Finally, we show that the DM cores do not possess photon rings, a key difference from the BH paradigm, which could help discriminate between the models.

True unicorns and false positives: Simulated probabilities of dark massive companions to bright stars

Abstract Many compact objects (black holes and neutron stars) exist in binaries. These binaries are normally discovered through their interactions, either from accretion as an X-ray binary or collisions as a gravitational wave source. However, the majority of compact objects in binaries should be non-interacting. Recently proposed discoveries have used radial velocities of a bright star (main sequence or evolved) that are indicative of a massive but dark companion, which is inferred to be a compact object. Unfortunately, this burgeoning new field has been hindered by false positives, including the “Unicorn” (V723 Mon) which was initially believed to be a red giant/black hole binary before being refuted. In this work, we investigate the evolution of stellar binary populations over time, using the binary evolution code COSMIC to simulate binary populations and determine the probability of a candidate object being either a “true Unicorn” (actual compact objects in binaries) or a false positive. We find that main sequence stars have a higher true Unicorn probability than red giants or naked helium stars (an exposed core of an evolved star), particularly if the companion is more massive and is ≥3 times less luminous than the MS star. We also find that a top-heavy initial mass function raises the true Unicorn probability further, that super-solar metallicity reduces the probability, and that most true Unicorns are found at periods ≤100 days. Finally, we find that a significant fraction of true Unicorns do not evolve into x-ray binaries during the age of the universe.

The Evolution of Massive Binary Stars

Massive stars play a major role in the evolution of their host galaxies and serve as important probes of the distant Universe. It has been established that the majority of massive stars reside in close binaries and interact with their companion stars during their lifetimes. Such interactions drastically alter their life cycles and complicate our understanding of their evolution, but are also responsible for the production of interesting and exotic interaction products. ▪Extensive observation campaigns with well-understood detection sensitivities have enabled the conversion of observed properties into intrinsic characteristics, facilitating a direct comparison to theory.▪Studies of large samples of massive stars in our Galaxy and the Magellanic Clouds have unveiled new types of interaction products, providing critical constraints on the mass transfer phase and the formation of compact objects.▪The direct detection of gravitational waves has revolutionized the study of stellar mass compact objects, providing a new window to study massive star evolution. Their formation processes are, however, still unclear. The known sample of compact object mergers will increase by orders of magnitude in the coming decade, which is vastly outgrowing the number of stellar-mass compact objects detected through electromagnetic radiation.

Finch–Skea quintessence models in non-conservative theory of gravity

This study is dedicated to presenting a new solution of the field equations in the Rastall theory with a quintessence field defined by the parameter ωq as −1<ωq<−13 by considering the isotropic matter content inside the sphere. The Finch–Skea ansatz (FS) is used in a static and spherically symmetric geometry to obtain the feasible relativistic solution. The results obtained in the physical evaluation are analyzed analytically and graphically. In the appropriate limit of the Rastall coupling parameter, one can regain the original results in the General Relativity. This complete analysis considers five different compact stars: HerX−1 with mass 0.88M⊙ and radius 7.7 km, VelaX−12 with mass 1.77M⊙ and radius 9.99 km, SAXJ1808−3658(SSI) with mass 1.435M⊙ and radius 7.07 km, 4U1608−52 with mass 1.74M⊙ and radius 9.30 km, 4U1538−52 with mass 0.87M⊙ and radius 7.86 km and PSRJ1416−2230 with mass 1.97M⊙ and radius 10.30 km. The physical validity of the obtained solution is verified by computing the necessary physical parameters like energy density and pressure, quintessence density, energy conditions, sound speed via the Herrera cracking concept, hydrostatic equilibrium of forces, mass function, compactness, Buchdahl limit, and surface redshift and analyze their behavior graphically. To investigate the demeanor of these parameters more closely, we computed the numerical values and manifested them in tabular form. We conclude that our presented mathematical model of compact stars in the Finch–Skea geometry with quintessence field fulfills all the requirements for a physically viable solution.

Suggested magnetic braking prescription derived from field complexity fails to reproduce the cataclysmic variable orbital period gap

Magnetic wind braking drives the spin-down of low-mass stars and the evolution of most interacting binary stars. A magnetic braking prescription that was claimed to reproduce both the period distribution of cataclysmic variables (CVs) and the evolution of the rotation rates of low-mass stars is based on a relation between the angular momentum loss rate and magnetic field complexity. The magnetic braking model based on field complexity has been claimed to predict a detached phase that could explain the observed period gap in the period distribution of but has never been tested in detailed models of evolution. Here we fill this gap. We incorporated the suggested magnetic braking law in MESA and simulated the evolution of for different initial stellar masses and initial orbital periods. We find that the prescription for magnetic braking based on field complexity fails to reproduce observations of The predicted secondary star radii are smaller than measured, and an extended detached phase that is required to explain the observed period gap (a dearth of non-magnetic with periods sim 2$ and $ sim 3$ hours ) is not predicted. Proposed magnetic braking prescriptions based on a relation between the angular momentum loss rate and field complexity are too weak to reproduce the bloating of donor stars in derived from observations and, in contrast to previous claims, do not provide an explanation for the observed period gap. The suggested steep decrease in the angular momentum loss rate does not lead to detachment. Stronger magnetic braking prescriptions and a discontinuity at the fully convective boundary are needed to explain the evolution of close binary stars that contain compact objects. The tension between braking laws derived from the spin-down of single stars and those required to explain and other close binaries containing compact objects remains.

Gravitational Wave Constraints on Planetary-Mass Primordial Black Holes Using LIGO O3a Data.

Gravitational waves from subsolar mass inspiraling compact objects would provide almost smoking-gun evidence for primordial black holes (PBHs). We perform the first search for inspiraling planetary-mass compact objects in equal-mass and highly asymmetric mass-ratio binaries using data from the first half of the LIGO-Virgo-KAGRA third observing run. Though we do not find any significant candidates, we determine the maximum luminosity distance reachable with our search to be of O(0.1-100) kpc, and corresponding model-independent upper limits on the merger rate densities to be O(10^{3}-10^{-7}) kpc^{-3} yr^{-1} for systems with chirp masses of O(10^{-4}-10^{-2})M_{⊙}, respectively. Furthermore, we interpret these rate densities as arising from PBH binaries and constrain the fraction of dark matter that such objects could comprise. For equal-mass PBH binaries, we find that these objects would compose less than 4%-100% of DM for PBH masses of 10^{-2}M_{⊙} to 2×10^{-3}M_{⊙}, respectively. For asymmetric binaries, assuming one black hole mass corresponds to a peak in the mass function at 2.5M_{⊙}, a PBH dark-matter fraction of 10% and a second, much lighter PBH, we constrain the mass function of the second PBH to be less than 1 for masses between 1.5×10^{-5}M_{⊙} and 2×10^{-4}M_{⊙}. Our constraints, recently released, are robust enough to be applied to any PBH or exotic compact object binary formation models, and complement existence microlensing results.

The abundance of clustered primordial black holes from quasar microlensing

While elementary particles are the favored candidate for the elusive dark matter, primordial black holes (PBHs) have also been considered to fill that role. Gravitational microlensing is a very well-suited tool to detect and measure the abundance of compact objects in galaxies. Previous studies based on quasar microlensing exclude a significant presence of substellar to intermediate-mass black holes (BHs; $ 100M_ However, these studies were based on a spatially uniform distribution of BHs while, according to current theories of PBH formation, they are expected to appear in clusters. We study the impact of clustering in microlensing flux magnification, finding that at large scales clusters act like giant pseudo-particles, strongly affecting the emission coming from the broad-line region, which can no longer be used to define the zero microlensing baseline. As an alternative, we set this baseline from the intrinsic magnification ratios of quasar images predicted by macro lens models and compared them with the observed flux ratios in emission lines, infrared, and radio. The (magnitude) differences are the flux-ratio anomalies attributable to microlensing, which we estimate for 35 image pairs corresponding to 12 lens systems. A Bayesian analysis indicates that the observed anomalies are incompatible with the existence of a significant population of clustered PBHs. Furthermore, we find that more compact clusters exhibit a stronger microlensing impact. Consequently, we conclude that clustering makes the existence of a significant population of BHs in the substellar to intermediate mass range even more unlikely.