Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2603.14576

IQP Born Machines under Data-dependent and Agnostic Initialization Strategies

Sacha Lerch, Joseph Bowles, Ricard Puig, Erik Armengol, Zoë Holmes, Supanut Thanasilp

quant-phcs.LGstat.MLSubmitted Mar 15, 2026• Announced Mar 17, 2026• 1 citation
arXivHTMLPDF

Abstract

Quantum circuit Born machines based on instantaneous quantum polynomial-time (IQP) circuits are natural candidates for quantum generative modeling, both because of their probabilistic structure and because IQP sampling is provably classically hard in certain regimes. Recent proposals focus on training IQP-QCBMs using Maximum Mean Discrepancy (MMD) losses built from low-body Pauli-$Z$ correlators, but the effect of initialization on the resulting optimization landscape remains poorly understood. In this work, we address this by first proving that the MMD loss landscape suffers from barren plateaus for random full-angle-range initializations of IQP circuits. We then establish lower bounds on the loss variance for identity and an unbiased data-agnostic initialization. We then additionally consider a data-dependent initialization that is better aligned with the target distribution and, under suitable assumptions, yields provable gradients and generally converges quicker to a good minimum (as indicated by our training of circuits with 150 qubits on genomic data). Finally, as a by-product, the developed variance lower bound framework is applicable to a general class of non-linear losses, offering a broader toolset for analyzing warm-starts in quantum machine learning.

Comments
16 + 35 pages, 3 + 4 figures

Related papers

Trainability of IQP Quantum Circuit Born Machines Under Gaussian Initialization

Gennaro De Luca, Vinayak Sharma, Aviral Shrivastava

Quantum Circuit Born Machines (QCBMs) offer a natural approach to generative machine learning by leveraging the Born rule. Recent work has provided a method to classically train QCBMs with Instantaneous Quantum Polynomial (IQP) circuits via the Maximum Mean Discrepancy (MMD) loss. Despite the assumed intractability of sampling from IQP circuits classically, their expectation values can be computed classically, enabling training of these IQP QCBMs. However, quantum machine learning (QML) models have various other challenges, including trainability issues caused by exponential concentration or barren plateaus. While these issues have been explored for parameters sampled from a uniform distribution, little work has been done to rigorously treat the use of arbitrary Gaussian initialization schemes. This work leverages Stein's lemma and Lipschitz concentration bounds for Gaussian random variables to provide an analytical lower bound of the variance of the gradient and a probabilistic concentration bound of the deviation of the gradient from its mean. It discusses strategies to either avoid or encourage exponential concentration, as well as the conditions under which barren plateaus are more likely to occur.

quant-pharXiv:2606.10179•Jun 10, 2026•0 citations

Characterizing Trainability of Instantaneous Quantum Polynomial Circuit Born Machines

Kevin Shen, Susanne Pielawa, Vedran Dunjko +1 more

Instantaneous quantum polynomial quantum circuit Born machines (IQP-QCBMs) have been proposed as quantum generative models with a classically tractable training objective based on the maximum mean discrepancy (MMD) and a potential quantum advantage motivated by sampling-complexity arguments, making them an exciting model worth deeper investigation. While recent works have further proven the universality of a (slightly generalized) model, the next immediate question pertains to its trainability, i.e., whether it suffers from the exponentially vanishing loss gradients, known as the barren plateau issue, preventing effective use, and how regimes of trainability overlap with regimes of possible quantum advantage. Here, we provide significant strides in these directions. To study the trainability at initialization, we analytically derive closed-form expressions for the variances of the partial derivatives of the MMD loss function and provide general upper and lower bounds. With uniform initialization, we show that barren plateaus depend on the generator set and the spectrum of the chosen kernel. We identify regimes in which low-weight-biased kernels avoid exponential gradient suppression in structured topologies. Also, we prove that a small-variance Gaussian initialization ensures polynomial scaling for the gradient under mild conditions. As for the potential quantum advantage, we further argue, based on previous complexity-theoretic arguments, that sparse IQP families can output a probability distribution family that is classically intractable, and that this distribution remains trainable at initialization at least at lower-weight frequencies.

quant-pharXiv:2602.11042•Feb 12, 2026•0 citations

Trainability and Mode Separation of Mixed IQP-QCBMs

Youngseok Lee, Hyunwoo Kim

Quantum circuit Born machines (QCBMs) based on instantaneous quantum polynomial-time (IQP) circuits are promising quantum generative models for their classical trainability. It is known that their ancilla-free form avoids barren plateaus under certain initializations, but remains non-universal. Although adding ancilla qubits raises the expressivity, whether the ancilla-extended model retains local trainability remains unknown. We propose the mixed IQP-QCBM, which generalizes the ancilla-extended circuit as a weighted mixture of ancilla-free IQP circuits, called branches. For a polynomial number of branches, we prove local barren-plateau avoidance from data-agnostic and, under certain assumptions, data-dependent initializations. We further show that the mixed IQP-QCBM can surpass the best ancilla-free IQP circuit only if its branches generate a number of distinct distributions. In particular, we focus on a behavior we call \emph{mode separation}, in which each branch captures a particular feature of the target. Mode separation is hard to attain from an initialization whose branches generate the same distribution: the gradients that would separate them are suppressed while the distributions they generate remain close. This motivates \emph{cluster initialization}, which assigns a different unsupervised data cluster to each branch and provides an initial degree of mode separation. Exact calculations on two 16-bit datasets support the barren-plateau and gradient-suppression claims. On four benchmarks, binary clusters, a two-dimensional Ising model, binarized MNIST, and a 484-spin glass, cluster initialization converges fastest and reaches the lowest mean test $\mathrm{MMD}^2$. We observe that, when achieving the lowest test $\mathrm{MMD}^2$, the mixed IQP-QCBM contains branches specialized to distinguishable data features such as blob patterns, magnetization sectors, or digit shapes.

quant-pharXiv:2607.27883•Jul 31, 2026•0 citations

Deployment of Large IQP Circuit Born Machines on Qubit-Limited Hardware Using Noise Injection

Ju-Young Ryu, Wooyeong Song, Kwangil Bae +2 more

Quantum circuit Born machines (QCBMs) built with instantaneous quantum polynomial (IQP) circuits have loss functions that can be efficiently estimated classically, but require quantum computers during inference. Thus, classical computers can train large IQP QCBMs that exceed the qubit numbers of available quantum hardware. This motivates research into how qubit-limited quantum computers can be used to sample from larger IQP QCBMs. Leveraging insights from classical simulation of noisy IQP circuits, we propose a deployment protocol for IQP QCBMs that injects noise at the logical level to probabilistically decompose the original circuit into smaller clusters that can be sampled on qubit-limited quantum processors. We also apply lightweight resampling of known noisy qubits to mitigate the impact of the introduced error. Numerical simulations of our protocol applied to shallow IQP QCBMs and IBM Heron experiments deploying a 316-qubit IQP circuit onto 156-qubit hardware show that mitigation improves the similarity of the empirical distribution to the target distribution. The proposed deployment protocol also crucially preserves correlation structures more accurately than a classical baseline that has accurate single-qubit marginals.

quant-pharXiv:2609.35163•Sep 29, 2026•0 citations

An IQP Born Machine for Calorimeter Image Generation at 64 Qubits with Compiled-IQP Deployment

Jamal Slim, Saverio Monaco, Florian Rehm +2 more

The challenge to scaling quantum generative models on near-term hardware is training. Variational circuit Born machines require repeated quantum sampling and are prone to barren plateaus. Instantaneous Quantum Polynomial-time (IQP) Born machines sidestep both, since their loss is built from low-order Pauli-Z correlators that admit an unbiased classical estimator, while sampling from worst-case circuits in the class is conjectured to be classically hard. We take this train-on-classical, deploy-on-quantum workflow to a real high-energy-physics generative task, learning calorimeter shower profiles at 64 qubits and running the trained model on an IBM Heron r2 superconducting processor at 67 physical qubits. Three ingredients make it work. A uniform mixture of IQP circuits (MoIQP) widens the model class at single-circuit training cost. The Pearson-Stabilized Correlation Kernel (PSCK) biases descent toward the pairwise correlations that carry the shower-development physics, which the standard heat kernel systematically compresses. An exact deferred-measurement compilation collapses the mixture into a single IQP circuit, realized on hardware as a constant-depth dynamic circuit with zero SWAP insertions on the device's native heavy-hex graph. The trained model reconstructs the correlation structure to within 0.016 of the floor imposed by the encoding itself. Raw device samples reproduce the per-cell energy spectra and the full pairwise correlation structure at Pearson r = 0.989, up to a single global amplitude compression of depolarizing origin. A Gaussian copula fitted to the same training split matches the pairwise target more accurately than the quantum model at negligible cost. The contribution is therefore the classical trainability, exact compilation, and hardware deployability of a quantum generative model at this scale, not superiority over classical surrogates.

quant-pharXiv:2605.27735•May 28, 2026•0 citations

Universality and kernel-adaptive training for classically trained, quantum-deployed generative models

Andrii Kurkin, Shen, Kevin, Susanne Pielawa +2 more

The instantaneous quantum polynomial (IQP) quantum circuit Born machine (QCBM) has been proposed as a promising quantum generative model over bitstrings. Recent works have shown that the training of IQP-QCBM is classically tractable w.r.t. the so-called Gaussian kernel maximum mean discrepancy (MMD) loss function, while maintaining the potential of a quantum advantage for sampling itself. Nonetheless, the model has a number of aspects where improvements would be important for more general utility: (1) the basic model is known to be not universal - i.e. it is not capable of representing arbitrary distributions, and it was not known whether it is possible to achieve universality by adding hidden (ancillary) qubits; (2) a fixed Gaussian kernel used in the MMD loss can cause training issues, e.g., vanishing gradients. In this paper, we resolve the first question and make decisive strides on the second. We prove that for an $n$-qubit IQP generator, adding $n + 1$ hidden qubits makes the model universal. For the latter, we propose a kernel-adaptive training method, where the kernel is adversarially trained. We show that in the kernel-adaptive method, the convergence of the MMD value implies weak convergence in distribution of the generator. We also analytically analyze the limitations of the MMD-based training method. Finally, we verify the performance benefits on the dataset crafted to spotlight improvements by the suggested method. The results show that kernel-adaptive training outperforms a fixed Gaussian kernel in total variation distance, and the gap increases with the dataset dimensionality. These modifications and analyses shed light on the limits and potential of these new quantum generative methods, which could offer the first truly scalable insights in the comparative capacities of classical versus quantum models, even without access to scalable quantum computers.

arXiv:2510.08476•Oct 9, 2025•0 citations
Quantum circuit Born machines based on instantaneous quantum polynomial-time (IQP) circuits are natural candidates for quantum generative modeling, both because of their probabilistic structure and because IQP sampling is provably classically hard in certain regimes. Recent proposals focus on training IQP-QCBMs using Maximum Mean Discrepancy (MMD) losses built from low-body Pauli-ZZZ correlators, but the effect of initialization on the resulting optimization landscape remains poorly understood. In this work, we address this by first proving that the MMD loss landscape suffers from barren plateaus for random full-angle-range initializations of IQP circuits. We then establish lower bounds on the loss variance for identity and an unbiased data-agnostic initialization. We then additionally consider a data-dependent initialization that is better aligned with the target distribution and, under suitable assumptions, yields provable gradients and generally converges quicker to a good minimum (as indicated by our training of circuits with 150 qubits on genomic data). Finally, as a by-product, the developed variance lower bound framework is applicable to a general class of non-linear losses, offering a broader toolset for analyzing warm-starts in quantum machine learning.
Quantum Circuit Born Machines (QCBMs) offer a natural approach to generative machine learning by leveraging the Born rule. Recent work has provided a method to classically train QCBMs with Instantaneous Quantum Polynomial (IQP) circuits via the Maximum Mean Discrepancy (MMD) loss. Despite the assumed intractability of sampling from IQP circuits classically, their expectation values can be computed classically, enabling training of these IQP QCBMs. However, quantum machine learning (QML) models have various other challenges, including trainability issues caused by exponential concentration or barren plateaus. While these issues have been explored for parameters sampled from a uniform distribution, little work has been done to rigorously treat the use of arbitrary Gaussian initialization schemes. This work leverages Stein's lemma and Lipschitz concentration bounds for Gaussian random variables to provide an analytical lower bound of the variance of the gradient and a probabilistic concentration bound of the deviation of the gradient from its mean. It discusses strategies to either avoid or encourage exponential concentration, as well as the conditions under which barren plateaus are more likely to occur.
Instantaneous quantum polynomial quantum circuit Born machines (IQP-QCBMs) have been proposed as quantum generative models with a classically tractable training objective based on the maximum mean discrepancy (MMD) and a potential quantum advantage motivated by sampling-complexity arguments, making them an exciting model worth deeper investigation. While recent works have further proven the universality of a (slightly generalized) model, the next immediate question pertains to its trainability, i.e., whether it suffers from the exponentially vanishing loss gradients, known as the barren plateau issue, preventing effective use, and how regimes of trainability overlap with regimes of possible quantum advantage. Here, we provide significant strides in these directions. To study the trainability at initialization, we analytically derive closed-form expressions for the variances of the partial derivatives of the MMD loss function and provide general upper and lower bounds. With uniform initialization, we show that barren plateaus depend on the generator set and the spectrum of the chosen kernel. We identify regimes in which low-weight-biased kernels avoid exponential gradient suppression in structured topologies. Also, we prove that a small-variance Gaussian initialization ensures polynomial scaling for the gradient under mild conditions. As for the potential quantum advantage, we further argue, based on previous complexity-theoretic arguments, that sparse IQP families can output a probability distribution family that is classically intractable, and that this distribution remains trainable at initialization at least at lower-weight frequencies.
Quantum circuit Born machines (QCBMs) based on instantaneous quantum polynomial-time (IQP) circuits are promising quantum generative models for their classical trainability. It is known that their ancilla-free form avoids barren plateaus under certain initializations, but remains non-universal. Although adding ancilla qubits raises the expressivity, whether the ancilla-extended model retains local trainability remains unknown. We propose the mixed IQP-QCBM, which generalizes the ancilla-extended circuit as a weighted mixture of ancilla-free IQP circuits, called branches. For a polynomial number of branches, we prove local barren-plateau avoidance from data-agnostic and, under certain assumptions, data-dependent initializations. We further show that the mixed IQP-QCBM can surpass the best ancilla-free IQP circuit only if its branches generate a number of distinct distributions. In particular, we focus on a behavior we call \emph{mode separation}, in which each branch captures a particular feature of the target. Mode separation is hard to attain from an initialization whose branches generate the same distribution: the gradients that would separate them are suppressed while the distributions they generate remain close. This motivates \emph{cluster initialization}, which assigns a different unsupervised data cluster to each branch and provides an initial degree of mode separation. Exact calculations on two 16-bit datasets support the barren-plateau and gradient-suppression claims. On four benchmarks, binary clusters, a two-dimensional Ising model, binarized MNIST, and a 484-spin glass, cluster initialization converges fastest and reaches the lowest mean test MMD2\mathrm{MMD}^2MMD2. We observe that, when achieving the lowest test MMD2\mathrm{MMD}^2MMD2, the mixed IQP-QCBM contains branches specialized to distinguishable data features such as blob patterns, magnetization sectors, or digit shapes.
Quantum circuit Born machines (QCBMs) built with instantaneous quantum polynomial (IQP) circuits have loss functions that can be efficiently estimated classically, but require quantum computers during inference. Thus, classical computers can train large IQP QCBMs that exceed the qubit numbers of available quantum hardware. This motivates research into how qubit-limited quantum computers can be used to sample from larger IQP QCBMs. Leveraging insights from classical simulation of noisy IQP circuits, we propose a deployment protocol for IQP QCBMs that injects noise at the logical level to probabilistically decompose the original circuit into smaller clusters that can be sampled on qubit-limited quantum processors. We also apply lightweight resampling of known noisy qubits to mitigate the impact of the introduced error. Numerical simulations of our protocol applied to shallow IQP QCBMs and IBM Heron experiments deploying a 316-qubit IQP circuit onto 156-qubit hardware show that mitigation improves the similarity of the empirical distribution to the target distribution. The proposed deployment protocol also crucially preserves correlation structures more accurately than a classical baseline that has accurate single-qubit marginals.
The challenge to scaling quantum generative models on near-term hardware is training. Variational circuit Born machines require repeated quantum sampling and are prone to barren plateaus. Instantaneous Quantum Polynomial-time (IQP) Born machines sidestep both, since their loss is built from low-order Pauli-Z correlators that admit an unbiased classical estimator, while sampling from worst-case circuits in the class is conjectured to be classically hard. We take this train-on-classical, deploy-on-quantum workflow to a real high-energy-physics generative task, learning calorimeter shower profiles at 64 qubits and running the trained model on an IBM Heron r2 superconducting processor at 67 physical qubits. Three ingredients make it work. A uniform mixture of IQP circuits (MoIQP) widens the model class at single-circuit training cost. The Pearson-Stabilized Correlation Kernel (PSCK) biases descent toward the pairwise correlations that carry the shower-development physics, which the standard heat kernel systematically compresses. An exact deferred-measurement compilation collapses the mixture into a single IQP circuit, realized on hardware as a constant-depth dynamic circuit with zero SWAP insertions on the device's native heavy-hex graph. The trained model reconstructs the correlation structure to within 0.016 of the floor imposed by the encoding itself. Raw device samples reproduce the per-cell energy spectra and the full pairwise correlation structure at Pearson r = 0.989, up to a single global amplitude compression of depolarizing origin. A Gaussian copula fitted to the same training split matches the pairwise target more accurately than the quantum model at negligible cost. The contribution is therefore the classical trainability, exact compilation, and hardware deployability of a quantum generative model at this scale, not superiority over classical surrogates.
The instantaneous quantum polynomial (IQP) quantum circuit Born machine (QCBM) has been proposed as a promising quantum generative model over bitstrings. Recent works have shown that the training of IQP-QCBM is classically tractable w.r.t. the so-called Gaussian kernel maximum mean discrepancy (MMD) loss function, while maintaining the potential of a quantum advantage for sampling itself. Nonetheless, the model has a number of aspects where improvements would be important for more general utility: (1) the basic model is known to be not universal - i.e. it is not capable of representing arbitrary distributions, and it was not known whether it is possible to achieve universality by adding hidden (ancillary) qubits; (2) a fixed Gaussian kernel used in the MMD loss can cause training issues, e.g., vanishing gradients. In this paper, we resolve the first question and make decisive strides on the second. We prove that for an nnn-qubit IQP generator, adding n+1n + 1n+1 hidden qubits makes the model universal. For the latter, we propose a kernel-adaptive training method, where the kernel is adversarially trained. We show that in the kernel-adaptive method, the convergence of the MMD value implies weak convergence in distribution of the generator. We also analytically analyze the limitations of the MMD-based training method. Finally, we verify the performance benefits on the dataset crafted to spotlight improvements by the suggested method. The results show that kernel-adaptive training outperforms a fixed Gaussian kernel in total variation distance, and the gap increases with the dataset dimensionality. These modifications and analyses shed light on the limits and potential of these new quantum generative methods, which could offer the first truly scalable insights in the comparative capacities of classical versus quantum models, even without access to scalable quantum computers.