Quantum Pebbles: Discrete Dynamics in Physics

• Quantum Pebbles are discrete, localized entities that model complex systems—from Planck-scale lattice sites in black holes to qubit guidance in anonymous graph search.
• They underpin quantum memory management in reversible computation, enabling efficient ancilla allocation with reported improvements up to 52.77% over classical strategies.
• In astrophysics, they metaphorically represent icy solids whose microphysical properties drive pebble drift, pile-up, and planetesimal formation in protoplanetary discs.

In the supplied literature, “Quantum Pebbles” denotes several distinct constructs unified by discreteness, locality, and information-bearing capacity rather than by a single formalism. The term covers Planck-scale lattice sites in a proposed gravitational crystal inside black holes, stationary qubit-emitting markers for treasure hunt in anonymous graphs, pebble-based abstractions for ancilla management in reversible quantum computation, and, more metaphorically, planet-forming pebbles whose macroscopic evolution is controlled by microphysical ice, binding, and surface-reaction physics (Nikolic, 2015, Gaur et al., 3 Sep 2025, Meuli et al., 2019, Topchieva et al., 2024).

1. Terminological scope and principal usages

In the available sources, “pebbles” are discrete localized units inserted into a larger dynamical system. What changes from field to field is the ontology of the pebble: a lattice site, a qubit source, a live intermediate state, or a drifting solid. “Quantum” likewise varies in meaning: literal quantum states in graph search, quantum-gravity microstructure in black-hole interiors, quantum-circuit state management in reversible computation, and microphysics-derived control parameters in planet-forming discs.

Usage	Pebble object	Representative result
Gravitational microstructure	Planck-scale lattice sites in a crystal phase	$M=\gamma M_0^{2/3}$ remnant scaling
Anonymous graph search	Source emitting identical pure-state qubits	Treasure found in exactly $D$ moves using $D$ quantum pebbles
Quantum memory management	Pebbles as live stored values on a DAG	SAT-based clean-up with average improvement of 52.77%
Protoplanetary discs	mm–cm solids with ice mantles	Pebble pile-up, chemical reprocessing, and planetesimal formation

This distribution of meanings is important because the same noun supports very different technical roles. In one case it labels microscopic degrees of freedom of spacetime; in another it is an advice-bearing device in a distributed algorithm; in another it is a combinatorial proxy for ancilla occupation; in astrophysical applications it is a mesoscale solid whose evolution transmits microscopic chemistry into disc-scale and planetary-scale structure (Nikolic, 2015, Gaur et al., 3 Sep 2025, Meuli et al., 2019, Eistrup et al., 2022).

2. Gravitational crystal and black-hole “quantum pebbles”

In Nikolić’s proposal, Einstein gravity is an effective description of a fluid phase of unknown microscopic degrees of freedom, and a phase transition to a crystal phase can occur under extreme conditions inside a black hole (Nikolic, 2015). The proposed crystal has lattice spacing of order the Planck length,

$$\ell_{\text{Pl}}=\frac{1}{m_{\text{Pl}}}\quad (\hbar=c=1),$$

and its entropy density is taken to be of order the Planckian entropy density,

$s\sim \frac{1}{\ell_{\text{Pl}}^3}.$

In the details supplied for the paper, the Planck-scale lattice sites are explicitly interpreted as the “quantum pebbles”: discrete grains of an underlying quantum geometry aggregated into a macroscopic crystalline core.

The phenomenological model is two-phase. Outside the core, gravity remains in the fluid phase and the geometry is Schwarzschild, with Bekenstein–Hawking entropy

$S_{\text{BH}}=\frac{A}{4}=4\pi M^2,\qquad R=2M.$

Inside the core, GR is assumed to be invalid; the effective metric is taken to be flat Minkowski, and the core entropy is modeled by volume scaling,

$$S_{\text{core}}=\alpha V_{\text{core}}=\alpha \frac{4\pi}{3}r_{\text{core}}^3,$$

with $\alpha\sim \mathcal{O}(1)$. No microscopic Lagrangian, quasi-particle content, or explicit degrees of freedom are specified.

The core grows during evaporation because the fluid-phase entropy decreases while ingoing Hawking partners transport entropy inward. Using Page’s estimate

$$S_{\text{radiation}}=\eta\left(4\pi M_0^2-4\pi M^2\right),\qquad \eta\approx 1.5,$$

the model yields

$$r_{\text{core}}^3=(1+\eta)\,\frac{3}{\alpha}\,\left(M_0^2-M^2\right),$$

schematically written in the paper as

$D$0

When the crystal core reaches the horizon, Hawking radiation is expected to shut off. The resulting remnant mass satisfies

$D$1

and for $D$2,

$D$3

Thus the remnant is macroscopic relative to the Planck scale but much lighter than the initial black hole.

The central informational claim is that the crystal phase stores the information that semiclassical evaporation would otherwise lose. Because $D$4, the storage capacity scales with volume rather than area, and the final remnant can in principle encode information from an arbitrarily large initial black hole (Nikolic, 2015). The paper contrasts this with fixed-mass remnant scenarios, fuzzballs, firewalls, energetic curtains, Planck stars, and holographic expectations. Its main limitation is equally explicit: the model is exploratory and phenomenological, with no microscopic derivation of the crystal phase, no detailed transition dynamics, and no explicit matching conditions at the phase boundary.

3. Quantum pebbles in anonymous graph search

In distributed computing, “quantum pebbles” are introduced as a formal guidance mechanism for an oblivious agent searching for a static treasure in an anonymous graph (Gaur et al., 3 Sep 2025). Vertices have no unique identifiers, edges have local port numbers, the agent has no persistent internal memory, and in each synchronous round it can only observe the current node degree, perform local computation, and move. The paper proves that classical pebbles are too weak in this setting: for any classical pebble placement strategy, no deterministic treasure-hunt algorithm for an oblivious agent always succeeds.

The paper’s definition is precise: a quantum pebble is a source that periodically emits qubits, and all emitted qubits are in the same pure quantum state $D$5. The oracle chooses the state per node so as to encode the outgoing port on a shortest path. For $D$6, the encoding uses the $D$7 and $D$8 bases,

$D$9

For general even $D$0, the construction uses $D$1 rotated bases

$D$2

with

$D$3

Two consecutive ports are encoded by the $D$4 and $D$5 states of a single basis.

Decoding exploits basis dependence rather than stored classical entropy. In the correct basis, outcomes are deterministic; in a wrong basis, the overlap is bounded by

$D$6

At each node, the agent measures $D$7 fresh qubits from the local pebble in each candidate basis, identifies the basis yielding a uniform output string, and maps the sign of that uniform string back to the corresponding port. With this procedure, an oblivious agent can locate the treasure in exactly $D$8 moves using $D$9 quantum pebbles, where $$\ell_{\text{Pl}}=\frac{1}{m_{\text{Pl}}}\quad (\hbar=c=1),$$0 is the length of a shortest path from source to treasure. The required number of measurements per node is

$$\ell_{\text{Pl}}=\frac{1}{m_{\text{Pl}}}\quad (\hbar=c=1),$$1

A frequent misconception, addressed directly in the paper, is that the advantage comes from a pure state somehow storing “more bits.” The paper explicitly rejects that reading: the emitted qubits have von Neumann entropy $$\ell_{\text{Pl}}=\frac{1}{m_{\text{Pl}}}\quad (\hbar=c=1),$$2, and the advantage comes from repeated sampling in multiple bases on identical copies, which makes $$\ell_{\text{Pl}}=\frac{1}{m_{\text{Pl}}}\quad (\hbar=c=1),$$3 possible ports statistically distinguishable to a memoryless agent (Gaur et al., 3 Sep 2025). This establishes a sharp separation between classical and quantum advice in anonymous-graph search.

4. Reversible pebbling and quantum memory management

In quantum compilation and reversible computing, pebbles denote currently stored intermediate values on a computation DAG, and the pebbling game becomes an exact model of ancilla allocation and clean-up (Meuli et al., 2019). Let $$\ell_{\text{Pl}}=\frac{1}{m_{\text{Pl}}}\quad (\hbar=c=1),$$4 be a DAG and $$\ell_{\text{Pl}}=\frac{1}{m_{\text{Pl}}}\quad (\hbar=c=1),$$5 its sinks. A reversible pebbling configuration is a set $$\ell_{\text{Pl}}=\frac{1}{m_{\text{Pl}}}\quad (\hbar=c=1),$$6 of pebbled vertices. A reversible pebbling strategy is a sequence $$\ell_{\text{Pl}}=\frac{1}{m_{\text{Pl}}}\quad (\hbar=c=1),$$7 such that $$\ell_{\text{Pl}}=\frac{1}{m_{\text{Pl}}}\quad (\hbar=c=1),$$8, $$\ell_{\text{Pl}}=\frac{1}{m_{\text{Pl}}}\quad (\hbar=c=1),$$9, each move changes exactly one vertex status, and pebbling or unpebbling a vertex requires all of its children to be pebbled.

The mapping to quantum memory management is direct. Placing a pebble corresponds to computing an intermediate value and storing it in ancilla qubits; removing a pebble corresponds to uncomputing that value and returning its ancillae to $s\sim \frac{1}{\ell_{\text{Pl}}^3}.$0. The requirement that all children be present when pebbling or unpebbling a node is the reversible-computation constraint that all needed inputs remain available for forward computation and inverse clean-up. Under the abstraction adopted in the paper, “the problem of finding a strategy to compute and uncompute intermediate states for a given fixed number of qubits corresponds to solving the reversible pebbling game” (Meuli et al., 2019).

The optimization problem is then encoded as SAT. For a step bound $s\sim \frac{1}{\ell_{\text{Pl}}^3}.$1 and pebble bound $s\sim \frac{1}{\ell_{\text{Pl}}^3}.$2, Boolean variables

$s\sim \frac{1}{\ell_{\text{Pl}}^3}.$3

indicate whether vertex $s\sim \frac{1}{\ell_{\text{Pl}}^3}.$4 is pebbled at time $s\sim \frac{1}{\ell_{\text{Pl}}^3}.$5. The encoding enforces initial and final conditions, move clauses,

$s\sim \frac{1}{\ell_{\text{Pl}}^3}.$6

and cardinality constraints

$s\sim \frac{1}{\ell_{\text{Pl}}^3}.$7

The outer optimization loop increments $s\sim \frac{1}{\ell_{\text{Pl}}^3}.$8 until satisfiable, giving a minimum-step strategy for the chosen qubit bound.

The reported empirical result is an average improvement of 52.77% in pebble count relative to the Bennett-style baseline, at an average multiplicative increase of 2.68× in steps (Meuli et al., 2019). The paper also presents a constrained-hardware example: a 9-input AND on a 16-qubit device. The Bennett strategy requires 17 qubits and thus does not fit; a Barenco decomposition uses 11 qubits but 48 gates; the SAT-derived pebbling strategy fits exactly into 16 qubits and uses 23 gates. The practical message is that pebbling is not merely a metaphor: it is a hardware-aware scheduling formalism for navigating the gate–qubit trade-off.

The main caveats are also explicit. Optimal reversible pebbling is PSPACE-complete; the SAT encoding scales with $s\sim \frac{1}{\ell_{\text{Pl}}^3}.$9; the tool assumes straight-line reversible computations; and connectivity, noise, and coherence are not encoded in the SAT instance (Meuli et al., 2019).

5. Pebbles in protoplanetary discs: drift, pile-up, and planetesimal belts

In planetary-disc research, pebbles are mm–cm solids with Stokes numbers typically in the range $S_{\text{BH}}=\frac{A}{4}=4\pi M^2,\qquad R=2M.$0, partially decoupled from the gas and therefore subject to substantial radial drift. The associated models are classical fluid and kinetic models, but several papers emphasize that the macroscopic behavior of these pebbles is controlled by microscopic material parameters, especially ice binding and sticking properties (Hyodo et al., 2020, Hyodo et al., 2020).

Hyodo, Ida, and Guillot describe a “no-drift” runaway pile-up in discs whose midplane turbulence increases with radius, as in the outer region of a dead zone (Hyodo et al., 2020). The basic feedback is that inward-drifting pebbles entering a low-$S_{\text{BH}}=\frac{A}{4}=4\pi M^2,\qquad R=2M.$1 region settle vertically, their scale height decreases, the midplane density ratio $S_{\text{BH}}=\frac{A}{4}=4\pi M^2,\qquad R=2M.$2 rises, and drag back-reaction on the gas suppresses the pebble drift velocity. At fixed pebble-to-gas mass flux $S_{\text{BH}}=\frac{A}{4}=4\pi M^2,\qquad R=2M.$3, slower drift forces $S_{\text{BH}}=\frac{A}{4}=4\pi M^2,\qquad R=2M.$4 upward, which raises $S_{\text{BH}}=\frac{A}{4}=4\pi M^2,\qquad R=2M.$5 further. This can drive a runaway toward a near “no-drift” state without invoking a pressure bump. The critical midplane turbulence in the turbulence-dominated regime is

$S_{\text{BH}}=\frac{A}{4}=4\pi M^2,\qquad R=2M.$6

and the paper identifies conditions under which this feedback dominates over Kelvin–Helmholtz regulation (Hyodo et al., 2020).

The follow-up study adds planetesimal formation through the streaming instability and shows that the no-drift state produces a finite planetesimal belt rather than an indefinitely growing pile-up (Hyodo et al., 2022). Planetesimals initially form in a narrow ring whose width expands through radial diffusion of accumulating pebbles and then saturates. With nominal parameters, more than one Earth mass of planetesimals forms for a disk having $S_{\text{BH}}=\frac{A}{4}=4\pi M^2,\qquad R=2M.$7 for more than $S_{\text{BH}}=\frac{A}{4}=4\pi M^2,\qquad R=2M.$8 kyr, with $S_{\text{BH}}=\frac{A}{4}=4\pi M^2,\qquad R=2M.$9, $$S_{\text{core}}=\alpha V_{\text{core}}=\alpha \frac{4\pi}{3}r_{\text{core}}^3,$$0, and $$S_{\text{core}}=\alpha V_{\text{core}}=\alpha \frac{4\pi}{3}r_{\text{core}}^3,$$1 at $$S_{\text{core}}=\alpha V_{\text{core}}=\alpha \frac{4\pi}{3}r_{\text{core}}^3,$$2 au (Hyodo et al., 2022). The model therefore converts a drag-feedback instability into a self-regulated belt of solids.

Around the water snow line, Hyodo et al. identify a second family of pile-up phenomena driven by sublimation, recondensation, and recycling between icy pebbles and silicate dust (Hyodo et al., 2020). The sublimation width $$S_{\text{core}}=\alpha V_{\text{core}}=\alpha \frac{4\pi}{3}r_{\text{core}}^3,$$3 is broad in the advection-dominated regime and narrow in the diffusion-dominated regime, which strongly affects the vertical thickness of the released dust layer. The paper distinguishes “runaway pile-up of silicate dust inside the snow line” for $$S_{\text{core}}=\alpha V_{\text{core}}=\alpha \frac{4\pi}{3}r_{\text{core}}^3,$$4 from “runaway pile-up of pebbles outside the snow line” for $$S_{\text{core}}=\alpha V_{\text{core}}=\alpha \frac{4\pi}{3}r_{\text{core}}^3,$$5. The former favors rocky planetesimals inside the snow line; the latter favors icy planetesimals outside it (Hyodo et al., 2020).

This astrophysical body of work uses “pebbles” in the ordinary planet-formation sense, but it also motivates the broader phrase “Quantum Pebbles” by stressing that microscopic parameters—binding energies, fragmentation velocities, adsorption and desorption kinetics—control disc-scale transport, pile-up, and planetesimal formation. This suggests a microphysics-to-macrophysics analogy rather than a formal quantum-information or quantum-gravity construction.

6. Chemical evolution of icy pebbles and planetary composition

The chemical composition carried by pebbles is itself time-dependent. In the FEOSAD-based study of a self-gravitating, viscous protoplanetary disc, pebbles form as early as 50 kyr after disc formation and persist to 500 kyr, with all pebbles covered by icy mantles (Topchieva et al., 2024). The total pebble mass reaches

$$S_{\text{core}}=\alpha V_{\text{core}}=\alpha \frac{4\pi}{3}r_{\text{core}}^3,$$6

and the ice mantles consist mainly of H$$S_{\text{core}}=\alpha V_{\text{core}}=\alpha \frac{4\pi}{3}r_{\text{core}}^3,$$7O and CO$$S_{\text{core}}=\alpha V_{\text{core}}=\alpha \frac{4\pi}{3}r_{\text{core}}^3,$$8, being carbon-depleted compared to gas and ices on small and grown dust. The paper therefore argues that planets formed from these pebbles are plausibly oxygen-rich under the modeled conditions.

A complementary local-chemistry calculation follows 0.6 mm icy pebbles drifting inward from 128 AU and 200 AU on timescales of 10 kyr, 100 kyr, and 1 Myr (Eistrup et al., 2022). For ice species with initial abundances relative to hydrogen $$S_{\text{core}}=\alpha V_{\text{core}}=\alpha \frac{4\pi}{3}r_{\text{core}}^3,$$9, including H$\alpha\sim \mathcal{O}(1)$0O, CO$\alpha\sim \mathcal{O}(1)$1, CH$\alpha\sim \mathcal{O}(1)$2OH, and NH$\alpha\sim \mathcal{O}(1)$3, the abundances change by less than 20% for both radii of origin and for the two smaller drift timescales. For less abundant species, and especially for the 1 Myr drift timescale, the changes are larger. The net trend is that pebble drift chemistry generally increases the ice abundances of CO$\alpha\sim \mathcal{O}(1)$4, HCN, and SO while decreasing other volatile molecules (Eistrup et al., 2022). Thus, for rapid drift the initial outer-disc composition is largely preserved, whereas long residence times permit substantial in-transit reprocessing.

These compositional effects propagate into planet formation calculations. In semi-analytical giant-planet models using the chemcomp code, atmospheric enrichment is dominated by gas that has been enriched by inward-drifting and evaporating pebbles (Danti et al., 2023). When planetesimal formation is included, pebbles locked into planetesimals no longer evaporate into the gas, and the accreted heavy-element content drops sharply. The paper therefore concludes that planetesimal formation needs to be inefficient in order to explain planets with high heavy-element content. Planetesimal accretion, when added, enhances the refractory component of the atmosphere and lowers volatile-to-refractory ratios, whereas pure pebble accretion tends to produce higher atmospheric C/H and O/H (Danti et al., 2023).

A similar compositional bifurcation appears in compact M-dwarf systems. In simulations tailored to TRAPPIST-1 analogues, pebble accretion and planetesimal accretion both reproduce observed masses, periods, and resonant structure, but their water outcomes differ strongly (Coleman et al., 2019). With ablation of icy pebbles in planetary envelopes and full recycling of that envelope with the disc, planets formed from pebbles are extremely dry; if water is not fully recycled, or if ablation is neglected, pebble-formed planets become extremely wet, similar to planets formed from planetesimals. The paper therefore identifies water content, rather than architecture, as the main discriminant between the two growth channels (Coleman et al., 2019).

7. Comparative themes, misconceptions, and open problems

Across these literatures, pebbles act as discrete carriers of state, entropy, memory, or composition. This suggests a shared structural motif: a continuum-level description fails or becomes incomplete, and the relevant correction enters through localized units whose collective behavior is decisive. In black-hole physics the correction is a phase change from Einstein-fluid to gravitational crystal; in anonymous-graph search it is repeated basis-dependent interrogation of identical pure states; in quantum compilation it is the explicit accounting of live intermediate states; in protoplanetary discs it is drag-mediated transport and chemistry of partially decoupled solids (Nikolic, 2015, Gaur et al., 3 Sep 2025, Meuli et al., 2019, Hyodo et al., 2020).

Several misconceptions are explicitly addressed in the source material. In the graph-search setting, the power of quantum pebbles does not come from larger static information content: each emitted qubit is in a pure state, and the advantage comes from measurement flexibility across multiple bases and many identical copies (Gaur et al., 3 Sep 2025). In the gravitational-crystal setting, “quantum pebbles” do not denote a derived microscopic theory; the lattice spacing and entropy density are postulated phenomenologically, with no explicit microdynamics (Nikolic, 2015). In the quantum-memory setting, pebbling is exact only at the chosen level of abstraction: it models straight-line reversible computation with ancilla constraints, not full architecture-specific compilation (Meuli et al., 2019). In the astrophysical literature, the phrase is best read as metaphorical shorthand for the fact that quantum-level or microphysical parameters such as binding energies and fragmentation thresholds regulate disc-scale solid transport and planetary composition (Topchieva et al., 2024, Eistrup et al., 2022).

The open problems are equally diverse. Nikolić’s model lacks a microscopic crystal theory, transition dynamics, and a derivation of remnant stability (Nikolic, 2015). Quantum-pebble search in anonymous graphs invites tighter distinguishability bounds, directed-graph generalizations, and adaptive stopping rules (Gaur et al., 3 Sep 2025). SAT-based reversible pebbling remains limited by PSPACE-hardness and solver scalability (Meuli et al., 2019). Planet-formation models still rely on reduced-dimensional or parcel approximations, prescribed pebble fluxes, simplified turbulence, and limited chemistry; more self-consistent coupling of drift, growth, reactions, and planet formation remains an open agenda (Hyodo et al., 2022, Eistrup et al., 2022, Danti et al., 2023).

Taken together, the literature does not define a single doctrine of “Quantum Pebbles.” It defines a family of high-resolution descriptions in which localized discrete entities—Planck-scale lattice grains, qubit sources, DAG pebbles, or drifting icy solids—supply the missing degrees of freedom required to explain memory, guidance, entropy storage, or compositional transport beyond a coarse-grained continuum picture.