We propose that the massless dark photons produced from the annihilation of keV dark fermions in the Galaxy can induce the excess of electron recoil events recently observed in the XENON1T experiment. The minimal model for this is the extension of a $U(1)_{X}$ gauge symmetry, under which the massless dark photon couples to both dark and visible matter currents. We find that the best-fit parameters of the dark sector are compatible with the most stringent constraints from stellar cooling. The dark fermion relic density in both freeze-in and freeze-out scenarios has been calculated and found to be able to account for the observed dark matter abundance.

It is known that the Higgs potential in the Standard Model can drive a successful inflation as long as the Higgs field couples non-minimally to gravity. It is then inevitable to take into account the loop corrections of the Standard Model particles to the Higgs potential in the Higgs inflation. In this paper, we discuss the one-loop corrections at finite temperature to the curvature perturbation generated during the Higgs inflation. We find that the thermal loop effects can suppress the power of the curvature perturbation at large scales, thus resulting in a low quadrupole of the temperature anisotropy in the cosmic microwave background.

We perform a-priori validation tests of subgrid-scale (SGS) models for the turbulent transport of momentum, energy and passive scalars. To this end, we conduct two sets of high-resolution hydrodynamical simulations with a Lagrangian code: an isothermal turbulent box with rms Mach number of 0.3, 2 and 8, and the classical wind tunnel where a cold cloud traveling through a hot medium gradually dissolves due to fluid instabilities. Two SGS models are examined: the eddy diffusivity (ED) model wildly adopted in astrophysical simulations and the "gradient model" due to Clark et al. (1979). We find that both models predict the magnitude of the SGS terms equally well (correlation coefficient > 0.8). However, the gradient model provides excellent predictions on the orientation and shape of the SGS terms while the ED model predicts poorly on both, indicating that isotropic diffusion is a poor approximation of the instantaneous turbulent transport. The best-fit coefficient of the gradient model is in the range of [0.16, 0.21] for the momentum transport, and the turbulent Schmidt number and Prandtl number are both close to unity, in the range of [0.92, 1.15].

We first present new Subaru narrow-band observations of the Ly$\alpha$ halo around the quasi-stellar object (QSO) CFHQ J232908$-$030158 at $z=6.42$, which appears the most luminous and extended halo at $z>5$ ($L_{Ly\alpha}=9.8\times10^{43}$ erg s$^{-1}$ within $37$ pkpc diameter). Then, combining these measurements with available data in the literature, we find two different evolutions of QSOs' Ly$\alpha$ halos. First is a possible short-term evolution with QSO age seen in four $z>6$ QSOs. We find the anti-correlation between the Ly$\alpha$ halo scales with QSOs' IR luminosity, with J2329-0301's halo being the brightest and largest. It indicates that ionizing photons escape more easily out to circum-galactic regions when host galaxies are less dusty. We also find a positive correlation between IR luminosity and black hole mass ($M_\text{BH}$). Given $M_\text{BH}$ as an indicator of QSO age, we propose a hypothesis that a large Ly$\alpha$ halo mainly exists around QSOs in the young phase of their activity due to a small amount of dust. The second is an evolution with cosmic time seen over $z\sim2-5$. We find the increase of surface brightness toward lower-redshift with a similar growth rate to that of dark matter halos (DHs) which evolve to $M_\text{DH}=10^{12}-10^{13}$ M$_\odot$ at $z=2$. The extent of Ly$\alpha$ halos is also found to increase at a rate scaling with the virial radius of growing DHs, $r_\text{vir} \propto M_\text{DH}^{1/3}$($1+z$)$^{-1}$. These increases are consistent with a scenario that the CGM around QSOs evolves in mass and size keeping pace with hosting DHs.

Understanding large-angular-scale galactic foregrounds is crucial for future CMB experiments aiming to detect $B$-mode polarization from primordial gravitational waves. Traditionally, the dust component has been separated using its different frequency dependence. However, using non-CMB observations has potential to increase fidelity and decrease the reconstruction noise. In this exploratory paper we investigate the capability of galactic 21 cm observations to predict the dust foreground in intensity. We train a neural network to predict the dust foreground as measured by the Planck Satellite from the full velocity data-cube of galactic 21 cm emission as measured by the HI4PI survey. We demonstrate that information in the velocity structure clearly improves the predictive power over both a simple integrated emission model and a simple linear model. The improvement is significant at arc-minute scales but more modest at degree scales. This proof of principle on temperature data indicates that it might also be possible to improve foreground polarization templates from the same input data.

This white paper submitted for the 2020 Canadian Long-Range Planning process (LRP2020) presents the prospects for Canada and the Square Kilometre Array (SKA) from 2020-2030, focussing on the first phase of the project (SKA1) scheduled to begin construction early in the next decade. SKA1 will make transformational advances in our understanding of the Universe across a wide range of fields, and Canadians are poised to play leadership roles in several. Canadian key SKA technologies will ensure a good return on capital investment in addition to strong scientific returns, positioning Canadian astronomy for future opportunities well beyond 2030. We therefore advocate for Canada's continued scientific and technological engagement in the SKA from 2020-2030 through participation in the construction and operations phases of SKA1.

The Quijote simulations are a set of 43100 full N-body simulations spanning more than 7000 cosmological models in the $\{\Omega_{\rm m}, \Omega_{\rm b}, h, n_s, \sigma_8, M_\nu, w \}$ hyperplane. At a single redshift the simulations contain more than 8.5 trillions of particles over a combined volume of 43100 $(h^{-1}{\rm Gpc})^3$. Billions of dark matter halos and cosmic voids have been identified in the simulations, whose runs required more than 35 million core hours. The Quijote simulations have been designed for two main purposes: 1) to quantify the information content on cosmological observables, and 2) to provide enough data to train machine learning algorithms. In this paper we describe the simulations and show a few of their applications. We also release the Petabyte of data generated, comprising hundreds of thousands of simulation snapshots at multiple redshifts, halo and void catalogs, together with millions of summary statistics such as power spectra, bispectra, correlation functions, marked power spectra, and estimated probability density functions.

Galaxies evolve from a blue star-forming phase into a red quiescent one by quenching their star formation activity. In high density environments, this galaxy evolution proceeds earlier and more efficiently. Therefore, local galaxy clusters are dominated by well-evolved red, elliptical galaxies. The fraction of blue galaxies in clusters monotonically declines with decreasing redshift, i.e., the Butcher-Oemler effect. In the local Universe, observed blue fractions of massive clusters are as small as $\lesssim$ 0.2. Here we report a discovery of a \lq \lq blue cluster\rq \rq, that is a local galaxy cluster with an unprecedentedly high fraction of blue star-forming galaxies yet hosted by a massive dark matter halo. The blue fraction is 0.57, which is 4.0 $\sigma$ higher than those of the other comparison clusters under the same selection and identification criteria. The velocity dispersion of the member galaxies is 510 km s$^{-1}$, which corresponds to a dark matter halo mass of 2.0$^{+1.9}_{-1.0}\times 10^{14}$ M$_{\odot}$. The blue fraction of the cluster is more than 4.7 $\sigma$ beyond the standard theoretical predictions including semi-analytic models of galaxy formation. The probability to find such a high blue fraction in an individual cluster is only 0.003\%, which challenges the current standard frameworks of the galaxy formation and evolution in the $\Lambda$CDM Universe. The spatial distribution of galaxies around the blue cluster suggests that filamentary cold gas streams can exist in massive halos even in the local Universe. However these cold streams have already disappeared in the theoretically simulated local universes.

Due to ionosphere absorption and the interference by natural and artificial radio emissions, ground observation of the sky at the decameter or longer is very difficult. This unexplored part of electromagnetic spectrum has the potential of great discoveries, notably in the study of cosmic dark ages and dawn, but also in heliophysics and space weather, planets, cosmic ray and neutrinos, pulsar and interstellar medium, extragalactic radio sources, and even SETI. At a forum organized by the International Space Science Institute-Beijing (ISSI-BJ), we discussed the prospect of opening up this window for astronomical observations by using a constellation of small or micro-satellites. We discussed the past experiments and the current ones such as the low frequency payload on Chang'e-4 mission lander, relay satellite and the Longjiang satellite, and also the future DSL mission, which is a linear array on lunar orbit which can make synthesized map of the whole sky as well as measure the global spectrum. We also discuss the synergy with other experiments, including ground global experiments such as EDGES, SARAS, SCI-HI and High-z, PRIZM/Albatros, ground imaging facillities such as LOFAR and MWA, and space experiments such as SUNRISE, DARE/DAPPER and PRATUSH. We also discussed some technical aspects of the DSL concept.

We have carried out a systematic analysis of the nearby (z=0.0279) active galaxy Zw 229.015 using multi-epoch, multi-instrument and deep pointed observations with XMM-Newton, Suzaku, Swift and NuSTAR. Spectral and temporal variability are examined in detail on both the long (weeks-to-years) and short (hours) timescales. A deep Suzaku observation of the source shows two distinct spectral states; a bright-soft state and a dim-hard state in which changes in the power law component account for the differences. Partial covering, blurred reflection and soft Comptonisation models describe the X-ray spectra comparably well, but the smooth, rather featureless, spectrum may be favouring the soft Comptonisation scenario. Moreover, independent of the spectral model, the observed spectral variability is ascribed to the changes in the power law continuum only and do not require changes in the properties of the absorber or blurred reflector incorporated in the other scenarios. The multi-epoch observations between 2009 and 2018 can be described in similar fashion. This could be understood if the primary emission is originating at a large distance from a standard accretion disc or if the disc is optically thin and geometrically thick as recently proposed for Zw 229.015. Our investigation shows that Zw 229.015 behaves similar to sources like Akn 120 and Mrk 530, that exhibit a strong soft-excess, but weak Compton hump and Fe K${\alpha}$ emission.

A linear correlation has been proposed between the CO luminosity ($\rm{L}^{\prime}_{\rm{CO}}$) and full-width at half maximum (FWHM) for high-redshift (z > 1) submillimeter galaxies. However, the controversy concerning the $\rm{L}^{\prime}_{\rm{CO}}$-FWHM correlation seems to have been caused by the use of heterogeneous samples (e.g., different transition lines) and/or data with large measurement uncertainties. In order to avoid the uncertainty caused by using different rotational transitions, in this work we make an extensive effort to select only CO($J = 1-0$) data from the literature. We separate these wide-ranging redshift data into two samples : the low-redshift (z < 1) and high-redshift (z > 1) samples. The samples are corrected for lensing magnification factors if gravitational-lensing effects appeared in the observations. The correlation analysis shows that there exists significant $\rm{L}^{\prime}_{\rm{CO}}$-FWHM correlations for both the low-redshift and high-redshift samples. A comparison of the low- and high-redshift $\rm{L}^{\prime}_{\rm{CO}}$-FWHM correlations does not show strong evolution with redshift. Assuming that there is no evolution, we can use this relation to determine the model independent distances of high-redshift galaxies. We then constrain cosmological models with the calibrated high-redshift CO data and the sample of Type Ia supernovae in the Union 2.1 compilation. In the constraint for wCDM with our samples, the derived values are w_{0} = -1.02 {\pm} 0.17, {\Omega}_{m0} = 0.30{\pm}0.02, and H_{0} = 70.00 {\pm}0.60 km\,s^{-1}\,Mpc^{-1}.

Revealing what fraction of galaxies harbor AGN is central in understanding black hole accretion history of the Universe. However, optical and soft X-ray surveys miss the most highly obscured AGNs. Infrared (IR), instead, is more robust against absorption. Previous IR photometric surveys, however, only had 4 or 5 filters in mid-IR. Our AKARI North Ecliptic Pole (NEP) wide field sample has 18 filters in mid-IR (9 from AKARI, 4 from WISE, and 5 from Spitzer), for the first time, allowing a sophisticated mid-IR SED fitting diagnosis for a statistical number of sources (89178 over 5.4 deg$^2$). By using a SED fitting technique, we investigate the evolution of AGN fraction as a function of redshift and IR (8-1000 $\mu$m) luminosity in an extinction-free way. We found that the AGN fraction (F$_{\rm AGN}$) shows no sign of strong redshift evolution. Instead, F$_{\rm AGN}$ increases with increasing IR luminosity in all redshifts bins ($0<z<2$).

Much of the star formation is obscured by dust. For the complete understanding of the cosmic star formation history (CSFH), infrared (IR) census is indispensable. AKARI carried out deep mid-infrared observations using its continuous 9-band filters in the North Ecliptic Pole (NEP) field (5.4 deg$^2$). This took significant amount of satellite's lifetime, $\sim$10\% of the entire pointed observations. By combining archival Spitzer (5 bands) and WISE (4 bands) mid-IR photometry, we have, in total, 18 band mid-IR photometry, which is the most comprehensive photometric coverage in mid-IR for thousands of galaxies. However previously, we only had shallow optical imaging ($\sim$25.9ABmag) in a small area of 1.0 deg$^2$. As a result, there remained thousands of AKARI's infrared sources undetected in optical. Using the new Hyper Suprime-Cam on Subaru telescope, we obtained deep enough optical images of the entire AKARI NEP field in 5 broad bands ($g\sim$27.5mag). These provided photometric redshift, and thereby IR luminosity for the previously undetected faint AKARI IR sources. Combined with the accurate mid-IR luminosity measurement, we constructed mid-IR LFs, and thereby performed a census of dust-obscured CSFH in the entire AKARI NEP field. We have measured restframe 8$\mu$m, 12$\mu$m luminosity functions (LFs), and estimated total infrared LFs at 0.35$<$z$<$2.2. Our results are consistent with our previous work, but with much reduced statistical errors thanks to the large area coverage of the new data. We have possibly witnessed the turnover of CSFH at $z\sim$2.

Evolutionary properties of infrared (IR) luminous galaxies are important keys to understand dust-obscured star formation history and galaxy evolution. Based on the near- to mid-IR imaging with 9 continuous filters of AKARI space telescope, we present the characteristics of dusty star-forming (SF) galalxies showing polycyclic aromatic hydrocarbon (PAH) features observed by the North Ecliptic Pole (NEP) wide field survey of AKARI and Herschel. All the sample galaxies from the AKARI/NEP-Wide data are selected based both on the Herschel/SPIRE 250 {\mu}m detection and optical spectroscopic redshift data. The physical modelling of spectral energy distribution (SED) using all available data points from u to sub-mm 500 {\mu}m band, including WISE and PACS data where available, takes unique advantages of the continuous near to mid-IR coverage, reliable constraint on far-IR peak, spectroscopically determined accurate redshifts, as well as energy balance principle by MAGPHYS. This enables us to derive physically meaningful and accurate total infrared luminosity and 8 {\mu}m (or PAH) luminosity consistently. Our sample galaxies are in the redshift range z <1, and majority of them appear to be normal SF/spiral populations showing PAH features near the 8 {\mu}m. These SF galaxies showing PAHs in the mid-IR include various types from quiescent to starbursts. Some of our sample show shortage of 8 {\mu}m luminosity compared to the total IR luminosity and this PAH deficit gets severe in more luminous IR galaxies, suggesting PAH molecules in these galaxies destroyed by strong radiation field from SF region or a large amount of cold dust in ISM. The specific SFR of our sample shows mass dependent time evolution which is consistent with downsizing evolutionary pattern.

Using $N$-body simulations with massive neutrino density perturbations, we detect the scale-dependent linear halo bias with high significance. This is the first time that this effect is detected in simulations containing neutrino density perturbations on all scales, confirming the same finding from separate universe simulations. The scale dependence is the result of the additional scale in the system, i.e. the massive neutrino free-streaming length, and it persists even if the bias is defined with respect to the cold dark matter plus baryon (instead of total matter) power spectrum. The separate universe approach provides a good model for the scale-dependent linear bias, and the effect is approximately $0.25f_\nu$ and $0.43f_\nu$ for halos with bias of 1.7 and 3.5, respectively. While the size of the effect is small, it is ${\it not}$ insignificant in terms of $f_\nu$ and should therefore be included to accurately constrain neutrino mass from clustering statistics of biased tracers. More importantly, this feature is a distinct signature of free-streaming particles and cannot be mimicked by other components of the standard cosmological model.

We present the analysis of a new near-infrared (NIR) spectrum of a recently discovered $z=6.621$ quasar PSO J006+39 in an attempt to explore the early growth of supermassive black holes (SMBHs). This NIR (rest-frame ultraviolet, UV) spectrum shows blue continuum slope and rich metal emission lines in addition to Ly$\alpha$ line. We utilize the MgII line width and the rest frame luminosity $L_\text{3000\AA}$ to find the mass of SMBH ($M_\text{BH}$) to be $\sim 10^8 M_\odot$, making this one of the lowest mass quasars at high redshift. The power-law slope index ($\alpha_\lambda$) of the continuum emission is $-2.94\pm0.03$, significantly bluer than the slope of $\alpha_\lambda=-7/3$ predicted from standard thin disc models. We fit the spectral energy distribution (SED) using a model which can fit local SMBHs, which includes warm and hot Comptonisation powered by the accretion flow as well as an outer standard disc. The result shows that the very blue slope is probably produced by a small radial ($\sim230$ gravitational radius, $R_\text{g}$) extent of the standard accretion disc. All plausible SED models require that the source is super-Eddington ($L_\text{bol}/L_\text{Edd} \gtrsim 9$), so the apparently small disc may simply be the inner funnel of a puffed up flow, and clearly the SMBH in this quasar is in a rapid growth phase. We also utilize the rest-frame UV emission lines to probe the chemical abundance in the broad line region (BLR) of this quasar. We find that this quasar has super solar metallicity through photoionization model calculations.

We investigate if the observed small and nearly scale-invariant primordial cosmic perturbation, i.e. the perturbation amplitude $P_\zeta\sim10^{-9}$ and the spectral index $n_s \simeq 0.965$, is typical in the landscape of vacua after imposing anthropic selections on them. We consider the situation where the universe begins from a metastable vacuum driving a precedent inflation, a curvature-dominated open universe is created by tunneling, and the curvature energy is inflated away by new inflation. We argue that the initial inflaton field value is homogeneous but typically non-zero because of the quantum fluctuation of long wavelength modes created during the precedent inflation, and only the universe which accidentally has a small inflaton field value is anthropically selected. We show that this bias, together with certain distributions of inflation model parameters that are physically well-motivated, makes the observed small and nearly scale-invariant spectrum typical.

We discuss the relations between swampland conjectures and observational constraints on both inflation and dark energy. Using the requirement $|\nabla V|\geq c V$, with $c$ as a universal constant whose value can be derived from inflation, there may be no observable distinction between constant and non-constant models of dark energy. However, the latest modification of the above conjecture, which utilizes the second derivative of the potential, opens up the opportunity for observations to determine if the dark energy equation of state deviates from that of a cosmological constant. We also comment on the observability of tensor fluctuations despite the conjecture that field excursions are smaller than the Planck scale.

The role of feedback in triggering or quenching star formation and hence driving galaxy evolution can be directly studied with high resolution integral field observations. The manifestation of feedback in shocks is particularly important to examine in galaxy mergers, where violent interactions of gas takes place in the interstellar medium during the course of the galactic collision. As part of our effort to systematically study the local population of luminous infrared galaxies within the Great Observatories All-Sky LIRG Survey, we undertook the Keck OSIRIS AO LIRG Analysis observing campaign to study the gas dynamics in the inner kiloparsec regions of these systems at spatial scales of a few 10s of parsecs. With high-resolution near-infrared adaptive optics-assisted integral-field observations taken with OSIRIS on the Keck Telescopes, we employ near-infrared diagnostics such as Brg and the ro-vibrationally excited H2 lines to quantify the nuclear star formation rate and identify feedback associated with shocked molecular gas seen in 21 nearby luminous infrared galaxies. Shocked molecular gas is preferentially found in the ultraluminous infrared systems, but may also be triggered at a lower luminosity, earlier merging stage. On circumnuclear scales, AGN have a strong effect on heating the surrounding molecular gas, though their coupling is not simply driven by AGN strength but rather is complicated by orientation, dust shielding, density, and other factors. We find that the nuclear star formation correlates with merger class and diminishing projected nuclear separations. These trends are largely consistent with the picture of merger-induced starbursts within the center of galaxy mergers.

Measurements of the Hubble parameter from the distance ladder are in tension with indirect measurements based on the cosmic microwave background (CMB) data and the inverse distance ladder measurements at 3-4 $\sigma$ level. We consider phenomenological modification to the timing and width of the recombination process and show that they can significantly affect this tension. This possibility is appealing, because such modification affects both the distance to the last scattering surface and the calibration of the baryon acoustic oscillations (BAO) ruler. Moreover, because only a very small fraction of the most energetic photons keep the early universe in the plasma state, it is possible that such modification could occur without affecting the energy density budget of the universe or being incompatible with the very tight limits on the departure from the black-body spectrum of CMB. In particular, we find that under this simplified model, with a conservative subset of Planck data alone, $H_0=73.44_{-6.77}^{+5.50}~{\rm km\ s}^{-1}\ {\rm Mpc}^{-1}$ and in combination with BAO data $H_0=68.86_{-1.35}^{+1.31}~{\rm km\ s}^{-1}\ {\rm Mpc}^{-1}$, decreasing the tension to $\sim 2\sigma$ level. However, when combined with Planck lensing reconstruction and high-$\ell$ polarization data, the tension climbs back to $\sim 2.7\sigma$, despite the uncertainty on non-ladder $H_0$ measurement more than doubling.

It was recently pointed out that the cosmological constant (even metastable one) belongs to the so-called "swampland" and hence cannot be obtained as the low-energy limit of string theory that requires $|\nabla V| > c V$. If true, the dark energy needs to be described by an evolving scalar field, i.e., quintessence with $w>-1$ within supergravity. However, the large hierarchy between the supersymmetry breaking scale and the energy scale of dark energy imposes a challenge on building quintessence models in supergravity as the quintessence field typically acquires a mass of order the gravitino mass. We investigate two approaches to circumvent this obstacle. One is imposing shift symmetry to the quintessence sector, and we demonstrate any quintessence potential can be embedded into supergravity and the fifth force constraint gives little limit on quintessence field displacement, leading to possible observational signature $w>-1$. The structure is stable against quantum corrections. A particular example can address the cosmic coincidence problem. The other approach is sequestered supergravity, and the stability requirement strongly constrains the form of the Kahler potential and superpotential, and the quintessence field displacement is typically much smaller than Planck mass. In addition, to satisfy the fifth force constraint, the quintessence field displacement is further restricted in the sequestered case, requiring $c \ll 1$.

We consider the power spectrum of a biased tracer observed in a finite volume in the presence of a large-scale overdensity and tidal fields. Expanding both the observed power spectrum and the source fields (linear power spectrum, scalar overdensity and tidal field tensor) in spherical harmonics, we explicitly confirm that each $(\ell,m)$ source generates just the corresponding $(\ell,m)$ modes of power spectrum in real space. In redshift space, each $(\ell,m)$ source additionally couples only to $(\ell+2n,m)$ modes of tracer power spectra. This generalizes the Kaiser formula for monopole, quadrupole and hexadecapole of the power spectrum to all $(\ell,m)$ modes generated to the second order in perturbation theory. This formalism can find applications in constraining the super-sample covariance and in the local power spectrum based bispectrum estimators. As an example application, we forecast the ability to measure these modes a survey with BOSS-like galaxy number densities.

The acoustic peaks in the angular power spectrum of cosmic microwave background (CMB) temperature and polarization anisotropies play an important role as a probe of the nature of new relativistic particles contributing to the radiation density in the early universe, parametrized by $\Delta N_{eff}$. The amplitude and phase of the acoustic oscillations provide information about whether the extra species are free-streaming particles, like neutrinos, or tightly-coupled, like the photons, during eras probed by the CMB. On the other hand, some extensions of the Standard Model produce new relativistic particles that decouple from their own non-gravitational interactions after neutrinos, but prior to photons. We study the signature of new relativistic species that decouple during this intermediate epoch. We argue that the decoupling species will cause a scale-dependent change in the amplitude and phase shift of the acoustic oscillations, different from the usual constant shifts on small scales. For intermediate decoupling times, the phase and amplitude shifts depend not only on $\Delta N_{eff}$ but the redshift $z_{dec,X}$ at which the new species decoupled. For $\Delta N_{eff} >0.334$, a Stage IV CMB experiment could determine $N_{eff}$ at the percent level and $z_{dec,X}$ at the $\sim 10\%$ level. For smaller values, $\Delta N_{eff}\sim 0.1$, constraints on $z_{dec,X}$ weaken but remain $\sim 20-50\%$ for $z_{dec,X} \sim \mathcal{O}(10^3-10^4)$. As an application, we study the contributions to $\Delta N_{eff}$ and determine the $z_{dec,X}$ values for simple implementations of the so-called $N$naturalness model.

Cosmic background neutrinos have a large velocity dispersion, which causes the evolution of long-wavelength density perturbations to depend on scale. This scale-dependent growth leads to the well-known suppression in the linear theory matter power spectrum that is used to probe neutrino mass. In this paper, we study the impact of long-wavelength density perturbations on small-scale structure formation. By performing separate universe simulations where the long-wavelength mode is absorbed into the local expansion, we measure the responses of the cold dark matter (CDM) power spectrum and halo mass function, which correspond to the squeezed-limit bispectrum and halo bias. We find that the scale-dependent evolution of the long-wavelength modes causes these quantities to depend on scale and provide simple expressions to model them in terms of scale and the amount of massive neutrinos. Importantly, this scale-dependent bias reduces the suppression in the linear halo power spectrum due to massive neutrinos by 13 and 26% for objects of bias $\bar{b}=2$ and $\bar{b} \gg1$, respectively. We demonstrate with high statistical significance that the scale-dependent halo bias ${\it cannot}$ be modeled by the CDM and neutrino density transfer functions at the time when the halos are identified. This reinforces the importance of the temporal nonlocality of structure formation, especially when the growth is scale dependent.

We present a detailed spectral analysis of the recent 1.5\,Ms XMM-Newton observing campaign on the narrow line Seyfert 1 galaxy IRAS~13224$-$3809, taken simultaneously with 500\,ks of NuSTAR data. The X-ray lightcurve shows three flux peaks, registering at about 100 times the minimum flux seen during the campaign, and rapid variability with a time scale of kiloseconds. The spectra are well fit with a primary powerlaw continuum, two relativistic-blurred reflection components from the inner accretion disk with very high iron abundance, and a simple blackbody-shaped model for the remaining soft excess. The spectral variability is dominated by the power law continuum from a corona region within a few gravitational radii from the black hole. Additionally, blueshifted Ne \textsc{x}, Mg \textsc{xii}, Si \textsc{xiv} and S \textsc{xvi} absorption lines are identified in the stacked low-flux spectrum, confirming the presence of a highly ionized outflow with velocity up to $v= 0.267$ and $0.225$\,c. We fit the absorption features with \texttt{xstar} models and find a relatively constant velocity outflow through the whole observation. Finally, we replace the \texttt{bbody} and supersolar abundance reflection models by fitting the soft excess successfully with the extended reflection model \texttt{relxillD}, which allows for higher densities than the standard \texttt{relxill} model. This returns a disk electron density $n_{\rm e}>10^{18.7}$\,cm$^{-3}$ and lowers the iron abundance from $Z_{\rm Fe}=24^{+3}_{-4}Z_\odot$ with $n_{\rm e}\equiv10^{15}$\,cm$^{-3}$ to $Z_{\rm Fe}=6.6^{+0.8}_{-2.1}Z_\odot$.

We explain the $e^+ e^-$ excess observed by the DAMPE Collaboration using a dark matter model based upon the Higgs triplet model and an additional hidden $SU(2)_X$ gauge symmetry. Two of the $SU(2)_X$ gauge bosons are stable due to a residual discrete symmetry and serve as the dark matter candidate. We search the parameter space for regions that can explain the observed relic abundance, and compute the flux of $e^+ e^-$ coming from a nearby dark matter subhalo. With the inclusion of background cosmic rays, we show that the model can render a good fit to the entire energy spectrum covering the AMS-02, Fermi-LAT and DAMPE data.

The Seyfert 1 galaxy NGC 4593 was monitored spectroscopically with the Hubble Space Telescope as part of a reverberation mapping campaign that also included Swift, Kepler and ground-based photometric monitoring. During 2016 July 12 - August 6, we obtained 26 spectra across a nearly continuous wavelength range of ~1150 - 10,000A. These were combined with Swift data to produce a UV/optical "lag spectrum", which shows the interband lag relative to the Swift UVW2 band as a function of wavelength. The broad shape of the lag spectrum appears to follow the $\tau \propto \lambda^{4/3}$ relation seen previously in photometric interband lag measurements of other active galactic nuclei (AGN). This shape is consistent with the standard thin disk model but the magnitude of the lags implies a disk that is a factor of ~3 larger than predicted, again consistent with what has been previously seen in other AGN. In all cases these large disk sizes, which are also implied by independent gravitational microlensing of higher-mass AGN, cannot be simply reconciled with the standard model. However the most striking feature in this higher resolution lag spectrum is a clear excess around the 3646A Balmer jump. This strongly suggests that diffuse emission from gas in the much larger broad-line region (BLR) must also contribute significantly to the interband lags. While the relative contributions of the disk and BLR cannot be uniquely determined in these initial measurements, it is clear that both will need to be considered in comprehensively modeling and understanding AGN lag spectra.

We investigate the three-point correlation between the Lyman-$\alpha$ forest and the CMB weak lensing ($\delta_F \delta_F \kappa$) expressed as the cross-correlation between the CMB weak lensing field and local variations in the forest power spectrum. In addition to the standard gravitational bispectrum term, we note the existence of a non-standard systematic term coming from mis-estimation of the mean flux over the finite length of Lyman-$\alpha$ skewers. We numerically calculate the angular cross-power spectrum and discuss its features. We integrate it into zero-lag correlation function and compare our predictions with recent results by Doux et al.. We find that our predictions are statistically consistent with the measurement, and including the systematic term improves the agreement with the measurement. We comment on the implication of the response of the Lyman-$\alpha$ forest power spectrum to the long-wavelength density perturbations.

We present results from the optical, ultraviolet and X-ray monitoring of the NLS1 galaxy IRAS 13224-3809 taken with Swift and XMM-Newton during 2016. IRAS 13224-3809 is the most variable bright AGN in the X-ray sky and shows strong X-ray reflection, implying that the X-rays strongly illuminate the inner disc. Therefore, it is a good candidate to study the relationship between coronal X-ray and disc UV emission. However, we find no correlation between the X-ray and UV flux over the available ~40 day monitoring, despite the presence of strong X-ray variability and the variable part of the UV spectrum being consistent with irradiation of a standard thin disc. This means either that the X-ray flux which irradiates the UV emitting outer disc does not correlate with the X-ray flux in our line of sight and/or that another process drives the majority of the UV variability. The former case may be due to changes in coronal geometry, absorption or scattering between the corona and the disc.

We present a public code to generate a mock galaxy catalog in redshift space assuming a log-normal probability density function (PDF) of galaxy and matter density fields. We draw galaxies by Poisson-sampling the log-normal field, and calculate the velocity field from the linearised continuity equation of matter fields, assuming zero vorticity. This procedure yields a PDF of the pairwise velocity fields that is qualitatively similar to that of N-body simulations. We check fidelity of the catalog, showing that the measured two-point correlation function and power spectrum in real space agree with the input precisely. We find that a linear bias relation in the power spectrum does not guarantee a linear bias relation in the density contrasts, leading to a cross-correlation coefficient of matter and galaxies deviating from unity on small scales. We also find that linearising the Jacobian of the real-to-redshift space mapping provides a poor model for the two-point statistics in redshift space. That is, non-linear redshift-space distortion is dominated by non-linearity in the Jacobian. The power spectrum in redshift space shows a damping on small scales that is qualitatively similar to that of the well-known Fingers-of-God (FoG) effect due to random velocities, except that the log-normal mock does not include random velocities. This damping is a consequence of non-linearity in the Jacobian, and thus attributing the damping of the power spectrum solely to FoG, as commonly done in the literature, is misleading.

X-ray reverberation lags from the vicinity of supermassive black holes have been detected in almost 30 AGN. The soft lag, which is the time delay between the hard and soft X-ray light curves, is usually interpreted as the time difference between the direct and reflected emission, but is alternatively suggested to arise from the direct and scattering emission from distant clouds. By analysing the archival Suzaku observations totalling an exposure time of ~ 770 ks, we discover a soft lag of $10\pm3.4$ ks at $9.58\times10^{-6}$ Hz in the luminous quasar PDS 456, which is the longest soft lag and lowest Fourier frequency reported to date. In this study, we use the maximum likelihood method to deal with non-continuous nature of the Suzaku light curves. The result follows the mass-scaling relation for soft lags, which further supports that soft lags originate from the innermost areas of AGN and hence are best interpreted by the reflection scenario. Spectral analysis has been performed in this work and we find no evidence of clumpy partial-covering absorbers. The spectrum can be explained by a self-consistent relativistic reflection model with warm absorbers, and spectral variations over epochs can be accounted for by the change of the continuum, and of column density and ionization states of the warm absorbers.

The squeezed-limit bispectrum, which is generated by nonlinear gravitational evolution as well as inflationary physics, measures the correlation of three wavenumbers, in the configuration where one wavenumber is much smaller than the other two. Since the squeezed-limit bispectrum encodes the impact of a large-scale fluctuation on the small-scale power spectrum, it can be understood as how the small-scale power spectrum "responds" to the large-scale fluctuation. Viewed in this way, the squeezed-limit bispectrum can be calculated using the response approach even in the cases which do not submit to perturbative treatment. To illustrate this point, we apply this approach to the cross-correlation between the large-scale quasar density field and small-scale Lyman-$\alpha$ forest flux power spectrum. In particular, using separate universe simulations which implement changes in the large-scale density, velocity gradient, and primordial power spectrum amplitude, we measure how the Lyman-$\alpha$ forest flux power spectrum responds to the local, long-wavelength quasar overdensity, and equivalently their squeezed-limit bispectrum. We perform a Fisher forecast for the ability of future experiments to constrain local non-Gaussianity using the bispectrum of quasars and the Lyman-$\alpha$ forest. Combining with quasar and Lyman-$\alpha$ forest power spectra to constrain the biases, we find that for DESI the expected $1-\sigma$ constraint is ${\rm err}[f_{\rm NL}]\sim60$. Ability for DESI to measure $f_{\rm NL}$ through this channel is limited primarily by the aliasing and instrumental noise of the Lyman-$\alpha$ forest flux power spectrum. The combination of response approach and separate universe simulations provides a novel technique to explore the constraints from the squeezed-limit bispectrum between different observables.

Modeling the nonlinearity of the halo bispectrum remains a major challenge in modern cosmology, in particular for ongoing and upcoming large-scale structure observations that are performed to study the inflationary physics. The "power spectrum response" offers a solution for bispectrum in the so-called squeezed limit, in which one wavenumber is much smaller than the other two. As a first step, we demonstrate that the halo squeezed-limit bispectrum computed from the second-order standard perturbation theory agrees precisely with the responses of linear halo power spectrum to large-scale density and potential fluctuations. Since the halo power spectrum responses to arbitrarily small scales can straightforwardly be obtained by separate universe simulations, the response approach provides a novel and powerful technique for modeling the nonlinear halo squeezed-limit bispectrum.

Active galactic nucleus (AGN) feedback is the process by which supermassive black holes in the centres of galaxies may moderate the growth of their hosts. Gas outflows from supermassive black holes release huge quantities of energy into the interstellar medium, clearing the surrounding gas. The most extreme of these, the ultra-fast outflows (UFOs), are the subset of X-ray detected outflows with velocities higher than 10,000 km/s, believed to originate in relativistic disc winds, a few hundred gravitational radii from the black hole. The absorption features produced by these outflows are variable, but no clear link has been found between the behaviour of the X-ray continuum and the energy or equivalent width of the outflow features due to the long time-scales of quasar variability. Here, we present the detection of multiple absorption lines from an extreme ultra-fast gas flow in the X-ray spectrum of the active galactic nucleus IRAS 13224-3809, at 0.236+/-0.006 times the speed of light (71,000 km/s), where the absorption is strongly anti-correlated with the emission from the inner regions of the accretion disk. If the gas flow is identified as a genuine outflow then it is in the fastest 5 per cent of such winds, and its variability is hundreds of times faster than in other variable winds, allowing us to observe in hours what would take months in a quasar. We find signatures of the wind simultaneously in both low and high energy detectors, which are consistent with a single ionized outflow, linking the two phenomena. The detection of the wind responding to the emission from the inner disk demonstrates a connection between accretion processes occurring on very different scales, with the X-rays from within a few gravitational radii of the black hole ionizing the relativistically outflowing gas as the flux rises.

By absorbing fluctuations into a local background, separate universe simulations provide a powerful technique to characterize the response of small-scale observables to the long-wavelength density fluctuations, for example those of the power spectrum and halo mass function which lead to the squeezed-limit $n$-point function and halo bias, respectively. Using quintessence dark energy as the paradigmatic example, we extend these simulation techniques to cases where non-gravitational forces in other sectors establish a Jeans scale across which the growth of density fluctuations becomes scale dependent. By characterizing the separate universes with matching background expansion histories, we show that the power spectrum and mass function responses depend on whether the long-wavelength mode is above or below the Jeans scale. Correspondingly, the squeezed bispectrum and halo bias also become scale dependent. Models of bias that are effectively local in the density field at a single epoch, initial or observed, cannot describe this effect which highlights the importance of temporal nonlocality in structure formation. Validated by these quintessence tests, our techniques are applicable to a wide range of models where the complex dynamics of additional fields affect the clustering of matter in the linear regime and it would otherwise be difficult to simulate their impact in the nonlinear regime.

We analyze the latest \emph{Suzaku} observation of the bright neutron star low-mass X-ray binary Serpens X-1 taken in 2013 October and 2014 April. The observation was taken using the burst mode and only suffered mild pile-up effects. A broad iron line is clearly detected in the X-ray spectrum. We test different models and find that the iron line is asymmetric and best interpreted by relativistic reflection. The relativistically broadened iron line is generally believed to originate from the innermost regions of the accretion disk, where strong gravity causes a series of special and general relativistic effects. The iron line profile indicates an inner radius of $\sim8$ $R_{\rm G}$, which gives an upper limit on the size of the neutron star. The asymmetric iron line has been observed in a number of previous observations, which gives several inner radius measurements at different flux states. We find that the inner radius of Serpens X-1 does not evolve significantly over the range of $L/L_{\rm Edd}\sim0.4-0.6$, and the lack of flux dependence of the inner radius implies that the accretion disk may be truncated outside the innermost stable circular orbit by the boundary layer rather than the stellar magnetic field.

The separate universe technique provides a means of establishing consistency relations between short wavelength observables and the long wavelength matter density fluctuations within which they evolve by absorbing the latter into the cosmological background. We extend it to cases where non-gravitational forces introduce a Jeans scale in other species like dynamical dark energy or massive neutrinos. The technique matches the synchronous gauge matter density fluctuations to the local expansion using the acceleration equation and accounts for the temporal nonlocality and scale dependence of the long wavelength response of small scale matter observables, e.g. the nonlinear power spectrum, halo abundance and the implied halo bias, and $N$-point correlation functions. Above the Jeans scale, the local Friedmann equation relates the expansion to real energy densities and a curvature that is constant in comoving coordinates. Below the Jeans scale, the curvature evolves and acts like a fake density component. In all cases, the matter evolution on small scales is correctly modeled as we illustrate using scalar field dark energy with adiabatic or isocurvature initial conditions across the Jeans scale set by its finite sound speed.

Broad Fe K emission lines have been widely observed in the X-ray spectra of black hole systems, and in neutron star systems as well. The intrinsically narrow Fe K fluorescent line is generally believed to be part of the reflection spectrum originating in an illuminated accretion disk, and broadened by strong relativistic effects. However, the nature of the lines in neutron star LMXBs has been under debate. We therefore obtained the longest, high-resolution X-ray spectrum of a neutron star LMXB to date with a 300 ks Chandra HETGS observation of Serpens X-1. The observation was taken under the "continuous clocking" mode and thus free of photon pile-up effects. We carry out a systematic analysis and find that the blurred reflection model fits the Fe line of Serpens X-1 significantly better than a broad Gaussian component does, implying that the relativistic reflection scenario is much preferred. Chandra HETGS also provides highest spectral resolution view of the Fe K region and we find no strong evidence for additional narrow lines.

We present a new observable, position-dependent power spectrum, to measure the large-scale structure bispectrum in the squeezed configuration, where one wavenumber is much smaller than the other two. The squeezed-limit bispectrum measures how the small-scale power spectrum is modulated by a long-wavelength overdensity, which is due to gravitational evolution and possibly inflationary physics. We divide a survey into small subvolumes, compute the local power spectrum and the mean overdensity in each subvolume, and measure the correlation between them. The correlation measures the integral of the bispectrum, which is dominated by squeezed configurations if the scale of the local power spectrum is much smaller than the subvolume size. We use the separate universe approach to model how the small-scale power spectrum is affected by a long-wavelength overdensity gravitationally. This models the nonlinearity of the bispectrum better than the perturbation theory approach. Not only the new observable is easy to interpret, but it sidesteps the complexity of the full bispectrum estimation as both power spectrum and mean overdensity are easier to estimate than the full bispectrum. We report on the first measurement of the position-dependent correlation function from the SDSS-III BOSS DR10 CMASS sample. We detect the bispectrum of the CMASS sample, and constrain their nonlinear bias combining with anisotropic clustering and weak lensing. We finally study the response of the small-scale power spectrum to 1-3 long-wavelength overdensities. We compare the separate universe approach to separate universe simulations to unprecedented accuracy. We test the standard perturbation theory (SPT) hypothesis that the nonlinear n-point function is fully predicted by the linear power spectrum at the same time. We find discrepancies on small scales, which suggest that SPT fails even if it is calculated to all orders.

We report on the first measurement of the three-point function with the position-dependent correlation function from the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 10 CMASS sample. This new observable measures the correlation between two-point functions of galaxy pairs within different subvolumes, $\hat{\xi}({\rm r},{\rm r}_L)$, where ${\rm r}_L$ is the location of a subvolume, and the corresponding mean overdensities, $\bar{\delta}({\rm r}_L)$. This correlation, which we call the "integrated three-point function", $i\zeta(r)=\langle\hat{\xi}({\rm r},{\rm r}_L)\bar{\delta}({\rm r}_L)\rangle$, measures a three-point function of two short- and one long-wavelength modes, and is generated by nonlinear gravitational evolution and possibly also by the physics of inflation. The $i\zeta(r)$ measured from the BOSS data lies within the scatter of those from the mock galaxy catalogs in redshift space, yielding a ten-percent-level determination of the amplitude of $i\zeta(r)$. The tree-level perturbation theory in redshift space predicts how this amplitude depends on the linear and quadratic nonlinear galaxy bias parameters ($b_1$ and $b_2$), as well as on the amplitude and linear growth rate of matter fluctuations ($\sigma_8$ and $f$). Combining $i\zeta(r)$ with the constraints on $b_1\sigma_8$ and $f\sigma_8$ from the global two-point correlation function and that on $\sigma_8$ from the weak lensing signal of BOSS galaxies, we measure $b_2=0.41\pm0.41$ (68% C.L.) assuming standard perturbation theory at the tree level and the local bias model.

We propose a dark matter explanation to simultaneously account for the excess of antiproton-to-proton and positron power spectra observed in the AMS-02 experiment while having the right dark matter relic abundance and satisfying the current direct search bounds. We extend the Higgs triplet model with a hidden gauge symmetry of $SU(2)_X$ that is broken to $Z_3$ by a quadruplet scalar field, rendering the associated gauge bosons stable weakly-interacting massive particle dark matter candidates. By coupling the complex Higgs triplet and the $SU(2)_X$ quadruplet, the dark matter candidates can annihilate into triplet Higgs bosons each of which in turn decays into lepton or gauge boson final states. Such a mechanism gives rise to correct excess of positrons and antiprotons with an appropriate choice of the triplet vacuum expectation value. Besides, the model provides a link between neutrino mass and dark matter phenomenology.

Galaxy proto-clusters at z >~ 2 provide a direct probe of the rapid mass assembly and galaxy growth of present day massive clusters. Because of the need of precise galaxy redshifts for density mapping and the prevalence of star formation before quenching, nearly all the proto-clusters known to date were confirmed by spectroscopy of galaxies with strong emission lines. Therefore, large emission-line galaxy surveys provide an efficient way to identify proto-clusters directly. Here we report the discovery of a large-scale structure at z = 2.44 in the HETDEX Pilot Survey. On a scale of a few tens of Mpc comoving, this structure shows a complex overdensity of Lya emitters (LAE), which coincides with broad-band selected galaxies in the COSMOS/UltraVISTA photometric and zCOSMOS spectroscopic catalogs, as well as overdensities of intergalactic gas revealed in the Lya absorption maps of Lee et al. (2014). We construct mock LAE catalogs to predict the cosmic evolution of this structure. We find that such an overdensity should have already broken away from the Hubble flow, and part of the structure will collapse to form a galaxy cluster with 10^14.5 +- 0.4 M_sun by z = 0. The structure contains a higher median stellar mass of broad-band selected galaxies, a boost of extended Lya nebulae, and a marginal excess of active galactic nuclei relative to the field, supporting a scenario of accelerated galaxy evolution in cluster progenitors. Based on the correlation between galaxy overdensity and the z = 0 descendant halo mass calibrated in the simulation, we predict that several hundred 1.9 < z < 3.5 proto-clusters with z = 0 mass of > 10^14.5 M_sun will be discovered in the 8.5 Gpc^3 of space surveyed by the Hobby Eberly Telescope Dark Energy Experiment.

We show that in a certain, angle-averaged squeezed limit, the $N$-point function of matter is related to the response of the matter power spectrum to a long-wavelength density perturbation, $P^{-1}d^nP(k|\delta_L)/d\delta_L^n|_{\delta_L=0}$, with $n=N-2$. By performing N-body simulations with a homogeneous overdensity superimposed on a flat Friedmann-Robertson-Lema\^itre-Walker (FRLW) universe using the \emph{separate universe} approach, we obtain measurements of the nonlinear matter power spectrum response up to $n=3$, which is equivalent to measuring the fully nonlinear matter $3-$ to $5-$point function in this squeezed limit. The sub-percent to few percent accuracy of those measurements is unprecedented. We then test the hypothesis that nonlinear $N$-point functions at a given time are a function of the linear power spectrum at that time, which is predicted by standard perturbation theory (SPT) and its variants that are based on the ideal pressureless fluid equations. Specifically, we compare the responses computed from the separate universe simulations and simulations with a rescaled initial (linear) power spectrum amplitude. We find discrepancies of 10\% at $k\simeq 0.2 - 0.5 \,h\,{\rm Mpc}^{-1}$ for $5-$ to $3-$point functions at $z=0$. The discrepancy occurs at higher wavenumbers at $z=2$. Thus, SPT and its variants, carried out to arbitrarily high order, are guaranteed to fail to describe matter $N$-point functions ($N>2$) around that scale.

The large-scale statistics of observables such as the galaxy density are chiefly determined by their dependence on the local coarse-grained matter density. This dependence can be measured directly and efficiently in N-body simulations by using the fact that a uniform density perturbation with respect to some fiducial background cosmology is equivalent to modifying the background and including curvature, i.e., by simulating a "separate universe". We derive this mapping to fully non-linear order, and provide a step-by-step description of how to perform and analyse the separate universe simulations. This technique can be applied to a wide range of observables. As an example, we calculate the response of the non-linear matter power spectrum to long-wavelength density perturbations, which corresponds to the angle-averaged squeezed limit of the matter bispectrum and higher $n$-point functions. Using only a modest simulation volume, we obtain results with percent-level precision over a wide range of scales.

We report on the deepest X-ray observation of the narrow-line Seyfert 1 galaxy Mrk 335 in the low-flux state obtained with Suzaku. The data are compared to a 2006 high-flux Suzaku observation when the source was ~10-times brighter. Describing the two flux levels self-consistently with partial covering models would require extreme circumstances, as the source would be subject to negligible absorption during the bright state and ninety-five per cent covering with near Compton-thick material when dim. Blurred reflection from an accretion disc around a nearly maximum spinning black hole (a>0.91, with preference for a spin parameter as high as ~ 0.995) appears more likely and is consistent with the long-term and rapid variability. Measurements of the emissivity profile and spectral modelling indicate the high-flux Suzaku observation of Mrk 335 is consistent with continuum-dominated, jet-like emission (i.e. beamed away from the disc). It can be argued that the ejecta must be confined to within ~25 rg if it does not escape the system. During the low-flux state the corona becomes compact and only extends to about 5 rg from the black hole, and the spectrum becomes reflection-dominated. The low-frequency lags measured at both epochs are comparable indicating that the accretion mechanism is not changing between the two flux levels. Various techniques to study the spectral variability (e.g. principal component analysis, fractional variability, difference spectra, and hardness ratio analysis) indicate that the low-state variability is dominated by changes in the power law flux and photon index, but that changes in the ionisation state of the reflector are also required. Most notably, the ionisation parameter becomes inversely correlated with the reflected flux after a long-duration flare-like event during the observation.

The extreme NLS1 galaxy IRAS 13224-3809 shows significant variability, frequency depended time lags, and strong Fe K line and Fe L features in the long 2011 XMM-Newton observation. In this work we study the spectral properties of IRAS 13224-3809 in detail, and carry out a series of analyses to probe the nature of the source, focusing in particular on the spectral variability exhibited. The RGS spectrum shows no obvious signatures of absorption by partially ionised material (warm absorbers). We fit the 0.3-10.0 keV spectra with a model that includes relativistic reflection from the inner accretion disc, a standard powerlaw AGN continuum, and a low-temperature (~0.1 keV) blackbody, which may originate in the accretion disc, either as direct or reprocessed thermal emission. We find that the reflection model explains the time-averaged spectrum well, and we also undertake flux-resolved and time-resolved spectral analyses, which provide evidence of gravitational light-bending effects. Additionally, the temperature and flux of the blackbody component are found to follow the $L\propto T^{4}$ relation expected for simple thermal blackbody emission from a constant emitting area, indicating a physical origin for this component.

We introduce the light-weight carbon fiber and aluminum gondola designed for the SPIDER balloon-borne telescope. SPIDER is designed to measure the polarization of the Cosmic Microwave Background radiation with unprecedented sensitivity and control of systematics in search of the imprint of inflation: a period of exponential expansion in the early Universe. The requirements of this balloon-borne instrument put tight constrains on the mass budget of the payload. The SPIDER gondola is designed to house the experiment and guarantee its operational and structural integrity during its balloon-borne flight, while using less than 10% of the total mass of the payload. We present a construction method for the gondola based on carbon fiber reinforced polymer tubes with aluminum inserts and aluminum multi-tube joints. We describe the validation of the model through Finite Element Analysis and mechanical tests.

The influence of large-scale density fluctuations on structure formation on small scales is described by the three-point correlation function (bispectrum) in the so-called "squeezed configurations," in which one wavenumber, say $k_3$, is much smaller than the other two, i.e., $k_3\ll k_1\approx k_2$. This bispectrum is generated by non-linear gravitational evolution and possibly also by inflationary physics. In this paper, we use this fact to show that the bispectrum in the squeezed configurations can be measured without employing three-point function estimators. Specifically, we use the "position-dependent power spectrum," i.e., the power spectrum measured in smaller subvolumes of the survey (or simulation box), and correlate it with the mean overdensity of the corresponding subvolume. This correlation directly measures an integral of the bispectrum dominated by the squeezed configurations. Measuring this correlation is only slightly more complex than measuring the power spectrum itself, and sidesteps the considerable complexity of the full bispectrum estimation. We use cosmological $N$-body simulations of collisionless particles with Gaussian initial conditions to show that the measured correlation between the position-dependent power spectrum and the long-wavelength overdensity agrees with the theoretical expectation. The position-dependent power spectrum thus provides a new, efficient, and promising way to measure the squeezed-limit bispectrum from large-scale structure observations such as galaxy redshift surveys.

We construct a model based on an extra gauge symmetry, SU(2)_X x U(1)_{B-L}, which can provide gauge bosons to serve as weakly-interacting massive particle dark matter. The stability of the dark matter is naturally guaranteed by a discrete Z_2 symmetry that is a subgroup of SU(2)_X. The dark matter interacts with standard model fermions by exchanging gauge bosons which are linear combinations of SU(2)_X x U(1)_{B-L} gauge bosons. With the appropriate choice of representation for the new scalar multiplet whose vacuum expectation value spontaneously breaks the SU(2)_X symmetry, the relation between the new gauge boson masses can naturally lead to resonant pair annihilation of the dark matter. After exploring the parameter space of the new gauge couplings subject to constraints from collider data and the observed relic density, we use the results to evaluate the cross section of the dark matter scattering off nucleons and compare it with data from the latest direct detection experiments. We find allowed parameter regions that can be probed by future direct searches for dark matter and LHC searches for new particles.

Survey observations of the three-dimensional locations of galaxies are a powerful approach to measure the distribution of matter in the universe, which can be used to learn about the nature of dark energy, physics of inflation, neutrino masses, etc. A competitive survey, however, requires a large volume (e.g., Vsurvey is roughly 10 Gpc3) to be covered, and thus tends to be expensive. A "sparse sampling" method offers a more affordable solution to this problem: within a survey footprint covering a given survey volume, Vsurvey, we observe only a fraction of the volume. The distribution of observed regions should be chosen such that their separation is smaller than the length scale corresponding to the wavenumber of interest. Then one can recover the power spectrum of galaxies with precision expected for a survey covering a volume of Vsurvey (rather than the volume of the sum of observed regions) with the number density of galaxies given by the total number of observed galaxies divided by Vsurvey (rather than the number density of galaxies within an observed region). We find that regularly-spaced sampling yields an unbiased power spectrum with no window function effect, and deviations from regularly-spaced sampling, which are unavoidable in realistic surveys, introduce calculable window function effects and increase the uncertainties of the recovered power spectrum. While we discuss the sparse sampling method within the context of the forthcoming Hobby-Eberly Telescope Dark Energy Experiment, the method is general and can be applied to other galaxy surveys.

In view of ongoing measurements of the Higgs-like boson at the LHC and direct searches for dark matter, we explore the possibility of accommodating the potential results in a simple new-physics model with discrete gauge symmetry as well as light neutrino masses. Specifically, we study collider and relic-density constraints on the new gauge coupling, predict the cross section of the dark matter scattering off nucleons, and compare it with current direct search data. We also discuss some of the implications if the dark matter is light. The new gauge sector of the model allows it to be compatible with the latest LHC information on the Higgs-like particle and simultaneously satisfy the requirements in its dark matter sector.

We present the first long \emph{Suzaku} observation of the hyperluminous infrared galaxy IRAS 09104+4109 which is dominated by a Type 2 AGN. The infrared to X-ray SED indicates that the source is an obscured quasar with a Compton-thin absorber. However, the 3$\sigma$ hard X-ray detection of the source with the \emph{BeppoSAX} PDS suggested a reflection-dominated, Compton-thick view. The high-energy detection was later found to be possibly contaminated by another Type 2 AGN, NGC 2785, which is only 17 arcmin away. Our new \emph{Suzaku} observation offers simultaneous soft and hard X-ray coverage and excludes contamination from NGC 2785. We find that the hard X-ray component is not detected by the \emph{Suzaku} HXD/PIN (effective energy band 14-45 keV). Both reflection and transmission models have been tested on the latest \emph{Suzaku} and \emph{Chandra} data. The 0.5-10 keV spectrum can be well modelled by the two scenarios. In addition, our analysis implied that the absorption column required in both models is $N_{\rm H} \sim 5 \times 10^{23}$ cm$^{-2}$. Unless IRAS 09104+4109 is a "changing-look" quasar, we confirm that it is a Compton-thin AGN. Although the lack of detection of X-ray emission above 10 keV seems to favour the transmission scenario, we found that the two models offer fairly similar flux predictions over the X-ray band below $\sim$ 40 keV. We also found that the strong iron line shown in the \emph{Suzaku} spectrum is in fact a blend of two emission lines, in which the 6.4 keV one is mostly contributed from the AGN and the 6.7 keV from the hot cluster gas. This implies that the neutral line is perhaps caused by disc reflection, and the reflection-dominated model is more likely the explanation. The transmission model should not be completely ruled out, but a deeper hard X-ray spectrum observation is needed to discriminate between the two scenarios.

We study non-linear massive gravity in the spherically symmetric context. Our main motivation is to investigate the effect of helicity-0 mode which remains elusive after analysis of cosmological perturbation around an open Friedmann-Lemaitre-Robertson-Walker (FLRW) universe. The non-linear form of the effective energy-momentum tensor stemming from the mass term is derived for the spherically symmetric case. Only in the special case where the area of the two sphere is not deviated away from the FLRW universe, the effective energy momentum tensor becomes completely the same as that of cosmological constant. This opens a window for discriminating the non-linear massive gravity from general relativity (GR). Indeed, by further solving these spherically symmetric gravitational equations of motion in vacuum to the linear order, we obtain a solution which has an arbitrary time-dependent parameter. In GR, this parameter is a constant and corresponds to the mass of a star. Our result means that Birkhoff's theorem no longer holds in the non-linear massive gravity and suggests that energy can probably be emitted superluminously (with infinite speed) on the self-accelerating background by the helicity-0 mode, which could be a potential plague of this theory.

Results are presented from a 500ks long XMM-Newton observation of the Narrow-Line Seyfert 1 galaxy IRAS13224-3809. The source is rapidly variable on timescales down to a few 100s. The spectrum shows strong broad Fe-K and L emission features which are interpreted as arising from reflection from the inner parts of an accretion disc around a rapidly spinning black hole. Assuming a power-law emissivity for the reflected flux and that the innermost radius corresponds to the innermost stable circular orbit, the black hole spin is measured to be 0.988 with a statistical precision better than one per cent. Systematic uncertainties are discussed. A soft X-ray lag of 100s confirms this scenario. The bulk of the power-law continuum source is located at a radius of 2-3 gravitational radii.

We consider the prospects to use polarized dark-matter detectors to discriminate between various dark-matter models. If WIMPs are fermions and participate in parity-violating interactions with ordinary matter, then the recoil-direction and recoil-energy distributions of nuclei in detectors will depend on the orientation of the initial nuclear spin with respect to the velocity of the detector through the Galactic halo. If, however, WIMPS are scalars, the only possible polarization-dependent interactions are extremely velocity-suppressed and, therefore, unobservable. Since the amplitude of this polarization modulation is fixed by the detector speed through the halo, in units of the speed of light, exposures several times larger than those of current experiments will be required to be probe this effect.

The XMM-Newton spectrum of the black hole candidate XTE J1652-453 shows a broad and strong Fe K-alpha emission line, generally believed to originate from reflection of the inner accretion disc. These data have been analysed by Hiemstra et al. (2011) using a variety of phenomenological models. We re-examine the spectrum with a self-consistent relativistic reflection model. A narrow absorption line near 7.2 keV may be present, which if real is likely the Fe XXVI absorption line arising from highly ionised, rapidly outflowing disc wind. The blue shift of this feature corresponds to a velocity of about 11100 km/s, which is much larger than the typical values seen in stellar-mass black holes. Given that we also find the source to have a low inclination (i < 32 degrees; close to face-on), we would therefore be seeing through the very base of outflow. This could be a possible explanation for the unusually high velocity. We use a reflection model combined with a relativistic convolution kernel which allows for both prograde and retrograde black hole spin, and treat the potential absorption feature with a physical model for a photo-ionised plasma. In this manner, assuming the disc is not truncated, we could only constrain the spin of the black hole in XTE J1652-453 to be less than ~ 0.5 Jc/GM^{2} at the 90% confidence limit.

We analysed new simultaneous \emph{Suzaku} and \emph{Swift} data of the Broad Line Seyfert 1 (BLS1) galaxy CBS 126. A clear Fe emission line and a strong soft excess are present in the source spectra. We fit the spectra with a relativistic reflection model and find the model tends to fit the data with a high iron abundance possibly due to the large soft excess present. By checking the difference and the RMS spectra, we find there is likely an absorption edge at $\sim$ 0.89 keV, which might be caused by oxygen or neon. We also produced an analysis of the time-resolved spectra in order to examine the existence of the edge. Although high iron abundance is required more in the time-resolved spectra, a model of solar iron abundance together with an absorption edge gives a more physical explanation. Variation of the ionisation parameter is an alternative, plausible explanation for the excess seen in the difference spectra. It is difficult to know if there are warm absorbers in this source from the current data. To further investigate the presence of possible warm absorbers, higher signal-to-noise low-energy data are needed. The work presented here tentatively suggests that the spectra of a BLS1 can be explained by a relativistic reflection model similar to that often seen in their narrow-line cousins.

Certain dark matter interactions with nuclei are mediated possibly by a scalar or pseudoscalar Higgs boson. The estimation of the corresponding cross sections requires a correct evaluation of the couplings between the scalar or pseudoscalar Higgs boson and the nucleons. Progress has been made in two aspects relevant to this study in the past few years. First, recent lattice calculations show that the strange-quark sigma term $\sigma_s$ and the strange-quark content in the nucleon are much smaller than what are expected previously. Second, lattice and model analyses imply sizable SU(3) breaking effects in the determination on the axial-vector coupling constant $g_A^8$ that in turn affect the extraction of the isosinglet coupling $g_A^0$ and the strange quark spin component $\Delta s$ from polarized deep inelastic scattering experiments. Based on these new developments, we re-evaluate the relevant nucleon matrix elements and compute the scalar and pseudoscalar couplings of the proton and neutron. We also find that the strange quark contribution in both types of couplings is smaller than previously thought.

The performance of the Planck instruments in space is enabled by their low operating temperatures, 20K for LFI and 0.1K for HFI, achieved through a combination of passive radiative cooling and three active mechanical coolers. The scientific requirement for very broad frequency coverage led to two detector technologies with widely different temperature and cooling needs. Active coolers could satisfy these needs; a helium cryostat, as used by previous cryogenic space missions (IRAS, COBE, ISO, Spitzer, AKARI), could not. Radiative cooling is provided by three V-groove radiators and a large telescope baffle. The active coolers are a hydrogen sorption cooler (<20K), a 4He Joule-Thomson cooler (4.7K), and a 3He-4He dilution cooler (1.4K and 0.1K). The flight system was at ambient temperature at launch and cooled in space to operating conditions. The HFI bolometer plate reached 93mK on 3 July 2009, 50 days after launch. The solar panel always faces the Sun, shadowing the rest of Planck, andoperates at a mean temperature of 384K. At the other end of the spacecraft, the telescope baffle operates at 42.3K and the telescope primary mirror operates at 35.9K. The temperatures of key parts of the instruments are stabilized by both active and passive methods. Temperature fluctuations are driven by changes in the distance from the Sun, sorption cooler cycling and fluctuations in gas-liquid flow, and fluctuations in cosmic ray flux on the dilution and bolometer plates. These fluctuations do not compromise the science data.

We describe the processing of the 336 billion raw data samples from the High Frequency Instrument (HFI) which we performed to produce six temperature maps from the first 295 days of Planck-HFI survey data. These maps provide an accurate rendition of the sky emission at 100, 143, 217, 353, 545 and 857 GHz with an angular resolution ranging from 9.9 to 4.4^2. The white noise level is around 1.5 {\mu}K degree or less in the 3 main CMB channels (100--217GHz). The photometric accuracy is better than 2% at frequencies between 100 and 353 GHz and around 7% at the two highest frequencies. The maps created by the HFI Data Processing Centre reach our goals in terms of sensitivity, resolution, and photometric accuracy. They are already sufficiently accurate and well-characterised to allow scientific analyses which are presented in an accompanying series of early papers. At this stage, HFI data appears to be of high quality and we expect that with further refinements of the data processing we should be able to achieve, or exceed, the science goals of the Planck project.

All-sky data from the Planck survey and the Meta-Catalogue of X-ray detected Clusters of galaxies (MCXC) are combined to investigate the relationship between the thermal Sunyaev-Zeldovich (SZ) signal and X-ray luminosity. The sample comprises ~ 1600 X-ray clusters with redshifts up to ~ 1 and spans a wide range in X-ray luminosity. The SZ signal is extracted for each object individually, and the statistical significance of the measurement is maximised by averaging the SZ signal in bins of X-ray luminosity, total mass, or redshift. The SZ signal is detected at very high significance over more than two decades in X-ray luminosity (10^43 erg/s < L_500 E(z)^-7/3 < 2 X 10^45 erg/s). The relation between intrinsic SZ signal and X-ray luminosity is investigated and the measured SZ signal is compared to values predicted from X-ray data. Planck measurements and X-ray based predictions are found to be in excellent agreement over the whole explored luminosity range. No significant deviation from standard evolution of the scaling relations is detected. For the first time the intrinsic scatter in the scaling relation between SZ signal and X-ray luminosity is measured and found to be consistent with the one in the luminosity -- mass relation from X-ray studies. There is no evidence of any deficit in SZ signal strength in Planck data relative to expectations from the X-ray properties of clusters, underlining the robustness and consistency of our overall view of intra-cluster medium properties.

(abridged) Planck has observed the entire sky from 30 GHz to 857GHz. The observed foreground emission contains contributions from different phases of the interstellar medium (ISM). We have separated the observed Galactic emission into the different gaseous components (atomic, molecular and ionised) in each of a number of Galactocentric rings. Templates are created for various Galactocentric radii using velocity information from atomic (neutral hydrogen) and molecular (12CO) observations. The ionised template is assumed to be traced by free-free emission as observed by WMAP, while 408 MHz emission is used to trace the synchrotron component. Gas emission not traced by the above templates, namely "ark gas", as evidenced using Planck data, is included as an additional template, the first time such a component has been used in this way. These templates are then correlated with each of the Planck frequency bands, as well as other ancillary data. The emission per column density of the gas templates allows us to create distinct spectral energy distributions (SEDs) per Galactocentric ring and in each of the gaseous tracers from 1.4 GHz to 25 THz (12\mu m). Apart from the thermal dust and free-free emission, we have probed the Galaxy for anomalous (e.g., spinning) dust as well as synchrotron emission. We show that anomalous dust emission is present in the atomic, molecular and dark gas phases throughout the Galactic disk. The derived dust propeties associated with the dark gas phase are derived but do not allow us to reveal the nature of this phase. For all environments, the anomalous emission is consistent with rotation from polycyclic aromatic hydrocarbons (PAHs) and, according to our simple model, accounts for $(25\pm5)%$ (statistical) of the total emission at 30 GHz.

Anomalous microwave emission (AME) has been observed by numerous experiments in the frequency range ~10-60 GHz. Using Planck maps and multi-frequency ancillary data, we have constructed spectra for two known AME regions: the Perseus and Rho Ophiuchi molecular clouds. The spectra are well fitted by a combination of free-free radiation, cosmic microwave background, thermal dust, and electric dipole radiation from small spinning dust grains. The spinning dust spectra are the most precisely measured to date, and show the high frequency side clearly for the first time. The spectra have a peak in the range 20-40 GHz and are detected at high significances of 17.1 sigma for Perseus and 8.4 sigma for Rho Ophiuchi. In Perseus, spinning dust in the dense molecular gas can account for most of the AME; the low density atomic gas appears to play a minor role. In Rho Ophiuchi, the ~30 GHz peak is dominated by dense molecular gas, but there is an indication of an extended tail at frequencies 50-100 GHz, which can be accounted for by irradiated low density atomic gas. The dust parameters are consistent with those derived from other measurements. We have also searched the Planck map at 28.5 GHz for candidate AME regions, by subtracting a simple model of the synchrotron, free-free, and thermal dust. We present spectra for two of the candidates; S140 and S235 are bright HII regions that show evidence for AME, and are well fitted by spinning dust models.

This paper presents the first results of comparison of Planck along with IRAS data with Green Bank Telescope 21-cm observations in 14 fields covering more than 800 square degrees at high Galactic latitude. Galactic dust emission for fields with average HI column density lower than 2 x 10^20 cm^-2 is well correlated with 21-cm emission. The residual emission in these fields, once the HI-correlated emission is removed, is consistent with the expected statistical properties of the cosmic infrared background fluctuations. Fields with larger column densities show significant excess dust emission compared to the HI column density. Regions of excess lie in organized structures that suggest the presence of hydrogen in molecular form, though they are not always correlated with CO emission. Dust emission from intermediate-velocity clouds is detected with high significance. Its spectral properties are consistent with, compared to the local ISM values, significantly hotter dust (T~20 K), lower sub-millimeter dust opacity, and a relative abundance of very small grains to large grains about four times higher. These results are compatible with expectations for clouds that are part of the Galactic fountain in which there is dust shattering and fragmentation. Correlated dust emission in HVCs is not detected; the average of the 99.9% confidence upper limits to the emissivity is 0.15 times the local ISM value at 857 and 3000 GHz, in accordance with gas phase evidence for lower metallicity and depletion in these clouds. Unexpected anti-correlated variations of the dust temperature and emission cross-section per H atom are identified in the local ISM and IVCs, a trend that continues into molecular environments. This suggests that dust growth through aggregation, seen in molecular clouds, is active much earlier in the cloud condensation and star formation processes.

We present the XMM-Newton follow-up for confirmation of Planck cluster candidates. Twenty-five candidates have been observed to date using snapshot (~10 ksec) exposures, ten as part of a pilot programme to sample a low range of signal-to-noise ratios (4<S/N<6), and a further 15 in a programme to observe a sample of S/N>5 candidates. The sensitivity and spatial resolution of XMM-Newton allows unambiguous discrimination between clusters and false candidates. The 4 false candidates have S/N <= 4.1. A total of 21 candidates are confirmed as extended X-ray sources. Seventeen are single clusters, the majority of which are found to have highly irregular and disturbed morphologies (about ~70%). The remaining four sources are multiple systems, including the unexpected discovery of a supercluster at z=0.45. For 20 sources we are able to derive a redshift estimate from the X-ray Fe K line (albeit of variable quality). The new clusters span the redshift range 0.09 <= z <= 0.54, with a median redshift of z~0.37. A first determination is made of their X-ray properties including the characteristic size, which is used to improve the estimate of the SZ Compton parameter, Y_SZ. The follow-up validation programme has helped to optimise the Planck candidate selection process. It has also provided a preview of the X-ray properties of these newly-discovered clusters, allowing comparison with their SZ properties, and to the X-ray and SZ properties of known clusters observed in the Planck survey. Our results suggest that Planck may have started to reveal a non-negligible population of massive dynamically perturbed objects that is under-represented in X-ray surveys. However, despite their particular properties, these new clusters appear to follow the Y_SZ-Y_X relation established for X-ray selected objects, where Y_X is the product of the gas mass and temperature.

We present precise Sunyaev-Zeldovich (SZ) effect measurements in the direction of 62 nearby galaxy clusters (z <0.5) detected at high signal-to-noise in the first Planck all-sky dataset. The sample spans approximately a decade in total mass, 10^14 < M_500 < 10^15, where M_500 is the mass corresponding to a total density contrast of 500. Combining these high quality Planck measurements with deep XMM-Newton X-ray data, we investigate the relations between D_A^2 Y_500, the integrated Compton parameter due to the SZ effect, and the X-ray-derived gas mass M_g,500, temperature T_X, luminosity L_X, SZ signal analogue Y_X,500 = M_g,500 * T_X, and total mass M_500. After correction for the effect of selection bias on the scaling relations, we find results that are in excellent agreement with both X-ray predictions and recently-published ground-based data derived from smaller samples. The present data yield an exceptionally robust, high-quality local reference, and illustrate Planck's unique capabilities for all-sky statistical studies of galaxy clusters.

Planck allows unbiased mapping of Galactic sub-millimetre and millimetre emission from the most diffuse regions to the densest parts of molecular clouds. We present an early analysis of the Taurus molecular complex, on line-of-sight-averaged data and without component separation. The emission spectrum measured by Planck and IRAS can be fitted pixel by pixel using a single modified blackbody. Some systematic residuals are detected at 353 GHz and 143 GHz, with amplitudes around -7 % and +13 %, respectively, indicating that the measured spectra are likely more complex than a simple modified blackbody. Significant positive residuals are also detected in the molecular regions and in the 217 GHz and 100 GHz bands, mainly caused by to the contribution of the J=2-1 and J=1-0 12CO and 13CO emission lines. We derive maps of the dust temperature T, the dust spectral emissivity index beta, and the dust optical depth at 250 microns tau. The temperature map illustrates the cooling of the dust particles in thermal equilibrium with the incident radiation field, from 16-17 K in the diffuse regions to 13-14 K in the dense parts. The distribution of spectral indices is centred at 1.78, with a standard deviation of 0.08 and a systematic error of 0.07. We detect a significant T-beta anti-correlation. The dust optical depth map reveals the spatial distribution of the column density of the molecular complex from the densest molecular regions to the faint diffuse regions. We use near-infrared extinction and HI data at 21-cm to perform a quantitative analysis of the spatial variations of the measured dust optical depth at 250 microns per hydrogen atom tau/NH. We report an increase of tau/NH by a factor of about 2 between the atomic phase and the molecular phase, which has a strong impact on the equilibrium temperature of the dust particles.

Using Planck maps of six regions of low Galactic dust emission with a total area of about 140 square degrees, we determine the angular power spectra of cosmic infrared background (CIB) anisotropies from multipole l = 200 to l = 2000 at 217, 353, 545 and 857 GHz. We use 21-cm observations of HI as a tracer of thermal dust emission to reduce the already low level of Galactic dust emission and use the 143 GHz Planck maps in these fields to clean out cosmic microwave background anisotropies. Both of these cleaning processes are necessary to avoid significant contamination of the CIB signal. We measure correlated CIB structure across frequencies. As expected, the correlation decreases with increasing frequency separation, because the contribution of high-redshift galaxies to CIB anisotropies increases with wavelengths. We find no significant difference between the frequency spectrum of the CIB anisotropies and the CIB mean, with Delta I/I=15% from 217 to 857 GHz. In terms of clustering properties, the Planck data alone rule out the linear scale- and redshift-independent bias model. Non-linear corrections are significant. Consequently, we develop an alternative model that couples a dusty galaxy, parametric evolution model with a simple halo-model approach. It provides an excellent fit to the measured anisotropy angular power spectra and suggests that a different halo occupation distribution is required at each frequency, which is consistent with our expectation that each frequency is dominated by contributions from different redshifts. In our best-fit model, half of the anisotropy power at l=2000 comes from redshifts z<0.8 at 857 GHz and z<1.5 at 545 GHz, while about 90% come from redshifts z>2 at 353 and 217 GHz, respectively.

We present first results on PLCK G266.6-27.3, a galaxy cluster candidate detected at a signal-to-noise ratio of 5 in the Planck All Sky survey. An XMM-Newton validation observation has allowed us to confirm that the candidate is a bona fide galaxy cluster. With these X-ray data we measure an accurate redshift, z = 0.94 +/- 0.02, and estimate the cluster mass to be M_500 = (7.8 +/- 0.8)e+14 solar masses. PLCK G266.6-27.3 is an exceptional system: its luminosity of L_X(0.5-2.0 keV)=(1.4 +/- 0.05)e+45 erg/s, equals that of the two most luminous known clusters in the z > 0.5 universe, and it is one of the most massive clusters at z~1. Moreover, unlike the majority of high-redshift clusters, PLCK G266.6-27.3 appears to be highly relaxed. This observation confirms Planck's capability of detecting high-redshift, high-mass clusters, and opens the way to the systematic study of population evolution in the exponential tail of the mass function.

(Abridged) A brief description of the methodology of construction, contents and usage of the Planck Early Release Compact Source Catalogue (ERCSC), including the Early Cold Cores (ECC) and the Early Sunyaev-Zeldovich (ESZ) cluster catalogue is provided. The catalogue is based on data that consist of mapping the entire sky once and 60% of the sky a second time by Planck, thereby comprising the first high sensitivity radio/submillimetre observations of the entire sky. A Monte-Carlo algorithm based on the injection and extraction of artificial sources into the Planck maps was implemented to select reliable sources among all extracted candidates such that the cumulative reliability of the catalogue is >=90%. The 10sigma photometric flux density limit of the catalogue at |b|>30 deg is 0.49, 1.0, 0.67, 0.5, 0.33, 0.28, 0.25, 0.47 and 0.82 Jy at each of the nine frequencies between 30 and 857 GHz. Sources which are up to a factor of ~2 fainter than this limit, and which are present in "clean" regions of the Galaxy where the sky background due to emission from the interstellar medium is low, are included in the ERCSC if they meet the high reliability criterion. The Planck ERCSC sources have known associations to stars with dust shells, stellar cores, radio galaxies, blazars, infrared luminous galaxies and Galactic interstellar medium features. A significant fraction of unclassified sources are also present in the catalogs. In addition, two early release catalogs that contain 915 cold molecular cloud core candidates and 189 SZ cluster candidates that have been generated using multi-frequency algorithms are presented. The entire source list, with more than 15000 unique sources, is ripe for follow-up characterisation with Herschel, ATCA, VLA, SOFIA, ALMA and other ground-based observing facilities.

We present the first all-sky sample of galaxy clusters detected blindly by the Planck satellite through the Sunyaev-Zeldovich (SZ) effect from its six highest frequencies. This early SZ (ESZ) sample is comprised of 189 candidates, which have a high signal-to-noise ratio ranging from 6 to 29. Its high reliability (purity above 95%) is further ensured by an extensive validation process based on Planck internal quality assessments and by external cross-identification and follow-up observations. Planck provides the first measured SZ signal for about 80% of the 169 previously-known ESZ clusters. Planck furthermore releases 30 new cluster candidates, amongst which 20 meet the ESZ signal-to-noise selection criterion. At the submission date, twelve of the 20 ESZ candidates were confirmed as new clusters, with eleven confirmed using XMM-Newton snapshot observations, most of them with disturbed morphologies and low luminosities. The ESZ clusters are mostly at moderate redshifts (86% with z below 0.3) and span more than a decade in mass, up to the rarest and most massive clusters with masses above 10^15 Msol.

The European Space Agency's Planck satellite was launched on 14 May 2009, and has been surveying the sky stably and continuously since 13 August 2009. Its performance is well in line with expectations, and it will continue to gather scientific data until the end of its cryogenic lifetime. We give an overview of the history of Planck in its first year of operations, and describe some of the key performance aspects of the satellite. This paper is part of a package submitted in conjunction with Planck's Early Release Compact Source Catalogue, the first data product based on Planck to be released publicly. The package describes the scientific performance of the Planck payload, and presents results on a variety of astrophysical topics related to the sources included in the Catalogue, as well as selected topics on diffuse emission.

The spectrum of the bright Seyfert I galaxy MCG-6-30-15 shows a broad Fe K-alpha emission line, generally interpreted to originate very close to the central black hole which must then have a high spin. The observed X-ray variability is driven by a powerlaw component, with little variability found in the iron line energy band. The disconnection of continuum and line variability may be a consequence of strong gravitational light bending. The X-ray spectrum of MCG-6-30-15 is however complex, strongly modified at soft X-ray energies by warm absorbers, which some workers have extended to build an absorption-dominated model of the whole source behaviour. The absorption model interprets the whole X-ray spectrum without a relativistic reflection component and attributes most variability to the warm absorbers. We re-examine the XMM-Newton, Chandra and Suzaku data taken in 2001, 2004 and 2006, respectively, and construct a model consisting of several warm absorbers confirmed by spectral analysis, together with a relativistically blurred reflection component that explains both the soft and hard excesses as well. The model works well on data from all epochs, demonstrating that the reflection model provides a consistent interpretation of the broadband spectrum of MCG-6-30-15.

The integrated Spectral Energy Distributions of the Large and Small Magellanic Cloud appear significantly flatter than expected from dust models based on their FIR and radio emission. The origin of this millimetre excess is still unexplained, and is here investigated using the Planck data. The background CMB contribution is subtracted using an ILC method performed locally around the galaxies. The foreground emission from the Milky Way is subtracted. After subtraction, the emission of both galaxies correlates closely with the gas emission of the LMC and SMC. The millimetre excess in the LMC can be explained by CMB fluctuations, but a significant excess is still present in the SMC SED. The Planck and IRIS data at 100 micron are combined to produce thermal dust temperature and optical depth maps of the two galaxies. The LMC temperature map shows the presence of a warm inner arm already found with the Spitzer data, but also shows the existence of a previously unidentified cold outer arm. Several cold regions are found along this arm, some of which are associated with known molecular clouds. The average emissivity spectral index is found to be consistent with beta=1.5 and beta=1.2 below 500 microns for the LMC and SMC respectively, significantly flatter than the values observed in the Milky Way. Furthermore, there is evidence in the SMC for a further flattening of the SED in the sub-mm. The spatial distribution of the millimetre dust excess in the SMC follows the gas and thermal dust distribution. Different models are explored in order to fit the dust emission in the SMC. It is concluded that the millimetre excess is unlikely to be caused by very cold dust emission and that it could be due to a combination of spinning dust. emission and thermal dust emission by more amorphous dust grains than those present in our Galaxy.

We construct an all-sky map of the apparent temperature and optical depth of thermal dust emission using the Planck-HFI and IRAS data. The optical depth maps are correlated to tracers of the atomic and molecular gas. The correlation is linear in the lowest column density regions at high galactic latitudes. At high NH, the correlation is consistent with that of the lowest NH. In the intermediate NH range, we observe departure from linearity, with the dust optical depth in excess to the correlation. We attribute this excess emission to thermal emission by dust associated with a Dark-Gas phase, undetected in the available HI and CO measurements. We show the 2D spatial distribution of the Dark-Gas in the solar neighborhood and show that it extends around known molecular regions traced by CO. The average dust emissivity in the HI phase in the solar neighborhood follows roughly a power law distribution with beta = 1.8 all the way down to 3 mm, although the SED flattens slightly in the millimetre. The threshold for the existence of the Dark-Gas is found at NH = (8.0\pm 0.58) 10^{20} Hcm-2. Assuming the same dust emissivity at high frequencies for the dust in the atomic and molecular phases leads to an average XCO = (2.54\pm0.13) 10^{20} H2cm-2/(K km s-1). The mass of Dark-Gas is found to be 28% of the atomic gas and 118% of the CO emitting gas in the solar neighborhood. A possible origin for the Dark-Gas is the existence of a dark molecular phase, where H2 survives photodissociation but CO does not. The observed transition for the onset of this phase in the solar neighborhood (AV = 0.4 mag) appears consistent with recent theoretical predictions. We also discuss the possibility that up to half of the Dark-Gas could be in atomic form, due to optical depth effects in the HI measurements.

We present the Sunyaev-Zeldovich (SZ) signal-to-richness scaling relation (Y500-N200) for the MaxBCG cluster catalogue. Employing a multi-frequency matched filter on the Planck sky maps, we measure the SZ signal for each cluster by adapting the filter according to weak-lensing calibrated mass-richness relations (N200-M500). We bin our individual measurements and detect the SZ signal down to the lowest richness systems (N200=10) with high significance, achieving a detection of the SZ signal in systems with mass as low as M500~5e13 Msolar. The observed Y500-N200 relation is well modeled by a power law over the full richness range. It has a lower normalisation at given N200 than predicted based on X-ray models and published mass-richness relations. An X-ray subsample, however, does conform to the predicted scaling, and model predictions do reproduce the relation between our measured bin-average SZ signal and measured bin-average X-ray luminosities. At fixed richness, we find an intrinsic dispersion in the Y500-N200 relation of 60% rising to of order 100% at low richness. Thanks to its all-sky coverage, Planck provides observations for more than 13,000 MaxBCG clusters and an unprecedented SZ/optical data set, extending the list of known cluster scaling laws to include SZ-optical properties. The data set offers essential clues for models of galaxy formation. Moreover, the lower normalisation of the SZ-mass relation implied by the observed SZ-richness scaling has important consequences for cluster physics and cosmological studies with SZ clusters.

The Planck High Frequency Instrument (HFI) is designed to measure the temperature and polarization anisotropies of the Cosmic Microwave Background and galactic foregrounds in six wide bands centered at 100, 143, 217, 353, 545 and 857 GHz at an angular resolution of 10' (100 GHz), 7' (143 GHz), and 5' (217 GHz and higher). HFI has been operating flawlessly since launch on 14 May 2009. The bolometers cooled to 100 mK as planned. The settings of the readout electronics, such as the bolometer bias current, that optimize HFI's noise performance on orbit are nearly the same as the ones chosen during ground testing. Observations of Mars, Jupiter, and Saturn verified both the optical system and the time response of the detection chains. The optical beams are close to predictions from physical optics modeling. The time response of the detection chains is close to pre-launch measurements. The detectors suffer from an unexpected high flux of cosmic rays related to low solar activity. Due to the redundancy of Planck's observations strategy, the removal of a few percent of data contaminated by glitches does not affect significantly the sensitivity. The cosmic rays heat up significantly the bolometer plate and the modulation on periods of days to months of the heat load creates a common drift of all bolometer signals which do not affect the scientific capabilities. Only the high energy cosmic rays showers induce inhomogeneous heating which is a probable source of low frequency noise.

In 2009 Ho\v{r}ava proposed a power-counting renormalizable quantum gravity theory. Afterwards a term in the action that softly violates the detailed balance condition has been considered with the attempt of obtaining a more realistic theory in its IR-limit. This term is proportional to $\omega R^{(3)}$, where $\omega$ is a constant parameter and $R^{(3)}$ is the spatial Ricci scalar. In this paper we derive constraints on this IR-modified Ho\v{r}ava theory using the late-time cosmic accelerating expansion observations. We obtain a lower bound of $|\omega|$ that is nontrivial and depends on $\Lambda_W$, the cosmological constant of the three dimensional spatial action in the Ho\v{r}ava gravity. We find that to preserve the detailed balance condition, one needs to fine-tune $\Lambda_W$ such that $- 2.29\times 10^{-4}< (c^2 \Lambda_W)/(H^2_0 \currentDE) - 2 < 0 $, where $H_0$ and $\currentDE$ are the Hubble parameter and dark energy density fraction in the present epoch, respectively. On the other hand, if we do not insist on the detailed balance condition, then the valid region for $\Lambda_W$ is much relaxed to $-0.39< (c^2 \Lambda_W)/(H^2_0 \currentDE) - 2 < 0.12$. We find that although the detailed balance condition cannot be ruled out, it is strongly disfavored.

We test the truncated disc models using multiwavelength (optical/UV/X-ray) data from the 2005 hard state outburst of the black hole SWIFT J1753.5-0127. This system is both fairly bright and has fairly low interstellar absorption, so gives one of the best datasets to study the weak, cool disc emission in this state. We fit these data using models of an X-ray illuminated disc to constrain the inner disc radius throughout the outburst. Close to the peak, the observed soft X-ray component is consistent with being produced by the inner disc, with its intrinsic emission enhanced in temperature and luminosity by reprocessing of hard X-ray illumination in an overlap region between the disc and corona. This disc emission provides the seed photons for Compton scattering to produce the hard X-ray spectrum, and these hard X-rays also illuminate the outer disc, producing the optical emission by reprocessing. However, the situation is very different as the outburst declines. The optical is probably cyclo-synchrotron radiation, self-generated by the flow, rather than tracing the outer disc. Similarly, limits from reprocessing make it unlikely that the soft X-rays are directly tracing the inner disc radius. This is seen more clearly in a similarly dim low/hard state spectrum from XTE J1118+480. The very small emitting area implied by the relatively high temperature soft X-ray component is completely inconsistent with the much larger, cooler, UV component which is well fit by a truncated disc. We speculate on the origin of this component, but its existence as a clearly separate spectral component from the truncated disc in XTE J1118+480 shows that it does not simply trace the inner disc radius, so cannot constrain the truncated disc models.