
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.

Among the components contributing to particle transport, geometry navigation
is an important consumer of CPU cycles. The tasks performed to get answers to
"basic" queries such as locating a point within a geometry hierarchy or
computing accurately the distance to the next boundary can become very
computing intensive for complex detector setups. So far, the existing geometry
algorithms employ mainly scalar optimisation strategies (voxelization, caching)
to reduce their CPU consumption. In this paper, we would like to take a
different approach and investigate how geometry navigation can benefit from the
vector instruction set extensions that are one of the primary source of
performance enhancements on current and future hardware. While on paper, this
form of microparallelism promises increasing performance opportunities,
applying this technology to the highly hierarchical and multiply branched
geometry code is a difficult challenge. We refer to the current work done to
vectorise an important part of the critical navigation algorithms in the ROOT
geometry library. Starting from a short critical discussion about the
programming model, we present the current status and first benchmark results of
the vectorisation of some elementary geometry shape algorithms. On the path
towards a full vectorbased geometry navigator, we also investigate the
performance benefits in connecting these elementary functions together to
develop algorithms which are entirely based on the flow of vectordata. To this
end, we discuss core components of a simple vector navigator that is tested and
evaluated on a toy detector setup.

We investigate the groundstate properties of the highly degenerate
noncoplanar phase of the classical bilinearbiquadratic Heisenberg model on
the triangular lattice with Monte Carlo simulations. For that purpose, we
introduce an Ising pseudospin representation of the ground states, and we use a
simple Metropolis algorithm with local updates, as well as a powerful cluster
algorithm. At sizes that can be sampled with local updates, the presence of
longrange order is surprisingly combined with an algebraic decay of
correlations and the complete disordering of the chirality. It is only thanks
to the investigation of unusually large systems (containing $\sim 10^8$ spins)
with cluster updates that the true asymptotic regime can be reached and that
the system can be proven to consist of equivalent (i.e., equally ordered)
sublattices. These largescale simulations also demonstrate that the scalar
chirality exhibits longrange order at zero temperature, implying that the
system has to undergo a finitetemperature phase transition. Finally, we show
that the average distance in the order parameter space, which has the structure
of an infinite Cayley tree, remains remarkably small between any pair of
points, even in the limit when the real space distance between them tends to
infinity.

Using extensive classical and quantum Monte Carlo simulations, we investigate
the groundstate phase diagram of the fully frustrated transverse field Ising
model on the square lattice. We show that pure columnar order develops in the
lowfield phase above a surprisingly large length scale, below which an
effective U(1) symmetry is present. The same conclusion applies to the Quantum
Dimer Model with purely kinetic energy, to which the model reduces in the
zerofield limit, as well as to the stacked classical version of the model. By
contrast, the 2D classical version of the model is shown to develop plaquette
order. Semiclassical arguments show that the transition from plaquette to
columnar order is a consequence of quantum fluctuations.

We report on large scale finitetemperature Monte Carlo simulations of the
classical $120^\circ$ or $e_g$ orbitalonly model on the simple cubic lattice
in three dimensions with a focus towards its critical properties. This model
displays a continuous phase transition to an orbitally ordered phase. While the
correlation length exponent $\nu\approx0.665$ is close to the 3D XY value, the
exponent $\eta \approx 0.15$ differs substantially from O(N) values. We also
introduce a discrete variant of the $e_g$ model, called $e_g$clock model,
which is found to display the same set of exponents. Further, an emergent U(1)
symmetry is found at the critical point $T_c$, which persists for $T<T_c$ below
a crossover length scaling as $\Lambda \sim \xi^a$, with an unusually small
$a\approx1.3$.

We perform large scale finitetemperature Monte Carlo simulations of the
classical $e_g$ and $t_{2g}$ orbital models on the simple cubic lattice in
three dimensions. The $e_g$ model displays a continuous phase transition to an
orbitally ordered phase. While the correlation length exponent
$\nu\approx0.66(1)$ is close to the 3D XY value, the exponent $\eta \approx
0.15(1)$ differs substantially from O(N) values. At $T_c$ a U(1) symmetry
emerges, which persists for $T<T_c$ below a crossover length scaling as
$\Lambda \sim \xi^a$, with an unusually small $a\approx1.3$. Finally, for the
$t_{2g}$ model we find a {\em first order} transition into a lowtemperature
latticenematic phase without orbital order.

We study the directionalordering transition in the twodimensional classical
and quantum compass models on the square lattice by means of Monte Carlo
simulations. An improved algorithm is presented which builds on the Wolff
cluster algorithm in onedimensional subspaces of the configuration space. This
improvement allows us to study classical systems up to $L=512$. Based on the
new algorithm we give evidence for the presence of strongly anomalous scaling
for periodic boundary conditions which is much worse than anticipated before.
We propose and study alternative boundary conditions for the compass model
which do not make use of extended configuration spaces and show that they
completely remove the problem with finitesize scaling. In the last part, we
apply these boundary conditions to the quantum problem and present a
considerably improved estimate for the critical temperature which should be of
interest for future studies on the compass model. Our investigation identifies
a strong onedimensional magnetic ordering tendency with a large correlation
length as the cause of the unusual scaling and moreover allows for a precise
quantification of the anomalous length scale involved.

We present a pseudospin model which should be experimentally accessible using
solidstate devices and, being a variation on the compass model, adds to the
toolbox for the protection of qubits in the area of quantum information. Using
Monte Carlo methods, we find for both classical and quantum spins in two and
three dimensions Isingtype Neel ordering of energy fluctuations at finite
temperatures without magnetic order. We also readdress the controversy
concerning the stability of the ordered state in the presence of quenched
impurities and present numerical results which are at clear variance with
earlier claims in the literature.

We study two planar square lattice Heisenberg models with explicit
dimerization or quadrumerization of the couplings in the form of ladder and
plaquette arrangements. We investigate the quantum critical points of those
models by means of (stochastic series expansion) quantum Monte Carlo
simulations as a function of the coupling ratio $\alpha = J^\prime/J$. The
critical point of the orderdisorder quantum phase transition in the ladder
model is determined as $\alpha_\mathrm{c} = 1.9096(2)$ improving on previous
studies. For the plaquette model we obtain $\alpha_\mathrm{c} = 1.8230(2)$
establishing a first benchmark for this model from quantum Monte Carlo
simulations. Based on those values we give further convincing evidence that the
models are in the threedimensional (3D) classical Heisenberg universality
class. The results of this contribution shall be useful as references for
future investigations on planar Heisenberg models such as concerning the
influence of nonmagnetic impurities at the quantum critical point.

The twodimensional $J$$J^\prime$ dimerized quantum Heisenberg model is
studied on the square lattice by means of (stochastic series expansion) quantum
Monte Carlo simulations as a function of the coupling ratio
\hbox{$\alpha=J^\prime/J$}. The critical point of the orderdisorder quantum
phase transition in the $J$$J^\prime$ model is determined as
\hbox{$\alpha_\mathrm{c}=2.5196(2)$} by finitesize scaling for up to
approximately $10 000$ quantum spins. By comparing six dimerized models we
show, contrary to the current belief, that the critical exponents of the
$J$$J^\prime$ model are not in agreement with the threedimensional classical
Heisenberg universality class. This lends support to the notion of nontrivial
critical excitations at the quantum critical point.

A comprehensive study of the twodimensional (2D) compass model on the square
lattice is performed for classical and quantum spin degrees of freedom using
Monte Carlo and quantum Monte Carlo methods. We employ stateoftheart
implementations using Metropolis, stochastic series expansion and parallel
tempering techniques to obtain the critical ordering temperatures and critical
exponents. In a preinvestigation we reconsider the classical compass model
where we study and contrast the finitesize scaling behavior of ordinary
periodic boundary conditions against annealed boundary conditions. It is shown
that periodic boundary conditions suffer from extreme finitesize effects which
might be caused by closed loop excitations on the torus. These excitations also
appear to have severe effects on the Binder parameter. On this footing we
report on a systematic Monte Carlo study of the quantum compass model. Our
numerical results are at odds with recent literature on the subject which we
trace back to neglecting the strong finitesize effects on periodic lattices.
The critical temperatures are obtained as $T_\mathrm{c}=0.1464(2)J$ and
$T_\mathrm{c}=0.055(1)J$ for the classical and quantum version, respectively,
and our data support a transition in the 2D Ising universality class for both
cases.

The compact Abelian Higgs model is simulated on a cubic lattice where it
possesses vortex lines and pointlike magnetic monopoles as topological defects.
The focus of this highprecision Monte Carlo study is on the vortex network,
which is investigated by means of percolation observables. In the region of the
phase diagram where the Higgs and confinement phases are separated by a
firstorder transition, it is shown that the vortices percolate right at the
phase boundary, and that the firstorder nature of the transition is reflected
by the network. In the crossover region, where the phase boundary ceases to be
first order, the vortices are shown to still percolate. In contrast to other
observables, the percolation observables show finitesize scaling. The
exponents characterizing the critical behavior of the vortices in this region
are shown to fall in the random percolation universality class.

It is argued that the phase diagram of the 3D Compact U(1) Lattice Higgs
Model is more refined than generally thought. The confined and Higgs phases are
separated by a welldefined phase boundary, marked by proliferating vortices.
It is shown that the confinement mechanism at work is precisely the dual
superconductor scenario.

The threedimensional lattice Higgs model with compact U(1) gauge symmetry
and unit charge is investigated by means of Monte Carlo simulations. The full
model with fluctuating Higgs amplitude is simulated, and both energy as well as
topological observables are measured. The data show a Higgs and a confined
phase separated by a welldefined phase boundary, which is argued to be caused
by proliferating vortices. For fixed gauge coupling, the phase boundary
consists of a line of firstorder phase transitions at small Higgs
selfcoupling, ending at a critical point. The phase boundary then continues as
a Kertesz line across which thermodynamic quantities are nonsingular. Symmetry
arguments are given to support these findings.