We propose a non-parametric method to denoise 1D stellar spectra based on wavelet shrinkage followed by adaptive Kalman thresholding. Wavelet shrinkage denoising involves applying the Discrete Wavelet Transform (DWT) to the input signal, `shrinking' certain frequency components in the transform domain, and then applying inverse DWT to the reduced components. The performance of this procedure is influenced by the choice of base wavelet, the number of decomposition levels, and the thresholding function. Typically, these parameters are chosen by `trial and error', which can be strongly dependent on the properties of the data being denoised. We here introduce an adaptive Kalman-filter-based thresholding method that eliminates the need for choosing the number of decomposition levels. We use the `Haar' wavelet basis, which we found to be the best-suited for 1D stellar spectra. We introduce various levels of Poisson noise into synthetic PHOENIX spectra, and test the performance of several common denoising methods against our own. It proves superior in terms of noise suppression and peak shape preservation. We expect it may also be of use in automatically and accurately filtering low signal-to-noise galaxy and quasar spectra obtained from surveys such as SDSS, Gaia, LSST, PESSTO, VANDELS, LEGA-C, and DESI.

A common problem in cosmology is to integrate the product of two or more spherical Bessel functions (sBFs) with different configuration-space arguments against the power spectrum or its square, weighted by powers of wavenumber. Naively computing them scales as $N_{\rm g}^{p+1}$ with $p$ the number of configuration space arguments and $N_{\rm g}$ the grid size, and they cannot be done with Fast Fourier Transforms (FFTs). Here we show that by rewriting the sBFs as sums of products of sine and cosine and then using the product to sum identities, these integrals can then be performed using 1-D FFTs with $N_{\rm g} \log N_{\rm g}$ scaling. This "rotation" method has the potential to accelerate significantly a number of calculations in cosmology, such as perturbation theory predictions of loop integrals, higher order correlation functions, and analytic templates for correlation function covariance matrices. We implement this approach numerically both in a free-standing, publicly-available \textsc{Python} code and within the larger, publicly-available package \texttt{mcfit}. The rotation method evaluated with direct integrations already offers a factor of 6-10$\times$ speed-up over the naive approach in our test cases. Using FFTs, which the rotation method enables, then further improves this to a speed-up of $\sim$$1000-3000\times$ over the naive approach. The rotation method should be useful in light of upcoming large datasets such as DESI or LSST. In analysing these datasets recomputation of these integrals a substantial number of times, for instance to update perturbation theory predictions or covariance matrices as the input linear power spectrum is changed, will be one piece in a Monte Carlo Markov Chain cosmological parameter search: thus the overall savings from our method should be significant.

The nonlinear mass is a characteristic scale in halo formation that has wide-ranging applications across cosmology. Naively, computing it requires repeated numerical integration to calculate the variance of the power spectrum on different scales and determine which scales exceed the threshold for nonlinear collapse. We accelerate this calculation by working in configuration space and approximating the correlation function as a polynomial at $r <= 5$ $h^{-1}$ Mpc. This enables an analytic rather than numerical solution, accurate across a variety of cosmologies to 0.1$-$1% (depending on redshift) and 10$-$20 times faster than the naive numerical method. We also present a further acceleration (40$-$80 times faster than the naive method) in which we determine the polynomial coefficients using a Taylor expansion in the cosmological parameters rather than re-fitting a polynomial to the correlation function. Our acceleration greatly reduces the cost of repeated calculation of the nonlinear mass. This will be useful for MCMC analyses to constrain cosmological parameters from the highly nonlinear regime, e.g. with data from upcoming surveys.

MegaMapper is a proposed ground-based experiment to measure Inflation parameters and Dark Energy from galaxy redshifts at 2<z<5. A 6.5-m Magellan telescope will be coupled with DESI spectrographs to achieve multiplexing of 20,000. MegaMapper would be located at Las Campanas Observatory to fully access LSST imaging for target selection.

Support for early career astronomers who are just beginning to explore astronomy research is imperative to increase retention of diverse practitioners in the field. Since 2010, Astrobites has played an instrumental role in engaging members of the community -- particularly undergraduate and graduate students -- in research. In this white paper, the Astrobites collaboration outlines our multi-faceted online education platform that both eases the transition into astronomy research and promotes inclusive professional development opportunities. We additionally offer recommendations for how the astronomy community can reduce barriers to entry to astronomy research in the coming decade.

We propose an experiment, the Cosmic Accelerometer, designed to yield velocity precision of $\leq 1$ cm/s with measurement stability over years to decades. The first-phase Cosmic Accelerometer, which is at the scale of the Astro2020 Small programs, will be ideal for precision radial velocity measurements of terrestrial exoplanets in the Habitable Zone of Sun-like stars. At the same time, this experiment will serve as the technical pathfinder and facility core for a second-phase larger facility at the Medium scale, which can provide a significant detection of cosmological redshift drift on a 6-year timescale. This larger facility will naturally provide further detection/study of Earth twin planet systems as part of its external calibration process. This experiment is fundamentally enabled by a novel low-cost telescope technology called PolyOculus, which harnesses recent advances in commercial off the shelf equipment (telescopes, CCD cameras, and control computers) combined with a novel optical architecture to produce telescope collecting areas equivalent to standard telescopes with large mirror diameters. Combining a PolyOculus array with an actively-stabilized high-precision radial velocity spectrograph provides a unique facility with novel calibration features to achieve the performance requirements for the Cosmic Accelerometer.

Nearly a century after the discovery that we live in an expanding Universe, and two decades after the discovery of accelerating cosmic expansion, there remains no direct detection of this acceleration via redshift drift - a change in the cosmological expansion velocity versus time. Because cosmological redshift drift directly determines the Hubble parameter H(z), it is arguably the cleanest possible measurement of the expansion history, and has the potential to constrain dark energy models (e.g. Kim et al. 2015). The challenge is that the signal is small - the best observational constraint presently has an uncertainty several orders of magnitude larger than the expected signal (Darling 2012). Nonetheless, direct detection of redshift drift is becoming feasible, with upcoming facilities such as the ESO-ELT and SKA projecting possible detection within two to three decades. This timescale is uncomfortably long given the potential of this cosmological test. With dedicated experiments it should be possible to rapidly accelerate progress and detect redshift drift with only a five-year observational baseline. Such a facility would also be ideal for precision radial velocity measurements of exoplanets, which could be obtained as a byproduct of the ongoing calibration measurements for the experiment.

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space.

The DESI Legacy Imaging Surveys are a combination of three public projects (the Dark Energy Camera Legacy Survey, the Beijing-Arizona Sky Survey, and the Mayall z-band Legacy Survey) that will jointly image approximately 14,000 deg^2 of the extragalactic sky visible from the northern hemisphere in three optical bands (g, r, and z) using telescopes at the Kitt Peak National Observatory and the Cerro Tololo Inter-American Observatory. The combined survey footprint is split into two contiguous areas by the Galactic plane. The optical imaging is conducted using a unique strategy of dynamically adjusting the exposure times and pointing selection during observing that results in a survey of nearly uniform depth. In addition to calibrated images, the project is delivering a catalog, constructed by using a probabilistic inference-based approach to estimate source shapes and brightnesses. The catalog includes photometry from the grz optical bands and from four mid-infrared bands (at 3.4, 4.6, 12 and 22 micorons) observed by the Wide-field Infrared Survey Explorer (WISE) satellite during its full operational lifetime. The project plans two public data releases each year. All the software used to generate the catalogs is also released with the data. This paper provides an overview of the Legacy Surveys project.

The galaxy bispectrum, measuring excess clustering of galaxy triplets, offers a probe of dark energy via baryon acoustic oscillations (BAOs). However up to now it has been severely underused due to the combinatorically explosive number of triangles. Here we exploit interference in the bispectrum to identify triangles that amplify BAOs. This approach reduces the computational cost of estimating covariance matrices, offers an improvement in BAO constraints equivalent to lengthening BOSS by 30%, and simplifies adding bispectrum BAO information to future large-scale redshift survey analyses.

Perturbation theory (PT) is often used to model statistical observables capturing the translation and rotation-invariant information in cosmological density fields. PT produces higher-order corrections by integration over linear statistics of the density fields weighted by kernels resulting from recursive solution of the fluid equations. These integrals quickly become high-dimensional and naively require increasing computational resources the higher the order of the corrections. Here we show how to decouple the integrands that often produce this issue, enabling PT corrections to be computed as a sum of products of independent 1-D integrals. Our approach is related to a commonly used method for calculating multi-loop Feynman integrals in Quantum Field Theory, the Gegenbauer Polynomial $x$-Space Technique (GPxT). We explicitly reduce the three terms entering the 2-loop power spectrum, formally requiring 9-D integrations, to sums over successive 1-D radial integrals. These 1-D integrals can further be performed as convolutions, rendering the scaling of this method $N_{\rm g} \log N_{\rm g}$ with $N_{\rm g}$ the number of grid points used for each Fast Fourier Transform. This method should be highly enabling for upcoming large-scale structure redshift surveys where model predictions at an enormous number of cosmological parameter combinations will be required by Monte Carlo Markov Chain searches for the best-fit values.

We present a picture of the matter bispectrum in a novel "interferometric" basis designed to highlight interference of the baryon acoustic oscillations (BAO) in the power spectra composing it. Triangles where constructive interference amplifies BAO provide stronger cosmic distance constraints than triangles with destructive interference. We show that the amplitude of the BAO feature in the full cyclically summed bispectrum can be decomposed into simpler contributions from single terms or pairs of terms in the perturbation theory bispectrum, and that across large swathes of our parameter space the full BAO amplitude is described well by the amplitude of BAO in a single term. The dominant term is determined largely by the $F^{(2)}$ kernel of Eulerian standard perturbation theory. We present a simple physical picture of the BAO amplitude in each term; the BAO signal is strongest in triangle configurations where two wavenumbers differ by a multiple of the BAO fundamental wavelength.

We present the first application of the angle-dependent 3-Point Correlation Function (3PCF) to the density fields magnetohydrodynamic (MHD) turbulence simulations intended to model interstellar (ISM) turbulence. Previous work has demonstrated that the angle-averaged bispectrum, the 3PCF's Fourier-space analog, is sensitive to the sonic and Alfv\'enic Mach numbers of turbulence. Here we show that introducing angular information via multipole moments with respect to the triangle opening angle offers considerable additional discriminatory power on these parameters. We exploit a fast, order $N_{\rm g} \log N_{\rm g}$ ($N_{\rm g}$ the number of grid cells used for a Fourier Transform) 3PCF algorithm to study a suite of MHD turbulence simulations with 10 different combinations of sonic and Alfv\'enic Mach numbers over a range from sub to super-sonic and sub to super-Alfv\'{e}nic. The 3PCF algorithm's speed for the first time enables full quantification of the time-variation of our signal: we study 9 timeslices for each condition, demonstrating that the 3PCF is sufficiently time-stable to be used as an ISM diagnostic. In future, applying this framework to 3-D dust maps will enable better treatment of dust as a cosmological foreground as well as reveal conditions in the ISM that shape star formation.

We present nbodykit, an open-source, massively parallel Python toolkit for analyzing large-scale structure (LSS) data. Using Python bindings of the Message Passing Interface (MPI), we provide parallel implementations of many commonly used algorithms in LSS. nbodykit is both an interactive and scalable piece of scientific software, performing well in a supercomputing environment while still taking advantage of the interactive tools provided by the Python ecosystem. Existing functionality includes estimators of the power spectrum, 2 and 3-point correlation functions, a Friends-of-Friends grouping algorithm, mock catalog creation via the halo occupation distribution technique, and approximate N-body simulations via the FastPM scheme. The package also provides a set of distributed data containers, insulated from the algorithms themselves, that enable nbodykit to provide a unified treatment of both simulation and observational data sets. nbodykit can be easily deployed in a high performance computing environment, overcoming some of the traditional difficulties of using Python on supercomputers. We provide performance benchmarks illustrating the scalability of the software. The modular, component-based approach of nbodykit allows researchers to easily build complex applications using its tools. The package is extensively documented at http://nbodykit.readthedocs.io, which also includes an interactive set of example recipes for new users to explore. As open-source software, we hope nbodykit provides a common framework for the community to use and develop in confronting the analysis challenges of future LSS surveys.

We obtain novel closed form solutions to the Friedmann equation for cosmological models containing a component whose equation of state is that of radiation $(w=1/3)$ at early times and that of cold pressureless matter $(w=0)$ at late times. The equation of state smoothly transitions from the early to late-time behavior and exactly describes the evolution of a species with a Dirac Delta function distribution in momentum magnitudes $|\vec{p}_0|$ (i.e. all particles have the same $|\vec{p}_0|$). Such a component, here termed "hot matter", is an approximate model for both neutrinos and warm dark matter. We consider it alone and in combination with cold matter and with radiation, also obtaining closed-form solutions for the growth of super-horizon perturbations in each case. The idealized model recovers $t(a)$ to better than $1.5\%$ accuracy for all $a$ relative to a Fermi-Dirac distribution (as describes neutrinos). We conclude by adding the second moment of the distribution to our exact solution and then generalizing to include all moments of an arbitrary momentum distribution in a closed form solution.

We present an algorithm enabling computation of the anisotropic redshift-space galaxy 3-point correlation function (3PCF) scaling as $N^2$, with $N$ the number of galaxies. Our previous work showed how to compute the isotropic 3PCF with this scaling by expanding the radially-binned density field around each galaxy in the survey into spherical harmonics and combining these coefficients to form multipole moments. The $N^2$ scaling occurred because this approach never explicitly required the relative angle between a galaxy pair about the primary galaxy. Here we generalize this work, demonstrating that in the presence of azimuthally-symmetric anisotropy produced by redshift-space distortions (RSD) the 3PCF can be described by two triangle side lengths, two independent total angular momenta, and a spin. This basis for the anisotropic 3PCF allows its computation with negligible additional work over the isotropic 3PCF. We also present the covariance matrix of the anisotropic 3PCF measured in this basis. Our algorithm tracks the full 5-D redshift-space 3PCF, uses an accurate line of sight to each triplet, is exact in angle, and easily handles edge correction. It will enable use of the anisotropic large-scale 3PCF as a probe of RSD in current and upcoming large-scale redshift surveys.

The nature of dark energy and the complete theory of gravity are two central questions currently facing cosmology. A vital tool for addressing them is the 3-point correlation function (3PCF), which probes deviations from a spatially random distribution of galaxies. However, the 3PCF's formidable computational expense has prevented its application to astronomical surveys comprising millions to billions of galaxies. We present Galactos, a high-performance implementation of a novel, O(N^2) algorithm that uses a load-balanced k-d tree and spherical harmonic expansions to compute the anisotropic 3PCF. Our implementation is optimized for the Intel Xeon Phi architecture, exploiting SIMD parallelism, instruction and thread concurrency, and significant L1 and L2 cache reuse, reaching 39% of peak performance on a single node. Galactos scales to the full Cori system, achieving 9.8PF (peak) and 5.06PF (sustained) across 9636 nodes, making the 3PCF easily computable for all galaxies in the observable universe.

Measurements of line-of-sight dependent clustering via the galaxy power spectrum's multipole moments constitute a powerful tool for testing theoretical models in large-scale structure. Recent work shows that this measurement, including a moving line-of-sight, can be accelerated using Fast Fourier Transforms (FFTs) by decomposing the Legendre polynomials into products of Cartesian vectors. Here, we present a faster, optimal means of using FFTs for this measurement. We avoid redundancy present in the Cartesian decomposition by using a spherical harmonic decomposition of the Legendre polynomials. Consequently, our method is substantially faster: a given multipole of order $\ell$ requires only $2\ell+1$ FFTs rather than the $(\ell+1)(\ell+2)/2$ FFTs of the Cartesian approach. For the hexadecapole ($\ell = 4$), this translates to $40\%$ fewer FFTs, with increased savings for higher $\ell$. The reduction in wall-clock time enables the calculation of finely-binned wedges in $P(k,\mu)$, obtained by computing multipoles up to a large $\ell_{\rm max}$ and combining them. This transformation has a number of advantages. We demonstrate that by using non-uniform bins in $\mu$, we can isolate plane-of-sky (angular) systematics to a narrow bin at $\mu \simeq 0$ while eliminating the contamination from all other bins. We also show that the covariance matrix of clustering wedges binned uniformly in $\mu$ becomes ill-conditioned when combining multipoles up to large values of $\ell_{\rm max}$, but that the problem can be avoided with non-uniform binning. As an example, we present results using $\ell_{\rm max}=16$, for which our procedure requires a factor of 3.4 fewer FFTs than the Cartesian method, while removing the first $\mu$ bin leads only to a 7% increase in statistical error on $f \sigma_8$, as compared to a 54% increase with $\ell_{\rm max}=4$.

We present the large-scale 3-point correlation function (3PCF) of the SDSS DR12 CMASS sample of $777,202$ Luminous Red Galaxies, the largest-ever sample used for a 3PCF or bispectrum measurement. We make the first high-significance ($4.5\sigma$) detection of Baryon Acoustic Oscillations (BAO) in the 3PCF. Using these acoustic features in the 3PCF as a standard ruler, we measure the distance to $z=0.57$ to $1.7\%$ precision (statistical plus systematic). We find $D_{\rm V}= 2024\pm29\;{\rm Mpc\;(stat)}\pm20\;{\rm Mpc\;(sys)}$ for our fiducial cosmology (consistent with Planck 2015) and bias model. This measurement extends the use of the BAO technique from the 2-point correlation function (2PCF) and power spectrum to the 3PCF and opens an avenue for deriving additional cosmological distance information from future large-scale structure redshift surveys such as DESI. Our measured distance scale from the 3PCF is fairly independent from that derived from the pre-reconstruction 2PCF and is equivalent to increasing the length of BOSS by roughly 10\%; reconstruction appears to lower the independence of the distance measurements. Fitting a model including tidal tensor bias yields a moderate significance ($2.6\sigma)$ detection of this bias with a value in agreement with the prediction from local Lagrangian biasing.

We search for a galaxy clustering bias due to a modulation of galaxy number with the baryon-dark matter relative velocity resulting from recombination-era physics. We find no detected signal and place the constraint $b_v < 0.01$ on the relative velocity bias for the CMASS galaxies. This bias is an important potential systematic of Baryon Acoustic Oscillation (BAO) method measurements of the cosmic distance scale using the 2-point clustering. Our limit on the relative velocity bias indicates a systematic shift of no more than $0.3\%$ rms in the distance scale inferred from the BAO feature in the BOSS 2-point clustering, well below the $1\%$ statistical error of this measurement. This constraint is the most stringent currently available and has important implications for the ability of upcoming large-scale structure surveys such as DESI to self-protect against the relative velocity as a possible systematic.

We present a configuration-space model of the large-scale galaxy 3-point correlation function (3PCF) based on leading-order perturbation theory and including redshift space distortions (RSD). This model should be useful in extracting distance-scale information from the 3PCF via the Baryon Acoustic Oscillation (BAO) method. We include the first redshift-space treatment of biasing by the baryon-dark matter relative velocity. Overall, on large scales the effect of RSD is primarily a renormalization of the 3PCF that is roughly independent of both physical scale and triangle opening angle; for our adopted $\Omega_{\rm m}$ and bias values, the rescaling is a factor of $\sim 1.8$. We also present an efficient scheme for computing 3PCF predictions from our model, important for allowing fast exploration of the space of cosmological parameters in future analyses.

We report a measurement of the large-scale 3-point correlation function of galaxies using the largest dataset for this purpose to date, 777, 202 Luminous Red Galaxies in the Sloan Digital Sky Survey Baryon Acoustic Oscillation Spectroscopic Survey (SDSS BOSS) DR12 CMASS sample. This work exploits the novel algorithm of Slepian & Eisenstein (2015b) to compute the multipole moments of the 3PCF in $\mathcal{O}(N^2)$ time, with $N$ the number of galaxies. Leading-order perturbation theory models the data well in a compressed basis where one triangle side is integrated out. We also present an accurate and computationally efficient means of estimating the covariance matrix. With these techniques the redshift-space linear and non-linear bias are measured, with 2.6% precision on the former if $\sigma_8$ is fixed. The data also indicates a $2.8\sigma$ preference for the BAO, confirming the presence of BAO in the 3-point function.

The anisotropic 2-point correlation function (2PCF) of galaxies measures pairwise clustering as a function of the pair separation's angle to the line of sight. The latter is often defined as either the angle bisector of the observer-galaxy-pair triangle or the vector from the observer to the separation midpoint. Here we show how to accelerate either of these measurements with Fourier Transforms, using a slight generalization of the Yamamoto et al. (2006) estimator in which each member of the pair is used successively as the line of sight. We also present perturbation theory predictions for our generalized estimator including wide-angle corrections.

In linear perturbation theory, all information about the growth of structure is contained in the Green's function, or equivalently, transfer function. These functions are generally computed using numerical codes or by phenomenological fitting formula anchored in accurate analytic results in the limits of large and small scale. Here we present a framework for analytically solving all scales, in particular the intermediate scales relevant for the baryon acoustic oscillations (BAO). We solve for the Green's function and transfer function using spherically-averaged overdensities and the approximation that the density of the coupled baryon-photon fluid is constant interior to the sound horizon.

Though Fourier Transforms (FTs) are a common technique for finding correlation functions, they are not typically used in computations of the anisotropy of the two-point correlation function (2PCF) about the line of sight in wide-angle surveys because the line-of-sight direction is not constant on the Cartesian grid. Here we show how FTs can be used to compute the multipole moments of the anisotropic 2PCF. We also show how FTs can be used to accelerate the 3PCF algorithm of Slepian & Eisenstein (2015). In both cases, these FT methods allow one to avoid the computational cost of pair counting, which scales as the square of the number density of objects in the survey. With the upcoming large datasets of DESI, Euclid, and LSST, FT techniques will therefore offer an important complement to simple pair or triplet counts.

We present an algorithm that computes the multipole coefficients of the galaxy three-point correlation function (3PCF) without explicitly considering triplets of galaxies. Rather, centering on each galaxy in the survey, it expands the radially-binned density field in spherical harmonics and combines these to form the multipoles without ever requiring the relative angle between a pair about the central. This approach scales with number and number density in the same way as the two-point correlation function, allowing runtimes that are comparable, and 500 times faster than a naive triplet count. It is exact in angle and easily handles edge correction. We demonstrate the algorithm on the LasDamas SDSS-DR7 mock catalogs, computing an edge corrected 3PCF out to $90\;{\rm Mpc}/h$ in under an hour on modest computing resources. We expect this algorithm will render it possible to obtain the large-scale 3PCF for upcoming surveys such as Euclid, LSST, and DESI.

We develop a configuration-space picture of the relative velocity between baryons and dark matter that clearly explains how it can shift the BAO scale in the galaxy-galaxy correlation function. The shift occurs because the relative velocity is non-zero only within the sound horizon and thus adds to the correlation function asymmetrically about the BAO peak. We further show that in configuration space the relative velocity has a localized, distinctive signature in the three-point galaxy correlation function (3PCF). In particular, we find that a multipole decomposition is a favorable way to isolate the relative velocity in the 3PCF, and that there is a strong signature in the l=1 multipole for triangles with 2 sides around the BAO scale. Finally, we investigate a further compression of the 3PCF to a function of only one triangle side that preserves the localized nature of the relative velocity signature while also nicely separating linear from non-linear bias. We expect that this scheme will substantially lessen the computational burden of finding the relative velocity in the 3PCF. The relative velocity's 3PCF signature can be used to correct the shift induced in the galaxy-galaxy correlation function so that no systematic error due to this effect is introduced into the BAO as used for precision cosmology.

Numerous upcoming observations, such as WFIRST, BOSS, BigBOSS, LSST, Euclid, and Planck, will constrain dark energy (DE)'s equation of state with great precision. They may well find the ratio of pressure to energy density, $w$, is -1, meaning DE is equivalent to a cosmological constant. However, many time-varying DE models have also been proposed. A single parametrization to test a broad class of them and that is itself motivated by a physical picture is therefore desirable. We suggest the simplest model of DE has the same mechanism as inflation, likely a scalar field slowly rolling down its potential. If this is so, DE will have a generic equation of state and the Universe will have a generic dependence of the Hubble constant on redshift independent of the potential's starting value and shape. This equation of state and expression for the Hubble constant offer the desired model-independent but physically motivated parametrization, because they will hold for most of the standard scalar-field models of DE such as quintessence and phantom DE. Up until now two-parameter descriptions of $w$ have been available, but this work finds an additional approximation that leads to a single-parameter model. Using it, we conduct a $\chi^2$ analysis and find that experiments in the next seven years should be able to distinguish any of these time-varying DE models on the one hand from a cosmological constant on the other to 73% confidence if $w$ today differs from -1 by 3.5%. In the limit of perfectly accurate measurements of $\Omega_m$ and $H_0$, this confidence would rise to 96%. We also include discussion of the current status of DE experiment, a table compiling the techniques each will use, and tables of the precisions of the experiments for which this information was available at the time of publication.

Self-interacting dark matter (SIDM), especially bosonic, has been considered a promising candidate to replace cold dark matter (CDM) as it resolves some of the problems associated with CDM. Here, we rule out the possibility that dark matter is a repulsive boson in thermal equilibrium. We develop the model first proposed by Goodman (2000) and derive the equation of state at finite temperature. Isothermal spherical halo models indicate a Bose-Einstein condensed core surrounded by a non-degenerate envelope, with an abrupt density drop marking the boundary between the two phases. Comparing this feature with observed rotation curves constrains the interaction strength of our model's DM particle, and Bullet Cluster measurements constrain the scattering cross section. Both ultimately can be cast as constraints on the particle's mass. We find these two constraints cannot be satisfied simultaneously in any realistic halo model---and hence dark matter cannot be a repulsive boson in thermal equilibrium. It is still left open that DM may be a repulsive boson provided it is not in thermal equilibrium; this requires that the mass of the particle be significantly less than a millivolt.

It is well known that in 3-d the average projected area of a convex solid is 1/4 the surface area. In this work, we generalize this theorem to higher dimensions by computing the analogous ratio as a function of the dimension. We prove a method for calculating this ratio in higher dimensions. We use this method to obtain both a recursion relation for this ratio from dimension d to d+1 and an explicit formula for it. We discuss the limiting behavior as the dimension becomes infinite and also mention possible application of the theorem as a test of the dimensionality of space.

We propose a simple model for dark energy useful for comparison with observations. It is based on the idea that dark energy and inflation should be caused by the same physical process. As motivation, we note that Linde's simple chaotic inflation $V=(1/2)m^{2}\phi^{2}$ produces values of $n_{s}=0.967$ and $r=0.13$, which are consistent with the WMAP 1-sigma error bars. We therefore propose $V=(1/2)m_{2}^{2}\phi_{2}^{2}+(1/2)m_{1}^{2}\phi_{1}^{2}$ with $m_{1}\sim10^{-5}$ and $m_{2}\leq10^{-60}$, where $c=1=\hbar$ and the reduced Planck mass is set to unity. The field $\phi_{1}$ drives inflation and has damped by now ($\phi_{1,0}=0$), while $\phi_{2}$ is currently rolling down its potential to produce dark energy. Using this model, we derive the formula $\delta w(z)\equiv w(z)+1=\delta w_{0}(H_0/H(z))^2$ via the slow-roll approximation. Our numerical results from exact and self-consistent solution of the equations of motion for $\phi_2$ and the Friedmann equations support this formula, and it should hold for any slow-roll dark energy. Our potential can be easily realized in N-flation models with many fields, and is easily falsifiable by upcoming experiments -- for example, if Linde's chaotic inflation is ruled out. But if $r$ values consistent with Linde's chaotic inflation are detected then one should take this model seriously indeed.

A few examples are given of Chandra's work on statistical and stochastic problems that relate to open questions in astrophysics, in particular his theory of dynamical relaxation in systems with inverse-square interparticle forces. The roles of chaos and integrability in this theory require clarification, especially for systems having a dominant central mass. After this prelude, a hypothetical form of repulsive bosonic dark matter is discussed. The repulsion leads to nontrivial thermodynamic behavior, including superfluidity, and would tend to suppress dynamical friction, greatly reducing the drag exerted on rotating galactic bars. However, this form of dark matter can probably be ruled out, at least for parameters that allow halos to reach thermal equilibria within a Hubble time. One combination of the particle mass and interparticle repulsion determines the minimum core radius of dark halos. Bounds on dark-matter collisionality inferred from the Bullet Cluster constrain a second combination. It is possible to satisfy both constraints only for parameters that predict unacceptable rotation curves outside the halo core.

We study the dynamics of a vertically thin, dispersion-dominated disk of planetesimals with eccentricities $e$ and inclinations $i$ (normalized in Hill units) satisfying $e >> 1$, $i << e^{-2} << 1$. This situation may be typical for e.g. a population of protoplanetary cores in the end of the oligarchic phase of planet formation. In this regime of orbital parameters planetesimal scattering has an anisotropic character and strongly differs from scattering in thick ($i ~ e$) disks. We derive analytical expressions for the planetesimal scattering coefficients and compare them with numerical calculations. We find significant discrepancies in the inclination scattering coefficients obtained by the two approaches and ascribe this difference to the effects not accounted for in the analytical calculation: multiple scattering events (temporary captures, which may be relevant for the production of distant planetary satellites outside the Hill sphere) and distant interaction of planetesimals prior to their close encounter. Our calculations show that the inclination of a thin, dispersion-dominated planetesimal disk grows exponentially on a very short time scale implying that (1) such disks must be very short-lived and (2) planetesimal accretion in this dynamical phase is insignificant. Our results are also applicable to the dynamics of shear-dominated disks switching to the dispersion-dominated regime.