We propose a local model for general fermionic systems, which we express in the Heisenberg picture. To this end, we shall use a recently proposed formalism, the so-called “Raymond-Robichaud” construction, which allows one to construct an explicitly local model for any dynamical theory that satisfies no-signaling, in terms of equivalence classes of transformations that can be attached to each individual subsystem. By following the rigorous use of the parity superselection rule for fermions, we show how this construction removes the usual difficulties that fermionic systems display in regard to the definition of local states and local transformations.

## I. INTRODUCTION

Since the well-known works of Bell,^{1,2} quantum theory has been labeled as “non-local.” This terminology arises from the fact that quantum systems have been proven, both theoretically and experimentally, to violate Bell's inequalities. As a consequence of this fact, Bell's theorem implies that quantum theory cannot admit a local hidden variable model in terms of a real-valued stochastic theory—hence the term non-local.

However, quantum theory has long been known to be local in a deeper, dynamical sense. First, it satisfies the so-called “no-signaling theorem”; moreover, it satisfies the principle of no-action at a distance, as formulated by Einstein.^{3}

The fact that the principle of no-action at a distance is satisfied is most explicit in the Heisenberg picture formulation. That is the case in both non-relativistic and relativistic quantum theory.^{4,5} Specifically, for a system of *N* qubits, it is possible to define elements of reality attached to each qubit, which fully specify the dynamics of the composite system of *N* qubits jointly. We refer to these elements of reality as “descriptors.”^{4} Such descriptors are not real-valued elements of a stochastic hidden-variable theory, precisely as Bell's theorem requires. They are q-numbers,^{6} i.e., elements of a non-commutative algebra, which generate the algebra of local observables of each qubit. One can recover the dynamics of the system of all qubits by simply considering the set of all local descriptors of each qubit.

This relation between no-signaling and no-action at a distance, which emerges in quantum theory, is not accidental. A recent work^{7} proved formally and in a general way that every no-signaling theory with reversible dynamics admits a so-called local-realistic model, which Raymond-Robichaud^{7} defines as a model where (i) each system has a mathematical object that describes its real factual situation; (ii) given the objects (say, *a* and *b*) fully describing one of two disjoint subsystems *A* and *B*, the objects describing the global system *AB* are entirely determined by the ordered pair (*a*, *b*); (iii) a dynamical transformation that affects system *A* only can only change the object *a* describing *A*, but cannot change the object *b* describing *B*. The objects *a* and *b* are called local elements of reality, or the ontic states of the subsystems *A* and *B*. Furthermore, the authors give a constructive method to define such local elements of reality of any theory that is no-signaling and dynamically reversible—which we shall refer to as the equivalence class formalism. Raymond-Robichaud^{8} uses this formalism to find the local elements of reality in qubit networks; such local elements of reality can be represented in terms of the standard local descriptors of qubits, i.e., the generators of the algebra of observables.^{9}

In this paper, we shall apply these recent results to fermionic quantum field theory. We will show how it admits local elements of reality as q-numbers by explicitly constructing them following the equivalence class formalism. As one would naturally expect, we will demonstrate that the familiar generators of the fermionic algebra can be used to represent the local elements of reality. Our work provides the essential ground to understand better the foundations of fermionic phenomena, such as the fermionic Aharonov–Bohm effect.^{10} Moreover, it helps to ease further the conceptual difficulties with the locality of fermionic quantum field theory and seek a unified description of classical and quantum theories.

The structure of this paper is as follows. We first define fermionic systems (Sec. II) and explain the formal notion of local realism in fermionic field theory (Sec. III). Then, we use the structure of fermionic theory to derive the fermionic descriptors using the equivalence class formalism (Sec. IV). Finally, we discuss the results, commenting on the subtleties of the formalism and suggesting generalizations of our results (Sec. V).

## II. FERMIONIC THEORY

Fermions can be regarded as a type of quantum system. In this work, we shall adopt a general information-theoretic approach to represent them, following previous works.^{11,12} This perspective adheres to the position of studying locality in fermionic systems in terms of modes,^{13} not particles.^{14,15} Within this approach, a fermionic system consists of a set of modes *I* and for each mode $i\u2208I$ there are fermionic creation and annihilation operators $f\u0302i,f\u0302i\u2020$. The ground state $|\Omega \u27e9$ is defined by the requirement that $f\u0302i|\Omega \u27e9=0$. The characteristic fermionic operator algebra is given by the creation and annihilation operators of all modes $i\u2208I$ obeying the following commutation relations:

We focus throughout the paper on systems with a finite set of modes for simplicity. In the Schrödinger picture, the allowed physical states are spanned by the states obtained by acting with the creation operators on the ground state $|\Omega \u27e9$, allowing linear combinations of states with the same parity only. With parity, we mean the parity of the number of fermions in the state. This restriction to allow only superpositions of states with the same parity is the parity superselection rule. Remarkably, this rule must be introduced for fermions to satisfy the no-signaling principle.

As an example, the density operator $\rho =|\Omega \u27e9\u27e81|+|1\u27e9\u27e8\Omega |$ is not allowed to physically exist, otherwise we could have faster than light communication between two parties. Formally, one can state the parity superselection rule as requiring that any physically allowed density operator *ρ* has to commute with the parity operator $\mathbb{P}$. $\mathbb{P}$ has as eigenstates of eigenvalue 1(–1) the states with an even (odd) number of fermions. This parity superselection rule is long known^{16} to restrict the set of physically allowed density operators in a fermionic system. However, it also restricts the set of unitaries and observables that one can apply to the fermionic system. It has been shown recently^{11} that physical unitaries and observables also need to commute with the parity operator $\mathbb{P}$ in order not to be able to signal.

More concretely, in a two-mode system, unitaries and observables with the form $f\u03021+f\u03021\u2020$ are not physically allowed. It follows from the parity superselection rule that if $U\u0302I$ is a parity superselection rule unitary that involves modes in the set *I*, and we consider a mode $j\u2209I$, then the action of $U\u0302I$ on $f\u0302j$ leaves it invariant: $U\u0302I\u2020f\u0302jU\u0302I=f\u0302j$. This result is a direct implication of $U\u0302I$ commuting with the parity operator. This fact is crucial for when we analyze in the Heisenberg picture the local mode structure. We will see that a local unitary in *A* does not affect the elements of reality of a disjoint region *B*. However, before going into detail about descriptors and the Heisenberg picture, let us first set the basis of what we mean by locality. We adhere to the position that locality is a structural requirement for a theory. We believe that the definition provided by Raymond-Robichaud^{7} formalizes perfectly this notion of locality at a structural level.

## III. LOCAL REALISM

After providing a general definition of local-realistic physical theories, Raymond-Robichaud^{7} proves that every physical theory satisfying the no-signaling principle and having reversible dynamics is a local-realistic theory. In this section, we apply this theorem to fermions. We first introduce the notion of a realistic theory, then clarify the notion of locality that we use. Finally, we show how fermionic theory fits into the construction of the equivalence class formalism.

### A. Realism

To define a realistic theory, we use a triad of sets $(P,R,T)$ on which we impose the following structures.

The set $P$ is the set of the elements of the physical theory that describe the phenomenal properties of the theory—i.e., those that are empirically accessible. In other words, an element $\rho \u2208P$ details all the properties that can be observed of a specific configuration of a physical system. For fermionic theories, this is the set of density operators *ρ* that satisfy the parity superselection rule. We call the set $P$ the phenomenal state space and its elements phenomenal states since they give us the observational properties of the system.

The set $T$ is the set of operations that describe the allowed physical transformations. In the fermionic case, this corresponds to the set of parity superselected unitary operators. The set of transformations for fermions forms a group under composition $\xb0$. This fact gives reversible dynamics in the theory. The transformations group $T$ acts on the phenomenal state space $P$, defining a group action. The group action in the fermionic case is given by $\rho \u2032=U\u0302\xb7\rho \xb7U\u0302\u2020$, where $U\u0302\u2208T$ and $\rho ,\rho \u2032\u2208P$.

The set $R$ is the set of states of a physical system describing all properties that entirely specify the configuration of that physical system, rather than merely describing the observational properties of it. We call the set $R$ the ontic state space. Since the term ontic appears more often in the literature^{17} than the terms real or noumenal, we will be using it throughout this work. The group of operations $T$ also acts on $R$, but not necessarily with the same action as in $P$. We denote the action of an element $U\u2208T$ to an element $r\u2208R$ as $r\u2032=U\u2009\u22c6\u2009r$. We will see in Sec. IV how the descriptors are a representation of these ontic states.

Let us consider a familiar example to exemplify the meaning of the ontic states. Consider a qubit composite system with the singlet state $|\psi \u27e9=12(|01\u27e9\u2212|10\u27e9)$, and we apply a local rotation in one of the qubits. We know that the Schrödinger state $|\psi \u27e9$ will remain unchanged. One can imagine in this case that the phenomenal state is $|\psi \u27e9$ since it gives the observational properties of the system. However, we can consider that the fact that we have applied a rotation is relevant for the constitutive properties of the system, even though such rotation does not produce any observational effect. In that case, we could consider that the realized rotation is part of our ontic state.

The phenomenal state must be completely characterized by the ontic state of the system: there must be an epimorphism $\phi $ from $R$ to $P$. In other words, there has to be a mapping $\phi :R\u2192P$ where for every $\rho \u2208P$ exists a $r\u2208R$ such that $\phi (r)=\rho $. Moreover, $\phi $ has to be faithful with respect to the two actions of the group of operations $T$ onto the sets $R$ and $P$. In the fermionic case, for all $U\u0302\u2208T,r\u2208R$, then $\phi (U\u0302\u2009\u22c6\u2009r)=U\u0302\xb7\phi (r)\xb7U\u0302\u2020$.

Thus, if we want to see that the fermionic theory is a realistic theory, we have to find a set $R$, an appropriate group action $\u22c6\u2009$ and an epimorphism $\phi $ that satisfy the conditions above. One could choose $R$ as being $P$, where the action of $T$ is the same as in $P$ and the epimorphism is the trivial identification $\phi (\rho )=\rho $. However, this trivial identification may not follow the conditions to be considered a local-realistic theory. In order to show that fermionic quantum field theory is a local-realistic theory, we need to find a suitable ontic state space for a fermionic theory that satisfies the locality properties described below.

### B. Locality

The notion of locality presented^{7} does not refer to position in any way. It refers to the structure of subsystems.^{18} So, when we talk about locality, we refer to being localized in a subsystem of the global system in consideration.

As we comment in Sec. II, we take the fermionic mode as the primary subsystem of a fermionic system. More concretely, if we have a physical system of *N* modes $N={1,\u2026,N}$, we consider that the systems where the set of modes is a subset of $N$ are the physical subsystems of the theory.

The notion of locality refers to the relation between different fermionic modes of the set $N$. We need to specify which systems our states refer to and how our transformations act in terms of mathematical structure. We will call an object local to a subsystem *A* if it only refers to or acts on that subsystem *A*. We are also interested in restricting objects that live or act in the global system to the particular subsystems.

In the phenomenal state space $P$, this so-called projection map $\pi P$ is given by the partial trace. If we have a phenomenal state of the whole system of *N* modes *ρ*, we say that the object that gives all the phenomenal properties of the system within the subsystem of modes $M={r1,\u2026,rM}\u2286N$ is given by the fermionic partial trace of *ρ* over the set of modes $N\M$ denoted by $\rho r1,\u2026,rM=trN\M(\rho )$. Indeed $\rho r1,\u2026,rM$ is a valid phenomenal state, since it is a parity superselected density operator. Moreover, it can be considered to be local to the subset of modes $M={r1,\u2026,rM}$ since $\rho r1,\u2026,rM$ can be given by combining creation and annihilation operators of the $M$ modes alone.

For operations, it is also necessary to have a notion of an operation being local to a subset of modes $M$. For fermions, an operation that acts on the global system of modes $N,\u2009U\u0302\u2208T$, is local to the subsystem given by the subset of modes $M$ iff $U\u0302$ can be given by a combination of creation and annihilation operators of the modes in $M$ alone.

### C. Fermions in the equivalence class formalism

We now analyze what we are required to obtain in order to find a local-realistic theory for fermions. First, we need a realistic theory; thus, we need an ontic state space with an action of the transformation space and the epimorphism $\phi $ to the fermionic phenomenal state space. In order for a realistic theory to be local for a lattice of systems, we need to be able to define two operations in the ontic state space $R$. First, we need a projection map $\pi SR:R\u2192R$. The projection map sends global ontic states to local ontic states. We require that the ontic projections and the phenomenal projections act in parallel. It is required that for any *ρ* being the phenomenal state of the whole system of fermionic modes $N$, then for the ontic state(s) $R\u2208R$ such that $\phi (R)=\rho $, $\phi (\pi SR(R))=\rho S=\pi SP(\rho )$, where *S* is any subset of the modes of $N$. So, we want a projection to subsystems that is faithful under the $\phi $ epimorphism. These properties give us a notion of ontic states belonging to a subsystem.

Last and crucially, we want a structure where knowing the ontic states of the parts, we know the ontic state of the whole. More specifically, we want to be able to define a map $\u2299:RA\xd7RB\u2192RAB$. Consider a phenomenal state $\rho \u2208P$ of the whole system of modes, such that it has an underlying ontic state $R\u2208R$ given by $\phi (R)=\rho $. Also, consider that then we split the global system of modes in two disjoint pieces by choosing two disjoint subsets of modes *A* and *B*, such that $A\u222aB=N$ and $A\u2229B=\u2205$. Then, we can define ontic states $RA=\pi AR(R),\u2009RB=\pi BR(R)$ by the previously defined projection operation. Now, we require that we can assign $RA\u2299RB=RA\u222aB=R$ uniquely. Observe that this is the condition of separability in the Einstein locality principle.^{19} Moreover, notice that with the usual understanding that density operators are the phenomenal states and the ontic states of the system, such join products cannot be defined. Given an entangled state *ρ* for the disjoint systems *A* and *B*, one can obtain $\pi AP(\rho )=trB(\rho )=\rho A$ and $\pi BP(\rho )=trA(\rho )=\rho B$, where the projection map is the partial trace. However, from the un-entangled product state $\rho \u2032=\rho A\u2227\rho B$, using the same map we also get to $\pi AP(\rho \u2032)=trB(\rho \u2032)=\rho \u2032A=\rho A$ and $\pi BP(\rho \u2032)=trA(\rho \u2032)=\rho \u2032B=\rho B$. Thus, we cannot define a join product $\rho A\u2299\rho B=\rho $ and $\rho A\u2299\rho B=\rho \u2032$ since $\rho \u2260\rho \u2032$. Thus, the usual realistic view of fermionic theory is not local-realistic. However, this does not stop us from discussing if it is possible to find a different ontic state space where the local realistic structure manifests. We see in Sec. IV how using the constructive proof of the equivalence class formalism theorem, we find such a structure.

In the equivalence class formalism, an operational no-signaling theory needs a phenomenal state space $P$, a set of transformations $T$ with an action $\xb7$ on $P$. Moreover, it needs a lattice of systems with associated faithful projectors $\pi SP$ and that the no-signaling condition is respected. We can understand the condition of $T$ having a group structure and giving a group action on $P$ as that the dynamics of the theory are reversible. This reversibility is the case in fermionic theory since the set of transformations $T$ corresponds to the set of fermionic parity superselected unitaries for a given complete set of modes. Such unitaries set, of course, conforms a group with the usual product. The phenomenal state space $P$ for fermionic theory is, as mentioned before, the set of parity superselected density operators for the given total set of modes. Given $U\u0302\u2208T,\u2009\rho \u2208P$ the group action is defined as $U\u0302\xb7\rho \xb7U\u0302\u2020$. Furthermore, the partial trace over modes in the complete set of modes gives the projections of the system lattice in the fermionic case. If we have a total set of modes $N$, we can define any subsystem by any subset $M\u2286N$ and we obtain $\pi MP(\u2212)=trN\M(\u2212)$.

Finally, we express the no-signaling principle presented in the equivalence class formalism as: given a system with a bipartition of disjoint subsystems *A* and *B*, any operation local in *A* $UA\u2208T$, any operation local in *B* $VB\u2208T$, and any phenomenal state of the global system *AB* $\rho \u2208P$, the following equation is satisfied:

which in fermionic theory corresponds to the usual no-signaling principle equation given by

Therefore, since the above no-signaling condition is satisfied for the fermionic system that considers the parity superselection rule for both states and unitaries,^{11} all the conditions of the equivalence class formalism are satisfied. Therefore, we can find with their construction the ontic states that dote fermionic theory of the structure of local realism.

The proof that any operational no-signaling theory is a local-realistic theory is constructive.^{7} The authors construct for a general operational theory with a group of transformations $T$, a lattice of systems, a phenomenal state space $P$, with the associated group action $\xb7$ and projection operators $\pi AP$ the ontic state space $R$, the ontic group action $\u22c6$, the epimorphism $\phi $, the ontic projection operators $\pi AR$, and the join product $\u2299$. They do so by defining an equivalence relation $\u223cA$ that states that two transformations *U*, *V* are equivalent in a subsystem *A* $U\u223cAV$ if and only if there exists a transformation $WA\xaf$ local on the subsystem $A\xaf$, complementary to *A*, such that $U=WA\xafV$. The equivalence classes $[U]A$ under this equivalence are the ontic states of the subsystem *A*, forming an ontic state space $R$ with all the desired properties. Some of them are

where *A* and *B* are systems, with $A\u2229B=\u2205$ and $U,V,WA$ are transformations with *W _{A}* being local in

*A*.

In order to find these equivalence classes in the fermionic setting, we utilize the equivalence between the descriptors and the equivalence class formalism in qubit networks.^{20} The paper aims to generalize it to the fermionic case and use the excellent properties and physical intuition that the descriptor picture provides to unveil the set of ontic states that give the fermionic theory a local realistic structure.

## IV. DESCRIPTORS

As we have explored in Sec. III, in order to provide a fermionic theory with a local realistic structure for an orthonormal set of modes, we need to find a set of ontic states $R$ with proper operations $\pi AR,\u2009\phi $ and $\u2299$. In order to construct these objects presented in the equivalence class formalism, we utilize the descriptor formalism presented by Deutsch and Hayden^{4} and then we use the result by Bédard^{20} that establishes the equivalence between the two formalisms.

Their work focuses on qubit networks, but we will see that one can apply their ideas to fermionic systems. Here, we will be applying their procedure directly to the fermionic case. Using the Heisenberg picture, the descriptor formalism consists of re-interpreting the real elements of quantum mechanics. The formalism considers pure states with unitary evolutions, assuming that one can always increase the dimension of the Hilbert space and use the Stinespring dilation theorem^{21} to treat any quantum channel as a unitary evolution. Thus, from here on, we consider the Heisenberg state being pure unless stated otherwise.

The derivation of descriptors starts with the fact that in quantum theory, the expectation value of an observable $O\u0302$ of a system in a particular *ρ* is given by $tr(O\u0302\xb7\rho )$. Now, if we apply a unitary transformation $U\u0302$ to *ρ*, the updated data will be given by $tr(O\u0302\xb7U\u0302\xb7\rho \xb7U\u0302\u2020)$. We use the Heisenberg picture. We consider that the state *ρ* remains invariant under the transformation $U\u0302$, but the observables $O\u0302$ evolve under the transformation given by the unitary operator $U\u0302$ as $O\u0302\u2032=U\u0302\u2020\xb7O\u0302\xb7U\u0302$. Since the observables of the system have evolved and *ρ* has not, it is reasonable to assume that we could find a framework where the observables are the ones that determine the state of the system and not *ρ*.

In order to do so, it is helpful to try to find which operators we can track to determine any $U\u0302\u2020\xb7O\u0302\xb7U\u0302$. We start by defining that the space that we work with consists of a fermionic space of *N* modes under the parity superselection rule. Therefore, the observables $O\u0302$ are Hermitian operators that commute with the parity operator. Using the anticommutation relations of the creation and annihilation operators, we can write any observable $O\u0302$ as

where $O\u0302i(j)$ are monomials of any degree of the creation and annihilation operators $f\u0302i,f\u0302i\u2020$. Notice that for each *j* in the sum, the global operator $O\u03021(j)\xb7\u2026\xb7O\u0302N(j)$ is a monomial of even degree of $(f\u0302i,f\u0302i\u2020)i=1N$. Now, we can see that for any $O\u0302$, we have this decomposition and, then,

Therefore, if we know how all the possible local operators $O\u0302i(j)$ transform under $U\u0302$, then we know how any observable $O\u0302$ evolves under $U\u0302$. However, we can reduce this further. Due to the fermionic algebra, it is easy to see that there are only four linearly independent monomials of $fi,fi\u2020$. Thus, the transformation under $U\u0302$ of $O\u0302i(j)$ will be given by

If we draw our attention to the second line, we can observe that $U\u0302\u2020\xb7(f\u0302i\xb7f\u0302i\u2020)\xb7U\u0302=(U\u0302\u2020\xb7f\u0302i\xb7U\u0302)\xb7(U\u0302\u2020\xb7f\u0302i\u2020\xb7U\u0302)$ and that $U\u0302\u2020\xb7(f\u0302i\u2020\xb7f\u0302i)\xb7U\u0302=(U\u0302\u2020\xb7f\u0302i\u2020\xb7U\u0302)\xb7(U\u0302\u2020\xb7f\u0302i\xb7U\u0302)$. Therefore, we can conclude that if we know the transformation under $U\u0302$ of the 2*N* operators $f\u0302i$ and $f\u0302i\u2020$ for $i=1,\u2026,N$, then we know the evolution under $U\u0302$ of any observable $O\u0302$. However, we can even reduce this further, since it can be observed that $(U\u0302\u2020\xb7f\u0302i\xb7U\u0302)\u2020=U\u0302\u2020\xb7f\u0302i\u2020\xb7U\u0302$.

This fact implies that we only need to know the evolution under $U\u0302$ of the *N* operators $f\u0302i$ in order to know the evolution under $U\u0302$ of any observable of the *N*-mode fermionic system. We call the set $(f\u03021,\u2026,f\u0302N)$ the set of fermionic descriptors of an *N*-mode fermionic system.

By knowing the initial state of the system $|\psi \u27e9\u27e8\psi |$ and knowing the expression of $(U\u0302\u2020\xb7f\u03021\xb7U\u0302,\u2026,U\u0302\u2020\xb7f\u0302N\xb7U\u0302)$, one can recover all the expectation values of any observable in the current state of the system. In the Schrodinger picture, such a phenomenal state would be given by $U\u0302|\psi \u27e9$. A few things to note is that each $U\u0302\u2020\xb7f\u0302i\xb7U\u0302$ can be a global *N*-mode operator, so it can have contributions in its expression from other modes $j\u2260i$. However, at the same time, the “updated” descriptors satisfy the fermionic algebra. In the sense that if we name $f\u0302\xafi=U\u0302\u2020\xb7f\u0302i\xb7U\u0302$, then ${f\u0302\xafi,f\u0302\xafj}=0$ and ${f\u0302\xafi,f\u0302\xafj\u2020}=\delta ijI\u0302$.

### A. Descriptors as ontic states

The curious reader may be wondering how these objects relate to the notion of local realism that we have introduced in Sec. III. The answer is inspired by the work of Bédard.^{20} The set of descriptors $U\u0302\u2020\xb7(f\u03021,\u2026,f\u0302N)\xb7U\u0302$ is a compact way to represent the equivalence classes $[U\u0302]i$ that according to the equivalence class formalism can be considered as the ontic states of the system. This process endows fermionic theory with a local-realistic structure.

We develop the specifics of this claim in the following lines, where we present the results that justify the connection between the fermionic descriptors that we have proposed and the equivalence class formalism of local-realistic theories.

**Theorem 1.** *Given a fermionic theory with a set of modes* $I={i1,\u2026,iN}$*. Given the equivalence relation on the group of parity superselected unitaries of the theory for each mode* $ij\u2208I$ *given by* $U\u0302\u223cijV\u0302$ *iff* $U\u0302=W\u0302I\{ij}\xb7V\u0302$*, then*

Thus, $[U\u0302]ij={V\u0302\u2208T|U\u0302\u2020\xb7f\u0302ij\xb7U\u0302=V\u0302\u2020\xb7f\u0302ij\xb7V\u0302}$.

With Theorem 1, we have obtained a direct connection between the fermionic descriptors and the ontic states of the equivalence class formalism. The proof follows easily from the algebraic properties of unitaries, and it is in the Appendix. We can define the equivalence classes that are the local elements of reality in terms of properties satisfied by the descriptors. Moreover, we see that we can define the equivalence class that gives the local ontic state in a subsystem in terms of only the descriptors associated with that subsystem.

#### 1. Ontic action

Using the result from Theorem 1, we can represent the ontic fermionic states of the equivalence class formalism using the fermionic descriptors. The action $\u2009\u22c6\u2009$ of the group of transformations $T$ on the ontic state space $R,\u2009W\u2009\u22c6\u2009[U]S$ that we have defined in Sec. III in the representation by descriptors is given by

In what follows, we express the join product $\u2299$, the ontic projectors $\pi AR$, and the epimorphism $\phi $ in respect of fermionic descriptors. Consider an *N*-mode fermionic system with set of modes $I={i1,\u2026,iN}$. We have an initial phenomenal Heisenberg state $\rho 0=|\psi 0\u27e9\u27e8\psi 0|$. The system undergoes a specific transformation $U\u0302$ given as a polynomial of the operators $f\u0302ij,f\u0302ik\u2020$. After such a transformation is applied, the pure state $|\psi \u27e9=U\u0302|\psi 0\u27e9$ is the phenomenal state of the system. In regard to the descriptor picture, we have that in respect of the lattice of systems given by the canonical mode set *I*, the set of descriptors $(f\u0302i1,\u2026,f\u0302iN)$ and the initial Heisenberg state $|\psi 0\u27e9$ represent the initial ontic state of the system. After the transformation $U\u0302$ is applied, the ontic state of the system can be represented by $(U\u0302\u2020f\u0302i1U\u0302,\u2026,U\u0302\u2020f\u0302iNU\u0302)$ and the Heisenberg state $|\psi 0\u27e9$.

#### 2. Descriptor epimorphism

A result that is important to understand the completeness of the descriptor picture and its relationship with Theorem 1, the group of transformations $T$, and its equivalence classes is the following Theorem 2. The following theorem helps to understand the definition of the epimorphism $\phi $ for a complete set of descriptors.

**Theorem 2.** *Using the complete set of descriptors* $(f\xaf\u03021,\u2026,f\xaf\u0302N)$, it *is possible to find uniquely the* $U\u0302$ *that acts on the N-mode system such that* $(f\xaf\u03021,\u2026,f\xaf\u0302N)=(U\u0302\u2020f\u03021U\u0302,\u2026,U\u0302\u2020f\u0302NU\u0302)$.

The proof of this result is entirely algebraic, but the idea is to use a decomposition of the unitary operators on a basis of the operators that act as descriptors. The complete proof is in the Appendix.

The construction of the epimorphism $\phi $ seems straightforward from this theorem. We choose the annihilation operators to describe the ontic state of our initial configuration, with a given Heisenberg state $|\psi 0\u27e9\u27e8\psi 0|$. We say that for such a configuration, the associated phenomenal state is the Heisenberg state. Then, if we have a different set of descriptors that is unitarily conjugated to the canonical one, one can obtain the unitary $U\u0302$ that connects them, and associate to the new configuration the phenomenal state $\phi ((U\u0302\u2020f\u0302i1U\u0302,\u2026,U\u0302\u2020f\u0302iNU\u0302),|\psi 0\u27e9\u27e8\psi 0|)=U\u0302|\psi 0\u27e9\u27e8\psi 0|U\u0302\u2020$.

This definition is consistent with what one might expect since we start with a phenomenal state $|\psi 0\u27e9\u27e8\psi 0|$ that we use as a Heisenberg state. The mentioned Heisenberg state and a corresponding set of canonical annihilation operators give the initial ontic state. Then, a unitary evolution makes the states evolve. The action of the unitary updates the descriptors. Moreover, since we can deduce the unitary, the new phenomenal state is the unitary acting on the initial phenomenal state. So, the definition of the epimorphism is consistent with our physical intuition on how transformations should act on the ontic and phenomenal states.

However, this only defines the epimorphism for a total set of descriptors, so for the global states of the system. We now present how from the ontic local state $((U\u0302\u2020f\u0302i1U\u0302\u2026,U\u0302\u2020f\u0302iMU\u0302),|\psi 0\u27e9\u27e8\psi 0|)$, we associate the phenomenal state *ρ _{M}*, which is a parity superselected density operator for the system of

*M*modes $i1,\u2026,iM$. The idea is the same as for the proof of Theorem 2. It is not difficult to see that

*ρ*can be decomposed in an operator basis $O\u0302M(j)$ given by monomials of creation and annihilation operators $f\u0302i1,f\u0302i1\u2020,\u2026,f\u0302iM,f\u0302iM\u2020$, then manipulating the expression we get

_{M}Since $U\u0302\u2020O\u0302M(j)U\u0302$ can be calculated from the descriptors $((U\u0302\u2020f\u0302i1U\u0302\u2026,U\u0302\u2020f\u0302iMU\u0302),|\psi 0\u27e9\u27e8\psi 0|)$ and $|\psi 0\u27e9\u27e8\psi 0|$ is part of the ontic state, we can assign from our ontic local states the phenomenal local states through the epimorphism $\phi $ given by

The complete details of the decomposition we use are in the proof of Theorem 2 in the Appendix.

### B. Separability

Once we have a representation of the epimorphism $\phi $ in terms of fermionic descriptors, we can discuss separability. Separability is the condition that gives the notion of locality to the ontic state-space construction. In the fermionic descriptor setting, we now discuss how we represent the ontic local states, how the join product operates, and how they relate to the local phenomenal states (the reduced density operators).

#### 1. Ontic projections

We know that given a subset *I _{A}* of the mode set $IA\u2286I$: ${a1,\u2026,aM}\u2286{i1,\u2026,iN}$ induces an

*M*mode subsystem

*A*of the

*N*mode system

*S*, where

*A*is the fermionic system of

*M*modes generated by the mode set $IA={i1,\u2026,iM}$. At the ontic state space in the equivalence class formalism to define the ontic projection mappings of the equivalence classes, one just defines $\pi AR([U]AB)=[U]A$. It is straightforward to see that the following simple operation defines the ontic projection mappings in the descriptor representation of fermionic systems by using Theorem 1.

**Definition 3.** *If we have an ontic global state of the system S represented by* $((U\u0302\u2020f\u0302i1U\u0302,\u2026,U\u0302\u2020f\u0302iNU\u0302),|\psi 0\u27e9\u27e8\psi 0|)$*. Then, we define*

#### 2. Joint ontic product

To define the ontic join product, we need to see that for any fermionic subsystem *A* defined in terms of modes subsets, we can consider a bipartition of the global system *S*. We can consider the fermionic system *B* of *N*–*M* modes spanned by the set of modes $IB=I/IA$. It is straightforward to see that *A* and *B* are disjoint systems and that the global system *S* can be considered the join space of both subsystems.

Now, we are in a position to define ontic states of systems *A* and *B* to be compatible. Consider an ontic state of system *A* represented with descriptors $((U\u0302\u2020f\u0302a1U\u0302,\u2026,U\u0302\u2020f\u0302aMU\u0302),|\psi \u27e9\u27e8\psi |)$ and an ontic state of system *B* represented with descriptors $((V\u0302\u2020f\u0302b1V\u0302,\u2026,V\u0302\u2020f\u0302bN\u2212MV\u0302),|\eta \u27e9\u27e8\eta |)$. Then, the state is said to be compatible if there exists a global state of the system $S=A\u222aB$, represented by descriptors as $((W\u0302\u2020f\u0302i1W\u0302,\u2026,W\u0302\u2020f\u0302iNW\u0302),|\varphi \u27e9\u27e8\varphi |)$ such that when projected to subsystems *A* and *B* equals the beforementioned states.

Note that this means that for two states to be compatible, they have to have the same Heisenberg state and that there must exist a unitary transformation *W* such that $[W]A=[U]A$ and $[W]B=[V]B$. Then, it is straightforward for compatible states to see how the join product can be defined and operates in the descriptor representation.

**Definition 4.** *Consider two compatible local ontic states of disjoint subsystems A and B,* $((U\u0302\u2020f\u0302a1U\u0302,\u2026,U\u0302\u2020f\u0302aMU\u0302),|\psi 0\u27e9\u27e8\psi 0|)$ *and* $((V\u0302\u2020f\u0302b1V\u0302,\u2026,V\u0302\u2020f\u0302bN\u2212MV\u0302),|\psi 0\u27e9\u27e8\psi 0|)$*. Then, there exists a unique global state* $((W\u0302\u2020f\u0302i1W\u0302,\u2026,W\u0302\u2020f\u0302iNW\u0302),|\psi \u27e9\u27e8\psi |)$ *such that the definition of the join product of the compatible states can be defined as*

Thus, due to the uniqueness property of the definition, knowing the ontic states of the subsystems is the same as knowing the ontic state of the global system. The projection and join structures are the same for the fermionic case to the qubit case.^{4,7} However, the critical result for the fermionic case is that for a global system where we know the unitaries that are applied, the projection and joining operations consist in separate and merge the descriptors of the relevant set of modes that conform to our subsystems of interest.

### C. Faithfulness of splitting operation

Of course, we can repeat the bipartition process until we reach the point of individual modes forming subsystems. We see how the local and global states are on the same footing, which gives us strong notions to believe in their physicality given a particular partition of the global system in terms of subsystems. A particular excellent result that one can obtain from the representation of ontic states with descriptors is that splitting into subsystems is a faithful operation because the diagram in Fig. 1 commutes. Theorem 5 condenses the information in Fig. 1.

**Theorem 5.** *Given a fermionic system with the mode set* $I={1,\u2026,N}$*. For any set of descriptors* $(f\xaf\u03021,\u2026,f\xaf\u0302N)$ *together with any Heisenberg state* $|\psi 0\u27e9\u27e8\psi 0|$ *and any non-empty subset of modes* $J\u2286I$*. The associated diagram of Fig. 1 commutes. In other words,*

The proof of the theorem is algebraic. We apply the form of the epimorphism $\phi $ for descriptors that we have on Eq. (11). We also use the definition of ontic projection $\pi AR$ that consists in forgetting the descriptors for the modes that are not in *A*. Furthermore, finally, we use the properties of $\pi AP$ being the fermionic partial trace. Moreover, we use some algebraic properties of fermionic systems under the parity superselection rule derived in previous works.^{11} The complete proof is in Appendix.

The definitions and theorems allow us to claim that fermions also have as qubits a local-realistic structure. We have found explicitly a neat representation of the ontic states using the widely known Heisenberg picture of quantum mechanics. Furthermore, we have seen that tracking the creation or annihilation operators is enough to have a complete description of the ontic states of the system. In Sec. V, we discuss the implications of our findings and the possible uses that this way of representing fermionic systems may have.

## V. DISCUSSION

Following the equivalence class formalism, this work shows how one can formally prove that fermionic systems are local-realistic. We have seen how a compact representation of the local structure can be given by extending the descriptor picture of qubit networks to fermionic modes. We have also seen that fermionic annihilation operators form a valid set of fermionic descriptors.

Readers familiar with the Heisenberg picture in quantum field theory may wonder how the analysis we have performed differs from it. We are indeed just using the usual Heisenberg picture. We have seen that the most efficient way to use the Heisenberg picture is to track the fermionic annihilation operators since, from their evolution, we can retrieve the evolution of any observable of the theory. However, the interpretation of these objects by the equivalence class formalism differentiates the usual Heisenberg picture from our treatment. First, we have stressed the relevance and necessity of having local objects associated with each fermionic mode, which are the subsystems in our theory and second, we propose to interpret the dynamical evolution of the annihilation operators as a representation of the equivalence classes of parity superselected unitaries that correspond to the ontic states in our system. This consistent interpretation is a novel element that introduces a possible new angle on the meaning of states and dynamics in fermionic systems.

Marletto *et al.*^{22} have analyzed how a phase in a fermionic Mach–Zender interferometer is acquired locally by using the Heisenberg picture and tracking some relevant observables. This paper is a generalization of that work, using the new construction of the equivalence class formalism. We specify the notion of a local element of reality in terms of local ontic states, and we can then describe any fermionic physical process with a local description.

We emphasize, as also pointed out by Chiribella,^{18} that the lattice of systems of a theory is an extra layer of the structure that the fermionic algebra alone does not determine. More concretely, in the fermionic case, the lattice of systems is determined by which notion of modes one chooses. Moreover, under that partition, we can assign the ontic local states in terms of subsystems of the theory. However, the choice of modes is not unique. We could perfectly define a different set of modes by unitary conjugation to the chosen one that would define a different partitioning of the global system in terms of subsystems. Furthermore, in this other partition, we would be able to identify the ontic local states. Therefore, even though we can uniquely identify the local elements of reality given a notion of locality, the locality is not unique for a group of physical transformations.

For a set of local operators to be a suitable set of descriptors, we need that the evolution of such operators suffices to determine the evolution of any observable. One can wonder if there is a straightforward way to understand why the annihilation operators are suitable descriptors. There is. We first observe that fermionic observables are fermionic operators, which form an algebra. We now notice that in this algebra, in particular, the 2*N* creation and annihilation operators are generators of the algebra of fermionic operators and since $U\u0302\u2020(\alpha A\u0302+B\u0302)U\u0302=\alpha U\u2020(A)U+U\u2020(B)U$ and $U\u0302\u2020(A\u0302B\u0302)U\u0302=U\u0302\u2020(A\u0302)U\u0302U\u0302\u2020(B\u0302)U\u0302$, tracking the generators of the algebra is enough to track any element of it. Using the fact that $(f\u0302)\u2020=f\u0302\u2020$, we can claim that tracking $f\u0302$ is enough to track $f\u0302\u2020$, so it is enough to track the annihilation operators.

This observation is useful in understanding why in the *N*-qubit network case, the descriptors can be the $\sigma \u0302jx,\sigma \u0302jz$ operators for each qubit *j*. These $\sigma jx,\sigma jz$ are the generators of the *N*-qubit network operator algebra. Following the same spirit as in this article, we can now also use the $\u2020$ operation properties. It is straightforward to see that if one defines $q\u0302j=12(\sigma \u0302jx+i\sigma \u0302jy)$, then $q\u0302j,q\u0302j\u2020$ are generators of the *N*-qubit algebra (see in Appendix). We can see that since $(q\u0302j)\u2020=q\u0302j\u2020$, it is enough to track all the $q\u0302j$ to track any *N*-qubit observable. Therefore, we only need to track one operator for each qubit, with the caveat of that operator not being an observable since it is a non-Hermitian operator.

These operators $q\u0302j$ behave like annihilation operators. They satisfy the condition that $q\u0302j|0\u27e9\u2026|0\u27e9=0,\u2009q\u0302j2=0,\u2009{q\u0302j,q\u0302j\u2020}=I\u0302$ and $[q\u0302j,q\u0302i]=0$. At the individual qubit/mode level, the anticommutation relations are the same as for fermionic modes, but operators of different modes/qubits commute. This realization is not related to the Jordan–Wigner transformation; it is just a property that arises from the defining commutation/anticommutation relations of generators of the qubit network algebra.

We can now postulate what the descriptors will be for any quantum system of particles. For the bosonic case, we can see that the direct candidate for bosonic descriptors in a system of *N* bosonic modes is the set of bosonic annihilation operators $b\u0302j$. We know that the set of $b\u0302j,b\u0302j\u2020$ is the generators of the bosonic operator algebra. We can simplify further, by tracking only the set of $b\u0302j$—this is enough to track any bosonic observable. In our future research, we would like to use the notion of algebra generators to study the locality and descriptors of anyonic systems,^{23} where the identifications are not so clear due to the lack of field theory creation and annihilation operators. Similarly, we would like to expand on these concepts focusing on well-known quantum algebras such as the Virasoro algebra used in string theory to find a compact local description of such systems.

A confusing question can arise in our main fermionic case. Due to the parity superselection rule, $f\u0302i$ are not fermionic observables, nor fermionic unitaries, nor a linear combination of them; how can the collection of them together with $|\psi 0\u27e9\u27e8\psi 0|$ be the ontic state of the parity super selected fermionic system? We answer this by regarding the set of fermionic annihilation operators and $|\psi 0\u27e9\u27e8\psi 0|$ not as the ontic state of the fermionic system, but as a representation of the ontic state. We use here representation as a broad term, not the precise meaning used in group theory. The ontic states are the equivalence classes of the parity superselected unitaries that we have described with the equivalence relation. Note that only the parity superselected unitaries (real, physically allowed) are required to define such equivalence classes as the ontic states. We use the descriptors and the Heisenberg picture to have a compact and convenient way of representing the ontic states of the parity superselected unitaries' equivalence classes and their structure.

Finally, in the fermionic and the general equivalence class formalism case, only the set of physically allowed transformations is necessary to define the ontic states of a physical system. Thus, we think these results suggest that we should study the transformation structures by themselves to see what range of structures and phenomena can offer. Identifying which transformations are physically possible seems crucial in understanding the constitutive properties and the locality features of a physical system. Further future work could involve extending the equivalence class formalism to nonreversible dynamics so that theories such as the objective collapse model^{24} could be studied in detail.

## ACKNOWLEDGMENTS

The authors thank David Deutsch, Paul Raymond-Robichaud, and Charles Alexandre Bédard for fruitful discussions. N.T.V. acknowledges financial support from “la Caixa” Foundation (ID 100010434, LCF/BQ/EU18/11650048). C.M. and V.V. thank the John Templeton Foundation and the Eutopia Foundation. V.V.'s research is supported by the National Research Foundation and the Ministry of Education in Singapore and administered by the Centre for Quantum Technologies, National University of Singapore. This publication was made possible through the support of the ID 61466 grant from the John Templeton Foundation, as part of the Quantum Information Structure of Spacetime (QISS) Project (qiss.fr). The opinions expressed in this publication are those of the authors and do not necessarily reflect the views of the John Templeton Foundation.

## AUTHOR DECLARATIONS

### Conflict of Interest

The authors have no conflicts to disclose.

## DATA AVAILABILITY

Data sharing is not applicable to this article as no new data were created or analyzed in this study.

### APPENDIX: PROOFS AND MATHEMATICAL DETAILS

In this section, we include mathematical details of the paper and laborious proofs of statements made in the main article.

(**Theorem 1**). *Given a fermionic theory with a set of modes* $I={i1,\u2026,iN}$*. Given the equivalence relation on the group of parity superselected unitaries of the theory for each mode* $ij\u2208I$ *given by* $U\u0302\u223cijV\u0302$ *iff* $U\u0302=W\u0302I\{ij}\xb7V\u0302$*, then*

*Thus,*$[U\u0302]ij={V\u0302\u2208T|U\u0302\u2020\xb7f\u0302ij\xb7U\u0302=V\u0302\u2020\xb7f\u0302ij\xb7V\u0302}$.

*Proof.* The last statement follows directly from the definition of an equivalence class, so the equation that needs to be proven is Eq. (8):

“$\u21d2$”: $U\u0302\u223cIjV\u0302$ implies $U\u0302=W\u0302I\{ij}\xb7V\u0302$ for some $W\u0302I\{ij}$ being a parity superselected unitary local on the set of modes that excludes the mode *i _{j}*. Thus, since $W\u0302I\{ij}$ is an even operator that does not contain any terms involving $f\u0302ij,f\u0302ij\u2020$, it is straightforward to check that $[W\u0302I\{ij},f\u0302ij]=0$. Thus, $U\u0302\u2020\xb7f\u0302ij\xb7U\u0302=V\u0302\u2020\xb7W\u0302I\{ij}\u2020\xb7f\u0302ij\xb7W\u0302I\{ij}\xb7V\u0302=V\u0302\u2020\xb7f\u0302ij\xb7W\u0302I\{ij}\u2020\xb7W\u0302I\{ij}\xb7V\u0302=V\u0302\u2020\xb7f\u0302ij\xb7V\u0302$.

“$\u21d0$”: We have that $U\u0302\u2020\xb7f\u0302ij\xb7U\u0302=V\u0302\u2020\xb7f\u0302ij\xb7V\u0302$, and to see that $U\u0302\u223cijV\u0302$ we need to see that $U\u0302=W\u0302I\{ij}\xb7V\u0302$. Or, equivalently, since we have a group structure, where transformations are unitaries, proving that $U\u0302\xb7V\u0302\u2020=W\u0302I\{ij}$ is enough to prove that $U\u0302\u223cijV\u0302$. From $U\u0302\u2020\xb7f\u0302ij\xb7U\u0302=V\u0302\u2020\xb7f\u0302ij\xb7V\u0302$, it is straightforward to deduce that $f\u0302ij\xb7(U\u0302\xb7V\u0302\u2020)=(U\u0302\xb7V\u0302\u2020)\xb7f\u0302ij$. Naming $U\u0302\xb7V\u0302\u2020=W\u0302$, noticing $W\u0302$ is a parity superselected unitary and taking the dagger of the previous equation, we have that the two following equalities hold:

Moreover, now since *W* is *a priori* a general parity superselected unitary, is not difficult to see that we can decompose it as follows: $W\u0302=O\u03020+f\u0302ijO\u03021+f\u0302ij\u2020O\u03022+f\u0302ijf\u0302ij\u2020O\u03023$, where $O\u03020,O\u03023,O\u03021,O\u03022$ are local operators on the set of modes $I\{ij}$, where *O*_{0}, *O*_{3} are even operators and $O\u03021,O\u03022$ are odd operators (since the operator $W\u0302$ is an even operator). Using this decomposition of $W\u0302$ in the first condition of Eq. (A1) and commuting/anticommuting the $f\u0302ij,f\u0302ij\u2020$ terms with the $O\u0302k$ operators, we obtain that

Is not difficult to see that this implies that $O\u03022=0\u0302$ and $O\u03023=0\u0302$. Then, using that $W\u0302=O\u03020+f\u0302ijO\u03021$ and replacing in the second condition of Eq. (A1), we obtain

It is also easy to see that it implies that $O\u03021=0\u0302$. Thus, we have seen that the conditions imply that $W\u0302=O\u03020$, thus being a local parity superselected unitary on the set of modes $I\{ij}$, thus allowing us to name $W\u0302=W\u0302I\{ij}$, and therefore we have proven that $U\u0302\u223cijV\u0302$. ◻

(**Theorem 2**). *Using the complete set of descriptors* $(f\xaf\u03021,\u2026,f\xaf\u0302N)$, it *is possible to find uniquely the* $U\u0302$ *that acts on the N-mode system such that* $(f\xaf\u03021,\u2026,f\xaf\u0302N)=(U\u0302\u2020f\u03021U\u0302,\u2026,U\u0302\u2020f\u0302NU\u0302)$.

*Proof*. First, note that the unitary conjugation action is up to a global phase. Keeping this in mind, considering the unitary $U\u0302$ as an operator, we can regard it as a vector on the vector space of operators. We can consider in the operator algebra an orthonormal basis, where the scalar product between two operators $A\u0302,B\u0302$ is given by $tr(A\u0302\u2020B\u0302)$. It is straightforward to see that for the fermionic operators, if one labels the orthonormal Fock basis as ${|k\u27e9}k=12N$, then $|k\u27e9\u27e8l|$ is an orthonormal basis of the operator vector space with scalar product given by $tr{A\u0302\u2020B\u0302}$.

Considering this, now $U\u0302$ can be written as $U\u0302=\u2211k,l=12Ntr(U\u0302|k\u27e9\u27e8l|)|k\u27e9\u27e8l|$. So, if we know $tr(U\u0302|k\u27e9\u27e8l|)$ we know the unitary. $|k\u27e9\u27e8l|$ is a product of creation and annihilation operators, since $|k\u27e9=f\u0302i1\u2020\u2026f\u0302in\u2020|\Omega \u27e9$ and $|\Omega \u27e9\u27e8\Omega |=f\u0302N\u2026f\u03021f\u03021\u2020\u2026f\u0302N\u2020$. Using this fact, having at our disposal $(f\u03021,\u2026,f\u0302N)$ and $(f\xaf\u03021,\u2026,f\xaf\u0302N)$, we can construct $|k\u27e9\u27e8l|$ and $|k\xaf\u27e9\u27e8l\xaf|$, where to construct the second we have replaced the $f\u0302i,f\u0302i\u2020$ in the decomposition of $|k\u27e9\u27e8l|$ by $f\xaf\u0302i,f\xaf\u0302i\u2020$. We can see easily that $|k\xaf\u27e9\u27e8l\xaf|=U\u0302\u2020|k\u27e9\u27e8l|U\u0302$.

We take the scalar product of any two of these objects. In other words, consider $tr(|k\xaf\u27e9\u27e8l\xaf||m\u27e9\u27e8n|)=tr(U\u0302\u2020|k\u27e9\u27e8l|U\u0302|m\u27e9\u27e8n|)$. Now, if we insert the decomposition of $U\u0302$ found above and we use the linearity properties of the trace, we obtain that

where we have used the orthonormality of the Fock basis and the properties of the Kronecker delta. Using the cyclic properties of the trace and complex conjugation, we obtain that $tr(U\u0302\u2020|k\u27e9\u27e8n|)tr(U\u0302|m\u27e9\u27e8l|)=tr(U\u0302|n\u27e9\u27e8k|)*tr(U\u0302|m\u27e9\u27e8l|)$. The question that now arises is that if knowing all the values of $tr(|k\xaf\u27e9\u27e8l\xaf||m\u27e9\u27e8n|)$ (that we can obtain since we only use $(f1,\u2026,fN)$ and $(f\xaf1,\u2026,f\xafN)$), we can retrieve $tr(U\u0302|m\u27e9\u27e8l|)$. We can. We need to notice that $tr(U\u0302|m\u27e9\u27e8l|)$ is a complex number, so knowing its polar form is enough. We see that $tr(|l\xaf\u27e9\u27e8l\xaf||m\u27e9\u27e8m|)=|tr(U\u0302|m\u27e9\u27e8l|)|2$. Thus, we obtain the modulus of the complex number. We can see now that $tr(U\u0302|m\u27e9\u27e8l|)=tr(|l\xaf\u27e9\u27e8l\xaf||m\u27e9\u27e8m|)ei\varphi m,l$. So, only the phases are up to determine. Here is where the issue of the overall phase freedom intervenes. Since $U\u0302$ is unitary, we know there must exist *m*_{0}, *l*_{0} such that $|tr(U\u0302|m0\u27e9\u27e8l0|)|2=tr(|l0\xaf\u27e9\u27e8l0\xaf||m0\u27e9\u27e8m0|)>0$. We have the freedom to fix the phase $\varphi m0,l0=0$ due to the overall phase redundancy. In other words, we could always choose a global phase to cancel the phase $\varphi m0,l0$ so it is set to 0. Now, we can see that if we consider $tr(|l0\xaf\u27e9\u27e8l\xaf||m\u27e9\u27e8m0|)=tr(U\u0302|m0\u27e9\u27e8l0|)*tr(U\u0302|m\u27e9\u27e8l|)=tr(|l0\xaf\u27e9\u27e8l0\xaf||m0\u27e9\u27e8m0|)tr(|l\xaf\u27e9\u27e8l\xaf||m\u27e9\u27e8m|)ei\varphi ml$. Thus, we obtain that

Therefore, indeed we can retrieve the unitaries that conjugate $(f\xaf\u03021,\u2026,f\xaf\u0302N)$ to $(f\u03021,\u2026,f\u0302N)$. □

(**Theorem 5**). *Given a fermionic system with the mode set* $I={i1,\u2026,iN}$*. For any set of descriptors* $(f\xaf\u03021,\u2026,f\xaf\u0302N)$ *together with any Heisenberg state* $|\psi 0\u27e9\u27e8\psi 0|$ *and any non-empty subset of modes* $J\u2286I$*. The associated diagram of Fig. 1 commutes. In other words,*

*Proof.* We begin by expanding the right-hand side of the equation, by applying the definition of the ontic projection operator, obtaining

Applying now the definition of $\phi $ of Eq. (11), we obtain

where we have chosen as orthonormal basis for the operator space of the fermionic system with modes $j1,\u2026,jn$ the basis given by ${f\u0302l1\u2020\u2026f\u0302lm\u2020|\Omega \u27e9\u27e8\Omega |f\u0302pr\u2026f\u0302p1}$. Where, $l1<\cdots <lm,\u2009p1<\cdots <pr$ with $l1,\u2026,lm,p1,\u2026,pr\u2208J$, where $|\Omega \u27e9$ is the vacuum for the set of modes *I*, and thus also *J*.

We now turn to expand the left-hand side of the initial equation by applying the definition of the epimorphism and choosing the same basis as before, but for the global set of modes {1,…, N}, we obtain

The main difference with the expanded version of the right-hand side is the values of the vectors $l\u2192,p\u2192$; in the first expression, they can have length up to *n* with the values being elements of *J*, while in the second case they can have lengths up to *N* with the values being elements of *I*.

Now, let us take the fermionic partial trace. In order to ease up on notation, let us reorder our set of modes such that *J* corresponds to the first *n* elements of our set of modes *I*. First, let us trivially observe that if *J *=* I* the equality holds, so, now we consider *J* a strict subset of *I*. Using the reordering, we can see how the expression of the left-hand side becomes

where $p\u2192,l\u2192$ are vectors of elements of *J* up to length *n* and the vectors $u\u2192,v\u2192$ are vectors of elements of $I\J$ up to length *N* − *n*. Once we have the expression of this form, we first notice that *ρ* will be a parity superselected state, since *ρ*_{0} is and the unitary is also, and the trace of an odd operator is always zero. With this notion, we can apply the properties of the fermionic partial trace^{11} where we need $u\u2192=v\u2192$ for the component not to vanish. Thus, we obtain that the left-hand side becomes

Moreover, we can focus now only on the coefficients to see that we can pull the sum over the vector $u\u2192$ inside. By using the cyclic property of the trace and then using its linearity, we obtain that the coefficients of the left-hand side are

Now, we can use the fact that $f\u0302l1\u2020\u2026f\u0302lm\u2020f\u0302u1\u2020\u2026f\u0302us\u2020|\Omega \u27e9$ can be written as $|l\u27e9\u2227|u\u27e9$, thus, obtaining that the coefficients of the left-hand side can be written as

we have that $\u2211u|u\u27e9\u27e8u|=I$ and $O\u0302\u2227I=O\u0302$, thus we obtain indeed that the left-hand side is $\u2211p\u2192,l\u2192tr(U\u0302\u2020(f\u0302l1\u2020\u2026f\u0302lm\u2020|\Omega \u27e9\u27e8\Omega |f\u0302pr\u2026f\u0302p1)U\u0302|\psi 0\u27e9\u27e8\psi 0|)f\u0302l1\u2020\u2026f\u0302lm\u2020|\Omega \u27e9\u27e8\Omega |f\u0302pr\u2026f\u0302p1$. Precisely as the right-hand side, proving the theorem. ◻

**Proposition 6.** *If one defines* $q\u0302j=12(\sigma \u0302jx+i\sigma \u0302jy)$*, then* $q\u0302j,q\u0302j\u2020$ *are generators of the N-qubit algebra.*

*Proof.* It is well established that $Ij,\sigma \u0302jx,\sigma \u0302jy,\sigma \u0302jz$ are generators of the *N*-qubit algebra. Moreover, since $(\sigma \u0302jx)2=Ij$ and $\u2212i\sigma \u0302jx\sigma \u0302jy=\sigma \u0302jz$ we can say that $\sigma \u0302jx,\sigma \u0302jy$ are generators of the *N*-qubit algebra. We will see how $q\u0302j,q\u0302j\u2020$ are enough to generate such operators that generate the algebra, thus, becoming generators. We first note that $q\u0302j+q\u0302j\u2020=\sigma \u0302jx$ and finally, $\u2212i(q\u0302j\u2212q\u0302j\u2020)=\sigma \u0302jy$. ◻