JP7426012B2 - エネルギーの微分を求めるための量子情報処理方法、古典コンピュータ、量子コンピュータ、量子情報処理プログラム、及びデータ構造 - Google Patents

エネルギーの微分を求めるための量子情報処理方法、古典コンピュータ、量子コンピュータ、量子情報処理プログラム、及びデータ構造 Download PDF

Info

Publication number
JP7426012B2
JP7426012B2 JP2020195516A JP2020195516A JP7426012B2 JP 7426012 B2 JP7426012 B2 JP 7426012B2 JP 2020195516 A JP2020195516 A JP 2020195516A JP 2020195516 A JP2020195516 A JP 2020195516A JP 7426012 B2 JP7426012 B2 JP 7426012B2
Authority
JP
Japan
Prior art keywords
quantum
energy
derivative
quantum circuit
computer
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2020195516A
Other languages
English (en)
Other versions
JP2021028847A (ja
Inventor
光祐 御手洗
裕也 中川
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
QunaSys Inc
Original Assignee
QunaSys Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Priority claimed from JP2019090332A external-priority patent/JP6804113B2/ja
Application filed by QunaSys Inc filed Critical QunaSys Inc
Priority to JP2020195516A priority Critical patent/JP7426012B2/ja
Publication of JP2021028847A publication Critical patent/JP2021028847A/ja
Application granted granted Critical
Publication of JP7426012B2 publication Critical patent/JP7426012B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Optical Modulation, Optical Deflection, Nonlinear Optics, Optical Demodulation, Optical Logic Elements (AREA)

Description

開示の技術は、エネルギーの微分を求めるための量子情報処理方法、古典コンピュータ、量子コンピュータ、量子情報処理プログラム、及びデータ構造に関する。
従来、変分量子固有値計算（VQE : Variational-Quantum-Eigensolver）（以下、単に「ＶＱＥ」と称する。）が知られている（例えば、非特許文献１を参照。）。ＶＱＥは、量子回路のパラメータを変分的に更新することにより、ハミルトニアンの最小の固有値を近似的に計算する。ハミルトニアンの最小の固有値に対応するエネルギーの情報は、例えば、量子化学計算の際の有用な情報である。
量子化学計算において、時間に依存しない多くの物理的性質又は化学的性質は、エネルギーの導関数を用いて定義されることが多い。この点、量子コンピュータによる量子計算を用いて量子位相推定を行うことにより、エネルギーの導関数を計算する技術が知られている（例えば、非特許文献２を参照。）
A. Peruzzo, J. McClean, P. Shadbolt, M.-H. Yung, X.-Q. Zhou, P. J. Love, A. Aspuru-Guzik and J. L.O’Brien, "A variational eigenvalue solver on a photonic quantum processor," Nature Communications,5, article number: 4213, 2014. I. Kassal and A. Aspuru-Guzik, "Quantum algorithm for molecular properties and geometry optimization", J. Chem. Phys. 131, 224102 (2009);
しかし、上記非特許文献１に開示されているＶＱＥは、エネルギーの導関数の導出については考慮されていない。また、上記非特許文献２に開示されている技術は、エネルギーの導関数を計算する際に量子位相推定を利用しており、ＶＱＥを用いてエネルギーを計算する際にエネルギーの導関数を計算することは開示されていない。
開示の技術は、上記の事情を鑑みてなされたものであり、ＶＱＥを用いて系のエネルギーを量子計算する際に、エネルギーの導関数を得ることができる、エネルギーの微分を求めるための量子情報処理方法、古典コンピュータ、量子コンピュータ、量子情報処理プログラム、及びデータ構造を提供することを目的とする。
上記の目的を達成するために本開示の第１態様のエネルギーの微分を求めるための量子情報処理方法は、古典コンピュータと量子コンピュータとを含むハイブリッドシステムが実行する量子情報処理方法であって、前記古典コンピュータが、ハミルトニアンと量子回路を表すパラメータの初期情報とを出力し、前記量子コンピュータが、前記古典コンピュータから出力された、前記ハミルトニアンと前記初期情報とに応じて、ＶＱＥ（Variational-Quantum-Eigensolver）を用いた量子計算を実行し、エネルギーを計算するための第１の量子回路を表すパラメータを生成し、前記第１の量子回路を表すパラメータを出力し、前記古典コンピュータが、前記量子コンピュータから出力された前記第１の量子回路を表すパラメータに応じて、回転ゲートを含む第２の量子回路を表すパラメータを生成し、前記第２の量子回路を表すパラメータを出力し、前記量子コンピュータが、前記古典コンピュータから出力された前記第２の量子回路を表すパラメータに応じた量子計算を実行し、前記量子計算の測定結果を出力し、前記古典コンピュータが、前記量子コンピュータから出力された前記測定結果と、前記ハミルトニアンと、前記ハミルトニアンの導関数とに応じて、前記ハミルトニアンに対応するエネルギーの導関数を計算し、前記エネルギーの導関数を出力する、処理を含む。
本開示の第２態様は、エネルギーの微分を求めるための量子情報処理方法であって、前記ハミルトニアンは、物質のハミルトニアンである。
本開示の第３態様は、エネルギーの微分を求めるための量子情報処理方法であって、前記古典コンピュータは、前記第１の量子回路に含まれる複数のユニタリゲートの間に、前記回転ゲートを挿入することにより、前記第２の量子回路を表すパラメータを生成する。
本開示の第４態様は、エネルギーの微分を求めるための量子情報処理方法であって、前記古典コンピュータは、基底状態のエネルギーの前記導関数を計算する場合に、以下の式（１）によって表される回転ゲートＲ± ａ，μを、前記第１の量子回路に含まれる複数のユニタリゲートの間に挿入する。






（１）
ここで、ａは量子回路のパラメータベクトルθの要素を識別するためのインデックスａ，ｂ，ｃ，・・・であり、μはパウリ演算子を識別するためのインデックスである。
本開示の第５態様は、エネルギーの微分を求めるための量子情報処理方法であって、前記古典コンピュータは、励起状態のエネルギーの前記導関数を計算する場合に、以下の式（２）によって表される回転ゲートＲ（ｓ），± ａ，μを、前記第１の量子回路に含まれる複数のユニタリゲートの間に挿入する。






（２）
ここで、ａは量子回路のパラメータベクトルθの要素を識別するためのインデックスａ，ｂ，ｃ，・・・であり、μはパウリ演算子を識別するためのインデックスであり、ｓは励起状態を表すインデックスである。
本開示の第６態様は、エネルギーの微分を求めるための量子情報処理方法であって、前記古典コンピュータと前記量子コンピュータとはコンピュータネットワークを介して接続されており、前記古典コンピュータと前記量子コンピュータとは、前記コンピュータネットワークを介して情報の送受信を行う。
本開示の第７態様は、エネルギーの微分を求めるための量子情報処理方法であって、古典コンピュータが、ハミルトニアンと量子回路を表すパラメータの初期情報とを出力し、量子コンピュータから出力された、前記ハミルトニアンと前記初期情報とに基づきＶＱＥ（Variational-Quantum-Eigensolver）を用いた量子計算により生成された第１の量子回路を表すパラメータに応じて、回転ゲートを含む第２の量子回路を表すパラメータを生成し、前記第２の量子回路を表すパラメータを出力し、前記量子コンピュータから出力された、前記第２の量子回路を表すパラメータに応じた量子計算の測定結果と、前記ハミルトニアンと、前記ハミルトニアンの導関数とに応じて、前記ハミルトニアンに対応するエネルギーの導関数を生成し、前記エネルギーの導関数を出力する、処理を実行する古典コンピュータによる量子情報処理方法である。
本開示の第８態様は、エネルギーの微分を求めるための量子情報処理方法であって、量子コンピュータが、古典コンピュータから出力された、ハミルトニアンＨと量子回路を表すパラメータの初期情報とに応じて、ＶＱＥ（Variational-Quantum-Eigensolver）を用いた量子計算を実行し、第１の量子回路を表すパラメータを生成し、前記第１の量子回路を表すパラメータを出力し、前記古典コンピュータから出力された、前記第１の量子回路を表すパラメータに応じた回転ゲートを含む第２の量子回路を表すパラメータに基づいて、量子計算を実行し、前記量子計算の測定結果を出力する、処理を実行する量子コンピュータによる量子情報処理方法である。
本開示の第９態様は、古典コンピュータであって、ハミルトニアンと量子回路を表すパラメータの初期情報とを出力し、量子コンピュータから出力された、前記ハミルトニアンと前記初期情報とに基づきＶＱＥ（Variational-Quantum-Eigensolver）を用いた量子計算により生成された第１の量子回路を表すパラメータに応じて、回転ゲートを含む第２の量子回路を表すパラメータを生成し、前記第２の量子回路を表すパラメータを出力し、前記量子コンピュータから出力された、前記第２の量子回路を表すパラメータに応じた量子計算の測定結果と、前記ハミルトニアンと、前記ハミルトニアンの導関数とに応じて、前記ハミルトニアンに対応するエネルギーの導関数を生成し、前記エネルギーの導関数を出力する、処理を実行する古典コンピュータである。
本開示の第１０態様は、量子コンピュータであって、古典コンピュータから出力された、ハミルトニアンと量子回路を表すパラメータの初期情報とに応じて、ＶＱＥ（Variational-Quantum-Eigensolver）を用いた量子計算を実行し、エネルギーを計算するための第１の量子回路を表すパラメータを生成し、前記第１の量子回路を表すパラメータを出力し、前記古典コンピュータから出力された、前記第１の量子回路を表すパラメータに応じた回転ゲートを含む第２の量子回路を表すパラメータに基づいて、量子計算を実行し、前記量子計算の測定結果を出力する、処理を実行する量子コンピュータである。
本開示の第１１態様は、量子情報処理プログラムであって、ハミルトニアンと量子回路を表すパラメータの初期情報とを出力し、量子コンピュータから出力された、前記ハミルトニアンと前記初期情報とに基づきＶＱＥ（Variational-Quantum-Eigensolver）を用いた量子計算により生成された、エネルギーを計算するための第１の量子回路を表すパラメータに応じて、回転ゲートを含む第２の量子回路を表すパラメータを生成し、前記第２の量子回路を表すパラメータを出力し、前記量子コンピュータから出力された、前記第２の量子回路を表すパラメータに応じた量子計算の測定結果と、前記ハミルトニアンと、前記ハミルトニアンの導関数とに応じて、前記ハミルトニアンに対応するエネルギーの導関数を生成し、前記エネルギーの導関数を出力する、処理を古典コンピュータに実行させるための量子情報処理プログラムである。
本開示の第１２態様は、量子情報処理プログラムであって、古典コンピュータから出力された、ハミルトニアンと量子回路を表すパラメータの初期情報とに応じて、ＶＱＥ（Variational-Quantum-Eigensolver）を用いた量子計算を実行し、エネルギーを計算するための第１の量子回路を表すパラメータを生成し、前記第１の量子回路を表すパラメータを出力し、前記古典コンピュータから出力された、前記第１の量子回路を表すパラメータに応じた回転ゲートを含む第２の量子回路を表すパラメータに基づいて、量子計算を実行し、前記量子計算の測定結果を出力する、処理を量子コンピュータに実行させるための量子情報プログラムである。
本開示の第１３態様は、量子コンピュータの量子計算に用いられるデータ構造であって、ハミルトニアンと、量子回路を表すパラメータの初期情報と、を含み、前記ハミルトニアンと前記初期情報とに基づくＶＱＥ（Variational-Quantum-Eigensolver）を用いた前記量子コンピュータによる量子計算によって、基底状態又は励起状態のエネルギーを計算するための第１の量子回路のパラメータを生成する処理に用いられる、データ構造である。
本開示の第１４態様は、量子コンピュータの量子計算に用いられる量子回路の構成に関するデータ構造であって、複数のユニタリゲートと、複数のユニタリゲートの間に挿入された回転ゲートと、を含み、前記複数のユニタリゲートと前記回転ゲートとを含む量子回路に基づく前記量子コンピュータによる量子計算によって、前記量子計算の測定結果を測定する処理に用いられる、データ構造である。
開示の技術によれば、ＶＱＥを用いて系のエネルギーを量子計算する際に、エネルギーの導関数を得ることができる、という効果が得られる。
本実施形態のハイブリッドシステム１００の概略構成の一例を示す図である。 古典コンピュータ１１０、制御装置１２１、及びユーザ端末１３０として機能するコンピュータの概略ブロック図である。 第２の量子回路の一例を模式的に示す図である。 第２の量子回路の一例を模式的に示す図である。 系の状態を表すパラメータによるエネルギーの導関数の計算処理の一例を示す図である。 シミュレーションに用いられる仮の量子回路を示す図である。 シミュレーション結果の一例を示す図である。
以下、図面を参照して開示の技術の実施形態を詳細に説明する。
＜第１実施形態に係るハイブリッドシステム１００＞
図１に、第１実施形態に係るハイブリッドシステム１００を示す。本実施形態のハイブリッドシステム１００は、古典コンピュータ１１０と量子コンピュータ１２０とユーザ端末１３０とを備える。古典コンピュータ１１０と量子コンピュータ１２０とユーザ端末１３０とは、図１に示されるように、一例としてＩＰネットワークなどのコンピュータネットワークを介して接続されている。
本実施形態のハイブリッドシステム１００においては、量子コンピュータ１２０が古典コンピュータ１１０からの要求に応じて所定の量子計算を行い、当該量子計算の計算結果を古典コンピュータ１１０へ出力する。古典コンピュータ１１０はユーザ端末１３０へ量子計算に応じた計算結果を出力する。これにより、ハイブリッドシステム１００全体として所定の計算処理が実行される。
古典コンピュータ１１０は、通信インターフェース等の通信部１１１と、プロセッサ、ＣＰＵ（Central processing unit）等の処理部１１２と、メモリ、ハードディスク等の記憶装置又は記憶媒体を含む情報記憶部１１３とを備え、各処理を行うためのプログラムを実行することによって構成されている。なお、古典コンピュータ１１０は１又は複数の装置ないしサーバを含むことがある。また、当該プログラムは１又は複数のプログラムを含むことがあり、また、コンピュータ読み取り可能な記憶媒体に記録して非一過性のプログラムプロダクトとすることできる。
量子コンピュータ１２０は、一例として、古典コンピュータ１１０から送信される情報に基づいて量子ビット群１２３のうちの少なくとも何れかの量子ビットへ照射するための電磁波を生成する。そして、量子コンピュータ１２０は、生成された電磁波を、量子ビット群１２３のうちの少なくとも何れかの量子ビットへ照射することにより、量子回路を実行する。
図１の例では、量子コンピュータ１２０は、古典コンピュータ１１０と通信を行う制御装置１２１と、制御装置１２１からの要求に応じて電磁波を生成する電磁波生成装置１２２と、電磁波生成装置１２２からの電磁波照射を受ける量子ビット群１２３とを備える。なお、本実施形態において「量子コンピュータ」とは、古典ビットによる演算を一切行わないことを意味するものではなく、量子ビットによる演算を含むコンピュータをいう。
制御装置１２１は、古典ビットにより演算を行う古典コンピュータであり、古典コンピュータ１１０において行うものとして本明細書にて説明する処理の一部又は全部を代替的に行う。例えば、制御装置１２１は、量子回路を予め記憶又は決定しておき、量子回路Ｕ（θ）のパラメータθを受信したことに応じて、量子ビット群１２３において量子回路Ｕ（θ）を実行するための量子ゲート情報を生成してもよい。
ユーザ端末１３０は、古典ビットにより演算を行う古典コンピュータである。ユーザ端末１３０は、ユーザから入力された情報を受け付け、当該情報に応じた処理を実行する。
古典コンピュータ１１０、制御装置１２１、及びユーザ端末１３０は、例えば、図２に示すコンピュータ５０で実現することができる。コンピュータ５０はＣＰＵ５１、一時記憶領域としてのメモリ５２、及び不揮発性の記憶部５３を備える。また、コンピュータ５０は、外部装置及び出力装置等が接続される入出力interface（Ｉ／Ｆ）５４、及び記録媒体５９に対するデータの読み込み及び書き込みを制御するread/write（Ｒ／Ｗ）部５５を備える。また、コンピュータ５０は、インターネット等のネットワークに接続されるネットワークＩ／Ｆ５６を備える。ＣＰＵ５１、メモリ５２、記憶部５３、入出力Ｉ／Ｆ５４、Ｒ／Ｗ部５５、及びネットワークＩ／Ｆ５６は、バス５７を介して互いに接続される。
第１実施形態のハイブリッドシステム１００は、基底状態のエネルギーの導関数を算出する。以下、前提となる事項について説明する。
［ＶＱＥについて］
まず、ＶＱＥのアルゴリズムについて簡単に説明する。
量子コンピュータは、量子回路Ｕ（θ）に基づいて量子計算を実行する。量子回路Ｕ（θ）と量子状態｜ψ（θ）＞との間の関係は、次式によって表される。なお、ψ（θ）は波動関数を表す。

なお、次式の項は、初期化されたｎ[qubit]の状態を表す。

θは、次式に示されるように、量子回路を表すＮ次元のパラメータベクトルである。また、以下において、θｐ（ｐ＝ａ，ｂ，ｃ，・・・）は、量子回路Ｕ（θ）のパラメータベクトルθの要素を表す。量子回路Ｕ（θ）のパラメータθは、量子回路の構成を表す情報である。このため、量子コンピュータは、量子回路Ｕ（θ）のパラメータθに応じて量子計算を実行する。

ＶＱＥは、与えられたハミルトニアンＨに対応するエネルギーＥ（θ）＝＜ψ（θ）｜Ｈ｜ψ（θ）＞が最小となるように、量子回路Ｕ（θ）のパラメータθを最適化する。なお、最適な量子回路Ｕ（θ）のパラメータθにおいては、全てのａについて、以下の式（１）が成立する。


（１）
ここで、量子回路Ｕ（θ）の最適なパラメータをθ*とする。以下の式（２）の表現が採用された場合、上記式（１）は、以下の式（３）によって表される。


（２）


（３）
また、量子回路Ｕ（θ）のパラメータθによる波動関数ψ（θ）の高次の偏導関数は、以下の式（４）によって表される。


（４）
時間に依存しない物理的特性又は化学的特性の多くは、系の状態を表すパラメータｘによるエネルギーの導関数に基づき計算される。系の状態を表すパラメータｘは、例えば、系における電場、磁場、及び原子核の位置座標等を表すパラメータである。本実施形態においては、以下の式に示されるように、系の状態を表すパラメータｘがＭ次元ベクトルによって表される。なお、ｘｑ（ｑ＝ｉ，ｊ，ｋ，・・・）は、系の状態を表すパラメータベクトルｘの要素を表す。

ハミルトニアンＨ及び量子回路Ｕ（θ）の最適なパラメータθ*は、系の状態を表すパラメータｘの関数である。そのため、ハミルトニアンＨは、Ｈ（ｘ）と表される。また、量子回路Ｕ（θ）の最適なパラメータθ*は、θ*（ｘ）と表される。このため、系のエネルギーＥは以下の式（５）によって表される。


（５）
ここで、基底状態のエネルギーをＥ*（ｘ）とする場合、Ｅ*（ｘ）＝Ｅ*（θ*（ｘ），ｘ）と表される。本実施形態では、次式に示されるような、系の状態を表すパラメータｘによるＥ*（θ*（ｘ），ｘ）の導関数を算出する。

,

,

［導関数の解析的な表現について］
次に、エネルギーの導関数の解析的な表現について以下説明する。基底状態のエネルギーの導関数は、以下の式（６）～（８）によって表される。


（６）


（７）


（８）
なお、上記式（６）～（８）においては、次式が仮定される。

また、最適なパラメータθ＊（ｘ）のパラメータｘによる偏導関数は、以下の式（９）～（１０）を解くことによって求められる。


（９）


（１０）
なお、上記式（１０）におけるγは、以下の式（１１）によって表される。


（１１）
［基底状態における導関数の計算と測定について］
本実施形態のハイブリッドシステム１００は、ｎ[qubit]で動作するものとする。また、ハミルトニアンＨは、次式に示されるパウリ演算子の集合に含まれるパウリ演算子Ｉ，Ｘ，Ｙ，Ｚの和へ分解される。

このため、ハミルトニアンＨは、以下の式（１２）によって表される。


（１２）
なお、上記式（１２）におけるｈＰ（ｘ）は、次式に示される係数である。

エネルギーの導関数を計算するためには、以下の式（１３）に示されるような、ハミルトニアンＨの導関数を計算する必要がある。

,

,


（１３）
上記式（１３）に示されるようなハミルトニアンＨの導関数は、古典コンピュータによって計算可能である。なお、上記式（１３）の計算は、上記式（１２）におけるｈＰ（ｘ）の導関数に相当する。
［量子回路の表現形式］
次に、量子回路の表現形式について説明する。なお、本実施形態では、量子回路Ｕ（θ）は、以下の式（１４）に示されるように、ユニタリ行列の積によって表されるとする。


（１４）
また、各ユニタリ行列Ｕａ（θａ）は、次式に示される生成子Ｇａによって生成される。なお、ｉは虚数を表す。

上記の生成子Ｇａは、以下の式（１５）によって表される。




（１５）
［２階偏導関数の測定］
次に、量子回路Ｕ（θ）のパラメータθによるエネルギーＥの２階偏導関数の測定について説明する。
系の状態を表すパラメータｘによるエネルギーＥの導関数を求めるためには、次式に示されるような、量子回路Ｕ（θ）のパラメータθによるエネルギーＥの導関数に関する情報を得る必要がある。



なお、次式に示されるような、量子回路Ｕ（θ）のパラメータθによるエネルギーＥの１階偏導関数は、参考文献１（Y. Li and S. C. Benjamin, "Efficient Variational Quantum Simulator Incorporating Active Error Minimization ", Phys. Rev. X 7, 021050 （2017）.）及び参考文献２（K Mitarai, M Negoro, M Kitagawa, K Fujii, "Quantum circuit learning", Physical Review A 98, 032309, 2018）に開示されている技術によって計算可能である。

そのため、次式に示されるような、量子回路Ｕ（θ）のパラメータθによるエネルギーＥの２階偏導関数の計算について以下説明する。

上記の２階偏導関数は、以下の式（１６）によって表される。


（１６）
上記式（１６）における｜∂ａψ（θ）＞は、以下の式（１７）によって表される。


（１７）
また、上記式（１６）における｜∂ａ∂ｂψ（θ）＞は、以下の式（１８）によって表される。


（１８）
なお、ここでは、以下の式（１９）に示されるような定義をおく。


（１９）
ここで、以下の式（２０）が成立する。


（２０）
図３に、本実施形態の量子回路を説明するための説明図を示す。図３に示される量子回路は、上記式（２０）に示される、量子回路Ｕ（θ）のパラメータθによるエネルギーＥの２階偏導関数を計算するための量子回路である。
なお、図３に示される量子回路のうちのＵについては、次式が成立する。

また、図３に示されるＲ± ａ，μ及びＲ± ｂ，νは回転ゲートを表す。回転ゲートＲ± ａ，μ及びＲ± ｂ，νは、以下の式（２１）によって表される。なお、次式において添え字として出現する「±」は、パリティに応じて決定される符号である。




（２１）
なお、上記式（２１）におけるａ，ｂは上記式（１５）等に示されるように、量子回路Ｕ（θ）のパラメータθの要素を識別するためのインデックスである。また、上記式（２１）におけるμは上記式（１５）等に示されるように、ユニタリＵａを生成する生成子に関するインデックスであり、予め設定される。
本実施形態では、量子コンピュータ１２０が、図３に示される量子回路に基づいて、次式に示される量を測定する。

このため、上記式（２０）に含まれる次式の項は、量子コンピュータ１２０によって量子計算され、その計算結果が測定される。


（２２）
したがって、上記式（２２）の量子計算の計算結果に基づき、上記式（２０）に示されるエネルギーＥの２階偏導関数が得られる。
［３階偏導関数の測定］
次に、量子回路Ｕ（θ）のパラメータθによるエネルギーＥの３階偏導関数の測定について説明する。
量子回路Ｕ（θ）のパラメータθによるエネルギーＥの３階偏導関数は、以下の式（２３）によって表される。


（２３）
図４に、量子回路Ｕ（θ）のパラメータθによるエネルギーＥの３階偏導関数を算出するための量子回路を示す。なお、図４に示されるＲ± ｃ，ρは、次式によって表される。

本実施形態では、量子コンピュータ１２０が、図４に示される量子回路に基づいて、次式に示される量を測定する。

このため、上記式（２３）に含まれる次式の項は、量子コンピュータ１２０によって量子計算され、その計算結果が測定される。


（２４）
したがって、上記式（２４）の量子計算の計算結果に基づき、上記式（２３）に示されるエネルギーＥの３階偏導関数が得られる。
［その他の偏導関数の測定］
次に、エネルギーＥのより高次な偏導関数の計算について説明する。次式（２５Ａ）に示されるエネルギーＥの導関数は、上述した手続きと同様の手続きによって計算される。


（２５Ａ）
上記式（２５Ａ）に示される導関数は、系の状態を表すパラメータｘによるエネルギーＥの導関数が含まれている。ここで、パラメータｘによるエネルギーＥの導関数については、例えば、上記式（２０）及び上記式（２３）に示されるｈＱに、次式に示されるような導関数を代入することにより、上記式（２５Ａ）の導関数を計算することができる。なお、次式は、ハミルトニアンの導関数に相当するため、古典コンピュータ１１０により計算可能である。

また、例えば、上記式（８）に出現するような次式（２５Ｂ）のような導関数も、同様の手続きによって計算することができる。


（２５Ｂ）
［本実施形態のハイブリッドシステム１００の動作の概要］
本実施形態のハイブリッドシステム１００は、古典コンピュータ１１０と量子コンピュータ１２０とによって、上記の各計算処理を実行する。具体的には、ハイブリッドシステム１００は、以下の手順（Ａ）～（Ｆ）に従って、エネルギーEの導関数を算出する。
（Ａ）ハイブリッドシステム１００の量子コンピュータ１２０が、ＶＱＥを実行し、量子回路Ｕ（θ）の最適なパラメータθ*（ｘ）を得る。
（Ｂ）ハイブリッドシステム１００の古典コンピュータ１１０が、上記式（１３）に示されるようなハミルトニアンＨの導関数を計算する。
（Ｃ）ハイブリッドシステム１００の古典コンピュータ１１０が、上記式（６）～（８）に従って、系の状態を表すパラメータｘについてのエネルギーＥの導関数を設定する。
（Ｄ）ハイブリッドシステム１００の量子コンピュータ１２０が、量子回路Ｕ（θ）のパラメータθによるエネルギーＥの導関数を求める。
（Ｅ）ハイブリッドシステム１００の古典コンピュータ１１０が、上記式（９）、（１０）、及び（１１）に従って、次式に示されるような、系の状態を表すパラメータｘによる量子回路Ｕ（θ）の最適なパラメータθ＊の導関数を求める。

なお、上記式（９）、（１０）、及び（１１）に含まれる次式の項は、上記（Ｄ）において求められる。このため、上記（Ｄ）において求められた次式の項を、上記式（９）、（１０）、及び（１１）へ代入することにより、系の状態を表すパラメータｘによる量子回路Ｕ（θ）のパラメータθ＊の導関数が計算される。

,

,

,

,

（Ｆ）ハイブリッドシステム１００の古典コンピュータ１１０が、上記（Ａ）～（Ｅ）において計算された各項を上記式（６）～（８）へ代入することにより、系の状態を表すパラメータｘによるエネルギーＥの導関数を得る。
［第１実施形態のハイブリッドシステム１００の動作］
次に、第１実施形態のハイブリッドシステム１００の具体的な動作について説明する。ハイブリッドシステム１００の各装置において、図５に示される各処理が実行される。
まず、ステップＳ１００において、ユーザ端末１３０は、ユーザから入力された問題情報を、古典コンピュータ１１０へ送信する。問題情報は、量子計算によって解かれる問題に関する情報であり、例えば、解析対象の物質に関する情報と当該物質の物性値に関する情報とが含まれている。解析対象の物質に関する情報の一例としては物質の分子パラメータ等が挙げられ、物性値に関する情報の一例としては解析対象の分子の誘電率等が挙げられる。
次に、ステップＳ１０２において、古典コンピュータ１１０は、ユーザ端末１３０から送信された問題情報を受信する。そして、ステップＳ１０２において、古典コンピュータ１１０は、受信した問題情報のうちの物質の分子に関する情報に基づいて、上記式（１２）に従って、当該分子の系のエネルギー状態を表すハミルトニアンＨを計算する。なお、上記式（１２）におけるｈP（ｘ）は、解析対象の物質の分子に関する情報に応じて設定される。
ステップＳ１０４において、古典コンピュータ１１０は、上記ステップＳ１０２で計算されたハミルトニアンＨに基づいて、ハミルトニアンＨの微分を計算する。具体的には、古典コンピュータ１１０は、既存の量子化学計算のライブラリソフトウェアにより、上記式（１３）に示されるハミルトニアンＨの微分を計算する。
ステップＳ１０６において、古典コンピュータ１１０は、上記ステップＳ１０２で計算されたハミルトニアンＨ及び量子回路を表すパラメータθの初期情報を出力する。具体的には、上記ステップＳ１０２で計算されたハミルトニアンＨ、量子回路Ｕ（θ）を表すパラメータの初期情報、及び最適化手法を、量子コンピュータ１２０へ送信する。最適化手法としては、例えば、Ｎｅｌｄｅｒ－Ｍｅａｄ法等が挙げられる。
ステップＳ１０６において、古典コンピュータ１１０から量子コンピュータ１２０へ送信されるデータは、量子コンピュータの量子計算に用いられるデータ構造であって、ハミルトニアンと量子回路を表すパラメータの初期情報とを含む。このデータ構造は、ハミルトニアンと初期情報とに基づくＶＱＥを用いた量子計算によって、基底状態又は励起状態のエネルギーを計算するための第１の量子回路のパラメータを生成する処理に用いられる。
ステップＳ１０８において、制御装置１２１は、上記ステップＳ１０６で古典コンピュータから送信された、ハミルトニアンＨと初期情報と最適化手法とを受信する。そして、制御装置１２１は、ハミルトニアンＨと初期情報と最適化手法とに応じて、ＶＱＥを用いた量子計算を量子コンピュータ１２０に実行させる。
具体的には、量子コンピュータ１２０は、制御装置１２１の制御に応じて、量子ビット群１２３のうちの少なくとも何れかの量子ビットへ照射するための電磁波を生成する。そして、量子コンピュータ１２０は、生成された電磁波を、量子ビット群１２３のうちの少なくとも何れかの量子ビットへ照射し、初期情報に応じた量子回路を実行することにより、最適な量子回路である第１の量子回路を表すパラメータθ＊を生成する。量子回路に含まれる各量子ゲートのゲート操作は対応する電磁波波形へと変換され、生成された電磁波が電磁波生成装置１２２によって量子ビット群１２３に照射される。そして、量子コンピュータ１２０は、第１の量子回路を表すパラメータθ＊を出力する。
ステップＳ１１０において、制御装置１２１は、上記ステップＳ１０８で得られた第１の量子回路を表すパラメータθ＊を、古典コンピュータ１１０へ送信する。
ステップＳ１１２において、古典コンピュータ１１０は、上記ステップＳ１１０で制御装置１２１から送信された第１の量子回路を表すパラメータθ＊を受信する。そして、古典コンピュータ１１０は、第１の量子回路を表すパラメータθ＊に応じて、回転ゲートを含む第２の量子回路を表すパラメータを生成する。
具体的には、古典コンピュータ１１０は、第１の量子回路を表すパラメータθ＊に応じて、図３及び図４に示されるような第２の量子回路を生成する。第２の量子回路は、最適なパラメータθ＊に対応する第１の量子回路に含まれる複数のユニタリゲートの間に、上記式（２１）に示される回転ゲートが挿入された量子回路である。
第１の量子回路に含まれる複数のユニタリゲートの間に回転ゲートが挿入されることにより、エネルギーＥの導関数を算出するための第２の量子回路が生成される。量子コンピュータ１２０によって、第２の量子回路が実行されることにより、上記式（２２）及び上記式（２４）に示される量が計算される。
ステップＳ１１４において、古典コンピュータ１１０は、上記ステップＳ１１２で生成された第２の量子回路を表すパラメータを量子コンピュータ１２０へ送信する。
ステップＳ１１４において、古典コンピュータ１１０から量子コンピュータ１２０へ送信されるデータは、量子コンピュータ１２０の量子計算に用いられる量子回路の構成に関するデータ構造であり、複数のユニタリゲートと、複数のユニタリゲートの間に挿入された回転ゲートとを含む。このデータ構造は、複数のユニタリゲートと回転ゲートとを含む量子回路に基づく量子計算によって、量子計算の測定結果を測定する処理に用いられる。
ステップＳ１１６において、制御装置１２１は、上記ステップＳ１１４で古典コンピュータから送信された、第２の量子回路を表すパラメータを受信する。そして、制御装置１２１は、第２の量子回路に応じた量子計算を量子コンピュータ１２０に実行させる。量子コンピュータ１２０は、第２の量子回路を表すパラメータに応じた量子計算を実行する。
具体的には、量子コンピュータ１２０は、制御装置１２１の制御に応じて、量子ビット群１２３のうちの少なくとも何れかの量子ビットへ照射するための電磁波を生成する。そして、量子コンピュータ１２０は、生成された電磁波を、量子ビット群１２３のうちの少なくとも何れかの量子ビットへ照射することにより、第２の量子回路を実行することにより観測される情報を測定する。そして、量子コンピュータ１２０は、量子計算により得られた測定結果を出力する。
例えば、以下の表に示すようなビット列が測定結果として得られる。以下の表では、ある量子回路が実行された際に得られるサンプル数が示されている。以下の表では、量子回路が実行された際に、「0000」の情報が「10」回、「0001」の情報が「50」回、「0010」の情報が「14」回、「0011」の情報が「12」回、「0100」の情報が「85」回得られていることが、一例として示されている。

ステップＳ１１８において、制御装置１２１は、上記ステップＳ１１６で得られた測定結果を、古典コンピュータ１１０へ送信する。
ステップＳ１２０において、古典コンピュータ１１０は、上記ステップＳ１１８で制御装置１２１から送信された測定結果を受信する。そして、古典コンピュータ１１０は、測定結果を統計処理することにより、上記式（２０）に示されるパラメータθによるエネルギーＥの２階偏導関数と、上記式（２３）に示されるパラメータθによるエネルギーＥの３階偏導関数とを計算する。具体的には、古典コンピュータ１１０は、得られたビット列を統計処理することにより偏導関数の期待値を計算する。また、古典コンピュータ１１０は、上記参考文献１に開示の技術を用いて、パラメータθによるエネルギーＥの１階偏導関数を計算する。また、古典コンピュータ１１０は、上記式（２５Ａ）及び上記式（２５Ｂ）に示されるような導関数を計算する。
ステップＳ１２２において、古典コンピュータ１１０は、上記ステップＳ１０４で計算されたハミルトニアンＨの導関数と、上記ステップＳ１２０で得られた、量子計算の測定結果に応じたパラメータθによるエネルギーＥの導関数の各々とに応じて、ハミルトニアンＨに対応する、系の状態を表すパラメータｘによるエネルギーＥの導関数を計算する。
具体的には、ステップＳ１２２において、古典コンピュータ１１０は、上記ステップＳ１０２で計算されたハミルトニアンＨと、上記ステップＳ１０４で計算されたハミルトニアンＨの導関数と、上記ステップＳ１２０で計算された各導関数とを、上記式（６）、上記式（７）、及び上記式（８）へ代入して、系の状態を表すパラメータｘによるエネルギーＥの導関数の各々を計算する。
ステップＳ１２４において、古典コンピュータ１１０は、上記ステップＳ１２２で得られた、系の状態を表すパラメータｘによるエネルギーＥの導関数の各々に基づいて、上記ステップＳ１０２で受信した問題情報のうちの物性値を計算する。これにより、ユーザ端末１３０から送信された問題情報に応じた物質の分子の物性値が得られる。
ステップＳ１２６において、古典コンピュータ１１０は、上記ステップＳ１２４で得られた物性値の計算結果をユーザ端末１３０へ送信する。
ステップＳ１２８において、ユーザ端末１３０は、上記ステップＳ１２４で古典コンピュータ１１０から送信された物性値の計算結果を受信する。
以上説明したように、第１実施形態のハイブリッドシステムは、古典コンピュータが、ハミルトニアンＨと量子回路を表すパラメータの初期情報とを出力する。そして、量子コンピュータが、古典コンピュータから出力された、ハミルトニアンと初期情報とに応じて、ＶＱＥを用いた量子計算を実行し、基底状態のエネルギーを計算するための第１の量子回路を表すパラメータを生成し、第１の量子回路を表すパラメータを出力する。そして、古典コンピュータが、量子コンピュータから出力された第１の量子回路を表すパラメータに応じて、回転ゲートを含む第２の量子回路を表すパラメータを生成し、第２の量子回路を表すパラメータを出力する。そして、量子コンピュータが、古典コンピュータから出力された第２の量子回路を表すパラメータに応じた量子計算を実行し、量子計算の測定結果を出力する。そして、古典コンピュータが、量子コンピュータから出力された測定結果と、ハミルトニアンの導関数とに応じて、ハミルトニアンに対応するエネルギーの導関数を計算し、エネルギーの導関数を出力する。これにより、ＶＱＥを用いて系のエネルギーを量子計算する際に、エネルギーの導関数を得ることができる。
また、古典コンピュータと量子コンピュータとの間の適切な役割分担により、ＶＱＥを用いて系のエネルギーを量子計算する際に、系の状態を表すパラメータｘによるエネルギーＥの導関数を効率的に得ることができる。
具体的には、古典コンピュータが、ハミルトニアンＨとハミルトニアンＨの導関数とを計算し、量子コンピュータがＶＱＥと量子回路のパラメータθによるエネルギーＥの導関数とを量子計算する。これにより、系の状態を表すパラメータｘによるエネルギーＥの導関数を効率的に得ることができる。
また、本実施形態では、図３及び図４に示されるような、浅い量子回路を用いて、系の状態を表すパラメータｘによるエネルギーＥの導関数を計算する。このような浅い量子回路は、ＮＩＳＱ（Noisy Intermediate-Scale Quantum Computer）デバイスを用いる際に有用である。このため、本実施形態によれば、量子計算の誤りを抑制しつつ、系の状態を表すパラメータｘによるエネルギーＥの導関数を効率的に得ることができる。
＜第２実施形態に係るハイブリッドシステム１００＞
次に、第２実施形態について説明する。なお、第２実施形態に係るハイブリッドシステムの構成は、第１実施形態と同様の構成となるため、同一符号を付して説明を省略する。
第２実施形態のハイブリッドシステム１００は、励起状態のエネルギーの導関数を算出する。
ＶＱＥを用いて励起状態のエネルギーを計算する手法は、以下の参考文献３～参考文献６に開示されている。
参考文献３：K. M. Nakanishi, K. Mitarai, and K. Fujii, “Subspace-search variational quantum eigensolver for excited states”, (2018), arXiv:1810.09434.
参考文献４：S. Endo, T. Jones, S. McArdle, X. Yuan, and S. Benjamin, “Variational quantum algorithms for discovering Hamiltonian spectra”, (2018), arXiv:1806.05707.
参考文献５：O. Higgott, D. Wang, and S. Brierley, “Variational Quantum Computation of Excited States”, (2018), arXiv:1805.08138.
参考文献６：J. R. McClean, M. E. Kimchi-Schwartz, J. Carter, and W. A. de Jong, “Hybrid quantum-classical hierarchy for mitigation of decoherence and determination of excited states”, Phys. Rev. A 95, 042308 (2017).
第２実施形態では、参考文献４又は参考文献５に開示されている技術を用いて、励起状態のエネルギーの導関数を計算する場合を例に説明する。
与えられたハミルトニアンの基底状態をＨ０（ｘ）とする場合、ｒ番目の励起状態のハミルトニアンＨｒ（ｘ）は、以下の式（２６）によって表される。なお、以下の式（２６）では、βｓが十分に大きい場合を想定する。また、励起状態を表すインデックスｒ＝１，２，・・・である。


（２６）
上記式（２６）において、基底状態は次式によって表される。

上記式におけるθ（０）（ｘ）は、基底状態における量子回路Ｕ（θ）の最適なパラメータを表す。第１実施形態では、基底状態における量子回路Ｕ（θ）の最適なパラメータはθ＊（ｘ）によって表されたが、第２実施形態においては、基底状態における量子回路Ｕ（θ）の最適なパラメータがθ（０）（ｘ）によって表される。また、第２実施形態では、次式が成立するものとする。なお、次式のＵ（ｒ）（θ）は、第１実施形態のＵ（θ）と同様の構造である。

第２実施形態のハイブリッドシステム１００の量子コンピュータ１２０は、２つの量子状態の内積を計算する必要がある。このため、第２実施形態の量子コンピュータ１２０は、２つの量子状態の内積を計算することが可能な量子コンピュータである。
ハミルトニアンＨｒ（ｘ）の状態｜ψ（ｒ）（θ）＞における期待値は、次式によって表される。

また、ｒ番目の励起状態における最適なエネルギーは次式によって定義される。

第２実施形態のハイブリッドシステム１００は、系の状態を表すパラメータｘによる励起状態のエネルギーＥ＊の導関数を計算する。なお、ｒ番目の励起状態のハミルトニアンＨｒにとっては、ｒ番目の励起状態のエネルギーＥｒ ＊は基底状態のエネルギーに対応する。このため、上記式（６）、（７）、及び（８）は、励起状態のエネルギーＥ＊の導関数を計算する際にも適用可能である。
上記式（６）は、系の状態を表すパラメータｘによるｒ番目の励起状態のエネルギーＥの１階偏導関数に対応する。また、上記式（７）は、系の状態を表すパラメータｘによるｒ番目の励起状態のエネルギーＥの２階偏導関数に対応する。ｒ番目の励起状態のエネルギーＥの１階偏導関数と、ｒ番目の励起状態のエネルギーＥの２階偏導関数とは、次式によって表される。

,

しかし、古典コンピュータ１１０は、励起状態のハミルトニアンＨｒの導関数を計算することができない。励起状態のハミルトニアンＨｒの導関数は、以下の式（２７）によって表される。


（２７）
上記式（６）に対して上記式（２７）式を代入した場合には、次式が出現する。第２実施形態のハイブリッドシステム１００の量子コンピュータ１２０は、次式を量子計算する。

［内積の測定］
次式に示されるような量子状態の内積は、以下の式（２８）に示される形式に展開される。



（２８）
上記式（２８）におけるφ（ｓ） （ａ，μ）は、上記式（１９）に従う。上記式（２８）におけるΣ内の各項は、以下の関係式によって計算される。


（２９）
第２実施形態の量子コンピュータ１２０は、上記式（２９）の左辺と、上記式（２９）の右辺における第１項及び第２項とを量子計算する。また、第２実施形態の古典コンピュータ１１０は、上記式（２９）の左辺と右辺における第１項及び第２項との量子計算の計算結果に基づき、上記式（２９）の右辺における第３項を計算する。
なお、上記式（２９）の左辺を計算するためには、次式の量子状態を計算する必要がある。

上記式によって表される量子状態は、以下の式（３０Ａ）に示される量子回路によって表され、簡易に生成することができる。


（３０Ａ）
上記式（３０Ａ）は、例えば、複数のユニタリゲートの間に回転ゲートを挿入することに相当する。例えば、次式に示される回転ゲートＲ（ｓ），＋ ａ，μは、ａ番目のユニタリゲートとａ－１番目のユニタリゲートとの間に挿入される。




（３０Ｂ）
量子コンピュータ１２０によって、複数のユニタリゲートの間に回転ゲートが挿入された量子回路が実行されることにより、上記式（２９）の左辺が測定される。
上記の方法は、他の項にも拡張可能である。例えば、上記式（６）に対応する以下の式（３１）を展開すると、以下の式（３２）のような項が現れる。


（３１）


（３２）
上記式（３２）を量子計算する際には、以下の式（３３）が用いられる。


（３３）
上記式（３３）の右辺の全ての項を量子コンピュータ１２０を用いて計算することにより、上記式（３２）の計算をすることができる。
上記各式の量子状態を計算するため、第２実施形態においては、次式に示される回転ゲートＲ（ｓ），± ａ，μ及びＲ（ｒ），± ａ，μがユニタリゲートの間に挿入される。

［第２実施形態のハイブリッドシステム１００の動作］
次に、第２実施形態のハイブリッドシステム１００の具体的な動作について説明する。第１実施形態と同様に、第２実施形態のハイブリッドシステム１００の各装置において、上記図５に示される各処理が実行される。
ステップＳ１００～ステップＳ１１０の各処理は、第１実施形態と同様に実行される。
ステップＳ１１２において、第２実施形態の古典コンピュータ１１０は、第１の量子回路を表すパラメータθ＊に応じて、回転ゲートを含む第２の量子回路を表すパラメータを生成する。
具体的には、第２実施形態の古典コンピュータ１１０は、最適なパラメータθ＊に対応する第１の量子回路に含まれる複数のユニタリゲートの間に、上記式（３０Ｂ）に示される回転ゲートを挿入して、第２の量子回路を生成する。
なお、ステップＳ１１４～ステップＳ１２６の各処理は、第１実施形態と同様に実行される。
以上説明したように、第２実施形態のハイブリッドシステムは、古典コンピュータが、励起状態のエネルギーの導関数を計算する場合に、励起状態ｓに応じた回転ゲートを、第１の量子回路に含まれる複数のユニタリゲートの間に挿入する。これにより、ＶＱＥを用いて系のエネルギーを量子計算する際に、励起状態のエネルギーの導関数を得ることができる。
次に、実施例を説明する。本実施例では、水素分子の電子ハミルトニアンを用いて数値シミュレーションを行った。本実施例では、既存のオープンソースライブラリであるPySCF（参考文献６（Q. Sun, T. C. Berkelbach, N. S. Blunt, G. H. Booth, S. Guo, Z. Li, J. Liu, J. D. McClain, E. R. Sayfutyarova, S. Sharma, S.Wouters, and G. K. Chan, Wiley Interdisciplinary Reviews: Computational Molecular Science 8, e1340 (2017).）を参照。）とOpenFermion（参考文献７（J. R. McClean, K. J. Sung, I. D. Kivlichan, Y. Cao, C. Dai, E. S. Fried, C. Gidney, B. Gimby, P. Gokhale, T. Hner, T. Hardikar, V. Havlek, O. Higgott, C. Huang, J. Izaac, Z. Jiang, X. Liu, S. McArdle, M. Neeley, T. O'Brien, B. O'Gorman, I. Ozdan, M. D. Radin, J. Romero, N. Rubin, N. P. D. Sawaya, K. Setia, S. Sim, D. S. Steiger, M. Steudtner, Q. Sun, W. Sun, D. Wang, F. Zhang, and R. Babbush, (2017), arXiv:1710.07629.を参照。）を用いてハミルトニアンを計算した。また、量子回路のシミュレーションは、Qulacs（Qulacs," https://github.com/qulacs/qulacs. "）を用いて実施した。
図６は数値シミュレーションに用いた量子回路である。Ｒｙ，Ｒｘはそれぞれｙ軸回転ゲート、ｘ軸回転ゲートを表している。図７は、図６の量子回路と水素分子のハミルトニアンを用いて、ハミルトニアンのパラメータｘとして水素原子間の距離として本手法を数値シミュレーションし、エネルギーの２次微分及び３次微分を求め、エネルギー曲線の２次の近似曲線及び３次の近似曲線を描いたものである。なお、図７の「Full CI」は理論値を表し、「Harmonic approx.」は２次の近似曲線を表し、「3rd-oder approx.」は３次の近似曲線を表し、「Enegy minimum」はエネルギーの最小値を表す。図７に示されるように、近似曲線が精度良く求められていることがわかる。
なお、本開示の技術は、上述した実施形態に限定されるものではなく、この発明の要旨を逸脱しない範囲内で様々な変形や応用が可能である。
例えば、上記各実施形態において、古典コンピュータ１１０と量子コンピュータ１２０との間の情報の送受信はどのようになされてもよい。例えば、古典コンピュータ１１０と量子コンピュータ１２０との間における、量子回路を表すパラメータの送受信及び測定結果の送受信等は、所定の計算が完了する毎に逐次送受信が行われてもよいし、全ての計算が完了した後に送受信が行われてもよい。
また、上記各実施形態では、ユーザ端末１３０から古典コンピュータ１１０へ問題情報が送信され、古典コンピュータ１１０が問題情報に応じたハミルトニアンＨを計算する場合を例に説明したが、これに限定されるものではない。例えば、ユーザ端末１３０を操作するユーザにおいて、ハミルトニアンとして当該問題を表現できる場合には、古典コンピュータ１１０は、ハミルトニアンＨを問題情報として受信してもよい。ユーザ端末１３０は、ＩＰネットワークなどのコンピュータネットワークを介して古典コンピュータ１１０又は古典コンピュータ１１０がアクセス可能な記憶媒体又は記憶装置に問題情報を送信してもよいが、記憶媒体又は記憶装置に記憶して古典コンピュータ１１０の運営者に渡し、当該運営者が古典コンピュータ１１０に当該記憶媒体又は記憶装置を用いて問題情報を入力するようにしてもよい。
また、上記各実施形態では、解析対象の物質に関する情報と当該物質の物性値に関する情報とが問題情報に含まれており、解析対象の物質に応じたエネルギーＥの導関数に基づき物性値が計算される場合を例に説明したが、これに限定されるものではない。例えば、所定の最適化問題に関する情報が問題情報として与えられても良い。この場合には、ハイブリッドシステム１００によって、与えられた最適化問題についての系の状態を表すパラメータｘによるエネルギーＥの導関数が計算される。例えば、最適化問題の一例である巡回セールスマン問題が問題情報として与えられた場合には、巡回対象の場所間の距離の総和を表すコスト関数がエネルギー関数Ｅとして設定され、巡回対象となる場所間の距離を表すパラメータｘによるコスト関数の導関数が計算される。
また、上記各実施形態では、上記式（２１）及び上記式（３０Ｂ）に示されるように、回転ゲートの角度がπ／４である場合を例に説明したが、これに限定されるものではなく、どのような角度であってもよい。例えば、回転ゲートの角度はπ／３であってもよい。なお、回転ゲートの角度を変更する場合には、回転ゲートの角度に応じた適切な形式に上記式（２１）及び上記式（３０Ｂ）が変形される。
また、上記各実施形態では、電磁波の照射によって量子回路が実行される場合を例に説明したが、これに限定されるものではなく、異なる方式によって量子回路が実行されてもよい。
また、上記各実施形態では、異なる組織によって古典コンピュータ１１０及び量子コンピュータ１２０が管理されている場合を想定しているが、古典コンピュータ１１０及び量子コンピュータ１２０は同一の組織によって一体として管理されていてもよい。この場合には、量子計算情報の古典コンピュータ１１０から量子コンピュータ１２０への送信及び量子コンピュータ１２０から古典コンピュータ１１０への測定結果の送信は不要となる。また、この場合には、量子コンピュータ１２０の制御装置１２１において上述の説明における古典コンピュータ１１０の役割を担うことが考えられる。
なお、上記各実施形態においては、「××のみに基づいて」、「××のみに応じて」、「××のみの場合」というように「のみ」との記載がなければ、本明細書においては、付加的な情報も考慮し得ることが想定されていることに留意されたい。一例として、「ａの場合にｂする」という記載は、明示した場合を除き、「ａの場合に常にｂする」ことを必ずしも意味しない。
また、何らかの方法、プログラム、端末、装置、サーバ又はシステム（以下「方法等」）において、本明細書で記述された動作と異なる動作を行う側面があるとしても、開示の技術の各態様は、本明細書で記述された動作のいずれかと同一の動作を対象とするものであり、本明細書で記述された動作と異なる動作が存在することは、当該方法等を本開示の技術の各態様の範囲外とするものではない。
また、本願明細書中において、プログラムが予めインストールされている実施形態として説明したが、当該プログラムを、コンピュータ読み取り可能な記録媒体に格納して提供することも可能である。
１００ ハイブリッドシステム
１１０ 古典コンピュータ
１１１ 通信部
１１２ 処理部
１１３ 情報記憶部
１２０ 量子コンピュータ
１２１ 制御装置
１２２ 電磁波生成装置
１２３ 量子ビット群
１３０ ユーザ端末

Claims (12)

  1. 古典コンピュータと量子コンピュータとを含むハイブリッドシステムが実行する量子情報処理方法であって、
    前記古典コンピュータが、ハミルトニアンと量子回路を表すパラメータの初期情報とを出力し、
    前記量子コンピュータが、前記古典コンピュータから出力された、前記ハミルトニアンと前記初期情報とに応じて、ＶＱＥ（Variational-Quantum-Eigensolver）を用いた量子計算を実行し、エネルギーを計算するための第１の量子回路を表すパラメータを生成し、前記第１の量子回路を表すパラメータを出力し、
    前記古典コンピュータが、前記量子コンピュータから出力された前記第１の量子回路を表すパラメータに応じて、回転ゲートを含む第２の量子回路を表すパラメータを生成し、前記第２の量子回路を表すパラメータを出力し、
    前記量子コンピュータが、前記古典コンピュータから出力された前記第２の量子回路を表すパラメータに応じた量子計算を実行し、前記量子計算の測定結果を出力し、
    前記古典コンピュータが、前記量子コンピュータから出力された前記測定結果と、前記ハミルトニアンと、前記ハミルトニアンの導関数とに応じて、前記ハミルトニアンに対応するエネルギーの導関数を計算し、前記エネルギーの導関数を出力し、
    前記エネルギーの導関数は、系の状態を表すパラメータに対するエネルギーの導関数であり、
    前記第２の量子回路は、量子回路のパラメータθに対する前記エネルギーの２階偏導関数を計算するための量子回路と、量子回路のパラメータθに対する前記エネルギーの３階偏導関数を計算するための量子回路とを含んで構成されており、
    前記エネルギーの導関数は、前記エネルギーの２階偏導関数の測定結果と前記エネルギーの３階偏導関数の測定結果とに基づいて計算される、
    処理を含むエネルギーの微分を求めるための量子情報処理方法。
  2. 前記ハミルトニアンは、物質のハミルトニアンである、
    請求項１に記載のエネルギーの微分を求めるための量子情報処理方法。
  3. 前記古典コンピュータは、前記第１の量子回路に含まれる複数のユニタリゲートの間に、前記回転ゲートを挿入することにより、前記第２の量子回路を表すパラメータを生成する、
    請求項１又は請求項２に記載のエネルギーの微分を求めるための量子情報処理方法。
  4. 前記古典コンピュータは、基底状態のエネルギーの前記導関数を計算する場合に、以下の式（１）によって表される回転ゲートＲ± ａ，μを、前記第１の量子回路に含まれる複数のユニタリゲートの間に挿入する、
    請求項３に記載のエネルギーの微分を求めるための量子情報処理方法。



    （１）
    ここで、ａは量子回路のパラメータベクトルθの要素を識別するためのインデックスａ，ｂ，ｃ，・・・であり、μはパウリ演算子を識別するためのインデックスである。
  5. 前記古典コンピュータは、励起状態のエネルギーの前記導関数を計算する場合に、以下の式（２）によって表される回転ゲートＲ（ｓ），± ａ，μを、前記第１の量子回路に含まれる複数のユニタリゲートの間に挿入する、
    請求項３に記載のエネルギーの微分を求めるための量子情報処理方法。



    （２）
    ここで、ａは量子回路のパラメータベクトルθの要素を識別するためのインデックスａ，ｂ，ｃ，・・・であり、μはパウリ演算子を識別するためのインデックスであり、ｓは励起状態を表すインデックスである。
  6. 前記古典コンピュータと前記量子コンピュータとはコンピュータネットワークを介して接続されており、
    前記古典コンピュータと前記量子コンピュータとは、前記コンピュータネットワークを介して情報の送受信を行う、
    請求項１～請求項５の何れか１項に記載のエネルギーの微分を求めるための量子情報処理方法。
  7. 古典コンピュータが、
    ハミルトニアンと量子回路を表すパラメータの初期情報とを出力し、
    量子コンピュータから出力された、前記ハミルトニアンと前記初期情報とに基づきＶＱＥ（Variational-Quantum-Eigensolver）を用いた量子計算により生成された第１の量子回路を表すパラメータに応じて、回転ゲートを含む第２の量子回路を表すパラメータを生成し、前記第２の量子回路を表すパラメータを出力し、
    前記量子コンピュータから出力された、前記第２の量子回路を表すパラメータに応じた量子計算の測定結果と、前記ハミルトニアンと、前記ハミルトニアンの導関数とに応じて、前記ハミルトニアンに対応するエネルギーの導関数を生成し、前記エネルギーの導関数を出力し、
    前記エネルギーの導関数は、系の状態を表すパラメータに対するエネルギーの導関数であり、
    前記第２の量子回路は、量子回路のパラメータθに対する前記エネルギーの２階偏導関数を計算するための量子回路と、量子回路のパラメータθに対する前記エネルギーの３階偏導関数を計算するための量子回路とを含んで構成されており、
    前記エネルギーの導関数は、前記エネルギーの２階偏導関数の測定結果と前記エネルギーの３階偏導関数の測定結果とに基づいて計算される、
    処理を実行するエネルギーの微分を求めるための量子情報処理方法。
  8. 量子コンピュータが、
    古典コンピュータから出力された、ハミルトニアンと量子回路を表すパラメータの初期情報とに応じて、ＶＱＥ（Variational-Quantum-Eigensolver）を用いた量子計算を実行し、第１の量子回路を表すパラメータを生成し、前記第１の量子回路を表すパラメータを出力し、
    前記古典コンピュータから出力された、前記第１の量子回路を表すパラメータに応じた回転ゲートを含む第２の量子回路を表すパラメータに基づいて、量子計算を実行し、前記量子計算の測定結果を出力する、
    処理を実行するエネルギーの微分を求めるための量子情報処理方法であって、
    前記古典コンピュータは、前記量子コンピュータから出力された前記測定結果と、前記ハミルトニアンと、前記ハミルトニアンの導関数とに応じて、前記ハミルトニアンに対応するエネルギーの導関数を計算し、前記エネルギーの導関数を出力し、
    前記エネルギーの導関数は、系の状態を表すパラメータに対するエネルギーの導関数であり、
    前記第２の量子回路は、量子回路のパラメータθに対する前記エネルギーの２階偏導関数を計算するための量子回路と、量子回路のパラメータθに対する前記エネルギーの３階偏導関数を計算するための量子回路とを含んで構成されており、
    前記エネルギーの導関数は、前記エネルギーの２階偏導関数の測定結果と前記エネルギーの３階偏導関数の測定結果とに基づいて計算される、
    量子情報処理方法。
  9. ハミルトニアンと量子回路を表すパラメータの初期情報とを出力し、
    量子コンピュータから出力された、前記ハミルトニアンと前記初期情報とに基づきＶＱＥ（Variational-Quantum-Eigensolver）を用いた量子計算により生成された第１の量子回路を表すパラメータに応じて、回転ゲートを含む第２の量子回路を表すパラメータを生成し、前記第２の量子回路を表すパラメータを出力し、
    前記量子コンピュータから出力された、前記第２の量子回路を表すパラメータに応じた量子計算の測定結果と、前記ハミルトニアンと、前記ハミルトニアンの導関数とに応じて、前記ハミルトニアンに対応するエネルギーの導関数を生成し、前記エネルギーの導関数を出力し、
    前記エネルギーの導関数は、系の状態を表すパラメータに対するエネルギーの導関数であり、
    前記第２の量子回路は、量子回路のパラメータθに対する前記エネルギーの２階偏導関数を計算するための量子回路と、量子回路のパラメータθに対する前記エネルギーの３階偏導関数を計算するための量子回路とを含んで構成されており、
    前記エネルギーの導関数は、前記エネルギーの２階偏導関数の測定結果と前記エネルギーの３階偏導関数の測定結果とに基づいて計算される、
    処理を実行する古典コンピュータ。
  10. 古典コンピュータから出力された、ハミルトニアンと量子回路を表すパラメータの初期情報とに応じて、ＶＱＥ（Variational-Quantum-Eigensolver）を用いた量子計算を実行し、エネルギーを計算するための第１の量子回路を表すパラメータを生成し、前記第１の量子回路を表すパラメータを出力し、
    前記古典コンピュータから出力された、前記第１の量子回路を表すパラメータに応じた回転ゲートを含む第２の量子回路を表すパラメータに基づいて、量子計算を実行し、前記量子計算の測定結果を出力する、
    処理を実行する量子コンピュータであって、
    前記古典コンピュータは、前記量子コンピュータから出力された前記測定結果と、前記ハミルトニアンと、前記ハミルトニアンの導関数とに応じて、前記ハミルトニアンに対応するエネルギーの導関数を計算し、前記エネルギーの導関数を出力し、
    前記エネルギーの導関数は、系の状態を表すパラメータに対するエネルギーの導関数であり、
    前記第２の量子回路は、量子回路のパラメータθに対する前記エネルギーの２階偏導関数を計算するための量子回路と、量子回路のパラメータθに対する前記エネルギーの３階偏導関数を計算するための量子回路とを含んで構成されており、
    前記エネルギーの導関数は、前記エネルギーの２階偏導関数の測定結果と前記エネルギーの３階偏導関数の測定結果とに基づいて計算される、
    量子コンピュータ。
  11. ハミルトニアンと量子回路を表すパラメータの初期情報とを出力し、
    量子コンピュータから出力された、前記ハミルトニアンと前記初期情報とに基づきＶＱＥ（Variational-Quantum-Eigensolver）を用いた量子計算により生成された、エネルギーを計算するための第１の量子回路を表すパラメータに応じて、回転ゲートを含む第２の量子回路を表すパラメータを生成し、前記第２の量子回路を表すパラメータを出力し、
    前記量子コンピュータから出力された、前記第２の量子回路を表すパラメータに応じた量子計算の測定結果と、前記ハミルトニアンと、前記ハミルトニアンの導関数とに応じて、前記ハミルトニアンに対応するエネルギーの導関数を生成し、前記エネルギーの導関数を出力し、
    前記エネルギーの導関数は、系の状態を表すパラメータに対するエネルギーの導関数であり、
    前記第２の量子回路は、量子回路のパラメータθに対する前記エネルギーの２階偏導関数を計算するための量子回路と、量子回路のパラメータθに対する前記エネルギーの３階偏導関数を計算するための量子回路とを含んで構成されており、
    前記エネルギーの導関数は、前記エネルギーの２階偏導関数の測定結果と前記エネルギーの３階偏導関数の測定結果とに基づいて計算される、
    処理を古典コンピュータに実行させるための量子情報処理プログラム。
  12. 古典コンピュータから出力された、ハミルトニアンと量子回路を表すパラメータの初期情報とに応じて、ＶＱＥ（Variational-Quantum-Eigensolver）を用いた量子計算を実行し、エネルギーを計算するための第１の量子回路を表すパラメータを生成し、前記第１の量子回路を表すパラメータを出力し、
    前記古典コンピュータから出力された、前記第１の量子回路を表すパラメータに応じた回転ゲートを含む第２の量子回路を表すパラメータに基づいて、量子計算を実行し、前記量子計算の測定結果を出力する
    処理を量子コンピュータに実行させるための量子情報プログラムであって、
    前記古典コンピュータは、前記量子コンピュータから出力された前記測定結果と、前記ハミルトニアンと、前記ハミルトニアンの導関数とに応じて、前記ハミルトニアンに対応するエネルギーの導関数を計算し、前記エネルギーの導関数を出力し、
    前記エネルギーの導関数は、系の状態を表すパラメータに対するエネルギーの導関数であり、
    前記第２の量子回路は、量子回路のパラメータθに対する前記エネルギーの２階偏導関数を計算するための量子回路と、量子回路のパラメータθに対する前記エネルギーの３階偏導関数を計算するための量子回路とを含んで構成されており、
    前記エネルギーの導関数は、前記エネルギーの２階偏導関数の測定結果と前記エネルギーの３階偏導関数の測定結果とに基づいて計算される、
    量子情報プログラム。
JP2020195516A 2019-05-13 2020-11-25 エネルギーの微分を求めるための量子情報処理方法、古典コンピュータ、量子コンピュータ、量子情報処理プログラム、及びデータ構造 Active JP7426012B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2020195516A JP7426012B2 (ja) 2019-05-13 2020-11-25 エネルギーの微分を求めるための量子情報処理方法、古典コンピュータ、量子コンピュータ、量子情報処理プログラム、及びデータ構造

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2019090332A JP6804113B2 (ja) 2019-05-13 2019-05-13 エネルギーの微分を求めるための量子情報処理方法、古典コンピュータ、量子コンピュータ、量子情報処理プログラム、及びハイブリッドシステム
JP2020195516A JP7426012B2 (ja) 2019-05-13 2020-11-25 エネルギーの微分を求めるための量子情報処理方法、古典コンピュータ、量子コンピュータ、量子情報処理プログラム、及びデータ構造

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
JP2019090332A Division JP6804113B2 (ja) 2019-05-13 2019-05-13 エネルギーの微分を求めるための量子情報処理方法、古典コンピュータ、量子コンピュータ、量子情報処理プログラム、及びハイブリッドシステム

Publications (2)

Publication Number Publication Date
JP2021028847A JP2021028847A (ja) 2021-02-25
JP7426012B2 true JP7426012B2 (ja) 2024-02-01

Family

ID=89718162

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2020195516A Active JP7426012B2 (ja) 2019-05-13 2020-11-25 エネルギーの微分を求めるための量子情報処理方法、古典コンピュータ、量子コンピュータ、量子情報処理プログラム、及びデータ構造

Country Status (1)

Country Link
JP (1) JP7426012B2 (ja)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114528996B (zh) * 2022-01-27 2023-08-04 本源量子计算科技（合肥）股份有限公司 一种目标体系试验态初始参数的确定方法、装置及介质

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2018080553A1 (en) 2016-10-24 2018-05-03 Google Llc Simulating materials using quantum computation
JP6804113B2 (ja) 2019-05-13 2020-12-23 株式会社ＱｕｎａＳｙｓ エネルギーの微分を求めるための量子情報処理方法、古典コンピュータ、量子コンピュータ、量子情報処理プログラム、及びハイブリッドシステム

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2018080553A1 (en) 2016-10-24 2018-05-03 Google Llc Simulating materials using quantum computation
JP6804113B2 (ja) 2019-05-13 2020-12-23 株式会社ＱｕｎａＳｙｓ エネルギーの微分を求めるための量子情報処理方法、古典コンピュータ、量子コンピュータ、量子情報処理プログラム、及びハイブリッドシステム

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
MITARAI, Kosuke et al.，"Methodology for replacing indirect measurements with direct measurements"，arXiv.org [online]，arXiv:1901.00015v2，Cornell University，2019年01月，pp.1-5，[令和２年８月７日検索], インターネット:<URL: https://arxiv.org/pdf/1901.00015v2>
ROMERO, Jonathan et al.，"Strategies for quantum computing molecular energies using the unitary coupled cluster ansatz"，arXiv.org [online]，arXiv:1701.02691v2，Cornell University，2018年，pp.1-18，[令和２年８月７日検索], インターネット:<URL: https://arxiv.org/pdf/1701.02691v2>

Also Published As

Publication number Publication date
JP2021028847A (ja) 2021-02-25

Similar Documents

Publication Publication Date Title
JP6804113B2 (ja) エネルギーの微分を求めるための量子情報処理方法、古典コンピュータ、量子コンピュータ、量子情報処理プログラム、及びハイブリッドシステム
McArdle et al. Digital quantum simulation of molecular vibrations
O’Brien et al. Calculating energy derivatives for quantum chemistry on a quantum computer
Tranter et al. A comparison of the Bravyi–Kitaev and Jordan–Wigner transformations for the quantum simulation of quantum chemistry
Kyriienko Quantum inverse iteration algorithm for programmable quantum simulators
Atia et al. Fast-forwarding of Hamiltonians and exponentially precise measurements
Daskin et al. Decomposition of unitary matrices for finding quantum circuits: application to molecular Hamiltonians
Schleife et al. Plane-wave pseudopotential implementation of explicit integrators for time-dependent Kohn-Sham equations in large-scale simulations
Ruiz-Serrano et al. A variational method for density functional theory calculations on metallic systems with thousands of atoms
Zelovich et al. Parameter-free driven Liouville-von Neumann approach for time-dependent electronic transport simulations in open quantum systems
Zuehlsdorff et al. Linear-scaling time-dependent density-functional theory beyond the Tamm-Dancoff approximation: obtaining efficiency and accuracy with in situ optimised local orbitals
McArdle et al. Improving the accuracy of quantum computational chemistry using the transcorrelated method
Kowalski et al. Quantum simulations employing connected moments expansions
Sanghi et al. Coarse-grained potential models for structural prediction of carbon dioxide (CO2) in confined environments
Whitfield Communication: Spin-free quantum computational simulations and symmetry adapted states
Man et al. Coherent states formulation of polymer field theory
JP7426012B2 (ja) エネルギーの微分を求めるための量子情報処理方法、古典コンピュータ、量子コンピュータ、量子情報処理プログラム、及びデータ構造
Huang et al. Protecting quantum correlation from correlated amplitude damping channel
Bian et al. Quantum computing for atomic and molecular resonances
Saitow et al. A local pair-natural orbital-based complete-active space perturbation theory using orthogonal localized virtual molecular orbitals
JP7515866B2 (ja) 遷移振幅を計算するための量子情報処理方法、古典コンピュータ、量子コンピュータ、ハイブリッドシステム、及び量子情報処理プログラム
Tsuchimochi et al. Multi-state quantum simulations via model-space quantum imaginary time evolution
Li Absorbing boundary conditions for time-dependent Schrödinger equations: A density-matrix formulation
Cao et al. Lindblad equation and its semiclassical limit of the Anderson-Holstein model
Allan et al. Taming the third order cumulant approximation to linear optical spectroscopy

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20220405

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20230627

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20230828

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20231019

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20231212

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20240110

R150 Certificate of patent or registration of utility model

Ref document number: 7426012

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150