Determination of nuclear quadrupole moments for $^{25}$Mg, $^{87}$Sr, and $^{135,137}$Ba via configuration-interaction combined with a coupled-cluster approach

本論文は、構成間相互作用および結合クラスター法を用いたアプローチにより、Mg、Sr、およびBaの低励起状態における電場勾配および磁気双極子超微細構造定数を計算しており、これにより$^{25}$Mg、$^{87}$Sr、および$^{135,137}$Baの核四重極モーメントの正確な決定を可能にし、それらはストロンチウムおよびバリウムに関して以前に採用されていた値との著しい相違を明らかにしている。

要約

原子核を、完璧な球体ではなく、少し押しつぶされたり引き伸ばされたりした風船のようなものだと想像してみてください。この形はランダムなものではなく、内部の陽子と中性子がどのように配置されているかを示す特定の「署名」なのです。科学者たちは、この形状を核四重極モーメントと呼んでいます。この「風船」の正確な形を知ることは、原子がどのように結合するかから材料の振る舞いに至るまで、物理学の根本的なルールを理解する上で極めて重要です。

しかし、この「つぶれ具合」を直接測定することは非常に困難です。それは、暗くて密閉された箱の中にある風船の正確な形を、箱を叩いた時に聞こえる音だけで推測しようとするようなものです。

実験：原子の「ハミング」を聞く

この論文の中で、著者であるYong-Bo Tang氏は、熟練のオーディオエンジニアとして振る舞っています。彼はマグネシウム（Mg）、ストロンチウム（Sr）、バリウム（Ba）という3つの特定の「原子の家族」に焦点を当てています。

これらの原子が励起されると（ギターの弦を弾くときのように）、超微細構造と呼ばれる非常に特定の「ハミング」や振動を放出します。このハミングは、外側で回転する電子と、内側にある押しつぶされた原子核との相互作用によって引き起こされます。

• 測定された部分： 科学者たちは、実験室でこのハミングのピッチ（音高）を非常に精密に測定しています。
• 失われた環（ミッシングリンク）： そのピッチから原子核の形状（四重極モーメント）を導き出すには、電子が核の周囲にどのように配置されているかを正確に知る必要があります。この配置は「電場勾配」（電子が転がり落ちる丘の傾斜のようなもの）を作り出します。

問題：計算するには丘が急すぎる

この「丘」の形状（電場勾配）を計算することは、コンピュータにとって悪夢です。電子はただ静止しているわけではありません。彼らは「電子相関」と呼ばれる複雑な方法で、互いに押し合い、引き合いながら踊っています。

• もしこれらの「ダンス」を無視すれば、丘の計算は間違ったものになります。
• もしすべてのダンスを計算しようとすれば、コンピュータはクラッシュしてしまいます。

これまでの試みは、ぼやけた衛星写真を使って山脈を地図に描き起こそうとするようなものでした。その結果は一貫していませんでした。ストロンチウムやバリウムについては、異なる研究によって結果が異なり、最大で10%もの差が出ていました。

解決策：ハイブリッドな「スイスアーミーナイフ」的アプローチ

これを解決するために、Tang氏は2つの強力な手法を組み合わせた新しい計算手法を開発しました。

1. 構成間相互作用法 (CI): これは、電子がどのように配置され得るかのあらゆる可能性を、一つずつ丁寧に見ていくようなものです。徹底していますが、時間がかかります。
2. 結合クラスター法 (CC): これは、電子がグループとしてどのように影響し合うかを予測する、洗練されたショートカット（近道）を使うようなものです。高速ですが、時として細部を見落とすことがあります。

Tang氏の手法であるCI+CCは、これら両方の良いとこ取りをしたものです。コア電子間の大きく重い相互作用を処理するために「ショートカット」を使用し、外殻電子の詳細を微調整するために「徹底的な」手法を使用します。これは、ドローンを使って森林の全体的な形をマッピングした後、特定の木一本一本の正確な高さを測るためにハイカーのチームを派遣するようなものです。

結果：混乱を晴らす

この高精度な「スイスアーミーナイフ」を用いて、Tang氏はMg、Sr、Baのいくつかの低エネルギー状態における電場勾配を計算しました。そして、これらの計算を既知の実験による「ハミング」と組み合わせることで、原子核の形状を決定しました。

得られた結果は以下の通りです：

• マグネシウム (25Mg): 結果はこれまでの実験と完璧に一致しました。ラジオのチューニングをして、放送局が極めてクリアに聞こえる状態を見つけたようなものです。計算された形状は、「ミューオンX線」実験（別のハイテクな測定方法）で発見されたものと一致しています。
• ストロンチウム (87Sr): ここで事態は複雑になります。Tang氏の結果は、原子核が現在教科書に記載されている値よりも約10%ほど押しつぶされていることを示唆しています。古い値は、ストロンチウムの「イオン」（電子を一つ失った原子）を観察したことから得られたものですが、Tang氏は「中性原子」を観察しました。この差は、従来の計算が何らかの微妙な電子のダンスを見落としていた可能性を示唆しています。
• バリウム (135,137Ba): ストロンチウムと同様に、Tang氏のバリウムに関する結果は、現在受け入れられている値（バリウムイオンから導出されたもの）と約4%の差があります。

まとめ

論文は、この手法がマグネシウムに対しては見事に機能する一方で、ストロンチウムやバリウムについては、科学者が現在使用している「ゴールドスタンダード（標準指標）」と比較して大きな不一致があることを結論づけています。

Tang氏は、その違いの原因として、現在の「ゴールドスタンダード」の計算が、三重励起（3つの電子が同時に相互作用すること）と呼ばれる特定の種類の電子相互作用を見逃していた可能性があると示唆しています。合唱団において、誰も予想していなかったハーモニーを3人の歌手が奏でれば曲の響きが変わるように、これらの三重の相互作用が、原子の形状の「ピッチ」を変化させているのかもしれません。

要約すると： 著者は、原子核の形状を測定するためのより優れたコンピュータモデルを構築しました。マグネシウムについては、そのモデルは既知の事実を裏付けました。ストロンチウムとバリウムについては、彼のモデルは現在の「公式」の測定値がわずかにずれている可能性を示唆しており、これらの原子核の真の形状を知るためには、3つの電子がどのように相互作用するかをもっと詳しく見る必要があることを示しています。

技術概要

技術要約：$^{25}\text{Mg}$、$^{87}\text{Sr}$、および $^{135,137}\text{Ba}$ の核四重極モーメントの決定

問題提起
核電気四重極モーメント（$Q$）は、核構造の対称性を記述する基本的なパラメータであり、原子、核、および凝縮系物理学において重要な意味を持つ。実験的な超微細構造定数（$B$）は高精度で測定可能であるが、正確な$Q$の値を抽出するには、電気電場勾配（$q$）の同等に精密な理論計算が必要となる。既存の文献では、$^{87}\text{Sr}$ および $^{137}\text{Ba}$ の推奨される核四重極モーメントの値に関して、特に中性原子とイオンから導出された値を比較した際に、顕著な不一致が見られる。例えば、中性Sr原子から導出された値は、イオンのデータから推奨される値よりも約8%大きく、Baについても同様の不整合が存在する。さらに、標準的な理論的手法は、重いアルカリ土類系における高精度な勾配計算に不可欠な、複雑な電子相関効果を説明することにしばしば苦慮する。

手法
著者らは、構成間相互作用（CI）法と結合クラスター（CC）法を組み合わせた、RCI+all-order法と呼ばれる相対論的ハイブリッド・アプローチを採用している。このフレームワークは、正確な電気電場勾配計算において決定的な役割を果たす、コア－コア、コア－価電子、および価電子－価電子の電子相関を包括的に説明するように設計されている。

主要な手法構成要素は以下の通りである：

• ハミルトニアンと相互作用： 超微細相互作用は、磁気双極子（$k=1$）および電気四重極（$k=2$）相互作用のための球面テンソル演算子を用いてモデル化されている。原子核は、一様に磁化され、かつ帯電した球体として扱われる。
• ハイブリッド計算スキーム： ディラック波動関数を展開するために、有限の基底関数セットとして、等間隔のガウス型関数を用いている。相関ポテンシャルは、不確実性を評価するために5つの異なる計算戦略を通じて導出される：
  1. CI + 二次多体摂動論（MBPT(2)）
  2. CI + 線形化結合クラスター単一および二重励起（LCCSD）
  3. CI + 完全結合クラスター単一および二重励起（CCSD）
  4. 再スケーリングパラメータを用いたCI + LCCSD（CI+LCCSDs）
  5. 再スケーリングパラメータを用いたCI + CCSD（CI+CCSDs）
• 補正： 遷移行列要素には、ランダム位相近似（RPA）、コア・ブリューケン効果、構造放射、および規格化補正（総称して高次（HO）補正と呼ばれる）に加え、二次までの二粒子（TP）相互作用が含まれる。
• $Q$ の抽出： 核四重極モーメントは、$Q = B / (234.9648867 \cdot q)$ という関係式を用いて抽出される。ここで、$B$ は実験的に測定された超微細構造定数であり、$q$ は理論的に計算された電気電場勾配である。

主な貢献と結果
本研究は、中性 $^{25}\text{Mg}$、$^{87}\text{Sr}$、および $^{135,137}\text{Ba}$ の低励起状態に対する電気電場勾配および磁気双極子超微細構造定数（$A$）を系統的に計算している。解析対象の特定の状態は以下の通りである：

• Mg: $3s3p \ ^3P_{1,2}$
• Sr: $5s4d \ ^1D_2$ および $5s5p \ ^3P_{1,2}$
• Ba: $6s5d \ ^3D_{1,2,3}$、$6s5d \ ^1D_2$、および $6s6p \ ^1P_1$

エネルギーおよび超微細定数：
計算されたエネルギーおよび磁気双極子定数（$A$）は、実験データ（NIST）および他の理論結果と極めて良好な一致を示している。CI+LCCSD法が最も信頼できる手法として特定されており、最終的な推奨値としての $A$ は、その不確実性の範囲内で（ほとんどの状態において通常2%以内）実験値と一致している。これは、アルカリ土類系に対するハイブリッドCI+CCアプローチの妥当性を検証するものである。

核四重極モーメント（$Q$）：
計算された勾配と実験的な $B$ 値を組み合わせることで、著者らは以下の核四重極モーメントを決定した：

• $^{25}\text{Mg}$: $0.203(2)$ b。この結果は、メソンX線実験による値（$0.201(3)$ b）および従来の超微細構造に基づく決定値と完全に一致している。
• $^{87}\text{Sr}$: $0.336(4)$ b。この値は、現在採用されている $0.305(2)$ b（Sr$^+$ の $4d_{5/2}$ 状態から導出）よりも約10%大きいが、Sr$^+$ の $5p_{3/2}$ 状態および中性Srの状態（$5s5p \ ^3P_{1,2}$）から導出された値とはよく一致している。
• $^{137}\text{Ba}$: $0.246(4)$ b。これは、推奨値である $0.236(3)$ b（Ba$^+$ の $5d_{3/2,5/2}$ 状態から導出）よりも約4%大きい。しかし、Ba$^+$ の $6p_{3/2}$ 状態から導出された値とは極めて良く一致している。
• $^{135}\text{Ba}$: $0.161(2)$ b.

意義と主張
本論文は、構成間相互作用（CI）と結合クラスター（CC）を組み合わせたアプローチが、特に直接測定が困難な不安定または重い原子核に対して、強固で正確な枠組みを提供することを主張している。

著者らは、SrおよびBaにおいて、中性原子から導出された値とイオンから導出された値の間の重大な不一致を強調している。具体的には、本研究における $^{87}\text{Sr}$ および $^{137}\text{Ba}$ の中性原子から導出された結果は、文献で現在採用されている値（イオンのデータに基づいているもの）に対し、それぞれ約10%および約4%異なっている。著者らは、これらの不一致は、標準的な結合クラスター計算における三重励起の寄与の欠落に起因している可能性があると示唆しており、近年のイオンに関する研究では、三重励起が超微細定数に1〜2%寄与することが示されていることを指摘している。

本研究は、$^{25}\text{Mg}$ に関する結果は確立されたデータと一致しているものの、$^{87}\text{Sr}$ および $^{137}\text{Ba}$ に関する知見は、現在採用されているイオン状態からの値を覆すものであると結論付けている。著者らは、磁気双極子定数の計算における高い精度によって裏付けられた、彼らの中性原子ベースの抽出が信頼できるものであると断言しているが、中性原子とイオンの決定の間の差異を完全に解決するためには、三重励起や三体相互作用を含むさらなる理論的研究が必要であることを認めている。