Special entangled fermionic systems and exceptional symmetries

Special fermionic systems entered the realm of quantum chemistry in the seventies in the work of Borland and Dennis in the form of a toy model. This work was leading to a detailed study of the N-representability problem by Klyachko. The topic then has been reconsidered in the light of entanglement theory boiling down to the notion of entanglement polytopes. Recently building on certain properties of such special fermionic systems, a connection between the coupled cluster method and entanglement has been established. In this paper we show that precisely such a special class of systems also provides an interesting physical realization for structures related to the Lie algebras of exceptional groups. This result draws such exotic symmetry structures under the umbrella of entangled systems of physical relevance.


Introduction
In recent developments of entanglement theory special fermionic systems have played an important role. The interest in studying such systems dates back to 1972, to the classical work of Borland and Dennis on the N -representability problem [1]. They were studying three fermion systems with six single particle states, and have shown that the eigenvalues of the one particle reduced density matrices are subject to special inequalities [1,2]. It turned out that these inequalities for the one particle marginals hold if and only if these marginals are coming from a pure three fermionic state. Hence these inequalities provide a sufficient and necessary condition for the pure state B Péter Lévay pplevay@gmail.com 1 MTA-BME Quantum Dynamics and Correlations Research Group, Department of Theoretical Physics, Budapest University of Technology and Economics, 1521 Budapest, Hungary 3-representability of a one-particle density matrix. After many years as a first new result in this field in 2006 Klyachko solved the pure state N -representability problem [3,4]. Thanks to this work it has turned out that the constraints for the eigenvalues of the one particle marginals define a convex polytope. Surprisingly the mathematics of such polytopes is inherently connected to quantum entanglement theory [5,6]. These observations are incapsulated in by now well-known terms terms like: entanglement polytopes and generalized Pauli principle.
Interestingly the system originally studied by Borland and Dennis also turned out to be very useful in another branch of entanglement theory: entanglement classification. The special role of three fermion systems with six single particle states in entanglement theory was first observed in Ref. [7]. It has been pointed out that for this system the so called SLOCC entanglement classes [8] can explicitely be given thanks to special mathematical structures, called Freudenthal triple systems [9]. These structures are familiar to mathematicians from studies concerning exceptional groups [10]. Later further connections between other entangled systems and Freudenthal systems have been discovered and the underlying correspondence has been elaborated [11]. This correspondence has recently been used to shed some light on connections between entanglement theory and the coupled cluster method [12]. In the meantime such special entangled systems have also been identified as embedded ones lurking within the realm of fermionic entanglement theory [13][14][15]. However, in these studies the precise mathematical connection between special entangled systems and exceptional symmetries have never been explored. The aim of the present paper is to fill in this gap, and to provide the first modest step in the direction of understanding exceptional Lie algebras generating exceptional symmetries, in fermionic entanglement terms.
As is well-known a Lie algebra is a mathematical structure describing infinitesimal symmetry transformations of a physical system which is exhibiting continuous symmetry properties described by the corresponding Lie group. The basic building blocks of Lie algebras are the simple Lie algebras. Over the field of complex numbers C the finite simple Lie algebras have been classified by Cartan, Killing and Dynkin [16]. The complete classification (neglecting exceptional isomorphisms) is given by the list of four infinite series A n , n ≥ 1, B n , n ≥ 2, C n , n ≥ 3, D n , n ≥ 4, n ∈ N, and five exceptional cases E 6 , E 7 , E 8 , F 4 , G 2 in terms of Dynkin diagrams [16]. Every simple complex Lie algebra has a unique real form whose corresponding centerless Lie-group is compact. These are also classified by Dynkin diagrams and usually in the same "ABC D E F G manner. For a physicist these compact groups of the four infinite series A n , B n , C n , D n are connected to well-known matrix Lie groups SU (n + 1), SO(2n + 1), Sp(2n), SO(2n) respectively. All of these matrix groups are leaving invariant some extra structure of physical relevance. SU (n + 1) is leaving invariant a Hermitian inner product and the associated norm of C n+1 familiar from quantum theory. SO(2n + 1) and SO(2n) being the rotation groups which are leaving invariant an ordinary scalar product and hence the associated length of a vector in R 2n+1 and R 2n . Finally Sp(2n) is leaving invariant a symplectic form (an alternating bilinear form) on R 2n regarded as a phase space. R 2n equipped with this structure is familiar as the symplectic geometry of classical Hamiltonian mechanics expressed in terms of Poisson brackets.
Then the question is: what kind of extra structures of physical relevance the exceptional symmetries of the corresponding groups are leaving invariant? Moreover, what kind of physical systems they are describing in a natural manner? There have been many possible answers to these question addressed in the physics literature. In this paper we would like to add an entanglement based physical system, which realizes its symmetries in an exceptional manner, to the list. We coin a term for such systems: special entangled embedded fermionic systems. We show that for them the extra structures are of the following kind. One of them is a well-known one: a symplectic structure which is already known from classical mechanics. The second one however, is quite unusual: a triple product. Indeed, multiplying two objects to obtain a third one is familiar to everyone, but in the case of a triple product one should multiply three object to get a fourth one! To cap all this we will see (Sect. 6.) that these structures are interrelated in an intricate manner. As one suspects these rather exotic structures are well-known to mathematicians [9,17]. However, to physicists having familiarity with fermionic systems they are definitely not. This paper is aimed to take the first step in the direction of drawing such symmetries under the umbrella of entangled systems of physical relevance.
The organization of this paper is as follows. In Sect. 2 we summarize some background material needed for our considerations. Fermionic systems with an even number of modes N = 2k are defined here. Our special entangled fermionic systems will be ones with six modes (k = 3). In Sect 3. we describe them as spinors in either the positive or negative chirality representation. There are special transformations acting on such fermionic states. These are called local stochastic entanglement transformations associated with classical communication [8]. The orbits of such transformations will be called SLOCC entanglement classes. Such transformations are considered in Sect. 4. Since we are interested in embedded systems, here we also describe how qubit systems can naturally be embedded into our fermionic ones via applying the so called single occupancy representation.
In Sect. 5. we list our five special entangled embedded femionic systems. They will be the ones corresponding to the groups D 4 , E 7 , E 6 , F 4 , G 2 , four of them 1 being exceptional. The mathematical basis of this correpondence will be elaborated in Sect. 6. As the main result of this paper here it will be shown that these systems exhibit a symplectic triple system structure. In proving this we explicitely give the form of the triple product and the symplectic form, and show that these are satisfying the basic axioms [17] of such systems. Curiously in order to establish this result an interesting mapping between states an operators shows up. Physically this means that the entangled states and the entanglement transformations acting on them merge into a unified mathematical structure. As a byproduct of these results in Sect. 7 we introduce a dual state, and an entanglement invariant. These objects are serving as physically relevant quantities which enable an explicit entanglement type identification.
In Sect. 8. we realize that by adjoining an extra single qubit SLOCC subgroup and an extra complex degree of freedom one can construct nonlinear realizations for our groups D 4 , E 7 , E 6 , F 4 , G 2 . These groups enlarge the SLOCC group substantially capable of incorporating transformations changing the value of the basic entanglement measure hence the entanglement type. Though the basic ideas of this realization are already available in the literature [18] we did not attempt here to give an explicit construction of it and explore its physical consequences. We postpone such elaborations for future work. The conclusions and some comments are left for Sect. 8. Here we also speculate on the role of the largest exceptional group E 8 for which we did not manage to give a simple fermionic entangled system interpretation. Hopefully in future investigations we will be able to add this item to our list.

Fermionic systems with an even number of modes
In the following we need the formalism as developed in [19] for fermionic systems with an even number of modes. Let k ∈ N and let V be an N = 2k dimensional complex vector space and V * its dual. We regard V = C N with {e I }, I = 1, 2, . . . N the canonical basis and {e I } the dual basis. Elements of V will be called single particle states. The number N = 2k is called the number of modes. Equipped with a Hermitian inner product ·|· our vector space (V , ·|· ) can be regarded as the Hilbert space of single particle states.
We also introduce the 2N dimensional vector space An element of V is of the form x = v + α where v is a vector and α is a linear form. According to the method of second quantization to any element x ∈ V one can associate a linear operator O x acting on the fermionic Fock space F N which is represented as a direct sum of m-particle subspaces F (m) , m = 0, 1, . . . N . Explicitly, to the basis vectors of V one associates fermionic creation and annihilation operators satisfying the usual fermionic anticommutation relations p I n J + n J p I = δ I J , p I p J + p J p I = n I n J + n J n I = 0.
Since a Hermitian inner product on V is at our disposal one can regard p I as the Hermitian conjugate of n I . In this way one arrives at the usual operators familiar from the literature: f I † ≡ p I and f I ≡ n I . However, for notational simplicity we refrain from this notation and we rather use the p I s and n I s in the text. Now the direct sum structure of F N in this formalism is made explicit as follows. First we define the vacuum state |vac ∈ F N by the property n I |vac = 0, Then F N is a vector space of dimension 2 N which is spanned by the basis vectors created from the vacuum |vac , p I 1 p I 2 · · · p I m |vac , Hence According to (6) an arbitrary state |ϕ ∈ F N can be written in the form Here the mth order totally antisymmetric tensors ϕ (m) I 1 I 2 ···I m encapsulate the complex amplitudes of the m-"particle" subspace and we introduced the new notation Before proceeding we should warn the reader that the formalism developed so far is not the canonical one familiar to physicists. This formalism is based on the work of Elie Cartan [20] and Claude Chevalley [21] called "The theory of spinors", which has been introduced to the physics literature by Trautman, Budinich and others [22]. Though it had turned out to be extremely useful for studying fermionic entanglement [14,15], our treatise is still not widespread among physicists working in the field of entanglement theory. We decided to use it due its mathematical convenience for exploring connections with exceptional groups our main concern here. In order to orient the reader for connecting our setup to a one which should be well-known here is an alternative way of looking at the structure of F N . In this more familiar picture one describes the 2 N basis vectors of F N as hence |00 . . . 0 ≡ |vac , |10 . . . 0 ≡ p 1 |vac etc. Let us now define 2N Majorana fermion operators as follows These operators are satisfying We will refer to the number 2N as the number of Majorana modes. One can represent the c μ via the N -fold tensor product 2 of Pauli matrices (with identity matrices in the remaining slots implicit) as In this way we arrive at the familiar Jordan-Wigner representation of Majorana operators.
Define now the chirality operator as For m even (odd) the states of Eq. (5) are eigenvectors of with eigenvalues +1(−1). The corresponding eigensubspaces of will be denoted by F ± N . One can also refer to their elements as spinors of positive or negative chirality [20][21][22].
Later it will be rewarding to introduce instead of the canonical labeling of single particle states as 1, 2, 3, 4, . . . 2k − 1, 2k, two other labelings. The first one will be called the occupied non-occupied labeling, the second the so called odd-even labeling. In the first case one defines a relabeling according to the formula (1, 2, . . . k, k + 1, k + 2, . . . , 2k) ↔ (1, 2, . . . , k, 1, 2, . . . , k) (14) in the case of the second one has The first labeling is useful for studying entanglement and the coupled cluster method [12]. In this context the 1, 2, 3, . . . states can be regarded as the occupied and the 1, 2, 3, . . . states as the non occupied ones. The second one is important when considering qubit systems embedded into fermionic ones. It is well-known that k-qubit systems [15,[23][24][25] can be embedded into fermionic ones with 2k modes. In this case the pairs j j will be referring to the two possible states of the jth qubit in the computational basis. We will give a more detailed description of these embedding procedures in the next sections.

Special fermionic systems with six modes
The states of our special fermionic systems we will be interested in are represented by the state vectors of the spaces F ± 6 . Hence in the following we will be exploring the special properties of fermionic systems with six modes (or twelve Majorana modes) with definite chirality (k = 3, N = 6, 2N = 12). Explicitly Notice that the dimension of F 6 is 2 6 = 64 and the fermionic Fock space splits into two 32 dimensional subspaces of positive (with representative |ϕ + ) and negative (with representative |ϕ − ) chirality. The 32 complex amplitudes of the representative cases split into two scalars η and ξ , two antisymmetric tensors x I J , y I J two vectors u I , v I , and a totally antisymmetric tensor of rank three z I J K where I , J , K = 1, 2, . . . 6.
Hence in the first case we have a split of dimensions 32 = 1 + 15 + 15 + 1, and in the second 32 = 6 + 20 + 6. Note that there exists an intertwiner between the positive and negative chirality representations such that This intertwiner relates the amplitudes via a bijection with explicit form [15] x where for later convenience the left hand side is labelled (x I J , y I J ) by I , J = 1, 2, . . . 6, and for the right hand side we applied the occupied-non occupied labeling of Eq. (14). Hence for example x 25 = z 231 = z 264 = −z 246 , y 54 = u 3 = u 6 etc.
One can also choose another intertwiner for connecting the positive and negative chirality representation. A one which is using the odd-even labeling is given by

Entanglement transformations, and embedded qubits
Until this point we have been discussing states. However, focusing on states is merely half of the physically interesting story, since in order to connect our mathematical considerations to physical phenomena we also have to discuss observables, and symmetry transformations, or even a special class of such transformations which generate interesting dynamics. In order to do this we have to describe special subsets of operators as infinitesimal generators of transformations acting on the fermionic Fock space. We will be interested in a special class of such transformations. Such transformations will form a group. The orbits of these transformations are states that can be converted to each other. Such orbits form equivalence classes that we will call entanglement classes.
A natural set of infinitesimal transformations that immediately springs into ones mind is the set of transformations generated by self-adjoint operators i.e. observables. The exponentials of such operators generate unitary symmetry transformations. The canonical example of that kind is the Hamiltonian as an observable associated to energy. The Hamiltonian then, via exponentialization, gives rise to the operator of time evolution. More generally one can consider a local set of unitary operators, that are acting merely on a chosen subsystem of our fermionic one. In this case the orbits regarded as entanglement classes are certain local unitary orbits represented by particular states. For an example (of a special system wich is a subset of our system with six modes showing up in the previous section) of such a classification of entanglement classes see Ref. [25]. Such local unitary transformations, combined with classical manipulations, forms the LOCC transformations [26,27] (local operations and classical comminucation).
However, in entanglement theory it has turned out that confining our attention to such subsets of transformations is too restrictive, and by relaxing the assumption of unitarity a more general set of transformations should be considered [8]. Such transformations are called SLOCC transformations (stochastic local operations and classical communication). In the fermionic context later even a further generalization has been introduced the GSLOCC transformations [15] (generalized SLOCC transformations). Now the set of operators we will be interested in is of the form Here A, B, C are N × N complex matrices with B and C skew-symmetric and for the repeated indices summation is understood. Our special choice is dictated by the convenience that via forming ≡ e s we obtain special group elements capable of representing GSLOCC transformations. Indeed, it can be shown [15] that = e s ∈ Spin(2N , C) taken together with a nonzero complex number λ is generating transformations of the form representing GSLOCC transformations. 3 First we will be interested in the Spin(2N , C) subgroup represented by elements of the form (1, ).
In order to understand what type of transformations are showing up in this mathematical setting let us consider the particle number conserving subgroup of the generalized SLOCC group which is obtained by setting B = C = 0 in Eq.(23). Then we get The action of on a state of the m-particle subspace where with G L(N , C) being the group of invertible N × N complex matrices. This group is known to represent SLOCC transformations on N state systems [8].
Since we have choosen N = 2k our formalism is also capable of describing embedded k-qubit systems [15]. In order to see this first we write the states of the k-fermion Here for the special case m = k we introduced a special notation z I 1 I 2 ...I k for the amplitudes of |z . Notice that for the k = 3 case, our main concern here, we get back to the amplitudes z I J K showing up in Eq. (17). Such states are precisely the Fock space representatives of the three fermionic ones with six single particle states studied by Borland and Dennis in their classical paper [1]. Under the subgroup of (25) the amplitudes of this state transform as where Hence in this special case the (25) subgroup of transformations coming from the group Spin(4k, C) with B I J = C I J = 0 will produce an SL(2k, C) subgroup comprising SLOCC transformations with unit determinant. Now we employ our odd-even labeling of Eq. (15).
Now the operators 4 the generators of n copies of the group SL(2, C). Clearly on the basis states (the computational basis for the j-th qubit) these operators give rise to the usual Pauli spin matrices. In order to see the physical meaning of our single occupancy representation [15] we express this unusual form of the Pauli operators in terms of the Majorana operators c μ , μ = 1, 2, . . . 4k introduced in Eqs. (10)- (11). On F 2k we have Notice now that on K s we also have c 4 j−3 c 4 j = −c 4 j−2 c 4 j−1 etc. hence on K s one can alternatively choose This result [19] connects the single occupancy embedding [15,25] of qubits into fermions to the influential work of Kitaev [28] on toric codes popular in possible implementations of quantum computing in solid state physics. For our considerations the value of this formalism is that we can now relate this action of the Pauli operators to a special subset of the infinitesimal fermionic SL(2k, C) transformations acting on F (k) in the (29) form. Indeed, these are transformations, characterized by a 2k × 2k matrix S of the (30) form, that leave the subspace K s invariant. Looking at Eq. (29) it is easy to see that such transformations in the (15) odd-even representation can be organized into a block diagonal matrix of the form where the 2 × 2 blocks of S ( j) of S are of the form with the Pauli operators, in the Kitaev representation discussed above, being the generators of this group action. Now the transformation |z ψ → |z ψ gives rise to the one which is just an SL(2, C) ×k action for k-qubits. Let us finally look at Eq. (24). Clearly after taking also into account transformations of the form (λ, ) where λ ∈ C × G L(1, C) gives rise to the G L(2, C) ×k action which represents SLOCC transformations for k-qubit systems [8]. Hence we can conclude that in this single occupancy Kitaev representation GSLOCC transformations for N = 2k fermionic systems incorporate also SLOCC transformations for k-qubit ones. In the following we will be interested only in the k = 3 case. We will see that such systems are very special ones, namely in this case one can define a mapping between states and operators representing entanglement transformations. Our aim is to show that such a mapping gives rise to structures well-known from mathematical constructions of exceptional groups. In this manner we will be able to establish a nice connection between embedded entangled fermionic systems and exceptional symmetries.

Embedded systems of the k = 3 case
As a starting point let us consider a fermionic system represented by a state in the negative chirality (17) form. Our choice for the negative chirality representation is dictated by convenience. Using the intertwiner (18) one can convert mathematical data to the (16) positive chirality one as well. Though such a conversion via an intertwiner obviously will yield a very different physical interpretation, it will turn out that the mathematical structures related to exceptional groups can be revealed in any of such representations. Now on the one hand elaborations are more easy in the (16) representation, on the other hand as far as physics is concerned the structure of embedded systems is manifesting itself nicer in the (17) one. But this is not a problem since using the intertwiner the physically relevant systems based on (17) will be also exhibiting the mathematical structures we are interested in. In the following first we give a list of the physically relevant embedded systems. They will be the main actors of our further considerations. Then after presenting some motivation for the structures we are intending to study, in the next section we turn to their mathematical elaboration.
1. According to Eq. (17) our basic system is a fermionic one with six modes. The representative state containing 1, 3, 5 "particle" states is encapsulated in the 32 = 6 + 20 + 6 complex amplitudes (u I , z I J K , v I ). The GSLOCC group is acting on this space in the manner as shown in Eq. (24) where (λ, ) ∈ G L(1, C) × Spin(12, C).
2. Three fermion states with six single particle states are naturally embedded into our system. The corresponding states are of the (28) form As we have already discussed, the GSLOCC group action is represented on this subspace by the action of the SLOCC group in the (24) and (29) way where (λ, ) ∈ G L(1, C) × SL(6, C). Explicitly, the amplitudes transform as 3. Our six dimensional vector space V of single particle states is equipped with a Hermitian form. One can put further structures on this space. A structure that has also appeared in the literature [29] is a sympectic one, which should be familiar from classical mechanics. As is well-known this can be represented by an alternating bilinear form ω, having the property ω(u, v) = −ω(v, u), where u, v ∈ V . The group which is leaving invariant this form is the symplectic group.
One can then also consider entanglement SLOCC transformations that also leave invariant ω. In this picture two states can be regarded to be equivalent iff they can be mutually converted to each other by transformations representing physical manipulations that are leaving invariant the extra structure of the single particle space. However, it can be shown [30] that on the space F

4.
One can also consider three qubit systems as embedded ones in the single occupancy representation. Here the 8 complex amplitudes ψ αβγ of three qubits are embedded into the three fermion state by the k = 3 version of (34), i.e. the relevant |ψ → |z ψ mapping. The transformation rule for the amplitudes is then with the group of SLOCC transformations being [8]

5.
There is still one more embedded system that can naturally be included in our list. This is the case of three bosonic qubits. In this case the three qubits are indistinguishable hence instead of the 8 amplitudes we are having merely 4. Explicitly a bosonic three qubit state |ψ b is having the form Clearly by using (34) this case can also be embedded in the |ψ b → |z ψ b manner in our basic (17) fermionic state with the prescription ψ 0 = ψ 000 , ψ 1 = ψ 001 = ψ 010 = ψ 100 etc. Since now the qubits are indistinguishable the SLOCC transformations form the group (G L(2, C) ×2 ) diag i.e. the diagonal action of three copies of G L(2, C). Explicitly we have i.e. due to indistinguishability now the same SL(2, C)s are acting on each qubit. Taking into account the arbitrary λ ∈ G L(1, C) degree of freedom we are left with the same G L(2, C)s acting on each qubit, i.e. the diagonal action. The entanglement classes (the SLOCC classes) of these systems are orbits under the action of the corresponding groups. The classification problem of entanglement types in this case boils down to finding representative states from each class. For the systems 1.-5. the structure of these classes (with representatives) follows a similar pattern. For example in the three qubit case 4. we have the fully separable class (|000 ), biseparable classes (|011 + |000 and permutations of the qubits), the so called W class (|001 + |010 + |100 ) and the famous GHZ class (|000 + |111 ). Hence we see that the number of nontrivial entanglement classes is four. This pattern survives for the remaining cases except for 5. where the biseparable class is missing.
For system 1. the classification of orbits under the group G L(1, C) × Spin(12, C) is also known as the classification of the corresponding spinors. This problem has been solved by Igusa [31]. The interpretation of the classification of spinors in terms of entanglement theory first appeared in [14,15]. The classification problem for three fermions with six single particle states has been solved by Reichel [32] in 1907, his classification has been interpreted as classification of entanglement types in [7]. It has also been shown that the different entanglement types can be related to the structure of cluster operators in the coupled cluster method [12]. The classification problem of three-qubit systems has appeared in [8] in 2000, however the corresponding mathematical problem named as the classification problem of trilinear forms is due to Le Paige [33] and dates back to 1881.
Notice that after applying the intertwiner one can embed our systems also into a state of the (16) form. Then looking at the corresponding amplitudes in Eqs.(20)-(21) the following interesting pattern emerges. We always have a single vacuum state and a single six fermion state. For the two (labeled by the matrix y) and four fermion states (labeled by the matrix x) we have the following numbers for the complex amplitudes of the corresponding cases: (1., 2., 3., 4., 5) ↔ (15, 9, 6, 3, 1). One can notice that the (15, 9, 6, 3, 1) numbers are the numbers of independent components of 6 × 6 antisymmetric matrices (15), 3 × 3 matrices (9), 3 × 3 symmetric matrices (6), 3 × 3 diagonal matrices (3), and 3 × 3 diagonal matrices with the same elements in the diagonal.
Let us now consider 3 × 3 Hermitian matrices with elements taken from R,C, H, i.e. the elements of these matrices are taken to real, komplex or quaternionic. Using overline for the corresponding conjugation one can give the structure to the matrices x and y as 3 × 3 matrices of the form These matrices have (15,9,6) real elements. Now we see that by using the complex extensions of these matrices, and using instead of (R, C, H) the triple (R × C, C × C, H × C we obtain matrices with (15,9,6) complex elements. For example in the quaternionic case one can consider the correspondence ⎛ Here on the left hand side α, β, γ ∈ C, a, b, c ∈ H ⊗ C, and overline means quaternionic conjugation. On the right hand side α, β, γ ∈ C, a, b, c ∈ Mat(2, C), i.e. on the right hand side these are 2 ×2 matrices.ã ≡ −εa T ε Here ε is the 2 ×2 matrix with components ε 12 = −ε 21 = 1, ε 11 = ε 22 = 0. One can then immediately verify that the 6 × 6 matrix on the right hand side is antisymmetric i.e. an element of Skew(6, C) having 15 independent complex components. On the other hand on the left hand side we have a Hermitian 3×3 matrix with elements being biquaternions (i.e. complexified quaternions). Hence we have the following sets of splits 32 = 1 + 15 + 15 + 1, 20 = 1 + 9 + 9 + 1, 14 = 1 + 6 + 6 + 1, (53) Then in the (53) case we get the following objects all of them described according to the similar (51) pattern! The remaining cases of course also fit into this scheme, since these are merely arising as restriction to the diagonal components of the relevant matrices. These considerations clearly indicate that we should be able to describe our systems in a unified framework. This framework is of course well-known in the mathematics literature. It is related to the construction of the so called Freudenthal triple systems [9,35] based on cubic Jordan algebras [34]. Note that the connection between entangled systems and the Freudenthal construction has already been studied in [11]. In the next section we would like to initiate a new way of looking at this construction based on embedded fermionic systems. In order to do this we have to introduce a quite unusual structure a triple product.

Special fermionic systems as symplectic triple systems
A symplectic triple system [17] is an R vector space over C equipped with a {·, ·} : R × R → C symplectic form 5  One can show that as a result of (57)-(58) we also have the property In order to obtain some insight into this rather abstract structure let us observe that for a fixed pair of vectors u, v the map is a linear one. Then in terms of this map the first axiom says that and the third This property shows that L uv is a derivation, since its action resembles the action of the ordinary derivative on a triple product. The (59) property shows that this linear map preserves the symplectic structure Clearly the second axiom of Eq.(57) of complicated appearance gives a very stringent constraint between the symplectic structure and the triple product. As a main result of this paper now we prove that the fermionic states (16) and (17) and their embedded subsystems define a symplectic triple system. In order to prove this it is enough to use the (16) positive chirality states, since this result can be transferred to the (17) negative chirality one via the use of the (18) intertwiner. In the following we use the notation for the basic state of Eq. we have We also introduce the notation To every fermionic state |ϕ let us now define an operator K ϕ which is of the (65) form as follows Using these expressions the explicit form of K ϕ is As we see this operator is quadratic in the amplitudes of the original state. This association of an operator to a state is a special case of a more general construction valid for N = 2k with k odd. This association is based on the so called moment map. For more information on the mathematical background of this map see Ref. [14]. One can generalize this construction further by noticing that even for a pair of states (|ϕ 1 , |ϕ 2 ) one can associate an operator of the above form. Indeed, in this case one can define Here I is the 6 × 6 identity matrix and the matrix indices are left implicit. One can also check that Tr A ϕ 1 ϕ 2 = −2ω 12 . Notice that Moreover one also has This property reminds us of Eq. (61). Moreover, since K ϕ 1 ϕ 2 is a linear operator, motivated by Eq. (60) it is worth considering the expression to define a new fermionic state associated to the three original ones (ϕ 1 , ϕ 2 , ϕ 3 ). This new |ϕ state has the (16) form where now the corresponding 32 amplitudes are of the form (η , y , x , ξ ). Similarly for the original states |ϕ α we have (η α , y α , x α , ξ α ) (α = 1, 2, 3). One can then calculate the explicit form of (η , y , x , ξ ) using the expressions (67)-(70) where now instead of the A, B and C matrices one has to use the (79)-(80) ones, and instead of (η, y, x, ξ) the set (η 3 , y 3 , x 3 , ξ 3 ) is to be used. With this process we have associated to the triplet (|ϕ 1 , |ϕ 2 , |ϕ 3 ) the new |ϕ fermionic state. Now we conjecture that the map is a triple product. In order to check this let us denote the 6×6 matrices of Eqs.(79)-(80) by A 12 ≡ A ϕ 1 ϕ 2 , B 12 ≡ B ϕ 1 ϕ 2 and C 12 ≡ C ϕ 1 ϕ 2 . Then the explicit form of our conjectured product is Let us refer to this product as Let us also define a symplectic form as where Now via a straightforward calculation one can check that with these definitions the axioms of Eqs.(56)-(58) are satisfied. Hence we have obtained the nice result that the complex 32 dimensional vector space F of fermionic states equipped with the (90) symplectic form and the triple product (89), i.e. the mathematical object (F, {·, ·}, [·, ·, ·]) is a symplectic triple system.
As we have already discussed using the (18) intertwiner one can convert our results obtained for the positive chirality representation to the negative chirality one which has a direct connection with the fermionic embedded systems related to physically relevant ones. For example if we are using the (17) negative chirality representation (now |ϕ = |ϕ − ) then instead of Eq. (78) for the components of K ϕ we obtain where Notice now that the (90) symplectic form is based on the polarization of the quantity −ξη − (x, y) hence by virtue of Eqs.(20)-(21) its new form can be rewritten as (now where Clearly after going through the same steps as in the positive chirality case we obtain again a symplectic triple system. However, since now we are in the realm of embedded fermionic systems we can also easily study symplectic triple subsystems corresponding to our physically relevant ones. Indeed, the triple products for these systems are just arising from the expressions above via doing suitable restrictions to the set of parameters: (v α , u α , z α ) where α = 1, 2, 3. For example in order to obtain a symplectic triple system for three fermion systems with six single particle states one just have to set u = v = 0. In this case instead of K ϕ one has to use K z and objects like K z 1 z 2 |z 3 to create the relevant triple product. The symplectic form in this case is 1 2 {z 1 , z 2 }. In this way we arrive at the nice result that all of our five embedded entangled systems discussed in Sect. 5.can be given the structure of a symplectic triple one. In the next section we elaborate further on the structure of such systems revealing their connection to exceptional groups.

Invariants, dual states and entanglement classes
Since symplectic triple systems have a general appearance in the following we refrain from the use of |ϕ for states and we simply refer to a state in the notation ϕ ∈ F where the notation F can be used for any of the prositive or negative chirality representations. However, for historic reasons connected to the literature on exceptional groups we opted to use the positive chirality one.
Let us introduce the following notation and the new quantity (a dual state cubic in the original amplitudes) (Recall Eq. (82) and (89).) Then the new form of (57) is with the remaining axioms left intact.
Let us now combineφ with the (90) symplectic form to create a polinomial which is quartic in the amplitudes of our fermionic state ϕ ∈ F.
It is useful to organize formally the (η, y, x, ξ) quantities into a matrix [10] then in this notation one has The explicit form of the quartic polynomial is where Pf(x) = is the Pfaffian of a 6 × 6 antisymmetric matrix. One can then show [14] that This means that our polynomial and the symplectic form is invariant under the nontrivial part of the SLOCC group G L(1, C) × Spin (12, C). Now the SLOCC classification of entanglement types for fermionic states (i.e. the classification of spinors [31]) can be described in the following elegant manner. Suppose we are given a state ϕ ∈ F. We read off the 32 complex amplitudes. Then we calculate the matrix of the form (100), the operator K ϕ of the form (77), and the quartic invariant of the form (102). Then one has to check which of the following five cases holds (see Table 1.).
Having determined the entanglement class one can be sure that via a SLOCC transformation our state ϕ can be transformed into some of the representatives showing up in the Table 1.
Notice that in the sixth column of Table 1. we have given the embedded qubit version of the unnormalized fermionic state. In arriving at this form we have used Eqs. (22) and (34). The I. class in this picture corresponds to the GHZ class of three qubits, since the representative state is the (X ⊗ X ⊗ X )(H ⊗ H ⊗ H ) (bit flip gates followed by Hadamard gates) transformed version (a SLOCC transformation) of the famous 6 GHZ state |000 + |111 . As far as the II. class is concerned its representative can be transformed by a bit flip gate I ⊗ X ⊗ I to the usual representative of the so called W -class i.e. |110 + |101 + |011 . The class III. representative is a biseparable state of the form (|11 + |00 ) ⊗ |1 . Finally |111 is representing the totally separable class. For a discussion of these classes see the original paper of Dür et. al. [8]. For a discussion on the geometric meaning of the fact that the same entanglement patterns appear for all of our special embedded entangled systems see Ref. [29].

Connecting fermionic systems to exceptional symmetries
Let us summarize what we have found in this paper in connection with special embedded fermionic systems. We have five such systems, all of them embedded in the fermionic system with 6 modes or 12 Majorana modes. In the positive chirality representation the number of complex amplitudes of their representative quantum states splits according to the pattern: 1 + σ + σ + 1 where σ = 1, 3, 6, 9, 15. The numbers 6,9,15 are related to the complexifications of 3 × 3 Hermitian matrices over R, C, H. The remaining numbers are arising as restriction to diagonal matrices. Such matrices are forming Jordan algebras [34] J F with F taken from the three possibilites R, C, H. There is also a construction of systems based on two copies of such Jordan algebras [34] and two scalars answering the 1 + σ + σ + 1 type split, called Freudenthal triple systems [9,10,35]. Such objects are featuring triple products and are used for constructions of exceptional groups. Since we have come accross these algebraic structures in connection with our fermionic systems, we are expecting that these structures can also be linked to exceptional groups in a natural manner.
According to Table 2. We observe that for four of our systems the number 2m+|k|+3 gives the dimensions of exceptional groups i.e. |g|. But we have to point out that there is an important caviat to this. Indeed, for the largest exceptional group e 8 the question mark shows that we did not manage to associate to this case any simple fermionic system. The reader familiar with Jordan algebras realizes that the last entry of our table should be associated with Jordan algebras coming from the complexification of octonion hermitian matrices, i.e. J O . Since the multiplication of octonions is non associative it is not at all surprising that a simple fermionic system is not capable of reproduccing this algebraic structure. In the last section of our conclusions we present some hints, how this case could also be incorporated into our formalism.
In closing this paper we will not present an explicit dictionary between the generators of the groups g 2 , so(8), f 4 , e 6 , e 7 and fermionic structures. In the following we only would like to reveal two different ways how the exceptional symmetries could be realized on our fermionic systems.
First we would only like to indicate that based on Ref. [18] in principle one should be able to obtain a nonlinear realization of exceptional symmetries in a following straightforward manner. We merely sketch the basic ideas details will be intended to be explored in future work.
The idea is to take our basic 32 dimensional fermionic system denoted by F and an extra copy of C. The direct sum of these spaces forms a 33 dimensional vector space. This will be the representation space for the exceptional group E 7 . Define now 32 + 32 extra generators acting on R ≡ F ⊕ C in a convenient manner [18]. Also take an extra copy of the one qubit SLOCC subgroup SL(2, C) with 3 generators E ± and H satisfying the well-known commutation relations and specify their action on R. An interesting feature of this realization is the action of the generator F on an element of R of the form (ϕ, ζ ) ∈ F ⊕ C. We have [18] These expressions show that the action of the generator E − is featuring two of the basic quantities needed for the classification of entanglement types showing up in Table 1. They are the quartic invariant D(ϕ) and the dual stateφ. Hence transformations generated by E − are clearly capable of changing the entanglement type of the state ϕ.
Hence the E − generator of the extra SLOCC-like transformation is representing some sort of global manipulation on our entangled system. It would be very interesting to explore the physical meaning of this transformation. Further observations on this structure reveal other interesting possibilites. Take for instance now the infinitesimal generators of the Spin(12, C) part of the SLOCC group (66 generators) and also specify their action on R. Then it is well-known [18,35] that one can decompose the adjoint representation of the exceptional group E 7 which has 133 generators with respect to the subgroup of Spin(12, C) × D where D refers to the subgroup generated by H in the following way Since D itself is a part of the extra SL(2, C) group, then the above decomposition of the adjoint of E 7 corresponds to a more general one with respect to the group Spin(12, C) × SL(2, C) of the form This gives the mathematical meaning of the formula 2m +|g|+3 giving the dimensions of the exceptional groups we have found in Table 2.
Indeed, similar decompositions also hold for the remaining groups E 6 , F 4 , D 4 and G 2 . The corresponding decompositions are of the form Moreover, though in this paper we did not manage to find a fermionic realization of the largest exceptional group a similar decomposition holds even in this case These results give rise to a second way of realizing exceptional symmetries. First notice that the common feature of these decompositions is the appearance of a new copy of SL(2, C) which also suggests a new role for our exceptional symmetries in a new type of system consisting of embedded fermionic ones and an extra qubit. The states of such physical systems are represented by elements of a vector space of the form F ⊗ C 2 .
As an example have a look at Eq. (111). Clearly the left hand side of SL(2, C) ×3 × SL(2, C) ⊂ D 4 shows that in this case F ⊗ C 2 contains the four-qubit state space on which the four qubit SLOCC subgroup SL(2, C) ×4 acts. Since D 4 SO(8, C) in this case the (111) decomposition suggests that the four qubit states space taken togeteher with the corresponding SLOCC transformations can alternatively be described via studying the adjoint action of the D 4 Lie algebra so (8)

Conclusions
In this paper we have been considering five special entangled quantum systems. They were embedded ones in the sense that they can be obtained by restriction from a basic one. This basic system was a fermionic one with six modes or twelve Majorana modes. The embedded systems contained in this system were: three fermions with six modes, three fermions with six symplectic modes, three qubits, and three bosonic qubits. These systems are called special because of their special symmetry properties. We observed that their symmetry properties are realized by simple Lie groups of exceptional type. We have established that four from our five systems is connected to the exceptional groups E 6 , E 7 , F 4 and G 2 . It is well-known in the mathematical literature that the exceptional symmetries are connected to the appearance of special structures left invariant by the correponding exceptional transformations. These structures are a symplectic form and a triple product connected to each other in a peculiar manner forming a symplectic triple system. As the main result of the paper we verified by a detailed explicit construction that our special embedded fermionic systems are form-ing such triple systems. We also suggested two possible ways how exceptional group actions should be realized on our embedded systems. Our results draw exceptional symmetry structures under the umbrella of entangled systems of physical relevance. We must point out however, that we did not manage to give a fermionic entanglement interpretation of the largest exceptional group E 8 . However, the results of Sect. 8. clearly show that this system should be considered as a basic one capable of incorporating all of the fermionic systems studied in this paper. Therefore the problem of finding the physical interpretation of the E 8 symplectic triple system is of basic importance. Clearly the construction of a fermionic Fock space representation of this E 8 case is not at all obvious. One solid piece of evidence is that if we use two copies of embedded fermionic systems one in the single and the other in the so called double occupancy representation [15] of qubits there is a fermionic code [19] featuring operators related to the Lie algebra e 8 . However, it is not at all clear how such constructs should be adjusted in a manner to be able to form a fermionic system containing our basic fermionic ones with six modes via a symplectic triple system structure.
There are two interesting possibilities left open by this paper. One of them is the explicit construction of exceptional symmetry generators acting on F ⊕ C, and the other is connected to F ⊗ C 2 .
In the first case what we have is just a fermionic entangled state ϕ ∈ F and an extra variable ζ ∈ C. ζ is related to the change in the quantities D(ϕ) andφ characterizing the entanglement type. This change is effected via the action of a special generator E − of the exceptional symmetry. See Eq. (106). Hence this exceptional symmetry generator is connected to global, i.e. entanglement type changing transformations. This is to be contrasted with the local nature of SLOCC transformations. This example already shows that in order to reveal the physical meaning of the exceptional symmetry transformations the explicit form of its group action on (F, ζ ) is of utmost importance.
In the second case we have a fermionic system and an extra qubit. Hence in this case we have a new type of entangled system with properties encapsulated in the m part of the exceptional Lie algebra. On the other hand, see Eq. (114), the h part is encapsulating the algebra of SLOCC entanglement transformations. However for realizing the exceptional symmetry transformations in the adjoint manner briefly discussed in Sect. 8. we are faced with the following problem. What is the physical meaning of commutators of the form [m, m] ⊂ h that have the strange interpretation of commutators of entangled states giving rise to SLOCC transformations? material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.