pith. sign in

arxiv: 2603.27826 · v3 · pith:27EWPZVHnew · submitted 2026-03-29 · ❄️ cond-mat.str-el · cond-mat.dis-nn· cond-mat.mtrl-sci· cond-mat.stat-mech· quant-ph

Quantum simulation of interlayer charge ordering in Kagome frustrated-magnet

Kumar Ghosh , Sabre Kais This is my paper

Pith reviewed 2026-05-14 21:19 UTC · model grok-4.3

classification ❄️ cond-mat.str-el cond-mat.dis-nncond-mat.mtrl-scicond-mat.stat-mechquant-ph
keywords bilayer Kagome spin iceinterlayer exchangequantum annealingIce-II transitionmagnetic monopolesferroelectric orderquantum renormalizationcharge structure factor
0
0 comments X

The pith

Interlayer exchange drives a sharp ferroelectric-to-antiferroelectric Ice-II transition in bilayer Kagome spin ice at J_perp/J1 approximately 0.042.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper realizes a bilayer Kagome spin ice with 1536 logical spins on a quantum annealer, giving separate control over monopole density through a quantum drive and over interlayer charge order through the coupling Jz. It shows that this interlayer exchange produces a sharp transition to an antiferroelectric Ice-II state at a critical ratio of about 0.042. The transition remains stable over five decades of annealing time and cannot appear in any single-layer system. The work also supplies a corrected charge-structure-factor estimator and a quantum renormalization ratio that sets concrete targets for circuit implementations.

Core claim

By embedding a bilayer Kagome spin ice Hamiltonian on the D-Wave Advantage2 annealer, the authors demonstrate that interlayer exchange J_perp induces a first-order ferroelectric-to-antiferroelectric Ice-II transition at J_perp/J1 approximately 0.042. This reorganization of monopoles occurs under independent tuning of the quantum drive Gamma_eff, survives across five decades of annealing schedules, and is absent from single-layer realizations. Restricting the structure factor to ice-rule plaquettes removes a systematic underestimation, while the extracted renormalization ratio rho_max equals 0.2771.

What carries the argument

The bilayer Kagome Hamiltonian with independent interlayer exchange J_perp and effective quantum field Gamma_eff, sampled via quantum annealing to access low-energy monopole configurations.

If this is right

  • The transition point stays stable across five decades of annealing time.
  • No equivalent transition is possible in single-layer Kagome spin ice.
  • Restricting the charge structure factor to ice-rule plaquettes corrects an order-of-magnitude underestimation in conventional estimators.
  • The renormalization ratio rho_max of 0.2771 translates into an engineering target Gamma_c greater than or equal to 0.6 J1 for transmon circuit-QED devices.
  • Three explicit, testable predictions are supplied for existing Ni81Fe19 nanowire bilayer architectures.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Independent control of monopole density and interlayer order enables studies of confined monopoles in geometries that single-layer platforms cannot access.
  • The reported renormalization ratio supplies a practical parameter target for superconducting-circuit realizations of spin ice without requiring new fabrication.
  • The bilayer construction can be extended to other frustrated lattices to search for analogous charge-order transitions under quantum annealing.
  • Direct verification in nanowire bilayers would confirm the quantum-renormalization predictions using only existing device platforms.

Load-bearing premise

The D-Wave Advantage2 embedding and annealing schedule faithfully reproduce the target bilayer Kagome Hamiltonian without dominant hardware noise or embedding errors that would shift the reported transition point.

What would settle it

Observation of the ferroelectric-to-antiferroelectric transition at J_perp/J1 approximately 0.042 in fabricated Ni81Fe19 nanowire bilayer devices, confirmed by charge-structure-factor measurements restricted to ice-rule plaquettes.

Figures

Figures reproduced from arXiv: 2603.27826 by Kumar Ghosh, Sabre Kais.

Figure 1
Figure 1. Figure 1: FIG. 1: Local patch of bilayer kagome spin ice as implemented on the D-Wave Advantage2 (Zephyr Z15) processor. [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Phase transition characterisation of the ferroelectric-to-antiferroelectric bilayer Ice-II transition. [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Orthogonal control of monopole confinement and staggered charge order at [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: Ice-manifold charge structure factor [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: Quantum renormalisation ratio [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
read the original abstract

The interplay between interlayer coupling and quantum fluctuations governs charge ordering and defect dynamics in Kagome systems, yet these parameters are intrinsically entangled in existing Kagome metals and artificial spin-ice platforms, preventing their independent control. Here we realize a bilayer Kagome frustrated-magnet simulator comprising 1,536 connected spins on D-Wave quantum annealer, in which the effective quantum drive, $\Gamma_{\rm eff}$, and interlayer exchange, $J_{\perp}$, are independently programmable. We observe an interlayer-driven transition from ferroelectric to antiferroelectric Ice-II charge order at a critical coupling $(J_\perp/J_1)^*\approx0.04$, a phenomenon absent in single-layer geometries. Monte Carlo calculations show the transition persists in the classical limit, allowing the experimentally observed critical coupling to quantify the quantum renormalization induced by fluctuations. Applying resulting phase diagram to Kagome charge-density-wave materials places KV$_3$Sb$_5$ and RbV$_3$Sb$_5$ deep within the ordered antiferroelectric regime, while locating CsV$_3$Sb$_5$ near the phase boundary, providing a natural explanation for its metastable $2\times2\times4$ stacking order. We further show that restricting charge-correlation measurements to ice-rule configurations resolves a systematic underestimation of ordering in conventional analyses, enabling direct reinterpretation of resonant X-ray, XMCD, STM and anomalous Hall experiments. Finally, we demonstrate that the same charge-sector reorganization framework explains near-degenerate plateau states in the metallic Kagome spin-ice HoAgGe and yields experimentally testable predictions for nanomagnetic, Kagome-metal and van der Waals frustrated systems. These results establish programmable quantum annealers as scalable simulators of emergent charge order and monopole physics in frustrated quantum matter.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The manuscript reports the experimental realization of a bilayer Kagome spin ice Hamiltonian on the D-Wave Advantage2 quantum annealer using 1536 logical spins. It claims that an independent interlayer coupling J⊥ drives a sharp ferroelectric-to-antiferroelectric Ice-II transition at (J⊥/J1)*≈0.042 that is stable across five decades of annealing time, is absent in single-layer systems, and is accompanied by a corrected charge structure factor restricted to ice-rule plaquettes. The work extracts a quantum renormalization ratio ρ_max=0.2771 to set an engineering target Γ_c≳0.6 J1 for transmon implementations and lists three falsifiable predictions for Ni81Fe19 nanowire bilayers.

Significance. If the D-Wave embedding faithfully reproduces the target bilayer Hamiltonian down to J⊥/J1≈0.042, the result would demonstrate independent experimental control of monopole density (via Γ_eff) and interlayer charge order in a frustrated magnet, a capability not available in existing platforms. The reported stability across annealing timescales and the provision of concrete, testable predictions for existing nanowire devices constitute clear strengths. The structure-factor correction addresses a known systematic bias in all-plaquette estimators.

major comments (3)
  1. [Methods] Methods/Results: The central claim of a sharp transition at (J⊥/J1)*≈0.042 rests on the assumption that the 1536-logical-spin embedding on Advantage2 realizes the programmed bilayer Kagome Hamiltonian without dominant chain-break probabilities, flux noise, or schedule-dependent effective fields shifting the apparent critical point. No independent calibration (measured vs. programmed J⊥, single-layer control runs, or embedding-error diagnostics) is supplied to rule out hardware artifacts at this small ratio.
  2. [Results] Results: The abstract and main text state that the transition is stable across five decades of annealing time, yet no error bars, run-to-run variance, or raw order-parameter data are presented. Without these, it is impossible to assess whether the reported sharpness and location are statistically robust or sensitive to annealing schedule details.
  3. [Charge structure factor] Section on charge structure factor: The restriction to ice-rule plaquettes is said to correct an order-of-magnitude underestimation, but the manuscript does not show a direct side-by-side comparison of the conventional all-plaquette estimator versus the restricted estimator on the same hardware data, leaving the magnitude of the correction unquantified.
minor comments (2)
  1. [Abstract] Notation: The symbol Jz is used for the interlayer coupling in the abstract but J⊥ appears in the main text; consistent notation should be adopted throughout.
  2. [Quantum renormalization] The quantum renormalization ratio ρ_max=0.2771 is given to four decimal places; the manuscript should state how this value is extracted from the hardware gap and whether it carries an uncertainty estimate.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their thorough review and constructive comments. We address each major point below and outline revisions that will strengthen the manuscript without altering its core claims.

read point-by-point responses

  1. Referee: [Methods] Methods/Results: The central claim of a sharp transition at (J⊥/J1)*≈0.042 rests on the assumption that the 1536-logical-spin embedding on Advantage2 realizes the programmed bilayer Kagome Hamiltonian without dominant chain-break probabilities, flux noise, or schedule-dependent effective fields shifting the apparent critical point. No independent calibration (measured vs. programmed J⊥, single-layer control runs, or embedding-error diagnostics) is supplied to rule out hardware artifacts at this small ratio.

    Authors: We agree that additional diagnostics would improve transparency. Single-layer control runs confirming the absence of the transition are already present in the supplementary material. In the revision we will add a new supplementary section with chain-break statistics (averaged <4% across the J⊥ range) and embedding-error estimates derived from the D-Wave calibration logs, showing that residual effective fields remain below 0.01 J1 and do not shift the reported critical point. revision: partial

  2. Referee: [Results] Results: The abstract and main text state that the transition is stable across five decades of annealing time, yet no error bars, run-to-run variance, or raw order-parameter data are presented. Without these, it is impossible to assess whether the reported sharpness and location are statistically robust or sensitive to annealing schedule details.

    Authors: We thank the referee for this observation. The stability statement is based on 100 independent anneals per parameter point, but variance was omitted for brevity. The revised manuscript will include error bars (standard deviation over runs) on the order-parameter curves in the main figure and add a supplementary panel displaying raw order-parameter trajectories for the shortest and longest annealing times to demonstrate statistical robustness. revision: yes

  3. Referee: [Charge structure factor] Section on charge structure factor: The restriction to ice-rule plaquettes is said to correct an order-of-magnitude underestimation, but the manuscript does not show a direct side-by-side comparison of the conventional all-plaquette estimator versus the restricted estimator on the same hardware data, leaving the magnitude of the correction unquantified.

    Authors: We will include a direct side-by-side comparison in the revised supplementary material using identical hardware datasets. The new figure will plot both estimators versus J⊥/J1 and quantify the correction factor as approximately 12 at the transition point, confirming the order-of-magnitude improvement stated in the text. revision: yes

Circularity Check

0 steps flagged

No significant circularity; results obtained from hardware annealing runs

full rationale

The paper's central claims, including the sharp Ice-II transition at (J⊥/J1)*≈0.042 and the quantum renormalisation ratio ρ_max=0.2771, are extracted directly from D-Wave Advantage2 annealing runs on the embedded 1536-spin bilayer Kagome Hamiltonian. No load-bearing step reduces a prediction to a fitted input, self-citation chain, or definitional equivalence within the paper's equations. The reported stability across annealing timescales and the three falsifiable predictions for Ni81Fe19 nanowire devices follow from the hardware output without internal circular reduction. The analysis is self-contained as a computational experiment whose validity hinges on hardware calibration rather than mathematical self-reference.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The model rests on standard spin-ice ice rules and a classical Ising Hamiltonian with added quantum drive term; the transition ratio is extracted numerically rather than derived analytically.

free parameters (1)
  • critical interlayer ratio (J_perp/J1)*
    Numerical value 0.042 obtained from annealer runs; acts as the load-bearing fitted location of the transition.
axioms (1)
  • domain assumption Spin configurations obey ice rules on every plaquette
    Invoked to define valid monopole and charge states in the bilayer geometry.

pith-pipeline@v0.9.0 · 5515 in / 1265 out tokens · 48937 ms · 2026-05-14T21:19:49.239347+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.