
The DESI Legacy Imaging Surveys are a combination of three public projects
(the Dark Energy Camera Legacy Survey, the BeijingArizona Sky Survey, and the
Mayall zband Legacy Survey) that will jointly image ~14,000 square degrees 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 InterAmerican 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 dynamic observing that results in a survey
of nearly uniform depth. In addition to calibrated images, the project is
delivering an inferencebased catalog which includes photometry from the grz
optical bands and from four midinfrared bands (at 3.4um, 4.6um, 12um and 22um)
observed by the Widefield 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.

Low redshift measurements of Baryon Acoustic Oscillations (BAO) test the late
time evolution of the Universe and are a vital probe of Dark Energy. Over the
past decade both the 6degree Field Galaxy Survey (6dFGS) and Sloan Digital Sky
Survey (SDSS) have provided important distance constraints at $z < 0.3$. In
this paper we reevaluate the cosmological information from the BAO detection
in 6dFGS making use of HOD populated COLA mocks for a robust covariance matrix
and taking advantage of the now commonly implemented technique of density field
reconstruction. For the 6dFGS data, we find consistency with the previous
analysis, and obtain an isotropic volume averaged distance measurement of
$D_{V}(z_{\mathrm{eff}}=0.097) =
372\pm17(r_{s}/r_{s}^{\mathrm{fid}})\,\mathrm{Mpc}$, which has a nonGaussian
likelihood outside the $1\sigma$ region. We combine our measurement from both
the postreconstruction clustering of 6dFGS and SDSS MGS offering the most
robust constraint to date in this redshift regime,
$D_{V}(z_{\mathrm{eff}}=0.122)=539\pm17(r_{s}/r^{\mathrm{fid}}_{s})\,\mathrm{Mpc}$.
These measurements are consistent with standard $\Lambda\mathrm{CDM}$ and after
fixing the standard ruler using a Planck prior on $\Omega_{m}h^{2}$, the joint
analysis gives $H_{0}=64.0\pm3.5\,\mathrm{kms}^{1}\mathrm{Mpc}^{1}$. In the
near future both the Taipan Galaxy Survey and the Dark Energy Spectroscopic
Instrument (DESI) will improve this measurement to $1\%$ at low redshift.

The peculiar velocities of galaxies cause their redshiftspace clustering to
depend on the angle to the lineofsight, providing a key test of gravitational
physics on cosmological scales. These effects may be described using a
multipole expansion of the clustering measurements. Focussing on Fourierspace
statistics, we present a new analysis of the effect of the survey window
function, and the variation of the lineofsight across a survey, on the
modelling of power spectrum multipoles. We determine the joint covariance of
the Fourierspace multipoles in a Gaussian approximation, and indicate how
these techniques may be extended to studies of overlapping galaxy populations
via multipole crosspower spectra. We apply our methodology to one of the
widestarea galaxy redshift surveys currently available, the 6degree Field
Galaxy Survey, deducing a normalized growth rate f*sigma_8 = 0.38 +/ 0.12 in
the lowredshift Universe, in agreement with previous analyses of this dataset
using different techniques. Our framework should be useful for processing
future wideangle galaxy redshift surveys.

We study invasion fronts in the FitzHughNagumo equation in the oscillatory
regime using singular perturbation techniques. Phenomenologically, localized
perturbations of the unstable steadystate grow and spread, creating temporal
oscillations whose phase is modulated spatially. The phase modulation appears
to be selected by an invasion front that describes the behavior in the leading
edge of the spreading process. We construct these invasion fronts for large
regions of parameter space using singular perturbation techniques. Key
ingredients are the construction of periodic orbits, their unstable manifolds,
and the analysis of pushed and pulled fronts in the fast system. Our results
predict the wavenumbers and frequencies of oscillations in the wake of the
front through a phase locking mechanism. We also identify a parameter regime
where nonlinear phase locked fronts are inaccessible in the singularly
perturbed geometry of the travelingwave equation. Direct simulations confirm
our predictions and point to interesting phase slip dynamics.

We present a new test of gravitational physics by comparing the growth rate
of cosmic structure measured around voids with that measured around galaxies in
the same largescale structure dataset, the lowredshift 6degree Field Galaxy
Survey. By fitting a selfconsistent Redshift Space Distortion model to the 2D
galaxygalaxy and voidgalaxy correlation functions, we recover growth rate
values f\sigma_8 = 0.42 \pm 0.06 and 0.39 \pm 0.11, respectively. The
environmentaldependence of cosmological statistics can potentially
discriminate between modifiedgravity scenarios which modulate the growth rate
as a function of scale or environment and test the underlying assumptions of
homogeneity and isotropy.

We investigate functions with the property that for every interval, the slope
at the midpoint of the interval is the same as the average slope. More
generally, we find functions whose average slopes over intervals are given by
the slope at a weighted average of the endpoints of those intervals. This is
equivalent to finding functions satisfying a weighted mean value property. In
the course of our exploration, we find connections to harmonic functions that
prompt us to explore multivariable analogues and the existence of "weighted
harmonic functions."