(57)【要約】【課題】  光格子時計の精度を高める。【解決手段】  本発明の実施形態では、光格子時計１００のための原子を収容しうる空洞領域を囲み、外界とつながる開口が少なくとも２つ設けられているシールド壁を備える輻射シールド１０が提供される。空洞領域に向かうシールド壁の内壁面の幾何形状は、内壁面が鏡面反射性である条件と拡散反射性である条件とにおいて、内壁面が放つ輻射と開口を通じ外界から浸入する輻射とこれら輻射の内壁面での反射成分とが原子２に生じさせるシュタルクシフト量を求め、両条件でのシュタルクシフト量の差分が光格子時計における時計遷移の動作をする原子の位置にわたって所定の値を超えないようにされている。本発明の実施形態では、そのような輻射シールドを備える光格子時計１００や、そのための設計方法も提供される。【選択図】  図２
2019129166.tif 000002
【特許請求の範囲】
【請求項１】
光格子時計のための原子を収容しうる空洞領域を囲み、外界とつながる開口が少なくとも２つ設けられているシールド壁を備えており、
前記空洞領域に向かう前記シールド壁の内壁面の幾何形状は、前記内壁面が鏡面反射性である条件と拡散反射性である条件とにおいて、前記内壁面が放つ輻射と前記開口を通じ外界から浸入する輻射とこれら輻射の前記内壁面での反射成分とが前記原子に生じさせるシュタルクシフト量を求め、両条件でのシュタルクシフト量の差分が光格子時計における時計遷移の動作をする前記原子の位置にわたって所定の値を超えないようにされている、輻射シールド。
【請求項２】
前記所定の値が０．８２×１０−１８である、
請求項１に記載の輻射シールド。
【請求項３】
前記幾何形状が前記直線回りの回転に対し実質的に重なる軸対称となっている、
請求項１に記載の輻射シールド。
【請求項４】
前記幾何形状が、
前記２つの開口の一方を頂点とする円錐面をなす第１のファネル部と、
前記２つの開口の他方を頂点とする円錐面をなす第２のファネル部と、
前記第１のファネル部と前記第２のファネル部とをつなぐ円筒面をなす中間部と
を含んでいるものである、
請求項３に記載の輻射シールド。
【請求項５】
前記第１のファネル部および前記第２のファネル部の前記２つの開口を結ぶ向きの長さが、前記中間部の同じ向きの長さと概ね等しくされている、
請求項４に記載の輻射シールド。
【請求項６】
前記幾何形状が前記直線に沿って反転させて実質的に重なる反転対称になっているものである、
請求項１に記載の輻射シールド。
【請求項７】
前記幾何形状が、
前記２つの開口の間の距離を半径としてもち前記２つの開口の一方を中心とする第１の部分球面と、
前記距離を半径としてもち前記２つの開口の他方を中心とする第２の部分球面と
を含んでいるものである、
請求項１に記載の輻射シールド。
【請求項８】
前記内壁面が０．１以下の反射率を示すものである、
請求項１に記載の輻射シールド。
【請求項９】
前記光格子時計のための原子がＳｒ、Ｙｂ、Ｍｇ、Ｃｄ、Ｈｇからなる群から選択されるいずれかである場合において、前記シールド壁は、前記内壁面が、前記開口を通る３００Ｋの外界からの輻射場が動作時の前記原子に与える鏡面反射黒体輻射シフトと拡散反射黒体輻射シフトとの差分が、１０−１８以下となるようなものである、
請求項１に記載の輻射シールド。
【請求項１０】
請求項１〜請求項９のいずれか１項に記載の輻射シールドを備える光格子時計。
【請求項１１】
光格子時計のための原子を収容しうる空洞領域を囲み、外界とつながる開口が少なくとも２つ設けられているシールド壁を備える輻射シールドの設計方法であって、
前記空洞領域に向かう内壁面が鏡面反射性である条件での前記内壁面が放つ輻射と前記開口を通じ外界から浸入する輻射とこれら輻射の前記内壁面での反射成分とが前記原子に生じさせるシュタルクシフト量である第１の黒体輻射シフト値を、前記空洞領域内で時計遷移分光を行う前記原子の各位置において求める鏡面反射計算ステップと、
前記内壁面が拡散反射性である条件での前記内壁面が放つ輻射と前記開口を通じ外界から浸入する輻射とこれら輻射の前記内壁面での反射成分とが前記原子に生じさせるシュタルクシフト量である第２の黒体輻射シフト値を、前記空洞領域内で時計遷移分光を行う前記原子の各位置において求める拡散反射計算ステップと、
前記第１の黒体輻射シフト値と前記第２の黒体輻射シフト値との差分を各位置において求めるステップと
を含み、
前記シールド壁の前記内壁面の幾何形状が、光格子時計における時計遷移分光を行う前記原子の位置にわたって前記差分が所定の値を超えないよう決定される、
輻射シールドの設計方法。
【請求項１２】
前記鏡面反射計算ステップおよび前記拡散反射計算ステップの少なくともいずれかにおいて光線追跡法によりシュタルクシフト量が算出される、
請求項１１に記載の輻射シールドの設計方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は輻射シールドおよびそれを利用する光格子時計に関する。さらに詳細には、本発明は、原子に作用する黒体輻射の影響を抑制した輻射シールドおよびそれを利用する光格子時計に関する。
【背景技術】
【０００２】
従来、高精度な計時のために各種の原子時計、特に光原子時計（optical atomic clocks）が開発されている。特に光の定在波である光格子を魔法波長の光で実現した光格子時計においては、時計レーザーの波長に比べ十分小さい空間的領域に光格子で原子を閉じ込める手法で、時計遷移の周波数についての不確かさを１０−１８レベルに抑制した動作が実証されている。この不確かさを決定する主要な要因の一つが原子を束縛している光格子による光シフトである。この光シフトについては例えば本発明者の改良により抑制が可能となっている（特許文献１）。不確かさを決定するもう一つの要因が黒体輻射（Blackbody Radiation, ＢＢＲ）、すなわちステファン＝ボルツマン則に従う輻射場に起源を持つ黒体輻射に起因する周波数シフトである。例えば、中性ストロンチウム原子（Sr）を用いる光格子時計では、室温の黒体輻射のために時計周波数が１０−１４でシフトすることが知られている。この黒体輻射の影響を抑制するためには、原子の周囲環境を極低温に冷却することにより温度Ｔの４乗（Ｔ４）に比例する輻射エネルギー自体を弱めることが有効である。原子周囲に冷却された周囲環境を実現するために、時計のための電子遷移（時計遷移）を行わせる原子を内部に配置する冷却可能な遮蔽部材（輻射シールド）が一定の成果を挙げている。例えば輻射シールドを９５Ｋに冷却し内部に生じる黒体輻射を抑制したもので、黒体輻射起因の不確かさが０．９×１０−１８、総計の不確かさが７．２×１０−１８となるようなＳｒ光格子時計が報告されている（非特許文献１）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】国際公開２０１６／１２２００１号公報
【非特許文献】
【０００４】
【非特許文献１】Ichiro Ushijima, Masao Takamoto, Manoj Das, Takuya Ohkubo, and Hidetoshi Katori, "Cryogenic optical lattice clocks," Nature Photonics 9, 185-189 (2015), doi:10.1038/nphoton.2015.5
【非特許文献２】T. Middelmann, C. Lisdat, S. Falke, J. Winfred, F. Riehle, and U. Sterr, "Tackling the blackbody shift in a strontium optical lattice clock," IEEE Transactions on 60, 2550 (2011), doi:10.1109/TIM.2010.2088470
【発明の概要】
【発明が解決しようとする課題】
【０００５】
時計遷移の周波数に現われる不確かさは１０−１８に到達しているものの、それをさらに凌駕する１０−１９レベルまたはそれ以上の精度での時間計測が求められている。そのため、冷却する輻射シールドにさらなる改良の余地が残されている。
【０００６】
光格子時計のための輻射シールドには、導入する原子や、光格子や時計遷移のためのレーザー光を内部空間に届けるための開口が少なくとも２つ開けられている。輻射シールドとなる壁（シールド壁）の外には、典型的には室温に応じた黒体輻射が空間に充満しているため、これら開口を通じ室温の放射が輻射シールド内へ浸入する。室温の放射エネルギーはＴ４に応じた大きさを持つため、開口自体のサイズが可能な限り小さくされていてもその影響は無視できない。シールド壁の内壁面に吸収性の材質を配置すれば浸入する黒体輻射の影響をある程度抑制できるが、現実に採用できる低反射コーティングの吸収は必ずしも完全ではなく、例えば１０％程度の反射率が残存することもある。内壁面で吸収されなかった電磁波は多重反射を生じるにもかかわらず、その寄与はこれまで考慮されていない。厄介なのが、実際の低反射コーティングでは反射には拡散性が伴っており、しかもその拡散反射の特性は多くの場合不明であり、正しい反射特性を設計に反映させることに手間がかかる。
【０００７】
本発明は上記問題の少なくともいくつかを解決することを課題とする。本発明は新規な着想に基づいて設計した輻射シールドやそれを用いる光格子時計を提供することにより、時間計測技術の一層の発展に寄与するものである。
【課題を解決するための手段】
【０００８】
本発明者は、開口を通って外界から侵入する放射（主に赤外域となる電磁波）を考慮した輻射シールドの新たな設計手法を見出し、またその設計手法に基づく輻射シールドが光格子時計の時計遷移の不確かさを一層低減させうることを確認し、本発明を完成させた。
【０００９】
すなわち、本発明のある態様においては、光格子時計のための原子を収容しうる空洞領域を囲み、外界とつながる開口が少なくとも２つ設けられているシールド壁を備えており、  前記空洞領域に向かう前記シールド壁の内壁面の幾何形状は、前記内壁面が鏡面反射性である条件と拡散反射性である条件とにおいて、前記内壁面が放つ輻射と前記開口を通じ外界から浸入する輻射とこれら輻射の前記内壁面での反射成分とが前記原子に生じさせるシュタルクシフト量を求め、両条件でのシュタルクシフト量の差分が光格子時計における時計遷移の動作をする前記原子の位置にわたって所定の値を超えないようにされている、輻射シールドが提供される。
【００１０】
加えて本発明では、上記態様の輻射シールドを備える光格子時計も提供される。さらに本発明では、輻射シールドの設計方法も提供される。
【００１１】
輻射シールドのシールド壁には、光格子時計として動作させるために一般には２つの開口を備えている。輻射シールド内の空洞領域には、シールド壁の内壁面が放つ輻射、開口を通じ外界から浸入する輻射、そしてこれらの輻射が内壁面で反射した成分の輻射が、熱に起因した輻射の成分として存在し、これらが原子の遷移周波数にシュタルクシフト（熱的光シフト）として影響を与える。本発明の上記態様では、内壁面が鏡面反射性である条件と、内壁面が拡散反射性である条件、という２つの条件における原子のシュタルクシフト量をそれぞれ求める。ここで、現実の内壁面を反映してシュタルクシフト量を求めるとしても、内壁面の反射の詳細な特性（反射時にどの程度拡散されたり、散乱されたりするか）がわかっていなくてもかまわず、例えば反射率（エネルギーとして反射される比率）だけがわかっていればよい。シュタルクシフト量を求める計算のための内壁面のモデルには、当該反射率の値をもつ完全拡散反射面と、当該反射率の値をもつ鏡面反射面、という仮想的な反射面を仮定するのである。その仮定で算出した２つのシュタルクシフト量から求まる差分を、時計遷移をする時点の原子の位置にわたって求めてそれを所定の値と比較する。この比較でシュタルクシフト量の差分が所定の値を超えていないようなシールド壁の内壁面の幾何形状が得られれば、その幾何学形状は、必ずしも反射特性の細部が判明していない現実の内壁面の場合であっても十分な性能を持っていることとなる。
【００１２】
本発明の上記態様において内壁面とは、シールド壁の表面のうちシールド壁が囲む空洞領域に向かう範囲のものである。この内壁面は、シールド壁に何らかの層や膜などがコーティングとして設けられていればそのコーティングの表面でありうる。シールド壁は、その材質が任意であり、例えば冷却が容易なように熱伝導率の高い金属（銅など）で作製される。また、本出願の説明において、不確かさを伴う数値を、最終の桁に含まれうる不確かさの幅を括弧書きで直後に記載する場合がある。この表記に従えば、例えば０．２０（１）とは、括弧内に記載した「１」程度の不確かさが、その直前の「０．２０」の値の最終の桁である１／１００の位に見込まれることを意味している。つまり、当該表記は、下限および上限をもつ範囲の形式で０．１９〜０．２１と示せることと同一の意味である。ただし、これらの数値表現において不確かさが明示されていることは、その表現対象となった量が伴う誤差やその量が本質的に統計的であるとの性質が表れていることのみを意図するものであり、対象となった量がその範囲に常に限定されることを意味するものではない。これらの用語および表記を含め、本出願においては、不明確にならない限り本発明の属する技術分野の慣用に従う用語法を採用することがある。例えば、赤外域でありうる熱に起因した輻射のビームを光線と呼んだり、必ずしも理想的な黒体輻射といえない現実の輻射に対しても黒体輻射の用語を用いることがある。
【発明の効果】
【００１３】
本発明のある態様では輻射による原子への影響を極力排除して不確かさを１０−１９程度またはそれ以上にまで抑制しうるような輻射シールドや、その輻射シールドを採用する光格子時計が提供される。
【図面の簡単な説明】
【００１４】
【図１】図１は本実施形態の計算にて考慮する要素を示す説明図であり、図１Ａは２つの開口をもつ輻射シールドのモデル、図１Ｂはモンテカルロ法にて計算するパストレースモデル、図１Ｃは拡散反射および鏡面反射の様子を示す模式図をそれぞれ示す。
【図２】図２は、輻射シールドを用いる光格子時計のための概略構成のための説明図である。
【図３】図３は、従来の光格子時計の光格子および原子と輻射シールドの配置および構成を示す構成図である。
【図４】図４は、従来の輻射シールドについて、熱的光シフトの位置依存性を計算した条件およびその結果を示す説明図である。図４Ａは例示の従来の輻射シールドを描画したものである。図４Ｂは、低反射コーティングを施し拡散反射および鏡面反射の二つの反射特性の内壁面を仮定して得られた熱的光シフトの計算結果を示すグラフである。図４Ｃは分光時の原子の想定される位置（z＝１０ｍｍ）からの熱的光シフトの差分を示すグラフである。
【図５】図５は、従来の輻射シールドにおいて計算により求めた熱的光シフトの反射率への依存性を示すグラフである。
【図６】図６は種々の輻射シールドの形状と輻射シールドの内部における熱的光シフトの位置依存性との関係を示す説明図である。
【図７】図７は、本実施形態における輻射シールドの設計方法のフローチャートである。
【図８】図８は、熱的光シフトの位置依存性を、種々の幾何形状で設計された輻射シールドについて、形状のメッシュモデルとともに示すグラフである。
【図９】図９は、熱的光シフトの位置依存性を、種々の幾何形状で設計された輻射シールドについて、形状のメッシュモデルとともに示すグラフである。
【図１０】図１０は、図１０Ａ〜Ｃの順に、図８Ｂ、８Ｃ、８Ｄを拡大したグラフである。
【発明を実施するための形態】
【００１５】
以下図面を参照し、本発明に係る輻射シールドと光格子時計の実施形態を説明する。全図を通じ当該説明に際し特に言及がない限り、共通する部分または要素には共通する参照符号が付される。また、図中、各実施形態の要素のそれぞれは、必ずしも互いの縮尺比を保って示されてはいない。
【００１６】
１．輻射シールドの設計
従来のＳｒ光格子時計において採用した輻射シールドの場合に概算値として見積もられた黒体輻射に起因する周波数シフト（ＢＢＲシフト）の不確かさは０．９×１０−１８であった（上述、非特許文献１）。この概算値は、開口が張る立体角と環境温度の測定とによって見積もられたものである。より精密な計算のためには、多重反射、および環境からシールドに浸入しその後にシールド壁の内壁面にて生じる散乱の両方も加えた上でのシフトを注意深く考慮しなくてはならない。非特許文献１におけるＢＢＲシフトの見積は、輻射シールドが単純な球面モデルである場合に基づいている（非特許文献２）。本実施形態では、多重反射の効果をより精密に扱うために、光線追跡法を採用してシールドの内壁面の幾何形状、反射率、および反射特性の影響も考慮することとする。
【００１７】
具体的には、本実施形態にて、輻射シールドの内部で熱輻射によって誘起されるＢＢＲシフトの光線追跡計算が説明される。当該計算モデルは、非特許文献１で採用した輻射シールドの形状設計に適用される。１０−１９のＢＢＲ不確かさをもつ光格子時計を実現するための輻射シールドの最適な幾何形状を説明する。光格子による光シフトの評価法の改良とも相まって、総計の不確かさが１０−１９となる極低温動作の光格子時計が実現することとなる。
【００１８】
１−１．計算モデル
図１は本実施形態の計算にて考慮する要素を示す説明図であり、図１Ａは２つの開口をもつ輻射シールドのモデルであり、光線追跡アルゴリズムにおいて原子に照射される輻射の３成分を示す。内壁面が放射する放射輝度Ｌｅ、入射する放射輝度Ｌｉ、および出射する放射輝度Ｌｏである。図１Ｂは、モンテカルロ法にて計算するパストレースモデルである。図１Ｃは、完全拡散反射および鏡面反射の様子を示す模式図である。
【００１９】
図１Ａに示すように、開口のあるシールド内の原子は、開口を通って漏れてくる光線と放射率ε（ν）＝１−Ｒ（ν）にてシールドの内壁面から放射される光線とに加え、Ｒ（ν）の反射率をもつ内壁面による熱輻射の反射光線および散乱光線を考慮する必要がある。ここでνは輻射の周波数を表す。そのような多重反射による熱的光シフトによる間接的な成分が、時計の不確かさを１０−１９前半のために重要な寄与をもつ要因である。
【００２０】
ＢＢＲシフトは、原子に照射される熱輻射のエネルギーを受けて生じるシュタルクシフトが、時計遷移の基底・励起それぞれの状態で差があるために生じる周波数シフトである。温度Ｔの熱源から生じる黒体輻射による電界Ｅ（ν，Ｔ）は、黒体輻射の周波数νにおける単位周波数範囲ｄνあたり、次のように定義することができる：
【数１】
2019129166.tif 000003
ここでｈはプランク定数、ε０が真空の誘電率、ｋＢはボルツマン定数である。これを用いてＢＢＲシフトは
【数２】
2019129166.tif 000004
のように書くことができる。ここでΔα（ν）は、時計遷移の基底・励起状態の分極率の差である。ここでは、磁気双極子および高次多重極の寄与を無視している。
【００２１】
しかしながら、ターゲットとなる原子の位置ｒａにおいて熱輻射の分光放射照度を式（１）により表現することは、温度の空間的不均一性のために不正確となる。ｍｍスケールのシールドが熱輻射の波長（〜μｍ）を大きく超えることから、真空において、直線（光線）に沿って伝播する分光放射輝度を考える。空間不均一性を説明する熱輻射のシュタルクシフト（つまり熱的光シフトνth(ｒa)）を正しく見積もるべく、原子を取り囲む球面にわたり放射による光線を積分する目的で、光線追跡法を採用する。
【００２２】
図１Ａに示すような開口を持つシールド内の原子にとっては、開口を通って漏れてきて入射する光線とε（ν）＝１−Ｒ（ν）の放射率にて輻射シールドの内壁面から放たれる光線とに加えて、反射率Ｒ（ν）をもつ内壁面により反射と散乱した熱輻射の光線を考えなくてはならない。多重反射による熱的光シフトのそのような間接的な成分が、時計の不確かさを１０−１９前半にまで減少させるために重要な寄与をもつファクターとなる。
【００２３】
１−２．モンテカルロ光線追跡計算法
モンテカルロ光線追跡計算法は、空間での光の伝播を解くための確率論的なアプローチである。実際の環境において光は、光源からから放たれ、表面で反射され、そしてターゲットとなる原子に到達する。しかし、光源から放たれる殆どの光線は、多重反射の後に原子にヒットすることはない。このため、計算効率の観点から、光線が、ターゲットとなる原子の位置から光源に向かって追跡されるようなバックワードアプローチを光線追跡アルゴリズムとして採用する。
【００２４】
本発明者は、光の反射と散乱の効果の計算のためにバックワードモンテカルロ光線追跡を適用するパストレーシングアルゴリズムを使用する。そのアルゴリズムでは任意の形状、材質、および任意の周波数の光を扱うことができる。本発明者は、pbrt- a physically based ray tracerと呼ばれるオープンソースの光線追跡ソフトを実際のパストレースのために採用する(http://www.pbrt.org/index.html)。モンテカルロ光線追跡の包括的理論およびその完全な実装の詳細が成書（Matt Pharr, Wenzel Jakob, and Greg Humphreys, "Physically Based Rendering, Third Edition: From Theory to Implementation," Morgan Kaufmann (2016), ISBN-13: 978-0128006450）に述べられている。
【００２５】
パストレーシングアルゴリズムは大略次のとおりである。図１Ｂは、当該アルゴリズムの説明のために示す。まず、ターゲット原子からランダムに選んだ方向に向かい１パス当たり１つの光線が追跡される。光線が表面にヒットした時点で、反射率に基づいて当該光線が吸収されるか反射されるかが内壁の反射率に従い確率的に決定される。光線が反射される場合、当該表面の双方向反射率分布関数（bidirectional reflectance distribution function, ＢＲＤＦ）に基づいて新たな方向が確率的に選択される。光線は、表面で完全に吸収されるか、輻射シールドの外に到達するまで追跡される。
【００２６】
レンダリング方程式は、光線追跡アルゴリズムにて光の伝播を記述する積分方程式である。この積分方程式は、光線がヒットする表面の全ての位置で値が求められるものであり、次式のように記述される。
【数３】
2019129166.tif 000005
ここで、Ｌｅは放出されるスペクトル放射輝度、Ｌｉは到来するスペクトル放射輝度、Ｌｏは出射するスペクトル放射輝度、ｆはＢＲＤＦ、Ωｒは位置ｒにおける半球、ｅｉは入射光の逆を向く方向、ｅｏは出射光の方向、ｎは反射する表面に垂直な方向の単位ベクトルである。
【００２７】
このパストレースアルゴリズムは、任意のＢＲＤＦを扱いうるものの、本発明者はｐｂｒｔの実装ではいくつかのＢＲＤＦについてエネルギー保存則を確認できなかった。そこで本発明者は、本願の輻射シールドの反射モデルのために、エネルギー保存則を確認した完全拡散反射モデルおよび完全鏡面反射モデルを使用する。
【００２８】
図１Ｃは、これら２つのモデルを図示している。完全拡散反射モデルは、反射面のためにランバート則に従う面を仮定するものである。光は、入射方向に無関係に半球方向に一様に反射される。他方の完全鏡面反射モデルでは、反射角が入射角と等しくなる。
【００２９】
熱的光シフトνth(ｒa)  が入射スベクトル放射輝度の積分によって計算できる：
【数４】
2019129166.tif 000006
本願において２つの別々の光源（熱源）を仮定する：一つは輻射シールド外部の室温環境で温度Ｔextのものであり、もう一つは輻射シールドの内壁面で温度Ｔwallのものである。
【００３０】
本発明者は、光源からの放射輝度を、
【数５】
2019129166.tif 000007
のように定義する。本願ではε（ν, Ｔ）はε（ν, Ｔext）＝１、また、ε（ν，Ｔwall）＝ε（ν）を表している。後者ε（ν）は、輻射シールドの内壁面の放射率である。熱的光シフトは、ＴextとＴwallの２つのスペクトル放射輝度積分値を合算することによって計算される。
【００３１】
本発明者は、全ての方向について光線が一様にサンプルされる全方向カメラモデルを光線追跡計算に導入した。このカメラをターゲット原子の位置に置き、モデルの輻射シールドにより囲まれるシーンを描画する（図示しない）。この輻射シールドの内壁面は、適当な数の画素（例えば球面を２００×１００の正方形を２つ並べた画像）に対応させることができる。出射スペクトル放射輝度Ｌｏ（式（３））の値が各画素値として対応する。その画素値によって表される反射した放射輝度の和をとることによって式４の積分方程式を計算する。
【００３２】
１−３．輻射シールドのための光線追跡計算
１−３−１．  低温動作光格子時計
図２は、輻射シールドを用いる光格子時計のための概略構成のための説明図である。また図３は、従来の光格子時計の光格子および原子と輻射シールドの配置および構成を示す構成図である。低温動作光格子時計１００では、Ｓｒ原子を極低温環境下のもとで時計遷移分光を行うために、分光時の原子周辺の領域を輻射シールド１０により囲むことにより実現される。Ｓｒ原子等の原子２は、原子供給部から供給されて輻射シールドの外部で三次元的に対向する冷却レーザー４の光線で冷却・捕獲される。その後、光格子レーザー６のつくる光格子にトラップされる。光格子は、輻射シールド１０を貫くように例えば一次元に配列する対向する同じ波長（例えば魔法波長）の光によって作られる定在波であり、波長を変調できる変調素子６２によって片方の光の波長を変調することにより、トラップした原子を輻射シールド１０内部に精密に搬送することができる（移動光格子）。原子２は、輻射シールド内部に例えば１０ｍｍだけ搬送され、光格子レーザー６による光格子にトラップされたまま輻射シールド１０内の適切な位置で時計レーザー８等の光学的手段によって時計遷移が観察される。なおこれらの光格子時計の動作は、適切な範囲をもつ真空槽２０内で実行される。ここに説明したもの以外も含め、光格子の動作の詳細は、非特許文献１において説明される。中性のＳｒ原子を用いる光格子時計では、８７Ｓｒ原子のレーザー冷却およびトラップが行われ、当該原子の一次元光格子への投入が行われ、移動光格子によって、輻射シールドの内部１０ｍｍに搬送される。図３においては、輻射シールド内の原子の位置を、輻射シールド１０の中心軸であり光軸でもあり、輻射シールド１０の一つの開口を原点にもつｚ軸に対して、ｍｍの単位で示している。光格子は、１０ｍｍのレイリー長、ビームウエストが５０μｍ半径程度である。輻射シールドの幾何形状や位置は、トラップや光格子のためのレーザービームとの干渉を避けるよう設計される。輻射シールドは２つの開口を持っており、従来のＳｒ光格子時計では、原子とレーザーが内部へ達しうるよう、それぞれ例えばφ１＝０．５ｍｍ、およびφ２＝１．０ｍｍとされて、両開口の原子に対する立体角はΩＡＰ＝８．９ｍｓｒである。
【００３３】
輻射シールドは冷却され、例えばスターリング冷凍機によって９５Ｋにされる。この輻射シールドの内壁面には、室温（例えばＴext＝２９６（５）Ｋ程度）の放射が多重反射するのを防止する目的で低反射コーティングがなされている。低反射コーティングの半球反射率Ｒ（ν）は、ν＝２７〜１５０ＴＨｚで約０．０２であり、これはν＝１０〜３００ＴＨｚで０．１未満であるため、ν＜１０ＴＨｚ、ν＞３００ＴＨｚについて０．１（１）の反射率を想定する。ここまでの一般的な計算手法の最初の解析例として、最初に、従来の輻射シールドで得られる結果を説明する。
【００３４】
１−３−２．熱的光シフトの位置依存性
最初に従来の輻射シールド内での熱的光シフトの位置依存性を計算する。図４は、従来の輻射シールドについて熱的光シフトの位置依存性を計算した条件およびその結果を示す説明図である。図４Ａは、例示の従来の輻射シールドを模し計算に用いた輻射シールドをメッシュで描画したものである。図４Ｂは、低反射コーティングを施し拡散反射および鏡面反射の二つの反射特性の内壁面を仮定して得られた熱的光シフトの計算結果を示すグラフである。図４Ｃは、ｚ＝１０ｍｍ付近の位置における熱的光シフトのｚ＝１０ｍｍからの差分を示すグラフである。
【００３５】
輻射シールドは図４Ａに示すように６３５０個の三角メッシュにてモデル化される。輻射シールドのモデルにおいて角ができるコーナー部分に２ｍｍ半径の丸み（隅肉半径（fillet radius））を与えることにより、エネルギー保存則が、拡散反射および鏡面反射それぞれの熱的光シフトにおいて、２×１０−２１および４×１０−２０の範囲で満たされた。
【００３６】
輻射シールド内壁面の各部の放射輝度Ｌｏ（式（３））を表す２００×１００画素の画像において、各画素について、光線追跡のために原子から逆方向に６４光線が追跡され、１本の光線の追跡において最大１０００回の反射がシミュレーションされた。隅肉半径ならびにメッシュ数および光線数を変更し、その値のバラつきから、拡散反射および鏡面反射に対する計算でのシフトの不確かさは、４×１０−２０および２×１０−１９と見積もることができる。位置依存性は、図３の輻射シールドにおける光線の始点をｚ軸で変更することにより評価される。
【００３７】
図４Ｂは、分数周波数単位で熱的光シフトνth(ｒa)  ／νＳｒの計算結果を示している。ここでＳｒの時計周波数およびＤＣ分極率差を、νＳｒ〜４２９ＴＨｚおよびΔα（０）＝４．０７８７３（１１）×１０−３９ｃｍ２／Ｖとした。黒四角および白円は、それぞれ、拡散反射および鏡面反射の場合の熱的光シフトを示している。拡散反射においては内壁面の形状に起因する位置依存性は内壁面による一様な反射のために消失している（後述）。
【００３８】
しかしながら、鏡面反射の場合においては、内壁面の形状がシフトに位置依存性をもたらしている。内壁面の角度のために、熱輻射がｚ＜９ｍｍの領域に集中しており、これが熱的光シフトの増強をもたらしている。つまりシフトの位置依存性は鏡面反射の場合に顕著である。
【００３９】
位置依存性は実験的にも計測される。輻射シールドの温度を９５Ｋに固定しながら原子の位置を変化しつつ周波数シフトを測定する。図４Ｃにおいて、白抜き三角が実験的に測定された熱的光シフトを示しており、採用した低反射コーティングにおいては、鏡面反射による寄与が支配的ではないことを示唆している。縦方向のエラーバーは、計測における統計的不確かさ示している。原子の位置の不確かさは１ｍｍである。しかしながら相対的な位置の不確かさは、移動光格子の周波数離調により決まるものであることから、無視することができる。室温の黒体輻射における波長全域について、表面のＢＲＤＦを特徴付けることは容易ではないため、熱的光シフトの位置依存性と、拡散反射および鏡面反射の間における周波数シフト差とを低減するような輻射シールドの設計を検討する。
【００４０】
１−３−３．熱的光シフトの反射率依存性
図５は、従来の輻射シールドにおいて、その中心位置（ｚ＝１０ｍｍ）で拡散反射および鏡面反射の場合について計算により求めた熱的光シフトの反射率への依存性を示すグラフである。黒四角および白円は、それぞれ拡散反射および鏡面反射についての光シフトの反射率依存性を示している。反射率Ｒ＝０についてのシフトが、室温放射の反射による間接的な寄与がない場合のシフト（図５における鎖線）に対応している。
【００４１】
間接寄与のシフトは、多重反射のために反射率とともに増加する。その増え方は、鏡面反射について拡散反射のものよりも顕著である。反射率の０．１の増分によって、拡散反射および鏡面反射の成分の間に３×１０−１８の周波数差がもたらされる。すなわち、原子の位置において、拡散反射と鏡面反射の間における周波数差が小さくなるように設計できるなら、そのような設計は有用なものとなるであろう。
【００４２】
１−３−４．熱輻射によるシフトの総計の不確かさ
温度Ｔextの熱輻射についての熱的光シフトは、輻射シールドの外界から原子に直接照射する放射輝度に起因する直接項νextｄ（ｒａ）と、内壁面によって反射した放射輝度に起因する間接項νextｉ（ｒａ）との２つの寄与に分類される。間接シフトは式（３）から求まる総シフトから、容易に計算が可能な直接シフトを差し引くことによって計算される。
【００４３】
低反射コーティングを持つ輻射シールドについて間接熱輻射シフトは、不確かさ６×１０−２０を伴う−０．２８×１０−１８と求められた。光線追跡を利用する間接熱的光シフトの正確な計算により、総計の熱的光シフトの不確かさが０．８２×１０−１８と求められる。
【００４４】
１−４．輻射シールドの形状設計における最適化例
以上の一連の説明では、従来採用されている輻射シールドを対象として、その性能を再検討したものであった。本実施形態では、輻射シールドの幾何形状を設計し最適化するために、光線追跡を用いた上記手法を採用する。特に、輻射シールドの内壁面の反射特性の細部が不明なままであっても、従来の輻射シールドで求まった上記０．８２×１０−１８の値を凌駕する輻射シールドが設計できることについて説明する。
【００４５】
１−４−１．チューブ状、球面、共焦点形状の輻射シールド
本節では、最初に、幾何形状のもたらす影響を調査するために、仮想的な輻射シールドの形状に対するシフトの依存性について説明する。構造体内部での反射の基本的な性質を理解するために、本発明者は、２つの開口を持つ、チューブ状、球面、および共焦点の輻射シールドという単純な幾何形状を取り上げる。内壁面について、輻射シールドの反射率Ｒ＝０．１と、Ｔwall＝９５Ｋの極低温と、輻射シールド外のＴext＝２９６Ｋの室温とを仮定する。両端にφ＝１ｍｍの２つの開口を持ち、長さ２０ｍｍの輻射シールドを想定する。
【００４６】
図６は種々の輻射シールドの形状と輻射シールドの内部における熱的光シフトの位置依存性との関係を示す説明図である。図６Ａは、チューブ状の設計に対し計算したシフトの位置依存性を示す。熱的光シフトは、拡散反射（黒四角）についてｚ＝１０ｍｍの位置で−５．７×１０−１７である。しかし、５．８×１０−１８の周波数差が拡散反射および鏡面反射（白円）の間に生じる。図６Ｂは、球面の設計に対し計算したシフトの位置依存性を示す。ｚ＝８ｍｍ〜１２ｍｍの範囲の位置で、拡散反射および鏡面反射（白円）の間で周波数差が、１．７×１０−１８以下となる。鏡面反射では、中心から半径を二分した距離だけ離れたところに近い２つの位置でディップが生じた。これらのディップは、球面の内壁面からの反射による熱輻射の集積量の増加に起因するものである。これは、２つの開口の位置にある点光源の場合について予想されるものである。
【００４７】
熱輻射が輻射シールド中心の近くに集中するのを防ぐため、２つの球面が焦点を共有して持ち、２つの球面のそれぞれの中心が開口に位置するような、共焦点の設計を想定する。計算結果が図６Ｃに示されている。球面の輻射シールドの鏡面反射の場合に存在したディップは消失し、拡散反射と鏡面反射のシフトは、ｚ＝１０ｍｍの±４ｍｍの範囲にわたって６．５×１０−１９の不確かさの範囲で互いに一致した。
【００４８】
１−４−２．現実の反射特性の下での輻射シールドの最適化指針
本実施形態における最適化指針は、輻射シールドの開口を結ぶ線上のうち原子の時計遷移分光を行う可能性のある各位置（原子をレーザーで捕獲する際の初期位置の不確かさも含む）に対し、熱的光シフトの位置依存性を拡散反射と鏡面反射とで計算し、それらの間の差をできるだけ低減するような輻射シールドを採用する、というものである。上述したように、現実の低反射のコーティングは複雑な反射特性を示し、反射特性を例えばＢＲＤＦ等で輻射の全波長域で特徴付けることは容易ではない。他方、ランバート則に従う完全拡散反射面の仮定で、および鏡面反射のみの完全鏡面反射の仮定で計算を行うことは、図４、図５の結果からわかるようにいずれも十分に現実的である。完全拡散反射と完全鏡面反射は、現実の反射現象を理想化した対照的な２つの条件といえるため、仮にそれらの条件の間における性能に差がない形状設計が可能であれば、その形状設計を採用し内壁面に現実の低反射コーティングを施した輻射シールドに期待しうる熱的光シフトの値にも不確定さは残らないはずである。反射率程度だけが知られている拡散性が特徴付けできないような低反射コーティングがあったとしても、拡散反射のみおよび鏡面反射のみという理想化した２つの場合で差がないのであれば、現実の低反射コーティングを施したものもその差の少ない範囲にあることが期待できる。この指針での最適化された形状設計で輻射シールドを作製すれば、現実の低反射コーティングなどの表面に位置する材質について、必ずしも明確とはいえない拡散反射特性を持っていても、予測された十分な特性が実現されるのである。
【００４９】
１−４−３．軸対称設計での輻射シールドの最適化
本発明者は、図３に示したレーザービームの配置と両立可能な幾何配置にまで上述した計算を拡張し、その上で、上記指針に従って最適化を行う。図７は、本実施形態における輻射シールドの設計方法のフローチャートである。黒体輻射場の反射の影響を弱めるために低反射コーティングを施した輻射シールド１０（図２）を設計するとしても、低反射コーティングを含めた内壁面の反射特性が、ＢＲＤＦのような詳細なレベルで判明していることは希である。本実施形態では、そのような一般的な状況においても、高性能な輻射シールドを設計しうる手法を採用する。その設計では、図７に示すように、まず初期設計として、内壁面が何らかの幾何形状を設定する（Ｓ１０２）。この設定は、光線追跡の計算を行う準備であるため、コンピュータ内にて形状を定義しうる任意の手法で行われる。次に、例えば内壁面の低反射コーティングとして実際に用いる材質の反射率Ｒを用いて、鏡面反射での黒体輻射シフトの値を計算する（Ｓ１０４）。この値を第１黒体輻射シフト値と呼ぶ。同様に、例えば上記反射率Ｒを用いて、拡散反射での黒体輻射シフトの値を計算する（Ｓ１０６）。こちらは第２黒体輻射シフト値と呼ぶ。黒体輻射シフト値は、いずれも輻射シールド１０内の空洞領域のどこでも計算は可能であるが、開口をつなぐ直線上など原子が時計遷移分光の際に存在しうる位置が輻射シールドの直接的な性能指標となるのでそのような位置で計算する。そして、各位置で第１黒体輻射シフト値と第２黒体輻射シフト値の差分を算出する（Ｓ１０８）。この差分の値が所定の値を超える場合には（Ｓ１１０、Ｙｅｓ）、再度、内壁面の幾何形状を設定する。差分の値が所定の値を超えなければ（Ｓ１１０、Ｎｏ）、その時点の幾何形状を採用する（Ｓ１１２）。ここでの所定の値の例として例えば０．８２×１０−１８を採用すれば、設計された幾何形状は、従来の輻射シールド（図４）よりも改善されたものといえる。なお、上記差分を求めるために想定するべき原子の位置範囲やこの所定の値は、動作条件や設計の基準であるため任意に設定することができる。また、図示しないが、追加の条件を課して同様の設計を行うことも有用である。例えば、冷却レーザー４（図２）のためのレーザーの配置を干渉しないように設計することがこの追加条件の例の一つである。さらに、所定の値として、固定された数値の比較ではなく、直前の設計における上記差分を、次の設計では所定の値として、性能の改善を継続的に行うこともできる。
【００５０】
現実の輻射シールドを採用する光格子時計において、レーザー冷却された原子を光格子によりトラップしたまま、輻射シールド内に搬送する必要がある。この間において原子は光格子をなすレーザーに適切にトラップされ続けなくてはならない。したがって、輻射シールドの一方の開口は、レーザー冷却のためのレーザーの照射が可能なようにファネル形状になってある程度突出しているものが好ましい。また、レーザー冷却の位置から輻射シールド内で時計遷移分光をする原子の位置までの距離（搬送距離）は、概ねレイリー長の２倍の２０ｍｍ以内に留めたい。原子を光格子にトラップするために光格子レーザーはある点でフォーカスさせてその光強度を上げているが、搬送距離を長くしフォーカスから離れると、光格子のレーザーは拡散し、トラップに十分な光強度を持たせにくくなるためである。本実施形態では、上述した計算手法をこれらの現実的な制約を持つ実の光格子時計のための輻射シールドの設計に適用する。その際輻射シールドの内壁面に反射率Ｒ＝０．１、２つの開口間距離Ｌ＝２０ｍｍ、温度９５Ｋが仮定されている。また、ファネル形状となっている部分（ファネル部）は、互いに円錐の開く向きが反対となっていて、二つの開口を結ぶ向きにそれぞれ長さＬ２をもつ。Ｌ１は、二つのファネル部をつなぐ円筒面をなす中間部の、同じ向きの長さである。本実施形態の最適化では、Ｌ＝Ｌ１＋２Ｌ２の関係を満たしながら、Ｌ１およびＬ２を変更する。黒四角および白円は、それぞれ拡散反射と鏡面反射の熱輻射シフトを示している。
【００５１】
図８および図９は、熱的光シフトの位置依存性を、種々の幾何形状で設計された輻射シールドについて、形状のメッシュモデルとともに示すグラフである。図８、９における設計の殆どにおいて、球面形状の輻射シールドについて説明したのと同様に、鏡面反射において２つのディップが生じる。拡散反射および鏡面反射の間でのシフトの差は、図８Ａから図８ＣにかけてＬ１を増加するのに応じて減少し、輻射シールドの中心付近の間接的熱輻射の集積のために、図８Ｄから図９Ｂにかけて増加し始める。Ｌ１がさらに増加すると、鏡面反射の場合に中心にディップが現われはじめ（図９Ｃ）、曲線形状はチューブ状の輻射シールドのものに類似し始める（図９Ｄ）。
【００５２】
図１０は、図１０Ａ〜Ｃの順に、図８Ｂ、８Ｃ、８Ｄの縦軸をｚ＝１０ｍｍの拡散反射のシフト値からの差に変更し、横軸および縦軸とも拡大したものである。シフトの位置依存性は、輻射シールドの中心近傍（ｚ＝９〜１１ｍｍ）にて２×１０−１９まで減少する（図１０Ｂ）。つまり、拡散反射および鏡面反射のシフトの差分は、この範囲にて２×１０−１９まで減少する。図１０Ａ〜１０Ｃ（図８Ｂ〜８Ｄ）において、ｚの範囲を広げた際に熱的光シフトの観点で最も有利なものは、図１０Ｂ（図８Ｃ）のＬ１＝７ｍｍ、Ｌ２＝６．５ｍｍの条件のものである。ただし、図１０Ａ〜１０Ｃを見る限りいずれも１０−１９前半というわずかな不確定さにとどまっており、Ｌ１＝７ｍｍ、Ｌ２＝６．５ｍｍで輻射シールドを作製する際に１ｍｍ程の許容範囲があることが示唆される。
【００５３】
１−４−４．他の原子への適用
ここまでの計算はＳｒ原子を対象とした計算であったが、光格子時計のための他の候補であるＹｂ、Ｍｇ、Ｃｄ、およびＨｇといった他の原子種に容易に拡張できる。プランクの法則による分光放射輝度の式および式（４）から実効立体角を次式により定義する。
【数６】
2019129166.tif 000008
νstat（静的項）をＴ0＝３００Ｋで時計遷移の基底と励起の２つの状態間での静的分極率の差から生じるＢＢＲシフトの寄与、νdyn（動的項）をＢＢＲスペクトルの周波数範囲に応じた分極率差の変動による温度Ｔ0でのＢＢＲシフトの寄与とする。補正のためにνstatおよびνdynを用いることによって、
熱的光シフトは、次式のように表現される：
【数７】
2019129166.tif 000009
ここで、ａ＝νstatΩeff (Ｔext)／４π、ｄ＝νdynΩeff (Ｔext)／４π、ｃ＝νstatΩeff(Ｔwall)／４π、およびｄ＝νdynΩeff (Ｔwall)／４πであり、図7Ｃの形状設計において、Ωeff (Ｔext）＝１８．４ｍｓｒ、Ωeff (Ｔwall）  ＝１２．５４８ｓｒである。
【００５４】
ＳｒおよびＹｂ原子について、νstatおよびνdynが実験と理論で調査されている。Ｍｇ、Ｃｄ、Ｈｇ原子についてνstatが理論で調査されている。表３は、温度Ｔ１、Ｔ２において各原子の寄与をまとめたものである。熱的光シフトは、これらの係数を利用して計算される。表１は輻射シールドについて図７Ｃの幾何設計を採用して光線追跡法を使用して計算された原子種についての熱的光シフト係数を示す。
【表１】
2019129166.tif 000010
【００５５】
２．まとめ
本願発明者は、光線追跡モデルに基づいて、極低温にした輻射シールドをもつ光格子時計について、熱輻射により誘起される光シフトを計算した。多数の輻射シールドの構造に対してシフトを計算し、熱輻射の多重反射によるＢＢＲシフトの空間的不均一性のために、２つのタイプの反射（完全拡散および鏡面）において１０−１７レベルの周波数差があることを見出した。
【００５６】
２つのタイプの異なる反射において、原子位置の不確かさが±１ｍｍのときにシフトの不確かさが２×１０−１９にまで削減されるような最適で実用的な輻射シールドを設計した。このデザインは、拡散反射および鏡面反射の間においてシフトの周波数差を２×１０−１９未満に減少させるものである。この設計は、他の原子種を用いる極低温の光格子時計に容易に適用可能である。この新規の設計は、１０−１９レベルの原子時計の可能性を拓くものといえる。
【００５７】
以上、本発明の実施形態を具体的に説明した。上述の実施形態、変形例および実施例は、本出願において開示される発明を説明するために記載されたものであり、本出願の発明の範囲は、特許請求の範囲の記載に基づき定められるべきものである。実施形態の他の組合せを含む本発明の範囲内に存在する変形例もまた、特許請求の範囲に含まれるものである。
【産業上の利用可能性】
【００５８】
本発明は、計時のための光格子時計の輻射シールドとして、またそのような輻射シールドを備える任意の原子時計のために使用可能である。
【符号の説明】
【００５９】
１００  光格子時計
２  原子  
４  原子冷却用レーザー
６  光格子レーザー
６２  変調素子
８  時計遷移分光用レーザー
１０  輻射シールド
２０  真空槽
【図１】
2019129166.tif 000011
【図２】
2019129166.tif 000012
【図３】
2019129166.tif 000013
【図４】
2019129166.tif 000014
【図５】
2019129166.tif 000015
【図６】
2019129166.tif 000016
【図７】
2019129166.tif 000017
【図８】
2019129166.tif 000018
【図９】
2019129166.tif 000019
【図１０】
2019129166.tif 000020