III. 物質生命研究部門
III-1 原子核理論グループ
１．メンバ教授

矢花 一浩

講師

橋本 幸男、中務 孝 （平成１９年７月３１日まで）

２．概要
当グループでは、強い相互作用により束縛した量子多体系である原子核、そしてフェルミ多粒子系として共通する理論や方法が有効な多電子系としての原子・分子・クラスター・固体などの物質科学を対象に、特に量子ダイナミクスに対する計算科学を特色として研究に取り組んでいる。原子核分野では、最近発展が著しい不安
定原子核の構造と反応の分野で、物質科学分野では、強レーザー場と物質の相互作用で起こる様々な電子ダイナミクスの解明を中心的な課題とし、時間依存密度汎関数理論を物質相にまたがる共通の理論として用い、研究の展開を図っている。
３．研究成果【１】原子核集団運動の理論、不安定核の構造（１）Skyrme 汎関数による光吸収断面積の実時間 TDDFT 計算(中務、矢花)時間依存密度汎関数理論(TDDFT)に基き、実験データが存在する安定核の光吸収断面積のエネルギー依存性を非経験的手法で計算した。光子エネルギーが中性子(陽子)分離エネルギーを超えると、原子核からの核子放出が可能になり、強度分布は連続分布となる。このような核子の連続状態の取扱には、これまで我々が開発してきた吸収境界条件法を用いた。ヘリウムからカルシウムまでの比較的軽い原子核について、実験データと比較した。その結果、断面積のエネルギー依存性について定性的に再現する一方で、強度のピーク位置が系統的に過小評価することが分った。重い原子核の光吸収が対称エネルギーで主に規定されているのに対して、軽い核では、それ以外の相関が重要であることを示唆していると考えられる。
（2） 断熱自己無撞着集団座標（ASCC）法による変形共存現象の解析（中務、日野原（京都大）、松尾（新潟大）、松柳（京都大））
昨年度、数値的不安定性が理論のゲージ不変性にあることをつきとめ、ゲージ固定によってこの不安定性が解消されることが分った。今年度、ASCC 法で求めた正準集団座標、運動量によって記述される集団ハミルトニアンを量子化し、励起スペクトルと集団波動関数を求めた。その結果、波動関数の広がりが角運動量の大きさに強く依存することが分った。例えば、プロレート・オブレートの変形状態が１つの原子核に共存する「変形共存現象」において、角運動量がゼロでは、プロレートとオブレートの変形が強く混ざっていても、わずかに回転を与えることで、これらの変形は独自性を高め、別々の回転バンドとして存在することが示された。実験データとの比較を実行し、この予言を確かめる観測量を現在模索中である。
（3）時間依存密度汎関数理論による原子核の応答関数の系統的計算（稲倉、矢花、中務(理研)）
重い短寿命の不安定核も生成可能となり、実験でアクセスできる核種が増えている。安定核、不安定核を分け隔てる事なく全核種の励起状態を計算し、不安定核まで含めた原子核の諸性質を議論する為の系統的計算を開始した。先行する計算では対象核が球形に限られるものが殆どであり、系統的計算に適用可能な変形自由度を考慮に入れた計算コードを作成、質量数 80 前後までの原子核の電気双極子励起(E1)モードについて計算を実施した。巨大共鳴のエネルギーの質量依存性、基底状態の変形度が巨大共鳴に及ぼす影響などについて議論した。また、数種の重い核についても計算を実施し、実験と比較する事で、我々の計算の有用性を示した。これ迄に幾つかの不安定核で観測された、低エネルギー領域での E1 遷移強度が、安定線から離れた不安定核一般で現わる事を示した。また、特定の一粒子軌道がその出現に重要な役割を果している事も指摘した。
（４）BCS 形式による時間依存密度汎関数理論の拡張（江幡、中務、稲倉、橋本、矢花）
原子核を超重領域まで系統的に計算するためには、対相関を取り入れることが重要である。最も良く知られた方法としては、超伝導の理論として開発された BCS(Bardeen-Cooper-Schriefer)の理論があり、原子核の核子超流動性を扱う理論として、古くから用いられてきた。我々は、この BCS 状態を時間依存として、一粒子軌道および(u,v)因子が従うべき方程式を導き、原子核の応答計算に応用した。これまでのところ、BKN 相互作用によるアイソスカラー型の強度分布関数の計算に応用し、変形により単極子強度分布に著しい変化が現れることを示した。これは、回転対称性が破れたことにより、巨大四重極共鳴の状態と結合し、エネルギーの低い領域にピークが現れることによる。現在、より現実的な Skyrme 汎関数へのプログラムの改訂を実行中である。
（５）有限振幅法の開発（中務(理研)、稲倉、矢花）
原子核の励起状態の計算には RPA がよく使われる。しかしながら、エネルギー密度汎関数を密度の 2 階微分する事で得られる残留相互作用は、その繁雑さから簡略化される事が殆どであった。核力の自己無撞着性を満たしつつ、計算が容易な手法として、有限振幅法を開発した。波動関数の時間発展をエネルギー表示の枠組みで線形化する事で RPA と同値な方程式を導出し、残留相互作用は容易に(数値的に)計算できる手法である。この手法は平均場の計算コードに若干の修正を施すだけで良いので、コードの開発自体も容易に行なえる利点もある。この新しい手法で計算を実施し、完全に自己無撞着な RPA の結果をよく再現する事を示した。
（６）原子核の三次元的回転運動の理論（橋本、堀端（青森大））
有限量子系である原子核は様々な集団運動のモードを持っている。回転運動は、その中でも最も顕著な集団運動の例である。原子核平均場は密度分布と緊密に関係している（nuclear self-consistency）ので、原子核平均場の回転運動は内部核子によるコヒーレントな運動の現れである。原子核の回転運動の研究は主に軸対称変形をした核の主軸まわりの定常回転を対象に行われてきた。クランキング模型に基づくこの考え方は、原子核の回転スペクトルの構造とその背景にある力学について説明する上で、定性的にも定量的にも成功してきた。一方、理論的な立場からは、主軸のまわりの回転運動だけではなく、より一般的な回転運動の存在が期待されている。原子核が軸対称から離れ、三軸非対称変形をすると、“主軸まわりの定常的な回転”という基礎の上に一種のフォノンが生じたような運動モードが起こることが Bohr と Mottelson の教科書でも指摘されている。ウォブリング(wobbling)と呼ばれるこの揺動運動は回転軸が平均場の主軸から離れて才差運動のように動き出すことを意味している。本研究では、ウォブリング運動のような三次元的な回転運動が原子核において生じる力学的機構を微視的に理解することを目的としている。
今年度は特にオスミウム 182Os の J=24 の場合に、ウォブリングをはじめとする励起モードの計算を生成座標法（GCM）に基づいて進めた。J=24 の場合に特徴的なことは、角運動量の拘束を主軸から離れた値にして行う傾斜角運動量ハートレーフォックボゴリュボフ（tilted axis HFB）の計算において、対称軸に垂直なｘ軸から対称軸（ｚ軸）へと拘束角運動量を倒す角度（チルト角）に依存して HFB 解の内部構造が３種類に変化することである。ひとつは、エネルギーの局所平衡点に相当する傾斜角回転（tilted axis rotation; TAR）と呼ばれるモードを含む分枝である。この、TAR モードは、主軸回転運動を仮定した計算では出てこないモードであり、P.M.Walker らの実験で K 量子数が 8 のバンドに相当すると考えられる。これらの K=8 のバンドはｘ軸まわりの回転運動に対応するイラスト線に対して、新しいタイプのイラスト線になっている。他のひとつは、角運動量が対称軸に近い向きのときに起きる歳差運動（precession）である。さらに、これらのモードをつなぐ過渡的な領域が存在することがわかった。GCM 計算による歳差運動の励起エネルギーは、実験値に近い値が得られている。これらのモードの内部構造を解析することが次の段階になる。
（７）Gogny 力を用いた時間依存 HFB コードの開発（橋本、野出木）不安定核の励起状態の研究は、主に Skyrme 力+ゼロレンジ対相関力に基づいた準粒子乱雑位相近似（QRPA）によって行われている。この組み合わせの場合には、Skyrme 力は粒子・空孔相関のチャネルに、また、ゼロレンジ対相関力は粒子・粒子相関のチャネルに働く。対相関は、適当なカットオフを導入して実用化されている。一方、Gogny 力を用いた QRPA の場合には粒子・空孔チャネルと粒子・粒子チャネルは共通の有効核力を用いて計算される。われわれは、この自己完結性に着目し、Gogny 力に基づいた時間依存 HFB（TDHFB）の方程式を解くためのコードを開発した。目的は、Skyrme 力+ゼロレンジ対相関の計算と比較しつつ軽い不安定核から始めて不安定核の系統的な励起状態の情報を蓄積することである。微小振幅の TDHFB 方程式は QRPAと同等であるので、これは、QRPA を解く方法のひとつでもある。応用は、酸素の同位体で四重極および双極型振動の強度関数を調べ、スペクトルの位置と全強度を確かめた。また、変形核を含むマグネシウム同位体で同様に計算し、変形の効果と対相関の効果が取り入れられていることを確認した。目下、クーロン力を取り入れた計算過程全体の時間的な効率を上げること、及び、基底空間を大きく取れるように改良することが課題となっている。
【２】強レーザー場と物質の相互作用に対する第一原理計算
原子核は核子からなる凝縮したフェルミ粒子多体系であるが、物質は同じくフェルミ粒子である電子の凝縮した多体系である。このため、電子の量子状態を記述する多体理論は原子核を記述する多体理論と共通する場合が多い。我々は、原子核物理学で長い歴史を持つ時間依存平均場理論を物質科学分野で応用する研究を進めている。特に強いレーザー場中で起こる電子のダイナミクスは、様々な非線形現象が現れ興味深い。
（１）誘電体の光絶縁破壊に対する第一原理シミュレーション（矢花、中務、岩田（計科セ）、乙部（関西原研）、山極（関西原研）、Ｂｅｒｔｓｃｈ（ワシントン大）
透明な絶縁体に可視領域の強いレーザー光を照射すると不可逆な損傷が生じ、レーザーを用いた加工に利用される。この過程の最初のステップでは、レーザー光による電子励起とそれに伴う光絶縁破壊が起こると考えられている。この光絶縁破壊に対し、時間依存密度汎関数理論を用いた第一原理シミュレーションを行っ

ミクスを記述する枠組みは、以前実時間計算から誘電関数を求める際に構築した。その基礎方程式は、外場と分極場を表す空間的に一様なベクトルポテンシャルを含む時間依存 Kohn-Sham 方程式である。外場のポテンシャルとして一様なパルスレーザーを表すものを用い、ブロッホ関数の時間変化を計算することにより、レーザー場に起因するダイナミクスを記述することができる。

-1

周期系に一様な外場が働く場合の電子ダイナ

Excitation energy [eV atom ]た。
110101010Frequency 3.1 eV-1-3TDDFT2ΔE∝CI-5Keldysh theory

10-710121013101415102

Irradiance [W/cm ]1610右図は、ダイアモンドに強レーザー場を照射した場合の計算結果を示す。横軸にレーザー場の最大強度、縦軸にレーザー場から電子へのエネルギー移行をプロットしている。レーザーの振動数は 3.1eV、パルス長は 16fs にとっている。レーザーの強度が 7x1014W/cm2 に達した段階で、電子へのエネルギー移行が急に増加している。この強度において、外場の振動数と励起によって伝導帯に生じたキャリアの作るプラズマ振動数が一致し、共鳴的にエネルギー移行が進んだものと考えられる。この強度が光絶縁破壊に相当するものと考えられる。
（２）強レーザー場中で起こる分子のクーロン爆発（川下、矢花、中務）
強いレーザー場中にある分子は、多重イオン化の後にクーロン爆発を起こす。この過程に対し、時間依存密度汎関数理論によるシミュレーションを行った。イオン化の記述のため必要な広範な空間領域で計算を行うた
め、超並列計算機ＰＡＣS-CS を利用した。空間領域を３次元的に分割し、各々の領域の計算を異なるＣＰＵが担うことにより、効率の良い並列計算の実施が可能となっている。
Ｎ２分子を対象に計算を行った。まず、交換相関ポテンシャルの選択に関する検討を行った。中性原子や分子のイオン化に対しては、最適化有効ポテンシャルなどの自己相互作用補正を行ったポテンシャルを用いることの重要性が知られている。しかし一方で、最適化有効ポテンシャルでは全エネルギー計算が困難な場合があり、イオンに働く力を求める際に問題が生じる。局所密度近似は最も単純だが、中性からのイオン化率を過大評価することが知られている。これに対して、我々が検討を行ったところ、クーロン爆発が問題となる４価程度のイオン化では、局所密度近似と最適化有効ポテンシャルは、同一のレーザー強度でほぼ等しいイオン化価数となることが確認された。
分子のクーロン爆発実験で測定される量に、爆発後のイオンの運動エネルギー分布と運動量分布がある。
運動エネルギーは中性分子の平衡位置におけるクーロンエネルギーよりも小さいことが知られており、イオン化が進む過程でイオン間距離が増大することに起因すると考えられている。また、イオンの運動量がレーザーの偏極方向に偏ることが見いだされており、これは分子がレーザー場中でレーザー偏極方向にトルクを受けて回転することに起因すると指摘されていた。我々は、第一原理的なシミュレーションからこれらの解釈の妥当性について検討した。その結果、運動エネルギーの抑制に関しては、実験で見出されている 4 割程度の抑制が計算でも見られた。ただし、実験では非常に短いパルスを用いた場合には抑制がほとんど見られないことが報告されているが、計算ではパルスの長さが短い場合にも抑制が見いだされた。レーザー偏極方向への運動量の整列についても、計算によりそのような効果があることが確認できた。
３．研究業績
（１）研究論文
１．G. E. Ballentine, G.F. Bertsch, N. Onishi, K. Yabana,
Moment distributions of clusters and molecules in the adiabatic rotor model
Comp. Phys. Comm. 178, 48-51 (2008).
２．K. Nobusada, K. Yabana
Electric currents in ring-shaped molecules induced by circularly polarized laser pulses
Comp. Phys. Comm. 177, 54-54 (2007).
３．T. Otobe, K. Yabana,Density-functional calculation for the tunnel ionization rate of hydrocarbon molecules
Phys. Rev. A75, 062507 (2007).

４．T. Otobe, M. Yamagiwa, J.-I. Iwata, K. Yabana, T. Nakatsukasa, G.F. Bertsch,
First-principles electron dynamics simulation for optical breakdown of dielectrics under an
intense laser field
Phys. Rev. B77, 165104 (2008).
５．T. Otobe, K. Yabana, J.-I. Iwata,First-principle calculation for high harmonic generation in diamond
Submitted to Journal of Computational and Theoretical Nanoscience
６．K. Yabana, T. Nakatsukasa, M. Ito,
Time-dependent description for nuclear reaction dynamics in the continuum
Submitted to Proceedings of the 20th European Few-Body Conference.
７．矢花一浩多電子ダイナミクスの量子シミュレーション－時間依存密度汎関数理論の最近の発展－
日本物理学会誌 Vol.62, No.6, pp.406-414 (2007).
８．N.Hinohara, T.Nakatsukasa, M.Matsuo, and K.Matsuyanagi,
Gauge-Invariant Formulation of Adiabatic Self-Consistent Collective Coordinate Method
Prog. Theor. Phys. 117 (2007) 451-478.
９．M.Ito, K.Yabana, T.Nakatsukasa, and M.Ueda,Fusion Reaction of Halo Nuclei: A Real-Time Wave-Packet Method for Three-Body Tunneling Dynamics
Nucl. Phys. A787 (2007) 267c-274c.
１０．T.Nakatsukasa and K.Yabana,
Real-time Skyrme TDHF dynamics of giant resonances
Nucl. Phys. A788 (2007) 349c-354c.
１１．T.Nakatsukasa, T.Inakura, and K.Yabana
Finite Amplitude Method for the Solution of the Random-Phase Approximation
Phys. Rev. C 76 (2007) 024318 (9 pages).
１２．N.Hinohara, T.Nakatsukasa, M.Matsuo, and K.Matsuyanagi,
Microscopic derivation of collective Hamiltonian by means of the adiabatic self-consistentcollective coordinate method
Prog. Theor. Phys. 119 (2008) 59-101.
１３．T.Nakatsukasa, K.Yabana, and M.Ito,
Time-dependent approaches for reaction and response in unstable nuclei
Eur. Phys. J. Special Topics 156 (2008) 249-256.
１４．T.Nakatsukasa, S.Shinohara, K.Yabana, H.Ohta,
Stochastic approach to correlation beyond the mean field with Skyrme interaction

to be published in Proceedings of the international nuclear physics
conference (INPC2007).
１５．T.Nakatsukasa, T.Inakura, K.Yabana
Linear response calculations with the time-dependent Skyrme density functional
to be published in Proceedings of the international symposium on Physics of Unstable Nuclei(ISPUN07).
１６．Lu Guo, J.A.Maruhn, P.-G.Reinhard, and Y.HashimotoConservation properties in the time-dependent Hartree-Fock theory
Phys. Rev. C77 (2008), 041301(R).
（２）学会発表（A）招待講演
１．K. Yabana,
TDDFT Description for Laser-Induced Breakdown in Dielectrics,
Minerva Gentner Symposium on Time-Dependent Density Functional Theory, Eilat, Israle, 16-21, Dec,
2007.
２．K. Yabana
Time-dependent description for nuclear reaction dynamics in the continuum – Reaction and responseof exotic nuclei –20th European Conference on Few-Body Problems in Physics, Pisa, Italy, 10-14 September 2007.
３．K. Yabana
Configuration mixing calculation with Skyrme Hamiltonian
UNEDF workshop, Pac Forest, Washington, U.S.A., 12-17, August, 2007.
４．K. Yabana,
Simulation for electron dynamics in solid under intense laser pulse
ISSP International Symposium on “Foundations and Applications of the Density Functional Theory”,1-3, August, 2007.
５．K. Yabana,Time-dependent description for nuclear reaction dynamics in the continuum – Fusion mechanismof halo nuclei –International Nuclear Physics Conference, Tokyo, Japan, June 3-8, 2007.
６．T.Nakatsukasa
Time-dependent approaches for reactions and responses of unstable nuclei
Invited talk at Theoretical Nuclear Physics School ``Exotic Nuclei: New Challenges'',(Les Houches, France, May 7-18, 2007).

７．T.Nakatsukasa
Time-dependent method for structure and reaction
Invited talk at International Workshop on

Many-body open quantum systems: From atomic nuclei

to quantum dots (Trento, Italy, May 14-18, 2007).
８．T.Nakatsukasa
Dipole response calculations with the Skyrme density functional
Invited talk at UNEDF Workshop (Pack Forest, Washington, USA, August 13-17, 2007).
９．T.Nakatsukasa
Nuclear response function calculated with the time-dependent Skyrme density functional
Invited talk at 3rd Japanese-German Workshop on Nuclear structure and astrophysics(Chiemsee, Germany, September 29-October 2, 2007).
１０．T.Nakatsukasa, S.Shinohara, K.Yabana, H.Ohta
Stochastic approach to correlation beyond the mean field with Skyrme interaction
Talk at the international nuclear physics conference (INPC2007) (Tokyo, Japan, June 3-8, 2007).
１１．T.Nakatsukasa, T.Inakura, K.Yabana
Linear response calculations with the time-dependent Skyrme density functional Talk at theinternational symposium on Physics of Unstable Nuclei (ISPUN07 (Hoi An, Vietnam, July 3-7, 2007).
１２．T.Nakatsukasa
Time-dependent method for nuclear dipole response
Talk at 4-th JUSTIPEN-EFES Workshop on shell structure of exotic nuclei (RIKEN, Wako, Japan,June 22, 2007).
１３．Y. Hashimoto and T. Horibata
Tilted-axis rotation and wobbling motion in the framework of three-dimensional cranked HFBInternational Nuclear Physiscs Conference (INPC 2007) 東京国際展示場, 2007 年 6 月 3-8 日
１４．Y. Hashimoto and K. Nodeki
Small-amplitude responses in light nuclei in TDHFB calculation with Gogny force
First FIDIPRO-JSPS Workshop on Energy Density Functionals in Nuclei Keurusselka,(Jyvaskyla),Finland. 2007 年 10 月 25-27 日
１５．T. Inakura, T. Nakatsukasa and K. Yabana
Fully self-consistent calculation of isovector dipole response: Systematic study up to A=50
First FIDIPRO-JSPS Workshop on Energy Density Functionals in Nuclei
Keurusselka,(Jyvaskyla), Finland. 2007 年 10 月 25-27 日
１６．Y. Kawashita, T. Nakatsukasa, K. Yabana
TDDFT calculation for electron-ion dynamics in molecules under intense laser pulse
IPPS International Workshop and Symposium on"Foundations and Applications the Density

Functional Theory" 東大物性研 2007.8.1-3
１７．Shuichiro Ebata, Takashi Nakatsukasa, Tsunenori Inakura, Yukio Hshimoto, Kazuhiro Yabana
Calculation of Response Function with TDHF+BCS in Real-TimeCNS International Summer School [理研 仁科ホール,2007.08.28 - 09.01]（B）その他の学会発表
１．矢花一浩光に起因する電子・イオンダイナミクスに対する実時間シミュレーション法の開発「次世代量子シミュレータ・量子デザイン手法開発」特定領域研究19年度成果報告会、国際高等研究所、2008
年 3 月 18－19 日
２．矢花一浩超短パルスレーザーと物質の相互作用－TDDFT による計算－「凝縮系の超高速現象とコヒーレント制御」、東京工業大学国際交流会館多目的ホール、2008 年 2 月 21－22日
３．矢花一浩電子励起ダイナミクスの第一原理シミュレーション
第 3 回「計算科学による新たな知の発見・統合・創出」シンポジウム「PACS-CS システムと計算科学」、筑波大学
計算科学研究センター、2007 年 9 月 3－4 日
４．矢花一浩光に起因する電子・イオンダイナミクスに対する実時間シミュレーション法の開発
特定領域研究「次世代量子シミュレータ・量子デザイン手法の開発」平成 19 年度顔合わせ会、国際高等研究所、2007 年 6 月 8－9 日
５．矢花一浩、はじめに
日本物理学会第 63 回年次大会、理論核物理領域シンポジウム「原子核物理学における大規模計算の現在と未来」、近畿大学本部キャンパス、2008 年 3 月 23 日６．中務孝
原子核密度汎関数理論の発展
研究会「原子核物理學の将来像」、京都大学、2007 年 8 月 8-9 日
７．中務孝
核構造大規模計算と密度汎関数法による数値シミュレーション
日本物理学会年次大会シンポジウム「原子核物理學における大規模計算の現在と未来」、
近畿大学、2008 年 3 月 23-26 日
８．橋本幸男、堀端孝俊
GCM による 182Os の TAR モードの計算
日本物理学会 第 62 回年次大会 北海道大学 2007 年 9 月 21-24 日

９．橋本幸男、堀端孝俊
182Os における傾斜角回転に対応するバンド状態
日本物理学会 第 63 回年次大会 近畿大学 2008 年 3 月 22-26 日
１０．稲倉恒法, 中務孝, 矢花一浩
有限振幅法を用いた Skyrme-RPA 計算
日本物理学会 第 62 回年次大会 北海道大学 2007 年 9 月 21 日
１１．稲倉恒法, 中務孝, 矢花一浩
Systematic study of isovector dipole mode up to A=50KEK 研究会『原子核・ハドロン物理 横断研究会』 KEK 2007 年 11 月 19-21 日
１２．稲倉恒法, 中務孝, 矢花一浩
Systematic study of photoabsorption spectra in nuclei -Fully selfconsistent RPA calculaitonシンポジウム「未来の素粒子・原子核数値シミュレーション」
つくば国際会議場「エポカルつくば」 2007 年 12 月 19-20 日
１３．稲倉恒法, 中務孝, 矢花一浩
平均場理論に基づいた電気双極子励起の応答関数に関する系統的計算
2008 年国立天文台研究会 プロセス元素合成の統合的理解 ー 量子ビームで探る宇宙進化の理解を目指して ー筑波大学 2008 年 3 月 13-14 日
１４．稲倉恒法, 中務孝, 矢花一浩
Electric dipole mode の系統的計算－中重核までの Skyrme-RPA 計算
日本物理学会 第 63 回年次大会 近畿大学 2008 年 3 月 26 日
１５．川下洋輔 中務孝 矢花一浩
「分子における強レーザー場中の電子・イオンダイナミクスの第一原理計算」
第 8 回光量子科学研究シンポジウム 関西原研 2007.6.4-5
１６．川下洋輔 中務孝 矢花一浩
「分子の多重イオン化過程に対する第一原理シミュレーション」
第 1 回分子科学討論会２００７仙台 東北大 2007.9.17-20
１７．川下洋輔 中務孝 矢花一浩
「分子の多重イオン化過程に対する第一原理シミュレーション」
日本物理学会第 62 回年次大会 北海道大2007.9.21-24

１８．川下洋輔 中務孝 矢花一浩
「実時間・実空間法による電子・イオンダイナミクスの大規模並列計算」
次世代スーパーコンピューティングシンポジウム２００７ MY PLAZA ホール 2007.10.3-4

１９．川下洋輔 中務孝 矢花一浩
「実時間・実空間 TDDFT による電子・イオンダイナミクスの大規模並列計算」
スーパーコンピューターワークショップ２００８"大規模計算と分子のダイナミクス"

自然科学研究機構 岡崎コンファレンスセンター 2008.2.18-19
２０．川下洋輔 中務孝 矢花一浩
「TDDFT を用いた分子の多重イオン化過程に対する第一原理シミュレーション」
凝縮系の超高速現象とコヒーレント制御 東工大 2008.2.21-22
２１．川下洋輔 中務孝 矢花一浩
「実時間・実空間 TDDFT による電子・イオンダイナミクスの大規模並列計算」
次世代ナノ統合シミュレーションソフトウェアの研究開発 第 2 回公開シンポジウム
自然科学研究機構 岡崎コンファレンスセンター 2008.3.4-5
２２．川下洋輔 中務孝 矢花一浩
「強レーザーパルスによる分子の多重イオン化過程に対する第一原理シミュレーション」
第 55 回応用物理学関係連合講演会

日本大学理工学部船橋キャンパス 2008.03.27-30

２３．江幡修一郎、中務 孝、稲倉恒法、橋本幸男、矢花一浩
TDHF-BCS を用いた実時間応答関数の計算 21aYE-7
日本物理学会 第 62 回年次大会 北海道大学, 2007.09.21-24
２４．江幡修一郎、中務 孝、稲倉恒法、橋本幸男、矢花一浩
TDHF+BCS を用いた実時間応答関数の計算 26pZH-5
日本物理学会 第 63 回年次大会 近畿大学, 2008.03.22-26
２５．野出木健一、橋本幸男、中務 孝、矢花一浩
Gogny 相互作用を用いた時間依存 HFB の計算
日本物理学会 第 62 回年次大会 北海道大学, 2007.09.21-24
２６．野出木健一、橋本幸男
Gogny 力による TDHFB 計算
KEK 研究会『原子核・ハドロン物理 横断研究会』 KEK 2007 年 11 月 19-21 日
２７．野出木健一、橋本幸男、中務 孝、矢花一浩
Gogny 力を用いた時間依存 HFB 計算による軽い変形核の応答
日本物理学会 第 63 回年次大会 近畿大学, 2008.03.22-26
２８．松山裕右、丸山敏毅、矢花一浩
クォークの分子動力学法による原子核衝突のシミュレーション
日本物理学会第６３回年次大会、近畿大学、2008 年 3 月 22-26 日