Concept 4 of 7 · How we make it fit

Active spaces

Molecules have thousands of orbitals; almost all of them are boring — full or empty all the time. The active space is the handful of orbitals where the interesting chemistry actually happens, and it’s the only part we put on the quantum computer.

energy Virtual orbitals always empty · frozen out LUMO HOMO Active space → encoded into qubits CAS(4, 4) Core orbitals always full · frozen out CAS(4, 4) → 8 spin-orbitals → 6 circuit qubits
Orbitals split into three energy bands. Core and virtual orbitals are frozen; only the active-space band — where the W–H bond forms — is put on the chip. For WH that is CAS(4, 4): 8 spin-orbitals, tapered to 6 circuit qubits with a Z2 symmetry.

What it actually means

A molecule’s electrons sit in molecular orbitals, conceptually like shells in a building. Real molecules have thousands of orbitals in any realistic basis set. But for almost every chemistry question you might ask:

So you carve out a band in the middle — the active space — containing the highest few occupied orbitals and the lowest few empty orbitals. Those are where electrons can actually be promoted, where bond formation and breaking happens, where electron correlation gets interesting. Only the active space goes onto the quantum computer.

The notation CAS(n, m) means Complete Active Space with n electrons distributed among m orbitals. Each spatial orbital is two qubits (one for spin-up, one for spin-down), so CAS(4, 4) is 8 qubits before optimization. Symmetry tricks (the parity mapper with Z2 tapering) shave that down to 6 qubits — which is what we use for WH.

This isn’t an approximation made up for quantum computers. Classical chemistry has used active spaces for decades — they’re how you make high-accuracy methods like CASSCF and CASPT2 tractable. The quantum adaptation is just: the same active space, encoded into qubits instead of into a classical wavefunction expansion.

Why it matters for our problem

The whole reason WH fits in 6 qubits is the active space.

The caveat

Automatic active-space selection drifts at larger bond distances. Our notebook’s PEC is capped at r = 3.0 Å because qiskit-nature’s ActiveSpaceTransformer starts picking different (non-bonding) orbitals beyond that range. Production-grade work uses natural-orbital or AO-character-based selection, which keeps a consistent active space across geometries.

Scaling up: from 6 qubits to 200+

The fusion-wall problem at scale — tungsten vacancy + embedded hydrogens — needs roughly CAS(20, 20) to capture the multi-centre correlation. In the canonical terminology from the README, that’s ~40 spin-orbital qubits before tapering, ~36 circuit qubits after a Z2 taper. Even bigger problems (a W8 vacancy cluster with multiple H sites, the real production target) push into the ~150+ circuit-qubit regime — which is where the error-correction story starts mattering, because at that circuit depth noise compounds faster than NISQ chips can tolerate. The “200+ logical qubits” figure you’ll see quoted elsewhere is the error-corrected count needed to actually run those circuits to chemical accuracy.

Problem Active space Qubits
WH binding this POC CAS(4, 4) 8 spin-orbitals → 6 circuit (after Z2 taper)
Single W vacancy + embedded H CAS(20, 20) ~40 spin-orbitals → ~36 circuit (after Z2 taper)
W8 vacancy cluster production target multi-centre ~150+ circuit
Fault-tolerant execution error-corrected 200+ logical