WO2021125189A1

WO2021125189A1 - プログラム、計算機システム、及び計算機

Info

Publication number: WO2021125189A1
Application number: PCT/JP2020/046819
Authority: WO
Inventors: 純太渕脇; 裕也大西; 嘉夫滝本
Original assignee: Jsr株式会社
Priority date: 2019-12-19
Filing date: 2020-12-15
Publication date: 2021-06-24
Also published as: CN114846479A; JPWO2021125189A1; US20220310212A1

Abstract

本発明に係るプログラムは、与えられた処理を、量子アルゴリズムにより、または古典アルゴリズムと量子アルゴリズムとの併用アルゴリズムにより行う場合に、当該与えられた処理を実行する計算機を選択する処理を古典計算機に実行させるプログラムであって、算出ステップと、選択ステップと、制御ステップとを当該古典計算機に実行させる。算出ステップは、与えられた処理を量子アルゴリズムを用いて実行するのに必要な、または与えられた処理を前記併用アルゴリズムを用いて実行する場合における量子アルゴリズムの部分における、必要な量子ビット又は要求量子体積の大きさを算出する。選択ステップは、算出ステップにより算出された量子ビット又は要求量子体積に基づいて、与えられた処理を実行する計算機を選択する。制御ステップは、選択ステップにより選択された計算機が量子計算機を含む場合に、量子計算機に送信する、量子計算機が量子アルゴリズムを開始する指示に関する制御信号を生成する。

Description

プログラム、計算機システム、及び計算機

　本発明は、プログラム、計算機システム及び計算機に関する。

　古典計算機とは異なる計算パラダイムの下で計算を行う計算機として量子計算機がある。量子計算機は、例えば量子もつれ状態等の初期生成された量子状態に対して、ショアの素因数分解アルゴリズムやグローバーの探索アルゴリズム等の所定の量子アルゴリズムのもとでユニタリ行列を作用させ、生成された量子状態を基に計算結果を生成する計算機である。

　量子計算のシミュレーション及び実行環境は、例えば、Ａｚｕｒｅ（登録商標）　Ｑｕａｎｔｕｍ、ＡＷＳ（登録商標）　Ｂｒａｋｅｔ、ＩＢＭ　Ｑ－Ｅｘｐｅｒｉｅｎｃｅ（登録商標）などの、クラウドサービスの形で提供されている。これらのクラウドサービスでは、例えば、量子アルゴリズムの古典コンピュータ上でのテスト及びシミュレーション環境が与えられ、また作成した量子アルゴリズムを、ゲート方式の超伝導コンピュータ、量子アニーリング超伝導コンピュータ、イオントラップコンピュータ等の量子ハードウェア上で実行することができる。

　ここで、量子計算機と古典計算機の違いを例えば分子軌道計算を行う場合で説明すると、量子計算機は、古典計算機と比較して、原子数Ｎに対する計算の増大速度がゆるやかな一方で、大きな量子ビット数の量子計算機を作成するのが難しいため、量子ビット数による制限が課せられる。従って、与えられた処理の大きさによって、古典計算機と量子計算機のどちらで計算を行うかを使い分けるのが望ましい場合がある。

　また、すべての計算を量子計算機で行うのではなく、古典計算機と量子計算機の両方を用いて計算を行う併用アルゴリズムが知られている。例えば、ＶＱＥ（Variational　Quantum　Eigensolver）は、変分法をベースにして、量子計算機により所定の試行波動関数を生成した上でハミルトニアンの期待値を測定し、測定結果を基に古典計算機で処理を行うことで、変分法ベースで、古典計算機の中に量子計算機を組み込むアルゴリズムである。

　また、量子ビットの異なる複数の量子計算機を用いて計算を行う場合が考えられる。

Kandala,　A.,　Mezzacapo,　A.,　Temme,　K.　et　al.,　Nature,　２０１７年、５４９巻、ｐ２４２－２４６

　本発明は、古典計算機と量子計算機とが、または複数の量子計算機が混在する場合において、処理の最適化を行うことを目的とする。

　上述した課題を解決し、目的を達成するために、本発明に係るプログラムは、与えられた処理を、量子アルゴリズムにより、または古典アルゴリズムと量子アルゴリズムとの併用アルゴリズムにより行う場合に、当該与えられた処理を実行する計算機を選択する処理を古典計算機に実行させるプログラムであって、算出ステップと、選択ステップと、制御ステップとを当該古典計算機に実行させる。算出ステップは、与えられた処理を量子アルゴリズムを用いて実行するのに必要な、または与えられた処理を前記併用アルゴリズムを用いて実行する場合における量子アルゴリズムの部分における、必要な量子ビット又は要求量子体積の大きさを算出する。選択ステップは、算出ステップにより算出された量子ビット又は要求量子体積に基づいて、与えられた処理を実行する計算機を選択する。制御ステップは、選択ステップにより選択された計算機が量子計算機を含む場合に、量子計算機に送信する、量子計算機が量子アルゴリズムを開始する指示に関する制御信号を生成する。

　また、本発明に係る計算機システムは、与えられた処理を、量子アルゴリズムにより、または古典アルゴリズムと量子アルゴリズムとの併用アルゴリズムにより行う場合に、当該与えられた処理を実行する計算機を複数の計算機の中から選択する処理を実行する古典計算機と、複数の計算機とからなる。この複数の計算機は、少なくとも１つの量子計算機を含む。古典計算機は、算出部と、振り分け部と、送信部とを備える。算出部は、与えられた処理を量子アルゴリズムを用いて実行するのに必要な、または与えられた処理を併用アルゴリズムを用いて実行する場合における量子アルゴリズムの部分における、必要な量子ビット又は要求量子体積の大きさを算出する。振り分け部は、算出部が算出した量子ビット又は要求量子体積に基づいて、与えられた処理を実行する計算機を選択する。送信部と、振り分け部が選択した計算機が量子計算機を含む場合に、量子計算機が量子アルゴリズムを開始する指示を量子計算機に送信する。

　また、本発明に係る計算機は、与えられた処理を、量子アルゴリズムにより、または古典アルゴリズムと量子アルゴリズムとの併用アルゴリズムにより行う場合に、与えられた処理を実行する計算機を複数の計算機の中から選択する処理を実行する古典計算機である。算出部は、与えられた処理を量子アルゴリズムを用いて実行するのに必要な、または与えられた処理を併用アルゴリズムを用いて実行する場合における量子アルゴリズムの部分における、必要な量子ビット又は要求量子体積の大きさを算出する。振分け部は、算出部が算出した量子ビット又は要求量子体積に基づいて、与えられた処理を実行する計算機を選択する。送信部は、振り分け部が選択した計算機が量子計算機を含む場合に、量子計算機が量子アルゴリズムを開始する指示を量子計算機に送信する。

　本発明によれば、古典計算機と量子計算機とが、または複数の量子計算機が混在する場合において、処理の最適化を行うことができる。

図１は、本発明の実施形態に係る計算機システムの概略を表す図である。図２Ａは、実施形態に係る背景について説明する図である。図２Ｂは、実施形態に係る背景について説明する図である。図３は、第１の実施形態に係るプログラムが実行する処理の流れを示したフローチャートである。図４は、第１の実施形態の第１の変形例に係るプログラムが実行する処理の流れを示したフローチャートである。図５は、第２の実施形態に係るプログラムが実行する処理の流れを示したフローチャートである。図６は、図５におけるステップＳ２００の処理について詳細に説明したフローチャートである。

　以下、添付図面を参照しながら、本発明の実施形態に係るプログラム、計算機システム及び計算機を説明する。

（第１の実施形態）
　図１は、本発明の実施形態に係る計算機システム１の概略を表す図である。計算機システム１は、他の計算機の制御を行うための従来型の古典計算機１００と、量子計算を用いず従来型の方法で計算を行う古典計算機２００ａ、２００ｂ、２００ｃ…と、量子計算を用いて計算を行う量子計算機３００ａ、３００ｂ、３００ｃ…を含む。ここで、古典計算機１００及び２００ａ、２００ｂ、２００ｃは、例えばチューリングマシンとして表現できる通常の意味での計算機、例えばノイマン型の計算機である。

　なお、計算機システム１の構成としては、古典計算機２００ａ、２００ｂ、２００ｃ…と、量子計算機３００ａ、３００ｂ、３００ｃ…とのうち２以上の計算機を有していればよく、量子計算機のみを有する場合、古典計算機と量子計算機とのうち両方を有する場合の、いずれの場合であってもよい。

　古典計算機１００は、他の計算機（古典計算機２００ａ、２００ｂ、２００ｃ…、量子計算機３００ａ、３００ｂ、３００ｃ）の制御を行うための計算機であり、処理回路１０と、記憶部２０と、ディスプレイ３０と、入力装置４０とを備える。処理回路１０は、算出機能１０ａと、選択機能１０ｂと、制御機能１０ｃと、分割決定機能１０ｄと、分割機能１０ｅと、表示制御機能１０ｆと、生成機能１０ｇを備える。

　処理回路１０は、プログラムを記憶部２０から読み出して実行することで、算出機能１０ａ、選択機能１０ｂ等の機能を実現するプロセッサ（例えば、ＣＰＵ（Ｃｅｎｔｒａｌ　Ｐｒｏｃｅｓｓｉｎｇ　Ｕｎｉｔ）、ＧＰＵ（Ｇｒａｐｈｉｃａｌ　Ｐｒｏｃｅｓｓｉｎｇ　Ｕｎｉｔ）、ＡＳＩＣ（Ａｐｐｌｉｃａｔｉｏｎ　Ｓｐｅｃｉｆｉｃ　Ｉｎｔｅｇｒａｔｅｄ　Ｃｉｒｃｕｉｔ）、ＰＬＤ（Ｐｒｏｇｒａｍｍａｂｌｅ　Ｌｏｇｉｃ　Ｄｅｖｉｃｅ）等）である。すなわち、算出機能１０ａ、選択機能１０ｂ、制御機能１０ｃ、分割決定機能１０ｄ、分割機能１０ｅ、表示制御機能１０ｆ、生成機能１０ｇ等による処理機能は、当該プロセッサにより実行可能なプログラムの形態で記憶部２０に記憶される。これにより、プログラムが読み出された状態の処理回路１０は、算出機能１０ａ、選択機能１０ｂ、制御機能１０ｃ、分割決定機能１０ｄ、分割機能１０ｅ、表示制御機能１０ｆ、生成機能１０ｇ等の機能を有し、算出機能１０ａにより算出ステップを、選択機能１０ｂにより選択ステップを、制御機能１０ｃにより制御ステップを、分割決定機能１０ｄにより分割決定ステップを、分割機能１０ｅにより分割ステップを、表示制御機能１０ｆにより表示ステップを、生成機能１０ｇにより生成ステップを、古典計算機１００に実行させる。すなわち、算出機能１０ａ、選択機能１０ｂ、制御機能１０ｃ、分割決定機能１０ｄ、分割機能１０ｅ、表示制御機能１０ｆ、生成機能１０ｇは、それぞれ算出部、振り分け部(選択部)、制御部、分割決定部、分割部、表示制御部、生成部に対応する。なお、処理回路１０が有する各機能の処理の詳細については後述する。また、制御機能１０ｃは、種々の情報を送信する機能を備える送信部の一例ともなる。

　記憶部２０は、ＲＡＭ（Ｒａｎｄｏｍ　Ａｃｃｅｓｓ　Ｍｅｍｏｒｙ）、ＲＯＭ（Ｒｅａｄ　Ｏｎｌｙ　Ｍｅｍｏｒｙ）、フラッシュメモリ、ハードディスク等の記憶領域であり、処理回路１０が実行する各種プログラム及び当該各種プログラムの実行結果等を記憶する。

　ディスプレイ３０は、各種情報を表示する。例えば、ディスプレイ３０は、ＣＰＵの処理結果や、ユーザから各種操作を受け付けるためのＧＵＩ（Graphical　User　Interface）等を表示する。ディスプレイ３０は、例えば液晶ディスプレイ、有機ＥＬ（Electro　Luminescence）ディスプレイ又はブラウン管ディスプレイで構成される。なお、例えばタッチパネルのような形態で、入力装置４０とディスプレイ３０とが一体に構成されても良い。　

　入力装置４０は、ユーザが各種の操作を行うためのデバイスである。入力装置４０は、例えばマウス、キーボード、タッチパネル又はハードウェアキーで構成される。

　なお、実施形態は上述の例に限られない。実施形態によっては、記憶部２０、ディスプレイ３０、入力装置４０等は、古典計算機１００の必須の構成要素ではなく、すなわち古典計算機１００は記憶部２０、ディスプレイ３０、入力装置４０等を含まなくてもよい。また、別の例として、記憶部２０、ディスプレイ３０、入力装置４０等は、例えば、ネットワークを介して古典計算機１００の外に配置され、処理回路１０とデータのやりとりが行われてもよい。また、古典計算機１００は、２以上の古典計算機から構成されてもよい。例えば、古典計算機は、算出機能１０ａを有する第１の古典計算機と、第１の古典計算機とネットワークを通じて接続され、振り分け機能（選択機能１０ｂ）及び送信部としての制御機能１０ｃを有する第２の古典計算機とを備えても良い。

　古典計算機１００が、他の計算機の制御を行うための計算を行う計算機であるのに対して、古典計算機２００ａ、古典計算機２００ｂ、古典計算機２００ｃ…は、量子化学計算等、与えられた処理そのものを実行するための計算機である。例えば古典計算機２００ａは、古典計算機１００と同様に、処理回路１０と同様の処理回路１０ｘと、記憶部２０と同様の記憶部２０ｘとを備える。なお、処理回路１０の有するすべての機能を、処理回路１０ｘが有している必要はない。以下同様に、古典計算機２００ｂ、２００ｃ…も、処理回路１０と同様の処理回路と、記憶部２０と同様の記憶部を備える。

　量子計算機３００ａ、３００ｂ、３００ｃは、量子状態を用いて計算を行う量子計算機である。例えば、量子計算機３００ａは、初期量子状態を生成する量子情報生成部５０と、複数の量子ゲート６０と、量子ゲート制御部７０と、測定部８０とを備える。

　量子計算機の典型的な動作の例を説明すると、処理回路１０は、送信部としての制御機能１０ｃにより、初期状態の量子状態を生成するための制御信号を量子計算機３００ａに送信する。初期状態の量子状態を生成するための制御信号を古典計算機１００から取得した量子計算機３００ａは、量子情報生成部５０により、例えば量子エンタングルメント状態等の初期量子状態を生成する。

　量子ゲート６０は、生成された量子状態に作用し、量子状態を変化させる。量子ゲート６０の量子状態に対する作用は、ユニタリ行列の形で表すことができる。また、量子ゲート制御部７０は、量子ゲート６０の動作、すなわち当該ユニタリ行列を変化させる。測定部８０は、量子ゲート６０を通過した量子状態を測定する。これにより、量子計算機３００ａは、測定部８０が測定した測定結果を出力データとして得ることができる。

　また、計算機システム１は、ユーザからの情報の入力を受け付け、またはユーザに情報の表示を行うユーザ端末４００を備えてもよい。一例として、ユーザ端末４００は、ディスプレイ３０や入力装置４０と同様に、ユーザ端末側に、ディスプレイ３０ｚ及び入力装置４０ｚを備えても良い。

　続いて、実施形態に係る背景について簡単に説明する。

　古典計算機とは異なる計算パラダイムの下で計算を行う計算機として量子計算機がある。量子計算機は、例えば、初期生成された量子もつれ状態等の量子状態に対して、ショアの素因数分解アルゴリズムやグローバーの探索アルゴリズム等の所定の量子アルゴリズムのもとでユニタリ行列を作用させ、生成された量子状態を測定して計算結果を生成する計算機である。

　ここで、量子計算機と古典計算機の違いについて、例えば、分子軌道計算等の計算を例に取ると、図２Ａに示されるように、古典計算機では、原子数Ｎが大きくなると計算時間が急激に増大する。例えば、原子数Ｎに対して、計算時間がＮ^６に比例として増大する。一方で、図２Ｂに示されるように、量子計算機では、原子数Ｎに対して計算時間がそれほど大きく増大しない。例えば、原子数Ｎに対して、計算時間がＮに比例する時間で抑えられる。ただし、大きな量子ビットの量子計算機を作成するのは現状難しく、量子ビット数による制限がある。従って、与えられた処理の大きさによって、古典計算機と量子計算機のどちらで計算を行うかを使い分けるのが望ましい場合がある。

　また、大きな量子ビットの量子計算機を作成する困難さを考慮して、すべての計算を量子計算機による量子アルゴリズムで行うのではなく、古典計算機により実行される古典アルゴリズムと量子計算機により実行される量子アルゴリズムの両方を用いて計算を行う併用アルゴリズムが知られている。併用アルゴリズムとしては、例えばＶＱＥ（Variational　Quantum　Eigensolver）が挙げられる。ＶＱＥは、変分法をベースにして、量子計算機により所定の試行波動関数を生成した上でハミルトニアンの期待値を測定し、測定結果を基に古典計算機で処理を行うことで、変分法ベースで、古典計算機の中に量子計算機を組み込むアルゴリズムである。

　また、古典計算機と量子計算機の両方を用いた併用アルゴリズムを用いる場合の他に、与えられた処理を、古典計算機が担当する処理と、量子計算機が担当する処理とに分割し、各処理をそれぞれの計算機が実行することにより並列計算を行う場合が考えられる。また、与えられた処理を量子ビット数の異なる複数の量子計算機で分担して行う場合、与えられた処理をそれぞれの量子計算機が担当する処理に分割し、各処理をそれぞれの計算機が実行することにより並列計算を行う場合が考えられる。

　実施形態は、かかる背景に基づくものであって、少なくとも１つの量子計算機を含む複数の計算機で処理を行う場合における処理を最適化しようとするものである。

　図３は、第１の実施形態に係るプログラムが実行する処理の流れを示したフローチャートである。なお、以下の例としては、与えられた計算処理の例として、化合物の分子軌道計算を行う場合について説明する。

　はじめに、ステップＳ１００において、処理回路１０は、化合物の構造等の情報を記憶部より取得し、与えられた処理、例えば与えられた化学構造に対して量子化学計算を行う処理を行う準備を行う。すなわち、処理回路１０は、化学構造のリストを記憶する記憶部２０から読み出された化学構造を、当該与えられた化学構造とし、量子化学計算を行う処理を行う準備を行う。

　続いて、ステップＳ１１０において、処理回路１０は、算出機能１０ａにより、必要な量子ビット(qubit)又は要求量子体積の大きさを算出する。すなわち、実施形態に係るプログラムは、ステップＳ１１０において、与えられた処理、例えば与えられた化学構造に対して量子化学計算を行う処理を、量子アルゴリズムを用いて実行するのに必要な、または与えられた処理を併用アルゴリズムを用いて実行する場合における量子アルゴリズムの部分における、必要な量子ビット又は要求量子体積の大きさを算出する算出ステップを古典計算機１００に実行させる。ここで、要求量子体積は、要求される量子体積(quantum　volume)の大きさ、すなわち、行うことのできる量子計算の規模、すなわち、空間的サイズ及び時間を表す。量子体積は、単純な量子ビット数に、発生するエラーの大きさを加味して算出される指標である。

　なお、ステップＳ１１０において処理回路１０が算出する必要な量子ビット又は要求量子体積の算出に当たっては、処理回路１０は、算出機能１０ａにより、若干の余裕度を持たせて必要な量子ビット又は要求量子体積の値を算出してもよく、また算出するのが望ましい。ここで、当該余裕度は、計算の種類や要求される計算速度、計算機の能力、キューの混雑状況などに応じて任意に設定可能である。

　必要な量子ビット又は要求量子体積の大きさを算出する方法として、例えば、処理回路１０は、算出機能１０ａにより、原子数と、当該原子数の分子軌道計算を行うのに必要な量子ビットまたは要求量子体積との関係を表すテーブルと、原子数とに基づいて、必要な量子ビット又は要求量子体積の大きさを算出する。当該テーブルは、例えば、計算量を正確に見積もることができる程度に原子数が小さい場合における原子数と必要な量子ビットまたは要求体積との関係式を、原子数が大きい場合に外挿することにより作成される。

　続いて、ステップＳ１２０において、与えられた処理を、量子アルゴリズムによりまたは古典アルゴリズムと量子アルゴリズムとの併用アルゴリズムにより、または古典アルゴリズムと量子アルゴリズムの両者を併用して処理を行う場合に、第１の実施形態に係るプログラムは、与えられた処理を実行する計算機を選択する選択ステップを、選択機能１０ｂにより古典計算機１００に実行させる。具体的には、処理回路１０は、選択機能１０ｂにより、ステップＳ１１０において算出された量子ビットまたは要求量子体積に基づいて、与えられた処理を実行する計算機を選択する。

　ここで、処理を行う候補となる計算機が複数の量子計算機である場合には、与えられた処理を実行する計算機を選択する選択ステップは、複数の量子計算機の中から、ステップＳ１１０において算出された量子ビットまたは要求量子体積に基づいて、与えられた処理を実行する量子計算機を選択するステップとなる。

　また、処理を行う候補となる計算機が１以上の量子計算機と１以上の古典計算機である場合には、与えられた処理を実行する計算機を選択する選択ステップは、古典計算機と量子計算機との両方の計算機の中から、ステップＳ１１０において算出された量子ビットまたは要求量子体積に基づいて、古典計算機と量子計算機とのうち、与えられた処理を実行する計算機を選択するステップとなる。この場合、一例として、処理回路１０は、選択機能１０ｂにより、量子ビットまたは要求量子体積が、候補となる量子計算機により実現できる量子ビットまたは要求量子体積の範囲内であれば、量子計算機を、与えられた処理を実行する計算機として選択し、量子ビットまたは要求量子体積が、候補となる量子計算機により実現できる量子ビットまたは要求量子体積の範囲を超えている場合、古典計算機を、与えられた処理を実行する計算機として選択する。

　例えば、量子計算機３００ａの量子ビット数が１００量子ビットであり、ステップＳ１１０において算出された量子ビット数が８０量子ビットである場合、ステップＳ１２０において、処理回路１０は、選択機能１０ｂにより、量子計算機３００ａを、与えられた処理を実行する計算機として選択する。量子計算機３００ａは、量子アルゴリズムにより、当該与えられた処理を実行する。一方、量子計算機３００ａの量子ビット数が１００量子ビットであり、ステップＳ１１０において算出された量子ビット数が１２０量子ビットである場合、ステップＳ１２０において、処理回路１０は、選択機能１０ｂにより、古典計算機２００ａを、与えられた処理を実行する計算機として選択する。与えられた処理が分子軌道計算の場合、古典計算機２００ａは、例えばＨａｒｔｒｅｅ－Ｆｏｃｋ法、ＤＦＴ法、ＣＣＳＤ（Ｔ）法等により、分子軌道計算を実行する。

　また、古典アルゴリズムと量子アルゴリズムとの併用アルゴリズムを用いる場合、与えられた処理を実行する計算機を選択する選択ステップは、ステップＳ１１０において算出された量子ビットまたは要求量子体積に基づいて、与えられた処理を実行する量子計算機を選択し、また、量子計算機の選択とは別途、古典アルゴリズムを実行する古典計算機を選択するステップとなる。

　これらの選択ステップにおいて、ステップＳ１１０で算出された量子ビットが、選択される計算機で利用可能な量子ビットを上回らないような、またはステップＳ１１０で算出された要求量子体積が、選択される計算機で実現可能な要求量子体積を上回らないように、処理回路１０は、選択機能１０ｂにより、量子計算機の選択を行う。

　なお、かかる選択ステップにおいて、システム中に量子ビット又は量子体積の異なる複数の量子計算機が存在する場合、選択ステップはその中から対象化合物の計算に必要な量子ビット又は量子体積を下回らず、かつ最も近い量子ビット又は体積を有する量子計算機を優先的に選択することが好ましい。

　なお、選択ステップにおいて、処理回路１０は、選択機能１０ｂにより、算出された量子ビットまたは要求量子体積の他に、利用可能な計算機リソースの能力、特性、空き具合などをも加味して、計算機の選択を行ってもよい。一例として、処理回路１０は、選択機能１０ｂにより、計算機リソースの能力が高く、特性が与えられた処理に適した特性であり、または比較的空いている計算機を処理を実行する計算機として選択する。

　また、選択ステップにおいて、処理回路１０は、選択機能１０ｂにより、計算機の選択のみならず、計算機を実行する実行タイミングや実行順序をあわせて選択してもよい。

　また、ステップＳ１２０の処理の実行とあわせて、処理回路１０は、表示制御機能１０ｆにより、算出機能１０ａにより算出された量子ビット若しくは要求量子体積の大きさ、選択機能１０ｂにより選択された計算機に係る情報及び、構造式又は原子座標又は化学構造を表す文字列または化合物を特定する情報を表示部に適宜表示する。ここで、表示部は、ディスプレイ３０であってもよいし、古典計算機１００が配置される端末であるプロバイダ側端末とは異なる端末であって、ユーザが操作可能な端末であるユーザ端末４００であってもよい。

　また、処理回路１０は、制御機能１０ｃにより、算出機能１０ａにより算出された量子ビット若しくは量子体積の大きさ、又は選択機能１０ｂにより選択された計算機に係る情報を、他の計算機やユーザ端末４００、計算機システム１の外部の端末等に送信する機能を備えても良い。

　続いて、ステップＳ１３０において、実施形態に係るプログラムは、ステップＳ１２０において選択された計算機が量子計算機を含む場合に、量子計算機に送信する制御信号を生成する制御ステップを、送信部としての選択機能１０ｂにより、古典計算機１００に実行させる。ここで、当該制御信号は、量子計算機が量子アルゴリズムを開始する指示に関する制御信号であり、例えば、量子計算機が量子アルゴリズムを開始時に生成する量子状態を生成するのに必要な情報である。例えば、量子アルゴリズムを開始時に生成する量子状態が、α_１|０＞+α_２｜１＞で表せるとすると、係数α_１、α_２が、量子計算機が量子アルゴリズムを開始時に生成する量子状態α_１|０＞+α_２｜１＞を生成するのに必要な情報となる。続いて、処理回路１０は、送信部としての制御機能１０ｃにより、量子計算機に、量子計算機が量子アルゴリズムの開始時に生成する量子状態を生成するのに必要な情報を制御信号として送信する。

　続いて、ステップＳ１４０において、古典計算機１００から取得した制御信号に基づいて、量子計算機は、量子計算を実行する。また、古典計算機により計算を実行する場合、古典計算機１００から取得した制御信号に基づいて、古典計算機は、通常の計算を実行する。また、古典計算機による計算と量子計算機による計算の併用アルゴリズムにより計算を実行する場合、古典計算機１００から取得した制御信号に基づいて、古典計算機と量子計算機とが、計算を行う。古典計算機１００は、送信部としての制御機能１０ｃにより、複数の計算機に対して制御信号を送信し、複数の計算機のうちそれぞれが、それぞれが受信した制御信号を基に、それぞれの処理を行う。

　以下の実施形態では、古典計算機２００ａ、２００ｂ、２００ｃ…及び量子計算機３００ａ、３００ｂ、３００ｃ…が、古典アルゴリズムと量子アルゴリズムの併用アルゴリズムであるＶＱＥ（Variational　Quantum　Eigensolver）法を用いて分子軌道計算を行う場合について説明する。

　ＶＱＥ法は、変分パラメータをθとして、系のハミルトニアンＨに対する試行波動関数ψ（θ）に係る期待値＜ψ（θ）｜Ｈ｜ψ（θ）＞が常に基底状態のエネルギーＥを上回るという変分原理を利用して、分子軌道計算を行うアルゴリズムの中に、量子計算機による計算を組み込んだ方法である。

　ここで、量子計算機３００ａ等は、ある変分パラメータθに対して、変分法の対象となる量子状態ψ（θ）を、量子アルゴリズムを用いて量子状態として生成する。また、量子計算機３００ａ等は、量子情報生成部５０により、分子軌道計算等を行う系のハミルトニアンＨを、量子計算機上に生成される量子状態ψ（θ）に対するハミルトニアンとして量子計算機上で表現する。しかる後、測定部８０は、量子計算機上に生成される量子状態ψ（θ）に対する当該ハミルトニアンの期待値＜ψ（θ）｜Ｈ｜ψ（θ）＞を測定する。これにより、量子計算機３００ａは、与えられた変分パラメータθのもとでの、系のハミルトニアンＨの、生成された量子状態ψ（θ）における期待値を、測定値として得ることができる。

　一方、古典計算機２００ａ等は、その他の部分の処理、例えば、変分パラメータθを選択する処理、ハミルトニアンの期待値が最小となる変分パラメータθを算出する処理、変分パラメータθの値に基づいて、どのような量子状態ψ（θ）を選択するかを決定する処理などの処理を担当し、変分パラメータθの最適化を行う。これにより、古典計算機２００ａ等と量子計算機３００ａ等が共同して、量子化学計算等を行うことができる。

　以上のように、第１の実施形態に係るプログラムによれば、古典計算機と量子計算機とが、または複数の量子計算機が混在する場合において、処理を行う計算機を適切に選択することができ、その結果、処理の最適化を行うことができる。

　（第１の実施形態の第１の変形例）

　第１の実施形態では、算出ステップにおいて算出された量子ビット又は要求量子体積の大きさに基づいて、与えられた処理を実行する計算機をプログラムが自動的に選択する場合について説明した。実施形態はこれに限られず、第１の実施形態の第１の変形例においては、プログラムが、算出した量子ビット又は要求量子体積の大きさをディスプレイ３０に表示させ、ディスプレイ３０に表示された量子ビットまたは要求量子体積の大きさを見たユーザからのフィードバックを受け付けてもよい。

　かかる状況が、図４に示されている。なお、ステップＳ１００、ステップＳ１１０、ステップＳ１３０及びステップＳ１４０については、図２ですでに説明した処理と同様の処理であるので、重ねての説明は省略する。

　ステップＳ１２１において、第１の実施形態の第１の変形例に係るプログラムは、処理回路１０の表示制御機能１０ｆにより、算出ステップＳ１１０により算出された量子ビット又は要求量子体積の大きさを表示部に表示させる表示制御ステップを、古典計算機１００に更に実行させる。

　続いて、ステップＳ１２２において、第１の実施形態の第１の変形例に係るプログラムは、古典計算機１００に、処理回路１０により、ステップＳ１２１において算出された量子ビット又は要求量子体積に基づいて、処理を実行可能な計算機のリストを生成する生成ステップを実行させる。

　続いて、ステップＳ１２３において、当該プログラムは、古典計算機１００に、処理回路１０の表示制御機能１０ｆにより、ステップＳ１２２において生成された処理を実行可能な計算機のリストを表示部に表示させる表示制御ステップを実行させる。

　続いて、ステップＳ１２４において、当該プログラムは、古典計算機１００に、処理回路１０の制御機能１０ｃにより、化合物名等の入力データを修正するか否かの入力をユーザから受け付ける。

　処理回路１０は、制御機能１０ｃにより、入力データを修正する旨の入力をユーザから受け付けた場合、（ステップＳ１２５　Ｙｅｓ）、処理はステップＳ１００に戻る。一方、入力データを修正しない場合（ステップＳ１２５　Ｎｏ）、処理はステップＳ１２３に移る。

　続いて、ステップＳ１２６において、第１の実施形態の第１の変形例に係るプログラムは、古典計算機１００に、処理回路１０の制御機能１０ｃにより、与えられた処理を実行する計算機の入力を手動で受け付けるステップを実行させる。

　続いて、ステップＳ１２７において、第１の実施形態の第１の変形例に係るプログラムは、古典計算機１００に、処理回路１０の選択機能１０ｂにより、ステップＳ１２６でユーザから受け付けた入力に基づいて、与えられた処理を実行する計算機を選択させるステップを実行させる。

　また、第１の実施形態と同様に、第１の実施形態の第１の変形例に係るプログラムは、古典計算機１００に、ステップＳ１３０及びステップＳ１４０と同様の処理を実行させる。

　以上のように、第１の実施形態に係るプログラムによれば、古典計算機と量子計算機とが、または複数の量子計算機が混在する場合において、必要に応じて、ユーザからの入力を受け付けることで、処理を行う計算機を適切に選択することができ、その結果、処理の最適化を行うことができる。

　（第２の実施形態）
　第１の実施形態では、与えられた処理を、複数の小処理に分割せず、一つの計算機においてすべての処理を行う場合について説明した。第２の実施形態では、与えられた処理を、複数の小処理に分割し、複数の小処理のそれぞれごとに選択された計算機で、古典計算機又は量子計算機が処理を行う場合について説明する。

　図５及び図６に、かかる状況が示されている。図５は、第２の実施形態に係るプログラムが古典計算機１００に実行させる処理の流れを示したフローチャートであり、図６は、図５におけるステップＳ２００の処理について詳細に説明したフローチャートである。

　まずはじめに、第２の実施形態に係るプログラムは、ステップＳ２００において、古典計算機１００に、処理回路１０の分割決定機能１５ｄにより、与えられた処理を、複数の小処理に分割する分割の仕方を決定する分割決定ステップを実行させる。

　図６を用いて、ステップＳ２００の分割決定ステップについて詳細に説明する。図６全体の処理が、図５のステップＳ２００に対応する。

　はじめに、第２の実施形態に係るプログラムは、ステップＳ２０１において、古典計算機１００に、各計算機の計算性能を取得させる。

　続いて、第２の実施形態に係るプログラムは、ステップＳ２０２において、古典計算機１００に、処理回路１０の分割決定機能１０ｄにより、ステップＳ２０１で取得した計算性能に基づいて、与えられた処理を複数の小処理に仮分割させるステップを実行させる。

　続いて、第２の実施形態に係るプログラムは、ステップＳ２０３において、古典計算機１００に、仮分割された複数の小処理それぞれについて、量子ビットまたは要求量子体積の大きさを算出させるステップを実行させる。

　続いて、仮分割が、所定の条件を満たしており適切な分割である場合（ステップＳ２０４　Ｙｅｓ）、処理はステップＳ２０６に進み、処理回路１０は、当該仮分割が、適切な分割方法であると決定し、処理回路１０は、当該分割を、適切な分割方法であると決定する。

　一方、仮分割が、適切な分割でない場合（ステップＳ２０４　Ｎｏ）、第２の実施形態に係るプログラムは、ステップＳ２０５において、古典計算機１００に、新たな分割方法で仮分割を行われるステップを実行させ、その後処理はステップＳ２０２へと戻る。

　このようなステップＳ２０１～ステップＳ２０６により、第２の実施形態に係るプログラムは、上述したステップＳ２００に係る処理を、古典計算機１００に実行させる。

　続いて、図５に戻り、第２の実施形態に係るプログラムは、ステップＳ２１０において、古典計算機１００に、処理回路１０の分割機能１０ｅにより、ステップＳ２００において決定された分割の仕方に基づいて、与えられた処理を、複数の小処理に分割するステップを実行させる。

　なお、実施形態は、ステップＳ２００及びステップＳ２１０の処理を経るまでもなく、与えられた処理が、複数の小処理からなるものである場合であってもよく、かかる場合は、ステップＳ２００及びステップＳ２１０の操作は必要に応じて省略される。

　続いて、第２の実施形態に係るプログラムは、ステップＳ２２０において、古典計算機１００に、処理回路１０の算出機能１０ａにより、複数の小処理それぞれについて、量子ビット又は要求量子体積の大きさを算出させる算出ステップを実行させる。

　続いて、第２の実施形態に係るプログラムは、ステップＳ２３０において、古典計算機１００に、複数の小処理それぞれについて、算出された量子ビット又は要求量子体積と、利用可能な計算機リソースの能力、特性、空き具合等に応じて、計算を行う計算機を選択させる選択ステップを実行させる。例えば、第２の実施形態に係るプログラムは、古典計算機１００に、計算能力の高い計算機、与えられた処理の性質と、計算機の特性がマッチする計算機、比較的空いている計算機を、計算を行う計算機として優先的に選択させる。

　続いて、第２の実施形態に係るプログラムは、ステップＳ２４０において、古典計算機１００に、処理回路１０の選択機能１０ｂにより、更に、複数の小処理それぞれについて、算出された量子ビット又は要求量子体積と、利用可能な計算機リソースの能力、特性、空き具合等に応じて、計算機を実行する実行タイミングまたは実行順序を選択させる。一例として、第２の実施形態に係るプログラムは、計算機の処理負荷が特定の計算機に集中しないように、計算機を実行する実行タイミングまたは実行順序を選択させる。

　続いて、第２の実施形態に係るプログラムは、ステップＳ２５０において、古典計算機１００に、処理回路１０の制御機能１０ｃにより、量子計算機が量子アルゴリズムの開始時に生成する量子状態を生成するのに必要な情報を制御信号として生成させる制御ステップを実行させる。

　処理回路１０は、送信部としての制御機能１０ｃにより、生成された制御信号を、ステップＳ２３０で選択した計算機に送信させる。

　続いて、ステップＳ２６０において、量子計算機が量子計算を実行する。

　以上のように、第２の実施形態に係るプログラムによれば、古典計算機と量子計算機とが、または複数の量子計算機が混在する場合において、与えられた処理が複数の処理からなる場合において、適切に処理が分割され、分割された処理ごとに処理を行う計算機が適切に選択されることで、処理の最適化を行うことができる。

　このように、本発明の実施形態によれば、古典計算機と量子計算機とが、または複数の量子計算機が混在する場合において、処理を最適化することができる。

　なお、本発明の実施形態は、例として提示されたものにすぎず、種々の置き換え、変更等を行うことができる。

　１０　処理回路
　１０ａ　算出機能
　１０ｂ　選択機能
　１０ｃ　制御機能
　１０ｄ　分割決定機能
　１０ｅ　分割機能
　１０ｆ　表示制御機能
　１０ｇ　生成機能
　２０　記憶部
　３０　ディスプレイ
　４０　入力装置
　２００ａ　古典計算機
　３００ａ　量子計算機
　４００　ユーザ端末

Claims

　与えられた処理を、量子アルゴリズムにより、または古典アルゴリズムと量子アルゴリズムとの併用アルゴリズムにより行う場合に、前記与えられた処理を実行する計算機を選択する処理を古典計算機に実行させるプログラムであって、
　前記与えられた処理を前記量子アルゴリズムを用いて実行するのに必要な、または前記与えられた処理を前記併用アルゴリズムを用いて実行する場合における前記量子アルゴリズムの部分における、必要な量子ビット又は要求量子体積の大きさを算出する算出ステップと、
　前記算出ステップにより算出された前記量子ビット又は前記要求量子体積に基づいて、前記与えられた処理を実行する前記計算機を選択する選択ステップと、
　前記選択ステップにより選択された前記計算機が量子計算機を含む場合に、前記量子計算機に送信する、前記量子計算機が前記量子アルゴリズムを開始する指示に関する制御信号を生成する制御ステップ
とを前記古典計算機に実行させるプログラム。
　前記制御ステップは、前記選択ステップにより選択された前記計算機が量子計算機を含む場合に、前記量子計算機が前記量子アルゴリズムの開始時に生成する量子状態を生成するのに必要な情報を前記制御信号として生成する、請求項１に記載のプログラム。
　前記与えられた処理は、与えられた化学構造に対して量子化学計算を行う処理である、請求項１に記載のプログラム。
　前記与えられた化学構造は、化学構造のリストを記憶する記憶部から読み出された化学構造である、請求項３に記載のプログラム。
　前記算出ステップは、前記与えられた処理が複数の小処理からなるものである場合に、前記複数の小処理のうちそれぞれについて前記量子ビット又は前記要求量子体積の大きさを算出し、
　前記選択ステップは、前記与えられた処理が複数の小処理からなるものである場合に、
前記複数の小処理のうちそれぞれについて、前記計算機を選択し、
　前記制御ステップは、前記与えられた処理が複数の小処理からなるものである場合に、
前記複数の小処理のうちそれぞれについて、前記制御信号を生成する、請求項１に記載のプログラム。
　前記与えられた処理を、複数の小処理に分割する分割のしかたを決定する分割決定ステップと、
　決定された前記分割のしかたに基づいて、前記与えられた処理を、前記複数の小処理に分割する分割ステップとを、前記古典計算機に実行させ、
　前記選択ステップは、前記分割ステップにより分割された前記複数の小処理のうちそれぞれについて、前記計算機を選択し、
　前記制御ステップは、前記分割ステップにより分割された前記複数の小処理のうちそれぞれについて、前記制御信号を生成する、請求項１に記載のプログラム。
　前記選択ステップは、前記算出ステップにより算出された前記量子ビットまたは前記要求量子体積と、利用可能な計算機リソースの能力、特性及び空き具合のうち少なくとも一つとに応じて、前記計算機を選択する、請求項５又は６に記載のプログラム。
前記選択ステップは、前記算出ステップにより算出された前記量子ビットまたは前記要求量子体積と、利用可能な計算機リソースの能力、特性及び空き具合のうち少なくとも一つとに応じて、更に前記計算機を実行する実行タイミングまたは実行順序を選択する、請求項７に記載のプログラム。
　前記選択ステップは、前記量子ビットが、前記計算機で利用可能な量子ビットを上回らないように、または前記要求量子体積が、前記計算機の要求量子体積を上回らないように、前記計算機を選択する、請求項１～８のいずれか一つに記載のプログラム。
　前記算出ステップにより算出された前記量子ビット又は前記要求量子体積の大きさを表示部に表示させる表示制御ステップを更に前記古典計算機に実行させる、請求項１～９のいずれか一つに記載のプログラム。
　前記算出ステップにより算出された前記量子ビット又は前記要求量子体積に基づいて、前記処理を実行可能な計算機のリストを生成する生成ステップと、
　前記生成ステップにより生成された前記処理を実行可能な計算機のリストを表示部に表示させる表示制御ステップとを更に前記古典計算機に実行させる、請求項１～９のいずれか一つに記載のプログラム。
　前記表示部は、前記古典計算機が配置される端末であるプロバイダ側端末とは異なる端末であって、ユーザが操作可能な端末であるユーザ端末に配置される、請求項１０又は１１に記載のプログラム。
　前記選択ステップは、古典計算機と量子計算機との両方の計算機の中から、前記量子ビット又は前記要求量子体積に基づいて、前記処理を実行する計算機を選択する、請求項１～１２のいずれか一つに記載のプログラム。
　前記併用アルゴリズムは、ＶＱＥ（variational　quantum　eigensolver）法である、請求項１～１３のいずれか一つに記載のプログラム。
　与えられた処理を、量子アルゴリズムにより、または古典アルゴリズムと量子アルゴリズムとの併用アルゴリズムにより行う場合に、前記与えられた処理を実行する計算機を複数の計算機の中から選択する処理を実行する古典計算機と、前記複数の計算機とからなる計算機システムであって、
　前記複数の計算機は、少なくとも１つの量子計算機を含み、
　前記古典計算機は、
　前記与えられた処理を前記量子アルゴリズムを用いて実行するのに必要な、または前記与えられた処理を前記併用アルゴリズムを用いて実行する場合における前記量子アルゴリズムの部分における、必要な量子ビット又は要求量子体積の大きさを算出する算出部と、
　前記算出部が算出した前記量子ビット又は前記要求量子体積に基づいて、前記与えられた処理を実行する前記計算機を選択する振り分け部と、
　前記振り分け部が選択した前記計算機が前記量子計算機を含む場合に、前記量子計算機が前記量子アルゴリズムを開始する指示を前記量子計算機に送信する送信部と
を備える、計算機システム。
　前記算出部が算出した前記量子ビット若しくは前記要求量子体積の大きさ、又は前記振り分け部が選択した前記計算機に係る情報を送信する機能を備える請求項１５に記載の計算機システム。
　前記算出部が算出した前記量子ビット若しくは前記要求量子体積の大きさ、又は前記振り分け部が選択した前記計算機に係る情報をユーザに表示する表示部を更に備える、請求項１５に記載の計算機システム。
　前記算出部が算出した前記量子ビット若しくは前記要求量子体積の大きさ、又は前記振り分け部が選択した前記計算機に係る情報及び、構造式又は原子座標又は化学構造を表す文字列または化合物を特定する情報をユーザに表示する表示部を更に備える、請求項１７に記載の計算機システム。
　前記古典計算機は、前記算出部の機能を有する第１の古典計算機と、前記第１の古典計算機とネットワークを通じて接続され、前記振り分け部及び前記送信部の機能を有する第２の古典計算機とを備える、請求項１５に記載の計算機システム。
　与えられた処理を、量子アルゴリズムにより、または古典アルゴリズムと量子アルゴリズムとの併用アルゴリズムにより行う場合に、前記与えられた処理を実行する計算機を複数の計算機の中から選択する処理を実行する古典計算機であって、
　前記与えられた処理を前記量子アルゴリズムを用いて実行するのに必要な、または前記与えられた処理を前記併用アルゴリズムを用いて実行する場合における前記量子アルゴリズムの部分における、必要な量子ビット又は要求量子体積の大きさを算出する算出部と、
　前記算出部が算出した前記量子ビット又は前記要求量子体積に基づいて、前記与えられた処理を実行する前記計算機を選択する振り分け部と、
　前記振り分け部が選択した前記計算機が量子計算機を含む場合に、前記量子計算機が前記量子アルゴリズムを開始する指示を前記量子計算機に送信する送信部と
を備える、計算機。