特許 7352916 イジングモデルの計算装置

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】

(24)【登録日】2023-09-21

(45)【発行日】2023-09-29

(54)【発明の名称】イジングモデルの計算装置

(51)【国際特許分類】

   G06N 99/00 20190101AFI20230922BHJP

   G06E 3/00 20060101ALI20230922BHJP

   G02F 3/00 20060101ALI20230922BHJP

   G02F 1/39 20060101ALI20230922BHJP

【ＦＩ】

G06N99/00 180

G06E3/00

G02F3/00 501

G02F1/39

【請求項の数】 8

(21)【出願番号】P 2022511161

(86)(22)【出願日】2021-04-02

(86)【国際出願番号】 JP2021014330

(87)【国際公開番号】W WO2021201279

(87)【国際公開日】2021-10-07

【審査請求日】2022-08-22

(31)【優先権主張番号】P 2020066491

(32)【優先日】2020-04-02

(33)【優先権主張国・地域又は機関】JP

【国等の委託研究の成果に係る記載事項】（出願人による申告）平成２６年度、国立研究開発法人 科学技術振興機構、革新的研究開発推進プログラム（ＩｍＰＡＣＴ）「大規模時分割多重光パラメトリック発振器」受託研究、産業技術力強化法第１７条の適用を受ける特許出願

(73)【特許権者】

【識別番号】000004226

【氏名又は名称】日本電信電話株式会社

(73)【特許権者】

【識別番号】504202472

【氏名又は名称】大学共同利用機関法人情報・システム研究機構

(74)【代理人】

【識別番号】110001243

【氏名又は名称】弁理士法人谷・阿部特許事務所

(72)【発明者】

【氏名】武居 弘樹

(72)【発明者】

【氏名】稲垣 卓弘

(72)【発明者】

【氏名】本庄 利守

(72)【発明者】

【氏名】稲葉 謙介

(72)【発明者】

【氏名】宇都宮 聖子

(72)【発明者】

【氏名】玉手 修平

【審査官】稲葉 崇

(56)【参考文献】

【文献】特開２０２０－３８３００（ＪＰ，Ａ）

【文献】国際公開第２０１７／０４７６６６（ＷＯ，Ａ１）

【文献】特開２０１７－７３１０６（ＪＰ，Ａ）

(58)【調査した分野】(Int.Cl.，ＤＢ名)

Ｇ０６Ｎ ３／００－９９／００

Ｇ０６Ｆ １８／００－１８／４０

Ｇ０６Ｅ １／００－３／００

Ｇ０２Ｆ １／１５－１／１９

Ｇ０２Ｆ １／００－１／１２５，１／２１－７／００

(57)【特許請求の範囲】

【請求項1】

イジングモデルの複数のスピンに対応し、同一の発振周波数を有する複数の光パルスを０またはπの位相でパラメトリック発振させる位相感応増幅器と、
前記位相感応増幅器から出力された前記複数の光パルスの位相および振幅を測定する、光パルス測定部と、
前記光パルス測定部において測定した光パルスの位相および振幅の情報を入力として、前記イジングモデルの結合係数と前記測定した光パルスとから決定される、ある光パルスに関わる相互作用に基づいて、フィードバック値を算出し、フィードバック信号を出力する、高速演算回路と、
前記高速演算回路において算出された前記フィードバック信号に基づいて、入力された前記複数の光パルスと同じ数の複数の光パルスの位相および振幅を変調することにより、前記ある光パルスに関わる相互作用を実装する変調器と、
を備え、
前記変調器は、ｐｕｓｈ－ｐｕｌｌ変調器であって、マッハツェンダ干渉計の両アーム導波路に配置された２つの位相変調器を絶対値が同じで符号が異なる信号で駆動し、前記変調器の動作点を透過する光強度が最小となる点で駆動するように構成され、
前記フィードバック信号を重畳する際に、前記動作点を、印加した磁場に相当する電圧分シフトさせることにより、光共振器のＤＯＰＯと同じ波長で、かつポンプ光との初期位相差が固定されているＤＯＰＯのパルスに定数の光が注入されるように構成されていることを特徴とする、イジングモデルの計算装置。

【請求項2】

前記動作点の電圧のシフトとして、定数電圧を用いるように構成される請求項１に記載のイジングモデルの計算装置。

【請求項3】

前記動作点の電圧のシフトとして、あらかじめ定められた関数に従って電圧を変化するように構成される請求項１に記載のイジングモデルの計算装置。

【請求項4】

前記電圧の変化は、線形増加である請求項３に記載のイジングモデルの計算装置。

【請求項5】

i番目（iは、自然数）のＤＯＰＯにおける前記動作点の電圧のシフトとして、i番目のパルスの測定の結果得られたＤＯＰＯ振幅c_iに比例する電圧を用いるように構成される請求項１に記載のイジングモデルの計算装置。

【請求項6】

i番目（iは、自然数）のＤＯＰＯにおける前記動作点の電圧のシフトとして、i番目のパルスの測定の結果得られたＤＯＰＯ振幅の絶対値c_iに比例する電圧を用いるように構成される請求項１に記載のイジングモデルの計算装置。

【請求項7】

i番目（ｎは、自然数）のＤＯＰＯにおける前記動作点の電圧のシフトとして、Ｎを全パルス数、ＤＯＰＯ振幅をc_iとすると、
全てのパルスの測定の結果得られたＤＯＰＯ振幅の絶対値の平均

【数1】

に比例する電圧を用いるように構成される請求項１に記載のイジングモデルの計算装置。

【請求項8】

前記動作点の電圧のシフトとして、定数電圧を用いる第１の構成、
前記動作点の電圧のシフトとして、あらかじめ定められた関数に従って電圧を変化する第２の構成、
前記動作点の電圧のシフトとして、あらかじめ定められた関数に従って電圧を変化し、前記電圧の変化は、線形増加である第３の構成、
i番目（iは、自然数）のＤＯＰＯにおける前記動作点の電圧のシフトとして、i番目のパルスの測定の結果得られたＤＯＰＯ振幅c_iに比例する電圧を用いるように構成される第４の構成、及び
i番目（iは、自然数）のＤＯＰＯにおける前記動作点の電圧のシフトとして、i番目のパルスの測定の結果得られたＤＯＰＯ振幅の絶対値c_iに比例する電圧を用いるように構成される第５の構成のうちの複数個の構成を同時に組み合わせて下式を満たすように、動作するように構成される請求項１に記載のイジングモデルの計算装置。

【数2】

A, B, C, D, 及びEは任意の実数

【発明の詳細な説明】

【技術分野】

【0001】

本開示は、計算機や光システムに関し、特に、組み合わせ最適化問題を解くイジングモデルの計算装置に関する技術に関する。

【背景技術】

【0002】

配送ルートの最適化や無線周波数割当など、現代社会の多くの重要な問題は組合せ最適化問題であり、大規模サイズの最適化問題をできるだけ短時間で、できるだけ高精度に求めることが急務となっている。２０１１年にレーザーネットワークを用いてイジング問題を解く「コヒーレント・イジングマシン」が提案され、その後の縮退光パラメトリック発振器（Ｄｅｇｅｎｅｒａｔｅ Ｏｐｔｉｃａｌ Ｐａｒａｍｅｔｒｉｃ Ｏｓｃｉｌｌａｔｏｒ：以下、ＤＯＰＯ）ネットワーク用いた時分割多重方式の提案により、大規模化が見込めるようになった。[非特許文献１]

【0003】

ＤＯＰＯパルスの位相０, πでスピン値を表現し、ＤＯＰＯパルス群を測定・フィードバックにより相互作用させることで相互作用するスピン群の理論モデル(イジングモデル)の基底状態探索問題を解くコヒーレントイジングマシン(Ｃｏｈｅｒｅｎｔ Ｉｓｉｎｇ Ｍａｃｈｉｎｅ)が開発されている［非特許文献２］。

【0004】

このコヒーレントイジングマシンに実装されているＤＯＰＯパルス間の相互作用を表現するハミルトニアンは式(１)で表わされる。J_ijはスピン間相互作用係数、

【0005】

【数1】

【0006】

イジングモデルでは｛１,-１｝の値をとるスピン σ_i を、コヒーレントイジングマシンでは、正負のアナログ値をとり、その絶対値が時間発展に伴い飽和する振幅の余弦成分 ｃｉ を用いて近似的に表現する。また、演算途中のスピンを測定して、その情報をもとにスピン間相互作用を入力するための信号をＦＰＧＡ(Ｆｉｅｌｄ-Ｐｒｏｇｒａｍｍａｂｌｅ Ｇａｔｅ Ａｒｒａｙ)を用いて計算し、その信号を光に乗せてフィードバックする仕組みによって、収束するまで繰り返し計算を行っている。このコヒーレントイジングマシンを用いて、最大カット問題とよばれる組合せ最適化問題の高速な解探索が可能であることが報告されている。

【0007】

ＤＯＰＯを用いて式（１）の形で実装するとき、i番目（iは自然数）のＤＯＰＯ振幅の時間発展は以下の式（１ｂ）で近似的に記述される[非特許文献３]。

【0008】

【数2】

【0009】

ここで、ｐは独立ＤＯＰＯの発振閾値での値で規格化したポンプ振幅、ｃ_iはｐ=２における振幅値で規格化したＤＯＰＯ振幅の余弦成分である。ＤＯＰＯの性質として振幅の正弦成分は減衰することから無視した。

【0010】

一方で、より一般化されたイジングモデルは、式（２）のハミルトニアンで表される。

【0011】

【数3】

【0012】

ここでJ_ijはスピン間相互作用係数、σ_iはサイトiにおけるスピン、ｈ_iはサイトiにおける局所磁場である。式（２）の右辺第一項は、スピン間相互作用項であり、右辺第二項は、各スピンに対する磁場項である。地図の四色塗り分け問題、巡回セールスマン問題などより多種類の組合せ最適化問題を式(２)のハミルトニアンに変換して解を求められることが知られている。

【先行技術文献】

【非特許文献】

【0013】

【文献】https://qistokyo.wordpress.com/research/coherent-ising-machine/

【文献】T. Inagaki et al., Science, 2016, Vol. 354, pp. 603-606

【文献】Z. Wang et al., Phys. Rev., 2013, A 88, 063853-1-9

【文献】A. Lucus, Frontiers in physics, 2014, ｖol.2, article 5, pp.1-15

【文献】W. H. Steier, IEEE J. Quant. Electron. 1967, ｖol.QE-3, pp.664-667

【発明の概要】

【0014】

従来のコヒーレントイジングマシンのハミルトニアン式（１）においては、式（２）左辺の第二項である各スピンに対する磁場項が無いため、地図の四色塗り分け問題、巡回セールスマン問題などの組み合わせ最適化問題が表現できない。よって、従来のコヒーレントイジングマシンで解ける問題が最大カット問題などに限られている、という課題があった。

【0015】

本開示はかかる従来の問題に鑑みなされたものであって、本開示の課題は、コヒーレントイジングマシンのフィードバック信号の生成に用いられてきたｐｕｓｈ－ｐｕｌｌ型の光変調器の動作点を調整することにより、局所磁場(hi)に相当する役割を実装することにある。

【0016】

上述の課題を解決するために、本開示のイジングモデルの計算装置の変調器は、従来のコヒーレントイジングマシンに、磁場項を模擬した機能を追加することにより、各サイトにおける磁場を表現する項を実装することを特徴とする。具体的には従来のコヒーレントイジングマシンにおいて、ｐｕｓｈ－ｐｕｌｌ型光変調器の動作点をずらすことにより、あらかじめ決まった符号の振幅を持つ光を入れる。動作点を＋方向にずらしておけば、その光を注入されたＤＯＰＯは０位相で発振する傾向が強くなり、動作点をマイナス方向にずらしておけば、その逆となる。

【0017】

本開示のイジングモデルの計算装置の態様は、イジングモデルの複数のスピンに対応し、同一の発振周波数を有する複数の光パルスを０またはπの位相でパラメトリック発振させる位相感応増幅器と、前記位相感応増幅器から出力された前記複数の光パルスの位相および振幅を測定する、光パルス測定部と、前記光パルス測定部において測定した光パルスの位相および振幅の情報を入力として、前記イジングモデルの結合係数と前記測定した光パルスとから決定される、ある光パルスに関わる相互作用に基づいて、フィードバック値を算出し、フィードバック信号を出力する、高速演算回路と、前記高速演算回路において算出された前記フィードバック信号に基づいて、入力された前記複数の光パルスと同じ数の複数の光パルスの位相および振幅を変調することにより、前記ある光パルスに関わる相互作用を実装する変調器と、を備え、前記変調器は、ｐｕｓｈ－ｐｕｌｌ変調器であって、マッハツェンダ干渉計の両アーム導波路に配置された２つの位相変調器を絶対値が同じで符号が異なる信号で駆動し、前記変調器の動作点を透過する光強度が最小となる点で駆動するように構成され、前記フィードバック信号を重畳する際に、前記動作点を、印加した磁場に相当する電圧分シフトさせることにより、光共振器のＤＯＰＯと同じ波長で、かつポンプ光との初期位相差が固定されているＤＯＰＯのパルスに定数の光が注入されるように構成されていることを特徴とする。

【0018】

本開示によって、コヒーレントイジングマシンにおいて、光変調器の動作点を調整することにより、局所磁場項を実装して、磁場項を含む一般的なイジングモデルを表現できるようにして、多様な組み合わせ最適化問題の解が求められるようになる。

【図面の簡単な説明】

【0019】

【図1】本発明の実施例にかかるイジングモデルの計算装置の構成を示す図である。

【図2】ｐｕｓｈ－ｐｕｌｌ変調器を示す図である。

【図3】光変調器への印加電圧と、光変調器の透過率との関係を示す図である。

【図4】（ａ）地図の４色塗り分け問題の解法を示す図である。（ｂ）地図の４色塗り分け問題の解法を示す図である。

【図5】本発明の実施例にかかる実験結果（地図の４色塗り分け問題）を示す図である。

【発明を実施するための形態】

【0020】

以下、本開示の光変調器を備えたイジングモデルの計算装置の形態は、実施例及び図面を用いて詳細に説明される。但し、本発明は以下に示す実施の形態の記載内容に限定されず、本明細書等において開示する発明の趣旨から逸脱することなく形態および詳細を様々に変更し得ることは当業者にとって自明である。

【0021】

以下に、コヒーレントイジングマシンが説明される。

【0022】

縮退光パラメトリック発振器ＤＯＰＯとは、(ある位相基準に対して)位相０またはπのいずれかでしか発振しないレーザである。コヒーレントイジングマシンとは、ＤＯＰＯ間の光の相互注入による同期現象を利用するものである。おなじ波長で発振するＤＯＰＯ１とＤＯＰＯ２の２つのＤＯＰＯがあるとする。ＤＯＰＯ１の光を、位相差０でＤＯＰＯ２に注入すると、光の注入同期現象によりＤＯＰＯ２はＤＯＰＯ１と同位相で発振しようとする。光を位相差πでＤＯＰＯ２に注入すると、逆位相(位相差π)でＤＯＰＯ２は発振しようとする。コヒーレントイジングマシンでは、このような光の注入同期現象を用いてスピン間の相互作用を模擬している。

【0023】

本発明の実施形態のコヒーレントイジングマシンは、上述したスピン間の相互作用を模擬することに加え、磁場項を模擬することにより、上述した式（２）で表されるイジングモデルを実装する。上述した従来のコヒーレントイジングマシンでは、問題としてＤＯＰＯの間の関係性(すなわち、この２つのＤＯＰＯの位相は揃った方が良いとか、逆の方が良いとか)が与えられ、その関係性(複数ある)をもっともよく満たすＤＯＰＯの位相の組み合せを求めるものである。すなわち、ＤＯＰＯ間の位相関係に注目し、各ＤＯＰＯがどちらの位相になるかということに対する条件は無い。これに対し、本実施形態は、各ＤＯＰＯに模擬的に局所磁場をかける。換言すれば、「磁場をかける」とは、各ＤＯＰＯに対して「どちらかというと０位相で発振せよ」、又は「どちらかというとπ位相で発振せよ」という条件を加えることに相当し、具体的には、各ＤＯＰＯに対し、位相基準に対して「位相差０の光」、又は「位相差πの光」を注入することに相当する。これらの注入される光は相互注入の他のＤＯＰＯからのものではなく、DOPOで用いるポンプ光によるレーザを「位相基準」としたものである。

【0024】

ＤＯＰＯパルスは位相感応増幅器と呼ばれる光増幅器を光共振器中に配置し、この位相感応増幅器にポンプ光を入力することで生成される。位相感応増幅器は非線形光学媒質を有して構成される。そして、位相感応増幅器は、その非線形光学媒質によって、光共振器に注入される光のうちポンプ光との初期位相差が０またはπの光をもっともよく増幅する。よって、ＤＯＰＯパルスは、利得の最も高い０位相差またはπ位相差で光が発振する。

【0025】

磁場項を模擬してＤＯＰＯに条件を付加した「測定・フィードバック法」を用いることで、全てのＤＯＰＯの間に任意の相互作用項および磁場項を実装可能なコヒーレントイジグマシンが実現される。これにより、多数の頂点や辺を持つ複雑なグラフも、そのまま解くことができる。

【0026】

［実施例１］
全長１ｋｍの長距離光ファイバーリング共振器中に、Ｎ個（Ｎは自然数）のＤＯＰＯの光パルス群３を一括発生させ、一部のＤＯＰＯパルスを取り出して測定，演算，変調を行ない、再び変調されたＤＯＰＯパルスを共振器に戻す「測定・フィードバック法」が採用される。

【0027】

図１は本実施形態のイジングモデルの計算装置１００の概略構成を示す図である。図１において、イジングモデルの計算装置は、リング状の光ファイバで構成された光共振器１と、光共振器１内に設けられた、位相感応増幅器２と、光共振器１から分岐された、フィードバックループの一部を構成する、光パルス測定部５と高速演算回路（ＦＰＧＡ）６と光変調器（変調器）８とを備えている。本実施形態のイジングモデルの計算装置では、光パルス測定部５と高速演算回路（ＦＰＧＡ）６と光変調器８とがフィードバックループを構成している。光共振器１は、カプラ４ａ, ４ｂを介して、フィードバックループと接続している。

【0028】

イジングモデルの計算装置１００は、図１に示すように、光共振器１として機能するリング状の光ファイバ内に設けられた位相感応増幅器（ＰＳＡ：Ｐｈａｓｅ Ｓｅｎｓｉｔｉｖｅ Ａｍｐｌｉｆｉｅｒ）２に対して、ポンプ光パルス（ｐｕｍｐ）を注入することによりイジングモデルのサイト数に対応する数のＤＯＰＯの光パルス群３を生成するように構成している（２値化ＯＰＯ：Ｏｐｔｉｃａｌ Ｐａｒａｍｅｔｒｉｃ Ｏｓｃｉｌｌａｔｉｏｎ：０またはπ位相の光パラメトリック発振）。光共振器１内のカプラ４ａを介して、ＤＯＰＯの光パルス群３から分岐された光パルス３ａは、光パルス測定部５で、その位相および振幅が測定される。光パルス測定部５において測定した光パルスの位相および振幅の情報を入力として、あらかじめ与えられるイジングモデルの結合係数と測定した光パルス３ａとから決定される、ある光パルスに関わる相互作用に基づいて、高速演算回路（ＦＰＧＡ）６でフィードバック信号７が算出される。測定・フィードバック法に基づくイジングモデルの計算装置（コヒーレントイジングマシン）１００では、高速演算回路（ＦＰＧＡ）６で計算されたフィードバック信号７を、光変調器８を用いて光共振器１中のＤＯＰＯと同じ波長で、かつポンプ光との初期位相差が固定されている光パルスに重畳して元のＤＯＰＯの光パルスに注入することで、ＤＯＰＯ間の結合を行う。光変調器８としては、図２に示すようなｐｕｓｈ-ｐｕｌｌ変調器と呼ばれる、マッハツェンダ干渉計の両アーム導波路に配置された２つの位相変調器を絶対値が同じで符号が異なる信号で駆動する光変調器[非特許文献５]を、その動作点を透過する光強度が最小となる点で駆動する。この時、入力信号が０なら透過率はゼロ、すなわち光は出ない。入力信号が＋なら、入力光に対して正の符号の振幅を持つ光が出力され、－なら負の符号を持つ光が出力される。これは、入力光に対して０またはπの位相変調を行うことに相当する。なお入力信号の大きさにより、出力信号の振幅の大きさも変えることができる。

【0029】

ここで、図３に示すように、横軸が変調器への印加電圧、縦軸が変調器の光透過率とすると、ｐｕｓｈ-ｐｕｌｌ変調器の動作点を矢印方向にずらす（シフトする）ことにより、あらかじめ決まった符号の振幅を持つ光（定数の光）を光共振器のＤＯＰＯと同じ波長で、かつポンプ光との初期位相差が固定されているＤＯＰＯのパルスに入れることができる。例えば動作点を＋方向にずらしておけば、その光を注入されたＤＯＰＯは０位相で発振する傾向が強くなる。動作点のシフト（磁場項）は、磁場印加後の動作点と磁場無の場合の動作点との差になる。すなわち、上述の「磁場」に相当する効果を実装することができる。

【0030】

カプラ４ｂを介して、変調されたＤＯＰＯの光パルス３ｂは、リング共振器１に注入される。以上の手順を繰り返し、Ｎ個のＤＯＰＯのパルスは、ＰＳＡに共振器１周時間時間あたりＮ個のポンプパルスを入力することにより、発生する。そうして、例えば、２０４８個のＤＯＰＯの任意のペアを結合でき、合計４００万の組み合わせが得られる。

【0031】

本実施例の測定・フィードバック法を用い、４色定理に基づく地図の塗分け問題の実装法が説明される。この問題は、隣接する県同士が同じ色にならないように、各都道府県を４色で塗分ける問題である。この問題を解くためによく用いられている手法は、各都道府県に４色に対応するスピンを割り当て、４つのスピンのうちのいずれか一つだけが上向きスピンとなり、他は下向きスピンとなるような条件を入れると同時に、隣接する県同士のスピンの値は逆となるような条件を入れるものである。ここでは、赤（Ｒ）,緑（Ｇ）,青（Ｂ）,及び黄（Ｙ）の４色を用いるとする。このとき、各県のＲ,Ｇ, Ｂ, 及びＹに対応するスピンを全結合し、全体にバイアス項を印加することにより、４つのスピンのうちいずれか１つが上向きになった時が最もエネルギーが下がるようにすることができる。この原理を利用して上に述べた「４つのスピンのうちのいずれか一つだけが上向きスピンとなる」、すなわち、各県の色が一色に定まるための条件を入れることができる（図４(ａ)）。さらに、隣接県の同色のスピン結合を反強磁性結合にすることで、隣接県が同じ色にならないための条件を入れることができる(図４(ｂ))。

【0032】

イジングモデルの複数のスピンに対応し、同一の発振周波数を有する複数の光パルスを０またはπの位相でパラメトリック発振させる位相感応増幅器２と、位相感応増幅器２から出力された複数の光パルスの位相および振幅を測定する、光パルス測定部５と、光パルス測定部５において測定した光パルスの位相および振幅の情報を入力として、イジングモデルの結合係数と測定した光パルスとから決定される、ある光パルスに関わる相互作用に基づいて、フィードバック値を算出し、フィードバック信号７を出力する、高速演算回路６と、高速演算回路６において算出されたフィードバック信号７に基づいて、入力された複数の光パルスと同じ数の複数の光パルスの位相および振幅を変調することにより、ある光パルスに関わる相互作用を実装する光変調器８と、を備え、光変調器８は、ｐｕｓｈ－ｐｕｌｌ変調器であって、マッハツェンダ干渉計の両アーム導波路に配置された２つの位相変調器を絶対値が同じで符号が異なる信号で駆動し、光変調器８の動作点を透過する光強度が最小となる点で駆動するように構成され、フィードバック信号７を重畳する際に、動作点を、印加した磁場に相当する電圧分シフトさせることにより、光共振器のＤＯＰＯと同じ波長で、かつポンプ光との初期位相差が固定されているＤＯＰＯのパルスに定数の光が注入されるように構成されていることを特徴とする、イジングモデルの計算装置を用いることにより、測定・フィードバック法に、模擬した磁場項を入れることで、より多彩な組合せ最適化問題を実装することが可能となった。

【0033】

［実施例２］
もっとも単純な方法で、動作点をずらす方法が説明される。本実施例では、動作点の電圧のシフトとして、定数電圧を用いる。課題のなっている局所磁場hiに応じて(式(１))、定数で動作点をずらす方法である。一般的にはi番目のＤＯＰＯに対してそれぞれ違う値が割り当てられるので、決まったパターンでパルス光を変調することになる。i番目のＤＯＰＯへフィードバックするために、変調器に印加する信号は、スピン間の相互作用に相当する式（１ｂ）の右辺第２項に定数の磁場項を加えた

【0034】

【数4】

【0035】

となる。ここでｃｊはｊ番目のパルス振幅の余弦成分の測定結果、ｒは比例係数であり、 (２)式の右辺第１項は、測定結果を演算してフィードバックするために用いられるＦＰＧＡからの出力信号である。第２項も同様に、高速演算回路（ＦＰＧＡ）６から出力がされても良いし、別の波形発生器から出力された信号がＦＰＧＡ出力信号と電気的に加算されても良い。

【0036】

［実施例３］
本実施例では、動作点の電圧のシフトとして、あらかじめ定められた関数に従って電圧を変化させる方法である。式(２)において、局所磁場hiが時間(ＤＯＰＯの周回数)と共にあらかじめ定められた関数で変化するものである。例えば、線形増加が用いられればよい。

【0037】

【数5】

【0038】

［実施例４］
本実施例では、i番目（iは、自然数）のＤＯＰＯにおける上記動作点の電圧のシフトとして、i番目のパルスの測定の結果得られたＤＯＰＯ振幅ｃｉに比例する電圧を用いる。

【0039】

ＤＯＰＯのパルスは、計算の開始時には光パワー０の状態で始まり、位相感応増幅器２から出力される雑音光が種となって、光共振器１中を周回し、位相感応増幅器２により徐々に増幅されていくため、計算の初期状態にはパルス毎にその振幅の大きさはまちまちである。それを保証するのが実施例４～６である。もっとも単純には、測定で得られたciに比例した信号を磁場項として追加する。

【0040】

【数6】

【0041】

［実施例５］
本実施例では、i番目（ｎは、自然数）のＤＯＰＯにおける動作点の電圧のシフトとして、i番目のパルスの測定の結果得られたＤＯＰＯ振幅の絶対値c_iに比例する電圧を用いる。c_iの絶対値に比例する信号が磁場項として追加される。

【0042】

【数7】

【0043】

［実施例６］
本実施例では、全てのパルスの測定結果の絶対値の平均値に比例する信号を磁場項として追加する。具体的には、i番目（iは、自然数）のＤＯＰＯにおける上記動作点の電圧のシフトとして、全てのパルスの測定の結果得られたＤＯＰＯ振幅の絶対値の平均

【0044】

【数8】

【0045】

【数9】

【0046】

［実施例７］
本実施例では、式(２)～(６)及び(８)のカッコ内第２項の組合せの信号を磁場項として追加する。上述の実施例２乃至６に記載された構成うちの複数個を同時に組み合わせて式 (９)を満たすように、動作がされる。

【0047】

【数10】

【0048】

なお、実施例３～７も実施例２と同様、高速演算回路（ＦＰＧＡ）６からフィードバック信号７が出力されても良い。また、別の波形発生器から出力された信号が高速演算回路（ＦＰＧＡ）６の出力信号と電気的に加算されても良い。

【産業上の利用可能性】

【0049】

本開示は、組み合わせ最適化問題を解くイジングモデルの計算装置に関する技術分野に適用できる。

【図1】

【図2】

【図3】

【図4】

【図5】