This paper presents a formal validation framework establishing equivalence between Quadratic Unconstrained Binary Optimization (QUBO) formulations and classical minimal cut set (MCS) identification for fault tree analysis in probabilistic risk assessment. The framework develops a 21-gate hierarchical verification system that proves QUBO ground states correspond exactly to minimal cut sets computed by the SCRAM fault tree solver. The validation methodology employs exhaustive enumeration for tractable models and conservative spot checks for larger instances, achieving PASS status across all verification gates. The theoretical contribution formalizes MCS identification as a constrained Boolean satisfiability problem amenable to quantum optimization. Fault trees with OR-of-ANDs structure are encoded into QUBO form with penalty terms enforcing logic constraints. The encoding preserves semantic equivalence: valid QUBO ground states map bijectively to classical minimal cut sets. Empirical validation demonstrates the framework on a canonical test case (TOP = A ∨ B ∨ (C ∧ D)) using QAOA simulation with plateau-based convergence methodology. The quantum circuit correctly identifies valid MCS configurations with p_opt = 0.01058 (95% CI: [0.00972, 0.01151]). Hardware feasibility is demonstrated through execution on IBM's ibm_torino 133-qubit Heron processor, where all circuits passed sanity and interpretation checks despite noise distributing probability across 110-120 unique states. The framework establishes foundational infrastructure for quantum PRA research by providing: (1) formal QUBO-MCS bijection verification against ground truth, (2) a reproducible classical baseline using the SCRAM solver, (3) plateau-based convergence criteria for stochastic quantum optimization, and (4) a hierarchical verification gate system ensuring scientific defensibility. No claim of quantum advantage is made; results are evaluated solely for semantic correctness, reproducibility, and NISQ feasibility.
This paper presents a formal validation framework establishing equivalence between Quadratic Unconstrained Binary Optimization (QUBO) formulations and classical minimal cut set (MCS) identification for fault tree analysis in probabilistic risk assessment. The framework develops a 21-gate hierarchical verification system that proves QUBO ground states correspond exactly to minimal cut sets computed by the SCRAM fault tree solver. The validation methodology employs exhaustive enumeration for tractable models and conservative spot checks for larger instances, achieving PASS status across all verification gates. The theoretical contribution formalizes MCS identification as a constrained Boolean satisfiability problem amenable to quantum optimization. Fault trees with OR-of-ANDs structure are encoded into QUBO form with penalty terms enforcing logic constraints. The encoding preserves semantic equivalence: valid QUBO ground states map bijectively to classical minimal cut sets. Empirical validation demonstrates the framework on a canonical test case (TOP = A ∨ B ∨ (C ∧ D)) using QAOA simulation with plateau-based convergence methodology. The quantum circuit correctly identifies valid MCS configurations with p_opt = 0.01058 (95% CI: [0.00972, 0.01151]). Hardware feasibility is demonstrated through execution on IBM's ibm_torino 133-qubit Heron processor, where all circuits passed sanity and interpretation checks despite noise distributing probability across 110-120 unique states. The framework establishes foundational infrastructure for quantum PRA research by providing: (1) formal QUBO-MCS bijection verification against ground truth, (2) a reproducible classical baseline using the SCRAM solver, (3) plateau-based convergence criteria for stochastic quantum optimization, and (4) a hierarchical verification gate system ensuring scientific defensibility. No claim of quantum advantage is made; results are evaluated solely for semantic correctness, reproducibility, and NISQ feasibility.
Quadratic Unconstrained Binary Optimization formulations for fault tree analysis encode gate logic through penalty terms, but prior work has demonstrated that these encodings fail to preserve fault tree semantics at scale: across 80 extended fault tree instances, only 7.5% of Rosenberg-based soft QUBO solutions match the constrained exact baseline, while 27.5% produce outright semantic violations. This paper develops two complementary methods to restore semantic correctness. A Hierarchical Semantic Validator and repair framework achieves 100% semantic validity across all 80 cohort instances through structure-aware bit-flip correction, improving the exact match rate from 7.5% to 28.75%. A constraint-level QUBO encoding eliminates penalty-based failures entirely by restricting optimization to the feasible subspace, reducing qubit requirements by a factor of 1.85x relative to Rosenberg encodings. Validation through Qiskit mixer circuits implementing a feasibility-preserving Quantum Alternating Operator Ansatz confirms zero probability leakage onto infeasible states for representative instances. These methods establish that semantic-preserving quantum optimization for fault tree minimal cut set identification is achievable with current algorithmic tools. Complete experimental artifacts with SHA-256 verification are provided under the DOI above.
Minimal cut set (MCS) identification is central to fault-tree-based probabilistic risk assessment (PRA), but quantum MCS methods have largely remained simulation-based or require Grover-style depths that are not executable on current hardware. This study reports an end-to-end IBM Quantum validation using shallow, fixed-parameter alternating-operator constraint-level QUBO (CL-QUBO) sampling circuits. The primary campaign executed 132 full-paired raw RL0 cells across 66 SCRAM-derived reactor fault-tree subtrees on two IBM devices, with 12 mapping-replicate controls held separate. Under an independently locked little-endian decoding convention, top-200 hardware outputs produced 2,596 exact-MCS state matches against 778 expected under a matched random-projection baseline, a 3.34× enrichment with 0.901 aggregate unique-MCS recovery. A 2,000-trial matched random-state null placed the all-campaign observed count outside the null distribution (empirical p < 5.0 × 10⁻⁴). Recovered cut sets preserved local PRA importance rankings with high fidelity, with mean Spearman correlations of 0.969 for Fussell-Vesely, 0.950 for Risk Reduction Worth, 0.934 for Birnbaum, and 0.920 for Risk Achievement Worth. A post-primary, pre-registered attribution ladder showed that removing the CL-QUBO cost phase collapsed the hardware result to the random/no-cost simulation scale (top-200 ratio 0.985), while the full p = 1 circuit remained at 3.34×. The results support cost-phase attribution for exact-MCS enrichment on current hardware, while making no claim of quantum advantage, plant-level risk quantification, or uniform behavior across all subtree classes.
This paper presents a reproducible, audit-ready quantum hardware experiment for minimal cut set (MCS) inference on reactor fault tree subtrees using IBM Quantum hardware. Eight two-qubit subtrees were selected from a Phase 2B reactor-scale fault tree corpus spanning four nuclear hazard families (EQK, LOCA, TRANS, SGTR). Classical baselines were computed with SCRAM 0.16.2, yielding two order-1 MCS per instance. Shallow p = 1 Quantum Alternating Operator Ansatz (QAOA+) circuits were constructed using the Constraint-Level QUBO (CL-QUBO) encoding from Peters (2026b), transpiled to ibm_torino, and executed as batch jobs totaling 57,344 shots across four experiments. Raw hardware results showed apparent infeasible leakage of 13.2% on the hw = 0 state. A dedicated four-circuit readout calibration protocol revealed that this leakage was almost entirely attributable to readout error: after confusion-matrix correction (κ = 1.56), the baseline experiment concentrates 98.3% of probability mass on the MCS manifold (hw = 1), with residual infeasible mass of 0.5%. Extension experiments substantially broaden the empirical foundation. Cross-pair validation on three additional qubit pairs (ibm_torino pairs 61–62 and 36–48; ibm_marrakesh pair 6–7) confirms that corrected MCS fidelity exceeds 95% across all tested physical qubit assignments (range 95.2%–98.3%), with the best pairs achieving 98.2% (ibm_marrakesh 6–7) and 97.9% (ibm_torino 61–62). A QAOA depth sweep (p ∈ {1, 2, 3}) on ibm_torino pair 0–1 quantifies noise–depth tradeoffs: because the mixer is inactive (β = 0), the observed redistribution from hw = 2 toward hw = 1 (7.8% at p = 1, 25.6% at p = 2, 40.9% at p = 3) is attributable to gate noise rather than productive mixing, with infeasible hw = 0 mass rising concurrently from 0.2% to 7.2%. A three-qubit exactly-one constraint experiment on ibm_torino triple 0–1–2 achieves 92.2% corrected hw = 1 mass from a feasible initialization, demonstrating that the CL-QUBO framework extends beyond two-qubit OR gates to richer constraint types on quantum hardware. Hard-start initialization experiments confirm that p = 1 circuits do not induce sufficient mixing to recover MCS states from non-minimal initializations, with this conclusion stable across independent replication and all tested qubit pairs. Complete experimental artifacts with SHA-256 verification are provided under the DOI above.
Quadratic Unconstrained Binary Optimization (QUBO) formulations have been proposed as a pathway for applying quantum optimization algorithms to fault tree analysis. These formulations encode fault tree logic through penalty terms, converting constrained minimal cut set identification into unconstrained optimization. Prior work has demonstrated QUBO-MCS equivalence for small, carefully structured instances. This paper investigates whether penalty-based encodings of the type commonly used in the literature reliably preserve fault tree semantics at realistic problem scales when solved exactly. A cohort of 80 extended fault tree instances with OR, AND, and NOT gates was evaluated using two exact solvers: a constrained baseline that enforces fault tree semantics directly, and a branch-and-bound solver operating on soft QUBO formulations using Rosenberg product reductions. Even with the TOP event variable fixed, penalty coefficients conservatively scaled to exceed the maximum objective contribution, and exact (not heuristic) optimization, the soft QUBO solutions matched the constrained baseline in only 6 of 80 instances (7.5%). More critically, 22 of 80 instances (27.5%) produced solutions that violate fault tree gate semantics entirely. These results demonstrate that Rosenberg-based soft QUBO formulations with dominance-safe penalties do not reliably preserve fault tree semantics at scale under exact optimization, indicating a structural limitation of this class of unconstrained penalty encodings rather than a tuning deficiency. Complete experimental artifacts with SHA-256 verification are provided under the DOI above.
This paper presents the first end-to-end evaluation of whether quantum-identified minimal cut sets (MCS) from physical IBM Quantum hardware can preserve standard probabilistic risk assessment (PRA) importance measures when compared against frozen SCRAM-validated classical baselines. It consumes frozen CL-QUBO hardware outputs from Paper 10 and propagates recovered MCS through two defensible probability-assignment tracks: a primary direct-component track based on NUREG/CR-6928 original component unreliability data and a secondary CCF extension track based on the 2020 NRC/INL common cause failure update using alpha factor methodology. Of 21 selected representative subtrees, 16 are promoted into the formal quantitative analysis set (10 primary, 6 extension), while five additional subtrees remain executed but quarantined in exploratory holdout lanes. Exact MCS recovery remains incomplete across both tracks. At the shared controlled operating point (ibm_marrakesh RL1 top-k = 10), the primary track achieves mean recovery fraction 0.225 with mean top-event absolute error 4.52 × 10⁻⁷, while the extension track achieves mean recovery fraction 0.292 with mean RAW Spearman 0.488, Birnbaum Spearman 0.438, FV Spearman 0.236, and mean top-event absolute error 7.50 × 10⁻⁷. Across the promoted set, tighter extraction windows generally outperform wider windows, set Jaccard is more informative than simple recovery fraction for downstream utility within this study, and the extension-track advantage survives common-support restriction but remains conditional and metric-dependent. The resulting pipeline therefore shows measurable but limited PRA utility: enough for screening-level assessment of relative component significance under favorable conditions, but not enough for standalone regulatory-grade quantification.
This study presents a bounded, auditable quantum execution framework for probabilistic risk assessment (PRA) using OpenPRA-derived fault tree encodings and provides the first platform-integrated empirical characterization of structural scaling limits within the OpenPRA quantum execution workflow. A tri-provider architecture integrating local structural validation, IBM gate-based quantum hardware, and a D-Wave-ready interface is applied to 132 constraint-level QUBO (CL-QUBO) export artifacts retained from the OpenPRA quantum-feasible candidate pool. A systematic scaling gap analysis reveals that all 132 cases occupy exactly two qubit levels—5 qubits (23 cases, 17.4%) and 8 qubits (109 cases, 82.6%)—with zero semantically valid same-root expansions beyond 8 qubits. IBM Quantum hardware validation on ibm_marrakesh confirms bounded MCS recovery at both qubit levels, including exact hardware recovery (3/3 reference MCS) at 5 qubits and partial recovery (3/4 reference MCS) at 8 qubits, with job-level provenance for every execution. A hybrid filtering framework generates 4,692 synthetic cross-root expansion candidates, of which zero preserve semantic PRA validity; 14 selected 16-qubit candidates are validated through a deterministic structural execution validator, and three are executed on IBM hardware as stress tests (156 physical qubits after transpilation, 205–220 unique bitstrings observed per job). These results establish strong evidence that the observed 8-qubit ceiling is a property of the current CL-QUBO export pipeline applied to real reactor-scale fault trees, that scaling beyond this ceiling requires structural decomposition or alternative execution paradigms, and that the integrated OpenPRA quantum platform can execute, validate, and bundle evidence across the full pipeline with auditable provenance. No claim of quantum advantage or production PRA quantification is made.