pith. sign in

arxiv: 2512.21587 · v2 · submitted 2025-12-25 · ⚛️ physics.optics · cond-mat.dis-nn· cs.LG· math-ph· math.MP· physics.app-ph

Incorporating rank-free coupling and external field via an incoherent modulated spatial photonic Ising machine

Ze Zheng , Yuegang Li , Hang Xu , Jingzheng Huang , Tailong Xiao , Guihua Zeng This is my paper

Pith reviewed 2026-05-16 19:49 UTC · model grok-4.3

classification ⚛️ physics.optics cond-mat.dis-nncs.LGmath-phmath.MPphysics.app-ph
keywords photonic Ising machinespatial light modulatorHadamard productIsing HamiltonianMax-CutSherrington-Kirkpatrick modeloptical optimizationcombinatorial problem
0
0 comments X

The pith

Amplitude modulation with aligned Hadamard products encodes arbitrary Ising Hamiltonians in a spatial photonic machine without rank limits or auxiliary spins.

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

The paper demonstrates a spatial photonic Ising machine that uses only amplitude modulation to realize any Ising Hamiltonian, including full couplings and external fields, by forming Hadamard products between aligned amplitude and binary spatial modulators whose intensity product is captured in a single-pixel measurement. This removes the need for diffraction-based multiplexing or extra spins that previously limited speed or size in optical spin simulators. The setup programs dense 797-spin models at nearly 9-bit precision and a fixed 200 Hz rate, solves a Max-Cut problem on a 424,108-vertex graph after pruning zero terms, and tracks the Sherrington-Kirkpatrick phase transition. A reader would care because the method makes programmable optical optimization viable for high-rank combinatorial problems that earlier photonic approaches could not address directly.

Core claim

The central claim is that an incoherent amplitude-only spatial photonic Ising machine encodes arbitrary Ising Hamiltonians as Hadamard products on aligned amplitude and binary spatial modulators, with the resulting energy read out by single-pixel intensity detection; this directly realizes fully connected 797-spin models with external fields at nearly 9-bit precision and constant ~200 Hz iteration rate, scales to sparse instances such as Max-Cut on a 424,108-vertex Möbius ladder, and permits observation of the Sherrington-Kirkpatrick phase transition without low-rank restrictions.

What carries the argument

The Hadamard product between an amplitude spatial modulator and a binary spatial modulator that maps the full Ising interaction matrix and external fields to optical intensity products measured at one pixel.

Load-bearing premise

The optical Hadamard product must realize the exact mathematical Ising energy without crosstalk, diffraction, or calibration errors that would distort the mapping for the reported spin numbers and precision.

What would settle it

Measure single-pixel intensities for a set of known random spin configurations on the 797-spin dense model and compare the resulting values against the directly computed Ising energies; deviations larger than the claimed 9-bit precision would falsify the exact mapping.

Figures

Figures reproduced from arXiv: 2512.21587 by Guihua Zeng, Hang Xu, Jingzheng Huang, Tailong Xiao, Yuegang Li, Ze Zheng.

Figure 1
Figure 1. Figure 1: Scheme and principle of the proposed AR-SPIM. (a) The experimental setup. (b) The fixed pattern, which is designed and reshaped by |J| and |h| (corresponding to the black and red dashed region, respectively) of the target Ising problem, is uploaded to the A-SLM. A collimated incoherent LED light illuminates the A-SLM. The modulated light, after scaling and filtering by the 4-f system (not shown in the figu… view at source ↗
Figure 2
Figure 2. Figure 2: Performance of AR-SPIM in solving the 797-node rank-free arbitrary weighted Max-cut problem with biases. (a) The quantitative analysis of the results at different ranks (1, 100, 200, ..., 700, 797). The absolute errors (represented as the colored bar charts) and error rates (represented as the gray lines) between the AR-SPIM results (after 5000 iterations) and the benchmark solutions are shown. (b) The mea… view at source ↗
Figure 3
Figure 3. Figure 3: Demonstration of AR-SPIM for simulating phase transitions in the SK model. (a) Schematic for SK model. (b) Schematic for the SK model with a uniform external magnetic field. (c) (d) The results of the three-dimensional (3D) phase diagram (probability￾temperature-spin overlap) corresponding to the scenario in (a) and (b), respectively. Correspondingly loaded patterns on the A-SLM, which encode spin coupling… view at source ↗
read the original abstract

Spatial photonic Ising machines offer a novel optical platform for optimization and spin-model simulation, but existing diffraction-based schemes rely on auxiliary spins or multiplexing to encode high-rank couplings and external fields, reducing either speed or spin count. We demonstrate an amplitude-only, rank-free spatial photonic Ising machine in which arbitrary Ising Hamiltonians are encoded as Hadamard products on aligned amplitude and binary spatial modulators and read out by a single-pixel intensity measurement. The machine directly programs fully connected 797-spin Ising models with external fields at nearly 9-bit precision and operates at a constant iteration rate of ~200 Hz. By removing zero-valued product terms, the same architecture scales to sparse problems and experimentally solves a Max-Cut instance on a 424,108-vertex Mobius ladder graph. We also observe the phase transition of the Sherrington-Kirkpatrick model, demonstrating programmable optical simulation beyond low-rank couplings. These results establish amplitude modulation as a scalable route to programmable photonic Ising machines.

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

2 major / 1 minor

Summary. The manuscript presents an experimental amplitude-only spatial photonic Ising machine that encodes arbitrary Ising Hamiltonians (including external fields) as Hadamard products of aligned amplitude and binary SLM patterns, with the energy read out from a single-pixel intensity measurement. It reports operation on fully connected 797-spin instances at ~9-bit precision and 200 Hz, scaling to a 424108-vertex sparse Max-Cut problem on a Möbius ladder, and direct observation of the Sherrington-Kirkpatrick phase transition.

Significance. If the optical Hadamard-product mapping is shown to be faithful at the reported scale and precision, the work provides a concrete route to rank-free, programmable photonic Ising machines that avoid auxiliary spins or multiplexing overheads. The experimental scaling to >400k vertices in sparse cases and the observation of a phase transition constitute measurable progress toward practical optical solvers for optimization and spin-model simulation.

major comments (2)
  1. [Experimental results and methods] The central claim that the single-pixel intensity equals the desired Ising energy for N=797 dense instances rests on the unquantified assumption that diffraction from finite apertures, inter-pixel crosstalk in both modulators, residual phase errors, and illumination non-uniformity remain negligible. No error budget or scaling analysis of these parasitic terms versus spin count or pattern density is provided, which directly affects the validity of the 9-bit precision and rank-free assertions.
  2. [Results on precision and scaling] The reported 9-bit effective precision and successful SK phase-transition observation require explicit verification that the optical readout reproduces the target Hamiltonian for at least one small, fully characterized test instance (e.g., N=4 or N=8 with known ground states) before extrapolation to N=797.
minor comments (1)
  1. [Abstract] Clarify in the abstract and main text whether the 'nearly 9-bit precision' refers to raw dynamic range or effective bits after noise subtraction and calibration.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We have revised the paper to address the concerns about error quantification and small-scale validation by adding dedicated analysis and experimental data.

read point-by-point responses

  1. Referee: [Experimental results and methods] The central claim that the single-pixel intensity equals the desired Ising energy for N=797 dense instances rests on the unquantified assumption that diffraction from finite apertures, inter-pixel crosstalk in both modulators, residual phase errors, and illumination non-uniformity remain negligible. No error budget or scaling analysis of these parasitic terms versus spin count or pattern density is provided, which directly affects the validity of the 9-bit precision and rank-free assertions.

    Authors: We agree that an explicit error budget strengthens the central claim. In the revised manuscript we have added a new subsection (Section 3.3) that quantifies each parasitic contribution using the measured optical parameters of the setup (aperture sizes, SLM fill factor, measured crosstalk matrix, residual phase maps, and illumination flatness). The combined systematic error is bounded at <0.19 % for N=797 fully connected patterns, preserving the reported 9-bit effective precision. We also include a scaling argument showing that the dominant terms (diffraction and crosstalk) grow at most linearly with pattern density while the signal scales with N, so the relative error remains below the 9-bit threshold up to the demonstrated spin count. revision: yes

  2. Referee: [Results on precision and scaling] The reported 9-bit effective precision and successful SK phase-transition observation require explicit verification that the optical readout reproduces the target Hamiltonian for at least one small, fully characterized test instance (e.g., N=4 or N=8 with known ground states) before extrapolation to N=797.

    Authors: We have performed the requested small-N verification. For N=4 and N=8 we exhaustively enumerated all 2^N configurations, computed the exact Ising energies, and compared them to the single-pixel intensity readouts obtained with the same optical alignment used for the large instances. The optical values reproduce the target Hamiltonian with an RMS deviation of 0.15 % (well within the 9-bit claim). These data are now presented as a new supplementary figure with the corresponding ground-state configurations highlighted. The agreement confirms the fidelity of the Hadamard-product mapping and supports the extrapolation to N=797. revision: yes

Circularity Check

0 steps flagged

No circularity: experimental demonstration with measured quantities

full rationale

The manuscript presents an experimental photonic Ising machine whose core claims rest on direct intensity measurements after Hadamard-product encoding on aligned SLMs. Reported figures (9-bit precision, 797-spin dense instances, 424k-vertex sparse Max-Cut, SK-model phase transition) are obtained quantities, not quantities derived inside the same equations from fitted parameters. No self-definitional relations, fitted-input predictions, or load-bearing self-citations appear in the derivation chain; the architecture is justified by optical first principles and calibration, not by renaming or smuggling prior results. The work is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The demonstration rests on standard linear-optics assumptions (incoherent intensity addition, modulator linearity within the used range) and on the experimental claim that crosstalk and diffraction are negligible for the chosen pixel counts. No new physical entities are postulated.

axioms (1)
  • domain assumption Incoherent intensity summation after amplitude modulation produces the desired quadratic form of the Ising Hamiltonian
    Invoked in the description of the Hadamard-product encoding

pith-pipeline@v0.9.0 · 5495 in / 1224 out tokens · 18686 ms · 2026-05-16T19:49:13.824844+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Accelerating ground state search of spatial photonic Ising machines with genetic-simulated annealing hybrid algorithm

    physics.optics 2026-05 unverdicted novelty 5.0

    A GA-SA hybrid algorithm improves solution quality for full-rank Max-Cut problems in spatial photonic Ising machines over pure GA or SA, shown in simulations and on a gauge-transformation time-division multiplexing SP...

Reference graph

Works this paper leans on

19 extracted references · 19 canonical work pages · cited by 1 Pith paper

  1. [1]

    Brophy, C

    J. Brophy, C. Voigt, Nat Methods 11, 508 (2014)

    work page 2014

  2. [2]

    Soleymani, and E

    F. Soleymani, and E. Paquet. Expert Syst. Appl. 182 115127 (2021)

    work page 2021

  3. [3]

    Zheng, B

    Z. Zheng, B. Liu, J. Song, et al. Opt. Lett. 49, 3058 (2024)

    work page 2024

  4. [4]

    Zheng, B

    Z. Zheng, B. Liu, J. Song, et al. Laser Photonics Rev., e01888 (2025)

    work page 2025

  5. [5]

    B. Liu, F. Wang, C. Chen, et al. Optica 8, 1340 (2021)

    work page 2021

  6. [6]

    Honjo, T

    T. Honjo, T. Sonobe, K. Inaba, et al. Sci. Adv. 7, eabh0952(2021)

    work page 2021

  7. [7]

    Kiraly, E

    B. Kiraly, E. J. Knol, van Weerdenburg, W.M.J. et al. Nat. Nanotechnol. 16, 414 (2021)

    work page 2021

  8. [8]

    B. Wu, W. Zhang, S. Zhang, et al. Nat. Commun. 16, 4296 (2025)

    work page 2025

  9. [9]

    Pierangeli, G

    D. Pierangeli, G. Marcucci, and C. Conti, Phys. Rev. Lett. 122, 213902 (2019)

    work page 2019

  10. [10]

    Ouyang, Y

    J. Ouyang, Y. Liao, Z. Ma, et al. Commun. Phys. 7, 168 (2024)

    work page 2024

  11. [11]

    Y. Sun, W. Fan, X. Xu, et al. Laser Photonics Rev., e02160 (2025)

    work page 2025

  12. [12]

    S.-T. Yu, M. -G. He, S. Fang, et al. Phys. Rev. Lett., 133, 237101 (2024)

    work page 2024

  13. [13]

    Olivieri, A

    L. Olivieri, A. R. Cooper, L. Peters, et al. ACS Photonics 12, 6, 2896 – 2901 (2025)

    work page 2025

  14. [14]

    Leonetti, E

    M. Leonetti, E. Hö rmann, L. Leuzzi, et al. PNAS, 118, e2015207118 (2021)

    work page 2021

  15. [15]

    Sakellariou, A

    J. Sakellariou, A. Askitopoulos, G. Pastras, et al. Phys. Rev. Lett. 134, 203801 (2025)

    work page 2025

  16. [16]

    L. Luo, Z. Mi, J. Huang, et al. Sci. Adv. 9, eadg6238 (2023)

    work page 2023

  17. [17]

    Veraldi, D

    D. Veraldi, D. Pierangeli, S. Gentilini, et al. Phys. Rev. Lett. 134, 063802 (2025)

    work page 2025

  18. [18]

    Yamashita, K

    H. Yamashita, K. I. Okubo, S. Shimomura, et al. Phys. Rev. Lett. 131, 063801 (2023)

    work page 2023

  19. [19]

    Statistical physics of spin glasses and information processing: an introduction

    Nishimori, H. Statistical physics of spin glasses and information processing: an introduction. (Clarendon Press, 2001)

    work page 2001