The Brans-Dicke theory of gravity is one of the oldest ideas to extend general relativity by introducing a non-minimal coupling between the scalar field and gravity. The Solar System tests put tight constraints on the theory. In order to evade these constraints, various screening mechanisms have been proposed. These screening mechanisms allow the scalar field to couple to matter as strongly as gravity in low density environments while suppressing it in the Solar System. The Vainshtein mechanism, which is found in various modified gravity models such as massive gravity, braneworld models and scalar tensor theories, suppresses the scalar field efficiently in the vicinity of a massive object. This makes it difficult to test these theories from gravitational wave observations. We point out that the recently found scalar gravitational wave memory effect, which is caused by a permanent change in spacetime geometry due to the collapse of a star to a back hole can be significantly enhanced in the Brans-Dicke theory of gravity with the Vainshtein mechanism. This provides a possibility to detect scalar gravitational waves by a network of three or more gravitational wave detectors.

We present a method to introduce relativistic corrections including linear dark energy perturbations in Horndeski theory into Newtonian simulations based on the N-body gauge approach. We assume that standard matter species (cold dark matter, baryons, photons and neutrinos) are only gravitationally-coupled with the scalar field and we then use the fact that one can include modified gravity effects as an effective dark energy fluid in the total energy-momentum tensor. In order to compute the scalar field perturbations, as well as the cosmological background and metric perturbations, we use the Einstein-Boltzmann code \hiclass. As an example, we study the impact of relativistic corrections on the matter power spectrum in k-essence, a subclass of Horndeski theory, including the effects of massless and massive neutrinos. For massive neutrinos with $\sum m_{\nu} = 0.1$ eV, the corrections due to relativistic species (photons, neutrinos and dark energy) can introduce a maximum deviation of approximately $8\%$ to the power spectrum at $k \sim 10^{-3} \ \textrm{Mpc}^{-1}$ at $z=0$, for a scalar field with sound speed $c_{s}^{2}\sim0.013$ during matter domination epoch. Our formalism makes it possible to test beyond \lcdm \ models probed by upcoming large-scale structure surveys on very large scales.

We provide further evidence that a massless cosmological scalar field with a non-minimal coupling to the Ricci curvature of the type $M^2_{\rm pl}(1+\xi \sigma^n/M_{\rm pl}^n) $ alleviates the existing tension between local measurements of the Hubble constant and its inference from CMB anisotropies and baryonic acoustic oscillations data in presence of a cosmological constant. In these models, the expansion history is modified compared to $\Lambda$CDM at early time, mimicking a change in the effective number of relativistic species, and gravity weakens after matter-radiation equality. Compared to $\Lambda$CDM, a quadratic ($n=2$) coupling increases the Hubble constant when {\em Planck} 2018 (alone or in combination with BAO and SH0ES) measurements data are used in the analysis. Negative values of the coupling, for which the scalar field decreases, seem favored and consistency with Solar System can be naturally achieved for a large portion of the parameter space without the need of any screening mechanism. We show that our results are robust to the choice of $n$, also presenting the analysis for $n=4$.

We propose new massive gravity theories with 5 dynamical degrees of freedom. We evade uniqueness theorems regarding the form of the kinetic and potential terms by adopting the "generalized massive gravity" framework, where a global translation invariance is broken. By exploiting the rotation symmetry in the field space, we determine two novel classes of theories. The first one is an extension of generalized massive gravity with a non-minimal coupling. On the other hand, the second theory produces a mass term that is different from de Rham, Gabadadze, Tolley construction and trivially has 5 degrees of freedom. Both theories allows for stable cosmological solutions without infinite strong coupling, which are free of ghost and gradient instabilities.

One of the most appealing approaches to ease the Hubble tension is the inclusion of an early dark energy (EDE) component that adds energy to the Universe in a narrow redshift window around the time of recombination and dilutes faster than radiation afterwards. In this paper, we analyze EDE in the framework of $\alpha$-attractor models. As well known, the success in alleviating the Hubble tension crucially depends on the shape of the energy injection. We show how different types of energy injections can be easily obtained, thanks to the freedom in choosing the functional form of the potential inspired by $\alpha$-attractor models. To confirm our intuition we perform an MCMC analysis for three representative cases and find indeed that $H_0$ is significantly larger than in $\Lambda$CDM like in other EDE models. Unlike axion-driven EDE models with super Planckian decay constant, the curvature of EDE models required by the data is natural in the context of recent theoretical developments in $\alpha$-attractors.

Generative adversarial networks (GANs) have been recently applied as a novel emulation technique for large scale structure simulations. Recent results show that GANs can be used as a fast, efficient and computationally cheap emulator for producing novel weak lensing convergence maps as well as cosmic web data in 2-D and 3-D. However, like any algorithm, the GAN approach comes with a set of limitations, such as an unstable training procedure and the inherent randomness of the produced outputs. In this work we employ a number of techniques commonly used in the machine learning literature to address the mentioned limitations. In particular, we train a GAN to produce both weak lensing convergence maps and dark matter overdensity field data for multiple redshifts, cosmological parameters and modified gravity models. In addition, we train a GAN using the newest Illustris data to emulate dark matter, gas and internal energy distribution data simultaneously. Finally, we apply the technique of latent space interpolation to control which outputs the algorithm produces. Our results indicate a 1-20% difference between the power spectra of the GAN-produced and the training data samples depending on the dataset used and whether Gaussian smoothing was applied. Finally, recent research on generative models suggests that such algorithms can be treated as mappings from a lower-dimensional input (latent) space to a higher dimensional (data) manifold. We explore such a theoretical description as a tool for better understanding the latent space interpolation procedure.

We present a thorough analysis of the cosmological perturbations in Generalised Massive Gravity. This is an extension of de Rham-Gabadadze-Tolley theory where the translation invariance in the St\"uckelberg field space is broken. This allows the mass parameters to be promoted to functions of the St\"uckelberg fields. We consider an exact cosmological background in this theory and study the stability of perturbations. We derive conditions to avoid ghost, gradient and tachyonic instability. The cosmology is an extension of the self-accelerating branch of the constant mass parameter theory, but now all five massive graviton polarisations propagate. For concreteness, we consider a minimal version of the theory where cosmology undergoes an accelerated expansion at late times and show that the perturbative stability is preserved for a range of parameters.

We derive the one-loop perturbative formula of the redshift-space matter power spectrum after density field reconstruction in the Zeldovich approximation. We find that the reconstruction reduces the amplitudes of nonlinear one-loop perturbative terms significantly by partially erasing the nonlinear mode-coupling between density and velocity fields. In comparison with N-body simulations, we find that both the monopole and quadrupole spectra of reconstructed matter density fields agree with the one-loop perturbation theory up to higher wavenumber than those before reconstruction. We also evaluate the impact on cosmic growth rate assuming the survey volume and the number density like the Baryon Oscillation Spectroscopic Survey and find that the total error, including statistical and systematic ones due to one-loop approximation, decreases by half.

We study relativistic stars in Hordenski theories that evade the gravitational wave constraints and exhibit the Vainshtein mechanism, focusing on a model based on the cubic Galileon Lagrangian. We derive the scalar field profile for static spherically symmetric objects in asymptotically de Sitter space-time with a linear time dependence. The exterior solution matches to the black hole solution found in the literature. Due to the Vainshtein mechanism, the stellar structure is indistinguishable from that of General Relativity with the same central density as long as the radius of the star is shorter than the Vainshtein radius. On the other hand, the scalar field is not suppressed beyond the Vainshtein radius. These solutions have an additional integration constant in addition to the mass of the star.

Accurate predictions for the non-linear matter power spectrum are needed to confront theory with observations in current and near future weak lensing and galaxy clustering surveys. We propose a computationally cheap method to create an emulator for modified gravity models by utilizing existing emulators for $\Lambda$CDM. Using a suite of $N$-body simulations we construct a fitting function for the enhancement of both the linear and non-linear matter power spectrum in the commonly studied Hu-Sawicki $f(R)$ gravity model valid for wave-numbers $k \lesssim 5-10\, h\text{Mpc}^{-1}$ and redshifts $z \lesssim 3$. We show that the cosmology dependence of this enhancement is relatively weak so that our fit, using simulations coming from only one cosmology, can be used to get accurate predictions for other cosmological parameters. We also show that the cosmology dependence can, if needed, be included by using linear theory, approximate $N$-body simulations (such as COLA) and semi-analytical tools like the halo model. Our final fit can easily be combined with any emulator or semi-analytical models for the non-linear $\Lambda$CDM power spectrum to accurately, and quickly, produce a non-linear power spectrum for this particular modified gravity model. The method we use can be applied to fairly cheaply construct an emulator for other modified gravity models. As an application of our fitting formula we use it to compute Fisher-forecasts for how well galaxy clustering and weak lensing in a Euclid-like survey will be at constraining modifications of gravity.

We present a systematic study of galaxy bias in the presence of primordial non-Gaussianity in General Relativity (GR) at second order in perturbation theory. The non-linearity of the Poisson equation in GR and primordial non-Gaussianity are consistently included. We show that the inclusion of non-local primordial non-Gaussianity in addition to local non-Gaussianity is important to show the absence of the modulation of small scale clustering by the long-wavelength mode in the single field slow-roll inflation. We study the bispectrum of the relativistic galaxy density in several gauges and identify the effect of primordial non-Gaussianity and GR corrections.

We develop perturbation theory approaches to model the marked correlation function constructed to up-weight low density regions of the Universe, which might help distinguish modified gravity models from the standard cosmology model based on general relativity. Working within Convolution Lagrangian Perturbation Theory, we obtain that weighted correlation functions are expressible as double convolution integrals that we approximate using a combination of Eulerian and Lagrangian schemes. We find that different approaches agree within 1$\%$ on quasi non-linear scales. Compared with {\it N}-body simulations, the perturbation theory is found to provide accurate predictions for the marked correlation function of dark matter fields, dark matter halos as well as Halo Occupation Distribution galaxies down to $30$ Mpc/h. These analytic approaches help to understand the degeneracy between the mark and the galaxy bias and find a way to maximize the differences among various cosmological models.

We provide in depth MCMC comparisons of two different models for the halo redshift space power spectrum, namely a variant of the commonly applied Taruya-Nishimichi-Saito (TNS) model and an effective field theory of large scale structure (EFTofLSS) inspired model. Using many simulation realisations and Stage IV survey-like specifications for the covariance matrix, we check each model's range of validity by testing for bias in the recovery of the fiducial growth rate of structure formation. The robustness of the determined range of validity is then tested by performing additional MCMC analyses using higher order multipoles, a larger survey volume and a more highly biased tracer catalogue. We find that under all tests, the TNS model's range of validity remains robust and is found to be much higher than previous estimates. The EFTofLSS model fails to capture the spectra for highly biased tracers as well as becoming biased at higher wavenumbers when considering a very large survey volume. Further, we find that the marginalised constraints on $f$ for all analyses are stronger when using the TNS model.

We introduce The Novel Probes Project, an initiative to advance the field of astrophysical tests of the dark sector by creating a forum that connects observers and theorists. This review focuses on tests of gravity and is intended to be of use primarily to observers, but also to theorists with interest in the development of experimental tests. It is twinned with a separate review on tests of dark matter self-interactions (Adhikari et al., in prep.). Our focus is on astrophysical probes of gravity in the weak-field regime, ranging from stars to quasilinear cosmological scales. These are complementary to both strong-field tests and background and linear probes in cosmology. In particular, the nonlinear screening mechanisms that are an integral part of viable modified gravity models lead to characteristic signals specifically on astrophysical scales. The constraining power of these signals is not limited by cosmic variance, but comes with the challenge of building robust theoretical models of the nonlinear dynamics of stars, galaxies, clusters and large scale structure. In this review we lay the groundwork for a thorough exploration of the astrophysical regime with an eye to using the current and next generation of observations for tests of gravity. We begin by setting the scene for how theories beyond General Relativity are expected to behave, focusing primarily on screened fifth forces. We describe the analytic and numerical techniques for exploring the pertinent astrophysical systems, as well as the signatures of modified gravity. With these in hand we present a range of observational tests, and discuss prospects for future measurements and theoretical developments.

There is a well known degeneracy between the enhancement of the growth of large-scale structure produced by modified gravity models and the suppression due to the free-streaming of massive neutrinos at late times. This makes the matter power-spectrum alone a poor probe to distinguish between modified gravity and the concordance $\Lambda$CDM model when neutrino masses are not strongly constrained. In this work, we investigate the potential of using redshift-space distortions (RSD) to break this degeneracy when the modification to gravity is scale-dependent in the form of Hu-Sawicki $f(R)$. We find that if the linear growth rate can be recovered from the RSD signal, the degeneracy can be broken at the level of the dark matter field. However, this requires accurate modelling of the non-linearities in the RSD signal, and we here present an extension of the standard perturbation theory-based model for non-linear RSD that includes both Hu-Sawicki $f(R)$ modified gravity and massive neutrinos.

Verlinde's theory of Emergent Gravity (EG) describes gravity as an emergent phenomenon rather than a fundamental force. Applying this reasoning in de Sitter space leads to gravity behaving differently on galaxy and galaxy cluster scales; this excess gravity might offer an alternative to dark matter. Here we test these ideas using the data from the Coma cluster and from 58 stacked galaxy clusters. The X-ray surface brightness measurements of the clusters at $0.1 < z < 1.2$ along with the weak lensing data are used to test the theory. We find that the simultaneous EG fits of the X-ray and weak lensing datasets are significantly worse than those provided by General Relativity (with cold dark matter). For the Coma cluster, the predictions from Emergent Gravity and General Relativity agree in the range of 250 - 700 kpc, while at around 1 Mpc scales, EG total mass predictions are larger by a factor of 2. For the cluster stack the predictions are only in good agreement at around the 1 - 2 Mpc scales, while for $r \gtrsim 10$ Mpc EG is in strong tension with the data. According to the Bayesian information criterion analysis, GR is preferred in all tested datasets; however, we also discuss possible modifications of EG that greatly relax the tension with the data.

The local galaxy bias formalism relies on the energy constraint equation at the formation time to relate the metric perturbation to the matter density contrast. In the Newtonian approximation, this relationship is linear, which allows us to specify the initial galaxy density as a function of local physical operators. In general relativity however, the relationship is intrinsically nonlinear and a modulation of the short-wavelength mode by the long-wavelength mode might be expected. We describe in detail how to obtain local coordinates where the coupling of the long- to the short-wavelength modes is removed through a change of coordinates (in the absence of the primordial non-Gaussianity). We derive the general-relativistic correction to the galaxy bias expansion at second order. The correction does not come from the modulation of small-scale clustering by the long-wavelength mode; instead, it arises from distortions of the volume element by the long-wavelength mode and it does not lead to new bias parameters.

We revisit the question of viability of bigravity cosmology as a candidate for dark energy. In the context of the low energy limit model, where matter couples to a single metric, we study linear perturbations around homogeneous and isotropic backgrounds to derive the Poisson's equation for the Newtonian potential. Extending to second order perturbations, we identify the Vainshtein radius below which non-linear scalar self interactions conspire to reproduce GR on local scales. We combine all of these results to determine the parameter space that allows a late time de-Sitter attractor compatible with observations and a successful Vainsthein mechanism. We find that the requirement on having a successful Vainsthein mechanism is not compatible with the existence of cosmological solutions at early times.

A recent article "An X-ray chimney extending hundreds of parsecs above and below the Galactic Centre" (2019, Nature, 567, 34) reported the detection of chimney-like X-ray-emitting features above and below the Galactic Center from XMM-Newton observations. We note here that these features were already reported by our Suzaku papers: Nakashima et al. (2013, ApJ, 773, 20, arXiv:1310.4236) for the southern feature and Nakashima et al. (2019, ApJ, in press, arXiv:1903.02571) for the northern feature. In particular, Nakashima et al. (2013) show that the ionization state of the southern feature is not in collisional ionization equilibrium and is most likely in a recombining or over-ionized state, which suggests its origin in the Galactic Center about 0.1 Myr ago.

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale.

We report the discovery of a $1^\circ$ scale X-ray plume in the northern Galactic Center (GC) region observed with Suzaku. The plume is located at ($l$, $b$) $\sim$ ($0\mbox{$.\!\!^\circ$}2$, $0\mbox{$.\!\!^\circ$}6$), east of the radio lobe reported by previous studies. No significant X-ray excesses are found inside or to the west of the radio lobe. The spectrum of the plume exhibits strong emission lines from highly ionized Mg, Si, and S that is reproduced by a thin thermal plasma model with $kT \sim 0.7$ keV and solar metallicity. There is no signature of non-equilibrium ionization. The unabsorbed surface brightness is $3\times10^{-14}$ erg cm$^{-2}$ s$^{-1}$ arcmin$^{-2}$ in the 1.5-3.0 keV band. Strong interstellar absorption in the soft X-ray band indicates that the plume is not a foreground source but is at the GC distance, giving a physical size of $\sim$100 pc, a density of 0.1 cm$^{-3}$, thermal pressure of $1\times10^{-10}$ erg cm$^{-3}$, mass of 600 $M_\odot$ and thermal energy of $7\times10^{50}$ erg. From the apparent association with a polarized radio emission, we propose that the X-ray plume is a magnetized hot gas outflow from the GC.

This paper reports the analysis procedure and results of simultaneous spectral fits of the Suzaku archive data for Sagittarius (Sgr) A East and the nearby Galactic center X-ray emission (GCXE). The results are that the mixed-morphology supernova remnant Sgr A East has a recombining plasma (RP) with Cr and Mn He$\alpha$ lines, and a power-law component (PL) with an Fe I K$\alpha$ line. The nearby GCXE has a $\sim$1.5-times larger surface brightness than the mean GCXE far from Sgr A East, although the spectral shape is almost identical. Based on these results, we interpret that the origins of the RP and the PL with the Fe I K$\alpha$ line are past big flares of Sgr A$^*$.

In quadratic-order degenerate higher-order scalar-tensor (DHOST) theories compatible with gravitational-wave constraints, we derive the most general Lagrangian allowing for tracker solutions characterized by $\dot{\phi}/H^p={\rm constant}$, where $\dot{\phi}$ is the time derivative of a scalar field $\phi$, $H$ is the Hubble expansion rate, and $p$ is a constant. While the tracker is present up to the cubic-order Horndeski Lagrangian $L=c_2X-c_3X^{(p-1)/(2p)} \square \phi$, where $c_2, c_3$ are constants and $X$ is the kinetic energy of $\phi$, the DHOST interaction breaks this structure for $p \neq 1$. Even in the latter case, however, there exists an approximate tracker solution in the early cosmological epoch with the nearly constant field equation of state $w_{\phi}=-1-2p\dot{H}/(3H^2)$. The scaling solution, which corresponds to $p=1$, is the unique case in which all the terms in the field density $\rho_{\phi}$ and the pressure $P_{\phi}$ obey the scaling relation $\rho_{\phi} \propto P_{\phi} \propto H^2$. Extending the analysis to the coupled DHOST theories with the field-dependent coupling $Q(\phi)$ between the scalar field and matter, we show that the scaling solution exists for $Q(\phi)=1/(\mu_1 \phi+\mu_2)$, where $\mu_1$ and $\mu_2$ are constants. For the constant $Q$, i.e., $\mu_1=0$, we derive fixed points of the dynamical system by using the general Lagrangian with scaling solutions. This result can be applied to the model construction of late-time cosmic acceleration preceded by the scaling $\phi$-matter-dominated epoch.

Ongoing and upcoming cosmological surveys will significantly improve our ability to probe the equation of state of dark energy, $w_{\rm DE}$, and the phenomenology of Large Scale Structure. They will allow us to constrain deviations from the $\Lambda$CDM predictions for the relations between the matter density contrast and the weak lensing and the Newtonian potential, described by the functions $\Sigma$ and $\mu$, respectively. In this work, we derive the theoretical prior for the joint covariance of $w_{\rm DE}$, $\Sigma$ and $\mu$, expected in general scalar-tensor theories with second order equations of motion (Horndeski gravity), focusing on their time-dependence at certain representative scales. We employ Monte-Carlo methods to generate large ensembles of statistically independent Horndeski models, focusing on those that are physically viable and in broad agreement with local tests of gravity, the observed cosmic expansion history and the measurement of the speed of gravitational waves from a binary neutron star merger. We identify several interesting features and trends in the distribution functions of $w_{\rm DE}$, $\Sigma$ and $\mu$, as well as in their covariances; we confirm the high degree of correlation between $\Sigma$ and $\mu$ in scalar-tensor theories. The derived prior covariance matrices will allow us to reconstruct jointly $w_{\rm DE}$, $\Sigma$ and $\mu$ in a non-parametric way.

We show how Newtonian cosmological simulations can be employed to investigate the non-linear evolution of two particle species in a relativistic context. We discuss the application for massive neutrinos and other multi-species systems such as Cold Dark Matter (CDM) plus baryons or Warm Dark Matter (WDM). We propose a method that allows us to perform simulations including massive neutrinos and general relativistic effects at almost the same computational cost as ordinary CDM only N-body simulations, employing tailor-made initial conditions and a dictionary for the interpretation of the simulation output.

Most of young and middle-aged supernova remnants (SNRs) exhibit an ionizing plasma (IP), an ionizing process following a shock heated SNR gas. On the other hand, significant fractions of SNRs exhibit a recombining plasma (RP). The origin and the mechanisms of the RP, however, are not yet well understood. This paper proposes a new model that the RP is followed after the IP process taken at the first epoch of the SNR evolution. Using the high quality and wide band (0.6-10 keV) spectrum of IC 443, we nicely fitted with a model of two RP (two-RP model) plus a power law (PL) with an Fe I Kalpha line component. The ionization temperature in one RP monotonously increases from Ne-Ca, while that in the other RP shows a drastic increase from Cr-Ni. Origin and mechanism of the two-RP and PL with an Fe I Kalpha line components are possibly due to a different evolution of two plasmas and ionization by the low-energy cosmic ray.

In the context of Degenerate Higher-Order Scalar-Tensor (DHOST) theories, we study cosmological solutions and their stability properties. In particular, we explicitly illustrate the crucial role of degeneracy by showing how the higher order homogeneous equations in the physical frame (where matter is minimally coupled) can be recast in a system of equations that do not involve higher order derivatives. We study the fixed points of the dynamics, finding the conditions for having a de Sitter attractor at late times. Then we consider the coupling to matter field (described for convenience by a k-essence Lagrangian) and find the conditions to avoid gradient and ghost instabilities at linear order in cosmological perturbations, extending previous work. Finally, we apply these results to a simple subclass of DHOST theories, showing that de Sitter attractor conditions, no ghost and no gradient instabilities conditions (both in the self-accelerating era and in the matter dominated era) can be compatible.

We present the results from the Hitomi Soft Gamma-ray Detector (SGD) observation of the Crab nebula. The main part of SGD is a Compton camera, which in addition to being a spectrometer, is capable of measuring polarization of gamma-ray photons. The Crab nebula is one of the brightest X-ray / gamma-ray sources on the sky, and, the only source from which polarized X-ray photons have been detected. SGD observed the Crab nebula during the initial test observation phase of Hitomi. We performed the data analysis of the SGD observation, the SGD background estimation and the SGD Monte Carlo simulations, and, successfully detected polarized gamma-ray emission from the Crab nebula with only about 5 ks exposure time. The obtained polarization fraction of the phase-integrated Crab emission (sum of pulsar and nebula emissions) is (22.1 $\pm$ 10.6)% and, the polarization angle is 110.7$^o$ + 13.2 / $-$13.0$^o$ in the energy range of 60--160 keV (The errors correspond to the 1 sigma deviation). The confidence level of the polarization detection was 99.3%. The polarization angle measured by SGD is about one sigma deviation with the projected spin axis of the pulsar, 124.0$^o$ $\pm$0.1$^o$.

Using the consistency relation in Fourier space, we derive the observed galaxy bispectrum from single-field inflation in the squeezed limit, in which one of the three modes has a wavelength much longer than the other two. This provides a non-trivial check of the full computation of the bispectrum based on second-order cosmological perturbation theory in this limit. We show that gauge modes need to be carefully removed in the second-order cosmological perturbations in order to calculate the observed galaxy bispectrum in the squeezed limit. We then give an estimate of the effective non-Gaussianity due to general relativistic lightcone effects that could mimic a primordial non-Gaussian signal.

We study the scalar-tensor-tensor non-Gaussian signal in an inflationary model comprising also an axion coupled with SU(2) gauge fields. In this set-up, metric fluctuations are sourced by the gauge fields already at the linear level providing an enhanced chiral gravitational waves spectrum. The same mechanism is at work in generating an amplitude for the three-point function that is parametrically larger than in standard single-field inflation.

Einstein's theory of gravity, General Relativity, has been precisely tested on Solar System scales, but the long-range nature of gravity is still poorly constrained. The nearby strong gravitational lens, ESO 325-G004, provides a laboratory to probe the weak-field regime of gravity and measure the spatial curvature generated per unit mass, $\gamma$. By reconstructing the observed light profile of the lensed arcs and the observed spatially resolved stellar kinematics with a single self-consistent model, we conclude that $\gamma = 0.97 \pm 0.09$ at 68% confidence. Our result is consistent with the prediction of 1 from General Relativity and provides a strong extragalactic constraint on the weak-field metric of gravity.

The Galactic Center X-ray Emission (GCXE) is composed of high temperature (7 keV) and low temperature (1 keV) plasmas (HTP and LTP, respectively). The global structure of the HTP is roughly uniform over the Galactic center (GC) region, and the origin of the HTP has been extensively studied. On the other hand, the LTP is more clumpy, and the origin has not been studied in detail. In the S XV He alpha line map, a pair of horn-like soft diffuse sources are seen at the symmetric positions with respect to Sagittarius A^{star}. The X-ray spectra of the pair are well represented by an absorbed thin thermal plasma model of a temperature and N_H of 0.6-0.7 keV and 4x10^{22} cm^{-2}, respectively. The N_H values indicate that the pair are located near at the GC. Then the dynamical time scales of the pair are 10^{5} yr. The Si and S abundances and the surface brightnesses in the S XV He alpha line band are 0.7-1.2 and 0.6-1.3 solar, and (2.0-2.4)x10^{-15} erg s^{-1} cm^{-2} arcmin^{-2}, respectively. The temperature, abundances, and surface brightness are similar to those of the LTP in the GCXE, while the abundances are far larger than those of known point sources, typically coronal active stars and RS CVn-type active binaries. Based on these results, possible origin of the LTP is discussed.

We present an X-ray study of the GeV gamma-ray supernova remnant (SNR) HB 21 with Suzaku. HB 21 is interacting with molecular clouds and the faintest in the GeV band among known GeV SNRs. We discovered strong radiative recombination continua of Si and S from the center of the remnant, which provide the direct evidence of a recombining plasma (RP). The total emission can be explained with the RP and ionizing plasma components. The electron temperature and recombination timescale of the RP component were estimated as 0.17 (0.15-0.18) keV and 3.2 (2.0-4.8) $\times$ 10$^{11}$ s cm$^{-3}$, respectively. The estimated age of the RP (RP age; $\sim$ 170 kyr) is the longest among known recombining GeV SNRs, because of very low density of electrons ($\sim$ 0.05 cm$^{-3}$). We have examined dependencies of GeV spectral indices on each of RP ages and SNR diameters for nine recombining GeV SNRs. Both showed possible positive correlations, indicating that both the parameters can be good indicators of properties of accelerated protons, for instance, degree of escape from the SNR shocks. A possible scenario for a process of proton escape is introduced; interaction with molecular clouds makes weaker magnetic turbulence and cosmic-ray protons escape, simultaneously cooling down the thermal electrons and generate an RP.

We compare analytical computations with numerical simulations for dark-matter clustering, in general relativity and in the normal branch of DGP gravity (nDGP). Our analytical frameword is the Effective Field Theory of Large-Scale Structure (EFTofLSS), which we use to compute the one-loop dark-matter power spectrum, including the resummation of infrared bulk displacement effects. We compare this to a set of 20 COLA simulations at redshifts $z = 0$, $z=0.5$, and $z =1$, and fit the free parameter of the EFTofLSS, called the speed of sound, in both $\Lambda$CDM and nDGP at each redshift. At one-loop at $z = 0$, the reach of the EFTofLSS is $k_{\rm reach}\approx 0.14 \, h { \rm Mpc^{-1}}$ for both $\Lambda$CDM and nDGP. Along the way, we compare two different infrared resummation schemes and two different treatments of the time dependence of the perturbative expansion, concluding that they agree to approximately $1\%$ over the scales of interest. Finally, we use the ratio of the COLA power spectra to make a precision measurement of the difference between the speeds of sound in $\Lambda$CDM and nDGP, and verify that this is proportional to the modification of the linear coupling constant of the Poisson equation.

The recent simultaneous detection of gravitational waves and a gamma ray burst from a neutron star merger significantly shrank the space of viable scalar-tensor theories by demanding that the speed of gravity is equal to that of light. The survived theories belong to the class of degenerate higher order scalar-tensor theories. We study whether these theories are suitable as dark energy candidates. We find scaling solutions in the matter dominated universe that lead to de Sitter solutions at late times without the cosmological constant, realising self-acceleration. We evaluate quasi-static perturbations around self-accelerating solutions and show that the stringent constraints coming from astrophysical objects and gravitational waves can be satisfied, leaving interesting possibilities to test these theories by cosmological observations.

Supernova remnants (SNRs) are the primary candidate of Galactic cosmic-ray accelerators. It is still an open issue when and how young SNRs, which typically exhibit strong synchrotron X-rays and GeV and TeV gamma-rays, undergo the state transition to middle-aged SNRs dominated by thermal X-rays and GeV gamma-rays. The SNR N132D in the Large Magellanic Cloud is an ideal target to study such a transition, exhibiting bright X-rays and gamma-rays, and with the expected age of ~2500 yrs. In this paper we present results of NuSTAR and Suzaku spectroscopy. We reveal that N132D has a nearly equilibrium plasma with a temperature of > 5 keV or a recombining plasma with a lower temperature (~1.5 keV) and a recombining timescale (net) of 8.8 (7.0--10.0)e12 cm^-3s. Together with the center filled morphology observed in the iron K line image, our results suggest that N132D is now at transition stage from a young SNR to middle-aged. We have constrained the tight upper-limit of nonthermal X-rays. Bright gamma-rays compared to faint nonthermal X-rays suggest that the gamma-rays are hadronic in origin. The spectral energy distribution from radio to gamma-rays shows a proton cut-off energy of ~30 TeV. These facts confirm that N132D is in the transition from young to middle-aged SNR. The large thermal energy of > 10^51 erg and accelerated proton energy of ~ 10^50 erg suggest the supernova explosion might have been very energetic.

We present results from the Hitomi X-ray observation of a young composite-type supernova remnant (SNR) G21.5$-$0.9, whose emission is dominated by the pulsar wind nebula (PWN) contribution. The X-ray spectra in the 0.8-80 keV range obtained with the Soft X-ray Spectrometer (SXS), Soft X-ray Imager (SXI) and Hard X-ray Imager (HXI) show a significant break in the continuum as previously found with the NuSTAR observation. After taking into account all known emissions from the SNR other than the PWN itself, we find that the Hitomi spectra can be fitted with a broken power law with photon indices of $\Gamma_1=1.74\pm0.02$ and $\Gamma_2=2.14\pm0.01$ below and above the break at $7.1\pm0.3$ keV, which is significantly lower than the NuSTAR result ($\sim9.0$ keV). The spectral break cannot be reproduced by time-dependent particle injection one-zone spectral energy distribution models, which strongly indicates that a more complex emission model is needed, as suggested by recent theoretical models. We also search for narrow emission or absorption lines with the SXS, and perform a timing analysis of PSR J1833$-$1034 with the HXI and SGD. No significant pulsation is found from the pulsar. However, unexpectedly, narrow absorption line features are detected in the SXS data at 4.2345 keV and 9.296 keV with a significance of 3.65 $\sigma$. While the origin of these features is not understood, their mere detection opens up a new field of research and was only possible with the high resolution, sensitivity and ability to measure extended sources provided by an X-ray microcalorimeter.

Phenomenological functions $\Sigma$ and $\mu$, also known as $G_{\rm light}/G$ and $G_{\rm matter}/G$, are commonly used to parameterize modifications of the growth of large-scale structure in alternative theories of gravity. We study the values these functions can take in Horndeski theories, i.e. the class of scalar-tensor theories with second order equations of motion. We restrict our attention to models that are in a broad agreement with tests of gravity and the observed cosmic expansion history. In particular, we require the speed of gravity to be equal to the speed of light today, as required by the recent detection of gravitational waves and electromagnetic emission from a binary neutron star merger. We examine the correlations between the values of $\Sigma$ and $\mu$ analytically within the quasi-static approximation, and numerically, by sampling the space of allowed solutions. We confirm that the conjecture made in [Pogosian:2016pwr], that $(\Sigma-1)(\mu -1) \ge 0$ in viable Horndeski theories, holds very well. Along with that, we check the validity of the quasi-static approximation within different corners of Horndeski theory. Our results show that, even with the tight bound on the present day speed of gravitational waves, there is room within Horndeski theories for non-trivial signatures of modified gravity at the level of linear perturbations.

Supernova remnants (SNRs) have been prime candidates for Galactic cosmic-ray accelerators. When low-energy cosmic-ray protons (LECRp) collide with interstellar gas, they ionize neutral iron atoms and emit the neutral iron line (Fe I K$\alpha$) at 6.40keV. We search for the iron K-shell line in seven SNRs from the Suzaku archive data of the Galactic plane in the $6^{\circ} \lesssim l \lesssim 40^{\circ}, |b| < 1^{\circ}$ region. All these SNRs interact with molecular clouds. We discover Fe I K$\alpha$ line emissions from five SNRs (W28, Kes 67, Kes 69, Kes 78, and W44). The spectra and morphologies suggest that the Fe I K$\alpha$ line is produced by interactions between LECRp and the adjacent cold gas. The proton energy density is estimated to be $\gtrsim$ 10-100 eV cm$^{-3}$, which is more than 10 times higher than that in the ambient interstellar medium.

In this paper, the Suzaku X-ray data of the Galactic Supernova Remnant (SNR) candidate G323.7$-$1.0 are analyzed to search for X-ray emission. Spatially-extended enhancements in the 6.4 keV line and in soft X-rays are found inside or on the radio shell. The soft X-ray enhancement would be the hottest part of the shell-like X-ray emission along the radio shell. The 6.4 keV line enhancement is detected at a significance level of $4.1 \sigma$. The lower limit of the equivalent width (EW) is 1.2 keV. The energy centroid of the 6.4 keV line is $6.40 \pm 0.04$ keV, indicating that the iron is less ionized than the Ne-like state. If the 6.4 keV line originates from ionizing non-equilibrium thermal plasma, presence of iron-rich ejecta in a low-ionization state is required, which is disfavored by the relatively old age of the SNR. The 6.4 keV line enhancement would be due to K-shell ionization of iron atoms in a dense interstellar medium by high-energy particles. Since there is no irradiating X-ray source, the origin of the 6.4 keV line enhancement is not likely the photoionization. The large EW can only be explained by K-shell ionization due to cosmic-ray protons with an energy of $\sim 10$ MeV, which might be generated by the shock acceleration in G323.7$-$1.0.

The Vainshtein mechanism, present in many models of gravity, is very effective at screening dark matter halos such that the fifth force is negligible and general relativity is recovered within their Vainshtein radii. Vainshtein screening is independent of halo mass and environment, in contrast to e.g. chameleon screening, making it difficult to test. However, our previous studies have found that the dark matter particles in filaments, walls, and voids are not screened by the Vainshtein mechanism. We therefore investigate whether cosmic voids, identified as local density minima using a watershed technique, can be used to test models of gravity that exhibit Vainshtein screening. We measure density, velocity, and screening profiles of stacked voids in cosmological $N$-body simulations using both dark matter particles and dark matter halos as tracers of the density field. We find that the voids are completely unscreened, and the tangential velocity and velocity dispersion profiles of stacked voids show a clear deviation from $\Lambda$CDM at all radii. Voids have the potential to provide a powerful test of gravity on cosmological scales.

The present paper investigates the temperature structure of the X-ray emitting plasma in the core of the Perseus cluster using the 1.8--20.0 keV data obtained with the Soft X-ray Spectrometer (SXS) onboard the Hitomi Observatory. A series of four observations were carried out, with a total effective exposure time of 338 ks and covering a central region $\sim7'$ in diameter. The SXS was operated with an energy resolution of $\sim$5 eV (full width at half maximum) at 5.9 keV. Not only fine structures of K-shell lines in He-like ions but also transitions from higher principal quantum numbers are clearly resolved from Si through Fe. This enables us to perform temperature diagnostics using the line ratios of Si, S, Ar, Ca, and Fe, and to provide the first direct measurement of the excitation temperature and ionization temperature in the Perseus cluster. The observed spectrum is roughly reproduced by a single temperature thermal plasma model in collisional ionization equilibrium, but detailed line ratio diagnostics reveal slight deviations from this approximation. In particular, the data exhibit an apparent trend of increasing ionization temperature with increasing atomic mass, as well as small differences between the ionization and excitation temperatures for Fe, the only element for which both temperatures can be measured. The best-fit two-temperature models suggest a combination of 3 and 5 keV gas, which is consistent with the idea that the observed small deviations from a single temperature approximation are due to the effects of projection of the known radial temperature gradient in the cluster core along the line of sight. Comparison with the Chandra/ACIS and the XMM-Newton/RGS results on the other hand suggests that additional lower-temperature components are present in the ICM but not detectable by Hitomi SXS given its 1.8--20 keV energy band.

The almost simultaneous detection of gravitational waves and a short gamma-ray burst from a neutron star merger has put a tight constraint on the difference between the speed of gravity and light. In the four-dimensional scalar-tensor theory with second order equations of motion, the Horndeski theory, this translates into a significant reduction of the viable parameter space of the theory. Recently, extensions of Horndeski theory, which are free from Ostrogradsky ghosts despite the presence of higher order derivatives in the equations of motion, have been identified and classified exploiting the degeneracy criterium. In these new theories, the fifth force mediated by the scalar field must be suppressed in order to evade the stringent Solar System constraints. We study the Vainshtein mechanism in the most general degenerate higher order scalar-tensor theory in which light and gravity propagate at the same speed. We find that the Vainshtein mechanism generally works outside a matter source but it is broken inside matter, similarly to beyond Horndeski theories. This leaves interesting possibilities to test these theories that are compatible with gravitational wave observations using astrophysical objects.

The Hitomi SXS spectrum of the Perseus cluster, with $\sim$5 eV resolution in the 2-9 keV band, offers an unprecedented benchmark of the atomic modeling and database for hot collisional plasmas. It reveals both successes and challenges of the current atomic codes. The latest versions of AtomDB/APEC (3.0.8), SPEX (3.03.00), and CHIANTI (8.0) all provide reasonable fits to the broad-band spectrum, and are in close agreement on best-fit temperature, emission measure, and abundances of a few elements such as Ni. For the Fe abundance, the APEC and SPEX measurements differ by 16%, which is 17 times higher than the statistical uncertainty. This is mostly attributed to the differences in adopted collisional excitation and dielectronic recombination rates of the strongest emission lines. We further investigate and compare the sensitivity of the derived physical parameters to the astrophysical source modeling and instrumental effects. The Hitomi results show that an accurate atomic code is as important as the astrophysical modeling and instrumental calibration aspects. Substantial updates of atomic databases and targeted laboratory measurements are needed to get the current codes ready for the data from the next Hitomi-level mission.

Ghost-free theories beyond the Horndeski class exhibit a partial breaking of the Vainshtein mechanism inside non-relativistic sources of finite extent. We exploit this breaking to identify new and novel astrophysical probes of these theories. Non-relativistic objects feel a gravitational force that is weaker than that predicted by general relativity. The new equation of hydrostatic equilibrium equation is derived and solved to predict the modified behaviour of stars. It is found that main-sequence stars are dimmer and cooler than their general relativity counterparts but the red giant phase is largely indistinguishable. The rotation curves and lensing potential of Milky Way-like galaxies are calculated. The circular velocities are smaller than predicted by general relativity at fixed radius and the lensing mass is larger than the dynamical mass. We discuss potential astrophysical probes of these theories and identify strong lensing as a particularly promising candidate.

We present Hitomi observations of N132D, a young, X-ray bright, O-rich core-collapse supernova remnant in the Large Magellanic Cloud (LMC). Despite a very short observation of only 3.7 ks, the Soft X-ray Spectrometer (SXS) easily detects the line complexes of highly ionized S K and Fe K with 16-17 counts in each. The Fe feature is measured for the first time at high spectral resolution. Based on the plausible assumption that the Fe K emission is dominated by He-like ions, we find that the material responsible for this Fe emission is highly redshifted at ~800 km/s compared to the local LMC interstellar medium (ISM), with a 90% credible interval of 50-1500 km/s if a weakly informative prior is placed on possible line broadening. This indicates (1) that the Fe emission arises from the supernova ejecta, and (2) that these ejecta are highly asymmetric, since no blue-shifted component is found. The S K velocity is consistent with the local LMC ISM, and is likely from swept-up ISM material. These results are consistent with spatial mapping that shows the He-like Fe concentrated in the interior of the remnant and the S tracing the outer shell. The results also show that even with a very small number of counts, direct velocity measurements from Doppler-shifted lines detected in extended objects like supernova remnants are now possible. Thanks to the very low SXS background of ~1 event per spectral resolution element per 100 ks, such results are obtainable during short pointed or slew observations with similar instruments. This highlights the power of high-spectral-resolution imaging observations, and demonstrates the new window that has been opened with Hitomi and will be greatly widened with future missions such as the X-ray Astronomy Recovery Mission (XARM) and Athena.

We show how standard Newtonian N-body simulations can be interpreted in terms of the weak-field limit of general relativity by employing the recently developed Newtonian motion gauge. Our framework allows the inclusion of radiation perturbations and the non-linear evolution of matter. We show how to construct the weak-field metric by combining Newtonian simulations with results from Einstein-Boltzmann codes. We discuss observational effects on weak lensing and ray tracing, identifying important relativistic corrections.

We report a Hitomi observation of IGR J16318-4848, a high-mass X-ray binary system with an extremely strong absorption of N_H~10^{24} cm^{-2}. Previous X-ray studies revealed that its spectrum is dominated by strong fluorescence lines of Fe as well as continuum emission. For physical and geometrical insight into the nature of the reprocessing material, we utilize the high spectroscopic resolving power of the X-ray microcalorimeter (the soft X-ray spectrometer; SXS) and the wide-band sensitivity by the soft and hard X-ray imager (SXI and HXI) aboard Hitomi. Even though photon counts are limited due to unintended off-axis pointing, the SXS spectrum resolves Fe K{\alpha_1} and K{\alpha_2} lines and puts strong constraints on the line centroid and width. The line width corresponds to the velocity of 160^{+300}_{-70} km s^{-1}. This represents the most accurate, and smallest, width measurement of this line made so far from any X-ray binary, much less than the Doppler broadening and shift expected from speeds which are characteristic of similar systems. Combined with the K-shell edge energy measured by the SXI and HXI spectra, the ionization state of Fe is estimated to be in the range of Fe I--IV. Considering the estimated ionization parameter and the distance between the X-ray source and the absorber, the density and thickness of the materials are estimated. The extraordinarily strong absorption and the absence of a Compton shoulder component is confirmed. These characteristics suggest reprocessing materials which are distributed in a narrow solid angle or scattering primarily with warm free electrons or neutral hydrogen.

The origin of the narrow Fe-K{\alpha} fluorescence line at 6.4 keV from active galactic nuclei has long been under debate; some of the possible sites are the outer accretion disk, the broad line region, a molecular torus, or interstellar/intracluster media. In February-March 2016, we performed the first X-ray microcalorimeter spectroscopy with the Soft X-ray Spectrometer (SXS) onboard the Hitomi satellite of the Fanaroff-Riley type I radio galaxy NGC 1275 at the center of the Perseus cluster of galaxies. With the high energy resolution of ~5 eV at 6 keV achieved by Hitomi/SXS, we detected the Fe-K{\alpha} line with ~5.4 {\sigma} significance. The velocity width is constrained to be 500-1600 km s$^{-1}$ (FWHM for Gaussian models) at 90% confidence. The SXS also constrains the continuum level from the NGC 1275 nucleus up to ~20 keV, giving an equivalent width ~20 eV of the 6.4 keV line. Because the velocity width is narrower than that of broad H{\alpha} line of ~2750 km s$^{-1}$, we can exclude a large contribution to the line flux from the accretion disk and the broad line region. Furthermore, we performed pixel map analyses on the Hitomi/SXS data and image analyses on the Chandra archival data, and revealed that the Fe-K{\alpha} line comes from a region within ~1.6 kpc from the NGC 1275 core, where an active galactic nucleus emission dominates, rather than that from intracluster media. Therefore, we suggest that the source of the Fe-K{\alpha} line from NGC 1275 is likely a low-covering fraction molecular torus or a rotating molecular disk which probably extends from a pc to hundreds pc scale in the active galactic nucleus system.

Extending the earlier measurements reported in Hitomi collaboration (2016, Nature, 535, 117), we examine the atmospheric gas motions within the central 100~kpc of the Perseus cluster using observations obtained with the Hitomi satellite. After correcting for the point spread function of the telescope and using optically thin emission lines, we find that the line-of-sight velocity dispersion of the hot gas is remarkably low and mostly uniform. The velocity dispersion reaches maxima of approximately 200~km~s$^{-1}$ toward the central active galactic nucleus (AGN) and toward the AGN inflated north-western `ghost' bubble. Elsewhere within the observed region, the velocity dispersion appears constant around 100~km~s$^{-1}$. We also detect a velocity gradient with a 100~km~s$^{-1}$ amplitude across the cluster core, consistent with large-scale sloshing of the core gas. If the observed gas motions are isotropic, the kinetic pressure support is less than 10\% of the thermal pressure support in the cluster core. The well-resolved optically thin emission lines have Gaussian shapes, indicating that the turbulent driving scale is likely below 100~kpc, which is consistent with the size of the AGN jet inflated bubbles. We also report the first measurement of the ion temperature in the intracluster medium, which we find to be consistent with the electron temperature. In addition, we present a new measurement of the redshift to the brightest cluster galaxy NGC~1275.

The effect of massive neutrinos on the growth of cold dark matter perturbations acts as a scale-dependent Newton's constant and leads to scale-dependent growth factors just as we often find in models of gravity beyond General Relativity. We show how to compute growth factors for $\Lambda$CDM and general modified gravity cosmologies combined with massive neutrinos in Lagrangian perturbation theory for use in COLA and extensions thereof. We implement this together with the grid-based massive neutrino method of Brandbyge and Hannestad in $\texttt{MG-PICOLA}$ and compare COLA simulations to full $\it N$-body simulations of $\Lambda$CDM and $f(R)$ gravity with massive neutrinos. Our implementation is computationally cheap if the underlying cosmology already has scale-dependent growth factors and it is shown to be able to produce results that match $\it N$-body to percent level accuracy for both the total and CDM matter power-spectra up to $k\lesssim 1 h/$Mpc.

Thanks to its high spectral resolution (~5 eV at 6 keV), the Soft X-ray Spectrometer (SXS) on board Hitomi enables us to measure the detailed structure of spatially resolved emission lines from highly ionized ions in galaxy clusters for the first time. In this series of papers, using the SXS we have measured the velocities of gas motions, metallicities and the multi-temperature structure of the gas in the core of the Perseus cluster. Here, we show that when inferring physical properties from line emissivities in systems like Perseus, the resonant scattering (RS) effect should be taken into account. In the Hitomi waveband, RS mostly affects the FeXXV He$\alpha$ line ($w$) - the strongest line in the spectrum. The flux measured by Hitomi in this line is suppressed by a factor ~1.3 in the inner ~30 kpc, compared to predictions for an optically thin plasma; the suppression decreases with the distance from the center. The $w$ line also appears slightly broader than other lines from the same ion. The observed distortions of the $w$ line flux, shape and distance dependence are all consistent with the expected effect of the resonant scattering in the Perseus core. By measuring the ratio of fluxes in optically thick ($w$) and thin (FeXXV forbidden, He$\beta$, Ly$\alpha$) lines, and comparing these ratios with predictions from Monte Carlo radiative transfer simulations, the velocities of gas motions have been obtained. The results are consistent with the direct measurements of gas velocities from line broadening described elsewhere in this series, although the systematic and statistical uncertainties remain significant. Further improvements in the predictions of line emissivities in plasma models, and deeper observations with future X-ray missions will enable RS measurements to provide powerful constraints on the amplitude and anisotropy of clusters gas motions.

Bayesian inference techniques are used to investigate situations where an additional light scalar field is present during inflation and reheating. This includes (but is not limited to) curvaton-type models. We design a numerical pipeline where $\simeq 200$ inflaton setups $\times\, 10$ reheating scenarios $= 2000$ models are implemented and we present the results for a few prototypical potentials. We find that single-field models are remarkably robust under the introduction of light scalar degrees of freedom. Models that are ruled out at the single-field level are not improved in general, because good values of the spectral index and the tensor-to-scalar ratio can only be obtained for very fine-tuned values of the extra field parameters and/or when large non-Gaussianities are produced. The only exception is quartic large-field inflation, so that the best models after Planck are of two kinds: plateau potentials, regardless of whether an extra field is added or not, and quartic large-field inflation with an extra light scalar field, in some specific reheating scenarios. Using Bayesian complexity, we also find that more parameters are constrained for the models we study than for their single-field versions. This is because the added parameters not only contribute to the reheating kinematics but also to the cosmological perturbations themselves, to which the added field contributes. The interplay between these two effects lead to a suppression of degeneracies that is responsible for having more constrained parameters.

The Galactic Diffuse X-ray Emission (GDXE) in the Milky Way Galaxy is spatially and spectrally decomposed into the Galactic Center X-ray Emission (GCXE), the Galactic Ridge X-ray Emission (GRXE) and the Galactic Bulge X-ray Emission (GBXE). The X-ray spectra of the GDXE are characterized by the strong K-shell lines of the highly ionized atoms, the brightest are the K-shell transition (principal quantum number transition of $n=2\rightarrow1$) of neutral iron (SXV-He$\alpha$), He-like iron (FeXXV-He$\alpha$) and He-like sulfur (SXV-He$\alpha$) lines. Accordingly, the GDXE is composed of a high-temperature plasma of $\sim$7 keV (HTP) and a low-temperature plasma of $\sim$1 keV (LTP), which emit the FeXXV-He$\alpha$ and SXV-He$\alpha$ lines, respectively. The FeI-K$\alpha$ line is emitted from nearly neutral irons, and hence the third component of the GDXE is a Cool Gas (CG). The FeI-K$\alpha$ distribution in the GCXE region is clumpy (FeI-K$\alpha$ clump), associated with giant molecular cloud (MC) complexes, Sagittarius A, B, C, D and E, in the central molecular zone. The origin of the FeI-K$\alpha$ clumps is the fluorescence and Thomson scattering from the MCs irradiated by past big flares of the super massive black hole Sagittarius A$^*$. The scale heights and equivalent widths of the FeI-K$\alpha$, FeXXV-He$\alpha$, and FeXXVI-Ly$\alpha$ ($n=2\rightarrow1$ transition of H-like iron) lines are different among the GCXE, GBXE and GRXE. Therefore the structure and origin are separately examined. This paper overviews the research history and the present understandings of the GDXE, in particular focus on the origin of the HTP and CG in the GCXE.

To search for giant X-ray pulses correlated with the giant radio pulses (GRPs) from the Crab pulsar, we performed a simultaneous observation of the Crab pulsar with the X-ray satellite Hitomi in the 2 -- 300 keV band and the Kashima NICT radio observatory in the 1.4 -- 1.7 GHz band with a net exposure of about 2 ks on 25 March 2016, just before the loss of the Hitomi mission.The timing performance of the Hitomi instruments was confirmed to meet the timing requirement and about 1,000 and 100 GRPs were simultaneously observed at the main and inter-pulse phases, respectively, and we found no apparent correlation between the giant radio pulses and the X-ray emission in either the main or inter-pulse phases.All variations are within the 2 sigma fluctuations of the X-ray fluxes at the pulse peaks, and the 3 sigma upper limits of variations of main- or inter- pulse GRPs are 22\% or 80\% of the peak flux in a 0.20 phase width, respectively, in the 2 -- 300 keV band.The values become 25\% or 110\% for main or inter-pulse GRPs, respectively, when the phase width is restricted into the 0.03 phase.Among the upper limits from the Hitomi satellite, those in the 4.5-10 keV and the 70-300 keV are obtained for the first time, and those in other bands are consistent with previous reports.Numerically, the upper limits of main- and inter-pulse GRPs in the 0.20 phase width are about (2.4 and 9.3) $\times 10^{-11}$ erg cm$^{-2}$, respectively. No significant variability in pulse profiles implies that the GRPs originated from a local place within the magnetosphere and the number of photon-emitting particles temporally increases.However, the results do not statistically rule out variations correlated with the GRPs, because the possible X-ray enhancement may appear due to a $>0.02$\% brightening of the pulse-peak flux under such conditions.

We extend our previous redshift space power spectrum code to the redshift space correlation function. Here we focus on the Gaussian Streaming Model (GSM). Again, the code accommodates a wide range of modified gravity and dark energy models. For the non-linear real space correlation function used in the GSM, we use the Fourier transform of the RegPT 1-loop matter power spectrum. We compare predictions of the GSM for a Vainshtein screened and Chameleon screened model as well as GR. These predictions are compared to the Fourier transform of the Taruya, Nishimichi and Saito (TNS) redshift space power spectrum model which is fit to N-Body data. We find very good agreement between the Fourier transform of the TNS model and the GSM predictions, with $\leq 6\%$ deviations in the first two correlation function multipoles for all models for separations in $50$Mpc$/h \leq s \leq 180$Mpc/$h$. Excellent agreement is found in the differences between the modified gravity and GR multipole predictions for both approaches to the redshift space correlation function, highlighting their matched ability in picking up deviations from GR. We elucidate the timeliness of such non-standard templates at the dawn of stage-IV surveys and discuss necessary preparations and extensions needed for upcoming high quality data.

We compute the power spectrum at one-loop order in standard perturbation theory for the matter density field to which a standard Lagrangian Baryonic acoustic oscillation (BAO) reconstruction technique is applied. The BAO reconstruction method corrects the bulk motion associated with the gravitational evolution using the inverse Zel'dovich approximation (ZA) for the smoothed density field. We find that the overall amplitude of one-loop contributions in the matter power spectrum substantially decrease after reconstruction. The reconstructed power spectrum thereby approaches the initial linear spectrum when the smoothed density field is close enough to linear, i.e., the smoothing scale $R_s$ larger than around 10$h^{-1}$Mpc. On smaller $R_s$,however, the deviation from the linear spectrum becomes significant on large scales ($k\lt R_s^{-1}$) due to the nonlinearity in the smoothed density field, and the reconstruction is inaccurate. Compared with N-body simulations, we show that the reconstructed power spectrum at one loop order agrees with simulations better than the unreconstructed power spectrum. We also calculate the tree-level bispectrum in standard perturbation theory to investigate non-Gaussianity in the reconstructed matter density field. We show that the amplitude of the bispectrum significantly decreases for small $k$ after reconstruction and that the tree-level bispectrum agrees well with N-body results in the weakly nonlinear regime.

We present a general parallelized and easy-to-use code to perform numerical simulations of structure formation using the COLA (COmoving Lagrangian Acceleration) method for cosmological models that exhibit scale-dependent growth at the level of first and second order Lagrangian perturbation theory. For modified gravity theories we also include screening using a fast approximate method that covers all the main examples of screening mechanisms in the literature. We test the code by comparing it to full simulations of two popular modified gravity models, namely $f(R)$ gravity and nDGP, and find good agreement in the modified gravity boost-factors relative to $\Lambda$CDM even when using a fairly small number of COLA time steps.

A flat Friedman-Roberson-Walker universe dominated by a cosmological constant ($\Lambda$) and cold dark matter (CDM) has been the working model preferred by cosmologists since the discovery of cosmic acceleration. However, tensions of various degrees of significance are known to be present among existing datasets within the $\Lambda$CDM framework. In particular, the Lyman-$\alpha$ forest measurement of the Baryon Acoustic Oscillations (BAO) by the Baryon Oscillation Spectroscopic Survey (BOSS) prefers a smaller value of the matter density fraction $\Omega_{\rm M}$ compared to the value preferred by cosmic microwave background (CMB). Also, the recently measured value of the Hubble constant, $H_0=73.24\pm1.74 \ {\rm km}\ {\rm s}^{-1} \ {\rm Mpc}^{-1}$, is $3.4\sigma$ higher than $66.93\pm0.62 \ {\rm km}\ {\rm s}^{-1} \ {\rm Mpc}^{-1}$ inferred from the Planck CMB data. In this work, we investigate if these tensions can be interpreted as evidence for a non-constant dynamical dark energy (DE). Using the Kullback-Leibler (KL) divergence to quantify the tension between datasets, we find that the tensions are relieved by an evolving DE, with the dynamical DE model preferred at a $3.5\sigma$ significance level based on the improvement in the fit alone. While, at present, the Bayesian evidence for the dynamical DE is insufficient to favour it over $\Lambda$CDM, we show that, if the current best fit DE happened to be the true model, it would be decisively detected by the upcoming DESI survey.

We present a novel cosmological solution in the framework of extended quasidilaton theory which underwent scrutiny recently. We only consider terms that do not generate the Boulware-Deser degree of freedom, hence the "ghost-free" quasidilaton theory, and show three new branches of cosmological evolution therein. One of the solutions passes the perturbative stability tests. This new solution exhibits a late time self-acceleration and all graviton polarizations acquire masses that converge to a constant in the asymptotic future. Moreover, all modes propagate at the speed of light. We propose that this solution can be used as a benchmark model for future phenomenological studies.

The Crab nebula originated from a core-collapse supernova (SN) explosion observed in 1054 A.D. When viewed as a supernova remnant (SNR), it has an anomalously low observed ejecta mass and kinetic energy for an Fe-core collapse SN. Intensive searches were made for a massive shell that solves this discrepancy, but none has been detected. An alternative idea is that the SN1054 is an electron-capture (EC) explosion with a lower explosion energy by an order of magnitude than Fe-core collapse SNe. In the X-rays, imaging searches were performed for the plasma emission from the shell in the Crab outskirts to set a stringent upper limit to the X-ray emitting mass. However, the extreme brightness of the source hampers access to its vicinity. We thus employed spectroscopic technique using the X-ray micro-calorimeter onboard the Hitomi satellite. By exploiting its superb energy resolution, we set an upper limit for emission or absorption features from yet undetected thermal plasma in the 2-12 keV range. We also re-evaluated the existing Chandra and XMM-Newton data. By assembling these results, a new upper limit was obtained for the X-ray plasma mass of <~ 1Mo for a wide range of assumed shell radius, size, and plasma temperature both in and out of the collisional equilibrium. To compare with the observation, we further performed hydrodynamic simulations of the Crab SNR for two SN models (Fe-core versus EC) under two SN environments (uniform ISM versus progenitor wind). We found that the observed mass limit can be compatible with both SN models if the SN environment has a low density of <~ 0.03 cm-3 (Fe core) or <~ 0.1 cm-3 (EC) for the uniform density, or a progenitor wind density somewhat less than that provided by a mass loss rate of 10-5 Mo yr-1 at 20 km s-1 for the wind environment.

Initial conditions for (Newtonian) cosmological N-body simulations are usually set by re-scaling the present-day power spectrum obtained from linear (relativistic) Boltzmann codes to the desired initial redshift of the simulation. This back-scaling method can account for the effect of inhomogeneous residual thermal radiation at early times, which is absent in the Newtonian simulations. We analyse this procedure from a fully relativistic perspective, employing the recently-proposed Newtonian motion gauge framework. We find that N-body simulations for LambdaCDM cosmology starting from back-scaled initial conditions can be self-consistently embedded in a relativistic space-time with first-order metric potentials calculated using a linear Boltzmann code. This space-time coincides with a simple "N-body gauge" for z<50 for all observable modes. Care must be taken, however, when simulating non-standard cosmologies. As an example, we analyse the back-scaling method in a cosmology with decaying dark matter, and show that metric perturbations become large at early times in the back-scaling approach, indicating a breakdown of the perturbative description. We suggest a suitable "forwards approach" for such cases.

We elucidate the importance of the consistent treatment of gravity-model specific non-linearities when estimating the growth of cosmological structures from redshift space distortions (RSD). Within the context of standard perturbation theory (SPT), we compare the predictions of two theoretical templates with redshift space data from COLA (COmoving Lagrangian Acceleration) simulations in the normal branch of DGP gravity (nDGP) and General Relativity (GR). Using COLA for these comparisons is validated using a suite of full N-body simulations for the same theories. The two theoretical templates correspond to the standard general relativistic perturbation equations and those same equations modelled within nDGP. Gravitational clustering non-linear effects are accounted for by modelling the power spectrum up to one loop order and redshift space clustering anisotropy is modelled using the Taruya, Nishimichi and Saito (TNS) RSD model. Using this approach, we attempt to recover the simulation's fiducial logarithmic growth parameter $f$. By assigning the simulation data with errors representing an idealised survey with a volume of $10\mbox{Gpc}^3/h^3$, we find the GR template is unable to recover fiducial $f$ to within 1$\sigma$ at $z=1$ when we match the data up to $k_{\rm max}=0.195h$/Mpc. On the other hand, the DGP template recovers the fiducial value within $1\sigma$. Further, we conduct the same analysis for sets of mock data generated for generalised models of modified gravity using SPT, where again we analyse the GR template's ability to recover the fiducial value. We find that for models with enhanced gravitational non-linearity, the theoretical bias of the GR template becomes significant for stage IV surveys. Thus, we show that for the future large data volume galaxy surveys, the self-consistent modelling of non-GR gravity scenarios will be crucial in constraining theory parameters.

We describe and demonstrate the potential of a new and very efficient method for simulating certain classes of modified gravity theories, such as the widely studied $f(R)$ gravity models. High resolution simulations for such models are currently very slow due to the highly nonlinear partial differential equation that needs to be solved exactly to predict the modified gravitational force. This nonlinearity is partly inherent, but is also exacerbated by the specific numerical algorithm used, which employs a variable redefinition to prevent numerical instabilities. The standard Newton-Gauss-Seidel iterative method used to tackle this problem has a poor convergence rate. Our new method not only avoids this, but also allows the discretised equation to be written in a form that is analytically solvable. We show that this new method greatly improves the performance and efficiency of $f(R)$ simulations. For example, a test simulation with $512^3$ particles in a box of size $512 \, \mathrm{Mpc}/h$ is now 5 times faster than before, while a Millennium-resolution simulation for $f(R)$ gravity is estimated to be more than 20 times faster than with the old method. Our new implementation will be particularly useful for running very high resolution, large-sized simulations which, to date, are only possible for the standard model, and also makes it feasible to run large numbers of lower resolution simulations for covariance analyses. We hope that the method will bring us to a new era for precision cosmological tests of gravity.

We use a set of N-body simulations employing a modified gravity (MG) model with Vainshtein screening to study matter and halo hierarchical clustering. As test-case scenarios we consider two normal branch Dvali-Gabadadze-Porrati (nDGP) gravity models with mild and strong growth rate enhancement. We study higher-order correlation functions $\xi_n(R)$ up to $n=9$ and associated hierarchical amplitudes $S_n(R)\equiv\xi_n(R)/\sigma(R)^{2n-2}$. We find that the matter PDFs are strongly affected by the fifth-force on scales up to $50h^{-1}$Mpc, and the deviations from GR are maximised at $z=0$. For reduced cumulants $S_n$, we find that at small scales $R\leq10h^{-1}$Mpc the MG is characterised by lower values, with the deviation growing from $7\%$ in the reduced skewness up to even $40\%$ in $S_5$. To study the halo clustering we use a simple abundance matching and divide haloes into thee fixed number density samples. The halo two-point functions are weakly affected, with a relative boost of the order of a few percent appearing only at the smallest pair separations ($r\leq 5h^{-1}$Mpc). In contrast, we find a strong MG signal in $S_n(R)$'s, which are enhanced compared to GR. The strong model exhibits a $>3\sigma$ level signal at various scales for all halo samples and in all cumulants. In this context, we find that the reduced kurtosis to be an especially promising cosmological probe of MG. Even the mild nDGP model leaves a $3\sigma$ imprint at small scales $R\leq3h^{-1}$Mpc, while the stronger model deviates from a GR-signature at nearly all scales with a significance of $>5\sigma$. Since the signal is persistent in all halo samples and over a range of scales, we advocate that the reduced kurtosis estimated from galaxy catalogues can potentially constitute a strong MG-model discriminatory as well as GR self-consistency test.

High-resolution X-ray spectroscopy with Hitomi was expected to resolve the origin of the faint unidentified E=3.5 keV emission line reported in several low-resolution studies of various massive systems, such as galaxies and clusters, including the Perseus cluster. We have analyzed the Hitomi first-light observation of the Perseus cluster. The emission line expected for Perseus based on the XMM-Newton signal from the large cluster sample under the dark matter decay scenario is too faint to be detectable in the Hitomi data. However, the previously reported 3.5 keV flux from Perseus was anomalously high compared to the sample-based prediction. We find no unidentified line at the reported high flux level. Taking into account the XMM measurement uncertainties for this region, the inconsistency with Hitomi is at a 99% significance for a broad dark-matter line and at 99.7% for a narrow line from the gas. We do not find anomalously high fluxes of the nearby faint K line or the Ar satellite line that were proposed as explanations for the earlier 3.5 keV detections. We do find a hint of a broad excess near the energies of high-n transitions of Sxvi (E=3.44 keV rest-frame) -- a possible signature of charge exchange in the molecular nebula and another proposed explanation for the unidentified line. While its energy is consistent with XMM pn detections, it is unlikely to explain the MOS signal. A confirmation of this interesting feature has to wait for a more sensitive observation with a future calorimeter experiment.

Theories of gravity in the beyond Horndeski class recover the predictions of general relativity in the solar system whilst admitting novel cosmologies, including late-time de Sitter solutions in the absence of a cosmological constant. Deviations from Newton's law are predicted inside astrophysical bodies, which allow for falsifiable, smoking-gun tests of the theory. In this work we study the pulsations of stars by deriving and solving the wave equation governing linear adiabatic oscillations to find the modified period of pulsation. Using both semi-analytic and numerical models, we perform a preliminary survey of the stellar zoo in an attempt to identify the best candidate objects for testing the theory. Brown dwarfs and Cepheid stars are found to be particularly sensitive objects and we discuss the possibility of using both to test the theory.

Theories of gravity in the beyond Horndeski class encompass a wide range of scalar-tensor theories that will be tested on cosmological scales over the coming decade. In this work, we investigate the possibility of testing them in the strong-field regime by looking at the properties of compact objects-neutron, hyperon, and quark stars-embedded in an asymptotically de Sitter space-time, for a specific subclass of theories. We extend previous works to include slow rotation and find a relation between the dimensionless moment of intertia, ($\bar{I}=Ic^2/G_{\rm N} M^3$), and the compactness, $\cal{C}=G_{\rm N} M/Rc^2$ (an $\bar{I}$-$\cal{C}$ relation), independent of the equation of state, that is reminiscent of but distinct from the general relativity prediction. Several of our equations of state contain hyperons and free quarks, allowing us to revisit the hyperon puzzle. We find that the maximum mass of hyperon stars can be larger than $2M_\odot$ for small values of the beyond Horndeski parameter, thus providing a resolution of the hyperon puzzle based on modified gravity. Moreover, stable quark stars exist when hyperonic stars are unstable, which means that the phase transition from hyperon to quark stars is predicted just as in general relativity, albeit with larger quark star masses. Two important and potentially observable consequences of some of the theories we consider are the existence of neutron stars in a range of masses significantly higher than in GR, and $\bar{I}$-$\mathcal{C}$ relations that differ from their GR counterparts. In the former case, we find objects that, if observed, could not be accounted for in GR because they violate the usual GR causality condition. We end by discussing several difficult technical issues that remain to be addressed in order to reach more realistic predictions that may be tested using gravitational wave searches or neutron star observations.

We derive the first systematic observational constraints on reheating in models of inflation where an additional light scalar field contributes to primordial density perturbations and affects the expansion history during reheating. This encompasses the original curvaton model but also covers a larger class of scenarios. We find that, compared to the single-field case, lower values of the energy density at the end of inflation and of the reheating temperature are preferred when an additional scalar field is introduced. For instance, if inflation is driven by a quartic potential, which is one of the most favoured models when a light scalar field is added, the upper bound $T_{\mathrm{reh}}<5\times 10^{4}\,\mathrm{GeV}$ on the reheating temperature $T_{\mathrm{reh}}$ is derived, and the implications of this value on post-inflationary physics are discussed. The information gained about reheating is also quantified and it is found that it remains modest in plateau inflation (though still larger than in the single-field version of the model) but can become substantial in quartic inflation. The role played by the vev of the additional scalar field at the end of inflation is highlighted, and opens interesting possibilities for exploring stochastic inflation effects that could determine its distribution.

This paper reports detailed K-shell line profiles of iron (Fe) and nickel (Ni) of the Galactic Center X-ray Emission (GCXE), Galactic Bulge X-ray Emission (GBXE), Galactic Ridge X-ray Emission (GRXE), magnetic Cataclysmic Variables (mCVs), non-magnetic Cataclysmic Variables (non-mCVs), and coronally Active Binaries (ABs). For the study of the origin of the GCXE, GBXE, and GRXE, the spectral analysis is focused on equivalent widths of the Fe K$\alpha$, Fe He$\alpha$, and Fe Ly$\alpha$ lines. The global spectrum of the GBXE is reproduced by a combination of the mCVs, non-mCVs, and ABs spectra. On the other hand, the GRXE spectrum shows significant data excesses at the Fe K$\alpha$ and Fe He$\alpha$ line energies. This means that additional components other than mCVs, non-mCVs, and ABs are required, which have symbiotic phenomena of cold gas and very high-temperature plasma. The GCXE spectrum shows larger excesses than those found in the GRXE spectrum at all the K-shell lines of iron and nickel. Among them the largest ones are the Fe K$\alpha$, Fe He$\alpha$, Fe Ly$\alpha$, and Fe Ly$\beta$ lines. Together with the fact that the scale heights of the Fe K$\alpha$, Fe He$\alpha$, and Fe Ly$\alpha$ lines are similar to that of the central molecular zone (CMZ), the excess components would be related to high-energy activity in the extreme envelopment of the CMZ.

We present all scalar-tensor Lagrangians that are cubic in second derivatives of a scalar field, and that are degenerate, hence avoiding Ostrogradsky instabilities. Thanks to the existence of constraints, they propagate no more than three degrees of freedom, despite having higher order equations of motion. We also determine the viable combinations of previously identified quadratic degenerate Lagrangians and the newly established cubic ones. Finally, we study whether the new theories are connected to known scalar-tensor theories such as Horndeski and beyond Horndeski, through conformal or disformal transformations.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

The extended quasidilaton theory is one of the simplest Lorentz-invariant massive gravity theories which can accommodate a stable self-accelerating vacuum solution. In this paper we revisit this theory and study the effect of matter fields. For a matter sector that couples minimally to the physical metric, we find hints of a Jeans type instability in the IR. In the analogue k-essence field set-up, this instability manifests itself as an IR ghost for the scalar field perturbation, but this can be interpreted as a classical instability that becomes relevant below some momentum scale in terms of matter density perturbations. We also consider the effect of the background evolution influenced by matter on the stability of the gravity sector perturbations. In particular, we address the previous claims of ghost instability in the IR around the late time attractor. We show that, although the matter-induced modification of the evolution potentially brings tension to the stability conditions, one goes beyond the regime of validity of the effective theory well before the solutions become unstable. We also draw attention to the fact that the IR stability conditions are also enforced by the existence requirements of consistent background solutions.

This work studies relativistic stars in beyond Horndeski scalar-tensor theories that exhibit a breaking of the Vainshtein mechanism inside matter, focusing on a model based on the quartic beyond Horndeski Lagrangian. We self-consistently derive the scalar field profile for static spherically symmetric objects in asymptotically de Sitter space-time and show that the Vainshtein breaking branch of the solutions is the physical branch thereby resolving several ambiguities with non-relativistic frameworks. The geometry outside the star is shown to be exactly Schwarzschild-de Sitter and therefore the PPN parameter $\beta_{\rm PPN}=1$, confirming that the external screening works at the post-Newtonian level. The Tolman-Oppenheimer-Volkoff (TOV) equations are derived and a new lower bound on the Vainshtein breaking parameter $\Upsilon_1>-4/9$ is found by requiring the existence of static spherically symmetric stars. Focusing on the unconstrained case where $\Upsilon_1<0$, we numerically solve the TOV equations for polytropic and realistic equations of state and find stars with larger radii at fixed mass. Furthermore, the maximum mass can increase dramatically and stars with masses in excess of $3M_\odot$ can be found for relatively small values of the Vainshtein breaking parameter. We re-examine white dwarf stars and show that post-Newtonian corrections are important in beyond Horndeski theories and therefore the bounds coming from previous analyses should be revisited.

The Beyond Horndeski class of alternative gravity theories allow for Self-accelerating de-Sitter cosmologies with no need for a cosmological constant. This makes them viable alternatives to $\Lambda$CDM and so testing their small-scale predictions against General Relativity is of paramount importance. These theories generically predict deviations in both the Newtonian force law and the gravitational lensing of light inside extended objects. Therefore, by simultaneously fitting the X-ray and lensing profiles of galaxy clusters new constraints can be obtained. In this work, we apply this methodology to the stacked profiles of 58 high-redshift ($ 0.1<z<1.2$) clusters using X-ray surface brightness profiles from the XMM Cluster Survey and weak lensing profiles from CFHTLenS. By performing a multi-parameter Markov chain Monte Carlo analysis, we are able to place new constraints on the parameters governing deviations from Newton's law $\Upsilon_{1}=-0.11^{+0.93}_{-0.67}$ and light bending $\Upsilon_{2}=-0.22^{+1.22}_{-1.19}$. Both constraints are consistent with General Relativity, for which $\Upsilon_{1}=\Upsilon_{2}=0$. We present here the first observational constraints on $\Upsilon_{2}$, as well as the first extragalactic measurement of both parameters.

The standard numerical tools for studying non-linear collapse of matter are Newtonian $N$-body simulations. Previous work has shown that these simulations are in accordance with General Relativity (GR) up to first order in perturbation theory, provided that the effects from radiation can be neglected. In this paper we show that the present day matter density receives more than 1$\%$ corrections from radiation on large scales if Newtonian simulations are initialised before $z=50$. We provide a relativistic framework in which \emph{unmodified} Newtonian simulations are compatible with linear GR even in the presence of radiation. Our idea is to use GR perturbation theory to keep track of the evolution of relativistic species and the relativistic space-time consistent with the Newtonian trajectories computed in $N$-body simulations. If metric potentials are sufficiently small, they can be computed using a first-order Einstein--Boltzmann code such as CLASS. We make this idea rigorous by defining a class of GR gauges, the \emph{Newtonian motion} gauges, which are defined such that matter particles follow Newtonian trajectories. We construct a simple example of a relativistic space-time within which unmodified Newtonian simulations can be interpreted.

We use two new hydrodynamical simulations of $\Lambda$CDM and $f(R)$ gravity to test the methodology used by Wilcox et al. 2015 (W15) in constraining the effects of a fifth force on the profiles of clusters of galaxies. We construct realistic simulated stacked weak lensing and X-ray surface brightness cluster profiles from these cosmological simulations, and then use these data projected along various lines-of-sight to test the spherical symmetry of our stacking procedure. We also test the applicability of the NFW profile to model weak lensing profiles of clusters in $f(R)$ gravity. Finally, we test the validity of the analytical model developed in W15 against the simulated profiles. Overall, we find our methodology is robust and broadly agrees with these simulated data. We also apply our full Markov Chain Monte Carlo (MCMC) analysis from W15 to our simulated X-ray and lensing profiles, providing consistent constraints on the modified gravity parameters as obtained from the real cluster data, e.g. for our $\Lambda$CDM simulation we obtain $|f_{\rm{R}0}| < 8.3 \times 10^{-5}$ (95% CL), which is in good agreement with the W15 measurement of $|f_{\rm{R}0}| < 6 \times 10^{-5}$. Overall, these tests confirm the power of our methodology which can now be applied to larger cluster samples available with the next generation surveys.

We present a fully relativistic calculation of the matter bispectrum at second order in cosmological perturbation theory assuming a Gaussian primordial curvature perturbation. For the first time we perform a full numerical integration of the bispectrum for both baryons and cold dark matter using the second-order Einstein-Boltzmann code, SONG. We review previous analytical results and provide an improved analytic approximation for the second-order kernel in Poisson gauge which incorporates Newtonian nonlinear evolution, relativistic initial conditions, the effect of radiation at early times and the cosmological constant at late times. Our improved kernel provides a percent level fit to the full numerical result at late times for most configurations, including both equilateral shapes and the squeezed limit. We show that baryon acoustic oscillations leave an imprint in the matter bispectrum, making a significant impact on squeezed shapes.

We develop a code to produce the power spectrum in redshift space based on standard perturbation theory (SPT) at 1-loop order. The code can be applied to a wide range of modified gravity and dark energy models using a recently proposed numerical method by A.Taruya. This includes Horndeski's theory with a general potential, which accommodates both chameleon and Vainshtein screening mechanisms and provides a non-linear extension of the effective theory of dark energy up to third order. Focus is on a recent non-linear model of the redshift space power spectrum which has been shown to model the anisotropy very well at relevant scales for the SPT framework, as well as capturing relevant non-linear effects typical of modified gravity theories. We provide consistency checks of the code against established results and elucidate it's application within the light of upcoming high precision RSD data.

Clusters of galaxies are the most massive gravitationally-bound objects in the Universe and are still forming. They are thus important probes of cosmological parameters and a host of astrophysical processes. Knowledge of the dynamics of the pervasive hot gas, which dominates in mass over stars in a cluster, is a crucial missing ingredient. It can enable new insights into mechanical energy injection by the central supermassive black hole and the use of hydrostatic equilibrium for the determination of cluster masses. X-rays from the core of the Perseus cluster are emitted by the 50 million K diffuse hot plasma filling its gravitational potential well. The Active Galactic Nucleus of the central galaxy NGC1275 is pumping jetted energy into the surrounding intracluster medium, creating buoyant bubbles filled with relativistic plasma. These likely induce motions in the intracluster medium and heat the inner gas preventing runaway radiative cooling; a process known as Active Galactic Nucleus Feedback. Here we report on Hitomi X-ray observations of the Perseus cluster core, which reveal a remarkably quiescent atmosphere where the gas has a line-of-sight velocity dispersion of 164+/-10 km/s in a region 30-60 kpc from the central nucleus. A gradient in the line-of-sight velocity of 150+/-70 km/s is found across the 60 kpc image of the cluster core. Turbulent pressure support in the gas is 4% or less of the thermodynamic pressure, with large scale shear at most doubling that estimate. We infer that total cluster masses determined from hydrostatic equilibrium in the central regions need little correction for turbulent pressure.

This paper reports the analysis of the X-ray spectra of the Galactic diffuse X-ray emission (GDXE) in the Suzaku archive. The fluxes of the Fe I K alpha (6.4 keV), Fe XXV,He alpha (6.7 keV) and Fe XXVI Ly alpha (6.97 keV) lines are separately determined. From the latitude distributions, we confirm that the GDXE is decomposed into the Galactic center (GCXE), the Galactic bulge (GBXE) and the Galactic ridge (GRXE) X-ray emissions. The scale heights (SHs) of the Fe XXV He alpha line of the GCXE, GBXE and GRXE are determined to be ~40, ~310 and ~140 pc, while those of the Fe I K alpha line are ~30, ~160 and ~70 pc, respectively. The mean equivalent widths (EWs) of the sum of the Fe XXV He alpha and Fe XXVI Ly alpha lines are ~750 eV, ~600 eV and ~550 eV, while those of the Fe I K alpha line are ~150~eV, ~60~eV and ~100~eV for the GCXE, GBXE and GRXE, respectively. The origin of the GBXE, GRXE and GCXE is separately discussed based on the new results of the SHs and EWs, in comparison with those of the Cataclysmic Variables (CVs), Active Binaries (ABs) and Coronal Active stars (CAs).

In General Relativity, the constraint equation relating metric and density perturbations is inherently nonlinear, leading to an effective non-Gaussianity in the dark matter density field on large scales - even if the primordial metric perturbation is Gaussian. Intrinsic non-Gaussianity in the large-scale dark matter overdensity in GR is real and physical. However, the variance smoothed on a local physical scale is not correlated with the large-scale curvature perturbation, so that there is no relativistic signature in the galaxy bias when using the simplest model of bias. It is an open question whether the observable mass proxies such as luminosity or weak lensing correspond directly to the physical mass in the simple halo bias model. If not, there may be observables that encode this relativistic signature.

We study new consistent scalar-tensor theories of gravity recently introduced by Langlois and Noui with potentially interesting cosmological applications. We derive the conditions for the existence of a primary constraint that prevents the propagation of an additional dangerous mode associated with higher order equations of motion. We then classify the most general, consistent scalar-tensor theories that are at most quadratic in the second derivatives of the scalar field. In addition, we investigate the possible connection between these theories and (beyond) Horndeski through conformal and disformal transformations. Finally, we point out that these theories can be associated with new operators in the effective field theory of dark energy, which might open up new possibilities to test dark energy models in future surveys.

We review recent progress in the construction of modified gravity models as alternatives to dark energy as well as the development of cosmological tests of gravity. Einstein's theory of General Relativity (GR) has been tested accurately within the local universe i.e. the Solar System, but this leaves the possibility open that it is not a good description of gravity at the largest scales in the Universe. This being said, the standard model of cosmology assumes GR on all scales. In 1998, astronomers made the surprising discovery that the expansion of the Universe is accelerating, not slowing down. This late-time acceleration of the Universe has become the most challenging problem in theoretical physics. Within the framework of GR, the acceleration would originate from an unknown dark energy. Alternatively, it could be that there is no dark energy and GR itself is in error on cosmological scales. In this review, we first give an overview of recent developments in modified gravity theories including $f(R)$ gravity, braneworld gravity, Horndeski theory and massive/bigravity theory. We then focus on common properties these models share, such as screening mechanisms they use to evade the stringent Solar System tests. Once armed with a theoretical knowledge of modified gravity models, we move on to discuss how we can test modifications of gravity on cosmological scales. We present tests of gravity using linear cosmological perturbations and review the latest constraints on deviations from the standard $\Lambda$CDM model. Since screening mechanisms leave distinct signatures in the non-linear structure formation, we also review novel astrophysical tests of gravity using clusters, dwarf galaxies and stars.

Determining the most general, consistent scalar tensor theory of gravity is important for building models of inflation and dark energy. In this work we investigate the number of degrees of freedom present in the theory of beyond Horndeski. We discuss how to construct the theory from the extrinsic curvature of the constant scalar field hypersurface, and find a simple expression for the action which guarantees the existence of the primary constraint necessary to avoid the Ostrogradsky instability. Our analysis is completely gauge-invariant. However we confirm that, mixing together beyond Horndeski with a different order of Horndeski, obstructs the construction of this primary constraint. Instead, when the mixing is between actions of the same order, the theory can be mapped to Horndeski through a generalised disformal transformation. This mapping however is impossible with beyond Horndeski alone, since we find that the theory is invariant under such a transformation. The picture that emerges is that beyond Horndeski is a healthy but isolated theory: combined with Horndeski, it either becomes Horndeski, or likely propagates a ghost.

We investigate the nonlocal gravity theory by deriving nonlocal equations of motion using the traditional variation principle in a homogeneous background. We focus on a class of models with a linear nonlocal modification term in the action. It is found that the resulting equations of motion contain the advanced Green's function, implying that there is an acausality problem. As a consequence, a divergence arises in the solutions due to contributions from the future infinity unless the Universe will go back to the radiation dominated era or become the Minkowski spacetime in the future. We also discuss the relation between the original nonlocal equations and its biscalar-tensor representation and identify the auxiliary fields with the corresponding original nonlocal terms. Finally, we show that the acusality problem cannot be avoided by any function of nonlocal terms in the action.

Abell 754 is a galaxy cluster in which an ongoing merger is evident on the plane of the sky, from the southeast to the northwest. We study the spatial variation of the X-ray spectra observed with Suzaku along the merging direction, centering on the Fe Ly-alpha / Fe He-alpha line ratio to search for possible deviation from ionization equilibrium. Fitting with a single temperature collisional non-equilibrium plasma model shows that the electron temperature increases from the southeast to the northwest. The ionization parameter is consistent with that in equilibrium (n_et>10^{13} s cm^{-3}) except for a specific region with the highest temperature (kT=13.3^{+1.4}_{-1.1} keV) where n_et=10^{11.6^{+0.6}_{-1.7}} s cm^{-3}. The elapsed time from the plasma heating estimated from the ionization parameter is 0.36-76 Myr at the 90% confidence level. This time scale is quite short but consistent with the traveling time of a shock to pass through that region. We thus interpret that the non-equilibrium ionization plasma in Abell 754 observed is a remnant of the shock heating in the merger process. We, however, note that the X-ray spectrum of the specific region where the non-equilibrium is found can also be fitted with a collisional ionization plasma model with two temperatures, low kT=4.2^{+4.2}_{-1.5} keV and very high kT > 19.3 keV. The very high temperature component is alternatively fitted with a power law model. Either of these spectral models is interpreted as a consequence of the ongoing merger process as in the case of that with the non-equilibrium ionization plasma.

Modified theories of gravity provide us with a unique opportunity to generate innovative tests of gravity. In Chameleon f(R) gravity, the gravitational potential differs from the weak-field limit of general relativity (GR) in a mass dependent way. We develop a probe of gravity which compares high mass clusters, where Chameleon effects are weak, to low mass clusters, where the effects can be strong. We utilize the escape velocity edges in the radius/velocity phase space to infer the gravitational potential profiles on scales of 0.3-1 virial radii. We show that the escape edges of low mass clusters are enhanced compared to GR, where the magnitude of the difference depends on the background field value |fR0|. We validate our probe using N-body simulations and simulated light cone galaxy data. For a DESI (Dark Energy Spectroscopic Instrument) Bright Galaxy Sample, including observational systematics, projection effects, and cosmic variance, our test can differentiate between GR and Chameleon f(R) gravity models, |fR0| = 4e-6 (2e-6) at > 5{\sigma} (> 2{\sigma}), more than an order of magnitude better than current cluster-scale constraints.

We present a science forecast for the eBOSS survey, part of the SDSS-IV project, which is a spectroscopic survey using multiple tracers of large-scale structure, including luminous red galaxies (LRGs), emission line galaxies (ELGs) and quasars (both as a direct probe of structure and through the Ly-$\alpha$ forest). Focusing on discrete tracers, we forecast the expected accuracy of the baryonic acoustic oscillation (BAO), the redshift-space distortion (RSD) measurements, the $f_{\rm NL}$ parameter quantifying the primordial non-Gaussianity, the dark energy and modified gravity parameters. We also use the line-of-sight clustering in the Ly-$\alpha$ forest to constrain the total neutrino mass. We find that eBOSS LRGs ($0.6<z<1.0$) (combined with the BOSS LRGs at $z>0.6$), ELGs ($0.6<z<1.2$) and Clustering Quasars (CQs) ($0.6<z<2.2$) can achieve a precision of 1%, 2.2% and 1.6% precisions, respectively, for spherically averaged BAO distance measurements. Using the same samples, the constraint on $f\sigma_8$ is expected to be 2.5%, 3.3% and 2.8% respectively. For primordial non-Gaussianity, eBOSS alone can reach an accuracy of $\sigma(f_{\rm NL})\sim10-15$, depending on the external measurement of the galaxy bias and our ability to model large-scale systematic errors. eBOSS can at most improve the dark energy Figure of Merit (FoM) by a factor of $3$ for the Chevallier-Polarski-Linder (CPL) parametrisation, and can well constrain three eigenmodes for the general equation-of-state parameter (Abridged).

The initial conditions for Newtonian $N$-body simulations are usually generated by applying the Zel'dovich approximation to the initial displacements of the particles using an initial power spectrum of density fluctuations generated by an Einstein-Boltzmann solver. We show that in most gauges the initial displacements generated in this way receive a first-order relativistic correction. We define a new gauge, the $N$-body gauge, in which this relativistic correction vanishes and show that a conventional Newtonian $N$-body simulation includes all first-order relativistic contributions (in the absence of radiation) if we identify the coordinates in Newtonian simulations with those in the relativistic $N$-body gauge.

This paper reports results of a Suzaku observation of the supernova remnant (SNR) Kes 79 (G33.6+0.1). The X-ray spectrum is best fitted by a two-temperature model: a non-equilibrium ionization (NEI) plasma and a collisional ionization equilibrium (CIE) plasma. The NEI plasma is spatially confined within the inner radio shell with kT~0.8 keV, while the CIE plasma is found in more spatially extended regions associated with the outer radio shell with kT~0.2 keV and solar abundance. Therefore, the NEI plasma is attributable to the SN ejecta and the CIE plasma is forward shocked interstellar medium. In the NEI plasma, we discovered K-shell line of Al, Ar and Ca for the first time. The abundance pattern and estimated mass of the ejecta are consistent with the core-collapse supernova explosion of a ~30-40 solar mass progenitor star. An Fe line with center energy of ~6.4 keV is also found in the southeast (SE) portion of the SNR, a close peripheral region around dense molecular clouds. One possibility is that the line is associated with the ejecta. However, the centroid energy of ~6.4 keV and the spatial distribution of enhancement near the SE peripheral do not favor this scenario. Since the ~6.4 keV emitting region coincides to the molecular clouds, we propose another possibility that the Fe line is due to K-shell ionization of neutral Fe by the interaction of locally accelerated protons (LECRp) with the surrounding molecular cloud. Both these possibilities, heated ejecta or LECRp origin, are discussed based on the observational facts.

We investigate whether the predictions of single-field models of inflation are robust under the introduction of additional scalar degrees of freedom, and whether these extra fields change the potentials for which the data show the strongest preference. We study the situation where an extra light scalar field contributes both to the total curvature perturbations and to the reheating kinematic properties. Ten reheating scenarios are identified, and all necessary formulas allowing a systematic computation of the predictions for this class of models are derived. They are implemented in the public library ASPIC, which contains more than 75 single-field potentials. This paves the way for a forthcoming full Bayesian analysis of the problem. A few representative examples are displayed and discussed.

Many modified gravity theories are under consideration in cosmology as the source of the accelerated expansion of the universe and linear perturbation theory, valid on the largest scales, has been examined in many of these models. However, smaller non-linear scales offer a richer phenomenology with which to constrain modified gravity theories. Here, we consider the Hu-Sawicki form of $f(R)$ gravity and apply the post-Friedmann approach to derive the leading order equations for non-linear scales, i.e. the equations valid in the Newtonian-like regime. We reproduce the standard equations for the scalar field, gravitational slip and the modified Poisson equation in a coherent framework. In addition, we derive the equation for the leading order correction to the Newtonian regime, the vector potential. We measure this vector potential from $f(R)$ N-body simulations at redshift zero and one, for two values of the $f_{R_0}$ parameter. We find that the vector potential at redshift zero in $f(R)$ gravity can be close to 50\% larger than in GR on small scales for $|f_{R_0}|=1.289\times10^{-5}$, although this is less for larger scales, earlier times and smaller values of the $f_{R_0}$ parameter. Similarly to in GR, the small amplitude of this vector potential suggests that the Newtonian approximation is highly accurate for $f(R)$ gravity, and also that the non-linear cosmological behaviour of $f(R)$ gravity can be completely described by just the scalar potentials and the $f(R)$ field.

The stellar distribution derived from an $H$ and $K_{\mathrm S}$-band survey of the central region of our Galaxy is compared with the Fe XXV K$\alpha$ (6.7 keV) line intensity observed with the Suzaku satellite. The survey is for the Galactic coordinates $|l| \lesssim 3^{\circ}.0$ and $|b| \lesssim 1^{\circ}.0$ (equivalent to 0.8 kpc $\times$ 0.3 kpc for $R_0 = 8$ kpc), and the number-density distribution $N(K_{\mathrm S,0}; l, b)$ of stars is derived using the extinction-corrected magnitude $K_{\mathrm S,0}=10.5$. This is deep enough to probe the old red giant population and in turn to estimate the ($l$, $b$) distribution of faint X-ray point sources such as coronally active binaries and cataclysmic variables. In the Galactic plane ($b=0^{\circ}$), $N(10.5; l, b)$ increases to the Galactic center as $|l|^{-0.30 \pm 0.03}$ in the range of $-0^{\circ}.1 \geq l \geq -0^{\circ}.7$, but this increase is significantly slower than the increase ($|l|^{-0.44 \pm 0.02}$ ) of the Fe XXV K$\alpha$ line intensity. If normalized with the ratios in the outer region $1^{\circ}.5 \leq |l| \leq 2^{\circ}.8$, where faint X-ray point sources are argued to dominate the diffuse Galactic X-ray ridge emission, the excess of the Fe XXV K$\alpha$ line intensity over the stellar number density is at least a factor of two at $|l| = 0^{\circ}.1$. This indicates that a significant part of the Galactic center diffuse emission arises from a truly diffuse optically-thin thermal plasma, and not from an unresolved collection of faint X-ray point sources related to the old stellar population.

We report re-analyses of the Suzaku observations of the Galactic supernova remnant (SNR), G272.2$-$3.2, for which the previous studies were limited below 3 keV. With careful data reduction and background subtraction, we discover the K-shell lines of Ar, Ca, and Fe above 3 keV. The X-ray spectrum of G272.2$-$3.2 consists of two components, a low-temperature collisional ionization equilibrium (CIE) plasma ($kT_{\rm e} \sim 0.2$ keV) and a high-temperature non-equilibrium ionization (NEI) plasma ($kT_{\rm e} = 0.6$-$3$ keV). The CIE plasma has solar abundances over the entire area, hence it would originate from the interstellar medium. On the other hand, the abundances of the NEI plasma increase toward the inner region, suggesting the ejecta origin. The line center energy of the Fe K-shell emission ($\sim 6.4$ keV) suggests that the ejecta are recently heated by the reverse shock, a common feature in Type Ia SNRs.

We report the first detection of thermal X-ray line emission from the supernova remnant (SNR) RX J1713.7-3946, the prototype of the small class of synchrotron dominated SNRs. A softness-ratio map generated using XMM-Newton data shows that faint interior regions are softer than bright shell regions. Using Suzaku and deep XMM-Newton observations, we have extracted X-ray spectra from the softest area, finding clear line features at 1 keV and 1.35 keV. These lines can be best explained as Ne Ly-alpha and Mg He-alpha from a thermal emission component. Since the abundance ratios of metals to Fe are much higher than solar values in the thermal component, we attribute the thermal emission to reverse-shocked SN ejecta. The measured Mg/Ne, Si/Ne, and Fe/Ne ratios of 2.0-2.6, 1.5-2.0, and <0.05 solar suggest that the progenitor star of RX J1713.7-3946 was a relatively low-mass star (<~20 M_sun), consistent with a previous inference based on the effect of stellar winds of the progenitor star on the surrounding medium. Since the mean blastwave speed of ~6000 km/s (the radius of 9.6 pc divided by the age of 1600 yr) is relatively fast compared with other core-collapse SNRs, we propose that RX J1713.7-3946 is a result of a Type Ib/c supernova whose progenitor was a member of an interacting binary. While our analysis provides strong evidence for X-ray line emission, our interpretation of its nature as thermal emission from SN ejecta requires further confirmation especially through future precision spectroscopic measurements using ASTRO-H.

We discuss the question of gauge choice when analysing relativistic density perturbations at second order. We compare Newtonian and General Relativistic approaches. Some misconceptions in the recent literature are addressed. We show that the comoving-synchronous gauge is the unique gauge in General Relativity that corresponds to the Lagrangian frame and is entirely appropriate to describe the matter overdensity at second order. The comoving-synchronous gauge is the simplest gauge in which to describe Lagrangian bias at second order.

Self-consistent ${\it N}$-body simulations of modified gravity models are a key ingredient to obtain rigorous constraints on deviations from General Relativity using large-scale structure observations. This paper provides the first detailed comparison of the results of different ${\it N}$-body codes for the $f(R)$, DGP, and Symmetron models, starting from the same initial conditions. We find that the fractional deviation of the matter power spectrum from $\Lambda$CDM agrees to better than $1\%$ up to $k \sim 5-10~h/{\rm Mpc}$ between the different codes. These codes are thus able to meet the stringent accuracy requirements of upcoming observational surveys. All codes are also in good agreement in their results for the velocity divergence power spectrum, halo abundances and halo profiles. We also test the quasi-static limit, which is employed in most modified gravity ${\it N}$-body codes, for the Symmetron model for which the most significant non-static effects among the models considered are expected. We conclude that this limit is a very good approximation for all of the observables considered here.

We develop the transport operator formalism, a new line-of-sight integration framework to calculate the anisotropies of the Cosmic Microwave Background (CMB) at the linear and non-linear level. This formalism utilises a transformation operator that removes all inhomogeneous propagation effects acting on the photon distribution function, thus achieving a split between perturbative collisional effects at recombination and non-perturbative line-of-sight effects at later times. The former can be computed in the framework of standard cosmological perturbation theory with a second-order Boltzmann code such as SONG, while the latter can be treated within a separate perturbative scheme allowing the use of non-linear Newtonian potentials. We thus provide a consistent framework to compute all physical effects contained in the Boltzmann equation and to combine the standard remapping approach with Boltzmann codes at any order in perturbation theory, without assuming that all sources are localised at recombination.

We analyse the clustering of cosmic large scale structure using a consistent modified gravity perturbation theory, accounting for anisotropic effects along and transverse to the line of sight. The growth factor has a particular scale dependence in f(R) gravity and we fit for the shape parameter f_{R0} simultaneously with the distance and the large scale (general relativity) limit of the growth function. Using more than 690,000 galaxies in the Baryon Oscillation Spectroscopy Survey Data Release 11, we find no evidence for extra scale dependence, with the 95\% confidence upper limit |f_{R0}| <8 \times 10^{-4}. Future clustering data, such as from the Dark Energy Spectroscopic Instrument, can use this consistent methodology to impose tighter constraints.

Vector theories with non-linear derivative self-interactions that break gauge symmetries have been shown to have interesting cosmological applications. In this paper we introduce a way to spontaneously break the gauge symmetry and construct these theories via a Higgs mechanism. In addition to the purely gauge field interactions, our method generates new ghost-free scalar-vector interactions between the Higgs field and the gauge boson. We show how these additional terms are found to reduce, in a suitable decoupling limit, to scalar bi-Galileon interactions between the Higgs field and Goldstone bosons. Our formalism is first developed in the context of abelian symmetry, which allows us to connect with earlier work on the extension of the Proca action. We then show how this formalism is straightforwardly generalised to generate theories with non-abelian symmetry.

Theories which modify general relativity to explain the accelerated expansion of the Universe often use screening mechanisms to satisfy constraints on Solar System scales. We investigate the effects of the cosmic web and the local environmental density of dark matter halos on the screening properties of the Vainshtein and chameleon screening mechanisms. We compare the cosmic web morphology of dark matter particles, mass functions of dark matter halos, mass and radial dependence of screening, velocity dispersions and peculiar velocities, and environmental dependence of screening mechanisms in $f(R)$ and nDGP models. Using the ORIGAMI cosmic web identification routine we find that the Vainshtein mechanism depends on the cosmic web morphology of dark matter particles, since these are defined according to the dimensionality of their collapse, while the chameleon mechanism shows no morphology dependence. The chameleon screening of halos and their velocity dispersions depend on halo mass, and small halos and subhalos can be environmentally screened in the chameleon mechanism. On the other hand, the screening of halos in the Vainshtein mechanism does not depend on mass nor environment, and their velocity dispersions are suppressed. The peculiar velocities of halos in the Vainshtein mechanism are enhanced because screened objects can still feel the fifth force generated by external fields, while peculiar velocities of chameleon halos are suppressed when the halo centers are screened.

The chameleon gravity model postulates the existence of a scalar field that couples with matter to mediate a fifth force. If it exists, this fifth force would influence the hot X-ray emitting gas filling the potential wells of galaxy clusters. However, it would not influence the clusters' weak lensing signal. Therefore, by comparing X-ray and weak lensing profiles, one can place upper limits on the strength of a fifth force. This technique has been attempted before using a single, nearby cluster (Coma, $z=0.02$). Here we apply the technique to the stacked profiles of 58 clusters at higher redshifts ($0.1<z<1.2$), including 12 new to the literature, using X-ray data from the XMM Cluster Survey (XCS) and weak lensing data from the Canada France Hawaii Telescope Lensing Survey (CFHTLenS). Using a multi-parameter MCMC analysis, we constrain the two chameleon gravity parameters ($\beta$ and $\phi_{\infty}$). Our fits are consistent with general relativity, not requiring a fifth force. In the special case of $f(R)$ gravity (where $\beta = \sqrt{1/6}$), we set an upper limit on the background field amplitude today of $|f_{\rm{R0}}| < 6 \times 10^{-5}$ (95% CL). This is one of the strongest constraints to date on $|f_{\rm{R0}}|$ on cosmological scales. We hope to improve this constraint in future by extending the study to hundreds of clusters using data from the Dark Energy Survey.

We have analyzed XMM-Newton, Chandra, and Suzaku observations of Kepler's supernova remnant (SNR) to investigate the properties of both the SN ejecta and the circumstellar medium (CSM). For comparison, we have also analyzed two similarly-aged, ejecta-dominated SNRs: Tycho's SNR, thought to be the remnant of a typical Type Ia SN, and SNR 0509-67.5 in the Large Magellanic Cloud, thought to be the remnant of an overluminous Type Ia SN. By simply comparing the X-ray spectra, we find that line intensity ratios of iron-group elements (IGE) to intermediate-mass elements (IME) for Kepler's SNR and SNR 0509-67.5 are much higher than those for Tycho's SNR. We therefore argue that Kepler is the product of an overluminous Type Ia SN. This inference is supported by our spectral modeling, which reveals the IGE and IME masses respectively to be ~0.95 M_sun and ~0.12 M_sun (Kepler's SNR), ~0.75 M_sun and ~0.34 M_sun (SNR 0509-67.5), and ~0.35 M_sun and ~0.70 M_sun (Tycho's SNR). We find that the CSM component in Kepler's SNR consists of tenuous diffuse gas (~0.3 M_sun) present throughout the entire remnant, plus dense knots (~0.035 M_sun). Both of these components have an elevated N abundance (N/H ~ 4 times the solar value), suggesting that they originate from CNO-processed material from the progenitor system. The mass of the diffuse CSM allows us to infer the pre-SN mass-loss rate to be ~1.5e-5 (V_w/10 km/s) M_sun/yr, in general agreement with results from recent hydrodynamical simulations. Since the dense knots have slow proper motions and relatively small ionization timescales, they were likely located a few pc away from the progenitor system. Therefore, we argue that Kepler's SN was an overluminous event that started to interact with massive CSM a few hundred years after the explosion. This supports the possible link between overluminous SNe and the so-called "Ia-CSM" SNe.

HESS J1841$-$055 is a diffuse unidentified gamma-ray source with the size of $\sim\,1.3^{\circ}\times\,1^{\circ}$. No conclusive counterpart in other wavelengths has so far detected. To search for X-rays responsible for the TeV emission, the Suzaku observations were conducted, which covered a half region of the HESS source. In the soft band (0.5-2.0 keV), we discovered a diffuse emission, Suzaku J1840.2$-$0552, with the size of $\sim10'$. Since its spectrum was fitted by an optically thin thermal plasma model, Suzaku J1840.2$-$0552 is likely to be a supernova remnant. We also discovered an extended source, Suzaku J1840.2$-$0544, in the hard band (2.0-8.0 keV) with an emission line at 6.1 keV. From the spectral feature and large interstellar absorption, this source is likely to be a cluster of galaxies behind the Galactic plane at the red-shift of $\sim$0.09. The other diffuse source spatially overlaps with the SNR candidate G26.6$-$0.2, which shows a non-thermal dominant spectrum. Since no other candidate is found in the hard X-ray band, we infer that these largely extended sources could be possible counterparts of HESS J1841$-$055.

A common idea for the origin of the Galactic diffuse X-ray emission, particularly that of the iron lines from neutral and highly ionized atoms, is a superposition of many cataclysmic variables and coronally active binaries. In this scenario, the flux should symmetrically distribute between the east and west on the plane with respect to Sagittarius A$^*$ because the stellar mass distribution determined by infrared observations is nearly symmetric. This symmetry is confirmed for the highly ionized iron line as well as the continuum emission. However, a clear excess of the neutral iron line in the near east of the Galactic center compared to the near-west side is found. The flux distribution of the excess emission well correlates with molecular column density. The X-ray spectrum of the excess emission is described by a power-law continuum plus a 6.4 keV line with the large equivalent width of $\sim1.3$ keV, which is hardly explained by the low-energy electron bombardment scenario. The longitudinal and latitudinal distribution of the excess emission disfavors the X-ray irradiation, neither by Sagittarius A$^*$ nor by nearby X-ray binaries. Then the low-energy proton bombardment is the most probable origin although the high energy density $\sim 80$ eV cm$^{-3}$ in 0.1-1000 MeV is required and there is no conventional proton source in the vicinity.

In the approach of the effective field theory of modified gravity, we derive the equations of motion for linear perturbations in the presence of a barotropic perfect fluid on the flat isotropic cosmological background. In a simple version of Gleyzes-Langlois-Piazza-Vernizzi (GLPV) theories, which is the minimum extension of Horndeski theories, we show that a slight deviation of the tensor propagation speed squared $c_{\rm t}^2$ from 1 generally leads to the large modification to the propagation speed squared $c_{\rm s}^2$ of a scalar degree of freedom $\phi$. This problem persists whenever the kinetic energy $\rho_X$ of the field $\phi$ is much smaller than the background energy density $\rho_m$, which is the case for most of dark energy models in the asymptotic past. Since the scaling solution characterized by the constant ratio $\rho_X/\rho_m$ is one way out for avoiding such a problem, we study the evolution of perturbations for a scaling dark energy model in the framework of GLPV theories in the Jordan frame. Provided the oscillating mode of scalar perturbations is fine-tuned so that it is initially suppressed, the anisotropic parameter $\eta=-\Phi/\Psi$ between the two gravitational potentials $\Psi$ and $\Phi$ significantly deviates from 1 for $c_{\rm t}^2$ away from 1. For other general initial conditions, the deviation of $c_{\rm t}^2$ from 1 gives rise to the large oscillation of $\Psi$ with the frequency related to $c_{\rm s}^2$. In both cases, the model can leave distinct imprints for the observations of CMB and weak lensing.

Recent discoveries of recombining plasmas (RPs) in supernova remnants (SNRs) have dramatically changed our understanding of SNR evolution. To date, a dozen of RP SNRs have been identified in the Galaxy. Here we present Suzaku deep observations of four SNRs in the Large Magellanic Cloud (LMC), N49, N49B, N23, and DEM L71, for accurate determination of their plasma state. Our uniform analysis reveals that only N49 is in the recombining state among them, which is the first robust discovery of a RP from an extra-galactic SNR. Given that RPs have been identified only in core-collapse SNRs, our result strongly suggests a massive star origin of this SNR. On the other hand, no clear evidence for a RP is confirmed in N23, from which detection of recombination lines and continua was previously claimed. Comparing the physical properties of the RP SNRs identified so far, we find that all of them are categorized into the "mixed-morphology" class and interacting with surrounding molecular clouds. This might be a key to solve formation mechanisms of the RPs.

The basic building block for Lorentz invariant and ghost free massive gravity is the square root of the combination $g^{-1}\eta\,$, where $g^{-1}$ is the inverse of the physical metric and $\eta$ is a reference metric. Since the square root of a matrix is not uniquely defined, it is possible to have physically inequivalent potentials corresponding to different branches. We show that around Minkowski background the only perturbatively well defined branch is the potential proposed by de Rham, Gabadadze and Tolley. On the other hand, if Lorentz symmetry is broken spontaneously, other potentials exist with a standard perturbative expansion. We show this explicitly building new Lorentz invariant, ghost-free massive gravity potentials for theories that in the background preserve rotational invariance, but break Lorentz boosts.

We study cosmology in the bigravity formulation of the dRGT model where matter couples to both metrics. At linear order in perturbation theory two mass scales emerge: an hard one from the dRGT potential, and an environmental dependent one from the coupling of bigravity with matter. At early time, the dynamics is dictated by the second mass scale which is of order of the Hubble scale. The set of gauge invariant perturbations that couples to matter follow closely the same behaviour as in GR. The remaining perturbations show no issue in the scalar sector, while problems arise in the tensor and vector sectors. During radiation domination, a tensor mode grows power-like at super-horizon scales. More dangerously, the only propagating vector mode features an exponential instability on sub-horizon scales. We discuss the consequences of such instabilities and speculate on possible ways to deal with them.

Despite decades of intense efforts, many fundamental aspects of Type Ia supernova (SNe Ia) remain elusive. One of the major open questions is whether the mass of the exploding white dwarf (WD) is close to the Chandrasekhar limit. Here we report the detection of strong K-shell emission from stable Fe-peak elements in the Suzaku X-ray spectrum of the Type Ia supernova remnant (SNR) 3C 397. The high Ni/Fe and Mn/Fe mass ratios (0.11-0.24 and 0.018-0.033, respectively) in the hot plasma component that dominates the K-shell emission lines indicate a degree of neutronization in the SN ejecta which can only be achieved by electron captures in the dense cores of exploding WDs with a near-Chandrasekhar mass. This suggests a single-degenerate origin for 3C 397, since Chandrasekhar mass progenitors are expected naturally if the WD accretes mass slowly from a companion. Together with other results supporting the double-degenerate scenario, our work adds to the mounting evidence that both progenitor channels make a significant contribution to the SN Ia rate in star-forming galaxies.

The joint JAXA/NASA ASTRO-H mission is the sixth in a series of highly successful X-ray missions developed by the Institute of Space and Astronautical Science (ISAS), with a planned launch in 2015. The ASTRO-H mission is equipped with a suite of sensitive instruments with the highest energy resolution ever achieved at E > 3 keV and a wide energy range spanning four decades in energy from soft X-rays to gamma-rays. The simultaneous broad band pass, coupled with the high spectral resolution of Delta E < 7 eV of the micro-calorimeter, will enable a wide variety of important science themes to be pursued. ASTRO-H is expected to provide breakthrough results in scientific areas as diverse as the large-scale structure of the Universe and its evolution, the behavior of matter in the gravitational strong field regime, the physical conditions in sites of cosmic-ray acceleration, and the distribution of dark matter in galaxy clusters at different redshifts.

The most characteristic high-energy phenomena in the Galactic center (GC) region is the presence of strong K-shell emission lines from highly ionized Si, S, Ar, Ca, Fe and Ni, which form the Galactic Center X-ray Emission (GCXE). These multiple lines suggest that the GCXE is composed of at least two plasmas with temperatures of ~1 and ~7 keV. The GCXE also exhibits the K-shell lines from neutral Si, S, Ar, Ca, Fe and Ni atoms. A debatable issue is the origin of the GCXE plasma; whether it is a diffuse plasma or integrated emission of many unresolved point sources such as cataclysmic variables and active binaries. Detailed spectroscopy for these lines may provide a reliable picture of the GCXE plasma. The origin of the K-shell lines from neutral atoms is most likely the fluorescence by X-rays from a putative past flare of Sgr A*. Therefore ASTRO-H may provide unprecedented data for the past light curve of Sgr A*. All these lines may provide key information for the dynamics of the GCXE, using possible Doppler shift and/or line broadening. This paper overviews these line features and the previous interpretation of their origin. We propose extended or revised science with the ASTRO-H observations of some select objects in the GC region.

Thanks to the unprecedented spectral resolution and sensitivity of the Soft X-ray Spectrometer (SXS) to soft thermal X-ray emission, ASTRO-H will open a new discovery window for understanding young, ejecta-dominated, supernova remnants (SNRs). In particular we study how ASTRO-H observations will address, comprehensively, three key topics in SNR research: (1) using abundance measurements to unveil SNR progenitors, (2) using spatial and velocity distribution of the ejecta to understand supernova explosion mechanisms, (3) revealing the link between the thermal plasma state of SNRs and the efficiency of their particle acceleration.

We report on the Suzaku results of the mixed-morphology supernova remnant (SNR) G290.1$-$0.8 (MSH 11-61A). The SNR has an asymmetric structure extended to the southeast and the northwest. In the X-ray spectra of the center and the northwest regions, we discover recombining plasma features with the strong Si Ly$\alpha$ and radiative recombination continuum at $\sim$ 2.7 keV. These features are the most significant in the northwest region, and the spectra are well-reproduced with a recombining plasma of $kT_{\rm e} = 0.5$ keV. Whereas the spectra of other regions are expressed by an ionizing plasma of $kT_{\rm e} = 0.6$ keV. The recombining plasma has over-solar abundances, while the ionizing plasma has roughly solar abundances. Hence they are likely ejecta and interstellar medium (ISM) origin, respectively. The recombining plasma in the northwest of G290.1$-$0.8 would be generated by a break-out of the supernova ejecta from a high density circumstellar medium to a low density ISM.

We analyse modelling techniques for the large-scale structure formed in scalar-tensor theories of constant Brans-Dicke parameter which match the concordance model background expansion history and produce a chameleon suppression of the gravitational modification in high-density regions. Thereby, we use a mass and environment dependent chameleon spherical collapse model, the Sheth-Tormen halo mass function and linear halo bias, the Navarro-Frenk-White halo density profile, and the halo model. Furthermore, using the spherical collapse model, we extrapolate a chameleon mass-concentration scaling relation from a LCDM prescription calibrated to N-body simulations. We also provide constraints on the model parameters to ensure viability on local scales. We test our description of the halo mass function and nonlinear matter power spectrum against the respective observables extracted from large-volume and high-resolution N-body simulations in the limiting case of f(R) gravity, corresponding to a vanishing Brans-Dicke parameter. We find good agreement between the two; the halo model provides a good qualitative description of the shape of the relative enhancement of the f(R) matter power spectrum with respect to LCDM caused by the extra attractive gravitational force but fails to recover the correct amplitude. Introducing an effective linear power spectrum in the computation of the two-halo term to account for an underestimation of the chameleon suppression at intermediate scales in our approach, we accurately reproduce the measurements from the N-body simulations.

In the context of a cubic Galileon model in which the Vainshtein mechanism suppresses the scalar field interactions with matter, we study low-density stars with slow rotation and static relativistic stars. We develop an expansion scheme to find approximated solutions inside the Vainshtein radius, and show that deviations from General Relativity (GR), while considering rotation, are also suppressed by the Vainshtein mechanism. In a quadratic coupling model, in which the scalarisation effect can significantly enhance deviations from GR in normal scalar tensor gravity, the Galileon term successfully suppresses the large deviations away from GR. Moreover, using a realistic equation of state, we construct solutions for a relativistic star, and show that deviations from GR are more suppressed for higher density objects. However, we found that the scalar field solution ceases to exist above a critical density, which roughly corresponds to the maximum mass of a neutron star. This indicates that, for a compact object described by a polytropic equation of state, the configuration that would collapse into a black hole cannot support a non-trivial scalar field.

Redshift-space distortions (RSD) offer an attractive method to measure the growth of cosmic structure on large scales, and combining with the measurement of the cosmic expansion history, it can be used as cosmological tests of gravity. With the advent of future galaxy redshift surveys aiming at precisely measuring the RSD, an accurate modeling of RSD going beyond linear theory is a critical issue in order to detect or disprove small deviations from general relativity (GR). While several improved models of RSD have been recently proposed based on the perturbation theory (PT), the framework of these models heavily relies on GR. Here, we put forward a new PT prescription for RSD in general modified gravity models. As a specific application, we present theoretical predictions of the redshift-space power spectra in f(R) gravity model, and compare them with N-body simulations. Using the PT template that takes into account the effects of both modifications of gravity and RSD properly, we successfully recover the fiducial model parameter in N-body simulations in an unbiased way. On the other hand, we found it difficult to detect the scale dependence of the growth rate in a model-independent way based on GR templates.

In this note, we study cosmic microwave background (CMB) constraints on primordial non-Gaussianity in DBI galileon models in which an induced gravity term is added to the Dirac-Born-Infeld (DBI) action. In this model, the non-Gaussianity of orthogonal shape can be generated. We provide a relation between theoretical parameters and orthogonal/equilateral non-linear parameters using the Fisher matrix approach for the CMB bispectrum. In doing so, we include the effect of the CMB transfer functions and experimental noise properties by employing the recently developed \textsf{SONG} code. The relation is also shown in the language of effective theory so that it can be applied to general single-field models. Using the bispectrum Fisher matrix and the central values for equilateral and orthogonal non-Gaussianities found by the Planck temperature survey, we provide forecasts on the theoretical parameters of the DBI galileon model. We consider the upcoming Planck polarisation data and the proposed post-Planck experiments COrE and PRISM. We find that Planck polarisation measurements may provide a hint for a non-canonical sound speed at the 68% confidence level. COrE and PRISM will not only confirm a non-canonical sound speed but also exclude the conventional DBI inflation model at more than the 95% and 99% confidence level respectively, assuming that the central values will not change. This indicates that improving constraints on non-Gaussianity further by future CMB experiments is invaluable to constrain the physics of the early universe.

Recent X-ray study of middle-aged supernova remnants (SNRs) reveals strong radiative recombination continua (RRCs) associated with overionized plasmas, of which the origin still remains uncertain. We report our discovery of an RRC in the middle-aged SNR 3C 391. If the X-ray spectrum is fitted with a two-temperature plasma model in collisional ionization equilibrium (CIE), residuals of Si XIV Ly alpha line at 2.006 keV, S XVI Ly alpha line at 2.623 keV and the edge of RRC of Si XIII at 2.666 keV are found. The X-ray spectrum is better described by a composite model consisting of a CIE plasma and a recombining plasma (RP). The abundance pattern suggests that the RP is associated to the ejecta from a core-collapse supernova with a progenitor star of 15 solar mass. There is no significant difference of the recombining plasma parameters between the southeast region and the northwest region surrounded by dense molecular clouds. We also find a hint of Fe I K alpha line at 6.4 keV (~2.4 sigma detection) from the southeast region of the SNR.

An elongated X-ray source with a strong K-shell line from He-like iron (Fe XXVI) is found at (RA, Dec)_{J2000.0}=(17h44m00s.0, -29D13'40''.9) in the Galactic center region. The position coincides with the X-ray thread, G359.55+0.16, which is aligned with the radio non-thermal filament. The X-ray spectrum is well fitted with an absorbed thin thermal plasma (apec) model. The best-fit temperature, metal abundance, and column density are 4.1^{+2.7}_{-1.8} keV, 0.58^{+0.41}_{-0.32} solar, and 6.1^{+2.5}_{-1.3}x10^{22} cm^{-2}, respectively. These values are similar to those of the largely extended Galactic center X-ray emission.

We study the behaviour of weak gravitational fields in the 6D Cascading DGP model using a bulk-based approach. To deal with the ambiguity in the thin limit of branes of codimension higher than one, we consider a specific regularization of the internal structure of the branes where the 5D brane can be considered thin with respect to the 4D one. We consider the solutions corresponding to pure tension sources on the 4D brane, and study perturbations at first order around these background solutions. We adopt a 4D scalar-vector-tensor decomposition, and focus on the scalar sector of perturbations. We show that, in a suitable 4D limit, the trace part of the 4D metric perturbations obeys a decoupled equation which suggests that it is a ghost for background tensions smaller than a critical tension, while it is a healthy field otherwise. We give a geometrical interpretation of the existence of the critical tension and of the reason why the relevant field is a ghost or not depending on the background tension. We however find a value of the critical tension which is different from the one already found in the literature. Differently from the results in the literature, our analysis implies that, choosing the background tension suitably, we can construct ghost-free models for any value of the free parameters of the theory. We suggest that the difference lies in the procedure used to evaluate the pillbox integration across the codimension-2 brane. We confirm the validity of our analysis by performing numerically the integration in a particular case where the solution inside the thick cod-2 brane is known exactly. We stress that the singular structure of the perturbation fields in the nested branes set-ups is very subtle, and that great care has to be taken when deriving the codimension-2 junction conditions.

Based on the multi-point propagator expansion, we present resummed perturbative calculations for cosmological power spectra and correlation functions in the context of modified gravity. In a wide class of modified gravity models that have a screening mechanism to recover general relativity (GR) on small scales, we apply the eikonal approximation to derive the governing equation for resummed propagator that partly includes the non-perturbative effect in the high-k limit. The resultant propagator in the high-k limit contains the new corrections arising from the screening mechanism as well as the standard exponential damping. We explicitly derive the expression for new high-k contributions in specific modified gravity models, and find that in the case of f(R) gravity for a currently constrained model parameter, the corrections are basically of the sub-leading order and can be neglected. Thus, in f(R) gravity, similarly to the GR case, we can analytically construct the regularized propagator that reproduces both the resummed high-k behavior and the low-k results computed with standard perturbation theory, consistently taking account of the nonlinear modification of gravity valid at large scales. With the regularized multi-point propagators, we give predictions for power spectrum and correlation function at one-loop order, and compare those with N-body simulations in f(R) gravity model. As an important application, we also discuss the redshift-space distortions and compute the anisotropic power spectra and correlation functions.

We investigate the dependence of the Vainshtein screening mechanism on the cosmic web morphology of both dark matter particles and halos as determined by ORIGAMI. Unlike chameleon and symmetron screening, which come into effect in regions of high density, Vainshtein screening instead depends on the dimensionality of the system, and screened bodies can still feel external fields. ORIGAMI is well-suited to this problem because it defines morphologies according to the dimensionality of the collapsing structure and does not depend on a smoothing scale or density threshold parameter. We find that halo particles are screened while filament, wall, and void particles are unscreened, and this is independent of the particle density. However, after separating halos according to their large scale morphological environment, we find no difference in the screening properties of halos in filaments versus halos in clusters. We find that the fifth force enhancement of dark matter particles in halos is greatest well outside the virial radius. We confirm the theoretical expectation that even if the internal field is suppressed by the Vainshtein mechanism, the object still feels the fifth force generated by the external fields, by measuring peculiar velocities and velocity dispersions of halos. Finally, we investigate the morphology and gravity model dependence of halo spins, concentrations, and shapes.

We report the first results from high-statistics observation of the 6.4-keV line in the region of $l= +1.5^\circ$ to $+3.5^\circ$ (hereafter referred to as GC East), with the goal to uncover the origin of the Galactic ridge X-ray emission (GRXE). By comparing this data with that from the previous observations in the region $l=-1.5^\circ$ to $-3.5^\circ$ (hereafter referred to as GC West), we discovered that the 6.4-keV line is asymmetrically distributed with respect to the Galactic center, whereas the 6.7-keV line is symmetrically distributed. The distribution of the 6.4-keV line follows that of $^{13}$CO and its flux is proportional to the column density of the molecular gas. This correlation agrees with that seen between the 6.4-keV line and the cold interstellar medium (ISM) (H$_{\rm I}$ $+$ H$_2$) in the region $|l|>4^\circ$. This result suggests that the 6.4-keV emission is diffuse fluorescence from the cold ISM not only in GC East and West but also in the entire Galactic plane. This observational result suggests that the surface brightness of the 6.4-keV line is proportional to the column density of the cold ISM in the entire Galactic plane. For the ionizing particles, we consider X-rays and low energy cosmic-ray protons and electrons .

We study static, spherically symmetric solutions in a recently proposed ghost-free model of non-linear massive gravity. We focus on a branch of solutions where the helicity-0 mode can be strongly coupled within certain radial regions, giving rise to the Vainshtein effect. We truncate the analysis to scales below the gravitational Compton wavelength, and consider the weak field limit for the gravitational potentials, while keeping all non-linearities of the helicity-0 mode. We determine analytically the number and properties of local solutions which exist asymptotically on large scales, and of local (inner) solutions which exist on small scales. We find two kinds of asymptotic solutions, one of which is asymptotically flat, while the other one is not, and also two types of inner solutions, one of which displays the Vainshtein mechanism, while the other exhibits a self-shielding behaviour of the gravitational field. We analyse in detail in which cases the solutions match in an intermediate region. The asymptotically flat solutions connect only to inner configurations displaying the Vainshtein mechanism, while the non asymptotically flat solutions can connect with both kinds of inner solutions. We show furthermore that there are some regions in the parameter space where global solutions do not exist, and characterise precisely in which regions of the phase space the Vainshtein mechanism takes place.

We present the results of a Suzaku observation of Abell 1674, an optically very rich (richness class 3) cluster cataloged as z = 0.1066. We discover the He-like Fe K-shell line from the central region for the first time, and find that the X-ray spectrum yields a high redshift of 0.215^{+0.007}_{-0.006}. On the other hand, the spectrum of the southwestern region is fitted with a redshift of 0.11 \pm 0.02 by the He-like Fe L-shell lines, consistent with the optically determined value. The gas temperature, metal abundance, and core radius of the X-ray emission in the central region are 3.8 \pm 0.2 keV, 0.20 \pm 0.05 Z_solar and 450 \pm 40 kpc, respectively, while those in the southwestern region are 2.0 \pm 0.2 keV, 0.41^{+0.17}_{-0.13} Z_solar and 220^{+90}_{-70} kpc, respectively. These parameters are typical for clusters. We thus conclude that Abell 1674 consists of two independent clusters, A1674-C at z ~0.22 and A1674-SW at z ~0.11, overlapping along the line of sight. The X-ray luminosities of A1674-C within r = 2 Mpc is 15.9 \pm 0.6 \times 10^{43} erg s^{-1} in the 0.1-2.4 keV energy band, while that for A1674-SW is 1.25 \pm 0.07 \times 10^{43} erg s^{-1}. Both are consistent with those expected from the L-T relation of clusters within a factor of 2. This is another support for our conclusion.

We compute the bispectrum induced in the cosmic microwave background (CMB) temperature and polarisation by the evolution of the primordial density perturbations using the second-order Boltzmann code SONG. We show that adding polarisation increases the signal-to-noise ratio by a factor four with respect to temperature alone and we estimate the observability of this intrinsic bispectrum and the bias it induces on measurements of primordial non-Gaussianity. When including all physical effects except the late-time non-linear evolution, we find for the intrinsic bispectrum a signal-to-noise of $S/N=3.8,\,2.9,\,1.6$ and $0.5$ for, respectively, an ideal experiment with an angular resolution of $\ell_\text{max}=3000$, the proposed CMB surveys PRISM and COrE, and Planck's polarised data; the bulk of this signal comes from the $E$-polarisation and from squeezed configurations. We discuss how CMB lensing is expected to reduce these estimates as it suppresses the bispectrum for squeezed configurations and contributes to the noise in the estimator. We find that the presence of the intrinsic bispectrum will bias a measurement of primordial non-Gaussianity of local type by $f_{NL}^\text{bias}=0.66$ for an ideal experiment with $\ell_\text{max}=3000$. Finally, we verify the robustness of our results by reproducing the analytical approximation for the squeezed-limit bispectrum in the general polarised case.

We study the behaviour of weak gravitational fields in models where a 4D brane is embedded inside a 5D brane equipped with induced gravity, which in turn is embedded in a 6D spacetime. We consider a specific regularization of the branes internal structures where the 5D brane can be considered thin with respect to the 4D one. We find exact solutions corresponding to pure tension source configurations on the thick 4D brane, and study perturbations at first order around these background solutions. To perform the perturbative analysis, we adopt a bulk-based approach and we express the equations in terms of gauge invariant and master variables using a 4D scalar-vector-tensor decomposition. We then propose an ansatz on the behaviour of the perturbation fields when the thickness of the 4D brane goes to zero, which corresponds to configurations where gravity remains finite everywhere in the thin limit of the 4D brane. We study the equations of motion using this ansatz, and show that they give rise to a consistent set of differential equations in the thin limit, from which the details of the internal structure of the 4D brane disappear. We conclude that the thin limit of the "ribbon" 4D brane inside the (already thin) 5D brane is well defined (at least when considering first order perturbations around pure tension configurations), and that the gravitational field on the 4D brane remains finite in the thin limit. We comment on the crucial role of the induced gravity term on the 5D brane.

We estimate the B-polarisation induced in the Cosmic Microwave Background by the non-linear evolution of density perturbations. Using the second-order Boltzmann code SONG, our analysis incorporates, for the first time, all physical effects at recombination. We also include novel contributions from the redshift part of the Boltzmann equation and from the bolometric definition of the temperature in the presence of polarisation. The remaining line-of-sight terms (lensing and time-delay) have previously been studied and must be calculated non-perturbatively. The intrinsic B-mode polarisation is present independent of the initial conditions and might contaminate the signal from primordial gravitational waves. We find this contamination to be comparable to a primordial tensor-to-scalar ratio of $r\simeq10^{-7}$ at the angular scale $\ell\simeq100\,$, where the primordial signal peaks, and $r\simeq 5 \cdot 10^{-5}$ at $\ell\simeq700\,$, where the intrinsic signal peaks. Therefore, we conclude that the intrinsic B-polarisation from second-order effects is not likely to contaminate future searches of primordial gravitational waves.

We present Suzaku results of the two Galactic supernova remnants (SNRs), G350.1-0.3 and G349.7+0.2. We find Al and Ni K alpha lines from both the SNRs for the first time, in addition to previously detected K-shell lines of Mg, Si, S, Ar, Ca and Fe. The spectra are well described by two optically thin thermal plasmas: a low-temperature (low-kT) plasma in collisional ionization equilibrium and a high-temperature (high-kT) plasma in non-equilibrium ionization. Since the low-kT plasma has solar metal abundances, it is thought to be of interstellar medium origin. The high-kT plasma has super-solar abundances, hence it is likely to be of ejecta origin. The abundance patterns of the ejecta components are similar to those of core-collapse supernovae with the progenitor mass of ~15-25 M_solar for G350.1-0.3 and ~35-40 M_solar for G349.7+0.2. We find extremely high abundances of Ni compared to Fe (Z_Ni/Z_Fe ~8). Based on the measured column densities between the SNRs and the near sky background, we propose that G350.1-0.3 and G349.7+0.2 are located at the distance of 9+/-3 kpc and 12+/-5 kpc, respectively. Then the ejecta masses are estimated to be ~13 M_solar and ~24 M_solar for G350.1-0.3 and G349.7+0.2, respectively. These values are consistent with the progenitor mass of ~15-25 M_solar and ~35-40 M_solar for G350.1-0.3 and G349.7+0.2, respectively.

We present deep observations of the Galactic supernova remnant IC 443 with the {\it Suzaku X-ray satellite}. We find prominent K-shell lines from iron and nickel, together with a triangle residual at 8--10~keV, which corresponds to the energy of the radiative recombination continuum (RRC) of He-like iron. In addition, the wavy residuals have been seen at $\sim$5.1 and $\sim$5.5~keV. We confirm that the residuals show the first enhanced RRCs of He- and H-like calcium found in supernova remnants. These facts provide robust evidence for the recombining plasma. We reproduce the plasma in the 3.7--10~keV band using a recombining plasma model at the electron temperature 0.65~keV. The recombination parameter $n_{\rm e}t$ ($n_{\rm e}$ is electron density and $t$ is elapsed time after formation of a recombining plasma) and abundances of iron and nickel are strongly correlated, and hence the errors are large. On the other hand, the ratio of nickel to iron relative to the solar abundances is well constrained to 11$^{+4}_{-3}$ (1$\sigma$). A possibility is that the large abundance ratio is a result of an asymmetric explosion of the progenitor star.

We propose a novel method to test the gravitational interactions in the outskirts of galaxy clusters. When gravity is modified, this is typically accompanied by the introduction of an additional scalar degree of freedom, which mediates an attractive fifth force. The presence of an extra gravitational coupling, however, is tightly constrained by local measurements. In chameleon modifications of gravity, local tests can be evaded by employing a screening mechanism that suppresses the fifth force in dense environments. While the chameleon field may be screened in the interior of the cluster, its outer region can still be affected by the extra force, introducing a deviation between the hydrostatic and lensing mass of the cluster. Thus, the chameleon modification can be tested by combining the gas and lensing measurements of the cluster. We demonstrate the operability of our method with the Coma cluster, for which both a lensing measurement and gas observations from the X-ray surface brightness, the X-ray temperature, and the Sunyaev-Zel'dovich effect are available. Using the joint observational data set, we perform a Markov chain Monte Carlo analysis of the parameter space describing the different profiles in both the Newtonian and chameleon scenarios. We report competitive constraints on the chameleon field amplitude and its coupling strength to matter. In the case of f(R) gravity, corresponding to a specific choice of the coupling, we find an upper bound on the background field amplitude of |f_{R0}|<6*10^{-5}, which is currently the tightest constraint on cosmological scales.

We calculate the shape correlations between trispectra in various equilateral non-Gaussian models, including DBI inflation, ghost inflation and Lifshitz scalars, using the full trispectrum as well as the reduced trispectum. We find that most theoretical models are distinguishable from the shapes of primordial trispectra except for several exceptions where it is difficult to discriminate between the models, such as single field DBI inflation and a Lifshitz scalar model. We introduce an estimator for the amplitude of the trispectrum, $g_{\rm NL} ^{equil}$ and relate it to model parameters in various models. Using constraints on $g_{\rm NL} ^{equil}$ from WMAP5, we give constraints on the model parameters.

We report the Suzaku/XIS results of the Galactic oxygen-rich supernova remnant (SNR), G292.0+1.8, a remnant of a core-collapse supernova. The X-ray spectrum of G292.0+1.8 consists of two type plasmas, one is in collisional ionization equilibrium (CIE) and the other is in non-equilibrium ionization (NEI). The CIE plasma has nearly solar abundances, and hence would be originated from the circumstellar and interstellar mediums. The NEI plasma has super-solar abundances, and the abundance pattern indicates that the plasma originates from the supernova ejecta with a main sequence of 30-35 Msolar. Iron K-shell line at energy of 6.6 keV is detected for the first time in the NEI plasma.

Gravitational vector degrees of freedom typically arise in many examples of modified gravity models. We start to systematically explore their role in these scenarios, studying the effects of coupling gravitational vector and scalar degrees of freedom. We focus on set-ups that enjoy a Galilean symmetry in the scalar sector and an Abelian gauge symmetry in the vector sector. These symmetries, together with the requirement that the equations of motion contain at most two space-time derivatives, only allow for a small number of operators in the Lagrangian for the gravitational fields. We investigate the role of gravitational vector fields for two broad classes of phenomena that characterize modified gravity scenarios. The first is self-acceleration: we analyze in general terms the behavior of vector fluctuations around self-accelerating solutions, and show that vanishing kinetic terms of vector fluctuations lead to instabilities on cosmological backgrounds. The second phenomenon is the screening of long range fifth forces by means of Vainshtein mechanism. We show that if gravitational vector fields are appropriately coupled to a spherically symmetric source, they can play an important role for defining the features of the background solution and the scale of the Vainshtein radius. Our general results can be applied to any concrete model of modified gravity, whose low-energy vector and scalar degrees of freedom satisfy the symmetry requirements that we impose.

Current constraints on f(R) gravity from the large-scale structure are at the verge of penetrating into a region where the modified forces become nonlinearly suppressed. For a consistent treatment of observables at these scales, we study cluster quantities produced in chameleon and linearized Hu-Sawicki f(R) gravity dark matter N-body simulations. We find that the standard Navarro-Frenk-White halo density profile and the radial power law for the pseudo phase-space density provide equally good fits for f(R) clusters as they do in the Newtonian scenario. We give qualitative arguments for why this should be the case. For practical applications, we derive analytic relations, e.g., for the f(R) scalar field, the gravitational potential, and the velocity dispersion as seen within the virialized clusters. These functions are based on three degrees of freedom fitted to simulations, i.e., the characteristic density, scale, and velocity dispersion. We further analyze predictions for these fitting parameters from the gravitational collapse and the Jeans equation, which are found to agree well with the simulations. Our analytic results can be used to consistently constrain chameleon f(R) gravity with future observations on virialized cluster scales without the necessity of running a large number of simulations.

We develop a numerical algorithm to solve the high-order nonlinear derivative-coupling equation associated with the quartic Galileon model, and implement it in a modified version of the RAMSES N-body code to study the effect of the Galileon field on the large-scale matter clustering. The algorithm is tested for several matter field configurations with different symmetries, and works very well. This enables us to perform the first simulations for a quartic Galileon model which provides a good fit to the cosmic microwave background (CMB) anisotropy, supernovae and baryonic acoustic oscillations (BAO) data. Our result shows that the Vainshtein mechanism in this model is very efficient in suppressing the spatial variations of the scalar field. However, the time variation of the effective Newtonian constant caused by the curvature coupling of the Galileon field cannot be suppressed by the Vainshtein mechanism. This leads to a significant weakening of the strength of gravity in high-density regions at late times, and therefore a weaker matter clustering on small scales. We also find that without the Vainshtein mechanism the model would have behaved in a completely different way, which shows the crucial role played by nonlinearities in modified gravity theories and the importance of performing self-consistent N-body simulations for these theories.

We explore the signatures of quintessence and modified gravity theories in the relativistic description of galaxy clustering within a parametrized post-Friedmann framework. For this purpose, we develop a calibration method to consistently account for horizon-scale effects in the linear parametrized Post-Friedmann perturbations of minimally and nonminimally coupled scalar-tensor theories and test it against the full model-specific fluctuations. We further study the relativistic effects in galaxy clustering for the normal and self-accelerating branches of the Dvali-Gabadadze-Porrati braneworld model as well as for phenomenological modifications of gravity. We quantify the impact of modified gravity and dark energy models on galaxy clustering by computing the velocity-to-matter density ratio F, the velocity contribution R, and the potential contribution P and give an estimate of their detectability in future galaxy surveys. Our results show that, in general, the relativistic correction contains additional information on gravity and dark energy, which needs to be taken into account in consistent horizon-scale tests of departures from LCDM using the galaxy-density field.

On cosmological scales, observations of the cluster abundance currently place the strongest constraints on f(R) gravity. These constraints lie in the large-field limit, where the modifications of general relativity can correctly be modeled by setting the Compton wavelength of the scalar field to its background value. These bounds are, however, at the verge of penetrating into a regime, where the modifications become nonlinearly suppressed due to the chameleon mechanism and cannot be described by this linearized approximation. For future constraints based on observations subjected to cluster abundance, it is therefore essential to consistently model the chameleon effect. We analyze descriptions of the halo mass function in chameleon f(R) gravity using a mass- and environment-dependent spherical collapse model in combination with excursion set theory and phenomenological fits to N-body simulations in the LCDM and f(R) gravity scenarios. Our halo mass functions consistently incorporate the chameleon suppression and cosmological parameter dependencies, improving upon previous formalisms and providing an important extension to N-body simulations for the application in consistent tests of gravity with observables sensitive to the abundance of clusters.

We report Suzaku results for soft X-ray emission to the south of the Galactic center (GC). The emission (hereafter "GC South") has an angular size of ~42' x 16' centered at (l, b) ~ (0.0, -1.4), and is located in the largely extended Galactic ridge X-ray emission (GRXE). The X-ray spectrum of GC South exhibits emission lines from highly ionized atoms. Although the X-ray spectrum of the GRXE can be well fitted with a plasma in collisional ionization equilibrium (CIE), that of GC South cannot be fitted with a plasma in CIE, leaving hump-like residuals at ~2.5 and 3.5 keV, which are attributable to the radiative recombination continua of the K-shells of Si and S, respectively. In fact, GC South spectrum is well fitted with a recombination-dominant plasma model; the electron temperature is 0.46 keV while atoms are highly ionized (kT = 1.6 keV) in the initial epoch, and the plasma is now in a recombining phase at a relaxation scale (plasma density x elapsed time) of 5.3 x 10^11 s cm^-3. The absorption column density of GC South is consistent with that toward the GC region. Thus GC South is likely to be located in the GC region (~8 kpc distance). The size of the plasma, the mean density, and the thermal energy are estimated to be 97 pc x 37 pc, 0.16 cm^-3, and 1.6 x 10^51 erg, respectively. We discuss possible origins of the recombination-dominant plasma as a relic of past activity in the GC region.

We present results of the Suzaku/XIS observation around the radio supernova remnant (SNR) G12.0-0.1. No significant diffuse emission extending in or along the radio shell was observed. Instead two compact X-ray sources, Suzaku J181205-1835 and Suzaku J181210-1842, were found in or near G12.0-0.1. Suzaku J181205-1835 is located at the northwest of the radio shell of G12.0-0.1. The X-ray profile is slightly extended over the point spread function of the Suzaku telescope. The X-ray spectrum has no line-like structure and is well represented by a power-law model with a photon index of 2.2 and an absorption column of N_{H}=4.9x10^{22} cm^{-2}. The distances of Suzaku J181205-1835 and G12.0-0.1 are estimated from the absorption column and the Sigma-D relation, respectively, and are nearly the same with each other. These results suggest that Suzaku J181205-1835 is a candidate of a pulsar wind nebula associated with G12.0-0.1. From its location, Suzaku J181210-1842 would be unrelated object to G12.0-0.1. The X-ray profile is point-like and the spectrum is thin thermal emission with Fe K-lines at 6.4, 6.7, and 6.97 keV, similar to those of cataclysmic variables.

The discovery of cosmic acceleration has stimulated theorists to consider dark energy or modifications to Einstein's General Relativity as possible explanations. The last decade has seen advances in theories that go beyond smooth dark energy -- modified gravity and interactions of dark energy. While the theoretical terrain is being actively explored, the generic presence of fifth forces and dark sector couplings suggests a set of distinct observational signatures. This report focuses on observations that differ from the conventional probes that map the expansion history or large-scale structure. Examples of such novel probes are: detection of scalar fields via lab experiments, tests of modified gravity using stars and galaxies in the nearby universe, comparison of lensing and dynamical masses of galaxies and clusters, and the measurements of fundamental constants at high redshift. The observational expertise involved is very broad as it spans laboratory experiments, high resolution astronomical imaging and spectroscopy and radio observations. In the coming decade, searches for these effects have the potential for discovering fundamental new physics. We discuss how the searches can be carried out using experiments that are already under way or with modest adaptations of existing telescopes or planned experiments. The accompanying paper on the Growth of Cosmic Structure describes complementary tests of gravity with observations of large-scale structure.

We report the \Suzaku/XIS and HXD and \Chandra/ACIS-I results on the X-ray spectra of the Phoenix cluster at the redshift $z=0.596$. The spectrum of the intracluster medium (ICM) is well-reproduced with the emissions from a low temperature ($\sim3.0$\,keV and $\sim0.76$\,solar) and a high temperature ($\sim11$\,keV and $\sim0.33$\,solar) plasmas; the former is localized at the cluster core, while the latter distributes over the cluster. In addition to these ICM emissions, a strongly absorbed power-law component is found, which is due to an active galactic nucleus (AGN) in the cluster center. The absorption column density and unobscured luminosity of the AGN are $\sim3.2\times10^{23}$\,cm$^{-2}$ and $\sim4.7\times10^{45}$\,ergs\,s$^{-1}$ (2-10\,keV), respectively. Furthermore, a neutral iron (\ion{Fe}{1}) K-shell line is discovered for the first time with the equivalent width (EW) of $\sim150$\,eV at the rest frame. The column density and the EW of the \ion{Fe}{1} line are exceptionally large for such a high luminosity AGN, and hence the AGN is classified as a type 2 quasi-stellar object (QSO). We speculate that the significant fraction of the ICM cooled gas would be consumed to maintain the torus and to activate the type 2 QSO. The Phoenix cluster has a massive starburst in the central galaxy, indicating suppression in the cooling flow is less effective. This may be because the onset of the latest AGN feedback has occurred recently and it has not yet been effective. Alternatively, the AGN feedback is predominantly in radiative-mode not in kinetic-mode and the torus may work as a shield to reduce its effect.

We perform a principal component analysis to assess ability of future observations to measure departures from General Relativity in predictions of the Poisson and anisotropy equations on linear scales. In particular, we focus on how the measurements of redshift-space distortions (RSD) observed from spectroscopic galaxy redshift surveys will improve the constraints when combined with lensing tomographic surveys. Assuming a Euclid-like galaxy imaging and redshift survey, we find that adding the 3D information decreases the statistical uncertainty by a factor between 3 and 10 compared to the case when only observables from lensing tomographic surveys are used. We also find that the number of well-constrained modes increases by a factor between 3 and 7. Our study indicates the importance of joint galaxy imaging and redshift surveys such as SuMIRe and Euclid to give more stringent tests of the {\Lambda}CDM model and to distinguish between various modified gravity and dark energy models.

We study linear perturbations around time dependent spherically symmetric solutions in the Lambda_3 massive gravity theory, which self-accelerate in the vacuum. We find that the dynamics of the scalar perturbations depend on the coordinate choice for the background solutions. For particular choices of coordinates there is a symmetry enhancement, leaving no propagating scalar degrees of freedom at linear order in perturbations. In contrast, any other coordinate choice propagates a single scalar mode. We find that the Hamiltonian of this scalar mode is unbounded from below for all self-accelerating solutions, signalling an instability.

We report on Suzaku results of the two distinct regions in the Galactic supernova remnant G348.5+0.1: extended thermal X-rays ("soft diffuse") at the northeast region and non-thermal X-rays (CXOU J171419.8-383023) at the northwest region. The X-ray spectrum of the soft diffuse is fitted with neither an ionization equilibrium nor a non-equilibrium (ionizing) plasma model, leaving saw-teeth residuals in the 1.5-3 keV energy band. The residual structures can be produced when free electrons are recombined to the K-shells of highly ionized Mg and Si ions. In fact, the X-ray spectrum is nicely fitted with a recombination-dominant plasma model. We propose a scenario that the plasma in nearly full ionized state at high temperature quickly changed to a recombining phase due to selective cooling of electrons to a lower temperature of ~0.5 keV. The spectrum of CXOU J171419.8-383023 is well explained by a simple power-law model with a photon index of 1.9, nearly equal to the typical value of pulsar wind nebulae. Since the distance is estimated to be the same as that of the soft diffuse, we infer that both the soft diffuse and CXOU J171419.8-383023 are associated with the same object, SNR G348.5+0.1.

Starting from the general Horndeski action, we derive the most general effective theory for scalar perturbations around flat space that allows us to screen fifth forces via the Vainshtein mechanism. The effective theory is described by a generalization of the Galileon Lagrangian, which we use to study the stability of spherically symmetric configurations exhibiting the Vainshtein effect. In particular, we discuss the phenomenological consequences of a scalar-tensor coupling that is absent in the standard Galileon Lagrangian. This coupling controls the superluminality and stability of fluctuations inside the Vainshtein radius in a way that depends on the density profile of a matter source. Particularly we find that the vacuum solution is unstable due to this coupling.

We present results of the Galactic supernova remnant (SNR) G355.6-0.0 observed with Suzaku. We resolved the diffuse emission detected with ASCA into two objects, G355.6-0.0 and a point-like source, Suzaku J173530-3236. The X-ray emission from G355.6-0.0 exhibits a center-filled morphology within the radio shell. The X-ray spectrum is well represented by a thin thermal plasma model with enhanced metal abundances. The spatial and spectral properties imply that G355.6-0.0 is a member of the mixed-morphology SNRs. The NH value of ~ 6 x 10^{22} cm^{-2} supports that G355.6-0.0 is a distant SNR. Suzaku J173530-3236 exhibits a hard X-ray spectrum with a strong Fe emission line, similar to those of cataclysmic variables. The NH value of ~ 3 x 10^{22} cm^{-2} is smaller than that of G355.6-0.0, and hence Suzaku J173530-3236 is located at the near side of G355.6-0.0.

We present the results from deep X-ray observations (~400 ks in total) of SN 1006 by the X-ray astronomy satellite Suzaku. The thermal spectrum from the entire supernova remnant (SNR) exhibits prominent emission lines of O, Ne, Mg, Si, S, Ar, Ca, and Fe. The observed abundance pattern in the ejecta components is in good agreement with that predicted by a standard model of Type Ia supernovae (SNe). The spatially resolved analysis reveals that the distribution of the O-burning and incomplete Si-burning products (Si, S, and Ar) is asymmetric, while that of the C-burning products (O, Ne, and Mg) is relatively uniform in the SNR interior. The peak position of the former is clearly shifted by 5' (~3.2 pc at a distance of 2.2 kpc) to the southeast from the SNR's geometric center. Using the SNR age of ~1000 yr, we constrain the velocity asymmetry (in projection) of ejecta to be ~3100 km/s. The abundance of Fe is also significantly higher in the southeast region than in the northwest region. Given that the non-uniformity is observed only among the heavier elements (Si through Fe), we argue that SN 1006 originates from an asymmetric explosion, as is expected from recent multi-dimensional simulations of Type Ia SNe, although we cannot eliminate the possibility that an inhomogeneous ambient medium induced the apparent non-uniformity. Possible evidence for the Cr K-shell line and line broadening in the Fe K-shell emission is also found.

Models of inflation in a gravitational background with an anisotropic space-time scaling are studied. The background is a higher-dimensional Lifshitz throat with the anisotropy scaling $z\neq 1$. After the dimensional reduction, the four-dimensional general covariance is explicitly broken to a three-dimensional spatial diffeomorphism. As a result the cosmological perturbation theory in this set up with less symmetries have to be formulated. We present the consistent cosmological perturbation theory for this set up. We find that the effective four-dimensional gravitational wave perturbations propagate with a different speed than the higher dimensional gravitational excitations. Depending on the model parameters, for an observer inside the throat, the four-dimensional gravitational wave propagation can be superluminal. We also find that the Bardeen potential and the Newtonian potential are different. This can have interesting observational consequences for lensing and CMB fluctuations. Furthermore, we show that at the linearized level the inflaton field excitations vanish.

In the context of cosmological perturbation theory, we derive the second order Boltzmann equation describing the evolution of the distribution function of radiation without a specific gauge choice. The essential steps in deriving the Boltzmann equation are revisited and extended given this more general framework: i) the polarisation of light is incorporated in this formalism by using a tensor-valued distribution function; ii) the importance of a choice of the tetrad field to define the local inertial frame in the description of the distribution function is emphasized; iii) we perform a separation between temperature and spectral distortion, both for the intensity and for polarisation for the first time; iv) the gauge dependence of all perturbed quantities that enter the Boltzmann equation is derived, and this enables us to check the correctness of the perturbed Boltzmann equation by explicitly showing its gauge-invariance for both intensity and polarization. We finally discuss several implications of the gauge dependence for the observed temperature.

There has been a lot of research interest in modified gravity theories which utilise the Vainshtein mechanism to recover standard general relativity in regions with high matter density, such as the Dvali-Gabadadze-Porrati and Galileon models. The strong nonlinearity in the field equations of these theories implies that accurate theoretical predictions could only be made using high-resolution cosmological simulations. Previously, such simulations were usually done on regular meshes, which limits both their performance and the accuracy. In this paper, we report the development of a new algorithm and code, based on ECOSMOG, that uses adaptive mesh refinements to improve the efficiency and precision in simulating the models with Vainshtein mechanism. We have made various code tests against the numerical reliability, and found consistency with previous simulations. We also studied the velocity field in the self-accelerating branch of the DGP model. The code, parallelised using MPI, is suitable for large cosmological simulations of Galileon-type modified gravity theories.

Massive gravity is a good theoretical laboratory to study modifications of General Relativity. The theory offers a concrete set-up to study models of dark energy, since it admits cosmological self-accelerating solutions in the vacuum, in which the size of the acceleration depends on the graviton mass. Moreover, non-linear gravitational self-interactions, in the proximity of a matter source, manage to mimic the predictions of linearised General Relativity, hence agreeing with solar-system precision measurements. In this article, we review our work in the subject, classifying, on one hand, static solutions, and on the other hand, self-accelerating backgrounds. For what respects static solutions we exhibit black hole configurations, together with other solutions that recover General Relativity near a source via the Vainshtein mechanism. For the self-accelerating solutions we describe a wide class of cosmological backgrounds, including an analysis of their stability.

In massive gravity, galileon, and braneworld explanations of cosmic acceleration, force modifications are screened by nonlinear derivative self-interactions of the scalar field mediating that force. Interactions between the field of a central body ("A") and an orbiting body ("B") imply that body B does not move as a test body in the field of body A if the orbit is smaller than the Vainshtein radius of body B. We find through numerical solutions of the joint field at the position of B that the A-field Laplacian is nearly perfectly screened by the B self-field, whereas first derivative or net forces are reduced in a manner that scales with the mass ratio of the bodies as (M_B/M_A)^{3/5}. The latter causes mass-dependent reductions in the universal perihelion precession rate, with deviations for the Earth-Moon system at the ~4% level. In spite of universal coupling, which preserves the microscopic equivalence principle, the motion of macroscopic screened bodies depends on their mass providing in principle a means for testing the Vainshtein mechanism.

We calculate synchrotron radio emission and gamma-ray emission due to bremsstrahlung, inverse-Compton scattering and $pi^0$-decay from the remnant of supernova which exploded in the circumstellar matter (CSM) formed by the progenitor's stellar wind. This sort of situation is a possible origin of mixed-morphology supernova remnants (SNRs) like W49B, which exhibit recombination-radiation spectra in X-ray emission. We assume that the CSM of 1.5 $M_{odot}$ exists at 0.07--3 pc away from the supernova in the interstellar medium (ISM) of density 0.016 cm$^{-3}$. When the blast wave breaks out of the CSM into the ISM, its velocity rapidly increases and hence particle acceleration is enhanced. The maximum energy of protons reaches $sim$ 1300 TeV just after the break-out with $sim$ 0.5% of the explosion energy. We consider the non-thermal emission from the blast-shocked ISM shell after the break-out. Synchrotron radio flux at 1 GHz is tens Jy, comparable to the observed from mixed-morphology SNRs. Because of low density, the gamma-ray luminosity is dominated by inverse-Compton scattering, which is higher than the $pi^0$-decay luminosity by an order of magnitude. The total gamma-ray luminosity including bremsstrahlung is of the order of $10^{33}$ erg s$^{-1}$ lower than the typical value $10^{35}$--$10^{36}$ erg s$^{-1}$ observed from mixed-morphology SNRs. However, if, e.g., $sim$10% of accelerated protons interact with some matter of density $sim$100 cm$^{-3}$, $pi^0$-decay gamma-ray luminosity would be enhanced to be comparable with the observed value.

Using the unprecedented spectral resolution of the reflection grating spectrometer (RGS) onboard XMM-Newton, we reveal dynamics of X-ray--emitting ejecta in the oxygen-rich supernova remnant Puppis A. The RGS spectrum shows prominent K-shell lines, including O VII He_alpha forbidden and resonance, O VIII Ly_alpha, O VIII Ly_beta, and Ne IX He_alpha resonance, from an ejecta knot positionally coincident with an optical oxygen-rich filament (the so-called Omega filament in the northeast of the remnant). We find that the line centroids are blueshifted by 1480+/-140+/-60 km/s (the first and second term errors are measurement and calibration uncertainties, respectively), which is fully consistent with that of the optical Omega filament. Line broadening at 654 eV (corresponding to O VIII Ly_alpha) is obtained to be less than ~0.9 eV, indicating an oxygen temperature of less than ~30 keV. Analysis of XMM-Newton MOS spectra shows an electron temperature of ~0.8 keV and an ionization timescale of ~2e10 cm^-3 s. We show that the oxygen and electron temperatures as well as the ionization timescale can be reconciled if the ejecta knot was heated by a collisionless shock whose velocity is ~600-1200 km/s and was subsequently equilibrated due to Coulomb interactions. The RGS spectrum also shows relatively weak K-shell lines of another ejecta feature located near the northeastern edge of the remnant, from which we measure redward Doppler velocities of 650+/-70+/-60 km/s.

We develop a new, efficient code for solving the second-order Einstein-Boltzmann equations, and use it to estimate the intrinsic CMB non-Gaussianity arising from the non-linear evolution of density perturbations. The full calculation involves contributions from recombination and less tractable contributions from terms integrated along the line of sight. We investigate the bias that this intrinsic bispectrum implies for searches of primordial non-Gaussianity. We find that the inclusion or omission of certain line of sight terms can make a large impact. When including all physical effects but lensing and time-delay, we find that the local-type f_nl would be biased by f_nl ~ 0.5, below the expected sensitivity of the Planck satellite. The speed of our code allows us to confirm the robustness of our results with respect to a number of numerical parameters.

Forthcoming radio continuum surveys will cover large volumes of the observable Universe and will reach to high redshifts, making them potentially powerful probes of dark energy, modified gravity and non-Gaussianity. Here we extend recent works by analyzing the general relativistic (GR) corrections to the angular power spectrum. These GR corrections to the standard Newtonian analysis of the power spectrum become significant on scales near and beyond the Hubble scale at each redshift. We consider the continuum surveys with LOFAR, WSRT and ASKAP, and examples of continuum surveys with the SKA. We find that the GR corrections are at most percent-level in LOFAR, WODAN and EMU surveys, but they can produce $O(10%)$ changes for high enough sensitivity SKA continuum surveys. The signal is however dominated by cosmic variance, and multiple-tracer observations will be needed to overcome this problem. The GR corrections are suppressed in continuum surveys because of the integration over redshift -- we expect that GR corrections will be enhanced for future SKA HI surveys in which the source redshifts will be known. We also provide predictions for the angular power spectra in the case where the primordial perturbations have local non-Gaussianity. We find that non-Gaussianity dominates over GR corrections, and rises above cosmic variance when $f_{\rm NL}\gtrsim5$ for SKA continuum surveys.

We study the observational viability of coupled quintessence models with their expansion and growth histories matched to modified gravity cosmologies. We find that for a Dvali-Gabadadze-Porrati model which has been fitted to observations, the matched interacting dark energy models are observationally disfavoured. We also study the distinguishability of interacting dark energy models matched to scalar-tensor theory cosmologies and show that it is not always possible to find a physical interacting dark energy model which shares their expansion and growth histories.

We study multi-field DBI inflation models with a waterfall phase transition. This transition happens for a D3 brane moving in the warped conifold if there is an instability along angular directions. The transition converts the angular perturbations into the curvature perturbation. Thanks to this conversion, multi-field models can evade the stringent constraints that strongly disfavour single field ultra-violet DBI inflation models in string theory. We explicitly demonstrate that our model satisfies current observational constraints on the spectral index and equilateral non-Gaussianity as well as the bound on the tensor to scalar ratio imposed in string theory models. In addition we show that large local type non-Gaussianity is generated together with equilateral non-Gaussianity in this model.

We investigate the effect of modified gravity on the specific angular momentum of galactic halos by analyzing the halo catalogs at z=0 from high-resolution N-body simulations for a f(R) gravity model that meets the solar-system constraint. It is shown that the galactic halos in the f(R) gravity model tend to acquire significantly higher specific angular momentum than those in the standard LCDM model. The largest difference in the specific angular momentum distribution between these two models occurs for the case of the isolated galactic halos with mass less than 10^{11}M_sun/h, which are likely least shielded by the chameleon screening mechanism. As the specific angular momentum of galactic halos is rather insensitive to the other cosmological parameters, it can in principle be an independent discriminator of modified gravity. We speculate a possibility of using the relative abundance of the low surface brightness galaxies (LSBGs) as a test of GR given that the formation of the LSBGs befalls in the fast spinning dark halos.

This paper presents the Suzaku results obtained for the Sagittarius (Sgr) C region using the concept of X-ray reflection nebulae (XRNe) as the echo of past flares from the super massive black hole, Sgr A*. The Sgr C complex is composed of several molecular clouds proximately located in projected distance. The X-ray spectra of Sgr C were analyzed on the basis of a view that XRNe are located inside the Galactic center plasma X-ray emission with an oval distribution around Sgr A*. We found that the XRNe are largely separated in the line-of-sight position, and are associated with molecular clouds in different velocity ranges detected by radio observations. We also applied the same analysis to the Sgr B XRNe and completed a long-term light curve for Sgr A* occurring in the past. As a new finding, we determined that Sgr A* was experiencing periods of high luminosity already 500 years ago, which is longer than the previously reported value. Our results are consistent with a scenario that Sgr A* was continuously active with sporadic flux variabilities of Lx = 1-3 x 10^39 erg s^-1 in the past 50 to 500 years. The average past luminosity was approximately 4-6 orders of magnitude higher than that presently observed. In addition, two short-term flares of 5-10 years are found. Thus, the past X-ray flare should not be a single short-term flare, but can be interpreted as multiple flares superposed on a long-term high state.

The joint JAXA/NASA ASTRO-H mission is the sixth in a series of highly successful X-ray missions initiated by the Institute of Space and Astronautical Science (ISAS). ASTRO-H will investigate the physics of the high-energy universe via a suite of four instruments, covering a very wide energy range, from 0.3 keV to 600 keV. These instruments include a high-resolution, high-throughput spectrometer sensitive over 0.3-2 keV with high spectral resolution of Delta E < 7 eV, enabled by a micro-calorimeter array located in the focal plane of thin-foil X-ray optics; hard X-ray imaging spectrometers covering 5-80 keV, located in the focal plane of multilayer-coated, focusing hard X-ray mirrors; a wide-field imaging spectrometer sensitive over 0.4-12 keV, with an X-ray CCD camera in the focal plane of a soft X-ray telescope; and a non-focusing Compton-camera type soft gamma-ray detector, sensitive in the 40-600 keV band. The simultaneous broad bandpass, coupled with high spectral resolution, will enable the pursuit of a wide variety of important science themes.

We investigate vector contributions to the Lagrangian of $\Lambda_3-$massive gravity in the decoupling limit, the less explored sector of this theory. The main purpose is to understand the stability of maximally symmetric %self-accelerating vacuum solutions. Around self-accelerating configurations, vector degrees of freedom become strongly coupled since their kinetic terms vanish, so their dynamics is controlled by higher order interactions. Even in the decoupling limit, the vector Lagrangian contains an infinite number of terms. We develop a systematic method to covariantly determine the vector Lagrangian at each order in perturbations, fully manifesting the symmetries of the system. We show that, around self-accelerating solutions, the structure of higher order $p$-form Galileons arise, avoiding the emergence of a sixth BD ghost mode. However, a careful analysis shows that there are directions along which the Hamiltonian is unbounded from below. This instability can be interpreted as one of the available fifth physical modes behaving as a ghost. Therefore, we conclude that self-accelerating configurations, in the decoupling limit of $\Lambda_3$-massive gravity, are generically unstable.

We report the global distribution of the intensities of the K-shell lines from the He-like and H-like ions of S, Ar, Ca and Fe along the Galactic plane. From the profiles, we clearly separate the Galactic center X-ray emission (GCXE) and the Galactic ridge X-ray emission (GRXE). The intensity profiles of the He-like K$\alpha$ lines of S, Ar, Ca and Fe along the Galactic plane are approximately similar with each other, while not for the H-like Ly$\alpha$ lines. In particular, the profiles of H-like Ly$\alpha$ of S and Fe show remarkable contrast; a large excess of Fe and almost no excess of S lines in the GCXE compared to the GRXE. Although the prominent K-shell lines are represented by $\sim$1 keV and $\sim$7 keV temperature plasmas, these two temperatures are not equal between the GCXE and GRXE. In fact, the spectral analysis of the GCXE and GRXE revealed that the $\sim$1 keV plasma in the GCXE has lower temperature than that in the GRXE, and vice versa for the $\sim$7 keV plasma.

We study the matter and velocity divergence power spectra in a f(R) gravity theory and their time evolution measured from several large-volume N-body simulations with varying box sizes and resolution. We find that accurate prediction of the matter power spectrum in f(R) gravity places stronger requirements on the simulation than is the case with LCDM, because of the nonlinear nature of the fifth force. Linear perturbation theory is shown to be a poor approximation for the f(R) models, except when the chameleon effect is very weak. We show that the relative differences from the fiducial LCDM model are much more pronounced in the nonlinear tail of the velocity divergence power spectrum than in the matter power spectrum, which suggests that future surveys which target the collection of peculiar velocity data will open new opportunities to constrain modified gravity theories. A close investigation of the time evolution of the power spectra shows that there is a pattern in the evolution history, which can be explained by the properties of the chameleon-type fifth force in f(R) gravity. Varying the model parameter |f_R0|, which quantifies the strength of the departure from standard gravity, mainly varies the epoch marking the onset of the fifth force, as a result of which the different f(R) models are in different stages of the same evolutionary path at any given time

We present the Suzaku results of the supernova remnant (SNR) G346.6-0.2. The X-ray emission has a center-filled morphology with the size of 6' x 8' within the radio shell. Neither an ionization equilibrium nor non-equilibrium (ionizing) plasma can reproduce the spectra remaining shoulder-like residuals in the 2-4 keV band. These structures are possibly due to recombination of free electrons to the K-shell of He-like Si and S. The X-ray spectra are well fitted with a plasma model in a recombination dominant phase. We propose that the plasma was in nearly full ionized state at high temperature of 5 keV, then the plasma changed to a recombining phase due to selective cooling of electrons to lower temperature of 0.3 keV. G346.6-0.2 would be in an epoch of the recombining phase.

We report new features of the typical mixed-morphology (MM) supernova remnant (SNR) W44. In the X-ray spectra obtained with Suzaku, radiative recombination continua (RRCs) of highly ionized atoms are detected for the first time. The spectra are well reproduced by a thermal plasma in a recombining phase. The best-fit parameters suggest that the electron temperature of the shock-heated matters cooled down rapidly from $\sim1$,keV to $\sim 0.5$,keV, possibly due to adiabatic expansion (rarefaction) occurred $\sim20,000$ years ago. We also discover hard X-ray emission which shows an arc-like structure spatially-correlated with a radio continuum filament. The surface brightness distribution shows a clear anti-correlation with $^{12}$CO (J=2-1) emission from a molecular cloud observed with NANTEN2. While the hard X-ray is most likely due to a synchrotron enhancement in the vicinity of the cloud, no current model can quantitatively predict the observed flux.

We use large volume N-body simulations to predict the clustering of dark matter in redshift space in f(R) modified gravity cosmologies. This is the first time that the nonlinear matter and velocity fields have been resolved to such a high level of accuracy over a broad range of scales in this class of models. We find significant deviations from the clustering signal in standard gravity, with an enhanced boost in power on large scales and stronger damping on small scales in the f(R) models compared to GR at redshifts z<1. We measure the velocity divergence (P_\theta \theta) and matter (P_\delta \delta) power spectra and find a large deviation in the ratios \sqrt{P_\theta \theta/P_\delta \delta} and P_\delta \theta/P_\delta\delta, between the f(R) models and GR for 0.03<k/(h/Mpc)<0.5. In linear theory these ratios equal the growth rate of structure on large scales. Our results show that the simulated ratios agree with the growth rate for each cosmology (which is scale dependent in the case of modified gravity) only for extremely large scales, k<0.06h/Mpc at z=0. The velocity power spectrum is substantially different in the f(R) models compared to GR, suggesting that this observable is a sensitive probe of modified gravity. We demonstrate how to extract the matter and velocity power spectra from the 2D redshift space power spectrum, P(k,\mu), and can recover the nonlinear matter power spectrum to within a few percent for k<0.1h/Mpc. However, the model fails to describe the shape of the 2D power spectrum demonstrating that an improved model is necessary in order to reconstruct the velocity power spectrum accurately. The same model can match the monopole moment to within 3% for GR and 10% for the f(R) cosmology at k<0.2 h/Mpc at z=1. Our results suggest that the extraction of the velocity power spectrum from future galaxy surveys is a promising method to constrain deviations from GR.

We present non-linear weak lensing predictions for coupled dark energy models using the CoDECS simulations. We calculate the shear correlation function and error covariance expected for these models, for forthcoming ground-based (such as DES) and space-based (Euclid) weak lensing surveys. We obtain predictions for the discriminatory power of a ground-based survey similar to DES and a space-based survey such as Euclid in distinguishing between $\Lambda$CDM and coupled dark energy models; we show that using the non-linear lensing signal we could discriminate between $\Lambda$CDM and exponential constant coupling models with $\beta_0\geq0.1$ at $4\sigma$ confidence level with a DES-like survey, and $\beta_0\geq0.05$ at $5\sigma$ confidence level with Euclid. We also demonstrate that estimating the coupled dark energy models' non-linear power spectrum, using the $\Lambda$CDM Halofit fitting formula, results in biases in the shear correlation function that exceed the survey errors.

We test Einstein gravity using cosmological observations of both expansion and structure growth, including the latest data from supernovae (Union2.1), CMB (WMAP7), weak lensing (CFHTLS) and peculiar velocity of galaxies (WiggleZ). We fit modified gravity parameters of the generalized Poisson equations simultaneously with the effective equation of state for the background evolution, exploring the covariances and model dependence. The results show that general relativity is a good fit to the combined data. Using a Pad{\'e} approximant form for the gravity deviations accurately captures the time and scale dependence for theories like $f(R)$ and DGP gravity, and weights high and low redshift probes fairly. For current observations, cosmic growth and expansion can be fit simultaneously with little degradation in accuracy, while removing the possibility of bias from holding one aspect fixed.

Astrophysical tests of modified modified gravity theories in the nearby universe have been emphasized recently by Hui, Nicolis and Stubbs (2009) and Jain and VanderPlas (2011). A key element of such tests is the screening mechanism whereby general relativity is restored in massive halos or high density environments like the Milky Way. In chameleon theories of gravity, including all f(R) models, field dwarf galaxies may be unscreened and therefore feel an extra force, as opposed to screened galaxies. The first step to study differences between screened and unscreened galaxies is to create a 3D screening map. We use N-body simulations to test and calibrate simple approximations to determine the level of screening in galaxy catalogs. Sources of systematic errors in the screening map due to observational inaccuracies are modeled and their contamination is estimated. We then apply our methods to create a map out to 200 Mpc in the Sloan Digital Sky Survey footprint using data from the Sloan survey and other sources. In two companion papers this map will be used to carry out new tests of gravity using distance indicators and the disks of dwarf galaxies. We also make our screening map publicly available.

We expose a new symmetry for linear perturbations around a solution of non-linear Fierz-Pauli massive gravity plus a bare cosmological constant. The cosmological constant is chosen such that the background metric is flat while the Stuckelberg fields have a non-trivial profile. Around this background, at linear order the new symmetry reduces the propagating degrees of freedom to those of General Relativity, namely the massless helicity 2 modes only. We discuss the physical consequences and possible applications of these findings.

Local non-Gaussianity, parametrized by $f_{\rm NL}$, introduces a scale-dependent bias that is strongest at large scales, precisely where General Relativistic (GR) effects also become significant. With future data, it should be possible to constrain $f_{\rm NL} = {\cal O}(1)$ with high redshift surveys. GR corrections to the power spectrum and ambiguities in the gauge used to define bias introduce effects similar to $f_{\rm NL}= {\cal O}(1)$, so it is essential to disentangle these effects. For the first time in studies of primordial non-Gaussianity, we include the consistent GR calculation of galaxy power spectra, highlighting the importance of a proper definition of bias. We present observable power spectra with and without GR corrections, showing that an incorrect definition of bias can mimic non-Gaussianity. However, these effects can be distinguished by their different redshift and scale dependence, so as to extract the true primordial non-Gaussianity.

We present the Suzaku results of the mixed-morphology supernova remnant W28. The X-ray spectra of the central region of W28 exhibit many bright emission lines from highly ionized atoms. An optically thin thermal plasma in collisional ionization equilibrium, either of single-temperature or multi-temperature failed to reproduce the data with line-like and bump-like residuals at the Si Lyman$\alpha$ energy and at 2.4--5.0 keV, respectively. The bumps probably correspond to radiative recombination continua from He-like Si and S. A simple recombining plasma model nicely fit the bump structures, but failed to fit low energy bands. The overall spectra can be fit with a multi-ionization temperature plasma with a common electron temperature. The multi-ionization temperatures are interpreted as elemental difference of ionization and recombination timescales. These results prefer the rarefaction scenario for the origin of the recombining plasma.

The next generation of weak lensing surveys will trace the evolution of matter perturbations and gravitational potentials from the matter dominated epoch until today. Along with constraining the dynamics of dark energy, they will probe the relations between matter overdensities, local curvature, and the Newtonian potential. We work with two functions of time and scale to account for any modifications of these relations in the linear regime from those in the LCDM model. We perform a Principal Component Analysis (PCA) to find the eigenmodes and eigenvalues of these functions for surveys like DES and LSST. This paper builds on and significantly extends the PCA analysis of Zhao et al. (2009) in several ways. In particular, we consider the impact of some of the systematic effects expected in weak lensing surveys. We also present the PCA in terms of other choices of the two functions needed to parameterize modified growth on linear scales, and discuss their merits. We analyze the degeneracy between the modified growth functions and other cosmological parameters, paying special attention to the effective equation of state w(z). Finally, we demonstrate the utility of the PCA as an efficient data compression stage which enables one to easily derive constraints on parameters of specific models without recalculating Fisher matrices from scratch.

We report on results of imaging and spectral analysis of the supernova remnant (SNR) RCW 86 observed with Suzaku. The SNR is known to exhibit K-shell emission of low ionized Fe, possibly originating from supernova ejecta. We revealed the global distribution of the Fe-rich plasma in the entire remnant, for the first time; the Fe-K emission was clearly detected from the west, north, and south regions, in addition to the X-ray brighter shells of southwest and northeast, where the presence of the Fe-rich ejecta has already been reported. The spectrum of each region is well represented by a three-component model consisting of low- and high-temperature thermal plasmas and a non-thermal emission. The lower-temperature component, with elemental abundances of near the solar values, likely originates from the forward shocked interstellar medium, while the Fe-rich ejecta is described by the hotter plasma. From the morphologies of the forward and reverse shocks in the west region, the total ejecta mass is estimated to be 1-2M_sun for the typical explosion energy of ~ 1 x 10^{51} erg. The integrated flux of the Fe-K emission from the entire SNR roughly corresponds to a total Fe mass of about 1M_sun. Both of these estimates suggest a Type Ia supernova origin of this SNR. We also find possible evidence of an Fe-rich clump located beyond the forward-shock front in the north rim, which is reminiscent of ejecta knots observed in the Tycho and Vela SNRs.

In standard cosmologies, dark energy interacts only gravitationally with dark matter. There could be a non-gravitational interaction in the dark sector, leading to changes in the effective DE equation of state, in the redshift dependence of the DM density and in structure formation. We use CMB, BAO and SNIa data to constrain a model where the energy transfer in the dark sector is proportional to the DE density. There are two subclasses, defined by the vanishing of momentum transfer either in the DM or the DE frame. We conduct a Markov-Chain Monte-Carlo analysis to obtain best-fit parameters. The background evolution allows large interaction strengths, and the constraints from CMB anisotropies are weak. The growth of DM density perturbations is much more sensitive to the interaction, and can deviate strongly from the standard case. However, the deviations are degenerate with galaxy bias and thus more difficult to constrain. Interestingly, the ISW signature is suppressed since the non-standard background evolution can compensate for high growth rates. We also discuss the partial degeneracy between interacting DE and modified gravity, and how this can be broken.

We study a cosmological scenario in which the DBI action governing the motion of a D3-brane in a higher-dimensional spacetime is supplemented with an induced gravity term. The latter reduces to the quartic Galileon Lagrangian when the motion of the brane is non-relativistic and we show that it tends to violate the null energy condition and to render cosmological fluctuations ghosts. There nonetheless exists an interesting parameter space in which a stable phase of quasi-exponential expansion can be achieved while the induced gravity leaves non trivial imprints. We derive the exact second-order action governing the dynamics of linear perturbations and we show that it can be simply understood through a bimetric perspective. In the relativistic regime, we also calculate the dominant contribution to the primordial bispectrum and demonstrate that large non-Gaussianities of orthogonal shape can be generated, for the first time in a concrete model. More generally, we find that the sign and the shape of the bispectrum offer powerful diagnostics of the precise strength of the induced gravity.

In this paper, we study the distribution of dark matter halos and voids using high resolution simulations in f(R) gravity models with the chameleon mechanism to screen the fifth force in dense environment. For dark matter halos, we show that the semi-analytic thin shell condition, with a suitably-defined environment, provides a good approximation to describe the mass and environmental dependence of the screening of the fifth force in halos. Due to stronger gravity, there are far more massive halos and large voids in f(R) models compared with the \Lambda CDM model. The numbers of voids with an effective radius of 15Mpc/h are twice and four times as many as those in \Lambda CDM for f(R) models with |f_{R0}|=1e-5 and 1e-4 respectively. This provides a new means to test the models using the upcoming observational data. We also find that halos inside voids are all unscreened in our simulations, which are ideal objects for the gravity test.

We investigate strong coupling effects in a covariant massive gravity model, which is a candidate for a ghost free non-linear completion of Fierz-Pauli. We analyse the conditions to recover general relativity via Vainshtein mechanism in the weak field limit, and find three main cases depending on the choice of parameters. In the first case, the potential is such that all non-linearities disappear and the vDVZ discontinuity cannot be avoided. In the second case, the Vainshtein mechanism allows to recover general relativity within a macroscopic radius from a source. In the last case, the strong coupling of the scalar graviton completely shields the massless graviton, and weakens gravity when approaching the source. In the second part of the paper, we explore new exact vacuum solutions, that asymptote de Sitter or anti de Sitter space depending on the parameter choice. The curvature of the space is proportional to the mass of the graviton, thus providing a cosmological background which may explain the present day acceleration in terms of the graviton mass. Moreover, by expressing the potential for non-linear massive gravity in a convenient form, we also suggest possible connections with a higher dimensional framework.

We carry out 3D-hydrodynamical calculations for the interaction of expanding supernova ejecta with the dense circumstellar matter (CSM) and the rarefied interstellar medium (ISM) outside. The CSM is composed of the stellar wind matter from the progenitor in its pre-supernova phase, and assumed to be axially symmetric: more matter around the equator than in the polar direction driven by rotation of the progenitor. Because of high density of the CSM, the ionization state of the shock-heated ejecta quickly becomes equilibrium with the electron temperature. When the blast wave breaks out of the CSM into the rarefied ISM, the shocked ejecta cools rapidly due to adiabatic expansion, and hence an over-ionized/recombining plasma would be left. The ejecta is reheated by the second reverse shock due to the interaction with the ISM. We calculate the emission measure of the supernova remnant (SNR) along the line of sight, and find that the over-ionized plasma appears to be bar-like with wings in the edge-on (equatorial view), while shell-like in the face-on (polar view) geometry with respect to the rotation axis. The hot gas heated by the blast wave exists in the outermost region of the SNR with a nearly complete shell, but the X-rays therefrom are too faint to be observable. Thus, depending on the viewing angle, the SNR of the over-ionized plasma would exhibit center-filled morphology in X-rays, like W49B, a mixed-morphology SNR. The bar-like structure is swept out by the second reverse shock and disappears eventually, and then the SNR becomes shell-like in both the equatorial and polar views in the later phase of the evolution.

We explore the possibility of realising self-accelerated expansion of the Universe taking into account the vector components of a massive graviton. The effective action in the decoupling limit contains an infinite number of terms, once the vector degrees of freedom are included. These can be re-summed in physically interesting situations, which result in non-polynomial couplings between the scalar and vector modes. We show there are self-accelerating background solutions for this effective action, with the possibility of having a non-trivial profile for the vector fields. We then study fluctuations around these solutions and show that there is always a ghost, if a background vector field is present. When the background vector field is switched off, the ghost can be avoided, at the price of entering into a strong coupling regime, in which the vector fluctuations have vanishing kinetic terms. Finally we show that the inclusion of a bare cosmological constant does not change the previous conclusions and it does not lead to a ghost mode in the absence of a background vector field.

We introduce a new code, ECOSMOG, to run N-body simulations for a wide class of modified gravity and dynamical dark energy theories. These theories generally have one or more new dynamical degrees of freedom, the dynamics of which are governed by their (usually rather nonlinear) equations of motion. Solving these non-linear equations has been a great challenge in cosmology. Our code is based on the RAMSES code, which solves the Poisson equation on adaptively refined meshes to gain high resolutions in the high-density regions. We have added a solver for the extra degree(s) of freedom and performed numerous tests for the f(R) gravity model as an example to show its reliability. We find that much higher efficiency could be achieved compared with other existing mesh/grid-based codes thanks to two new features of the present code: (1) the efficient parallelisation and (2) the usage of the multigrid relaxation to solve the extra equation(s) on both the regular domain grid and refinements, giving much faster convergence even under much more stringent convergence criteria. This code is designed for performing high-accuracy, high-resolution and large-volume cosmological simulations for modified gravity and general dark energy theories, which can be utilised to test gravity and the dark energy hypothesis using the upcoming and future deep and high-resolution galaxy surveys.

We analyze the diffuse Fe I K-alpha line generated in the diffuse interstellar molecular hydrogen by primary photons or subrelativistic protons injected by Sagittarius (Sgr) A*. We showed that unlike emission from compact molecular clouds, this emission can be permanently observed in the directions of the Galactic center. We conclude that the diffuse emission of 6.4 keV line observed at present is probably due to Fe I K-alpha vacancy production by primary photons if the X-ray luminosity of Sgr A* was about Lx ~ 10^39-10^40 erg/s. In principle these data can also be described in the framework of the model when the 6.4 keV line emission is generated by subrelativistic protons generated by accretion onto the central black hole but in this case extreme parameters of this model are necessary.

Prominent K-shell emission lines of neutral iron (hereafter, FeI-K) and hard-continuum X-rays were found from molecular clouds (MCs) in the Sagittarius B (Sgr B) region with the two separate Suzaku observations in 2005 and 2009. The X-ray flux of FeI-K decreased in correlation to the hard-continuum flux by factor of 0.4-0.5 in 4 years, which is nearly equal to the light-travelling across the MCs. The rapid and correlated time-variability, the equivalent width of FeI-K, and the K-edge absorption depth of FeI are consistently explained by "X-ray echoes" due to the fluorescent and Thomson-scattering of an X-ray flare from an external source. The required flux of the X-ray flare depends on the distance to the MCs and the duration time. The flux, even in the minimum case, is larger than those of the brightest Galactic X-ray sources. Based on these facts, we conclude that the super-massive black hole, Sgr A*, exhibited a big-flare about a few hundred years ago and the luminosity of higher than 4x10^39 erg s^{-1}. The "X-ray echo" from Sgr B, located at a few hundred light-years from Sgr A*, now arrived at the Earth.

We study spherically symmetric solutions in a covariant massive gravity model, which is a candidate for a ghost-free non-linear completion of the Fierz-Pauli theory. There is a branch of solutions that exhibits the Vainshtein mechanism, recovering General Relativity below a Vainshtein radius given by $(r_g m^2)^{1/3}$, where $m$ is the graviton mass and $r_g$ is the Schwarzschild radius of a matter source. Another branch of exact solutions exists, corresponding to Schwarzschild-de Sitter spacetimes where the curvature scale of de Sitter space is proportional to the mass squared of the graviton.

In this Letter, we investigate the environmental dependence of dark matter halos in theories that attempt to explain the accelerated expansion of the Universe by modifying general relativity (GR). Using high-resolution N-body simulations in f(R) gravity models which recover GR in dense environments by virtue of the chameleon mechanism, we find a strong environmentally-dependent difference between the lensing mass and dynamical mass estimates of dark matter halos. This environmental dependence of the halo properties can be used as a smoking gun to test GR observationally.

We have surveyed spatial profiles of the Fe K$\alpha$ lines in the Galactic center diffuse X-rays (GCDX), including the transient region from the GCDX to the Galactic ridge X-ray emission (GRXE), with the Suzaku satellite. We resolved Fe K$\alpha$ line complex into three lines of Fe \emissiontype{I}, Fe \emissiontype{XXV} and Fe \emissiontype{XXVI} K$\alpha$, and obtained their spatial intensity profiles with the resolution of $\sim \timeform{0D.1}$. We compared the Fe \emissiontype{XXV} K$\alpha$ profile with a stellar mass distribution (SMD) model made from near infrared observations. The intensity profile of Fe \emissiontype{XXV} K$\alpha$ is nicely fitted with the SMD model in the GRXE region, while that in the GCDX region shows $3.8\pm0.3$ $(\timeform{0D.2}<|l|<\timeform{1D.5})$ or $19\pm6$ $(|l|<\timeform{0D.2})$ times excess over the best-fit SMD model in the GRXE region. Thus Fe \emissiontype{XXV} K$\alpha$ in the GCDX is hardly explained by the same origin of the GRXE. In the case of point source origin, a new population with the extremely strong Fe \emissiontype{XXV} K$\alpha$ line is required. An alternative possibility is that the majority of the GCDX is truly diffuse optically thin thermal plasma.

The Tornado nebula (G357.7-0.1) is a mysterious radio source with bright "head" and faint "tail" located in the direction of the Galactic center (GC) region. We here report the discovery of two diffuse X-ray sources at the head and tail of the Tornado with the Suzaku satellite. We found emission lines from highly ionized atoms in the two sources. The spectra are reproduced by an optically thin thermal plasma with a common temperature of 0.6-0.7 keV. The interstellar absorption (NH) of these sources are the same and are slightly larger than that of the GC distance. Since the estimated distance using the NH value is consistent with the radio observation of the Tornado, these X-ray sources are likely associated with the Tornado nebula. The twin-plasma morphology at the both ends of the Tornado suggests that the system is a bipolar/outflow source.

With the Suzaku satellite, we observed an unidentified TeV gamma-ray source HESS J1741$-$302 and its surroundings. No diffuse or point-like X-ray sources are detected from the bright southern emission peak of HESS J1741$-$302. From its neighborhood, we found a new intermediate polar candidate at the position of $(\alpha, \delta)_{\rm J2000.0} = (\timeform{17h40m35.6s}, \timeform{-30D14m16s})$, which is designated as Suzaku J174035.6$-$301416. The spectrum of Suzaku J174035.6$-$301416 exhibits emission lines at the energy of 6.4, 6.7 and 7.0 keV, which can be assigned as the K$\alpha$ lines from neutral, He-like and H-like iron, respectively. A coherent pulsation is found at a period of 432.1 $\pm$ 0.1 s. The pulse profile is quasi-sinusoidal in the hard X-ray band (4$-$8 keV), but is more complicated in the soft X-ray band (1$-$3 keV). The moderate period of pulsation, the energy flux, and the presence of the iron K$\alpha$ lines indicate that Suzaku J174035.6$-$301416 is likely an intermediate polar, a subclass of magnetized white dwarf binaries (cataclysmic variables). Based on these discoveries, we give some implications on the origin of GCDX and brief comments on HESS J1741$-$302 and PSR B1737$-$30.

The next generation of weak lensing surveys will trace the growth of large scale perturbations through a sequence of epochs, offering an opportunity to test General Relativity (GR) on cosmological scales. We review in detail the parametrization used in MGCAMB to describe the modified growth expected in alternative theories of gravity and generalized dark energy models. We highlight its advantages and examine several theoretical aspects. In particular, we show that the same set of equations can be consistently used on super-horizon and sub-horizon linear scales. We also emphasize the sensitivity of data to scale-dependent features in the growth pattern, and propose using Principal Component Analysis to converge on a practical set of parameters which is most likely to detect departures from GR. The connection with other parametrizations is also discussed.

We investigated the emission of K$\alpha$ iron line from the massive molecular clouds in the Galactic center (GC). We assume that at present the total flux of this emission consists of time variable component generated by primary X-ray photons ejected by Sagittarius A$^\ast$ (Sgr A$^\ast$) in the past and a relatively weak quasi-stationary component excited by impact of protons which were generated by star accretion onto the central black hole. The level of background emission was estimated from a rise of the 6.4 keV line intensity in the direction of several molecular clouds, that we interpreted as a stage when the X-ray front ejected by Sgr A$^\ast$ entered into these clouds. The 6.4 keV emission before this intensity jump we interpreted as emission generated by subrelativistic cosmic rays there. The cross-section of K$\alpha$ vacancies produced by protons differs from that of electrons or X-rays. Therefore, we expect that this processes can be distinguished from the analysis of the equivalent width of the iron line and time variations of the width can be predicted. The line intensity from the clouds depends on their distance from Sgr A$^\ast$ and the coefficient of spacial diffusion near the Galactic center. We expect that in a few years the line intensity for the cloud G\,0.11$-$0.11 which is relatively close to Sgr A$^\ast$ will decreases to the level $\lesssim$ 10% from its present value. For the cloud Sagittarius B2 (Sgr B2) the situation is more intricate. If the diffusion coefficient $D\gtrsim 10^{27}$ cm$^2$ s$^{-1}$ then the expected stationary flux should be about 10% of its level in 2000. In the opposite case the line intensity from Sgr B2 should drop down to zero because the protons do not reach the cloud.

We present the Suzaku results of a supernova remnant (SNR), G359.1-0.5 in the direction of the Galactic center region. From the SNR, we find prominent K-shell lines of highly ionized Si and S ions, together with unusual structures at 2.5-3.0 and 3.1-3.6 keV. No canonical SNR plasma model, in either ionization equilibrium or under-ionization, can explain the structures. The energies and shapes of the structures are similar to those of the radiative transitions of free electrons to the K-shell of He-like Si and S ions (radiative recombination continuum: RRC). The presence of the strong RRC structures indicates that the plasma is in over-ionization. In fact, the observed spectrum is well fitted with an over-ionized plasma model. The best-fit electron temperature of 0.29 keV is far smaller than the ionization temperature of 0.77 keV, which means that G359.1-0.5 is in extreme condition of over-ionization. We report some cautions on the physical parameters, and comment possible origins for the over-ionized plasma.

We perform high resolution N-body simulations for f(R) gravity based on a self-adaptive particle- mesh code MLAPM. The Chameleon mechanism that recovers General Relativity on small scales is fully taken into account by self-consistently solving the non-linear equation for the scalar field. We independently confirm the previous simulation results, including the matter power spectrum, halo mass function and density profiles, obtained by Oyaizu et al. (Phys.Rev.D 78, 123524, 2008) and Schmidt et al. (Phys.Rev.D 79, 083518, 2009), and extend the resolution up to k~20 h/Mpc for the measurement of the matter power spectrum. Based on our simulation results, we discuss how the Chameleon mechanism affects the clustering of dark matter and halos on full non-linear scales.

We explore the complementarity of weak lensing and galaxy peculiar velocity measurements to better constrain modifications to General Relativity. We find no evidence for deviations from GR on cosmological scales from a combination of peculiar velocity measurements (for Luminous Red Galaxies in the Sloan Digital Sky Survey) with weak lensing measurements (from the CFHT Legacy Survey). We provide a Fisher error forecast for a Euclid-like space-based survey including both lensing and peculiar velocity measurements, and show that the expected constraints on modified gravity will be at least an order of magnitude better than with present data, i.e. we will obtain 5% errors on the modified gravity parametrization described here. We also present a model--independent method for constraining modified gravity parameters using tomographic peculiar velocity information, and apply this methodology to the present dataset.

We study the inflationary dynamics in a model of slow-roll inflation in warped throat. Inflation is realized by the motion of a D-brane along the radial direction of the throat, and at later stages instabilities develop in the angular directions. We closely investigate both the single field potential relevant for the slow-roll phase, and the full multi-field one including the angular modes which becomes important at later stages. We study the main features of the instability process, discussing its possible consequences and identifying the vacua towards which the angular modes are driven.

The joint JAXA/NASA ASTRO-H mission is the sixth in a series of highly successful X-ray missions initiated by the Institute of Space and Astronautical Science (ISAS). ASTRO-H will investigate the physics of the high-energy universe by performing high-resolution, high-throughput spectroscopy with moderate angular resolution. ASTRO-H covers very wide energy range from 0.3 keV to 600 keV. ASTRO-H allows a combination of wide band X-ray spectroscopy (5-80 keV) provided by multilayer coating, focusing hard X-ray mirrors and hard X-ray imaging detectors, and high energy-resolution soft X-ray spectroscopy (0.3-12 keV) provided by thin-foil X-ray optics and a micro-calorimeter array. The mission will also carry an X-ray CCD camera as a focal plane detector for a soft X-ray telescope (0.4-12 keV) and a non-focusing soft gamma-ray detector (40-600 keV) . The micro-calorimeter system is developed by an international collaboration led by ISAS/JAXA and NASA. The simultaneous broad bandpass, coupled with high spectral resolution of Delta E ~7 eV provided by the micro-calorimeter will enable a wide variety of important science themes to be pursued.

We study primordial fluctuations generated during inflation in a class of models motivated by the DBI Galileons, which are extensions of the DBI action that yield second order field equations. This class of models generalises the DBI Galileons in a similar way with K-inflation. We calculate the primordial non-Gaussianity from the bispectrum of the curvature perturbations at leading order in the slow-varying approximations. We show that the estimator for the equilateral-type non-Gaussianity, $f_{\rm NL} ^{equil}$, can be applied to measure the amplitude of the primordial bispectrum even in the presence of the Galileon-like term although it gives a slightly different momentum dependence from K-inflation models. For the DBI Galileons, we find $-0.32 /c_s^2 < f_{\rm NL} ^{equil} < -0.16/c_s^2$ and large primordial non-Gaussianities can be obtained when $c_s$ is much smaller than 1 as in the usual DBI inflation. In G-inflation models, where a de Sitter solution is obtained without any potentials, the non-linear parameter is given by $f_{\rm NL}^{equil} = 4.62 r^{-2/3}$ where $r$ is the tensor to scalar ratio, giving a stringent constraint on the model.

The observable universe could be a 1+3-surface (the "brane") embedded in a 1+3+\textit{d}-dimensional spacetime (the "bulk"), with Standard Model particles and fields trapped on the brane while gravity is free to access the bulk. At least one of the \textit{d} extra spatial dimensions could be very large relative to the Planck scale, which lowers the fundamental gravity scale, possibly even down to the electroweak ($\sim$ TeV) level. This revolutionary picture arises in the framework of recent developments in M theory. The 1+10-dimensional M theory encompasses the known 1+9-dimensional superstring theories, and is widely considered to be a promising potential route to quantum gravity. At low energies, gravity is localized at the brane and general relativity is recovered, but at high energies gravity "leaks" into the bulk, behaving in a truly higher-dimensional way. This introduces significant changes to gravitational dynamics and perturbations, with interesting and potentially testable implications for high-energy astrophysics, black holes, and cosmology. Brane-world models offer a phenomenological way to test some of the novel predictions and corrections to general relativity that are implied by M theory. This review analyzes the geometry, dynamics and perturbations of simple brane-world models for cosmology and astrophysics, mainly focusing on warped 5-dimensional brane-worlds based on the Randall--Sundrum models. We also cover the simplest brane-world models in which 4-dimensional gravity on the brane is modified at \emph{low} energies -- the 5-dimensional Dvali--Gabadadze--Porrati models. Then we discuss co-dimension two branes in 6-dimensional models.

We investigate an estimator to measure the primordial trispectrum in equilateral type non-Gaussian models such as k-inflation, single field DBI inflation and multi-field DBI inflation models from Cosmic Microwave Background (CMB) anisotropies. The shape of the trispectrum whose amplitude is not constrained by the bispectrum in the context of effective theory of inflation and k-inflation is known to admit a separable form of the estimator for CMB anisotropies. We show that this shape is $87 \%$ correlated with the full quantum trispectrum in single field DBI inflation, while it is $33 \%$ correlated with the one in multi-field DBI inflation when curvature perturbation is originated from purely entropic contribution. This suggests that $g_{\rm NL} ^{equil}$, the amplitude of this particular shape, provides a reasonable measure of the non-Gaussianity from the trispectrum in equilateral non-Gaussian models. We relate model parameters such as the sound speed, $c_s$ and the transfer coefficient from entropy perturbations to the curvature perturbation, $T_{\mathcal{R} S}$ with $g_{\rm NL} ^{equil}$, which enables us to constrain model parameters in these models once $g_{\rm NL}^{equil}$ is measured in WMAP and Planck.

Diffuse X-rays from the Galactic center (GC) region were found to exhibit many K-shell lines from iron and nickel atoms in the 6--9 keV band. The strong emission lines seen in the spectrum are neutral iron K$\alpha$ at 6.4~keV, He-like iron K$\alpha$ at 6.7~keV, H-like iron Ly$\alpha$ at 6.9~keV, and He-like iron K$\beta$ at 7.8~keV. Among them, the 6.4~keV emission line is a probe of non-thermal phenomena. We have detected strong 6.4~keV emission in several giant molecular clouds, some of which were newly discovered by Suzaku. All the spectra exhibit large equivalent widths of 1-2~keV and absorption columns of $2-10\times 10^{23}\{\rm H\ cm}^{-2}$. We found time variability of diffuse 6.4~keV emission in the Sgr B2 region comparing the maps and spectra obtained from 1994 to 2005 with ASCA, Chandra, XMM-Newton and Suzaku. We also report discovery of K$\alpha$ lines of neutral argon, calcium, chrome, and manganese atoms in the Sgr~A region. We show that the equivalent width of the 6.4~keV emission line detected in X-ray faint region against the 6.4 keV-associated continuum (power-law component) is $\sim 800\ {\rm eV}$. These features are naturally explained by the X-ray reflection nebula scenario rather than the low energy cosmic-ray electrons scenario. On the other hand, a 6.4~keV clump, G~0.162$-$0.217, discovered at the south end of the Radio Arc has a small equivalent width of 6.4~keV emission line of $\sim200\ {\rm eV}$. The Radio Arc is a site of relativistic electrons. Thus, it is conceivable that the X-rays of G~0.162$-$0.217 are due to low energy cosmic-ray electrons

We study the halo bispectrum from non-Gaussian initial conditions. Based on a set of large $N$-body simulations starting from initial density fields with local type non-Gaussianity, we find that the halo bispectrum exhibits a strong dependence on the shape and scale of Fourier space triangles near squeezed configurations at large scales. The amplitude of the halo bispectrum roughly scales as $f_nl^2$. The resultant scaling on the triangular shape is consistent with that predicted by Jeong & Komatsu based on perturbation theory. We systematically investigate this dependence with varying redshifts and halo mass thresholds. It is shown that the $f_nl$ dependence of the halo bispectrum is stronger for more massive haloes at higher redshifts. This feature can be a useful discriminator of inflation scenarios in future deep and wide galaxy redshift surveys.

We test General Relativity (GR) using current cosmological data: the cosmic microwave background (CMB) from WMAP5 (Komatsu et al. 2009), the integrated Sachs-Wolfe (ISW) effect from the cross-correlation of the CMB with six galaxy catalogs (Giannantonio et al. 2008), a compilation of supernovae Type Ia (SNe) including the latest SDSS SNe (Kessler et al. 2009), and part of the weak lensing (WL) data from CFHTLS (Fu et al. 2008, Kilbinger et al. 2009) that probe linear and mildly non-linear scales. We first test a model where the effective Newton's constant, mu, and the ratio of the two gravitational potentials, eta, transit from the GR value to another constant at late times; in this case, we find that standard GR is fully consistent with the combined data. The strongest constraint comes from the ISW effect which would arise from this gravitational transition; the observed ISW signal imposes a tight constraint on a combination of mu and eta that characterizes the lensing potential. Next, we consider four pixels in time and space for each function mu and eta, and perform a Principal Component Analysis (PCA) finding that seven of the resulting eight eigenmodes are consistent with GR within the errors. Only one eigenmode shows a 2-sigma deviation from the GR prediction, which is likely to be due to a systematic effect. However, the detection of such a deviation demonstrates the power of our time- and scale-dependent PCA methodology when combining observations of structure formation and expansion history to test GR.

We report the Suzaku observation of 1E 1740.7-2942, a black hole candidate called the "Great Annihilator" (GA). The high-quality spectrum of Suzaku provides the severest constraints on the parameters of the GA. Two clumpy structures are found around the GA in the line images of FeI Kalpha at 6.4 keV and SXV Kalpha at 2.45 keV. One clump named M359.23-0.04 exhibits the 6.4-keV line with an equivalent width of ~ 1.2 keV, and is associated with a molecular cloud in the radio CS(J=1-0) map. Thus the 6.4-keV line from M359.23-0.04 is likely due to X-ray fluorescence irradiated by an external X-ray source. The irradiating X-rays would be either the past flare of Sagittarius A* or the bright nearby source, the GA. The other clump named G359.12-0.05 is associated with the radio supernova remnant candidate G359.07-0.02. We therefore propose that G359.12-0.05 is an X-ray counterpart of G359.07-0.02. G359.12-0.05 has a thin thermal plasma spectrum with a temperature of kT ~ 0.9 keV. The plasma parameters of G359.12-0.05 are consistent with those of a single supernova remnant in the Galactic center region.

The K-shell emission line of neutral irons from the Galactic center (GC) region is one of the key for the structure and activity of the GC. The origin is still open question, but possibly due either to X-ray radiation or to electron bombarding to neutral atoms. To address this issue, we analyzed the Suzaku X-ray spectrum from the GC region of intense neutral iron line emission, and report on the discovery of Kalpha lines of neutral argon, calcium, chrome, and manganese atoms. The equivalent widths of these Kalpha lines indicate that the metal abundances in the GC region should be ~1.6 and ~4 of solar value, depending on the X-ray and the electron origins, respectively. On the other hand, the metal abundances in the hot plasma in the GC region are found to be ~1-2 solar. These results favor that the origin of the neutral Kalpha lines are due to X-ray irradiation.

In this review, we discuss how non-Gaussianity of cosmological perturbations arises from inflation. After introducing the in-in formalism to calculate the $n$-point correlation function of quantum fields, we present the computation of the bispectrum of the curvature perturbation generated in general single field inflation models. The shapes of the bispectrum are compared with the local-type non-Gaussianity that arises from non-linear dynamics on super-horizon scales.

We confirm the recent claims that, in the infrared limit of Ho\v{r}ava-Lifshitz gravity, the scalar graviton becomes a ghost if the sound speed squared is positive on the flat de Sitter and Minkowski background. In order to avoid the ghost and tame the instability, the sound speed squared should be negative and very small, which means that the flow parameter $\lambda$ should be very close to its General Relativity (GR) value. We calculate the cubic interactions for the scalar graviton which are shown to have a similar structure with those of the curvature perturbation in k-inflation models. The higher order interactions become increasing important for a smaller sound speed squared, that is, when the theory approaches GR. This invalidates any linearized analysis and any predictability is lost in this limit as quantum corrections are not controllable. This pathological behaviour of the scalar graviton casts doubt on the validity of the projectable version of the theory.

Next generation surveys will observe the large-scale structure of the Universe with unprecedented accuracy. This will enable us to test the relationships between matter over-densities, the curvature perturbation and the Newtonian potential. Any large-distance modification of gravity or exotic nature of dark energy modifies these relationships as compared to those predicted in the standard smooth dark energy model based on General Relativity. In linear theory of structure growth such modifications are often parameterised by virtue of two functions of space and time that enter the relation of the curvature perturbation to, first, the matter over-density, and second, the Newtonian potential. We investigate the predictions for these functions in Brans-Dicke theory, clustering dark energy models and interacting dark energy models. We find that each theory has a distinct path in the parameter space of modified growth. Understanding these theoretical priors on the parameterisations of modified growth is essential to reveal the nature of cosmic acceleration with the help of upcoming observations of structure formations.

We compute the second-order matching conditions for tensor metric perturbations at an abrupt change in the equation of state. For adiabatic perturbations on large scales the matching hypersurface coincides with a uniform-density hypersurface. We show that in the uniform-density gauge both the tensor perturbation and its time-derivative are continuous in this case. For non-adiabatic perturbations, the matching hypersurface need not coincide with a uniform-density hypersurface and the tensor perturbation in the uniform-density gauge may be discontinuous. However, we show that in the Poisson gauge both the tensor perturbation and its time-derivative are continuous for adiabatic or non-adiabatic perturbations. As an application we solve the evolution equation for second-order tensor perturbations on large scales for a constant equation of state and we use the matching conditions to evolve the solutions through the transition from an inflationary era to a radiation era. We show that in the radiation era the resulting free part of the large-scale tensor perturbation (constant mode) is slow-roll suppressed in both the uniform-density and Poisson gauges. Thus, we conclude that second-order gravitational waves from slow-roll inflation are suppressed.

We present a set of predictions for weak lensing correlation functions in the context of modified gravity models, including a prescription for the impact of the nonlinear power spectrum regime in these models. We consider the DGP and f(R) models, together with dark energy models with the same expansion history. We use the requirement that gravity is close to GR on small scales to estimate the non-linear power for these models. We then calculate weak lensing statistics, showing their behaviour as a function of scale and redshift, and present predictions for measurement accuracy with future lensing surveys, taking into account cosmic variance and galaxy shape noise. We demonstrate the improved discriminatory power of weak lensing for testing modified gravities once the nonlinear power spectrum contribution has been included. We also examine the ability of future lensing surveys to constrain a parameterisation of the non-linear power spectrum, including sensitivity to the growth factor.

We present a cosmological model with a solution that self-accelerates at late-times without signs of ghost instabilities on small scales. The model is a natural extension of the Brans-Dicke (BD) theory including a non-linear derivative interaction, which appears in a theory with the Galilean shift symmetry. The existence of the self-accelerating universe requires a negative BD parameter but, thanks to the non-linear term, small fluctuations around the solution are stable on small scales. General relativity is recovered at early times and on small scales by this non-linear interaction via the Vainshtein mechanism. At late time, gravity is strongly modified and the background cosmology shows a phantom-like behaviour and the growth rate of structure formation is enhanced. Thus this model leaves distinct signatures in cosmological observations and it can be distinguished from standard $\Lambda$CDM cosmology.

Dark energy in General Relativity is typically non-interacting with other matter. However, it is possible that the dark energy interacts with the dark matter, and in this case, the dark matter can violate the universality of free fall (the weak equivalence principle). We show that some forms of the dark sector interaction do not violate weak equivalence. For those interactions that do violate weak equivalence, there are no available laboratory experiments to probe this violation for dark matter. But cosmology provides a test for violations of the equivalence principle between dark matter and baryons -- via a test for consistency of the observed galaxy velocities with the Euler equation.

We compute the leading order connected four-point function of the primordial curvature perturbation coming from the four-point function of the fields in multi-field DBI inflation models. We confirm that the consistency relations in the squeezed limit and in the counter-collinear limit hold as in single field models thanks to special properties of the DBI action. We also study the momentum dependence of the trispectra coming from the adiabatic, mixed and purely entropic contributions separately and we find that they have different momentum dependence. This means that if the amount of the transfer from the entropy perturbations to the curvature perturbation is significantly large, the trispectrum can distinguish multi-field DBI inflation models from single field DBI inflation models. A large amount of transfer $T_{\mathcal{RS}} \gg 1 $ suppresses the tensor to scalar ratio $r \propto T_{\mathcal{RS}}^{-2}$ and the amplitude of the bispectrum $f_{NL}^{equi} \propto T_{\mathcal{RS}}^{-2}$ and so it can ease the severe observational constraints on the DBI inflation model based on string theory. On the other hand, it enhances the amplitude of the trispectrum $\tau_{NL}^{equi} \propto T_{\mathcal{RS}}^2 f_{NL}^{equi 2}$ for a given amplitude of the bispectrum.

We present a hard X-ray spectrum of unprecedented quality of the Galactic supernova remnant W49B obtained with the Suzaku satellite. The spectrum exhibits an unusual structure consisting of a saw-edged bump above 8 keV. This bump cannot be explained by any combination of high-temperature plasmas in ionization equilibrium. We firmly conclude that this bump is caused by the strong radiative recombination continuum (RRC) of iron, detected for the first time in a supernova remnant. The electron temperature derived from the bremsstrahlung continuum shape and the slope of the RRC is 1.5 keV. On the other hand, the ionization temperature derived from the observed intensity ratios between the RRC and K-alpha lines of iron is 2.7 keV. These results indicate that the plasma is in a highly overionized state. Volume emission measures independently determined from the fluxes of the thermal and RRC components are consistent with each other, suggesting the same origin of these components.

We present the Suzaku spectroscopic study of the Galactic middle-aged supernova remnant (SNR) IC 443. The X-ray spectrum in the 1.75-6.0 keV band is described by an optically-thin thermal plasma with the electron temperature of 0.6 keV and several additional Lyman lines. We robustly detect, for the first time, strong radiative recombination continua (RRC) of H-like Si and S around at 2.7 and 3.5 keV. The ionization temperatures of Si and S determined from the intensity ratios of the RRC to He-like K-alpha line are 1.0 keV and 1.2 keV, respectively. We thus find firm evidence for an extremely-overionized (recombining) plasma. As the origin of the overionization, a thermal conduction scenario argued in previous work is not favored in our new results. We propose that the highly-ionized gas were made at the initial phase of the SNR evolution in dense regions around a massive progenitor, and the low electron temperature is due to a rapid cooling by an adiabatic expansion.

We show that higher-order actions for cosmological perturbations in the multi-field DBI-inflation model are obtained by a Lorentz boost from the rest frame of the brane to the frame where the brane is moving. We confirm that this simple method provides the same third- and fourth- order actions at leading order in slow-roll and in the small sound speed limit as those obtained by the usual ADM formalism. As an application, we compute the leading order connected four-point function of the primordial curvature perturbation coming from the intrinsic fourth-order contact interaction in the multi-field DBI-inflation model. At the third order, the interaction Hamiltonian arises purely by the boost from the second-order action in the rest frame of the brane. The boost acts on the adiabatic and entropy modes in the same way thus there exists a symmetry between the adiabatic and entropy modes. But at fourth order this symmetry is broken due to the intrinsic fourth-order action in the rest frame and the difference between the Lagrangian and the interaction Hamiltonian. Therefore, contrary to the three-point function, the momentum dependence of the purely adiabatic component and the components including the entropic contributions are different in the four-point function. This suggests that the trispectrum can distinguish the multi-field DBI-inflation model from the single field DBI-inflation model.

We study the effect of primordial isocurvature perturbations on non-Gaussian properties of CMB temperature anisotropies. We consider generic forms of the non-linearity of isocurvature perturbations which can be applied to a wide range of theoretical models. We derive analytical expressions for the bispectrum and the Minkowski Functionals for CMB temperature fluctuations to describe the non-Gaussianity from isocurvature perturbations. We find that the isocurvature non-Gaussianity in the quadratic isocurvature model, where the isocurvature perturbation S is written as a quadratic function of the Gaussian variable sigma, S=sigma^2-<sigma^2>, can give the same signal-to-noise as f_NL=30 even if we impose the current observational limit on the fraction of isocurvature perturbations contained in the primordial power spectrum alpha. We give constraints on isocurvature non-Gaussianity from Minkowski Functionals using WMAP 5-year data. We do not find a significant signal of the isocurvature non-Gaussianity. For the quadratic isocurvature model, we obtain a stringent upper limit on the isocurvature fraction alpha<0.070 (95% CL) for a scale invariant spectrum which is comparable to the limit obtained from the power spectrum.

We present a formalism to calculate the non-linear matter power spectrum in modified gravity models that explain the late-time acceleration of the Universe without dark energy. Any successful modified gravity models should contain a mechanism to recover General Relativity (GR) on small scales in order to avoid the stringent constrains on deviations from GR at solar system scales. Based on our formalism, the quasi non-linear power spectrum in the Dvali-Gabadadze-Porratti (DGP) braneworld models and $f(R)$ gravity models are derived by taking into account the mechanism to recover GR properly. We also extrapolate our predictions to fully non-linear scales using the Parametrized Post Friedmann (PPF) framework. In $f(R)$ gravity models, the predicted non-linear power spectrum is shown to reproduce N-body results. We find that the mechanism to recover GR suppresses the difference between the modified gravity models and dark energy models with the same expansion history, but the difference remains large at weakly non-linear regime in these models. Our formalism is applicable to a wide variety of modified gravity models and it is ready to use once consistent models for modified gravity are developed.

We compute the tree-level connected four-point function of the primordial curvature perturbation for a fairly general minimally coupled single field inflationary model, where the inflaton's Lagrangian is a general function of the scalar field and its first derivatives. This model includes K-inflation and DBI-inflation as particular cases. We show that, at the leading order in the slow-roll expansion and in the small sound speed limit, there are two important tree-level diagrams for the trispectrum. One is a diagram where a scalar mode is exchanged and the other is a diagram where the interaction occurs at a point, i.e. a contact interaction diagram. The scalar exchange contribution is comparable to the contact interaction contribution. For the DBI-inflation model, in the so-called equilateral configuration, the scalar exchange trispectrum is maximized when the angles between the four momentum vectors are equal and in this case the amplitude of the trispectrum from the scalar exchange is one order of magnitude higher than the contact interaction trispectrum.

A 20' \times 16' elliptical ring-like structure has been found near the Galactic center in the narrow energy band corresponding to the K$\alpha$ line from He-like sulfur. In the ring, two diffuse sources are found, a supernova remnant candidate G359.79-0.26 and an unidentified source G359.77-0.09. The X-ray spectrum of G359.77-0.09 is similar to that of G359.79-0.26, which can be explained by an absorbed thin thermal plasma model with temperatures of 0.7 and 1.0 keV. The absorption column densities of these two sources are large (N_H = 6.9 \times 10^{22} and 4.5 \times 10^{22} cm^{-2}) and are consistent with that of the Galactic center distance. The X-ray spectrum extracted from the ring-like structure is also represented by an absorbed thin thermal plasma model (kT_e \sim 0.9 keV). The thermal energy of the plasma filling the ring-like structure is estimated to be 1.0 \times 10^{51} erg. We therefore propose that the two sources comprise a single ring-like object, which is possibly a super bubble with a size of 49 pc \times 40 pc in the Galactic center region.

We present a new methodology to derive the positions of the Sgr B molecular clouds (MCs) along the line of sight, as an application study of the Galactic center diffuse X-rays (GCDX). The GCDX is composed of hot plasma emission of about 7 keV and 1 keV temperatures, and non-thermal continuum emission including the 6.4 keV line from neutral irons. The former, the Galactic center plasma emission (GCPE), is uniformly distributed over 1 degree in longitude, while the latter is clumpy emission produced by Thomson scattering and fluorescence from MCs irradiated by external X-rays (the X-ray reflection nebula emission: XRNE). We examined the Suzaku X-ray spectra of the GCPE and XRNE near to the Sgr B MC complex, and found that the spectra suffer from two different absorptions of N_H (Abs1) >10^23 H cm^-2 and N_H (Abs2) ~6 x 10^22 H cm^-2. Abs1 is proportional to the 6.4 keV-line flux, and hence is due to the MCs, while Abs2 is typical of interstellar absorption toward the Galactic center. Assuming that the GCPE plasma is spherically-extended around Sgr A* with a uniform density and the same angular distribution of the two temperature components, we quantitatively estimated the line-of-sight positions of the MCs from the flux ratio the GCPE spectrum suffered by Abs1 and that with no Abs1. The results suggest that the Sgr B MCs are located at the near side of Sgr A* in the GCPE.

We compute the non-Gaussianity of the curvature perturbation generated by ekpyrotic collapse with multiple fields. The transition from the multi-field scaling solution to a single-field dominated regime converts initial isocurvature field perturbations to an almost scale-invariant comoving curvature perturbation. In the specific model of two fields, $\phi_1$ and $\phi_2$, with exponential potentials, $-V_i \exp (-c_i \phi_i)$, we calculate the bispectrum of the resulting curvature perturbation. We find that the non-Gaussianity is dominated by non-linear evolution on super-Hubble scales and hence is of the local form. The non-linear parameter of the curvature perturbation is given by $f_{NL} = 5 c_j^2 /12$, where $c_j$ is the exponent of the potential for the field which becomes sub-dominant at late times. Since $c_j^2$ must be large, in order to generate an almost scale invariant spectrum, the non-Gaussianity is inevitably large. By combining the present observational constraints on $f_{\rm NL}$ and the scalar spectral index, the specific model studied in this paper is thus ruled out by current observational data.

A new and powerful probe of the origin and evolution of structures in the Universe has emerged and been actively developed over the last decade. In the coming decade, non-Gaussianity, i.e., the study of non-Gaussian contributions to the correlations of cosmological fluctuations, will become an important probe of both the early and the late Universe. Specifically, it will play a leading role in furthering our understanding of two fundamental aspects of cosmology and astrophysics: (i) the physics of the very early universe that created the primordial seeds for large-scale structures, and (ii) the subsequent growth of structures via gravitational instability and gas physics at later times. To date, observations of fluctuations in the Cosmic Microwave Background (CMB) and the Large-Scale Structure of the Universe (LSS) have focused largely on the Gaussian contribution as measured by the two-point correlations (or the power spectrum) of density fluctuations. However, an even greater amount of information is contained in non-Gaussianity and a large discovery space therefore still remains to be explored. Many observational probes can be used to measure non-Gaussianity, including CMB, LSS, gravitational lensing, Lyman-alpha forest, 21-cm fluctuations, and the abundance of rare objects such as clusters of galaxies and high-redshift galaxies. Not only does the study of non-Gaussianity maximize the science return from a plethora of present and future cosmological experiments and observations, but it also carries great potential for important discoveries in the coming decade.

We report on observations of SN 1006 with the X-ray Imaging Spectrometers (XIS) on board Suzaku. We firmly detected K-shell emission from Fe, for the first time, and find that the Fe ionization state is quite low. The broad band spectrum extracted from the southeast of the remnant is well fitted with a model consisting of three optically thin thermal non-equilibrium ionization plasmas and a power-law component. Two of the thermal models are highly overabundant in heavy elements and, hence, are likely due to ejecta. These components have different ionization parameters: $n_et \sim 1.4\times 10^{10}$ cm$^{-3}$ s and $n_et \sim 7.7\times 10^8$ cm$^{-3}$ s and it is the later one that produces the Fe-K emission. This suggests that Fe has been heated by the reverse shock more recently than the other elements, consistent with a picture where the ejecta are stratified by composition with Fe in the interior. On the other hand, the third thermal component is assumed to be solar abundance, and we associate it with emission from the interstellar medium (ISM). The electron temperature and ionization parameter are $kT_e \sim $0.5 keV and $n_et \sim 5.8\times 10^9$ cm$^{-3}$ s. The electron temperature is lower than that expected from the shock velocity which suggests a lack of collisionless electron heating at the forward shock. The extremely low ionization parameter and extreme non-equilibrium state are due to the low density of the ambient medium.

This paper reports the Suzaku results on the northeast shell of RCW 86. With the spatial and spectral analysis, we separated the X-rays into three distinct components; low (kT_e~0.3keV) and high (kT_e~1.8keV) temperature plasmas and a non-thermal component, and discovered their spatial distributions are different from each other. The low temperature plasma is dominated at the east rim, whereas the non-thermal emission is the brightest at the northeast rim which is spatially connected from the east rim. The high temperature plasma, found to contain the ~6.42keV line (K alpha of low-ionized iron), is enhanced at the inward region with respect to the east rim and has no spatial correlation with the non-thermal X-ray (the northeast). The Fe-Kalpha line, therefore, is not related to the non-thermal emission but originates from Fe-rich ejecta heated to the high temperatures by the reverse shock. Since the metal abundances of the low temperature plasma are sub-solar, the most possible origin of this component is interstellar medium heated by a blast wave. The non-thermal X-ray, which has a power-law index of ~2.8, is likely to be synchrotron emission. A possible scenario to explain these morphologies and spectra is: A fast moving blast wave in a thin cavity of OB association collided with a dense interstellar medium or cloud at the east region very recently. As the result, the reverse shock in this interior decelerated, and arrived at the Fe-rich region of the ejecta and heated it. In the northeast rim, on the other hand, the blast wave is still moving fast, and accelerated high energy electrons to emit synchrotron X-rays.

This paper reports on the Suzaku results of thermal and non-thermal features of 30 Dor C, a supernova remnant (SNR) in a superbubble of the Large Magellanic Cloud (LMC). The west rim exhibits a non-thermal X-ray spectrum with no thermal component. A single power-law model is rejected but a power-law model with spectral cutoff is accepted. The cutoff frequency of $(3-7)\times 10^{17}$ Hz is the highest among the shell type SNRs like SN 1006 ($\sim 6\times 10^{16}$ Hz), and hence 30 Dor C would be the site of the highest energy accelerator of the SNR shock. The southeast (SE) and northeast (NE) rims have both the thermal and non-thermal components. The thin-thermal plasmas in the both rims are in collisional ionization equilibrium state. The electron temperature of the plasma in the SE rim ($kT_e \sim 0.7$ keV) is found to be higher than the previously reported value. The power-law index from SE is nearly the same as, while that from the NE is larger than that of the West rim. The SNR age would be in the range of $(4-20)\times 10^3$ yr. Thus, 30 Dor C is likely to be the oldest shell-like SNR with non-thermal emission.

Suzaku observed the Galactic center region near the Radio Arc at ~20' southeast of Sagittarius A*.In the 18'x18' field of view, we found four distinct X-ray sources: a bright star and a diffuse source associated with the star clusters in the soft band (0.5-2.0 keV), a small clump in a higher energy band (4-6 keV), and a peculiar clump in the 6.4 keV line band.The latter two clumps are located at the south end of the Radio Arc. This paper reports on the results, and discusses the origin of these X-ray sources, with a particular emphasis on small clumps.

We present the Suzaku results on a new candidate of a supernova remnant (SNR) in the Sagittarius C region. We detected diffuse X-rays of an elliptical shape (G359.41-0.12) and a chimney-like structure (the Chimney), both of which were fitted with a thin thermal the model of kBT ~1 keV temperature. The absorption columns are same between these two structures, indicating that both are located at the same distance in the same line of sight. The narrow band image and one-dimensional profile of S XV Kalpha at 2.45 keV show that the Chimney is emanating from G359.41-0.12. Therefore, these two sources are physically connected with each other. The sum of the thermal energies of the Chimney and G359.41-0.12 is estimated to be 1.4x10^50 erg, typical for a galactic SNR. G359.41-0.12 is likely to be a new SNR candidate and the Chimney is an associated outflow.

A scale-invariant spectrum of isocurvature perturbations is generated during collapse in the ekpyrotic scaling solution in models where multiple fields have steep negative exponential potentials. The scale invariance of the spectrum is realized by a tachyonic instability in the isocurvature field. This instability drives the scaling solution to the late time attractor that is the old ekpyrotic collapse dominated by a single field. We show that the transition from the scaling solution to the single field dominated ekpyrotic collapse automatically converts the initial isocurvature perturbations about the scaling solution to comoving curvature perturbations about the late-time attractor. The final amplitude of the comoving curvature perturbation is determined by the Hubble scale at the transition.

We construct a consistency test of General Relativity (GR) on cosmological scales. This test enables us to distinguish between the two alternatives to explain the late-time accelerated expansion of the universe, that is, dark energy models based on GR and modified gravity models without dark energy. We derive the consistency relation in GR which is written only in terms of observables - the Hubble parameter, the density perturbations, the peculiar velocities and the lensing potential. The breakdown of this consistency relation implies that the Newton constant which governs large-scale structure is different from that in the background cosmology, which is a typical feature in modified gravity models. We propose a method to perform this test by reconstructing the weak lensing spectrum from measured density perturbations and peculiar velocities. This reconstruction relies on Poisson's equation in GR to convert the density perturbations to the lensing potential. Hence any inconsistency between the reconstructed lensing spectrum and the measured lensing spectrum indicates the failure of GR on cosmological scales. The difficulties in performing this test using actual observations are discussed.

We study the spectrum of gravitational perturbations around a vacuum de Sitter brane in a 5D asymmetric braneworld model, with induced curvature on the brane. This generalises the stealth acceleration model proposed by Charmousis, Gregory and Padilla (CGP) which realises the Cardassian cosmology in which power law cosmic acceleration can be driven by ordinary matter. Whenever the bulk has infinite volume we find that there is always a perturbative ghost propagating on the de Sitter brane, in contrast to the Minkowski brane case analysed by CGP. We discuss the implication of this ghost for the stealth acceleration model, and identify a limiting case where the ghost decouples as the de Sitter curvature vanishes.