The spatial photonic Ising machine has achieved remarkable advancements in solving combinatorial optimization problems. However, it still remains a huge challenge to flexibly mapping an arbitrary problem to Ising model. In this paper, we propose a general spatial photonic Ising machine based on interaction matrix eigendecomposition method. Arbitrary interaction matrix can be configured in the two-dimensional Fourier transformation based spatial photonic Ising model by using values generated by matrix eigendecomposition. The error in the structural representation of the Hamiltonian decreases substantially with the growing number of eigenvalues utilized to form the Ising machine. In combination with the optimization algorithm, as low as 65% of the eigenvalues is required by intensity modulation to guarantee the best probability of optimal solution for a 20-vertex graph Max-cut problem, and this probability decreases to below 20% for zero best chance. Our work provides a viable approach for spatial photonic Ising machines to solve arbitrary combinatorial optimization problems with the help of multi-dimensional optical property.The spatial photonic Ising machine has achieved remarkable advancements in solving combinatorial optimization problems. However, it still remains a huge challenge to flexibly mapping an arbitrary problem to Ising model. In this paper, we propose a general spatial photonic Ising machine based on interaction matrix eigendecomposition method. Arbitrary interaction matrix can be configured in the two-dimensional Fourier transformation based spatial photonic Ising model by using values generated by matrix eigendecomposition. The error in the structural representation of the Hamiltonian decreases substantially with the growing number of eigenvalues utilized to form the Ising machine. In combination with the optimization algorithm, as low as 65% of the eigenvalues is required by intensity modulation to guarantee the best probability of optimal solution for a 20-vertex graph Max-cut problem, and this probability decreases to below 20% for zero best chance. Our work provides a viable approach for spatial photonic Ising machines to solve arbitrary combinatorial optimization problems with the help of multi-dimensional optical property.

arXiv:2306.10076Jun 16, 20230 citations

Parallel spatial photonic Ising machine using spatial multiplexing for accelerating combinatorial optimization

Suguru Shimomura, Jun Tanida, Yusuke Ogura

A spatial photonic Ising machine (SPIM) handles large-scale combinatorial optimization problems due to optical processing with spatial parallelism. However, iterative feedback in the search for an optimal solution limits processing speed, even though the Ising Hamiltonian is computed optically. We propose a parallel spatial photonic Ising machine (pSPIM) utilizing spatial multiplexing to search for an optimal solution efficiently. By employing grating patterns and encoding multiple configurations of Ising spins in a phase distribution, several Ising Hamiltonians are computed simultaneously. We demonstrated that the optimal solutions of max-cut problems with 100 Ising spins were obtained more frequently as the number of processing units increased. In addition, combining the multicomponent model with parallel processing provides an efficient search for optimal solutions to problems represented by using interaction matrices with a rank greater than one. The spin-update strategy of pSPIM offers an effective approach for efficiently solving large-scale combinatorial optimization problems.A spatial photonic Ising machine (SPIM) handles large-scale combinatorial optimization problems due to optical processing with spatial parallelism. However, iterative feedback in the search for an optimal solution limits processing speed, even though the Ising Hamiltonian is computed optically. We propose a parallel spatial photonic Ising machine (pSPIM) utilizing spatial multiplexing to search for an optimal solution efficiently. By employing grating patterns and encoding multiple configurations of Ising spins in a phase distribution, several Ising Hamiltonians are computed simultaneously. We demonstrated that the optimal solutions of max-cut problems with 100 Ising spins were obtained more frequently as the number of processing units increased. In addition, combining the multicomponent model with parallel processing provides an efficient search for optimal solutions to problems represented by using interaction matrices with a rank greater than one. The spin-update strategy of pSPIM offers an effective approach for efficiently solving large-scale combinatorial optimization problems.

arXiv:2502.18918Optics LettersOct 7, 20254 citations

Quadrature Photonic Spatial Ising Machine

Wenjia Zhang, Wencheng Sun, Yuanyuan Liu +3 more

Abstract The mining in physics and biology for accelerating the hardcore algorithm to solve non-deterministic polynomial (NP) hard problems has inspired a great amount of special-purpose ma-chine models. Ising machine has become an efficient solver for various combinatorial optimization problems. As a computing accelerator, large-scale photonic spatial Ising machine have great advantages and potentials due to excellent scalability and compact system. However, current fundamental limitation of photonic spatial Ising machine is the configuration flexibility of problem implementation in the accelerator model. Arbitrary spin interaction is highly desired for solving various NP hard problems. Moreover, the absence of external magnetic field in the proposed photonic Ising machine will further narrow the freedom to map the optimization applications. In this paper, we propose a novel quadrature photonic spatial Ising machine to break through the limitation of photonic Ising accelerator by synchronous phase manipulation in two and three sections. Max-cut problem solution with graph order of 100 and density from 0.5 to 1 is experimentally demonstrated after almost 100 iterations. We derive and verify using simulation the solution for Max-cut problem with more than 1600 nodes and the system tolerance for light misalignment. Moreover, vertex cover problem, modeled as an Ising model with external magnetic field, has been successfully implemented to achieve the optimal solution. Our work suggests flexible problem solution by large-scale photonic spatial Ising machine.Abstract The mining in physics and biology for accelerating the hardcore algorithm to solve non-deterministic polynomial (NP) hard problems has inspired a great amount of special-purpose ma-chine models. Ising machine has become an efficient solver for various combinatorial optimization problems. As a computing accelerator, large-scale photonic spatial Ising machine have great advantages and potentials due to excellent scalability and compact system. However, current fundamental limitation of photonic spatial Ising machine is the configuration flexibility of problem implementation in the accelerator model. Arbitrary spin interaction is highly desired for solving various NP hard problems. Moreover, the absence of external magnetic field in the proposed photonic Ising machine will further narrow the freedom to map the optimization applications. In this paper, we propose a novel quadrature photonic spatial Ising machine to break through the limitation of photonic Ising accelerator by synchronous phase manipulation in two and three sections. Max-cut problem solution with graph order of 100 and density from 0.5 to 1 is experimentally demonstrated after almost 100 iterations. We derive and verify using simulation the solution for Max-cut problem with more than 1600 nodes and the system tolerance for light misalignment. Moreover, vertex cover problem, modeled as an Ising model with external magnetic field, has been successfully implemented to achieve the optimal solution. Our work suggests flexible problem solution by large-scale photonic spatial Ising machine.

Jun 16, 20212 citations

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.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.09161Jul 12, 20241 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.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 LettersMay 20, 202512 citations

Quadrature photonic spatial Ising machine

Wenchen Sun, Wenjia Zhang, Yuanyuan Liu +2 more

As a computing accelerator, a large-scale photonic spatial Ising machine has great advantages and potential due to its excellent scalability and compactness. However, the current fundamental limitation of a photonic spatial Ising machine is the configuration flexibility for problem implementation in the accelerator model. Arbitrary spin interactions are highly desired for solving various non-deterministic polynomial (NP)-hard problems. In this paper, we propose a novel quadrature photonic spatial Ising machine to break through the limitation of the photonic Ising accelerator by synchronous phase manipulation in two sections. The max-cut problem solution with a graph order of 100 and density from 0.5 to 1 is experimentally demonstrated after almost 100 iterations. Our work suggests flexible problem solving by the large-scale photonic spatial Ising machine.As a computing accelerator, a large-scale photonic spatial Ising machine has great advantages and potential due to its excellent scalability and compactness. However, the current fundamental limitation of a photonic spatial Ising machine is the configuration flexibility for problem implementation in the accelerator model. Arbitrary spin interactions are highly desired for solving various non-deterministic polynomial (NP)-hard problems. In this paper, we propose a novel quadrature photonic spatial Ising machine to break through the limitation of the photonic Ising accelerator by synchronous phase manipulation in two sections. The max-cut problem solution with a graph order of 100 and density from 0.5 to 1 is experimentally demonstrated after almost 100 iterations. Our work suggests flexible problem solving by the large-scale photonic spatial Ising machine.

arXiv:2105.09307Optics LettersFeb 9, 202235 citations

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.