Source author record

Changbom Park

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

76works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

76 published item(s)

preprint2026arXiv

A New Collision Avoidance Fiber Assignment Algorithm for Robotic Fiber Positioners in Multi-Object Spectroscopy

We present a new fiber assignment algorithm for a robotic fiber positioner system in multi-object spectroscopy. Modern fiber positioner systems typically have overlapping patrol regions, resulting in the number of observable targets being highly dependent on the fiber assignment scheme. To maximize observable targets without fiber collisions, the algorithm proceeds in three steps. First, it assigns the maximum number of targets for a given field of view without considering any collisions between fiber positioners. Then, the fibers in collision are grouped, and the algorithm finds the optimal solution resolving the collision problem within each group. We compare the results from this new algorithm with those from a simple algorithm that assigns targets in descending order of their rank by considering collisions. As a result, we could increase the overall completeness of target assignments by 10% with this new algorithm in comparison with the case using the simple algorithm in a field with 150 fibers. Our new algorithm is designed for the All-sky SPECtroscopic survey of nearby galaxies (A-SPEC) based on the K-SPEC spectrograph system, but can also be applied to similar fiber-based systems with heavily overlapping fiber positioners.

preprint2026arXiv

ODIN: Clustering Properties of Ly$α$ Blobs at $z$ $\sim$ 2.4 and 3.1

Spatially extended Ly$α$ nebulae, known as Ly$α$ blobs (LABs), are a rare population at $z > 2$ that are thought to trace proto-groups or the progenitors of massive galaxies in the present-day universe. However, their dark matter halo properties (e.g., halo mass) are still uncertain due to their rarity and strong field-to-field variation. The One-hundred-deg$^2$ DECam Imaging in Narrowbands (ODIN) survey has discovered 103 and 112 LABs in the extended ($\sim$9~\sqdeg) COSMOS field at $z\sim2.4$ and 3.1, respectively, enabling estimation of their bias and host halo masses through clustering analysis. We measure the angular auto-correlation functions (ACFs) of LABs and derive galaxy bias factors of $b$ = $4.0\pm0.8$ and $3.8\pm0.7$, corresponding to minimum halo masses of $2.8^{+3.0}_{-1.8}$ and $0.7^{+0.8}_{-0.5}\times10^{12}~M_\odot$ and median halo masses of $4.2^{+3.8}_{-2.5}$ and $1.1^{+1.1}_{-0.7}\times10^{12}~M_\odot$ at $z\sim2.4$ and 3.1, respectively. LABs occupy $\sim$11$^{+39}_{-8}$\% and $\sim$3$^{+9}_{-2}$\% of all dark matter halos above these minimum halo masses. These findings suggest that LABs inhabit massive dark matter halos, likely tracing proto-group environments that evolve into present-day massive halos ($\sim$10$^{13}~M_\odot$), where massive elliptical galaxies or galaxy groups reside, by $z=0$.

preprint2022arXiv

A novel approach for calculating galaxy rotation curves using spaxel cross-correlation and iterative smoothing

Precise measurements of the internal dynamics of galaxies have proven of great importance for understanding the internal dark matter distribution of galaxies. We present a novel method for measuring the line-of-sight (LOS) velocities across the face of galaxies by cross-correlation of spectral pixels (spaxels) and an iterative method of smoothing. On simulated data the method can accurately recover the input LOS velocities for different types of spectra (absorption line dominated, emission line dominated, and differing shapes of the continuum), and can handle stellar population radial gradients. Most important of all, it continues to provide reliable measurements of LOS velocities with reasonable uncertainties even when the spectra are very low signal-to-noise (approaching $\sim 1$), which is a challenge for traditional template-fitting approaches. We apply our method to data from a real MaNGA galaxy as a demonstration and find promising results with good precision. This novel approach can be complementary to existing methods primarily based on template fitting.

preprint2022arXiv

The Fundamental Plane Is Not a Plane: Warped Nature of the Fundamental Plane of Early-type Galaxies and Its Implication for Galaxy Formation

Based on $16,283$ early-type galaxies (ETGs) in $0.025\le z_\mathrm{spec}<0.055$ from Sloan Digital Sky Survey data, we show that the fundamental plane (FP) of ETGs is not a plane in the strict sense but is a curved surface with a twisted shape whose orthogonal direction to the surface is shifted as the central velocity dispersion ($σ_0$) or mean surface brightness within the half-light radius ($μ_e$) changes. When ETGs are divided into subsamples according to $σ_0$, the coefficient of $μ_e$ of the FP increases, whereas the zero-point of the FP decreases at higher $σ_0$. Taking the $z$ band as an example, the coefficient of $μ_e$ rises from $0.28$ to $0.36$ as $σ_0$ increases from $\sim100$ to $\sim300$ km s$^{-1}$. At the same time, the zero-point of the FP falls from $-7.5$ to $-9.0$ in the same $σ_0$ range. The consistent picture on the curved nature of the FP is also reached by inspecting changes in the FP coefficients for ETG subsamples with different $μ_e$. By examining scaling relations that are projections of the FP, we suggest that the warped nature of the FP may originate from dry merger effects that are imprinted more prominently in ETGs with higher masses.

preprint2021arXiv

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

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.

preprint2021arXiv

Evidence for Impact of Galaxy Mergers on Stellar Kinematics of Early-type Galaxies

We provide observational evidence that galaxy mergers significantly affect stellar kinematics of early-type galaxies (ETGs) such as specific stellar angular momentum within the half-light radius ($λ_{R_e}$) and kinematic misalignment ($ψ_\mathrm{mis}$), using MaNGA integral field unit spectroscopic data that are in the Stripe 82 region of the Sloan Digital Sky Survey. In this study, tidal features around ETGs, which are detected in deep coadded images, are used as direct evidence for mergers that occurred recently. In the case of ETGs that do not have dust lanes, $λ_{R_e}$ is lower in ETGs with tidal features than in those without tidal features (median $λ_{R_e}$: $0.21$ versus $0.39$) in all stellar mass and Sérsic index ranges except the most massive bin, so that the fraction of ETGs with tidal features in slow rotators is more than twice as large as that in fast rotators ($42\%$ versus $18\%$). Moreover, ETGs with tidal features have larger $ψ_\mathrm{mis}$ than those without tidal features (mean $ψ_\mathrm{mis}$: $28^\circ$ versus $15^\circ$). By contrast, ETGs with dust lanes are fast rotators, and ETGs with both dust lanes and tidal features have the highest $λ_{R_e}$ (median $λ_{R_e}$: $0.59$) among all ETG categories. In addition, ETGs with dust lanes have small $ψ_\mathrm{mis}$ regardless of the existence of tidal features ($ψ_\mathrm{mis}<7.5^\circ$). Our results can be explained if mergers with different gas fractions generate merger remnants that have different kinematic properties.

preprint2021arXiv

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

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$.

preprint2020arXiv

Cosmological Information from the Small-scale Redshift Space Distortions

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.

preprint2020arXiv

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

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$.

preprint2020arXiv

Dependence of the Fundamental Plane of Early-type Galaxies on Age and Internal Structure

We investigate the scatter in the fundamental plane (FP) of early-type galaxies (ETGs) and its dependence on age and internal structure of ETGs, using $16,283$ ETGs with $M_r\le-19.5$ and $0.025\le z<0.055$ in Sloan Digital Sky Survey data. We use the relation between the age of ETGs and photometric parameters such as color, absolute magnitude, and central velocity dispersion of ETGs and find that the scatter in the FP depends on age. The FP of old ETGs with age $\gtrsim9$ Gyrs has a smaller scatter of $\sim0.06$ dex ($\sim14\%$) while that of young ETGs with age $\lesssim6$ Gyrs has a larger scatter of $\sim0.075$ dex ($\sim17\%$). In the case of young ETGs, less compact ETGs have a smaller scatter in the FP ($\sim0.065$ dex; $\sim15\%$) than more compact ones ($\sim0.10$ dex; $\sim23\%$). On the other hand, the scatter in the FP of old ETGs does not depend on the compactness of galaxy structure. Thus, among the subpopulations of ETGs, compact young ETGs have the largest scatter in the FP. This large scatter in compact young ETGs is caused by ETGs that have low dynamical mass-to-light ratio ($M_\mathrm{dyn}/L$) and blue color in the central regions. By comparing with a simple model of the galaxy that has experienced a gas-rich major merger, we find that the scenario of recent gas-rich major merger can reasonably explain the properties of the compact young ETGs with excessive light for a given mass (low $M_\mathrm{dyn}/L$) and blue central color.

preprint2019arXiv

Constraining Cosmology with Big Data Statistics of Cosmological Graphs

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.

preprint2019arXiv

Ensemble Average of Three-Dimensional Minkowski Tensors of a Gaussian Random Field in Redshift Space

We present the ensemble expectation values for the translation invariant, rank-2 Minkowski tensors in three-dimensions, for a linearly redshift space distorted Gaussian random field. The Minkowski tensors $W^{0,2}_{1}$, $W^{0,2}_{2}$ are sensitive to global anisotropic signals present within a field, and by extracting these statistics from the low redshift matter density one can place constraints on the redshift space distortion parameter $β= f/b$. We begin by reviewing the calculation of the ensemble expectation values $\langle W^{0,2}_{1} \rangle$, $\langle W^{0,2}_{2} \rangle $ for isotropic, Gaussian random fields, then consider how these results are modified by the presence of a linearly anisotropic signal. Under the assumption that all fields remain Gaussian, we calculate the anisotropic correction due to redshift space distortion in a coordinate system aligned with the line of sight, finding inequality between the diagonal elements of $\langle W^{0,2}_{1} \rangle $, $\langle W^{0,2}_{2} \rangle $. The ratio of diagonal elements of these matrices provides a set of statistics that are sensitive only to the redshift space distortion parameter $β$. We estimate the Fisher information that can be extracted from the Minkowski tensors, and find $W^{0,2}_{1}$ is more sensitive to $β$ than $W^{0,2}_{2}$, and a measurement of $W^{0,2}_{1}$ accurate to $\sim 1\%$ can yield a $\sim 4\%$ constraint on $β$. Finally, we discuss the difference between using the matrix elements of the Minkowski tensors directly against measuring the eigenvalues. For the purposes of cosmological parameter estimation we advocate the use of the matrix elements, to avoid spurious anisotropic signals that can be generated by the eigenvalue decomposition.

preprint2019arXiv

Isolated dark matter deprived galaxies in hydrodynamical simulations: real objects or artefacts?

We searched for isolated dark matter deprived galaxies within several state-of-the-art hydrodynamical simulations: Illustris, IllustrisTNG, EAGLE, and Horizon-AGN and found a handful of promising objects in all except Horizon-AGN. While our initial goal was to study their properties and evolution, we quickly noticed that all of them were located at the edge of their respective simulation boxes. After carefully investigating these objects using the full particle data, we concluded that they are not merely caused by a problem with the algorithm identifying bound structures. We provide strong evidence that these oddballs were created from regular galaxies that get torn apart due to unphysical processes when crossing the edge of the simulation box. We show that these objects are smoking guns indicating an issue with the implementation of the periodic boundary conditions of the particle data in Illustris, IllustrisTNG, and EAGLE, which was eventually traced down to be a minor bug occurring for a very rare set of conditions.

preprint2016arXiv

HectoMAP and Horizon Run 4: Dense Structures and Voids in the Real and Simulated Universe

HectoMAP is a dense redshift survey of red galaxies covering a 53 $deg^{2}$ strip of the northern sky. HectoMAP is 97\% complete for galaxies with $r<20.5$, $(g-r)>1.0$, and $(r-i)>0.5$. The survey enables tests of the physical properties of large-scale structure at intermediate redshift against cosmological models. We use the Horizon Run 4, one of the densest and largest cosmological simulations based on the standard $Λ$ Cold Dark Matter ($Λ$CDM) model, to compare the physical properties of observed large-scale structures with simulated ones in a volume-limited sample covering 8$\times10^6$ $h^{-3}$ Mpc$^3$ in the redshift range $0.22<z<0.44$. We apply the same criteria to the observations and simulations to identify over- and under-dense large-scale features of the galaxy distribution. The richness and size distributions of observed over-dense structures agree well with the simulated ones. Observations and simulations also agree for the volume and size distributions of under-dense structures, voids. The properties of the largest over-dense structure and the largest void in HectoMAP are well within the distributions for the largest structures drawn from 300 Horizon Run 4 mock surveys. Overall the size, richness and volume distributions of observed large-scale structures in the redshift range $0.22<z<0.44$ are remarkably consistent with predictions of the standard $Λ$CDM model.

preprint2016arXiv

Probing scalar tensor theories for gravity in redshift space

We present measurements of the spatial clustering statistics in redshift space of various scalar field modified gravity simulations. We utilise the two-point and the three-point correlation functions to quantify the spatial distribution of dark matter halos within these simulations and thus discern between the models. We compare $Λ$CDM simulations to various modified gravity scenarios and find consistency with previous work in terms of 2-point statistics in real and redshift-space. However using higher order statistics such as the three-point correlation function in redshift space we find significant deviations from $Λ$CDM hinting that higher order statistics may prove to be a useful tool in the hunt for deviations from General Relativity.

preprint2016arXiv

Quasars as a tracer of large-scale structures in the distant universe

We study the dependence of the number density and properties of quasars on the background galaxy density using the currently largest spectroscopic datasets of quasars and galaxies. We construct a galaxy number density field smoothed over the variable smoothing scale of between approximately 10 and $20\,h^{-1}$Mpc over the redshift range of $0.46<z<0.59$ using the Sloan Digital Sky Survey (SDSS) Data Release 12 (DR12) Constant MASS (CMASS) galaxies. The quasar sample is prepared from the SDSS I/II DR7. We examine the correlation of incidence of quasars with the large-scale background density and dependence of quasar properties such as bolometric luminosity, black hole mass, and Eddington ratio on the large-scale density. We find a monotonic correlation between the quasar number density and large-scale galaxy number density, which is fitted well with a power law relation, $n_Q\proptoρ_G^{0.618}$. We detect weak dependences of quasar properties on the large-scale density such as a positive correlation between black hole mass and density, and a negative correlation between luminosity and density. We discuss the possibility of using quasars as a tracer of large-scale structures at high redshifts, which may be useful for studies of growth of structures in the high redshift universe.

preprint2016arXiv

The DESI Experiment Part I: Science,Targeting, and Survey Design

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$α$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

preprint2016arXiv

The DESI Experiment Part II: Instrument Design

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= λ/Δλ$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

preprint2016arXiv

The Detailed Science Case for the Maunakea Spectroscopic Explorer: the Composition and Dynamics of the Faint Universe

MSE is an 11.25m aperture observatory with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. More than 3200 fibres will feed spectrographs operating at low (R ~ 2000 - 3500) and moderate (R ~ 6000) spectral resolution, and approximately 1000 fibers will feed spectrographs operating at high (R ~ 40000) resolution. MSE is designed to enable transformational science in areas as diverse as tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of thick disk and halo stars; connecting galaxies to their large scale structure; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars; next generation cosmological surveys using redshift space distortions and peculiar velocities. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, WFIRST, PLATO, and the SKA, and is designed to complement and go beyond the science goals of other planned and current spectroscopic capabilities like VISTA/4MOST, WHT/WEAVE, AAT/HERMES and Subaru/PFS. It is an ideal feeder facility for E-ELT, TMT and GMT, and provides the missing link between wide field imaging and small field precision astronomy. MSE is optimized for high throughput, high signal-to-noise observations of the faintest sources in the Universe with high quality calibration and stability being ensured through the dedicated operational mode of the observatory. (abridged)

preprint2016arXiv

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

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σ$.

preprint2016arXiv

Two-dimensional Topology of Cosmological Reionization

We study the two-dimensional topology of the 21-cm differential brightness temperature for two hydrodynamic radiative transfer simulations and two semi-numerical models. In each model, we calculate the two dimensional genus curve for the early, middle and late epochs of reionization. It is found that the genus curve depends strongly on the ionized fraction of hydrogen in each model. The genus curves are significantly different for different reionization scenarios even when the ionized faction is the same. We find that the two-dimensional topology analysis method is a useful tool to constrain the reionization models. Our method can be applied to the future observations such as those of the Square Kilometer Array.

preprint2015arXiv

Cosmological constraints from the redshift dependence of the Alcock-Paczynski test and volume effect: galaxy two-point correlation function

We propose a method using the redshift dependence of the Alcock-Paczynski (AP) test and volume effect to measure the cosmic expansion history. The galaxy two-point correlation function as a function of angle, $ξ(μ)$, is measured at different redshifts. Assuming an incorrect cosmological model to convert galaxy redshifts to distances, the shape of $ξ(μ)$ appears anisotropic due to the AP effect, and the amplitude shifted by the change in comoving volume. Due to the redshift dependence of the AP and volume effect, both the shape and amplitude of $ξ(μ)$ exhibit redshift dependence. Similar to Li et.al (2014), we find the redshift-space distortions (RSD) caused by galaxy peculiar velocities, although significantly distorts $ξ(μ)$, exhibit much less redshift evolution compared to the AP and volume effects. By focusing on the redshift dependence of $ξ(μ)$, we can correctly recover the cosmological parameters despite the contamination of RSD. The method is tested by using the Horizon Run 3 N-body simulation, from which we made a series of $1/8$-sky mock surveys having 8 million physically self-bound halos and sampled to have roughly a uniform number density in $z=0-1.5$. We find the AP effect results in tight, unbiased constraints on the density parameter and dark energy equation of state, with 68.3% CL intervals $δΩ_m\sim0.03$ and $δw\sim0.1$, and the volume effect leads to much tighter constraints of $δΩ_m\sim0.007$ and $δw\sim0.035$.

preprint2015arXiv

Cosmological constraints from the redshift dependence of the Alcock-Paczynski test: galaxy density gradient field

We propose a method based on the redshift dependence of the Alcock-Paczynski (AP) test to measure the expansion history of the Universe. It uses the isotropy of the galaxy density gradient field to constrain cosmological parameters. If the density parameter $Ω_m$ or the dark energy equation of state $w$ are incorrectly chosen, the gradient field appears to be anisotropic with the degree of anisotropy varying with redshift. We use this effect to constrain the cosmological parameters governing the expansion history of the Universe. Although redshift-space distortions (RSD) induced by galaxy peculiar velocities also produce anisotropies in the gradient field, these effects are close to uniform in magnitude over a large range of redshift. This makes the redshift variation of the gradient field anisotropy relatively insensitive to the RSD. By testing the method on mock surveys drawn from the Horizon Run 3 cosmological N-body simulations, we demonstrate that the cosmological parameters can be estimated without bias. Our method is complementary to the baryon acoustic oscillation or topology methods as it depends on $D_AH$, the product of the angular diameter distance and the Hubble parameter.

preprint2015arXiv

Dark matter halos of barred disk galaxies

We use a large volume-limited sample of disk galaxies drawn from the Sloan Digital Sky Survey Data Release 7 to study the dependence of the bar fraction on the stellar-to-halo mass ratio, making use of a group catalog, we identify central and satellite galaxies in our sample. For the central galaxies in the sample we estimate the stellar-to-halo mass ratio (M$_{\mathrm{*}}/$M$_{\mathrm{h}}$) and find that the fraction of barred galaxies is a strong function of this ratio, especially for the case of strong bars. Bars are more common in galaxies with high M$_{\mathrm{*}}/$M$_{\mathrm{h}}$ values, as expected from early theoretical works that showed that systems with massive dark matter halos are more stable against bar instabilities. We find that the change of the bar fraction with M$_{\mathrm{h}}$ and M$_{\mathrm{*}}$ is stronger if we consider a relation with the form $f_{\mathrm{bar}}=f_{\mathrm{bar}}$(M$_{\mathrm{*}}^α/$M$_{\mathrm{h}}$) with $α=1.5$, and that the bar fraction is largely independent of other physical properties such as color and spin parameter when M$_{\mathrm{*}}^{3/2}/$M$_{\mathrm{h}}$ is fixed. With our sample of galaxies segregated into centrals and satellites, we also compare the fraction of barred galaxies in each group, finding a slightly higher bar fraction for satellites when compared with centrals at fixed stellar mass, but at fixed color this difference becomes very weak. This result, in agreement with previous studies, confirms that the bar fraction does not directly depend on the group/cluster environment, but the dependence exists through its dependence on internal morphology.

preprint2015arXiv

Effects of hot halo gas on the star formation and mass transfer during distant galaxy-galaxy encounters

We use $N$-body/smoothed particle hydrodynamics simulations of encounters between an early-type galaxy (ETG) and a late-type galaxy (LTG) to study the effects of hot halo gas on the evolution for a case with the mass ratio of the ETG to LTG of 2:1 and the closest approach distance of $\sim$100 kpc. We find that the dynamics of the cold disk gas in the tidal bridge and the amount of the newly formed stars depend strongly on the existence of a gas halo. In the run of interacting galaxies not having a hot gas halo, the gas and stars accreted into the ETG do not include newly formed stars. However, in the run using the ETG with a gas halo and the LTG without a gas halo, a shock forms along the disk gas tidal bridge and induces star formation near the closest approach. The shock front is parallel to a channel along which the cold gas flows toward the center of the ETG. As a result, the ETG can accrete star-forming cold gas and newly born stars at and near its center. When both galaxies have hot gas halos, a shock is formed between the two gas halos somewhat before the closest approach. The shock hinders the growth of the cold gas bridge to the ETG and also ionizes it. Only some of the disk stars transfer through the stellar bridge. We conclude that the hot halo gas can give significant hydrodynamic effects during distant encounters.

preprint2015arXiv

The Eleventh and Twelfth Data Releases of the Sloan Digital Sky Survey: Final Data from SDSS-III

The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra.

preprint2015arXiv

The environment of barred galaxies revisited

We present a study of the environment of barred galaxies using galaxies drawn from the SDSS. We use several different statistics to quantify the environment: the projected two-point cross-correlation function, the background-subtracted number count of neighbor galaxies, the overdensity of the local environment, the membership of our galaxies to galaxy groups to segregate central and satellite systems, and for central galaxies we estimate the stellar to halo mass ratio (M$_{\mathrm{*}}/$M$_{\mathrm{h}}$) . When we split our sample into early- and late-type galaxies, we see a weak but significant trend for early-type galaxies with a bar to be more strongly clustered on scales from a few 100 kpc to 1 Mpc when compared to unbarred early-type galaxies. This indicates that the presence of a bar in early-type galaxies depends on the location within their host dark matter halos. This is confirmed by the group catalog in the sense that for early-types, the fraction of central galaxies is smaller if they have a bar. For late-type galaxies, we find fewer neighbors within $\sim50$ kpc around the barred galaxies when compared to unbarred galaxies from the control sample, suggesting that tidal forces from close companions suppress the formation/growth of bars. For central late-type galaxies, bars are more common on galaxies with high M$_{\mathrm{*}}/$M$_{\mathrm{h}}$ values, as expected from early theoretical works which showed that systems with massive dark matter halos are more stable against bar instabilities. Finally, we find no obvious correlation between overdensity and the bars in our sample, showing that galactic bars are not obviously linked to the large-scale structure of the universe.

preprint2015arXiv

Topology of neutral hydrogen distribution with the Square Kilometer Array

Morphology of the complex HI gas distribution can be quantified by statistics like the Minkowski functionals, and can provide a way to statistically study the large scale structure in the HI maps both at low redshifts, and during the epoch of reionization (EoR). At low redshifts, the 21cm emission traces the underlying matter distribution. Topology of the HI gas distribution, as measured by the genus, could be used as a "standard ruler". This enables the determination of distance-redshift relation and also the discrimination of various models of dark energy and of modified gravity. The topological analysis is also sensitive to certain primordial non-Gaussian features. Compared with two-point statistics, the topological statistics are more robust against the nonlinear gravitational evolution, bias, and redshift-space distortion. The HI intensity map observation naturally avoids the sparse sampling distortion, which is an important systematic in optical galaxy survey. The large cosmic volume accessible to SKA would provide unprecedented accuracy using such a measurement... [abridged]

preprint2015arXiv

Unusual A2142 supercluster with a collapsing core: distribution of light and mass

We study the distribution, masses, and dynamical properties of galaxy groups in the A2142 supercluster. We analyse the global luminosity density distribution in the supercluster and divide the supercluster into the high-density core and the low-density outskirts regions. We find galaxy groups and filaments in the regions of different global density, calculate their masses and mass-to-light ratios and analyse their dynamical state with several 1D and 3D statistics. We use the spherical collapse model to study the dynamical state of the supercluster. We show that in A2142 supercluster groups and clusters with at least ten member galaxies lie along an almost straight line forming a 50 Mpc/h long main body of the supercluster. The A2142 supercluster has a very high density core surrounded by lower-density outskirt regions. The total estimated mass of the supercluster is M_est = 6.2 10^{15}M_sun. More than a half of groups with at least ten member galaxies in the supercluster lie in the high-density core of the supercluster, centered at the rich X-ray cluster A2142. Most of the galaxy groups in the core region are multimodal. In the outskirts of the supercluster, the number of groups is larger than in the core, and groups are poorer. The orientation of the cluster A2142 axis follows the orientations of its X-ray substructures and radio halo, and is aligned along the supercluster axis. The high-density core of the supercluster with the global density D8 > 17 and perhaps with D8 > 13 may have reached the turnaround radius and started to collapse. A2142 supercluster with luminous, collapsing core and straight body is an unusual object among galaxy superclusters. In the course of the future evolution the supercluster may be split into several separate systems.

preprint2014arXiv

Exact Third-Order Density Perturbation and One-Loop Power Spectrum in General Dark Energy Models

Under the standard perturbation theory (SPT), we obtain the fully consistent third-order density fluctuation and kernels for the general dark energy models without using the Einstein-de Sitter (EdS) universe assumption for the first time. We also show that even though the temporal and spatial components of the SPT solutions can not be separable, one can find the exact solutions to any order in general dark energy models. With these exact solutions, we obtain the less than \% error correction of one-loop matter power spectrum compared to that obtained from the EdS assumption for $k = 0.1 {\rm h\, Mpc}^{-1}$ mode at $z = 0$ (1, 1.5). Thus, the EdS assumption works very well at this scale. However, if one considers the correction for $P_{13}$, the error is about 6 (9, 11) \% for the same mode at $z = 0$ (1, 1.5). One absorbs $P_{13}$ into the linear power spectrum in the renormalized perturbation theory (RPT) and thus one should use the exact solution instead of the approximation one. The error on the resummed propagator $N$ of RPT is about 14 (8, 6) \% at $z =0$ (1, 1.5) for $k = 0.4 {\rm h\, Mpc}^{-1}$. For $k = 1 {\rm h\, Mpc}^{-1}$, the error correction of the total matter power spectrum is about 3.6 (4.6, 4.5) \% at $z = 0$ (1, 1.5). Upcoming observation is required to archive the sub-percent accuracy to provide the strong constraint on the dark energy and this consistent solution is prerequisite for the model comparison.

preprint2014arXiv

Primordial non-Gaussian signatures in CMB polarization

We study the signatures of local type primordial non-Gaussianity, parametrized by $\fnl$, of scalar perturbations in CMB polarization using the probability distribution functions, Minkowski Functionals and Betti numbers. We show that the lowest order non-Gaussian deviation of the PDF of the total polarization intensity is at order $(\fnlσ)^2$. We calculate the non-Gaussian deviations of Minkowski Functionals and Betti numbers from simulated polarization maps. We find that $E$ mode polarization provides independent and equally strong constraint on $\fnl$ as temperature fluctuations. The non-Gaussian signal in the total polarization intensity, however, is much weaker and has a relatively large cosmic variance and hence may not be useful for detecting local type non-Gaussianity.

preprint2014arXiv

Tracing high redshift cosmic web with quasar systems

We trace the cosmic web at redshifts 1.0 <= z <= 1.8 using the quasar data from the SDSS DR7 QSO catalogue (Schneider et al. 2010). We apply a friend-of-friend (FoF) algorithm to the quasar and random catalogues to determine systems at a series of linking lengths, and analyse richness and sizes of these systems. At the linking lengths l <= 30 Mpc/h the number of quasar systems is larger than the number of systems detected in random catalogues, and systems themselves have smaller diameters than random systems. The diameters of quasar systems are comparable to the sizes of poor galaxy superclusters in the local Universe, the richest quasar systems have four members. The mean space density of quasar systems is close to the mean space density of local rich superclusters. At intermediate linking lengths (40 <= l <= 70 Mpc/h) the richness and length of quasar systems are similar to those derived from random catalogues. Quasar system diameters are similar to the sizes of rich superclusters and supercluster chains in the local Universe. At the linking length 70 Mpc/h the richest systems of quasars have diameters exceeding 500 Mpc/h. The percolating system which penetrate the whole sample volume appears in quasar sample at smaller linking length than in random samples (85 Mpc/h). Quasar luminosities in systems are not correlated with the system richness. Quasar system catalogues at our web pages http://www.aai.ee/~maret/QSOsystems.html serve as a database to search for superclusters of galaxies and to trace the cosmic web at high redshifts.

preprint2013arXiv

Alpha, Betti and the Megaparsec Universe: on the Topology of the Cosmic Web

We study the topology of the Megaparsec Cosmic Web in terms of the scale-dependent Betti numbers, which formalize the topological information content of the cosmic mass distribution. While the Betti numbers do not fully quantify topology, they extend the information beyond conventional cosmological studies of topology in terms of genus and Euler characteristic. The richer information content of Betti numbers goes along the availability of fast algorithms to compute them. For continuous density fields, we determine the scale-dependence of Betti numbers by invoking the cosmologically familiar filtration of sublevel or superlevel sets defined by density thresholds. For the discrete galaxy distribution, however, the analysis is based on the alpha shapes of the particles. These simplicial complexes constitute an ordered sequence of nested subsets of the Delaunay tessellation, a filtration defined by the scale parameter, $α$. As they are homotopy equivalent to the sublevel sets of the distance field, they are an excellent tool for assessing the topological structure of a discrete point distribution. In order to develop an intuitive understanding for the behavior of Betti numbers as a function of $α$, and their relation to the morphological patterns in the Cosmic Web, we first study them within the context of simple heuristic Voronoi clustering models. Subsequently, we address the topology of structures emerging in the standard LCDM scenario and in cosmological scenarios with alternative dark energy content. The evolution and scale-dependence of the Betti numbers is shown to reflect the hierarchical evolution of the Cosmic Web and yields a promising measure of cosmological parameters. We also discuss the expected Betti numbers as a function of the density threshold for superlevel sets of a Gaussian random field.

preprint2013arXiv

Betti numbers of Gaussian fields

We present the relation between the genus in cosmology and the Betti numbers for excursion sets of three- and two-dimensional smooth Gaussian random fields, and numerically investigate the Betti numbers as a function of threshold level. Betti numbers are topological invariants of figures that can be used to distinguish topological spaces. In the case of the excursion sets of a three-dimensional field there are three possibly non-zero Betti numbers; $β_0$ is the number of connected regions, $β_1$ is the number of circular holes, and $β_2$ is the number of three-dimensional voids. Their sum with alternating signs is the genus of the surface of excursion regions. It is found that each Betti number has a dominant contribution to the genus in a specific threshold range. $β_0$ dominates the high-threshold part of the genus curve measuring the abundance of high density regions (clusters). $β_1$ dominates the genus near the median thresholds which measures the topology of negatively curved iso-density surfaces, and $β_2$ corresponds to the low-threshold part measuring the void abundance. We average the Betti number curves (the Betti numbers as a function of the threshold level) over many realizations of Gaussian fields and find that both the amplitude and shape of the Betti number curves depend on the slope of the power spectrum $n$ in such a way that their shape becomes broader and their amplitude drops less steeply than the genus as $n$ decreases. This behaviour contrasts with the fact that the shape of the genus curve is fixed for all Gaussian fields regardless of the power spectrum. Even though the Gaussian Betti number curves should be calculated for each given power spectrum, we propose to use the Betti numbers for better specification of the topology of large scale structures in the universe.

preprint2013arXiv

Horizon Run 3: Topology as a Standard Ruler

We study the Physically Self Bound Cold Dark Matter Halo distribution which we associate with the massive galaxies within the Horizon Run 3 to estimate the accuracy in determination of the cosmological distance scale measured by the topology analysis. We apply the routine "Contour 3D" to 108 Mock Survey of $π$ steradians out to redshift z = 0.6, which effectively correspond to the SDSS-III BOSS survey, and compare the topology with that of a Gaussian Random Phase Field. We find that given three separate smoothing lengths $λ=$ 15, 21, and 34 $h^{-1}{\rm Mpc}$, the least $χ^2$ fit genus per unit volume g yields a 1.7 % fractional uncertainty in smoothing length and angular diameter distance to $z = 0.6$. This is an improvement upon former calibrations of and presents a competitive error estimate with next BAO scale techniques. We also present three dimensional graphics of the Horizon Run 3 spherical mock survey to show a wealth of large-scale structures of the universe that are predicted in surveys like BOSS.

preprint2013arXiv

On the galactic spin of barred disk galaxies

We present a study of the connection between the galactic spin parameter $λ_{d}$ and the bar fraction in a volume-limited sample of 10,674 disk galaxies drawn from the Sloan Digital Sky Survey Data Release 7. The galaxies in our sample are visually classified into galaxies hosting long or short bars, and non-barred galaxies. We find that the spin distributions of these three classes are statistically different, with galaxies hosting long bars with the lowest $λ_{d}$ values, followed by non-barred galaxies, while galaxies with short bars present typically high spin parameters. The bar fraction presents its maximum at low to intermediate $λ_{d}$ values for the case of long bars, while the maximum for short bars is at high $λ_{d}$. This bi-modality is in good agreement with previous studies finding longer bars hosted by luminous, massive, red galaxies with low content of cold gas, while short bars are found in low luminosity, low mass, blue galaxies, usually gas rich. In addition, the rise and fall of the bar fraction as a function of $λ_{d}$, within the long-bar sample, shown in our results, can be explained as a result of two competing factors: the self-gravity of the disk that enhances bar instabilities, and the support by random motions instead of ordered rotational motion, that prevents the formation/growth of bars.

preprint2012arXiv

A Second-order bias model for the Logarithmic Halo Mass Density

We present an analytic model for the local bias of dark matter halos in a LCDM universe. The model uses the halo mass density instead of the halo number density and is searched for various halo mass cuts, smoothing lengths, and redshift epoches. We find that, when the logarithmic density is used, the second-order polynomial can fit the numerical relation between the halo mass distribution and the underlying matter distribution extremely well. In this model the logarithm of the dark matter density is expanded in terms of log halo mass density to the second order. The model remains excellent for all halo mass cuts (from M_{cut}=3\times10^{11}$ to $3\times10^{12}h^{-1}M_{\odot}$), smoothing scales (from $R=5h^{-1}$Mpc to $50h^{-1}$Mpc), and redshift ranges (from z=0 to 1.0) considered in this study. The stochastic term in the relation is found not entirely random, but a part of the term can be determined by the magnitude of the shear tensor.

preprint2012arXiv

Hot and cold spots counts as probes of non-Gaussianity in the CMB

We introduce the numbers of hot and cold spots, $n_h$ and $n_c$, of excursion sets of the CMB temperature anisotropy maps as statistical observables that can discriminate different non-Gaussian models. We numerically compute them from simulations of non-Gaussian CMB temperature fluctuation maps. The first kind of non-Gaussian model we study is the local type primordial non-Gaussianity. The second kind of models have some specific form of the probability distribution function from which the temperature fluctuation value at each pixel is drawn, obtained using HEALPIX. We find the characteristic non-Gaussian deviation shapes of $n_h$ and $n_c$, which is distinct for each of the models under consideration. We further demonstrate that $n_h$ and $n_c$ carry additional information compared to the genus, which is just their linear combination, making them valuable additions to the Minkowski Functionals in constraining non-Gaussianity.

preprint2012arXiv

Measures of Galaxy Environment - I. What is "Environment"?

The influence of a galaxy's environment on its evolution has been studied and compared extensively in the literature, although differing techniques are often used to define environment. Most methods fall into two broad groups: those that use nearest neighbours to probe the underlying density field and those that use fixed apertures. The differences between the two inhibit a clean comparison between analyses and leave open the possibility that, even with the same data, different properties are actually being measured. In this work we apply twenty published environment definitions to a common mock galaxy catalogue constrained to look like the local Universe. We find that nearest neighbour-based measures best probe the internal densities of high-mass haloes, while at low masses the inter-halo separation dominates and acts to smooth out local density variations. The resulting correlation also shows that nearest neighbour galaxy environment is largely independent of dark matter halo mass. Conversely, aperture-based methods that probe super-halo scales accurately identify high-density regions corresponding to high mass haloes. Both methods show how galaxies in dense environments tend to be redder, with the exception of the largest apertures, but these are the strongest at recovering the background dark matter environment. We also warn against using photometric redshifts to define environment in all but the densest regions. When considering environment there are two regimes: the 'local environment' internal to a halo best measured with nearest neighbour and 'large-scale environment' external to a halo best measured with apertures. This leads to the conclusion that there is no universal environment measure and the most suitable method depends on the scale being probed.

preprint2012arXiv

Residual foreground contamination in the WMAP data and bias in non-Gaussianity estimation

We analyze whether there is any residual foreground contamination in the cleaned WMAP 7 years data for the differential assemblies, Q, V and W. We calculate the correlation between the foreground map, from which long wavelength correlations have been subtracted, and the foreground reduced map for each differential assembly after applying the Galaxy and point sources masks. We find positive correlations for all the differential assemblies, with high statistical significance. For Q and V, we find that a large fraction of the contamination comes from pixels where the foreground maps have positive values larger than three times the rms values. These findings imply the presence of residual contamination from Galactic emissions and unresolved point sources. We redo the analysis after masking the extended point sources cataloque of Scodeller et al. [7] and find a drop in the correlation and corresponding significance values. To quantify the effect of the residual contamination on the search for primordial non-Gaussianity in the CMB we add estimated contaminant fraction to simulated Gaussian CMB maps and calculate the characteristic non-Gaussian deviation shapes of Minkowski Functionals that arise due to the contamination. We find remarkable agreement of these deviation shapes with those measured from WMAP data, which imply that a major fraction of the observed non-Gaussian deviation comes from residual foreground contamination. We also compute non-Gaussian deviations of Minkowski Functionals after applying the point sources mask of Scodeller et al. and find a decrease in the overall amplitudes of the deviations which is consistent with a decrease in the level of contamination.

preprint2012arXiv

The Challenge of the Largest Structures in the Universe to Cosmology

Large galaxy redshift surveys have long been used to constrain cosmological models and structure formation scenarios. In particular, the largest structures discovered observationally are thought to carry critical information on the amplitude of large-scale density fluctuations or homogeneity of the universe, and have often challenged the standard cosmological framework. The Sloan Great Wall (SGW) recently found in the Sloan Digital Sky Survey (SDSS) region casts doubt on the concordance cosmological model with a cosmological constant (i.e. the flat LCDM model). Here we show that the existence of the SGW is perfectly consistent with the LCDM model, a result that only our very large cosmological N-body simulation (the Horizon Run 2, HR2) could supply. In addition, we report on the discovery of a void complex in the SDSS much larger than the SGW, and show that such size of the largest void is also predicted in the LCDM paradigm. Our results demonstrate that an initially homogeneous isotropic universe with primordial Gaussian random phase density fluctuations growing in accordance with the General Relativity, can explain the richness and size of the observed large-scale structures in the SDSS. Using the HR2 simulation we predict that a future galaxy redshift survey about four times deeper or with 3 magnitude fainter limit than the SDSS should reveal a largest structure of bright galaxies about twice as big as the SGW.

preprint2012arXiv

The New Horizon Run Cosmological N-Body Simulations

We present two large cosmological N-body simulations, called Horizon Run 2 (HR2) and Horizon Run 3 (HR3), made using 6000^3 = 216 billions and 7210^3 = 374 billion particles, spanning a volume of (7.200 Gpc/h)^3 and (10.815 Gpc/h)^3, respectively. These simulations improve on our previous Horizon Run 1 (HR1) up to a factor of 4.4 in volume, and range from 2600 to over 8800 times the volume of the Millennium Run. In addition, they achieve a considerably finer mass resolution, down to 1.25x10^11 M_sun/h, allowing to resolve galaxy-size halos with mean particle separations of 1.2 Mpc/h and 1.5 Mpc/h, respectively. We have measured the power spectrum, correlation function, mass function and basic halo properties with percent level accuracy, and verified that they correctly reproduce the LCDM theoretical expectations, in excellent agreement with linear perturbation theory. Our unprecedentedly large-volume N-body simulations can be used for a variety of studies in cosmology and astrophysics, ranging from large-scale structure topology, baryon acoustic oscillations, dark energy and the characterization of the expansion history of the Universe, till galaxy formation science - in connection with the new SDSS-III. To this end, we made a total of 35 all-sky mock surveys along the past light cone out to z=0.7 (8 from the HR2 and 27 from the HR3), to simulate the BOSS geometry. The simulations and mock surveys are already publicly available at http://astro.kias.re.kr/Horizon-Run23/.

preprint2012arXiv

The spin of late-type galaxies at redshifts z < 1.2

We study the evolution of the galactic spin using data of high redshift galaxies in the fields of the Great Observatories Origins Deep Survey (GOODS). Through simple dynamical considerations we estimate the spin for the disc galaxies in our sample and find that its distribution is consistent with that found for nearby galaxies. Defining a dimensionless angular momentum parameter for the disc component of the galaxies ($λ_{d}$), we do not find signs of evolution in the redshift range $0.4 \leq z \leq 1.2$. We find that the mass and environmental dependence of the spin of our high redshift galaxies are similar to that of low-$z$ galaxies; showing a strong dependence on mass, in the sense that low-mass systems present higher $λ_{d}$ values than high-mass galaxies, with no significant dependence on the environmental density. These results lead us to conclude that, although individual disc galaxies might occasionally suffer from strong evolution, they evolve in such a way that the overall spin distribution of the galactic population remains constant from $z\sim1$ to the present epoch.

preprint2011arXiv

Activity in galactic nuclei of cluster and field galaxies in the local universe

We study the environmental effects on the activity in galactic nuclei by comparing galaxies in clusters and in the field. Using a spectroscopic sample of galaxies in Abell clusters from the SDSS DR7, we investigate the dependence of nuclear activity on the physical parameters of clusters as well as the nearest neighbor galaxy. We also compare galaxy properties between AGN hosts and non-AGN galaxies. We find that the AGN fraction of early-type galaxies starts to decrease around one virial radius of clusters (r_{200,cl}) as decreasing clustercentric radius, while that of late types starts to decrease close to the cluster center (R~0.1-0.5r_{200,cl}). The AGN fractions of early-type cluster galaxies on average are found to be lower than those of early-type field galaxies by a factor ~3. However, the mean AGN fractions of late-type cluster galaxies are similar to those of late-type field galaxies. The AGN fraction of early-type BCGs lies between those of other early-type, cluster and field galaxies with similar luminosities. In the field, the AGN fraction strongly depends on the morphology of and the distance to the nearest neighbor galaxy. We find an anti-correlation between the AGN fraction and the velocity dispersion of clusters for all subsamples divided by morphology and luminosity of host galaxies. The AGN power indicated by L_[OIII]/M_BH is found to depend strongly on the mass of host galaxies rather than the clustercentric radius. The difference in physical parameters such as luminosity, (u-r) colors, SFRs, and (g-i) color gradients between AGN hosts and non-AGN galaxies is seen for both early and late types at all clustercentric radii, while the difference in structure parameters between the two is significant only for late types. These results support the idea that the nuclear activity is triggered through galaxy-galaxy interactions and mergers when gas supply for AGN is available.

preprint2011arXiv

Cosmic Voids in Sloan Digital Sky Survey Data Release 7

We study the distribution of cosmic voids and void galaxies using Sloan Digital Sky Survey Data Release 7 (SDSS DR7). Using the VoidFinder algorithm as described by Hoyle 2002, we identify 1054 statistically significant voids in the northern galactic hemisphere with radii > 10 h^{-1} Mpc. The filling factor of voids in the sample volume is 62%. The largest void is just over 30 h^{-1} Mpc in effective radius. The median effective radius is 17 h^{-1} Mpc. The voids are found to be significantly underdense, with density contrast δ< -0.85 at the edges of the voids. The radial density profiles of these voids are similar to predictions of dynamically distinct underdensities in gravitational theory. We find 8,046 galaxies brighter than M_r = -20.09 within the voids, accounting for 7% of the galaxies. We compare the results of VoidFinder on SDSS DR7 to mock catalogs generated from a SPH halo model simulation as well as other Λ-CDM simulations and find similar void fractions and void sizes in the data and simulations. This catalog is made publicly available at http://www.physics.drexel.edu/~pan/voidcatalog.html for download.

preprint2011arXiv

Dependence of Barred Galaxy Fraction on Galaxy Properties and Environment

We investigate the dependence of occurrence of bars in galaxies on galaxy properties and environment. We use a volume-limited sample of 33,391 galaxies brighter than $M_{r}=-19.5+5$log$h$ at $0.02\le z\le0.05489$, drawn from the SDSS DR 7. We classify the galaxies into early and late types, and identify bars by visual inspection. Among 10,674 late-type galaxies with axis ratio $b/a>0.60$, we find 3,240 barred galaxies ($f_{bar}=30.4%$) which divide into 2,542 strong bars ($f_{SB1}=23.8%$) and 698 weak bars ($f_{SB2}=6.5%$). We find that $f_{SB1}$ increases as $u-r$ color becomes redder, and that it has a maximum value at intermediate velocity dispersion ($σ\simeq$150 km s$^{-1}$). This trend suggests that strong bars are dominantly hosted by intermediate-mass systems. Weak bars prefer bluer galaxies with lower mass and lower concentration. In the case of strong bars, their dependence on the concentration index appears only for massive galaxies with $σ>150$ km s$^{-1}$. We also find that $f_{bar}$ does not directly depend on the large-scale background density when other physical parameters ($u-r$ color or $σ$) are fixed. We discover that $f_{SB1}$ decreases as the separation to the nearest neighbor galaxy becomes smaller than 0.1 times the virial radius of the neighbor regardless of neighbor's morphology. These results imply that strong bars are likely to be destroyed during strong tidal interactions, and that the mechanism for this phenomenon is gravitational and not hydrodynamical. The fraction of weak bars has no correlation with environmental parameters. We do not find any direct evidence for environmental stimulation of bar formation.

preprint2011arXiv

Probing Dark Energy with Alpha Shapes and Betti Numbers

We introduce a new descriptor of the weblike pattern in the distribution of galaxies and matter: the scale dependent Betti numbers which formalize the topological information content of the cosmic mass distribution. While the Betti numbers do not fully quantify topology, they extend the information beyond conventional cosmological studies of topology in terms of genus and Euler characteristic used in earlier analyses of cosmological models. The richer information content of Betti numbers goes along with the availability of fast algorithms to compute them. When measured as a function of scale they provide a "Betti signature" for a point distribution that is a sensitive yet robust discriminator of structure. The signature is highly effective in revealing differences in structure arising in different cosmological models, and is exploited towards distinguishing between different dark energy models and may likewise be used to trace primordial non-Gaussianities. In this study we demonstrate the potential of Betti numbers by studying their behaviour in simulations of cosmologies differing in the nature of their dark energy.

preprint2011arXiv

Quantifying galactic morphological transformations in the cluster environment

We study the effects of the cluster environment on galactic morphology by defining a dimensionless angular momentum parameter $λ_{d}$, to obtain a quantitative and objective measure of galaxy type. The use of this physical parameter allows us to take the study of morphological transformations in clusters beyond the measurements of merely qualitative parameters, e.g. S/E ratios, to a more physical footing. To this end, we employ an extensive Sloan Digital Sky Survey sample (Data Release 7), with galaxies associated with Abell galaxy clusters. The sample contains 121 relaxed Abell clusters and over 51,000 individual galaxies, which guarantees a thorough statistical coverage over a wide range of physical parameters. We find that the median $λ_{d}$ value tends to decrease as we approach the cluster center, with different dependences according to the mass of the galaxies and the hosting cluster; low and intermediate mass galaxies showing a strong dependence, while massive galaxies seems to show, at all radii, low $λ_{d}$ values. By analysing trends in $λ_{d}$ as functions of the nearest neighbour environment, clustercentric radius and velocity dispersion of clusters, we can identify clearly the leading physical processes at work. We find that in massive clusters ($σ>700$ km/s), the interaction with the cluster central region dominates, whilst in smaller clusters galaxy-galaxy interactions are chiefly responsible for driving galactic morphological transformations.

preprint2011arXiv

The galactic spin of AGN galaxies

Using an extensive sample of galaxies selected from the Sloan Digital Sky Survey Data Release 5, we compare the angular momentum distribution of active galactic nucleus (AGN) with non-AGN hosting late-type galaxies. To this end we characterize galactic spin through the dimensionless angular momentum parameter $λ$, which we estimate approximately through simple dynamical considerations. Using a volume limited sample, we find a considerable difference when comparing the empirical distributions of $λ$ for AGNs and non-AGN galaxies, the AGNs showing typically low $λ$ values and associated dispersions, while non-AGNs present higher $λ$ values and a broader distribution. A more striking difference is found when looking at $λ$ distributions in thin $M_{r}$ cuts, while the spin of non-AGN galaxies presents an anti-correlation with $M_{r}$, with bright (massive) galaxies having low spins, AGN host galaxies present uniform values of $λ$ at all magnitudes, a behavior probably imposed by the fact that most late-type AGN galaxies present a narrow range in color, with a typical constant $λ$ value. We also find that the fraction of AGN hosting galaxies in our sample strongly depends on galactic spin, increasing dramatically for decreasing $λ$. For AGN host galaxies we compute the mass of their supermassive black holes and find that this value tends to be higher for low spin galaxies, even at fixed luminosity, a result that could account, to a certain extent, for the spread on the luminosity black-hole mass relation.

preprint2011arXiv

Topology of large scale structure as test of modified gravity

The genus of the iso-density contours is a robust measure of the topology of large scale structure, and it is relatively insensitive to nonlinear gravitational evolution, galaxy bias and redshift-space distortion. We show that the growth of density fluctuations is scale-dependent even in the linear regime in some modified gravity theories, which opens a new possibility of testing the theories observationally. We propose to use the genus of the iso-density contours, an intrinsic measure of the topology of large scale structure, as a statistic to be used in such tests. In Einstein's general theory of relativity, density fluctuations are growing at the same rate on all scales in the linear regime, and the genus per comoving volume is almost conserved as structures are growing homologously, so we expect that the genus-smoothing scale relation is basically time-independent. However, in some modified gravity models where structures grow with different rates on different scales, the genus-smoothing scale relation should change over time. This can be used to test the gravity models with large scale structure observations. We studied the case of the f(R) theory, DGP braneworld theory as well as the parameterized post-Friedmann (PPF) models. We also forecast how the modified gravity models can be constrained with optical/IR or redshifted 21cm radio surveys in the near future.

preprint2010arXiv

Distribution of satellite galaxies in high redshift groups

We use galaxy groups at redshifts between 0.4 and 1.0 selected from the Great Observatories Origins Deep Survey (GOODS) to study the color-morphological properties of satellite galaxies, and investigate possible alignment between the distribution of the satellites and the orientation of their central galaxy. We confirm the bimodal color and morphological type distribution for satellite galaxies at this redshift range: the red and blue classes corresponds to the early and late morphological types respectively, and the early-type satellites are on average brighter than the late-type ones. Furthermore, there is a {\it morphological conformity} between the central and satellite galaxies: the fraction of early-type satellites in groups with an early-type central is higher than those with a late-type central galaxy. This effect is stronger at smaller separations from the central galaxy. We find a marginally significant signal of alignment between the major axis of the early-type central galaxy and its satellite system, while for the late-type centrals no significant alignment signal is found. We discuss the alignment signal in the context of shape evolution of groups.

preprint2010arXiv

Environmental Dependence of Local Luminous Infrared Galaxies

We study the environmental dependence of local luminous infrared galaxies (LIRGs) and ultraluminous infrared galaxies (ULIRGs) found in the Sloan Digital Sky Survey (SDSS) data. The LIRG and ULIRG samples are constructed by cross-correlating spectroscopic catalogs of galaxies of the SDSS Data Release 7 and the Infrared Astronomical Satellite Faint Source Catalog. We examine the effects of the large-scale background density (Sigma_5), galaxy clusters, and the nearest neighbor galaxy on the properties of infrared galaxies (IRGs). We find that the fraction of LIRGs plus ULIRGs among IRGs (f_(U)LIRGs) and the infrared luminosity (L_IR) of IRGs strongly depend on the morphology of and the distance to the nearest neighbor galaxy: the probability for an IRG to be a (U)LIRG (f_(U)LIRGs) and its L_IR both increase as it approaches a late-type galaxy, but decrease as it approaches an early-type galaxy (within half the virial radius of its neighbor). We find no dependence of f_(U)LIRGs on the background density (surface galaxy number density) at fixed stellar mass of galaxies. The dependence of f_(U)LIRGs on the distance to galaxy clusters is also found to be very weak, but in highest-density regions such as the center of galaxy clusters, few (U)LIRGs are found. These environmental dependence of LIRGs and ULIRGs and the evolution of star formation rate (SFR)-environment relation from high redshifts to low redshifts seem to support the idea that galaxy-galaxy interactions/merging play a critical role in triggering the star formation activity of LIRGs and ULIRGs.

preprint2010arXiv

Galaxy Clustering Topology in the Sloan Digital Sky Survey Main Galaxy Sample: a Test for Galaxy Formation Models

We measure the topology of the main galaxy distribution using the Seventh Data Release of the Sloan Digital Sky Survey, examining the dependence of galaxy clustering topology on galaxy properties. The observational results are used to test galaxy formation models. A volume-limited sample defined by $M_r<-20.19$ enables us to measure the genus curve with amplitude of $G=378$ at $6h^{-1}$Mpc smoothing scale, with 4.8\% uncertainty including all systematics and cosmic variance. The clustering topology over the smoothing length interval from 6 to $10 h^{-1}$Mpc reveals a mild scale-dependence for the shift ($Δν$) and void abundance ($A_V$) parameters of the genus curve. We find substantial bias in the topology of galaxy clustering with respect to the predicted topology of the matter distribution, which varies with luminosity, morphology, color, and the smoothing scale of the density field. The distribution of relatively brighter galaxies shows a greater prevalence of isolated clusters and more percolated voids. Even though early (late)-type galaxies show topology similar to that of red (blue) galaxies, the morphology dependence of topology is not identical to the color dependence. In particular, the void abundance parameter $A_V$ depends on morphology more strongly than on color. We test five galaxy assignment schemes applied to cosmological N-body simulations of a $Λ$CDM universe to generate mock galaxies: the Halo-Galaxy one-to-one Correspondence model, the Halo Occupation Distribution model, and three implementations of Semi-Analytic Models (SAMs). None of the models reproduces all aspects of the observed clustering topology; the deviations vary from one model to another but include statistically significant discrepancies in the abundance of isolated voids or isolated clusters and the amplitude and overall shift of the genus curve. (Abridged)

preprint2010arXiv

Gravitational Potential Environment of Galaxies: I. Simulation

We extend the concept of galaxy environment from the local galaxy number density to the gravitational potential and its functions like the shear tensor. For this purpose we examine whether or not one can make an accurate estimation of the gravitational potential from an observational sample which is finite in volume, biased due to galaxy biasing, and subject to redshift space distortion. Dark halos in a $Λ$CDM simulation are used in this test. We find that one needs to stay away from the sample boundaries by more than 30$h^{-1}$Mpc to reduce the error within 20% of the root mean square values of the potential or the shear tensor. The error due to the galaxy biasing can be significantly reduced by using the galaxy mass density field instead of the galaxy number density field. The error caused by the redshift space distortion can be effectively removed by correcting galaxy positions for the peculiar velocity effects. We inspect the dependence of dark matter halo properties on four environmental parameters; local density, gravitational potential, and the ellipticity and prolateness of the shear tensor. We find the local density has the strongest correlation with halo properties. This is evidence that the internal physical properties of dark halos are mainly controlled by small-scale physics. In high density regions dark halos are on average more massive and spherical, and have higher spin parameter and velocity dispersion. In high density regions dark halos are on average more massive and spherical, and have higher spin parameter and velocity dispersion. We also study the relation between the environmental parameters and the subtypes of dark halos. The spin parameter of satellite halos depends only weakly on the local density for all mass ranges studied while that of isloated or central halos depends more sensitively on the local density. (abridged)

preprint2010arXiv

Statistics of the excursion sets in models with local primordial non-Gaussianity

We use the statistics of regions above or below a temperature threshold (excursion sets) to study the cosmic microwave background (CMB) anisotropy in models with primordial non-Gaussianity of the local type. By computing the full-sky spatial distribution and clustering of pixels above/below threshold from a large set of simulated maps with different levels of non-Gaussianity, we find that a positive value of the dimensionless non-linearity parameter f_NL enhances the number density of the cold CMB excursion sets along with their clustering strength, and reduces that of the hot ones. We quantify the robustness of this effect, which may be important to discriminate between the simpler Gaussian hypothesis and non-Gaussian scenarios, arising either from non-standard inflation or alternative early-universe models. The clustering of hot and cold pixels exhibits distinct non-Gaussian signatures, particularly at angular scales of about 75 arcmin (i.e. around the Doppler peak), which increase linearly with f_NL. Moreover, the clustering changes strongly as a function of the smoothing angle. We propose several statistical tests to maximize the detection of a local primordial non-Gaussian signal, and provide some theoretical insights within this framework, including an optimal selection of the threshold level. We also describe a procedure which aims at minimizing the cosmic variance effect, the main limit within this statistical framework.

preprint2010arXiv

The Properties of Type Ia Supernova Host Galaxies from the Sloan Digital Sky Survey

We investigate the properties and environments of Type Ia Supernova (SN Ia) host galaxies in the Stripe 82 of the Sloan Digital Sky Survey-II Supernova Survey centered on the celestial equator. Host galaxies are defined as the galaxy nearest to the supernova (SN) in terms of angular distance whose velocity difference from the SN is less than 1000 km s^{-1}. Eighty seven SN Ia host galaxies are selected from the SDSS Main galaxy sample with the apparent r-band magnitude m_r < 17.77, and compared with the SDSS Main galaxies. The SN Ia rates for early and late-type galaxies are 0.81 +- 0.19 SN (100yr)^{-1} and 0.99 +- 0.21 SN (100yr)^{-1}, respectively. We find that the host galaxies have a color distribution consistent with that of the Main galaxies, regardless of their morphology. However, host galaxies are on average brighter than the Main galaxies by ~ 0.3 mag over the range of -18.3 > M_r > -21.3. But the brighter ends of their luminosity distributions are similar. The distribution of the distance to the nearest neighbor galaxy shows that SNe Ia are more likely to occur in isolated galaxies without close neighbors. We also find that the SN Ia host galaxies are preferentially located in a region close to massive galaxy clusters compared to the Main galaxies.

preprint2010arXiv

The UNAM-KIAS Catalog of Isolated Galaxies

A new catalog of isolated galaxies from The Sloan Digital Sky Survey (DR5) is presented. 1520 isolated galaxies were found in 1.4 steradians of sky. The selection criteria in this so called UNAM-KIAS catalog was implemented from a variation on the criteria developed by Karachentseva 1973 including full redshift information. Through an image processing pipeline that takes advantage from the high resolution (~ 0.4 ''/pix) and high dynamic range of the SDSS images, a uniform g band morphological classification for all these galaxies is presented. We identify 80% (SaSm) spirals (50% later than Sbc types) on one hand, and a scarce population of early-type E(6.5%) and S0(8%) galaxies amounting to 14.5% on the other hand. This magnitude-limited catalog is ~ 80% complete at 16.5, 15.6, 15.0, 14.6 and 14.4 magnitudes in the ugriz bands respectively. Some representative physical properties including SDSS magnitudes and color distributions, color-color diagrams, absolute magnitude-color, and concentration-color diagrams as a function of morphological type are presented. The UNAM-KIAS Morphological Atlas is also released along with this paper. For each galaxy of type later than Sa, a mosaic is presented that includes: (1) a g-band logarithmic image, (2) a g band filtered-enhanced image where a Gaussian kernel of various sizes was applied and (3) an RGB color image from the SDSS database. For E/S0/Sa galaxies, in addition to the images in (1), (2) and (3), plots of r band surface brightness and geometric profiles (ellipticity, Position Angle PA and A4/B4 coefficients of the Fourier series expansions of deviations of a pure ellipse) are provided...

preprint2009arXiv

Clues on the origin of galactic angular momentum from looking at galaxy pairs

We search for correlations between the spin in pairs of spiral galaxies, to study if the angular momentum gain for each galaxy was the result of tidal torques imprint by the same tidal field. To perform our study we made use of a sample of galaxy pairs identified using the Sloan Digital Sky Survey. We find a weak, but statistically significant correlation between the spin magnitude of neighbouring galaxies, but no clear alignment between their orientation. We show that events such as interactions with close neighbours play an important role in the value of the spin for the final configuration, as we find these interactions tend to reduce the value of the $λ$ spin parameter of late-type galaxies considerably, with dependence on the morphology of the neighbour. This implies that the original tidal field for each pair could have been similar, but the redistribution of angular momentum at later stages of evolution is important.

preprint2009arXiv

Interactions of Galaxies in the Galaxy Cluster Environment

We study the dependence of galaxy properties on the clustercentric radius and the environment attributed to the nearest neighbor galaxy using the SDSS galaxies associated with the Abell galaxy clusters. We find that there exists a characteristic scale where the properties of galaxies suddenly start to depend on the clustercentric radius at fixed neighbor environment. The characteristic scale is $1\sim 3$ times the cluster virial radius depending on galaxy luminosity. Existence of the characteristic scale means that the local galaxy number density is not directly responsible for the morphology-density relation in clusters because the local density varies smoothly with the clustercentric radius and has no discontinuity in general. What is really working in clusters is the morphology-clustercentric radius-neighbor environment relation, where the neighbor environment means both neighbor morphology and the local mass density attributed to the neighbor. The morphology-density relation appears working only because of the statistical correlation between the nearest neighbor distance and the local galaxy number density. We find strong evidence that the hydrodynamic interactions with nearby early-type galaxies is the main drive to quenching star formation activity of late-type galaxies in clusters. The hot cluster gas seems to play at most a minor role down to one tenth of the cluster virial radius. We also find that the viable mechanisms which can account for the clustercentric radius dependence of the structural and internal kinematics parameters are harassment and interaction of galaxies with the cluster potential. The morphology transformation of the late-type galaxies in clusters seems to have taken place through both galaxy-galaxy hydrodynamic interactions and galaxy-cluster/galaxy-galaxy gravitational interactions.

preprint2009arXiv

Large-Scale Structure of the Universe as a Cosmic Standard Ruler

We propose to use the large-scale structure of the universe as a cosmic standard ruler, based on the fact that the pattern of galaxy distribution should be maintained in the course of time on large scales. By examining the scale-dependence of the pattern in different redshift intervals it is possible to reconstruct the expansion history of the universe, and thus to measure the cosmological parameters governing the expansion of the universe. The features in the galaxy distribution that can be used as standard rulers include the topology of large-scale structure and the overall shapes of galaxy power spectrum and correlation function. The genus, being an intrinsic topology measure, is resistant against the non-linear gravitational evolution, galaxy biasing, and redshift-space distortion effects, and thus is ideal for quantifying the primordial topology of the large-scale structure. The expansion history of the universe can be constrained by comparing among the genus measured at different redshifts. In the case of initially Gaussian fluctuations the genus accurately recovers the slope of the primordial power spectrum near the smoothing scale, and the expansion history can be constrained by comparing between the predicted and measured genus.

preprint2009arXiv

Origin and evolution of galactic spin from looking at galaxy pairs

To study if the angular momentum gain for each member of a galaxy pair was the result of tidal torques imprinted by the same tidal field, we search for correlations between the spin in pairs of spiral galaxies identified using the Sloan Digital Sky Survey. We find a weak, but statistically significant correlation between the spin magnitude of neighbouring galaxies. We show that events such as close interactions with neighbours play an important role in the value of the spin for the final configuration, as we find these interactions tend to reduce the value of the $λ$ spin parameter of late-type galaxies considerably. This implies that the original tidal field for each pair could have been similar, but the redistribution of angular momentum at later stages of evolution is important

preprint2009arXiv

Reconstructing galaxy fundamental distributions and scaling relations from photometric redshift surveys. Applications to the SDSS early-type sample

Noisy distance estimates associated with photometric rather than spectroscopic redshifts lead to a mis-estimate of the luminosities, and produce a correlated mis-estimate of the sizes. We consider a sample of early-type galaxies from the SDSS DR6 for which both spectroscopic and photometric information is available, and apply the generalization of the V_max method to correct for these biases. We show that our technique recovers the true redshift, magnitude and size distributions, as well as the true size-luminosity relation. We find that using only 10% of the spectroscopic information randomly spaced in our catalog is sufficient for the reconstructions to be accurate within about 3%, when the photometric redshift error is dz = 0.038. We then address the problem of extending our method to deep redshift catalogs, where only photometric information is available. In addition to the specific applications outlined here, our technique impacts a broader range of studies, when at least one distance-dependent quantity is involved. It is particularly relevant for the next generation of surveys, some of which will only have photometric information.

preprint2009arXiv

Statistical nature of non-Gaussianity from cubic order primordial perturbations: CMB map simulations and genus statistic

We simulate CMB maps including non-Gaussianity arising from cubic order perturbations of the primordial gravitational potential, characterized by the non-linearity parameter $g_{NL}$. The maps are used to study the characteristic nature of the resulting non-Gaussian temperature fluctuations. We measure the genus and investigate how it deviates from Gaussian shape as a function of $g_{NL}$ and smoothing scale. We find that the deviation of the non-Gaussian genus curve from the Gaussian one has an antisymmetric, sine function like shape, implying more hot and more cold spots for $g_{NL}>0$ and less of both for $g_{NL}<0$. The deviation increases linearly with $g_{NL}$ and also exhibits mild increase as the smoothing scale increases. We further study other statistics derived from the genus, namely, the number of hot spots, the number of cold spots, combined number of hot and cold spots and the slope of the genus curve at mean temperature fluctuation. We find that these observables carry signatures of $g_{NL}$ that are clearly distinct from the quadratic order perturbations, encoded in the parameter $f_{NL}$. Hence they can be very useful tools for distinguishing not only between non-Gaussian temperature fluctuations and Gaussian ones but also between $g_{NL}$ and $f_{NL}$ type non-Gaussianities.

preprint2009arXiv

The Horizon Run N-body Simulation: Baryon Acoustic Oscillations and Topology of Large Scale Structure of the Universe

In support of the new Sloan III survey, which will measure the baryon oscillation scale using the luminous red galaxies (LRGs), we have run the largest N-body simulation to date using $4120^3 = 69.9$ billion particles, and covering a volume of $(6.592 h^{-1} {\rm Gpc})^3$. This is over 2000 times the volume of the Millennium Run, and corner-to-corner stretches all the way to the horizon of the visible universe. LRG galaxies are selected by finding the most massive gravitationally bound, cold dark matter subhalos, not subject to tidal disruption, a technique that correctly reproduces the 3D topology of the LRG galaxies in the Sloan Survey.We have measured the covariance function, power spectrum, and the 3D topology of the LRG galaxy distribution in our simulation and made 32 mock surveys along the past lightcone to simulate the Sloan III survey. Our large N-body simulation is used to accurately measure the non-linear systematic effects such as gravitational evolution, redshift space distortion, past light cone space gradient, and galaxy biasing, and to calibrate the baryon oscillation scale and the genus topology. For example, we predict from our mock surveys that the baryon acoustic oscillation peak scale can be measured with the cosmic variance-dominated uncertainty of about 5% when the SDSS-III sample is divided into three equal volume shells, or about 2.6% when a thicker shell with $0.4<z<0.6$ is used. We find that one needs to correct the scale for the systematic effects amounting up to 5.2% to use it to constrain the linear theories. And the uncertainty in the amplitude of the genus curve is expected to be about 1% at 15 $h^{-1}$Mpc scale. We are making the simulation and mock surveys publicly available.

preprint2009arXiv

The nature of the SDSS galaxies in various classes based on morphology, colour and spectral features - II. Multi-wavelength properties

We present a multi-wavelength study of the nature of the SDSS galaxies divided into fine classes based on their morphology, colour and spectral features. The SDSS galaxies are classified into early-type and late-type; red and blue; passive, HII, Seyfert and LINER, which returns a total of 16 fine classes of galaxies. The properties of galaxies in each fine class are investigated from radio to X-ray, using 2MASS, IRAS, FIRST, NVSS, GALEX and ROSAT data. The UV - optical - NIR colours of blue early-type galaxies (BEGs) seem to result from the combination of old stellar population and recent star formation (SF). Non-passive red early-type galaxies (REGs) have larger metallicity and younger age than passive REGs, which implies that non-passive REGs have suffered recent SF adding young and metal-rich stars to them. The radio detection fraction of REGs strongly depends on their optical absolute magnitudes, while that of most late-type galaxies does not, implying the difference in their radio sources: AGN and SF. The UV - optical colours and the radio detection fraction of passive RLGs show that they have properties similar to REGs rather than non-passive RLGs. Dust extinction may not be a dominant factor making RLGs red, because RLGs are detected in the mid- and far-infrared bands less efficiently than blue late-type galaxies (BLGs). The passive BLGs have very blue UV - optical - NIR colours, implying either recent SF quenching or current SF in their outskirts. Including star formation rate, other multi-wavelength properties in each fine class are investigated, and their implication on the identity of each fine class is discussed (abridged).

preprint2009arXiv

The nature of the Sloan Digital Sky Survey galaxies in various classes based on morphology, colour and spectral features - III. Environments

We present a study on the environments of the SDSS galaxies divided into fine classes based on their morphology, colour and spectral features. The SDSS galaxies are classified into early-type and late-type; red and blue; passive, HII, Seyfert and LINER, which returns a total of 16 fine classes of galaxies. We estimate the local number density, target-excluded local luminosity density, local colour, close pair fraction and the luminosity and colour of the brightest neighbour, which are compared between the fine classes comprehensively. The morphology-colour class of galaxies strongly depends on the local density, with the approximate order of high-density preference: red early-type galaxies (REGs) -- red late-type galaxies (RLGs) -- blue early-type galaxies (BEGs) -- blue late-type galaxies (BLGs). We find that high-density environments (like cluster environments) seem to suppress AGN activity. The pair fraction of HII REGs does not show statistically significant difference from that of passive REGs, while the pair fraction of HII BLGs is smaller than that of non-HII BLGs. HII BLGs show obvious double (red + blue) peaks in the distribution of the brightest neighbour colour, while red galaxies show a single red peak. The brightest neighbours of Seyfert BLGs tend to be blue, while those of LINER BLGs tend to be red, which implies that the difference between Seyfert and LINER may be related to the pair interaction. Other various environments of the fine classes are investigated, and their implication on galaxy evolution is discussed.

preprint2008arXiv

Combined Effects of Galaxy Interactions and Large-Scale Environment on Galaxy Properties

We inspect the coupled dependence of physical parameters of the SDSS galaxies on the small-scale (distance to and morphology of the nearest neighbor galaxy) and the large-scale (background density smoothed over 20 nearby galaxies) environments. The impacts of interaction on galaxy properties are detected at least out to the neighbor separation corresponding to the virial radius of galaxies, which is typically between 200 and 400 h^{-1}kpc for the galaxies in our sample. To detect these long-range interaction effects it is crucial to divide galaxy interactions into four cases dividing the morphology of target and neighbor galaxies into early and late types. We show that there are two characteristic neighbor-separation scales where the galaxy interactions cause abrupt changes in the properties of galaxies. The first scale is the virial radius of the nearest neighbor galaxy r_{vir,nei}. Many physical parameters start to deviate from those of extremely isolated galaxies at the projected neighbor separation r_p of about r_{vir,nei}. The second scale is at r_p \approx 0.05 r_{vir,nei} = 10--20 h^{-1} kpc, and is the scale at which the galaxies in pairs start to merge. We find that late-type neighbors enhance the star formation activity of galaxies while early-type neighbors reduce it, and that these effects occur within r_{vir,nei}. The hot halo gas and cold disk gas must be participating in the interactions at separations less than the virial radius of the galaxy plus dark halo system. Our results also show that the role of the large-scale density in determining galaxy properties is minimal once luminosity and morphology are fixed. We propose that the weak residual dependence of galaxy properties on the large-scale density is due to the dependence of the halo gas property on the large-scale density.

preprint2008arXiv

Internal Extinction in the SDSS Late-Type Galaxies

We study internal extinction of late-type galaxies in the Sloan Digital Sky Survey. We find that the degree of internal extinction depends on both the concentration index c and K_s-band absolute magnitude M_K. We give simple fitting functions for internal extinction. In particular, we present analytic formulae giving the extinction-corrected magnitudes from the observed optical parameters. For example, the extinction-corrected r-band absolute magnitude can be obtained by M_{r,0}=-20.77 +(-1+\sqrt{1+4Δ(M_{r,obs}+20.77+4.93Δ)})/2Δ, where Δ= 0.236\{1.35(c-2.48)^2-1.14\}\log(a/b), c=R_{90}/R_{50} is the the concentration index, and a/b is the isophotal axis ratio of the 25 mag/arcsec^2 isophote in the i-band. The 1σerror in M_{r,0} is 0.21{\rm log}(a/b). The late-type galaxies with very different inclinations are found to trace almost the same sequence in the (u-r)-$M_r$ diagram when our prescriptions for extinction correction are applied. We also find that (u-r) color can be a third independent parameter that determines the degree of internal extinction.

preprint2007arXiv

Internal and Collective Properties of Galaxies in the Sloan Digital Sky Survey

We examine volume-limited samples from the SDSS galaxies to look for relations among internal and collective physical parameters of galaxies as faint as M_r=- 17.5, which include morphology, luminosity, color, color gradient, concentration, size, velocity dispersion, equivalent width (EW) of H_alpha line,axis ratio, the luminosity and velocity dispersion functions. At fixed morphology and luminosity, we find that bright (M_r<-20) early-types show very small dispersions in color, color gradient, concentration, size, and velocity dispersion. These dispersions increase at fainter magnitudes, where the fraction of blue early-types increases. Late-types show wider dispersions in all physical parameters compared to early types at the same luminosity. Concentration indices of early-types are well-correlated with velocity dispersion, but are insensitive to luminosity and color for bright galaxies. The slope of the Faber-Jackson relation continuously changes from 4.6 +- 0.4 to 2.7+- 0.2 when luminosity changes from M_r = -22 to -20. The size of early- types is well-correlated with stellar velocity dispersion (for >100 km/s). We find that passive spirals are well-separated from star-forming late-types at EW (H_alpha) of about 4. An interesting finding is that many physical parameters of galaxies manifest different behaviors across the absolute magnitude of about M_* +- 1. The morphology fraction as a function of luminosity depends less sensitively on large scale structure than the luminosity function (LF) does, and thus seems to be more universal. The effects of internal extinction in late-types on the completeness of volume limited samples and on the LF and morphology fraction are found to be very important. (abridged)

preprint2005arXiv

Fourier Phase Analysis of SDSS Galaxies

We present a first analysis of the clustering of SDSS galaxies using the distribution function of the sum of Fourier phases. This statistic was recently proposed by one of authors as a new method to probe phase correlations of cosmological density fields. Since the Fourier phases are statistically independent of Fourier amplitudes, the phase statistic plays a complementary role to the conventional two-point statistics of galaxy clustering. In particular, we focus on the distribution functions of phase sum over three closed wavevectors as a function of the configurations of triangle wavevectors. We find that the observed distribution functions of phase sum are in good agreement with the lowest-order approximation from perturbation theory. For direct comparison with observations, we construct mock catalogs from N-body simulations taking account of several observational effects such as the survey geometry, the redshift distortion, and the discreteness due to the limited number of data. Indeed the observed phase correlations for the galaxies in the range of absolute magnitude -22<Mr<-18 agree well with those for Lambda CDM predictions with sigma8=0.9 evolved from Gaussian initial condition. This agreement implies that the galaxy biasing is approximately linear in redshift space. Instead, assuming that the galaxy biasing is described by a quadratic deterministic function at scale larger than 30Mpc/h, we can constrain the ratio of the quadratic biasing parameter b2 to the linear biasing parameter b1. We find that resulting b2/b1 is mainly dependent on sigma8 and that it is well fitted by b2/b1=0.54(+-0.06)-0.62(+-0.08)sigma8 in CDM models. Indeed, b2/b1 is nearly zero when sigma8=0.9.

preprint2003arXiv

GOTPM: A Parallel Hybrid Particle-Mesh Treecode

We describe a parallel, cosmological N-body code based on a hybrid scheme using the particle-mesh (PM) and Barnes-Hut (BH) oct-tree algorithm. We call the algorithm GOTPM for Grid-of-Oct-Trees-Particle-Mesh. The code is parallelized using the Message Passing Interface (MPI) library and is optimized to run on Beowulf clusters as well as symmetric multi-processors. The gravitational potential is determined on a mesh using a standard PM method with particle forces determined through interpolation. The softened PM force is corrected for short range interactions using a grid of localized BH trees throughout the entire simulation volume in a completely analogous way to P$^3$M methods. This method makes no assumptions about the local density for short range force corrections and so is consistent with the results of the P$^3$M method in the limit that the treecode opening angle parameter, $θ\to 0$. (abridged)

preprint2002arXiv

Correlation Functions of the Abell, APM, and X-ray Clusters of Galaxies

We have measured the correlation functions of the optically selected clusters of galaxies in the Abell and the APM catalogs, and of the X-ray clusters in the X-ray-Brightest Abell-type Clusters of galaxies (XBACs) catalog and the Brightest Clusters Sample (BCS). The same analysis method and the same method of characterizing the resulting correlation functions are applied to all observational samples. We have found that the amplitude of the correlation function of the APM clusters is much higher than what has been previously claimed, in particular for richer subsamples. The correlation length of the APM clusters with the richness R >= 70 (as defined by the APM team) is found to be r_0 = 25.4(+3.1,-3.0) h^(-1)Mpc. The amplitude of correlation function is about 2.4 times higher than that of Croft et al. (1997). The correlation lengths of the Abell clusters with the richness class RC >= 0 and 1 are measured to be r_0 = 17.4(+1.2,-1.1) and 21.0(+2.8,-2.8) h^(-1)Mpc, respectively, which is consistent with our results for the APM sample at the similar level of richness. The richness dependence of cluster correlations is found to be r_0 = 0.40*d_c + 3.2 where d_c is the mean intercluster separation. This is identical in slope with the Bahcall & West (1992)'s estimate, but is inconsistent with the weak dependence of Croft et al. (1997). The X-ray bright Abell clusters in the XBACs catalog and the X-ray selected clusters in the BCS catalog show strong clustering. The correlation length of the XBACs clusters with L_x >= 0.65*10^44 h^(-2)erg/s is 30.3(+8.2,-6.5) h^(-1)Mpc, and that of the BCS clusters with L_x >= 0.70*10^44 h^(-2)erg/s is 30.2(+9.8,-8.9) h^(-1)Mpc. The clustering strength of the X-ray clusters is much weaker than what is expected from the optical clusters.

preprint2001arXiv

Effects of Foreground Contamination on the Cosmic Microwave Background Anisotropy Measured by MAP

We study the effects of diffuse Galactic, far-infrared extragalactic source, and radio point source emission on the cosmic microwave background (CMB) anisotropy data anticipated from the MAP experiment. We focus on the correlation function and genus statistics measured from mock MAP foreground-contaminated CMB anisotropy maps generated in a spatially-flat cosmological constant dominated cosmological model. Analyses of the simulated MAP data at 90 GHz (0.3 deg FWHM resolution smoothed) show that foreground effects on the correlation function are small compared with cosmic variance. However, the Galactic emission, even just from the region with |b| > 20 deg, significantly affects the topology of CMB anisotropy, causing a negative genus shift non-Gaussianity signal. Given the expected level of cosmic variance, this effect can be effectively reduced by subtracting existing Galactic foreground emission models from the observed data. IRAS and DIRBE far-infrared extragalactic sources have little effect on the CMB anisotropy. Radio point sources raise the amplitude of the correlation function considerably on scales below 0.5 deg. Removal of bright radio sources above a 5 σdetection limit effectively eliminates this effect. Radio sources also result in a positive genus curve asymmetry (significant at 2 σ) on 0.5 deg scales. Accurate radio point source data is essential for an unambiguous detection of CMB anisotropy non-Gaussianity on these scales. Non-Gaussianity of cosmological origin can be detected from the foreground-subtracted CMB anisotropy map at the 2 σlevel if the measured genus shift parameter |Δν| >= 0.02 (0.04) or if the measured genus asymmetry parameter |Δg| >= 0.03 (0.08) on a 0.3 (1.0) deg FWHM scale.

preprint2000arXiv

Topology from the Simulated Sloan Digital Sky Survey

We measure the topology (genus curve) of the galaxy distribution in a mock redshift catalog designed to resemble the upcoming Sloan Digital Sky Survey (SDSS). The catalog, drawn from a large N-body simulation of a Lambda-CDM cos- mological model, mimics the anticipated spectroscopic selection procedures of the SDSS in some detail. Sky maps, redshift slices, and 3-D contour maps of the mock survey reveal a rich and complex structure, including networks of voids and superclusters that resemble the patterns seen in the CfA redshift survey and the Las Campanas Redshift Survey (LCRS). The 3-D genus curve can be measured from the simulated catalog with superb precision; this curve has the general shape predicted for Gaussian, random phase initial conditions, but the error bars are small enough to demonstrate with high significance the subtle departures from this shape caused by non-linear gravitational evolution. These distortions have the form predicted by Matsubara's (1994) perturbative anal- ysis, but they are much smaller in amplitude. We also measure the 3-D genus curve of the radial peculiar velocity field measured by applying distance- indicator relations (with realistic errors) to the mock catalog. This genus curve is consistent with the Gaussian random phase prediction, though it is of relatively low precision because of the large smoothing length required to overcome noise in the measured velocity field. Finally, we measure the 2-D topology in redshift slices, similar to early slices from the SDSS and to slices already observed in the LCRS. The genus curves of these slices are consistent with the observed genus curves of the LCRS, providing further evidence in favor of the inflationary CDM model with Omega_M~0.4. The catalog is publicly available at http://www.astronomy.ohio-state.edu/~dhw/sdss.html.

preprint2000arXiv

Topology of the Galaxy Distribution in the Hubble Deep Fields

We have studied topology of the distribution of the high redshift galaxies identified in the Hubble Deep Field (HDF) North and South. The two-dimensional genus is measured from the projected distributions of the HDF galaxies at angular scales from $3.8''$ to $ 6.1''$. We have also divided the samples into three redshift slices with roughly equal number of galaxies using photometric redshifts to see possible evolutionary effects on the topology. The genus curve of the HDF North clearly indicates clustering of galaxies over the Poisson distribution while the clustering is somewhat weaker in the HDF South. This clustering is mainly due to the nearer galaxies in the samples. We have also found that the genus curve of galaxies in the HDF is consistent with the Gaussian random phase distribution with no significant redshift dependence.