Encoding a black-hole mass model in qubits

A black hole as a computational problem

A quantum computer cannot hold an unlimited mathematical description of a black hole. Researchers must first choose a model, reduce it to a manageable size and identify what its calculations actually mean.

In an arXiv preprint, listed as forthcoming in Zeitschrift für Naturforschung A, physicist Babak Vakili develops that translation for the mass observable of a simplified Schwarzschild black-hole interior. Schwarzschild geometry describes an idealized, nonrotating, uncharged black hole.

The work combines mathematical analysis with a small noiseless computer simulation. It does not measure an astronomical black hole, create one in a processor or experimentally test quantum gravity. Its question is narrower: how can this model’s mass operator be encoded and measured using qubits?

Keeping only allowed states

The starting point is a symmetry-reduced description, called minisuperspace, that retains only a few gravitational variables rather than every possible variation of spacetime. Quantization leads to the Wheeler–DeWitt equation, a constraint specifying which quantum states the model permits.

With the paper’s chosen variables, that constraint becomes especially simple. The system resembles two mathematical oscillators, and allowed states have matching oscillator occupation numbers: zero paired with zero, one with one, and so on. These oscillators describe the reduced model; they are not two objects vibrating inside a black hole.

The mass operator connects neighboring matched pairs. Consequently, a state with one definite occupation number is generally not a state with one definite mass. Finding the possible mass values requires solving the operator’s spectral problem, not merely reading its diagonal entries.

Vakili’s analytical solution gives a continuous mass spectrum extending from zero upward within this quantization. That qualification matters: the result concerns this particular model, not every proposed theory of quantum gravity.

A cutoff is not fundamental mass quantization

To obtain a finite computation, the paper retains only the first N allowed basis states. This produces an N-by-N matrix, which has a discrete set of eigenvalues.

Those discrete values are consequences of the computational cutoff. They do not establish that real black holes can possess only particular masses. As the cutoff increases, the lowest matrix eigenvalue approaches zero, consistent with the continuous spectrum derived analytically.

When N is a power of two, the retained physical basis fits into log₂ N qubits. Four basis states therefore need two qubits; 64 need six. This compact encoding stores the shared occupation label rather than maintaining a separate register for each oscillator.

The matrix is then expanded into Pauli strings: products of elementary qubit operators from which its expectation value can be reconstructed.

Compact storage, but several kinds of cost

A general operator on an N-dimensional space can require N² Pauli terms. For this mass matrix, Vakili derives exactly 2N − 1 nonzero terms: seven for four retained states and 127 for 64.

The paper also derives an optimal partition into log₂ N + 1 mutually commuting groups, whose members can be measured together in a suitable common basis. For 64 states, that means seven groups rather than treating all 127 terms separately.

These counts describe different resources. The number of Pauli terms remains exponential in the number of encoding qubits, even though it grows only linearly with N. And a small number of measurement groups is not a complete runtime estimate: repeated measurements, the circuits needed to access their common bases and state preparation also matter. The paper notes that preparing a generic real state can involve N − 1 rotation parameters.

What the benchmark establishes

For the four-state matrix, Vakili uses a variational quantum eigensolver. An adjustable two-qubit circuit prepares trial states, while a classical optimizer minimizes their mass expectation value.

In a noiseless state-vector simulation, the result agrees with exact matrix diagonalization at approximately 0.1612738448 in the model’s units. This is the lowest eigenvalue of that finite matrix—not a measured black-hole mass or a fundamental minimum mass.

The variational benchmark is not repeated at larger cutoffs, and it does not establish performance under realistic gate errors or finite measurement sampling.

Separately, the paper proposes checking the gravitational constraint using two oscillator registers. Measuring both occupation numbers reveals mismatches; averaging their squared difference avoids cancellations that could conceal violations. This requires twice as many qubits as the compact encoding.

The contribution is thus a concrete encoding and resource analysis for one gravitational observable. Black-hole evaporation, thermodynamic entropy and experimentally validated quantum-gravity predictions remain outside its demonstrated scope.