こんにちは、最新の量子コンピュータを使ってみます。Amazon BraketのAria1プロセッサを使いました。最新のイオントラップ25量子ビットです。ツールを読み込みます。 <pre class="ql-syntax" spellcheck="false">import numpy as np
import matplotlib.pyplot as plt

%matplotlib inline
import string
import time

from braket.circuits import Circuit, Gate, Observable, FreeParameter
from braket.devices import LocalSimulator
from braket.aws import AwsDevice, AwsQuantumTask
</pre> まずは簡単な量子もつれ回路を。<pre class="ql-syntax" spellcheck="false">bell = Circuit().h(0).cnot(0, 1)
print(bell)
</pre> 回路を描画します。<pre class="ql-syntax" spellcheck="false">T : |0|1|
 
q0 : -H-C-
 | 
q1 : ---X-

T : |0|1|
</pre> シミュレーターを使ってまずは結果を見ます。<pre class="ql-syntax" spellcheck="false">device = LocalSimulator()
result = device.run(bell, shots=1000).result()
counts = result.measurement_counts
print(counts)

Counter({&#39;11&#39;: 505, &#39;00&#39;: 495})
</pre> プロットしてみます。<pre class="ql-syntax" spellcheck="false">plt.bar(counts.keys(), counts.values())
plt.xlabel(&#39;bitstrings&#39;)
plt.ylabel(&#39;counts&#39;)
</pre> <img src="https://assets.blueqat.com/public/uploads/us-east-2:4805ff4b-c3cc-4344-b165-86544c34d0bf/2023/10/13/result1.png"/> 次にシミュレーターではなく実機でAria1を指定します。<pre class="ql-syntax" spellcheck="false">aria = AwsDevice(&#34;arn:aws:braket:us-east-1::device/qpu/ionq/Aria-1&#34;)
aria.queue_depth().quantum_tasks
</pre> 上記の回路を投げます。<pre class="ql-syntax" spellcheck="false">aria_task = aria.run(bell, shots=1000)
aria_task_id = aria_task.id
aria_status = aria_task.state()
print(&#39;Status of quantum task:&#39;, aria_status)

Status of quantum task: CREATED
</pre> 以前と違って待ち番号が見れます。<pre class="ql-syntax" spellcheck="false">aria_task.queue_position().queue_position
</pre> 終わると結果が見れます。 <pre class="ql-syntax" spellcheck="false">task_load_aria = AwsQuantumTask(arn=aria_task_id)

status_aria = task_load_aria.state()
print(&#39;Status of (reconstructed) quantum task:&#39;, status_aria)

if status_aria == &#39;COMPLETED&#39;:
  results = task_load_aria.result()
   
  metadata = task_load_aria.metadata()
  shots = metadata[&#39;shots&#39;]
  machine = metadata[&#39;deviceArn&#39;]
  print(&#34;{} shots taken on machine {}.&#34;.format(shots, machine))
   
  counts = results.measurement_counts
  print(&#39;Measurement counts:&#39;, counts)

  plt.bar(counts.keys(), counts.values())
  plt.xlabel(&#39;bitstrings&#39;)
  plt.ylabel(&#39;counts&#39;)
  plt.tight_layout()
  plt.savefig(&#39;bell_ionq.png&#39;, dpi=700)
   
elif status in [&#39;FAILED&#39;, &#39;CANCELLED&#39;]:
  print(&#39;Your quantum task is in terminal status, but has not completed.&#39;)

else:
  print(&#39;Sorry, your quantum task is still being processed and has not been finalized yet.&#39;)
</pre> しばらく待つと結果が出ました。だいたいシミュレータと同じです。エラーが少ないですね。<pre class="ql-syntax" spellcheck="false">Status of (reconstructed) quantum task: COMPLETED
1000 shots taken on machine arn:aws:braket:us-east-1::device/qpu/ionq/Aria-1.
Measurement counts: Counter({&#39;11&#39;: 496, &#39;00&#39;: 478, &#39;01&#39;: 17, &#39;10&#39;: 9})
</pre> <img src="https://assets.blueqat.com/public/uploads/us-east-2:4805ff4b-c3cc-4344-b165-86544c34d0bf/2023/10/13/ionq-ipynb-JupyterLab.png"/> 次にもちょい大きな回路でやってみます。 <pre class="ql-syntax" spellcheck="false">ghz = Circuit().h(0)

for i in range(24):
  ghz.cnot(i, i+1)
   
print(ghz)
</pre> <pre class="ql-syntax" spellcheck="false">T : |0|1|2|3|4|5|6|7|8|9|10|11|12|13|14|15|16|17|18|19|20|21|22|23|24|
 
q0 : -H-C--------------------------------------------------------------
 | 
q1 : ---X-C------------------------------------------------------------
 | 
q2 : -----X-C----------------------------------------------------------
 | 
q3 : -------X-C--------------------------------------------------------
 | 
q4 : ---------X-C------------------------------------------------------
 | 
q5 : -----------X-C----------------------------------------------------
 | 
q6 : -------------X-C--------------------------------------------------
 | 
q7 : ---------------X-C------------------------------------------------
 | 
q8 : -----------------X-C----------------------------------------------
 | 
q9 : -------------------X-C--------------------------------------------
 | 
q10 : ---------------------X--C-----------------------------------------
 | 
q11 : ------------------------X--C--------------------------------------
 | 
q12 : ---------------------------X--C-----------------------------------
 | 
q13 : ------------------------------X--C--------------------------------
 | 
q14 : ---------------------------------X--C-----------------------------
 | 
q15 : ------------------------------------X--C--------------------------
 | 
q16 : ---------------------------------------X--C-----------------------
 | 
q17 : ------------------------------------------X--C--------------------
 | 
q18 : ---------------------------------------------X--C-----------------
 | 
q19 : ------------------------------------------------X--C--------------
 | 
q20 : ---------------------------------------------------X--C-----------
 | 
q21 : ------------------------------------------------------X--C--------
 | 
q22 : ---------------------------------------------------------X--C-----
 | 
q23 : ------------------------------------------------------------X--C--
 | 
q24 : ---------------------------------------------------------------X--

T : |0|1|2|3|4|5|6|7|8|9|10|11|12|13|14|15|16|17|18|19|20|21|22|23|24|
</pre> まずはローカルシミュレータで。<pre class="ql-syntax" spellcheck="false">result = device.run(ghz, shots=1000).result()
counts = result.measurement_counts
print(counts)

Counter({&#39;0000000000000000000000000&#39;: 507, &#39;1111111111111111111111111&#39;: 493})
</pre> <img src="https://assets.blueqat.com/public/uploads/us-east-2:4805ff4b-c3cc-4344-b165-86544c34d0bf/2023/10/13/ghz.png"/> 早速実機で<pre class="ql-syntax" spellcheck="false">ionq_task_ghz = ionq.run(ghz, shots=1000)

ionq_task_id_ghz = ionq_task_ghz.id
ionq_status_ghz = ionq_task_ghz.state()
print(&#39;Status of quantum task:&#39;, ionq_status_ghz)
</pre> 結構25量子ビットにしては衝撃的な結果が。。。かなり精度いいですね！<pre class="ql-syntax" spellcheck="false">Status of (reconstructed) quantum task: COMPLETED
1000 shots taken on machine arn:aws:braket:us-east-1::device/qpu/ionq/Aria-1.
Measurement counts: Counter({&#39;0000000000000000000000000&#39;: 285, &#39;1111111111111111111111111&#39;: 276, &#39;0000000000000000010000000&#39;: 11, &#39;1111111111111111110000000&#39;: 11, &#39;0000000000000000000000010&#39;: 10, &#39;1111111111110000000000000&#39;: 9, &#39;0111111111111111111111111&#39;: 9, &#39;0000000000000000000111111&#39;: 8, &#39;1111111111111111101111111&#39;: 8, &#39;0000000000111111111111111&#39;: 8, &#39;1111111111111111100000000&#39;: 7, &#39;1111111111111111111000000&#39;: 7, &#39;0000000000000000000001000&#39;: 7, &#39;1111111111111111111111110&#39;: 7, &#39;1111111111111000000000000&#39;: 6, &#39;0000000000000000000000011&#39;: 6, &#39;1111111111111111111011111&#39;: 6, &#39;0000000000000000001111111&#39;: 6, &#39;1111111111110111111111111&#39;: 6, &#39;0000000000001111111111111&#39;: 6, &#39;0000001000000000000000000&#39;: 5, &#39;0000000000001000000000000&#39;: 5, &#39;0000000000000000100000000&#39;: 5, &#39;0000000000000000000000001&#39;: 5, &#39;1111111111111111110111111&#39;: 5, &#39;1111111111111111011111111&#39;: 5, &#39;1111111111111110111111111&#39;: 5, &#39;0000100000000000000000000&#39;: 4, &#39;0000010000000000000000000&#39;: 4, &#39;0000000000100000000000000&#39;: 4, &#39;0000000000000000001000000&#39;: 4, &#39;0000000000000000000100000&#39;: 4, &#39;0000000000000000000010000&#39;: 4, &#39;0000000000000000000000111&#39;: 4, &#39;1111111111111111111110111&#39;: 4, &#39;1111111111111111111101111&#39;: 4, &#39;1111111111111101111111111&#39;: 4, &#39;1111111111111011111111111&#39;: 4, &#39;0000000000000111111111111&#39;: 4, &#39;1111111111011111111111111&#39;: 4, &#39;1101111111111111111111111&#39;: 4, &#39;0011111111111111111111111&#39;: 4, &#39;0100000000000000000000000&#39;: 3, &#39;1110000000000000000000000&#39;: 3, &#39;0001000000000000000000000&#39;: 3, &#39;0000000010000000000000000&#39;: 3, &#39;1111111111000000000000000&#39;: 3, &#39;1111111111100000000000000&#39;: 3, &#39;0000000000000100000000000&#39;: 3, &#39;0000000000000001000000000&#39;: 3, &#39;1111111111111111111111000&#39;: 3, &#39;1111111111111111111111100&#39;: 3, &#39;0000000000000000000011111&#39;: 3, &#39;0000000000000000011111111&#39;: 3, &#39;1111111111101111111111111&#39;: 3, &#39;0000000000011111111111111&#39;: 3, &#39;0000000001111111111111111&#39;: 3, &#39;1111110111111111111111111&#39;: 3, &#39;1110111111111111111111111&#39;: 3, &#39;1011111111111111111111111&#39;: 3, &#39;1100000000000000000000000&#39;: 2, &#39;1111110000000000000000000&#39;: 2, &#39;0000000100000000000000000&#39;: 2, &#39;0000000000010000000000000&#39;: 2, &#39;0000000000000010000000000&#39;: 2, &#39;1111111111111111000000000&#39;: 2, &#39;0000000000100000010000000&#39;: 2, &#39;0111111111111111111000000&#39;: 2, &#39;0000000000000000000000100&#39;: 2, &#39;0111111111111111111111101&#39;: 2, &#39;1111111111111111111111101&#39;: 2, &#39;0000000000000000010000011&#39;: 2, &#39;1111111111110000000000111&#39;: 2, &#39;1111111111111111011011111&#39;: 2, &#39;0000000000000001111111111&#39;: 2, &#39;1111111101111111111111111&#39;: 2, &#39;1111111011111111111111111&#39;: 2, &#39;1111011111111111111111111&#39;: 2, &#39;1000000000000000000000000&#39;: 1, &#39;0010000000000000000000000&#39;: 1, &#39;1111111100000000000000000&#39;: 1, &#39;1000000010000000000000000&#39;: 1, &#39;1110111110000000000000000&#39;: 1, &#39;1111111110000000000000000&#39;: 1, &#39;0000000001000000000000000&#39;: 1, &#39;1111101111000000000000000&#39;: 1, &#39;0001000000100000000000000&#39;: 1, &#39;1101111111100000000000000&#39;: 1, &#39;0010000000010000000000000&#39;: 1, &#39;0000000001111100000000000&#39;: 1, &#39;1111111111111100000000000&#39;: 1, &#39;1000000000000010000000000&#39;: 1, &#39;0100000000000010000000000&#39;: 1, &#39;0000000001000110000000000&#39;: 1, &#39;0000001000000001000000000&#39;: 1, &#39;0000000000001111000000000&#39;: 1, &#39;1000000000000000100000000&#39;: 1, &#39;1000000010000000100000000&#39;: 1, &#39;1111111111111101100000000&#39;: 1, &#39;0000000000011111100000000&#39;: 1, &#39;0000001000000000010000000&#39;: 1, &#39;0000000010000000010000000&#39;: 1, &#39;0000000000010000010000000&#39;: 1, &#39;1111111111111010110000000&#39;: 1, &#39;0000000000001111110000000&#39;: 1, &#39;0001111111101111110000000&#39;: 1, &#39;1111111111101111110000000&#39;: 1, &#39;0000000000011111110000000&#39;: 1, &#39;0000000000000001001000000&#39;: 1, &#39;0000000000000000011000000&#39;: 1, &#39;0000000000001000011000000&#39;: 1, &#39;1111111111111111011000000&#39;: 1, &#39;1111111111101111111000000&#39;: 1, &#39;0001111111111111111000000&#39;: 1, &#39;0000000000000100000100000&#39;: 1, &#39;0000000000000000100100000&#39;: 1, &#39;1111111111111111111100000&#39;: 1, &#39;1111111111111111111110000&#39;: 1, &#39;1100000000001000000001000&#39;: 1, &#39;0000000000000010000001000&#39;: 1, &#39;0000111111111111110001000&#39;: 1, &#39;0101000011001011111101000&#39;: 1, &#39;0000000000000000000011000&#39;: 1, &#39;0000000000000000111111000&#39;: 1, &#39;0000000000111111111111000&#39;: 1, &#39;0111111111111111111111000&#39;: 1, &#39;0000000000000100000000100&#39;: 1, &#39;1111111111111111111101100&#39;: 1, &#39;1111111111110110111111100&#39;: 1, &#39;1111111111111110111111100&#39;: 1, &#39;1000000000000000000000010&#39;: 1, &#39;0000000001000000000000010&#39;: 1, &#39;1011111111111110100000010&#39;: 1, &#39;1111111111111111110000010&#39;: 1, &#39;0000000001111111111111010&#39;: 1, &#39;0000000000011000000000110&#39;: 1, &#39;0000000000000010000001110&#39;: 1, &#39;0000000010000000000011110&#39;: 1, &#39;0000000000000001000011110&#39;: 1, &#39;0000000000000000011111110&#39;: 1, &#39;0000000000000111111111110&#39;: 1, &#39;0001000000000000000000001&#39;: 1, &#39;1111111111100000000000001&#39;: 1, &#39;1101111111110000000000001&#39;: 1, &#39;1111111111111010000000001&#39;: 1, &#39;1111111111111111111001101&#39;: 1, &#39;1111111111111111110111101&#39;: 1, &#39;0000000000001111111111101&#39;: 1, &#39;1111111111111000000000011&#39;: 1, &#39;1111111111111111111000011&#39;: 1, &#39;0000000000011111111110011&#39;: 1, &#39;1111111111111111010111011&#39;: 1, &#39;0000000111111111101111011&#39;: 1, &#39;1111111111111111111111011&#39;: 1, &#39;0010000000000000000000111&#39;: 1, &#39;0000000010000000000000111&#39;: 1, &#39;1111111111100000010110111&#39;: 1, &#39;0000000000000000011110111&#39;: 1, &#39;1111111111101111111110111&#39;: 1, &#39;1011111111111111111110111&#39;: 1, &#39;0000000000000000000001111&#39;: 1, &#39;0010000000000000000001111&#39;: 1, &#39;0000000000000000000101111&#39;: 1, &#39;1111111110111111111101111&#39;: 1, &#39;1111111111111111110011111&#39;: 1, &#39;1101111111111111111011111&#39;: 1, &#39;1000000000000000000111111&#39;: 1, &#39;0001100000000000000111111&#39;: 1, &#39;1111111111111000000111111&#39;: 1, &#39;1100000000000000001111111&#39;: 1, &#39;1111111111011100001111111&#39;: 1, &#39;1111111111111110101111111&#39;: 1, &#39;0000000000000111101111111&#39;: 1, &#39;0000000001000000011111111&#39;: 1, &#39;1111111111000000011111111&#39;: 1, &#39;0000000000001000111111111&#39;: 1, &#39;0000000000011101111111111&#39;: 1, &#39;0000000000011011111111111&#39;: 1, &#39;1111110000000111111111111&#39;: 1, &#39;0000000001110111111111111&#39;: 1, &#39;0000001111110111111111111&#39;: 1, &#39;1111000000001111111111111&#39;: 1, &#39;0000000000101111111111111&#39;: 1, &#39;1111110111101111111111111&#39;: 1, &#39;0000000100011111111111111&#39;: 1, &#39;1111110111011111111111111&#39;: 1, &#39;1001111111011111111111111&#39;: 1, &#39;1000000000111111111111111&#39;: 1, &#39;1111111110111111111111111&#39;: 1, &#39;0111111101111111111111111&#39;: 1, &#39;0000000111111111111111111&#39;: 1, &#39;1111101111111111111111111&#39;: 1, &#39;0000011111111111111111111&#39;: 1, &#39;0001111111111111111111111&#39;: 1})
</pre> <img src="https://assets.blueqat.com/public/uploads/us-east-2:4805ff4b-c3cc-4344-b165-86544c34d0bf/2023/10/13/ghz2.png"/> 両側にきちんとGHZが実行できてるのが見えます。結構革命的にすごいです。よだれが出ますね。。。量子コンピュータは着実に開発が進んでいるようです。。。

最新版のamazon braketで、IonQのAria1プロセッサ25量子ビットを稼働させてみる。

Yuichiro Minato