Fast radio bursts (FRBs) are bright, millisecond-duration radio transients originating from extragalactic distances. Their origin is unknown. Some FRB sources emit repeat bursts, ruling out cataclysmic origins for those events. Despite searches for periodicity in repeat burst arrival times on time scales from milliseconds to many days, these bursts have hitherto been observed to appear sporadically, and though clustered, without a regular pattern. Here we report the detection of a $16.35\pm0.15$ day periodicity (or possibly a higher-frequency alias of that periodicity) from a repeating FRB 180916.J0158+65 detected by the Canadian Hydrogen Intensity Mapping Experiment Fast Radio Burst Project (CHIME/FRB). In 38 bursts recorded from September 16th, 2018 through February 4th, 2020, we find that all bursts arrive in a 5-day phase window, and 50% of the bursts arrive in a 0.6-day phase window. Our results suggest a mechanism for periodic modulation either of the burst emission itself, or through external amplification or absorption, and disfavour models invoking purely sporadic processes.

Magnetars are highly magnetized young neutron stars that occasionally produce enormous bursts and flares of X-rays and gamma-rays. Of the approximately thirty magnetars currently known in our Galaxy and Magellanic Clouds, five have exhibited transient radio pulsations. Fast radio bursts (FRBs) are millisecond-duration bursts of radio waves arriving from cosmological distances. Some have been seen to repeat. A leading model for repeating FRBs is that they are extragalactic magnetars, powered by their intense magnetic fields. However, a challenge to this model has been that FRBs must have radio luminosities many orders of magnitude larger than those seen from known Galactic magnetars. Here we report the detection of an extremely intense radio burst from the Galactic magnetar SGR 1935+2154 using the Canadian Hydrogen Intensity Mapping Experiment (CHIME) FRB project. The fluence of this two-component bright radio burst and the estimated distance to SGR 1935+2154 together imply a 400-800 MHz burst energy of $\sim 3 \times 10^{34}$ erg, which is three orders of magnitude brighter than those of any radio-emitting magnetar detected thus far. Such a burst coming from a nearby galaxy would be indistinguishable from a typical FRB. This event thus bridges a large fraction of the radio energy gap between the population of Galactic magnetars and FRBs, strongly supporting the notion that magnetars are the origin of at least some FRBs.

The Packed Ultra-wideband Mapping Array (PUMA) is a proposed low-resolution transit interferometric radio telescope operating over the frequency range 200 - 1100MHz. Its rich science portfolio will include measuring structure in the universe from redshift z = 0.3 to 6 using 21cm intensity mapping, detecting one million fast radio bursts, and monitoring thousands of pulsars. It will allow PUMA to advance science in three different areas of physics (the physics of dark energy, the physics of cosmic inflation and time-domain astrophysics). This document is a response to a request for information (RFI) by the Panel on Radio, Millimeter, and Submillimeter Observations from the Ground (RMS) of the Decadal Survey on Astronomy and Astrophysics 2020. We present the science case of PUMA, the development path and major risks to the project.

The redshifted 21 cm line is an emerging tool in cosmology, in principle permitting three-dimensional surveys of our Universe that reach unprecedentedly large volumes, previously inaccessible length scales, and hitherto unexplored epochs of our cosmic timeline. Large radio telescopes have been constructed for this purpose, and in recent years there has been considerable progress in transforming 21 cm cosmology from a field of considerable theoretical promise to one of observational reality. Increasingly, practitioners in the field are coming to the realization that the success of observational 21cm cosmology will hinge on software algorithms and analysis pipelines just as much as it does on careful hardware design and telescope construction. This review provides a pedagogical introduction to state-of-the-art ideas in 21 cm data analysis, covering a wide variety of steps in a typical analysis pipeline, from calibration to foreground subtraction to mapmaking to power spectrum estimation to parameter estimation.

This white paper envisions a revolutionary post-DESI, post-LSST dark energy program based on intensity mapping of the redshifted 21cm emission line from neutral hydrogen at radio frequencies. The proposed intensity mapping survey has the unique capability to quadruple the volume of the Universe surveyed by optical programs, provide a percent-level measurement of the expansion history to $z \sim 6$, open a window to explore physics beyond the concordance $\Lambda$CDM model, and to significantly improve the precision on standard cosmological parameters. In addition, characterization of dark energy and new physics will be powerfully enhanced by cross-correlations with optical surveys and cosmic microwave background measurements. The rich dataset obtained by the proposed intensity mapping instrument will be simultaneously useful in exploring the time-domain physics of fast radio transients and pulsars, potentially in live "multi-messenger" coincidence with other observatories. The core dark energy/inflation science advances enabled by this program are the following: (i) Measure the expansion history of the universe over $z=0.3-6$ with a single instrument, extending the range deep into the pre-acceleration era, providing an unexplored window for new physics; (ii) Measure the growth rate of structure in the universe over the same redshift range; (iii) Observe, or constrain, the presence of inflationary relics in the primordial power spectrum, improving existing constraints by an order of magnitude; (iv) Observe, or constrain, primordial non-Gaussianity with unprecedented precision, improving constraints on several key numbers by an order of magnitude. Detailed mapping of the enormous, and still largely unexplored, volume of cosmic space will thus provide unprecedented information on fundamental questions of the vacuum energy and early-universe physics.

Development of the hardware, data analysis, and simulation techniques for large compact radio arrays dedicated to mapping the 21 cm line of neutral hydrogen gas has proven to be more difficult than imagined twenty years ago when such telescopes were first proposed. Despite tremendous technical and methodological advances, there are several outstanding questions on how to optimally calibrate and analyze such data. On the positive side, it has become clear that the outstanding issues are purely technical in nature and can be solved with sufficient development activity. Such activity will enable science across redshifts, from early galaxy evolution in the pre-reionization era to dark energy evolution at low redshift.

PUMA is a proposal for an ultra-wideband, low-resolution and transit interferometric radio telescope operating at $200-1100\,\mathrm{MHz}$. Its design is driven by six science goals which span three science themes: the physics of dark energy (measuring the expansion history and growth of the universe up to $z=6$), the physics of inflation (constraining primordial non-Gaussianity and primordial features) and the transient radio sky (detecting one million fast radio bursts and following up SKA-discovered pulsars). We propose two array configurations composed of hexagonally close-packed 6m dish arrangements with 50% fill factor. The initial 5,000 element 'petite array' is scientifically compelling, and can act as a demonstrator and a stepping stone to the full 32,000 element 'full array'. Viewed as a 21cm intensity mapping telescope, the program has the noise equivalent of a traditional spectroscopic galaxy survey comprised of 0.6 and 2.5 billion galaxies at a comoving wavenumber of $k=0.5\,h\mathrm{Mpc}^{-1}$ spanning the redshift range $z = 0.3 - 6$ for the petite and full configurations, respectively. At redshifts beyond $z=2$, the 21cm technique is a uniquely powerful way of mapping the universe, while the low-redshift range will allow for numerous cross-correlations with existing and upcoming surveys. This program is enabled by the development of ultra-wideband radio feeds, cost-effective dish construction methods, commodity radio-frequency electronics driven by the telecommunication industry and the emergence of sufficient computing power to facilitate real-time signal processing that exploits the full potential of massive radio arrays. The project has an estimated construction cost of 55 and 330 million FY19 USD for the petite and full array configurations. Including R&D, design, operations and science analysis, the cost rises to 125 and 600 million FY19 USD, respectively.

Radio interferometers consisting of identical antennas arranged on a regular lattice permit fast Fourier transform beamforming, which reduces the correlation cost from $\mathcal{O}(n^2)$ in the number of antennas to $\mathcal{O}(n\log n)$. We develop a formalism for describing this process and apply this formalism to derive a number of algorithms with a range of observational applications. These include algorithms for forming arbitrarily pointed tied-array beams from the regularly spaced Fourier-transform formed beams, sculpting the beams to suppress sidelobes while only losing percent-level sensitivity, and optimally estimating the position of a detected source from its observed brightness in the set of beams. We also discuss the effect that correlations in the visibility-space noise, due to cross-talk and sky contributions, have on the optimality of Fourier transform beamforming, showing that it does not strictly preserve the sky information of the $n^2$ correlation, even for an idealized array. Our results have applications to a number of upcoming interferometers, in particular the Canadian Hydrogen Intensity Mapping Experiment--Fast Radio Burst (CHIME/FRB) project.

Fast radio bursts (FRBs) are highly dispersed millisecond-duration radio flashes likely arriving from far outside the Milky Way galaxy. This phenomenon was discovered at radio frequencies near 1.4 GHz and to date has been observed in one case at as high as 8 GHz, but not below 700 MHz in spite of significant searches at low frequencies. Here we report detections of FRBs at radio frequencies as low as 400 MHz, on the Canadian Hydrogen Intensity Mapping Experiment (CHIME) using the CHIME/FRB instrument. We present 13 FRBs detected during a telescope pre-commissioning phase, when our sensitivity and field-of-view were not yet at design specifications. Emission in multiple events is seen down to 400 MHz, the lowest radio frequency to which we are sensitive. The FRBs show a variety of temporal scattering behaviours, with the majority significantly scattered, and some apparently unscattered to within measurement uncertainty even at our lowest frequencies. Of the 13 reported here, one event has the lowest dispersion measure yet reported, implying it is among the closest yet known, and another has shown multiple repeat bursts, as described in a companion paper. Our low-scattering events suggest that efforts to detect FRBs at radio frequencies below 400 MHz will eventually be successful. The overall scattering properties of our sample suggest that FRBs as a class are preferentially located in environments that scatter radio waves more strongly than the diffuse interstellar medium (ISM) in the Milky Way.

The discovery of a repeating Fast Radio Burst (FRB) source, FRB 121102, eliminated models involving cataclysmic events for this source. No other repeating FRB has yet been detected in spite of many recent FRB discoveries and follow-ups, suggesting repeaters may be rare in the FRB population. Here we report the detection of six repeat bursts from FRB 180814.J0422+73, one of the 13 FRBs detected by the Canadian Hydrogen Intensity Mapping Experiment (CHIME) FRB project during its pre-commissioning phase in July and August 2018. These repeat bursts are consistent with originating from a single position on the sky, with the same dispersion measure (DM), ~189 pc cm-3. This DM is approximately twice the expected Milky Way column density, and implies an upper limit on the source redshift of 0.1, at least a factor of ~2 closer than FRB 121102. In some of the repeat bursts, we observe sub-pulse frequency structure, drifting, and spectral variation reminiscent of that seen in FRB 121102, suggesting similar emission mechanisms and/or propagation effects. This second repeater, found among the first few CHIME/FRB discoveries, suggests that there exists -- and that CHIME/FRB and other wide-field, sensitive radio telescopes will find -- a substantial population of repeating FRBs.

In radio interferometry, the quantization process introduces a bias in the magnitude and phase of the measured correlations which translates into errors in the measurement of source brightness and position in the sky, affecting both the system calibration and image reconstruction. In this paper we investigate the biasing effect of quantization in the measured correlation between complex-valued inputs with a circularly symmetric Gaussian probability density function (PDF), which is the typical case for radio astronomy applications. We start by calculating the correlation between the input and quantization error and its effect on the quantized variance, first in the case of a real-valued quantizer with a zero mean Gaussian input and then in the case of a complex-valued quantizer with a circularly symmetric Gaussian input. We demonstrate that this input-error correlation is always negative for a quantizer with an odd number of levels, while for an even number of levels this correlation is positive in the low signal level regime. In both cases there is an optimal interval for the input signal level for which this input-error correlation is very weak and the model of additive uncorrelated quantization noise provides a very accurate approximation. We determine the conditions under which the magnitude and phase of the measured correlation have negligible bias with respect to the unquantized values: we demonstrate that the magnitude bias is negligible only if both unquantized inputs are optimally quantized (i.e., when the uncorrelated quantization error model is valid), while the phase bias is negligible when 1) at least one of the inputs is optimally quantized, or when 2) the correlation coefficient between the unquantized inputs is small. Finally, we determine the implications of these results for radio interferometry.

A new generation of wide-field radio interferometers designed for 21-cm surveys is being built as drift scan instruments allowing them to observe large fractions of the sky. With large numbers of antennas and frequency channels the enormous instantaneous data rates of these telescopes require novel, efficient, data management and analysis techniques. The $m$-mode formalism exploits the periodicity of such data with the sidereal day, combined with the assumption of statistical isotropy of the sky, to achieve large computational savings and render optimal analysis methods computationally tractable. We present an extension to that work that allows us to adopt a more realistic sky model and treat objects such as bright point sources. We develop a linear procedure for deconvolving maps, using a Wiener filter reconstruction technique, which simultaneously allows filtering of these unwanted components. We construct an algorithm, based on the Sherman-Morrison-Woodbury formula, to efficiently invert the data covariance matrix, as required for any optimal signal-to-noise weighting. The performance of our algorithm is demonstrated using simulations of a cylindrical transit telescope.

We present results from a new incoherent-beam Fast Radio Burst (FRB) search on the Canadian Hydrogen Intensity Mapping Experiment (CHIME) Pathfinder. Its large instantaneous field of view (FoV) and relative thermal insensitivity allow us to probe the ultra-bright tail of the FRB distribution, and to test a recent claim that this distribution's slope, $\alpha\equiv-\frac{\partial \log N}{\partial \log S}$, is quite small. A 256-input incoherent beamformer was deployed on the CHIME Pathfinder for this purpose. If the FRB distribution were described by a single power-law with $\alpha=0.7$, we would expect an FRB detection every few days, making this the fastest survey on sky at present. We collected 1268 hours of data, amounting to one of the largest exposures of any FRB survey, with over 2.4\,$\times$\,10$^5$\,deg$^2$\,hrs. Having seen no bursts, we have constrained the rate of extremely bright events to $<\!13$\,sky$^{-1}$\,day$^{-1}$ above $\sim$\,220$\sqrt{(\tau/\rm ms)}$ Jy\,ms for $\tau$ between 1.3 and 100\,ms, at 400--800\,MHz. The non-detection also allows us to rule out $\alpha\lesssim0.9$ with 95$\%$ confidence, after marginalizing over uncertainties in the GBT rate at 700--900\,MHz, though we show that for a cosmological population and a large dynamic range in flux density, $\alpha$ is brightness-dependent. Since FRBs now extend to large enough distances that non-Euclidean effects are significant, there is still expected to be a dearth of faint events and relative excess of bright events. Nevertheless we have constrained the allowed number of ultra-intense FRBs. While this does not have significant implications for deeper, large-FoV surveys like full CHIME and APERTIF, it does have important consequences for other wide-field, small dish experiments.

We have developed FFT beamforming techniques for the CHIME radio telescope, to search for and localize the astrophysical signals from Fast Radio Bursts (FRBs) over a large instantaneous field-of-view (FOV) while maintaining the full angular resolution of CHIME. We implement a hybrid beamforming pipeline in a GPU correlator, synthesizing 256 FFT-formed beams in the North-South direction by four formed beams along East-West via exact phasing, tiling a sky area of ~250 square degrees. A zero-padding approximation is employed to improve chromatic beam alignment across the wide bandwidth of 400 to 800 MHz. We up-channelize the data in order to achieve fine spectral resolution of $\Delta\nu$=24 kHz and time cadence of 0.983 ms, desirable for detecting transient and dispersed signals such as those from FRBs.

The Canadian Hydrogen Intensity Mapping Experiment (CHIME) Pathfinder radio telescope is currently surveying the northern hemisphere between 400 and 800 MHz. By mapping the large scale structure of neutral hydrogen through its redshifted 21 cm line emission between $z \sim 0.8-2.5$ CHIME will contribute to our understanding of Dark Energy. Bright astrophysical foregrounds must be separated from the neutral hydrogen signal, a task which requires precise characterization of the polarized telescope beams. Using the DRAO John A. Galt 26 m telescope, we have developed a holography instrument and technique for mapping the CHIME Pathfinder beams. We report the status of the instrument and initial results of this effort.

Close tidal encounters among large planetesimals and moons should have been more common than grazing or normal impacts. Using a mass spring model within an N-body simulation, we simulate the deformation of the surface of an elastic spherical body caused by a close parabolic tidal encounter with a body that has similar mass as that of the primary body. Such an encounter can induce sufficient stress on the surface to cause brittle failure of an icy crust and simulated fractures can extend a large fraction of the radius of body. Strong tidal encounters may be responsible for the formation of long graben complexes and chasmata in ancient terrain of icy moons such as Dione, Tethys, Ariel and Charon.

We present a procedure for efficiently compressing astronomical radio data for high performance applications. Integrated, post-correlation data are first passed through a nearly lossless rounding step which compares the precision of the data to a generalized and calibration-independent form of the radiometer equation. This allows the precision of the data to be reduced in a way that has an insignificant impact on the data. The newly developed Bitshuffle lossless compression algorithm is subsequently applied. When the algorithm is used in conjunction with the HDF5 library and data format, data produced by the CHIME Pathfinder telescope is compressed to 28% of its original size and decompression throughputs in excess of 1 GB/s are obtained on a single core.

We present the design and implementation of a custom GPU-based compute cluster that provides the correlation X-engine of the CHIME Pathfinder radio telescope. It is among the largest such systems in operation, correlating 32,896 baselines (256 inputs) over 400MHz of radio bandwidth. Making heavy use of consumer-grade parts and a custom software stack, the system was developed at a small fraction of the cost of comparable installations. Unlike existing GPU backends, this system is built around OpenCL kernels running on consumer-level AMD GPUs, taking advantage of low-cost hardware and leveraging packed integer operations to double algorithmic efficiency. The system achieves the required 105TOPS in a 10kW power envelope, making it among the most power-efficient X-engines in use today.

The Canadian Hydrogen Intensity Mapping Experiment (CHIME) is an ambitious new radio telescope project for measuring cosmic expansion and investigating dark energy. Keeping good records of both physical configuration of its 1280 antennas and their analogue signal chains as well as the ~100 TB of data produced daily from its correlator will be essential to the success of CHIME. In these proceedings we describe the database-driven software we have developed to manage this complexity.

A pathfinder version of CHIME (the Canadian Hydrogen Intensity Mapping Experiment) is currently being commissioned at the Dominion Radio Astrophysical Observatory (DRAO) in Penticton, BC. The instrument is a hybrid cylindrical interferometer designed to measure the large scale neutral hydrogen power spectrum across the redshift range 0.8 to 2.5. The power spectrum will be used to measure the baryon acoustic oscillation (BAO) scale across this poorly probed redshift range where dark energy becomes a significant contributor to the evolution of the Universe. The instrument revives the cylinder design in radio astronomy with a wide field survey as a primary goal. Modern low-noise amplifiers and digital processing remove the necessity for the analog beamforming that characterized previous designs. The Pathfinder consists of two cylinders 37\,m long by 20\,m wide oriented north-south for a total collecting area of 1,500 square meters. The cylinders are stationary with no moving parts, and form a transit instrument with an instantaneous field of view of $\sim$100\,degrees by 1-2\,degrees. Each CHIME Pathfinder cylinder has a feedline with 64 dual polarization feeds placed every $\sim$30\,cm which Nyquist sample the north-south sky over much of the frequency band. The signals from each dual-polarization feed are independently amplified, filtered to 400-800\,MHz, and directly sampled at 800\,MSps using 8 bits. The correlator is an FX design, where the Fourier transform channelization is performed in FPGAs, which are interfaced to a set of GPUs that compute the correlation matrix. The CHIME Pathfinder is a 1/10th scale prototype version of CHIME and is designed to detect the BAO feature and constrain the distance-redshift relation.

The Canadian Hydrogen Intensity Mapping Experiment (CHIME) is a transit interferometer currently being built at the Dominion Radio Astrophysical Observatory (DRAO) in Penticton, BC, Canada. We will use CHIME to map neutral hydrogen in the frequency range 400 -- 800\,MHz over half of the sky, producing a measurement of baryon acoustic oscillations (BAO) at redshifts between 0.8 -- 2.5 to probe dark energy. We have deployed a pathfinder version of CHIME that will yield constraints on the BAO power spectrum and provide a test-bed for our calibration scheme. I will discuss the CHIME calibration requirements and describe instrumentation we are developing to meet these requirements.

In this paper we continue to develop the m-mode formalism, a technique for efficient and optimal analysis of wide-field transit radio telescopes, targeted at 21 cm cosmology. We extend this formalism to give an accurate treatment of the polarised sky, fully accounting for the effects of polarisation leakage and cross-polarisation. We use the geometry of the measured set of visibilities to project down to pure temperature modes on the sky, serving as a significant compression, and an effective first filter of polarised contaminants. We use the m-mode formalism with the Karhunen-Loeve transform to give a highly efficient method for foreground cleaning, and demonstrate its success in cleaning realistic polarised skies observed with an instrument suffering from substantial off axis polarisation leakage. We develop an optimal quadratic estimator in the m-mode formalism, which can be efficiently calculated using a Monte-Carlo technique. This is used to assess the implications of foreground removal for power spectrum constraints where we find that our method can clean foregrounds well below the foreground wedge, rendering only scales $k_\parallel < 0.02 h \,\mathrm{Mpc}^{-1}$ inaccessible. As this approach assumes perfect knowledge of the telescope, we perform a conservative test of how essential this is by simulating and analysing datasets with deviations about our assumed telescope. Assuming no other techniques to mitigate bias are applied, we recover unbiased power spectra when the per-feed beam width to be measured to 0.1%, and amplifier gains to be known to 1% within each minute. Finally, as an example application, we extend our forecasts to a wideband 400-800 MHz cosmological observation and consider the implications for probing dark energy, finding a medium-sized cylinder telescope improves the DETF Figure of Merit by around 70% over Planck and Stage II experiments alone.

We model a 21 cm intensity mapping survey in the redshift range 0.01<z<1.5 designed to simulate the skies as seen by future radio telescopes such as the Square Kilometre Array (SKA), including instrumental noise and Galactic foregrounds. In our pipeline, we remove the introduced Galactic foregrounds with a fast independent component analysis (fastica) technique. We present the power spectrum of the large-scale matter distribution, C(l), before and after the application of this foreground removal method and calculate the resulting systematic errors. We attempt to reduce systematics in the foreground subtraction by optimally masking the maps to remove high foregrounds in the Galactic plane. Our simulations show a certain level of bias remains in the power spectrum at all scales l<400. At large-scales l<30 this bias is particularly significant. We measure the impact of these systematic effects in two different ways: firstly we fit cosmological parameters to the broadband shape of the power spectrum and secondly we extract the position of the Baryon Acoustic Oscillations (BAO). In the first analysis, we find that the systematics introduce an significant shift in the best fit cosmological parameters at the 2 to 3 sigma level which depends on the masking and noise levels. However, cosmic distances can be recovered in an unbiased way after foreground removal at all simulated redshifts by fitting the BAOs in the power spectrum. We conclude that further advances in foreground removal are needed in order to recover unbiased information from the broadband shape of the power spectrum, however, intensity mapping experiments will be a powerful tool for mapping cosmic distances across a wide redshift range.

The large-scale distribution of neutral hydrogen in the Universe will be luminous through its 21 cm emission. Here, for the first time, we use the auto-power spectrum of 21 cm intensity fluctuations to constrain neutral hydrogen fluctuations at z~0.8. Our data were acquired with the Green Bank Telescope and span the redshift range 0.6 < z < 1 over two fields totalling ~41 deg. sq. and 190 h of radio integration time. The dominant synchrotron foregrounds exceed the signal by ~10^3, but have fewer degrees of freedom and can be removed efficiently. Even in the presence of residual foregrounds, the auto-power can still be interpreted as an upper bound on the 21 cm signal. Our previous measurements of the cross-correlation of 21 cm intensity and the WiggleZ galaxy survey provide a lower bound. Through a Bayesian treatment of signal and foregrounds, we can combine both fields in auto- and cross-power into a measurement of Omega_HI b_HI = [0.62^{+0.23}_{-0.15}] * 10^{-3} at 68% confidence with 9% systematic calibration uncertainty, where Omega_HI is the neutral hydrogen (HI) fraction and b_HI is the HI bias parameter. We describe observational challenges with the present data set and plans to overcome them.

In this paper we describe the spherical harmonic transit telescope, a novel formalism for the analysis of transit radio telescopes. This all-sky approach bypasses the curved sky complications of traditional interferometry and so is particularly well suited to the analysis of wide-field radio interferometers. It enables compact and computationally efficient representations of the data and its statistics that allow new ways of approaching important problems like map-making and foreground removal. In particular, we show how it enables the use of the Karhunen-Loeve transform as a highly effective foreground filter, suppressing realistic foreground residuals for our fiducial example by at least a factor twenty below the 21cm signal even in highly contaminated regions of the sky. This is despite the presence of the angle-frequency mode mixing inherent in real-world instruments with frequency-dependent beams. We show, using Fisher forecasting, that foreground cleaning has little effect on power spectrum constraints compared to hypothetical foreground-free measurements. Beyond providing a natural real-world data analysis framework for 21cm telescopes now under construction and future experiments, this formalism allows accurate power spectrum forecasts to be made that include the interplay of design constraints and realistic experimental systematics with twenty-first century 21cm science.

In this letter, 21 cm intensity maps acquired at the Green Bank Telescope are cross-correlated with large-scale structure traced by galaxies in the WiggleZ Dark Energy Survey. The data span the redshift range 0.6 < z < 1 over two fields totaling ~41 deg. sq. and 190 hours of radio integration time. The cross-correlation constrains Omega_HI b_HI r = [0.43 \pm 0.07 (stat.) \pm 0.04(sys.)] x 10^-3, where Omega_HI is the neutral hydrogen HI fraction, r is the galaxy-hydrogen correlation coefficient, and b_HI is the HI bias parameter. This is the most precise constraint on neutral hydrogen density fluctuations in a challenging redshift range. Our measurement improves the previous 21 cm cross-correlation at z ~ 0.8 both in its precision and in the range of scales probed.

Primordial magnetic fields could provide an explanation for the galactic magnetic fields observed today, in which case they may also leave interesting signals in the CMB and the small-scale matter power spectrum. We discuss how to approximately calculate the important non-linear magnetic effects within the guise of linear perturbation theory, and calculate the matter and CMB power spectra including the SZ contribution. We then use various cosmological datasets to constrain the form of the magnetic field power spectrum. Using solely large-scale CMB data (WMAP7, QUaD and ACBAR) we find a 95% CL on the variance of the magnetic field at 1 Mpc of B_\lambda < 6.4 nG. When we include SPT data to constrain the SZ effect, we find a revised limit of B_\lambda < 4.1 nG. The addition of SDSS Lyman-alpha data lowers this limit even further, roughly constraining the magnetic field to B_\lambda < 1.3 nG.

We report the discovery of an unusually red brown dwarf found in a search for high proper motion objects using WISE and 2MASS data. WISEP J004701.06+680352.1 is moving at 0.44$ arcsec/yr and lies relatively close to the Galactic Plane (b=5.2 degrees). Near-infrared photometry and spectroscopy reveals that this is one of the reddest (2MASS J-K_s = 2.55 +/- 0.08 mag) field L dwarfs yet detected, making this object an important member of the class of unusually red L dwarfs. We discuss evidence for thick condensate clouds and speculate on the age of the object. Although models by different research groups agree that thick clouds can explain the red spectrum, they predict dramatically different effective temperatures, ranging from 1100K to 1600K. This brown dwarf is well suited for additional studies of extremely dusty substellar atmospheres because it is relatively bright (K_s = 13.05 +/- 0.03 mag), which should also contribute to an improved understanding of young gas-giant planets and the transition between L and T brown dwarfs.

We use the new cosmological recombination code, CosmoRec, for parameter estimation in the context of (future) precise measurements of the CMB temperature and polarization anisotropies. We address the question of how previously neglected physical processes in the recombination model of Recfast affect the determination of key cosmological parameters, for the first time performing a model-by-model computation of the recombination problem. In particular we ask how the biases depend on different combinations of parameters, e.g. when varying the helium abundance or the effective number of neutrino species in addition to the standard six parameters. We also forecast how important the recombination corrections are for a combined Planck, ACTPol and SPTpol data analysis. Furthermore, we ask which recombination corrections are really crucial for CMB parameter estimation, and whether an approach based on a redshift-dependent correction function to Recfast is sufficient in this context.

Primordial magnetic fields and massive neutrinos can leave an interesting signal in the CMB temperature and polarization. We perform a systematic analysis of general perturbations in the radiation-dominated universe, accounting for any primordial magnetic field and including leading- order effects of the neutrino mass. We show that massive neutrinos qualitatively change the large- scale perturbations sourced by magnetic fields, but that the effect is much smaller than previously claimed. We calculate the CMB power spectra sourced by inhomogeneous primordial magnetic fields, from before and after neutrino decoupling, including scalar, vector and tensor modes, and consistently modelling the correlation between the density and anisotropic stress sources. In an appendix we present general series solutions for the possible regular primordial perturbations.

We present new photometric and spectroscopic observations for 2M 1533+3759 (= NSVS 07826147). It has an orbital period of 0.16177042 day, significantly longer than the 2.3--3.0 hour periods of the other known eclipsing sdB+dM systems. Spectroscopic analysis of the hot primary yields Teff = 29230 +/- 125 K, log g = 5.58 +/- 0.03 and log N(He)/N(H) = -2.37 +/- 0.05. The sdB velocity amplitude is K1 = 71.1 +/- 1.0 km/s. The only detectable light contribution from the secondary is due to the surprisingly strong reflection effect. Light curve modeling produced several solutions corresponding to different values of the system mass ratio, q(M2/M1), but only one is consistent with a core helium burning star, q=0.301. The orbital inclination is 86.6 degree. The sdB primary mass is M1 = 0.376 +/- 0.055 Msun and its radius is R1 = 0.166 +/- 0.007 Rsun. 2M1533+3759 joins PG0911+456 (and possibly also HS2333+3927) in having an unusually low mass for an sdB star. SdB stars with masses significantly lower than the canonical value of 0.48 Msun, down to as low as 0.30 Msun, were theoretically predicted by Han et al. (2002, 2003), but observational evidence has only recently begun to confirm the existence of such stars. The existence of core helium burning stars with masses lower than 0.40--0.43 Msun implies that at least some sdB progenitors have initial main sequence masses of 1.8--2.0 Msun or more, i.e. they are at least main sequence A stars. The secondary is a main sequence M5 star.

Distances in cosmology are usually inferred from observed redshifts - an estimate that is dependent on the local peculiar motion - giving a distorted view of the three dimensional structure and affecting basic observables such as the correlation function and power spectrum. We calculate the full non-linear redshift-space power spectrum for Gaussian fields, giving results for both the standard flat sky approximation and the directly-observable angular correlation function and angular power spectrum. Coupling between large and small scale modes boosts the power on small scales when the perturbations are small. On larger scales power is slightly suppressed by the velocities perturbations on smaller scales. The analysis is general, but we comment specifically on the implications for future high-redshift observations, and show that the non-linear spectrum has significantly more complicated angular structure than in linear theory. We comment on the implications for using the angular structure to separate cosmological and astrophysical components of 21 cm observations.

Bayesian model selection provides the cosmologist with an exacting tool to distinguish between competing models based purely on the data, via the Bayesian evidence. Previous methods to calculate this quantity either lacked general applicability or were computationally demanding. However, nested sampling (Skilling 2004), which was recently applied successfully to cosmology by Muhkerjee et al. 2006, overcomes both of these impediments. Their implementation restricts the parameter space sampled, and thus improves the efficiency, using a decreasing ellipsoidal bound in the n-dimensional parameter space centred on the maximum likelihood point. However, if the likelihood function contains any multi-modality, then the ellipse is prevented from constraining the sampling region efficiently. In this paper we introduce a method of clustered ellipsoidal nested sampling which can form multiple ellipses around each individual peak in the likelihood. In addition we have implemented a method for determining the expectation and variance of the final evidence value without the need to use sampling error from repetitions of the algorithm ; this further reduces the computational load by at least an order of magnitude. We have applied our algorithm to a pair of toy models and one cosmological example where we demonstrate that the number of likelihood evaluations required is ~ 4% of that necessary for using previous algorithms. We have produced a fortran library containing our routines which can be called from any sampling code, in addition for convenience we have incorporated it into the popular CosmoMC code as CosmoClust. Both are available for download at http://www.mrao.cam.ac.uk/software/cosmoclust .

We report on the discovery of J0644+3344, a bright deeply eclipsing cataclysmic variable (CV) binary. Spectral signatures of both binary components and an accretion disk can be seen at optical wavelengths. The optical spectrum shows broad H I, He I, and He II accretion disk emission lines with deep narrow absorption components from H I, He I, Mg II and Ca II. The absorption lines are seen throughout the orbital period, disappearing only during primary eclipse. These absorption lines are either the the result of an optically-thick inner accretion disk or from the photosphere of the primary star. Radial velocity measurements show that the H I, He I, and Mg II absorption lines phase with the the primary star, while weak absorption features in the continuum phase with the secondary star. Radial velocity solutions give a 150+/-4 km/s semi-amplitude for the primary star and 192.8+/-5.6 km/s for the secondary. The individual stellar masses are 0.63-0.69 Mdot for the primary and 0.49-0.54 Mdot for the secondary. The bright eclipsing nature of this binary has helped provide masses for both components with an accuracy rarely achieved for CVs. This binary most closely resembles a nova-like UX UMa or SW Sex type of CV. J0644+3344, however, has a longer orbital period than most UX UMa or SW Sex stars. Assuming an evolution toward shorter orbital periods, J0644+3344 is therefore likely to be a young interacting binary. The secondary star is consistent with the size and spectral type of a K8 star, but has an M0 mass.

We report the discovery of a new detached eclipsing binary known as NSVS01031772 with component masses M_1= 0.5428 +/- 0.0027Msun, M_2= 0.4982 +/- 0.0025Msun, and radii R_1= 0.5260 +/- 0.0028Rsun, R_2= 0.5088 +/- 0.0030Rsun. The system has an orbital period of 0.3681414(3) days and an apparent magnitude V ~ 12.6. The estimated effective temperatures of the stars are T_eff1 = 3615 +/- 72K and T_eff2 = 3513 +/- 31K. The space velocities of the system suggest that it is a main sequence binary and with a metallicity that is approximately solar. The two stars in this binary are located in a region of the Mass-Radius relation where no accurate observational data were previously available. Similarly to other low-mass binaries recently studied, the radius of each star in NSVS01031772 exceeds the best evolutionary model predictions by about 8.5% on average.

We present new observations of cool <6000K and low mass <1Msun binary systems that have been discovered by searching several modern stellar photometric databases. The search has led to a factor of 10 increase in the number of known cool close eclipsing binary systems.

Full tests to stellar models below 1 Msun have been hindered until now by the scarce number of precise measurements of the stars' most fundamental parameters: their masses and radii. With the current observational techniques, the required precision to distinguish between different models (errors < 2-3 %) can only be achieved using detached eclipsing binaries where 1) both stars are similar in mass, i.e. q = M1/M2 ~ 1.0, and 2) each star is a main sequence object below 1 Msun. Until 2003 only three such binaries had been found and analyzed in detail. Two new systems were published in 2005 (Creevey et al.; Lopez-Morales & Ribas), almost doubling the previous number of data points. Here we present preliminary results for 3 new low-mass detached eclipsing binaries. These are the first studied systems from our sample of over 41 new binaries (Shaw & Lopez-Morales, this proceedings). We also provide an updated comparison between the Mass-Radius and the Mass-Teff relations predicted by the models and the observational data from detached eclipsing binaries.

We report 323 hours of nearly uninterrupted time series photometric observations of the DBV star GD 358 acquired with the Whole Earth Telescope (WET) during May 23rd to June 8th, 2000. We acquired more than 232 000 independent measurements. We also report on 48 hours of time-series photometric observations in Aug 1996. We detected the non-radial g-modes consistent with degree l=1 and radial order 8 to 20 and their linear combinations up to 6th order.We also detect, for the first time, a high amplitude l=2 mode, with a period of 796s. In the 2000 WET data, the largest amplitude modes are similar to those detected with the WET observations of 1990 and 1994, but the highest combination order previously detected was 4th order. At one point during the 1996 observations, most of the pulsation energy was transferred into the radial order k=8 mode, which displayed a sinusoidal pulse shape in spite of the large amplitude. The multiplet structure of the individual modes changes from year to year, and during the 2000 observations only the k=9 mode displays clear normal triplet structure. Even though the pulsation amplitudes change on timescales of days and years, the eigenfrequencies remain essentially the same, showing the stellar structure is not changing on any dynamical timescale.