Inflation predicts that the Universe is spatially flat. The Planck 2018 measurements of the cosmic microwave background anisotropy favour a spatially closed universe at more than 2$\sigma$ confidence level. We use model independent methods to study the issue of cosmic curvature. The method reconstructs the Hubble parameter $H(z)$ from cosmic chronometers data with the Gaussian process method. The distance modulus is then calculated with the reconstructed function $H(z)$ and fitted by type Ia supernovae data. Combining the cosmic chronometers and type Ia supernovae data, we obtain $\Omega_{k0}h^2=0.102\pm 0.066$ which shows a spatially open universe at more than 1$\sigma$ confidence level. By adding the redshift space distortions data to the type Ia supernovae data, we obtain $\Omega_{k0}h^2=0.117^{+0.058}_{-0.045}$ and no interesting deviation from $\Lambda$CDM model is found. The almost $2\sigma$ evidence of a spatially open universe is independent of cosmological model and gravitational theory.

Directly comparing the 6 expansion rate measured by type Ia supernovae data and the lower bound on the expansion rate set by the strong energy conditions or the null hypothesis that there never exists cosmic acceleration, we see $3\sigma$ direct evidence of cosmic acceleration and the $Rh=ct$ model is strongly excluded by the type Ia supernovae data. We also use Gaussian process method to reconstruct the expansion rate and the deceleration parameter from the 31 cosmic chronometers data and the 6 data points on the expansion rate measured from type Ia supernoave data, the direct evidence of cosmic acceleration is more than $3\sigma$ and we find that the transition redshift $z_t=0.60_{-0.12}^{+0.21}$ at which the expansion of the Universe underwent the transition from acceleration to deceleration. The Hubble constant inferred from the cosmic chronometers data with the Gaussian process method is $H_0=67.46\pm4.75$ Km/s/Mpc. To understand the properties of cosmic acceleration and dark energy, we fit two different two-parameter models to the observational data, and we find that the constraints on the model parameters from either the full distance modulus data by the Pantheon compilation or the compressed expansion rate data are very similar, and the derived Hubble constants are consistent with the Planck 2018 result. Our results confirm that the 6 compressed expansion rate data can replace the full 1048 distance modulus data from the Pantheon compilation. We derive the transition redshift $z_t=0.61_{-0.16}^{+0.24}$ by fitting a simple $q(z)$ model to the combination of cosmic chronometers data and the Pantheon compilation, the result is consistent with that obtained from the reconstruction with Gaussian process. Abridged.

Glycolamide is a glycine isomer and also one of the simplest derivatives of acetamide (e.g., one hydrogen atom is replaced with a hydroxyl group), which is a known interstellar molecule. Using a battery of state of the art rotational spectroscopic techniques in the frequency and time domain, around 1500 transitions have been newly assigned. Based on the reliable frequency predictions, we report a radioastronomical search for glycolamide in the well known high-mass star forming region Sgr B2(N) using the ALMA imaging spectral line survey ReMoCA. We also searched for glycolamide toward Sgr B2(N) with the Effelsberg radio telescope. We report the nondetection of glycolamide toward this source with an abundance at least six and five times lower than that of acetamide and glycolaldehyde, respectively. Our astrochemical model suggests that glycolamide may be present in this source at a level just below the upper limit, which was derived from the EMoCA survey. We could also not detect the molecule in the region's extended molecular envelope, which was probed with the Effelsberg telescope. We find an upper limit to its column density that is similar to the column densities obtained earlier for acetamide and glycolaldehyde with the Green Bank Telescope.

The discovery of gravitational waves, which are ripples of space-time itself, opened a new window to test general relativity, because it predicts that there are only plus and cross polarizations for gravitational waves. For alternative theories of gravity, there may be up to six polarizations. The measurement of the polarization is one of the major scientific goals for future gravitational wave detectors. To evaluate the capability of the detector, we need to use the frequency dependent response functions averaged over the source direction and polarization angle. We derive the full analytical formulas of the averaged response functions for all six possible polarizations and present their asymptotic behaviors based on these analytical formulas. Compared with the numerical simulation, the full analytical formulas are more efficient and valid for any equal-arm interferometric gravitational wave detector without optical cavities in the arms and for a time-delay-interferometry Michelson combination.

We present the results of a ~250 square degrees CO mapping (+26d<l<+50d and -5d<b<+5d) toward the Aquila Rift region at a spatial resolution of ~50" and a grid spacing of 30". The high dynamic range CO maps with a spectral resolution of ~0.2km/s display highly structured molecular cloud (MC) morphologies with valuable velocity information, revealing complex spatial and dynamical features of the local molecular gas. In combination with the MWISP CO data and the Gaia DR2, distances of the main MC structures in the local ISM are well determined toward the Aquila Rift. We find that the total MC mass within 1 kpc is about >4.1x10^5 Msun in the whole region. In fact, the mass of the molecular gas is dominated by the W40 giant molecular cloud (GMC) at ~474 pc (~1.4x10^5 Msun) and the GMC complex G036.0+01.0 at ~560-670 pc (~2.0x10^5 Msun), while the MCs at ~220-260pc have gas masses of 10^2-10^3 Msun. Interestingly, an ~80pc long filamentary MC G044.0-02.5 at a distance of ~404 pc shows a systematic velocity gradient along and perpendicular to the major axis of the filament. The HI gas with the enhanced emission has the similar spatial morphologies and velocity features compared to the corresponding CO structure, indicating that the large-scale converging HI flows are probably responsible for the formation of the MC. Meanwhile, the long filamentary MC consists of many sub-filaments with the lengths ranging from ~0.5 pc to several pc, as well as prevalent networks of filaments in other large-scale local MCs.

The possibility that in the mass range around $10^{-12}\ M_\odot$ most of dark matter constitutes of primordial black holes (PBHs) is a very interesting topic. To produce PBHs with this mass, the primordial scalar power spectrum needs to be enhanced to the order of 0.01 at the scale $k\sim 10^{12}\ \text{Mpc}^{-1}$. The enhanced power spectrum also produces large secondary gravitational waves at the mHz band. A phenomenological delta function power spectrum is usually used to discuss the production of PBHs and secondary gravitational waves. Based on G and k inflations, we propose a new mechanism to enhance the power spectrum at small scales by introducing a non-canonical kinetic term $[1-2G(\phi)]X$ with the function $G(\phi)$ having a peak. Away from the peak, $G(\phi)$ is negligible and we recover the usual slow-roll inflation which is constrained by the cosmic microwave background anisotrpy observations. Around the peak, the slow-roll inflation transiently turns to ultra slow-roll inflation. The enhancement of the power spectrum can be obtained with generic potentials, and there is no need to fine tune the parameters in $G(\phi)$. The energy spectrum $\Omega_{GW}(f)$ of secondary gravitational waves have the characteristic power law behaviour $\Omega_{GW}(f)\sim f^{n}$ and is testable by pulsar timing array and space based gravitational wave detectors.

Understanding the formation and evolution of the Universe is crucial for cosmological studies, and the line intensity mapping provides a powerful tool for this kind of study. We propose to make use of multipole moments of redshift-space line intensity power spectrum to constrain the cosmological and astrophysical parameters, such as the equation of state of dark energy, massive neutrinos, primordial non-Gaussianity, and star formation rate density. As an example, we generate mock data of multipole power spectra for H-alpha 6563AA, [OIII] 5007AA and [OII] 3727AA measured by SPHEREx experiment at z=1 considering contaminations from interloper lines, and use Markov Chain Monte Carlo (MCMC) method to constrain the parameters in the model. We find a good fitting result of the parameters compared to their fiducial values, which means that the multipole power spectrum can effectively distinguish signal and interloper lines, and break the degeneracies between parameters, such as line mean intensity and bias. We also explore the cross power spectrum with CSST (Chinese Space Station Telescope) spectroscopic galaxy survey in the constraints. Since more accurate fitting results can be obtained by including measurements of the emission lines at higher redshifts out to z=3 at least and cross-correlations between emission lines can be involved, the line intensity mapping is expected to offer excellent results in future cosmological and astrophysical studies.

The inclination distribution of circumbinary planets (CBPs) is an important scientific issue. It is of great significance in estimating the occurrence rate of CBPs and studying their formation and evolution. Although the CBPs currently discovered by the transit method are nearly coplanar, the true inclination distribution of CBPs is still unknown. Previous researches on CBPs mostly regarded them as an isolated binary-planet system, without considering the birth environment of their host binaries. It is generally believed that almost all stars are born in clusters. Therefore, it is necessary to consider the impact of the close encounters of stars on the CBP systems. This article discusses how the close encounters of fly-by stars affect the inclination of CBP. Based on extensive numerical simulations, we found that CBPs in close binary with a spacing of about 0.2 au are almost unaffected by star fly-bys. Their orbits remain coplanar. However, when the spacing of the binary stars is greater than 1 au, the 2-3 fly-bys of the intruding star can excite a considerable inclination even for the CBP near the unstable boundary of the binary. For the planets in the outer region, a single star fly-by can excite an inclination to more than 5 degrees. Especially, CBPs in near polar or retrograde orbits can naturally form through bianry-star encounters. If close binaries are born in open clusters, our simulations suggest that there may be high-inclination CBPs in binaries with spacing > 1 au.

The cosmic infrared background (CIB) anisotropies and cosmic microwave background (CMB) lensing are powerful measurements for exploring cosmological and astrophysical problems. In this work, we measure the auto-correlation power spectrum of the CIB anisotropies in the Herschel-SPIRE HerMES Large Mode Survey (HeLMS) field, which covers more than 280 square degrees of the sky, and the cross power spectrum between HeLMS CIB and the CMB lensing from Planck. We use the Herschel Level 1 timeline data to merge the CIB maps at 250, 350 and 500 um bands, and mask the areas where the flux is greater than 3-sigma (~50 mJy/beam) or no measured data. We obtain the final CIB power spectra at 100<ell<20,000 by considering several effects, such as beam function, mode coupling, transfer function, and so on. We also calculate the theoretical CIB auto- and cross-power spectra of CIB and CMB lensing by assuming that the CIB emissivity follows Gaussian distribution in redshift. The Markov chain Monte Carlo (MCMC) method is used to perform a joint fit of the CIB auto-power spectra and cross-power spectra of CIB and CMB lensing data. We find that our model can fit the power spectrum data very well, and obtain basically consistent results with higher accuracy comparing to previous studies.

Ultra-light axions (ULAs) with mass less than 10^-20 eV have interesting behaviors that may contribute to either dark energy or dark matter at different epochs of the Universe. Its properties can be explored by cosmological observations, such as expansion history of the Universe, cosmic large-scale structure, cosmic microwave background, etc. In this work, we study the ULAs with a mass around 10^-33 eV, which means the ULA field still rolls slowly at present with the equation of state w=-1 as dark energy. In order to investigate the mass and other properties of this kind of ULA field, we adopt the measurements of Type Ia supernova (SN Ia), baryon acoustic oscillation (BAO), and Hubble parameter H(z). The Markov Chain Monte Carlo (MCMC) technique is employed to perform the constraints on the parameters. Finally, by exploring four cases of the model, we find that the mass of this ULA field is about 3x10^-33 eV if assuming the initial axion field phi_i=M_pl. We also investigate a general case by assuming phi_i< M_pl and find that the fitting results of phi_i/M_pl are consistent with or close to 1 for the datasets we use.

Primordial black holes and secondary gravitational waves can be used to probe the small scale physics at very early time. For secondary gravitational waves produced after the horizon reentry, we derive an analytical formula for the time integral of the source and analytical behavior of the time dependence of the energy density of induced gravitational waves is obtained. By proposing a piecewise power law parametrization for the power spectrum of primordial curvature perturbations, we use the observational constraints on primordial black hole dark matter to obtain an upper bound on the power spectrum, and discuss the test of the model with future space based gravitational wave antenna.

Methanol (CH3OH) is one of the most abundant interstellar molecules, offering a vast number of transitions to be studied, including many maser lines. While the strongest Galactic CH3OH lines, the so-called class II masers, show no indications for the presence of superluminous counterparts in external galaxies, the less luminous Galactic class I sources appear to be different. Here we report class I 36GHz CH3OH 4(-1) - 3(0) E line emission from the nearby galaxies Maffei2 and IC342, measured with the 100-m telescope at Effelsberg at three different epochs within a time span of about five weeks. The 36GHz methanol line of Maffei2 is the second most luminous among the sources detected with certainty outside the Local Group of galaxies. This is not matched by the moderate infrared luminosity of Maffei2. Higher resolution data are required to check whether this is related to its prominent bar and associated shocks. Upper limits for M82, NGC4388, NGC5728 and Arp220 are also presented. The previously reported detection of 36GHz maser emission in Arp220 is not confirmed. Non-detections are reported from the related class I 44GHz methanol transition towards Maffei2 and IC342, indicating that this line is not stronger than its 36GHz counterpart. In contrast to the previously detected 36GHz CH3OH emission in NGC253 and NGC4945, our 36GHz profiles towards Maffei2 and IC342 are similar to those of previously detected non-masing lines from other molecular species. However, by analogy to our Galactic center region, it may well be possible that the 36GHz methanol lines in Maffei~2 and IC~342 are composed of a large number of faint and narrow maser features that remain spatially unresolved. In view of this, a search for a weak broad 36GHz line component would also be desirable in NGC253 and NGC4945.

In order to search for further observational evidence of cloud-cloud collisions in one of the promising candidates, L1188, we carried out observations of multiple molecular lines toward the intersection region of the two nearly orthogonal filamentary molecular clouds in L1188. Based on these observations, we find two parallel filamentary structures, both of which have at least two velocity components being connected with broad bridging features. We also found a spatially complementary distribution between the two molecular clouds, as well as enhanced $^{13}$CO emission and $^{12}$CO self-absorption toward their abutting regions. At the most blueshifted velocities, we unveil a 1~pc-long arc ubiquitously showing $^{12}$CO line wings. We discover two 22 GHz water masers, which are the first maser detections in L1188. An analysis of line ratios at a linear resolution of 0.2 pc suggests that L1188 is characterised by kinetic temperatures of 13--23~K and H$_{2}$ number densities of 10$^{3}$--10$^{3.6}$ cm$^{-3}$. On the basis of previous theoretical predictions and simulations, we suggest that these observational features can be naturally explained by the scenario of a cloud-cloud collision in L1188, although an additional contribution of stellar feedback from low-mass young stellar objects cannot be ruled out.

The precise measurement of the spectrum of protons, the most abundant component of the cosmic radiation, is necessary to understand the source and acceleration of cosmic rays in the Milky Way. This work reports the measurement of the cosmic ray proton fluxes with kinetic energies from 40 GeV to 100 TeV, with two and a half years of data recorded by the DArk Matter Particle Explorer (DAMPE). This is the first time an experiment directly measures the cosmic ray protons up to ~100 TeV with a high statistics. The measured spectrum confirms the spectral hardening found by previous experiments and reveals a softening at ~13.6 TeV, with the spectral index changing from ~2.60 to ~2.85. Our result suggests the existence of a new spectral feature of cosmic rays at energies lower than the so-called knee, and sheds new light on the origin of Galactic cosmic rays.

Space-based gravitational wave detectors cannot keep rigid structures and precise arm length equality, so the precise equality of detector arms which is required in a ground-based interferometer to cancel the overwhelming laser noise is impossible. The time-delay interferometry method is applied to unequal arm lengths to cancel the laser frequency noise. We give analytical formulas of the averaged response functions for tensor, vector, breathing and longitudinal polarizations in different TDI combinations, and obtain their asymptotic behaviors. At low frequencies, $f\ll f_*$, the averaged response functions of all TDI combinations increase as $f^2$ for all six polarizations. The one exception is the averaged response functions of $\zeta$ for all six polarizations increase as $f^4$ in the equilateral-triangle case. At high frequencies, $f\gg f_*$, the averaged response functions of all TDI combinations for the tensor and breathing modes fall off as $1/f^2$, the averaged response functions of all TDI combinations for the vector mode fall off as $\ln(f)/f^2$ , and the averaged response functions of all TDI combinations for the longitudinal mode fall as $1/f$. We also give LISA and TianQin sensitivity curves in different TDI combinations for tensor, vector, breathing and longitudinal polarizations.

The Chinese Space Station Optical Survey (CSS-OS) is a planned full sky survey operated by the Chinese Space Station Telescope (CSST). It can simultaneously perform the photometric imaging and spectroscopic slitless surveys, and will probe weak and strong gravitational lensing, galaxy clustering, individual galaxies and galaxy clusters, active galactic nucleus (AGNs), and so on. It aims to explore the properties of dark matter and dark energy and other important cosmological problems. In this work, we focus on two main CSS-OS scientific goals, i.e. the weak gravitational lensing (WL) and galaxy clustering surveys. We generate the mock CSS-OS data based on the observational COSMOS and zCOSMOS catalogs. We investigate the constraints on the cosmological parameters from the CSS-OS using the Markov Chain Monte Carlo (MCMC) method. The intrinsic alignments, galaxy bias, velocity dispersion, and systematics from instrumental effects in the CSST WL and galaxy clustering surveys are also included, and their impacts on the constraint results are discussed. We find that the CSS-OS can improve the constraints on the cosmological parameters by a factor of a few (even one order of magnitude in the optimistic case), compared to the current WL and galaxy clustering surveys. The constraints can be further enhanced when performing joint analysis with the WL, galaxy clustering, and galaxy-galaxy lensing data. Therefore, the CSS-OS is expected to be a powerful survey for exploring the Universe. Since some assumptions may be still optimistic and simple, it is possible that the results from the real survey could be worse. We will study these issues in details with the help of simulations in the future.

For non-spinning, charged (Reissner-Nordstr\"om) black holes, the particles with an opposite sign of charge with respect to that of the black hole will be pulled into the black hole by the extra electromagnetic force. Such a hole will be quickly neutralized so that there should not exist significantly charged, non-spinning black holes in the universe. The case of spinning, charged (Kerr-Newmann) black holes is more complicated. For a given initial position and initial velocity of the particle, an oppositely charged particle does not always more easily fall into the black hole than a neutral particle. The possible existence of a magnetosphere further complicate the picture. One therefore cannot straightforwardly conclude that a charged spinning black hole will be neutralized. In this paper, we make the first step to investigate the neutralization of Kerr-Newmann black holes without introducing a magnetosphere. We track the particle trajectories under the influence of the curved spacetime and the electromagnetic field carried by the spinning, charged black hole. A statistical method is used to investigate the neutralization problem. We find a universal dependence of the falling probability into the black hole on the charge of the test particle, with the oppositely charged particles having a higher probability of falling. We therefore conclude that charged, spinning black holes without a magnetosphere should be quickly neutralized, consistent with people's intuition. The neutralization problem of Kerr-Newmann black holes with a co-rotating force-free magnetosphere is subject to further studies.

The DArk Matter Particle Explorer (DAMPE), a satellite-based cosmic ray and gamma-ray detector, was launched on December 17, 2015, and began its on-orbit operation on December 24, 2015. In this work we document the on-orbit calibration procedures used by DAMPE and report the calibration results of the Plastic Scintillator strip Detector (PSD), the Silicon-Tungsten tracKer-converter (STK), the BGO imaging calorimeter (BGO), and the Neutron Detector (NUD). The results are obtained using Galactic cosmic rays, bright known GeV gamma-ray sources, and charge injection into the front-end electronics of each sub-detector. The determination of the boundary of the South Atlantic Anomaly (SAA), the measurement of the live time, and the alignments of the detectors are also introduced. The calibration results demonstrate the stability of the detectors in almost two years of the on-orbit operation.

We derive molecular-gas-phase $^{12}$C/$^{13}$C isotope ratios for the central few 100 pc of the three nearby starburst galaxies NGC 253, NGC 1068, and NGC 4945 making use of the $\lambda$ $\sim$ 3 mm $^{12}$CN and $^{13}$CN $N$ = 1--0 lines in the ALMA Band 3. The $^{12}$C/$^{13}$C isotopic ratios derived from the ratios of these lines range from 30 to 67 with an average of 41.6 $\pm$ 0.2 in NGC 253, from 24 to 62 with an average of 38.3 $\pm$ 0.4 in NGC 1068, and from 6 to 44 with an average of 16.9 $\pm$ 0.3 in NGC 4945. The highest $^{12}$C/$^{13}$C isotopic ratios are determined in some of the outskirts of the nuclear regions of the three starburst galaxies. The lowest ratios are associated with the northeastern and southwestern molecular peaks of NGC 253, the northeastern and southwestern edge of the mapped region in NGC 1068, and the very center of NGC 4945. In case of NGC 1068, the measured ratios suggest inflow from the outer part of NGC 1068 into the circum-nuclear disk through both the halo and the bar. Low $^{12}$C/$^{13}$C isotopic ratios in the central regions of these starburst galaxies indicate the presence of highly processed material.

The GLOSTAR survey will study the star formation in the Galactic plane between $-2^{\circ}< \textit{l}< 85^{\circ}$ and $|b| < 1^{\circ}$ with unprecedented sensitivity in both, flux density ($\sim$ 40 $\mu Jy$ beam$^{-1}$) and range of angular scales ($\sim$ 1."5 to the largest radio structures in the Galaxy). In this paper we present the first results obtained from a radio continuum map of a 16 square degree sized region of the Galactic plane centered on $\textit{l} = 32^{\circ}$ and $b = 0^{\circ}$ ($28^{\circ} < \textit{l} < 36^{\circ}$ and $|b| < 1^{\circ}$). This map has a resolution of 18" and sensitivity of $\sim$ 60-150 $\mu Jy$ beam$^{-1}$. We present data acquired with the VLA in D-configuration. Two 1 GHz wide sub-bands were observed simultaneously and centred at 4.7 and 6.9 GHz. These data were calibrated and imaged using the $\textit{Obit}$ software package. The source extraction has been performed using the BLOBCAT software package and verified through a combination of visual inspection and cross-matching with other radio and mid-infrared surveys. The final catalog consists of 1575 discrete radio sources and 27 large scale structures (including W43 and W44). By cross-matching with other catalogs and calculating the spectral indices ($S(\nu) \propto \nu^\alpha$), we have classified 231 continuum sources as HII regions, 37 as ionization fronts, and 46 as planetary nebulae. The longitude and latitude distribution and negative spectral indices are all consistent with the vast majority of the unclassified sources being extragalactic background sources. We present a catalog of 1575 radio continuum sources and discuss their physical properties, emission nature and relation with previously reported. These first GLOSTAR results have increased the number of reliable HII regions in this part of the Galaxy by a factor of four.

We use the cross-correlation between the thermal Sunyaev-Zeldovich (tSZ) signal measured by the Planck satellite and the luminous red galaxy (LRG) samples provided by the SDSS DR7 to study the properties of galaxy cluster and intra-cluster gas. We separate the samples into three redshift bins z_1=(0.16, 0.26), z_2=(0.26, 0.36), z_3=(0.36, 0.47), and stack the Planck y-map against LRGs to derive the averaged y-profile for each redshift bin. We then fit the stacked profile with the theoretical prediction from the universal pressure profile (UPP) by using the Markov-Chain Monte-Carlo method. We find that the best-fit values of the UPP parameters for the three bins are generally consistent with the previous studies, except for the noticeable evolution of the parameters in the three redshift bins. We simultaneously fit the data in the three redshift bins together, and find that the original UPP model cannot fit the data at small angular scales very well in the first and third redshift bins. The joint fits can be improved by including an additional parameter eta to change the redshift-dependence of the model (i.e. E(z)^{8/3} -> E(z)^{8/3+eta}) with best-fit value as eta=-3.11^{+1.09}_{-1.13}. This suggests that the original UPP model with less redshift-dependence may provide a better fit to the stacked thermal Sunyaev-Zeldovich profile.

Gravitational waves are perturbations of the metric of space-time. Six polarizations are possible, although general relativity predicts that only two such polarizations, tensor plus and tensor cross are present for gravitational waves. We give the analytical formulas for the antenna response functions for the six polarizations which are valid for any equal-arm interferometric gravitational-wave detectors without optical cavities in the arms.The response function averaged over the source direction and polarization angle decreases at high frequencies which deteriorates the signal-to-noise ratio registered in the detector. At high frequencies, the averaged response functions for the tensor and breathing modes fall of as $1/f^2$, the averaged response function for the longitudinal mode falls off as $1/f$ and the averaged response function for the vector mode falls off as $\ln(f)/f^2$.

Different inflationary models predict oscillatory features in the primordial power spectrum. These can leave an imprint on both the cosmic microwave background (CMB) and the large-scale structure (LSS) of our Universe, that can be searched for with current data. Inspired by the axion-monodromy model of inflation, we search for primordial oscillations that are logarithmic in wavenumber, using both CMB data from the Planck satellite and LSS data from the WiggleZ galaxy survey. We find that, within our search range for the new oscillation parameters (amplitude, frequency and phase), both CMB-only and CMB+LSS data yield the same best-fit oscillation frequency of $\log_{10}\omega=1.5$, improving the fit over $\Lambda$CDM by $\Delta\chi^2\! =\! -9$ and $\Delta\chi^2\! =\! -13$ (roughly corresponding to $2\sigma$ and $2.8\sigma$ significance), respectively. Bayesian evidence for the log-oscillation model versus $\Lambda$CDM indicates a very slight preference for the latter. Future CMB and LSS data will further probe this scenario.

The Milky Way Imaging Scroll Painting (MWISP) project is an unbiased Galactic plane CO survey for mapping regions of l=-10deg to +250deg and |b|<5.2deg with the 13.7 m telescope of the Purple Mountain Observatory. The legacy survey aims to observe the 12CO, 13CO, and C18O (J=1-0) lines simultaneously with full sampling using the nine-beam Superconducting SpectroScopic Array Receiver (SSAR) system with an instantaneous bandwidth of 1 GHz. In this paper, the completed 250 deg^2 data from l=+25.8deg to +49.7deg are presented with a grid spacing of 30" and a typical rms noise level of ~0.5 K for 12CO at the channel width of 0.16 km/s and ~0.3 K for 13CO and C18O at 0.17 km/s. The high quality data with moderate resolution (50"), uniform sensitivity, and high spatial dynamic range, allow us to investigate the details of molecular clouds (MCs) traced by the three CO isotope lines. Three interesting examples are briefly investigated, including distant Galactic spiral arms traced by CO emission with VLSR<0 km/s, the bubble-like dense gas structure near the HII region W40, and the MCs distribution perpendicular to the Galactic plane.

We give the procedure to reconstruct the extended inflationary potentials for attractors and use the $\alpha$ attractor and the constant-roll model as examples to show how to reconstruct the class of extended inflationary potentials in the strong coupling limit. We also derive the strong coupling condition which the coupling constant $\xi$ satisfies. The class of extended inflationary potentials has the same attractor, and the reconstructed extended inflationary potentials are consistent with the observational constraints. The existence of the attractors from extended inflationary models further challenges the model discrimination just by the observables $n_s$ and $r$.

The two swampland criteria are generically in tension with the single field slow-roll inflation because the first swampland criterion requires small tensor to scalar ratio while the second swampland criterion requires large tensor to scalar ratio. The challenge to the single field slow-roll inflation imposed by the swampland criteria can be avoided by modifying the relationship between the tensor to scalar ratio and the slow-roll parameter. We show that the Gauss-Bonnet inflation with the coupling function inversely proportional to the potential overcomes the challenge by adding a constant factor in the relationship between the tensor to scalar ratio and the slow-roll parameter. For the Gauss-Bonnet inflation, while the swampland criteria are satisfied, the slow-roll conditions are also fulfilled, so the scalar spectral tilt and the tensor to scalar ratio are consistent with the observations. We use the potentials for chaotic inflation and the E-model as examples to show that the models pass all the constraints. The swampland criteria may imply Gauss-Bonnet coupling.

The Serpens filament, as one of the nearest infrared dark clouds, is regarded as a pristine filament at a very early evolutionary stage of star formation. In order to study its molecular content and dynamical state, we mapped this filament in seven species. Among them, HCO$^{+}$, HNC, HCN, and CS show self-absorption, while C$^{18}$O is most sensitive to the filamentary structure. A kinematic analysis demonstrates that this filament forms a velocity-coherent (trans-)sonic structure, a large part of which is one of the most quiescent regions in the Serpens cloud. Widespread C$^{18}$O depletion is found throughout the Serpens filament. Based on the Herschel dust-derived H$_{2}$ column density map, the line mass of the filament is 36--41~M$_{\odot}$~pc$^{-1}$, and its full width at half maximum is 0.17$\pm$0.01~pc, while its length is ~1.6~pc. The inner radial column density profile of this filament can be well fitted with a Plummer profile with an exponent of 2.2$\pm$0.1, a scale radius of $0.018\pm 0.003$ pc, and a central density of $(4.0\pm 0.8)\times 10^{4}$~cm$^{-3}$. The Serpens filament appears to be slightly supercritical. The widespread blue-skewed HNC and CS line profiles and HCN hyperfine line anomalies across this filament indicate radial infall in parts of the Serpens filament. C$^{18}$O velocity gradients also indicate accretion flows along the filament. The velocity and density structures suggest that such accretion flows are likely due to a longitudinal collapse parallel to the filament's long axis. Both the radial infall rate and the longitudinal accretion rate along the Serpens filament are lower than all previously reported values in other filaments. This indicates that the Serpens filament lies at an early evolutionary stage when collapse has just begun, or that thermal and non-thermal support are effective in providing support against gravity.

The speed of gravitational waves provides us a new tool to test alternative theories of gravity. The constraint on the speed of gravitational waves from GW170817 and GRB170817A is used to test some classes of Horndeski theory. In particular, we consider the coupling of a scalar field to Einstein tensor and the coupling of the Gauss-Bonnet term to a scalar field. The coupling strength of the Gauss-Bonnet coupling is constrained to be in the order of $10^{-15}$. In the Horndeski theory we show that in order for this theory to satisfy the stringent constraint on the speed of GWs the mass scale $M$ introduced in the non-minimally derivative coupling is constrained to be in the range $10^{15}\text{GeV}\gg M \gtrsim 2\times 10^{-35}$GeV taking also under consideration the early times upper bound for the mass scale $M$. The large mass ranges require no fine-tuning because the effect of non-minimally derivative coupling is negligible at late times.

The Chinese Space Station Optical Survey (CSS-OS) is a major science project of the Space Application System of the China Manned Space Program. This survey is planned to perform both photometric imaging and slitless spectroscopic observations, and it will focus on different cosmological and astronomical goals. Most of these goals are tightly dependent on the accuracy of photometric redshift (photo-z) measurement, especially for the weak gravitational lensing survey as a main science driver. In this work, we assess if the current filter definition can provide accurate photo-z measurement to meet the science requirement. We use the COSMOS galaxy catalog to create a mock catalog for the CSS-OS. We compare different photo-z codes and fitting methods that using the spectral energy distribution (SED) template-fitting technique, and choose to use a modified LePhare code in photo-z fitting process. Then we investigate the CSS-OS photo-z accuracy in certain ranges of filter parameters, such as band position, width, and slope. We find that the current CSS-OS filter definition can achieve reasonably good photo-z results with sigma_z~0.02 and outlier fraction ~3%.

The next generation of space-borne gravitational wave detectors may detect gravitational waves from extreme mass-ratio inspirals with primordial black holes. To produce primordial black holes which contribute a non-negligible abundance of dark matter and are consistent with the observations, a large enhancement in the primordial curvature power spectrum is needed. For a single field slow-roll inflation, the enhancement requires a very flat potential for the inflaton, and this will increase the number of $e$-folds. To avoid the problem, an ultra-slow-roll inflation at the near inflection point is required. We elaborate the conditions to successfully produce primordial black hole dark matter from single field inflation and propose a toy model with polynomial potential to realize the big enhancement of the curvature power spectrum at small scales while maintaining the consistency with the observations at large scales. The power spectrum for the second order gravitational waves generated by the large density perturbations at small scales is consistent with the current pulsar timing array observations.

Although several S-type and P-type planets in binary systems were discovered in past years, S-type planets have not yet been found in close binaries with an orbital separation not more than 5 au. Recent studies suggest that S-type planets in close binaries may be detected through high-accuracy observations. However, nowadays planet formation theories imply that it is difficult for S-type planets in close binaries systems to form in situ. In this work, we extensively perform numerical simulations to explore scenarios of planet-planet scattering among circumbinary planets and subsequent tidal capture in various binary configurations, to examine whether the mechanism can play a part in producing such kind of planets. Our results show that this mechanism is robust. The maximum capture probability is $\sim 10\%$, which can be comparable to the tidal capture probability of hot Jupiters in single star systems. The capture probability is related to binary configurations, where a smaller eccentricity or a low mass ratio of the binary will lead to a larger probability of capture, and vice versa. Furthermore, we find that S-type planets with retrograde orbits can be naturally produced via capture process. These planets on retrograde orbits can help us distinguish in situ formation and post-capture origin for S-type planet in close binaries systems. The forthcoming missions (PLATO) will provide the opportunity and feasibility to detect such planets. Our work provides several suggestions for selecting target binaries in search for S-type planets in the near future.

We analyze the polarization content of gravitational waves in Horndeski theory. Besides the familiar plus and cross polarizations in Einstein's General Relativity, there is one more polarization state which is the mixture of the transverse breathing and longitudinal polarizations.The additional mode is excited by the massive scalar field. In the massless limit, the longitudinal polarization disappears, while the breathing one persists. The upper bound on the graviton mass severely constrains the amplitude of the longitudinal polarization, which makes its detection highly unlikely by the ground-based or space-borne interferometers in the near future. However, pulsar timing arrays might be able to detect the polarization excited by the massive scalar field. Since additional polarization states appear in alternative theories of gravity, the measurement of the polarizations of gravitational waves can be used to probe the nature of gravity. In addition to the plus and cross states, the detection of the breathing polarization means that gravitation is mediated by massless spin 2 and spin 0 fields, and the detection of both the breathing and longitudinal states means that gravitation is propagated by the massless spin 2 and massive spin 0 fields.

For the constant-roll tachyon inflation, we derive the analytical expressions for the scalar and tensor power spectra, the scalar and tensor spectral tilts and the tensor to scalar ratio up to the first order by using the method of Bessel function approximation. The derived $n_s-r$ results for the constant-roll inflation are also compared with the observations, we find that only one constant-roll inflation is consistent with the observations and observations constrain the constant-roll inflation to be slow-roll inflation. The tachyon potential is also reconstructed for the constant-roll inflation which is consistent with the observations.

The primordial power spectra of scalar and tensor perturbations during slow-roll inflation are usually calculated with the method of Bessel function approximation. For constant-roll or ultra slow-roll inflation, the method of Bessel function approximation may be invalid. We compare the numerical results with the analytical results derived from the Bessel function approximation, and we find that they differ significantly on super-horizon scales if the constant slow-roll parameter $\eta_H$ is not small. More accurate method is needed for calculating the primordial power spectrum for constant-roll inflation.

NGC 4945 is one of the nearest (~3.8 Mpc; 1" ~ 19 pc) starburst galaxies. ALMA band 3 (3--4\,mm) observations of HCN, HCO+, CS, C3H2, SiO, HCO, and CH3C2H were carried out with ~2" resolution. The lines reveal a rotating nuclear disk of projected size 10" x 2" with position angle ~45 deg, inclination ~75 deg and an unresolved bright central core of size <2.5". The continuum source (mostly free-free radiation) is more compact than the nuclear disk by a linear factor of two but shows the same position angle and is centered 0.39" +_ 0.14" northeast of the nuclear accretion disk defined by H2O maser emission. Outside the nuclear disk, both HCN and CS delineate molecular arms on opposite sides of the dynamical center. These are connected by a (deprojected) 0.6 kpc sized molecular bridge, likely a dense gaseous bar seen almost ends-on, shifting gas from the front and back side into the nuclear disk. Modeling this nuclear disk located farther inside <100 pc) with tilted rings indicates a coplanar outflow reaching a characteristic deprojectd velocity of ~50 km/s. All our molecular lines, with the notable exception of CH3C2H, show significant absorption near the systemic velocity (~571 km/s), within a range of ~500-660 km/s. Apparently, only molecular transitions with low critical H2-density do not show absorption. The velocity field of the nuclear disk, derived from CH3C2H, provides evidence for rigid rotation in the inner few arcseconds and a dynamical mass of M = (2.1+_0.2) x 10^8 Mo inside a galactocentric radius of 2.45", with a significantly flattened rotation curve farther out. Velocity integrated line intensity maps with most pronounced absorption show molecular peak positions up to 1.5" southwest of the continuum peak, presumably due to absorption, which appears to be most severe slightly northeast of the nuclear maser disk.

High energy cosmic ray electrons plus positrons (CREs), which lose energy quickly during their propagation, provide an ideal probe of Galactic high-energy processes and may enable the observation of phenomena such as dark-matter particle annihilation or decay. The CRE spectrum has been directly measured up to $\sim 2$ TeV in previous balloon- or space-borne experiments, and indirectly up to $\sim 5$ TeV by ground-based Cherenkov $\gamma$-ray telescope arrays. Evidence for a spectral break in the TeV energy range has been provided by indirect measurements of H.E.S.S., although the results were qualified by sizeable systematic uncertainties. Here we report a direct measurement of CREs in the energy range $25~{\rm GeV}-4.6~{\rm TeV}$ by the DArk Matter Particle Explorer (DAMPE) with unprecedentedly high energy resolution and low background. The majority of the spectrum can be properly fitted by a smoothly broken power-law model rather than a single power-law model. The direct detection of a spectral break at $E \sim0.9$ TeV confirms the evidence found by H.E.S.S., clarifies the behavior of the CRE spectrum at energies above 1 TeV and sheds light on the physical origin of the sub-TeV CREs.

We mapped the kinetic temperature structure of the Orion molecular cloud 1 with para-H2CO(303-202, 322-221, and 321-220) using the APEX 12m telescope. This is compared with the temperatures derived from the ratio of the NH3(2,2)/(1,1) inversion lines and the dust emission. Using the RADEX non-LTE model, we derive the gas kinetic temperature modeling the measured averaged line ratios of para-H2CO 322-221/303-202 and 321-220/303-202. The gas kinetic temperatures derived from the para-H2CO line ratios are warm, ranging from 30 to >200 K with an average of 62 K at a spatial density of 10$^5$ cm$^{-3}$. These temperatures are higher than those obtained from NH3(2,2)/(1,1) and CH3CCH(6-5) in the OMC-1 region. The gas kinetic temperatures derived from para-H2CO agree with those obtained from warm dust components measured in the mid infrared (MIR), which indicates that the para-H2CO(3-2) ratios trace dense and warm gas. The cold dust components measured in the far infrared (FIR) are consistent with those measured with NH3(2,2)/(1,1) and the CH3CCH(6-5) line series. With dust at MIR wavelengths and para-H2CO(3-2) on one side and dust at FIR wavelengths, NH3(2,2)/(1,1), and CH3CCH(6-5) on the other, dust and gas temperatures appear to be equivalent in the dense gas of the OMC-1 region, but provide a bimodal distribution, one more directly related to star formation than the other. The non-thermal velocity dispersions of para-H2CO are positively correlated with the gas kinetic temperatures in regions of strong non-thermal motion (Mach number >2.5) of the OMC-1, implying that the higher temperature traced by para-H2CO is related to turbulence on a 0.06 pc scale. Combining the temperature measurements with para-H2CO and NH3(2,2)/(1,1) line ratios, we find direct evidence for the dense gas along the northern part of the OMC-1 10 km s$^{-1}$ filament heated by radiation from the central Orion nebula.

The DArk Matter Particle Explorer (DAMPE), one of the four scientific space science missions within the framework of the Strategic Pioneer Program on Space Science of the Chinese Academy of Sciences, is a general purpose high energy cosmic-ray and gamma-ray observatory, which was successfully launched on December 17th, 2015 from the Jiuquan Satellite Launch Center. The DAMPE scientific objectives include the study of galactic cosmic rays up to $\sim 10$ TeV and hundreds of TeV for electrons/gammas and nuclei respectively, and the search for dark matter signatures in their spectra. In this paper we illustrate the layout of the DAMPE instrument, and discuss the results of beam tests and calibrations performed on ground. Finally we present the expected performance in space and give an overview of the mission key scientific goals.

Recent studies reveal that the free eccentricities of Kepler-34b and Kepler-413b are much larger than their forced eccentricities, implying that the scattering events may take place in their formation. The observed orbital configuration of Kepler-34b cannot be well reproduced in disk-driven migration models, whereas a two-planet scattering scenario can play a significant role of shaping the planetary configuration. These studies indicate that circumbinary planets discovered by Kepler may have experienced scattering process. In this work, we extensively investigate the scattering outcomes of circumbinary planets focusing on the effects of planet mass ratio. We find that the planetary mass ratio and the the initial relative locations of planets act as two important parameters which affect the eccentricity distribution of the surviving planets. As an application of our model, we discuss the observed orbital configurations of Kepler-34b and Kepler-413b. We first adopt the results from the disk-driven models as the initial conditions, then simulate the scattering process occurred in the late evolution stage of circumbinary planets. We show that the present orbital configurations of Kepler-34b and Kepler-413b can be well reproduced when considering two unequal-mass planet ejection model. Our work further suggests that some of the currently discovered circumbinary single-planet systems may be the survivals of original multiple-planet systems. The disk-driven migration and the scattering events occurring in the late stage both play an irreplaceable role in sculpting the final systems.

The detection of gravitational waves by the Laser Interferometer Gravitational-Wave Observatory opens a new era to use gravitational waves to test alternative theories of gravity. We investigate the polarizations of gravitational waves in $f(R)$ gravity and Horndeski theory, both containing scalar modes. These theories predict that in addition to the familiar $+$ and $\times$ polarizations, there are transverse breathing and longitudinal polarizations excited by the massive scalar mode and the new polarization is a single mixed state. It would be very difficult to detect the longitudinal polarization by interferometers, while pulsar timing array may be the better tool to detect the longitudinal polarization.

We derive a lower bound on the field excursion for the tachyon inflation, which is determined by the amplitude of the scalar perturbation and the number of $e$-folds before the end of inflation. Using the relation between the observables like $n_s$ and $r$ with the slow-roll parameters, we reconstruct three classes of tachyon potentials. The model parameters are determined from the observations before the potentials are reconstructed, and the observations prefer the concave potential. We also discuss the constraints from the reheating phase preceding the radiation domination for the three classes of models by assuming the equation of state parameter $w_{re}$ during reheating is a constant. Depending on the model parameters and the value of $w_{re}$, the constraints on $N_{re}$ and $T_{re}$ are different. As $n_s$ increases, the allowed reheating epoch becomes longer for $w_{re}=-1/3$, 0 and $1/6$ while the allowed reheating epoch becomes shorter for $w_{re}=2/3$.

HESS J1912+101 is a shell-like TeV source that has no clear counterpart in multiwavelength. Using CO and H i data, we reveal that VLSR~+60 km/s molecular clouds (MCs), together with shocked molecular gas and high-velocity neutral atomic shells, are concentrated toward HESS J1912+101. The prominent wing profiles up to VLSR~+80 km/s seen in 12CO (J=1-0 and J=3-2) data, as well as the high-velocity expanding H i shells up to VLSR~ +100 km/s, exhibit striking redshifted-broadening relative to the quiescent gas. These features provide compelling evidences for large-scale perturbation in the region. We argue that the shocked MCs and the high-velocity Hi shells may originate from an old supernova remnant (SNR). The distance to the SNR is estimated to be ~4.1 kpc based on the Hi self-absorption method, which leads to a physical radius of 29.0 pc for the ~(0.7-2.0)x10e5 years old remnant with an expansion velocity of >40 km/s. The +60 km/s MCs and the disturbed gas are indeed found to coincide with the bright TeV emission, supporting the physical association between them. Naturally, the shell-like TeV emission comes from the decay of neutral pions produced by interactions between the accelerated hadrons from the SNR and the surrounding high-density molecular gas.

We present new Effelsberg-100 m, ATCA, and VLA observations of rotational SiS transitions in the circumstellar envelope (CSE) of IRC +10216. Thanks to the high angular resolution achieved by the ATCA observations, we unambiguously confirm that the molecule's J=1-0 transition exhibits maser action in this CSE, as first suggested more than thirty years ago. The maser emission's radial velocity peaking at a local standard of rest velocity of -39.862$\pm$0.065 km/s indicates that it arises from an almost fully accelerated shell. Monitoring observations show time variability of the SiS (1-0) maser. The two lowest-$J$ SiS quasi-thermal emission lines trace a much more extended emitting region than previous high-J SiS observations. Their distributions show that the SiS quasi-thermal emission consists of two components: one is very compact (radius<1.5", corresponding to <3$\times 10^{15}$ cm), and the other extends out to a radius >11". An incomplete shell-like structure is found in the north-east, which is indicative of existing SiS shells. Clumpy structures are also revealed in this CSE. The gain of the SiS (1-0) maser (optical depths of about -5 at the blue-shifted side and, assuming inversion throughout the entire line's velocity range, about -2 at the red-shifted side) suggests that it is unsaturated. The SiS (1-0) maser can be explained in terms of ro-vibrational excitation caused by infrared pumping, and we propose that infrared continuum emission is the main pumping source.

We point out that there are only three polarizations for gravitational waves in $f(R)$ gravity, and the polarization due to the massive scalar mode is a mix of the pure longitudinal and transverse breathing polarization. The classification of the six polarizations by the Newman-Penrose quantities is based on weak, plane and null gravitational waves, so it is not applicable to the massive mode.

The recently measured 3.5 keV line in a number of galaxy clusters, the Andromeda galaxy (M31) and the Milky Way (MW) center can be well accounted for by a scenario in which dark matter decays to axion-like particles (ALPs) and subsequently convert to 3.5 keV photons in magnetic fields of galaxy clusters or galaxies. We propose to test this hypothesis by performing X-ray polarization measurements. Since ALPs can only couple to photons with polarization orientation parallel to magnetic field, we can confirm or reject this model by measuring the polarization of 3.5 keV line and comparing it to the orientation of magnetic field. We discuss luminosity and polarization measurements for both galaxy cluster and spiral galaxy, and provide a general relation between polarization and galaxy inclination angle. This effect is marginally detectable with X-ray polarimetry detectors currently under development, such as the enhanced X-ray Timing and Polarization (eXTP), the Imaging X-ray Polarimetry Explorer (IXPE) and the X-ray Imaging Polarimetry Explorer (XIPE). The sensitivity can be further improved in the future with detectors of larger effective area or better energy resolutions.

We present a new large-scale (4 square degrees) simultaneous $^{12}$CO, $^{13}$CO, and C$^{18}$O ($J$=1$-$0) mapping of L1188 with the PMO 13.7-m telescope. Our observations have revealed that L1188 consists of two nearly orthogonal filamentary molecular clouds at two clearly separated velocities. Toward the intersection showing large velocity spreads, we find several bridging features connecting the two clouds in velocity, and an open arc structure which exhibits high excitation temperatures, enhanced $^{12}$CO and $^{13}$CO emission, and broad $^{12}$CO line wings. This agrees with the scenario that the two clouds are colliding with each other. The distribution of young stellar object (YSO) candidates implies an enhancement of star formation in the intersection of the two clouds. We suggest that a cloud-cloud collision happened in L1188 about 1~Myr ago, possibly triggering the formation of low- and intermediate-mass YSOs in the intersection.

Intensity mapping is now becoming a useful tool to study the large-scale structure of the universe through spatial variations in the integrated emission from galaxies and the intergalactic medium. We study intensity mapping of the H-alpha 6563AA, [OIII]5007AA, [OII]3727AA and H-beta 4861AA lines at 0.8<z<5.2. The mean intensities of these four emission lines are estimated using the observed luminosity functions (LFs), cosmological simulations, and the star formation rate density (SFRD) derived from observations at z<5. We calculate the intensity power spectra and consider the foreground contamination of other lines at lower redshifts. We use the proposed NASA small explorer SPHEREx (the Spectro-Photometer for the History of the Universe, Epoch of Reionization, and Ices Explorer) as a case study for the detectability of the intensity power spectra of the four emission lines. We also investigate the cross correlation with the 21-cm line probed by CHIME (the Canadian Hydrogen Intensity Mapping Experiment), Tianlai experiment and SKA (the Square Kilometer Array) at 0.8<z<2.4. We find both the auto and cross power spectra can be well measured for the H-alpha, [OIII] and [OII] lines at z<3, while it is more challenging for the H-beta line. Finally, we estimate the constraint on the SFRD from intensity mapping, and find we can reach accuracy higher than 7% at z<4, which is better than usual measurements using the LFs of galaxies.

We show explicitly how the T model, E model, and Hilltop inflations are obtained from the general scalar-tensor theory of gravity with arbitrary conformal factors in the strong coupling limit. We argue that $\xi$ attractors can give any observables $n_s$ and $r$ by this method. The existence of attractors imposes a challenge to distinguish different models.

We derive the general formulae for the the scalar and tensor spectral tilts to the second order for the inflationary models with non-minimally derivative coupling without taking the high friction limit. The non-minimally kinetic coupling to Einstein tensor brings the energy scale in the inflationary models down to be sub-Planckian. In the high friction limit, the Lyth bound is modified with an extra suppression factor, so that the field excursion of the inflaton is sub-Planckian. The inflationary models with non-minimally derivative coupling are more consistent with observations in the high friction limit. In particular, with the help of the non-minimally derivative coupling, the quartic power law potential is consistent with the observational constraint at 95\% CL.

The nonminimal coupling of the kinetic term to Einstein's tensor helps the implementation of inflationary models due to the gravitationally enhanced friction. We calculate the parametrized post-Newtonian (PPN) parameters for the scalar--tensor theory of gravity with nonminimally derivative coupling. We find that under experimental constraint from the orbits of millisecond pulsars in our galaxy, the theory deviates from Einstein's general relativity in the order of $10^{-20}$, and the effect of the nonminimal coupling is negligible if we take the scalar field as dynamical dark energy. With the assumed conditions that the background scalar field is spatially homogeneous and evolves only on cosmological timescales and the contribution to stress--energy in the solar system from the background scalar field is subdominant, the scalar field is required to be massless.

We perform a phase space analysis of a generalized modified gravity theory with nonminimally coupling between geometry and matter. We apply the dynamical system approach to this generalized model and find that in the cosmological context, different choices of Lagrangian density will apparently result in different phases of the Universe. By carefully choosing the variables, we prove that there is an attractor solution to describe the late time accelerating universe when the modified gravity is chosen in a simple power-law form of the curvature scalar. We further examine the temperature evolution based on the thermodynamic understanding of the model. Confronting the model with supernova type Ia data sets, we find that the nonminimally coupled theory of gravity is a viable model to describe the late time Universe acceleration.

The observational data on the anisotropy of the cosmic microwave background constraints the scalar spectral tilt $n_s$ and the tensor to scalar ratio $r$ which depend on the first and second derivatives of the inflaton potential. The information can be used to reconstruct the inflaton potential in the polynomial form up to some orders. However, for some classes of potentials, $n_s$ and $r$ behave as $n_s(N)$ and $r(N)$ universally in terms of the number of e-folds $N$. The universal behaviour of $n_s(N)$ can be used to reconstruct a class of inflaton potentials. By parametrizing one of the parameters $n_s(N)$, $\epsilon(N)$ and $\phi(N)$, and fitting the parameters in the models to the observational data, we obtain the constraints on the parameters and reconstruct the classes of the inflationary models which include the chaotic inflation, T-model, hilltop inflation, s-dual inflation, natural inflation and $R^2$ inflation.

The extragalactic background light (EBL) is comprised of the cumulative radiation from all galaxies and active galactic nuclei over the cosmic history. In addition to point sources, EBL also contains information from diffuse sources of radiation. The angular power spectra of the near-infrared intensities could contain additional signals and a complete understanding of the nature of the IR background is still lacking in the literature. Here we explore the constraints that can be placed on particle decays, especially candidate dark matter models involving axions that trace dark matter halos of galaxies. Axions with a mass around a few eV will decay via two photons with wavelengths in the near-IR band, and will leave a signature in the IR background intensity power spectrum. Using recent power spectra measurements from the Hubble Space Telescope (HST) and Cosmic Infrared Background Experiment (CIBER), we find that the 0.6 to 1.6 micron power spectra can be explained by axions with masses around 4 eV. The total axion abundance Omega_a~0.05, and it is comparable to the baryon density of the Universe. The suggested mean axion mass and abundance are not ruled out by existing cosmological observations. Interestingly, the axion model with a mass distribution is preferred by the data, which cannot be explained by the standard quantum chromodynamics (QCD) theory and needs further discussion.

Cosmic Dawn Intensity Mapper is a "Probe Class" mission concept for reionization studies of the universe. It will be capable of spectroscopic imaging observations between 0.7 to 6-7 microns in the near-Infrared. The primary observational objective is pioneering observations of spectral emission lines of interest throughout the cosmic history, but especially from the first generation of distant, faint galaxies when the universe was less than 800 million years old. With spectro-imaging capabilities, using a set of linear variable filters (LVFs), CDIM will produce a three-dimensional tomographic view of the epoch of reionization (EoR). CDIM will also study galaxy formation over more than 90% of the cosmic history and will move the astronomical community from broad-band astronomical imaging to low-resolution (R=200-300) spectro-imaging of the universe.

With the usual definitions for the entropy and the temperature associated with the apparent horizon, we discuss the first law of the thermodynamics on the apparent in the general scalar-tensor theory of gravity with the kinetic term of the scalar field non-minimally coupling to Einstein tensor. We show the equivalence between the first law of thermodynamics on the apparent horizon and Friedmann equation for the general models, by using a mass-like function which is equal to the Misner-Sharp mass on the apparent horizon. The results further support the universal relationship between the first law of thermodynamics and Friedmann equation.

As a part of the Milky Way Imaging Scroll Painting (MWISP) survey, we performed a simultaneous 12CO(1-0), 13CO(1-0), C18O(1-0) mapping toward molecular clouds in a region encompassing 3.75 square degrees. We reveal three molecular clouds, the 15 km/s cloud, the 27 km/s cloud, and the 50 km/s cloud, in the surveyed region. The 50 km/s cloud is resolved with an angular resolution of around 1' for the first time. Investigating their morphology and velocity structures, we find that the 27 km/s cloud is likely affected by feedback from the stellar association Mon OB3 and the 50 km/s cloud is characterised by three large expanding molecular shells. The region is mapped in C18O (1-0) for the first time. We discover seven C18O clumps, which are likely to form massive stars, and 15 dust clumps based on the BGPS archival data. Using infrared color-color diagrams, we find 56 Class I and 107 Class II young stellar object (YSO) candidates. According to the distribution of YSO candidates, an overdensity is found around the HII region S287 and the intersection of two shells, indicative of triggering. The star formation efficiency and rate of the 27~km/s cloud are discussed. Comparing the observed values of the filament S287-main with the models of fragmentation, we suggest that turbulence controls the large scale of fragmentation in the filament while gravitational fragmentation plays an important role in the formation of YSOs on the small scale. We find that star-forming gas tends to have higher excitation temperature, higher 13CO opacity, and higher column density than non-star-forming gas, which is consistent with the point that star formation happens in denser gas and star-forming gas is heated by YSOs. Using the 1.1 mm dust emission to trace dense gas, we obtain a dense gas fraction of 2.7%-10.4% for the 27 km/s cloud.

OB type stars have strong ionizing radiation, and drive energetic winds. The ultraviolet (UV) radiation from ionizing stars may heat dust and ionize gas to sweep up an expanding bubble shell. This shell may be the result of feedback leading to a new generation of stars. N131 is an infrared dust bubble residing in a molecular filament. We study the formation and fragmentation of this bubble with multi-wavelength dust and gas observations. Towards the bubble N131, we analyzed archival multi-wavelength observations including 3.6, 4.5, 5.8, 8.0, 24, 70, 160, 250, 350, 500 $\mu$m, 1.1 mm, and 21 cm. In addition, we performed new observations of CO (2-1), CO (1-0), and $^{13}$CO (1-0) with the IRAM 30-m telescope. Multi-wavelength dust and gas observations reveal a ringlike shell with compact fragments, two filamentary structures, and a secondary bubble N131-A. The bubble N131 is a rare object with a large hole at 24 $\mu$m and 21 cm in the direction of its center. The dust and gas clumps are compact and might have been compressed at the inner edge of the ringlike shell, while they are extended and might be pre-existing at the outer edge. The column density, excitation temperature, and velocity show a potentially hierarchical distribution from the inner to outer edge of the ringlike shell. We also detected the front and back sides of the secondary bubble N131-A in the direction of its center. The derived Lyman-continuum ionizing photon flux within N131-A is equivalent to an O9.5 star. Based on the above, we suggest that the bubble N131 might be triggered by the strong stellar winds from a group of massive stars inside the bubble. We propose a scenario in which the bubble N131 forms from the disruption of a gas filament by expansion of HII region, strong stellar winds, and fragments under self-gravity.

TianQin is a proposal for a space-borne detector of gravitational waves in the millihertz frequencies. The experiment relies on a constellation of three drag-free spacecraft orbiting the Earth. Inter-spacecraft laser interferometry is used to monitor the distances between the test masses. The experiment is designed to be capable of detecting a signal with high confidence from a single source of gravitational waves within a few months of observing time. We describe the preliminary mission concept for TianQin, including the candidate source and experimental designs. We present estimates for the major constituents of the experiment's error budget and discuss the project's overall feasibility. Given the current level of technology readiness, we expect TianQin to be flown in the second half of the next decade.

We combine new Cosmic Microwave Background (CMB) data from Planck with Baryon Acoustic Oscillation (BAO) data to constrain the Brans-Dicke (BD) theory, in which the gravitational constant $G$ evolves with time. Observations of type Ia supernovae (SNeIa) provide another important set of cosmological data, as they may be regarded as standard candles after some empirical corrections. However, in theories that include modified gravity like the BD theory, there is some risk and complication when using the SNIa data because their luminosity may depend on $G$. In this paper, we assume a power law relation between the SNIa luminosity and $G$, but treat the power index as a free parameter. We then test whether the difference in distances measured with SNIa data and BAO data can be reduced in such a model. We also constrain the BD theory and cosmological parameters by making a global fit with the CMB, BAO and SNIa data set. For the CMB+BAO+SNIa data set, we find $0.08\times10^{-2} < \zeta <0.33\times10^{-2} $ at the 68\% confidence level (CL) and $-0.01\times10^{-2} <\zeta <0.43\times 10^{-2} $ at the 95\% CL, where $\zeta$ is related to the {BD} parameter $\omega$ by $\zeta=\ln(1+1/\omega)$.

The spatial fluctuations of the extragalactic background light trace the total emission from all stars and galaxies in the Universe. A multi-wavelength study can be used to measure the integrated emission from first galaxies during reionization when the Universe was about 500 million years old. Here we report arcminute-scale spatial fluctuations in one of the deepest sky surveys with the Hubble Space Telescope in five wavebands between 0.6 and 1.6 $\mu$m. We model-fit the angular power spectra of intensity fluctuation measurements to find the ultraviolet luminosity density of galaxies at $z$ > 8 to be $\log \rho_{\rm UV} = 27.4^{+0.2}_{-1.2}$ erg s$^{-1}$ Hz$^{-1}$ Mpc$^{-3}$ $(1\sigma)$. This level of integrated light emission allows for a significant surface density of fainter primeval galaxies that are below the point source detection level in current surveys.

We propose a new and robust method to test the consistency of the cosmic evolution given by a cosmological model. It is realized by comparing the combined quantity r_d^CMB/D_V^SN, which is derived from the comoving sound horizon r_d from cosmic microwave background (CMB) measurements and the effective distance D_V derived from low-redshift Type-Ia supernovae (SNe Ia) data, with direct and independent r_d/D_V obtained by baryon acoustic oscillation (BAO) measurements at median redshifts. We apply this test method for the Lambda-CDM and wCDM models, and investigate the consistency of the derived value of r_d/D_V from Planck 2015 and the SN Ia data sets of Union2.1 and JLA (z<1.5), and the r_d/D_V directly given by BAO data from six-degree-field galaxy survey (6dFGS), Sloan Digital Sky Survey Data Release 7 Main Galaxy Survey (SDSS-DR7 MGS), DR11 of SDSS-III, WiggleZ and Ly-alpha forecast surveys from Baryon Oscillation Spectroscopic Data (BOSS) DR-11 over 0.1<z<2.36. We find that r_d^CMB/D_V^SN for both non-flat Lambda-CDM and flat wCDM models with Union2.1 and JLA data are well consistent with the BAO and CMB measurements within 1-sigma CL. Future surveys will further tight up the constraints significantly, and provide stronger test on the consistency.

Orion KL has served as a benchmark for spectral line searches throughout the (sub)millimeter regime. The main goal is to systematically study spectral characteristics of Orion KL in the 1.3 cm band. We carried out a spectral line survey (17.9 GHz to 26.2 GHz) with the Effelsberg-100 m telescope towards Orion KL. We find 261 spectral lines, yielding an average line density of about 32 spectral features per GHz above 3$\sigma$. The identified lines include 164 radio recombination lines (RRLs) and 97 molecular lines. A total of 23 molecular transitions from species known to exist in Orion KL are detected for the first time in the interstellar medium. Non-metastable 15NH3 transitions are detected in Orion KL for the first time. Based on the velocity information of detected lines and the ALMA images, the spatial origins of molecular emission are constrained and discussed. A narrow feature is found in SO2 ($8_{1,7}-7_{2,6}$), possibly suggesting the presence of a maser line. Column densities and fractional abundances relative to H2 are estimated for 12 molecules with LTE methods. Rotational diagrams of non-metastable 14NH3 transitions with J=K+1 to J=K+4 yield different results; metastable 15NH3 is found to have a higher excitation temperature than non-metastable 15NH3, indicating that they may trace different regions. Elemental and isotopic abundance ratios are estimated: 12C/13C=63+-17, 14N/15N=100+-51, D/H=0.0083+-0.0045. The dispersion of the He/H ratios derived from H$\alpha$/He$\alpha$ pairs to H$\delta$/He$\delta$ pairs is very small, which is consistent with theoretical predictions that the departure coefficients bn factors for hydrogen and helium are nearly identical. Based on a non-LTE code neglecting excitation by the infrared radiation field and a likelihood analysis, we find that the denser regions have lower kinetic temperature, which favors an external heating of the Hot Core.

In quantum theory of gravity, we expect the Lorentz Invariance Violation (LIV) and the modification of the dispersion relation between energy and momentum for photons. The effect of the energy-dependent velocity due to the modified dispersion relation for photons was studied in the standard cosmological context by using a sample of Gamma Ray Bursts (GRBs). In this paper we mainly discuss the possible LIV effect by using different cosmological models for the accelerating universe. Due to the degeneracies among model parameters, the GRBs' time delay data are combined with the cosmic microwave background data from the Planck first year release, the baryon acoustic oscillation data at six different redshifts, as well as Union2 type Ia supernovae data, to constrain both the model parameters and the LIV effect. We find no evidence of LIV.

We produce simulations of emission of the atomic CII line in large sky fields in order to determine the current prospects for mapping this line during the high redshift Epoch of Reionization. We estimate the CII line intensity, redshift evolution and spatial fluctuations using observational relations between CII emission and the SFR in a galaxy for the frequency range of 200 GHz to 300 GHz. We obtained a frequency averaged intensity of CII emission of ${\rm I_{\rm CII}=(4 \pm 2)\times10^{2}\, Jy\, \rm sr^{-1}}$ in the redshift range $z\, \sim\, 5.3\, -\, 8.5$. Observations of CII emission in this frequency range will suffer contamination from emission lines at lower redshifts, in particular from the CO rotation lines. For the relevant frequency range we estimated the CO contamination (originated in emission from galaxies at $z\, <\, 2.5$), using simulations, to be ${\rm I_{\rm CO} \approx 10^{3}\, Jy \, sr^{-1}}$ and independently confirmed the result based in observational relations. We generated maps as a function of angle and frequency using detailed simulations of the CII and CO emission across several redshifts in order to properly take into account the observational pipeline and light cone effects. In order to reduce the foreground contamination we found that we should mask galaxies below redshifts $\sim 2.5$ with a CO flux in one of the CO(J:2-1) to CO(J:6-5) lines higher than ${\rm 5\times 10^{-22}\, W\ m^{-2}}$ or a AB magnitude lower than ${\rm m_{\rm K}\, =\, 22}$. We estimate that the additional continuum contamination is of the order of ${\rm 10^{5}\, Jy\, sr^{-1}}$. It is also considered the possibility of cross correlating foreground lines with galaxies in order to probe the intensity of the foregrounds.

SPHEREx (Spectro-Photometer for the History of the Universe, Epoch of Reionization, and Ices Explorer) ( http://spherex.caltech.edu ) is a proposed all-sky spectroscopic survey satellite designed to address all three science goals in NASA's Astrophysics Division: probe the origin and destiny of our Universe; explore whether planets around other stars could harbor life; and explore the origin and evolution of galaxies. SPHEREx will scan a series of Linear Variable Filters systematically across the entire sky. The SPHEREx data set will contain R=40 spectra fir 0.75$<\lambda<$4.1$\mu$m and R=150 spectra for 4.1$<\lambda<$4.8$\mu$m for every 6.2 arc second pixel over the entire-sky. In this paper, we detail the extra-galactic and cosmological studies SPHEREx will enable and present detailed systematic effect evaluations. We also outline the Ice and Galaxy Evolution Investigations.

To reconcile the BICEP2 measurement on the tensor-to-scalar ratio $r$ with Planck constraint, a large negative running of scalar spectral index $n_s$ is needed. So the inflationary observable such as $n_s$ should be expanded at least to the second-order slow-roll parameters for single-field inflationary models. The large value of $r$ and the Lyth bound indicate that it is impossible to obtain the sub-Planckian excursion for the inflaton. However, we derive an absolutely minimal bound $\Delta\phi/M_{\rm Pl}>\sqrt{r/2}$ on the inflaton excursion for single-field inflationary models, which can be applied to non-slow-roll inflationary models as well. This bound excludes the possibility of the small-field inflation with $\Delta\phi<0.1 M_{\rm Pl}$ if the BICEP2 result on $r$ stands and it opens the window of sub-Planckian excursion with $\Delta\phi<M_{\rm Pl}$ even if $r$ is as large as $0.1$. To get the sub-Planckian excursion with $\Delta\phi<0.1 M_{\rm Pl}$, our modified bound requires $r<0.02$. Using a fifth-order polynomial potential as an explicit example, we show that it not only agrees with the observational results, but also violates the Lyth bound.

We present the science enabled by cross-correlations of the SKA1-LOW 21-cm EoR surveys with other line mapping programs. In particular, we identify and investigate potential synergies with planned programs, such as the line intensity mapping of redshifted CO rotational lines, [CII] and Ly-$\alpha$ emissions during reionization. We briefly describe how these tracers of the star-formation rate at $z \sim 8$ can be modeled jointly before forecasting their auto- and cross-power spectra measurements with the nominal 21cm EoR survey. The use of multiple line tracers would be invaluable to validate and enrich our understanding of the EoR.

We present the cross-correlation between the far-infrared background fluctuations as measured with the Herschel Space Observatory at 250, 350, and 500 {\mu}m and the near-infrared background fluctuations with Spitzer Space Telescope at 3.6 {\mu}m. The cross-correlation between far and near-IR background anisotropies are detected such that the correlation coefficient at a few to ten arcminute angular scales decreases from 0.3 to 0.1 when the far-IR wavelength increases from 250 {\mu}m to 500 {\mu}m. We model the cross-correlation using a halo model with three components: (a) far-IR bright or dusty star-forming galaxies below the masking depth in Herschel maps, (b) near-IR faint galaxies below the masking depth at 3.6 {\mu}m, and (c) intra-halo light, or diffuse stars in dark matter halos, that likely dominates fluctuations at 3.6 {\mu}m. The model is able to reasonably reproduce the auto correlations at each of the far-IR wavelengths and at 3.6 {\mu}m and their corresponding cross-correlations. While the far and near-IR auto-correlations are dominated by faint dusty, star-forming galaxies and intra-halo light, respectively, we find that roughly half of the cross-correlation between near and far-IR backgrounds is due to the same galaxies that remain unmasked at 3.6 {\mu}m. The remaining signal in the cross-correlation is due to intra-halo light present in the same dark matter halos as those hosting the same faint and unmasked galaxies. In this model, the decrease in the cross-correlation signal from 250 {\mu}m to 500 {\mu}m comes from the fact that the galaxies that are primarily contributing to 500 {\mu}m fluctuations peak at a higher redshift than those at 250 {\mu}m.

IRC +10216 is the prototypical carbon star exhibiting an extended molecular circumstellar envelope. Its spectral properties are therefore the template for an entire class of objects. The main goal is to systematically study the $\lambda$ $\sim$1.3 cm spectral line characteristics of IRC +10216. We carried out a spectral line survey with the Effelsberg-100 m telescope toward IRC +10216. It covers the frequency range between 17.8 GHz and 26.3 GHz (K-band). In the circumstellar shell of IRC +10216, we find 78 spectral lines, among which 12 remain unidentified. The identified lines are assigned to 18 different molecules and radicals. A total of 23 lines from species known to exist in this envelope are detected for the first time outside the Solar System and there are additional 20 lines first detected in IRC +10216. The potential orgin of "U" lines is also discussed. Assuming local thermodynamic equilibrium (LTE), we then determine rotational temperatures and column densities of 17 detected molecules. Molecular abundances relative to H$_{2}$ are also estimated. A non-LTE analysis of NH$_{3}$ shows that the bulk of its emission arises from the inner envelope with a kinetic temperature of 70$\pm$20 K. Evidence for NH$_{3}$ emitting gas with higher kinetic temperature is also obtained, and potential abundance differences between various $^{13}$C-bearing isotopologues of HC$_{5}$N are evaluated. Overall, the isotopic $^{12}$C/$^{13}$C ratio is estimated to be 49$\pm$9. Finally, a comparison of detected molecules in the $\lambda$ $\sim$1.3 cm range with the dark cloud TMC-1 indicates that silicate-bearing molecules are more predominant in IRC +10216.

Extragalactic background light (EBL) anisotropy traces variations in the total production of photons over cosmic history, and may contain faint, extended components missed in galaxy point source surveys. Infrared EBL fluctuations have been attributed to primordial galaxies and black holes at the epoch of reionization (EOR), or alternately, intra-halo light (IHL) from stars tidally stripped from their parent galaxies at low redshift. We report new EBL anisotropy measurements from a specialized sounding rocket experiment at 1.1 and 1.6 micrometers. The observed fluctuations exceed the amplitude from known galaxy populations, are inconsistent with EOR galaxies and black holes, and are largely explained by IHL emission. The measured fluctuations are associated with an EBL intensity that is comparable to the background from known galaxies measured through number counts, and therefore a substantial contribution to the energy contained in photons in the cosmos.

The BICEP2 experiment confirms the existence of primordial gravitational wave with the tensor-to-scalar ratio $r=0$ ruled out at $7\sigma$ level. The consistency of this large value of $r$ with the {\em Planck} data requires a large negative running $n'_s$ of the scalar spectral index. Herein we propose two types of the single field inflation models with simple potentials to study the possibility of the consistency of the models with the BICEP2 and {\em Planck} observations. One type of model suggested herein is realized in the supergravity model building. These models fail to provide the needed $n'_s$ even though both can fit the tensor-to-scalar ratio and spectral index.

The detection of B-mode power spectrum by the BICEP2 collaboration constrains the tensor-to-scalar ratio $r=0.20^{+0.07}_{-0.05}$ for the lensed-$\Lambda$CDM model. The consistency of this big value with the {\em Planck} results requires a large running of the spectral index. The large values of the tensor-to-scalar ratio and the running of the spectral index put a challenge to single field inflation. For the chaotic inflation, the larger the value of the tensor-to-scalar ratio is, the smaller the value of the running of the spectral index is. For the natural inflation, the absolute value of the running of the spectral index has an upper limit.

Planck measurements of the cosmic microwave background power spectra find a lower value of the Hubble constant $H_0$ and a higher value of the fractional matter energy density $\Omega_{m0}$ for the concordance $\Lambda$CDM model, and these results are in tension with other measurements. The {\em Planck} group argued that the tension came either from some sources of unknown systematic errors in some astrophysical measurements or the wrong $\Lambda$CDM model applied in fitting the data. We studied the reason for the tension on $H_0$ from different measurements by considering two dynamical dark energy models. We found that there is no tension between different data, the constraint on $H_0$ is almost unchanged for different dark energy models and the tension with the local measurements remains when the error bar on $H_0$ is tightened to be around 1. We argue that the tension on $H_0$ is not caused by the fitting model.

We discuss the general dynamical behaviors of quintessence field, in particular, the general conditions for tracking and thawing solutions are discussed. We explain what the tracking solutions mean and in what sense the results depend on the initial conditions. Based on the definition of tracking solution, we give a simple explanation on the existence of a general relation between $w_\phi$ and $\Omega_\phi$ which is independent of the initial conditions for the tracking solution. A more general tracker theorem which requires large initial values of the roll parameter is then proposed. To get thawing solutions, the initial value of the roll parameter needs to be small. The power-law and pseudo-Nambu Goldstone boson potentials are used to discuss the tracking and thawing solutions. A more general $w_\phi-\Omega_\phi$ relation is derived for the thawing solutions. Based on the asymptotical behavior of the $w_\phi-\Omega_\phi$ relation, the flow parameter is used to give an upper limit on $w_\phi'$ for the thawing solutions. If we use the observational constraint $w_{\phi 0}<-0.8$ and $0.2<\Omega_{m0}<0.4$, then we require $n\lesssim 1$ for the inverse power-law potential $V(\phi)=V_0(\phi/m_{pl})^{-n}$ with tracking solutions and the initial value of the roll parameter $|\lambda_i|<1.3$ for the potentials with the thawing solutions.

The intensity mapping of Ly-alpha emission during the epoch of reionization (EoR) will be contaminated by foreground emission lines from lower redshifts. We calculate the mean intensity and power spectrum of Ly-alpha emission at z~7, and estimate the uncertainties according to the relevant astrophysical processes. We find that the low-redshift emission lines from 6563 A H-alpha, 5007 A [OIII] and 3727 A [OII] will be the strong contaminants on the observed Ly-alpha power spectrum. We make use of both the star formation rate (SFR) and luminosity functions (LF) to estimate the mean intensity and power spectra of the three foreground lines at z~0.5 for H-alpha, z~0.9 for [OIII] and z~1.6 for [OII], as they will contaminate the Ly-alpha emission at z~7. The [OII] line is found to be the strongest. We analyze the masking of the bright survey pixels with a foreground line above some line intensity threshold as a way to reduce the contamination in the intensity mapping survey. We find that the foreground contamination can be neglected if we remove the pixels with fluxes above 1.4x10^-20 W/m^2.

The dynamics of scalar fields as dark energy is well approximated by some general relations between the equation of state parameter $w(z)$ and the fraction energy density $\Omega_\phi$. Based on the approximation, for slowly-rolling scalar fields, we derived the analytical expressions of $w(z)$ which reduce to the popular Chevallier-Polarski-Linder parametrization with explicit degeneracy relation between $w_0$ and $w_a$. The models approximate the dynamics of scalar fields well and help eliminate the degeneracies among $w_a$, $w_0$ and $\Omega_{m0}$. With the explicit degeneracy relations, we test their effects on the constraints of cosmological parameters. We find that: (1) The analytical relations between $w_0$ and $w_a$ for the two models are consistent with observational data; (2) The degeneracies have little effect on $\Omega_{m0}$; (3) The $1\sigma$ error of $w_0$ was reduced about 30% with the degeneracy relations.

The transient acceleration which the current cosmic acceleration is not eternal is possible by introducing the interaction between dark matter and dark energy. If the energy transfer is from dark energy to dark matter, then it is possible to realize the transient acceleration. We study the possibility of transient acceleration by considering two time-dependent phenomenological interaction forms so that the energy transfer increases as the universe evolves. Starting from a simple and extending to a more complicated ansatz, we obtain analytical expressions for the evolutions of the deceleration and the various energy density parameters. We find the ranges of the parameters in the models for a transient acceleration.

In this letter, the effects of dark radiation (DR) are tested. Theoretically, the phase-space analysis method is applied to check whether the model is consist with the history of our universe which shows positive results. Observationally, by using the observational data ($SNLS$ (SuperNovae Legacy Survey), $WMAP9$(Wilkinson Microwave Anisotropy Probe 9 Years Result), $PLANCK$ (Planck First Data Release), $BAO$ (Baryon Acoustic Oscillations), $H(z)$ (Hubble Parameter Data) and $BBN$ (Big Bang Nucleosynthesis)), the dark radiation is found to have the effect of wiping out the tension between the $SNLS$ data and the other data in flat $\Lambda CDM$ model. The effects of dark radiation also make the best fit value of $N_{eff}$ slightly larger than 3.04.

The attenuation of high-energy gamma-ray spectrum due to the electron-positron pair production against the extragalactic background light (EBL) provides an indirect method to measure the EBL of the universe. We use the measurements of the absorption features of the gamma-rays from blazars as seen by Fermi Gamma-ray Space Telescope to explore the EBL flux density and constrain the EBL spectrum, star formation rate density (SFRD) and photon escape fraction from galaxies out to z=6. Our results are basically consistent with the existing determinations of the quantities. We find a larger photon escape fraction at high redshifts, especially at z=3, compared to the result from the recent Ly-alpha measurements. Our SFRD result is consistent with the data from both gamma-ray burst and UV observations in 1-sigma level. However, the average SFRD we obtain at z>~3 matches the gamma-ray data better than the UV data. Thus our SFRD result at z>~6 favors that it is sufficiently high enough to reionize the universe.

Compared to the giant planets in the solar system, exoplanets have many remarkable properties such as the prevalence of giant planets on eccentric orbits and the presence of hot Jupiters. Planet-planet scattering (PPS) between giant planets is a possible mechanism in interpreting above and other observed properties. If the observed giant planet architectures are indeed the outcomes of PPS, such drastic dynamical process must affect their primordial moon systems. In this Letter, we discuss the effect of the PPS on the survival of their regular moons. From the viewpoint of observations, some preliminary conclusions are drawn from the simulations. 1. PPS is a destructive process to the moon systems, single planets on eccentric orbits are not the ideal moon-search targets. 2. If hot Jupiters formed through PPS, their original moons have little chance to survive. 3. Planets in multiple systems with small eccentricities are more likely holding their primordial moons. 4. Compared to the lower-mass planets, the massive ones in multiple systems may not be the preferred moon-search targets if the system underwent a PPS history.

We consider a dark energy model with a relation between the equation of state parameter $w$ and the energy density parameter $\Omega_\phi$ derived from thawing scalar field models. Assuming the variation of the fine structure constant is caused by dark energy, we use the observational data of the variation of the fine structure constant to constrain the current value of $w_0$ and $\Omega_{\phi 0}$ for the dark energy model. At the $1\sigma$ level, the observational data excluded some areas around $w_0=-1$, which explains the positive detection of the variation of the fine structure constant at the $1\sigma$ level, but $\Lambda$CDM model is consistent with the data at the $2\sigma$ level.

The constraints on the $\Lambda$CDM model from type Ia supernova (SNe Ia) data alone and BAO data alone are similar, so it is worthwhile to compare their constraints on the property of dark energy. We apply the SNLS3 compilation of 472 SNe Ia data, the baryon acoustic oscillation measurement of distance, the cosmic microwave background radiation data from the WMAP7, and the Hubble parameter data to study the effect of their different combinations on the fittings of cosmological parameters in the modified holographic dark energy model and the Chevallier-Polarski-Linder model. Neither BAO nor WMAP7 data alone give good constraint on the equation of state parameter of dark energy, but both WMAP7 data and BAO data help SNe Ia data break the degeneracies among the model parameters, hence tighten the constraint on the variation of equation of state parameter $w_a$, and WMAP7 data do the job a little better. Although BAO and WMAP7 data provide reasonably good constraints on $\Omega_m$ and $\Omega_k$, they are not able to constrain the dynamics of dark energy. On the other hand, SNe Ia data do not provide good constraints on $\Omega_m$ and $\Omega_k$, but they provide good constraint on the dynamics of dark energy, especially the variation of the equation of state parameter of dark energy, so we need to combine SNe Ia with BAO and WMAP7 data to probe the property of dark energy. The addition of $H(z)$ data helps better constrain the geometry of the Universe $\Omega_k$ and the property of dark energy. For the SNLS SNe Ia data, the nuisance parameters $\alpha$ and $\beta$ are consistent for all different combinations of the above data, and their impacts on the fittings of cosmological parameters are minimal. By fitting the data to different models, $\Lambda$CDM model is still consistent with all the observational data.

We calculate the absolute intensity and anisotropies of the Lyman-alpha radiation field present during the epoch of reionization. We consider emission from both galaxies and the intergalactic medium (IGM) and take into account the main contributions to the production of Lyman-alpha photons: recombinations, collisions, continuum emission from galaxies and scattering of Lyman-n photons in the IGM. We find that the emission from individual galaxies dominates over the IGM with a total Lyman-alpha intensity (times frequency) of about (1.43-3.57)x10^{-8} erg s^{-1} cm^{-2} sr^{-1} at a redshift of 7. This intensity level is low so it is unlikely that the Lyman-\alpha background during reionization can be established by an experiment aiming at an absolute background light measurement. Instead we consider Lyman-\alpha intensity mapping with the aim of measuring the anisotropy power spectrum which has rms fluctuations at the level of 1 x 10^{-16} [erg s^[-1} cm^{-2} sr^{-1}]^2 at a few Mpc scales. These anisotropies could be measured with a spectrometer at near-IR wavelengths from 0.9 to 1.4 \mu m with fields in the order of 0.5 to 1 sq. degrees. We recommend that existing ground-based programs using narrow band filters also pursue intensity fluctuations to study statistics on the spatial distribution of faint Lyman-\alpha emitters. We also discuss the cross-correlation signal with 21 cm experiments that probe HI in the IGM during reionization. A dedicated sub-orbital or space-based Lyman-\alpha intensity mapping experiment could provide a viable complimentary approach to probe reionization, when compared to 21 cm experiments, and is likely within experimental reach.

Molecular hydrogen is now understood to be the main coolant of the primordial gas clouds leading to the formation of the very first stars and galaxies. The line emissions associated with molecular hydrogen should then be a good tracer of the matter distribution at the onset of reionization of the universe. Here we propose intensity mapping of H2 line emission in rest-frame mid-infrared wavelengths to map out the spatial distribution of gas at redshifts z > 10. We calculate the expected mean intensity and clustering power spectrum for several H2 lines. We find that the 0-0S(3) rotational line at a rest wavelength of 9.66 microns is the brightest line over the redshift range of 10 to 30 with an intensity of about 5 to 10 Jy/sr at z~15. To reduce astrophysical and instrumental systematics, we propose the cross-correlation between multiple lines of the H2 rotational and vibrational line emission spectrum. Our estimates of the intensity can be used as a guidance in planning instruments for future mid-IR spectroscopy missions such as SPICA.

According to the core accretion theory, circumbinary embryos can form only beyond a critical semimajor axis (CSMA). However, due to the relatively high density of solid materials in the inner disk, significant amount of small planetesimals must exist in the inner zone when embryos were forming outside this CSMA. So embryos migration induced by the planetesimal swarm is possible after the gas disk depletion. Through numerical simulations, we found (i) the scattering-driven inward migration of embryos is robust, planets can form in the habitable zone if we adopt a mass distribution of MMSN-like disk; (ii) the total mass of the planetesimals in the inner region and continuous embryo-embryo scattering are two key factors that cause significant embryo migrations; (iii) the scattering-driven migration of embryos is a natural water-deliver mechanism. We propose that planet detections should focus on the close binary with its habitable zone near CSMA.

We use the new gamma-ray bursts (GRBs) data, combined with the baryon acoustic oscillation(BAO) observation from the spectroscopic Sloan Digital Sky Survey (SDSS) data release, the newly obtained $A$ parameter at $z=0.6$ from the WiggleZ Dark Energy Survey, the cosmic microwave background (CMB) observations from the 7-Year Wilkinson Microwave Anisotropy Probe (WMAP7) results, and the type Ia supernovae (SNeIa) from Union2 set, to constrain a phenomenological model describing possible interactions between dark energy and dark matter, which was proposed to alleviate the coincidence problem of the standard $\Lambda$CDM model. By using the Markov Chain Monte Carlo (MCMC) method, we obtain the marginalized $1\sigma$ constraints $\Omega_{m}=0.2886\pm{0.0135}$, $r_m=-0.0047\pm{0.0046}$, and $w_X=-1.0658\pm{0.0564}$. We also consider other combinations of these data for comparison. These results show that: (1) the energy of dark matter is slightly transferring to that of dark energy; (2) even though the GRBs+BAO+CMB data present less stringent constraints than SNe+BAO+CMB data do, the GRBs can help eliminate the degeneracies among parameters.

Using several cosmological observations, i.e. the cosmic microwave background anisotropies (WMAP), the weak gravitational lensing (CFHTLS), the measurements of baryon acoustic oscillations (SDSS+WiggleZ), the most recent observational Hubble parameter data, the Union2.1 compilation of type Ia supernovae, and the HST prior, we impose constraints on the sum of neutrino masses ($\mnu$), the effective number of neutrino species ($\neff$) and dark energy equation of state ($w$), individually and collectively. We find that a tight upper limit on $\mnu$ can be extracted from the full data combination, if $\neff$ and $w$ are fixed. However this upper bound is severely weakened if $\neff$ and $w$ are allowed to vary. This result naturally raises questions on the robustness of previous strict upper bounds on $\mnu$, ever reported in the literature. The best-fit values from our most generalized constraint read $\mnu=0.556^{+0.231}_{-0.288}\rm eV$, $\neff=3.839\pm0.452$, and $w=-1.058\pm0.088$ at 68% confidence level, which shows a firm lower limit on total neutrino mass, favors an extra light degree of freedom, and supports the cosmological constant model. The current weak lensing data are already helpful in constraining cosmological model parameters for fixed $w$. The dataset of Hubble parameter gains numerous advantages over supernovae when $w=-1$, particularly its illuminating power in constraining $\neff$. As long as $w$ is included as a free parameter, it is still the standardizable candles of type Ia supernovae that play the most dominant role in the parameter constraints.

Unresolved near-infrared background anisotropies are expected to have contributions from the earliest galaxies during reionization and faint, dwarf galaxies at intermediate redshifts. Previous measurements were unable to conclusively pinpoint the dominant origin because they did not sample spatial scales that were sufficiently large to distinguish between these two possibilities. Here we report a measurement of the anisotropy power spectrum from sub-arcminute to one degree angular scales and find the clustering amplitude to be larger than the model predictions involving the two existing explanations. As the shot-noise level of the power spectrum is consistent with that expected from faint galaxies, a new source population on the sky is not necessary to explain the observations. A physical mechanism that increases the clustering amplitude, however, is needed. Motivated by recent results related to the extended stellar light profile in dark matter halos, we consider the possibility that the fluctuations originate from diffuse intrahalo stars of all galaxies. We find that the measured power spectrum can be explained by an intrahalo light fraction of 0.07 to 0.2 % relative to the total luminosity in dark matter halos of masses log(M/M_Sun) ~ 9 to 12 at redshifts of ~ 1 to 4.

The ghost free massive gravity modified Friedmann equations at cosmic scale and provided an explanation of cosmic acceleration without dark energy. We analyzed the cosmological solutions of the massive gravity in detail and confronted the cosmological model with current observational data. We found that the model parameters $\alpha_3$ and $\alpha_4$ which are the coefficients of the third and fourth order nonlinear interactions cannot be constrained by current data at the background level. The mass of graviton is found to be the order of current Hubble constant if $\alpha_3=\alpha_4=0$, and the mass of graviton can be as small as possible in the most general case.

A fraction of the extragalactic near-infrared (near-IR) background light involves redshifted photons from the ultraviolet (UV) emission from galaxies present during reionization at redshifts above 6. The absolute intensity and the anisotropies of the near-IR background provide an observational probe of the first-light galaxies and their spatial distribution. We estimate the extragalactic background light intensity during reionization by accounting for the stellar and nebular emission from first-light galaxies. We require the UV photon density from these galaxies to generate a reionization history that is consistent with the optical depth to electron scattering from cosmic microwave background measurements. We also require the bright-end luminosity function of galaxies in our models to reproduce the measured Lyman drop-out luminosity functions at redshifts of 6 to 8. The absolute intensity is about 0.1 to 0.3 nW m$^{-2}$ sr$^{-1}$ at the peak of its spectrum at $\sim$ 1.1 $\mu$m. We also discuss the anisotropy power spectrum of the near-IR background using a halo model to describe the galaxy distribution. We compare our predictions for the anisotropy power spectrum to existing measurements from deep near-IR imaging data from {\it Spitzer}/IRAC, {\it Hubble}/NICMOS, and {\it AKARI}. The predicted rms fluctuations at tens of arcminute angular scales are roughly an order of magnitude smaller than the existing measurements. While strong arguments have been made that the measured fluctuations do not have an origin involving faint low-redshift galaxies, we find that the existing measurements are also incompatible with an origin during the era of reionization. The measured near-IR background anisotropies remain unexplained and could be associated with an unidentified non-astrophysical origin.

Possible bulk compositions of the super-Earth exoplanets, CoRoT-7b, Kepler-9d, and Kepler-10b are investigated by applying a commonly used silicate and a non-standard carbon model. Their internal structures are deduced using the suitable equation of state of the materials. The degeneracy problems of their compositions can be partly overcome, based on the fact that all three planets are extremely close to their host stars. By analyzing the numerical results, we conclude: 1) The iron core of CoRoT-7b is not more than 27% of its total mass within 1 $\sigma$ mass-radius error bars, so an Earth-like composition is less likely, but its carbon rich model can be compatible with an Earth-like core/mantle mass fraction; 2) Kepler-10b is more likely with a Mercury-like composition, its old age implies that its high iron content may be a result of strong solar wind or giant impact; 3) the transiting-only super-Earth Kepler-9d is also discussed. Combining its possible composition with the formation theory, we can place some constraints on its mass and bulk composition.

The wide-area imaging surveys with the {\it Herschel} Space Observatory at sub-mm wavelengths have now resulted in catalogs of order one hundred thousand dusty, star-burst galaxies. We make a statistical estimate of $N(z)$ using a clustering analysis of sub-mm galaxies detected at each of 250, 350 and 500 $\mu$m from the Herschel Multi-tiered Extragalactic Survey (HerMES) centered on the Bo\"{o}tes field. We cross-correlate {\it Herschel} galaxies against galaxy samples at optical and near-IR wavelengths from the Sloan Digital Sky Survey (SDSS), the NOAO Deep Wide Field Survey (NDWFS) and the Spitzer Deep Wide Field Survey (SDWFS). We create optical and near-IR galaxy samples based on their photometric or spectroscopic redshift distributions and test the accuracy of those redshift distributions with similar galaxy samples defined with catalogs of the Cosmological Evolution Survey (COSMOS), as the COSMOS field has superior spectroscopy coverage. We model-fit the clustering auto and cross-correlations of {\it Herschel} and optical/IR galaxy samples to estimate $N(z)$ and clustering bias factors. The $S_{350} > 20$ mJy galaxies have a bias factor varying with redshift as $b(z)=1.0^{+1.0}_{-0.5}(1+z)^{1.2^{+0.3}_{-0.7}}$. This bias and the redshift dependence is broadly in agreement with galaxies that occupy dark matter halos of mass in the range of 10$^{12}$ to 10$^{13}$ M$_{\sun}$. We find that the redshift distribution peaks around $z \sim 0.5$ to 1 for galaxies selected at 250 $\mu$m with an average redshift of $< z > = 1.8 \pm 0.2$. For 350 and 500 $\mu$m-selected SPIRE samples the peak shifts to higher redshift, but the average redshift remains the same with a value of $1.9 \pm 0.2$.

The Herschel Multi-tiered Extragalactic Survey, HerMES, is a legacy program designed to map a set of nested fields totalling ~380 deg^2. Fields range in size from 0.01 to ~20 deg^2, using Herschel-SPIRE (at 250, 350 and 500 \mu m), and Herschel-PACS (at 100 and 160 \mu m), with an additional wider component of 270 deg^2 with SPIRE alone. These bands cover the peak of the redshifted thermal spectral energy distribution from interstellar dust and thus capture the re-processed optical and ultra-violet radiation from star formation that has been absorbed by dust, and are critical for forming a complete multi-wavelength understanding of galaxy formation and evolution. The survey will detect of order 100,000 galaxies at 5\sigma in some of the best studied fields in the sky. Additionally, HerMES is closely coordinated with the PACS Evolutionary Probe survey. Making maximum use of the full spectrum of ancillary data, from radio to X-ray wavelengths, it is designed to: facilitate redshift determination; rapidly identify unusual objects; and understand the relationships between thermal emission from dust and other processes. Scientific questions HerMES will be used to answer include: the total infrared emission of galaxies; the evolution of the luminosity function; the clustering properties of dusty galaxies; and the properties of populations of galaxies which lie below the confusion limit through lensing and statistical techniques. This paper defines the survey observations and data products, outlines the primary scientific goals of the HerMES team, and reviews some of the early results.

The atomic CII fine-structure line is one of the brightest lines in a typical star-forming galaxy spectrum with a luminosity ~ 0.1% to 1% of the bolometric luminosity. It is potentially a reliable tracer of the dense gas distribution at high redshifts and could provide an additional probe to the era of reionization. By taking into account of the spontaneous, stimulated and collisional emission of the CII line, we calculate the spin temperature and the mean intensity as a function of the redshift. When averaged over a cosmologically large volume, we find that the CII emission from ionized carbon in individual galaxies is larger than the signal generated by carbon in the intergalactic medium (IGM). Assuming that the CII luminosity is proportional to the carbon mass in dark matter halos, we also compute the power spectrum of the CII line intensity at various redshifts. In order to avoid the contamination from CO rotational lines at low redshift when targeting a CII survey at high redshifts, we propose the cross-correlation of CII and 21-cm line emission from high redshifts. To explore the detectability of the CII signal from reionization, we also evaluate the expected errors on the CII power spectrum and CII-21 cm cross power spectrum based on the design of the future milimeter surveys. We note that the CII-21 cm cross power spectrum contains interesting features that captures physics during reionization, including the ionized bubble sizes and the mean ionization fraction, which are challenging to measure from 21-cm data alone. We propose an instrumental concept for the reionization CII experiment targeting the frequency range of $\sim$ 200 to 300 GHz with 1, 3 and 10 meter apertures and a bolometric spectrometer array with 64 independent spectral pixels with about 20,000 bolometers.

The large-scale structure of the Universe can be mapped with unresolved intensity fluctuations of the 21 cm line. The power spectrum of the intensity fluctuations has been proposed as a probe of the baryon acoustic oscillations at low to moderate redshifts with interferometric experiments now under consideration. We discuss the contamination to the low-redshift 21 cm intensity power spectrum generated by the 18 cm OH line since the intensity fluctuations of the OH line generated at a slightly higher redshift contribute to the intensity fluctuations observed in an experiment. We assume the OH megamaser luminosity is correlated with the star formation rate, and use the simulation to estimate the OH signal and the spatial anisotropies. We also use a semi-analytic simulation to predict the 21 cm power spectrum. At z=1 to 3, we find that the OH contamination could reach 0.1% to 1% of the 21 cm rms fluctuations at the scale of the first peak of the baryon acoustic oscillation. When z>3 the OH signal declines quickly, so that the contamination on the 21-cm becomes negligible at high redshifts.

A new conception is proposed in Ref.\cite{Verlinde:2010hp,Padmanabhan:2009kr} that gravity is one kind of entropic force. In this letter, we try to discuss its applications to the modified gravities by using three different corrections to the area law of entropy which are derived from the quantum effects and extra dimensions. According to the assumption of holographic principle, the number of bits $N$ which is related to the equipartition law of energy is modified. Then, the modified law of Newton's gravity and the modified Friedmann equations are obtained by using the new notion. By choosing suitable parameters, the modified area law of entropy leads to de-Sitter solutions which can be used to explain the accelerating expansion of our universe. It suggests that the accelerating phase in our universe may be an emergent phenomenon based on holographic principle and thermodynamics.

The cosmological models based on teleparallel gravity with nonzero torsion are considered. To investigate the evolution of this theory, we consider the phase-space analysis of the $f(T)$ theory. It shows when the tension scalar can be written as an inverse function of $x$ where $x=\rho_{e}/(3m_{pl}^{2}H^{2})$ and $T=g(x)$, the system is an autonomous one. Furthermore,the $\omega_{e}-\omega'_{e}$ phase analysis is given out. We perform the dynamical analysis for the models $f(T)=\beta T\ln(T/T_{0})$ and $f(T)=\alpha m_{pl}^{2}(-T/m_{pl}^{2})^{n}$ particularly. We find that the universe will settle into de-Sitter phase for both models. And we have examined the evolution behavior of the power law form in the $\omega_{ep}-\omega'_{ep}$ plane.

We apply the Union2 compilation of 557 supernova Ia data, the baryon acoustic oscillation measurement of distance, the cosmic microwave background radiation data from the seven year Wilkinson Microwave Anisotropy Probe, and the Hubble parameter data to study the geometry of the Universe and the property of dark energy by using models and parametrizations with different high redshift behaviours of $w(z)$. We find that $\Lambda$CDM model is consistent with current data, that the Dvali-Gabadadze-Porrati model is excluded by the data at more than $3\sigma$ level, that the Universe is almost flat, and that the current data is unable to distinguish models with different behaviours of $w(z)$ at high redshift. We also add the growth factor data to constrain the growth index of Dvali-Gabadadze-Porrati model and find that it is more than $1\sigma$ away from its theoretical value.

We first consider the Einstein-aether theory with a gravitational coupling and a Lagrange multiplier field, and then consider the non-minimally coupled quintessence field theory with Lagrange multiplier field. We study the influence of the Lagrange multiplier field on these models. We show that the energy density evolution of the Einstein-aether field and the quintessence field are significantly modified. The energy density of the Einstein-aether is nearly a constant during the entire history of the Universe. The energy density of the quintessence field can also be kept nearly constant in the matter dominated Universe, or even exhibit a phantom-like behavior for some models. This suggests a possible dynamical origin of the cosmological constant or dark energy. Further more, for the canonical quintessence in the absence of gravitational coupling, we find that the quintessence scalar field can play the role of cold dark matter with the introduction of a Lagrange multiplier field. We conclude that the Lagrange multiplier field could play a very interesting and important role in the construction of cosmological models.

We perform a full dynamical analysis by considering the interactions between dark energy and radiation, and dark energy and dark matter. We find that the interaction helps alleviate the coincidence problem for the quintessence model.

We apply the holographic principle and the equipartition law of energy to the apparent horizon of a Friedmann-Robertson-Walker universe and derive the Friedmann equation describing the dynamics of the universe. We also show that the equipartition law of energy can be interpreted as the first law of thermodynamics at the apparent horizon.

Cosmic magnification is due to the weak gravitational lensing of sources in the distant Universe by foreground large-scale structure leading to coherent changes in the observed number density of the background sources. Depending on the slope of the background source number counts, cosmic magnification causes a correlation between the background and foreground galaxies, which is unexpected in the absence of lensing if the two populations are spatially disjoint. Previous attempts using submillimetre (sub-mm) sources have been hampered by small number statistics. The large number of sources detected in the {\it Herschel} Multi-tiered Extra-galactic Survey (HerMES) Lockman-SWIRE field enables us to carry out the first robust study of the cross-correlation between sub-mm sources and sources at lower redshifts. Using ancillary data we compile two low-redshift samples from SDSS and SWIRE with <z> ~ 0.2 and 0.4, respectively, and cross-correlate with two sub-mm samples based on flux density and colour criteria, selecting galaxies preferentially at z ~ 2. We detect cross-correlation on angular scales between ~1 and 50 arcmin and find clear evidence that this is primarily due to cosmic magnification. A small, but non-negligible signal from intrinsic clustering is likely to be present due to the tails of the redshift distribution of the sub-mm sources overlapping with those of the foreground samples.

We propose observations of the molecular gas distribution during the era of reionization. At $z \sim $ 6 to 8, $^{12}$CO$(J=1-0)$ line intensity results in a mean brightness temperature of about 0.5 $\mu$K with rms fluctuations of 0.1 $\mu$K at 1 to 10 Mpc spatial scales, corresponding to 30 arcminute angular scales. This intensity fluctuations can be mapped with an interferometer, similar to existing and planned 21-cm background experiments, but operating at $\sim$ 12 to 17 GHz. We discuss the feasibility of detecting the cross-correlation between HI and CO molecular gas since such a cross-correlation has the advantage that it will be independent of systematics and most foregrounds in each of the 21-cm and CO$(1-0)$ line experiments. Additional instruments tuned to higher-order transitions of the CO molecule or an instrument operating with high spectral resolution at mm-wavelengths targeting 158 $\mu$m CII could further improve the reionization studies with molecular gas. The combined 21-cm and CO line observations has the potential to establish the relative distribution of gas in the inter-galactic medium and molecular gas that are clumped in individual first-light galaxies that are closely connected to the formation of massive stars in these galaxies.

We examine the effectiveness of the weak gravity conjecture in constraining the dark energy by comparing with observations. For general dark energy models with plausible phenomenological interactions between dark sectors, we find that although the weak gravity conjecture can constrain the dark energy, the constraint is looser than that from the observations.

We consider the interaction between dark matter and dark energy in the framework of holographic dark energy, and propose a natural and physically plausible form of interaction, in which the interacting term is proportional to the product of the powers of the dark matter and dark energy densities. We investigate the cosmic evolution in such models. The impact of the coupling on the dark matter and dark energy components may be asymmetric. While the dark energy decouples from the dark matter at late time, just as other components of the cosmic fluid become decoupled as the universe expands, interestingly, the dark matter may actually become coupled to the dark energy at late time. We shall name such a phenomenon as "incoupling". We use the latest type Ia supernovae data from the SCP team, baryon acoustics oscillation data from SDSS and 2dF surveys, and the position of the first peak of the CMB angular power spectrum to constrain the model. We find that the interaction term which is proportional to the the first power product of the dark energy and dark matter densities gives excellent fit to the current data.

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

We study the radion stability and radion mass in the framework of the Horava-Witten (HW) heterotic M-Theory on $S^{1}/Z_{2}$, and find that the radion is stable and its mass can be of the order of GeV. The gravity is localized on the visible brane, and the spectrum of the gravitational Kaluza-Klein (KK) modes is discrete and can have a mass gap of TeV. The corrections to the 4D Newtonian potential from the higher order gravitational KK modes are exponentially suppressed. Applying such a setup to cosmology, we find the generalized Friedmann-like equations on each of the two orbifold branes.

We apply the Constitution compilation of 397 supernova Ia, the baryon acoustic oscillation measurements including the $A$ parameter, the distance ratio and the radial data, the five-year Wilkinson microwave anisotropy probe and the Hubble parameter data to study the geometry of the universe and the property of dark energy by using the popular Chevallier-Polarski-Linder and Jassal-Bagla-Padmanabhan parameterizations. We compare the simple $\chi^2$ method of joined contour estimation and the Monte Carlo Markov chain method, and find that it is necessary to make the marginalized analysis on the error estimation. The probabilities of $\Omega_k$ and $w_a$ in the Chevallier-Polarski-Linder model are skew distributions, and the marginalized $1\sigma$ errors are $\Omega_m=0.279^{+0.015}_{-0.008}$, $\Omega_k=0.005^{+0.006}_{-0.011}$, $w_0=-1.05^{+0.23}_{-0.06}$, and $w_a=0.5^{+0.3}_{-1.5}$. For the Jassal-Bagla-Padmanabhan model, the marginalized $1\sigma$ errors are $\Omega_m=0.281^{+0.015}_{-0.01}$, $\Omega_k=0.000^{+0.007}_{-0.006}$, $w_0=-0.96^{+0.25}_{-0.18}$, and $w_a=-0.6^{+1.9}_{-1.6}$. The equation of state parameter $w(z)$ of dark energy is negative in the redshift range $0\le z\le 2$ at more than $3\sigma$ level. The flat $\Lambda$CDM model is consistent with the current observational data at the $1\sigma$ level.

We use the redshift Hubble parameter $H(z)$ data derived from relative galaxy ages, distant type Ia supernovae (SNe Ia), the Baryonic Acoustic Oscillation (BAO) peak, and the Cosmic Microwave Background (CMB) shift parameter data, to constrain cosmological parameters in the Undulant Universe. We marginalize the likelihood functions over $h$ by integrating the probability density $P\propto e^{-\chi^2/2}$. By using the Markov Chain Monte Carlo (MCMC) technique, we obtain the best fitting results and give the confidence regions on the $b-\Omega_{\rm m0}$ plane. Then we compare their constraints. Our results show that the $H(z)$ data play a similar role with the SNe Ia data in cosmological study. By presenting the independent and joint constraints, we find that the BAO and CMB data play very important roles in breaking the degeneracy compared with the $H(z)$ and SNe Ia data alone. Combined with the BAO or CMB data, one can improve the constraints remarkably. The SNe Ia data sets constrain $\Omega_{\rm m0}$ much tighter than the $H(z)$ data sets, but the $H(z)$ data sets constrain $b$ much tighter than the SNe Ia data sets. All these results show that the Undulant Universe approaches the $\Lambda \rm$CDM model. We expect more $H(z)$ data to constrain cosmological parameters in future.

To realize the accelerations in the early and late periods of our universe, we need to specify potentials for the dominant fields. In this paper, by using the Noether symmetry approach, we try to find suitable potentials in the "cosmic triad" vector field scenario. Because the equation of state parameter of dark energy has been constrained in the range of $-1.21\leq \omega\leq -0.89$ by observations, we derive the Noether conditions for the vector field in quintessence, phantom and quintom models, respectively. In the first two cases, constant potential solutions have been obtained. What is more, a fast decaying point-like solution with power-law potential is also found for the vector field in quintessence model. For the quintom case, we find an interesting constraint $\tilde{C}V_{p}'=-CV_{q}'$ on the field potentials, where $C$ and $\tilde{C}$ are constants related to the Noether symmetry.

In this Letter, we find suitable potentials in the multiple scalar fields scenario by using the Noether symmetry approach. We discussed three models with multiple scalar fields: N-quintessence with positive kinetic terms, N-phantom with negative kinetic terms and N-quintom with both positive and negative kinetic terms. In the N-quintessence case, the exponential potential which could be derived from several theoretic models is obtained from the Noether conditions. In the N-phantom case, the potential $\frac{V_{0}}{2}(1-\cos(\sqrt{\frac{3N}{2}}\frac{\phi}{m_{pl}}))$, which could be derived from the Pseudo Nambu-Goldstone boson model, is chosen as the Noether conditions required. In the N-quintom case, we derive a relation $DV'_{\phi q}=-\tilde{D}V'_{\phi p}$ between the potential forms for the quintessence-like fields and the phantom-like fields by using the Noether symmetry.

The extended holographic dark energy model with the Hubble horizon as the infrared cutoff avoids the problem of the circular reasoning of the holographic dark energy model. Unfortunately, it is hit with the no-go theorem. In this paper, we consider the extended holographic dark energy model with a potential, $V(\phi)$, for the Brans-Dicke scalar field. With the addition of a potential for the Brans-Dicke scalar field, the extended holographic dark energy model using the Hubble horizon as the infrared cutoff is a viable dark energy model, and the model has the dark energy dominated attractor solution.

We use current and future simulated data of the growth rate of large scale structure in combination with data from supernova, BAO, and CMB surface measurements, in order to put constraints on the growth index parameters. We use a recently proposed parameterization of the growth index that interpolates between a constant value at high redshifts and a form that accounts for redshift dependencies at small redshifts. We also suggest here another exponential parameterization with a similar behaviour. The redshift dependent parametrizations provide a sub-percent precision level to the numerical growth function, for the full redshift range. Using these redshift parameterizations or a constant growth index, we find that current available data from galaxy redshift distortions and Lyman-alpha forests is unable to put significant constraints on any of the growth parameters. For example both $\Lambda$CDM and flat DGP are allowed by current growth data. We use an MCMC analysis to study constraints from future growth data, and simulate pessimistic and moderate scenarios for the uncertainties. In both scenarios, the redshift parameterizations discussed are able to provide significant constraints and rule out models when incorrectly assumed in the analysis. The values taken by the constant part of the parameterizations as well as the redshift slopes are all found to significantly rule out an incorrect background. We also find that, for our pessimistic scenario, an assumed constant growth index over the full redshift range is unable to rule out incorrect models in all cases. This is due to the fact that the slope acts as a second discriminator at smaller redshifts and therefore provide a significant test to identify the underlying gravity theory.

We examine in some detail the influence of the systematics in different data sets including type Ia supernova sample, baryon acoustic oscillation data and the cosmic microwave background information on the fitting results of the Chevallier-Polarski-Linder parametrization. We find that the systematics in the data sets does influence the fitting results and leads to different evolutional behavior of dark energy. To check the versatility of Chevallier-Polarski-Linder parametrization, we also perform the analysis on the Wetterich parametrization of dark energy. The results show that both the parametrization of dark energy and the systematics in data sets influence the evolutional behavior of dark energy.

By assuming a phenomenological form for the ratio of the DE and DM densities $\rho_X\propto \rho_{m} a^{\xi}$ (Dalal et al.2001), we discuss the cosmic coincidence problem in light of current observational data. $\xi$ is a key parameter to denote the severity of the coincidence problem. In this scenario, $\xi=3$ and $\xi=0$ correspond to $\Lambda$CDM and the self-similar solution without coincidence problem, respectively. Hence, any solution with a scaling parameter $0<\xi<3$ makes the coincidence problem less severe. In addition, the standard cosmology without interaction between DE and DM is characterized by $\xi+3\omega_X=0$, whereas the inequality $\xi+3\omega_X\neq 0$ represents non-standard cosmology. We place observational constraints on the parameters $(\Omega_{X,0},\omega_{X},\xi)$ of this model, by using the Constitution Set, the CMB shift parameter from the 5-year WMAP and the SDSS baryon acoustic peak. Combining the three samples, we get $\Omega_{X,0}=0.72\pm0.02$, $\omega_{X}=-0.98\pm0.07$ and $\xi=3.06\pm0.35$ at 68.3% confidence level. The result shows that the $\Lambda$CDM model still remains a good fit to the recent observational data, as well as, the coincidence problem indeed exists and is quite severe, in the framework of this simple phenomenological model. We further constrain the model with the transition redshift (deceleration / acceleration). It shows that if the transition from deceleration to acceleration happens at the redshit $z > 0.73$, within the framework of this model, we can conclude that the interaction between DE and DM is necessary.

We use the Constitution supernova, the baryon acoustic oscillation, the cosmic microwave background, and the Hubble parameter data to analyze the evolution property of dark energy. We obtain different results when we fit different baryon acoustic oscillation data combined with the Constitution supernova data to the Chevallier-Polarski-Linder model. We find that the difference stems from the different values of $\Omega_{m0}$. We also fit the observational data to the model independent piecewise constant parametrization. Four redshift bins with boundaries at $z=0.22$, 0.53, 0.85 and 1.8 were chosen for the piecewise constant parametrization of the equation of state parameter $w(z)$ of dark energy. We find no significant evidence for evolving $w(z)$. With the addition of the Hubble parameter, the constraint on the equation of state parameter at high redshift isimproved by 70%. The marginalization of the nuisance parameter connected to the supernova distance modulus is discussed.

We discuss the extent to which photometric measurements alone can be used to identify Type Ia supernovae (SNIa) and to determine redshift and other parameters of interest for cosmological studies. We fit the light curve data of the type expected from a survey such as the one planned with Large Synoptic Survey Telescope (LSST) and also to remove the contamination from the core-collapse supernovae to SNIa samples. We generate 1000 SNIa mock flux data for each of the LSST filters based on existing design parameters, then use a Markov Chain Monte-Carlo (MCMC) analysis to fit for the redshift, apparent magnitude, stretch factor and the phase of the SNIa. We find that the model fitting works adequately well when the true SNe redshift is below 0.5, while at $z < 0.2$ the accuracy of the photometric data is almost comparable with spectroscopic measurements of the same sample. We discuss the contamination of Type Ib/c (SNIb/c) and Type II supernova (SNII) on the SNIa data set. We find it is easy to distinguish the SNII through the large $\chi^2$ mismatch when fitting to photometric data with Ia light curves. This is not the case for SNIb/c. We implement a statistical method based on the Bayesian estimation in order to statistically reduce the contamination from SNIb/c for cosmological parameter measurements from the whole SNe sample. The proposed statistical method also evaluate the fraction of the SNIa in the total SNe data set, which provides a valuable guide to establish the degree of contamination.

The primordial non-Gaussianity (PNG) in matter density perturbation is a very powerful probe of the physics of the very early Universe. The local PNG can induce a distinct scale-dependent bias on the large scale structure distribution of galaxies and quasars, which could be used for constraining it. We study the detection limits on PNG from the surveys of the LAMOST telescope. The cases of the main galaxy survey, the luminous red galaxy (LRG) survey, and the quasar survey of different magnitude limits are considered. We find that the MAIN1 sample (i.e. the main galaxy survey with one magnitude deeper than the SDSS main galaxy survey, or r<18.8) could only provide very weak constraint on PNG. For the MAIN2 sample (r<19.8) and the LRG survey, the 2\sigma (95.5%) limit on the PNG parameter f_{NL} are |f_{NL}|<145 and |f_{NL}|<114 respectively, comparable to the current limit from cosmic microwave background (CMB) data. The quasar survey could provide much more stringent constraint, and we find that the 2\sigma limit for |f_{NL}| is between 50 and 103, depending on the magnitude limit of the survey. With Planck-like priors on cosmological parameters, the quasar survey with g<21.65 would improve the constraints to |f_{NL}|<43 (2\sigma). We also discuss the possibility of further tightening the constraint by using the relative bias method proposed by Seljak(2008).

The five dimensional Brans-Dicke theory naturally provides two scalar fields by the Killing reduction mechanism. These two scalar fields could account for the accelerated expansion of the universe. We test this model and constrain its parameter by using the type Ia supernova (SN Ia) data. We find that the best fit value of the 5-dimensional Brans-Dicke coupling contant is $\omega = -1.9$. This result is also consistent with other observations such as the baryon acoustic oscillation (BAO).

The growth rate of matter perturbation and the expansion rate of the Universe can be used to distinguish modified gravity and dark energy models in explaining cosmic acceleration. We explore here the inclusion of spatial curvature into the growth factor. We expand previous results using the approximation $\Omega_{m}^\gamma$ and then suggest a new form, $f_a=\Omega_m^\gamma+(\gamma-4/7)\Omega_k$, as an approximation for the growth factor when the curvature $\Omega_k$ is not negligible, and where the growth index $\gamma$ is usually model dependent. The expression recovers the standard results for the curved and flat $\Lambda$CDM and Dvali-Gabadadze-Porrati models. Using the best fit values of $\Omega_{m0}$ and $\Omega_{k0}$ to the expansion/distance measurements from Type Ia supernovae, baryon acoustic oscillation, WMAP5, and $H(z)$ data, we fit the growth index parameter to current growth factor data and obtain $\gamma_{\Lambda}(\Omega_{k} \not= 0) = 0.65^{+0.17}_{-0.15}$ and $\gamma_{DGP}(\Omega_{k} \not= 0) = 0.53^{+0.14}_{-0.12}$. For the $\Lambda$CDM model, the 1-$\sigma$ observational bounds are found consistent with theoretical value, unlike the case for the Dvali-Gabadadze-Porrati model. We also find that the current data we used is not enough to put significant constraints when the 3 parameters in $f_a$ are fit simultaneously. Importantly, we find that, in the presence of curvature, the analytical expression proposed for $f_a$ provides a better fit to the growth factor than other forms and should be useful for future high precision missions and studies.

The cosmic coincidence problem is a serious challenge to dark energy model. We suggest a quantitative criteria for judging the severity of the coincidence problem. Applying this criteria to three different interacting models, including the interacting quintessence, interacting phantom, and interacting Chaplygin gas models, we find that the interacting Chaplygin gas model has a better chance to solve the coincidence problem. Quantitatively, we find that the coincidence index C for the interacting Chaplygin gas model is smaller than that for the interacting quintessence and phantom models by six orders of magnitude.

The dynamical behaviors of two interacting dark energy models are considered. In addition to the scaling attractors found in the non-interacting quintessence model with exponential potential, new accelerated scaling attractors are also found in the interacting dark energy models. The coincidence problem is reduced to the choice of parameters in the interacting dark energy models.

We perform a detailed phase-space analysis of various phantom cosmological models, where the dark energy sector interacts with the dark matter one. We examine whether there exist late-time scaling attractors, corresponding to an accelerating universe and possessing dark energy and dark matter densities of the same order. We find that all the examined models, although accepting stable late-time accelerated solutions, cannot alleviate the coincidence problem, unless one imposes a form of fine-tuning in the model parameters. It seems that interacting phantom cosmology cannot fulfill the basic requirement that led to its construction.

We investigate the observational signatures of the holographic dark energy models in this paper, including both the original model and a model with an interaction term between the dark energy and dark matter. We first delineate the dynamical behavior of such models, especially whether they would have "Big Rip" for different parameters, then we use several recent observations, including 182 high-quality type Ia supernovae data observed with the Hubble Space Telescope, the SNLS and ESSENCE surveys, 42 latest Chandra X-ray cluster gas mass fraction, 27 high-redshift gamma-ray burst samples, the baryon acoustic oscillation measurement from the Sloan Digital Sky Survey, and the CMB shift parameter from WMAP three years result to give more reliable and tighter constraints on the holographic dark energy models. The results of our constraints for the holographic dark energy model without interaction is $c=0.748^{+0.108}_{-0.009}$, $\Omega_{\mathrm{m0}}=0.276^{+0.017}_{-0.016}$, and for model with interaction ($c=0.692^{+0.135}_{-0.107}$, $\Omega_{\mathrm{m0}}=0.281^{+0.017}_{-0.017}$ ,$\alpha=-0.006 ^{+0.021}_{-0.024}$, where $\alpha$ is an interacting parameter). As these models have more parameters than the $\Lambda$CDM model, we use the Bayesian evidence as a model selection criterion to make comparison. We found that the holographic dark energy models are mildly favored by the observations compared with the $\mathrm{\Lambda CDM}$ model.

The growth rate of matter perturbation and the expansion rate of the Universe can be used to distinguish modified gravity and dark energy models in explaining the cosmic acceleration. The growth rate is parametrized by the growth index $\gamma$. We discuss the dependence of $\gamma$ on the matter energy density $\Omega$ and its current value $\Omega_0$ for a more accurate approximation of the growth factor. The observational data, including the data of the growth rate, are used to fit different models. The data strongly disfavor the Dvali-Gabadadze-Porrati model. For the dark energy model with a constant equation of state, we find that $\Omega_0=0.27\pm 0.02$ and $w=-0.97\pm 0.09$. For the $\Lambda$CDM model, we find that $\gamma=0.64^{+0.17}_{-0.15}$. For the Dvali-Gabadadze-Porrati model, we find that $\gamma=0.55^{+0.14}_{-0.13}$.

We use the current weak lensing data to constrain the equation of state of dark energy $w$ and the total mass of massive neutrinos $\sum m_{\nu}$. The constraint on $w$ would be weak if only the current weak lensing data are used. With the addition of other observational data such as the type Ia supernovae, baryon acoustic oscillation, and the high redshift Hubble parameter data $H(z)$ derived from relative galaxy ages to break the degeneracy, the result is significantly improved. For the pure $w$CDM model without massive neutrinos, we find $w=-1.00^{+0.10}_{-0.12}$. For the $w$CDM model with the massive neutrino component, we show that the constraint on $w$ is almost unchanged, there is very little degeneracy between $w$ and $\sum m_{\nu}$. After marginalizing over other parameters, we obtain the probability distribution function of $\sum m_{\nu}$, and find that the upper limit is $\sum m_{\nu} \leq 0.8$ eV at 95.5% confidence level for the combined data sets. Our constraints of $w$ and $\sum m_{\nu}$ are both compatible and comparable with the constraints obtained from the WMAP 5-year data.

The extended holographic dark energy model with the Hubble horizon as the infrared cutoff avoids the problem of the circular reasoning of the holographic dark energy model. We show that the infrared cutoff of the extended holographic dark energy model cannot be the Hubble horizon provided that the Brans-Dicke parameter $\omega$ satisfies the experimental constraint $\omega> 10^4$, and this is proved as a no-go theorem. The no-go theorem also applies to the case in which the dark matter interacts with the dark energy.

In unparticle dark matter (unmatter) models the equation of state of the unmatter is given by $p=\rho/(2d_U+1)$, where $d_U$ is the scaling factor. Unmatter with such equations of state would have a significant impact on the expansion history of the universe. Using type Ia supernovae (SNIa), the baryon acoustic oscillation (BAO) measurements and the shift parameter of the cosmic microwave background (CMB) to place constraints on such unmatter models we find that if only the SNIa data is used the constraints are weak. However, with the BAO and CMB shift parameter data added strong constraints can be obtained. For the $\Lambda$UDM model, in which unmatter is the sole dark matter, we find that $d_U > 60$ at 95% C.L. For comparison, in most unparticle physics models it is assumed $d_U<2$. For the $\Lambda$CUDM model, in which unmatter co-exists with cold dark matter, we found that the unmatter can at most make up a few percent of the total cosmic density if $d_U<10$, thus it can not be the major component of dark matter.

We use the Monte-Carlo Markov Chain method to explore the dark energy property and the cosmic curvature by fitting two popular dark energy parameterizations to the observational data. The new 182 gold supernova Ia data and the ESSENCE data both give good constraint on the DE parameters and the cosmic curvature for the dark energy model $w_0+w_a z/(1+z)$. The cosmic curvature is found to be $|\Omega_k|\la 0.03$. For the dark energy model $w_0+w_a z/(1+z)^2$, the ESSENCE data gives better constraint on the cosmic curvature and we get $|\Omega_k|\leq 0.02$.

Orbifold branes are studied in the framework of the 11-dimensional Horava-Witten heterotic M-Theory. It is found that the effective cosmological constant can be easily lowered to its current observational value by the mechanism of large extra dimensions. The domination of this constant over the evolution of the universe is only temporarily. Due to the interaction of the bulk and the branes, the universe will be in its decelerating expansion phase again in the future, whereby all problems connected with a far future de Sitter universe are resolved.

The cold dark matter may be in a meta-stable state and decays to other particles with a very long lifetime. If the decaying products of the dark matter are weakly interacting, e.g. neutrinos, then it would have little impact on astrophysical processes and is therefore difficult to observe. However, such a decay would affect the expansion history of the Universe because of the change of the equation of state. We utilize a high-quality type Ia supernovae (SN Ia) data set selected from several resent observations and the position of the first peak of the Cosmic Microwave Background (CMB) angular spectrum given by the WMAP three-year data to constrain the dark matter decay-to-neutrino rate $\Gamma=\alpha \Gamma_{\chi}$, where $\alpha$ is the fraction of the rest mass which gets converted to neutrinos, and $\Gamma_{\chi}$ is the decay width. We find that $\Gamma^{-1} > 0.7\times10^3$ Gyr at 95.5% confidence level.

Although there is mounting observational evidence that the cosmic expansion is undergoing a late-time acceleration, the physical mechanism behind such a phenomenon is yet unknown. In this paper, we investigate a holographic dark energy (HDE) model with interaction between the components of the dark sector in the light of current cosmological observations. We use both the new \emph{gold} sample of 182 type Ia supernovae (SNe Ia) and the 192 SNe Ia ESSENCE data, the baryon acoustic oscillation measurement from the Sloan Digital Sky Survey and the shift parameter from the three-year Wilkinson Microwave Anisotropy Probe data. In agreement with previous results, we show that these observations suggest a very weak coupling or even a noninteracting HDE. The phantom crossing behavior in the context of these scenarios is also briefly discussed.

With the help of a masslike function which has dimension of energy and equals to the Misner-Sharp mass at the apparent horizon, we show that the first law of thermodynamics of the apparent horizon $dE=T_AdS_A$ can be derived from the Friedmann equation in various theories of gravity, including the Einstein, Lovelock, nonlinear, and scalar-tensor theories. This result strongly suggests that the relationship between the first law of thermodynamics of the apparent horizon and the Friedmann equation is not just a simple coincidence, but rather a more profound physical connection.

We consider the possibility of observing the onset of the late time inflation of our patch of the Universe. The Hubble size criterion and the event horizon criterion are applied to several dark energy models to discuss the problem of future inflation of the Universe. We find that the acceleration has not lasted long enough to confirm the onset of inflation by present observations for the dark energy model with constant equation of state, the holographic dark energy model and the generalized Chaplygin gas (GCG) model. For the flat $\Lambda$CDM model with $\Omega_{m0}=0.3$, we find that if we use the Hubble size criterion, we need to wait until the $a_v$ which is the scale factor at the time when the onset of inflation is observed reaches 3.59 times of the scale factor $a_T$ when the Universe started acceleration, and we need to wait until $a_v=2.3 a_T$ to see the onset of inflation if we use the event horizon criterion. For the flat holographic dark energy model with $d=1$, we find that $a_v=3.46 a_T$ with the Hubble horizon and $a_v=2.34 a_T$ with the event horizon, respectively. For the flat GCG model with the best supernova fitting parameter $\alpha=1.2$, we find that $a_v=5.50 a_T$ with the Hubble horizon and $a_v=2.08 a_T$ with the event horizon, respectively.

We consider the possibility that the dark energy is made up of two or more independent components, each having a different equation of state. We fit the model with supernova and gamma-ray burst (GRB) data from resent observations, and use the Markov Chain Monte Carlo (MCMC) technique to estimate the allowed parameter regions. We also use various model selection criteria to compare the two component model with the LCDM, one component dark energy model with static or variable w(XCDM), and with other multi-component models. We find that the two component models can give reasonably good fit to the current data. For some data sets, and depending somewhat on the model selection criteria, the two component model can give better fit to the data than XCDM with static w and XCDM with variable w parameterized by w = w_0 + w_az/(1+z).

Using the data coming from the new 182 Gold type Ia supernova samples, the shift parameter of the Cosmic Microwave Background given by the three-year Wilkinson Microwave Anisotropy Probe observations, and the baryon acoustic oscillation measurement from the Sloan Digital Sky Survey, $H(z)$ and lookback time measurements, we have performed a statistical joint analysis of the interacting holographic dark energy model. Consistent parameter estimations show us that the interacting holographic dark energy model is a viable candidate to explain the observed acceleration of our universe.

The energy conditions give upper bounds on the luminosity distance. We apply these upper bounds to the 192 essence supernova Ia data to show that the Universe had experienced accelerated expansion. This conclusion is drawn directly from the distance modulus-reshift graph. In addition to be a very simple method, this method is also totally independent of any cosmological model. From the degeneracy of the distance modulus at low redshift, we argue that the choice of $w_0$ for probing the property of dark energy is misleading. One explicit example is used to support this argument.

The energy conditions provide a very promising model-independent study of the current acceleration of the universe. However, in order to connect these conditions with observations, one often needs first to integrate them, and then find the corresponding constraints on some observational variables, such as the distance modulus. Those integral forms can be misleading, and great caution is needed when one interprets them physically. A typical example is that the transition point of the deceleration parameter $q(z)$ is at about $z \simeq 0.76$ in the $\Lambda$CDM model. However, with the same model when we consider the dimensionless Hubble parameter $E(z)$, which involves the integration of $q(z)$, we find that $E(z)$ does not cross the line of $q(z) = 0$ before $z = 2$. Therefore, to get the correct result, we cannot use the latter to determine the transition point. With these in mind, we carefully study the constraints from the energy conditions, and find that, among other things, the current observational data indeed strongly indicate that our universe has ocne experienced an accelerating expansion phase between the epoch of galaxy formation and the present.

We have investigated the thermodynamical properties of dark energy. Assuming that the dark energy temperature $T\sim a^{-n}$ and considering that the volume of the Universe enveloped by the apparent horizon relates to the temperature, we have derived the dark energy entropy. For dark energy with constant equation of state $w>-1$ and the generalized Chaplygin gas, the derived entropy can be positive and satisfy the entropy bound. The total entropy, including those of dark energy, the thermal radiation and the apparent horizon, satisfies the generalized second law of thermodynamics. However, for the phantom with constant equation of state, the positivity of entropy, the entropy bound, and the generalized second law cannot be satisfied simultaneously.

The new 182 gold supernova Ia data, the baryon acoustic oscillation measurement and the shift parameter determined from the Sloan Digital Sky Survey and the three-year Wilkinson Microwave Anisotropy Probe data are combined to reconstruct the dark energy equation of state parameter $w(z)$ and the deceleration parameter $q(z)$. We find that the strongest evidence of acceleration happens around the redshift $z\sim 0.2$ and the stringent constraints on $w(z)$ lie in the redshift range $z\sim 0.2-0.5$. At the sweet spot, $-1.2<w(z)<-0.6$ for the dark energy parametrization $w(z)=w_0+w_a z/(1+z)^2$ at the $3\sigma$ confidence level. The transition redshift $z_t$ when the Universe underwent the transition from deceleration to acceleration is derived to be $z_t=0.36^{+0.23}_{-0.08}$. The combined data is also applied to find out the geometry of the Universe, and we find that at the $3\sigma$ confidence level, $|\Omega_k|\alt 0.05$ for the simple one parameter dark energy model, and $-0.064<\Omega_k<0.028$ for the $\Lambda$CDM model.

We have investigated the thermodynamical properties of the Universe with dark energy. Adopting the usual assumption in deriving the constant co-moving entropy density that the physical volume and the temperature are independent, we observed some strange thermodynamical behaviors. However, these strange behaviors disappeared if we consider the realistic situation that the physical volume and the temperature of the Universe are related. Based on the well known correspondence between the Friedmann equation and the first law of thermodynamics of the apparent horizon, we argued that the apparent horizon is the physical horizon in dealing with thermodynamics problems. We have concentrated on the volume of the Universe within the apparent horizon and considered that the Universe is in thermal equilibrium with the Hawking temperature on the apparent horizon. For dark energy with $w\ge -1$, the holographic principle and the generalized second law are always respected.

We discuss the phenomenological model in which the potential energy of the quintessence field depends linearly on the energy density of the spatial curvature. We find that the pressure of the scalar field takes a different form when the potential of the scalar field also depends on the scale factor and the energy momentum tensor of the scalar field can be expressed as the form of a perfect fluid. A general coupling was proposed to explain the current accelerating expansion of the Universe and solve the fine-tuning problem.

The compactification of M theory with time dependent hyperbolic internal space gives an effective scalar field with exponential potential which provides a transient acceleration in Einstein frame in four dimensions. Ordinary matter and radiation are present in addition to the scalar field coming from compactification. We find that we have to fine-tune the initial conditions of the scalar field so that our Universe experiences acceleration now. During the evolution of our Universe, the volume of the internal space increases about 12 times. The time variation of the internal space results in a large time variation of the fine structure constant which violates the observational constraint on the variation of the fine structure constant. The large variation of the fine structure constant is a generic feature of transient acceleration models.

Two new one-parameter tracking behavior dark energy representations $\omega=\omega_0/(1+z)$ and $\omega=\omega_0 e^{z/(1+z)}/(1+z)$ are used to probe the geometry of the Universe and the property of dark energy. The combined type Ia supernova (SN Ia), Sloan Digital Sky Survey (SDSS), and Wilkinson Microwave Anisotropy Probe (WMAP) data indicate that the Universe is almost spatially flat and that dark energy contributes about 72% of the matter content of the present universe. The observational data also tell us that $\omega(0)\sim -1$. It is argued that the current observational data can hardly distinguish different dark energy models to the zeroth order. The transition redshift when the expansion of the Universe changed from deceleration phase to acceleration phase is around $z_{\rm T}\sim 0.6$ by using our one-parameter dark energy models.

We show that the most general dark energy model that possesses a scaling solution $\rho_\phi\propto a^n$ is the k-essence model, which includes both of the quintessence and tachyon models. The exact scaling solutions are then derived. The potential that gives the tracking solution in which dark energy exactly tracks the background matter field is the inverse squared potential. The quintessence field with exponential potential can be obtained from the k-essence field with the inverse squared potential. We also find the fixed points and study their main properties, whereby the scalar field dominant fixed point is identified.

We propose a new parametrization of the deceleration parameter to study its time-variation behavior. The advantage of parameterizing the deceleration parameter is that we do not need to assume any underlying theory of gravity. By fitting the model to the 157 gold sample supernova Ia data, we find strong evidence that the Universe is currently accelerating and it accelerated in the past. By fitting the model to the 115 nearby and Supernova Legacy Survey supernova Ia data, the evidence that the Universe is currently accelerating is weak, although there is still a strong evidence that the Universe once accelerated in the past. The results obtained from the 157 gold sample supernova Ia data and those from the 115 supernova Ia data are not directly comparable because the two different data sets measure the luminosity distance up to different redshifts. We then use the Friedmann equation and a dark energy parametrization to discuss the same problem. When we fit the model to the supernova Ia data alone, we find weak evidence that the Universe is accelerating and the current matter density is higher than that measured from other experiments. After we add the Sloan Digital Sky Survey data to constrain the dark energy model, we find that the behavior of the deceleration parameter is almost the same as that obtained from parameterizing the deceleration parameter.

We consider a brane cosmological model with energy exchange between brane and bulk. Parameterizing the energy exchange term by the scale factor and Hubble parameter, we are able to exactly solve the modified Friedmann equation on the brane. In this model, the equation of state for the effective dark energy has a transition behavior changing from $w_{de}^{eff}>-1$ to $w_{de}^{eff}<-1$, while the equation of state for the dark energy on the brane has $w>-1$. Fitting data from type Ia supernova, Sloan Digital Sky Survey and Wilkinson Microwave Anisotropy Probe, our universe is predicted now in the state of super-acceleration with $w_{de0}^{eff}=-1.21$.

The holographic principle is used to discuss the holographic dark energy model. We find that the Bekenstein-Hawking entropy bound is far from saturation under certain conditions. A more general constraint on the parameter of the holographic dark energy model is also derived.

We have reexamined the holographic dark energy model by considering the spatial curvature. We have refined the model parameter and observed that the holographic dark energy model does not behave as phantom model. Comparing the holographic dark energy model to the supernova observation alone, we found that the closed universe is favored. Combining with the Wilkinson Microwave Anisotropy Probe (WMAP) data, we obtained the reasonable value of the spatial curvature of our universe.

A model of holographic dark energy with an interaction with matter fields has been investigated. Choosing the future event horizon as an IR cutoff, we have shown that the ratio of energy densities can vary with time. With the interaction between the two different constituents of the universe, we observed the evolution of the universe, from early deceleration to late time acceleration. In addition, we have found that such an interacting dark energy model can accommodate a transition of the dark energy from a normal state where $w_D>-1$ to $w_D<-1$ phantom regimes. Implications of interacting dark energy model for the observation of dark energy transition has been discussed.

The recent observations of type Ia supernovae strongly support that the universe is accelerating now and decelerated in the recent past. This may be the evidence of the breakdown of the standard Friemann equation. We consider a general modified Friedmann equation. Three different models are analyzed in detail. The current supernovae data and the Wilkinson microwave anisotropy probe data are used to constrain these models. A detailed analysis of the transition from the deceleration phase to the acceleration phase is also performed.

The nature of dark energy is a mystery to us. This paper uses the supernova data to explore the property of dark energy by some model independent methods. We first Talyor expanded the scale factor $a(t)$ to find out the deceleration parameter $q_0<0$. This result just invokes the Robertson-Walker metric. Then we discuss several different parameterizations used in the literature. We find that $\Omega_{\rm DE0}$ is almost less than -1 at $1\sigma$ level. We also find that the transition redshift from deceleration phase to acceleration phase is $z_{\rm T}\sim 0.3$.

The recent observations support that our universe is flat and expanding with acceleration. A quintessence model with a general relation between the quintessence potential and the quintessence kinetic energy was proposed to explain the phenomenon. The dark energy potential includes both the hyperbolic and the double exponential potentials. We analyze this model in detail by using the recent supernova and the first year Wilkinson Microwave Anisotropy Probe (WMAP) observations. For a flat universe dominated by a dark energy with constant $\omega$ which is a special case of the general model, we find that $\Omega_{\rm m0}=0.30^{+0.06}_{-0.08}$ and $\omega_{\rm Q}\le -0.82$, and the turnaround redshift $z_{\rm T}$ when the universe switched from the deceleration phase to the acceleration phase is $z_{\rm T}=0.65$. For the general model, we find that $\Omega_{\rm m0}\sim 0.3$, $\omega_{\rm Q0}\sim -1.0$, $\beta\sim 0.5$ and $z_{\rm T}\sim 0.67$. A model independent polynomial parameterization is also considered, the best fit model gives $\Omega_{\rm m0}=0.40\pm 0.14$, $\omega_{\rm Q0}=-1.4$ and $z_{\rm T}=0.37$.

The recent type Ia supernova data suggest that the universe is accelerating now and decelerated in recent past. This may provide the evidence that the standard Friedmann equation needs to be modified. We analyze in detail a new model in the context of modified Friedmann equation using the supernova data published by the High-$z$ Supernova Search Team and the Supernova Cosmology Project. The new model explains recent acceleration and past deceleration. Furthermore, the new model also gives a decelerated universe in the future.

The generalized Chaplygin gas model with parameter space $\alpha>-1$ is studied in this paper. Some reasonable physical constraints are added to justify the use of the larger parameter space. The Type Ia supernova data and age data of some clusters are then used to fit the model. We find that the parameters have bimodal distributions. For the generalized Chaplygin gas model, we also find that less free parameters fit the data better. The best fit model is the spatially flat model with baryons. The best fit parameters are: $\Omega_{m0}=0.044$, $w_{c0}=-0.881$ and $\alpha=1.57$. The transition redshift is $z_T=0.395$.

We study the possibility to probe the spatial geometry of the Universe by supernova measurement of the cubic correction to the luminosity distance. We illustrate with an accelerating universe model with infinite-volume extra dimensions, for which the 1$\sigma$ level supernova results indicate that the Universe is closed.

The idea of relating the infrared and ultraviolet cutoffs is applied to Brans-Dicke theory of gravitation. We find that extended holographic dark energy from the Hubble scale or the particle horizon as the infrared cutoff will not give accelerating expansion. The dynamical cosmological constant with the event horizon as the infrared cutoff is a viable dark energy model.

In this paper, we use the type Ia supernova data to constrain the model of holographic dark energy. For $d=1$, the best fit result is $\Omega^0_m = 0.25$, the equation of the state of the holographic dark energy $w^0_\d = -0.91$ and the transition between the decelerating expansion and the accelerating expansion happened when the cosmological red-shift was $z_T=0.72$. If we set $d$ as a free parameter, the best fit results are $d=0.21$, $\Omega_m^0=0.46$, $w_\d^0=-2.67$, which sounds like a phantom today, and the transition redshift is $z_T=0.28$.

The recent observations of type Ia supernovae strongly support that the universe is accelerating now and decelerated in the recent past. This may be the evidence of the breakdown of the standard Friedmann equation. Instead of a linear function of the matter density, we consider a general function of the matter density to modify the Freidmann equation. We propose a new model which explains the recent acceleration and the past deceleration. Furthermore, the new model also gives a decelerated universe in the future. The new model gives $\Omega_{m0}=0.46$ and $z_T=0.44$.

The recent observations of type Ia supernovae suggest that the universe is accelerating now and decelerated in the recent past. This may be the evidence of the breakdown of the standard Friedmann equation. The Friedmann equation $H^2\sim \rho$ is modified to be a general form $H^2=g(\rho)$. Three models with particular form of $g(\rho)$ are considered in detail. The supernova data published by \citet{tonry}, \citet{ddsg} and \citet{raknop03} are used to analyze the models. After the best fit parameters are obtained, we then find out the transition redshift $z_{\rm T}$ when the universe switched from the deceleration phase to the acceleration phase.