不確定な世界

科学の話題を中心に、勉強したことや考えたことを残していきたいと思います

「世界を変えた書物」展＠上野の森美術館

その他

３連休の目玉イベントとして、金沢工業大学が主催する「世界を変えた書物」展（上野の森美術館）に行ってきた。
www.kanazawa-it.ac.jp

歴史上の天才科学者たちが後世に残した偉大な書物の数々。その実物をこの目で見ることができるという、科学ファンなら必見のイベント。しかも、入場は無料で会場内は写真撮影OK（ただしフラッシュ禁止）、SNSへの投稿もOK（他の来場者の映り込みに注意）なのである。なんとも太っ腹なことだ。２時間半ほど、じっくり鑑賞させていただいた。
というわけで、さっそく私の撮った本の写真をいくつか紹介しよう。

エウクレイデス（ユークリッド）「原論」、１４８２年、初版

f:id:quanta087:20180923205447j:plain

コペルニクス「天球の回転について」、１５４３年、初版

f:id:quanta087:20180923210035j:plain

ガリレオ・ガリレイ「星界の報告」、１６１０年、初版

f:id:quanta087:20180923210450j:plain

ニュートン「自然哲学の数学的原理（プリンキピア）」、１６８７年、初版*1

f:id:quanta087:20180923211204j:plain

ネーピア「驚くべき対数法則の記述」、１６１４年、初版

f:id:quanta087:20180923211813j:plain

ライプニッツ「極大と極小に関する新しい方法」、１６８４年、初版

f:id:quanta087:20180923212434j:plain

オイラー「無限解析入門」、１７４８年、初版

f:id:quanta087:20180923215556j:plain

デカルト「方法序説」、１６３７年、初版

f:id:quanta087:20180923215826j:plain

ホイヘンス「光についての論考」、１６９０年、初版

f:id:quanta087:20180923220101j:plain

ダニエル・ベルヌーイ「流体力学」、１７３８年、初版

f:id:quanta087:20180923220234j:plain

アグリコラ「金属について（デ・レ・メタリカ）」、１５５６年、初版

f:id:quanta087:20180923215103j:plain

ファラデー「電気の実験的研究 第Ⅰ、Ⅱ、Ⅲ巻」、１８３９－４４、１８５５年、初版

f:id:quanta087:20180923213232j:plain

マイケルソン & モーリ「地球と発光性エーテルの相対運動について」、１８８７年、初版

f:id:quanta087:20180923214652j:plain

アインシュタイン「特殊相対性理論及び一般相対性理論」、１９１７年、初版

f:id:quanta087:20180923214024j:plain

ダーウィン「種の起源」、１８５９年、初版

f:id:quanta087:20180923214225j:plain

ここに載せたのは展示のほんの一例である。
この他にもウィトルウィウス、アルキメデス、アポロニウス、ケプラー、カルダーノ、フック、ヤング、ゲーテ、ボイル、ラヴォアジェ、クーロン、エジソン、ガルヴァーニ、ヴォルタ、アンペール、オーム、グラハム・ベル、ヘルツ、マルコーニ、ローレンツ、マクスウェル、レントゲン、キュリー、シュレディンガー、ラザフォード、湯川秀樹、プランク、リーマン、ミンコフスキー、パスカル、ヘルムホルツ、メンデルなどなどなどなど。
およそ科学のあらゆる分野の偉人による書籍や論文が目白押しであった。古今東西の天才科学者が生み出した知性の結晶に囲まれるという、一生の思い出になった。

最後に、本企画展の関連書籍として、以下の本を紹介しよう*2。

図説 世界を変えた書物 科学知の系譜
posted with ヨメレバ
竺 覚暁 グラフィック社 2017-11-17
Amazonで購入
Kindleで購入

*1:本展示会で最も楽しみにしていた本の一つ。言うまでもなく、物理学を学んだ人間にとって聖書のような存在である。

*2:会場の物販に置いてあるだろうからそこで買おうと思っていたのだが、置いていなかったので帰宅後にポチった。

「nature 科学」シリーズ；竹内薫 監修　読書感想

読書

今回紹介するのは、
「nature 科学　系譜の知」
「nature 科学　未踏の知」
「nature 科学　深層の知」
の３冊。
「nature」は有名な科学論文雑誌なので、ご存じの方も多いだろう。そのnature誌の一般向け解説コラム「News & Views」過去１０年分の中からイチオシの内容が選定され、分野別に収録されているのが本書である。

本書は分野別に３冊に分かれているのだが、それぞれの巻の話をする前にまずは全体の感想を述べておこう。
本書の内容はかなり難しい。ブルーバックスのような本と違って特定の論文をダイレクトに解説しているため、細分化されたニッチな話が多いし、専門用語もバンバン出てくる。正直に言うと、自分の専門分野に比較的近い部分を除き、書いてあることの殆どは理解できなかった。
それでも私が本書を読んだのは、科学の「現場」に近い情報が欲しかったからだ。素人向けの雑学本やニュース記事から学ぶことも多くあるが、やはりそのような知識と最先端の現場にはある程度のギャップがある。常識的な手段*1でそのギャップを埋める最も手っ取り早い手段は学術論文を読むことだが、大抵の人にとってそれはハードルが高すぎる。本書を読めば、学術論文を読んだり、あるいはそれを理解するための体系的な勉強をすっ飛ばして、現場を覗き見ることができるというわけだ。万人のおすすめできる本ではないが、私と同じように分からないなりに現場を知りたい、という人には一読の価値がある。

さて、それでは、個々の巻の中で、興味を持ったテーマをそれぞれ１つずつピックアップしてみよう。

系譜の知

nature科学 系譜の知 バイオ(生命科学)、医学、進化(古生物)
posted with ヨメレバ
竹内 薫 実業之日本社 2015-01-31
Amazonで購入
Kindleで購入
この巻では主にバイオ・医学・生物系の内容が扱われている。私の生物学の知識は高校１年生の時に必修だった生物Ⅰで止まっており、特にタンパク質やら薬品やらの名前が次々に出てくる分子生物学系の話には苦手意識さえある。一方で、脳科学には人並みに興味があるし、進化論・古生物学の話は割と好きだ。
そういうわけで、この巻で最も興味を持った話題は「進化（古生物）」カテゴリの、「始祖の地位から墜ちた始祖鳥」（オリジナルの論文は文献[1]を参照）。始祖鳥といえば恐竜から鳥への進化を表す化石として有名であるが、なんと、その始祖鳥が鳥ではないかもしれないというのだ。もしこれが本当なら、冥王星が惑星から墜ちたときと同じくらいの衝撃である*2。
おおざっぱに要約すると、中国で新たに発見された新種を含めて分析すると、始祖鳥は鳥群（現在の鳥類に繋がる分類）から外れて別系統の種として分類される、ということらしい。ただし、重要なことは、その新種を除外して再分析すると、始祖鳥は鳥群に復活するということだ。たった1種の有無だけで分類が変わってしまうわけなのだから、今後の発見次第でこの研究の結果も覆る可能性は十分にありうるということだ。なんとも難しい分野である。

参考文献

[1] Xing Xu et al., An Archaeopteryx-like theropod from China and the origin of Avialae. Nature 475, 465–470 (2011)

未踏の知

nature科学 未踏の知 地球(地球文明)・環境・宇宙
posted with ヨメレバ
竹内 薫 実業之日本社 2015-01-31
Amazonで購入
Kindleで購入
地球・環境・宇宙系の内容が扱われている巻ではあるが、宇宙に関する内容がおよそ半分を占めている。一般に宇宙の話が一番需要があるのだろうが、ここではあえて「地球」カテゴリから、「シート状の対流が作る地球ダイナモ」を紹介したい。
そもそもダイナモ理論とは、地球などの天体が磁場を生み出す仕組みを、天体内部の導電性流体の対流に求める理論である。ここで紹介されている研究は、従来のダイナモ理論のシミュレーションでは流体の対流が円柱状になるのに対し、より現実に近い条件でシミュレーションを行ったところ、対流が薄いシート状になったという内容だ。
私の印象に残った部分は二つ。一つは条件が現実に近いにもかかわらず、棒磁石のような双極子型の磁場形状を説明するという観点からは、従来の理論よりも悪くなってしまったらしいということ。もう一つは、コンピュータの性能限界から、その「現実に近い条件」でさえ非現実的にせざるを得ないということだ。例えば、粘性率に関わるエクマン数という数値は、現実の10^{-15}ではなく、10^{-6}\sim10^{-7}程度の数値を使っているというのだ。なんと現実値と８桁も違う！ほかにも、磁気プラントル数という数値は現実値10^{-6}に対して、1となっているらしい。この点について記事の中では、「モデル作成の秘訣は「適当な割合で間違っている」変数を用いることにあるのかもしれない」と指摘されている。つまり、模型飛行機で風洞実験を行うのと同じように、適切なスケーリングで計算しなくては現実は再現できない、ということなのだろうと理解した。
この記事の元になった論文（文献[2]）は２００８年のもの。この１０年で進展があったのかどうか気になって少し調べてみた。文献[3]（２０１５年）に目を通す限りそこまで大きなブレイクスルーはないようだが、地道に進んでいる様子がうかがえる。今後の研究に期待したい。

参考文献

[2] Kageyama, A., Miyagoshi, T. & Sato, T., Formation of current coils in geodynamo simulations. Nature 454, 1106–1109 (2008)
[3] 宮腰剛広, 陰山聡 「地球ダイナモ研究のこれまでとこれから」, J. Plasma Fusion Res. 91, 676-682 (2015)

深層の知

nature科学 深層の知 物理数学・物理化学・工学・ロボット
posted with ヨメレバ
竹内 薫 実業之日本社 2015-01-31
Amazonで購入
Kindleで購入
この巻は学生時代の専門に最も近い物理・化学・工学分野を扱っており、上の２冊に比べて比較的理解しやすく、興味を持った記事も多かった。その中から一つ選ぶとしたら、やはり量子コンピュータ関連の記事「ピンクのキュービット」だろう。
ダイヤモンドと言えばもちろん宝石としての用途が最も一般的だが、その硬さを活かして研磨剤やダイヤモンド・カッターなどの加工機械に使われていることも有名である。しかし、その輝きの中に量子コンピュータとしての機能も備わっていることは、あまり知られていないのではないだろうか。
ダイヤモンド結晶の中で窒素不純物と空孔が隣り合った構造を窒素-空孔（NV）中心という。NV中心はダイヤモンドにピンク色を帯びさせるという、宝石の観点からも重要な特性を持つ。しかし今重要なのは、NV中心が天然の量子ドットとなって電子を捕獲し、その電子のスピンをキュービットとして利用することができるという点だ。
ダイヤモンド中のスピンは驚異的なコヒーレンス時間を持つという優れた性能を持つが、一方で拡張性、つまり多キュービット系を構築することに対して課題がある。その課題にチャレンジしたのがこの研究で、NV中心と、その近くに存在する別の電子スピン（窒素中心）が２キュービット系として扱えることを明らかにしている。
ただし、この記事で紹介されている論文（文献[4]）は２００６年のものであり、この分野の進展スピードを考えるとかなり古い。現在、ダイヤモンドキュービットでは既に量子テレポーテーションや量子エラー訂正アルゴリズムが実装されている（文献[5][6]）*3。また、ダイヤモンドは超伝導回路などの他のハードウェアと比べて「光」との相性がよいため、量子通信関連の文脈でよく登場する。実をいうと、この記事を書いている途中でまさにそのようなニュースが飛び込んできた（文献[7]）*4。
最先端の科学というは、今、この瞬間にも進んでいるのである。

参考文献

[4] R. Hanson et al., Polarization and Readout of Coupled Single Spins in Diamond. Phys. Rev. Lett. 97, 087601 (2006).
[5] W. Pfaff et al., Unconditional quantum teleportation between distant solid-state qubits. Science 345, 532–535 (2014).
[6] T. H. Taminiau et al., Universal control and error correction in multi-qubit spin registers in diamond. Nat. Nanotechnol. 9, 171-176 (2014).
[7] Peter C. Humphreys et al., Deterministic delivery of remote entanglement on a quantum network. Nature 558, 268–273 (2018).

最後に

最後にもう一つ、本書の内容そのものとは直接は関係ないのだが、本書のおすすめポイントとしてぜひ挙げておきたいのが、「帯」のデザインである。
f:id:quanta087:20180611230412j:plain

光の３原色をモチーフにした色合いが美しい。特に私のお気に入りは「未踏の知」の青色だ。ぼぉっと眺めていると、宇宙を漂いながら世界の深淵を覗いているような、神秘的な気分に浸ることができる。私は普段はあまり本の帯の有無を気にする性格ではないのだが、本書に限って言えば内容以上に帯に惹かれて手に取ったと言っても過言ではない。もし本書をこれから読む方がいれば、ぜひ３冊そろえて並べてみて欲しい。

*1:研究室に連絡を取ってラボ見学させてもらうとか、twitterをやってる大学教授をフォローして現場の情報を流してもらうとかいう方法も考えられるが、そんなコミュ力お化けのような真似は私にはできない。

*2:冥王星の話はかつて「子供の科学」を読んでいた時に知ったのだが、私はその時に「科学的知識は固定されたものではなく、定説が変わることもありうる」ということを学んだ。

*3:少し専門的になるので脚注とするが、現在では多キュービット化のためにダイヤモンドに一定確率で含まれる炭素同位体の核スピンを用いるのがトレンドである。

*4:gigazineの記事ではダイヤモンドとは書いていないが、論文や英語記事の方を読めばダイヤモンドを用いていることがわかるだろう

「ディープラーニングがわかる数学入門」；涌井良幸・涌井貞美 著　読書感想

人工知能 読書

本日紹介するのは、涌井良幸・涌井貞美 著「ディープラーニングがわかる数学入門」。

ディープラーニングに関してはこれまでに「ゼロからわかるDeep Learning （感想記事）」「機械学習と深層学習-C言語によるシミュレーション- （感想記事）」など、実装中心の本を読んでいたが、そろそろ理論的な理解も深めたいと考えていた。ディープラーニングの理論的教科書と言えば青本が有名だが、書店で立ち読みして自分にはまだ早いと感じ、代わりに手に取ったのが本書である。

著者のお二人は人工知能の専門的研究者ではなく、数学の（元）高校教員だ。退職後、数学や統計学を中心に数多くの入門書や解説書を執筆されており、元高校教員という経歴からか初心者への分かりやすさに定評がある。本書でも、ベクトルや微分、行列のような基本的な数学から始まり、ディープラーニングの本質がわかりやすく解説されている。

特に初心者向けだと思った本書の特徴は、記号の定義にページが割かれていることだ。ニューラルネットワークの理論では、扱わなければならないパラメータの数が非常に多い。しかも、一つの変数にこれでもかというくらい添え字がくっついていて、一目見ただけでは意味が分からない。学術書と呼ばれるレベルの本では記号の定義は定義式を一行示すだけで済ますことが多いが、本書では実に１０ページ近くを割いて変数と添え字が何を指しているのかを図と一緒に解説してくれている。私も添え字が何を指しているのかすぐに混乱してしまうので、このような丁寧な説明は非常にありがたい。
また、「注がしつこい」ことも挙げておこう。誤差関数の定義に1/2がついている理由、フィルターサイズを一般的な5x5ではなく3x3にしている理由など、細かい点に関する注が繰り返し書いてある。わかっている人には鬱陶しいと感じられるかもしれないが、初心者は大枠についていくのに必死で細かい注のことなど何回聞いても忘れてしまうものだ。このようなしつこさも教員の性なのだろうか。

私自身が学んだこととして挙げられるのは、まず勾配降下法が「最急」であるという事実を再認識できたこと。数学を避けたニューラルネットワークの解説では、勾配降下法の更新量を決める公式
\displaystyle \Delta\omega = -\eta\nabla E
は天下り的に与えられることが多いが、本書ではそもそもの勾配ベクトルの導出までさかのぼり、この公式が最も効率よく誤差を減らす（最急である）ことを説明している。勾配ベクトルが「最急」であること自体は知識として知っていたが、”ベクトルの内積が最も小さくなるのはベクトルが逆向きのときである”という性質から勾配ベクトルを導出することで、その意味がストンと腑に落ちた。もっとも、改めて線形代数やベクトル解析の教科書を開いてみると、導出過程まできちんと載っていたわけで、単に自分が不勉強だっただけなのだが…。

もう一つ、誤差逆伝播法についての理解が深まったことも大きい。「ゼロから作るDeep Learning」では、誤差逆伝播の説明として数学を用いるものと計算グラフを用いるのもを挙げ、後者を採用していた。本書ではもちろん前者、数学を用いた説明である。
以前計算グラフを学んだときにも、数値微分をするより計算効率がいいことに一応の納得はできた。しかし、色々なところで「誤差逆伝播といえばチェーンルール」というような説明ばかり聞いていたため、「合成関数を微分するのにチェーンルールを使うのは当たり前ではないか。何がそんなに凄いんだ」という思いは拭い切れていなかった。本書を読んだことでその疑問に自分なりの答えが出た。凄いのはチェーンルールそのものではなく、チェーンルールを微分を含まない漸化式に変形したところだったのだ。微分というコストの高い計算を避けて積和のみで計算できること、そして後ろの計算をするために前の計算結果をそのまま使うことができること。これこそが、単純にチェーンルールに則って数値微分を行うのに比べて誤差逆伝播法の計算効率が高い理由なのだ。
チェーンルールから漸化式を導く計算過程はそれほど自明ではない（少なくとも私には思いつかない）。きっと、式変形が得意な数学の天才が試行錯誤の末にたどり着いたのだろう。「誤差逆伝播がなぜ凄いのか」、この点を理解できたのが本書を読んでの最大の収穫である。

最後に、注意点もいくつか挙げておこう。
まず、本書で扱われている最適化手法は、用意した全学習用データを使って誤差を計算するバッチ勾配降下法（最急勾配法）である。実際には、全データの中からランダムにサンプリングしたデータに対して最適化を行う確率的勾配降下法が用いられる*1。また、畳み込みニューラルネットワークの説明では、ストライドが変数になっていないなど、やや一般性に欠ける部分がある。本書でも初心者がディープラーニングについて知っているべき本質は十分に学ぶことができるが、実践レベルに到達するためには、もう少し高度な教科書を読む必要があるだろう。私も、いずれはあの青いシリーズを理解したいものである。


ディープラーニングがわかる数学入門
posted with ヨメレバ
涌井 良幸,涌井 貞美 技術評論社 2017-03-28
Amazonで購入
Kindleで購入

*1:確率的勾配法の中でも、サンプリング数が１つであるものをオンライン学習、ある程度まとまった数をサンプリングする場合をミニバッチ学習と呼ぶ、らしい

「機械学習と深層学習-C言語によるシミュレーション-」；小高知宏 著　読書感想

読書 人工知能

年始に買って積読していた、小高知宏著「機械学習と深層学習-C言語によるシミュレーション-」を読み終わった。
本書の特徴は、事実上Pythonに支配されている機械学習や深層学習のプログラムをC言語で学べること。しかも、特殊なライブラリに頼ることなく、きわめて基本的な文法だけで１ファイルに収まるようにコーディングされている。例えば、畳み込み処理などもすべてforループで愚直に行っている。以前読んだ「ゼロから作るDeep Learning」でさえ、numpyはフル活用しているし、後半は著者作成のライブラリを導入してしまっていたことから考えると、この方針を一切ぶれずに貫いているのは素晴らしいと思う。ただ、その分本書で扱っている例題はごく簡単なものが多く、実際の問題にすぐに応用できるものではないのも事実ではある。

サンプルコードとデータセットは本書の公式ページからダウンロードすることができる*1。単にサンプルを実行したり写経するだけでは頭に入ってこない人は、自分の得意な言語に移植してみるのも一興だろう。私もとりあえずCで実行して動くことを確認してからPythonで書いてみる、という進め方をした。特別うまく書かなくてもコード量が半分程度になり、改めてPythonの便利さを実感した次第である。

さて、本書で学べる手法は以下の通り：
・帰納的学習（ランダムパターン生成）
・強化学習（Q学習）
・蟻コロニー最適化法
・遺伝的アルゴリズム
・パーセプトロン
・深層学習（畳み込みニューラルネットワーク）

見ての通り、機械学習の中でも生物学的手法を多く扱っている。この中で私にとって予想外の収穫だったのは遺伝的アルゴリズムが扱われていたことである。遺伝的アルゴリズムのことは以前から知っていた。きっかけは確かYoutubeの動画で、ヒト型の３Dモデルに二足歩行を教えるものと、マリオをクリアさせるものであった。

マリオの動画は残念ながら削除されていた…。まあ、ともかくそれ以来遺伝的アルゴリズムにはずっと興味を持ちつつも専門的に勉強する機会はなかった。今回思いがけず遺伝的アルゴリズムを実装ベースで勉強することができたのは私にとって大きな一歩であった。

話が逸れてしまった。本書は以上のように幅広い手法を扱っており、機械学習の素人である私にとってはどれも新鮮に学ぶことができた。しかし、ニューラルネットワークや深層学習をメインに学びたい人にとって本書の内容は決して満足できるものではないだろう。深層学習メインの方には、正直言っておすすめできない。だが、畳み込みなどを本当に一から実装している点などには何か得るものがあるのではないだろうか。
もう一つマイナス面として、理論的な面ではほとんど参考にならない、ということも注意しておこう。もともと約２００ページという比較的薄目の本であるが、その半分はサンプルコードで占められている。扱っている内容の多さもあり、一テーマあたりの理論的解説は本当に微々たるものである。本書はあくまでもコードを通して理解するタイプのものであることを忘れてはならない。

結論。深層学習に限らずに生物学的機械学習を幅広く学びたい方、C言語は使えるのにわざわざPython覚えなくちゃいけないの…？という方、上級者様の「エレガントな」プログラミングテクニックについていけず、もっとシンプルな文法だけで機械学習を勉強したい方 etc...　にはぜひ本書をお勧めしたい。

機械学習と深層学習 ―C言語によるシミュレーション―
posted with ヨメレバ
小高知宏 オーム社 2016-05-21
Amazonで購入
Kindleで購入

*1:このサンプルコードを利用できるのは本書を購入した方のみ、ということが規約に書いてあるので注意して欲しい

量子コンピュータの基本素子・量子ビットのハードウェア実装（超伝導磁束編その５～データの読み出し（測定）～）

量子テクノロジー

その１～素子構造～
その２～超伝導リング詳細～
その２.５～ノイズ耐性～
その３～初期化～
その４～データの書き込み・演算


では最後に、量子ビットに書き込んだデータを読み出す技術、つまり「量子を観測する」とは実際にはどういうことなのか、を見ていこう。かなり難しく・長くなるので注意して欲しい。
今回キーとなるのは図１９で量子ビット本体を囲んでいる②の回路、SQUIDだ。SQUIDは必ずしも量子コンピュータ専用の部品というわけではなく、一般的に超高感度磁場センサとして使われているので、すでに知っている方もいるかもしれない。SQUIDの動作原理を数式で理解しようとすると結構大変だが、ここでは量子ビット読み出しの仕組みを理解するのに必要な部分だけをおおざっぱに説明する。

f:id:quanta087:20171222074806p:plain
図１９　回路全体図（文献[1]より転載）

SQUIDの仕組み

SQUIDの構造

SQUIDの構造自体は量子ビット本体とよく似ており、ジョセフソン接合を含む超伝導リングである。量子ビットと違うのは、ジョセフソン接合が２つであること、外部から電流を注入するための導線がくっついていることである（図２０）。

f:id:quanta087:20171225232026p:plain
図２０　SQUIDの回路図

ジョセフソン接合の限界

ジョセフソン接合は、薄い絶縁体を超伝導体でサンドイッチすると、絶縁体にも超伝導電流が流れるというものであった。しかし、このような状態にはやはり限界がある。ジョセフソン接合に流す電流値がある限界を超えると、超えた分の電流は普通の（トンネル）電流として流れてしまうのだ*1*2。つまり、ジョセフソン接合に電気抵抗が生じ、電圧降下が発生する。このことがSQUIDの動作原理のキーポイントである。

SQUIDによる磁場検出

では、一般的にSQUIDでどのように磁場を検出するのかを見ていこう。なお、以下の説明の大部分は文献[2]（こちらのウェブページ）を参考にさせていただいた。
まず具体的にどういう操作をするのか結論だけ述べておく。SQUIDに外部から電流を注入し、それを強くしていって、どのくらいの電流でジョセフソン接合が限界に達して電圧が生じるのかを観察するのである。

（１）磁場がない（ゼロ）場合

最も簡単な場合として、検出すべき磁場がゼロである場合のことを考えよう。図２１のようにSQUIDに生えている端子から電流\displaystyle I_{in}を注入すると、SQUID内を\displaystyle \frac{I_{in}}{2}ずつに分かれて流れ、もう片方の端子で合流して出ていく。\displaystyle I_{in}をどんどん大きくしていって、ジョセフソン接合を通る電流\displaystyle \frac{I_{{in}}}{2}が限界値を超えるとジョセフソン接合の両端に電圧が生じるので、あらかじめSQUIDに電圧計をつけておいて電圧計が反応した時の\displaystyle I_{in}を記録しておく。ただし、これはあくまで一般論である。量子ビットに応用するときには電圧計は別のところについているなど、ちょっとした違いがあるので、その点については後述する。

f:id:quanta087:20171226080744p:plain
図２１　磁場がない場合

（２）弱い磁場がある場合

次に、SQUIDリングを弱い磁場が下から上に貫いている場合を考えよう*3。その２で述べたように、超伝導リング内に中途半端な磁場が存在すると、それを相殺する磁場を生み出すために、リングに（右回りの）電流が流れる。この状態で電流\displaystyle I_{in}を注入してみよう（図２２）。

f:id:quanta087:20171226205740p:plain
図２２　弱い磁場がある場合

ここで、ジョセフソン接合に流れる電流に注目して欲しい。右側では、元々SQUIDに流れている電流と外部から注入された電流\displaystyle \frac{I_{{in}}}{2}が、同じ向きに流れているのだ。キルヒホッフの法則により、ジョセフソン接合にはこの２つの電流の和が流れ込むことになる。そのため、\displaystyle I_{in}を強くしていくと、磁場がなかった時よりも小さい\displaystyle I_{in}でジョセフソン接合が限界を迎えるのだ*4。

（３）そこそこ強い磁場がある場合

続いて、もう少し強い（具体的には\displaystyle \frac{\Phi _0}{2}を超えない程度の）磁場がSQUIDを貫いているとしよう。加わっている磁場が大きいと、それを相殺するために流れる電流も大きくなる。その結果、ジョセフソン接合に電圧降下を起こすために必要な外部電流はさらに小さくなるのがわかるだろう（図２３）。

f:id:quanta087:20171226220313p:plain
図２３　そこそこ強い磁場がある場合

このように、SQUIDを貫く磁場が大きくなると、電圧降下を起こすのに必要な外部電流が小さくなっていく。さらに正確に言えば、磁場が\displaystyle \frac{\Phi _0}{2}を超えると再び外部電流の限界値が大きくなり、磁場が\displaystyle \Phi _0になると磁場がないときの状態に戻る、という周期性がある（磁場が\displaystyle \Phi _0であれば磁束量子条件を満たしているので、SQUIDに円電流は流れない）。この効果を利用して磁場の強さを測定するのが、SQUIDの基本的な原理である。ただし、文献[2]として挙げたウェブページでは、注入する外部電流は常に一定で、常伝導成分の変化によって生じる電圧値の変化を読み取る、という説明になっている。ここでは、量子ビットの説明に都合がいいように説明を改変したことをお断りしておく。

ここまでのまとめ

SQUIDについて、ここから先に必要になる部分をまとめておこう。最低でも、以下の２点だけは押さえて欲しい。
・ジョセフソン接合に限界を超えた電流を流すと、抵抗および電圧降下が生じる。
・SQUIDを貫く磁場が大きいほど、より小さな外部電流を注入するだけで電圧降下が起きるようになる（正確には、\displaystyle \frac{\Phi _0}{2}を境にして\displaystyle \Phi _0で元に戻るという周期性がある）。


量子ビット読み出し回路としてのSQUID

では、SQUIDがどのように量子ビットに応用されているのかを見ていこう。ただ、ここの部分の説明は比較的優先度が低いので、どうしても我慢できなければ「ここまでのまとめ」を読んでいただくだけで構わない。

量子ビットが生み出す磁場

再び図１９の写真を見ると、SQUIDの内部に量子ビットが配置されていることがわかる。量子ビットから発生する磁場をSQUIDで測定するわけだ。ところで、素朴に写真を眺めていると、SQUIDは量子ビットの内側部分の磁場を検出のだろう、となんとなく直感的に考えると思う（図２４）。

f:id:quanta087:20171228145028p:plain
図２４　内側の磁場を検出？

しかしよく考えてみると、量子ビットの内側の磁場は、まさに量子ビットが生み出す磁場によってゼロまたは\displaystyle \Phi _0に保たれている。上述したように、\displaystyle \Phi _0の磁場はSQUIDにとってはゼロと同じことであり、ゼロと区別することができない。ではどうするのかというと、実は、SQUIDが検出するのは量子ビットの内側ではなく、外側の磁場なのである（図２５）*5。

f:id:quanta087:20171228152514p:plain
図２５　外側の磁場を検出

一応、量子ビットの外側に磁場が生じる仕組みも説明しておこう。量子ビットリングには、加えられている磁場を相殺するための右回りの電流が流れている。ここで、リングの一辺だけを考えて直線状の電流に対する「右ねじの法則」を適用する。すると、リングの内側では磁場が下向きであるのに対し、リングの外側では磁場が上向きであることがわかるだろう（図２６）。

f:id:quanta087:20171228173823p:plain
図２６　リングの外側に磁場が生じる

このように、量子ビットはリングの内側だけでなく、外側にも磁場をつくる。SQUIDはこの外側の磁場を検出するわけだ*6。ただし、量子ビットが内側につくる磁場が、例えば\displaystyle 0.1 \Phi _0だったとして、外側も\displaystyle 0.1 \Phi _0だとは限らないことは注意しておく。磁場（正確には磁束）の値は、領域の面積に比例するからだ。図１９の写真を見ると、どう考えても明らかに量子ビットの外側の面積の方が小さいので、具体的な磁場の数値は外側では小さくなるはずだ。とは言っても、「リングに流れる電流が弱ければ生じる磁場も弱いし、リングに流れる電流が強ければ生じる磁場も強い」という定性論は成り立つ。ここで思い出してもらいたいのが、その４の後半で述べた通り、データ読み出しの段階では量子ビットは「右回り」と「左回り」が基準になる状態に変換されており、右回り状態が電流が弱く（＝エネルギーが低く）、左回り状態は電流が強い（＝エネルギーが高い）ということだ（その２の図６を参照）。つまり、量子ビットが外側につくる磁場にも強弱がある。量子ビットが右回りなのか左回りなのかによってSQUIDを貫く磁場の強さが変わり、その違いが検出されるのだ。

ここまでのまとめ

・SQUIDは、量子ビットが（リングの外側に）つくる磁場の強弱を読み取る。
・量子ビットが右回り状態であればSQUIDを貫く磁場は弱く、左回り状態であればSQUIDを貫く磁場は強い。

量子ビットのデータ読み出し

ゴールまであともう少しだ。
ここで、入出力線まで含めた量子ビット読み出し回路の全体図を確認しておこう（図２７）。ただし、図は文献[1]を参考に説明しやすいように改変した*7。

f:id:quanta087:20171228201040p:plain
図２７　量子ビット読み出し回路全体図

SQUIDの説明でもちょっと触れたが、測定器（電圧計＋カウンタ・パソコンなど）がSQUID本体とは少し離れて設置されている。そのため、信号の流れがやや複雑になるのだが、それを理解するためにあと一つだけ勉強していただくことがある。中学校の理科（あれ、もしかしたら高校の物理かも？）で習う回路理論だ。

並列抵抗回路

問題：以下の図２８のように、抵抗\displaystyle R_1と抵抗\displaystyle R_2が並列に接続されている。そこに電流\displaystyle Iを流すと、それぞれの経路を流れる電流はいくらか？
答え：\displaystyle I_1 = \frac{R_2}{R_1 + R_2}I\qquad \displaystyle I_2 = \frac{R_1}{R_1 + R_2}I

f:id:quanta087:20171228210957p:plain
図２8　並列抵抗

一気に答えまで出してしまったが、並列回路では抵抗値に応じて電流が分配されるということが分かればよい。例えば、\displaystyle R_1=1 \Omega、\displaystyle R_2=2 \Omegaとすると、\displaystyle I_1=\frac{2}{3}I、\displaystyle I_2=\frac{1}{3}Iである。両方の経路に電流が流れてはいるが、その大きさは異なる。抵抗とは電流の流れにくさであるのだから当然だが、電流は抵抗が低い経路を優先して流れるのだ。このことをもう少し極端な状況で試してみよう。\displaystyle R_1=0 \Omegaにしてみたらどうだろうか。上記の式に当てはめると、\displaystyle I_1=I、\displaystyle I_2=0となる。つまり、抵抗ゼロの経路がある場合、電流はその経路にすべて流れ込み、抵抗がある経路には一切流れなくなる。この状況を短絡（ショート）と呼ぶ（図２９）。

f:id:quanta087:20171228212646p:plain
図２９　短絡（ショート）

ここで理解しておいて欲しいことをまとめると、
・並列回路では電流は（抵抗値に応じて）それぞれの経路に分配される。ただし、抵抗ゼロの経路がある場合、電流は全てそちらに流れる（ショート）
ということになる。

全ての話がつながる

さあ、これで準備は整った。この記事で説明してきた知識が、ここで全てつながる。ここまでの知識を、改めてまとめ直してみよう。
・SQUIDの一般論
　・ジョセフソン接合に限界を超えた電流を流すと、抵抗が生じる。
　・SQUIDを貫く磁場が強いほど、より小さな外部電流を注入するだけで抵抗が生じる。

・量子ビットとSQUIDの関係
　・SQUIDは、量子ビットが（リングの外側に）つくる磁場の強弱を読み取る。
　・量子ビットが右回り状態であればSQUIDを貫く磁場は弱く、左回り状態であればSQUIDを貫く磁場は強い。

・回路理論
　・並列回路では電流は（抵抗値に応じて）それぞれの経路に分配される。ただし、抵抗ゼロの経路がある場合、電流は全てそちらに流れる（ショート）。


では、最後のピースをはめよう。量子ビットの状態によってSQUIDを貫く磁場の強さが変わり、それに応じてSQUIDの抵抗ゼロを保ったまま入力できる電流が大きくなったり小さくなったりする。では、「中くらい」の電流を、トリガー信号として入力してみよう。何が起こるだろうか。

量子ビットが「右回り」状態だった場合

もし量子ビットが「右回り」状態だったならば、SQUIDを貫く磁場は小さい。ここに「中くらいの電流」が注入されても、ジョセフソン接合の限界まで余裕がある。抵抗はゼロのままだ。このとき、入力されたトリガー信号は全て抵抗ゼロのSQUIDに流れ、わざわざ「抵抗」がある出力線には流れない。当然、出力線に設置された電圧計も反応しない。トリガー信号を入力してからどれくらい待てば判断できるかは事前に調節してあるとすると、「電圧計に反応がない」ことをもって、量子ビットが「右回り」状態だったのだと確定する（図３０）。

f:id:quanta087:20171228222617p:plain
図３０　量子ビットが右回り状態のとき

量子ビットが「左回り」状態だった場合

もし量子ビットが「左回り」状態だったならば、SQUIDを貫く磁場は大きい。ここに「中くらいの電流」が注入されると、ジョセフソン接合が限界に達して抵抗が生じる。このとき、入力されたトリガー信号が、SQUID側と出力線側に分配される（具体的な比率までは分からないが）。すると、「抵抗」に電流が流れて電圧が生じ、電圧計が反応する。これが、量子ビットが「左回り」に確定したという信号だ（図３１）。このイベントはカウンタやパソコンに記録され、実験データが十分に集まれば統計処理などが行われる。例えば、最適動作点でのラビ振動によって「８０％\displaystyle |+>、２０％\displaystyle |->」という状態をつくり、「８０％\displaystyle |右>、２０％\displaystyle |左>」に変換して測定する。実験を繰り返せば、統計的に電圧計が反応する確率は２０％だという結果が得られるだろう*8*9。

f:id:quanta087:20171228225526p:plain
図３１　量子ビットが左回り状態のとき


以上が、超伝導量子ビットのデータ読み出し、ひいては"量子を「観測」する"ということの正体である。

参考文献

[1] K. Kakuyanagi et al., Dephasing of a Superconducting Flux Qubit. Phys. Rev. Lett. 98, 047004 (2007)
[2]SQUIDの動作原理についてはこちらのウェブページを参考にした
[3] 仙場浩一, 超伝導磁束量子ビットの単一回読み出し, NTT技術ジャーナル2004年1月号
[4] S. Saito et al., Parametric Control of a Superconducting Flux Qubit. Phys. Rev. Lett. 96, 107001 (2006)

最後に

なんとか年内に終わらせることができた。実をいうと磁束量子ビットについてもう一回分書きたいことがあるのだが、まだ勉強していないので、文献が理解できるかどうか次第である。
それにしてもここ数年で本当に量子コンピュータはメジャーな存在になった。（おそらく将来高級言語として用いることを想定した）量子シミュレータ言語などもかなり増えてきた。しかし理論やソフトの情報に比べてまだまだハードの情報が少ないと思う。例えば、「量子コンピュータ Advent Calendar 2017」の記事を見ると、ほとんどが理論やアルゴリズム、シミュレータの話である。もちろんこれらの記事のレベルは大変高い。私自身は理論やシミュレータには弱いので本当に勉強になる。だが、これだけの量の記事がありながらハードに全く触れられていないのは、Qiitaがそもそもプログラミングやソフトウェアのためのサイトだということを差し引いても、少し寂しい。みんなは「”観測”、”観測”って、言葉では簡単にいうけど、実際のチップからはどんな信号が出てくるんだ？」とは疑問に思わないのだろうか？ぜひ、ハード関連の情報がもっと増えて欲しい。そうすれば私も勉強する手間が省けるのに。
というわけで、今年の更新はこれで終わる。皆さん、よいお年を。

*1:そもそも超伝導体そのものにも流せる電流の限界値が存在するが、ジョセフソン接合はそれよりも弱い

*2:一度ジョセフソン接合の限界を超えても、流す電流を弱めれば再び超伝導状態に戻る

*3:向きはあまり関係ないので、上から下に貫いていてもよい

*4:電圧降下は片方のジョセフソン接合が限界を迎えた時点で発生するので、もう片方では電流が相殺しているということは気にしなくてもよい。 また、SQUIDに流れる電流が左回りであるならば左側のジョセフソン接合が電圧降下を起こすというだけの話なので、磁場や電流の向きもあまり気にする必要はない。

*5:正直にいうと、この説明には１００％の自信はない。文献中にはっきりとした記述がないからだ。例えば文献[1]には「量子ビットとSQUIDはinductiveに結合している」とある。また、文献[3]には「量子ビットとSQUIDの２つのループの面積比は最適化されています」とある。これらの情報を総合して自分なりに噛み砕いた結果が、本文中の説明である。

*6:元々超伝導リングに意図的に加えている磁場の一部もSQUIDに入っているはずだが、それは考慮に入れなくてもおおざっぱな理解には問題ないだろう

*7:例えば、原図では回路安定化のためのキャパシタやローパスフィルタなどがついているが、全部省いた。また、測定器は原図には書かれておらず、おそらく「抵抗」よりも奥に設置されているはずだが、私の図では説明の都合上「抵抗」にくっつけた。本質的な部分は変わらないはずなので許してほしい。

*8:最初はどのくらいの時間のラビ振動でどのくらいの確率の重ね合わせになるかは分からないので、試行錯誤で調節する。そのあたりの話はシリコン編おまけを参照のこと。

*9:もともと最適動作点の状態でいうと、\displaystyle |->であれば信号が検知されるということである。ところで実は、最適動作点からずらすときに磁場を強めて「左回り」のエネルギーが低いような状態にすると、測定信号が反転する。すなわち、最適動作点で\displaystyle |+>であれば信号が検出されるようになる。この仕組みをうまく説明できるほど理解できなかったので、磁場を弱める方を採用した。信号が反転するということ自体は、文献[4]を参照のこと。

量子コンピュータの基本素子・量子ビットのハードウェア実装（超伝導磁束編その４～データの書き込み・演算～）

量子テクノロジー

その１～素子構造～
その２～超伝導リング詳細～
その２.５～ノイズ耐性～
その３～初期化～


では、初期化した量子ビットにデータを書き込んで演算を行う工程、つまり量子ゲートがいかにハードウェアで実現されるかを見ていこう。基本的には以前紹介したシリコン量子ビットと同じく「磁気共鳴」という物理現象を利用する。磁気共鳴や量子ゲートの基本的な説明は既にシリコン編その４で書いているので、まずはこちらを読んでいただくのがよいと思う。あまり新しく説明することはないが、シリコン量子ビットと超伝導量子ビットの状況の違いには注意していただきたい。


磁束量子ビットの制御回路

まず、量子ビット周りの回路の全体図をもう一度確認しよう（図１５）。

f:id:quanta087:20171222074806p:plain
図１５　回路全体図（文献[1]より転載）

量子ビットの制御に用いるのは図の⑤、「Microwave line」と書かれた導線である。シリコン量子ビットの時と同じく、この導線に交流電流を流すと超伝導リングを垂直に貫く交流磁場が発生する（図１６）。

f:id:quanta087:20171222221924p:plain
図１６　交流磁場

この交流磁場の周波数（≒電磁波のエネルギー）が量子ビットのエネルギー差に一致しているとき、「磁気共鳴（ラビ振動）」という物理現象が起き、量子ビットというハードウェアが駆動されるのだった。なお、実際に使われる電磁波の周波数は大体数ギガヘルツ程度のオーダーである（例えば文献[1]では3.9GHz）。この共鳴周波数はジョセフソン接合の接合面積に依存する。すなわち、工場でチップを製造した時点で決まってしまう。これでは実用上都合が悪いので、後から現場でエネルギーを調節できるように工夫した「ギャップ可変」量子ビットという技術も存在する（文献[2]）。

ところで、シリコン編では磁気共鳴に必要な条件として、「静磁場とそれに垂直な交流磁場」という言い方をした。しかし実は、大事なのは「静磁場と交流磁場」ではなくて「量子ビットと交流磁場」が垂直なことである*1。スピンの場合は静磁場と平行な向きを向いていたが、超伝導量子ビットは最適動作点では静磁場に対して垂直であり*2、交流磁場と静磁場は平行である*3。つまり超伝導量子ビットが回転する様子をブロッホ球で表すと、図１７のように赤道面を回っているように見えるのだ（文献[3]）。

f:id:quanta087:20171223015014p:plain
図１７　超伝導量子ビットの「回転」

このように、量子ビットを適当な時間ラビ振動させることによって「プラスの重ね合わせ状態」と「マイナスの重ね合わせ状態」が重ね合わせ状態になる*4。すなわち、１量子ビットに対する「θ-回転ゲート」という量子論理演算が行われ、超伝導量子ビットに \displaystyle \alpha |+> + \beta |->というデータが書き込まれたことになる。

なお、この記事では「（データの）書き込み」と「演算」をあまり区別していない。普通のコンピュータなら「論理回路で演算した結果をレジスタに書き込む」などといった使い分けをするだろうが、量子ビットは演算結果をそのまま自分で持ち続けるため、論理回路とメモリを兼ねているのだ。将来技術が進んで、演算用量子ビットと保存用量子ビットが明確に分かれたらまた話が違ってくるだろう。

データ読み出し（測定）の準備

量子ビットへのデータ書き込みの話はこれで終わる。ここからの話は本来は次回に回すべき内容である。しかし次回は話が少し難しくなるので量的に分散しておきたいという事情に加え、頭の中に \displaystyle \alpha |+> + \beta |->という式が残っているうちに説明しておきたいことがあるのだ。そういうわけで、少しだけフライングさせてもらいたい。

端的に言おう。実は、最適動作点で動作している量子ビットからはデータを取り出すことができないのだ。その理由はデータ読み出しの仕組みにある。詳細は次回の記事で説明するが、磁束量子ビット読み出し用回路（SQUID）は、右回りの電流が生み出す磁場と、左回りの電流が生み出す磁場の強さの違いを読み取る。しかし最適動作点ではどちらの磁場の強さも同じであるため、この手法を使うことができないのである（文献[1]）。このため、データを読み出すときだけは超伝導リングに与えている磁場を\displaystyle \frac{\Phi_0}{2}から弱め、「右回り」と「左回り」がビットとなるような状態に戻してあげなくてはいけないのだ*5*6。
ここで重要なのは、この操作の前後で、量子ビットに書き込んだデータはきちんと保存されているということだ。すなわち、 \displaystyle \alpha |+> + \beta |->という重ね合わせ状態は、そのまま \displaystyle \alpha |右> + \beta |左>という状態に変換されるのである（図１８）。

f:id:quanta087:20171223194507p:plain
図１８　データ読み出し（測定）の準備

とりあえずこれでデータ読み出しのための準備は完了だ。この続きは次回の記事で説明することにする。

今回はここまで

次回は量子ビットに書き込んだデータを読み出す技術、言い換えれば量子ビットを「観測」する方法について説明する。
その５～データの読み出し（測定）～

参考文献

[1] K. Kakuyanagi et al., Dephasing of a Superconducting Flux Qubit. Phys. Rev. Lett. 98, 047004 (2007)
[2] 齊藤志郎 他, 超伝導量子ビットとスピン集団のコヒーレント結合, NTT技術ジャーナル2012年6月号
[3] Irinel Chiorescu, Non-linear driving and Entanglement of a quantum bit with a quantum readout. （スライド）

*1:なお、ここでいう向きとはブロッホ球上での話であり、必ずしも実空間と対応するとは限らない

*2:専門用語では「システムハミルトニアンが\displaystyle \sigma _xである」という

*3:同じく「場のハミルトニアンが\displaystyle \sigma _zである」という

*4:ややこしいので、「右回り」と「左回り」が基準になるような見方に直して、重ね合わせの「位相」が変化していると考えれば見通しがいいかもしれない

*5:磁場は強めてもよくて、実際強めている文献の方が多いのだが、ここでは説明の都合上弱める方を採用する

*6:実際の論文では、例えば「初めから磁場を少し弱くしておいて、交流磁場に直流オフセットをつけることで演算時だけ最適動作点に移し、演算が終わると自動的に磁場が弱い領域に戻る」といった手法が使われている。このような工夫は面白いが話がややこしくなるので、この記事では採用していない

量子コンピュータの基本素子・量子ビットのハードウェア実装（超伝導磁束編その３～初期化～）

量子テクノロジー

その１～素子構造～
その２～超伝導リング詳細～
その２.５～ノイズ耐性～


では、具体的に超伝導磁束量子ビットの制御方法を見ていこう。まずは、量子ビットにデータを書き込む前に強制的にリセットを行う過程である「初期化」について説明する。
なお、量子ビットにおける初期化の意義など、細かい部分は「シリコン編その３」で述べている。この記事でも以前の記事に言及する部分が多いので、適宜参照して欲しい。

量子ビットの初期化

さて、肝心の初期化法だが…実は、それほど難しくはない。ここまで引っ張っておいて何だが、基本的にただ待つだけでよいのだ。
例えば文献[1]には、

The qubit is initialized to the ground state simply by allowing it to relax.
（訳：量子ビットは単に緩和させることによって基底状態へと初期化される）

とある。ground state とは「エネルギーが最も低い状態」のことであり、relaxとは「エネルギーが高い方から低い方に落ちる」ことを指している*1。
また、文献[2]では、

The qubit is thermally initialized to be in \displaystyle |g> by waiting for 300 \displaystyle \mu s.
（訳：量子ビットは３００マイクロ秒間待つことによって熱的に \displaystyle |g> へと初期化される）

と記述されている。こちらは具体的に待ち時間にまで言及している。言うまでもないが、 \displaystyle |g>というのはground state（エネルギーが低い状態）のことである。thermally（熱的に）というのは、エネルギー緩和という物理現象が熱に関係していることを示している。

ところで、具体的に「エネルギーが最も低い状態」が右回り状態なのか左回り状態なのかは、超伝導リングに加えている磁場によって変わる。この記事では基本的に「最適動作点」での量子ビットを扱うので、エネルギーが低いのは「プラスの重ね合わせ」である。なので、もし初期状態が \displaystyle |->状態ならエネルギーが落ちて \displaystyle |+>状態になるし、初期状態が元々 \displaystyle |+>状態であったのなら、何も起こらない。以上のことから、結局どのように量子ビットを初期化するのかをまとめると、
「３００マイクロ秒ほど待てば自然に \displaystyle |+>状態に落ち着く」
ということになる（図１３）*2。

f:id:quanta087:20171215205902p:plain
図１３　磁束量子ビットの初期化


シリコン量子ビットの初期化には量子ドットのエネルギーを底上げするためのトリガー信号が必要であったが*3、超伝導磁束量子ビットではそのような明示的な操作は必要なく、勝手にエネルギーが落ちるのを待つ。単純な仕組みで初期化可能と言えば聞こえはいいが、これは要するに短時間で勝手にデータが消えるということを意味しており、量子ビットの性質としては決して望ましくない。できれば、普段はデータは消えにくいがリセット信号を入れると即座に初期化されるような能動的な仕組みの方が望ましいと、個人的には思う。

エネルギー緩和

さすがに「待つだけ」では内容が薄すぎるので、この機会にエネルギー緩和についても軽く触れておこう。 前回の記事では位相緩和について言及したので、バランスを取る意味合いもある。
上述したように、量子ビットのエネルギーが高い方から低い方に自然に落ちる現象のことをエネルギー緩和という*4。もっと情報科学寄りの言葉で「ビット反転エラー」と呼ぶこともある。0が1になったり、1が0になったりするエラーのことだ。当然、データの保存や演算の最中にエネルギー緩和が起きると、それはデコヒーレンスの原因となる。ちょうどデータが消えて欲しいタイミングでエネルギー緩和が起これば「初期化」と呼ばれるだけの話である。
位相緩和の原因としては磁場ノイズを挙げたが、エネルギー緩和は熱ノイズが原因となって起こる。要するに、量子ビットの持っているエネルギーが熱になって逃げてしまう、言い換えれば「冷めて」しまうのだ（図１４）。量子力学的な観点からは、二準位系が結晶の熱振動の量子である「フォノン」という形でエネルギーを放出する現象であるとも言える。

f:id:quanta087:20171215212733p:plain
図１４　熱ノイズによるエネルギー緩和

ところで、超伝導磁束量子ビットではどのくらいの時間オーダーでエネルギー緩和が起こるのかというと、実験環境によっても違うので一概には言えないが、数マイクロ秒程度のオーダーのようだ。例えば、文献[3]では１.９マイクロ秒、文献[4]では４マイクロ秒という測定値が報告されている*5。先ほど初期化のためには３００マイクロ秒待つ必要があると述べたが、エネルギー緩和は数マイクロ秒で起きるのに、なぜそんなに待たなければいけないのだろうか。おそらく、エネルギー緩和も量子力学的現象である以上確率的に起こるため、「１００％確実に緩和が起こった。初期化は完了した。」と言い切るためには十分な時間的マージンが必要なのであろう*6。まあ、古い文献値をそこまで気にする必要はない。具体的な数値はともかく、「いわゆるデコヒーレンスには位相緩和とエネルギー緩和の二種類がある」という点を抑えておくだけでも、量子コンピュータを理解する上で役に立つと思う。

まとめ

最後にもう一度肝心な部分をまとめておこう。
量子コンピュータで計算を行う前に、量子ビットのデータを強制的にリセットする「初期化」と呼ばれる工程が必要になる。超伝導量子ビットの場合は具体的な操作は必要なく、待っているだけで量子ビットは自然に低エネルギー状態（ここでは \displaystyle |+>状態）へとエネルギー緩和を起こすので、これをもって初期化とみなす。エネルギー緩和自体の時定数は数マイクロ秒程度であるが、実際の待ち時間はもう少し長く、この記事では３００マイクロ秒という例を挙げた。

今回はここまで

次回は、超伝導量子ビットにデータを書き込み、演算を行う手法について説明する。とは言っても、シリコン量子ビットとほとんど一緒なので、書くことはそれほどないが…。

その４～データの書き込み・演算～

参考文献

[1] I. Chiorescu et al., Coherent Quantum Dynamics of a Superconducting Flux Qubit. Science 299, 1869 (2003)
[2] S. Saito et al., Parametric Control of a Superconducting Flux Qubit. Phys. Rev. Lett. 96, 107001 (2006)
[3] F. Yoshihara et al., Decoherence of Flux Qubits due to 1/f Flux Noise. Phys. Rev. Lett. 97, 167001 (2006)
[4] P. Bertet et al., Dephasing of a Superconducting Qubit Induced by Photon Noise., Phys. Rev. Lett. 95, 257002 (2005) (arXiv版はこちら)

*1:ちなみに、エネルギーが低い方から高い方に上がることはexciteといい、エネルギーが高い状態はexcited state（励起状態）という。

*2:具体的に何秒待つかはそれぞれの実験環境によっても違うだろうから、３００マイクロ秒というのは単なる１例だと思ってほしい。

*3:シリコン量子ビットの記事でも、「量子ドットのエネルギーを引き上げた後は待つだけ」と表現した。しかし、「貯蔵電極」を経由する緩和は普通のエネルギー緩和よりはるかに短時間で起こりやすい。これは緩和しやすい状況を意図的に作り出す強制リセット機能であり、何もせずにただ待っているのとはちょっと違う。

*4:高温環境では低い方から高い方に行くことも含めることが多いが、今回のような低温環境では無視できる。

*5:古い文献ばかり読んでいるが、最近のものは改良されて性能が向上しているかもしれない

*6:そもそもエネルギー緩和時間が数マイクロ秒というのは、量子ビットが高エネルギーを保っている確率が約1/3まで減ってしまうまでの時間である。完全に緩和するまでには短くても数十マイクロ秒かかる。