{
"cells": [
{
"cell_type": "markdown",
"id": "d2c31ae8",
"metadata": {},
"source": [
"---\n",
"title: \"期待値推定のためのワイヤー切断\"\n",
"description: \"ワイヤカットを用いて回路を多くの小さなサブ回路に分割する。\"\n",
"---\n",
"\n",
"{/* cspell:ignore edgecolor Cutqc forall fontsize */}\n",
"\n",
"\n",
"\n",
"# 期待値推定のためのワイヤー切断\n",
"\n",
"*推定実行時間：Heronプロセッサで22秒（注：これはあくまで推定値です。 （実行時間は環境によって異なる場合があります。）*\n",
"\n"
]
},
{
"cell_type": "markdown",
"id": "8bf80006",
"metadata": {},
"source": [
"\n",
"\n",
"## 学習成果\n",
"\n",
"このチュートリアルを終えた後、ユーザーは以下の点を理解できるようになります：\n",
"\n",
"* 大きな回路を小さなサブ回路に分割し、それによってノイズの影響を低減する方法 [`qiskit-addon-cutting`](https://github.com/Qiskit/qiskit-addon-cutting)\n",
"\n",
"\n",
"\n",
"## 前提条件\n",
"\n",
"このチュートリアルを進める前に、以下のトピックについてあらかじめ理解しておいていただくことをお勧めします：\n",
"\n",
"* このワークフローで使用される「 [Sampler](/docs/api/qiskit-ibm-runtime/sampler-v2) 」プリミティブを使用する\n",
"\n",
"\n",
"\n",
"## 背景\n",
"\n",
"回路編みとは、回路を、より少ないゲートや量子ビットで構成される複数の小さなサブ回路に分割するさまざまな手法を総称する用語である。 各サブ回路は独立して実行することができ、最終結果は各サブ回路の出力に対して古典的な後処理を行うことで得られる。 この手法は、 [Qiskitの](https://qiskit.github.io/qiskit-addon-cutting/index.html) 「Circuit cutting」アドオンで利用可能です。詳細については、 [ドキュメント](https://qiskit.github.io/qiskit-addon-cutting/explanation/index.html)やその他の[入門資料](https://qiskit.github.io/qiskit-addon-cutting/tutorials/index.html)をご参照ください。\n",
"\n",
"このチュートリアルでは、配線に沿って回路を分割する 「**ワイヤーカッティング** 」と呼ばれる手法に焦点を当てています [\\[1\\], \\[2\\]](#references)。 なお、従来の回路では、分割点での出力が決定論的に決定でき、0か1のいずれかとなるため、分割は簡単である。 しかし、カットの時点における量子ビットの状態は、一般に混合状態である。 したがって、各サブ回路については、異なる基底（通常はパウリ基底 [\\[3\\]、\\[4\\]](#references) のようなトモグラフィー的に完全な基底）を用いて複数回測定を行い、それに応じて固有状態に準備する必要がある。 下図（出典： [\\[7\\]](#references) ）は、4量子ビットのGHZ状態を3つのサブ回路に分割するためのワイヤーカットの例を示している。 ここで、 $M_j$ は基底の集合（通常は Pauli X, Y, Z）を表し、 $P_i$ は固有状態の集合（通常は $|0\\rangle$, $|1\\rangle$, $|+\\rangle$, $|+i\\rangle$ ）を表す。\n",
"\n",
"\n",
"\n",
"\n",
"各サブ回路は量子ビットやゲートの数が少ないため、ノイズの影響を受けにくいと予想される。 このチュートリアルでは、この手法を用いてシステムのノイズを効果的に低減できる例を紹介します。\n",
"\n"
]
},
{
"cell_type": "markdown",
"id": "55b94021",
"metadata": {},
"source": [
"\n",
"\n",
"## 要件\n",
"\n",
"このチュートリアルを始める前に、以下のものがインストールされていることを確認してください：\n",
"\n",
"* Qiskit SDK v2.0 またはそれ以降、 [可視化](/docs/api/qiskit/visualization)サポート付き\n",
"* Qiskit Runtime v0.22 またはそれ以降 ( `pip install qiskit-ibm-runtime` )\n",
"* 回路作成用Qiskitアドオン v0.10.0 以降 (`pip install qiskit-addon-cutting`)\n",
"* Qiskit アドオン utils 0.3 以降 (`pip install qiskit-addon-utils`)\n",
"* Qiskit Aer (`pip install qiskit-aer` )\n",
"\n"
]
},
{
"attachments": {},
"cell_type": "markdown",
"id": "7db2e559",
"metadata": {},
"source": [
"\n",
"\n",
"## セットアップ\n",
"\n"
]
},
{
"cell_type": "code",
"execution_count": 1,
"id": "bc380c46",
"metadata": {},
"outputs": [],
"source": [
"import numpy as np\n",
"import matplotlib.pyplot as plt\n",
"\n",
"from qiskit.circuit import Parameter, ParameterVector, QuantumCircuit\n",
"from qiskit.quantum_info import PauliList, SparsePauliOp\n",
"from qiskit.transpiler import generate_preset_pass_manager\n",
"from qiskit_aer import AerSimulator\n",
"from qiskit.result import sampled_expectation_value\n",
"\n",
"from qiskit_addon_cutting.instructions import CutWire\n",
"from qiskit_addon_cutting import (\n",
" cut_wires,\n",
" expand_observables,\n",
" partition_problem,\n",
" generate_cutting_experiments,\n",
" reconstruct_expectation_values,\n",
")\n",
"\n",
"from qiskit_ibm_runtime import QiskitRuntimeService\n",
"from qiskit_ibm_runtime import SamplerV2, Batch"
]
},
{
"cell_type": "markdown",
"id": "431a5bd2-e6ed-471b-ad9e-c4edd27784a8",
"metadata": {},
"source": [
"\n",
"\n",
"## 小規模シミュレータの例\n",
"\n",
"このチュートリアルでは、1次元（ 1D ）の多体局在（MBL）回路を[シミュ](/docs/guides/intro-to-patterns)レートするためのQiskitパターンを実装します。 MBL回路はハードウェア効率に優れた回路であり、 $\\theta$ および $\\vec{\\phi}$ の2つのパラメータによって定数化されます。 $\\theta$ を $0$ に設定し、 $|0\\rangle$ ですべての量子ビットに対して初期状態を準備した場合、各量子ビットサイト $i$ において、 $\\langle Z_i \\rangle$ の理想的な期待値は、 $\\vec{\\phi}$ の値にかかわらず、 $+1$ となります。この回路に関する詳細は、こちらの[記事](https://www.nature.com/articles/s41467-025-57623-x)をご覧ください。\n",
"\n",
"なお、ノイズのないシミュレータでは、回路を切断した場合としない場合で得られる期待値は同じになります。\n",
"\n"
]
},
{
"attachments": {},
"cell_type": "markdown",
"id": "988ee237",
"metadata": {},
"source": [
"\n",
"\n",
"### ステップ1：古典的な入力を量子問題にマッピングする\n",
"\n",
"\n",
"\n",
"#### 1D のMBL回路を組み立てる\n",
"\n",
"まず、 1D のMBL回路を構成するための関数を紹介します。\n",
"\n"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "d2aa988a",
"metadata": {},
"outputs": [],
"source": [
"class MBLChainCircuit(QuantumCircuit):\n",
" def __init__(\n",
" self, num_qubits: int, depth: int, use_cut: bool = False\n",
" ) -> None:\n",
" super().__init__(\n",
" num_qubits, name=f\"MBLChainCircuit<{num_qubits}, {depth}>\"\n",
" )\n",
" evolution = MBLChainEvolution(num_qubits, depth, use_cut)\n",
" self.compose(evolution, inplace=True)\n",
"\n",
"\n",
"class MBLChainEvolution(QuantumCircuit):\n",
" def __init__(self, num_qubits: int, depth: int, use_cut) -> None:\n",
" super().__init__(\n",
" num_qubits, name=f\"MBLChainEvolution<{num_qubits}, {depth}>\"\n",
" )\n",
"\n",
" theta = Parameter(\"θ\")\n",
" phis = ParameterVector(\"φ\", num_qubits)\n",
"\n",
" for layer in range(depth):\n",
" layer_parity = layer % 2\n",
" # print(\"layer parity\", layer_parity)\n",
" for qubit in range(layer_parity, num_qubits - 1, 2):\n",
" # print(qubit)\n",
" self.cz(qubit, qubit + 1)\n",
" self.u(theta, 0, np.pi, qubit)\n",
" self.u(theta, 0, np.pi, qubit + 1)\n",
" if (\n",
" use_cut\n",
" and layer_parity == 0\n",
" and (\n",
" qubit == num_qubits // 2 - 1\n",
" or qubit == num_qubits // 2\n",
" )\n",
" ):\n",
" self.append(CutWire(), [num_qubits // 2])\n",
" if use_cut and layer < depth - 1 and layer_parity == 1:\n",
" if qubit == num_qubits // 2:\n",
" self.append(CutWire(), [qubit])\n",
" for qubit in range(num_qubits):\n",
" self.p(phis[qubit], qubit)"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "a3debf65-06df-4277-933e-14b6f6170756",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
""
]
},
"execution_count": 3,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"num_qubits = 10\n",
"depth = 2\n",
"mbl = MBLChainCircuit(num_qubits, depth)\n",
"mbl.draw(\"mpl\", fold=-1)"
]
},
{
"cell_type": "markdown",
"id": "ef13f367",
"metadata": {},
"source": [
"$\\theta = 0$ について、すべての量子ビットにわたる平均期待値 $O = \\frac{1}{n} \\sum_i Z_i$ を計算する。 $\\langle Z_i \\rangle = 1$ の理想的な期待値は $\\forall$ であり、 $i$ であるため、 $O$ の理想的な期待値も $1$ となる。パラメータ $\\phi$ はランダムに選択される。\n",
"\n"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "b0f1e5fd",
"metadata": {},
"outputs": [],
"source": [
"np.random.seed(42)\n",
"phis = list(np.random.rand(mbl.num_parameters - 1))\n",
"theta = [0]\n",
"params = theta + phis"
]
},
{
"cell_type": "markdown",
"id": "a81f0346",
"metadata": {},
"source": [
"回路を分割するには、分割したい箇所に **CutWire** を挿入して注釈を付ける必要があります。 このチュートリアルでは、等分分割を採用します。 MBL回路は、関数内の設定 `use_cut=True` により、 $\\frac{n}{2}$ 個の量子ビットの後に注釈が適切に挿入されるように設計されています。ここで、 $n$ は元の回路に含まれる量子ビットの数です。 また、ランダムに生成されたパラメータを回路に割り当てました。\n",
"\n"
]
},
{
"cell_type": "code",
"execution_count": 5,
"id": "5208e0a8",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
""
]
},
"execution_count": 5,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"mbl_cut = MBLChainCircuit(num_qubits, depth, use_cut=True)\n",
"mbl_cut.assign_parameters(params, inplace=True)\n",
"mbl_cut.draw(\"mpl\", fold=-1)"
]
},
{
"cell_type": "markdown",
"id": "ac6f36e3",
"metadata": {},
"source": [
"\n",
"\n",
"### ステップ2：量子ハードウェア実行に向けた問題の最適化\n",
"\n",
"\n",
"\n",
"#### 回路を小さなサブ回路に分割する\n",
"\n",
"[`qiskit-addon-cutting`](https://qiskit.github.io/qiskit-addon-cutting/)ここで、を使用して回路を2つの小さな部分回路に分割します。 `qiskit-addon-cutting` クビットの数を適切に調整することで、ワイヤの切断位置を分割する仮想 `Move` ゲートを追加する。 それでは、この仮想ゲートを使って回路を作成しましょう。 1本の配線が切断されているため、関連する量子ビットの数は1つ増えます。\n",
"\n"
]
},
{
"cell_type": "code",
"execution_count": 6,
"id": "1834cb22",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
""
]
},
"execution_count": 6,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"mbl_move = cut_wires(mbl_cut)\n",
"mbl_move.draw(\"mpl\", fold=-1)"
]
},
{
"cell_type": "markdown",
"id": "1df16e6b",
"metadata": {},
"source": [
"\n",
"\n",
"#### オブザーバブルを構築および拡張する\n",
"\n",
"前述の定義に従えば、観測量は各量子ビットにおける $Z$ の平均値となる。 しかし、仮想 `Move` ゲートを挿入すると、回路における有効な量子ビット数は増加する。 また、量子ビット数のこの変化を考慮に入れるため、観測量をそれに応じて展開する必要があります。 仮想 `Move` ゲートのために追加された余分な量子ビットに対して、観測演算子は常に自明な作用（ $I$ のように）しか及ぼさないことに注意してください。\n",
"\n"
]
},
{
"cell_type": "code",
"execution_count": 7,
"id": "3074b173",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"PauliList(['ZIIIIIIIII', 'IZIIIIIIII', 'IIZIIIIIII', 'IIIZIIIIII',\n",
" 'IIIIZIIIII', 'IIIIIZIIII', 'IIIIIIZIII', 'IIIIIIIZII',\n",
" 'IIIIIIIIZI', 'IIIIIIIIIZ'])"
]
},
"execution_count": 7,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"observable = PauliList(\n",
" [\"I\" * i + \"Z\" + \"I\" * (num_qubits - i - 1) for i in range(num_qubits)]\n",
")\n",
"observable"
]
},
{
"cell_type": "code",
"execution_count": 8,
"id": "32b7081b",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"PauliList(['ZIIIIIIIIII', 'IZIIIIIIIII', 'IIZIIIIIIII', 'IIIZIIIIIII',\n",
" 'IIIIZIIIIII', 'IIIIIIZIIII', 'IIIIIIIZIII', 'IIIIIIIIZII',\n",
" 'IIIIIIIIIZI', 'IIIIIIIIIIZ'])"
]
},
"execution_count": 8,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"new_obs = expand_observables(observable, mbl, mbl_move)\n",
"new_obs"
]
},
{
"cell_type": "markdown",
"id": "0ad0886b",
"metadata": {},
"source": [
"これで、回路をゲート `Move` に沿って分割することができ、各サブ回路に対応する元の観測量の部分であるサブ観測量とともに、サブ回路が得られる。\n",
"\n"
]
},
{
"cell_type": "code",
"execution_count": 12,
"id": "c8894a06",
"metadata": {},
"outputs": [],
"source": [
"partitioned_problem = partition_problem(circuit=mbl_move, observables=new_obs)\n",
"subcircuits = partitioned_problem.subcircuits\n",
"subobservables = partitioned_problem.subobservables"
]
},
{
"cell_type": "markdown",
"id": "8f47ab33",
"metadata": {},
"source": [
"ここでは、2つのサブ回路を図示します：\n",
"\n"
]
},
{
"cell_type": "code",
"execution_count": 13,
"id": "1b0e779f",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
""
]
},
"execution_count": 13,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"subcircuits[0].draw(\"mpl\", fold=-1)"
]
},
{
"cell_type": "code",
"execution_count": 14,
"id": "3c802f28",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
""
]
},
"execution_count": 14,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"subcircuits[1].draw(\"mpl\", fold=-1)"
]
},
{
"cell_type": "markdown",
"id": "0e4c03fb",
"metadata": {},
"source": [
"演算 `Move` を用いて観測可能領域を拡張するには、データ構造が必要 `PauliList` となる。 元の回路の期待値を再構築するには、観測量が の `SparsePauliOp` 形式である必要があります。\n",
"\n"
]
},
{
"cell_type": "code",
"execution_count": 15,
"id": "7dc51220",
"metadata": {},
"outputs": [],
"source": [
"M_z = SparsePauliOp(\n",
" [\"I\" * i + \"Z\" + \"I\" * (num_qubits - i - 1) for i in range(num_qubits)],\n",
" coeffs=[1 / num_qubits] * num_qubits,\n",
")"
]
},
{
"cell_type": "markdown",
"id": "2b5a5746",
"metadata": {},
"source": [
"前述の通り、各カットにおいて、上流回路はパウリ基底で測定されなければならず、下流回路はその基底の固有状態へと準備されなければならない。 この関 `generate_cutting_experiments` 数は、再構成に必要なすべての回路と、各回路に関連付けられた係数を生成します。 詳細については[、本論文](https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.125.150504)をご覧ください。\n",
"\n"
]
},
{
"cell_type": "code",
"execution_count": 16,
"id": "648519f0",
"metadata": {},
"outputs": [],
"source": [
"subexperiments, coefficients = generate_cutting_experiments(\n",
" circuits=subcircuits,\n",
" observables=subobservables,\n",
" num_samples=np.inf,\n",
")"
]
},
{
"cell_type": "markdown",
"id": "43f58cfb",
"metadata": {},
"source": [
"\n",
"\n",
"#### 回路をバックエンドにトランスパイルする\n",
"\n",
"シミュレーションのみを行う最初の例として、回路をバックエンドの基本ゲートセットに変換します：\n",
"\n"
]
},
{
"cell_type": "code",
"execution_count": 17,
"id": "29d71cd3",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"\n"
]
}
],
"source": [
"service = QiskitRuntimeService()\n",
"backend = service.least_busy(\n",
" operational=True, simulator=False, min_num_qubits=133\n",
")\n",
"\n",
"print(backend)"
]
},
{
"attachments": {},
"cell_type": "markdown",
"id": "b4d480b3",
"metadata": {},
"source": [
"\n",
"\n",
"### ステップ3: `Qiskit primitives`を使用して実行する\n",
"\n",
"それでは、各サブ実験を実行してください：\n",
"\n"
]
},
{
"cell_type": "code",
"execution_count": 20,
"id": "05fabb27",
"metadata": {},
"outputs": [],
"source": [
"pm_basis = generate_preset_pass_manager(\n",
" optimization_level=2, basis_gates=backend.configuration().basis_gates\n",
")\n",
"basis_subexperiments = {\n",
" label: pm_basis.run(partition_subexpts)\n",
" for label, partition_subexpts in subexperiments.items()\n",
"}"
]
},
{
"cell_type": "code",
"execution_count": 21,
"id": "b22a1b00",
"metadata": {},
"outputs": [],
"source": [
"sampler = SamplerV2(mode=AerSimulator())\n",
"jobs = {\n",
" label: sampler.run(subsystem_subexpts, shots=2**12)\n",
" for label, subsystem_subexpts in basis_subexperiments.items()\n",
"}"
]
},
{
"cell_type": "markdown",
"id": "50b94af2",
"metadata": {},
"source": [
"\n",
"\n",
"### ステップ4：後処理を行い、結果を希望の古典形式で返す\n",
"\n",
"ここで、各サブ実験の実行結果を取得し、未カット回路の期待値を再構築します：\n",
"\n"
]
},
{
"cell_type": "code",
"execution_count": 22,
"id": "31dc35ea-6554-4ca7-9c3b-0b5394c46e4e",
"metadata": {},
"outputs": [],
"source": [
"# Retrieve results\n",
"results = {label: job.result() for label, job in jobs.items()}"
]
},
{
"cell_type": "code",
"execution_count": 23,
"id": "4972e4e4",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"np.float64(0.9953821063041687)"
]
},
"execution_count": 23,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"reconstructed_expval_terms = reconstruct_expectation_values(\n",
" results,\n",
" coefficients,\n",
" subobservables,\n",
")\n",
"reconstructed_expval = np.dot(reconstructed_expval_terms, M_z.coeffs).real\n",
"reconstructed_expval"
]
},
{
"cell_type": "code",
"execution_count": 32,
"id": "fb5f955a",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"Text(0, 0.5, '$M_Z$')"
]
},
"execution_count": 32,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"methods = [\n",
" \"Uncut\",\n",
" \"Wire cut\",\n",
"]\n",
"values = [\n",
" 1,\n",
" reconstructed_expval,\n",
"] # since the ideal expectation value in noiseless simulation is +1\n",
"\n",
"ax = plt.gca()\n",
"plt.bar(methods, values, color=\"#a56eff\", width=0.4, edgecolor=\"#8a3ffc\")\n",
"ax.set_ylabel(r\"$M_Z$\", fontsize=12)"
]
},
{
"cell_type": "markdown",
"id": "0d6db390-e7a8-4efe-902c-8d9a312170c6",
"metadata": {},
"source": [
"\n",
"\n",
"## 大規模なハードウェアの例\n",
"\n"
]
},
{
"cell_type": "markdown",
"id": "ae69c5e0-32b1-4f03-ab13-7b95a9acfd25",
"metadata": {},
"source": [
"ここでは、60キュービットのMBL回路におけるワイヤーカッティングを実演します。 未カットの回路およびカットされた回路は、 IBM Quantum® ハードウェア上で実行されます：\n",
"\n"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "37834c72",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"num_qubits = 60\n",
"depth = 2\n",
"\n",
"# construct the circuit\n",
"mbl = MBLChainCircuit(num_qubits, depth)\n",
"\n",
"# create parameters\n",
"phis = list(np.random.rand(mbl.num_parameters - 1))\n",
"theta = [0]\n",
"params = theta + phis\n",
"\n",
"# construct the cut circuit\n",
"mbl_cut = MBLChainCircuit(num_qubits, depth, use_cut=True)\n",
"mbl_cut.assign_parameters(params, inplace=True)\n",
"mbl_move = cut_wires(mbl_cut)\n",
"\n",
"# Define observable and expand to account for the wire cut\n",
"observable = PauliList(\n",
" [\"I\" * i + \"Z\" + \"I\" * (num_qubits - i - 1) for i in range(num_qubits)]\n",
")\n",
"new_obs = expand_observables(observable, mbl, mbl_move)\n",
"\n",
"# Construct a SparsePauliOp version of the observable for later use in reconstruction\n",
"M_z = SparsePauliOp(\n",
" [\"I\" * i + \"Z\" + \"I\" * (num_qubits - i - 1) for i in range(num_qubits)],\n",
" coeffs=[1 / num_qubits] * num_qubits,\n",
")\n",
"\n",
"# Partition the circuit and get subcircuits and subobservables\n",
"partitioned_problem = partition_problem(circuit=mbl_move, observables=new_obs)\n",
"subcircuits = partitioned_problem.subcircuits\n",
"subobservables = partitioned_problem.subobservables\n",
"\n",
"# Obtain subexperiments and coefficients\n",
"subexperiments, coefficients = generate_cutting_experiments(\n",
" circuits=subcircuits,\n",
" observables=subobservables,\n",
" num_samples=np.inf,\n",
")\n",
"\n",
"# Transpile the subexperiments to the backend\n",
"pm = generate_preset_pass_manager(optimization_level=2, backend=backend)\n",
"isa_subexperiments = {\n",
" label: pm.run(partition_subexpts)\n",
" for label, partition_subexpts in subexperiments.items()\n",
"}\n",
"\n",
"# Execute the subexperiments and retrieve results\n",
"with Batch(backend=backend) as batch:\n",
" sampler = SamplerV2(mode=batch)\n",
" sampler.options.environment.job_tags = [\"TUT_WC\"]\n",
" jobs = {\n",
" label: sampler.run(subsystem_subexpts, shots=2**12)\n",
" for label, subsystem_subexpts in isa_subexperiments.items()\n",
" }\n",
"results = {label: job.result() for label, job in jobs.items()}\n",
"\n",
"# Reconstruct the expectation value of the original observable\n",
"reconstructed_expval_terms = reconstruct_expectation_values(\n",
" results,\n",
" coefficients,\n",
" subobservables,\n",
")\n",
"reconstructed_expval = np.dot(reconstructed_expval_terms, M_z.coeffs).real\n",
"\n",
"# Compute the uncut circuit to obtain the noisy expectation value for comparison\n",
"sampler = SamplerV2(mode=backend)\n",
"sampler.options.environment.job_tags = [\"TUT_WC\"]\n",
"\n",
"if mbl.num_clbits == 0:\n",
" mbl.measure_all()\n",
"isa_mbl = pm.run(mbl)\n",
"\n",
"pub = (isa_mbl, params)\n",
"uncut_job = sampler.run([pub])\n",
"\n",
"uncut_counts = uncut_job.result()[0].data.meas.get_counts()\n",
"uncut_expval = sampled_expectation_value(uncut_counts, M_z)\n",
"\n",
"# visualize the results\n",
"ax = plt.gca()\n",
"methods = [\"uncut\", \"cut\"]\n",
"values = [uncut_expval, reconstructed_expval]\n",
"\n",
"plt.bar(methods, values, color=\"#a56eff\", width=0.4, edgecolor=\"#8a3ffc\")\n",
"plt.axhline(y=1, color=\"k\", linestyle=\"--\")\n",
"plt.text(0.3, 0.95, \"Exact result\")\n",
"plt.show()"
]
},
{
"cell_type": "code",
"execution_count": 29,
"id": "1445d099",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"0.9202473958333336"
]
},
"execution_count": 29,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"uncut_expval"
]
},
{
"cell_type": "markdown",
"id": "75f48e6a-c7e4-46f3-9d39-a7a877427a04",
"metadata": {},
"source": [
"\n",
"\n",
"## 次のステップ\n",
"\n",
"\n",
" この作品に興味を持たれた方は、以下の資料もご参照ください：\n",
"\n",
" * [回路切断を伴う周期的境界条件](/docs/tutorials/periodic-boundary-conditions-with-circuit-cutting)\n",
" * [深さ低減のための回路切断](/docs/tutorials/depth-reduction-with-circuit-cutting)\n",
"\n",
"\n"
]
},
{
"cell_type": "markdown",
"id": "23cd3042",
"metadata": {},
"source": [
"\n",
"\n",
"## 参照\n",
"\n",
"\\[1] Peng, T、 ハロー、A． Ozols, M., & Wu, X. (2020). 小型量子コンピュータで大規模量子回路をシミュレートする。 Physical review letters, 125(15), 150504.\n",
"\n",
"\\[2] Tang, W、 トメッシュ、T、 スチャラ、M、 Larson, J., & Martonosi, M. (2021, April). Cutqc：大規模量子回路評価のための小型量子コンピュータの使用。 プログラミング言語とオペレーティングシステムのアーキテクチャサポートに関する第26回ACM国際会議予稿集 (pp. 473-486).\n",
"\n",
"\\[3] パーリン、M. A、 サリーム、Z. H、 スチャラ、M.、オズボーン、J. C. (2021). 最尤トモグラフィによる量子回路切断。 量子情報, 7(1), 64.\n",
"\n",
"\\[4] Majumdar, R., & Wood, C. J. (2022). エラーを軽減した量子回路切断 arXiv preprint arXiv:2211.13431.\n",
"\n",
"\\[5] カレ、T、 マジュムダル、R、 サングル、R、 レイ、A、 Seshadri, P. V., & Simmhan, Y. (2023). 量子古典ワークロードの並列化：分割技術の影響のプロファイリング。 2023 IEEE International Conference on Quantum Computing and Engineering (QCE) (Vol. 1, pp. 990-1000). IEEE。\n",
"\n",
"\\[6] Bhoumik, D、 マジュムダル、R、 Saha, A., & Sur-Kolay, S. (2023). ノイズと時間の最適化を伴う量子回路の分散スケジューリング. arXiv preprint arXiv:2309.06005.\n",
"\n",
"\\[7] Majumdar, R. (2024). 「離散量子計算回路におけるリソースとノイズの効率的な低減」（博士論文、インド統計研究所－コルカタ）。 [https://www.proquest.com/openview/b481def90b1cc80e6b58a77c99e8385c/1?pq-origsite=gscholar\\&cbl=2026366\\&diss=y](https://www.proquest.com/openview/b481def90b1cc80e6b58a77c99e8385c/1?pq-origsite=gscholar\\&cbl=2026366\\&diss=y)\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"
},
"hours": 1.5,
"qpuSeconds": 22
},
"nbformat": 4,
"nbformat_minor": 5
}