Case Study: Is This Problem Actually Hard Classically?
Executive Summary
Quantum simulation is the most defensible quantum application, and "defensible" is not "automatic." A great many proposed simulation targets are comfortably solved by existing classical methods, and proposing them wastes the credibility that the genuinely hard cases deserve.
This case study builds the triage: given a candidate system, determine whether DMRG, coupled cluster, or quantum Monte Carlo already handles it, and — if not — identify which property defeats them. Getting this right is what separates a serious quantum simulation proposal from a press release.
Skills applied
- Assessing entanglement scaling and area laws (§17.14).
- Recognizing the fermionic sign problem and where it bites (§17.15).
- Evaluating coupled-cluster applicability from correlation strength.
- Selecting the correct classical baseline for a quantum advantage claim.
Phase 1: The three classical champions
DMRG / tensor networks. Essentially exact for 1D systems with area-law entanglement. Cost scales with bond dimension $\chi$, which must grow as $e^{S}$ where $S$ is the entanglement entropy across a cut.
- Wins: 1D chains, ladders, quasi-1D, gapped systems, many 2D systems of modest width.
- Loses: genuinely 2D/3D systems with volume-law entanglement, where $\chi$ becomes exponential.
Coupled cluster (CCSD(T)). The "gold standard" of quantum chemistry. Polynomial cost, $O(N^7)$ for CCSD(T).
- Wins: weakly correlated molecules — most organic chemistry, most of the periodic table's main group, equilibrium geometries. Accuracy routinely ~1 kcal/mol.
- Loses: strong (static) correlation — stretched bonds, bond breaking, transition metals with near-degenerate d orbitals, open-shell systems, excited states. It can fail catastrophically and without warning, which is the dangerous part.
Quantum Monte Carlo. Statistically exact where applicable, polynomial cost.
- Wins: bosonic systems, spin systems on bipartite lattices, sign-problem-free fermionic models.
- Loses: systems with a fermionic sign problem — frustrated magnets, doped Hubbard models, most realistic fermionic systems away from half filling.
Phase 2: The triage flowchart
For a candidate system, ask in order:
- Is it 1D or quasi-1D with an entanglement area law? → DMRG solves it. No quantum advantage.
- Is it weakly correlated (single-reference, small $T_1$ diagnostic)? → CCSD(T) solves it. No quantum advantage.
- Is it sign-problem-free? → QMC solves it. No quantum advantage.
- Does it fail all three? → Candidate for genuine quantum advantage.
The systems surviving step 4 share a profile: 2D or 3D, strongly correlated, fermionic, and frustrated or doped. That is a narrow class — and it contains several of the most scientifically important open problems in physics and chemistry.
Phase 3: Working examples
| System | DMRG | CCSD(T) | QMC | Verdict |
|---|---|---|---|---|
| H₂O equilibrium geometry | ✓ | ✓ | ✓ | Classically solved. A common but meaningless quantum demo |
| 1D Heisenberg chain | ✓ | — | ✓ | Classically solved |
| Benzene ground state | ~ | ✓ | ✓ | Classically solved |
| N₂ dissociation curve | ✓ (small) | ✗ (breaks down) | ~ | Classical methods exist (CASPT2, DMRG); marginal |
| 2D Hubbard model, doped | ✗ | ✗ | ✗ (sign) | Genuine candidate — and a central open problem |
| FeMoco active space | ✗ | ✗ | ✗ | Genuine candidate |
| Frustrated 2D spin liquid | ✗ | — | ✗ (sign) | Genuine candidate |
| Protein folding energetics | — | ✓ (force fields adequate) | — | Not a quantum problem — classical MD suffices |
The last row is worth dwelling on: protein folding appears on nearly every "quantum applications" slide, and it is dominated by classical force fields and conformational sampling, not by electronic-structure accuracy. It is a sampling problem, not a correlated-electron problem.
Phase 4: The doped 2D Hubbard model
The strongest single case, worth understanding in detail:
$$H = -t\sum_{\langle ij\rangle\sigma}\left(c^\dagger_{i\sigma}c_{j\sigma} + \text{h.c.}\right) + U\sum_i n_{i\uparrow}n_{i\downarrow}$$
Why each classical method fails:
- DMRG: the model is 2D, so entanglement across a cut grows with system width (volume law in the relevant regime); bond dimension becomes exponential.
- CCSD(T): at intermediate $U/t$ the system is strongly correlated and multireference; the single-reference expansion diverges.
- QMC: away from half filling the determinant weights change sign, and the variance grows exponentially with system size and inverse temperature.
And the stakes are real: the doped 2D Hubbard model is widely believed to contain the essential physics of high-temperature superconductivity, unresolved after nearly four decades.
Resource estimate: a $10\times10$ lattice needs ~200 qubits for the fermionic modes plus ancillas, with time evolution costing $10^{6}$–$10^{8}$ gates depending on method and accuracy — error-corrected, but far cheaper than a chemistry-scale QPE calculation. This is among the most plausible early fault-tolerant applications.
Phase 5: How to state a claim honestly
Weak: "We simulated H₂ on a quantum computer." (Solvable by hand; useful only as a hardware demonstration, and should be labelled as such.)
Weak: "Quantum computers will revolutionize drug discovery." (Most of drug discovery is classical MD, docking, and statistics; the electronic-structure fraction is small and often adequately handled.)
Strong: "The doped 2D Hubbard model at $U/t = 8$ resists DMRG (volume-law entanglement in 2D), coupled cluster (multireference character), and QMC (sign problem). A fault-tolerant quantum simulation with ~200 logical qubits could resolve the ground-state phase diagram. We estimate X logical qubits and Y Toffoli gates."
The strong version names the classical methods, states why each fails, and gives a resource estimate. That is the form a defensible claim takes.
Phase 6: The meta-lesson
The best classical algorithms are extremely good, and they improve. DMRG's reach in 2D has expanded considerably; QMC sign-problem mitigations continue to advance; machine-learned wavefunctions are a fast-moving area. A quantum advantage claim is a claim about the classical state of the art, and it therefore requires tracking that state of the art.
This is the same lesson as the Sycamore result from Chapter 1: the classical baseline moved and the claimed margin shrank. Assume yours will too, and build the margin accordingly.
Discussion Questions
- Why does entanglement scaling determine whether DMRG succeeds? Connect area law to bond dimension.
- The sign problem is NP-hard in general. Does that mean it can never be circumvented for specific models?
- Protein folding appears on most quantum applications slides. Construct the strongest argument for including it, then rebut it.
- Classical methods keep improving. How large a margin should a quantum advantage claim require to be credible over a 15-year horizon?
Your Turn: Extensions
- Compute entanglement entropy across a cut for 1D and 2D Heisenberg models and observe the scaling difference.
- Run CCSD(T) on N₂ at equilibrium and at 3× bond length; observe the breakdown.
- Find a published quantum-chemistry quantum-advantage claim and check its classical baseline against the triage.
- Estimate the bond dimension DMRG would need for a $10\times10$ Hubbard lattice.
Key Takeaways
- Most proposed quantum simulation targets are already solved by DMRG, coupled cluster, or QMC; proposing them is a credibility cost.
- Run the triage: area-law/1D → DMRG; weakly correlated → CCSD(T); sign-problem-free → QMC. Only systems failing all three are candidates.
- The surviving class is 2D/3D, strongly correlated, fermionic, frustrated or doped — narrow, but containing genuinely important open problems.
- The doped 2D Hubbard model is the strongest case: it defeats all three classical methods for identifiable reasons and is within reach of early fault-tolerant machines.
- A quantum advantage claim is a claim about the classical state of the art, which keeps improving. Name the methods, say why each fails, and leave margin.