Source author record

Zunli Yuan

Zunli Yuan appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

8works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

8 published item(s)

preprint2022arXiv

A flexible method for estimating luminosity functions via Kernel Density Estimation -- II. Generalization and Python implementation

We propose a generalization of our previous KDE (kernel density estimation) method for estimating luminosity functions (LFs). This new upgrade further extend the application scope of our KDE method, making it a very flexible approach which is suitable to deal with most of bivariate LF calculation problems. From the mathematical point of view, usually the LF calculation can be abstracted as a density estimation problem in the bounded domain of $\{Z_1<z<Z_2,~ L>f_{\mathrm{lim}}(z) \}$. We use the transformation-reflection KDE method ($\hatϕ$) to solve the problem, and introduce an approximate method ($\hatϕ_{\mathrm{1}}$) based on one-dimensional KDE to deal with the small sample size case. In practical applications, the different versions of LF estimators can be flexibly chosen according to the Kolmogorov-Smirnov test criterion. Based on 200 simulated samples, we find that for both cases of dividing or not dividing redshift bins, especially for the latter, our method performs significantly better than the traditional binning method $\hatϕ_{\mathrm{bin}}$. Moreover, with the increase of sample size $n$, our LF estimator converges to the true LF remarkably faster than $\hatϕ_{\mathrm{bin}}$. To implement our method, we have developed a public, open-source Python Toolkit, called \texttt{kdeLF}. With the support of \texttt{kdeLF}, our KDE method is expected to be a competitive alternative to existing nonparametric estimators, due to its high accuracy and excellent stability. \texttt{kdeLF} is available at \url{http://github.com/yuanzunli/kdeLF} with extensive documentation available at \url{http://kdelf.readthedocs.org/en/latest~}.

preprint2022arXiv

Gamma ray burst constraints on cosmological models from the improved Amati correlation

An improved Amati correlation was constructed in (ApJ 931 (2022) 50) by us recently. In this paper, we further study constraints on the $Λ$CDM and $w$CDM models from the gamma ray bursts (GRBs) standardized with the standard and improved Amati correlations, respectively. By using the Pantheon type Ia supernova sample to calibrate the latest A220 GRB data set, the GRB Hubble diagram is obtained model-independently. We find that at the high redshift region ($z>1.4$) the GRB distance modulus from the improved Amati correlation is larger apparently than that from the standard Amati one. The GRB data from the standard Amati correlation only give a lower bound limit on the present matter density parameter $Ω_{\mathrm{m0}}$, while the GRBs from the improved Amati correlation constrain the $Ω_{\mathrm{m0}}$ with the $68\%$ confidence level to be $0.308^{+0.066}_{-0.230}$ and $0.307^{+0.057}_{-0.290}$ in the $Λ$CDM and $w$CDM models, respectively, which are consistent very well with those given by other current popular observational data including BAO, CMB and so on. Once the $H(z)$ data are added in our analysis, the constraint on the Hubble constant $H_0$ can be achieved. We find that two different correlations provide slightly different $H_0$ results but the marginalized mean values seem to be close to that from the Planck 2018 CMB observations.

preprint2022arXiv

The improved Amati correlations from Gaussian copula

In this paper, we obtain two improved Amati correlations of the Gamma-Ray burst (GRB) data via a powerful statistical tool called copula. After calibrating, with the low-redshift GRB data, the improved Amati correlations based on a fiducial $Λ$CDM model with $Ω_\mathrm{m0}=0.3$ and $H_0=70~\mathrm{km~s^{-1}Mpc^{-1}}$, and extrapolating the results to the high-redshift GRB data, we obtain the Hubble diagram of GRB data points. Applying these GRB data to constrain the $Λ$CDM model, we find that the improved Amati correlation from copula can give a result well consistent with $Ω_\mathrm{m0}=0.3$, while the standard Amati and extended Amati correlations do not. This results suggest that when the improved Amati correlation from copula is used in the low-redshift calibration method, the GRB data can be regarded as a viable cosmological explorer. However, the Bayesian information criterion indicates that the standard Amati correlation remains to be favored mildly since it has the least model parameters. Furthermore, once the simultaneous fitting method rather than the low-redshift calibration one is used, there is no apparent evidence that the improved Amati correlation is better than the standard one. Thus, more works need to be done in the future in order to compare different Amati correlations.

preprint2020arXiv

A flexible method for estimating luminosity functions via Kernel Density Estimation

We propose a flexible method for estimating luminosity functions (LFs) based on kernel density estimation (KDE), the most popular nonparametric density estimation approach developed in modern statistics, to overcome issues surrounding binning of LFs. One challenge in applying KDE to LFs is how to treat the boundary bias problem, since astronomical surveys usually obtain truncated samples predominantly due to the flux-density limits of surveys. We use two solutions, the transformation KDE method ($\hatϕ_{\mathrm{t}}$), and the transformation-reflection KDE method ($\hatϕ_{\mathrm{tr}}$) to reduce the boundary bias. We develop a new likelihood cross-validation criterion for selecting optimal bandwidths, based on which, the posterior probability distribution of bandwidth and transformation parameters for $\hatϕ_{\mathrm{t}}$ and $\hatϕ_{\mathrm{tr}}$ are derived within a Markov chain Monte Carlo (MCMC) sampling procedure. The simulation result shows that $\hatϕ_{\mathrm{t}}$ and $\hatϕ_{\mathrm{tr}}$ perform better than the traditional binned method, especially in the sparse data regime around the flux-limit of a survey or at the bright-end of the LF. To further improve the performance of our KDE methods, we develop the transformation-reflection adaptive KDE approach ($\hatϕ_{\mathrm{tra}}$). Monte Carlo simulations suggest that it has a good stability and reliability in performance, and is around an order of magnitude more accurate than using the binned method. By applying our adaptive KDE method to a quasar sample, we find that it achieves estimates comparable to the rigorous determination by a previous work, while making far fewer assumptions about the LF. The KDE method we develop has the advantages of both parametric and non-parametric methods.

preprint2016arXiv

A mixture evolution scenario of the AGN radio luminosity function

We propose a mixture evolution scenario to model the evolution of the radio luminosity function (RLF) of steep spectrum AGNs (active galactic nuclei) based on a Bayesian method. In this scenario, the shape of the RLF is determined by both the density and luminosity evolution. Our models indicate that the density evolution is positive until a redshift of $\thicksim 0.9$, at which point it becomes negative, while the luminosity evolution is positive to a higher redshift ($z \thicksim 5$ for model B and $z \thicksim 3.5$ for model C), where it becomes negative. Our mixture evolution model works well, and the modeled RLFs are in good agreement with previous determinations. The mixture evolution scenario can naturally explain the luminosity-dependent evolution of the RLFs.

preprint2016arXiv

Effect of spectral index distribution on estimating the AGN radio luminosity function

In this paper, we scrutinize the effect of spectral index distribution on estimating the AGN (active galactic nucleus) radio luminosity function (RLF) by a Monte Carlo method. We find that the traditional bivariate RLF estimators can cause bias in varying degree. The bias is especially pronounced for the flat-spectrum radio sources whose spectral index distribution is more scattered. We believe that the bias is caused because the $K$-corrections complicate the truncation boundary on the $L-z$ plane of the sample, but the traditional bivariate RLF estimators have difficulty in dealing with this boundary condition properly. We suggest that the spectral index distribution should be incorporated into the RLF analysis process to obtain a robust estimation. This drives the need for a trivariate function of the form $Φ(α,z,L)$ which we show provides an accurate basis for measuring the RLF.

preprint2013arXiv

A graphical analysis of the systematic error of classical binned methods in constructing luminosity functions

The classical 1/Va and PC methods of constructing binned luminosity functions (LFs) are revisited and compared by graphical analysis. Using both theoretical analysis and illustration with an example, we show why the two methods give different results for the bins which are crossed by the flux limit curves $L=L_{lim}(z)$. Based on a combined sample simulated by a Monte Carlo method, the estimate $ϕ$ of two methods are compared with the input model LFs. The two methods give identical and ideal estimate for the high luminosity points of each redshift interval. However, for the low luminosity bins of all the redshift intervals both methods give smaller estimate than the input model. We conclude that once the LF is evolving with redshift, the classical binned methods will unlikely give an ideal estimate over the total luminosity range. Page & Carrera (2000) noticed that for objects close to the flux limit $ϕ_{1/V_{a}}$ nearly always to be too small. We believe this is due to the arbitrary choosing of redshift and luminosity intervals. Because $ϕ_{1/V_{a}}$ is more sensitive to how the binning are chosen than $ϕ_{PC}$. We suggest a new binning method, which can improve the LFs produced by the 1/Va method significantly, and also improve the LFs produced by the PC methods. Our simulations show that after adopting this new binning, both the 1/Va and PC methods have comparable results.

preprint2011arXiv

On the Evolution of the Cores of Radio Sources and Their Extended Radio Emission

The work in this paper aims at determining the evolution and possible co-evolution of radio-loud active galactic nuclei (AGNs) and their cores via their radio luminosity functions (i.e., total and core RLF respectively). Using a large combined sample of 1063 radio-loud AGNs selected at low radio frequency, we investigate the radio luminosity function (RLF) at 408 MHz of steep-spectrum radio sources. Our results support a luminosity-dependent evolution. Using core flux density data of the complete sample 3CRR, we investigate the core RLF at 5.0 GHz. Based on the combined sample with incomplete core flux data, we also estimate the core RLF using a modified factor of completeness. Both results are consistent and show that the comoving number density of radio cores displays a persistent decline with redshift, implying a negative density evolution. We find that the core RLF is obviously different from the total RLF at 408 MHz band which is mainly contributed by extended lobes, implying that the cores and extended lobes could not be co-evolving at radio emission.