{ "cells": [ { "cell_type": "markdown", "id": "f767af9e-cdd6-41f4-9279-9e92fa2dde06", "metadata": {}, "source": [ "---\n", "title: \"Qiskitを用いた不確実性の探求\"\n", "description: \"物理的な測定値が互いにどのように影響し合うかを学び、実際の量子コンピュータを用いて量子力学における基本不確定性を探求する\"\n", "---\n", "\n", "{/* cspell:ignore prodxz delx delz dely prodxy uparrow downarrow */}\n", "\n", "\n", "\n", "# 不確実性の探求\n", "\n" ] }, { "cell_type": "markdown", "id": "132cdf03-164e-4334-aab6-206aacacc25b", "metadata": {}, "source": [ "このQiskit in Classroomsモジュールでは、以下のパッケージがインストールされた Python 環境が必要です：\n", "\n", "* `qiskit` v2.1.0 または新しい\n", "* `qiskit-ibm-runtime` v0.40.1 または新しい\n", "* `qiskit-aer` v0.17.0 または新しい\n", "* `qiskit.visualization`\n", "* `numpy`\n", "* `pylatexenc`\n", "\n", "上記のパッケージをセットアップしてインストールするには、 [Qiskitのインストール](/docs/guides/install-qiskit)ガイドをご覧ください。\n", "実際の量子コンピュータでジョブを実行するには、 IBM Quantum® のアカウントを設定する必要があります。 [IBM Cloud アカウントの設定](/docs/guides/cloud-setup)ガイドの手順に従ってください。\n", "\n", "このモジュールはテストされ、8分のQPU時間を使用した。これはあくまでも目安である。実際の使用状況とは異なる場合があります。時間のかかる2つの計算は、ヘッダーコメントにそのように記されており、QPUの時間が足りない場合はシミュレーターで行うことができる。これらを取り除くと、このモジュールに必要なQPU時間はわずか30秒になります。\n", "\n" ] }, { "cell_type": "code", "execution_count": null, "id": "dbda80f2-62c5-4f25-95cf-c5495bc171b1", "metadata": {}, "outputs": [], "source": [ "# Uncomment and modify this line as needed to install dependencies\n", "#!pip install 'qiskit>=2.1.0' 'qiskit-ibm-runtime>=0.40.1' 'qiskit-aer>=0.17.0' 'numpy' 'pylatexenc'" ] }, { "cell_type": "markdown", "id": "4283087e-cc02-467d-b02c-c1c131ad15f3", "metadata": {}, "source": [ "ケイティ博士（ McCormick ）によるモジュールのウォークスルーを以下でご覧いただくか、 [こちらを](https://www.youtube.com/watch?v=3h3pwrECbb8\\&list=PLOFEBzvs-Vvrs2fuvsuT039ariYPsua3d\\&index=5)クリックして YouTube でご覧ください。\n", "\n", "***\n", "\n", "\n", "\n" ] }, { "cell_type": "markdown", "id": "be385aa2-9080-4c39-aba0-bc35087ced78", "metadata": {}, "source": [ "\n", "\n", "## 概要\n", "\n", "物理の授業以外でも、不確定性原理という言葉を聞いたことがあるだろう。不確実性についての一般的な口語表現は、\"何かを見ることによって、それに影響を与える \"というものだ確かにその通りだ。しかし、より物理的な不確かさの表現方法としては、ある種の物理的な観測値があり、それらは同時に任意の精度で知ることができない非互換性を持っているということである。多くの生徒が最初に出会うのは、 $x$ と $p_x$ という相容れない変数のペアである。これはそれぞれ、 $x$ -軸と呼ばれる1つの軸に沿った位置と、その方向に沿った直線運動量を意味する。これらの変数に対する不確実性の制約は次のようになる。 $\\Delta x \\Delta p_x \\geq \\frac{\\hbar}{2}.$ ここで、 $\\Delta x$ は「 $x$ の不確実性」と呼ばれ、統計学における標準偏差と同じ定義を持っており、次のように定義できる。 $\\Delta x = \\sqrt{\\langle x^2 \\rangle - \\langle x \\rangle^2}.$ $\\Delta p_x$ も同様に定義される。\n", "ここでは、この不確定性関係を導出するのではなく、それが古典的波動に対する我々の理解と一致していることを指摘する。つまり、完全な周波数（ $f$ ）と波長（ $\\lambda$ ）を持つ波は、完全な正弦波として永遠に続くことになる。量子力学的には、これはド・ブロイの仮説に従って運動量を完全に知ることに相当する。 $\\lambda = h/p$。しかし、波のような粒子の位置（ $where$ ）を知るためには、その粒子を記述する波が、例えば非常に狭いガウスのように、空間においてより鋭く尖ったものにならなければならない。私たちは、このような鋭く尖った波動関数を含むあらゆる連続関数を、波長の異なる正弦波関数のフーリエ級数として表現できることを知っている。しかし、波動関数がより鋭くピークを持つようになると（そして位置がよりよく知られるようになると）、フーリエ級数により多くの項が必要となり、より多くの波長の混合が必要となる（したがって、量子力学的には、より多くの運動量の値が必要となる）。\n", "\n", "もっと簡単に言えば、運動量が明確に定義された状態（空間上の完全な正弦波）は、位置が非常に不確かである。明確な位置を持つ状態（ディラック・デルタ分布のような）は、非常に不確かな運動量を持つ。\n", "\n", "このような相性の悪さを示す変数は他にもある。例えば、粒子のスピンはある軸に沿って明確に投影されるかもしれないが、直交する軸の投影については何もわからない。例えば、状態 $|0\\rangle \\sim |\\uparrow\\rangle$ （量子ビットまたは spin-1/2 粒子の場合）は、 $z$ （量子ビットの場合は1、 spin-1/2 粒子の場合は $\\hbar/2$ ）軸に沿った明確な投影を持つ。しかし、この状態は2つの状態の重ね合わせとして書くことができ、それぞれの状態は $x$ ： $|0\\rangle = \\frac{1}{\\sqrt{2}}(|+\\rangle_x+|-\\rangle_x)$ あるいは等価的に $\\begin{pmatrix} 1 \\\\ 0\\end{pmatrix} = \\frac{1}{\\sqrt{2}}\\left[\\frac{1}{\\sqrt{2}}\\begin{pmatrix} 1 \\\\ 1\\end{pmatrix}+\\frac{1}{\\sqrt{2}}\\begin{pmatrix} 1 \\\\ -1\\end{pmatrix}\\right].$ $|+\\rangle_x$ は、 $|-\\rangle_x$ と同様に、 $x$ へのよく定義された射影を持つ。したがって、 $x$ 軸に沿った状態の投影を指定しても、 $z$ 軸に沿った投影はわからない。また、 $z$ 軸への投影を指定した場合、 $x$ に沿った投影はわかりません。スピンの文脈と量子ビットの文脈でこれを議論する場合、細かい違いがあります。しかし、一般的に言えば、パウリ行列の固有状態には興味深い関係があり、それを探ることができる。\n", "このレッスンでは、これらの相容れない変数の不確定性についての直観を実験的にチェックし、 IBM® 量子コンピュータ上で不確定性関係が成り立つことを検証する。\n", "\n" ] }, { "cell_type": "markdown", "id": "eb8c7a0c-c962-4efd-ab82-905ba7415284", "metadata": {}, "source": [ "\n", "\n", "## 直感の簡単な確認\n", "\n", "この最初の実験とモジュール全体を通して、私たちは「Qiskitパターン」として知られる量子コンピューティングのフレームワークを使用する：\n", "\n", "* ステップ1：古典的入力を量子問題にマップする\n", "* ステップ2：量子実行のための問題の最適化\n", "* ステップ 3: `IBM Quantum` プリミティブを使用して実行する\n", "* ステップ4：後処理と古典的分析\n", "\n", "私たちは通常、これらのステップに従うが、必ずしも明示的にラベルを貼るとは限らない。\n", "\n", "まず、ランタイム・プリミティブを含む必要なパッケージをロードすることから始めよう。また、利用可能な量子コンピュータの中で最も暇なものを選ぶ。\n", "\n", "初回使用時に認証情報を保存するためのコードが以下にあります。ノートブックを自分の環境に保存した後、必ずこの情報をノートブックから削除してください。そうすれば、ノートブックを共有するときにあなたの認証情報が誤って共有されることはありません。詳しいガイダンスについては、 [IBM Cloud アカウントの設定](/docs/guides/initialize-account)および[信頼できない環境でのサービスの初期](/docs/guides/cloud-setup-untrusted)化を参照してください。\n", "\n" ] }, { "cell_type": "code", "execution_count": null, "id": "c26c8f6f-1ed7-4079-ba21-4747d318996f", "metadata": {}, "outputs": [], "source": [ "from numpy import pi\n", "\n", "# Load IBM Quantum Compute Service\n", "from qiskit_ibm_runtime import QiskitRuntimeService\n", "\n", "# Syntax for first saving your token. Delete these lines after saving your credentials.\n", "# QiskitRuntimeService.save_account(channel='ibm_quantum_platform',\n", "# instance = '', token='', overwrite=True, set_as_default=True)\n", "# service = QiskitRuntimeService(channel='ibm_quantum_platform')\n", "\n", "# Load saved credentials\n", "service = QiskitRuntimeService()" ] }, { "cell_type": "code", "execution_count": 2, "id": "7d34071a-6fb1-4a0a-9e66-c9ee40c3767b", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "ibm_sherbrooke\n" ] } ], "source": [ "# Load the Runtime primitive and session\n", "from qiskit_ibm_runtime import (\n", " Batch,\n", " SamplerV2 as Sampler,\n", " EstimatorV2 as Estimator,\n", ")\n", "\n", "# Use the least busy backend\n", "backend = service.least_busy(min_num_qubits=127)\n", "print(backend.name)" ] }, { "cell_type": "markdown", "id": "fcbf932e-45b7-4a50-90e1-59dcb7577dfc", "metadata": {}, "source": [ "もし生徒がレッスン中に量子コンピュータの計算時間を使い果たした場合、以下の行はコメントなしで、上記で選択した量子コンピュータのノイズ動作を部分的に模倣したシミュレータをセットアップするために使用することができます。\n", "\n" ] }, { "cell_type": "code", "execution_count": null, "id": "7265c432-19eb-41ba-940d-a80541fe0218", "metadata": {}, "outputs": [], "source": [ "# Import an estimator, this time from qiskit (we will import from Runtime for real hardware)\n", "from qiskit_aer.primitives import SamplerV2, EstimatorV2\n", "from qiskit_aer.noise import NoiseModel\n", "\n", "# Generate the noise model from the backend properties\n", "noise_model = NoiseModel.from_backend(backend)\n", "\n", "noisy_sampler = SamplerV2(options={\"backend_options\": {\"noise_model\": noise_model}})\n", "noisy_estimator = EstimatorV2(options={\"backend_options\": {\"noise_model\": noise_model}})" ] }, { "cell_type": "markdown", "id": "5f072c75-efb7-4596-8bc6-ccae553684df", "metadata": {}, "source": [ "ある作用素Zの固有状態は、別の作用素Xの固有状態ではないことを思い出してほしい。 $x$、 $z$ の軸に沿って測定を行い、実験的にそれを観察する。 IBM 量子コンピュータは、 $z$ に沿って測定する構造になっているため、 $z$ に沿った測定には、単純に` qc.measure ()` を使用する。しかし、 $x$ に沿って測定するためには、システムを回転させ、 $x$ の軸を測定する方向まで効果的に移動させる必要がある。これはハダマードゲートで達成される。 $y$ に沿った測定にも同様のステップが必要である。ここでは便宜上、必要なステップをまとめた：\n", "\n", "* $z$ に沿って測定する： `qc.measure()`\n", "* $x$ : `qc.h()` に沿って測定する。 `qc.measure()`\n", "* $y$ ： `qc.sdg()`、 `qc.h()`、 `qc.s`。 `qc.measure()`\n", "\n", "\n", "\n", "### ステップ1：古典的な入力を量子問題にマッピングする\n", "\n", "この場合、マッピングのステップは、上述の測定と回転を量子回路で表現するだけである：\n", "\n" ] }, { "cell_type": "code", "execution_count": null, "id": "8baef7f9-78bf-48c8-8359-efd65f943d96", "metadata": {}, "outputs": [ { "data": { "text/plain": [ " $\"Output$ " ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# Step 1: Map\n", "\n", "# Import some general packages\n", "from qiskit import ClassicalRegister, QuantumCircuit, QuantumRegister\n", "\n", "# Define registers\n", "qr = QuantumRegister(1, \"q\")\n", "cr = ClassicalRegister(2, \"c\")\n", "qc = QuantumCircuit(qr, cr)\n", "\n", "# Add a first measurement\n", "qc.measure(qr, cr[0])\n", "qc.barrier()\n", "\n", "# Change basis so that measurements made on quantum computer which normally tell us about z,\n", "# now tell us about x.\n", "qc.h(qr)\n", "\n", "# Add a second measurement\n", "qc.measure(qr, cr[1])\n", "\n", "qc.draw(\"mpl\")" ] }, { "cell_type": "markdown", "id": "4d931ec8-6f8f-405f-a9ce-41724435f3a1", "metadata": {}, "source": [ "\n", "\n", "### ステップ2：量子実行のための問題最適化\n", "\n", "このステップでは、私たちが実行したいオペレーションを、特定の量子コンピュータの機能の観点から表現する。また、我々の問題を量子コンピューターのレイアウトにマッピングすることもできる。\n", "\n" ] }, { "cell_type": "code", "execution_count": 5, "id": "d613ac18-61ff-4e91-9f7d-73887bd2306a", "metadata": {}, "outputs": [], "source": [ "# Step 2: Transpile\n", "from qiskit.transpiler.preset_passmanagers import generate_preset_pass_manager\n", "\n", "target = backend.target\n", "pm = generate_preset_pass_manager(target=target, optimization_level=3)\n", "\n", "qc_isa = pm.run(qc)" ] }, { "cell_type": "markdown", "id": "2f31c2b1-08be-43c7-8af3-5d9ee8a4dab0", "metadata": {}, "source": [ "\n", "\n", "### ステップ 3: `IBM Quantum` プリミティブを使用して実行する\n", "\n", "サンプラーを使用して、測定値に関する統計データを収集することができます。 . `mode = backend`を使用して、実際の量子コンピュータ上で動作するSamplerプリミティブを構築します。ワークフローに応じて他にもモードがありますが、以下ではそのうちの1つを使用します。 Samplerは、「pubs」（プリミティブ・ユニファイド・ブロック）のリストを引数としてrun()メソッドを呼び出すことで使用されます。各パブには最大3つの値が含まれており、これらが組み合わさって、推定器が完了すべき計算単位を定義します。具体的には、回路、観測可能変数、パラメータです。回路の一覧、観測可能変数の一覧、およびパラメータの一覧を指定することもできます。詳細については、「 [PUBの概要](/docs/guides/primitive-input-output#pubs) 」をご覧ください\n", "\n", "私たちは本物の量子コンピューター上で動作させ、本物の量子物理学の実験を実施したい。実際の量子コンピューターで割り当てられた時間を使い果たした場合は、以下の量子コンピューター用のコードをコメントアウトし、シミュレーター上で実行するためのコードをコメントアウトしないようにすればよい。\n", "\n" ] }, { "cell_type": "code", "execution_count": null, "id": "8996bcf0-1848-4eda-b5a1-47190fede630", "metadata": {}, "outputs": [], "source": [ "# Step 3: Run the job on a real quantum computer\n", "\n", "\n", "sampler = Sampler(mode=backend)\n", "pubs = [qc_isa]\n", "job = sampler.run(pubs)\n", "res = job.result()\n", "\n", "counts = res[0].data.c.get_counts()\n", "\n", "# Run the job on the Aer simulator with noise model from real backend\n", "\n", "# job = noisy_sampler.run([qc_isa])\n", "# res=job.result()\n", "# counts=res[0].data.c.get_counts()" ] }, { "cell_type": "markdown", "id": "01dedf77-2ea7-45f4-9a84-8cd1caa26c14", "metadata": {}, "source": [ "\n", "\n", "### ステップ4：後処理\n", "\n", "これは特にシンプルな後処理のケースで、単純にカウント数を視覚化する。\n", "\n", "Qiskitは、量子ビット、測定値、その他のものを、番号の小さいものを最後/右側に並べるという、「リトルエンディアン」と呼ばれる慣例があることに注意してください。つまり、以下の「10」と書かれた列は、1回目の測定で「0」、2回目の測定で「1」が得られたカウントを指している。\n", "\n" ] }, { "cell_type": "code", "execution_count": 7, "id": "7cfeafbf-6a40-4829-9e1d-6c8dbdab9f23", "metadata": {}, "outputs": [ { "data": { "text/plain": [ " $\"Output$ " ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# Step 4: Post-process\n", "\n", "from qiskit.visualization import plot_histogram\n", "\n", "plot_histogram(counts)" ] }, { "cell_type": "markdown", "id": "2d21379b-b27f-4448-adb5-91f55e79b0e2", "metadata": {}, "source": [ "この慣例に魅力を感じない場合は、 `marginal_counts` 、各測定結果を個別に視覚化することができる：\n", "\n" ] }, { "cell_type": "code", "execution_count": 8, "id": "8c94042b-0e52-41f3-8231-ede9fab3ee51", "metadata": {}, "outputs": [ { "data": { "text/plain": [ " $\"Output$ " ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "from qiskit.result import marginal_counts\n", "\n", "plot_histogram(\n", " marginal_counts(counts, indices=[0]), title=\"Counts after first measurement\"\n", ")" ] }, { "cell_type": "code", "execution_count": 9, "id": "d0ad384b-c2f1-4bc6-8730-ec8a4b593a50", "metadata": {}, "outputs": [ { "data": { "text/plain": [ " $\"Output$ " ] }, "execution_count": 9, "metadata": {}, "output_type": "execute_result" } ], "source": [ "plot_histogram(\n", " marginal_counts(counts, indices=[1]), title=\"Counts after second measurement\"\n", ")" ] }, { "cell_type": "markdown", "id": "258d66a0-ddfc-4c69-8943-37efb7c09b40", "metadata": {}, "source": [ "デフォルトでは、Qiskitの状態は $|0\\rangle$ 状態に初期化されます。だから、最初の測定のほとんどすべてが $|0\\rangle$。しかし、2回目の測定（ $x$ への状態の投影に関する情報を与えるもの）では、ほぼ均等に分かれたことに注意してほしい。この状態は、 $z$ に沿った測定では非常に予測可能な結果を与えるが、 $x$ に沿った測定では非常に予測不可能な結果を与えるようだ。これを探ってみよう。\n", "\n", "逆の順序で測定したらどうなるか？まず、ハダマードゲートを使って、 $|0\\rangle$ が $|\\pm\\rangle_x$ で測定される確率の統計を取ることから始める。次に、2回目の測定では、2回目のハダマードゲートを使って $z$。\n", "\n" ] }, { "cell_type": "code", "execution_count": 10, "id": "4029873f-1260-4fc3-bfd2-e5a3f579a09c", "metadata": {}, "outputs": [ { "data": { "text/plain": [ " $\"Output$ " ] }, "execution_count": 10, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# Step 1:\n", "\n", "# Define registers\n", "qr = QuantumRegister(1, \"q\")\n", "cr = ClassicalRegister(2, \"c\")\n", "qc = QuantumCircuit(qr, cr)\n", "\n", "# Change basis to measure along x.\n", "qc.h(qr)\n", "qc.measure(qr, cr[0])\n", "qc.barrier()\n", "\n", "# Change our basis back to z and make a second measurement\n", "qc.h(qr)\n", "qc.measure(qr, cr[1])\n", "\n", "qc.draw(\"mpl\")" ] }, { "cell_type": "code", "execution_count": 11, "id": "c56d1649-de18-4417-8033-2be2a04a2582", "metadata": {}, "outputs": [], "source": [ "# Step 2: Transpile the circuit for running on a quantum computer\n", "\n", "pm = generate_preset_pass_manager(target=target, optimization_level=3)\n", "qc_isa = pm.run(qc)" ] }, { "cell_type": "code", "execution_count": null, "id": "782f01c9-f9c2-4e5e-8638-7c9ba3347bfb", "metadata": {}, "outputs": [], "source": [ "# Step 3: Run the job on a real quantum computer\n", "\n", "sampler = Sampler(mode=backend)\n", "pubs = [qc_isa]\n", "job = sampler.run(pubs)\n", "res = job.result()\n", "counts = res[0].data.c.get_counts()\n", "\n", "# Run the job on the Aer simulator with noise model from real backend\n", "# job = noisy_sampler.run([qc_isa])\n", "# res=job.result()\n", "# counts=res[0].data.c.get_counts()" ] }, { "cell_type": "code", "execution_count": 13, "id": "ce3fb0fd-a859-4c91-b2e5-c09343a397c2", "metadata": {}, "outputs": [ { "data": { "text/plain": [ " $\"Output$ " ] }, "execution_count": 13, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# Step 4: Post-process\n", "from qiskit.visualization import plot_histogram\n", "\n", "plot_histogram(counts)" ] }, { "cell_type": "markdown", "id": "d7862ea1-a9fd-4b23-9198-cfbfb05cf3dc", "metadata": {}, "source": [ "ここでは、予測可能性はさらに低いようだ！以前は、少なくとも最初の測定でどのような結果が出るかわかっていた。なぜこのようなことが起きたのか、その理由を理解するのはそれほど難しいことではない。 $|0\\rangle$ $|+\\rangle_x$ とを半々で混ぜたものである。したがって、最初の測定で＋または-の状態（チャートでは0と1にマップされている）を得る確率は、明らかに等しいはずである。 $|-\\rangle_x$ $|0\\rangle=\\frac{1}{\\sqrt{2}}(|+\\rangle_x+|-\\rangle_x)$ $x$ に沿った測定は、状態を固有状態 $|+\\rangle_x$ か固有状態 $|-\\rangle_x$ のどちらかに崩壊させる。これらの状態はそれぞれ、 $|0\\rangle$ と $|1\\rangle$ の50対50の混合物である。 $|+\\rangle_x = \\frac{1}{\\sqrt{2}}(|0\\rangle+|1\\rangle)$ $|-\\rangle_x = \\frac{1}{\\sqrt{2}}(|0\\rangle-|1\\rangle)$ したがって、系が $x$ の固有状態にあるとき、 $z$ に沿って測定すると、 $|0\\rangle$ と $|1\\rangle$ の両方が得られ、ほぼ等しい確率で得られることは明らかである。\n", "つまり、最初の例では、州によっては、ある測定結果については非常に予測しやすいが、他の測定結果については予測不可能であることが示された。今回の例は、それよりももっと悪いことができることを示している。測定の順番を入れ替えただけでも、両方の測定で予測不可能な結果が得られる状態がある。ある量について、ある状態がどの程度確かなのか、あるいは不確かなのかを調べてみよう。\n", "\n" ] }, { "cell_type": "markdown", "id": "292b878f-1792-4f55-8991-49e37c7db58b", "metadata": {}, "source": [ "\n", "\n", "## 不確実性の計算\n", "\n", "これを定量化するには、不確実性（分散）を用いる。不確かさ」はしばしば分布の「分散」の平方根と定義される。すなわち、ある観測値 $S$ の不確かさを $\\Delta S$ と表し、次式で与えられる\n", "\n", "$$\n", "\\begin{aligned}\n", "(\\Delta S)^2 & \\equiv \\langle (S - \\langle S \\rangle)^2 \\rangle\\\\\n", "\n", "(\\Delta S)^2 & = \\langle S^2 - 2 S \\langle S \\rangle +\\langle S \\rangle^2 \\rangle\\\\\n", "\n", "(\\Delta S)^2 & = \\langle S^2 \\rangle - \\langle S \\rangle^2\n", "\\end{aligned}\n", "$$\n", "\n", "パウリ行列の場合、 $S^2 = I$、これは次のようになる\n", "\n", "$$\n", "(\\Delta S)^2 = 1 - \\langle S \\rangle^2\n", "$$\n", "\n", "これを具体的な例に当てはめてみよう。状態 $|\\psi\\rangle = |+\\rangle_y = \\frac{1}{\\sqrt{2}}\\begin{pmatrix} 1 \\\\ i \\end{pmatrix},$、その状態における観測可能な $X$ の不確かさを決定することから始めよう。\n", "\n", "\n", "\n", "#### 理解度チェック\n", "\n", "$|+\\rangle_y = |+i\\rangle$ の状態 $X$ の不確かさを手計算する。\n", "\n", "\n", " \n", " $$\n", " \\Delta X =\\sqrt{\\langle+i| X^2 |+i\\rangle - \\langle+i| X |+i\\rangle^2}\n", " $$\n", "\n", " 与えられた状態では、こうなる：\n", "\n", " $$\n", " \\begin{aligned}\n", " \\Delta X & =\\sqrt{\\frac{1}{\\sqrt{2}}\\begin{pmatrix}1 & -i\\end{pmatrix} \\begin{pmatrix}0 & 1 \\\\ 1 & 0\\end{pmatrix} \\begin{pmatrix}0 & 1 \\\\ 1 & 0\\end{pmatrix} \\frac{1}{\\sqrt{2}}\\begin{pmatrix}1 \\\\ i\\end{pmatrix} - \\left(\\frac{1}{\\sqrt{2}}\\begin{pmatrix}1 & -i\\end{pmatrix} \\begin{pmatrix}0 & 1 \\\\ 1 & 0\\end{pmatrix} \\frac{1}{\\sqrt{2}}\\begin{pmatrix}1 \\\\ i\\end{pmatrix}\\right)^2}\\\\\n", "\n", " \\Delta X & =\\sqrt{\\frac{1}{2}\\begin{pmatrix}1 & -i\\end{pmatrix} \\begin{pmatrix}1 & 0 \\\\ 0 & 1\\end{pmatrix} \\begin{pmatrix}1 \\\\ i\\end{pmatrix} - \\left(\\frac{1}{2}\\begin{pmatrix}1 & -i\\end{pmatrix} \\begin{pmatrix}i \\\\ 1\\end{pmatrix}\\right)^2}\\\\\n", "\n", " \\Delta X & =\\sqrt{\\frac{1}{2}\\begin{pmatrix}1 & -i\\end{pmatrix} \\begin{pmatrix}1 \\\\ i\\end{pmatrix} - \\left(\\frac{1}{2}(0))\\right)^2}\\\\\n", "\n", " \\Delta X & =\\sqrt{\\frac{1}{2}(2)} = 1\n", " \\end{aligned}\n", " $$\n", " \n", "\n", "\n", "`qc.initialize()` を使って任意の初期状態を作ることができる。なお、ここでの虚数単位のシンタックスは $1j$ である。\n", "\n" ] }, { "cell_type": "code", "execution_count": null, "id": "bf16b90b-d77a-477b-b432-963f8b81feed", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "-0.02408454165642664\n" ] } ], "source": [ "# Step 1: Map the problem into a quantum circuit\n", "\n", "from qiskit.quantum_info import SparsePauliOp\n", "import numpy as np\n", "\n", "obs = SparsePauliOp(\"X\")\n", "\n", "# Define registers\n", "qr = QuantumRegister(1, \"q\")\n", "cr = ClassicalRegister(1, \"c\")\n", "qc = QuantumCircuit(qr, cr)\n", "\n", "# Initialize the state\n", "qc.initialize([1, 1j] / np.sqrt(2))\n", "\n", "# Step 2: Transpile the circuit\n", "\n", "pm = generate_preset_pass_manager(target=target, optimization_level=3)\n", "qc_isa = pm.run(qc)\n", "obs_isa = obs.apply_layout(layout=qc_isa.layout)\n", "\n", "# Step 3: Run the circuit on a real quantum computer\n", "\n", "\n", "estimator = Estimator(mode=backend)\n", "pubs = [(qc_isa, obs_isa)]\n", "job = estimator.run([[qc_isa, obs_isa]])\n", "res = job.result()\n", "\n", "# Run the job on the Aer simulator with noise model from real backend\n", "# job = noisy_estimator.run([[qc_isa,obs_isa]])\n", "# res=job.result()\n", "\n", "# Step 4: Return the result in classical form, and analyze.\n", "\n", "print(res[0].data.evs)" ] }, { "cell_type": "markdown", "id": "e162abde-f346-4f2c-8358-d95c818a1021", "metadata": {}, "source": [ "上の式によれば $(\\Delta X)^2 = 1 - \\langle X \\rangle^2 = 1-(0.0015...)^2 \\rightarrow \\Delta X = 0.999...$ 同じ状態にこだわって、 $Z$ の期待値を求めてみよう：\n", "\n" ] }, { "cell_type": "code", "execution_count": null, "id": "c1918196-a26e-4419-8a25-488a3d1d4cc1", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "0.04958271968581247\n" ] } ], "source": [ "# Step 1: Map the problem into a quantum circuit\n", "\n", "obs = SparsePauliOp(\"Z\")\n", "\n", "# Define registers\n", "qr = QuantumRegister(1, \"q\")\n", "cr = ClassicalRegister(1, \"c\")\n", "qc = QuantumCircuit(qr, cr)\n", "\n", "# Initialize the state to |+>_y\n", "qc.initialize([1, 1j] / np.sqrt(2))\n", "\n", "\n", "# Step 2: Transpile the circuit\n", "\n", "pm = generate_preset_pass_manager(target=target, optimization_level=3)\n", "qc_isa = pm.run(qc)\n", "obs_isa = obs.apply_layout(layout=qc_isa.layout)\n", "\n", "# Step 3: Run the circuit on a real quantum computer\n", "\n", "estimator = Estimator(mode=backend)\n", "pubs = [(qc_isa, obs_isa)]\n", "job = estimator.run(pubs)\n", "res = job.result()\n", "\n", "# Run the job on the Aer simulator with noise model from real backend\n", "# job = noisy_estimator.run([[qc_isa,obs_isa]])\n", "# res=job.result()\n", "\n", "# Step 4: Return the result in classical form, and analyze.\n", "\n", "print(res[0].data.evs)" ] }, { "cell_type": "markdown", "id": "8ba18eed-5982-45ca-b0ed-5708b7d4aa17", "metadata": {}, "source": [ "先ほどと同じ計算ができるが、分散は再び 1.0 に非常に近いことがわかるだろう。 $\\Delta X \\Delta Z \\approx 1.0$ と結論づけることができる。実際、我々が選んだ状態では、これはほぼ正しい。しかし、我々はもっとうまくやれるだろうか？それとももっと悪い？\n", "\n", "一方向に沿った位置、 $x,$ と同じ方向に沿った運動量、 $p_x.$ の間に不確定性関係があることを思い出してほしい。 $\\Delta x \\Delta p_x \\geq \\hbar/2$ これしか覚えていないと、 $\\Delta X$ と $\\Delta Z$ も不確かさについてこのような基本的な限界を持っているのではないかと考えたくなるかもしれない。もしかしたら、 $\\Delta X \\Delta Z$、ゼロになることは不可能なのだろうか？別の州を試してみて、これが成り立つかどうか見てみよう。今回は $|\\psi\\rangle = \\frac{1}{\\sqrt{2}}\\begin{pmatrix} 1 \\\\ 1 \\end{pmatrix}.$ を使ってみよう。以下のコードでは、estimatorは同じジョブ投入で2組の回路と観測値を受け取ることができる。\n", "\n" ] }, { "cell_type": "code", "execution_count": null, "id": "ef107fd6-03c2-42f2-a51d-e52957cfb82e", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "The expectation value of the first observable is: 1.0011036174126302\n", "The expectation value of the second observable is: 0.0029429797670141016\n" ] } ], "source": [ "# Step 1: Map the problem into a quantum circuit\n", "\n", "obs1 = SparsePauliOp(\"X\")\n", "obs2 = SparsePauliOp(\"Z\")\n", "\n", "# Define registers\n", "\n", "qr = QuantumRegister(1, \"q\")\n", "cr = ClassicalRegister(1, \"c\")\n", "qc = QuantumCircuit(qr, cr)\n", "\n", "# Initialize the state\n", "qc.initialize([1, 1] / np.sqrt(2))\n", "\n", "# Step 2: Transpile the circuit\n", "\n", "pm = generate_preset_pass_manager(target=target, optimization_level=3)\n", "qc_isa = pm.run(qc)\n", "obs1_isa = obs1.apply_layout(layout=qc_isa.layout)\n", "obs2_isa = obs2.apply_layout(layout=qc_isa.layout)\n", "\n", "# Step 3: Run the circuit on a real quantum computer\n", "\n", "with Batch(backend=backend) as batch:\n", " estimator = Estimator(mode=batch)\n", " pubs = [(qc_isa, obs1_isa), (qc_isa, obs2_isa)]\n", " job = estimator.run(pubs)\n", " res = job.result()\n", "batch.close()\n", "\n", "# Run the job on the Aer simulator with noise model from real backend\n", "\n", "# job = noisy_estimator.run([[qc,obs1],[qc,obs2]])\n", "# res=job.result()\n", "\n", "# Step 4: Return the result in classical form, and analyze.\n", "\n", "print(\"The expectation value of the first observable is: \", res[0].data.evs)\n", "print(\"The expectation value of the second observable is: \", res[1].data.evs)" ] }, { "cell_type": "markdown", "id": "489472e8-eefc-4740-aa7e-eba8e71d0d72", "metadata": {}, "source": [ "$X$ の期待値は 1.0 に近いが、 1.0 を超えてはならない。 1.0 をごく少量上回ったとしても気にする必要はない。これはノイズや読み出しエラーなどの要因によるものである。これは非常に重要なテーマだが、今は無視してもいいだろう。\n", "\n", "我々は、 1.0 に非常に近い $X$ の期待値を得た（ $X$ の分散が非常に低いことに対応する）。これにより、2つの分散の積はかなり低くなる：\n", "\n", "$$\n", "\\Delta X \\Delta Z = \\sqrt{1-(0.9853)^2} \\times \\sqrt{1-(-0.00195)^2} = 0.171.\n", "$$\n", "\n", "これは正確にはゼロではないが、この値はパウリ作用素の固有値（ $\\pm 1$ ）に比べて小さくなっている。さて、線形位置と運動量の間の不確定性関係は、演算子 $x$ と $p_x$ の間の転置関係を用いて、別の書き方ができることを思い出してほしい：\n", "\n", "$$\n", "\\Delta x \\Delta p_x \\geq \\frac{1}{2}|\\langle [x,p_x] \\rangle|\n", "$$\n", "\n", "ここで\n", "\n", "$$\n", "[x,p_x] = xp_x-p_xx\n", "$$\n", "\n", "は、 $x$ と $p_x$ のコミュテーターである。\n", "\n", "これはパウリ作用素に最も簡単に拡張できる形式である。一般に、2つの演算子 $A$ と $B$、\n", "\n", "$$\n", "\\Delta A \\Delta B \\geq \\frac{1}{2}|\\langle [A,B] \\rangle|.\n", "$$\n", "\n", "また、パウリ行列の場合、 $X$ と $Z$ を計算するためには、 $[X,Z]$ が必要である\n", "\n", "$$\n", "\\Delta X \\Delta Z \\geq \\frac{1}{2}|\\langle [X,Z] \\rangle|.\n", "$$\n", "\n", "ここではこれを示し、同様の計算は練習として読者に委ねる：\n", "\n", "$$\n", "[X,Z] = XZ-ZX = \\begin{pmatrix}0 & 1 \\\\ 1 & 0\\end{pmatrix}\\begin{pmatrix}1 & 0 \\\\ 0 & -1\\end{pmatrix}-\\begin{pmatrix}1 & 0 \\\\ 0 & -1\\end{pmatrix}\\begin{pmatrix}0 & 1 \\\\ 1 & 0\\end{pmatrix}\n", "$$\n", "\n", "$$\n", "[X,Z] = \\begin{pmatrix}0 & -1 \\\\ 1 & 0\\end{pmatrix}-\\begin{pmatrix}0 & 1 \\\\ -1 & 0\\end{pmatrix} = 2\\begin{pmatrix}0 & -1 \\\\ 1 & 0\\end{pmatrix}\n", "$$\n", "\n", "これは完全に受け入れられる答えだが、もう一歩踏み込むとこうなる\n", "\n", "$$\n", "[X,Z] = -2i\\begin{pmatrix}0 & -i \\\\ i & 0\\end{pmatrix}=-2iY\n", "$$\n", "\n", "したがって、不確定性関係は次のようになる\n", "\n", "$$\n", "\\Delta X \\Delta Z \\geq |\\langle Y \\rangle|.\n", "$$\n", "\n", "\n", "\n", "#### 理解度チェック\n", "\n", "$[X,Y]$ と $[Y,Z]$ を決定する。これを用いて、 $X$ & $Y$、 $Y$ & $Z$ の間の不確定性関係を書き出す。\n", "\n", "\n", " \n", " $$\n", " [X,Y] = XY-YX = \\begin{pmatrix}0 & 1 \\\\ 1 & 0\\end{pmatrix}\\begin{pmatrix}0 & -i \\\\ i & 0\\end{pmatrix}-\\begin{pmatrix}0 & -i \\\\ i & 0\\end{pmatrix}\\begin{pmatrix}0 & 1 \\\\ 1 & 0\\end{pmatrix} = 2 \\begin{pmatrix}i & 0 \\\\ 0 & -i\\end{pmatrix}=2iZ\n", " $$\n", "\n", " $$\n", " [Y,Z] = YZ-ZY = \\begin{pmatrix}0 & -i \\\\ i & 0\\end{pmatrix}\\begin{pmatrix}1 & 0 \\\\ 0 & -1\\end{pmatrix}-\\begin{pmatrix}1 & 0 \\\\ 0 & -1\\end{pmatrix}\\begin{pmatrix}0 & -i \\\\ i & 0\\end{pmatrix} = 2 \\begin{pmatrix}0 & i \\\\ i & 0\\end{pmatrix}=2iX\n", " $$\n", "\n", " 一般的な不確定性関係と組み合わせると、次のようになる\n", "\n", " $$\n", " \\Delta X \\Delta Y \\geq |\\langle Z \\rangle|,\n", " $$\n", "\n", " $$\n", " \\Delta Y \\Delta Z \\geq |\\langle X \\rangle|.\n", " $$\n", " \n", "\n", "\n", "\n", "\n", "### 一貫性検査\n", "\n", "次に進む前に、これが前回の発見と一致していることを確認しておこう。 $|\\psi\\rangle = \\frac{1}{\\sqrt{2}}\\begin{pmatrix} 1 \\\\ 1 \\end{pmatrix}.$ そして、 $\\Delta X \\Delta Z = 0.171.$、この積が以下の値以上であることがわかった\n", "\n", "$$\n", "|\\langle Y \\rangle|=\\frac{1}{\\sqrt{2}}\\begin{pmatrix} 1 & 1 \\end{pmatrix}\\begin{pmatrix}0 & -i \\\\ i & 0\\end{pmatrix}\\frac{1}{\\sqrt{2}}\\begin{pmatrix} 1 \\\\ 1 \\end{pmatrix}\n", "$$\n", "\n", "$$\n", "|\\langle Y \\rangle| = \\frac{1}{2}\\begin{pmatrix} 1 & 1 \\end{pmatrix}\\begin{pmatrix} -i \\\\ i \\end{pmatrix} = \\frac{1}{2}(-i+i) = 0.\n", "$$\n", "\n", "確かにそうだ、 $\\Delta X \\Delta Z = 0.171 \\geq |\\langle Y \\rangle|=0$\n", "\n" ] }, { "cell_type": "markdown", "id": "0ea50490-36e9-4ab6-b9e1-6840f6dedfec", "metadata": {}, "source": [ "以下の質問を使って、これらの調査結果に対する直感を養ってください：\n", "\n", "\n", "\n", "#### 理解度チェック\n", "\n", "以下の項目をセットで答えなさい：\n", "\n", "(a) $X$、不確かさがゼロになると予想される状態は？\n", "\n", "(b) $Z$、不確かさがゼロになると予想される状態は？\n", "\n", "(c) どのような状態でゼロ期待値 $\\langle Y \\rangle$ が得られるか？\n", "\n", "(d) 上記の質問に対する回答は、ケース（ $\\Delta X \\Delta Z \\geq |\\langle Y \\rangle|$ ）と一致しているか？\n", "\n", "(e) 推定子を用いてこれを明示的にチェックするコードを書く。\n", "\n", "\n", " \n", " (a)我々は、 $X$ オペレータの固有状態が、 $X$ の不確かさをゼロにすることを期待するかもしれない。実際、 $|\\psi\\rangle = |+\\rangle_x,$ を使うと、次のようになる。 $\\Delta X = \\sqrt{1-\\langle X \\rangle^2} = \\sqrt{1-1^2} = 0.$\n", "\n", " (b)私たちは、 $Z$ オペレータの固有状態が、 $Z$ の不確かさをゼロにすることを期待するかもしれない。実際、 $|\\psi\\rangle = |1\\rangle,$ を使うと、次のようになる。 $\\Delta Z = \\sqrt{1-\\langle Z \\rangle^2} = \\sqrt{1-(-1)^2} = 0.$\n", "\n", " (c)どのような状態であれ、測定したときに、 $y$、負の投影と同じくらい頻繁に、 $\\langle Y \\rangle=0$、正の投影が得られると予想される。これには $X$ と $Z$ の固有状態が含まれる。\n", "\n", " (d) はい。 $X$、 $Z$ の固有状態では、不確かさの積 $\\Delta X \\Delta Z $ の値が非常に小さくなると予想される： $\\Delta X \\Delta Z \\approx 0.$ これは、同じ状態に対して $\\langle Y \\rangle=0$。つまり、不確定性関係を満たすことができる。\n", "\n", " (e)以下のようなコードで確認できる：\n", "\n", " ```python\n", " obs1 = SparsePauliOp.from_list(\n", " [(\"X\", 1.000)]\n", " )\n", " obs2 = SparsePauliOp.from_list(\n", " [(\"Y\", 1.000)]\n", " )\n", " obs3 = SparsePauliOp.from_list(\n", " [(\"Z\", 1.000)]\n", " )\n", " qc = QuantumCircuit(1,1)\n", " qc.ry(pi/2,0)\n", "\n", " job = estimator.run([(qc, [[obs1], [obs2], [obs3]])], precision=0.001)\n", " res=job.result()\n", " ```\n", "\n", " 結果はすべての期待値を返す。すべての期待値を取得し、不確実性を計算するには、次のようにする：\n", "\n", " ```python\n", " xs=res[0].data.evs[0]\n", " ys=abs(res[0].data.evs[1])\n", " zs=res[0].data.evs[2]\n", "\n", " import math\n", " prodxz=((1-xs[i]*xs[i])**0.5)*(1-zs[i]*zs[i])**0.5\n", "\n", " ```\n", " \n", "\n", "\n", "以下の項目をセットで答えなさい：\n", "\n", "(a)大きな期待値 $\\langle Y \\rangle$ を持つような状態が考えられるか？\n", "\n", "(b)同じ状態であれば、 $X$ の不確かさは大きいか小さいか？\n", "\n", "(c)同じ状態であれば、 $Z$ の不確かさは大きいか小さいか？\n", "\n", "(d) 上記の質問に対する回答は、ケース（ $\\Delta X \\Delta Z \\geq |\\langle Y \\rangle|$ ）と一致しているか？\n", "\n", "(e) 推定子を用いてこれを明示的にチェックするコードを書く。\n", "\n", "\n", " \n", " (a) $Y$ : $|+\\rangle_y$ の固有状態について、 $\\langle Y \\rangle\\approx 1$ が見つかると予想される。\n", "\n", " (b) $|+\\rangle_y,$ の状態で $X$ を測定すると、同じ頻度／確率でプラスとマイナスの結果が得られるため、 $X$ の状態には大きな不確かさがあると予想される。\n", "\n", " (c) $|+\\rangle_y,$ の状態で $Z$ を測定すると、同じ頻度／確率でプラスとマイナスの結果が得られるため、 $Z$ の状態には大きな不確かさがあると予想される。\n", "\n", " (d) はい。 $Y,$、特に $|+\\rangle_y,$ の固有状態については、不確定性の積 $\\Delta X \\Delta Z $ の値が大きいと予想される。また、 $\\langle Y \\rangle\\approx 1$。つまり、この状態では $\\langle Y \\rangle$、 $\\Delta X \\Delta Z $、どちらもかなり大きく、不確定性関係が再び満たされる可能性は十分にある。\n", "\n", " (e)以下のようなコードで確認できる：\n", "\n", " ```python\n", " obs1 = SparsePauliOp.from_list(\n", " [(\"X\", 1.000)]\n", " )\n", " obs2 = SparsePauliOp.from_list(\n", " [(\"Y\", 1.000)]\n", " )\n", " obs3 = SparsePauliOp.from_list(\n", " [(\"Z\", 1.000)]\n", " )\n", " qc = QuantumCircuit(1,1)\n", " qc.rx(-pi/2,0)\n", "\n", " job = estimator.run([(qc, [[obs1], [obs2], [obs3]])], precision=0.001)\n", " res=job.result()\n", " ```\n", "\n", " 結果はすべての期待値を返す。すべての期待値を取得し、不確実性を計算するには、次のようにする：\n", "\n", " ```python\n", " xs=res[0].data.evs[0]\n", " ys=abs(res[0].data.evs[1])\n", " zs=res[0].data.evs[2]\n", "\n", " import math\n", " prodxz=((1-xs[i]*xs[i])**0.5)*(1-zs[i]*zs[i])**0.5\n", "\n", " ```\n", " \n", "\n", "\n" ] }, { "cell_type": "markdown", "id": "26f11700-d6be-4ce1-9f97-b2d47ea933ec", "metadata": {}, "source": [ "\n", "\n", "## 不確定性関係の検証\n", "\n", "上記のテストは、状態ベクトル $|\\psi\\rangle = |+\\rangle_x$ の単一選択に対する不確定性関係の妥当性を示したに過ぎない。これが一般的に実験と一致していることを確信するためには、状態ベクトルの多くの選択肢について推定量を用いて同様の計算を行う必要がある。 RY\\`ゲートを使って、パラメータ $\\theta$ を使って異なる初期状態を生成するために、状態ベクトルを $z$ の軸から回転させることから始めよう。\n", "\n" ] }, { "cell_type": "code", "execution_count": null, "id": "c434a14d-c0ce-40bd-8f5e-7d35722e97bb", "metadata": {}, "outputs": [], "source": [ "# The calculation below uses approximately 3-4 minutes of QPU time.\n", "# Step 1: Map the problem into a quantum circuit\n", "\n", "from qiskit.circuit import Parameter\n", "import numpy as np\n", "\n", "# Specify observables\n", "obs1 = SparsePauliOp(\"X\")\n", "obs2 = SparsePauliOp(\"Y\")\n", "obs3 = SparsePauliOp(\"Z\")\n", "\n", "# Define registers\n", "qr = QuantumRegister(1, \"q\")\n", "cr = ClassicalRegister(1, \"c\")\n", "qc = QuantumCircuit(qr, cr)\n", "\n", "# Rotate away from |0>\n", "theta = Parameter(\"θ\")\n", "qc.ry(theta, 0)\n", "\n", "params = np.linspace(0, 2, num=21)\n", "\n", "# Step 2: Transpile the circuit\n", "\n", "pm = generate_preset_pass_manager(target=target, optimization_level=3)\n", "qc_isa = pm.run(qc)\n", "obs1_isa = obs1.apply_layout(layout=qc_isa.layout)\n", "obs2_isa = obs2.apply_layout(layout=qc_isa.layout)\n", "obs3_isa = obs3.apply_layout(layout=qc_isa.layout)\n", "\n", "# Step 3: Run the circuit on a real quantum computer\n", "\n", "with Batch(backend=backend) as batch:\n", " estimator = Estimator(mode=batch)\n", " pubs = [(qc_isa, [[obs1_isa], [obs2_isa], [obs3_isa]], [params])]\n", " job = estimator.run(pubs, precision=0.01)\n", " res = job.result()\n", "\n", "batch.close()\n", "\n", "# Run the job on the Aer simulator with noise model from real backend\n", "\n", "# job = noisy_estimator.run([(qc, [[obs1], [obs2], [obs3]], [params])])\n", "# res=job.result()" ] }, { "cell_type": "code", "execution_count": 18, "id": "41d974bb-da9e-459e-8145-294f71cab4c7", "metadata": {}, "outputs": [], "source": [ "# Step 4: Post-processing and classical analysis.\n", "xs = res[0].data.evs[0]\n", "ys = abs(res[0].data.evs[1])\n", "zs = res[0].data.evs[2]\n", "\n", "# Calculate uncertainties\n", "\n", "delx = []\n", "delz = []\n", "prodxz = []\n", "for i in range(len(xs)):\n", " delx.append(abs((1 - xs[i] * xs[i])) ** 0.5)\n", " delz.append(abs((1 - zs[i] * zs[i])) ** 0.5)\n", " prodxz.append(delx[i] * delz[i])" ] }, { "cell_type": "code", "execution_count": 19, "id": "1599d4df-38f8-4171-95ae-6c3c61409201", "metadata": {}, "outputs": [ { "data": { "text/plain": [ " $\"Output$ " ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# Here we can plot the results from this simulation.\n", "import matplotlib.pyplot as plt\n", "\n", "plt.plot(params, delx, label=r\"$\\Delta$ X\")\n", "plt.plot(params, ys, label=r\"$\\langle$ Y $\\rangle$\")\n", "plt.plot(params, delz, label=r\"$\\Delta$ Z\")\n", "plt.plot(params, prodxz, label=r\"$\\Delta$X $\\Delta$Z\")\n", "plt.xlabel(r\"$\\theta$\")\n", "plt.ylabel(\"Expectation/Uncertainty Values\")\n", "plt.legend()\n", "plt.show()" ] }, { "cell_type": "markdown", "id": "1694ab4f-b5e2-46be-b453-2599f73c582d", "metadata": {}, "source": [ "赤い曲線 $(\\Delta X \\Delta Z)$ は、常にオレンジの曲線 $\\langle Y \\rangle.$ よりも大きいことに注意。不確実性積が沈んで限界にいくらか近づくこともあれば、上昇して限界から遠ざかることもあるが、常に不確定性関係に従う。\n", "\n", "もちろん、これは不確定性関係の最良のテストではないかもしれない。私たちの極限 $\\langle Y \\rangle$ は常にゼロに非常に近いからだ。 $Y$。具体的には、 $|0\\rangle$ を $z$ の軸から様々な角度で回転させ、の固有状態により大きな投影を持つ量子状態を使ってみよう。しかし今度は、その結果として得られる状態を、 $z$ の周りにある角度、おそらく $\\pi/4$、回転させて、何が起こるか見てみよう。\n", "\n" ] }, { "cell_type": "code", "execution_count": null, "id": "c9700dd1-d116-4f94-8c81-70383ce38671", "metadata": {}, "outputs": [], "source": [ "# The calculation below uses approximately 3-4 minutes of QPU time.\n", "from qiskit.circuit import Parameter\n", "import numpy as np\n", "\n", "# Step 1: Map the problem to a quantum circuit\n", "\n", "# Specify observables\n", "obs1 = SparsePauliOp(\"X\")\n", "obs2 = SparsePauliOp(\"Y\")\n", "obs3 = SparsePauliOp(\"Z\")\n", "\n", "# Define registers\n", "qr = QuantumRegister(1, \"q\")\n", "cr = ClassicalRegister(1, \"c\")\n", "qc = QuantumCircuit(qr, cr)\n", "\n", "# Rotate away from |0> along one plane, and then along a transverse direction.\n", "theta = Parameter(\"θ\")\n", "qc.ry(theta, 0)\n", "qc.rz(pi / 4, 0)\n", "\n", "params = np.linspace(0, 2, num=21)\n", "\n", "# Step 2: Transpile the circuit\n", "\n", "pm = generate_preset_pass_manager(target=target, optimization_level=3)\n", "qc_isa = pm.run(qc)\n", "\n", "obs1_isa = obs1.apply_layout(layout=qc_isa.layout)\n", "obs2_isa = obs2.apply_layout(layout=qc_isa.layout)\n", "obs3_isa = obs3.apply_layout(layout=qc_isa.layout)\n", "\n", "# Step 3: Run the circuit on a real quantum computer\n", "\n", "with Batch(backend=backend) as batch:\n", " estimator = Estimator(mode=batch)\n", " pubs = [(qc_isa, [[obs1_isa], [obs2_isa], [obs3_isa]], [params])]\n", " job = estimator.run(pubs, precision=0.01)\n", " res = job.result()\n", "\n", "batch.close()\n", "\n", "# Run the job on the Aer simulator with noise model from real backend\n", "\n", "# job = noisy_estimator.run([(qc, [[obs1], [obs2], [obs3]], [params])])\n", "# res=job.result()" ] }, { "cell_type": "code", "execution_count": 21, "id": "7e85e759-33a7-4539-80f7-05c769c46144", "metadata": {}, "outputs": [], "source": [ "# Step 4: Post-processing and classical analysis.\n", "xs = res[0].data.evs[0]\n", "ys = abs(res[0].data.evs[1])\n", "zs = res[0].data.evs[2]\n", "\n", "# Calculate uncertainties\n", "\n", "delx = []\n", "delz = []\n", "prodxz = []\n", "for i in range(len(xs)):\n", " delx.append(abs((1 - xs[i] * xs[i])) ** 0.5)\n", " delz.append(abs((1 - zs[i] * zs[i])) ** 0.5)\n", " prodxz.append(delx[i] * delz[i])" ] }, { "cell_type": "code", "execution_count": 22, "id": "6f129a2d-e029-45b7-bb77-bfa83babd8c9", "metadata": {}, "outputs": [ { "data": { "text/plain": [ " $\"Output$ " ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# Here we can plot the results from this simulation.\n", "import matplotlib.pyplot as plt\n", "\n", "plt.plot(params, delx, label=r\"$\\Delta$ X\")\n", "plt.plot(params, ys, label=r\"$\\langle$ Y $\\rangle$\")\n", "plt.plot(params, delz, label=r\"$\\Delta$ Z\")\n", "plt.plot(params, prodxz, label=r\"$\\Delta$X $\\Delta$Z\")\n", "plt.xlabel(r\"$\\theta$\")\n", "plt.ylabel(\"Expectation/Uncertainty Values\")\n", "plt.legend()\n", "plt.show()" ] }, { "cell_type": "markdown", "id": "7b90bcd5-f561-4f7f-a5a3-7882b3ddb464", "metadata": {}, "source": [ "今や不確定性の限界 $(\\Delta X \\Delta Z)$ が試されていることが明らかになった！赤い曲線は以前よりもオレンジ色の曲線にずっと近づいている。実際、ノイズが存在しない場合、不確定性関係は一点で完全に飽和する（ $(\\Delta X \\Delta Z) = \\langle Y \\rangle$ ）。ノイズや読み出し誤差が存在する場合、測定結果が時折わずかに大きい値 $(\\Delta X \\Delta Z)$ を示すことは驚くべきことではない。 $\\langle Y \\rangle.$ これは不確定性の真の違反ではなく、単に誤差がゼロでないことに起因する現象である。\n", "\n", "\n", "\n", "#### 理解度チェック\n", "\n", "$\\langle Y \\rangle$ を可能な限り大きくし、これを限界まで追求する方法を説明してください。\n", "\n", "\n", " \n", " 現在のコードには、デフォルトの初期状態 $|0\\rangle$ を、 $z$ の軸から、パラメトリックな角度 $\\theta$ で回転させ、さらに、 $z$ の軸を中心に、状態ベクトルを $y$ の軸まで回転させる角度 $\\pi/4,$ で回転させる行がある。\n", "\n", " `qc.ry(theta,0)`\n", "\n", " `qc.rz(pi/4,0)`\n", "\n", " $z$ の回転を $\\pi/4$ から $\\pi/2$ に変え、 $Y$ の固有状態まで回転させることができる：\n", "\n", " `qc.ry(theta,0)`\n", "\n", " `qc.rz(pi/2,0)`\n", "\n", " それ以外の変更は必要ない。\n", " \n", "\n", "\n", "コードを変更するかコピーして、Yの期待値を最大化した不確かさ関係のチェックを実装する。不確定性関係は成り立つのか？\n", "\n", "\n", " \n", " 上の例とまったく同じコードを、次のように使います\n", "\n", " `qc.rz(pi/2,0)`\n", "\n", " を交換した\n", "\n", " `qc.rz(pi/4,0).`\n", "\n", " その結果、下の図のようになるはずで、そう、不確定性原理はまだ有効であるはずだ。\n", "\n", " ![不確かさとY演算子の最大期待値を比較したプロット。](https://quantum.cloud.ibm.com/learning/images/modules/quantum-mechanics/exploring-uncertainty-with-qiskit/checkin.avif)\n", " \n", "\n", "\n" ] }, { "cell_type": "markdown", "id": "6dd8c749-00d6-4a08-8d3c-2f5b5f3a0902", "metadata": {}, "source": [ "上のコードを修正して同様の絵を作り、量子コンピューターでの測定から、積 $\\Delta X \\Delta Y$ があるべき振る舞いをすることを示す。好きな州を選んでください。\n", "\n", "\n", " \n", " 上の例とまったく同じコードを使うことになるが、実際には、異なる不確かさを計算するために期待値を使うだけで、上と同じ結果を使うことができる。例えば\n", "\n", " ```python\n", " xs=res[0].data.evs[0]\n", " ys=res[0].data.evs[1]\n", " zs=abs(res[0].data.evs[2])\n", " import math\n", " delx = []\n", " dely = []\n", " prodxy=[]\n", " for i in range(len(xs)):\n", " delx.append((1-xs[i]*xs[i])**0.5)\n", " dely.append((1-ys[i]*ys[i])**0.5)\n", " prodxy.append(((1-xs[i]*xs[i])**0.5)*(1-ys[i]*ys[i])**0.5)\n", " ```\n", "\n", " をプロットすることができる\n", "\n", " ```python\n", " import matplotlib.pyplot as plt\n", " plt.plot(params, delx, label=r'$\\Delta$ X')\n", " plt.plot(params, dely, label=r'$\\langle$ Y $\\rangle$')\n", " plt.plot(params, zs, label=r'$\\Delta$ Z')\n", " plt.plot(params, prodxy, label=r'$\\Delta$X $\\Delta$Z')\n", " plt.xlabel(r'$\\theta$')\n", " plt.ylabel('Expectation/Uncertainty Values')\n", " plt.legend()\n", " plt.show()\n", " ```\n", " \n", "\n", "\n" ] }, { "cell_type": "markdown", "id": "571a0120-a9d3-44c5-9dce-1b021090f9d1", "metadata": {}, "source": [ "課題： $\\theta$ の多くの値をスキャンしたのと同じように、 $\\phi$ の多くの値をスキャンし、不確定性関係が決して違反しないことを示す 3 次元プロットを作成するコードを書く。好きな観測点を選んでください。\n", "\n" ] }, { "cell_type": "markdown", "id": "b89eaae3-6ea9-486c-b93c-55431839d5f1", "metadata": {}, "source": [ "\n", "\n", "## 質問\n", "\n", "指導者は、このノートがどのように使用されているかについての[簡単なアンケートに](https://ibm.biz/classrooms_instructor_key_request)答えることで、解答と一般的なカリキュラムにおける配置についてのガイダンスが付いたバージョンのノートを要求することができる。\n", "\n", "\n", "\n", "### 重要な概念：\n", "\n", "* 位置や直線運動量、スピンの成分など、多くの物理的観測値の間には不確定性関係がある。\n", "* パウリ行列は交わらない。これは、スピンのすべての成分を同時に知る／決定することはできないという事実を数学的に反映したものである。\n", "* 量子コンピュータはパウリ作用素/行列を多用するので、パウリ作用素の不確定性関係を知っておくことは、密接に関連するスピン作用素と同様に有用である。\n", "* 2つの演算子 $A$ と $B$ の不確かさを表す一般式は以下の通りである。 $\\Delta A \\Delta B \\geq \\frac{1}{2}|\\langle [A,B] \\rangle|.$\n", "* ある演算子 $A$ の固有状態 $|a\\rangle$ は、その演算子に関連する物理的観測量の不確定性をゼロにする。実験的に見ても、 $\\Delta A \\approx 0.$\n", "* ある演算子 $A$ の固有状態 $|a\\rangle$ は、 $A$ と交わらない演算子 $B$ に対して、より大きな不確かさをもたらす。\n", "* 実際の量子コンピュータを用いた実験結果は、物理作用素の行列表現から得られる直感を裏付けている。\n", "\n", "\n", "\n", "### 正誤問題：\n", "\n", "1. T/F $X$ と $Y$ を同時に測定することはできるが、 $Z$ を測定することはできない。\n", "2. T/F $X$ と $Z$ を同時に測定することはできるが、 $Y$ を測定することはできない。\n", "3. T/F 線形の位置演算子と線形の運動量演算子は交わらない。\n", "4. T/F IBM 量子コンピュータはデフォルトで $Z$ に沿って測定するので、他の方向に沿って測定するには回転を行わなければならない。\n", "5. T/F 以下の回路は $Z$ を効果的に測定し、次に $X$ を測定する。\n", "\n", "![測定、ハダマードゲート、そして別の測定を示す回路図。](https://quantum.cloud.ibm.com/learning/images/modules/quantum-mechanics/exploring-uncertainty-with-qiskit/tf-q5.avif)\n", "\n", "\n", "\n", "### MC質問：\n", "\n", "1. 下の図は、次の不確かさの関係のどれを示しているか？\n", "\n", " * a. $\\Delta X \\Delta Y \\geq |\\langle Z \\rangle|$\n", " * b. $\\Delta Y \\Delta Z \\geq |\\langle X \\rangle|$\n", " * c. $\\Delta Z \\Delta X \\geq |\\langle Y \\rangle|$\n", " * d. 上記のいずれでもない\n", "\n", "![不確かさとY演算子の最大期待値を比較したプロット。](https://quantum.cloud.ibm.com/learning/images/modules/quantum-mechanics/exploring-uncertainty-with-qiskit/checkin.avif)\n", "\n", "2. $x$ に沿って測定を行う標準的なシーケンスは次のうちどれですか？\n", " * a. のみ `qc.measure()`\n", " * b. `qc.h()` その後 `qc.measure()`\n", " * c. `qc.h()`, `qc.h()` その後 `qc.measure()`\n", " * d. `qc.h()` `qc.s` `qc.h()` `qc.measure()`\n", " * e. `qc.sdg()` `qc.h()`, `qc.s` `qc.measure()`\n", " * f. `qc.sdg()` `qc.h()`, `qc.s`, `qc.h()` それから `qc.measure()`\n", "\n", "3. 次のうち、最大の期待値（ $\\langle X \\rangle$ ）をもたらす状態はどれか？\n", " * a. $|+\\rangle_x$\n", " * b. $|-\\rangle_x$\n", " * c. $|+\\rangle_y$ とも呼ばれる。 $|+i\\rangle$\n", " * d. $|-\\rangle_y$ とも呼ばれる $|-i\\rangle$\n", " * e. $|0\\rangle$ とも呼ばれる $|\\uparrow\\rangle$\n", " * f. $|1\\rangle$ とも呼ばれる。 $|\\downarrow\\rangle$\n", "\n", "4. 次のうち、最大の不確定性（ $\\Delta X$ ）をもたらす状態はどれか？\n", " * a. $|+\\rangle_x$\n", " * b. $|+\\rangle_y$ とも呼ばれる $|+i\\rangle$\n", " * c. $|0\\rangle$ とも呼ばれる。 $|\\uparrow\\rangle$\n", " * d. a と b が同点\n", " * e. bとcは同点\n", " * f. a、b、cが同点\n", "\n", "\n", "\n", "### ディスカッションの質問：\n", "\n", "1. この不確定性の概念は、直交空間におけるベクトルの矢印としてのスピンの概念と相反するものですか？ブロッホ球ではどうですか？\n", "\n", "2. $x$ 軸と $y$ 軸の中間の方向に沿って測定装置を方向付けるとする。何が起こりますか。この方向に沿って測定できるか？ $X$、 $Y$ における不確実性とどのような関係があるのだろうか？\n", "\n", "3. ここで得られた結果に納得するために、さらにどのような実験をしたいですか？\n", "\n" ] }, { "cell_type": "markdown", "metadata": {}, "id": "a1b8767d", "source": "© IBM Corp., 2017-2026" } ], "metadata": { "in_page_toc_max_heading_level": 2, "in_page_toc_min_heading_level": 2, "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" }, "widgets": { "application/vnd.jupyter.widget-state+json": { "state": {}, "version_major": 2, "version_minor": 0 } } }, "nbformat": 4, "nbformat_minor": 4 }