Search

In this study, we investigate gravitational lensing within the framework of more realistic dark matter halo models, transcending the limitations of spherical-collapse approximations. Through analytical computations utilizing diverse mass functions, we address critical factors typically overlooked in the standard Press-Schechter formalism, including ellipsoidal-collapse conditions, angular momentum dynamics, dynamical friction, and the cosmological constant. Our analysis incorporates two widely recognized halo density profiles, the Navarro-Frenk-White and Einasto profiles, considering both spherical and ellipsoidal-collapse scenarios. We present relevant calculations of pivotal gravitational lensing observables, such as Einstein radii, lensing optical depths, and time delays, spanning a wide range of redshifts and masses across two distinct lensing models: the point mass and singular isothermal sphere (SIS) lens models. Our findings illuminate that adopting more realistic dark matter halo models leads to heightened lensing effects compared to their spherical-collapse counterparts. Furthermore, our analyses of lensing optical depths and time delays reveal distinct characteristics between point mass and SIS lens models. These outcomes highlight the need for more realistic halo descriptions instead of simple approximations when modeling gravitational lensing, as this approach can potentially better reveal the complex structures of dark matter., Comment: 13 pages; 5 figures; Refs included

In this work we investigate the merger rate of primordial black hole-neutron star (PBH-NS) binaries in two widely-studied modified gravity (MG) models: Hu-Sawicki $f(R)$ gravity and the normal branch of Dvali-Gabadadze-Porrati (nDGP) gravity. In our analysis, we take into account the effects of MG on the halo properties including halo mass function, halo concentration parameter, halo density profile, and velocity dispersion of dark matter particles. We find that these MG models, due to their stronger gravitational field induced by an effective fifth force, predict enhanced merger rates compared to general relativity. This enhancement is found to be redshift-dependent and sensitive to model parameters, PBH mass and fraction. Assuming PBH mass range of $5-50 M_{\odot}$, we compare the predicted merger rate of PBH-NS binaries with those inferred from LIGO-Virgo-KAGRA observations of gravitational waves (GWs). We find that the merger rates obtained from MG models will be consistent with the GW observations, if the abundance of PBHs is relatively large, with the exact amount depending on the MG model and its parameter values, as well as PBH mass. We also establish upper limits on the abundance of PBHs in these MG frameworks while comparing with the existing non-GW constraints, which can potentially impose even more stringent constraints., Comment: 15 pages, 6 figures

Ultradense dark matter halos (UDMHs) are high concentrations of dark matter, assumed to have formed deep in the radiation-dominated era from amplified primordial perturbations. In this work we improve the previous works for the calculation of UDMH abundance by elaborating on the formation process of these halos by including various physical and geometrical modifications in the analysis. In particular, we investigate the impact of angular momentum, dynamical friction and triaxial collapse on the predicted mass functions for UDMHs. We perform the calculations for four primordial power spectra with different amplified features that allow for primordial black hole and UDHM formation in wide and narrow mass ranges. We also apply this analysis in the context of two possible scenarios for dark matter: the single-component and the multi-component. We find that the abundance of UDMHs is prominently enhanced in the presence of these more realistic mass functions., Comment: 11 pages. 6 figures

Primordial black holes (PBHs) are known as one of the potential candidates for dark matter. They are expected to have formed due to the direct gravitational collapse of density fluctuations in the early Universe. Therefore, the study of the merger rate of PBHs in modified theories of gravity can provide more detailed information about their abundance. In this work, we delve into the calculation of the merger rate of PBHs within the theoretical framework of $f(R)$ gravity. Our analysis reveals an enhancement in the merger rate of PBHs compared to that obtained from general relativity (GR). Additionally, modulating the field strength $f_{R0}$ induces shifts in the PBH merger rate, presenting a potential observational signature of modified gravity. We also find that the total merger rate of PBHs will be consistent with the merger rate of black holes estimated by the Laser Interferometer Gravitational-Wave Observatory (LIGO)-Virgo-KAGRA detectors if $f_{PBH}\gtrsim 0.1$. While further improvements might be required, relative enhancement of the merger rate of PBHs in the framework of $f(R)$ gravity and its consistency with gravitational wave data underscore the importance of employing modified theories of gravity to examine diverse scenarios related to the formation of black holes., Comment: 17 pages, 11 figures, references added

In this work, we extend the model proposed by White concerning the post-collapse evolution of density peaks while considering the role of angular momentum. On a timescale smaller than the peak collapse, $t_{0}$, the inner regions of the peak reach the equilibrium forming a cuspy profile, as in White's paper, but the power-law density profile is flatter, namely $\rho \propto r^{-1.52}$, using the specific angular momentum $J$ obtained in theoretical models of how it evolves in CDM universes, namely $J \propto M^{2/3}$. The previous result shows how angular momentum influences the slope of the density profile, and how a slightly flatter profile obtained in high-resolution numerical simulations, namely $\rho \propto r^{\alpha}$, $(\alpha \simeq -1.5)$ can be reobtained. Similarly to simulations, in our model adiabatic contraction was not taken into account. This means that more comprehensive simulations could give different values for the slope of the density profile, similar to an improvement of our model., Comment: 7 pages, 1 figure, refs. added

Nowadays, the existence of supermassive black holes (SMBHs) in the center of galactic halos is almost confirmed. An extremely dense region referred to as dark-matter spike is expected to form around central SMBHs as they grow and evolve adiabatically. In this work, we calculate the merger rate of compact binaries in dark-matter spikes while considering halo models with spherical and ellipsoidal collapses. Our findings exhibit that ellipsoidal-collapse dark matter halo models can potentially yield the enhancement of the merger rate of compact binaries. Finally, our results confirm that the merger rate of primordial black hole binaries is consistent with the results estimated by the LIGO-Virgo detectors, while such results can not be realized for primordial black hole-neutron star binaries., Comment: 13 pages; 5 figures; references added, typos were fixed

In this work, we study the effect of a high-precision semianalytical mass function on the merger rate of primordial black holes (PBHs) in dark matter halos. For this purpose, we first explain a theoretical framework for dark matter halo models and introduce relevant quantities such as halo density profile, concentration parameter, and a high-precision semianalytical function namely Del Popolo (DP) mass function. In the following, we calculate the merger rate of PBHs in the framework of ellipsoidal-collapse dark matter halo models while considering the DP mass function, and compare it with our previous study for the Sheth-Tormen (ST) mass function. The results show that by taking the mass of PBHs as $M_{\rm PBH} = 30M_{\odot}$, the DP mass function can potentially amplify the merger rate of PBHs. Moreover, we calculate the merger rate of PBHs for the DP mass function as a function of their mass and fraction and compare it with the black hole mergers recorded by the LIGO-Virgo detectors during the latest observing run. Our findings show that the merger rate of PBHs will fall within the LIGO-Virgo band if $f_{\rm PBH} \gtrsim \mathcal{O}(10^{-1})$. This implies that the DP mass function can be used to strengthen constraints on the fraction of PBHs., Comment: 11 pages, 4 figures, 1 table, accepted version

Cosmic voids are known as underdense substructures of the cosmic web that cover a large volume of the Universe. It is known that cosmic voids contain a small number of dark matter halos, so the existence of primordial black holes (PBHs) in these secluded regions of the Universe is not unlikely. In this work, we calculate the merger rate of PBHs in dark matter halos structured in cosmic voids and determine their contribution to gravitational wave events resulting from black hole mergers recorded by the Advanced Laser Interferometer Gravitational-Wave Observatory (aLIGO)-Advanced Virgo (aVirgo) detectors. Relying on the PBH scenario, the results of our analysis indicate that about $2 \sim 3$ annual events of binary black hole mergers out of all those recorded by the aLIGO-aVirgo detectors should belong to cosmic voids. We also calculate the redshift evolution of the merger rate of PBHs in cosmic voids. The results show that the evolution of the merger rate of PBHs has minimum sensitivity to the redshift changes, which seems reasonable while considering the evolution of cosmic voids. Finally, we specify the behavior of the merger rate of PBHs as a function of their mass and fraction in cosmic voids and we estimate $\mathcal{R} (M_{PBH}, f_{PBH})$ relation, which is well compatible with our findings., Comment: 11 pages; 4 figures; References are included. Physics of the Dark Universe (2023)

Primordial black holes (PBHs), as a potential macroscopic candidate for dark matter, can encounter other compact objects in dark matter halos because of their random distribution. Besides, the detection of gravitational waves (GWs) related to the stellar-mass black hole-neutron star (BH-NS) mergers raises the possibility that the BHs involved in such events may have a primordial origin. In this work, we calculate the merger rate of PBH-NS binaries within the framework of ellipsoidal-collapse dark matter halo models and compare it with the corresponding results derived from spherical-collapse dark matter halo models. Our results exhibit that ellipsoidal-collapse dark matter halo models can potentially amplify the merger rate of PBH-NS binaries in such a way that it is very close to the range estimated by the LIGO-Virgo observations. While spherical-collapse dark matter halo models cannot justify PBH-NS merger events as consistent results with the latest GW data reported by the LIGO-Virgo collaboration. In addition, we calculate the merger rate of PBH-NS binaries as a function of PBH mass and fraction within the context of ellipsoidal-collapse dark matter halo models. The results indicate that PBH-NS merger events with the mass of $(M_{PBH}\le 5 M_{\odot}, M_{NS}\simeq 1.4 M_{\odot})$ will be consistent with the LIGO-Virgo observations if $f_{PBH}\simeq 1$. We also show that to have at least on $(M_{PBH}\simeq 5 M_{\odot}, M_{NS}\simeq 1.4 M_{\odot})$ event in the comoving volume $1 Gpc^{3}$ annually, ellipsoidal-collapse dark matter halo models constrain the abundance of PBHs as $f_{PBH} \geq 0.1$., Comment: 10 pages, 4 figures, 1 table, references added, accepted in The Astrophysical Journal (APJ)

Interacting dark matter-dark energy (IDMDE) models can be taken to account as one of the present challenges that may affect the cosmic structures. In this work, we study the integrated Sachs-Wolfe (ISW) effect in IDMDE models. To this end, we initially introduce a theoretical framework for IDMDE models. Moreover, we briefly discuss the stability conditions of IDMDE models and by specifying a simple functional form for the energy density transfer rate, we calculate the perturbation equations. In the following, we calculate the amplitude of the matter power spectrum for the IDMDE model and compare it with the corresponding result obtained from the $\Lambda$CDM model. Furthermore, we calculate the amplitude of the ISW auto-power spectrum as a function of multipole order l for the IDMDE model. The results indicate that the amplitude of the ISW auto-power spectrum in the IDMDE model for different phantom dark energy equations of state behaves similar to the one for the $\Lambda$CDM model, whereas, for the quintessence dark energy equations of state, the amplitude of the ISW-auto power spectrum for the IDMDE model should be higher than the one for the $\Lambda$CDM model. Also, it turns out that the corresponding results by different values of the coupling parameter demonstrate that $\xi$ is inversely proportional to the amplitude of the ISW-auto power spectrum in the IDMDE model. Finally, by employing four different surveys, we calculate the amplitude of the ISW-cross power spectrum as a function of multipole order $l$ for the IDMDE model. The results exhibit that the amplitude of the ISW-cross power spectrum for the IDMDE model for all values of $\omega_{\rm x}$ is higher than the one obtained for the $\Lambda$CDM model. Also, it turns out that the amplitude of the ISW-cross power spectrum in the IDMDE model changes inversely with the value of coupling parameter $\xi$., Comment: 13 pages, 2 tables, 8 figures, references added

We study the cosmological inflation within the context of f(Q, T) gravity, wherein Q is the nonmetricity scalar and T is the trace of the matter energy-momentum tensor. By choosing a linear combination of Q and T, we first analyze the realization of an inflationary scenario driven via the geometrical effects of the linear f(Q, T) gravity and then, we obtain the modified slow-roll parameters, the scalar and the tensor spectral indices, and the tensor-to-scalar ratio for the proposed model. In addition, by choosing three inflationary potentials, i.e. the power-law, hyperbolic and natural potentials, and by applying the slow-roll approximations, we calculate these inflationary observables in the presence of an inflaton scalar field. The results indicate that by properly restricting the free parameters, the proposed model provides appropriate predictions that are consistent with the observational data obtained from the Planck 2018. Also, we specify that the contribution of linear model of f(Q, T) gravity with the hyperbolic and natural potentials can impose different restrictions on the parameters of these potentials. Furthermore, the predictions of natural inflation in this model are in good agreement with the joint Planck, BK15 and BAO data, justifying the use of the f(Q, T) gravity., Comment: revised 18 pages, a few figures, references added

A novel 4-dimensional Einstein-Gauss-Bonnet (4D EGB) gravity has been proposed that asserts to bypass the Lovelock's theorem and to result in a non-trivial contribution to the gravitational dynamics in four-dimensional spacetime. In this work, we study the integrated Sachs-Wolfe (ISW) effect in the 4D EGB model. For this purpose, we calculate the evolution of the gravitational potential and the linear growth factor as a function of redshift for the 4D EGB model and compare it with the corresponding result obtained from the $\Lambda$-cold dark matter ($\Lambda$CDM) model. We also calculate the ISW-auto power spectrum and the ISW-cross power spectrum as functions of cosmic microwave background (CMB) multipoles for the 4D EGB model and compare those with the one obtained from the $\Lambda$CDM model. To do this, we use the strongest constraint on the coupling parameter proposed for the 4D EGB model. Additionally, to calculate the ISW effect for the 4D EGB model, we employ three large-scale structure surveys from different wavelengths. The results exhibit that the ISW effect in the 4D EGB model is higher than the one obtained from the $\Lambda$CDM model. Hence, we show that the 4D EGB model can amplify the amplitude of the ISW power spectrum, which can be considered as a relative advantage of the 4D EGB model comparing the $\Lambda$CDM one. Also, we indicate that the deviation from the $\Lambda$CDM model is directly proportional to the value of the dimensionless coupling parameter $\beta$., Comment: 10 pages, 2 tables, fixed some typos, added some reference

We study the merger rate of primordial black holes (PBHs) in self-interacting dark matter (SIDM) halo models. To explore a numerical description for the density profile of SIDM halo models, we use the result of a previously performed simulation for SIDM halo models with $\sigma/m=10~{\rm cm^{2}g^{-1}}$. We also propose a concentration-mass-time relation that can explain the evolution of the halo density profile related to SIDM models. Furthermore, we investigate the encounter condition of PBHs that may have been randomly distributed in the medium of dark matter halos. Under these assumptions, we calculate the merger rate of PBHs within each halo considering SIDM halo models and compare the results with that obtained for cold dark matter (CDM) halo models. To do this, we employ the definition of the time after halo virialization as a function of halo mass. We indicate that SIDM halo models for $f_{\rm PBH}>0.32$ can generate sufficient PBH mergers in such a way that those exceed the one resulted from CDM halo models. By considering the spherical-collapse halo mass function, we obtain similar results for the cumulative merger rate of PBHs. Moreover, we calculate the redshift evolution of the PBH total merger rate. To determine a constraint on the PBH abundance, we study the merger rate of PBHs in terms of their fraction and masses and compare those with the black hole merger rate estimated by the Advanced LIGO (aLIGO)-Advanced Virgo (aVirgo) detectors during the third observing run. The results demonstrate that within the context of SIDM halo models, the merger rate of $10~M_{\odot}-10~M_{\odot}$ events can potentially fall within the aLIGO-aVirgo window. We also estimate a relation between the fraction of PBHs and their masses, which is well consistent with our findings., Comment: 17 pages, 9 figures, refs. added, fixed some typos

We study numerically the nonlinear stability of {\it excited} fermion-boson stars in spherical symmetry. Such compound hypothetical stars, composed by fermions and bosons, are gravitationally bound, regular, and static configurations described within the coupled Einstein-Klein-Gordon-Euler theoretical framework. The excited configurations are characterized by the presence in the radial profile of the (complex, massive) scalar field -- the bosonic piece -- of at least one node across the star. The dynamical emergence of one such configuration from the accretion of a cloud of scalar field onto an already-formed neutron star, was numerically revealed in our previous investigation. Prompted by that finding we construct here equilibrium configurations of excited fermion-boson stars and study their stability properties using numerical-relativity simulations. In addition, we also analyze their dynamical formation from generic, constraint-satisfying initial data. Contrary to purely boson stars in the excited state, which are known to be generically unstable, our study reveals the appearance of a cooperative stabilization mechanism between the fermionic and bosonic constituents of those excited-state mixed stars. While similar examples of stabilization mechanisms have been recently discussed in the context of $\ell-$boson stars and multi-field, multi-frequency boson stars, our results seem to indicate that the stabilization mechanism is a purely gravitational effect and does not depend on the type of matter of the companion star., Comment: 11 pages, 16 figures

We have studied the merger rate of primordial black holes (PBHs) in the ellipsoidal-collapse model of halo to explain the dark matter abundance by the PBH merger estimated from the gravitational waves detections via the Advanced LIGO (aLIGO) detectors. We have indicated that the PBH merger rate within each halo for the ellipsoidal models is more significant than for the spherical models. We have specified that the PBH merger rate per unit time and per unit volume for the ellipsoidal-collapse halo models is about one order of magnitude higher than the corresponding spherical models. Moreover, we have calculated the evolution of the PBH total merger rate as a function of redshift. The results indicate that the evolution for the ellipsoidal halo models is more sensitive than spherical halo models, as expected from the models. Finally, we have presented a constraint on the PBH abundance within the context of ellipsoidal and spherical models. By comparing the results with the aLIGO mergers during the third observing run (O3), we have shown that the merger rate in the ellipsoidal-collapse halo models falls within the aLIGO window, while the same result is not valid for the spherical-collapse ones. Furthermore, we have compared the total merger rate of PBHs in terms of their fraction in the ellipsoidal-collapse halo models for several masses of PBHs. The results suggest that the total merger rate of PBHs changes inversely with their masses. We have also estimated the relation between the fraction of PBHs and their masses in the ellipsoidal-collapse halo model and have shown it for a narrow mass distribution of PBHs. The outcome shows that the constraint inferred from the PBH merger rate for the ellipsoidal-collapse halo models can be potentially stronger than the corresponding result obtained for the spherical-collapse ones., Comment: 13 pages, 6 figures, refs. added

We obtain analytical approximate black hole solutions for higher derivative gravity in the presence of Maxwell electromagnetic source. We construct near horizon and asymptotic solutions and then use these to obtain an approximate analytic solution using a continued fraction method to get a complete solution. We compute the thermodynamic quantities and check the first law and Smarr formula. Finally, we investigate the null and time-like geodesics of this black hole., Comment: 17 pages, 47 figures, Accepted in PRD

Gravitationally bound structures composed by fermions and scalar particles known as fermion-boson stars are regular and static configurations obtained by solving the coupled Einstein-Klein-Gordon-Euler (EKGE) system. In this work, we discuss one possible scenario through which these fermion-boson stars may form by solving numerically the EKGE system under the simplifying assumption of spherical symmetry. Our initial configurations assume an already existing neutron star surrounded by an accreting cloud of a massive and complex scalar field. The results of our simulations show that once part of the initial scalar field is expelled via gravitational cooling the system gradually oscillates around an equilibrium configuration that is asymptotically consistent with a static solution of the system. The formation of fermion-boson stars for large positive values of the coupling constant in the self-interaction term of the scalar-field potential reveal the presence of a node in the scalar field. This suggests that a fermionic core may help stabilize configurations with nodes in the bosonic sector, as happens for purely boson stars in which the ground state and the first excited state coexist., Comment: 13 pages, 16 figures