
Quantumclassical hybrid algorithms are emerging as promising candidates for
nearterm practical applications of quantum information processors in a wide
variety of fields ranging from chemistry to physics and materials science. We
report on the experimental implementation of such an algorithm to solve a
quantum chemistry problem, using a digital quantum simulator based on trapped
ions. Specifically, we implement the variational quantum eigensolver algorithm
to calculate the molecular ground state energies of two simple molecules and
experimentally demonstrate and compare different encoding methods using up to
four qubits. Furthermore, we discuss the impact of measurement noise as well as
mitigation strategies and indicate the potential for adaptive implementations
focused on reaching chemical accuracy, which may serve as a crossplatform
benchmark for multiqubit quantum simulators.

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.

The variational quantum eigensolver (VQE) algorithm combines the ability of
quantum computers to efficiently compute expectation values with a classical
optimization routine in order to approximate ground state energies of quantum
systems. In this paper, we study the application of VQE to the simulation of
molecular energies using the unitary coupled cluster (UCC) ansatz. We introduce
new strategies to reduce the circuit depth for the implementation of UCC and
improve the optimization of the wavefunction based on efficient classical
approximations of the cluster amplitudes. Additionally, we propose an
analytical method to compute the energy gradient that reduces the sampling cost
for gradient estimation by several orders of magnitude compared to numerical
gradients. We illustrate our methodology with numerical simulations for a
system of four hydrogen atoms that exhibit strong correlation and show that the
circuit depth of VQE using a UCC ansatz can be reduced without introducing
significant loss of accuracy in the final wavefunctions and energies.

We show that qubit stabilizer states can be represented by nonnegative
quasiprobability distributions associated with a WignerWeylMoyal formalism
where Clifford gates are positive stateindependent maps. This is accomplished
by generalizing the WignerWeylMoyal formalism to three generators instead of
twoproducing an exterior, or Grassmann, algebrawhich results in Clifford
group gates for qubits that act as a permutation on the finite Weyl phase space
points naturally associated with stabilizer states. As a result, a nonnegative
probability distribution can be associated with each stabilizer state's
threegenerator Wigner function, and these distributions evolve
deterministically to one another under Clifford gates. This corresponds to a
hidden variable theory that is noncontextual and local for qubit Clifford
gates while Clifford (Pauli) measurements have a contextdependent
representation. Equivalently, we show that qubit Clifford gates can be
expressed as propagators within the threegenerator WignerWeylMoyal formalism
whose semiclassical expansion is truncated at order $\hbar^0$ with a finite
number of terms. The $T$gate, which extends the Clifford gate set to one
capable of universal quantum computation, require a semiclassical expansion of
the propagator to order $\hbar^1$. We compare this approach to previous
quasiprobability descriptions of qubits that relied on the twogenerator
WignerWeylMoyal formalism and find that the twogenerator Weyl symbols of
stabilizer states result in a description of evolution under Clifford gates
that is statedependent, in contrast to the threegenerator formalism. We have
thus extended Wigner nonnegative quasiprobability distributions from the odd
$d$dimensional case to $d=2$ qubits, which describe the noncontextuality of
Clifford gates and contextuality of Pauli measurements on qubit stabilizer
states.

The GottesmanKnill theorem established that stabilizer states and operations
can be efficiently simulated classically. For qudits with dimension three and
greater, stabilizer states and Clifford operations have been found to
correspond to positive discrete Wigner functions and dynamics. We present a
discrete Wigner functionbased simulation algorithm for odd$d$ qudits that has
the same time and space complexity as the AaronsonGottesman algorithm. We show
that the efficiency of both algorithms is due to the harmonic evolution in the
symplectic structure of discrete phase space. The differences between the
Wigner function algorithm and AaronsonGottesman are likely due only to the
fact that the WeylHeisenberg group is not in $SU(d)$ for $d=2$ and that qubits
have stateindependent contextuality. This may provide a guide for extending
the discrete Wigner function approach to qubits.

We give a path integral formulation of the time evolution of qudits of odd
dimension. This allows us to consider semiclassical evolution of discrete
systems in terms of an expansion of the propagator in powers of $\hbar$. The
largest power of $\hbar$ required to describe the evolution is a traditional
measure of classicality. We show that the action of the Clifford operators on
stabilizer states can be fully described by a single contribution of a path
integral truncated at order $\hbar^0$ and so are "classical," just like
propagation of Gaussians under harmonic Hamiltonians in the continuous case.
Such operations have no dependence on phase or quantum interference.
Conversely, we show that supplementing the Clifford group with gates necessary
for universal quantum computation results in a propagator consisting of a
finite number of semiclassical path integral contributions truncated at order
$\hbar^1$ , a number that nevertheless scales exponentially with the number of
qudits. The same sum in continuous systems has an infinite number of terms at
order $\hbar^1$.