In this paper we consider quantum walks whose evolution converges to the Dirac equation in the limit of small wave-vectors. We show exact Fast Fourier implementation of the Dirac quantum walks in one, two, and three space dimensions. The behaviour of particle states—defined as states smoothly peaked in some wave-vector eigenstate of the walk—is described by an approximated dispersive differential equation that for small wave-vectors gives the usual Dirac particle and antiparticle kinematics. The accuracy of the approximation is provided in terms of a lower bound on the fidelity between the exactly evolved state and the approximated one. The jittering of the position operator expectation value for states having both a particle and an antiparticle component is analytically derived and observed in the numerical implementations.
In this paper we consider quantum walks whose evolution converges to the Dirac equation in the limit of small wave-vectors. We show exact Fast Fourier implementation of the Dirac quantum walks in one, two, and three space dimensions. The behaviour of particle states—defined as states smoothly peaked in some wave-vector eigenstate of the walk—is described by an approximated dispersive differential equation that for small wave-vectors gives the usual Dirac particle and antiparticle kinematics. The accuracy of the approximation is provided in terms of a lower bound on the fidelity between the exactly evolved state and the approximated one. The jittering of the position operator expectation value for states having both a particle and an antiparticle component is analytically derived and observed in the numerical implementations.
The Dirac equation can be modelled as a quantum walk (QW), whose main features are being: discrete in time and space (i.e. a unitary evolution of the wave-function of a particle on a lattice); homogeneous (i.e. translation-invariant and time-independent) and causal (i.e. information propagates at a bounded speed, in a strict sense). This link, which was proposed already by Succi and Benzi, Bialynicki-Birula and Meyer, is shown to hold for Bargmann–Wigner equations and symmetric hyperbolic systems in general. We then analytically prove the convergence of the solution of the QW to the solution of the Cauchy problem for the Dirac equation. We do so by adapting a powerful method from standard numerical analysis, which is of general interest to the field of quantum simulation. At the practical level, this result provides precise error bounds and convergence rates, thereby validating the QW as a quantum simulation scheme. At the theoretical level, it reinforces the status of this QW as a simple, discrete toy model of relativistic particles.
A quantum walk describes the discrete unitary evolution of a quantum particle on a discrete graph. Some quantum walks, referred to as the Weyl and Dirac walks, provide a description of the free evolution of relativistic quantum fields in the small wave-vector regime. The clash between the intrinsic discreteness of quantum walks and the continuous symmetries of special relativity is resolved by giving a definition of change of inertial frame in terms of a change of values of the constants of motion, which leaves the walk operator unchanged. Starting from the family of 1 + 1 dimensional Dirac walks with all possible values of the mass parameter, we introduce a unique walk encompassing the latter as an extra degree of freedom, and we derive its group of changes of inertial frames. This symmetry group contains a non-linear realization of SO+(2,1)⋉R3; since one of the two space-like dimensions does not correspond to an actual spatial degree of freedom but rather the mass, we interpret it as a 2 + 1 dimensional de Sitter group. This group also contains a non-linear realization of the proper orthochronous Poincaré group SO+(1,1)⋉R2 in 1 + 1 dimension, as the ones considered within the framework of doubly special relativity, which recovers the usual relativistic symmetry in the limit of small wave-vectors and masses. Surprisingly, for the Dirac walk with a fixed value of the mass parameter, the group of allowed changes of reference frame does not have a consistent interpretation in the limit of small wave-vectors.
Discrete-time quantum walks (QWs) are transportation models of single quantum particles over a lattice. Their evolution is driven through causal and local unitary operators. QWs are a powerful tool for quantum simulation of fundamental physics, as some of them have a continuum limit converging to well-known physics partial differential equations, such as the Dirac or the Schr\"odinger equation. In this paper, we show how to recover the Dirac equation in (3+1) dimensions with a QW evolving in a tetrahedral space. This paves the way to simulate the Dirac equation on a curved space-time. This also suggests an ordered scheme for propagating matter over a spin network, of interest in loop quantum gravity, where matter propagation has remained an open problem.
From the unitary operator used for implementing two-state discrete-time quantum walk on one-, two- and three- dimensional lattice we obtain a two-component Dirac-like Hamiltonian. In particular, using different pairs of Pauli basis as position translation states we obtain three different form of Hamiltonians for evolution on one-dimensional lattice. We extend this to two- and three-dimensional lattices using different Pauli basis states as position translation states for each dimension and show that the external coin operation, which is necessary for one-dimensional walk is not a necessary requirement for a walk on higher dimensions but can serve as an additional resource to control the dynamics. The two-component Hamiltonian we present here for quantum walk on different lattices can serve as a general framework to simulate, control, and study the dynamics of quantum systems governed by Dirac-like Hamiltonian.
We introduce a quantum algorithm for simulating the time-dependent Dirac equation in 3+1 dimensions using discrete-time quantum walks. Thus far, promising quantum algorithms have been proposed to simulate quantum dynamics in non-relativistic regimes efficiently. However, only some studies have attempted to simulate relativistic dynamics due to its theoretical and computational difficulty. By leveraging the convergence of discrete-time quantum walks to the Dirac equation, we develop a quantum spectral method that approximates smooth solutions with exponential convergence. This mitigates errors in implementing potential functions and reduces the overall gate complexity that depends on errors. We demonstrate that our approach does not require additional operations compared to the asymptotic gate complexity of non-relativistic real-space algorithms. Our findings indicate that simulating relativistic dynamics is achievable with quantum computers and can provide insights into relativistic quantum physics and chemistry.
The weak convergence theorems of the one- and two-dimensional simple quantum walks, SQW(d),d=1,2, show a striking contrast to the classical counterparts, the simple random walks, SRW(d). In the SRW(d), the distribution of position X(t) of the particle starting from the origin converges to the Gaussian distribution in the diffusion scaling limit, in which the time scale T and spatial scale L both go to infinity as the ratio L/T is kept finite. On the other hand, in the SQW(d), the ratio L/T is kept to define the pseudovelocity V(t)=X(t)/t, and then all joint moments of the components ${V}_{j}(t),1\ensuremath{\leqslant}j\ensuremath{\leqslant}d$, of V(t) converge in the $T=L\ensuremath{\rightarrow}\ensuremath{\infty}$ limit. The limit distributions have novel structures such that they are inverted-bell shaped and their supports are bounded. In the present paper we claim that these properties of the SQW(d) can be explained by the theory of relativistic quantum mechanics. We show that the Dirac equation with a proper ultraviolet cutoff can provide a quantum walk model in three dimensions, where the walker has a four-component qubit. We clarify that the pseudovelocity V(t) of the quantum walker, which solves the Dirac equation, is identified with the relativistic velocity. Since the quantum walker should be a tardyon, not a tachyon, ∣V(t)∣<c, where c is the speed of light, and this restriction (the causality) is the origin of the finiteness of supports of the limit distributions universally found in quantum walk models. By reducing the number of components of momentum in the Dirac equation, we obtain the limit distributions of pseudovelocities for the lower dimensional quantum walks. We show that the obtained limit distributions for the one- and two-dimensional systems have common features with those of SQW(1) and SQW(2).