Mutually Orthogonal Quantum Latin Squares

• Mutually orthogonal quantum Latin squares (MOQLS) are quantum analogues of classical Latin squares, replacing symbols with unit vectors in Hilbert space to form orthonormal bases.
• They facilitate novel constructions in quantum error correction, quantum codes, and entanglement by ensuring pairwise orthogonality via tensor-product conditions.
• The theory mirrors classical MOLS with stringent existence bounds, yet introduces critical distinctions such as nonexistence in the product-state framework for order six.

Searching arXiv for recent and foundational papers on mutually orthogonal quantum Latin squares. Mutually orthogonal quantum Latin squares are quantum-combinatorial analogues of classical mutually orthogonal Latin squares in which symbols are replaced by vectors in Hilbert space and orthogonality is imposed through tensor-product basis conditions rather than symbol-pair uniqueness. A quantum Latin square of order $n$ is an $n\times n$ array of unit vectors in $\mathbb C^n$ such that every row and every column forms an orthonormal basis of $\mathbb C^n$; two such squares are orthogonal when the $n^2$ vectors obtained by tensoring corresponding entries form an orthonormal basis of $\mathbb C^n\otimes\mathbb C^n$ (Musto et al., 2018). In this sense, mutually orthogonal quantum Latin squares (MOQLS) extend the classical theory of MOLS while retaining a stringent product-state geometry. Their study connects combinatorial design theory, unitary error bases, mutually unbiased maximally entangled bases, quantum orthogonal arrays, and multipartite entanglement (Musto et al., 2018, Zang et al., 2021).

1. Definition and basic formalism

A quantum Latin square (QLS) of dimension $n$ is an $n\times n$ array

$$\bigl\{\,|\Psi_{ij}\rangle \in \mathbb C^n \mid i,j\in [n]\,\bigr\},$$

with $[n]=\{0,1,\dots,n-1\}$, such that each row and each column is an orthonormal basis of $n\times n$0 (Musto et al., 2018). Equivalently, for each fixed row index $n\times n$1, the set $n\times n$2 is an orthonormal basis, and for each fixed column index $n\times n$3, the set $n\times n$4 is an orthonormal basis. Classical Latin squares embed into this framework by replacing each symbol with a computational-basis vector (Musto et al., 2018).

For two QLS $n\times n$5 and $n\times n$6 of dimension $n\times n$7, the standard orthogonality notion is the tensor-basis criterion

$n\times n$8

forming an orthonormal basis of $n\times n$9 (Musto et al., 2018). Equivalently,

$\mathbb C^n$0

(Musto et al., 2018). A family of $\mathbb C^n$1 QLS is mutually orthogonal if it is pairwise orthogonal in this sense (Musto et al., 2018).

The literature also contains a distinct notion based on bipartite or multipartite array entries, where a “pair of orthogonal quantum Latin squares” is itself an array of bipartite states and row and column sums satisfy 1-uniformity conditions (Goyeneche et al., 2017). This entangled notion is not identical to the product-state MOQLS formalism above. The distinction is mathematically consequential, especially at order $\mathbb C^n$2, where the entangled setting behaves differently from the non-entangled one (Ball et al., 2 Mar 2026).

2. Orthogonality notions and their unification

A central development in the theory was the simplification of orthogonality for QLS. Earlier work of Goyeneche, Raissi, Di Martino, and Życzkowski used a larger tensor expression with partial traces. A later formulation proved that the essential content reduces exactly to the tensor-basis condition above, and that the extra partial-trace conditions are redundant (Musto et al., 2018). In particular, the simplified criterion is equivalent to the earlier GRMZ notion, but only one trace condition carries substantive information (Musto et al., 2018).

The same simplification extends from pairs to families. Pairwise orthogonality is sufficient; no genuinely higher-arity trace condition is needed for MOQLS families (Musto et al., 2018). This aligns the subject with classical MOLS, where mutual orthogonality is likewise pairwise.

A separate line of work developed “weak orthogonality” for QLS in connection with mutually unbiased maximally entangled bases. There, for QLS $\mathbb C^n$3, weak orthogonality requires that for all $\mathbb C^n$4, there exists a unique $\mathbb C^n$5 such that

$\mathbb C^n$6

(Musto, 2016). In the classical case this is equivalent to left orthogonality of Latin squares, and hence to ordinary orthogonality up to left conjugation (Musto, 2016). This notion underlies a different MUB-oriented construction and should not be conflated with the product-state MOQLS definition used in the modern QLS literature (Musto, 2016, Musto et al., 2018).

The distinction between product-state orthogonality and entangled orthogonality is equally important. In the entangled framework of quantum combinatorial designs, entries may be bipartite or multipartite states that do not factor into separate QLS, and orthogonality is encoded by global basis orthogonality together with 1-uniformity of line sums (Goyeneche et al., 2017). A plausible implication is that several statements in the literature that use the phrase “orthogonal quantum Latin squares” are only comparable after fixing which orthogonality notion is intended.

3. Structural results and bounds

The quantum theory preserves a key extremal feature of classical MOLS: any family of MOQLS of dimension $\mathbb C^n$7 has size at most $\mathbb C^n$8 (Musto et al., 2018). The proof normalizes first rows via equivalence transformations and then extracts $\mathbb C^n$9 linearly independent vectors in $\mathbb C^n$0, forcing $\mathbb C^n$1 (Musto et al., 2018). Thus quantization does not enlarge the maximal family size beyond the classical bound.

This upper bound has stronger consequences in certain regimes. It is known that if equality holds, the family is classical (Ball et al., 2 Mar 2026). The order-six nonexistence work further notes that if there exist $\mathbb C^n$2 MOQLS$\mathbb C^n$3, then there exist $\mathbb C^n$4 MOLS$\mathbb C^n$5, hence $\mathbb C^n$6 MOLS$\mathbb C^n$7, and therefore a projective plane of order $\mathbb C^n$8 (Ball et al., 2 Mar 2026). This suggests that near-maximal MOQLS families are forced back into classical incidence geometry.

Equivalence transformations also play a structural role. Orthogonality is preserved under common row and column permutations together with individual unitaries and phase multipliers on each QLS (Musto et al., 2018). In the non-classical existence theory, two QLS $\mathbb C^n$9 and $n^2$0 are equivalent if there exist a unitary $n^2$1, phases $n^2$2, and permutations $n^2$3 such that

$n^2$4

for all $n^2$5 (Han et al., 27 Jul 2025). A QLS equivalent to one arising from a classical Latin square is called classical; otherwise it is non-classical (Han et al., 27 Jul 2025).

A useful diagnostic for non-classicality is that classical QLS have only inner products of modulus $n^2$6 or $n^2$7 among entries. Hence any QLS exhibiting an entry overlap of modulus strictly between $n^2$8 and $n^2$9 is genuinely quantum (Zang et al., 2021). This criterion is used repeatedly in explicit constructions of non-classical orthogonal pairs (Zang et al., 2021).

4. Existence theory and explicit constructions

The first explicit pair of orthogonal QLS not equivalent to any pair of orthogonal classical Latin squares was constructed in dimension $\mathbb C^n\otimes\mathbb C^n$0 (Musto et al., 2018). That work also proved that the example is genuinely nonclassical by showing that one square cannot be transformed into a computational-basis Latin square via a single global unitary, phases, and permutations (Musto et al., 2018). The same paper’s constructive novelty is at the level of pairs rather than larger nonclassical MOQLS families (Musto et al., 2018).

A more systematic existence theory was later developed for non-classical MOQLS using combinatorial design methods (Han et al., 27 Jul 2025). That work introduced idempotent QLS, self-orthogonal QLS, holey quantum Latin squares, and orthogonality on these variants, then applied PBD constructions and filling-in-holes constructions to obtain broad asymptotic existence theorems (Han et al., 27 Jul 2025). The principal results include:

Structure	Existence claim	Exception range stated
non-classical $\mathbb C^n\otimes\mathbb C^n$1-idempotent MOQLS$\mathbb C^n\otimes\mathbb C^n$2	exists for $\mathbb C^n\otimes\mathbb C^n$3	except possibly $\mathbb C^n\otimes\mathbb C^n$4
non-classical $\mathbb C^n\otimes\mathbb C^n$5-MOQLS$\mathbb C^n\otimes\mathbb C^n$6	exists for $\mathbb C^n\otimes\mathbb C^n$7	except possibly $\mathbb C^n\otimes\mathbb C^n$8
non-classical $\mathbb C^n\otimes\mathbb C^n$9-MOQLS$n$0	exists for $n$1	except possibly $n$2
non-classical SOQLS$n$3	exists for $n$4	none stated beyond threshold

These results place MOQLS existence into the recursive-combinatorial tradition of classical MOLS, but with local non-computational bases injected into chosen subspaces to force non-classicality (Han et al., 27 Jul 2025).

An earlier construction paper already provided direct product methods for MOQLS and lower bounds for the maximal number $n$5 of mutually orthogonal non-classical QLS of dimension $n$6 (Zang et al., 2021). If there exists a $n$7-MOQLS$n$8 and a $n$9-MOQLS$n\times n$0, then there exists a $n\times n$1-MOQLS$n\times n$2 by tensoring corresponding entries (Zang et al., 2021). More generally, if there exists a $n\times n$3-MOQLS$n\times n$4 for each $n\times n$5, then there exists a $n\times n$6-MOQLS$n\times n$7 with $n\times n$8 (Zang et al., 2021).

The same work transferred classical MOLS existence into quantum existence. If $n\times n$9, $$\bigl\{\,|\Psi_{ij}\rangle \in \mathbb C^n \mid i,j\in [n]\,\bigr\},$$0, and $$\bigl\{\,|\Psi_{ij}\rangle \in \mathbb C^n \mid i,j\in [n]\,\bigr\},$$1 for all $$\bigl\{\,|\Psi_{ij}\rangle \in \mathbb C^n \mid i,j\in [n]\,\bigr\},$$2, then there exists $$\bigl\{\,|\Psi_{ij}\rangle \in \mathbb C^n \mid i,j\in [n]\,\bigr\},$$3-MOQLS$$\bigl\{\,|\Psi_{ij}\rangle \in \mathbb C^n \mid i,j\in [n]\,\bigr\},$$4 with $$\bigl\{\,|\Psi_{ij}\rangle \in \mathbb C^n \mid i,j\in [n]\,\bigr\},$$5 (Zang et al., 2021). It also proved lower bounds such as

$$\bigl\{\,|\Psi_{ij}\rangle \in \mathbb C^n \mid i,j\in [n]\,\bigr\},$$6

when $$\bigl\{\,|\Psi_{ij}\rangle \in \mathbb C^n \mid i,j\in [n]\,\bigr\},$$7 with $$\bigl\{\,|\Psi_{ij}\rangle \in \mathbb C^n \mid i,j\in [n]\,\bigr\},$$8, together with broad existence statements outside explicit exceptional sets $$\bigl\{\,|\Psi_{ij}\rangle \in \mathbb C^n \mid i,j\in [n]\,\bigr\},$$9 (Zang et al., 2021).

5. Incomplete, generalized, and entangled variants

Incomplete quantum Latin squares (IQLS) were introduced as quantum analogues of incomplete Latin squares with holes indexed by mutually orthogonal subspaces $[n]=\{0,1,\dots,n-1\}$0 (Zang et al., 2021). Orthogonality for two IQLS is defined by requiring their superimposed tensor states to form an orthonormal basis of

$[n]=\{0,1,\dots,n-1\}$1

(Zang et al., 2021). These incomplete objects support hole-filling constructions that produce full QLS, MOQLS, and self-orthogonal QLS (Zang et al., 2021, Han et al., 27 Jul 2025).

Generalized mutually orthogonal quantum Latin squares (GMOQLS) enlarge the framework further by allowing each cell to contain a $[n]=\{0,1,\dots,n-1\}$2-partite pure state

$[n]=\{0,1,\dots,n-1\}$3

rather than a product of $[n]=\{0,1,\dots,n-1\}$4 local vectors (Zang et al., 2021). The defining conditions require orthogonality of all cells, row and column quantum-Latin conditions after tracing out $[n]=\{0,1,\dots,n-1\}$5 parties, and a global support condition after tracing out $[n]=\{0,1,\dots,n-1\}$6 parties (Zang et al., 2021). If every entry is fully separable, GMOQLS reduce to ordinary MOQLS (Zang et al., 2021).

This generalized theory is equivalent to quantum orthogonal arrays of size $[n]=\{0,1,\dots,n-1\}$7. Specifically,

$[n]=\{0,1,\dots,n-1\}$8

(Zang et al., 2021). The corresponding multipartite state

$[n]=\{0,1,\dots,n-1\}$9

is then 2-uniform (Zang et al., 2021). This places MOQLS within a broader entanglement-theoretic architecture in which ordinary MOQLS appear as the fully separable corner of a larger family of quantum combinatorial designs.

A different but related entangled framework from quantum combinatorial designs gives, for every $n\times n$00, a triple of entangled MOQLS of size $n\times n$01 derived from generalized Bell states

$n\times n$02

and equivalent to a $n\times n$03 (Goyeneche et al., 2017). In that setting the resulting designs are “essentially quantum,” with entangled entries that need not decompose into separate QLS (Goyeneche et al., 2017). This suggests that the entangled and non-entangled theories should be read as parallel, not interchangeable, generalizations.

6. Order six, Euler’s problem, and current interpretation

Order $n\times n$04 is the critical case in the subject because it sharpens the distinction between non-entangled MOQLS and entangled quantum Latin structures. A recent theorem proves: $n\times n$05 (Ball et al., 2 Mar 2026). This is the product-state quantum analogue of the classical nonexistence of orthogonal Latin squares of order $n\times n$06, and it establishes that “quantization without entanglement does not defeat Euler’s obstruction” (Ball et al., 2 Mar 2026).

The proof combines pattern analysis, local-unitary normalization, graph-theoretic reformulation via orthonormal representations of complements of Latin square graphs, classification of Latin squares of order six up to paratopy, and targeted computer-assisted elimination (Ball et al., 2 Mar 2026). A key reduction shows that if two MOQLS$n\times n$07 existed, then one could assume one of them classical (Ball et al., 2 Mar 2026). This reduction allows the remaining problem to be phrased in terms of orthonormal representations of complements of Latin square graphs in $n\times n$08 (Ball et al., 2 Mar 2026).

The same paper emphasizes the sharp contrast with entangled quantum Latin squares of order $n\times n$09, which do exist and are equivalent to an $n\times n$10 state (Ball et al., 2 Mar 2026). Thus the slogan that “the thirty-six quantum officers, if they exist, must be entangled” is literal in this context (Ball et al., 2 Mar 2026).

Single QLS of order $n\times n$11 nevertheless remain relevant as potential test cases. Two explicit QLS of order $n\times n$12 with cardinalities $n\times n$13 and $n\times n$14 were constructed in 2026 (Xu, 15 May 2026). The first uses a direct-sum decomposition

$n\times n$15

while the second is built from two-dimensional Hadamard pairs supported on coordinate planes (Xu, 15 May 2026). The paper does not prove any orthogonality theorem, but it exposes structural features—support patterns, repeated vectors up to phase, and direct-sum or coordinate-plane geometry—that are directly relevant to MOQLS searches (Xu, 15 May 2026).

For the cardinality-$n\times n$16 square, the repeated-use pattern is severe: only $n\times n$17 phase-classes occur across $n\times n$18 cells (Xu, 15 May 2026). This implies that any orthogonal partner would face strong packing constraints, because repeated vectors in one square force orthogonality relations among the matching entries of the partner (Xu, 15 May 2026). For the cardinality-$n\times n$19 square, the geometry is local to coordinate $n\times n$20-planes, suggesting a different search strategy based on support partitions and local Hadamard blocks (Xu, 15 May 2026). A plausible implication is that explicit low-cardinality QLS can serve more as obstruction witnesses and structured ansätze than as direct evidence toward existence of MOQLS at order $n\times n$21.

7. Applications and related frameworks

Orthogonal QLS and MOQLS have several direct applications in quantum information. Orthogonal QLS yield quantum codes, and the same principle extends to quantum Latin isometry squares, from which one obtains error-detecting codes and a characterization of unitary error bases (Musto et al., 2018). In that generalized setting, orthogonality of isometry squares leads to an encoding tensor

$n\times n$22

that detects a single error (Musto et al., 2018).

Weakly orthogonal QLS also support the construction of mutually unbiased maximally entangled bases in square dimension (Musto, 2016). Given two families of Hadamards $n\times n$23 and $n\times n$24, and a pair of weak orthogonal QLS $n\times n$25, the bases $n\times n$26 and $n\times n$27 are mutually unbiased (Musto, 2016). The construction strictly generalizes the Beth–Wocjan construction because it allows quantum Latin squares rather than classical Latin squares and allows indexed families of Hadamards rather than a single fixed Hadamard (Musto, 2016).

The broader quantum-combinatorial design framework connects MOQLS, quantum orthogonal arrays, and $n\times n$28-uniform states. In particular, a family of $n\times n$29 MOQLH gives a $n\times n$30-uniform state

$n\times n$31

on $n\times n$32 parties (Goyeneche et al., 2017). For the square case $n\times n$33, this relationship underlies the role of MOQLS and GMOQLS in constructing 2-uniform and absolutely maximally entangled states (Goyeneche et al., 2017, Zang et al., 2021).

A structural perspective from quantum magic squares is also relevant. Quantum Latin squares correspond to rank-one quantum magic squares, and the semiclassical quantum Latin squares are precisely those constructed from classical Latin squares (Cuevas et al., 2022). This clarifies that the classical-to-quantum embedding is exact at the semiclassical level, while the full set of quantum Latin squares is strictly larger (Cuevas et al., 2022). This suggests that any genuinely quantum MOQLS theory must extend beyond the semiclassical subclass.

8. Conceptual landscape and common misunderstandings

A common misconception is that any “quantum solution” of a Latin-square orthogonality problem automatically belongs to the MOQLS framework. Order $n\times n$34 disproves this: entangled quantum Latin squares of order $n\times n$35 exist, but non-entangled MOQLS of order $n\times n$36 do not (Ball et al., 2 Mar 2026). The resource that changes the answer is entanglement, not merely superposition.

Another misconception is that orthogonality of quantum Latin squares is intrinsically multipartite or requires higher-order trace constraints. In the non-entangled QLS framework, orthogonality is simply the pairwise tensor-basis condition, and the more complicated GRMZ-style family conditions collapse to ordinary pairwise orthogonality (Musto et al., 2018).

It is also inaccurate to treat all uses of “orthogonal quantum Latin squares” as equivalent across the literature. At least three notions coexist: the product-state MOQLS notion used in modern QLS theory (Musto et al., 2018), the weak-orthogonality notion developed for MUB constructions (Musto, 2016), and the entangled orthogonality notion for quantum Latin arrangements and quantum orthogonal arrays (Goyeneche et al., 2017). These notions agree only in restricted classical or separable settings.

Finally, the existence theory for non-classical MOQLS should not be read as implying arbitrary abundance in every order. There are strong positive results outside explicit exceptional sets (Han et al., 27 Jul 2025, Zang et al., 2021), but order $n\times n$37 is now known to be a sharp nonexistence point for non-entangled pairs (Ball et al., 2 Mar 2026). The subject therefore exhibits both asymptotic abundance and low-dimensional rigidity.

9. Outlook

The modern theory of mutually orthogonal quantum Latin squares has settled several foundational issues. The product-state orthogonality notion is now conceptually minimal and equivalent to earlier cumbersome formulations (Musto et al., 2018). Explicit non-classical orthogonal pairs exist (Musto et al., 2018). Broad recursive existence results for non-classical $n\times n$38- and $n\times n$39-MOQLS are available through PBD and filling-in-holes constructions (Han et al., 27 Jul 2025). The generalized entangled theory is tied precisely to quantum orthogonal arrays and 2-uniform states (Zang et al., 2021), and order $n\times n$40 has been sharply separated into impossible non-entangled MOQLS versus possible entangled analogues (Ball et al., 2 Mar 2026).

Open structure remains. The order-$n\times n$41 question “Are 2 MOQLS(7) classical?” is explicitly identified as unresolved (Ball et al., 2 Mar 2026). More broadly, the theory still lacks a full classification of genuinely quantum MOQLS in low orders, and explicit constructions beyond pairs remain comparatively sparse in the non-entangled framework (Musto et al., 2018). The recent appearance of highly structured QLS of order $n\times n$42 with controlled cardinalities and transparent geometric support patterns suggests that future progress may rely on combining combinatorial recursion with fine-grained Hilbert-space geometry, especially support partitions, Hadamard rotations, and low-cardinality obstruction analysis (Xu, 15 May 2026).