arXiv量子コンピュータニュース（2026年10月15日）

キーワード: 明示的相関とフィードバック不要の変分法 / 任意ユニタリのTゲート数指数8/7 / フォック状態による位相非依存の機械センシング / 多モード損失・フェージング通信路の有限エネルギー容量 / EFIペアと単一コピー擬似ランダム状態の同値性

オープニング：2026年10月15日 — arxiv量子コンピュータ論文解説

2026年10月14日にarXivへ掲載された量子物理の新着から、これまでの回で扱っていない5本を選んだ。ニューメキシコ大学によるクオンティニュアムH2を使った2次元ハバード模型の明示的相関計算、ローレンス・バークレー国立研究所・マイクロソフト・デューク大学による任意ユニタリ合成のTゲート数の改善、チューリッヒ工科大学と北京大学による微小グラム級音響モードでの量子増強センシング、ヨーク大学による多モード損失・フェージング通信路の有限エネルギー容量、スイス連邦工科大学ローザンヌ校による量子暗号の最小仮定の同値性である。前半の2本は、量子計算の資源をどこで節約するかを、測定後の古典処理への移し替えと、データ読み込みの分割という別々の手段で扱う。後半の3本は、実機の非古典状態、通信路の情報理論的上限、計算量的な暗号の基礎という順に、量子情報の「上限」と「仮定」を問う。

論文1: 運動量空間の行列積状態と実空間の相関因子を量子プロセッサでつなぐ

出典: Explicitly correlated ground states with a feedback-free variational hybrid algorithm. arXiv:2610.18681（2026年10月14日掲載）。著者はニューメキシコ大学クオンタム・ニューメキシコ研究所のBenjamin Corbett氏とAkimasa Miyake氏。

古典計算ではどちらかを制限しないと両立しなかった二つの表現、すなわち運動量空間で圧縮される参照状態と実空間で働く対角な相関因子を、量子プロセッサの基底回転でつなぎ、その組合せを32量子ビットの実機で示した研究なので選んだ。量子と古典の最適化ループを持たない点も、変分量子アルゴリズムの運用上の弱点への直接の回答になっている。

周期境界の2次元フェルミ・ハバード模型では、同じボンド次元なら運動量空間の行列積状態（MPS）が実空間MPSよりはるかに精度が高い。一方、グッツヴィラー因子のような物理的に動機づけられた対角相関因子は実空間で自然に働く。変分モンテカルロで両者を組み合わせるには、運動量空間MPSとスレーター行列式の重なりとして実空間振幅を計算する必要があり、これは一般に#P困難である。トランスコリレーテッド法は相関因子を相似変換ハミルトニアンに取り込んでこの衝突を避けるが、得られるエネルギーは変分的でなく、物理量には非エルミート演算子の左右固有ベクトルが要る。著者らの手順は四段階である。まず古典計算のトランスコリレーテッドDMRGで小さいボンド次元の運動量空間MPSを作る。次に量子プロセッサでそのMPSを準備し、フェルミオンのフーリエ変換で実空間へ回転して、占有数基底と各ホッピング対を対角化する基底で測定する。測定基底はハミルトニアンだけで決まり相関因子に依存しないので、相関因子 exp(Σθ_r D_r) の形とパラメータは測定後に古典処理だけで選び直し、最適化できる。最後に、MPSから古典的に分かる低次モーメント（二重占有数の期待値など）を再現するよう、サンプルの重みを最小ノルムの乗法的変化で補正する。

$$E_{\mathrm{ec}}(\boldsymbol\theta)=\frac{\langle\phi|C_{\boldsymbol\theta}HC_{\boldsymbol\theta}|\phi\rangle}{\langle\phi|C_{\boldsymbol\theta}^{2}|\phi\rangle},\qquad C_{\boldsymbol\theta}=\exp\Big(\sum_r\theta_rD_r\Big)$$

半充填・U/t=4の4×4格子の古典シミュレーションでは、ボンド次元16のMPSに103パラメータの相関因子を組み合わせると−13.609±0.005（単位t）に達した。同じ−13.6に届くには、トランスコリレーションなしの運動量空間DMRGでボンド次元1024、実空間なら4096が要り、パラメータ数は141分の1になる。実機実験はクオンティニュアムの56量子ビット機H2-2のうち32量子ビットを使い、ボンド次元4のMPSを2層の階段型回路で準備した。コンパイル後の対角測定回路は2量子ビットゲート555個、1量子ビットゲート1,095個で、回路は一次元の最近接結合だけで済む。対角12,250、ホッピング23,616のサンプルのうち、上下スピンが各8個の対称性セクターに残ったのは7,526と14,509だった。グッツヴィラー因子で最適化したエネルギーは、生データで−9.21±0.49、事後選択で−11.49±0.49、較正後に−13.14±0.07となり、雑音なしの基準値−13.2914±0.0028に近づいた。基底エネルギーに対する相対誤差は16±4%から3.5±0.5%へ下がり、較正は標準誤差を6.6分の1、基準値からのずれを12.2分の1にした。比較では、量子選択配置間相互作用法（QSCI）が厳密な基底状態から雑音なしにサンプルできても同じ精度に約2×10⁸サンプルを要し、本実験の3.6×10⁴より4桁近く多い。量子クリロフ法は、ボンド次元64の実空間MPSから始めて雑音を無視しても、2量子ビットゲート544個（本手法の回転回路は464個）を要する。

限界も具体的である。比の推定量には有限サンプルでの厳密な変分上界がなく、状態準備に誤りがあると較正によって上界の保証も失われる。実機の推定値は、同じサンプル数の雑音なし試行1,000回のうち95.2%より高く、事後選択と較正の後にも上方へのずれが残った。半充填では必要サンプル数が格子サイト数Nに対して exp(cN)、c<0.1で増えると著者らは仮説を立てており、6×6格子では4×4に比べて生データで約3.4倍、較正後では約30倍のサンプルが要った。著者ら自身も原理実証と位置づけ、量子優位には補助場量子モンテカルロが符号問題を抱える半充填から外れた大きな格子が必要だとしている。スピン相関関数では、実機の推定値がトランスコリレーテッド法の双直交推定より厳密値に近く、非変分的なエネルギーの良さが物理量の正確さを意味しないことも示された。

論文2: 任意ユニタリ合成のTゲート数を指数8/7まで下げる

出典: T-efficient unitary synthesis via block-sparse Hamiltonian simulation. arXiv:2610.18589（2026年10月14日掲載）。著者はローレンス・バークレー国立研究所のDiyi Liu氏、Wibe Albert de Jong氏、Chao Yang氏、マイクロソフト・クオンタムのShuchen Zhu氏、David Williams-Young氏、デューク大学のDi Fang氏。

任意のn量子ビットユニタリをクリフォード+T回路で作るときの最悪Tゲート数について、上界の指数を5/4から8/7へ下げ、その指数が特定の計算モデル内では最適であることも同時に示した理論結果なので選んだ。誤り耐性計算ではTゲートがマジック状態蒸留で供給されるため、データ読み込みを含めた回路全体のT数が直接のコスト指標になる。

N=2ⁿとすると、一般のユニタリは N² 個の実パラメータを持ち、数え上げによる下界は Ω(N) である。上界はLow氏らの N^{3/2}、タン氏の再帰的コサイン・サイン分解による4/3、ユアン氏らのランダム・ブロック平坦化と特異値変換による5/4と更新されてきた。ユアン氏らの構成では、各クエリが N² 個すべての行列要素を含むルックアップ表を参照し、特異値増幅のたびにその読み込みコストを払う。本論文は、目標ユニタリのエルミート対数Hをとり、トレースを除いた部分にアダマール変換とランダム符号行列を2層重ねた相似変換 P=F S_τ F S_σ F を施して、b×bブロックすべてのノルムを ‖H‖/D（D=N/b）程度に平坦化する。符号行列1層では、Hをアダマール基底で対角化したものと可換になって効かない入力があるため2層が要る。次に、エルミート性を保つよう非対角ブロックの対を一緒にして、s個の部分にランダムに振り分ける。各部分のルックアップ表は D²/s 程度に縮み、1クエリのTコストは Õ(N/s+b²) になる。これらの部分の時間発展をqDRIFTで組み合わせ、各区間をQSVTで実装すると、総クエリ数は Õ(D+s) となる。s を D 程度に取って b=Θ(N^{3/7}) で二つの項を釣り合わせると、Tゲート数が Õ(N^{8/7})、クリーンな補助量子ビットが Õ(N^{6/7}) になる。部分を増やすと表は小さくなるが必要な区間数が増える、という交換条件がこの指数を決めている。

精度への依存は、残差を逐次補正する「残差精緻化」で改善した。誤差 ε に対して L=n+log₂(1/ε) とおくと、L≤N の範囲で Tゲート数 Õ(N^{8/7}L^{3/7}+NL)、補助量子ビット Õ(N) となり、単一のqDRIFT近似で生じる ε^{−6/7} の依存を避ける。N^{1/4}≤L≤N では NL の項が支配的になる。下界は「逐次ブロックSELECT模型」の中で示された。この模型は目標をハミルトニアン発展の列として表し、ルックアップ、ブロック内処理、シミュレーションに費用を課す一方、発展の間に挟むユニタリが運べる情報量を N²log(1/ε)/4 ビット以下に制限する。その中では、一定精度で Ω̃(N^{8/7}) が最悪ケースで必要になり、構成と一致する。限界として、Õ と Ω̃ は n の多項式因子を隠しており、定理4の明示形では n³ がかかる。結果は十分大きな n での漸近評価で、具体的な小規模ユニタリでの定数因子は示されていない。最適性は模型内の主張であり、模型の外では数え上げの Ω(N) との間に N^{1/7} の隙間が残る。サンプリングはすべて古典の前処理で行うため、回路自体は決定的に固定される。

論文3: 7フォノンまでのフォック状態で微小グラム級の音響モードを位相非依存のセンサーにする

出典: Quantum-enhanced mechanical sensing with machine-learning estimation. arXiv:2610.18776（2026年10月14日掲載）。著者はチューリッヒ工科大学のFinn Müßigbrodt氏、Yiwen Chu氏、北京大学とチューリッヒ工科大学のMatteo Fadel氏ら6名。

スクイーズド状態による機械系の感度向上は既知の直交位相に限られ、信号の位相が分からない場合に効く非ガウス状態の計測上の利点は、これまでイオンのような軽い系でしか示されていなかった。それを有効質量約16マイクログラムの音響モードで実証し、古典的な探針の上限を明確に超えた実験なので選んだ。10月12日回のスタンフォード大学による量子音響の量子ジャンプ実験は測定の量子非破壊性を扱ったが、本論文はフォノン数分布そのものを計測資源として使う。

素子はトランズモン量子ビットを高次倍音バルク音響波共振器（HBAR）にフリップチップ接合した回路量子音響力学系である。音響モードは5.088ギガヘルツ、エネルギー緩和時間196マイクロ秒、ラムゼー位相緩和時間401マイクロ秒で、量子ビットとの結合は2π×269キロヘルツ、量子ビットのT1は27マイクロ秒である。最適制御パルスで真空からフォック状態 |k⟩ を準備し、忠実度は k=1 で93.6%、k=4 で92.0%、k=7 で67.4%だった。検出対象の信号は音響モードに共鳴するマイクロ波駆動で模擬し、位相空間での変位の振幅 d を推定する。真空を探針にして変位を較正し、0から0.832までの31点の格子を作った。どの方向の変位にも同じ応答を示すフォック状態は、平均フォノン数を固定したとき位相非依存の変位推定で最適であり、理想的には量子フィッシャー情報が 4(2k+1) となる。一方、コヒーレント状態の混合で書ける古典的な探針は、エネルギーによらずフィッシャー情報が4を超えない。

フィッシャー情報は二つの方法で独立に求めた。近接した変位でのフォノン数分布のヘリンジャー距離の曲率からは、|4⟩ で d≃0.3 において F=21.2(1.3) を得た。滑らかな応答モデルを当てはめたベイズ推定からは、3≤k≤7 の探針で F≃20 が得られ、両者は一致した。古典的上限4の5倍であり、7.0デシベルの感度向上、同じ試行回数での変位不確かさ約2.2分の1に当たる。力の感度に換算すると、5ギガヘルツの力に対して0.77ピコニュートン毎ルートヘルツになる。ただし F は |4⟩ で頭打ちになり、それ以上のフォノン数では約20にとどまった。損失だけのモデルでは透過率 η≃0.97 から飽和値は約67で、それに達するのは k≫50 なので、観測された飽和は説明できない。準備誤りの分を |k−1⟩ に割り当てるモデルは平均フォノン数4以上で約25に飽和し、測定とよく合う。つまり現状の上限は寿命ではなく状態準備の忠実度で決まっている。

後半では、解析モデルを使わない推定器として、測定したフォノン数1個を入力し31個の較正変位への事後確率を出力するニューラルネットワーク（20ニューロンの隠れ層2層）を較正データで学習した。100ショットの記録から得た事後分布は、モデルに基づく推定と同じ位置・同じ幅を示し、同じ感度向上を再現した。較正範囲全体で平均した推定バイアスの絶対値は、モデル推定の6.1×10⁻³から3.2×10⁻³へ約半分になった。著者らは、入力がフォノン数1個のスカラーなのでネットワークは条件付き確率の表以上の表現力を持たず、その役割は希な出力での標本ゼロや大きな相対揺らぎを、関数形を仮定せずに平滑化することだと明記している。応用として高周波重力波、暗黒物質、波動関数の収縮模型の検証を挙げるが、本論文で示されたのは共鳴駆動で模擬した信号の推定までである。

論文4: 多モード損失とフェージング通信路の有限エネルギー双方向容量を厳密に求める

出典: Finite-energy two-way assisted capacities: Exact formulas for multimode loss and fading. arXiv:2610.18731（2026年10月14日掲載）。著者はヨーク大学のStefano Pirandola氏と、ロンドンのnodeQのGaetana Spedalieri氏。

純損失通信路の双方向容量の上限として広く使われるPLOB限界は、入力エネルギーに制約がない場合の値である。実際の光通信は有限の光子数予算を多数のモードへ配分し、透過率も時間的に揺らぐ。その状況での量子通信・エンタングルメント配送・秘密鍵の双方向容量を、任意の適応的プロトコルを許したうえで厳密な式として与えた結果なので選んだ。

透過率 η の純損失通信路について、無制約の双方向容量は −log₂(1−η) であり、平均光子数 N のもとでは逆コヒーレント情報 g(N)−g((1−η)N) が達成可能な速度として知られていた（g はボソンのエントロピー関数）。先行研究（論文中の参考文献8）は、セクター・テレポーテーションによる模擬で、この速度が平均光子数制約下の任意の適応的プロトコルに対して最適であることを示していた。本論文はこれを m モードに拡張し、容量が各モードへの光子配分 N_j の最適化 max Σ c_{η_j}(N_j)、ΣN_j≤N で与えられ、量子通信、エンタングルメント配送、秘密通信、秘密鍵の4つの容量がすべて一致することを示した。制約は全分岐・全試行を含む無条件の平均であり、各モードに N 個ずつの予算ではない。最適配分は一意で、透過するモードにはすべて正の光子数が割り当てられ、透過率の高いモードほど多くの光子を受け取る。条件は各モードの限界的な容量増分が等しいことである。例えば η=(0.8, 0.4, 0.1)、N=3 では、最適配分は約 (2.13, 0.76, 0.11)、容量は3モード1回あたり約2.04である。受動的なモード混合を含む通信路では、容量は透過行列の特異値の2乗だけで決まる。

周波数に依存する損失と入力パワー制約へ拡張すると、最適な占有数は c′(ν)=λħω という1個のラグランジュ乗数の条件で決まる。損失が一定で帯域が無制限なら容量は閉じた式で書け、入力パワーの平方根に比例する。フェージング通信路では、透過率の実現値が送信前に公開される場合と送信後に公開される場合を区別して厳密な容量を求め、観測されない場合には上下界を与えた。確率 p で透過率 η₀、それ以外は完全遮断の2値フェージングでは、送信前の情報があれば容量は p·c(N/p)、送信後なら p·c(N) となる。前者は「通る」と分かっている時刻にエネルギーを集中できる分だけ大きい。p=1/2、N=0.1 のオン・オフ通信路では、送信前の情報で約0.390、送信後で約0.242、観測されない場合は約0.072から約0.138の間に挟まれる。

強い逆定理のしきい値、すなわちそれを超えると誤りが確実に1へ近づく速度も求めた。平均エネルギー制約のもとでは、このしきい値は無制約のPLOB値の和に等しく、有限エネルギー容量より一般に大きい。平均だけを縛る制約では、公開の乱数で低確率に高エネルギーのプロトコルを走らせる戦略が許されるためである。帯域無制限・損失一定の場合、容量は有限なのにしきい値は無限大になる。著者らは、観測されないフェージングでエンタングルメント配送と秘密鍵の容量が一致するかは未解決とし、熱雑音やランダム位相、連続するフェージングの相関（ブロック・フェージング）はモデルに含まれないと明記している。連続周波数の定理は独立に扱えるスペクトルモードを仮定しており、有限時間の光学散乱問題への適用には時間境界での結合の制御が別に要る。なお論文末尾には、候補となる数学的導出と証明の生成に生成AIを使い、すべて著者が検証・修正したとの記載がある。

論文5: EFIペアと伸長型の単一コピー擬似ランダム状態生成器の同値性

出典: EFI⇔1PRS. arXiv:2610.18799（2026年10月14日掲載）。著者はスイス連邦工科大学ローザンヌ校のRishabh Batra氏の単著。

量子暗号の「最小の仮定」とされるEFIペアについて、2026年のCavalar氏らの同値性証明が残していた非一様性、すなわち対数長の古典アドバイスを取り除き、量子版の「一方向関数・コミットメント・擬似乱数生成器の同値性」の図式を完成させた結果なので選んだ。

古典暗号では、一方向関数、擬似乱数生成器、計算量的に安全なコミットメント、効率的に標本できて統計的には遠いが計算量的に区別できない分布の対（EFID）が、存在の仮定として同値である。量子では、擬似ランダム状態はBQP=QMAとなる量子オラクルのもとでも存在しうること、単一コピー版はP=NPとなる古典オラクルのもとでも存在しうることが知られ、一方向関数より弱い仮定で暗号が組める可能性が議論されてきた。EFIペアは、効率的に生成でき、統計的には遠いが計算量的には区別できない量子状態の対で、量子コミットメントと同値である。伸長型の単一コピー擬似ランダム状態生成器（1PRS）は、短い古典鍵から純粋状態を作り、鍵について平均した状態が1コピーでは最大混合状態と区別できず、しかも出力の量子ビット数が鍵のビット数より多いものを指す。伸長を課さないと、鍵kを基底状態|k⟩に写すだけで自明に条件を満たしてしまう。Cavalar氏らは、抽出長を指定する O(log λ) ビットのアドバイスつきでEFIと1PRSの同値性を示していた。アドバイスが要ったのは、EFI状態のエントロピーの基準値が未知で、抽出する長さを決められなかったためである。

本論文の工夫は、EFI生成器を純粋化し、その純粋化レジスタを手元に残すことにある。こうすると、1ビットのフラグと純粋化した状態からなる源のエントロピーはちょうど1ビットになり、個々のEFI状態のエントロピーを推定せずに済む。そのうえで、各接頭部分について、実際の周辺状態と計算量的に区別できず、直前までの周辺状態は厳密に一致する局所的な「証人」状態を用意し、条件付きエントロピーの総和が 1+χ≥2−negl になることを示した（χは2状態の統計的分離から来る量）。個々の項は未知で負にもなりうるので、Haitner氏らのランダムスタートの手法を量子側情報つきに拡張し、独立な鍵の源を L=4M 個連結して開始位置を一様に選ぶことで、各座標の条件付きエントロピーを平均値で下から抑えた。漸近等分配の後、すべての列で同じ抽出長を使って抽出し、残りの量子ビットを秘密のパウリで覆う。出力長と鍵長の差は (M−3/2)t>0 となり、伸長が得られる。帰結として、EFI、統計的に検証可能な一方向状態生成器、量子コミットメント、1PRSが存在の仮定として同値になり、Goldin氏らが「ナノクリプト」に分類した原始要素もすべて同値になる。さらに、ガーチ氏の量子アルゴリズム的エントロピーに基づくギャップ問題の平均困難性によって、EFIの一様なメタ計算量的特徴づけが初めて得られた。結果は存在の同値性であり、具体的な候補構成の安全性や、構成の効率に関する数値は与えていない。

まとめ

前半の2本は、量子計算のコストをどこへ移すかの例だった。Corbett氏とMiyake氏の手法は、相関因子の最適化を測定後の古典処理へ移し、32量子ビットの実機で相対誤差3.5%に達したが、サンプル数はサイト数に対して指数的に増える。Liu氏らは、データ読み込みをブロック疎な部分へ分けてTゲート数の指数を8/7に下げ、特定の模型内ではそれが最適であることを示した。後半の3本は、上限と仮定を扱った。チューリッヒ工科大学の実験は、古典的探針の上限の5倍のフィッシャー情報を微小グラム級の音響モードで示し、頭打ちの原因が状態準備にあることを定量化した。Pirandola氏らは、有限の光子数予算での双方向容量を多モードとフェージングへ広げ、送信前の情報が持つ価値を数式で示した。Batra氏は、量子暗号の最小仮定をめぐる同値性からアドバイスを取り除いた。

参考ソース

• 論文1: Explicitly correlated ground states with a feedback-free variational hybrid algorithm. arXiv:2610.18681
• 論文2: T-efficient unitary synthesis via block-sparse Hamiltonian simulation. arXiv:2610.18589
• 論文3: Quantum-enhanced mechanical sensing with machine-learning estimation. arXiv:2610.18776
• 論文4: Finite-energy two-way assisted capacities: Exact formulas for multimode loss and fading. arXiv:2610.18731
• 論文5: EFI⇔1PRS. arXiv:2610.18799