---
title: Wave Function Cloning Algorithms
url: https://www.emergentmind.com/topics/wave-function-cloning-algorithm
type: topic
---

# Wave Function Cloning Algorithms

“Wave function cloning algorithm” is not a single standardized algorithmic primitive in quantum information theory. In the arXiv literature, the term refers to several distinct but related constructions: hybrid variational procedures that learn an approximate cloner for a restricted ensemble of states, symmetry-assisted exact cloning on a specially structured two-state family, physical phase-covariant optical cloners and the numerical methods used to model them, and broader channel-cloning frameworks in which state cloning appears as a special case. Across these usages, the common feature is that the objective is never universal exact copying of an arbitrary unknown quantum state; rather, each construction either respects the no-cloning theorem by restricting the input family or by allowing only approximate cloning, or else recasts cloning as a computational or higher-order process problem [2407.06026].

## 1. Scope of the term in the literature

The literature uses “wave function cloning” for heterogeneous tasks that share only a family resemblance. In one line of work, a parameterized quantum circuit or programmable interferometer is trained to maximize cloning fidelity for a chosen state ensemble. In another, a structured memory state is “perfectly cloned up to global mirroring,” which is operationally sufficient only because the downstream associative-retrieval task identifies mirrored and original memories up to a final inversion. In yet another, phase-covariant optical amplification is treated as a physical cloning transformation for equatorial polarization qubits, while the algorithmic novelty lies in numerically evaluating the resulting multiphoton states. A more recent abstraction embeds state cloning into deterministic cloning of quantum channels via trash-and-replace channels, whereas complexity-theoretic work addresses the opposite question: whether cloning may remain computationally intractable even when a verifier description is known [2012.11424].

| Usage | Object cloned | Characterization |
|---|---|---|
| Variational quantum cloning | Restricted state family | Hybrid quantum-classical approximate cloner |
| Mirror modular cloning | \(|M\rangle\) or \(|\bar M\rangle\) | Exact only on a two-state, symmetry-related family |
| Phase-covariant optical cloning | Equatorial polarization qubits | Physical restricted cloner, not universal copying |
| Channel cloning reformulation | Trash-and-replace channels \(\mathcal T_\rho\) | State cloning as a special case of channel cloning |
| Witness-cloning hardness | Accepted states of a verifier | Conditional impossibility framework, not a constructive cloner |

This multiplicity of meanings matters because apparently similar claims can refer to very different operational tasks. “Optimal,” “perfect,” and even “cloning” are task-relative terms in this area: optimality is usually fidelity-optimality for a specified ensemble; perfection may hold only up to a symmetry; and a cloning algorithm may mean a trainable search procedure, a fixed physical transformation, or an optimization framework rather than a closed-form circuit.

## 2. Information-theoretic constraints and fidelity benchmarks

The foundational constraint is the no-cloning theorem. In the formulation recalled in experimental variational cloning, the forbidden ideal process would require a unitary satisfying
\[
U\ket{\psi}_1 \ket{0}_2\ket{\mathcal{M}}_\mathcal{M} \rightarrow \ket{\psi}_1\ket{\psi}_2\ket{\mathcal{M}'}_\mathcal{M}
\quad \forall \ket{\psi},
\]
which cannot exist for all unknown \(\ket{\psi}\). A generic qubit is written as
\[
\ket{\psi} = \cos(\theta)\ket{0} + \sin(\theta)e^{i\phi}\ket{1}.
\]
Accordingly, legitimate cloning problems replace universal exact copying by approximate deterministic cloning or by exact cloning on a restricted state set [2407.06026].

The canonical universal symmetric \(1\rightarrow 2\) qubit benchmark is the Bužek–Hillery cloner, for which the reduced clone states satisfy
\[
\rho_1=\rho_2=\frac{5}{6}\dyad{\psi}+\frac{1}{6}\dyad{\psi_\bot},
\]
with fidelity
\[
F_{1,2}= \Tr\!\left(\sqrt{\sqrt{\sigma}\rho_{1,2}\sqrt{\sigma}}\right)^2 = \frac{5}{6},
\qquad \sigma=\dyad{\psi}.
\]
The value \(5/6\) is the optimal universal symmetric \(1\to 2\) qubit fidelity benchmark. Restricted ensembles admit higher fidelities. For symmetric phase-covariant \(1\to 2\) cloning on the equator of the Bloch sphere,
\[
\ket{\psi_\phi}=\frac{\ket{0}+e^{i\phi}\ket{1}}{\sqrt{2}},
\]
the quoted optimal fidelity is
\[
F_O = 0.853,
\]
to be compared with the semiclassical measure-and-prepare benchmark
\[
F_{SC}=0.750.
\]
For state-dependent cloning of only two known nonorthogonal states, the achievable fidelity can exceed the phase-covariant value because the ensemble is smaller [2407.06026].

Variational quantum cloning work formalizes several figures of merit. For an \(M\to N\) task with output state \(\rho_{\boldsymbol\theta}\), the reduced state of clone \(j\) is \(\rho^j_{\boldsymbol\theta}=\operatorname{Tr}_{\bar j}(\rho_{\boldsymbol\theta})\). Local clone fidelity is
\[
F^j_{\mathsf L}(\boldsymbol\theta):=
F\!\left(|\psi\rangle\langle\psi|,\rho^j_{\boldsymbol\theta}\right),
\]
while global fidelity is
\[
F_{\mathsf G}(\boldsymbol\theta):=
F\!\left(|\psi\rangle\langle\psi|^{\otimes N},\rho_{\boldsymbol\theta}\right).
\]
The paper distinguishes local, squared-local, global, and asymmetric objectives, and proves the general inequality
\[
\mathcal C_{\mathsf L}(\boldsymbol\theta)\le \mathcal C_{\mathsf G}(\boldsymbol\theta)\le N\,\mathcal C_{\mathsf L}(\boldsymbol\theta),
\]
while stressing that this does not imply global and local optima always coincide [2012.11424].

At the asymptotic level, the channel-cloning formulation gives a sharp restriction for deterministic state cloning of any continuous set of states: the replication rate is at most linear. In that framework, deterministic super-linear replication of continuous state families obeys a strong converse; if one attempts \(M=N^{1+\delta}\) with \(\delta>0\), then asymptotically
\[
\mathsf D^{\mathcal S}\ge \frac{\pi}{4},
\qquad
F^{\mathcal S}\le \frac{1}{\sqrt 2}.
\]
This situates approximate wave-function cloning within a stricter asymptotic landscape than the finite-\(N\) fidelity benchmarks usually emphasized in circuit-design work [2509.08059].

## 3. Variational quantum cloning as a hybrid learning algorithm

Variational quantum cloning (VQC) is a hybrid quantum-classical algorithm that learns a short-depth parameterized circuit implementing an approximate cloner for a specified family of states. The central object is not a universal exact cloner, but a family-specific circuit \(U_{\boldsymbol\theta,\boldsymbol g}\) acting on \(M\) input copies, \(N-M\) blanks, and optional ancilla. The optimization problem is written as
\[
(\boldsymbol\theta^*,\boldsymbol g^*)=
\arg\min_{\boldsymbol\theta,\ \boldsymbol g\in\mathcal G}
\mathcal C(\boldsymbol\theta,\boldsymbol g),
\]
where \(\boldsymbol\theta\) are continuous parameters and \(\boldsymbol g\) is a discrete circuit structure chosen from a gate pool. The algorithm samples training states from a target family \(\mathcal S\), estimates local or global fidelities on quantum hardware, and updates both parameters and, optionally, circuit structure on a classical processor. The continuous optimization uses analytic gradients via the parameter-shift rule, while the structure-learning component perturbs a gate sequence and reoptimizes it. The paper further gives a sample-complexity estimate
\[
\mathrm{Samp}=L\times K=
\mathcal O\!\left(\gamma^{-2}\log(2/\delta)\right),
\]
and argues that sufficiently shallow alternating layered ansätze avoid barren plateaus for the local cloning cost [2012.11424].

This framework is explicitly motivated by quantum cryptanalysis. The 2020 VQC study treats approximate cloning as a resource for attacks on quantum coin-flipping protocols and emphasizes hardware-efficient short-depth circuits rather than analytically optimal but hardware-unfriendly constructions. For phase-covariant \(1\to 2\) cloning, it recalls the optimal local fidelity
\[
F_{\mathsf L}^{\mathrm{opt}}=
\frac12\left(1+\frac1{\sqrt2}\right)\approx 0.853,
\]
and the optimal global fidelity
\[
F_{\mathsf G}^{\mathrm{opt}}=
\frac18(1+\sqrt2)^2\approx 0.72.
\]
A learned circuit achieves nearly optimal local fidelity,
\[
F_B^{(a)}\approx F_E^{(a)}\approx 0.85,
\]
but with lower global fidelity,
\[
F_{\mathsf G}^{(a)}\approx 0.638.
\]
On Rigetti Aspen-8, the learned hardware-adapted circuit outperforms the textbook ideal circuit by about \(15\%\) in observed cloning fidelity because it uses fewer entangling gates and fits the device topology better. The same study reports application-specific attack biases, including \(\epsilon\approx 0.29\) for a VQC-based attack on \(\mathcal P_1\), and \(\epsilon^{\mathrm I}_{\mathcal P_2,\mathrm{VQC}}\approx 0.345\), \(\epsilon^{\mathrm{II}}_{\mathcal P_2,\mathrm{VQC}}=0.241\) for two attack models on \(\mathcal P_2\) [2012.11424].

A 2024 photonic implementation realizes the same general idea in integrated optics. There, two dual-rail photonic qubits are used: one carries the state to be cloned, and the other is a blank target register initialized to \(\ket{0}\). The input state is encoded in modes \(2\) and \(3\), the blank in modes \(1\) and \(4\), and photon pairs are injected into modes \(2\) and \(4\). The trainable ansatz is a \(4\)-mode universal interferometric mesh embedded in a \(6\)-mode universal photonic integrated circuit; the orange variational region optimizes \(12\) internal phases. Output clones are identified as dual-rail qubits on modes \((1,2)\) and \((3,4)\). Rather than performing full tomography at each iteration, the chip is configured to project each clone onto the target input-state basis, and fidelities are estimated from post-selected two-fold coincidence events [2407.06026].

For phase-covariant cloning in that device, the cost function is
\[
\mathcal{C}_{PC}=
\mathds{E}\big[(1-F_{1,\phi})^2+(1-F_{2,\phi})^2+(F_{1,\phi}-F_{2,\phi})^2\big].
\]
Because the cost contains at most second moments of the fidelities, the continuous equatorial average is replaced exactly by a minimal planar \(2\)-design with four states,
\[
\mathcal F=\left[0,\frac{\pi}{2},\pi,\frac{3\pi}{2}\right],
\]
namely the eigenstates of Pauli \(X\) and \(Y\). For state-dependent cloning, an additional regularization term involving post-selection probabilities \(P_A,P_B\) is added,
\[
\lambda\left[(1-P_A)^2+(1-P_B)^2+(P_A-P_B)^2\right],
\qquad \lambda=1,
\]
to suppress pathological solutions with high conditional fidelity but vanishing success probability. The classical optimizer is the gradient-free Nelder–Mead simplex method over \(\vec\theta\in[0,2\pi]^{12}\), augmented by a reboot strategy that enlarges the simplex around the best point found so far if noise and finite sampling induce premature collapse [2407.06026].

The photonic experiment demonstrates near-optimal rather than exact learning. For phase-covariant cloning, training on the four-state design and validating on \(50\) evenly spaced equatorial test states leads to convergence after roughly \(700\) iterations; the paper states that the machine reaches an almost optimal circuit in under \(1000\) iterations. The best measured average fidelities on the \(50\)-state test set are
\[
F_1^{opt}=0.802\pm 0.001,\qquad
F_2^{opt}=0.838\pm 0.001,
\]
which are below \(F_O=0.853\) but above the semiclassical threshold \(F_{SC}=0.750\). The reported asymmetry indicates only approximate symmetry of the learned solution, despite the symmetry penalty in the cost. The authors attribute the residual gap mainly to incomplete convergence induced by experimental noise and finite counting statistics [2407.06026].

## 4. Symmetry-restricted and physical cloning constructions

A conceptually different use of “wave function cloning” appears in mirror modular cloning for quantum associative memory. The memory state is
\[
|M\rangle=\frac{1}{\sqrt p}\sum_{i=1}^{p}|p^i\rangle,
\]
and its mirrored version is
\[
|\bar M\rangle=\frac{1}{\sqrt p}\sum_{i=1}^{p}|\bar p^i\rangle
= X^{\otimes n}|M\rangle,
\qquad
\bar p^i_j=1-p^i_j.
\]
The proposed state-dependent unitary acts only on the two-state set spanned by \(|M\rangle\) and \(|\bar M\rangle\), with amplitudes determined by the overlap \(\langle M|\bar M\rangle\). On input \(|M\rangle|\Sigma\rangle|0\rangle\), the unitary produces a flagged superposition whose ancilla indicates whether the output branch carries \( |M\rangle \) or \( |\bar M\rangle \). Because retrieval from mirrored memory is equivalent to retrieval from the original memory followed by a global bitwise inversion, the paper describes this as “perfectly cloned up to global mirroring.” The caveat is essential: this is not a universal wave-function cloner, but a restricted symmetry-assisted transformation defined on a highly structured family of states [2206.01644].

The same paper couples this cloning step to a content-addressable retrieval algorithm. Using an input register \(|i\rangle\), a cloned memory register, and \(b\) control qubits, it computes Hamming distances and amplifies patterns according to
\[
P(i,p^k)=\frac{1}{Z}\cos^{2b}\!\left(\frac{\pi}{2n}d_H(i,p^k)\right).
\]
The number of amplitude-amplification iterations is
\[
C=\left[
\frac{1}{\frac{1}{p}\sum_{k=1}^{p}\cos^{2b}\!\left(\frac{\pi}{2n}d_H(i,p^k)\right)}
\right]^{1/2},
\]
and under an approximately uniform-pattern assumption the paper estimates
\[
C\sim (\pi b)^{1/4}.
\]
Its claim of “perfect cloning” is therefore inseparable from associative-memory semantics: the mirrored branch is acceptable only because the downstream task treats \( |M\rangle \) and \( |\bar M\rangle \) as operationally equivalent up to a final inversion [2206.01644].

An older but still relevant strand concerns phase-covariant optical cloning of polarization qubits. There the physical transformation is generated by
\[
H=i\hbar\Gamma\left((a_{\varphi}^{\dagger})^2+(a_{\varphi^\perp}^{\dagger})^2\right)+\mathrm{h.c.},
\]
with
\[
g=\int dt\,\Gamma,\qquad C_g=\cosh g,\qquad T_g=\tanh g.
\]
The input family is the equator of the Poincaré sphere, expressed through
\[
a_{\varphi}^{\dagger}=\frac{1}{\sqrt2}\left(e^{i\varphi}a_H^\dagger+e^{-i\varphi}a_V^\dagger\right),
\qquad
a_{\varphi^\perp}^{\dagger}=
\frac{i}{\sqrt2}\left(e^{i\varphi}a_H^\dagger-e^{-i\varphi}a_V^\dagger\right).
\]
Amplification of one photon of a polarization singlet yields macroscopic states
\[
|\Phi\rangle=\sum_{i,j=0}^{\infty}\gamma_{ij}\bigl|(2i+1)_\varphi,(2j)_{\varphi^\perp}\bigr\rangle,
\qquad
|\Phi_\perp\rangle=\sum_{i,j=0}^{\infty}\gamma_{ij}\bigl|(2j)_\varphi,(2i+1)_{\varphi^\perp}\bigr\rangle,
\]
and a micro-macro singlet
\[
|\Psi^-\rangle=
\frac{1}{\sqrt2}
\Bigl(
|1_\varphi\rangle_a|\Phi_\perp\rangle_b-
|1_{\varphi^\perp}\rangle_a|\Phi\rangle_b
\Bigr).
\]
The associated numerical paper makes clear that the algorithmic contribution is primarily the simulation of these macrostates, not a general cloning algorithm [1112.0632].

That numerical work addresses slowly convergent hypergeometric sums governing probabilities, mean photon number, variance, and distinguishability. Its methods include logarithmic evaluation of factorial expressions, dynamic cutoff estimation, reordering of nested sums, switching between complementary summation formulas depending on threshold size, and tile-based parallelization with OpenMP. The reported implementation uses about \(200\) MB RAM for a \(1000\times1000\) tile and requires \(1\)–\(2\) hours runtime per such tile. In this context, “wave function cloning algorithm” can denote the computational machinery used to model the output of a restricted optical cloner rather than the cloning map itself [1112.0632].

## 5. Reformulation as channel cloning

A more general and conceptually unifying reformulation treats state cloning as a special case of deterministic cloning of quantum channels. Given a state family \(\mathcal S\), define the corresponding trash-and-replace channels
\[
\mathcal T_\rho[\cdot]=\rho,
\qquad \rho\in\mathcal S.
\]
The central equivalence statement is
\[
F^{\mathcal C}_{\rm CJ}=F^{\mathcal C}_{\diamond}=F^{\mathcal S},
\]
where \(\mathcal C=\{\mathcal T_\rho\}\), and any state-cloning map \(P\) induces a channel-cloning process \(\mathbf P\) through
\[
\mathbf P[\mathcal T_\rho^{\times N}]
=
\mathcal T_{P[\rho^{\otimes N}]}.
\]
For trash-and-replace channels, the paper further proves that parallel, sequential, and non-causal higher-order processes collapse to the same power as an ordinary CPTP state-cloning map. In this sense, wave-function cloning is exactly the trash-and-replace corner of a broader theory of channel replication [2509.08059].

The channel framework replaces ordinary circuits by higher-order quantum operations, or processes, that take \(N\) queries of an unknown channel \(E\) and return a channel approximating \(E^{\otimes M}\). The quality criteria are worst-case Choi–Jamiołkowski fidelity,
\[
F_{\rm CJ}(A,B)=F(\mathtt{CJ}[A],\mathtt{CJ}[B]),
\]
and worst-case diamond fidelity,
\[
F_\diamond(A,B)=
\min_{\Psi}
F\bigl(({\rm id}\otimes A)[\Psi],({\rm id}\otimes B)[\Psi]\bigr).
\]
The optimal worst-case channel-cloning fidelity is then
\[
F^{\mathcal C}_\bullet(N,M)=
\max_{\mathbf P}\inf_{E\in\mathcal C}
F_\bullet\big(\mathbf P[E^{\times N}],E^{\otimes M}\big),
\qquad
\bullet\in\{{\rm CJ},\diamond\}.
\]
This gives a single formalism for state cloning, unitary-gate replication, noisy-channel cloning, and super-replication questions [2509.08059].

Within that framework, the paper provides an SDP-based method to search for optimal cloning processes. Fidelity is represented by the semidefinite characterization
\[
F(\rho,\sigma)=
\max_Y \ \operatorname{tr}\frac{1}{2}(Y+Y^\dagger)
\quad\text{subject to}\quad
\begin{pmatrix}
\rho & Y\\
Y^\dagger & \sigma
\end{pmatrix}\ge 0.
\]
A finite net of target channels is chosen, a process Choi operator is optimized under complete-positivity and process constraints, and measure-and-prepare protocols can be approximately isolated by adding a PPT condition. The paper also draws a direct connection between optimal channel cloning and Bayesian channel estimation, and shows that measure-and-prepare super-replication can already achieve the optimal quadratic scaling for some unitary-channel families, while deterministic state cloning remains limited to linear rate [2509.08059].

## 6. Computational hardness, limitations, and recurrent misconceptions

A distinct line of work asks whether cloning may remain hard even when the cloner is given classical side information specifying a verifier circuit. This is not the standard no-cloning setting. The task is: given one copy of a witness state \(\ket{\psi}\) and a classical description of a polynomial-size verification circuit \(V\), output two valid copies, or more generally map the accepting space \(\mathcal L\) of \(V\) to \(\mathcal L\otimes\mathcal L\). The paper develops a conditional hardness framework using hidden maximally entangled states over representation-theoretically defined subspaces. For a verifier \(V\) uniquely accepting a hidden-EPR state \(\ket{\Phi_\Pi}\) with completeness \(c=1\) and soundness \(s\le 8/9\), it conjectures that if a uniform polynomial-time algorithm performs
\[
\ket{\Phi_\Pi}\ket{0}\mapsto \ket{\Phi_\Pi}^{\otimes 2},
\]
then there must also exist a polynomial-time algorithm that outputs some state in \(\Pi\). From this, and assuming \(\BQP\not\supseteq \NP\), the paper derives a conditional theorem excluding efficient cloning for a representation-theoretic family of unique-witness verifiers in the case \(m_{\mu\nu\lambda}=1\) [2411.11805].

This hardness result is computational rather than information-theoretic. It does not say that exact cloning is newly forbidden by linearity; instead, it argues that even when a classical verifier description weakens the usual “unknown state” premise, efficient cloning may still be impossible because witness generation or hidden-subspace reconstruction is computationally hard. The construction relies on projectors
\[
\Xi_\lambda=\frac{d_\lambda}{|G|}\sum_{g\in G}\chi^\lambda(g)^*\,\sigma(g),
\]
whose dimensions encode Kronecker multiplicities \(m_{\mu\nu\lambda}\), and on hidden maximally entangled states
\[
\ket{\Phi_\Pi}=
\frac{1}{\sqrt{d_1}}\sum_{i=1}^{d_1}\ket{b_i}\otimes \ket{b_i^*}.
\]
The article’s relevance to “wave function cloning algorithms” is therefore negative but important: it delineates a regime in which algorithmic cloning may fail for reasons of complexity even after the strict information-theoretic unknown-state barrier has been softened [2411.11805].

Several misconceptions recur across the literature. First, no paper surveyed here provides a loophole around the no-cloning theorem for arbitrary unknown states. Second, “perfect cloning” can mean exactness only within a restricted equivalence class, as in mirror modular cloning up to global bitwise complementation. Third, “optimal cloning” usually means attaining the maximum fidelity allowed by quantum mechanics for a chosen ensemble, not fidelity \(1\). Fourth, some papers concern the search for a cloner or the simulation of a physical cloner rather than a standalone closed-form cloning circuit. A plausible implication is that “wave function cloning algorithm” is best understood as an umbrella label for restricted, task-specific, and often hardware- or application-adapted procedures, rather than as the name of a single universally accepted algorithmic object.

Source: https://www.emergentmind.com/topics/wave-function-cloning-algorithm