
Particle physics has an ambitious and broad experimental programme for the
coming decades. This programme requires large investments in detector hardware,
either to build new facilities and experiments, or to upgrade existing ones.
Similarly, it requires commensurate investment in the R&D of software to
acquire, manage, process, and analyse the shear amounts of data to be recorded.
In planning for the HLLHC in particular, it is critical that all of the
collaborating stakeholders agree on the software goals and priorities, and that
the efforts complement each other. In this spirit, this white paper describes
the R&D activities required to prepare for this software upgrade.

We present MCFR, a multicasting concurrent face routing algorithm that uses
geometric routing to deliver a message from source to multiple targets. We
describe the algorithm's operation, prove it correct, estimate its performance
bounds and evaluate its performance using simulation. Our estimate shows that
MCFR is the first geometric multicast routing algorithm whose message delivery
latency is independent of network size and only proportional to the distance
between the source and the targets. Our simulation indicates that MCFR has
significantly better reliability than existing algorithms.

We show how observations of the evolution of the galaxy cluster number
abundance can be used to constrain primordial nonGaussianity in the universe.
We carry out a maximum likelihood analysis incorporating a number of current
datasets and accounting for a wide range of sources of systematic error. Under
the assumption of Gaussianity, the current data prefer a universe with matter
density $\Omega_m\simeq 0.3$ and are inconsistent with $\Omega_m=1$ at the
$2\sigma$ level. If we assume $\Omega_m=1$, the predicted degree of cluster
evolution is consistent with the data for nonGaussian models where the
primordial fluctuations have at least two times as many peaks of height
$3\sigma$ or more as a Gaussian distribution does. These results are robust to
almost all sources of systematic error considered: in particular, the
$\Omega_m=1$ Gaussian case can only be reconciled with the data if a number of
systematic effects conspire to modify the analysis in the right direction.
Given an independent measurement of $\Omega_m$, the techniques described here
represent a powerful tool with which to constrain nonGaussianity in the
primordial universe, independent of specific details of the nonGaussian
physics. We discuss the prospects and strategies for improving the constraints
with future observations.

I compare various semianalytic models for the bias of dark matter halos with
halo clustering properties observed in recent numerical simulations. The best
fitting model is one based on the collapse of ellipsoidal perturbations
proposed by Sheth, Mo & Tormen (1999), which fits the halo correlation length
to better than 8 per cent accuracy. Using this model, I confirm that the
correlation length of clusters of a given separation depends primarily on the
shape and amplitude of mass fluctuations in the universe, and is almost
independent of other cosmological parameters. Current observational
uncertainties make it difficult to draw robust conclusions, but for
illustrative purposes I discuss the constraints on the mass power spectrum
which are implied by recent analyses of the APM cluster sample. I also discuss
the prospects for improving these constraints using future surveys such as the
Sloan Digital Sky Survey. Finally, I show how these constraints can be combined
with observations of the cluster number abundance to place strong limits on the
matter density of the universe.

We carry out Nbody simulations of several nonGaussian structure formation
models, including Peebles' isocurvature cold dark matter model, cosmic string
models, and a model with primordial voids. We compare the evolution of the
cluster mass function in these simulations with that predicted by a modified
version of the PressSchechter formalism. We find that the PressSchechter
formula can accurately fit the cluster evolution over a wide range of redshifts
for all of the models considered, with typical errors in the mass function of
less than 25%, considerably smaller than the amount by which predictions for
different models may differ. This work demonstrates that the PressSchechter
formalism can be used to place strong model independent constraints on
nonGaussianity in the universe.

The WavelengthOriented Microwave Background Analysis Team (WOMBAT) is
constructing microwave maps which will be more realistic than previous
simulations. Our foreground models represent a considerable improvement: where
spatial templates are available for a given foreground, we predict the flux and
spectral index of that component at each place on the sky and estimate
uncertainties. We will produce maps containing simulated CMB anisotropy
combined with expected foregrounds. The simulated maps will be provided to the
community as the WOMBAT Challenge, so such maps can be analyzed to extract
cosmological parameters by scientists who are unaware of their input values.
This will test the efficacy of foreground subtraction, power spectrum analysis,
and parameter estimation techniques and help identify the areas most in need of
progress. These maps are also part of the FORECAST project, which allows
webbased access to the known foreground maps for the planning of CMB missions.

The WavelengthOriented Microwave Background Analysis Team (WOMBAT) is
constructing microwave skymaps which will be more realistic than previous
simulations. Our foreground models represent a considerable improvement: where
spatial templates are available for a given foreground, we predict the flux and
spectral index of that component at each place on the sky and estimate the
uncertainties in these quantities. We will produce maps containing simulated
Cosmic Microwave Background anisotropies combined with all major expected
foreground components. The simulated maps will be provided to the cosmology
community as the WOMBAT Challenge, a "hounds and hares" exercise where such
maps can be analyzed to extract cosmological parameters by scientists who are
unaware of their input values. This exercise will test the efficacy of current
foreground subtraction, power spectrum analysis, and parameter estimation
techniques and will help identify the areas most in need of progress.

We consider constraints on the amplitude of mass fluctuations in the
universe, sigma_8, derived from two simple observations: the present number
density of clusters and the amplitude of their correlation function. Allowing
for the possibility that the primordial fluctuations are nongaussian
introduces a degeneracy in the value of sigma_8 preferred by each of these
constraints. However, when the constraints are taken together this degeneracy
is broken, yielding a precise determination of sigma_8 and the degree of
nongaussianity for a given background cosmology. For a flat, Omega_m=1
universe with a power spectrum parameterized by a CDM shape parameter
Gamma=0.2, we find that the perturbations are consistent with a gaussian
distribution with sigma_8=0.49(+0.080.07) (95% limits). For some popular
choices of background model, including the favored low matter density models,
the hypothesis that the primordial fluctuations are gaussian is ruled out with
a high degree of confidence.

We calculate predictions from defect models of structure formation for both
the matter and Cosmic Microwave Background (CMB) over all observable scales.
Our results point to a serious problem reconciling the observed largescale
galaxy distribution with the COBE normalization, a result which is robust for a
wide range of defect parameters. We conclude that standard scaling defect
models are in conflict with the data, and show how attempts to resolve the
problem by considering nonscaling defects would require radical departures
from the standard scaling picture.

We examine the statistical properties of defects formed by the breaking of a
U(1) symmetry when the Higgs field has a power spectrum $P(k) \propto k^n$. We
find a marked dependence of the amount of infinite string on the spectral index
$n$ and empirically identify an analytic form for this quantity. We also
confirm that this result is robust to changes in the definition of infinite
string. It is possible that this result could account for the apparent absence
of infinite string in recent latticefree simulations.

We describe an implementation of the structure functions of Babul \&
Starkman, in order to quantify the ``sheetlike'' nature of a distribution of
matter. We test this statistic on a toy model describing cosmic string wakes,
and show that it does a better job than other statistics which have been
proposed for distinguishing nonGaussianity in the form of sheets. We conclude
that the most favoured cosmic string model is unlikely to produce a significant
increase in the sheetlike nature of the matter distribution beyond that which
occurs in Gaussian models (with the same power spectrum) due to the formation
of Zeldovich pancakes. Although the statistic was developed in the context of
cosmic string wake formation, we expect it to be useful for comparing the
observed galaxy distribution with a wide range of theoretical models with
different power spectra.

We find constraints on the generation of supercausalhorizon energy
perturbations from a smooth initial state, under a simple physical scheme. We
quantify these constraints by placing the upper limit $\lambda_c = 3.0 d_H$ on
the wavelength at which the power spectrum turns over to $k^4$ behavior. This
means that subhorizon processes can generate significant power on scales
further outside the horizon than one might naively expect. The existence of
this limit may have important implications for the interpretation of the small
scale power spectrum of the Cosmic Microwave Background.