It is commonly assumed that the energy density of the Universe was dominated by radiation between reheating after inflation and the onset of matter domination 54,000 years later. While the abundance of light elements indicates that the Universe was radiation dominated during Big Bang Nucleosynthesis (BBN), there is scant evidence that the Universe was radiation dominated prior to BBN. It is therefore possible that the cosmological history was more complicated, with deviations from the standard radiation domination during the earliest epochs. Indeed, several interesting proposals regarding various topics such as the generation of dark matter, matter-antimatter asymmetry, gravitational waves, primordial black holes, or microhalos during a nonstandard expansion phase have been recently made. In this paper, we review various possible causes and consequences of deviations from radiation domination in the early Universe - taking place either before or after BBN - and the constraints on them, as they have been discussed in the literature during the recent years.

LCDM cosmological models with Early Dark Energy (EDE) have been proposed to resolve tensions between the Hubble constant H0 = 100h km/s/Mpc measured locally, giving h ~ 0.73, and H0 deduced from Planck cosmic microwave background (CMB) and other early universe measurements plus LCDM, giving h ~ 0.67. EDE models do this by adding a scalar field that temporarily adds dark energy equal to about 10% of the cosmological energy density at the end of the radiation-dominated era at redshift z ~ 3500. Here we compare linear and nonlinear predictions of a Planck-normalized LCDM model including EDE giving h = 0.728 with those of standard Planck-normalized LCDM with h = 0.678. We find that nonlinear evolution reduces the differences between power spectra of fluctuations at low redshifts. As a result, at z = 0 the halo mass functions on galactic scales are nearly the same, with differences only 1-2%. However, the differences dramatically increase at high redshifts. The EDE model predicts 50% more massive clusters at z = 1 and twice more galaxy-mass halos at z = 4. Even greater increases in abundances of galaxy-mass halos at higher redshifts may make it easier to reionize the universe with EDE. Predicted galaxy abundances and clustering will soon be tested by JWST observations. Positions of baryonic acoustic oscillations (BAOs) and correlation functions differ by about 2% between the models -- an effect that is not washed out by nonlinearities. Both standard LCDM and the EDE model studied here agree well with presently available acoustic-scale observations, but DESI and Euclid measurements will provide stringent new tests.

Pulsar timing arrays (PTAs) are currently the only experiments directly sensitive to gravitational waves with decade-long periods. Within the next five to ten years, PTAs are expected to detect the stochastic gravitational-wave background (SGWB) collectively sourced by inspiralling supermassive black hole binaries. It is expected that this background is mostly isotropic, and current searches focus on the monopole part of the SGWB. Looking ahead, anisotropies in the SGWB may provide a trove of additional information both on known and unknown astrophysical and cosmological sources. In this paper, we build a simple yet realistic Fisher formalism for anisotropic SGWB searches with PTAs. Our formalism is able to accommodate realistic properties of PTAs, and allows simple and accurate forecasts. We illustrate our approach with an idealized PTA consisting of identical, isotropically distributed pulsars. In a companion paper, we apply our formalism to current PTAs and show that it can be a powerful tool to guide and optimize real data analysis.

Baryon acoustic oscillations (BAO) provide a robust standard ruler, and can be used to constrain the expansion history of the Universe at low redshift. Standard BAO analyses return a model-independent measurement of the expansion rate and the angular diameter distance as function of redshift, normalized by the sound horizon at radiation drag. However, this methodology relies on anisotropic distance distortions of a fixed, pre-computed template (obtained in a given fiducial cosmology) in order to fit the observations. Therefore, it may be possible that extensions to the consensus $\Lambda$CDM add contributions to the BAO feature that cannot be captured by the template fitting. We perform mock BAO fits to power spectra computed assuming cosmological models which modify the growth of perturbations prior to recombination in order to test the robustness of the standard BAO analysis. We find no significant bias in the BAO analysis for the models under study ($\Lambda$CDM with a free effective number of relativistic species, early dark energy, and a model with interactions between neutrinos and a fraction of the dark matter), even for cases which do not provide a good fit to Planck measurements of the CMB power spectra. This result supports the use of the standard BAO analysis and its measurements to perform cosmological parameter inference and to constrain exotic models. In addition, we provide a methodology to reproduce our study for different models and surveys, as well as discuss different options to handle eventual biases in the BAO measurements.

Ultra-light axions (ULAs) are a promising dark-matter candidate. ULAs may have implications for small-scale challenges to the $\Lambda$CDM model, and arise in string scenarios. ULAs are already constrained by cosmic microwave background (CMB) experiments and large-scale structure surveys, and will be probed with much greater sensitivity by future efforts. It is challenging to compute observables in ULA scenarios with sufficient speed and accuracy for cosmological data analysis because the ULA field oscillates rapidly. In past work, an effective fluid approximation has been used to make these computations feasible. Here this approximation is tested against an exact solution of the ULA equations, comparing the induced error of CMB observables with the sensitivity of current and future experiments. In the most constrained mass range for a ULA dark matter component ($10^{-27}~{\rm eV}\leq m_{\rm ax}\leq 10^{-25}~{\rm eV}$), the induced bias on the allowed ULA fraction of dark matter from Planck data is less than $1\sigma$. In the cosmic-variance limit (including temperature and polarization data), the bias is $\lesssim 2\sigma$ for primary CMB anisotropies, with more severe biases (as high as $\sim 4\sigma$) resulting for less reliable versions of the effective fluid approximation. If all of the standard cosmological parameters are fixed by other measurements, the expected bias rises to $4-20\sigma$ (well beyond the validity of the Fisher approximation), though the required level of degeneracy breaking will not be achieved by any planned surveys.

We present the steps to forecast the sensitivity of the Laser Interferometer Space Antenna (LISA) to both a stochastic gravitational wave background and deterministic wave sources. We show how to use these expressions to estimate the precision with which LISA can determine parameters associated with these sources. Tools are included to enable easy calculation of the signal-to-noise ratio and draw sensitivity curves. Benchmark values are given for easy comparison and checking of methods in the case of two worked examples. The first benchmark is the threshold stochastic gravitational wave background $\Omega_{GW} h^2$ that LISA can observe. The second is the signal-to-noise ratio that LISA would observe for a binary black hole system identical to GW150914, radiating 4 years before merger.

We present a detailed investigation of a sub-dominant oscillating scalar field ('early dark energy', EDE) in the context of resolving the Hubble tension. Consistent with earlier work, but without relying on fluid approximations, we find that a scalar field frozen due to Hubble friction until ${\rm log}_{10}(z_c)\sim3.5$, reaching $\rho_{\rm EDE}(z_c)/\rho_{\rm tot}\sim10$%, and diluting faster than matter afterwards can bring cosmic microwave background (CMB), baryonic acoustic oscillations, supernovae luminosity distances, and the late-time estimate of the Hubble constant from the SH0ES collaboration into agreement. A scalar field potential which scales as $V(\phi) \propto \phi^{2n}$ with $2\lesssim n\lesssim 3.4$ around the minimum is preferred at the 68% confidence level, and the {\em Planck} polarization places additional constraints on the dynamics of perturbations in the scalar field. In particular, the data prefers a potential which flattens at large field displacements. An MCMC analysis of mock data shows that the next-generation CMB observations (i.e., CMB-S4) can unambiguously detect the presence of the EDE at very high significance. This projected sensitivity to the EDE dynamics is mainly driven by improved measurements of the $E$-mode polarization. We also explore new observational signatures of EDE scalar field dynamics: (i) We find that depending on the strength of the tensor-to-scalar ratio, the presence of the EDE might imply the existence of isocurvature perturbations in the CMB. (ii) We show that a strikingly rapid, scale-dependent growth of EDE field perturbations can result from parametric resonance driven by the anharmonic oscillating field for $n\approx 2$. This instability and ensuing potentially nonlinear, spatially inhomogenoues, dynamics may provide unique signatures of this scenario.

Pulsar timing arrays (PTAs) will enable the detection of nanohertz gravitational waves (GWs) from a population of supermassive binary black holes (SMBBHs) in the next $\sim 3-7$ years. In addition, PTAs provide a rare opportunity to probe exotic physics. Potential sources of GWs in the nanohertz band include, cosmic strings and cosmic superstrings, inflation and, phase transitions in the early universe. GW observations will also make possible tests of gravitational theories that, by modifying Einstein's theory of general relativity, attempt to explain the origin of cosmic acceleration and reconcile quantum mechanics and gravity, two of the most profound challenges facing fundamental physics today. Finally, PTAs also provide a new means to probe certain dark matter models.

We construct realistic sensitivity curves for pulsar timing array searches for gravitational waves, incorporating both red and white noise contributions to individual pulsar noise spectra, as well as the effect of fitting to a pulsar timing model. We demonstrate the method on both simulated pulsars and a realistic array consisting of a subset of NANOGrav pulsars used in recent analyses. A comparison between the results presented here and exact sensitivity curves shows agreement to tens of percent. The resulting sensitivity curves can be used to assess the detectability of predicted gravitational-wave signals in the nanohertz frequency band in a fraction of the time that it would take to compute the exact sensitivity curves.

Early dark energy (EDE) that behaves like a cosmological constant at early times (redshifts $z\gtrsim3000$) and then dilutes away like radiation or faster at later times can solve the Hubble tension. In these models, the sound horizon at decoupling is reduced resulting in a larger value of the Hubble parameter $H_0$ inferred from the cosmic microwave background (CMB). We consider two physical models for this EDE, one involving an oscillating scalar field and another a slowly-rolling field. We perform a detailed calculation of the evolution of perturbations in these models. A Markov Chain Monte Carlo search of the parameter space for the EDE parameters, in conjunction with the standard cosmological parameters, identifies regions in which $H_0$ inferred from {\it Planck} CMB data agrees with the SH0ES local measurement. In these cosmologies, current baryon acoustic oscillation and supernova data are described as successfully as in \LCDM, while the fit to {\it Planck} data is slightly improved. Future CMB and large-scale-structure surveys will further probe this scenario.

Pulsar timing array (PTA) collaborations in North America, Australia, and Europe, have been exploiting the exquisite timing precision of millisecond pulsars over decades of observations to search for correlated timing deviations induced by gravitational waves (GWs). PTAs are sensitive to the frequency band ranging just below 1 nanohertz to a few tens of microhertz. The discovery space of this band is potentially rich with populations of inspiraling supermassive black-holes binaries, decaying cosmic string networks, relic post-inflation GWs, and even non-GW imprints of axionic dark matter. This article aims to provide an understanding of the exciting open science questions in cosmology, galaxy evolution, and fundamental physics that will be addressed by the detection and study of GWs through PTAs. The focus of the article is on providing an understanding of the mechanisms by which PTAs can address specific questions in these fields, and to outline some of the subtleties and difficulties in each case. The material included is weighted most heavily towards the questions which we expect will be answered in the near-term with PTAs; however, we have made efforts to include most currently anticipated applications of nanohertz GWs.

The axion is a hypothetical, well-motivated dark-matter particle whose existence would explain the lack of charge-parity violation in the strong interaction. In addition to this original motivation, an `axiverse' of ultra-light axions (ULAs) with masses $10^{-33}\,{\rm eV}\lesssim m_{\rm a}\lesssim 10^{-10}\,{\rm eV}$ also emerges from string theory. Depending on the mass, such a ULA contributes to the dark-matter density, or alternatively, behaves like dark energy. At these masses, ULAs' classical wave-like properties are astronomically manifested, potentially mitigating observational tensions within the $\Lambda$CDM paradigm on local-group scales. ULAs also provide signatures on small scales such as suppression of structure, interference patterns and solitons to distinguish them from heavier dark matter candidates. Through their gravitational imprint, ULAs in the presently allowed parameter space furnish a host of observational tests to target in the next decade, altering standard predictions for microwave background anisotropies, galaxy clustering, Lyman-$\alpha$ absorption by neutral hydrogen along quasar sightlines, pulsar timing, and the black-hole mass spectrum.

There are two big questions cosmologists would like to answer -- How does the Universe work, and what are its origin and destiny? A long wavelength gravitational wave detector -- with million km interferometer arms, achievable only from space -- gives a unique opportunity to address both of these questions. A sensitive, mHz frequency observatory could use the inspiral and merger of massive black hole binaries as standard sirens, extending our ability to characterize the expansion history of the Universe from the onset of dark energy-domination out to a redshift z ~ 10. A low-frequency detector, furthermore, offers the best chance for discovery of exotic gravitational wave sources, including a primordial stochastic background, that could reveal clues to the origin of our Universe.

A primordial spectrum of gravitational waves serves as a backlight to the relativistic degrees of freedom of the cosmological fluid. Any change in the particle physics content, due to a change of phase or freeze-out of a species, will leave a characteristic imprint on an otherwise featureless primordial spectrum of gravitational waves and indicate its early-Universe provenance. We show that a gravitational wave detector such as the Laser Interferometer Space Antenna would be sensitive to physics near 100 TeV in the presence of a sufficiently strong primordial spectrum. Such a detection could complement searches at newly proposed 100 km circumference accelerators such as the Future Circular Collider at CERN and the Super Proton-Proton Collider in China, thereby providing insight into a host of beyond Standard Model issues, including the hierarchy problem, dark matter, and baryogenesis.

The fine-structure constant, $\alpha$, controls the strength of the electromagnetic interaction. There are extensions of the standard model in which $\alpha$ is dynamical on cosmological length and time-scales. The physics of the cosmic microwave background (CMB) depends on the value of $\alpha$. The effects of spatial variation in $\alpha$ on the CMB are similar to those produced by weak lensing: smoothing of the power spectrum, and generation of non-Gaussian features. These would induce a bias to estimates of the weak-lensing potential power spectrum of the CMB. Using this effect, Planck measurements of the temperature and polarization power spectrum, as well as estimates of CMB lensing, are used to place limits (95% C. L.) on the amplitude of a scale-invariant angular power spectrum of $\alpha$ fluctuations relative to the mean value ($C_L^{\alpha}=A^\alpha_{\rm SI}/[L(L+1)]$) of $A^\alpha_{\rm SI}\leq 1.6 \times 10^{-5}$. The limits depend on the assumed shape of the $\alpha$-fluctuation power spectrum. For example, for a white noise angular power spectrum ($C_L^{\alpha}=A^\alpha_{\rm WN}$), the limit is $A^\alpha_{\rm WN}\leq 2.3 \times 10^{-8}$. It is found that the response of the CMB to $\alpha$ fluctuations depends on a separate-universe approximation, such that theoretical predictions are only reliable for $\alpha$ multipoles with $L\lesssim 100$ . An optimal trispectrum estimator can be constructed and it is found that it is only marginally more sensitive than lensing techniques for Planck but significantly more sensitive when considering the next generation of experiments. For a future CMB experiment with cosmic-variance limited polarization sensitivity (e.g., CMB-S4), the optimal estimator could detect $\alpha$ fluctuations with $A^\alpha_{\rm SI}>1.9 \times 10^{-6}$ and $A^\alpha_{\rm WN} > 1.4 \times 10^{-9}$.

Cosmological observations are used to test for imprints of an ultra-light axion-like field (ULA), with a range of potentials $V(\phi)\propto[1-\cos(\phi/f)]^n$ set by the axion-field value $\phi$ and decay constant $f$. Scalar field dynamics dictate that the field is initially frozen and then begins to oscillate around its minimum when the Hubble parameter drops below some critical value. For $n\!=\!1$, once dynamical, the axion energy density dilutes as matter; for $n\!=\!2$ it dilutes as radiation and for $n\!=\!3$ it dilutes faster than radiation. Both the homogeneous evolution of the ULA and the dynamics of its linear perturbations are included, using an effective fluid approximation generalized from the usual $n=1$ case. ULA models are parameterized by the redshift $z_c$ when the field becomes dynamical, the fractional energy density $f_{z_c} \equiv \Omega_a(z_c)/\Omega_{\rm tot}(z_c)$ in the axion field at $z_c$, and the effective sound speed $c_s^2$. Using Planck, BAO and JLA data, constraints on $f_{z_c}$ are obtained. ULAs are degenerate with dark energy for all three potentials if $1+z_c \lesssim 10$. When $3\times10^4 \gtrsim 1+z_c \gtrsim 10 $, $f_{z_c}$ is constrained to be $ \lesssim 0.004 $ for $n=1$ and $f_{z_c} \lesssim 0.02 $ for the other two potentials. The constraints then relax with increasing $z_c$. These results strongly constrain ULAs as a resolution to cosmological tensions, such as discrepant measurements of the Hubble constant, or the EDGES measurement of the global 21 cm signal.

We present astroplan - an open source, open development, Astropy affiliated package for ground-based observation planning and scheduling in Python. astroplan is designed to provide efficient access to common observational quantities such as celestial rise, set, and meridian transit times and simple transformations from sky coordinates to altitude-azimuth coordinates without requiring a detailed understanding of astropy's implementation of coordinate systems. astroplan provides convenience functions to generate common observational plots such as airmass and parallactic angle as a function of time, along with basic sky (finder) charts. Users can determine whether or not a target is observable given a variety of observing constraints, such as airmass limits, time ranges, Moon illumination/separation ranges, and more. A selection of observation schedulers are included which divide observing time among a list of targets, given observing constraints on those targets. Contributions to the source code from the community are welcome.

Compensated isocurvature perturbations (CIPs) are primordial fluctuations that balance baryon and dark-matter isocurvature to leave the total matter density unperturbed. The effects of CIPs on the cosmic microwave background (CMB) anisotropies are similar to those produced by weak lensing of the CMB: smoothing of the power spectrum, and generation of non-Gaussian features. Previous work considered the CIP effects on the CMB power-spectrum but neglected to include the CIP effects on estimates of the lensing potential power spectrum (though its contribution to the non-Gaussian, connected, part of the CMB trispectrum). Here, the CIP contribution to the standard estimator for the lensing potential power-spectrum is derived, and along with the CIP contributions to the CMB power-spectrum, Planck data is used to place limits on the root-mean-square CIP fluctuations on CMB scales, $\Delta_{\rm rms}^2(R_{\rm CMB})$. The resulting constraint of $\Delta_{\rm rms}^2(R_{\rm CMB}) < 4.3 \times 10^{-3}$ using this new technique improves on past work by a factor of $\sim 3$. We find that for Planck data our constraints almost reach the sensitivity of the optimal CIP estimator. The method presented here is currently the most sensitive probe of the amplitude of a scale-invariant CIP power spectrum placing an upper limit of $A_{\rm CIP}< 0.017$ at 95% CL. Future measurements of the large-scale CMB lensing potential power spectrum could probe CIP amplitudes as low as $\Delta_{\rm rms}^2(R_{\rm CMB}) = 8 \times 10^{-5}$ ($A_{\rm CIP} = 3.2 \times 10^{-4}$).

This book lays out the scientific goals to be addressed by the next-generation ground-based cosmic microwave background experiment, CMB-S4, envisioned to consist of dedicated telescopes at the South Pole, the high Chilean Atacama plateau and possibly a northern hemisphere site, all equipped with new superconducting cameras. CMB-S4 will dramatically advance cosmological studies by crossing critical thresholds in the search for the B-mode polarization signature of primordial gravitational waves, in the determination of the number and masses of the neutrinos, in the search for evidence of new light relics, in constraining the nature of dark energy, and in testing general relativity on large scales.

We calculate the sensitivity to a circular polarization of an isotropic stochastic gravitational wave background (ISGWB) as a function of frequency for ground- and space-based interferometers and observations of the cosmic microwave background. The origin of a circularly polarized ISGWB may be due to exotic primordial physics (i.e., parity violation in the early universe) and may be strongly frequency dependent. We present calculations within a coherent framework which clarifies the basic requirements for sensitivity to circular polarization, in distinction from previous work which focused on each of these techniques separately. We find that the addition of an interferometer with the sensitivity of the Einstein Telescope in the southern hemisphere improves the sensitivity of the ground-based network to circular polarization by about a factor of two. The sensitivity curves presented in this paper make clear that the wide range in frequencies of current and planned observations ($10^{-18}\ {\rm Hz} \lesssim f \lesssim 100\ {\rm Hz}$) will be critical to determining the physics that underlies any positive detection of circular polarization in the ISGWB. We also identify a desert in circular polarization sensitivity for frequencies between $10^{-15}\ {\rm Hz} \lesssim f \lesssim 10^{-3}\ {\rm Hz}$, given the inability for pulsar timing arrays and indirect-detection methods to distinguish the gravitational wave polarization.

Quantum fluctuations of the gravitational field in the early Universe, amplified by inflation, produce a primordial gravitational-wave background across a broad frequency band. We derive constraints on the spectrum of this gravitational radiation, and hence on theories of the early Universe, by combining experiments that cover 29 orders of magnitude in frequency. These include Planck observations of cosmic microwave background temperature and polarization power spectra and lensing, together with baryon acoustic oscillations and big bang nucleosynthesis measurements, as well as new pulsar timing array and ground-based interferometer limits. While individual experiments constrain the gravitational-wave energy density in specific frequency bands, the combination of experiments allows us to constrain cosmological parameters, including the inflationary spectral index, $n_t$, and the tensor-to-scalar ratio, $r$. Results from individual experiments include the most stringent nanohertz limit of the primordial background to date from the Parkes Pulsar Timing Array, $\Omega_{\rm gw}(f)<2.3\times10^{-10}$. Observations of the cosmic microwave background alone limit the gravitational-wave spectral index at 95\% confidence to $n_t\lesssim5$ for a tensor-to-scalar ratio of $r = 0.11$. However, the combination of all the above experiments limits $n_t<0.36$. Future Advanced LIGO observations are expected to further constrain $n_t<0.34$ by 2020. When cosmic microwave background experiments detect a non-zero $r$, our results will imply even more stringent constraints on $n_t$ and hence theories of the early Universe.

In the curvaton scenario, primordial curvature perturbations are produced by a second field that is sub-dominant during inflation. Depending on how the curvaton decays [possibly producing baryon number, lepton number, or cold dark matter (CDM)], mixtures of correlated isocurvature perturbations are produced, allowing the curvaton scenario to be tested using cosmic microwave background (CMB) data. Here, a full range of 27 curvaton-decay scenarios is compared with CMB data, placing limits on the curvaton fraction at decay, $r_D$, and the lepton asymmetry, $\xi_{\rm lep}$. If baryon number is generated by curvaton decay and CDM before (or vice-versa), these limits imply specific predictions for non-Gaussian signatures testable by future CMB experiments and upcoming large-scale-structure surveys.

We provide APASS photometry in the Landolt BV and Sloan g'r'i' bands for all the 425,743 stars included in the latest 4th RAVE Data Release. The internal accuracy of the APASS photometry of RAVE stars, expressed as error of the mean of data obtained and separately calibrated over a median of 4 distinct observing epochs and distributed between 2009 and 2013, is 0.013, 0.012, 0.012, 0.014 and 0.021 mag for B, V, g', r' and i' band, respectively. The equally high external accuracy of APASS photometry has been verified on secondary Landolt and Sloan photometric standard stars not involved in the APASS calibration process, and on a large body of literature data on field and cluster stars, confirming the absence of offsets and trends. Compared with the Carlsberg Meridian Catalog (CMC-15), APASS astrometry of RAVE stars is accurate to a median value of 0.098 arcsec. Brightness distribution functions for the RAVE stars have been derived in all bands. APASS photometry of RAVE stars, augmented by 2MASS JHK infrared data, has been chi2 fitted to a densely populated synthetic photometric library designed to widely explore in temperature, surface gravity, metallicity and reddening. Resulting Teff and E(B-V), computed over a range of options, are provided and discussed, and will be kept updated in response to future APASS and RAVE data releases. In the process it is found that the reddening caused by an homogeneous slab of dust, extending for 140 pc on either side of the Galactic plane and responsible for E(B-V,poles)=0.036 +/- 0.002 at the galactic poles, is a suitable approximation of the actual reddening encountered at Galactic latitudes |b|>=25 deg.

An improved estimator for the amplitude fnl of local-type non-Gaussianity from the cosmic microwave background (CMB) bispectrum is discussed. The standard estimator is constructed to be optimal in the zero-signal (i.e., Gaussian) limit. When applied to CMB maps which have a detectable level of non-Gaussianity the standard estimator is no longer optimal, possibly limiting the sensitivity of future observations to a non-Gaussian signal. Previous studies have proposed an improved estimator by using a realization-dependent normalization. Under the approximations of a flat sky and a vanishingly thin last-scattering surface, these studies showed that the variance of this improved estimator can be significantly smaller than the variance of the standard estimator when applied to non-Gaussian CMB maps. Here this technique is generalized to the full sky and to include the full radiation transfer function, yielding expressions for the improved estimator that can be directly applied to CMB maps. The ability of this estimator to reduce the variance as compared to the standard estimator in the face of a significant non-Gaussian signal is re-assessed using the full CMB transfer function. As a result of the late time integrated Sachs-Wolfe effect, the performance of the improved estimator is degraded. If CMB maps are first cleaned of the late-time ISW effect using a tracer of foreground structure, such as a galaxy survey or a measurement of CMB weak lensing, the new estimator does remove a majority of the excess variance, allowing a higher significance detection of fnl.

We examine model independent constraints on the high redshift and prerecombination expansion history from cosmic microwave background observations, using a combination of principal component analysis and other techniques. This can be translated to model independent limits on early dark energy and the number of relativistic species $N_{\rm eff}$. Models such as scaling (Doran-Robbers), dark radiation ($\Delta N_{\rm eff}$), and barotropic aether fall into distinct regions of eigenspace and can be easily distinguished from each other. Incoming CMB data will map the expansion history from $z=0$--$10^5$, achieving subpercent precision around recombination, and enable determination of the amount of early dark energy and valuable guidance to its nature.

Considerable recent attention has focussed on the prospects to use the cosmic microwave background (CMB) trispectrum to probe the physics of the early universe. Here we evaluate the probability distribution function (PDF) for the standard estimator tau_nle for the amplitude tau_nl of the CMB trispectrum both for the null-hypothesis (i.e., for Gaussian maps with tau_nl = 0) and for maps with a non-vanishing trispectrum (|tau_nl|>0). We find these PDFs to be highly non-Gaussian in both cases. We also evaluate the variance with which the trispectrum amplitude can be measured, <tau_nle^2>, as a function of its underlying value, tau_nl. We find a strong dependence of this variance on tau_nl. We also find that the variance does not, given the highly non-Gaussian nature of the PDF, effectively characterize the distribution. Detailed knowledge of these PDFs will therefore be imperative in order to properly interpret the implications of any given trispectrum measurement. For example, if a CMB experiment with a maximum multipole of lmax = 1500 (such as the Planck satellite) measures tau_nle = 0 then at the 95% confidence our calculations show that we can conclude tau_nl < 1005; assuming a Gaussian PDF but with the correct tau_nl-dependent variance we would incorrectly conclude tau_nl < 4225; further neglecting the tau_nl-dependence in the variance we would incorrectly conclude tau_nl < 361.

The cosmic microwave background (CMB) is affected by the total radiation density around the time of decoupling. At that epoch, neutrinos comprised a significant fraction of the radiative energy, but there could also be a contribution from primordial gravitational waves with frequencies greater than ~ 10^-15 Hz. If this cosmological gravitational wave background (CGWB) were produced under adiabatic initial conditions, its effects on the CMB and matter power spectrum would mimic massless non-interacting neutrinos. However, with homogenous initial conditions, as one might expect from certain models of inflation, pre big-bang models, phase transitions and other scenarios, the effect on the CMB would be distinct. We present updated observational bounds for both initial conditions using the latest CMB data at small scales from the South Pole Telescope (SPT) in combination with Wilkinson Microwave Anisotropy Probe (WMAP), current measurements of the baryon acoustic oscillations, and the Hubble parameter. With the inclusion of the data from SPT the adiabatic bound on the CGWB density is improved by a factor of 1.7 to 10^6 Omega_gw < 8.7 at the 95% confidence level (C.L.), with weak evidence in favor of an additional radiation component consistent with previous analyses. The constraint can be converted into an upper limit on the tension of horizon-sized cosmic strings that could generate this gravitational wave component, with Gmu < 2 10^-7 at 95% C.L., for string tension Gmu. The homogeneous bound improves by a factor of 3.5 to 10^6 Omega_gw < 1.0 at 95% C.L., with no evidence for such a component from current data.

Observations of the cosmic microwave background (CMB) and large-scale structure (LSS) provide a unique opportunity to explore the fundamental properties of the constituents that compose the cosmic dark radiation background (CDRB), of which the three standard neutrinos are thought to be the dominant component. We report on the first constraint to the CDRB rest-frame sound speed, ceff^2, using the most recent CMB and LSS data. Additionally, we report improved constraints to the CDRB viscosity parameter, cvis^2. For a non-interacting species, these parameters both equal 1/3. Using current data we find that a standard CDRB, composed entirely of three non-interacting neutrino species, is ruled out at the 99% confidence level (C.L.) with ceff^2 = 0.30 +0.027 -0.026 and cvis^2 = 0.44 +0.27 -0.21 (95% C.L.). We also discuss how constraints to these parameters from current and future observations (such as the Planck satellite) allow us to explore the fundamental properties of any anomalous radiative energy density beyond the standard three neutrinos.

We report the first detections of Blue Straggler Stars (BSS) in the bulge of the Milky Way galaxy. Proper motions from extensive space-based observations along a single sight-line allow us to separate a sufficiently clean and well-characterized bulge sample that we are able to detect a small population of bulge objects in the region of the color-magnitude diagram commonly occupied young objects and blue strgglers. However, variability measurements of these objects clearly establish that a fraction of them are blue stragglers. Out of the 42 objects found in this region of the color-magnitude diagram, we estimate that at least 18 are genuine BSS. We normalize the BSS population by our estimate of the number of horizontal branch stars in the bulge in order to compare the bulge to other stellar systems. The BSS fraction is clearly discrepant from that found in stellar clusters. The blue straggler population of dwarf spheroidals remains a subject of debate; some authors claim an anticorrelation between the normalised blue straggler fraction and integrated light. If this trend is real, then the bulge may extend it by three orders of magnitude in mass. Conversely, we find that the genuinely young (~5Gy or younger) population in the bulge, must be at most 3.4% under the most conservative scenario for the BSS population.

One of the principle efforts in cosmic microwave background (CMB) research is measurement of the parameter fnl that quantifies the departure from Gaussianity in a large class of non-minimal inflationary (and other) models. Estimators for fnl are composed of a sum of products of the temperatures in three different pixels in the CMB map. Since the number ~Npix^2 of terms in this sum exceeds the number Npix of measurements, these ~Npix^2 terms cannot be statistically independent. Therefore, the central-limit theorem does not necessarily apply, and the probability distribution function (PDF) for the fnl estimator does not necessarily approach a Gaussian distribution for N_pix >> 1. Although the variance of the estimators is known, the significance of a measurement of fnl depends on knowledge of the full shape of its PDF. Here we use Monte Carlo realizations of CMB maps to determine the PDF for two minimum-variance estimators: the standard estimator, constructed under the null hypothesis (fnl=0), and an improved estimator with a smaller variance for |fnl| > 0. While the PDF for the null-hypothesis estimator is very nearly Gaussian when the true value of fnl is zero, the PDF becomes significantly non-Gaussian when |fnl| > 0. In this case we find that the PDF for the null-hypothesis estimator fnl_hat is skewed, with a long non-Gaussian tail at fnl_hat > |fnl| and less probability at fnl_hat < |fnl| than in the Gaussian case. We provide an analytic fit to these PDFs. On the other hand, we find that the PDF for the improved estimator is nearly Gaussian for observationally allowed values of fnl. We discuss briefly the implications for trispectrum (and other higher-order correlation) estimators.

Using data from the COSMOS survey, we perform the first joint analysis of galaxy-galaxy weak lensing, galaxy spatial clustering, and galaxy number densities. Carefully accounting for sample variance and for scatter between stellar and halo mass, we model all three observables simultaneously using a novel and self-consistent theoretical framework. Our results provide strong constraints on the shape and redshift evolution of the stellar-to-halo mass relation (SHMR) from z=0.2 to z=1. At low stellar mass, we find that halo mass scales as Mh M*^0.46 and that this scaling does not evolve significantly with redshift to z=1. We show that the dark-to-stellar ratio, Mh/M*, varies from low to high masses, reaching a minimum of Mh/M*~27 at M*=4.5x10^10 Msun and Mh=1.2x10^12 Msun. This minimum is important for models of galaxy formation because it marks the mass at which the accumulated stellar growth of the central galaxy has been the most efficient. We describe the SHMR at this minimum in terms of the "pivot stellar mass", M*piv, the "pivot halo mass", Mhpiv, and the "pivot ratio", (Mh/M*)piv. Thanks to a homogeneous analysis of a single data set, we report the first detection of mass downsizing trends for both Mhpiv and M*piv. The pivot stellar mass decreases from M*piv=5.75+-0.13x10^10 Msun at z=0.88 to M*piv=3.55+-0.17x10^10 Msun at z=0.37. Intriguingly, however, the corresponding evolution of Mhpiv leaves the pivot ratio constant with redshift at (Mh/M*)piv~27. We use simple arguments to show how this result raises the possibility that star formation quenching may ultimately depend on Mh/M* and not simply Mh, as is commonly assumed. We show that simple models with such a dependence naturally lead to downsizing in the sites of star formation. Finally, we discuss the implications of our results in the context of popular quenching models, including disk instabilities and AGN feedback.

The Cosmic Origins Spectrograph (COS) was installed in the Hubble Space Telescope in May, 2009 as part of Servicing Mission 4 to provide high sensitivity, medium and low resolution spectroscopy at far- and near-ultraviolet wavelengths (FUV, NUV). COS is the most sensitive FUV/NUV spectrograph flown to date, spanning the wavelength range from 900{\AA} to 3200{\AA} with peak effective area approaching 3000 cm^2. This paper describes instrument design, the results of the Servicing Mission Orbital Verification (SMOV), and the ongoing performance monitoring program.

Minimum-variance estimators for the parameter fnl that quantifies local-model non-Gaussianity can be constructed from the cosmic microwave background (CMB) bispectrum (three-point function) and also from the trispectrum (four-point function). Some have suggested that a comparison between the estimates for the values of fnl from the bispectrum and trispectrum allow a consistency test for the model. But others argue that the saturation of the Cramer-Rao bound by the bispectrum estimator implies that no further information on fnl can be obtained from the trispectrum. Here we elaborate the nature of the correlation between the bispectrum and trispectrum estimators for fnl. We show that the two estimators become statistically independent in the limit of large number of CMB pixels and thus that the trispectrum estimator does indeed provide additional information on fnl beyond that obtained from the bispectrum. We explain how this conclusion is consistent with the Cramer-Rao bound. Our discussion of the Cramer-Rao bound may be of interest to those doing Fisher-matrix parameter-estimation forecasts or data analysis in other areas of physics as well.

The cosmic expansion history proceeds in broad terms from a radiation dominated epoch to matter domination to an accelerated, dark energy dominated epoch. We investigate whether intermittent periods of acceleration are possible in the early universe -- between Big Bang nucleosynthesis (BBN) and recombination and beyond. We establish that the standard picture is remarkably robust: observations of anisotropies in the cosmic microwave background exclude any extra period of accelerated expansion between 1 \leq z \lesssim 10^5 (corresponding to 5\times10^{-4}\ {\rm eV} \leq T \lesssim 25\ {\rm eV}).

Deviations from general relativity, such as could be responsible for the cosmic acceleration, would influence the growth of large scale structure and the deflection of light by that structure. We clarify the relations between several different model independent approaches to deviations from general relativity appearing in the literature, devising a translation table. We examine current constraints on such deviations, using weak gravitational lensing data of the CFHTLS and COSMOS surveys, cosmic microwave background radiation data of WMAP5, and supernova distance data of Union2. Markov Chain Monte Carlo likelihood analysis of the parameters over various redshift ranges yields consistency with general relativity at the 95% confidence level.

Modifications to GR generically predict time and scale-dependent effects which may be probed by observations of strong lensing by galaxies. Measurements of the stellar velocity dispersion determine the dynamical mass whereas measurements of the Einstein radius determine the lensing mass. In GR these two masses are equal; in alternative gravity theories they may not be. Using measurements of the stellar velocity dispersion and strong lensing around galaxies from the Sloan Lens ACS (SLACS) survey we place constraints on lensing in modified gravity theories and extend previous studies by applying this data to explore its dependence on various properties of the lens such as the lens redshift or mass and thereby constrain scalar-tensor, f(R) gravity theories, and generic parameterizations of deviations from GR. Besides applying the observations to these specific gravity theories, the data places a constraint on a generic dependence of modifications to GR on the lens mass and redshift. At the 68% confidence level we find that the ratio between the lensing and dynamical masses can only vary by less then 50% over a mass range for the lens galaxies of 1E12 < M/Msun < 1E14 and less than 40% over the redshift range 0.06 < z < 0.36.

Several anomalies have been identified which may imply a breakdown of the statistical isotropy of the cosmic microwave background (CMB). In particular, an anomalous alignment of the quadrupole and octopole and a hemispherical power asymmetry have increased in significance as the data have improved. There have been several attempts to explain these observations which explore isotropy breaking mechanisms within the early universe, but little attention has been given to the possibility that these anomalies have their origin within the local universe. We explore such a mechanism by considering the kinetic Sunyaev-Zel'dovich effect due to a gaseous halo associated with the Milky Way. Considering several physical models of an anisotropic free electron optical depth contributed by such a halo, we find that the associated screening maps of the primordial anisotropies have the necessary orientations to affect the anomaly statistics very significantly, but only if the column density of free electrons in the halo is at least an order of magnitude higher than indicated by current observations.

There is no direct evidence for radiation domination prior to big-bang nucleosynthesis, and so it is useful to consider how constraints to thermally-produced axions change in non-standard thermal histories. In the low-temperature-reheating scenario, radiation domination begins at temperatures as low as 1 MeV, and is preceded by significant entropy generation. Axion abundances are then suppressed, and cosmological limits to axions are significantly loosened. In a kination scenario, a more modest change to axion constraints occurs. Future possible constraints to axions and low-temperature reheating are discussed.

We summarize the utility of precise cosmic microwave background (CMB) polarization measurements as probes of the physics of inflation. We focus on the prospects for using CMB measurements to differentiate various inflationary mechanisms. In particular, a detection of primordial B-mode polarization would demonstrate that inflation occurred at a very high energy scale, and that the inflaton traversed a super-Planckian distance in field space. We explain how such a detection or constraint would illuminate aspects of physics at the Planck scale. Moreover, CMB measurements can constrain the scale-dependence and non-Gaussianity of the primordial fluctuations and limit the possibility of a significant isocurvature contribution. Each such limit provides crucial information on the underlying inflationary dynamics. Finally, we quantify these considerations by presenting forecasts for the sensitivities of a future satellite experiment to the inflationary parameters.

We present stellar proper motions in the Galactic bulge from the Sagittarius Window Eclipsing Extrasolar Search (SWEEPS) project using ACS/WFC on HST. Proper motions are extracted for more than 180,000 objects, with >81,000 measured to accuracy better than 0.3 mas/yr in both coordinates. We report several results based on these measurements: 1. Kinematic separation of bulge from disk allows a sample of >15,000 bulge objects to be extracted based on >6-sigma detections of proper motion, with <0.2% contamination from the disk. This includes the first detection of a candidate bulge Blue Straggler population. 2. Armed with a photometric distance modulus on a star by star basis, and using the large number of stars with high-quality proper motion measurements to overcome intrinsic scatter, we dissect the kinematic properties of the bulge as a function of distance along the line of sight. This allows us to extract the stellar circular speed curve from proper motions alone, which we compare with the circular speed curve obtained from radial velocities. 3. We trace the variation of the {l,b} velocity ellipse as a function of depth. 4. Finally, we use the density-weighted {l,b} proper motion ellipse produced from the tracer stars to assess the kinematic membership of the sixteen transiting planet candidates discovered in the Sagittarius Window; the kinematic distribution of the planet candidates is consistent with that of the disk and bulge stellar populations.

It is usually assumed that dark matter is produced during the radiation dominated era. There is, however, no direct evidence for radiation domination prior to big-bang nucleosynthesis. Two non-standard thermal histories are considered. In one, the low-temperature-reheating scenario, radiation domination begins as late as 1 MeV, and is preceded by significant entropy generation. Thermal axion relic abundances are then suppressed, and cosmological limits to axions are loosened. For reheating temperatures less than 35 MeV, the large-scale structure limit to the axion mass is lifted. The remaining constraint from the total density of matter is significantly relaxed. Constraints are also relaxed for higher reheating temperatures. In a kination scenario, a more modest change to cosmological axion constraints is obtained. Future possible constraints to axions and low-temperature reheating from the helium abundance and next-generation large-scale-structure surveys are discussed.

Inflation predicts a stochastic background of gravitational waves over a broad range of frequencies, from those accessible with cosmic microwave background (CMB) measurements, to those accessible directly with gravitational-wave detectors, like NASA's Big-Bang Observer (BBO), currently under study. In a previous paper [Phys. Rev. D73, 023504 (2006)] we connected CMB constraints to the amplitude and tensor spectral index of the inflationary gravitational-wave background (IGWB) at BBO frequencies for four classes of models of inflation by directly solving the inflationary equations of motion. Here we extend that analysis by including results obtained in the WMAP third-year data release as well as by considering two additional classes of inflationary models. As often noted in the literature, the recent indication that the primordial density power-spectrum has a red spectral index implies (with some caveats) that the amplitude of the IGWB may be large enough to be observable in the CMB polarization. Here we also explore the implications for the direct detection of the IGWB.

One of the possible low-energy consequences of string theory is the addition of a Chern-Simons term to the standard Einstein-Hilbert action of general relativity. It can be argued that the quintessence field should couple to this Chern-Simons term, and if so, it drives in the linearized theory a parity-violating interaction between the gravito-electric and gravitomagnetic fields. In this paper, the linearized spacetime for Chern-Simons gravity around a massive spinning body is found to include new modifications to the gravitomagnetic field that have not appeared in previous work. The orbits of test bodies and the precession of gyroscopes in this spacetime are calculated, leading to new constraints on the Chern-Simons parameter space due to current satellite experiments.

The exoplanets discovered so far have been mostly around relatively nearby and bright stars. As a result, the host stars are mostly (i) in the Galactic disk, (ii) relatively massive, and (iii) relatively metal rich. The aim of the SWEEPS project is to extend our knowledge to stars which (i) are in a different part of the Galaxy, (ii) have lower masses, and (iii) have a large range of metallicities. To achieve this goal, we used the Hubble Space Telescope to search for transiting planets around F, G, K, and M dwarfs in the Galactic bulge. We photometrically monitored 180,000 stars in a dense bulge field continuously for 7 days. We discovered 16 candidate transiting extrasolar planets with periods of 0.6 to 4.2 days, including a new class of ultra-short period planets (USPPs) with P < 1.2 days. Radial-velocity observations of the two brightest candidates support their planetary nature. These results suggest that planets are as abundant in the Galactic bulge as they are in the solar neighborhood, and they are equally abundant around low-mass stars (within a factor 2). The planet frequency increases with metallicity even for the stars in the Galactic bulge. All the USPP hosts are low-mass stars, suggesting either that close-in planets around higher-mass stars are irradiatively evaporated, or that the planets can migrate to close-in orbits only around such old and low-mass stars.

We report results of a recently-completed pre-Formulation Phase study of SPIRIT, a candidate NASA Origins Probe mission. SPIRIT is a spatial and spectral interferometer with an operating wavelength range 25 - 400 microns. SPIRIT will provide sub-arcsecond resolution images and spectra with resolution R = 3000 in a 1 arcmin field of view to accomplish three primary scientific objectives: (1) Learn how planetary systems form from protostellar disks, and how they acquire their inhomogeneous composition; (2) characterize the family of extrasolar planetary systems by imaging the structure in debris disks to understand how and where planets of different types form; and (3) learn how high-redshift galaxies formed and merged to form the present-day population of galaxies. Observations with SPIRIT will be complementary to those of the James Webb Space Telescope and the ground-based Atacama Large Millimeter Array. All three observatories could be operational contemporaneously.

It has been proposed that cosmic acceleration or inflation can be driven by replacing the Einstein-Hilbert action of general relativity with a function f(R) of the Ricci scalar R. Such f(R) gravity theories have been shown to be equivalent to scalar-tensor theories of gravity that are incompatible with Solar System tests of general relativity, as long as the scalar field propagates over Solar System scales. Specifically, the PPN parameter in the equivalent scalar-tensor theory is gamma=1/2, which is far outside the range allowed by observations. In response to a flurry of papers that questioned the equivalence of f(R) theory to scalar-tensor theories, it was recently shown explicitly, without resorting to the scalar-tensor equivalence, that the vacuum field equations for 1/R gravity around a spherically symmetric mass also yield gamma= 1/2. Here we generalize this analysis to f(R) gravity and enumerate the conditions that, when satisfied by the function f(R), lead to the prediction that gamma=1/2.

Shortly after the addition of a 1/R term to the Einstein-Hilbert action was proposed as a solution to the cosmic-acceleration puzzle, Chiba showed that such a theory violates Solar System tests of gravity. A flurry of recent papers have called Chiba's result into question. They argue that the spherically-symmetric vacuum spacetime in this theory is the Schwarzschild-de Sitter solution, making this theory consistent with Solar System tests. We point out that although the Schwarzschild-de Sitter solution exists in this theory, it is not the unique spherically-symmetric vacuum solution, and it is not the solution that describes the spacetime in the Solar System. The solution that correctly matches onto the stellar-interior solution differs from Schwarzschild-de Sitter in a way consistent with Chiba's claims. Thus, 1/R gravity is ruled out by Solar System tests.

More than 200 extrasolar planets have been discovered around relatively nearby stars, primarily through the Doppler line shifts owing to the reflex motions of their host stars, and more recently through transits of some planets across the face of the host stars. The detection of planets with the shortest known periods, 1.2 to 2.5 days, has mainly resulted from transit surveys which have generally targeted stars more massive than 0.75 M_sun. Here we report the results from a planetary transit search performed in a rich stellar field towards the Galactic bulge. We discovered 16 candidates with orbital periods between 0.4 and 4.2 days, five of which orbit stars of 0.44 to 0.75 M_sun. In two cases, radial-velocity measurements support the planetary nature of the companions. Five candidates have orbital periods below 1.0 day, constituting a new class of ultra-short-period planets (USPPs), which occur only around stars of less than 0.88 M_sun. This indicates that those orbiting very close to more luminous stars might be evaporatively destroyed, or that jovian planets around lower-mass stars might migrate to smaller radii.

The B-mode polarization lensing signal is a useful probe of the neutrino mass and to a lesser extent the dark energy equation of state as the signal depends on the integrated mass power spectrum between us and the last scattering surface. This lensing B-mode signal, however, is non-Gaussian and the resulting non-Gaussian covariance to the power spectrum cannot be ignored as correlations between B-mode bins are at a level of 0.1. For temperature and E-mode polarization power spectra, the non-Gasussian covariance is not significant, where we find correlations at the 10^{-5} level even for adjacent bins. The resulting degradation on neutrino mass and dark energy equation of state is about a factor of 2 to 3 when compared to the case where statistics are simply considered to be Gaussian. We also discuss parameter uncertainties achievable in upcoming experiments and show that at a given angular resolution for polarization observations, increasing the sensitivity beyond a certain noise value does not lead to an improved measurement of the neutrino mass and dark energy equation of state with B-mode power spectrum. For Planck, the resulting constraints on the sum of the neutrino masses is ~ 0.2 eV and on the dark energy equation of state parameter we find, sigma_w ~ 0.5.

A detection of the primordial gravitational wave background is considered to be the ``smoking-gun '' evidence for inflation. While super-horizon waves are probed with cosmic microwave background (CMB) polarization, the relic background will be studied with laser interferometers. The long lever arm spanned by the two techniques improves constraints on the inflationary potential and validation of consistency relations expected under inflation. If gravitational waves with a tensor-to-scalar amplitude ratio greater than 0.01 are detected by the CMB, then a direct detection experiment with a sensitivity consistent with current concept studies should be pursued vigorously. If no primordial tensors are detected by the CMB, a direct detection experiment to understand the simplest form of inflation must have a sensitivity improved by two to three orders of magnitude over current plans.

Primordial gravitational waves (GWs) with frequencies > 10^{-15} Hz contribute to the radiation density of the Universe at the time of decoupling of the cosmic microwave background (CMB). The effects of this GW background on the CMB and matter power spectra are identical to those due to massless neutrinos, unless the initial density-perturbation amplitude for the gravitational-wave gas is non-adiabatic, as may occur if such GWs are produced during inflation or some post-inflation phase transition. In either case, current observations provide a constraint to the GW amplitude that competes with that from big-bang nucleosynthesis (BBN), although it extends to much lower frequencies (~10^{-15} Hz rather than the ~10^{-10} Hz lower limit from BBN): at 95% confidence-level, Omega_gw h^2 < 6.9 x 10^{-6} for homogeneous (i.e., non-adiabatic) initial conditions. Future CMB experiments, like Planck and CMBPol, should allow sensitivities to Omega_gw h^2 < 1.4 x 10^{-6} and Omega_gw h^2 < 5 x 10^{-7}, respectively.

Inflation generically predicts a stochastic background of gravitational waves over a broad range of frequencies, from those accessible with cosmic microwave background (CMB) measurements, to those accessible directly with gravitational-wave detectors, like NASA's Big-Bang Observer (BBO) or Japan's Deci-Hertz Interferometer Gravitational-wave Observer (DECIGO), both currently under study. Here we investigate the detectability of the inflationary gravitational-wave background at BBO/DECIGO frequencies. To do so, we survey a range of slow-roll inflationary models consistent with constraints from the CMB and large-scale structure (LSS). We go beyond the usual assumption of power-law power spectra, which may break down given the 16 orders of magnitude in frequency between the CMB and direct detection, and solve instead the inflationary dynamics for four classes of inflaton potentials. Direct detection is possible in a variety of inflationary models, although probably not in any in which the gravitational-wave signal does not appear in the CMB polarization. However, direct detection by BBO/DECIGO can help discriminate between inflationary models that have the same slow-roll parameters at CMB/LSS scales.

We have studied a sample of 809 Mg II absorption systems with 1.0 < z_abs < 1.86 in the spectra of SDSS QSOs, with the aim of understanding the nature and abundance of the dust and the chemical abundances in the intervening absorbers. Normalized, composite spectra were derived, for abundance measurements, for the full sample and several sub-samples, chosen on the basis of the line strengths and other absorber and QSO properties. Average extinction curves were obtained for the sub-samples by comparing their geometric mean spectra with those of matching samples of QSOs without absorbers in their spectra. There is clear evidence for the presence of dust in the intervening absorbers. The 2175 A feature is not present in the extinction curves, for any of the sub-samples. The extinction curves are similar to the SMC extinction curve with a rising UV extinction below 2200 A. The absorber rest frame colour excess, E(B-V), derived from the extinction curves, depends on the absorber properties and ranges from < 0.001 to 0.085 for various sub-samples. The column densities of several ions do not show such a correspondingly large variation. The depletion pattern is similar to halo clouds in the Galaxy. Assuming an SMC gas-to-dust ratio we find a trend of increasing abundance with decreasing extinction; systems with N_H I \~ 10^{20} cm^{-2} show solar abundance of Zn. The large velocity spread of strong Mg II systems seems to be mimicked by weak lines of other elements. The ionization of the absorbers, in general appears to be low. QSOs with absorbers are, in general, at least three times as likely to have highly reddened spectra as compared to QSOs without any absorption systems in their spectra.

Using the Hubble Space Telescope, we have resolved individual red-giant branch stars in the halos of eight nearby spiral galaxies. The fields lie at projected distances between 2 and 13 kpc along the galaxies' minor axes. The data set allows a first look at the systematic trends in halo stellar populations. We have found that bright galaxies tend to have broad red-giant branch star color distributions with redder mean colors, suggesting that the heavy element abundance spread increases with the parent galaxy luminosity. The mean metallicity of the stellar halo, estimated using the mean colors of red-giant branch stars, correlates with the parent galaxy luminosity. The metallicity of the Milky Way halo falls nearly 1 dex below this luminosity-metallicity relation, suggesting that the halo of the Galaxy is more the exception than the rule for spiral galaxies; i.e., massive spirals with metal-poor halos are unusual. The luminosity-halo stellar abundance relation is consistent with the scaling relation expected for stellar systems embedded in dominant halos, suggesting that the bulk of the halo stellar population may have formed in situ.

(Abriged) We report results of a campaign to image the stellar populations in the halos of highly inclined spiral galaxies, with the fields roughly 10 kpc (projected) from the nuclei. We use the F814W (I) and F606W (V) filters in the Wide Field Planetary Camera 2, on board the Hubble Space telescope. Extended halo populations are detected in all galaxies. The color-magnitude diagrams appear to be completely dominated by giant-branch stars, with no evidence for the presence of young stellar populations in any of the fields. We find that the metallicity distribution functions are dominated by metal-rich populations, with a tail extending toward the metal poor end. To first order, the overall shapes of the metallicity distribution functions are similar to what is predicted by simple, single-component model of chemical evolution with the effective yields increasing with galaxy luminosity. However, metallicity distributions significantly narrower than the simple model are observed for a few of the most luminous galaxies in the sample. It appears clear that more luminous spiral galaxies also have more metal-rich stellar halos. The increasingly significant departures from the closed-box model for the more luminous galaxies indicate that a parameter in addition to a single yield is required to describe chemical evolution. This parameter, which could be related to gas infall or outflow either in situ or in progenitor dwarf galaxies that later merge to form the stellar halo, tends to act to make the metallicity distributions narrower at high metallicity.

We have imaged the halo populations of a sample of nearby spiral galaxies using the Wide Field Planetary Camera 2 on broad the Hubble Space Telescope with the aim of studying the stellar population properties and relating them to those of the host galaxies. In four galaxies, the red-giant branch is sufficiently well populated to measure the magnitude of the tip of the red-giant branch (TRGB), a well-known distance indicator. Using both the Sobel edge-detection technique and maximum-likelihood analysis to measure the $I$-band magnitude of the red giant branch tip, we determine distances to four nearby galaxies: NGC 253, NGC 4244, and NGC 4945, NGC 4258. For the first three galaxies, the TRGB distance determined here is more direct, and likely to be more accurate, than previous distance estimates. In the case of NGC 4258, our TRGB distance is in good agreement with the the geometrical maser distance, supporting the the Large Magellanic Cloud distance modulus $(m-M)_0 = 18.50$ that is generally adopted in recent estimates of the Hubble constant.

The carrier of the dust-associated photoluminescence process causing the extended red emission (ERE) in many dusty interstellar environments remains unidentified. Several competing models are more or less able to match the observed broad, unstructured ERE band. We now constrain the character of the ERE carrier further by determining the wavelengths of the radiation that initiates the ERE. Using the imaging capabilities of the Hubble Space Telescope, we have resolved the width of narrow ERE filaments appearing on the surfaces of externally illuminated molecular clouds in the bright reflection nebula NGC 7023 and compared them with the depth of penetration of radiation of known wavelengths into the same cloud surfaces. We identify photons with wavelengths shortward of 118 nm as the source of ERE initiation, not to be confused with ERE excitation, however. There are strong indications from the well-studied ERE in the Red Rectangle nebula and in the high-|b| Galactic cirrus that the photon flux with wavelengths shortward of 118 nm is too small to actually excite the observed ERE, even with 100% quantum efficiency. We conclude, therefore, that ERE excitation results from a two-step process. While none of the previously proposed ERE models can match these new constraints, we note that under interstellar conditions most polycyclic aromatic hydrocarbon (PAH) molecules are ionized to the di-cation stage by photons with E > 10.5 eV and that the electronic energy level structure of PAH di-cations is consistent with fluorescence in the wavelength band of the ERE. Therefore, PAH di-cations deserve further study as potential carriers of the ERE. (abridged)

The unidentified infrared bands (UIBs) have been observed in virtually every dusty astrophysical environment investigated. The UIB carrier must be abundant and ubiquitous. Strong evidence points to polycyclic aromatic hydrocarbons (PAHs) as likely candidates, but the identification is not complete. Additional diagnostics are needed to further constrain the UIB carrier, such as probing excitation sources which range from UV-strong to UV-weak, in order to determine the ``band gap'' of the UIB carrier. Observations and models suggest that the UIBs can be found in sources with weak UV fields. To that end, we present new results of observing the 3.3 micron spectral region in six stars embedded in reflection nebulae, and in six Vega-like stars. These objects have effective temperatures ranging from 3500 to 12,000 K. Their environments include dust that should be relatively unprocessed (reflection nebulae) and dust that has most likely undergone significant processing (Vega-like) by the embedded illumination source. Together with data from the literature, we have a sample of 27 sightlines. Our analysis suggests that neither the strength of the UV field impinging on the dust nor the effective temperature of the star is the determining factor in whether the 3.3 micron UIB emission is present in an object. We found three detections of the 3.3 micron-emission band, one in a Vega-type object, one in a Herbig Ae/Be object and one in a reflection nebula, and all with disks. The role of disk geometry is likely to be important in revealing or obscuring the photo-dissociation regions from which the UIB emission arises.

It is well known that the sightline toward HD 204827 in the cluster Trumpler 37 shows an UV extinction curve that does not follow CCM. However, when a dust component, foreground to the cluster, is removed, the residual extinction curve is identical to that found in the SMC within the uncertainties. The curve is very steep and has little or no 2175 A bump. The position of HD 204827 in the sky is projected onto the edge of the Cepheus IRAS Bubble. In addition, HD 204827 has an IRAS bowshock indicating that it may be embedded in dust swept up by the supernova that created the IRAS Bubble. Shocks due to the supernova may have led to substantial processing of this dust. The HD 204827 cloud is dense and rich in carbon molecules. The 3.4 $\mu$m feature indicating a C-H grain mantle is present in the dust toward HD 204827. The environment of the HD 204827 cloud dust may be similar to the dust associated with HD 62542 which lies on the edge of a stellar wind bubble and is also dense and rich in molecules. This sightline may be a Rosetta Stone if its environment can be related to those in the SMC having similar dust.

We have obtained new STIS/HST spectra to search for structure in the ultraviolet interstellar extinction curve, with particular emphasis on a search for absorption features produced by polycyclic aromatic hydrocarbons (PAHs). The presence of these molecules in the interstellar medium has been postulated to explain the infrared emission features seen in the 3-13 $\mu$m spectra of numerous sources. UV spectra are uniquely capable of identifying specific PAH molecules. We obtained high S/N UV spectra of stars which are significantly more reddened than those observed in previous studies. These data put limits on the role of small (30-50 carbon atoms) PAHs in UV extinction and call for further observations to probe the role of larger PAHs. PAHs are of importance because of their ubiquity and high abundance inferred from the infrared data and also because they may link the molecular and dust phases of the interstellar medium. A presence or absence of ultraviolet absorption bands due to PAHs could be a definitive test of this hypothesis. We should be able to detect a 20 \AA wide feature down to a 3$\sigma$ limit of $\sim$0.02 A$_V$. No such absorption features are seen other than the well-known 2175 \AA bump.

As part of the Hubble Deep Field South program, a set of shorter 2-orbit observations were obtained of the area adjacent to the deep fields. The WFPC2 flanking fields cover a contiguous solid angle of 48 square arcminutes. Parallel observations with the STIS and NICMOS instruments produce a patchwork of additional fields with optical and near-infrared (1.6 micron) response. Deeper parallel exposures with WFPC2 and NICMOS were obtained when STIS observed the NICMOS deep field. These deeper fields are offset from the rest, and an extended low surface brightness object is visible in the deeper WFPC2 flanking field. In this data paper, which serves as an archival record of the project, we discuss the observations and data reduction, and present SExtractor source catalogs and number counts derived from the data. Number counts are broadly consistent with previous surveys from both ground and space. Among other things, these flanking field observations are useful for defining slit masks for spectroscopic follow-up over a wider area around the deep fields, for studying large-scale structure that extends beyond the deep fields, for future supernova searches, and for number counts and morphological studies, but their ultimate utility will be defined by the astronomical community.

Space Telescope Imaging Spectrograph (STIS) and archival International Ultraviolet Explorer (IUE) data were used to search for structure in the ultraviolet interstellar extinction curve. Between 1250 and 3200 A, in addition to the well-known 2175 A bump, there is only one fairly small (10%) absorption feature at 2700 A. The feature is very broad (>400 A). It coincides extremely well in position and shape with laboratory measurements of presolar diamond grains found in meteorites. As this feature appears in both the new STIS and archival IUE data, it is unlikely to be an instrumental effect. This feature appears to be ubiquitous in the ISM as it is seen along sightlines toward various clouds with reddenings in the range of E(B-V) = 0.35 to 1.73 mag. The identification with diamonds remains tentative since it is based on a single, weak feature. If the feature is due to interstellar diamonds then they require an order of magnitude more absorption strength per gram/cm^2 than the presolar meteoritic diamonds measured in the lab.

Interstellar dust contains a component which reveals its presence by emitting a broad, unstructured band of light in the 540 to 950 nm wavelength range, referred to as Extended Red Emission (ERE). The presence of interstellar dust and ultraviolet photons are two necessary conditions for ERE to occur. This is the basis for suggestions which attribute ERE to an interstellar dust component capable of photoluminescence. In this study, we have collected all published ERE observations with absolute-calibrated spectra for interstellar environments, where the density of ultraviolet photons can be estimated reliably. In each case, we determined the band-integrated ERE intensity, the wavelength of peak emission in the ERE band, and the efficiency with which absorbed ultraviolet photons are contributing to the ERE. The data show that radiation is not only driving the ERE, as expected for a photoluminescence process, but is modifying the ERE carrier as manifested by a systematic increase in the ERE band's peak wavelength and a general decrease in the photon conversion efficiency with increasing densities of the prevailing exciting radiation. The overall spectral characteristics of the ERE and the observed high quantum efficiency of the ERE process are currently best matched by the recently proposed silicon nanoparticle (SNP) model. Using the experimentally established fact that ionization of semiconductor nanoparticles quenches their photoluminescence, we proceeded to test the SNP model by developing a quantitative model for the excitation and ionization equilibrium of SNPs under interstellar conditions for a wide range of radiation field densities.

We have made high-resolution, high-sensitivity dynamic spectra of a sample of strong pulsars at 430 MHz with the Arecibo radiotelescope. For 4 pulsars we find faint but sharply delineated features in the secondary spectra. These are examples of the previously observed ``crisscross'' or ``multiple drift slope'' phenomenon presumed to be due to multiple imaging of the pulsar by the interstellar medium. The unprecedented resolution and dynamic range of our observations allow a deeper level of analysis. Distances to the dominant scattering screen along the line of sight are determined and are shown to agree well with those inferred from other scintillation phenomena. Multiple imaging of the pulsar by the ISM is required. A compact central image surrounded by a faint scattering halo, roughly circularly symmetric, is consistent with the data. Scattering from filaments may also be consistent. The angular extent of the scattering material parallel to the direction of the pulsar velocity is roughly 5 mas, corresponding to a linear extent of about 2 AU. Further observations of these features should allow better discrimination between models and an identification of the scattering structures.

The Hubble Deep Field-South observations targeted a high-galactic-latitude field near QSO J2233-606. We present WFPC2 observations of the field in four wide bandpasses centered at roughly 300, 450, 606, and 814 nm. Observations, data reduction procedures, and noise properties of the final images are discussed in detail. A catalog of sources is presented, and the number counts and color distributions of the galaxies are compared to a new catalog of the HDF-N that has been constructed in an identical manner. The two fields are qualitatively similar, with the galaxy number counts for the two fields agreeing to within 20%. The HDF-S has more candidate Lyman-break galaxies at z > 2 than the HDF-N. The star-formation rate per unit volume computed from the HDF-S, based on the UV luminosity of high-redshift candidates, is a factor of 1.9 higher than from the HDF-N at z ~ 2.7, and a factor of 1.3 higher at z ~ 4.

We present 0.35 to 2.5 micron spectra of the south and northwest filaments in the reflection nebula NGC 7023. These spectra were used to test the theory of Seahra & Duley that carbon nanoparticles are responsible for Extended Red Emission (ERE). Our spectra fail to show their predicted second emission band at 1.0 micron even though both filaments exhibit strong emission in the familiar 0.7 micron ERE band. The northwest filament spectrum does show one, and possibly two, new dust emission features in the near-infrared. We clearly detect a strong emission band at 1.5 micron which we tentatively attribute to beta-FeSi_2 grains. We tentatively detect a weaker emission band at 1.15 micron which coincides with the location expected for transitions from the conduction band to mid-gap defect states of silicon nanoparticles. This is added evidence that silicon nanoparticles are responsible for ERE as they already can explain the observed behavior of the main visible ERE band.

We present the imaging observations made with the Space Telescope Imaging Spectrograph of the Hubble Deep Field - South. The field was imaged in 4 bandpasses: a clear CCD bandpass for 156 ksec, a long-pass filter for 22-25 ksec per pixel typical exposure, a near-UV bandpass for 23 ksec, and a far-UV bandpass for 52 ksec. The clear visible image is the deepest observation ever made in the UV-optical wavelength region, reaching a 10 sigma AB magnitude of 29.4 for an object of area 0.2 square arcseconds. The field contains QSO J2233-606, the target of the STIS spectroscopy, and extends 50"x50" for the visible images, and 25"x25" for the ultraviolet images. We present the images, catalog of objects, and galaxy counts obtained in the field.

We have obtained V and I images of the lone globular cluster that belongs to the dwarf Local Group irregular galaxy known as WLM. The color-magnitude diagram of the cluster shows that it is a normal old globular cluster with a well-defined giant branch reaching to M_V=-2.5, a horizontal branch at M_V=+0.5, and a sub-giant branch extending to our photometry limit of M_V=+2.0. A best fit to theoretical isochrones indicates that this cluster has a metallicity of [Fe/H]=-1.52\pm0.08 and an age of 14.8\pm0.6 Gyr, thus indicating that it is similar to normal old halo globulars in our Galaxy. From the fit we also find that the distance modulus of the cluster is 24.73\pm0.07 and the extinction is A_V=0.07\pm0.06, both values that agree within the errors with data obtained for the galaxy itself by others. We conclude that this normal massive cluster was able to form during the formation of WLM, despite the parent galaxy's very small intrinsic mass and size.

To examine a recently proposed hypothesis that silicon nanoparticles are the source of extended red emission (ERE) in the interstellar medium, we performed a detailed modeling of the mean Galactic extinction in the presence of silicon nanoparticles. For this goal we used the appropriate optical constants of nanosized Si, essentially different from those of bulk Si due to quantum confinement. It was found that a dust mixture of silicon nanoparticles, bare graphite grains, silicate core-organic refractory mantle grains and three-layer silicate-water ice-organic refractory grains works well in explaining the extinction and, in addition, results in the acceptable fractions of UV/visible photons absorbed by silicon nanoparticles: 0.071-0.081. Since these fractions barely agree with the fraction of UV/visible photons needed to excite the observed ERE, we conclude that the intrinsic photon conversion efficiency of the photoluminescence by silicon nanoparticles must be near 100%, if they are the source of the ERE.

We have used spectra taken between 1992 and 1997 to derive the spectroscopic orbit of the eclipsing double-lined spectroscopic binary HD 197770. This binary has a period of 99.69 +/- 0.02 days and K amplitudes of 31.2 +/- 0.8 and 47.1 +/- 0.4 km s^{-1} for components A & B, respectively. The msin^{3}i values for A & B are 2.9 and 1.9, respectively, and are close to the actual masses due to the eclipsing nature of this binary. Both components of HD 197770 have spectral types near B2 III. This means both components are undermassive by about a factor of five and, thus, evolved stars. Additional evidence of the evolved nature of HD 197770 is found in 25, 60, and 100 micron IRAS images of HD 197770. These images show 2 apparent shells centered on HD 197770; a bright 60 micron shell with a 14' diameter and a larger (1.2 degeree diameter) bubble-like feature. At least one of the components of HD 197770 is likely to be a post-AGB star.