知財求人 - 知財ポータルサイト「IP Force」
▶ マイクロソフト テクノロジー ライセンシング，エルエルシーの特許一覧
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公表特許公報(A)
(11)【公表番号】
(43)【公表日】2024-02-21
(54)【発明の名称】ハードウェアソルバアーキテクチャ
(51)【国際特許分類】
G06N 99/00 20190101AFI20240214BHJP
G06E 3/00 20060101ALI20240214BHJP
【ＦＩ】
G06N99/00 180
G06E3/00
【審査請求】未請求
【予備審査請求】未請求
(21)【出願番号】P 2023547509
(86)(22)【出願日】2022-01-28
(85)【翻訳文提出日】2023-08-29
(86)【国際出願番号】 US2022014172
(87)【国際公開番号】W WO2022169670
(87)【国際公開日】2022-08-11
(31)【優先権主張番号】21155447.2
(32)【優先日】2021-02-05
(33)【優先権主張国・地域又は機関】EP
(81)【指定国・地域】
(71)【出願人】
【識別番号】314015767
【氏名又は名称】マイクロソフト テクノロジー ライセンシング，エルエルシー
(74)【代理人】
【識別番号】100118902
【弁理士】
【氏名又は名称】山本 修
(74)【代理人】
【識別番号】100106208
【弁理士】
【氏名又は名称】宮前 徹
(74)【代理人】
【識別番号】100196508
【弁理士】
【氏名又は名称】松尾 淳一
(72)【発明者】
【氏名】パルミジャーニ，フランチェスカ
(72)【発明者】
【氏名】ガンツシディス，クリストス
(72)【発明者】
【氏名】ハラー，イシュトヴァン
(72)【発明者】
【氏名】バッラーニ，ヒテシュ
(72)【発明者】
【氏名】ローストロン，アントニー・イアン・テイラー
(72)【発明者】
【氏名】リバルチェンコ，アンドレイ
(57)【要約】
変数間の相互作用の重み付き和を含む関数を最適化する変数のベクトルの値を推定するシステムであって、各々が関数に対するそれぞれの変数の寄与をモデル化するように配置された複数の並列ハードウェアチャネルを含み、各チャネルは、それぞれの変数の値をモデル化する信号を生成するように構成された信号生成器と、信号をチャネルのそれぞれに供給するように配置されたスプリッタと、信号のベクトルに、それぞれの変数と変数のベクトルとの間の相互作用をモデル化する重みのベクトルを乗算し、それによってそれぞれの変数の寄与を表すフィードバック信号を生成するように構成された相互作用ロジックと、フィードバック信号に応じて信号を適応させるように構成された信号生成器にフィードバック信号を戻すためのフィードバック経路とを含み、各チャネルは、光学コンポーネントおよび／またはアナログ電子部品のみを使用して実装される、システム。
【特許請求の範囲】
【請求項1】
関数を最適化する変数のベクトルの値を推定するためのシステムであって、前記関数は、複数の項の重み付き和を含み、各項は、前記ベクトルからの前記変数の対応する部分集合の積を含み、各項は、前記変数間の相互作用をモデル化する重みの行列からの対応する重みによって重み付けされ、前記システムは、互いに同時に動作するように配置された複数の並列ハードウェアチャネルを含み、各チャネルは、前記関数に対する前記変数のそれぞれの１つの寄与をモデル化するように配置され、前記並列チャネルのそれぞれは、それぞれの前記変数の値をモデル化する変調された特性を有するそれぞれのモデリング信号を生成するように構成されたそれぞれの信号生成器と、
それぞれの前記モデリング信号のインスタンスを前記並列チャネルのそれぞれに供給するように配置されたそれぞれのスプリッタであって、各チャネルは、したがって、変数の前記ベクトルをモデル化した信号のベクトルを受信する、スプリッタと、
それぞれのベクトル乗算器を含むそれぞれの相互作用ロジックであって、それぞれの前記ベクトル乗算器は、受信した信号の前記ベクトルに、それぞれの前記変数と変数の前記ベクトルとの間の相互作用をモデル化する重みの前記行列からの重みのそれぞれのベクトルを乗算するように構成され、前記相互作用ロジックは、それによって、それぞれの前記チャネルによってモデル化されたそれぞれの前記変数の前記寄与を表すそれぞれのフィードバック信号を生成する、それぞれの相互作用ロジックと、
前記フィードバック信号をそれぞれの前記信号生成器に戻すように配置されたそれぞれのフィードバック経路であって、それぞれの前記信号生成器は、前記フィードバック信号に応じて、それぞれの前記モデリング信号を適応させるように構成される、フィードバック経路と
を含み、
各チャネルにおいてそれぞれの前記信号生成器、スプリッタ、相互作用ロジックおよびフィードバック経路を含む各チャネルは、光学コンポーネントおよびアナログ電子部品の一方または双方のみを使用して実装される、
システム。
【請求項2】
請求項１に記載のシステムであって、前記変数のそれぞれは、二値である、システム。
【請求項3】
請求項１または２に記載のシステムであって、前記変数のそれぞれは、正または負の値を取ることができ、前記関数は【数1】
σ１．．．Ｎは、前記ベクトルの前記変数であり、Ｊｉｊは重みの前記行列であり、ｉとｊは１からＮの間の前記変数のインスタンスを列挙するためのインデックスであり、Σｉｊは前記関数に含まれる変数のすべての部分集合（σｉ，σｊ）にわたる和を表し、
前記システムは、各ｉに対する前記チャネルのそれぞれの１つを含み、それぞれの前記ベクトル乗算器は、乗算Σｊσｊ・Ｊｉｊを実行するように構成され、したがって、前記関数に対するσｉのそれぞれの前記寄与をモデル化し、
前記信号生成器は、前記ハミルトニアンのエネルギーを最小化するように前記モデリング信号の前記適応を実行することによって前記関数を最適化するように構成される、
システム。
【請求項4】
請求項３に記載のシステムであって、前記正および負の値は、＋１／２および－１／２、または、＋１および－１である、システム。
【請求項5】
請求項２に記載のシステムであって、前記変数のそれぞれはブールであり、前記関数は、【数2】
ν１．．．Ｎは、前記ベクトルの前記変数であり、Ｑｉｊは重みの前記行列であり、ｉとｊは１からＮの間の前記変数のインスタンスを列挙するためのインデックスであり、Σｉｊは前記関数に含まれる変数のすべての部分集合（νｉ，νｊ）にわたる和を表し、
前記システムは、各ｉに対する前記チャネルのそれぞれの１つを含み、それぞれの前記ベクトル行列乗算器は、乗算Σｊνｊ・Ｑｉｊを実行するように構成され、したがって、前記関数に対するνｉのそれぞれの前記寄与をモデル化し、
前記信号生成器は、前記関数を最小化するように、前記モデリング信号の前記適応を実行することによって、前記関数を最適化するように構成される、
システム。
【請求項6】
請求項１から５のうちのいずれか一項に記載のシステムであって、各チャネルにおいて、それぞれの前記フィードバック経路は、それぞれの前記信号生成器に戻る前に、それぞれの前記フィードバック信号にそれぞれのノイズ成分を導入するように配置される、システム。
【請求項7】
請求項１から６のうちのいずれか一項に記載のシステムであって、それぞれの前記変数の前記値をモデル化するために使用される各モデリング信号の変調された前記特性は、前記信号の振幅と、
前記信号の位相と
のうちの１つを含む、システム。
【請求項8】
請求項１から７のうちのいずれか一項に記載のシステムであって、各チャネルにおいて、それぞれの前記信号生成器によって生成されたそれぞれの前記モデリング信号は、光信号を含み、それぞれの前記ベクトル乗算器は、光領域においてそのそれぞれの乗算を実行するように構成された光乗算器を含む、システム。
【請求項9】
請求項８に記載のシステムであって、それぞれの信号生成器によって生成される前記光信号は、異なるそれぞれの光波長を有する、システム。
【請求項10】
請求項８または９に記載のシステムであって、各ベクトル乗算器の前記光乗算器は、空間光変調器と、
波長選択スイッチと、
リング共振器と、
マッハツェンダー干渉計と
のうちの１つを含む、
システム。
【請求項11】
請求項８、９または１０に記載のシステムであって、各相互作用ロジックは、アナログ電子形式でそれぞれの前記信号生成器への前記フィードバック信号を発生するために、それぞれの前記光乗算器の光出力を検出するように配置されたそれぞれの光検出器をさらに含む、システム。
【請求項12】
請求項１１に記載のシステムであって、それぞれの信号生成器は、それぞれの光源と、それぞれの前記フィードバック信号に応じて、アナログ電子信号の形式でそれぞれのスピンを生成するように構成されたそれぞれのスピン生成器と、それぞれの前記スピンをそれぞれの前記光信号の変調された前記特性に変調するように配置されたそれぞれの変調器とを含む、システム。
【請求項13】
請求項１２に記載のシステムであって、各チャネルにおけるそれぞれの前記スピン生成器は、さらなる光源と、それぞれの前記フィードバック信号に応じて前記さらなる光源からの光を変調するように配置されたさらなる変調器と、前記さらなる変調器からの変調された前記光を検出し、それに応じて前記スピンを生成するように配置されたさらなる光検出器とを含む、システム。
【請求項14】
請求項１２または１３に記載のシステムであって、各チャネルにおいて、それぞれの前記変数の前記値をモデル化するために使用される各信号の前記特性は、前記信号の振幅を含み、
それぞれのスピン生成器は、前記スピンを表すために正レベルと負レベルとの間で振幅変調されたアナログ電子信号の形式でそれぞれの前記スピンを生成するように配置され、
それぞれの前記変調器は、前記スピンをそれぞれの前記光信号の前記振幅に変調するように構成され、前記光信号の前記振幅およびそれぞれの前記光乗算器の前記出力は、正のレベルのみを取ることができ、負のレベルは取ることができず、
各チャネルにおける前記相互作用ロジックのそれぞれの前記光検出器は、インコヒーレント検出によってそれぞれの前記光乗算器の前記出力を検出し、それによってそれぞれの読み取り値を生成するように構成され、
各チャネルにおけるそれぞれの前記相互作用ロジックは、正レベルと負レベルの間で変調されたそれぞれの前記フィードバック信号を発生するために、それぞれの前記光検出器からのそれぞれの前記読み取り値にＤＣオフセットを加えるように構成される、
システム。
【請求項15】
関数を最適化する変数のベクトルの値を推定するための方法であって、前記関数は、複数の項の重み付き和を含み、各項は、前記ベクトルからの前記変数の対応する部分集合の積を含み、各項は、前記変数間の相互作用をモデル化する重みの行列からの対応する重みによって重み付けされ、前記方法は、複数の並列ハードウェアチャネルのそれぞれにおいて、それぞれの信号生成器によって、それぞれの変数の値をモデル化する変調された特性を有するそれぞれのモデリング信号を生成するステップと、
それぞれの前記モデリング信号のインスタンスを前記並列チャネルのそれぞれに供給するステップであって、各チャネルは、したがって、変数の前記ベクトルをモデル化した信号のベクトルを受信する、ステップと、
それぞれの相互作用ロジックにおいて、受信した信号の前記ベクトルに、それぞれの前記変数と変数の前記ベクトルとの間の相互作用をモデル化する重みの前記行列からの重みのそれぞれのベクトルを乗算し、それによって、それぞれの前記チャネルによってモデル化されたそれぞれの前記変数の前記寄与を表すそれぞれのフィードバック信号を生成するステップと、
前記フィードバック信号をそれぞれの前記信号生成器に戻すステップと、
それぞれの前記信号生成器によって、前記フィードバック信号に応じて、それぞれの前記モデリング信号を適応させるステップと
を含み、
前記方法は、光学ハードウェアおよびアナログ電子ハードウェアの一方または双方のみを使用して実施される、
方法。
【発明の詳細な説明】
【技術分野】
【0001】
[001]ロジスティクス、財務ポートフォリオ管理、創薬、およびその他の応用領域における多くの問題は、目的を最適化することを目標として、入力（通常、変数と呼ばれる）への値の割り当てを見つけることを必要とする。
【背景技術】
【0002】
例えば、このような問題には「組合せ最適化問題」が含まれる。他の最適化分野とは異なり、組合せ最適化は、変数が有限集合から値を取る問題に関連する。例えば、有効な代入は、限定された集合（例えば、３つの利用可能なルートから１つを選ぶ）、または一般的には整数の有限部分集合からの二値（例えば、投資をするかしないか）であり得る。このような問題では、各変数の値の組合せ方は有限である。原理的には、すべての可能な組合せを列挙し、最適な割り当てを見つけることが可能である。しかし、実際には、組合せの集合が非常に大きい（変数の数が指数関数的）ため、このような網羅的な探索は、中程度の大きさの問題では実行不可能である。
【0003】
[002]このような問題の構造を理解するための広範な研究が行われてきた。組合せ最適化問題の部分集合は、ＮＰ完全（ＮＰは非決定性多項式を意味する）として知られる問題のカテゴリに属する。ＮＰ完全性は計算複雑さ理論で知られている概念であり、すべてのＮＰ完全問題は他のＮＰ完全問題に変換することができる。任意のＮＰ完全問題に対する効率的なソルバは、任意のＮＰ完全問題を効率的に解くことができることを意味する。すべてのＮＰ完全問題はまた、「ＮＰ困難」として知られるより大きな問題の部分集合に属し、すべてのＮＰ困難問題はまた、他のすべてのＮＰ困難問題に変換することができる。
【0004】
[003]この設定における「効率的」という用語は、すべての可能性を列挙することなく問題の解を見つけることを意味する。具体的には、本明細書で説明するグラフ最適化に対する効率的な解とは、解を求めるのに要する時間が問題の変数（グラフの頂点など）の数に対して多項式になるような解であり、一方、すべての可能な解を列挙することは問題の変数の数に対して指数関数的になるような解である。しかし、そのような効率的なソルバは存在しないと広く受け入れられている。その代わりに、この分野での研究は、「十分に良い」解を見つけるアルゴリズムを考案することに集中している。多くの場合、そのような近似アルゴリズムが、厳密解に十分に近い解を本当に提供するという保証はない。
【0005】
[004]組合せ最適化問題には様々なものが存在し、上述したように、ＮＰ完全問題は他のＮＰ完全問題に変換することができる。例えば、巡回セールスマン問題は次のように定義される：与えられた都市集合と都市間の対距離に対して、すべての都市を経由する経路を求める問題であり、各都市はちょうど１回訪問され、経路の全長は最短である。
【0006】
[005]組合せ最適化問題の一般的な形式は、二次無制約二値最適化（ＱＵＢＯ）問題と呼ばれ、これは、二値変数Ｖ＝｛ν１，ν２，．．．νＮ｝の集合によって定義され、それぞれが０または１のいずれかの値、および式ΣｉΣｊＱｉｊ・νｉ・νｊを取り、ここで係数Ｑｉｊは変数νｉとνｊの間の相互作用を定義する。巡回セールスマン問題は、変数を訪問可能な各都市の各都市間の経路における位置として定義し（例えば、最初の変数はロンドンが最初に訪問する都市であるか否かを示し得る）、すべての都市をちょうど１回訪問するという制約の下で総距離が最小になるように都市間の距離を行列Ｑに符号化することで、ＱＵＢＯ問題として定式化することができる。ＱＵＢＯは多項式無制約二値最適化（ＰＵＢＯ）問題の一種であり、式Σｖ’⊆ｖＱｖ’・Πν∈Ｖ’νを最小化するように二値変数Ｖ＝｛ν１，ν２，．．．νＮ｝の集合に値を割り当て、この場合、Ｑは、任意の数の変数間の相互作用を符号化し得る。上述したように、ＮＰ困難問題の異なる定式化間の変換は可能である。式に補助変数および補助項を導入することで、ＰＵＢＯ問題をＱＵＢＯに変換することができる。
【0007】
[006]物理学で強磁性および他の物理過程をモデル化するために使用されるイジングモデルと呼ばれる定式化は、上で定義したＱＵＢＯ問題と等価である。イジングモデルは、２つの離散的な状態に存在し得る変数を持つ物理システムで記述され、これらの変数は互いに相互作用することができ、システムの全エネルギーはＨ（σ）＝－Σｉ，ｊＪｉｊσｉσｊ－μΣｉｈｉσｉで与えられる。イジング定式化では、二値変数（「スピン」と呼ばれることもある）は通常、１／０または他の二値割り当てではなく、＋１／－１のいずれかに割り当てられる。ただし、イジング定式化は、σｉ＝２νｉ－１の式を適用することで、ブール変数のＱＵＢＯ定式化に簡単にマッピングすることができる。
【0008】
[007]上記で使用された表記法は、ＱＵＢＯとイジング定式化の間でわずかに異なっており、割り当てられる変数はσで表され、相互作用係数はイジングモデルのＪで表されることに注意されたい。簡単にするため、このアプリケーションでは、値が割り当てられる変数にはσという表記を使用し、イジングソルバの文脈では、相互作用係数の配列を表すためにＪを使用する。しかし、本明細書では、重みの行列および変数をそれぞれ示すために、Ｑおよびνを一般的に使用することもできる。
【0009】
[008]総エネルギーに関する上記の式の第２項－μΣｉｈｉσｉは、モデル化されているシステムに対する外部「場」または何らかの外部効果の影響を表す。例えば、強磁性体の場合、第１項は磁気双極子間の相互作用によるエネルギー寄与を表し、第２項は外部磁場によるシステムのエネルギーを表す。外部磁場がない場合のイジング問題を解く方がはるかに単純であるため、多くの問題は外部磁場なしのイジング問題としてモデル化される。しかし、余分なスピンを導入し、重みを慎重に選択したエッジを追加することで、外部磁場のある問題を外部磁場のない問題に変換することは可能である。以下の説明でイジング問題と呼ばれる問題またはモデルは、外部磁場のない問題、またはすでに外部磁場のない問題に変換された問題を想定する。
【0010】
[009]最近まで、組合せ最適化のアルゴリズムは、一般に、汎用ＣＰＵ、ＦＰＧＡ、ＧＰＵ、ＡＳＩＣなどのデジタルハードウェアで実装されてきた。デジタルハードウェアは、柔軟性（すなわち、異なるアルゴリズムをプログラムする能力）および信頼性に関して大きな利点を有する。しかし、デジタルソリューションには、実行速度と消費電力による制限もある。かつては、デジタルハードウェアの世代が進むごとに、計算能力の向上と消費電力の削減が達成されてきた。しかし、基本的な物理的限界に近づくにつれて、デジタルハードウェアの性能向上はますます困難になると広く予測されている。組合せ最適化問題に対するより良い解答を探したり、より大きなインスタンスに取り組んだりすることは、より大きなハードウェアコストを伴うことになる。
【0011】
[0010]しかし、最近では、非デジタル物理プロセスに基づくハードウェアを用いて、このような問題を解決しようとする試みがなされている。物理プロセスによるイジングモデルの一般的な実現には、量子アニーラが使用されている。既存のシステムでは、問題の変数は量子ビットで表現され、通常「スピン」と呼ばれる＋１と－１の値を取る。しかし、このトポロジーでは完全な接続性が得られない。その代わりに、量子ビットは、４つの水平量子ビットがカプラを介して４つの垂直量子ビットに接続されたユニットセルの集合を含むアーキテクチャで相互接続される。ユニットセルは垂直方向と水平方向にタイル状に配置され、隣接する量子ビットが接続され、まばらに接続された量子ビットの格子が形成される。このアーキテクチャの接続性が制限されているため、望ましくない影響があり、その結果、スピンへの問題変数の非効率な表現が生じる、つまり、問題を表現するために物理システムに必要な量子ビットの数は、元の変数の数よりもはるかに多くなる。
【0012】
[0011]量子アニーラハードウェアのこの固有の物理的制限のため、古典的ハードウェアで実行でき、量子の物理的特性に触発されたアルゴリズムが開発されてきた。例えば、Ｍｉｃｒｏｓｏｆｔ Ａｚｕｒｅは、量子インスパイアされた最適化（ＱＩＯ）アルゴリズムを開発し、ＰＵＢＯ問題の近似に有望であることを示している。
【0013】
[0012]光ソルバでは、入力変数（例えば、イジング問題のσｉ＝１．．．Ｎ）を表すために光信号が使用され、変数間の相互作用（例えば、イジング問題の行列Ｊ）をモデル化する方法で信号を結合するために光学素子が使用される。液晶ディスプレイまたはリング共振器など、光学領域でベクトル同士の乗算を行う光学素子は当技術分野で知られている。加算（Σ）は、光検出器の光検出素子に入射する信号のコヒーレント加算またはインコヒーレント加算を行うことができる光検出器を用いて実施することができる。
【0014】
[0013]ソルバへの入力（すなわち、値が決定される変数）が、－１と＋１、または－１／２と＋１／２のような正負の二値を取り得る場合、これらは、単にスピンの量子特性との類推によって「スピン」と呼ばれることがある。しかし、このような文脈では、実際にはスピンの量子特性を意味しない。代わりに、２つの可能な「スピン」値は、単に二値変数の２つの可能な値を指し、例えば、光の振幅または位相の２つの異なる値を使用して表される可能性がある。
【0015】
[0014]光学に基づく、あるいは光学に着想を得た最先端のソリューションは、デジタルアプローチのみ（東芝のＳＢＭ［模擬分岐マシン］による巡回セールスマン問題の解決、Ｉｋｕｋｏ Ｈａｓｕｍｉ ｈｔｔｐｓ：／／ｍｅｄｉｕｍ．ｃｏｍ／ｔｏｓｈｉｂａ－ｓｂｍ／ｓｏｌｖｉｎｇ－ｔｒａｖｅｌｉｎｇ－ｓａｌｅｓｍａｎ－ｐｒｏｂｌｅｍ－ｗｉｔｈ－ｓｂｍ－ｓｉｍｕｌａｔｅｄ－ｂｉｆｕｒｃａｔｉｏｎ－ｍａｃｈｉｎｅ－８９７４０ｃ８３ｅｄ３７を参照）、あるいはＢｏｈｍ，Ｆａｂｉａｎ，Ｇｕｙ Ｖｅｒｓｃｈａｆｆｅｌｔ，Ｇｕｙ Ｖａｎ ｄｅｒ Ｓａｎｄｅの「Ａ ｐｏｏｒ ｍａｎ’ｓ ｃｏｈｅｒｅｎｔ Ｉｓｉｎｇ ｍａｃｈｉｎｅ ｂａｓｅｄ ｏｎ ｏｐｔｏ－ｅｌｅｃｔｒｏｎｉｃ ｆｅｅｄｂａｃｋ ｓｙｓｔｅｍｓ ｆｏｒ ｓｏｌｖｉｎｇ ｏｐｔｉｍｉｚａｔｉｏｎ Ｐｒｏｂｌｅｍｓ（最適化問題を解決するための光電子フィードバックシステムに基づく貧乏人のコヒーレントイジングマシン）」Ｎａｔｕｒｅ ｃｏｍｍｕｎｉｃａｔｉｏｎｓ １０．１（２０１９）：１－９（ｈｔｔｐｓ：／／ｗｗｗ．ｎａｔｕｒｅ．ｃｏｍ／ａｒｔｉｃｌｅｓ／ｓ４１４６７－０１９－１１４８４－３）およびＩｎａｇａｋｉ，Ｔａｋａｈｉｒｏらの「Ａ ｃｏｈｅｒｅｎｔ Ｉｓｉｎｇ ｍａｃｈｉｎｅ ｆｏｒ ２０００－ｎｏｄｅ ｏｐｔｉｍｉｚａｔｉｏｎ Ｐｒｏｂｌｅｍｓ（２０００ノードの最適化問題に対するコヒーレントイジングマシン）」Ｓｃｉｅｎｃｅ ３５４．６３１２（２０１６）：６０３－６０６（ｈｔｔｐｓ：／／ｓｃｉｅｎｃｅ．ｓｃｉｅｎｃｅｍａｇ．ｏｒｇ／ｃｏｎｔｅｎｔ／３５４／６３１２／６０３）のようなハイブリッドアプローチのいずれかを提案している。
【0016】
[0015]ハイブリッドアプローチでは、変数値を表す信号を生成するためのビルディングブロックは、通常、光ハードウェアで実装されるが、変数相互作用を計算するロジックは、デジタルハードウェアと、光領域とデジタル領域とを変換するハードウェアとを用いて実装される。これとは対照的に、「全アナログ」ソルバでは、代わりに非デジタルハードウェアが、光（すなわち光信号）とアナログ電子領域の間で信号を変換するために使用される。全アナログソルバの利点は、光信号とアナログ電子信号の伝送速度が速いことである（デジタルエレクトロニクスは、フリップフロップを介してビット列をクロックする必要があるため、本質的にはるかに遅い）。一方、反復の一部をデジタル領域で実装することは、デジタルエレクトロニクスと比較して伝送速度が速いという全アナログソルバの利点を失うことになる。システムの速度は最も遅い部分によって制限されるため、デジタルエレクトロニクスを含めると、光ソルバの利点が否定されることになる。
【0017】
[0016]４つのスピン／変数に対する全対全接続性と１６のスピン／変数に対する部分的接続性を持つ全光ソリューションが提案され、実証されている。（Ｍａｒａｎｄｉ，Ａ．，Ｗａｎｇ，Ｚ．，Ｔａｋａｔａ，Ｋ．，Ｂｙｅｒ，Ｒ．Ｌ．＆ Ｙａｍａｍｏｔｏ，Ｙ．による「Ｎｅｔｗｏｒｋ ｏｆ ｔｉｍｅ－ｍｕｌｔｉｐｌｅｘｅｄ ｏｐｔｉｃａｌ ｐａｒａｍｅｔｒｉｃ ｏｓｃｉｌｌａｔｏｒｓ ａｓ ａ ｃｏｈｅｒｅｎｔ Ｉｓｉｎｇ ｍａｃｈｉｎｅ（コヒーレントイジングマシンとしての時間多重光パラメトリック発振器のネットワーク）」Ｎａｔ．Ｐｈｏｔｏｎｉｃｓ ８，９３７－９４２（２０１４），Ｋ．Ｔａｋａｔａらによる「Ａ １６－ｂｉｔ Ｃｏｈｅｒｅｎｔ Ｉｓｉｎｇ Ｍａｃｈｉｎｅ ｆｏｒ Ｏｎｅ－Ｄｉｍｅｎｓｉｏｎａｌ Ｒｉｎｇ ａｎｄ Ｃｕｂｉｃ Ｇｒａｐｈ Ｐｒｏｂｌｅｍｓ（１次元リングおよび三次グラフ問題用の１６ビットコヒーレントイジングマシン）」，Ｓｃｉｅｎｔｉｆｉｃ Ｒｅｐｏｒｔ ２０１６．「Ａ １６－ｂｉｔ Ｃｏｈｅｒｅｎｔ Ｉｓｉｎｇ Ｍａｃｈｉｎｅ ｆｏｒ Ｏｎｅ－Ｄｉｍｅｎｓｉｏｎａｌ Ｒｉｎｇ ａｎｄ Ｃｕｂｉｃ Ｇｒａｐｈ Ｐｒｏｂｌｅｍｓ（１次元リングおよび三次グラフ問題用の１６ビットコヒーレントイジングマシン）」（ｅｕｒｏｐｅｐｍｃ．ｏｒｇ））。
【0018】
[0017]最先端の全光ソルバは、時分割多重アーキテクチャの光信号を用いて変数を生成する。すなわち、信号は同じ光ビームに直列に多重化され、変数間の相互作用をモデル化するために信号を結合できるように、変数ごとに異なる遅延経路が導入される。しかし、時分割多重の場合、スピン生成が直列に行われるため、ソルバの時間的複雑さはモデル化される変数の数に比例する。図５Ａは、時分割多重化アーキテクチャの概略ブロック図である。時分割アーキテクチャの場合、各「スピン」は直列に生成され、それらの間の相互作用を実装できるように、異なるパス長だけ遅延される。これは、スピン間の相互作用を計算し、スピン生成にフィードバックを渡す前に、１回の反復ですべてのスピンが生成されるのを待たなければならないことを意味する。
【0019】
[0018]完全にアナログ領域で実装されるソルバは、イジングシステムの「スピン」相互作用をモデル化するために光学的またはアナログ電子ベクトル乗算器を使用する。光学的なベクトル乗算器を実装することは、光通信の速度を利用するという、上で指摘した利点がある。
【発明の概要】
【0020】
[0019]本開示は、（光領域、アナログ電子領域、またはその組合せのいずれかの）全アナログシステムにおいて実装されるソルバアーキテクチャに関し、これは、イジング問題のようなベクトル行列乗算を含む問題の近似解を決定するように構成される。
【0021】
[0020]上述したように、既存の全光ソルバは、問題の異なる変数を表す信号（すなわち、「スピン信号」）を直列に生成するために時分割多重化を使用し、次に、それらを結合して変数間の相互作用をモデル化するために異なる経路長の遅延線を使用する（例えば、図５Ａを参照）。これは、時間の複雑さがモデル化される変数の数に対して線形であることを意味する。
【0022】
[0021]一方、本開示は、空間分割アーキテクチャを採用する、光および／またはアナログ電子領域で実装された全アナログソルバを開示する（例えば、後でさらに詳細に説明する図５Ｂを参照）。言い換えれば、スピン生成は並列に行うことができ、反復あたりの時間はスピン（またはより一般的には変数）の数で準一定であることを意味する。
【0023】
[0022]アーキテクチャは、モデル化される各変数が異なるチャネルに関連付けられるように構成され、光学またはアナログ電子ハードウェアは、最適化される全体的な関数の更新に対する各変数の寄与を計算するように構成される。光ソルバまたはアナログ電子ソルバを使用する利点は、計算が光または電気信号の速度を利用するため、計算をより迅速に実行できることであり、これらの信号を検出するのにかかる時間によってのみ制限されるが、デジタルハードウェアの電力および速度の考慮事項の一部（ラッチを通してビットをシフトする必要性など）によって制限されることはない。
【0024】
[0023]さらに、変数を複数の並列チャネルに分割すること（すなわち、既存のソルバアーキテクチャで使用される時分割多重化ではなく、空間分割多重化）は、問題の時間的複雑性を大幅に増加させることなく、より多くの変数の数に拡張できるため有利である。既存の時分割ベースの光ソルバは、変数間の相互作用を計算するために各変数に遅延経路を導入する必要があるため、モデル化できる変数の数に制限があり、時分割全光ソルバの反復ごとにかかる時間は、問題の変数の数に対して線形に拡張する。さらに、既存のハイブリッドソルバは、すべてのデジタルソルバと同じ制約によって、最終的に速度が制限される。
【0025】
[0024]本開示の第１の態様は、関数を最適化する変数のベクトルの値を推定するシステムを提供する。関数は、複数の項の重み付き和を含み、各項は、ベクトルからの変数の対応する部分集合の積を含み、各項は、変数間の相互作用をモデル化する重みの行列からの対応する重みによって重み付けされる。システムは、互いに同時に動作するように配置された複数の並列ハードウェアチャネルを含み、各チャネルは、関数に対する変数のそれぞれの１つの寄与をモデル化するように配置される。並列チャネルのそれぞれは、それぞれの変数の値をモデル化する変調された特性を有するそれぞれのモデリング信号を生成するように構成されたそれぞれの信号生成器と、それぞれのモデリング信号のインスタンスを並列チャネルのそれぞれに供給するように配置されたそれぞれのスプリッタであって、各チャネルは、変数のベクトルをモデル化した信号のベクトルを受信する、スプリッタと、受信した信号のベクトルに、それぞれの変数と変数のベクトルとの間の相互作用をモデル化する重みの行列からの重みのそれぞれのベクトルを乗算するように構成されたそれぞれのベクトル乗算器を含み、それによって、それぞれのチャネルによってモデル化されたそれぞれの変数の寄与を表すそれぞれのフィードバック信号を生成するそれぞれの相互作用ロジックと、フィードバック信号をそれぞれの信号生成器に戻すように配置されたそれぞれのフィードバック経路であって、それぞれの信号生成器は、フィードバック信号に応じて、それぞれのモデリング信号を適応させるように構成される、フィードバック経路と、を含み、各チャネルにおいてそれぞれの信号生成器、スプリッタ、相互作用ロジックおよびフィードバック経路を含む各チャネルは、光学および／またはアナログ電子部品のみを使用して実装される。
【0026】
[0025]本開示の別の態様は、関数を最適化する変数のベクトルの値を推定する方法を提供する。関数は、複数の項の重み付き和を含み、各項は、ベクトルからの変数の対応する部分集合の積を含み、各項は、変数間の相互作用をモデル化する重みの行列からの対応する重みによって重み付けされる。方法は、複数の並列ハードウェアチャネルのそれぞれにおいて、それぞれの信号生成器によって、それぞれの変数の値をモデル化する変調された特性を有するそれぞれのモデリング信号を生成するステップと、それぞれのモデリング信号のインスタンスを並列チャネルのそれぞれに供給するステップであって、各チャネルは、変数のベクトルをモデル化した信号のベクトルを受信する、ステップと、それぞれの相互作用ロジックにおいて、受信した信号のベクトルに、それぞれの変数と変数のベクトルとの間の相互作用をモデル化する重みの行列からの重みのそれぞれのベクトルを乗算し、それによって、それぞれのチャネルによってモデル化されたそれぞれの変数の寄与を表すそれぞれのフィードバック信号を生成するステップと、フィードバック信号をそれぞれの信号生成器に戻すステップと、それぞれの信号生成器によって、フィードバック信号に応じて、それぞれのモデリング信号を適応させるステップと、を含み、方法は、光学および／またはアナログ電子ハードウェアのみを使用して実施される。
【0027】
[0026]本開示をより良く理解するために、また、本開示の実施形態がどのように実施され得るかを示すために、添付の図を参照する。
【図面の簡単な説明】
【0028】
【図1】[0027]２次制約なし二値最適化問題のための例示的な光ソルバアーキテクチャの概略ブロック図である。
【図2】[0028]光ソルバアーキテクチャの１チャネルのハードウェア実装の概略ブロック図である。
【図3】[0029]スピン生成ハードウェアの概略ブロック図である。
【図4】[0030]ソルバの動作中のアナログ領域と光領域との間の信号変換の概略ブロック図である。
【図7】[0034]光ベクトル乗算の概念を示す概略図である。
【図8】[0035]光ベクトル間乗算を実行するための波長選択スイッチの動作を示す図である。
【図9】[0036]光ベクトルごと行列乗算を実行するための修正波長選択スイッチアーキテクチャの動作を示す図である。
【発明を実施するための形態】
【0029】
[0037]図１は、組合せ最適化問題を解くように構成された光ソルバアーキテクチャの概略ブロック図である。本明細書において最適化問題を「解く」とは、近似解を求める可能性を含むことに留意されたい。以下の説明では、上述のイジング問題のような２変数の相互作用をモデル化する問題に焦点を当てる。しかし、３つ以上の変数間の相互作用をモデル化する高次のＰＵＢＯ問題は、変数の数が多いとはいえ、ＱＵＢＯ問題に変換することができ、同じ光学ハードウェアを使用して解くことができることに留意されたい。上述したように、このような問題の目標は、問題に特有の特定の関数を最小化する変数の割り当てを見つけることである。最適化において最小化される関数は、本明細書では、強磁性に対するイジングモデルのような特定の物理システムにおいてこの形式を取る全エネルギーを反映するために、「エネルギー」と呼ばれることがある。
【0030】
[0038]一般的な組合せ最適化問題は、まず、当該問題をＱＵＢＯ問題にマッピングし、これをイジング問題にマッピングすることによって解くことができる。多くの問題については、ＱＵＢＯ定式化へのマッピングが知られている。それ以外の問題については、マッピングを導出する必要があり得る。一般的なＮＰ困難問題のＱＵＢＯまたはイジング定式化へのマッピングは、それ自体、数学の専門家であれば理解できる話題である。例えば、ＰＵＢＯ形式で表現された問題、例えば、式ΣｉｊｋＱｉｊｋνｉνｊνｋを有する３次制約なし二値最適化問題は、余分な変数および項を導入することにより、ＱＵＢＯ問題として表現することができ、したがって、本明細書に開示されるイジングソルバによって解くことができる。本明細書で開示されるソルバは、イジング問題の解を提供し、その問題へのマッピングが見出され得る任意のＮＰ困難問題を解くのに使用することができる。
【0031】
[0039]この場合、問題は、全エネルギーがイジングハミルトニアン、すなわち－Σｉ，ｊＪｉｊσｉσｊ（外部磁場なしと仮定）によって与えられる物理システムにマッピングされる。与えられた問題をイジング系にマッピングするには、全エネルギー－Σｉ，ｊＪｉｊ・σｉ・σｊを最小化する（すなわち、Σｉ，ｊＪｉｊ・σｉ・σｊを最大化する）ことが問題を最適化することと等価になるように行列Ｊを決定する必要がある。
【0032】
[0040]したがって、ハミルトニアンは複数の項Ｊｉｊ・σｉ・σｊの和であり、各項は変数σｉ，σｊのそれぞれの部分集合と対応する重みＪｉｊとの積である。（したがって、第１項は１つの変数σ１にそれ自身と重みＪ１１を乗じた部分集合であり、第２項は部分集合σ１、σ２と重みＪ１２を乗じたものである）。
【0033】
[0041]見てわかるように、この和は、σｉを和の左側に取り出すことによって、一連のベクトルごとの（ドット積）乗算に分解することができる：
【0034】
【数1】
【0035】
[0043]応用例として、巡回セールスマン問題がある。簡単な例では、営業担当者が訪問する必要のある３つの都市があるとする：ロンドン、エディンバラ、カーディフである。これらはＱＵＢＯ問題の９つの変数でモデル化できる：ν１は最初にロンドンを訪問する、ν２は２番目にロンドンを訪問する、ν３は３番目にロンドンを訪問する、ν４は最初にエディンバラを訪問する、ν５は２番目にエディンバラを訪問する、ν６は３番目にエディンバラを訪問する、ν７は最初にカーディフを訪問する、ν８は３番目にカーディフを訪問する、ν９は３番目にカーディフを訪問することを表す。行列Ｑｉｊの要素は、対応する都市のペア間を移動する際のペナルティを表す。つまり、Ｑ１５（最初にロンドンを訪問し、２番目にエディンバラを訪問する）は、ロンドンからエディンバラまでの距離のペナルティを表している。Ｑ１９（最初にロンドンを訪問し、３番目にカーディフを訪問する）のような重みは、連続する都市間の距離だけで総距離が決まるため、この問題では意味がないため、０に設定されていることに注意されたい。Ｑ１２またはＱ１３（最初にロンドンを訪問し、２番目にロンドンを訪問する、または、最初にロンドンを訪問し、３番目にロンドンを訪問する）のような他の重みは、各都市を１回ずつ訪問するという制約を課すために、大きなペナルティ値に設定することができる。ＱＵＢＯ問題（ΣｉΣｊＱｉｊ・νｉ・νｊを最小化する）は、その後、イジング問題（ハミルトニアン項－Σｉ，ｊＪｉｊσｉσｊのエネルギーを最小化する）に変換され、図１のソルバシステムを使用して解くことができる。
【0036】
[0044]他の例として、２つの分子間の分子類似性を推定する分子類似性問題がある。これは例えば、ある分子が別の分子を阻害する可能性が高いことを推定し、薬剤に使用するために用いることができる。分子類似性をＱＵＢＯ問題としてモデル化すること自体は、当技術分野で知られている。
【0037】
[0045]問題をモデル化するイジングシステムのハミルトニアンを最小化する方向に、システムによって生成される信号を適応させるための更新規則が導出され得る。イジングモデルの可能な更新方程式は、以下のように記述される：
【0038】
【数2】
【0039】
【数3】
[0047]ここで、上記の括弧内は、２ΣｊＪｉｊ＜ｘｊ＞として評価することができる。更新の項は、各反復における各スピンの更新の大きさを制御して、システム全体が最小になるように適合させるために、定数αとβを乗算することができる。上式のζ項はガウシアンノイズで、非最適解で「スタック」するのを避けるために、各反復でシステムを摂動させるために適用される。最後に、式１のｃｏｓ２（）項は、
【0040】
【数4】
【0041】
[0048]図１に示すように、より詳細に後述するが、それぞれの変数をモデル化するそれぞれのモデリング信号ｘｉは、それぞれのチャネル１０２内のそれぞれの信号生成器１００（例えば光信号生成器）によって生成される。行列Ｊに対応する変数間の相互作用は、相互作用ロジック１０４によってモデル化される。
【0042】
[0049]上記の式は、後でさらに詳細に説明する。他の定式化も可能である。どのような定式化が使用されるにせよ、更新方程式の基本的な特性は、モデル化される物理システムが最小エネルギー（すなわち、上記で与えられたハミルトニアンの最小値）に向かうように、信号をプッシュまたは適応させることである。これは、システム全体のエネルギーに対する与えられた信号の寄与を表し、その符号が更新の方向を決定する項βΣｊＪｉｊｘｊ［ｋ］によって駆動される。言い換えれば、この項は、それぞれのチャネル１０２のそれぞれのモデリング信号ｘｉを適応させるために、信号生成器１０２にフィードバックを提供する。このフィードバックの符号は、それぞれのモデリング信号（ｘ）の適応を引き起こし、その信号を、システムのハミルトニアンにおける全体的なエネルギーを減少させる方向に駆動する。このフィードバックの値は、適応の程度を決定する（任意に、係数αとβによって信号ｘに対して減衰される）。
【0043】
[0050]ソルバが直接イジング変数σｉを表す信号を決定する一方で、これはＱＵＢＯ変数νｉの最適マッピングを見つけることと等価であり、元の問題の形で別の変数集合に変換することができることに留意されたい。ただし、ソルバが決定できるイジング変数（スピン）の集合と、元の問題を最適化する変数の集合との間にマッピングが存在することが重要である。以下の説明では、νｉまたはσｉのいずれかが、ソルバによってモデル化された二値変数を示すために使用される場合があることに留意されたい。
【0044】
[0051]二値変数σｉの値を光学的または電気的アナログ信号としてモデル化し、非デジタルハードウェアコンポーネントの組合せを使用してモデル化された各変数に対して上記の更新を実行する全アナログソルバを実装することができる。ソルバは、与えられた変数の割り当てを表す信号の初期集合を生成し、アナログ電子または光学ハードウェアに実装された相互作用ロジックを使用して計算されたフィードバック信号に基づいて、一連の反復ステップで新しい信号を生成する。イジング問題のためのソルバアーキテクチャの実装例については、後で詳しく説明する。
【0045】
[0052]本開示に従って配置され得るソルバ構成の多くの選択肢があり、各構成は、時間とともに生成される信号を「イジング」システムの全エネルギーを最小化する信号の集合に促すフィードバック信号を生成し、これは、所与の問題定義に対する変数の最適な割り当てにマッピングされ得る。
【0046】
[0053]本開示は、Ｎ個の変数のイジング問題（「スピン」と呼ばれることもある）にマッピングされ得る組合せ最適化問題を解くための新規なアーキテクチャを提供し、問題の変数は、Ｎ個の別個のハードウェアチャネルの集合によってモデル化され、所与の問題定義に従って変数の相互作用をモデル化する信号相互作用ロジックによって提供されるフィードバックに基づいて反復的に更新される。このシステムは、光学およびアナログ電子領域でのみ発生し、信号相互作用ロジックは、光学またはアナログ電子ハードウェアのいずれかによってモデル化され得る。これは、図１を参照して、さらに詳細に説明される。
【0047】
[0054]図１は、組合せ最適化問題のためのアナログソルバの例示的な概略ブロック図である。このアーキテクチャは、Ｎ個のチャネル１０２を含み、各チャネルｉは、各チャネルについて導出されたフィードバック方程式に従ってチャネルフィードバック信号を計算するように構成され、各チャネルのフィードバックは、問題を定義する「エネルギー」関数を最小化する変数の集合になりやすい。
【0048】
[0055]第１のチャネル１０２は、正または負の「スピン」値のいずれかを取るイジング変数σ１に対応するモデリング信号ｘ１を計算するように構成され、モデリング変数ｘ１は、最適化の各反復で受信されたフィードバックに基づいて更新される。モデル化される変数σは二値である可能性があるが、モデリング信号ｘは、変数の２つの可能な二進値の間で変化し得るソフトな値を取り得ることに留意されたい。次に、各チャネルによる寄与を決定するプロセスについて説明する。各チャネルは、関数へのそれぞれの寄与を計算するために同じステップを実行するハードウェアコンポーネントを含むことに留意されたい。
【0049】
[0056]各チャネル１０２は、信号生成器１００、スプリッタ１０６、および信号相互作用ロジック１０４を含み、これらの各々は、１つまたは複数のハードウェアコンポーネントを含み得る。本明細書においてこの文脈で使用される「ロジック」は、デジタルロジックを指すのではなく、むしろアナログまたは光ハードウェアを使用して実行される信号演算を指すことに留意されたい。信号生成器１００は、σｉのモデリング信号を生成し、信号の測定可能な特性は変数σｉの二進値を表す。信号は、例えば、レーザのような光源によって生成される光信号であってもよい。光変調器は、変数σｉをモデル化するために、光信号の特性を変調するために使用することができる。二値変数σｉが振幅のような特性の値で符号化されるためには、可能な変調された特性値（振幅）と二進値（例えば１と－１）との間のマッピングが定義されなければならない。例えば、ｘｉは［－ａ，＋ａ］の間の範囲にあり、ａはある定数であり、正の振幅はイジング変数σ＝１にマッピングされ、負の振幅はイジング変数σｉ＝－１にマッピングされる。変数σｉをモデリングする変調信号（これを本明細書では「モデリング信号」ｘｉと呼ぶ場合がある）が生成されると、この信号は、スプリッタ１０６を適用することによってコピーすることができ、同じ変数νｉを符号化するモデリング信号ｘｉの複数のインスタンスを生成し、図１の矢印で示すように、他のチャネルに伝達することができる。
【0050】
[0057]信号相互作用ロジック１０４は、変数のベクトルを表す複数のモデリング信号を受信し、各信号は、それぞれのチャネルｊのスプリッタ１０６から受信される。相互作用ロジック１０４は、モデリング信号ｘｉを、モデリングされた変数の重み付き和を表す信号に結合するベクトルごとの乗算器を含み、重みは、イジング問題のスピン相互作用を定義する行列Ｊの関連要素に対応する。ベクトルごとの乗算を実行するために使用され得る様々な可能なハードウェア構成が存在する。本明細書で開示する一例は、波長選択スイッチ（ＷＳＳ）である。これについては後でさらに詳しく説明する。光ベクトルごとの乗算は、代替的に、空間光変調器（ＳＬＭ）、リング共振器またはマッハツェンダー干渉計（ＭＺＩ）、またはそのような技術や他の適切な光学コンポーネントのいくつかの組合せを含む他の既知の光技術によって実施されてもよい。別の選択肢として、ベクトルごとの乗算演算は、例えば、メモリスタを使用することによって、アナログ電子領域で（すなわち、電気信号を使用して）実施することもできる。
【0051】
[0058]図１は、各チャネルのための別個の相互作用ロジック１０４を示すが、この相互作用ロジックは、各チャネルのための信号相互作用を計算するための別個のハードウェア構成要素を必ずしも含まないことに留意されたい。いくつかの実施形態では、システム全体の相互作用ロジックは、グローバルなベクトルごとの行列乗算器を含む。この場合、所与のチャネルのためのベクトルごとの乗算を実行するベクトルごとの乗算器の部分は、図１に示されるそのチャネルのための相互作用ロジックに対応する。例えば、各相互作用ロジックは、図９の例を参照して後で詳細に説明するように、２Ｄ空間光変調器の異なる行に対応してもよい。より一般的には、異なる個々の相互作用ロジックブロック１０４を実装するグローバルなベクトルごとの行列乗算器の物理的部分は、互いに重なっても重ならなくてもよいことにも留意されたい。一方、他の実施形態では、別個のハードウェア構成要素が、各チャネルについてベクトルごとの乗算を実行するように構成される。両方のアーキテクチャの効果は同じである。これらの代替アーキテクチャの両方について、図８および図９に示すベクトルごとの（またはベクトル行列）乗算器の波長選択スイッチ実装の文脈で、以下により詳細に説明する。
【0052】
[0059]フィードバック信号は、フィードバック経路１０８に沿って信号生成器１０４に戻され、信号生成器１０４は、システムのハードウェアに従って新しい信号を決定する。更新された信号は、例えば、フィードバック信号を変調器に渡して光源からの入力信号を変調し、その結果生じる光信号をフォトダイオードで検出することによって生成することができる。あるいは、一部の実施形態では、フィードバック信号を符号化したアナログ電子信号を、アナログ電子部品を用いて、例えば、メモリスタを用いて直接生成してもよい。いずれにせよ、システムは、時間の経過とともに、与えられた問題式のエネルギー関数を最小化する変数の最適割り当てにマップされる安定状態になるように設計される。
【0053】
[0060]各チャネルは、変数の値への特定の割り当てに対応するすべての信号が安定状態に達するまで、上述したのと同じ方式に従って信号を更新する。任意の数の変数σ１，．．．，σＮは、Ｎ個のチャネルを設定し、各信号をＮ個の同一コピー信号に分割し、各チャネルに１つずつ送信することによって、この方法でモデル化することができる。
【0054】
[0061]図２は、光ハードウェアに実装されたＮ個の変数のイジング問題のための光ソルバの例を示す。ソルバは、Ｎ個のチャネル１０２を含み、各チャネルｉは、変数ｉに割り当てられた二進値を表す信号を生成するために、信号生成器１００、スプリッタ１０６および相互作用ロジック１０４を実装するハードウェアを含む。わかりやすくするために、図２は単一のチャネル１０２のみの概略図を示しているが、これは各チャネル１０２で重複してもよい。モデリング信号は、図２のスプリッタ１０６によって他のすべてのチャネルに出力され、他の各チャネルから相互作用ロジック１０４で受信されることに留意されたい。
【0055】
[0062]各チャネル１０２は、システムが、解決されるべき最適化問題に従って変数の最適な割り当てを表す、安定した状態セットに落ち着くまで、フィードバック信号に従って更新されたモデリング信号ｘ１を反復的に生成する。上述したように、信号の更新は、例えば、更新方程式によって与えられる：
【0056】
【数5】
【0057】
[0064]初期信号は、所定のチャネルによってモデル化された変数σｉの初期二進値を表すスピン生成ハードウェア３００によって生成される。「スピン」は、本明細書では、イジングシステムの二値変数を表す信号を指すために使用され、スピンの量子力学的定義と混同されるべきではないことに留意されたい。スピン生成ハードウェア３００のハードウェアコンポーネントの実装例を、図３を参照して以下に詳細に説明する。スピン生成器によって出力された信号は、スピン生成ハードウェア３００から破線の矢印で示すように、電気信号としてスプリッタ３２６に送られ、スプリッタ３２６は信号を２つの異なる経路で送る。
【0058】
[0065]この実施形態では、スピン生成器３００は、図１の信号生成器１００の一部のみを含むことに留意されたい。図２のアーキテクチャ例では、信号生成器１００は、スピン生成器３００、光源３０２、変調器３０４、スプリッタ３２６、増幅器３２４、および電気信号を合計するアナログハードウェア３２２を含む。上述したように、信号生成器１００の出力はモデリング信号ｘ１である。
【0059】
[0066]第１の経路に沿って、信号は、変調器３０４において、特定の波長のレーザである光源３０２の出力と組み合わされ、レーザビームを変調し、それによって、図１を参照して上述したように、モデリング信号ｘｉを生成する。変調器３０４は、変数をモデル化する電気信号を符号化するために光信号を変調し、例えば、光信号の振幅、その位相、またはそれらの組合せを使用することができる。上述したように、検出された光信号が正負両方の値を符号化できるように、光信号の測定可能な特性と二値変数との間にマッピングが存在すべきである。これは、受信信号の位相と周波数情報、およびその強度を測定するコヒーレント検出方式を使用して行うことができる。信号の振幅とその符号のみが検出の対象である場合、直接検出の形式を使用することができ、これは本明細書において「差動検出」と呼ばれることがある。これについては以下に詳述する。
【0060】
[0067]変調器３０４は、変調信号を１からＮのスプリッタ３０６に送り、このスプリッタ３０６は、同一の光信号をシステムのＮ個のチャネルの各々のベクトルごとの乗算器３１４（ＶＶＭ）に伝達する。図２の例では、ベクトルごとの乗算器３１３は波長選択スイッチとして実装されている。しかしながら、他の実施形態では、これは、リング共振器、空間光変調器またはマッハツェンダー干渉計（ＭＺＩ）、またはアナログ電子部品によって実装されてもよい。ＷＳＳを採用するような実施形態では、各チャネルで生成される信号は、Ｎ個の信号間の干渉を引き起こすことなくベクトルごとの乗算における異なる信号の相互作用を可能にするために、異なる波長（すなわち色）で生成されてもよいことに留意されたい。ＶＶＭ３１４は、変数｛σ１，．．．，σＮ｝のベクトルを表す信号｛ｘ１，．．．，ｘＮ｝の入力集合に、行列Ｊの要素の対応する部分集合を乗算するように構成されている。所与のチャネルｉについて、波長選択スイッチは、各反復ｋについて式１のベクトルごとの項を表す信号、すなわちΣｊＪｉｊｘｊ［ｋ］を計算し、ここでｘｊ［ｋ］はチャネルｊについて信号生成ハードウェアで生成された信号である。ベクトル乗算を計算する波長選択スイッチの動作については、後でさらに詳細に説明する。
【0061】
[0068]実施形態では、ＶＶＭ３１４によって出力された信号は、図２の破線で示されているように、光学領域に留まり、この信号は、光検出器で検出することによって電気信号に変換される。図示の実施形態では、直接検出方式が使用され、和を実行するように構成されたアナログハードウェア３１０を使用して差動検出を達成するために、適応直流（ＤＣ）項３１２が追加される。増幅器３１６は、適応された和に定数βを乗算し、フィードバック経路１０８に出力する。式１の
【0062】
【数6】
【0063】
[0069]第２の経路では、信号ｉは増幅器に出力され、増幅器は、式（１）に示す定数αによる変数σｉの乗算を表す電気信号を増幅する。この信号は、アナログ加算ハードウェア３１８からフィードバック経路を経由して伝達される和
【0064】
【数7】
【0065】
【数8】
【0066】
[0070]スピン生成器３００による更新方程式の評価の一例を、図３を参照して以下に説明する。このプロセスは各反復で繰り返され、各チャネルの更新された信号を生成し、Ｎ個の信号の相互作用を決定し、新しいフィードバック信号を計算し、システムが安定するまで繰り返される。
【0067】
[0071]図２において一緒に動作する複数の構成要素は、図１を参照して上述した一般的な信号生成器１００に対応することに留意されたい。この実施態様における信号生成器１００は、スピン生成ハードウェア３００、レーザ３０２、アナログ加算ハードウェア３２２、スプリッタ３２６および変調器３０４のすべてを含む。同様に、図１を参照して上述した相互作用ロジック１０４は、波長選択スイッチに実装されたベクトルごとの乗算器と、フィードバック信号の加算および増幅を実行するアナログハードウェアコンポーネントの両方を含む。他の実施形態は、例えば、光検出器および変調器を含む、光またはアナログ信号に対する演算を実行するためのさらなるハードウェアを含み得る。
【0068】
[0072]各チャネルｉは、そのチャネルの信号の更新を並列に計算するハードウェアで実装される。更新は、例えば、Ｍ回の反復の所定の停止点の後にシステムが停止されるまで継続される。あるいは、信号を定期的に測定し、その後の測定で変化が見られない場合にシステムを停止することもできる。近似解は、システムが安定したとき、すなわち、生成された信号によってモデル化された変数の集合が、ある反復から次の反復まで一定のままであるときに発見される。この安定した信号の集合は、与えられたイジング問題の解を近似するＮ個の変数の割り当てに直接マッピングすることができる。
【0069】
[0073]図２に示す例示的な実施形態は、空間分割多重アーキテクチャ（すなわち、各変数が別々のハードウェアで生成される）を使用し、二値変数は信号振幅に符号化され、ＶＶＭは複数の並列ハードウェアチャネルの各々において波長選択スイッチによって実装される。しかし、代替の実施形態では、位相または周波数など、信号の異なる測定可能な特性を使用して変数をモデル化することができ、リング共振器、マッハツェンダー干渉計、またはメモリスタを使用して実装された電子ＶＶＭなど、他のタイプのＶＶＭを使用することもできる。また、時分割多重信号を結合するために遅延線を使用して、空間分割の代わりに時分割を使用することも可能である。
【0070】
[0074]図３は、スピン生成ハードウェア構成要素３００の一例の概略ブロック図である。指定波長の連続波レーザなどの光源４００は、光信号を生成するために使用され、光信号はスピンを生成するためにスピン生成変調器４０２に渡される。前の反復で生成された信号に基づく電気フィードバック信号も、相互作用ロジック１０４のＶＶＭ３１５から変調器４０２に受け取られる。変調器４０２は、電気フィードバック信号に従って光源４００からの光信号を変調する。変調器４０２は、マッハツェンダー変調器であってもよく、干渉される入力信号をそれ自体で分割する。フィードバック信号に対応する値は、変調器によってアームの１つ内に設定され、変調器４０２の出力は、干渉電界の同相成分、言い換えれば、フィードバック信号に適用される式１のｃｏｓ（）関数であり、光検出器４０４で光強度として検出されるとｃｏｓ２（）となる。信号生成光検出器４０４での検出は、得られた信号をアナログ領域に戻す。光検出器４０４を用いて光強度を直接検出する場合、信号は常に正であり、ｃｏｓ２（）関数によって与えられる［０，１］の範囲内にある。方程式の加法項、すなわち式１の－１／２項を表す適応ＤＣ項４０６は、イジング変数を適切にモデル化するために、信号の範囲を［－１／２，１／２］に変換するために正信号に加えられることがある。次に、出力信号は、図２を参照して上述したように、他のチャネルのＶＶＭに入力されるモデリング信号ｘｉを生成するためにスプリッタ３０６に渡される。
【0071】
[0075]図３に示す例のソルバの各反復の間、信号は光領域からアナログ領域へ、およびその逆へ２回変換される。図３に示すように、フィードバック信号は最初に信号生成器３００によって光信号に変換され、次に光検出器で検出されてアナログ信号に戻される。変調器は、アナログ領域から光領域に信号を変換して光信号相互作用ロジックで処理するために光源とともに使用され、光検出器は相互作用ロジックの出力信号を光領域からアナログ領域に戻すために使用される。
【0072】
[0076]図４は、Ｎ個のハードウェアチャネルを含む、上述した図２に示される例示的なソルバにおける１つのチャネルについての、光領域とアナログ領域との間の信号の変換の概略ブロック図である。図３を参照して上述したスピン生成コンポーネント３００は、図４において「信号１」と示される電気信号を出力する。本例のソルバでは、信号相互作用ロジックは、光信号に作用する、例えば波長選択スイッチの形態の光ベクトルごとの乗算器３１４を含む。しかしながら、上述したように、任意の所与のチャネルの相互作用ロジックは、複数のチャネルのためのフィードバック信号のベクトルを生成するためのベクトル行列乗算を実装する全体的な信号相互作用ハードウェアの一部を含み得る。このアーキテクチャの一例を、図９を参照して以下に詳細に説明する。最初のアナログ光信号変換５００は、光信号を信号間相互作用ロジック５０２に渡すために信号が生成された後に発生する。ベクトルごとの乗算結果は、光検出器３０８での検出によって評価され、図２に示すように、信号が光領域に変換され、さらなる算術演算がアナログ領域で実行される。
【0073】
[0077]図４の信号間相互作用ロジック５０２は、光ＶＶＭ演算と電気フィードバック信号への変換の両方を含む。出力された電気的フィードバック信号は、次に第２のアナログ光信号変換５０４に渡され、ここで電気的フィードバック信号は変調器によって光信号に戻される。図３と図４では、これはスピン生成ハードウェア３００内で発生し、スピン生成光源４００と変調器４０２が電気フィードバック信号を光信号に変換する。図４の信号生成５０６は、図３で説明した完全な信号生成プロセスではなく、信号を光領域からアナログ領域に変換する信号生成ハードウェア３００によって生成された信号の検出のみに対応する。
【0074】
[0078]上述した例示的な実施形態では、信号相互作用ロジック５０４のベクトルごとの乗算動作は、例えば、後述する波長選択スイッチによって、光領域で実装されることに留意されたい。しかしながら、他の実施形態では、信号相互作用はアナログ電子領域で実施されてもよい。同様に、いくつかの実施形態では、信号の加算のような他の算術演算は、アナログ電子領域ではなく光学領域で実行されてもよい。したがって、図４に示すプロセスは、図２に示すソルバを実装するために使用される特定のハードウェア構成に固有のものである。
【0075】
[0079]上述したように、本明細書で説明するアーキテクチャの利点は、「空間分割」多重化アーキテクチャを使用することであり、これは、Ｎ個の変数のシステムが、各変数に対して別々のハードウェアを使用してモデル化されることを意味する。対照的に、いくつかの最先端のソルバは、時分割多重化アーキテクチャを使用する。図５Ａは、時分割多重化を使用するアーキテクチャの例を示している。時分割多重化アーキテクチャは、システムのすべての変数を表す信号を生成するための１組の信号生成ハードウェア５１０ａと、信号相互作用ロジックを実装する１つのハードウェアと、を使用する。信号生成ハードウェア５１２ａは、時間間隔で信号を生成し、相互作用ロジックは、相互作用ハードウェアで異なる時間に受信された信号に遅延を適用することにより、時間分割された信号の相互作用を可能にする。ソルバの時間的複雑性は、モデル化される変数の数によって線形に増加するので、このアーキテクチャは、より大きなシステムに対してより遅い。
【0076】
[0080]対照的に、図５Ｂに示す空間分割多重化アーキテクチャは、各変数に対して個別の信号生成器５１０ｂを含み、すべての信号は信号相互作用ハードウェア５１２ｂで相互作用する。別個の物理的生成により、すべての信号を同時に生成することができ、システムは、反復ごとに準一定の時間を維持しながら、より多くのハードウェアを追加するだけで、多数の変数にスケールアップすることができる。図５Ｂの完全ソルバシステム用の信号相互作用ロジックハードウェア５１２ｂは、単純化のために単一のブロックとして示されているが、相互作用ロジックは、代替的に、上述された図２に示されるように、各チャネル用の別個のハードウェアＶＶＭを含むように実装されてもよい。
直接検出
[0081]図４を参照して説明したように、実施形態では、信号は、ソルバの動作中に光領域とアナログ領域との間で変換される。例えば、式１を表す電気フィードバック信号は、光検出器で検出することによって完全な式を評価する前に、ｃｏｓ（）関数を適用するマッハツェンダー変調器を使用して光領域に変換される場合がある。光強度の直接検出は、信号の２乗を測定するので正のみである。しかし、本明細書でイジング変数をモデル化する信号のように、正または負の値を取る信号を扱う場合、アナログ電子領域と光学領域との間の変換は、信号の符号情報を保持する必要があり、または、信号が例えば正の値に制限されるような変換は、関連する領域での後続の操作によって修正されるべきである。「実数値」は、本明細書では、複素平面の実軸に沿った値を有する光信号、すなわち、正または負の実数値を取り得る信号を指すために使用され得る。
直接検出
[0081]図４を参照して説明したように、実施形態では、信号は、ソルバの動作中に光領域とアナログ領域との間で変換される。例えば、式１を表す電気フィードバック信号は、光検出器で検出することによって完全な式を評価する前に、ｃｏｓ（）関数を適用するマッハツェンダー変調器を使用して光領域に変換される場合がある。光強度の直接検出は、信号の２乗を測定するので正のみである。しかし、本明細書でイジング変数をモデル化する信号のように、正または負の値を取る信号を扱う場合、アナログ電子領域と光学領域との間の変換は、信号の符号情報を保持する必要があり、または、信号が例えば正の値に制限されるような変換は、関連する領域での後続の操作によって修正されるべきである。「実数値」は、本明細書では、複素平面の実軸に沿った値を有する光信号、すなわち、正または負の実数値を取り得る信号を指すために使用され得る。
【0077】
[0082]正負の値の検出を可能にする１つの可能な検出方式は、コヒーレント検出であり、これは、正または負のいずれかである受信光信号の振幅および位相情報を測定する。しかし、コヒーレント検出の欠点は、光強度の直接検出よりも実装が複雑であることである。コヒーレント検出方式は多くの場合、デジタル信号処理を必要とする。コヒーレント検出を実行するためにデジタル領域に変換し直すと、信号の伝送速度など、光領域およびアナログ電子領域で信号を処理する利点の一部が失われるか、または低下する。
【0078】
[0083]別の検出方法では、コヒーレント検出のようなシステムの複雑さを必要としない直接検出、すなわち光強度の検出が使用される。直接検出は、アナログ電子領域で正のみの信号を測定し、これは、その後、信号の範囲を正または負の値を許容するように補正する適応項を加算または減算することによって、アナログ電子領域でオフセットされ得る。これは「差動検出」とも呼ばれ得る。同様の検出方式は、実数値である二値位相シフトキーイング信号を検出するために電気通信で使用される。
【0079】
[0084]この直接検出方式の概略図を図６に示すが、これは光ベクトル乗算器の出力のような実数値の光信号を検出するためのものである。これを、上述した図２のソルバの適用例について説明する。正または負の値を取り得る実信号７００、例えば式１を表すためにスピン生成器３００によってチャネルｉで生成されたモデリング信号ｘｉ（［－１／２、＋１／２］の範囲の値を取り得る）は、例えば変調器３０４によって、少なくとも１つのアナログ光変換によって最初に光領域に変換される。チャネルのモデリング信号および他のチャネルから受信したモデリング信号に対して光ベクトルごとの乗算が実行され、入力ベクトルにチャネルｉの相互作用に対応する関連する重みが乗算され、入力信号の重み付き和を表すチャネルの光信号が出力される。
【0080】
[0085]この信号は、光検出器３０８で検出することによりアナログ信号に変換される。しかし、光検出器４０４は光強度を測定するため、検出される信号は正の値のみに制限され、負の値を持つことはできない。これを修正するために、ＶＶＭ動作の出力信号は、図２に３１２として示されているＤＣオフセット項を追加することによって、正または負の値を許容するように補正される。このＤＣオフセット項は各スピンに固有であり、システムの初期化時に評価および設定される。
【0081】
[0086]これにより、アナログ領域で信号を調整することにより、ソルバが必要とする正と負の信号の測定が可能になる。この差動検出方式はコヒーレントなものよりも単純であり、信号を光領域からアナログ電子領域に直接変換するために簡単に実装することができる。しかし、ＶＶＭ出力では、与えられた入力信号の波長が異なる場合、すべての信号の経路長が一致するように注意を払う必要がある。信号のインコヒーレントな加算については、波長選択スイッチの動作との関連で以下に詳しく説明する。
【0082】
[0087]この差動検出方式は、本ソルバアーキテクチャに関連して上述されているが、適応オフセット項を有する直接検出は、実の正および負の値を取る光ベクトルごとの乗算演算が実装され得る任意のアプリケーションに使用され得る。例えば、これは、入力ベクトルがネットワーク重みによって乗算され得るディープニューラルネットワークのような機械学習アプリケーションにおいて使用され得る。この差動検出方式は、以下に詳述する空間光変調器（ＳＬＭ）、リング共振器、または波長選択スイッチなど、様々なタイプの光ＶＶＭを使用するアプリケーションに適用することができる。この差動検出方法は、光領域で演算を行うことができるため、デジタル演算に比べて大幅な高速化を実現できる一方、コヒーレント検出方式の難しい実装を必要とせずに、所望の範囲の実数値信号をモデル化できるという利点がある。このような差動検出方式は、波長選択スイッチまたはリング共振器ＶＶＭのように、ＯＶＭの入力信号に異なる波長が使用される場合、システムの位相感度を必要とせずに実装することができる。
【0083】
[0088]図６Ａは、差動検出方式の概略ブロック図であり、一定のアナログ電子ＤＣ項７１２が光検出器によって検出された信号に加えられる。光信号は光検出器７１０で受信され、この光検出器は、上述したように、実光信号７００を正の値のみを取り得るアナログ信号７０４に変換する。別個のアナログ定数ＤＣ項７１２が、図６に示される減算器７０８の可能な実装である微分器７１４において正のアナログ信号から減算される。アナログＤＣ項７１２は、図６Ｂを参照して後述するように、１つまたは複数のアナログ電子部品を使用して、または光信号を使用して生成することができる。アナログ微分器は、電子工学の技術分野において周知であり、本明細書ではこれ以上説明しない。微分器の出力は、正のみの信号７０４から一定のＤＣ項を減算することによって得られる実アナログ信号７１８であり、これにより、アナログ信号は、正または負の値を取り得る所望の範囲に入る。前述したように、各ＤＣ項は光検出器に固有であり、システムが初期化されるときに評価され、設定されてもよい。ＤＣ項は、等価的に、正のアナログ信号７０４に加算される負の値のアナログ信号として、または正のアナログ信号７０４から減算される正の信号として見ることができることに留意されたい。
【0084】
[0089]図６Ｂは、光オフセット差動検出方式の概略ブロック図であり、この方式では、任意の２つの光信号７２０ａ、７２０ｂが一対の光検出器７２２ａ、７２２ｂで検出され、アナログ微分器７２６を使用して、それぞれの検出器７２２ａ、７２２ｂで生成された電気信号７２４ａ、７２４ｂの間で差が取られ得る。これは既知の検出構成であり、上述の差動検出方法を実施するために使用されてもよく、実数値信号、例えば光ＶＶＭの出力が第１の光信号７２０ａに符号化され、定数値が第２のオフセット光信号７２０ｂに符号化され、これは定数値によって変調された別個の光源によって提供されてもよい。各信号は、正の値のみを取るアナログ電子信号７２４ａ、７２４ｂに変換されてもよい。これらは、評価すべき実数値信号を検出して得られたアナログ信号１から一定のオフセットアナログ信号２を減算する微分器７２６で結合される。この差動検出方式の効果は、図６Ａの差動検出と同じであり、この場合、ＤＣオフセット項は、一定の光信号７２０ｂの検出によって得られる。ただし、この差動検出の実施では、アナログ電子部品を使用してアナログＤＣオフセット信号を直接生成する代わりに、第２の光オフセット信号を生成する必要がある。
【0085】
[0090]説明した実施形態では、ソルバは、イジング問題に容易にマッピングされるＱＵＢＯ問題を解く方法として、イジング変数を表す正負の「スピン」信号、例えば－１／１の形でｘおよび対応するフィードバック信号をモデル化することに留意されたい。フィードバック信号の符号は、ハミルトニアンのエネルギーを減少させるためにモデリング信号ｘを駆動する方向を表す。しかし、他の実施形態では、純粋に正の信号を使用することを排除するものではない。代わりに、行列Ｊは、正および負の重みを含み得る。このような実施形態では、ＤＣオフセット３１０、３２０は必ずしも必要ではない。例えば、ＱＵＢＯ変数１／０を直接モデル化してもよい。この場合、直接検出によって生成された正信号を補正する必要はない場合がある。
光ＶＶＭ用波長選択スイッチ
[0091]上述のように、各チャネルは、相互作用ロジック１０４の一部として、それぞれのベクトルごとの乗算器を実装することができる。様々な可能なベクトルごとの乗算器構成が、本明細書に開示されるソルバアーキテクチャとともに使用され得る。一部のＶＶＭは、空間光変調器、リング共振器、マッハツェンダー干渉計など、完全に光学領域で実装されてもよい。他のＶＶＭは、例えば、電気信号の重み付き和を計算するためにメモリスタを使用するなど、アナログ電子領域で実装されてもよい。
光ＶＶＭ用波長選択スイッチ
[0091]上述のように、各チャネルは、相互作用ロジック１０４の一部として、それぞれのベクトルごとの乗算器を実装することができる。様々な可能なベクトルごとの乗算器構成が、本明細書に開示されるソルバアーキテクチャとともに使用され得る。一部のＶＶＭは、空間光変調器、リング共振器、マッハツェンダー干渉計など、完全に光学領域で実装されてもよい。他のＶＶＭは、例えば、電気信号の重み付き和を計算するためにメモリスタを使用するなど、アナログ電子領域で実装されてもよい。
【0086】
[0092]本明細書で開示されるソルバアーキテクチャのいくつかの実施形態で使用するために本明細書で開示される光ＶＶＭ（ＯＶＶＭ）の一例は、波長選択スイッチである。（ＷＳＳ）である。ＷＳＳは、電気通信アプリケーションで使用され、すべての信号が同じ電力で伝送されることを保証するために、異なる波長の信号が独立して最適化されることを可能にし、また、異なる波長の信号が単一の光ファイバ内で一緒に結合されることを可能にし、またはその逆に、伝送ノードにおける加算機能またはドロップ機能のために結合されることを可能にする。
【0087】
[0093]光ベクトル乗算のためのＷＳＳの実装は、ＷＳＳが個々の波長を減衰させる（重くする）ので積関数をエミュレートする能力を有し、異なる波長を単一のファイバに結合する能力によって達成される加算関数が、その後少なくとも１つの光検出器によって検出されるという事実に基づいている。
【0088】
[0094]図７は、光信号に一定の係数（重み）を乗算する原理を示す概略ブロック図である。図７の例では、光源８０２ａ、８０２ｂから２つの別個の光信号が供給される。これらの信号のパワーまたは強度は、入射光の強度に依存する電気信号を生成する光検出器８０６で測定することができる。したがって、光信号によって符号化された値の乗算は、減衰器または損失誘発コンポーネント８０４ａ、８０４ｂを適用することによって実施することができ、減衰器または損失誘発コンポーネント８０４ａ、８０４ｂはそれぞれ、設定可能な定数係数によってそれぞれの光信号の強度を低下させる、つまり、光検出器またはカメラ８０６で検出されたときの得られるアナログ信号は、入力信号に適用される「重み」と解釈され得る、適用される減衰器の損失係数でスケーリングされる。２つの別々の光信号が同じ波長を共有する場合、光検出器でのそれらの加算はコヒーレント（すなわち、２つの光信号は電界で加算される）になるが、２つの別々の光信号が異なる波長を持ち、その差が光検出器の帯域幅よりもはるかに大きい場合、光検出器でのそれらの加算はインコヒーレント（すなわち、２つの光信号は電力で加算される）になる。
【0089】
[0095]ベクトル行列の乗算は、以下の形式の一連のベクトルごとの積に分解することができる：
【0090】
【数9】
【0091】
[0097]図７の構成は、入力ベクトルを表す入力信号の集合に減衰器８０４ａ、８０４ｂを適用することによって、光ベクトルごとの乗算の重み付き和を実装するために使用することができ、減衰器は、重み行列の適切な要素に対応するように設定される。重み付けされた入力の合計は、光信号のコヒーレントまたはインコヒーレントな加算を実行できるカメラまたは他の検出システム８０６で信号を結合することによって実行される。コヒーレント加算は、光信号が同じ波長を共有する場合に必要であるが、正しい出力信号を計算するためには、すべての信号を正確に位相整合させる必要がある。これは光波長では難しい。インコヒーレント加算の場合、光検出器は、入力信号が、所定の検出器の帯域幅によって与えられる許容範囲内で、同じ経路長を有することを必要とする。光検出器での直接検出によるインコヒーレント加算では、入力光信号の合計パワーまたは強度が計算される。これについては、波長選択スイッチにおけるインコヒーレント検出について後述するが、この場合、加算されるそれぞれの信号は異なる波長である。
【0092】
[0098]次に、上記の原理に基づいてベクトルごとの乗算またはベクトル行列乗算を実行する波長選択スイッチの動作を、図８を参照して説明する。
【0093】
[0099]入力ベクトルｖは、異なる波長の光信号８００の集合によって表され、これは、例えば、図２に示され上述されたようなソルバのＮ個のチャネルから受信されたモデリング信号｛ｘ１，．．．，ｘＮ｝の集合である。入力ベクトル用の別個の光信号８００を運ぶ入力ファイバ８０８は、図８の右下に示されている。
【0094】
[00100]注：図示を簡単にするために、ソルバの１つのチャネル１０２のみのファイバ８０８、８１８が図８に示されている。しかしながら、このようなソルバにおいてＷＳＳが適用される実施形態では、各チャネル１０２に対してファイバ８０８、８１８の別個のセットが存在し得る。入力ファイバの各セットは、所定のチャネル１００のスプリッタ１０６からの入力に対応する。出力ファイバ８１８は、それぞれのベクトルごとの乗算器の出力を提供する（したがって、全体としてベクトル行列乗算を実行する）。この場合、ＳＬＭ８１０の要素の異なる対応する部分集合は、異なるチャネルの相互作用ロジック１０４を実装する。実施形態では、これらの要素は、ＳＬＭ８１０の同じ物理的プレート（例えば、同じガラス片またはプラスチック片）上の異なる位置の異なる画素として実装することができる。
【0095】
[00101]重み行列Ｑの対応する要素は、空間光変調器（ＳＬＭ）８１０に実装され、その一例は、シリコン上の液晶空間光変調器（ＬＣｏＳ－ＳＬＭ）であり、上述したように、各入力光信号８００を特定の係数で変調する。この場合、信号は入力の波長に依存する係数で変調され、ＳＬＭ８１０の各列は異なる入射波長に対応する。入力信号８０８は、各信号がそれぞれの波長に対して正しい水平位置でＳＬＭに到達するようにレンズを通る。
【0096】
[00102]所与のチャネルの出力信号は、変調された光信号を光検出器８２０で検出することによって取得され、単一ビーム８１８に結合され、その後光検出器８２０で検出される。各々が異なる波長を有する様々な光信号を光検出器８２０で単一のビームに結合することは、波長分割多重（ＷＤＭ）と呼ばれることがある。これは、１つまたは複数のレンズ８１６および／または分散素子８１４（例えば、プリズムまたは回折格子などの回折素子）の配置によって促進され、一方、ＳＬＭは、各個々の波長に対する独立した重みを保証する。
【0097】
[00103]光検出器８２０は、異なる波長の様々な構成の光信号のインコヒーレント加算を実行する。インコヒーレント検出が構成信号の強度の和を計算するためには、結合されるそれぞれの信号の周波数の差が光検出器の周波数帯域幅よりもはるかに大きいことが保証されなければならず、これは、光検出器が信号間の相互作用から交差項を検出しないことを意味する。異なる波長の信号をインコヒーレント検出する場合、信号の位相を合わせる必要はない。これとは対照的に、同じ波長の光源を入力とするＶＶＭアーキテクチャを使用する場合は、検出器でコヒーレント加算を実行する必要があり、すべての信号が位相整合されている必要があるという困難な要件がある。
【0098】
[00104]図８に示すものと同様のアーキテクチャは、電気通信アプリケーション用の波長選択スイッチデバイスを構築するために知られており、通常、逆方向に実行され、入力は単一の光ファイバであり、出力は異なる光ファイバ内の異なる波長の光信号の集合である。そのシナリオでは、波長選択スイッチは、特定の伝送ノードで特定の波長を追加またはドロップするために使用されるだけでなく、信号が伝送スペクトル全体で平坦であることを保証する。
光ベクトル行列乗算のための拡張ＷＳＳアーキテクチャ
[00105]上述のように、本明細書で説明するイジング問題のためのソルバは、２つのアーキテクチャのうちの１つで実装することができる。第１では、図２に示すように、信号相互作用ロジック１０４は、波長スイッチのような別個のベクトル行列乗算器３１４として実装され、重み行列Ｊの１行によって与えられる重みのベクトルによって入力ベクトルを乗算するために上述のように動作する。各チャネルによって生成された信号は、１対Ｎスプリッタによって、他のすべてのチャネルでベクトルごとの乗算器に向けられる。
光ベクトル行列乗算のための拡張ＷＳＳアーキテクチャ
[00105]上述のように、本明細書で説明するイジング問題のためのソルバは、２つのアーキテクチャのうちの１つで実装することができる。第１では、図２に示すように、信号相互作用ロジック１０４は、波長スイッチのような別個のベクトル行列乗算器３１４として実装され、重み行列Ｊの１行によって与えられる重みのベクトルによって入力ベクトルを乗算するために上述のように動作する。各チャネルによって生成された信号は、１対Ｎスプリッタによって、他のすべてのチャネルでベクトルごとの乗算器に向けられる。
【0099】
[00106]しかし、第２のアーキテクチャでは、グローバルなベクトル行列乗算器（ＶＭＭ）が実装されてもよく、ソルバのチャネルはそれぞれ、入力ベクトルを形成するために、そのモデリング信号ｘｉをＶＭＭに提供し、完全な行列はこのＶＭＭに実装される。図１に示すソルバアーキテクチャはこのアーキテクチャに対応し、各チャネルの相互作用ロジック１０４は、各チャネルの個別の乗算器ではなく、グローバルベクトル行列乗算器の要素のそれぞれの部分集合（図示の例では行など）によって少なくとも部分的に実装される。実施形態では、これらの要素は、異なるピクセル位置で同じＳＬＭプレート（例えば、ガラスまたはプラスチックの同じ部分）に統合されてもよい。
【0100】
[00107]次に、ＷＳＳベクトルごとの乗算器のアーキテクチャを拡張してベクトル行列演算を実行するＷＳＳアーキテクチャの例を説明する。このアーキテクチャは、上述のベクトルごとのＷＳＳよりも多くのスピンを同時に処理できるという利点がある。
【0101】
[00108]図９は、第２のソルバアーキテクチャを実装するために使用され得る、ベクトル行列乗算を計算するための波長選択スイッチ構成例を示す図である。この波長選択スイッチは、異なる波長の信号の集合を生成する光源の入力アレイ９０８と、レンズレットアレイ９００と、修正ＳＬＭ９０２と、分散素子８１４（例えば回折素子）と、図９には示されていないが光検出器のアレイで検出される光信号の出力アレイと、を含む。
【0102】
[00109]この例のソルバアーキテクチャにおいてベクトル行列乗算に空間光変調器を使用するには、ベクトル行列乗算の全機能が達成されるように、ＳＬＭの垂直軸は、同じ波長であっても異なる重みを提供する必要がある。これは、行列乗算の場合、完全な出力ベクトルを生成するために、入力ベクトルを行列Ｑの各行と乗算する必要があるためである。ＳＬＭ９０８は、図８に示すものの修正バージョンであり、変調器のアレイがアレイ状に配列され、各変調器によって適用される損失は、入力に適用される行列の重みを反映する、すなわち、修正ＳＬＭの行は、行列Ｊの行の重みを符号化する。上述したように、各チャネルは、入力信号と行列Ｊの単一行とのベクトル乗算を含む出力を計算する。したがって、入力信号の各々は、一連のベクトルごとの乗算に対応するＳＬＭ９０２の各行に当たるように、垂直方向に広がるように処理される必要がある。
【0103】
[00110]グローバルＶＭＭを有する例示的なソルバアーキテクチャでは、単一の入力アレイ９０８が、各チャネルで生成されたモデリング信号ｘｉを含む。このベクトルは、信号の波長に対応するＳＬＭ９０２の水平方向にビームをコリメートしながら、信号を垂直方向に広げる特定の形状を有するレンズレットアレイ９００に通される。これにより、異なる波長のより多くの入力信号を単一のＳＬＭで処理することができる。図８のような単一レンズから図９のようなレンズレットアレイに移行することで、より多くの波長へのスケーリングが可能になる。レンズレットアレイは、両方向のビームのコリメーション特性を改善する。
【0104】
【0105】
[00112]ＳＬＭ９０２は、変調器の２Ｄアレイを含み、アレイの各要素は、受信した入力信号にそれぞれの重みを適用し、これは、図８のベクトルごとの乗算器について説明したＳＬＭとは対照的に、同じ「列」にあるが異なる垂直位置の要素が異なる値を有することを必要としない。特定の波長で、しかし異なる重みで変調された各重み付け信号、すなわちＳＬＭ９０２の列の各信号は、ＷＤＭ信号の組合せが光検出器アレイ９０４内の正しい光検出器で発生することを保証するために、異なる垂直位置で分散素子８１４から跳ね返る必要がある。分散素子は、回折格子またはプリズムなどの回折素子として実装することができる。
【0106】
[00113]実施形態では、出力信号は、そのビームによって表される出力ベクトル要素に対応する光検出器においてインコヒーレント加算を使用して検出されるように、信号を正しい垂直高さのビームに向けるために、１つまたは複数のレンズを介して要素８１４から向けられてもよい。例えば、分散素子８１４とシステムの末端の複数のチャネル（潜在的にファイバ）との間に別のレンズレットアレイを含めることもできる。
【0107】
[00114]光検出器アレイ９０４は、垂直アレイの光検出器の集合として配置され、分散素子８１４から向けられた各結合信号は、異なるチャネルの出力信号に対応する。
【0108】
[00115]上述したベクトル行列乗算器アーキテクチャを使用するソルバは、入力をコリメートするためにレンズレットアレイを使用する一方で、図９に示すような単一のハードウェア配置を使用して、すべてのチャネルに対するスピンの相互作用の同時処理を可能にする。これは、各ビームを複数のＳＬＭ９０２の構成に向けられる複数のビームに分割することによって、さらに多数の入力を可能にするようにさらに拡大することができる。
【0109】
[00116]波長分割多重を使用しない空間光変調器、リング共振器、およびマッハツェンダー干渉計など、多くの既存技術によって光ベクトル乗算も実装されている。このような技術は、例えば、Ｋ．Ｋｉｔａｙａｍａらによる、「Ｎｏｖｅｌ ｆｒｏｎｔｉｅｒ ｏｆ ｐｈｏｔｏｎｉｃｓ ｆｏｒ ｄａｔａ ｐｒｏｃｅｓｓｉｎｇ－Ｐｈｏｔｏｎｉｃ ａｃｃｅｌｅｒａｔｏｒ（データ処理のためのフォトニクスの新たなフロンティア－フォトニックアクセラレータ）」、ＡＰＬ Ｐｈｏｔｏｎｉｃｓ ２０１９、ｈｔｔｐｓ：／／ｄｏｉ．ｏｒｇ／１０．１０６３／１．５１０８９１２に詳細に記載されており、参照によりその全体が本明細書に組み込まれる。波長選択スイッチの実装は、ＳＬＭの空間変調とリング共振器の波長分割を組み合わせたものであるが、リング共振器の実装では、入力信号を一連のリング共振器に通す必要があるのに対し、ＳＬＭでは、各信号を１つの変調器に通すだけでよく、システム損失の点で有利である。ＳＬＭ ＶＭＭ実装では、波長分割は行わず、代わりに、単一の光源を使用し、光検出器でのコヒーレント加算を用いて出力アレイの各要素の重み付き和を計算する。波長選択スイッチは、これら両方の技術の利点を兼ね備えている。
【0110】
[00117]波長選択スイッチに関する上記の説明は、本明細書で説明するようなソルバアーキテクチャにおける実装に言及している。しかしながら、ベクトル行列乗算は、特に機械学習において、例えばニューラルネットワークの重みを入力ベクトルに適用するなど、多くの用途がある。本明細書で説明する波長選択スイッチは、このような用途に使用することができる。同様に、波長選択スイッチＶＭＭは、図５Ａに示される時分割多重アーキテクチャなどの他のソルバアーキテクチャに適用されてもよい。
【0111】
[00118]本明細書に開示された技術は、広範囲のアプリケーションに適用することができ、特に、本明細書に開示されたソルバの実装は、イジング定式化への既知の変換が存在する任意のＮＰ困難問題を解くために使用することができる。そのような問題のよく知られた例は、巡回セールスマン問題である。これは、他の分野の問題、例えば、分子の類似性を決定する問題にも使用することができ、分子のグラフ表現のグラフ類似性問題のＱＵＢＯ定式化への変換を見出す作業が行われている。この研究は、Ｈｅｒｎａｎｄｅｚ，Ｍａｒｉｔｚａらの「Ａ ｑｕａｎｔｕｍ－ｉｎｓｐｉｒｅｄ ｍｅｔｈｏｄ ｆｏｒ ｔｈｒｅｅ－ｄｉｍｅｎｓｉｏｎａｌ ｌｉｇａｎｄ－ｂａｓｅｄ ｖｉｒｔｕａｌ ｓｃｒｅｅｎｉｎｇ（３次元リガンドに基づく仮想スクリーニングのための量子的手法）」に記載されている。Ｊｏｕｒｎａｌ ｏｆ Ｃｈｅｍｉｃａｌ Ｉｎｆｏｒｍａｔｉｏｎ ａｎｄ Ｍｏｄｅｌｉｎｇ ５９．１０（２０１９）：４４７５－４４８５
[00119]上記の実施形態は、例として説明されたことが理解されよう。開示された技術の他の変形および適用は、本明細書における概念の開示が与えられれば、当業者には明らかになるであろう。
[00119]上記の実施形態は、例として説明されたことが理解されよう。開示された技術の他の変形および適用は、本明細書における概念の開示が与えられれば、当業者には明らかになるであろう。
【0112】
[00120]より一般的には、本明細書に開示される一態様によれば、関数を最適化する変数のベクトルの値を推定するシステムが提供される。関数は、複数の項の重み付き和を含み、各項は、ベクトルからの変数の対応する部分集合の積を含み、各項は、変数間の相互作用をモデル化する重みの行列からの対応する重みによって重み付けされる。システムは、互いに同時に動作するように配置された複数の並列ハードウェアチャネルを含み、各チャネルは、関数に対する変数のそれぞれの１つの寄与をモデル化するように配置される。並列チャネルのそれぞれは、それぞれの変数の値をモデル化する変調された特性を有するそれぞれのモデリング信号を生成するように構成されたそれぞれの信号生成器と、それぞれのモデリング信号のインスタンスを並列チャネルのそれぞれに供給するように配置されたそれぞれのスプリッタであって、各チャネルは、変数のベクトルをモデル化した信号のベクトルを受信する、スプリッタと、受信した信号のベクトルに、それぞれの変数と変数のベクトルとの間の相互作用をモデル化する重みの行列からの重みのそれぞれのベクトルを乗算するように構成されたそれぞれのベクトル乗算器を含み、それによって、それぞれのチャネルによってモデル化されたそれぞれの変数の寄与を表すそれぞれのフィードバック信号を生成するそれぞれの相互作用ロジックと、フィードバック信号をそれぞれの信号生成器に戻すように配置されたそれぞれのフィードバック経路であって、それぞれの信号生成器は、フィードバック信号に応じて、それぞれのモデリング信号を適応させるように構成される、フィードバック経路と、を含み、各チャネルにおいてそれぞれの信号生成器、スプリッタ、相互作用ロジックおよびフィードバック経路を含む各チャネルは、光学コンポーネントおよび／またはアナログ電子部品のみを使用して実装される。
【0113】
[00121]実施形態では、変数のそれぞれは、二値である。
【0114】
[00122]実施形態において、変数のそれぞれは、正または負の値を取り、関数は
【0115】
【数10】
【0116】
[00123]実施形態において、正および負の値は、＋１／２および－１／２、または＋１および－１である。
【0117】
[00124]実施形態において、変数のそれぞれはブールであり、関数は、
【0118】
【数11】
【0119】
[00125]実施形態において、各チャネルにおいて、それぞれのフィードバック経路は、それぞれの信号生成器に戻る前に、それぞれのフィードバック信号にそれぞれのノイズ成分を導入するように配置される。
【0120】
[00126]実施形態において、それぞれの変数の値をモデル化するために使用される各モデリング信号の変調された特性は、信号の振幅、または信号の位相、のいずれか一方を含む。
【0121】
[00127]実施形態において、各チャネルにおいて、それぞれの信号生成器によって生成されたそれぞれのモデリング信号は、光信号を含み、それぞれのベクトル乗算器は、光領域においてそれぞれの乗算を実行するように構成された光乗算器を含む。
【0122】
[00128]実施形態において、それぞれの信号生成器によって生成される光信号は、異なるそれぞれの光波長を有する。
【0123】
[00129]実施形態において、各ベクトル乗算器の光乗算器は、
－ 空間光変調器、
－ 波長選択スイッチ、
－ リング共振器、または
－ マッハツェンダー干渉計、
のいずれか１つを含む。
－ 空間光変調器、
－ 波長選択スイッチ、
－ リング共振器、または
－ マッハツェンダー干渉計、
のいずれか１つを含む。
【0124】
[00130]実施形態において、各相互作用ロジックは、アナログ電子形式でそれぞれの信号生成器へのフィードバック信号を生成するために、それぞれの光乗算器の光出力を検出するように配置されたそれぞれの光検出器をさらに含む。
【0125】
[00131]実施形態において、それぞれの信号生成器は、それぞれの光源と、それぞれのフィードバック信号に応じて、アナログ電子信号の形式でそれぞれのスピンを生成するように構成されたそれぞれのスピン生成器と、それぞれのスピンをそれぞれの光信号の変調された特性に変調するように配置されたそれぞれの変調器と、を含む。
【0126】
[00132]実施形態において、各チャネルにおけるそれぞれのスピン生成器は、さらなる光源と、それぞれのフィードバック信号に応じてさらなる光源からの光を変調するように配置されたさらなる変調器と、さらなる変調器からの変調された光を検出し、それに応じてスピンを生成するように配置されたさらなる光検出器と、を含む。
【0127】
[00133]実施形態では、各チャネルにおいて、それぞれの変数の値をモデル化するために使用される各信号の特性は、信号の振幅を含み、それぞれのスピン生成器は、スピンを表すために正レベルと負レベルとの間で振幅変調されたアナログ電子信号の形式でそれぞれのスピンを生成するように配置され、それぞれの変調器は、スピンをそれぞれの光信号の振幅に変調するように構成され、光信号の振幅およびそれぞれの光乗算器の出力は、正のレベルのみを取ることができ、負のレベルは取ることができず、各チャネルにおける相互作用ロジックのそれぞれの光検出器は、インコヒーレント検出によってそれぞれの光乗算器の出力を検出し、それによってそれぞれの読み取り値を生成するように構成される。
【0128】
[00134]本明細書に開示される別の態様によれば、関数を最適化する変数のベクトルの値を推定する方法が提供される。関数は、複数の項の重み付き和を含み、各項は、ベクトルからの変数の対応する部分集合の積を含み、各項は、変数間の相互作用をモデル化する重みの行列からの対応する重みによって重み付けされる。方法は、複数の並列ハードウェアチャネルのそれぞれにおいて、それぞれの信号生成器によって、それぞれの変数の値をモデル化する変調された特性を有するそれぞれのモデリング信号を生成するステップと、それぞれのモデリング信号のインスタンスを並列チャネルのそれぞれに供給するステップであって、各チャネルは、変数のベクトルをモデル化した信号のベクトルを受信する、ステップと、それぞれの相互作用ロジックにおいて、受信した信号のベクトルに、それぞれの変数と変数のベクトルとの間の相互作用をモデル化する重みの行列からの重みのそれぞれのベクトルを乗算し、それによって、それぞれのチャネルによってモデル化されたそれぞれの変数の寄与を表すそれぞれのフィードバック信号を生成するステップと、フィードバック信号をそれぞれの信号生成器に戻すステップと、それぞれの信号生成器によって、フィードバック信号に応じて、それぞれのモデリング信号を適応させるステップと、を含み、方法は、光学および／またはアナログ電子ハードウェアのみを使用して実施される。
【0129】
[00135]実施形態において、本方法は、本明細書に開示されたシステムの特徴のいずれかに従ったステップをさらに含み得る。
【図1】
【図2】
【図3】
【図4】
【図5A】
【図5B】
【図6】
【図6A】
【図6B】
【図7】
【図8】
【図9】
【国際調査報告】