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.

As a physical and sufficient compression of the full CMB data, the CMB distance priors, or shift parameters, have been widely used and provide a convenient way to include CMB data when obtaining cosmological constraints. In this paper, we revisit this data vector and examine its stability under different cosmological models. We find that the CMB distance priors are an accurate substitute for the full CMB data when probing dark energy dynamics. This is true when the primordial power spectrum model is directly generalized from the power spectrum of the model used in the derivation of the distance priors from the CMB data. We discover a difference when a non-flat model with the untilted primordial inflation power spectrum is used to measure the distance priors. This power spectrum is a radical change from the more conventional tilted primordial power spectrum and violates fundamental assumptions for the reliability of the CMB shift parameters. We also investigate the performance of CMB distance priors when the sum of neutrino masses $\sum m_{\nu}$ and the effective number of relativistic species $N_{\text{eff}}$ are allowed to vary. Our findings are consistent with earlier results: the neutrino parameters can change the measurement of the sound horizon from CMB data, and thus the CMB distance priors. We find that when the neutrino model is allowed to vary, the cold dark matter density $\omega_{c}$ and $N_{\text{eff}}$ need to be included in the set of parameters that summarize CMB data, in order to reproduce the constraints from the full CMB data. We present an updated and expanded set of CMB distance priors which can reproduce constraints from the full CMB data within $1\sigma$, and are applicable to models with massive neutrinos, as well as non-standard cosmologies.

Horizon Run 5 (HR5) is a cosmological hydrodynamics simulation which captures the properties of the Universe on a Gpc scale while achieving a resolution of 1 kpc. This enormous dynamic range allows us to simultaneously capture the physics of the cosmic web on very large scales and account for the formation and evolution of dwarf galaxies on much smaller scales. Inside the simulation box we zoom-in on a high-resolution cuboid region with a volume of 1049x114x114 Mpc^3. The sub-grid physics chosen to model galaxy formation includes radiative heating/cooling, reionization, star formation, supernova feedback, chemical evolution tracking the enrichment of oxygen and iron, the growth of supermassive black holes and active galactic nuclei (AGN) feedback in the form of a dual jet-heating mode. For this simulation we implemented a hybrid MPI-OpenMP version of the RAMSES code, specifically targeted for modern many-core many thread parallel architectures. For the post-processing, we extended the Friends-of-Friend (FoF) algorithm and developed a new galaxy finder to analyse the large outputs of HR5. The simulation successfully reproduces many observations, such as the cosmic star formation history, connectivity of galaxy distribution and stellar mass functions. The simulation also indicates that hydrodynamical effects on small scales impact galaxy clustering up to very large scales near and beyond the baryonic acoustic oscillation (BAO) scale. Hence, caution should be taken when using that scale as a cosmic standard ruler: one should carefully understand the corresponding biases. The simulation is expected to be an invaluable asset for the interpretation of upcoming deep surveys of the Universe.

Horizon Run 5 (HR5) is a cosmological hydrodynamics simulation which captures the properties of the Universe on a Gpc scale while achieving a resolution of 1 kpc. This enormous dynamic range allows us to simultaneously capture the physics of the cosmic web on very large scales and account for the formation and evolution of dwarf galaxies on much smaller scales. Inside the simulation box we zoom-in on a high-resolution cuboid region with a volume of 1049x114x114 Mpc^3. The sub-grid physics chosen to model galaxy formation includes radiative heating/cooling, reionization, star formation, supernova feedback, chemical evolution tracking the enrichment of oxygen and iron, the growth of supermassive black holes and active galactic nuclei (AGN) feedback in the form of a dual jet-heating mode. For this simulation we implemented a hybrid MPI-OpenMP version of the RAMSES code, specifically targeted for modern many-core many thread parallel architectures. For the post-processing, we extended the Friends-of-Friend (FoF) algorithm and developed a new galaxy finder to analyse the large outputs of HR5. The simulation successfully reproduces many observations, such as the cosmic star formation history, connectivity of galaxy distribution and stellar mass functions. The simulation also indicates that hydrodynamical effects on small scales impact galaxy clustering up to very large scales near and beyond the baryonic acoustic oscillation (BAO) scale. Hence, caution should be taken when using that scale as a cosmic standard ruler: one should carefully understand the corresponding biases. The simulation is expected to be an invaluable asset for the interpretation of upcoming deep surveys of the Universe.

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 $\xi_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 $\Omega_m$, the satellite velocity bias with respect to MBPs, $b_v^s$, and the merger-time-scale parameter $\alpha$. We find that $\alpha=1.5$ gives an excellent fit to the $w(r_p)$ and $\xi_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 $(\Omega_m, b_v^s, \alpha)$ by breaking the degeneracy among them. The resulting parameter values from all measurements are $(\Omega_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.

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 $\Omega_{\rm c}h^{2}$ and scalar spectral index $n_{\rm s}$. A strong negative degeneracy between $\Omega_{\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} \Omega_{\rm c}h^{2}$ -- we find $n_{\rm s}^{3/2} \Omega_{\rm c}h^{2} = 0.1121 \pm 0.0043$ after combining the LOWZ and CMASS data sets, assuming a flat $\Lambda$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}\Omega_{\rm c}h^{2} = 0.1139 \pm 0.0009$.

We use six tilted spatially-flat and untilted non-flat dark energy cosmological models in analyses of South Pole Telescope polarization (SPTpol) cosmic microwave background (CMB) data, alone and in combination with Planck 2015 CMB data and non-CMB data, namely, the Pantheon Type Ia supernovae apparent magnitudes, a collection of baryon acoustic oscillation data points, Hubble parameter measurements, and growth rates. Although the cosmological models that best-fit the Planck CMB and non-CMB data do not provide good fits to the SPTpol data, with the $\chi^2$'s exceeding the expected value, given the uncertainties, in each model the cosmological parameter constraints from the SPTpol data and from the Planck CMB and non-CMB data are largely mutually consistent. When the smaller angular scale SPTpol data are used jointly with either the Planck data alone or with the Planck CMB and the non-CMB data to constrain untilted non-flat models, spatially-closed models remain favored over their corresponding flat limits. When used in conjunction with Planck data, non-CMB data (baryon acoustic oscillation measurements in particular, from six experiments) have significantly more constraining power than the SPTpol data.

We have developed observation control software for the Immersion GRating INfrared Spectrometer (IGRINS) slit-viewing camera module, which maintains the position of an astronomical target on the spectroscopic slit. It is composed of several packages that monitor and control the system, acquire the images, and compensate for the tracking error by sending tracking feedback information to the telescope control system. For efficient development and maintenance of each software package, we have applied software engineering methods, i.e., a spiral software development with model-based design. It is not trivial to define the shape and center of astronomical object point spread functions (PSFs), which do not have symmetric Gaussian profiles in short exposure (<4 s) guiding images. Efforts to determine the PSF centroid are additionally complicated by the core saturation of bright guide stars. We have applied both a two-dimensional Gaussian fitting algorithm (2DGA) and center balancing algorithm (CBA) to identify an appropriate method for IGRINS in the near-infrared K-band. The CBA derives the expected center position along the slit-width by referencing the spillover flux ratio of the PSF wings on both sides of the slit. In this research, we have compared the accuracy and reliability of the CBA to the 2DGA by using data from IGRINS commissioning observations at McDonald Observatory. We find that the performance of each algorithm depends on the brightness of the targets and the seeing conditions, with the CBA performing better in typical observing scenarios. The algorithms and test results we present can be utilized with future spectroscopic slit observations in various observing conditions and for a variety of spectrograph designs.

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.

Cosmic microwave background (CMB) anisotropy (spatial inhomogeneity) data provide the tightest constraints on the Hubble constant, matter density, spatial curvature, and dark energy dynamics. Other data, sensitive to the evolution of only the spatially homogeneous part of the cosmological model, such as Type Ia supernova apparent magnitude, baryon acoustic oscillation distance, and Hubble parameter measurements, can be used in conjunction with the CMB data to more tightly constrain parameters. Recent joint analyses of CMB and such non-CMB data indicate that slightly closed spatial hypersurfaces are favored in nonflat untilted inflation models and that dark energy dynamics cannot be ruled out, and favor a smaller Hubble constant. We show that the constraints that follow from these non-CMB data alone are consistent with those that follow from the CMB data alone and so also consistent with, but weaker than, those that follow from the joint analyses of the CMB and non-CMB data.

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 $\beta = 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 $\beta$. We estimate the Fisher information that can be extracted from the Minkowski tensors, and find $W^{0,2}_{1}$ is more sensitive to $\beta$ than $W^{0,2}_{2}$, and a measurement of $W^{0,2}_{1}$ accurate to $\sim 1\%$ can yield a $\sim 4\%$ constraint on $\beta$. 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.

We use the physically-consistent tilted spatially-flat and untilted non-flat $\Lambda$CDM inflation models to constrain cosmological parameter values with the Planck 2015 cosmic microwave background (CMB) anisotropy data and recent Type Ia supernovae measurements, baryonic acoustic oscillations (BAO) data, growth rate observations, and Hubble parameter measurements. The most dramatic consequence of including the four non-CMB data sets is the significant strengthening of the evidence for non-flatness in the non-flat $\Lambda$CDM model, from 1.8$\sigma$ for the CMB data alone to 5.1$\sigma$ for the full data combination. The BAO data is the most powerful of the non-CMB data sets in more tightly constraining model parameter values and in favoring a spatially-closed Universe in which spatial curvature contributes about a percent to the current cosmological energy budget. The untilted non-flat $\Lambda$CDM model better fits the large-angle CMB temperature anisotropy angular spectrum and is more consistent with the Dark Energy Survey constraints on the current value of the rms amplitude of mass fluctuations ($\sigma_8$) as a function of the current value of the nonrelativistic matter density parameter ($\Omega_m$) but does not provide as good a fit to the smaller-angle CMB temperature anisotropy data as does the tilted flat-$\Lambda$CDM model. Some measured cosmological parameter values differ significantly between the two models, including the reionization optical depth and the baryonic matter density parameter, both of whose 2$\sigma$ ranges (in the two models) are disjoint or almost so.

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

We examine the effects of dark energy dynamics and spatial curvature on cosmic reionization by studying reionization in tilted spatially-flat and untilted non-flat XCDM and $\phi$CDM dynamical dark energy inflation models that best fit the Planck 2015 cosmic microwave background (CMB) anisotropy and a large compilation of non-CMB data. We carry out a detailed statistical study, based on a principal component analysis and a Markov chain Monte Carlo analysis of a compilation of lower-redshift reionization data, to estimate the uncertainties in the cosmological model reionization histories. We find that, irrespective of the nature of dark energy, there are significant differences between the reionization histories of the spatially-flat and non-flat models. Although both the flat and non-flat models can accurately match the low-redshift ($z\lesssim6$) reionization observations, there is a clear discrepancy between high-redshift ($z>7$) Lyman-$\alpha$ emitter data and the predictions from non-flat models. This is solely due to the fact that the non-flat models have a significantly larger electron scattering optical depth, $\tau_{\rm el}$, compared to the flat models, which requires an extended and much earlier reionization scenario supported by more high-redshift ionizing sources in the non-flat models. Non-flat models also require strong redshift evolution in the photon escape fraction, that can become unrealistically high ($\gtrsim1$) at some redshifts. However, $\tau_{\rm el}$ is about 0.9-$\sigma$ lower in the tilted flat $\Lambda$CDM model when the new Planck 2018 data are used and this reduction will partially alleviate the tension between the non-flat model predictions and the data.

Using SDSS DR15 to its full extent, we derived fundamental plane distances to over 317 000 early-type galaxies up to a redshift of 0.4. In addition to providing the largest sample of fundamental plane distances ever calculated, as well as a well calibrated group catalogue covering the entire SDSS spectroscopic footprint as far a redshift of 0.5, we present several improvements reaching beyond the traditional definition of the fundamental plane. In one approach, we adjusted the distances by removing systematic biases and selection effects in redshift-magnitude space, thereby greatly improving the quality of measurements. Alternatively, by expanding the traditional fundamental plane by additional terms, we managed to remove systematic biases caused by the selection of our SDSS spectroscopic galaxy sample as well as notably reducing its scatter. We discuss the advantages and caveats of these various methods and calibrations in detail. We found that improving the fundamental plane distance estimates beyond the established methods requires a delicate balancing act between various systematic biases and gains, but managed to reduce the uncertainty of our distance measurements by about a factor of two compared to the traditional fundamental plane.

We constrain tilted spatially-flat and untilted nonflat XCDM dynamical dark energy inflation parameterizations using Planck 2015 cosmic microwave background (CMB) anisotropy data and recent baryonic acoustic oscillations distance measurements, Type Ia supernovae data, Hubble parameter observations, and growth rate measurements. Inclusion of the four non-CMB data sets results in a significant strengthening of the evidence for nonflatness in the nonflat XCDM model from 1.1$\sigma$ for the CMB data alone to 3.4$\sigma$ for the full data combination. In this untilted nonflat XCDM case the data favor a spatially-closed model in which spatial curvature contributes a little less than a percent of the current cosmological energy budget; they also mildly favor dynamical dark energy over a cosmological constant at 1.2$\sigma$. These data are also better fit by the flat-XCDM parameterization than by the standard $\Lambda$CDM model, but only at 0.3$\sigma$ significance. Current data is unable to rule out dark energy dynamics. The nonflat XCDM parameterization is compatible with the Dark Energy Survey limits on the present value of the rms mass fluctuations amplitude ($\sigma_8$) as a function of the present value of the nonrelativistic matter density parameter ($\Omega_m$), however it does not provide as good a fit to the higher multipole CMB temperature anisotropy data as does the standard tilted flat-$\Lambda$CDM model. A number of measured cosmological parameter values differ significantly when determined using the tilted flat-XCDM and the nonflat XCDM parameterizations, including the baryonic matter density parameter and the reionization optical depth.

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 $\Lambda$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) $\alpha$, (2) closed-to-connected triple fraction (the so-called transitivity) $\tau_\Delta$, 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.

This study presents a numerical analysis of the topology of a set of cosmologically interesting three-dimensional Gaussian random fields in terms of their Betti numbers $\beta_0$, $\beta_1$ and $\beta_2$. We show that Betti numbers entail a considerably richer characterization of the topology of the primordial density field. Of particular interest is that Betti numbers specify which topological features - islands, cavities or tunnels - define its spatial structure. A principal characteristic of Gaussian fields is that the three Betti numbers dominate the topology at different density ranges. At extreme density levels, the topology is dominated by a single class of features. At low levels this is a \emph{Swiss-cheeselike} topology, dominated by isolated cavities, at high levels a predominantly \emph{Meatball-like} topology of isolated objects. At moderate density levels, two Betti number define a more \emph{Sponge-like} topology. At mean density, the topology even needs three Betti numbers, quantifying a field consisting of several disconnected complexes, not of one connected and percolating overdensity. A {\it second} important aspect of Betti number statistics is that they are sensitive to the power spectrum. It reveals a monotonic trend in which at a moderate density range a lower spectral index corresponds to a considerably higher (relative) population of cavities and islands. We also assess the level of complementary information that Betti numbers represent, in addition to conventional measures such as Minkowski functionals. To this end, we include an extensive description of the Gaussian Kinematic Formula (GKF), which represents a major theoretical underpinning for this discussion.

Ultrarefractory (UR) phases in CAIs could have formed at higher T compared to common CAI minerals and thus they potentially provide constraints on very high-T processes in the solar nebula. We report a detailed characterization of an UR phase davisite bearing CAI from a reduced type CV chondrite. Absence of secondary iron- and/or alkali-rich phases and occurrence of LIME olivine indicate that primitive chemical and isotopic compositions are preserved in the CAI. Davisite occur only in one lithological unit that consists of three chemically and isotopically distinct parts: i) $^{16}$O-poor regions with reversely-zoned melilite and davisite; ii) $^{16}$O-rich regions consisting of unzoned, gehlenitic melilite, diopside and spinel; and iii) spinel framboids composed of $^{16}$O-rich spinel and $^{16}$O-poor melilite. Random distribution of chemical and isotopic heterogeneities with sharp boundaries in the CAI indicates its formation by an aggregation of mineral assemblages formed and processed separately at different time and/or space. Although isotope exchange prior to the final agglomeration of the CAI cannot be ruled out, we suggest that modification of chemical and isotopic composition of porous CAI precursors or aggregation of isotopically distinct mineral assemblages are alternative scenarios for the origin of O-isotopic heterogeneity in CAIs. In either case, coexistence of spatially and/or temporally distinct $^{16}$O-rich and -poor gaseous reservoirs at the earliest stage of the solar system formation is required. The grain-scale oxygen isotopic disequilibrium in the CAI indicate that post-formation heating of the CAI was short, which can be achieved by rapid outward transport of the CAI. High Ti$^{3+}$/Ti$^{\mathrm{tot}}$ ratios of pyroxene and presence of LIME olivine document that the entire CAI formation process took place under highly reducing conditions.

We study the rotational motion of the galaxy cluster Abell 2107 at redshift z = 0.04 and its connection to nearby large-scale structures using a large amount of spectroscopic redshift data. By combining 978 new redshifts from the MMT/Hectospec observations with data in the literature, we construct a large sample of 1968 galaxies with measured redshifts at clustercentric radius $R < 60^\prime$. Our sample has high (80%) and spatially uniform completeness at r-band apparent magnitude $m_{r,\textrm{Petro},0}$ < 19.1. We first apply the caustic method to the sample and identify 285 member galaxies of Abell 2107 at $R < 60^\prime$. Then the rotation amplitude and the position angle of rotation axis are measured. We find that the member galaxies show strong global rotation at $R < 20^\prime$ ($v_\textrm{rot}/\sigma_v$ ~0.6) with a significance of > 3.8$\sigma$, which is confirmed by two independent methods. The rotation becomes weaker in outer regions. There are at least five filamentary structures that are connected to the cluster and that consist of known galaxy groups. These structures are smoothly connected to the cluster, which seem to be inducing the global rotation of the cluster through inflow of galaxies.

We present hydrodynamic equations with relativistic pressure and velocity in the presence of weak gravity, in a cosmological context. Previously we consistently derived special relativistic hydrodynamic equations with weak gravity in Minkowski background. With the relativistic pressure and velocity one cannot derive the cosmological counterpart by a simple transformation from equations in the Minkowski background. Here we present a proper derivation. We point out the potential importance of relativistic pressure and velocity in gravitational lensing.

We study the dynamical state and properties of galaxies and groups in the supercluster SClA2142 that has a collapsing core, to understand its possible formation and evolution. We find the substructure of galaxy groups using normal mixture modelling. We have used the projected phase space (PPS) diagram, spherical collapse model, clustercentric distances, and magnitude gap between the brightest galaxies in groups to study the dynamical state of groups and to analyse group and galaxy properties. We compared the alignments of groups and their brightest galaxies with the supercluster axis. The supercluster core has a radius of about $8 h^{-1}$Mpc and total mass $M_{\mathrm{tot}} \approx 2.3\times~10^{15}h^{-1}M_\odot$ and is collapsing. Galaxies in groups on the supercluster axis have older stellar populations than off-axis groups, with median stellar ages $4 - 6$ and $< 4$Gyr, correspondingly. The cluster A2142 and the group Gr8 both host galaxies with the oldest stellar populations among groups in SClA2142 having the median stellar age $t > 8$Gyr. Recently quenched galaxies and active galactic nuclei (AGNs) are mostly located at virial radii or in merging regions of groups, and at clustercentric distances $D_c \approx 6 h^{-1}$ Mpc. The most elongated groups lie along the supercluster axis and are aligned with it. Magnitude gaps between the brightest galaxies of groups are less than one magnitude, suggesting that groups in SClA2142 are dynamically young. The collapsing core of the supercluster, infall of galaxies and groups, and possible merging groups, which affect galaxy properties and may trigger the activity of AGNs, show how the whole supercluster is evolving.

We present interesting aspects of a modified theory of gravity called Ho\v{r}ava-Lifshitz (HL) gravity that can be constrained by the phenomenological and observational prospects in terms of neutron star (NS) structure and equation-of-state models. The deformation in HL gravity from general relativity (GR) can change typical features of the NS structure. In this Letter, we investigate the NS structure by deriving Tolman-Oppenheimer-Volkoff equation in HL gravity. We find that a NS in HL gravity with larger radius and heavier mass than a NS in GR remains stable without collapsing into a black hole.

We constrain spatially-flat tilted and nonflat untilted scalar field ($\phi$) dynamical dark energy inflation ($\phi$CDM) models by using Planck 2015 cosmic microwave background (CMB) anisotropy measurements and recent baryonic acoustic oscillation distance observations, Type Ia supernovae apparent magnitude data, Hubble parameter measurements, and growth rate data. We assume an inverse power-law scalar field potential energy density $V(\phi)=V_0 \phi^{-\alpha}$. We find that the combination of the CMB data with the four non-CMB data sets significantly improves parameter constraints and strengthens the evidence for nonflatness in the nonflat untilted $\phi$CDM case from $1.8\sigma$ for the CMB measurements only to more than $3.1\sigma$ for the combined data. In the nonflat untilted $\phi$CDM model current observations favor a spatially closed universe with spatial curvature contributing about two-thirds of a percent of the present cosmological energy budget. The flat tilted $\phi$CDM model is a 0.4$\sigma$ better fit to the data than is the standard flat tilted $\Lambda$CDM model: current data allow for the possibility that dark energy is dynamical. The nonflat tilted $\phi$CDM model is in better accord with the Dark Energy Survey bounds on the rms amplitude of mass fluctuations now ($\sigma_8$) as a function of the nonrelativistic matter density parameter now ($\Omega_m$) but it does not provide as good a fit to the larger-multipole Planck 2015 CMB anisotropy data as does the standard flat tilted $\Lambda$CDM model. A few cosmological parameter value measurements differ significantly when determined using the tilted flat and the untilted nonflat $\phi$CDM models, including the cold dark matter density parameter and the reionization optical depth.

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

We present fore-optics and calibration unit design of Devasthal Optical Telescope Integral Field Spectrograph (DOTIFS). DOTIFS fore-optics is designed to modify the focal ratio of the light and to match its plate scale to the physical size of Integral Field Units (IFUs). The fore-optics also delivers a telecentric beam to the IFUs on the telescope focal plane. There is a calibration unit part of which is combined with the fore-optics to have a light and compact system. We use Xenon-arc lamp as a continuum source and Krypton/Mercury-Neon lamps as wavelength calibration sources. Fore-optics and calibration unit shares two optical lenses to maintain compactness of the overall subsystem. Here we present optical and opto-mechanical design of the calibration unit and fore-optics as well as calibration scheme of DOTIFS.

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

JINGLE is a new JCMT legacy survey designed to systematically study the cold interstellar medium of galaxies in the local Universe. As part of the survey we perform 850um continuum measurements with SCUBA-2 for a representative sample of 193 Herschel-selected galaxies with M*>10^9Msun, as well as integrated CO(2-1) line fluxes with RxA3m for a subset of 90 of these galaxies. The sample is selected from fields covered by the Herschel-ATLAS survey that are also targeted by the MaNGA optical integral-field spectroscopic survey. The new JCMT observations combined with the multi-wavelength ancillary data will allow for the robust characterization of the properties of dust in the nearby Universe, and the benchmarking of scaling relations between dust, gas, and global galaxy properties. In this paper we give an overview of the survey objectives and details about the sample selection and JCMT observations, present a consistent 30 band UV-to-FIR photometric catalog with derived properties, and introduce the JINGLE Main Data Release (MDR). Science highlights include the non-linearity of the relation between 850um luminosity and CO line luminosity, and the serendipitous discovery of candidate z>6 galaxies.

We present perturbation initial conditions for two types of scalar field potential often used in the dark energy study: one inverse power-law and the other exponential. The solutions are presented in the presence of the $w =$ constant fluid ($w = 1/3$ for radiation fluid), a minimally coupled scalar field and a sub-dominating zero-pressure fluid (cold dark matter and baryon dust). We consider two gauge conditions, the $w$-fluid comoving gauge and the cold-dark-matter comoving gauge; solutions in the latter gauge are derived by the gauge transformation and this method can be applied to derive solutions in any other gauge condition.

We present a library of high-resolution (R $\equiv$ $\lambda$/$\Delta$$\lambda$ $\sim$ 45,000) and high signal-to-noise ratio (S/N $\geq$ 200) near-infrared spectra for stars of a wide range of spectral types and luminosity classes. The spectra were obtained with the Immersion GRating INfrared Spectrograph (IGRINS) covering the full range of the H (1.496-1.780 $\mu$m) and K (2.080-2.460 $\mu$m) atmospheric windows. The targets were primarily selected for being MK standard stars covering a wide range of effective temperatures and surface gravities with metallicities close to the Solar value. Currently, the library includes flux-calibrated and telluric-absorption-corrected spectra of 84 stars, with prospects for expansion to provide denser coverage of the parametric space. Throughout the H and K atmospheric windows, we identified spectral lines that are sensitive to $T_\mathrm{eff}$ or $\log g$ and defined corresponding spectral indices. We also provide their equivalent widths. For those indices, we derive empirical relations between the measured equivalent widths and the stellar atmospheric parameters. Therefore, the derived empirical equations can be used to calculate $T_\mathrm{eff}$ and $\log g$ of a star without requiring stellar atmospheric models.

We apply the Minkowski tensor statistics to three dimensional Gaussian random fields. Minkowski tensors contain information regarding the orientation and shape of excursion sets, that is not present in the scalar Minkowski functionals. They can be used to quantify globally preferred directions, and additionally provide information on the mean shape of subsets of a field. This makes them ideal statistics to measure the anisotropic signal generated by redshift space distortion in the low redshift matter density field. We review the definition of the Minkowski tensor statistics in three dimensions, focusing on two coordinate invariant quantities $W^{0,2}_{1}$ and $W^{0,2}_{2}$. We calculate the ensemble average of these $3 \times 3$ matrices for an isotropic Gaussian random field, finding that they are proportional to products of the identity matrix and a corresponding scalar Minkowski functional. We show how to numerically reconstruct $W^{0,2}_{1}$ and $W^{0,2}_{2}$ from discretely sampled fields and apply our algorithm to isotropic Gaussian fields generated from a linear $\Lambda$CDM matter power spectrum. We then introduce anisotropy by applying a linear redshift space distortion operator to the matter density field, and find that both $W^{0,2}_{1}$ and $W^{0,2}_{2}$ exhibit a distinct signal characterised by inequality between their diagonal components. We discuss the physical origin of this signal and how it can be used to constrain the redshift space distortion parameter $\Upsilon \equiv f/b$.

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

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS); the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data driven machine learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS website (www.sdss.org) has been updated for this release, and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020, and will be followed by SDSS-V.

We present a high-resolution, near-IR spectroscopic study of multiple outflows in the LkH$\alpha$ 234 star formation region using the Immersion GRating INfrared Spectrometer (IGRINS). Spectral mapping over the blueshifted emission of HH 167 allowed us to distinguish at least three separate, spatially overlapped, outflows in H${_2}$ and [Fe II] emission. We show that the H${_2}$ emission represents not a single jet, but complex multiple outflows driven by three known embedded sources: MM1, VLA 2, and VLA 3. There is a redshifted H${_2}$ outflow at a low velocity, $\VLSR$ $<$ $+$50 {\kms}, with respect to the systemic velocity of $\VLSR$ $=$ $-$11.5 {\kms}, that coincides with the H${_2}$O masers seen in earlier radio observations two arcseconds southwest of VLA 2. We found that the previously detected [Fe II] jet with $|$$\VLSR$$|$ $>$ 100 {\kms} driven by VLA 3B is also detected in H${_2}$ emission, and confirm that this jet has a position angle about 240$\degree$. Spectra of the redshifted knots at 14$\arcsec$$-$65$\arcsec$ northeast of LkH$\alpha$ 234 are presented for the first time. These spectra also provide clues to the existence of multiple outflows. We detected high-velocity (50$-$120 {\kms}) H${_2}$ gas in the multiple outflows around LkH$\alpha$ 234. Since these gases move at speeds well over the dissociation velocity ($>$ 40 {\kms}), the emission must originate from the jet itself rather than H${_2}$ gas in the ambient medium. Also, position-velocity diagrams and excitation diagram indicate that emission from knot C in HH 167 come from two different phenomena, shocks and photodissociation.

Late-type galaxies falling into a cluster would evolve being influenced by the interactions with both the cluster and the nearby cluster member galaxies. Most numerical studies, however, tend to focus on the effects of the former with little work done on those of the latter. We thus perform a numerical study on the evolution of a late-type galaxy interacting with neighboring early-type galaxies at high speed, using hydrodynamic simulations. Based on the information obtained from the Coma cluster, we set up the simulations for the case where a Milky Way-like late-type galaxy experiences six consecutive collisions with twice as massive early-type galaxies having hot gas in their halos at the closest approach distances of 15-65 kpc/h at the relative velocities of 1500-1600 km/s. Our simulations show that the evolution of the late-type galaxy can be significantly affected by the accumulated effects of the high-speed multiple collisions with the early-type galaxies, such as on cold gas content and star formation activity of the late-type galaxy, particularly through the hydrodynamic interactions between cold disk and hot gas halos. We find that the late-type galaxy can lose most of its cold gas after the six collisions and have more star formation activity during the collisions. By comparing our simulation results with those of galaxy-cluster interactions, we claim that the role of the galaxy-galaxy interactions on the evolution of late-type galaxies in clusters could be comparable with that of the galaxy-cluster interactions, depending on the dynamical history.

We have explored the evolution of gas distributions from cosmological simulations carried out using the RAMSES adaptive mesh refinement (AMR) code, to explore the effects of resolution on cosmological hydrodynamical simulations. It is vital to understand the effect of both the resolution of initial conditions and the final resolution of the simulation. Lower initial resolution simulations tend to produce smaller numbers of low mass structures. This will strongly affect the assembly history of objects, and has the same effect of simulating different cosmologies. The resolution of initial conditions is an important factor in simulations, even with a fixed maximum spatial resolution. The power spectrum of gas in simulations using AMR diverges strongly from the fixed grid approach - with more power on small scales in the AMR simulations - even at fixed physical resolution and also produces offsets in the star formation at specific epochs. This is because before certain times the upper grid levels are held back to maintain approximately fixed physical resolution, and to mimic the natural evolution of dark matter only simulations. Although the impact of hold back falls with increasing spatial and initial-condition resolutions, the offsets in the star formation remain down to a spatial resolution of 1 kpc. These offsets are of order of 10-20%, which is below the uncertainty in the implemented physics but are expected to affect the detailed properties of galaxies. We have implemented a new grid-hold-back approach to minimize the impact of hold back on the star formation rate.

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

We consider five dimensional conformal gravity theory which describes an anisotropic extra dimension. Reducing the theory to four dimensions yields Brans-Dicke theory with a potential and a hidden parameter $z$ which implements the anisotropy between the four dimensional spacetime and the extra dimension. We find that a range of value of the parameter $z$ can address the current dark energy density compared to the Planck energy density. Constraining the parameter $z$ and the other cosmological model parameters using the recent observational data consisting of the Hubble parameters, type Ia supernovae, and baryon acoustic oscillations, together with the Planck or WMAP 9-year data of the cosmic microwave background radiation, we find $z>-2.05$ for Planck data and $z>-2.09$ for WMAP 9-year data at 95\% confidence level. We also obtained constraints on the rate of change of the effective Newtonian constant~($G_{\rm eff}$) at present and the variation of $G_{\rm eff}$ since the epoch of recombination to be consistent with observation.

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

The hierarchical galaxy formation scenario in the Cold Dark Matter cosmogony with a non-vanishing cosmological constant and geometrically flat space has been very successful in explaining the large-scale distribution of galaxies. However, there have been claims that the scenario predicts too many satellite galaxies associated with massive galaxies compared to observations, called the missing satellite galaxy problem. Isolated groups of galaxies hosted by passively evolving massive early-type galaxies are ideal laboratories for finding the missing physics in the current theory. Here we report from a deep spectroscopic survey of such satellite systems that isolated massive early-type galaxies with no recent star formation through wet mergers or accretion have almost no satellite galaxies fainter than the r-band absolute magnitude of about Mr =-14. If only early-type satellites are used, the cutoff is at somewhat brighter magnitude of about Mr =-15. Such a cutoff has not been found in other nearby satellite galaxy systems hosted by late-type galaxies or those with merger features. Various physical properties of satellites depend strongly on the host-centric distance. Our observation indicates that the satellite galaxy luminosity function is largely determined by the interaction of satellites with the environment provided by their host, which sheds light on the missing satellite galaxy problem.

We employ the rank-2 {\em contour} Minkowski Tensor in two dimensions to probe length and time scales of ionized bubbles during the epoch of reionization. We demonstrate that the eigenvalues of this tensor provide excellent probes of the distribution of the sizes of ionized bubbles, and from it the characteristic bubble sizes, at different redshifts. We show that ionized bubbles are not circular, and hence not spherical in three dimensions, as is often assumed for simplified analytic arguments. We quantify their shape anisotropy by using the ratio of the two eigenvalues. The shape parameter provides the characteristic time epochs when bubble mergers begin and end. Our method will be very useful to reconstruct the reionization history using data of the brightness temperature field.

We generalize the translation invariant tensor-valued Minkowski Functionals which are defined on two-dimensional flat space to the unit sphere. We apply them to level sets of random fields. The contours enclosing boundaries of level sets of random fields give a spatial distribution of random smooth closed curves. We obtain analytic expressions for the ensemble expectation values for the matrix elements of the tensor-valued Minkowski Functionals for isotropic Gaussian and Rayleigh fields. We elucidate the way in which the elements of the tensor Minkowski Functionals encode information about the nature and statistical isotropy (or departure from isotropy) of the field. We then implement our method to compute the tensor-valued Minkowski Functionals numerically and demonstrate how they encode statistical anisotropy and departure from Gaussianity by applying the method to maps of the Galactic foreground emissions from the PLANCK data.

We study galaxy populations and search for possible merging substructures in the rich galaxy cluster A2142. Normal mixture modelling revealed in A2142 several infalling galaxy groups and subclusters. The projected phase space diagram was used to analyse the dynamics of the cluster and study the distribution of various galaxy populations in the cluster and subclusters. The cluster, supercluster, BCGs, and one infalling subcluster are aligned. Their orientation is correlated with the alignment of the radio and X-ray haloes of the cluster. Galaxies in the centre of the main cluster at the clustercentric distances $0.5~h^{-1}Mpc$ have older stellar populations (with the median age of $10 - 11$~Gyrs) than galaxies at larger clustercentric distances. Star-forming and recently quenched galaxies are located mostly in the infall region at the clustercentric distances $D_{\mathrm{c}} \approx 1.8~h^{-1}Mpc$, where the median age of stellar populations of galaxies is about $2$~Gyrs. Galaxies in A2142 have higher stellar masses, lower star formation rates, and redder colours than galaxies in other rich groups. The total mass in infalling groups and subclusters is $M \approx 6\times10^{14}h^{-1}M_\odot$, approximately half of the mass of the cluster, sufficient for the mass growth of the cluster from redshift $z = 0.5$ (half-mass epoch) to the present. The cluster A2142 may have formed as a result of past and present mergers and infallen groups, predominantly along the supercluster axis. Mergers cause complex radio and X-ray structure of the cluster and affect the properties of galaxies in the cluster, especially in the infall region. Explaining the differences between galaxy populations, mass, and richness of A2142, and other groups and clusters may lead to better insight about the formation and evolution of rich galaxy clusters.

Large deviation statistics is implemented to predict the statistics of cosmic densities in cylinders applicable to photometric surveys. It yields few percent accurate analytical predictions for the one-point probability distribution function (PDF) of densities in concentric or compensated cylinders; and also captures the density-dependence of their angular clustering (cylinder bias). All predictions are found to be in excellent agreement with the cosmological simulation Horizon Run 4 in the quasi-linear regime where standard perturbation theory normally breaks down. These results are combined with a simple local bias model that relates dark matter and tracer densities in cylinders and validated on simulated halo catalogues. This formalism can be used to probe cosmology with existing and upcoming photometric surveys like DES, Euclid or WFIRST containing billions of galaxies.

The mergers of galaxy clusters are the most energetic events in the universe after the Big Bang. With the increased availability of multi-object spectroscopy and X-ray data an ever increasing fraction of local clusters are recognised as exhibiting signs of recent or past merging events on various scales. Our goal is to probe how these mergers affect the evolution and content of their member galaxies. We specifically aim to answer the following questions: Is the quenching of star formation in merging clusters enhanced when compared with relaxed clusters? Is the quenching preceded by a (short lived) burst of star formation? We obtained optical spectroscopy of >400 galaxies in the field of the merging cluster Abell 520. We combine these observations with archival data to get a comprehensive picture of the state of star formation in the members of this merging cluster. Finally, we compare these observations with a control sample of 10 non-merging clusters at the same redshift from The Arizona Cluster Redshift Survey (ACReS). We split the member galaxies in passive, star forming or recently quenched depending on their spectra. The core of the merger shows a decreased fraction of star forming galaxies compared to clusters in the non-merging sample. This region, dominated by passive galaxies, is extended along the axis of the merger. We find evidence of rapid quenching of the galaxies during the core passage with no signs of a star burst on the time scales of the merger (~0.4Gyr). Additionally, we report the tentative discovery of an infalling group along the main filament feeding the merger, currently at ~2.5 Mpc from the merger centre. This group contains a high fraction of star forming galaxies as well as ~2/3 of all the recently quenched galaxies in our survey.

We study the star formation activity of nearby galaxies with bars using a sample of late-type galaxies at 0.02$\leq z \leq$ 0.05489 and $M_r <-19.5$ from the Sloan Digital Sky Survey. We compare the physical properties of strongly and weakly barred galaxies with those of non-barred galaxies that have stellar mass and redshift distributions similar to barred galaxies. We find that the star formation activity of strongly barred galaxies probed by starburstiness, $\it{g-r}$, NUV$-r$, and mid-infrared [3.4]$-$[12] colors is, on average, lower than that of non-barred galaxies. However, weakly barred galaxies do not show such a difference between barred and non-barred galaxies. The amounts of atomic and molecular gas in strongly barred galaxies are smaller than those of non-barred galaxies, and the gas metallicity is higher in strongly barred galaxies than in non-barred galaxies. The gas properties of weakly barred galaxies again show no difference from those of non-barred galaxies. We stack the optical spectra of barred and non-barred galaxies in several mass bins and fit to the stacked spectra with a spectral fitting code, STARLIGHT. We find no significant difference in stellar populations between barred and non-barred galaxies for both strongly and weakly barred galaxies. Our results are consistent with the idea that the star formation activity of barred galaxies is enhanced in the past along with significant gas consumption, and is currently lower than or similar to that of non-barred galaxies. The past star formation enhancement depends on the strength of bars.

We present analysis of near-infrared, high-resolution spectroscopy towards the Flat-spectrum YSO (Young Stellar Object) ESO H$\alpha$ 279a (1.5 solar mass) in the Serpens star forming region, at the distance of 429 pc. Using the Immersion GRating INfrared Spectrometer (IGRINS, R=45,000), we detect emission lines originating from the accretion channel flow, jet, and inner disk. Specifically, we identify hydrogen Brackett series recombination, [Fe II], [Fe III], [Fe IV], Ca I, Na I, H2, H2O and CO overtone emission lines. By modeling five bands of CO overtone emission lines, and the symmetric double-peaked line profile for Na I emission lines, we find that ESO H$\alpha$ 279a has an actively accreting Keplerian disk. From our Keplerian disk model, we find that Na I emission lines originate between 0.04 AU and 1.00 AU, while CO overtone emission lines are from the outer part of disk, in the range between 0.22 AU and 3.00 AU. It reveals that the neutral atomic Na gas is a good tracer of the innermost region of the actively accreting disk. We derive a mass accretion rate of 2-10x10^{-7} M_solar/yr from the measured Br_gamma emission luminosity of 1.78(+-0.31)x10^{31} erg/s.

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

We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016.

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

We have measured the probability distribution function (PDF) of cosmic matter density field from a suite of N-body simulations. We propose the generalized normal distribution of version 2 (Nv2) as an alternative fitting formula to the well-known log-normal distribution. We find that Nv2 provides significantly better fit than the log-normal distribution for all smoothing radii (2, 5, 10, 25 [Mpc/h]) that we studied. The improvement is substantial in the underdense regions. The development of non- Gaissianities in the cosmic matter density field is captured by continuous evolution of the skewness and shifts parameters of the Nv2 distribution. We present the redshift evolution of these parameters for aforementioned smoothing radii and various background cosmology models. All the PDFs measured from large and high-resolution N-body simulations that we use in this study can be obtained from a Web site at https://astro.kias.re.kr/jhshin.

Simple parameter-free analytic bias functions for the two-point correlation of densities in spheres at large separation are presented. These bias functions generalize the so-called Kaiser bias to the mildly non-linear regime for arbitrary density contrasts. The derivation is carried out in the context of large deviation statistics while relying on the spherical collapse model. A logarithmic transformation provides a saddle approximation which is valid for the whole range of densities and shown to be accurate against the 30 Gpc cube state-of-the-art Horizon Run 4 simulation. Special configurations of two concentric spheres that allow to identify peaks are employed to obtain the conditional bias and a proxy to BBKS extrema correlation functions. These analytic bias functions should be used jointly with extended perturbation theory to predict two-point clustering statistics as they capture the non-linear regime of structure formation at the percent level down to scales of about 10 Mpc/h at redshift 0. Conversely, the joint statistics also provide us with optimal dark matter two-point correlation estimates which can be applied either universally to all spheres or to a restricted set of biased (over- or underdense) pairs. Based on a simple fiducial survey, this estimator is shown to perform five times better than usual two-point function estimators. Extracting more information from correlations of different types of objects should prove essential in the context of upcoming surveys like Euclid, DESI, PFS or LSST.

CMB polarization data is usually analyzed using $E$ and $B$ modes because they are scalars quantities under rotations along the lines of sight and have distinct physical origins. We explore the possibility of using the Stokes parameters $Q$ and $U$ for complementary analysis and consistency checks in the context of searches for non-Gaussianity. We show that the Minkowski Functionals (MFs) of $Q,U$ are invariant under local rotations along the lines of sight even though $Q,U$ are spin-2 variables, for full sky analysis. The invariance does not hold for incomplete sky. For local type primordial non-Gaussianity, when we compare the non-Gaussian deviations of MFs for $Q,U$ to what is obtained for $E$ mode or temperature fluctuations, we find that the amplitude is about an order of magnitude lower and the shapes of the deviations are different. This finding can be useful in distinguishing local type non-Gaussianity from other origins of non-Gaussianity in the observed data. Lastly, we analyze the sensitivity of the amplitudes of the MFs for $Q$, $U$ and the number density of singularities of the total polarization intensity to the tensor-to-scalar ratio, $r$, and find that all of them decrease as $r$ increases.

We present the results from an extensive spectroscopic survey of the central region of the nearby galaxy cluster Abell 2199 at $z=0.03$. By combining 775 new redshifts from the MMT/Hectospec observations with the data in the literature, we construct a large sample of 1624 galaxies with measured redshifts at $R<30^\prime$, which results in high spectroscopic completeness at $r_{\rm petro,0}<20.5$ (77%). We use these data to study the kinematics and clustering of galaxies focusing on the comparison with those of the intracluster medium (ICM) from Suzaku X-ray observations. We identify 406 member galaxies of A2199 at $R<30^\prime$ using the caustic technique. The velocity dispersion profile of cluster members appears smoothly connected to the stellar velocity dispersion profile of the cD galaxy. The luminosity function is well fitted with a Schechter function at $M_r<-15$. The radial velocities of cluster galaxies generally agree well with those of the ICM, but there are some regions where the velocity difference between the two is about a few hundred kilometer per second. The cluster galaxies show a hint of global rotation at $R<5^\prime$ with $v_{\rm rot}=300{-}600\,\textrm{km s}^{-1}$, but the ICM in the same region do not show such rotation. We apply a friends-of-friends algorithm to the cluster galaxy sample at $R<60^\prime$ and identify 32 group candidates, and examine the spatial correlation between the galaxy groups and X-ray emission. This extensive survey in the central region of A2199 provides an important basis for future studies of interplay among the galaxies, the ICM and the dark matter in the cluster.

The homogeneity of matter distribution at large scales, known as the cosmological principle, is a central assumption in the standard cosmological model. The case is testable though, thus no longer needs to be a principle. Here we perform a test for spatial homogeneity using the Sloan Digital Sky Survey Luminous Red Galaxies (LRG) sample by counting galaxies within a specified volume with the radius scale varying up to 300 Mpc/h. We directly confront the large-scale structure data with the definition of spatial homogeneity by comparing the averages and dispersions of galaxy number counts with allowed ranges of the random distribution with homogeneity. The LRG sample shows significantly larger dispersions of number counts than the random catalogues up to 300 Mpc/h scale, and even the average is located far outside the range allowed in the random distribution; the deviations are statistically impossible to be realized in the random distribution. This implies that the cosmological principle does not hold even at such large scales. The same analysis of mock galaxies derived from the N-body simulation, however, suggests that the LRG sample is consistent with the current paradigm of cosmology, thus the simulation is also not homogeneous in that scale. We conclude that the cosmological principle is not in the observed sky and nor is demanded to be there by the standard cosmological world model. This reveals the nature of the cosmological principle adopted in the modern cosmology paradigm, and opens new field of research in theoretical cosmology.

We present a deep near-infrared spectrum of the Orion Bar Photodissociation Region (PDR) taken with the Immersion Grating INfrared Spectrometer (IGRINS) on the 2.7 m telescope at the McDonald Observatory. IGRINS has high spectral resolution (R~45000) and instantaneous broad wavelength coverage (1.45-2.45 microns), enabling us to detect 87 emission lines from rovibrationally excited molecular hydrogen (H_2) that arise from transitions out of 69 upper rovibration levels of the electronic ground state. These levels cover a large range of rotational and vibrational quantum numbers and excitation energies, making them an excellent probe of the excitation mechanisms of H_2 and physical conditions within the PDR. The Orion Bar PDR is thought to consist of cooler high density clumps or filaments (T=50-250 K, n_H = 10^5 - 10^7 cm^-3) embedded in a warmer lower density medium (T=250-1000 K, n_H=10^4 - 10^5 cm^-3). We fit a grid of simple constant-temperature and constant-density Cloudy models, which recreate the observed H_2 level populations well, to constrain the temperature to a range of 600 to 650 K and the density to n_H = 2.5 x 10^3 to 10^4 cm^-3. The best fit model gives T = 625 K and n_H = 5x10^3 cm^-3. This well constrained warm temperature is consistent with kinetic temperatures found by other studies for the Orion Bar's lower density medium. However, the range of densities well fit by the model grid is marginally lower than those reported by other studies. We could be observing lower density gas than the surrounding medium, or perhaps a density-sensitive parameter in our models is not properly estimated.

We have analyzed the temperature, velocity and density of H2 gas in NGC 7023 with a high-resolution near-infrared spectrum of the northwestern filament of the reflection nebula. By observing NGC 7023 in the H and K bands at R ~ 45,000 with the Immersion GRating INfrared Spectrograph (IGRINS), we detected 68 H2 emission lines within the 1" x 15" slit. The diagnostic ratios of 2-1 S(1)/1-0 S(1) is 0.41-0.56. In addition, the estimated ortho-to-para ratios (OPR) is 1.63-1.82, indicating that the H2 emission transitions in the observed region arises mostly from gas excited by UV fluorescence. Gradients in the temperature, velocity, and OPR within the observed area imply motion of the photodissociation region (PDR) relative to the molecular cloud. In addition, we derive the column density of H2 from the observed emission lines and compare these results with PDR models in the literature covering a range of densities and incident UV field intensities. The notable difference between PDR model predictions and the observed data, in high rotational J levels of v = 1, is that the predicted formation temperature for newly-formed H2 should be lower than that of the model predictions. To investigate the density distribution, we combine pixels in 1" x 1" areas and derive the density distribution at the 0.002 pc scale. The derived gradient of density suggests that NGC 7023 has a clumpy structure, including a high clump density of ~10^5 cm^-3 with a size smaller than ~5 x 10^-3 pc embedded in lower density regions of 10^3-10^4 cm^-3.

We study the properties of dark matter haloes in a wide range of modified gravity models, namely, $f(R)$, DGP, and interacting dark energy models. We study the effects of modified gravity and dark energy on the internal properties of haloes, such as the spin and the structural parameters. We find that $f(R)$ gravity enhance the median value of the Bullock spin parameter, but could not detect such effects for DGP and coupled dark energy. $f(R)$ also yields a lower median sphericity and oblateness, while coupled dark energy has the opposite effect. However, these effects are very small. We then study the interaction rate of haloes in different gravity, and find that only strongly coupled dark energy models enhance the interaction rate. We then quantify the enhancement of the alignment of the spins of interacting halo pairs by modified gravity. Finally, we study the alignment of the major axes of haloes with the large-scale structures. The alignment of the spins of interacting pairs of haloes in DGP and coupled dark energy models show no discrepancy with GR, while $f(R)$ shows a weaker alignment. Strongly coupled dark energy shows a stronger alignment of the halo shape with the large-scale structures.

The variations of galaxy stellar masses and colour-types with the distance to projected cosmic filaments are quantified using the precise photometric redshifts of the COSMOS2015 catalogue extracted from COSMOS field (2 deg$^{2}$). Realistic mock catalogues are also extracted from the lightcone of the cosmological hydrodynamical simulation Horizon-AGN. They show that the photometric redshift accuracy of the observed catalogue ($\sigma_z<0.015$ at $M_*>10^{10}{\rm M}_{\odot}$ and $z<0.9$) is sufficient to provide 2D filaments that closely match their projected 3D counterparts. Transverse stellar mass gradients are measured in projected slices of thickness 75 Mpc between $0.5< z <0.9$, showing that the most massive galaxies are statistically closer to their neighbouring filament. At fixed stellar mass, passive galaxies are also found closer to their filament while active star-forming galaxies statistically lie further away. The contributions of nodes and local density are removed from these gradients to highlight the specific role played by the geometry of the filaments. We find that the measured signal does persist after this removal, clearly demonstrating that proximity to a filament is not equivalent to proximity to an over-density. These findings are in agreement with gradients measured both in 2D or 3D in the Horizon-AGN simulation and those observed in the spectroscopic VIPERS survey (which rely on the identification of 3D filaments). They are consistent with a picture in which the influence of the geometry of the large-scale environment drives anisotropic tides which impact the assembly history of galaxies, and hence their observed properties.

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

We use the Horizon Run 4 cosmological N -body simulation to study the effects of distant and close interactions on the alignments of the shapes, spins, and orbits of targets haloes with their neighbours, and their dependence on the local density environment and neighbour separation. Interacting targets have a significantly lower spin and higher sphericity and oblateness than all targets. Interacting pairs initially have anti-parallel spins, but the spins develop parallel alignment as time goes on. Neighbours tend to evolve in the plane of rotation of the target, and in the direction of the major axis of prolate haloes. Moreover, interactions are preferentially radial, while pairs with non-radial orbits are preferentially prograde. The alignment signals are stronger at high-mass and for close separations, and independent on the large-scale density. Positive alignment signals are found at redshifts up to 4, and increase with decreasing redshifts. Moreover, the orbits tend to become prograde at low redshift, while no alignment is found at high redshift (z = 4).

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= \lambda/\Delta\lambda$ 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.

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

The missing historical record of the Cassiopeia A (Cas A) supernova (SN) event implies a large extinction to the SN, possibly greater than the interstellar extinction to the current SN remnant. Here we investigate the possibility that the guest star that appeared near Cas A in 1592--1593 in Korean history books could have been an `impostor' of the Cas A SN, i.e., a luminous transient that appeared to be a SN but did not destroy the progenitor star, with strong mass loss to have provided extra circumstellar extinction. We first review the Korean records and show that a spatial coincidence between the guest star and Cas A cannot be ruled out, as opposed to previous studies. Based on modern astrophysical findings on core-collapse SN, we argue that Cas A could have had an impostor and derive its anticipated properties. It turned out that the Cas A SN impostor must have been bright ($M_V =-14.7 \pm 2.2$ mag) and an amount of dust with visual extinction of $\ge 2.8\pm 2.2 $ mag should have formed in the ejected envelope and/or in a strong wind afterwards. The mass loss needs to have been spherically asymmetric in order to see the light echo from the SN event but not the one from the impostor event.

We report a study of the three-dimensional (3D) outflow structure of a 15$\arcsec$ $\times$ 13$\arcsec$ area around H$_{2}$ peak 1 in Orion KL with slit-scan observations (13 slits) using the Immersion Grating Infrared Spectrograph. The datacubes, with high velocity-resolution ($\sim$ 7.5 {\kms}) provide high contrast imaging within ultra-narrow bands, and enable the detection of the main stream of the previously reported H$_{2}$ outflow fingers. We identified 31 distinct fingers in H$_{2}$ 1$-$0 S(1) $\lambda$2.122 $\micron$ emission. The line profile at each finger shows multiple-velocity peaks with a strong low-velocity component around the systemic velocity at ${\VLSR}$ = $+$8 {\kms} and high velocity emission ($|$${\VLSR}$$|$ = 45$-$135 {\kms}) indicating a typical bow-shock. The observed radial velocity gradients of $\sim$ 4 {\kms} arcsec$^{-1}$ agree well with the velocities inferred from large-scale proper motions, where the projected motion is proportional to distance from a common origin. We construct a conceptual 3D map of the fingers with the estimated inclination angles of 57$\degree$$-$74$\degree$. The extinction difference ($\Delta$$A_{\rm v}$ $>$ 10 mag) between blueshifted and redshifted fingers indicates high internal extinction. The extinction, the overall angular spread and scale of the flow argue for an ambient medium with very high density (10$^{5}$$-$10$^{6}$ cm$^{-3}$), consistent with molecular line observations of the OMC core. The radial velocity gradients and the 3D distributions of the fingers together support the hypothesis of simultaneous, radial explosion of the Orion KL outflow.

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

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\rho_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.

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 $\Lambda$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 $\Lambda$CDM hinting that higher order statistics may prove to be a useful tool in the hunt for deviations from General Relativity.

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

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)

We investigate cosmology of massive electrodynamics and explore the possibility whether massive photon could provide an explanation of the dark energy. The action is given by the scalar-vector-tensor theory of gravity which is obtained by non-minimal coupling of the massive Stueckelberg QED with gravity and its cosmological consequences are studied by paying a particular attention to the role of photon mass. We find that the theory allows cosmological evolution where the radiation- and matter-dominated epochs are followed by a long period of virtually constant dark energy that closely mimics $\Lambda$CDM model and the main source of the current acceleration is provided by the nonvanishing photon mass governed by the relation $\Lambda\sim m^2$. A detailed numerical analysis shows that the nonvanishing photon mass of the order of $\sim 10^{-34}$ eV is consistent with the current observations. This magnitude is far less than the most stringent limit on the photon mass available so far, which is of the order of $m \leq 10^{-27}$eV.

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 $\Lambda$ Cold Dark Matter ($\Lambda$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 $\Lambda$CDM model.

From high resolution (R = 45,000), high signal-to-noise (S/N > 400) spectra gathered with the Immersion Grating Infrared Spectrograph (IGRINS) in the H and K photometric bands, we have derived elemental abundances of two bright, well-known metal-poor halo stars: the red giant HD 122563 and the subgiant HD 140283. Since these stars have metallicities approaching [Fe/H] = -3, their absorption features are generally very weak. Neutral-species lines of Mg, Si, S and Ca are detectable, as well as those of the light odd-Z elements Na and Al. The derived IR-based abundances agree with those obtained from optical-wavelength spectra. For Mg and Si the abundances from the infrared transitions are improvements to those derived from shorter wavelength data. Many useful OH and CO lines can be detected in the IGRINS HD 122563 spectrum, from which derived O and C abundances are consistent to those obtained from the traditional [O I] and CH features. IGRINS high resolutions H- and K-band spectroscopy offers promising ways to determine more reliable abundances for additional metal-poor stars whose optical features are either not detectable, or too weak, or are based on lines with analytical difficulties.

We present the results of high-resolution near-IR spectroscopy toward the multiple outflows around the Herbig Be star Lk{\Ha} 234 using the Immersion Grating Infrared Spectrograph (IGRINS). Previous studies indicate that the region around Lk{\Ha} 234 is complex, with several embedded YSOs and the outflows associated with them. In simultaneous H$-$ and K$-$band spectra from HH 167, we detected 5 {\FeII} and 14 H$_{2}$ emission lines. We revealed a new {\FeII} jet driven by radio continuum source VLA 3B. Position-velocity diagrams of H$_{2}$ 1$-$0 S(1) $\lambda$2.122 $\micron$ line show multiple velocity peaks. The kinematics may be explained by a geometrical bow shock model. We detected a component of H$_{2}$ emission at the systemic velocity (V$_{LSR}$ $=$ $-$10.2 {\kms}) along the whole slit in all slit positions, which may arise from the ambient photodissociation region. Low-velocity gas dominates the molecular hydrogen emission from knots A and B in HH 167, which is close to the systemic velocity, {\FeII} emission lines are detected at farther from the systemic velocity, at V$_{LSR}$ $=$ $-$100 $-$ $-$130 {\kms}. We infer that the H$_{2}$ emission arises from shocked gas entrained by a high-velocity outflow. Population diagrams of H$_{2}$ lines imply that the gas is thermalized at a temperature of 2,500 $-$ 3,000 K and the emission results from shock excitation.

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.

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

We present fully nonlinear and exact cosmological perturbation equations in the presence of multiple components of fluids and minimally coupled scalar fields. We ignore the tensor-type perturbation. The equations are presented without taking the temporal gauge condition in the Friedmann background with general curvature and the cosmological constant. For each fluid component we ignore the anisotropic stress. The multiple component nature, however, introduces the anisotropic stress in the collective fluid quantities. We prove the Newtonian limit of multiple fluids in the zero-shear gauge and the uniform-expansion gauge conditions, present the Newtonian hydrodynamic equations in the presence of general relativistic pressure in the zero-shear gauge, and present the fully nonlinear equations and the third-order perturbation equations of the nonrelativistic pressure fluids in the CDM-comoving gauge.

We investigate aspects of axion as a coherently oscillating massive classical scalar field by analyzing third order perturbations in Einstein's gravity in the axion-comoving gauge. The axion fluid has its characteristic pressure term leading to an axion Jeans scale which is cosmologically negligible for a canonical axion mass. Our classically derived axion pressure term in Einstein's gravity is identical to the one derived in the non-relativistic quantum mechanical context in the literature. We show that except for the axion pressure term, the axion fluid equations are exactly the same as the general relativistic continuity and Euler equations of a zero-pressure fluid up to third order perturbation. The general relativistic density and velocity perturbations of the CDM in the CDM-comoving gauge are exactly the same as the Newtonian perturbations to the second order (in all scales), and the pure general relativistic corrections appearing from the third order are numerically negligible (in all scales as well) in the current paradigm of concordance cosmology. Therefore, here we prove that, in the super-Jeans scale, the classical axion can be handled as the Newtonian CDM fluid up to third order perturbation. We also show that the axion fluid supports the vector-type (rotational) perturbation from the third order. Our analysis includes the cosmological constant.

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

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.

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{*}}^{\alpha}/$M$_{\mathrm{h}}$) with $\alpha=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.

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, $\xi(\mu)$, is measured at different redshifts. Assuming an incorrect cosmological model to convert galaxy redshifts to distances, the shape of $\xi(\mu)$ 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 $\xi(\mu)$ exhibit redshift dependence. Similar to Li et.al (2014), we find the redshift-space distortions (RSD) caused by galaxy peculiar velocities, although significantly distorts $\xi(\mu)$, exhibit much less redshift evolution compared to the AP and volume effects. By focusing on the redshift dependence of $\xi(\mu)$, 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 $\delta \Omega_m\sim0.03$ and $\delta w\sim0.1$, and the volume effect leads to much tighter constraints of $\delta \Omega_m\sim0.007$ and $\delta w\sim0.035$.

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.

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.

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.

Interactions such as mergers and flybys play a fundamental role in shaping galaxy morphology. Using the Horizon Run 4 cosmological N-body simulation, we studied the frequency and type of halo interactions, and their redshift evolution as a function of the environment defined by the large-scale density, pair separation, mass ratio, and target halo mass. Most interactions happen at large-scale density contrast $\delta \approx 20$, regardless of the redshift, corresponding to groups and relatively dense part of filaments. However, the fraction of interacting target is maximum at $\delta \approx 1000$. We provide a new empirical fitting form for the interaction rate as a function of the halo mass, large-scale density, and redshift. We also report the existence of two modes of interactions from the distributions of mass ratio and relative distance, implying two different physical origins of the interaction. Satellite targets lose their mass as they proceed deeper into the host halo. The relative importance of these two trends strongly depends on the large-scale density, target mass, and redshift.

We investigate observational consequences of the early episodically dominating dark energy on the evolution of cosmological structures. For this aim, we introduce the minimally coupled scalar field dark energy model with the Albrecht-Skordis potential which allows a sudden ephemeral domination of dark energy component during the radiation or early matter era. The conventional cosmological parameters in the presence of such an early dark energy are constrained with WMAP and Planck cosmic microwave background radiation data including other external data sets. It is shown that in the presence of such an early dark energy the estimated cosmological parameters can deviate substantially from the currently known $\Lambda\textrm{CDM}$-based parameters, with best-fit values differing by several percents for WMAP and by a percent level for Planck data. For the latter case, only a limited amount of dark energy with episodic nature is allowed since the Planck data strongly favors the $\Lambda\textrm{CDM}$ model. Compared with the conventional dark energy model, the early dark energy dominating near radiation-matter equality or at the early matter era results in the shorter cosmic age or the presence of tensor-type perturbation, respectively. Our analysis demonstrates that the alternative cosmological parameter estimation is allowed based on the same observations even in Einstein's gravity.

The $\Lambda$ cold dark matter ($\Lambda\textrm{CDM}$) model is currently known as the simplest cosmology model that best describes observations with minimal number of parameters. Here we introduce a cosmology model that is preferred over the conventional $\Lambda\textrm{CDM}$ one by constructing dark energy as the sum of the cosmological constant $\Lambda$ and the additional fluid that is designed to have an extremely short transient spike in energy density during the radiation-matter equality era and the early scaling behavior with radiation and matter densities. The density parameter of the additional fluid is defined as a Gaussian function plus a constant in logarithmic scale-factor space. Searching for the best-fit cosmological parameters in the presence of such a dark energy spike gives a far smaller chi-square value by about five times the number of additional parameters introduced and narrower constraints on matter density and Hubble constant compared with the best-fit $\Lambda\textrm{CDM}$ model. The significant improvement in reducing chi-square mainly comes from the better fitting of Planck temperature power spectrum around the third ($\ell \approx 800$) and sixth ($\ell \approx 1800$) acoustic peaks. The likelihood ratio test and the Akaike information criterion suggest that the model of dark energy spike is strongly favored by the current cosmological observations over the conventional $\Lambda\textrm{CDM}$ model. However, based on the Bayesian information criterion which penalizes models with more parameters, the strong evidence supporting the presence of dark energy spike disappears. Our result emphasizes that the alternative cosmological parameter estimation with even better fitting of the same observational data is allowed in the Einstein's gravity.

We use a volume-limited sample of quasars in the Sloan Digital Sky Survey (SDSS) DR7 quasar catalog to identify quasar groups and address their statistical significance. This quasar sample has a uniform selection function on the sky and nearly a maximum possible contiguous volume that can be drawn from the DR7 catalog. Quasar groups are identified by using the Friend-of-Friend algorithm with a set of fixed comoving linking lengths. We find that the richness distribution of the richest 100 quasar groups or the size distribution of the largest 100 groups are statistically equivalent with those of randomly-distributed points with the same number density and sky coverage when groups are identified with the linking length of 70 h-1Mpc. It is shown that the large-scale structures like the huge Large Quasar Group (U1.27) reported by Clowes et al. (2013) can be found with high probability even if quasars have no physical clustering, and does not challenge the initially homogeneous cosmological models. Our results are statistically more reliable than those of Nadathur (2013), where the test was made only for the largest quasar group. It is shown that the linking length should be smaller than 50 h-1Mpc in order for the quasar groups identified in the DR7 catalog not to be dominated by associations of quasars grouped by chance. We present 20 richest quasar groups identified with the linking length of 70 h-1Mpc for further analyses.

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]

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\sigma)^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.

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.

We have developed a parallel cosmological hydrodynamic simulation code designed for the study of formation and evolution of cosmological structures. The gravitational force is calculated using the TreePM method and the hydrodynamics is implemented based on the smoothed particle hydrodynamics. The initial displacement and velocity of simulation particles are calculated according to second-order linear perturbation theory using the power spectra of dark matter and baryonic matter. The initial background temperature is given by Recfast and the temperature uctuations at the initial particle position are determined by the adiabatic model. We use a time-limiter scheme over the individual time steps to capture shock-fronts and to ease the time-step tension between the shock and preshock particles. We also include the astrophysical gas processes of radiative heating/cooling, star formation, metal enrichment, and supernova feedback. We have tested the code in several standard cases such as one-dimensional Riemann problems, Kelvin-Helmholtz, and Sedov blast wave instability. Star formation on the galactic disk is investigated to check whether the Schmidt-Kennicutt relation is properly recovered. We also study global star formation histories at different simulation resolutions and compare them with observations.

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.

A novel method to characterize the topology of the early-universe intergalactic medium during the epoch of cosmic reionization is presented. The 21-cm radiation background from high redshift is analyzed through the calculation of the 2-dimensional (2D) genus. The radiative transfer of hydrogen-ionizing photons and ionization-rate equations are calculated in a suite of numerical simulations under various input parameters. The 2D genus is calculated from the mock 21-cm images of the high-redshift Universe. We construct the 2D genus curve by varying the threshold differential brightness temperature, and compare this to the 2D genus curve of the underlying density field. We find that (1) the 2D genus curve reflects the evolutionary track of cosmic reionization and (2) the 2D genus curve can discriminate between certain reionization scenarios and thus indirectly probe the properties of radiation-sources. Choosing the right beam shape of a radio antenna is crucial for this analysis. To this end, the Square Kilometre Array (SKA) is found to be a suitable apparatus for this analysis in terms of sensitivity, even though some deterioration of the data for this purpose is unavoidable under the planned size of the antenna core.

We investigate the cosmology of massive spinor electrodynamics when torsion is non-vanishing. A non-minimal interaction is introduced between the torsion and the vector field and the coupling constant between them plays an important role in subsequential cosmology. It is shown that the mass of the vector field and torsion conspire to generate dark energy and pressureless dark matter, and for generic values of the coupling constant, the theory effectively provides an interacting model between them with an additional energy density of the form $\sim 1/a^6$. The evolution equations mimic $\Lambda$CDM behavior up to $1/a^3$ term and the additional term represents a deviation from $\Lambda$CDM. We show that the deviation is compatible with the observational data, if it is very small. We find that the non-minimal interaction is responsible for generating an effective cosmological constant which is directly proportional to the mass squared of the vector field and the mass of the photon within its current observational limit could be the source of the dark energy.

Cosmology is entering an era of percent level precision due to current large observational surveys. This precision in observation is now demanding more accuracy from numerical methods and cosmological simulations. In this paper, we study the accuracy of $N$-body numerical simulations and their dependence on changes in the initial conditions and in the simulation algorithms. For this purpose, we use a series of cosmological $N$-body simulations with varying initial conditions. We test the influence of the initial conditions, namely the pre-initial configuration (preIC), the order of the Lagrangian perturbation theory (LPT), and the initial redshift, on the statistics associated with the large scale structures of the universe such as the halo mass function, the density power spectrum, and the maximal extent of the large scale structures. We find that glass or grid pre-initial conditions give similar results at $z\lesssim 2$. However, the initial excess of power in the glass initial conditions yields a subtle difference in the power spectra and the mass function at high redshifts. The LPT order used to generate the ICs of the simulations is found to play a crucial role. First-order LPT (1LPT) simulations underestimate the number of massive haloes with respect to second-order (2LPT) ones, typically by 2% at $10^{14} h^{-1} M_\odot$ for an initial redshift of 23, and the small-scale power with an underestimation of 6% near the Nyquist frequency for $z_\mathrm{ini} = 23$. Moreover, at higher redshifts, the high-mass end of the mass function is significantly underestimated in 1LPT simulations. On the other hand, when the LPT order is fixed, the starting redshift has a systematic impact on the low-mass end of the halo mass function.

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 $\pi$ 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 $\lambda =$ 15, 21, and 34 $h^{-1}{\rm Mpc}$, the least $\chi^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.

We prove that the axion as a coherently oscillating scalar field acts as a cold dark matter (CDM) to the second-order perturbations in all cosmological scales including the super-horizon scale. The proof is made in the axion-comoving gauge. For a canonical mass, the axion pressure term causes deviation from the CDM only on scales smaller than the Solar System size. Beyond such a small scale the equations of the axion fluid are the same as the ones of the CDM based on the CDM-comoving gauge which are exactly identical to the Newtonian equations to the second order. We also show that the axion fluid does not generate the rotational (vector-type) perturbation even to the second order. Thus, in the case of axion fluid, we have the relativistic/Newtonian correspondence to the second order, even considering the rotational perturbation. Our analysis is made in the presence of the cosmological constant, and can be easily extended to the realistic situation including other components of fluids and fields.

Large sets of objects with spectroscopic redshift measurements will be needed for imaging dark energy experiments to achieve their full potential, serving two goals:_training_, i.e., the use of objects with known redshift to develop and optimize photometric redshift algorithms; and_calibration_, i.e., the characterization of moments of redshift (or photo-z error) distributions. Better training makes cosmological constraints from a given experiment stronger, while highly-accurate calibration is needed for photo-z systematics not to dominate errors. In this white paper, we investigate the required scope of spectroscopic datasets which can serve both these purposes for ongoing and next-generation dark energy experiments, as well as the time required to obtain such data with instruments available in the next decade. Large time allocations on kilo-object spectrographs will be necessary, ideally augmented by infrared spectroscopy from space. Alternatively, precision calibrations could be obtained by measuring cross-correlation statistics using samples of bright objects from a large baryon acoustic oscillation experiment such as DESI. We also summarize the additional work on photometric redshift methods needed to prepare for ongoing and future dark energy experiments.

We present a study of the connection between the galactic spin parameter $\lambda_{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 $\lambda_{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 $\lambda_{d}$ values for the case of long bars, while the maximum for short bars is at high $\lambda_{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 $\lambda_{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.

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; $\beta_0$ is the number of connected regions, $\beta_1$ is the number of circular holes, and $\beta_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. $\beta_0$ dominates the high-threshold part of the genus curve measuring the abundance of high density regions (clusters). $\beta_1$ dominates the genus near the median thresholds which measures the topology of negatively curved iso-density surfaces, and $\beta_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.

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, $\alpha$. 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 $\alpha$, 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.

We estimate the power spectra of CMB temperature anisotropy in localized regions on the sky using the WMAP 7-year data. Here, we report that the north hat and the south hat regions at the high Galactic latitude (|b| >= 30 deg) show anomaly in the power spectrum amplitude around the third peak, which is statistically significant up to 3 sigma. We try to figure out the cause of the observed anomaly by analyzing the low Galactic latitude (|b|< 30 deg) regions where the galaxy contamination is expected to be stronger, and regions that are weakly or strongly dominated by the WMAP instrument noise. We also consider the possible effect of unresolved radio point sources. We found another but less statistically significant anomaly in the low Galactic latitude north and south regions whose behavior is opposite to the one at the high latitude. Our analysis shows that the observed north-south anomaly at high latitude becomes weaker on the regions with high number of observations (weak instrument noise), suggesting that the anomaly is significant at sky regions that are dominated by the WMAP instrument noise. We have checked that the observed north-south anomaly has weak dependences on the bin-width used in the power spectrum estimation and the Galactic latitude cut. We have also discussed the possibility that the detected anomaly may hinge on the particular choice of the multipole bin around the third peak. We anticipate that the issue of whether the anomaly is intrinsic one or due to the WMAP instrument noise will be resolved by the forthcoming Planck data.

We construct several Milky Way-like galaxy models containing a gas halo (as well as gaseous and stellar disks, a dark matter halo, and a stellar bulge) following either an isothermal or an NFW density profile with varying mass and initial spin. In addition, galactic winds associated with star formation are tested in some of the simulations. We evolve these isolated galaxy models using the GADGET-3 $N$-body/hydrodynamic simulation code, paying particular attention to the effects of the gas halo on the evolution. We find that the evolution of the models is strongly affected by the adopted gas halo component. The model without a gas halo shows an increasing star formation rate (SFR) at the beginning of the simulation for some hundreds of millions of years and then a continuously decreasing rate to the end of the run at 3 Gyr. On the other hand, the SFRs in the models with a gas halo emerge to be either relatively flat throughout the simulations or increasing over a gigayear and then decreasing to the end. The models with the more centrally concentrated NFW gas halo show overall higher SFRs than those with the isothermal gas halo of the equal mass. The gas accretion from the halo onto the disk also occurs more in the models with the NFW gas halo, however, this is shown to take place mostly in the inner part of the disk and not to contribute significantly to the star formation unless the gas halo has very high density at the central part. The rotation of a gas halo is found to make SFR lower in the model. The SFRs in the runs including galactic winds are found to be lower than the same runs but without winds. We conclude that the effects of a hot gaseous halo on the evolution of galaxies are generally too significant to be simply ignored, and expect that more hydrodynamical processes in galaxies could be understood through numerical simulations employing both gas disk and gas halo components.

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.

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.

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 ($\lambda_{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 $\lambda_{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.

Axion as a coherently oscillating scalar field is known to behave as a cold dark matter in all cosmologically relevant scales. For conventional axion mass with 10^{-5} eV, the axion reveals a characteristic damping behavior in the evolution of density perturbations on scales smaller than the solar system size. The damping scale is inversely proportional to the square-root of the axion mass. We show that the axion mass smaller than 10^{-24} eV induces a significant damping in the baryonic density power spectrum in cosmologically relevant scales, thus deviating from the cold dark matter in the scale smaller than the axion Jeans scale. With such a small mass, however, our basic assumption about the coherently oscillating scalar field is broken in the early universe. This problem is shared by other dark matter models based on the Bose-Einstein condensate and the ultra-light scalar field. We introduce a simple model to avoid this problem by introducing evolving axion mass in the early universe, and present observational effects of present-day low-mass axion on the baryon density power spectrum, the cosmic microwave background radiation (CMB) temperature power spectrum, and the growth rate of baryon density perturbation. In our low-mass axion model we have a characteristic small-scale cutoff in the baryon density power spectrum below the axion Jeans scale. The small-scale deviations from the cold dark matter model in both matter and CMB power spectra clearly differ from the ones expected in the cold dark matter model mixed with the massive neutrinos as a hot dark matter component.

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.

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.

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.

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

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.

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.

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.

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 ($\sigma\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 $\sigma>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 $\sigma$) 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.

We study characters of recent type Ia supernova (SNIa) data using evolving dark energy models with changing equation of state parameter w. We consider sudden-jump approximation of w for some chosen redshift spans with double transitions, and constrain these models based on Markov Chain Monte Carlo (MCMC) method using the SNIa data (Constitution, Union, Union2) together with baryon acoustic oscillation A parameter and cosmic microwave background shift parameter in a flat background. In the double-transition model the Constitution data shows deviation outside 1 sigma from LCDM model at low (z < 0.2) and middle (0.2 < z < 0.4) redshift bins whereas no such deviations are noticeable in the Union and Union2 data. By analyzing the Union members in the Constitution set, however, we show that the same difference is actually due to different calibration of the same Union sample in the Constitution set, and is not due to new data added in the Constitution set. All detected deviations are within 2 sigma from the LCDM world model. From the LCDM mock data analysis, we quantify biases in the dark energy equation of state parameters induced by insufficient data with inhomogeneous distribution of data points in the redshift space and distance modulus errors. We demonstrate that location of peak in the distribution of arithmetic means (computed from the MCMC chain for each mock data) behaves as an unbiased estimator for the average bias, which is valid even for non-symmetric likelihood distributions.

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 $\lambda$, which we estimate approximately through simple dynamical considerations. Using a volume limited sample, we find a considerable difference when comparing the empirical distributions of $\lambda$ for AGNs and non-AGN galaxies, the AGNs showing typically low $\lambda$ values and associated dispersions, while non-AGNs present higher $\lambda$ values and a broader distribution. A more striking difference is found when looking at $\lambda$ 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 $\lambda$ 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 $\lambda$ value. We also find that the fraction of AGN hosting galaxies in our sample strongly depends on galactic spin, increasing dramatically for decreasing $\lambda$. 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.

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 \delta < -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 \Lambda -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.

We study the effects of the cluster environment on galactic morphology by defining a dimensionless angular momentum parameter $\lambda_{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 $\lambda_{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 $\lambda_{d}$ values. By analysing trends in $\lambda_{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 ($\sigma>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.

The modified gravity with $f(R)=R^{1+\epsilon}$ ($\epsilon>0$) allows a scaling solution where the density of gravity sector follows the density of the dominant fluid. We present initial conditions of background and perturbation variables during the scaling evolution regime in the modified gravity. As a possible dark energy model we consider a gravity with a form $f(R)=R^{1+\epsilon}+qR^{-n}$ ($-1<n \le 0$) where the second term drives the late-time acceleration. We show that our $f(R)$ gravity parameters are very sensitive to the baryon perturbation growth and baryon density power spectrum, and present observational constraints on the model parameters. Our analysis suggests that only the parameter space extremely close to the $\Lambda\textrm{CDM}$ model is allowed.

Versions of parameterized pseudo-Newtonian gravity theories specially designed for cosmology have been introduced in recent cosmology literature. The modifications demand a zero-pressure fluid in the context of versions of modified Poisson-like equation with two different gravitational potentials. We consider such modifications in the context of relativistic gravity theories where the action is a general algebraic function of the scalar curvature, the scalar field, and the kinetic term of the field. In general it is not possible to isolate the zero-pressure fluid component simultaneously demanding a modification in the Poisson-like equation. Only in the small-scale limit we can realize some special forms of the attempted modifications. We address some loopholes in the possibility of showing non-Einstein gravity nature based on pseudo-Newtonian modifications in the cosmological context. We point out that future observations of gravitational weak lensing together with velocity perturbation can potentially test the validity of Einstein's gravity in cosmology context.

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.

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.

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 ($\Delta\nu$) 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 $\Lambda$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)

We study the dependence of satellite galaxy properties on the distance to the host galaxy and the orbital motion using the SDSS data. From SDSS DR7 we find 3515 isolated satellite systems of galaxies at z<0.03 that contain 8904 satellite galaxies. Using this sample we construct a catalog of 635 satellites associated with 215 host galaxies whose spin directions are determined by our inspection of the SDSS color images and/or by spectroscopic observations in the literature. We divide satellite galaxies into prograde and retrograde orbit subsamples depending on their orbital motion respect to the spin direction of the host. We find that the number of galaxies in prograde orbit is nearly equal to that of retrograde orbit galaxies: the fraction of satellites in prograde orbit is 50+/-2%. The velocity distribution of satellites with respect to their hosts is found almost symmetric: the median bulk rotation of satellites is -1+/- 8 km/s. It is found that the radial distribution of early-type satellites in prograde orbit is strongly concentrated toward the host while that of retrograde ones shows much less concentration. We also find the orbital speed of late-type satellites in prograde orbit increases as the projected distance to the host (R) decreases while the speed decreases for those in retrograde orbit. At R less than 0.1 times the host virial radius the orbital speed decreases in both prograde and retrograde orbit cases. Prograde satellites are on average fainter than retrograde satellites for both early and late morphological types. The u-r color becomes redder as R decreases for both prograde and retrograde orbit late-type satellites. The differences between prograde and retrograde orbit satellites may be attributed to their different origin or the different strength of physical processes that they have experienced through hydrodynamic interactions with their host 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.

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.

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

Despite the continuing importance of ground-based parallax measurements, few active programs remain. Because new members of the solar neighborhood tend towards later spectral types, infrared parallax programs are particularly desirable. Consequently, the astrometric quality of the new infrared camera, FanCam, developed by the Virginia Astronomical Instrumentation Laboratory (VAIL) for the 31-in (0.8-m) Tinsley reflector at Fan Mountain Observatory was assessed using 68 J-band exposures of an open cluster, NGC 2420, over a range of hour angles during 2005. Positions of 16 astrometric evaluation stars were measured and the repeatability of those positions was evaluated using the mean error in a single observation of unit weight. Overall, a precision of 1.3 +/- 0.7 microns in x (RA) and 1.3 +/- 0.8 microns in y (Dec) was attained, which corresponds to 0.04" +/- 0.02" in each axis. Although greater precision is expected from CCDs in the visual and near-infrared, this instrument can achieve precision similar to that of the ESO NTT infrared parallax program. Therefore, measuring parallaxes in the infrared would be feasible using this equipment. If initiated, such a program could provide essential distances for brown dwarfs and very low mass stars that would contribute significantly to the solar neighborhood census.

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 $\Lambda$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)

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

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 $\lambda$ 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

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.

We study a generalized version of Chaplygin gas as unified model of dark matter and dark energy. Using realistic theoretical models and the currently available observational data from the age of the universe, the expansion history based on the type Ia supernovae, the matter power spectrum, the cosmic microwave background radiation anisotropy power spectra, and the perturbation growth factor we put the unified model under observational test. As the model has only two free parameters in the flat Friedmann background [$\Lambda$CDM (cold dark matter) model has only one free parameter] we show that the model is already tightly constrained by currently available observations. The only parameter space extremely close to the $\Lambda$CDM model is allowed in this unified model.

We show the importance of properly including the perturbations of the dark energy component in the dynamical dark energy models based on a scalar field and modified gravity theories in order to meet with present and future observational precisions. Based on a simple scaling scalar field dark energy model, we show that observationally distinguishable substantial differences appear by ignoring the dark energy perturbation. By ignoring it the perturbed system of equations becomes inconsistent and deviations in (gauge-invariant) power spectra depend on the gauge choice.

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.

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

We search for isolated galaxies using a volume-limited sample of galaxies with 0.02< z < 0.04742 from SDSS DR7 supplemented by bright galaxies. We devise a diagnostic tool to select isolated galaxies in different environments using the projected separation (r_p) normalized by the virial radius of the nearest neighbor (r_{vir,nei}) and the local background density. We find that the isolation condition of r_p>r_{vir,nei} and \rho <\bar{\rho} well segregates the CIG galaxies. We confirm the morphology conformity between the host and their satellites, which suggests importance of hydrodynamic interaction among galaxies within their virial radii in galaxy evolution.

Using a sample of galaxy groups found in the Sloan Digital Sky Survey Data Release 4, we measure the following four types of alignment signals: (1) the alignment between the distributions of the satellites of each group relative to the direction of the nearest neighbor group (NNG); (2) the alignment between the major axis direction of the central galaxy of the host group (HG) and the direction of the NNG; (3) the alignment between the major axes of the central galaxies of the HG and the NNG; and (4) the alignment between the major axes of the satellites of the HG and the direction of the NNG. We find strong signal of alignment between the satellite distribution and the orientation of central galaxy relative to the direction of the NNG, even when the NNG is located beyond $3r_{\rm vir}$ of the host group. The major axis of the central galaxy of the HG is aligned with the direction of the NNG. The alignment signals are more prominent for groups that are more massive and with early type central galaxies. We also find that there is a preference for the two major axis of the central galaxies of the HG and NNG to be parallel for the system with both early central galaxies, however not for the systems with both late type central galaxies. For the orientation of satellite galaxies, we do not find any significant alignment signals relative to the direction of the NNG. From these four types of alignment measurements, we conclude that the large scale environment traced by the nearby group affects primarily the shape of the host dark matter halo, and hence also affects the distribution of satellite galaxies and the orientation of central galaxies. (Shortened version)

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.

We present measurements of the clustering of hot and cold patches in the microwave background sky as measured from the Wilkinson Microwave Anisotropy Probe (WMAP) five-year data. These measurements are compared with theoretical predictions which assume that the cosmological signal obeys Gaussian statistics. We find significant differences from the simplest Gaussian-based prediction. However, the measurements are sensitive to the fact that the noise is spatially inhomogeneous (e.g., because different parts of the sky were observed for different lengths of time). We show how to account for this spatial inhomogeneity when making predictions. Differences from the Gaussian-based expectation remain even after this more careful accounting of the noise. In particular, we note that hot and cold pixels cluster differently within the same temperature thresholds at few-degree scales. While these findings may indicate primordial non-Gaussianity, we discuss other plausible explanations for these discrepancies. In addition, we find some deviations from Gaussianity at sub-degree scales, especially in the W band, whose origin may be associated with extragalactic dust emission.

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.

We study the galaxy morphology-luminosity-environmental relation and its redshift evolution using a spectroscopic sample of galaxies in the Great Observatories Origins Deep Survey (GOODS). In the redshift range of $0.4\leq z\leq1.0$ we detect conformity in morphology between neighboring galaxies. The realm of conformity is confined within the virialized region associated with each galaxy plus dark matter halo system. When a galaxy is located within the virial radius of its nearest neighbor galaxy, its morphology strongly depends on the neighbor's distance and morphology: the probability for a galaxy to be an early type ($f_E$) strongly increases as it approaches an early-type neighbor, but tends to decrease as it approaches a late-type neighbor. We find that $f_E$ evolves much faster in high density regions than in low density regions, and that the morphology-density relation becomes significantly weaker at $z\approx 1$. This may be because the rate of galaxy-galaxy interactions is higher in high density regions, and a series of interactions and mergers over the course of galaxy life eventually transform late types into early types. We find more isolated galaxies are more luminous, which supports luminosity transformation through mergers at these redshifts. Our results are consistent with those from nearby galaxies, and demonstrate that galaxy-galaxy interactions have been strongly affecting the galaxy evolution over a long period of time.

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.

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.

Using a large sample of local galaxies (144,940) with -17.5< Mr <-22 and 0.025< z <0.107, selected from SDSS DR5, we compare AGN host galaxies with non-AGN galaxies at matched luminosity, velocity dispersion, color, color gradient, or concentration index, to investigate how AGN activity is related with galaxy properties. The AGN sample is composed of Type II AGNs identified with flux ratios of narrow-emission lines with S/N > 6. We find that the fraction of galaxies hosting an AGN (f_AGN) depends strongly on morphology together with color, and very weakly on luminosity or velocity dispersion of host galaxies. In particular, f_AGN of early-type galaxies is almost independent of luminosity nor velocity dispersion when color is fixed. The host galaxy color preferred by AGNs is u-r ~2.0 for early-type hosts and u-r=2.0-2.4 for late-type hosts. This trend suggests that AGNs are dominantly hosted by intermediate-mass late-type galaxies. We also investigate how the accretion power varies with galaxy properties. While the Eddington ratio ([OIII] line luminosity normalized by black hole mass) ranges over three orders of magnitude for both morphological types, late-type galaxies are the dominant hosts over all AGN power. Among late-type galaxies, bluer color galaxies host higher power AGNs. These results are consistent with a scenario that more massive and redder galaxies are harder to host AGNs since these galaxiesalready consumed gas at the center or do not have sufficient gas supply to feed the black hole. In contrast, intermediate-mass, intermediate-color, and more concentrated late-type galaxies are more likely to host AGNs, implying that perhaps some fraction of low-mass, blue, and less concentrated lat e-type galaxiesmay not host massive black holes or may host very low-power AGNs.

We measure the 3D genus topology of large scale structure using Luminous Red Galaxies (LRGs) in the Sloan Digital Sky Survey and find it consistent with the Gaussian random phase initial conditions expected from the simplest scenarios of inflation. This studies 3D topology on the largest scales ever obtained. The topology is sponge-like. We measure topology in two volume-limited samples: a dense shallow sample studied with smoothing length of 21h^{-1}Mpc, and a sparse deep sample studied with a smoothing length of 34h^{-1}Mpc. The amplitude of the genus curve is measured with 4% uncertainty. Small distortions in the genus curve expected from non-linear biasing and gravitational effects are well explained (to about 1-\sigma accuracy) by N-body simulations using a subhalo-finding technique to locate LRGs. This suggests the formation of LRGs is a clean problem that can be modeled well without any free fitting parameters. This bodes well for using LRGs to measure the characteristic scales such as the baryon oscillation scale in future deep redshift surveys.

The Yuan-Tseh Lee Array for Microwave Background Anisotropy (AMiBA) is the first interferometer dedicated to studying the cosmic microwave background (CMB) radiation at 3mm wavelength. The choice of 3mm was made to minimize the contributions from foreground synchrotron radiation and Galactic dust emission. The initial configuration of seven 0.6m telescopes mounted on a 6-m hexapod platform was dedicated in October 2006 on Mauna Loa, Hawaii. Scientific operations began with the detection of a number of clusters of galaxies via the thermal Sunyaev-Zel'dovich effect. We compare our data with Subaru weak lensing data in order to study the structure of dark matter. We also compare our data with X-ray data in order to derive the Hubble constant.

Near-infrared photometric observations of the Type IIn SN 2005ip in NGC 2906 reveal large fluxes (>1.3 mJy) in the K_s-band over more than 900 days. While warm dust can explain the late-time K_s-band emission of SN 2005ip, the nature of the dust heating source is ambiguous. Shock heating of pre-existing dust by post-shocked gas is unlikely because the forward shock is moving too slowly to have traversed the expected dust-free cavity by the time observations first reveal the K_s emission. While an infrared light echo model correctly predicts a near-infrared luminosity plateau, heating dust to the observed temperatures of ~1400-1600 K at a relatively large distance from the supernova (> 10^{18} cm) requires an extraordinarily high early supernova luminosity (~1 X 10^{11} L_solar). The evidence instead favors condensing dust in the cool, dense shell between the forward and reverse shocks. Both the initial dust temperature and the evolutionary trend towards lower temperatures are consistent with this scenario. We infer that radiation from the circumstellar interaction heats the dust. While this paper includes no spectroscopic confirmation, the photometry is comparable to other SNe that do show spectroscopic evidence for dust formation. Observations of dust formation in SNe are sparse, so these results provide a rare opportunity to consider SNe Type IIn as dust sources.

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\Delta (M_{r,obs}+20.77+4.93\Delta)})/2\Delta, where \Delta = 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\sigma 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.

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.

We use a sample of galaxies from the Sloan Digital Sky Survey (SDSS) to search for correlations between the $\lambda$ spin parameter and the environment and mass of galaxies. In order to calculate the total value of $\lambda$ for each observed galaxy, we employed a simple model of the dynamical structure of the galaxies, which allows a rough estimate of the value of $\lambda$ using only readily obtainable observables from the luminous galaxies. Use of a large volume-limited sample (upwards of 11,000) allows reliable inferences of mean values and dispersions of $\lambda$ distributions. We find, in agreement with some N-body cosmological simulations, no significant dependence of $\lambda$ on the environmental density of the galaxies. For the case of mass, our results show a marked correlation with $\lambda$, in the sense that low-mass galaxies present both higher mean values of $\lambda$ and associated dispersions, than high-mass galaxies. These results provide interesting constrain on the mechanisms of galaxy formation and acquisition of angular momentum, a valuable test for cosmological models.

We present a comprehensive study of the nature of the SDSS galaxies divided into various 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, returning a total of 16 fine classes of galaxies. We examine the luminosity dependence of seven physical parameters of galaxies in each class. We find that more than half of red early-type galaxies (REGs) have star formation or AGN activity, and that these active REGs have smaller axis ratio and bluer outside compared to the passive REGs. Blue early-type galaxies (BEGs) show structural features similar to those of REGs, but their centres are bluer than REGs. HII BEGs are found to have bluer centres than passive BEGs, but HII REGs have bluer outside than passive REGs. Bulge-dominated late-type galaxies have red colours. Passive red late-types are similar to REGs in several aspects. Most blue late-type galaxies (BLGs) have forming stars, but a small fraction of BLGs do not show evidence for current star formation activity. Differences of other physical parameters among different classes are inspected, and their implication on galaxy evolution is discussed.

We use a sample of galaxies from the Sloan Digital Sky Survey (SDSS) to search for correlations between the $\lambda$ spin parameter and the environment and mass of galaxies. In order to calculate the total value of $\lambda$ for each observed galaxy, we employed a simple model of the dynamical structure of the galaxies which allows a rough estimate of the value of $\lambda$ using only readily obtainable observables from the luminous galaxies. Use of a large volume limited sample (upwards of 11,000) allows reliable inferences of mean values and dispersions of $\lambda$ distributions. We find, in agreement with some N-body cosmological simulations, no significant dependence of $\lambda$ on the environmental density of the galaxies. For the case of mass, our results show a marked correlation with $\lambda$, in the sense that low mass galaxies present both higher mean values of $\lambda$ and associated dispersions, than high mass galaxies. This last direct empirical result, at odds with expectations from N-body cosmological simulations, provides interesting constrain on the mechanisms of galaxy formation and acquisition of angular momentum, a valuable test for cosmological models.

We have studied the radial distribution of the early (E/S0) and late (S/Irr) types of satellites around bright host galaxies. We made a volume-limited sample of 4,986 satellites brighter than M_r = -18.0 associated with 2,254 hosts brighter than M_r =-19.0 from the SDSS DR5 sample. The morphology of satellites is determined by an automated morphology classifier, but the host galaxies are visually classified. We found segregation of satellite morphology as a function of the projected distance from the host galaxy. The amplitude and shape of the early-type satellite fraction profile are found to depend on the host luminosity. This is the morphology-radius/density relation at the galactic scale. There is a strong tendency for morphology conformity between the host galaxy and its satellites. The early-type fraction of satellites hosted by early-type galaxies is systematically larger than that of late-type hosts, and is a strong function of the distance from the host galaxies. Fainter satellites are more vulnerable to the morphology transformation effects of hosts. Dependence of satellite morphology on the large-scale background density was detected. The fraction of early-type satellites increases in high density regions for both early and late-type hosts. It is argued that the conformity in morphology of galactic satellite system is mainly originated by the hydrodynamical and radiative effects of hosts on satellites.

We present results from the first phase of the Seoul National University Bright Quasar Survey in Optical (SNUQSO) as well as its basic observational setup. Previous and current large-area surveys have been successful in identifying many quasars, but they could have missed bright quasars due to their survey design. In order to help complete the census of bright quasars, we have performed spectroscopic observations of new bright quasar candidates selected from various methods based on optical colors, near-infrared colors, radio, and X-ray data. In 2005/2006, we observed 55 bright quasar candidates using the Bohyunsan Optical Echelle Spectrograph (BOES) on the 1.8 m telescope at the Bohyunsan Optical Astronomy Observatory in Korea. We identify 14 quasars/Seyferts from our observation, including an optically bright quasar with i=14.98 mag at z=0.092 (SDSS J003236.59-091026.2). Non-quasar/Seyfert objects are found to be mostly stars, among which there are five M-type stars and one cataclysmic variable. Our result shows that there still exist bright quasars to be discovered. However, at the same time, we conclude that finding new bright quasars in high Galactic latitude regions is very challenging and that the existing compilation of optically bright quasars is nearly complete in the northern hemisphere.

The general world model for homogeneous and isotropic universe has been roposed. For this purpose, we introduce a global and fiducial system of reference (world reference frame) constructed on a 5-dimensional space-time that is embedding the universe, and define the line element as the separation between two neighboring events that are distinct in space and time, as viewed in the world reference frame. The effect of cosmic expansion on the measurement of physical distance has been correctly included in the new metric, which differs from the Friedmann-Robertson-Walker metric where the spatial separation is measured for events on the hypersurface at a constant time while the temporal separation is measured for events at different time epochs. The Einstein's field equations with the new metric imply that closed, flat, and open universes are filled with positive, zero, and negative energy, respectively. We have demonstrated that the flat universe is empty and stationary, equivalent to the Minkowski space-time, and that the universe with positive energy density is always spatially closed and finite. In the closed universe, the proper time of a comoving observer does not elapse uniformly as judged in the world reference frame, in which both cosmic expansion and time-varying light speeds cannot exceed the limiting speed of the special relativity. We have also reconstructed cosmic evolution histories of the closed world models that are consistent with recent astronomical observations, and derived useful formulas such as energy-momentum relation of particles, redshift, total energy in the universe, cosmic distance and time scales, and so forth. It has also been shown that the inflation with positive acceleration at the earliest epoch is improbable.

We propose a model of allocating galaxies in cosmological N-body simulations. we identify each subhalo with a galaxy, and assign luminosity and morphological type assuming that the galaxy luminosity is a monotonic function of its host subhalo mass. The morphology assignment is made by using two simple relations between subhalo mass and galaxy luminosity of different types. One is using a constant ratio in luminosity of early (E/SO) and late (S/Irr) type galaxies at a fixed subhalo mass. And the other assumes that galaxies of different morphological types but having an equal luminosity have a constant ratio in their subhalo masses. We made a series of comparisons of the properties of these simulated galaxies with those of the SDSS galaxies. The resulting simulated galaxy sample is found to successfully reproduce the observed local number density distribution except for in high density regions. The luminosity function is studied as a function of local density. It was found that the observed luminosity functions in different local density environments are overall well-reproduced by the simulated galaxies. Discrepancy is found at the bright end of the luminosity function of early types in the underdense regions and at the faint end of both morphological types in very high density regions. A significant fraction of the observed early type galaxies in voids seems to have undergone a relatively recent star formation and became brighter. The lack of faint simulated galaxies in dense regions may be due to the strong tidal force of the central halo which destroys less massive satellite subhalos around in the simulation. The mass-to-light ratio is found to depend on the local density in the way similar to that observed in the SDSS sample. (abridged)

We present a unified picture on the evolution of galaxy luminosity and morphology. Galaxy morphology is found to depend critically on the local environment set up by the nearest neighbor galaxy in addition to luminosity and the large scale density. When a galaxy is located farther than the virial radius from its closest neighbor, the probability for the galaxy to have an early morphological type is an increasing function only of luminosity and the local density due to the nearest neighbor ($\rho_n$). The tide produced by the nearest neighbor is thought to be responsible for the morphology transformation toward the early type at these separations. When the separation is less than the virial radius, i.e. when $\rho_n > \rho_{\rm virial}$, its morphology depends also on the neighbor's morphology and the large-scale background density over a few Mpc scales ($\rho_{20}$) in addition to luminosity and $\rho_n$. The early type probability keeps increasing as $\rho_n$ increases if its neighbor is an early type. But the probability decreases as $\rho_n$ increases when the neighbor is a late type. The cold gas streaming from the late type neighbor can be the reason for the morphology transformation toward late type. The overall early-type fraction increases as $\rho_{20}$ increases when $\rho_n > \rho_{\rm virial}$. This can be attributed to the hot halo gas of the neighbor which is confined by the pressure of the ambient medium held by the background mass. We have also found that galaxy luminosity depends on $\rho_n$, and that the isolated bright galaxies are more likely to be recent merger products. We propose a scenario that a series of morphology and luminosity transformation occur through distant interactions and mergers, which results in the morphology--luminosity--local density relation.

We find an exact coordinate transformation rule from the $AdS_5$ Schwarzschild black hole in the Poincare and the global patch to the Fefferman-Graham coordinate system. Using these results, we evaluate the corresponding holographic stress tensor and trace anomaly of the boundary theory as a function of the radial coordinate. Following the AdS/CFT correspondence, we reinterpret the radial coordinate dependence of the trace anomaly as the Wilsonian renormalization group(RG) flow of the boundary theory.

We present a study of the connection between star-forming galaxies, AGN host galaxies, and normal early-type galaxies in the Sloan Digital Sky Survey (SDSS). Using the SDSS DR5 and DR4plus data, we select our early-type galaxy sample in the color versus color-gradient space, and we classify the spectral types of the selected early-type galaxies into normal, star-forming, Seyfert, and LINER classes, using several spectral line flux ratios. We investigate the slope in the fundamental space for each class of early-type galaxies and find that there are obvious differences in the slopes of the fundamental planes (FPs) among the different classes of early-type galaxies, in the sense that the slopes for Seyferts and star-forming galaxies are flatter than those for normal galaxies and LINERs. This may be the first identification of the systematic variation of the FP slope among the subclasses of early-type galaxies. The difference in the FP slope might be caused by the difference in the degree of nonhomology among different classes or by the difference of gas contents in their merging progenitors. One possible scenario is that the AGN host galaxies and star-forming galaxies are formed by gas-rich merging and that they may evolve into normal early-type galaxies after finishing their star formation or AGN activities.

We have analyzed the HI aperture synthesis image of the Large Magellanic Cloud (LMC), using an objective and quantitative measure of topology to understand the HI distribution hosting a number of holes and clumps of various sizes in the medium. The HI distribution shows different topology at four different chosen scales. At the smallest scales explored (19-29 pc), the HI mass is distributed in such a way that numerous clumps are embedded on top of a low density background. At the larger scales from 73 to 194 pc, it shows a generic hole topology. These holes might have been formed mainly by stellar winds from hot stars. At the scales from 240 to 340 pc, slightly above the disk scale-height of the gaseous disk, major clumps in the HI map change the distribution to have a slight clump topology. These clumps include the giant cloud associations in the spiral arms and the thick filaments surrounding superholes. At the largest scales studied (390-485 pc), the hole topology is present again. Responsible to the hole topology at this scale are a few superholes which seem mainly associated with supernova explosions in the outer disk. The gaps between the bar and spiral arms have a minor effect on the topology at this scale.

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)

We investigate the dependence of physical properties of galaxies brighter than M_r=-18.0 in the SDSS on environment, as measured by local density using an adaptive smoothing kernel. We find that variations of galaxy properties with environment are almost entirely due to the dependence of morphology and luminosity on environment. Because galaxy properties depend not only on luminosity but also on morphology, it is clear that galaxy properties cannot be determined solely by dark halo mass. When morphology and luminosity are fixed, other physical properties, such as color, color-gradient, concentration, size, velocity dispersion, and star formation rate, are nearly independent of local density. The only feature is the sharp decrease of late type fraction above the critical luminosity of about M_r=-21.3 in the morphology versus luminosity relation. Weak residual dependences on environment include that of the color of late-types (bluer at lower density) and of the L-sigma relation of early-types (larger dispersion at higher density for bright galaxies). The early-type fraction is a monotonically increasing function of local density and luminosity. The morphology-density- luminosity relation should be a key constraint on galaxy formation models. We demonstrate that the dependence on environment of the morphology of galaxie originates from variations in density on effective Gaussian smoothing scales much smaller than 12 h^{-1}Mpc. We find that galaxy morphology varies both with density measured on an effective Gaussian smoothing scale of 4.7 h^{-1} Mpc and on distance to the nearest bright galaxy (particularly, for about 0.2 h^{-1}Mpc). We propose a mechanism that the morphology of galaxies in galaxy systems is transformed by the tidal force.

Using simple dimensional arguments for both spiral and elliptical galaxies, we present formulas to derive an estimate of the halo spin parameter $\lambda$ for any real galaxy, in terms of common observational parameters. This allows a rough estimate of $\lambda$, which we apply to a large volume limited sample of galaxies taken from the SDSS data base. The large numbers involved (11,597) allow the derivation of reliable $\lambda$ distributions, as signal adds up significantly in spite of the errors in the inferences for particular galaxies. We find that if the observed distribution of $\lambda$ is modeled with a log-normal function, as often done for this distribution in dark matter halos that appear in cosmological simulations, we obtain parameters $\lambda_{0}=0.04 \pm 0.005$ and $\sigma_{\lambda}=0.51 \pm 0.05$, interestingly consistent with values derived from simulations. For spirals, we find a good correlation between empirical values of $\lambda$ and visually assigned Hubble types, highlighting the potential of this physical parameter as an objective classification tool.

We measure the large-scale real-space power spectrum P(k) using luminous red galaxies (LRGs) in the Sloan Digital Sky Survey (SDSS) and use this measurement to sharpen constraints on cosmological parameters from the Wilkinson Microwave Anisotropy Probe (WMAP). We employ a matrix-based power spectrum estimation method using Pseudo-Karhunen-Loeve eigenmodes, producing uncorrelated minimum-variance measurements in 20 k-bands of both the clustering power and its anisotropy due to redshift-space distortions, with narrow and well-behaved window functions in the range 0.01h/Mpc < k < 0.2h/Mpc. Results from the LRG and main galaxy samples are consistent, with the former providing higher signal-to-noise. Our results are robust to omitting angular and radial density fluctuations and are consistent between different parts of the sky. They provide a striking confirmation of the predicted large-scale LCDM power spectrum. Combining only SDSS LRG and WMAP data places robust constraints on many cosmological parameters that complement prior analyses of multiple data sets. The LRGs provide independent cross-checks on Om and the baryon fraction in good agreement with WMAP. Within the context of flat LCDM models, our LRG measurements complement WMAP by sharpening the constraints on the matter density, the neutrino density and the tensor amplitude by about a factor of two, giving Omega_m=0.24+-0.02 (1 sigma), sum m_nu < 0.9 eV (95%) and r<0.3 (95%). Baryon oscillations are clearly detected and provide a robust measurement of the comoving distance to the median survey redshift z=0.35 independent of curvature and dark energy properties. Within the LCDM framework, our power spectrum measurement improves the evidence for spatial flatness, sharpening the curvature constraint Omega_tot=1.05+-0.05 from WMAP alone to Omega_tot=1.003+-0.010. Assuming Omega_tot=1, the equation of state parameter is constrained to w=-0.94+-0.09, indicating the potential for more ambitious future LRG measurements to provide precision tests of the nature of dark energy. All these constraints are essentially independent of scales k>0.1h/Mpc and associated nonlinear complications, yet agree well with more aggressive published analyses where nonlinear modeling is crucial.

We have independently measured the genus topology of the temperature fluctuations in the cosmic microwave background seen in the Wilkinson Microwave Anisotropy Probe (WMAP) 3-year data. A genus analysis of the WMAP data indicates consistency with Gaussian random-phase initial conditions, as predicted by standard inflation. We set 95% confidence limits on non-linearities of -101 < f_{nl} < 107. We also find that the observed low l (l <= 8) modes show a slight anti-correlation with the Galactic foreground, but not exceeding 95% confidence, and that the topology defined by these modes is consistent with that of a Gaussian random-phase distribution (within 95% confidence).

We measure the three-dimensional topology of large-scale structure in the Sloan Digital Sky Survey (SDSS). This allows the genus statistic to be measured with unprecedented statistical accuracy. The sample size is now sufficiently large to allow the topology to be an important tool for testing galaxy formation models. For comparison, we make mock SDSS samples using several state-of-the-art N-body simulations: the Millennium run of Springel et al. (2005)(10 billion particles), Kim & Park (2006) CDM models (1.1 billion particles), and Cen & Ostriker (2006) hydrodynamic code models (8.6 billion cell hydro mesh). Each of these simulations uses a different method for modeling galaxy formation. The SDSS data show a genus curve that is broadly characteristic of that produced by Gaussian random phase initial conditions. Thus the data strongly support the standard model of inflation where Gaussian random phase initial conditions are produced by random quantum fluctuations in the early universe. But on top of this general shape there are measurable differences produced by non-linear gravitational effects (cf. Matsubara 1994), and biasing connected with galaxy formation. The N-body simulations have been tuned to reproduce the power spectrum and multiplicity function but not topology, so topology is an acid test for these models. The data show a ``meatball'' shift (only partly due to the Sloan Great Wall of Galaxies; this shift also appears in a sub-sample not containing the Wall) which differs at the 2.5\sigma level from the results of the Millennium run and the Kim & Park dark halo models, even including the effects of cosmic variance.

We have independently measured the genus topology of the temperature fluctuations in the cosmic microwave background seen in the Wilkinson Microwave Anisotropy Probe (WMAP) 3-year data. A genus analysis of the WMAP data indicates consistency with Gaussian random-phase initial conditions, as predicted by standard inflation. We set 95% confidence limits on non-linearities of -101 < f_{nl} < 107. We also find that the observed low l (l <= 8) modes show a slight anti-correlation with the Galactic foreground, but not exceeding 95% confidence, and that the topology defined by these modes is consistent with that of a Gaussian random-phase distribution (within 95% confidence).

[Abridged] We have produced a cleaned map of the Wilkinson Microwave Anisotropy Probe (WMAP) 3-year data using an improved foreground subtraction technique. We perform an internal linear combination (ILC) to subtract the Galactic foreground emission from the temperature fluctuations observed by the WMAP. We divide the whole sky into hundreds of pixel groups with similar foreground spectral indices over a range of WMAP frequencies, apply the ILC for each group, and obtain a CMB map with foreground emission effectively reduced. With the resulting foreground-reduced ILC map, we have investigated the known anomalies in CMB maps at large scales. The quadrupole and the octopole powers measured from our ILC map are \delta T_2^2 = 276_{-126}^{+94} uK^2 and \delta T_3^2 = 952_{-83}^{+64} uK^2, respectively. The measured quadrupole power is lower than the value expected in the concordance LCDM model (1250 uK^2), in which the probability of finding a quadrupole power lower than the measured value is 5.7%. We have confirmed that the quadrupole and the octopole are strongly aligned with angle \theta_{23} = 11.8_{-8.0}^{+6.4} degree, and are planar with high planarity parameters t=0.98_{-0.02}^{+0.02} for l=2 and t=0.91_{-0.03}^{+0.02} for l=3. The observed angular separation \theta_{23} is marginally statistically significant because the probability of finding the angular separation as low as the observed value is 4.3%. However, the observed planarity is not statistically significant. The probability of observing such a planarity as high as the measured t values is over 18%. The ILC simulations show that the residual foreground emission in the ILC map does not affect the estimated values significantly.

[Abridged] We have produced a cleaned map of the Wilkinson Microwave Anisotropy Probe (WMAP) 3-year data using an improved foreground subtraction technique. We perform an internal linear combination (ILC) to subtract the Galactic foreground emission from the temperature fluctuations observed by the WMAP. We divide the whole sky into hundreds of pixel groups with similar foreground spectral indices over a range of WMAP frequencies, apply the ILC for each group, and obtain a CMB map with foreground emission effectively reduced. With the resulting foreground-reduced ILC map, we have investigated the known anomalies in CMB maps at large scales. The quadrupole and the octopole powers measured from our ILC map are \delta T_2^2 = 276_{-126}^{+94} uK^2 and \delta T_3^2 = 952_{-83}^{+64} uK^2, respectively. The measured quadrupole power is lower than the value expected in the concordance LCDM model (1250 uK^2), in which the probability of finding a quadrupole power lower than the measured value is 5.7%. We have confirmed that the quadrupole and the octopole are strongly aligned with angle \theta_{23} = 11.8_{-8.0}^{+6.4} degree, and are planar with high planarity parameters t=0.98_{-0.02}^{+0.02} for l=2 and t=0.91_{-0.03}^{+0.02} for l=3. The observed angular separation \theta_{23} is marginally statistically significant because the probability of finding the angular separation as low as the observed value is 4.3%. However, the observed planarity is not statistically significant. The probability of observing such a planarity as high as the measured t values is over 18%. The ILC simulations show that the residual foreground emission in the ILC map does not affect the estimated values significantly.

We have found the u -r color versus g -i color gradient space can be used for highly successful morphology classification of galaxies in the Sloan Digital Sky Survey. In this space galaxies form early and late type branches well-separated from each other. The location of galaxies along the branches reflects the degree and locality of star formation activity, and monotonically corresponds to the sequence of morphological subclasses. When the concentration index is used together, the completeness and reliability of classification reaches about 91% for a training set of SDSS galaxies brighter than r~15.9. At faintest magnitudes (r~17.5) of the SDSS spectroscopic sample, the performance still remains at about 88%. The new classification scheme will help us find accurate relations of galaxy morphology with spatial and temporal environments, and help one understand the origin of morphology of galaxies.

We have measured the cosmic momentum power spectrum from the peculiar velocities of galaxies in the SFI sample. The SFI catalog contains field spiral galaxies with radial peculiar velocities derived from the I-band Tully-Fisher relation. As a natural measure of the large-scale peculiar velocity field, we use the cosmic momentum field that is defined as the peculiar velocity field weighted by local number of galaxies. We have shown that the momentum power spectrum can be derived from the density power spectrum for the constant linear biasing of galaxy formation, which makes it possible to estimate \beta_S = \Omega_m^{0.6} / b_S parameter precisely where \Omega_m is the matter density parameter and b_S is the bias factor for optical spiral galaxies. At each wavenumber k we estimate \beta_S(k) as the ratio of the measured to the derived momentum power over a wide range of scales (0.026 h^{-1}Mpc <~ k <~ 0.157 h^{-1}Mpc) that spans the linear to the quasi-linear regimes. The estimated \beta_S(k)'s have stable values around 0.5, which demonstrates the constancy of \beta_S parameter at scales down to 40 h^{-1}Mpc. We have obtained \beta_S=0.49_{-0.05}^{+0.08} or \Omega_m = 0.30_{-0.05}^{+0.09} b_S^{5/3}, and the amplitude of mass fluctuation as \sigma_8\Omega_m^{0.6}=0.56_{-0.21}^{+0.27}. The 68% confidence limits include the cosmic variance. We have also estimated the mass density power spectrum. For example, at k=0.1047 h Mpc^{-1} (\lambda=60 h^{-1}Mpc) we measure \Omega_m^{1.2} P_{\delta}(k)=(2.51_{-0.94}^{+0.91})\times 10^3 (h^{-1}Mpc)^3, which is lower compared to the high-amplitude power spectra found from the previous maximum likelihood analyses of peculiar velocity samples like Mark III, SFI, and ENEAR.

We have measured the cosmic momentum power spectrum from the peculiar velocities of galaxies in the SFI sample. The SFI catalog contains field spiral galaxies with radial peculiar velocities derived from the I-band Tully-Fisher relation. As a natural measure of the large-scale peculiar velocity field, we use the cosmic momentum field that is defined as the peculiar velocity field weighted by local number of galaxies. We have shown that the momentum power spectrum can be derived from the density power spectrum for the constant linear biasing of galaxy formation, which makes it possible to estimate \beta_S = \Omega_m^{0.6} / b_S parameter precisely where \Omega_m is the matter density parameter and b_S is the bias factor for optical spiral galaxies. At each wavenumber k we estimate \beta_S(k) as the ratio of the measured to the derived momentum power over a wide range of scales (0.026 h^{-1}Mpc <~ k <~ 0.157 h^{-1}Mpc) that spans the linear to the quasi-linear regimes. The estimated \beta_S(k)'s have stable values around 0.5, which demonstrates the constancy of \beta_S parameter at scales down to 40 h^{-1}Mpc. We have obtained \beta_S=0.49_{-0.05}^{+0.08} or \Omega_m = 0.30_{-0.05}^{+0.09} b_S^{5/3}, and the amplitude of mass fluctuation as \sigma_8\Omega_m^{0.6}=0.56_{-0.21}^{+0.27}. The 68% confidence limits include the cosmic variance. We have also estimated the mass density power spectrum. For example, at k=0.1047 h Mpc^{-1} (\lambda=60 h^{-1}Mpc) we measure \Omega_m^{1.2} P_{\delta}(k)=(2.51_{-0.94}^{+0.91})\times 10^3 (h^{-1}Mpc)^3, which is lower compared to the high-amplitude power spectra found from the previous maximum likelihood analyses of peculiar velocity samples like Mark III, SFI, and ENEAR.

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.

We measure the topology of volume-limited galaxy samples selected from a parent sample of 314,050 galaxies in the Sloan Digital Sky Survey(SDSS), which is now complete enough to describe the fully three-dimensional topology and its dependence on galaxy properties. We compare the observed genus statistic G\nu_f) to predictions for a Gaussian random field and to the genus measured for mock surveys constructed from new large-volume simulations of the LambdaCDM cosmology. In this analysis we carefully examine the dependence of the observed genus statistic on the Gaussian smoothing scale R_G from 3.5 to 11 h^-1Mpc and on the luminosity of galaxies over the range -22.50<M_r<-18.5. The void multiplicity A_V is less than unity at all smoothing scales. Because A_V cannot become less than 1 through gravitational evolution, this result provides strong evidence for biased galaxy formation in low density environments. We also find clear evidence of luminosity bias of topology within the volume-limited sub-samples. The shift parameter Delta nu indicates that the genus of brighter galaxies shows a negative shift toward a ``meatball'' (i.e. cluster-dominated) topology, while faint galaxies show a positive shift toward a ``bubble'' (i.e. void-dominated) topology. The transition from negative to positive shift occurs approximately at the characteristic absolute magnitude M_{r*}=-20.4. Even in this analysis of the largest galaxy sample to date, we detect the influence of individual large-scale structures, as the shift parameter Delta nu and cluster multiplicity A_C reflect (at ~3 sigma) the presence of the Sloan Great Wall and a X-shaped structure which runs for several hundred Mpc across the survey volume.

We have studied the dependence of topology of large scale structure on tracer, gravitational evolution, redshift space distortion, and cosmology. A series of large N-body simulations of the $\Lambda$CDM and SCDM models that have evolved 1.1 or 8.6 billion particles, are used in the study. Evolution of the genus statistic, used as a topology measure, from redshift 8 to 0 is accurately calculated over a wide range of smoothing scales using the simulations. The tracers of large scale structure considered are the CDM matter, biased peaks in the initial density field, dark halos, and `galaxies' populating the dark halos in accordance with a Halo Occupation Distribution (HOD) model. We have found that the effects of biasing, gravitational evolution, and initial conditions on topology of large scale structure are all comparable. The redshift space distortion effects are relatively small down to about 5 $h^{-1}$Mpc for all tracers except for the high threshold part of the genus curve. The gravitational effects are found to be well-modeled by the analytic perturbation theory when the CDM matter distribution is considered. But the direction of gravitational evolution of topology can be even reversed for different tracers. For example, the shift parameter of the genus curve evolves in opposite directions for matter and HOD `galaxies' at large scales. At small scales there are interesting deviations of the genus curve of dark halos and `galaxies' from that of matter in our initially Gaussian simulations. The deviations should be understood as due to combined effects of gravitational evolution and biasing.

R Aqr is known to be a symbiotic binary system with an associated extended emission nebula, possibly produced by a historic outburst. To find the associated historic records, we searched for and compiled all 'Guest Star' and 'Peculiar Star' records in three Korean 'official' history books that cover almost two thousand years, Samguksagi, Goryeosa, Joseonwangjosillok. In addition to the record of A.D. 1073, previously noted by Li (1985), we have found in Goryeosa another candidate record of A.D. 1074, which has the same positional description as that of A.D. 1073 with an additional brightness description. We examined various aspects of the two records and conclude that they both are likely to be the records of outburst of R Aqr. This means that there were two successive outbursts in A.D. 1073 and in A.D. 1074, separated by approximately one year. Based on these records, we estimate the distance to R Aqr to be 273 pc if the expansion of the nebula has been at a constant rate. The brightness record of A.D. 1074 corresponds to the absolute magnitude at outburst of M_(outburst) = -6.2 mag. ~ -5.2 mag. at this distance. The two Korean records associated with R Aqr may provide astronomically meaningful constraints to the outburst model of R Aqr and the formative process of its nebulosity.

We have compiled and analyzed historical Korean meteor and meteor shower records in three Korean official history books, Samguksagi which covers the three Kingdoms period (57 B.C -- A.D. 935), Goryeosa of Goryeo dynasty (A.D. 918 -- 1392), and Joseonwangjosillok of Joseon dynasty (A.D. 1392 -- 1910). We have found 3861 meteor and 31 meteor shower records. We have confirmed the peaks of Perseids and an excess due to the mixture of Orionids, north-Taurids, or Leonids through the Monte-Carlo test. The peaks persist from the period of Goryeo dynasty to that of Joseon dynasty, for almost one thousand years. Korean records show a decrease of Perseids activity and an increase of Orionids/north-Taurids/Leonids activity. We have also analyzed seasonal variation of sporadic meteors from Korean records. We confirm the seasonal variation of sporadic meteors from the records of Joseon dynasty with the maximum number of events being roughly 1.7 times the minimum. The Korean records are compared with Chinese and Japanese records for the same periods. Major features in Chinese meteor shower records are quite consistent with those of Korean records, particularly for the last millennium. Japanese records also show Perseids feature and Orionids/north-Taurids/Leonids feature, although they are less prominent compared to those of Korean or Chinese records.

We present results from a test for the Gaussianity of the whole sky sub-degree scale CMB temperature anisotropy measured by the Wilkinson Microwave Anisotropy Probe (WMAP). We calculate the genus from the foreground-subtracted and Kp0-masked WMAP maps and measure the genus shift parameters defined at negative and positive threshold levels and the asymmetry parameter to quantify the deviation from the Gaussian relation. At WMAP Q, V, and W bands, the genus and genus-related statistics imply that the observed CMB sky is consistent with Gaussian random phase field. However, from the genus measurement on the Galactic northern and southern hemispheres, we have found two non-Gaussian signatures at the W band resolution (0.35 degree scale), i.e., the large difference of genus amplitudes between the north and the south and the positive genus asymmetry in the south, which are statistically significant at 2.6 sigma and 2.4 sigma levels, respectively. The large genus amplitude difference also appears in the WMAP Q and V band maps, deviating the Gaussian prediction with a significance level of about 2 sigma. The probability that the genus curves show such a large genus amplitude difference exceeding the observed values at all Q, V, and W bands in a Gaussian sky is only 1.4%. Such non-Gaussian features are reduced as the higher Galactic cut is applied, but their dependence on the Galactic cut is weak. We discuss possible sources that can induce such non-Gaussian features, and conclude that the CMB data with higher signal-to-noise ratio and the accurate foreground model are needed to understand the non-Gaussian signatures.

We study the problem of separating E and B modes in interferometric observations of the polarization of the cosmic microwave background. The E and B band powers and their mixings are measured from both single-dish and interferometric mock observations using the quadratic estimator of the maximum likelihood analysis. We find that the interferometer can separate E and B modes in a single-pointing measurement and is thus well suited for detecting the faint lensing induced and gravity-wave induced B modes. In mosaicking observation, compared to the single dish, the interferometer is in general more efficient in separating E and B modes, and for high signal-to-noise per pixel it needs about three times fewer pixels to measure extremely blue polarization power spectra.

We have developed a new halo finding method, Physically Self-Bound (PSB) group finding algorithm, which can efficiently identify halos located even at crowded regions. This method combines two physical criteria such as the tidal radius of a halo and the total energy of each particle to find member particles. Two hierarchical meshes are used to increase the speed and the power of halo identification in the parallel computing environments. First, a coarse mesh with cell size equal to the mean particle separation $l_{\rm mean}$ is used to obtain the density field over the whole simulation box. Mesh cells having density contrast higher than a local cutoff threshold $\delta_{\rm LOC}$ are extracted and linked together for those adjacent to each other. This produces local-cell groups. Second, a finer mesh is used to obtain density field within each local-cell group and to identify halos. If a density shell contains only one density peak, its particles are assigned to the density peak. But in the case of a density shell surrounding at least two density peaks, we use both the tidal radii of halo candidates enclosed by the shell and the total energy criterion to find physically bound particles with respect to each halo. Similar to DENMAX and HOP, the \hfind method can efficiently identify small halos embedded in a large halo, while the FoF and the SO do not resolve such small halos. We apply our new halo finding method to a 1-Giga particle simulation of the $\Lambda$CDM model and compare the resulting mass function with those of previous studies. The abundance of physically self-bound halos is larger at the low mass scale and smaller at the high mass scale than proposed by the Jenkins et al. (2001) who used the FoF and SO methods. (abridged)

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, $\theta \to 0$. (abridged)

We have simulated the interferometric observation of the Cosmic Microwave Background (CMB) temperature and polarization fluctuations. We have constructed data pipelines from the time-ordered raw visibility samples to the CMB power spectra which utilize the methods of data compression, maximum likelihood analysis, and optimal subspace filtering. They are customized for three observational strategies, such as the single pointing, the mosaicking, and the drift-scanning. For each strategy, derived are the optimal strategy parameters that yield band power estimates with minimum uncertainty. The results are general and can be applied to any close-packed array on a single platform such as the CBI and the forthcoming AMiBA experiments. We have also studied the effect of rotation of the array platform on the band power correlation by simulating the CBI single pointing observation. It is found that the band power anti-correlations can be reduced by rotating the platform and thus densely sampling the visibility plane. This enables us to increase the resolution of the power spectrum in the l-space down to the limit of the sampling theorem (Delta l = 226 = pi / theta), which is narrower by a factor of about sqrt{2} than the resolution limit (Delta l = 300) used in the recent CBI single pointing observation. The validity of this idea is demonstrated for a two-element interferometer that samples visibilities uniformly in the uv-annulus. From the fact that the visibilities are the Fourier modes of the CMB field convolved with the beam, a fast unbiased estimator (FUE) of the CMB power spectra is developed and tested. It is shown that the FUE gives results very close to those from the quadratic estimator method without requiring large computer resources even though uncertainties in the results increase.

We have simulated the interferometric observation of the Cosmic Microwave Background (CMB) temperature and polarization fluctuations. We have constructed data pipelines from the time-ordered raw visibility samples to the CMB power spectra which utilize the methods of data compression, maximum likelihood analysis, and optimal subspace filtering. They are customized for three observational strategies, such as the single pointing, the mosaicking, and the drift-scanning. For each strategy, derived are the optimal strategy parameters that yield band power estimates with minimum uncertainty. The results are general and can be applied to any close-packed array on a single platform such as the CBI and the forthcoming AMiBA experiments. We have also studied the effect of rotation of the array platform on the band power correlation by simulating the CBI single pointing observation. It is found that the band power anti-correlations can be reduced by rotating the platform and thus densely sampling the visibility plane. This enables us to increase the resolution of the power spectrum in the l-space down to the limit of the sampling theorem (Delta l = 226 = pi / theta), which is narrower by a factor of about sqrt{2} than the resolution limit (Delta l = 300) used in the recent CBI single pointing observation. The validity of this idea is demonstrated for a two-element interferometer that samples visibilities uniformly in the uv-annulus. From the fact that the visibilities are the Fourier modes of the CMB field convolved with the beam, a fast unbiased estimator (FUE) of the CMB power spectra is developed and tested. It is shown that the FUE gives results very close to those from the quadratic estimator method without requiring large computer resources even though uncertainties in the results increase.

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.

We have made a topological study of cosmic microwave background (CMB) polarization maps by simulating the AMiBA experiment results. A $\Lambda$CDM CMB sky is adopted to make mock interferometric observations designed for the AMiBA experiment. CMB polarization fields are reconstructed from the AMiBA mock visibility data using the maximum entropy method. We have also considered effects of Galactic foregrounds on the CMB polarization fields. The genus statistic is calculated from the simulated $Q$ and $U$ polarization maps, where $Q$ and $U$ are Stokes parameters. Our study shows that the Galactic foreground emission, even at low Galactic latitude, is expected to have small effects on the CMB polarization field. Increasing survey area and integration time is essential to detect non-Gaussian signals of cosmological origin through genus measurement.

We have made a topological study of cosmic microwave background (CMB) polarization maps by simulating the AMiBA experiment results. A $\Lambda$CDM CMB sky is adopted to make mock interferometric observations designed for the AMiBA experiment. CMB polarization fields are reconstructed from the AMiBA mock visibility data using the maximum entropy method. We have also considered effects of Galactic foregrounds on the CMB polarization fields. The genus statistic is calculated from the simulated $Q$ and $U$ polarization maps, where $Q$ and $U$ are Stokes parameters. Our study shows that the Galactic foreground emission, even at low Galactic latitude, is expected to have small effects on the CMB polarization field. Increasing survey area and integration time is essential to detect non-Gaussian signals of cosmological origin through genus measurement.

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 \sigma detection limit effectively eliminates this effect. Radio sources also result in a positive genus curve asymmetry (significant at 2 \sigma) 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 \sigma level if the measured genus shift parameter |\Delta\nu| >= 0.02 (0.04) or if the measured genus asymmetry parameter |\Delta g| >= 0.03 (0.08) on a 0.3 (1.0) deg FWHM scale.

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 \sigma detection limit effectively eliminates this effect. Radio sources also result in a positive genus curve asymmetry (significant at 2 \sigma) 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 \sigma level if the measured genus shift parameter |\Delta\nu| >= 0.02 (0.04) or if the measured genus asymmetry parameter |\Delta g| >= 0.03 (0.08) on a 0.3 (1.0) deg FWHM scale.

We present results from a first test of the Gaussianity of degree-scale cosmic microwave background (CMB) anisotropy. We investigate Gaussianity of the CMB anisotropy by studying the topology of CMB anisotropy maps from the QMAP and Saskatoon experiments. We also study the QMASK map, a combination map of the QMAP and Saskatoon data. We measure the genus from noise-suppressed Wiener-filtered maps at an angular scale of about 1.5 degrees. To test the Gaussianity of the observed anisotropy, we compare these results to those derived from a collection of simulated maps for each experiment in a Gaussian spatially-flat cosmological constant dominated cold dark matter model. The genus-threshold level relations of the QMAP and Saskatoon maps are consistent with Gaussianity. While the combination QMASK map has a mildly non-Gaussian genus curve which is not a consequence of known foreground contamination, this result is not statistically significant at the 2 sigma level. These results extend previous upper limits on the non-Gaussianity of the large angular scale (> 10 degrees) CMB anisotropy (measured by the COBE DMR experiment) down to degree angular scales.

We present results from a first test of the Gaussianity of degree-scale cosmic microwave background (CMB) anisotropy. We investigate Gaussianity of the CMB anisotropy by studying the topology of CMB anisotropy maps from the QMAP and Saskatoon experiments. We also study the QMASK map, a combination map of the QMAP and Saskatoon data. We measure the genus from noise-suppressed Wiener-filtered maps at an angular scale of about 1.5 degrees. To test the Gaussianity of the observed anisotropy, we compare these results to those derived from a collection of simulated maps for each experiment in a Gaussian spatially-flat cosmological constant dominated cold dark matter model. The genus-threshold level relations of the QMAP and Saskatoon maps are consistent with Gaussianity. While the combination QMASK map has a mildly non-Gaussian genus curve which is not a consequence of known foreground contamination, this result is not statistically significant at the 2 sigma level. These results extend previous upper limits on the non-Gaussianity of the large angular scale (> 10 degrees) CMB anisotropy (measured by the COBE DMR experiment) down to degree angular scales.

We introduce the cosmic momentum field as a new measure of the large-scale peculiar velocity and matter fluctuation fields. The momentum field is defined as the peculiar velocity field traced and weighted by galaxies, and is equal to the velocity field in the linear regime. We show that the radial component of the momentum field can be considered as a scalar field with the power spectrum which is practically 1/3 of that of the total momentum field. We present a formula for the power spectrum directly calculable from the observed radial peculiar velocity data. The momentum power spectrum is measured for the MAT sample in the Mark III catalog of peculiar velocities of galaxies. Using the momentum power spectrum we find the amplitude of the matter power spectrum is 6400^{+2800}_{-1800} and 4500^{+2000}_{-1300} \Omega^{-1.2} (h^{-1} Mpc)^{-3} at the wavenumbers 0.049 and 0.074 h Mpc^{-1}, respectively, where \Omega is the density parameter. The 68% confidence limits include the cosmic variance. The measured momentum and density power spectra together indicate that the parameter \beta_O = \Omega^{0.6}/b_O = 0.51^{+0.13}_{-0.08} or \Omega = 0.33^{+0.15}_{-0.09} b_O^{5/3} where b_O is the bias factor for optical galaxies.

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.

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.

We have applied an effective numerical scheme for cosmic-ray transport to 3D MHD simulations of jet flow in radio galaxies (see the companion paper by Jones et al. 1999). The marriage of relativistic particle and 3D magnetic field information allows us to construct a rich set of ``synthetic observations'' of our simulated objects. The information is sufficient to calculate the ``true'' synchrotron emissivity at a given frequency using explicit information about the relativistic electrons. This enables us to produce synchrotron surface-brightness maps, including polarization. Inverse-Compton X-ray surface-brightness maps may also be produced. First results intended to explore the connection between jet dynamics and electron transport in radio lobes are discussed. We infer lobe magnetic field values by comparison of synthetically observed X-ray and synchrotron fluxes, and find these ``inverse-Compton'' fields to be quite consistent with the actual RMS field averaged over the lobe. The simplest minimum energy calculation from the synthetic observations also seems to agree with the actual simulated source properties.

We have found that the two-point correlation function of the APM clusters of galaxies has an amplitude much higher than that claimed by the APM group. As the richness limit increases from R = 53 to 80, the correlation length increases from 17.5 to 28.9 h^{-1}Mpc. This indicates that the richness dependence of the APM cluster correlation function is also much stronger than what the APM group has reported. The richness dependence can be represented by a fitting formula r_0 = 0.53 d_c + 0.01, which is consistent with the Bahcall's formula r_0 = 0.4 d_c. We have tried to find the possible reason for discrepancies. However, our estimates for the APM cluster correlation function are found to be robust against variation of the method of calculation and of sample definition.

In the Hubble Deep Field (HDF) we have identified candidate regions where primordial galaxies might be forming. These regions are identified from negative or positive peaks in the difference maps obtained from the HDF maps smoothed over 0.8'' and 4''. They have apparent V magnitudes typically between 29 and 31 (missing flux below the local average level for the dark objects). The identified objects are shown to be real by two ways. First, the cross-correlations of these peaks detected in different filters are strong. Second, their auto-correlation functions indicate that these faint diffuse objects are self-clustered. The subset of objects dark in the F450W and F606W bandpasses, but bright in F814W, also shows stronger correlation compared to the whole dark sample. This further supports that our samples are indeed physical objects. The amplitude and slope of the angular correlation function of the bright objects indicates that these objects are ancestors of the present nearby bright galaxies. We have inspected individual bright objects and noted that they have several tiny spots embedded in extended backgrounds. They are likely to be the primordial galaxies at high redshifts in the process of active star formation and merging. Our subset of dark objects is thought to be the `intergalactic dark clouds' he blocking the background far UV light (at the rest frame) at high redshifts instead of empty spaces between the first galaxies at the edge of the universe of galaxies.

We have simulated cosmic microwave background (CMB) anisotropy maps for several COBE-DMR-normalized cold dark matter (CDM) cosmogonies, to make predictions for the upcoming MAP experiment. We have studied the sensitivity of the simulated MAP data to cosmology, sky coverage, and instrumental noise. With accurate knowledge of instrumental noise, MAP data will discriminate among the cosmogonies considered, and determine the topology of the initial fluctuations. A correlation function analysis of the simulated MAP data results in a very accurate measurement of the acoustic Hubble radius at decoupling. A low-density open CDM model with Omega_0=0.4 can be distinguished from the Omega_0=1 fiducial CDM model or a Lambda CDM model with > 99% confidence from the location of the acoustic "valley" in the correlation function. A genus analysis of the simulated MAP data indicates that in cosmogonies with Gaussian random-phase initial conditions, a shift of the zero-crossing point of the genus curve near the mean temperature threshold level nu=0 should not exceed |Delta nu| = 0.01 (0.04) when the total effective FWHM smoothing is 0.3 deg (1.0 deg). The asymmetry of the genus curve at the positive and negative threshold levels should not exceed |Delta g/g(nu=1,-1)| = 0.8% (4%) at 0.3 deg (1.0 deg) FWHM smoothing. Deviations of the observed MAP data in excess of these small values will be evidence for non-Gaussian behavior. The amplitude of the genus curve is a measure of the shape of the power spectrum at the smoothing scale. Even with the expected amount of instrumental noise and partial sky coverage (due to the Galaxy), the MAP data should allow discrimination amongst the cosmogonies considered at more than 99% confidence solely from a genus amplitude analysis.

We have measured the topology (genus) of the fluctuations in the cosmic microwave background seen in the recently completed (four-year) data set produced by the COBE satellite. We find that the genus is consistent with that expected from a random-phase Gaussian distribution, as might be produced naturally in inflationary models.

We compute the power spectrum of galaxy density fluctuations in a recently completed redshift survey of optically-selected galaxies in the southern hemisphere (SSRS2). The amplitude and shape of the SSRS2 power spectrum are consistent with results of the Center for Astrophysics redshift survey of the northern hemisphere (CfA2), including the abrupt change of slope on a scale of 30-50Mpc/h; these results are reproducible for independent volumes of space and variations are consistent with the errors estimated from mock surveys. Taken together, the SSRS2 and CfA2 form a complete sample of 14,383 galaxies which covers one-third of the sky. The power spectrum of this larger sample continues to rise on scales up to ~ 200Mpc/h, with weak evidence for flattening on the largest scales. The SSRS2+CfA2 power spectrum and the power spectrum constraints implied by COBE are well-matched by an Omega*h ~ 0.2, Omega+lambda_0=1 CDM model with minimal biasing of optically-selected galaxies.