Coherent Ising Machines (CIMs) have emerged as promising hardware solvers for NP-hard combinatorial optimization problems, potentially offering energy-efficiency, accuracy, speed, and scalability. These problems span applications including scheduling, routing, resource allocation, and machine learning. Spin amplitude inhomogeneity (an artifact of its analog nature) can hinder CIM convergence to the ground state. Several CIM modifications have been proposed to address this issue. However, no prior work has systematically compared the performance of these concepts. This paper presents the first side-by-side evaluation of six promising CIM architectures using a unified simulation model. We develop stochastic differential equations, find optimal hyperparameters via Bayesian optimization per problem-architecture pair, and assess the performance on multiple metrics using MaxCut instances from the Gset and BiqMac datasets. CIM-SNN outperforms the rest on most metrics, with a 54.21% average success probability (SP) improvement and 1.36× less time-to-solution (TTS) at 96% GS energy compared to the Standard CIM on BiqMac graphs. Other CIMs show some advantages on certain Gset graphs. CIM-QA is usually the worst-performing, with an 18.55% average SP loss and a 1.30× higher TTS at 96% GS compared to the Standard CIM on BiqMac graphs. CIM-SNN demonstrates the greatest noise robustness and promise for further tuning via a linear pump schedule. CIM-QA benefits from high noise levels. This study makes two key contributions: (1) identifying promising CIM architectures to guide targeted research and accelerate practical adoption, and (2) providing a unified framework for evaluating future architectural modifications against existing ones.
Coherent Ising Machines (CIMs) have emerged as promising hardware solvers for NP-hard combinatorial optimization problems, potentially offering energy-efficiency, accuracy, speed, and scalability. These problems span applications including scheduling, routing, resource allocation, and machine learning. Spin amplitude inhomogeneity (an artifact of its analog nature) can hinder CIM convergence to the ground state. Several CIM modifications have been proposed to address this issue. However, no prior work has systematically compared the performance of these concepts. This paper presents the first side-by-side evaluation of six promising CIM architectures using a unified simulation model. We develop stochastic differential equations, find optimal hyperparameters via Bayesian optimization per problem-architecture pair, and assess the performance on multiple metrics using MaxCut instances from the Gset and BiqMac datasets. CIM-SNN outperforms the rest on most metrics, with a 54.21% average success probability (SP) improvement and 1.36× less time-to-solution (TTS) at 96% GS energy compared to the Standard CIM on BiqMac graphs. Other CIMs show some advantages on certain Gset graphs. CIM-QA is usually the worst-performing, with an 18.55% average SP loss and a 1.30× higher TTS at 96% GS compared to the Standard CIM on BiqMac graphs. CIM-SNN demonstrates the greatest noise robustness and promise for further tuning via a linear pump schedule. CIM-QA benefits from high noise levels. This study makes two key contributions: (1) identifying promising CIM architectures to guide targeted research and accelerate practical adoption, and (2) providing a unified framework for evaluating future architectural modifications against existing ones.
Coherent Ising Machines (CIMs) are emerging hardware solvers for approximating solutions to combinatorial optimization problems, and may offer advantages in energy-efficiency, speed, and scalability across applications. However, CIM computing is a nascent and methodologically fragmented field, and substantial research is required to determine its utility for real-world combinatorial optimization problems. The core argument of this thesis is that several complementary efforts are required: (1) developing a unified simulations framework that systematically evaluates CIM variations by isolating implementation and methodological differences from architectural ones; (2) introducing layers of abstraction, analogous to those in mature CMOS-based computing, to enable effort (1); and (3), by extension of efforts (1) and (2), potentially enabling CIM-inspired dynamics -- which, when simulated, are simply classical algorithms that run on conventional computers -- for efficient in-silico approximation of combinatorial optimization problems.As a concrete example of the utility of the unified simulations framework, this thesis compares several CIM modifications proposed to address spin amplitude inhomogeneity (an artifact of the CIM's analog nature) that can hinder convergence to the ground state. These variations have been analyzed under varying assumptions and methodologies, and no prior work has systematically compared their performance by isolating architectural differences. This thesis presents the first side-by-side evaluation of six promising CIM architectures using a unified simulation model. We develop implementation-agnostic stochastic differential equations, find optimal hyperparameters via Bayesian optimization per problem-architecture pair, and assess the performance on multiple metrics using MaxCut instances from the Gset and BiqMac datasets. CIM-SNN outperforms the rest on most metrics, with a 54.21% average success probability (SP) improvement and 1.36× faster time-to-solution (TTS) at 96% GS energy compared to the Standard CIM on BiqMac graphs. Other CIMs show some advantages on certain Gset graphs. CIM-QA is usually the worst-performing, with an 18.55% average SP loss and a 1.30× slower TTS at 96% GS compared to the Standard CIM on BiqMac graphs. CIM-SNN demonstrates the greatest noise robustness and promise for further tuning via a linear pump schedule. CIM-QA benefits from high noise levels. These analyses were made possible only by the unified simulations framework developed in this thesis, which can also be used to evaluate future architectural modifications against existing ones. Lastly, we present a potential abstraction stack that further extends the unified simulations framework, which can benefit future research and introduce a new direction in the field.
Coherent Ising machines (CIMs) have emerged as a hybrid form of quantum computing devices designed to solve NP-complete problems, offering an exciting opportunity for discovering optimal solutions. Despite challenges such as susceptibility to noise-induced local minima, we achieved notable advantages in improving the computational accuracy and stability of CIMs. We conducted a successful experimental demonstration of CIM via femtosecond laser pumping that integrates optimization strategies across optical and structural dimensions, resulting in significant performance enhancements. The results are particularly promising. An average success rate of 55% was achieved to identify optimal solutions within a Möbius Ladder graph comprising 100 vertices. Compared with other alternatives, the femtosecond pulse results in significantly higher peak power, leading to more pronounced quantum effects and lower pump power in optical fiber-based CIMs. In addition, we have maintained an impressive success rate for a continuous period of 8 hours, emphasizing the practical applicability of CIMs in real-world scenarios. Furthermore, our research extends to the application of these principles in practical applications such as molecular docking and credit scoring. The results presented substantiate the theoretical promise of CIMs, paving the way for their integration into large-scale practical applications.
Coherent Ising Machines (CIMs) have emerged as a hybrid form of quantum computing devices designed to solve NP-complete problems, offering an exciting opportunity for discovering optimal solutions. Despite challenges such as susceptibility to noise-induced local minima, we achieved notable advantages in improving the computational accuracy and stability of CIMs. We conducted a successful experimental demonstration of CIM via femto-second laser pumping that integrates optimization strategies across optical and structural dimensions, resulting in significant performance enhancements. The results are particularly promising. An average success rate of 55% was achieved to identify optimal solutions within a Mobius Ladder graph comprising 100 vertices. Compared with other alternatives, the femto-second pulse results in significantly higher peak power, leading to more pronounced quantum effects and lower pump power in optical fiber based CIMs. In addition, we have maintained an impressive success rate for a continuous period of 8 hours, emphasizing the practical applicability of CIMs in real-world scenarios. Furthermore, our research extends to the application of these principles in practical applications such as molecular docking and credit scoring. The results presented substantiate the theoretical promise of CIMs, paving the way for their integration into large-scale practical applications.
Coherent Ising machines (CIMs), leveraging the bistable physical properties of coherent light to emulate Ising spins, exhibit great potential as hardware accelerators for tackling complex combinatorial optimization problems. Recent advances have demonstrated that the performance of CIMs can be enhanced either by incorporating large random noise or higher-order nonlinearities, yet their combined effects on CIM performance remain mainly unexplored. In this work, we develop a numerical CIM model that utilizes a tunable fifth-order polynomial nonlinear dynamic function under large noise levels, which has the potential to be implemented in all-optical platforms. We propose a normal form of a CIM model that allows for both supercritical and subcritical pitchfork bifurcation operational regimes, with fifth-order nonlinearity and tunable hyperparameters to control the Ising spin dynamics. In the benchmark studies, we simulate various sets of MaxCut problems using our fifth-order polynomial CIM model. The results show a significant performance improvement, achieving an average of 59.5% improvement in median time-to-solution (TTS) and an average of 6 times improvement in median success rate (SR) for dense Maxcut problems in the BiqMac library, compared to the commonly used third-order polynomial CIM model with low noise. The fifth-order polynomial CIM model in the large-noise regime also shows better performance trends as the problem size scales up. These findings reveal the enhancements on the computational performance of Ising machines in the large-nose regime from fifth-order nonlinearity, showing important implications for both simulation and hardware perspectives.
The coherent Ising machine (CIM) is a network of optical parametric oscillators (OPOs) that solves for the ground state of Ising problems through OPO bifurcation dynamics. Here, we present experimental results comparing the performance of the CIM to quantum annealers (QAs) on two classes of NP-hard optimization problems: ground state calculation of the Sherrington-Kirkpatrick (SK) model and MAX-CUT. While the two machines perform comparably on sparsely-connected problems such as cubic MAX-CUT, on problems with dense connectivity, the QA shows an exponential performance penalty relative to CIMs. We attribute this to the embedding overhead required to map dense problems onto the sparse hardware architecture of the QA, a problem that can be overcome in photonic architectures such as the CIM.
Coherent Ising Machines (CIMs) have recently gained attention as a promising computing model for solving combinatorial optimization problems. In particular, the Chaotic Amplitude Control (CAC) algorithm has demonstrated high solution quality, but its performance is highly sensitive to a large number of hyperparameters, making efficient tuning essential. In this study, we present an algorithm portfolio approach for hyperparameter tuning in CIMs employing Chaotic Amplitude Control with momentum (CACm) algorithm. Our method incorporates multiple search strategies, enabling flexible and effective adaptation to the characteristics of the hyperparameter space. Specifically, we propose two representative tuning methods, Method A and Method B. Method A optimizes each hyperparameter sequentially with a fixed total number of trials, while Method B prioritizes hyperparameters based on initial evaluations before applying Method A in order. Performance evaluations were conducted on the Supercomputer “Flow” at Nagoya University, using planted Wishart instances and Time to Solution (TTS) as the evaluation metric. Compared to the baseline performance with best-known hyperparameters, Method A achieved up to 1.47× improvement, and Method B achieved up to 1.65× improvement. These results demonstrate the effectiveness of the algorithm portfolio approach in enhancing the tuning process for CIMs.