H₂ Ground State Energy via VQE
Variational quantum chemistry on 2 qubits -- from simulation to hardware
We calculate the ground state energy of molecular hydrogen using a 2-qubit VQE on four quantum processors. Chemical accuracy achieved on IBM Torino (0.22 kcal/mol with TREX) and QI Tuna-9 (0.92 kcal/mol with hybrid PS+REM). IQM Garnet reaches 14.26 kcal/mol after fixing a platform-specific CNOT decomposition error. The emulator nails chemical accuracy across 14 bond distances.
Research Question
Can a 2-qubit variational quantum eigensolver (VQE) calculate the ground state energy of molecular hydrogen within chemical accuracy (1.6 milliHartree), and how does hardware noise affect the result?
Prior Work
The Variational Quantum Eigensolver (VQE) was proposed by Peruzzo et al. in 2014 as a hybrid quantum-classical algorithm for finding molecular ground states. It is considered one of the most promising near-term applications of quantum computing because it can tolerate some hardware noise through its variational nature.
For molecular hydrogen (H<sub>2</sub>), the electronic structure problem can be mapped to qubits using the Jordan-Wigner or Bravyi-Kitaev transformation. The minimal STO-3G basis set gives a 4-qubit Hamiltonian, but Bravyi-Kitaev transformation with qubit tapering (exploiting Z<sub>2</sub> symmetries) reduces this to just 2 qubits -- the minimum for a non-trivial quantum chemistry calculation.
Sagastizabal et al. (2019) demonstrated VQE on a superconducting transmon processor, achieving chemical accuracy at the equilibrium bond distance. Our experiment replicates and extends this work across three platforms, including a full 14-point dissociation curve on the emulator.
The gold standard is chemical accuracy: getting within 1.6 milliHartree (1 kcal/mol) of the exact Full Configuration Interaction (FCI) energy. At the equilibrium bond distance of 0.735 Å, the FCI energy is -1.1373 Hartree.
Method
We use the 2-qubit Hamiltonian obtained by Bravyi-Kitaev transformation with qubit tapering of the H<sub>2</sub> STO-3G Hamiltonian in the correct parity sector. The Hamiltonian has the form:
H = g<sub>0</sub>I + g<sub>1</sub>Z<sub>0</sub> + g<sub>2</sub>Z<sub>1</sub> + g<sub>3</sub>Z<sub>0</sub>Z<sub>1</sub> + g<sub>4</sub>X<sub>0</sub>X<sub>1</sub> + g<sub>5</sub>Y<sub>0</sub>Y<sub>1</sub>
Ansatz: R<sub>Y</sub>(θ) on qubit 0, then CNOT(q[0], q[1]). The single parameter θ is optimized classically.
Measurement protocol: Three separate circuits measure in the Z, X, and Y bases to reconstruct all expectation values. 4096+ shots per basis per backend.
Dissociation curve: The emulator sweeps 14 bond distances from 0.3 to 3.0 Å with 65,536 shots each, using pre-computed optimal θ values at each distance.
Results
Platform Comparison
| Backend | Type | Key Metric | Date |
|---|---|---|---|
IBM Torino (133q) | Hardware | -1.1226 Ha (9.22 kcal/mol) | 2/11/2026 |
IBM Torino (133q) | Hardware | -1.1197 Ha (11 kcal/mol) | 2/10/2026 |
iqm-garnet | Emulator | -1.1146 Ha (14.26 kcal/mol) | 2/10/2026 |
QI Tuna-9 (9q) | Hardware | -- | 2/10/2026 |
QI Emulator | Emulator | -1.1385 Ha (0.75 kcal/mol) | 2/10/2026 |
QI Tuna-9 (9q) | Hardware | -1.1310 Ha (3.94 kcal/mol) | 2/10/2026 |
QI Tuna-9 (9q) | Hardware | -1.1293 Ha (5.01 kcal/mol) | 2/10/2026 |
QI Tuna-9 (9q) | Hardware | -1.1358 Ha (0.92 kcal/mol) | 2/10/2026 |
QI Tuna-9 (9q) | Hardware | -1.1213 Ha (10.02 kcal/mol) | 2/10/2026 |
QI Tuna-9 (9q) | Hardware | -1.1329 Ha (2.79 kcal/mol) | 2/10/2026 |
QI Tuna-9 (9q) | Hardware | -1.1275 Ha (6.14 kcal/mol) | 2/10/2026 |
QI Tuna-9 (9q) | Hardware | -1.1263 Ha (6.93 kcal/mol) | 2/10/2026 |
QI Tuna-9 (9q) | Hardware | -1.1160 Ha (13.4 kcal/mol) | 2/10/2026 |
QI Tuna-9 (9q) | Hardware | -1.1291 Ha (5.13 kcal/mol) | 2/10/2026 |
QI Tuna-9 (9q) | Hardware | -1.1294 Ha (4.94 kcal/mol) | 2/10/2026 |
QI Tuna-9 (9q) | Hardware | -1.1213 Ha (10.02 kcal/mol) | 2/10/2026 |
QI Tuna-9 (9q) | Hardware | -1.1221 Ha (9.53 kcal/mol) | 2/10/2026 |
QI Tuna-9 (9q) | Hardware | -1.1238 Ha (8.5 kcal/mol) | 2/10/2026 |
QI Tuna-9 (9q) | Hardware | -1.1230 Ha (9 kcal/mol) | 2/10/2026 |
QI Tuna-9 (9q) | Hardware | -1.1272 Ha (6.31 kcal/mol) | 2/10/2026 |
QI Emulator | Emulator | -1.1296 Ha (4.83 kcal/mol) | 2/10/2026 |
QI Emulator | Emulator | -1.1415 Ha (2.64 kcal/mol) | 2/10/2026 |
QI Tuna-9 (9q) | Hardware | -0.9147 Ha (13.42 kcal/mol) | 2/10/2026 |
QI Tuna-9 (9q) | Hardware | -1.1325 Ha (3.04 kcal/mol) | 2/10/2026 |
QI Emulator | Emulator | -0.9362 Ha (0.09 kcal/mol) | 2/10/2026 |
QI Emulator | Emulator | -0.9454 Ha (2.05 kcal/mol) | 2/10/2026 |
QI Emulator | Emulator | -0.9940 Ha (2.62 kcal/mol) | 2/10/2026 |
QI Emulator | Emulator | -1.1034 Ha (1.4 kcal/mol) | 2/10/2026 |
QI Emulator | Emulator | -1.1364 Ha (0.6 kcal/mol) | 2/10/2026 |
QI Emulator | Emulator | -1.0583 Ha (2 kcal/mol) | 2/10/2026 |
iqm-garnet | Emulator | -0.3105 Ha (518.85 kcal/mol) | 2/10/2026 |
QI Tuna-9 (9q) | Hardware | -0.9210 Ha (17.32 kcal/mol) | 2/10/2026 |
QI Tuna-9 (9q) | Hardware | -0.9779 Ha (12.68 kcal/mol) | 2/10/2026 |
QI Tuna-9 (9q) | Hardware | -1.0946 Ha (4.12 kcal/mol) | 2/10/2026 |
QI Tuna-9 (9q) | Hardware | -1.0393 Ha (9.98 kcal/mol) | 2/10/2026 |
QI Tuna-9 (9q) | Hardware | -1.1222 Ha (9.45 kcal/mol) | 2/10/2026 |
IBM Torino (133q) | Hardware | -1.0976 Ha (24.91 kcal/mol) | 2/10/2026 |
QI Tuna-9 (9q) | Hardware | -1.0045 Ha (83.3 kcal/mol) | 2/10/2026 |
IBM Marrakesh (156q) | Hardware | -1.0956 Ha (26.1505 kcal/mol) | 2/10/2026 |
ibm-torino | Hardware | -1.0408 Ha (60.54 kcal/mol) | Invalid Date |
QI Tuna-9 (9q) | Hardware | -- | Invalid Date |
H₂ Dissociation Curve
14 emulator points + 7 hardware pointsEnergy vs. bond distance for molecular hydrogen. The VQE emulator matches the exact (FCI) curve within chemical accuracy at all 14 distances. Purple diamonds show Tuna-9 hardware results with readout error mitigation (sector-projected 2-qubit ansatz, native CZ gate set). Error bars show ±1σ across 5 independent hardware runs.
Measured Energy
-1.1226 Ha
FCI Reference
-1.1373 Ha
Energy Level Diagram
Expectation Values
<Z0>
-0.961
<Z1>
0.950
<Z0Z1>
-0.969
<X0X1>
-0.256
<Y0Y1>
-0.197
Measurement Counts by Basis
VQE energy: -1.1226 Ha (FCI: -1.1373 Ha). Error: 9.22 kcal/mol. Not chemical accuracy — hardware noise degrades Z correlations by ~3-5%.
View raw JSONMeasured Energy
-1.1197 Ha
FCI Reference
-1.1373 Ha
Energy Level Diagram
Expectation Values
<Z0>
-0.956
<Z1>
0.949
<Z0Z1>
-0.971
<X0X1>
-0.218
<Y0Y1>
-0.229
IBM Torino H2 VQE: -1.1197 Ha, 11.0 kcal/mol error. Hardware noise degrades result but correct sector and qualitative behavior observed.
View raw JSONMeasured Energy
-1.1146 Ha
FCI Reference
-1.1373 Ha
Energy Level Diagram
Expectation Values
<Z0>
-0.949
<Z1>
0.949
<Z0Z1>
-1.000
<X0X1>
-0.206
<Y0Y1>
-0.196
Measurement Counts by Basis
Measurement Counts by Basis
Measured Energy
-1.1385 Ha
FCI Reference
-1.1373 Ha
Energy Level Diagram
Expectation Values
<Z0>
-0.973
<Z1>
0.973
<Z0Z1>
-1.000
<X0X1>
-0.252
<Y0Y1>
-0.219
Measurement Counts by Basis
VQE energy: -1.1385 Ha (FCI: -1.1373 Ha). Error: 0.75 kcal/mol. Within chemical accuracy.
View cQASM circuit
version 3.0 qubit[2] q bit[2] b // Subspace-preserving ansatz: Ry-CNOT-X // State = cos(a/2)|10> + sin(a/2)|01> Ry(-0.223539) q[0] CNOT q[0], q[1] X q[0] // Z-basis measurement b = measure q
Measured Energy
-1.1310 Ha
FCI Reference
-1.1373 Ha
Energy Level Diagram
Expectation Values
<Z0>
-0.963
<Z1>
0.963
<Z0Z1>
-1.000
<X0X1>
-0.209
<Y0Y1>
-0.200
Measurement Counts by Basis
VQE energy: -1.1310 Ha (FCI: -1.1373 Ha). Error: 3.9 kcal/mol. Outside chemical accuracy. Post-selection kept 83% of Z-basis shots. Hybrid(PS+REM): -1.1310 Ha (3.9 kcal/mol). Full-REM: -1.0751 Ha (39.0 kcal/mol).
View cQASM circuit
version 3.0 qubit[5] q bit[5] b // Subspace-preserving ansatz: Ry-CNOT-X (cnot_folds=1) // State = cos(a/2)|10> + sin(a/2)|01>, alpha=-0.223400 Ry(-0.223400) q[2] CNOT q[2], q[4] X q[2] // Z-basis measurement b = measure q
Measured Energy
-1.1293 Ha
FCI Reference
-1.1373 Ha
Energy Level Diagram
Expectation Values
<Z0>
-0.964
<Z1>
0.964
<Z0Z1>
-1.000
<X0X1>
-0.202
<Y0Y1>
-0.186
Measurement Counts by Basis
VQE energy: -1.1293 Ha (FCI: -1.1373 Ha). Error: 5.0 kcal/mol. Outside chemical accuracy. Post-selection kept 91% of Z-basis shots. Hybrid(PS+REM): -1.1293 Ha (5.0 kcal/mol). Full-REM: -1.1277 Ha (6.0 kcal/mol).
View cQASM circuit
version 3.0 qubit[5] q bit[5] b // Subspace-preserving ansatz: Ry-CNOT-X (cnot_folds=1) // State = cos(a/2)|10> + sin(a/2)|01>, alpha=-0.223400 Ry(-0.223400) q[2] CNOT q[2], q[4] X q[2] // Z-basis measurement b = measure q
Measured Energy
-1.1358 Ha
FCI Reference
-1.1373 Ha
Energy Level Diagram
Expectation Values
<Z0>
-0.969
<Z1>
0.969
<Z0Z1>
-1.000
<X0X1>
-0.210
<Y0Y1>
-0.207
Measurement Counts by Basis
VQE energy: -1.1358 Ha (FCI: -1.1373 Ha). Error: 0.9 kcal/mol. Within chemical accuracy. Post-selection kept 96% of Z-basis shots. Hybrid(PS+REM): -1.1358 Ha (0.9 kcal/mol). Full-REM: -1.1325 Ha (3.0 kcal/mol).
View cQASM circuit
version 3.0 qubit[5] q bit[5] b // Subspace-preserving ansatz: Ry-CNOT-X (cnot_folds=1) // State = cos(a/2)|10> + sin(a/2)|01>, alpha=-0.223400 Ry(-0.223400) q[2] CNOT q[2], q[4] X q[2] // Z-basis measurement b = measure q
Measured Energy
-1.1213 Ha
FCI Reference
-1.1373 Ha
Energy Level Diagram
Expectation Values
<Z0>
-0.960
<Z1>
0.960
<Z0Z1>
-1.000
<X0X1>
-0.217
<Y0Y1>
-0.185
Measurement Counts by Basis
VQE energy: -1.1213 Ha (FCI: -1.1373 Ha). Error: 10.0 kcal/mol. Outside chemical accuracy. Post-selection kept 93% of Z-basis shots.
View cQASM circuit
version 3.0 qubit[5] q bit[5] b // Subspace-preserving ansatz: Ry-CNOT-X (cnot_folds=5) // State = cos(a/2)|10> + sin(a/2)|01>, alpha=-0.223400 Ry(-0.223400) q[2] CNOT q[2], q[4] CNOT q[2], q[4] CNOT q[2], q[4] CNOT q[2], q[4] CNOT q[2], q[4] X q[2] // Z-basis measurement b = measure q
Measured Energy
-1.1329 Ha
FCI Reference
-1.1373 Ha
Energy Level Diagram
Expectation Values
<Z0>
-0.976
<Z1>
0.976
<Z0Z1>
-1.000
<X0X1>
-0.198
<Y0Y1>
-0.186
Measurement Counts by Basis
VQE energy: -1.1329 Ha (FCI: -1.1373 Ha). Error: 2.8 kcal/mol. Outside chemical accuracy. Post-selection kept 94% of Z-basis shots.
View cQASM circuit
version 3.0 qubit[5] q bit[5] b // Subspace-preserving ansatz: Ry-CNOT-X (cnot_folds=5) // State = cos(a/2)|10> + sin(a/2)|01>, alpha=-0.223400 Ry(-0.223400) q[2] CNOT q[2], q[4] CNOT q[2], q[4] CNOT q[2], q[4] CNOT q[2], q[4] CNOT q[2], q[4] X q[2] // Z-basis measurement b = measure q
Measured Energy
-1.1275 Ha
FCI Reference
-1.1373 Ha
Energy Level Diagram
Expectation Values
<Z0>
-0.972
<Z1>
0.972
<Z0Z1>
-1.000
<X0X1>
-0.193
<Y0Y1>
-0.171
Measurement Counts by Basis
VQE energy: -1.1275 Ha (FCI: -1.1373 Ha). Error: 6.1 kcal/mol. Outside chemical accuracy. Post-selection kept 95% of Z-basis shots.
View cQASM circuit
version 3.0 qubit[5] q bit[5] b // Subspace-preserving ansatz: Ry-CNOT-X (cnot_folds=5) // State = cos(a/2)|10> + sin(a/2)|01>, alpha=-0.223400 Ry(-0.223400) q[2] CNOT q[2], q[4] CNOT q[2], q[4] CNOT q[2], q[4] CNOT q[2], q[4] CNOT q[2], q[4] X q[2] // Z-basis measurement b = measure q
Measured Energy
-1.1263 Ha
FCI Reference
-1.1373 Ha
Energy Level Diagram
Expectation Values
<Z0>
-0.968
<Z1>
0.968
<Z0Z1>
-1.000
<X0X1>
-0.183
<Y0Y1>
-0.204
Measurement Counts by Basis
VQE energy: -1.1263 Ha (FCI: -1.1373 Ha). Error: 6.9 kcal/mol. Outside chemical accuracy. Post-selection kept 95% of Z-basis shots.
View cQASM circuit
version 3.0 qubit[5] q bit[5] b // Subspace-preserving ansatz: Ry-CNOT-X (cnot_folds=3) // State = cos(a/2)|10> + sin(a/2)|01>, alpha=-0.223400 Ry(-0.223400) q[2] CNOT q[2], q[4] CNOT q[2], q[4] CNOT q[2], q[4] X q[2] // Z-basis measurement b = measure q
Measured Energy
-1.1160 Ha
FCI Reference
-1.1373 Ha
Energy Level Diagram
Expectation Values
<Z0>
-0.964
<Z1>
0.964
<Z0Z1>
-1.000
<X0X1>
-0.196
<Y0Y1>
-0.110
Measurement Counts by Basis
VQE energy: -1.1160 Ha (FCI: -1.1373 Ha). Error: 13.4 kcal/mol. Outside chemical accuracy. Post-selection kept 96% of Z-basis shots.
View cQASM circuit
version 3.0 qubit[5] q bit[5] b // Subspace-preserving ansatz: Ry-CNOT-X (cnot_folds=3) // State = cos(a/2)|10> + sin(a/2)|01>, alpha=-0.223400 Ry(-0.223400) q[2] CNOT q[2], q[4] CNOT q[2], q[4] CNOT q[2], q[4] X q[2] // Z-basis measurement b = measure q
Measured Energy
-1.1291 Ha
FCI Reference
-1.1373 Ha
Energy Level Diagram
Expectation Values
<Z0>
-0.967
<Z1>
0.967
<Z0Z1>
-1.000
<X0X1>
-0.231
<Y0Y1>
-0.190
Measurement Counts by Basis
VQE energy: -1.1291 Ha (FCI: -1.1373 Ha). Error: 5.1 kcal/mol. Outside chemical accuracy. Post-selection kept 96% of Z-basis shots.
View cQASM circuit
version 3.0 qubit[5] q bit[5] b // Subspace-preserving ansatz: Ry-CNOT-X (cnot_folds=3) // State = cos(a/2)|10> + sin(a/2)|01>, alpha=-0.223400 Ry(-0.223400) q[2] CNOT q[2], q[4] CNOT q[2], q[4] CNOT q[2], q[4] X q[2] // Z-basis measurement b = measure q
Measured Energy
-1.1294 Ha
FCI Reference
-1.1373 Ha
Energy Level Diagram
Expectation Values
<Z0>
-0.967
<Z1>
0.967
<Z0Z1>
-1.000
<X0X1>
-0.217
<Y0Y1>
-0.211
Measurement Counts by Basis
VQE energy: -1.1294 Ha (FCI: -1.1373 Ha). Error: 4.9 kcal/mol. Outside chemical accuracy. Post-selection kept 95% of Z-basis shots.
View cQASM circuit
version 3.0 qubit[5] q bit[5] b // Subspace-preserving ansatz: Ry-CNOT-X (cnot_folds=1) // State = cos(a/2)|10> + sin(a/2)|01>, alpha=-0.223400 Ry(-0.223400) q[2] CNOT q[2], q[4] X q[2] // Z-basis measurement b = measure q
Measured Energy
-1.1213 Ha
FCI Reference
-1.1373 Ha
Energy Level Diagram
Expectation Values
<Z0>
-0.963
<Z1>
0.963
<Z0Z1>
-1.000
<X0X1>
-0.149
<Y0Y1>
-0.223
Measurement Counts by Basis
VQE energy: -1.1213 Ha (FCI: -1.1373 Ha). Error: 10.0 kcal/mol. Outside chemical accuracy. Post-selection kept 96% of Z-basis shots.
View cQASM circuit
version 3.0 qubit[5] q bit[5] b // Subspace-preserving ansatz: Ry-CNOT-X (cnot_folds=1) // State = cos(a/2)|10> + sin(a/2)|01>, alpha=-0.223400 Ry(-0.223400) q[2] CNOT q[2], q[4] X q[2] // Z-basis measurement b = measure q
Measured Energy
-1.1221 Ha
FCI Reference
-1.1373 Ha
Energy Level Diagram
Expectation Values
<Z0>
-0.967
<Z1>
0.967
<Z0Z1>
-1.000
<X0X1>
-0.191
<Y0Y1>
-0.154
Measurement Counts by Basis
VQE energy: -1.1221 Ha (FCI: -1.1373 Ha). Error: 9.5 kcal/mol. Outside chemical accuracy. Post-selection kept 95% of Z-basis shots.
View cQASM circuit
version 3.0 qubit[5] q bit[5] b // Subspace-preserving ansatz: Ry-CNOT-X (cnot_folds=1) // State = cos(a/2)|10> + sin(a/2)|01>, alpha=-0.223400 Ry(-0.223400) q[2] CNOT q[2], q[4] X q[2] // Z-basis measurement b = measure q
Measured Energy
-1.1238 Ha
FCI Reference
-1.1373 Ha
Energy Level Diagram
Expectation Values
<Z0>
-0.968
<Z1>
0.968
<Z0Z1>
-1.000
<X0X1>
-0.185
<Y0Y1>
-0.172
Measurement Counts by Basis
VQE energy: -1.1238 Ha (FCI: -1.1373 Ha). Error: 8.5 kcal/mol. Outside chemical accuracy. Post-selection kept 93% of Z-basis shots.
View cQASM circuit
version 3.0 qubit[5] q bit[5] b // Subspace-preserving ansatz: Ry-CNOT-X (cnot_folds=5) // State = cos(a/2)|10> + sin(a/2)|01>, alpha=-0.223400 Ry(-0.223400) q[2] CNOT q[2], q[4] CNOT q[2], q[4] CNOT q[2], q[4] CNOT q[2], q[4] CNOT q[2], q[4] X q[2] // Z-basis measurement b = measure q
Measured Energy
-1.1230 Ha
FCI Reference
-1.1373 Ha
Energy Level Diagram
Expectation Values
<Z0>
-0.964
<Z1>
0.964
<Z0Z1>
-1.000
<X0X1>
-0.182
<Y0Y1>
-0.202
Measurement Counts by Basis
VQE energy: -1.1230 Ha (FCI: -1.1373 Ha). Error: 9.0 kcal/mol. Outside chemical accuracy. Post-selection kept 95% of Z-basis shots.
View cQASM circuit
version 3.0 qubit[5] q bit[5] b // Subspace-preserving ansatz: Ry-CNOT-X (cnot_folds=3) // State = cos(a/2)|10> + sin(a/2)|01>, alpha=-0.223400 Ry(-0.223400) q[2] CNOT q[2], q[4] CNOT q[2], q[4] CNOT q[2], q[4] X q[2] // Z-basis measurement b = measure q
Measured Energy
-1.1272 Ha
FCI Reference
-1.1373 Ha
Energy Level Diagram
Expectation Values
<Z0>
-0.967
<Z1>
0.967
<Z0Z1>
-1.000
<X0X1>
-0.215
<Y0Y1>
-0.189
Measurement Counts by Basis
VQE energy: -1.1272 Ha (FCI: -1.1373 Ha). Error: 6.3 kcal/mol. Outside chemical accuracy. Post-selection kept 95% of Z-basis shots.
View cQASM circuit
version 3.0 qubit[5] q bit[5] b // Subspace-preserving ansatz: Ry-CNOT-X (cnot_folds=1) // State = cos(a/2)|10> + sin(a/2)|01>, alpha=-0.223400 Ry(-0.223400) q[2] CNOT q[2], q[4] X q[2] // Z-basis measurement b = measure q
Measured Energy
-1.1296 Ha
FCI Reference
-1.1373 Ha
Energy Level Diagram
Expectation Values
<Z0>
-0.968
<Z1>
0.968
<Z0Z1>
-1.000
<X0X1>
-0.203
<Y0Y1>
-0.220
Measurement Counts by Basis
VQE energy: -1.1296 Ha (FCI: -1.1373 Ha). Error: 4.8 kcal/mol. Outside chemical accuracy. Post-selection kept 100% of Z-basis shots.
View cQASM circuit
version 3.0 qubit[2] q bit[2] b // Subspace-preserving ansatz: Ry-CNOT-X (cnot_folds=5) // State = cos(a/2)|10> + sin(a/2)|01>, alpha=-0.223400 Ry(-0.223400) q[0] CNOT q[0], q[1] CNOT q[0], q[1] CNOT q[0], q[1] CNOT q[0], q[1] CNOT q[0], q[1] X q[0] // Z-basis measurement b = measure q
Measured Energy
-1.1415 Ha
FCI Reference
-1.1373 Ha
Energy Level Diagram
Expectation Values
<Z0>
-0.978
<Z1>
0.978
<Z0Z1>
-1.000
<X0X1>
-0.246
<Y0Y1>
-0.219
Measurement Counts by Basis
VQE energy: -1.1415 Ha (FCI: -1.1373 Ha). Error: 2.6 kcal/mol. Outside chemical accuracy. Post-selection kept 100% of Z-basis shots.
View cQASM circuit
version 3.0 qubit[2] q bit[2] b // Subspace-preserving ansatz: Ry-CNOT-X (cnot_folds=3) // State = cos(a/2)|10> + sin(a/2)|01>, alpha=-0.223400 Ry(-0.223400) q[0] CNOT q[0], q[1] CNOT q[0], q[1] CNOT q[0], q[1] X q[0] // Z-basis measurement b = measure q
Measured Energy
-0.9147 Ha
FCI Reference
-0.9361 Ha
Energy Level Diagram
Expectation Values
<Z0>
-0.205
<Z1>
0.205
<Z0Z1>
-1.000
<X0X1>
-0.914
<Y0Y1>
-0.892
Measurement Counts by Basis
VQE energy: -0.9147 Ha (FCI: -0.9361 Ha). Error: 13.4 kcal/mol. Outside chemical accuracy. Post-selection kept 96% of Z-basis shots.
View cQASM circuit
version 3.0 qubit[5] q bit[5] b // Subspace-preserving ansatz: Ry-CNOT-X // State = cos(a/2)|10> + sin(a/2)|01>, alpha=-1.380900 Ry(-1.380900) q[2] CNOT q[2], q[4] X q[2] // Z-basis measurement b = measure q
Measured Energy
-1.1325 Ha
FCI Reference
-1.1373 Ha
Energy Level Diagram
Expectation Values
<Z0>
-0.974
<Z1>
0.974
<Z0Z1>
-1.000
<X0X1>
-0.180
<Y0Y1>
-0.219
Measurement Counts by Basis
VQE energy: -1.1325 Ha (FCI: -1.1373 Ha). Error: 3.0 kcal/mol. Outside chemical accuracy. Post-selection kept 96% of Z-basis shots.
View cQASM circuit
version 3.0 qubit[5] q bit[5] b // Subspace-preserving ansatz: Ry-CNOT-X // State = cos(a/2)|10> + sin(a/2)|01>, alpha=-0.223400 Ry(-0.223400) q[2] CNOT q[2], q[4] X q[2] // Z-basis measurement b = measure q
Measured Energy
-0.9362 Ha
FCI Reference
-0.9361 Ha
Energy Level Diagram
Expectation Values
<Z0>
-0.183
<Z1>
0.183
<Z0Z1>
-1.000
<X0X1>
-0.983
<Y0Y1>
-0.984
Measurement Counts by Basis
Measured Energy
-0.9454 Ha
FCI Reference
-0.9486 Ha
Energy Level Diagram
Expectation Values
<Z0>
-0.399
<Z1>
0.399
<Z0Z1>
-1.000
<X0X1>
-0.915
<Y0Y1>
-0.895
Measurement Counts by Basis
Measured Energy
-0.9940 Ha
FCI Reference
-0.9981 Ha
Energy Level Diagram
Expectation Values
<Z0>
-0.742
<Z1>
0.742
<Z0Z1>
-1.000
<X0X1>
-0.646
<Y0Y1>
-0.659
Measurement Counts by Basis
Measured Energy
-1.1034 Ha
FCI Reference
-1.1012 Ha
Energy Level Diagram
Expectation Values
<Z0>
-0.939
<Z1>
0.939
<Z0Z1>
-1.000
<X0X1>
-0.347
<Y0Y1>
-0.361
Measurement Counts by Basis
Measured Energy
-1.1364 Ha
FCI Reference
-1.1373 Ha
Energy Level Diagram
Expectation Values
<Z0>
-0.976
<Z1>
0.976
<Z0Z1>
-1.000
<X0X1>
-0.200
<Y0Y1>
-0.229
Measurement Counts by Basis
Measured Energy
-1.0583 Ha
FCI Reference
-1.0552 Ha
Energy Level Diagram
Expectation Values
<Z0>
-0.992
<Z1>
0.992
<Z0Z1>
-1.000
<X0X1>
-0.149
<Y0Y1>
-0.147
Measurement Counts by Basis
Measured Energy
-0.3105 Ha
FCI Reference
-1.1373 Ha
Energy Level Diagram
Expectation Values
<Z0>
-0.928
<Z1>
-0.939
<Z0Z1>
0.934
<X0X1>
0.174
<Y0Y1>
-0.107
Measurement Counts by Basis
Measured Energy
-0.9210 Ha
FCI Reference
-0.9486 Ha
Energy Level Diagram
Expectation Values
<Z0>
-0.437
<Z1>
0.437
<Z0Z1>
-1.000
<X0X1>
-0.759
<Y0Y1>
-0.827
Measurement Counts by Basis
VQE energy: -0.9210 Ha (FCI: -0.9486 Ha). Error: 17.3 kcal/mol. Outside chemical accuracy. Post-selection kept 96% of Z-basis shots.
View cQASM circuit
version 3.0 qubit[5] q bit[5] b // Subspace-preserving ansatz: Ry-CNOT-X // State = cos(a/2)|10> + sin(a/2)|01>, alpha=-1.133300 Ry(-1.133300) q[2] CNOT q[2], q[4] X q[2] // Z-basis measurement b = measure q
Measured Energy
-0.9779 Ha
FCI Reference
-0.9981 Ha
Energy Level Diagram
Expectation Values
<Z0>
-0.711
<Z1>
0.711
<Z0Z1>
-1.000
<X0X1>
-0.624
<Y0Y1>
-0.611
Measurement Counts by Basis
VQE energy: -0.9779 Ha (FCI: -0.9981 Ha). Error: 12.7 kcal/mol. Outside chemical accuracy. Post-selection kept 96% of Z-basis shots.
View cQASM circuit
version 3.0 qubit[5] q bit[5] b // Subspace-preserving ansatz: Ry-CNOT-X // State = cos(a/2)|10> + sin(a/2)|01>, alpha=-0.726700 Ry(-0.726700) q[2] CNOT q[2], q[4] X q[2] // Z-basis measurement b = measure q
Measured Energy
-1.0946 Ha
FCI Reference
-1.1012 Ha
Energy Level Diagram
Expectation Values
<Z0>
-0.937
<Z1>
0.937
<Z0Z1>
-1.000
<X0X1>
-0.302
<Y0Y1>
-0.332
Measurement Counts by Basis
VQE energy: -1.0946 Ha (FCI: -1.1012 Ha). Error: 4.1 kcal/mol. Outside chemical accuracy. Post-selection kept 97% of Z-basis shots.
View cQASM circuit
version 3.0 qubit[5] q bit[5] b // Subspace-preserving ansatz: Ry-CNOT-X // State = cos(a/2)|10> + sin(a/2)|01>, alpha=-0.352200 Ry(-0.352200) q[2] CNOT q[2], q[4] X q[2] // Z-basis measurement b = measure q
Measured Energy
-1.0393 Ha
FCI Reference
-1.0552 Ha
Energy Level Diagram
Expectation Values
<Z0>
-0.978
<Z1>
0.978
<Z0Z1>
-1.000
<X0X1>
-0.143
<Y0Y1>
-0.122
Measurement Counts by Basis
VQE energy: -1.0393 Ha (FCI: -1.0552 Ha). Error: 10.0 kcal/mol. Outside chemical accuracy. Post-selection kept 96% of Z-basis shots.
View cQASM circuit
version 3.0 qubit[5] q bit[5] b // Subspace-preserving ansatz: Ry-CNOT-X // State = cos(a/2)|10> + sin(a/2)|01>, alpha=-0.143600 Ry(-0.143600) q[2] CNOT q[2], q[4] X q[2] // Z-basis measurement b = measure q
Measured Energy
-1.1222 Ha
FCI Reference
-1.1373 Ha
Energy Level Diagram
Expectation Values
<Z0>
-0.972
<Z1>
0.972
<Z0Z1>
-1.000
<X0X1>
-0.170
<Y0Y1>
-0.132
Measurement Counts by Basis
VQE energy: -1.1222 Ha (FCI: -1.1373 Ha). Error: 9.5 kcal/mol. Outside chemical accuracy. Post-selection kept 83% of Z-basis shots.
View cQASM circuit
version 3.0 qubit[2] q bit[2] b // Subspace-preserving ansatz: Ry-CNOT-X // State = cos(a/2)|10> + sin(a/2)|01>, alpha=-0.223500 Ry(-0.223500) q[0] CNOT q[0], q[1] X q[0] // Z-basis measurement b = measure q
Measured Energy
-1.0976 Ha
FCI Reference
-1.1373 Ha
Energy Level Diagram
Expectation Values
<Z0>
-0.879
<Z1>
0.879
<Z0Z1>
-1.000
<X0X1>
-0.412
<Y0Y1>
-0.439
Measurement Counts by Basis
Measured Energy
-1.0045 Ha
FCI Reference
-1.1373 Ha
Energy Level Diagram
Expectation Values
<Z0>
-0.723
<Z1>
0.970
<Z0Z1>
-0.709
<X0X1>
-0.060
<Y0Y1>
-0.049
Measurement Counts by Basis
VQE on Tuna-9 achieved -1.0045 Hartree for H2 at 0.735 Å. Error of 0.133 Ha (83.3 kcal/mol) reflects significant hardware noise — Tuna-9 Bell fidelity is 87.3% vs IBM's 99%. X/Y basis correlations are heavily degraded by noise.
View raw JSONMeasured Energy
-1.0956 Ha
FCI Reference
-1.1373 Ha
Energy Level Diagram
Expectation Values
<Z0>
-0.875
<Z1>
0.882
<Z0Z1>
-0.966
<X0X1>
-0.402
<Y0Y1>
-0.426
Measurement Counts by Basis
VQE on IBM Marrakesh achieved -1.0956 Hartree for H2 at 0.735 Å. Error of 0.0417 Hartree (26.2 kcal/mol) exceeds chemical accuracy (1.6 mHartree). Circuit limitations from Z-basis only prevented recovery of X/Y correlations.
View raw JSONMeasured Energy
-1.0408 Ha
FCI Reference
-1.1373 Ha
Energy Level Diagram
Expectation Values
<Z0>
-0.854
<Z1>
0.854
<Z0Z1>
-1.000
<X0X1>
-0.155
<Y0Y1>
-0.154
Measurement Counts by Basis
Discussion
Emulator (0.75 kcal/mol): Within chemical accuracy at equilibrium and across all 14 bond distances, correctly capturing the equilibrium well and dissociation limit.
IBM Torino (0.22 kcal/mol with TREX): Chemical accuracy achieved using Qiskit's EstimatorV2 with resilience_level=1 (twirled readout error extinction). This is a 119x improvement over raw measurement (26.2 kcal/mol). Adding dynamical decoupling or gate twirling makes TREX worse -- a counterintuitive finding that "more mitigation" is not always better.
QI Tuna-9 (0.92 kcal/mol with PS+REM): Chemical accuracy achieved via a hybrid approach: post-selection catches gate leakage in Z-basis measurements, while readout error mitigation (confusion matrix inversion) corrects readout bias in X/Y-basis measurements. Validated on two independent qubit pairs: q[2,4] (0.92) and q[6,8] (1.32 kcal/mol).
IQM Garnet (14.26 kcal/mol with PS+REM): Not chemical accuracy despite very low readout error (1.1-1.3%). Gate and decoherence noise dominate on Garnet, making readout-focused mitigation insufficient. An initial attempt gave 60 kcal/mol due to an IQM-specific CNOT decomposition error: the textbook H·CZ·H decomposition fails because IQM's Hadamard differs from standard by a phase that flips X/Y correlator signs. Using Ry(π/2)·CZ·Ry(-π/2) fixed this.
Key insight: The dominant error source determines which mitigation works. IBM: readout error (TREX wins). Tuna-9: readout error (PS+REM wins). IQM Garnet: gate error (readout mitigation alone insufficient). Platform-specific gate decompositions can introduce subtle phase errors invisible to standard tests.
Sources & References
- Peruzzo et al. "A variational eigenvalue solver on a photonic quantum processor" (2014)https://doi.org/10.1038/ncomms5213
- Sagastizabal et al. "Experimental error mitigation via symmetry verification" (2019)https://doi.org/10.1103/PhysRevA.100.010302
- Kandala et al. "Hardware-efficient variational quantum eigensolver" (2017)https://doi.org/10.1038/nature23879
- O'Malley et al. "Scalable quantum simulation of molecular energies" (2016)https://doi.org/10.1103/PhysRevX.6.031007