Weyl Dephasing Maps in Quantum Channels

• Weyl dephasing maps are finite-dimensional quantum channels built from discrete Weyl operators that provide a basis-selective dephasing mechanism, generalizing Pauli dynamics for qudits.
• They employ a Fourier-positivity criterion and Choi–Jamiołkowski framework to rigorously ensure complete positivity and trace preservation in their operation.
• Their structure encompasses both static and dynamical aspects with anisotropic and isotropic versions, revealing intricate subgroup geometries and memory effects in convex mixtures.

Searching arXiv for papers on Weyl dephasing maps and related Weyl channels. Weyl dephasing maps are finite-dimensional quantum channels built from discrete Weyl operators and characterized by a basis-selective action on Weyl components of a density operator. In the multipartite setting, they arise as a direct generalization of Pauli-diagonal channels to tensor-product systems of qudits, while in the one-parameter dynamical setting they appear as random-unitary evolutions supported on cyclic subgroups of the discrete phase space $\mathbb Z_d\times \mathbb Z_d$. Their defining feature is that they preserve a distinguished sector of Weyl coordinates, possibly with roots-of-unity phase factors, while suppressing or averaging over complementary sectors; in this sense they form a generalized dephasing or phase-erasing family in Weyl coordinates (Basile et al., 2023, Xu et al., 22 May 2026).

1. Weyl-operator framework

For a single $d$-dimensional system, the Weyl operators are

$$U(m,n)=\sum_{k=0}^{d-1}\omega^{mk}\,|k\rangle\langle k+n|, \qquad \omega=e^{2\pi i/d},$$

with all arithmetic modulo $d$. For an $N$-particle system of $N$ qudits, the tensor-product Weyl operators are

$$U(\vec m,\vec n)=\bigotimes_{\alpha=1}^N U(m_\alpha,n_\alpha),$$

equivalently,

$$U(\vec m,\vec n)=\sum_{\vec k}\omega^{\vec m\cdot \vec k}\,|\vec k\rangle\langle \vec k+\vec n|.$$

The same construction extends to different local dimensions $d_\alpha$, with $\omega_\alpha=e^{2\pi i/d_\alpha}$ (Basile et al., 2023).

These operators form an orthogonal unitary basis of Hilbert–Schmidt space. Their algebra is governed by

$d$0

$d$1

$d$2

and

$d$3

In the single-qudit dynamical notation, one often writes $d$4 for the Weyl basis elements indexed by $d$5; the corresponding projective commutation law is

$d$6

This indexing makes explicit that Weyl dephasing maps are organized by the discrete phase space $d$7 (Xu et al., 22 May 2026).

2. Diagonal Weyl maps and complete positivity

Any density matrix on $d$8 admits the expansion

$d$9

with hermiticity condition

$$U(m,n)=\sum_{k=0}^{d-1}\omega^{mk}\,|k\rangle\langle k+n|, \qquad \omega=e^{2\pi i/d},$$0

A Weyl map is diagonal in this basis: $$U(m,n)=\sum_{k=0}^{d-1}\omega^{mk}\,|k\rangle\langle k+n|, \qquad \omega=e^{2\pi i/d},$$1 The coefficients $$U(m,n)=\sum_{k=0}^{d-1}\omega^{mk}\,|k\rangle\langle k+n|, \qquad \omega=e^{2\pi i/d},$$2 are therefore the eigenvalues of the superoperator in the Weyl basis. Trace preservation and hermiticity preservation require

$$U(m,n)=\sum_{k=0}^{d-1}\omega^{mk}\,|k\rangle\langle k+n|, \qquad \omega=e^{2\pi i/d},$$3

(Basile et al., 2023).

Complete positivity is encoded by the Choi–Jamiołkowski matrix

$$U(m,n)=\sum_{k=0}^{d-1}\omega^{mk}\,|k\rangle\langle k+n|, \qquad \omega=e^{2\pi i/d},$$4

which must satisfy $$U(m,n)=\sum_{k=0}^{d-1}\omega^{mk}\,|k\rangle\langle k+n|, \qquad \omega=e^{2\pi i/d},$$5. Because the operators $$U(m,n)=\sum_{k=0}^{d-1}\omega^{mk}\,|k\rangle\langle k+n|, \qquad \omega=e^{2\pi i/d},$$6 commute, $$U(m,n)=\sum_{k=0}^{d-1}\omega^{mk}\,|k\rangle\langle k+n|, \qquad \omega=e^{2\pi i/d},$$7 is diagonalized by a Fourier-transform basis, yielding Choi eigenvalues

$$U(m,n)=\sum_{k=0}^{d-1}\omega^{mk}\,|k\rangle\langle k+n|, \qquad \omega=e^{2\pi i/d},$$8

together with the inverse relation

$$U(m,n)=\sum_{k=0}^{d-1}\omega^{mk}\,|k\rangle\langle k+n|, \qquad \omega=e^{2\pi i/d},$$9

and normalization

$d$0

Thus valid Weyl channels are exactly those whose Fourier-transformed coefficients are nonnegative and normalized. The same condition may be written as a linear transform involving quantum Fourier matrices,

$d$1

This Fourier-positivity criterion is the basic structural constraint behind Weyl dephasing maps (Basile et al., 2023).

3. Weyl erasing channels as generalized dephasing maps

The dephasing interpretation becomes explicit in the subclass of Weyl erasing channels, defined by

$d$2

These channels “completely erase, preserve or introduce specific phases to the projections of the density matrix onto the Weyl basis.” They are the direct analogue of component-erasing maps, but in the Weyl basis rather than the Pauli basis (Basile et al., 2023).

Weyl coefficient	Action on the corresponding Weyl component
$d$3	preserved
$d$4	erased
$d$5 with phase	preserved up to phase

Operationally, these maps selectively suppress some Weyl-basis coherences while leaving a subgroup of components untouched up to phase. The paper does not use “phase-damping channel” as a formal definition, but the operational meaning is exactly a generalized dephasing or phase-erasing channel in Weyl coordinates (Basile et al., 2023).

The nonzero sector has a precise algebraic description. The set of indices with $d$6 forms a subgroup $d$7, where

$d$8

or, for different local dimensions,

$d$9

On that subgroup, $N$0 acts as a homomorphism into roots of unity,

$N$1

The coefficients outside $N$2 are set to zero. In Kraus form, the erasing channel admits a canonical interpretation in which the Kraus operators are a subset of Weyl matrices $N$3, specifically those satisfying

$N$4

with probabilities proportional to

$N$5

This identifies Weyl dephasing maps with subgroup-controlled random-unitary channels whose preserved sector is phase-twisted by a homomorphism and whose erased sector is annihilated in Weyl coordinates (Basile et al., 2023).

4. Convex geometry and subgroup classification

Because the Choi eigenvalues $N$6 are nonnegative and normalized, the set of Weyl channels forms a simplex of dimension $N$7. Its extreme points occur when exactly one Choi eigenvalue is nonzero,

$N$8

which gives

$N$9

These extreme channels are unitary conjugations because

$N$0

Hence every Weyl channel is a convex combination of unitary Weyl conjugations, i.e. a random-unitary channel. The extremality criterion can be stated equivalently as

$N$1

(Basile et al., 2023).

The group-theoretic classification of the dephasing sector depends on the arithmetic of the dimension. For prime $N$2, the index group $N$3 is a vector space over $N$4, so subgroup analysis reduces to linear algebra. For prime powers and composite dimensions, the construction proceeds through decomposition into cyclic $N$5-power groups. The classification is organized by

$N$6

together with partitions $N$7 satisfying

$N$8

Automorphisms are represented by constrained matrices $N$9, and this finite-abelian-group machinery yields an algorithmic classification of Weyl erasing channels (Basile et al., 2023).

For single-qudit dynamical maps, the subgroup structure of $$U(\vec m,\vec n)=\bigotimes_{\alpha=1}^N U(m_\alpha,n_\alpha),$$0 is classified completely via the Hermite normal form. Every subgroup $$U(\vec m,\vec n)=\bigotimes_{\alpha=1}^N U(m_\alpha,n_\alpha),$$1 has a unique HNF representation generated by

$$U(\vec m,\vec n)=\bigotimes_{\alpha=1}^N U(m_\alpha,n_\alpha),$$2

with

$$U(\vec m,\vec n)=\bigotimes_{\alpha=1}^N U(m_\alpha,n_\alpha),$$3

so that

$$U(\vec m,\vec n)=\bigotimes_{\alpha=1}^N U(m_\alpha,n_\alpha),$$4

The classification distinguishes cyclic, split rank-2, and non-split rank-2 subgroups. This subgroup geometry is the algebraic skeleton for the dynamical theory of Weyl dephasing maps (Xu et al., 22 May 2026).

5. Dynamical Weyl dephasing maps and memory effects

A Weyl dynamical map on $$U(\vec m,\vec n)=\bigotimes_{\alpha=1}^N U(m_\alpha,n_\alpha),$$5 is

$$U(\vec m,\vec n)=\bigotimes_{\alpha=1}^N U(m_\alpha,n_\alpha),$$6

where $$U(\vec m,\vec n)=\bigotimes_{\alpha=1}^N U(m_\alpha,n_\alpha),$$7 is a probability distribution. Each Weyl operator is an eigenoperator,

$$U(\vec m,\vec n)=\bigotimes_{\alpha=1}^N U(m_\alpha,n_\alpha),$$8

with eigenvalues

$$U(\vec m,\vec n)=\bigotimes_{\alpha=1}^N U(m_\alpha,n_\alpha),$$9

The dynamics is therefore completely encoded by the time-dependent weights $$U(\vec m,\vec n)=\sum_{\vec k}\omega^{\vec m\cdot \vec k}\,|\vec k\rangle\langle \vec k+\vec n|.$$0, or equivalently by the spectral data $$U(\vec m,\vec n)=\sum_{\vec k}\omega^{\vec m\cdot \vec k}\,|\vec k\rangle\langle \vec k+\vec n|.$$1 (Xu et al., 22 May 2026).

The time-local description uses

$$U(\vec m,\vec n)=\sum_{\vec k}\omega^{\vec m\cdot \vec k}\,|\vec k\rangle\langle \vec k+\vec n|.$$2

with GKLS form

$$U(\vec m,\vec n)=\sum_{\vec k}\omega^{\vec m\cdot \vec k}\,|\vec k\rangle\langle \vec k+\vec n|.$$3

For Weyl maps one sets $$U(\vec m,\vec n)=\sum_{\vec k}\omega^{\vec m\cdot \vec k}\,|\vec k\rangle\langle \vec k+\vec n|.$$4 and typically $$U(\vec m,\vec n)=\sum_{\vec k}\omega^{\vec m\cdot \vec k}\,|\vec k\rangle\langle \vec k+\vec n|.$$5, obtaining

$$U(\vec m,\vec n)=\sum_{\vec k}\omega^{\vec m\cdot \vec k}\,|\vec k\rangle\langle \vec k+\vec n|.$$6

The RHP criterion yields the standard dynamical trichotomy: CP-divisible or Markovian evolution when all $$U(\vec m,\vec n)=\sum_{\vec k}\omega^{\vec m\cdot \vec k}\,|\vec k\rangle\langle \vec k+\vec n|.$$7, semigroup evolution when the rates are constant and nonnegative, non-Markovianity when some rate becomes negative, and eternal non-Markovianity (ENM) when at least one rate satisfies

$$U(\vec m,\vec n)=\sum_{\vec k}\omega^{\vec m\cdot \vec k}\,|\vec k\rangle\langle \vec k+\vec n|.$$8

(Xu et al., 22 May 2026).

The isotropic/anisotropic distinction is central. An anisotropic Weyl map has the form

$$U(\vec m,\vec n)=\sum_{\vec k}\omega^{\vec m\cdot \vec k}\,|\vec k\rangle\langle \vec k+\vec n|.$$9

with fixed nonuniform weights $d_\alpha$0, whereas an isotropic Weyl map is

$d_\alpha$1

A key structural result is that anisotropic Weyl maps with nonuniform weights and at least two distinct nontrivial eigenvalues cannot form a semigroup. By contrast, isotropic maps have only two spectral values,

$d_\alpha$2

and they generate a Markovian semigroup precisely when

$d_\alpha$3

Equivalently,

$d_\alpha$4

with zero rates otherwise (Xu et al., 22 May 2026).

For $d_\alpha$5, the isotropic map reduces to the Weyl dephasing map

$d_\alpha$6

Let

$d_\alpha$7

for $d_\alpha$8, so that the cyclic subgroup generated by $d_\alpha$9 has order $\omega_\alpha=e^{2\pi i/d_\alpha}$0. If

$\omega_\alpha=e^{2\pi i/d_\alpha}$1

then the map is eternally non-Markovian if either $\omega_\alpha=e^{2\pi i/d_\alpha}$2 is odd, or $\omega_\alpha=e^{2\pi i/d_\alpha}$3 is even and $\omega_\alpha=e^{2\pi i/d_\alpha}$4. In higher-dimensional Weyl systems, therefore, a single dephasing map can already be irreducibly eternally non-Markovian (Xu et al., 22 May 2026).

6. Convex mixtures, qutrit structure, and relation to Pauli maps

Convexity plays a nontrivial role in Weyl dephasing dynamics. The isotropic Weyl map can be written as an equal-weight convex mixture of dephasing maps,

$\omega_\alpha=e^{2\pi i/d_\alpha}$5

A striking consequence is that a convex combination of eternally non-Markovian Weyl dephasing maps can nevertheless generate a Markovian semigroup. More precisely, when the constituent dephasing maps are ENM with odd $\omega_\alpha=e^{2\pi i/d_\alpha}$6, the equal-weight isotropic mixture can still yield a semigroup provided

$\omega_\alpha=e^{2\pi i/d_\alpha}$7

Thus non-Markovianity is not additive under mixing (Xu et al., 22 May 2026).

The converse direction is also established for mixtures of $\omega_\alpha=e^{2\pi i/d_\alpha}$8 distinct isotropic Weyl semigroups,

$\omega_\alpha=e^{2\pi i/d_\alpha}$9

If all semigroups have the same subgroup order

$d$00

then the mixture is eternally non-Markovian whenever

$d$01

where

$d$02

The underlying mechanism is geometric: if the union of subgroup supports leaves at least one nontrivial phase-space point uncovered, at least one decay rate remains negative for all $d$03 (Xu et al., 22 May 2026).

The qutrit case $d$04 makes these statements explicit. Here $d$05 has $d$06 elements, there are $d$07 subgroups of order $d$08, and distinct order-$d$09 subgroups intersect only at $d$10. For $d$11 or $d$12, the union of chosen subgroups does not cover the entire phase space, so the mixtures are eternally non-Markovian. For $d$13, the union can cover the whole phase space: uniform mixing $d$14 yields a Markovian map, whereas nonuniform mixing can produce non-Markovian or even eternally non-Markovian dynamics (Xu et al., 22 May 2026).

Weyl dephasing maps reduce to familiar Pauli structures in the binary case. For $d$15, Weyl operators reduce to Pauli operators, and the extreme Weyl channels correspond to the vertices of the qubit Pauli-channel tetrahedron (Basile et al., 2023). The Weyl framework is nevertheless strictly broader: generalized Pauli maps can be embedded into it, and for $d$16 they may be rewritten using Weyl operators with grouped weights over commuting subsets (Xu et al., 22 May 2026). This places Weyl dephasing maps as the higher-dimensional extension of Pauli dephasing channels, with additional subgroup geometry and convex phenomena that do not appear in the qubit setting.