-
PRIYA: A New Suite of Lyman-alpha Forest Simulations for Cosmology
Authors:
Simeon Bird,
Martin Fernandez,
Ming-Feng Ho,
Mahdi Qezlou,
Reza Monadi,
Yueying Ni,
Nianyi Chen,
Rupert Croft,
Tiziana Di Matteo
Abstract:
We present the PRIYA suite of cosmological simulations, based on the code and hydrodynamic model of the ASTRID simulation, and designed for cosmological analyses of the Lyman-$α$ forest. Our simulation suite spans a $9$-dimensional parameter space, including $4$ cosmological parameters and $5$ astrophysical/thermal parameters. We have run $48$ low fidelity simulations with $1536^3$ particles in a…
▽ More
We present the PRIYA suite of cosmological simulations, based on the code and hydrodynamic model of the ASTRID simulation, and designed for cosmological analyses of the Lyman-$α$ forest. Our simulation suite spans a $9$-dimensional parameter space, including $4$ cosmological parameters and $5$ astrophysical/thermal parameters. We have run $48$ low fidelity simulations with $1536^3$ particles in a $120$ Mpc/h box and $3$ high fidelity simulations with $3072^3$ particles in a $120$ Mpc/h box. All our simulations include a full physics model for galaxy formation, including supernova and AGN feedback, and thus also contain a realistic population of DLAs. We advance on earlier simulations suites by larger particle loads, by incorporating new physical models for patchy hydrogen and helium reionization, and by self-consistently incorporating a model for AGN feedback. We show that patchy helium reionization imprints an excess in the 1D flux power spectrum on large scales, which may allow future measurements of helium reionization bubble sizes. Simulation parameters are chosen based on a Latin hypercube design and a Gaussian process is used to interpolate to arbitrary parameter combinations. We build a multi-fidelity emulator for the 1D flux power spectrum and the mean IGM temperature. We show that our final interpolation error is $< 1\%$ and that our simulations produce a flux power spectrum converged at the percent level for $z=5.4$ - $2.2$. Our simulation suite will be used to interpret Lyman-$α$ forest 1D flux power spectra from SDSS and future DESI data releases.
△ Less
Submitted 2 January, 2024; v1 submitted 8 June, 2023;
originally announced June 2023.
-
Machine Learning Uncovers the Universe's Hidden Gems: A Comprehensive Catalogue of CIV Absorption Lines in SDSS DR12
Authors:
Reza Monadi,
Ming-Feng Ho,
Kathy L. Cooksey,
Simeon Bird
Abstract:
We assemble the largest CIV absorption line catalogue to date, leveraging machine learning, specifically Gaussian processes, to remove the need for visual inspection for detecting CIV absorbers. The catalogue contains probabilities classifying the reliability of the absorption system within a quasar spectrum. Our training set was a sub-sample of DR7 spectra that had no detectable CIV absorption in…
▽ More
We assemble the largest CIV absorption line catalogue to date, leveraging machine learning, specifically Gaussian processes, to remove the need for visual inspection for detecting CIV absorbers. The catalogue contains probabilities classifying the reliability of the absorption system within a quasar spectrum. Our training set was a sub-sample of DR7 spectra that had no detectable CIV absorption in a large visually inspected catalogue. We used Bayesian model selection to decide between our continuum model and our absorption-line models. Using a random hold-out sample of 1301 spectra from all of the 26,030 investigated spectra in DR7 CIV catalogue, we validated our pipeline and obtained an 87% classification performance score. We found good purity and completeness values, both ~80%, when a probability of ~95% is used as the threshold. Our pipeline obtained similar CIV redshifts and rest equivalent widths to our training set. Applying our algorithm to 185,425 selected quasar spectra from SDSS DR12, we produce a catalogue of 113,775 CIV doublets with at least 95% confidence. Our catalogue provides maximum a posteriori values and credible intervals for CIV redshift, column density, and Doppler velocity dispersion. We detect CIV absorption systems with a redshift range of 1.37 $\!-\!$ 5.1, including 33 systems with a redshift larger than 5 and 549 absorbers systems with a rest equivalent width greater than 2 A at more than 95% confidence. Our catalogue can be used to investigate the physical properties of the circumgalactic and intergalactic media.
△ Less
Submitted 23 September, 2023; v1 submitted 28 April, 2023;
originally announced May 2023.
-
Improved selection of extremely red quasars with boxy CIV lines in BOSS
Authors:
Reza Monadi,
Simeon Bird
Abstract:
Extremely red quasars (ERQs) are an interesting sample of quasars in the Baryon Oscillation Spectroscopic Sample (BOSS) in the redshift range of $2.0 - 3.4$ and have extreme red colours of $i-W3\ge4.6$. Core ERQs have strong CIV emission lines with rest equivalent width of $\ge100$Å. Many core ERQs also have CIV line profiles with peculiar boxy shapes which distinguish them from normal blue quasar…
▽ More
Extremely red quasars (ERQs) are an interesting sample of quasars in the Baryon Oscillation Spectroscopic Sample (BOSS) in the redshift range of $2.0 - 3.4$ and have extreme red colours of $i-W3\ge4.6$. Core ERQs have strong CIV emission lines with rest equivalent width of $\ge100$Å. Many core ERQs also have CIV line profiles with peculiar boxy shapes which distinguish them from normal blue quasars. We show, using a combination of kernel density estimation and local outlier factor analyses on a space of the $i-W3$ colour, CIV rest equivalent width and line kurtosis, that core ERQs likely represent a
separate population rather than a smooth transition between normal
blue quasars and the quasars in the tail of the colour-REW distribution. We apply our analyses to find new criteria for selecting ERQs in this 3D parameter space. Our final selection produces $133$ quasars, which are \emph{three} times more likely to have a visually verified CIV broad absorption line feature than the previous core ERQ sample. We further show that our newly selected sample are extreme objects in the intersection of the WISE AGN catalogue with the MILLIQUAS quasar catalogue in the colour-colour space of ($W1-W2$, $W2-W3$). This paper validates an improved selection method for red quasars which can be applied to future datasets such as the quasar catalogue from the Dark Energy Spectroscopic Instrument (DESI).
△ Less
Submitted 29 January, 2022; v1 submitted 25 July, 2021;
originally announced July 2021.
-
Statistical Analysis of I Stokes Parameter of Millisecond Pulsars
Authors:
Hossein Panahi,
Issa Eghdami,
Reza Monadi
Abstract:
Using Detrended Fluctuation Analysis (DFA) and box counting method, we test spacial correlation and fractality of Polarization Pulse Profiles (PPPs) of 24 millisecond pulsars (MSPs) which were observed in Parkes Pulsar Timing Array (PPTA) project. DFA analysis indicates that MSPs' PPPs are persistent and the results of box counting method confirm the fractality in the majority of PPPs. A Kolmogoro…
▽ More
Using Detrended Fluctuation Analysis (DFA) and box counting method, we test spacial correlation and fractality of Polarization Pulse Profiles (PPPs) of 24 millisecond pulsars (MSPs) which were observed in Parkes Pulsar Timing Array (PPTA) project. DFA analysis indicates that MSPs' PPPs are persistent and the results of box counting method confirm the fractality in the majority of PPPs. A Kolmogorov-Smirnov test indicates that isolated MSPs have more complex PPPs than binary ones. Then we apply our analysis on a random sample of normal pulsars. Comparing the results of our analysis on MSPs and normal pulsars shows that MSPs have more complex PPPs which is resulted from smaller angular half-width of the emission cone and more peaks in MSPs PPPs. On the other hand, high values of Hurst exponent in MSPs confirm compact emission regions in these pulsars.
△ Less
Submitted 28 March, 2016;
originally announced April 2016.
-
Gaussian Anisotropy In Strange Quark Stars
Authors:
H. Panahi,
R. Monadi,
I. Eghdami
Abstract:
In this paper for studying the anisotropic strange quark stars, we assume that the radial pressure inside the anisotropic star is a superposition of pressure in an isotropic case plus a Gaussian perturbation term. Considering a proportionality between electric charge density and the density of matter, we solve the TOV equation for different cases numerically. Our results indicate that anisotropy i…
▽ More
In this paper for studying the anisotropic strange quark stars, we assume that the radial pressure inside the anisotropic star is a superposition of pressure in an isotropic case plus a Gaussian perturbation term. Considering a proportionality between electric charge density and the density of matter, we solve the TOV equation for different cases numerically. Our results indicate that anisotropy increases the maximum mass $M_{max}$ and also its corresponding radius $R$ for a typical strange quark star. According to our calculations, an anisotropy amplitude of $A=3\times10^{33}Nm^{-2}$ with a standard deviation of $σ=3\times10^{3}m$ leads to a neutron star of 1.97$M_{\odot}$. Furthermore, electric charge not only increases the maximum mass and its corresponding radius, but also raises up the anisotropy factor. We can see that the tangential pressure $p_{t}$ and anisotropy factor $Δ$ unlike the radial pressure $p_{r}$ have a maximum on the surface and this maximum increases by adding electric charge effect. However, we show that anisotropy can be more effective than electric charge in rasing maximum mass of strange quark stars.
△ Less
Submitted 26 April, 2015;
originally announced April 2015.