I. 宇宙・原子核物理研究部門
I-1. 原子核物理分野
1. メンバー
教授

矢花一浩、中務	
  孝

准教授

寺崎	
  順（HPCI 戦略プログラム）

講師

橋本幸男

助教

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

研究員

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

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

2. 概要
	  核子（陽子・中性子）の多体系である原子核の構造・反応・応答などの多核子量子ダイナミ
クスの研究を推進している。安定線（ハイゼンベルグの谷）から離れた放射性アイソトープの原子核の構造と反応、エキゾチックな励起状態の性質、様々な集団運動の発現機構など、未解
決の謎の解明に取り組んでいる。原子核の研究は、フェルミ粒子の量子多体系計算という観点
で、物質科学や光科学、冷却原子系の物理と密接なつながりをもつ。また、クォーク・グルー
オンのダイナミクスを記述する格子 QCD に基づく核力の計算、軽い原子核の直接計算などが進展する中、素粒子物理学との連携も重要性が増している。ニュートリノの解明に向けたニュー
トリノレス二重ベータ崩壊の観測実験や、素粒子標準模型のテストに関わる実験などにも原子
核理論の精密計算が不可欠とされている。また、元素の起源や星の構造にも原子核の性質は深く関わり、宇宙物理学とも密接に関係している。本部門・原子核物理分野のメンバーは、このような幅広い課題に取り組み、分野の枠を超えた研究を推進している。	
 	
 3. 研究成果【	 1	
 】時間依存ハートレー・フォック法による核子移行過程の記述（矢花、関澤）	
 多核子移行反応は、核子多体系に発現する非平衡量子輸送過程としてその反応機構に興味が持たれており、また近年、多核子移行反応を用いこれまでに生成されたことのない不安定原子
核を生成する可能性が盛んに議論されている。我々は、多核子移行反応の微視的反応機構を明らかにすると共に、不安定核の生成方法を予言することを目指し研究を進めてきた。昨年度までに、TDHF 法に粒子数射影法を組み合わせることにより、起こる確率の大きい移行反応過程の断面積を定量的に記述できることを明らかにした。しかし、移行する核子数が大きく、起こる
確率の小さい移行反応過程については、実験値を過小評価するという問題に直面した。そこで本年度には、この実験値と食い違いを解消することを試みた。相違の原因として、励起した原子核からの核子放出による脱励起（核子蒸発）の効果が考えられる。核子蒸発過程は、我々の数値シミュレーションで扱う時間スケールより遥かに長い時間スケールで起こる現象であり、TDHF 計算とは別に統計模型を用いた計算を行う必要がある。統計模型を利用するために、我々は粒子数射影法を応用し、生成された原子核の励起エネルギーを評価する方法を開発した。そして、核子蒸発の効果を取り入れた断面積を評価し、実験値との比較を行った。その結果、核fPDand short principal axes of deformation. On the other hand,for the low-lying E1 mode, the geometrical effect is not2
prominent. The authors of Ref. [74] have reported that the
deformation hinders the PDR strength for the prolate neutron1
rich Sn isotopes. They also found that the K = 0 componentof the calculated low-lying E1 strength is larger than the0
K = 1 component, which was interpreted by the effect of40
anisotropic neutron skins. Namely, the neutron skin thickness
子蒸発の効果が実験値との相違を改善することを確認したが、効果は不十分であった。このこin the z (K = 0) direction is expected to be thicker than in
FIG. 13. (Color onとは、TDHF 計算で扱えていない相関の影響を示唆している。
the x-y (K = 1) directions. In this section, we investigate thisarrows indicate the ob
deformation effect for selected isotopes: strongly deformed(Sr and Zr) and weakly deformed nuclei (Pd and Cd).
state, the PDR stre【	 2	 】正準基底時間依存ハートレー・フォック・ボゴリューボフ理論による多核子ダseen in N = 72, theA. Deformation and PDR strength
イナミクス計算（ 中務、江幡（北大）、稲倉（京大））	
 going back to the sp
In 3Fig.
12, fPDR are shown for Zr isotopes. In the present
In weakly deform	  原子核理論分野で発展が目覚ましい
次元実空間表示の時間依存ハートレーフォック法に
calculation, Zr isotopes with N = 60–72 are deformed with a
PDR is much mild
対して、計算コストの増加をミニマムに抑えた上で、核子対凝縮と超流動ダイナミクスを扱え
prolate shape, and the one with N = 74 has a small triaxial
ground state of Pd
shape. The other nuclei are calculated to be spherical in the
るように拡張した理論を、2010 年に我々が提唱した。このときの論文に対して、第一著者でat N = 58–74; howground state. Here, we can clearly identify the deformation
ある江幡氏（当時筑波大学大学院生・理化学研究所 JRA）が、本年度の日本物理学会若手奨励 significantly smalle
effects. The total PDR strength suddenly reduces at N = 58 →gradually increases
賞を受賞した。この方法を線形応答領域に応用し、原子核の応答関数の系統的計算を行い、軽60 at the onset of deformation (β2 ≈ 0.37). This is consistentthat. In Cd isotopes
with Ref. [74], in that the deformation
the PDR
Pd.
い原子核から希土類領域の重い原子核まで、中性子分離エネルギーが2MeVhinders
以上のアイソトー
strength. However, at N = 72 → 76 where the deformation reWe have foundプについて電気双極子(E1)強度分布を調査した。その結果、中性子過剰核の広い領域においてduces and disappears, the PDR strength decreases again. Thus,
strength is consisten
低エネルギーの E1 ピーク（ピグミー共鳴）が現れ、その出現のメカニズムに中性子の殻効果
the deformation does not always hinder the PDR strength.
potentials. In Figs.
Very similar強度のほとんどは巨大共鳴状態に集中するた
behaviors can be found for Sr isotopes as well,−λn are shown foが重要であることが分かった。また、通常、E1
in Fig. 13. In Sr isotopes, the ground state has a prolate shapeisotopes, the sudd
め、低エネルギー領域の E1 強度は大幅に弱められていることが知られているが、非常に中性with β2 ≈ 0.37 for N = 60–72, and an oblate shape at N = 38
causes the change
and 74. Whenthe prolate deformation
develops in the ground
子過剰な重い原子核の励起エネルギー5	
 MeV
以下には、このE1 抑制が働かない新種のピグミ|λn | monotonically
exceptfor two placeー共鳴が現れることを予言した。さらに、これらのピグミー共鳴の性質と中性子星の構造を決
place at N = 58 →定する状態方程式との関係を
9numbers perfectly精査した。	 PDR strength in n
8
Zr isotopes
7

た核融合反応計算を開始し、6軽い原子核に限定的ではある5

味深い結果を得た。衝突エネ

6
5

4

4

3

ルギーを同じにして、対相関2を入れた計算と無視した計算

1とを比べると、対相関によっ

0て核融合が阻害されるケース

Total
z-axisy-axisx-axisfPDR [%]

が、対相関の影響に関して興fPDR [%]

	  また、同じ計算手法を用い

40

50

60

70

Neutron Number

82

90があることが分かった。さら

図１：Zr
アイソトープの低エネルギーE1強度の割合。特定のFIG. 12.
(Color online) The PDR strength in
Zr isotopes decomposed
into
x,
y,
and
z
components
which
are
denotedby red, blue,
中性子数 N=50, 60, 72, 82 で傾きが突然変化する。
研究を続行中である。	
 and green bars. In spherical nuclei, the three components are equal
to each other. In this figure, we include isotopes with the neutron	 chemical potentials smaller than 2 MeV (N > 82).
【	 3	 】大振幅集団運動理論を用いた核反応ダイナミクスの記述（温、中務）	
 にこの効果を確認するため、024303-8	  線形領域を超える大振幅集団運動を扱う理論として、断熱近似型の理論が提案されているが、そのうちの一つが、2000 年に提案された断熱自己無撞着集団座標法(Adiabatic	 Self-consistent	
 Collective	 Coordinate	
 Method:	
 ASCC 法)である。この理論では、少数自由度の集団空間（座標）の自己無撞着な抽出が可能であり、特に、過去の断熱時間依存平均場理論では不可能であった一意的な抽出ができる点が優れている。この理論に基づいて、低エネル
ギーの多核子反応ダイナミクスを記述することを目指した研究を開始した。まず、テストとして、比較的簡単な BKN 相互作用と呼ばれる有効相互作用を用いて、3 次元座標表示の TDHF を基

Pd isoTotal
z-axisy-axisx-axis

3
2
1
0

50

FIG. 14. (Color o礎にしたコード開発を行い、2 つのアルファ粒子の散乱・融合のミクロな機構を明らかにすることにした。拘束演算子を局所調和方程式と呼ばれる式から求め、それを拘束条件付き平均場計算に応用する。これを自己無撞着に計算し、核反応の集団運動を支配するポテンシャルと質量パラメータを求める。これまでに、このようにして自己無撞着に決めた質量が遠方で換算質量に一致すること、また 2 つのアルファ粒子が接触した後には換算質量が大きくなり、反応速
度が弱められることなどが分かっている。	
 
	
 
【	
 4	
 】アイソスピン不変なエネルギー汎関数による原子核密度汎関数計算（中務、佐藤（理研）、Dobaczewski（ワルシャワ大）、	
 Satula	 (ワルシャワ大））	
 	
  現在主流となっている原子核のエネルギー密度汎関数は、Skyrme 形式、Gogny 形式、共変形()

式（相対的）の 3 つに大別されるが、どれも陽子と中性子の密度 ρ p , ρ n の汎関数としてエネ
ルギーが与えられている。
しかし、陽子や中性子はアイソスピンの第 3 成分の固有状態であり、アイソスピン空間における回転に対して不変ではなく、一般にはアイソスピンが任意の方向を向いた状態、すなわち陽子と中性子が混合した状態に拡張する必要がある。これを実行するため、陽子・中性子を区別せずに「核子」として扱う新しい Kohn-Sham 方程式と、それに対応す()

る非対角要素 ρ pn , ρ np を含むエネルギー汎関数を構築し、その計算コード開発を実施した。
アイソスピン不変に拡張することによる新しいパラメータの導入はなく、これまでのエネルギー汎関数のパラメータセットをそのまま利用することができる。計算コードは、ポーランド・ワルシャワ大学のグループを中心に開発されている HFODD にこの拡張を適用することにした。
この研究はワルシャワ大学との共同研究で進められている。	
 	
  この拡張による大きなメリットとして、荷電類似状態（アイソバリックアナログ状態）と呼ばれる状態が、自然かつ簡便な形で表現できることがあげられる。図には、質量数 A=54のアイソスピン T=1 の状態が示されている。左端の点(T=−1)
が 54Ni に、右端(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)を用いて和則を計算すると、大次元の行列の対角化とそのすべての解が必要となり、さらに計算量による制限のため行列の次元を減らす近似が必須である。全自由度を扱うことが出来る有限振幅法を用い、複素エネルギー平面で励起状態に対応する極をすべて囲む複素積分を行うことによって、和則を効率よく計算する
定式化を行った。全自由度が扱える小さな模型空間において QRPA の行列対角化による解と一致することを数値的に確認し、さらに和則を与える Thouless 定理や絶縁定理の密度汎関数理論での有効性を、現実的な場合において数値的に示した。有限振幅法による本定式化はこれらの定理が使えない場合においても有効であり、複素積分の並列化によって和則を効率よく正確に計算することが可能となった。	
 
	
 
【	
 7	
 】二重ベータ崩壊の核行列要素（陽子−中性子対振幅の量子ゆらぎを取り込んだ核
行列要素の計算	 (日野原、Engel(Univ.	
 North	 Carolina))	
 	
  二重ベータ崩壊の核行列要素はアイソスカラー型の陽子−中性子対相関力によって抑制されることが乱雑位相近似(QRPA)計算によって知られているが、この陽子−中性子対相関力は基底状態に質的な影響を与えないため、その強さについてはよくわかっていない。QRPA 計算は相転移を記述出来ないため、基底状態の四重極変形ゆらぎが大きく変形共存や変形転移状態にある
場合や、アイソスカラー対凝縮相への転移状態近傍にある場合は近似が破綻する。この QRPAの問題を回避し核行列要素の相互作用依存性を明らかにするために、四重極変形度と陽子−中性子対振幅を拘束した変分によってこれらの量子ゆらぎの大きい状態をも含んだ基底を生成し、量子数射影後に重ねあわせることで、崩壊の始状態と終状態を構築した(生成座標法)。特に陽子−中性子対を自由度に取り込んだ生成座標法はこれまでに行われたことのないユニークなものである。第一の計算例として、二主殻模型空間で 76Ge→76Se の二重ベータ崩壊の核行列
要素を評価した。生成座標法では殻模型計算より広い一粒子模型空間を用いることが出来、さらに二準粒子励起に制限される QRPA よりも豊富な多体相関を取り込むことが出来るため、適切な有効相互作用を用いることで、将来的には殻模型計算と QRPA のどちらにも含まれていな
い効果を議論することが期待できる。もう一つの将来計画は Skyrme 型の原子核密度汎関数に基づいて陽子−中性子対の自由度を入れた生成座標法計算を行うことで、β崩壊などの関連する物理量も含めた荷電交換過程の系統的な計算へと発展させることである。	
 
	
 
【	
 8	
 】Gogny-TDHFB による ２０Ｏ- ２０Ｏ正面衝突の計算	  （橋本）	
 	  原子核の振る舞いを研究するうえで、時間依存平均場の方法は定性的にも定量的にも極めて有用な枠組みとなっている。原子核の基底状態近傍の状態においては対相関が重要な働きをすることが知られているので、中性子過剰核の性質を理解する目的で、TDHF の枠組みを拡張して対相関を扱えるようにした時間依存ハートレーフォックボゴリュボフ（TDHFB）の実用的な数値計算が行われるようになってきた。橋本は、調和振動子基底と空間格子としてのLagrange 格子点を組み合わせた基底を利用し、Gogny 力を用いた TDHFB を解く、という数
値計算の方法を整備してきた。今年度は、２０Ｏ+２０Ｏの衝突に応用した。ちなみに、TDHF+BCSの形式での原子核衝突の計算は存在するが、TDHFB 方程式を直接解いて原子核衝突を計算する試みは、現時点では極めて少ない。計算に必要な格子空間の選定は初期条件の設定と関係し、そのための試行錯誤に計算時間と労力を割く必要があった。実際の計算は、クーロン障壁頂上付近の軌道を３本（３個の初期エネルギー）設定し、２１fm 離れた点から計算を開始した。
この一連の計算により、① 障壁より上のエネルギーの軌道では融合が起こり、核子の移行量の期待値を計算できた。② 障壁より下の軌道において、対相関の働く中性子側では、陽子側の約１０倍の移行数があり、これは２０Ｏの中性子数・陽子数比である１．５よりずっと多いことがわかる。さらに、③ 障壁よりも上のエネルギーを持つ２本の軌道から、巨視的な摩擦係数を計算することができた。摩擦係数は、融合の初期に急速に増加することが示された。今後、これらの量の初期エネルギー依存性を系統的に調べていく。
【 9 】Energy	 dependence	
 of	
 the	
 nucleus-nucleus	
 potential	
 and	
 the	
 friction	 parameter	 in	
 fusion	 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, the friction parameter, and the random
force characterizing a one-dimensional Langevin-type description of the heavy-ion fusion process areinvestigated. Systematic calculations with the ImQMDmodel show that the fluctuation-dissipation
relation found in symmetric head-on 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 is transferred
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 the dissipative 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 fusion potential andexplains 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 年度採択、	
 
交付額（直接経費）1,300 千円、課題名	
  QRPA を用いたニュートリノレス二重ベータ崩壊の原子核行列要素計算	
 6. 研究業績(1) 研究論文

A)

査読付き論文1. K. Sekizawa, K. Yabana,"Particle-number projection method in time-dependent Hartree-Fock theory: Properties of reaction
products", Phys. Rev. C90, 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, 0 64316 (2014).
3. S. Ebata, T. Nakatsukasa, and T. Inakura,“Systematic investigation of low-lying dipole modes using the canonical-basis time-dependent
Hartree-Fock-Bogoliubov 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 large-amplitude 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, K.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-phase approximation in nuclear matrix elementsof neutrinoless double-β decay”, Physical Review 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).

Ｂ） 査読無し論文
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", 2ndConf. on Advances in Radioactive Isotope Science, JPS Conference Proceedings Vol. 6 (inpress).
4.	
 K. Sekizawa, K. Yabana,

"Time-dependent Hartree-Fock calculations for multi-nucleon transfer processes: Effects ofparticle evaporation on production cross sections", 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 andnuclear physics, Pohang, Korea, May 12-13, 2014.
5. T. Nakatsukasa,“Nuclear structure studies with energy density functionals”, International workshop on Progress innuclear shell-model calculations in CNS-RIKEN collaboration, Wako, Japan, Nov. 26-28, 2014.
6. T. Nakatsukasa,“Isospin invariant energy density functional and isobaric analogue states”, International Symposiumon Physics of Unstable Nuclei 2014, Ho Chi Minh City, Vietnam, Nov. 3-8, 2014.
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 of neutrinoless double-beta decay throughoverlap matrix”, International Conference on 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 and challenge”, 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 of neutrinoless 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 meeting of the APS division of nuclear physics andthe 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 inRadioactive Isotope Science (ARIS2014), Univ. of Tokyo, 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 inRadioactive Isotope Science (ARIS2014), Univ. of Tokyo, Jun. 1-6, 2014.
7. N. Hinohara,“QRPA calculations using the complex-energy finite amplitude method”, NUCLEISciDAC collaboration meeting, Santa Fe, NM, USA, Jun. 15-19, 2014.
8. N. Hinohara, M. Kortelainen, W. Nazarewicz, and E. Olsen, “Finite-amplitude method for discrete
collective excited states and sum rules”, Fourth Joint Meeting of the Nuclear Physics Divisions of theAmerican 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 from particle-number
projection method", ECT* Workshop on "From nuclear structure to particle-transfer reactions andback 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.
12. K. Sekizawa, K. Yabana, "Time-dependent Hartree-Fock Calculation for Multinucleon Transfer
Processes", The 2nd conference on Advances in Radioactive 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 multi-nucleon transfer andquasi-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 日、早稲田大学早稲田キャンパス）	
 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	
 and	
 challenges	
 (RIKEN,	
 June	
 9-13,	
 2014)の組織委員長を務めた（中務）。	
 3. HPCI 戦略プログラム 5 共催の国際ワークショップ	
 Advances	
 and	
 perspectives	
 in	
 computational	
 nuclear	 physics	
 (Waikoloa,	
 HI,	
 USA,	
 Oct.	
 5-7,	 2014)の組織委員長を務めた（中務）。	
 4. 第 13 回国際サマースクール	
 CNSSS14	
 (RIKEN,	
 Aug.	
 21-27,	 2014)	
 の組織委員を務めた（中
務）。	
 
9.	
 管理・運営	 矢花一浩
矢花は、センターの共同研究担当主幹として、当センターの全国共同利用業務である
学際共同利用プログラムの運営を統括した。また、数理物質系物理学域長・数理物質科学研究科物理学専攻長を務めた。	
 
中務孝	
 
	
  	
  	
  	
  	
  計算科学研究センター	  宇宙・原子核物理研究部門	
  原子核物理分野リーダー	
 
	
  	
  	
  	
  	
  計算科学研究センター	  運営委員会委員	
 
	
  	
  	
  	
  	
  計算科学研究センター	
  共同研究委員会委員	
 
	
  	
 	
 	
 	
 	
 	
 	
 	
 数理物質系物理学域	
  運営委員会委員	
 
	
  	
 	
 	
 	
 	
 	
 	
 	
 数理物質系物理学域	
  原子核理論グループ長	
 
	
  	
 	
 	
 	
 	
 	
 	
 	
 数理物質系物理学域	
  図書委員長	
 
	
  	
 	
 	
 	
 	
 	
 	
 	
 理化学研究所・中務原子核理論研究室	
  准主任研究員（兼務）	
 
	
  	
 	
 	
 	
 	
 	
 	
 	
 京都大学・基礎物理学研究所	
  運営協議会委員	
 
	
  	
 	
 	
 	
 	
 	
 	
 	
 核理論委員会委員	
 
	
  	
 	
 	
 	
 	
 	
 	
 	
 日本物理学会	
  第 67-70 期代議員	
 
	
  	
 	
 	
 	
 	
 	
 	
 	
 理化学研究所 RI ビームファクトリー	
  国際プログラム諮問委員会(NP-PAC)委員	
 

	
  	
  	
 	
 	
 	
 	
 	
 	
 	
 雑誌「原子核研究」編集委員	
 
	
 	
 	
 	
 	
 	
 	
 寺崎	
  順	
 
	
  	
  	
 	
 	
 	
 	
 	
 	
 	
 文部科学省 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.	
 その他	
 
海外長期滞在	
 
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 日(次年度に継続)