Standard siren cosmology of gravitational wave (GW) merger events relies on the identification of host galaxies and their redshifts. But this can be highly challenging due to numerous candidates of galaxies in the GW localization area. We point out that the number of candidates can be reduced by orders of magnitude for strongly lensed GW events, due to extra observational constraints. For the next-generation GW detectors like Einstein Telescope (ET), we estimate that this number is usually significantly less than one, as long as the GW localization uncertainty is better than $\sim 10\, \rm deg^2$. This implies that the unique identification of the host galaxy of lensed GW event detected by ET and Cosmic Explorer (CE) is possible. This provides us a promising opportunity to measure the redshift of the GW event and facilitate the standard siren cosmology. We also discuss its potential applications in understanding the evolution process and environment of the GW event.

The integrated Sachs-Wolfe (ISW) effect is caused by the decay of cosmological gravitational potential, and is therefore a unique probe of dark energy. However, its robust detection is still problematic. Various tensions between different data sets, different large scale structure (LSS) tracers, and between data and the $\Lambda$CDM theory prediction, exist. We propose a novel method of ISW measurement by cross correlating CMB and the LSS traced by "low-density-position" (LDP, \citet{2019ApJ...874....7D}). It isolates the ISW effect generated by low-density regions of the universe, but insensitive to selection effects associated with voids. We apply it to the DR8 galaxy catalogue of the DESI Legacy imaging surveys, and obtain the LDPs at $z\leq 0.6$ over $\sim$ 20000 $deg^2$ sky coverage. We then cross correlate with the Planck temperature map, and detect the ISW effect at $3.2\sigma$. We further compare the measurement with numerical simulations of the concordance $\Lambda$CDM cosmology, and find the ISW amplitude parameter $A_{ISW}=1.14\pm0.38$ when we adopt a LDP definition radius $R_s=3^{'}$, fully consistent with the prediction of the standard $\Lambda$CDM cosmology ($A_{ISW}=1$). This agreement with $\Lambda$CDM cosmology holds for all the galaxy samples and $R_s$ that we have investigated. Furthermore, the S/N is comparable to that of galaxy ISW measurement. These results demonstrate the LDP method as a competitive alternative to existing ISW measurement methods, and provide independent checks to existing tensions.

The disagreement between direct late-time measurements of the Hubble constant from the SH0ES collaboration, and early-universe measurements based on the $\Lambda$CDM model from the Planck collaboration might, at least in principle, be explained by new physics in the early universe. Recently, the application of the Effective Field Theory of Large-Scale Structure to the full shape of the power spectrum of the SDSS/BOSS data has revealed a new, rather powerful, way to measure the Hubble constant and the other cosmological parameters from Large-Scale Structure surveys. In light of this, we analyze two models for early universe physics, Early Dark Energy and Rock 'n' Roll, that were designed to significantly ameliorate the Hubble tension. Upon including the information from the full shape to the Planck, BAO, and Supernovae measurements, we find that the degeneracies in the cosmological parameters that were introduced by these models are well broken by the data, so that these two models do not significantly ameliorate the tension.

The radioheliograph image is essential for the study of solar short term activities and long term variations, while the continuity and granularity of radioheliograph data is not so ideal, due to the short visible time of the sun and the complex electron-magnetic environment near the ground-based radio telescope. In this work, we develop a multi-channel input single-channel output neural network, which can generate radioheliograph image in microwave band from the Extreme Ultra-violet (EUV) observation of the Atmospheric Imaging Assembly (AIA) on-board the Solar Dynamic Observatory (SDO). The neural network is trained with nearly 8 years of data of Nobeyama Radioheliograph (NoRH) at 17 GHz and SDO/AIA from January 2011 to September 2018. The generated radioheliograph image is in good consistency with the well-calibrated NoRH observation. SDO/AIA provides solar atmosphere images in multiple EUV wavelengths every 12 seconds from space, so the present model can fill the vacancy of limited observation time of microwave radioheliograph, and support further study of the relationship between the microwave and EUV emission.

Non-Gaussianity of the cosmological matter density field can be largely reduced by a local Gaussianization transformation (and its approximations such as the logrithmic transformation). Such behavior can be recasted as the Gaussian copula hypothesis, and has been verified to very high accuracy at two-point level. On the other hand, statistically significant non-Gaussianities in the Gaussianized field have been detected in simulations. We point out that, this apparent inconsistency is caused by the very limited degrees of freedom in the copula function, which make it misleading as a diagnosis of residual non-Gaussianity in the Gaussianized field. Using the copula density, we highlight the departure from Gaussianity. We further quantify its impact in the predicted n-point correlation functions. We explore a remedy of the Gaussian copula hypothesis, which alleviates but not completely solves the above problems.

The Jiangmen Underground Neutrino Observatory is proposed to determine neutrino mass hierarchy using a 20~ktonne liquid scintillator detector. Strict radio-purity requirements have been put forward for all the components of the detector. According to the MC simulation results, the radon dissolved in the water Cherenkov detector should be below 200~mBq/m$^3$. Radium, the progenitor of radon, should also be taken seriously into account. In order to measure the radium concentration in water, a radium measurement system, which consists of a radium extraction system, a radon emanation chamber and a radon concentration measurement system, has been developed. In this paper, the updated radon concentration in gas measurement system with a one-day-measurement sensitivity of $\sim$5~mBq/m$^3$, the detail of the development of the radium concentration in water measurement system with a sensitivity of $\sim$23~mBq/m$^3$ as well as the measurement results of Daya Bay water samples will be presented.

Context. Solar radio bursts originate mainly from high energy electrons accelerated in solar eruptions like solar flares, jets, and coronal mass ejections. A subcategory of solar radio bursts with short time duration may be used as a proxy to understand the wave generation and propagation within the corona. Aims. Complete case studies of the source size, position and kinematics of short term bursts are very rare due to instrumental limitations. A comprehensive multi-frequency spectroscopic and imaging study was carried out of a clear example of a solar type IIIb-III pair. Methods. In this work, the source of the radio burst was imaged with the interferometric mode, using the remote baselines of the LOw Frequency ARray (LOFAR). A detailed analysis of the fine structures in the spectrum and of the radio source motion with imaging was conducted. Results. The study shows how the fundamental and harmonic components have a significantly different source motion. The apparent source of the fundamental emission at 26.56MHz displaces away from the solar disk center at about 4 times the speed of light, while the apparent source of the harmonic emission at the same frequency shows a speed of < 0.02c. The source size of the harmonic emission, observed in this case, is smaller than that in previous studies, indicating the importance of the use of the remote baselines.

To reach the full potential for the next generation of weak lensing surveys, it is necessary to mitigate the contamination of intrinsic alignments (IA) of galaxies in the observed cosmic shear signal. The self calibration (SC) of intrinsic alignments provides an independent method to measure the IA signal from the survey data and the photometric redshift information. It operates differently from the marginalization method based on IA modeling. In this work, we present the first application of SC to the KiDS450 data and the KV450 data, to split directly the intrinsic shape - galaxy density (Ig) correlation signal and the gravitational shear - galaxy density (Gg) correlation signal, using the information from photometric redshift (photo-z). We achieved a clear separation of the two signals and performed several validation tests. Our measured signals are found to be in general agreement with the KiDS450 cosmic shear best-fit cosmology, for both lensing and IA measurements. For KV450, we use partial (high-z) data, and our lensing measurements are also in good agreement with KV450 cosmic shear best-fit, however, our IA signal suggests a larger IA amplitude. We discussed the impact of photo-z quality on IA detection and several other potential systematic biases. Finally, we discuss the potential application of the information extracted for both the lensing signal and the IA signal in future surveys.

The most commonly used estimators of the anisotropic galaxy power spectrum employ Fast Fourier transforms, and rely on a specific choice of the line-of-sight that breaks the symmetry between the galaxy pair. This leads to wide-angle effects, including the presence of odd power spectrum multipoles like the dipole ($\ell = 1$) and octopole ($\ell = 3$). In Fourier-space these wide-angle effects also couple to the survey window function. We present a self-consistent framework extending the commonly used window function treatment to include the wide-angle effects. We show that our framework can successfully model the wide-angle effects in the BOSS DR12 dataset. We present estimators for the odd power spectrum multipoles and, detect these multipoles in BOSS DR12 with high significance. Understanding the impact of the wide-angle effects on the power spectrum multipoles is essential for many cosmological observables like primordial non-Gaussianity and the detection of General Relativistic effects and represents a potential systematic for measurements of Baryon Acoustic Oscillations and redshift-space distortions.

The Five-Hundred-Meter Aperture Spherical Radio Telescope(FAST) is the largest single-dish radio telescope in the world. In this paper, we make forecast on the FAST HI large scale structure survey by mock observations. We consider a drift scan survey with the L-band 19 beam receiver, which may be commensal with the pulsar search and Galactic HI survey. We also consider surveys at lower frequency, either using the current single feed wide band receiver, or a future multi-beam phased array feed (PAF) in the UHF band. We estimate the number density of detected HI galaxies and the measurement error in positions, the precision of the surveys are evaluated using both Fisher matrix and simulated observations. The measurement error in the HI galaxy power spectrum is estimated, and we find that the error is relatively large even at moderate redshifts, as the number of positively detected galaxies drops drastically with increasing redshift. However, good cosmological measurement could be obtained with the intensity mapping technique where the large scale HI distribution is measured without resolving individual galaxies. The figure of merit (FoM) for the dark energy equation of state with different observation times are estimated, we find that with the existing L-band multi-beam receiver, a good measurement of low redshift large scale structure can be obtained, which complements the existing optical surveys. With a PAF in the UHF band, the constraint can be much stronger, reaching the level of a dark energy task force (DETF) stage IV experiment.

In this work, we propose a powerful probe of neutrino effects on the large-scale structure (LSS) of the Universe, i.e., Minkowski functionals (MFs). The morphology of LSS can be fully described by four MFs. This tool, with strong statistical power, is robust to various systematics and can comprehensively probe all orders of N-point statistics. By using a pair of high-resolution N-body simulations, for the first time, we comprehensively studied the subtle neutrino effects on the morphology of LSS. For an ideal LSS survey of volume $\sim1.73$ Gpc$^3$/$h^3$, neutrino signals are mainly detected from void regions with a significant level up to $\thicksim 10\sigma$ and $\thicksim 300\sigma$ for CDM and total matter density fields, respectively. This demonstrates its enormous potential for much improving the neutrino mass constraint in the data analysis of up-coming ambitious LSS surveys.

An accurate theoretical template for the galaxy power spectrum is a key for the success of ongoing and future spectroscopic surveys. We examine to what extent the Effective Field Theory of Large Scale Structure is able to provide such a template and correctly estimate cosmological parameters. To that end, we initiate a blinded challenge to infer cosmological parameters from the redshift-space power spectrum of high-resolution mock catalogs mimicking the BOSS galaxy sample but covering a hundred times larger cumulative volume. This gigantic simulation volume allows us to separate systematic bias due to theoretical modeling from the statistical error due to sample variance. The challenge task was to measure three unknown input parameters used in the simulation: the Hubble constant, the matter density fraction, and the clustering amplitude. We present analyses done by two independent teams, who have fitted the mock simulation data generated by yet another independent group. This allows us to avoid any confirmation bias by analyzers and pin down possible tuning of the specific EFT implementations. Both independent teams have recovered the true values of the input parameters within sub-percent statistical errors corresponding to the total simulation volume.

We apply the Effective Field Theory of Large-Scale Structure to analyze the $w$CDM cosmological model. By using the full shape of the power spectrum and the BAO post-reconstruction measurements from BOSS, the Supernovae from Pantheon, and a prior from BBN, we set the competitive CMB-independent limit $w=-1.046_{-0.052}^{+0.055}$ at $68\%$ C.L.. After adding the Planck CMB data, we find $w=-1.023_{-0.030}^{+0.033}$ at $68\%$ C.L.. Our results are obtained using PyBird, a new, fast Python-based code which we make publicly available.

Galaxy intrinsic alignment (IA) is both a source of systematic contamination to cosmic shear measurement and its cosmological applications, and a source of valuable information on the large scale structure of the universe and galaxy formation. The self-calibration (SC) method \citep{SC2008} was designed to separate IA from cosmic shear, free of IA modeling. It was first successfully applied to the KiDS450 and KV450 data \citep{Yao2019}. We apply the SC method to the DECaLS DR3 shear + photo-z catalog and significantly improve the IA detection to $\sim 20\sigma$. We find a strong dependence of IA on galaxy color, with strong IA signal for red galaxies, and weak/vanishing IA signal for blue galaxies. The detected IA for red galaxies is in reasonable agreement with the non-linear tidal alignment model and the inferred IA amplitude increases with redshift. We discuss various caveats and possible improvements in the measurement, theory and parameter fitting that will be addressed in future works.

We report on our analysis of the 2019 outburst of the X-ray accreting pulsar 4U 1901+03 observed with Insight-HXMT and NICER. Both spectra and pulse profiles evolve significantly in the decaying phase of the outburst. Dozens of flares are observed throughout the outburst. They are more frequent and brighter at the outburst peak. We find that the flares, which have a duration from tens to hundreds of seconds, are generally brighter than the persistent emission by a factor of $\sim$ 1.5. The pulse profile shape during the flares can be significantly different than that of the persistent emission. In particular, a phase shift is clearly observed in many cases. We interpret these findings as direct evidence of changes of the pulsed beam pattern, due to transitions between the sub- and super-critical accretion regimes on a short time scale. We also observe that at comparable luminosities the flares' pulse profiles are rather similar to those of the persistent emission. This indicates that the accretion on the polar cap of the neutron star is mainly determined by the luminosity, i.e., the mass accretion rate.

Solar S-bursts are short duration ($<1$ s at decameter wavelengths) radio bursts that have been observed during periods of moderate solar activity, where S stands for short. The frequency drift of S-bursts can reflect the density variation and the motion state of the electron beams. In this work, we investigate the frequency drift and the fine structure of the S-bursts with the LOw Frequency ARray (LOFAR). We find that the average frequency drift rate of the S-bursts within 110-180MHz could be described by $df/dt=-0.0077f^{1.59}$. With the high time and frequency resolution of LOFAR, we can resolve the fine structures of the observed solar S-bursts. A fine drift variation pattern was found in the structure of S-bursts (referred to as solar Sb-bursts in this paper) during the type-III storm on 2019 April 13, in the frequency band of 120-240 MHz. The Sb-bursts have a quasi-periodic segmented pattern, and the relative flux intensity tends to be large when the frequency drift rate is relatively large. This kind of structure exists in about 20\% of the solar S-burst events within the observed frequency range. We propose that the fine structure is due to the density fluctuations of the background coronal density. We performed a simulation based on this theory which can reproduce the shape and relative flux intensity of the Sb-bursts. This work shows that the fine structure of solar radio bursts can be used to diagnose the coronal plasma.

We perform a systematic search of quasi-periodic variabilities in $\gamma$-ray active galactic nuclei (AGNs) in the third \emph{Fermi} Large Area Telescope source catalog (3FGL). We employe two techniques, Lomb-Scargle Periodogram (LSP) and Weighted Wavelet Z-transform (WWZ), to obtain power spectra of $\gamma$-ray light curves covering from 2008 August to 2016 December. The results show that besides several objects that have been reported in previous works, an additional source, the FSRQ PKS 0601-70 has a possible quasi-periodic variability of 450 days in its $\gamma$-ray light curves with the significance of $>3\sigma$. The physical implications of our findings are discussed.

A second peak in extreme ultraviolet sometimes appears during the gradual phase of solar flares, which is known as EUV late phase (ELP). Stereotypically ELP is associated with two separated sets of flaring loops with distinct sizes, and it has been debated whether ELP is caused by additional heating or extended plasma cooling in the longer loop system. Here we carry out a survey of 55 M-and-above GOES-class flares with ELP during 2010--2014. Based on the flare-ribbon morphology, these flares are categorized as circular-ribbon (19 events), two-ribbon (23 events), and complex-ribbon (13 events) flares. Among them, 22 events (40%) are associated with coronal mass ejections, while the rest are confined. An extreme ELP, with the late-phase peak exceeding the main-phase peak, is found in 48% of two-ribbon flares, 37% of circular-ribbon flares, and 31% of complex-ribbon flares, suggesting that additional heating is more likely present during ELP in two-ribbon than in circular-ribbon flares. Overall, cooling may be the dominant factor causing the delay of the ELP peak relative to the main-phase peak, because the loop system responsible for the ELP emission is generally larger than, and well separated from, that responsible for the main-phase emission. All but one of the circular-ribbon flares can be well explained by a composite "dome-plate" quasi-separatrix layer (QSL). Only half of these show a magnetic null point, with its fan and spine embedded in the dome and plate, respectively. The dome-plate QSL, therefore, is a general and robust structure characterizing circular-ribbon flares.

A new black hole X-ray binary (BHXRB) MAXI J1535-571 was discovered by MAXI during its outburst in 2017. Using observations taken by the first Chinese X-ray satellite, the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT), we perform a joint spectral analysis (2-150 keV) in both energy and time domains. The energy spectra provide the essential input for probing the intrinsic Quasi-Periodic Oscillation (QPO) fractional rms spectra (FRS). Our results show that during the intermediate state, the energy spectra are in general consistent with those reported by Swift/XRT and NuSTAR. However, the QPO FRS become harder and the FRS residuals may suggest the presence of either an additional power-law component in the energy spectrum or a turn-over in the intrinsic QPO FRS at high energies.

Weak gravitational lensing is one of the most promising cosmological probes to constrain dark matter, dark energy and the nature of gravity at cosmic scales. Intrinsic alignments (IA) of galaxies have been recognized as one of the most serious systematic effects facing gravitational lensing. Such alignments must be isolated and removed to obtain a pure lensing signal. Furthermore, the alignments are related to the processes of galaxy formation, so their extracted signal can help in understanding such formation processes and improving their theoretical modeling. We report in this letter the first detection of the gravitational shear--intrinsic shape (GI) correlation and the intrinsic shape--galaxy density (Ig) correlation in a photometric redshift survey using the self-calibration method. These direct measurements are made from the KiDS-450 photometric galaxy survey with a significance of 2.74$\sigma$ in the third bin for the Ig correlation, and 2.73$\sigma$ for the GI cross-correlation between the third and fourth bins. The self-calibration method uses the information available from photometric surveys without needing to specify an IA model and will play an important role in validating IA models and IA mitigation in future surveys such as LSST, Euclid and WFIRST.

The second observational campaign of gravitational waves organized by the LIGO/Virgo Collaborations has led to several breakthroughs such as the detection of gravitational wave signals from merger systems involving black holes or neutrons stars. During O2,14 gravitational wave alerts were sent to the astronomical community with sky regions covering mostly over hundreds of square degrees. Among them, 6 have been finally confirmed as real astrophysical events. Since 2013, a new set of ground-based robotic telescopes called GWAC and its pathfinder mini-GWAC have been developed to contribute to the various challenges of themulti-messenger and time domain astronomy. The GWAC system is built up in the framework of the ground-segment system of the SVOM mission that will be devoted to the study of the multi-wavelength transient sky in the next decade. During O2, only the mini-GWAC telescopenetwork was fully operational. Due to the wide field of view and fast automatic follow-up capabilities of the mini-GWAC telescopes, they were well adapted to efficiently cover the sky localization areas of the gravitational wave event candidates. In this paper, we present the mini-GWAC pipeline we have set up to respond to the GW alerts and we report our optical follow-up observations of 8 GW alerts detected during the O2 run. Our observations provided the largest coverage of the GW localization areas in a short latency made by any optical facility. We found tens of optical transient candidates in our images, but none of those could be securely associated with any confirmed black hole-black hole merger event. Based on this first experience and the near future technical improvements of our network system, we will be more competitive to detect the optical counterparts from some gravitational wave events that will be detected during the upcoming O3 run, especially those emerging from binary neutron star mergers.

The precision of the cosmological data allows us to accurately approximate the predictions for cosmological observables by Taylor expanding up to a low order the dependence on the cosmological parameters around a reference cosmology. By applying this observation to the redshift-space one-loop galaxy power spectrum of the Effective Field Theory of Large-Scale Structure, we analyze the BOSS DR12 data by scanning over all the parameters of $\Lambda$CDM cosmology with massive neutrinos. We impose several sets of priors, the widest of which is just a Big Bang Nucleosynthesis prior on the current fractional energy density of baryons, $\Omega_b h^2$, and a bound on the sum of neutrino masses to be less than 0.9 eV. In this case we measure the primordial amplitude of the power spectrum, $A_s$, the abundance of matter, $\Omega_m$, the Hubble parameter, $H_0$, and the tilt of the primordial power spectrum, $n_s$, to about $19\%$, $5.7\%$, $2.2\%$ and $7.3\%$ respectively, obtaining $\ln ( 10^{10} A_s) =2.91\pm 0.19$, $\Omega_m=0.314\pm 0.018$, $H_0=68.7\pm 1.5$ km/(s Mpc) and $n_s=0.979\pm 0.071$ at $68\%$ confidence level. A public code is released with this preprint.

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.

Solar type III radio bursts are excited by electron beams propagating outward from the Sun. The flux of type III radio burst has a time profile of rising and decay phase at a given frequency, which has been actively studied since 1970s. Several factors that may influence the duration of a type III radio burst has been proposed. However, the major cause of the duration is still an open question. In this work, to study the dominant cause of the duration, we investigate the source positions of the front edge, the peak, and the tail edge in the dynamic spectrum of a single and clear type III radio burst. The duration of this type III burst at a given frequency is about 3 second for decameter wave. The beam-formed observations by the LOw-Frequency ARray (LOFAR) are used, which can provide the radio source positions and the dynamic spectra at the same time. We find that, for this burst, the source positions of the front edge, the peak, and the tail edge split with each other spatially. The radial speed of the front edge, the peak, and the tail edge is 0.42 c, 0.25 c, and 0.16 c, respectively. We estimate the influences of the corona density fluctuation and the electron-velocity dispersion on the duration, and the scattering effect by comparison with a few short-duration bursts from the same region. The analysis yields that, in the frequency range of 30 - 41 MHz, the electron-velocity dispersion is the dominant factor that determines the time duration of type III radio bursts with long duration, while scattering may play important role in the duration of short bursts.

Cosmological tensions can arise within $\Lambda$CDM scenario amongst different observational windows, such as measurements of the $H_0$ and $\sigma_8$ parameters. These tensions, if finally confirmed by measurements, may indicate new physics beyond the standard paradigm. In this Letter, we report how to alleviate both the $H_0$ and $\sigma_8$ tensions simultaneously within torsional gravity from the perspective of effective field theory (EFT). We apply the EFT approach, which allows to investigate the evolution equations at the background and perturbation levels in a systematic way, and examine the conditions followed by the coefficients of various possibly involved operators such that cosmological tensions can be relaxed. Following these observations we construct concrete models of Lagrangians of torsional gravity. Specifically, we consider the parametrization $f(T)=-T -2\Lambda/M_P^2 +\alpha T^\beta$, where two out of the three parameters are independent (in which an additional term of the form $cT^{1/2}$ can be added). This model can efficiently fit observations solving all statistical tensions. To our knowledge, this is the first time where a modified gravity theory can alleviate both $H_0$ and $\sigma_8$ tensions simultaneously, hence offering an additional argument in favor of gravitational modification.

The Analytical Blind Separation (ABS) method is dedicated to extract the CMB signal from its contaminants without any prior information about foregrounds. In this study, we apply the ABS method to simulated polarization sky maps in order to test its ability in the recovery of the CMB $E$- and $B$-mode power spectra. We simulate multi-frequency $Q$ and $U$ microwave sky maps in the presence of two polarized foreground components (synchrotron and thermal dust), together with an uncorrelated Gaussian distributed instrumental noise, measured with a hypothetical PICO-like experiment operating in 10-frequency bands in the range of 30 to 321 GHz. Both full and partial sky surveys are taken into account in our analysis. For the full sky case with the tensor-to-scalar ratios of $r=0$ and $r=0.05$, the ABS can accurately recover the $E$-mode and $B$-mode power spectra with a relative error below $20\%$ with respect to the input CMB power spectra for the full multipole range ($2 \leq \ell \leq 1050$), close to theoretical expectations. For the partial sky analysis with a E/B separation proposed by \cite{ebmixture10} to avoid "$E$-to-$B$" leakage, the results are comparable to the full sky case. We find that the ABS is able to estimate both CMB $E$-mode and $B$-mode power spectra within 1-$\sigma$ at all scales for full sky, and for most of the scales ($50 \leq \ell \leq 1050$) for partial sky. For low-$\ell$ multipoles ($\ell \leq 50$), a noticeable divergence occurs in the partial sky analysis, which might be due to uncertainties originated from both the ABS and the E/B separation. A more detailed analysis for the partial sky case is thus needed for estimating primordial $B$-modes unbiasedly and accurately.

The Effective Field Theory of Large-Scale Structure (EFTofLSS) is a formalism that allows us to predict the clustering of Cosmological Large-Scale Structure in the mildly non-linear regime in an accurate and reliable way. After validating our technique against several sets of numerical simulations, we perform the analysis for the cosmological parameters of the DR12 BOSS data. We assume $\Lambda$CDM, a fixed value of the baryon/dark-matter ratio, $\Omega_b/\Omega_c$, and of the tilt of the primordial power spectrum, $n_s$, and no significant input from numerical simulations. By using the one-loop power spectrum multipoles, we measure the primordial amplitude of the power spectrum, $A_s$, the abundance of matter, $\Omega_m$, and the Hubble parameter, $H_0$, to about $13\%$, $3.2\%$ and $3.2\%$ respectively, obtaining $\ln(10^{10}As)=2.72\pm 0.13$, $\Omega_m=0.309\pm 0.010$, $H_0=68.5\pm 2.2$ km/(s Mpc) at 68\% confidence level. If we then add a CMB prior on the sound horizon, the error bar on $H_0$ is reduced to $1.6\%$. These results are a substantial qualitative and quantitative improvement with respect to former analyses, and suggest that the EFTofLSS is a powerful instrument to extract cosmological information from Large-Scale Structure.

Considering that the existence of relativistic particles in the protostellar jet has been confirmed by the detection of linearly polarized radio emission from the HH 80-81 jet, we search for $\gamma$-rays (which could be produced by the relativistic particles) from the HH 80-81 system using ten-year {\it Fermi}-LAT observations. A significant $\gamma$-ray excess is found in the direction of the HH 80-81 system with Test-Statistic (TS) value$>$100 , which is likely produced in the HH 80-81 jet. The $\gamma$-ray spectrum extends only to 1 GeV with a photon index of 3.5. No significant variability is found in the $\gamma$-ray emission. Combining the $\gamma$-ray and previous radio observations, we may exclude the relativistic Bremsstrahlung process for the production of the $\gamma$-rays under the assumption of a magnetic field strength of 0.1 mG. $\gamma$-rays can be produced through inverse-Compton (IC) scattering of relativistic electrons or inelastic proton-proton ($pp$) interaction. The maximum energies of electrons and protons are constrained to be $<0.1\ $TeV and $\sim\ $20 GeV, respectively.

The statistical properties of X-ray flares from two separate locations (nucleus and HST-1) in the M87 jet are investigated to reveal the physical origin of the flares. We analyse the archival \textit{Chandra} data for M87, and identify 14 flares in the nucleus and 9 flares in HST-1. The peak intensity ($I_{\rm{P}}$) and the flaring duration time ($T_{\rm{fl}}$) for each flare are obtained. It is found that the distributions of both $I_{\rm{P}}$ and $T_{\rm{fl}}$ for the nucleus obey a power-law form with a similar index. A similar result is also obtained for HST-1, and no significant inconsistency between the nucleus and HST-1 is found for the indices. Similar to solar X-ray flares, the power-law distributions of the flare event parameters can be well explained by a self-organized criticality (SOC) system, which are triggered by magnetic reconnection. Our results suggest that the flares from nucleus and HST-1 are possibly triggered by magnetic reconnection process. The consistent indices for the distributions of $I_{\rm{P}}$ and $T_{\rm{fl}}$ in the CORE and HST-1 indicate that the dimensions of the energy dissipation of the magnetic reconnection is identical in the two regions. A strong correlation between the flares in the two regions also suggests a similar physical origin for the flares.

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.

In this work, we propose a forward-modeling method to study the trajectory and speed of the interplanetary (IP) Type-III radio burst exciter. The model assumes that the source of an IP Type-III radio burst moves outward from the Sun following the Parker spiral field line. Using the arrival time of the radio waves at multiple spacecraft, we are able to determine the trajectory of the radio source in the Ecliptic plane, and its outward speed, as well as the injection time and longitude of the associated electron beam near the solar surface that triggers the Type-III radio burst. For the application of this method, we design a system to gather the arrival time of the radio wave from the radio dynamic spectra observed by Solar Terrestrial Relations Observatory (STEREO)/WAVES and Wind/WAVES. Then the system forward models the trajectory and speed of the radio burst exciter iteratively according to an evaluation function. Finally, we present a survey of four Type-III radio bursts that are well discussed in the literature. The modeled trajectories of the radio source are consistent with the previous radio-triangulation results, the longitude of the associated active region, or the location of Langmuir waves excited by the electron beam.

The Eisenhart lift of a Paul Trap used to store ions in molecular physics is a linearly polarized periodic gravitational wave. A modified version of Dehmelt's Penning Trap is in turn related to circularly polarized periodic gravitational waves, sought for in inflationary models. Similar equations rule also the Lagrange points in Celestial Mechanics. The explanation is provided by anisotropic oscillators.

Both multi-streaming (random motion) and bulk motion cause the Finger-of-God (FoG) effect in redshift space distortion (RSD). We apply a direct measurement of the multi-streaming effect in RSD from simulations, proving that it induces an additional, non-negligible FoG damping to the redshift space density power spectrum. We show that, including the multi-streaming effect, the RSD modelling is significantly improved. We also provide a theoretical explanation based on halo model for the measured effect, including a fitting formula with one to two free parameters. The improved understanding of FoG helps break the $f\sigma_8-\sigma_v$ degeneracy in RSD cosmology, and has the potential of significantly improving cosmological constraints.

Extracting CMB B-mode polarization from complicated foregrounds is a challenging task in searching for inflationary gravitational waves. We propose the ABS method as a blind and analytical solution to this problem. It applies to the measured cross bandpower between different frequency bands and obtains the CMB B-mode bandpower analytically. It does not rely on assumptions of foregrounds and does not require multiple parameter fitting. Testing against a variety of foregrounds, survey frequency configurations and instrument noise, we verify its applicability and numerical stability. The ABS method also applies to CMB temperature, E-mode polarization, the thermal Sunyaev Zel'dovich effect, spectral distortion, and even significantly different problems such as cosmic magnification.

Recently Chen et al. (2018, ApJ, 861, 58) accurately determined the volume weighted halo velocity bias in simulations, and found that the deviation of velocity bias from unity is much weaker than the peak model prediction. Here we present a possible explanation of this vanishing velocity bias. The starting point is that, halos are peaks in the low redshift {\it non-Gaussian} density field with smoothing scale $R_{\Delta}$ (virial radius), instead of peaks in the high redshift initial {\it Gaussian} density field with a factor of $\mathcal{O}(\Delta^{1/3})$ larger smoothing scale. Based on the approximation that the density field can be Gaussianized by a local and monotonic transformation, we extend the peak model to the non-Gaussian density field and derive the analytical expression of velocity dispersion and velocity power spectrum of these halos. The predicted deviation of velocity bias from unity is indeed much weaker than the previous prediction, and the agreement with the simulation results is significantly improved.

We point out a new possibility to determine the average redshift distribution of a large sample of gravitational wave standard sirens, without spectroscopic follow-ups. It is based on the cross correlation between the luminosity-distance space large scale structure (LSS) traced by standard sirens, and the redshift space LSS traced by galaxies in preexisting electromagnetic wave observations. We construct an unbiased and model independent estimator $E_z$ to realize this possibility. We demonstrate with BBO and Euclid that, $0.1\%$ accuracy in redshift determination can be achieved. This method can significantly alleviate the need of spectroscopic follow-up of standard sirens, and enhance their cosmological applications.

In this study, we apply the Analytical method of Blind Separation (ABS) of the cosmic microwave background (CMB) from foregrounds to estimate the CMB temperature power spectrum from multi-frequency microwave maps. We test the robustness of the ABS estimator and assess the accuracy of the power spectrum recovery by using realistic simulations based on the seven-frequency Planck data, including various frequency-dependent and spatially-varying foreground components (synchrotron, free-free, thermal dust and anomalous microwave emission), as well as an uncorrelated Gaussian-distributed instrumental noise. Considering no prior information about the foregrounds, the ABS estimator can analytically recover the CMB power spectrum over almost all scales with less than $0.5\%$ error for maps where the Galactic plane region ($|b|<10^{\circ}$) is masked out. To further test the flexibility and effectiveness of the ABS approach in a variety of situations, we apply the ABS to the simulated Planck maps in three cases: (1) without any mask, (2) imposing a two-times-stronger synchrotron emission and (3) including only the Galactic plane region ($|b|<10^{\circ}$) in the analysis. In such extreme cases, the ABS approach can still provide an unbiased estimate of band powers at the level of 1 $\mu\rm{K}^2$ on average over all $\ell$ range, and the recovered powers are consistent with the input values within 1-$\sigma$ for most $\ell$ bins.

Since 2016 October, the active galaxy PKS 2247-131 has undergone a gamma-ray outburst, which we studied using data obtained with the Fermi Gamma-ray Space Telescope. The emission arises from a relativistic jet in PKS 2247-131, as an optical spectrum only shows a few weak absorption lines, typical of the BL Lacertae sub-class of the blazar class of active galactic nuclei. Here we report a ~34.5 day quasi-periodic oscillation (QPO) in the emission after the initial flux peak of the outburst. Compared to one-year time-scale QPOs, previously identified in blazars in Fermi energies, PKS 2247-131 exhibits the first clear case of a relatively short, month-like oscillation. We show that this QPO can be explained in terms of a helical structure in the jet, where the viewing angle to the dominant emission region in the jet undergoes periodic changes. The time scale of the QPO suggests the presence of binary supermassive black holes in PKS 2247-131.

A relativistic electron-positron pair beam can be produced in the interaction of TeV photons from a blazar with the extragalactic background light (EBL). The relativistic $e^{\pm}$ pairs would loss energy through inverse-Compton scattering (ICS) photons of cosmic microwave background (CMB) or plasma instabilities. The dominant energy-loss process is under debate. Based on the assumption that the dominant energy-loss process is ICS, the resulted cascade GeV radiation is usually used to constrain the intergalactic magnetic field (IGMF). Here, we include the energy-loss due to plasma oblique instability in the calculation of cascade gamma-ray flux, and investigate the impact of the plasma instability on the constraint of IGMF. The up-to-date GeV data and archival TeV data of the blazar 1ES 0229+200 are used. The results indicate that even if the oblique instability cooling is dominating over ICS cooling, the cascade flux could be still used to constrain the IGMF. It is found that with the ratio between the cooling rates of the oblique instability and the ICS varying from 0.1, 1 to 10, the lower limit of the IGMF putted by the cascade flux and the gamma-ray data changes from $8\times10^{-18}\ $G, $5\times10^{-18}\ $G to $10^{-18}\ $G. If the ratio between the two cooling rates is 30, the estimate of IGMF based on the cascade flux is invalid.

Future gravitational wave (GW) observations are capable of detecting millions of compact star binary mergers in extragalactic galaxies, with $1\%$ luminosity-distance ($D_L$) measurement accuracy and better than arcminute positioning accuracy. This will open a new window of the large scale structure (LSS) of the universe, in the 3D {\bf luminosity-distance space (LDS)}, instead of the 3D redshift space of galaxy spectroscopic surveys. The baryon acoustic oscillation and the AP test encoded in the LDS LSS constrain the $D_L$-$D^{\rm co}_A$ (comoving angular diameter distance) relation and therefore the expansion history of the universe. Peculiar velocity induces the LDS distortion, analogous to the redshift space distortion, and allows for a new structure growth measure $f_L\sigma_8$. When the distance duality is enforced ($1+z=D_L/D^{\rm co}_A$), the LDS LSS by itself determines the redshift to $\sim 1\%$ level accuracy, and alleviates the need of spectroscopic follow-up of GW events.But a more valuable application is to test the distance duality to $1\%$ level accuracy, in combination with conventional BAO and supernovae measurements. This will put stringent constraints on modified gravity models in which the gravitational wave $D^{\rm GW}_L$ deviates from the electromagnetic wave $D^{EM}_L$. All these applications require no spectroscopic follow-ups.

The Globular Cluster 47 Tucanae was firstly detected in gamma-rays by the Large Area Telescope (LAT) onboard the \emph{Fermi} Gamma-ray Space Telescope, and the gamma-ray emission has been widely attributed to the millisecond pulsars. In this work, we analyze the Fermi-LAT pass 8 data ranging from 2008 August to 2017 May and report the detection of a modulation with a period of $18.416\pm0.008$ hours at a significance level of $\sim 4.8\sigma$. This is the first time to detect a significant modulation with a period much longer than that of millisecond pulsars in gamma-rays from Globular Clusters. The periodic modulation signal displays in the {\it Swift}-BAT data as well. The phase-folded Chandra X-ray light curve of a point source may have provided an additional clue. These facts indicate the detection of the first young pulsar binary in gamma-rays in Globular Clusters.

We design an event recognition-analysis system that can automatically detect solar type III radio burst and can mine information of the burst from the dynamic spectra observed by Nancay Decameter Array (NDA). We investigate the frequency drift rate of type III bursts and the speed of electron beams responsible for the generation of the bursts. Several computer vision methods are used in this automatic analysis system. The Hough transform is performed to recognize the line segment associated with type III bursts in the dynamic spectra. A modified active contour model is used to track the backbone of the burst and estimate the frequency drift rate at different frequency channels. We run this system on the NDA data from 2012 to 2017, and give a statistical survey of the event number distribution, the starting and stopping frequencies of bursts, the frequency dependence of the drift rate, and the exciter speed using three corona density models. The median value of the average frequency drift rates is about 6.94MHz/s for 1389 simple well-isolated type III bursts detected in the frequency range 10--80 MHz of NDA observation. The frequency drift rate changes with frequency as $df/dt = -0.0672 f^{1.23}$ from a least-squares fitting. The average exciter speed is about 0.2c based the density models. We do not find any significant dependence of the drift rate and the exciter speed on the solar activity cycle.

The mechanism for producing gamma-ray quasi-periodic oscillation (QPO) in blazar is unknown. One possibility is the geometric model, in which without the need for intrinsic quasi-periodic variation, the relativistic Doppler factor changes periodically, resulting in observed gamma-ray QPO. We propose a method to test this geometric model. We analyze the {\it Fermi}-LAT data of PG 1553+113 spanning from 2008 August until 2018 February. According to 29 four-month average spectral energy distributions (SEDs) in the energy range of 0.1-300 GeV, we split the {\it Fermi}-LAT energy range into three bands: 0.1-1 GeV, 1-10 GeV, and 10-300 GeV. The spectrum in each energy range can be successfully fitted by a power-law. The light curves and photon indices in the three energy ranges are obtained. Then, light curves in three narrow energy ranges, i.e., 0.2-0.5 GeV, 2-5 GeV and 20- 40 GeV, are constructed, and the relative variability amplitudes in the three narrow energy ranges are calculated. A discrete-correlation analysis is performed for the light curves. Our results indicate that (i) the light curves in the different energy ranges follow the same pattern showed in the light curve above 0.1 GeV; (ii) the three groups of photon indices in the energy ranges of 0.1-1 GeV, 1-10 GeV, and 10-300 GeV keep nearly constant; (iii) the ratio between relative variability amplitudes in different narrow energy ranges are equal (within their errors) to the prediction by the Doppler effect. Our results support the scenario of the relativistic boost producing the gamma-ray QPO for PG 1553+113.

The galaxy intrinsic alignment (IA) is a major challenge of weak lensing cosmology. To alleviate this problem, Zhang (2010, MNRAS, 406, L95) proposed a self-calibration method, independent of IA modeling. This proposal relies on several scaling relations between two-point clustering of IA and matter/galaxy fields, which were previously only tested with analytical IA models. In this paper, these relations are tested comprehensively with an $N$-body simulation of $3072^3$ simulation particles and boxsize 600 $h^{-1} \, \mathrm{Mpc}$. They are verified at the accuracy level of $\mathcal{O}(1)\%$ over angular scales and source redshifts of interest. We further confirm that these scaling relations are generic, insensitive to halo mass, weighting in defining halo ellipticities, photo-$z$ error, and misalignment between galaxy ellipticities and halo ellipticities. We also present and verify three new scaling relations on the B-mode IA. These results consolidate and complete the theory side of the proposed self-calibration technique.

The energy dissipation mechanism in blazar jet is unknown. Blazar's flares could provide insights on this problem. Here we report statistical results of XMM-Newton X-ray flares of Mrk 421. We analyze all public XMM-Newton X-ray observations for Mrk 421, and construct the light curves. Through fitting light curves, we obtain the parameters of flare-profiles, such as peak flux ($F_{\rm p}$) and flaring time duration ($T_{\rm fl}$). It is found that both the distributions of $F_{\rm p}$ and $T_{\rm fl}$ obey a power-law form, with the same index of $\alpha_{\rm F}=\alpha_{\rm T}\approx1$. The statistical properties are consistent with the predictions by a self-organized criticality (SOC) system with energy dissipation in one-dimensional space. This is similar to solar flare, but with different space dimensions of the energy dissipation domain. This suggests that X-ray flaers of Mrk 421 are possibly driven by a magnetic reconnection mechanism. Moreover, in the analysis, we find that variability on timescale of $\sim1000\ $s frequently appears. Such rapid variability indicates a magnetic field of $\geq 2.1\delta_{\rm D}^{-1/3}$ G ($\delta_{\rm D}$ is the Doppler factor) in emission region.

A long-standing issue in peculiar velocity cosmology is whether the halo/galaxy velocity bias $b_v=1$ at large scale. The resolution of this important issue must resort to high precision cosmological simulations. However, this is hampered by another long-standing `sampling artifact' problem in volume weighted velocity measurement. We circumvent this problem with a hybrid approach. We first measure statistics free of sampling artifact, then link them to volume weighted statistics in theory, finally solve for the velocity bias. $b_v$ determined by our method is not only free of sampling artifact, but also free of cosmic variance. We apply this method to a $\Lambda$CDM N-body simulation of $3072^3$ particles and $1200 Mpc/{\rm h}$ box size. For the first time, we determine the halo velocity bias to $0.1\%$-$1\%$ accuracy. Our major findings are as follows: (1) $b_v\neq 1$ at $k>0.1 h/{\rm Mpc}$. The deviation from unity ($|b_v-1|$) increases with $k$. Depending on halo mass and redshift, it may reach $\mathcal{O}(0.01)$ at $k=0.2 h/{\rm Mpc}$ and $\mathcal{O}(0.05)$ at $k\sim 0.3 h/{\rm Mpc}$. The discovered $b_v\neq 1$ has statistically significant impact on structure growth rate measurement by spectroscopic redshift surveys, including DESI, Euclid and SKA. (2) Both the sign and the amplitude of $b_v-1$ depend on mass and redshift. These results disagree with the peak model prediction in that $b_v$ has much weaker deviation from unity, varies with redshift, and can be bigger than unity. (3) Most of the mass and redshift dependences can be compressed into a single dependence on the halo density bias. Based on this finding, we provide an approximate two-parameter fitting formula.

In paper I of this series (Yang et al. 2017, ApJ), we proposed an analytical method of blind separation ({\bf ABS}) to extract the cosmic magnification signal in galaxy number distribution and reconstruct the weak lensing power spectrum. Here we report a new version of the ABS method, with significantly improved performance. This version is characterized by a shift parameter $\mathcal{S}$, with the special case of $\mathcal{S}=0$ corresponding to the original ABS method. We have tested this new version, compared to the previous one, and confirmed its supreme performance in all investigated situations. Therefore it supersedes the previous version. The proof of concept studies presented in this paper demonstrate that it may enable surveys such as LSST and SKA to reconstruct the lensing power spectrum at $z\simeq 1$ with $1\%$ accuracy. We will test with more realistic simulations to verify its applicability in real data.

POLAR is a compact wide-field space-borne detector for precise measurements of the linear polarisation of hard X-rays emitted by transient sources in the energy range from 50 keV to 500 keV. It consists of a 40$\times$40 array of plastic scintillator bars used as a detection material. The bars are grouped in 25 detector modules. The energy range sensitivity of POLAR is optimized to match with the prompt emission photons from the gamma-ray bursts (GRBs). Polarization measurements of the prompt emission would probe source geometries, emission mechanisms and magnetic structures in GRB jets. The instrument can also detect hard X-rays from solar flares and be used for precise measurement of their polarisation. POLAR was launched into a low Earth orbit on-board the Chinese space-lab TG-2 on September 15th, 2016. To achieve high accuracies in polarisation measurements it is essential to assure both before and during the flight a precise energy calibration. Such calibrations are performed with four low activity $^{22}$Na radioactive sources placed inside the instrument. Energy conversion factors are related to Compton edge positions from the collinear annihilation photons from the sources. This paper presents main principles of the in-flight calibration, describes studies of the method based on Monte Carlo simulations and its laboratory verification and finally provides some observation results based on the in-flight data analysis.

Circularly polarized gravitational sandwich waves exhibit, as do their linearly polarized counterparts, the Velocity Memory Effect: freely falling test particles in the flat after-zone fly apart along straight lines with constant velocity. In the inside zone their trajectories combine oscillatory and rotational motions in a complicated way. For circularly polarized periodic gravitational waves some trajectories remain bounded, while others spiral outward. These waves admit an additional "screw" isometry beyond the usual five. The consequences of this extra symmetry are explored.

Using isochronous mass spectrometry at the experimental storage ring CSRe in Lanzhou, the masses of $^{82}$Zr and $^{84}$Nb were measured for the first time with an uncertainty of $\sim 10$ keV, and the masses of $^{79}$Y, $^{81}$Zr, and $^{83}$Nb were re-determined with a higher precision. %The latter differ significantly from their literature values. The latter are significantly less bound than their literature values. Our new and accurate masses remove the irregularities of the mass surface in this region of the nuclear chart. Our results do not support the predicted island of pronounced low $\alpha$ separation energies for neutron-deficient Mo and Tc isotopes, making the formation of Zr-Nb cycle in the $rp$-process unlikely. The new proton separation energy of $^{83}$Nb was determined to be 490(400)~keV smaller than that in the Atomic Mass Evaluation 2012. This partly removes the overproduction of the $p$-nucleus $^{84}$Sr relative to the neutron-deficient molybdenum isotopes in the previous $\nu p$-process simulations.

We report the first full-dimensional potential energy surface (PES) and quantum mechanical close-coupling calculations for scattering of SiO due to H$_2$. The full-dimensional interaction potential surface was computed using the explicitly correlated coupled-cluster (CCSD(T)-F12b) method and fitted using an invariant polynomial approach. Pure rotational quenching cross sections from initial states $v_1=0$, $j_1$=1-5 of SiO in collision with H$_2$ are calculated for collision energies between 1.0 and 5000 cm$^{-1}$. State-to-state rotational rate coefficients are calculated at temperatures between 5 and 1000 K. The rotational rate coefficients of SiO with para-H$_2$ are compared with previous approximate results which were obtained using SiO-He PESs or scaled from SiO-He rate coefficients. Rovibrational state-to-state and total quenching cross sections and rate coefficients for initially excited SiO($v_1=1, j_1$=0 and 1) in collisions with para-H$_2$($v_2=0,j_2=0$) and ortho-H$_2$($v_2=0,j_2=1$) were also obtained. The application of the current collisional rate coefficients to astrophysics is briefly discussed.

In the narrow-line Seyfert 1 galaxy 1H 0707-495, recently a transient quasi-periodic oscillation (QPO) signal with a frequency of $\sim 2.6\times 10^{-4}$ Hz has been detected at a high statistical significance. Here, we reanalyze the same set of XMM-Newton data observed on 2008 February 4 with the Weighted-Wavelet Z-transform (WWZ) method. In addition to confirming the previous findings, we also find another QPO signal with a frequency of $\sim 1.2\times 10^{-4}$ Hz in a separated X-ray emission phase at the significance level of $\sim 3.7\sigma$. The signal is also found fitting an auto-regressive model though at a lower significance. The frequency ratio between these two signals is $\sim 2:1$. The analysis of other XMM-Newton measurements of 1H 0707-495 also reveals the presence of the $\sim 2.6\times 10^{-4}$ Hz ($\sim 1.2\times 10^{-4}$ Hz) QPO signal on 2007 May 14 (2010 September 17) at the significance level of $\sim 4.2\sigma$ ($\sim 3.5\sigma$). The QPO frequency found in this work follows the $f_{QPO}-M_{BH}$ relation reported in previous works spanning from stellar-mass to supermassive black holes. This is the first time to observe two separated transient X-ray QPO signals in active galactic nuclei (AGNs), which sheds new light on the physics of accreting supermassive black holes.

In the narrow-line Seyfert 1 galaxy 1H 0707-495, recently a transient quasi-periodic oscillation (QPO) signal with a frequency of $\sim 2.6\times 10^{-4}$ Hz has been detected at a high statistical significance. Here, we reanalyze the same set of XMM-Newton data observed on 2008 February 4 with the Weighted-Wavelet Z-transform (WWZ) method. In addition to confirming the previous findings, we also find another QPO signal with a frequency of $\sim 1.2\times 10^{-4}$ Hz in a separated X-ray emission phase at the significance level of $\sim 3.7\sigma$. The signal is also found fitting an auto-regressive model though at a lower significance. The frequency ratio between these two signals is $\sim 2:1$. The analysis of other XMM-Newton measurements of 1H 0707-495 also reveals the presence of the $\sim 2.6\times 10^{-4}$ Hz ($\sim 1.2\times 10^{-4}$ Hz) QPO signal on 2007 May 14 (2010 September 17) at the significance level of $\sim 4.2\sigma$ ($\sim 3.5\sigma$). The QPO frequency found in this work follows the $f_{QPO}-M_{BH}$ relation reported in previous works spanning from stellar-mass to supermassive black holes. This is the first time to observe two separated transient X-ray QPO signals in active galactic nuclei (AGNs), which sheds new light on the physics of accreting supermassive black holes.

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.

Photo-z error is one of the major sources of systematics degrading the accuracy of weak lensing cosmological inferences. Zhang et al. (2010) proposed a self-calibration method combining galaxy-galaxy correlations and galaxy-shear correlations between different photo-z bins. Fisher matrix analysis shows that it can determine the rate of photo-z outliers at a level of 0.01-1% merely using photometric data and do not rely on any prior knowledge. In this paper, we develop a new algorithm to implement this method by solving a constrained nonlinear optimization problem arising in the self-calibration process. Based on the techniques of fixed-point iteration and non-negative matrix factorization, the proposed algorithm can efficiently and robustly reconstruct the scattering probabilities between the true-z and photo-z bins. The algorithm has been tested extensively by applying it to mock data from simulated stage IV weak lensing projects. We find that the algorithm provides a successful recovery of the scatter rates at the level of 0.01-1%, and the true mean redshifts of photo-z bins at the level of 0.001, which may satisfy the requirements in future lensing surveys.

The "gravitational memory effect" due to an exact plane wave provides us with an elementary description of the diffeomorphisms associated with soft gravitons. It is explained how the presence of the latter may be detected by observing the motion of freely falling particles or other forms of gravitational wave detection. Numerical calculations confirm the relevance of the first, second and third time integrals of the Riemann tensor pointed out earlier. Solutions for various profiles are constructed. It is also shown how to extend our treatment to Einstein-Maxwell plane waves and a midi-superspace quantization is given.

In the narrow-line Seyfert 1 galaxy Mrk 766, a Quasi-Periodic Oscillation (QPO) signal with a period of $\sim6450$ s is detected in the \emph{XMM-Newton} data collected on 2005 May 31. This QPO signal is highly statistical significant at the confidence level at $\sim5\sigma$ with the quality factor of $Q=f/\Delta f>13.6$. The X-ray intensity changed by a factor of 3 with root mean square fractional variability of $14.3\%$. Furthermore, this QPO signal presents in the data of all three EPIC detectors and two RGS cameras and its frequency follows the $f_{\rm QPO}$-$M_{\rm BH}$ relation spanning from stellar-mass to supermassive black holes. Interestingly, a possible QPO signal with a period of $\sim4200$ s had been reported in the literature. The frequency ratio of these two QPO signals is $\sim$ 3:2. Our result is also in support of the hypothesis that the QPO signals can be just transient. The spectral analysis reveals that the contribution of the soft excess component below $\sim$ 1 keV is different between epochs with and without QPO, this property as well as the former frequency-ratio are well detected in X-ray BH binaries, which may have shed some lights on the physical origin of our event.

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.

In the narrow-line Seyfert 1 galaxy Mrk 766, a Quasi-Periodic Oscillation (QPO) signal with a period of $\sim6450$ s is detected in the \emph{XMM-Newton} data collected on 2005 May 31. This QPO signal is highly statistical significant at the confidence level at $\sim5\sigma$ with the quality factor of $Q=f/\Delta f>13.6$. The X-ray intensity changed by a factor of 3 with root mean square fractional variability of $14.3\%$. Furthermore, this QPO signal presents in the data of all three EPIC detectors and two RGS cameras and its frequency follows the $f_{\rm QPO}$-$M_{\rm BH}$ relation spanning from stellar-mass to supermassive black holes. Interestingly, a possible QPO signal with a period of $\sim4200$ s had been reported in the literature. The frequency ratio of these two QPO signals is $\sim$ 3:2. Our result is also in support of the hypothesis that the QPO signals can be just transient. The spectral analysis reveals that the contribution of the soft excess component below $\sim$ 1 keV is different between epochs with and without QPO, this property as well as the former frequency-ratio are well detected in X-ray BH binaries, which may have shed some lights on the physical origin of our event.

We give an account of the gravitational memory effect in the presence of the exact plane wave solution of Einstein's vacuum equations. This allows an elementary but exact description of the soft gravitons and how their presence may be detected by observing the motion of freely falling particles. The theorem of Bondi and Pirani on caustics (for which we present a new proof) implies that the asymptotic relative velocity is constant but not zero, in contradiction with the permanent displacement claimed by Zel'dovich and Polnarev. A non-vanishing asymptotic relative velocity might be used to detect gravitational waves through the "velocity memory effect", considered by Braginsky, Thorne, Grishchuk, and Polnarev.

We propose an Analytical method of Blind Separation (ABS) of cosmic magnification from the intrinsic fluctuations of galaxy number density in the observed galaxy number density distribution. The ABS method utilizes the different dependences of the signal (cosmic magnification) and contamination (galaxy intrinsic clustering) on galaxy flux, to separate the two. It works directly on the measured cross galaxy angular power spectra between different flux bins. It determines/reconstructs the lensing power spectrum analytically, without assumptions of galaxy intrinsic clustering and cosmology. It is unbiased in the limit of infinite number of galaxies. In reality the lensing reconstruction accuracy depends on survey configurations, galaxy biases, and other complexities, due to finite number of galaxies and the resulting shot noise fluctuations in the cross galaxy power spectra. We estimate its performance (systematic and statistical errors) in various cases. We find that, stage IV dark energy surveys such as SKA and LSST are capable of reconstructing the lensing power spectrum at $z\simeq 1$ and $\ell\la 5000$ accurately. This lensing reconstruction only requires counting galaxies, and is therefore highly complementary to the cosmic shear measurement by the same surveys.

We report a nominally high-confidence $\gamma$-ray quasi-periodic modulation in the blazar PKS 0301$-$243. For this target, we analyze its \emph{Fermi}-LAT Pass 8 data covering from 2008 August to 2017 May. Two techniques, i.e., the maximum likelihood optimization and the exposure-weighted aperture photometry, are used to build the $\gamma$-ray light curves. Then both the Lomb-Scargle Periodogram and the Weighted Wavelet Z-transform are applied to the light curves to search for period signals. A quasi-periodicity with a period of $2.1\pm0.3$ yr appears at the significance level of $\sim5\sigma$, although it should be noted that this putative quasi-period variability is seen in a data set barely four times longer. We speculate that this $\gamma$-ray quasi-periodic modulation might be evidence of a binary supermassive black hole.

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

We search for $\gamma$-ray and optical periodic modulations in a distant flat spectrum radio quasar (FSRQ) PKS 0426-380 (the redshift $z=1.1$). Using two techniques (i.e., the maximum likelihood optimization and the exposure-weighted aperture photometry), we obtain $\gamma$-ray light curves from \emph{Fermi}-LAT Pass 8 data covering from 2008 August to 2016 December. We then analyze the light curves with the Lomb-Scargle Periodogram (LSP) and the Weighted Wavelet Z-transform (WWZ). A $\gamma$-ray quasi-periodicity with a period of 3.35 $\pm$ 0.68 years is found at the significance-level of $\simeq3.6\ \sigma$. The optical-UV flux covering from 2005 August to 2013 April provided by ASI SCIENCE DATA CENTER is also analyzed, but no significant quasi-periodicity is found. It should be pointed out that the result of the optical-UV data could be tentative because of the incomplete of the data. Further long-term multiwavelength monitoring of this FSRQ is needed to confirm its quasi-periodicity.

Weak lensing statistics is typically measured as weighted sum of shear estimators or their products (shear-shear correlation). The weighting schemes are designed in the hope of minimizing the statistical error without introducing systematic errors. It would be ideal to approach the Cramer-Rao bound (the lower bound of the statistical uncertainty) in shear statistics, though it is generally difficult to do so in practice. The reasons may include: difficulties in galaxy shape measurement, inaccurate knowledge of the probability-distribution-function (PDF) of the shear estimator, misidentification of point sources as galaxies, etc.. Using the shear estimators defined in Zhang et al. (2015), we show that one can overcome all these problems, and allow shear measurement accuracy to approach the Cramer-Rao bound. This can be achieved by symmetrizing the PDF of the shear estimator, or the joint PDF of shear estimator pairs (for shear-shear correlation), without any prior knowledge of the PDF. Using simulated galaxy images, we demonstrate that under general observing conditions, this idea works as expected: it minimizes the statistical uncertainty without introducing systematic error.

As a space-borne detector POLAR is designed to conduct hard X-ray polarization measurements of gamma-ray bursts on the statistically significant sample of events and with an unprecedented accuracy. During its development phase a number of tests, calibrations runs and verification measurements were carried out in order to validate instrument functionality and optimize operational parameters. In this article we present results on gain optimization togeter with verification data obtained in the course of broad laboratory and environmental tests. In particular we focus on exposures to the $^{137}$Cs radioactive source and determination of the gain dependence on the high voltage for all 1600 detection channels of the polarimeter. Performance of the instrument is described in detail with respect to the dynamic range, energy resolution and temperature dependence. Gain optimization algorithms and response non-uniformity studies are also broadly discussed. Results presented below constitute important parts for development of the POLAR calibration and operation database.

Measuring the volume-weighted peculiar velocity statistics from inhomogeneously and sparsely distributed galaxies/halos, by existing velocity assignment methods, suffers from a significant sampling artifact. As an alternative, the Kriging interpolation based on Gaussian processes was introduced and evaluated [Y. Yu, J. Zhang, Y. Jing, and P. Zhang, Phys. Rev. D 92, 083527 (2015)]. Unfortunately, the most straightforward application of Kriging does not perform better than the existing methods in the literature. In this work, we investigate two physically motivated extensions. The first takes into account of the anisotropic velocity correlations. The second introduces the nugget effect, on account of multi-streaming of the velocity field. We find that the performance is indeed improved. For sparsely sampled data [$n_P\lesssim 6\times 10^{-3}(h^{-1} {\rm Mpc})^{-3}$] where the sampling artifact is the most severe, the improvement is significant and is two-fold: 1) The scale of reliable measurement of the velocity power spectrum is extended by a factor $\sim 1.6$, and 2) the dependence on the velocity correlation prior is weakened by a factor of $\sim 2$. We conclude that such extensions are desirable for accurate velocity assignment by Kriging.

We examine the gamma-ray quasi-periodic variability of PKS 2155-304 with the latest publicly available \emph{Fermi}-LAT Pass 8 data which covers the years from 2008 August to 2016 October. We produce the light curves in two ways, i.e., the exposure-weighted aperture photometry and the maximum likelihood optimization. The light curves are then analyzed by using Lomb-Scargle Periodogram (LSP) and Weighted Wavelet Z-transform (WWZ); and the results reveal a significant quasi-periodicity with a period of $1.74\pm0.13$ years and a significance of $\sim$ 4.9 $\sigma$. The constraint of multifrequencies quasi-periodic variabilities on blazar emission model is discussed.

To take full advantage of the unprecedented power of upcoming weak lensing surveys, understanding the noise, such as cosmic variance and geometry/mask effects, is as important as understanding the signal itself. Accurately quantifying the noise requires a large number of statistically independent mocks for a variety of cosmologies. This is impractical for weak lensing simulations, which are costly for simultaneous requirements of large box size (to cover a significant fraction of the past light cone) and high resolution (to robustly probe the small scale where most lensing signal resides). Therefore fast mock generation methods are desired and are under intensive investigation. We propose a new fast weak lensing map generation method, named the inverse-Gaussianization method, based on the finding that a lensing convergence field can be Gaussianized to excellent accuracy by a local transformation [Yu et al, Phys. Rev. D 84, 023523 (2011)]. Given a simulation, it enables us to produce as many as infinite statistically independent lensing maps as fast as producing the simulation initial conditions. The proposed method is tested against simulations for each tomography bin centered at lens redshift $z \sim 0.5$, 1, and 2, with various statistics. We find that the lensing maps generated by our method have reasonably accurate power spectra, bispectra, and power spectrum covariance matrix. Therefore, it will be useful for weak lensing surveys to generate realistic mocks. As an example of application, we measure the probability distribution function of the lensing power spectrum, from 16384 lensing maps produced by the inverse-Gaussianization method.

We test one-zone synchrotron self-Compton (SSC) models with high-quality multiwavelength spectral energy distribution (SED) data of Mrk 421. We use Markov chain Monte Carlo (MCMC) technique to fit twelve day-scale SEDs of Mrk 421 with one-zone SSC models. Three types of electron energy distribution (EED), a log-parabola (LP) EED, a power-law log-parabola (PLLP) EED and a broken power-law (BPL) EED, are assumed in fits. We find that the one-zone SSC model with the PLLP EED provides successful fits to all the twelve SEDs. However, the one-zone SSC model with the LP and BPL EEDs fail to provide acceptable fits to the highest energy X-ray data or GeV data in several states. We therefore conclude that the one-zone SSC model works well in explaining the SEDs of Mrk 421, and the PLLP EED is preferred over the LP and BPL EEDs for Mrk 421 during the flare in March 2010. We derive magnetic field $B'\sim0.01$\ G, Doppler factor $\delta_{\rm D}\sim$30--50, and the curvature parameter of EED $r\sim1$--$10$ in the model with the PLLP EED. The evolutions of model parameters are explored. The physical implications of our results are discussed.

We present a measurement of two-dimensional (2D) redshift-space power spectrum for the Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 11 CMASS galaxies in the North Galactic Cap (NGC) based on the method developed by Jing & Borner (2001). In this method, we first measure the 2D redshift-space correlation function for the CMASS galaxies, and obtain the 2D power spectrum based on Fourier Transform of the correlation function. The method is tested with an N-body mock galaxy catalog, which demonstrates that the method can yield an accurate and unbiased measurement of the redshift-space power spectrum given the input 2D correlation function is correct. Compared with previous measurements in literature that are usually based on direct Fourier Transform in redshift space, our method has the advantages that the window function and shot-noise are fully corrected. In fact, our 2D power spectrum, by its construction, can accurately reproduce the 2D correlation function, and in the meanwhile can reproduce, for example, the 2D power spectrum of Beutler et al. (2014) accurately if ours is convolved with the window function they provided. Thus, our measurement can facilitate a direct comparison with the theoretical predictions. With this accurate measurement of the 2D power spectrum, we then develop a method to measure the structure growth rate, by separating the anisotropic redshift-space power spectrum from the isotropic real-space power spectrum. We have also carefully corrected for the nonlinearities in the mapping from real space to redshift space, according to the theoretical model of Zhang et al. (2013). Finally, we obtain f(zeff)sigma_8(zeff)=0.438\pm0.037 at the effective redshift zeff=0.57, where f(zeff) is the linear growth rate at redshift zeff. The result is useful for constraining cosmological parameters. The measurements of 2D power spectrum will be released soon.

Power-law frequency distributions of the peak flux of solar flare X-ray emission have been studied extensively and attributed to a system of self-organized criticality (SOC). In this paper, we first show that, so long as the shape of the normalized light curve is not correlated with the peak flux, the flux histogram of solar flares also follows a power-law distribution with the same spectral index as the power-law frequency distribution of the peak flux, which may partially explain why power-law distributions are ubiquitous in the Universe. We then show that the spectral indexes of the histograms of soft X-ray fluxes observed by GOES satellites in two different energy channels are different: the higher energy channel has a harder distribution than the lower energy channel, which challenges the universal power-law distribution predicted by SOC models and implies a very soft distribution of thermal energy content of plasmas probed by the GOES. The temperature ($T$) distribution, on the other hand, approaches a power-law distribution with an index of 2 for high values of $T$. Application of SOC models to statistical properties of solar flares needs to be revisited.

The $Chandra$ archival data is a valuable resource for various studies on different topics of X-ray astronomy. In this paper, we utilize this wealth and present a uniformly processed data set, which can be used to address a wide range of scientific questions. The data analysis procedures are applied to 10,029 ACIS observations, which produces 363,530 source detections, belonging to 217,828 distinct X-ray sources. This number is twice the size of the $Chandra$ Source Catalog (Version 1.1). The catalogs in this paper provide abundant estimates of the detected X-ray source properties, including source positions, counts, colors, fluxes, luminosities, variability statistics, etc. Cross-correlation of these objects with galaxies shows 17,828 sources are located within the $D_{25}$ isophotes of 1110 galaxies, and 7504 sources are located between the $D_{25}$ and 2$D_{25}$ isophotes of 910 galaxies. Contamination analysis with the log$N$--log$S$ relation indicates that 51.3\% of objects within 2$D_{25}$ isophotes are truly relevant to galaxies, and the "net" source fraction increases to 58.9\%, 67.3\%, and 69.1\% for sources with luminosities above $10^{37}$, $10^{38}$, and $10^{39}$ erg s$^{-1}$. Among the possible scientific uses of this catalog, we discuss the possibility to study intra-observation variability, inter-observation variability, and supersoft sources.

In the third Fermi catalogue (3FGL) there are sixteen gamma-ray globular clusters. After analyzing the recent released Pass 8 data of Fermi Large Area Telescope (LAT), we report the discovery of significant gamma-ray emission from M 15 and NGC 6397, confirm that NGC 5904 is a gamma-ray emitter and find evidence of gamma-ray emission from NGC 6218 and NGC 6139. Inter- estingly, in the globular clusters M 15, NGC 6397 and NGC 5904, millisecond pulsars (MSPs) have been found in radio or X-rays, which are strongly in support of the MSP origin of the gamma-ray emission. However, due to the relative low luminosity of the gamma-ray emission we do not find any evidence for the gamma-ray pulsation or flux variability of these sources.

We present the spectroscopic and photometric observations on the Be/X-ray binary RX~J0440.9+4431 from 2001 to 2014. The short-term and long-term variability of the H$\alpha$ line profile indicates that one-armed global oscillations existed in the circumstellar disk. Several positive and negative correlations between the $V$-band brightness and the H$\alpha$ intensity \textbf{were} found from the long-term photometric and spectroscopic observations. We suggest that the monotonic increase of the $V$-band brightness and the H$\alpha$ brightness between our 2005-2007 observations might be the result of a continuous mass ejection from the central Be star, while the negative correlation in 2007-2010 should be caused by the cessation of mass loss from the Be star just before the decline in $V$-band brightness began (around our 2007 observations). With the extension of the ejection material, the largest circumstellar disk during the last two decades has been observed in our 2010 observations with an equivalent width of $\sim$$-12.88$ \AA, which corresponds to a circumstellar disk with a size of 12.9 times the radius of the central Be star. Three consecutive X-ray outbursts peaking around MJD 55293, 55444, and 55591 might be connected with the largest circumstellar disk around the Be star. We also use the orbital motion of the neutron star as a probe to constrain the structure of the circumstellar disk and estimate the eccentricity of the binary system to be $\ge$ 0.4. After three years of the H$\alpha$ intensity decline after the X-ray outbursts, a new circumstellar disk was being formed around the Be star after our 2013 observations.

[abridged] Volume-weighted statistics of large scale peculiar velocity is preferred by peculiar velocity cosmology, since it is free of uncertainties of galaxy density bias entangled in mass-weighted statistics. However, measuring the volume-weighted velocity statistics from galaxy (halo/simulation particle) velocity data is challenging. For the first time, we apply the Kriging interpolation to obtain the volume-weighted velocity field. Kriging is a minimum variance estimator. It predicts the most likely velocity for each place based on the velocity at other places. We test the performance of Kriging quantified by the E-mode velocity power spectrum from simulations. Dependences on the variogram prior used in Kriging, the number $n_k$ of the nearby particles to interpolate and the density $n_P$ of the observed sample are investigated. First, we find that Kriging induces $1\%$ and $3\%$ systematics at $k\sim 0.1h{\rm Mpc}^{-1}$ when $n_P\sim 6\times 10^{-2} ({\rm Mpc}/h)^{-3}$ and $n_P\sim 6\times 10^{-3} ({\rm Mpc}/h)^{-3}$, respectively. The deviation increases for decreasing $n_P$ and increasing $k$. When $n_P\lesssim 6\times 10^{-4} ({\rm Mpc}/h)^{-3}$, a smoothing effect dominates small scales, causing significant underestimation of the velocity power spectrum. Second, increasing $n_k$ helps to recover small scale power. However, for $n_P\lesssim 6\times 10^{-4} ({\rm Mpc}/h)^{-3}$ cases, the recovery is limited. Finally, Kriging is more sensitive to the variogram prior for lower sample density. The most straightforward application of Kriging on the cosmic velocity field does not show obvious advantages over the nearest-particle method (Zheng et al. 2013) and could not be directly applied to cosmology so far. However, whether potential improvements may be achieved by more delicate versions of Kriging is worth further investigation.

The tensor-vector-scalar (TeVeS) model is considered a viable theory of gravity. It produces the Milgrom's modified Newtonian dynamics in the nonrelativistic weak field limit and is free from ghosts. This model has been tested against various cosmological observations. Here we investigate whether new observations such as the galaxy velocity power spectrum measured by 6dF and the kinetic Sunyaev Zel'dovich effect power spectrum measured by ACT/SPT can put further constraints on the TeVeS model. Furthermore, we perform the test of TeVeS cosmology with a sterile neutrino by confronting to Planck data, and find that it is ruled out by cosmic microwave background measurements from the Planck mission.

We systematically analyze the 6-year {\it Fermi}/LAT data of the lobe-dominated quasars (LDQs) in the complete LDQ sample from 3CRR survey and report the discovery of high-energy $\gamma$-ray emission from 3C 275.1. The $\gamma$-ray emission of 3C 207 is confirmed and significant variability of the lightcurve is identified. We do not find statistically significant $\gamma$-ray emission from other LDQs. 3C 275.1 is the known $\gamma$-ray quasar with the lowest core dominance parameter (i.e., $R=0.11$). We also show that both the northern radio hotspot and parsec jet models can reasonably reproduce the $\gamma$-ray data. The parsec jet model, however, is favored by the potential $\gamma$-ray variability at the timescale of months. We suggest that some dimmer $\gamma$-ray LDQs will be detected in the future and LDQs could contribute non-negligibly to the extragalactic $\gamma$-ray background.

A profound assumption in peculiar velocity cosmology is $b_v=1$ at sufficiently large scales, where $b_v$ is the volume weighted halo(galaxy) velocity bias with respect to the matter velocity field. However, this fundamental assumption has not been robustly verified in numerical simulations. Furthermore, it is challenged by structure formation theory (BBKS, 1986, ApJ; Desjacques and Sheth, 2010, PRD), which predicts the existence of velocity bias (at least for proto-halos) due to the fact that halos reside in special regions (local density peaks). The major obstacle to measure the volume weighted velocity from N-body simulations is an unphysical sampling artifact. It is entangled in the measured velocity statistics and becomes significant for sparse populations. With recently improved understanding of the sampling artifact (Zhang, Zheng and Jing, 2015, PRD; Zheng, Zhang and Jing, 2015, PRD), for the first time we are able to {\it appropriately correct this sampling artifact and then robustly measure the volume weighted halo velocity bias}. (1) We verify $b_v=1$ within $2\%$ model uncertainty at $k\lesssim 0.1h/$Mpc and $z=0$-$2$ for halos of mass $\sim 10^{12}$-$10^{13} h^{-1} M_\odot$, and, therefore, consolidates a foundation of the peculiar velocity cosmology. (2) We also find statistically significant signs of $b_v\neq 1$ at $k\gtrsim 0.1h/$Mpc. Unfortunately, whether this is real or caused by residual sampling artifact requires further investigation. Nevertheless, cosmology based on $k\gtrsim 0.1h/$Mpc velocity data shall be careful this potential velocity bias.

Perhaps the most controversial idea in modern cosmology is that our observable universe is contained within one bubble among many, all inhabiting the eternally inflating multiverse. One of the few way to test this idea is to look for evidence of the relic inhomogeneities left by the collisions between other bubbles and our own. Such relic inhomogeneities induces a coherent bulk flow over gigaparsec scales. Therefore, bubble collisions leave unique imprints in the cosmic microwave background (CMB) through the kinetic Sunyaev Zel'dovich (kSZ) effect, temperature anisotropies induced by the scattering of photons from coherently moving free electrons in the diffuse intergalactic medium. The kSZ signature produced by bubble collisions has a unique directional dependence and is tightly correlated with the galaxy distribution; it can therefore be distinguished from other contributions to the CMB anisotropies. An important advantage of the kSZ signature is that it peaks on arcminute angular scales, where the limiting factors in making a detection are instrumental noise and foreground subtraction. This is in contrast to the collision signature in the primary CMB, which peaks on angular scales much larger than one degree, and whose detection is therefore limited by cosmic variance. In this paper, we examine the prospects for probing the inhomogeneities left by bubble collisions using the kSZ effect. We provide a forecast for detection using cross-correlations between CMB and galaxy surveys, finding that the detectability using the kSZ effect can be competitive with constraints from CMB temperature and polarization data.

Cosmic shear can only be measured where there are galaxies. This source-lens clustering (SLC) effect has two sources, intrinsic source clustering and cosmic magnification (magnification/size bias). Lensing tomography can suppress the former. However, this reduction is limited by the existence of photo-z error and nonzero redshift bin width. Furthermore, SLC induced by cosmic magnification cannot be reduced by lensing tomography. Through N-body simulations, we quantify the impact of SLC on the lensing power spectrum in the context of lensing tomography. We consider both the standard estimator and the pixel-based estimator. We find that none of them can satisfactorily handle both sources of SLC. (1) For the standard estimator, SLC induced by both sources can bias the lensing power spectrum by O(1)-O(10)%. Intrinsic source clustering also increases statistical uncertainties in the measured lensing power spectrum. However, the standard estimator suppresses intrinsic source clustering in the cross-spectrum. (2) In contrast, the pixel-based estimator suppresses SLC through cosmic magnification. However, it fails to suppress SLC through intrinsic source clustering and the measured lensing power spectrum can be biased low by O(1)-O(10)%. In short, for typical photo-z errors (sigma_z/(1+z)=0.05) and photo-z bin sizes (Delta_z^P=0.2), SLC alters the lensing E-mode power spectrum by 1-10%, with ell~10^3$ and z_s~1 being of particular interest to weak lensing cosmology. Therefore the SLC is a severe systematic for cosmology in Stage-IV lensing surveys. We present useful scaling relations to self-calibrate the SLC effect.

Weak lensing causes spatially coherent fluctuations in flux of type Ia supernovae (SNe Ia). This lensing magnification allows for weak lensing measurement independent of cosmic shear. It is free of shape measurement errors associated with cosmic shear and can therefore be used to diagnose and calibrate multiplicative error. Although this lensing magnification is difficult to measure accurately in auto correlation, its cross correlation with cosmic shear and galaxy distribution in overlapping area can be measured to significantly higher accuracy. Therefore these cross correlations can put useful constraint on multiplicative error, and the obtained constraint is free of cosmic variance in weak lensing field. We present two methods implementing this idea and estimate their performances. We find that, with $\sim 1$ million SNe Ia that can be achieved by the proposed D2k survey with the LSST telescope (Zhan et al. 2008), multiplicative error of $\sim 0.5\%$ for source galaxies at $z_s\sim 1$ can be detected and larger multiplicative error can be corrected to the level of $0.5\%$. It is therefore a promising approach to control the multiplicative to the sub-percent level required for stage IV projects. The combination of the two methods even has the potential to diagnose and calibrate galaxy intrinsic alignment, which is another major systematic error in cosmic shear cosmology.

Cosmology based on large scale peculiar velocity preferes volume weighted velocity statistics. However, measuring the volume weighted velocity statistics from inhomogeneously distributed galaxies (simulation particles/halos) suffer from an inevitable and significant sampling artifact. We study this sampling artifact in the velocity power spectrum measured by the nearest-particle (NP) velocity assignment method(Zheng et al. 2013, PRD). We derive the analytical expression of leading and higher order terms. We find that the sampling artifact suppresses the $z=0$ E-mode velocity power spectrum by $\sim 10\%$ at $k=0.1h/$Mpc , for samples with number density $10^{-3}({\rm Mpc}/h)^{-3}$. This suppression becomes larger for larger $k$ and for sparser samples. We argue that, this source of systematic errors in peculiar velocity cosmology, albeit severe, can be self-calibrated in the framework of our theoretical modelling. We also work out the sampling artifact in the density-velocity cross power spectrum measurement. More robust evaluation of related statistics through simulations will be presented in a companion paper (Zheng, Zhang & Jing, 2015, PRD). We also argue that similar sampling artifact exists in other velocity assignment methods and hence must be carefully corrected to avoid systematic bias in peculiar velocity cosmology.

Measuring the volume weighted velocity power spectrum suffers from a severe systematic error, due to imperfect sampling of the velocity field from inhomogeneous distribution of dark matter particles/halos in simulations or galaxies with velocity measurement. This "sampling artifact" depends on both the mean particle number density $\bar{n}_P$ and the intrinsic large scale structure (LSS) fluctuation in the particle distribution. (1) We report robust detection of this sampling artifact in N-body simulations. It causes $\sim 12$% underestimation of the velocity power spectrum at $k=0.1$h/Mpc for samples with $\bar{n}_P=6\times10^{-3}$ (Mpc/h)$^{-3}$. This systematic underestimation increases with decreasing $\bar{n}_P$ and increasing $k$. Its dependence on the intrinsic LSS fluctuations is also robustly detected. (2) All these findings are expected by our theoretical modelling in paper I \cite{Zhang14}. In particular, the leading order theoretical approximation agrees quantitatively well with simulation result for $\bar{n}_P\gtrsim6\times 10^{-4}$(Mpc/h)$^{-3}$. Furthermore, we provide an ansatz to take high order terms into account. It improves the model accuracy to $\lesssim1$% at $k\lesssim0.1$h/Mpc over 3 orders of magnitude in $\bar{n}_P$ and over typical LSS clustering from $z=0$ to $z=2$. (3) The sampling artifact is determined by the deflection ${\bf D}$ field, which is straightforwardly available in both simulations and data of galaxy velocity. Hence the sampling artifact in the velocity power spectrum measurement can be self-calibrated within our framework. By applying such self-calibration in simulations, it becomes promising to determine the {\it real} large scale velocity bias of $10^{13}M_\odot$ halos with $\sim 1$% accuracy, and that of lower mass halos by better accuracy. ...[abridged]

Accurate rate coefficients for molecular vibrational transitions due to collisions with H$_2$, critical for interpreting infrared astronomical observations, are lacking for most molecules. Quantum calculations are the primary source of such data, but reliable values that consider all internal degrees of freedom of the collision complex have only been reported for H$_2$-H$_2$ due to the difficulty of the computations. Here we present essentially exact full-dimensional dynamics computations for rovibrational quenching of CO due to H$_2$ impact. Using a high-level six-dimensional potential surface, time-independent scattering calculations, within a full angular-momentum-coupling formulation, were performed for the deexcitation of vibrationally excited CO. Agreement with experimentally-determined results confirms the accuracy of the potential and scattering computations, representing the largest of such calculations performed to date. This investigation advances computational quantum dynamics studies representing initial steps toward obtaining CO-H$_2$ rovibrational quenching data needed for astrophysical modeling.

Globular clusters (GCs) are the oldest stellar system in the Galaxy, in which the millisecond pulsars are widely believed to be the only steady {\gamma}-ray emitters. So far 9 {\gamma}-ray GCs have been identified and a few candidates such as 2MS-GC01 and IC 1257 have been suggested. In this work, after analyzing the publicly-available Fermi-LAT data we confirm the significant {\gamma}-ray emission from 2MS-GC01 and IC 1257 and report the discovery for {\gamma}-ray emission from NGC 5904 and NGC 6656 within their tidal radii. Also a strong evidence of significant {\gamma}-ray emission is found from FSR 1735. From the observed {\gamma}-ray luminosities, the numbers of MSPs that are expected to be present in these GCs are estimated.

The extended TeV gamma-ray source ARGO J2031+4157 (or MGRO J2031+41) is positionally consistent with the Cygnus Cocoon discovered by $Fermi$-LAT at GeV energies in the Cygnus superbubble. Reanalyzing the ARGO-YBJ data collected from November 2007 to January 2013, the angular extension and energy spectrum of ARGO J2031+4157 are evaluated. After subtracting the contribution of the overlapping TeV sources, the ARGO-YBJ excess map is fitted with a two-dimensional Gaussian function in a square region of $10^{\circ}\times 10^{\circ}$, finding a source extension $\sigma_{ext}$= 1$^{\circ}$.8$\pm$0$^{\circ}$.5. The observed differential energy spectrum is $dN/dE =(2.5\pm0.4) \times 10^{-11}(E/1 TeV)^{-2.6\pm0.3}$ photons cm$^{-2}$ s$^{-1}$ TeV$^{-1}$, in the energy range 0.2-10 TeV. The angular extension is consistent with that of the Cygnus Cocoon as measured by $Fermi$-LAT, and the spectrum also shows a good connection with the one measured in the 1-100 GeV energy range. These features suggest to identify ARGO J2031+4157 as the counterpart of the Cygnus Cocoon at TeV energies. The Cygnus Cocoon, located in the star-forming region of Cygnus X, is interpreted as a cocoon of freshly accelerated cosmic rays related to the Cygnus superbubble. The spectral similarity with Supernova Remnants indicates that the particle acceleration inside a superbubble is similar to that in a SNR. The spectral measurements from 1 GeV to 10 TeV allows for the first time to determine the possible spectrum slope of the underlying particle distribution. A hadronic model is adopted to explain the spectral energy distribution.

In this work, we aim to show the possibilities of the variable speed of light (VSL) theory in explaining the type Ia supernovae observations without introducing dark energy. The speed of light is assumed to be scale factor dependent, which is the most popular assumption in VSL theory. We show the modified calculation of the distance modulus, and the validity of the redshift-scale factor relation in VSL theory. Three different models of VSL are tested SNe data-sets with proper constraints on the model parameters. The comparison of the three models and flat $\Lambda$CDM in distance modulus is showed. Some basic problems and the difficulties of the confirmation of the VSL theory are also discussed

The ARGO-YBJ detector is an extensive air shower array that has been used to monitor the northern $\gamma$-ray sky at energies above 0.3 TeV from 2007 November to 2013 January. In this paper, we present the results of a sky survey in the declination band from $-10^{\circ}$ to $70^{\circ}$, using data recorded over the past five years. With an integrated sensitivity ranging from 0.24 to $\sim$1 Crab units depending on the declination, six sources have been detected with a statistical significance greater than 5 standard deviations. Several excesses are also reported as potential $\gamma$-ray emitters. The features of each source are presented and discussed. Additionally, $95\%$ confidence level upper limits of the flux from the investigated sky region are shown. Specific upper limits for 663 GeV $\gamma$-ray AGNs inside the ARGO-YBJ field of view are reported. The effect of the absorption of $\gamma$-rays due to the interaction with extragalactic background light is estimated.

The interaction between Dark Matter and Dark Energy has been proposed as a mechanism to alleviate the coincidence problem. We analyze the effect of the interaction on the evolution of the gravitational field and propose two new observables based on its effect on the matter peculiar velocity field. We find that for different model parameters the matter peculiar velocity could be factor 2 times larger or 5 times smaller than the amplitude of velocity perturbations in the concordance LCDM cosmological model at the same scales. We compare the effect on the peculiar velocities with those on the Integrated Sachs-Wolfe effect and we show that velocities can potentially provide constraints on the strength of the interaction stronger than those currently available. We show that the current upper limits on the amplitude of the kinetic Sunyaev-Zel'dovich power spectrum provide constraints on the coupling within the dark sectors that are consistent with those obtained previously from the CMB and galaxy clusters. In particular, we show that current upper limits from the Atacama Cosmology Telescope and the South Pole Telescope favor Dark Energy decaying into Dark Matter, as required to solve the coincidence problem.

The existence of galaxy intrinsic clustering severely hampers the weak lensing reconstruction from cosmic magnification. In paper I \citep{Yang2011}, we proposed a minimal variance estimator to overcome this problem. By utilizing the different dependences of cosmic magnification and galaxy intrinsic clustering on galaxy flux, we demonstrated that the otherwise overwhelming galaxy intrinsic clustering can be significantly suppressed such that lensing maps can be reconstructed with promising accuracy. This procedure relies heavily on the accuracy of determining the galaxy bias from the same data. Paper I adopts an iterative approach, which degrades toward high redshift. The current paper presents an alternative method, improving over paper I. We prove that the measured galaxy clustering between flux bins allows for simultaneous determination of the lensing power spectrum and the flux dependence of galaxy bias, at this redshift bin. Comparing to paper I, the new approach is not only more straightforward, but also more robust. It identifies an ambiguity in determining the galaxy bias and further discovers a mathematically robust way to suppress this ambiguity to non-negligible level ($\sim 0.1%$). The accurately determined galaxy bias can then be applied to the minimal variance estimator proposed in paper I to improve the lensing map-making. The gain at high redshift is significant. These maps can be used to measure other statistics, such as cluster finding and peak statistics. Furthermore, by including galaxy clustering measurement between different redshift bins, we can also determine the lensing cross power spectrum between these bins, up to a small and correctable multiplicative factor.

Wide range equation of state (EOS) for liquid hydrogen is ultimately built by combining two kinds of density functional theory (DFT) molecular dynamics simulations, namely, first-principles molecular dynamics simulations and orbital-free molecular dynamics simulations. Specially, the present introduction of short cutoff radius pseudopotentials enables the hydrogen EOS to be available in the range $9.82\times10^{-4}$ to $1.347\times10^{3}$ g/cm$^{3}$ and up to $5\times10^{7}$ K. By comprehensively comparing with various attainable experimental and theoretical data, we derive the conclusion that our DFT-EOS can be readily and reliably conducted to hydrodynamic simulations of the inertial confinement fusion.

The dynamics of collisionless galaxy can be described by the Vlasov-Poisson system. By the Jean's theorem, all the spherically symmetric steady galaxy models are given by a distribution of {\Phi}(E,L), where E is the particle energy and L the angular momentum. In a celebrated Doremus-Feix-Baumann Theorem, the galaxy model {\Phi}(E,L) is stable if the distribution {\Phi} is monotonically decreasing with respect to the particle energy E. On the other hand, the stability of {\Phi}(E,L) remains largely open otherwise. Based on a recent abstract instability criterion of Guo-Lin, we constuct examples of unstable galaxy models of f(E,L) and f(E) in which f fails to be monotone in E.

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

We report the observation of a very high energy \gamma-ray source, whose position is coincident with HESS J1841-055. This source has been observed for 4.5 years by the ARGO-YBJ experiment from November 2007 to July 2012. Its emission is detected with a statistical significance of 5.3 standard deviations. Parameterizing the source shape with a two-dimensional Gaussian function we estimate an extension \sigma=(0.40(+0.32,-0.22}) degree, consistent with the HESS measurement. The observed energy spectrum is dN/dE =(9.0-+1.6) x 10^{-13}(E/5 TeV)^{-2.32-+0.23} photons cm^{-2} s^{-1} TeV^{-1}, in the energy range 0.9-50 TeV. The integral \gamma-ray flux above 1 TeV is 1.3-+0.4 Crab units, which is 3.2-+1.0 times the flux derived by HESS. The differences in the flux determination between HESS and ARGO-YBJ, and possible counterparts at other wavelengths are discussed.

The statistics of dark matter halos is an essential component of understanding the nonlinear evolution in modified gravity cosmology. Based on a series of modified gravity N-body simulations, we investigate the halo mass function, concentration and bias. We model the impact of modified gravity by a single parameter \zeta, which determines the enhancement of particle acceleration with respect to GR, given the identical mass distribution (\zeta=1 in GR). We select snapshot redshifts such that the linear matter power spectra of different gravity models are identical, in order to isolate the impact of gravity beyond modifying the linear growth rate. At the baseline redshift corresponding to z_S=1.2 in the standard \Lambda CDM, for a 10% deviation from GR(|\zeta-1|=0.1), the measured halo mass function can differ by about 5-10%, the halo concentration by about 10-20%, while the halo bias differs significantly less. These results demonstrate that the halo mass function and/or the halo concentration are sensitive to the nature of gravity and may be used to make interesting constraints along this line.

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

We search for the dark flow induced diffuse kinetic Sunyaev Zel'dovich (kSZ) effect through CMB-galaxy cross correlation. Such angular correlation is anisotropic, with a unique $cos(theta_DF)$ angular dependence and hence can be distinguished from other components. Here, $theta_DF$ is the angle between the opposite dark flow direction and the direction of the sky where the correlation is measured. We analyze the KIAS-VAGC galaxy catalog of SDSS-DR7 and the WMAP seven-year temperature maps, applying an unbiased optimal weighting scheme to eliminate any statistically isotropic components and to enhance the dark flow detection signal. Non-zero weighted cross correlations are detected at 3.5 sigma for the redshift bin z<0.1 and at 3 sigma for the bin 0.1<z<0.2, implying the existence of statistically anisotropic components in CMB. However, further analysis does not support the dark flow explanation. The observed directional dependence deviates from the $\propto cos(theta_DF)$ relation expected, and hence can not be explained by the presence of a single dark flow, and if the observed cross correlation is generated by the dark flow induced kSZ effect, the velocity would be too high (> 6000 km/s). We report this work as the first attempt to search for dark flow through weighted CMB-galaxy cross correlation and to draw the attention on the sources of the detected anomalous CMB-galaxy cross correlation.

As one of the brightest active blazars in both X-ray and very high energy $\gamma$-ray bands, Mrk 501 is very useful for physics associated with jets from AGNs. The ARGO-YBJ experiment is monitoring it for $\gamma$-rays above 0.3 TeV since November 2007. Starting from October 2011 the largest flare since 2005 is observed, which lasts to about April 2012. In this paper, a detailed analysis is reported. During the brightest $\gamma$-rays flaring episodes from October 17 to November 22, 2011, an excess of the event rate over 6 $\sigma$ is detected by ARGO-YBJ in the direction of Mrk 501, corresponding to an increase of the $\gamma$-ray flux above 1 TeV by a factor of 6.6$\pm$2.2 from its steady emission. In particular, the $\gamma$-ray flux above 8 TeV is detected with a significance better than 4 $\sigma$. Based on time-dependent synchrotron self-Compton (SSC) processes, the broad-band energy spectrum is interpreted as the emission from an electron energy distribution parameterized with a single power-law function with an exponential cutoff at its high energy end. The average spectral energy distribution for the steady emission is well described by this simple one-zone SSC model. However, the detection of $\gamma$-rays above 8 TeV during the flare challenges this model due to the hardness of the spectra. Correlations between X-rays and $\gamma$-rays are also investigated.

In paper I (Yu et al. [1]), we show through N-body simulation that a local monotonic Gaussian transformation can significantly reduce non-Gaussianity in a noise-free lensing convergence field. This makes the Gaussianization a promising theoretical tool to understand high-order lensing statistics. Here we present a study of its applicability in lensing data analysis, in particular when shape measurement noise is presented in lensing convergence maps. (i) We find that shape measurement noise significantly degrades the Gaussianization performance and the degradation increases for shallower surveys. (ii) The Wiener filter is efficient in reducing the impact of shape measurement noise. The Gaussianization of the Wiener-filtered lensing maps is able to suppress skewness, kurtosis, and the 5th- and 6th-order cumulants by a factor of 10 or more. It also works efficiently to reduce the bispectrum to zero.

Hill's equations, which first arose in the study of the Earth-Moon-Sun system, admit the two-parameter centrally extended Newton-Hooke symmetry without rotations. This symmetry allows for extending Kohn's theorem about the center-of-mass decomposition. Particular light is shed on the problem using Duval's "Bargmann" framework. The separation of the center-of-mass motion into that of a guiding center and relative motion is derived by a generalized chiral decomposition.

We present a concept study on weak lensing map reconstruction through the cosmic magnification effect in galaxy number density distribution. We propose a minimal variance linear estimator to minimize both the dominant systematical and statistical errors in the map reconstruction. It utilizes the distinctively different flux dependences to separate the cosmic magnification signal from the overwhelming galaxy intrinsic clustering noise. It also minimizes the shot noise error by an optimal weighting scheme on the galaxy number density in each flux bin. Our method is in principle applicable to all galaxy surveys with reasonable redshift information. We demonstrate its applicability against the planned Square Kilometer Array survey, under simplified conditions. Weak lensing maps reconstructed through our method are complementary to that from cosmic shear and CMB and 21cm lensing. They are useful for cross checking over systematical errors in weak lensing reconstruction and for improving cosmological constraints.

Consistent descriptions of the equation of states, and information about transport coefficients of deuterium-tritium mixture are demonstrated through quantum molecular dynamic (QMD) simulations (up to a density of 600 g/cm$^{3}$ and a temperature of $10^{4}$ eV). Diffusion coefficients and viscosity are compared with one component plasma model in different regimes from the strong coupled to the kinetic one. Electronic and radiative transport coefficients, which are compared with models currently used in hydrodynamic simulations of inertial confinement fusion, are evaluated up to 800 eV. The Lorentz number is also discussed from the highly degenerate to the intermediate region.

We consider the generic scenario of dark energy which arises through the latent heat of a hidden sector first order cosmological phase transition. This field could account for the extra radiation degree of freedom suggested by the CMB. We present the bubble nucleation solution for the viscous limit. The decay rate of the field is constrained by published KSZ data, and may be an explanation of current excess ISW correlations. Cross correlation of current and future surveys can further constrain or test the parameter space. The decay model is plausibly in the observable range, and avoids anthropic problems. This class of models is not well constrained by the popular dark energy figure of merit.

We report the observation of TeV gamma-rays from the Cygnus region using the ARGO-YBJ data collected from 2007 November to 2011 August. Several TeV sources are located in this region including the two bright extended MGRO J2019+37 and MGRO J2031+41. According to the Milagro data set, at 20 TeV MGRO J2019+37 is the most significant source apart from the Crab Nebula. No signal from MGRO J2019+37 is detected by the ARGO-YBJ experiment, and the derived flux upper limits at 90% confidence level for all the events above 600 GeV with medium energy of 3 TeV are lower than the Milagro flux, implying that the source might be variable and hard to be identified as a pulsar wind nebula. The only statistically significant (6.4 standard deviations) gamma-ray signal is found from MGRO J2031+41, with a flux consistent with the measurement by Milagro.

We present a detailed theoretical analysis of a non-thermal escape of molecular hydrogen from Mars induced by collisions with hot atomic oxygen from martian corona. To accurately describe the energy transfer in O + H$_2(v,j)$ collisions, we performed extensive quantum-mechanical calculations of state-to-state elastic, inelastic, and reactive cross sections. The escape flux of H$_2$ molecules was evaluated using a simplified 1D column model of the martian atmosphere with realistic densities of atmospheric gases and hot oxygen production rates for the low solar activity conditions. An average density of the non-thermal escape flux of H$_2$ of $1.9\times10^5$ cm$^{-2}$s$^{-1}$ was obtained considering energetic O atoms produced in dissociative recombinations of O$_{2}^{+}$ ions. Predicted rovibrational distribution of the escaping H$_2$ was found to contain a significant fraction of higher rotational states. While the non-thermal escape rate was found to be lower than Jeans flux for H$_2$ molecules, the non-thermal escape rates of HD and D$_2$ are significantly higher than their respective Jeans rates. The accurate values of non-thermal escape fluxes of different molecular isotopes of H$_2$ may be important in analyses of evolution of the martian atmosphere. The described molecular ejection mechanism is general and expected to contribute to atmospheric escape of H$_2$ and other light molecules from planets, satellites, and exoplanetary bodies.

Motivated by recent works of Neyrinck et al. 2009 and Scherrer et al. 2010, we proposed a Gaussianization transform to Gaussianize the non-Gaussian lensing convergence field $\kappa$. It performs a local monotonic transformation $\kappa\rightarrow y$ pixel by pixel to make the unsmoothed one-point probability distribution function of the new variable $y$ Gaussian. We tested whether the whole $y$ field is Gaussian against N-body simulations. (1) We found that the proposed Gaussianization suppresses the non-Gaussianity by orders of magnitude, in measures of the skewness, the kurtosis, the 5th- and 6th-order cumulants of the $y$ field smoothed over various angular scales relative to that of the corresponding smoothed $\kappa$ field. The residual non-Gaussianities are often consistent with zero within the statistical errors. (2) The Gaussianization significantly suppresses the bispectrum. Furthermore, the residual scatters around zero, depending on the configuration in the Fourier space. (3) The Gaussianization works with even better performance for the 2D fields of the matter density projected over $\sim 300 \mpch$ distance interval centered at $z\in(0,2)$, which can be reconstructed from the weak lensing tomography. (4) We identified imperfectness and complexities of the proposed Gaussianization. We noticed weak residual non-Gaussianity in the $y$ field. We verified the widely used logarithmic transformation as a good approximation to the Gaussianization transformation. However, we also found noticeable deviations.

The Copernican principle, a cornerstone of modern cosmology, remains largely unproven at Gpc radial scale and above. Here we will show that, violations of this type will inevitably cause a first order anisotropic kinetic Sunyaev Zel'dovich (kSZ) effect. If large scale radial inhomogeneities have amplitude large enough to explain the "dark energy" phenomena, the induced kSZ power spectrum will be much larger than the ACT/SPT upper limit. This single test confirms the Copernican principle and rules out the adiabatic void model as a viable alternative to dark energy.

BigBOSS is a Stage IV ground-based dark energy experiment to study baryon acoustic oscillations (BAO) and the growth of structure with a wide-area galaxy and quasar redshift survey over 14,000 square degrees. It has been conditionally accepted by NOAO in response to a call for major new instrumentation and a high-impact science program for the 4-m Mayall telescope at Kitt Peak. The BigBOSS instrument is a robotically-actuated, fiber-fed spectrograph capable of taking 5000 simultaneous spectra over a wavelength range from 340 nm to 1060 nm, with a resolution R = 3000-4800. Using data from imaging surveys that are already underway, spectroscopic targets are selected that trace the underlying dark matter distribution. In particular, targets include luminous red galaxies (LRGs) up to z = 1.0, extending the BOSS LRG survey in both redshift and survey area. To probe the universe out to even higher redshift, BigBOSS will target bright [OII] emission line galaxies (ELGs) up to z = 1.7. In total, 20 million galaxy redshifts are obtained to measure the BAO feature, trace the matter power spectrum at smaller scales, and detect redshift space distortions. BigBOSS will provide additional constraints on early dark energy and on the curvature of the universe by measuring the Ly-alpha forest in the spectra of over 600,000 2.2 < z < 3.5 quasars. BigBOSS galaxy BAO measurements combined with an analysis of the broadband power, including the Ly-alpha forest in BigBOSS quasar spectra, achieves a FOM of 395 with Planck plus Stage III priors. This FOM is based on conservative assumptions for the analysis of broad band power (kmax = 0.15), and could grow to over 600 if current work allows us to push the analysis to higher wave numbers (kmax = 0.3). BigBOSS will also place constraints on theories of modified gravity and inflation, and will measure the sum of neutrino masses to 0.024 eV accuracy.

ARGO-YBJ is an air shower detector array with a fully covered layer of resistive plate chambers. It is operated with a high duty cycle and a large field of view. It continuously monitors the northern sky at energies above 0.3 TeV. In this paper, we report a long-term monitoring of Mrk 421 over the period from 2007 November to 2010 February. This source was observed by the satellite-borne experiments Rossi X-ray Timing Explorer and Swift in the X-ray band. Mrk 421 was especially active in the first half of 2008. Many flares are observed in both X-ray and gamma-ray bands simultaneously. The gamma-ray flux observed by ARGO-YBJ has a clear correlation with the X-ray flux. No lag between the X-ray and gamma-ray photons longer than 1 day is found. The evolution of the spectral energy distribution is investigated by measuring spectral indices at four different flux levels. Hardening of the spectra is observed in both X-ray and gamma-ray bands. The gamma-ray flux increases quadratically with the simultaneously measured X-ray flux. All these observational results strongly favor the synchrotron self-Compton process as the underlying radiative mechanism.

Cosmic magnification has been detected through cross correlation between foreground and background populations (galaxies or quasars). It has been shown that weighing each background object by its $\alpha-1$ can significantly improve the cosmic magnification measurement \citep{Menard02,Scranton05}. Here, $\alpha$ is the logarithmic slope of the luminosity function of background populations. However, we find that this weighting function is optimal only for sparse background populations in which intrinsic clustering is negligible with respect to shot noise. We derive the optimal weighting function for general case including scale independent and scale dependent weights. The optimal weighting function improves the S/N (signal to noise ratio) by $\sim 20%$ for a BigBOSS-like survey and the improvement can reach a factor of $\sim 2$ for surveys with much denser background populations.

We study the thermophysical properties of dense helium plasmas by using quantum molecular dynamics and orbital-free molecular dynamics simulations, where densities are considered from 400 to 800 g/cm$^{3}$ and temperatures up to 800 eV. Results are presented for the equation of state. From the Kubo-Greenwood formula, we derive the electrical conductivity and electronic thermal conductivity. In particular, with the increase in temperature, we discuss the change in the Lorenz number, which indicates a transition from strong coupling and degenerate state to moderate coupling and partial degeneracy regime for dense helium.

We use 1-dimensional numerical simulations to study spherical collapse in the f(R) gravity models. We include the nonlinear coupling of the gravitational potential to the scalar field in the theory and use a relaxation scheme to follow the collapse. We find an unusual enhancement in density near the virial radius which may provide observable tests of gravity. We also use the estimated collapse time to calculate the critical overdensity $\delta_c$ used in calculating the mass function and bias of halos. We find that analytical approximations previously used in the literature do not capture the complexity of nonlinear spherical collapse.

We derive the exact analytical solution of the linear structure growth rate in LCDM cosmology with flat or curved geometry, under the Newtonian gauge. Unlike the well known solution under the Newtonian limit (Heath 1977), our solution takes all general relativistic corrections into account and is hence valid at both the sub- and super-horizon scales. With this exact solution, we evaluate cosmological impacts induced by these relativistic corrections. (1) General relativistic corrections alter the density growth from z=100 to z=0 by 10% at k=0.01h/Mpc and the impact becomes stronger toward larger scales. We caution the readers that the overdensity is not gauge invariant and the above statement is restrained to the Newtonian gauge. (2) Relativistic corrections introduce a k^{-2} scale dependence in the density fluctuation. It mimics a primordial non-Gaussianity of the local type with $f^{\rm local}_{\rm NL}\sim 1$. This systematical error may become non-negligible for future all sky deep galaxy surveys. (3) Cosmological simulations with box size greater than 1Gpc are also affected by these relativistic corrections. We provide a post-processing recipe to correct for these effects. (4) These relativistic corrections affect the redshift distortion. However, at redshifts and scales relevant to redshift distortion measurements, such effect is negligible.

The thermal Sunyaev-Zel'dovich (tSZ) effect directly measures the thermal pressure of free electrons integrated along the line of sight and thus contains valuable information on the thermal history of the universe. However, the redshift information is entangled in the projection along the line of sight. This projection effect severely degrades the power of the tSZ effect to reconstruct the thermal history. We investigate the tSZ tomography technique to recover this otherwise lost redshift information by cross correlating the tSZ effect with galaxies of known redshifts, or alternatively with matter distribution reconstructed from weak lensing tomography. We investigate in detail the 3D distribution of the gas thermal pressure and its relation with the matter distribution, through our adiabatic hydrodynamic simulation and the one with additional gastrophysics including radiative cooling, star formation and supernova feedback. (1) We find a strong correlation between the gas pressure and matter distribution, with a typical cross correlation coefficient r ~ 0.7 at k . 3h/Mpc and z < 2. This tight correlation will enable robust cross correlation measurement between SZ surveys such as Planck, ACT and SPT and lensing surveys such as DES and LSST, at ~20-100{\sigma} level. (2) We propose a tomography technique to convert the measured cross correlation into the contribution from gas in each redshift bin to the tSZ power spectrum. Uncertainties in gastrophysics may affect the reconstruction at ~ 2% level, due to the ~ 1% impact of gastrophysics on r, found in our simulations. However, we find that the same gastrophysics affects the tSZ power spectrum at ~ 40% level, so it is robust to infer the gastrophysics from the reconstructed redshift resolved contribution.

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

The sun blocks cosmic ray particles from outside the solar system, forming a detectable shadow in the sky map of cosmic rays detected by the ARGO-YBJ experiment in Tibet. Because the cosmic ray particles are positive charged, the magnetic field between the sun and the earth deflects them from straight trajectories and results in a shift of the shadow from the true location of the sun. Here we show that the shift measures the intensity of the field which is transported by the solar wind from the sun to the earth.

Weak lensing surveys exploit measurements of galaxy ellipticities. These measurements are subject to errors which degrade the cosmological information that can be extracted from the surveys. Here we propose a way of using the galaxy data themselves to calibrate the measurement errors. In particular, the cosmic shear field, which causes the galaxies to appear elliptical, also changes their sizes and fluxes. Information about the sizes and fluxes of the galaxies can be added to the shape information to obtain more robust information about the cosmic shear field. The net result will be tighter constraints on cosmological parameters such as those which describe dark energy.

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

The galaxy intrinsic alignment is a severe challenge to precision cosmic shear measurement. We propose to self-calibrate the induced gravitational shear-galaxy intrinsic ellipticity correlation (the GI correlation, \citealt{Hirata04b}) in weak lensing surveys with photometric redshift measurement. (1) We propose a method to extract the intrinsic ellipticity-galaxy density cross correlation (I-g) from the galaxy ellipticity-density measurement in the same redshift bin. (2) We also find a generic scaling relation to convert the extracted I-g correlation to the demanded GI correlation. We perform concept study under simplified conditions and demonstrate its capability to significantly reduce the GI contamination. We discuss the impact of various complexities on the two key ingredients of the self-calibration technique, namely the method to extract the I-g correlation and the scaling relation between the I-g and the GI correlation. We expect none of them is likely able to completely invalidate the proposed self-calibration technique.

The galaxy intrinsic alignment causes the galaxy ellipticity-ellipticity power spectrum between two photometric redshifts to decrease faster with respect to the redshift separation $\Delta z^P$, for fixed mean redshift. This offers a valuable diagnosis on the intrinsic alignment. We show that the distinctive dependences of the GG, II and GI correlations on $\Delta z^P$ over the range $|\Delta z^P|\la 0.2$ can be understood robustly without strong assumptions on the intrinsic alignment. This allows us to measure the intrinsic alignment within each conventional photo-z bin of typical size $\ga 0.2$, through lensing tomography of photo-z bin size $\sim 0.01$. Both the statistical and systematical errors in the lensing cosmology can be reduced by this self-calibration technique.

Recently Kashlinsky et al. 2008, 2010 reported a discovery of a $\sim 10^3$ km/$s$ bulk flow of the universe out to $z\simeq 0.3$, through the dark flow induced CMB dipole in directions of clusters. We point out that, if this dark flow exists, it will also induce observable CMB temperature fluctuations at multipole $\ell\sim 10^3$, through modulation of the inhomogeneous electron distribution on the uniform dark flow. The induced small scale kinetic Sunyaev Zel'dovich (SZ) effect will reach $\sim 1\muk^2$ at multipole $10^3\la \ell\la 10^4$, only a factor of $\sim 2$ smaller than the conventional kinetic SZ effect. Furthermore, it will be correlated with the large scale structure (LSS) and its correlation with 2MASS galaxy distribution reaches $0.3 \mu$K at $\ell=10^3$, under a directional dependent optimal weighting scheme. We estimate that, WMAP plus 2MASS should already be able to detect this dark flow induced small scale kinetic SZ effect with $\sim 6\sigma$ confidence. Deeper galaxy surveys such as SDSS can further improve the measurement. Planck plus existing galaxy surveys can reach $\ga 14\sigma$ detection. Existing CMB-LSS cross correlation measurements shall be reanalyzed to test the existence of the dark flow and, if it exists, shall be used to eliminate possible bias on the integrated Sachs-Wolfe effect measurement through the CMB-LSS cross correlation.

Correlated imaging through atmospheric turbulence is studied, and the analytical expressions describing turbulence effects on image resolution are derived. Compared with direct imaging, correlated imaging can reduce the influence of turbulence to a certain extent and reconstruct high-resolution images. The result is backed up by numerical simulations, in which turbulence-induced phase perturbations are simulated by random phase screens inserting propagation paths.

A large fraction of cosmological information on dark energy and gravity is encoded in the nonlinear regime. Precision cosmology thus requires precision modeling of nonlinearities in general dark energy and modified gravity models. We modify the Gadget-2 code and run a series of N-body simulations on modified gravity cosmology to study the nonlinearities. The modified gravity model that we investigate in the present paper is characterized by a single parameter \zeta, which determines the enhancement of particle acceleration with respect to general relativity (GR), given the identical mass distribution (\zeta = 1 in GR). The first nonlinear statistics we investigate is the nonlinear matter power spectrum at k < 3h/Mpc, which is the relevant range for robust weak lensing power spectrum modeling at l < 2000. In this study, we focus on the relative difference in the nonlinear power spectra at corresponding redshifts where different gravity models have the same linear power spectra. This particular statistics highlights the imprint of modified gravity in the nonlinear regime and the importance to include the nonlinear regime in testing GR. By design, it is less susceptible to the sample variance and numerical artifacts. We adopt a mass assignment method based on wavelet to improve the power spectrum measurement. We run a series of tests to determine the suitable simulation specifications (particle number, box size and initial redshift). We find that, the nonlinear power spectra can differ by ~30% for 10% deviation from GR (|\zeta-1| = 0.1) where the rms density fluctuations reach 10. This large difference, on one hand, shows the richness of information on gravity in the corresponding scales, and on the other hand, invalidates simple extrapolations of some existing fitting formulae to modified gravity cosmology.

We present the X-ray variability properties of the X-ray and TeV bright blazar Mrk 421 with a ~60 ks long XMM-Newton observation performed on November 9-10, 2005. The source experienced a pronounced flare, of which the inter-band time lags were determined with a very high confidence level. The soft (0.6-0.8 keV) X-ray variations lagged the hard (4-10 keV) ones by 1.09 ks, and the soft lag increases with increasing difference in the photon energy. The energy-dependent soft lags can be well fitted with the difference of the energy-dependent cooling timescales of the relativistic electron distribution responsible for the observed X-ray emission, which constrains the magnetic field strength and Doppler factor of the emitting region to be B\delta^(1/3)~1.78 Gauss.

Photo-z errors, especially catastrophic errors, are a major uncertainty for precision weak lensing cosmology. We find that the shear-(galaxy number) density and density-density cross correlation measurements between photo-z bins, available from the same lensing surveys, contain valuable information for self-calibration of the scattering probabilities between the true-z and photo-z bins. The self-calibration technique we propose does not rely on cosmological priors nor parameterization of the photo-z probability distribution function, and preserves all of the cosmological information available from shear-shear measurement. We estimate the calibration accuracy through the Fisher matrix formalism. We find that, for advanced lensing surveys such as the planned stage IV surveys, the rate of photo-z outliers can be determined with statistical uncertainties of 0.01-1% for $z<2$ galaxies. Among the several sources of calibration error that we identify and investigate, the {\it galaxy distribution bias} is likely the most dominant systematic error, whereby photo-z outliers have different redshift distributions and/or bias than non-outliers from the same bin. This bias affects all photo-z calibration techniques based on correlation measurements. Galaxy bias variations of $O(0.1)$ produce biases in photo-z outlier rates similar to the statistical errors of our method, so this galaxy distribution bias may bias the reconstructed scatters at several-$\sigma$ level, but is unlikely to completely invalidate the self-calibration technique.

Dark energy interacting with dark matter is a promising model to solve the cosmic coincidence problem. We study the signature of such interaction on large scale cosmic microwave background (CMB) temperature anisotropies. Based on the detail analysis in perturbation equations of dark energy and dark matter when they are in interaction, we find that the large scale CMB, especially the late Integrated Sachs Wolfe effect, is a useful tool to measure the coupling between dark sectors. We also discuss the possibility to detect the coupling by cross-correlating CMB maps with tracers of the large scale structure. We finally perform the global fitting to constrain the coupling by using the CMB power spectrum data together with other observational data. We find that in the $1\sigma$ range, the constrained coupling between dark sectors can solve the coincidence problem.

We reconstruct the interaction rate between the dark matter and the holographic dark energy with the parameterized equation of states and the future event horizon as the infrared cut-off length. It is shown that the observational constraints from the 192 SNIa and BAO measurement permit the negative interaction in the wide region. Moreover, the usual phenomenological descriptions can not describe the reconstructed interaction well for many cases. The other possible interaction is also discussed.

The Large Area Multi-Object Spectroscopic Telescope (LAMOST) is a dedicated spectroscopic survey telescope being built in China, with an effective aperture of 4 meters and equiped with 4000 fibers. Using the LAMOST telescope, one could make redshift survey of the large scale structure (LSS). The baryon acoustic oscillation (BAO) features in the LSS power spectrum provide standard rulers for measuring dark energy and other cosmological parameters. In this paper we investigate the meaurement precision achievable for a few possible surveys: (1) a magnitude limited survey of all galaxies, (2) a survey of color selected red luminous galaxies (LRG), and (3) a magnitude limited, high density survey of z<2 quasars. For each survey, we use the halo model to estimate the bias of the sample, and calculate the effective volume. We then use the Fisher matrix method to forecast the error on the dark energy equation of state and other cosmological parameters for different survey parameters. In a few cases we also use the Markov Chain Monte Carlo (MCMC) method to make the same forecast as a comparison. The fiber time required for each of these surveys is also estimated. These results would be useful in designing the surveys for LAMOST.

A generic prediction of the standard cosmology, based on general relativity (GR), dark matter and the cosmological constant (and more generally, smooth dark energy), is that, the two gravitational potentials describing the spatial and temporal scalar perturbations of the universe are equivalent. Modifications in GR or dark energy clustering in general violate this relation. Thus this ratio serves as a smoking gun of the dark universe. We propose a method to extract this ratio at various cosmological scales and redshifts from a set of measurements, in a model independent way. The ratio measured by future surveys has strong discriminating power for a variety of dark universe scenarios.

The application of large scale peculiar velocity (LSPV), as a crucial probe of dark matter, dark energy and gravity, is severely limited by measurement obstacles. We show that fluctuations in type Ia supernovae (SNe Ia) fluxes induced by LSPV offer a promising approach to measure LSPV at intermediate redshifts. In the 3D Fourier space, gravitational lensing, the dominant systematical error, is well suppressed, localized and can be further corrected effectively. Advance in SN observations can further significantly reduce shot noise induced by SN intrinsic fluctuations, which is the dominant statistical error. Robust mapping on the motion of the dark universe through SNe Ia is thus feasible to $z\sim 0.5$.

We propose to use spatial correlations of the kinetic Sunyaev-Zeldovich (KSZ) flux as an estimator of the peculiar velocity power spectrum. In contrast with conventional techniques, our new method does not require measurements of the thermal SZ signal or the X-ray temperature. Moreover, this method has the special advantage that the expected systematic errors are always sub-dominant to statistical errors on all scales and redshifts of interest. We show that future large sky coverage KSZ surveys may allow a peculiar velocity power spectrum estimates of an accuracy reaching ~10%.

Free electrons deplete photons from type Ia supernovae through the (inverse) Compton scattering. This Compton dimming increases with redshift and reaches 0.004 mag at $z=1$ and 0.01 mag at $z=2$. Although far from sufficient to invalidate the existence of dark energy, it can bias constraint on dark energy at a level non-negligible for future supernova surveys. This effect is correctable and should be incorporated in supernova analysis. The Compton dimming has similar impact on cosmology based on gamma ray bursts as standard candles.

Baryon acoustic oscillation (BAO) is a powerful probe on the expansion of the universe, shedding light on elusive dark energy and gravity at cosmological scales. BAO measurements through biased tracers of the underlying matter density field, as most proposals do, can reach high statistical accuracy. However, possible scale dependence in bias may induce non-negligible systematical errors, especially for the most ambitious spectroscopic surveys proposed. We show that precision spectroscopic redshift information available in these surveys allows for {\it self calibration} of the galaxy bias and its stochasticity, as function of scale and redshift. Through the effect of redshift distortion, one can simultaneously measure the real space power spectra of galaxies, galaxy-velocity and velocity, respectively. At relevant scales of BAO, galaxy velocity faithfully traces that of the underlying matter. This valuable feature enables a rather model independent way to measure the galaxy bias and its stochasticity by comparing the three power spectra. For the square kilometer array (SKA), this self calibration is statistically accurate to correct for 1% level shift in BAO peak positions induced by bias scale dependence. Furthermore, we find that SKA is able to detect BAO in the velocity power spectrum, opening a new window for BAO cosmology.

The standard cosmological model is based on general relativity and includes dark matter and dark energy. An important prediction of this model is a fixed relationship between the gravitational potentials responsible for gravitational lensing and the matter overdensity. Alternative theories of gravity often make different predictions for this relationship. We propose a set of measurements which can test the lensing/matter relationship, thereby distinguishing between dark energy/matter models and models in which gravity differs from general relativity. Planned optical, infrared and radio galaxy and lensing surveys will be able to measure $E_G$, an observational quantity whose expectation value is equal to the ratio of the Laplacian of the Newtonian potentials to the peculiar velocity divergence, to percent accuracy. We show that this will easily separate alternatives such as $\Lambda$CDM, DGP, TeVeS and $f(R)$ gravity.

We consider the Randall-Sundrum brane-world model with bulk-brane energy transfer where the Einstein-Hilbert action is modified by curvature correction terms: a four-dimensional scalar curvature from induced gravity on the brane, and a five-dimensional Gauss-Bonnet curvature term. It is remarkable that these curvature terms will not change the dynamics of the brane universe at low energy. Parameterizing the energy transfer and taking the dark radiation term into account, we find that the phantom divide of the equation of state of effective dark energy could be crossed, without the need of any new dark energy components. Fitting the two most reliable and robust SNIa datasets, the 182 Gold dataset and the Supernova Legacy Survey (SNLS), our model indeed has a small tendency of phantom divide crossing for the Gold dataset, but not for the SNLS dataset. Furthermore, combining the recent detection of the SDSS baryon acoustic oscillations peak (BAO) with lower matter density parameter prior, we find that the SNLS dataset also mildly favors phantom divide crossing.

Modified gravity theories have richer observational consequences for large-scale structure than conventional dark energy models, in that different observables are not described by a single growth factor even in the linear regime. We examine the relationships between perturbations in the metric potentials, density and velocity fields, and discuss strategies for measuring them using gravitational lensing, galaxy cluster abundances, galaxy clustering/dynamics and the ISW effect. We show how a broad class of gravity theories can be tested by combining these probes. A robust way to interpret observations is by constraining two key functions: the ratio of the two metric potentials, and the ratio of the Gravitational ``constant'' in the Poisson equation to Newton's constant. We also discuss quasilinear effects that carry signatures of gravity, such as through induced three-point correlations. Clustering of dark energy can mimic features of modified gravity theories and thus confuse the search for distinct signatures of such theories. It can produce pressure perturbations and anisotropic stresses, which breaks the equality between the two metric potentials even in general relativity. With these two extra degrees of freedom, can a clustered dark energy model mimic modified gravity models in all observational tests? We show with specific examples that observational constraints on both the metric potentials and density perturbations can in principle distinguish modifications of gravity from dark energy models. We compare our result with other recent studies that have slightly different assumptions (and apparently contradictory conclusions).

A recent modeling study of brightness ratios for CO rotational transitions in gas typical of the diffuse ISM by Liszt found the role of H collisions to be more important than previously assumed. This conclusion was based on quantum scattering calculations using the so-called WKS potential energy surface (PES) which reported a large cross section for the important 0->1 rotational transition. New close-coupling (CC) rigid-rotor calculations for CO(v=0,J=0) excitation by H are performed on four different PESs. Two of the PESs are obtained in this work using state-of-the-art quantum chemistry techniques at the CCSD(T) and MRCI levels of theory. Cross sections for the J=0->1, as well as other odd Delta J, transitions are significantly suppressed compared to even Delta J transitions in thermal energy CC calculations using the CCSD(T) and MRCI surfaces. This is consistent with a expected even Delta J propensity and in contrast to CC calculations using the earlier WKS PES which predict a dominating 0->1 transition. The current results suggest that the original astrophysical assumption that excitation of CO by H_2 dominates the kinetics of CO in diffuse ISM gas is likely to remain valid.

Lyman-$\alpha$ spheres, i.e. regions around the first stars which are illuminated by Lyman-$\alpha$ photons and show 21cm absorption feature against the CMB, are smoking guns at the dawn of the reionization epoch. Though overwhelming radio foreground makes their detections extremely difficult, we pointed out that, strong gravitational lensing can significantly improve their observational feasibility. Since Lyman-$\alpha$ spheres have ~10" sizes, comparable to the caustic size of galaxy clusters, individual images of each strongly lensed Lyman-$\alpha$ sphere often merge together and form single structures in the 21cm sky with irregular shapes. Using high-resolution N-body LCDM simulations, we found that the lensing probability to have magnification bigger than 10 is ~10^{-5}. This results in $\ga 10^6$ strongly lensed Lyman-$\alpha$ spheres across the sky, which should be the primary targets for first detections of Lyman-$\alpha$ spheres. Although the required total radio array collecting area for their detection is large (~100 km^2), the design of long fixed cylindrical reflectors can significantly reduce the total cost of such array to the level of the square kilometer array (SKA) and makes the detection of these very first objects feasible.

The Sunyaev Zel'dovich (SZ) signal is highly non-Gaussian, so the SZ power spectrum (along with the mean $y$ parameter) does not provide a complete description of the SZ effect. Therefore, SZ-based constraints on cosmological parameters and on cluster gastrophysics which assume Gaussianity will be biased. We derive an analytic expression for the $n$-point joint PDF of the SZ power spectrum. Our derivation, which is based on the halo model, has several advantages: it is expressed in an integral form which allows quick computation; it is applicable to any given survey and any given angular scale; it is straightforward to incorporate many of the complexities which arise when modeling the SZ signal. To illustrate, we use our expression to estimate $p(C_\ell)$, the one-point PDF of the SZ power spectrum. For small sky coverage (applicable to BIMA/CBI and the Sunyaev Zel'dovich Array experiments), our analysis shows that $p(C_\ell)$ on the several arc-minute scale is expected to be strongly skewed, peaking at a value well below the mean and with a long tail which extends to tail high $C_\ell$ values. In the limit of large sky coverage (applicable to the South Pole Telescope and Planck), $p(C_\ell)$ approaches a Gaussian form. However, even in this limit, the variance of the power spectrum is very different from the naive Gaussian-based estimate. This is because different $\ell$ models are strongly correlated, making the cosmic variance of the SZ band-power much larger than the naive estimate. Our analysis should also be useful for modeling the PDF of the power spectrum induced by gravitational lensing at large $\ell$.

We study the impact of astrophysical processes on the gamma-ray background produced by the annihilation of dark matter particles in cosmological halos, with particular attention to the consequences of the formation of supermassive black holes. In scenarios where these objects form adiabatically from the accretion of matter on small seeds, dark matter is first compressed into very dense ``spikes'', then its density progressively decreases due to annihilations and scattering off of stellar cusps. With respect to previous analyses, based on non-evolving halos, the predicted annihilation signal is higher and significantly distorted at low energies, reflecting the large contribution to the total flux from unevolved spikes at high redshifts. The peculiar spectral feature arising from the specific redshift distribution of the signal, would discriminate the proposed scenario from more conventional astrophysical explanations. We discuss how this affects the prospects for detection and demonstrate that the gamma-ray background from DM annihilations might be detectable even in absence of a signal from the Galactic center.

For cosmologically interesting $f(R)$ gravity models, we derive the complete set of the linearized field equations in the Newtonian gauge, under environments of the solar system, galaxies and clusters respectively. Based on these equations, we confirmed previous $\gamma=1/2$ solution in the solar system. However, $f(R)$ gravity models can be strongly environment-dependent and the high density (comparing to the cosmological mean) solar system environment can excite a viable $\gamma=1$ solution for some $f(R)$ gravity models. Although for $f(R)\propto -1/R$, it is not the case; for $f(R)\propto -\exp(-R/\lambda_2H_0^2)$, such $\gamma=1$ solution does exist. This solution is virtually indistinguishable from that in general relativity (GR) and the value of the associated curvature approaches the GR limit, which is much higher than value in the $\gamma=1/2$ solution. We show that for some forms of $f(R)$ gravity, this solution is physically stable in the solar system and can smoothly connect to the surface of the Sun. The derived field equations can be applied directly to gravitational lensing of galaxies and clusters. We find that, despite significant difference in the environments of galaxies and clusters comparing to that of the solar system, gravitational lensing of galaxies and clusters can be virtually identical to that in GR, for some forms of $f(R)$ gravity. Fortunately, galaxy rotation curve and intra-cluster gas pressure profile may contain valuable information to distinguish these $f(R)$ gravity models from GR.

Extended 21cm absorption regions (dubbed ``21cm absorption halos'') around first galaxies at $z\sim 30$ are likely the first distinctive structures accessible to radio observations. Though the radio array capable of detecting and resolving them must have $\sim 200$ km$^2$ total collecting area, given the great impact of such detections to the understanding of the reionization process and cosmology, such radio survey would be extremely profitable. As an example, we point out a potentially useful byproduct of such survey. The resolved 21cm absorption ``halos'', likely close to spherical, can serve as (almost) ideal sources for measuring the {\it cosmic shear} and mapping the matter distribution to $z\sim 30$. We investigate the expected lensing signal and consider a variety of noise contributions on the shear measurement. We find that S/N $\sim 1$ can be achieved for individual ``halos''. Given millions of 21cm absorption ``halos'' across the sky, the total S/N will be comparable to traditional shear measurement of $\sim$$10^9$ galaxies at $z\sim 1$.

Cosmic dust extinction alters the flux of type Ia supernovae. Inhomogeneities in the dust distribution induce correlated fluctuations of the SN fluxes. We find that such correlation can be up to 60% of the signal caused by gravitational lensing magnification, with an opposite sign. Therefore if not corrected, cosmic dust extinction is the dominant source of systematic uncertainty for future SNe Ia lensing measurement limiting the overall S/N to be $\la 10$. On the other hand, SN flux correlation measurements can be used in combination with other lensing data to infer the level of dust extinction. This will provide a viable method to eliminate gray dust contamination from the SN Ia Hubble diagram.

A generic prediction of general relativity is that the cosmological linear density growth factor $D$ is scale independent. But in general, modified gravities do not preserve this signature. A scale dependent $D$ can cause time variation in gravitational potential at high redshifts and provides a new cosmological test of gravity, through early time integrated Sachs-Wolfe (ISW) effect-large scale structure (LSS) cross correlation. We demonstrate the power of this test for a class of $f(R)$ gravity, with the form $f(R)=-\lambda_1 H_0^2\exp(-R/\lambda_2H_0^2)$. Such $f(R)$ gravity, even with degenerate expansion history to $\Lambda$CDM, can produce detectable ISW effect at $z\ga 3$ and $l\ga 20$. Null-detection of such effect would constrain $\lambda_2$ to be $\lambda_2>1000$ at $>95%$ confidence level. On the other hand, robust detection of ISW-LSS cross correlation at high $z$ will severely challenge general relativity.

The decay rate of cosmological gravitational potential measures the deviation from Einstein-de Sitter universe and can put strong constraints on the nature of dark energy and gravity. Usual method to measure this decay rate is through the integrated Sachs-Wolfe (ISW) effect-large scale structure (LSS) cross correlation. However, the interpretation of the measured correlation signal is complicated by the galaxy bias and matter power spectrum. This could bias and/or degrade its constraints to the nature of dark energy and gravity. But, combining the lensing-LSS cross correlation measurements, the decay rate of gravitational potential can be isolated. For any given narrow redshift bin of LSS, the ratio of the two cross correlations directly measures $[d\ln D_{\phi}/d\ln a]H(z)/W(\chi,\chi_s)$, where $D_{\phi}$ is the linear growth factor of the gravitational potential, $H$ is the Hubble constant at redshift $z$, $W(\chi,\chi_s)$ is the lensing kernel and $\chi$ and $\chi_s$ are the comoving angular diameter distance to lens and source, respectively. This method is optimal in the sense that (1) the measured quantity is essentially free of systematic errors and is only limited by cosmic variance and (2) the measured quantity only depends on several cosmological parameters and can be predicted from first principles unambiguously. Though fundamentally limited by inevitably large cosmic variance associated with the ISW measurements, it can still put useful independent constraints on the amount of dark energy and its equation of state. It can also provide a powerful test of modified gravity and can distinguish the Dvali-Gabadadze-Porrati model from $\Lambda$CDM at $>2.5\sigma$ confidence level.

Future weak lensing measurements of cosmic shear will reach such high accuracy that second order effects in weak lensing modeling, like the influence of baryons on structure formation, become important. We use a controlled set of state of the art cosmological simulations to quantify this effect by comparing pure N-body dark matter runs with corresponding hydrodynamical simulations, carried out both in non-radiative, and in dissipative form with cooling and star formation. In both hydrodynamical simulations, the clustering of the gas is suppressed while that of dark matter is boosted at scales k>1 h/Mpc. Despite this counterbalance between dark matter and gas, the clustering of the total matter is suppressed by up to 1 percent at 1<k<10 h/Mpc, while for k ~ 20 h/Mpc it is boosted, up to 2 percent in the non-radiative run and 10 percent in the run with star formation. The stellar mass formed in the latter is highly biased relative to the dark matter in the pure N-body simulation. Using our power spectrum measurements to predict the effect of baryons on the weak lensing signal at 100<l<10000, we find that baryons may change the lensing power spectrum by less than 0.5 percent at l<1000, but by 1 to 10 percent at 1000<l<10000. The size of the effect exceeds the predicted accuracy of future lensing power spectrum measurements and will likely be detected. Precise determinations of cosmological parameters with weak lensing, and studies of small-scale fluctuations and clustering, therefore rely on properly including baryonic physics.

We show how precision lensing measurements can be obtained through the lensing magnification effect in high redshift 21cm emission from galaxies. Normally, cosmic magnification measurements have been seriously complicated by galaxy clustering. With precise redshifts obtained from 21cm emission line wavelength, one can correlate galaxies at different source planes, or exclude close pairs to eliminate such contaminations. We provide forecasts for future surveys, specifically the SKA and CLAR. SKA can achieve percent precision on the dark matter power spectrum and the galaxy dark matter cross correlation power spectrum, while CLAR can measure an accurate cross correlation power spectrum. The neutral hydrogen fraction was most likely significantly higher at high redshifts, which improves the number of observed galaxies significantly, such that also CLAR can measure the dark matter lensing power spectrum. SKA can also allow precise measurement of lensing bispectrum.

We develop a new tool to generate statistically precise dark matter maps from the cosmic magnification of galaxies with distance estimates. We show how to overcome the intrinsic clustering problem using the slope of the luminosity function, because magnificability changes strongly over the luminosity function, while intrinsic clustering only changes weakly. This may allow precision cosmology beyond most current systematic limitations. SKA is able to reconstruct projected matter density map at smoothing scale $\sim 10^{'}$ with S/N$\geq 1$, at the rate of 200-4000 deg$^2$ per year, depending on the abundance and evolution of 21cm emitting galaxies. This power of mapping dark matter is comparable to, or even better than that of cosmic shear from deep optical surveys or 21cm surveys.

In the mean field approximation of the relativistic $\sigma$-$\omega$-$\rho$ model, the magnetic fields are incorporated, and its influence on the properties of n-p-e neutron star matter are studied. When the strength of the magnetic field is weaker than $\sim 10^{18}$G, the particle fractions and chemical potentials, matter energy density and pressure hardly change with the magnetic field; when the strength of the magnetic field is stronger than $\sim 10^{20}$G, the above quantities change with the magnetic field evidently. Furthermore, the pressure is studied in both thermodynamics and hydrodynamics. The difference between these two ways exits in the high density region, that is, the thermal self-consistency may not be satisfied in this region if the magnetic field is considered.

We constrain the Cosmic Star Formation Rate (CSFR) by requiring that massive stars produce the observed UV, optical, and IR light while at the same time not overproduce the Diffuse Supernova Neutrino Background as bounded by Super-Kamiokande. With the massive star component so constrained we then show that a reasonable choice of stellar Initial Mass Function and other parameters results in SNIa rates and iron yields in good agreement with data. In this way we define a `concordance' CSFR that predicts the optical SNII rate and the SNIa contribution to the MeV Cosmic Gamma-Ray Background. The CSFR constrained to reproduce these and other proxies of intermediate and massive star formation is more clearly delineated than if it were measured by any one technique and has the following testable consequences: (1) SNIa contribute only a small fraction of the MeV Cosmic Gamma-Ray Background, (2) massive star core-collapse is nearly always accompanied by a successful optical SNII, and (3) the Diffuse Supernova Neutrino Background is tantalizingly close to detectability.

We consider the impact of second order corrections to the geodesic equation governing gravitational lensing. We start from the full second order metric, including scalar, vector and tensor perturbations, and retain all relevant contributions to the cosmic shear corrections that are second order in the gravitational potential. The relevant terms are: the nonlinear evolution of the scalar gravitational potential, the Born correction, and lens-lens coupling. No other second order terms contribute appreciably to the lensing signal. Since ray-tracing algorithims currently include these three effects, this derivation serves as rigorous justification for the numerical predictions.

We discuss the nature of the possible high-l excess in the Cosmic Microwave Background (CMB) anisotropy power spectrum observed by the Cosmic Background Imager (CBI). We probe the angular structure of the excess in the CBI deep fields and investigate whether it could be due to the scattering of CMB photons by hot electrons within clusters, the Sunyaev-Zeldovich (SZ) effect. We estimate the density fluctuation parameters for amplitude, sigma_8, and shape, Gamma, from CMB primary anisotropy data and other cosmological data. We use the results of two separate hydrodynamical codes for Lambda-CDM cosmologies, consistent with the allowed sigma_8 and Gamma values, to quantify the expected contribution from the SZ effect to the bandpowers of the CBI experiment and pass simulated SZ effect maps through our CBI analysis pipeline. The result is very sensitive to the value of sigma_8, and is roughly consistent with the observed power if sigma_8 ~ 1. We conclude that the CBI anomaly could be a result of the SZ effect for the class of Lambda-CDM concordance models if sigma_8 is in the upper range of values allowed by current CMB and Large Scale Structure (LSS) data.

Weak gravitational lensing of background galaxies offers an excellent opportunity to study the intervening distribution of matter. While much attention to date has focused on the two-point function of the cosmic shear, the three-point function, the {\it bispectrum}, also contains very useful cosmological information. Here, we compute three corrections to the bispectrum which are nominally of the same order as the leading term. We show that the corrections are small, so they can be ignored when analyzing present surveys. However, they will eventually have to be included for accurate parameter estimates from future surveys.

We show that by directly correlating the cluster kinetic Sunyaev Zeldovich (KSZ) flux, the cluster peculiar velocity power spectrum can be measured to $\sim 10%$ accuracy by future large sky coverage KSZ surveys. This method is almost free of systemics entangled in the usual velocity inversion method. The direct correlation brings extra information of density and velocity clustering. We utilize these information to construct two indicators of the Hubble constant and comoving angular distance and propose a novel method to constrain cosmology.

The measured cosmic gamma ray background (CGB) spectrum at MeV energies is in reasonable agreement with the predicted contribution from type Ia supernovae (SNIa). But the characteristic features in the SNIa gamma ray spectrum, weakened by integration over source redshifts, are hard to measure, and additionally the contributions from other sources in the MeV range are uncertain, so that the SNIa origin of the MeV CGB remains unproven. Since different CGB sources have different clustering properties and redshift distributions, by combining the CGB spectrum and angular correlation measurements, the contributions to the CGB could be identified and separated. The SNIa CGB large-scale structure follows that of galaxies. Its rms fluctuation at degree scales has a characteristic energy dependence, ranging from $\sim 1%$ to order of unity and can be measured to several percent precision by proposed future satellites such as the Advanced Compton Telescope. With the identification of the SNIa contribution, the SNIa rate could be measured unambiguously as a function of redshift up to $z \sim 1$, by combining both the spectrum and angular correlation measurements, yielding new constraints on the star formation rate to even higher redshifts. Finally, we show that the gamma ray and neutrino backgrounds from supernovae should be closely connected, allowing an important consistency test from the measured data. Identification of the astrophysical contributions to the CGB would allow much greater sensitivity to an isotropic high-redshift CGB contribution arising in extra dimension or dark matter models.

(Abridged) The precision modeling of the thermal Sunyaev Zeldovich (SZ) effect and future high accuracy measurements will provide a powerful way to constrain the thermal history of the universe. In this paper, we combine high resolution adiabatic hydro simulations and analytical models to provide a precision modeling of the gas density weighted temperature $\bar{T}_g$ and the mean SZ Compton $y$ parameter. Simulation artifacts are quantified and corrected and our analytical model passes tests of various simulations and is able to predict $\bar{T}_g$ and $\bar{y}$ to several percent accuracy. For low matter density $\Lambda$CDM cosmology, the present $\bar{T}_g$ is $0.32 (\sigma_8/0.84)^{3.05-0.15\Omega_m}(\Omega_m/0.268)^{1.28-0.2\sigma_8} {\rm keV}$. The mean $y$ parameter is $2.6\times 10^{-6} (\sigma_8/0.84)^{4.1-2\Omega_m}(\Omega_m/0.268)^{1.28-0.2\sigma_8}$. The current upper limit of $y<1.5\times 10^{-5}$ measured by FIRAS has already ruled out combinations of high $\sigma_8\ga 1.1$ and high $\Omega_m\ga 0.5$.

The cosmic infrared background (CIB) is a sensitive measure of the star formation history. But this background is overwhelmed by foregrounds, which bias the CIB mean flux and auto correlation measurement severely. Since dominant foregrounds do not correlate with the large scale structure, the cross correlation of CIB with galaxies is free of such foregrounds and presents as an unbiased measure of the star formation history. This also allows the measurement of CIB based on integrated intensity, whose theoretical prediction is based on energy conservation, thus is fairly model independent and robust. The redshift information of CIB sources can be directly recovered with the aid of galaxy photometric redshifts. In the far IR region, correlation signal is around 10% at degree scale. Combining FIRAS and SDSS, the cross correlation can be measured with 20% accuracy (statistical and systematics errors) at $l\sim 40$. Such measurement would allow the measurement of star formation rate with comparable accuracies to $z\sim 1$. Future CIB experiments with $1^{\circ}$ resolution are able to measure the star formation rate with 10% accuracy to $z\sim 1.5$. Secondary CMB anisotropies arising from the large scale structure correlate with CIB. The correlation strength is about $0.5\mu$K at $\sim 10^{\circ}$ (the integrated Sachs-Wolfe effect) and $\sim -0.3\mu$K at $\sim 1^{\circ}$ (the Sunyaev Zeldovich effect, Rayleigh-Jeans regime). FIRAS+WMAP would allow the measurement of the ISW effect with $\ga 5\sigma$ confidence. Its detection helps to constrain the amount of dark energy and far IR emissivity bias.

Precision measurements of the Sunyaev Zeldovich (SZ) effect in upcoming blank sky surveys require theoretical understanding of all physical processes with $\ga 10%$ effects on the SZ power spectrum. We show that, observed cluster magnetic field could reduce the SZ power spectrum by $\sim 20%$ at $l\sim 4000$, where the SZ power spectrum will be precisely measured by the Sunyaev Zeldovich array (SZA) and the Atacama cosmology telescope (ACT). At smaller scale, this effect is larger and could reach a factor of several. Such effect must be considered for an unbiased interpretation of the SZ data. Though the magnetic effect on the SZ power spectrum is very similar to that of radiative cooling, it is measurable by multi-band CMB polarization measurement.

The kinetic SZ effect, which is the dominant CMB source at arc-minute scales and $\nu \sim 217$ Ghz, probes the ionized gas peculiar momentum up to the epoch of reionization and is a sensitive measure of the reionization history. We ran high resolution self-similar and $\Lambda$CDM hydro simulations and built an analytical model to study this effect. Our model reproduces the $\Lambda$CDM simulation results to several percent accuracy, passes various tests against self-similar simulations, and shows a wider range of applicability than previous analytical models. Our model in its continuous version is free of simulation limitations such as finite simulation box and finite resolution and allows an accurate prediction of the kinetic SZ power spectrum $C_l$. For the WMAP cosmology, we find $l^2C_l/(2\pi)\simeq 0.91 \times 10^{-12} [(1+z_{\rm reion})/10]^{0.34}(l/5000)^{0.23-0.015(z_{\rm reion}-9)}$ for the reionization redshift $6<z_{\rm reion}<20$ and $3000<l<9000$. The corresponding temperature fluctuation is several $\mu$K at these ranges. The dependence of $C_l$ on the reionization history allows an accurate measurement of the reionization epoch. For the Atacama cosmology telescope experiment, $C_l$ can be measured with $\sim 1%$ accuracy. $C_l$ scales as $(\Omega_b h)^2 \sigma_8^{4\sim 6}$. Given cosmological parameters, ACT would be able to constrain $z_{\rm reion}$ with several percent accuracy. Some multi-reionization scenarios degenerate in the primary CMB temperature and TE measurement can be distinguished with $\sim 10 \sigma$ confidence.

Weak lensing measurements are entering a precision era to statistically map the distribution of matter in the universe. The most common measurement has been of the variance of the projected surface density of matter, which corresponds to the induced correlation in alignments of background galaxies. This measurement of the fluctuations is insensitive to the total mass content, like using waves on the ocean to measure its depths. But when the depth is shallow as happens near a beach, waves become skewed. Similarly, a measurement of skewness in the projected matter distribution directly measures the total matter content of the universe. While skewness has already been convincingly detected, its constraint on cosmology is still weak. We address optimal analyses for the CFHT Legacy Survey in the presence of noise. We show that a compensated Gaussian filter with a width of 2.5 arc minutes optimizes the cosmological constraint, yielding $\Delta \Omega_m/\Omega_m\sim 10%$. This is significantly better than other filters which have been considered in the literature. This can be further improved with tomography and other sophisticated analyses.

Weak gravitational lensing provides a direct statistical measure of the dark matter distribution. The variance is easiest to measure, which constrains the degenerate product \sigma_8\Omega^0.6. The degeneracy is broken by measuring the skewness arising from the fact that densities must remain positive, which is not possible when the initially symmetric perturbations become non-linear. Skewness measures the non-linear mass scale, which in combination with the variance measures \Omega directly. We present the first detection of dark matter skewness from the Virmos-Decart survey. We have measured the full three point function, and its projections onto windowed skewness. We separate the lensing mode and the B mode. The lensing skewness is detected for a compensated Gaussian on scales of 5.37 arc minutes to be \kappa^3=1.06+/-0.06x10^-6. The B-modes are consistent with zero at this scale. The variance for the same window function is \kappa^2= 5.32+/-0.62+/-0.98x10^-5, resulting in S_3=375^{+342}_{-124}. Comparing to N-body simulations, we find \Omega_0<0.5 at 90% confidence. The Canada-France-Hawaii-Telescope legacy survey and newer simulations should be able to improve significantly on the constraint.

We explore the intergalactic medium (IGM) as a potential source of the unresolved soft X-ray background (XRB) and the feasibility to extract the IGM state and evolution from XRB observations. We build two analytical models, the continuum field model and the halo model, to calculate the IGM XRB mean flux, angular auto correlation and cross correlation with galaxies. Our results suggest that the IGM may contribute a significant fraction to the unresolved soft XRB flux and correlations. We calibrated non-Gaussian errors estimated against our $512^3$ moving mesh hydro simulation and estimate that the ROSAT all sky survey plus Sloan galaxy photometric redshift survey would allow a $\sim 10%$ accuracy in the IGM XRB-galaxy cross correlation power spectrum measurement for $800<l<5000$ and a $\sim 20%$ accuracy in the redshift resolved X-ray emissivity-galaxy cross correlation power spectrum measurement for $z\lesssim 0.5$. At small scales, non-gravitational heating, e.g. feedback, dominates over gravity and leaves unique signatures in the IGM XRB, which allows a comparable accuracy in the measurement of the amount of non-gravitational heating and the length scales where non-gravitational energy balances gravity.

We review the recently published results from the CBI's first season of observations. Angular power spectra of the CMB were obtained from deep integrations of 3 single fields covering a total of 3 deg^2 and 3 shallower surveys of overlapping (mosaiced) fields covering a total of 40 deg^2. The observations show a damping of the anisotropies at high-l as expected from the standard scenarios of recombination. We present parameter estimates obtained from the data and discuss the significance of an excess at l>2000 observed in the deep fields.

We present a new approach to doing Eulerian computational fluid dynamics that is designed to work at high Mach numbers encountered in hydrodynamical simulations of the IGM. In conventional Eulerian CFD, the thermal energy is poorly tracked in supersonic bulk flows where local fluid variables cannot be accurately separated from the much larger bulk flow components. We described a method in which local fluid quantities can be directly tracked and the Eulerian fluid equations solved in a local frame moving with the flow. The new algorithm has been used to run large hydrodynamical simulations on a 1024^3 grid to study the kinetic SZ effect. The KSZ power spectrum is broadly peaked at l~10^4 with temperature fluctuations on micro Kelvin levels.

We derive cosmological parameters from the CBI measurements of the Cosmic Microwave Background (CMB) angular power spectrum. Our results provide an independent confirmation of the standard $\Omega_{\rm tot}=1$ $\Lambda$CDM model within the adiabatic, inflationary paradigm. Above $\ell=2000$ the observations show evidence of power in excess of that expected in the standard models. We use hydrodynamical simulations to show how Sunyaev-Zeldovich Effect (SZE) may account for the excess power for models with fluctuation amplitude $\sigma_8\sim 1$ which is in the high end of the range allowed by the primary CMB observations.

The Sunyaev Zel'dovich effect (SZ effect) is a complete probe of ionized baryons, the majority of which are likely hiding in the intergalactic medium. We ran a $512^3$ $\Lambda$CDM simulation using a moving mesh hydro code to compute the statistics of the thermal and kinetic SZ effect such as the power spectra and measures of non-Gaussianity. The thermal SZ power spectrum has a very broad peak at multipole $l\sim 2000-10^4$ with temperature fluctuations $\Delta T \sim 15\mu$K. The power spectrum is consistent with available observations and suggests a high $\sigma_8\simeq 1.0$ and a possible role of non-gravitational heating. The non-Gaussianity is significant and increases the cosmic variance of the power spectrum by a factor of $\sim 5$ for $l<6000$. We explore optimal driftscan survey strategies for the AMIBA CMB interferometer and their dependence on cosmology. For SZ power spectrum estimation, we find that the optimal sky coverage for a 1000 hours of integration time is several hundred square degrees. One achieves an accuracy better than 40% in the SZ measurement of power spectrum and an accuracy better than 20% in the cross correlation with Sloan galaxies for $2000<l<5000$. For cluster searches, the optimal scan rate is around 280 hours per square degree with a cluster detection rate 1 every 7 hours, allowing for a false positive rate of 20% and better than 30% accuracy in the cluster SZ distribution function measurement.

About 90% of baryons in the universe have thus far escaped direct observation. This is known as the {\it missing baryon problem}. The Sunyaev Zel'dovich effect (SZ effect) has the potential to directly measure the state of the majority of these missing baryons. The next generation CMB experiments such as AMIBA will provide an unbiased sample of the intergalactic medium through the SZ effect. The existing and upcoming simulations and analytical studies provide a quantitative understanding of the SZ effect. All these make the SZ effect a powerful baryon probe. We present an overview of this probe from both phenomelogical and theoretical aspects.

We apply the hierarchical clustering model and non-linear perturbation theory to the cosmological density and temperature fields. This allows us to calculate the intergalactic gas pressure power spectrum, SZ anisotropy power spectrum, skewness and related statistics. Then we show the effect of the non-gravitational heating. Our model confirms recent simulations yielding mass weighted gas temperature \bar{T}_g \sim 0.35 keV and reproduces the power spectra found in these simulations. While the SZ effect contains only angular information, we show that it is possible to extract the full time resolved gas pressure power spectrum when combined with galaxy photometric redshift surveys by using a variation on the cross-correlation. This method further allows the disentanglement of the gravitational and the non-gravitational heating.