Solving large-scale quadratic unconstrained binary optimization (QUBO) problems is critical in various fields, including physics, finance, and engineering. However, these problems remain intractable on conventional computing architectures. Alternative solvers, such as Ising machines (IMs) based on networks of coupled electronic, mechanical, or photonic parametric oscillators (POs), have recently been developed. PO-based IMs aim to find the ground state of an Ising Hamiltonian, which encodes the solution to a QUBO problem. However, their analog nature and their energy-minimization process based on gradient descent make PO-based IMs inherently susceptible to identifying inaccurate solutions. In this work, we introduce and validate a QUBO solver—the analog Floquet solver (AFS)—which enhances the dynamics of PO-based IMs by leveraging Floquet states that emerge spontaneously in POs coupled to high-quality-factor resonances. These states enable the AFS to embed periodic time modulation into its energy-minimization process, allowing it to escape local minima during the search for QUBO problem solutions. As a result, the AFS significantly increases the likelihood of identifying accurate solutions compared to conventional PO-based IMs. More generally, this work defines a new paradigm in analog computing—spanning both classical and quantum realms—that can be physically realized with existing technologies across diverse physical domains.
Solving large-scale quadratic unconstrained binary optimization (QUBO) problems is critical in various fields, including physics, finance, and engineering. However, these problems remain intractable on conventional computing architectures. Alternative solvers, such as Ising machines (IMs) based on networks of coupled electronic, mechanical, or photonic parametric oscillators (POs), have recently been developed. PO-based IMs aim to find the ground state of an Ising Hamiltonian, which encodes the solution to a QUBO problem. However, their analog nature and their energy-minimization process based on gradient descent make PO-based IMs inherently susceptible to identifying inaccurate solutions. In this work, we introduce and validate a QUBO solver—the analog Floquet solver (AFS)—which enhances the dynamics of PO-based IMs by leveraging Floquet states that emerge spontaneously in POs coupled to high-quality-factor resonances. These states enable the AFS to embed periodic time modulation into its energy-minimization process, allowing it to escape local minima during the search for QUBO problem solutions. As a result, the AFS significantly increases the likelihood of identifying accurate solutions compared to conventional PO-based IMs. More generally, this work defines a new paradigm in analog computing—spanning both classical and quantum realms—that can be physically realized with existing technologies across diverse physical domains.
Abstract Solving large-scale Quadratic Unconstrained Binary Optimization (QUBO) problems is critical in various fields, including finance, physics, biology and engineering. However, these problems remain intractable on conventional computing architectures. As a result, alternative solvers—such as Ising Machines (IMs) based on networks of coupled electronic, mechanical or photonic parametric oscillators (POs)—have recently been developed. PO-based Ising Machines (IMs) aim to find the ground state of an Ising Hamiltonian, which encodes the solution to a QUBO problem. However, they rely on an energy minimization process based on gradient descent, making them inherently susceptible to getting trapped in local minima and, as a result, to identifying inaccurate solutions. In this work, we introduce and validate a QUBO solver—the Analog Floquet Solver (AFS)—which enhances the dynamics of PO-based Ising machines (IMs) by leveraging Floquet states that emerge spontaneously in POs coupled to high quality-factor resonances. These states enable the AFS to embed periodic time modulation into its energy minimization process, allowing it to escape local minima during the search for QUBO problems’ solutions. As a result, the AFS significantly increases the likelihood of identifying accurate solutions compared to conventional PO-based IMs. More generally, this work establishes a new paradigm in analog computing that can be physically implemented using existing technologies across a variety of physical domains.
With the growing demand for efficient computing to tackle complex real-world challenges like combinatorial optimization, both academia and industry are exploring alternative computing paradigms. The Ising machine (IM) is one such approach, solving optimization problems by mapping their cost functions onto the Ising model. In this framework, finding the optimal solution equates to determining the Ising spin configuration corresponding to the ground state of the Ising Hamiltonian. Analog IMs, relying on artificial spins that can take analogue rather than binary values, have shown great potential in reaching this ground at high speeds. Optical implementations, such as the opto-electronic IM developed within our research group, are particularly promising due to their higher bandwidth, energy efficiency, and lower latency compared to electronic IMs. However, scaling remains a challenge for large problems. To overcome this, we propose a novel analog optical IM utilizing a spatial light modulator (SLM). Our design differs from other IMs relying on SLMs by employing analog spin variables instead of binary ones and leveraging gradient descent-based feedback rather than a computer-driven Metropolis–Hastings sampling process. By exploiting the inherent parallelism of light, this SLM-based IM is expected to efficiently handle optimization problems with millions of spins while still maintaining high energy efficiency. At the conference, we will also offer more insight into the performance of these analogue IMs and will delve deeper in the bit resolution required for the different optical components such as the employed SLM. Based on the dynamical evolution of the IM, we will discuss strategies to further improve the problem-solving capacity and speed of these optical IMs.
An Ising machine (IM), as a type of analog computer tailored for tackling intractable combinatorial optimization problems, has attracted remarkable attention in recent years. In contrast to the blossoming field of bespoke IM hardware, developing metaheuristics from IMs remains largely uninvestigated. Here, we propose a physics-inspired evolutionary computation paradigm, termed the Ising optimization framework (IOF); it comprises a unique Ising algorithm and a hybrid annealing scheme, which together are well-suited for solving quadratic unconstrained binary optimization (QUBO) problems embedded in Ising system energy. The Ising algorithm leverages a set of iterated self-mapping functions to evolve an Ising-spin swarm, enabling efficient energy minimization in artificial Ising systems while mitigating detrimental chaos. Complementing the algorithm, a hybrid annealing scheme integrating singular value dropout, bifurcation control, and a nudging strategy, is devised to augment the overall optimization capacity. The effectiveness of the IOF is validated on various Ising and Max-cut problems with decision variables ranging from 625 to 5000 in number. In comparison to four major types of methods for solving QUBOs, including IM simulations, nature-inspired algorithms, a state-of-the-art heuristic, and the commercial solver Gurobi, the IOF consistently demonstrates notable optimization quality and computational efficiency. This paper provides a theoretical foundation and practical guidelines for bridging Ising-inspired approaches with evolutionary computation, offering an evolutionary perspective on Ising optimizations and suggesting a fertile avenue for future research and application.
We present a heuristic algorithm designed to solve Quadratic Unconstrained Binary Optimization (QUBO) problems efficiently. The algorithm, referred to as IC-D2S, leverages a hybrid approach using Ising and classical machines to address very large problem sizes. Considering the practical limitation on the size of the Ising machine(IM), our algorithm partitions the QUBO problem into a collection of QUBO subproblems (called subQUBOs) and utilizes the IM to solve each subQUBO. Our proposed heuristic algorithm uses a set of control parameters to generate the subQUBOs and explore the search space. Also, it utilizes an annealer based on cosine waveform and applies a mutation operator at each step of the search to diversify the solution space and facilitate the process of finding the global minimum of the problem. We have evaluated the effectiveness of our IC-D2S algorithm on three large-sized problem sets and compared its efficiency in finding the (near-)optimal solution with three QUBO solvers. One of the solvers is a software-based algorithm (D2TS), while the other one (D-Wave) employs a similar approach to ours, utilizing both classical and Ising machines. The results demonstrate that for large-sized problems (>= 5000) the proposed algorithm identifies superior solutions. Additionally, for smaller-sized problems (= 2500), IC-D2S efficiently finds the optimal solution in a significantly faster manner.
We present a heuristic algorithm designed to solve Quadratic Unconstrained Binary Optimization (QUBO) problems efficiently. The algorithm, referred to as IC-D2S, leverages a hybrid approach using Ising and classical machines to address very large problem sizes. Considering the practical limitation on the size of the Ising machine (IM), our algorithm partitions the QUBO problem into a collection of QUBO subproblems (called subQUBOs) and utilizes the IM to solve each subQUBO. Our proposed heuristic algorithm uses a set of control parameters to generate the subQUBOs and explore the search space. Also, it utilizes an annealer based on cosine waveform and applies a mutation operator at each step of the search to diversify the solution space and facilitate the process of finding the global minimum of the problem. We have evaluated the effectiveness of our IC-D2S algorithm on three large-sized problem sets and compared its efficiency in finding the (near-)optimal solution with three QUBO solvers. One of the solvers is a software-based algorithm (D2TS), while the other one (D-Wave) employs a similar approach to ours, utilizing both classical and Ising machines. The results demonstrate that for large-sized problems (≥ 5000) the proposed algorithm identifies superior solutions. Additionally, for smaller-sized problems (= 2500), IC-D2S efficiently finds the optimal solution in a significantly faster manner.
Solving optimization problems is challenging for existing digital computers and even for future quantum hardware. The practical importance of diverse problems, from healthcare to financial optimization, has driven the emergence of specialised hardware over the past decade. However, their support for problems with only binary variables severely restricts the scope of practical problems that can be efficiently embedded. We build analog iterative machine (AIM), the first instance of an opto-electronic solver that natively implements a wider class of quadratic unconstrained mixed optimization (QUMO) problems and supports all-to-all connectivity of both continuous and binary variables.Beyond synthetic 7-bit problems at small-scale, AIM solves the financial transaction settlement problem entirely in analog domain with higher accuracy than quantum hardware and at room temperature. With compute-in-memory operation and spatial-division multiplexed representation of variables, the design of AIM paves the path to chip-scale architecture with 100 times speed-up per unit-power over the latest GPUs for solving problems with 10,000 variables. The robustness of the AIM algorithm at such scale is further demonstrated by comparing it with commercial production solvers across multiple benchmarks, where for several problems we report new best solutions. By combining the superior QUMO abstraction, sophisticated gradient descent methods inspired by machine learning, and commodity hardware, AIM introduces a novel platform with a step change in expressiveness, performance, and scalability, for optimization in the post-Moores law era.