Gravitational waves from coalescing neutron stars encode information about nuclear matter at extreme densities, inaccessible by laboratory experiments. The late inspiral is influenced by the presence of tides, which depend on the neutron star equation of state. Neutron star mergers are expected to often produce rapidly-rotating remnant neutron stars that emit gravitational waves. These will provide clues to the extremely hot post-merger environment. This signature of nuclear matter in gravitational waves contains most information in the 2-4 kHz frequency band, which is outside of the most sensitive band of current detectors. We present the design concept and science case for a neutron star extreme matter observatory (NEMO): a gravitational-wave interferometer optimized to study nuclear physics with merging neutron stars. The concept uses high circulating laser power, quantum squeezing and a detector topology specifically designed to achieve the high-frequency sensitivity necessary to probe nuclear matter using gravitational waves. Above one kHz, the proposed strain sensitivity is comparable to full third-generation detectors at a fraction of the cost. Such sensitivity changes expected event rates for detection of post-merger remnants from approximately one per few decades with two A+ detectors to a few per year, and potentially allows for the first gravitational-wave observations of supernovae, isolated neutron stars, and other exotica.

Fast Fourier transform based estimators are formulated for measuring momentum power spectra, including the auto power spectra of the momentum, the momentum divergence, and the cross spectrum of density fluctuation and momentum divergence. Algorithms using the third order Bettle-Lemari\'{e} scaling function to assign discrete objects to regular grids for fast Fourier transform are proposed to clean alias effects. Numerical experiments prove that the implementation can achieve sub-percent precision till close to the Nyquist frequency. Impact of removing bulk flow on estimation of momentum power spectra are derived theoretically and verified numerically, subtracting bulk flow has little effects at large scales but might induce meaningful differences in nonlinear regime, and probably it is not necessary to subtract bulk flow for samples which peculiar velocities are exact or sufficiently accurate. Momentum power spectra of dark matter samples from N-body simulation are measured and discussed. As expected, prediction of the one loop Eulerian perturbation theory agrees with simulation only slightly better than the linear theory at z=0, but can be applied to higher redshift with improved accuracy. Measurements of simulation data and the one loop Eulerian theory both reveal that the momentum field contains strong rotational part, and there is a large stochastic component in the divergence of momentum which is not correlated with the density field. The three kinds of momentum power spectra have their own characteristics.

Kilonovae are optical flashes produced in the aftermath of neutron star-neutron star mergers (NNMs) or neutron star-black hole mergers (NBMs). In this work, we use the Millennium Simulation, combined with a semi-analytic galaxy formation model--GABE (Galaxy Assembly with Binary Evolution) which adopts binary stellar population synthesis models, to explore the cosmic event rate of kilonovae, and the properties of their host galaxies in a cosmological context. We find that model with supernova kick velocity of 0 km/s fits the observation best, in agreement with the exception of some formation channels of binary neutron star. This indicates that NNMs prefer to originate from binary systems with low kick velocities. With V$_{\rm kick}$=0 km/s, the cosmic event rate of NNMs and NBMs at z=0 are 283 Gpc$^{-3}$yr$^{-1}$ and 91 Gpc$^{-3}$yr$^{-1}$, respectively, marginally consistent with the constraint from LVC GWTC-1. For Milky Way-mass galaxies, we predict the NNM rate is $25.7^{+59.6}_{-7.1}$ Myr$^{-1}$, which is also in good agreement with the observed properties of binary neutron stars in the Milky Way. Taking all the NNMs into account in the history of Milky Way-mass galaxies, we find that the averaged r-process elements yield with A>79 in a NNM and NBM event should be 0.01 M$_{\odot}$ to be consistent with observation. We conclude that NGC 4993, the host galaxy of GW170817, is a typical host galaxy for NNMs. However, generally NNMs and NBMs tend to reside in young, blue, star-forming, late-type galaxies, with stellar mass and gaseous metallicity distribution peaking at $10^{10.65}$ M$_{\odot}$ and 8.72-8.85, respectively. By studying kilonovae host galaxies in the cosmological background, it is promising to constrain model details better when we have more events in the forthcoming future. (abridged)

We report an experiment to test quantum interference, entanglement and nonlocality using two dissimilar photon sources, the Sun and a semiconductor quantum dot on the Earth, which are separated by 150 million kilometers. By making the otherwise vastly distinct photons indistinguishable all degrees of freedom, we observe time-resolved two-photon quantum interference with a raw visibility of 0.796(17), well above the 0.5 classical limit, providing the first evidence of quantum nature of thermal light. Further, using the photons with no common history, we demonstrate post-selected two-photon entanglement with a state fidelity of 0.826(24), and a violation of Bell's inequality by 2.20(6). The experiment can be further extended to a larger scale using photons from distant stars, and open a new route to quantum optics experiments at an astronomical scale.

We developed a new semi-analytic galaxy formation model: Galaxy Assembly with Binary Evolution (GABE). For the first time we introduce binary evolution into semi-analytic models of galaxy formation by using Yunnan-II stellar population synthesis model, which includes various binary interactions. When implementing our galaxy formation model onto the merger trees extracted from the Millennium simulation, it can reproduce a large body of observational results. We find that in the local universe the model including binary evolution reduces the luminosity at optical and infrared wavelengths slightly, while increases the luminosity at ultraviolet wavelength significantly, especially in $F_{\rm UV}$ band. The resulting luminosity function does not change very much over SDSS optical bands and infrared band, but the predicted colors are bluer, especially when $F_{\rm UV}$ band is under consideration. The new model allows us to explore the physics of various high energy events related to the remnants of binary stars, e.g. type Ia supernovae, short gamma-ray bursts and gravitational wave events, and their relation with host galaxies in a cosmological context.

We compare galaxy morphology predicted by the Illustris hydrodynamical simulation and a semi-analytic model (SAM) grafted in the halo merger trees from the Illustris-Dark matter simulation. Morphology is classified according to the luminous profile and the kinematic bulge-to-total ratio for Illustris galaxy, and the bulge-to-total stellar mass ratio for SAM galaxy. For late-type galaxies in the Illustris catalogue, most of their counterparts in the SAM model have the same type, and the consistency between two models is higher for lower mass galaxies. For early-type (ET) galaxies in Illustris, the consistency is quite low for the counterparts except for most massive galaxies. By comparing in detail the growth histories of some matched galaxy pairs of Milky Way mass in Illustris and the SAM model, we notice two aspects of differences in determining galaxy morphology between the two galaxy formation implementations. First, in the SAM, major merger and frequent minor mergers result in the growth of bulges and turn the galaxy into ET, while bulge formation is not connected to mergers as tightly as in SAM for the Illustris galaxies. In addition, the satellite stellar mass can decrease significantly due to tidal stripping before merging into the central galaxy in Illustris, while it does not decrease in the SAM model. This results in less mergers with large (stellar) mass ratios in the Illustris simulation, and less effect of mergers on shaping galaxy morphology.

We use a volume-limited galaxy sample from the SDSS Data Release 7 to explore the dependence of galactic conformity on the large-scale environment, measured on $\sim$ 4 Mpc scales. We find that the star formation activity of neighbour galaxies depends more strongly on the environment than on the activity of their primary galaxies. In under-dense regions most neighbour galaxies tend to be active, while in over-dense regions neighbour galaxies are mostly passive, regardless of the activity of their primary galaxies. At a given stellar mass, passive primary galaxies reside in higher density regions than active primary galaxies, leading to the apparently strong conformity signal. The dependence of the activity of neighbour galaxies on environment can be explained by the corresponding dependence of the fraction of satellite galaxies. Similar results are found for galaxies in a semi-analytical model, suggesting that no new physics is required to explain the observed large-scale conformity.

We study the radial acceleration relation (RAR) for early-type galaxies (ETGs) in the SDSS MaNGA MPL5 dataset. The complete ETG sample show a slightly offset RAR from the relation reported by McGaugh et al. (2016) at the low-acceleration end; we find that the deviation is due to the fact that the slow rotators show a systematically higher acceleration relation than the McGaugh's RAR, while the fast rotators show a consistent acceleration relation to McGaugh's RAR. There is a 1\sigma significant difference between the acceleration relations of the fast and slow rotators, suggesting that the acceleration relation correlates with the galactic spins, and that the slow rotators may have a different mass distribution compared with fast rotators and late-type galaxies. We suspect that the acceleration relation deviation of slow rotators may be attributed to more galaxy merger events, which would disrupt the original spins and correlated distributions of baryons and dark matter orbits in galaxies.

We take advantage of the statistical power of the large-volume dark-matter-only Millennium simulation, combined with a sophisticated semi-analytic galaxy formation model, to explore whether the recently reported $z=3.7$ quiescent galaxy ZF-COSMOS-20115 (ZF; Glazebrook et al. 2017) can be accommodated in current galaxy formation models. In our model, a population of quiescent galaxies (QGs) with stellar masses and star formation rates comparable to those of ZF naturally emerges at redshifts $z<4$. There are two and five ZF analogues at the redshift $3.86$ and $3.58$ in the Millennium simulation volume, respectively. We demonstrate that, while the $z>3.5$ massive QGs are rare (about 2\% of the galaxies with the similar stellar masses), the existing AGN feedback model implemented in the semi-analytic galaxy formation model can successfully explain the formation of the high-redshift QGs as it does on their lower redshift counterparts.

We present a fully integrated InGaAs/InP negative feedback avalanche diode (NFAD) based free-running single-photon detector (SPD) designed for accurate lidar applications. A free-piston Stirling cooler is used to cool down the NFAD with a large temperature range, and an active hold-off circuit implemented in a field programmable gate array is applied to further suppress the afterpulsing contribution. The key parameters of the free-running SPD including photon detection efficiency (PDE), dark count rate (DCR), afterpulse probability, and maximum count rate (MCR) are dedicatedly optimized for lidar application in practice. We then perform a field experiment using a Mie lidar system with 20 kHz pulse repetition frequency to compare the performance between the free-running InGaAs/InP SPD and a commercial superconducting nanowire single-photon detector (SNSPD). Our detector exhibits good performance with 1.6 Mcps MCR (0.6 {\mu}s hold-off time), 10% PDE, 950 cps DCR, and 18% afterpulse probability over 50 {\mu}s period. Such performance is worse than the SNSPD with 60% PDE and 300 cps DCR. However, after performing a specific algorithm that we have developed for afterpulse and count rate corrections, the lidar system performance in terms of range-corrected signal (Pr2) distribution using our SPD agrees very well with the result using the SNSPD, with only a relative error of ~2%. Due to the advantages of low-cost and small size of InGaAs/InP NFADs, such detector provides a practical solution for accurate lidar applications.

A particular population of galaxies have drawn much interest recently, which are as faint as typical dwarf galaxies but have the sizes as large as $L^*$ galaxies, the so called "ultra-diffuse galaxie" (UDGs). The lack of tidal features of UDGs in dense environments suggests that their host halos are perhaps as massive as that of the Milky Way. On the other hand, galaxy formation efficiency should be much higher in the halos of such masses. Here we use the model galaxy catalog generated by populating two large simulations: the Millennium-II cosmological simulation and Phoenix simulations of 9 big clusters with the semi-analytic galaxy formation model. This model reproduces remarkably well the observed properties of UDGs in the nearby clusters, including the abundance, profile, color, and morphology, etc. We search for UDG candidates using the public data and find 2 UDG candidates in our Local Group and 23 in our Local Volume, in excellent agreement with the model predictions. We demonstrate that UDGs are genuine dwarf galaxies, formed in the halos of $\sim 10^{10}M_{\odot}$. It is the combination of the late formation time and high-spins of the host halos that results in the spatially extended feature of this particular population. The lack of tidal disruption features of UDGs in clusters can also be explained by their late infall-time.

Massive spectroscopic redshift surveys open a promising window to accurately measure peculiar velocity at cosmological distances through redshift space distortion (RSD). In paper I of this series of work we proposed to decompose peculiar velocity into three eigen-modes (v_\delta, v_S and v_B) in order to facilitate the RSD modeling and peculiar velocity reconstruction. In the current paper we measure the dark matter RSD related statistics of the velocity eigen-modes through a set of N-body simulations, including the velocity power spectra, correlation functions, one-point probability distribution functions, cumulants and the damping functions describing the Finger of God effect. (1) The power spectrum measurement shows that these velocity components have distinctly different spatial distribution and redshift evolution. In particular, we measure the window function \tilde{W}(k,z), which describes the impact of nonlinear evolution on the v_\delta-density relation. We confirm that it can induce a significant systematic error of O(10%) in RSD cosmology. We demonstrate that \tilde{W} can be accurately described by a simple fitting formula with one or two free parameters. (2) The correlation function measurement shows that the correlation length is O(100), O(10) and O(1) Mpc for v_\delta, v_S and v_B respectively. These correlation lengths determine where we can treat the velocity fields as spatially uncorrelated. (3) The velocity PDFs and cumulants quantify non-Gaussianities of the velocity fields. We confirm speculation in paper I that v_\delta is largely Gaussian, nevertheless with non-negligible non-Gaussianity, v_B is significantly non-Gaussian. We also measure the damping functions. Despite the observed non-Gaussianities, the damping functions and hence the FOG effect are all well approximated as Gaussian ones at scales of interest.

An all-fiber, micro-pulse and eye-safe high spectral resolution wind lidar (HSRWL) at 1550nm is proposed and demonstrated by using a pair of upconversion single-photon detectors and a fiber Fabry-Perot scanning interferometer (FFP-SI). In order to improve the optical detection efficiency, both the transmission spectrum and the reflection spectrum of the FFP-SI are used for spectral analyses of the aerosol backscatter and the reference laser pulse. The reference signal is tapped from the outgoing laser and served as a zero velocity indicator. The Doppler shift is retrieved from a frequency response function Q, which is defined as the ratio of difference of the transmitted signal and the reflected signal to their sum. Taking advantages of high signal-to-noise ratio of the detectors and high spectral resolution of the FFP-SI, the Q spectra of the aerosol backscatter are reconstructed along the line-of-sight (LOS) of the telescope. By applying a least squares fit procedure to the measured Q spectra, the center frequencies and the bandwidths are obtained simultaneously. And then the Doppler shifts are determined relative to the center frequency of the reference signal. To eliminate the influence of temperature fluctuations on the FFP-SI, the FFP-SI is cased in a chamber with temperature stability of 0.001 during the measurement. Continuous LOS wind observations are carried out on two days at Hefei (31.843 N, 117.265 E), China. In the meantime, LOS wind measurements from the HSRWL show good agreement with the results from an ultrasonic wind sensor (Vaisala windcap WMT52). Due to the computational expensive of the convolution operation of the Q function, an empirical method is adopted to evaluate the quality of the measurements. The standard deviation of the wind speed is 0.76 m/s at the 1.8 km. The standard deviation of the retrieved bandwidth variation is 2.07 MHz at the 1.8 km.

For the first time, a versatile, eyesafe, compact and direct detection Doppler lidar is developed using upconversion single-photon detection method. An all-fiber and polarization maintaining architecture is realized to guarantee the high optical coupling efficiency and the system stability. Using integrated-optic components, the conservation of etendue of the optical receiver is achieved by manufacturing a fiber-coupled periodically poled Lithium niobate waveguide and an all-fiber Fabry-Perot interferometer (FPI). The so-called double-edge direct detection is implemented using a single-channel FPI and a single upconversion detector, incorporating time-division multiplexing method. The relative error of the system is lower than 0.1% over 9 weeks. To show the robust of the system, atmospheric wind and visibility over 48 hours are detected in the boundary layer. In the intercomparison experiments, lidar shows good agreement with the ultrasonic wind sensor (Vaisala windcap WMT52), with standard deviation of 1.04 m/s in speed and 12.3{\deg} in direction.

A maximum-likelihood method, tested as an unbiased estimator from numerical simulations, is used to estimate cosmic bulk flow from peculiar velocity surveys. The likelihood function is applied to four observational catalogues (ENEAR, SFI++, A1SN and SC) constructed from galaxy peculiar velocity surveys and Type-Ia supernovae data at low redshift ($z \leq 0.03$). We find that the Spiral Field {\it I}-band catalogue constrains the bulk flow to be $V=290 \pm 30 {\,{\rm km}\,{\rm s}^{-1}}$ towards $l=281^{\circ} \pm 7^{\circ}$, $b=8^{\circ +6^{\circ}}_{-5^{\circ}}$ on effective scales of $ 58 {\,{h^{-1} {\rm Mpc}}}$, which is the tightest constraints achievable at the present time. By comparing the amplitudes of our estimated bulk flows with theoretical prediction, we find excellent agreement between the two. In addition, directions of estimated bulk flows are also consistent with measurements in other studies.

We present a spectroscopic catalog of 67082 M dwarfs from the LAMOST pilot survey. For each spectrum of the catalog, spectral subtype, radial velocity, equivalent width of H${\alpha}$, a number of prominent molecular band indices and the metal sensitive parameter $\zeta$ are provided . Spectral subtype have been estimated by a remedied Hammer program (Original Hammer: Covey et al. 2007), in which indices are reselected to obtain more accurate auto-classified spectral subtypes. All spectra in this catalog have been visually inspected to confirm the spectral subtypes. Radial velocities have been well measured by our developed program which uses cross-correlation method and estimates uncertainty of radial velocity as well. We also examine the magnetic activity properties of M dwarfs traced by H${\alpha}$ emission line. The molecular band indices included in this catalog are temperature or metallicity sensitive and can be used for future analysis of the physical properties of M dwarfs. The catalog is available on the website \url{http://sciwiki.lamost.org/MCatalogPilot/}.

Massive spectroscopic surveys will measure the redshift space distortion (RSD) induced by galaxy peculiar velocity to unprecedented accuracy. We develop a new method to improve the RSD modeling and to carry out robust reconstruction of the 3D large scale peculiar velocity through galaxy redshift surveys, in light of RSD. (1) We propose a mathematically unique and physically motivated decomposition of peculiar velocity into three eigen-components of different origins, different scale dependences and different impacts on RSD. (2) This decomposition has the potential to simplify and improve the RSD modeling. We derive a new formula for the redshift space power spectrum. Under the velocity decomposition scheme, all high order Gaussian corrections and non-Gaussian correction of order $\delta^3$ can be taken into account without introducing extra model uncertainties. We find that the leading order contribution $\propto (1+f\tilde{W}(k)u^2)^2$ differs from the Kaiser formula by a function $\tilde{W}(k)\leq 1$. We hence identify a significant systematic error causing underestimation of the structure growth parameter $f$ by as much as $O(10%)$ even at relatively large scale $k=0.1h/$Mpc. The velocity decomposition also reveals the three origins of the finger of God (FOG) effect. (3) The velocity decomposition clarifies issues in peculiar velocity reconstruction through 3D galaxy distribution. We discuss two possible ways to carry out the velocity reconstruction. Both use the otherwise troublesome RSD in velocity reconstruction as a valuable source of information. Both have the advantage to render the reconstruction of a stochastic 3D field into the reconstruction of a deterministic window function $W^s(k,u)$ of limited degrees of freedom...[abridged]

Analysis of the Pangu N-body simulation validates that the bulk flow of halos follows a Maxwellian distribution which variance is consistent with the prediction of the linear theory of structure formation. We propose that the consistency between the observed bulk velocity and theories should be examined at the effective scale of the radius of a spherical top-hat window function yielding the same smoothed velocity variance in linear theory as the sample window function does. We compared some recently estimated bulk flows from observational samples with the prediction of the \Lambda CDM model we used; some results deviate from expectation at a level of ~ 3\sigma but the discrepancy is not as severe as previously claimed. We show that bulk flow is only weakly correlated with the dipole of the internal mass distribution, the alignment angle between the mass dipole and the bulk flow has a broad distribution peaked at ~ 30-50 deg., and also that the bulk flow shows little dependence on the mass of the halos used in the estimation. In a simulation of box size 1Gpc/h, for a cell of radius 100 Mpc/h the maximal bulk velocity is >500 km/s, dipoles of the environmental mass outside the cell are not tightly aligned with the bulk flow, but are rather located randomly around it with separation angles ~ 20-40 deg. In the fastest cell there is a slightly smaller number of low-mass halos; however halos inside are clustered more strongly at scales > ~ 20 Mpc/h, which might be a significant feature since the correlation between bulk flow and halo clustering actually increases in significance beyond such scales.

Baryonic physical processes could leave non-negligible imprint on cosmic matter distribution pattern. Series of high precision simulation data sets with identical initial condition are employed for count-in-cell (CIC) analysis, including one N-body dark matter run, one with adiabatic gas only and one with dissipative processes. Variances and higher order correlation functions of dark matter and gas are estimated. It is found that baryon physical processes mainly affected dark matter distribution at scales less than $1h^{-1}$Mpc. In comparison with the pure dark matter run, adiabatic process alone strengthens variance of dark matter by \sim 10% at scale $0.1h^{-1}$Mpc, while $S_n$s of dark matter deviate from pure dark matter case only mildly at a few percentages. Dissipative gas run does not differ much to the adiabatic run in dark matter variance, but renders significantly different $S_n$ parameters of dark matter, bringing about more than 10% enhancement to $S_3$ at $0.1h^{-1}$Mpc and $z=0$. Distribution patterns of gas in two hydrodynamical simulations are prominently different. Variance of gas at $z=0$ decreases by $\sim 30%$ in adiabatic simulation while by $\sim 60%$ in non-adiabatic simulation at $0.1h^{-1}$Mpc, the attenuation is weaker at larger scales but still obvious at $\sim 10h^{-1}$Mpc. $S_n$ parameters of gas are biased upward at scales $< \sim 4h^{-1}$Mpc, dissipative processes give $\sim 84%$ promotion at $z=0$ to $S_3$ at $0.1h^{-1}$Mpc against the moderate $\sim 7%$ in adiabatic simulation. The clustering segregation we observed between gas and dark matter could have intricate implication on modeling galaxy distribution and relevant cosmological application demanding fine details of matter distribution in strongly nonlinear regime.

Flux-limited and volume-limited galaxy samples are constructed from SDSS data releases DR4, DR6 and DR7 for statistical analysis. The two-point correlation functions $\xi(s)$, monopole of three-point correlation functions $\zeta_0$, projected two-point correlation function $w_p$ and pairwise velocity dispersion $\sigma_{12}$ are measured to test if galaxy samples are fair for these statistics. We find that with increment of sky coverage of SDSS, $\xi(s)$ of flux-limited sample is extremely robust and insensitive to local structures at low redshift. But for volume-limited samples fainter than $L^*$ at large scales $s>\sim 10\hmpc$, deviation of $\xi(s)$ and $\zeta_0$ of DR7 to those of DR4 and DR6 increases with larger absolute magnitude. In the weakly nonlinear regime, there is no agreement between $\zeta_0$ of different data releases in all luminosity bins. Furthermore, $w_p$ of volume-limited samples of DR7 in luminosity bins fainter than $-M_{r,0.1}=[18.5,19.5]$ are significantly larger, and $\sigma_{12}$ of the two faintest volume-limited samples of DR7 display very different scale dependence than results of DR4 and DR6. Our findings call for cautions in understanding clustering analysis results of SDSS faint galaxy samples, and higher order statistics of SDSS volume-limited samples in the weakly nonlinear regime. The first zero-crossing points of $\xi(s)$ of volume-limited samples are also investigated and discussed.

The kinetic Sunyaev Zel'dovich effect (kSZ) effect is a potentially powerful probe to the missing baryons. However, the kSZ signal is overwhelmed by various contaminations and the cosmological application is hampered by loss of redshift information due to the projection effect. We propose a kSZ tomography method to alleviate these problems, with the aid of galaxy spectroscopic redshift surveys. We propose to estimate the large scale peculiar velocity through the 3D galaxy distribution, weigh it by the 3D galaxy density and adopt the product projected along the line of sight with a proper weighting as an estimator of the true kSZ temperature fluctuation $\Theta$. We thus propose to measure the kSZ signal through the $\Hat{\Theta}$-$\Theta$ cross correlation. This approach has a number of advantages (see details in the abstract of the paper). We test the proposed kSZ tomography against non-adiabatic and adiabatic hydrodynamical simulations. We confirm that $\hat{\Theta}$ is indeed tightly correlated with $\Theta$ at $k\la 1h/$Mpc, although nonlinearities in the density and velocity fields and nonlinear redshift distortion do weaken the tightness of the $\hat{\Theta}$-$\Theta$ correlation. We further quantify the reconstruction noise in $\Hat{\Theta}$ from galaxy distribution shot noise. Based on these results, we quantify the applicability of the proposed kSZ tomography for future surveys. We find that, in combination with the BigBOSS-N spectroscopic redshift survey, the PLANCK CMB experiment will be able to detect the kSZ with an overall significance of $\sim 50\sigma$ and further measure its redshift distribution at many redshift bins over $0<z<2$.

The zeroth-order component of the cosine expansion of the projected three-point correlation function is proposed for clustering analysis of cosmic large scale structure. These functions are third order statistics but can be measured similarly to the projected two-point correlations. Numerical experiments with N-body simulations indicate that the advocated statistics are redshift distortion free within 10% in the non-linear regime on scales ~0.2-10Mpc/h. Halo model prediction of the zeroth-order component of the projected three-point correlation function agrees with simulations within ~10%. This lays the ground work for using these functions to perform joint analyses with the projected two-point correlation functions, exploring galaxy clustering properties in the framework of the halo model and relevant extensions.

The spatial distribution of galaxies we observed is subject to the given condition that we, human beings are sitting right in a galaxy -- the Milky Way. Thus the ergodicity assumption is questionable in interpretation of the observed galaxy distribution. The resultant difference between observed statistics (volume average) and the true cosmic value (ensemble average) is termed as the ergodicity bias. We perform explicit numerical investigation of the effect for a set of galaxy survey depths and near-end distance cuts. It is found that the ergodicity bias in observed two- and three-point correlation functions in most cases is insignificant for modern analysis of samples from galaxy surveys and thus close a loophole in precision cosmology. However, it may become non-negligible in certain circumstances, such as those applications involving three-point correlation function at large scales of local galaxy samples. Thus one is reminded to take extra care in galaxy sample construction and interpretation of the statistics of the sample, especially when the characteristic redshift is low.

We present an analysis of timing residual (noise) of 54 pulsars obtained from 25-m radio telescope at Urumqi Observatory with a time span of 5~8 years, dealing with statistics of the Hurst parameter. The majority of these pulsars were selected to have timing noise that look like white noise rather than smooth curves. The results are compared with artificial series of different constant pairwise covariances. Despite the noise like appearance, many timing residual series showed Hurst parameters significantly deviated from that of independent series. We concluded that Hurst parameter may be capable of detecting dependence in timing residual and of distinguishing chaotic behavior from random processes.

We present detailed clustering analysis of a large K-band selected local galaxy sample, which is constructed from the 2MASS and the SDSS and consists of $82,486$ galaxies with $10 < K < 13.5$ and $0.01 < z < 0.1$. The two-point correlation function of the magnitude-limited sample in real space at small scales is well described by a power law $\xi(r)=(r/6.44\pm0.23)^{-1.81\pm0.02}$. The pairwise velocity dispersion is derived from the anisotropic two-point correlation function and find the dispersion $\sigma_{12}=685\pm 17\kms$ if its scale invariance is assumed, which is larger than values measured in optical bands selected galaxy samples. We further investigate the dependence of the two-point correlation function and the $\sigma_{12}$ on the $g-r$ color and the $K$-band luminosity, obtain similar results to previous works in optical bands. Comparing a mock galaxy sample with our real data indicates that the semi-analytical model can not mimic the $\sigma_{12}$ in observation albeit it can approximate the two-point correlation function within measurement uncertainties.

We present the current status of the QUIJOTE (Q-U-I JOint TEnerife) CMB Experiment, a new instrument which will start operations early 2009 at Teide Observatory, with the aim of characterizing the polarization of the CMB and other processes of galactic and extragalactic emission in the frequency range 10-30 GHz and at large angular scales. QUIJOTE will be a valuable complement at low frequencies for the PLANCK mission, and will have the required sensitivity to detect a primordial gravitational-wave component if the tensor-to-scalar ratio is larger than r=0.05.

The excursion set theory of halo formation is modified by adopting the fractional Brownian motion, to account for possible correlation between merging steps. We worked out analytically the conditional mass function, halo merging rate and formation time distribution in the spherical collapse model. We also developed an approximation for the ellipsoidal collapse model and applied it to the calculation of the conditional mass function and the halo formation time distribution. For models in which the steps are positively correlated, the halo merger rate is enhanced when the accreted mass is less than $\sim 25M^*$, while for the negatively correlated case this rate is reduced. Compared with the standard model in which the steps are uncorrelated, the models with positively correlated steps produce more aged population in small mass halos and more younger population in large mass halos, while for the models with negatively correlated steps the opposite is true. An examination of simulation results shows that a weakly positive correlation between successive merging steps appears to fit best. We have also found a systematic effect in the measured mass function due to the finite volume of simulations. In future work, this will be included in the halo model to accurately predict the three point correlation function estimated from simulations.

Scale transformations have played an extremely successful role in studies of cosmological large-scale structure by relating the non-linear spectrum of cosmological density fluctuations to the linear primordial power at longer wavelengths. Here we generalize this approach to investigate the usefulness of scale transformations for nonlinear higher-order statistics, specifically the bispectrum. We find that the bispectrum predicted by perturbation theory at tree-level can be rescaled to match the results of full numerical simulations in the weakly and intermediately nonlinear regimes, especially at high redshifts, with an accuracy that is surprising given the simplicity of the procedure used. This discovery not only offers a simple practical way of calculating the matter bispectrum, but also suggests that scale transformations may yet yield even deeper insights into the physics of hierarchical clustering.

The fractional Brownian motion with index $\alpha$ is introduced to construct the fractional excursion set model. A new mass function with single parameter $\alpha$ is derived within the formalism, of which the Press-Schechter mass function (PS) is a special case when $\alpha=1/2$. Although the new mass function is computed assuming spherical collapse, comparison with the Sheth-Tormen fitting function (ST) shows that the new mass function of $\alpha\approx 0.435$ agrees with ST remarkably well in high mass regime, while predicts more small mass halos than the ST but less than the PS. The index $\alpha$ is the Hurst exponent, which exact value in context of structure formation is modulated by properties of the smoothing window function and the shape of power spectrum. It is conjectured that halo merging rate and merging history in the fractional set theory might be imprinted with the interplay between halos at small scales and their large scale environment. And the mass function in high mass regime can be a good tool to detect the non-Gaussianity of the initial density fluctuation.

The power spectrum of the two-degree Field Galaxy Redshift Survey (2dFGRS) sample is estimated with the discrete wavelet transform (DWT) method. The DWT power spectra within $0.04 <k< 2.3 h$Mpc$^{-1}$ are measured for three volume-limited samples defined in connective absolute magnitude bins $-19 \sim -18$, $-20 \sim -19$ and $-21 \sim -20$. We show that the DWT power spectrum can effectively distinguish $\Lambda$CDM models of $\sigma_8=0.84$ and $\sigma_8=0.74$. We adopt maximum likelihood method to perform three-parameter fitting with bias parameter $b$, pairwise velocity dispersion $\sigma_{pv}$ and redshift distortion parameter $\beta=\Omega_m^{0.6}/b$ to the measured DWT power spectrum. Fitting results denotes that in a $\sigma_8=0.84$ universe the best fitted $\Omega_m$ given by the three samples are consistent in the range $0.28 \sim 0.36$, and the best fitted $\sigma_{pv}$ are $398^{+35}_{-27}$, $475^{+37}_{-29}$ and $550 \pm 20$km/s for the three samples, respectively. However in the model of $\sigma_8=0.74$, our three samples give very different values of $\Omega_m$. We repeat the fitting by using empirical formula of redshift distortion. The result of the model of low $\sigma_8$ is still poor, especially, one of the best value $\sigma_{pv}$ is as large as $10^3$km/s. The power spectrum of 2dFGRS seems in disfavor of models with low amplitude of density fluctuations.

Multifractal measures are applied to three IRAS galaxy subsamples selected by colour from the PSCz catalogue. As shown by generalised dimension spectrum, hot IRAS galaxies are found less clustered than cold galaxies, but the difference is very weak. An alternative tool, the conditional multifractal dimension spectrum reveals apparent peculiarity of the distribution of hot galaxies, especially at high orders. The conditional multifractal measure is basically measuring the environment of selected galaxies. A detailed analysis of the distribution of galaxies with their number of neighbours turns out that cold galaxies are more likely living in over dense regions than hot galaxies. Further studies suggest that, since the colour of IRAS galaxies is a good trace of star formation rate, it is possible that we have statistical evidence here for enhanced star formation rate due to galaxy interactions

We explore features of redshift distortion in Fourier analysis of N-body simulations. The phases of the Fourier modes of the dark matter density fluctuation are generally shifted by the peculiar motion along the line of sight, the induced phase shift is stochastic and has probability distribution function (PDF) symmetric to the peak at zero shift while the exact shape depends on the wave vector, except on very large scales where phases are invariant by linear perturbation theory. Analysis of the phase shifts motivates our phenomenological models for the bispectrum in redshift space. Comparison with simulations shows that our toy models are very successful in modeling bispectrum of equilateral and isosceles triangles at large scales. In the second part we compare the monopole of the power spectrum and bispectrum in the radial and plane-parallel distortion to test the plane-parallel approximation. We confirm the results of Scoccimarro (2000) that difference of power spectrum is at the level of 10%, in the reduced bispectrum such difference is as small as a few percents. However, on the plane perpendicular to the line of sight of k_z=0, the difference in power spectrum between the radial and plane-parallel approximation can be more than 10%, and even worse on very small scales. Such difference is prominent for bispectrum, especially for those configurations of tilted triangles. The non-Gaussian signals under radial distortion on small scales are systematically biased downside than that in plane-parallel approximation, while amplitudes of differences depend on the opening angle of the sample to the observer. The observation gives warning to the practice of using the power spectrum and bispectrum measured on the k_z=0 plane as estimation of the real space statistics.

[Abridged] We use mock catalogues based on the GALICS model (Hatton et al. 03) to explore the nature of galaxy clustering observed in the SDSS. We measure low and high order angular clustering statistic from these mock catalogues, after selecting galaxies the same way as for observations, and compare them directly to estimates from SDSS data. Note that we also present measurements of S3-S5 on the SDSS DR1. We find that our model is in general good agreement with observations in the scale/luminosity range where we can trust the predictions. This range is found to be limited (i) by the size of the dark matter simulation used -- which introduces finite volume effects at large scales -- and by the mass resolution of this simulation -- which introduces incompleteness at apparent magnitudes fainter than $r\sim 20$. We then focus on the small scale clustering properties of galaxies and investigate the behaviour of three different prescriptions for positioning galaxies within haloes of dark matter. We show that galaxies are poor tracers both of DM particles or DM sub-structures, within groups and clusters. Instead, SDSS data tells us that the distribution of galaxies lies somewhat in between these two populations. This confirms the general theoretical expectation from numerical simulations and semi-analytic modelling.

We perform detailed comparison of the semi-analytic halo model predictions with measurements in numerical simulations of the two and three point correlation functions (3PCF), as well as power spectrum and bispectrum. We discuss the accuracy and self-consistency of the halo model description of gravitational clustering in the non-linear regime and constrain halo model parameters. We exploit the recently proposed multipole expansion of three point statistics that expresses rotation invariance in the most natural way. This not only offers technical advantages by reducing the integrals required for the halo model predictions, but amounts to a convenient way of compressing the information contained in the 3PCF. We find that, with an appropriate choice of the halo boundary and mass function cut-off, halo model predictions are in good agreement with the bispectrum measured in numerical simulations. However, the halo model predicts less than the observed configuration dependence of the 3PCF on ~ Mpc scales. This effect is mainly due to quadrupole moment deficit, possibly related to the assumption of spherical halo geometry. Our analysis shows that using its harmonic decomposition, the full configuration dependence of the 3PCF in the non-linear regime can be compressed into just a few numbers, the lowest multipoles. Moreover, these multipoles are closely related to the highest signal to noise eigenmodes of the 3PCF. Therefore this estimator may simplify future analyses aimed at constraining cosmological and halo model parameters from observational data.

We measure the monopole moment of the three-point correlation function on scales $1\mpc-70\mpc$ in the Two degree Field Galaxy Redshift Survey (2dFGRS). Volume limited samples are constructed using a series of integral magnitudes bins between $M = -18 ... -22$. Our measurements with a novel edge-corrected estimator represent most, if not all, three-point level monopole or angular averaged information in the catalogue. We fit a perturbative non-linear bias model to a joint data vector formed from the estimated two- and three-point correlation functions. Two different models are used: an analytic model based on Eulerian perturbation theory including bias and redshift distortions, and a phenomenological bias model based on the direct redshift space measurements in the large Virgo simulations. To interpret the clustering results, we perform a three parameter Gaussian maximum likelihood analysis. In the canonical $-21 \sim -20$ volume limited sample we find $\sigma_8 = 0.93^{+0.06}_{-0.2}$, $b = 1.04^{+0.23}_{-0.09}$, and $b_2 = -0.06^{+0.003}_{-0.001}$. Our estimate of $\sigma_8$, is robust across the different volume limited samples constructed. These results, based solely on the large scale clustering of galaxies, are in excellent agreement with previous analyses using the Wilkinson Anisotropy Probe: this is a spectacular success of the concordance model. We also present two-parameter fits for the bias parameters, which are in excellent agreement with previous findings of the bias evolution in the 2dFGRS.

We present two related techniques to measure the two-point correlation function and the power spectrum with edge correction in any spatial dimensions. The underlying algorithm uses fast Fourier transforms for calculating the two-point function with an heuristically weighted edge corrected estimator. Once the correlation function is measured, we estimate the power spectrum via numerical integration of the Hankel transform connecting the two. We introduce an efficient numerical technique based on Gauss-Bessel-quadrature and double exponential transformation. This, combined with our, or any similar, two-point function estimator leads to a novel edge corrected estimator for power spectra. The pair of techniques presented are the Euclidean analogs of those developed and widely used in cosmic microwave background research for spherical maps.

We present measurements of the redshift--space three-point correlation function of galaxies in the Sloan Digital Sky Survey (SDSS). For the first time, we analyze the dependence of this statistic on galaxy morphology, color and luminosity. In order to control systematics due to selection effects, we used $r$--band, volume-limited samples of galaxies, constructed from the magnitude-limited SDSS data ($14.5<r<17.5$), and further divided the samples into two morphological types (early and late) or two color populations (red and blue). The three-point correlation function of SDSS galaxies follow the hierarchical relation well and the reduced three-point amplitudes in redshift--space are almost scale-independent ($Q_z=0.5\sim1.0$). In addition, their dependence on the morphology, color and luminosity is not statistically significant. Given the robust morphological, color and luminosity dependences of the two-point correlation function, this implies that galaxy biasing is complex on weakly non-linear to non-linear scales. We show that simple deterministic linear relation with the underlying mass could not explain our measurements on these scales.

Conditional cumulants form a set of unique statistics which represent a sensible compromise between N-point correlation functions and cumulants measured from moments of counts in cells. They share accurate edge corrected estimators with $N$-point correlation functions, yet, they are as straightforward to measure and interpret as counts in cells. The conditional cumulants have three equivalent views as i) degenerate N-point correlation functions ii) or integrated monopole moments of the bispectrum iii) they are closely related to neighbour counts. We compute the predictions of weakly non-linear perturbation theory for conditional cumulants and compare them with measurements in simulations, both in real and redshift space. We find excellent agreement between theory and simulations, especially on scales >~20Mpc. Due to their advantageous statistical properties and well understood dynamics, we propose conditional cumulants as tools for high precision cosmology. Potential applications include constraining bias and redshift distortions from galaxy redshift surveys

We present a simple and accurate method to constrain galaxy bias based on the distribution of counts in cells. The most unique feature of our technique is that it is applicable to non-linear scales, where both dark matter statistics and the nature of galaxy bias are fairly complex. First, we estimate the underlying continuous distribution function from precise counts-in-cells measurements assuming local Poisson sampling. Then a robust, non-parametric inversion of the bias function is recovered from the comparison of the cumulative distributions in simulated dark matter and galaxy catalogs. Obtaining continuous statistics from the discrete counts is the most delicate novel part of our recipe. It corresponds to a deconvolution of a (Poisson) kernel. For this we present two alternatives: a model independent algorithm based on Richardson-Lucy iteration, and a solution using a parametric skewed lognormal model. We find that the latter is an excellent approximation for the dark matter distribution, but the model independent iterative procedure is more suitable for galaxies. Tests based on high resolution dark matter simulations and corresponding mock galaxy catalogs show that we can reconstruct the non-linear bias function down to highly non-linear scales with high precision in the range of $-1 \le \delta \le 5$. As far as the stochasticity of the bias, we have found a remarkably simple and accurate formula based on Poisson noise, which provides an excellent approximation for the scatter around the mean non-linear bias function. In addition we have found that redshift distortions have a negligible effect on our bias reconstruction, therefore our recipe can be safely applied to redshift surveys.

The analysis of redshift surveys with fractal tools requires one to apply some form of statistical correction for galaxies lying near the geometric boundary of the sample. In this paper we compare three different methods of performing such a correction upon estimates of the correlation integral in order to assess the extent to which estimates may be biased by boundary terms. We apply the corrections illustrative examples, including a simple fractal set (L\'{e}vy flight), a random $\beta$-model, and a subset of the CfA2 Southern Cap survey. This study shows that the new ``angular'' correction method we present is more generally applicable than the other methods used to date: the conventional ``capacity'' correction imposes a bias towards homogeneity, and the ``deflation'' method discards large-scale information, and consequently reduces the statistical usefulness of data sets. The ``angular'' correction method is effective at recovering true fractal dimensions, although the extent to which boundary corrections are important depends on the form of fractal distribution assumed as well as the details of the survey geometry. We also show that the CfA2 Southern sample does not show any real evidence of a transition to homogeneity. We then revisit the IRAS PSCz survey and ``mock'' PSCz catalogues made using N-body simulations of two different cosmologies. The results we obtain from the PSCz survey are not significantly affected by the form of boundary correction used, confirming that the transition from fractal to homogeneous behaviour reported by Pan & Coles (2000) is real.

We investigate the behaviour of galaxy clustering on large scales using the PSCz catalogue. In particular, we ask whether there is any evidence of large-scale fractal behaviour in this catalogue. We find the correlation dimension in this survey varies with scale, consistent with other analyses. For example, our results on small and intermediate scales are consistent those obtained from the QDOT sample, but the larger PSCz sample allows us to extend the analysis out to much larger scales. We find firm evidence that the sample becomes homogeneous at large scales; the correlation dimension of the sample is D_2=2.992(3) for r>30 h^{-1} Mpc. This provides strong evidence in favour of a Universe which obeys the Cosmological Principle.