Cosmological N-body simulations represent an excellent tool to study the formation and evolution of dark matter (DM) halos and the mechanisms that have originated the universal profile at the largest mass scales in the Universe. In particular, the combination of the velocity dispersion $\sigma_\mathrm{v}$ with the density $\rho$ can be used to define the pseudo-entropy $S(r)=\sigma_\mathrm{v}^2/\rho^{2/3}$, whose profile is well-described by a simple power-law $S\propto r^{\alpha}$. We analyze a set of cosmological hydrodynamical re-simulations of massive galaxy clusters and study the pseudo-entropy profiles as traced by different collisionless components in simulated galaxy clusters: DM, stars, and substructures. We analyze four sets of simulations, exploring different resolution and physics (N-body and full hydrodynamical simulations) to investigate convergence and the impact of baryons. We find that baryons significantly affect the inner region of pseudo-entropy profile as traced by substructures, while pseudo-entropy profiles as traced by DM particles are characterized by an almost universal behavior and can be used as a low-scatter ($\sim6-7\%$) mass-proxy. We compare observed and simulated pseudo-entropy profiles and find good agreement in both normalization and slope. We demonstrate, however, that the method used to derive observed pseudo-entropy profiles could introduce biases and underestimate the impact of mergers. Finally, we investigate the pseudo-entropy traced by the stars focusing our interest in the dynamical distinction between intracluster light (ICL) and the stars bound to the brightest cluster galaxy (BCG): the combination of these two pseudo-entropy profiles is well-described by a single power-law out to almost the entire cluster virial radius.

SN 2018hti is a Type I superluminous supernova (SLSN~I) with an absolute $g$-band magnitude of $-22.2$ at maximum brightness, discovered in a metal-poor galaxy at a redshift of 0.0612. We present extensive photometric and spectroscopic observations of this supernova, covering the phases from $\sim -35$ days to more than +340 days from the $r$-band maximum. Combining our $BVgri$-band photometry with {\it Swift} UVOT optical/ultraviolet photometry, we calculated the peak luminosity as $\sim 3.5\times10^{44}$ erg s$^{-1}$. Modeling the observed light curve reveals that the luminosity evolution of SN 2018hti can be produced by an ejecta mass of 5.8 $M_\odot$ and a magnetar with a magnetic field of $B=1.8\times10^{13}$~G having an initial spin period of $P_0=1.8$ ms. Based on such a magnetar-powered scenario and a larger sample, a correlation between the spin of the magnetar and the kinetic energy of the ejecta can be inferred for most SLSNe~I, suggesting a self-consistent scenario. Like for other SLSNe~I, the host galaxy of SN 2018hti is found to be relatively faint ($M_{g} = -17.75$ mag) and of low metallicity ($Z=0.3~Z_\odot$), with a star-formation rate of 0.3 $M_\odot$ yr$^{-1}$. According to simulation results of single-star evolution, SN 2018hti could originate from a massive, metal-poor star with a zero-age main sequence (ZAMS) mass of 25--40 $M_\odot$, or from a less massive rotating star with $M_\mathrm{ZAMS} \approx 16$--25 $M_\odot$. For the case of a binary system, its progenitor could also be a star with $M_\mathrm{ZAMS} \gtrsim 25$ $M_\odot$.

GRANDMA is a network of 25 telescopes of different sizes, including both photometric and spectroscopic facilities. The network aims to coordinate follow-up observations of gravitational-wave candidate alerts, especially those with large localisation uncertainties, to reduce the delay between the initial detection and the optical confirmation. In this paper, we detail GRANDMA's observational performance during Advanced LIGO/Advanced Virgo Observing Run 3 (O3), focusing on the second part of O3; this includes summary statistics pertaining to coverage and possible astrophysical origin of the candidates. To do so, we quantify our observation efficiency in terms of delay between gravitational-wave candidate trigger time, observations, and the total coverage. Using an optimised and robust coordination system, GRANDMA followed-up about 90 % of the gravitational-wave candidate alerts, i.e. 49 out of 56 candidates. This led to coverage of over 9000 deg2 during O3. The delay between the gravitational-wave candidate trigger and the first observation was below 1.5 hour for 50 % of the alerts. We did not detect any electromagnetic counterparts to the gravitational-wave candidates during O3, likely due to the very large localisation areas (on average thousands of degrees squares) and relatively large distance of the candidates (above 200 Mpc for 60 % of BNS candidates). We derive constraints on potential kilonova properties for two potential binary neutron star coalescences (GW190425 and S200213t), assuming that the events' locations were imaged.

We propose a more general version of the Trans-Planckian Censorship Conjecture (TCC) which can apply to models of inflation with varying speed of sound. We find that inflation models with $c_S < 1$ are in general more strongly constrained by censorship of trans-Planckian modes than canonical inflation models, with the upper bound on the tensor/scalar ratio reduced by as much as three orders of magnitude for sound speeds consistent with bounds from data. In particular, models which satisfy the TCC, and therefore the de Sitter Swampland Conjecture, can still violate the more general condition for non-classicality of trans-Planckian modes. As a concrete example, we apply the constraint to Dirac-Born-Infeld inflation models motivated by string theory.

The first measurements from the CIBER experiment of extragalactic background light (EBL) in near-infrared (NIR) band exhibit a higher intensity than those inferred through $\gamma$-ray observations. Recent theoretical-EBL intensities are typically consistent with the very high energy (VHE) $\gamma$-ray observations. Yet, it is possible that the excess NIR radiation is a new component of EBL and not in tension with the TeV spectra of distant blazars, since the hypothetical axion-like particle (ALP) may lead to a reduced opacity of the Universe for VHE $\gamma$-rays. In order to probe whether the excess component arises mainly from EBL, thirteen observed spectra in high energy and VHE ranges from ten distant TeV BL Lac objects are fitted by four theoretical spectra which involve theoretical EBL (Gilmore ), Gilmore's EBL model including photon/ALP coupling, Gilmore's EBL with CIBER excess and the latter including photon/ALP coupling respectively. We find the goodness of fit for the model with CIBER excess can be improved with a significance of $7.6~\sigma$ after including the photon/ALP coupling; Thus, the ALP/photon mixing mechanism can effectively alleviate the tension; However, the Gilmore EBL model, on the whole, is more compatible with the observed spectra compared to those with ALP, although individual blazars such as PKS 1424+240 and 1ES 1101-232 prefer the ALP-model. Our results suggest that the recent EBL models can solely explain the VHE $\gamma$-ray observations, and assuming the existence of the ALP to alleviate the tension is not required in a statistical sense, thus the excess over the EBL models is less likely to be a new EBL component.

After extensively explored, broad agreement between observations and theories has been reached that satellites are preferentially aligned with major axes of their host centrals. There are still some issues unsolved on this topic. In this paper, we present studies on satellite spatial distribution. To fairly compare with observations, we develop a novel galaxy finder and reconstruction algorithm in hydrodynamical simulation, which is based on the projected mock image, taking into account the full consideration of the point spread function, pixel size, surface brightness limit, resolution and redshift dimming effects. With galaxy samples constructed using such an algorithm, the satellite alignment is examined by comparing with observational results. It is found that the observational alignment can be reproduced for red galaxies, which dominate the sample in this study, but not for blue galaxies. Satellites' radial distribution is also investigated. It exhibits that outer satellites within host halos show stronger alignment signal than satellites in the inner regions, especially for red satellites, which is in contrast with previous studies. The disagreement is mainly due to extra galaxies identified by our new galaxy finder, which are mainly located in the inner region of host halos. Our study illustrates that at lower redshift, the alignment strength becomes stronger, while radial distribution curve becomes flatter. This suggests differences in the evolution of the angular distribution between satellites residing in the inner and outer halos, and implies that the post-infall evolution reduces the original alignment signal, that the impact decreases for satellites with later infall times.

We construct a SCUBA-2 450-$\mu$m map in the COSMOS field that covers an area of 300 arcmin$^{2}$ and reaches a 1$\sigma$ noise level of 0.65 mJy in the deepest region. We extract 256 sources detected at 450 $\mu$m with signal-to-noise ratio $>$ 4.0 and analyze the physical properties of their multi-wavelength counterparts. We find that most of the sources are at $z\lesssim3$, with a median of $z = 1.79^{+0.03}_{-0.15}$. About $35^{+32}_{-25}$% of our sources are classified as starburst galaxies based on their total star-formation rates (SFRs) and stellar masses ($M_{\ast}$). By fitting the far-infrared spectral energy distributions, we find that our 450-$\mu$m-selected sample has a wide range of dust temperatures (20 K $ \lesssim T_{\rm d} \lesssim$ 60 K), with a median of ${T}_{\rm d} = 38.3^{+0.4}_{-0.9}$ K. We do not find a redshift evolution in dust temperature for sources with $L_{\rm IR}$ > $10^{12}$ $\rm L_\odot$ at $z<3$. However, we find a moderate correlation where dust temperature increases with the deviation from the SFR-$M_{\ast}$ relation. The increase in dust temperature also correlates with optical morphology, which is consistent with merger-triggered starbursts in sub-millimeter galaxies. Our galaxies do not show the tight IRX-$\beta_{\rm UV}$ correlation that has been observed in the local Universe. We construct the infrared luminosity functions of our 450-$\mu$m sources and measure their comoving SFR densities. The contribution of the $L_{\rm IR}$ > $10^{12}$ $\rm L_\odot$ population to the SFR density rises dramatically from $z$ = 0 to 2 ($\propto$ ($1+z$)$^{3.9\pm1.1}$) and dominates the total SFR density at $z \gtrsim 2$.

The current Hubble constant tension is usually presented by comparing constraints on $H_0$ only. However, the post-recombination background cosmic evolution is determined by two parameters in the standard $\Lambda$CDM model, the Hubble constant ($H_0$) and today's matter energy fraction ($\Omega_{\rm{m}}$). If we therefore compare all constraints individually in the $H_0$-$\Omega_{\rm{m}}$ plane, (1) various constraints can be treated as independently as possible, (2) single-sided constraints are easier to consider, (3) compatibility among different constraints can be viewed in a more robust way, and (4) whether or not a nonstandard model is able to reconcile all constraints in tension can be seen more effectively. We perform a systematic comparison of independent constraints in the $H_0$-$\Omega_{\rm{m}}$ space based on a flat $\Lambda$CDM model. Constraints along different degeneracy directions consistently overlap in one region of the space, with the exception of the local measurement from Cepheid variable-calibrated supernovae. Due to the different responses of individual constraints to a modified model, it is difficult for nonstandard models with modifications at high-, mid- or low-redshifts to reconcile all constraints if none of them have unaccounted-for systematic effects. Based on our analysis, the local measurement is the most outlying and therefore drives the bulk of the tension. This may suggest that the most likely solution to the tension is an alteration to the local result, either due to some previously unseen feature in our local cosmic environment, or some other unknown systematic effect.

We present the Global Rapid Advanced Network Devoted to the Multi-messenger Addicts (GRANDMA). The network consists of 21 telescopes with both photometric and spectroscopic facilities. They are connected together thanks to a dedicated infrastructure. The network aims at coordinating the observations of large sky position estimates of transient events to enhance their follow-up and reduce the delay between the initial detection and the optical confirmation. The GRANDMA program mainly focuses on follow-up of gravitational-wave alerts to find and characterise the electromagnetic counterpart during the third observational campaign of the Advanced LIGO and Advanced Virgo detectors. But it allows for any follow-up of transient alerts involving neutrinos or gamma-ray bursts, even with poor spatial localisation. We present the different facilities, tools, and methods we developed for this network, and show its efficiency using observations of LIGO/Virgo S190425z, a binary neutron star merger candidate. We furthermore report on all GRANDMA follow-up observations performed during the first six months of the LIGO-Virgo observational campaign, and we derive constraints on the kilonova properties assuming that the events' locations were imaged by our telescopes.

Recent studies show that there is tension between the de Sitter swampland conjectures proposed by Obeid, et al. and inflationary cosmology. In this paper, we consider an alternative to inflation, `tachyacoustic' cosmology, in light of swampland conjectures. In tachyacoustic models, primordial perturbations are generated by a period of superluminal sound speed instead of accelerating expansion. We show that realizations of tachyacoustic Lagrangians can be consistent with the de Sitter swampland conjectures, and therefore can in principle be consistent with a UV-complete theory. We derive a general condition for models with $c_S > 1$ to be consistent with swampland conjectures.

Consistency between cosmological data sets is essential for ongoing and future cosmological analyses. We first investigate the questions of stability and applicability of some moment-based inconsistency measures to multiple data sets. We show that the recently introduced index of inconsistency (IOI) is numerically stable while it can be applied to multiple data sets. We use an illustrative construction of constraints as well as an example with real data sets (i.e. WMAP versus Planck) to show some limitations of the application of the Karhunen-Loeve decomposition to discordance measures. Second, we perform various consistency analyzes using IOI between multiple current data sets while \textit{working with the entire common parameter spaces}. We find current Large-Scale-Structure (LSS) data sets (Planck CMB lensing, DES lensing-clustering and SDSS RSD) all to be consistent with one another. This is found to be not the case for Planck temperature (TT) versus polarization (TE,EE) data, where moderate inconsistencies are present. Noteworthy, we find a strong inconsistency between joint LSS probes and Planck with IOI=5.27, and a moderate tension between DES and Planck with IOI=3.14. Next, using the IOI metric, we compare the Hubble constant from five independent probes. We confirm previous strong tensions between local measurement (SH0ES) and Planck as well as between H0LiCOW and Planck, but also find new strong tensions between SH0ES measurement and the joint LSS probes with IOI=6.73 (i.e. 3.7-$\sigma$ in 1D) as well as between joint LSS and combined probes SH0ES+H0LiCOW with IOI=8.59 (i.e. 4.1-$\sigma$ in 1D). Whether due to systematic effects in the data sets or problems with the underlying model, sources of these old and new tensions need to be identified and dealt with.

Some inconsistencies between cosmological observations continue to persist. The currently reported inconsistencies can be due to problems with the underlying model or systematic effects associated with some particular datasets. In an era of precision cosmology, it is important to develop proper tools to quantify the degree of these inconsistencies and to look for methods to identify their causes. This turns out to be a task that is not free of ambiguities, especially when the model is multi-dimensional, which is usually the case in cosmology. Measures currently-proposed in the literature disagree even in Gaussian cases. We discuss, with illustrative examples, some requirements that should be fulfilled in Gaussian cases by inconsistency measures and suggest a guiding definition of (in)consistency. As an example, we show that the recently-proposed index of inconsistency (IOI) meets those requirements and is in line with the (in)consistency definition proposed. Next, we examine the common practice to convert some measures, including IOI and other similar quantities, to a probability to exceed or a significance level that depends on the number of parameters. In the context of quantifying inconsistencies, we show that such a procedure can underestimate inconsistencies when there is more than one model parameter. We also discuss multiple-dataset comparisons and introduce a new tool based on the multi-dataset IOI that can identify outlying constraints when present. Comparison of constraints from various datasets can help identify the source of inconsistencies and will become important as more independent constraints will become available from ongoing and future surveys.

There has been considerable recent interest in a new class of non-slow roll inflationary solutions known as \textit{constant roll} inflation. Constant roll solutions are a generalization of the ultra-slow roll (USR) solution, where the first Hubble slow roll parameter $\epsilon$ is small, but the second Hubble slow roll parameter $\eta$ is not. While it is known that the USR solutions represent dynamical transients, there has been some disagreement in literature about whether or not large-$\eta$ constant roll solutions are attractors or are also a class of transient solutions. In this paper we show that the large-$\eta$ constant roll solutions do in fact represent transient solutions by performing stability analysis on the exact analytic (large-$\eta$) constant roll solutions.

One of the main science goals of the Large Synoptic Survey Telescope (LSST) is to uncover the nature of cosmic acceleration. In the base analysis, possible deviations from the Lambda-Cold-Dark-Matter ($\Lambda$CDM) background evolution will be probed by fitting a $w(z)$CDM model, which allows for a redshift-dependent dark energy equation of state with $w(z)$, within general relativity (GR). A rich array of other phenomena can arise due to deviations from the standard $\Lambda$CDM+GR model though, including modifications to the growth rate of structure and lensing, and novel screening effects on non-linear scales. Concrete physical models are needed to provide consistent predictions for these (potentially small) effects, to give us the best chance of detecting them and separating them from astrophysical systematics. A complex plethora of possible models has been constructed over the past few decades, with none emerging as a particular favorite. This document prioritizes a subset of these models along with rationales for further study and inclusion into the LSST Dark Energy Science Collaboration (DESC) data analysis pipelines, based on their observational viability, theoretical plausibility, and level of theoretical development. We provide references and theoretical expressions to aid the integration of these models into DESC software and simulations, and give justifications for why other models were not prioritized. While DESC efforts are free to pursue other models, we provide here guidelines on which theories appear to have higher priority for collaboration efforts due to their perceived promise and greater instructional value.

The Gamma-Ray Integrated Detectors (GRID) is a space mission concept dedicated to monitoring the transient gamma-ray sky in the energy range from 10 keV to 2 MeV using scintillation detectors onboard CubeSats in low Earth orbits. The primary targets of GRID are the gamma-ray bursts (GRBs) in the local universe. The scientific goal of GRID is, in synergy with ground-based gravitational wave (GW) detectors such as LIGO and VIRGO, to accumulate a sample of GRBs associated with the merger of two compact stars and study jets and related physics of those objects. It also involves observing and studying other gamma-ray transients such as long GRBs, soft gamma-ray repeaters, terrestrial gamma-ray flashes, and solar flares. With multiple CubeSats in various orbits, GRID is unaffected by the Earth occultation and serves as a full-time and all-sky monitor. Assuming a horizon of 200 Mpc for ground-based GW detectors, we expect to see a few associated GW-GRB events per year. With about 10 CubeSats in operation, GRID is capable of localizing a faint GRB like 170817A with a 90% error radius of about 10 degrees, through triangulation and flux modulation. GRID is proposed and developed by students, with considerable contribution from undergraduate students, and will remain operated as a student project in the future. The current GRID collaboration involves more than 20 institutes and keeps growing. On August 29th, the first GRID detector onboard a CubeSat was launched into a Sun-synchronous orbit and is currently under test.

We explain dark energy with equipartition theorem in string landscape.

According to the traditional scenario for core-collapse supernovae, the core of the collapsing star forms a neutron star and its gravitational energy release sends out a shockwave into the stellar envelope. However, in a significant number of numerical simulations, the shock stalls and the star cannot be exploded successfully, especially for a massive, compact star. We consider an alternative scenario that with mass fallback, the collapsing star forms a black hole in the center, surrounded by a dense, hot accretion disk, which blows out an intense outflow (wind). The kinetic energy of the wind may result in a successful stellar explosion. With an improved version of the formulism in Kohri et al. (2005) who studied neutron star accretion of minor fallback, we study this disk wind-driven explosion by calculating the accretion history for a suite of pre-SN stellar models with different initial surface rotational velocities, masses and metallicities, and by comparing the disk wind energy with the binding energy of the infalling stellar envelope. We show that the most promising models to be exploded successfully by this new channel are those relatively compact pre-SN stars with relatively low metallicities and not too low specific angular momenta. The total energies of the explosions are $\sim 10^{51-52}$ergs, and a more massive progenitor may produce a more energetic explosion.

The intrinsic alignment of galaxies is an important systematic effect in weak-lensing surveys, which can affect the derived cosmological parameters. One direct way to distinguish different alignment models and quantify their effects on the measurement is to produce mocked weak-lensing surveys. In this work, we use full-sky ray-tracing technique to produce mock images of galaxies from the ELUCID $N$-body simulation run with the WMAP9 cosmology. In our model we assume that the shape of central elliptical galaxy follows that of the dark matter halo, and spiral galaxy follows the halo spin. Using the mocked galaxy images, a combination of galaxy intrinsic shape and the gravitational shear, we compare the predicted tomographic shear correlations to the results of KiDS and DLS. It is found that our predictions stay between the KiDS and DLS results. We rule out a model in which the satellite galaxies are radially aligned with the center galaxy, otherwise the shear-correlations on small scales are too high. Most important, we find that although the intrinsic alignment of spiral galaxies is very weak, they induce a positive correlation between the gravitational shear signal and the intrinsic galaxy orientation (GI). This is because the spiral galaxy is tangentially aligned with the nearby large-scale overdensity, contrary to the radial alignment of elliptical galaxy. Our results explain the origin of detected positive GI term from the weak-lensing surveys. We conclude that in future analysis, the GI model must include the dependence on galaxy types in more detail.

The intra-cluster light (ICL) in observations is usually identified through the surface brightness limit method. In this paper, for the first time we produce the mock images of galaxy groups and clusters using a cosmological hydro- dynamical simulation, to investigate the ICL fraction and focus on its dependence on observational parameters, e.g., the surface brightness limit (SBL), the effects of cosmological redshift dimming, point spread function and CCD pixel size. Detailed analyses suggest that the width of point spread function has a significant effect on the measured ICL fraction, while the relatively small pixel size shows almost no influence. It is found that the measured ICL fraction depends strongly on the SBL. At a fixed SBL and redshift, the measured ICL fraction decreases with increasing halo mass, while with a much faint SBL, it does not depend on halo mass at low redshifts. In our work, the measured ICL fraction shows clear dependence on the cosmological redshift dimming effect. It is found that there are more mass locked in ICL component than light, suggesting that the use of a constant mass-to-light ratio at high surface brightness levels will lead to an underestimate of ICL mass. Furthermore, it is found that the radial profile of ICL shows a characteristic radius which is almost independent of halo mass. The current measurement of ICL from observations has a large dispersion due to different methods, and we emphasize the importance of using the same definition when observational results are compared with the theoretical predictions.

A subset of short gamma-ray bursts (sGRBs) have been found to be characterized by near-infrared/optical bumps at $\sim 1$ days, some of which exhibit almost concurrent X-ray flares. Although the near-infrared/optical bumps may be a signature of kilonovae, the X-ray flares are not consistent with kilonovae. It is widely believed that sGRBs are produced by the mergers of double compact objects, during which sub-relativistic ejecta are launched. In this paper we propose that the above optical/X-ray features are indicative of the formation of long-lived magnetars following the mergers of double neutron stars. Observations and theoretical works imply that the spin-down power of the magnetars is injected into the ejecta as ultra-relativistic electron-positron pairs, i.e. pulsar wind nebulae (PWNe). Here we suggest such a PWN model and find that the optical bump and X-ray flare observed in GRB 050724 can be well understood in this PWN model. We show that the optical bump and X-ray flare may have different origins. Our results strengthen the evidence for the formation of magnetars in double neutron star mergers and justify the validity of the PWN model.

We examine the quenched fraction of central and satellite galaxies as a function of galaxy stellar mass, halo mass, and the matter density of their large scale environment. Matter densities are inferred from our ELUCID simulation, a constrained simulation of local Universe sampled by SDSS, while halo masses and central/satellite classification are taken from the galaxy group catalog of Yang et al. The quenched fraction for the total population increases systematically with the three quantities. We find that the `environmental quenching efficiency', which quantifies the quenched fraction as function of halo mass, is independent of stellar mass. And this independence is the origin of the stellar mass-independence of density-based quenching efficiency, found in previous studies. Considering centrals and satellites separately, we find that the two populations follow similar correlations of quenching efficiency with halo mass and stellar mass, suggesting that they have experienced similar quenching processes in their host halo. We demonstrate that satellite quenching alone cannot account for the environmental quenching efficiency of the total galaxy population and the difference between the two populations found previously mainly arises from the fact that centrals and satellites of the same stellar mass reside, on average, in halos of different mass. After removing these halo-mass and stellar-mass effects, there remains a weak, but significant, residual dependence on environmental density, which is eliminated when halo assembly bias is taken into account. Our results therefore indicate that halo mass is the prime environmental parameter that regulates the quenching of both centrals and satellites.

Intrinsic alignments (IA) of galaxies have been recognized as one of the most serious contaminants to weak lensing. These systematics need to be isolated and mitigated in order for ongoing and future lensing surveys to reach their full potential. The IA self-calibration (SC) method was shown in previous studies to be able to reduce the GI contamination by up to a factor of 10 for the 2-point and 3-point correlations. The SC method does not require to assume an IA model in its working and can extract the GI signal from the same photo-z survey offering the possibility to test and understand structure formation scenarios and their relationship to IA models. In this paper, we study the effects of the IA SC mitigation method on the precision and accuracy of cosmological parameter constraints from future cosmic shear surveys LSST, WFIRST and Euclid. We perform analytical and numerical calculations to estimate the loss of precision and the residual bias in the best fit cosmological parameters after the self-calibration is performed. We take into account uncertainties from photometric redshifts and the galaxy bias. We find that the confidence contours are slightly inflated from applying the SC method itself while a significant increase is due to the inclusion of the photo-z uncertainties. The bias of cosmological parameters is reduced from several-$\sigma$, when IA is not corrected for, to below 1-$\sigma$ after SC is applied. These numbers are comparable to those resulting from applying the method of marginalizing over IA model parameters despite the fact that the two methods operate very differently. We conclude that implementing the SC for these future cosmic-shear surveys will not only allow one to efficiently mitigate the GI contaminant but also help to understand their modeling and link to structure formation.

We examine systematically the (in)consistency between cosmological constraints as obtained from various current data sets of the expansion history, Large Scale Structure (LSS), and Cosmic Microwave Background (CMB) from Planck. We run (dis)concordance tests within each set and across the sets using a recently introduced index of inconsistency (IOI) capable of dissecting inconsistencies between two or more data sets. First, we compare the constraints on $H_0$ from five different methods and find that the IOI drops from 2.85 to 0.88 (on Jeffreys' scales) when the local $H_0$ measurements is removed. This seems to indicate that the local measurement is an outlier, thus favoring a systematics-based explanation. We find a moderate inconsistency (IOI=2.61) between Planck temperature and polarization. We find that current LSS data sets including WiggleZ, SDSS RSD, CFHTLenS, CMB lensing and SZ cluster count, are consistent one with another and when all combined. However, we find a persistent moderate inconsistency between Planck and individual or combined LSS probes. For Planck TT+lowTEB versus individual LSS probes, the IOI spans the range 2.92--3.72 and increases to 3.44--4.20 when the polarization data is added in. The joint LSS versus the combined Planck temperature and polarization has an IOI of 2.83 in the most conservative case. But if Planck lowTEB is added to the joint LSS to constrain $\tau$ and break degeneracies, the inconsistency between Planck and joint LSS data increases to the high-end of the moderate range with IOI=4.81. Whether due to systematic effects in the data or to the underlying model, these inconsistencies need to be resolved. Finally, we perform forecast calculations using LSST and find that the discordance between Planck and future LSS data, if it persists as present, can rise up to a high IOI of 17, thus falling in the strong range of inconsistency. (Abridged).

The SCUBA-2 Ultra Deep Imaging EAO Survey (STUDIES) is a three-year JCMT Large Program aiming at reaching the 450 $\mu$m confusion limit in the COSMOS-CANDELS region, to study a representative sample of the high-redshift far-infrared galaxy population that gives rise to the bulk of the far-infrared background. We present the first-year data from STUDIES. We have reached a 450 $\mu$m noise level of 0.91~mJy for point sources at the map center, covered an area of 151 arcmin$^2$, and detected 98 and 141 sources at 4.0 and 3.5 $\sigma$, respectively. Our derived counts are best constrained in the 3.5-25 mJy regime using directly detected sources. Below the detection limits, our fluctuation analysis further constrains the slope of the counts down to 1 mJy. The resulting counts at 1-25 mJy are consistent with a power law having a slope of $-2.59$ ($\pm0.10$ for 3.5-25 mJy, and $^{+0.4}_{-0.7}$ for 1-3.5 mJy). There is no evidence of a faint-end termination or turn-over of the counts in this flux density range. Our counts are also consistent with previous SCUBA-2 blank-field and lensing cluster surveys. The integrated surface brightness from our counts down to 1 mJy is $90.0\pm17.2$ Jy deg$^{-2}$, which can account for up to $83^{+15}_{-16}\%$ of the COBE 450 $\mu$m background. We show that Herschel counts at 350 and 500 $\mu$m are significantly higher than our 450 $\mu$m counts, likely caused by its large beam and source clustering. High-angular resolution instruments like SCUBA-2 at 450 $\mu$m are therefore highly beneficial for measuring the luminosity and spatial density of high-redshift dusty galaxies.

It is unclear whether bulge growth is responsible for the flattening of the star formation main sequence (MS) at the high mass end. To investigate the role of bulges in shaping the MS, we compare the NUV$-r$ color between the central ($r<R_{50}$) and outer regions for a sample of 6401 local star-forming galaxies. The NUV$-r$ color is a good specific star formation rate indicator. We find that at $M_{\ast}<10^{10.2}M_{\sun}$, the central NUV$-r$ is on average only $\sim$ 0.25 mag redder than the outer NUV$-r$. Above $M_{\ast}=10^{10.2}M_{\sun}$, the central NUV$-r$ becomes systematically much redder than the outer NUV$-r$ for more massive galaxies, indicating that the central bulge is more evolved at the massive end. When dividing the galaxies according to their S\'ersic index $n$, we find that galaxies with $n$>2.0 tend to be redder in the central NUV$-r$ color than those with $n$<2.0, even at fixed B/T and $M_{\ast}$. This suggests that star formation in bulges is more strongly dependent on $n$ (or central mass density) than on B/T. Finally, we find that the fraction of galaxies with $n$>2.0 rapidly increases with $M_{\ast}$ at $M_{\ast}>10^{10.2}M_{\sun}$, which is consistent with the turning over of the MS at the same transition mass. We conclude that the increasing fraction of low-sSFR dense bulges in $M_{\ast}>10^{10.2}M_{\sun}$ galaxies, rather than increasing B/T, is responsible for the flattened slope of the $M_{\ast}$$-$SFR relation at high masses.

The continuous progress toward more precise cosmological surveys and experiments has galvanized recent interest into consistency tests on cosmological parameters and models. At the heart of this effort is quantifying the degree of inconsistency between two or more cosmological data sets. We introduce an intuitive moment-based measure we call the \textit{index of inconsistency} (IOI) and show that it is sensitive to the separation of the means, the size of the constraint ellipsoids, and their orientations in the parameter space. We find that it tracks accurately the inconsistencies when present. Next, we show that parameter marginalization can cause a loss of information on the inconsistency between two experiments and we quantify such a loss using the drop in IOI. In order to zoom on a given parameter, we define the relative residual IOI and the relative drop in IOI. While these two quantities can provide insights on the parameters that are most responsible for inconsistencies, we find that the full IOI applied to the whole parameter spaces is what must be used to correctly reflect the degree of inconsistency between two experiments. We discuss various properties of IOI, provide its eigenmode decomposition, and compare it to other measures of discordance. Finally, we apply IOI to current geometry data sets (i.e. an improved Supernovae Type Ia compilation, baryon acoustic oscillations from 6dF, SDSS MGS and Lyman-$\alpha$ forest, and high-$\ell$ CMB temperature data from Planck-2015) versus growth data sets (i.e. Redshift Space Distortions from WiggleZ and SDSS, Weak Lensing from CFHTLenS, CMB Lensing, Sunyav-Zeldovich effect, and low-$\ell$ CMB temperature and polarization data from Planck-2015). We find that a persistent inconsistency is present between the two data sets. This could reflect the presence of systematics in the data or inconsistencies in the underlying model.

The mixed Hydrogen and Helium (H + He) spectrum with a clear steepening at $\sim 700$ TeV has been detected by ARGO-YBJ experiments. In this paper, we demonstrate that the observed H + He spectrum can be well reproduced with the model of cosmic rays escaping from the supernova remnants (SNRs) in our Galaxy. In this model, particles are accelerated in a SNR through a non-linear diffusive shock acceleration mechanism and three components of high energy light nuclei escaped from the SNR are considered. It should be noted that the proton spectrum observed by KASCADE can be also explained by this model given a higher acceleration efficiency.

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.

Primordial gravitational waves constitute a promising probe of the very early Universe and the laws of gravity. We study in this work changes to tensor-mode perturbations (TMPs) that can arise in various proposed modified gravity (MG) theories. These include additional friction effects, nonstandard dispersion relations involving a massive graviton, a modified speed, and a small-scale modification. We introduce a physically motivated parameterization of these effects and use current available data to obtain exclusion regions in the parameter spaces. Taking into account the foreground subtraction, we then perform a forecast analysis focusing on tensor-mode MG parameters as constrained by future experiments COrE, Stage-IV and PIXIE. For a fiducial value of the tensor-to-scalar ratio $r=0.01$, we find that an additional friction of $3.5\sim4.5\%$ compared to GR will be detected at $3$-$\sigma$ by these experiments, while a decrease in friction will be more difficult to detect. The speed of gravitational waves needs to be by $5\sim15\%$ different from the speed of light for detection. We find that the minimum detectable graviton mass is about $7.8\sim9.7\times 10^{-33}$ eV, which is of the same order of magnitude as the graviton mass that allows massive gravity theories to produce cosmic acceleration. Finally, we study TMPs in MG during inflation using our parameterization. We find that the tensor spectral index is related to $r$ and the friction parameter $\nu_0$ by $n_T=-3\nu_0-r/8$. Assuming that the friction parameter is unchanged throughout the history of the Universe and is much larger than $r$, future experiments will be able to distinguish this MG consistency relation from the standard inflation one, and thus can be used as a further test of MG. In summary, TMPs and cosmic-microwave-background B-mode polarization provide a complementary avenue to test gravity theories. (Abridged)

One of the outstanding problems in general relativistic cosmology is that of the averaging. That is, how the lumpy universe that we observe at small scales averages out to a smooth Friedmann-Lemaitre-Robertson-Walker (FLRW) model. The root of the problem is that averaging does not commute with the Einstein equations that govern the dynamics of the model. This leads to the well-know question of backreaction in cosmology. In this work, we approach the problem using the covariant framework of Macroscopic Gravity (MG). We use its cosmological solution with a flat FLRW macroscopic background where the result of averaging cosmic inhomogeneities has been encapsulated into a backreaction density parameter denoted $\Omega_\mathcal{A}$. We constrain this averaged universe using available cosmological data sets of expansion and growth including, for the first time, a full CMB analysis from Planck temperature anisotropy and polarization data, the supernovae data from Union 2.1, the galaxy power spectrum from WiggleZ, the weak lensing tomography shear-shear cross correlations from the CFHTLenS survey and the baryonic acoustic oscillation data from 6Df, SDSS DR7 and BOSS DR9. We find that $-0.0155 \le \Omega_\mathcal{A} \le 0$ (at the 68\% CL) thus providing a tight upper-bound on the backreaction term. We also find that the term is strongly correlated with cosmological parameters such $\Omega_\Lambda$, $\sigma_8$ and $H_0$. While small, a backreaction density parameter of a few percent should be kept in consideration along with other systematics for precision cosmology.

The scenario consistent with a wealth of observations for the missing mass problem is that of weakly interacting dark matter particles. However, arguments or proposals for a Newtonian or relativistic modified gravity scenario continue to be made. A distinguishing characteristic between the two scenarios is that dark matter particles can produce a gravitational effect, in principle, without the need of baryons while this is not the case for the modified gravity scenario where such an effect must be correlated with the amount of baryonic matter. We consider here ultra-faint dwarf (UFD) galaxies as a promising arena to test the two scenarios based on the above assertion. We compare the correlation of the luminosity with the velocity dispersion between samples of UFD and non-UFD galaxies, finding a significant loss of correlation for UFD galaxies. For example, we find for 28 non-UFD galaxies a strong correlation coefficient of -0.688 which drops to -0.077 for the 23 UFD galaxies. Incoming and future data will determine whether the observed stochasticity for UFD galaxies is physical or due to systematics in the data. Such a loss of correlation (if it is to persist) is possible and consistent with the dark matter scenario for UFD galaxies but would constitute a new challenge for the modified gravity scenario.

In this work, we mainly study the magnification relations of quad lens models for cusp, fold and cross configurations. By dividing and ray-tracing in different image regions, we numerically derive the positions and magnifications of the four images for a point source lying inside of the astroid caustic. Then, based on the magnifications, we calculate the signed cusp and fold relations for the singular isothermal elliptical lenses. The signed fold relation map has positive and negative regions, and the positive region is usually larger than the negative region as has been confirmed before. It can also explain that for many observed fold image pairs, the fluxes of the Fermat minimum images are apt to be larger than those of the saddle images. We define a new quantity cross relation which describes the magnification discrepancy between two minimum images and two saddle images. Distance ratio is also defined as the ratio of the distance of two saddle images to that of two minimum images. We calculate the cross relations and distance ratios for nine observed Einstein crosses. In theory, for most of the quad lens models, the cross relations decrease as the distance ratios increase. In observation, the cross relations of the nine samples do not agree with the quad lens models very well, nevertheless, the cross relations of the nine samples do not give obvious evidence for anomalous flux ratio as the cusp and fold types do. Then, we discuss several reasons for the disagreement, and expect good consistencies for more precise observations and better lens models in the future.

Merger trees are routinely used to follow the growth and merging history of dark matter haloes and subhaloes in simulations of cosmic structure formation. Srisawat et al. (2013) compared a wide range of merger-tree-building codes. Here we test the influence of output strategies and mass resolution on tree-building. We find that, somewhat surprisingly, building the tree from more snapshots does not generally produce more complete trees; instead, it tends to short- en them. Significant improvements are seen for patching schemes which attempt to bridge over occasional dropouts in the underlying halo catalogues or schemes which combine the halo-finding and tree-building steps seamlessly. The adopted output strategy does not affec- t the average number of branches (bushiness) of the resultant merger trees. However, mass resolution has an influence on both main branch length and the bushiness. As the resolution increases, a halo with the same mass can be traced back further in time and will encounter more small progenitors during its evolutionary history. Given these results, we recommend that, for simulations intended as precursors for galaxy formation models where of order 100 or more snapshots are analysed, the tree-building routine should be integrated with the halo finder, or at the very least be able to patch over multiple adjacent snapshots.

Using a sample of ~6,000 local face-on star-forming galaxies (SFGs), we examine the correlations between the NUV-r colors both inside and outside the half-light radius, stellar mass M* and S\'{e}rsic index n in order to understand how the quenching of star formation is linked to galaxy structure. For these less dust-attenuated galaxies, NUV-r is found to be linearly correlated with Dn4000, supporting that NUV-r is a good photometric indicator of stellar age (or specific star formation rate). We find that: (1) At M*<10^{10.2}M_{\sun}, the central NUV-r is on average only~ 0.25 mag redder than the outer NUV-r. The intrinsic value would be even smaller after accounting for dust correction. However, the central NUV-r becomes systematically much redder than the outer NUV-r for more massive galaxies at M*>10^{10.2}M_{\sun}. (2) The central NUV-r shows no dependence on S\'{e}rsic index n at M*<10^{10.2}M_{\sun}, while above this mass galaxies with a higher n tend to be redder in the central NUV-r color. These results suggest that galaxies with M*<10^{10.2}M_{\sun} exhibit similar star formation activity from the inner R<R_{50} region to the R>R_{50} region. In contrast, a considerable fraction of the M*>10^{10.2}M_{\sun} galaxies, especially those with a high n, have harbored a relatively inactive bulge component.

The Extended Baryon Oscillation Spectroscopic Survey (eBOSS) will conduct novel cosmological observations using the BOSS spectrograph at Apache Point Observatory. Observations will be simultaneous with the Time Domain Spectroscopic Survey (TDSS) designed for variability studies and the Spectroscopic Identification of eROSITA Sources (SPIDERS) program designed for studies of X-ray sources. eBOSS will use four different tracers to measure the distance-redshift relation with baryon acoustic oscillations (BAO). Using more than 250,000 new, spectroscopically confirmed luminous red galaxies at a median redshift z=0.72, we project that eBOSS will yield measurements of $d_A(z)$ to an accuracy of 1.2% and measurements of H(z) to 2.1% when combined with the z>0.6 sample of BOSS galaxies. With ~195,000 new emission line galaxy redshifts, we expect BAO measurements of $d_A(z)$ to an accuracy of 3.1% and H(z) to 4.7% at an effective redshift of z= 0.87. A sample of more than 500,000 spectroscopically-confirmed quasars will provide the first BAO distance measurements over the redshift range 0.9<z<2.2, with expected precision of 2.8% and 4.2% on $d_A(z)$ and H(z), respectively. Finally, with 60,000 new quasars and re-observation of 60,000 quasars known from BOSS, we will obtain new Lyman-alpha forest measurements at redshifts z>2.1; these new data will enhance the precision of $d_A(z)$ and H(z) by a factor of 1.44 relative to BOSS. Furthermore, eBOSS will provide improved tests of General Relativity on cosmological scales through redshift-space distortion measurements, improved tests for non-Gaussianity in the primordial density field, and new constraints on the summed mass of all neutrino species. Here, we provide an overview of the cosmological goals, spectroscopic target sample, demonstration of spectral quality from early data, and projected cosmological constraints from eBOSS.

Eternal inflation arising from a potential landscape predicts that our universe is one realization of many possible cosmological histories. One way to access different cosmological histories is via the nucleation of bubble universes from a metastable false vacuum. Another way to sample different cosmological histories is via classical transitions, the creation of pocket universes through the collision between bubbles. Using relativistic numerical simulations, we examine the possibility of observationally determining if our observable universe resulted from a classical transition. We find that classical transitions produce spatially infinite, approximately open Friedman-Robertson-Walker universes. The leading set of observables in the aftermath of a classical transition are negative spatial curvature and a contribution to the Cosmic Microwave Background temperature quadrupole. The level of curvature and magnitude of the quadrupole are dependent on the position of the observer, and we determine the possible range of observables for two classes of single-scalar field models. For the first class, where the inflationary phase has a lower energy than the vacuum preceding the classical transition, the magnitude of the observed quadrupole generally falls to zero with distance from the collision while the spatial curvature grows to a constant. For the second class, where the inflationary phase has a higher energy than the vacuum preceding the classical transition, the magnitude of the observed quadrupole generically falls to zero with distance from the collision while the spatial curvature grows without bound. We find that the magnitude of the quadrupole and curvature grow with increasing centre of mass energy of the collision, and explore variations of the parameters in the scalar field lagrangian.

In this Letter, we investigate how galaxy mass assembly mode depends on stellar mass $M_{\ast}$, using a large sample of $\sim$10, 000 low redshift galaxies. Our galaxy sample is selected to have SDSS $R_{90}>5\arcsec.0$, which allows the measures of both the integrated and the central NUV$-r$ color indices. We find that: in the $M_{\ast}-($ NUV$-r$) green valley, the $M_{\ast}<10^{10}~M_{\sun}$ galaxies mostly have positive or flat color gradients, while most of the $M_{\ast}>10^{10.5}~M_{\sun}$ galaxies have negative color gradients. When their central $D_{n}4000$ index values exceed 1.6, the $M_{\ast}<10^{10.0}~M_{\sun}$ galaxies have moved to the UV red sequence, whereas a large fraction of the $M_{\ast}>10^{10.5}~M_{\sun}$ galaxies still lie on the UV blue cloud or the green valley region. We conclude that the main galaxy assembly mode is transiting from "the outside-in" mode to "the inside-out" mode at $M_{\ast}< 10^{10}~M_{\sun}$ and at $M_{\ast}> 10^{10.5}~M_{\sun}$. We argue that the physical origin of this is the compromise between the internal and the external process that driving the star formation quenching in galaxies. These results can be checked with the upcoming large data produced by the on-going IFS survey projects, such as CALIFA, MaNGA and SAMI in the near future.

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.

We review five often used quad lens models, each of which has analytical solutions and can produce four images at most. Each lens model has two parameters, including one that describes the intensity of non-dimensional mass density, and the other one that describes the deviation from the circular lens. In our recent work, we have found that the cusp and the fold summations are not equal to 0, when a point source infinitely approaches a cusp or a fold from inner side of the caustic. Based on the magnification invariant theory, which states that the sum of signed magnifications of the total images of a given source is a constant, we calculate the cusp summations for the five lens models. We find that the cusp summations are always larger than 0 for source on the major cusps, while can be larger or smaller than 0 for source on the minor cusps. We also find that if these lenses tend to the circular lens, the major and minor cusp summations will have infinite values, and with positive and negative signs respectively. The cusp summations do not change significantly if the sources are slightly deviated from the cusps. In addition, through the magnification invariants, we also derive the analytical signed cusp relations on the axes for three lens models. We find that both on the major and the minor axes the larger the lenses deviated from the circular lens, the larger the signed cusp relations. The major cusp relations are usually larger than the absolute minor cusp relations, but for some lens models with very large deviation from circular lens, the minor cusp relations can be larger than the major cusp relations.

Investigating the spin parameter distribution of subhaloes in two high resolution isolated halo simulations, re- cent work by Onions et al. suggested that typical subhalo spins are consistently lower than the spin distribution found for field haloes. To further examine this puzzle, we have analyzed simulations of a cosmological volume with sufficient resolution to resolve a significant subhalo population. We confirm the result of Onions et al. and show that the typical spin of a subhalo decreases with decreasing mass and increasing proximity to the host halo center. We interpret this as the growing influence of tidal stripping in removing the outer layers, and hence the higher angular momentum particles, of the subhaloes as they move within the host potential. Investigating the redshift dependence of this effect, we find that the typical subhalo spin is smaller with decreasing redshift. This indicates a temporal evolution as expected in the tidal stripping scenario.

Galaxies are not randomly distributed across the universe but showing different kinds of alignment on different scales. On small scales satellite galaxies have a tendency to distribute along the major axis of the central galaxy, with dependence on galaxy properties that both red satellites and centrals have stronger alignment than their blue counterparts. On large scales, it is found that the major axes of Luminous Red Galaxies (LRGs) have correlation up to 30Mpc/h. Using hydro-dynamical simulation with star formation, we investigate the origin of galaxy alignment on different scales. It is found that most red satellite galaxies stay in the inner region of dark matter halo inside which the shape of central galaxy is well aligned with the dark matter distribution. Red centrals have stronger alignment than blue ones as they live in massive haloes and the central galaxy-halo alignment increases with halo mass. On large scales, the alignment of LRGs is also from the galaxy-halo shape correlation, but with some extent of mis-alignment. The massive haloes have stronger alignment than haloes in filament which connect massive haloes. This is contrary to the naive expectation that cosmic filament is the cause of halo alignment.

Simulating the evolution of the local universe is important for studying galaxies and the intergalactic medium in a way free of cosmic variance. Here we present a method to reconstruct the initial linear density field from an input non-linear density field, employing the Hamiltonian Markov Chain Monte Carlo (HMC) algorithm combined with Particle Mesh (PM) dynamics. The HMC+PM method is applied to cosmological simulations, and the reconstructed linear density fields are then evolved to the present day with N-body simulations. The constrained simulations so obtained accurately reproduce both the amplitudes and phases of the input simulations at various $z$. Using a PM model with a grid cell size of 0.75 Mpc/h and 40 time-steps in the HMC can recover more than half of the phase information down to a scale k~0.85 h/Mpc at high z and to k~3.4 h/Mpc at z=0, which represents a significant improvement over similar reconstruction models in the literature, and indicates that our model can reconstruct the formation histories of cosmic structures over a large dynamical range. Adopting PM models with higher spatial and temporal resolutions yields even better reconstructions, suggesting that our method is limited more by the availability of computer resource than by principle. Dynamic models of structure evolution adopted in many earlier investigations can induce non-Gaussianity in the reconstructed linear density field, which in turn can cause large systematic deviations in the predicted halo mass function. Such deviations are greatly reduced or absent in our reconstruction.

Observations have shown that the spatial distribution of satellite galaxies is not random, but rather is aligned with the major axes of central galaxies (CGs). The strength of the alignment is dependent on the properties of both the satellites and centrals. Theoretical studies using dissipationless N-body simulations are limited by their inability to directly predict the shape of CGs. Using hydrodynamical simulations including gas cooling, star formation, and feedback, we carry out a study of galaxy alignment and its dependence on the galaxy properties predicted directly from the simulations.We found that the observed alignment signal is well produced, as is the color dependence: red satellites and red centrals both show stronger alignments than their blue counterparts. The reason for the stronger alignment of red satellites is that most of them stay in the inner region of the dark matter halo where the shape of the CG better traces the dark matter distribution. The dependence of alignment on the color of CGs arises from the halo mass dependence, since the alignment between the shape of the central stellar component and the inner halo increases with halo mass. We also find that the alignment of satellites is most strongly dependent on their metallicity, suggesting that the metallicity of satellites, rather than color, is a better tracer of galaxy alignment on small scales. This could be tested in future observational studies.

In this Letter, we analyse the radial UV-optical color distributions in a sample of low redshift green valley (GV) galaxies, with the Galaxy Evolution Explorer (GALEX)+Sloan Digital Sky Survey (SDSS) images, to investigate how the residual recent star formation distribute in these galaxies. We find that the dust-corrected $u-r$ colors of early-type galaxies (ETGs) are flat out to $R_{90}$, while the colors turn blue monotonously when $r>0.5R_{50}$ for late-type galaxies (LTGs). More than a half of the ETGs are blue-cored and have remarkable positive NUV$-r$ color gradients, suggesting that their star formation are centrally concentrated; the rest have flat color distributions out to $R_{90}$. The centrally concentrated star formation activity in a large portion of ETGs is confirmed by the SDSS spectroscopy, showing that $\sim$50 % ETGs have EW(H$\rm \alpha$)$>6.0$ \AA. For the LTGs, 95% of them show uniform radial color profiles, which can be interpreted as a red bulge plus an extended blue disk. The links between the two kinds of ETGs, e.g., those objects having remarkable "blue-cored" and those having flat color gradients, are less known and require future investigations. It is suggested that the LTGs follow a general picture that quenching first occur in the core regions, and then finally extend to the rest of the galaxy. Our results can be re-examined and have important implications for the IFU surveys, such as MaNGA and SAMI.

Observations have shown that the spatial distribution of satellite galaxies is not random, but aligned with the major axes of central galaxies. This alignment is dependent on galaxy properties, such that red satellites are more strongly aligned than blue satellites. Theoretical work done to interpret this phenomena has found that it is due to the non-spherical nature of dark matter halos. However, most studies over-predict the alignment signal under the assumption that the central galaxy shape follows the shape of the host halo. It is also not clear whether the color dependence of alignment is due to an assembly bias or an evolution effect. In this paper we study these problems using a cosmological N-body simulation. Subhalos are used to trace the positions of satellite galaxies. It is found that the shape of dark matter halos are mis-aligned at different radii. If the central galaxy shares the same shape as the inner host halo, then the alignment effect is weaker and agrees with observational data. However, it predicts almost no dependence of alignment on the color of satellite galaxies, though the late accreted subhalos show stronger alignment with the outer layer of the host halo than their early accreted counterparts. We find that this is due to the limitation of pure N-body simulations that satellites galaxies without associated subhalos ('orphan galaxies') are not resolved. These orphan (mostly red) satellites often reside in the inner region of host halos and should follow the shape of the host halo in the inner region.

A thorough understanding of many astrophysical phenomena associated with compact objects requires reliable knowledge about both the equation of state (EOS)of super-dense nuclear matter and the theory of strong-field gravity simultaneously because of the EOS-gravity degeneracy. Currently, variations of the neutron star (NS) mass-radius correlation from using alternative gravity theories are much larger than those from changing the NS matter EOS within known constraints. At least two independent observables are required to break the EOS-gravity degeneracy. Using model EOSs for hybrid stars and a Yukawa-type non-Newtonian gravity, we investigate both the mass-radius correlation and pulsating frequencies of NSs. While the maximum mass of NSs increases, the frequencies of the $f$, $p_1$, $p_2$, and $w_I$ pulsating modes are found to decrease with the increasing strength of the Yukawa-type non-Newtonian gravity, providing a useful reference for future determination simultaneously of both the strong-field gravity and the supranuclear EOS by combining data of x-ray and gravitational wave emissions of neutron stars.

Using analytical method, we study the Singular Isothermal Quadrupole (SIQ) lens system, which is the simplest lens model that can produce four images. In this case, the radial mass distribution is in accord with the profile of the Singular Isothermal Sphere (SIS) lens, and the tangential distribution is given by adding a quadrupole on the monopole component. The basic properties of the SIQ lens have been studied in this paper, including deflection potential, deflection angle, magnification, critical curve, caustic, pseudo-caustic and transition locus. Analytical solutions of the image positions and magnifications for the source on axes are derived. As have been found, naked cusps will appear when the relative intensity $k$ of quadrupole to monopole is larger than 0.6. According to the magnification invariant theory of the SIQ lens, the sum of the signed magnifications of the four images should be equal to unity \citep{dal98}. However, if a source lies in the naked cusp, the summed magnification of the left three images is smaller than the invariant 1. With this simple lens system, we study the situations that a point source infinitely approaches a cusp or a fold. The sum of magnifications of cusp image triplet is usually not equal to 0, and it is usually positive for major cusp while negative for minor cusp. Similarly, the sum of magnifications of fold image pair is usually neither equal to 0. Nevertheless, the cusp and fold relations are still equal to 0, in that the sum values are divided by infinite absolute magnifications by definition.

An arbitrary surface mass density of gravitational lens can be decomposed into multipole components. We simulate the ray-tracing for the multipolar mass distribution of generalized SIS (Singular Isothermal Sphere) model, based on the deflection angles which are analytically calculated. The magnification patterns in the source plane are then derived from inverse shooting technique. As have been found, the caustics of odd mode lenses are composed of two overlapping layers for some lens models. When a point source traverses such kind of overlapping caustics, the image numbers change by \pm 4, rather than \pm 2. There are two kinds of images for the caustics. One is the critical curve and the other is the transition locus. It is found that the image number of the fold is exactly the average value of image numbers on two sides of the fold, while the image number of the cusp is equal to the smaller one. We also focus on the magnification patterns of the quadrupole (m = 2) lenses under the perturbations of m = 3, 4 and 5 mode components, and found that one, two, and three butterfly or swallowtail singularities can be produced respectively. With the increasing intensity of the high-order perturbations, the singularities grow up to bring sixfold image regions. If these perturbations are large enough to let two or three of the butterflies or swallowtails contact, eightfold or tenfold image regions can be produced as well. The possible astronomical applications are discussed.

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.

This document summarizes the results of a community-based discussion of the potential science impact of the Mayall+BigBOSS highly multiplexed multi-object spectroscopic capability. The KPNO Mayall 4m telescope equipped with the DOE- and internationally-funded BigBOSS spectrograph offers one of the most cost-efficient ways of accomplishing many of the pressing scientific goals identified for this decade by the "New Worlds, New Horizons" report. The BigBOSS Key Project will place unprecedented constraints on cosmological parameters related to the expansion history of the universe. With the addition of an open (publicly funded) community access component, the scientific impact of BigBOSS can be extended to many important astrophysical questions related to the origin and evolution of galaxies, stars, and the IGM. Massive spectroscopy is the critical missing ingredient in numerous ongoing and planned ground- and space-based surveys, and BigBOSS is unique in its ability to provide this to the US community. BigBOSS data from community-led projects will play a vital role in the education and training of students and in maintaining US leadership in these fields of astrophysics. We urge the NSF-AST division to support community science with the BigBOSS multi-object spectrograph through the period of the BigBOSS survey in order to ensure public access to the extraordinary spectroscopic capability.

Brown dwarfs are important objects because they may provide a missing link between stars and planets, two populations that have dramatically different formation history. In this paper, we present the candidate binaries with brown dwarf companions that are found by analyzing binary microlensing events discovered during 2004 - 2011 observation seasons. Based on the low mass ratio criterion of q < 0.2, we found 7 candidate events, including OGLE-2004-BLG-035, OGLE-2004-BLG-039, OGLE-2007-BLG-006, OGLE-2007-BLG-399/MOA-2007-BLG-334, MOA-2011-BLG-104/OGLE-2011-BLG-0172, MOA-2011-BLG-149, and MOA-201-BLG-278/OGLE-2011-BLG-012N. Among them, we are able to confirm that the companions of the lenses of MOA-2011-BLG-104/OGLE-2011-BLG-0172 and MOA-2011-BLG-149 are brown dwarfs by determining the mass of the lens based on the simultaneous measurement of the Einstein radius and the lens parallax. The measured mass of the brown dwarf companions are (0.02 +/- 0.01) M_Sun and (0.019 +/- 0.002) M_Sun for MOA-2011-BLG-104/OGLE-2011-BLG-0172 and MOA-2011-BLG-149, respectively, and both companions are orbiting low mass M dwarf host stars. More microlensing brown dwarfs are expected to be detected as the number of lensing events with well covered light curves increases with new generation searches.

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.

We present the analysis of the light curves of 9 high-magnification single-lens gravitational microlensing events with lenses passing over source stars, including OGLE-2004-BLG-254, MOA-2007-BLG-176, MOA-2007-BLG-233/OGLE-2007-BLG-302, MOA-2009-BLG-174, MOA-2010-BLG-436, MOA-2011-BLG-093, MOA-2011-BLG-274, OGLE-2011-BLG-0990/MOA-2011-BLG-300, and OGLE-2011-BLG-1101/MOA-2011-BLG-325. For all events, we measure the linear limb-darkening coefficients of the surface brightness profile of source stars by measuring the deviation of the light curves near the peak affected by the finite-source effect. For 7 events, we measure the Einstein radii and the lens-source relative proper motions. Among them, 5 events are found to have Einstein radii less than 0.2 mas, making the lenses candidates of very low-mass stars or brown dwarfs. For MOA-2011-BLG-274, especially, the small Einstein radius of $\theta_{\rm E}\sim 0.08$ mas combined with the short time scale of $t_{\rm E}\sim 2.7$ days suggests the possibility that the lens is a free-floating planet. For MOA-2009-BLG-174, we measure the lens parallax and thus uniquely determine the physical parameters of the lens. We also find that the measured lens mass of $\sim 0.84\ M_\odot$ is consistent with that of a star blended with the source, suggesting that the blend is likely to be the lens. Although we find planetary signals for none of events, we provide exclusion diagrams showing the confidence levels excluding the existence of a planet as a function of the separation and mass ratio.

The microlensing event OGLE-2008-BLG-510 is characterised by an evident asymmetric shape of the peak, promptly detected by the ARTEMiS system in real time. The skewness of the light curve appears to be compatible both with binary-lens and binary-source models, including the possibility that the lens system consists of an M dwarf orbited by a brown dwarf. The detection of this microlensing anomaly and our analysis demonstrates that: 1) automated real-time detection of weak microlensing anomalies with immediate feedback is feasible, efficient, and sensitive, 2) rather common weak features intrinsically come with ambiguities that are not easily resolved from photometric light curves, 3) a modelling approach that finds all features of parameter space rather than just the `favourite model' is required, and 4) the data quality is most crucial, where systematics can be confused with real features, in particular small higher-order effects such as orbital motion signatures. It moreover becomes apparent that events with weak signatures are a silver mine for statistical studies, although not easy to exploit. Clues about the apparent paucity of both brown-dwarf companions and binary-source microlensing events might hide here.

Using N-body simulations and galaxy formation models, we study the galaxy stellar mass correlation and the two-point auto-correlation. The simulations are run with cosmological parameters from the WMAP first, third and seven year results, which mainly differ in the perturbation amplitude of \sigma_{8}. The stellar mass of galaxies are determined using either a semi-analytical galaxy formation model or a simple empirical abundance matching method. Compared to the SDSS DR7 data at z=0 and the DEEP2 results at z=1, we find that the predicted galaxy clusterings from the semi-analytical model are higher than the data at small scales, regardless of the adopted cosmology. Conversely, the abundance matching method predicts good agreement with the data at both z=0 and z=1 for high \sigma_8 cosmologies (WMAP1 & WMAP7), but the predictions from a low \sigma_8 cosmology (WMAP3) are significantly lower than the data at z=0. We find that the excess clustering at small-scales in the semi-analytical model mainly arises from satellites in massive haloes, indicating that either the star formation is too efficient in low-mass haloes or tidal stripping is too inefficient at high redshift. Our results show that galaxy clustering is strongly affected by the models for galaxy formation, thus can be used to constrain the baryonic physics. The weak dependence of galaxy clustering on cosmological parameters makes it difficult to constrain the WMAP1 and WMAP7 cosmologies.

Microlensing can provide a useful tool to probe binary distributions down to low-mass limits of binary companions. In this paper, we analyze the light curves of 8 binary lensing events detected through the channel of high-magnification events during the seasons from 2007 to 2010. The perturbations, which are confined near the peak of the light curves, can be easily distinguished from the central perturbations caused by planets. However, the degeneracy between close and wide binary solutions cannot be resolved with a $3\sigma$ confidence level for 3 events, implying that the degeneracy would be an important obstacle in studying binary distributions. The dependence of the degeneracy on the lensing parameters is consistent with a theoretic prediction that the degeneracy becomes severe as the binary separation and the mass ratio deviate from the values of resonant caustics. The measured mass ratio of the event OGLE-2008-BLG-510/MOA-2008-BLG-369 is $q\sim 0.1$, making the companion of the lens a strong brown-dwarf candidate.

We investigate the solar-system constraint on the f(R) theory of modified gravity with chameleon mechanism, where f(R) represents the deviation from general relativity in the gravity action. We obtain a stringent bound to a general, non-constant deviation function f(R): -10^{-15} < df/dR < 0 when R ~ 3*10^5*H0^2, and a loose bound: 0 < R*d(df/dR)/dR < 2/5 when R > 3*10^5*H0^2, by requiring the thin-shell condition in the solar system, particularly in the atmosphere of the Earth. These bounds can be conveniently utilized to test the f(R) models with given functional forms of f(R) and to obtain the constraints on the parameters therein. For demonstration we apply these bounds to several widely considered f(R) models. (H0: Hubble constant)

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.

We present the discovery and mass measurement of the cold, low-mass planet MOA-2009-BLG-266Lb, made with the gravitational microlensing method. This planet has a mass of m_p = 10.4 +- 1.7 Earth masses and orbits a star of mass M_* = 0.56 +- 0.09 Solar masses at a semi-major axis of a = 3.2 (+1.9 -0.5) AU and an orbital period of P = 7.6 (+7.7 -1.5} yrs. The planet and host star mass measurements are enabled by the measurement of the microlensing parallax effect, which is seen primarily in the light curve distortion due to the orbital motion of the Earth. But, the analysis also demonstrates the capability to measure microlensing parallax with the Deep Impact (or EPOXI) spacecraft in a Heliocentric orbit. The planet mass and orbital distance are similar to predictions for the critical core mass needed to accrete a substantial gaseous envelope, and thus may indicate that this planet is a "failed" gas giant. This and future microlensing detections will test planet formation theory predictions regarding the prevalence and masses of such planets.

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.

Since 1995, more than 500 exoplanets have been detected using different techniques, of which 11 were detected with gravitational microlensing. Most of these are gravitationally bound to their host stars. There is some evidence of free-floating planetary mass objects in young star-forming regions, but these objects are limited to massive objects of 3 to 15 Jupiter masses with large uncertainties in photometric mass estimates and their abundance. Here, we report the discovery of a population of unbound or distant Jupiter-mass objects, which are almost twice (1.8_{-0.8}^{+1.7}) as common as main-sequence stars, based on two years of gravitational microlensing survey observations toward the Galactic Bulge. These planetary-mass objects have no host stars that can be detected within about ten astronomical units by gravitational microlensing. However a comparison with constraints from direct imaging suggests that most of these planetary-mass objects are not bound to any host star. An abrupt change in the mass function at about a Jupiter mass favours the idea that their formation process is different from that of stars and brown dwarfs. They may have formed in proto-planetary disks and subsequently scattered into unbound or very distant orbits.

We report the result of the analysis of the light curve of a caustic-crossing binary-lens microlensing event OGLE-2009-BLG-023/MOA-2009-BLG-028. Even though the event was observed solely by survey experiments, we could uniquely determine the mass of the lens and distance to it by simultaneously measuring the Einstein radius and lens parallax. From this, we find that the lens system is composed of M-type dwarfs with masses $(0.50\pm 0.07) \ M_\odot$ and $(0.15\pm 0.02)\ M_\odot$ located in the Galactic disk with a distance of $\sim 1.8$ kpc toward the Galactic bulge direction. The event demonstrates that physical lens parameters of binary-lens events can be routinely determined from future high-cadence lensing surveys and thus microlensing can provide a new way to study Galactic binaries.

Significant progress has been made in recent years in constraining the density dependence of nuclear symmetry energy using terrestrial nuclear laboratory data. Around and below the nuclear matter saturation density, the experimental constraints start to merge in a relatively narrow region. At supra-saturation densities, there are, however, still large uncertainties. After summarizing the latest experimental constraints on the density dependence of nuclear symmetry energy, we highlight a few recent studies examining imprints of nuclear symmetry energy on the binding energy, energy release during hadron-quark phase transitions as well as the $w$-mode frequency and damping time of gravitational wave emission of neutron stars.

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

We present the first example of binary microlensing for which the parameter measurements can be verified (or contradicted) by future Doppler observations. This test is made possible by a confluence of two relatively unusual circumstances. First, the binary lens is bright enough (I=15.6) to permit Doppler measurements. Second, we measure not only the usual 7 binary-lens parameters, but also the 'microlens parallax' (which yields the binary mass) and two components of the instantaneous orbital velocity. Thus we measure, effectively, 6 'Kepler+1' parameters (two instantaneous positions, two instantaneous velocities, the binary total mass, and the mass ratio). Since Doppler observations of the brighter binary component determine 5 Kepler parameters (period, velocity amplitude, eccentricity, phase, and position of periapsis), while the same spectroscopy yields the mass of the primary, the combined Doppler + microlensing observations would be overconstrained by 6 + (5 + 1) - (7 + 1) = 4 degrees of freedom. This makes possible an extremely strong test of the microlensing solution. We also introduce a uniform microlensing notation for single and binary lenses, we define conventions, summarize all known microlensing degeneracies and extend a set of parameters to describe full Keplerian motion of the binary lenses.

We report the gravitational microlensing discovery of a sub-Saturn mass planet, MOA-2009-BLG-319Lb, orbiting a K or M-dwarf star in the inner Galactic disk or Galactic bulge. The high cadence observations of the MOA-II survey discovered this microlensing event and enabled its identification as a high magnification event approximately 24 hours prior to peak magnification. As a result, the planetary signal at the peak of this light curve was observed by 20 different telescopes, which is the largest number of telescopes to contribute to a planetary discovery to date. The microlensing model for this event indicates a planet-star mass ratio of q = (3.95 +/- 0.02) x 10^{-4} and a separation of d = 0.97537 +/- 0.00007 in units of the Einstein radius. A Bayesian analysis based on the measured Einstein radius crossing time, t_E, and angular Einstein radius, \theta_E, along with a standard Galactic model indicates a host star mass of M_L = 0.38^{+0.34}_{-0.18} M_{Sun} and a planet mass of M_p = 50^{+44}_{-24} M_{Earth}, which is half the mass of Saturn. This analysis also yields a planet-star three-dimensional separation of a = 2.4^{+1.2}_{-0.6} AU and a distance to the planetary system of D_L = 6.1^{+1.1}_{-1.2} kpc. This separation is ~ 2 times the distance of the snow line, a separation similar to most of the other planets discovered by microlensing.

We have observed 7 new transits of the `hot Jupiter' WASP-5b using a 61 cm telescope located in New Zealand, in order to search for transit timing variations (TTVs) which can be induced by additional bodies existing in the system. When combined with other available photometric and radial velocity (RV) data, we find that its transit timings do not match a linear ephemeris; the best fit \chi^2 values is 32.2 with 9 degrees of freedom which corresponds to a confidence level of 99.982 % or 3.7 \sigma. This result indicates that excess variations of transit timings has been observed, due either to unknown systematic effects or possibly to real TTVs. The TTV amplitude is as large as 50 s, and if this is real, it cannot be explained by other effects than that due to an additional body or bodies. From the RV data, we put an upper limit on the RV amplitude caused by the possible secondary body (planet) as 21 m s^{-1}, which corresponds to its mass of 22-70 M_{Earth} over the orbital period ratio of the two planets from 0.2 to 5.0. From the TTVs data, using the numerical simulations, we place more stringent limits down to 2 M_{Earth} near 1:2 and 2:1 mean motion resonances (MMRs) with WASP-5b at the 3 \sigma level, assuming that the two planets are co-planer. We also put an upper limit on excess of Trojan mass as 43 M_{Earth} (3 \sigma) using both RV and photometric data. We also find that if the possible secondary planet has non- or a small eccentricity, its orbit would likely be near low-order MMRs. Further follow-up photometric and spectroscopic observations will be required to confirm the reality of the TTV signal, and results such as these will provide important information for the migration mechanisms of planetary systems.

Describing the hyperonic and quark phases of neutron stars with an isospin- and momentum-dependent effective interaction for the baryon octet and the MIT bag model, respectively, and using the Gibbs conditions to construct the mixed phase, we study the energy release due to the hadron-quark phase transition. Moreover, the frequency and damping time of the first axial $w$-mode of gravitational waves are studied for both hyperonic and hybrid stars. We find that the energy release is much more sensitive to the bag constant than the density dependence of the nuclear symmetry energy. Also, the frequency of the $w$-mode is found to be significantly different with or without the hadron-quark phase transition and depends strongly on the value of the bag constant. Effects of the density dependence of the nuclear symmetry energy become, however, important for large values of the bag constant that lead to higher hadron-quark transition densities.

We investigate the cosmological and the local tests of the f(R) theory of modified gravity via the observations of (1) the cosmic expansion and (2) the cosmic structures and via (3) the solar-system experiments. To fit the possible cosmic expansion histories under consideration, for each of them we reconstruct f(R), known as "designer f(R)". We then test the designer f(R) via the cosmic-structure constraints on the metric perturbation ratio Psi/Phi and the effective gravitational coupling G_eff and via the solar-system constraints on the Brans-Dicke theory with the chameleon mechanism. We find that among the designer f(R) models specified by the CPL effective equation of state w_eff, only the model closely mimicking general relativity with a cosmological constant (LambdaCDM) can survive all the tests. Accordingly, these tests rule out the frequently studied "w_eff = -1" designer f(R) models which are distinct in cosmic structures from LambdaCDM. When considering only the cosmological tests, we find that the surviving designer f(R) models, although exist for a variety of w_eff, entail fine-tuning.

We report the result of the analysis of a dramatic repeating gravitational microlensing event OGLE-2009-BLG-092/MOA-2009-BLG-137, for which the light curve is characterized by two distinct peaks with perturbations near both peaks. We find that the event is produced by the passage of the source trajectory over the central perturbation regions associated with the individual components of a wide-separation binary. The event is special in the sense that the second perturbation, occurring $\sim 100$ days after the first, was predicted by the real-time analysis conducted after the first peak, demonstrating that real-time modeling can be routinely done for binary and planetary events. With the data obtained from follow-up observations covering the second peak, we are able to uniquely determine the physical parameters of the lens system. We find that the event occurred on a bulge clump giant and it was produced by a binary lens composed of a K and M-type main-sequence stars. The estimated masses of the binary components are $M_1=0.69 \pm 0.11\ M_\odot$ and $M_2=0.36\pm 0.06\ M_\odot$, respectively, and they are separated in projection by $r_\perp=10.9\pm 1.3\ {\rm AU}$. The measured distance to the lens is $D_{\rm L}=5.6 \pm 0.7\ {\rm kpc}$. We also detect the orbital motion of the lens system.

We study the dependence of the galaxy contents within halos on the halo formation time using two galaxy formation models, one being a semianalytic model utilizing the halo assembly history from a high resolution N-body simulation and the other being a smoothed particle hydrodynamics simulation including radiative cooling, star formation, and energy feedback from galactic winds. We confirm the finding by Gao et al. that at fixed mass, the clustering of halos depends on the halo formation time, especially for low-mass halos. This age dependence of halo clustering makes it desirable to study the correlation between the occupation of galaxies within halos and the halo age. We find that, in halos of fixed mass, the number of satellite galaxies has a strong dependence on halo age, with fewer satellites in older halos. The youngest one-third of the halos can have an order of magnitude more satellites than the oldest one-third. For central galaxies, in halos that form earlier, they tend to have more stars and thus appear to be more luminous, and the dependence of their luminosity on halo age is not as strong as that of stellar mass. The results can be understood through the star formation history in halos and the merging of satellites onto central galaxies. The age dependence of the galaxy contents within halos would constitute an important ingredient in a more accurate halo-based model of galaxy clustering.

We report the result of the analysis of the light curve of the microlensing event MOA-2009-BLG-016. The light curve is characterized by a short-duration anomaly near the peak and an overall asymmetry. We find that the peak anomaly is due to a binary companion to the primary lens and the asymmetry of the light curve is explained by the parallax effect caused by the acceleration of the observer over the course of the event due to the orbital motion of the Earth around the Sun. In addition, we detect evidence for the effect of the finite size of the source near the peak of the event, which allows us to measure the angular Einstein radius of the lens system. The Einstein radius combined with the microlens parallax allows us to determine the total mass of the lens and the distance to the lens. We identify three distinct classes of degenerate solutions for the binary lens parameters, where two are manifestations of the previously identified degeneracies of close/wide binaries and positive/negative impact parameters, while the third class is caused by the symmetric cycloid shape of the caustic. We find that, for the best-fit solution, the estimated mass of the lower-mass component of the binary is (0.04 +- 0.01) M_sun, implying a brown-dwarf companion. However, there exists a solution that is worse only by \Delta\chi^2 ~ 3 for which the mass of the secondary is above the hydrogen-burning limit. Unfortunately, resolving these two degenerate solutions will be difficult as the relative lens-source proper motions for both are similar and small (~ 1 mas/yr) and thus the lens will remain blended with the source for the next several decades.

We present a new analysis of the Jupiter+Saturn analog system, OGLE-2006-BLG-109Lb,c, which was the first double planet system discovered with the gravitational microlensing method. This is the only multi-planet system discovered by any method with measured masses for the star and both planets. In addition to the signatures of two planets, this event also exhibits a microlensing parallax signature and finite source effects that provide a direct measure of the masses of the star and planets, and the expected brightness of the host star is confirmed by Keck AO imaging, yielding masses of M_* = 0.51(+0.05-0.04) M_sun, M_b = 231+-19 M_earth, M_c = 86+-7 M_earth. The Saturn-analog planet in this system had a planetary light curve deviation that lasted for 11 days, and as a result, the effects of the orbital motion are visible in the microlensing light curve. We find that four of the six orbital parameters are tightly constrained and that a fifth parameter, the orbital acceleration, is weakly constrained. No orbital information is available for the Jupiter-analog planet, but its presence helps to constrain the orbital motion of the Saturn-analog planet. Assuming co-planar orbits, we find an orbital eccentricity of eccentricity = 0.15 (+0.17-0.10) and an orbital inclination of i = 64 (+4-7) deg. The 95% confidence level lower limit on the inclination of i > 49 deg. implies that this planetary system can be detected and studied via radial velocity measurements using a telescope of >30m aperture.

Gravitational microlensing is not only a successful tool for discovering distant exoplanets, but it also enables characterization of the lens and source stars involved in the lensing event. In high magnification events, the lens caustic may cross over the source disk, which allows a determination of the angular size of the source and additionally a measurement of its limb darkening. When such extended-source effects appear close to maximum magnification, the resulting light curve differs from the characteristic Paczynski point-source curve. The exact shape of the light curve close to the peak depends on the limb darkening of the source. Dense photometric coverage permits measurement of the respective limb-darkening coefficients. In the case of microlensing event OGLE 2008-BLG-290, the K giant source star reached a peak magnification of about 100. Thirteen different telescopes have covered this event in eight different photometric bands. Subsequent light-curve analysis yielded measurements of linear limb-darkening coefficients of the source in six photometric bands. The best-measured coefficients lead to an estimate of the source effective temperature of about 4700 +100-200 K. However, the photometric estimate from colour-magnitude diagrams favours a cooler temperature of 4200 +-100 K. As the limb-darkening measurements, at least in the CTIO/SMARTS2 V and I bands, are among the most accurate obtained, the above disagreement needs to be understood. A solution is proposed, which may apply to previous events where such a discrepancy also appeared.

The flux and line shape of the fine-structure transitions of \NeII\ and \NeIII\ at 12.8 and 15.55\,$\mu$m and of the forbidden transitions of \OI\ $\lambda6300$ are calculated for young stellar objects with a range of mass-loss rates and X-ray luminosities using the X-wind model of jets and the associated wide-angle winds. For moderate and high accretion rates, the calculated \NeII\ line luminosity is comparable to or much larger than produced in X-ray irradiated disk models. All of the line luminosities correlate well with the main parameter in the X-wind model, the mass-loss rate, and also with the assumed X-ray luminosity --- and with one another. The line shapes of an approaching jet are broad and have strong blue-shifted peaks near the effective terminal velocity of the jet. They serve as a characteristic and testable aspect of jet production of the neon fine-structure lines and the \OI\ forbidden transitions.

The evolution of the galaxy stellar mass function is especially useful to test the current model of galaxy formation. Observational data have revealed a few inconsistencies with predictions from the $\Lambda {\rm CDM}$ model. For example, most massive galaxies have already been observed at very high redshifts, and they have experienced only mild evolution since then. In conflict with this, semi-analytical models of galaxy formation predict an insufficient number of massive galaxies at high redshift and a rapid evolution between redshift 1 and 0 . In addition, there is a strong correlation between star formation rate and stellar mass for star-forming galaxies, which can be roughly reproduced with the model, but with a normalization that is too low at high redshift. Furthermore, the stellar mass density obtained from the integral of the cosmic star formation history is higher than the measured one by a factor of 2. In this paper, we study these issues using a semi-analytical model that includes: 1) cold gas accretion in massive halos at high redshift; 2) tidal stripping of stellar mass from satellite galaxies; and 3) an evolving stellar initial mass function (bottom-light) with a higher gas recycle fraction. Our results show that the combined effects from 1) and 2) can predict sufficiently massive galaxies at high redshifts and reproduce their mild evolution at low redshift, While the combined effects of 1) and 3) can reproduce the correlation between star formation rate and stellar mass for star-forming galaxies across wide range of redshifts. A bottom-light/top-heavy stellar IMF could partly resolve the conflict between the stellar mass density and cosmic star formation history.

We present the discovery of a Neptune-mass planet OGLE-2007-BLG-368Lb with a planet-star mass ratio of q=[9.5 +/- 2.1] x 10^{-5} via gravitational microlensing. The planetary deviation was detected in real-time thanks to the high cadence of the MOA survey, real-time light curve monitoring and intensive follow-up observations. A Bayesian analysis returns the stellar mass and distance at M_l = 0.64_{-0.26}^{+0.21} M_\sun and D_l = 5.9_{-1.4}^{+0.9} kpc, respectively, so the mass and separation of the planet are M_p = 20_{-8}^{+7} M_\oplus and a = 3.3_{-0.8}^{+1.4} AU, respectively. This discovery adds another cold Neptune-mass planet to the planetary sample discovered by microlensing, which now comprise four cold Neptune/Super-Earths, five gas giant planets, and another sub-Saturn mass planet whose nature is unclear. The discovery of these ten cold exoplanets by the microlensing method implies that the mass ratio function of cold exoplanets scales as dN_{\rm pl}/d\log q \propto q^{-0.7 +/- 0.2} with a 95% confidence level upper limit of n < -0.35 (where dN_{\rm pl}/d\log q \propto q^n). As microlensing is most sensitive to planets beyond the snow-line, this implies that Neptune-mass planets are at least three times more common than Jupiters in this region at the 95% confidence level.

To improve the planet detection efficiency, current planetary microlensing experiments are focused on high-magnification events searching for planetary signals near the peak of lensing light curves. However, it is known that central perturbations can also be produced by binary companions and thus it is important to distinguish planetary signals from those induced by binary companions. In this paper, we analyze the light curves of microlensing events OGLE-2007-BLG-137/MOA-2007-BLG-091, OGLE-2007-BLG-355/MOA-2007-BLG-278, and MOA-2007-BLG-199/OGLE-2007-BLG-419, for all of which exhibit short-term perturbations near the peaks of the light curves. From detailed modeling of the light curves, we find that the perturbations of the events are caused by binary companions rather than planets. From close examination of the light curves combined with the underlying physical geometry of the lens system obtained from modeling, we find that the short time-scale caustic-crossing feature occurring at a low or a moderate base magnification with an additional secondary perturbation is a typical feature of binary-lens events and thus can be used for the discrimination between the binary and planetary interpretations.

Clusters, filaments, sheets and voids are the building blocks of the cosmic web. In this study, we present and compare two distinct algorithms for finding cosmic filaments and sheets, a task which is far less well established than the identification of dark matter halos or voids. One method is based on the smoothed dark matter density field, the other uses the halo distributions directly. We apply both techniques to one high resolution N-body simulation and reconstruct the filamentary/sheet like network of the dark matter density field. We focus on investigating the properties of the dark matter halos inside these structures, in particular on the directions of their spins and the orientation of their shapes with respect to the directions of the filaments and sheets. We find that both the spin and the major axes of filament-halos with masses <= 10^{13} M_sun/h are preferentially aligned with the direction of the filaments. The spins and major axes of halos in sheets tend to lie parallel to the sheets. There is an opposite mass dependence of the alignment strengths for the spin (negative) and major (positive) axes, i.e. with increasing halo mass the major axis tends to be more strongly aligned with the direction of the filament whereas the alignment between halo spin and filament becomes weaker with increasing halo mass. The alignment strengths as a function of distance to the most massive node halo indicate that there is a transit large scale environment impact: from the 2-D collapse phase of the filament to the 3-D collapse phase of the cluster/node halo at small separation. Overall, the two algorithms for filament/sheet identification investigated here agree well with each other. The method based on halos alone can be easily adapted for use with observational data sets.

We analyze OGLE-2007-BLG-050, a high magnification microlensing event (A ~ 432) whose peak occurred on 2 May, 2007, with pronounced finite-source and parallax effects. We compute planet detection efficiencies for this event in order to determine its sensitivity to the presence of planets around the lens star. Both finite-source and parallax effects permit a measurement of the angular Einstein radius \theta_E = 0.48 +/- 0.01 mas and the parallax \pi_E = 0.12 +/- 0.03, leading to an estimate of the lens mass M = 0.50 +/- 0.14 M_Sun and its distance to the observer D_L = 5.5 +/- 0.4 kpc. This is only the second determination of a reasonably precise (<30%) mass estimate for an isolated unseen object, using any method. This allows us to calculate the planetary detection efficiency in physical units (r_\perp, m_p), where r_\perp is the projected planet-star separation and m_p is the planet mass. When computing planet detection efficiency, we did not find any planetary signature and our detection efficiency results reveal significant sensitivity to Neptune-mass planets, and to a lesser extent Earth-mass planets in some configurations. Indeed, Jupiter and Neptune-mass planets are excluded with a high confidence for a large projected separation range between the planet and the lens star, respectively [0.6 - 10] and [1.4 - 4] AU, and Earth-mass planets are excluded with a 10% confidence in the lensing zone, i.e. [1.8 - 3.1] AU.

Following our previous work(Jiang et al.(2008)), in which we studied the merger time-scale of galaxies in a high-resolution cosmological hydro/N-body simulation, we investigate the potential influence of uncertainties in the numerical implementation of baryonic physics on the merger time-scale. The simulation used in the previous work suffers from the overcooling problem which causes the central galaxies of large halos too massive. This may result in a shortened merger time-scale compared to that in the real universe. We run a similar simulation, but the stellar mass is significantly reduced to model another extreme case of low stellar mass. Our result shows that the merger time-scale is little affected by the star formation recipes, except for the satellites in nearly radial orbits which show a 22 percent higher time-scale in the lower stellar mass case. Since the radial orbits only account for a small part of the satellites' orbits, the fitting formula in Jiang et al.(2008) is still applicable to a reasonable accuracy, nearly immune to the uncertainty in the baryonic physics.

We report the detection of the cool, Jovian-mass planet MOA-2007-BLG-400Lb. The planet was detected in a high-magnification microlensing event (with peak magnification A_max = 628) in which the primary lens transited the source, resulting in a dramatic smoothing of the peak of the event. The angular extent of the region of perturbation due to the planet is significantly smaller than the angular size of the source, and as a result the planetary signature is also smoothed out by the finite source size. Thus the deviation from a single-lens fit is broad and relatively weak (~ few percent). Nevertheless, we demonstrate that the planetary nature of the deviation can be unambiguously ascertained from the gross features of the residuals, and detailed analysis yields a fairly precise planet/star mass ratio of q = 0.0026+/-0.0004, in accord with the large significance (\Delta\chi^2=1070) of the detection. The planet/star projected separation is subject to a strong close/wide degeneracy, leading to two indistinguishable solutions that differ in separation by a factor of ~8.5. Upper limits on flux from the lens constrain its mass to be M < 0.75 M_Sun (assuming it is a main-sequence star). A Bayesian analysis that includes all available observational constraints indicates a primary in the Galactic bulge with a mass of ~0.2-0.5 M_Sun and thus a planet mass of ~ 0.5-1.3 M_Jupiter. The separation and equilibrium temperature are ~0.6-1.1AU (~5.3-9.7AU) and ~103K (~34K) for the close (wide) solution. If the primary is a main-sequence star, follow-up observations would enable the detection of its light and so a measurement of its mass and distance.

We report the detection of an extrasolar planet of mass ratio q ~ 2 x 10^(-4) in microlensing event MOA-2007-BLG-192. The best fit microlensing model shows both the microlensing parallax and finite source effects, and these can be combined to obtain the lens masses of M = 0.060 (+0.028 -0.021) M_sun for the primary and m = 3.3 (+4.9 -1.6) M_earth for the planet. However, the observational coverage of the planetary deviation is sparse and incomplete, and the radius of the source was estimated without the benefit of a source star color measurement. As a result, the 2-sigma limits on the mass ratio and finite source measurements are weak. Nevertheless, the microlensing parallax signal clearly favors a sub-stellar mass planetary host, and the measurement of finite source effects in the light curve supports this conclusion. Adaptive optics images taken with the Very Large Telescope (VLT) NACO instrument are consistent with a lens star that is either a brown dwarf or a star at the bottom of the main sequence. Follow-up VLT and/or Hubble Space Telescope (HST) observations will either confirm that the primary is a brown dwarf or detect the low-mass lens star and enable a precise determination of its mass. In either case, the lens star, MOA-2007-BLG-192L, is the lowest mass primary known to have a companion with a planetary mass ratio, and the planet, MOA-2007-BLG-192Lb, is probably the lowest mass exoplanet found to date, aside from the lowest mass pulsar planet.

Searches for extrasolar planets have uncovered an astonishing diversity of planetary systems, yet the frequency of solar system analogs remains unknown. The gravitational microlensing planet search method is potentially sensitive to multiple-planet systems containing analogs of all the solar system planets except Mercury. We report the detection of a multiple-planet system with microlensing. We identify two planets with masses of ~0.71 and ~0.27 times the mass of Jupiter and orbital separations of ~2.3 and ~4.6 astronomical units orbiting a primary star of mass ~0.50 solar masses at a distance of ~1.5 kiloparsecs. This system resembles a scaled version of our solar system in that the mass ratio, separation ratio, and equilibrium temperatures of the planets are similar to those of Jupiter and Saturn. These planets could not have been detected with other techniques; their discovery from only six confirmed microlensing planet detections suggests that solar system analogs may be common.

We study galaxy mergers using a high-resolution cosmological hydro/N-body simulation with star formation, and compare the measured merger timescales with theoretical predictions based on the Chandrasekhar formula. In contrast to Navarro et al., our numerical results indicate, that the commonly used equation for the merger timescale given by Lacey and Cole, systematically underestimates the merger timescales for minor mergers and overestimates those for major mergers. This behavior is partly explained by the poor performance of their expression for the Coulomb logarithm, \ln (m_pri/m_sat). The two alternative forms \ln (1+m_pri/m_sat) and 1/2\ln [1+(m_pri/m_sat)^2] for the Coulomb logarithm can account for the mass dependence of merger timescale successfully, but both of them underestimate the merger time scale by a factor 2. Since \ln (1+m_pri/m_sat) represents the mass dependence slightly better we adopt this expression for the Coulomb logarithm. Furthermore, we find that the dependence of the merger timescale on the circularity parameter \epsilon is much weaker than the widely adopted power-law \epsilon^{0.78}, whereas 0.94*{\epsilon}^{0.60}+0.60 provides a good match to the data. Based on these findings, we present an accurate and convenient fitting formula for the merger timescale of galaxies in cold dark matter models.

We use high-resolution N-body numerical simulations to study the number of predicted large-separation multiply-imaged systems produced by clusters of galaxies in the SDSS photometric and spectroscopic quasar samples. We incorporate the condensation of baryons at the centre of clusters by (artificially) adding a brightest central galaxy (BCG) as a truncated isothermal sphere. We make predictions in two flat cosmological models: a LCDM model with a matter density $\Omega_0=0.3$, and $\sigma_8=0.9$ (LCDM0), and a model favoured by the WMAP three-year data with $\Omega_0=0.238$, and $\sigma_8=0.74$ (WMAP3). We found that the predicted multiply-imaged quasars with separation >10" is about 6.2 and 2.6 for the SDSS photometric (with an effective area 8000 deg$^2$) and spectroscopic (with an effective area 5000 deg$^2$) quasar samples respectively in the LCDM0 model; the predicted numbers of large-separation lensed quasars agree well with the observations. These numbers are reduced by a factor of 7 or more in the WMAP3 model, and are consistent with data at < 8% level. The predicted cluster lens redshift peaks around redshift 0.5, and 90% are between 0.3 and 1. We find that the BCG creates a central circular region, comparable to the Einstein ring of the BCG, where the central image disappears in the usual three-image and five-image configurations. If we include four image systems as an extreme case of five-image systems (with an infinitely demagnified central image), we find that 68% of the central images are fainter by a factor of 100 than the brightest image, and about 80% are within 1.5"of the BCG.

(abridged) We study the mass assembly history (MAH) of dark matter haloes. We compare MAHs obtained using (i) merger trees constructed with the extended Press-Schechter (EPS) formalism, (ii) numerical simulations, and (iii) the Lagrangian perturbation code PINOCCHIO. We show that the PINOCCHIO MAHs are in excellent agreement with those obtained using numerical simulations. Using a suite of 55 PINOCCHIO simulations, with 256^3 particles each, we study the MAHs of 12,924 cold dark matter haloes in a \LambdaCDM concordance cosmology. We show that haloes less massive than the characteristic non-linear mass scale establish their potential wells much before they acquire most of their mass. The time when a halo reaches its maximum virial velocity roughly divides its mass assembly into two phases, a fast accretion phase which is dominated by major mergers, and a slow accretion phase dominated by minor mergers. Each halo experiences about 3 \pm 2 major mergers since its main progenitor had a mass equal to one percent of the final halo mass. This major merger statistic is found to be virtually independent of halo mass. However, the average redshift at which these major mergers occur, is strongly mass dependent, with more massive haloes experiencing their major mergers later.

Using a set of high-resolution N-body/SPH cosmological simulations with identical initial conditions but run with different numerical setups, we investigate the influence of baryonic matter on the mass distribution of dark halos when radiative cooling is NOT included. We compare the concentration parameters of about 400 massive halos with virial mass from $10^{13}$ \Msun to $7.1 \times 10^{14}$ \Msun. We find that the concentration parameters for the total mass and dark matter distributions in non radiative simulations are on average larger by ~3% and 10% than those in a pure dark matter simulation. Our results indicate that the total mass density profile is little affected by a hot gas component in the simulations. After carefully excluding the effects of resolutions and spurious two-body heating between dark matter and gas particles, we conclude that the increase of the dark matter concentration parameters is due to interactions between baryons and dark matter. We demonstrate this with the aid of idealized simulations of two-body mergers. The results of individual halos simulated with different mass resolutions show that the gas profiles of densities, temperature and entropy are subjects of mass resolution of SPH particles. In particular, we find that in the inner parts of halos, as the SPH resolution increases the gas density becomes higher but both the entropy and temperature decrease.

We use high-resolution $N$-body simulations to investigate the optical depth of giant arcs with length-to-width ratio larger than 7.5 and 10 in the `standard' $\LCDM$ model with $\sigma_8=0.9$ and $\Omega_{\rm m,0}=0.3$ and a model based on three-year Wilkinson Microwave Anisotropy Probe (WMAP) data. We find that, in dark-matter only simulations, the lensing probability in the three-year WMAP model (with $\sigma_8=0.74$ and $\Omega_{\rm m,0}=0.238)$ decreases by a factor of $\sim 6$ compared with that in the `standard' $\LCDM$ model. The effects of baryonic cooling, star formation and feedbacks are uncertain, but we argue that baryons will only increase the the lensing cross-section by a moderate factor, $\sim 2$. We conclude that the low central value of $\sigma_8$ and $\Omega_{\rm m,0}$ preferred by the WMAP three-year data may be too low to be compatible with observations if conventional assumptions of the background source population are correct.

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 present 40 fully hydrodynamical numerical simulations of the intergalactic gas that gives rise to the Ly-alpha forest. We make artificial spectra from each and measure three output parameters: the mean flux, a measure of the most common Ly-alpha line width (b-value), and the 1D power spectrum of the flux on scales from 0.01 - 0.1 s/km. We investigate how these output parameters are correlated with the input parameters that describe the simulations. These input parameters include the main cosmological parameters; and two astrophysical parameters gammah and gammahe. The parameter gammah controls the rate of ionization of H I, He I and He II while gammahe controls the rate of heating. We compare the output from the simulations to data from HIRES and UVES spectra and we give a new measurement of the flux power from these spectra at z=1.95. We find that simulations with a wide variety of sigma_8 values, from at least 0.8 - 1.1, can fit the small scale flux power and b-values when we adjust gammahe to compensate for the sigma_8 change. We can also use gammah to adjust the H I ionization rate to simultaneously match the mean flux. We derive scaling relations that give the output parameter values expected for a variety of input parameters with errors of 2% for optical depth and 1.4% for b-balues. We can break the degeneracy between gammahe and sigma_8 using large scale power or other data to fix sigma_8.

We identify a concordant model for the intergalactic medium (IGM) at redshift z=1.9 that uses popular values for cosmological and astrophysical parameters and accounts for all baryons with an uncertainty of 6%. We have measured the amount of absorption, DA, in the Ly-alpha forest at redshift 1.9 in spectra of 77 QSO from the Kast spectrograph. We calibrated the continuum fits with realistic artificial spectra, and we found that averaged over all 77 QSOs the mean continuum level is within 1-2% of the correct value. Absorption from all lines in the Ly-alpha forest at z=1.9 removes DA=15.1 +/- 0.7% of the flux between 1070 and 1170 (rest) Angstroms. This is the first measurement using many QSOs at this z, and the first calibrated measurement at any redshift. Metal lines absorb 2.3 +/- 0.5%, and LLS absorb 1.0 +/- 0.4% leaving 11.8 +/- 1.0% from the lower density bulk of the IGM. Averaging over Delta z=0.1 or 154 Mpc, the dispersion is 6.1 +/- 0.3% including LLS and metal lines, or 3.9 (+0.5, -0.7)% for the lower density IGM alone, consistent with the usual description of large scale structure. LLS and metal lines are major contributors to the variation in the mean flux, and they make the flux field significantly non-Gaussian. We find that a hydrodynamic simulation on a 1024 cubed grid in a 75.7 Mpc box reproduces the observed DA from the low density IGM with parameters values H_o=71 km/s/Mpc, Omega_Lambda=0.73, Omega_m=0.27, Omega_b=0.044, sigma_8=0.9 and a UV background that has an ionization rate that is 1.08 +/- 0.08 times the prediction by Madau, Haardt & Rees (1999).

This paper uses numerical simulations to test the formation time distribution of dark matter haloes predicted by the analytic excursion set approaches. The formation time distribution is closely linked to the conditional mass function and this test is therefore an indirect probe of this distribution. The excursion set models tested are the extended Press-Schechter (EPS) model, the ellipsoidal collapse (EC) model, and the non-spherical collapse boundary (NCB) model. Three sets of simulations (6 realizations) have been used to investigate the halo formation time distribution for halo masses ranging from dwarf-galaxy like haloes ($M=10^{-3} M_*$, where $M_*$ is the characteristic non-linear mass scale) to massive haloes of $M=8.7 M_*$. None of the models can match the simulation results at both high and low redshift. In particular, dark matter haloes formed generally earlier in our simulations than predicted by the EPS model. This discrepancy might help explain why semi-analytic models of galaxy formation, based on EPS merger trees, under-predict the number of high redshift galaxies compared with recent observations.

In a first of a series of studies of the H-alpha + [NII] emission from nearby spiral galaxies, we present measurements of H-alpha + [NII] emission from HII regions in M81. Our method uses large-field-CCD images and long-slit spectra, and is part of the ongoing Beijing-Arizona-Taipei-Connecticut Sky Survey. The CCD images are taken with the NAOC 0.6/0.9m f/3 Schmidt telescope at the Xinglong Observing Station, using a multicolor filter set. Spectra of 10 of the brightest HII regions are obtained using the NAOC 2.16m telescope with a Tek 1024 X 1024 CCD. The continua of the spectra are calibrated by flux-calibrated images taken from the Schmidt observations. We determine the continuum component of our H-alpha + [NII] image via interpolation from the more accurately-measured backgrounds (M81 starlight) obtained from the two neighboring (in wavelength) BATC filter images. We use the calibrated fluxes of H-alpha + [NII] emission from the spectra to normalize this interpolated, continuum-subtracted H-alpha + [NII] image. We estimate the zero point uncertainty of the measured H-alpha + [NII] emission flux to be $\sim$ 8%. A catalogue of H-alpha + [NII] fluxes for 456 HII regions is provided, with those fluxes are on a more consistent linear scale than previously available. The logarithmically-binned H-alpha + [NII] luminosity function of HII regions is found to have slope $\alpha$ = -0.70, consistent with previous results (which allowed $\alpha=-0.5 \sim -0.8$). From the overall H-alpha + [NII] luminosity of the HII regions, the star formation rate of M81 is found to be $\sim 0.68 M_{\odot} {\rm yr}^{-1}$, modulo uncertainty with extinction corrections.

Measurements of the Hubble constant to distant galaxy clusters using the Sunyaev-Zeldovich effect are systematically low in comparison to values obtained by other means. These measurements usually assume a spherical isothermal $\beta$ model for the intracluster medium (ICM). We present the results of a statistical analysis of {\em temperature bias} in $H_0$ determinations in a sample of 27 numerically simulated X-ray clusters drawn from a $\Lambda$CDM model at z=0.5. We employ adaptive mesh refinement to provide high resolution (15.6 h$^{-1} kpc$) in cluster cores which dominate the X-ray and radio signals. Fitting synthetic X-ray and y-parameter maps to the standard isothermal beta model, we find a broad, skewed distribution in $f\equiv H_0(SZ)/H_0(actual)$ with a mean, median, and standard deviation of 0.89, 0.83 and 0.32 respectively, where $H_0(SZ)$ is the value of $H_0$ derived by using Sunyaev-Zeldovich effect method and $H_0(actual)$ is the value used in the cosmological simulation. We find that the clusters' declining temperature profiles systematically lower estimates of $H_0$ by 10% to 20%. The declining temperature profile of our adiabatic system is consistent with the result including radiative cooling and supernovae feedback (Loken et al. 2002). We thereby introduce a non-isothermal $\beta$ model as an improvement. Applying the non-isothermal $\beta$ model to the refined sample with well-fitted temperature profile, the value of f improves 9% relative to the actual value. The study of the morphology and the clumping effects conclude that these two factors combine to induce scatter in f of $\pm$ 30 %.

In recent years, the possibility of measuring the cosmological constant $\Omega_\Lambda$ through the application of the Alcock-Paczynski test to the Lyman Alpha (Ly$\alpha$) forest has been suggested (McDonald et al. 1999; Hui et al. 1999). Despite the theoretical uncertainties due to a few other cosmological parameters, some of the greatest difficulties we encounter concern the huge uncertainties due to cosmic variance and noise. In this paper, we propose a maximum likelihood estimation (MLE) method to deal with cosmic variance and noise using synthetic spectra of quasistellar objects (QSOs) from our cosmological hydrodynamic simulations. We demonstrate that the MLE method can overcome the cosmic variance problem. Applying the MLE method, we find that we have more than 90% probability to determine $\Omega_\Lambda$ within 20% error and approximately of 66% probability to determine $\Omega_\Lambda$ within 10% error by using 30 pairs QSO spectra when other cosmological parameters are assumed. Another important source of error is from noise in the flux spectra, and we have modeled the corresponding effect by studying artificial spectra with different kinds of noise added. We discover that the noise distribution does not have significant effect on the final cross-correlation functions as long as the signal-to-noise ratio (S/N) is fixed. Finally, a preliminary test and discussion about the sensitivities to other cosmological parameters are included in this paper as well.

We study the effects of gravitational back-reaction in models of Quintessence. The effective energy-momentum tensor with which cosmological fluctuations back-react on the background metric will in some cases lead to a termination of the phase of acceleration. The fluctuations we make use of are the perturbations in our present Universe. Their amplitude is normalized by recent measurements of anisotropies in the cosmic microwave background, their slope is taken to be either scale-invariant, or characterized by a slightly blue tilt. In the latter case, we find that the back-reaction effect of fluctuations whose present wavelength is smaller than the Hubble radius but which are stretched beyond the Hubble radius by the accelerated expansion during the era of Quintessence domination can become large. Since the back-reaction effects of these modes oppose the acceleration, back-reaction will lead to a truncation of the period of Quintessence domination. This result impacts on the recent discussions of the potential incompatibility between string theory and Quintessence.

It is shown that the power spectrum defined in the Synchronous Gauge can not be directly used to calculate the predictions of cosmological models on the large-scale structure of universe, which should be calculated directly by a suitable gauge-invariant power spectrum or the power spectrum defined in the Newtonian Gauge.

There is increasing evidence that conventional cold dark matter (CDM) models lead to conflicts between observations and numerical simulations of dark matter halos on sub-galactic scales. Spergel and Steinhardt showed that if the CDM is strongly self-interacting, then the conflicts disappear. However, the assumption of strong self-interaction would rule out the favored candidates for CDM, namely weakly interacting massive particles (WIMPs), such as the neutralino. In this paper we propose a mechanism of non-thermal production of WIMPs and study its implications on the power spectrum. We find that the non-vanishing velocity of the WIMPs suppresses the power spectrum on small scales compared to what it obtained in the conventional CDM model. Our results show that, in this context, WIMPs as candidates for dark matter can work well both on large scales and on sub-galactic scales.

In this paper we present Monte-Carlo simulations of Lyman alpha absorption systems which originate in galactic haloes, galaxy discs and dark matter (DM) satellites around big central haloes. It is found that for strong Lyman alpha absorption lines galactic haloes and satellites can explain ~20% and 40% of the line number density of QSO absorption line key project respectively. If big galaxies indeed possess such large numbers of DM satellites and they possess gas, these satellites may play an important role for strong Lyman alpha lines. However the predicted number density of Lyman-limit systems by satellites is \~0.1 (per unit redshift), which is four times smaller than that by halo clouds. Including galactic haloes, satellites and HI discs of spirals, the predicted number density of strong lines can be as much as 60% of the HST result. The models can also predict all of the observed Lyman-limit systems. The average covering factor within 250 kpc/h is estimated to be ~0.36. And the effective absorption radius of a galaxy is estimated to be ~150 kpc/h. The models predict W_r propto rho^{-0.5} L_B^{0.15} (1+z)^{-0.5}. We study the selection effects of selection criteria similar to the imaging and spectroscopic surveys. We simulate mock observations through known QSO lines-of-sight and find that selection effects can statistically tighten the dependence of line width on projected distance. (abridged)

We re-examine approximations in the analytical calculation of the primordial spectrum of cosmological perturbation produced during inflation. Taking two inflation models (chaotic inflation and natural inflation) as examples, we numerically verify the accuracy of these approximations.

We calculate the thermal flux of low-energy solar neutrinos and antineutrinos of all flavors arising from a variety of neutrino pair processes: Compton production (including plasmon-pole diagrams), neutral current decay of thermally populated nuclear states, plasmon decay, and electron transitions from free to atomic bound states. The resulting flux density per flavor is significant (10E8-10E9/cm2/sec/MeV) below about 5 keV, and the distributions fill much of the valley between the high-energy edge of the cosmic background neutrino spectrum and the low energy tails of the pp-chain electron neutrino and terrestrial electron antineutrino spectra. Thermal neutrinos carry information on the solar core temperature distribution and on heavy flavor masses in the range of 1 keV. The detection of these neutrinos is a daunting but interesting challenge.

We study the possibilities of parametric amplification of the gravitational perturbation during reheating in single-field inflation models. Our result shows that there is no additional growth of the super-horizon modes beyond the usual predictions.

The edge-on, nearby spiral galaxy NGC5907 has long been used as the prototype of a ``non-interacting'' warped galaxy. We report here the discovery of two interactions with companion dwarf galaxies that substantially change this picture. First, a faint ring structure is discovered around this galaxy that is likely due to the tidal disruption of a companion dwarf spheroidal galaxy. The ring is elliptical in shape with the center of NGC5907 close to one of the ring's foci. This suggests the ring material is in orbit around NGC5907. No gaseous component to the ring has been detected either with deep H$\alpha$ images or in Very Large Array (VLA) HI 21-cm line maps. The visible material in the ring has an integrated luminosity <= 10^8 L_sun and its brightest part has a color R-I ~ 0.9. All of these properties are consistent with the ring being a tidally-disrupted dwarf spheroidal galaxy. Second, we find that NGC5907 has a dwarf companion galaxy, PGC54419, projected to be only 36.9 kpc from the center of NGC5907, close in radial velocity (\Delta V = 45 km s^{-1}) to the giant spiral galaxy. This dwarf is seen at the tip of the HI warp and in the direction of the warp. Hence, NGC5907 can no longer be considered ``non-interacting,'' but is obviously interacting with its dwarf companions much as the Milky Way interacts with its dwarf galaxies. These results, coupled with the finding by others that dwarf galaxies tend to be found around giant galaxies, suggest that tidal interaction with companions, even if containing a mere 1% of the mass of the parent galaxy, might be sufficient to excite the warps found in the disks of many large spiral galaxies.

Core and conal emission beams of radio pulsars have been identified observationally (Rankin 1983,1993; Lyne and Manchester 1988), In the inner gap model (Ruderman and Sutherland 1975, hereafter RS75), the gap continually breaks down (sparking) by forming electron-positron pairs on a time scale of a few microseconds (RS75). This makes a large amplitude low frequency wave to be produced and which would be scattered by relativistic particles moving out from the gap. Under this assumption, Qiao (1988a,1992) presented an Inverse Compton Scattering (ICS) model for both core and conal emissions. This paper presented a development of the model. Retardation and aberration effects due to different radiation components emitted from different heights are considered. The luminosity of pulsar radio emission via the ICS process is discussed. Coherent emission by bunches of particles is adopted and which is adequate to explain pulsar radiation. The theoretical results are in agreement with observations very well.