In this paper, we investigate whether galaxy assembly bias can reconcile the 20-40% disagreement between the observed galaxy projected clustering signal and the galaxy-galaxy lensing signal in the BOSS CMASS galaxy sample reported in Leauthaud et al. (2017). We use the suite of AbacusCosmos Lamda-CDM simulations at Planck best-fit cosmology and two flexible implementations of extended halo occupation distribution (HOD) models that incorporate galaxy assembly bias to build forward models and produce joint fits of the observed galaxy clustering signal and the galaxy-galaxy lensing signal. We find that our models using the standard HODs without any assembly bias generalizations continue to show a 20-40% over-prediction of the observed galaxy-galaxy lensing signal. We find that our implementations of galaxy assembly bias do not reconcile the two measurements at Planck best-fit cosmology. In fact, despite incorporating galaxy assembly bias, the satellite distribution parameter, and the satellite velocity bias parameter into our extended HOD model, our fits still strongly suggest a 31-34% discrepancy between the observed projected clustering and galaxy-galaxy lensing measurements. It remains to be seen whether a combination of other galaxy assembly bias models, alternative cosmological parameters, or baryonic effects can explain the amplitude difference between the two signals.

We present the first images of a coordinated campaign to follow active region NOAA 12709 on 2018 May 13 as part of a joint effort between three observatories (China-Europe). The active region was close to disk center and enclosed a small pore, a tight polarity inversion line and a filament in the chromosphere. The active region was observed with the 1.5-meter GREGOR solar telescope on Tenerife (Spain) with spectropolarimetry using GRIS in the He I 10830 \r{A} spectral range and with HiFI using two broad-band filter channels. In addition, the Lomnicky Stit Observatory (LSO, Slovakia) recorded the same active region with the new Solar Chromospheric Detector (SCD) in spectroscopic mode at H$\alpha$ 6562 \r{A}. The third ground-based telescope was located at the Fuxian Solar Observatory (China), where the active region was observed with the 1-meter New Vacuum Solar Telescope (NVST), using the Multi-Channel High Resolution Imaging System at H$\alpha$ 6562 \r{A}. Overlapping images of the active region from all three telescopes will be shown as well as preliminary Doppler line-of-sight (LOS) velocities. The potential of such observations are discussed.

This paper is devoted to the primary spectro-polarimetric observation performed at the New Vacuum Solar Telescope of China since 2017, and our aim is to precisely evaluate the real polarimetric accuracy and sensitivity of this polarimetry by using full Stokes spectro-polarimetric observations of the photospheric line Fe I 532.4 nm. In the work, we briefly describe the salient characteristic of the NVST as a polarimeter in technology and then characterize its instrumental polarization based on the operation in 2017 and 2019. It is verified that the calibration method making use of the instrumental polarization calibration unit (ICU) is stable and credible. The calibration accuracy can reach up to 3$\times 10^{-3}$ . Based on the scientific observation of the NOAA 12645 on April 5th, 2017, we estimate that the residual cross-talk from Stokes $I$ to Stokes $Q$, $U$ and $V$, after the instrumental polarization calibration, is about 4$\times10^{-3}$ on average, which is consistent with the calibration accuracy and close to the photon noise. The polarimetric sensitivity (i.e., the detection limit) for polarized light is of the order of $10^{-3}$ with an integration time over 20 seconds. Slow modulation rate is indeed an issue for the present system. The present NVST polarimeter is expected to be integrated with an high-order adaptive optics system and a field scanner to realize 2D magnetic field vector measurements in the following instrumentation update.

We present a deep machine learning (ML)-based technique for accurately determining $\sigma_8$ and $\Omega_m$ from mock 3D galaxy surveys. The mock surveys are built from the AbacusCosmos suite of $N$-body simulations, which comprises 40 cosmological volume simulations spanning a range of cosmological models, and we account for uncertainties in galaxy formation scenarios through the use of generalized halo occupation distributions (HODs). We explore a trio of ML models: a 3D convolutional neural network (CNN), a power-spectrum-based fully connected network, and a hybrid approach that merges the two to combine physically motivated summary statistics with flexible CNNs. We describe best practices for training a deep model on a suite of matched-phase simulations and we test our model on a completely independent sample that uses previously unseen initial conditions, cosmological parameters, and HOD parameters. Despite the fact that the mock observations are quite small ($\sim0.07h^{-3}\,\mathrm{Gpc}^3$) and the training data span a large parameter space (6 cosmological and 6 HOD parameters), the CNN and hybrid CNN can constrain $\sigma_8$ and $\Omega_m$ to $\sim3\%$ and $\sim4\%$, respectively.

Weak lensing peak abundance analyses have been applied in different surveys and demonstrated to be a powerful statistics in extracting cosmological information complementary to cosmic shear two-point correlation studies. Future large surveys with high number densities of galaxies enable tomographic peak analyses. Focusing on high peaks, we investigate quantitatively how the tomographic redshift binning can enhance the cosmological gains. We also perform detailed studies about the degradation of cosmological information due to photometric redshift (photo-z) errors. We show that for surveys with the number density of galaxies $\sim40\,{\rm arcmin^{-2}}$, the median redshift $\sim1$, and the survey area of $\sim15000\,{\rm deg^{2}}$, the 4-bin tomographic peak analyses can reduce the error contours of $(\Omega_{{\rm m}},\sigma_{8})$ by a factor of $5$ comparing to 2-D peak analyses in the ideal case of photo-z error being absent. More redshift bins can hardly lead to significantly better constraints. The photo-z error model here is parametrized by $z_{{\rm bias}}$ and $\sigma_{{\rm ph}}$ and the fiducial values of $z_{{\rm bias}}=0.003$ and $\sigma_{{\rm ph}}=0.02$ is taken. We find that using tomographic peak analyses can constrain the photo-z errors simultaneously with cosmological parameters. For 4-bin analyses, we can obtain $\sigma(z_{{\rm bias}})/z_{{\rm bias}}\sim10\%$ and $\sigma(\sigma_{{\rm ph}})/\sigma_{{\rm ph}}\sim5\%$ without assuming priors on them. Accordingly, the cosmological constraints on $\Omega_{{\rm m}}$ and $\sigma_{8}$ degrade by a factor of $\sim2.2$ and $\sim1.8$, respectively, with respect to zero uncertainties on photo-z parameters. We find that the uncertainty of $z_{{\rm bias}}$ plays more significant roles in degrading the cosmological constraints than that of $\sigma_{{\rm ph}}$.

Covariance matrix estimation is a persistent challenge for cosmology, often requiring a large number of synthetic mock catalogues. The off-diagonal components of the covariance matrix also make it difficult to show representative error bars on the 2-point correlation function (2PCF), since errors computed from the diagonal values of the covariance matrix greatly underestimate the uncertainties. We develop a routine for decorrelating the projected and anisotropic 2PCF with simple and scale-compact transformations on the 2PCF. These transformation matrices are modeled after the Cholesky decomposition and the symmetric square root of the Fisher matrix. Using mock catalogues, we show that the transformed projected and anisotropic 2PCF recover the same structure as the original 2PCF, while producing largely decorrelated error bars. Specifically, we propose simple Cholesky based transformation matrices that suppress the off-diagonal covariances on the projected 2PCF by ~95% and that on the anisotropic 2PCF by ~87%. These transformations also serve as highly regularized models of the Fisher matrix, compressing the degrees of freedom so that one can fit for the Fisher matrix with a much smaller number of mocks.

We present our deep $Y$-band imaging data of a two square degree field within the F22 region of the VIMOS VLT Deep Survey. The observations were conducted using the WIRCam instrument mounted at the Canada--France--Hawaii Telescope (CFHT). The total on-sky time was 9 hours, distributed uniformly over 18 tiles. The scientific goals of the project are to select faint quasar candidates at redshift $z>2.2$, and constrain the photometric redshifts for quasars and galaxies. In this paper, we present the observation and the image reduction, as well as the photometric redshifts that we derived by combining our $Y$-band data with the CFHTLenS $u^*g'r'i'z'$ optical data and UKIDSS DXS $JHK$ near-infrared data. With $J$-band image as reference total $\sim$80,000 galaxies are detected in the final mosaic down to $Y$-band $5\sigma$ point source limiting depth of 22.86 mag. Compared with the $\sim$3500 spectroscopic redshifts, our photometric redshifts for galaxies with $z<1.5$ and $i'\lesssim24.0$ mag have a small systematic offset of $|\Delta{z}|\lesssim0.2$, 1$\sigma$ scatter $0.03<\sigma_{\Delta z} < 0.06$, and less than 4.0% of catastrophic failures. We also compare to the CFHTLenS photometric redshifts, and find that ours are more reliable at $z\gtrsim0.6$ because of the inclusion of the near-infrared bands. In particular, including the $Y$-band data can improve the accuracy at $z\sim 1.0-2.0$ because the location of the 4000\AA-break is better constrained. The $Y$-band images, the multi-band photometry catalog and the photometric redshifts are released at \url{http://astro.pku.edu.cn/astro/data/DYI.html}.

We examine the possible influence of Earth-directed coronal mass ejections (ECMEs) on the Sun's shadow in the 3~TeV cosmic-ray intensity observed by the Tibet-III air shower (AS) array. We confirm a clear solar-cycle variation of the intensity deficit in the Sun's shadow during ten years between 2000 and 2009. This solar-cycle variation is overall reproduced by our Monte Carlo (MC) simulations of the Sun's shadow based on the potential field model of the solar magnetic field averaged over each solar rotation period. We find, however, that the magnitude of the observed intensity deficit in the Sun's shadow is significantly less than that predicted by MC simulations, particularly during the period around solar maximum when a significant number of ECMEs is recorded. The $\chi^2$ tests of the agreement between the observations and the MC simulations show that the difference is larger during the periods when the ECMEs occur, and the difference is reduced if the periods of ECMEs are excluded from the analysis. This suggests the first experimental evidence of the ECMEs affecting the Sun's shadow observed in the 3~TeV cosmic-ray intensity.

The convergence peaks, constructed from galaxy shape measurement in weak lensing, is a powerful probe of cosmology as the peaks can be connected with the underlined dark matter haloes. However the capability of convergence peak statistic is affected by the noise in galaxy shape measurement, signal to noise ratio as well as the contribution from the projected mass distribution from the large-scale structures along the line of sight (LOS). In this paper we use the ray-tracing simulation on a curved sky to investigate the correspondence between the convergence peak and the dark matter haloes at the LOS. We find that, in case of no noise and for source galaxies at $z_{\rm s}=1$, more than $65\%$ peaks with $\text{SNR} \geq 3$ (signal to noise ratio) are related to more than one massive haloes with mass larger than $10^{13} {\rm M}_{\odot}$. Those massive haloes contribute $87.2\%$ to high peaks ($\text{SNR} \geq 5$) with the remaining contributions are from the large-scale structures. On the other hand, the peaks distribution is skewed by the noise in galaxy shape measurement, especially for lower SNR peaks. In the noisy field where the shape noise is modelled as a Gaussian distribution, about $60\%$ high peaks ($\text{SNR} \geq 5$) are true peaks and the fraction decreases to $20\%$ for lower peaks ($ 3 \leq \text{SNR} < 5$). Furthermore, we find that high peaks ($\text{SNR} \geq 5$) are dominated by very massive haloes larger than $10^{14} {\rm M}_{\odot}$.

We present the GeneRalized ANd Differentiable Halo Occupation Distribution (GRAND-HOD) routine that generalizes the standard 5 parameter halo occupation distribution model (HOD) with various halo-scale physics and assembly bias. We describe the methodology of 4 different generalizations: satellite distribution generalization, velocity bias, closest approach distance generalization, and assembly bias. We showcase the signatures of these generalizations in the 2-point correlation function (2PCF) and the squeezed 3-point correlation function (squeezed 3PCF). We identify generalized HOD prescriptions that are nearly degenerate in the projected 2PCF and demonstrate that these degeneracies are broken in the redshift-space anisotropic 2PCF and the squeezed 3PCF. We also discuss the possibility of identifying degeneracies in the anisotropic 2PCF and further demonstrate the extra constraining power of the squeezed 3PCF on galaxy-halo connection models. We find that within our current HOD framework, the anisotropic 2PCF can predict the squeezed 3PCF better than its statistical error. This implies that a discordant squeezed 3PCF measurement could falsify the particular HOD model space. Alternatively, it is possible that further generalizations of the HOD model would open opportunities for the squeezed 3PCF to provide novel parameter measurements. The GRAND-HOD Python package is publicly available at https://github.com/SandyYuan/GRAND-HOD .

High peaks in weak lensing (WL) maps originate dominantly from the lensing effects of single massive halos. Their abundance is therefore closely related to the halo mass function and thus a powerful cosmological probe. On the other hand, however, besides individual massive halos, large-scale structures (LSS) along lines of sight also contribute to the peak signals. In this paper, with ray tracing simulations, we investigate the LSS projection effects. We show that for current surveys with a large shape noise, the stochastic LSS effects are subdominant. For future WL surveys with source galaxies having a median redshift $z_{\mathrm{med}}\sim1$ or higher, however, they are significant. For the cosmological constraints derived from observed WL high peak counts, severe biases can occur if the LSS effects are not taken into account properly. We extend the model of \citet{Fan2010} by incorporating the LSS projection effects into the theoretical considerations. By comparing with simulation results, we demonstrate the good performance of the improved model and its applicability in cosmological studies.

We analyze the Sun's shadow observed with the Tibet-III air shower array and find that the shadow's center deviates northward (southward) from the optical solar disc center in the "Away" ("Toward") IMF sector. By comparing with numerical simulations based on the solar magnetic field model, we find that the average IMF strength in the "Away" ("Toward") sector is $1.54 \pm 0.21_{\rm stat} \pm 0.20_{\rm syst}$ ($1.62 \pm 0.15_{\rm stat} \pm 0.22_{\rm syst}$) times larger than the model prediction. These demonstrate that the observed Sun's shadow is a useful tool for the quantitative evaluation of the average solar magnetic field.

We report the result of a faint quasar survey in a one square degree field. The aim is to test the Y-K/g-z and J-K/i-Y color selection criteria for quasars at faint magnitude, to obtain a complete sample of quasars based on deep optical and near-infrared color-color selection, and to measure the faint end of quasar luminosity function (QLF) over a wide redshift range. We carried out a quasar survey based on the Y-K/g-z and J-K/i-Y quasar selection criteria, using the deep Y-band data obtained from our CFHT/WIRCam Y-band images in a two-degree field within the F22 field of the VIMOS VLT deep survey, optical co-added data from Sloan Digital Sky Survey Stripe 82 and deep near-infrared data from the UKIDSS Deep Extragalactic Survey in the same field. We discovered 25 new quasars at 0.5 < z < 4.5 and i < 22.5 mag within one square degree field. The survey significantly increases the number of faint quasars in this field, especially at z ~ 2-3. It confirms that our color selections are highly complete in a wide redshift range (z < 4.5), especially over the quasar number density peak at z ~ 2-3, even for faint quasars. Combining all previous known quasars and new discoveries, we construct a sample with 109 quasars, and measure the binned QLF and parametric QLF. Although the sample is small, our results agree with a pure luminosity evolution at lower redshift and luminosity evolution and density evolution model at redshift z > 2.5.

We investigate the three-point correlation function (3PCF) in the squeezed limit by considering galaxy pairs as discrete objects and cross-correlating them with the galaxy field. We develop an efficient algorithm using Fast Fourier Transforms to compute such cross-correlations and their associated pair-galaxy bias bpg and the squeezed 3PCF coefficient Qeff. We implement our method using N-body cosmological simulations and a fiducial Halo Occupation Distribution (HOD) and present the results in both the real space and redshift space. In real space, we observe a peak in bpg and Qeff at pair separation of ~ 2 Mpc, attributed to the fact that galaxy pairs at 2 Mpc separation trace the most massive dark matter halos. We also see strong anisotropy in the bpg and Qeff signals that track the large-scale filamentary structure. In redshift space, both the 2 Mpc peak and the anisotropy are significantly smeared out along the line-of-sight due to Finger-of-God effect. In both the real space and redshift space, the squeezed 3PCF shows a factor of 2 variation, contradicting the hierarchical ansatz but offering rich information on the galaxy-halo connection. Thus, we explore the possibility of using the squeezed 3PCF to constrain the HOD. When we compare two simple HOD models that are closely matched in their projected two-point correlation function (2PCF), we do not yet see a strong variation in the 3PCF that is clearly disentangled from variations in the projected 2PCF. Nevertheless, we propose that more complicated HOD models, e.g. those incorporating assembly bias, can break degeneracies in the 2PCF and show a distinguishable squeezed 3PCF signal.

Constraining neutrino mass remains an elusive challenge in modern physics. Precision measurements are expected from several upcoming cosmological probes of large-scale structure. Achieving this goal relies on an equal level of precision from theoretical predictions of neutrino clustering. Numerical simulations of the non-linear evolution of cold dark matter and neutrinos play a pivotal role in this process. We incorporate neutrinos into the cosmological N-body code CUBEP3M and discuss the challenges associated with pushing to the extreme scales demanded by the neutrino problem. We highlight code optimizations made to exploit modern high performance computing architectures and present a novel method of data compression that reduces the phase-space particle footprint from 24 bytes in single precision to roughly 9 bytes. We scale the neutrino problem to the Tianhe-2 supercomputer and provide details of our production run, named TianNu, which uses 86% of the machine (13,824 compute nodes). With a total of 2.97 trillion particles, TianNu is currently the world's largest cosmological N-body simulation and improves upon previous neutrino simulations by two orders of magnitude in scale. We finish with a discussion of the unanticipated computational challenges that were encountered during the TianNu runtime.

The influence of adiabatic focusing on particle diffusion is an important topic in astrophysics and plasma physics. In the past several authors have explored the influence of along-field adiabatic focusing on parallel diffusion of charged energetic particles. In this paper by using the Unified NonLinear Transport (UNLT) theory developed by Shalchi (SH2010) and the method of He and Schlickeiser (HS2014) we derive a new nonlinear perpendicular diffusion coefficient for non-uniform background magnetic field. This formula demonstrates that particle perpendicular diffusion coefficient is modified by along-field adiabatic focusing. For isotropic pitch-angle scattering and weak adiabatic focusing limit the derived perpendicular diffusion coefficient is independent of the sign of adiabatic focusing characteristic length. For two-component model we simplify the perpendicular diffusion coefficient up to second order of the power series of adiabatic focusing characteristic quantity. We find that the first order modifying factor is equal to zero and the sign of the second one is determined by the energy of particles.

The field line random walk (FLRW) of magnetic turbulence is one of the important topics in plasma physics and astrophysics. In this article by using the field line tracing method mean square displacements (MSD) of FLRW is calculated in all possible length scales for pure two-dimensional turbulence with damping dynamical model. We demonstrate that in order to describe FLRW with damping dynamical model a new dimensionless quantity $R$ is needed to be introduced. In different length scales dimensionless MSD shows different relationship with the dimensionless quantity $R$. Although temporal effect impacts MSD of FLRW and even changes regimes of FLRW, it dose not affect the relationship between the dimensionless MSD and dimensionless quantity $R$ in all possible length scales.

The physics and demographics of type 2 quasars remain poorly understood, and new samples of such objects selected in a variety of ways can give insight into their physical properties, evolution, and relationship to their host galaxies. We present a sample of 2758 type 2 quasars at z $\leq$ 1 from the SDSS-III/BOSS spectroscopic database, selected on the basis of their emission-line properties. We probe the luminous end of the population by requiring the rest-frame equivalent width of [OIII] to be > 100 {\AA}. We distinguish our objects from star-forming galaxies and type 1 quasars using line widths, standard emission line ratio diagnostic diagrams at z < 0.52 and detection of [Ne V]{\lambda}3426{\AA} at z > 0.52. The majority of our objects have [OIII] luminosities in the range 10^8.5-10^10 L$_{\odot}$ and redshifts between 0.4 and 0.65. Our sample includes over 400 type 2 quasars with incorrectly measured redshifts in the BOSS database; such objects often show kinematic substructure or outflows in the [OIII] line. The majority of the sample has counterparts in the WISE survey, with median infrared luminosity {\nu}L{\nu}[12{\mu}m] = 4.2 x 10^44 erg/sec. Only 34 per cent of the newly identified type 2 quasars would be selected by infrared color cuts designed to identify obscured active nuclei, highlighting the difficulty of identifying complete samples of type 2 quasars. We make public the multi-Gaussian decompositions of all [OIII] profiles for the new sample and for 568 type 2 quasars from SDSS I/II, together with non-parametric measures of line profile shapes and identify over 600 candidate double-peaked [OIII] profiles.

Measurements of neutrino mass in cosmological observations rely on two point statistics that are hindered by significant degeneracies with the optical depth and galaxy bias. The relative velocity effect between cold dark matter and neutrinos induces a large scale dipole into the matter density field and may be able to provide orthogonal constraints to standard techniques. We numerically investigate this dipole in the TianNu Simulation, which contains cold dark matter and 50 meV neutrinos. We first compute the dipole using a new linear response technique where we treat the displacement caused by the relative velocity as a phase in Fourier space and then integrate the matter power spectrum over redshift. Then, we compute the dipole numerically in real space using the simulation density and velocity fields. We find excellent agreement between the linear response and N-body methods. Utilizing the dipole as an observational tool will require two tracers of the matter distribution that are differently biased with respect to the neutrino density.

Astrophysical techniques have pioneered the discovery of neutrino mass properties. Current cosmological observations give an upper bound on neutrino masses by attempting to disentangle the small neutrino contribution from the sum of all matter using precise theoretical models. We discover the differential neutrino condensation effect in our TianNu N-body simulation. Neutrino masses can be inferred using this effect by comparing galaxy properties in regions of the universe with different neutrino relative abundance (i.e. the local neutrino to cold dark matter density ratio). In "neutrino-rich"' regions, more neutrinos can be captured by massive halos compared to "neutrino-poor" regions. This effect differentially skews the halo mass function and opens up the path to independent neutrino mass measurements in current or future galaxy surveys.

In this Letter, we report the observational constraints on the Hu-Sawicki $f(R)$ theory derived from weak lensing peak abundances, which are closely related to the mass function of massive halos. In comparison with studies using optical or x-ray clusters of galaxies, weak lensing peak analyses have the advantages of not relying on mass-baryonic observable calibrations. With observations from the Canada-France-Hawaii-Telescope Lensing Survey, our peak analyses give rise to a tight constraint on the model parameter $|f_{R0}|$ for $n=1$. The $95\%$ CL limit is $\log_{10}|f_{R0}| < -4.82$ given WMAP9 priors on $(\Omega_{\rm m}, A_{\rm s})$. With Planck15 priors, the corresponding result is $\log_{10}|f_{R0}| < -5.16$.

Utilizing the high-resolution, large-scale LAOZI cosmological simulations we investigate the nature of the metal-poor (${\rm [Z/H]<-2}$) damped Lyman alpha systems (mpDLA) at $z=3$. The following physical picture of mpDLAs emerges. The majority of mpDLAs inhabit regions $\ge 20$~kpc from the host galaxy center on infalling cold gas streams originating from the intergalactic medium, with infall velocity of $\sim 100$ km/s and temperature of $\sim 10^{4}$ K. For each host galaxy, on average, about $1\%$ of the area within a radius $150$~kpc is covered by mpDLAs. The mpDLAs are relatively diffuse ($n_{\rm{gas}} \sim 10^{-2}$ cm$^{-3}$), Jeans quasi-stable, and have very low star formation rate ($\dot{\Sigma} \le 10^{-4} M_{\odot} \rm{\ yr}^{-1} \rm{\ kpc}^{-2}$). As mpDLAs migrate inward to the galaxy center, they mix with high metallicity gas and stellar outflows in the process, removing themselves from the metal-poor category and rendering the central ($\le 5$ kpc) regions of galaxies devoid of mpDLAs. Thus, the central regions of the host galaxies are populated by mostly metal-rich DLAs instead of mpDLAs. All observables of the simulated mpDLAs are in excellent agreement with observations, except the gas density, which is about a factor of ten lower than the value inferred observationally. However, the observationally inferred value is based on simplified assumptions that are not borne out in the simulations.

We propose a valid scheme to measure the Hubble parameter $H(z)$ at high redshifts by detecting the Sandage-Loeb signal (SL signal) which can be realized by the next generation extremely large telescope. It will largely extend the current observational Hubble parameter data (OHD) towards the redshift region of $z \in [2.0,5.0]$, the so-called "redshift desert", where other dark energy probes are hard to provide useful information of the cosmic expansion. Quantifying the ability of this future measurement by simulating observational data for a CODEX (COsmic Dynamics and EXo-earth experiment)-like survey and constraining various cosmological models, we find that the SL signal scheme brings the redshift upper-limit of OHD from $z_\mathrm{max}=2.3$ to $z_\mathrm{max}\simeq 5.0$, provides more accurate constraints on different dark energy models, and greatly changes the degeneracy direction of the parameters. For the $\Lambda$CDM case, the accuracy of $\Omega_m$ is improved by $58\%$ and the degeneracy between $\Omega_m$ and $\Omega_ {\Lambda}$ is rotated to the vertical direction of $\Omega_k = 0$ line strongly; for the $w$CDM case, the accuracy of $w$ is improved by $15\%$. The Fisher matrix forecast on different time-dependent $w(z)$ is also performed.

By adopting the differential age method, we utilize selected 17832 luminous red galaxies (LRGs) from Sloan Digital Sky Survey Data Release Seven (SDSS DR7) covering redshift $0<z<0.4$ to measure Hubble parameters. Using a full spectrum fitting package UlySS, these spectra are reduced with single stellar population (SSP) models and optimal age information of our selected sample are derived. With the decreasing age-redshift relation, four new observational $H(z)$ data (OHD) points are obtained, which are $H(z)=69.0\pm19.6$ km s$^{-1}$ Mpc$^{-1}$ at $z=0.07$, $H(z)=68.6\pm26.2$ km s$^{-1}$ Mpc$^{-1}$ at $z=0.12$, $H(z)$=$72.9\pm29.6$ km s$^{-1}$ Mpc$^{-1}$ at $z=0.2$ and $H(z)$=$88.8\pm36.6$ km s$^{-1}$ Mpc$^{-1}$ at $z=0.28$, respectively. Combined with other 21 available OHD data points, a performance of constraint on both flat and non-flat $\Lambda$CDM model is presented.

Fluctuations of the 21 cm brightness temperature before the formation of the first stars hold the promise of becoming a high-precision cosmological probe in the future. The growth of over densities is very well described by perturbation theory at that epoch and the signal can in principle be predicted to arbitrary accuracy for given cosmological parameters. Recently, Tseliakhovich and Hirata pointed out a previously neglected and important physical effect, due to the fact that baryons and cold dark matter (CDM) have supersonic relative velocities after recombination. This relative velocity suppresses the growth of matter fluctuations on scales $k \sim 10-10^3$ Mpc$^{-1}$. In addition, the amplitude of the small-scale power spectrum is modulated on the large scales over which the relative velocity varies, corresponding to $k \sim 0.005-1$ Mpc$^{-1}$. In this paper, the effect of the relative velocity on 21 cm brightness temperature fluctuations from redshifts $z \geq 30$ is computed. We show that the 21 cm power spectrum is affected on ${\it most}$ scales. On small scales, the signal is typically suppressed several tens of percent, except for extremely small scales ($k \gtrsim 2000$ Mpc$^{-1}$) for which the fluctuations are boosted by resonant excitation of acoustic waves. On large scales, 21 cm fluctuations are enhanced due to the non-linear dependence of the brightness temperature on the underlying gas density and temperature. The enhancement of the 21 cm power spectrum is of a few percent at $k \sim 0.1$ Mpc$^{-1}$ and up to tens of percent at $k \lesssim 0.005$ Mpc$^{-1}$, for standard $\Lambda$CDM cosmology. In principle this effect allows to probe the small-scale matter power spectrum not only through a measurement of small angular scales but also through its effect on large angular scales.

We study the power of current and future observational Hubble parameter data (OHD) on non-parametric estimations of the dark energy equation of state, $w(z)$. We propose a new method by conjunction of principal component analysis (PCA) and the criterion of goodness of fit (GoF) criterion to reconstruct $w(z)$, ensuring the sensitivity and reliability of the extraction of features in the EoS. We also give an new error model to simulate future OHD data, to forecast the power of future OHD on the EoS reconstruction. The result shows that current OHD, despite in less quantity, give not only a similar power of reconstruction of dark energy compared to the result given by type Ia supernovae, but also extend the constraint on $w(z)$ up to redshift $z\simeq2$. Additionally, a reasonable forecast of future data in more quantity and better quality greatly enhances the reconstruction of dark energy.

Many schemes have been proposed to perform a model-independent constraint on cosmological dynamics, such as nonparametric dark energy equation of state (EoS) \omega(z) or the deceleration parameter q(z). These methods usually contain derivative processes with respect to observational data with noise. However, it still remains remarkably uncertain when one estimates the numerical differentiation, especially the corresponding truncation errors. In this work, we introduce a global numerical differentiation method, first formulated by Reinsch(1967), which is smoothed by cubic spline functions. The optimal solution is obtained by minimizing the functional \Phi(f). To investigate the potential of the algorithm further, we apply it to the estimation of the transition redshift z_{t} with simulated expansion rate E(z) based on observational Hubble parameter data(OHD). An effective method to determine the free parameter S appearing in Reinsch Splines is provided.

The recent spectacular progress in the experimental and theoretical understanding of graphene, the basic constituent of graphite, is applied here to compute, from first principles, the UV extinction of nano-particles made of stacks of graphene layers. The theory also covers cases where graphene is affected by structural, chemical or orientation disorder, each disorder type being quantitatively defined by a single parameter. The extinction bumps carried by such model materials are found to have positions and widths falling in the same range as the known astronomical 2175 \AA features: as the disorder parameter increases, the bump width increases from 0.85 to 2.5 $\mu$m$^{-1}$, while its peak position shifts from 4.65 to 4.75 $\mu$m$^{-1}$. Moderate degrees of disorder are enough to cover the range of widths of the vast majority of observed bumps (0.75 to 1.3 $\mu$m$^{-1}$). Higher degrees account for outliers, also observed in the sky. The introduction of structural or chemical disorder amounts to changing the initial $sp^{2}$ bondings into $sp^{3}$ or $sp^{1}$, so the optical properties of the model material become similar to those of the more or less amorphous carbon-rich materials studied in the laboratory: a-C, a-C:H, HAC, ACH, coals etc. The present treatment thus bridges gaps between physically different model materials.