原子衝突学会誌 2012 年第 9 巻第 5 号
Journal of atomic collision research, vol. 9, issue 5, 2012.

しょうとつ

原子衝突学会 2012 年 9 月 15 日発行
http://www.atomiccollision.jp/

原子衝突学会賛助会員(五十音順)
アイオーピー・パブリッシング・リミテッド（IOP英国物理学会出版局）

http://journals.iop.org/

アステック株式会社

http://www.astechcorp.co.jp/

アドキャップバキュームテクノロジー株式会

http://www.adcap-vacuum.com

有限会社

イーオーアール

http://www.eor.jp/

Electronics Optics Research Ltd.
株式会社 オプティマ

http://www.optimacorp.co.jp/

カクタス・コミュニケーションズ株式会社

http://www.editage.jp
http://www.cactus.co.jp

キャンベラジャパン株式会社

http://www.canberra.com/jp/

クリムゾン インタラクティブ プライベート リミテッド

株式会社 サイエンス ラボラトリーズ

http://www.enago.jp/
http://ulatus.jp/
http://www.voxtab.jp /http://www.scilab.co.jp/

1真空光学株式会社

http://www.shinku-kogaku.co.jp/

スペクトラ・フィジックス株式会社

http://www.spectra-physics.jp/

ソーラボジャパン株式会社

http://www.thorlabs.jp/

ツジ電子株式会社

http://www.tsujicon.jp/

株式会社東京インスツルメンツ

http://www.tokyoinst.co.jp/

株式会社東和計測

http://www.touwakeisoku.co.jp/

株式会社トヤマhttp://www.toyama-jp.com/

株式会社

ナバテック

http://www.navatec.co.jp/

2仁木工芸株式会社

http://www.nikiglass.co.jp/

伯東株式会社

http://www.g5-hakuto.jp/

丸菱実業株式会社

http://www.ec-marubishi.co.jp/

株式会社 ラボラトリ・イクイップメント・コーポレーション

3

http://www.labo-eq.co.jp/

しょうとつ
第9巻 第5号
目 次

市川 行和

… 5

水野 智也

… 8

（総説） 低速多価イオン・原子衝突における電子捕獲過程と共鳴現象 島倉 紀之
その 1：He2+ + Li, Na, K 衝突系… 16（解説） 強い短距離斥力相関の下での弱結合少数多体系
上村正康，肥山詠美子－ エフィモフ物理と 4He 原子の 3,4 クラスター系－(後編)… 31（原子衝突のキーワード） 2 電子性再結合

中村 信行

… 41

（原子衝突のキーワード） 放射性電子再結合と放射性電子捕獲

神原 正

… 42

（原子衝突の新しい風）

後藤 基

… 43

原子衝突学会改称記念式典報告

髙橋正彦，森下 亨… 44

原子衝突学会第 37 回年会報告

行事委員会

… 44

学会財政の現状と今後の見通し

髙橋 正彦

… 45

原子衝突学会への改称に寄せて

山﨑 優一

… 47

原子衝突学会への改称に寄せて

羽馬 哲也

… 47

（シリーズ） 宇宙と原子
第三回
多価イオン
－宇宙には極端に熱いところがある－（若手奨励賞受賞研究）
高速重イオン荷電変換衝突における分子分解過程の研究

原子衝突若手の会 第 33 回秋の学校開催のお知らせ

第 33 回 秋の学校開催事務局… 48

2012 年度第 3 回運営委員会報告

庶務幹事

… 48

第 39 回定期総会報告

庶務幹事

… 48

第 13 回若手奨励賞選考理由報告

庶務幹事

… 49

国際会議発表奨励事業に関するお知らせ

庶務幹事

… 50

編集委員会からのお知らせ：「しょうとつ」の配信方法変更について

編集委員長

… 50

「しょうとつ」原稿募集

編集委員会

… 51
… 51

今月のユーザー名とパスワード

4

シリーズ

「宇宙と原子」
第三回 多価イオン－宇宙には極端に熱いところがある－
市川行和
yukitikawa@nifty.com
平成 24 年 7 月 6 日原稿受付「世の中に多価イオンというものがあり，原子物
理学的にも応用上も興味深い」ということが言われて久しくなる [1]．多価イオンの研究は現在では目新しさがなくなかったが，その分着実に続けられている．2 年に 1 度開かれている国際会議も盛んのようである [2]．
気体を熱すると衝突が盛んになり，やがて原子から電子がはぎ取られてイオンができる．つまりプラズマ状態になる．衝突がはげしくなると多数の電子がはぎ取られる．すなわち，高温のプラ

図 2: 超新星残骸 IC443 のＸ線スペクトル．人工

ズマには必ず多価イオンが含まれる．多価イオン

衛星「すざく」による観測（文献 [4] よりアメリカ天は同じ数の束縛電子をもつ中性原子と比べると，文学会の許可を得て転載）．

必ずしも同じ電子構造をもつわけではない．むし太陽コロナの分光スペクトルの中であった．

ろ「新しい原子」として振る舞う．価数の大きなイ（1869 年）ただしその実体は不明で，当初は新しオンはきわめて強いクーロン引力を周囲の物質

い元素「コロニウム」と考えられた．その後，地上に及ぼし，反応性に富んでいる．固体表面の改の実験室で多価イオンが見つかり（1920 年代）コ質など応用に役立てようという研究も進んでいロニウムは鉄の多価イオンと判明した．

る．

前回説明したように，高温のプラズマからの放地上で高温プラズマを作るには特別な装置が射は基本的に波長の短いものになる（励起状態必要だが，宇宙にはいたるところに高温プラズマ間の遷移に伴って放出される可視光が観測されが存在する．その典型的な例が太陽コロナであることもある）．図 1 は太陽観測用人工衛星「ひのる．そもそも多価イオンが初めて見つかったのはで」搭載の極端紫外線撮像分光装置（EIS）で観測した太陽コロナのスペクトルの一例である [3]．
さまざまな多価イオンの輝線が見られる．この例
では，鉄の 11 価 Fe XII の 195.12 Å の線が最も顕著である．イオン化と再結合のつり合いからイオンの存在割合をもとめる，いわゆる電離平衡の計算によると，Fe XII の割合が最大になるのは

図 1:

プラズマの温度が 126 万度のときである．このス

太陽コロナの紫外線スペクトル．横軸は波

ペクトルは太陽活動が静かなときのものである．

長（Å），縦軸は強度（相対値）．破線は感度曲線

すなわち，太陽コロナは常時その程度の高温にに相当するもの．人工衛星「ひので」による観測なっている．ところで太陽はきわめて普通の星で（文献 [3] より転載）．

Copyrigh t© 2012 The Atom ic Colli sion Soci ety of Japan, All ri gh ts reserv ed.

5

あり，そのことから高温プラズマが宇宙のどこにでもあることが推察される．
高温プラズマのもう一つの例は超新星残骸（SNR）である．IC443 と名付けられている SNRの X 線スペクトルを図 2 に示す [4]．この天体は
今から 4000 年ほどまえに爆発した超新星の残骸で地球から 5000 光年離れたところにある．このス
ペクトルは X 線観測用人工衛星「すざく」により観測された．波長分解能は 130 eV である．ここにはさまざまな原子の H 様あるいは He 様イオンのKα 線（2p-1s, 1s2p-1s2 ）が示されている．このス

図 3: 電子衝突による Fe11+ のイオン化断面積．

ペクトルの解析から約 700 万度のプラズマからの黒丸がイオン蓄積リングを用いた実験値，その上放射であることがわかる．宇宙からの X 線の観測下にある小さな黒点は誤差の範囲を表す．図のはロケットの開発とともにはじめられ，現在では大上方にある誤差棒つきの点は交差ビーム法を用型の人工衛星を用いて精力的に行われている．

いた従来の実験値 [8]．それに重なっている破線観測された X 線のうち，線スペクトルは基本的にはこれまでの推奨値．実線は理論値である．詳細

多価イオンからの放射である．

は文献[7] を参照のこと．（文献 [7] よりアメリカ

高温プラズマからの紫外線や X 線の観測には天文学会の許可を得て転載．）

多価イオンの分光学的データが必要である． またそのスペクトルの解析には多価イオンの生成

れ，2012 年現在第 7 版 [6] が使用可能である．

消滅や励起に関係する原子衝突素過程の知識

ところで，Chianti というのは，イタリアのトスカーが不可欠である [5]．特に，電子衝突によるイオ

ナ地方にある地名で，そこで生産されるワインのンの励起，イオン化，再結合，光吸収によるイオ

名前ともなっている．関係者にきいたところ，どうンの励起，イオン化，イオンと原子の電荷移行衝

してこういう名前がデータベースについたかは不突などが重要である．

明であるが，おそらく最初に名前を考える際にワ紫外線や X 線の観測の発展に伴って，それにインを飲みながら議論していて誰かが思いつき

必要な多価イオンに関する分光学的データや衝で提案したのがそのままになったのであろう，と突断面積データの整備が進められてきた．そののことである．その後よくあるように，何かの略語

成果はいくつかのデータベースとしてまとめられにならないかと努力したがうまくいかなかったそうている．その中でも特に充実しているのが，である．

CHIANTI と呼ばれるものである．これは主として多価イオンの関与する衝突過程の断面積デ太陽コロナのスペクトルの解析を目的として作らータはこのようにかなり整備されてきているが，まれたもので，水素（Z=1）から亜鉛（Z=30）までのだ不足しているものも多くあり，既存のデータもそほとんどあらゆるイオン（一部中性原子も含む）にの信頼度が十分でない場合がある．そこで現在

ついても多価イオン衝突過程の研究は精力的に行われ・エネルギー準位ている．その一例を紹介しよう．

・輝線の波長，遷移確率

多価イオンのイオン化断面積を求める実験は・電子・イオン衝突による励起・イオン化・比較的容易なので古くからおこなわれている [5]．

再結合の速度係数

その多くは標的イオンビームに電子ビームを当・その他関連する素過程データ

てる交差ビーム型実験である．しかし，通常のイが集められており，ネットを通じてだれでも利用オン源で生成されるイオンには準安定状態にあできる．1997 年に第 1 版が発表になって以来，るものが混ざっていることが多い．何らかの方法

多価イオン研究の進展に伴い随時改訂が行わでそれらの励起イオンを取り除かないと，基底状

Copyrigh t© 2012 The Atom ic Colli sion Soci ety of Japan, All ri gh ts reserv ed.

6

態にあるイオンについての断面積を得ることがで[4]

きない．もし準安定状態にあるイオンを取り除け
なければ，その分だけ得られた断面積の誤差が[5]大きくなる．
準安定状態にあるイオンを除く一つの方法は，

[6]
[7]

イオン蓄積リングを使うことである．リングの中を回っているうちに，準安定状態にあるイオンは脱励起し，基底状態にあるもののみが残る．そうな[8]ってから電子と衝突させる．その例が Hahn らの実験 [7] である．ここでは59, S857 (2007)．
H. Yamaguchi et al., Astrophys. J. 705, L6(2009)．
T. R. Kallman and P. Palmeri, Rev. Mod.
Phys. 79, 79 (2007)．
E. Landi et al., Astrophys. J. 744, 99 (2012)．
M. Hahn et al., Astrophys. J. 729. 76(2011)．
D. C. Gregory et al., Phys. Rev. A 35, 3256
(1987)．

e + Fe11+ (3s23p3 4S3/2 ) → 2e + Fe12+の断面積を求めた．その結果を図 3 に示す．
実験はハイデルベルクにあるイオン蓄積リング
TSR で行われた．イオンはパルス状に注入され，入射後数秒たってから電子と衝突させる．準安定状態（今の場合は 3s23p3 2D5/2）の寿命は 0.5秒とされており，衝突の際にはほとんどなくなってしまっている．なお，イオン化のしきい値は330.39 eV である．結果をみると，これまでに発表されていた交差ビーム法による値 [8] は今回の値より 30 %ほど大きい．これは準安定状態の寄与である．現在コロナのスペクトルの解析に用いられているのはこの古い値なので，今後変更
が 必 要 で あ る ． Hahn ら の 論 文 の タ イ ト ル（ Storage ring cross section measurements forelectron impact ionization of Fe11+ forming Fe12+
and Fe13+）をみると，純粋に原子衝突の研究であり天文学とは無関係にみえる．実際，内容も基本
的には原子物理である．しかし，それが天体物
理の雑誌（Astrophysical Journal）に掲載されていることは，宇宙と原子の研究が密接な関係にあることを示す証拠となる．

参考文献[1] 多価イオンに関する解説は多数あるが，最近のものではプラズマ・核融合学会誌 83 巻 8 号（2007 年）「小特集：多価イオン原子過程の基礎と拡がる応用研究」．
[2] 2010 年に開かれた第 15 回多価イオン国際会議（HCI）の報告は以下にある：
Phys. Scr. T144 (2011) ．
[3] P. R. Young et al., Publ. Astron. Soc. Japan
Copyrigh t© 2012 The Atom ic Colli sion Soci ety of Japan, All ri gh ts reserv ed.

7

若手奨励賞受賞研究

高速重イオン荷電変換衝突における分子分解過程の研究
水野智也
高エネルギー加速器研究機構 〒 305-0801 茨城県つくば市大穂 1-1
tmizuno@post.kek.jp
平成 24 年 7 月 18 日原稿受付
入射イオンの衝突後の価数を選別し，荷電変換衝突毎の反応過程を調べると，衝突過程をより詳細に議論できるようになる．本稿では CO と C2H2 分子の電離分解過程をそれぞれ紹介し，分解分岐比や分解イオン片の運動エネルギーの分布が衝突径数によって整理できる事を示す．また荷電変換衝突
による電離分解過程の分子配向依存性について紹介し，最後に C2H2 分子の衝突誘起による構造変形過程について議論を行う．

1．はじめにGas jet

SSDp

近年の測定技術の発展により高速重イオンとdeflector

skimmer

原子分子の相互作用の研究は盛んに行われている．特に運動量画像観測法[1,2]を用いた，衝Delay Line Detector

突時の分子配向依存，分解イオン片の運動エネ

ent
gm
fra

ルギー(KER=kinetic energy release)の分布から励起状態の同定の研究など，数多く報告されてionhole slitAq+charge selector
X1,X2,Y1,Y2
X1
Y1

いる．
高速重イオン衝突の研究は純粋に原子分子物理としての興味だけでなく，イオンビームを用

hole slit

SSDp

いた分析技術，重粒子線癌治療の原子分子過X1,X2,Y1,Y2
Y2
Digital Oscilloscope
X2
X, Y anode wires in Delay Line Detector

程の解明などの観点からも関心がもたれている．
高速重イオンが気相分子に衝突すると，電離・図 1: 実験セットアップ．

電子励起が引き起こされ，その後分解に到る．そ出したイオンビームと衝突させた．そのとき生成しれと同時に入射イオンの荷電変換も起こる．入射た分解片イオンを衝突領域にかけてある静電場

イオンの荷電変換と標的分子の電離分解過程にによって引き出し Delay Line Detector (DLD) には強い相関があり，衝突過程を十分に理解しよう

よって二次元検出位置と飛行時間を計測し，分と思うと，荷電変換毎に衝突過程を分けた研究

解イオン片の三次元運動量分布を得た．

が求められる．
本研究では衝突後の入射イオン価数を特定し，それと同時に標的分子の分解イオンを飛行時間

また同時に衝突後の入射イオンを静電偏向器
によって価数選別し PIPS (Passivated Implanted
Planar Silicon) タイプの半導体検出器によって法と画像観測法を用いて三次元運動量分布を検出した．

測定する装置の開発を行い，高速重イオンと気相分子の衝突過程を詳細に調べた．

3．CO 分子に対する高速 Si2+の衝突過程
MeV エネルギー領域の高速重イオンと二原子

2．実験セットアップ

分子の衝突過程は衝突径数，衝突時の分子配図 1 に本研究の実験セットアップを示す．実験向に強く依存する．入射イオンの荷電変換過程

槽上部から標的ガスを導入し，加速器から引きCopyrigh t© 2012 The Atom ic Col li sion Soci ety of Japan, All ri gh ts reserv ed.

8

1.2

+

CO

1.0

1e-capture

0.8
0.6

0.2

2000

2+

+

C

O

20.00

+

2+

飛行時間 2 (nsec)

CO

0.4

0.0

相対 強度 (任 意 単位 )

1200

1600

2000

2400

2800

1.2

1e-loss

1.0

C +O

1800

C +O

+
5.000

+

C +O

2+
4.000

1600
2+

C +O

1400

0.8

+

+

3.000

2+

2.000

0.6
2+

0.4

C

O

0.2

1200

2+

1.000

1200

0.0
1200

1600

2000

2400

2800

1400

1600

1800

2000

0

飛行時間 1 (nsec)

2e-loss1.2
1.0

図 3: 2 MeV Si2+衝突(1e-loss)における CO 分子の0.8

分解イオン片質量相関図．(飛行時間 1 に飛行時間の短いほうをプロット)

0.6
0.4
0.2
0.0
1200

1600

2000

2400

ンが二重電離するためには一重電離するよりも2800

衝突径数が小さくなければならない．このため飛行時間 (nsec)

2e-loss の方が分解が促進される事になる．
2+

図 2: 2 MeV Si 衝突における CO 分子の解離イオ
ン片 TOF スペクトル．

イオン片質量相関を示す．標的分子が三重電離

は衝突径数に強く依存するため荷電変換衝突が起こった際には以下の分解過程が存在する．

毎に分けた測定結果は定性的には衝突径数依(CO)3+* → C3+ + O

図 3 に 1e-loss でのイオン対解離における分解存性として理解できると期待される．
2++→ C +O

(2)

+

2+

(3)

本研究では 2 MeV Si -CO 衝突による CO 分

→ C +O

子電離分解過程を詳細に調べる事により，高速

3+

→ C+O
2+

+

(1)

2+

(4)
+

2+

重イオン衝突がどのように整理できるか検討した(C , O )イオン対解離のほうが(C , O )イオン[3]．図 2 に CO 分子の衝突による分解イオン片

対解離より収量が多い事がわかる．これは O のの飛行時間(TOF)質量スペクトルを荷電変換衝

第二イオン化ポテンシャルは 35.1 eV であるのに突毎に示す．荷電変換衝突はそれぞれ次のよう

対し C は 24.4 eV であり，(C+,O2+)イオン対解離のに 表 記 す る 事 に す る ． {1e-loss(Si2+ → Si3+) ，

ほうがよりエネルギー的に高い状態にあるためで

2+

4+

2+

+

ある．

2e-loss(Si →Si )， 1e-capture(Si →Si )}TOF スペクトルは荷電変換衝突毎に大きく異

表 1 に荷電変換衝突毎の相対電離断面積となっている事がわかる．1e-capture では e-loss と各電離度(r, COr+*)毎の分解分岐比を示した．相比べて親イオン(CO+)の収量が多く，分解イオン

対電離断面積は荷電変換衝突毎に 100 に規格の収量が少ない．また e-loss だけで比較すると化しており，四重電離以上の断面積の割合は1e-loss より 2e-loss の方が分解イオンの収量が増

1 %を切っているため無視している．

2+

えているのが見て取れる．2 MeV Si (速度 v ～先行研究として 865 MeV Xe44+(速度 v = 16

1.69 a.u.)の速度領域では e-loss に比べ e-capturea.u.) [5]，19 MeV F4+ (速度 v = 6.3 a.u.) [6]の結の方が衝突径数が大きい衝突である事が知られ果も示した．先行研究では入射イオンの衝突後ている [4]．衝突径数が大きいため標的分子にの価数選別は行われていない．この結果をみる

対するエネルギー付与が小さく，分解が抑制さと先行研究である 865 MeV Xe44+と 19 MeV F4+れると考えられる．1e-loss と 2e-loss では入射イオの結果は非常に近い値になっている．衝突を特Copyrigh t© 2012 The Atom ic Colli sion Soci ety of Japan, All ri gh ts reserv ed.

9表 1: 相対断面積と分解分岐比
2 MeVSi 荷電変換衝突(v = 1.69 a.u.) [3], 865 MeV Xe44+(v = 16 a.u.) [5], 19 MeV F4+(v = 6.3 a.u.) [6]2+

徴づけるパラメーターの一つとして Sommerfeld

1

parameter k = q/v (価数 q, 速度 v)が知られており，

いる．また k > 1 となる衝突過程においては分解
分岐比，分解イオン片の KER が入射イオン速度，（ 任 意 単位 ）した．先行研究では k > 1 の領域においては，k

相対 強度

は衝突過程を良く表すパラメーターにはならず
q/vb (b は衝突径数)をパラメーターとして用いな
いといけない事が指摘されている [7]．本研究の結果は衝突径数によって分解パターンが異なることを意味しており，k > 1 においては q/vb が衝突

+

10

20

30

40

50

60
3+*

0

10

20

30

40

50

60
3+*

10

20

30

40

50

60
4+*

(KER)の分布を荷電変換毎に示した．KER においても荷電変換毎に大きく変化している事がわ

0

かる．1e-capture において KER 分布は最も小さく，
1e-loss，2e-loss になるにつれて KER 分布が大き

80
2+

70
2+

CO ->C +O

1.2
1.0
0.8
0.6
0.4
0.2
0.0

次に図 4 に分解イオン片の運動エネルギー

80
+

70
+

CO ->C +O

1.2
1.0
0.8
0.6
0.4
0.2
0.0
0

る．

70
2+

CO ->C +O

1.2
1.0
0.8
0.6
0.4
0.2
0.0

を特徴づけるパラメーターである事を支持してい+CO ->C +O
1e-capture
1e-loss
2e-loss0

価数，イオン種に殆ど依存しない事が報告され領域においても分解過程が大きく異なる事を示

ク ーロ ン 爆発モ デル2+*

k > 1 では非摂動論的な取り扱いが必要とされて換毎に分けて分解分岐比を測定すると，k > 1 の

1 

 2

1.2
1.0
0.8
0.6
0.4
0.2
0.0

k < 1 では摂動論を用いて衝突過程を表現でき，ている [7, 8]．本研究(k = 1.18)において荷電変

3 

20

40

60

80

100

80
2+

120

KER (eV)

図 4: 2 MeV Si2+衝突における CO 分子 KER スペク
トル(縦線はポテンシャル曲線から予想される値)．

くなっている．KER の大きさは励起エネルギーの大きさと関係しており，衝突径数の小さい，エネ

Copyrigh t© 2012 The Atom ic Colli sion Soci ety of Japan, All ri gh ts reserv ed.

10

ルギー付与の大きな衝突の場合は同じ電離度においても標的分子をエネルギー的に高い状態に励起していると言える．

C

KER の値は単純に純粋なクーロンポテンシャ
ルによって記述すると以下の式で近似できる(クーロン爆発モデル)．
KER = q1q2/R

(5)

入射イオン

O

q1 は炭素イオンの電荷，q2 は酸素イオンの電荷，R は炭素原子と酸素原子の核間距離である．

図 5: 分子配向 の定義．(入射イオンの進行方
向と炭素イオンの運動量ベクトルのなす角で定義
している)

KER の分布は基底状態の波動関数を用いて核間距離分布を計算し，(5)式に代入する事により求めている．KER 分布は純粋なクーロン爆発モ
デルでは再現できない構造を有していることが知られており，イオン化状態ごとにポテンシャル曲

1e-loss線を計算し，フランク・コンドン原理により基底状

態から電離状態に垂直遷移する事を用いて解離状態との差として KER を計算する事で再現する

dN/d (任 意 単位 )事が出来る [9]．図 4 に示した縦線は存在するイ
オン化状態の KER の値である．(CO)2+* → C+
+ O+イオン対解離に関しては対応する電子状態も示しておいた．
本研究において k > 1 のエネルギー領域においても荷電変換衝突を分けて測定すると，分解

0

20

40

60

80 100 120 140 160 180

1e-capture

分岐比，KER が大きく変化していることがわかった．衝突過程をより良く理解するためには衝突径0数を分けた測定が必要であり，荷電変換衝突を

20

40

60

80 100 120 140 160 180

 (degree)

分けて測定すると定性的ではあるが衝突径数の議論が出来る事が示された．

図 6: (C2+，O+)分解過程の配向依存性

4．入射イオン荷電変換衝突の CO 分子
配向依存性

O4+-CO 衝突における 1e-loss と 1e-capture 毎の(C2+, O+)分解過程における配向依存性を示した．

標的分子の電離過程は衝突時の分子配向に高速領域でのイオン衝突での電離微分断面積

強く依存し，二原子分子の場合では，多重電離の配向依存性は以下の式で近似できる事が知ら

過程は分子軸と入射イオンの進行方向が平行にれている [11]．

なるような衝突において優勢に起こる事が報告さ

d/dP2(cos))

れている [7,8]．標的分子の電離過程の配向依(6)

ここで N0 は規格化定数，は異方性パラメータ

存性は，標的分子の電子密度の異方性によってー，P2(x)=(3x2-1)/2 でありルジャンドルの二次の入射イオンの進行方向と分子軸が平行な衝突の多項式である．立体角 dは sinddであるため

方がエネルギー付与が大きくなる事から説明さ相対強度 N の微分は以下の式になる．

れている [10,11]．

dN ndP2(cos))sindd(7)

本研究では，入射イオンの荷電変換過程が分I は入射イオン強度であり，n は標的数である．N

子配向から受ける影響について調べた [12]．図の配向角  での微分は以下の式になる．(方向

5 に分子配向 の定義を示す．図 6 に 6 MeVで積分することにより 2の係数が掛かっている)

Copyrigh t© 2012 The Atom ic Colli sion Soci ety of Japan, All ri gh ts reserv ed.

11

dN/d=2P2(cos))sin(8)
電離微分断面積に配向依存性が無い時(=0)，2500

C2H+ ( m = 0 - 2 )1e-capture

相対強度の角度微分 dN/dは sinの依存性を

m

2000

示す事になる．

H

カウ ン ト

配向依存性は荷電変換衝突によって大きく異なり 1e-loss と 1e-capture では正反対の依存性を示すことが明らかになった．1e-capture において+

1500
CH+ ( m = 0 - 2 )
m

1000
H+

は sin分布に近く分子軸と入射イオンの進行方

500

2

C
x10

2000

2400

向が垂直のとき起こりやすい傾向を示している．
01600

相対強度分布が sinに比較的近いことから電離
微分断面積の角度依存性は小さいと言える．

2+

x10

x10

2800

3200

飛 行時 間 ( n s e c )他方で 1e-loss においては分子軸と入射イオン

図 7: 2 MeV Si2+ 1e-capture(Si2+ →Si3+)におけるC2H2 分子の解離イオン片 TOF スペクトル

の進行方向が平行な場合において優勢に起こっており sin分布から大きくずれている．1e-loss では O4+から O5+に電離するために 114 eV もの大きなエネルギー付与が中性分子から与えられる必

二つのピークを持つ構造をしている．放射光を要がある．そのため必然的にエネルギー付与の用いた先行研究で，C2H2 分子の C 1s に空孔を大きな分子軸に平行な衝突が優勢に起こると考あけ Auger 過程によって二価の C2H2 を生成させ，えられる．逆に 1e-capture では O4+から O3+への(C2H+，H+)に分解した過程の KER が報告されて

電子移行反応は発熱過程であるため衝突径数

いる [15]．その結果では KER は 3.5 eV と 5.3 eVが大きな場合でも起こる．衝突径数が大きい場の二つのピークを持つ構造をしている事が報告

合,標的分子の電子密度の異方性は小さく分子されている．また分解イオン片の KER と同時に軸に大きく依存しなくなる．そのため 1e-capture

Auger 電子の運動エネルギーと角度分布も測定

衝突においては配向依存性が小さくなると言えされており，3.5 eV の KER の場合は二価の C2H2分子の電子状態は 1g+と同定され，KER が 5.3

る．

eV の場合は 1u と同定されている．よって本研究で観測された 3 eV と 5 eV の KER はそれぞれ

5．C2H2 分子の電離・分解・構造変形過程

g+と 1u と考えられる．2e-loss において 1u が1

顕著に観測される理由は，より高い励起状態で標的分子が多原子分子となった場合，衝突誘ある 1u は一番近接衝突でエネルギー付与の大起による緩和過程には二原子分子標的では存

きい 2e-loss で最も生成されやすいためである．

在しない結合の組み換えや構造変形などの過程

図 9 に(CH+,CH+)分解過程の KER 分布を示しが存在する．分解イオン片の運動エネルギー分

た．2e-loss は収量が十分でないために載せてい布(KER)や運動量ベクトル相関から分解経路のない．3～4eV 付近にピークを持つ構造をしてい同定と構造変形過程について調べた [10,11]．

る事がわかる．先ほどの放射光の結果では，6.5

図 7 に 2 MeV Si2+-C2H2 における 1e-capture

eV にピークを持ち，この状態として 1g+と同定さ衝突における分解イオン片の TOF 質量スペクトれている [15]．明らかに本研究での KER は放ルを示した．H2+や CH2+などが観測され，組み換射光の結果と比較して小さい値を取っている．重え反応がイオン衝突によって誘起されている事

水素化アセチレン C2D2 を用いた，800nm の波長がわかる．

でパルス幅 9 fs，強度 0.13 PWcm-2 の超短パル

図 8 に(C2H+, H+)分解過程における分解イオ

スレーザーでの実験[16]でも放射光の結果同様

ン片 KER 分布を示した．1e-capture と 1e-loss にに 6.5 eV 近辺にピークを有している．他方で 1.2

おいては 3 eV にピークを持つ分布となり，非常MeV Ar8+衝突の実験結果も報告されており，2.5に良く一致した．2e-loss において 3 eV と 5 eV に

eV 近 辺 に ピ ー ク を 有 し て い る [17] ． 更 にCopyrigh t© 2012 The Atom ic Colli sion Soci ety of Japan, All ri gh ts reserv ed.

12

1.2

1.2

相対 強度 (任 意 単位 )0.8

相対 強度 (任 意 単位 )1e-capture
1e-loss
2e-loss

1.0

0.6
0.4

1.0

1e-capture
1e-loss

0.8
0.6
0.4

0.2

0.2

0.0

0.0
0

2

4

6

8

10

12

14

16

18

20

0

2

4

6

8

12

14

16

18

20

図 8: 2 MeV Si2+ 衝突における(C2H+,H+)分解過程
における KER 分布

図 9: 2 MeV Si2+ 衝突における(CH+,CH+)分解過程における KER 分布

1e-capture 同様に KER 分布がダブルピークにな

態からの分解を考慮する必要があり，禁制遷移ることも報告されている．つまりイオン衝突では光状態のポテンシャル曲面の計算がまたれる．

照射とは異なる分解過程が存在する事を示して

次に C2H2 分子のイオン衝突による構造変形

いる．

過程の研究を行った．6 MeV O4+-C2H2 衝突におイオン衝突では光では誘起されない禁制遷移

ける分解過程の C+の運動量ベクトル P(C+) と H+も起こるため，光照射では生成されない電子状

2+

C C

+

+

2+

C +C +H

H

0.20

C

+

H

0.18
0.16
0.14
0.12
0.10
0.08

の運動量ベクトル P(H+)のなす角 をそれぞれ

+

2+

C C

+

+

C

2+

C +C +H

H

C

+

H

0.16

+

+

+

C

0.14

2+

0.12
0.10
0.08
0.06
0.04
0.02

20

40

60

80

100

120

140
+

0.14

160
+

C +C +H

0.08
0.06
0.04
0.02
0.00

40

60

80

100

120

0.18

140

160

+

+

C +CH +H

180
+

0.14
0.12
0.10
0.08
0.06
0.04
0.02
0.00

20

40
+

60

80

100

120

140

160

40

60

80

100

120

180

140
+

160
+

C +C +H

0.16

180
+

0.14
0.12
0.10
0.08
0.06
0.04
0.02
0.00

0

20

40

60

80

100

120

0.20
0.18
0.16
0.14
0.12
0.10

0.16

0

20

0.18

0.10

20

0

+

0.12

0

0.00

180

相対 強度 (任 意 単位 )0

相対 強度 (任 意 単位 )

+

+

0.18

2+

0.06
0.04
0.02
0.00

140

160

+

+

C +CH +H

+

180

140

160

180

0.08
0.06
0.04
0.02
0.00

0

20

40

60

80

100

120

+

C と H の な す角

10KER (eV)

KER (eV)C +と H +の な す 角

図 10: 6 MeV O4+ 衝突における C+の運動量ベクトルと H+の運動量ベクトルのなす角．
左:1e-loss (O4+→O5+)，右:1e-capture(O4+→O3+)
Copyrigh t© 2012 The Atom ic Colli sion Soci ety of Japan, All ri gh ts reserv ed.

13

1e-loss, 1e-capture ごとに(H+，C+，C2+)，(H+，C+，
+

+

+

れる．

+

C )，(H ，C ，CH )分解過程にわけて図 10 に示した．1e-loss と 1e-capture で大きな違いは観測さ6．まとめれていない．次に 4 体解離(H+， C+， C2+)，(H+，高速重イオン荷電変換衝突に伴う気相分子のC+ ， C+)においてピークが二つ観測されている

電離分解過程を分解イオン片運動量画像観測のは H+がどちらの C と結合していたかを見ている

法を用い，分解イオンの運動エネルギー(KER)，

ことに対応している．すなわち H と結合している

入射イオンに対する分子軸依存性，分解イオン

ほうの C 原子とは同じ方向の運動量ベクトルを持

片運動量ベクトル相関から議論してきた．

ち，結合していない方の C 原子とは逆方向の運

Sommerfeld parameter k が 1 以上の値を持つ衝

動量ベクトルを持つ．4 体解離の場合 0 度と 180

突過程においても荷電変換毎に分けて測定す

度が直線分子に対応し，20 度と 160 度のピーク

れば分解分岐比，KER は大きく異なることを明ら

は C-C 結合からどちらとも 20 度折れ曲がって H+

かにし，衝突過程を特徴づけるパラメーターとしが放出されていることを意味している．他方で 3て衝突径数 b を含めた q/vb を用いる必要がある

体解離の場合は H+は必ず C+と結合していたこと

ことを改めて指摘した．

が自明であることから(もう一方が CH+であるため)，

分子軸依存性では荷電変換の過程に必要なエH+と C+の運動量のなす角はそのまま C-C 結合か

ネルギー付与量によって分子配向依存性が議らの折れ曲がりを表す．4 体解離では 20 度の折論できることを明らかにした．

れ曲がりでピークを持つのに対し，3 体解離では

C2H2 を標的に用いた実験では，分解過程によ40～50 度にピークを持つ．つまり 3 体解離の方っては光照射とは異なる状態を経由した解離過が 4 体解離より分解前に構造が直線構造から屈程がイオン衝突では存在することを明らかにし，

曲構造へ変形していることが言える．

禁制遷移を経由した解離過程をイオン衝突では

以下では(H+ ，C+ ，CH+)分解過程において構

考慮しなければならない事を改めて示した．イオ

造変形が起こるメカニズムを考察する．放射光をン衝突による構造変形を議論し，そのメカニズム



用 い た 先 行 研 究 で C 1s →  遷 移 に よ るとして価電子励起，振動励起が考えられることをRenner-Teller 効果によって直線分子から屈曲し示した．

た分子への構造変形が報告されている [18]．本
研究の速度領域では C 1s 内殻電離や励起過程

7．謝辞

は全断面面積において 2～3 %と非常に小さく，

本研究は京都大学大学院工学研究科 原子

入射イオン速度が 1 a.u. 近辺の衝突過程では核工学専攻に在学中に伊藤秋男教授の御指導無視することができる．二価の C2H2 分子においの下行ったものです．ここに改めて感謝いたしまては価電子励起状態において屈曲構造に構造

す．また土田秀次准教授(京都大学)，間嶋拓也

変形することが予想されており [19]，実際に超助教(京都大学)には本研究の遂行に対して数多

短パルスレーザーを用いたポンプ-プローブによくのご助言を頂いたことを感謝いたします．共同

り二価イオンのビニリデン構造(CCH2)への構造研究者として研究に参加して頂いた中井陽一博

変形が実時間で観測されている [16]．三価イオ

士(理研)には研究全般にわたり数多くのご助力

ンがどのようなポテンシャル曲面を持っているの頂いた事感謝いたします．

か報告されていないため確定的な事は言えない

また共に実験を行った京都大学原子核工学

が三価イオンにおいても構造変形を起こす価電専攻量子システム工学研究室の皆様方に深く御子励起状態が存在する可能性は高い．他方で礼を申し上げます．

120 keV Ar8+-CS2 衝突では振動励起によって屈曲が誘起されていると結論づけられている [20]．
これらのことを総合すると構造変形は電子励起も参考文献[1] J. Ullrich et al., Rep. Prog. Phys. 66, 1463

しくは振動励起によって誘起されていると考えら(2003).

Copyrigh t© 2012 The Atom ic Colli sion Soci ety of Japan, All ri gh ts reserv ed.

14

[2] J. Ullrich et. al., J. Phys. B 30, 2917 (1997).
[3] T. Mizuno, T. Yamada, H. Tsuchida, Y. Nakai,
and A. Itoh, Phys. Rev. A 81, 012704 (2010).
[4] S. De, P. N. Ghosh, A. Roy, C. P. Safvan,
Nucl. Instr. and Meth. B 243, 435 (2006).
[5] L. Adoui et al., J. Phys. B 32, 631 (1999).
[6] Ben-Itzhak et al., Phys. Rev. A 47, 2827
(1993).
[7] V, Krishnamurthi, I. Ben-Itzhak and K. D.
Carnes, J. Phys. B 29, 287 (1996).
[8] B. Siegmann et. al., Phys. Rev. A 62, 022718(2000).
[9] M. Lundqivst et. al, Phys. Rev. Lett. 75,
1058 (1995).
[10] U. Werner et al. Phys. Rev. Lett. 79, 1662
(1997).
[11] B. Siegmann et al., Phys. Rev. A 66, 052701(2002).
[12] T. Mizuno, T. Majima, H. Tsuchida, Y. Nakai,
and A. Itoh, J. Phys. Conf. Sires 58, 173(2007).
[13] T. Mizuno, T. Yamada, H. Tsuchida, Y. Nakai,
and A. Itoh, J. Phys. Conf. Sires 163, 12039(2009).
[14] A. Itoh, T. Mizuno, T. Majima, J. Phys. Conf.
Ser. 288, 012027 (2011).
[15] T. Osipov et al., J. Phys. B 41, 091001(2008).
[16] A. Matsuda, M. Fushitani, E. Takahashi, and
A. Hishikawa, Phys. Chem,. Chem. Phys. 13,
8697 (2011).
[17] S. De, J. Rajput, A. Roy, P. N. Ghosh, and C.
P. Safvan, Phys. Rev. A 77, 022708 (2008).
[18] N. Sato et al., Chem. Phys. Lett. 393, 295(2004).
[19] T. S. Zyubina et al., J. Chem. Phys. 123,
134320 (2005).
[20] F. A. Rajgara et al., Phys. Rev. A 64, 032712(1999).

Copyrigh t© 2012 The Atom ic Colli sion Soci ety of Japan, All ri gh ts reserv ed.

15

総

説

低速多価イオン・原子衝突における電子捕獲過程と共鳴現象
その 1：He2+ + Li, Na, K 衝突系
島倉紀之
新潟大学理学部
shima@shimakura.gs.niigata-u.ac.jp
平成 24 年 5 月 30 日原稿受付

重粒子衝突非弾性散乱過程において形状共鳴はいまだ観測されていない．我々は，これまで，共鳴観測に供するため，量子論的緊密結合法，波束法を用いて，共鳴現象を系統的に調べる理論的研究を行ってきた．今回は，その中から特に He2+イオンとアルカリ金属原子(Li, Na, K)という衝突系を取り上げ，これまでの理論的研究成果をもとに，形状共鳴について解説する．次回は，これ以外の衝突対を取り上げ，入射イオンの電荷依存性，原子番号依存性などについて解説する予定である．非弾性過程を利用して形状共鳴状態を作るという要求は最近特に強くなっている．
なり以前からさまざまな衝突系において観測され1．はじめに共鳴は，物理，化学，生物の分野で広く観測ており，我々も Arq+ + H2 衝突系の電子捕獲断されている普遍的な現象であり，何故そのような面積に現れたオービティング共鳴を報告した [1]．

現象が起こるか，その現象はどのような観測に現

また，オービティング共鳴の断面積への関与をれるか研究されてきた．原子衝突の分野では，明確にするため，N4+ + H 衝突系に対して，実共鳴は主として電子衝突の分野で良く知られた

験，理論の共同研究を行った [2]．もう 1 つは現象である．例えば，電子と分子の衝突における0.1 eV u-1 以下の極低衝突エネルギー領域にお共鳴の重要性は，弾性散乱，振動励起，電子付

いてある特定の衝突エネルギーで断面積がピー

着のような過程において長いこと認識されてきた．

ク状になる形状共鳴である．古い論文で形状共ある振動励起過程は，共鳴のため，狭いエネル

鳴をオービティング共鳴と呼んでいるものもある

ギー範囲ではあるが 2，3 桁もその断面積の値がので注意が必要である．この総説で扱うのは 2 つ大きくなる．電子付着により生成された負の電荷

目の形状共鳴である．

を持つ分子イオンは，不安定であることが多く，

形状共鳴は古くは 1898 年にラジカル—ラジカ

短い寿命で解離反応を起こし，ラジカルやイオン

ル再結合反応を説明するためボルツマンによっのようないろいろな種類の解離フラグメントを生じて提案されたが [3]，孤立衝突系における形状

る．それ故に，共鳴の知識は電子付着を経由し共鳴の実験による観測は，それから 75 年後のカた解離過程において，その生成物をコントロール

ナダの Waterloo 大学の Scols のグループによする上で重要である．

る Hg 原子からの H 原子の弾性散乱断面積の

電子衝突の場合に比べ重粒子衝突の場合，測定まで待たなければならなかった [4]．ほぼ時共鳴現象の研究は十分に行われているとは言えを同じくして， Toennies のグループは水素原子，ない．重粒子衝突における共鳴現象は 2 種類に水素分子と希ガス（Ne, Ar, Kr, Xe）の弾性衝突分類される．1 つは 10 eV u

-1

以下の低衝突エ過程で現れる重粒子共鳴現象の観測結果を報

ネルギー領域において衝突エネルギーが減少

告した [5]．これ以降も，重粒子 2 体衝突においするにつれ断面積が増加するいわゆるオービテていくつかの系の弾性衝突過程で形状共鳴がィング共鳴である．このオービティング共鳴はか観測されている．一方，理論的な研究としては

Copyrigh t© 2012 The Atom ic Colli sion Soci ety of Japan, All ri gh ts reserv ed.

16

1978 年のカナダの Alberta 大学の Davis とThorson の報告が初めてである [6]．彼らは H+

He2+ + アルカリ金属原子(ns)→ He+(3l) + 1 価のアルカリ金属イオン (1)+ H(1s), D(1s) 衝突系における弾性散乱，共鳴型電荷移行過程を理論的に扱った．状態間の(アルカリ金属原子 = Li, Na, K)遷移を含む非弾性散乱過程における共鳴は，
1984 年に，(N3+ + H) 衝突系における電子捕獲を取り上げ解説する．扱ったエネルギー領域で過程において，Rittby らによって初めて報告さは電子捕獲に寄与するのは He+(n = 3) チャン

れた [7]．非弾性散乱過程に対して形状共鳴をネルであることが知られているので，電子捕獲チ扱った 2 番目の報告は我々の (N5++ H) 衝突ャンネルとしてそれらのみ考慮した．Scols らのグ系における電子捕獲過程の研究である [8]．こループ，Toennis らのグループ，Davis と Thorson

れ以降，我々の研究も含め，形状共鳴に対してによって扱われた衝突系はいずれも共鳴は弱い

多くの計算結果が報告されるようになった．

化学結合力によって生じる．これに対して，多価

我々は，N

5+

+ H 衝突系における形状共鳴の報告以後 [8]，C
2+

系，Be

5+

+ H 衝突系，B+ Li 衝突系，N4+3+

イオンによる中性粒子からの電子捕獲過程では，+ Li 衝突3++ H 衝突系，Be

共鳴現象は，入射イオンの電荷により中性標的+

原子が分極し，ポテンシャルエネルギー曲線に+ H 衝突系 [10]に対して，引力部分（ポテンシャルの井戸）が現れ，遠心力

量子論的緊密結合法にもとづき断面積の計算を（斥力）を考慮した有効ポテンシャルに障壁が生行うとともに波束法を適用することにより，共鳴現

じることによる．注意したいのは，弱い化学結合

象解明の研究を続けてきた．その結果，(1) 共

力に比べて多価イオンによる中性粒子の分極力

鳴は衝突の途中で形成される準分子の振動回はさらに弱いということである．多価イオンによる

転エネルギーと衝突エネルギーが一致すると起電子捕獲過程の場合，標的原子の分極率が大こること，(2) ある特定の部分波の散乱行列が大

きいと，大きい軌道角運動量で，従って大きい衝

きくなることにより断面積がピーク状に大きな値を突エネルギーで共鳴が観測できる可能性がある．

持つようになることを見出した．また，(3) 共鳴が

我々が He2+とアルカリ金属原子という衝突系をと高い衝突エネルギーで観測されるには入射イオりあげたのは，これらがンの荷数，標的の分極率が大きい必要があるが，(1) 実験で取り扱い易い衝突系である，

その反面，そのような場合，多くの状態が複雑に(2) 大きな双極子分極率を持つ標的と多価入絡み合い，共鳴が観測しづらくなることもあること

射イオンの組み合わせである，を明らかにした．これらの結果は日本物理学会(3) 電子捕獲断面積の値が大きい，

2002 年秋季大会のシンポジューム（主題：低速(4) 共鳴ピークとそのバックグラウンドの差が著し

多価イオンの理論）で筆者が報告済みである．こい，の頃から，世界中で，電子捕獲過程に現れる共と考えたからである．得られた結果は，実験によ

鳴ピークを実験によって観測しようという試みがる共鳴現象の観測が予想以上に難しいことを示始まったが，現在のところすべて徒労に終わってしたが，それでもこれまで扱われた系に比べれいる．これは，共鳴現象が極端に低い衝突エネ

ば，共鳴が観測される可能性を示していると思わルギーで現れること，更に入射イオンの衝突エネ

れる．

He 衝突系 [9] ，N

5+

ルギーを揃えて衝突させなければならないことに非弾性過程を利用して形状共鳴状態を作るとよる．

いう要求は最近特に強くなっている．これは，極そこで我々は，共鳴の観測を行ってもらうため，

端に冷えた分子に対して様々な応用が考えられ

電子捕獲過程に現れる共鳴の特徴を系統的にるからであり，特定の量子状態だけを作る，また調べる理論的研究を行った．今回はその中から，余分なエネルギーを持たせないという観点から

2+

特に，He イオンとアルカリ金属原子という衝突

いうと弾性散乱過程よりも非弾性過程を用いた方対の組み合わせが有利だからである．また，最近の冷たい原子，

Copyrigh t© 2012 The Atom ic Colli sion Soci ety of Japan, All ri gh ts reserv ed.

17

分子を作る技術の進展がこのような研究に拍車で，非常に重要である．しかし，量子論的緊密結をかけている [11]．

合法では，断面積の数値を得ることはできるが，本論文では，断りがない限り原子単位系を用

共鳴のダイナミクスのイメージは把握しにくい．こいる．

のような理由から，我々は，原子核の運動も量子的に扱いながら，状態間の遷移も含めて衝突の2．理論と計算法

追跡ができ，共鳴のイメージを描くことのできる多2.1 電子状態の計算

チャンネル核波束法を開発した [10]．1 つのポ共鳴現象が現れる極低衝突エネルギー領域

テンシャル曲線上に核波束を走らせる核波束法

では衝突系全体を分子として扱う分子基底が良はポピュラーなものであったが，状態間遷移が起い近似となる．また，状態間の遷移を引き起こす

こり複数の状態を考慮しなければならない衝突

結合項としては，動径結合のみ考慮すれば良く，系を扱える核波束法の開発はコンピュータの進

回転結合は遷移にほとんど影響を与えない．更歩（計算速度と記憶容量）した最近になってように，遷移は断熱ポテンシャル曲線の擬交差点付

やく可能となった．我々の開発した多チャンネル

近でのみ起こると考えて良い．

核波束法では，系のハミルトニアンと初期波束を

この総説で扱う衝突系では，内殻電子は衝突与えることで，任意の時間における波束を得るこ

動力学に影響を与えないと考えることができるのとができる．初期波束としてはガウス型波束を用で，原子核と内殻電子をまとめてコアとしてガウス

いるが，それに含まれるパラメーターの具体的な型の擬ポテンシャルで表し，基底関数としてはス

決定法については論文[10]を参考にして欲しい．

レーター型軌道を用い，配置間相互作用法を使

アルゴリズムとして，より正確な結果を得るためって，準分子の固有状態（波動関数）と固有値をsplit-operator 法を，断熱基底—透熱基底間の変求めた [12]．

換を効率的かつ正確に行うため高速フーリェ変換法を用いている．核波束法の原理的な事柄に

2.2 量子論的緊密結合方程式

関しては Balakrishnan らの論文[15]を，また，いく

断熱基底の量子論的緊密結合方程式は 2 次

つかの点で我々の方法と異なるが，Vaeck ら[16]微分の結合項を含むこと，結合項が核間距離のが多チャンネル核波束法の開発を行っているの狭い範囲で大きな値となることから，数値的に解で，それらの論文も参照して欲しい．

くのは不便である．一方，透熱基底では，核間距

時間発展の結果得られる波束を用いて，スペ離の変化に対して状態が比較的ゆっくりと変化

クトル法に基づき固有値，固有関数を得る方法もするので，数値積分の際に刻み幅をそれほど細開発した [10]．この方法では，数値積分を効率かく取らなくても断面積が精度良く求まる [8, 13,

的に行うためウインドウ関数を用いるなど工夫を

14]．緊密結合方程式は部分波（全角運動量 K）

しているが，基本的には Feit らの方法[17]を用いごとに数値的に解かれ，i 状態から f 状態へのている．量子論的緊密結合法によって得られる S遷移に対応する散乱行列 S Kfi が得られる．全散行列を調べることにより特定の衝突エネルギーで乱断面積は得られた S Kfi 行列を用いて準分子の回転量子数を知ることができる．一方で，

 fi E 


ki

2

 2K 1S

K
fi

2

現れた共鳴ピークに寄与する部分波を，つまりはスペクトル法によって，有効ポテンシャルの井戸（2）

中の準分子の振動回転状態のエネルギーが求K

まる．この両者を比較することによって，2 つの方と書ける．ここで ki は初期波数である．

法の計算精度の確認がある程度できる．また，2つの方法を比較するもう 1 つの利点は共鳴がど

2.3 多チャンネル核波束法とスペクトル法

量子的な現象である共鳴の物理的イメージをのポテンシャル井戸によるものなのか知ることができる点にある．しかし，量子論的緊密結合法に描くことは，そのダイナミクスを理解する上で，更よるピークの位置は電子捕獲に関与する複数のにはどんな系で共鳴が観測し易いか予測する上Copyrigh t© 2012 The Atom ic Colli sion Soci ety of Japan, All ri gh ts reserv ed.

18

状態を同時に考慮して得られた結果であり，スペクトル法に基づく同定はあくまでも 1 つの擬状態（断熱と透熱を組み合わせた）のポテンシャルに基づく結果であるという違いに注意して欲しい．

3．形状共鳴の性質
この章では，まず 3.1 節で量子論的緊密結合法に基づく計算結果において共鳴がどのような形で現れるか解説し，3.2 節で共鳴が何故生じるのかを有効ポテンシャルの立場から解説した後，

図 1: 散乱行列の衝突エネルギーおよび全角運
動量量子数依存性．(B3+ + Li)衝突系の衝突エネ
ルギーが 0.04548 eV で現れる振動，回転量子数
が(0,189)に対応した共鳴の場合を例として示した．

3.3 節では「量子論的緊密結合法とスペクトル法
を組み合わせて共鳴ピークを同定する」方法について述べる．これにより我々の得た計算結果の理解の一助にされたい．
3.1 量子論的緊密結合法の計算結果から見た形状共鳴式(2)で角運動量 K に関する和は 0 から散乱行列が断面積に寄与しなくなる十分大きな全角
運動量 K max までとる．形状共鳴が起こらない場
合， K max の値としては運動量 p と電子捕獲が起こる衝突径数 b の積が 1 つの目安となる．つまり，

 衝突エネルギーが大きいと， K max は大きくなり，
 大きな K max まで和を取る必要が生じる．しかし形

図 2: B3+ + Li → B2+(4d,4f) + Li+電子捕獲断面積の衝突エネルギー依存性．カッコ内の 2 つの数
値はそれぞれ準分子の振動量子数と回転量子数を表している．

状共鳴が起こる場合注意が必要である．一例として，図 1 に，B3+ +
Li 衝突系の衝突エネルギー が 0.04548 eV の共鳴ピーク付近における，散乱

これが，図 2 に示した B3+ + Li → B2+(4d, 4f) +行列の衝突エネルギーおよび全角運動量量子Li+ 電子捕獲過程において，衝突エネルギー数依存性を示した．K = 0−185 では散乱行列の0.04548 eV で，断面積がピークになる原因であ

角運動量依存性は衝突エネルギーが変わってもる．量子論的緊密結合方程式を解いて断面積をあまり変化しないなめらかな振動関数になってい求めるとき，我々は，K = pb の関係式から K max をる．またその振動は，K が小さい領域ではゆった

求めそこまでの和を取るとか，更には K に関するりしているが，K が大きくなるにつれ短い周期の振幅の大きい振動に変わって行く．少々見にく

和を 1 つおきにとるとかしがちであるが，そのようなことをしてしまうと共鳴現象が起こっていること

いが， K  186で振動はほとんどなくなり，振動のを見逃してしまうことがある．

外側の K > 186 では，K = 189 においてのみ散乱

行列が大きな値を持っている．この K = 189 のと3.2 有効ポテンシャルの立場からみた形状共きの散乱行列の異常さはその大きさにある．K =鳴0−185 では散乱行列の二乗の値は最大で 0.02

断熱ポテンシャルに遠心力ポテンシャルを加程度であるが，K = 189 では約 0.85 と 50 倍にもえた有効ポテンシャル

なる（図 1 では縦軸は 0.2 で切ってある）．このようVeff ( R)  V ( R) に，共鳴現象が起こる場合，散乱行列はある特定のそれも 1 つの角運動量で大きな値を持つ．

K (K 1)
2R2(3)

Copyrigh t© 2012 The Atom ic Colli sion Soci ety of Japan, All ri gh ts reserv ed.



19

4 標的原子依存性
共鳴は非常に低い衝突エネルギーで起こる現象なので，できるだけ高い衝突エネルギーで起
こる衝突系を探す必要がある．一方で，電子捕獲過程は反応過程であり，共鳴が観測されるにはその反応断面積が大きいということも必要である．Li, Na, K はいずれも双極子分極率が大きい（Li, Na, K の双極子分極率はそれぞれ 162 au3,
162 au3, 290 au3）ので，He2+ + Li, Na, K 系のポ
テンシャル井戸は深くなる．ポテンシャル井戸が
図 3: 断熱ポテンシャル V に遠心力ポテンシャル

深くなると，図 3 からわかるように，角運動量が高K (K 1)を加えた有効ポテンシャル Veff.
2R2

くなっても，有効ポテンシャルの井戸が埋まらず，
大きな振動回転エネルギーをもつ不連続状態が存在することになり，その結果，高い衝突エネルを考える．モデル図を図 3 に示した．ポテンシャ



ギーで共鳴を観測できる可能性が生じる．最初

ルエネルギー曲線が一定の条件を満たす引力に，反応が起こり易いかどうか予測する目的で，型で，角運動量が 0 でない場合，有効ポテンシャ

断熱ポテンシャルエネルギーについて，特に擬ルはポテンシャル障壁を持ち，それよりも内側に

交差点の位置およびその点でのエネルギー差とポテンシャル井戸を持つ．この井戸がある程度

電子捕獲断面積の大きさとの関係に注目して，深いと井戸の中に不連続状態である振動状態が解説する．引き続き，共鳴の同定とメカニズムにできる．衝突対は，トンネル効果によりポテンシャ

ついて解説する．

ル障壁にしみ込み，この振動回転状態のエネル
ギーと衝突エネルギーが一致すると，一時的に

4.1 断熱ポテンシャルエネルギー曲線と擬交

井戸中に振動回転状態として留まることになる．

差点

角運動量がある値以上に大きくなると，井戸は消HeLi2+系，HeNa2+系, HeK2+系の断熱ポテンシ

滅する．

ャルエネルギー曲線を，それぞれ，図 4，図 5，図6 に示す．また，擬似交差点でのエネルギー差を

3.3 共鳴ピークの同定法

表 1 に，ポテンシャル井戸の深さを表 2 にまとめ

電子捕獲断面積において現れた共鳴ピーク

た．断熱状態 4Σ から 8Σ はそれぞれ核間距離無の衝突エネルギー付近で，散乱行列の部分波

限遠で(He+(3s) + Li+, Na+, K+)，(He+(3p) + Li+,

依存性を調べることにより，共鳴に寄与するただNa+, K+)，(He+(3d) + Li+, Na+, K+)，(He2+ +1 つの部分波の角運動量がわかる．更に，この角Li(2s), Na(3s), K(4s))，(He2+ + Li(2p), Na(3p),運動量をもとに遠心力ポテンシャルを求め，その有効ポテンシャル上で波束を走らせることにより，振動状態のエネルギーと固有関数を求めること

K(4p)) 状態に繋がっている．核間距離が大きい
領域では，4Σ から 6Σ のポテンシャルエネルギー曲線は核間距離が大きくなるにつれ右下がりにができる．スペクトルパワーのエネルギーを衝突なっているが，これはそれらの状態では衝突対対の距離が無限大のときのエネルギーを基準にの間にクーロン反発が働いているためである．こすることにより，電子捕獲断面積において現れたれに対して，7Σ，8Σ 状態では，核間距離が大き共鳴ピークの衝突エネルギーと比較できることにい領域では分極力しか働かないため，ポテンシ

なる．これらの作業から電子捕獲断面積に現れャルエネルギー曲線はほぼ水平となっている．

た共鳴ピークが準分子のどんな振動回転量子数

Massey の判別条件[18]や Landau-Zener 公式に対応しているか同定できる．

[18]により，擬交差点でのエネルギー差と衝突エネルギーから，その点で遷移が起こり易いかどう

Copyrigh t© 2012 The Atom ic Colli sion Soci ety of Japan, All ri gh ts reserv ed.

20

小さく，共鳴が現れる極低衝突エネルギーでも電子捕獲断面積は小さいことが予測される．表 2からわかるようにポテンシャル井戸(A)の深さは0.025 au とかなり大きいので，この井戸による共鳴は比較的高い衝突エネルギーで現れると考えられる．一方，井戸(B)の深さは 2 × 10-5 au と浅い．小さい核間距離に現れる井戸(A)はアルカリ金属原子の分極率がそのまま反映されたものであるが，外側の井戸(B)は擬交差点で入射チャン
ネル He2+ + Li と電子捕獲チャンネル He+(3s) +図 4: 準分子 HeLi2+の断熱ポテンシャルエネルギ
ー曲線．

Li+のポテンシャル曲線の交差によるものである．
この擬交差点 R が大きい領域にあり，この点で分極によるポテンシャルが小さいことも井戸が浅い理由である．
②He2+ + Na 衝突系：擬交差点は He2+ + Li 系と比べ R の小さい方に移動している．これは無限遠で He2+ + Na(3s) 状態が He2+ + Li(2s) よりも，
エネルギー的に，それぞれのイオン化状態（He2++ Na+ + e ， He2+ + Li+ + e 状態）に近いことによる．この衝突系でも 6Σ−7Σ 間，5Σ−6Σ 間の擬交差点はほぼ透熱的に通り抜け，電子捕獲は主と

図 5: 準分子 HeNa2+の断熱ポテンシャルエネル
ギー曲線．

して 4Σ−5Σ 間の擬交差点で起こり， He+(3s) +Na+ 状態への電子捕獲が支配的であると予想できる．4Σ−5Σ 間の擬交差点のエネルギー差は，
He2+ + Li 衝突系と比べて大きくなっており，共鳴が現れる極低衝突エネルギー領域の遷移にちょうど良い大きさなので，電子捕獲断面積は大きくなると予測される．井戸(B)の深さも He2+ + Li系のように浅くはない．いろいろな点を総合すると，この系はこれまで扱った系の中では最も高い
衝突エネルギーで共鳴現象が現れると考えられる．
③He2+ + K 衝突系：無限遠での始状態 He2+ +図 6: 準分子 HeK2+の断熱ポテンシャルエネルギ
ー曲線．

K(4s)のエネルギーが前述の 2 つの系よりも更に高いため，擬交差点は 2 つの系よりも更に R の小さい方へ移動している．これに伴い 4Σ−5Σ 間のか見積もることができる．

擬交差点のエネルギー差は大きくなりすぎ，低い①He2+ + Li 衝突系：表 1 によると核間距離 R =衝突エネルギーではこの点で遷移は起こらなく

35−45 au に擬交差点が 3 つあるが， R  45 au

なってしまう．一方，6Σ−7Σ 間の擬交差点ではにある 6Σ−7Σ 間の擬交差点， R  42 au にあるエネルギー差が小さすぎ透熱的に通過すること

5Σ−6Σ 間の擬交差点は透熱的に通過し，おもか ら ， 5Σ−6Σ 間 の 擬 交 差 点 で の 遷 移 に よ るに R  36 au の 4Σ−5Σ 間 の 擬 交 差 点 で++He+(3p) + K+ 状態への電子捕獲断面積だけが



He (3s) + Li 状態への電子捕獲が起こると考え

ある程度の大きさを持つと思われる．共鳴はポテられる．しかし，この点のエネルギー差はかなり

Copyrigh t© 2012 The Atom ic Colli sion Soci ety of Japan, All ri gh ts reserv ed.

21

表 1: 疑交差点の位置とその点でのエネルギー差．
系
状態間
位置(au)ΔE(eV)
He2+ + Li

He2+ + Na

He2+ + K

4−5

35.62

6.0×10-5

5−6

41.56

2.7×10-7

6−7

45.07

2.7×10-8

4−5

24.1

1.6×10-2

5−6

30.64

2.7×10-6

6−7

33.83

1.8×10-7

4−5

18

6.8×10-1

5−6

17.2

2.7×10-3

6−7

20.06

4.9×10-5

表 2: 断熱ポテンシャルの井戸の深さ．R = ∞ での始状態のエネルギーとの差の大きさである．
ポテンシャル井戸の深さ
系
記号
(au)
He2+ + Li

(A)

0.025

(B)

0.00002

He2+ + Na(A)

0.0251(B)

0.0009(A)

0.0397(B)

0.0061

(C)

0.004

2+

He + K

ンシャル井戸(B)，(C)によるものが現れると考え部分波の角運動量量子数が得られるだけなので，られるが，断面積が小さいことからこの系で共鳴

例えば，振動回転準位(55,34)と(4,34)は同じエ

は観測しにくいと予想される．

ネルギーとして示してある．図では井戸(A)による共鳴には*をつけて井戸(B)による共鳴と区別してある（He2+ + Na 系の場合も同様）．井戸(B)より

4.2 共鳴の同定とそのメカニズム
He2+ + Li，Na，K 系に対して，4Σ−7Σ の 4 状

も井戸(A)の場合の方がスペクトル法と量子論的態を用いて量子論的緊密結合法により電子捕獲

緊密結合法による帰属結果が一致している準位

断面積を求め，スペクトル法により振動回転準位の数は多い．同じ回転量子数では，井戸(A)のの同定を行うことにより得られた共鳴に関する知

方が振動量子数の大きい状態になっている．ま見は以下のとおりである．

た井戸(A)が関係するピークは比較的鋭く，井戸

2+

①He

+ Li 衝突系：得られた電子捕獲断面積(B)の場合には幅広いものが多い．これは井戸+のうち，He (3s)状態への電子捕獲断面積を図+(A)の方が井戸(B)の場合よりポテンシャル障壁+

7(a), (b) に示す．He (3p) 状態，He (3d) 状態が高く，幅広いため，井戸に捕まった準分子のへの電子捕獲断面積は計算したすべての衝突

寿命が長いためである．計算結果によると共鳴+エネルギー領域で He (3s) 状態の場合に比べ 2が現れる最大の衝突エネルギーは約 0.29 eV で桁程小さかったので割愛した．また，各共鳴ピーある．図 7(b)に現れている共鳴ピークは，いずれクに対しスペクトル法を用いて行った振動回転準も井戸(A)によるものであり，回転量子数が高い．

位の帰属の結果を本文末表 3 に示した（帰属方共鳴が高い衝突エネルギーで現れるには，回転2+法の説明は He+ Na 衝突系の所で行う）．この量子数が高くなっても有効ポテンシャルの井戸

表には，比較のため，量子論的緊密結合法によ

中に不連続状態が生き続ける必要があるからでって得られた断面積のピークの位置も示してある．

ある．

量子論的緊密結合法によるピークの位置からは
Copyrigh t© 2012 The Atom ic Colli sion Soci ety of Japan, All ri gh ts reserv ed.

22図 7: He2+ + Li → He+(3s) + Li+電子捕獲断面積に現れる共鳴ピーク．括弧内の数字は共鳴状態の振動，回転量子数．*がついているものは 4Σ
の，*がついていないものは 5Σ の井戸による共鳴を表す．(上図，a)，(下図，b)．

②He2+ + Na 衝突系：得られた He+(3s)状態への
電 子 捕 獲 断 面 積 を 図 8(a) ， (b) ， (c) に 示 す ．
He+(3p) 状態，He+(3d) 状態への電子捕獲断面

図 8: He2+ + Na → He+(3s) + Na+において現れる
共鳴ピーク．括弧内の数字は共鳴状態の振動，
回転量子数．*がついているものは 4Σ の，*がついていないものは 4Σ と 5Σ の擬交差点における井戸による共鳴を表す．(上図，a)，(中図，b)，(下図，c)．

積は He+(3s) 状態への電子捕獲断面積より 2 桁
以上小さいので割愛した．図 8(a)，(b)から，衝突エネルギーが 0.05 meV から 3 meV の領域では，
He+(3s) 状 態 へ の 電 子 捕 獲 断 面 積 は(2000−8500)×10-16 cm2 という大きな値になることがわかる．一方，図 8(c)は，共鳴が現れる最も高8(c)のピークのように stückelberg 振動と重なって

い衝突エネルギーでも電子捕獲断面積は 250 ×-16

いないことによる．先に示した B3+ + Li 衝突系の

2

10 cm と大きいが，そのピークは stückelberg 振場合の図 1 ほどはっきりはしていないが，K =動と重なってしまい，観測が難しい可能性を示し

0−65 の領域では|S|2 の K 依存性は衝突エネルギている．各共鳴ピークに対して，散乱行列の 2 乗ーが変わってもほとんど差がないこと，K = 66 のの角運動量量子数依存性および衝突エネルギ

場合だけ|S|2 の K 依存性が衝突エネルギーに大ー依存性を調べた．一例として，図 8(b) 中の衝

きく依存していることがわかった．次に，この共鳴

突エネルギーE = 2.5229 meV のピーク (1) を取状態の振動量子数をスペクトル法によって求めり上げ，簡単に説明する．このピークを選んだ理た．まず HeNa2+系の断熱ポテンシャルの 6Σ−7Σ

由は，比較的高い衝突エネルギーで現れ，図Copyrigh t© 2012 The Atom ic Colli sion Soci ety of Japan, All ri gh ts reserv ed.

23図 9: 準 分 子 HeNa2+ の 断 熱 ポ テ ン シ ャ ル の6Σ−7Σ 間，5Σ−6Σ 間の擬似交差点を透熱的に通過させたポテンシャルエネルギーに K = 66 の遠心力ポテンシャルを加えた有効ポテンシャル中の振動状態のスペクトルパワー．

図 10: 準分子 HeNa2+ の 断熱ポ テン シ ャ ルの6Σ−7Σ 間，5Σ−6Σ 間の擬交差点を透熱的に通過させたポテンシャルエネルギーに K = 66 の遠心力ポテンシャルを加えた有効ポテンシャル中の振動状態  = 5 の固有関数の絶対値の二乗．

間，5Σ−6Σ 間の擬交差点を透熱的に通過させたポテンシャルエネルギーに K = 66 の遠心力ポ断面積を図 11 に示す．He+(3s) 状態，He+(3d)テンシャルを加えた有効ポテンシャルエネルギ

状態への電子捕獲断面積は He+(3p) 状態への電子捕獲過程に比べ 2 桁程小さいので割愛した．
断熱ポテンシャルの 6Σ−7Σ 間の擬交差点を透
熱的に通過させた場合の有効ポテンシャル（ポテンシャル I）および 6Σ−7Σ 間に加え 5Σ−6Σ 間の擬交差点も透熱的に通過させた有効ポテンシ
ャル（ポテンシャル II）の 2 つを用いての共鳴ピー
クの帰属結果を本文末表 5 にまとめるとともに図11 の各ピークに記した．図中の*印はポテンシャ
ル II を用いての帰属であることを示す．共鳴のピークは，He2+ + Li 衝突系と同様鋭く，比較的等間隔に現れている．共鳴ピークの位置，形が単調なことから，共鳴ピークについていろいろな知
見を得ることができる．その特徴は(1) 衝突エネ
ルギーが増加するにつれベースとなっている断面積の値は増加している．(2) 振動回転量子数
が(30，63)から(23，84)までの 8 個の鋭いピークと振動回転量子数が(30，64)から(23，82)までの 7個の少し幅広いピークが，ともに，ほぼ等間隔で現れている．両グループとも，それぞれ，順に振動量子数が 1 減少，回転量子数が 3 増加している．振動エネルギー1 個分と回転エネルギー3 個分にそれほど大きい差がないため，有効ポテン
シャルの形，井戸の底から計ったそれらの状態のエネルギーなどが似ているものと思われる．(3)振動量子数が同じで回転量子数が 1 異なる(24，81)のピークと(24，82)のピークを比べると，(24，ー（ポテンシャル I）を用いて「スペクトルパワー」を求めた．その結果得られた 4 つのピークを図 9に示した．そのうち E = 0.00009179 au (= 2.4977
meV)のピークのエネルギーが図 8(b)の E =2.5229 meV のピークのエネルギーに近い．そこで，更にスペクトル法を用いてポテンシャル I が作る井戸中の，E = 0.00009179 au のピークに対応する，準安定状態の固有関数を求めた．得られた固有関数の二乗を有効ポテンシャルとともに
図 10 に示す．有効ポテンシャルには R = 50 au
付近にポテンシャル障壁があり，固有関数はこの障壁の内側に捕まっている振動量子数 5 の状態であることがわかる．他の共鳴ピークに対しても同様な方法を用いて振動回転量子数を求めた．
得られた共鳴ピークの帰属結果を本文末表 4 にまとめるとともに図 8(a)，(b)，(c) の各ピークに記した．但し，*印を付けたピークは断熱ポテンシャ
ルの 6Σ−7Σ 間，5Σ−6Σ 間に加え，4Σ−5Σ 間の擬交差点も透熱的に通過させた有効ポテンシャ
ル（ポテンシャル II）を用いた結果である． He2+
+ Li 衝突系の場合とは異なり，スペクトル法によって同定できたピークには井戸(B)により生じるピークが多い．計算結果によるとこの衝突系で共鳴が現れる最大の衝突エネルギーは約 0.3 eVである．
③He2+ + K 衝突系：衝突エネルギーが 0.5 meVから 3 meV の領域の He+(3p)状態への電子捕獲

Copyrigh t© 2012 The Atom ic Colli sion Soci ety of Japan, All ri gh ts reserv ed.

24

図 11: He2+ + K → He+(3p) + K+において現れる
共鳴ピーク．括弧内の数字は共鳴状態の振動，
回転量子数．*は 4Σ の井戸による共鳴を表す．

82)のピークの方が高さは低いが幅広い．これは(24，82)のピークの方が有効ポテンシャル障壁の厚さが薄く，寿命が短いことによる．計算結果に

図 12: He2+ + Na 衝突系に対して，ポテンシャル
II 上を衝突エネルギー2.4977meV で波束を走らせた場合の核波束の時間変化．

よると共鳴が起こる最大の衝突エネルギーは約0.044 eV と小さく，この系は He2+ + Na 系と比べ共鳴は観測しにくいと思われる．しかし，この系は数十 meV 以上の衝突エネルギーでの電子捕から，この波束はポテンシャル井戸にトラップされ獲断面積の計算が難しいため，実際はより高いた振動量子数が 5 の振動回転状態に対応してい衝突エネルギーで共鳴が現れているにもかかわることがわかる．波束の確率の減少は，波束がポらず見逃している可能性もある．

テンシャル障壁を通り抜けて外に逃げ出していることに対応している．
波束法により共鳴状態にトラップされた波束を

4.3 多チャンネル核波束法による共鳴ダイナ

可視化することができた．可視化により，電子捕

ミクスの解析
共鳴を観測するのに最適である可能性が高い

獲断面積のピークの位置で共鳴が起こっている

He2+ + Na 系に対して，スペクトル法により振動回こと，共鳴のメカニズムと動力学をよりはっきりと転状態(5, 66)であると帰属されたピークを選び，理解することが可能になった．

衝突エネルギーE = 2.4977 meV で核波束を走ら
4.4 まとめせることにより，共鳴のダイナミクスを解析した．
得られた計算結果を図 12 に示す．この図より，入Li, Na, K と原子番号が大きくなるにつれ，イ

射波束は，衝突により，3 つに分かれることがわ

オン化エネルギーが小さくなる（Li, Na, K のイオ

かる．1 つ目は R =50 au 付近にあるポテンシャ

ン化エネルギーはそれぞれ 5.39 eV, 5.14 eV,

ル障壁によって跳ね返されたものであり，2 つ目

4.34 eV）とともに双極子分極率は大きくなる．そはより内側の R ≦ 25 au にある He2+イオンとの結果，始状態（He2+ + Li, Na, K）と電子捕獲状

Li+ イオンの反発ポテンシャルによって跳ね返さ態（He+(3l) + Li+, Na+, K+）のポテンシャルエネル

れたものである．2 つの波束の大きさを比べると，

ギー曲線の疑交差点は内側に移動し井戸が深く

後者の波束の確率の方が大きいことがわかる．3

なるので，この順番で共鳴が高い衝突エネルギ

つ目は，R = 25 au まで侵入し，時間が経過してーで現れると予想していた．しかし，計算結果は，も同じ核間距離にとどまっている波束である．こ

最も高い衝突エネルギーで共鳴が現れるのはの波束は時間と共に徐々にその確率は減少してHe2+ + Na 衝突系であることを示した．これにはいるが同じ形を保っており，6 つの山からできて非断熱遷移確率という要因がからんでいる．疑いる．図 10 に示した固有関数と同じ形であること

交差点の位置とその位置での断熱状態間のエ

Copyrigh t© 2012 The Atom ic Colli sion Soci ety of Japan, All ri gh ts reserv ed.

25

ネルギー差にはある程度関連があり，小さい核

れらから，より高い衝突エネルギーで共鳴を観測

間距離にある擬交差点ほどエネルギー差は大きするにはどのような系を選べば良いか予測が可い．エネルギー差が大きいと低エネルギーでは能になる．今回述べた衝突系の中では He2+ +遷移が起こりにくくなり，またエネルギー差が小さ

Na 衝突系が最も高い衝突エネルギーで共鳴がいとほとんどの衝突エネルギーで非断熱遷移確

観測されうる系であることがわかったが，(1) 標的

率が 1 になってしまう．衝突対が擬交差点を通過原子の原子番号，(2) 入射イオンの原子番号，する際にほとんど遷移してしまうと，衝突では「行(3) 入射イオンの電荷をうまく調整することで，よき」と「帰り」の 2 度必ず交差点を通過するため，

り適切な衝突系を見つけることも可能になる．
理論的に，十分正確に，共鳴現象が起こる衝

衝突が終わった後反応が起こっていないというこ
2++ Na 衝突系では電子突エネルギーがわかったとする．そのとき共鳴を捕獲過程の擬交差点が，高い衝突エネルギー

観測するための実験的条件としてどのようなこと

で共鳴が生じる適当な位置にたまたまあった．

が要求されるだろうか．これまで報告されているとになってしまう． He

これまでは，1 番内側にある 1 番深いポテンシ

共鳴はすべて衝突エネルギー1 eV 以下に現れ

ャルの井戸もしくはその外側の擬交差点での井ている．この衝突エネルギーで多価イオンを標的

戸に関連させて共鳴について議論してきた．そに衝突させるという条件は十分クリヤーできる．

の他に，もっと大きい核間距離にも疑交差点があ次に要求される条件はビームのエネルギー幅をる．これらの交差点において生じるポテンシャル

十分狭くするということである．共鳴ピークを観測

井戸は内側のポテンシャル井戸よりは浅いので，するにはビームのエネルギー幅が共鳴ピークのそれらに起因する共鳴は衝突エネルギーの低い

エネルギー幅よりも十分狭い必要がある．我々のところに現れるはずである．しかし，交差点が内計算結果によると，高い衝突エネルギーで現れ側に移動するとそれらの井戸も深くなるので，共る場合でも，共鳴ピークのエネルギー幅は 10-3

鳴が高い衝突エネルギーで現れる可能性は捨

eV 以下である．このような場合，少なくともビームてきれない．

のエネルギー幅を 10-3 eV 以下に，更に衝突エネルギーを 10-3 eV 以下の間隔で変化させなけ

5. 共鳴観測に向けてればならない．

より高い衝突エネルギーで，電子捕獲過程に共鳴状態の寿命が衝突時間よりも短いと共鳴

おいて，共鳴を観測するには，標的原子の分極状態とは言えない．仮に共鳴が衝突エネルギー率が大きく，入射イオンの電荷が大きいことが好10-1 eV で現れ，その共鳴ピークのエネルギー幅

ましい．しかし，疑交差点の位置，その点でのエが 10-3 eV で，ビームのエネルギー幅を 10-4 eV

ネルギー差に注意する必要があることがわかっにすることができ，衝突エネルギーを 10-4 eV のた．また，一般的にはポテンシャル井戸をもつ状

間隔で変えることができたとする．衝突対の換算

態よりもエネルギー的に下に他の状態が存在す

質量を 1 u とすると，距離 1 Å を通過するのに 2

るので，ポテンシャル井戸が内側に移動し深くな

× 10-14 sec かかる．共鳴ピークのエネルギー幅かると，共鳴を引き起こす状態とそれら下の状態とら寿命は 6 × 10-13 sec であるので，このような状が相互作用して井戸が浅くなることにも注意が必況下で，共鳴は十分観測できることがわかる．実要である．

験的に測定できるのであれば，共鳴ピークの幅

疑交差点の位置は，アルカリ金属原子が標的が狭ければ狭い程寿命が長いので測定には有の場合，その原子番号が大きくなるにつれ小さ利と言える．

い方へ移動する．一方，「低速多価イオン・原子

共鳴状態の寿命を求める理論的方法は 2 つあ

衝突における電子捕獲過程と共鳴現象 その 2」

る．1 つは共鳴状態にある波束の確率の減少かで扱うように，入射イオンの原子番号が大きくな

ら求める方法であり，もう 1 つは断面積のエネル

るとともにまた入射イオンの電荷が小さくなるとと

ギー依存性に現れるピークの半値幅から求める

もに疑交差点の位置は大きい方へ移動する．こ方法である．2 つの方法を用いて計算を行ったとCopyrigh t© 2012 The Atom ic Colli sion Soci ety of Japan, All ri gh ts reserv ed.

26

ころ，これまで扱った系では，共鳴状態の寿命は[5] J. P. Toennies, W. Welz, and G. Wolf, J.

0.5−18 psec の範囲にあることがわかった．衝突

Chem. Phys. 72, 614 (1979).

領域 1 au，衝突エネルギー1 meV − 1 eV の場合[6] J. P. Davis and W. R. Thorson, Can. J. Phys.

の衝突時間は 4×10-3 − 0.1 psec であるので，共56, 996 (1978).

鳴状態の寿命は衝突時間より長く，共鳴状態は[7] M. Rittby, N. Elander, E. Brandas, and A.

観測可能なはずである．

Barany, J. Phys. B17, L677 (1984).
[8] N. Shimakura and M. Kimura, Phys. Rev.

本総説は九州大学理学研究院の故季村夆生

A44, 1659 (1991).

教授との共同研究の思い出として引き受けた．

[9] S. Suzuki, N. Shimakura, T. Shirai, and M.

二人で形状共鳴についていくつか論文を書き，Kimura, J. Phys. B17, 1741 (1998).

そのあと総説をまとめるという計画を立てていた[10] N. Suzuki, N. Shimakura, and H. Kono,

が，それも今は夢となった．この総説は二人で書Physica Scripta T92, 435 (2001).

きかけていた論文の一部をまとめ，総説用に書き[11] D. W. Chandler, J. Chem. Phys. 132, 110901

直したものである．研究遂行にあたり，新潟大学(2010).

自然科学研究科の 2 人の修了生，桜井智君，鈴木伸啓君に協力していただいた．謝意を表した[12] J. N. Bardsley, Case Stud. At. Phys. 4, 299

い．

(1974).
[13] H. G. Heil, S. E. Butler, and A. Dalgarno,

参考文献

Phys. Rev. A44, 1659 (1991).

[1] S. Kravis, H. Saitoh, K. Okuno, K. Soejima,[14] R. Boyd, T-S Ho, H. Rabitz, D. A.

M. Kimura, I. Shimamura, Y. Awaya, Y.

Padmavathi, and M. K. Mishra, J. Chem.

Kaneko, M. Oura, and N. Shimakura, Phys.

Phys. 101, 2023 (1994).

Rev. A52 (2), 1206 (1995).

[15] N. Balakrishnan, C. Kalyanaraman, and N.

[2] L. Folkerts, M. A. Haque, C. C. Havener, N.

Sathyamurth, Phys. Rep. 280, 79 (1997).

Shimakura, and M. Kimura, Phys. Rev. A51[16] N. Vaeck, M. D. Lecomte, and J. Liévin, J.

(5), 3685 (1995).

Phys. B32, 409 (1999).

[3] L. Bolzmann, Vorlesung Uber Gastheorie II[17] M. D. Feit, J. A. Fleck, and A. Steiger, J.

(J. A. Barth, Leipzig, 1989)

Comput. Phys. 47, 412 (1982).
[18] 例えば，高柳和夫，電子・原子・分子の衝突[4] A. Schutte, G. Scols, F. Tommasini, and D.

（培風館，1972，改訂版，1996）

Bassi, Phys. Rev. Lett. 29, 979 (1972).

Copyrigh t© 2012 The Atom ic Colli sion Soci ety of Japan, All ri gh ts reserv ed.

27

表 3: He2+ + Li 衝突系における振動回転準位の同定．
井戸(A)による共鳴
井戸(B)による共鳴
振動回転準位

緊密結合法
(meV)

スペクトル法(meV)

振動回転準位

緊密結合法
(meV)

スペクトル法
(meV)

55,34

0.370520.171434,340.37052

0.3646254,36

0.393080.29932

3,36

0.39308

0.39728

53,38

0.463

0.473474,35

0.425110.3210952,40

0.533420.51701

3,38

0.463

0.3646353,39

0.570450.378232,40

0.533420.5687152,41

0.67119

0.666673,39

0.570450.4653151,43

0.776290.764632,41

0.671190.6693950,45

0.880280.868032,430.776290.8789251,44

0.92880.661230,450.880280.5986449,47

0.980910.97962,440.92880.7401449,48

1.21371

1.208171,470.98091

1.1591948,50

1.35727

1.34151,48

1.21371

1.2979649,49

1.41903

1.417690,50

1.35727

1.1700747,52

1.49397

1.480281,49

1.41903

1.4476348,51

1.60627

1.61633

0,52

1.49397

1.5646347,53

1.83082

1.79593

1,51

1.60627

1.3659946,55

1.98696

1.975520,531.83082

1.7632745,57

2.16127

2.14967

0,55

1.98696

2.1632746,56

2.31476

2.30205

0,57

2.16127

2.5877745,58

2.55926

2.544230,56

2.31476

2.35375

44,60

2.78932.778240,58

2.55926

2.808181,202

265.144265.933

0,60

2.78933.341510,204

270.548271.402K=202

265.144−0,205

276.796277.661K=204

270.548−0,206

282.94283.784K=205276.796−0,207

289.005289.852

K=206282.94−

Copyrigh t© 2012 The Atom ic Colli sion Soci ety of Japan, All ri gh ts reserv ed.

28

表 4: He2+ + Na 衝突系における振動回転準位の同定．
井戸(A)による共鳴
井戸(B)による共鳴
振動回転準位

緊密結合法
(meV)

スペクトル法(meV)

振動回転準位

緊密結合法
(meV)

スペクトル法
(meV)

K= 19

0.05216−

16,19

0.052160.03265

K = 220.05733

0.0761915,220.05733

0.04082

K = 200.06959−

16,20

0.06959

0.0462671,23

0.07798

0.0952415,230.07798

0.06531

K = 210.08856

0.1088416,21

0.08856

0.06258

K = 26

0.097530.05442

K = 26

0.09753−

K = 270.12188−15,27

0.121880.13061

K = 280.14726−13,280.14726

0.15238

K = 290.17340.155113,29

0.17340.17922

K = 340.18494

0.103413,34

0.184940.13877

K = 330.18938

0.141514,33

0.189380.20136

K = 350.21382

0.1632713,35

0.21382

0.1741566,360.225610.2231313,360.225610.22585

K = 350.24517−14,35

0.245170.26122

K = 36

0.273110.1142914,36

0.273110.23946

K = 390.31746−12,39

0.317460.31565

K = 400.357450.29932

12,40

0.357450.244964,410.390260.38095

11,41

0.39026

0.3238162,45

0.52253−

10,45

0.52253

0.5088561,48

0.67838−9,48

0.67838

0.631361,50

0.863960.906139,50

0.863960.851760,51

0.93069

0.906129,51

0.93069

0.9659959,53

1.0475−8,53

1.04751.0313

K = 56

1.42177−

7,56

1.42177

1.3823257,58

1.53453−7,58

1.53453

1.5428655,60

1.64834−

6,60

1.64834

1.6217856,62

2.01335

2.048997,622.01335

1.9972955,632.18766

1.983686,632.18766

2.1932155,65

2.31685

2.141515,652.31685

2.2911753,66

2.523962.440825,662.52292.497753,69

2.92992

2.966

4,69

2.92992

2.82994

1,244

273.332

273.226

K = 244273.332−

1,245

278.216278.069

K = 245278.216−1,246

283.24282.831

K = 246283.24−0,248

288.304

288.164

K = 248288.304−

0,249

293.528293.171

K = 249

293.528−0,250

298.205298151

K = 250

298.205−

0,251

302.695

302.994

K = 251302.695−

0,252

307.882308.056

K = 252

307.882−0,253

312.168312.763

K = 253

312.168−

Copyrigh t© 2012 The Atom ic Colli sion Soci ety of Japan, All ri gh ts reserv ed.

29

表 5: He2+ + K 衝突系における振動回転準位の同定．
井戸(B)による共鳴
井戸(C)による共鳴
振動回転準位

緊密結合法
(meV)

スペクトル法(meV)

振動回転準位

緊密結合法
(meV)

スペクトル法
(meV)

K=46
K=45

0.32223

0.2585K=460.32223

0.34286

0.337850.31837K=450.33785

0.29116K=480.35320.35646K=480.35320.2884436,47

0.36169

0.4

36,47

0.361690.3078435,49

0.39333

0.2911635,49

0.39333

0.41361K=510.43206

0.31837

K=510.43206

0.4381K=530.47698

0.42993K=53

0.47698

0.42993K=550.52433

0.46803

K=55

0.52433

0.59592

32,57

0.56958

0.503432,57

0.56958

0.644933,56

0.59298

0.55234

K=560.59298

0.55238

32,58

0.65883

0.65306

K=58

0.65883

0.6013631,60

0.71671

0.71293

K=600.716710.644931,61

0.81666

0.8136131,61

0.81666

0.7564730,63

0.88686

0.8843629,63

0.886860.8054529,66

1.08241.08028

29,66

1.08240.9877629,67

1.20561

1.2081729,671.20561

1.1292628,69

1.302211.3034128,69

1.302211.19456

29,70

1.439791.4503528,70

1.43979

1.3523927,72

1.54768

1.5564727,72

1.54768

1.65443

26,75

1.82384

1.8367427,75

1.82384

1.9374226,76

1.99495

2.0108926,76

1.99495

1.8802825,78

2.13498

2.1496725,78

2.134982.2585125,79

2.322192.3455925,79

2.322192.46804

24,81

2.481372.49253

24,81

2.481372.6049324,82

2.687522.7102224,82

2.68752

2.8272223,842.86015

2.8680423,84

2.86015

2.98233

3,164

38.3655

38.3539

K=164

38.2655−

3,165

39.570239.6603

K=165

39.5702−2,170

44.0152

44.1226

K=170

44.0152−

Copyrigh t© 2012 The Atom ic Colli sion Soci ety of Japan, All ri gh ts reserv ed.

30

解

説

強い短距離斥力相関の下での弱結合少数多体系− エフィモフ物理と 4 He 原子の 3,4 クラスター系 −(後編)
上村正康, 肥山詠美子
理化学研究所 仁科センター

ストレンジネス核物理研究室

mkamimura@riken.jp, hiyama@riken.jp
平成 24 年 7 月 23 日原稿受付

前編（9 巻 4 号 20 ページ）に引き続き, 本稿では表題に関する解説を行う. 4 He 原子のクラスターの理論研究は エフィモフ物理や冷却気体原子の物理の基礎研究に重要な知見を提供してきた．最近，

4 体系のエフィモフ物理も注目されている．本稿では, 4 He 原子の 4 体系 (tetramer) の基底・励起状態について, 4 He-4 He の realistic potential を用いた精密計算結果を基に, エネルギーレベル, 空間
構造, 短距離相関, 3 体・4 体の束縛エネルギーの相関 (generalized Tjon lines) などについて議論する. 「大きな散乱長で短距離相互作用をしている少数粒子系」の持つ universality という観点から,
原子核の弱結合 3 体・4 体系との類似性にも着目する.

1. はじめに1.1 背後に何があるのか?
筆者の 1 人 (上村) は, 核物理の弱結合系に関する研究会 (2012 年 7 月) において, 本解説 (前・後編) に相当する内容の講演を行った. その冒頭において, 以下の図 1-3 の 3 枚のスライドを示して, 原子核の弱結合系と 4 He 原子のクラス
ター系との間に有る共通性について注意を喚起
した.
図 1 は, 「このエネルギー レベル (模式図) は,何という粒子 (X) の 2 体, 3 体, 4 体系のエネル
ギーレベルか?」という質問形であるが, 答は,

図 1:

X= He 原子である. しかし, 核物理研究者から
4

見ると, この図は, X=α 粒子（4 He 原子核）として, 2 体系=8 Be 原子核, 3 体系＝ 12 C 原子核,

4 体系＝ 16 O 原子核のレベルに特徴がそっくり

4

He 原子の 2, 3, 4 体系のレベル構造の模式図. 4 He 原子核 (α 粒子）の 2, 3, 4 体系
(8 Be, 12 C, 16 O 核) のレベル構造と特徴がそっくりである. また, 核子の 2, 3, 4 体系
（d, t, α) ともよく似ている（但し, t の励起状態を除く）.

である. また, X=核子として, 2 体系=重陽子,

3 体系＝三重陽子, 4 体系＝ α 粒子とも似ている

図 2 は, 3 体系と 4 体系の波動関数の overlap（但し, 三重陽子の励起状態は束縛状態, 共鳴状function と呼ばれるもののグラフである. 左が態としては見つかっていないが）. この類似性

4の背後に何があるのだろうか.

He 原子核の基底・励起状態 (0+ ) と三重陽子との overlap [1] であり, 右が 4 He 原子の tetramerc 2012 The Atomic Collision Society of Japan, All rights reserved.
Copyright⃝31図 3: 左は, 核物理で Tjon line と呼ばれるもので,様々な 2 体核力による三重陽子と α 粒子の結合エネルギーが, ほぼ直線の相関になっている [3, 4]. 右図では, 4 He 原子の trimerと tetramer の基底状態の結合エネルギーの相関 (Atomic Tjon line) が様々な realistic
4
He-4 He potentials に関して, ほぼ完璧に直線になっていることを表す [2].

図 2: 3 体系と 4 体系の波動関数の overlap function. 左が 4 He 原子核の基底・励起状態 (0+ )と三重陽子との overlap [1], 右が 4 He 原子のtetramer の基底・励起状態 (0+ ) と trimerとの overlap [2]. 両者のスケールは数十万倍違うが, 空間構造は非常によく似ている.

の基底・励起状態 (0+ ) と trimer との overlap [2]である. 両者のスケールは数十万倍の違いがあa large scattering length have universal low-

るが, この空間構造の類似性の背後に何がある

energy properties that do not depend on theのだろうか.

details of their intrinsic structure or their in-

図 3 は, 3 体系と 4 体系の結合エネルギーの相

teractions at short distances.” ということであ関を表す. 左図は 核物理で Tjon line [3, 4] とろう.

呼ばれるもので, 様々な 2 体核力による三重陽

1.3 Eﬀective theory

子と α 粒子の結合エネルギーが, ほぼ直線の相エフィモフ物理, 言い換えれば Universal Few-

関になっていることを示している. 右図は 4 He

body Physics, における 3 体系, 4 体系の理論計

原子の trimer と tetramer の基底状態の結合エ

算には 2 つのタイプがある.

ネルギーの相関を表す. 筆者等の論文 [2] におい(a) Eﬀective theory の leading order 計算,て, それが 様々な realistic 4 He-4 He potentials(b) 原子間の realistic potential を用いる計算.

に関して, ほぼ完璧に直線に乗っていることが筆者等の勝手な理解では, 前者 (a) の立場は次明らかになった（§2.5). 直線の傾き (slope) は,のように見える— “散乱長が十分大きければ,左が 4.8, 右が 4.778 である. この一致の背後に短距離相互作用の詳細に依らない” のだから難何があるのだろうか.

しい realistic potential を使う必要は無く, 簡

1.2 背後に universality

単な eﬀective potential でよい. i) 2 体の con-

上記研究会で述べたことだが, 背後にあるもtact (point-like) inteteraction で散乱長を出し,の（即ち, エフィモフ物理の核心）を一言でいうii) 3-body parameter (3 体の contact inteter-と（著名なレビュー論文 [5] の表題を取って）,action または 3 体の内部領域境界条件) によっ“Universality in few-body systems with largeて 3 体・4 体のエネルギー・波動関数を出す—

scattering length” であろう. もう少し引き伸ばという様に, realistic-potential 3 体・4 体問題

すと (この論文の abstract の冒頭文から取って）を renormalize できるはずである [5]. (本稿で“Particles with short-range interactions andは, 相互作用の中身について, これ以上立ち入ら

c 2012 The Atomic Collision Society of Japan, All rights reserved.
Copyright⃝

32ない). universal property は leading-order ef-

一 方,realistic

potential を 用 い た 4 He

fective theory で理解でき, それを超えることは

tetramer の計算は, 極度に強い短距離斥力コア

non-universal eﬀect である.

の下での非常に浅い束縛励起状態を非断熱的に

上記 (b) の realistic potential は, 目下のとこ

解く, という量子力学的 4 体問題の難しさによろ He- He potential しかない. その他のケー

り, 進展が遅れていた. 即ち, 2006 年の時点で,ス (Li, K, Rb, Cs など）は, 作るのが難しい. し基底状態のエネルギー (LM2M2 potential) は,かし, He 3 体・4 体の結果は, scaling によって数件の計算がほぼ一致した値を出していたが大きな原子にも敷衍できるであろう. Eﬀective（有効数字 3 桁目が僅かに異なる）, 励起状態

theory による分析の基礎を固める, または, 十については, trimer+atom の threshold energy

分記述できない部分 (non-universal eﬀect) を明(−126.40 mK) から測った束縛エネルギーは,らかにするという役割もある.

1.0 mK,

4

4

4

6.6 mK,

51.6 mK というように,バラついていた（文献 [12] の Table VIII 参

1.4 4 体のエフィモフ状態照）. 特に, 1.0 mK という値は, Faddeev 法が前編で述べたように, 2006 年に 3 体のエフィ

低エネルギー散乱計算の S 行列の振る舞いかモフ状態 [6, 7] は, Cs の超冷却原子気体においてら求めたものであり, 同法としては, “束縛状態観測された（その後, Rb, K, Li などで）. 散乱長計算” によって このように極度に強い短距離を変化させて, エフィモフ状態が 3 体 breakup

斥力コアの下での非常に浅い束縛励起状態を

threshold のすぐ上の共鳴状態になる散乱長の解く事は現在の計算技術では不可能である, と値（leading-order eﬀective theory が予測する

give-up 宣言をしていた [14].

値）で, 3-body recombination と呼ばれる現象

しかし, 上述したように, leading-order eﬀec-(B+B+B → B2 +B) が観測され（B はボソン),

tive theory 計算が予測した通りに 4 体のエフィ

エフィモフ状態の証拠とみなされた. その後,モフ状態の証拠が見つかったという刺激もあり,2009 年に 4 体ボソンのエフィモフ状態についrealistic potential を用いた 4 He tetramer 計算て, leading-order eﬀective theory の計算がなが渇望されていた. また, eﬀective theory 計算

され [8], それが予測する散乱長で, 4-body re-が, 「ボソン 3 体の基底状態および各励起状態のcombination (B+B+B+B→B3 +B, B2 +B+B,

エネルギーのすぐ下に, ボソン 4 体の束縛状態

B2 +B2 ) が観測され [9], 4 体のエフィモフ状態

が 2 個付随している」という明解な予言 [8] をの証拠とされた.

しており, これを 4 He-4 He realistic potential を

1.5 Realistic-potential 4 体計算

使って確認することも急務であった. このこと

realistic potential を用いた 4 He trimer の計を筆者等が知って, 4 He trimer, tetramer の精密

算は, 前編で述べたように多くの文献で実行さ計算を行い, 文献 [12] に発表した次第である.

れ, 十分一致する結果が得られた. 最近の重要

1.6 3 体・4 体の束縛エネルギーの相関な新しい計算は, P. Naidon, 上田正仁と筆者の散乱長の大きい 2 体短距離相互作用をしてい1 人（肥山） [10] による次の計算であろう —4る 3 体系・4 体系が持つ universality の 1 つに,He-4 He の LM2M2 realistic potential [11] をそれぞれの束縛エネルギーの間の相関がある.

用いて, 「4 He trimer の基底状態が 3 体 breakup

図 3 で見たように, 原子核物理では, 3 He 原子核

threshold のすぐ上の共鳴状態になる散乱長のと 4 He 原子核の基底状態の束縛エネルギー（B3値」を求め, 実験で見つかった Cs, Rb, K, Li のと B4 ) を様々な realistic 2 体核力 potential で計

3-body recombination を起こす散乱長と同じで算して, B3 -B4 平面に plot したもの（ほぼ直線あることを示した（potential の eﬀective rangeになる, 傾き=4.8）が, Tjon line として知られで normalize する）.

c 2012 The Atomic Collision Society of Japan, All rights reserved.
Copyright⃝

33

ている. 4 He 原子の 3 体系と 4 体系には, どちら

2. ヘリウム原子 4 体系 4 He tetramer

も基底状態と励起状態があるので, その束縛エ
ネルギー (B3(0) , B3(1) と B4(0) , B4(1) ) の間には, 計 6

少 数 粒 子 系 の 精 密 解 法「 ガ ウ ス 関 数 展 開通りの相関がある（generalized Tjon lines と呼法（GEM）」によってヘリウム原子の 4 体系ばれる）.

tetramer を解く. まず, 解き方を説明する（詳それらの相関は, leading-order eﬀective the-

しくは, 文献 [12] の Sec.IIIA）. ヤコビ座標を

ory を用いて, 文献 [13] で計算され (2004 年）, 妥

図 4 に示す.

当なエネルギー領域でどれも直線の相関になっていることが報告されていた（B3(0) -B4(0) 相関の傾き 4.08 と核力の場合の傾き 4.8 との違いは, 前者に next-to-leading order eﬀect=non-universal

eﬀect が入っていないためか？ §2.5 で明らかになる）. しかし, 4 He の realistic potential を用いた計算とどのような相違があるか（両者の差として, non-universal eﬀect が見えるか）は, 不明であった — 上述したように, realistic potentialによる 4 体計算が不十分であったためである.

図 4: 4 粒子系のヤコビ座標. K-type と H-typeがある. 4 粒子を cyclic に入れ替えたセッ
ト i = 1, ..., 12 (K-type) と i = 13, ..., 18(H-type) とがある. 18 個のあからさまな図
は, 文献 [16] の Fig.18 にある.

筆者等は, この問題を解くため, 多数の 4 He4

He realistic potential を用いて, エネルギー相関の計算を行い, 6 通りの相関が全てほぼ完璧な直線の相関であることを示した [2]. 同時に,

eﬀective theory の leading-order 計算では取り
この 4 体系の Schrödiner 方程式

込めない non-universal eﬀect を明らかにした.

(H − E)Ψ4 = 0,

本稿（後編) では, まず, 4 He 原子の 4 体系

(tetramer) の基底状態と励起状態の 2 つの束縛

H =−状態の解き方と計算結果の説明を行う. また, 前編で紹介した筆者等の「直感モデル：dimerlike-

(1)

ℏ2 2
ℏ2 2
ℏ2 2
∇x −
∇y −
∇
2µx
2µy
2µz z
4
∑

+

(2)

V (rij ),1=i<jpair model」による励起状態エネルギーの予言を満たす波動関数 Ψ4 を求める.

が 精密に的中していることを示す. 使用する

µx =

1
2 m,

µy = 23 m, µz = 34 m, m は 4 He atom の mass.

realistic 4 He potential は, 前編では広く使われV (rij ) は 粒子 i, j 間の中心力 potential であり,てきた LM2M2 potential (1991 年）だけであっ通常, realistic 4 He-4 He potential が使われる.

たが, 後編では, 最近 (2010 年）発表された最も”realistic”とは, 合計 4 個の電子の自由度を然

精密な potential (PCKLJS potential [15] と呼るべき方法で消去して, 4 He 原子核間のみの関ばれる. 相対論補正, QED 補正あり）を含む多数

数として, 4 He 原子間相互作用を表すことであの realistic potentials を用いて, 3 体, 4 体系の束る.

縛エネルギーの間の 6 通りの相関 (Genelarized

4

He trimer, teramer の計算に最も頻繁に使われるのは, 前編の図 2a の赤点線で示した

Tjon lines) を計算し, 相関の universality, non-LM2M2 potential [11] である. V (rij ) の遠方はuniversal eﬀect について議論する.

−6
van der Waals 型 (∝ rij) であるが, 原点近傍

(rij ≲ 2Å) では, 極度に強い斥力になっている(V (0) ∼ 106 K).

c 2012 The Atomic Collision Society of Japan, All rights reserved.
Copyright⃝34Λ = lz = even である.

2.1 4 体系の基底関数

基底関数の動径部分は, 3 体系の場合と同様に,ヘリウム原子は, spinless boson であるから,
Ψ4 を, 対称化された L2 積分可能の基底関数 (K -次のように表す.

type と H -type がある）で展開する.

ϕnsin(x) = xlx e−(x/xnx ) ×
x lx

αmax
K

Ψ4 =

∑

(cos )

αmax
H(sym;K)A(K)+αK ΦαKαK =1

Φ(sym;K)
=
αK

∑

(10)

ψny ly (y) = y ly e−(y/yny ) ,(11)

φnz lz (z) = z lz e−(z/znz ) .

(12)

2

(4)

i=1

Φ(sym;H)=αHsin ω(x/xnx )2 ,ガウスのサイズは等比数列に取る：

Φ(K)αK (xi , yi , zi ),18∑{ cos ω(x/xnx )22(sym;H)A(H), (3)αH Φ αH

αH =1
12
∑

2

Φ(H)
αH (xi , yi , zi ).

(5)

i=13

xnx = x1 anx x −1

(nx = 1, ..., nmax
),
x

(13)

yny = y1 any y −1
znz = z1 anz z −1

(ny = 1, ..., nmax),ymax
(nz = 1, ..., nz ) .

(14)
(15)

ここで, (xi , yi , zi ) は, i 番目のヤコビ座標セット

4 体系の基底関数の総数は 29 056 個である.

である. K -type と H -type の基底関数で展開す基底関数の non-linear paramers は, 文献 [12] のることは重要必須であり, 関数空間を大きく広Table VI に載せてある.

げている.

2.2

(H)
固有エネルギー E と展開係数 A(K)αK (AαH ) は,4

He tetramer のエネルギー準位

表 1 に,4 He tetramer の基底・励起状態の束縛

Raileigh-Ritz の変分法で解く：

エネルギー, 2 粒子間平均距離の計算結果を示す.

⟨ Φ(sym;K)| H − E | Ψ4 ⟩ = 0,αK(6)励起状態の位置, B4(1) = 127.33 mK, は trimer⟨ Φ(sym;H)| H − E | Ψ4 ⟩ = 0.
αH

(7)の束縛エネルギー B3(0) = 126.40 mK （前編のmax
max
ここで, αK = 1, ..., αKおよび αH = 1, ..., αH.

表 2）から見て僅かに 0.93 mK である. trimer,

式 (6),(7) は, 前編 (4)–(6) の形の一般化行列固

tetramer について, 図 1 の模式図で示したレベル

有値問題に行き着く.

構造が見られる. 確かに, §1.5 で触れた eﬀective(H)
基底関数 Φ(K)αK (ΦαH ) を次の様に角運動量表示

theory 計算による予言「ボソン 3 体の基底状態にする：

および各励起状態のエネルギーのすぐ下に, ボソ
ン 4 体の束縛状態が 2 個付随している」が, 4 He

(cos
sin )
n x lxΦ(K)(xi ) ψny ly (yi ) φnz lz (zi )αK (xi , yi , zi ) = ϕ[[]]
× Ylx (bxi )Yly (byi ) Λ Ylz (bzi ),trimer の基底状態については, realistic potential
計算でも成り立っていることが分かった.

JM

(8)

(i = 1, ..., 12)

空間分布は, 大雑把に言うと, tetramer の基底状態は, trimer の基底状態とほぼ同じサイズで

(cos
sin )
n x lxΦ(H)(xi ) ψny ly (yi ) φnz lz (zi )αH (xi , yi , zi ) = ϕ[[]]
× Ylx (bxi )Yly (byi ) Λ Ylz (bzi ).

ある. tetramer の励起状態においては, 3 粒子がほぼ trimer の基底状態にあり, 第 4 の粒子が他

JM(i = 13, ..., 18)(9)の粒子から遠くそれぞれ dimer の粒子間距離程

ここで αK は, 次の量を代表する（H -type につ

度離れている（図 2 参照）.

いても同様）：

前編で, dimerlike-pair model が予言した

αK ≡ {cos or sin, ω, nx lx , ny ly , nz lz , Λ, JM }.

tetramer の励起状態の束縛エネルギーは B4(1) =J, M は, 全角運動量とその z 成分. 本計算では,126.27 mK（LM2M2 potential) であったから,

J = 0 のみを扱うので（J > 0 には束縛状態は非常によく的中している. より一般的な予言式無い）, 4 粒子の対称化から, i) K -type に対し[前編の (34) 式, B2 は dimer の束縛エネルギー]て, lx = even, ly + lz = even および Λ = lz , ii)

(1)

(0)

B4
B
2
= 3 + ,
B2
B2
3

H -type に対して, lx = even, ly = even および(16)

c 2012 The Atomic Collision Society of Japan, All rights reserved.
Copyright⃝

35

0.008(v)

4

He tetramer の 4 体計算結果. B4(v =0, 1) は束縛エネルギー. rij は 2 体間の距
離. LM2M2 potential 使用. ⟨T ⟩, ⟨V ⟩ は, 運動エネルギー, ポテンシャル エネルギーの(1)
期待値. Ref. [14] の B4 は, 散乱 S 行列から外挿したもの.

Tetramer(0)

B4 (mK)⟨T ⟩ (mK)⟨V ⟩ (mK)⟨rij ⟩ (Å)

Tetramer(1)

B4 (mK)⟨T ⟩ (mK)⟨V ⟩ (mK)⟨rij ⟩ (Å)

基底状態

−3

(Å )

tetramer(v)

P4 (x)

表 1:

0.004(trimer v=1)

v=0

0.002

v=0

GEM [12]

Faddeev [14]558.984282.2557.7
4107−4841.2
8.43

−4665

v=1

(trimer v=0)

0
0

5

10

x (Å)

図 5:

8.40

励起状態

0.006

v=1

GEM [12]

Faddeev [14]127.331639.2

127.5−1766.5

4

He tetramer における 2 粒子相関関数(v)
P4 (x) [12]. trimer, dimer のケース（前編の図 4a) も描き入れてある. dimer 基底状態の線は trimer 励起状態にほぼ重なっている. 各曲線はピークで一致するよう normalizes してある. x < 4 Å でどの曲線も同じ形（同じ短距離相関）をしている → 前
編 (27) 式の 4 体系のケースの基礎付けにもなっている.

35.8

は, §2.5 の atomic Tjon lines の内の B3(0) -B4(1)相関において更なる威力を発揮する.

2.3

4

He tetramer の短距離相関

前編で述べた大きな心配 (1) 短距離相関の記
述と (2) 漸近形の記述は, 4 体系でも同じであり, さらに, 計算の困難さは増しているので, 心配はより強い. しかし, GEM の枠組みのよさと基底関数（ガウス関数のレンジ, 個数）の適切な採用 (文献 [12] の Table VI) により, これらの
図 6: 4 核子系（4 He 原子核）の基底状態における 2
粒子相関関数. 同状態に関して文献 [17] において行われた 7 グループよるベンチマー
クテスト計算の結果. 4 He tetramer における 2 粒子相関関数（図 5）の方が 4 He 原子
核の場合より遥かに強いことが分かる.

心配が見事に解決されている.
まず, tetramer での短距離相関を見るために,

trimer のとき, 前編 (29) 式, と同様に 2 体相関関数 P4(v) (x)(v = 0, 1) を定義する.
(v)

(v)

(v)

P4 (x1 ) Y00 (bx1 ) = ⟨ Ψ4 | Ψ4 ⟩y1 ,z1 .

(17)

図 5 に tetramer の基底 (v = 0)・励起 (v = 1) 状態

短距離相関）をしている—これには著者自身驚の結果を示してある [12]. trimer, dimer のケー

いた. このことは, 4 体計算において, 強い短距ス（前編の図 4a) も描き入れてある. dimer 基

離相関を直接扱うことを避けるために, Jastrow

底状態の線は trimer 励起状態にほぼ重なっていtype などの 2 体の correlation factor を 4 体系

る. 各曲線はピークで一致するよう normalizes

波動関数に a priori に掛けておく — という仮してある. x ≲ 4 Å でどの曲線も同じ形（同じ

定に基礎を与えることにもなっている [3 体系のc 2012 The Atomic Collision Society of Japan, All rights reserved.
Copyright⃝36

前編 (29) 式参照].

010

参考までに, 核力による短距離相関と比較し

−1/2

)

tetramer(z) (Å

てみる. 図 6 は, 4 核子系（4 He 原子核）の基底(v =0 ,v)状態における 2 粒子相関（図 5 と同じ物理量）の ∼ 10−2

z O4 3を示す [17]. この図の原点付近の値はピーク値程度であるが, 図 5 における He 原子
4

10

−2

10

−4

(dimer)v=0

tetramer での相関では ∼ 10−6 であり, 後者の短10距離相関の方が遥かに強いことを示している.

2.4

4(trimer v=1 )0

He tetramer 波動関数の漸近形

v=1

−6

500

1000

z (Å)

「求まった 4 体波動関数の十分遠方 (zi → ∞)が正しく e(v)−k4 ziziの形になっているか」を見るた
図 7:

めには，その形と overlap function
(v)

(0)

(v)

O4 (z1 ) Y00 (bz1 ) = ⟨ Ψ3 | Ψ4 ⟩x1 ,y1 ,(18)とを比べればよい．図 7 に z × O4(v) (z) を示す．
(v)

open circle は 遠方の正解 C4(v) e−k4 zi を示す（C4(v) の値は [12])．十分過ぎる遠方まで満足できる精度で求まっている. 但し, 図中の trimer,

4

He tetramer における overlap function(v)
O4 (z)(v = 0, 1) [12](z を掛けてある）の遠方での振舞い. trimer 励起状態と dimer のケースも描き入れてある. 丸印は, 結合エネ
ルギーから計算した漸近形. 遠方の 3 本の線がほぼ平行になっている＝ dimerlike-pairmodel の予言.

献 [2] の中の [32, 35, 47]).

dimer の線に比べると, さすがに 4 体計算の難4

しさが見えてしまうが.

He 原子の 3 体・4 体系には, それぞれ基底状

態と励起状態があるので, その束縛エネルギー

特筆すべきは, 遠方での 3 本の線がほぼ平行(B3(0) , B3(1) と B4(0) , B4(1) ) の間には, 計 6 通りの相になっていることである（k4(1) ≃ k3(1) , k4(1) ≃ k2 ).

関がある（generalized Tjon lines と呼ばれる）.

これは, 前編の trimer 計算で提唱した dimerlike-

文献 [13] で, これらの相関図 (universal scaling

pair model が, tetramer でも成り立つことを保curve と呼ばれる）が leading-order eﬀetive the-

証している.

ory を用いて計算されて, 4 He クラスターの妥当

2.5

4

な束縛エネルギーの範囲で, どれも直線になってHe の generalized Tjon lines「S 波散乱長が大きい短距離相互作用をする

いることが示されている（図 8 と図 9 の中の黒実

3 体・4 体系」が見せる universality の内, 束縛

線として転載した）. しかし, 当時, 4 He tetramer

エネルギーの間の相関は好んで議論される. 原(特に励起状態）の信頼できる計算がそろってい子核物理では, 3 He 原子核, 4 He 原子核の基底状なかったので, 相関の直線性や non-universality

態の束縛エネルギー（B3(0) , B4(0) ) 間の相関が著

eﬀect の有無については議論が進まなかった.

名であり, 図 3 のように, 様々な 2 体の realistic

筆者等は, これらを解明すべく, 文献 [2] にお

potential で計算した値が, ほぼ直線（Tjon line,

いて, 発表されている realistic potential を 14 通

傾斜=4.8) に乗ることが知られている（そうなり持ち出し（本稿ではその中身に触れないが, 文

る理由は, 核物理側では深くは議論されていな

献 [2] に詳述されている）, その結果を, 6 通り

かった）.

の相関図に plot した. 図 8 と図 9 の中の closed

3 体・4 体ボソンの leading-order eﬀective the-circles と open squares がそれを表す. 赤点線はory の計算では, 散乱長が無限大である場合 (uni-

最小 2 乗法による直線での ﬁt である. どの相(0)(0)
B4 /B3 ≈ 4.6,

関もきれいに直線に乗っているので驚いた（異≈ 1.01 であることが知られていた（文なる potential を使った計算であり, ある特定のtary limit と呼ばれる）には,
(1)

(0)

B4 /B3potential のパラメタを変化させた計算ではない
c 2012 The Atomic Collision Society of Japan, All rights reserved.
Copyright⃝

37

100

100
(1)

(1)

(0)

a) B3 −B4

90

(1)

(1)

B4 / B290

B4 / B2(1)b) B3 −B480

Dimerlike−pairmodel, Eq.(16)80
Linear fit

Linear fit

Ref. [13]

70
1.55

1.60

1.65

1.70

Ref. [13]

70
70

1.75

80

90

100

440

c)

(1)
(0)
B3 −B4(0)(0)d) B3 −B4
B4 / B2

400

400
Linear fit(0)Linear fit(0)

B4 / B2100

B3 / B2

440

Ref. [13]360

Ref. [13]

360

320

320
1.55

1.60

1.65

1.70

70

1.75

80
(0)
B3 / B2(1)
B3 / B2

図 8:

90
(0)

(1)

B3 / B2

(1)

4

(1)

He trimer と tetramer の基底状態・励起状態のエネルギー相関（generalized Tjon line）[2]. a) B3 -B4
(0)
(1)
(1)
(0)
(0)
(0)
相関, b) B3 -B4 相関, c) B3 -B4 相関 および d) B3 -B4 相関を表す. 各エネルギーは dimer のエネルギー B2 で normalize してある. 14 個の closed circles, open squares は, 各種の 4 He-4 He realistic
potential を用いた計算結果を表す. potential の詳細については文献 Ref. [2] 参照. 赤点線は, 計算結果
を最小 2 乗法で直線近似したもの. b) の中の青鎖線は, dimelike-pair model による予言, (16) 式. 黒実
線は, 4 He trimer, tetramer に対する leading-order eﬀective theory による計算 [13].

100

1.75
(0)

(0)

(1)

(1)

b) B4 −B4a) B3 −B3

90

(1)

B4 / B2(1)

B3 / B2

1.70

1.65

80

Linear fit

1.60

Linear fit

Ref. [13]Ref. [13]

1.55
70

80

90

70

100

(0)
B3 / B2

図 9:

320

360

400

440

(0)
B4 / B2

(1)

4

(1)

He trimer と tetramer の基底状態・励起状態のエネルギー相関（generalized Tjon line）[2]. a) B3 -B3(0)(1)
相関, b) B4 -B4 相関を表す. その他は, 図 8 と同じ.

ことは特筆すべきである）. 最小 2 乗法による2-3%である）. しかし, 図 8(c),(d) および図 9(b)直線の係数, 計算値と直線のづれの具合などは,においては, 無視できないずれが見える. 後者

文献 [2] に記載されている.

のずれは, eﬀectcive theory が, B4(0) (tetramer

黒実線で示された leading-order eﬀective the-

基底状態のエネルギー) の値を 12-14% under-

ory 計算 [13] は, 図 8(a),(b) および図 9(a) にお

estimate していることが原因である. それがいては計算値とかなりよく一致している（ずれはnon-universal eﬀect であり, next-leading-orderc 2012 The Atomic Collision Society of Japan, All rights reserved.
Copyright⃝

38

eﬀect を持ち出すか, 何らかの 4 体パラメタを導本研究は, 筆者等が自身の計算法を携えて, 原

入する必要があろう. 今後, これらの研究が進子核分野から冷却原子分野に進出して行った例

むものと思われる.

であるが, 一方で, 「逆向き」に, 冷却原子分野

図 8(b) の中の青鎖線は, dimerlike-pair modelの主力理論研究者が, 同分野の基幹の論理を展の予言式 (16) を plot したものである. 計算値

開・適用すべく, 原子核の弱結合系（クラスター

および赤点線とほぼ重なっており, 予言の正確

物理, 不安定核物理など）に進出し始めているさを表している.

例（その種の研究のレビュー論文は文献 [18]）を横目で見て, 複雑な気持ちで居る. 確かに, 広い

3. おわりに概念 (universal few-body physics) で分析・整理

筆者等が, He 原子の trimer, tetramer の理が行われようとしていることは, 納得できると論研究に参入した動機は, 前編の §1 に述べたよ

同時に, 内心脅威を感じている. 両分野の交流

うに, 原子核で培った少数粒子系計算法「ガウが深く行われ, 両者の発展に寄与して行くことス関数展開法, GEM」を適用するためであった.

を祈って筆を措く.

4しかし, 研究を進めて行くうちに, 原子核の弱結合系の物理と 4 He 原子クラスターの物理（広参考文献

く冷却原子気体の物理）との間の共通性に関し[1] E. Hiyama, B.F. Gibson and M.
Kamimura, Phys. Rev. C 70, 031001(R)(2004).
[2] E. Hiyama and M. Kamimura, Phys. Rev.
A 85, 062505 (2012).
[3] J.A. Tjon, Phys. Lett. B 56, 217 (1975).
[4] A. Nogga, S.K. Bogner and A. Schwenk,
Phys. Rev. C 70, 061002 (2004).
[5] E. Braaten and H-W. Hammer, Phys. Reports 428, 259 (2006).
[6] V. Eﬁmov, Yad. Fiz. 12, 1080 (1970)[Sov. J. Nucl. Phys. 12, 589 (1971)].
[7] 数納広哉,「しょうとつ」第 8 巻第 6 号,
[原子衝突のキーワード—エフィモフ状態],
p.36 (2011).
[8] J. von Stecher, J.P. D’Incao and C.H.
Greene, Nature Phys. 5, 417 (2009).
[9] F. Ferlaino et al., Phys. Rev. Lett. 102,
140401 (2009).
[10] P. Naidon, E. Hiyama and M. Ueda,
Phys. Rev. A 86, 012502 (2012).
[11] R.A. Aziz and M.J. Slaman, J. Chem.
Phys. 94, 8047 (1991).
[12] E. Hiyama and M. Kamimura, Phys. Rev.
A 85, 022502 (2012).
[13] L. Platter, H.-W. Hammer and UlfG. Meissner, Phys. Rev. A 70, 052101(2004).
[14] R. Lazauskas and J. Carbonell, Phys.
Rev. A 73, 062717 (2006).

て, 自らの計算結果を通して, 認識が深くなってきた. 両分野に通底しているものは, universal

few-body physics (または, universality in fewbody systems with large scattering length) という言葉で括れる.
原子核の世界では, 2 体相互作用の大きさを外から変化させることはもちろん不可能だが, しかし, 冷却原子気体の世界では, 外から磁場をかけ, Feshbach resonance を利用して, 2 体短距離相互作用の散乱長を微調整できる（つまり, 相互作用の大きさを変化させられる）ことは驚異である. µK, 10−10 eV のレベル世界のことだから,
原子核の研究とは関係ない, と見過ごすことはできない. なぜなら, universality を通して, 同様なことが原子核の世界でも見られるはずだか
らである (例を §1 の図 1-3 で示した）.
冷却原子気体の世界で, 少数粒子系の構造や反応メカニズムを, 相互作用を変化させて調べることは, 原子核の弱結合系の世界で仮想的に相互作用を変化させて起こるダイナミクスを,

universality を通して, simmulate している可能性がある. そうであれば, これは原子核の研究に大いに役立つかも知れない — と想像（空想?）している.

c 2012 The Atomic Collision Society of Japan, All rights reserved.
Copyright⃝39

[15] M. Przybytek et al., Phys. Rev. Lett.
104, 183003 (2010).
[16] E. Hiyama, Y. Kino and M. Kamimura,
Prog. Part. Nucl. Phys. 51, 223 (2003).
[17] H. Kamada et al., Phys. Rev. C 64,
044001 (2001).
[18] H.-W. Hammer and L. Platter, Ann. Rev
Part. Sci. 60, 207 (2010).

Few-body systems interacting withstrong short-range repulsion – Eﬁmov physics and 4 He-atom trimer andtetramer
Masayasu KAMIMURA and Emiko HIYAMA
Theoretical study of few-body systems interacting with strong short-range repulsion is reported. Especially discussed isthe energy
level, spatial structure, correlation between thebinding energies and short-range correlationsof the 4 He trimer and tetramer that are calculated with the use of the Gaussian expansion
method (GEM) developed by the authors. Relation to Eﬁmov physics is also discussed.

c 2012 The Atomic Collision Society of Japan, All rights reserved.
Copyright⃝

40

「原子衝突のキーワード」の他，蓄積リングを用いた合流ビーム法や，電子ビームイオントラップ（Electron Beam Ion

2 電子性再結合(Dielectronic Recombination)

Trap: EBIT）を用いる方法などがあるが，重元
素多価イオンの場合には，後出の 2 つが主な方法となる [3]．

原子あるいは分子から電子を剥ぎ取ることを

重元素多価イオンの DR の研究は，プラズマ「電離」と言う．電離で生成される正の電荷を素過程として重要であるばかりでなく，「原子持ったイオンが，再び電子を奪い返すことを「再衝突学」的にも興味深い話題を多く含んでいる．

結合」(recombination) と言う．つまり再結合は

反応式 (1) および (2) から分かるように，RR と電離の逆過程と言える．電離にはエネルギーのDR は始状態と終状態が同じであるため，原理的付与が必要であるのに対して，再結合ではエネルにそれぞれ干渉し合い，ヘリウムの光電離スペ

ギーが放出される．このエネルギーが分子の解クトルで得られるような非対称な共鳴形状（い離に供される過程が，第 9 巻第 3 号で紹介された

わゆる Fano プロファイル）を示す．一般に多「解離性再結合」である．一方，エネルギーを光子

価イオンでは DR が RR に比べ桁違いに大きいとして放出する過程は「放射性再結合」(radiative

ため干渉効果は小さいが，その定性的な理解に

recombination, 以下 RR) と呼ばれる．

反して明瞭な非対称形状が観測されることもあ

→A

(1)

る．相対論効果が強い重元素多価イオンにおけここで h はプランク定数，ν は放出される光のる干渉効果を計算するのは容易ではないが，理振動数である．上式からも明らかなように放射論的にも研究が進められている [4]．また，とき性再結合は光電離の逆過程である．光電離におに相対論効果が支配的に現れることがある．電

いては，光子が直接電子をはじき出す直接電離通大の Tokyo-EBIT を用いた Au や Bi など重元過程と，2 電子励起状態を介した間接電離過程

素の Li 様イオンの DR の測定では，ブライト相があるが，後者の逆過程が本稿の主題である「2互作用（電子間相互作用における相対論効果）が電子性再結合」(dielectronic recombination, 以

共鳴強度を倍増させたり，放射 X 線の角度分布

下 DR) である．

を支配したりすることが示された [5]．「ブライRR : e + A

q+ DC

DR :e+Aq+−−→ A(q−1)+(q−1)+∗∗→A+ hν(q−1)+

ト相互作用」については，続号のこのコーナー+hν (2)

DR は上式に記したように，自動電離の逆過程

で改めて紹介される予定である．

である 2 電子性電子捕獲 (dielectronic capture:(電気通信大学 中村信行)

DC) と，それに伴い生成される 2 電子励起状態

参考文献の光放射脱励起から成る．

DR は特定の衝突エネルギーにおいて極めて[1] H. S. W. Massey and D.R.Bates, Rept.

大きな断面積を持つ共鳴過程であり，プラズマ

Prog. Phys. 9, 62 (1943).

中イオンの電離平衡状態での価数分布に強く影[2] A. Burgess, ApJ 139, 776 (1964).

響する．その重要性は 1940 年代に Massey と[3] 渡辺裕文，加藤太治，プラズマ核融合学会

Bates[1] によって指摘された．1960 年代には誌 83 巻 660 (2007).

Burgess[2] が，10 K 程度の高温プラズマであ
6[4] X. M. Tong et al., Phys. Rev. A 80,

る太陽コロナにおいて，鉄多価イオンの DR が042502 (2009).

重要であることを示した．以降，天体プラズマ[5] N. Nakamura et al., Phys. Rev. Lett. 100,や核融合プラズマに関係する多価イオンの DR

073203 (2008); Z. Hu et al., Phys. Rev.

の研究が精力的に進められてきた．

Lett. 108, 073002 (2012).

DR を実験的に研究するには，交差ビーム法
c 2012 The Atomic Collision Society of Japan, All rights reserved.
Copyright⃝

41

「原子衝突のキーワード」(2) の値を中心に分布し，その幅は標的内電子の運動量分布できまる．断面積や光子の角分布・放射性電子再結合と放射性電子捕獲(Radiative Recombination and Radiative
Electron Capture)

偏光も RR と同様の特徴を持つ．ただし実験室
系ではほぼ静止した電子に高速のイオンが衝突するので，光子の放出角度・エネルギーや微分
断面積はイオン静止系からローレンツ変換する

放射性電子再結合 (RR) は，イオンが自由電必要がある．なお 10 MeV/核子以上のイオンと子を捕獲すると同時に光子を放出する過程
e + Aq+ → A(q−1)+ + ℏω

原子分子の衝突による電子捕獲では REC の寄与が重要となる．

(1)

で，連続状態から離散状態への 1 電子遷移過程

REC の実験は 1970 年代から重イオン加速器

すなわち光電効果の逆過程である．プラズマでを用いて数多く行われている．一方 RR の実験

重要な原子過程で，また重イオン蓄積リングのには自由電子とイオンの衝突装置が必要で，重電子冷却によるビーム損失にも関与する．

イオン蓄積リングと電子ビームイオントラップ

静止したイオンに対して運動エネルギー T で(EBIT) の開発にともない 1990 年代以降に研究入射した電子が RR により束縛エネルギー EB

が大きく進展した．イオン蓄積リングでは周回の状態に捕獲される場合，光子のエネルギーはするイオンビームが同じ速度の電子ビームと合ℏω = T + EB

(2)

流して起こるきわめて低い相対速度の RR が，となる．終状態がネオンより重いイオンの 1s 軌

また EBIT では蓄積された高電離イオンと高速道であれば，光子エネルギーは 1keV 以上の Xの電子ビームの衝突による RR が観測されてい線の領域である．

る．REC と RR を合わせると，実験に用いられ

速度 v の電子が原子番号 Z の裸核の 1s 軌道にたイオンは水素からウランまで，またイオンにRR で捕獲される全断面積は，非相対論（Zα ≪ 1，

対する電子の速度はゼロから 0.99998c にわたり，v ≪ c，ただし α は微細構造定数）と電気双極子

全断面積や角分布，終状態分布だけでなく最近遷移を仮定すれば，ν = Ze /ℏv として

では X 線の偏光の測定も行われている．実験結2

σRR = 9.2×10−21(3ν
1 + ν2)2果によると，電気双極子近似による全断面積の−4ν arctan(1/ν)

e

1 − e−2πν(3)式（3）は広範囲の衝突系と衝突速度でそこそことなる（単位 cm2 ）[1]．σRR は v が同じならば成り立つが，Zα ∼ 1 あるいは v ∼ c では相対論

Z とともに増加し，Z が同じならば大きな v で的効果・多重極遷移の寄与が顕著である．また−5に従って減少する。また終状態の主量子数

終状態が励起状態の場合の磁気量子数分布，他n に対して近似的に Z R /nT (n T + Z R)（たの電子捕獲過程との干渉などについても，理論

だし R はリュードベリ定数）に比例する．終状的な研究と相まって多彩な研究が行われている．

態が 1s の場合、電子入射と光子放出の間の角度

くわしくは参考文献 [2] を参照．

v

4

2

2

2

θ に対して，光子の角分布は sin2 θ に比例し，電(理化学研究所 神原 正)

子光子の面内に 100%直線偏光する．

参考文献

放射性電子捕獲（REC）は，原子・分子と衝突したイオンが標的から電子を捕獲して光子を[1] H. A. Bethe and E. E. Salpeter,

放出する過程q+

Quan-tum Mechanics of One- and Two-Electron+B→A(q−1)++ B + ℏω+

(4)

Atoms , (Plenum, New York) 1977, p320.

である．標的にゆるく束縛された電子を自由電[2] J. Eichler and Th. Stöhlker, Phys. Rep.

A

子とみなせば，REC は運動量分布を持つ自由電

439, (2007) 1-99.

子の RR で近似できる．光子のエネルギーは式
c 2012 The Atomic Collision Society of Japan, All rights reserved.
Copyright⃝

42

2012 年度 役員・委員会等会長
髙橋正彦（東北大学）
幹事
渡部直樹（北海道大学）（副会長）森下亨（電気通信大学）

足立純一（高エネルギー加速器研究機構） 星野正光（上智大学）
運営委員
石井邦和（奈良女子大学）

高口博志（広島大学）

星野正光（上智大学）

間嶋拓也（京都大学）

美齊津文典（東北大学）

本橋健次（東洋大学）

森下 亨（電気通信大学）

渡辺信一（電気通信大学）

足立純一（高エネルギー加速器研究機構） 岸本直樹（東北大学）
小島隆夫（理化学研究所）

冨田成夫（筑波大学）

日高 宏（北海道大学）

渡部直樹（北海道大学）

渡辺 昇（東北大学）
会計監事
城丸春夫（首都大学東京）
中村義春
常置委員会等
編集委員会

委員長：

渡部直樹（北海道大学）

行事委員会

委員長：

森下

広報渉外委員会委員長：

足立純一（高エネルギー加速器研究機構）

若手奨励賞選考委員会委員長：

大野公一（豊田理化学研究所）

国際会議発表奨励者選考委員会

委員長：

髙橋正彦（東北大学）

学会事務局

担当幹事：星野正光（上智大学）

亨（電気通信大学）

編集委員会
足立純一，岸本直樹，長嶋泰之，中井陽一，羽馬哲也，早川滋雄，日高宏

森林健悟，渡部直樹

しょうとつ

第９巻 第５号（通巻 48 号）

Journal of Atomic Collision Researchc 原子衝突学会 2012⃝
http://www.atomiccollision.jp/
発行: 2012 年 9 月 15 日
配信: 原子衝突学会 事務局<acr-post@bunken.co.jp>