© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

arXiv:2410.10689 · Physical Review Letters

Fully Programmable Spatial Photonic Ising Machine by Focal Plane Division

Daniele Veraldi, Davide Pierangeli, Silvia Gentilini, Marcello Calvanese Strinati, Jason Sakellariou, James Cummins, Airat Kamaletdinov, Marvin Syed, Richard Zhipeng Wang, Natalia G. Berloff, Dimitris Karanikolopoulos, Pavlos G. Savvidis, Claudio Conti

Published Feb 14, 2025• 24 citations

Abstract

Ising machines are an emerging class of hardware that promises ultrafast and energy-efficient solutions to NP-hard combinatorial optimization problems. Spatial photonic Ising machines (SPIMs) exploit optical computing in free space to accelerate the computation, showcasing parallelism, scalability, and low power consumption. However, current SPIMs can implement only a restricted class of problems. This partial programmability is a critical limitation that hampers their benchmark. Achieving full programmability of the device while preserving its scalability is an open challenge. Here, we report a fully programmable SPIM achieved through a novel operation method based on the division of the focal plane. In our scheme, a general Ising problem is decomposed into a set of Mattis Hamiltonians, whose energies are simultaneously computed optically by measuring the intensity on different regions of the camera sensor. Exploiting this concept, we experimentally demonstrate the computation with high success probability of ground-state solutions of up to 32-spin Ising models on unweighted maximum cut graphs with and without ferromagnetic bias. Simulations of the hardware prove a favorable scaling of the accuracy with the number of spin. Our fully programmable SPIM enables the implementation of many quadratic unconstrained binary optimization problems, further establishing SPIMs as a leading paradigm in non-von Neumann hardware.

DOI: 10.1103/physrevlett.134.063802

Related papers

Encoding arbitrary Ising Hamiltonians on Spatial Photonic Ising Machines

Jason Sakellariou, Alexis Askitopoulos, Georgios Pastras +1 more

Photonic Ising Machines constitute an emergent new paradigm of computation, geared towards tackling combinatorial optimization problems that can be reduced to the problem of finding the ground state of an Ising model. Spatial Photonic Ising Machines have proven to be advantageous for simulating fully connected large-scale spin systems. However, fine control of a general interaction matrix $J$ has so far only been accomplished through eigenvalue decomposition methods that either limit the scalability or increase the execution time of the optimization process. We introduce and experimentally validate a SPIM instance that enables direct control over the full interaction matrix, enabling the encoding of Ising Hamiltonians with arbitrary couplings and connectivity. We demonstrate the conformity of the experimentally measured Ising energy with the theoretically expected values and then proceed to solve both the unweighted and weighted graph partitioning problems, showcasing a systematic convergence to an optimal solution via simulated annealing. Our approach greatly expands the applicability of SPIMs for real-world applications without sacrificing any of the inherent advantages of the system, and paves the way to encoding the full range of NP problems that are known to be equivalent to Ising models, on SPIM devices.

arXiv:2407.09161•Jul 12, 2024•1 citation

Encoding Arbitrary Ising Hamiltonians on Spatial Photonic Ising Machines

Jason Sakellariou, Alexis Askitopoulos, Georgios Pastras +1 more

Photonic Ising machines constitute an emergent new paradigm of computation geared toward tackling combinatorial optimization problems that can be reduced to the problem of finding the ground state of an Ising model. Spatial photonic Ising machines (SPIMs) have proven advantageous for simulating fully connected large-scale spin systems. Fine control of a general interaction matrix J has been accomplished so far only through matrix decomposition methods. We introduce and experimentally validate a SPIM instance that enables direct control over the full interaction matrix, allowing the encoding of Ising Hamiltonians with arbitrary couplings and connectivity. We demonstrate the conformity of the experimentally measured Ising energy with the theoretically expected value and then proceed to solve both the unweighted and weighted graph partitioning problems, showcasing a systematic convergence to an optimal solution via simulated annealing. Our approach significantly expands the applicability of SPIMs for real-world applications, as it is more efficient than matrix decomposition methods in the case of sparse problems. It paves the way to encoding the full range of NP problems that are known to be equivalent to Ising models on SPIM devices.

Physical Review Letters•May 20, 2025•12 citations

Multi-channel spatial photonic Ising machine for 102400-spin 0-1 knapsack problem solving

Junze Yao, Junjie Yu

Combinatorial optimization problems are challenging for classical computing architectures, motivating the development of specialized hardware like Ising machines. The Spatial Photonic Ising Machine (SPIM) and its enhanced versions have shown promise in solving such problems by mapping them to the Ising model and leveraging optical parallelism to compute Hamiltonian energies efficiently. However, the 0-1 knapsack problem, as a representative combinatorial optimization problem, has only been reported to be solved effectively in small-scale SPIM. In this work, we analyze the factors that limit the spin scalability and propose a novel hardware-aware Hamiltonian mapping scheme to overcome this limitation. Based on this enhanced SPIM, the scheme is validated through simulated annealing experiments on a 102400-spin 0-1 knapsack problem, achieving high solution accuracy. Our method not only enables efficient solutions to large-scale 0-1 knapsack problems but also offers generalizability to other combinatorial optimization problems with similar coupling configurations. This work demonstrates the potential of SPIMs as scalable and adaptable solvers for real-world optimization challenges.

Optics Letters•Jun 1, 2026•0 citations

Wavelength-division multiplexing optical Ising simulator enabling fully programmable spin couplings and external magnetic fields

Li Luo, Zhiyi Mi, Junyi Huang +1 more

Recently various physical systems have been proposed for modeling Ising spin Hamiltonians appealing to solve combinatorial optimization problems with remarkable performance. However, how to implement arbitrary spin-spin interactions is a critical and challenging problem in unconventional Ising machines. Here, we propose a general gauge transformation scheme to enable arbitrary spin-spin interactions and external magnetic fields as well, by decomposing an Ising Hamiltonian into multiple Mattis-type interactions. With this scheme, a wavelength-division multiplexing spatial photonic Ising machine (SPIM) is developed to show the programmable capability of general spin coupling interactions. We exploit the wavelength-division multiplexing SPIM to simulate three spin systems: ± J models, Sherrington-Kirkpatrick models, and only locally connected J 1 - J 2 models and observe the phase transitions. We also demonstrate the ground-state search for solving Max-Cut problem with the wavelength-division multiplexing SPIM. These results promise the realization of ultrafast-speed and high–power efficiency Boltzmann sampling to a generalized large-scale Ising model.

arXiv:2303.11565•Science Advances•Dec 1, 2023•38 citations

Ising accelerator with a reconfigurable interferometric photonic processor

Rausell-Campo, José Roberto, Kayed, Nayem Al, Pérez-Lppez, Daniel +3 more

The general-purpose programmable photonic processors offer a scalable and reconfigurable solution for a wide range of RF and optical applications. Therefore, implementing photonic Ising machines using programmable processors leverages the advantages of high speed and parallelism, enabling efficient hardware acceleration for finding ground-state solutions to combinatorial optimization problems. In this work, we demonstrate a novel programmable photonic Ising solver based on a hexagonal mesh general-purpose programmable photonic platform. The integrated system allows reconfigurable matrix multiplication and computes the Hamiltonian iteratively using an annealing algorithm that facilitates spin updates and effectively searches for the ground state. As a proof of concept, we experimentally solve two benchmark optimization problems, a fundamental three-node ferromagnetic coupling problem with external bias that demonstrates nontrivial spin interactions, and a four-node Max-Cut problem with arbitrary coupling matrices. Furthermore, to establish a large-scale capability, we emulated Ising problems with sizes up to N = 50, achieving success probabilities exceeding 80\%. Additionally, we examined the impact of errors, such as phase and coupling, on the performance of the programmable photonic Ising machine. Our general-purpose photonic Ising machine paves the way for implementing large-scale, programmable architectures for solving optimization problems.

arXiv:2511.13284•Nov 17, 2025•0 citations

Enhanced, fully connected 360 000-spin spatial photonic Ising machine

Junze Yao, Rongwei Zhu, Junjie Yu

The spatial photonic Ising machine (SPIM) is a promising architecture that leverages the scalability and parallelism of photons to solve Ising problems in combinatorial optimization. SPIMs have been held back fundamentally by limited computational accuracy. This study overcomes that limitation through a hybrid approach, integrating aberration correction with a dynamic algorithm to update simultaneously a number of flipped spins, rather than just one. The method demonstrably improves computational accuracy and presents a way for SPIMs to scale to significantly larger systems (order 10${}^{5}$ spins), a milestone on the path to practical application.

Physical Review Applied•Sep 25, 2025•5 citations

Ising machines are an emerging class of hardware that promises ultrafast and energy-efficient solutions to NP-hard combinatorial optimization problems. Spatial photonic Ising machines (SPIMs) exploit optical computing in free space to accelerate the computation, showcasing parallelism, scalability, and low power consumption. However, current SPIMs can implement only a restricted class of problems. This partial programmability is a critical limitation that hampers their benchmark. Achieving full programmability of the device while preserving its scalability is an open challenge. Here, we report a fully programmable SPIM achieved through a novel operation method based on the division of the focal plane. In our scheme, a general Ising problem is decomposed into a set of Mattis Hamiltonians, whose energies are simultaneously computed optically by measuring the intensity on different regions of the camera sensor. Exploiting this concept, we experimentally demonstrate the computation with high success probability of ground-state solutions of up to 32-spin Ising models on unweighted maximum cut graphs with and without ferromagnetic bias. Simulations of the hardware prove a favorable scaling of the accuracy with the number of spin. Our fully programmable SPIM enables the implementation of many quadratic unconstrained binary optimization problems, further establishing SPIMs as a leading paradigm in non-von Neumann hardware.

Photonic Ising Machines constitute an emergent new paradigm of computation, geared towards tackling combinatorial optimization problems that can be reduced to the problem of finding the ground state of an Ising model. Spatial Photonic Ising Machines have proven to be advantageous for simulating fully connected large-scale spin systems. However, fine control of a general interaction matrix

J

has so far only been accomplished through eigenvalue decomposition methods that either limit the scalability or increase the execution time of the optimization process. We introduce and experimentally validate a SPIM instance that enables direct control over the full interaction matrix, enabling the encoding of Ising Hamiltonians with arbitrary couplings and connectivity. We demonstrate the conformity of the experimentally measured Ising energy with the theoretically expected values and then proceed to solve both the unweighted and weighted graph partitioning problems, showcasing a systematic convergence to an optimal solution via simulated annealing. Our approach greatly expands the applicability of SPIMs for real-world applications without sacrificing any of the inherent advantages of the system, and paves the way to encoding the full range of NP problems that are known to be equivalent to Ising models, on SPIM devices.

Photonic Ising machines constitute an emergent new paradigm of computation geared toward tackling combinatorial optimization problems that can be reduced to the problem of finding the ground state of an Ising model. Spatial photonic Ising machines (SPIMs) have proven advantageous for simulating fully connected large-scale spin systems. Fine control of a general interaction matrix J has been accomplished so far only through matrix decomposition methods. We introduce and experimentally validate a SPIM instance that enables direct control over the full interaction matrix, allowing the encoding of Ising Hamiltonians with arbitrary couplings and connectivity. We demonstrate the conformity of the experimentally measured Ising energy with the theoretically expected value and then proceed to solve both the unweighted and weighted graph partitioning problems, showcasing a systematic convergence to an optimal solution via simulated annealing. Our approach significantly expands the applicability of SPIMs for real-world applications, as it is more efficient than matrix decomposition methods in the case of sparse problems. It paves the way to encoding the full range of NP problems that are known to be equivalent to Ising models on SPIM devices.

Combinatorial optimization problems are challenging for classical computing architectures, motivating the development of specialized hardware like Ising machines. The Spatial Photonic Ising Machine (SPIM) and its enhanced versions have shown promise in solving such problems by mapping them to the Ising model and leveraging optical parallelism to compute Hamiltonian energies efficiently. However, the 0-1 knapsack problem, as a representative combinatorial optimization problem, has only been reported to be solved effectively in small-scale SPIM. In this work, we analyze the factors that limit the spin scalability and propose a novel hardware-aware Hamiltonian mapping scheme to overcome this limitation. Based on this enhanced SPIM, the scheme is validated through simulated annealing experiments on a 102400-spin 0-1 knapsack problem, achieving high solution accuracy. Our method not only enables efficient solutions to large-scale 0-1 knapsack problems but also offers generalizability to other combinatorial optimization problems with similar coupling configurations. This work demonstrates the potential of SPIMs as scalable and adaptable solvers for real-world optimization challenges.

Recently various physical systems have been proposed for modeling Ising spin Hamiltonians appealing to solve combinatorial optimization problems with remarkable performance. However, how to implement arbitrary spin-spin interactions is a critical and challenging problem in unconventional Ising machines. Here, we propose a general gauge transformation scheme to enable arbitrary spin-spin interactions and external magnetic fields as well, by decomposing an Ising Hamiltonian into multiple Mattis-type interactions. With this scheme, a wavelength-division multiplexing spatial photonic Ising machine (SPIM) is developed to show the programmable capability of general spin coupling interactions. We exploit the wavelength-division multiplexing SPIM to simulate three spin systems: ± J models, Sherrington-Kirkpatrick models, and only locally connected J 1 - J 2 models and observe the phase transitions. We also demonstrate the ground-state search for solving Max-Cut problem with the wavelength-division multiplexing SPIM. These results promise the realization of ultrafast-speed and high–power efficiency Boltzmann sampling to a generalized large-scale Ising model.

The general-purpose programmable photonic processors offer a scalable and reconfigurable solution for a wide range of RF and optical applications. Therefore, implementing photonic Ising machines using programmable processors leverages the advantages of high speed and parallelism, enabling efficient hardware acceleration for finding ground-state solutions to combinatorial optimization problems. In this work, we demonstrate a novel programmable photonic Ising solver based on a hexagonal mesh general-purpose programmable photonic platform. The integrated system allows reconfigurable matrix multiplication and computes the Hamiltonian iteratively using an annealing algorithm that facilitates spin updates and effectively searches for the ground state. As a proof of concept, we experimentally solve two benchmark optimization problems, a fundamental three-node ferromagnetic coupling problem with external bias that demonstrates nontrivial spin interactions, and a four-node Max-Cut problem with arbitrary coupling matrices. Furthermore, to establish a large-scale capability, we emulated Ising problems with sizes up to N = 50, achieving success probabilities exceeding 80\%. Additionally, we examined the impact of errors, such as phase and coupling, on the performance of the programmable photonic Ising machine. Our general-purpose photonic Ising machine paves the way for implementing large-scale, programmable architectures for solving optimization problems.

The spatial photonic Ising machine (SPIM) is a promising architecture that leverages the scalability and parallelism of photons to solve Ising problems in combinatorial optimization. SPIMs have been held back fundamentally by limited computational accuracy. This study overcomes that limitation through a hybrid approach, integrating aberration correction with a dynamic algorithm to update simultaneously a number of flipped spins, rather than just one. The method demonstrably improves computational accuracy and presents a way for SPIMs to scale to significantly larger systems (order 10

{}^{5}

spins), a milestone on the path to practical application.