Quantum computers have advanced rapidly, and quantum chemistry is now regarded as one of the most promising applications of quantum computing. The considerable attention to this field originates from the fact that the quantum phase estimation (QPE) algorithm [1,2] can determine eigenstates of the electronic structure Hamiltonian in polynomial time, whereas the best-known classical algorithms scale exponentially [3]. Consequently, reducing the resource cost of fault-tolerant quantum computation has become a central goal of the field [4-9]. Ongoing progress in hardware and algorithms promises to extend quantum-chemical accuracy to a wide range of problems in chemistry and materials science. 量子コンピュータは急速に発展しており、量子化学は量子コンピュータの最も有望な応用分野の一つと考えられている。この分野に大きな注目が集まっているのは、量子位相推定(QPE)アルゴリズム[1,2]が電子構造ハミルトニアンの固有状態を多項式時間で決定できるのに対し、最もよく知られた古典的アルゴリズムは指数関数的にスケーリングする[3]という事実に端を発している。その結果、フォールト・トレラント量子計算のリソース・コストを削減することが、この分野の中心的な目標となっている[4-9]。ハードウェアやアルゴリズムの進歩により、量子化学の精度は化学や材料科学における様々な問題にまで拡大することが期待されている。
In contrast to the long-term, fault-tolerant scenario outlined above, recent demonstrations of quantumchemical calculations using current quantum hardware have steadily approached the accuracy and accessible system sizes of advanced conventional electronic-structure methods. These advances have been enabled by hybrid approaches that embed quantum computationscapturing static (strong) correlation through coherent superpositions - within classical post-Hartree-Fock procedures that recover dynamical (weak) correlation, thus enabling quantitative accuracy. The integration of quantum computation with auxiliary-field quantum Monte Carlo (AFQMC) [10-12] and coupled cluster (CC) [13, 14] exemplifies this strategy. 上記のような長期的でフォールトトレラントなシナリオとは対照的に、現在の量子ハードウェアを用いた最近の量子化学計算の実証実験では、従来の高度な電子構造法の精度と利用可能なシステムサイズに着実に近づいている。このような進歩は、動的相関(弱い相関)を回復する古典的なポストハートリーフォック法の中に、コヒーレントな重ね合わせによって静的相関(強い相関)をとらえる量子計算を組み込んだハイブリッド・アプローチによって可能になった。量子計算と補助場量子モンテカルロ法(AFQMC)[10-12]および結合クラスタ法(CC)[13, 14]の統合は、この戦略の一例である。
CC theory is one of the most successful methods in modern quantum chemistry [15,16][15,16]. Widely imple- CC 理論は、現代の量子化学 [15,16][15,16] で最も成功した手法の一つである。広く実装されている
mented in quantum chemistry software, CC excels at describing dynamical correlation and typically yields high accuracy for single-reference (SR) systems - those lacking significant orbital quasi-degeneracy. Its polynomial computational cost is attractive: CC singles and doubles (CCSD) scales as O(N^(6))O\left(N^{6}\right), while the perturbative triples extension, CCSD(T)\operatorname{CCSD}(\mathrm{T}), scales as O(N^(7))O\left(N^{7}\right). Further developments, such as linear-scaling CC methods [17-23], have extended CC to relatively small proteins. In addition, CC data now serve as valuable training data for machine learning [24,25][24,25]. Therefore, CC theory is significant for SR systems and has widespread applicability. CCは量子化学ソフトウェアに組み込まれており、動的相関の記述に優れ、通常、単一参照(SR)系(軌道の準-縮退が顕著でない系)で高い精度が得られます。計算コストが多項式であることも魅力です:CCシングル&ダブルス(CCSD)は O(N^(6))O\left(N^{6}\right) としてスケールするが,摂動トリプルの拡張である CCSD(T)\operatorname{CCSD}(\mathrm{T}) は O(N^(7))O\left(N^{7}\right) としてスケールする。線形スケーリングCC 法[17-23]などのさらなる発展により、CCは比較的小さなタンパク質にまで拡張された。さらに、CCデータは現在、機械学習 [24,25][24,25] の貴重な学習データとして役立っている。したがって、CC 理論はSRシステムにとって重要であり、広く応用可能である。
However, conventional CC cannot describe multireference (MR) systems, where quasi-degenerate orbitals yield strong correlation. The tailored CC (TCC) method [26,27] addresses this by embedding amplitudes obtained from an MR wave function into the CC ansatz. While early implementations employed complete active space configuration interaction (CASCI) [26, 27], recent studies have incorporated other MR solvers, such as density matrix renormalization group [28, 29] and full configuration interaction quantum Monte Carlo [30, 31], achieving a balanced treatment of strong and weak correlations. しかし、従来のCCでは、準縮退軌道が強い相関をもたらす多参照(MR)系を記述することができない。テーラードCC(TCC)法[26,27]は、MR 波動関数から得られる振幅をCCアサッツに埋め込むことでこの問題に対処する。初期の実装では完全活性空間配置相互作用(CASCI)[26,27]が採用されていましたが、最近の研究では密度行列繰り込み群[28,29]や完全配置相互作用量子モンテカルロ[30,31]などの他のMRソルバーが取り入れられ、強相関と弱相関のバランスの取れた取り扱いが実現されています。
The integration of quantum computation into TCC has recently been proposed [13, 14]. This provides a significant advantage because it avoids estimating higher-order reduced density matrices (RDMs), which are expensive to obtain on a quantum device [32]. By contrast, the quantum state realized on a quantum computer is reconstructed via state-tomography techniques and converted directly into CC amplitudes. Two initial studies employed computational basis tomography [13] and matchgate shadows tomography [14], demonstrating balanced 最近、量子計算をTCCに統合することが提案されている[13, 14]。これは、量子デバイスで得るには高価な高次の縮小密度行列(RDM)を推定する必要がないため、大きな利点となる[32]。対照的に、量子コンピュータ上で実現される量子状態は、状態トモグラフィ技術によって再構成され、直接 CC 振幅に変換される。計算基底トモグラフィー[13]とマッチゲートシャドウトモグラフィー[14]を用いた2つの初期研究があり、バランスの取れた
accuracy but highlighting the need for efficient sampling from quantum states. しかし、量子状態からの効率的なサンプリングの必要性を強調している。
Quantum-selected configuration interaction (QSCI) [33] is a scalable state tomography protocol that can work with relatively small shot budgets and is highly flexible. It can be combined with various quantum algorithms, such as the variational quantum eigensolver (VQE) [34] and QPE, to determine the eigenstate of the electronic structure Hamiltonian. QSCI has expanded quantum-classical hybrid calculations to 77 qubits [35], a scale that already enables practical quantum-chemical applications. Note that QSCI has been sometimes referred to as sampling-based quantum diagonalization (SQD) [35]; however, this study treats them as essentially identical. In QSCI, an exact quantum state is not always necessary, which reduces some workloads of quantum devices [33]. Furthermore, the hardness of classical sampling from unitary cluster Jastrow circuits suggests a potential quantum advantage in QSCI [36]. 量子選択配置相互作用(QSCI)[33]は、スケーラブルな状態トモグラフィープロトコルであり、比較的小さなショットバジェットで動作し、非常に柔軟である。これは、変分量子固有値解法(VQE)[34]やQPEなどの様々な量子アルゴリズムと組み合わせて、電子構造ハミルトニアンの固有状態を決定することができます。QSCIは、量子-古典ハイブリッド計算を77 量子ビット[35]まで拡張しており、これはすでに実用的な量子化学応用を可能にする規模である。なお、QSCIはサンプリングに基づく量子対角化(SQD)と呼ばれることもあります[35]。QSCIでは、厳密な量子状態は必ずしも必要ではないため、量子デバイスの負荷を軽減することができます[33]。さらに、ユニタリー・クラスター・ジャストロー回路からの古典的なサンプリングの硬さは、QSCIにおける量子的な利点の可能性を示唆しています[36]。
Moreover, QSCI offers a significant advantage: its configuration-interaction (CI) coefficients are free from the additive shot noise that plagues many other tomography schemes. In QSCI, computational-basis states are first sampled from the prepared quantum state; although this sampling is stochastic, it affects only which determinants enter the subspace. The subsequent stepsconstructing the effective Hamiltonian within that subspace and diagonalizing it-are performed exactly. Consequently, the sole error in the CI coefficients stems from subspace truncation, not from statistical noise. This feature contrasts sharply with computational basis tomography (sampling) [37] and classical shadows methods [3840], both of which propagate shot noise directly into the reconstructed state. A noise-free wave function is particularly beneficial for downstream classical algorithms. For example, a quantum-classical hybrid quantum Monte Carlo method was performed stably [11], as opposed to additive errors that had prevented stable computation. Therefore, it makes QSCI a powerful partner for advanced correlated wave-function methods. さらに、QSCIには、他の多くのトモグラフィ方式を悩ませている加法性ショットノイズがない、という大きな利点がある。QSCIでは、まず準備された量子状態から計算基底状態をサンプリングする。このサンプリングは確率的であるが、どの行列式が部分空間に入るかだけに影響する。このサンプリングは確率的であるが、どの行列式が部分空間に入るかだけに影響する。その後のステップは、部分空間内で有効ハミルトニアンを構築し、それを対角化することである。その結果、CI 係数の唯一の誤差は、統計的ノイズではなく、部分空間の切り捨てに起因する。この特徴は、計算基底トモグラフィ(サンプリング)[37]や古典的なシャドウ法[3840]とは対照的で、どちらもショットノイズを直接再構成状態に伝播させる。ノイズのない波動関数は、下流の古典的アルゴリズムにとって特に有益である。例えば、量子-古典ハイブリッド量子モンテカルロ法は[11]、安定な計算を妨げていた加法誤差とは対照的に、安定に実行された。従って、QSCIは高度な相関波動関数法の強力なパートナーとなる。
In this study, we introduce a TCC variant tailored with QSCI (denoted as QSCI-TCC). The QSCI procedure efficiently samples the electronic configurations realized on the quantum device and reconstructs the coherent wave function on a classical computer. This wave function, which captures the strong correlation, is mapped to CC amplitudes; subsequent CC iterations then recover the remaining dynamical correlation. We validate QSCI-TCC on several MR systems, such as the simultaneous bond dissociation in H_(2)O\mathrm{H}_{2} \mathrm{O} and the triple bond dissociation in N_(2)\mathrm{N}_{2}. We further investigate how finite shot measurements affect their accuracy. 本研究では、QSCIを用いて調整されたTCCの変形(QSCI-TCCと呼ぶ)を紹介する。QSCI 法は、量子デバイス上で実現される電子配置を効率的にサンプリングし、古典計算機上でコヒーレント波動関数を再構成する。この波動関数は強い相関を捉え、CC 振幅にマッピングされます。その後のCC 反復により、残りの動的相関が回復されます。われわれは、 H_(2)O\mathrm{H}_{2} \mathrm{O} の同時結合解離や N_(2)\mathrm{N}_{2} の三重結合解離など、いくつかのMR 系でQSCI-TCCを検証した。さらに、有限ショット測定がその精度にどのような影響を与えるかを調べる。
The remainder of this paper is organized as follows. Sec II reviews the TCC and QSCI formalisms and presents an overview of QSCI-TCC. Computational details are given in Sec III. Results are discussed in Sec IV, and Sec V summarizes our contributions and outlines the 本稿の残りの部分は以下のように構成されている。第 II 章では、TCCとQSCIの形式論をレビューし、QSCI-TCCの概要を示す。計算の詳細は第 III 節で述べる。結果は第 IV 節で議論し、第 V 節で我々の貢献を要約し、以下を概説する。
current limitations. 現在の限界
II. COMPUTATIONAL METHODS II.計算方法
A. Review of tailored coupled cluster A.テーラード・カップルド・クラスターのレビュー
TCC [26,27][26,27] factorizes the cluster exponential e^( hat(T))e^{\hat{T}} as TCC [26,27][26,27] はクラスタ指数 e^( hat(T))e^{\hat{T}} を次のように因数分解する。
where |Psi_(0):)\left|\Psi_{0}\right\rangle is a reference Slater determinant. The cluster operator is partitioned into an active-space component hat(T)^("active ")\hat{T}^{\text {active }}, which acts entirely within the active space, and a remainder hat(T)^("rest ")\hat{T}^{\text {rest }}. The splitting of the cluster exponential into static and dynamical parts in TCC is similar to that in the mixed-exponentially generated four (MEG4) method of Nakatsuji et al. [41, 42]; however, the subsequent formalism is substantially different. ここで |Psi_(0):)\left|\Psi_{0}\right\rangle は参照スレーター行列式である。クラスターオペレーターは、活性空間内で完全に作用する活性空間成分 hat(T)^("active ")\hat{T}^{\text {active }} と余り hat(T)^("rest ")\hat{T}^{\text {rest }} に分割される。TCCにおけるクラスターエクスポネンシャルの静的部分と動的部分への分割は、Nakatsujiら[41, 42]のMEG4(mixed-exponentially generated four)法と同様であるが、その後のフォーマリズムは大きく異なる。
The split cluster operators are explicitly written as, スプリットクラスター演算子は、明示的に次のように書かれる、
{:[ hat(T)^("active ")= hat(T)_(1)^("active ")+ hat(T)_(2)^("active ")],[=sum_(i,a)t_(i)^(a) hat(a)_(a)^(†) hat(a)_(i)+sum_(i,j,a,b)t_(ij)^(ab) hat(a)_(a)^(†) hat(a)_(b)^(†) hat(a)_(i) hat(a)_(j)","],[i","j","a","b in" active space "","]:}\begin{aligned}
\hat{T}^{\text {active }}= & \hat{T}_{1}^{\text {active }}+\hat{T}_{2}^{\text {active }} \\
= & \sum_{i, a} t_{i}^{a} \hat{a}_{a}^{\dagger} \hat{a}_{i}+\sum_{i, j, a, b} t_{i j}^{a b} \hat{a}_{a}^{\dagger} \hat{a}_{b}^{\dagger} \hat{a}_{i} \hat{a}_{j}, \\
& i, j, a, b \in \text { active space },
\end{aligned}
and そして
{:[ hat(T)^(rest)= hat(T)_(1)^(rest)+ hat(T)_(2)^(rest)],[=sum_(i,a)t_(i)^(a) hat(a)_(a)^(†) hat(a)_(i)+sum_(i,j,a,b)t_(ij)^(ab) hat(a)_(a)^(†) hat(a)_(b)^(†) hat(a)_(i) hat(a)_(j)","],[{i","j","a","b}⊄" active space "","]:}\begin{aligned}
\hat{T}^{\mathrm{rest}}= & \hat{T}_{1}^{\mathrm{rest}}+\hat{T}_{2}^{\mathrm{rest}} \\
= & \sum_{i, a} t_{i}^{a} \hat{a}_{a}^{\dagger} \hat{a}_{i}+\sum_{i, j, a, b} t_{i j}^{a b} \hat{a}_{a}^{\dagger} \hat{a}_{b}^{\dagger} \hat{a}_{i} \hat{a}_{j}, \\
& \{i, j, a, b\} \not \subset \text { active space },
\end{aligned}
where hat(a)_(a)^(†)\hat{a}_{a}^{\dagger} and hat(a)_(i)\hat{a}_{i} are the creation and annihilation operators for orbitals aa and ii, respectively, and t_(i)^(a),t_(ij)^(ab)t_{i}^{a}, t_{i j}^{a b} are the CC amplitudes. ここで、 hat(a)_(a)^(†)\hat{a}_{a}^{\dagger} と hat(a)_(i)\hat{a}_{i} はそれぞれ軌道 aa と ii の生成・消滅演算子であり、 t_(i)^(a),t_(ij)^(ab)t_{i}^{a}, t_{i j}^{a b} はCC 振幅である。
In conventional CC theory, the amplitudes are obtained by solving the amplitude equations [43]. In TCC, the amplitudes within the active space are fixed by CI coefficients via 従来のCC 理論では、振幅は振幅方程式を解くことによって得られる[43]。TCCでは、アクティブ空間内の振幅は、CI 係数によって固定される。
where hat(C)_(1)\hat{C}_{1} and hat(C)_(2)\hat{C}_{2} are the CI singles and doubles operators. With a constant hat(T)^("active ")\hat{T}^{\text {active }}, the remaining amplitudes are solved in the usual CC manner. ここで、 hat(C)_(1)\hat{C}_{1} と hat(C)_(2)\hat{C}_{2} はCIのシングルスとダブルスの演算子である。定数 hat(T)^("active ")\hat{T}^{\text {active }} で、残りの振幅は通常のCCの方法で解かれます。
In addition to this tailored coupled-cluster singles-anddoubles approach, a perturbative triples ( T ) correction can be introduced [44]. This correction is introduced by setting the active-space amplitudes to zero, thereby preventing the double counting of static correlation. この調整されたクラスターシングルアンドダブルスアプローチに加えて,摂動的トリプル( T )補正を導入することができる[44]。この補正は活性空間の振幅をゼロにすることで導入され、静的相関の二重計数を防ぐことができます。
B. Review of quantum selected configuration interaction B.量子選択配置相互作用のレビュー
1. Basic formulation of QSCI 1.QSCIの基本定式化
The QSCI method [33] samples the quantum state prepared on a quantum computer. The resulting configurations {|Phi_(i):)}\left\{\left|\Phi_{i}\right\rangle\right\} span a subspace of the Fock space, and an effective Hamiltonian QSCI 法[33]は、量子コンピュータ上で準備された量子状態をサンプリングします。得られた配置 {|Phi_(i):)}\left\{\left|\Phi_{i}\right\rangle\right\} はフォック空間の部分空間にまたがり、有効ハミルトニアン
where c_(i)c_{i} is the expansion coefficient of the ii th basis state. RR is the number of retained configurations, which can be truncated based on the sampling frequency. Thus, the quantum state is faithfully reconstructed on a classical computer. ここで、 c_(i)c_{i} は ii 番目の基底状態の展開係数である。 RR は保持される状態の数であり、サンプリング周波数に基づいて切り捨てることができる。このように、量子状態は古典コンピュータ上で忠実に再構成される。
2. Cartesian product of bitstrings 2.ビット列のデカルト積
Finite shot noise can lead to spin-symmetry breaking in the sampled space. To mitigate this, each sampled determinant can be separated into its alpha\alpha - and beta\beta-spin parts, 有限のショットノイズは、サンプリングされた空間におけるスピン対称性の破れにつながる可能性がある。これを緩和するために、サンプリングされた各行列式を alpha\alpha -と beta\beta -スピン部分に分離することができます、
and the Cartesian product can be formed そしてデカルト積は次のようになる。
{| tilde(Phi)_(k):)}={|Phi_(i)^(alpha):)|Phi_(j)^(beta):),^(AA)i,j <= R},\left\{\left|\tilde{\Phi}_{k}\right\rangle\right\}=\left\{\left|\Phi_{i}^{\alpha}\right\rangle\left|\Phi_{j}^{\beta}\right\rangle,{ }^{\forall} i, j \leq R\right\},
thereby enlarging the subspace spanned by {| tilde(Phi)_(i):)}\left\{\left|\tilde{\Phi}_{i}\right\rangle\right\} instead of {|Phi_(i):)}\left\{\left|\Phi_{i}\right\rangle\right\} and increasing the likelihood that all relevant spin eigenfunctions are included [11, 35]. Although this increases the dimension of hat(H)^("eff ")\hat{H}^{\text {eff }}, it often improves accuracy. それによって、 {|Phi_(i):)}\left\{\left|\Phi_{i}\right\rangle\right\} の代わりに {| tilde(Phi)_(i):)}\left\{\left|\tilde{\Phi}_{i}\right\rangle\right\} によってスパンされる部分空間を拡大し、すべての関連するスピン固有関数が含まれる可能性を高めます[11, 35]。これは hat(H)^("eff ")\hat{H}^{\text {eff }} の次元を増加させますが、しばしば精度を向上させます。
To ensure spin adaptation, we treat the union of the determinants sampled for the alpha\alpha - and beta\beta-spin sectors as the configuration set for each spin. Whenever a determinant of the form Eq. (9) is sampled, the subspace is enlarged by adding its spin-swapped partner |Phi_(i^(')):)\left|\Phi_{i^{\prime}}\right\rangle, defined as スピン適応を確実にするために、 alpha\alpha -と beta\beta -スピンセクターのためにサンプリングされた行列式の和を各スピンに対するコンフィギュレーションセットとして扱います。式(9)の行列式がサンプリングされるたびに、部分空間はスピンスワップされたパートナー |Phi_(i^(')):)\left|\Phi_{i^{\prime}}\right\rangle を追加することによって拡大されます。
Because TCC at the singles-and-doubles level includes only up to second-order excitations, mapping CI coefficients to CC amplitudes introduces an error in the QSCITCC energy E_(QSCI-TCC)E_{\mathrm{QSCI}-\mathrm{TCC}}. Two main sources contribute to this error: (i) the neglect of triples or higher excitations and (ii) an incomplete set of determinants in QSCI. To compensate, a corrected energy 一重二重励起レベルのTCCは二次励起までしか含まないので、CI 係数をCC 振幅にマッピングするとQSCITCCのエネルギー E_(QSCI-TCC)E_{\mathrm{QSCI}-\mathrm{TCC}} に誤差が生じる。この誤差には2つの主な原因がある:(i)3 倍以上の励起の無視、(ii)QSCIの不完全な行列式のセット。これを補うために、補正されたエネルギー
is introduced [13]. Here, E_(QC)^("active ")E_{Q C}^{\text {active }} denotes the total energy obtained within the active-space approximation via quantum computing. In this study, we substitute the QSCI analog E_(QSCI)^("active ")E_{\mathrm{QSCI}}^{\text {active }} for E_(QC)^("active ")E_{\mathrm{QC}}^{\text {active }}. The last term, が導入されている[13]。ここで、 E_(QC)^("active ")E_{Q C}^{\text {active }} は、量子計算によって活性空間近似内で得られる全エネルギーを示す。本研究では、QSCIアナログ E_(QSCI)^("active ")E_{\mathrm{QSCI}}^{\text {active }} を E_(QC)^("active ")E_{\mathrm{QC}}^{\text {active }} に置き換える。最後の項は
corresponds to the energy obtained immediately after tailoring, i.e., before the amplitudes outside the active space are optimized. は、テーラリング直後、つまりアクティブスペース外の振幅が最適化される前に得られるエネルギーに相当する。
III. COMPUTATIONAL DETAILS III.計算の詳細
This section summarizes the settings used for the numerical evaluation of QSCI-TCC. このセクションでは、QSCI-TCCの数値評価に使用した設定をまとめる。
We used PySCF 2.2.1 [45, 46] for quantum chemical calculations, including Hartree-Fock (HF), CCSD, CCSD(T), and full configuration interaction (FCI). The basis set used for each calculation is specified individually in the corresponding subsection of Sec. IV. 量子化学計算にはPySCF 2.2.1 [45, 46]を使用し、Hartree-Fock (HF)、CCSD、CCSD(T)、full configuration interaction (FCI)を含む。各計算に使用した基底セットは、Sec.IVの対応するサブセクションに個別に記載されている。
Simulations of quantum computing were performed with Quri-Parts 0.20.3. VQE simulation runs relied on Chemqulacs [47] and used two ansatzes: the disentangled unitary coupled cluster singles and doubles (UCCSD) ansatz [34, 48] and the GateFabric ansatz [49]. Fermionic operators were mapped to qubits using the Jordan-Wigner transformation. 量子コンピューティングのシミュレーションは、Quri-Parts 0.20.3で行った。VQEシミュレーションは、Chemqulacs [47]を使用し、2つのアサッツを使用しました。核融合作用素はJordan-Wigner 変換を用いて量子ビットにマッピングした。
We performed the 10^(7)10^{7} shot measurements for sampling quantum states in the investigation of energy curves in Sec. IV A. The influence of shot count is analyzed in detail in Sec. IV B. Effective Hamiltonians generated by QSCI were diagonalized using the kernel_fixed_space function in PySCF. QSCIで生成した有効ハミルトニアンはPySCFのkernel_fixed_space 関数を用いて対角化した。
IV. RESULTS AND DISCUSSION IV.結果と考察
A. Bond dissociation A.結合解離
We benchmark QSCI-TCC on two prototypical strongcorrelation problems: the simultaneous dissociation of both O-H\mathrm{O}-\mathrm{H} bonds in H_(2)O\mathrm{H}_{2} \mathrm{O} and the triple bond dissociation in N_(2)\mathrm{N}_{2}. These systems are widely used to test methods designed for MR situations. H_(2)O\mathrm{H}_{2} \mathrm{O} における O-H\mathrm{O}-\mathrm{H} 両結合の同時解離と N_(2)\mathrm{N}_{2} における三重結合の解離である。これらの系はMRの状況を想定した手法のテストに広く使われている。
1. Simultaneous dissociation of the OH bonds in H_(2)O\mathrm{H}_{2} \mathrm{O} 1. H_(2)O\mathrm{H}_{2} \mathrm{O} のOH 結合が同時に解離する。
Figure 2 (a) shows the potential energy curves for the simultaneous stretching of the two O-H\mathrm{O}-\mathrm{H} bonds. Starting from the HF baseline, static correlation within the (8e, 6o) active space is introduced and captured. Both active-space QSCI and VQE lower the energy relative to HF, and the resulting gap widens as the bond lengths rr become larger. 図 2(a)は、2つの O-H\mathrm{O}-\mathrm{H} 結合の同時伸長に対するポテンシャルエネルギー曲線を示している。HFベースラインから出発して、(8e, 6o)活性空間内の静的相関が導入され、捕捉される。活性空間 QSCIとVQEの両方がHFに比べてエネルギーを下げ、その結果、結合長 rr が大きくなるにつれてギャップが広がる。
The subsequent TCC calculations incorporate dynamical correlation from the remaining orbitals, further lowering the energy across the entire region. その後のTCC 計算では、残りの軌道からの動的相関が組み込まれ、領域全体のエネルギーがさらに低下する。
We also plot CCSD, CCSD(T), and exact FCI energies for comparison. Their curves are nearly indistinguishable from the QSCI-TCC curves; however, CCSD(T)\operatorname{CCSD}(\mathrm{T}) begins to depart from the others at approximately r=2.0"Å"r=2.0 \AAÅ. 比較のためにCCSD、CCSD(T)、正確なFCIエネルギーもプロットした。これらの曲線はQSCI-TCCの曲線とほとんど区別がつかない。しかし、 CCSD(T)\operatorname{CCSD}(\mathrm{T}) は約 r=2.0"Å"r=2.0 \AAÅ で他の曲線から離れ始める。
2. Triple-bond dissociation in N_(2)N_{2} 2. N_(2)N_{2} における三重結合解離。
Fig. 3 (a) presents the potential energy curves for the dissociation of the triple bond in N_(2)\mathrm{N}_{2}. Similar to H_(2)O\mathrm{H}_{2} \mathrm{O}, static correlation recovered within the ( 6e,6o6 \mathrm{e}, 6 \mathrm{o} ) active space is essential for a qualitatively correct description of the dissociation region. Our VQE calculation, which employs a three-layer GateFabric ansatz, fails to produce a smooth energy curve; nevertheless, QSCI delivers a stable result even when the underlying VQE wave function is imperfect. This robustness can be attributed to the fact that the effective Hamiltonian is built from computational basis states that belong to the exact wave function, regardless of the quality of the trial state. All methods that include correlation over the full set of orbitals - such as QSCI-TCC and its variants, CCSD(T), and the semistochastic heat-bath CI (SHCI) reference - converge to an essentially common curve on the scale of Fig. 3 (a). The only exception is that CCSD underestimates the energy. The referential SHCI data are taken from Ref. [50]. 図 3(a)は N_(2)\mathrm{N}_{2} における三重結合の解離のポテンシャルエネルギー曲線である。 H_(2)O\mathrm{H}_{2} \mathrm{O} と同様に、( 6e,6o6 \mathrm{e}, 6 \mathrm{o} )活性空間内で回収される静的相関は、解離領域を定性的に正しく記述するために不可欠である。しかし、QSCIはVQEの波動関数が不完全であっても安定した結果を得ることができる。この頑健性は、試行状態の質に関係なく、正確な波動関数に属する計算基底状態から有効ハミルトニアンを構築していることに起因しています。QSCI-TCCとその変種、CCSD(T)、セミストカスティック・ヒートバスCI(SHCI)など、軌道の全セットの相関を含むすべての方法は、図 3(a)のスケールで本質的に共通の曲線に収束します。唯一の例外は、CCSDがエネルギーを過小評価していることである。参考にしたSHCIデータは文献[50]から引用しました。[50].
For further evaluation of accuracy, Fig. 3 (b) shows the energy deviation from SHCI. The CCSD error grows monotonically over the scanned range. In the region from 1.8 to 2.0"Å"2.0 \AAÅ, the CCSD(T)\operatorname{CCSD}(\mathrm{T}) error also starts to increase, reflecting the well-known breakdown of CCSD(T)\operatorname{CCSD}(\mathrm{T}) for N_(2)\mathrm{N}_{2} dissociation. By contrast, the QSCI-TCC and QSCITCC(T) error curves remain nearly flat, mirroring the behavior already observed for H_(2)O\mathrm{H}_{2} \mathrm{O}. This stability indicates that the CC amplitudes not only describe dynamical correlation but also successfully incorporate static correlation by embedding the QSCI wave function. さらに精度を評価するために、図 3(b)にSHCIからのエネルギー偏差を示す。CCSDの誤差はスキャンされた範囲で単調に増加している。1.8から 2.0"Å"2.0 \AAÅ の領域では、 CCSD(T)\operatorname{CCSD}(\mathrm{T}) 誤差も増加し始め、これは N_(2)\mathrm{N}_{2} 解離に対する CCSD(T)\operatorname{CCSD}(\mathrm{T}) のよく知られた破綻を反映している。対照的に、QSCI-TCCとQSCITCC(T)の誤差曲線はほぼ平坦で、 H_(2)O\mathrm{H}_{2} \mathrm{O} で観測された振る舞いを反映している。この安定性は、CC 振幅が動的相関を記述するだけでなく、QSCI 波動関数を埋め込むことによって静的相関をうまく取り込んでいることを示している。
FIG. 2: (a) Energy curve and (b) energy difference compared to FCI for simultaneous O-H\mathrm{O}-\mathrm{H} bond dissociation in H_(2)O\mathrm{H}_{2} \mathrm{O} at the cc-pVDZ basis set level. The active space is (8e, 6o), and the ansatz for VQE is the UCCSD ansatz. 図 2: (a) cc-pVDZ 基底セットレベルでの H_(2)O\mathrm{H}_{2} \mathrm{O} における O-H\mathrm{O}-\mathrm{H} 結合同時解離のエネルギー曲線と(b)FCIと比較したエネルギー差。活性空間は(8e, 6o)で、VQEのアサッツはUCCSDアサッツである。
FIG. 3: (a) Energy curve and (b) energy difference compared to SHCI for N_(2)\mathrm{N}_{2} dissociation at the cc-pVTZ basis set level. The active space is (6e,6o). The GateFabric ansatz with three layers is used for VQE. 図 3:(a)cc-PVTZ 基底セットレベルでの N_(2)\mathrm{N}_{2} 解離のエネルギー曲線と(b)SHCIとのエネルギー差。活性空間は(6e,6o)。VQEには3 層のGateFabric 解析法を用いた。
In the cases of less than 10^(5)10^{5} shots, the observed energy distribution is discrete rather than continuous, reflecting the binary occurrence of key determinants in the sampled set. The energies can be lower than the converged value because CC theory is not variational in nature. 10^(5)10^{5} ショット未満の場合、観測されたエネルギー分布は連続的ではなく離散的であり、これはサンプリングされた集合における主要な決定基の二者択一的な出現を反映している。CC 理論は本質的に変分ではないため、エネルギーは収束値よりも低くなる可能性があります。
From the above observations, we conclude that the 10^(5)10^{5} shot measurements achieve an error below 1kcal//mol1 \mathrm{kcal} / \mathrm{mol} ( ∼1.6 xx10^(-3)E_(h)\sim 1.6 \times 10^{-3} E_{\mathrm{h}} ) in this setting. No apparent difference is observed between the union and non-union treatments in this study; in both cases, the number of measurements required for convergence remains of the same order. 以上の観察から、 10^(5)10^{5} ショット測定は、この設定において 1kcal//mol1 \mathrm{kcal} / \mathrm{mol} ( ∼1.6 xx10^(-3)E_(h)\sim 1.6 \times 10^{-3} E_{\mathrm{h}} )以下の誤差を達成すると結論づけられる。この研究では、ユニオン処理と非ユニオン処理の間に明らかな違いは見られず、どちらの場合も収束に必要な測定回数は同じオーダーのままである。
The present requirement of 10^(5)10^{5} shots outperforms earlier quantum-classical hybrid TCC studies. Computational basis tomography required 3xx10^(7)3 \times 10^{7} shots to achieve the standard deviation of 0.61 xx10^(-3)E_(h)0.61 \times 10^{-3} E_{\mathrm{h}} in the r=5.355r=5.355 Bohr ( ∼2.8"Å"\sim 2.8 \AAÅ ) case at the cc-pVDZ level [13]. Although this previous study was conducted under a challenging bond-distance setting, the relative phases need to be determined for computational basis tomography [37]; therefore, QSCI is considered superior. In addition, at the cc-pVDZ level, the estimated result of the matchgate shadows tomography demonstrates that 1.0 xx10^(5)1.0 \times 10^{5} and 2.2 xx10^(6)2.2 \times 10^{6} shots were required to obtain chemical precision at r=1.1r=1.1 and 2.2"Å"2.2 \AAÅ, respectively [14], where the latter bond length is more strongly correlated. These comparisons underscore the shot efficiency of QSCI-TCC. 現在の 10^(5)10^{5} ショットの必要性は、以前の量子-古典ハイブリッドTCC 研究よりも優れている。計算基底トモグラフィーでは、cc-pVDZレベルで r=5.355r=5.355 Bohr ( ∼2.8"Å"\sim 2.8 \AAÅ )の場合の標準偏差 0.61 xx10^(-3)E_(h)0.61 \times 10^{-3} E_{\mathrm{h}} を達成するために 3xx10^(7)3 \times 10^{7} ショットが必要であった[13]。この先行研究は難しい結合距離の設定の下で行われたが、計算基底トモグラフィー[37]のためには相対相を決定する必要があるため、QSCIが優れていると考えられる。さらに、cc-pVDZレベルでは、マッチゲート・シャドー・トモグラフィーによる推定結果は、 r=1.1r=1.1 と 2.2"Å"2.2 \AAÅ で化学的精度を得るために、それぞれ 1.0 xx10^(5)1.0 \times 10^{5} と 2.2 xx10^(6)2.2 \times 10^{6} ショットが必要であることを示している[14]。これらの比較は、QSCI-TCCのショット効率を強調している。
V. CONCLUSIONS V.結論
We introduced QSCI-TCC, a quantum-classical hybrid method that integrates the QSCI wave function into the TCC framework. QSCI provides a scalable, shot-efficient way to reconstruct a quantum-prepared state on a classical computer. The wave function reconstructed by QSCI contains no statistical error, only model error, providing a stable state for classical post-processing. Embedding the resulting CI coefficients as fixed active space amplitudes, the subsequent CC calculation supplies the missing dynamical correlation, yielding a balanced treatment of static (strong) and dynamical (weak) correlation effects. 我々は、QSCI 波動関数をTCCの枠組みに統合した量子-古典ハイブリッド法であるQSCI-TCCを導入した。QSCIは、古典計算機上で量子準備状態を再構成するスケーラブルでショット効率の高い方法を提供する。QSCIによって再構築された波動関数は統計誤差を含まず、モデル誤差のみを含むため、古典的な後処理に対して安定した状態を提供します。得られたCI 係数を固定された活性空間振幅として埋め込むことで、その後のCC 計算で欠落していた動的相関が供給され、静的相関効果(強相関)と動的相関効果(弱相関)をバランスよく取り扱うことができる。
Benchmark calculations for simultaneous O-H\mathrm{O}-\mathrm{H} bond dissociation in H_(2)O\mathrm{H}_{2} \mathrm{O} and triple-bond dissociation in N_(2)\mathrm{N}_{2} show that, even in regions of the dissociation curves where CCSD or CCSD(T) begins to deteriorate, QSCITCC and QSCI-TCC(T) continue to deliver accurate results. Thus, the static correlation captured by quantum computation, although a simulator has been used in this H_(2)O\mathrm{H}_{2} \mathrm{O} の O-H\mathrm{O}-\mathrm{H} 結合解離と N_(2)\mathrm{N}_{2} の三重結合解離を同時に行うベンチマーク計算では、解離曲線のCCSDやCCSD(T)が劣化し始める領域でも、QSCITCCとQSCI-TCC(T)は正確な結果を出し続けている。このように、量子計算によって捉えられる静的相関は、シミュレーターを用いたとはいえ、この
study, is successfully embedded within the coupled cluster description. の研究は、連成クラスターの記述の中にうまく組み込まれている。
Efficient quantum state tomography-using QSCI as well as other tomography methods-remains a key research topic. Most recently, Lenihan et al. have introduced a related study on another shadows-based protocol combined with a surrogate CC model [51]. QSCIや他のトモグラフィ手法を用いた効率的な量子状態トモグラフィは、依然として重要な研究テーマである。最近では、LenihanらがサロゲートCCモデルと組み合わせた別のシャドウ・ベース・プロトコルに関する関連研究を紹介している[51]。
It is worth mentioning a critical issue with QSCI and the proposed improvements in recent years. In strongly correlated systems, the number of Slater de- QSCIと近年提案されている改良について、重大な問題に言及する価値がある。強い相関を持つシステムでは、スレーターデ
terminants with non-trivial coefficients can grow exponentially with system size. To mitigate this, enhanced sampling schemes that utilize real-time evolution have been proposed [52-54]. Besides, Reinholdt et al. have compared QSCI with a conventional selected CI intensively [55]. Improving sampling efficiency, therefore, remains an important avenue for future work. 非自明な係数を持つターミナントは、システムサイズとともに指数関数的に増大する可能性がある。これを緩和するために、リアルタイム進化を利用した拡張サンプリング方式が提案されている[52-54]。また、Reinholdtらは、QSCIと従来の選択 CIを集中的に比較している[55]。従って、サンプリング効率の改善は、今後の重要な課題である。
ACKOWLEDGEMENTS 謝辞
The authors are grateful to Takuma Murokoshi for his technical support. This project was supported by funding from the MEXT Quantum Leap Flagship Program (MEXTQLEAP) through Grant No. JPMXS0120319794. This study was conducted as part of a joint research in the Quantum Software Research Hub (Grant No. JPMJPF2014), and further supported by TOPPAN Digital Inc. 室越琢磨氏の技術支援に感謝する。本プロジェクトは、文部科学省量子リープフラッグシップ研究プログラム(MEXTQLEAP)の助成金(助成番号:JPMXS0120319794)の支援を受けて実施した。本研究は、量子ソフトウェア研究ハブにおける共同研究(助成金番号 JPMJPF2014)の一環として実施され、さらに株式会社トッパンデジタルによる支援を受けた。
[1] D. S. Abrams and S. Lloyd, Simulation of many-body fermi systems on a universal quantum computer, Phys. Rev. Lett. 79, 2586 (1997). [1] D. S. Abrams and S. Lloyd, Simulation of many-body fermi systems on a universal quantum computer, Phys.79, 2586 (1997).
[2] D. S. Abrams and S. Lloyd, Quantum algorithm providing exponential speed increase for finding eigenvalues and eigenvectors, Phys. Rev. Lett. 83, 5162 (1999). [2] D. S. Abrams and S. Lloyd, Quantum algorithm providing exponential speed increase for finding eigenvalues and eigenvectors, Phys.83, 5162 (1999).
[3] A. Aspuru-Guzik, A. D. Dutoi, P. J. Love, and M. HeadGordon, Simulated quantum computation of molecular energies, Science 309, 1704 (2005).
[4] M. Reiher, N. Wiebe, K. M. Svore, D. Wecker, and M. Troyer, Elucidating reaction mechanisms on quantum computers, Proc. Natl. Acad. Sci. USA 114, 7555 (2017). [4] M. Reiher, N. Wiebe, K. M. Svore, D. Wecker, and M. Troyer, Elucidating reaction mechanisms on quantum computers, Proc.Natl. Acad.Sci. USA 114, 7555 (2017).
[5] D. W. Berry, C. Gidney, M. Motta, J. R. McClean, and R. Babbush, Qubitization of Arbitrary Basis Quantum Chemistry Leveraging Sparsity and Low Rank Factorization, Quantum 3, 208 (2019).
[6] V. von Burg, G. H. Low, T. Häner, D. S. Steiger, M. Reiher, M. Roetteler, and M. Troyer, Quantum computing enhanced computational catalysis, Phys. Rev. Res. 3, 033055 (2021).
[7] J. Lee, D. W. Berry, C. Gidney, W. J. Huggins, J. R. McClean, N. Wiebe, and R. Babbush, Even more efficient quantum computations of chemistry through tensor hypercontraction, PRX Quantum 2, 030305 (2021).
[8] D. Rocca, C. L. Cortes, J. F. Gonthier, P. J. Ollitrault, R. M. Parrish, G.-L. Anselmetti, M. Degroote, N. Moll, R. Santagati, and M. Streif, Reducing the Runtime of Fault-Tolerant Quantum Simulations in Chemistry through Symmetry-Compressed Double Factorization, J. Chem. Theory Comput. 20, 4639 (2024). [8] D. Rocca, C. L. Cortes, J. F. Gonthier, P. J. Ollitrault, R. M. Parrish, G.-L. Anselmetti, M. Degroote, N. Moll, R. Santagati, and M. Streif, Reducing Runtime of Fault-Tolerant Quantum Simulations in Chemistry through Symmetry-Compressed Double Factorization, J. Chem.Theory Comput.20, 4639 (2024).
[9] G. Hao Low, R. King, D. W. Berry, Q. Han, A. E. DePrince, III, A. White, R. Babbush, R. D. Somma, and N. C. Rubin, Fast quantum simulation of electronic structure by spectrum amplification, arXiv e-prints , arXiv:2502.15882 (2025), arXiv:2502.15882 [quant-ph].
[10] W. J. Huggins, B. A. O’Gorman, N. C. Rubin, D. R. Reichman, R. Babbush, and J. Lee, Unbiasing fermionic quantum monte carlo with a quantum computer, Nature 603, 416 (2022). [10] W. J. Huggins, B. A. O'Gorman, N. C. Rubin, D. R. Reichman, R. Babbush, and J. Lee, Unbiasing fermionic quantum monte carlo with a quantum computer, Nature 603, 416 (2022).
[11] Y. Yoshida, L. Erhart, T. Murokoshi, R. Nakagawa, C. Mori, T. Miyanaga, T. Mori, and W. Mizukami, Auxiliary-field quantum Monte Carlo method with quantum selected configuration interaction, arXiv e-prints , arXiv:2502.21081 (2025), arXiv:2502.21081 [quant-ph]. [11] Yoshida, L. Erhart, T. Murokoshi, R. Nakagawa, C. Mori, T. Miyanaga, T. Mori, and W. Mizukami, Auxiliary-field quantum Monte Carlo method with quantum selected configuration interaction, arXiv e-prints , arXiv:2502.21081 (2025), arXiv:2502.21081 [quant-ph].
[12] D. Danilov, J. Robledo-Moreno, K. J. Sung, M. Motta, and J. Shee, Enhancing the accuracy and efficiency of sample-based quantum diagonalization with phaseless auxiliary-field quantum Monte Carlo, arXiv e-prints , arXiv:2503.05967 (2025), arXiv:2503.05967 [quant-ph]. (12)D.Danilov、J.Robledo-Moreno、K.J.Sung、M.Motta、J.Shee, Enhancing accuracy and efficiency of sample-based quantum diagonalization with phaseless auxiliary-field quantum Monte Carlo, arXiv e-prints , arXiv:2503.05967 (2025), arXiv:2503.05967 [quant-ph].
[13] L. Erhart, Y. Yoshida, V. Khinevich, and W. Mizukami, Coupled cluster method tailored with quantum computing, Phys. Rev. Res. 6, 023230 (2024).
[14] M. Scheurer, G.-L. R. Anselmetti, O. Oumarou, C. Gogolin, and N. C. Rubin, Tailored and Externally Corrected Coupled Cluster with Quantum Inputs, J. Chem. Theory Comput. 20, 5068 (2024). [14] M. Scheurer, G.-L. R. Anselmetti, O. Oumarou, C. Gogolin, and N. C. Rubin, Tailored and Externally Corrected Coupled Cluster with Quantum Input, J. Chem.Theory Comput.20, 5068 (2024).
[15] R. J. Bartlett and M. Musial, Coupled-cluster theory in quantum chemistry, Rev. Mod. Phys. 79, 291 (2007). [15] R. J. Bartlett and M. Musial, Coupled-cluster theory in quantum chemistry, Rev. Mod.Phys. 79, 291 (2007).
[16] R. J. Bartlett, Perspective on coupled-cluster theory. the evolution toward simplicity in quantum chemistry, Phys. Chem. Chem. Phys. 26, 8013 (2024). [16] R. J. Bartlett, Perspective on coupled-cluster theory. The evolution towards simplicity in quantum chemistry, Phys.Chem.Phys. 26, 8013 (2024).
[17] C. Riplinger and F. Neese, An efficient and near linear scaling pair natural orbital based local coupled cluster method, J. Chem. Phys. 138, 034106 (2013). [17] C. Riplinger and F. Neese, An efficient and near linear scaling pair natural orbital based local coupled cluster method, J. Chem.Phys. 138, 034106 (2013).
[18] C. Riplinger, B. Sandhoefer, A. Hansen, and F. Neese, Natural triple excitations in local coupled cluster calculations with pair natural orbitals, J. Chem. Phys. 139, 134101 (2013). [18] C. Riplinger, B. Sandhoefer, A. Hansen, and F. Neese, Natural triple excitations in local coupled cluster calculations with pair natural orbitals, J. Chem.Phys. 139, 134101 (2013).
[19] C. Riplinger, P. Pinski, U. Becker, E. F. Valeev, and F. Neese, Sparse maps-A systematic infrastructure for reduced-scaling electronic structure methods. II. Linear scaling domain based pair natural orbital coupled cluster theory, J. Chem. Phys. 144, 024109 (2016). [19] C. Riplinger, P. Pinski, U. Becker, E. F. Valeev, and F. Neese, Sparse maps-A systematic infrastructure for reduced-scaling electronic structure methods.II.Phys. 144, 024109 (2016).
[20] M. Saitow, U. Becker, C. Riplinger, E. F. Valeev, and F. Neese, A new near-linear scaling, efficient and accurate, open-shell domain-based local pair natural orbital coupled cluster singles and doubles theory, J. Chem. Phys. 146, 164105 (2017). [20] M. Saitow, U. Becker, C. Riplinger, E. F. Valeev, and F. Neese, A new near-linear scaling, efficient and accurate, open-shell domain-based local pair natural orbital coupled cluster singles and doubles theory, J. Chem.Phys. 146, 164105 (2017).
[21] Y. Guo, C. Riplinger, U. Becker, D. G. Liakos, Y. Minenkov, L. Cavallo, and F. Neese, Communication: An improved linear scaling perturbative triples correction for the domain based local pair-natural orbital based singles and doubles coupled cluster method [DLPNO-CCSD(T)], J. Chem. Phys. 148, 011101 (2018). [21] Y. Guo, C. Riplinger, U. Becker, D. G. Liakos, Y. Minenkov, L. Cavallo, and F. Neese, Communication:(1)局所対自然軌道に基づくシングル・ダブルス結合クラスター法[DLPNO-CCSD(T)]の改良された線形スケーリング摂動トリプル補正, J. Chem.Phys. 148, 011101 (2018). [22][22] D. G. Liakos, Y. Guo, and F. Neese, Comprehensive Benchmark Results for the Domain Based Local Pair Natural Orbital Coupled Cluster Method (DLPNOCCSD(T)) for Closed- and Open-Shell Systems, J. Phys. Chem. A 124, 90 (2020). [22][22] D. G. Liakos, Y. Guo, and F. Neese, Comprehensive Benchmark Results for the Domain Based Local Pair Natural Orbital Coupled Cluster Method (DLPNOCCSD(T)) for Closed- and Open-Shell Systems, J. Phys. Chem.A 124, 90 (2020).
[23] Y. Guo, C. Riplinger, D. G. Liakos, U. Becker, M. Saitow, and F. Neese, Linear scaling perturbative triples correction approximations for open-shell domain-based local pair natural orbital coupled cluster singles and doubles theory [DLPNO-CCSD(T/T)], J. Chem. Phys. 152, 024116 (2020). (23)Y.Guo,C.Riplinger,D.G.Liakos,U.Becker,M.Saitow,F.Neese,Linear scaling perturbative triples correction approximations for open-shell domain-based local pair natural orbital coupled cluster singles and doubles theory [DLPNO-CCSD(T/T)], J. Chem.Phys. 152, 024116 (2020).
[24] D. M. Wilkins, A. Grisafi, Y. Yang, K. U. Lao, R. A. DiStasio, and M. Ceriotti, Accurate molecular polarizabilities with coupled cluster theory and machine learning, Proc. Natl. Acad. Sci. USA 116, 3401 (2019). [24] D. M. Wilkins, A. Grisafi, Y. Yang, K. U. Lao, R. A. DiStasio, and M. Ceriotti, Accurate molecular polarizabilities with coupled cluster theory and machine learning, Proc.Natl. Acad.Sci. USA 116, 3401 (2019).
[25] J. S. Smith, B. T. Nebgen, R. Zubatyuk, N. Lubbers, C. Devereux, K. Barros, S. Tretiak, O. Isayev, and A. E. Roitberg, Approaching coupled cluster accuracy with a general-purpose neural network potential through transfer learning, Nat. Commun. 10, 2903 (2019). [25] J. S. Smith, B. T. Nebgen, R. Zubatyuk, N. Lubbers, C. Devereux, K. Barros, S. Tretiak, O. Isayev, and A. E. Roitberg, Approaching coupled cluster accuracy with general-purpose neural network potential through transfer learning, Nat.Commun.10, 2903 (2019).
[26] T. Kinoshita, O. Hino, and R. J. Bartlett, Coupledcluster method tailored by configuration interaction, J. Chem. Phys. 123, 074106 (2005). このような場合、[26] 木下、日野、R. J. Bartlett, Coupledcluster method tailored by configuration interaction, J. Chem.Phys. 123, 074106 (2005).
[27] O. Hino, T. Kinoshita, G. K.-L. Chan, and R. J. Bartlett, Tailored coupled cluster singles and doubles method applied to calculations on molecular structure and harmonic vibrational frequencies of ozone, J. Chem. Phys. 124, 114311 (2006). [このため、オゾンの分子構造と調和振動周波数の計算には、クラスター一重・二重結合法を適用した。Phys. 124, 114311 (2006).
[28] F. M. Faulstich, M. Máté, A. Laestadius, M. A. Csirik, L. Veis, A. Antalik, J. Brabec, R. Schneider, J. Pittner, S. Kvaal, and Ö. Legeza, Numerical and Theoretical Aspects of the DMRG-TCC Method Exemplified by the Nitrogen Dimer, J. Chem. Theory Comput. 15, 2206 (2019). [28] F. M. Faulstich, M. Máté, A. Laestadius, M. A. Csirik, L. Veis, A. Antalik, J. Brabec, R. Schneider, J. Pittner, S. Kvaal, and Ö.Legeza, Nitrogen Dimerに代表されるDMRG-TCC 法の数値的および理論的側面, J. Chem.Theory Comput.15, 2206 (2019).
[29] J. Višňák, J. Brandejs, M. Máté, L. Visscher, Ö. Legeza, and J. Pittner, DMRG-Tailored Coupled Cluster Method in the 4c-Relativistic Domain: General Implementation and Application to the NUHFI and NUF3 Molecules, J. Chem. Theory Comput. 20, 8862 (2024). [29] J. Višňák, J. Brandejs, M. Máté, L. Visscher, Ö.Legeza, and J. Pittner, DMRG-Tailored Coupled Cluster Method in the 4c-Relativistic Domain:一般的な実装とNUHFIおよびNUF3 分子への応用、J. Chem.Theory Comput.20, 8862 (2024).
[30] E. Vitale, A. Alavi, and D. Kats, FCIQMC-Tailored Distinguishable Cluster Approach, J. Chem. Theory Comput. 16, 5621 (2020). [30] E. Vitale, A. Alavi, and D. Kats, FCIQMC-Tailored Distinguishable Cluster Approach, J. Chem.Theory Comput.16, 5621 (2020).
[31] E. Vitale, G. Li Manni, A. Alavi, and D. Kats, FCIQMCTailored Distinguishable Cluster Approach: Open-Shell Systems, J. Chem. Theory Comput. 18, 3427 (2022). [31] E. Vitale, G. Li Manni, A. Alavi, and D. Kats, FCIQMCTailored Distinguishable Cluster Approach:Open-Shell Systems, J. Chem.Theory Comput.18, 3427 (2022).
[32] S. Nishio, Y. Oba, and Y. Kurashige, Statistical errors in reduced density matrices sampled from quantum circuit simulation and the impact on multireference perturbation theory, Phys. Chem. Chem. Phys. 25, 30525 (2023). (32)西尾慎一郎, 大場義人, 倉重義人, 量子回路シミュレーションからサンプリングした縮小密度行列の統計誤差と多参照摂動論への影響, Phys.Chem.Phys. 25, 30525 (2023).
[33] K. Kanno, M. Kohda, R. Imai, S. Koh, K. Mitarai, W. Mizukami, and Y. O. Nakagawa, QuantumSelected Configuration Interaction: classical diagonalization of Hamiltonians in subspaces selected by quantum computers, arXiv e-prints , arXiv:2302.11320 (2023), arXiv:2302.11320 [quant-ph].
[34] A. Peruzzo, J. McClean, P. Shadbolt, M.-H. Yung, X.-Q. Zhou, P. J. Love, A. Aspuru-Guzik, and J. L. O’Brien, [34] 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, Nat. Commun. 5, 4213 (2014). フォトニック量子プロセッサーによる変分固有値ソルバー, Nat.Commun.5, 4213 (2014).
[35] J. Robledo-Moreno, M. Motta, H. Haas, A. JavadiAbhari, P. Jurcevic, W. Kirby, S. Martiel, K. Sharma, S. Sharma, T. Shirakawa, I. Sitdikov, R.-Y. Sun, K. J. Sung, M. Takita, M. C. Tran, S. Yunoki, and A. Mezzacapo, Chemistry Beyond Exact Solutions on a Quantum-Centric Supercomputer, arXiv e-prints , arXiv:2405.05068 (2024), arXiv:2405.05068 [quant-ph]. [35] J. Robledo-Moreno, M. Motta, H. Haas, A. JavadiAbhari, P. Jurcevic, W. Kirby, S. Martiel, K. Sharma, S. Sharma, T. Shirakawa, I. Sitdikov, R.-Y. Sun, K. J. Sung, M. Takita, M. C. Tran, S. Yunoki, and A. Mezzacapo, Chemistry Beyond Exact Solutions on Quantumentric Supercomputer.
[36] A. Hafid, H. Iwakiri, K. Tsubouchi, N. Yoshioka, and M. Kohda, Hardness of classically sampling quantum chemistry circuits, arXiv e-prints , arXiv:2504.12893 (2025), arXiv:2504.12893 [quant-ph].
[37] M. Kohda, R. Imai, K. Kanno, K. Mitarai, W. Mizukami, and Y. O. Nakagawa, Quantum expectation-value estimation by computational basis sampling, Phys. Rev. Res. 4, 033173 (2022). (37)今井亮太, 菅野和彦, 御手洗邦彦, 水上和彦, 中川洋一, 計算基底サンプリングによる量子期待値推定, Phys. Rev. Res. 4, 033173 (2022).
[38] H.-Y. Huang, R. Kueng, and J. Preskill, Predicting many properties of a quantum system from very few measurements, Nat. Phys. 16, 1050 (2020). [38] H.-Y.Huang, R. Kueng, and J. Preskill, Predicting many properties of a quantum system from very few measurements, Nat.Phys. 16, 1050 (2020).
[39] K. Wan, W. J. Huggins, J. Lee, and R. Babbush, Matchgate shadows for fermionic quantum simulation, Commun. Math. Phys. 404, 629 (2023). [39] K. Wan, W. J. Huggins, J. Lee, and R. Babbush, Matchgate shadows for fermionic quantum simulation, Commun.Math.Phys. 404, 629 (2023).
[40] V. Heyraud, H. Chomet, and J. Tilly, Unified framework for matchgate classical shadows, npj Quantum Inf. 11, 65 (2025). [40] V. Heyraud, H. Chomet, and J. Tilly, Unified framework for matchgate classical shadows, npj Quantum Inf.11, 65 (2025).
[41] H. Nakatsuji, Exponentially generated wave functions, J. Chem. Phys. 83, 5743 (1985). [41] 中辻秀夫,指数関数的に生成される波動関数,J. Chem.Phys. 83, 5743 (1985).
[42] H. Nakatsuji, Mixed-exponentially generated wave function method for ground, excited, ionized, and electron attached states of a molecule, J. Chem. Phys. 95, 4296 (1991). [42] 中辻秀之, 分子の基底状態,励起状態,電離状態,電子付着状態に対する混合指数発生波動関数法, J. Chem.Phys. 95, 4296 (1991).
[43] T. D. Crawford and H. F. Schaefer III, An introduction to coupled cluster theory for computational chemists, in Reviews in Computational Chemistry (John Wiley & Sons, Ltd, 2000) pp. 33-136. [44][44] D. I. Lyakh, V. F. Lotrich, and R. J. Bartlett, The ‘tailored’ CCSD(T)\operatorname{CCSD}(\mathrm{T}) description of the automerization of cyclobutadiene, Chem. Phys. Lett. 501, 166 (2011). [44][44] D. I.Lyakh, V. F. Lotrich, and R. J. Bartlett, The 'tailored' CCSD(T)\operatorname{CCSD}(\mathrm{T}) description of the automerization of cyclobutadiene, Chem.Phys. Lett.501, 166 (2011).
[45] Q. Sun, X. Zhang, S. Banerjee, P. Bao, M. Barbry, N. S. Blunt, N. A. Bogdanov, G. H. Booth, J. Chen, Z.-H. Cui, J. J. Eriksen, Y. Gao, S. Guo, J. Hermann, M. R. Hermes, K. Koh, P. Koval, S. Lehtola, Z. Li, J. Liu, N. Mardirossian, J. D. McClain, M. Motta, B. Mussard, H. Q. Pham, A. Pulkin, W. Purwanto, P. J. Robinson, E. Ronca, E. R. Sayfutyarova, M. Scheurer, H. F. Schurkus, J. E. T. Smith, C. Sun, S.-N. Sun, S. Upadhyay, L. K. Wagner, X. Wang, A. White, J. D. Whitfield, M. J. Williamson, S. Wouters, J. Yang, J. M. Yu, T. Zhu, T. C. Berkelbach, S. Sharma, A. Y. Sokolov, and G. K.L. Chan, Recent developments in the PySCF program package, J. Chem. Phys. 153, 024109 (2020). [45] Q. Sun, X. Zhang, S. Banerjee, P. Bao, M. Barbry, N. S. Blunt, N. A. Bogdanov, G. H. Booth, J. Chen, Z.-H. Cui, J. J. Eriksen, Y. Gao, S. Guo, J. Hermann, M. R. Hermes, Koh, P. Koval, S. Le.Cui、J. J. Eriksen、Y. Gao、S. Guo、J. Hermann、M. R. Hermes、Koh、P. Koval、S. Lehtola、Z. Li、J. Liu、N. Mardirossian、J. D. McClain、M. Motta、B.Mussard, H. Q. Pham, A. Pulkin, W. Purwanto, P. J. Robinson, E. Ronca, E. R. Sayfutyarova, M. Scheurer, H. F. Schurkus, J. E. T. Smith, C. Sun, S.-.N.このように、このような研究開発には、様々な工夫が凝らされています。Phys. 153, 024109 (2020).
[46] 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.-L. Chan, Pyscf: the python-based simulations of chemistry framework, WIREs Comput. Mol. Sci. 8, e1340 (2018). このような場合、[46] 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.-L. Chan, Pyscf: the python-based simulations of chemistry framework, WIREs Comput.Mol.Sci. 8, e1340 (2018).
[47] Chemqulacs, https://wmizukami.github.io/ chemqulacs (2023).
[48] W. Kutzelnigg, Quantum chemistry in Fock space. I. The universal wave and energy operators, J. Chem. Phys. 77, 3081 (1982); W. Kutzelnigg and S. Koch, Quantum chemistry in Fock space. II. Effective Hamiltonians in Fock space, ibid. 79, 4315 (1983); W. Kutzelnigg, [48] W. Kutzelnigg, Quantum Chemistry in Fock space.I. The universal wave and energy operators, J. Chem.Phys. 77, 3081 (1982); W. Kutzelnigg and S. Koch, Quantum chemistry in Fock space.II.Effective Hamiltonians in Fock space, ibid. 79, 4315 (1983); W. Kutzelnigg、
Quantum chemistry in Fock space. IV. The treatment of permutational symmetry. Spin-free diagrams with symmetrized vertices, ibid. 82, 4166 (1985); R. J. Bartlett, S. A. Kucharski, and J. Noga, Alternative coupled-cluster ansätze II. The unitary coupled-cluster method, Chem. Phys. Lett. 155, 133 (1989); W. Kutzelnigg, Error analysis and improvements of coupled-cluster theory, Theoret. Chim. Acta 80, 349 (1991); A. G. Taube and R. J. Bartlett, New perspectives on unitary coupled-cluster theory, Int. J. Quantum Chem. 106, 3393 (2006); F. A. Evangelista, G. K.-L. Chan, and G. E. Scuseria, Exact parameterization of fermionic wave functions via unitary coupled cluster theory, J. Chem. Phys. 151, 244112 (2019). フォック空間における量子化学。IV.順列対称性の扱い。(1985)、R. J. Bartlett, S. A. Kucharski, and J. Noga, Alternative coupled-cluster ansätze II.The unitary coupled-cluster method, Chem.Phys. Lett.155, 133 (1989); W. Kutzelnigg, Error analysis and improvements of coupled-cluster theory, Theoret.Chim.A. G. Taube and R. J. Bartlett, New Perspectives on unitary coupled-cluster theory, Int.J. Quantum Chem.106, 3393 (2006); F. A. Evangelista, G. K.-L. Chan, and G. E. Scuseria, Exact parameterization of fermionic wave functions via unitary coupled cluster theory, J. Chem.Phys. 151, 244112 (2019).
[49] G.-L. R. Anselmetti, D. Wierichs, C. Gogolin, and R. M. Parrish, Local, expressive, quantum-number-preserving vqe ansätze for fermionic systems, New J. Phys. 23, 113010 (2021).
[50] J. Lee, H. Q. Pham, and D. R. Reichman, Data repository for “Twenty Years of Auxiliary-Field Quantum Monte Carlo in Quantum Chemistry: An Overview and Assessment on Main Group Chemistry and Bond-Breaking” (2022), https://doi.org/10.5281/zenodo. 6816236. [50] J. Lee, H. Q. Pham, and D. R. Reichman, Data repository for "Twenty Years of Auxiliary-Field Quantum Monte Carlo in Quantum Chemistry:An Overview and Assessment on Main Group Chemistry and Bond-Breaking" (2022), https://doi.org/10.5281/zenodo.6816236.
[51] C. Lenihan, O. J. Backhouse, T. W. A. Montgomery, P. Lolur, M. J. Bhaseen, and G. H. Booth, Excita- [51] C. Lenihan, O. J. Backhouse, T. W. A. Montgomery, P. Lolur, M. J. Bhaseen, and G. H. Booth, Excita-.
tion Amplitude Sampling for Low Variance Electronic Structure on Quantum Computers, arXiv e-prints , arXiv:2506.15438 (2025), arXiv:2506.15438 [quant-ph].
[52] K. Sugisaki, S. Kanno, T. Itoko, R. Sakuma, and N. Yamamoto, Hamiltonian simulation-based quantumselected configuration interaction for large-scale electronic structure calculations with a quantum computer, arXiv e-prints , arXiv:2412.07218 (2024), arXiv:2412.07218 [quant-ph]. [52]杉崎健一郎、菅野慎一郎、糸子忠昭、佐久間亮、山本直樹、量子コンピュータを用いた大規模電子構造計算のためのハミルトニアンシミュレーションに基づく量子選択配置相互作用、arXiv e-prints , arXiv:2412.07218 (2024), arXiv:2412.07218 [quant-ph].
[53] M. Mikkelsen and Y. O. Nakagawa, Quantum-selected configuration interaction with time-evolved state, arXiv e-prints , arXiv:2412.13839 (2024), arXiv:2412.13839 [quant-ph]. このような場合、[53] M. Mikkelsen and Y. O. Nakagawa, Quantum-selected configuration interaction with time-evolved state, arXiv e-prints , arXiv:2412.13839 (2024), arXiv:2412.13839 [quant-ph].
[54] J. Yu, J. Robledo Moreno, J. T. Iosue, L. Bertels, D. Claudino, B. Fuller, P. Groszkowski, T. S. Humble, P. Jurcevic, W. Kirby, T. A. Maier, M. Motta, B. Pokharel, A. Seif, A. Shehata, K. J. Sung, M. C. Tran, V. Tripathi, A. Mezzacapo, and K. Sharma, QuantumCentric Algorithm for Sample-Based Krylov Diagonalization, arXiv e-prints , arXiv:2501.09702 (2025), arXiv:2501.09702 [quant-ph]. [54] J. Yu, J. Robledo Moreno, J. T. Iosue, L. Bertels, D. Claudino, B. Fuller, P. Groszkowski, T. S. Humble, P. Jurcevic, W. Kirby, T. A. Maier, M. Motta, B. Pokharel, A. Seif, A. Shehata, K. J. Sung, M. C. Tran, A. Mezzacapo, K. Sharma, QuantumCentric Algorithm for Sample-Based Krylov Diagonalal, V. Tripathi, A. Mezzacapo, K. Sharma.このような場合、"QuantumCentric Algorithm for Sample-Based Krylov Diagonalization", arXiv e-prints , arXiv:2501.09702 (2025), arXiv:2501.09702 [quant-ph].
[55] P. Reinholdt, K. M. Ziems, E. Rosendahl Kjellgren, S. Coriani, S. P. A. Sauer, and J. Kongsted, Critical Limitations in Quantum-Selected Configuration Interaction Methods, arXiv e-prints , arXiv:2501.07231 (2025), arXiv:2501.07231 [physics.chem-ph].