Arxiv:2012.04237V4 [Quant-Ph] 1 Aug 2021

Quantum-to-classical transition via quantum cellular automata

Pedro C.S. Costa1, ∗ 1Department of Physics and Astronomy, Macquarie University, Sydney, NSW 2109, Australia (Dated: August 3, 2021) A quantum cellular automaton (QCA) is an abstract model consisting of an array of ﬁnite-dimensional quantum systems that evolves in discrete time by local unitary operations. Here we propose a simple coarse-graining map, where the spatial structure of the QCA is merged into effective ones. Starting with a QCA that simulates the Dirac equation, we apply this coarse-graining map recursively until we get its effective dynamics in the semiclassical limit, which can be described by a classical cellular automaton. We show that the emergent-effective result of the former microscopic discrete model converges to the diffusion equation and to a classical transport equation under a speciﬁc initial condition. Therefore, QCA is a good model to validate the quantum-to-classical transition.

I. INTRODUCTION the CA, likewise done in [8, 14]. That is, in order to an- alyze the emergent dynamics in the semi-classical limit, In a realistic scenario, quantum systems are never we assume that we cannot resolve neighboring cells and completely isolated from the environment and its quan- then some information has to be discarded. tum effects naturally start to become suppressed. Such a Although in this work is the indistinguishability be- phenomenon is known as the decoherence process, and tween the states the source of decoherence to study the it is the main ingredient used to explore the quantum- quantum-to-classical transition dynamics, the procedure to-classical transition [23]. Completely positive trace- here, in analogy to [5], can also be helpful to track the preserving (CPTP) maps are powerful tools to study such nonclassicality that emerges due to other sources of de- transition. These maps, which can be seen as a unitary in- coherence in the new quantum technologies. teraction among three systems and further discard two of The paper is organized as follows: In sectionII we de- them [29], are therefore not reversible and their applica- fined the quantum model of computation used to simulate tion implies a loss of information about the system. Sim- the Dirac equation in the continuous limit. In section III ilarly with the classical picture, where coarse-graining we introduce the coarse-graining map used to study the (CG) procedure via the projective operators [13, 17] decoherence process for the massless Dirac equation. In helps us to extract effective dynamics that have lower sectionIV the CG map is applied and our main results computational cost, the CPTP maps also provide a data are established. Finally, in sectionV, we conclude with compressor for the quantum system. Thus, these maps a discussion of these results and we point out possible when applied in quantum systems can also be seen as a directions to go for next. coarse-graining procedure. Rather than verify the necessary conditions to extract the quantum channels at each coarse-grained level [11, II. A PUQCA THAT SIMULATES THE (1+1) DIRAC EQUATION 12], here we concentrate on just one specific quantum dynamics, one map and one CG level. The Dirac equation is the quantum dynamics chosen, and the coarse- In one spatial and one temporal dimension, the Dirac graining level is the one where the coherence disappears. equation of a free spin-1/2 particle of mass m may be Given that the quantum model of computation, known written as follows: as quantum cellular automata (proposed in [7]) is able to h¯ arXiv:2012.04237v4 [quant-ph] 1 Aug 2021 simulate a large class of quantum walks, in particular the i ∂t + ih¯σz∂x − mcσx ψ = 0, (1) c one that simulates the Dirac equation, we use it to study the semiclassical limit of the Dirac equation [9]. Our where c is the velocity of light, h¯ is Planck’s constant procedure respects the same constraints between differ- divided by 2π, and σ j, with j ∈ {x,z}, are the usual ent CG levels for the time evolution shown in [12, 16], Pauli matrices. We want to construct a partitioned uni- which is represented in figure (3). However, our formal- tary quantum cellular automata (PUQCA) whose dynam- ism takes advantage of the space structure provided by ics, in the continuous limit, is the same as the Dirac equation (1). A PUQCA is a discrete model of quantum computation, which is most easily understood in its graph per- ∗ [email protected] spective. We start by discretizing the space coordinate, The restriction to the single excitation subspace guarantees that the number of particles is fixed to one. To inter- act the cells, we employ a simple SWAP gate between the subcells connected by an edge. That means, we take:

W1 = |00ih00| + |01ih10| + |10ih01| + |11ih11|. (4)

Given the above construction, we can now see the evolution on the one-particle sector of this QCA. To do that we write the system state at time t, |Ψ(t)i as: " FIG. 1. a Graph representation of the space coordinates. La- 0 ∑ ψ (k∆x,t)|...(00)(k−1)∆x(10)k∆x(00)(k+1)∆x ···i beling the vertex positions by k their positions in space can k∈Z be calculated, i.e., x = k∆x. ) The vertices are related to the # PUQCA cells, with two subcells each. 1 + ψ (k∆x,t)|···(00)(k−1)∆x(01)k∆x(00)(k+1)∆x ···i . (5) and thus constructing a 1D linear graph–see figure1 a. Each vertex in this graph represents a point in space, and Applying the evolution ε to this state the following re- we will assume that the distance between two vertices currence relation for the amplitudes can be derived: is ∆x. In this way, the vertices are located at positions x = k ∆x, with k ∈ Z. ψ0 (x,t + ∆t) = cosθψ0 (x − ∆x,t) − isinθψ1 (x − ∆x,t) In order to describe the dynamics of the particle, an ex- ψ1 (x,t + ∆t) = cosθψ1 (x + ∆x,t) − isinθψ0(x + ∆x,t). citation moving in the space, each vertex is split into two (6) subcells and in each subcell, we place a qubit–see fig- ure1 b. More concretely, to the vertex at position x we Associating ψ0 and ψ1 with the spinor components of ψ, associate a Hilbert space (H )x = (H0)x ⊗ (H1)x with the last step is just to take the continuous limit. Making 2 2 H0 ' H1 ' C , for the left and right subcells, respec- θ = mc ∆t/h¯, and taking the limit ∆x → 0 and ∆t → 0, tively. The state |0i will be associated with no particles, such that ∆x/∆t → c, in first order we get the massive while the state |1i will be associated with the presence of Dirac equation in (1+1)-dimensional spacetime (1). We a particle/excitation. call such quantum automaton as the Dirac PUQCA. In this model, the evolution is also discrete. Defin- ing the time steps size as ∆t, the values for time are t = m∆t, with m ∈ N. For the cell structure defined, III. COARSE-GRAINING A 1D PUQCA the transition function E , i.e., the operator that updates the PUQCA state | (t)i = | (t + t)i, is composed E Ψ Ψ ∆ Now that we have established the quantum cellular by two parts. First, each vertex is individually and si- automaton that simulates the 1+1 Dirac equation, we multaneously updated by the application of a unitary use the automaton spatial structure to define a coarse- (W ) ( ) ⊗ ( ) → ( ) ⊗ ( ) . Then, to al- 0 x : H0 x H1 x H0 x H1 x graining map. The physical idea is that an observer with low for the particle propagation a second unitary is ap- low resolution will not be able to discriminate between plied, but now acting on the graph edges, i.e., (W )x+∆x 1 x : different automaton cells. ( ) ⊗ ( ) → ( ) ⊗ ( ) . The transi- H1 x H0 x+∆x H1 x H0 x+∆x To incorporate this idea in our approach, we employ tion function is then completely positive and trace preserving (CPTP) maps, O (k+1)∆x O 0 a.k.a. quantum channels, to throw away inaccessible in- E = (W1)k∆x (W0)k ∆x. (2) k∈Z k0∈Z formation in a general way. This approach has been de- veloped and used in other scenarios [1, 10, 12, 15, 16, Following the construction of a quantum walk that 18, 19, 21, 24, 25, 27, 28]. In general, a coarse-graining simulates the 1+1 Dirac equation [20], and the translation channel will be defined as a quantum channel which is between quantum walks and the PUQCA [7], we choose a trace preserving and completely positive map that re- W the unitary transformations as follows. For 0 we allow, duces the description of the system: Λ : L (HD) → a rotation in the single-excitation subspace for each cell: L (Hd ) with D > d, where L (H ) is the set of all linear operators acting on . W = |00ih00| − isinθ|01ih01| + cosθ|01ih10| H 0 Here, we will apply a coarse-graining channel in each + cosθ|10ih01| − isinθ|10ih10| + |11ih11|. (3) cell of the automaton, reducing their two subcells into a

2 single subcell. Locally the coarse-graining map acts as one-particle sector, we can write the automaton state at follows: level L of description as

2 2 2 1 C ⊗ C → C k . k∆x,a 0 0 Λk∆x : L 0 k∆x 1 k∆x L k mod 2 b c∆x | k a ih k a | 2 ρL = ∑ ∑ (ρL)k0∆x,a0 ( , ) ( , ) , (10) (7) k,k0∈Z a,a0=0 The application of the coarse-graining map in the whole automaton is then given by: where |(k,a)i represents the state with an excitation at position k∆x and subcell a. The state at level L + 1, i.e., after we apply the coarse-graining map ΛCG in the whole O automaton, is given by: ΛCG = Λk∆x. (8) k∈Z 1 Notice that the application of the coarse-graining map (1 − δk,k0 ) ρL+1 = ∑ ∑ δk,k0 δa,a0 + preserves the automaton local structure. It is also impor- k,k0∈Z a,a=0 3 tant to realize that after “compressing” the two subcells k k0 × k∆x,a k k0 at a position x into an effective subcell, this subcell is (ρL)k0∆x,a0 , mod 2 , mod 2 . merged with a neighboring effective subcell to form an 2 2 effective cell. The distance between two effective cells is (11) still ∆x. By this procedure, we establish a new level of The derivation of this expression is done in Ap- description of the automaton, with a new 1-d lattice. pendixVIA. To finalize our procedure to obtain an effective de- From the above expression, we can derive, by recur- scription, we must define the action of Λx. To model sion, the density matrix coefficients for level L in terms the lack of resolution, we will employ a coarse-graining of the density matrix coefficients of level 0: map that cannot distinguish if an excitation is in the left L−1− or right subcell of a point x. The action of Λx, which was 2 1 1 L L−1 ( )x,b = ( )2 x+2 b∆x+s∆x,a , (12) first suggested in [12] in the context of atoms in a optical ρL x,b ∑ ∑ ρ0 2Lx+2L−1b∆x+s∆x,a lattice, is as follows: s=0 a=0 and √ Λx (|01ih01|) = |1ih1| Λx (|01ih00|) = |1ih0|/√3 2L−1−1 1 Λx (|10ih10|) = |1ih1| Λx (|10ih00|) = |1ih0|/ 3 x,b 1 2Lx+2L−1b∆x+s∆x,a (ρ ) 0 0 = (ρ ) , (|00ih00|) = |0ih0| (|10ih01|) = 0 L x ,b L ∑ ∑ 0 2Lx0+2L−1b0∆x+s0∆x,a0 Λx Λx 3 0 0 (9) s,s =0 a,a =0 (13) The first column contains the diagonal elements and for x 6= x0 and b, b0 ∈ {0,1}. From the recursive rela- expresses the inability to distinguish whether the exci- tion above for the off-diagonal elements of ρL, we show tation is on the left or the right subcell. If no excita- √ x,b ≤ L tion is present, then effectively no excitation is observed. in AppendixVIB that (ρL)x0,b0 (2 2/3) , i.e. the It is clearly observed that this map is not equivalent to coherence elements converge exponentially fast to zero a partial trace. The second column contains the off- with the coarse-graining level. To illustrate this exponen- diagonal terms. As the the base vectors in the excited tial decay, consider the case where at L = 0 we start√ the subspace cannot be discriminated, then there can be no automaton in the state |Ψ(0)i = (|(0,0)i+|(0,1)i)/ 2 coherence among them. The only possible coherence and we evolve it with the Dirac PUQCA, setting θ = π/4 terms are in between the no excitation subspace√ and the in equation (3), for 200 time steps. Figure2 shows how single-excitation subspace. The factor 1/ 3 guarantees the sum of all absolute values of the off-diagonal terms the complete positivity of Λx. As quantum channels pre- decays with level of coarse-graining. serve the Hermiticity, the action on the Hermitian conju- gated terms can be apprehended from (9). Further discussion on this map and its action on two-particle subspace IV. COARSE-GRAINED DIRAC PUQCA can be found in [12, 25]. Given the map Λx, we are now ready to find the In the previous section, we computed the coarse- coarse-grained 1D PUQCA. As the coarse-graining grained density operator for an arbitrary level L and we channel maps the 1-d lattice to another similar 1-d lat- wrote the recursive equations over different L. In this tice, we can apply this map L times on the underlying section, we will compute the coarse-grained QCA. Be- automation to obtain its Lth level of description. In the fore the coarse-graining map, the QCA is governed by

3 (1) (0) L−1 (0) (1) where αl ,αl ∈ R and ∑l=0 (αl + αl ) = 1. Now, for a given time t 2L from the massless dynamics, equation (6), since any superposition is created anymore when θ = 0, all amplitudes that started on the left(right) of the center remain in the subcells 0(1) that is

L−1 2 −1 L (t) x,0 = (t) 2 x+s∆x,0 (ρL )x,0 ∑ (ρ0 )2Lx+s∆x,0 (15) s=0 L 2 −1 L (t) x,1 = (t) 2 x+s∆x,1 (ρL )x,1 ∑ (ρ0 )2Lx+s∆x,1 . (16) FIG. 2. This plot shows the semilog sum of the off-diagonal s=2L−1 elements in different levels of CG for a fixed time. Starting with a QCA state at t = 0 in L = 0, the sum of the absolute Now, let us coarse-graining the density operator ρ at time values of the coherence elements is taken from L = 0 to L = 6. t + ∆t and then use equation (6) for θ = 0: An exponential decay behavior can be noticed. L−1 2 −1 L (0)(t + t x) = (t + t) 2 x+s∆x,0 ρL ∆ , ∑ (ρ0 ∆ )2Lx+s∆x,0 (17) s=0 L−1 2 −1 L = (t) 2 x−∆x+s∆x,0 ∑ (ρ0 )2Lx−∆x+s∆x,0 , s=0 L−1 2 −1 L (1)(t + t x) = (t + t) 2 x+s∆x,1 ρL ∆ , ∑ (ρ0 ∆ )2Lx+s∆x,1 (18) s=0 L 2 −1 L = (t) 2 x+∆x+s∆x,1 ∑ (ρ0 )2Lx+∆x+s∆x,1 , s=2L−1 FIG. 3. The coarse-graining procedure for an arbitrary level of (i) where we adopted the notation ρL (t + ∆t,x) = L and time step h∆t. While at the level L = 0 the time evolution x,i operator, E , is given by a unitary map, at the level L the time (ρL(t))x,i for i = 0,1 to describe the the density ma- evolution is given by a quantum channel E˜ derived from E . trix elements in the level L. The challenge here is to see whether in the right-hand side of equations (17) and (18), there is a density matrix in the level L at time t which is derived from the a density matrix of level 0 at equation (5) and it is well known that such evolution time t. It turns out that the diagonal elements of the den- models a free evolution of Dirac particles, as we recalled sity matrix of level 0 after we coarse-grained L times are in the first section. Now, what can we say about the (i) i L dynamics at a generic level L? We certainly know that ρL (t,x − (−1 )∆x/2 ), and therefore, the QCA will lose coherence at each application of the (0) (0) L coarse-graining map. For sake of simplicity, we start ρL (t + ∆t,x) = ρL (t,x − ∆x/2 ) (19) (1) (1) restricting ourselves to θ = 0, the massless case. This ρ (t + ∆t,x) = ρ (t,x + ∆x/2L). means, from equation (6), that ψ0 and ψ1 are fully de- L L coupled: left(right)-moving modes at level L = 0 maps Now, Taylor expands the above equation around x and to left(right)-moving modes at level L ≥ 1. Thus, let us t up to a first order, and the following continuous differ- focus on the population terms. Before we coarse-grain ential equations can be established the density operator at time t, we add a new information (0) c (0) about the initial condition. Here we assume that we start ∂t ρ (t,x) = − ∂xρ (t,x) (20) L 2L L with a superposed state in the center of the lattice of size ( ) c ( ) L ∂ ρ 1 (t,x) = ∂ ρ 1 (t,x). 2 + 1, that is for l = 0,1,...,L − 1 t L 2L x L

2L−1−l ,0 (0) where c = ∆x/∆t. From the results above, it is clear that L (ρ0(0))2L−1−l ,0 = αl , (14) the factor 1/2 starts to create a stationary state as we L−1+l ( ) increase the level L, see ﬁgures (4). In these plots, we (ρ (0))2 ,1 = α 1 , 0 2L−1+l ,1 l can see the contrast between the light cones after the CG

4 ever, in the level L we rescaled the time from t + h∆t to t + ∆t when h 6= 1. As explained in [8], we do not ac- cept h > s since we want to preserve the local structure in the emergent dynamics which comes from the light cone existent in the QCA evolution. Let us go back to our previous problem, given by equation (20), where the localization starts to appear. Com- bining the temporal with the spatial CG, this no physical situation generated by the particle localization can be avoided. In fact, we can keep with the same light cone structure if before each CG operation we compensate it by applying the transition function twice. Thus, since we are moving from the level 0 straight to the level L, which is equivalent of we group 2L cells, the velocity can be kept constant if the compensation by h = 2L is done be- (a) L L fore ΛCG. In other words, we have t + 2 ∆t in the level 0 equivalent to t + ∆t in the level L,

2L−1−1 2Lx+s x,0 ( (t + t))x,0 = (t + 2L t) ∆ ρL ∆ x,0 ∑ ρ0 ∆ 2Lx+s∆x,0 s=0 2L−1 2Lx+s x,1 ( (t + t))x,1 = (t + 2L t) ∆ . ρL ∆ x,1 ∑ ρ0 ∆ 2Lx+s∆x,1 s=2L−1

Therefore, x ± ∆x/2L, goes to x ± ∆x, when we have t + 2L∆t in the level 0 instead t + ∆t, and consequently the factor 1/2L in equation (20) is removed, that is

(b) (0) (0) ∂t ρ (t,x) = −c∂xρ (t,x) (21) FIG. 4. Light cone propagation before (a) and after (b) two L L (1) (1) space coarse-grained transformation. The particle propagation ∂t ρL (t,x) = c∂xρL (t,x). becomes slower at each CG transformation, where the speed of the particle is rescaled by the factor 1/2L. Before we move to the case where m 6= 0, here we can map generated from a massless particle which started in already point out that the emergent dynamics derived in a superposition state in x = 0 at t = 0. In other words, by the classical limit can be easily described by the parti- doing only CG in space, trivial dynamics start to appear tioned cellular automata [8, 26, 30] with two subcells, fast. Inspired by the results given in [8] more interesting which is equivalent to two bits at each cell in the clas- cases can be derived when the coarse-graining in time is sical case. The transition function is composed of two permutation operators. Both permutations are the SWAP also done. In this procedure, the transition function, E , can be applied h times before each CG application. The operator, analogous to the PUQCA used to describe the massless case, that is, the first one acting locally inside number of E applications satisfies the following upper bound h ≤ s, where s is the amount of grouped subcells all vertices, moving the excitation from the subcell 0 (1) to the he subcell 1 (0), and the second interacting the where the local maps Λx are applied. In our case, from equation (7), s = 2. vertices, sending the excitation from the vertex x, in the subcell 1 (0) to x + 1 (x − 1) in the subcell 0 (1). In figure (3) we elucidated the procedure with an arbitrary value for h. Notice that in the new level the update The starting point for the massive case is the recur- state is achieved by only one application of E˜, i.e., the rence relation for the probabilities. From equation (6) quantum channel emerged from the CG process. How- when sinθ ≈ θ, in terms of the density matrix elements,

5 x,1 we have (ρL(t + ∆t))x,1 2L−1 x,0 x−∆x,0 h 2Lx−∆x+s∆x,0 2Lx+∆x+s∆x,1 ρ(t + ∆t)x,0 ≈ ρ(t)x− x,0 (22) ≈ (t) + (t) ∆ ∑ (ρ0 )2Lx−∆x+s∆x,0 (ρ0 )2Lx+∆x+s∆x,1 2 s=2L−1 mc ∆t x−∆x,0 x−∆x,1 + i ρ(t)x− x − ρ(t)x− x 2 h¯ ∆ ,1 ∆ ,0 mc ∆t 2Lx−∆x+s∆x,0 2Lx−∆x+s∆x,1 + i (ρ0(t)) L − (ρ0(t)) L x,1 x+∆x,1 h¯ 2 x−∆x+s∆x,1 2 x−∆x+s∆x,0 ρ(t + ∆t)x,1 ≈ ρ(t)x+∆x,1 2Lx+∆x+s∆x,1 2Lx+∆x+s∆x,0i 2 + (ρ (t)) − (ρ (t)) . mc ∆t 0 2Lx+∆x+s∆x,0 0 2Lx+∆x+s∆x,1 + i ρ(t)x+∆x,1 − ρ(t)x+∆x,0 , h¯ x+∆x,0 x+∆x,1 (25)

Notice that in the pair of equations above the second- In this stage, we can already do a Taylor expansion in order terms in time, i.e., ∆t2, are not included. Now, we both sides. On the left side, we expand until the first will not assume anything about the initial condition, and order in time. However, on the right side there is a sub- thus, in the level L we have tlety that should be noticed. First, the expansion is done around the point 2Lx + p, where p = ±s∆x. Second, 2L−1−1 h L L i while for the terms without the mass the expansion is (t) x,0 = (t) 2 x+s∆x,0 + (t) 2 x+s∆x,1 (ρL )x,0 ∑ (ρ0 )2Lx+s∆x,0 (ρ0 )2Lx+s∆x,1 done until the first order, that is s=0 L− 2 1−1 2Lx−∆x+s∆x,0 2Lx+∆x+s∆x,1 2Lx+s∆x (t) + (t) (ρ0 )2Lx−∆x+s∆x,0 (ρ0 )2Lx+∆x+s∆x,1 (26) = ∑ (ρ0(t))2Lx+s∆x , s=0 2Lx+s∆x ∆x 2Lx+s∆x,0 ≈ (ρ0(t)) L − ∂x (ρ0(t)) L 2L−1 2 x+s∆x L 2 x+s∆x,0 h L L i 2 ( (t))x,1 = ( (t))2 x+s∆x,0 + ( (t))2 x+s∆x,1 ρL x,1 ∑ ρ0 2Lx+s∆x,0 ρ0 2Lx+s∆x,1 , ∆x 2Lx+s∆x,1 L−1 + ∂ (ρ (t)) , s=2 2L x 0 2Lx+s∆x,1 2L−1 2Lx+s∆x = ∑ (ρ0(t))2Lx+s∆x , the ones where the mass appears we only consider the s=2L−1 zero order in the expansion. This difference in the order considered during the Taylor expansion appears because the ones which includes the mass are already first wherein the dependence of the internal state of the par- order terms in ∆t. However, it turns out that the result ticle was eliminated in the right part of the equality of adding all zero orders terms eliminates the mass in the above. This elimination can be done since the prob- emergent dynamics, since ability of finding a particle in the cell x is given by ρ(t)x,0 + ρ(t)x,1. If we also add both equations above, 2Lx−∆x+s∆x,0 2Lx+∆x+s∆x,0 ∆x x,0 x,1 (ρ0(t)) L − (ρ0(t)) L ≈ O . we get 2 x−∆x+s∆x,1 2 x+∆x+s∆x,1 2L (27) L Thus, the following dynamics is derived for the massive 2 −1 L x 2 x+s∆x case (ρL(t))x = ∑ (ρ0(t))2Lx+s∆x . (23) s=0 x ∂t (ρL(t)) Likewise the massless case, our next step is to look up x 2L−1 our density matrix at time t + ∆t and then making use of c 2Lx+s∆x,1 2Lx+s∆x,0 = ∂ (ρ (t)) − (ρ (t)) . equations (22). Thus, the update state at level L writing ∑ L x 0 2Lx+s∆x,1 0 2Lx+s∆x,0 s=0 2 in terms of the level 0 is (28) x,0 (ρL(t + ∆t))x,0 Therefore, all the interference terms that would appear as 2L−1−1 a consequence of the mass are not detectable in the level h L L ≈ (t) 2 x−∆x+s∆x,0 + (t) 2 x+∆x+s∆x,1 L. Moreover, since in equation (28) the dynamics of the ∑ (ρ0 )2Lx−∆x+s∆x,0 (ρ0 )2Lx+∆x+s∆x,1 s=0 level L is level 0 dependent, there is not a model in the 2 mc ∆t L L level L which is derived from the level 0 when the mass + i (ρ (t))2 x−∆x+s∆x,0 − (ρ (t))2 x−∆x+s∆x,1 h¯ 0 2Lx−∆x+s∆x,1 0 2Lx−∆x+s∆x,0 of the particle is included. L L i However, in the level L the internal degree of the par- + ( (t))2 x+∆x+s∆x,1 − ( (t))2 x+∆x+s∆x,0 , ρ0 2Lx+∆x+s∆x,0 ρ0 2Lx+∆x+s∆x,1 ticle is not perceived anymore, and then, since the states (24) 0 and 1 belong to the same cell but located in different

6 [H,ρ]. We used the partitioned unitary quantum cellular automata [7] as the computational model to derive exp(iHD). Inspired by [8], we provided an analogous coarse-graining scheme to explore emergent dynamics, but now for its quantum version, namely the PUQCA. We presented a numerical instance and analytical ev- idence that the coarse-graining map employed here gen- erates decoherence in our quantum system. As a consequence of multiple CG map applications, a semiclassical system can be achieved. Furthermore, we explicitly show that under a particular initial condition for the massless case, the classical transport given by ∂t ρL + vσz∂xρL = 0 FIG. 5. Variance against time from the QCA at level 0 and the is the emergent dynamics by computing the continuous emergent dynamics derived at level 3. In both cases, the same limit. In this emergent system, the velocity v can either value θ = 0.2 is used. go exponentially faster to zero, i.e., v ∼ 1/2L, or remain the same if the temporal CG is also included. On the other hand, when the mass is included and no assump- points, as shown in ﬁgure (1b), we can say that the dis- tions about the initial conditions are done, the diffusion placement between them is ∆x as well. Therefore, by equation is recovered. It is well known that the diffusion setting equation is achieved by computing the continuous limit

x,1 x of the random walk. Therefore, the CG procedure devel- (ρ0(t))x,1 = (ρ0(t))x , (29) oped here provided the right tools to see how the random walk can be derived from the quantum walk. the following equality holds Differently from [8] that shows a general procedure to study emergent dynamics, here just one map and one t x,0 t x−∆x (ρ0( ))x,0 = (ρ0( ))x−∆x . (30) type unitary operators are used. However, our results are pointing out that this procedure can be generalized for Using this equivalence into equation (28), and consider- different complete positive trace-preserving maps as well ing that we have x2 t = λ, we get ∆ /∆ as for different dynamics, including the cases where we 2L−1 2 have interacting quantum walks [3,6, 22], and we in- ∆x 2 2Lx+s∆x tend to address these possibilities in future publications. ∂t ρL(t,x) = ∑ L ∂x (ρ0(t))2Lx+s∆x (31) s=0 2 ∆t Another question that immediately raises from our pro- 2 2L−1 ! cedure is whether we can derive emergent local quan- ∆x L 2 2 x+s∆x tum channels E˜ arising from QCA, following the proce- = L ∂x ∑ (ρ0(t))2Lx+s∆x 2 ∆t s=0 dure shown in [11, 12]. If such quantum channels can λ be achieved, we will provide a powerful tool to generate = ∂ 2ρ (t,x). 2L x L transition functions such that noisy quantum cellular automata, as described in [2], naturally appear during the which is the diffusion equation. decoherence process. In ﬁgure (5) we contrast the variance between the level Therefore, this formalism might be an excellent can- 0, that simulates the Dirac equation, and the emerged didate to study open quantum systems and also to shed one. In the level 0, we can see the character linear of the some light on more fundamental questions in the quan- variance, as expected from the quantum walks dynam- tum theory, e.g., a tool to explore quantum-to-classical ics. However, in the level 3 we can notice a nonlinear transition via quantum cellular automata. Finally, an- behavior, similar to the dispersion relation provided by other possible direction to go is contrasting the emergent the random walk. Such result shown in ﬁgure (5) sup- dynamics derived here against alternative formulation to ports the continuous limit derived in equation (31). study different coarse levels like the use of the wavelets basis [4].

V. CONCLUSION ACKNOWLEDGEMENTS In this work, we studied emergent dynamics starting from the Dirac like Hamiltonian, i.e., HD = σx ⊗ I am really grateful to Fernando de Melo who in- px + mσz ⊗ I, where the dynamics is driven by ∂t ρ = troduced the subject of coarse-graining and proposed

7 the topic of emergent dynamics via quantum cellular manuscript, and Gavin Brennen for helpful comments. automata during my PhD I also would like to thank This work was partially supported by ARC Discovery Giuseppe Di Molfetta for helpful discussion during this Projects. project, Yuval R. Sanders for careful reading of the

[1] R. Alicki, M. Fannes, and M. Pogorzelska. Quantum gen- 2014. eralized subsystems. Phys. Rev. A, 79:052111, May 2009. [16] Oleg Kabernik. Quantum coarse graining, symmetries, [2] Michele Avalle and Alessio Serafini. Noisy quantum and reducibility of dynamics. Phys. Rev. A, 97:052130, cellular automata for quantum versus classical excitation May 2018. transfer. Phys. Rev. Lett., 112:170403, May 2014. [17] Leo P. Kadanoff. Scaling laws for ising models near. [3] Scott D. Berry and Jingbo B. Wang. Two-particle quan- Physics Physique Fizika, 2:263–272, Jun 1966. tum walks: Entanglement and graph isomorphism testing. [18] Johannes Kofler and Caslavˇ Brukner. Conditions for Phys. Rev. A, 83:042317, Apr 2011. quantum violation of macroscopic realism. Phys. Rev. [4] Gavin K. Brennen, Peter Rohde, Barry C. Sanders, and Lett., 101:090403, Aug 2008. Sukhwinder Singh. Multiscale quantum simulation of [19] N. D. Mermin. Quantum mechanics vs local realism near quantum field theory using wavelets. Phys. Rev. A, the classical limit: A bell inequality for spin s. Phys. Rev. 92:032315, Sep 2015. D, 22:356–361, Jul 1980. [5] Hong-Bin Chen, Ping-Yuan Lo, Clemens Gneiting, Joon- [20] Leonard Mlodinow and Todd A. Brun. Discrete space- woo Bae, Yueh-Nan Chen, and Franco Nori. Quantifying time, quantum walks, and relativistic wave equations. the nonclassicality of pure dephasing. Nature communi- Phys. Rev. A, 97:042131, Apr 2018. cations, 10(1):1–9, 2019. [21] David Poulin. Macroscopic observables. Phys. Rev. A, [6] Pedro C. S. Costa, Fernando de Melo, and Renato Portu- 71:022102, Feb 2005. gal. Multiparticle quantum walk with a gaslike interac- [22] Xizhou Qin, Yongguan Ke, Xiwen Guan, Zhibing Li, tion. Phys. Rev. A, 100:042320, Oct 2019. Natan Andrei, and Chaohong Lee. Quantum walks of [7] Pedro C. S. Costa, Renato Portugal, and Fernando two interacting particles in one dimension. arXiv preprint de Melo. Quantum walks via quantum cellular automata. arXiv:1402.3349, 2014. Quantum Information Processing, 17(9):226, Jul 2018. [23] Maximilian Schlosshauer. The quantum-to-classical tran- [8] Pedro C.S. Costa and Fernando de Melo. Coarse grain- sition and decoherence. arXiv preprint arXiv:1404.2635, ing of partitioned cellular automata. Journal of Cellular 2014. Automata, 15:305–331. [24] Pedro Silva Correia, Paola Concha Obando, Raúl O Valle- [9] Paul Adrien Maurice Dirac. The quantum theory of the jos, and Fernando de Melo. Macro-to-micro quantum electron. Proceedings of the Royal Society of London. Se- mapping and the emergence of nonlinearity. arXiv e- ries A, Containing Papers of a Mathematical and Physical prints, pages arXiv–2007, 2020. Character, 117(778):610–624, 1928. [25] Pedro Silva Correia and Fernando de Melo. Spin- [10] Cristhiano Duarte. Compatibility between agents as a tool entanglement wave in a coarse-grained optical lattice. for coarse-grained descriptions of quantum systems. J. Phys. Rev. A, 100:022334, Aug 2019. Phys. A, 53(39):395301, aug 2020. [26] Tommaso Toffoli and Norman Margolous. Cellular Au- [11] Cristhiano Duarte, Barbara Amaral, Marcelo Terra tomata Machines. MIT Press Series in Scientific Compu- Cunha, and Matthew Leifer. Investigating coarse- tation, 1985. grainings and emergent quantum dynamics with [27] Isadora Veeren and Fernando de Melo. Entropic uncer- four mathematical perspectives. arXiv preprint tainty relations and the quantum-to-classical transition. arXiv:2011.10349, 2020. Phys. Rev. A, 102:022205, Aug 2020. [12] Cristhiano Duarte, Gabriel Dias Carvalho, Nadja K. [28] Tian Wang, Roohollah Ghobadi, Sadegh Raeisi, and Bernardes, and Fernando de Melo. Emerging dynamics Christoph Simon. Precision requirements for observing arising from coarse-grained quantum systems. Phys. Rev. macroscopic quantum effects. Phys. Rev. A, 88:062114, A, 96:032113, Sep 2017. Dec 2013. [13] Paul Ehrenfest and Tatiana Ehrenfest. The conceptual [29] Michael M Wolf. Quantum channels & operations: foundations of the statistical approach in mechanics. Guided tour. Lecture notes available at http://www-m5. Courier Corporation, 1990. ma. tum. de/foswiki/pub M, 5, 2012. [14] Navot Israeli and Nigel Goldenfeld. Coarse-graining of [30] Dieter A. Wolf-Gladrow. Lattice-Gas Cellular Au- cellular automata, emergence, and the predictability of tomata and Lattice Boltzmann Models: An Introduction. complex systems. Physical Review E, 2006. Springer, Berlin, 2000. [15] Hyunseok Jeong, Youngrong Lim, and M. S. Kim. Coarsening measurement references and the quantum-to- classical transition. Phys. Rev. Lett., 112:010402, Jan

8 VI. APPENDICES and then, from equation (9), we see that the CPTP map covering the coherence terms leads us to

A. Effective state at level L + 1 from the state at level L 0 0 ΛCG |k,ai k ,a 0 1 k k 0 In the following, we apply the coarse-graining local = ,k mod 2 ,k mod 2 . (34) map to the operator (10), getting a general expression for 3 2 2 the coarse-grained density operator in level L + 1. To get this general form, note that we can distinguish three Finally, recast together equations (32), (33) and (34), we different cases: (i) x = x0 and a = a0, (ii) a 6= a0, for all get the ﬁnal expression shown in equation (11), which is for x = x, and (iii) x 6= x . In case (i), the state at level the coarse grained density operator at level L + 1. L + 1 reads: B. Decoherence by coarse-graining ΛCG (|k,aihk,a|)

= ···(0)k−1 (1)k (0)k+1 ··· ···(0)k−1 (1)k (0)k+1 ··· . Lemma 1. Let p ∈ [0,1]D, with D ≥ 2, be a probability D distribution, i.e., ∑ pi = 1. Then the inequality Note that the three conditions for x shown above imply i=1 the same for k. The number of subcells for each position D √ √ is reduced to one, and we have to space group subcells at ∑ pi ≤ D (35) different spatial positions, to build a new cell hosting two i=1 qubits. Consequently, the choice of the space grouping of the subcells, at a higher level, might be determinant holds, and it is saturated when pi = 1/D for all i ∈ [D]. and may depend on the spatial position at level L. In Proof. The proof goes by the way of induction. For the fact, applying the map at even position x leads to: base case, D = 2, we deﬁne the function f2 : [0,1] → R as ED k k ···(1) k ··· ···(1) k ··· = ,0 ,0 . √ p 2 2 2 2 f2(p1) = p1 + 1 − p1;

where the one particle excitation takes now this place at where we already employed the normalization property. a Setting d f dp = 0 is simple to obtain that the max- = 0, the left subcell. While if the local map is applied 2/ 1 √ at odd position x, the single-particle excitation will be imum is obtained for p = 1/2. Therefore, p + √ √ 1 1 placed in the right subcell a = 1: 1 − p1 ≤ 2. D−1 √ For the inductive step, we assume that ∑ pi ≤ √ i=0 ED k − 1 k − 1 D − 1. For D we have: ···(1) k−1 ··· ···(1) k−1 ··· = ,1 ,1 . 2 2 2 2 D √ D−1 √ √ Notice once again that, without lack of generality, we ∑ pi = ∑ pi + pD space-group from the left to the right. i=0 i=0 v   Both the above cases can be recast as follows: uD−1 D−1 s p √ = u p i + p . (36) t ∑ j  ∑ D−1  D k k j=0 i=0 ∑ j=0 p j ΛCG (|k,aihk,a|) = ,k mod 2ih ,k mod 2 (32) 2 2 0 q D−1 0 Noting that deﬁning p = pi/∑ p j, then p is a Case (ii), i.e., a 6= a0 and x = x0 leads directly to i j=0 probability distribution in [0,1]D−1, for which the in- 0 D−1 ΛCG |k,ai k,a = 0, (33) duction hypothesis can be applied, and that ∑ j=0 p j = 1 − pD we have: because Λx (|01ih10|) = 0. D Finally, for x 6= x0, we follow the same procedure, and √ p √ √ ∑ pi ≤ 1 − pD D − 1 + pD. (37) we apply the local map for each spatial position x: i=0 √ √ √ |k,ai k0,a0 = f (p ) = − p D − + p ΛCG Finally, the function D √D 1 D 1 D E has its maximum value, D, at pD = 1/D, what ﬁnishes ...(1) ,...(0) 0 ,...... ,(0) ,...(1) 0 ,... , k xL k k the proof.

9 n/2L n/2L Before we proof that the off-diagonal elements of the Theorem 2. Let (ρL) ∈ C × C be the density density matrix go exponentially faster to zero, the sim- matrix in the coarse level L, with n 2L, where their plest case of the upper bound value from the level 0 to off-diagonal elements in terms of the level 0 are given by the level 1 is explicitly computed. equation (13). Then, the absolute value of the coherence 0 From equation (13), for b = b = 0, we see that when terms of (ρL) is upper bounded by we move from the level 0 to 1 we have, √ !L x,0 1 h 2x,0 2x,1 2 2 (ρ ) = (ρ ) 0 + (ρ ) 0 (38) , (44) 1 x,0 3 0 2x ,0 0 2x ,1 3 2x,1 2x,0 i + (ρ0)2x0,0 + (ρ0)2x0,1 . for x 6= x0 and b,b0 = 0,1. Now, from the triangle inequality, i.e., |a + b| ≤ |a|+|b|, we can write Proof. From equation (13), the following inequality can be written as 1 h x,0 ≤ 2x,0 2x,1 (ρ1)x0,0 (ρ0)2x0,0 + (ρ0)2x0,1 (39) 3 2L−1−1 L L−1 2x,1 2x,0 i x,b 1 2 x+2 b∆x+s∆x,0 (ρ ) 0 0 ≤ (ρ ) + (ρ0)2x0,0 + (ρ0)2x0,1 . L x ,b L ∑ 0 2Lx0+2L−1b0∆x+s0∆x,0 3 s,s0=0 L L−1 L L−1 Using the following inequality for the density matrix + (ρ )2 x+2 b∆x+s∆x,1 (ρ )2 x+2 b∆x+s∆x,0 0 2Lx0+2L−1b0∆x+s0∆x,0 0 2Lx0+2L−1b0∆x+s0∆x,1 2 Lx+ L−1b x+s x r,a r,a r0,a0 2 2 ∆ ∆ ,1 ≤ + (ρ0) L 0 L−1 0 0 . (45) ρr0,a0 ρr,a ρr0,a0 , (40) 2 x +2 b ∆x+s ∆x,1

equation (39) can be written as Now, using the inequality given by equation (40)

1 q 0 q 0 (ρ )x,0 ≤ (ρ )2x,0 (ρ )2x ,0 + (ρ )2x,1 (ρ )2x ,1 1 x0,0 3 0 2x,0 0 2x0,0 0 2x,1 0 2x0,1 2L−1−1 x,b 1 q q (ρ ) 0 0 ≤ 2x,1 2x0,0 2x,0 2x0,1 L x ,b L ∑ + (ρ0)2x,1 (ρ0)2x0,0 + (ρ0)2x,0 (ρ0)2x0,1 . 3 s=0 " 2L−1−1 q L L−1 q L 0 L−1 0 0 0 After a simple algebraic manipulation in the equation × ( )2 x+2 b∆x+s∆x,0 ( )2 x +2 b ∆x +s∆x ,0 ρ0 2Lx+2L−1b∆x+s∆x,0 ∑ ρ0 2Lx0+2L−1b0∆x0+s0∆x,0 above, we have s0=0 q q q q 2Lx0+2L−1b0 x0+s x0,1 x,0 1 2x,0 2x0,0 2x0,1 + (ρ ) ∆ ∆ (ρ ) 0 ≤ (ρ ) (ρ ) 0 + (ρ ) 0 0 2Lx0+2L−1b0 x0+s0 x,1 1 x ,0 3 0 2x,0 0 2x ,0 0 2x ,1 ∆ ∆ q q q q 2L−1−1 q 1 2x,1 2x0,0 2x0,1 2Lx+2L−1b∆x+s∆x,1 2Lx0+2L−1b0∆x0+s∆x0,0 + (ρ ) (ρ ) 0 + (ρ ) 0 . + (ρ ) (ρ ) 3 0 2x,1 0 2x ,0 0 2x ,1 0 2Lx+2L−1b∆x+s∆x,1 ∑ 0 2Lx0+2L−1b0∆x0+s0∆x,0 s0=0

q L 0 L−1 0 0 0 Noticing that + ( )2 x +2 b ∆x +s∆x ,1 . ρ0 2Lx0+2L−1b0∆x0+s0∆x,1 q 2x,a q 2x0,0 q 2x0,1 q 2x0,0 q 2x0,1 (ρ0)2x,a (ρ0)2x0,0 + (ρ0)2x0,1 ≤ (ρ0)2x0,0 + (ρ0)2x0,1, (41) Now we can apply the inequality used in equation (41) 2x,a since (ρ0)2x,a ≤ 1 for a = 0,1 we have that to get

q q L−1 x,0 2 2x0,0 2x0,1 L−1 "2 −1 q (ρ1) 0 ≤ (ρ0) 0 + (ρ0) 0 . (42) x,b 2 2Lx0+2L−1b0∆x0+s∆x0,0 x ,0 2x ,0 2x ,1 (ρ ) 0 0 ≤ (ρ ) 3 L x ,b L ∑ 0 2Lx0+2L−1b0∆x0+s0∆x,0 3 s0=0

Now, from lemma1 q L 0 L−1 0 0 0 + ( )2 x +2 b ∆x +s∆x ,1 ρ0 2Lx0+2L−1b0∆x0+s0∆x,1 q 0 q 0 √ 2x ,0 2x ,1 ≤ (ρ0)2x0,0 + (ρ0)2x0,1 2, 2L−1−1 q L 0 L−1 0 0 0 + ( )2 x +2 b ∆x +s∆x ,0 ∑ ρ0 2Lx0+2L−1b0∆x0+s0∆x,0 therefore, s0=0 √ q 2Lx0+2L−1b0∆x0+s∆x0,1 x,0 2 2 + (ρ0)2Lx0+2L−1b0 x0+s0 x,1 , (ρ ) 0 ≤ < 1. (43) ∆ ∆ 1 x ,0 3

10 or The proof can be conclude by using the lemma1 in the inequality above for D = 2L L 2L−1−1 q x,b 2 2Lx0+2L−1b0∆x0+s∆x0,0 (ρ ) 0 0 ≤ (ρ ) L x ,b L ∑ 0 2Lx0+2L−1b0∆x0+s0∆x,0 3 s0=0

q L 0 L−1 0 0 0 + ( )2 x +2 b ∆x +s∆x ,1 . (46) ρ0 2Lx0+2L−1b0∆x0+s0∆x,1