Skip to content
Back to Academy

What Is Quantum Machine Learning?

Quantum machine learning, explained for the model owner: can quantum feature extraction move a model that has stopped improving, and do you need a quantum computer in production? The machine-learning session of the Kipu Academy business track.

TutorialIntermediate~20 minHands-on lab
Get in touch

Educational disclaimer. The method is described by analogy; the papers cited under each claim carry the scientific detail. Service limits, free tiers and hardware names are snapshots from the day this was written and will drift. Treat the Kipu Quantum Hub and the linked papers as the single sources of truth.

1. Place the quantum part

Why the processor belongs in the training path only, and what that means for serving.

2. See which models gain

The controlled experiment: fed the expert's features it beats them, denied them it finds nothing.

3. Judge the evidence

The single-split trap, out-of-fold scoring with correction, and what the published results support.

The model you have already tuned flat

Somewhere in your organisation there is a machine learning model that someone defends in a monthly review. It scores credit applications, or flags a machine before it fails, or ranks compounds before the lab spends a week to synthesise them. It works. It has worked for two years. And over the last four quarters nobody was able to improve its accuracy.

The obvious things don't work anymore: more capacity, more trees, a longer hyperparameter search, class weights, calibration. Each pass costs a week for marginal gains. That is a signal that the constraint is no longer the model and rather sits in how the data is presented. Your plateau is a re-presentation limit, not a tuning limit.

A learning curve rising steeply at first, then going flat. The early stretch is labelled easy early gains. The flat stretch, hatched and marked YOUR PLATEAU, is labelled more effort, no more accuracy. The horizontal axis is your own effort: more data, more capacity, longer search. The vertical axis is model performance. A note beside the flat stretch reads FLAT, tuning searches only the features you already have.

This session introduces an instrument built for those cases: Kipu's Rimay service, which computes extra input columns for a model you already tuned flat. However, it does not remove the tuning work, it can help ML models digest the data better. Over the last year we have seen Rimay work particularly well for problems where training data is rare.

For example four hundred failures recorded in nine years, the ninety toxic compounds, the few customers with high churn intent, the hundred defaults in a portfolio of overwhelmingly good loans. You can subscribe to it on the Kipu Quantum Hub Marketplace.

How does it work, where's the quantum computer?

First up: Using Rimay, your production system will not depend on a quantum computer. It is needed only for the training phase. Inference stays classical afterwards. Quantum machine learning is the use of a quantum processor inside a machine-learning pipeline. Whether to represent the data, to fit the model, or to make the prediction. Rimay, also called digitized quantum feature extraction, falls in the first category. It writes existing data columns onto qubits and returns quantum fit transform measurements for ordinary machine learning models.

quan·tum ma·chine learn·ing

/ˈkwɒn.təm məˈʃiːn ˈlɜː.nɪŋ/noun

  1. 1

    the use of a quantum processor somewhere inside a machine-learning pipeline, whether to represent the data, to fit the model, or to make the prediction.

End-to-end quantum machine learning would put a processor in the pipeline. So every training, every inference waits in a hardware queue on the cloud. That is unusable in industry: latency measured in minutes rather than milliseconds, on-premise policy prohibits cloud usage in production. Feature extraction can put the quantum processor in the training path only. Slow, scheduled, expensive compute during training is something organisations already tolerate. If required, a surrogate, a small classical model trained to reproduce the quantum output, extracts the performance gains and makes them usable. More on that later.

Two serving architectures drawn as flow diagrams. Path A, struck through with a heavy diagonal cross, runs data into a model and then out through a QPU on the way to every prediction, putting the processor in the request path. Path B leaves the request path classical, data to model to prediction in a straight line, with the QPU sitting above the model and looping back into it, used during training only.
Path A makes every prediction dependent on hardware availability. Path B keeps quantum out of inference.

Here's how Rimay works: raw data, scarce and expensive to label, goes into the quantum feature extraction service, which returns enriched columns. Your existing scikit-learn, XGBoost or PyTorch models consume them alongside existing data. Your model and stack remain mostly unchanged with one additional API call in between.

One row through three stages labelled input data, Rimay extraction and output data. Input data holds nine grey column strips, labelled nine feature columns, one per sensor reading. Rimay extraction holds a three by three lattice of nine grey ringed sites joined by faint grey links, with a single gold path snaking through every site in a shape like a two. Output data holds seventeen columns: nine grey ones labelled one per feature, then eight gold ones labelled feature pairs.
What goes in: nine feature columns, one qubit each. Where they sit: a schematic lattice, grey sites with a single gold path winding through them. What comes back: seventeen columns, nine per-feature and eight feature-pair.

The encoding uses one qubit per feature. Nine columns (we call them features) of a predictive-maintenance tabular dataset are mapped onto nine qubits. Each feature's value sets that qubit's bias. Quantum columns come back: one measured value per qubit, plus feature-pair fit transforms. The free simulator allows usage for up to fifteen features on the input dataset. Paid services allow usage on real quantum computers with up to 156 features.

The quantum stage of the published pipeline and the classical classifier it feeds. On the left, a bracketed block labelled digitized quantum feature mapping: your feature vector sets a Hamiltonian encoding, a short circuit of gates runs over four qubit wires and is measured, and the measured expectation values are collected into a vector. On the right, a bracketed block labelled classical classifier: that vector becomes one row of an ordinary table, which support-vector machines and ensemble models consume.
The Digitized Quantum Feature Extraction protocol as drawn by the Kipu authors in published work. DQFE, is the algorithm Rimay implements as a service. Figure 1 of Zhang, Q. et al., arXiv:2602.18350.

In a nutshell, Rimay re-presents structure in a form models can digest better if mixed with the original columns. It does not engineer data, or recover missing information. Overall these fit transforms help mostly with feature separation. Whether a model can use the re-presentation is a separate question and is one that we discover in the following paragraphs.

Which models improve in performance?

An important question, and perhaps the biggest misconception about Rimay is that it replaces tuning work of a data engineer. To investigate whether that is the case, we start two Rimay extractions on the same 3,000 rows of AI4I 2020, a public predictive-maintenance benchmark. All other things being equal (same split), only one of Rimay's inputs contained the three additional features a domain expert had hand-built: temperature difference, power and overstrain features.

Two arms of one experiment, drawn as two rows running left to right through three stages labelled input data, Rimay extraction and output data. In the first row the input data holds six solid grey raw columns and then three visibly empty dashed slots where the expert features would be; in the second row those same three slots are filled by outlined columns, the expert's three hand-built features, so the two input tables differ on sight. The Rimay extraction stage is drawn identically in both rows as a single gold chain of nine qubits snaking through three rows. The output data holds three labelled families side by side: plain grey classical columns, outlined expert columns and gold quantum columns. The denied row gets eleven quantum columns back from its six inputs, the fed row seventeen from its nine, so the fed row's output table is visibly wider, with the expert family showing as three empty dashed slots in the denied row.
All things being equal (same 3,000 rows, same split, same models). The schematics lay out two pipelines, one with only the raw sensor data (top) and one with the expert's engineered features added (bottom).

Denied the expert features, no model shows a significant gain. Fed them, the re-presentation improves the expert's data further for Naive Bayes. Gains: +0.0394 average precision over the expert-augmented table itself, at p 6.1e-05 after Holm correction, that means the result is statistically significant (p<0.05). This improvement does not hold for all models. Logistic regression: +0.0180. Linear SVM: no effect. RBF-kernel SVM: -0.0506. Gradient boosting: -0.0476. We can recover the losses with intelligent feature selection later.

Bar chart, one model, four arms, naive Bayes on AI4I. Raw sensors alone reach 0.549 average precision. Raw plus Rimay denied the expert physics reaches 0.547, a change of zero, annotated Rimay finds nothing. Raw plus the expert physics reaches 0.665. The expert table plus Rimay re-presented reaches 0.704, annotated plus 0.039 on the expert table at p 6.1e-05 Holm.
Denied the expert's columns, the extraction returns nothing usable for this model; fed them, it re-presents them so a model can read them better. 3,000 rows, results taken from the free Rimay Simulator on hub.kipu-quantum.com, 2,000 shots.

And this comes at no surprise: a well known effect that predates quantum computers. For interested readers, the following study generated sixteen families of engineered features and measured how well four regression models deal with them; random forests and gradient boosting work well on their own, having most benefit from only three (arXiv:1701.07852, v2, 2020). That also means that models require different tuning. If you want to investigate the tradeoffs in more detail you will look at a so-called PR-curve. This has nothing to do with marketing and everything with how well a model predicts outcomes under different settings.

Precision-recall curves for naive Bayes on the fed extraction, pooled out-of-fold predictions. The curve with Rimay's re-encoded columns, in gold, sits above the expert-table curve, in grey, across the range, a clean visible separation rather than curves trading places.
Validated curve separation from pooled out-of-fold estimates; paired per-split Average Precision +0.0394 at p 6.1e-05 Holm. AI4I 3,000 rows, free Rimay Simulator on hub.kipu-quantum.com, 2,000 shots paired with Naive Bayes.

Separation between the grey and yellow line is benefit or loss from the quantum features. Upper right means higher performance. Captured here is Recall (don't overlook anything) against Precision (avoid false alarms). ML models can run at different thresholds. If you want absolutely no false negatives you will need to run at high Recall thresholds which incentivise the model to label far more data as positives, at the cost of a lot of false alarms. At high precision thresholds you want the system to only ring the alarm when it is 100% certain. This comes at the expense of recall. The right setting sits somewhere in between. Depending on your threshold needs, the curve separation can be larger or smaller even for the same model and the same dataset.

Decision Logic: In this instance the cost of false positives and false negatives is highly asymmetric and the positive class is rare. Therefore the PR-curve is the appropriate metric to use. For a balanced dataset you'd go with a ROC-curve analysis (read more here).

Which datasets work well with Rimay?

To illustrate both gains and losses, this section discusses three datasets: AI4I (predictive maintenance), Pima (diabetes), UNSW-NB15 (network intrusion). After closer inspection it becomes clear that AI4I is the only table of the four that carried structure for the extraction to re-present.

Dot plot, change in average precision when Rimay columns are added, three datasets. AI4I predictive maintenance, where an expert hand-built the physics, shows naive Bayes plus 0.054, logistic regression plus 0.024, linear SVM plus 0.022. Pima diabetes shows naive Bayes minus 0.032, logistic regression minus 0.018, linear SVM minus 0.020. UNSW network traffic shows naive Bayes minus 0.006, logistic regression minus 0.051, linear SVM minus 0.026.

As you can see, Rimay is not a magic wand, but rather based on traceable data transformations. Measured gains are largest on a table that already carries structure, and losses appear on ones that do not. We just completed step 1 of 4 to determine whether we got a hit:

  1. Measure performance on raw Classical+Quantum features.
  2. Apply a Feature Selection Protocol to repair the losses.
  3. Perform a PR or ROC-Curve Analysis.
  4. Test for Statistical Significance.

Extracting gains with the Rimay simulator is considerably more difficult than with its larger quantum hardware counterpart. Limited size and inherent limited complexity of the small feature counts make it hard, yet not prohibitive to find gains.

A schematic tradeoff with no numbers on either axis. Columns added runs to the right along the horizontal axis, effect on the metric runs up the vertical axis. A gold curve, the signal the new columns add, rises steeply at first and then flattens. A grey curve, the variance cost of the extra width, rises steadily and keeps rising. The two cross, and from the crossing rightwards a hatched band labelled cost greater than impact marks the region past the crossing. A note at the right reads Balance: the impact is real, but it can be less than the columns cost.
Losses are caused by dilution can outweigh the benefits from the added signal. Every added features costs variance and needs investigation whether it is worth its cost.

Feature Selection (Step 2)

Using only training data an inner cross-validation decides which columns to keep. The so called Fisher score measures the variance of features depending on their class. High values are achieved if the within-class variance is low and the cross-class difference is large. That means features that separate classes well achieve higher scores than those that blur the lines. Other feature selection methods exist.

Three parts. On the left, nine grey column strips standing side by side, labelled all your columns, kept, with none removed. In the middle, seventeen thin gold strips labelled seventeen quantum columns pass down through a gold funnel, and a smaller set of gold strips labelled the best few, by Fisher, comes out below it. A dashed gold arrow carries those few across to the right, where the same nine grey strips appear unchanged with the few gold strips appended beside them, labelled the table the model sees. A note at the right reads keep everything you have, add only the best of the new.
The Fisher selection protocol. Ranking is computed on the training fold only, on inner cross-validation.
Grouped bar chart, naive Bayes, change in average precision against the raw columns alone, four datasets, two bars each. Grey bars, all quantum columns added unselected: AI4I plus 0.058, Pima minus 0.035, UNSW minus 0.017, Cleveland minus 0.033. Gold bars, keep every raw column and add only the top-Fisher quantum columns: AI4I plus 0.063, Pima plus 0.018, UNSW plus 0.008, Cleveland minus 0.002. Title reads: with rule-based selection, the losses disappear. Subtitle: zero is no change, not a win.
Same data from the free Rimay Simulator (hub.kipu-quantum.com, 500 shots). The losses disappear after selection.

Performance Curve Analysis (Step 3)

Once you have applied steps 1 and maybe 2, sometimes the perfomance score will improve and it is right to ask whether what you just observed is quantum advantage. That argument requires a little more nuance and background checks. Quantum Advantage carries a lot of importance and meaning in the world of computation and as soon as you will announce it, many will try to dispute it. The following two sections Step 3 and 4 will serve as an initial guide that help you determine whether a result is significant for your application. But at the same time, keep in mind that the ultimate end goal is industrial usefulness. So if it works for you and saves money effort or time, use it, and don't waste too much thought on artificial debates.

The following data comes from one held-out split of 140 carrying 50 failures of the AI4I training dataset, using logistic regression. From initial analysis failure-class F1 moves 4pp from 0.740 to 0.780. Accuracy improves by 2.9 points. Looking at the PR-plot: both curves sit nearly on top of each other on average. So where does the lift come from? An important lesson: Metrics you care about can still differ usefully for your application.

Precision-recall curves for one held-out split of the AI4I benchmark, 140 test rows containing 50 failures. The classical curve reaches average precision 0.843 and the curve with Rimay columns reaches 0.851. The two curves cross repeatedly and each step is a handful of rows changing side.
The two curves trade places all the way along, and the gap between the two average-precision numbers is small. 700 rows, 140 test samples computed with the free Rimay Simulator on hub.kipu-quantum.com, 500 shots.

Specific applications do not consume the average of a PR or ROC-curve. F1 and accuracy measure after you have chosen a threshold. Rimay's extracted features might help with separation in a certain threshold region, while providing less benefit in another. Depending on your application you can tune your thresholds: high recall settings apply when a false negative costs several times more than a false alarm. You then fix recall at, say, 0.8 and tune for maximum precision at that setting.

A schematic precision-recall diagram with no numbers on either axis. Precision runs up the vertical axis as a labelled direction, recall runs to the right. Two curves, one grey for the arm you run today and one gold for the candidate arm, fall together and sit almost on top of each other across most of the range. A hatched vertical band covers the high-recall end and is labelled your operating region. Inside that band the two curves separate, and a short gold bar marks the vertical distance between them, annotated as the gap you pay for. A note at the right reads: at the recall you must run at, the precision difference is the whole story.
How to read a curve. Schematic, no axis scale, and not the measured single-split figure above: on that campaign the region-first reading was not significant either way.

So determine your curve region-first: fix the recall your cost structure demands, read off the precision there, and price the difference in. While on the run above that reading was not significant either way, the point is the method. A whole-curve average can undersell or oversell a model for your specific decision.

Fair comparison of models

Machine Learning models are statistical and you can get a high performing model from structural advantage or just by luck. To determine whether we got structural advantage we analyse statistical significance with a McNemar test. Here classical and quantum models go head to head. Datapoints only the classical model got right award points to that model and inversely those only the hybrid did get right add to its points. Rows both models got right or both got wrong drop out; the two remaining scores provide the p-value.

In general:

  • p>0.05 not significant.
  • p<0.05 significant results.

A lot of times, if you use an 80/20 split, you only ever test on 20 percent of data at low statistical significance. Most of your data goes into training and is never scored and can provide the necessary statistics for more accurate p-values. Predicting every row exactly once, out of fold, from a model that never saw it, unlocks the entire dataset. By doing so, you also have to correct for the fact that you tried several models with a Holm correction. To give you a perspective, with five candidates, the chance that one of them clears p = 0.05 by luck alone is roughly one in four. Holm correction prevents luck reported as a result.

Five rows, each one fold of the same data split into five equal blocks. In the first row the leftmost block is gold, held out and scored, and the other four are hatched, trained on. Down the rows the gold block moves one place to the right, so across the five folds every block of the data is scored exactly once, always by a model that never saw it.
What every row measured actually means. Coverage goes from a fifth of your rows to all of them, and it costs compute rather than data. Nothing is selected on what it is scored on, and the move happens once.

Throughout your training and test splits you will see that the p-value can rise and fall with increasing coverage. Dataset coverage sharpens our understanding. If it is bouncing around, the results are likely insignificant, on the other hand insignificant results can turn out to show real impact and vice versa. Overall more data is usually better, but not always available. Many situations will require you to make decisions with only 70% of the understanding rather than 90%. It becomes important that you differentiate between those that are easily reversible with high upside and those that might need further investigation before forming a recommendation.

Three archetype trajectories drawn against accumulating out-of-fold coverage, with no numbers on either axis. The vertical axis is the p-value, the horizontal axis is coverage as more of the data gets scored. A dashed horizontal line marks p equals 0.05, with the hatched region below it marked significant. Every trajectory starts wide and jittery at low coverage. One grey trajectory settles high above the line, a null that holds. A second grey trajectory starts below the line and rises back above it, a lucky split whose significance fades. A gold trajectory settles below the line, a real effect that crystallises. A note at the right reads one split cannot tell them apart; coverage can: it sharpens whatever is true, and it manufactures nothing.
Which of these is behind your one good split? One split cannot say. Score more folds and whatever is true settles. Schematic; the measured trajectories are in the paragraph above.
Why the resampling is over folds

The obvious way to put an interval on an out-of-fold score is to bootstrap the rows: resample the predictions with replacement a few thousand times and read off the spread. It is wrong here, and Kipu's own evaluation pipeline states why: out-of-fold rows within a fold share a fitted model, so they are not independent draws, and resampling them as if they were inflates the apparent sample size and therefore the significance. Row-level tests on paired per-row losses show the same failure, returning vanishingly small p-values for comparisons that a paired test at the fold level shows to have little if no gap at all.

The fix is to resample whole folds, so the unit of resampling is the unit that carries an independently fitted model, and to use a variance-corrected paired test built for repeated cross-validation, where the training sets of different folds overlap. Holm-Bonferroni then runs across the whole family of comparisons rather than each one alone. If a vendor cannot tell you what their resampling unit was, the metrics they are quoting needs to be taken with a grain of salt.

A successful and durable pilot looks roughly like this: enough positives that fold-to-fold variance is smaller than the effect you want, out-of-fold prediction rather than one single split, and a correction that accounts for luck and the number of tests.

Published evidence for further reading

Result

Breast tumour detection, MedMNIST ultrasound, 702 samples, images reduced classically to 202 variables and then to 156, IBM ibm_kingston hardware at 156 qubits, support-vector classifier on the 156 SHAP-selected variables of which 109 are quantum-derived

Reported

AUC 0.937 on the benchmark's predefined test split

Baseline compared against

  • the same classifier on the classical variables alone 0.887
  • Google AutoML Vision 0.919
  • ResNet-18 0.891
  • AutoKeras 0.871
  • ResNet-50 0.866
  • Auto-sklearn 0.836

Result

Molecular toxicity classification, 171 molecules by 156 descriptors, one qubit per descriptor on the same 156-qubit hardware, gradient boosting under 5x5 repeated stratified cross-validation, 25 test sets

Reported

+121% in precision for the combined classical plus quantum feature set, relative, with no absolute precision published either side

Baseline compared against

  • the same gradient boosting on the 156 classical descriptors alone, which is how the paper's abstract and conclusion frame the comparison
  • no third-party model reported for this dataset
  • no absolute precision published either side

Result

Aerial tree-genus classification, balanced five-class TreeSatAI subset, 1,000 training and 200 test images, random forest on features reduced to each backend's qubit width, three IBM Heron processors plus the ibm_aer simulator

Reported

best 87.0%, from the quantum-only feature set on ibm_pittsburgh run at 120 of its 156 qubits; the hybrid set reached 86.5% on both 120-qubit backends and 81.5% on the 156-qubit one

Baseline compared against

  • ResNet-50 features plus a random forest at 84.0% on 120 variables, the strongest classical configuration
  • bare ResNet-50 83.0% on 120
  • both fall to 79.5% at the matched 156-variable width

Result

Motor-imagery EEG classification, 51 participants, 600 task cues each over three sessions, four-second 58-lead waveforms reduced by a spatial-temporal deep-learning stage to a 24-dimension feature set, no hardware named

Reported

88.8% accuracy (533 of 600) and AUROC 0.962 on the one held-out participant; 89.8% (27,478 of 30,590) and AUROC 0.970 across all 51, a figure that includes the 50 training participants and is therefore not out of sample

Baseline compared against

  • none published; the deep-learning stage feeds the final classifier rather than acting as a comparator, and no quantum-free ablation is reported
  • the 89.8% across all 51 includes the 50 training participants and is therefore not out of sample

Result

Off-line quantum feature extraction on the TreeSatAI tree-genus benchmark, the same 1,000 training and 200 test images as the card above, the quantum step run once on a 200-image subsample, one fifth of the training set, on an IBM processor at 120 active qubits of 156, a ridge surrogate then reproducing those features for the remaining data at classical cost

Reported

87% accuracy for the ridge surrogate, matching the 87.0% of the full quantum pipeline, at a fivefold reduction in the number of quantum executions

Baseline compared against

  • the full quantum pipeline, DQFE at 120 qubits, 87.0%
  • the classical ResNet-50 plus random forest configuration at 120 features, 84.0%
  • both figures are carried over from arXiv:2602.18350 rather than re-run for this paper
  • the 87.0% and 84.0% are the same measurements as the tree-genus card, not an independent confirmation of them; what is new here is the cost, not the accuracy
  • the surrogate figure is published as 87%, without a decimal place

Result

Toxic-activity classification on the UCI Toxicity dataset, the top 200 of its 1,203 chemical descriptors kept by mutual information and one qubit used per retained descriptor, features read from a 20-nanosecond quench on D-Wave Advantage2-prototype2.6, gradient boosting then trained on the quantum-enhanced features. This is the annealer route behind the Huk service, not the IBM feature extraction the rest of this page describes

Reported

median AUC 0.88, accuracy 0.75, precision 0.78, recall 0.52 and F1 0.77 for gradient boosting on the quantum-enhanced features, medians of the 50 scores from a stratified 10-fold, 5-repeat cross-validation

Baseline compared against

  • the same gradient boosting on the original features, AUC 0.62, accuracy 0.54, precision 0.40, recall 0.33, F1 0.53
  • random forest 0.63, support-vector machine 0.62 and XGBoost 0.56 on AUC, all on the original features
  • no third-party model reported for this dataset
  • the paper does not state how many molecules the dataset holds, only that the test fold is around 40 samples
  • a different platform from the IBM results above, so the two are not evidence for each other

Result

Prediction of atrial-fibrillation episodes on the UCI myocardial infarction complications dataset, 111 qubits, the same quench feature map on D-Wave Advantage2-prototype2.6 and the same stratified 10-fold, 5-repeat cross-validation, gradient boosting trained on the quantum-enhanced features

Reported

median AUC 0.78 and accuracy 0.63 for gradient boosting on the quantum-enhanced features; the support-vector variant reaches accuracy 0.70 with AUC 0.77

Baseline compared against

  • the same gradient boosting on the original features, AUC 0.70, accuracy 0.62
  • random forest 0.71, XGBoost 0.68 and support-vector machine 0.46 on AUC, all on the original features
  • no third-party model reported for this dataset
  • not every metric moves the same way: recall falls from 0.41 to 0.29 for gradient boosting on this task
  • the margin here is far smaller than on the toxicity dataset in the same paper

Open a card for its full result, its baselines and its source.

Feature width is the axis that separates our examples from the published wins. Our runs sit at nine to fifteen features; the published wins sit at 120 to 156 with similar sample counts. Rimay Simulator (free, hub.kipu-quantum.com) reaches 15 features and Rimay DQFE 133 to 156 features on IBM Heron processors with up to 20,000 samples. Those can then train surrogate models that help extract the gains and offer performance improvements on classical models and on premise. A functional map of features on a Heron QPU might look like this:

Coupling map of an IBM heavy-hexagonal quantum processor. Around a hundred and fifty small circles, each a numbered qubit, are joined by short links into a lattice of hexagons tiled edge to edge across a diagonal band. A label strip beneath the lattice reads 156 nodes, 120 active qubits, heavy-hex topology.
TreeSatAI benchmark: Each circle represents one qubit; 120 of the 156 are active, one feature each. Figure 2 of Flores-Garrigos, C. et al., arXiv:2605.19801.

Two boundaries of this session. The first is the hardware. Presented results are all based on the free Rimay simulator from hub.kipu-quantum.com, and never a quantum processor. The published wins above come from research pipelines. The feature count limitation also applies to quantum computers, at a different scale: up to 156 features, and the measurements also depend on the hardware architecture. An example of how features might be distributed on the chip is shown above.

The second limitation is publicly available benchmarks. Due to the sensitivity of its applications, there are few published Rimay numbers on enterprise business records: like churn, payment delay, supplier risk. If your question is what this does on your dataset an experiment might have to be run, before coming to a conclusion. The simulator can give hints, but is not the full story.

A horizontal feature-count axis. A marked point near the left end sits at 15 features, free simulator, labelled everything measured here. A hatched band far to the right, labelled where the published wins live, covers 120 to 156 features on hardware. The stretch between the two carries the note little structure here to expose.

The surrogate, and why serving stays classical

Regarding the promise of classical inference made in section two. Kipu's off-line framework runs the quantum processor over a representative sample of the training data during training. It is important that its distribution represents the full set, to then train a lightweight classical surrogate on the quantum data transformations. Once extracted, everything after that is classical: no quantum hardware queue at prediction time. Along with all its benefits: ordinary inference latency, ordinary MLOps cadence, on premise capability.

The Path B architecture extended by the surrogate, drawn as two lanes. The training lane runs once, on a schedule: the training data goes through the quantum pipeline and a classical surrogate learns the transformation from what comes back, then is handed down to the inference lane. The inference lane runs on every prediction and is marked classical prediction, no quantum step: data flows through the surrogate into your unchanged model and out as a prediction, with no quantum processor and no hardware queue anywhere in the path. A large callout reads runs classically on-premise: the surrogate is a classical model, deployed on your servers, preserving up to 99 percent of the fully quantum pipeline's accuracy.
The architecture from section 2, finished, and the point is the end product: the surrogate is a classical model you deploy and run in production yourself. The quantum pipeline appears once, in the training lane. The retention figure is the paper's breast-ultrasound case, AUC 0.932 against the full pipeline's 0.937 (arXiv:2605.19801, preprint); the satellite case in the text is reported separately by the paper itself.

On the published satellite case, the classical surrogate reached the same 87 percent accuracy as the fully quantum pipeline with fivefold fewer quantum executions, from a 200-of-1,000 subsample (arXiv:2605.19801, preprint). Separately, on the breast ultrasound case the surrogate retained AUC 0.932 against the full pipeline's 0.937. Not all improvements can be captured at all times. What can go wrong during the subsample selection is creating bias or sparse, non-representative data that extrapolates badly. Inputs falling outside the covered region need to be flagged and accounted for.

Conclusion

Rimay offers an interesting choice for models that have stopped improving despite existing tuning efforts. At one additional tool-call it offers a lightweight and cheap addition to Machine Learning pipelines. As a ready-made service it allows Machine Learning practitioners easy and high-velocity experiments with quantum computing. Rimay's product page covers the service boundaries and details. If you are looking for optimization, session two, quantum optimization, covers the evidence for Kipu's quantum optimizers. The full curriculum is available via the Academy page.

quan·tum ma·chine learn·ing

/ˈkwɒn.təm məˈʃiːn ˈlɜː.nɪŋ/noun

  1. 1

    the use of a quantum processor somewhere inside a machine-learning pipeline, whether to represent the data, to fit the model, or to make the prediction.

  2. 2

    Rimay is the tool that computes the re-presented fit-transform structure of your data so that a model can act on it.

How was this session?

Email is optional. If provided, it’s only used to follow up on your feedback.

Would it move your model?

Bring a table that has gone flat: the positives you have, the columns you feed, and the model you would keep. Get in touch and we will score whether the experiment is worth running.