新学術領域研究「重い電子の形成と秩序化」のニュースレターに最近の研究内容の解説文が掲載されました。PDFでご覧いただけます。（柳澤の記事はpp.25より）

新学術領域研究「重い電子の形成と秩序化」のニュースレターに最近の研究内容の解説文が掲載されました。PDFでご覧にいただけます。（柳澤の記事はpp.12より）

第１章　はじめに

わたしたちは物性を調べるとき，物質に様々な外場をかけ，それに対する物質の持つ様々な自由度の応答を観測する．例えば磁性を調べる場合，最も簡単な方法は物質に磁場を加え，物質中のスピン自由度の応答をその周りに巻いたコイル等で観測する方法である．それでは，わたしたちが扱う強相関電子系に於いて「超音波測定」とはどのような外場を加え，何の応答を観測できるのだろうか？その質問に大雑把に答えるとしたら

「超音波測定は物質中に歪み場を加え，電気四極子の応答を観測する手法である．」

と言えるだろう．

超音波は弾性波として固体中を伝搬する．局所的にその弾性波をみると，結晶中に然るべき対称性を持った歪み場が作り出されている．電子系はその歪み場を，ポテンシャルの変化として感じる．もし固体が完全結晶（※１） で周期的な格子を持つ場合，そのポテンシャルは結晶対称性によって周期函数で表される．超音波には縦波と横波が存在するため，様々な対称性の歪み場を加えることができる．超音波計測とは，いわば系のポテンシャルを外から直接揺さぶり，電子系（ならびに格子系）の応答を四極子感受率（あるいは歪み感受率）として観測する手法であり，磁気モーメントの応答に対応する帯磁率，エントロピーに対応する比熱とともに物性物理学における有効な測定手段の一つである．

そのため，超音波を用いて得られる「弾性定数」という基本的物理量は固体物理学の教科書には必ずと言っていいほど登場する．しかしながら，「磁性」や「誘電性」を観測する実験手法に比べ，超音波実験とそこから得られる物理には，正直なところ馴染みが少ないという学生諸君が多いのではなかろうか．確かにキッテル先生のIntroduction to Solid State Physicsでは一時期，弾性定数の章が割愛されていたし(※２），物性物理学の門を叩いた学生は自分の測定した比熱や磁化を手っ取り早く計算したいので，固体物理の教科書の「弾性」の章は読み飛ばしている可能性が高い．

そこで，本稿では学生の皆さんに「超音波で固体の電子状態を観る」ことにもっと馴染んでもらうべく，超音波で得られる物理量の持つ意味と，それを重い電子系や多極子秩序，ラットリング等を示す系に適用したときに得られるデータの解釈の仕方について，実験屋の視点から基礎的な解説をし，最後に最近著者が行っている充填スクッテルダイト化合物のラットリングの研究について，超音波実験から得られた知見を紹介する．

本稿が学会や論文で超音波の実験結果をみる際の手助けになれば幸いである.

(第二章に続く）

注釈
※１　 人類が手にすることができる最も「完全」に近い結晶はシリコンの単結晶である．しかし最近，極低温弾性定数測定によって１モル当り10¹⁴個程度の単原子空孔の存在が明らかになった. [1]
※２　7th Edition以降で復活

第２章　超音波実験の測定手法

2.1 パルス法

先ず，超音波実験の測定手法と測定の勘所を至極簡単に述べよう．超音波測定では，固体中を伝搬する弾性波の音速と吸収を測定する．パルス法と呼ばれる方法は固体の音速を測定する最もポピュラーな方法である．良く研磨された測定試料面に各種の接着剤によりトランスデューサを接着し，トランスデューサの共振周波数に合った連続波発信器を用意する．そこから出力された電気信号をゲートに通して充分に短い幅（200-500 [ns]）に切り，ドライブ・パルスとしてトランスデューサに入力する．トランスデューサによってパルス信号は超音波パルスに変換され，試料中を弾性波（歪み波）として伝搬する．試料の両端面で反射を繰り返した超音波パルスはトランスデューサにより再度電気信号に変換され，検出器へ送られる．オシロスコープでその信号の時間変化を監視すると，図1のような「超音波パルスエコー」が得られる．「エコー」とは山びこ現象で次々に聞こえる「こだま」そのものである．理想的なパルスエコーならばその振幅は exp(-βt)に比例して減衰する．t [s]は入射信号からの遅延時間である.

超音波が伝搬する経路長を，試料端で反射したエコーが届く時間間隔で割れば，超音波の音速$v$の絶対値が得られる．βは単位長さ当たりの超音波吸収係数α [dB m^-1]、音速v [m s^-1]との間にβ = α/vの関係を持つ．ここで，ドライブ・パルスはほぼ矩形のエンベロープを持っているが，電気音響変換の過程で歪みを生じたり、トランスデューサを試料に貼付ける為の接着剤による影響でエンベロープのエッジが丸みを帯びてくる。特に電気機械結合定数の高いLiNbO₃を圧電素子として用いた場合，エコー信号の立ち上がりを正確に知ることが難しくなるため，正確な絶対値の測定には直接試料表面に圧電結晶をスパッタリングする方法や圧電高分子トランスデューサをスピンコート法で塗布するなどの工夫が必要である.いずれにしても音速の絶対値の測定精度には限界があり，うまく測定しても数%の誤差が生じてしまう．

Fig. 01 超音波エコーの観測例：URu₂Si₂の横波弾性定数C₄₄モードの超音波エコーのスナップ写真．
（厚さ100 μmのLiNbO₃ウエーハをトランスデューサに用い，３倍高調波の105 MHzで観測した．）

Fig. 02 位相比較法で扱う信号のタイムチャート（概念図）

Fig. 03 超音波位相比較法のブロックダイアグラム

2.2 位相比較法

量子系の状態を反映したより微小な音速の変化を検出するためには，エコー信号を検出するのではなく，位相信号を用いることで高い分解能を達成することができる．これは位相比較法と呼ばれる手法である．位相比較法については様々なところで解説があるので[2,3] ，その詳細については省くが，ここでの勘所は音速の変化を基準信号からの遅延（位相差）として観測し（図２)，位相差をゼロ検出し，先程の発信器に負帰還をかけることで，音速の相対変化を信号発生器の周波数の相対変化に読み替えるところにある．即ち，原理的には測定周波数が高ければ高いほど，分解能が向上することになる．また，ゼロ検出法を用いることで，測定系における非線形性の影響を受けにくく，Δ v/v 〜 10^-7のきわめて高い分解能の測定が可能となる．この高い分解能は，特にSi単結晶内に存在する単原子空孔の量子状態の測定[4]や，音響ドハース効果などの極微小な量子振動の検出に有効である．我々が用いている測定系のブロックダイアグラムを図3に示す．これとは別に周波数を固定して生の位相差信号をデジタルストレージに保存し，そこから音速の相対変化と吸収係数を同時に算出する方法もある[5]．この方法は若干分解能で劣るものの，積分器による負帰還を行わないため，音速の変化が著しい場合やパルス磁場下など短い時定数で測定を行う必要がある場合に有効である．

(第３章に続く）

第３章　超音波実験の測定手法

3.5 多極子秩序における弾性異常

まずは歪みと線形結合する電気四極子が秩序を起こす典型物質の弾性応答を紹介しよう．反強四極子(AFQ)秩序を起こす物質の典型例といえば先ほど紹介したCeB₆である[24]．この物質はT_Q = 3.3 KでΓ₅型の四極子が秩序波数　k = [1/2 1/2 1/2]で整列するAFQ秩序を起こし，さらにT_N = 2.3 K以下で反強磁性(AFM)秩序を示すことが中性子散乱実験，NMR，共鳴X線散乱実験で確認されている．図13にCeB₆の弾性定数C₄₄の低温部拡大図を示す．Γ₈四重項結晶場基底状態はAFQ秩序の四極子-歪み相互作用によって２つのクラマース二重項に分裂する．そのためT_QでC₄₄のソフト化は止まり，低温でゆるやかに上昇する．T_Q以上の常磁性相では1/Tに比例したキュリー項があるため，Γ₈四重項が反強四極子秩序によって2つのクラマース2重項に分裂することで，四極子感受率-χ_Γ5が逆カスプ状に折れ曲がった結果としてこの弾性異常が理解できる．これは反強磁性磁化率において副格子磁化に平行に磁場を加えたときの磁化率χ_||に現れるカスプ状の応答を垂直方向にひっくり返したものと類推することができる．

同じような振る舞いは正方晶DyB₂C₂の弾性応答にも観られる．DyB₂C₂は正方晶系で初めて反強四極子秩序が報告された物質である[25]．T_Q = 24.7 KでO_xy型のAFQ秩序を示し，T_NN = 15.3 KでAFM秩序を示すことが共鳴X線散乱実験や磁場下における中性子散乱実験で検証された [26, 27]．Dy³⁺はJ = 15/2のクラマースイオンでありDyのサイトシンメトリーC_4hの下で4E_1/2(+)4E_3/2と８つのクラマース二重項に既約分解される．比熱で見積もられた磁気エントロピーがT_QでRln4に達することから，結晶場基底状態は2組のクラマース二重項が擬縮退した擬四重項基底状態になっていると考えられる．これらを踏まえて弾性定数を観ると，全ての弾性定数にソフト化が現れていることから [28]，選択則より基底状態は対称性の異なる二つのクラマース二重項の組み合わせに限定される．さらにT_Qにおいて秩序変数O_xyの応答を観る弾性定数C₆₆にはCeB₆で弾性定数C₄₄に観られたような）逆カスプ状の折れ曲がりが観測される．常磁性相におけるソフト化の変化量を比較すると四極子O_yz, O_zxの応答に対応する弾性定数C₄₄が最も大きい．これは四極子-歪み相互作用の結合定数の大小によって決まっており，AFQ 秩序の秩序変数の対称性と常磁性相におけるソフト化の大きさの間に相関は無い.

図13 CeB₆の弾性定数C₄₄の低温部拡大図 [18]

図14 DyB₂C₂の弾性定数の温度変化 [39]

一方，強四極子（FQ)秩序を示す物質の場合は，秩序変数とそれに対応する弾性定数に現れるソフト化の大小に相関がある．図14にHoB₆の弾性定数の温度依存性を示す．この物質はT_Q = 6.1 KでΓ₅型の四極子O_yz, O_zx, O_xyが<111>方向に整列するFQ秩序を示す．対応する弾性定数C₄₄は協力的ヤーン・テラー効果による三方晶への構造変化に伴い，T_Qに向かって弾性定数が発散していることがわかる．これは強磁性転移において自発磁化が生じ容易軸方向の磁化率が発散することと類推できる．図15にCe₃Pd₂₀Ge₆ の弾性定数の温度変化を示す．この物質はT_Q = 1.3 Kで正方晶もしくは斜方晶への構造相転移が起こっていることが熱膨張と中性子回折実験で確認されている [29,30]．弾性定数(C₁₁-C₁₂)/2$に50%を超えるソフト化が観測されることからΓ₃対称性の四極子O₂⁰あるいはO₂²がFQ秩序していることを強く示唆する.

図15 HoB₆の弾性定数C₁₁, C_L, (C₁₁-C₁₂)/2, C₄₄の温度変化 [18]

図16 Ce$_3$Pd$_{20}$Ge$_6$の弾性定数$C_{11}$, $C_{¥rm L}$, $C_{¥rm B}$, $(C_{11}-C_{12})/2$, $C_{44}$の温度変化[57]

さて，四極子以上のランクの多極子と直接結合する共役場のプローブは今のところ見つかっていない．そのため，磁気八極子秩序や電気十六極子秩序が起こっていると考えられている物質に対しては，磁場や歪み場によって誘起される下位の電気四極子や磁気双極子を観測することで傍証を集めるしかない．現在，八極子秩序が実験的に立証されている系で超音波の報告があるのはCe_xLa_1-xB₆だけである．CeB₆をLaで希釈していくと磁気相互作用と四極子相互作用が拮抗し，x = 0.75でAFQ転移温度とAFM転移温度が逆転した温度領域にO_yz, O_zx, O_xyの強四極子モーメントが発生する非磁性のIV相と呼ばれる新たな相が現れる [31,32]．超音波測定によって弾性定数C₄₄は強四極子相関により31%の巨大なソフト化を示すことがわかり [33]，赤津らによる熱膨張実験によって結晶がバルクで三方晶に歪んでいることが明らかとなった [34]．これらの実験事実はIV相で強四極子モーメントが発達していることを強く示唆する．後に，この強四極子モーメントは久保・倉本によってΓ_5u型の反強八極子秩序で二次的に誘起され得ることが理論的に示された [35,36]．その後，共鳴X線散乱実験や磁場中中性子散乱実験が行われ，それぞれ反強八極子モーメントの誘起する電気四極子と磁気双極子の超格子反射が観測され，微視的に反強八極子秩序が実証された[37,38]．

図17 Ce_0.25La_0.75B₆の弾性定数C₁₁, C_B, (C₁₁-C₁₂)/2, C₄₄の温度依存性 [31,33]

(第３章3.6節に続く）