{ "cells": [ { "attachments": {}, "cell_type": "markdown", "id": "7dc25056-2c3c-435a-a827-1930f8bb768f", "metadata": { "gloss": { "ansatz": { "text": "Combination of a reference operator and a variational form, to describe the search space that we're exploring.", "title": "Ansatz" }, "variational": { "text": "Layers of parameterized gates that are repeated a certain number of times, and the parameters of the gates are optimized during the algorithm to minimize the cost function.", "title": "Variational form" } } }, "source": [ "---\n", "title: \"アンザッツと変分形式\"\n", "description: \"このレッスンでは、アンザッツとは何か、そしてそれが基準状態とどのように異なるかを説明します。 アンザッツの例が示される。\"\n", "---\n", "\n", "{/* cspell:ignore ZZII IZZI ZIIZ IZIZ IIZZ ansä */}\n", "\n", "\n", "\n", "# アプローチと変分形式\n", "\n", "すべての変分アルゴリズムの核心は、状態間の差異を分析するという重要なアイデアにある。これらの状態は、パラメータや変数の集合から、何らかの適切な写像（例えば、連続的で微分可能なもの）を通じて関連付けられるのが便利である。\n", "\n", "まず、パラメータ化された回路を手作業で構築する方法を探る。 これらの回路を使って、変分アルゴリズムが探索するパラメータ化された状態のコレクションを表す へんけいけいしき。 次に、この変分形式を基準状態に適用して、 アナザッツ。\n", "\n", "また、この探索空間を探索する際に、スピードと正確さをどのようにトレードオフするかについても探求する。\n", "\n", "![ヒューリスティック・アンサゼッツと問題特異的アンサゼッツを含む、アンサゼッツ議論の主要な構成要素を示す図。](https://quantum.cloud.ibm.com/learning/images/courses/variational-algorithm-design/ansaetze-and-variational-forms/ansatz-workflow.svg)\n", "\n" ] }, { "attachments": {}, "cell_type": "markdown", "id": "5ebac8e8-6451-46e5-8a01-099f84296a17", "metadata": {}, "source": [ "\n", "\n", "## パラメータ化量子回路\n", "\n", "変分アルゴリズムは、 $|\\psi(\\vec{\\theta})\\rangle$ の量子状態を探索し、比較することで動作する。 の量子状態は、 $k$ の有限パラメータセット $\\vec{\\theta} = (\\theta^0, \\ldots, \\theta^{k-1})$ に依存する。これらの状態は、パラメトライズされた量子回路を用いて準備することができ、ゲートは調整可能なパラメータで定義される。 このパラメトリックな回路は、特定の角度を縛ることなく作ることができる：\n", "\n" ] }, { "cell_type": "code", "execution_count": 1, "id": "4c564fab-ed65-4949-b27f-658b0f1ed621", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "\"Output" ] }, "execution_count": 1, "metadata": {}, "output_type": "execute_result" } ], "source": [ "from qiskit.circuit import QuantumCircuit, Parameter\n", "\n", "theta = Parameter(\"θ\")\n", "\n", "qc = QuantumCircuit(3)\n", "qc.rx(theta, 0)\n", "qc.cx(0, 1)\n", "qc.x(2)\n", "\n", "qc.draw(\"mpl\")" ] }, { "cell_type": "code", "execution_count": 2, "id": "4096d161-2916-472b-8abc-d4a6b78955cf", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "\"Output" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/plain": [ "\"Output" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "from math import pi\n", "\n", "angle_list = [pi / 3, pi / 2]\n", "circuits = [qc.assign_parameters({theta: angle}) for angle in angle_list]\n", "\n", "for circuit in circuits:\n", " display(circuit.draw(\"mpl\"))" ] }, { "attachments": {}, "cell_type": "markdown", "id": "f01221c4-f7ee-496a-a664-4b3abc0748c2", "metadata": {}, "source": [ "\n", "\n", "## 変分形式とアンザッツ\n", "\n", "参照状態 $|\\rho\\rangle$ から目標状態 $|\\psi(\\vec\\theta)\\rangle$ へ反復的に最適化するためには、変分形式 $U_V(\\vec{\\theta})$ を定義する必要がある。これは、変分アルゴリズムが探索するパラメトライズされた状態のコレクションを表す：\n", "\n", "$$\n", "\\begin{aligned}\n", "|0\\rangle \\xrightarrow{U_R} U_R|0\\rangle\n", "\n", "& = |\\rho\\rangle \\xrightarrow{U_V(\\vec{\\theta})} U_A(\\vec{\\theta})|0\\rangle \\\\[1mm]\n", "\n", "& = U_V(\\vec{\\theta})U_R|0\\rangle \\\\[1mm]\n", "\n", "& = U_V(\\vec{\\theta})|\\rho\\rangle \\\\[1mm]\n", "\n", "& = |\\psi(\\vec{\\theta})\\rangle \\\\[1mm]\n", "\n", "\\end{aligned}\n", "$$\n", "\n", "パラメータ化された状態は、いかなるパラメータにも依存しない参照状態 $|\\rho\\rangle$ と、常にパラメータに依存する変分形式 $U_V(\\vec{\\theta})$ の両方に依存することに注意。 この2つの半分の組み合わせをアナザッツと呼ぶ。 $U_A(\\vec\\theta) := U_V(\\vec\\theta)U_R$。\n", "\n", "変分アルゴリズムが探索するパラメトライズされた状態のコレクションを表現するためにアナザッツを構築するとき、次元性という重要な問題に気づく。 $n$ -qubit系（つまりヒルベルト空間）は、配置空間内に膨大な数の異なる量子状態を持つ。 それを完全に探求するには、扱いにくい数のパラメーターが必要になる。 量的には、その次元数は $D = 2^{2n}$ である。さらに悪いことに、探索アルゴリズムなどの実行複雑度は、この次元数に応じて指数関数的に増大する。\n", "\n", "この欠点に対処するために、最も関連性の高い状態のみを探索するように、変分形式に合理的な制約を課すのが一般的である。 効率的なトランケート・アンサッツを見つけることは活発な研究分野であるが、ここでは2つの一般的なデザインを取り上げる。\n", "\n" ] }, { "attachments": {}, "cell_type": "markdown", "id": "6d6caf9e-4eb0-4ae1-b9be-6fcfacd23b39", "metadata": {}, "source": [ "\n", "\n", "## ヒューリスティックなアプローチとトレードオフ\n", "\n", "次元数を制限するのに役立つような、あなたの特定の問題に関する情報がない場合は、 $2^{2n}$ より少ないパラメータで、パラメータ化された回路の任意のファミリーを試すことができます。 しかし、考慮すべきトレードオフもある：\n", "\n", "* **スピード** ：探索空間を縮小することで、アルゴリズムをより高速に実行できる。\n", "* **正確さ** ：スペースを小さくすると、問題の実際の解が除外され、最適な解が得られなくなる恐れがある。\n", "* **ノイズ** ：より深い回路はノイズの影響を受けるので、アンサッツの接続性、ゲート、ゲートの忠実度を実験する必要がある。\n", "\n", "パラメータが多ければ多いほど、正確な結果が得られる可能性は高くなるが、アルゴリズムの実行には時間がかかる。\n", "\n" ] }, { "attachments": {}, "cell_type": "markdown", "id": "136b5311-f52b-47a0-a7d2-3f1255b0270e", "metadata": {}, "source": [ "\n", "\n", "### N-局所回路\n", "\n", "ヒューリスティック・アンセッツェの最も広く使われている例のひとつが [Nローカル回路](/docs/api/qiskit/qiskit.circuit.library.NLocal)である：\n", "\n", "* **効率的な実装** ：N-localアナザッツは、通常、少数の物理量子ビットを用いて量子コンピュータ上で効率的に実装できる単純なローカルゲートで構成される。 これにより、量子回路の構築と最適化が容易になる。\n", "* **重要な相関を捉える** N-localアナザッツは、ゲート数が少なくても、量子系における量子ビット間の重要な相関を捉えることができる。 これは、ローカルゲートが隣接する量子ビットに作用し、量子ビット間にエンタングルメントを生じさせることができるためで、複雑な量子系をシミュレートする上で重要になる。\n", "\n", "これらの回路は、以下のように交互に1回以上繰り返される回転層とエンタングルメント層で構成されている：\n", "\n", "* 各層は、最大でも $N$ のサイズのゲートで形成される。 $N$ は量子ビット数より小さくなければならない。\n", "* 回転レイヤーの場合、ゲートは互いに積み重ねられる。 のような標準的な回転操作を使うことができる。 [`RX`](/docs/api/qiskit/qiskit.circuit.library.RXGate) または [`CRZ`](/docs/api/qiskit/qiskit.circuit.library.CRZGate).\n", "* エンタングルメント層には、 [`Toffoli` ゲート](/docs/api/qiskit/qiskit.circuit.library.CCXGate)や [`CX`](/docs/api/qiskit/qiskit.circuit.library.CXGate) のようなゲートとエンタングルメント戦略を使うことができる。\n", "* どちらのタイプのレイヤーもパラメータ化することもしないこともできるが、少なくともどちらか一方はパラメータを含む必要がある。 そうでなければ、少なくとも1つのパラメータがなければ、バリエーションは生まれない！\n", "* オプションで、回路の最後に回転層を追加することもできる。\n", "\n", "例えば、5量子ビットの回路を作ってみよう。 [`NLocal`](/docs/api/qiskit/qiskit.circuit.library.NLocal) と [`RX`](/docs/api/qiskit/qiskit.circuit.library.RXGate) と [`CRZ`](/docs/api/qiskit/qiskit.circuit.library.CRZGate) ゲート、量子ビットに作用する [`Toffoli` ゲートで](/docs/api/qiskit/qiskit.circuit.library.CCXGate)形成されるエンタングルメント・ブロック $[0,1,2]$、 $[0,2,3]$、 $[4,2,1]$、 $[3,1,0]$、 $2$ の各層の繰り返し。\n", "\n" ] }, { "cell_type": "code", "execution_count": 3, "id": "c28d7b53-fbd7-42f5-a619-5a6b0b44e9fc", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "\"Output" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "from qiskit.circuit.library import NLocal, CCXGate, CRZGate, RXGate\n", "from qiskit.circuit import Parameter\n", "\n", "theta = Parameter(\"θ\")\n", "ansatz = NLocal(\n", " num_qubits=5,\n", " rotation_blocks=[RXGate(theta), CRZGate(theta)],\n", " entanglement_blocks=CCXGate(),\n", " entanglement=[[0, 1, 2], [0, 2, 3], [4, 2, 1], [3, 1, 0]],\n", " reps=2,\n", " insert_barriers=True,\n", ")\n", "ansatz.decompose().draw(\"mpl\")" ] }, { "attachments": {}, "cell_type": "markdown", "id": "d293b8cc-1512-4a65-a3f7-94a979d1736b", "metadata": {}, "source": [ "上記の例では、最大のゲートはトッフォリ・ゲートで、これは3つの量子ビットに作用し、回路を $3$ -localにする。 $N$ -ローカル回路の最も一般的なタイプは、シングルqubitローテーションゲートと $2$ -qubitエンタングルメントゲートを持つ [$2$](/docs/api/qiskit/qiskit.circuit.library.TwoLocal) -ローカル回路である。\n", "\n", "Qiskitのクラスを使って $2$ -local回路を作ってみよう。 [`TwoLocal`](/docs/api/qiskit/qiskit.circuit.library.TwoLocal) クラスを使ってみましょう。 構文は `NLocal` と同じだが、いくつか違いがある。 例えば、 `RX`、 `RZ`、 `CNOT` などのほとんどのゲートは、ゲートをインポートしたり、 `Parameter` インスタンスを作成したりすることなく、文字列として渡すことができる。\n", "\n" ] }, { "cell_type": "code", "execution_count": 4, "id": "c48521ee-db12-449e-af9f-7d8ec0fa7890", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "\"Output" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "from qiskit.circuit.library import TwoLocal\n", "\n", "ansatz = TwoLocal(\n", " num_qubits=5,\n", " rotation_blocks=[\"rx\", \"rz\"],\n", " entanglement_blocks=\"cx\",\n", " entanglement=\"linear\",\n", " reps=2,\n", " insert_barriers=True,\n", ")\n", "ansatz.decompose().draw(\"mpl\")" ] }, { "attachments": {}, "cell_type": "markdown", "id": "a6e80571-9bb4-4e33-9581-2c48cad1a566", "metadata": {}, "source": [ "この場合、それぞれの量子ビットが次の量子ビットとエンタングルされる線形エンタングルメント分布を用いた。 その他の戦略については、 [`TwoLocal` のドキュメントを](/docs/api/qiskit/qiskit.circuit.library.TwoLocal)参照のこと。\n", "\n" ] }, { "attachments": {}, "cell_type": "markdown", "id": "65d39e58-e54b-4ea3-95bf-ee68d9b89832", "metadata": {}, "source": [ "\n", "\n", "### 効率的な SU2\n", "\n", "[`efficient_su2`](/docs/api/qiskit/qiskit.circuit.library.EfficientSU2) は、SU(2)と `CX` エンタングルメントにまたがる1量子ビット演算のレイヤーで構成されるハードウェア効率の高い回路である。 これは、変分量子アルゴリズムのための試行波動関数の準備や、機械学習のための分類回路として使用できる発見的パターンである。\n", "\n" ] }, { "cell_type": "code", "execution_count": null, "id": "5abc6f15-d7f5-466a-9659-78a11a2c5270", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "\"Output" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "from qiskit.circuit.library import efficient_su2\n", "\n", "ansatz = efficient_su2(4, su2_gates=[\"rx\", \"y\"], entanglement=\"linear\", reps=1)\n", "ansatz.decompose().draw(\"mpl\")" ] }, { "attachments": {}, "cell_type": "markdown", "id": "14d01a25-2e2a-47e5-9c50-910599c59502", "metadata": {}, "source": [ "\n", "\n", "## 問題特異的なアプローチ\n", "\n", "発見的でハードウェア的に効率的な解法は、素朴な方法で問題を解決するのに役立つが、問題固有の知識を利用することで、回路の探索空間を特定のタイプに限定することができる。 これにより、検索プロセスの精度を落とすことなく、スピードを上げることができる。\n", "\n" ] }, { "attachments": {}, "cell_type": "markdown", "id": "ae376a5f-004c-44af-a6ba-4e4419e350c3", "metadata": {}, "source": [ "\n", "\n", "### 最適化\n", "\n", "max-cut問題では、異なるグループのノード間の辺の数が最大になるようにグラフのノードを分割したい。 左の0番目のノードは、右の残りのノードからカットによって分離されるべきである。\n", "\n" ] }, { "cell_type": "code", "execution_count": 6, "id": "ef9a5e1d-f70a-4f59-9112-2d8f71b253ba", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "\"Output" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "import rustworkx as rx\n", "from rustworkx.visualization import mpl_draw\n", "\n", "n = 4\n", "G = rx.PyGraph()\n", "G.add_nodes_from(range(n))\n", "# The edge syntax is (start, end, weight)\n", "edges = [(0, 1, 1.0), (0, 2, 1.0), (0, 3, 1.0), (1, 2, 1.0), (2, 3, 1.0)]\n", "G.add_edges_from(edges)\n", "\n", "mpl_draw(\n", " G, pos=rx.shell_layout(G), with_labels=True, edge_labels=str, node_color=\"#1192E8\"\n", ")" ] }, { "cell_type": "markdown", "id": "65d72527-77b1-4eb4-93e5-762a5f55b838", "metadata": {}, "source": [ "QAOAアルゴリズムをマックスカット問題に利用するためには、オペレータの最小期待値が2つの異なるグループのノード間の最大エッジ数に対応するような方法でコストをエンコードするパウリ・ハミルトニアンが必要である。\n", "\n", "この単純な例では、演算子はエッジで結ばれたノード上のZ演算子を持つ項の線形結合である（0番目の量子ビットが最も右にあることを思い出してほしい）： $ZZII + IZZI + ZIIZ + IZIZ + IIZZ$.いったん演算子が構築されれば、QAOAアルゴリズムのansatzは、Qiskit回路ライブラリの `QAOAAnsatz` 回路を使うことで簡単に構築できます。\n", "\n" ] }, { "cell_type": "code", "execution_count": 7, "id": "69599589-689e-4e0f-8135-937aedd9430d", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "\"Output" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# Pre-defined ansatz circuit, operator class and visualization tools\n", "from qiskit.circuit.library import QAOAAnsatz\n", "from qiskit.quantum_info import SparsePauliOp\n", "\n", "# Problem to Hamiltonian operator\n", "hamiltonian = SparsePauliOp.from_list(\n", " [(\"ZZII\", 1), (\"IZZI\", 1), (\"ZIIZ\", 1), (\"IZIZ\", 1), (\"IIZZ\", 1)]\n", ")\n", "# QAOA ansatz circuit\n", "ansatz = QAOAAnsatz(hamiltonian, reps=2)\n", "# Draw\n", "ansatz.decompose(reps=3).draw(\"mpl\")" ] }, { "cell_type": "markdown", "id": "23728bd9-ac92-4e3a-a886-2238fc5287ad", "metadata": {}, "source": [ "前の画像は、わかりやすくするために、基本的なゲートのansatzを示している。 しかし、 `reps` の引数を変えるか、decompose メソッドを使わずに回路を描くことで、複数レベルの分解で表現することができる。 例えば、以下の表現は、デフォルトのレップス値（ `reps=1` ）を持つQAOA構造を直接示している。\n", "\n" ] }, { "cell_type": "code", "execution_count": 8, "id": "9724a36b-aae6-4d98-aad2-1090f0e1a941", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "\"Output" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "ansatz.decompose(reps=2).draw(\"mpl\")" ] }, { "attachments": {}, "cell_type": "markdown", "id": "6a8aa7fc-3f80-46b2-b520-b6006d75a404", "metadata": {}, "source": [ "\n", "\n", "### 量子 Machine Learning\n", "\n", "機械学習では、データを2つ以上のカテゴリーに分類することが一般的なアプリケーションである。 これには、古典的な特徴ベクトルを量子ヒルベルト空間にマッピングする特徴マップにデータポイントをエンコードすることが含まれる。 古典的なシミュレーションが困難なパラメータ化された量子回路に基づいて量子特徴マップを構築することは、古典的な機械学習アプローチに対する潜在的な優位性を得るための重要なステップであり、現在活発な研究分野である。\n", "\n", "`zz_feature_map` を使ってパラメータ化された回路を作ることができる。 データ点を特徴マップに渡す ( $x$ )ことができ、また別の変分形式をパラメータとして重みを渡す ( $\\theta$ )ことができる。\n", "\n" ] }, { "cell_type": "code", "execution_count": null, "id": "9d3a813b-40cf-4399-b912-714893987c30", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "\"Output" ] }, "execution_count": 9, "metadata": {}, "output_type": "execute_result" } ], "source": [ "from qiskit.circuit.library import zz_feature_map, TwoLocal\n", "\n", "data = [0.1, 0.2]\n", "\n", "zz_feature_map_reference = zz_feature_map(feature_dimension=2, reps=2)\n", "zz_feature_map_reference = zz_feature_map_reference.assign_parameters(data)\n", "\n", "variation_form = TwoLocal(2, [\"ry\", \"rz\"], \"cz\", reps=2)\n", "vqc_ansatz = zz_feature_map_reference.compose(variation_form)\n", "vqc_ansatz.decompose().draw(\"mpl\")" ] }, { "attachments": {}, "cell_type": "markdown", "id": "4a1b37bd-8de9-479b-8453-1d8942492731", "metadata": {}, "source": [ "\n", "\n", "## サマリー\n", "\n", "このレッスンでは、変分形式で探索空間を定義する方法を学んだ：\n", "\n", "* *パラメトライズされた*量子回路で状態を準備。ゲートは調整可能なパラメータで定義される\n", "* スピードと正確さをトレードオフにする解答の組み立て方\n", "* ヒューリスティック分析\n", "* 問題特有の解決策\n", "\n", "高レベルの変分作業負荷は次のようになる：\n", "\n", "![参照状態を準備するユニタリーとアンサッツを準備するユニタリーの2つを示す回路図。](https://quantum.cloud.ibm.com/learning/images/courses/variational-algorithm-design/ansaetze-and-variational-forms/ansatz-circuit.svg)\n", "\n", "各変分パラメータ $\\vec\\theta$、異なる量子状態が生成される。 最適なパラメータを見つけるためには、問題に特化した*コスト関数を*定義して、アサッツのパラメータを繰り返し更新する必要がある。\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": 2 }