通常の水素である軽水素は、原子核が陽子 1 個だけで構成される。その質量は他の元素の原子核に比べて小さいため、金属結晶の中を移動するとき、熱エネルギーによって障壁を越える経路だけでなく、量子トンネルによって障壁を通り抜ける経路も取りうる。ただし、水素が軽いという事実だけでは、実際の移動速度も、量子トンネルが支配的になる温度も決まらない。東京大学生産技術研究所などの研究グループは、バナジウム薄膜中の水素について、低水素濃度の相では 35 K という低温でも高速な移動が起こる一方、高水素濃度の相では量子トンネルが抑制され、熱エネルギーを受け取って格子間を移る拡散が支配的になることを示した[1][2]。
発表資料は、この差を、水素が「波」として振る舞う場合と「粒子」として振る舞う場合として紹介している。この表現は現象の違いを直感的に伝えるが、水素原子核が量子力学に従う状態と、量子力学に従わない状態へ切り替わるわけではない。どちらの結晶相でも、水素原子核の状態は波動関数によって記述される。変化するのは、その波動関数の空間的な広がりである。低水素濃度の結晶では、波動関数が複数の格子間サイトにまたがって広がる。水素濃度が増えて結晶が一軸方向へひずむと、波動関数は特定のサイト付近へ集中する。
ここでいう局在化は、水素原子核の位置が測定によって一つに確定したという意味ではない。結晶が作るポテンシャルエネルギーの地形により、低いエネルギーを持つ状態の確率分布が、限られた格子間サイトの周囲へ偏ることを指す。測定による状態の確定を扱う量子測定問題との違いは、既稿で整理した[3]。結晶構造によって波動関数の広がりが制限される局在化とは、原因も論点も異なる。
研究結果をつなぐ因果関係は、水素の軽さから拡散速度へ直接進むものではない。第一段階では、水素濃度の上昇がバナジウム格子を一軸方向へ変形させる。第二段階では、その変形が結晶の対称性を低下させ、それまで等価だった格子間サイトにエネルギー差を生じさせる。さらに、移動元と移動先の量子状態がそろわなくなることで、波動関数は複数サイトへ広がりにくくなり、サイト間の量子トンネルが抑えられる。最終的に観測されるのが、低温でも高速に移動する相と、熱を受け取るまで移動しにくい相との差である。
本稿の焦点は、量子トンネルを水素原子核が単独で所有する固定的な能力として捉えるのではなく、結晶がどのような移動先と経路を用意しているかによって発現条件が変わる現象として読むことにある。水素の質量は結晶相が変わっても同じである。それにもかかわらず移動機構が変わる以上、説明すべき対象は水素原子核だけではない。水素濃度、格子ひずみ、サイトの等価性、エネルギー準位、波動関数の広がりを、一つの連続した構造として追う必要がある。
1. 水素の軽さだけでは量子拡散を説明できない
金属結晶では、金属原子が一定の規則に従って並んでいる。水素は、多くの場合、金属原子を押し出してその位置を占めるのではなく、金属原子の間に形成された格子間サイトへ入る。格子間サイトは単なる空洞ではない。周囲の金属原子との距離や配置によって、水素がそこに存在するときのエネルギーが決まり、隣のサイトへ移るときの経路と障壁も決まる。
水素が格子間サイトへ入ると、周囲の金属原子は水素との相互作用によって元の位置から変位する。水素が少ない段階では、局所的なひずみとして吸収できる場合がある。濃度が上がると、個々のひずみが独立した局所変形では済まなくなり、結晶軸の長さや安定な水素位置が系統的に変わる。さらに濃度が増えれば、結晶全体が別の相へ移る。金属中水素の拡散を考えるとき、水素原子核の質量だけでなく、格子ひずみ、欠陥、占有サイト、相構造が必要になるのは、このように水素の存在自体が移動環境を変えるためである[4]。
熱による拡散では、水素原子核が格子振動からエネルギーを受け取り、現在いるサイトと隣のサイトの間にある障壁を越える。温度が下がれば格子振動から受け取れるエネルギーが減るため、障壁越えの頻度も急速に低下する。古典的な熱活性化型拡散だけを考えれば、十分に低い温度では水素の移動はほぼ止まるはずである。
ところが、軽い水素原子核では、熱によって障壁の頂上まで上がらなくても、波動関数が障壁の反対側まで達する場合がある。タンタル中の水素を扱った回復実験では、低温での移動を量子拡散として解釈する根拠が示され、軽い格子間原子の拡散理論では、格子との相互作用を含めたトンネル過程が検討されてきた[5][6]。水素を、質量が約 2 倍の重水素へ置き換えると拡散挙動が変わることも、原子核の質量と波動関数の広がりが移動へ関与する証拠となる[7]。
これらの研究から、水素原子核の小さな質量が、低温で観測可能な量子拡散を生じやすくする主要因であることが分かる。ただし、質量だけでは実際の拡散機構は決まらない。水素原子核が軽くても、移動先までの距離が長ければ波動関数の重なりは小さくなる。障壁が高いだけでなく幅も広ければ、反対側へ到達する振幅は急速に減衰する。移動元と移動先の量子準位に大きな差があれば、二つの状態は有効に混成できない。水素が移る際に周囲の金属原子を大きく動かす必要があれば、その格子変形も移動を妨げる。
| 条件 | トンネルを促す状態 | トンネルを抑える状態 |
|---|---|---|
| 原子核の質量 | 水素のように軽い原子核では、波動関数が空間的に広がりやすい。 | 重い同位体や元素では、同じ障壁に対する透過確率が低下する。 |
| 移動距離 | 隣接サイトが近ければ、両側の波動関数が重なりやすい。 | サイト間距離が長ければ、障壁内で波動関数が大きく減衰する。 |
| 障壁の形状 | 低く狭い障壁では、反対側へ到達する確率が高くなる。 | 高い障壁や幅の広い障壁では、トンネル確率が低下する。 |
| 量子準位 | 移動元と移動先のエネルギーが近ければ、状態が混成しやすい。 | サイト間にエネルギー差があれば、共鳴的な移動が成立しにくい。 |
| 格子の応答 | 小さな格子変形で移動できれば、トンネル経路が保たれやすい。 | 周囲の金属原子を大きく動かす必要があれば、移動に追加のエネルギーが必要になる。 |
パラジウム中の水素では、八面体サイトと四面体サイトの量子準位が近づく条件で、両サイト間の共鳴トンネルが生じることが示された[8]。この結果は、障壁の高さだけでなく、始点と終点の量子状態が整合していることが移動を左右することを具体的に示している。続く研究では、水素の移動速度と格子振動、伝導電子との相関が測定され、量子トンネルが水素原子核だけで完結する運動ではなく、金属格子と電子系を含む環境との相互作用の中で成立することが確認された[9]。
バナジウム研究が加えたのは、材料の種類を変えた比較ではない。同じバナジウム薄膜の中で、水素濃度を変えると結晶相が変わり、その相変化に対応して量子トンネルが現れる条件と抑えられる条件が分かれることを示した。水素濃度が母体格子を変形させ、変形した格子が水素の安定サイトと量子準位を変える。その結果として、同じ質量を持つ水素原子核の移動機構が、量子トンネルから熱活性化型拡散へ切り替わる。
水素の軽さは、量子トンネルが起こりうる理由を説明する。しかし、α相では実際に高速移動が現れ、β相では抑えられる理由までは説明できない。その差を生むのは、バナジウム格子がどのような対称性を持ち、水素にどの格子間サイトを与えるかである。次に確認すべきなのは、低水素濃度の α相と高水素濃度の β相で、水素の居場所と結晶構造がどのように変わるかである。
2. α相と β相では水素の居場所が変わる
バナジウムは、常温では体心立方格子を取る。立方体の 8 個の頂点と中心にバナジウム原子が配置された構造であり、金属原子の間には水素が入りうる格子間サイトが存在する。代表的なのが、4 個のバナジウム原子に囲まれた四面体サイトと、6 個のバナジウム原子との位置関係で定まる八面体サイトである。名称は空隙の外形を示すだけではない。周囲の原子数と距離が異なるため、水素が感じるポテンシャルエネルギー、原子核の振動状態、隣接サイトまでの移動距離も異なる。
水素がどちらのサイトを選ぶかは、バナジウムに固有の一つの答えとして決まっているわけではない。水素濃度、軽水素と重水素の違い、格子の伸び方によって、四面体サイトと八面体サイトの相対的な安定性が変わる。先行研究でも、バナジウム中の水素と重水素について、占有サイトと量子状態が格子ひずみに応じて変化することが調べられてきた[10]。
水素濃度が低い α相では、バナジウム格子は水素を固溶させながら、ほぼ体心立方構造を保つ。この構造では、複数の四面体サイトが立方対称性によって等価になり、水素は主として四面体サイトを占有する。ある四面体サイトだけが特別に深い安定点になるのではなく、結晶の回転や反転によって互いに対応する複数のサイトが、ほぼ同じエネルギーを持つ。
水素濃度が上がると、格子は局所的に膨張するだけでは済まなくなる。水素を含む領域全体で原子配置が組み替わり、バナジウム水素化物は別の結晶相を形成する。β相では、体心立方格子の 1 本の軸が他の 2 本より長くなり、体心正方構造へ変わる。バナジウム水素系とバナジウム重水素系の相図は、形成される相が水素濃度と温度の双方に依存することを示している[11]。金属水素系で起こる相転移は、水素が既存の空隙へ受動的に蓄積する現象ではなく、水素と母体金属の相互作用によって結晶構造そのものが変わる過程である[12]。
| 観点 | α相 | β相 |
|---|---|---|
| 水素濃度 | 水素を低濃度で固溶した状態として形成される。 | 水素濃度の上昇に伴って形成される水素化物相である。 |
| 結晶構造 | 立方対称性を持つ体心立方構造をほぼ保つ。 | 1 本の結晶軸が伸びた体心正方構造になる。 |
| サイトの等価性 | 対称操作で対応する複数の四面体サイトが、ほぼ同じエネルギーを持つ。 | 一軸ひずみによってサイト間の等価性が失われ、エネルギー差が生じる。 |
| 主な水素位置 | 四面体サイトが安定位置になる。 | 伸長軸方向に位置する八面体サイトが安定化する。 |
| 低温での移動 | 複数サイトにまたがる量子状態が形成され、量子トンネルによる高速移動が可能になる。 | 波動関数が特定サイト付近へ局在し、熱活性化型の移動が支配的になる。 |
α相から β相への変化では、水素の居場所が四面体サイトから八面体サイトへ移るだけではない。第一段階として、一軸ひずみによって立方対称性が失われる。第二段階として、それまで等価だった格子間サイトのエネルギーが分裂し、伸長軸方向の八面体サイトが選択的に安定化する。水素原子核は、同じ条件を持つ複数の位置から移動先を選べる状態ではなく、特定の位置から別の非等価な位置へ移らなければならない状態へ置かれる。
占有サイトは、水素が静止している位置を示すだけの構造情報ではない。あるサイトから次のサイトまでの距離、途中にあるポテンシャル障壁、移動先とのエネルギー差は、出発地点の種類によって変わる。四面体サイトが連続して作る移動網と、八面体サイトを起点とする移動網では、水素原子核が通過しなければならない経路が異なる。サイト占有の変化は、そのまま拡散経路の組み替えを意味する。
バナジウム中では、四面体サイトと八面体サイトの占有率を変えると拡散速度も変化することが、計算と実験の比較から示されている[13]。中性子分光では、水素原子核の振動と拡散に対応する運動が調べられ、バナジウム、ニオブ、タンタルの間で異なる低温挙動が観測された[14]。核磁気共鳴による研究でも、バナジウム水素系の相転移に伴って、水素の移動速度と温度依存性が通常の熱拡散から外れる領域が確認されている[15]。これらの先行研究は、結晶相、占有サイト、拡散が独立した物性ではなく、同じ格子構造から生じる相互に結びついた現象であることを示している。
α相と β相の速度差を、水素原子数が増えたことで混雑が生じた結果と考えるだけでは、低温で移動機構そのものが変わる理由を説明できない。水素濃度の上昇は、まずバナジウム格子を一軸方向へ変形させる。その変形が水素の安定サイトとサイト間のエネルギー関係を変え、利用可能な移動経路を組み替える。濃度と拡散の間にある結晶構造の変化を追うことで、同じ水素原子核が α相では量子トンネルによって移動し、β相では熱を受け取るまで局在する理由が見えてくる。
3. 一軸ひずみが等価だった移動先を分ける
体心立方構造では、結晶の x、y、z 方向は回転対称性によって互いに入れ替えられる。ある方向だけが長い、あるいは特定の面だけが特別であるという区別がないため、対称操作によって対応づけられる四面体サイトは、同じ局所構造とほぼ同じエネルギーを持つ。水素原子核から見れば、複数の移動先が同じ条件で並び、どのサイトへ移っても結晶から受ける作用が変わらない状態である。
サイトが等価であれば、それぞれの位置に局在した量子状態も同じエネルギーを持つ。隣接サイト間の距離が短く、途中の障壁を通じて波動関数が重なれば、各サイトに対応する状態は独立したままではなく、複数サイトにまたがる固有状態へ組み替えられる。水素原子核が一つのサイトから別のサイトへ古典的に飛び移る前から、その存在確率が複数の位置へ広がるのは、この等価な状態同士の混成による。
β相では、この条件が一軸ひずみによって崩れる。格子の 1 本の軸が伸びると、伸長軸に沿った方向と、それに垂直な面内方向は対称操作で入れ替えられなくなる。周囲のバナジウム原子までの距離が方向ごとに変わるため、同じ種類に見えた格子間サイトでも、水素が受けるポテンシャルエネルギーは一致しない。
密度汎関数理論による先行研究では、バナジウム格子の軸比 c/a が増えるにつれて、四面体サイトと八面体サイトの相対的な安定性が逆転し、伸長軸方向の八面体サイトが水素の安定位置になることが示されている[16]。格子の伸長は、まず周囲のバナジウム原子との距離を変える。その距離変化によって各サイトのポテンシャルエネルギーと原子核の振動準位が変わり、幾何学的なひずみが量子状態のエネルギー差へ変換される。
基板に固定されたバナジウム薄膜では、この構造変化をバルク材料と同じようには緩和できない。膜面内の格子間隔は基板との整合によって拘束されるため、水素吸収による膨張は主として膜厚方向へ現れる。水素濃度が増えるほど一軸的なひずみが強まり、相構造と拡散係数が連動して変化することが報告されている[17]。今回用いられた MgO 基板上のエピタキシャル薄膜でも、水素濃度は単なる粒子数ではなく、格子の変形条件を決める変数になる。
| 段階 | 立方対称性を保つ状態 | 一軸ひずみを持つ状態 |
|---|---|---|
| 原子配置 | x、y、z 方向が対称操作によって入れ替えられる。 | 伸長軸と面内方向を同じものとして扱えなくなる。 |
| 格子間サイト | 複数の四面体サイトがほぼ同じ局所構造を持つ。 | 八面体サイトと四面体サイト、方向の異なるサイトに構造差が生じる。 |
| ポテンシャルエネルギー | 対称性で対応するサイトが同じ深さの安定点になる。 | 伸長軸方向の八面体サイトが選択的に低い安定点になる。 |
| 量子準位 | 複数サイトの低エネルギー状態がそろい、相互に混成できる。 | サイトごとの準位が分裂し、異なる位置の状態が混成しにくくなる。 |
| 波動関数 | 複数の四面体サイトへ広がる非局在化状態を形成できる。 | 最も安定な八面体サイトの周囲へ確率分布が偏る。 |
| 移動への帰結 | 等価なサイト網を通じて量子トンネルが拡散へ寄与できる。 | 非等価なサイト間の移動には熱による励起が必要になる。 |
対称性の低下は、同じ形の空隙にわずかな差を付けるだけではない。複数の移動先を同じエネルギーで接続していた条件を失わせる。あるサイトが他より深い安定点になると、水素原子核の低エネルギー状態はその位置へ偏る。別のサイトへ移るには、空間的な障壁を通過するだけでなく、異なるエネルギーを持つ状態へ移らなければならない。
この変化を「対称性が下がると波動関数が局在する」とだけまとめると、間の機構が抜け落ちる。直接の原因は、ひずみによってサイト間のエネルギー差が生じることである。その背後には、結晶対称性が複数の原子配置を等価に保ち、対応する量子準位をそろえていた構造がある。対称性が破れると準位の縮退が解け、状態の混成が弱まり、波動関数が特定サイトの周囲へ集中する。
水素濃度と拡散速度の関係にも、この中間過程が入る。水素が増えると格子が一軸方向へひずむ。ひずみによって格子間サイトの等価性が失われる。非等価になったサイト間では波動関数が連続して広がれず、トンネルによる移動が抑えられる。高濃度で拡散が遅くなるのは、水素同士が空間的に混雑するためだけではなく、水素の蓄積が移動経路の量子力学的な接続条件を変えるためである。
移動する水素と、その水素を受け入れる結晶は、固定された主体と環境に分けられない。水素濃度が結晶構造を変え、変化した結晶構造が次の水素の安定位置と移動可能性を決める。この循環によって、α相で開かれていた等価なサイト間の経路は、β相ではエネルギー差を持つ経路へ組み替えられる。次に確認すべきなのは、こうして非等価になったサイト間で、障壁の高さだけでは量子トンネルの有無を説明できない理由である。
4. 障壁の高さだけを見てもトンネル確率は決まらない
結晶中の水素拡散は、安定な格子間サイトから隣のサイトへ移る過程として記述できる。各サイトは水素にとってポテンシャルエネルギーの谷に相当し、サイト間にはエネルギーの高い領域がある。熱活性化型の拡散では、水素原子核が格子振動からエネルギーを受け取り、この障壁を上から越える。量子トンネルでは、障壁の頂上まで到達しなくても、波動関数が障壁内部で減衰しながら反対側へ達する。
バナジウム中の水素については、第一原理計算により、金属内部への水素の溶解、格子間サイトを結ぶ拡散経路、表面で解離した水素が内部へ入るまでの障壁が調べられてきた[18][19]。これらの計算は、水素の移動が始点と終点だけで決まらず、その間で周囲のバナジウム原子がどのように変位し、どの格子間位置を通過するかによって変わることを示している。
β相に相当する体心正方構造では、水素は伸長軸方向の八面体サイトに安定して存在する。別の八面体サイトへ移る経路には、両サイトを直接結ぶ経路だけでなく、四面体サイトなどの中間位置を経由する経路がある。体心正方構造の一水素化物相を扱った先行研究でも、八面体サイトの占有と、四面体サイトを介した拡散過程が検討されてきた[20]。
今回の研究では、α相と β相を模した構造について、水素の移動経路に沿ったポテンシャルエネルギーが計算された。水素だけを動かして周囲のバナジウム原子を固定した場合と、移動位置ごとに格子を緩和させた場合が比較されている。格子緩和を含めるのは、水素が移動するとき、周囲の金属原子も新しい水素位置に応じて変位するためである。原子配置を固定した計算だけでは、水素と格子が協調して移動障壁を下げる効果を見落とす。
立方対称性を保つ α相の模型では、隣接する四面体サイト間を直接結ぶ経路の障壁は、格子緩和を含めると 0.122 eV となった。電子ボルトは原子や電子のエネルギーを表す単位であり、1 eV は電子 1 個を 1 V の電位差で動かすときのエネルギーに相当する。この数値だけでトンネル確率を定量できるわけではないが、β相の経路と比べる基準になる[1]。
一軸ひずみを持つ β相の模型では、安定な八面体サイト間を直接結ぶ経路の障壁が、格子緩和後も 1.025 eV に達した。α相の四面体サイト間障壁の 8 倍を超える。八面体サイト間の距離も長いため、障壁は高いだけでなく、水素原子核が通過しなければならない領域も広い。波動関数は障壁内部で距離とともに減衰するため、高さと幅の両方が増える直接経路では、低温でのトンネル確率が著しく小さくなる[1]。
β相には、四面体サイトを経由することで、直接経路より低い障壁を通る間接経路も存在する。計算された最大障壁は経路によって 0.228 eV または 0.376 eV であり、1.025 eV の直接経路より低い。この結果だけを見ると、水素は八面体サイトから四面体サイトへトンネルし、そこから次の八面体サイトへ移れるようにも見える。しかし、間接経路が存在することと、経路全体を量子トンネルで通過できることは同じではない[1]。
| 移動経路 | 結晶状態 | 計算された障壁 | トンネルを左右する条件 | 低温移動への帰結 |
|---|---|---|---|---|
| 四面体サイト間の直接経路 | 立方対称性を持つ α相を模した構造である。 | 格子緩和を含めると 0.122 eV である。 | サイト間距離が短く、始点と終点が対称性によって等価である。 | 両サイトの状態が混成し、波動関数が複数サイトへ広がりやすい。 |
| 八面体サイト間の直接経路 | 一軸ひずみを持つ β相を模した構造である。 | 格子緩和を含めても 1.025 eV である。 | 障壁が高く、サイト間距離も長いため、障壁内部で波動関数が強く減衰する。 | 低温で直接トンネルする確率は極めて小さくなる。 |
| 四面体サイトを経由する間接経路 | β相を模した構造内で八面体サイトと中間サイトを結ぶ。 | 経路によって最大 0.228 eV または 0.376 eV である。 | 障壁は直接経路より低いが、始点、経由地、終点の量子準位が一致しない。 | 経路全体を覆う非局在化状態が形成されず、連続的なトンネル移動になりにくい。 |
量子トンネルを理解するには、古典的な障壁の高さと、各サイトに形成される量子状態を分けて見る必要がある。水素原子核は、ポテンシャルエネルギーの谷の底に静止しているわけではなく、最も低い状態でも零点運動を持つ。各サイトには、その局所的な形状に応じた離散的な振動準位が形成される。隣接する 2 個のサイトが同じ構造とエネルギーを持てば、両側の低エネルギー状態は混成し、二つのサイトにまたがる固有状態を作れる。
α相の四面体サイトは、立方対称性によって互いに等価である。始点と終点のポテンシャルの深さ、局所構造、低い振動準位がそろうため、障壁を挟んだ波動関数が同じエネルギーで接続される。第一段階として、短いサイト間距離と比較的低い障壁が、両側の波動関数に有限の重なりを与える。第二段階として、等価なサイトの準位整合が、その重なりを複数サイトに広がる固有状態へ変える。低温で観測される高速移動は、この二つの条件が同時に成立した結果である。
β相の間接経路では、障壁の一部が低くなっても、八面体サイトと経由する四面体サイトのポテンシャルの深さが異なる。一軸ひずみによって両サイトの局所構造が分かれ、各サイトに形成される水素原子核の準位もそろわない。始点の低エネルギー状態から経由地の状態へ移るには、空間的な障壁を透過するだけでなく、エネルギー差も埋めなければならない。
温度が十分に低ければ、水素原子核はそのエネルギー差を格子振動から受け取れない。始点と経由地の波動関数に空間的な重なりが残っていても、異なるエネルギーを持つ状態は強く混成しない。結果として、八面体サイト、四面体サイト、次の八面体サイトを一続きに覆う固有状態は形成されず、水素は最初の八面体サイト付近へ局在する。
| 判断要素 | 物理的な役割 | α相 | β相 |
|---|---|---|---|
| 障壁の高さ | 障壁内部で波動関数が減衰する強さを決める。 | 四面体サイト間の障壁は比較的低い。 | 八面体サイト間の直接障壁は高い。 |
| 障壁の幅 | 波動関数が減衰しながら進む距離を決める。 | 隣接する四面体サイト間の距離が短い。 | 安定な八面体サイト間の距離が長い。 |
| サイトの等価性 | 始点と終点に同じ量子状態を形成できるかを決める。 | 立方対称性によって複数サイトが等価になる。 | 一軸ひずみによってサイトが非等価になる。 |
| 量子準位の整合 | 異なるサイトの状態が一つの固有状態へ混成できるかを決める。 | 低い準位がそろい、非局在化状態を形成できる。 | 準位差が生じ、経路全体の状態が混成しにくい。 |
| 格子緩和 | 水素移動に伴う金属原子の変位が障壁をどこまで下げるかを決める。 | 緩和によって低い移動障壁が形成される。 | 緩和後も直接経路の高い障壁とサイト間の非等価性が残る。 |
α相の熱拡散係数は先行研究で測定されており、今回の研究では、その高温領域で得られた温度依存性を 35 K まで外挿した移動距離と、薄膜内で実際に観測された水素分布の変化が比較された[21]。熱活性化型の拡散だけなら、35 K における移動距離は極めて短くなる。それにもかかわらず、水素は数秒以内に厚さ 40 nm の膜内へ広がった。低い障壁、短い距離、等価なサイト、そろった量子準位という α相の構造が、この熱拡散からの大幅な逸脱を説明する。
β相でトンネルが抑えられる原因を、1.025 eV という高い直接障壁だけに求めると、より低い間接経路が存在する事実を説明できない。間接経路が低温で機能しないのは、一軸ひずみによって経由地のエネルギーが始点と異なり、量子準位の連続性が失われるためである。障壁の低下は移動の必要条件の一つにすぎず、経路を構成する状態が同じエネルギーで接続されなければ、量子トンネルは実際の拡散へ発展しない。
α相から β相への変化は、移動経路に高い壁が追加された現象ではない。立方対称性の低下によって、複数のサイトを同じ量子状態で結んでいた条件が失われた変化である。α相では、等価な四面体サイトが一つの量子力学的な移動網を作る。β相では、八面体サイトと中間サイトの準位が分かれ、経路は個別の熱活性化過程へ分解される。量子トンネルを決めているのは、障壁の高さだけでなく、結晶が経路全体にわたって量子状態の連続性を保てるかどうかである。
5. 水素の位置と移動を別々に測る
バナジウム中の水素拡散を明らかにするには、水素がどこに存在するかと、どの速さで移動するかを分けて測る必要がある。水素の深さ方向分布が変化しても、それだけでは占有している格子間サイトまでは分からない。反対に、水素が八面体サイトまたは四面体サイトにいることを確認しても、隣接サイトへ移る時間や温度依存性は決まらない。本研究では、位置、分布、速度を異なる測定で取得し、量子状態計算によって相互の関係を接続した。
水素の位置測定が難しい理由は、水素原子が電子を 1 個しか持たないことにある。X 線回折では、電子による散乱を通じて原子配置を調べる。多数の電子を持つバナジウム原子は強い信号を与えるが、水素の寄与は小さい。重い金属原子に囲まれた格子間水素について、通常の X 線回折だけから深さ方向の分布や占有サイトを詳しく特定することは難しい。
研究グループは、この制約に対して、共鳴核反応法とイオンチャネリングを組み合わせた。共鳴核反応法は電子ではなく水素原子核との核反応を利用するため、金属原子からの強い電子散乱に埋もれず、水素を選択的に検出できる。さらに、イオンの入射方向を結晶軸に合わせることで、水素が結晶内のどの位置にあるかを判別する。この方法は、エピタキシャル成長させたチタン水素化物薄膜について、水素と重水素のサイト占有差を調べる研究でも用いられている[22]。
5.1 共鳴核反応法で深さ方向の分布を測る
共鳴核反応法では、高エネルギーの窒素イオンを試料へ入射し、水素原子核との核反応によって放出される γ線を測定する。核反応は、窒素イオンが特定のエネルギーを持つときに強く起こる。入射したイオンは試料内部を進むにつれてエネルギーを失うため、入射時のエネルギーを変えると、共鳴条件に達する深さも変わる。
この性質を利用すると、表面からどの深さに水素が存在するかを測定できる。ある入射エネルギーで強い γ線が検出されれば、そのエネルギーに対応する深さに水素が多いことを意味する。入射エネルギーを順に変えて測定すれば、水素濃度の深さ方向分布が得られる。
本研究では、厚さ 40 nm のバナジウム薄膜へ 35 K で水素イオンを注入した。注入直後の水素は膜全体へ均一に分布せず、入射した表面に近い領域へ偏る。試料温度を上げた後に深さ分布を再測定すると、水素が膜の内部へ移り、分布が均一化したことを確認できる。共鳴核反応法が直接示すのは、この注入位置からの再分布である[1]。
5.2 イオンチャネリングで格子間サイトを識別する
深さ分布だけでは、水素が四面体サイトと八面体サイトのどちらを占有しているかは分からない。そこで、窒素イオンを結晶の特定方向から入射するイオンチャネリングを組み合わせる。
完全に不規則な方向からイオンを入射すると、イオンは多数のバナジウム原子と衝突しながら進む。結晶軸や結晶面に沿って入射すると、規則的に並んだ金属原子列の間に通路が生じ、イオンは原子との衝突を避けながら深く進む。この現象がチャネリングである。
格子間サイトの位置によって、チャネリング中の窒素イオンが通る軌道との重なり方は異なる。水素が特定のサイトを占有すると、入射方向ごとの核反応強度に固有の変化が現れる。実測した角度依存性を、四面体サイトや八面体サイトを仮定した計算結果と比較することで、主な占有位置を判定できる。
高水素濃度の β相では、この測定から、水素が一軸伸長方向に対応する八面体サイトを主として占有することが確認された。これは、結晶の一軸ひずみによって八面体サイトが安定化するという第一原理計算と一致する。格子変形と占有サイトの関係は、計算上の予測だけでなく、位置測定によって裏づけられた。
5.3 電気抵抗の時間変化から拡散速度を求める
共鳴核反応法は、水素が再分布した前後の状態を測定できるが、分布が変化する時間過程を連続的に追う用途には制約がある。拡散速度の測定には、バナジウム薄膜の電気抵抗が水素濃度と分布に応じて変化する性質を利用した。
35 K で水素を注入した直後には、膜の表面側に水素濃度の高い領域が形成され、深部には水素の少ないバナジウム領域が残る。電流は膜面に沿って流れるため、表面側の水素化された層と、深部の水素濃度が低い層は、電気的には並列に接続された伝導経路として働く。
温度を上げると、水素は濃度の高い表面側から深部へ拡散する。膜厚方向の濃度差が小さくなるにつれて、局所的に異なる抵抗を持つ層が並列に電流を流す状態から、膜全体がほぼ同じ水素濃度を持つ状態へ移る。この再分布に伴って試料全体の抵抗が増加する。
抵抗値そのものが水素原子の移動距離を直接表すわけではない。研究グループは、水素濃度と局所抵抗の関係を用い、膜厚方向の水素移動を 1 次元の無作為歩行として計算した。仮定した拡散係数ごとに抵抗の時間変化を求め、実測された緩和曲線と比較することで、β相の拡散係数を推定した。
この解析には、膜面内では水素濃度が一様であり、主な濃度勾配が膜厚方向にあるという近似が入る。抵抗緩和は拡散速度を反映するが、拡散係数は測定値から直接読み取った量ではなく、濃度分布と電気伝導を結ぶ模型を介して得られた値である。
5.4 測定結果と量子状態計算を接続する
| 測定または解析 | 直接得られる情報 | 模型または比較を介して得られる情報 | 研究全体で支える主張 |
|---|---|---|---|
| X 線回折 | バナジウム薄膜の格子定数と一軸方向の変形を測る。 | 水素濃度の上昇に伴う α相から β相への構造変化を判定する。 | 高水素濃度では立方対称性が低下する。 |
| 共鳴核反応法 | 水素濃度の深さ方向分布を測る。 | 注入直後の偏った分布と、加熱後の均一化を比較する。 | 水素が膜厚方向へ実際に移動したことを示す。 |
| イオンチャネリング | 入射方向に対する核反応強度の変化を測る。 | 候補サイトごとの計算結果と比較し、占有位置を判定する。 | β相では伸長軸方向の八面体サイトが安定化する。 |
| 電気抵抗緩和 | 水素再分布に伴う抵抗の時間変化を測る。 | 1 次元拡散模型との比較から拡散係数を推定する。 | β相の拡散が温度依存の熱活性化過程であることを示す。 |
| 量子状態計算 | 仮定した結晶構造に対するポテンシャルエネルギー面と固有状態を求める。 | 波動関数が複数サイトへ広がるか、特定サイトへ局在するかを判定する。 | α相の高速移動と β相のトンネル抑制を結晶対称性から説明する。 |
各手法が示す範囲を区別すると、研究の論理が明確になる。共鳴核反応法は水素の再分布を示すが、移動機構までは決めない。イオンチャネリングは占有サイトを示すが、移動速度までは測らない。電気抵抗緩和は速度の温度依存性を与えるが、原子核の波動関数を直接観測するものではない。量子状態計算は非局在化と局在化を説明するが、計算だけでは実際の薄膜中で水素が移動したことを証明できない。
X 線回折とイオンチャネリングによって、水素濃度の上昇が結晶をひずませ、安定な格子間サイトを四面体サイトから八面体サイトへ変えることを確認する。共鳴核反応法と電気抵抗緩和では、その構造変化に対応して水素の再分布速度と温度依存性が変わることを調べる。さらに量子状態計算を重ねることで、サイトの非等価化が波動関数の局在化を生み、量子トンネルを抑える機構へ接続する。
本研究の強さは、一つの測定値を量子トンネルの直接証拠として扱ったことにはない。結晶構造、水素の格子位置、深さ方向分布、拡散速度、波動関数の形を、それぞれ異なる方法で調べた点にある。個々の測定には直接示せる範囲と模型に依存する範囲がある。それでも、α相では等価な四面体サイトと高速な低温移動が対応し、β相では一軸ひずみ、八面体サイト占有、熱活性化型拡散が対応する。独立した証拠が同じ因果関係へ収束することで、結晶対称性が水素トンネルを左右するという結論が支えられている。
6. β相では熱活性化型の拡散が現れた
高水素濃度の β相では、水素の移動速度が温度に強く依存した。平均水素濃度が H/V = 0.29 の試料を 65 K 未満に保つと、12 時間測定しても有意な抵抗緩和は確認されなかった。表面側に偏っていた水素分布が、測定時間内には膜厚方向へほとんど均一化しなかったことを意味する。
温度を 65 K から 100 K の範囲へ上げると、抵抗は時間とともに増加し、温度が高いほど短時間で一定値へ近づいた。抵抗変化が速くなるのは、水素が表面側から膜の深部へ移動し、濃度分布が均一化するまでの時間が短くなったためである。温度上昇によって拡散が加速するという観測結果は、熱エネルギーを受け取った水素が格子間サイトの障壁を越える機構と整合する。
熱活性化型の拡散では、温度が上がると、移動障壁を越えられる水素原子核の割合が指数関数的に増える。拡散係数と温度の間に現れるこの関係をアレニウス則という。活性化エネルギーは、あるサイトから次のサイトへ移るために必要な代表的エネルギーを表す。値が大きいほど、同じ温度では移動できる水素の割合が少なくなり、拡散は遅くなる。
研究グループは、測定された抵抗緩和を、膜厚方向の水素移動を表す 1 次元の無作為歩行模型と比較した。無作為歩行模型では、水素が一定の時間間隔で隣接する層へ確率的に移ると仮定し、その移動を繰り返したときの濃度分布と電気抵抗を計算する。拡散係数を変えながら実測曲線を再現する値を求めることで、各温度における β相の拡散係数を推定した。
得られた拡散係数は 65 K から 100 K の範囲でアレニウス則に従い、活性化エネルギーは 148 ± 20 meV となった[1]。測定点は、一つの熱活性化過程を表す温度依存性の上に並んだ。低温側で拡散係数が一定値へ近づく挙動や、アレニウス則から高速側へ外れる変化は確認されなかった。
| 観測条件 | 測定結果 | 機構に対する意味 |
|---|---|---|
| β相、65 K 未満 | 12 時間測定しても有意な抵抗緩和が現れなかった。 | 低温で温度に依存せず残る高速なトンネル拡散は確認されなかった。 |
| β相、65 K から 100 K | 温度が高いほど抵抗緩和が速くなった。 | 水素が熱エネルギーを得て移動障壁を越える過程と整合する。 |
| β相の拡散係数 | アレニウス則に従い、活性化エネルギーは 148 ± 20 meV となった。 | 測定範囲では熱活性化型のサイト間移動が支配的である。 |
| α相、35 K | 測定開始前に抵抗緩和が完了するほど移動が速かった。 | 既知の熱拡散を低温へ外挿した予測では説明できない。 |
この結果から、β相で量子力学的な効果がすべて失われたとは言えない。水素原子核は β相でも量子状態を持ち、零点振動も存在する。実験から判定できるのは、65 K から 100 K における巨視的な水素再分布が、量子トンネルを主因とする温度非依存の拡散ではなく、熱エネルギーを必要とするサイト間移動として説明できることである。
低水素濃度の α相では、同じ測定系から対照的な結果が得られた。平均水素濃度が H/V = 0.03 の試料では、35 K で水素を注入した後の抵抗緩和を時間分解できなかった。測定を開始した時点で抵抗はすでに均一分布に対応する値へ達しており、水素は数秒以内に厚さ 40 nm の薄膜を移動したと判断された。
この高速移動が通常の熱拡散で説明できるかを確かめるため、先行研究で高温側から得られた α相の熱拡散係数を 35 K まで外挿した。外挿値に基づけば、水素が 100 秒間に移動する代表的な距離は約 0.1 nm となる。これはバナジウム原子間の距離よりも短い規模であり、膜厚 40 nm にわたる再分布を生じさせるには足りない。
実際には、水素分布は 100 秒より短い数秒のうちに膜全体へ広がった。移動時間は予測より短く、到達距離は熱活性化型拡散の外挿値より数百倍長い。単に拡散係数の数値が多少ずれたのではなく、低温へ向かうほど急速に遅くなるはずの熱拡散と、35 K でも測定限界を超えて速い実際の移動との間に、機構の違いを必要とする隔たりがある。
α相については、拡散が速すぎたため、β相と同じ方法で詳細な抵抗緩和曲線を取得できなかった。測定から得られるのは、数秒以内に膜厚 40 nm を移動したという時間と距離の上限、およびそこから求められる拡散係数の下限値である。温度ごとの拡散係数を並べ、低温で一定値へ近づく曲線を直接観測したわけではない。
α相の量子トンネルという判断は、この測定限界を含めて組み立てられている。実験は、35 K での実際の移動が熱活性化型拡散の外挿値を大幅に上回ることを示す。立方対称性を持つ α相の計算からは、複数の四面体サイトへ広がる量子状態が得られる。熱拡散では説明できない速度と、トンネルを許す波動関数の構造が一致することで、量子トンネルによる移動という解釈が支えられる。
β相では、温度を上げるほど移動が速くなり、拡散係数は一つの活性化エネルギーで記述できた。α相では、熱拡散ならほぼ停止するはずの 35 K で、測定時間より速く水素が移動した。両相の差は、同じ水素が温度だけによって異なる振る舞いを見せた結果ではない。水素濃度によって結晶構造と安定サイトが変わり、その構造差が、低温で利用できる移動機構を量子トンネルと熱活性化過程に分けた結果である。
7. 立方晶では波動関数が複数サイトへ広がる
実験では、低水素濃度の α相と高水素濃度の β相で、低温における水素の移動速度が大きく異なった。抵抗緩和と水素分布の測定からは速度差を確認できるが、その差が結晶対称性によって生じたのか、欠陥密度や界面、水素濃度に伴う別の効果によるのかまでは決まらない。研究グループは、結晶構造から水素原子核の量子状態までを計算し、立方対称性の有無が波動関数の広がりへどう反映されるかを調べた。
計算模型には、バナジウム原子 54 個と水素原子 1 個からなる周期構造を用いた。水素原子を 1 個に限定することで、水素同士の直接的な相互作用を避け、結晶ひずみだけが水素の安定位置と量子状態へ与える効果を取り出している。この模型は低濃度の α相に近い条件を表す一方、多数の水素を含む実際の β相を完全に再現するものではない。β相について計算から直接判定できるのは、一軸ひずみが単独でも水素の局在化を生じさせるかという点である[23]。
比較したのは、立方対称性を保つ構造と、1 本の結晶軸を伸ばした正方構造である。正方構造では、伸長方向の格子定数 c と、それに垂直な方向の格子定数 a の比を c/a = 1.10 とした。c/a = 1 なら 3 方向の長さが等しい立方構造であり、1.10 では 1 本の軸が面内方向より 10％長い。一軸ひずみを数値として与えることで、α相に近い対称性の高い構造と、β相に近い対称性の低い構造を同じ原子数で比較できる[23]。
7.1 結晶構造から水素のポテンシャルエネルギー面を求める
最初の計算では、水素を結晶内のさまざまな位置へ置き、それぞれの位置で系全体が持つエネルギーを密度汎関数理論によって求めた。水素の位置を x、y、z の 3 方向へ動かしながらエネルギーを計算すると、水素原子核が結晶中で感じる 3 次元のポテンシャルエネルギー面が得られる。
ポテンシャルエネルギー面の低い場所は、水素が安定して存在しやすい格子間サイトに対応する。隣接する低い場所の間にある高い領域は、サイト間を移動するときの障壁になる。立方構造では、対称性で対応する四面体サイトが同じ深さの谷になる。正方構造では、格子の伸長によって谷の深さが方向ごとに変わり、伸長軸方向の八面体サイトが選択的に低くなる。
ポテンシャルエネルギー面だけから分かるのは、水素を点状の粒子として置いたときの安定位置と移動障壁である。量子力学的な水素原子核は、谷の底の一点に静止せず、零点運動によって有限の空間へ広がる。谷が複数並び、各谷の量子状態が重なれば、原子核の波動関数は一つの格子間サイトを越えて広がる可能性がある。
7.2 シュレーディンガー方程式から固有状態を求める
研究グループは、計算した 3 次元のポテンシャルエネルギー面を使い、水素原子核についてシュレーディンガー方程式を数値的に解いた。この計算で得られる固有状態は、原子核が取りうる定常的な量子状態であり、それぞれに固有のエネルギーと波動関数が対応する。
波動関数の絶対値を 2 乗した量は、水素原子核を各位置で見つける確率密度を表す。確率密度が一つの格子間サイト周辺だけに集中していれば、その状態は局在している。複数のサイトに同程度の確率が分布していれば、水素原子核はそれらのサイトにまたがる非局在化状態を形成している。
固有状態のエネルギー差も、サイト間の結合を判断する材料になる。互いに独立した同じ形の谷があり、谷の間に波動関数の重なりがまったくなければ、各サイトに局在した状態は同じエネルギーを保つ。有限のトンネル結合が生じると、局在状態は複数サイトに広がる状態へ組み替えられ、もともと同じだったエネルギー準位が複数に分裂する。準位分裂は、サイト間に量子力学的な接続が存在することを示す。
7.3 立方構造では 12 個の四面体サイトが一つの状態網を作る
立方構造の計算では、最も低い 12 個の固有状態が、縮退度 1、2、6、2、1 を持つ 5 群に分かれた[23]。縮退度とは、同じエネルギーを持つ独立した固有状態の数である。たとえば縮退度 6 の群では、異なる形を持つ 6 個の波動関数が同じエネルギーに属する。
この 1、2、6、2、1 という並びは、単に 12 個の低い状態が偶然近いエネルギーへ集まった結果ではない。立方格子内にある 12 個の等価な四面体サイトを、有限のトンネル結合で接続した模型から得られる準位構造と一致する。第一段階では、立方対称性が 12 個のサイトを同じエネルギーに保つ。第二段階では、隣接サイト間で波動関数が重なることにより、12 個の局在状態が結晶対称性に従った複数の非局在化状態へ分裂する。
計算された波動関数も、この解釈を支持した。最も低い状態の確率密度は、特定の四面体サイトだけに集中せず、対称性で結ばれた複数のサイトへ分布していた。水素原子核が時間とともに一つのサイトから順番に飛び移った様子を示しているのではない。固有状態そのものが複数サイトにまたがり、水素原子核の位置を一つのサイトだけへ割り当てられない状態になっている。
この非局在化は、α相で観測された高速な低温移動と対応する。複数の四面体サイトが一つの量子状態によって接続されていれば、水素原子核は各サイト間の障壁を熱によって一つずつ越える必要がない。格子中に形成された等価なサイト網を通じて、量子トンネルが実際の拡散へ寄与できる。
7.4 一軸ひずみは低エネルギー状態を八面体サイトへ集める
c/a = 1.10 の正方構造では、低エネルギー状態の構成が変わった。立方構造で等価だった四面体サイトは、一軸ひずみによって異なる局所構造とポテンシャルエネルギーを持つようになる。伸長軸方向の八面体サイトが他の格子間サイトより深い安定点となり、最も低い状態の確率密度はその周辺へ集中した。
局在化は、直接にはポテンシャルの深さの差から生じる。八面体サイトが深い谷になれば、そこに形成される低エネルギー状態と、四面体サイトや別方向のサイトに形成される状態との間に準位差が生じる。準位がそろわなければ、各サイトの波動関数に空間的な重なりがあっても、一つの固有状態へ強く混成しない。
背後にあるのは、立方対称性によって保たれていた縮退の解除である。立方構造では、対称操作で移り合うサイトが同じエネルギーを持つため、複数サイトを同等に含む状態が成立する。一軸ひずみが入ると、伸長軸方向と面内方向の区別が生じ、対応する量子準位が分裂する。最低エネルギーの状態は、最も深くなった八面体サイトを中心に形成され、他のサイトへの広がりを失う。
| 計算上の観点 | 立方構造 | 一軸ひずみを持つ正方構造 | 拡散への帰結 |
|---|---|---|---|
| 模型の対称性 | x、y、z 方向が等価である。 | 伸長軸と面内方向が非等価になる。 | 移動先を同じ条件で接続できるかが変わる。 |
| 安定な格子間サイト | 複数の四面体サイトが同じ深さの安定点になる。 | 伸長軸方向の八面体サイトが選択的に安定化する。 | 水素が利用できる低エネルギー経路が組み替わる。 |
| 量子準位 | 12 個の等価なサイト間の結合を反映して 1、2、6、2、1 の群へ分裂する。 | サイト間の等価性が失われ、八面体サイトに対応する状態が低エネルギー側へ分かれる。 | 等価な準位間の混成が、非等価な準位間では弱くなる。 |
| 波動関数 | 複数の四面体サイトへ非局在化する。 | 特定の八面体サイト周辺へ局在する。 | 低温でサイト間トンネルが拡散へ寄与できる範囲が変わる。 |
| 支配的な移動機構 | 熱による障壁越えが困難な低温でも、トンネル経路が残る。 | 別サイトへ移るには、主として熱による励起が必要になる。 | α相の高速移動と β相の熱活性化型拡散の差を説明する。 |
7.5 計算結果は実験の速度差をどこまで説明するか
量子状態計算だけから、実際の薄膜中における拡散係数を直接求めたわけではない。計算で示されたのは、立方構造では水素原子核の低エネルギー状態が複数サイトへ広がり、一軸ひずみを持つ構造では特定サイトへ局在するという差である。拡散速度との接続には、共鳴核反応法による水素分布と、電気抵抗緩和による時間変化の測定が必要になる。
α相では、35 K において熱活性化型拡散から予測される距離を大幅に超えて水素が移動した。計算では、α相に対応する立方構造で、複数の四面体サイトにまたがる固有状態が得られた。実験が示す異常に速い移動と、計算が示す非局在化状態が、量子トンネルという同じ機構で接続される。
β相では、65 K から 100 K の拡散係数がアレニウス則に従った。計算では、一軸ひずみによって八面体サイトが安定化し、低エネルギー状態がその周囲へ局在した。実験で見られた熱を必要とする移動と、計算された非等価なサイト構造が対応する。
この対応は、立方対称性だけが拡散速度を完全に決めることを意味しない。実際の β相には多数の水素が存在し、水素間相互作用、占有サイトの秩序化、局所ひずみ、欠陥、界面の影響も加わる。水素原子 1 個の模型が示すのは、一軸ひずみだけでもサイトの等価性を失わせ、波動関数を局在させる方向へ作用することである。
β相で「量子性が消えた」と表現すると、この計算結果を誤って一般化することになる。β相の水素原子核もシュレーディンガー方程式に従い、局在した固有状態と零点運動を持つ。失われたのは量子力学そのものではなく、複数の格子間サイトを同じ低エネルギー状態で結ぶ条件である。
立方対称性を持つ α相では、等価な四面体サイト、低い移動障壁、短いサイト間距離、そろった量子準位が一つの移動網を作る。一軸ひずみを持つ β相では、八面体サイトが選択的に安定化し、経由するサイトとの準位差が生じ、その網が分断される。結晶対称性は水素原子核が量子系であるかを決めるのではなく、その量子状態が複数サイトへ広がり、巨視的な拡散として現れるかを決めている。
8. 実験で示したことと模型に依存すること
本研究の結論は、一つの測定によって量子トンネルを直接捉えた結果ではない。結晶構造、水素の格子間位置、深さ方向分布、抵抗の時間変化、ポテンシャルエネルギー面、原子核の波動関数を別々に調べ、それらの対応関係から移動機構を特定している。個々の結果が示す範囲を区別することで、結晶対称性が水素トンネルを左右するという結論が、どの事実に支えられ、どの部分で模型を必要とするかが見えてくる。
β相で水素が伸長軸方向の八面体サイトを占有したことは、チャネリング共鳴核反応法による位置測定から得られた。65 K から 100 K における抵抗緩和が温度とともに速くなったことも、実験で観測された事実である。一方、抵抗曲線から拡散係数を得るには、水素濃度分布と電気伝導を結ぶ模型が必要になる。α相の量子トンネルについては、35 K でも測定限界を超える高速移動が起きた事実と、立方構造で波動関数が複数の四面体サイトへ広がる計算結果を統合して判断している。
8.1 直接測定された量と、模型から求めた量を分ける
実験装置が出力する値と、論文が最終的に議論する物理量は必ずしも同じではない。共鳴核反応法が直接測るのは、窒素イオンの入射条件に対する γ線の強度である。水素の深さ方向分布は、イオンが試料内で失うエネルギーと核反応の共鳴条件を使って、その強度を位置情報へ変換した結果である。イオンチャネリングによるサイト占有も、入射角度ごとの反応強度を、候補となる格子間サイトの計算結果と比較して判定する。
電気抵抗測定が直接記録するのは、試料全体の抵抗が時間とともにどう変わったかである。そこから水素の拡散係数を求めるには、膜厚方向の濃度分布を仮定し、各層の水素濃度を局所抵抗へ変換し、薄膜全体の電気伝導を計算しなければならない。拡散係数は抵抗計から直接読み取った値ではなく、実測曲線を最もよく再現する模型の変数として決められる。量子状態計算も、薄膜内の水素原子核を直接観察した結果ではない。密度汎関数理論から得たポテンシャルエネルギー面を使い、特定の結晶構造を仮定してシュレーディンガー方程式を解いた結果であり、実験で観測された α相と β相の拡散挙動を説明する機構として使われる。
8.2 1 次元拡散模型は何を仮定しているか
β相の拡散係数を求める模型では、厚さ 40 nm の薄膜を複数の層へ分割し、水素が隣接層へ確率的に移る過程を計算している。初期状態では水素が表面側へ偏り、時間の経過とともに膜厚方向へ均一化すると仮定する。各層の水素濃度から局所的な電気抵抗を求め、それらの層が膜面に沿った電流に対して並列に接続されるとして、試料全体の抵抗を計算する。
この模型では、膜面内の水素濃度は一様であり、主要な濃度勾配は膜厚方向にだけ存在するとみなす。実際の薄膜に面内方向の濃度むらや欠陥に沿った高速経路があれば、抵抗緩和には膜厚方向以外の移動も含まれる可能性がある。測定された抵抗曲線だけから、面内拡散と膜厚方向の拡散を独立に分離したわけではない。
公開査読記録では、模型を構成する層数への依存性が検討された[24]。薄膜を 40 層に分割した場合と 120 層に分割した場合を比較しても、同じ拡散係数に対する抵抗緩和曲線はほぼ変わらなかった。この結果は、数値計算の空間分割が粗すぎるために拡散係数が決まったのではないことを示す。
一方で、層数を増やしても 1 次元という前提自体は変わらない。空間分割に対する数値的な安定性と、実際の拡散が 1 次元で近似できるかという物理的な妥当性は別の論点である。前者は確認されたが、後者には薄膜面内の均一性という仮定が残る。
抵抗から水素濃度を求める係数についても、温度依存性が留保になる。水素濃度と抵抗の関係は一定温度で評価され、拡散測定を行った温度範囲でも大きく変わらないと仮定された。係数が温度によって変化すれば、抵抗緩和の一部を水素移動ではなく電気伝導特性の変化が生んだ可能性が生じる。著者らは既知の温度依存性と実測結果を基に影響が限定的であると説明しているが、拡散係数にはこの変換過程に由来する不確かさが含まれる。
8.3 イオンビームによる加熱はどの測定へ影響するか
共鳴核反応法では、高エネルギーの窒素イオンを薄膜へ照射する。イオンが試料中でエネルギーを失うと、その一部は熱へ変わるため、照射位置の温度は試料台で測った温度より高くなる可能性がある。水素拡散は温度に敏感であり、数 K の差でも低温での分布変化へ影響しうる。
公開査読記録では、20 秒間の照射による温度上昇が 10 K 程度の規模になりうると見積もられた[24]。この値はすべての測定条件に共通する厳密な補正値ではなく、ビーム電流、照射時間、熱伝導、照射面積に依存する推定である。照射点の実温度をその場で直接測っていないため、共鳴核反応測定で示された温度境界には不確かさが残る。
この留保が影響するのは、主としてイオンビームを使用した水素分布とサイト占有の測定である。電気抵抗緩和による拡散速度の測定は、窒素イオンを照射し続けながら実施したものではない。抵抗測定で得られた 65 K から 100 K の温度依存性へ、ビーム加熱の見積もりをそのまま適用することはできない。
測定ごとの条件を分けると、ビーム加熱が研究全体を無効にするという結論にはならない。β相で八面体サイトが占有されるという位置情報には温度評価の幅が生じるが、抵抗緩和がアレニウス則に従うことは別の測定系から得られている。複数の手法を使った研究では、一つの装置に固有の誤差がどの主張まで伝わるかを限定して読む必要がある。
8.4 水素原子 1 個の模型で高濃度 β相をどこまで説明できるか
量子状態計算では、バナジウム原子 54 個に対して水素原子を 1 個だけ置いている。この模型は、水素同士の直接的な相互作用を避け、結晶ひずみが単独で水素の安定位置と波動関数へ与える影響を調べるには適している。低濃度の α相を表す模型としては、実験条件との対応も比較的明確である。
高水素濃度の β相では、実際には多数の水素が格子間サイトを占有する。隣接する水素同士の反発、占有サイトの秩序化、局所的な格子変形、利用可能な移動先の減少が生じる。水素がすでに占有しているサイトは移動先として使えず、計算上は開いている経路が実材料では部分的に遮られる可能性もある。
水素原子 1 個の模型は、これらの高濃度効果を再現しない。計算結果から直接言えるのは、一軸ひずみだけを与えた場合でも、伸長軸方向の八面体サイトが安定化し、波動関数が局在するということである。実際の β相で観測された拡散の遅さを、単一水素模型だけから定量的に再現したわけではない。
査読への回答で著者らは、多数の水素を含む計算が実際の β相をより忠実に表すことを認めている[24]。そのうえで、一軸ひずみ模型が、実験で確認された八面体サイト占有と、先行研究で知られた拡散障壁の増加を再現することを示した。水素間相互作用や経路遮蔽は、局在化を弱めるよりも、サイトの非等価性を増し、移動可能な経路をさらに制限する方向へ働く可能性がある。
この対応から、結晶対称性の低下を中心機構として位置づけることはできる。ただし、β相の拡散係数、相安定性、温度依存性のすべてを対称性だけへ還元することはできない。対称性低下は量子トンネルを抑える必要な構造変化を与えるが、高濃度相で観測される数値には、水素間相互作用、欠陥、基板拘束、局所ひずみも加わる。
8.5 α相の量子トンネルはどこまで直接示されたか
α相では、35 K における抵抗緩和が速すぎて時間分解できなかった。実験から確認できたのは、水素が数秒以内に厚さ 40 nm の薄膜内へ広がったことと、その速度が既知の熱拡散係数を低温へ外挿した予測を大幅に上回ったことである。拡散係数の詳細な温度依存性や、低温で一定値へ近づく曲線を測定したわけではない。
量子トンネルという解釈は、二つの独立した差から組み立てられている。実験側では、熱活性化型の拡散では到達できない距離を低温で移動した。計算側では、立方構造の複数の四面体サイトへ広がる低エネルギー状態が得られた。異常な速度だけなら、欠陥や界面に沿った高速拡散という別の説明も考えられる。非局在化状態の計算だけなら、実際の試料でその経路が使われたことを証明できない。両者が対応することで、量子トンネルが最も整合的な機構となる。
同位体を変えて拡散速度の質量依存性を詳細に比較する測定や、α相の高速緩和を捉えられる時間分解能を持つ測定が加われば、量子トンネルの実験的な直接性はさらに高まる。現段階では、熱拡散からの大幅な逸脱と量子状態計算の一致が、α相のトンネル拡散を支える中心的な証拠である。
| 主張 | 直接測定された事実 | 模型または計算による接続 | 残る留保 |
|---|---|---|---|
| β相では八面体サイトが占有される | チャネリング条件に対する核反応強度の角度依存性を測定した。 | 八面体サイトと四面体サイトを仮定した計算曲線との比較から占有位置を判定した。 | イオンビームによる照射点の温度上昇に不確かさがある。 |
| β相の拡散は熱活性化型である | 65 K から 100 K で温度上昇とともに抵抗緩和が速くなった。 | 1 次元無作為歩行模型から得た拡散係数がアレニウス則に従った。 | 膜面内を一様とする近似と、濃度から抵抗への変換係数を含む。 |
| α相では低温でも高速移動が起こる | 35 K で抵抗緩和を時間分解できず、水素分布が数秒以内に均一化した。 | 既知の熱拡散係数を低温へ外挿した移動距離と比較した。 | 拡散係数の詳細な温度依存性は測定できず、下限値だけが得られた。 |
| α相の高速移動は量子トンネルによる | 熱拡散では説明できない移動距離が観測された。 | 立方構造で複数の四面体サイトへ広がる固有状態が得られた。 | 移動中の波動関数を直接観測した結果ではない。 |
| 対称性低下がトンネルを抑える | 高濃度相で一軸ひずみ、八面体サイト占有、熱活性化型拡散が対応した。 | 一軸ひずみ模型でサイトの非等価化と波動関数の局在化が得られた。 | 高濃度相の水素間相互作用と経路遮蔽を完全には計算していない。 |
各留保は、すべての結論へ同じ強さで影響するわけではない。1 次元拡散模型の仮定は β相の拡散係数の定量値に関係するが、抵抗緩和が温度とともに速くなった観測事実までは消さない。ビーム加熱は共鳴核反応測定の温度評価に影響するが、別に測定された抵抗緩和のアレニウス型温度依存性を説明しない。単一水素模型は高濃度 β相の多体効果を省くが、一軸ひずみだけで八面体サイトの安定化と局在化が起こるという機構は示している。
証拠の強さは、結論ごとに異なる。β相の八面体サイト占有と熱活性化型の抵抗緩和は、比較的直接性の高い実験結果である。α相の量子トンネルは、熱拡散では説明できない高速移動と非局在化状態の計算を統合した機構解釈である。結晶対称性の低下がトンネルを抑えるという中心命題は、これらの位置、速度、構造、量子状態が同じ相変化に対応して変わることで支えられる。
外部ひずみによって結晶対称性を変え、量子トンネルを可逆に切り替える操作はまだ実証されていない。水素貯蔵容量、吸放出速度、動作温度、耐久性が改善された結果も得られていない。本研究が与えたのは、完成した水素制御技術ではなく、格子間サイトの等価性と量子準位の整合を材料設計の対象にできるという機構上の指針である。
9. 量子トンネルは結晶が許した経路として現れる
α相と β相で変わるのは、水素原子核の質量ではない。どちらの相でも、通常の軽水素であれば原子核は陽子 1 個からなり、量子力学に従う。変化するのは、水素濃度によって形成されるバナジウム格子の構造と、その格子が水素原子核へ与えるポテンシャルエネルギー面である。
9.1 水素の軽さは可能性を与え、結晶構造が発現条件を決める
「水素は軽いから量子トンネルする」という説明は、トンネルが起こりうる理由を示すが、α相と β相の差までは説明しない。水素原子核の小さな質量は、障壁を透過できる波動関数の広がりを与える。しかし、その広がりが隣接サイトへ届き、実際の輸送へ結びつくには、結晶側にも条件が必要になる。
本研究でトンネル拡散を支えたのは、短いサイト間距離、十分に低く狭い障壁、等価な移動先、整合した量子準位である。これらの条件は独立して働くのではない。結晶対称性がサイトを等価に保ち、等価なサイトが準位をそろえ、短い距離と有限の障壁が波動関数の重なりを作る。複数の条件が同時に成立したとき、局在状態は非局在化状態へ変わり、トンネルが拡散として観測される。
β相では、水素原子核の軽さそのものは残っている。それでも、結晶の一軸ひずみによってサイト網が非等価な状態へ組み替えられるため、軽さが高速拡散として発現する経路が閉じる。能力の有無ではなく、能力を実際の運動へ変換できる構造条件の有無が、両相を分けている。
9.2 水素濃度は移動する粒子の数だけを表していない
本研究では、水素濃度が拡散機構を変える。しかし、高濃度になるほど水素同士が混雑し、移動が遅くなったという説明だけでは不十分である。濃度の上昇は、まずバナジウム格子を一軸方向へ変形させる。その変形が安定な格子間サイトを変え、サイト間のエネルギー差を生み、低エネルギーの移動経路を組み替える。
水素は、あらかじめ固定された結晶内を受動的に移動しているわけではない。水素の蓄積が母体格子を変え、変化した格子が後から移動する水素の条件を変える。移動主体の量が、移動環境そのものを更新している。
この循環があるため、拡散係数を水素とバナジウムの組み合わせに固有の一つの値として扱うことはできない。同じ元素の組み合わせでも、水素濃度、結晶相、基板拘束、ひずみ方向が変われば、安定サイトと移動機構が変わる。材料の輸送特性は、組成だけでなく、その組成が作る構造状態に依存する。
9.3 量子的な機能は構成要素単独では説明できない
既稿では、量子計算機が壊れやすい物理量子ビットをそのまま論理情報として扱うのではなく、多数の量子ビットが作る制約と関係の中へ情報を置く構造を扱った[25]。量子誤り訂正と金属中水素のトンネルは対象も機構も異なるが、構成要素が量子的であるだけでは観測可能な機能が成立せず、要素間の関係が必要になる点に限れば共通する。本研究では、その関係がバナジウム格子の対称性、格子間サイト、移動障壁、量子準位として具体化されている。
9.4 材料設計ではポテンシャル地形全体を扱う必要がある
水素輸送を制御する材料設計では、元素の種類や平均的な活性化エネルギーだけを比較しても、低温で量子トンネルが現れる条件を十分に捉えられない。α相と β相の比較が示したのは、同じ元素組成でも、結晶対称性が変われば安定サイト、移動距離、障壁、量子準位、波動関数が連動して変わることである。
設計対象には、結晶軸の比、格子間サイトの等価性、サイト間距離、移動経路上の障壁、始点と終点の準位整合を含める必要がある。薄膜では、基板が膜面内の膨張を拘束するため、水素吸収によるひずみが特定方向へ集中する。基板材料、膜厚、界面、外部応力を変えれば、同じバナジウム水素系でも異なるポテンシャル地形を形成できる可能性がある。
ただし、本研究は外部応力を操作してトンネルを可逆に切り替えた実験ではない。水素貯蔵容量、吸放出速度、耐久性、常温での動作性能を改善した結果も示していない。得られたのは、結晶対称性の低下が量子トンネルを抑える因果経路と、その因果を利用して移動機構を設計できる可能性である。
9.5 結晶が変わると、同じ水素の移動原理が変わる
α相では、立方対称性が複数の四面体サイトを等価に保ち、低い障壁とそろった準位によって波動関数を非局在化させる。β相では、水素濃度の上昇が格子を一軸方向へひずませ、八面体サイトを選択的に安定化し、経由地との準位差によって波動関数を局在させる。前者では量子トンネルが低温拡散へ寄与し、後者では熱活性化型のサイト間移動が支配的になる。
この差から導かれる帰結は、水素原子核が「波」から「粒子」へ変化したということではない。同じ量子力学に従う水素原子核について、結晶が複数サイトへ広がれる状態を許すか、特定サイトへ閉じ込める状態を作るかが変わったのである。
水素原子核の軽さは、量子トンネルの可能性を与える。バナジウム結晶の対称性は、その可能性を低温の高速拡散として発現させる経路を作る。水素が増えて格子がひずむと、等価なサイト網と準位整合が失われ、その経路は閉じる。量子トンネルを決めていたのは粒子の性質だけではなく、粒子が移動できる状態空間を作る結晶構造だった。
参考文献
- Sudhansu Sekhar Das, Takahiro Ozawa, Taizo Kawauchi, Hiroshi Nakanishi and Katsuyuki Fukutani, Impact of crystal symmetry lowering on proton tunneling, Nature Communications, Vol. 17, Article 5928（2026）. https://doi.org/10.1038/s41467-026-75020-w
- 東京大学生産技術研究所, 水素はいつ「粒子」になり、いつ「波」になるのか？――結晶対称性が量子トンネル現象を支配することを解明――（2026-07-16）. https://www.iis.u-tokyo.ac.jp/ja/news/5104/
- id774, シュレーディンガー方程式と量子測定問題（2026-04-30）. https://blog.id774.net/entry/2026/04/30/4644/
- Anke Pundt and Reiner Kirchheim, Hydrogen in Metals: Microstructural Aspects, Annual Review of Materials Research, Vol. 36, pp. 555–608（2006）. https://doi.org/10.1146/annurev.matsci.36.090804.094451
- Yuh Fukai, Evidence for quantum diffusion of hydrogen in Ta: quench-recovery experiments revisited, Japanese Journal of Applied Physics, Vol. 23, L596（1984）. https://doi.org/10.1143/JJAP.23.L596
- Jun Kondo, Diffusion of light interstitials in metals, Physica B+C, Vol. 126, pp. 377–384（1984）. https://doi.org/10.1016/0378-4363(84)90191-8
- Gerhard Sicking, Isotope effects in metal-hydrogen systems, Journal of the Less-Common Metals, Vol. 101, pp. 169–190（1984）. https://doi.org/10.1016/0022-5088(84)90093-6
- Takahiro Ozawa et al., Observation of resonant tunneling of proton from octahedral to tetrahedral sites in Pd, Journal of Physics and Chemistry of Solids, Vol. 185, Article 111741（2024）. https://doi.org/10.1016/j.jpcs.2023.111741
- Takahiro Ozawa, Ryota Shimizu, Isao Harayama, Taro Hitosugi, Daiichiro Sekiba and Katsuyuki Fukutani, Observation of proton tunneling correlated with phonons and electrons in Pd, Science Advances, Vol. 11, eady8495（2025）. https://doi.org/10.1126/sciadv.ady8495
- Xue Xin, Robert Johansson, Martin Wolff and Björgvin Hjörvarsson, Hydrogen in vanadium: site occupancy and isotope effects, Physical Review B, Vol. 93, Article 134107（2016）. https://doi.org/10.1103/PhysRevB.93.134107
- Wolfgang Pesch, Thomas Schober and Helmut Wenzl, A TEM study of the phase diagrams VH and VD, Scripta Metallurgica, Vol. 16, pp. 307–312（1982）. https://doi.org/10.1016/0036-9748(82)90358-1
- V. A. Somenkov and S. S. Shil’stein, Phase transitions of hydrogen in metals, Progress in Materials Science, Vol. 24, pp. 267–335（1980）. https://doi.org/10.1016/0079-6425(79)90006-9
- L. P. A. Mooij et al., Influence of site occupancy on diffusion of hydrogen in vanadium, Physical Review B, Vol. 95, Article 064310（2017）. https://doi.org/10.1103/PhysRevB.95.064310
- Volker Lottner, Andreas Heim and Thomas Springer, The Debye-Waller factor and the diffusion process for hydrogen in Nb, Ta, and V single crystals investigated by neutron spectroscopy, Zeitschrift für Physik B, Vol. 32, pp. 157–165（1979）. https://doi.org/10.1007/BF01320112
- Yuh Fukai and Shinji Kazama, NMR studies of anomalous diffusion of hydrogen and phase transition in vanadium-hydrogen alloys, Acta Metallurgica, Vol. 25, pp. 59–70（1977）. https://doi.org/10.1016/0001-6160(77)90246-2
- Robert Johansson, Rajeev Ahuja, Olle Eriksson, Björgvin Hjörvarsson and Richard H. Scheicher, Effect of uniaxial strain on the site occupancy of hydrogen in vanadium from density-functional calculations, Scientific Reports, Vol. 5, Article 10301（2015）. https://doi.org/10.1038/srep10301
- Wen Huang et al., Concentration dependence of hydrogen diffusion in clamped vanadium (001) films, Journal of Physics: Condensed Matter, Vol. 29, Article 045402（2017）. https://doi.org/10.1088/1361-648X/29/4/045402
- Jian Luo et al., Dissolution, diffusion and permeation behavior of hydrogen in vanadium: a first-principles investigation, Journal of Physics: Condensed Matter, Vol. 23, Article 135501（2011）. https://doi.org/10.1088/0953-8984/23/13/135501
- Jian Qin et al., Investigation of adsorption, dissociation, and diffusion properties of hydrogen on the V(100) surface and in the bulk: a first-principles calculation, Journal of Advanced Research, Vol. 21, pp. 25–34（2020）. https://doi.org/10.1016/j.jare.2019.09.003
- Kazuhiro Asano, Satoshi Hayashi, Yasushi Nakamura and Etsuo Akiba, Effect of substitutional Mo on diffusion and site occupation of hydrogen in the BCT monohydride phase of V–H system studied by 1H NMR, Journal of Alloys and Compounds, Vol. 507, pp. 399–404（2010）. https://doi.org/10.1016/j.jallcom.2010.08.018
- J. E. Kleiner, E. H. Sevilla and R. M. Cotts, Diffusion of hydrogen in α′-VHx, Physical Review B, Vol. 33, pp. 6662–6666（1986）. https://doi.org/10.1103/PhysRevB.33.6662
- Takahiro Ozawa et al., Isotope-dependent site occupation of hydrogen in epitaxial titanium hydride nanofilms, Nature Communications, Vol. 15, Article 9558（2024）. https://doi.org/10.1038/s41467-024-53838-6
- Sudhansu Sekhar Das et al., Supplementary Information for Impact of crystal symmetry lowering on proton tunneling（2026）. https://static-content.springer.com/esm/art%3A10.1038%2Fs41467-026-75020-w/MediaObjects/41467_2026_75020_MOESM1_ESM.pdf
- Nature Communications, Transparent Peer Review file for Impact of crystal symmetry lowering on proton tunneling（2026）. https://static-content.springer.com/esm/art%3A10.1038%2Fs41467-026-75020-w/MediaObjects/41467_2026_75020_MOESM2_ESM.pdf
- id774, 量子計算機はなぜ壊れやすい情報で計算できるのか（2026-06-11）. https://blog.id774.net/entry/2026/06/11/4874/