
2
various quantum algorithms have been proposed, such
as quantum Monte Carlo, Hamiltonian simulation, and
imaginary-time evolution algorithms. The Hamiltonian
simulation, similar to TN approaches, allows quantum
algorithms to avoid the sign problem mentioned above.
However, it is essential to note that these techniques have
their own limitations. Due to the limited number of
“noise free” qubits or qubits to control the error, the em-
bedding of the Hamiltonian is significantly restricted. For
scenarios with infinite Hilbert space dimensions, Hamil-
tonian truncation is especially crucial [35,36].
In particular, nonlinear sigma models have been in-
vestigated in spin-lattice systems [37] and with quantum
computing techniques [38,39]. In the latter case, it is
demonstrated that nonlinear sigma models can have a
qubit-efficient description through a fuzzy-sphere repre-
sentation. Although this allows for a Hamiltonian trun-
cation that is suitable for quantum time evolution algo-
rithms, it is not prone to a straightforward generalisation
to nonlinear sigma models that feature a richer structure.
In this work, we go beyond this limitation and augment a
two-dimensional nonlinear sigma model with a topologi-
cal θ-term. In particular, using quantum-gate simulators,
we aim to study the entanglement entropy of the vacuum
to investigate the critical behaviour and determine the
mass gap of the quantum theory.
In addition to the simulation of QFTs, both TN-
inspired quantum circuits and conventional TN methods
have also been used in various machine-learning applica-
tions. It has been shown that their relation with quantum
many-body systems can be used to achieve more inter-
pretable networks [40–43].
This work is organised as follows. In Section II, we re-
view the Hamiltonian formulation of the O(3) nonlinear
sigma model at θ=πin terms of angular momentum
variables on a one-dimensional spin chain. Section III
shows how these can be embedded on a quantum com-
puter. In particular, we use this approach to compute
the bipartite entanglement entropy associated with the
half chain and we confirm the vanishing mass gap of the
theory in Section IV. Finally, we briefly summarise our
results and conclude in Section V.
II. HAMILTONIAN FORMULATION OF
NONLINEAR SIGMA MODELS
The field content of a general O(3) nonlinear sigma
model is given by a real vector field nthat takes values
on a sphere, S2, i.e. it is normalised to n2= 1. In a two-
dimensional Euclidean spacetime, the associated action
is commonly written as
S=1
2g2Zd2x(∂µn)2.(1)
Here, gis the dimensionless coupling constant, and we
consider Euclidean coordinates τand x. At the classical
level, the vector field is massless. This remains true to
all orders in perturbation theory. Nevertheless, it can be
shown that the quantum theory is gapped [7–11].
As the two-dimensional O(3) nonlinear sigma model
admits instanton solutions (see, e.g., [2,3]), it is feasible
to consider an additional topological term in this theory.
Along these lines, we can distinguish finite-action field
configurations by their topological charge, Q=Rd2x ρQ,
where
ρQ=1
4πabcna∂xnb∂τnc(2)
is the topological charge density. These field configura-
tions, in turn, contribute a finite θ-term to the action,
Sθ=1
2g2Zd2x(∂µn)2+iθQ . (3)
Since the topological charge is an integer, the action is
2π-periodic with respect to θ. In the following, we aim to
investigate the topological term’s effect on the quantum
theory’s mass gap. For concreteness, let us focus on the
case θ=πin the following. As we will see momentarily,
this choice allows for a simple Hamiltonian formulation
of the QFT. We will comment on the general case later
in this work.
For a numerical investigation of the two-dimensional
O(3) nonlinear sigma model using quantum algorithms,
a suitable Hamiltonian formulation of the former is
needed [7,44,45]. Somewhat fortunately, this allows us
to treat the two-dimensional theory from an effectively
one-dimensional perspective as follows. First, one con-
siders the theory on a discrete spatial axis while keeping
the time coordinate continuous at the same time. In this
scenario, the time derivative of the kinetic term can be
identified with an angular momentum per lattice site.
Therefore, a one-dimensional chain of coupled quantum
rotors can describe the two-dimensional field theory. In
particular, for the nonlinear sigma model at θ=π, the
Hamiltonian can be written as [12]
H=1
2β
N
X
k=1
L2
k+β
N−1
X
k=1
nknk+1 .(4)
Here, we set the lattice spacing in the spatial direction
to a= 1 and make the replacement β= 1/g2. In
this setup, Lkdenotes a (modified) quantum mechani-
cal angular momentum operator acting on the k-th site
of the spin chain.3While this operator acts on the local
Hilbert space at each site, the second term of the Hamil-
tonian corresponds to the interactions of neighbouring
sites, which we will characterise momentarily. Note that,
at this stage, we impose open boundary conditions to
keep the notation simple. This is why the summation
3We again remark that this spin chain fully characterises the two-
dimensional nonlinear sigma model. This is due to the choice of
Hamiltonian variables (see, e.g., [45]).