Providing a possible connection between neutrino emission and gravitational-wave (GW) bursts is important to our understanding of the physical processes that occur when black holes or neutron stars merge. In the Daya Bay experiment, using data collected from December 2011 to August 2017, a search has been performed for electron-antineutrino signals coinciding with detected GW events, including GW150914, GW151012, GW151226, GW170104, GW170608, GW170814, and GW170817. We used three time windows of $\mathrm{\pm 10~s}$, $\mathrm{\pm 500~s}$, and $\mathrm{\pm 1000~s}$ relative to the occurrence of the GW events, and a neutrino energy range of 1.8 to 100 MeV to search for correlated neutrino candidates. The detected electron-antineutrino candidates are consistent with the expected background rates for all the three time windows. Assuming monochromatic spectra, we found upper limits (90% confidence level) on electron-antineutrino fluence of $(1.13~-~2.44) \times 10^{11}~\rm{cm^{-2}}$ at 5 MeV to $8.0 \times 10^{7}~\rm{cm^{-2}}$ at 100 MeV for the three time windows. Under the assumption of a Fermi-Dirac spectrum, the upper limits were found to be $(5.4~-~7.0)\times 10^{9}~\rm{cm^{-2}}$ for the three time windows.

Strong gravitational lensing has been a powerful probe of cosmological models and gravity. To date, constraints in either domain have been obtained separately. We propose a new methodology through which the cosmological model, specifically the Hubble constant, and post-Newtonian parameter can be simultaneously constrained. Using the time-delay cosmography from strong lensing combined with the stellar kinematics of the deflector lens, we demonstrate the Hubble constant and post-Newtonian parameter are incorporated in two distance ratios which reflect the lensing mass and dynamical mass, respectively. Through the reanalysis of the four publicly released lenses distance posteriors from the H0LiCOW collaboration, the simultaneous constraints of Hubble constant and post-Newtonian parameter are obtained. Our results suggests no deviation from the General Relativity, $\gamma_{\texttt{PPN}}=0.87^{+0.19}_{-0.17}$ with a Hubble constant favors the local Universe value, $H_0=73.65^{+1.95}_{-2.26}$ km s$^{-1}$ Mpc$^{-1}$. Finally, we forecast the robustness of gravity tests by using the time-delay strong lensing for constraints we expect in the next few years. We find that the joint constraint from 40 lenses are able to reach the order of $7.7\%$ for the post-Newtonian parameter and $1.4\%$ for Hubble constant.

We examine the light curves of two quasars, motivated by recent suggestions that a supermassive black hole binary (SMBHB) can exhibit sharp lensing spikes. We model the variability of each light curve as due to a combination of two relativistic effects: the orbital relativistic Doppler boost and gravitational binary self-lensing. In order to model each system we extend previous Doppler plus self-lensing models to include eccentricity. The first quasar is identified in optical data as a binary candidate with a 20-yr period (Ark 120), and shows a prominent spike. For this source, we rule out the lensing hypothesis and disfavor the Doppler-boost hypothesis due to discrepancies in the measured vs. recovered values of the binary mass and optical spectral slope. The second source, which we nickname Spikey, is the rare case of an active galactic nucleus (AGN) identified in Kepler's high-quality, high-cadence photometric data. For this source, we find a model, consisting of a combination of a Doppler modulation and a narrow symmetric lensing spike, consistent with an eccentric SMBHB with a total mass of approximately 30 million solar masses, rest-frame orbital period T=418 days, eccentricity e=0.5, and seen at an inclination of 8 degrees from edge-on. This interpretation can be tested by monitoring Spikey for periodic behavior and recurring flares in the next few years. In preparation for such monitoring we present the first X-ray observations of this object taken by the Neil Gehrels Swift observatory.

We present results of continuum and spectral line observations with ALMA and 22 GHz water (H$_2$O) maser observations using KaVA and VERA toward a high-mass star-forming region, G25.82-0.17. Multiple 1.3 mm continuum sources are revealed, indicating the presence of young stellar objects (YSOs) at different evolutionary stages, namely an ultra-compact HII region, G25.82-E, a high-mass young stellar object (HM-YSO), G25.82-W1, and starless cores, G25.82-W2 and G25.82-W3. Two SiO outflows, at N-S and SE-NW orientations, are identified. The CH$_3$OH 8$_{-1}$-7$_{0}$ E line, known to be a class I CH$_3$OH maser at 229 GHz is also detected showing a mixture of thermal and maser emission. Moreover, the H$_2$O masers are distributed in a region ~0.25" shifted from G25.82-W1. The CH$_3$OH 22$_{4}$-21$_{5}$ E line shows a compact ring-like structure at the position of G25.82-W1 with a velocity gradient, indicating a rotating disk or envelope. Assuming Keplerian rotation, the dynamical mass of G25.82-W1 is estimated to be $>$25 M$_{\odot}$ and the total mass of 20 M$_\odot$-84 M$_\odot$ is derived from the 1.3 mm continuum emission. The driving source of the N-S SiO outflow is G25.82-W1 while that of the SE-NW SiO outflow is uncertain. Detection of multiple high-mass starless$/$protostellar cores and candidates without low-mass cores implies that HM-YSOs could form in individual high-mass cores as predicted by the turbulent core accretion model. If this is the case, the high-mass star formation process in G25.82 would be consistent with a scaled-up version of low-mass star formation.

Gravitational waves (GWs) are produced by colliding particles through the gravitational analogue of electromagnetic bremsstrahlung. We calculate the contribution of free-free emission in the radiation-dominated Universe to the stochastic GW background. We find that the energy density of the resulting GW radiation is heavily dependent on the number of elementary particles, $N_{\mathrm{tot}}$, and the maximum initial temperature, $T_{\mathrm{max}}$. We rule out $N_{\mathrm{tot}}\gtrsim N_{\mathrm{SM}}$ for $T_{\mathrm{max}}\sim T_{\mathrm{Planck}}\approx10^{19}$ GeV and $N_{\mathrm{tot}}\gtrsim10^{13}\times N_{\mathrm{SM}}$ for $T_{\mathrm{max}}\sim10^{16}$ GeV, where $N_{\mathrm{SM}}$ is the number of particles in the Standard Model. In the case of inflation, existing cosmological data constrain $T_{\mathrm{max}}\lesssim10^{16}$ GeV. However, alternative models to inflation such as bouncing cosmologies allow for $T_{\mathrm{max}}$ near $T_{\mathrm{Planck}}$. At the energy scales we are considering, the extra number of particles arise naturally in models of extra dimensions.

Strong gravitational lensing by galaxies provides us with a unique opportunity to understand the nature of gravity on galactic and extra-galactic scales. In this paper, we propose a new multimessenger approach using data from both gravitational wave (GW) and the corresponding electromagnetic (EM) counterpart to infer the constraint of the modified gravity (MG) theory denoted by the scale dependent phenomenological parameter. To demonstrate the robustness of this approach, we calculate the time-delay predictions by choosing various values of the phenomenological parameters and then compare them with that from the general relativity (GR). For the third generation ground-based GW observatory, with one typical strongly lensed GW+EM event, and assuming that the dominated error from the stellar velocity dispersions is 5\%, the GW time-delay data can distinguish an 18\% MG effect on a scale of tens of kiloparsecs with a $68\%$ confidence level. Assuming GR and a Singular Isothermal Sphere mass model, there exists a simplified consistency relationship between time-delay and imaging data. This relationship does not require for the velocity dispersion measurement, and hence can avoid major uncertainties. By using this relationship, the multimessenger approach is able to distinguish an $8\%$ MG effect. Our results show that the GW multimessenger approach can play an important role in revealing the nature of gravity on the galactic and extra-galactic scales.

Halo model is a physically intuitive method for modelling the non-linear power spectrum, especially for the alternatives to the standard $\Lambda$CDM models. In this paper, we exam the Sheth-Tormen barrier formula adopted in the previous \texttt{CHAM} method \citep{2018MNRAS.476L..65H}. As an example, we model the ellipsoidal collapse of top-hat dark matter haloes in $f(R)$ gravity. A good agreement between Sheth-Tormen formula and our result is achieved. The relative difference in the ellipsoidal collapse barrier is less than or equal to $1.6\%$. Furthermore, we verify that, for F4 and F5 cases of Hu-Sawicki $f(R)$ gravity, the screening mechanism do not play a crucial role in the non-linear power spectrum modelling up to $k\sim1[h/{\rm Mpc}]$. We compare two versions of modified gravity modelling, namely with/without screening. We find that by treating the effective Newton constant as constant number ($G_{\rm eff}=4/3G_N$) is acceptable. The scale dependence of the gravitational coupling is sub-relevant. The resulting spectra in F4 and F5, are in $0.1\%$ agreement with the previous \texttt{CHAM} results. The published code is accelerated significantly. Finally, we compare our halo model prediction with N-body simulation. We find that the general spectrum profile agree, qualitatively. However, via the halo model approach, there exists a systematic under-estimation of the matter power spectrum in the co-moving wavenumber range between $0.3 h/{\rm Mpc}$ and $3 h/{\rm Mpc}$. These scales are overlapping with the transition scales from two halo term dominated regimes to those of one halo term dominated.

We compile and analyze ~200 trigonometric parallaxes and proper motions of molecular masers associated with very young high-mass stars. These measurements strongly suggest that the Milky Way is a four-arm spiral. Fitting log-periodic spirals to the locations of the masers, allows us to significantly expand our view of the structure of the Milky Way. We present an updated model for its spiral structure and incorporate it into our previously published parallax-based distance-estimation program for sources associated with spiral arms. Modeling the three-dimensional space motions yields estimates of the distance to the Galactic center, Ro = 8.15 +/- 0.15 kpc, the circular rotation speed at the Sun's position, To = 236 +/- 7 km/s, and the nature of the rotation curve. Our data strongly constrain the full circular velocity of the Sun, To + Vsun = 247 +/- 4 km/s, and its angular velocity, (To + Vsun)/Ro = 30.32 +/- 0.27 km/s/kpc. Transforming the measured space motions to a Galactocentric frame which rotates with the Galaxy, we find non-circular velocity components typically about 10 km/s. However, near the Galactic bar and in a portion of the Perseus arm, we find significantly larger non-circular motions. Young high-mass stars within 7 kpc of the Galactic center have a scale height of only 19 pc and, thus, are well suited to define the Galactic plane. We find that the orientation of the plane is consistent with the IAU-defined plane to within +/-0.1 deg., and that the Sun is offset toward the north Galactic pole by Zsun = 5.5 +/- 5.8 pc. Accounting for this offset places the central supermassive black hole, Sgr A*, in the midplane of the Galaxy. Using our improved Galactic parameters, we predict the Hulse-Taylor binary pulsar to be at a distance of 6.54 +/- 0.24 kpc, assuming its orbital decay from gravitational radiation follows general relativity.

We report trigonometric parallax and proper motion measurements of 6.7-GHz CH3OH and 22-GHz H2O masers in eight high-mass star-forming regions (HMSFRs) based on VLBA observations as part of the BeSSeL Survey. The distances of these HMSFRs combined with their Galactic coordinates, radial velocities, and proper motions, allow us to assign them to a segment of the Perseus arm with ~< 70 deg. These HMSFRs are clustered in Galactic longitude from ~30 deg to ~50, neighboring a dirth of such sources between longitudes ~50 deg to ~90 deg.

The cosmic distance duality relation (CDDR) has been test through several astronomical observations in the last years. This relation establishes a simple equation relating the angular diameter ($D_A$) and luminosity ($D_L$) distances at a redshift $z$, $D_LD_A^{-1}(1+z)^{-2}=\eta=1$. However, only very recently this relation has been observationally tested at high redshifts ($z \approx 3.6$) by using luminosity distances from type Ia supernovae (SNe Ia) and gamma ray bursts (GRBs) plus angular diameter distances from strong gravitational lensing (SGL) observations. The results show that no significant deviation from the CDDR validity has been verified. In this work, we test the potentialities of future luminosity distances from gravitational waves (GWs) sources to impose limit on possible departures of CDDR jointly with current SGL observations. The basic advantage of $D_L$ from GWs is being insensitive to non-conservation of the number of photons. By simulating 600, 900 and 1200 data of GWs using the Einstein Telescope (ET) as reference, we derive limits on $\eta(z)$ function and obtain that the results will be at least competitive with current limits from the SNe Ia $+$ GRBs $+$ SGLs analyses.

Using infrared to (sub)millimeter data from Spitzer, Herschel, the James Clerk Maxwell Telescope, and the IRAM 30-m telescope, we conducted an unbiased survey of the massive dense cores (MDCs) in the Cygnus X molecular cloud complex, aimed at characterizing the physical conditions of high-mass star formation (HMSF) at ~0.1 pc scales. We created 5 deg times 6 deg images of the 70-1200 um dust continuum, gas column density, and dust temperature of Cygnus X. A spatial relation between the dense regions (Av > 15) and the developed HII regions was found, indicating the impact of the latter on the global structures of Cygnus X. With a 35-Msun mass threshold implied by HMSF signposts, we identified 151 MDCs with sizes of ~0.1 pc, masses of 35-1762 Msun, and temperatures of 8-35 K. Our MDC sample is statistically complete in Cygnus X and is three times larger than that in Motte et al. (2007). The MDCs were classified into IR-bright/IR-quiet ones based on their mid-infrared fluxes and a large "IR-quiet" proportion (90%) was found in our sample. Two possible scenarios were proposed to interpret this: accelerated HMSF and the incapability of HMSF of the IR-quiet MDCs. We also found 26 starless MDCs by their lack of compact emissions at 21-70 um wavelengths, of which the most massive ones are probably the best candidates of initial HMSF sites in Cygnus X.

We study the co-evolution of supermassive black holes (SMBHs) with galaxies by means of semi-analytic model (SAM) of galaxy formation based on sub-halo merger trees built from Millennium and Millennium-II simulation. We utilize the simulation results from Guo 2013 and Henriques 2015 to study two aspects of the co-evolution, \emph{i.e.} the stochastic gravitational wave (GW) background generated by SMBH merger and the SMBH/galaxy clustering. The characteristic strain amplitude of GW background predicted by Guo 2013 and Henriques 2015 models are $A_{yr^{-1}}=5.00\times10^{-16}$ and $A_{yr^{-1}}=9.42\times10^{-17}$, respectively. We find the GW amplitude is very sensitive to the galaxy merger rate. The difference in the galaxy merger rate between Guo 2013 and Henriques 2015, results in a factor $5$ deviation in the GW strain amplitude. For clusterings, we calculate the spatially isotropic two point auto- and cross-correlation functions (2PCFs) for both SMBHs and galaxies by using the mock catalogs generated from Guo 2013 model. We find that all 2PCFs have positive dependence on both SMBH and galaxy mass. And there exist a significant time evolution in 2PCFs, namely, the clustering effect is enhanced at lower redshifts. Interestingly, this result is not reported in the active galactic nuclei samples in SDSS. Our analysis also shows that, roughly, SMBHs and galaxies, with galaxy mass $10^2\sim10^3$ larger than SMBH mass, have similar pattern of clustering, which is a reflection of the co-evolution of SMBH and galaxy. Finally, we calculate the first ten multiples of the angular power spectrum of the energy density of GW background. We find the amplitude of angular power spectrum of the first ten multiples is about $10\%$ to $60\%$ of the monopole component in the whole frequency range.

The Weak Equivalence Principle states that the geodesics of a test particle in a gravitational field are independent of the particle's constitution. To constrain violations of the Weak Equivalence Principle, we use the one-meter telescope at Table Mountain Observatory near Los Angeles to monitor the relative arrival times of pulses from the Crab Pulsar in the optical ($\lambda \approx 585$ nm) and near-infrared ($\lambda \approx 814$ nm) using an instrument which detects single photons with nanosecond-timing resolution in those two bands. The infrared pulse arrives slightly before the visible pulse. Our three analysis methods give delays with statistical errors of $\Delta t_{obs} = 7.41 \pm 0.58$, $0.4 \pm 3.6$, and $7.35 \pm 4.48$ microseconds (at most 1/4000 of the pulsar period). We attribute this discrepancy to systematic error from the fact that the visible and infrared pulses have slightly different shapes. Whether this delay emerges from the pulsar, is caused by passing through wavelength-dependent media, or is caused by a violation of the equivalence principle, unless there is a fine-tuned cancellation among these, we set the first upper limit on the differential post-Newtonian parameter at these wavelengths of $\Delta \gamma < 1.07 \times 10^{-10}~(3\sigma)$. This result falls in an unexplored region of parameter space and complements existing limits on equivalence-principle violation from fast radio bursts, gamma ray bursts, as well as previous limits from the Crab.

We revisit the cosmology of Covariant Galileon gravity in view of the most recent cosmological data sets, including weak lensing. As a higher derivative theory, Covariant Galileon models do not have a $\Lambda$CDM limit and predict a very different structure formation pattern compared with the standard $\Lambda$CDM scenario. Previous cosmological analyses suggest that this model is marginally disfavoured, yet can not be completely ruled out. In this work we use a more recent and extended combination of data, and we allow for more freedom in the cosmology, by including a massive neutrino sector with three different mass hierarchies. We use the Planck measurements of Cosmic Microwave Background temperature and polarization; Baryonic Acoustic Oscillations measurements by BOSS DR12; local measurements of $H_0$; the joint light-curve analysis supernovae sample; and, for the first time, weak gravitational lensing from the KiDS collaboration. We find, that in order to provide a reasonable fit, a non-zero neutrino mass is indeed necessary, but we do not report any sizable difference among the three neutrino hierarchies. Finally, the comparison of the Bayesian Evidence to the $\Lambda$CDM one shows that in all the cases considered, Covariant Galileon models are statistically ruled out by cosmological data.

We present a fast numerical screened halo model algorithm (CHAM) for modeling non-linear power spectrum for the alternative models to LCDM. This method has three obvious advantages. First of all, it is not being restricted to a specific dark energy/modified gravity model. In principle, all of the screened scalar-tensor theories can be applied. Second, the least assumptions are made in the calculation. Hence, the physical picture is very easily understandable. Third, it is very predictable and does not rely on the calibration from N-body simulation. As an example, we show the case of Hu-Sawicki f(R) gravity. In this case, the typical CPU time with the current parallel Python script (8 threads) is roughly within $10$ minutes. The resulting spectra are in a good agreement with N-body data within a few percentage accuracy up to k~1 h/Mpc.

We compare Einstein-Boltzmann solvers that include modifications to General Relativity and find that, for a wide range of models and parameters, they agree to a high level of precision. We look at three general purpose codes that primarily model general scalar-tensor theories, three codes that model Jordan-Brans-Dicke (JBD) gravity, a code that models f(R) gravity, a code that models covariant Galileons, a code that models Ho\v{r}ava-Lifschitz gravity and two codes that model non-local models of gravity. Comparing predictions of the angular power spectrum of the cosmic microwave background and the power spectrum of dark matter for a suite of different models, we find agreement at the sub-percent level. This means that this suite of Einstein-Boltzmann solvers is now sufficiently accurate for precision constraints on cosmological and gravitational parameters.

We update the search for features in the Cosmic Microwave Background (CMB) power spectrum due to transient reductions in the speed of sound, using Planck 2015 CMB temperature and polarisation data. We enlarge the parameter space to much higher oscillatory frequencies of the feature, and define a robust prior independent of the ansatz for the reduction, guaranteed to reproduce the assumptions of the theoretical model and exhaustive in the regime in which the feature is easily distinguishable from the baseline cosmology. We find a fit to the $\ell\approx20$--$40$ minus/plus structure in Planck TT power spectrum, as well as features spanning along the higher $\ell$'s ($\ell\approx100$--$1500$). For the last ones, we compute the correlated features that we expect to find in the CMB bispectrum, and asses their signal-to-noise and correlation to the ISW-lensing secondary bispectrum. We compare our findings to the shape-agnostic oscillatory template tested in Planck 2015, and we comment on some tantalising coincidences with some of the traits described in Planck's 2015 bispectrum data.

EFTCAMB/EFTCosmoMC are publicly available patches to the CAMB/CosmoMC codes implementing the effective field theory approach to single scalar field dark energy and modified gravity models. With the present numerical notes we provide a guide to the technical details of the code. Moreover we reproduce, as they appear in the code, the complete set of the modified equations and the expressions for all the other relevant quantities used to construct these patches. We submit these notes to the arXiv to grant full and permanent access to this material which provides very useful guidance to the numerical implementation of the EFT framework. We will update this set of notes when relevant modifications to the EFTCAMB/EFTCosmoMC codes will be released. The present version is based on the version of EFTCAMB/EFTCosmoMC Sep17.

Future observations of cosmic microwave background (CMB) polarisation have the potential to answer some of the most fundamental questions of modern physics and cosmology. In this paper, we list the requirements for a future CMB polarisation survey addressing these scientific objectives, and discuss the design drivers of the CORE space mission proposed to ESA in answer to the "M5" call for a medium-sized mission. The rationale and options, and the methodologies used to assess the mission's performance, are of interest to other future CMB mission design studies. CORE is designed as a near-ultimate CMB polarisation mission which, for optimal complementarity with ground-based observations, will perform the observations that are known to be essential to CMB polarisation scienceand cannot be obtained by any other means than a dedicated space mission.

The BeSSeL Survey is mapping the spiral structure of the Milky Way by measuring trigonometric parallaxes of hundreds of maser sources associated with high-mass star formation. While parallax techniques for water masers at high frequency (22 GHz) have been well documented, recent observations of methanol masers at lower frequency (6.7 GHz) have revealed astrometric issues associated with signal propagation through the ionosphere that could significantly limit parallax accuracy. These problems displayed as a "parallax gradient" on the sky when measured against different background quasars. We present an analysis method in which we generate position data relative to an "artificial quasar" at the target maser position at each epoch. Fitting parallax to these data can significantly mitigate the problems and improve parallax accuracy.

Environment, such as the accretion disk, could modify the signal of the gravitational wave from the astrophysical black hole binaries. In this article, we model the matter field around the intermediate-mass binary black holes by means of an axion-like scalar field and investigate their joint evolution. In details, we consider the equal mass binary black holes surrounded by a shell of axion-like scalar field both in spherical symmetric and non-spherical symmetric cases, and with different strength of the scalar field. Our result shows that the environmental scalar field could essentially modify the dynamics. Firstly, in the spherical symmetric case, with increasing of the scalar field strength, the number of circular orbit of the binary black hole is reduced. It means that the scalar field could significantly accelerate the merger process. Secondly, once the scalar field strength exceeds certain critical value, the scalar field could collapse into a third black hole with its mass being larger than the binary. Consequently, the new black hole collapsed from the environmental scalar field could accrete the binary promptly and the binary collides head-on between each other. In this process, there is almost no any quadrupole signal produced, namely the gravitational wave is greatly suppressed. Thirdly, when the scalar field strength is relatively smaller than the critical value, the black hole orbit could develop eccentricity through the accretion of the scalar field. Fourthly, during the initial stage of the inspire, the gravitational attractive force from the axion-like scalar field could induce a sudden turn in the binary orbits, hence result in a transient wiggle in the gravitational waveform. Finally, in the non-spherical case, the scalar field could gravitationally attract the binary moving toward the mass center of the scalar field and slow down the merger process.

By using the effective field theory approach, we investigate the role of initial condition for the dark energy or modified gravity models. In details, we consider the constant and linear parametrization of the effective Newton constant models. Firstly, under the adiabatic assumption, the correction from the extra scalar degree of freedom in the beyond $\Lambda$CDM model is found to be negligible. The dominant ingredient in this setup is the primordial curvature perturbation originated from inflation mechanism, and the energy budget of the matter components is not very crucial. Secondly, the iso-curvature perturbation sourced by the extra scalar field is studied. For the constant and linear model of the effective Newton constant, there is no such kind of scalar mode exist. For the quadratic model, there is a non-trivial one. However, the amplitude of the scalar field is damped away very fast on all scales. Consequently, it could not support a reasonable structure formation. Finally, we study the importance of the setup of the scalar field starting time. By setting different turn-on time, namely $a=10^{-2} $ and $a=10^{-7} $, we compare the cosmic microwave background radiation temperature, lensing deflection angle auto-correlation function as well as the matter power spectrum in the constant and linear model. We find there is an order of $\mathcal{O}(1\%)$ difference in the observable spectra for constant model, while for the linear model, it is smaller than $\mathcal{O}(0.1\%)$.

We forecast the scientific capabilities to improve our understanding of cosmic inflation of CORE, a proposed CMB space satellite submitted in response to the ESA fifth call for a medium-size mission opportunity. The CORE satellite will map the CMB anisotropies in temperature and polarization in 19 frequency channels spanning the range 60-600 GHz. CORE will have an aggregate noise sensitivity of $1.7 \mu$K$\cdot \,$arcmin and an angular resolution of 5' at 200 GHz. We explore the impact of telescope size and noise sensitivity on the inflation science return by making forecasts for several instrumental configurations. This study assumes that the lower and higher frequency channels suffice to remove foreground contaminations and complements other related studies of component separation and systematic effects, which will be reported in other papers of the series "Exploring Cosmic Origins with CORE." We forecast the capability to determine key inflationary parameters, to lower the detection limit for the tensor-to-scalar ratio down to the $10^{-3}$ level, to chart the landscape of single field slow-roll inflationary models, to constrain the epoch of reheating, thus connecting inflation to the standard radiation-matter dominated Big Bang era, to reconstruct the primordial power spectrum, to constrain the contribution from isocurvature perturbations to the $10^{-3}$ level, to improve constraints on the cosmic string tension to a level below the presumptive GUT scale, and to improve the current measurements of primordial non-Gaussianities down to the $f_{NL}^{\rm local} < 1$ level. For all the models explored, CORE alone will improve significantly on the present constraints on the physics of inflation. Its capabilities will be further enhanced by combining with complementary future cosmological observations.

Cosmological observables show a dependence with the neutrino mass, which is partially degenerate with parameters of extended models of gravity. We study and explore this degeneracy in Horndeski generalized scalar-tensor theories of gravity. Using forecasted cosmic microwave background and galaxy power spectrum datasets, we find that a single parameter in the linear regime of the effective theory dominates the correlation with the total neutrino mass. For any given mass, a particular value of this parameter approximately cancels the power suppression due to the neutrino mass at a given redshift. The extent of the cancellation of this degeneracy depends on the cosmological large-scale structure data used at different redshifts. We constrain the parameters and functions of the effective gravity theory and determine the influence of gravity on the determination of the neutrino mass from present and future surveys.

We conducted VLA C-configuration observations to measure positions and luminosities of Galactic Class II 6.7 GHz methanol masers and their associated ultra-compact H II regions. The spectral resolution was 3.90625 kHz and the continuum sensitivity reached 45 \uJypb. We mapped 372 methanol masers with peak flux densities of more than 2 Jy selected from the literature, 367 of them were detected. Absolute positions have nominal uncertainties of 0.3 arcsec. In this first paper on the data analysis, we present three catalogs, the first gives information on the strongest feature of 367 methanol maser sources, and the second on all detected maser spots. The third catalog present derived data of the 279 radio continuum sources found in the vicinity of maser sources. Among them, 140 show evidence of physical association with maser sources. Our catalogs list properties including distance, flux density, radial velocity and the distribution of masers on the Galactic plane is then provided as well. We found no significant relationship between luminosities of masers and their associated radio continuum counterparts.

HERD is the High Energy cosmic-Radiation Detection instrument proposed to operate onboard China's space station in the 2020s. It is designed to detect energetic cosmic ray nuclei, leptons and photons with a high energy resolution ($\sim1\%$ for electrons and photons and $20\%$ for nuclei) and a large geometry factor ($>3\,{ m^2\,sr}$ for electrons and diffuse photons and $>2\,{ m^2\,sr}$ for nuclei). In this work we discuss the capability of HERD to detect monochromatic $\gamma$-ray lines, based on simulations of the detector performance. It is shown that HERD will be one of the most sensitive instruments for monochromatic $\gamma$-ray searches at energies between $\sim10$ to a few hundred GeV. Above hundreds of GeV, Cherenkov telescopes will be more sensitive due to their large effective area. As a specific example, we show that a good portion of the parameter space of a supersymmetric dark matter model can be probed with HERD.

We consider Horava gravity within the framework of the effective field theory (EFT) of dark energy and modified gravity. We work out a complete mapping of the theory into the EFT language for an action including all the operators which are relevant for linear perturbations with up to sixth order spatial derivatives. We then employ an updated version of the EFTCAMB/EFTCosmoMC package to study the cosmology of the low-energy limit of Horava gravity and place constraints on its parameters using several cosmological data sets. In particular we use cosmic microwave background (CMB) temperature-temperature and lensing power spectra by Planck 2013, WMAP low-l polarization spectra, WiggleZ galaxy power spectrum, local Hubble measurements, Supernovae data from SNLS, SDSS and HST and the baryon acoustic oscillations measurements from BOSS, SDSS and 6dFGS. We get improved upper bounds, with respect to those from Big Bang Nucleosynthesis, on the deviation of the cosmological gravitational constant from the local Newtonian one. At the level of the background phenomenology, we find a relevant rescaling of the Hubble rate at all epoch, which has a strong impact on the cosmological observables; at the level of perturbations, we discuss in details all the relevant effects on the observables and find that in general the quasi-static approximation is not safe to describe the evolution of perturbations. Overall we find that the effects of the modifications induced by the low-energy Horava gravity action are quite dramatic and current data place tight bounds on the theory parameters.

Despite its continued observational successes, there is a persistent (and growing) interest in extending cosmology beyond the standard model, $\Lambda$CDM. This is motivated by a range of apparently serious theoretical issues, involving such questions as the cosmological constant problem, the particle nature of dark matter, the validity of general relativity on large scales, the existence of anomalies in the CMB and on small scales, and the predictivity and testability of the inflationary paradigm. In this paper, we summarize the current status of $\Lambda$CDM as a physical theory, and review investigations into possible alternatives along a number of different lines, with a particular focus on highlighting the most promising directions. While the fundamental problems are proving reluctant to yield, the study of alternative cosmologies has led to considerable progress, with much more to come if hopes about forthcoming high-precision observations and new theoretical ideas are fulfilled.

Lensing of the Cosmic Microwave Background (CMB) changes the morphology of pattern of temperature fluctuations, so topological descriptors such as Minkowski Functionals can probe the gravity model responsible for the lensing. We show how the recently introduced two-to-two and three-to-one kurt-spectra (and their associated correlation functions), which depend on the power spectrum of the lensing potential, can be used to probe modified gravity theories such as $f({R})$ theories of gravity and quintessence models. We also investigate models based on effective field theory, which include the constant-$\Omega$ model, and low-energy Ho\vrava theories. Estimates of the cumulative signal-to-noise for detection of lensing-induced morphology changes, reaches ${\cal O}(10^3)$ for the future planned CMB polarization mission COrE$^{+}$. Assuming foreground removal is possible to $\ell_{max}=3000$, we show that many modified gravity theories can be rejected with a high level of significance, making this technique comparable in power to galaxy weak lensing or redshift surveys. These topological estimators are also useful in distinguishing {\em lensing} from other scattering secondaries at the level of the four-point function or trispectrum. Examples include the kinetic Sunyaev-Zel'dovich (kSZ) effect which shares, with lensing, a lack of spectral distortion. We also discuss the complication of foreground contamination from unsubtracted point sources.

We show how to fully map a specific model of modified gravity into the Einstein-Boltzmann solver EFTCAMB. This approach consists in few steps and allows to obtain the cosmological phenomenology of a model with minimal effort. We discuss all these steps, from the solution of the dynamical equations for the cosmological background of the model to the use of the mapping relations to cast the model into the effective field theory language and use the latter to solve for perturbations. We choose the Hu-Sawicki f(R) model of gravity as our working example. After solving the background and performing the mapping, we interface the algorithm with EFTCAMB and take advantage of the effective field theory framework to integrate the full dynamics of linear perturbations, returning all quantities needed to accurately compare the model with observations. We discuss some observational signatures of this model, focusing on the linear growth of cosmic structures. In particular we present the behavior of $f\sigma_8$ and $E_G$ that, unlike the $\Lambda$CDM scenario, are generally scale dependent in addition to redshift dependent. Finally, we study the observational implications of the model by comparing its cosmological predictions to the Planck 2015 data, including CMB lensing, the WiggleZ galaxy survey and the CFHTLenS weak lensing survey measurements. We find that while WiggleZ data favor a non-vanishing value of the Hu-Sawicki model parameter, $\log_{10}(-f^0_{R})$, and consequently a large value of $\sigma_8$, CFHTLenS drags the estimate of $\log_{10}(-f^0_{R})$ back to the $\Lambda$CDM limit.

The latest Planck results show a power deficit in the temperature anisotropies near $\ell \approx 20$ in the cosmic microwave background (CMB). This observation can hardly be explained within the standard inflationary $\Lambda$-cold-dark-matter ($\Lambda$CDM) scenario. In this paper we consider a string theory inspired inflationary model (axion monodromy inflation) with a step-like modulation in the potential which gives rise to observable signatures in the primordial perturbations. One interesting phenomenon is that the primordial scalar modes experience a sudden suppression at a critical scale when the modulation occurs. By fitting to the CMB data, we find that the model can nicely explain the $\ell \approx 20$ power deficit anomaly as well as predict specific patterns in the temperature-polarization correlation and polarization autocorrelation spectra. Though the significance of the result is not sufficient to claim a detection, our analysis reveals that fundamental physics at extremely high energy scales, namely, some effects inspired by string theory, may be observationally testable in forthcoming cosmological experiments.

We continue the study of mild transient reductions in the speed of sound of the adiabatic mode during inflation, of their effect on the primordial power spectrum and bispectrum, and of their detectability in the Cosmic Microwave Background (CMB). We focus on the regime of \emph{moderately sharp} mild reductions in the speed of sound during uninterrupted slow-roll inflation, a theoretically well motivated and self-consistent regime that admits an effective single-field description. The signatures on the power spectrum and bispectrum were previously computed using a slow-roll Fourier transform (SRFT) approximation, and here we compare it with generalized slow-roll (GSR) and in-in methods, for which we derive new formulas that account for moderately sharp features. The agreement between them is excellent, and also with the power spectrum obtained from the numerical solution to the equation of motion. We show that, in this regime, the SRFT approximation correctly captures with simplicity the effect of higher derivatives of the speed of sound in the mode equation, and makes manifest the correlations between power spectrum and bispectrum features. In a previous paper we reported hints of these correlations in the Planck data and here we perform several consistency checks and further analyses of the best fits, such as polarization and local significance at different angular scales. For the data analysis, we show the excellent agreement between the CLASS and CAMB Boltzmann codes. Our results confirm that the theoretical framework is consistent, and they suggest that the predicted correlations are robust enough to be searched for in CMB and Large Scale Structure (LSS) surveys.

Inspired by the study of mild transient reductions in the speed of sound of the adiabatic mode during inflation, we search for a primordial localized feature imprinted in cosmic microwave background and large-scale structure formation observables. We find some common oscillatory patterns both in the Planck CMB temperature-temperature power spectrum and the WiggleZ galaxy spectrum. By performing independent searches with these two data sets, we find a coincidence in the most significant mode previously found by Ach\'ucarro et al. 2013 by using only Planck data. Furthermore, the joint data analysis shows that the oscillation frequency of the feature gets better constrained, and the amplitude marginally deviates from zero, unlike what was observed using only Planck data. Besides the parameter estimation, we also discuss the Bayesian evidence. The addition of WiggleZ data mildly enhances the significance of the best mode found in the Planck data.

We present a description for setting initial particle displacements and field values for simulations of arbitrary metric theories of gravity, for perfect and imperfect fluids with arbitrary characteristics. We extend the Zel'dovich Approximation to nontrivial theories of gravity, and show how scale dependence implies curved particle paths, even in the entirely linear regime of perturbations. For a viable choice of Effective Field Theory of Modified Gravity, initial conditions set at high redshifts are affected at the level of up to 5% at Mpc scales, which exemplifies the importance of going beyond {\Lambda}-Cold Dark Matter initial conditions for modifications of gravity outside of the quasi-static approximation. In addition, we show initial conditions for a simulation where a scalar modification of gravity is modelled in a Lagrangian particle-like description. Our description paves the way for simulations and mock galaxy catalogs under theories of gravity beyond the standard model, crucial for progress towards precision tests of gravity and cosmology.

Planck-2015 data seem to favour a large value of the lensing amplitude parameter, $A_{\rm L}=1.22\pm0.10$, in CMB spectra. This result is in $2\sigma$ tension with the lensing reconstruction result, $A_{\rm L}^{\phi\phi}=0.95\pm0.04$. In this paper, we simulate several CMB anisotropy and CMB lensing spectra based on Planck-2015 best-fit cosmological parameter values and Planck blue book beam and noise specifications. We analyse several modified gravity models within the effective field theory framework against these simulations and find that models whose effective Newton constant is enhanced can modulate the CMB anisotropy spectra in a way similar to that of the $A_{\rm L}$ parameter. However, in order to lens the CMB anisotropies sufficiently, like in the Planck-2015 results, the growth of matter perturbations is substantially enhanced and gives a high $\sigma_8$ value. This in turn proves to be problematic when combining these data to other probes, like weak lensing from CFHTLenS, that favour a smaller amplitude of matter fluctuations.

We revisit the degeneracy between massive neutrinos and generalized theories of gravity in the framework of effective field theory of cosmic acceleration. In particular we consider f(R) theories and a class of non-minimally coupled models parametrized via a coupling to gravity which is linear in the scale factor. In the former case, we find a slightly lower degeneracy with respect to what found so far in the literature, due to the fact that we implement exact designer f(R) models and evolve the full linear dynamics of perturbations. As a consequence, our bounds are slightly tighter on the f(R) parameter but looser on the summed neutrino mass. We also set a new upper bound on the Compton wavelength parameter ${\rm Log}_{10}B_0<-4.1$ at 95% C.L. with fixed summed neutrino mass ($\sum m_{\nu}=0.06$ eV) in f(R) gravity with the combined data sets from cosmic microwave background temperature and lensing power spectra of Planck collaboration as well as galaxy power spectrum from WiggleZ dark energy survey. We do not observe a sizable degeneracy between massive neutrinos and modified gravity in the linear parametrization of non-minimally gravitational coupling model. The analysis is performed with an updated version of the EFTCAMB/EFTCosmoMC package, which is now publicly available and extends the first version of the code with the consistent inclusion of massive neutrinos, tensor modes, several alternative background histories and designer quintessence models.

We introduce EFTCAMB/EFTCosmoMC as publicly available patches to the commonly used CAMB/CosmoMC codes. We briefly describe the structure of the codes, their applicability and main features. To illustrate the use of these patches, we obtain constraints on parametrized pure EFT and designer $f(R)$ models, both on $\Lambda$CDM and $w$CDM background expansion histories, using data from Planck temperature and lensing potential spectra, WMAP low-$\ell$ polarization spectra (WP), and baryon acoustic oscillation (BAO). Upon inspecting theoretical stability of the models on the given background, we find non-trivial parameter spaces that we translate into viability priors. We use different combinations of data sets to show their individual effects on cosmological and model parameters. Our data analysis results show that, depending on the adopted data sets, in the $w$CDM background case this viability priors could dominate the marginalized posterior distributions. Interestingly, with Planck+WP+BAO+lensing data, in $f(R)$ gravity models, we get very strong constraints on the constant dark energy equation of state, $w_0\in(-1,-0.9997)\,\,(95% {\rm C.L.})$.

We use the optimised skew-spectrum as well as the skew-spectra associated with the Minkowski Functionals (MFs) to test the possibility of using the cross-correlation of the Integrated Sachs-Wolfe effect (ISW) and lensing of the cosmic microwave background (CMB) radiation to detect deviations in the theory of gravity away from General Relativity (GR). We find that the although both statistics can put constraints on modified gravity, the optimised skew-spectra are especially sensitive to the parameter $\rm B_0$ that denotes the the {\em Compton wavelength} of the scalaron at the present epoch. We investigate three modified gravity theories, namely: the Post-Parametrised Friedmanian (PPF) formalism; the Hu-Sawicki (HS) model; and the Bertschinger-Zukin (BZ) formalism. Employing a likelihood analysis for an experimental setup similar to ESA's Planck mission, we find that, assuming GR to be the correct model, we expect the constraints from the first two skew-spectra, $S_{\ell}^{(0)}$ and $S_{\ell}^{(1)}$, to be the same: $\rm B_0<0.45$ at $95%$ confidence level (CL), and $\rm B_0<0.67$ at $99%$ CL in the BZ model. The third skew-spectrum does not give any meaningful constraint. We find that the optimal skew-spectrum provides much more powerful constraint, giving $\rm B_0<0.071$ at $95%$ CL and $\rm B_0<0.15$ at $99%$ CL, which is essentially identical to what can be achieved using the full bispectrum.

We implement the effective field theory (EFT) approach to dark energy and modified gravity in the public Einstein-Boltzmann solver CAMB. The resulting code, which we dub EFTCAMB, is a powerful and versatile tool that can be used for several objectives. It can be employed to evolve the full dynamics of linear scalar perturbations in any given single field dark energy or modified gravity model, once the latter is mapped into the EFT formalism. It offers a numerical implementation of EFT as a model-independent framework to test gravity on cosmological scales. EFTCAMB has a built-in check for the fulfillment of general stability conditions such as the absence of ghost and superluminal propagation of perturbations. It handles phantom-divide crossing models and does not contain any quasi-static approximation, but rather evolves the full dynamics of perturbations on linear scales. As we will show, the latter is an important feature in view of the accuracy and scale range of upcoming surveys. We show the reliability and applicability of our code by evolving the dynamics of linear perturbations and extracting predictions for power spectra in several models. In particular we perform a thorough analysis of f(R) theories, comparing our outputs with those of an existing code for LCDM backgrounds, and finding an agreement that can reach 0.1% for models with a Compton wavelength consistent with current cosmological data. We then showcase the flexibility of our code studying two different scenarios. First we produce new results for designer f(R) models with a time-varying dark energy equation of state. Second, we extract predictions for linear observables in some parametrized EFT models with a phantom-divide crossing equation of state for dark energy.

By using the cubic spline interpolation method, we reconstruct the shape of the primordial scalar and tensor power spectra from the recently released {\it Planck} temperature and BICEP2 polarization cosmic microwave background data. We find that the vanishing scalar index running ($\dd n_s/\dd\ln k$) model is strongly disfavored at more than $3\sigma$ confidence level on the $k=0.0002$ Mpc$^{-1}$ scale. Furthermore, the power-law parameterization gives a blue-tilt tensor spectrum, no matter using only the first 5 bandpowers $n_t = 1.20^{+0.56}_{-0.64} (95% {\rm CL})$ or the full 9 bandpowers $n_t = 1.24^{+0.51}_{-0.58} (95% {\rm CL})$ of BICEP2 data sets. Unlike the large tensor-to-scalar ratio value ($r\sim0.20$) under the scale-invariant tensor spectrum assumption, our interpolation approach gives $r_{0.002} < 0.060 (95% {\rm CL})$ by using the first 5 bandpowers of BICEP2 data. After comparing the results with/without BICEP2 data, we find that {\it Planck} temperature with small tensor amplitude signals and BICEP2 polarization data with large tensor amplitude signals dominate the tensor spectrum reconstruction on the large and small scales, respectively. Hence, the resulting blue tensor tilt actually reflects the tension between {\it Planck} and BICEP2 data.

In order to explain cosmic acceleration without invoking "dark" physics, we consider $f(R)$ modified gravity models, which replace the standard Einstein-Hilbert action in General Relativity with a higher derivative theory. We use data from the WiggleZ Dark Energy survey to probe the formation of structure on large scales which can place tight constraints on these models. We combine the large-scale structure data with measurements of the cosmic microwave background from the Planck surveyor. After parameterising the modification of the action using the Compton wavelength parameter $B_0$, we constrain this parameter using ISiTGR, assuming an initial non-informative log prior probability distribution of this cross-over scale. We find that the addition of the WiggleZ power spectrum provides the tightest constraints to date on $B_0$ by an order of magnitude, giving ${\rm log}_{10}(B_0) < -4.07$ at 95% confidence limit. Finally, we test whether the effect of adding the lensing amplitude $A_{\rm Lens}$ and the sum of the neutrino mass $\sum m_\nu$ is able to reconcile current tensions present in these parameters, but find $f(R)$ gravity an inadequate explanation.

We update the cosmological parameter estimation for three non-vanilla models by a joint analysis of \CCCP\ X-ray cluster, the newly released \Planck\ CMB data as well as some external data sets, such as baryon acoustic oscillation measurements from the 6dFGS, SDSS DR7 and BOSS DR9 surveys, and Hubble Space Telescope $H_0$ measurement. First of all, we find that X-ray cluster data sets strongly favor a non-zero summed neutrino mass at more than 3$\sigma$ confidence level in these non-vanilla models. And then, we reveal some tensions between X-ray cluster and {\it Planck} data in some cosmological parameters. For the matter power spectrum amplitude $\sigma_8$, X-ray cluster data favor a lower value compared with {\it Planck}. Because of the strong $\sigma_8-\sum m_{\nu}$ degeneracy, this tension could beyond 2$\sigma$ confidence level when the summed neutrino mass $\sum m_{\nu}$ is allowed to vary. For the CMB lensing amplitude $A_L$, the addition of X-ray cluster data results in a 3$\sigma$ deviation from the vanilla model. Furthermore, {\it Planck}+X-ray data prefer a large Hubble constant and phantom-like dark energy equation of state, which are in $2\sigma$ tension with those from WMAP7+X-ray data. Finally, we find that these tensions/descrepencies could be relaxed in some sense by adding a $9\%$ systematic shift in the cluster mass functions.

We constrain $f(R)$ and chameleon-type modified gravity in the framework of the Berstchinger-Zukin parametrization using the recent released Planck data, including both CMB temperature power spectrum and lensing potential power spectrum. Some other external data sets are included, such as BAO measurements from the 6dFGS, SDSS DR7 and BOSS DR9 surveys, HST $H_0$ measurement and supernovae from Union2.1 compilation. We also use WMAP9yr data for consistency check and comparison. For $f(R)$ gravity, WMAP9yr results can only give quite a loose constraint on the modified gravity parameter $B_0$, which is related to the present value of the Compton wavelength of the extra scalar degree of freedom, $B_0<3.37$ at $95\% {\rm C.L.}$ We demonstrate that this constraint mainly comes from the late ISW effect. With only Planck CMB temperature power-spectrum data, we can improve the WMAP9yr result by a factor $3.7$ ($B_0<0.91$ at $95\% {\rm C.L.}$). If the Planck lensing potential power-spectrum data are also taken into account, the constraint can be further strenghtened by a factor $5.1$ ($B_0<0.18$ at $95\% {\rm C.L.}$). This major improvement mainly comes from the small-scale lensing signal. Furthermore, BAO, HST and supernovae data could slightly improve the $B_0$ bound ($B_0<0.12$ at $95\% {\rm C.L.}$).For the chameleon-type model, we find that the data set which we used cannot constrain the Compton wavelength $B_0$ and the potential index $s$ of chameleon field, but can give a tight constraint on the parameter $\beta_1=1.043^{+0.163}_{-0.104}$ at $95\% {\rm C.L.}$ ($\beta_1=1$ in general relativity), which accounts for the non-minimal coupling between the chameleon field and the matter component. In addition, we find that both modified gravity models we considered favor a relatively higher Hubble parameter than the concordance LCDM model in general relativity.

We forecast the constraints on both Hu-Sawicki model and Bertschinger-Zukin model of modified gravity within the Parameterized Post-Friedmann (PPF) formalism for the Planck satellite experiment by performing the joint analysis of ISW-Lensing bispectrum and CMB power spectrum. We find that, even considering the temperature-temperature mode of CMB power spectrum only, Planck data are expected to reduce the error bars on the modified gravity parameter $B_0$ (related to the present value of Compton wavelength of the extra scalar degree of freedom) at least one order magnitude compared with WMAP. The spectrum-bispectrum joint analysis can further improve the results by a factor ranging from 1.14 to 5.32 depending on the specific modified gravity model. One of our main results is that the cross-correlation between ISW-Lensing bispectrum and power spectrum can be safely neglected when performing the joint analysis. For simplicity, we only investigate the likelihood of one parameter ({$B_0$}) and fix all other cosmological parameters to their best-fit values in WMAP7yr results.

In this paper we study the non-Gaussian features of the primordial fluctuations in loop quantum cosmology with the inverse volume corrections. The detailed analysis is performed in the single field slow-roll inflationary models. However, our results reflect the universal characteristics of bispectrum in loop quantum cosmology. The main corrections to the scalar bispectrum come from two aspects: one is the modifications to the standard Bunch-Davies vacuum, the other is the corrections to the background dependent variables, such as slow-roll parameters. Our calculations show that the loop quantum corrections make $f_{{\rm NL}}$ of the inflationary models increase 0.1%. Moreover, we find that two new shapes arise, namely $\mathcal F_{1}$ and $\mathcal F_{2}$. The former gives a unique loop quantum feature which is less correlated with the local, equilateral and single types, while the latter is highly correlated with the local one.

The definition of the Galactic coordinate system was announced by the IAU Sub-Commission 33b on behalf of the IAU in 1958. For more than 50 years the definition of the Galactic coordinate system has remained unchanged from this IAU1958 version. On the basis of deep and all-sky catalogs, the position of the Galactic plane can be revised and updated definitions of the Galactic coordinate systems can be proposed. We re-determine the position of the Galactic plane based on modern large catalogs, such as the Two Micron All-Sky Survey (2MASS) and the SPECFIND v2.0. This paper also aims to propose a possible definition of the optimal Galactic coordinate system by adopting the ICRS position of the Sgr A* at the Galactic center. The near-infrared 2MASS point source catalog and the SPECFIND v2.0 catalog of radio continuum spectra are used to determine the mean position of the Galactic plane on the celestial sphere. By fitting the data to an ideal Galactic equator, the parameters defining the Galactic coordinate system are obtained. We find that the obliquity of the Galactic equator on the ICRS principal plane is about $0.4^\circ$ (2MASS) and $0.6^\circ$ (SPECFIND v2.0) larger than the J2000.0 value, which is widely used in coordinate transformations between the equatorial $(\alpha, \delta)$ and the Galactic $(\ell, b)$. Depending on the adopted parameters, data, and methods, the largest difference between the resulting Galactic coordinate systems is several arcminutes. We derive revised transformation matrices and parameters describing the orientation of the Galactic coordinate systems in the ICRS at the 1 milli-arcsecond level to match the precision of modern observations. For practical applications, we propose that a revised definition of the Galactic coordinate system should be required in the near future.

In this note we study the linear dynamics of scalar graviton in a de Sitter background in the infrared limit of the healthy extension of Ho\v{r}ava-Lifshitz gravity with the dynamical critical exponent $z=3$. Both our analytical and numerical results show that the non-zero Fourier modes of scalar graviton oscillate with an exponentially damping amplitude on the sub-horizon scale, while on the super-horizon scale, the phases are frozen and they approach to some asymptotic values. In addition, as the case of the non-zero modes on super-horizon scale, the zero mode also initially decays exponentially and then approaches to an asymptotic constant value.

Using the full radiation transfer function, we numerically calculate the CMB angular bispectrum seeded by the compensated magnetic scalar density mode. We find that, for the string inspired primordial magnetic fields characterized by index $n_B=-2.9$ and mean-field amplitude $B_{\lam}=9{\rm nG}$, the angular bispectrum is dominated by two primordial magnetic shapes. The first magnetic shape looks similar to the one from local-type primordial curvature perturbations, so both the amplitude and profile of the Komatsu-Spergel estimator (reduced bispectrum) seeded by this shape are almost the same as those of the primary CMB anisotropies. However, for different parameter sets ($l_1,l_2$), this "local-type" reduced bispectrum oscillates around different asymptotic values in the high-$l_3$ regime because of the effect of the Lorentz force, which is exerted by the primordial magnetic fields on the charged baryons. This feature is different from the standard case where all modes approach to zero asymptotically in the high-$l$ limit. On the other hand, the second magnetic shape appears only in the primordial magnetic field model. The amplitude of the Komatsu-Spergel estimator sourced by the second shape diverges in the low-$l$ regime because of the negative slope of shape. In the high-$l$ regime, this amplitude is approximately equal to that of the first estimator, but with a reversal phase.

In this paper we study the entropy perturbations in N-flation by using the $\d\ma{N}$ formalism. We calculate the entropy corrections to the power spectrum of the overall curvature perturbation $P_{\z}$. We obtain an analytic form of the transfer coefficient $T^2_{\ma{R}\ma{S}}$, which describes the correlation between the curvature and entropy perturbations, and investigate its behavior numerically. It turns out that the entropy perturbations cannot be neglected in N-flation, because the amplitude of entropy components is approximately in the same order as the adiabatic one at the end of inflation $T^2_{\ma{R}\ma{S}}\sim\ma{O}(1)$. The spectral index $n_S$ is calculated and it becomes smaller after the entropy modes are taken into account, i.e., the spectrum becomes redder, compared to the pure adiabatic case. Finally we study the modified consistency relation of N-flation, and find that the tensor-to-scalar ratio ($r\simeq0.006$) is greatly suppressed by the entropy modes, compared to the pure adiabatic one ($r\simeq0.017$) at the end of inflation.

We study the dynamics of a timelike vector field which violates Lorentz invariance when the background spacetime is in an accelerating phase in the early universe. It is shown that a timelike vector field is difficult to realize an inflationary phase, so we investigate the evolution of a vector field within a scalar field driven inflation model. And we calculate the power spectrum of the vector field without considering the metric perturbations. While the time component of the vector field perturbations provides a scale invariant spectrum when $\xi = 0$, where $\xi$ is a nonminimal coupling parameter, both the longitudinal and transverse perturbations give a scale invariant spectrum when $\xi = 1/6$.

We compute the large scalar four-point correlation functions in general single field inflation models, where the inflaton Lagrangian is an arbitrary function of the inflaton and its first derivative. We find that the leading order trispectra have four different shapes determined by three parameters. We study features in these shapes that can be used to distinguish among themselves, and between them and the trispectra of the local form. For the purpose of data analyses, we give two simple representative forms for these "equilateral trispectra". We also study the effects on the trispectra if the initial state of inflation deviates from the standard Bunch-Davies vacuum.

We investigate the primordial trispectra of the general multifield DBI inflationary model. In contrast with the single field model, the entropic modes can source the curvature perturbations on the super horizon scales, so we calculate the contributions from the interaction of four entropic modes mediating one adiabatic mode to the trispectra, at the large transfer limit ($T_{RS}\gg1$). We obtained the general form of the 4-point correlation functions, plotted the shape diagrams in two specific momenta configurations, "equilateral configuration" and "specialized configuration". Our figures showed that we can easily distinguish the two different momenta configurations.

We investigate the linear cosmological perturbations of Ho\v{r}ava-Lifshitz gravity in a FRW universe without any matter. Our results show that a new gauge invariant dynamical scalar mode emerges, due to the gauge transformation under the "foliation-preserving" diffeomorphism and "projectability condition", and it can produce a scale invariant power spectrum. In the infrared regime with $\lambda=1$, the dynamical scalar degree of freedom turns to be a non-dynamical one at the linear order level.

The constraint on the total energy in a given spatial region is given from holography by the mass of a black hole which just fits in that region, which leads to an UV/IR relation: the maximal energy density in that region is proportional to $M_p^2/L^2$, where $M_p$ is the Planck mass and $L$ is the spatial scale of that region under consideration. Assuming the maximal black hole in the universe is formed through gravitational collapse of perturbations in the universe, then the "Jeans" scale of the perturbations gives a causal connection scale $R_{\rm CC}$. For gravitational perturbations, $R^{-2}_{\rm CC}={\rm Max}(\dot H+2H^2, -\dot H)$ for a flat universe. We study the cosmological dynamics of the corresponding vacuum energy density by choosing the causal connection scale as the IR cutoff in the UV/IR relation, in the cases of the vacuum energy density as an independently conserved energy component and an effective dynamical cosmological constant, respectively. It turns out that only the case with the choice $R_{\rm CC}^{-2}= \dot H +2H^2$, could be consistent with the current cosmological observations when the vacuum density appears as an independently conserved energy component. In this case, the model is called holographic Ricci scalar dark energy model in the literature.

In this series we construct an effective field theory (EFT) in curved spacetime to study gravitational radiation and backreaction effects. We begin in this paper with a derivation of the self-force on a compact object moving in the background spacetime of a supermassive black hole. The EFT approach utilizes the disparity between two length scales, which in this problem are the size of the compact object and the radius of curvature of the background spacetime, to treat the orbital dynamics of the compact object, described as an effective point particle, separately from its tidal deformations. Ultraviolet divergences are regularized using Hadamard's {\it partie finie} to isolate the non-local finite part from the quasi-local divergent part. The latter is constructed from a momentum space representation for the graviton retarded propagator and is evaluated using dimensional regularization in which only logarithmic divergences are relevant for renormalizing the parameters of the theory. As a first important application of this framework we explicitly derive the first order self-force given by Mino, Sasaki, Tanaka, Quinn and Wald. Going beyond the point particle approximation, to account for the finite size of the object, we demonstrate that for extreme mass ratio inspirals the motion of a compact object is affected by tidally induced moments at $O(\epsilon^4)$, in the form of an Effacement Principle. The relatively large radius-to-mass ratio of a white dwarf star allows for these effects to be enhanced until the white dwarf becomes tidally disrupted, a potentially $O(\epsilon^2)$ process, or plunges into the supermassive black hole. This work provides a new foundation for further exploration of higher order self force corrections, gravitational radiation and spinning compact objects.

We analyze the phenomena of Bose Novae, as described by Donley et al [Nature 412, 295 (2001)], by focusing on the behavior of excitations or fluctuations above the condensate, as driven by the dynamics of the condensate (rather than the dynamics of the condensate alone or the kinetics of the atoms). The dynamics of the condensate squeezes and amplifies the quantum excitations, mixing the positive and negative frequency components of their wave functions thereby creating particles which appear as bursts and jets. By analyzing the changing amplitude and particle content of these excitations, our simple physical picture (based on a test field approximation) explains well the overall features of the Bose Novae phenomena and provide excellent quantitative fits with experimental data on several aspects, such as the scaling behavior of the collapse time and the amount of particles in the jet. The predictions of the bursts at this level of approximation is less than satisfactory but may be improved on by including the backreaction of the excitations on the condensate. The mechanism behind the dominant effect -- parametric amplification of vacuum fluctuations and freezing of modes outside of horizon -- is similar to that of cosmological particle creation and structure formation in a rapid quench (which is fundamentally different from Hawking radiation in black holes). This shows that BEC dynamics is a promising venue for doing `laboratory cosmology'.

Phenomena associated with the controlled collapse of a Bose-Einstein condensate described in the experiment of Donley et al (E. Donley et. al., Nature 412, 295 (2001); N. Claussen, Ph. D. Thesis, U. of Colorado (2003)) are explained here as a consequence of the squeezing and amplification of quantum fluctuations above the condensate by the condensate dynamics. In analyzing the changing amplitude and particle contents of these excitations, our simple physical picture provides excellent quantitative fits with experimental data on the scaling behavior of the collapse time and the amount of particles emitted in the jets.

We examine the recently observed deviation of the muon g - 2 from the Standard Model prediction within the framework of gravity mediated SUGRA models with R parity invariance. Universal soft breaking (mSUGRA) models, and models with non-universal Higgs and third generation squark/slepton masses at M_G are considered. All relic density constraints from stau-neutralino co-annihilation and large \tan\beta NLO corrections for b \to s\gamma decay are included, and we consider two possibilities for the light Higgs: m_h > 114 GeV and m_h > 120 GeV. The combined m_h, b \to s\gamma and a_{\mu} bounds give rise to lower bounds on \tan\beta and m_{1/2}, while the lower bound on a_{\mu} gives rise to an upper bounds on m_{1/2}. These bounds are sensitive to A_0, e.g. for m_h > 114 GeV, the 95% C.L. is \tan\beta > 7(5) for A_0 = 0(-4m_{1/2}), and for m_h > 120 GeV, \tan\beta > 15(10). The positive sign of the a_{\mu} deviation implies \mu > 0, eliminating the extreme cancellations in the dark matter neutralino-proton detection cross section so that almost all the SUSY parameter space should be accessible to future planned detectors. Most of the allowed parts of parameter space occur in the co-annihilation region where m_0 is strongly correlated with m_{1/2}. The lower bound on a_{\mu} then greatly reduces the allowed parameter space. Thus using 90% C. L. bounds on a_{\mu} we find for A_0 = 0 that \tan\beta \geq 10 and for \tan\beta \leq 40 that m_{1/2} = (290 - 550) GeV and m_0 = (70 - 300) GeV. Then the tri-lepton signal and other SUSY signals would be beyond the Tevatron Run II (except for the light Higgs), only the \tilde{\tau}_1 and h and (and for part of the parameter space) the \tilde{e}_1 will be accessible to a 500 GeV NLC, while the LHC would be able to see the full SUSY mass spectrum.

Using the concept of open systems where the classical geometry is treated as the system and the quantum matter field as the environment, we derive a fluctuation-dissipation theorem for semiclassical cosmology. This theorem which exists under very general conditions for dissipations in the dynamics of the system, and the noise and fluctuations in the environment, can be traced to the formal mathematical relation between the dissipation and noise kernels of the influence functional depicting the open system, and is ultimately a consequence of the unitarity of the closed system. In particular, for semiclassical gravity, it embodies the backreaction effect of matter fields on the dynamics of spacetime. The backreaction equation derivable from the influence action is in the form of a Einstein-Langevin equation. It contains a dissipative term in the equation of motion for the dynamics of spacetime and a noise term related to the fluctuations of particle creation in the matter field. Using the well-studied model of a quantum scalar field in a Bianchi Type-I universe we illustrate how this Langevin equation and the noise term are derived and show how the creation of particles and the dissipation of anisotropy during the expansion of the universe can be understood as a manifestation of this fluctuation-dissipation relation.

The quantum Brownian motion paradigm provides a unified framework where one can see the interconnection of some basic quantum statistical processes like decoherence, dissipation, particle creation, noise and fluctuation. We treat the case where the Brownian particle is coupled linearly to a bath of time dependent quadratic oscillators. While the bath mimics a scalar field, the motion of the Brownian particle modeled by a single oscillator could be used to depict the behavior of a particle detector, a quantum field mode or the scale factor of the universe. An important result of this paper is the derivation of the influence functional encompassing the noise and dissipation kernels in terms of the Bogolubov coefficients. This method enables one to trace the source of statistical processes like decoherence and dissipation to vacuum fluctuations and particle creation, and in turn impart a statistical mechanical interpretation of quantum field processes. With this result we discuss the statistical mechanical origin of quantum noise and thermal radiance from black holes and from uniformly- accelerated observers in Minkowski space as well as from the de Sitter universe discovered by Hawking, Unruh and Gibbons-Hawking. We also derive the exact evolution operator and master equation for the reduced density matrix of the system interacting with a parametric oscillator bath in an initial squeezed thermal state. These results are useful for decoherence and backreaction studies for systems and processes of interest in semiclassical cosmology and gravity. Our model and results are also expected to be useful for related problems in quantum optics. %\pacs {05.40.+j,03.65.Sq,98.80.Cq,97.60.Lf}

We begin by enumerating the many processes in gravitation and cosmology where quantum noise and fluctuations play an active role such as particle creation, galaxy formation and entropy generation. Using the influence functional we first explain the origin and nature of noise in quantum systems interacting with an environment at a finite temperature. With linear coupling to nonohmic baths or at low temperatures, colored noise and nonlocal dissipation would appear and for nonlinear coupling multiplicative noise is generally expected. We derive a generalized fluctuation- dissipation relation for these systems. Then using a model of quantum Brownian motion in a bath of parametric oscillators, we show how noise and dissipation can be related to the Bogolubov coefficients of parametric amplification, which in the second-quantized sense, depicts cosmological particle creation in a dynamic background. We then calculate the influence functional and study the noise characteristics of quantum fields as probed by a particle detector. As examples, we show that an uniformly- accelerated observer in flat space or an inertial observer in an exponentially expanding (de Sitter) universe would see a thermal particle spectrum, recovering the well-known results of Unruh and Gibbons and Hawking , as inspired by the Hawking effect in black holes. We show how this method can be effectively used for treating the backreaction of particle creation and other quantum field processes on the