2016/04/18

-An important step toward photonic computation that solves intractable mathematical problems-

=> Press Release
=> Quantum Optical State Control Research Group

Background

In our society, systems, such as the internet and traffic and social networks, are becoming more and more complex. Therefore, optimizing them is important for efficient use of natural and social resources. Many such problems are reduced to mathematical problems called combinatorial optimization problems, where we find an optimum combination from a finite set of combination. Such problems are known to be very hard to solve with modern digital computers, because the computation time ‘explodes’ as the number of combinations increases. On the other hand, combinatorial optimization problems can be mapped onto the ground-state-search problems of the Ising model, which is a theoretical model for interacting spins. Today, several institutions are trying to solve the Ising model by using artificial spin networks.

Among 'Ising machines', the CIM is now drawing attention as a way to solve the Ising model efficiently by utilizing artificial spins realized with photonics technologies (Fig. 1and Fig. 2). In the CIM, a group of OPOs is used as spins. An OPO is a special kind of laser oscillator that takes only the 0 or π phase and can therefore be used as a stable artificial spin. The spin-spin interactions can be implemented by mutual injection of OPO lights. With sufficiently slow pumping of the oscillator network, the OPO network tends to start oscillating with the combination of phases that gives the loss-minimum of the whole network, which is the ground state of the corresponding spin network.

In 2014, a group from Stanford University performed a proof-of-principle CIM experiment using four OPOs. However, in order to solve complex combinatorial optimization problems, we need to increase the number of artificial spins.

Achievements

We have succeeded in generating more than10,000 time-multiplexed OPOs using an optical fiber cavity as long as 1 km (Fig. 3). This was accomplished by utilizing technologies developed through the R&D of optical fiber communications. The OPOs can be used as artificial spins for a large-scale CIM to solve combinatorial optimization problems in the real world. In addition, we have realized a one-dimensional Ising spin network by introducing nearest-neighbor coupling between OPOs, and found that our room-temperature OPOs can simulate the behavior of low-temperature spins.

Generation of large-scale time-multiplexed OPOs

1. We realized a phase-sensitive amplifier (PSA), which is an optical amplifier that amplifies only light with the 0 or π phase, using four-wave mixing in a 1-km HNLF.
2. By placing the PSA in a fiber loop, we constructed a 1-km fiber cavity. We then operated the PSA with 2-GHz clock frequency to obtain a total 10,320 OPOs from the single cavity (Fig. 3).
3. We experimentally confirmed that the OPO phases were discretized to 0 or π, which suggests that they can be used as stable artificial spins.

Simulation of one-dimensional spin network

1. We realized one-dimensional spin coupling by inserting a 1-bit delay interferometer into the fiber cavity (Fig. 5).
2. By adjusting the OPO coupling phase, we observed that the spin network can behave like ferro- and anti-ferromagnetic spins. In addition, we confirmed that we could decrease the total network energy by changing the condition of the pump for the OPOs. As a result, we were able to confirm that the OPOs can be used as high-quality artificial spins.

Technical Features

1. Large-scale time-multiplexed OPOs
   Figure 3 shows the experimental setup. The 100-ps width, 2-GHz clock laser pulses with wavelengths of 1531 and 1551 nm were injected into a fiber cavity through a wavelength-division multiplexing (WDM) filter. The cavity contained a WDM filter, a 1-km HNLF, a narrowband wavelength filter, a polarization controller, and a 99:1 coupler. Through the four-wave mixing in the HNLF, we can obtain phase-sensitive amplification at1541-nm wavelength, which corresponds to the midpoint between the two pump wavelengths. The noise lights generated from the PSA circulate in the cavity, and eventually lead to the generation of an OPO in which only the 0- or π-phase light can oscillate. Since the cavity round-trip time was around 5 µs and the pulse interval was 500 ps, we were able to generate more than 10,000 OPOs using a single cavity. A portion of the OPOs were output from the cavity through the 99:1 coupler, and the phase differences between adjacent OPOs were measured using a 1-bit delay interferometer followed by photodetectors. A portion of phase difference measurement results is shown in Fig. 4(a), where the cosine components of the phase differences were discretized to +1 or -1, implying that the OPOs had only the 0 or π phase. A histogram of the phase difference measurement results for all of the OPOs (total number: 10,320), is shown in Fig. 4(b). We observed a clear bifurcation of OPO phases to 0 and π, which confirmed the generation of more than 10,000 stable artificial spins.
2. Simulation of one-dimensional spin network
   We performed an experiment to simulate the one-dimensional Ising model using the OPOs. We inserted another 1-bit delay interferometer into the fiber cavity so that we could inject a portion of an OPO light to an adjacent OPO to implement nearest neighbor spin couplings (Fig. 5). We can realize ferro- and anti-ferromagnetic spin coupling by adjusting the phase difference between the two arms of the interferometer to 0 or π. The phase difference measurement results are shown in Fig. 6 (a) (coupling phase 0) and (b) (coupling phase π). From these results, we confirmed that the spin network showed ferro- and anti-ferromagnetic behavior depending on the coupling phase. In addition, we found that not all of the spins were aligned. Instead, we observed probabilistic spin flips. As a result, domain formation was observed as shown in Fig. 6 (a), for example. According to the theory of the one-dimensional Ising model, the ground state is not obtained and domains will be formed at finite spin temperature. This characteristic was well simulated by our experiment. We also investigated how the domain lengths changed as we varied the pump amplitude in the case of ferromagnetic coupling. As shown in Fig. 7, the closer the pump amplitude was set to oscillation threshold, the larger the domains became. This suggests that the spin state will get closer to the ground state (i.e. all the spins will be aligned in the same direction) when the pump amplitude is set near the threshold.

Future Plans

In the present experiment, since we implemented only nearest-neighbor coupling, only a one-dimensional spin network could be simulated. We are now planning to implement more flexible coupling between OPOs to simulate more complex spin networks, so that we can realize a CIM that solves large-scale combinatorial optimization problems.

Publication

"Large-scale Ising spin network based on degenerate optical parametric oscillators"
Takahiro Inagaki, Kensuke Inaba, Ryan Hamerly, Kyo Inoue, Yoshihisa Yamamoto & Hiroki Takesue
Nature Photonics 2016