{ "cells": [ { "cell_type": "markdown", "id": "b1267e87-e0b3-4934-b4b8-0265fe84adfd", "metadata": {}, "source": [ "---\n", "title: \"エラー軽減および抑制技術\"\n", "description: \"IBM Quantum Compute Service で利用可能なエラー軽減および抑制手法の説明\"\n", "---\n", "\n", "{/* cspell:ignore extrapolators */}\n", "\n", "\n", "\n", "# エラー軽減および抑制技術\n", "\n" ] }, { "cell_type": "markdown", "id": "349a0e50", "metadata": {}, "source": [ "\n", " 新しい実行モデルのベータ版がリリースされた。 有向実行モデルは、エラー軽減のワークフローをカスタマイズする際に、より柔軟性を提供する。 詳しくは [Directed実行モデル](/docs/guides/directed-execution-model)ガイドを参照。\n", "\n", "\n" ] }, { "cell_type": "markdown", "id": "7fb66fc9-c519-4904-a6a1-33adc3f8df13", "metadata": { "tags": [ "version-info" ] }, "source": [ "{/*\n", " DO NOT EDIT THIS CELL!!!\n", " This cell's content is generated automatically by a script. Anything you add\n", " here will be removed next time the notebook is run. To add new content, create\n", " a new cell before or after this one.\n", " */}\n", "\n", "\n", " \n", " このページのコードは、以下の要件に基づいて開発された。\n", " これらのバージョンまたは新しいバージョンの使用をお勧めします。\n", "\n", " ```\n", " qiskit-ibm-runtime~=0.47.0\n", " ```\n", " \n", "\n", "\n" ] }, { "cell_type": "markdown", "id": "076872c0-f8de-4aef-ae13-29078ee8564e", "metadata": {}, "source": [ "エラーの軽減および抑制技術は、より大規模なワークロードへスケールアップする際に、結果の品質を向上させるために用いられます。 このページでは、 IBM Quantum Compute で利用可能なエラー抑制およびエラー軽減の手法について、概要を説明します。\n", "\n", "次のセルでは、Estimatorプリミティブをインポートし、後のコードセルでEstimatorの初期化に使用するバックエンドを作成しています。\n", "\n" ] }, { "cell_type": "code", "execution_count": 1, "id": "39ef65e5", "metadata": {}, "outputs": [], "source": [ "from qiskit_ibm_runtime import EstimatorV2 as Estimator\n", "from qiskit_ibm_runtime import QiskitRuntimeService\n", "\n", "service = QiskitRuntimeService()\n", "backend = service.least_busy()" ] }, { "cell_type": "markdown", "id": "ad548ef5-f1c1-4def-b551-b3ece2032065", "metadata": {}, "source": [ "\n", "\n", "## 動的分離\n", "\n", "量子回路は、 IBM®、マイクロ波パルスのシーケンスとしてハードウェア上で実行され、正確な時間間隔でスケジュールして実行する必要がある。\n", "残念ながら、量子ビット間の不要な相互作用は、アイドリング中の量子ビットのコヒーレントエラーにつながる可能性がある。 動的デカップリングは、アイドリング中の量子ビットにパルス列を挿入して、これらの誤差の影響をほぼ打ち消すことによって機能する。 挿入された各パルス・シーケンスはID演算に相当するが、パルスの物理的な存在はエラーを抑制する効果がある。\n", "パルスシーケンスには多くの選択肢があり、特定の症例ごとにどのシーケンスが良いかは、依然として[活発な研究分野](https://journals.aps.org/prapplied/abstract/10.1103/PhysRevApplied.20.064027)である。\n", "\n", "ダイナミック・デカップリングは、主に、いくつかの量子ビットが何も操作されずにアイドル状態になっているギャップを含む回路に有効であることに注意されたい。 回路内の演算が非常に高密度に詰め込まれており、すべての量子ビットがほとんどの時間ビジー状態である場合、動的デカップリング・パルスを追加しても性能が向上しない可能性がある。 それどころか、パルス自体の不完全性によって性能が悪化する可能性さえある。\n", "\n", "下図はXXパルスシーケンスによるダイナミックデカップリングを示している。 左の抽象的な回路は、右上のマイクロ波パルス・スケジュールにマッピングされている。 右下は同じスケジュールを示しているが、最初の量子ビットのアイドル期間に2つのXパルスが挿入されている。\n", "\n", "![ダイナミック・デカップリングの描写](https://quantum.cloud.ibm.com/docs/images/guides/error-mitigation-explanation/dd.avif)\n", "\n", "動的デカップリングは、 [動的デカップリング・オプションで](/docs/api/qiskit-ibm-runtime/options-dynamical-decoupling-options) `enable` を `True` に設定することで有効にできる。 `sequence_type` オプションを使用すると、複数の異なるパルスシーケンスから選択できる。 デフォルトのシーケンスタイプは `\"XX\"`。\n", "\n", "以下のコードセルでは、Estimator の動的デカップリングを有効にし、動的デカップリングシーケンスを選択する方法を示しています。\n", "\n" ] }, { "cell_type": "code", "execution_count": 2, "id": "cc057e3f-f592-434e-8bb2-96f901564318", "metadata": { "jp-MarkdownHeadingCollapsed": true, "vscode": { "languageId": "plaintext" } }, "outputs": [], "source": [ "estimator = Estimator(mode=backend)\n", "estimator.options.dynamical_decoupling.enable = True\n", "estimator.options.dynamical_decoupling.sequence_type = \"XpXm\"" ] }, { "attachments": {}, "cell_type": "markdown", "id": "56ed5785-e3a8-4dc5-8b55-ef876fccd592", "metadata": { "editable": true, "slideshow": { "slide_type": "" }, "tags": [] }, "source": [ "\n", "\n", "## パウリの回転\n", "\n", "[ランダム化コンパイルとしても](https://journals.aps.org/pra/abstract/10.1103/PhysRevA.94.052325)知られるトワーリングは、任意のノイズ・チャンネルをより特殊な構造を持つノイズ・チャンネルに変換するために広く使われている技術である。\n", "\n", "パウリ・ツワーリングは、パウリ演算を用いた特殊なツワーリングである。 量子チャンネルをパウリチャンネルに変換する効果がある。 パウリノイズは線形に蓄積されるのに対し、コヒーレントノイズは演算回数に応じて二次関数的に蓄積される傾向があるためだ。 パウリ・ツワーリングはしばしば、任意のノイズよりもパウリ・ノイズの方が効果的な他のエラー緩和技術と組み合わされる。\n", "\n", "パウリ・ツワーリングは、ゲートの理想的な効果が変わらないように、ランダムに選ばれた1量子ビットのパウリ・ゲートで、選ばれたゲートのセットを挟むことで実装される。 その結果、単一の回路が、すべて同じ理想的な効果を持つランダムな回路のアンサンブルに置き換えられる。 回路をサンプリングする場合、サンプルは1つだけではなく、複数のランダムなインスタンスから抽出される。\n", "\n", "![回転するパウリの描写](https://quantum.cloud.ibm.com/docs/images/guides/error-mitigation-explanation/pauli_twirling.svg)\n", "\n", "現在の量子ハードウェアのエラーのほとんどは2量子ビットゲートに起因するため、この技術は（ネイティブの）2量子ビットゲートにのみ適用されることが多い。 次の図は、CNOTゲートとECRゲートのパウリ・ツワールを示している。 列内のすべての回路が同じ理想的な効果を持つ。\n", "\n", "![ゲート回転の描写](https://quantum.cloud.ibm.com/docs/images/guides/error-mitigation-explanation/gate_twirls.avif)\n", "\n", "パウリ・ツワーリングは、 [ツワーリング・オプションで](/docs/api/qiskit-ibm-runtime/options-twirling-options) `enable_gates` を `True` に設定することで有効になる。 その他の注目すべきオプションは以下の通り：\n", "\n", "* `num_randomizations`:撚り合わされた回路のアンサンブルから引き出される回路インスタンスの数。\n", "* `shots_per_randomization`:各回路インスタンスからサンプリングするショット数。\n", "\n", "以下のコードセルでは、パウリ・トゥワーリングを有効にし、Estimatorに対してこれらのオプションを設定する方法を示しています。 これらのオプションはいずれも、明示的に設定する必要はありません。\n", "\n" ] }, { "cell_type": "code", "execution_count": 3, "id": "c3e12d44", "metadata": {}, "outputs": [], "source": [ "estimator = Estimator(mode=backend)\n", "estimator.options.twirling.enable_gates = True\n", "estimator.options.twirling.num_randomizations = 32\n", "estimator.options.twirling.shots_per_randomization = 100" ] }, { "attachments": {}, "cell_type": "markdown", "id": "ad9d9159-4277-4aec-bf26-da748e240f38", "metadata": { "editable": true, "slideshow": { "slide_type": "" }, "tags": [] }, "source": [ "\n", "\n", "\n", "\n", "## TREX（ eXtinction ）\n", "\n", "[Twirled Readout Error（TREX： eXtinction ）](https://journals.aps.org/pra/abstract/10.1103/PhysRevA.105.032620) は、パウリ観測量の期待値を推定する際の測定誤差の影響を軽減する。\n", "これは「ツイールド測定」という概念に基づいており、測定ゲートを、(1) パウリXゲート、(2) 測定、(3) 古典的なビット反転、という一連の操作でランダムに置き換えることによって実現される。 標準的なゲートツイリングと同様に、このシーケンスは、次の図に示すように、ノイズがない場合の単純な測定に相当します\n", "\n", "![測定時の回転の様子](https://quantum.cloud.ibm.com/docs/images/guides/error-mitigation-explanation/measurement_twirling.svg)\n", "\n", "読み出し誤差が存在する場合、測定の回転は読み出し誤差伝達行列を対角化するという効果があり、これにより行列の逆行列を求めることが容易になる。 ゼロ状態で初期化されたランダム回路のベンチマークを行うことにより、\n", "対角ノイズチャネルからのスケーリング項が学習される。 これにより、\n", "本サービスは、読み出しノイズに起因する\n", "期待値からバイアスを除去することができる。 読み出し誤差伝達行列を推定するには、追加の校正回路を実行する必要があり、それによってわずかなオーバーヘッドが生じる。\n", "\n", "\n", "\n", "### オプションを使用してTREXを `measure_mitigation` 有効にします\n", "\n", "TREXを有効にするには、 [IBM Quantum](/docs/api/qiskit-ibm-runtime/options-resilience-options-v2) の「Estimatorの計算耐障害性オプション」で `True` を に `measure_mitigation` 設定します。 [ここでは、](/docs/api/qiskit-ibm-runtime/options-measure-noise-learning-options) 測定ノイズの学習に関するオプションについて説明します。 ゲートの回転と同様に、回路のランダム化回数や、1回のランダム化あたりのショット数を設定できます。\n", "\n", "以下のコードセルでは、TREXを有効にし、Estimatorに対してこれらのオプションを設定する方法を示しています。 これらのオプションはいずれも、明示的に設定する必要はありません。\n", "\n" ] }, { "cell_type": "code", "execution_count": 4, "id": "80c858e3", "metadata": {}, "outputs": [], "source": [ "estimator = Estimator(mode=backend)\n", "estimator.options.resilience.measure_mitigation = True\n", "estimator.options.resilience.measure_noise_learning.num_randomizations = 32\n", "estimator.options.resilience.measure_noise_learning.shots_per_randomization = 100" ] }, { "cell_type": "markdown", "id": "5a0c3f58", "metadata": {}, "source": [ "\n", "\n", "### オプション `resilience_level` を使用してTREXを有効にします\n", "\n", "次の例に示すように、 [レジリエンスレベル](/docs/guides/estimator-noise-management#resilience)を に設定することで `1`、TREXを有効にできます。\n", "\n" ] }, { "cell_type": "code", "execution_count": 5, "id": "de92bffd", "metadata": {}, "outputs": [], "source": [ "from qiskit_ibm_runtime import QiskitRuntimeService\n", "from qiskit_ibm_runtime import EstimatorV2 as Estimator\n", "\n", "service = QiskitRuntimeService()\n", "backend = service.least_busy(operational=True, simulator=False)\n", "\n", "# Setting options during primitive initialization\n", "estimator = Estimator(backend, options={\"resilience_level\": 1})" ] }, { "attachments": {}, "cell_type": "markdown", "id": "7f345691-8330-40f5-a5c6-d6c5fb93409a", "metadata": { "editable": true, "slideshow": { "slide_type": "" }, "tags": [] }, "source": [ "\n", "\n", "\n", "\n", "## ゼロノイズ外挿法（ZNE）\n", "\n", "ゼロノイズ外挿（ZNE）は、観測値の期待値を推定する際の誤差を軽減するための手法である。 それによって結果が改善されることも多いが、偏りのない結果が得られる保証はない。\n", "\n", "ZNEは2つの段階からなる：\n", "\n", "1. *ノイズ増幅* ：元の量子回路を異なるノイズ率で複数回実行する。\n", "2. *外挿* ：理想的な結果は、ノイズの多い期待値の結果をゼロノイズ限界まで外挿することによって推定される。\n", "\n", "ノイズ増幅段階と外挿段階は、いずれもさまざまな方法で実装することができます。 IBM Quantum Computeは、「デジタルゲートフォールディング」によってノイズ増幅を実現しています。これは、2量子ビットゲートを、そのゲートとその逆ゲートからなる等価なシーケンスに置き換えることを意味します。 たとえば、単一の $U$ を $U U^\\dagger U$ に置き換えると、ノイズ増幅係数は3になります。 外挿を行う際には、線形近似や指数近似など、いくつかの関数形式の中から選択することができます。\n", "下の図は、左側がデジタルゲートの折り畳み、右側が外挿の手順を示しています。\n", "\n", "![ZNEの描写](https://quantum.cloud.ibm.com/docs/images/guides/error-mitigation-explanation/zne.avif)\n", "\n", "この手法のオーバーヘッドは、ノイズ因子の数に比例して増加する。 デフォルト設定では、3つのノイズ係数で期待値をサンプリングするため、\n", "おおよそ3xのオーバーヘッドが生じます。\n", "\n", "\n", "\n", "### オプション `zne_mitigation` を使用してZNEを有効にする\n", "\n", "ZNE を有効にするには、 [IBM Quantum](/docs/api/qiskit-ibm-runtime/options-resilience-options-v2) の「Estimator」の「Compute resilience」オプション `True` で を に `zne_mitigation` 設定します。\n", "ZNEのオプションについては、 [こちらで](/docs/api/qiskit-ibm-runtime/options-zne-options)説明しています。 特に注目すべきオプションは以下の通りです：\n", "\n", "* `noise_factors`:ノイズ増幅に使用するノイズ係数。\n", "* `extrapolator`:外挿に使用する関数形式。\n", "\n", "以下のコードセルでは、ZNEを有効にし、Estimatorのこれらのオプションを設定する方法を示しています。 これらのオプションはいずれも、明示的に設定する必要はありません。\n", "\n" ] }, { "cell_type": "code", "execution_count": 6, "id": "06de4431", "metadata": {}, "outputs": [], "source": [ "estimator = Estimator(mode=backend)\n", "estimator.options.resilience.zne_mitigation = True\n", "estimator.options.resilience.zne.noise_factors = (1, 3, 5)\n", "estimator.options.resilience.zne.extrapolator = \"exponential\"" ] }, { "cell_type": "markdown", "id": "600ee406", "metadata": {}, "source": [ "\n", "\n", "### オプション `resilience_level` を使用してZNEを有効にする\n", "\n", "次の例に示すように、 [レジリエンスレベル](/docs/guides/estimator-noise-management#resilience)を に設定することで、ZNE（およびTREX、 `2`ゲート・トゥワーリング）を有効にできます。\n", "\n" ] }, { "cell_type": "code", "execution_count": 7, "id": "7b644752", "metadata": {}, "outputs": [], "source": [ "from qiskit_ibm_runtime import QiskitRuntimeService\n", "from qiskit_ibm_runtime import EstimatorV2 as Estimator\n", "\n", "service = QiskitRuntimeService()\n", "backend = service.least_busy(operational=True, simulator=False)\n", "\n", "# Setting options during primitive initialization\n", "estimator = Estimator(backend, options={\"resilience_level\": 2})" ] }, { "cell_type": "markdown", "id": "e5dcab7b", "metadata": {}, "source": [ "\n", "\n", "\n", "\n", "## 確率的誤差増幅（PEA）\n", "\n", "ZNEにおける主な課題の1つは、ターゲット回路に影響を与えるノイズを正確に増幅することである。 ゲートフォールディングは、この増幅を簡単に行う方法を提供するが、不正確な可能性があり、誤った結果につながる可能性がある。 詳細は、論文 [\"Scalable error mitigation for noisy quantum circuits produces competitive expectation values\"](https://arxiv.org/pdf/2108.09197)、特に補足情報の4ページを参照のこと。 確率的誤差増幅は、ノイズ学習による誤差増幅に対してより正確なアプローチを提供する。\n", "\n", "PEAはより洗練された技術であり、ノイズを再構築するために予備実験を行い、その情報を使って正確な増幅を行う。 これは、回路の各レイヤーのエンタングルゲートが実行される前に、そのエンタングルゲートの渦巻きノイズモデルを学習することから始まる（関連する学習オプションを参照）。 [LayerNoiseLearningOptions](/docs/api/qiskit-ibm-runtime/options-layer-noise-learning-options) を参照)。 学習段階の後、回路は各ノイズ係数で実行され、回路の各エンタングル層は、対応する学習ノイズモデルに比例した1量子ビットノイズを確率的に注入することで増幅される。 詳しくは [「フォールト・トレランス以前の量子コンピューティングの有用性の証拠」という](https://www.nature.com/articles/s41586-023-06096-3)記事を参照されたい。\n", "\n", "PEAは3つの段階からなる：\n", "\n", "1. *学習* ：回路の各層のエンタングルゲートの渦巻きノイズモデルを学習する。\n", "2. *ノイズ増幅* ：元の量子回路を異なるノイズ係数で複数回実行する。\n", "3. *外挿* ：理想的な結果は、ノイズの多い期待値の結果をゼロノイズ限界まで外挿することによって推定される。\n", "\n", "実用規模の実験では、PEAが最良の選択となることが多い。\n", "\n", "PEAはZNEノイズ増幅技術であるため、 `resilience.zne_mitigation = True` を設定してZNEを有効にする必要もある。 その他の [`resilience.zne`](/docs/api/qiskit-ibm-runtime/options-zne-options) オプションはさらに、エクストラポレーターや増幅レベルなどを設定するために使うことができる。 PEAは、プリミティブを使用する際に自動的に生成されるノイズモデルを必要とする。\n", "\n", "次のスニペットは、Estimatorジョブの結果を軽減するためにPEAを使用した例です：\n", "\n" ] }, { "cell_type": "code", "execution_count": 8, "id": "4947f442", "metadata": {}, "outputs": [], "source": [ "estimator = Estimator(mode=backend)\n", "estimator.options.resilience.zne_mitigation = True\n", "estimator.options.resilience.zne.amplifier = \"pea\"" ] }, { "cell_type": "markdown", "id": "fcc40cb7-45b7-4248-84c4-aa7f26ba50ff", "metadata": {}, "source": [ "\n", "\n", "\n", "\n", "## 確率的誤差補償（PEC）\n", "\n", "確率的エラーキャンセル(PEC)は、観測値の期待値を推定する際のエラーを軽減する技術である。 ZNEとは異なり、期待値の不偏推定値を返す。 しかし、一般的にはオーバーヘッドが大きくなる。\n", "\n", "PECでは、理想的なターゲット回路の効果は、実際に実装可能なノイズ回路の線形結合として表現される：\n", "\n", "$$\n", "\\mathcal{O}_{\\text{ideal}} = \\sum_{i} \\eta_i \\mathcal{O}_{noisy, i}\n", "$$\n", "\n", "理想回路の出力は、線形結合で定義されたランダムアンサンブルから抽出された異なるノイズ回路インスタンスを実行することで再現できる。 係数 $\\eta_i$ が確率分布を形成していれば、アンサンブルの確率として直接使うことができる。 実際には、係数のいくつかは負なので、代わりに準確率分布を形成する。 それでもランダム・アンサンブルを定義するために使用することはできるが、準確率分布の否定性に関連するサンプリング・オーバーヘッドがあり、これは次の量によって特徴づけられる\n", "\n", "$$\n", "\\gamma = \\sum_{i} \\lvert \\eta_i \\rvert \\geq 1.\n", "$$\n", "\n", "サンプリングオーバーヘッドとは、期待値を所定の精度で推定するために必要なショット数に対する乗算係数であり、理想的な回路から必要とされるショット数と比較される。 これは、 $\\gamma$ に対して2次関数的にスケールする。\n", "\n", "PEC を有効にするには、 [IBM Quantum](/docs/api/qiskit-ibm-runtime/options-resilience-options-v2) の「Estimator の計算耐障害性オプション」で `True` を に `pec_mitigation` 設定します。\n", "PEC のオプションについては、 [こちらで](/docs/api/qiskit-ibm-runtime/options-pec-options)説明しています。 サンプリングのオーバーヘッドの上限は、オプション `max_overhead` を使用して設定できます。 サンプリングのオーバーヘッドを制限すると、結果の精度が要求された精度を超えてしまう可能性があることに注意してください。 のデフォルト値は 100 `max_overhead` です。\n", "\n", "次のコードセルでは、PECを有効にし、Estimatorのオプションを設定 `max_overhead` する方法を示しています。\n", "\n" ] }, { "cell_type": "code", "execution_count": 9, "id": "f4ae33f2-3a4d-4869-952b-e5ebd9e69efc", "metadata": {}, "outputs": [], "source": [ "estimator = Estimator(mode=backend)\n", "estimator.options.resilience.pec_mitigation = True\n", "estimator.options.resilience.pec.max_overhead = 100" ] }, { "cell_type": "markdown", "id": "b9b598a2-03c0-4ff9-ac47-d187693ac6e2", "metadata": {}, "source": [ "\n", "\n", "## 次のステップ\n", "\n", "* エラー軽減オプションとEstimatorプリミティブの組み合わせについては、 [チュートリアルを](/docs/tutorials/combine-error-mitigation-techniques)ご覧ください。\n", "* [Estimator を使用してノイズ管理を設定します](/docs/guides/estimator-noise-management)。\n", "* [Sampler を使用してノイズ管理を設定します](/docs/guides/sampler-noise-management)。\n", "* IBM Quantum プリミティブのその他の[オプション](/docs/guides/runtime-options-overview)を確認してください。\n", "* ジョブ[を実行](/docs/guides/execution-modes)するモードを決定します。\n", "\n" ] }, { "cell_type": "markdown", "metadata": {}, "id": "a1b8767d", "source": "© IBM Corp., 2017-2026" } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3" } }, "nbformat": 4, "nbformat_minor": 5 }