Wigner formalism for a particle on an infinite lattice: dynamics and spin

The recently proposed Wigner function for a particle in an infinite lattice [NJP 14, 103009 (2012)] is extended here to include an internal degree of freedom, as spin. The formalism is developed to account for dynamical processes, with or without decoherence. We show explicit solutions for the case of Hamiltonian evolution under a position-dependent potential, and for evolution governed by a master equation under some simple models of decoherence. Discrete processes are also discussed. Finally we discuss the possibility of introducing a negativity concept for the Wigner function in the case in which the spin degree of freedom is included.


I. INTRODUCTION
Since its introduction, the Wigner function (WF) [1] has played an important role in physics. Quantum mechanics can be entirely formulated using this tool, therefore providing an alternative description of quantum phenomena, along with their dynamics. Also from a more experimental perspective, the WF has proven instrumental for tomographic reconstruction of the states prepared in the lab. The WF is in hence completely equivalent to the standard quantum mechanical formalism. Nevertheless, the particular features of the phase space description make it advantageous in some situations, for instance recognizing the quantum features of states, or dealing with decoherence scenarios. In the WF interference effects manifest in a clear way [2][3][4][5]. Another interesting property that manifests in the visualization of the WF of some states is the appearance of negative values over the phase space. This fact has been considered as a direct manifestation of the quantum nature of such states, and used to characterize their quantumness [6][7][8]. The relativistic extension of the Wigner function [9] has also found applications to a wide variety of problems, ranging from general relativistic kinetic theory and statistical mechanics [10,11], nuclear matter at high densities and temperatures [12], electrons in magnetic fields [13,14], the quark-gluon plasma [15], to neutrino propagation in astrophysical or cosmological scenarios [16,17].
The applications mentioned above make use of a Wigner function defined in continuous space. It is nevertheless possible to introduce also a sensible Wigner function for systems on a discrete space. The definition for the case of a finite dimensional Hilbert space can be traced back to Stratonovich and Agarwal [18,19] (see also [20]), who introduced a spherical, continuous phase space for a spin particle. A possible generalization was proposed by Wootters in 1987 [21] for prime dimensional systems, and later generalized to any power of primes in [22]. A different construction was followed in [23][24][25] which could cope with any dimension of the Hilbert space at the expense of enlarging the size of the phase space grid (see [26,27] for a review). The discrete WF for a finite dimensional system is furthermore related to quantum information problems [7,25,[27][28][29][30][31].
If the discrete Hilbert space is infinite dimensional, a different extension of the WF is required. In [32] we proposed a definition of the WF that can be used for such systems, having the correct marginal properties and with the advantage that a closed form can be obtain in some cases, such as the Gaussian states. Notice that, in contrast to the continuous case, where the axiomatic definition of the WF uniquely determines its functional form [33], in the discrete case different definitions are possible that respect the mathematical conditions enumerated above (see also [34,35] for alternative, related definitions, motivated by the study of the angle and angular momentum phase space).
Many of the problems where the continuous WF has found application concern particles with spin, or with spinor descriptions of quantum fields. In order to use the phase space formalism in this scenario, a generalization has to be introduced which combines the spin and spatial degrees of freedom (dof). One of the most common prescriptions in the literature is the use of a matrix valued WF [36], where the spinor or spin indices give rise to various matrix elements. Indeed, other possibilities exist, such as introducing a phase space for the spin degrees of freedom, which correspond to another discrete, finite dimensional Hilbert space, and construct a real valued WF for the cartesian product of spin and space phase spaces. In the matrix-valued WF, the treatment of space and spin dof is not symmetric. The spatial part is described in terms of a phase space, while the spin is unchanged. Although the treatment is asymmetric, such description has some advantages when dealing with a particle subject to a spin-dependent force, since some effects like the spin precession, or motion that depends on the spin component, are better visualized with respect to a fixed spin basis. Examples of this description are the analysis of the Stern-Gerlach experiment [37], the study of entangled vibronic quantum states of a trapped atom [36], or the reconstruction of the full entangled quantum state for the cyclotron and spin degrees of freedom of an electron in a Penning trap [38].
In this paper, we have extended the definition of the Wigner function introduced in [32] to incorporate the spin of a particle, using the Wigner matrix formalism for the spin degrees of freedom, and we illustrate the consequences of this definition by analyzing some simple physical situations, such as states involving spatial and spin entanglement or dynamical evolution, as it appears for a particle subject to a spin-dependent force. The rest of this paper is organized as follows. In Section II we introduce a definition for the Wigner matrix (WM) that incorporates the spin of the particle, and summarize the main properties that are satisfied by this object. To illustrate the structure of this representation, we consider some simple cases in Section III. Section IV contains the main results of our paper, concerning the dynamics obeyed by the WM under the influence of an interacting Hamiltonian that may depend or not on the spin. First, we study the time evolution in continuous time, by deriving the equation of motion for the WM and solving this equation in some simple cases. The situation without spin serves us to consider the special case of a particle on a lattice interacting with a linear potential. We also investigate the interaction that appears for a spin-dependent force to visualize the main differences with the spinless case. Finally, we study the effect of decoherence for the system under consideration. Also in this Section, we show how one can make use of the WM to investigate the dynamics that appear in some discrete-time problems, and consider the particular example of the quantum walk. As before, we show the effect that decoherence may have on such problems.
One of the advantages of a WF description of continuous variable systems is the access to a negativity that measures the non-classicality of states. Although the relation of the negativity to non-classicality is well established, this quantity does not correspond to a physical observable. With a more general definition as the WM and the occurrence of (nonclassical) spin degrees of freedom, we may wonder if there is a generalized negativity quantity and whether it retains some physical information. This is discussed in Section V. Section VI presents our main conclusions. The derivation of some formulae has been relegated to the Appendix in order to make our presentation more transparent.

II. PARTICLE WITH SPIN ON A ONE-DIMENSIONAL LATTICE
We are interested in the phase space description of a spin 1/2 particle that is allowed to move on an infinite 1D lattice. A paradigmatic example is the quantum walk (QW) on the line, where a particle moves along the sites of a 1D lattice. In its discrete-time version [39], the direction of motion is dictated by the state of an extra two-dimensional Hilbert space (the coin), that can correspond to the internal spin of the moving particle. In fact, during the process the spatial and internal states become entangled, even if the initial state was separable, thus making clear the need for a joint description of both degrees of freedom. Another example is the study of spin dependent transport properties of single atoms in a 1D optical lattice [40].
We will start with the definition of the WF for a (spinless) particle on a 1D lattice already introduced in [32]. We consider a lattice with sites {na/n ∈ Z}, where a is the lattice spacing. To these sites one can associate a basis {|n }, with n ∈ Z. By a Fourier transformation we define a quasi-momentum basis, |q = a 2π n e iqna |n , which can be restricted to the first Brillouin zone, q ∈ [− π a , π a [. The phase space is defined by points (m, k), where m ∈ Z, whereas k is continuous and periodic, taking values in [−π, π[. With these notations, we define the WF as where ρ is the density operator corresponding to the state of the system, and A(m, k) = 1 2π n |m − n n|e −i(2n−m)k are the phase point operators for the lattice. It can be checked that the above definition fulfills the necessary requirements to be considered a valid WF. We refer the reader to the above reference for more information about the properties obeyed by (1).
We now would like to incorporate the additional degree of freedom arising from the spin of the particle. As discussed in the Introduction, there are different approaches in the literature to describe finite dimensional Hilbert spaces, such as the spin of a particle. One can combine both degrees of freedom (spin and lattice) by a tensor multiplication of the corresponding point operators, as done in [41] for angular momentum and spin states.
As discussed in the introduction, here we opt for a prescription with ample acceptance in the continuous applications, namely a matrix-valued WF. A similar choice has been used in relativistic and non-relativistic setups with continuous spatial dof. Among the latter we can mention the study of Stern-Gerlach experiment [37], the analysis of entangled vibronic quantum states of a trapped atom [36], or the reconstruction of the full entangled quantum state for the cyclotron and spin degrees of freedom of an electron in a Penning trap [38]. The Wigner function defined in this way combines the following properties: -It keeps a close analogy with the definition of the relativistic Wigner function [9][10][11], thus allowing to describe the transition from the relativistic to the non relativistic regime.
-It appears as a simple and convenient choice to describe the spin motion in some particular cases, like the Stern-Gerlach experiment in continuous space [37], or the dynamics of a spin 1/2 particle on a lattice under the effect of a spin-dependent force, as described in Sect. IV.
We consider the Hilbert space H = H l ⊗ H s , where H l stands for the motion on the lattice, and H s describes the spin states. The composed Hilbert space is spanned by the basis {|n, α ≡ |n ⊗ |α } with n ∈ Z and |α /α = 0, 1 designate the eigenvectors of the σ z Pauli matrix (these states might also correspond to the computational basis of a qubit, or to the levels of a two level system). According to the above discussion, we propose the following definition for the WM We then have a set of four functions W αβ (m, k), α, β = 0, 1 forming a 2 × 2 matrix. Each function, as before, is defined on the phase space of points (m, k), with m ∈ Z, and k takes values in [−π, π[. A similar definition can be made for any operator O acting on H: Unlike the spatial variables, where the relationship with phase space points is non trivial, there is a direct correspondence between spin indices in the state of the system and indices in the matrix WF. This implies that operations on the spin space, such as rotations, change of basis or interactions with a spin-dependent force, as studied below, become more transparent using the matrix WF than other kind of representations for the spin. Moreover, the definition Eq.
(2) keeps a closer analogy, for pure states, to the relativistic WF used in Quantum Field Theory. For such states one has ρ = |Ψ Ψ| and we can write with Ψ α (n) ≡ n, α |Ψ . In the continuum limit, the functions Ψ α (n) can be interpreted as the components of a Pauli spinor or a Dirac spinor. In this case, Eq. (4) can be related to the relativistic WF already mentioned in the Introduction.
Some of the properties discussed in [32] can be easily generalized for the matrix WF. 1) We have which implies that the matrix WF is Hermitian. The normalization condition becomes 2) Also, 3) Given two operators C, D and their corresponding Wigner 4) A complete knowledge of the WF can be used to reconstruct the density operator ρ: 5) The marginal distributions of (2) are related to matrix elements of the density operator As already discussed in [32], these equations reflect the distinction between the coordinates of the phase space points, m ∈ Z, k ∈ [−π, π[, and the position and quasimomentum bases, n, q. The k coordinate is adimensional and does not directly represent a momentum value, but is connected to q = k/a. The spatial label m in phase-space is only connected to a discrete position, s, for even values, m = 2s, while the odd values of m are analogous to the odd half-integer phase space grid points in [24,25].

III. PARTICULAR CASES
In order to obtain some insight about the characteristics of the matrix WF Eq. (2), we will give the explicit form it takes for some particular cases.

• Product state
We start by considering a product state of spatial and spin degrees of freedom where ρ L represents a general state on the lattice, and ρ S is an arbitrary spin state. In this case, we readily obtain with • Superposition of two deltas Let us consider the WM for the state formed by a superposition of two localized states at lattice sites |n 1 and |n 2 with n 1 = n 2 ∈ Z where α is an arbitrary complex number that represents the relative weight of the state |n 2 . For α = 1 we obtain a Schrödinger-cat state. The corresponding WF can be easily calculated. Written in matrix form in the above spin basis, In this case, the WM is zero everywhere except for three particular values of the space-like phase coordinate, m = 2n 1 , 2n 2 , n 1 + n 2 . It is interesting to compare the structure provided by Eq. (16) with the corresponding superposition of two localized states without spin [32], given by In that case, the WF is a scalar function where φ is the phase of the complex coefficient α, and ∆n = n 2 − n 1 . One observes the different terms in (18) appear distributed on different matrix positions in Eq. (16). In particular, the out of diagonal term in (16) corresponds to the interference, oscillating term in (18). This term plays an interesting role related to the non positivity of the WF. We will return to this point later.
• Superposition of two Gaussian states The superposition of two discretized pure Gaussian states with orthogonal spin components is another interesting state for which the WM defined in this work can be computed analytically. Such a state is defined as for arbitrary a, b ∈ Z, σ ∈ R + . For this state, the WF can be expressed as a matrix in the same {|0 , |1 } basis where with N = θ 3 (0, e − 1 σ 2 ) the normalization constant. The Jacobi theta function θ 3 (z, q) is defined as θ 3 (z, q) ≡ n q n 2 e 2izn for complex arguments q, z, with |q| < 1 [42]. As in the previous example, we find an important difference with the WF for the case without spin [32], since the components in the scalar function appear here distributed as the components of the matrix WF. In the limit a = −b ≫ σ with σ → 0 we recover the result for the two deltas (16) corresponding to the case n 1 = −n 2 = a and α = 1. Figure 1 shows the four components of the WM for a two-Gaussian state, as given by Eqs. (20)(21)(22). One can immediately observe on each component the presence of a secondary image that reflects the property Eq. (7). In [32] we discussed with some detail, for the spinless case, the peculiarities related to this duplicate.

IV. DYNAMICS
The WF formalism can be used, not only to allow for a description of a given state, but also to analyze the dynamics, and to visualize it in phase space. Our purpose is to study the motion of a particle on a lattice in terms of the corresponding WF. We start from the simplest case, which corresponds to the spinless particle, and then move to a more general situation, where the particle interacts with a spin-dependent term. The time evolution will be first considered within continuous time, a situation that can be applied to most problems in physics, and can be described by the Schrödinger equation.
A. Continuous time

Particle without spin
Let us consider a spinless particle moving on a lattice under the influence of a potential V that depends on the lattice site. We concentrate on the following Hamiltonian that appears as a consequence of the tight-binding approximation in crystals, where the parameter J is a characteristic of the system which is related to the hopping probability of an electron to the nearest neighbor, and the displacement operators T ± are defined by T ± |n = |n ± 1 . Notice that the Hamiltonian (23) can also be considered as a discretized version of (with M the mass of the particle) if one defines J = − 1 2Ma 2 . The wave function can be written as ψ(n, t), with t the time, so that the Schrödinger equation with V n ≡ n|V |n . The last term inside the brackets in Eq. (1) can be easily reabsorbed into the definition of the coefficients V n (it can be also understood as a term proportional to the identity in the Hamiltonian, thus contributing only as a position-independent phase as time evolves). Therefore we omit that term. It is straightforward to derive an evolution equation satisfied by the WF for the above problem. We begin with the von Neumann equation for the density operator Making use of (1) one arrives to where we have explicitly showed the time dependence of ρ and W (m, k) for the sake of clarity. Let us consider that V (x) is a continuous and infinitely derivable function. In this case, one can obtain a closed form of the above expression for the WF, as showed in the Appendix. As a result, one arrives to the following expression (28) It must be noticed that Eq. (28) also holds for the WM (2) if we introduce the spin of the particle, by simply replacing W (m, k, t) −→ W αβ (m, k, t), since none of the spatial operations in this equation can affect the spin indices.
Before we go on, we will consider the continuous limit (a → 0) of Eq. (28). In this limit, our WF has to be replaced by the corresponding function W c (x, q, t) following the prescription [32] By replacing J = − 1 2Ma 2 and substituting (29) in (28), and taking the limit (a → 0), one obtains the equation Eq. (30) is the equation of motion for the WF under the effect of an external potential V (x) in continuous space, where q represents the momentum of the particle (ranging from −∞ to ∞) (see, for example [4]).
As an interesting particular case, we will study the case of a linear potential, i.e. V (x) = λx, with λ a real constant. Eq. (28) adopts a simple form To solve this equation, we perform a Fourier transformation on the variable m by introducing the functioñ the new variable q taking values on the interval [−π, π[. With the help of this function, we can rewrite Eq. (32) as ∂ ∂tW (q, k, t) = −4iJ sin k sin qW (q, k, t) + λa ∂ ∂kW (q, k, t).
The change of functionW leads to the following equation for f (q, k, t): which implies that f (q, k, t) must be of the form f (q, k, t) = g(q, k + λat), with g(q, k) an unknown function that can be determined by the initial (t = 0) condition in Eq. (34), giving g(q, k) = e 4iJ cos k sin q λaW (q, k, 0).
We finally obtain, after some algebrã To derive an expression for the WF, we need the inverse relation of Eq. (32), given by and make use of the formula [43] J n (z) = 1 2πˆπ −π dqe −inq e iz sin q , where n ∈ Z, z ∈ C, and J n (z) are the Bessel functions of the first kind. After substituting Eq. (37) into (38) we arrive to the final expression Notice that, in the latter equation, the argument k + λat is to be understood modulo 2π. Using this fact, one can readily obtain that the above solution exhibits a time periodicity which corresponds to the well known phenomenon of Bloch oscillations, that can be observed for electrons confined in a periodic potential (the lattice) subject to a constant force, as for example a constant electric field. The corresponding frequency ω B = |λ| a is precisely what is expected for our linear potential V (x) = λx. Directly related to the above treatment, it appears quite natural to attempt a parallelism with a situation that describes the dynamics of a particle under the effect of a constant gravitational field, V (x) = m g gx, where m g is the gravitational mass and g the acceleration of gravity. Notice that, for the following discussion to make sense, one should design a physical system that is described by this potential, and that Eq. (25) can be considered as a discretized approximation to (24), with J = − 1 2Ma 2 . We will return to this discussion later. We find it convenient to use the symbol m i instead of M to represent the inertial mass, and to recover the Planck constant. We observe that the argument of the Bessel functions in Eq. (40) depends upon the combination where k g ≡ ( 2mimg g 2 ) 1/3 is a characteristic wave vector that modulates the spatial dependence of energy eigenstates in a gravitational field in continuous space [44]. As the authors of this work discuss, this is one of the possible effects for quantum particles under the effect of gravity, where various combinations of (powers of) m g and m i may appear depending on the problem under consideration, thus paving the way to measuring these two quantities independently.
The dynamics on the lattice we just considered offers a similar perspective. The time evolution in Eq. (40) is governed by the product k g a, which involves the lattice spacing as a new parameter, thus allowing an extra degree of freedom in the design of experiments, if they are performed on a lattice instead of in continuous space. However, one has to be careful about this point: Only if the design of the experiment is such that J and V (x) correspond to the above hypothesis, the previous discussion can make sense.
To illustrate the behavior of the WF, we plotted in Fig. 2 several snapshots obtained by evolving an initial Gaussian state of the form (21). The time evolution is governed by Eq. (40). One observes several features on this plot. First, the position of the maximum shows oscillations for the variable m, as corresponding to the Bloch oscillations discussed above, while variable k evolves linearly (and periodically) with time. During the evolution, the WF also experiences a distortion that is similar to the one observed in continuous space [44]. One also observes the presence of a secondary image which manifests as vertical strips.

Particle with spin
We return to the description of a particle with spin 1/2. Our purpose is to analyze the dynamics for such a system, and compare it with the spinless case. To do so, we need to introduce some spin-dependent potential, otherwise the different components in the WM will evolve exactly in the same way, and the results of the previous subsection apply. In order to make this comparison as close as possible, we will consider the time evolution under the effect of a Hamiltonian of the form where V is, as before, a site-dependent scalar potential. It is possible to obtain an evolution equation, similar to (28), when the particle is subject to the above Hamiltonian in the lattice. This derivation is made in the Appendix, the main difference with the spinless case being that the diagonal and off-diagonal components of the WM evolve differently. In what follows, we concentrate on the particular example of a discretized linear potential V n = λan, with λ a real constant. Then, Eq. (72) particularizes to and (valid for α = β). The first equation can be easily solved by comparison to (31). We only have to perform the replacement λ −→ (−1) α λ. Therefore, we can write the solution using the same procedure as in the case with no spin, to obtain The same comments made in the previous section hold here: W αα (m, k, t) is periodic in time, with frequency given by ω B = |λ| a . Eq. (45) can be solved by introducing a Fourier transform, as made with (31). We arrive, after some algebra, at (valid when α = β). To illustrate the evolution of the WM elements under the effect of the Hamiltonian (43) with a linear potential, we followed this evolution for an initial separable state of the form (13), with ρ S defined by the pure state 1 √ 2 (|0 + |1 ) and W L (m, k) corresponding to a Gaussian state, given by (c.f. Eq. (21)) and N = θ 3 (0, e − 1 σ 2 ) the normalization constant. The results are presented in Fig. 3, which shows different snapshots of the diagonal components W 00 (m, k, t) and W 11 (m, k, t) of the WM. We observe that both components present similar features to the case without spin, plotted in Fig. 2. However, they evolve differently on the m axis: Initially, the W 00 (m, k, t) component moves to the left, while the W 11 (m, k, t) component moves to the right, as a consequence of the different time dependence (−1) α λat in (46), a phenomenon which is reminiscent of the splitting into two beams on the Stern-Gerlach experiment, where the basic piece of the interaction is analogous to (43).

Decoherence
Another dynamical scenario of great relevance for the study of quantum systems is the presence of decoherence, which can be caused by interaction with the environment. In the following we show how the WF formalism we are discussing accommodates also such situation. In particular, we explore some typical cases, in which the spin structure of the WM allows a simple visualization of the decoherence effects.
We consider the case where the interaction with the environment can be described by a Lindblad-type equation [45] where A k are the Lindblad operators, and γ k represent the corresponding coupling constants. If these operators act only on the spin space, the Lindblad (noise) term Eq. (49) immediate translates in an analogous equation for the WM. In other words, under this hypothesis we can write for the WM In the latter equation, ∂ ∂t W (m, k, t)| H denotes the contribution of the Hamiltonian to the dynamics (without decoherence), and we used a matrix notation, so that spin indices are omitted.
As a simple example, let us consider the case when we only have a Lindblad operator A 1 = σ z with γ 1 ≡ γ. We then have which solution can be readily obtained, and expressed as In other words, in this example decoherence leaves the diagonal terms unaltered, while the off-diagonal terms are exponentially damped with time.
Our second example is provided by the Lindblad operator A 1 = σ x with γ 1 ≡ γ. In this case, Eq. (50) becomes This set of equations can be solved by elementary operations. We concentrate on the diagonal terms, for which the final solution reads Similar equations can be obtained involving W 01 (m, k, t) and W 10 (m, k, t). As a result, in the limit t −→ ∞ both W 00 (m, k, t) and W 11 (m, k, t) become an equally weighted mixture (the same happens with the off-diagonal terms).

Quantum walk
The examples studied in the previous Section arise as a consequence of the continuous interaction of a particle with an external potential acting on the lattice. However, we can envisage some situations in which we act on the particle with subsequent short pulses, or via some actions that appear suddenly, but regularly in time. A paradigmatic example of this kind is provided by the quantum walk [39,46], which has received a lot of interest in recent years. In the discrete quantum walk, a quantum particle moves on an (1D) lattice subject to the periodic influence of a displacement operator, that propagates the particle to the right or to the left, according to the state of a two-level system (the coin). The total Hilbert space has precisely the structure H = H l ⊗ H s , defined in Sect. II and, in fact, we can associate the states of the coin to the spin of the particle, without loss of generality. It is customary to use the basis states |L and |R in H s (instead of |0 and |1 ) and associate them to the left and right propagation, respectively. We consider the successive application of the unitary transformation where C(θ) = I ⊗ σ z e −iθσy , θ ∈ [0, π/2] is a parameter defining the bias of the coin toss, I is the identity operator in H l , and σ y and σ z are Pauli matrices acting on H s . The QW dynamics can be described entirely in terms of the WM [47], via a recursion formula that relates W (m, k, t + 1) to other components of this function at time t. Using Eq. (56) one obtains, after some algebra: where M L = (|L L|)C(θ) and M L = (|R R|)C(θ). A complete analysis of the time evolution in phase space with the help of the WF can be found in [47]. Notice that a different definition of the WF was used in [48] for the reduced density matrix of the walker (after tracing the coin) to study of the evolution and the effects of decoherence for the quantum walk.

Decoherence in discrete time
The WF formalism can easily accommodate the description of the general transformation of the quantum state via a completely positive (CP) map. In particular, we consider here trace preserving maps. These could, for instance, represent a decoherent QW process, with Kraus operators modeling the interaction of the system with the environment. The discrete evolution is represented by where E i are Kraus operators with the property i E † i E i = I. As an example, we analyze two simple models of decoherence which are applied as projective measurements in the different degrees of freedom of the system. The first model is defined as projectors in spin space, while the second model is defined by projecting in the lattice sites. We use the notation Π i to designate the different projectors, which satisfy Π † i = Π i and Π i Π j = δ ij Π i . With probability p, the system is projected onto the spin (or space) basis, so that Eq. (58) will be rewritten as By iteration of the above equation and making use of the properties of projectors, one can derive the following formula relating the final and initial density operators of the system, We start from a state consisting of superposition of two deltas with orthogonal spin components, Eq. (15) with α = 1.
For the first projective model we apply the spin projectors Π i = |i i|, i = 0, 1, while for the site projection they are given by Π n = |n n| , n ∈ Z. The iterated density operator ρ(t) that is obtained from Eq. (60) is the same in both cases, the reason being the spin and position entanglement structure in Eq. (15). The result is The corresponding WM becomes Thus, as a consequence of the projective measurements, the non-diagonal components in the WM (62) tend to zero with time. This was expected from the intuitive idea that these components appear from interference between the two spin states in Eq. (15) (or, correspondingly, between the two occupied positions): Once decoherence acts, this kind of interference is reduced and the responsible terms are consequently diminished. Qualitatively similar results are found if one starts from the superposition of two Gaussian states (19), and introduces projective measurements on the lattice states. Interestingly, these interference terms are non positive and tend to disappear as decoherence is acting. We will discuss the consequences of this idea in the next Section with more detail.

V. NEGATIVITY
In the context of continuous variables, it is well known that the Wigner function may present some zones in phase space where it is negative. This is interpreted as an indication of quantumness, in the sense that the state would not have a classical analogue. In order to quantify this quantum feature, the negative volume of the Wigner function has been defined as a measure on non-classicality [6] and has been applied to distinguish quantum states from classical ones [8]. The only pure states with non-negative Wigner function are Gaussian states [49], however the classification is not complete for mixed states.
For the continuous phase space, the negativity of a state ρ becomes The positive character of the Wigner function has also been studied for discrete systems. In the finite dimensional case, and for odd dimension, Gross showed [50] that the only pure states with positive Wigner function are stabilizer states. The presence of negative values in the Wigner function has been in this case connected to a quantum resource, related to a possible quantum speedup [28,51] or the non-simulability of certain quantum computations involving states with non-positive Wigner function [7,31].
In the case of spin 1 2 , the Wigner function defined by Wooters [21] has been used to establish a separability criterion for a system of two particles [52]. A connection between entanglement and negative Wigner functions was established also in [8] for two particles in a continuous space, when the state is a hyperradial s-wave.
Even without the additional degree of freedom, the discreteness of the Hilbert space causes the appearance of spurious negative terms in the Wigner function, which do not correspond directly to non-classical features of the state, but are due to the structure of the discrete phase space itself. Nevertheless, for the case of a spinless particle we showed in [32] that it is possible to introduce a modified negativity measure which excludes such negative contributions and contains information about the quantumness of the states, consistent with the continuum limit.
Not being a true quantum observable, the meaning of such a negativity measure will depend strongly on the definition used for the Wigner function and on the characteristics of the particular system, as the discussion above illustrates. It is then reasonable to ask what natural extension corresponds to the system we are discussing, and what information it maintains about the characteristics of the states.
It is possible to think of several extensions of the Wigner function. If we start with our definition (2) and trace out the spin, we are left with a scalar Wigner function representing the state of the spatial degree of freedom, which in general will be mixed. To this function we can immediately apply the definition of negativity discussed in [32]. It might be more interesting to think of a negativity definition η(ρ) that preserves some spin information.
One possibility is to define a negativity for the Wigner matrix, as in [47], where || A || 1 ≡ T r √ A † A is the trace norm of matrix A, and the second equality follows from normalization. We can easily check that this quantity fulfills the following desirable properties: 1. It reduces to Eq. (63) for product states in the continuum limit, with W c (x, p) obtained from W L (m, k) (see Eq. (13); 2. It is invariant under rotations in spin space.
The first property is also satisfied by the negativity computed after tracing out the spin. The second property, on the other hand, can be illustrated with the following example. We consider an electron, subject to an external magnetic field. To simplify, the electron is confined to a site on the lattice, so that its state is factorizable. The effect of the magnetic field manifests on the precession of the spin, which continuously changes the spin state of the electron. This property ensures that the value of the negativity is not influenced by the precession. In other words, simply changing the spin direction will not alter the negativity properties of the Wigner matrix. Notice that, for some alternative definitions of the Wigner function for a particle with spin [21], the function can contain negative values in the phase space for some states, while being completely positive for other states. We can further explore the significance of the definition (64) by considering different examples. We may then investigate, as in [52], whether this quantity holds information about the entanglement in the state.
We start by analyzing the cat state, | ψ = Another scenario where the emergence of classical behavior is often discussed is that of decoherent dynamics. It is thus reasonable to study how this quantity, η, changes under decoherence. To do so, we consider a very simple situation, in which the initial state subject to decoherence is the double delta considered in Sect. III. To simplify, we restrict ourselves to the discrete time dynamics already studied in Sect. IV, with decoherence arising from projections on spin or lattice sites. Similar qualitative conclusions can be drawn if we allow for a continuous time dynamics, or if we consider a double Gaussian state (19), although calculations are more involved. A simple application of (64) to Eq. (62) leads to the result for the negativity as a function of time. This simple result can be interpreted as the damping of the out-of-diagonal terms in (62). As time goes on, these interference terms tend to fade away, and one is left with an incoherent state with a positive Wigner function. This transition from a coherent superposition to an incoherent one is, of course, a well known phenomenon in the theory of open quantum systems which shows a change in the nature of the Wigner function that is monitored by our definition of the corresponding negativity.
Although it is obvious from this discussion that in the presence of spin the negativity does not have the clear unique physical meaning it had in the purely spatial case (either continuous or discrete), the quantity η introduced here may be useful to characterize some features of the quantum state or the dynamics when the study is restricted to particular families of states. The topic is nevertheless far from being closed, and could be the subject of further debate.

VI. CONCLUSIONS
In this paper, we have elaborated the previously introduced Wigner formalism for a particle in an infinite 1D lattice, in order to account for dynamics and for the presence of an additional, finite-dimensional, degree of freedom. Our goal was to describe the dynamics on the phase space associated to this problem. Although we have concentrated, for simplicity, on the case where such additional degree of freedom corresponds to a spin 1/2, one can envisage more general situations where higher spins, or different properties, such as the polarization of a photon, are considered. As we have showed, the matrix formalism is specially well suited to describe the interaction of the particle with a spin-dependent Hamiltonian on a fixed basis, and keeps a close resemblance to the relativistic WF formalism [9,11], a fact that might be useful in the investigation of the non relativistic limit of a given problem. We have illustrated the construction of the WF by analyzing first some simple static examples, like the "Schrödinger cat" double delta or two-Gaussian states. For these states, the position and spin variables are entangled, and this entanglement manifests in a particular structure of the WM.
We have studied the time evolution of the WM for some simple cases. We have explicitly shown the equation governing the evolution of the WF for a general space-dependent potential. This equation, however, can only be exactly solved for some special cases, as we have done for the case of a linear potential, where one recovers the well known phenomenon of Bloch oscillations. A similar statement is valid for a Hamiltonian that can be factored as a scalar part and a spin operator. We have obtained the equation of motion for a general scalar term, and solved it in the linear case, what allows us to compare with the dynamics in the spinless case. The presence of a "spin dependent force" introduces new features on the dynamics that manifest in phase space. To complete the above description, we have incorporated the role of decoherence which, for some simple examples, can be implemented for the WM in a closed form.
In some physical situations, the interaction appears as short pulses acting on the particle, a paradigmatic example being the Quantum Walk. It is possible to analyze the role of decoherence also in this case, and we have analyzed a simple example for the double delta state, when decoherence appears as projections either on the spatial or in the original spin basis. We have showed that both kind of mechanisms produce the same effect, which translates into a damping of the off-diagonal matrix components.
Finally, we have explored a possible extension of the concept of negativity, as defined for the scalar WF, to the spin 1/2 case. While it is not evident what the physical meaning of such negativity might have once the spin is incorporated to the particle, we have proposed the minimum requirements that, in our opinion, this magnitude should obey, and we have suggested a definition of negativity that fulfills these requirements. Following this proposal, we analyzed how decoherence translates into a decreasing of negativity in the above decoherence model. We also showed that our definition of negativity has not trivial correspondence with entanglement, as clearly indicated by an analysis of the Werner state. We think, however, that it is worth studying further the relationship of the Wigner description to the quantum properties of general states in a lattice.