Showing 1–50 of 60 results for author: Hong, E

The Environmental Dependence of the Stellar Mass - Gas Metallicity Relation in Horizon Run 5

Authors: Aaron R. Rowntree, Ankit Singh, Fiorenzo Vincenzo, Brad K. Gibson, Céline Gouin, Daniela Galárraga-Espinosa, Jaehyun Lee, Juhan Kim, Clotilde Laigle, Changbom Park, Christophe Pichon, Gareth Few, Sungwook E. Hong, Yonghwi Kim

Abstract: Metallicity offers a unique window into the baryonic history of the cosmos, being instrumental in probing evolutionary processes in galaxies between different cosmic environments. We aim to quantify the contribution of these environments to the scatter in the mass-metallicity relation (MZR) of galaxies. By analysing the galaxy distribution within the cosmic skeleton of the Horizon Run 5 cosmologic… ▽ More Metallicity offers a unique window into the baryonic history of the cosmos, being instrumental in probing evolutionary processes in galaxies between different cosmic environments. We aim to quantify the contribution of these environments to the scatter in the mass-metallicity relation (MZR) of galaxies. By analysing the galaxy distribution within the cosmic skeleton of the Horizon Run 5 cosmological hydrodynamical simulation at redshift $z = 0.625$, computed using a careful calibration of the T-ReX filament finder, we identify galaxies within three main environments: nodes, filaments and voids. We also classify galaxies based on the dynamical state of the clusters and the length of the filaments in which they reside. We find that the cosmic environment significantly contributes to the scatter in the MZR; in particular, both the gas metallicity and its average relative standard deviation increase when considering denser large-scale environments. The difference in the average metallicity between galaxies within relaxed and unrelaxed clusters is $\approx 0.1 \text{ dex}$, with both populations displaying positive residuals, $δZ_{g}$, from the averaged MZR. Moreover, the difference in metallicity between node and void galaxies accounts for $\approx 0.14 \, \text{dex}$ in the scatter of the MZR at stellar mass $M_{\star} \approx 10^{9.35}\,\text{M}_{\odot}$. Finally, both the average [O/Fe] in the gas and the galaxy gas fraction decrease when moving to higher large-scale densities in the simulation, suggesting that the cores of cosmic environments host, on average, older and more massive galaxies, whose enrichment is affected by a larger number of Type Ia Supernova events. △ Less

Submitted 3 June, 2024; v1 submitted 15 April, 2024; originally announced April 2024.

Comments: 17 pages, 17 figures

Maunakea Spectroscopic Explorer exposure time calculator for end-to-end simulator: to optimizing spectrograph design and observing simulation

Authors: Tae-Geun Ji, Jennifer Sobeck, Changgon Kim, Hojae Ahn, Mingyeong Yang, Taeeun Kim, Sungwook E. Hong, Kei Szeto, Jennifer L. Marshall, Christian Surace, Soojong Pak

Abstract: The Maunakea Spectroscopic Explorer (MSE) project will provide multi-object spectroscopy in the optical and near-infrared bands using an 11.25-m aperture telescope, repurposing the original Canada-France-Hawaii Telescope (CFHT) site. MSE will observe 4,332 objects per single exposure with a field of view of 1.5 square degrees, utilizing two spectrographs with low-moderate (R$\sim$3,000, 6,000) and… ▽ More The Maunakea Spectroscopic Explorer (MSE) project will provide multi-object spectroscopy in the optical and near-infrared bands using an 11.25-m aperture telescope, repurposing the original Canada-France-Hawaii Telescope (CFHT) site. MSE will observe 4,332 objects per single exposure with a field of view of 1.5 square degrees, utilizing two spectrographs with low-moderate (R$\sim$3,000, 6,000) and high (R$\approx$30,000) spectral resolution. In general, an exposure time calculator (ETC) is used to estimate the performance of an observing system by calculating a signal-to-noise ratio (S/N) and exposure time. We present the design of the MSE exposure time calculator (ETC), which has four calculation modes (S/N, exposure time, S/N trend with wavelength, and S/N trend with magnitude) and incorporates the MSE system requirements as specified in the Conceptual Design. The MSE ETC currently allows for user-defined inputs of target AB magnitude, water vapor, airmass, and sky brightness AB magnitude (additional user inputs can be provided depending on computational mode). The ETC is built using Python 3.7 and features a graphical user interface that allows for cross-platform use. The development process of the ETC software follows an Agile methodology and utilizes the Unified Modeling Language (UML) diagrams to visualize the software architecture. We also describe the testing and verification of the MSE ETC. △ Less

Submitted 16 January, 2024; originally announced January 2024.

Comments: 30 pages, 13 figures

Journal ref: J. Astron. Telesc. Instrum. Syst., Vol. 10, Issue 01, 018001 (2024)

Merger Tree-based Galaxy Matching: A Comparative Study Across Different Resolutions

Authors: Minyong Jung, Ji-hoon Kim, Boon Kiat Oh, Sungwook E. Hong, Jaehyun Lee, Juhan Kim

Abstract: We introduce a novel halo/galaxy matching technique between two cosmological simulations with different resolutions, which utilizes the positions and masses of halos along their subhalo merger tree. With this tool, we conduct a study of resolution biases through the {\it galaxy-by-galaxy} inspection of a pair of simulations that have the same simulation configuration but different mass resolutions… ▽ More We introduce a novel halo/galaxy matching technique between two cosmological simulations with different resolutions, which utilizes the positions and masses of halos along their subhalo merger tree. With this tool, we conduct a study of resolution biases through the {\it galaxy-by-galaxy} inspection of a pair of simulations that have the same simulation configuration but different mass resolutions, utilizing a suite of {\sc IllustrisTNG} simulations to assess the impact on galaxy properties. We find that, with the subgrid physics model calibrated for TNG100-1, subhalos in TNG100-1 (high resolution) have $\lesssim0.5$ dex higher stellar masses than their counterparts in the TNG100-2 (low-resolution). It is also discovered that the subhalos with $M_{\mathrm{gas}}\sim10^{8.5}\,{\rm M}_\odot$ in TNG100-1 have $\sim0.5$ dex higher gas mass than those in TNG100-2. The mass profiles of the subhalos reveal that the dark matter masses of subhalos in TNG100-2 converge well with those from TNG100-1, except within 4 kpc of the resolution limit. The differences in stellar mass and hot gas mass are most pronounced in the central region. We exploit machine learning to build a correction mapping for the physical quantities of subhalos from low- to high-resolution simulations (TNG300-1 and TNG100-1), which enables us to find an efficient way to compile a high-resolution galaxy catalog even from a low-resolution simulation. Our tools can easily be applied to other large cosmological simulations, testing and mitigating the resolution biases of their numerical codes and subgrid physics models. △ Less

Submitted 19 March, 2024; v1 submitted 4 December, 2023; originally announced December 2023.

Comments: 24 pages, 14 figures; Accepted for publication in to ApJ

Journal ref: ApJ 965 156 (2024)

Final parsec problem of black hole mergers and ultralight dark matter

Authors: Hyeonmo Koo, Dongsu Bak, Inkyu Park, Sungwook E. Hong, Jae-Weon Lee

Abstract: When two galaxies merge, they often produce a supermassive black hole binary (SMBHB) at their center. Numerical simulations with cold dark matter show that SMBHBs typically stall out at a distance of a few parsecs apart, and take billions of years to coalesce. This is known as the final parsec problem. We suggest that ultralight dark matter (ULDM) halos around SMBHBs can generate dark matter waves… ▽ More When two galaxies merge, they often produce a supermassive black hole binary (SMBHB) at their center. Numerical simulations with cold dark matter show that SMBHBs typically stall out at a distance of a few parsecs apart, and take billions of years to coalesce. This is known as the final parsec problem. We suggest that ultralight dark matter (ULDM) halos around SMBHBs can generate dark matter waves due to gravitational cooling. These waves can effectively carry away orbital energy from the black holes, rapidly driving them together. To test this hypothesis, we performed numerical simulations of black hole binaries inside ULDM halos. Our results imply that ULDM waves can lead to the rapid orbital decay of black hole binaries. △ Less

Submitted 6 November, 2023; originally announced November 2023.

The effects of non-linearity on the growth rate constraint from velocity correlation functions

Authors: Motonari Tonegawa, Stephen Appleby, Changbom Park, Sungwook E. Hong, Juhan Kim

Abstract: The two-point statistics of the cosmic velocity field, measured from galaxy peculiar velocity (PV) surveys, can be used as a dynamical probe to constrain the growth rate of large-scale structures in the universe. Most works use the statistics on scales down to a few tens of Megaparsecs, while using a theoretical template based on the linear theory. In addition, while the cosmic velocity is volume-… ▽ More The two-point statistics of the cosmic velocity field, measured from galaxy peculiar velocity (PV) surveys, can be used as a dynamical probe to constrain the growth rate of large-scale structures in the universe. Most works use the statistics on scales down to a few tens of Megaparsecs, while using a theoretical template based on the linear theory. In addition, while the cosmic velocity is volume-weighted, the observable line-of-sight velocity two-point correlation is density-weighted, as sampled by galaxies, and therefore the density-velocity correlation term also contributes, which has often been neglected. These effects are fourth order in powers of the linear density fluctuation $δ_{\rm L}^4$, compared to $δ_{\rm L}^2$ of the linear velocity correlation function, and have the opposite sign. We present these terms up to $δ_{\rm L}^4$ in real space based on the standard perturbation theory, and investigate the effect of non-linearity and the density-velocity contribution on the inferred growth rate $fσ_8$, using $N$-body simulations. We find that for a next-generation PV survey of volume $\sim {\cal O}(500 \, h^{-1} \, {\rm Mpc})^3$, these effects amount to a shift of $fσ_8$ by $\sim 10$ per cent and is comparable to the forecasted statistical error when the minimum scale used for parameter estimation is $r_{\rm min} = 20 \, h^{-1} \, {\rm Mpc}$. △ Less

Submitted 6 March, 2024; v1 submitted 25 September, 2023; originally announced September 2023.

Comments: 17 pages, 6 figures. Minor corrections to non-linear calculation in section 2.2, conclusions unchanged. Accepted for publication in MNRAS

Cluster-counterpart Voids: Void Identification from Galaxy Density Field

Authors: Junsup Shim, Changbom Park, Juhan Kim, Sungwook E. Hong

Abstract: We identify cosmic voids from galaxy density fields under the theory of void-cluster correspondence. We extend the previous novel void-identification method developed for the matter density field to the galaxy density field for practical applications. From cosmological N-body simulations, we construct galaxy number- and mass-weighted density fields to identify cosmic voids that are counterparts of… ▽ More We identify cosmic voids from galaxy density fields under the theory of void-cluster correspondence. We extend the previous novel void-identification method developed for the matter density field to the galaxy density field for practical applications. From cosmological N-body simulations, we construct galaxy number- and mass-weighted density fields to identify cosmic voids that are counterparts of galaxy clusters of specific mass. The parameters for the cluster-counterpart void identification such as Gaussian smoothing scale, density threshold, and core volume fraction are found for galaxy density fields. We achieve about $60$--$67\%$ of completeness and reliability for identifying the voids of corresponding cluster mass above $3\times10^{14}h^{-1}M_{\odot}$ from a galaxy sample with the mean number density, $\bar{n}=4.4\times10^{-3} (h^{-1}{\rm Mpc})^{-3}$. When the mean density is increased to $\bar{n}=10^{-2} (h^{-1}{\rm Mpc})^{-3}$, the detection rate is enhanced by $\sim2$--$7\%$ depending on the `mass scale' of voids. We find that the detectability is insensitive to the density weighting scheme applied to generate the density field. Our result demonstrates that we can apply this method to the galaxy redshift survey data to identify cosmic voids corresponding statistically to the galaxy clusters in a given mass range. △ Less

Submitted 16 May, 2023; originally announced May 2023.

Comments: 10 pages, 5 figures, 1 table. Accepted for publication in ApJ

Cosmological Parameter Constraints from the SDSS Density and Momentum Power Spectra

Authors: Stephen Appleby, Motonari Tonegawa, Changbom Park, Sungwook E. Hong, Juhan Kim, Yongmin Yoon

Abstract: We extract the galaxy density and momentum power spectra from a subset of early-type galaxies in the SDSS DR7 main galaxy catalog. Using galaxy distance information inferred from the improved fundamental plane described in \citet{Yoon_2020}, we reconstruct the peculiar velocities of the galaxies and generate number density and density-weighted velocity fields, from which we extract the galaxy dens… ▽ More We extract the galaxy density and momentum power spectra from a subset of early-type galaxies in the SDSS DR7 main galaxy catalog. Using galaxy distance information inferred from the improved fundamental plane described in \citet{Yoon_2020}, we reconstruct the peculiar velocities of the galaxies and generate number density and density-weighted velocity fields, from which we extract the galaxy density and momentum power spectra. We compare the measured values to the theoretical expectation of the same statistics, assuming an input $Λ$CDM model and using a third-order perturbative expansion. After validating our analysis pipeline with a series of mock data sets, we apply our methodology to the SDSS data and arrive at constraints $fσ_{8} = 0.471_{-0.080}^{+0.077}$ and $b_{1}σ_{8} = 0.920_{-0.070}^{+0.070}$ at a mean redshift $\bar{z} = 0.04$. Our result is consistent with the Planck cosmological best fit parameters for the $Λ$CDM model. The momentum power spectrum is found to be strongly contaminated by small scale velocity dispersion, which suppresses power by $\sim {\cal O}(30\%)$ on intermediate scales $k \sim 0.05 \, h \, {\rm Mpc}^{-1}$. △ Less

Submitted 2 October, 2023; v1 submitted 3 May, 2023; originally announced May 2023.

Comments: 15 figures, 4 tables, accepted for ApJ

Tomographic Alcock-Paczynski Test with Redshift-Space Correlation Function: Evidence for the Dark Energy Equation of State Parameter w>-1

Authors: Fuyu Dong, Changbom Park, Sungwook E. Hong, Juhan Kim, Ho Seong Hwang, Hyunbae Park, Stephen Appleby

Abstract: The apparent shape of galaxy clustering depends on the adopted cosmology used to convert observed redshift to comoving distance, the $r(z)$ relation, as it changes the line elements along and across the line of sight differently. The Alcock-Paczyński (AP) test exploits this property to constrain the expansion history of the universe. We present an extensive review of past studies on the AP test. W… ▽ More The apparent shape of galaxy clustering depends on the adopted cosmology used to convert observed redshift to comoving distance, the $r(z)$ relation, as it changes the line elements along and across the line of sight differently. The Alcock-Paczyński (AP) test exploits this property to constrain the expansion history of the universe. We present an extensive review of past studies on the AP test. We adopt an extended AP test method introduced by Park et al. (2019), which uses the full shape of redshift-space two-point correlation function (CF) as the standard shape, and apply it to the SDSS DR7, BOSS, and eBOSS LRG samples covering the redshift range up to $z=0.8$.We calibrate the test against the nonlinear cosmology-dependent systematic evolution of the CF shape using the Multiverse simulations. We focus on examining whether or not the flat $Λ$CDM `concordance' model is consistent with observation. We constrain the flat $w$CDM model to have $w=-0.892_{-0.050}^{+0.045}$ and $Ω_m=0.282_{-0.023}^{+0.024}$ from our AP test alone, which is significantly tighter than the constraints from the BAO or SNe I$a$ methods by a factor of 3 - 6. When the AP test result is combined with the recent BAO and SNe I$a$ results, we obtain $w=-0.903_{-0.023}^{+0.023}$ and $Ω_m=0.285_{-0.009}^{+0.014}$. This puts a strong tension with the flat $Λ$CDM model with $w=-1$ at $4.2σ$ level. Consistency with $w=-1$ is obtained only when the Planck CMB observation is combined. It remains to see if this tension between observations of galaxy distribution at low redshifts and CMB anisotropy at the decoupling epoch becomes greater in the future studies and leads us to a new paradigm of cosmology. △ Less

Submitted 29 April, 2023; originally announced May 2023.

Comments: 21 pages, 11 figures, accepted by ApJ

The Universe is worth $64^3$ pixels: Convolution Neural Network and Vision Transformers for Cosmology

Authors: Se Yeon Hwang, Cristiano G. Sabiu, Inkyu Park, Sungwook E. Hong

Abstract: We present a novel approach for estimating cosmological parameters, $Ω_m$, $σ_8$, $w_0$, and one derived parameter, $S_8$, from 3D lightcone data of dark matter halos in redshift space covering a sky area of $40^\circ \times 40^\circ$ and redshift range of $0.3 < z < 0.8$, binned to $64^3$ voxels. Using two deep learning algorithms, Convolutional Neural Network (CNN) and Vision Transformer (ViT),… ▽ More We present a novel approach for estimating cosmological parameters, $Ω_m$, $σ_8$, $w_0$, and one derived parameter, $S_8$, from 3D lightcone data of dark matter halos in redshift space covering a sky area of $40^\circ \times 40^\circ$ and redshift range of $0.3 < z < 0.8$, binned to $64^3$ voxels. Using two deep learning algorithms, Convolutional Neural Network (CNN) and Vision Transformer (ViT), we compare their performance with the standard two-point correlation (2pcf) function. Our results indicate that CNN yields the best performance, while ViT also demonstrates significant potential in predicting cosmological parameters. By combining the outcomes of Vision Transformer, Convolution Neural Network, and 2pcf, we achieved a substantial reduction in error compared to the 2pcf alone. To better understand the inner workings of the machine learning algorithms, we employed the Grad-CAM method to investigate the sources of essential information in activation maps of the CNN and ViT. Our findings suggest that the algorithms focus on different parts of the density field and redshift depending on which parameter they are predicting. This proof-of-concept work paves the way for incorporating deep learning methods to estimate cosmological parameters from large-scale structures, potentially leading to tighter constraints and improved understanding of the Universe. △ Less

Submitted 2 November, 2023; v1 submitted 17 April, 2023; originally announced April 2023.

Comments: 23 pages, 10 figures

Cosmic gas highways in C-EAGLE simulations

Authors: I. Vurm, J. Nevalainen, S. E. Hong, Y. M. Bahé, C. Dalla Vecchia, P. Heinämäki

Abstract: A substantial fraction of the cosmic baryons is expected to hide in the form of diffuse warm-hot intergalactic medium (WHIM), the majority of which resides in the filaments of the Cosmic Web and has proven very difficult to detect due to its low density. Close to galaxy clusters, the filament gas is affected by the cluster's gravitational potential and attains substantial infall velocities, eventu… ▽ More A substantial fraction of the cosmic baryons is expected to hide in the form of diffuse warm-hot intergalactic medium (WHIM), the majority of which resides in the filaments of the Cosmic Web and has proven very difficult to detect due to its low density. Close to galaxy clusters, the filament gas is affected by the cluster's gravitational potential and attains substantial infall velocities, eventually undergoing a termination shock which may boost its X-ray signal. We aim to identify the optimal locations of the enhanced X-ray emission and absorption arising from cluster-filament interactions, as well as improve our understanding of the various physical processes affecting the WHIM as it approaches the cluster. We applied the DisPerSE filament finder to the galaxy distribution in the surroundings of a Coma-like ($M_{200} \sim 10^{15.4}~{\rm M}_{\odot}$) simulated C-EAGLE galaxy cluster. We characterised the thermodynamic properties of the gas in such filaments as well as their dependence on distance from the cluster, and provided a physical interpretation for the results. The identified filaments account for $\sim 50$% of the hot WHIM ($T > 10^{5.5}$ K) in the cluster vicinity. The filament gas is in approximate free-fall all the way down to $\sim 2 \ r_{200}$ from the cluster, at which stage it begins to slow down due to the increasing pressure of the ambient gas. The deceleration is accompanied by the conversion of gas bulk kinetic energy into heat and the increase of density and temperature from the general Cosmic Web level of $ρ\sim 10ρ_{\rm av}$ and $T = 10^5-10^6$ K towards $ρ\sim 100ρ_{\rm av}$ and $T = 10^7-10^8$ K near the cluster boundary. We conclude that the detection of the cosmic filaments of galaxies around clusters may provide a practical observational avenue for locating the densest and hottest phase of the missing baryons. △ Less

Submitted 6 March, 2023; originally announced March 2023.

Comments: 20 pages, 16 figures, accepted for publication in A&A

Reconstructing the cosmological density and velocity fields from redshifted galaxy distributions using V-net

Authors: Fei Qin, David Parkinson, Sungwook E. Hong, Cristiano G. Sabiu

Abstract: The distribution of matter that is measured through galaxy redshift and peculiar velocity surveys can be harnessed to learn about the physics of dark matter, dark energy, and the nature of gravity. To improve our understanding of the matter of the Universe, we can reconstruct the full density and velocity fields from the galaxies that act as tracer particles. In this paper, we use the simulated ha… ▽ More The distribution of matter that is measured through galaxy redshift and peculiar velocity surveys can be harnessed to learn about the physics of dark matter, dark energy, and the nature of gravity. To improve our understanding of the matter of the Universe, we can reconstruct the full density and velocity fields from the galaxies that act as tracer particles. In this paper, we use the simulated halos as proxies for the galaxies. We use a convolutional neural network, a V-net, trained on numerical simulations of structure formation to reconstruct the density and velocity fields. We find that, with detailed tuning of the loss function, the V-net could produce better fits to the density field in the high-density and low-density regions, and improved predictions for the probability distribution of the amplitudes of the velocities. However, the weights will reduce the precision of the estimated $β$ parameter. We also find that the redshift-space distortions of the halo catalogue do not significantly contaminate the reconstructed real-space density and velocity field. We estimate the velocity field $β$ parameter by comparing the peculiar velocities of halo catalogues to the reconstructed velocity fields, and find the estimated $β$ values agree with the fiducial value at the 68\% confidence level. △ Less

Submitted 19 June, 2023; v1 submitted 3 February, 2023; originally announced February 2023.

Comments: 25 pages, 18 figures, 2 tables. Published JCAP

Journal ref: JCAP06(2023)062

MulGuisin, a Topological Network Finder and its Performance on Galaxy Clustering

Authors: Young Ju, Inkyu Park, Cristiano G. Sabiu, Sungwook E. Hong

Abstract: We introduce a new clustering algorithm, MulGuisin (MGS), that can identify distinct galaxy over-densities using topological information from the galaxy distribution. This algorithm was first introduced in an LHC experiment as a Jet Finder software, which looks for particles that clump together in close proximity. The algorithm preferentially considers particles with high energies and merges them… ▽ More We introduce a new clustering algorithm, MulGuisin (MGS), that can identify distinct galaxy over-densities using topological information from the galaxy distribution. This algorithm was first introduced in an LHC experiment as a Jet Finder software, which looks for particles that clump together in close proximity. The algorithm preferentially considers particles with high energies and merges them only when they are closer than a certain distance to create a jet. MGS shares some similarities with the minimum spanning tree (MST) since it provides both clustering and network-based topology information. Also, similar to the density-based spatial clustering of applications with noise (DBSCAN), MGS uses the ranking or the local density of each particle to construct clustering. In this paper, we compare the performances of clustering algorithms using controlled data and some realistic simulation data as well as the SDSS observation data, and we demonstrate that our new algorithm finds networks most efficiently and defines galaxy networks in a way that most closely resembles human vision. △ Less

Submitted 18 February, 2024; v1 submitted 9 January, 2023; originally announced January 2023.

Comments: 26 pages, 13 figures

Panel Discussion: Practical Problem Solving for Machine Learning

Authors: Guillermo Cabrera, Sungwook E. Hong, Lilianne Nakazono, David Parkinson, Yuan-Sen Ting

Abstract: Machine Learning is a powerful tool for astrophysicists, which has already had significant uptake in the community. But there remain some barriers to entry, relating to proper understanding, the difficulty of interpretability, and the lack of cohesive training. In this discussion session we addressed some of these questions, and suggest how the field may move forward. Machine Learning is a powerful tool for astrophysicists, which has already had significant uptake in the community. But there remain some barriers to entry, relating to proper understanding, the difficulty of interpretability, and the lack of cohesive training. In this discussion session we addressed some of these questions, and suggest how the field may move forward. △ Less

Submitted 30 November, 2022; originally announced November 2022.

Comments: 6 pages. Prepared for the proceedings of the International Astronomical Union Symposium 368 "Machine Learning in Astronomy: Possibilities and Pitfalls"

Modeling Cosmological Perturbations of Thermal Inflation

Authors: Jeong-Myeong Bae, Sungwook E. Hong, Heeseung Zoe

Abstract: We consider a simple system consisting of matter, radiation and vacuum components to model the impact of thermal inflation on the evolution of primordial perturbations. The vacuum energy magnifies the primordial modes entering the horizon before its domination, making them potentially observable, and the resulting transfer function reflects the phase changes and energy contents. To determine the t… ▽ More We consider a simple system consisting of matter, radiation and vacuum components to model the impact of thermal inflation on the evolution of primordial perturbations. The vacuum energy magnifies the primordial modes entering the horizon before its domination, making them potentially observable, and the resulting transfer function reflects the phase changes and energy contents. To determine the transfer function, we follow the curvature perturbation from well outside the horizon during radiation domination to well outside the horizon during vacuum domination and evaluate it on a constant radiation density hypersurface, as is appropriate for the case of thermal inflation. The shape of the transfer function is determined by the ratio of vacuum energy to radiation at matter-radiation equality, which we denote by $\upsilon$, and has two characteristic scales, $k_{\rm a}$ and $k_{\rm b}$, corresponding to the horizon sizes at matter radiation equality and the beginning of the inflation, respectively. If $\upsilon \ll 1$, the universe experiences radiation, matter and vacuum domination eras and the transfer function is flat for $k \ll k_{\rm b}$, oscillates with amplitude $1/5$ for $ k_{\rm b} \ll k \ll k_{\rm a}$ and oscillates with amplitude $1$ for $k \gg k_{\rm a}$. For $\upsilon \gg 1$, the matter domination era disappears, and the transfer function reduces to being flat for $k \ll k_{\rm b}$ and oscillating with amplitude $1$ for $k \gg k_{\rm b}$. △ Less

Submitted 8 December, 2023; v1 submitted 19 April, 2022; originally announced April 2022.

Comments: 22 pages, 6 figures, accepted to the Class. Quant. Grav

Machine-guided Exploration and Calibration of Astrophysical Simulations

Authors: Boon Kiat Oh, Hongjun An, Eun-jin Shin, Ji-hoon Kim, Sungwook E. Hong

Abstract: We apply a novel method with machine learning to calibrate sub-grid models within numerical simulation codes to achieve convergence with observations and between different codes. It utilizes active learning and neural density estimators. The hyper parameters of the machine are calibrated with a well-defined projectile motion problem. Then, using a set of 22 cosmological zoom simulations, we tune t… ▽ More We apply a novel method with machine learning to calibrate sub-grid models within numerical simulation codes to achieve convergence with observations and between different codes. It utilizes active learning and neural density estimators. The hyper parameters of the machine are calibrated with a well-defined projectile motion problem. Then, using a set of 22 cosmological zoom simulations, we tune the parameters of a popular star formation and feedback model within Enzo to match simulations. The parameters that are adjusted include the star formation efficiency, coupling of thermal energy from stellar feedback, and volume into which the energy is deposited. This number translates to a factor of more than three improvements over manual calibration. Despite using fewer simulations, we obtain a better agreement to the observed baryon makeup of a Milky-Way (MW) sized halo. Switching to a different strategy, we improve the consistency of the recommended parameters from the machine. Given the success of the calibration, we then apply the technique to reconcile metal transport between grid-based and particle-based simulation codes using an isolated galaxy. It is an improvement over manual exploration while hinting at a less known relation between the diffusion coefficient and the metal mass in the halo region. The exploration and calibration of the parameters of the sub-grid models with a machine learning approach is concluded to be versatile and directly applicable to different problems. △ Less

Submitted 14 March, 2022; originally announced March 2022.

A low-threshold ultrahigh-energy neutrino search with the Askaryan Radio Array

Authors: P. Allison, S. Archambault, J. J. Beatty, D. Z. Besson, A. Bishop, C. C. Chen, C. H. Chen, P. Chen, Y. C. Chen, B. A. Clark, W. Clay, A. Connolly, L. Cremonesi, P. Dasgupta, J. Davies, S. de Kockere, K. D. de Vries, C. Deaconu, M. A. DuVernois, J. Flaherty, E. Friedman, R. Gaior, J. Hanson, N. Harty, B. Hendricks , et al. (55 additional authors not shown)

Abstract: In the pursuit of the measurement of the still-elusive ultrahigh-energy (UHE) neutrino flux at energies of order EeV, detectors using the in-ice Askaryan radio technique have increasingly targeted lower trigger thresholds. This has led to improved trigger-level sensitivity to UHE neutrinos. Working with data collected by the Askaryan Radio Array (ARA), we search for neutrino candidates at the lowe… ▽ More In the pursuit of the measurement of the still-elusive ultrahigh-energy (UHE) neutrino flux at energies of order EeV, detectors using the in-ice Askaryan radio technique have increasingly targeted lower trigger thresholds. This has led to improved trigger-level sensitivity to UHE neutrinos. Working with data collected by the Askaryan Radio Array (ARA), we search for neutrino candidates at the lowest threshold achieved to date, leading to improved analysis-level sensitivities. A neutrino search on a data set with 208.7~days of livetime from the reduced-threshold fifth ARA station is performed, achieving a 68\% analysis efficiency over all energies on a simulated mixed-composition neutrino flux with an expected background of $0.10_{-0.04}^{+0.06}$ events passing the analysis. We observe one event passing our analysis and proceed to set a neutrino flux limit using a Feldman-Cousins construction. We show that the improved trigger-level sensitivity can be carried through an analysis, motivating the Phased Array triggering technique for use in future radio-detection experiments. We also include a projection using all available data from this detector. Finally, we find that future analyses will benefit from studies of events near the surface to fully understand the background expected for a large-scale detector. △ Less

Submitted 14 February, 2022; originally announced February 2022.

Comments: 15 pages, 8 figures

Minkowski Functionals of SDSS-III BOSS : Hints of Possible Anisotropy in the Density Field?

Authors: Stephen Appleby, Changbom Park, Pratyush Pranav, Sungwook E. Hong, Ho Seong Hwang, Juhan Kim, Thomas Buchert

Abstract: We present measurements of the Minkowski functionals extracted from the SDSS-III BOSS catalogs. After defining the Minkowski functionals, we describe how an unbiased reconstruction of these statistics can be obtained from a field with masked regions and survey boundaries, validating our methodology with Gaussian random fields and mock galaxy snapshot data. From the BOSS galaxy data we generate a s… ▽ More We present measurements of the Minkowski functionals extracted from the SDSS-III BOSS catalogs. After defining the Minkowski functionals, we describe how an unbiased reconstruction of these statistics can be obtained from a field with masked regions and survey boundaries, validating our methodology with Gaussian random fields and mock galaxy snapshot data. From the BOSS galaxy data we generate a set of four density fields in three dimensions corresponding to the northern and southern skies of LOWZ and CMASS catalogs, smoothing over large scales such that the field is perturbatively non-Gaussian. We extract the Minkowski functionals from each data set separately, and measure their shapes and amplitudes by fitting a Hermite polynomial expansion. For the shape parameter of the Minkowski functional curves $a_0$, that is related to the bispectrum of the field, we find that the LOWZ-South data presents a systematically lower value of $a_0 = -0.080 \pm 0.040$ than its northern sky counterpart $a_0 = 0.032 \pm 0.024$. Although the significance of this discrepancy is low, it potentially indicates some systematics in the data or that the matter density field exhibits anisotropy at low redshift. By assuming a standard isotropic flat $Λ$CDM cosmology, the amplitudes of Minkowski functionals from the combination of northern and southern sky data give the constraints $Ω_{\rm c} h^2 n_{\rm s} = 0.110 \pm 0.006$ and $0.111 \pm 0.008$ for CMASS and LOWZ, respectively, which is in agreement with the Planck $Λ$CDM best-fit $Ω_{\rm c}h^{2} n_{\rm s} = 0.116 \pm 0.001$. △ Less

Submitted 12 October, 2021; originally announced October 2021.

Panspermia in a Milky Way-like Galaxy

Authors: Raphael Gobat, Sungwook E. Hong, Owain Snaith, Sungryong Hong

Abstract: We study the process of panspermia in Milky Way-like galaxies by modeling the probability of successful travel of organic compounds between stars harboring potentially habitable planets. To this end, we apply the modified habitability recipe of Gobat & Hong (2016) to a model galaxy from the MUGS suite of zoom-in cosmological simulations. We find that, unlike habitability, which only occupies narro… ▽ More We study the process of panspermia in Milky Way-like galaxies by modeling the probability of successful travel of organic compounds between stars harboring potentially habitable planets. To this end, we apply the modified habitability recipe of Gobat & Hong (2016) to a model galaxy from the MUGS suite of zoom-in cosmological simulations. We find that, unlike habitability, which only occupies narrow dynamic range over the entire galaxy, the panspermia probability can vary be orders of magnitude between the inner ($R, b = 1-4 {\rm kpc}$) and outer disk. However, only a small fraction of star particles have very large values of panspermia probability and, consequently, the fraction of star particles where the panspermia process is more effective than prebiotic evolution is much lower than from naïve expectations based on the ratio between panspermia probability and natural habitability. △ Less

Submitted 18 September, 2021; originally announced September 2021.

Comments: 22 pages, 18 figures, ApJ accepted

Weak-lensing Mass Reconstruction of Galaxy Clusters with Convolutional Neural Network

Authors: Sungwook E. Hong, Sangnam Park, M. James Jee, Dongsu Bak, Sangjun Cha

Abstract: We introduce a novel method for reconstructing the projected matter distributions of galaxy clusters with weak-lensing (WL) data based on convolutional neural network (CNN). Training datasets are generated with ray-tracing through cosmological simulations. We control the noise level of the galaxy shear catalog such that it mimics the typical properties of the existing ground-based WL observations… ▽ More We introduce a novel method for reconstructing the projected matter distributions of galaxy clusters with weak-lensing (WL) data based on convolutional neural network (CNN). Training datasets are generated with ray-tracing through cosmological simulations. We control the noise level of the galaxy shear catalog such that it mimics the typical properties of the existing ground-based WL observations of galaxy clusters. We find that the mass reconstruction by our multi-layered CNN with the architecture of alternating convolution and trans-convolution filters significantly outperforms the traditional reconstruction methods. The CNN method provides better pixel-to-pixel correlations with the truth, restores more accurate positions of the mass peaks, and more efficiently suppresses artifacts near the field edges. In addition, the CNN mass reconstruction lifts the mass-sheet degeneracy when applied to our projected cluster mass estimation from sufficiently large fields. This implies that this CNN algorithm can be used to measure cluster masses in a model-independent way for future wide-field WL surveys. △ Less

Submitted 14 December, 2021; v1 submitted 10 February, 2021; originally announced February 2021.

Comments: 18 pages, 13 figures, ApJ accepted

Journal ref: ApJ 923, 266 (2021)

Cosmological Parameter Estimation from the Two-Dimensional Genus Topology -- Measuring the Expansion History using the Genus Amplitude as a Standard Ruler

Authors: Stephen Appleby, Changbom Park, Sungwook E. Hong, Ho Seong Hwang, Juhan Kim, Motonari Tonegawa

Abstract: We measure the genus of the galaxy distribution in two-dimensional slices of the SDSS-III BOSS catalog to constrain the cosmological parameters governing the expansion history of the Universe. The BOSS catalogs are divided into twelve concentric shells over the redshift range $0.25 < z < 0.6$ and we repeatedly measure the genus from the two-dimensional galaxy density fields, each time varying the… ▽ More We measure the genus of the galaxy distribution in two-dimensional slices of the SDSS-III BOSS catalog to constrain the cosmological parameters governing the expansion history of the Universe. The BOSS catalogs are divided into twelve concentric shells over the redshift range $0.25 < z < 0.6$ and we repeatedly measure the genus from the two-dimensional galaxy density fields, each time varying the cosmological parameters used to infer the distance-redshift relation to the shells. We also indirectly reconstruct the two-dimensional genus amplitude using the three-dimensional genus measured from SDSS Main Galaxy Sample with galaxies at low redshift $z < 0.12$. We combine the low- and high-redshift measurements, finding the cosmological model which minimizes the redshift evolution of the genus amplitude, using the fact that this quantity should be conserved. Being a distance measure, the test is sensitive to the matter density parameter ($Ω_{\rm m}$) and equation of state of dark energy ($w_{\rm de}$). We find a constraint of $w_{\rm de} = -1.05^{+0.13}_{-0.12}$, $Ω_{\rm m} = 0.303 \pm 0.036$ after combining the high- and low-redshift measurements and combining with Planck CMB data. Higher redshift data and combining data sets at low redshift will allow for stronger constraints. △ Less

Submitted 2 February, 2021; originally announced February 2021.

Comments: Accepted for publication in the Astrophysical Journal. 19 pages, 13 figures

Revealing the Local Cosmic Web from Galaxies by Deep Learning

Authors: Sungwook E. Hong, Donghui Jeong, Ho Seong Hwang, Juhan Kim

Abstract: The 80% of the matter in the Universe is in the form of dark matter that comprises the skeleton of the large-scale structure called the Cosmic Web. As the Cosmic Web dictates the motion of all matter in galaxies and inter-galactic media through gravity, knowing the distribution of dark matter is essential for studying the large-scale structure. However, the Cosmic Web's detailed structure is unkno… ▽ More The 80% of the matter in the Universe is in the form of dark matter that comprises the skeleton of the large-scale structure called the Cosmic Web. As the Cosmic Web dictates the motion of all matter in galaxies and inter-galactic media through gravity, knowing the distribution of dark matter is essential for studying the large-scale structure. However, the Cosmic Web's detailed structure is unknown because it is dominated by dark matter and warm-hot inter-galactic media, both of which are hard to trace. Here we show that we can reconstruct the Cosmic Web from the galaxy distribution using the convolutional-neural-network-based deep-learning algorithm. We find the mapping between the position and velocity of galaxies and the Cosmic Web using the results of the state-of-the-art cosmological galaxy simulations, Illustris-TNG. We confirm the mapping by applying it to the EAGLE simulation. Finally, using the local galaxy sample from Cosmicflows-3, we find the dark-matter map in the local Universe. We anticipate that the local dark-matter map will illuminate the studies of the nature of dark matter and the formation and evolution of the Local Group. High-resolution simulations and precise distance measurements to local galaxies will improve the accuracy of the dark-matter map. △ Less

Submitted 17 March, 2021; v1 submitted 4 August, 2020; originally announced August 2020.

Comments: 18 pages, 9 figures, ApJ accepted, added more data analyses

Deep-Learning Study of the 21cm Differential Brightness Temperature During the Epoch of Reionization

Authors: Yungi Kwon, Sungwook E. Hong, Inkyu Park

Abstract: We propose a deep learning analyzing technique with convolutional neural network (CNN) to predict the evolutionary track of the Epoch of Reionization (EoR) from the 21-cm differential brightness temperature tomography images. We use 21cmFAST, a fast semi-numerical cosmological 21-cm signal simulator, to produce mock 21-cm maps between $z=6 \sim 13$. We then apply two observational effects into tho… ▽ More We propose a deep learning analyzing technique with convolutional neural network (CNN) to predict the evolutionary track of the Epoch of Reionization (EoR) from the 21-cm differential brightness temperature tomography images. We use 21cmFAST, a fast semi-numerical cosmological 21-cm signal simulator, to produce mock 21-cm maps between $z=6 \sim 13$. We then apply two observational effects into those 21-cm maps, such as instrumental noise and limit of (spatial and depth) resolution somewhat suitable for realistic choices of the Square Kilometre Array (SKA). We design our deep learning model with CNN to predict the sliced-averaged neutral hydrogen fraction from the given 21-cm map. The estimated neutral fraction from our CNN model has a great agreement with its true value even after coarsely smoothing with broad beamsize and frequency bandwidth, and also heavily covered by noise with narrow. Our results have shown that deep learning analyzing method has a large potential to efficiently reconstruct the EoR history from the 21-cm tomography surveys in future. △ Less

Submitted 14 July, 2020; v1 submitted 11 June, 2020; originally announced June 2020.

Comments: 25 pages, 8 figures, added references

Journal ref: Journal of the Korean Physical Society 77, 49-59 (2020)

The Horizon Run 5 Cosmological Hydrodynamic Simulation: Probing Galaxy Formation from Kilo- to Giga-parsec Scales

Authors: Jaehyun Lee, Jihye Shin, Owain N. Snaith, Yonghwi Kim, C. Gareth Few, Julien Devriendt, Yohan Dubois, Leah M. Cox, Sungwook E. Hong, Oh-Kyoung Kwon, Chan Park, Christophe Pichon, Juhan Kim, Brad K. Gibson, Changbom Park

Abstract: Horizon Run 5 (HR5) is a cosmological hydrodynamical simulation which captures the properties of the Universe on a Gpc scale while achieving a resolution of 1kpc. Inside the simulation box we zoom-in on a high-resolution cuboid region with a volume of $1049\times119\times127\,{\rm cMpc}^3$.The sub-grid physics chosen to model galaxy formation includes radiative heating/cooling, UV background, star… ▽ More Horizon Run 5 (HR5) is a cosmological hydrodynamical simulation which captures the properties of the Universe on a Gpc scale while achieving a resolution of 1kpc. Inside the simulation box we zoom-in on a high-resolution cuboid region with a volume of $1049\times119\times127\,{\rm cMpc}^3$.The sub-grid physics chosen to model galaxy formation includes radiative heating/cooling, UV background, star formation, supernova feedback, chemical evolution tracking the enrichment of oxygen and iron, the growth of supermassive black holes and feedback from active galactic nuclei (AGN) in the form of a dual jet-heating mode. For this simulation we implemented a hybrid MPI-OMP version of RAMSES, specifically targeted for modern many-core many thread parallel architectures. In addition to the traditional simulation snapshots, light-cone data was generated on the fly. For the post-processing, we extended the Friends-of-Friend (FoF) algorithm and developed a new galaxy finder PGalF to analyse the outputs of HR5. The simulation successfully reproduces observations, such as the cosmic star formation history and connectivity of galaxy distribution, We identify cosmological structures at a wide range of scales, from filaments with a length of several cMpc, to voids with a radius of ~100 cMpc. The simulation also indicates that hydrodynamical effects on small scales impact galaxy clustering up to very large scales near and beyond the baryonic acoustic oscillation (BAO) scale. Hence, caution should be taken when using that scale as a cosmic standard ruler: one needs to carefully understand the corresponding biases. The simulation is expected to be an invaluable asset for the interpretation of upcoming deep surveys of the Universe. △ Less

Submitted 5 October, 2023; v1 submitted 1 June, 2020; originally announced June 2020.

Comments: The dimension of the zoom-in region has been corrected

Cosmological Information from the Small-scale Redshift Space Distortions

Authors: Motonari Tonegawa, Changbom Park, Yi Zheng, Hyunbae Park, Sungwook E. Hong, Ho Seong Hwang, Juhan Kim

Abstract: The redshift-space distortion (RSD) in the observed distribution of galaxies is known as a powerful probe of cosmology. Observations of large-scale RSD have given tight constraints on the linear growth rate of the large-scale structures in the universe. On the other hand, the small-scale RSD, caused by galaxy random motions inside clusters, has not been much used in cosmology, but also has cosmolo… ▽ More The redshift-space distortion (RSD) in the observed distribution of galaxies is known as a powerful probe of cosmology. Observations of large-scale RSD have given tight constraints on the linear growth rate of the large-scale structures in the universe. On the other hand, the small-scale RSD, caused by galaxy random motions inside clusters, has not been much used in cosmology, but also has cosmological information because universes with different cosmological parameters have different halo mass functions and virialized velocities. We focus on the projected correlation function $w(r_p)$ and the multipole moments $ξ_l$ on small scales ($1.4$ to $30\ h^{-1}\rm{Mpc}$). Using simulated galaxy samples generated from a physically motivated most bound particle (MBP)-galaxy correspondence scheme in the Multiverse Simulation, we examine the dependence of the small-scale RSD on the cosmological matter density parameter $��_m$, the satellite velocity bias with respect to MBPs, $b_v^s$, and the merger-time-scale parameter $α$. We find that $α=1.5$ gives an excellent fit to the $w(r_p)$ and $ξ_l$ measured from the SDSS-KIAS value added galaxy catalog. We also define the ``strength'' of Fingers-of-God as the ratio of the parallel and perpendicular size of the contour in the two-point correlation function set by a specific threshold value and show that the strength parameter helps constraining $(Ω_m, b_v^s, α)$ by breaking the degeneracy among them. The resulting parameter values from all measurements are $(Ω_m,b_v^s)=(0.272\pm0.013,0.982\pm0.040)$, indicating a slight reduction of satellite galaxy velocity relative to the MBP. However, considering that the average MBP speed inside haloes is $0.94$ times the dark matter velocity dispersion, the main drivers behind the galaxy velocity bias are gravitational interactions, rather than baryonic effects. △ Less

Submitted 25 May, 2020; originally announced May 2020.

Comments: 18 pages, 12 figures, accepted for publication in ApJ

SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES) IV: Spatial clustering and halo masses of 450-$μ$m-selected sub-millimeter galaxies

Authors: Chen-Fatt Lim, Chian-Chou Chen, Ian Smail, Wei-Hao Wang, Wei-Leong Tee, Yen-Ting Lin, Douglas Scott, Yoshiki Toba, Yu-Yen Chang, YiPing Ao, Arif Babul, Andy Bunker, Scott C. Chapman, David L Clements, Christopher J. Conselice, Yu Gao, Thomas R. Greve, Luis C. Ho, Sungwook E. Hong, Ho Seong Hwang, Maciej Koprowski, Michał J. Michałowski, Hyunjin Shim, Xinwen Shu, James M. Simpson

Abstract: We analyze an extremely deep 450-$μ$m image ($1σ=0.56$\,mJy\,beam$^{-1}$) of a $\simeq 300$\,arcmin$^{2}$ area in the CANDELS/COSMOS field as part of the SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES). We select a robust (signal-to-noise ratio $\geqslant 4$) and flux-limited ($\geqslant 4$\,mJy) sample of 164 sub-millimeter galaxies (SMGs) at 450-$μ$m that have $K$-band counterparts in the COSMOS… ▽ More We analyze an extremely deep 450-$μ$m image ($1σ=0.56$\,mJy\,beam$^{-1}$) of a $\simeq 300$\,arcmin$^{2}$ area in the CANDELS/COSMOS field as part of the SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES). We select a robust (signal-to-noise ratio $\geqslant 4$) and flux-limited ($\geqslant 4$\,mJy) sample of 164 sub-millimeter galaxies (SMGs) at 450-$μ$m that have $K$-band counterparts in the COSMOS2015 catalog identified from radio or mid-infrared imaging. Utilizing this SMG sample and the 4705 $K$-band-selected non-SMGs that reside within the noise level $\leqslant 1$\,mJy\,beam$^{-1}$ region of the 450-$μ$m image as a training set, we develop a machine-learning classifier using $K$-band magnitude and color-color pairs based on the thirteen-band photometry available in this field. We apply the trained machine-learning classifier to the wider COSMOS field (1.6\,deg$^{2}$) using the same COSMOS2015 catalog and identify a sample of 6182 450-$μ$m SMG candidates with similar colors. The number density, radio and/or mid-infrared detection rates, redshift and stellar mass distributions, and the stacked 450-$μ$m fluxes of these SMG candidates, from the S2COSMOS observations of the wide field, agree with the measurements made in the much smaller CANDELS field, supporting the effectiveness of the classifier. Using this 450-$μ$m SMG candidate sample, we measure the two-point autocorrelation functions from $z=3$ down to $z=0.5$. We find that the 450-$μ$m SMG candidates reside in halos with masses of $\simeq (2.0\pm0.5) \times10^{13}\,h^{-1}\,\rm M_{\odot}$ across this redshift range. We do not find evidence of downsizing that has been suggested by other recent observational studies. △ Less

Submitted 2 June, 2020; v1 submitted 28 April, 2020; originally announced April 2020.

Comments: ApJ accepted

Cosmological Parameter Estimation from the Two-Dimensional Genus Topology -- Measuring the Shape of the Matter Power Spectrum

Authors: Stephen A. Appleby, Changbom Park, Sungwook E. Hong, Ho Seong Hwang, Juhan Kim

Abstract: We present measurements of the two-dimensional genus of the SDSS-III BOSS catalogs to constrain cosmological parameters governing the shape of the matter power spectrum. The BOSS data are divided into twelve concentric shells over the redshift range $0.2 < z < 0.6$, and we extract the genus from the projected two-dimensional galaxy density fields. We compare the genus amplitudes to their Gaussian… ▽ More We present measurements of the two-dimensional genus of the SDSS-III BOSS catalogs to constrain cosmological parameters governing the shape of the matter power spectrum. The BOSS data are divided into twelve concentric shells over the redshift range $0.2 < z < 0.6$, and we extract the genus from the projected two-dimensional galaxy density fields. We compare the genus amplitudes to their Gaussian expectation values, exploiting the fact that this quantity is relatively insensitive to non-linear gravitational collapse. The genus amplitude provides a measure of the shape of the linear matter power spectrum, and is principally sensitive to $Ω_{\rm c}h^{2}$ and scalar spectral index $n_{\rm s}$. A strong negative degeneracy between $Ω_{\rm c}h^{2}$ and $n_{\rm s}$ is observed, as both can increase small scale power by shifting the peak and tilting the power spectrum respectively. We place a constraint on the particular combination $n_{\rm s}^{3/2} Ω_{\rm c}h^{2}$ -- we find $n_{\rm s}^{3/2} Ω_{\rm c}h^{2} = 0.1121 \pm 0.0043$ after combining the LOWZ and CMASS data sets, assuming a flat $Λ$CDM cosmology. This result is practically insensitive to reasonable variations of the power spectrum amplitude and linear galaxy bias. Our results are consistent with the Planck best fit $n_{\rm s}^{3/2}Ω_{\rm c}h^{2} = 0.1139 \pm 0.0009$. △ Less

Submitted 3 April, 2020; originally announced April 2020.

Comments: ApJ submitted, comments welcome

Constraints on the Diffuse Flux of Ultra-High Energy Neutrinos from Four Years of Askaryan Radio Array Data in Two Stations

Authors: ARA Collaboration, P. Allison, S. Archambault, J. J. Beatty, M. Beheler-Amass, D. Z. Besson, M. Beydler, C. C. Chen, C. H. Chen, P. Chen, B. A. Clark, W. Clay, A. Connolly, L. Cremonesi, J. Davies, S. de Kockere, K. D. de Vries, C. Deaconu, M. Duvernois, E. Friedman, R. Gaior, J. Hanson, K. Hanson, K. D. Hoffman, B. Hokanson-Fasig , et al. (49 additional authors not shown)

Abstract: The Askaryan Radio Array (ARA) is an ultra-high energy (UHE, $>10^{17}$ eV) neutrino detector designed to observe neutrinos by searching for the radio waves emitted by the relativistic products of neutrino-nucleon interactions in Antarctic ice. In this paper, we present constraints on the diffuse flux of ultra-high energy neutrinos between $10^{16}-10^{21}$ eV resulting from a search for neutrinos… ▽ More The Askaryan Radio Array (ARA) is an ultra-high energy (UHE, $>10^{17}$ eV) neutrino detector designed to observe neutrinos by searching for the radio waves emitted by the relativistic products of neutrino-nucleon interactions in Antarctic ice. In this paper, we present constraints on the diffuse flux of ultra-high energy neutrinos between $10^{16}-10^{21}$ eV resulting from a search for neutrinos in two complementary analyses, both analyzing four years of data (2013-2016) from the two deep stations (A2, A3) operating at that time. We place a 90 % CL upper limit on the diffuse all flavor neutrino flux at $10^{18}$ eV of $EF(E)=5.6\times10^{-16}$ $\textrm{cm}^{-2}$$\textrm{s}^{-1}$$\textrm{sr}^{-1}$. This analysis includes four times the exposure of the previous ARA result, and represents approximately 1/5 the exposure expected from operating ARA until the end of 2022. △ Less

Submitted 20 July, 2020; v1 submitted 2 December, 2019; originally announced December 2019.

Comments: 12 pages, 7 figures; Accepted to PRD

Journal ref: Phys. Rev. D 102, 043021 (2020)

Long-baseline horizontal radio-frequency transmission through polar ice

Authors: P. Allison, S. Archambault, J. J. Beatty, D. Z. Besson, C. C. Chen, C. H. Chen, P. Chen, A. Christenson, B. A. Clark, W. Clay, A. Connolly, L. Cremonesi, C. Deaconu, M. Duvernois, L. Friedman, R. Gaior, J. Hanson, K. Hanson, J. Haugen, K. D. Hoffman, E. Hong, S. Y. Hsu, L. Hu, J. J. Huang, A. M. -H. Huang , et al. (41 additional authors not shown)

Abstract: We report on analysis of englacial radio-frequency (RF) pulser data received over horizontal baselines of 1--5 km, based on broadcasts from two sets of transmitters deployed to depths of up to 1500 meters at the South Pole. First, we analyze data collected usingtwo RF bicone transmitters 1400 meters below the ice surface, and frozen into boreholes drilled for the IceCube experiment in 2011. Additi… ▽ More We report on analysis of englacial radio-frequency (RF) pulser data received over horizontal baselines of 1--5 km, based on broadcasts from two sets of transmitters deployed to depths of up to 1500 meters at the South Pole. First, we analyze data collected usingtwo RF bicone transmitters 1400 meters below the ice surface, and frozen into boreholes drilled for the IceCube experiment in 2011. Additionally, in Dec., 2018, a fat-dipole antenna, fed by one of three high-voltage (~1 kV), fast (~(1-5 ns)) signal generators was lowered into the 1700-m deep icehole drilled for the South Pole Ice Core Experiment (SPICE), approximately 3 km from the geographic South Pole. Signals from transmitters were recorded on the five englacial multi-receiver ARA stations, with receiver depths between 60--200 m. We confirm the long, >1 km RF electric field attenuation length, test our observed signal arrival timing distributions against models, and measure birefringent asymmetries at the 0.15% level. △ Less

Submitted 19 January, 2021; v1 submitted 28 August, 2019; originally announced August 2019.

Recent Results from The Askaryan Radio Array

Authors: ARA Collaboration, P. Allison, S. Archambault, R. Bard, J. J. Beatty, M. Beheler-Amass, D. Z. Besson, M. Beydler, C. -C. Chen, C. -H. Chen, P. Chen, B. Clark, A. Clough, A. Connolly, J. Davies, C. Deaconu, M. A. DuVernois, C. Fender, E. Friedman, J. Hanson, K. Hanson, J. Haugen, K. D. Hoffman, E. Hong, S. -Y. Hsu , et al. (39 additional authors not shown)

Abstract: The Askaryan Radio Array (ARA) is an ultra-high energy (UHE) neutrino telescope at the South Pole consisting of an array of radio antennas aimed at detecting the Askaryan radiation produced by neutrino interactions in the ice. Currently, the experiment has five stations in operation that have been deployed in stages since 2012. This contribution focuses on the development of a search for a diffuse… ▽ More The Askaryan Radio Array (ARA) is an ultra-high energy (UHE) neutrino telescope at the South Pole consisting of an array of radio antennas aimed at detecting the Askaryan radiation produced by neutrino interactions in the ice. Currently, the experiment has five stations in operation that have been deployed in stages since 2012. This contribution focuses on the development of a search for a diffuse flux of neutrinos in two ARA stations (A2 and A3) from 2013-2016. A background of $\sim 0.01-0.02$ events is expected in one station in each of two search channels in horizontal- and vertical-polarizations. The expected new constraints on the flux of ultra-high energy neutrinos based on four years of analysis with two stations improve on the previous limits set by ARA by a factor of about two. The projected sensitivity of ARA's five-station dataset is beginning to be competitive with other neutrino telescopes at high energies near $10^{10.5}\,$GeV. △ Less

Submitted 24 July, 2019; originally announced July 2019.

Comments: 7 pages, 5 figures, prepared for proceedings for the International Cosmic Ray Conference, Madison, WI, 2019

Alcock-Paczynski Test with the Evolution of Redshift-Space Galaxy Clustering Anisotropy

Authors: Hyunbae Park, Changbom Park, Cristiano G. Sabiu, Xiao-dong Li, Sungwook E. Hong, Juhan Kim, Motonari Tonegawa, Yi Zheng

Abstract: We develop an improved Alcock-Paczynski (AP) test method that uses the redshift-space two-point correlation function (2pCF) of galaxies. Cosmological constraints can be obtained by examining the redshift dependence of the normalized 2pCF, which should not change apart from the expected small non-linear evolution. An incorrect choice of cosmology used to convert redshift to comoving distance will m… ▽ More We develop an improved Alcock-Paczynski (AP) test method that uses the redshift-space two-point correlation function (2pCF) of galaxies. Cosmological constraints can be obtained by examining the redshift dependence of the normalized 2pCF, which should not change apart from the expected small non-linear evolution. An incorrect choice of cosmology used to convert redshift to comoving distance will manifest itself as redshift-dependent 2pCF. Our method decomposes the redshift difference of the two-dimensional correlation function into the Legendre polynomials whose amplitudes are modeled by radial fitting functions. Our likelihood analysis with this 2-D fitting scheme tightens the constraints on $Ω_m$ and ${w}$ by $\sim 40\%$ compared to the method of Li et al. (2016, 2017, 2018) that uses one dimensional angular dependence only. We also find that the correction for the non-linear evolution in the 2pCF has a non-negligible cosmology dependence, which has been neglected in previous similar studies by Li et al.. With an accurate accounting for the non-linear systematics and use of full two-dimensional shape information of the 2pCF down to scales as small as $5~h^{-1}{\rm Mpc}$ it is expected that the AP test with redshift-space galaxy clustering anisotropy can be a powerful method to constrain the expansion history of the universe. △ Less

Submitted 8 July, 2019; v1 submitted 10 April, 2019; originally announced April 2019.

Comments: 11 pages, 8 figures, accepted for publication

Constraining Cosmology with Big Data Statistics of Cosmological Graphs

Authors: Sungryong Hong, Donghui Jeong, Ho Seong Hwang, Juhan Kim, Sungwook E. Hong, Changbom Park, Arjun Dey, Milos Milosavljevic, Karl Gebhardt, Kyoung-Soo Lee

Abstract: By utilizing large-scale graph analytic tools implemented in the modern Big Data platform, Apache Spark, we investigate the topological structure of gravitational clustering in five different universes produced by cosmological $N$-body simulations with varying parameters: (1) a WMAP 5-year compatible $Λ$CDM cosmology, (2) two different dark energy equation of state variants, and (3) two different… ▽ More By utilizing large-scale graph analytic tools implemented in the modern Big Data platform, Apache Spark, we investigate the topological structure of gravitational clustering in five different universes produced by cosmological $N$-body simulations with varying parameters: (1) a WMAP 5-year compatible $Λ$CDM cosmology, (2) two different dark energy equation of state variants, and (3) two different cosmic matter density variants. For the Big Data calculations, we use a custom build of stand-alone Spark/Hadoop cluster at Korea Institute for Advanced Study (KIAS) and Dataproc Compute Engine in Google Cloud Platform (GCP) with the sample size ranging from 7 millions to 200 millions. We find that among the many possible graph-topological measures, three simple ones: (1) the average of number of neighbors (the so-called average vertex degree) $α$, (2) closed-to-connected triple fraction (the so-called transitivity) $τ_Δ$, and (3) the cumulative number density $n_{s\ge5}$ of subcomponents with connected component size $s \ge 5$, can effectively discriminate among the five model universes. Since these graph-topological measures are in direct relation with the usual $n$-points correlation functions of the cosmic density field, graph-topological statistics powered by Big Data computational infrastructure opens a new, intuitive, and computationally efficient window into the dark Universe. △ Less

Submitted 18 March, 2019; originally announced March 2019.

Comments: 16 pages, 11 figures, submitted to MNRAS

DOTIFS: spectrograph optical and opto-mechanical design

Authors: Haeun Chung, A. N. Ramaprakash, Pravin Khodade, Deepa Modi, Chaitanya V. Rajarshi, Sabyasachi Chattopadhyay, Pravin A. Chordia, Vishal Joshi, Sungwook E. Hong, Amitesh Omar, Swara Ravindranath, Yong-Sun Park, Changbom Park

Abstract: Devasthal Optical Telescope Integral Field Spectrograph (DOTIFS) is a new multi-Integral Field Unit (IFU) instrument, planned to be mounted on the 3.6m Devasthal optical telescope in Nainital, India. It has eight identical, fiber-fed spectrographs to disperse light coming from 16 IFUs. The spectrographs produce 2,304 spectra over a 370-740nm wavelength range simultaneously with a spectral resoluti… ▽ More Devasthal Optical Telescope Integral Field Spectrograph (DOTIFS) is a new multi-Integral Field Unit (IFU) instrument, planned to be mounted on the 3.6m Devasthal optical telescope in Nainital, India. It has eight identical, fiber-fed spectrographs to disperse light coming from 16 IFUs. The spectrographs produce 2,304 spectra over a 370-740nm wavelength range simultaneously with a spectral resolution of R=1200-2400. It is composed of all-refractive, all spherical optics designed to achieve on average 26.0% throughput from the telescope to the CCD with the help of high transmission spectrograph optics, volume phase holographic grating, and graded coated e2v 2K by 4K CCD. We present the optical and opto-mechanical design of the spectrograph as well as current development status. Optics and optomechanical components for the spectrographs are being fabricated. △ Less

Submitted 8 October, 2018; originally announced October 2018.

Comments: 14 pages, 10 figures, Conference poster at SPIE Astronomical Telescopes + Instrumentation, 2018, Austin, Texas, United States

Journal ref: Proc. SPIE 10702, Ground-based and Airborne Instrumentation for Astronomy VII, 107027A (8 July 2018)

Noise Characterization of IUCAA Digital Sampling Array Controller (IDSAC)

Authors: Sabyasachi Chattopadhyay, A. N. Ramaprakash, Bhushan Joshi, Pravin A. Chordia, Mahesh P. Burse, Kalpesh Chillal, Sakya Sinha, Sujit Punnadi, Ketan Rikame, Sungwook E. Hong, Dhruv Paranjpye, Haeun Chung, Changbom Park, Amitesh Omar

Abstract: IUCAA Digital Sampling Array Controller (IDSAC) is a flexible and generic yet powerful CCD controller which can handle a wide range of scientific detectors. Based on an easily scalable modular backplane architecture consisting of Single Board Controllers (SBC), IDSAC can control large detector arrays and mosaics. Each of the SBCs offers the full functionality required to control a CCD independentl… ▽ More IUCAA Digital Sampling Array Controller (IDSAC) is a flexible and generic yet powerful CCD controller which can handle a wide range of scientific detectors. Based on an easily scalable modular backplane architecture consisting of Single Board Controllers (SBC), IDSAC can control large detector arrays and mosaics. Each of the SBCs offers the full functionality required to control a CCD independently. The SBCs can be cold swapped without the need to reconfigure them. Each SBC can handle data from up to four video channels with or without dummy outputs at speeds up to 500 kilo Pixels Per Second (kPPS) Per Channel with a resolution of 16 bits. Communication with Linux based host computer is through a USB3.0 interface, with the option of using copper or optical fibers. A Field Programmable Gate Array (FPGA) is used as the master controller in each SBC which allows great flexibility in optimizing performance by adjusting gain, timing signals, bias levels, etc. using user-editable configuration files without altering the circuit topology. Elimination of thermal kTC noise is achieved via Digital Correlated Double Sampling (DCDS). We present the results of noise performance characterization of IDSAC through simulation, theoretical modeling, and actual measurements. The contribution of different types of noise sources is modeled using a tool to predict noise of a generic DCDS signal chain analytically. The analytical model predicts the net input referenced noise of the signal chain to be 5 electrons for 200k pixels per second per channel readout rate with 3 samples per pixel. Using a cryogenic test set up in the lab, the noise is measured to be 5.4 e (24.3 \muV), for the same readout configuration. △ Less

Submitted 3 October, 2018; originally announced October 2018.

Comments: 21 Pages, 13 Figures, 6 Tables

Journal ref: Sabyasachi Chattopadhyay et al., Noise Characterization of IUCAA Digital Sampling Array Controller (IDSAC), J. Astron. Telesc. Instrum. Syst. 4(3), 036002 (2018)

Multiple rooks of chess - a generic integral field unit deployment technique

Authors: Sabyasachi Chattopadhyay, A. N. Ramaprakash, Pravin Khodade, Kabir Chakrabarty, Shabbir Shaikh, Haeun Chung, Sungwook E. Hong

Abstract: A new field re-configuration technique, Multiple Rooks of Chess (MRC), for multiple deployable Integral Field Spectrographs has been developed. The method involves mechanical geometry as well as an optimized deployment algorithm. The geometry is found to be simple for mechanical implementation. The algorithm initially assigns the IFUs to the target objects and then devises the movement sequence ba… ▽ More A new field re-configuration technique, Multiple Rooks of Chess (MRC), for multiple deployable Integral Field Spectrographs has been developed. The method involves mechanical geometry as well as an optimized deployment algorithm. The geometry is found to be simple for mechanical implementation. The algorithm initially assigns the IFUs to the target objects and then devises the movement sequence based on the current and the desired IFU positions. The reconfiguration time using the suitable actuators which runs at 20 cm/s is found to be a maximum of 25 seconds for the circular DOTIFS focal plane (180 mm diameter). The Geometry Algorithm Combination (GAC) has been tested on several million mock target configurations with object-to-IFU (τ ) ratio varying from 0.25 to 16. The MRC method is found to-be efficient in target acquisition in terms of field revisit and deployment time without any collision or entanglement of the fiber bundles. The efficiency of the technique does not get affected by the increase in number density of target objects. The technique is compared with other available methods based on sky coverage, flexibility and overhead time. The proposed geometry and algorithm combination is found to have an advantage in all of the aspects. △ Less

Submitted 2 October, 2018; originally announced October 2018.

Comments: 18 Pages, 13 Figures, 1 Table

Journal ref: Proc. SPIE 10706, Advances in Optical and Mechanical Technologies for Telescopes and Instrumentation III, 107065S (10 July 2018)

Design and Performance of an Interferometric Trigger Array for Radio Detection of High-Energy Neutrinos

Authors: P. Allison, S. Archambault, R. Bard, J. J. Beatty, M. Beheler-Amass, D. Z. Besson, M. Beydler, M. Bogdan, C. -C. Chen, C. -H. Chen, P. Chen, B. A. Clark, A. Clough, A. Connolly, L. Cremonesi, J. Davies, C. Deaconu, M. A. DuVernois, E. Friedman, J. Hanson, K. Hanson, J. Haugen, K. D. Hoffman, B. Hokanson-Fasig, E. Hong , et al. (47 additional authors not shown)

Abstract: Ultra-high energy neutrinos are detectable through impulsive radio signals generated through interactions in dense media, such as ice. Subsurface in-ice radio arrays are a promising way to advance the observation and measurement of astrophysical high-energy neutrinos with energies above those discovered by the IceCube detector ($\geq$1 PeV) as well as cosmogenic neutrinos created in the GZK proces… ▽ More Ultra-high energy neutrinos are detectable through impulsive radio signals generated through interactions in dense media, such as ice. Subsurface in-ice radio arrays are a promising way to advance the observation and measurement of astrophysical high-energy neutrinos with energies above those discovered by the IceCube detector ($\geq$1 PeV) as well as cosmogenic neutrinos created in the GZK process ($\geq$100 PeV). Here we describe the $\textit{NuPhase}$ detector, which is a compact receiving array of low-gain antennas deployed 185 m deep in glacial ice near the South Pole. Signals from the antennas are digitized and coherently summed into multiple beams to form a low-threshold interferometric phased array trigger for radio impulses. The NuPhase detector was installed at an Askaryan Radio Array (ARA) station during the 2017/18 Austral summer season. $\textit{In situ}$ measurements with an impulsive, point-source calibration instrument show a 50% trigger efficiency on impulses with voltage signal-to-noise ratios (SNR) of $\le$2.0, a factor of $\sim$1.8 improvement in SNR over the standard ARA combinatoric trigger. Hardware-level simulations, validated with $\textit{in situ}$ measurements, predict a trigger threshold of an SNR as low as 1.6 for neutrino interactions that are in the far field of the array. With the already-achieved NuPhase trigger performance included in ARASim, a detector simulation for the ARA experiment, we find the trigger-level effective detector volume is increased by a factor of 1.8 at neutrino energies between 10 and 100 PeV compared to the currently used ARA combinatoric trigger. We also discuss an achievable near term path toward lowering the trigger threshold further to an SNR of 1.0, which would increase the effective single-station volume by more than a factor of 3 in the same range of neutrino energies. △ Less

Submitted 21 October, 2018; v1 submitted 12 September, 2018; originally announced September 2018.

Comments: 19 pages, 19 figures, submitted to Nuclear Instruments and Methods in Physics Research A

Observation of Reconstructable Radio Emission Coincident with an X-Class Solar Flare in the Askaryan Radio Array Prototype Station

Authors: P. Allison, S. Archambault, J. Auffenberg, R. Bard, J. J. Beatty, M. Beheler-Amass, D. Z. Besson, M. Beydler, C. Bora, C. -C. Chen, C. -H. Chen, P. Chen, B. A. Clark, A. Clough, A. Connolly, J. Davies, C. Deaconu, M. A. DuVernois, E. Friedman, B. Fox, P. W. Gorham, J. Hanson, K. Hanson, J. Haugen, B. Hill , et al. (52 additional authors not shown)

Abstract: The Askaryan Radio Array (ARA) reports an observation of radio emission coincident with the "Valentine's Day" solar flare on Feb. 15$^{\rm{th}}$, 2011 in the prototype "Testbed" station. We find $\sim2000$ events that passed our neutrino search criteria during the 70 minute period of the flare, all of which reconstruct to the location of the sun. A signal analysis of the events reveals them to be… ▽ More The Askaryan Radio Array (ARA) reports an observation of radio emission coincident with the "Valentine's Day" solar flare on Feb. 15$^{\rm{th}}$, 2011 in the prototype "Testbed" station. We find $\sim2000$ events that passed our neutrino search criteria during the 70 minute period of the flare, all of which reconstruct to the location of the sun. A signal analysis of the events reveals them to be consistent with that of bright thermal noise correlated across antennas. This is the first natural source of radio emission reported by ARA that is tightly reconstructable on an event-by-event basis. The observation is also the first for ARA to point radio from individual events to an extraterrestrial source on the sky. We comment on how the solar flares, coupled with improved systematic uncertainties in reconstruction algorithms, could aid in a mapping of any above-ice radio emission, such as that from cosmic-ray air showers, to astronomical locations on the sky. △ Less

Submitted 9 July, 2018; originally announced July 2018.

Comments: 18 pages, 22 figures

Cosmological constraints from the redshift dependence of the Alcock-Paczynski effect: Dynamical dark energy

Authors: Xiao-Dong Li, Cristiano G. Sabiu, Changbom Park, Yuting Wang, Gong-bo Zhao, Hyunbae Park, Arman Shafieloo, Juhan Kim, Sungwook E. Hong

Abstract: We perform an anisotropic clustering analysis of 1,133,326 galaxies from the Sloan Digital Sky Survey (SDSS-III) Baryon Oscillation Spectroscopic Survey (BOSS) Data Release (DR) 12 covering the redshift range $0.15<z<0.69$. The geometrical distortions of the galaxy positions, caused by incorrect cosmological model assumptions, are captured in the anisotropic two-point correlation function on scale… ▽ More We perform an anisotropic clustering analysis of 1,133,326 galaxies from the Sloan Digital Sky Survey (SDSS-III) Baryon Oscillation Spectroscopic Survey (BOSS) Data Release (DR) 12 covering the redshift range $0.15<z<0.69$. The geometrical distortions of the galaxy positions, caused by incorrect cosmological model assumptions, are captured in the anisotropic two-point correlation function on scales 6 -- 40 $h^{-1}\rm Mpc$. The redshift evolution of this anisotropic clustering is used to place constraints on the cosmological parameters. We improve the methodology of Li et al. 2016, to enable efficient exploration of high dimensional cosmological parameter spaces, and apply it to the Chevallier-Polarski-Linder parametrization of dark energy, $w=w_0+w_a{z}/({1+z})$. In combination with the CMB, BAO, SNIa and $H_0$ from Cepheid data, we obtain $Ω_m = 0.301 \pm 0.008,\ w_0 = -1.042 \pm 0.067,\ $ and $w_a = -0.07 \pm 0.29$ (68.3\% CL). Adding our new AP measurements to the aforementioned results reduces the error bars by $\sim$30 -- 40\% and improves the dark energy figure of merit by a factor of $\sim$2. We check the robustness of the results using realistic mock galaxy catalogues. △ Less

Submitted 5 March, 2018; originally announced March 2018.

Comments: 12 pages, 9 figures, accepted to ApJ

Journal ref: ApJ, 856, 2 (2018)

Cosmological Parameter Estimation Using the Genus Amplitude - Application to Mock Galaxy Catalogs

Authors: Stephen Appleby, Changbom Park, Sungwook E. Hong, Juhan Kim

Abstract: We study the topology of the matter density field in two dimensional slices, and consider how we can use the amplitude $A$ of the genus for cosmological parameter estimation. Using the latest Horizon Run 4 simulation data, we calculate the genus of the smoothed density field constructed from lightcone mock galaxy catalogs. Information can be extracted from the amplitude of the genus by considering… ▽ More We study the topology of the matter density field in two dimensional slices, and consider how we can use the amplitude $A$ of the genus for cosmological parameter estimation. Using the latest Horizon Run 4 simulation data, we calculate the genus of the smoothed density field constructed from lightcone mock galaxy catalogs. Information can be extracted from the amplitude of the genus by considering both its redshift evolution and magnitude. The constancy of the genus amplitude with redshift can be used as a standard population, from which we derive constraints on the equation of state of dark energy $w_{\rm de}$ - by measuring $A$ at $z \sim 0.1$ and $z \sim 1$, we can place an order $Δw_{\rm de} \sim {\cal O}(15\%)$ constraint on $w_{\rm de}$. By comparing $A$ to its Gaussian expectation value we can potentially derive an additional stringent constraint on the matter density $ΔΩ_{\rm mat} \sim 0.01$. We discuss the primary sources of contamination associated with the two measurements - redshift space distortion and shot noise. With accurate knowledge of galaxy bias, we can successfully remove the effect of redshift space distortion, and the combined effect of shot noise and non-linear gravitational evolution is suppressed by smoothing over suitably large scales $R_{\rm G} \ge 15 {\rm Mpc}/h$. Without knowledge of the bias, we discuss how joint measurements of the two and three dimensional genus can be used to constrain the growth factor $β= f/b$. The method can be applied optimally to redshift slices of a galaxy distribution generated using the drop-off technique. △ Less

Submitted 19 January, 2018; originally announced January 2018.

Comments: Published 01/2018. 15 pages, 10 figures

Journal ref: ApJ 853 17 (2018)

Minkowski Tensors in Two Dimensions - Probing the Morphology and Isotropy of the Matter and Galaxy Density Fields

Authors: Stephen Appleby, Pravabati Chingangbam, Changbom Park, Sungwook E. Hong, Juhan Kim, Vidhya Ganesan

Abstract: We apply the Minkowski Tensor statistics to two dimensional slices of the three dimensional density field. The Minkowski Tensors are a set of functions that are sensitive to directionally dependent signals in the data, and furthermore can be used to quantify the mean shape of density peaks. We begin by introducing our algorithm for constructing bounding perimeters around subsets of a two dimension… ▽ More We apply the Minkowski Tensor statistics to two dimensional slices of the three dimensional density field. The Minkowski Tensors are a set of functions that are sensitive to directionally dependent signals in the data, and furthermore can be used to quantify the mean shape of density peaks. We begin by introducing our algorithm for constructing bounding perimeters around subsets of a two dimensional field, and reviewing the definition of Minkowski Tensors. Focusing on the translational invariant statistic $W^{1,1}_{2}$ - a $2 \times 2$ matrix - we calculate its eigenvalues for both the entire excursion set ($Λ_{1},Λ_{2}$) and for individual connected regions and holes within the set ($λ_{1},λ_{2}$). The ratio of eigenvalues $Λ_{2}/Λ_{1}$ informs us of the presence of global anisotropies in the data, and $\langle λ_{2}/λ_{1} \rangle$ is a measure of the mean shape of peaks and troughs in the density field. We study these quantities for a Gaussian field, then consider how they are modified by the effect of gravitational collapse using the latest Horizon Run 4 cosmological simulation. We find $Λ_{1,2}$ are essentially independent of gravitational collapse, as the process maintains statistical isotropy. However, the mean shape of peaks is modified significantly - overdensities become relatively more circular compared to underdensities of the same area. When applying the statistic to a redshift space distorted density field, we find a significant signal in the eigenvalues $Λ_{1,2}$, suggesting that they can be used to probe the large-scale velocity field. △ Less

Submitted 22 May, 2018; v1 submitted 20 December, 2017; originally announced December 2017.

Comments: 17 pages, accepted for publication in APJ

Measurement of the real dielectric permittivity epsilon_r of glacial ice

Authors: P. Allison, S. Archambault, J. Auffenberg, R. Bard, J. J. Beatty, M. Beheler-Amass, D. Z. Besson, M. Beydler, C. Brabec, C. -C. Chen, C. -H. Chen, P. Chen, A. Christenson, B. A. Clark, A. Connolly, L. Cremonesi, C. Deaconu, M. Duvernois, L. Friedman, R. Gaior, P. W. Gorham, J. Hanson, K. Hanson, J. Haugen, K. D. Hoffman , et al. (44 additional authors not shown)

Abstract: Using data collected by the Askaryan Radio Array (ARA) experiment at the South Pole, we have used long-baseline propagation of radio-frequency signals to extract information on the radio-frequency index-of-refraction in South Polar ice. Owing to the increasing ice density over the upper 150--200 meters, rays are observed along two, nearly parallel paths, one of which is direct and a second which r… ▽ More Using data collected by the Askaryan Radio Array (ARA) experiment at the South Pole, we have used long-baseline propagation of radio-frequency signals to extract information on the radio-frequency index-of-refraction in South Polar ice. Owing to the increasing ice density over the upper 150--200 meters, rays are observed along two, nearly parallel paths, one of which is direct and a second which refracts through an inflection point, with differences in both arrival time and arrival angle that can be used to constrain the neutrino properties. We also observe indications, for the first time, of radio-frequency ice birefringence for signals propagating along predominantly horizontal trajectories, corresponding to an asymmetry of order 0.1% between the ordinary and extra-ordinary paths, numerically compatible with previous measurements of birefringent asymmetries for vertically-propagating radio-frequency signals at South Pole. Taken together, these effects offer the possibility of redundantly measuring the range from receiver to a neutrino interaction in Antarctic ice, if receiver antennas are deployed at shallow (25 m<z<100 m) depths. Such range information is essential in determining both the neutrino energy, as well as the incident neutrino direction. △ Less

Submitted 14 January, 2019; v1 submitted 8 December, 2017; originally announced December 2017.

Cosmological constraints from the redshift dependence of the volume effect using the galaxy 2-point correlation function across the line-of-sight

Authors: Xiao-Dong Li, Changbom Park, Cristiano G. Sabiu, Hyunbae Park, Cheng Cheng, Juhan Kim, Sungwook E. Hong

Abstract: We develop a methodology to use the redshift dependence of the galaxy 2-point correlation function (2pCF) across the line-of-sight, $ξ(r_{\bot})$, as a probe of cosmological parameters. The positions of galaxies in comoving Cartesian space varies under different cosmological parameter choices, inducing a {\it redshift-dependent scaling} in the galaxy distribution. This geometrical distortion can b… ▽ More We develop a methodology to use the redshift dependence of the galaxy 2-point correlation function (2pCF) across the line-of-sight, $ξ(r_{\bot})$, as a probe of cosmological parameters. The positions of galaxies in comoving Cartesian space varies under different cosmological parameter choices, inducing a {\it redshift-dependent scaling} in the galaxy distribution. This geometrical distortion can be observed as a redshift-dependent rescaling in the measured $ξ(r_{\bot})$. We test this methodology using a sample of 1.75 billion mock galaxies at redshifts 0, 0.5, 1, 1.5, 2, drawn from the Horizon Run 4 N-body simulation. The shape of $ξ(r_{\bot})$ can exhibit a significant redshift evolution when the galaxy sample is analyzed under a cosmology differing from the true, simulated one. Other contributions, including the gravitational growth of structure, galaxy bias, and the redshift space distortions, do not produce large redshift evolution in the shape. We show that one can make use of this geometrical distortion to constrain the values of cosmological parameters governing the expansion history of the universe. This method could be applicable to future large scale structure surveys, especially photometric surveys such as DES, LSST, to derive tight cosmological constraints. This work is a continuation of our previous works as a strategy to constrain cosmological parameters using redshift-invariant physical quantities. △ Less

Submitted 29 June, 2017; originally announced June 2017.

Comments: 12 pages, 6 figures, Accepted for publication in ApJ

Small-scale Effects of Thermal Inflation on Halo Abundance at High-$z$, Galaxy Substructure Abundance and 21-cm Power Spectrum

Authors: Sungwook E. Hong, Heeseung Zoe, Kyungjin Ahn

Abstract: We study the impact of thermal inflation on the formation of cosmological structures and present astrophysical observables which can be used to constrain and possibly probe the thermal inflation scenario. These are dark matter halo abundance at high redshifts, satellite galaxy abundance in the Milky Way, and fluctuation in the 21-cm radiation background before the epoch of reionization. The therma… ▽ More We study the impact of thermal inflation on the formation of cosmological structures and present astrophysical observables which can be used to constrain and possibly probe the thermal inflation scenario. These are dark matter halo abundance at high redshifts, satellite galaxy abundance in the Milky Way, and fluctuation in the 21-cm radiation background before the epoch of reionization. The thermal inflation scenario leaves a characteristic signature on the matter power spectrum by boosting the amplitude at a specific wavenumber determined by the number of e-foldings during thermal inflation ($N_{\rm bc}$), and strongly suppressing the amplitude for modes at smaller scales. For a reasonable range of parameter space, one of the consequences is the suppression of minihalo formation at high redshifts and that of satellite galaxies in the Milky Way. While this effect is substantial, it is degenerate with other cosmological or astrophysical effects. The power spectrum of the 21-cm background probes this impact more directly, and its observation may be the best way to constrain the thermal inflation scenario due to the characteristic signature in the power spectrum. The Square Kilometre Array (SKA) in phase 1 (SKA1) has sensitivity large enough to achieve this goal for models with $N_{\rm bc}\gtrsim 26$ if a 10000-hr observation is performed. The final phase SKA, with anticipated sensitivity about an order of magnitude higher, seems more promising and will cover a wider parameter space. △ Less

Submitted 14 November, 2017; v1 submitted 25 June, 2017; originally announced June 2017.

Comments: 28 pages, 8 figures

Journal ref: Phys. Rev. D 96, 103515 (2017)

A question of separation: disentangling tracer bias and gravitational nonlinearity with counts-in-cells statistics

Authors: Cora Uhlemann, Martin Feix, Sandrine Codis, Christophe Pichon, Francis Bernardeau, Benjamin L'Huillier, Juhan Kim, Sungwook E. Hong, Clotilde Laigle, Changbom Park, Jihye Shin, Dmitri Pogosyan

Abstract: Starting from a very accurate model for density-in-cells statistics of dark matter based on large deviation theory, a bias model for the tracer density in spheres is formulated. It adopts a mean bias relation based on a quadratic bias model to relate the log-densities of dark matter to those of mass-weighted dark haloes in real and redshift space. The validity of the parametrised bias model is est… ▽ More Starting from a very accurate model for density-in-cells statistics of dark matter based on large deviation theory, a bias model for the tracer density in spheres is formulated. It adopts a mean bias relation based on a quadratic bias model to relate the log-densities of dark matter to those of mass-weighted dark haloes in real and redshift space. The validity of the parametrised bias model is established using a parametrisation-independent extraction of the bias function. This average bias model is then combined with the dark matter PDF, neglecting any scatter around it: it nevertheless yields an excellent model for densities-in-cells statistics of mass tracers that is parametrised in terms of the underlying dark matter variance and three bias parameters. The procedure is validated on measurements of both the one and two point statistics of subhalo densities in the state-of-the-art Horizon Run 4 simulation showing excellent agreement for measured dark matter variance and bias parameters. Finally, it is demonstrated that this formalism allows for a joint estimation of the nonlinear dark matter variance and the bias parameters using solely the statistics of subhaloes. Having verified that galaxy counts in hydrodynamical simulations sampled on a scale of 10 Mpc/h closely resemble those of subhaloes, this work provides important steps towards making theoretical predictions for density-in-cells statistics applicable to upcoming galaxy surveys like Euclid or WFIRST. △ Less

Submitted 24 May, 2017; originally announced May 2017.

Comments: 14 pages, 11 figures

CMB Spectral Distortion Constraints on Thermal Inflation

Authors: Kihyun Cho, Sungwook E. Hong, Ewan D. Stewart, Heeseung Zoe

Abstract: Thermal inflation is a second epoch of exponential expansion at typical energy scales $V^{1/4} \sim 10^{6 \sim 8} \mathrm{GeV}$. If the usual primordial inflation is followed by thermal inflation, the primordial power spectrum is only modestly redshifted on large scales, but strongly suppressed on scales smaller than the horizon size at the beginning of thermal inflation,… ▽ More Thermal inflation is a second epoch of exponential expansion at typical energy scales $V^{1/4} \sim 10^{6 \sim 8} \mathrm{GeV}$. If the usual primordial inflation is followed by thermal inflation, the primordial power spectrum is only modestly redshifted on large scales, but strongly suppressed on scales smaller than the horizon size at the beginning of thermal inflation, $k > k_{\rm b} = a_{\rm b} H_{\rm b}$. We calculate the spectral distortion of the cosmic microwave background generated by the dissipation of acoustic waves in this context. For $k_{\rm b} \ll 10^3 \mathrm{Mpc}^{-1}$, thermal inflation results in a large suppression of the $μ$-distortion amplitude, predicting that it falls well below the standard value of $μ\simeq 2\times 10^{-8}$. Thus, future spectral distortion experiments, similar to PIXIE, can place new limits on the thermal inflation scenario, constraining $k_{\rm b} \gtrsim 10^3 \mathrm{Mpc}^{-1}$ if $μ\simeq 2\times 10^{-8}$ were found. △ Less

Submitted 7 August, 2017; v1 submitted 8 May, 2017; originally announced May 2017.

Comments: 18 pages, 7 figures

Topology of Large-Scale Structures of Galaxies in Two Dimensions - Systematic Effects

Authors: Stephen Appleby, Changbom Park, Sungwook E. Hong, Juhan Kim

Abstract: We study the two-dimensional topology of the galactic distribution when projected onto two-dimensional spherical shells. Using the latest Horizon Run 4 simulation data, we construct the genus of the two-dimensional field and consider how this statistic is affected by late-time nonlinear effects -- principally gravitational collapse and redshift space distortion (RSD). We also consider systematic a… ▽ More We study the two-dimensional topology of the galactic distribution when projected onto two-dimensional spherical shells. Using the latest Horizon Run 4 simulation data, we construct the genus of the two-dimensional field and consider how this statistic is affected by late-time nonlinear effects -- principally gravitational collapse and redshift space distortion (RSD). We also consider systematic and numerical artifacts such as shot noise, galaxy bias, and finite pixel effects. We model the systematics using a Hermite polynomial expansion and perform a comprehensive analysis of known effects on the two-dimensional genus, with a view toward using the statistic for cosmological parameter estimation. We find that the finite pixel effect is dominated by an amplitude drop and can be made less than $1\%$ by adopting pixels smaller than $1/3$ of the angular smoothing length. Nonlinear gravitational evolution introduces time-dependent coefficients of the zeroth, first, and second Hermite polynomials, but the genus amplitude changes by less than $1\%$ between $z=1$ and $z=0$ for smoothing scales $R_{\rm G} > 9 {\rm Mpc/h}$. Non-zero terms are measured up to third order in the Hermite polynomial expansion when studying RSD. Differences in shapes of the genus curves in real and redshift space are small when we adopt thick redshift shells, but the amplitude change remains a significant $\sim {\cal O}(10\%)$ effect. The combined effects of galaxy biasing and shot noise produce systematic effects up to the second Hermite polynomial. It is shown that, when sampling, the use of galaxy mass cuts significantly reduces the effect of shot noise relative to random sampling. △ Less

Submitted 15 February, 2017; originally announced February 2017.

Comments: Accepted for publication in ApJ

Journal ref: ApJ 836 45 (2017)

Cosmological constraints from the redshift dependence of the Alcock-Paczynski effect: application to the SDSS-III BOSS DR12 galaxies

Authors: Xiao-Dong Li, Changbom Park, Cristiano G. Sabiu, Hyunbae Park, David H. Weinberg, Donald P. Schneider, Juhan Kim, Sungwook E. Hong

Abstract: We apply the methodology developed in \cite{Li2014,Li2015} to BOSS DR12 galaxies and derive cosmological constraints from the redshift dependence of the Alcock-Paczynski (AP) effect. The apparent anisotropy in the distribution of observed galaxies arise from two main sources, the redshift-space distortion (RSD) effect due to the galaxy peculiar velocities, and the geometric distortion when incorre… ▽ More We apply the methodology developed in \cite{Li2014,Li2015} to BOSS DR12 galaxies and derive cosmological constraints from the redshift dependence of the Alcock-Paczynski (AP) effect. The apparent anisotropy in the distribution of observed galaxies arise from two main sources, the redshift-space distortion (RSD) effect due to the galaxy peculiar velocities, and the geometric distortion when incorrect cosmological models are assumed for transforming redshift to comoving distance, known as the AP effect. Anisotropies produced by the RSD effect are, although large, maintaining a nearly uniform magnitude over a large range of redshift, while the degree of anisotropies from the AP effect varies with redshift by much larger magnitude. We split the DR12 galaxies into six redshift bins, measure the 2-point correlation function in each bin, and assess the redshift evolution of anisotropies. We obtain constraints of $Ω_m=0.290 \pm 0.053,\ \ w = -1.07 \pm 0.15$, which are comparable with the current constraints from other cosmological probes such as type Ia supernovae, cosmic microwave background, and baryon acoustic oscillation (BAO). Combining these cosmological probes with our method yield tight constraints of $ Ω_m = 0.301 \pm 0.006,\ w=-1.054 \pm 0.025$. Our method is complementary to the other large scale structure probes like BAO and topology. We expect this technique will play an important role in deriving cosmological constraints from large scale structure surveys. △ Less

Submitted 18 September, 2016; originally announced September 2016.

Comments: 21 pages, 14 figures. Accepted for publication in ApJ

AAO Starbugs: software control and associated algorithms

Authors: Nuria P. F. Lorente, Minh V. Vuong, Keith Shortridge, Tony J. Farrell, Scott Smedley, Sungwook E. Hong, Carlos Bacigalupo, Michael Goodwin, Kyler Kuehn, Christophe Satorre

Abstract: The Australian Astronomical Observatory's TAIPAN instrument deploys 150 Starbug robots to position optical fibres to accuracies of 0.3 arcsec, on a 32 cm glass field plate on the focal plane of the 1.2 m UK-Schmidt telescope. This paper describes the software system developed to control and monitor the Starbugs, with particular emphasis on the automated path-finding algorithms, and the metrology s… ▽ More The Australian Astronomical Observatory's TAIPAN instrument deploys 150 Starbug robots to position optical fibres to accuracies of 0.3 arcsec, on a 32 cm glass field plate on the focal plane of the 1.2 m UK-Schmidt telescope. This paper describes the software system developed to control and monitor the Starbugs, with particular emphasis on the automated path-finding algorithms, and the metrology software which keeps track of the position and motion of individual Starbugs as they independently move in a crowded field. The software employs a tiered approach to find a collision-free path for every Starbug, from its current position to its target location. This consists of three path-finding stages of increasing complexity and computational cost. For each Starbug a path is attempted using a simple method. If unsuccessful, subsequently more complex (and expensive) methods are tried until a valid path is found or the target is flagged as unreachable. △ Less

Submitted 8 August, 2016; originally announced August 2016.

Comments: 10 pages, to be published in Proc. SPIE 9913, Software and Cyberinfrastructure for Astronomy IV; 2016

The Most Bound Halo Particle-Galaxy Correspondence Model: Comparison between Models with Different Merger Timescales

Authors: Sungwook E. Hong, Changbom Park, Juhan Kim

Abstract: We develop a galaxy assignment scheme that populates dark matter halos with galaxies by tracing the most bound member particles (MBPs) of simulated halos. Several merger-timescale models based on analytic calculations and numerical simulations are adopted as the survival time of mock satellite galaxies. We build mock galaxy samples from halo merger data of the Horizon Run 4 $N$-body simulation fro… ▽ More We develop a galaxy assignment scheme that populates dark matter halos with galaxies by tracing the most bound member particles (MBPs) of simulated halos. Several merger-timescale models based on analytic calculations and numerical simulations are adopted as the survival time of mock satellite galaxies. We build mock galaxy samples from halo merger data of the Horizon Run 4 $N$-body simulation from $z = 12-0$. We compare group properties and two-point correlation functions (2pCFs) of mock galaxies with those of volume-limited SDSS galaxies, with $r$-band absolute magnitudes of $\mathcal{M}_r - 5 \log h < -21$ and $-20$ at $z=0$. It is found that the MBP-galaxy correspondence scheme reproduces the observed population of SDSS galaxies in massive galaxy groups ($M > 10^{14} h^{-1} M_{\odot}$) and the small-scale 2pCF ($r_{\rm p} < 10 h^{-1} {\rm Mpc}$) quite well for the majority of the merger timescale models adopted. The new scheme outperforms the previous subhalo-galaxy correspondence scheme by more than $2σ$. △ Less

Submitted 1 June, 2016; originally announced June 2016.

Comments: 23 pages, 8 figures

Journal ref: ApJ 823, 103 (2016)

Evolution of galaxy habitability

Authors: R. Gobat, S. E. Hong

Abstract: We combine a semi-analytic model of galaxy evolution with constraints on circumstellar habitable zones and the distribution of terrestrial planets to probe the suitability of galaxies of different mass and type to host habitable planets, and how it evolves with time. We find that the fraction of stars with terrestrial planets in their habitable zone (known as habitability) depends only weakly on g… ▽ More We combine a semi-analytic model of galaxy evolution with constraints on circumstellar habitable zones and the distribution of terrestrial planets to probe the suitability of galaxies of different mass and type to host habitable planets, and how it evolves with time. We find that the fraction of stars with terrestrial planets in their habitable zone (known as habitability) depends only weakly on galaxy mass, with a maximum around 4e10 Msun. We estimate that 0.7% of all stars in Milky Way type galaxies to host a terrestrial planet within their habitable zone, consistent with the value derived from Kepler observations. On the other hand, the habitability of passive galaxies is slightly but systematically higher, unless we assume an unrealistically high sensitivity of planets to supernovae. We find that the overall habitability of galaxies has not changed significantly in the last ~8 Gyr, with most of the habitable planets in local disk galaxies having formed ~1.5 Gyr before our own solar system. Finally, we expect that ~1.4e9 planets similar to present-day Earth have existed so far in our galaxy. △ Less

Submitted 7 June, 2016; v1 submitted 21 May, 2016; originally announced May 2016.

Comments: 10 pages, 5 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 592, A96 (2016)

Horizon Run 4 Simulation: Coupled Evolution of Galaxies and Large-scale Structures of the Universe

Authors: Juhan Kim, Changbom Park, Benjamin L'Huillier, Sungwook E. Hong

Abstract: The Horizon Run 4 is a cosmological $N$-body simulation designed for the study of coupled evolution between galaxies and large-scale structures of the Universe, and for the test of galaxy formation models. Using $6300^3$ gravitating particles in a cubic box of $L_{\rm box} = 3150 ~h^{-1}{\rm Mpc}$, we build a dense forest of halo merger trees to trace the halo merger history with a halo mass resol… ▽ More The Horizon Run 4 is a cosmological $N$-body simulation designed for the study of coupled evolution between galaxies and large-scale structures of the Universe, and for the test of galaxy formation models. Using $6300^3$ gravitating particles in a cubic box of $L_{\rm box} = 3150 ~h^{-1}{\rm Mpc}$, we build a dense forest of halo merger trees to trace the halo merger history with a halo mass resolution scale down to $M_s = 2.7 \times 10^{11} h^{-1}{\rm M_\odot}$. We build a set of particle and halo data, which can serve as testbeds for comparison of cosmological models and gravitational theories with observations. We find that the FoF halo mass function shows a substantial deviation from the universal form with tangible redshift evolution of amplitude and shape. At higher redshifts, the amplitude of the mass function is lower, and the functional form is shifted toward larger values of $\ln (1/σ)$. We also find that the baryonic acoustic oscillation feature in the two-point correlation function of mock galaxies becomes broader with a peak position moving to smaller scales and the peak amplitude decreasing for increasing directional cosine $μ$ compared to the linear predictions. From the halo merger trees built from halo data at 75 redshifts, we measure the half-mass epoch of halos and find that less massive halos tend to reach half of their current mass at higher redshifts. Simulation outputs including snapshot data, past lightcone space data, and halo merger data are available at http://sdss.kias.re.kr/astro/Horizon-Run4/. △ Less

Submitted 23 August, 2015; v1 submitted 20 August, 2015; originally announced August 2015.

Comments: 17 pages, 17 Figures, Accepted for Publication in JKAS. Data available at http://sdss.kias.re.kr/astro/Horizon-Run4/

Journal ref: JKAS (2015) 48, 213