Next Article in Journal
Quantum Value Valuation Continuum
Previous Article in Journal
An Ultra-Energy-Efficient Reversible Quantum-Dot Cellular Automata 8:1 Multiplexer Circuit
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Tomographic Universality of the Discrete Wigner Function

1
Departamento de Física, Universidad de Guadalajara, Revolución 1500, Guadalajara 44420, Mexico
2
Departamento de Matemáticas, Universidad de Guadalajara, Revolución 1500, Guadalajara 44420, Mexico
*
Author to whom correspondence should be addressed.
Submission received: 6 November 2023 / Revised: 26 December 2023 / Accepted: 16 January 2024 / Published: 19 January 2024

Abstract

:
We observe that the discrete Wigner functions (DWFs) of n-partite systems with odd local dimensions are tomographically universal, as reflected in the delta function form of the DWF for any stabilizer. However, in the n-qubit case, this property does not hold due to the non-factorization of the mapping kernel, the explicit form of which depends on a particular partition of the discrete phase space. Nonetheless, it turns out that the DWF for some specific stabilizers, not included in the set used for the construction of the Wigner map, takes on the form of a delta function. This implies that the possibility of classical simulations of Pauli measurements in a given stabilizer state for qubit systems is closely tied to the experimental setup.

1. Introduction

The concept of quasidistribution functions, which emerged as a means to incorporate quantum corrections into statistical mechanics within two distinct contexts [1,2], has been successfully used to reformulate quantum mechanics in position and momentum (plane) phase spaces [3,4]. Subsequently, these concepts were generalized to more sophisticated geometries. In the framework of the phase space approach, operators suitable for a given quantum system are represented as real functions obtained from an invertible and covariant linear map w ^ ( Ω ) :
f ^ W f ( Ω ) = Tr f ^ w ^ ( Ω ) ,
where Ω is a point in the corresponding phase space M .
The most natural and widely used representation [5,6,7], commonly referred to as the Wigner correspondence [1], is self-dual. This self-duality allows for the treatment of both the states and observables in the same manner, making it particularly useful for studying the quantum classical correspondence [8]. However, the main drawback of the Wigner map lies in its representation of some states in the form of negative (quasi)distributions. In fact, it is noteworthy that non-negative Wigner functions (WFs) correspond only to the so-called Gaussian states, thereby restricting the type of positivity-preserving operations to Gaussian ones.
A discrete analog of the Wigner map, which is applicable for describing a finite number of particles with a local dimension that is a prime number, was introduced much later [9,10,11,12]. This discrete analog has recently found remarkable applications in the analysis of the classicality and the (classical) simulability of n-partite quantum systems [13,14,15,16,17,18,19,20]. In this context, the discrete phase space (DPS) M takes the form of a p n × p n grid, constituting an affine plane [21]. By labelling it with elements from the Galois field F p n , this DPS inherits the same geometric properties as the ordinary plane; e.g., parallel lines (not necessarily straight) do not intersect.
The construction of discrete quasidistributions is closely tied to the concept of mutually unbiased bases (MUBs) [22,23,24], which are eigenstates of disjoint sets of p n commuting monomials (stabilizers) formed by products of generalized Pauli operators. In power-of-prime dimensions, one can always arrange the operational basis formed by p 2 n 1 monomials (excluding the identity operator) into p n + 1 disjoint stabilizers [25]. In addition, it is possible to associate lines in the DPS with projectors onto elements of a certain (stabilizer) basis in a way that parallel lines correspond to different components of the same basis, while (single) intersecting lines are related to states belonging to the mutually unbiased bases [11]. It is noteworthy that in the prime-dimensional case, only straight lines can be related to MUBs. However, in power-of-prime dimensions, such an association exists with more complex geometrical structures, the so-called commutative curves [26,27]. This leads to a fundamental property of the discrete Wigner function: the tomographic condition. This condition means that summing the Wigner function over a line results in the probability that the system is in the state associated with that particular line [11,12].
The mapping kernel from the Hilbert space into the DPS can be recast as the sum of projectors onto elements of an appropriate complete set of MUBs. However, there are several inequivalent ways to form p n + 1 disjoint commuting sets, which correspond to different partitions of the DPS into non-intersecting curves. Such sets can be roughly characterized by their factorization properties, and both unitary equivalent and (globally) unitary inequivalent sets of stabilizers can be found. Consequently, the form of the discrete self-dual maps fundamentally depends on the chosen complete set of stabilizers [28]. For each state of the MUBs used in the construction of the mapping kernel, the Wigner function is represented as a delta function. The latter is equivalent to the tomographic condition [11,12,14,15,29].
On the other hand, it is known that the Wigner map for odd and prime local dimensions exhibits a property akin to the continuous position-momentum Wigner function: only stabilizer states correspond to non-negative Wigner functions [16]. This characteristic arises from the factorization of the self-dual mapping kernel for any partition of the DPS and results in the tomographic universality of the DWFs, which acquire a delta function form of any stabilizer state of any partition. This enables the use of the DWF as an indicator of quantumness, allowing the classification of quantum states based on their suitability for quantum computation speedup and estimations of the cost of classically simulating quantum circuits [18,19,20,30,31].
However, this property breaks down in the case of qubit systems [19,20,30]. The main reason for the loss of tomographic universality in n-qubit systems is the inequivalence among Wigner maps based on different sets of stabilizers [28,32]. This leads to the following question: is it possible for the DWF of an element of a stabilizer basis to take on the form of a delta function if that basis is not used for the construction of a Wigner map? In other words, are there states, aside from the elements of the MUBs fixing the map, that satisfy the tomographic condition?
In this paper, we show that in the n-qubit case, when a partition of the discrete phase space into a complete fixed set of stabilizers is given, it is possible to identify stabilizer states whose Wigner functions take on the form of delta functions along commutative curves belonging to other partitions. We analyze explicit results for the three-qubit case and discuss the implications of this observed property of n-qubit DWFs.
The structure of this paper is as follows: In Section 2, we recall the basic properties of n-qudit operations in the Hilbert space, particularly the construction of Abelian displacement operators and the relation to different types of MUBs. In Section 3, we discuss geometric structures related to MUBs in the discrete phase space which are employed for the construction of the Wigner map. In Section 4, we focus on the tomographic universality of n-qudit Wigner functions for odd local dimensions and on identifying the stabilizer states with non-negative delta-like Wigner functions in the n-qubit case.

2. The Generalized Pauli Group and Displacement Operators

Let us consider a system of n qudits, each with a local dimension that is a prime number p. It is convenient to relabel vectors on an orthonormal basis | k 1 , , k n , k i Z p in the corresponding Hilbert space H p n = H p n with elements of a finite field, | κ , κ F p n according to
| k 1 , , k n | κ , κ | κ = δ κ , κ , κ = i = 1 n k i θ i F p n ,
where { θ 1 , , θ n } is an abstract basis in the field F p n , considered as a linear space; the components k i , i = 1 , , n are obtained through the trace operation with a dual basis, { θ ˜ 1 , , θ ˜ n } [i.e., tr ( θ i θ ˜ j ) = δ i j ], with tr ( α ) = α + α 2 + + α p n 1 here for α F p n [33].
In even dimensions, there always exists a self-dual basis, i.e., a basis such that tr ( θ i θ j ) = δ i j , whereas in odd dimensions, there are almost self-dual bases, such that tr ( θ i θ j ) = q j δ i j , where q j is equal to 1 with one possible exception, e.g., q j = ( q 1 ) δ j n + 1 , for some q Z p .
The generators of the Pauli group P n can be expressed as:
Z ^ α = κ χ α κ | κ κ | , X ^ β = κ | κ + β κ | ,
Z ^ α = Z ^ α , X ^ β = X ^ β , Z ^ α p = X ^ β p = I ^ ,
where
χ α = ω tr α , ω = e 2 π i / p ,
are additive characters of the finite field,
χ α + β = χ α χ β , κ χ κ α = p n δ α , 0 .
The operators (2) satisfy the discrete counterpart of the Weyl form of the commutation relations:
Z α X β = χ α β X β Z α .
In an almost self-dual basis, the character of a product of field elements is factorized:
χ α β = Π j = 1 n ω α j β j q j ,
where a i = tr ( α θ i ) and b i = tr ( β θ i ) . This implies the factorization of Z ^ α and X ^ β into a product of cyclic operators [34,35]:
Z ^ = n = 0 p 1 ω n | n n | , X ^ = n = 0 p 1 | n + 1 n | , X ^ p = Z ^ p = I ^ , Z ^ X ^ = ω X ^ Z ^ .
according to
Z ^ α = Π j Z ^ q j a j , X ^ β = Π j X ^ b j .
The monomials { Z ^ α X ^ β , α , β F p n } form an operational basis in O p [ H p N ] and can be arranged in commuting sets of unitary displacement operators D ^ ( α , β ) obtained by Clifford transformations U ^ α , β of the set { Z ^ τ , τ F p n } :
D ^ α ( τ ) , β ( τ ) = Φ α , β τ Z ^ α ( τ ) X ^ β ( τ ) = U ^ α , β Z ^ τ U ^ α , β = D ^ α ( τ ) , β ( τ ) ,
where Φ α , β τ Φ α , β * τ = 1 is a phase factor. The operators (8) from each commuting set λ = { α , β } form an Abelian group,
D ^ α λ ( τ ) , β λ ( τ ) D ^ α λ ( τ ) , β λ ( τ ) = D ^ α λ ( τ + τ ) , β λ ( τ + τ ) ,
as a consequence of the relation Z ^ τ Z ^ τ = Z ^ τ + τ . This implies that the phases Φ α , β τ , where ( α , β ) λ , Φ α , β 0 = 1 , satisfy the equation:
Φ α , β τ Φ α , β τ = Φ α , β τ + τ χ α λ ( τ ) β λ ( τ ) = Φ α , β τ + τ χ α λ ( τ ) β λ ( τ ) ,
leading to the following condition for the labels of elements of λ :
tr α λ ( τ ) β λ ( τ ) = tr α λ ( τ ) β λ ( τ ) , τ , τ F p n .
The eigenstates | ψ κ λ of each commuting set { D ^ α λ ( τ ) , β λ ( τ ) D ^ λ ( τ ) , τ F p n } ,
D ^ λ ( τ ) | ψ κ λ = χ ( κ τ ) | ψ κ λ , κ F p n ,
form an orthonormal basis in H p n . The disjoint sets of displacement operators satisfy the relation:
Tr D ^ λ ( τ ) D ^ λ ( τ ) = p n δ λ , λ δ τ , τ + p n ( 1 δ λ , λ ) δ τ , 0 δ τ , 0 ,
and so the corresponding bases (12) are mutually unbiased,
| ψ κ λ | ψ κ λ | 2 = δ λ , λ δ κ , κ + p n ( 1 δ λ , λ ) .
Commuting sets of displacement operators { D ^ λ ( τ ) , τ F p n } , containing p n 1 monomials (excluding the identity operator), are usually called stabilizers and the respective eigenstates are the stabilizer states. The entire set of p 2 n 1 monomials can be partitioned into p n + 1 disjoint stabilizers in several locally inequivalent forms; i.e., they cannot be reduced to each other through local transformations. This is related to the possibility of constructing stabilizers with different factorization structures [36,37], which are related to the commutation condition between blocks of single-qudit operators (7), constituting commutative sets ( λ ) of direct products of n-generalized Pauli operators. Each monomial (up to a phase),
Z ^ α ( τ ) X ^ β ( τ ) = Π j Z ^ q j a j ( τ ) X ^ b j ( τ ) , τ F p n , α ( τ ) , β ( τ ) λ ,
can be divided into two parts, so that the first part contains k operators, corresponding to the ( j 1 , , j k ) qudits, and the second part contains the n k operators of the rest of the qudits. The stabilizer is factorized into at least two subsets if the products of the Pauli operators from the first blocks of all monomials commute between each other (the operators from all the second blocks satisfy the same condition). Some mutually commuting sub-blocks may exist inside the first or second blocks, etc. Thus, one can represent any stabilizer λ in the following form:
λ = { m 1 , m 2 , , m n } ,
where 0 < m 1 m 2 m n and j m j = n is the number of qudits in the j-th block that cannot be factorized into commuting sub-blocks any more. The partition λ = { 1 , 1 , , 1 } n corresponds to a completely factorized stabilizer and λ = { n } to a completely non-factorized one. It is worth noting that different partitions may have the same factorization structure.
The simplest complete set of stabilizers is labelled by linear functions:
α = τ , β = μ τ , μ F p n ; α = 0 ,
where the commuting sets α = τ , β = τ , α = 0 , β = 0 are always completely factorized. However, the factorization structure for the other values of μ depends on the number of qudits and the local dimensions.
More sophisticated dependences of pairs α λ ( τ ) , β λ ( τ ) over the parameter τ have the form [26,27]:
α λ ( τ ) = m = 0 n 1 α m λ τ p m , β λ ( κ ) = m = 0 n 1 β m λ τ p m ,
where the coefficients α m , β m F p n should satisfy conditions derived from (11), as, for instance,
m k tr ( α m β k ) = 0 , k = 0 , . . . , n .
It is worth noting that the displacement operators belonging to a given stabilizer λ can be expanded into the projectors of the corresponding basis vectors { | ψ κ λ , κ F p n } :
D ^ λ ( τ ) = κ χ ( κ τ ) | ψ κ λ ψ κ λ | .
The above equation allows us to represent any state in terms of the probabilities p κ ( λ ) = ψ κ λ | ρ ^ | ψ κ λ of projections onto elements of MUBs, i.e.,to tomographically reconstruct the state [11,12,14,15,22,29,38]:
ρ ^ = λ κ p κ ( λ ) | ψ κ λ ψ κ λ | I ^ = λ κ p κ ( λ ) 1 p n + 1 | ψ κ λ ψ κ λ | .

3. Phase Space Construction and the Wigner Map

The discrete phase space (DPS) [11,12,14,15,29] is a p n × p n grid, where the points ( α , β ) , α , β F p n label elements of monomials Z ^ α X ^ β . This DPS is endowed with a finite geometry [11,12,14,15,29,33] and admits a set of discrete symplectic transformations [39,40,41,42,43]. Similarly to the continuous case, the axes of the DPS are associated with the observables Z ^ α and X ^ β .
A separation of O p [ H p n ] into disjoint stabilizers corresponds to a partition of the DPS into non-intersecting (except at the origin), non-degenerate (containing p n different points including the origin) commutative curves { Γ = α ( τ ) , β ( τ ) , τ F p n } [26,27]. In other words, the points of every commutative curve Γ label a set of mutually commutative monomials:
{ Z ^ α ( τ ) X ^ β ( τ ) , [ Z ^ α ( τ ) X ^ β ( τ ) , Z ^ α ( τ ) X ^ β ( τ ) ] = 0 , α ( τ ) , β ( τ ) Γ } ,
such that condition (11) is satisfied. There exist p n + 1 commutative curves Γ λ in every DPS partition, so that the index λ , labelling curves within a given partition, takes p n + 1 distinct values. In addition, there are p n 1 parallel curves { Γ κ λ , κ F p n } for each curve passing through the origin. The whole bundle of p n parallel curves (called a striation) covers the DPS, and at every phase space point, p n + 1 curves (i.e., one curve from each striation) intersect. It is worth stressing that only the monomials, labelled by the points of the curves crossing the origin, commute among each other.
The set of rays (15) is the simplest example of such a partition, where the sets of parallel lines are explicitly defined as:
β = μ α + κ , μ , κ F p n ( p n sets ) α = κ , κ F p n , ( a single set ) .
Another way to establish a connection between the discrete geometry with the algebraic structures is to link the commutative curves with the eigenstates (12) of the corresponding stabilizers. This approach associates each basis { | ψ κ λ , κ F p n } with a striation of curves { Γ κ λ = α λ ( τ , κ ) , β λ ( τ , κ ) , τ F p n } , and the bases corresponding to different striations are mutually unbiased. In the conventional association, the curves that include the origin are put in correspondence with the states | ψ κ = 0 λ with unit eigenvalues (later called the unit stabilizer states). Such curves, Γ κ = 0 λ Γ λ , have the form (16) and will be called unit commutative curves.
Both n-qudit states and observables can be represented in the form of distributions in the DPS through a one-to-one map [16,39,40,41,42]:
W f α , β = Tr f ^ w ^ α , β ,
f ^ = p n α , β W f ( α , β ) w ^ ( α , β ) ,
where the Hermitian mapping kernel w ^ α , β = w ^ ( α , β ) is covariant under action of the discrete displacements (8), and is defined as
w ^ α , β = D ^ α , β w ^ 0 , 0 D ^ α , β ,
w ^ 0 , 0 = p n γ , δ D ^ ( γ , δ ) ,
so that
w ^ α , β = p n γ , δ χ ( γ β δ α ) D ^ ( γ , δ ) .
The kernel (23) is normalized according to
α , β w ^ α , β = p n I ^ ,
and satisfies the overlap relation:
Tr w ^ ( α , β ) w ^ ( α , β ) = p n δ α α δ β β .
It is convenient to recast Equation (24) as a sum of stabilizers corresponding to a given phase space partition,
w ^ 0 , 0 = p n λ τ D ^ λ ( τ ) I ^ ,
or, equivalently, as a plane superposition of projectors onto the elements of MUBs:
w ^ 0 , 0 = λ | ψ 0 λ ψ 0 λ | I ^ ,
where the index λ labels the stabilizers/projectors (12) within this partition. Geometrically, Equations (28) and (29) represent a summation along the curves that pass through the origin, i.e., of the form (16). In other words, each point ( α , β ) of the DPS is assigned a specific phase Φ α , β τ according to the unit commutative curve Γ λ that passes through that point and corresponding to the displacement operator D α λ ( τ ) , β λ ( τ ) . In Figure 1, we plot all of the unit commutative curves { Γ λ } and the distribution of the corresponding phases { Φ α , β , ( α , β ) Γ λ } along each curve in the case of three qubits for two different partitions.
By shifting w ^ 0 , 0 through the action of the displacement operators (23) to an arbitrary point ( α , β ) , one obtains
w ^ ( α , β ) = λ ( α , β ) , κ ( α , β ) | ψ κ λ ψ κ λ | I ^ ,
where the sum is performed over all the curves Γ κ λ passing through the point ( α , β ) .

4. Tomographic Universality of the Discrete Wigner Function

4.1. Tomographic Property for a Given DPS Partition

The Wigner function of any stabilizer state | ψ κ λ within a partition used for the construction of w ^ 0 , 0 is a convolution of characters computed at the points of the unit commutative curve Γ λ = γ λ τ , δ λ ( τ ) , which labels the elements of a stabilizer D ^ λ ( τ ) with the eigenvalues χ ( κ τ ) , κ F p n and fixes a particular stabilizer state according to (12),
ψ κ λ | w ^ α , β | ψ κ λ = p n τ χ γ λ ( τ ) β α δ λ ( τ ) χ ( κ τ ) ,
as it immediately follows from spectral decomposition (17) and Equation (25).
On the other hand, it immediately follows from decomposition (30) and the unbiasement condition (13) that
ψ κ λ | w ^ α , β | ψ κ λ = λ ( α , β ) , κ ( α , β ) | ψ κ λ | ψ κ λ | 2 1 = λ ( α , β ) , κ ( α , β ) δ λ , λ δ κ , κ ,
which is just a delta function along the curve Γ κ λ = { ( α λ ( τ , κ ) , β λ ( τ , κ ) ) , τ F p n } associated with the state | ψ κ λ :
ψ κ λ | w ^ α , β | ψ κ λ = τ δ α λ ( τ , κ ) , β λ ( τ , κ ) .
Equation (33) represents a flip side of the tomographic condition,
p n τ Tr w ^ α λ ( τ , κ ) , β λ ( τ , κ ) ρ ^ = p κ ( λ ) = ψ κ λ | ρ ^ | ψ κ λ .
In other words, the sum of the Wigner function of a state ρ ^ along a curve corresponding to a stabilizer state | ψ κ λ gives the probability of detecting the state ρ ^ in | ψ κ λ .
The properties (33) and (34) are automatically satisfied for the stabilizer states of the partition employed for the construction of the Wigner map. However, for stabilizer states from different partitions, Equations (33) and (34) do not necessarily hold. We will analyze the validity of these tomographic relations for both odd and even local dimensions.

4.2. Odd Local Dimensions

In the case of odd local dimensions, the kernel at the origin, Equation (28), takes on a simple form because the equation for phase (10):
Φ α , β ( τ ) Φ α , β ( τ ) = χ α λ ( τ ) β λ ( τ ) = χ α λ ( τ ) β λ ( τ ) , ( α , β ) λ ,
has an explicit solution for every stabilizer of any partition:
Φ α , β ( τ ) = χ 2 1 α λ ( τ ) β λ ( τ ) ,
where α ( τ ) = α ( τ ) and β ( τ ) = β ( τ ) due to (16). The primary property of phase (35) is its factorization, as follows from (5). As a consequence, the kernel w ^ 0 , 0 and thus w ^ α , β are factorized in an almost self-dual basis for any partition:
w ^ α , β = j = 1 n w ^ q j α j , β j ,
where
w ^ ( α j , β j ) = p 1 k , l ω ( β j k α j l 2 1 k l ) Z ^ l X ^ k .
This implies that the mapping kernel w ^ α , β does not depend on the particular partition of the DPS in the case of odd local dimensions. Such a statement may not be immediately obvious from the representation (29). This is because different partitions involve stabilizer states with different factorization properties (14) [28].
As a result, the tomographic conditions (33) and (34) are satisfied for any stabilizer state | ψ κ λ ( π ) belonging to an arbitrary phase space partition π , even if they are unrelated to the one employed for the construction of the Wigner map, i.e.,
ψ κ λ ( π ) | w ^ α , β | π | ψ κ λ ( π ) = τ δ α λ ( τ , κ ) , β λ ( τ , κ ) ,
where w ^ α , β | π is a Wigner kernel corresponding to some other partition π . Actually, the factorization of kernel (36) implies that the Wigner function of an arbitrary n-qudit state has exactly the same form for any DPS partition in the case of odd local dimensions.

4.3. Even Local Dimensions

The situation is significantly different for even local dimensions, p = 2 . In this case, none of the solutions of Equation (10) for the phase of the displacement operators for stabilizers labelled by points of commutative curves (16), and for n > 2 , admit a factorized form. A solution of the recurrence (10),
Φ α , β τ + τ = Φ α , β ( τ ) Φ α , β ( τ ) χ α λ ( τ ) β λ ( τ ) ,
under the condition,
Φ α , β 2 ( τ ) = χ α λ ( τ ) β λ ( τ ) , α λ ( τ ) , β λ ( τ ) Γ λ ,
is fixed by choosing signs of the phases corresponding to the elements of the self-dual basis { θ k , k = 1 , , n } ,
Φ α , β ( θ k ) = ± χ α λ ( θ k ) β λ ( θ k ) ,
where χ ( α ) = ( 1 ) t r ( α ) . For instance, for all positive signs, we have [43]:
Φ α , β τ = k = 1 n t k θ k = χ m = 1 n 1 α λ ( t m θ m ) β λ k = m + 1 n t k θ k k = 1 n c κ α , β ,
where c κ α , β = χ α λ ( t k θ k ) β λ t k θ k . The non-factorized form of the phase Φ α , β implies that the map w ^ ( α , β ) essentially depends on the choice of a complete set of disjoint stabilizers. In other words, it depends on the partition of the DPS into sets of commutative curves Γ λ .
In general, the Wigner function of a stabilizer state | ψ κ λ ( π ) from an arrangement π in the (core) partition π 0 has the form (see Appendix A):
ψ κ λ ( π ) | w ^ α , β | π 0 | ψ κ λ ( π ) = 2 n τ χ γ λ ( τ ) β + α δ λ ( τ ) χ ( κ τ ) Φ γ , δ π 0 ( τ ) Φ γ , δ π * ( τ ) ,
where γ λ ( τ ) , δ λ ( τ ) Γ λ ( π ) , Φ γ , δ π ( τ ) are the phases along Γ λ ( π ) , and Φ γ , δ π 0 ( τ ) are the phases assigned to each point γ λ ( τ ) , δ λ ( τ ) of the DPS according to the partition π 0 . Therefore, Equation (40) can be reduced to the form (33) if the phases Φ γ , δ π 0 ( τ ) and Φ γ , δ π ( τ ) are the same along the curve Γ λ ( π ) . Such curves will be referred as to Abelian curves. Obviously, all the curves that constitute a given partition are Abelian. However, there are Abelian curves outside of this partition, leading to a DWF of the corresponding stabilizer states that are delta functions:
ψ κ λ ( π ) | w ^ α , β | π 0 | ψ κ λ ( π ) = τ δ α λ ( τ , κ ) , β λ ( τ , κ ) , α λ ( τ , κ ) , β λ ( τ , κ ) Γ κ λ ( π ) ,
where Γ κ λ ( π ) is a curve associated with the state | ψ κ λ ( π ) in the partition π . It is worth noting that these Abelian curves may possess a factorization structure (14) not present in the original partition.
The density matrix of a stabilizer state corresponding to an Abelian curve not belonging to the core partition acquires a simple form:
| ψ κ λ ( π ) ψ κ λ ( π ) | = 2 n τ w ^ ( α λ ( τ , κ ) , β λ ( τ , κ ) | π 0 ) , α λ ( τ , κ ) , β λ ( τ , κ ) Γ κ λ ( π )
In other word, it is an equally weighted sum of the Wigner kernels in the core partition π 0 along the curve Γ κ λ ( π ) associated with the state | ψ κ λ ( π ) . However, the reconstruction expression (18) for a stabilizer state associated with an Abelian curve does not have the form of a convex sum of the core MUB projectors.
In general, searching for Abelian curves that do not belong to a given partition typically requires an involved analytical or numerical procedure. In addition, there is the same freedom (38) in the election of the phases of the displacement operators along such curves from a partition π to those of the core partition π 0 . In what follows, we give examples of Abelian curves in the case of three qubit systems.
Let us choose a partition of the DPS for the field F 2 3 corresponding to rays (15), and fix the phases of the stabilizers according to (39) in the same way for all slopes μ 0 . We do this by choosing the all positive solutions of the recurrence relation (37) for the basis elements, i.e., Φ α , μ α ( θ k ) = + χ μ θ k 2 , where it is always true that Φ α , 0 ( τ ) = 1 . This results in the following distribution of phases (where σ is the field generator):
Φ α , β = α β 0 σ σ 2 σ 3 σ 4 σ 5 σ 6 σ 7 0 1 1 1 1 1 1 1 1 σ 1 1 i 1 i i i 1 σ 2 1 i 1 i i i 1 1 σ 3 1 1 i i i 1 1 i σ 4 1 i i i 1 1 i 1 σ 5 1 i i 1 1 i 1 i σ 6 1 i 1 1 i 1 i i σ 7 1 1 1 i 1 i i i .
The partition corresponding to the rays has the factorization structure ( 3 , 0 , 6 ) , which means that there are three completely factorized bases and six completely non-factorized bases. In the three-qubit case, other admissible partitions include ( 1 , 6 , 2 ) , ( 2 , 3 , 4 ) and ( 0 , 9 , 0 ) [36,37]. Apart from the rays, there are two types of commutative curves constituting these partitions, specifically regular curves of the following form [26] (see also Appendix A):
β = f ( α ) = ϕ 0 α + ϕ 2 α 2 + ϕ α 4 , α = g ( β ) = ψ 0 β + ψ 2 β 2 + ψ β 4
These are characterized by non-degenerate values along the α or β axes and exceptional curves which are degenerated in both directions. These curves can be represented either in a parametric form or as relations of the type:
f ( α ) = g ( β ) ,
where α and β do not take all 2 n values of the field, but instead run only through an admissible set of points restricted by relations tr ( α σ 1 ) = 0 and tr ( β σ 2 ) = 0 , where σ 1 , 2 are some fixed elements of F 2 n .
We have found that for the primitive polynomial σ 3 + σ + 1 = 0 and the self-dual basis { σ 3 , σ 5 , σ 6 } , the regular curves (44) are Abelian under the positive sign selection (38) if the following relations between parameters ϕ 0 and ϕ are fulfilled:
1.
If ϕ = σ , σ 2 , σ 4 , then ϕ 0 = ϕ , ϕ 2 , ϕ 3 , ϕ 6 ;
2.
If ϕ = σ 3 , σ 5 , σ 6 , then ϕ 0 = 0 , 1 ;
3.
If ϕ = 1 , there are no ϕ 0 values producing Abelian curves, as can be observed in the example in Figure 1b, where the only Abelian curve is the ray α = 0 .
For instance, the bi-factorized curves in the partition ( 0 , 9 , 0 ) defined in Appendix B,
α = β + σ 6 β 2 + σ 3 β 4 ,
β = α + σ 3 α 2 + σ 5 α 4 ,
are Abelian, and the corresponding DWFs are
ψ κ ( α ; 1 , σ 3 ) ( 0 , 9 , 0 ) | w ^ α , β | ( 3 , 0 , 6 ) | ψ κ ( α ; 1 , σ 3 ) ( 0 , 9 , 0 ) = δ α , β + σ 6 β 2 + σ 3 β 4 + κ ,
ψ κ ( β ; 1 , σ 5 ) ( 0 , 9 , 0 ) | w ^ α , β | ( 3 , 0 , 6 ) | ψ κ ( β ; 1 , σ 5 ) ( 0 , 9 , 0 ) = δ β , α + σ 3 α 2 + σ 5 α 4 + κ ,
where | ψ κ ( α ; 1 , σ 3 ) ( 0 , 9 , 0 ) and | ψ κ ( β ; 1 , σ 5 ) ( 0 , 9 , 0 ) are eigenstates of the stabilizers
{ Z ^ β + σ 6 β 2 + σ 3 β 4 X ^ β , β F 2 3 } , { Z ^ α X ^ α + σ 3 α 2 + σ 5 α 4 , α F 2 3 }
respectively. Interestingly, both exceptional curves belonging to this partition are also Abelian.
On the other hand, the curves from the same partition ( 0 , 9 , 0 ) ,
α = σ 6 β + σ 3 β 2 + σ 5 β 4 ,
β = σ 2 α + σ 3 α 2 + σ 5 α 4 ,
are not Abelian; thus, the DFWs corresponding to the stabilizer states { | ψ κ ( β ; σ 2 , σ 5 ) ( 0 , 9 , 0 ) , κ F 2 3 } are not delta functions. In Figure 2, we plot the Wigner function for the unit stabilizer states corresponding to curves (46) and (50), respectively. Reconstruction of the stabilizer states | ψ 0 ( β ; 1 , σ 5 ) ( 0 , 9 , 0 ) and | ψ 0 ( β ; σ 2 , σ 5 ) ( 0 , 9 , 0 ) associated with curves (46) and (50) in the form of an expansion of the ray-related (15) MUB projectors fixed by the phases in (43) (and in terms of the corresponding probabilities (18)) is given in Appendix C. As it is expected, the obtained representations (A11) and (A14) are not convex sums of the projectors (A17) and (A18).
The above-mentioned curves become Abelian for a different assignment of signs to the phases (38). For instance, curve (50) is Abelian by choosing Φ α , μ α ( θ k ) = χ μ θ k 2 for any μ 0 . Under this assignment, the following regular curves (44) are Abelian:
1.
If ϕ = σ , σ 2 , σ 4 , then ϕ 0 = ϕ + ϕ 2 ;
2.
If ϕ = σ 3 , σ 5 , σ 6 , then ϕ 0 = ϕ + ϕ 2 ;
3.
If ϕ = 1 , ϕ 0 = 0 .
Additionally, curve (49) is Abelian under an assignment where the signs of Φ α , μ α ( θ k ) depend of the value of the slope μ .
In general, in the three-qubit case and the specific ( 3 , 0 , 6 ) DPS partition, certain choices of phase signs Φ α , μ α ( θ k ) can make all possible commutative curves Abelian. The inverse statement is however not true: for a given curve, it is not always possible to establish a choice of signs of Φ α , β ( θ k ) so that this curve becomes Abelian for fixed phases Φ α , μ α ( τ ) .

5. Conclusions

The Wigner mapping kernel can be constructed as a sum of projectors onto elements of a complete set of MUBs, which are eigenstates of disjoint stabilizers corresponding to a given partition of the discrete phase space into non-intersecting commutative curves. Different partitions lead to different factorization properties of MUBs, which are required for the mapping. For qudit systems with odd local dimensions, the Wigner kernel is factorized in the same form for all possible partitions. This results in tomographic universality, reflected in the delta function form of Wigner functions of any stabilizer state corresponding to any partition.
However, in the case of n-qubit systems, the mapping kernel is not factorizable, and its form depends on the chosen discrete phase space partition [28] and on the selected set of stabilizer states in this core partition, particularly on their phases. Consequently, the qubit discrete Wigner function is not tomographically universal for an arbitrary election of phases of the stabilizer states used for Wigner map construction. Nonetheless, this property is not completely lost. It is shown that for a given partition of the DPS, there are stabilizer states corresponding to commutative curves which do not participate in the construction of the mapping kernel, such that their DWFs are delta functions. This means that the probability of detecting these states by measuring in an arbitrary state is obtained by summing the DWF along the corresponding curve. This property is directly related to the feasibility of the classical simulation of Pauli observable measurements in such n-qubit states (see Theorem 2 in [18]). However, since the measurement procedure is codified in the partition of the DPS, only some specific stabilizer states, beyond those used in the measurement scheme, are classically simulable in a given experimental setup (for a fixed DPS partition). In other words, different experimental configurations may lead to different classical simulability outcomes for stabilizer states. This highlights the role of experimental design and setup in determining the classical or quantum nature of measurement outcomes. On the other hand, it is always possible to adjust a detection scheme so that the Pauli measurements in any particular stabilizer state can be described by a non-contextual hidden variable model. This supports a previously discussed finding [18,30] that the stabilizer states cannot be considered as a resource for quantum computation [13].
Although we have only analyzed a three-qubit system, the present approach is extendable to a higher number of qubits. In particular, one can expect that for any given partition, there are sign assignments such that the DWF of an arbitrary stabilizer state acquires a delta function form; i.e., it confirms the tomographic universality of n-qubit DWFs under the freedom of the phase choice. This conjecture has been verified by extensive numerical simulations and it can be justified by the same functional form of the recurrence Equation (37) for the phase (10) of any Abelian curve.

Author Contributions

Conceptualization, I.S. and A.B.K.; Software, E.C. and A.G.; Validation, A.B.K.; Formal analysis, I.S., E.C., A.G. and A.B.K.; Investigation, I.S., E.C. and A.G.; Data curation, E.C. and A.G.; Writing—original draft, A.B.K.; Writing—review and editing, I.S., E.C. and A.B.K.; Supervision, A.B.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

In this Appendix, we derive Equation (40) in the main text. Taking into account the expansion inverse of (17):
| ψ κ λ ψ κ λ | = 2 n τ χ ( κ τ ) D ^ λ ( τ )
and representation (25) of the Wigner kernel, one obtains
ψ κ λ ( π ) | w ^ α , β | π 0 | ψ κ λ ( π ) = 2 2 n τ γ , δ χ γ β + δ α χ ( κ τ ) Tr D γ , δ | π 0 D α λ ( τ ) , β λ ( τ ) | π ,
where the sum over τ denotes a sum over the points of the curve α λ ( τ ) , β λ ( τ ) Γ λ ( π ) . The displacement operators D γ , δ | π 0 and D α λ ( τ ) , β λ ( τ ) | π from the partitions π 0 and π , respectively, according to (8), (10) carry the corresponding phases Φ γ , δ π 0 and Φ α ( τ ) , β ( τ ) π * that guarantee the Abelian property (9). Thus, one gets
ψ κ λ ( π ) | w ^ α , β | π 0 | ψ κ λ ( π ) = 2 n τ γ , δ χ γ β + δ α χ ( κ τ ) Φ γ , δ π 0 Φ α ( τ ) , β ( τ ) π * δ α ( τ ) , γ δ β ( τ ) , δ ,
which after performing a sum over γ and δ , is reduced to Equation (40).

Appendix B

In this Appendix, we recall some basic properties of commutative curves in even local dimensions. Commuting sets constituted by 2 n different monomials (stabilizers) { Z ^ α λ ( τ ) X ^ β λ ( τ ) , τ F 2 n } are labelled by points of a discrete grid belonging to a non-singular curve (i.e., with no self-intersection) Γ λ that passes through the origin ( α ( 0 ) , β ( 0 ) ) = ( 0 , 0 ) Γ λ and satisfies
tr α λ ( τ ) β λ ( τ ) = tr α λ ( τ ) β λ ( τ ) , α λ ( τ ) , β λ ( τ ) Γ λ .
A general parametric form of such a curve is given by
α ( τ ) = k = 0 n 1 α k τ 2 k , β ( τ ) = k = 0 n 1 β k τ 2 k , α k , β k F 2 n ,
where the coefficients α k and β k are such that m k tr ( α m β k ) = 0 , which implies
Z ^ α ( τ ) X ^ β ( τ ) , Z ^ α ( τ ) X ^ β ( τ ) = 0 , α λ ( τ ) , β λ ( τ ) Γ λ .
It is worth noting that n appropriately chosen points of the curve, e.g., τ = { θ 1 , , θ n } , are enough to generate the entire curve. This is because the set of points of these curves forms an Abelian group where each element is of order two, and therefore the group is isomorphic to Z 2 n .
The regular curves, which are non-degenerate in at least one of the directions α or β (i.e., they take on all the values of the field in such a direction), can be represented in the explicit form [26]:
β = f ( α ) = k = 0 n 1 ϕ k α 2 k or α = g ( β ) = k = 0 n 1 ψ k β 2 k ,
where the coefficients ϕ k , ψ k F 2 n satisfy the following commutativity restrictions:
ϕ k = ϕ n k 2 k , ψ k = ψ n k 2 k , k = 1 , , n 1 2 ,
where . denotes the integer part. In particular, for even n values, the coefficients satisfy the additional restriction:
ϕ n / 2 = ϕ n / 2 2 n / 2 , ψ n / 2 = ψ n / 2 2 n / 2 .
The degenerate (or exceptional) curves are characterized by multiple appearances of the admissible points in both directions α and β . In other words, for every point ( α j , β j ) of such a curve, α j and β j take on only 2 n r β and 2 n r α different values, respectively, where r α and r β are the degrees of degenerations along the corresponding axes. The admissible points of such curves are fixed by the relations tr ( σ j β ) = 0 and tr ( σ k α ) = 0 , where σ j , σ k are some given elements of F 2 n .
The partitions of the DPS constituted by 2 n + 1 commutative curves are classified by their factorization structures (14), λ = { m 1 , m 2 , , m n } , where k m k = 2 n + 1 , which indicates the number and the lengths of the commuting sub-blocks of the stabilizers labelled by the points of a curve. Locally equivalent stabilizers can be labelled by points of different curves, but have the same factorization structure. On the other hand, curves with different factorization structures are not locally equivalent.
In particular, in the partition given by Equation (15), there are three completely factorized rays ( β = 0 , α = 0 , and β = α ), i.e., of the structure { 1 , 1 , 1 } , and the other six rays have the structure { 3 } , so that the whole partition is ( 3 , 0 , 6 ) .
The ( 0 , 9 , 0 ) partitions include only curves with the factorization { 1 , 2 } and always contain exceptional curves. One example of these partitions is given by the following seven regular and two exceptional curves:
(a)
Regular curves
α = σ 2 β + σ 3 β 2 + σ 5 β 4 , α = σ 6 β + σ 3 β 2 + σ 5 β 4 , β = σ 2 α + σ 3 α 2 + σ 5 α 4 , β = σ 6 α 2 + σ 3 α 4 , α = β + σ 6 β 2 + σ 3 β 4 , β = α + σ 3 α 2 + σ 5 α 4 , α = σ 3 β 2 + σ 5 β 4 ,
(b)
Exceptional curves
β 2 + σ 5 β = σ 2 α 2 + σ 6 α , tr ( σ 4 β ) = 0 , tr ( σ 5 α ) = 0 ; β 2 + σ 2 β = σ 6 α 2 + σ 5 α , tr ( σ 6 β ) = 0 , tr ( σ 2 α ) = 0 .

Appendix C

In this Appendix, we provide the explicit expressions for the reconstruction of the states | ψ 0 ( β ; 1 , σ 5 ) ( 0 , 9 , 0 ) and | ψ 0 ( β ; σ 2 , σ 5 ) ( 0 , 9 , 0 ) associated with the curves (46) and (50) in terms of ray-related MUB projectors with the corresponding probabilities (18):
| ψ 0 ( β ; 1 , σ 5 ) ( 0 , 9 , 0 ) ψ 0 ( β ; 1 , σ 5 ) ( 0 , 9 , 0 ) | = λ , κ F 2 3 p κ ( λ ) ( 1 , σ 5 ) 1 9 | ψ κ λ ψ κ λ | + p ˜ κ ( 1 , σ 5 ) 1 9 | ψ ˜ κ ψ ˜ κ | ,
p κ λ ( 1 , σ 5 ) = | ψ κ λ | ψ 0 ( β ; 1 , σ 5 ) ( 0 , 9 , 0 ) | 2 = λ κ 0 σ σ 2 σ 3 σ 4 σ 5 σ 6 σ 7 0 1 4 0 1 4 0 0 0 1 4 1 4 σ 1 8 1 8 1 8 1 8 1 8 1 8 1 8 1 8 σ 2 1 8 1 8 1 8 1 8 1 8 1 8 1 8 1 8 σ 3 1 4 0 0 0 1 4 1 4 0 1 4 σ 4 1 4 0 0 1 4 1 4 0 1 4 0 σ 5 1 8 1 8 1 8 1 8 1 8 1 8 1 8 1 8 σ 6 1 2 1 2 0 0 0 0 0 0 σ 7 1 4 0 1 4 1 4 0 1 4 0 0 ,
p ˜ κ ( 1 , σ 5 ) = | ψ ˜ κ | ψ 0 ( β ; 1 , σ 5 ) ( 0 , 9 , 0 ) | 2 = 1 8 ,
and
| ψ 0 ( β ; σ 2 , σ 5 ) ( 0 , 9 , 0 ) ψ 0 ( β ; σ 2 , σ 5 ) ( 0 , 9 , 0 ) | = λ κ p κ ( λ ) σ 2 , σ 5 1 9 | ψ κ λ ψ κ λ | + p ˜ κ ( σ 2 , σ 5 ) 1 9 | ψ ˜ κ ψ ˜ κ | ,
p κ λ σ 2 , σ 5 = | ψ κ λ | ψ 0 ( β ; σ 2 , σ 5 ) ( 0 , 9 , 0 ) | 2 = λ κ 0 σ σ 2 σ 3 σ 4 σ 5 σ 6 σ 7 0 1 2 1 2 0 0 0 0 0 0 σ 1 8 1 8 1 8 1 8 1 8 1 8 1 8 1 8 σ 2 1 4 0 1 4 1 4 0 1 4 0 0 σ 3 1 4 0 0 1 4 1 4 0 1 4 0 σ 4 0 1 4 1 4 1 4 0 0 1 4 0 σ 5 1 8 1 8 1 8 1 8 1 8 1 8 1 8 1 8 σ 6 0 1 4 0 1 4 1 4 1 4 0 0 σ 7 1 8 1 8 1 8 1 8 1 8 1 8 1 8 1 8 ,
p ˜ κ ( σ 2 , σ 5 ) = | ψ ˜ κ | ψ 0 ( β ; σ 2 , σ 5 ) ( 0 , 9 , 0 ) | 2 = 1 8 ,
where the MUBs associated with the partition of the DPS in rays (15) have the form:
| ψ κ λ = 2 3 α , γ F 2 3 Φ α , λ α χ α ( κ + γ ) | γ ,
| ψ ˜ κ = 2 3 / 2 γ F 2 3 χ κ γ | γ ,
and the phases Φ α , β are defined in (43).

References

  1. Wigner, E. On the quantum correction for thermodynamic equilibrium. Phys. Rev. A 1932, 40, 749–759. [Google Scholar] [CrossRef]
  2. Bloch, F. Zur Theorie des Austauschproblems und der Remanenzerscheinung der Ferromagnetika. Z. für Phys. 1932, 74, 295–335. [Google Scholar] [CrossRef]
  3. Groenewold, H.J. On the principles of elementary quantum mechanics. Physica 1949, 12, 405–460. [Google Scholar] [CrossRef]
  4. Moyal, J.E. Quantum mechanics as a statistical theory. Math. Proc. Camb. Philos. Soc. 1947, 45, 99–124. [Google Scholar] [CrossRef]
  5. Mandel, L.; Wolf, E. Coherence Properties of Optical Fields. Rev. Mod. Phys. 1965, 37, 231–287. [Google Scholar] [CrossRef]
  6. Barker, J.R.; Murray, S. A quasi-classical formulation of the Wigner function approach to quantum ballistic transport. Phys. Lett. A 1983, 93, 271–274. [Google Scholar] [CrossRef]
  7. Lin, J.; Chiu, L.C. Quantum theory of electron transport in the Wigner formalism. J. Appl. Phys. 1985, 57, 1373–1376. [Google Scholar] [CrossRef]
  8. Berry, M.V. Semi-classical mechanics in phase space: A study of Wigner’s function. Philos. Trans. R. Soc. A 1977, 287, 237–271. [Google Scholar] [CrossRef]
  9. O’Connell, R.F.; Wigner, E.P. Manifestations of Bose and Fermi statistics on the quantum distribution functionfor systems of spin-0 and spin-1/2 particles. Phys. Rev. A 1984, 30, 2613–2618. [Google Scholar] [CrossRef]
  10. Cohen, M.; Scully, M.O. Joint Wigner distribution for spin-1/2 particles. Foud. Phys. 1986, 16, 295–310. [Google Scholar] [CrossRef]
  11. Wootters, W.K. A Wigner-function formulation of finite-state quantum mechanics. Ann. Phys. 1987, 176, 1–21. [Google Scholar] [CrossRef]
  12. Gibbons, K.S.; Hoffman, M.J.; Wootters, W.K. Discrete phase space based on finite fields. Phys. Rev. A 2004, 70, 062101. [Google Scholar] [CrossRef]
  13. Gottesman, D. The Heisenberg representation of quantum computers. In Group 22: International Colloquium on Group Theoretical Methods in Physics, Proceedings of 22nd International Colloquium, Group22, ICGTMP’98, Hobart, Australia, 13–17 July 1998; Corney, S.P., Delbourgo, R., Jarvis, P.D., Eds.; International Press: Cambridge, MA, USA, 1999; pp. 32–43. [Google Scholar]
  14. Galvão, E.F. Discrete Wigner functions and quantum computational speedup. Phys. Rev. A 2005, 71, 042302. [Google Scholar] [CrossRef]
  15. Cormick, C.; Galvão, E.F.; Gottesman, D.; Paz, J.P.; Pittenger, A.O. Classicality in discrete Wigner functions. Phys. Rev. A 2006, 73, 012301. [Google Scholar] [CrossRef]
  16. Gross, D. Hudson’s theorem for finite-dimensional quantum systems. J. Math. Phys. A 2006, 47, 122107. [Google Scholar] [CrossRef]
  17. Veitch, V.C.; Ferrie, C.; Gross, D.; Emerson, J. Negative quasi-probability as a resource for quantum computation. New J. Phys. 2012, 14, 113011. [Google Scholar] [CrossRef]
  18. Raussendorf, R.; Browne, D.E.; Delfosse, N.; Okay, C.; Bermejo-Vega, J. Contextuality and Wigner-function negativity in qubit quantum computation. Phys. Rev. A 2017, 95. [Google Scholar] [CrossRef]
  19. Schmid, D.; Du, H.; Shelby, J.H.; Pusey, M.F. Uniqueness of Noncontextual Models for Stabilizer Subtheories. Phys. Rev. Lett. 2022, 129, 120403. [Google Scholar] [CrossRef]
  20. Raussendorf, R.; Okay, C.; Zurel, M.; Feldmann, P. The role of cohomology in quantum computation with magic states. Quantum 2023, 7, 979. [Google Scholar] [CrossRef]
  21. Saniga, M.; Planat, M.; Rosu, H. Mutually unbiased bases and finite projective planes. J. Opt. B Quantum Semiclass. Opt. 2004, 6, L19–L20. [Google Scholar] [CrossRef]
  22. Ivanovic, I.D. Geometrical description of quantal state determination. J. Phys. A 1984, 14, 3241–3245. [Google Scholar] [CrossRef]
  23. Wootters, W.K.; Fields, B.D. Optimal State-Determination by Mutually Unbiased Measurements. Ann. Phys. 1989, 191, 363–381. [Google Scholar] [CrossRef]
  24. Klappenecker, A.; Rötteler, M. Constructions of mutually unbiased bases. In Lecture Notes in Computer Science Vol. 2948: Finite Fields and Applications, Procceedings of 7th International Conference, Fq7, Toulouse, France, 5–9 May 2003; Mullen, G., Poli, A., Stichtenoth, H., Eds.; Springer: Berlin/Heidelberg, Germany, 2003; pp. 137–144. [Google Scholar]
  25. Bandyopadhyay, S.; Boykin, P.O.; Roychowdhury, V.; Vatan, F. A new proof for the existence of mutually unbiased bases. Algorithmica 2002, 34, 512–528. [Google Scholar] [CrossRef]
  26. Klimov, A.B.; Romero, J.L.; Björk, G.; Sánchez-Soto, L.L. Discrete phase-space structure of n-qubit mutually unbiased bases. Ann. Phys. 2009, 324, 53–72. [Google Scholar] [CrossRef]
  27. Klimov, A.B.; Romero, J.L.; Björk, G.; Sánchez-Soto, L.L. Geometrical approach to mutually unbiased bases. J. Phys. A Math. Gen. 2007, 40, 9177. [Google Scholar] [CrossRef]
  28. Pittenger, A.O.; Rubin, M.H. Wigner functions and separability for finite systems. J. Phys. A Math. Gen. 2005, 38, 6005–6036. [Google Scholar] [CrossRef]
  29. Delfosse, N.; Guerin, P.A.; Bian, J.; Raussendorf, R. Wigner Function Negativity and Contextuality in Quantum Computation on Rebits. Phys. Rev. X 2015, 5, 021003. [Google Scholar] [CrossRef]
  30. Howard, M.; Wallman, J.; Vietch, V.; Emerson, J. Contextuality supplies the ‘magic’ for quantum computation. Nature 2014, 510, 351–355. [Google Scholar] [CrossRef] [PubMed]
  31. Delfosse, N.; Okay, C.; Bermejo-Vega, J.; Browne, D.E.; Raussendorf, R. Equivalence between contextuality and negativity of the Wigner function for qudits. New J. Phys. 2017, 19, 123024. [Google Scholar] [CrossRef]
  32. Muñoz, C.; Klimov, A.B.; Sanchez-Soto, L.L. Discrete phase-space structures and Wigner functions for N qubits. Quantum Inf. Process. 2017, 16, 158. [Google Scholar] [CrossRef]
  33. Lidl, R.; Niederreiter, H. Introduction to Finite Fields and Their Applications, 2nd ed.; Cambridge University Press: Cambridge, UK, 1994; ISBN 978-052-146-094-1. [Google Scholar]
  34. Schwinger, J. Unitary Operator Bases. Proc. Natl. Acad. Sci. USA 1960, 46, 570–579. [Google Scholar] [CrossRef] [PubMed]
  35. Schwinger, J. Unitary Transformations and the Action Principle. Proc. Natl. Acad. Sci. USA 1960, 46, 883–897. [Google Scholar] [CrossRef]
  36. Lawrence, J.; Brukner, Č.; Zeilinger, A. Mutually unbiased binary observable sets on N qubits. Phys. Rev. A 2002, 65, 032320. [Google Scholar] [CrossRef]
  37. Romero, J.L.; Björk, G.; Klimov, A.B.; Sánchez-Soto, L.L. Structure of the sets of mutually unbiased bases for N qubits. Phys. Rev. A 2005, 72, 062310. [Google Scholar] [CrossRef]
  38. Durt, T. About Weyl and Wigner tomography in finite-dimensional Hilbert spaces. Open Sys. Inf. Dyn. 2006, 13, 403–413. [Google Scholar] [CrossRef]
  39. Vourdas, A. Quantum systems with finite Hilbert space. Rep. Prog. Phys. 2004, 67, 267. [Google Scholar] [CrossRef]
  40. Vourdas, A. Factorization in finite quantum systems. J. Phys. A Math. Gen. 2003, 36, 5645. [Google Scholar] [CrossRef]
  41. Paz, J.P.; Roncaglia, A.J.; Saraceno, M. Qubits in phase space: Wigner-function approach to quantum-error correction and the mean-king problem. Phys. Rev. A 2005, 72, 012309. [Google Scholar] [CrossRef]
  42. Bjork, G.; Klimov, A.B.; Sanchez-Soto, L.L. The discrete Wigner function. In Progress in Optics, 1st ed.; Wolf, E., Ed.; Elsevier: Amsterdam, The Netherlands, 2008; Volume 51, pp. 469–516. [Google Scholar] [CrossRef]
  43. Klimov, A.B.; Muñoz, C.; Romero, J.L. Geometrical approach to the discrete Wigner function in prime power dimensions. J. Phys. A Math. Gen. 2006, 39, 14471. [Google Scholar] [CrossRef]
Figure 1. Unit commutative curves along with the distribution of the respective phases (8) satisfying condition (10) for the partitions: (a), ( 3 , 0 , 6 ) described by the rays β = μ α , α = 0 , μ F 2 3 ; (b) ( 1 , 6 , 2 ) described by the curves β = μ α + α 2 + α 4 , μ F 2 3 .
Figure 1. Unit commutative curves along with the distribution of the respective phases (8) satisfying condition (10) for the partitions: (a), ( 3 , 0 , 6 ) described by the rays β = μ α , α = 0 , μ F 2 3 ; (b) ( 1 , 6 , 2 ) described by the curves β = μ α + α 2 + α 4 , μ F 2 3 .
Quantumrep 06 00005 g001
Figure 2. The Wigner function for the unit stabilizer states corresponding to (a) the Abelian curve β = α + σ 3 α 2 + σ 5 α 4 (46) and (b) the non-Abelian curve β = σ 2 α + σ 3 α 2 + σ 5 α 4 (50) for the partition of the DPS according to Equations (15) and (43). Note that the first Wigner function has a delta function form (48) corresponding to the Abelian curve, while the second function corresponding to a non-Abelian curve takes on negative values.
Figure 2. The Wigner function for the unit stabilizer states corresponding to (a) the Abelian curve β = α + σ 3 α 2 + σ 5 α 4 (46) and (b) the non-Abelian curve β = σ 2 α + σ 3 α 2 + σ 5 α 4 (50) for the partition of the DPS according to Equations (15) and (43). Note that the first Wigner function has a delta function form (48) corresponding to the Abelian curve, while the second function corresponding to a non-Abelian curve takes on negative values.
Quantumrep 06 00005 g002
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Sainz, I.; Camacho, E.; García, A.; Klimov, A.B. Tomographic Universality of the Discrete Wigner Function. Quantum Rep. 2024, 6, 58-73. https://0-doi-org.brum.beds.ac.uk/10.3390/quantum6010005

AMA Style

Sainz I, Camacho E, García A, Klimov AB. Tomographic Universality of the Discrete Wigner Function. Quantum Reports. 2024; 6(1):58-73. https://0-doi-org.brum.beds.ac.uk/10.3390/quantum6010005

Chicago/Turabian Style

Sainz, Isabel, Ernesto Camacho, Andrés García, and Andrei B. Klimov. 2024. "Tomographic Universality of the Discrete Wigner Function" Quantum Reports 6, no. 1: 58-73. https://0-doi-org.brum.beds.ac.uk/10.3390/quantum6010005

Article Metrics

Back to TopTop