We study the structure of neutron stars in f(R) gravity theories with perturbative constraints. We derive the modified Tolman-Oppenheimer-Volkov equations and solve them for a polytropic equation of state. We investigate the resulting modifications to the masses and radii of neutron stars and show that observations of surface phenomena alone cannot break the degeneracy between altering the theory of gravity versus choosing a different equation of state of neutron-star matter. On the other hand, observations of neutron-star cooling, which depends on the density of matter at the stellar interior, can place significant constraints on the parameters of the theory.

Energy-dependent speeds of light have been considered an observable signature of quantum gravity effects. The two simplest dispersion relationships produce either linear or quadratic corrections, in particle energy, to the photon speed. The macroscopic limits of these theories -- how objects with small energy per particle, but with large mass, behave -- are not fully understood. We here briefly discuss some features of the macroscopic limit, that are necessary for understanding how detectors and emitters interact with the high-energy photons that probe spacetime.

We present a new method for detecting the missing baryons by generating a template for the kinematic Sunyaev-Zel'dovich effect. The template is computed from the product of a reconstructed velocity field with a galaxy field; we find that the combination of a galaxy redshift survey such as SDSS and a CMB survey such as ACT and PLANCK can detect the kSZ, and thus the ionized gas, at significant signal-to-noise. Unlike other techniques that look for hot gas or metals, this approach directly detects the electrons in the IGM through their signature on the CMB. The estimated signal-to-noise for detecting the galaxy-momentum kSZ cross-correlation is 4, 9, and 12 for ACT (with survey area of 2000 $\mathrm{deg}^2$) with SDSS-DR4, SDSS3 and ADEPT respectively. The estimated signal-to-noise for PLANCK with SDSS-DR4, SDSS3 and ADEPT is 11, 23, and 32. Our method provides a new mean for determining properties of the ionized gas in the Universe. We provide galaxy momentum templates constructed from Sloan Digital Sky Survey online at our website at http://www.astro.princeton.edu/~shirley/SZ/SZ.html. The predicted correlation coefficients are provided along with the momentum maps. One can download the momentum templates and cross-correlate directly with CMB maps from ACT and PLANCK to detect the missing baryons.

Major observational efforts in the coming decade are designed to probe the equation of state of dark energy. Measuring a deviation of the equation-of-state parameter w from -1 would indicate a dark energy that cannot be represented solely by a cosmological constant. While it is commonly assumed that any implied modification to the LambdaCDM model amounts to the addition of new dynamical fields, we propose here a framework for investigating whether or not such new fields are required when cosmological observations are combined with a set of minimal assumptions about the nature of gravitational physics. In our approach, we treat the additional degrees of freedom as perturbatively constrained and calculate a number of observable quantities, such as the Hubble expansion rate and the cosmic acceleration, for a homogeneous Universe. We show that current observations place our Universe within the perturbative validity of our framework and allow for the presence of non-dynamical gravitational degrees of freedom at cosmological scales.

Modifications to gravity that add additional functions of the Ricci curvature to the Einstein-Hilbert action -- collectively known as $f(R)$ theories -- have been studied in great detail. When considered as complete theories of gravity they can generate non-perturbative deviations from the general relativistic predictions in the solar system, and the simplest models show instabilites on cosmological scales. Here we show that it is possible to treat $f(R)=R\pm\mu^4/R$ gravity in a perturbative fashion such that it shows no instabilities on cosmological scales and, in the solar system, is consistent with measurements of the PPN parameters. We show that such a theory produces a spatially flat, accelerating universe, even in the absence of dark energy and when the matter density is too small to close the universe in the general relativistic case.

Upcoming Sunyaev-Zel'dovich surveys are expected to return ~10^4 intermediate mass clusters at high redshift. Their average masses must be known to same accuracy as desired for the dark energy properties. Internal to the surveys, the CMB potentially provides a source for lensing mass measurements whose distance is precisely known and behind all clusters. We develop statistical mass estimators from 6 quadratic combinations of CMB temperature and polarization fields that can simultaneously recover large-scale structure and cluster mass profiles. The performance of these estimators on idealized NFW clusters suggests that surveys with a ~1' beam and 10uK' noise in uncontaminated temperature maps can make a ~10sigma detection, or equivalently a ~10% mass measurement for each 10^3 set of clusters. With internal or external acoustic scale E-polarization measurements, the ET cross correlation estimator can provide a stringent test for contaminants on a first detection at \~1/3 the significance. For surveys that reach below 3muK', the EB cross correlation estimator should provide the most precise measurements and potentially the strongest control over contaminants.

We explore the physics of a gas of particles interacting with a condensate that spontaneously breaks Lorentz invariance. The equation of state of this gas varies from 1/3 to less than -1 and can lead to the observed cosmic acceleration. The particles are always stable. In our particular class of models these particles are fermions with a chiral coupling to the condensate. They may behave as relativistic matter at early times, produce a brief period where they dominate the expansion with w<0 today, and behave as matter at late time. There are no small parameters in our models, which generically lead to dark energy clustering and, depending on the choice of parameters, smoothing of small scale power.

Upcoming observatories will be able to detect the kinetic Sunyaev-Zel'dovitch (kSZ) effect with unprecendented signal-to-noise, and cross-correlations with foreground signals such as galaxy counts are a promising way to extract additional cosmological information. We consider how well a tomographic galaxy-count cross-correlation experiment, using data from WMAP, ACT and SALT, can significantly constrain the properties of dark energy. We include the need to model a wide range of effects, including those associated with complicated baryonic physics, in our analysis. We demonstrate how much of the cosmological information contained in the kSZ comes from larger scales than that in the galaxy power spectrum, and thus how use of the kSZ can help avoid difficult systematics associated with non-linear and scale-dependent bias at k>1h Mpc^{-1}.

The frequencies of quasi-periodic oscillations around neutron stars are believed to be related to characteristic frequencies in the gravitational fields of the compact objects. In different variability models, these include the Keplerian, epicyclic, and Lense-Thirring frequencies, which depend mostly on the properties of the stellar spacetimes. We argue that quasi-periodic oscillations in the X-ray flux of neutron stars can be used to map the external spacetimes of the compact objects and, therefore, lead to direct tests of general relativity in the strong-field regime. In particular, we show that particular extensions of General Relativity, in which the gravitational force felt by matter is mediated by both a rank-two tensor and a scalar field, can be constrained by current observations.

We consider how quintessence models in which the sound speed differs from the speed of light and varies with time affect the cosmic microwave background and the fluctuation power spectrum. Significant modifications occur on length scales related to the Hubble radius during epochs in which the sound speed is near zero and the quintessence contributes a non-negligible fraction of the total energy density. For the microwave background, we find that the usual enhancement of the lowest multipole moments by the integrated Sachs-Wolfe effect can be modified, resulting in suppression or bumps instead. Also, the sound speed can produce oscillations and other effects at wavenumbers $k > 10^{-2}$ h/Mpc in the fluctuation power spectrum.

In contrast to gravity in the weak-field regime, which has been subject to numerous experimental tests, gravity in the strong-field regime is largely unconstrained by observations. We show that gravity theories that cannot be rejected by solar system tests but that diverge from general relativity in the strong-field regime predict neutron stars with significantly different properties than their general relativistic counterparts. In particular, the range of redshifts of surface atomic lines predicted by such gravity theories is significantly larger than the uncertainty introduced by our lack of knowledge of the equation of state of ultra-dense matter. Measurements of such redshifted lines with current X-ray observatories such as Chandra and XMM-Newton can thus provide interesting new constraints on strong-field gravity.

Anomalous X-ray pulsars, compact non-pulsing X-ray sources in supernova remnants, and X-ray bursters are three distinct types of sources for which there are viable models that attribute their X-ray emission to thermal emission from the surface of a neutron star. Inferring the surface area of the emitting regions in such systems is crucial in assessing the viability of different models and in providing bounds on the radii of neutron stars. We show that the inferred areas of the emitting regions may be over- or under-estimated by a factor of <=2, because of the geometry of the system and general relativistic light deflection, combined with the effects of phase averaging. Such effects make the determination of neutron-star radii uncertain, especially when compared to the ~5% level required for constraining the equation of state of neutron-star matter. We also note that, for a given spectral shape, the inferred source luminosities and pulse fractions are anticorrelated because they depend on the same properties of the emitting regions, namely their sizes and orientations, i.e., brighter sources have on average weaker pulsation amplitudes than fainter sources. We argue that this property can be used as a diagnostic tool in distinguishing between different spectral models. As an example, we show that the high inferred pulse fraction and brightness of the pulsar RXS J1708-40 are inconsistent with isotropic thermal emission from a neutron-star surface. Finally, we discuss the implication of our results for surveys in the soft X-rays for young, cooling neutron stars in supernova remnants and show that the absence of detectable pulsations from the compact source at the center of Cas A (at a level of >=30%) is not a strong argument againts its identification with a spinning neutron star.

Most models of anomalous X-ray pulsars (AXPs) account for the observed X-ray spectra and pulsations by means of radiation processes that occur on the surfaces of neutron stars. For any such model, general relativistic deflection of light severely suppresses the amplitude of the observed pulsations. We calculate the expected pulsation amplitudes of AXPs according to various models and compare the results with observations. We show that the high (<= 70%) pulse amplitudes observed in some AXPs can be accounted for only if the surface emission is localized (spot radius <40 degrees) and strongly beamed (cos^n[theta'] with n>2, where theta' is the angle to the normal). These constraints are incompatible with those cooling and magnetar models in which the observed X-rays originate as thermal emission from the neutron-star surface. Accretion models, on the other hand, are compatible with observations for a wide range of parameters. Finally, definitive conclusions cannot be reached on magnetospheric models, since their localization and beaming properties are not well understood.