2016/04/18

～光を使って難問を解くコンピュータの実現に道～

⇒ニュースリリース
⇒量子光制御研究グループ

研究の背景と経緯

インターネットや交通網、ソーシャルネットワークなど、社会を構成する様々なシステムが大規模化、複雑化する現在、それらのシステムをいかに効率よく運用するかは重要な課題となっています。これらの課題の多くは、組合せ最適化問題と呼ばれる数学的問題に帰着します。組合せ最適化問題とは、数多くの選択肢の組合せの中から最も良いものを見つけ出す問題で、選択肢が多くなると計算時間が爆発的に増大するため、現代のコンピュータでは解くことが大変困難であることが知られています。

一方、組合せ最適化問題の多くは相互作用するスピン群のモデルである「イジングモデル」のエネルギー最小状態(基底状態)を求める問題に帰着可能です。最近、人工的に作製したスピンを用いてイジングモデルを模擬し、そのエネルギー最小状態を求めることで複雑な組合せ最適化問題を高速に解く試みが多くの研究機関で盛んになってきました。

中でも、コヒーレントイジングマシン(coherent Ising machine)は、光を用いて実現した人工スピン群により高速にイジングモデルを解く可能性がある計算機として現在注目を集めています(図１、図２)。本方式では、光パラメトリック発振器(optical parametric oscillator: OPO)をスピンとして用います。OPOは、0またはπの位相しかとらない特殊なレーザ発振器で、位相0, πをそれぞれ上向き、下向きのスピンに対応させることができます。各OPOから出力される光を、光伝送路を介して互いに注入することで、スピン間の相互作用を実現します。光伝送路によりネットワーク化されたOPO群は、多くの場合ネットワーク全体の損失が最小となる位相の組合せで発振するため、イジングモデルの基底状態を与えるスピンの組合せを高い確率で得ることができます。

2014年にスタンフォード大学のグループが4つのOPOを人工スピンとして用いてコヒーレントイジングマシンの原理確認実験に成功しています。しかし、現実社会で課題となっている複雑な組合せ最適化問題に適用するためには、スピン数を飛躍的に増大する必要がありました。

研究の内容

今回、長さ1 kmという長距離光ファイバ共振器を用いて、時間多重された10,000個を超えるOPOを一括発生することに成功しました(図３)。光ファイバ通信の研究開発で培った技術を用いることで、現実社会で課題となっている複雑な問題を解くことが可能な大規模コヒーレントイジングマシンを実現するための多数の人工スピンの生成が可能となりました。 　また、隣接するOPO間に光結合を導入することで、最もシンプルなネットワークである1次元のイジングモデルを模擬する実験を行いました。この実験により、室温で動作するOPO群が、低温下のスピンの振る舞いをよく模擬する、高品質な人工スピンとして動作することを確認しました。

大規模時分割多重OPOの発生

1. 長さ1 kmの高非線形光ファイバ中の四光波混合より、位相0とπの光だけを増幅する位相感応増幅器を実現
2. 位相感応増幅器を光ファイバループに挿入することで、長さ1 kmの長距離光ファイバ共振器を構築し、位相感応増幅器を2 GHzの繰り返しで動作することにより、10,320個のOPOを発生(図３)
3. 発生したOPOの位相が0またはπに二値化しており、安定な人工スピンとして使用可能であることを実験で確認(図４)

1次元スピンネットワークの模擬

1. 1ビット遅延干渉計を共振器中に挿入することで、各OPOの光の一部を隣接するOPOに注入し、1次元の結合を実現(図５)
2. OPO結合の位相を調整することで、ネットワーク化された人工スピンが強磁性的・反強磁性的に振る舞うことを観測(図６)。また、位相感応増幅器の駆動条件を変えることで、ネットワークのエネルギーが低下することを観測(図７)。この結果、室温で動作するOPO群が、高品質な人工スピンとして使用可能であることを確認。

技術のポイント

1. 大規模時分割多重OPOの発生
   実験系を図3に示します。パルス幅100 ps、繰り返し周波数 2 GHz (パルス間隔500 ps)、波長1531 nmおよび1551 nmの2波長のポンプパルス列を、波長多重フィルタを介して光共振器に入力します。光ファイバ共振器には、長さ1 kmの高非線形光ファイバ、波長フィルタ、偏波コントローラ、および99:1光カプラが含まれています。高非線形光ファイバ中の2波長四光波混合により、2つのポンプ波長の中間点にあたる1541 nmにおいてポンプ光の位相に対して0またはπの位相成分だけが増幅される位相感応増幅と呼ばれる現象が得られます。ポンプ光入力直後に位相感応増幅により発生した雑音が種光となり、波長フィルタを透過する1541 nmの光のみが光ファイバ共振器中を周回するようにすることで、位相が0またはπのみの光が発振するOPOを実現できます。光の共振器1周時間が約5 μs、ポンプパルスの時間間隔が500 psなので、1つの共振器だけを用いて10,000を超えるOPOを発生することができます。発生したOPOの一部は99:1カプラにより共振器外に取り出し、1ビット遅延干渉計を用いて隣接パルス間の位相差を測定します。位相差測定の結果を図4(a)に示します。このように、位相差のコサイン成分が1または-1に二極化、すなわち各OPOの位相が0またはπのみをとっています。発生した全OPOに相当する10,320パルスの位相差測定結果のヒストグラムを図4(b)に示します。このように、全てのOPOの位相は0またはπに離散化しており、1万個以上の安定な人工スピンが発生したことが確認できました。
2. 1次元スピンネットワークの模擬実験
   生成したOPO群を用いて最もシンプルなネットワークである1次元イジングモデルを模擬する実験を行いました。図３に示す実験系において、位相差測定に用いたものと同様の1ビット遅延干渉計を共振器中に挿入することで、各OPOの光の一部を隣接するOPOに注入し、1次元の結合を実現できます(図５)。干渉計の2光路の位相差を0またはπとすることで、強磁性結合(スピンの向きが同じとなる)および反強磁性結合(逆となる)を実現可能です。位相差測定の結果を図6(a)(干渉計位相差0)、(b)(π)に示します。このように、干渉計の位相差に応じて、スピン群が強磁性的／反強磁性的な振る舞いをすることが確認できました。また、図6(a)中に示すように、強磁性結合においても全てのスピンは揃わずドメイン構造を持つことが観測されました。これは、理論的に示されている「有限温度下では基底状態にならない」という1次元イジングモデルの特性をよく再現しています。さらに、強磁性結合においてポンプ光振幅を変えた場合にドメイン長の分布がどのように変化するかを調べた結果を図７に示します。ポンプ光振幅が発振閾値に近いほど、長いドメインが形成され、基底状態(ドメインが無くなり、全てのスピンが揃う状態)に近づくことがわかります。

今後の展開

今回の実験では隣接OPO間の結合のみを用いたため、模擬したネットワーク構造は1次元にとどまっていました。今後は、任意のOPO間に結合を実装することで、より複雑なスピンネットワークを実現し、大規模な組合せ最適化問題の解探索が可能なコヒーレントイジングマシンの実現を目指します。

論文情報

"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

用語解説

※１ ... スピン
電子などの素粒子は小さな磁石としての性質をもつ。この磁石としての性質を古典的な球の自転になぞらえて「スピン」と呼ぶ。

※２ ... 四光波混合
入力した光に対して非線形な応答を示す媒体に異なる波長をもつ2つまたは3つの光を入力したとき、新たな波長の光が発生する現象。