Development of energy demand models for cooperative distributed EMS

It is a project adopted for Research Proposal Based Use in FY 2019.

Numerical Approach to the particle and nuclear physics

核物理研究センターでは、サブアトミック世界の強い相互作用の基本原理である量子色力学（ＱＣＤ）を中心に、素粒子・ハドロン・原子核物理学とその関連分野にわたる研究を、全国の共同利用研究者とともに行っている。スーパーコンピューターは、大規模計算を必要とする理論的な研究を中心に使用されている。その内容は、格子ＱＣＤによる数値シミュレーション、原子核構造及びハドロン物理現象論の数値的な解法を始め、深層学習の加速器実験への応用を含む。

(i) 格子QCD
ハドロン分光学（質量などの計算）とハドロン間相互作用、クォーク間相互作用、有限温度・密度における真空の性質・状態方程式、閉じ込め・カイラル対称性の問題

(ii) 少数多体系・原子核構造
３-４核子系の精密計算、原子核密度汎関数法

(iii) 機械学習
加速器実験への応用（加速器の運転制御、Belle実験・ILC実験での信号識別・フレーバー識別、ILC SiDカロリメータのエネルギー較正）

(iv) その他
２次元格子フェルミオン系、超新星爆発、弦理論

Transport properties in Nafion

ナフィオン中の水および電解質イオンの輸送係数を求めることを目的とする。水溶液電解質を含むナフィオンを、圧力場下および電場下に置いた場合の、水分子および電解質イオンの輸送係数を、速度相関関数から線形応答理論を用いることによって求めた。数種類の水溶液電解質を含むナフィオンについて、Darcy係数および電導係数を求めた。イオン種および含水率によるこれらの係数の変化を得ることができた。

Simulation of radical reaction throughout plasma-chemical-vapor deposition process

目的：銅表面上のグラフェンプラズマ気相成長メカニズムを解明する。
内容：メタン分子より発生するラジカル分子を予想し、銅基板との相互作用をシミュレーョンで調べた
結果：シミュレーョンにて CH ラジカルを銅基板に運動エネルギー10eVで衝突させた。 C-H結合分解は容易に起きず、さらに複雑な過程の存在が示唆された。

Reaction path search for the asymmetric hydrogenation catalyzed by rhodium(III) monohydride complex

当研究室でこれまでに開発してきた3つの塩素によって架橋されたRh(III)モノヒドリド錯体による不斉水素化反応について、新たに単離した単核錯体をモデルとして、錯体の立体の評価および反応経路探索を行った。

Electronic properties of low-dimensional nanoscale materials

本プロジェクトでは、各種異種物質ならびに外部電界と広義複合構造を形成した低次元ナノスケール物質に着目し、複合構造形成下における電子物性の変調の可能性を密度汎関数理論に基づく第一原理計算の手法と有効遮蔽媒質法を組み合わせることにより明らかにした。
　本年度は、立体形状を有する炭化水素分子、トリプチセンとテトラフェニルメタンに着目し、それらの重合構造体を考えることで、それぞれ２次元と３次元の共有結合炭化水素ネットワーク物質の物質設計を行った。これ等の重合構造体は極めて安定な炭化水素結晶系であることが明らかになった。さらに、sp2炭素領域の共有結合ネットワーク中の特異な配置から、いずれの構造もディラックバンドと平坦バンドの組からな特異な電子状態を価電子帯と伝導体に有する半導体となることを明らかにした。

Numerical Approach to the particle and nuclear physics

核物理研究センターでは、サブアトミック世界の強い相互作用の基本原理である量子色力学（ＱＣＤ）を中心に、素粒子・ハドロン・原子核物理学とその関連分野にわたる研究を、全国の共同利用研究者とともに行っている。スーパーコンピューターは、大規模計算を必要とする理論的な研究を中心に使用されている。その内容は、格子ＱＣＤによる数値シミュレーション、原子核構造及びハドロン物理現象論の数値的な解法を始め、深層学習の加速器実験への応用を含む。

(i) 格子QCD
ハドロン分光学（質量などの計算）とハドロン間相互作用、クォーク間相互作用、有限温度・密度における真空の性質・状態方程式、閉じ込め・カイラル対称性の問題

(ii) 少数多体系・原子核構造
３-４核子系の精密計算、原子核密度汎関数法

(iii) 機械学習
加速器実験への応用（加速器の運転制御、Belle実験・ILC実験での信号識別・フレーバー識別、ILC SiDカロリメータのエネルギー較正）

(iv) その他
２次元格子フェルミオン系、超新星爆発、弦理論

Numerical Simulation on Gas-Liquid Two-phase Flow in Large-scale Porous Media

大規模多孔質場の構造形成過程について，気液二相流挙動の解析が可能な格子ボルツマン法を適用し，ファンデルワールス状態方程式パラメータに対する計算安定性の評価を行った．その上で，大規模多孔質系における二相流挙動について検討した．

Numerical Approach to the particle and nuclear physics

核物理研究センターでは、サブアトミック世界の強い相互作用の基本原理である量子色力学（ＱＣＤ）を中心に、素粒子・ハドロン・原子核物理学とその関連分野にわたる研究を、全国の共同利用研究者とともに行っている。スーパーコンピューターは、大規模計算を必要とする理論的な研究を中心に使用されている。その内容は、格子ＱＣＤによる数値シミュレーション、原子核構造及びハドロン物理現象論の数値的な解法を始め、深層学習の加速器実験への応用を含む。

(i) 格子QCD
ハドロン分光学（質量などの計算）とハドロン間相互作用、クォーク間相互作用、有限温度・密度における真空の性質・状態方程式、閉じ込め・カイラル対称性の問題

(ii) 少数多体系・原子核構造
３-４核子系の精密計算、原子核密度汎関数法

(iii) 機械学習
加速器実験への応用（加速器の運転制御、Belle実験・ILC実験での信号識別・フレーバー識別、ILC SiDカロリメータのエネルギー較正）

(iv) その他
２次元格子フェルミオン系、超新星爆発、弦理論

Resurgence structure of the SU(3) gauge theory

It is a project adopted for Research Proposal Based Use in FY 2019.

Numerical Approach to the particle and nuclear physics

核物理研究センターでは、サブアトミック世界の強い相互作用の基本原理である量子色力学（ＱＣＤ）を中心に、素粒子・ハドロン・原子核物理学とその関連分野にわたる研究を、全国の共同利用研究者とともに行っている。スーパーコンピューターは、大規模計算を必要とする理論的な研究を中心に使用されている。その内容は、格子ＱＣＤによる数値シミュレーション、原子核構造及びハドロン物理現象論の数値的な解法を始め、深層学習の加速器実験への応用を含む。

(i) 格子QCD
ハドロン分光学（質量などの計算）とハドロン間相互作用、クォーク間相互作用、有限温度・密度における真空の性質・状態方程式、閉じ込め・カイラル対称性の問題

(ii) 少数多体系・原子核構造
３-４核子系の精密計算、原子核密度汎関数法

(iii) 機械学習
加速器実験への応用（加速器の運転制御、Belle実験・ILC実験での信号識別・フレーバー識別、ILC SiDカロリメータのエネルギー較正）

(iv) その他
２次元格子フェルミオン系、超新星爆発、弦理論

Stabilization of spatially dependent chiral superconductivity by magnetic field

金属/超伝導接合面から距離変化した超伝導では、自発磁化軸に平行な外部磁場は自発磁化状態を安定化させる。この磁場中の自発磁化状態の安定化により、常磁性電流とchirality変化が起こることを、Ginzburg-Landau理論により説明した。結果として、D4h対称性のEu、Eg既約での自発磁化状態が安定化する磁場による自由エネルギー減少は、常磁性電流と距離におけるone-nodalなchirality変化をもたらす。このone-nodalなchirality変化は、自発磁化状態の安定化のための磁場によるエネルギー減少が全距離において起こることで生じる。この機構について、秩序変数と自由エネルギーの距離依存性を数値計算から示して、更にone-nodalなchirality変化に対応した秩序変数の距離における符号変化を自由エネルギー磁場依存項に関係した近似式として導出することで説明した。

Theoretical Calculations on Cobalt Complexes Bearing a Sterically Hindered Phenanthroline Ligand

新規に合成を行った2,9-位に嵩高い置換基を有する0価のコバルト錯体の電子構造、および末端アルキンとの反応における反応機構に関する知見を得ることを目的に、Gaussian 09 を利用してMulliken の電子スピン計算と末端アルキンとの反応時のエネルギー計算を行った。Mullikenの電子スピン計算の結果、不対電子がコバルト中心上に局在化していることを確認した。また、末端アルキンのC-H結合がコバルト中心にσ配位した中間体の存在を理論的に明らかにした。

Efficient Implementation by Vector Computer for Graph Algorithm

2018年度の研究で問題になっていたバンクコンフリクトの確認を行なった．その結果，OpenMPをしないときにバンクコンフリクトが発生する状況が確認された．同じコードでOpenMP利用では発生していない．なお，この後，SX-Aurora環境で確認したところ，発生が見られず，原因究明はできていない状況である．

Numerical Simulations of Turbulent Mixing in Oceans

海面からの熱をエネルギー源とする台風は，海洋中に活発な乱流混合を引き起こす。この混合は海面水温を低下させるため，台風へのエネルギー供給も低下し，台風の勢力が弱まる。これまでの研究では，海洋中の混合を鉛直一次元過程で近似することがしばしば行われてきた。しかしながら，海洋中の3次元的な流れが，海面水温の分布や混合強度を変え，それが台風の発達に影響を与える可能性があることが示唆される。そこで，海洋中の3次元的な流れも再現し，1次元過程のみ場合と比較することで，海洋中の3次元過程が果たす役割を調べた。その結果，3次元過程は，海面水温の平均値を低下させるだけでなく，水へ分布の非対称性を著しく増加させることが明らかとなった。この非対称性は，台風の進路に影響することから，海洋の３次元過程が台風の進路に影響する可能性が示唆された。

湖沼～河川～海域における流動・水質・物質輸送の数値解析

三次元数値モデルSCHISM，FVCOMを用いて，感潮河川，淡水湖，汽水湖，閉鎖性内湾，内海～外洋域における流動・水質・物質輸送のシミュレーションを行った．

Numerical simulation of statistical structure in three-dimensional incompressible flow over cavities.

開いたキャビティを過ぎる三次元自励振動流の渦構造を解析することを目的とする．キャビティアスペクト比2.0について直接数値シミュレーションを実行し，自励振動と渦構造を可視化した．せん断層とキャビティ内の速度と圧力の各種統計量の分布構造を可視化し調べた．

Numerical Approach to the particle and nuclear physics

核物理研究センターでは、サブアトミック世界の強い相互作用の基本原理である量子色力学（ＱＣＤ）を中心に、素粒子・ハドロン・原子核物理学とその関連分野にわたる研究を、全国の共同利用研究者とともに行っている。スーパーコンピューターは、大規模計算を必要とする理論的な研究を中心に使用されている。その内容は、格子ＱＣＤによる数値シミュレーション、原子核構造及びハドロン物理現象論の数値的な解法を始め、深層学習の加速器実験への応用を含む。

(i) 格子QCD
ハドロン分光学（質量などの計算）とハドロン間相互作用、クォーク間相互作用、有限温度・密度における真空の性質・状態方程式、閉じ込め・カイラル対称性の問題

(ii) 少数多体系・原子核構造
３-４核子系の精密計算、原子核密度汎関数法

(iii) 機械学習
加速器実験への応用（加速器の運転制御、Belle実験・ILC実験での信号識別・フレーバー識別、ILC SiDカロリメータのエネルギー較正）

(iv) その他
２次元格子フェルミオン系、超新星爆発、弦理論

On the eplica exchange with solute tempering method applied to charged systems

レプリカ交換溶質加熱第二(REST2)法は分子動力学計算のサンプリング効率を上げ、高精度の計算を行うための計算手法である。本研究では電荷を持つ系のREST2法の挙動を調査することを目的とした。REST2法を実装したGROMACSをOCTOPUSに移植し、動作確認と検証をまで進めることができたが、当該年度はtool chain, job managerや認証系の問題等で研究が難航し、OCTOPUSでのプロダクトランに至っておらず、有意な結果を得ていない。

Controlled growth of single-crystalline metal nanowires via thermomigration across a nanoscale junction

流体においてよく知られている温度勾配下における質量輸送がナノスケールの金属材料においても発生することが観察されている。本研究はこの現象を利用したナノワイヤ生成手法を実験的、数値的に解析し、その制御手法を提案するものである。利用者は実験的に観察されているこの現象の詳細を明らかにするために、分子動力学法を用いて基材表面及びその直上にあるナノワイヤおける原子の拡散プロセスの解析を行い、結晶粒界における原子拡散が粒界面に垂直な引張荷重により促進されることを確認した。その結果、発表論文で行った実験及びFEM解析により推測された、引張荷重によるナノワイヤの成長の促進が裏付けられる結果を得ることができた。

Electronic state analyses of ionic liquid / electrode interfaces by ATR-FUV spectroscopy and quantum chemical calculation

It is a project adopted for Research Proposal Based Use in FY 2019.

Atomistic evaluation of temperature dependence of critical stress of Mg{10-12} twinning deformation

目的　六方最密構造Mgにおける{10-12}双晶変形はMgにおける塑性変形の主な担い手として重要であるがその変形の温度依存性評価は未だ明らかにされていない．本研究ではその臨界応力の温度依存性評価について原子論に基づいて明らかにする
内容　加速分子動力学法を用いて温度条件を変更してMg{10-12}双晶変形の発生する臨界応力を評価した．その結果を用いてオロワンの式に基づいて温度ごとの臨界応力を予測した．
結果 　得られた臨界応力の温度依存性は常温での実験結果と良い一致を示し，さらに低温における臨界応力の温度依存性の変化が顕著であることがわかった．

Design of a photonic crystal wavelength division multiplex device

帯域幅が従来よりも広い、新たなフォトニック結晶直行格子導波路(OLW)を提案し、共振器からOLWへの光取り出し効率を計算した。計算には3次元FDTD法を用い、導波路有無によるQ値変化から結合効率を算出した。計算の結果、従来の導波路よりも帯域幅が2倍となり、光取り出し効率も良好な値となることがわかった。

Computational Research on a Factor that Governs Light-Emission Efficiency of Fluorescent Glasses

It is a project adopted for Research Proposal Based Use in FY 2019.

Comprehensive understanding of mechanism of laser-produced magnetic field generation.

目的　実験結果とシミュレーション結果との比較をするために、大型計算機が必要なため利用する

　　　　 　内容　金属平板に高強度レーザーを照射し，膨張するプラズマをシミュレーション計算する．
　　　 　　用いるコードはFLASHと呼ばれる理想MHDシミュレーションコードである．
　　　　　　 これによって得られたパラメーターとモデルを組み合わせ，物理現象の解明に努めている

　　 　結果　現在検証段階ではあるが，計算結果を論文に掲載する予定

Numerical Approach to the particle and nuclear physics

核物理研究センターでは、サブアトミック世界の強い相互作用の基本原理である量子色力学（ＱＣＤ）を中心に、素粒子・ハドロン・原子核物理学とその関連分野にわたる研究を、全国の共同利用研究者とともに行っている。スーパーコンピューターは、大規模計算を必要とする理論的な研究を中心に使用されている。その内容は、格子ＱＣＤによる数値シミュレーション、原子核構造及びハドロン物理現象論の数値的な解法を始め、深層学習の加速器実験への応用を含む。

(i) 格子QCD
ハドロン分光学（質量などの計算）とハドロン間相互作用、クォーク間相互作用、有限温度・密度における真空の性質・状態方程式、閉じ込め・カイラル対称性の問題

(ii) 少数多体系・原子核構造
３-４核子系の精密計算、原子核密度汎関数法

(iii) 機械学習
加速器実験への応用（加速器の運転制御、Belle実験・ILC実験での信号識別・フレーバー識別、ILC SiDカロリメータのエネルギー較正）

(iv) その他
２次元格子フェルミオン系、超新星爆発、弦理論

Numerical study on a conical shaped target for Laser Fusion Rocket

円錐ターゲットの推進性能を評価することを目的とし、円錐ターゲットと球状ターゲットの2つの場合において核融合反応が起きた直後を仮定して高温の核融合プラズマが推進剤を加熱し膨張する過程を二次元輻射流体コードStar2Dで計算し、両者を比較した。
数値計算の結果、円錐ターゲットの円錐推進剤の81 %は加熱されない事が分かった。また、輻射流体計算の終状態となる時間の推進剤部分の運動エネルギーと内部エネルギーの合計は計算体系全体のエネルギーの0.33 %であり円錐推進剤はほとんど推力に寄与しない事が分かった。球状及び円錐ターゲットの力積はそれぞれ12 Ns、0.75 Nsであった。運動量変換効率はそれぞれ47 %，55 %であった。

Formation of plasma channel in a high-density plasma using 3-dimentional particle simulation with dynamic load balance technique

高速点火レーザー核融合における実スケール計算を可能とするため、ノードあたりの粒子数を動的に再割当てしてノード間負荷を均等にする手法を粒子コードに採用し、最大10倍程度の高速化に成功した。このコードを用い、大きさ数百ミクロンｘ数十ミクロンの計算領域において、空間分解能数ナノメートル、時間分解能数フェムト秒で、最大１０ピコ秒まで計算を行い、高強度レーザーが高密度プラズマ中を侵入する際形成する密度チャンネル形成の実験結果再現に成功した。

Ab initio approach to light nuclear systems based on quantum chromodynamics

Comprehensive study of generalized baryon-baryon interaction including strangeness is one of the important subject of nuclear physics. In order to obtain a complete set of isospin-base baryon interactions, we perform a large scale lattice QCD calculation with almost physical quark masses corresponding to $(m_{\pi}, m_{K})\approx(146,525)$~MeV and large volume $(La)^4=(96a)^4\approx(8.1~{\rm fm})^4$. A large number of Nambu-Bethe-Salpeter correlation functions from nucleon-nucleon to $\Xi\Xi$ are calculated simultaneously. For the moment, we focus on the strangeness S=-1 channels of the hyperon interactions by means of HAL QCD method. Strong repulsive forces are found in Sigma-Nucleon systems.

Study on Fundamental Phenomena of Detonation and Its Application

水素/空気予混合気に対する詳細化学反応モデルを使用して，2次元及び3次元の非定常圧縮性粘性解析を行った．今年度は粘性の影響を考慮した，高次精度スキームを使用した2次元，3次元デトネーションについて解析を実施し，粘性の影響を評価した．

Topological analysis of neural network mitigating catastrophic forgetting

単一のゴールのみを学習させる代わりに、複数のゴールを切り替えながら学習させた場合、ニューラルネットワークの破滅的忘却が軽減できることが論理生物学の分野で示唆されている。本研究では、この知見を実データに応用し、複数のゴールを切り替えながら学習させた場合、破滅的忘却が軽減できることを確認した上で、得られたニューラルネットワークの構造的特徴を明らかにする。

Effects of high intensity Laser Plasma Interaction on hydrodynamic simulations

目的　プラズマの密度スケールが長い領域に高強度レーザーを照射した際に顕著に起こるレーザープラズマ相互作用(LPI)を考慮した流体シミュレーションコードの開発を行うために、その影響のモデル化に取り組む。特に、電磁粒子(PIC)シミュレーションを用いて現象の数値解析を行い、データ収集、および現象の詳細な解析を行う。
内容　 PICシミュレーションを用いて、長いスケール長のプラズマ中　でのLPIで発生する誘導ラマン散乱(SRS)、誘導ブリリアン散乱(SBS)をからの反射光スペクトル、および高速電子スペクトル等のデータ収集を行った。
結果　SRSとSBSによる反射光スペクトル、および発生する高速電子スペクトルのデータを収集した。従来研究に比べ、様々なスケール長、レーザー強度に対する現象について、長時間のシミュレーションを行うことができた。これによって、キャビティの発生、absolute SRS, convective SRSなどの成長に伴う反射光や高速電子の発生に強い影響を及ぼすことが分かり、今後のモデル化を進める上で、非常に有効なデータを蓄積することができた。

Exploring the origin of dielectric relaxation spectrum for liquid water with molecular dynamics simulation

It is a project adopted for Research Proposal Based Use in FY 2019.

Correlation between viscous shear relaxation and local stresses in liquids

本研究の目的は，液体金属の粘度の局所構造依存性を解明することである．この問題を解決するために，液体の局所せん断応力時間相関関数の初期応力状態依存性を調べた．結果，液体の局所構造パラメタの一つとして局所圧力を計算し，その初期圧力状態とせん断応力緩和の間に強い相関が存在することを見出した．圧縮環境下にある原子はせん断に対して速く緩和し，引っ張り環境下にある原子は遅く緩和することがわかった．さらに，材料によらない普遍性をもつことが示唆された．

Monopole effects on quark confinement and chiral symmetry breaking

核物理研究センターでは、サブアトミック世界の強い相互作用の基本原理である量子色力学（ＱＣＤ）を中心に、素粒子・ハドロン・原子核物理学とその関連分野にわたる研究を、全国の共同利用研究者とともに行っている。スーパーコンピューターは、大規模計算を必要とする理論的な研究を中心に使用されている。その内容は、格子ＱＣＤによる数値シミュレーション、原子核構造及びハドロン物理現象論の数値的な解法を始め、深層学習の加速器実験への応用を含む。

(i) 格子QCD
ハドロン分光学（質量などの計算）とハドロン間相互作用、クォーク間相互作用、有限温度・密度における真空の性質・状態方程式、閉じ込め・カイラル対称性の問題

(ii) 少数多体系・原子核構造
３-４核子系の精密計算、原子核密度汎関数法

(iii) 機械学習
加速器実験への応用（加速器の運転制御、Belle実験・ILC実験での信号識別・フレーバー識別、ILC SiDカロリメータのエネルギー較正）

(iv) その他
２次元格子フェルミオン系、超新星爆発、弦理論

Materials exploration by using materials informatics (MI)

新規無機化合物の化学的安定性を見積もるためのモデルを構築するために，Materials Projectデータベースから15万件のデータをダウンロードし，深層学習により任意の組成の1原子あたりのエネルギーを予測するモデルを構築した．平均絶対誤差が0.19eVのモデルを構築できた．これにより新規化合物の化学的安定性をある程度予測することができるようになった．

Numerical studies of mysteries in thermos-fluid physics

スーパーコンピュータを利用して様々な問題の熱流体数値解析が
おこなわれているが、例えば、直管内で乱流遷移する位置と入口乱れ強さの関係を解明で
きる数値解析や理論は存在してこなかった。当方では、確率論的Navier-Stokes方程式とそ
の境界条件の理論を新たに提案し、それに基づいて、現象の解明を可能にしてきている。
また、それを用いて、高効率ロケットエンジンの性能検討を行ってきている。

Prediction of turbulent modulation and erosion induced by turbulent vortex cavitation

クエット乱流において，気液相変化による密度変化と，流れ場とキャビテーションの相互作用を考慮した圧縮性気液二相流の直接数値シミュレーションを実現した．そして，乱流中に発生する渦キャビテーションを観察し，乱流渦とキャビテーションの相互作用のメカニズムを明らかにした．さらに，乱流渦キャビテーションによる圧力波から壊食を予測することで，乱流渦に起因するキャビテーション壊食の特徴を調べた．その結果，乱流渦の直下とその周囲において，キャビテーション壊食が発生する可能性が高いことを示した．

Numerical simulation of flows including dense solid particles

高濃度に固体粒子を含む流れでは，粒子スケールで起こる固体粒子間の衝突や，粒子－流体間の直接的な流体力学相互作用が，気泡やクラスターなどのメゾスケール構造の形成を引き起こし，さらにこれらの自発構造がより大きいスケールでの対流運動を誘起するなど，典型的なマルチスケール構造を持ちます．実際の工業的な応用では，多くの場合二相間の熱・物質輸送を伴い，さらには粒子サイズの不均一性や，非球形性が問題となることが多く，高精度な予測を行うには多くの物理因子を考慮する必要があります．本研究では，これらの素過程となる種々の現象のモデル開発に取り組んでいます．

Development of Large-Scale Simulation and Visualization Program for Air Traffic Flow Toward Post-K Supercomputer

本研究は複数の相互作用する航空機モデルで構成されるマルチエージェントシステムをベースに、安全性、効率性、ロバスト性を備えた最適なフライトスケジュールの議論に利用しやすい大規模航空交通流シミュレーションツールの開発を目的としている。その開発段階では、大量の航空機軌道データを扱うためのビジュアルアナリティクスツールを開発して活用している。シミュレーションと実際の航空交通流データに対応して、データ同化ステップと最適化結果の有効性を直感的に検証できるインタラクティブな可視化プログラムを提案した。

Development of energy demand models for cooperative distributed EMS

本研究は国勢調査等の社会調査データから日本全国の世帯を乱数生成し、生成した模擬人口を用いた生活行動、エネルギー需要シミュレーションを行い、住宅のエネルギー需要を高い時空間解像度で推計するアプローチを確立した。この結果、世帯構成、気象条件、住宅・設備の集積状況が重要な因子となりエネルギー需要が形成されていることが明らかになった。

QCDにおけるグルーオン励起とアーベリアン・ドミナンスの研究

核物理研究センターでは、サブアトミック世界の強い相互作用の基本原理である量子色力学（ＱＣＤ）を中心に、素粒子・ハドロン・原子核物理学とその関連分野にわたる研究を、全国の共同利用研究者とともに行っている。スーパーコンピューターは、大規模計算を必要とする理論的な研究を中心に使用されている。その内容は、格子ＱＣＤによる数値シミュレーション、原子核構造及びハドロン物理現象論の数値的な解法を始め、深層学習の加速器実験への応用を含む。

(i) 格子QCD
ハドロン分光学（質量などの計算）とハドロン間相互作用、クォーク間相互作用、有限温度・密度における真空の性質・状態方程式、閉じ込め・カイラル対称性の問題

(ii) 少数多体系・原子核構造
３-４核子系の精密計算、原子核密度汎関数法

(iii) 機械学習
加速器実験への応用（加速器の運転制御、Belle実験・ILC実験での信号識別・フレーバー識別、ILC SiDカロリメータのエネルギー較正）

(iv) その他
２次元格子フェルミオン系、超新星爆発、弦理論

Large eddy simulation of two-dimensional turbulent jet using the structure function model.

構造関数モデルを用いたLarge Eddy Simulationで二次元乱流噴流の発達を調べることを目的とする．レイノルズ数6000, 10000の二次元乱流噴流をＬＥＳで計算し，統計量を実験値と比較し解析した．構造関数モデルは数値的に安定で，時間平均速度や乱れ強さは実験値と一致する結果が得られた．

Development of CO2-free combustor for high intensity combustion and NOx reduction with simultaneous function

It is a project adopted for Research Proposal Based Use in FY 2019.

Making inflow database for coupling simulation of LES and thermal radiation analysis

本研究では街区の温度場・熱放射場の非定常的な特性を把握するための熱放射・LES連成解析を実施するために，準周期境界により作成された流入変動風にメソ気象モデルWRFにより作成された温度プロファイルをドライバ領域で与えることで風速と温度場の流入条件を作成した。その結果，ドライバ領域によって得られた流入条件では，WRFにより作成された温度の鉛直プロファイルを維持しながら，微細な変動を含んだ風速・温度場を持つ流入条件が作成されていることを確認した。

First-principles molecular dynamics simulation of water interfaces using hybrid functions

It is a project adopted for Research Proposal Based Use in FY 2019.

Study on turbulent heat transfer mechanism inner high-speed lithium flow irradiated by ion beam

核融合炉材料試験を実現するための核融合中性子源では、液体リチウム（Li）ターゲットに重陽子ビームを照射することで高エネルギー・高フラックス中性子場を形成するが、ビーム照射に伴う熱を除去するためにLiは高速で流される。施設の安全性や効率化を考える上で、液体Liターゲット内部の熱輸送現象を詳細に把握することは重要な意味を持つ。そこで本研究では熱負荷を受けるLi噴流内の乱流熱輸送のメカニズムの解明を目的とし、これまでなされていなかった大規模LES計算を実施した。結果として、底面壁近傍における乱流熱拡散と渦構造との相関を捉えることができた。

Free-energy analysis of sorption of small molecule to polymer studied by atomistic molecular simulation combined with the theory of solutions in the energy representation

目的：ポリマー材料への低分子成分の収着性を全原子モデルで解析する
内容：分子動力学(MD)シミュレーションとエネルギー表示溶液理論を用い、セルロースアセテートに低分子成分(H2O、CO2、N2、CH4)が収着する際の自由エネルギーを計算した。セルロースアセテートはアセチル置換度の異なる3水準を対象とした。
結果：H2Oは他成分より収着自由エネルギー変化（ΔG）が大きく、極性分子が収着しやすいことが示された。CH4のΔGはおおよそ0であり、セルロースアセテートへ収着しにくいと考えられる。また、アセチル置換度が増加すると各成分においてΔGは減少しており、置換度2が最も低分子成分を吸収しやすいことが示された。

Fully self-consistent Relativistic Brueckner-Hartree-Fock theory for nuclear matter

核物理研究センターでは、サブアトミック世界の強い相互作用の基本原理である量子色力学（ＱＣＤ）を中心に、素粒子・ハドロン・原子核物理学とその関連分野にわたる研究を、全国の共同利用研究者とともに行っている。スーパーコンピューターは、大規模計算を必要とする理論的な研究を中心に使用されている。その内容は、格子ＱＣＤによる数値シミュレーション、原子核構造及びハドロン物理現象論の数値的な解法を始め、深層学習の加速器実験への応用を含む。

(i) 格子QCD
ハドロン分光学（質量などの計算）とハドロン間相互作用、クォーク間相互作用、有限温度・密度における真空の性質・状態方程式、閉じ込め・カイラル対称性の問題

(ii) 少数多体系・原子核構造
３-４核子系の精密計算、原子核密度汎関数法

(iii) 機械学習
加速器実験への応用（加速器の運転制御、Belle実験・ILC実験での信号識別・フレーバー識別、ILC SiDカロリメータのエネルギー較正）

(iv) その他
２次元格子フェルミオン系、超新星爆発、弦理論

study on dynamics of viscoelastic fluid in wall-bounded flow

水などのニュートン流体に比べ，界面活性剤溶液にみられる粘弾性流体の流動は複雑である．粘弾性流体では，せん断流中で形成されるミセル構造などによって生じる付加的な応力が乱れを誘起したり，抑制したりする．しかし，そのメカニズムは十分な理解がされていない．本研究では，直接数値解析によって複雑な粘弾性流体の流動現象を理解する．

Reaction path search for the asymmetric hydrogenation catalyzed by rhodium(III) monohydride complex

当研究室でこれまでに開発してきた3つの塩素によって架橋されたRh(III)モノヒドリド錯体による不斉水素化反応について、新たに単離した単核錯体をモデルとして、錯体の立体の評価および反応経路探索を行った。

Accelerating In-Transit Co-Processing for Scientific Simulations Using Region-Based Data-Driven Adaptive Compression

It is a project adopted for Research Proposal Based Use in FY 2019.

Two-dimensional radiation-hydrodynamics simulation for Laser-produced plasmas

これまでに開発した2次元輻射流体シミュレーションコードSTAR2dに対して、本年度は流体ソルバーの改良や、相対論電子ビーム輸送ルーティンの組み込みを行った。また、開発したコードを以下の解析を行った。
(a) レーザー核融合爆縮過程をシミュレートするため、２次元輻射流体コードSTAR2Dの開発を進めた。本年度は流体ソルバー空間精度向上と境界条件設定の改良を行うとともに、保存量保存性の改良を進めた。また、開発したコードを爆縮過程のシミュレーションに供して、解析を行った。

(b) Ｘ線顕微鏡などの応用を目的とした軟Ｘ線源開発として、テーブルトップのジュール級レーザーを平板金ターゲットに照射し、生成される金アブレーションプラズマからの軟X線生成過程をシミュレートした。2次元輻射流体コードSTAR2Dでレーザーアブレーション過程をシミュレートし、得られたアブレーションプラズマの時間・空間プロファイルと輻射放出・吸収係数テーブルを用いたX線光線追跡シミュレーションにより放出輻射特性を評価した。シミュレーション結果と先行実験悔過とは、軟X線放出時間・空間プロファイル、照射レーザーから軟X線への変換効率等においてよい一致を示した。これより、シミュレーションの妥当性が確認された。今後、レーザー照射条件・ターゲット厚さを変えたシミュレーションを進め、軟Ｘ線放射の高効率化を目指す。

Application of ReaxFF to Large-Scale Chemical Reaction Systems in Materials Science

ReaxFF1) は、結合の生成・開裂を記述可能な反応力場を搭載した反応分子動力学計算プログラムである。この反応力場は、量子化学計算データに基づいて決定されているため、高い精度で化学反応を記述できる一方で、量子化学計算手法と比較して、数オーダー速い計算速度を有している2)。この特徴から、燃焼反応3,4)、触媒反応5)、および高分子合成反応6)など様々な分野で応用されている。
　我々は、この特異なプログラムを活用して、材料科学分野における複雑な化学反応が関与する現象の理解を進めたいと考えている。そこで、まずは ReaxFF の基本的な使い方を知ると共に、計算精度と速い計算速度の確認を目的に様々な先行文献の追試を行うこととした。例えば、Strachan 等によって報告された3)、広く知られた高エネルギー物質である RDX ([cyclic-[CH2N(NO2)]3) の熱分解挙動に関する研究を再現計算した。各反応条件における二酸化炭素や水などの低分子生成量の時間挙動は、先行文献と良い一致を示した。他にも様々な分野での ReaxFF 応用事例の再現計算を行ったが、いずれも先行報告と類似の結果を得た。よって、基本的使用法を修得すると共に、ReaxFF の反応ポテンシャルの精密さ、および高速計算に対応していることを再確認した。
　今後は、ReaxFF の特徴を生かして、今まで適用が困難であった材料科学分野の大規模な化学反応系へ適用することにより、複雑な反応ダイナミクスが現象を支配している新たな領域での研究を開始する予定である。

参考文献
1) https://www.scm.com/product/reaxff
2) van Duin, A. C. T. et al., J. Phys. Chem. A 2001, 105, 9396.
3) Strachan, A. et al., J. Chem. Phys. 2005, 122, 054502.
4) Ashraf, C. et al., J. Chem. Phys. 2017, 121,1051.
5) Senftle, T. P. et al., NPJ Comput. Mater. 2016, 2 15011.
6) Vashisth, A. et al., J. Phys. Chem. A 2018, 122, 6633.

Novel Method of PPI Process by Multi-Scale MD and Rigid-Docking

It is a project adopted for Research Proposal Based Use in FY 2019.

Atomistic study on deformation behavior of mille-feuille structures

種々の金属系および高分子系ミルフィーユ構造のミルフィーユ構造の形成条件および変形特性を解析するために必要となる電子論的・原子論的シミュレーション手法の開発・整備を進めた．具体的には，(1)第一原理計算に基づく合金の熱力学的安定性評価，(2)分子動力学法に基づく変形素過程の抽出，(3)粗視化・加速分子動力学法に基づく高分子凝集体の状態探索，を行うための手法の整備・開発を実施した．

Numerical Approach to the particle and nuclear physics

核物理研究センターでは、サブアトミック世界の強い相互作用の基本原理である量子色力学（ＱＣＤ）を中心に、素粒子・ハドロン・原子核物理学とその関連分野にわたる研究を、全国の共同利用研究者とともに行っている。スーパーコンピューターは、大規模計算を必要とする理論的な研究を中心に使用されている。その内容は、格子ＱＣＤによる数値シミュレーション、原子核構造及びハドロン物理現象論の数値的な解法を始め、深層学習の加速器実験への応用を含む。

(i) 格子QCD
ハドロン分光学（質量などの計算）とハドロン間相互作用、クォーク間相互作用、有限温度・密度における真空の性質・状態方程式、閉じ込め・カイラル対称性の問題

(ii) 少数多体系・原子核構造
３-４核子系の精密計算、原子核密度汎関数法

(iii) 機械学習
加速器実験への応用（加速器の運転制御、Belle実験・ILC実験での信号識別・フレーバー識別、ILC SiDカロリメータのエネルギー較正）

(iv) その他
２次元格子フェルミオン系、超新星爆発、弦理論

Influences of Nanoparticle on Ionics Fluids

陽陰イオンから構成される高温液体である溶融塩にナノスケールの微粒子を混合した際の熱物性の変化を調査することを目的とした．分子動力学法を用いて熱伝導率に与える影響を調査した．また，固体表面上でのイオン液体の濡れ性に関する調査を行った．数値解析にはOCTＯPUSを用いた．

Molecular dynamics simulation for thermal engineering at microscale

ナノ・マイクロメートルスケールのエネルギー輸送現象を原理的に理解して制御することを目的として，ナノ構造が凝縮過程や凝縮時の熱抵抗に与える影響を，分子動力学法を用いて調査した．また壁面での凝固現象が微粒子に及ぼす影響に関しても分子動力学解析を用いて現象解明を行った．
そのために大規模可視化対応PCクラスタを用いた．

Numerical simulation of self-sustained oscillations in two-dimensional jet – edge system.

噴流がエッジに衝突すると噴流が上下に発振する．噴流エッジ間距離を変え振動の変化を調べることを目的とした．噴流－エッジ間距離5から12まで変化させた計算を行い，流れ場や渦構造を可視化した．エッジ先端近傍の速度波形が変化し、渦の数が変化することによって振動モードの変化が起きることを明らかにした．

Numerical Approach to the particle and nuclear physics

核物理研究センターでは、サブアトミック世界の強い相互作用の基本原理である量子色力学（ＱＣＤ）を中心に、素粒子・ハドロン・原子核物理学とその関連分野にわたる研究を、全国の共同利用研究者とともに行っている。スーパーコンピューターは、大規模計算を必要とする理論的な研究を中心に使用されている。その内容は、格子ＱＣＤによる数値シミュレーション、原子核構造及びハドロン物理現象論の数値的な解法を始め、深層学習の加速器実験への応用を含む。

(i) 格子QCD
ハドロン分光学（質量などの計算）とハドロン間相互作用、クォーク間相互作用、有限温度・密度における真空の性質・状態方程式、閉じ込め・カイラル対称性の問題

(ii) 少数多体系・原子核構造
３-４核子系の精密計算、原子核密度汎関数法

(iii) 機械学習
加速器実験への応用（加速器の運転制御、Belle実験・ILC実験での信号識別・フレーバー識別、ILC SiDカロリメータのエネルギー較正）

(iv) その他
２次元格子フェルミオン系、超新星爆発、弦理論

Simulation of LEMP Propagation in the Earth-Ionosphere Waveguide Using the FDTD Method in the 2-D Spherical Coordinate System

二次元円筒座標系FDTD法(Finite Difference Time Domain)は二次元空間で三次元的に電磁界解析を行うことができるため，メモリの節約が可能になる。しかし，セルの形状が長方形のため地球の曲面を無視できる範囲での電磁界解析にしか適さない。一方で，二次元球座標系FDTD法は地球の曲面を考慮できるため，300 kmを超えて伝播する雷電磁波の計算にも有効である。本研究では，二次元球座標系FDTD法に必要な電磁界更新式の定式化を行い，垂直電界の実測波形と比較することで二次元球座標系FDTD法の適用可能性の検討を行った。

Accessing the Accuracy of Density Functional Theory through Structure and Dynamics of the Water–Air Interface

シュレディンガー方程式を解いて原子に働く力を計算する密度汎関数分子動力学法（density functional theory molecular dynamics, DFTMD）を用い、水/空気界面における水分子の性質を様々な汎関数を用いて比較した。試した計算手法の中では、revPBE0-D3(0)が水の物性をよく記述することがわかった。

Optical spectrum in the stipe-type pattern charge ordering of two-dimensional electron systems

2次元電子系においては、電子が規則的に並んだ電荷秩序が発生するが、様々な電荷秩序の光学励起状態に
予想される多電子励起性が興味深い。これまで電荷自由度に注目したスピンレスフェルミオンモデルを
6×6の格子で扱い、光学伝導度、および､各励起エネルギーにおける多電子励起性を求めてきたが、本年度は特に
ストライプ型電荷秩序を対象に解析を行った。スペクトル（σ）と電子励起数（Nex）を各振動数ごとに
求めたところ、スペクトル形状は比較的鋭いが、Nexを見ると、励起数は最大で６程度と比較的大きい。
これは１電子励起と多電子励起が混在していると大まかに理解される。

Consideration of numerical PIV measurement on turbulent premixed flames using DNS

乱流予混合火炎において重要な特性のうち，火炎変位速度を求めるために必要な局所流速の計測精度を評価するため，DNS（直接数値シミュレーション）データを用いて現実のPIV（粒子画像流速測定法）により近い数値計測手法（数値PIV計測）を構築する．
未燃混合気と既燃生成物の密度比を変化させた3次元DNSによって構築された乱流予混合火炎のデータに対して数値PIV計測を行う際に，トレーサ粒子の粒径が計測された局所流速に及ぼす影響について，DNSによって算出される局所流速と数値PIV計測によって計測された局所流速の値を比較して検討した．
本研究の条件下においては，未燃混合気と既燃生成物の密度比が高い場合，トレーサ粒子の粒径が大きい方が，DNSと数値PIV計測の局所流速の間で全体的に高い相関が見られた．

Study of energy-momentum tensor for Nf=2+1 full QCD using gradient flow

It is a project adopted for Research Proposal Based Use in FY 2019.

Large scale atomistic simulation of crack propagation in amorphous metals

大規模分子動力学シミュレーションにより金属アモルファスを作成し、作成した原子モデルに対して同シミュレーション手法を用いてせん断負荷を印加することにより、金属アモルファス材料中のき裂成長過程を解析した。これまでより大規模なモデル（650nm x 82nm x 4nm) を用いることにより、切り欠き部から原子ひずみの高い領域が帯状に形成され、せん断き裂が成長する過程をより長時間観察することができた。これによりこれまで明らかになっていないせん断帯の幅やせん断部での温度上昇に関する詳細な知見が得られた。

Conformational sampling of an mSin3 complex using multidimensional virtual-system coupled canonical MD

It is a project adopted for Research Proposal Based Use in FY 2019.

LES study of turbulent flow fields over hilly terrains ? Comparisons of inflow turbulence generation methods and SGS models

Large eddy simulation (LES) is performed to numerically simulate the turbulent flows over hilly terrains with smooth surface. Two types of inflow turbulence generation methods, wind tunnel replication method and consistent discrete random flow generation method, are used to simulate the approaching flows for LES study of turbulent flow fields over simplified topography. The advantages and limitations of these methods are discussed in the context of flow patterns and turbulence characteristics. The performance of three different subgrid-scale models are investigated: standard Smagorinsky model (SM), (b) coherent structure Smagorinsky model (CSM) and dynamic Lagrangian Smagorinsky model (DLSM). The numerical results obtained from the CSM and DLSM are in relatively poor agreement with those of the experiment. The size of the recirculation region predicted by the CSM and DLSM are smaller than that of the SM. It is found that the closed-wake region is formed behind the 2-D ridge, while the open-wake region is established on the leeward side of 3-D hill due to the existence of the secondary flow. The open-wake region is well reproduced by the SM rather than the CSM and DLSM. The non-Gaussian characteristics of wake flow and the topography-driven twist effects are also examined.

Computational fluid dynamics simulation of human left atrium based on 4D-CT images

左心房内血栓症の発症メカニズムに対する流体力学的観点からの検討にむけ，患者の4D-CT画像計測に基づく左心房内血流の数値流体力学シミュレーションを行った．左心房の壁面運動が，左心房内部の血流の攪拌，ひいては内部血流の滞留を抑制する役割を持つことが示唆された．

Exploring an approach using persistent homology for the confinement phase transition of quark-gluon matter

クォーク・グルーオン物質の非閉じ込め相転移に現れる物理現象を解明するため、パーシステントホモロジー解析を用いた有効なアプローチを模索している。今回はその一環として、QCDの閉じ込め機構を説明するとされているセンターボルテックスの時空間配置がどのように温度変化するかを調べた。数値シミュレーションによってセンターボルテックス構造を計算し、パーシステントホモロジー解析を行ってその「穴」を表す位相幾何学的特徴（1次のホモロジー）を調べた結果、転移点において閉じ込めからsQGPの状態への遷移を説明する顕著な変化が確認された。

Measuring Φ100 mm aspheric mirror by Non-Contact Three Dimensional Nano-profiler Using Normal Vector Tracing Method

我々は、非接触かつ自由曲面の測定が可能な形状測定法として、法線ベクトル追跡型非接触ナノ形状測定法の開発に取り組んでいる。本研究では、実際にΦ100 mm凹非球面ミラーの形状測定を行った。その結果、試料作製時の切削痕由来と思われる同心円状の形状誤差があることが分かった。

Development of transition-metal complex-mediated C-F bond cleavage reaction

炭素(sp3)－フッ素結合は非常に強固な化学結合の一つであり、その切断反応は困難な反応の一つとして知られている。当研究室では、電子豊富なニッケル錯体を用いることで、炭素(sp3)－フッ素結合切断が進行することを見出した。本研究では、このニッケル錯体による炭素(sp3)－フッ素結合切断過程を理論計算ソフトであるGaussian 09を使用した密度汎関数理論（DFT）計算によって解析した。その結果、ベンゼン環への電子供与を経由して炭素(sp3)－フッ素結合が切断されるといった、これまで提唱されたことのない反応機構によって本反応が進行していることが示唆された。

Simulation of polarization in ionosphere plasma

スーパーコンピューターを用いた2 流体シミュレーションを実施し、電離層プラズマ中の分極過程を評価した。計算機で与えた条件の下で、時間とともに分極が大きくなっている事が示された。今年度は、空間的、時間的に広大な範囲におけるプラズマ模擬実験のための、方法の構築をさらに進めた。

Lattice simulations in Z3-symmetric QCD

Z3対称な格子QCDの計算とそれに関連する計算を行った。
特に有限温度・有限密度領域における計算を行った。
位相クエンチ近似を使用しているので、現在の段階の結果はアイソスピン化学ポテンシャル存在する場合と同等の状態に対応する。
温度が一定のとき、化学ポテンシャルが増加するとポリヤコフループの絶対値は増加する傾向が見られた。

Numerical simulations of the spin-orbital glass transitions in a pyrochlore magnet

不純物などの外的乱れなしで起こるスピングラス転移が、Y2Mo2O7というパイロクロア酸化物において30年ほど前から実験的に知られており、興味が持たれてきた。そのメカニズムは、理論的には未解明であった。我々はこの系におけるガラス転移に関する微視的理論模型を導出し、その性質を大規模モンテカルロシミュレーションによって明らかにした。具体的には、OCTOPUS上で、MPIを用いたノード内(48コア)並列計算を行った。用いた手法はレプリカ交換モンテカルロ法で、プロセス毎に温度の異なるレプリカのシミュレーションを行いつつ、プロセス間で温度の交換を行うというものである。本研究の結果、Y2Mo2O7においてスピンと軌道自由度が同時にガラス転移を示すことを発見した。

Numerical Approach to the particle and nuclear physics

核物理研究センターでは、サブアトミック世界の強い相互作用の基本原理である量子色力学（ＱＣＤ）を中心に、素粒子・ハドロン・原子核物理学とその関連分野にわたる研究を、全国の共同利用研究者とともに行っている。スーパーコンピューターは、大規模計算を必要とする理論的な研究を中心に使用されている。その内容は、格子ＱＣＤによる数値シミュレーション、原子核構造及びハドロン物理現象論の数値的な解法を始め、深層学習の加速器実験への応用を含む。

(i) 格子QCD
ハドロン分光学（質量などの計算）とハドロン間相互作用、クォーク間相互作用、有限温度・密度における真空の性質・状態方程式、閉じ込め・カイラル対称性の問題

(ii) 少数多体系・原子核構造
３-４核子系の精密計算、原子核密度汎関数法

(iii) 機械学習
加速器実験への応用（加速器の運転制御、Belle実験・ILC実験での信号識別・フレーバー識別、ILC SiDカロリメータのエネルギー較正）

(iv) その他
２次元格子フェルミオン系、超新星爆発、弦理論

Numerical Approach to the particle and nuclear physics

核物理研究センターでは、サブアトミック世界の強い相互作用の基本原理である量子色力学（ＱＣＤ）を中心に、素粒子・ハドロン・原子核物理学とその関連分野にわたる研究を、全国の共同利用研究者とともに行っている。スーパーコンピューターは、大規模計算を必要とする理論的な研究を中心に使用されている。その内容は、格子ＱＣＤによる数値シミュレーション、原子核構造及びハドロン物理現象論の数値的な解法を始め、深層学習の加速器実験への応用を含む。

(i) 格子QCD
ハドロン分光学（質量などの計算）とハドロン間相互作用、クォーク間相互作用、有限温度・密度における真空の性質・状態方程式、閉じ込め・カイラル対称性の問題

(ii) 少数多体系・原子核構造
３-４核子系の精密計算、原子核密度汎関数法

(iii) 機械学習
加速器実験への応用（加速器の運転制御、Belle実験・ILC実験での信号識別・フレーバー識別、ILC SiDカロリメータのエネルギー較正）

(iv) その他
２次元格子フェルミオン系、超新星爆発、弦理論

Simulation study of action potential propagation in the heart

心臓不整脈の実体は、心臓組織上での旋回性興奮の発生にあると考えられている。遺伝性QT延長症候群や、薬物の副作用としての薬物誘発性QT延長症候群における致死性不整脈の発症は、心筋細胞レベルでの再分極異常（早期後脱分極：early afterdepolarization）の発生が引き金となって、旋回性興奮（リエントリー）の発生によるものと考えられている。しかしながら、そのリエントリー生成メカニズムは未だ明らかでない。
　心筋細胞の活動電位（細胞）モデルユニットを電気的に結合した6 cm × 6 cmの心室組織（シート）モデルに、再分極異常細胞ユニット（EADユニット）をクラスター状に配置し、興奮伝播シミュレーションを様々な条件下で繰り返し実施した（図）。EADクラスターのサイズ、シート内での位置、複数のEADユニットクラスターが存在した場合のクラスター間の相互作用の強度といったパラメータが、リエントリー生成に与える影響を検討した。
EADが発生することでEADユニットの興奮の終了（活動電位の再分極）が遅れる。EADユニットクラスター以外の非EADユニット（no-EAD units）は、既に興奮が終了しており、EADユニットとの間に大きな電位差が生まれる。この電位差によって、電流が発生し、すでに興奮が終了した非EADユニットを再び興奮させる。つまり、リエントリー（興奮の再入）を引き起こす。今回、EADユニットクラスターのサイズ（ユニット数）が心室組織全体の約20%以上を占め、かつ、そのEADユニットクラスターが組織辺縁部に近い場所にあるときに、旋回性興奮波の発生を見た。今後は、EADの発生からリエントリー誘発（致死性不整脈発症）までの条件を更に詳細に検討する予定である。

Dissimilar Transport of Momentum and Heat Transfer in Pulsating Turbulent Flows: Effect of Themal Boundary Condition

エンジン筒内の熱損失を高精度に予測する乱流伝熱モデルの開発を目標に，その基礎的な乱流場である脈動を伴う乱流熱伝達場を対象として，運動量輸送と熱輸送の非相似性が生じる機構について検討を行った．2019年度は熱的境界条件が非相似性に及ぼす影響を検討した結果、脈動乱流における非相似性は熱的境界条件が運動量輸送と相似な境界条件においても生じることを示した。

Numerical Simulations of Multiphase Flows

VOF法/MTHINC法，Boundary Data Immersion法に基づく気液二相，流体構造連成コードの整備を進め，MPI並列版の動作確認を行なった．OCTOPUSを使用してのプロダクトランは行なわなかった．

Numerical Approach to the particle and nuclear physics

核物理研究センターでは、サブアトミック世界の強い相互作用の基本原理である量子色力学（ＱＣＤ）を中心に、素粒子・ハドロン・原子核物理学とその関連分野にわたる研究を、全国の共同利用研究者とともに行っている。スーパーコンピューターは、大規模計算を必要とする理論的な研究を中心に使用されている。その内容は、格子ＱＣＤによる数値シミュレーション、原子核構造及びハドロン物理現象論の数値的な解法を始め、深層学習の加速器実験への応用を含む。

(i) 格子QCD
ハドロン分光学（質量などの計算）とハドロン間相互作用、クォーク間相互作用、有限温度・密度における真空の性質・状態方程式、閉じ込め・カイラル対称性の問題

(ii) 少数多体系・原子核構造
３-４核子系の精密計算、原子核密度汎関数法

(iii) 機械学習
加速器実験への応用（加速器の運転制御、Belle実験・ILC実験での信号識別・フレーバー識別、ILC SiDカロリメータのエネルギー較正）

(iv) その他
２次元格子フェルミオン系、超新星爆発、弦理論

Interglueball potential in SU(2) lattice gauge theory

SU(N) Yang-Mills理論におけるグルーボールは暗黒物質の良い候補であり、宇宙の中規模構造の形成において重要である暗黒物質間の散乱を研究する必要性がある。
本研究では格子ゲージ理論におけるHALQCD法を用いてSU(N) Yang-Mills理論におけるグルーボール間相互作用を計算し、近距離における斥力的な振る舞いを示した [1,2,3]。
また、グルーボール間の散乱断面積を計算し、銀河の衝突の観測などからSU(2) Yang-Mills理論のスケールパラメータに対する制限を導出した。

Helicity instability effects of supersonic Batchelor vortices

It is a project adopted for Research Proposal Based Use in FY 2019.

Numerical Approach to the particle and nuclear physics

核物理研究センターでは、サブアトミック世界の強い相互作用の基本原理である量子色力学（ＱＣＤ）を中心に、素粒子・ハドロン・原子核物理学とその関連分野にわたる研究を、全国の共同利用研究者とともに行っている。スーパーコンピューターは、大規模計算を必要とする理論的な研究を中心に使用されている。その内容は、格子ＱＣＤによる数値シミュレーション、原子核構造及びハドロン物理現象論の数値的な解法を始め、深層学習の加速器実験への応用を含む。

(i) 格子QCD
ハドロン分光学（質量などの計算）とハドロン間相互作用、クォーク間相互作用、有限温度・密度における真空の性質・状態方程式、閉じ込め・カイラル対称性の問題

(ii) 少数多体系・原子核構造
３-４核子系の精密計算、原子核密度汎関数法

(iii) 機械学習
加速器実験への応用（加速器の運転制御、Belle実験・ILC実験での信号識別・フレーバー識別、ILC SiDカロリメータのエネルギー較正）

(iv) その他
２次元格子フェルミオン系、超新星爆発、弦理論

The Development of Turbulent Mixing Parameterization Scheme in the Surface Ocean under the Stratification and Earth’s Rotation

世界規模の海洋循環を再現する海洋大循環モデルでは大気海洋相互作用などに重要な混合層深度の再現性が不十分である。海洋大循環モデルでは混合層深度の再現するために乱流混合をパラメタリゼーションで表現しているが、このパラメタリゼーションには海洋で大きな役割を果たす密度成層の影響が十分に考慮されていないため、乱流混合を陽に再現可能な数値モデルを用いて乱流混合に対して密度成層が与える影響を明らかにし、その効果を定式化することで新たなパラメタリゼーションを開発することを目的とした。数値実験を行った結果、密度成層が強い深さにおいて従来のパラメタリゼーションから大きくずれることが明らかとなった。従来のパラメタリゼーションに成層の効果を導入することで、表層流速の鉛直分布の再現性が向上した。

Simulation study of cardiac electrical excitation using virtual heart model

スーパーコンピュータ上に仮想心臓モデルを構築し，心臓電気生理学的シミュレーションを行うことで致死性不整脈のメカニズム解明や，予防・診断に役立てることを目的している．内容：心筋細胞の電気興奮に伴う電位変化（活動電位）を再現することが可能なユニット約2,000万個を組み合わせて心室形状モデルを構築した．モデルには，心外膜から心内膜にかけての電気生理学的性質の不均一性を組み込んだ．本研究では，遺伝性の心臓疾患の一つであるブルガダ症候群を想定し，心筋組織内の興奮伝導障害の領域を設定し，伝導障害と心室性不整脈の誘発性および持続性との関係について検討した．心臓電気生理学的検査を想定し，右室流出路の心内膜側を電気的に刺激し，不整脈が誘発されるときの応答をシミュレートした．伝導障害領域の大きさや，領域内におけるユニット間の電気的結合の強さなどの伝導障害の程度や，電気刺激の条件を変えながらシミュレーション実験を行った．伝導障害の領域が直径3 cm，ユニット間の電気的結合を正常の10%まで減少させると不整脈が誘発される場合があったが，電気刺激のタイミングが同一であっても，不整脈の誘発性は電気刺激の部位に大きく影響された．

Formation of plasma channel in a high-density plasma using 3-dimentional particle simulation with dynamic load balance technique

It is a project adopted for Research Proposal Based Use in FY 2019.

Large scale atomistic simulation of defects and hydrogen behavior in iron

大規模分子動力学シミュレーションによりbcc鉄に点欠陥（vacancyおよびinterstitial）を導入した原子モデルを作成し、作成した原子モデルに対して水素原子を導入した場合と導入しない場合とでこれらの点欠陥の拡散挙動がどのように変化するかを解析した。その結果　水素の有無にかかわらずinterstitialの拡散はvacancyに比べて速いことが確かめられた。また、水素の有無によるこれらの拡散速度の変化については、温度や水素濃度によって水素の影響が変化することがわかったが、その変化のばらつきが大きく、系統的な結果を得るためには今後さらに慎重な解析が必要であることがわかった。

Formation of planar structures with rollable and jumpable cubic modular robots

This research focuses on the control algorithm of an innovative modular robot, which is a cube-shaped robot capable of performing rolling and jumping actions in multiple directions, and coupling with the other robots to form closely-packed structures. The configuration of the robot referenced to the actuator mechanisms developed in our previous research for allowing rolling and jumping actions to be performed by cubic robots in multiple directions, and the bonding mechanism included in a robot called M-block in previous research for allowing robots to form closely-packed structures of different 3D shapes. Because of the difficulties in building a swarm of new robots and conducting experiment under micro gravity environment, a computer simulation was conducted to investigate the interaction of robots with the proposed control algorithm under micro gravity environment. Due to the numerous magnets involved in the interaction between robots, massive computation on magnetic torque and force are handled in parallel using large-scaled computer systems OCTOPUS. The result has confirmed the effectiveness of the proposed algorithm in the construction of a planar structure under micro gravity environment.

Radiation-hydrodynamic numerical simulation for laser processing, material processing, and laser fusion

多次元輻射流体数値シミュレーションを駆使し、レーザー駆動超高圧を利用したレーザー加工プロセス、新物質創生、及びレーザー核融合エネルギー創生に関わるレーザープラズマ相互作用の学理を深めることを目的としている。本年度は、レーザー核融合研究として、大阪大学レーザー科学研究所で実施したGXII/LFEX実験の対向ビーム爆縮シミュレーションを行い、爆縮密度温度等、物理パラメータを解析した。得られた結果の概要は以下のとおりである。

(a) 二次元輻射流体シミュレーションコードSTAR2Dを用いて、実験に対応したレーザーおよびターゲット条件で爆縮シミュレーションを行った。レーザー条件は２種類の配位、ターゲット条件は２種類である。計算の結果、対向ビーム爆縮では非球対称な楕円体状の爆縮コアが形成されることが見出された。
(b) STAR2Dコードによる爆縮シミュレーションの結果をさらにポストプロセス処理し、爆縮コアからのX線発光を実験で使用したX線ストリークカメラに対応した感度分布を用いて計算機上で再現し、実験計測結果と比較した。その結果、概ね実験との良好な一致がみられたものの、爆縮タイミングについて、実験よりも計算が100ps程度遅れていることが見出された。この実験結果とシミュレーション結果の差の原因としてはレーザー吸収のモデリングにおける不正確さなどが考えられ、今後のコード開発及びプラズマモデリングの方向性が示唆された。

Visualization of technical solution by self-organizing maps for accelerating innovation

特許文献データを用いて技術課題の解決手段の可視化を行い，新規アイデアの着想を支援する．単語レベルの共起確率ベクトル・相関係数ベクトル、そしてマルチビューを適用し自己組織化マップを構築することで，技術要素の直接的な関連性及び潜在的な関連性の視覚化を実現した．

Numerical Study on Compressible Viscous Flow including Supersonic Combustion

マッハ数3程度の超音速空気流中に音速の水素が噴射した場合の衝撃波/境界層干渉と水素拡散の挙動を明らかにするために，化学種の質量保存を含む3次元非定常圧縮性粘性解析を行った．大規模な並列を行うことで，効率的に計算結果を得ることが可能となった．そして，噴流の上流側のせん断層や噴流背後の後流域において非定常な渦が形成され，空気と混合していることが分かった．

Development of a comparative method for time-course single-cell data

遺伝子の発現量などから細胞状態の遷移過程を捉える疑似時間解析において、実験条件の異なるデータから得られる疑似時間経路を理論的に比較する計算手法が求められている。既存手法では線形構造の比較しか行えないため、経路に分岐が含まれる場合、その選択に解析者によるバイアスがかかっていた。本研究では、比較したい疑似時間経路に分岐が含まれるような場合に、客観的に比較可能な手法CAPITALを開発した。公開されているヒトの筋細胞データの組を用いたところ、マーカー遺伝子の発現ダイナミクスを正しく捉え、既報の結果と概ね合致したため、開発手法の妥当性がある程度示されたと考える。

Numerical Simulations of Wall-Turbulence Transition and Multi-phase Flows of Newtonian and Viscoelastic Fluids

壁面せん断流の亜臨界乱流遷移において，その遷移過程と発現構造は複雑で，理論的アプローチが困難である．本研究では，ニュートン流体流れの後退翼平板境界層やGiesekus粘弾性モデル流体流れのバックステップ流などの直接数値シミュレーションを行い，低レイノルズ数下での粘弾性による流動構造の変化や（ニュートン流体を含む）乱流遷移過程の変化を調査した．今年度は他に，テイラークエット流や，分散性粒子を含む固液二相流，OpenFOAMを用いた気液界面を含むミリスケール液滴・液膜の解析などを行った．

Hydrodynamic boundary condition at the atomic scale

分子動力学法で作り出したサブナノスケール流れが、固液摩擦係数の周波数依存性つまり粘弾性を考慮すればマクロな Navier 境界条件と Stokes 方程式で記述できるということを初めて示した．

The Development of Turbulent Mixing Parameterization Scheme in the Surface Ocean under the Stratification and Earth’s Rotation

It is a project adopted for Research Proposal Based Use in FY 2019.

Numerical Approach to the particle and nuclear physics

核物理研究センターでは、サブアトミック世界の強い相互作用の基本原理である量子色力学（ＱＣＤ）を中心に、素粒子・ハドロン・原子核物理学とその関
連分野にわたる研究を、全国の共同利用研究者とともに行っている。スーパーコンピューターは、大規模計算を必要とする理論的な研究を中心に使用されて
いる。その内容は、格子ＱＣＤによる数値シミュレーション、原子核構造及びハドロン物理現象論の数値的な解法を始め、深層学習の加速器実験への応用を
含む。

(i) 格子QCD
ハドロン分光学（質量などの計算）とハドロン間相互作用、クォーク間相互作
用、有限温度・密度における真空の性質・状態方程式、
閉じ込め・カイラル対称性の問題

(ii) 少数多体系・原子核構造
３-４核子系の精密計算、原子核密度汎関数法

(iii) 機械学習
加速器実験への応用（加速器の運転制御、Belle実験・ILC実験での信号識別・フレーバー識別、ILC SiDカロリメータのエネルギー較正）

(iv) その他
２次元格子フェルミオン系、超新星爆発、弦理論

Atomistic prediction of the temperature- and loading-rate-dependent first pop-in load in nanoindentation

材料のナノインデンテーション試験中に観察され、局所領域での材料の機械的強度に根本的に関連する最初の変位バースト現象(第一pop-in)の発生荷重の温度・押し込み速度依存性を予測することを目的とし、896000原子を含むBCC解析モデル結晶をEAMポテンシャルで駆動してナノインデンテーションの分子動力学解析を実施した。解析中にて圧子直下に生じる原子応力状態を獲得し、その結果を基にNudged Elastic Band法を用い、第一pop-inに起因する転位核生成の活性化エネルギーを求めた。そして確率モデルから異なる温度・押し込み速度下での第一pop-inの確率分布を求めることで予測を行った。予測した第一pop-in発生荷重は温度が高くなるほど、また押し込み速度が遅くなるほど低くなり、実験での温度・押し込み速度依存性の傾向と一致していることから、本手法によって第一pop-in発生荷重が予測可能であることを確認した。

Laser micro machining using a photonic nanojet

フォトニックナノジェット(Photonic nanojet: PNJ)は誘電体マイクロ球にレーザを照射することで発生するビームである．PNJの高分解能な強度分布制御によって，高分解能に加工幅や加工深さを制御可能なレーザ微細加工技術の確立を目的とする．加工対象であるシリコンにPNJを照射した際の強度分布をFinite-Difference Time-Domain法を用いた電磁場シミュレーションによって解析した．これによって，PNJを用いることでサブマイクロメートルスケールの加工が可能であることを解析的に明らかにした．

Synthesizing Real-Scale Populations with Location and Income Atributes of Whole Japan using Birth Cohort

It is a project adopted for Research Proposal Based Use in FY 2019.

Helicity instability effects of supersonic Batchelor vortices

超音速流では変動の成長が圧縮性効果によって抑制されるが，線形安定性理論を用いてヘリシティ分布による不安定効果を見出した。その効果を非線形発達も含めた高次精度3次元数値計算でも確認した。

Simulation Study of Collective Phenomena of Plasmas via Molecular Dynamics/N-body Approach

クーロン相互作用に支配されるプラズマの集団的ダイナミクスとその背後にある階層構造を明らかにすることを目的として，ミクロな運動方程式からプラズマのマクロ現象を再現するシミュレーションプログラムを開発した．大規模並列粒子シミュレーションコード開発フレームワークFDPSを用いてコード開発を行い，プラズマ中の基本的集団現象であるデバイ遮蔽のシミュレーションによってコードの動作確認を行った．周期境界条件下で系内の電子の運動を解き，イオンのポテンシャルが遮蔽される現象を再現した．周期境界条件の実装のために，Particle-Mesh 法と近傍領域のイメージ粒子を考慮する方法を実装し比較した．

Numerical Approach to the particle and nuclear physics

核物理研究センターでは、サブアトミック世界の強い相互作用の基本原理である量子色力学（ＱＣＤ）を中心に、素粒子・ハドロン・原子核物理学とその関連分野にわたる研究を、全国の共同利用研究者とともに行っている。スーパーコンピューターは、大規模計算を必要とする理論的な研究を中心に使用されている。その内容は、格子ＱＣＤによる数値シミュレーション、原子核構造及びハドロン物理現象論の数値的な解法を始め、深層学習の加速器実験への応用を含む。

(i) 格子QCD
ハドロン分光学（質量などの計算）とハドロン間相互作用、クォーク間相互作用、有限温度・密度における真空の性質・状態方程式、閉じ込め・カイラル対称性の問題

(ii) 少数多体系・原子核構造
３-４核子系の精密計算、原子核密度汎関数法

(iii) 機械学習
加速器実験への応用（加速器の運転制御、Belle実験・ILC実験での信号識別・フレーバー識別、ILC SiDカロリメータのエネルギー較正）

(iv) その他
２次元格子フェルミオン系、超新星爆発、弦理論

Lattice Boltzmann Simulation on Reaction and Transport Processes in Large-scale Porous Media

大規模多孔質場における反応輸送解析を対象として，格子ボルツマン解析の適用可能性ならびに並列化による計算性能向上について評価した．その上で，反応流動解析を実施し，多孔質構造内局所反応分布の時系列変化について検討を行った．

Numerical Simulation on Mixing Control of Free Jet

DNS(Direct Numerical Simulation)により，複数の自由噴流を配置し壁面に衝突させる衝突自由噴流について,多数配置した場合の流動解析,噴流自体を周期的に移動（振動）させる制御,開花噴流による制御を行い,伝熱特性を評価した．19本の多重衝突自由噴流を解析した結果,多重衝突噴流の伝熱性能の向上のためには吹き上げの発生を抑制することが必要であること,１列に並んだ4本の噴流を列方向に周期的に振動させることで比較的一様な流れ場と伝熱性能が得られること,開花噴流を特定の周波数と衝突距離の下で壁面に衝突させることで,ほぼ一定の均一な伝熱特性が得られることを明らかにした．