Matter Density Parameter Research Articles

Testing the standardizability of, and deriving cosmological constraints from, a new Amati-correlated gamma-ray burst data compilation

By using gamma-ray burst (GRB) data to simultaneously constrain Amati correlation parameters and cosmological parameters in six spatially flat and nonflat dark energy cosmological models, we show that an updated 220 GRB version of the Jia et al. [1] GRB data compilation are standardizable through the Amati correlation and so can be used for cosmological analyses. However, the resulting GRB data constraints on the current value of the nonrelativistic matter density parameter, Ω m 0, are in > 2σ tension with those from a joint analysis of better-established Hubble parameter [H(z)] and baryon acoustic oscillation (BAO) data for most of the cosmological models we consider, indicating that these GRB data cannot be jointly used with better-established H(z) + BAO data to constrain cosmological parameters.

Late-time constraints on interacting dark energy: Analysis independent of H0, rd, and MB

We investigated a possible interaction between cold dark matter and dark energy, corresponding to a well-known interacting dark energy model discussed in the literature within the context of resolving the Hubble tension. We put constraints on it in a novel way, by creating new likelihoods with an analytical marginalization over the Hubble parameter H0, the sound horizon rd, and the supernova absolute magnitude MB. Our aim is to investigate the impacts on the coupling parameter of the interacting model, ξ, and the equation of state of dark energy w and the matter density parameter Ωm, 0. The late-time cosmological probes used in our analysis include the PantheonPlus (calibrated and uncalibrated), cosmic chronometers, and baryon acoustic oscillation samples and the Pantheon for comparison. Through various combinations of these datasets, we demonstrate hints of an up to 2σ deviation from the standard Λ cold dark matter model.

Prospects for weighing neutrinos in interacting dark energy models using joint observations of gravitational waves and γ-ray bursts* *Supported by the National Natural Science Foundation of China (12305069, 11947022, 11975072, 11875102, 11835009), the National SKA Program of China (2022SKA0110200, 2022SKA0110203), the National 111 Project (B16009), and the Program of the Education Department of Liaoning Province, China (JYTMS20231695)

Cosmological observations can be used to weigh neutrinos, but this method is model-dependent, with results relying on the cosmological model considered. If we consider interactions between dark energy and dark matter, the neutrino mass constraints differ from those derived under the standard model. On the contrary, gravitational wave (GW) standard siren observations can measure absolute cosmological distances, helping to break parameter degeneracies inherent in traditional cosmological observations, thereby improving constraints on neutrino mass. This paper examines the constraints on neutrino mass within interacting dark energy (IDE) models and explores how future GW standard siren observations could enhance these results. For multi-messenger GW observations, we consider the joint observations of binary neutron star mergers by third-generation ground-based GW detectors and short γ-ray burst observations by missions similar to the THESEUS satellite project. Using current cosmological observations (CMB+BAO+SN), we obtain an upper limit on the neutrino mass in the IDE models of 0.15 (or 0.16) eV. With the inclusion of GW data, the upper limit on the neutrino mass improves to 0.14 eV. This indicates that in the context of IDE models, the improvement in neutrino mass constraints from GW observations is relatively limited. However, GW observations significantly enhance the constraints on other cosmological parameters, such as matter density parameter, the Hubble constant, and coupling strength between dark energy and dark matter.

KiDS-1000 cosmology: Combined second- and third-order shear statistics

Aims. In this work, we perform the first cosmological parameter analysis of the fourth release of Kilo Degree Survey (KiDS-1000) data with second- and third-order shear statistics. This paper builds on a series of studies aimed at describing the roadmap to third-order shear statistics. Methods. We derived and tested a combined model of the second-order shear statistic, namely, the COSEBIs and the third-order aperture mass statistics 〈ℳap3〉 in a tomographic set-up. We validated our pipeline with N-body mock simulations of the KiDS-1000 data release. To model the second- and third-order statistics, we used the latest version of HMCODE2020 for the power spectrum and BIHALOFIT for the bispectrum. Furthermore, we used an analytic description to model intrinsic alignments and hydro-dynamical simulations to model the effect of baryonic feedback processes. Lastly, we decreased the dimension of the data vector significantly by considering only equal smoothing radii for the 〈ℳap3〉 part of the data vector. This makes it possible to carry out a data analysis of the KiDS-1000 data release using a combined analysis of COSEBIs and third-order shear statistics. Results. We first validated the accuracy of our modelling by analysing a noise-free mock data vector, assuming the KiDS-1000 error budget, finding a shift in the maximum of the posterior distribution of the matter density parameter, ΔΩm < 0.02 σΩm, and of the structure growth parameter, ΔS8 < 0.05 σS8. Lastly, we performed the first KiDS-1000 cosmological analysis using a combined analysis of second- and third-order shear statistics, where we constrained Ωm = 0.248−0.055+0.062 and S8 = σ8√(Ωm/0.3 )= 0.772 ± 0.022. The geometric average on the errors of Ωm and S8 of the combined statistics decreases, compared to the second-order statistic, by a factor of 2.2.

A new binning method to choose a standard set of Quasars

Although the Lambda Cold Dark Matter model is the most accredited cosmological model, information at intermediate redshifts (z) between type Ia Supernovae (z = 2.26) and the Cosmic Microwave Background (z = 1100) is crucial to validate this model further. Here, we present a detailed and reliable methodology for binning the quasars (QSO) data that allows the identification of a golden sample of QSOs to be used as standard candles. This procedure has the advantage of being very general. Thus, it can be applied to any astrophysical sources at cosmological distances. This methodology allows us to avoid the circularity problem since it involves a flux–flux relation and includes the analysis of removing selection biases and the redshift evolution. With this method, we have discovered a sample of 1253 quasars up to z = 7.54 with reduced intrinsic dispersion of the relation between Ultraviolet and X-ray fluxes, with δint=0.096±0.003 (56% less than the original sample where δint=0.22). Once the luminosities are corrected for selection biases and redshift evolution, this ‘gold’ sample allows us to determine the matter density parameter to be ΩM=0.240±0.064. This value is aligned with the results of the ΛCDM model obtained with SNe Ia.

The Scavenger Hunt for Quasar Samples to Be Used as Cosmological Tools

Although the Λ Cold Dark Matter model is the most accredited cosmological model, information at high redshifts (z) between type Ia supernovae (z=2.26) and the Cosmic Microwave Background (z=1100) is crucial to validate this model further. To this end, we have discovered a sample of 1132 quasars up to z=7.54 exhibiting a reduced intrinsic dispersion of the relation between ultraviolet and X-ray fluxes, δF=0.22 vs. δF=0.29 (24% less), than the original sample. This gold sample, once we correct the luminosities for selection biases and redshift evolution, enables us to determine the matter density parameter ΩM with a precision of 0.09. Unprecedentedly, this quasar sample is the only one that, as a standalone cosmological probe, yields such tight constraints on ΩM while being drawn from the same parent population of the initial sample.

Constraining cosmology with the Gaia-unWISE Quasar Catalog and CMB lensing: structure growth

We study the angular clustering of Quaia, a Gaia- and unWISE-based catalog of over a million quasars with an exceptionally well-defined selection function. With it, we derive cosmology constraints from the amplitude and growth of structure across cosmic time. We divide the sample into two redshift bins, centered at z = 1.0 and z = 2.1, and measure both overdensity auto-correlations and cross-correlations with maps of the Cosmic Microwave Background convergence measured by Planck. From these data, and including a prior from measurements of the baryon acoustic oscillations scale, we place constraints on the amplitude of the matter power spectrum σ 8 = 0.766 ± 0.034, and on the matter density parameter Ω m = 0.343+0.017 -0.019. These measurements are in reasonable agreement with Planck at the ∼ 1.4σ level, and are found to be robust with respect to observational and theoretical uncertainties. We find that our slightly lower value of σ 8 is driven by the higher-redshift sample, which favours a low amplitude of matter fluctuations. We present plausible arguments showing that this could be driven by contamination of the CMB lensing map by high-redshift extragalactic foregrounds, which should also affect other cross-correlations with tracers of large-scale structure beyond z ∼ 1.5. Our constraints are competitive with those from state-of-the-art 3×2-point analyses, but arise from a range of scales and redshifts that is highly complementary to those covered by cosmic shear data and most galaxy clustering samples. This, coupled with the unprecedented combination of volume and redshift precision achieved by Quaia, allows us to break the usual degeneracy between Ω m and σ 8.

Euclid: Cosmology forecasts from the void-galaxy cross-correlation function with reconstruction

We have investigated the cosmological constraints that can be expected from measurement of the cross-correlation of galaxies with cosmic voids identified in the Euclid spectroscopic survey, which will include spectroscopic information for tens of millions of galaxies over 15 000 deg2 of the sky in the redshift range 0.9 ≤ z < 1.8. We have done this using simulated measurements obtained from the Flagship mock catalogue, the official Euclid mock that closely matches the expected properties of the spectroscopic dataset. To mitigate anisotropic selection-bias effects, we have used a velocity field reconstruction method to remove large-scale redshift-space distortions from the galaxy field before void-finding. This allowed us to accurately model contributions to the observed anisotropy of the cross-correlation function arising from galaxy velocities around voids as well as from the Alcock–Paczynski effect, and we studied the dependence of constraints on the efficiency of reconstruction. We find that Euclid voids will be able to constrain the ratio of the transverse comoving distance DM and Hubble distance DH to a relative precision of about 0.3%, and the growth rate fσ8 to a precision of between 5% and 8% in each of the four redshift bins covering the full redshift range. In the standard cosmological model, this translates to a statistical uncertainty ΔΩm = ±0.0028 on the matter density parameter from voids, which is better than what can be achieved from either Euclid galaxy clustering and weak lensing individually. We also find that voids alone can measure the dark energy equation of state to a 6% precision.

Cosmological Constraints from the BOSS DR12 Void Size Function

We present the first cosmological constraints derived from the analysis of the void size function. This work relies on the final Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 12 (DR12) data set, a large spectroscopic galaxy catalog, ideal for the identification of cosmic voids. We extract a sample of voids from the distribution of galaxies, and we apply a cleaning procedure aimed at reaching high levels of purity and completeness. We model the void size function by means of an extension of the popular volume-conserving model, based on two additional nuisance parameters. Relying on mock catalogs specifically designed to reproduce the BOSS DR12 galaxy sample, we calibrate the extended size function model parameters and validate the methodology. We then apply a Bayesian analysis to constrain the Lambda cold dark matter (ΛCDM) model and one of its simplest extensions, featuring a constant dark energy equation of state parameter, w. Following a conservative approach, we put constraints on the total matter density parameter and the amplitude of density fluctuations, finding Ωm = 0.29 ± 0.06 and . Testing the alternative scenario, we derive w = −1.1 ± 0.2, in agreement with the ΛCDM model. These results are independent and complementary to those derived from standard cosmological probes, opening up new ways to identify the origin of potential tensions in the current cosmological paradigm.

Observational constraints on Yukawa cosmology and connection with black hole shadows

We confront Yukawa modified cosmology, proposed in Jusufi et al. [arXiv:2304.11492], with data from Supernovae Type Ia (SNe Ia) and Hubble parameter (OHD) observations. Yukawa cosmology is obtained from a Yukawa-like gravitational potential, with coupling parameter α and wavelength parameter λ, which gives rise to modified Friedmann equations. We show that the agreement with observations is very efficient, and within 1σ confidence level we find the best-fit parameters λ=2693−1262+1191Mpc, α=0.416−0.326+1.137, and a graviton mass of mg=2.374−0.728+2.095×10−42GeV. Additionally, we establish a connection between the effective dark matter and dark energy density parameters and the angular radius of the black hole shadow of the SgrA and M87 black holes in the low-redshift limit, consistent with the Event Horizon Telescope findings.

KiDS-1000: Combined halo-model cosmology constraints from galaxy abundance, galaxy clustering, and galaxy-galaxy lensing

We present constraints on the flat Λ cold dark matter cosmological model through a joint analysis of galaxy abundance, galaxy clustering, and galaxy-galaxy lensing observables with the Kilo-Degree Survey. Our theoretical model combines a flexible conditional stellar mass function, which describes the galaxy-halo connection, with a cosmological N-body simulation-calibrated halo model, which describes the non-linear matter field. Our magnitude-limited bright galaxy sample combines nine-band optical-to-near-infrared photometry with an extensive and complete spectroscopic training sample to provide accurate redshift and stellar mass estimates. Our faint galaxy sample provides a background of accurately calibrated lensing measurements. We constrain the structure growth parameter to S8 = σ8√Ωm/0.3 =√0.773−0.030+0.028 and the matter density parameter to Ωm = 0.290−0.017+0.021. The galaxy-halo connection model adopted in the work is shown to be in agreement with previous studies. Our constraints on cosmological parameters are comparable to, and consistent with, joint ‘3 × 2pt’ clustering-lensing analyses that additionally include a cosmic shear observable. This analysis therefore brings attention to the significant constraining power in the often excluded non-linear scales for galaxy clustering and galaxy-galaxy lensing observables. By adopting a theoretical model that accounts for non-linear halo bias, halo exclusion, scale-dependent galaxy bias, and the impact of baryon feedback, this work demonstrates the potential for, and a way towards, including non-linear scales in cosmological analyses. Varying the width of the satellite galaxy distribution with an additional parameter yields a strong preference for sub-Poissonian variance, improving the goodness of fit by 0.18 in terms of the reduced χ2 value (and increasing the p-value by 0.25) compared to a fixed Poisson distribution.

Quasars: Standard Candles up to z = 7.5 with the Precision of Supernovae Ia

Currently, the Λ cold dark matter model, which relies on the existence of cold dark matter and a cosmological constant Λ, best describes the universe. However, we lack information in the high-redshift (z) region between Type Ia supernovae (SNe Ia; up to z = 2.26) and the cosmic microwave background (z = 1100), an interval crucial to test cosmological models and their possible evolution. We have defined a sample of 983 quasars up to z = 7.54 with a reduced intrinsic dispersion δ = 0.007, which determines the matter density parameter Ω M with the same precision of SNe Ia. Although previous analysis have used quasars as cosmological tools, this is the first time that high-redshift sources, in this case quasars, as standalone cosmological probes yield such tight constraints on Ω M . Our results show the importance of correcting cosmological relationships for selection biases and redshift evolution and how the choice of a golden sample reduces considerably the intrinsic scatter. This proves the reliability of quasars as standard cosmological candles.

Dark matter profiles of SPARC galaxies: a challenge to fuzzy dark matter

ABSTRACT Stellar and gas kinematics of galaxies are a sensitive probe of the dark matter distribution in the halo. The popular fuzzy dark matter models predict the peculiar shape of density distribution in galaxies: specific dense core with sharp transition to the halo. Moreover, fuzzy dark matter predicts scaling relations between the dark matter particle mass and density parameters. In this work, we use a Bayesian framework and several dark matter halo models to analyse the stellar kinematics of galaxies using the Spitzer Photometry and Accurate Rotation Curves data base. We then employ a Bayesian model comparison to select the best halo density model. We find that more than half of the galaxies prefer the fuzzy dark model against standard dark matter profiles (NFW, Burkert, and cored NFW). While this seems like a success for fuzzy dark matter, we also find that there is no single value for the particle mass that provides a good fit for all galaxies.

H0=69.8±1.3 km s−1 Mpc−1 , Ωm0=0.288±0.017 , and other constraints from lower-redshift, non-CMB, expansion-rate data

We use updated Type Ia Pantheon+ supernova, baryon acoustic oscillation, and Hubble parameter (now also accounting for correlations) data, as well as new reverberation-measured C IV quasar data, and quasar angular size, H II starburst galaxy, reverberation-measured Mg II quasar, and Amati-correlated gamma-ray burst data to constrain cosmological parameters. We show that these data sets result in mutually consistent constraints and jointly use them to constrain cosmological parameters in six different spatially-flat and non-flat cosmological models. Our analysis provides summary model-independent determinations of two key cosmological parameters: the Hubble constant, $H_0=69.8\pm1.3$ $\rm{km \ s^{-1} \ Mpc^{-1}}$, and the current non-relativistic matter density parameter, $\Omega_{m0}=0.288\pm0.017$. Our summary error bars are 2.4 and 2.3 times those obtained using the flat $\Lambda$CDM model and Planck TT,TE,EE+lowE+lensing cosmic microwave background (CMB) anisotropy data. Our $H_0$ value is very consistent with that from the local expansion rate based on the Tip of the Red Giant Branch data, is 2$\sigma$ lower than that from the local expansion rate based on Type Ia supernova and Cepheid data, and is 2$\sigma$ higher than that in the flat $\Lambda$CDM model based on Planck TT,TE,EE+lowE+lensing CMB data. Our data compilation shows at most mild evidence for non-flat spatial hypersurfaces, but more significant evidence for dark energy dynamics, 2$\sigma$ or larger in the spatially-flat dynamical dark energy models we study.

Exact Solutions and Cosmological Constraints in Fractional Cosmology

This paper investigates exact solutions of cosmological interest in fractional cosmology. Given μ, the order of Caputo’s fractional derivative, and w, the matter equation of state, we present specific exact power-law solutions. We discuss the exact general solution of the Riccati Equation, where the solution for the scale factor is a combination of power laws. Using cosmological data, we estimate the free parameters. An analysis of type Ia supernovae (SNe Ia) data and the observational Hubble parameter data (OHD), also known as cosmic chronometers, and a joint analysis with data from SNe Ia + OHD leads to best-fit values for the free parameters calculated at 1σ, 2σ and 3σ confidence levels (CLs). On the other hand, these best-fit values are used to calculate the age of the Universe, the current deceleration parameter (both at 3σ CL) and the current matter density parameter at 1σ CL. Finding a Universe roughly twice as old as the one of ΛCDM is a distinction of fractional cosmology. Focusing our analysis on these results, we can conclude that the region in which μ>2 is not ruled out by observations. This parameter region is relevant because fractional cosmology gives a power-law solution without matter, which is accelerated for μ>2. We present a fractional origin model that leads to an accelerated state without appealing to Λ or dark energy.

Quasar UV/X-ray relation luminosity distances are shorter than reverberation-measured radius–luminosity relation luminosity distances

ABSTRACT We use measurements of 59/58 quasars (QSOs), over a redshift range 0.0041 ≤ z ≤ 1.686, to do a comparative study of the radius–luminosity (R − L) and X-ray−UV luminosity (LX − LUV) relations and the implication of these relations for cosmological parameter estimation. By simultaneously determining R − L or LX − LUV relation parameters and cosmological parameters in six different cosmological models, we find that both R − L and LX − LUV relations are standardizable but provide only weak cosmological parameter constraints, with LX − LUV relation data favouring larger current non-relativistic matter density parameter Ωm0 values than R − L relation data and most other available data. We derive LX − LUV and R − L luminosity distances for each of the sources in the six cosmological models and find that LX − LUV relation luminosity distances are shorter than R − L relation luminosity distances as well as standard flat ΛCDM model luminosity distances. This explains why LX − LUV relation QSO data favour larger Ωm0 values than do R − L relation QSO data or most other cosmological measurements. While our sample size is small and only spans a small z range, these results indicate that more work is needed to determine whether the LX − LUV relation can be used as a cosmological probe.

Starlet higher order statistics for galaxy clustering and weak lensing

We present a first application to photometric galaxy clustering and weak lensing of wavelet-based multi-scale (beyond two points) summary statistics: starlet peak counts and starlet ℓ1-norm. Peak counts are the local maxima in the map, and ℓ1-norm is computed via the sum of the absolute values of the starlet (wavelet) decomposition coefficients of a map, providing a fast multi-scale calculation of the pixel distribution, encoding the information of all pixels in the map. We employ the cosmo-SLICS simulations sources and lens catalogues, and we compute wavelet-based non-Gaussian statistics in the context of combined probes and their potential when applied to the weak-lensing convergence maps and galaxy maps. We obtain forecasts on the matter density parameter Ωm, the reduced Hubble constant h, the matter fluctuation amplitude σ8, and the dark energy equation of state parameter w0. In our setting for this first application, we consider the two probes to be independent. We find that the starlet peaks and the ℓ1-norm represent interesting summary statistics that can improve the constraints with respect to the power spectrum, even in the case of photometric galaxy clustering and when the two probes are combined.

Gamma-ray bursts, quasars, baryonic acoustic oscillations, and supernovae Ia: new statistical insights and cosmological constraints

ABSTRACTThe recent $\sim 4 \, \sigma$ Hubble constant, H0, tension is observed between the value of H0 from the cosmic microwave background (CMB) and Type Ia supernovae (SNe Ia). It is a decade since this tension is excruciating the modern astrophysical community. To shed light on this problem is key to consider probes at intermediate redshifts between SNe Ia and CMB and reduce the uncertainty on H0. Toward these goals, we fill the redshift gap by employing gamma-ray bursts (GRBs) and quasars (QSOs), reaching z = 9.4 and z = 7.6, respectively, combined with baryonic acoustic oscillations (BAO), and SNe Ia. To this end, we employ the ‘Dainotti GRB 3D relation’ among the rest-frame end time of the X-ray plateau emission, its corresponding luminosity, and the peak prompt luminosity, and the ‘Risaliti–Lusso’ QSO relation between ultraviolet and X-ray luminosities. We inquire the commonly adopted Gaussianity assumption on GRBs, QSOs, and BAO. With the joint sample, we fit the flat Λ Cold Dark Matter model with both the Gaussian and the newly discovered likelihoods. We also investigate the impact of the calibration assumed for Pantheon and Pantheon + SNe Ia on this analysis. Remarkably, we show that only GRBs fulfil the Gaussianity assumption. We achieve small uncertainties on the matter-density parameter ΩM and H0. We find H0 values compatible within 2σ with the one from the Tip of the Red Giant Branch. Finally, we show that the cosmological results are heavily biased against the arbitrary calibration choice for SNe Ia.

UNIONS: The impact of systematic errors on weak-lensing peak counts

Context. The Ultraviolet Near-Infrared Optical Northern Survey (UNIONS) is an ongoing deep photometric multiband survey of the northern sky. As part of UNIONS, the Canada-France Imaging Survey (CFIS) provides r-band data, which we use to study weak-lensing peak counts for cosmological inference. Aims. We assess systematic effects for weak-lensing peak counts and their impact on cosmological parameters for the UNIONS survey. In particular, we present results on local calibration, metacalibration shear bias, baryonic feedback, the source galaxy redshift estimate, intrinsic alignment, and cluster member dilution. Methods. For each uncertainty and systematic effect, we describe our mitigation scheme and the impact on cosmological parameter constraints. We obtain constraints on cosmological parameters from Monte Carlo Markov chains using CFIS data and MassiveNuS N-body simulations as a model for peak counts statistics. Results. Depending on the calibration (local versus global, and the inclusion or not of the residual multiplicative shear bias), the mean matter density parameter, Ωm, can shift by up to −0.024 (−0.5σ). We also see that including baryonic corrections can shift Ωm by +0.027 (+0.5σ) with respect to the dark-matter-only simulations. Reducing the impact of the intrinsic alignment and cluster member dilution through signal-to-noise cuts leads to larger constraints. Finally, with a mean redshift uncertainty of Δz̄ = 0.03, we see that the shift in Ωm (+0.001, which corresponds to +0.02σ) is not significant. Conclusions. This paper investigates, for the first time with UNIONS weak-lensing data and peak counts, the impact of systematic effects. The value of Ωm is the most impacted and can shift by up to ∼0.03, which corresponds to 0.5σ depending on the choices for each systematics. We expect constraints to become more reliable with future (larger) data catalogs, for which the current pipeline will provide a starting point. The code used to obtain the results is available on GitHub.

Parameter estimation of eccentric gravitational waves with a decihertz observatory and its cosmological implications

Eccentricity of compact binaries can improve the parameter estimation of gravitational waves (GWs), which is due to the fact that the multiple harmonics induced by eccentricity can provide more information and break the degeneracy between waveform parameters. In this paper, we first investigate the parameter estimation of eccentric GWs with decihertz observatory. We consider two scenarios for the configuration of DECIGO, i.e., the one cluster of DECIGO with its design sensitivity and B-DECIGO which also has one cluster but with inferior sensitivity as a comparison. We adopt the Fisher matrix to estimate the parameter errors. By mocking up the typical binaries in GWTC-3, we find a nonvanishing eccentricity can significantly improve the estimation for almost all waveform parameters. In particular, the localization of typical binary black holes (BBH) can achieve $\mathcal{O}(10--{10}^{3.5})$ factors of improvement when the initial eccentricity ${e}_{0}=0.4$ at 0.1 Hz. The precise localization of binary neutron stars (BNS) and neutron star--black hole binaries (NSBH), together with the large improvement of localization of BBH from eccentricity in the midband, inspire us to construct the catalogs of golden dark sirens whose host galaxies can be uniquely identified. We find that with only one cluster of DECIGO running 1 year in its design sensitivity, hundreds of golden dark BNS, NSBH, and tens of golden dark BBH can be observed. Eccentricity can greatly increase the population of golden dark BBH from $\ensuremath{\sim}7({e}_{0}=0)$ to $\ensuremath{\sim}65({e}_{0}=0.2)$. Such an increase of population of golden dark BBH events can improve the precision of Hubble constant measurement from 2.06% to 0.68%, matter density parameter from 64% to 16% in $\mathrm{\ensuremath{\Lambda}}\mathrm{CDM}$ model. Through the phenomenological parametrization of GW propagation, the constraints of modified gravity can be improved from 6.2% to 1.6%. Our results show the remarkable significance of eccentricity for the detection and parameter estimation of GW events, allowing us to probe the Universe precisely.