筑波大学 計算科学研究センター 平成 26 年度 年次報告書

Ⅱ-2. 原子核物理分野
1. メンバー
教授

矢花一浩、中務 孝

准教授

寺崎

講師

橋本幸男

助教

日野原伸生（国際テニュアトラック）

研究員

温

学生

大学院生６名、学類生２名

順（HPCI 戦略プログラム）

凱（数理物質系物理学域）

2. 概要
核子（陽子・中性子）の多体系である原子核の構造・反応・応答などの多核子量子ダイナミ
クスの研究を推進している。安定線（ハイゼンベルグの谷）から離れた放射性アイソトープの原子核の構造と反応、エキゾチックな励起状態の性質、様々な集団運動の発現機構など、未解決の謎の解明に取り組んでいる。原子核の研究は、フェルミ粒子の量子多体系計算という観点で、物質科学や光科学、冷却原子系の物理と密接なつながりをもつ。また、クォーク・グルーオンのダイナミクスを記述する格子 QCD に基づく核力の計算、軽い原子核の直接計算などが進展する中、素粒子物理学との連携も重要性が増している。ニュートリノの解明に向けたニュートリノレス二重ベータ崩壊の観測実験や、素粒子標準模型のテストに関わる実験などにも原子核理論の精密計算が不可欠とされている。また、元素の起源や星の構造にも原子核の性質は深く関わり、宇宙物理学とも密接に関係している。本部門・原子核物理分野のメンバーは、このような幅広い課題に取り組み、分野の枠を超えた研究を推進している。

3. 研究成果【 1 】時間依存ハートレー・フォック法による核子移行過程の記述（矢花、関澤）
多核子移行反応は、核子多体系に発現する非平衡量子輸送過程としてその反応機構に興味が持たれており、また近年、多核子移行反応を用いこれまでに生成されたことのない不安定
原子核を生成する可能性が盛んに議論されている。我々は、多核子移行反応の微視的反応機構を明らかにすると共に、不安定核の生成方法を予言することを目指し研究を進めてきた。
昨年度までに、TDHF 法に粒子数射影法を組み合わせることにより、起こる確率の大きい移行
反応過程の断面積を定量的に記述できることを明らかにした。しかし、移行する核子数が大きく、起こる確率の小さい移行反応過程については、実験値を過小評価するという問題に直面した。そこで本年度には、この実験値と食い違いを解消することを試みた。相違の原因として、励起した原子核からの核子放出による脱励起（核子蒸発）の効果が考えられる。核子- 63 -

筑波大学 計算科学研究センター 平成 26 年度 年次報告書

蒸発過程は、我々の数値シミュレーションで扱う時間スケールより遥かに長い時間スケールで起こる現象であり、TDHF 計算とは別に統計模型を用いた計算を行う必要がある。統計模型を利用するために、我々は粒子数射影法を応用し、生成された原子核の励起エネルギーを評価する方法を開発した。そして、核子蒸発の効果を取り入れた断面積を評価し、実験値との比較を行った。その結果、核子蒸発の効果が実験値との相違を改善することを確認したが、効果は不十分であった。このことは、TDHF 計算で扱えていない相関の影響を示唆している。

【 2 】正準基底時間依存ハートレー・フォック・ボゴリューボフ理論による多核子ダイナ
ミクス計算（中務、江幡（北大）、稲倉（京大））
原子核理論分野で発展が目覚ましい 3 次元実空間表示の時間依存ハートレーフォック法に対して、計算コストの増加をミニマムに抑えた上で、核子対凝縮と超流動ダイナミクスを扱えるように拡張した理論を、2010 年に我々が提唱した。このときの論文に対して、第一著者である江幡氏（当時筑波大学大学院生・理化学研究所 JRA）が、本年度の日本物理学会若手奨励賞を受賞した。この方法を線形応答領域に応用し、原子核の応答関数の系統的計算を行い、軽い原子核から希土類領域の重い原子核まで、中性子分離エネルギーが 2MeV 以上のアイ
ソトープについて電気双極子(E1)強度分布を調査した。その結果、中性子過剰核の広い領域
において低エネルギーの E1 ピーク（ピグミー共鳴）が現れ、その出現のメカニズムに中性子の殻効果が重要であることが分かった。また、通常、E1 強度のほとんどは巨大共鳴状態に集中するため、低エネルギー領域の E1 強度は大幅に弱められていることが知られているが、非常に中性子過剰な重い原子核の励起エネルギー5 MeV 以下には、この E1 抑制が働かない新種のピグミー共鳴が現れることを予言した。さらに、これらのピグミー共鳴の性質と中性子星の構造を決定する状態方程式との関係を精査した。
また、同じ計算手法を用いた核融合反応計算を開始し、軽い原子核に限定的ではあるが、対相関の影響に関して興味深い結果を得た。

図１：Zr アイソトープの低エネルギーE1 強度の割合。

衝突エネルギーを同じにして、対特定の中性子数 N=50, 60, 72, 82 で傾きが突然変化す

相関を入れた計算と無視した計算- 64 -

筑波大学 計算科学研究センター 平成 26 年度 年次報告書とを比べると、対相関によって核融合が阻害されるケースがあることが分かった。さらに
この効果を確認するため、研究を続行中である。

【 3 】大振幅集団運動理論を用いた核反応ダイナミクスの記述（温、中務）
線形領域を超える大振幅集団運動を扱う理論として、断熱近似型の理論が提案されているが、そのうちの一つが、2000 年に提案された断熱自己無撞着集団座標法(Adiabatic Selfconsistent Collective Coordinate Method: ASCC 法)である。この理論では、少数自由度の集団空間（座標）の自己無撞着な抽出が可能であり、特に、過去の断熱時間依存平均場理論では不可能であった一意的な抽出ができる点が優れている。この理論に基づいて、低エネ
ルギーの多核子反応ダイナミクスを記述することを目指した研究を開始した。まず、テストとして、比較的簡単な BKN 相互作用と呼ばれる有効相互作用を用いて、3 次元座標表示のTDHF を基礎にしたコード開発を行い、2 つのアルファ粒子の散乱・融合のミクロな機構を明らかにすることにした。拘束演算子を局所調和方程式と呼ばれる式から求め、それを拘束条
件付き平均場計算に応用する。これを自己無撞着に計算し、核反応の集団運動を支配するポテンシャルと質量パラメータを求める。これまでに、このようにして自己無撞着に決めた質量が遠方で換算質量に一致すること、また 2 つのアルファ粒子が接触した後には換算質量が大きくなり、反応速度が弱められることなどが分かっている。

【 4 】アイソスピン不変なエネルギー汎関数による原子核密度汎関数計算（中務、佐藤（理研）、Dobaczewski（ワルシャワ大）、 Satula (ワルシャワ大））
現在主流となっている原子核のエネルギー密度汎関数は、Skyrme 形式、Gogny 形式、共変()

形式（相対的）の 3 つに大別されるが、どれも陽子と中性子の密度 ρ p , ρ n の汎関数としてエネルギーが与えられている。しかし、陽子や中性子はアイソスピンの第 3 成分の固有状態であり、アイソスピン空間における回転に対して不変ではなく、一般にはアイソスピンが任意の方向を向いた状態、すなわち陽子と中性子が混合した状態に拡張する必要がある。これを実行するため、陽子・中性子を区別せずに「核子」として扱う新しい Kohn-Sham 方程式()と、それに対応する非対角要素 ρ pn , ρ np を含むエネルギー汎関数を構築し、その計算コード開発を実施した。アイソスピン不変に拡張することによる新しいパラメータの導入はな
く、これまでのエネルギー汎関数のパラメータセットをそのまま利用することができる。計
算コードは、ポーランド・ワルシャワ大学のグループを中心に開発されている HFODD にこの拡張を適用することにした。この研究はワルシャワ大学との共同研究で進められている。

- 65 -

筑波大学 計算科学研究センター 平成 26 年度 年次報告書

この拡張による大きなメリット
として、荷電類似状態（アイソバリ
ックアナログ状態）と呼ばれる状態
が、自然かつ簡便な形で表現できることがあげられる。図には、質量数
A=54 のアイソスピン T=1 の状態が示されている。左端の点(T=−1)が54

Ni に、右端(T=+1)が 54Fe、真ん中(T=0)が 54Co に対応する。この真ん中の 54Co の状態は従来の密度汎関数で記述することが困難であった。さらに、パラメータフリーの計算で実 図２：A=54 の同重体におけるアイソスピン T=1 三重項状
験と非常によい一致を得ている。

態のエネルギー（左から Ni, Co, Fe に対応）。

【 5 】原子核の二重ベータ崩壊の原子核行列要素（寺崎）寺崎は四年前に赴任して以来ニュートリノレス二重ベータ崩壊の理論的研究を実施している。その中心的内容は、準粒子乱雑位相近似を用いたその崩壊の原子核行列要素の計算である。その崩壊が観測された場合にニュートリノの質量スケールを決定するために必要な原子
核行列要素の信頼できる数値を出すことが目標である。2014 年度は、150Nd から 150Sm へのニュートリノレス二重ベータ崩壊の原子核行列要素を準粒子乱雑位相近似（QRPA）を用いて求め、その計算に関する論文を出版した。その計算過程で、始状態と終状態から求めた QRPA 励
起状態の重複を求める必要があるが、ここで QRPA 基底状態の規格化因子の積が約 2 であり、原子核行列要素を 1/2 倍する効果があることを初めて見出した。このことは、今まで原子核
行列要素を小さくするために用いられていた現象論的相互作用の方法の見直しの必要性を示しており、同じ課題を研究している研究者への影響は大きい。今までより理論的整合性の高い計算を提示したことで、信頼性の高い原子核行列要素を求めるという目標に近づいた。

【 6 】有限振幅法を用いた和則計算 (日野原、Kortelainen (Univ. Jyväskylä), Nazarewicz(Michigan State Univ.), Olsen (Michigan State Univ.))
和則は巨大共鳴などの線形応答に対する情報と関連する扱いやすい量であるため、原子核
エネルギー密度汎関数の未知の結合定数を決定する際に、ダイナミクスを代表する量として有用であると考えられる。しかしながら乱雑位相近似(QRPA)を用いて和則を計算すると、大- 66 -

筑波大学 計算科学研究センター 平成 26 年度 年次報告書

次元の行列の対角化とそのすべての解が必要となり、さらに計算量による制限のため行列の次元を減らす近似が必須である。全自由度を扱うことが出来る有限振幅法を用い、複素エネ
ルギー平面で励起状態に対応する極をすべて囲む複素積分を行うことによって、和則を効率よく計算する定式化を行った。全自由度が扱える小さな模型空間において QRPA の行列対角化による解と一致することを数値的に確認し、さらに和則を与える Thouless 定理や絶縁定理の密度汎関数理論での有効性を、現実的な場合において数値的に示した。有限振幅法による本定式化はこれらの定理が使えない場合においても有効であり、複素積分の並列化によって和則を効率よく正確に計算することが可能となった。

【 7 】二重ベータ崩壊の核行列要素（陽子−中性子対振幅の量子ゆらぎを取り込んだ核行列
要素の計算 (日野原、Engel(Univ. North Carolina))
二重ベータ崩壊の核行列要素はアイソスカラー型の陽子−中性子対相関力によって抑制されることが乱雑位相近似(QRPA)計算によって知られているが、この陽子−中性子対相関力は基底状態に質的な影響を与えないため、その強さについてはよくわかっていない。QRPA 計
算は相転移を記述出来ないため、基底状態の四重極変形ゆらぎが大きく変形共存や変形転移状態にある場合や、アイソスカラー対凝縮相への転移状態近傍にある場合は近似が破綻する。この QRPA の問題を回避し核行列要素の相互作用依存性を明らかにするために、四重極変形度と陽子−中性子対振幅を拘束した変分によってこれらの量子ゆらぎの大きい状態をも含んだ基底を生成し、量子数射影後に重ねあわせることで、崩壊の始状態と終状態を構築した(生成座標法)。特に陽子−中性子対を自由度に取り込んだ生成座標法はこれまでに行われたことのないユニークなものである。第一の計算例として、二主殻模型空間で 76Ge→76Se の二重ベータ崩壊の核行列要素を評価した。生成座標法では殻模型計算より広い一粒子模型空間を用いることが出来、さらに二準粒子励起に制限される QRPA よりも豊富な多体相関を取り込むことが出来るため、適切な有効相互作用を用いることで、将来的には殻模型計算と
QRPA のどちらにも含まれていない効果を議論することが期待できる。もう一つの将来計画
は Skyrme 型の原子核密度汎関数に基づいて陽子−中性子対の自由度を入れた生成座標法計算
を行うことで、β崩壊などの関連する物理量も含めた荷電交換過程の系統的な計算へと発展させることである。

【 8 】Gogny-TDHFB による２０Ｏ-２０Ｏ正面衝突の計算 （橋本）
原子核の振る舞いを研究するうえで、時間依存平均場の方法は定性的にも定量的にも極めて有用な枠組みとなっている。原子核の基底状態近傍の状態においては対相関が重要な働きをすることが知られているので、中性子過剰核の性質を理解する目的で、TDHF の枠組みを- 67 -

筑波大学 計算科学研究センター 平成 26 年度 年次報告書

拡張して対相関を扱えるようにした時間依存ハートレーフォックボゴリュボフ（TDHFB）の実用的な数値計算が行われるようになってきた。橋本は、調和振動子基底と空間格子としての Lagrange 格子点を組み合わせた基底を利用し、Gogny 力を用いた TDHFB を解く、という数値計算の方法を整備してきた。今年度は、２０Ｏ+２０Ｏの衝突に応用した。ちなみに、
TDHF+BCS の形式での原子核衝突の計算は存在するが、TDHFB 方程式を直接解いて原子核衝突を計算する試みは、現時点では極めて少ない。計算に必要な格子空間の選定は初期条件の設定と関係し、そのための試行錯誤に計算時間と労力を割く必要があった。実際の計算は、クーロン障壁頂上付近の軌道を３本（３個の初期エネルギー）設定し、２１fm 離れた点から計算を開始した。この一連の計算により、① 障壁より上のエネルギーの軌道では融合が起
こり、核子の移行量の期待値を計算できた。② 障壁より下の軌道において、対相関の働く中性子側では、陽子側の約１０倍の移行数があり、これは２０Ｏの中性子数・陽子数比である
１．５よりずっと多いことがわかる。さらに、③ 障壁よりも上のエネルギーを持つ２本の軌道から、巨視的な摩擦係数を計算することができた。摩擦係数は、融合の初期に急速に増加することが示された。今後、これらの量の初期エネルギー依存性を系統的に調べていく。
【 9 】Energy dependence of the nucleus-nucleus potential and the friction parameter infusion reactions （温、坂田(茨城大)、Li(CIAE)、Wu(CIAE)、Zhang(CIAE)、Zhou(CAS))
Applying a macroscopic reduction procedure to the improved quantum molecular
dynamics (ImQMD) model, the energy dependences of the nucleus-nucleus potential, thefriction parameter, and the random force characterizing a one-dimensional Langevin-type
description of the heavy-ion fusion process are investigated. Systematic calculations with
the ImQMDmodel show that the fluctuation-dissipation relation found in symmetric headon fusion reactions at energies just above the Coulomb barrier fades out when the incident
energy increases. It turns out that this dynamical change with increasing incident energy
is caused by a specific behavior of the friction parameter which directly depends on the
microscopic dynamical process, i.e., on how the collective energy of the relative motion istransferred into the intrinsic excitation energy. It is shown microscopically that the energy
dissipation in the fusion process is governed by two mechanisms: One is caused by the
nucleon exchanges between two fusing nuclei, and the other is due to a rearrangement of
nucleons in the intrinsic system. The former mechanism monotonically increases thedissipative energy and shows a weak dependence on the incident energy, while the latter
depends on both the relative distance between two fusing nuclei and the incident energy.
It is shown that the latter mechanism is responsible for the energy dependence of the- 68 -

筑波大学 計算科学研究センター 平成 26 年度 年次報告書

fusion potential and explains the fading out of the fluctuation-dissipation relation.

4. 教育

1．五嶌祐希、修士（理学）、“Gogny 力を用いた HFB による質量数Ａ≦５０領域における中性子過剰原子核の構造”
2 ． 関 澤 一 之 、 博 士 （ 理 学 ）、“ Multinucleon transfer reactions and quasifission
processes in time-dependent Hartree-Fock theory”

5. 受賞、外部資金、知的財産権等受賞
1. 日本物理学会若手奨励賞、江幡修一郎、“Canonical-basis time-dependent
Hartree-Fock-Bogoliubov theory and linear response calculations”、2015 年 3 月
21 日
2. 筑波大学数理物質科学研究科長賞、関澤一之、“Multinucleon transfer reactions
and quasifission processes in time-dependent Hartree-Fock theory”、2015 年 3月 25 日3. The Award for Best Project of the 2014 TALENT Course #5 “Theory for
Exploring Nuclear Structure Experiments”、関澤一之、2015 年 12 月
4. 「京」を中核とする HPCI システム利用課題 優秀成果賞, 橋本幸男、関澤一之、
矢花一浩、2014 年 10 月

外部資金
1. 日本学術振興会科学研究費・基盤研究(B)、中務孝、代表、2013 年より継続、
5,500,000 円、「原子核の低エネルギー集団励起と核融合・核分裂機構の解明」
2. 科研費・新学術領域研究（研究領域提案型）、中務孝、分担、2012 年より継続、
1,000,000 円、「冷却原子を用いた中性子過剰な低密度核物質の状態方程式」
3. 学振・二国間交流事業（オープンパートナー）、中務孝、代表、2014 年、2,484,000 円、セミナー「計算核物理学の進展と展望」
4. JST ImPACT「核変換による高レベル放射性廃棄物の大幅な低減・資源化」、中務孝（研究開発責任者）、矢花一浩、橋本幸男（研究開発参加者）、2014 年、3,000,000 円、「核構造計算による核反応モデルの高精度化」
5. 科研費 平成 26 年度基盤研究 C、寺崎順、研究代表者、2014 年度採択、- 69 -

筑波大学 計算科学研究センター 平成 26 年度 年次報告書

交付額（直接経費）1,300 千円、課題名 QRPA を用いたニュートリノレス二重ベー
タ崩壊の原子核行列要素計算

6. 研究業績(1) 研究論文
A)

査読付き論文1. K. Sekizawa, K. Yabana,“Particle-number projection method in timedependent Hartree-Fock theory: Properties of reaction products”, Phys. Rev.
C 90, 064614 (2014).
2. T. Inakura, W. Horiuchi, Y. Suzuki, and T. Nakatsukasa, “Mean-field analysis
of ground state and low-lying electric dipole strength in 22C”, Phys. Rev. C 89, 064316 (2014).
3. S. Ebata, T. Nakatsukasa, and T. Inakura, “Systematic investigation of lowlying dipole modes using the canonical-basis time-dependent Hartree-FockBogoliubov theory”, Phys. Rev. C 90, 024303 (2014).
4. M. Matsuo, N. Hinohara, K. Sato, K. Matsuyanagi, T. Nakatsukasa, and K.
Yoshida, “Quadrupole shape dynamics from the viewpoint of a theory of largeamplitude collective motion”, Phys. Scr. 89, 054020 (2014).
5. T. Nakatsukasa, “Finite amplitude method in linear response TDDFT
calculations”, J. Phys. Conf. Ser. 533, 012054 (2014).
6. J.Sheikh, N.Hinohara, J.Dobaczewski, T.Nakatsukasa, W.Nazarewicz, andK.Sato, “Isospin invariant Skyrme density functional approach with axial
symmetry”, Phys. Rev. C 89, 054317 (2014)
7. J. Terasaki, “Many-body correlations of quasiparticle random-phaseapproximation in nuclear matrix elements of neutrinoless double-β decay”,
Phys. Rev. C 91, 034318 (2015).
8. J. Terasaki, “Relation between pairing gaps and transition probabilities in
132,136Te”, Journal of Physics: Conference Series 533, 012059 (2014).

9. Nobuo Hinohara and Jonathan Engel, “Proton-neutron pairing amplitudes as a
generator coordinate for double-beta decay”, Phys. Rev. C 90, 031301(R) (2014).

- 70 -

筑波大学 計算科学研究センター 平成 26 年度 年次報告書Ｂ） 査読無し論文
1. H.Z.Liang, J. Meng, T. Nakatsukasa, Z. M. Niu, P. Ring, X. Roca-Maza, N. Van
Giai, P. W. Zhao, “Nuclear charge-exchange excitations in localized covariant
de sity functional theory”,EPJ Web Conf. 66, 02064 (2014)
2. S.Ebata, T.Nakatsukasa, “Pairing effects in nuclear fusion reaction”, J. Phys.
Soc. Conf. Proc. 1, 013038 (2014)
3. K. Sekizawa, K. Yabana, “Strong Orientation Dependence of Multinucleon
Transfer Processes in 238U+124Sn Reaction”, 2nd Conf. on Advances inRadioactive Isotope Science, JPS Conference Proceedings Vol. 6 (in press).
4. K. Sekizawa, K. Yabana, “Time-dependent Hartree-Fock calculations for multinucleon transfer processes: Effects of particle evaporation on production crosssections”, VI Int. Conf. Fusion14, EPJ Web of Conference 86, 00043 (2015).
(2) 国際会議発表
A)

招待講演

1. K. Yabana, “Imaginary time formalism of triple-alpha reaction”, PKU-CUSTIPEN
Nuclear Reaction Workshop “Reactions and Spectroscopy of Unstable Nuclei”,
Peking Univ., China, Aug. 11-14, 2014.
2. K. Yabana, “Cluster Structure of Light Nuclei Superposing Multiple Slater
Determinants”, Int. Conf. Nuclear Theory in the Supercomputing Era -2014(NTSE-2014), Pacific National Univ. Khabarovsk, Russia, June 23-27, 2014.
3. T. Nakatsukasa,“Mean-field calculations for IoI”, RIBF Discussion Plus!: Island of
Inversion, Wako, Japan, Apr. 25, 2014.
4. T. Nakatsukasa,“Nuclear response and equation of state”, APCTP Workshop on
the intersection of cold-atomic and nuclear physics, Pohang, Korea, May 12-13,
2014.
5. T. Nakatsukasa,“Nuclear structure studies with energy density functionals”,
International workshop on Progress in nuclear shell-model calculations in CNSRIKEN collaboration, Wako, Japan, Nov. 26-28, 2014.
6. T. Nakatsukasa,“Isospin invariant energy density functional and isobaric
analogue states”, International Symposium on Physics of Unstable Nuclei 2014,
Ho Chi Minh City, Vietnam, Nov. 3-8, 2014.
- 71 -

筑波大学 計算科学研究センター 平成 26 年度 年次報告書

7. T. Nakatsukasa, “Time-dependent density functional calculation of nuclear
response functions”, International Conference: Nuclear Theory in the
Supercomputing Era – 2014, Khabarovsk, Russia, June 23-27, 2014.
8. T. Nakatsukasa, “Time-dependent approaches to nuclear many-body dynamics”,
International Workshop on New Frontier of Numerical Methods for Many-Body
Correlations, Tokyo, Japan, Feb. 18-21, 2015.
9. J. Terasaki, “Effects of QRPA correlations on nuclear matrix elements ofneutrinoless double-beta decay through overlap matrix”, International Conferenceon Nuclear Theory in the Supercomputing Era (NTSE2014), Khabarovsk, Russia,
June 23‒27, 2014.
10. J. Terasaki, “Many-body correlations of QRPA in nuclear matrix elements of 0vββ
decay”, EMMI Rapid Reaction Task Force on Nuclear Matrix Elements,
Darmstadt, Germany, November 10‒21, 2014.
11. J. Terasaki, “Many-body correlations of QRPA in nuclear matrix elements of 0vββ
decay”, Double-beta-decay meeting, Chapel Hill, USA, March 9 and 10, 2015.
12. N. Hinohara, “Finite-amplitude method for low-lying collective modes and QRPA
sum rules”, ICNT workshop “Physics of exotic nuclei: Theoretical advances andchallenge”, RIKEN, Wako, Japan, Jun. 9-13, 2014.
13. N. Hinohara, “Role of Fluctuations of collective coordinates in nuclear matrix
elements”, Mini-collaboration meeting on double-beta decay, Univ. of North
Carolina, Chapel Hill, NC, USA, Mar. 9-10, 2015.
B)

一般講演

1. K. Yabana, T. Akahori, Y. Funaki, “Imaginary time approach for reaction rate of
triple-alpha process”, Fourth Joint Meeting of the Nuclear Physics Divisions of
the American Physical Society and the Physical Society of Japan, Hilton Waikoloa
Village, Hawaii, USA, Oct. 7-11, 2014.
2. J. Terasaki, “Effects of QRPA correlations on nuclear matrix elements ofneutrinoless double-beta decay through overlap matrix”, Neutrino Nuclear
Responses for Neutrino Studies in Nuclei (NNR14), Ibaraki, Japan, November 5‒
6, 2014.
3. J. Terasaki, “Effect of QRPA correlations on nuclear matrix element of
neutrinoless double-beta decay through overlap of QRPA states”, 4th joint- 72 -

筑波大学 計算科学研究センター 平成 26 年度 年次報告書

meeting of the APS division of nuclear physics and the PSJ, Waikoloa, USA,
Sep.7‒11, 2014.
4. N. Hinohara and J. Engel, “Effect of Fluctuations of Quadrupole Deformation and
Neutron-Proton Correlations on Double Beta Decay Nuclear Matrix Element”,
2nd Conference on Advances in Radioactive Isotope Science (ARIS2014), Univ. ofTokyo, Jun. 1-6, 2014.
5. N. Hinohara, M. Kortelainen, W. Nazarewicz, and E. Olsen, “Low-energy
Collective Modes and Sum Rules with the Complex-energy Finite-amplitude
Method”, (poster presentation), 2nd Conference on Advances in Radioactive
Isotope Science (ARIS2014), Univ. of Tokyo, Jun. 1-6, 2014
6. N. Hinohara and J. Engel, “Effect of Fluctuations of Quadrupole Deformation and
Neutron-Proton Correlations on Double Beta Decay Nuclear Matrix Element”,
2nd Conference on Advances in Radioactive Isotope Science (ARIS2014), Univ. ofTokyo, Jun. 1-6, 2014.
7. N. Hinohara, “QRPA calculations using the complex-energy finite amplitude
method”, NUCLEI SciDAC collaboration meeting, Santa Fe, NM, USA, Jun. 1519, 2014.
8. N. Hinohara, M. Kortelainen, W. Nazarewicz, and E. Olsen, “Finite-amplitude
method for discrete collective excited states and sum rules”, Fourth Joint Meetingof the Nuclear Physics Divisions of the American Physical Society and The
Physical Society of Japan (HAWAII2014),

Waikoloa Village, HI, USA, Oct. 7-11,

2014.
9. T. Akahori, Y. Funaki, K. Yabana, “3α→12C reaction rate in imaginary time
theory”, 3rd Int. Workshop on “State of the Art in Nuclear Cluster Physics”
SOTANCP3, Kanto Gakuin Univ, Yokohama,, May 26-30, 2014.
10. K. Sekizawa, K. Yabana, “Transfer dynamics in the TDHF theory deduced fromparticle-number projection method”, ECT* Workshop on “From nuclear structure
to particle-transfer reactions and back II”, ECT*, Trento, Italy, Nov. 10-14, 2014.
11. K. Sekizawa, K. Yabana, “Strong orientation dependence of multinucleon
transfer processes in 238U+124Sn reaction”, The 4th Joint Meeting of the Nuclear
Physics Division of the APS and the JPS, Hilton Waikoloa Village, Hawaii, USA,
Oct. 7-11, 2014.

- 73 -

筑波大学 計算科学研究センター 平成 26 年度 年次報告書

12. K. Sekizawa, K. Yabana, “Time-dependent Hartree-Fock Calculation for
Multinucleon Transfer Processes”, The 2nd conference on Advances inRadioactive Isotope Science “ARIS2014”, ITO International Research Center,
Tokyo, Japan, 1-6 June, 2014
(3) 国内学会・研究会発表
A)

招待講演

1. 矢花一浩、“トリプルアルファ反応率の量子力学計算”、第２７回理論懇シンポジウ
ム「理論天文学・宇宙物理学と境界領域」、国立天文台三鷹キャンパス、2014 年 12
月 24-26 日
2. 中務

孝、“原子核密度汎関数理論の概要”、サマースクール「クォークから超新星爆発まで」、京都大学基礎物理学研究所、京都、2014 年 7 月 24 日

B)

その他の発表

1. 矢花一浩、“量子多体問題と密度汎関数理論の考え方”、サマースクール「クォークから超新星爆まで」−基礎物理の理想への挑戦−、京大基研、京都市左京区、2014 年 7
月 24 日2. 中務

孝、“核構造計算による核反応モデルの高精度化”、ImPACT 藤田プログラム

全体会議、JＳＴ別館、東京都千代田区、2015 年 3 月 26 日
3. K. Sekizawa, K. Yabana, “Time-dependent Hartree-Fock calculations for multinucleon transfer and quasi-fission processes”, A seminar at the ANU,
Australian National University, Canberra, Australia, July 21, 2014.
4. K. Sekizawa, K. Yabana, “Time-dependent Hartree-Fock calculations for
multinucleon transfer processes”, 262th Sendai Nuclear Science Colloquium,
Tohoku University, Sendai, Miyagi, Japan, April 15, 2014.
5. 関澤一之、矢花一浩、“射影演算子を用いた核子移行 TDHF 波動関数の分析 II”、日本物理学会第 70 回年会、早稲田大学、東京都新宿区、2015 年 3 月 21 日-24 日
6. 関澤一之、矢花一浩、“原子核ダイナミクスの微視的シミュレーション”、第 6 回「学際計算科学による新たな知の発見・統合・創出」シンポジウム －HA-PACS とCOMA による計算科学の発展と、分野融合への取り組み－、筑波大学大学会館国際会
議室、茨城県つくば市、2014 年 10 月 21 日- 22 日
7. 橋本幸男、“ラグランジュ格子を用いた Gogny-TDHFB による２０O＋２０O の計算”、日本物理学会第 70 回年次大会、東京都新宿区、2015 年 3 月 21 日-24 日- 74 -

筑波大学 計算科学研究センター 平成 26 年度 年次報告書

8. 橋本幸男、矢花一浩、関澤一之、“密度汎関数理論に基づく原子核ダイナミクスの研究”、第 1 回「京」を中核とする HPCI システム利用研究課題 成果報告会、コクヨホ
ール、東京都港区、2014 年 10 月 31 日

(4) 著書、解説記事等
1. 日野原伸生、“ノースカロライナ大学チャペルヒル校 (海外通信)”、原子核研究、第 59 巻 1 号(2014 年 9 月発行) pp. 14-15.

7.

異分野間連携・国際連携・国際活動等【国際連携】
1. 日本学術振興会二国間交流事業オープンパートナー・セミナー（代表：中務）の支援により、米国ロスアラモス国立研究所の Joe Carlson 氏を中心とする理論核物理研究者等とハワイにおいて共同セミナーを開催した（2014 年 10 月）（中務）。
2. ポーランド・ワルシャワ大学の原子核理論グループと共同で、オープンソースコード
HFODD をアイソスピン不変なエネルギー密度汎関数へ適用する拡張に取り組んでいる（中務）。
3. 韓国・APCTP およびイタリア・トレント大学の冷却原子系理論グループと、1 次元周期ポテンシャル中のフェルミ粒子系の超流動状態に関する共同研究を行っている（中務）。
4. 米国ノースカロライナ大学の Engel 教授と二重ベータ崩壊の核行列要素に関する共同
研究 (日野原)。
5. 米国ミシガン州立大学 Nazarewicz 教授およびフィンランド・ユバスキュラ大学のKortelainen 研究員と原子核密度汎関数の諸問題に関する共同研究 (日野原)
8. シンポジウム、研究会、スクール等の開催実績
1. 国際会議 2nd International Conference on Advances in Radioactive Isotope
Science (Tokyo, June 1-6, 2014)の組織委員を務めた（中務）。
2. 国際 ICNT ワークショップ Physics of exotic nuclei: Theoretical advances andchallenges (RIKEN, June 9-13, 2014)の組織委員長を務めた（中務）。
3. HPCI 戦略プログラム 5 共催の国際ワークショップ Advances and perspectives incomputational nuclear physics (Waikoloa, HI, USA, Oct. 5-7, 2014)の組織委員長を務めた（中務）。
4. 第 13 回国際サマースクール CNSSS14 (RIKEN, Aug. 21-27, 2014) の組織委員を務めた（中務）。
- 75 -

筑波大学 計算科学研究センター 平成 26 年度 年次報告書

9. 管理・運営
矢花一浩
矢花は、センターの共同研究担当主幹として、当センターの全国共同利用業務である学際共同利用プログラムの運営を統括した。また、数理物質系物理学域長・
数理物質科学研究科物理学専攻長を務めた。
中務孝
計算科学研究センター 宇宙・原子核物理研究部門
原子核物理分野リーダー
計算科学研究センター 運営委員会委員
計算科学研究センター 共同研究委員会委員
数理物質系物理学域 運営委員会委員
数理物質系物理学域 原子核理論グループ長
数理物質系物理学域 図書委員長
理化学研究所・中務原子核理論研究室 准主任研究員（兼務）
京都大学・基礎物理学研究所 運営協議会委員
核理論委員会委員
日本物理学会 第 67-70 期代議員
理化学研究所 RI ビームファクトリー 国際プログラム諮問委員会(NPPAC)委員
雑誌「原子核研究」編集委員

寺崎 順
文部科学省 HPCI 戦略プログラム分野 5「物質と宇宙の起源と構造」中の計算機ユーザーのプログラミング支援活動において連絡係を務めた。
10. 社会貢献・国際貢献
中務孝
理化学研究所一般公開(2014.4.19)において講演「量子力学と原子核」
カナダ・TRIUMF

国際プログラム諮問委員会(SAP-EEC)委員

Editor for Journal of Physical Society of Japan
Editor for International Journal of Modern Physics E

11. その他- 76 -

筑波大学 計算科学研究センター 平成 26 年度 年次報告書

海外長期滞在1) 日野原 伸生, Univ. North Carolina at Chapel Hill, NC, USA,2014 年 6 月

20 日〜2014 年 8 月 25 日
2) 日野原 伸生, National Superconducting Cyclotron Laboratory,
Michigan State Univ., East Lansing, MI, USA, 2014 年 8 月 27 日〜2015 年
3 月 31 日(次年度に継続)- 77 -