学際大規模情報基盤共同利用・共同研究拠点 平成 23 年度共同研究 最終報告書 2012 年 5 月

課題番号 11-NA17

原子衝突による材料科学のための大規模シミュレーション基盤
青木学聡（京都大学）佐竹信一（東京理科大学）剣持 貴弘（同志社大学）
岩下 武史（京都大学） 松尾二郎（京都大学） 瀬木利夫（京都大学）
概要
イオンビーム, プラズマ, 超音速ジェットによる材料加工プロセスの詳細な解析を行う
ため, 並列計算機上で原子衝突現象を扱うための分子動力学シミュレーションプログラ
ム, そしてこれに付随する多量の原子座標データを扱う解析環境の開発を行っている．
本年度は，原子衝突現象の分子動力学シミュレーションの周辺を支える，プリ/ポストプ
ロセスの並行実行，および衝突領域が重なりあう逐次衝突現象の検証の際に必要とされるタスクの同時実行処理を，Xcrypt を用いたスクリプトにより実装した．また，これらを用い，多数の原子の衝突の繰り返しによるスパッタリング現象のシミュレーションを実施し，得られた原子座標のデータから原子衝突に起因する欠陥の形成過程および表面
加工プロセスの検証を行った．

1.

最大 1 億(108)原子の標的材料を扱う，あるいは 1

研究の目的と意義

万サンプル以上の衝突現象を容易に実現するため(1) 研究の目的の原子衝突シミュレーション基盤を開発すること

本課題では，数 eV から MeV という幅広い範囲が求められる．また，大規模分子動力学シミュレ

での粒子線と固体材料表面との衝突素過程を解明

ーションにより得られる，極めて多量の原子位置するための，大規模分子動力学法を中心とした原情報に対する多面的な解析を可能とし，さらにマ子シミュレーション手法の技術展開を図ることをルチスケール，マルチフィジクスシミュレーショ

目的とする．入射エネルギーが数 keV から MeV 領ンとの連携を行うための，並列・分散処理，他階域での原子・固体衝突では，高効率の二次イオン層の記憶領域等を前提とするソフトウェア基盤の脱離や局所的な材料変性など特異な現象が見られ整備を進める．原子の脱離過程や表面構造等，原る．この現象の解明には，常温に比べて 5~8 桁と子座標を一次情報とする特徴的なパラメータの空

いう高いエネルギーの局所領域への導入とその散間，時間的変化の抽出により，粒子衝突シミュレ

逸，そしてこれに伴う大規模な標的原子の移動がーションデータに対する解析方法を提示し，ナノサブナノメートル～マイクロメートルの空間スケ

スケール材料プロセス，計測技術の進展に資する．

ール，サブピコ～ナノ秒の時間スケールで追跡す(2) 研究の意義と必要性ることが必要となる．一方，運動エネルギー1keV

イオンビーム，プラズマ，超音速ジェット技術

以下の低エネルギーイオンビーム，プラズマプロ

は，常温以上の運動エネルギーをもった原子や分セス，超音速分子流の領域では，材料極表面にお子を精密に制御し，標的材料に衝突させる技術でける超微細加工，表面クリーニング等，多数の原あり，これらは半導体製造プロセスと筆頭とする

子衝突に繰り返しの後に生じる表面脱離反応や，

ナノテクノロジーにおける製造，計測技術に必須

原子種に応じた表面反応特性に関する，より詳細の手法である．ナノテクノロジーの発展に伴い，な理解が求められる．

これらの粒子線技術も更なる高度化，多様化が求これらの衝突現象を原子スケールで再現・解析

められるようになっている．本課題で実施する分するため，分子動力学法，モンテカルロ法による，

子動力学シミュレーションをはじめとする大規模1

学際大規模情報基盤共同利用・共同研究拠点 平成 23 年度共同研究 最終報告書 2012 年 5 月

原子衝突シミュレーションにより，エネルギー粒分野を主な研究領域としている．原子衝突に関す

子と固体表面の相互作用を原子レベルで詳細に調るシミュレーションプログラムの開発，データ解べることで，材料表面衝突素過程の解明と応用に

析を中心として，実際の実験および大規模計算手

関する知見が得られることが期待される．

法に関する助言を受けながら研究を進めている．
(2) 共同研究分野

分子動力学法は，系内に存在する原子間の相互作用を求めた後，ニュートンの運動方程式に従い，入射エネルギーが数 keV から MeV 領域での原全ての粒子の運動の時間発展を追跡する手法であ子・固体衝突では，高効率の二次イオン脱離や局る．この方法は，近年の計算機資源(計算速度，記所的な材料変性など特異な現象が見られる．この憶容量)の増大とともに，ナノスケール物理シミュ現象の解明には，常温に比べて 5～8 桁という高い

レーションの手法として，生命科学，材料科学分

エネルギーの局所領域への導入とその散逸，そし野において広く利用されるようになっている．現てこれに伴う大規模な標的原子の移動がサブナノ状においても，分子動力学シミュレーションソフ

メートル～マイクロメートルの空間スケール，サトウェアは市販品およびオープンソースソフトウ

ブピコ～ナノ秒の時間スケールで追跡することがェアとして多数存在するが，その多くは比較的小必要となる．このようなスケールでの大規模原子規模，静的な状態を対象としたものに限られてい集団の運動を解析する手法として，分子動力学法

る．本課題が対象とする原子衝突過程の解析には， が挙げられる．本研究課題では，大規模並列計算
最大 108 個という巨大な標的材料，最大数桁に及ぶ

環境上における原子衝突現象の分子動力学シミュ

非熱平衡状態からの緩和過程等，原子衝突特有のレーションを主な対象とし，核となるシミュレー

現象に対応させる必要がある．また，分子動力学

ションプログラム及び，ポスト・プリプロセスをシミュレーションより得られる原子座標データを含めた，シミュレーション環境全体の高度化に取多量に収集することにより，電子顕微鏡イメージ， り組んでいる．
X 線回折，RBS 測定などの測定原理を直接再現す(3) 当公募型共同研究ならではという事項など

る他，マルチスケール，マルチフィジクスシミュ

本研究制度への応募を機に，シミュレーション

レーションモデルとのシームレスな連携が期待で技術を中心とし，実験，計算機理論の分野が相互きる．このような巨大原子座標データを対象としに連携できる研究コミュニティの運営を進めていた多面的な解析手法を容易に実現するためには，る．また，シミュレーションプログラムの開発及並列，分散処理，複数の記憶階層を前提とした原びその実施に際し，京都大学学術情報メディアセ子座標データへのアクセスモデルを設計する必要

ンターのスーパーコンピュータシステムの提供をがある．これらの背景より，本研究が対象とする受けた．最大 4 ノード 64CPU からなる計算機資源

シミュレーション技術，データ解析技術は，加速を年間を通じて利用できるという支援により，単器物理学，表面材料科学，計算物理工学，計算機

独の PC では時間的(繰り返し数)，空間的(扱う原科学の連携により実現される，超大規模数値計算子数)に困難な規模の原子衝突シミュレーション

系応用分野の課題に相当するものと考える．

を実施している．

2. 当拠点公募型共同研究として実施した意義

3.

(1) 共同研究を実施した大学名と研究体制(1)研究成果の詳細について

研究成果の詳細と当初計画の達成状況

本研究では京都大学，東京理科大学，同志社大

平成 22 年度での本研究課題において，分子動力学から 6 名の研究員による共同研究体制を立ち上学シミュレーションプログラムの大規模並列化にげた．このうち 3 名が原子衝突シミュレーション，

取り組み，108 原子系からなる原子衝突シミュレー

2 名が原子衝突現象の基礎実験と応用，
1 名が HPC

ションを実施することができた．そこで本年度は，
2

学際大規模情報基盤共同利用・共同研究拠点 平成 23 年度共同研究 最終報告書 2012 年 5 月

リスト 1: Xcrypt スクリプトの例

ポストプロセスにおける並行分散処理，およびジ

use base qw (core);

ョブ単位での並行実行の手法について検討を進め# リモートで実行する場合の環境を追加(オプション)
my $remote_kyoto_env = &add_host({'host'=>‘t123@thin.kudpc.kyoto-u.ac.jp','sched' => 't2k_kyoto','wd' => '/home/t/t123/multislice.xy','xd' => '/home/t/t123/local/opt/xcrypt'});た．
分子動力学シミュレーションにおいてシミュレ
ーションプログラムから得られる情報は，多数の# テンプレートの記述
my %template = ('id' => 'job',# ジョブ毎に変化させるパラメータ‘RANGE0’ => [0..59],# 実行するユーザー定義スクリプト# $VALUE[0] が 上記RANGE0の値で置換される'exe0@' => sub { "./asi_auto_byid_200keV.sh$VALUE[0] || exit 1"; },# リモートジョブのオプション'JS_cpu' => 1,'JS_node' => 1,'JS_memory' => '1000MB','JS_limit_time' => '30:00','JS_queue' => 'qh123,'JS_group' => 'qh123','env' => $remote_kyoto_env,);

原子の位置と速度のみである．一方，種々の実験より観察される物理量は，表面形状，格子周期，
原子結合状態といった，原子座標からみれば間接的な情報である．シミュレーション結果の妥当性に関する検討や，多面的な解析を行うためには，核となるシミュレーションプログラムの高速化は
もとより，その出力結果の収集，および解析過程# テンプレートから個別のジョブを生成
my @jobs = prepare(%template);の迅速化，効率化も合わせて重要となる．
様々なプラットホーム上で複数タスクの並行実# ジョブの実行、 終了を待機
submit_sync(@jobs);

行を実現するツールとして，京都大学学術情報メディアセンターにおいて開発されている Xcrypt 1

マルチスライス法 3 を用いて透過型電子顕微鏡がある．Xcrypt は Perl による DSL であり，その記(Transmitting Electron Microscope, TEM)イメージを

載内容はリスト 1 に示すように，並行に実施する生成する．マルチスライス法の概要は以下の通り

命令をテンプレートとして記述する．Xcrypt はテである．

ンプレートから逐次的な変数(リスト 1 中では

1. 標的材料をオングストロームレベルの薄い層

RANGE0)を展開，個別の命令(ユーザー定義部は毎に分解する

exe0@で記載される)を生成し，これらを個々の環

2. 各層の原子配置から，電子線の散乱の伝達関境に対応した形式で実行する．ここでいう環境と

数を求める

は
「ローカルまたは ssh によるリモート実行」「shell3. 各層の重ね合わせにより，試料全体の伝達関または NQS, PBS 等へのバッチ投入」を指し，こ数を求める

れらを任意の組み合わせで設定できる．

4. 試料透過後の電子線と対象とする電子顕微鏡

Xcrypt を利用した MD シミュレーションのポスの光学定数より，結像したイメージを作成す

トプロセスとして，シリコン(Si)の集束イオンビー

る

ム(Focused Ion Beam)加工前後の電子顕微鏡像の

ここで 1~3 の処理は，対象となる原子数に従い増生成を取り上げた．

大する．このようなポストプロセス処理について

まず，FIB 衝突過程の分子動力学シミュレーシ

も，プログラムレベルでの並列化により高速化がョンを京都大学 HX600 上で実施し，原子座標デー期待されるものの，ポストプロセスの多様性，即タを生成した．21.72nm×21.72nm×25.9nm，644400

応性等の問題から対応することは困難である．こ個の Si 原子からなる標的材料に対し，
Ar 原子を逐のような場合，Xcrypt に代表されるような簡便な次的に衝突させる. Ar イオンの加速エネルギーは手法により，大規模バッチシステム上での並行実

50keV とし，これを基板表面中心の 4.0nm×4.0nm

行環境が利用できることは，解析時間の大幅な短の領域に連続的に衝突させる2．1 回の衝突毎に原縮につながる．なお，Xcrypt 導入時に行った解析

子座標データを収集した後，原子座標データより，
1 http://super.para.media.kyoto-u.ac.jp/xcrypt/index.html
23 E. J. Kirkland, "Advanced Computing in Electron

詳細は平成 23 年度中間報告を参照

Microscopy" (Plenum 1998, Springer 2010)3

学際大規模情報基盤共同利用・共同研究拠点 平成 23 年度共同研究 最終報告書 2012 年 5 月

(a)

0

(b)

1
(6.25×1012/cm2)(c)10(6.25×1013/cm2)

TEM Image Simulation (E=200keV)

Snap

Dose

図 1: 50keV Ar イオンの Si への逐次衝突シミュレーションの結果. 上段のスナップショットでは第 4 象限を切り欠いている．

プログラム側の変更は，Xcrypt スクリプトのテン

本質的に含まれることから，TEM イメージ像は原

プレート中で生成される通し番号(リスト 1 中の子配置の乱れに対し極めて敏感な応答を示す．図

$VALUE[0])に対し，これを入出力ファイル名とい1(b)に示すように，1 回の原子衝突時において入射った実行パラメータに変換する部分をシェルスク

原子および入射原子との衝突により高速ではじき

リプトのレベルで追加したのみであった．

出された原子の軌跡に沿って形成される欠陥構造

50keV の Ar 原子を逐次的に Si(001)基板に衝突が明確に観察される．また，衝突回数の増加に伴させたときの原子配置スナップショット(図中上い，基板中心部の結晶構造が乱されていき，図 1(g)段)，および TEM イメージシミュレーションの結の 60 回衝突(約 3.75×1014/cm2 相当)ではロッド状果(中及び下段)を図 1 に示す．図からわかるようのアモルファス領域が形成されるとともに，衝突に，原子座標をそのまま 3 次元投影しただけでは，箇所の表面が大きく盛り上がることも観察された．

原子間の結合状況や，本来の結晶構造からの歪み

Xcrypt によるタスク並行実行のもう一つの応用を見積もることは困難である．これに対し，電子として，巨大材料に対する原子衝突現象の解析を顕微鏡の撮像原理には原子配置のフーリエ変換が検討した．模式図を図 2 に示す．巨大材料系に対4

学際大規模情報基盤共同利用・共同研究拠点 平成 23 年度共同研究 最終報告書 2012 年 5 月(d)20(1.25×1014/cm2)(e)40(2.5×1014/cm2)(f)60(3.75×1014/cm2)

TEM Image Simulation (E=200keV)

Snap

Dose

(図 1 続き)する原子衝突シミュレーションを計算機資源の許3. 取り出した原子座標データに対し，境界条件

す範囲で実現するには，局所的な時刻の同時性にを適切に設定し，原子衝突シミュレーション

関する厳密性を緩和，すなわち，個々の原子衝突を実施

による損傷形成・緩和過程については完全に独立

4. 更新された原子座標データを標的全体の座標であると仮定する．この場合，大きな空間スケー

データベースに書き戻す(チェックイン)

ルに渡って生じるような，極めて大規模な原子移

これらの手続きは，原子衝突による影響が重なり動については再現することはできないが，原子衝合わない範囲で並行に実施できる．逆に原子衝突

突の繰り返しによる逐次的なノックオン現象によが重なりあう条件については，衝突の順番を維持

る原子移動については検証することが可能である． するために，後続のシミュレーション実行をブロ
具体的な手続きは以下のようになる．

ックする必要がある．さらに，何らかの理由で衝

1. ランダムに衝突箇所を選択

突シミュレーションが失敗した場合，後続のシミ

2. 標的材料全体を構成する原子座標データベーュレーションの実行をスキップする必要も生じる．

スから，衝突箇所近傍の原子座標データのみ

このような依存性のあるジョブを制御するためを取り出す(チェックアウト)に，開発者である京都大学学術情報メディアセン
5学際大規模情報基盤共同利用・共同研究拠点 平成 23 年度共同研究 最終報告書 2012 年 5 月

リスト 2: 依存関係のあるジョブ制御のための Xcrypt
拡張モジュール

巨大標的材料

package dependency_abortable;
use strict;use NEXT;
use common;use builtin;

衝突近傍
のみを
取り出し
並行計算

衝突領域が
重なる場合はブロック#テンプレートに ‘depend_on’キーワードを追加&add_key('depend_on');# job をsubmitする前に実行するルーチン
sub before {my $self = shift;
my @dep_jobs0 = @{mkarray ($self->{depend_on})};my @dep_jobs = ();foreach my $j0 (@dep_jobs0) {if ( ref $j0 eq '' ) {my $j = jobsched::find_job_by_id ($j0);
if ($j) {push (@dep_jobs, $j);}} elsif ( ref $j0 eq 'user' ) {push (@dep_jobs, $j0);}}
sync (@dep_jobs); # 依存するジョブによって実行をブロックする# abort myself if any of @dep_jobs is aborted
foreach (@dep_jobs) {if ( jobsched::get_job_status($_) eq 'aborted' )core::abort($self);
last;
}
}
}

図 2: 巨大標的材料に対する原子衝突シミュレーショ
ンの模式図

ターにより，下記の機能を付与する Xcrypt の拡張モジュールが追加された．
・ ジョブの依存性を登録するための’depend_on’
キーワードの導入
・ ジョブを submit する直前に呼び出されるサブ
ルーチン before にて，depend_on キーワードに{1;登録されたジョブ ID の終了を待機させる機能
・ 同じく before ルーチンにて，依存するいずれ

リスト 3: 依存関係のあるジョブ制御のための Xcryptの例

かのジョブの終了ステータスが aborted である

use base qw (dependency_abortable core);

場合, 自分自身を abort する機能

# 依存関係(兼ジョブパラメータ)ファイル# structure# [ { 'param' : { 'x':f, 'y':f, 'in_rad':f, 'out_rad':f },#'graph" : { 'linkfrom': [d], 'linkto': [d] } } ]use YAML;
my $y = YAML::LoadFile(“depend_graph.yaml");

拡張モジュールとこれを使った Xcrypt スクリプトの例をリスト 2, 3 に示す．リストに示すように，
極めて簡潔な記述により，ジョブ間依存関係の定

my $params = $#$y + 1;
my $ids = [0..$params-1];

義，アボートメッセージの伝達が実現される．

my %template = ('id' => 'job','RANGE0' => $ids,'exe0@' => sub { "./mysh.sh $VALUE[0] || exit 1"; },);Xcrypt は，拡張元である Perl の特徴を生かし，モ
ジュール毎の独立性が高く設計されていることか

my @jobs = prepare(%template);

ら，ユーザーは特定の機能を取捨選択することがの膨大なライブラリを活用できる．例えばリスト# ジョブ依存関係を登録
foreach my $i ( @{$ids} ) {$jobs[$i]->{depend_on}= [map {$jobs[$_];} @{$y->[$i]->{graph}->{linkfrom}}];}3 ではジョブ間の依存性は YAML 形式のファイル

submit_sync(@jobs);可能である．また，これまでに開発されてきた Perl

で与えるようになっている．このように，汎用的衝突シミュレーションを並行実施した．それぞれなプログラミング言語およびデータ構造基盤を最のシミュレーションジョブでは，「高速 Xe 原子の大限活用することで，スクリプトの作成，ユーザ

衝突・緩和(10 万ステップ，4ps 相当)」「冷却(5 万ー教育のコストを大幅に削減できる．

ステップ，40ps 相当)」と，2 つのシミュレーショ

実際に，リスト 2, 3 に示すスクリプトを利用し，ンを続けて実施した．この時 Si 或いは Xe 原子が実際に単一標的に材料の広領域への連続衝突シミ

該当する円筒領域外に出た場合，これらの原子はュレーションを実施した．表面が 54nm×54nm，「無効」のフラグを付与した状態でチェックイン

深さ 16nm，2,400,000 原子からなる Si(100)表面にを行い，以後チェックアウトの対象とならないよ

対し，1keV に加速した Xe 原子を垂直に複数回逐うにした．

次的に衝突させた．この際，衝突近傍より半径 7nm

今回，衝突箇所のパラメータ，および衝突領域の円柱領域を図 2 に示すように個別に取りだし，6

学際大規模情報基盤共同利用・共同研究拠点 平成 23 年度共同研究 最終報告書 2012 年 5 月

シミュレーション前

140回 (8.6×1012/cm2相当)衝突後

第1層

第2層

第3層

第4層

図 3: 1keV の Xe 原子を Si(100)表面に対し逐次的に衝突させる分子動力学シミュレーションの結果．

の依存関係グラフを作成には，以下のような単純される．また衝突の状態によっては，上部 4 原子な方法をとった．

層に相当する領域の Si 原子が移動することで，そ・j 番目のジョブ担当領域が i 番目のジョブ(ただしれより深い領域の Si 原子が露出する，といった深i<j)と重なる場合，
j の linkfrom の要素に i を付加．

いクレーター状のダメージが形成されることも示

・リスト 3 の要領で Xcrypt スクリプト内においてされた．同時に，Xe の衝突の影響により表面近傍

依存関係が構築される．

の Si 原子が真空中へ脱離する現象(スパッタリン

この方法では，依存関係のリンクは単純にジョグ)も観察された．このようなシミュレーションを

ブ数の 2 乗で増加するため，ジョブ数の増加に従通して，Xe 原子 1 回あたりの衝突による脱離数のい Xcrypt は多量のメモリが要求される結果となっ他，脱離する原子が直前にいた位置情報を得るこた．そこで，今回のシミュレーションではジョブとにより，スパッタリングモデルの検討を行うこを 20 毎に分割し Xcrypt のスクリプトを実行するとが可能となった．

こととした．より多量のジョブを一度に取り扱う(2) 当初計画の達成状況についてためには，依存関係のリンクを直近のもののみに限定する，といった簡略化の手続きをあらかじめ

分子動力学シミュレーションを原子衝突現象に実施する必要がある．

適用する場合の問題として，原子の衝突箇所近傍

Si(100)表面への Xe 原子の衝突シミュレーショ

にのみ極めて高速の粒子が多数存在し，その影響ンを 140 回実行した結果を図 3 に示す．この照射により系全体の時間刻みを短く取らざるを得ない，

量はおよそ 8.6×1012/cm2 に相当する．図 3 では Siことが挙げられる．この点を改善すべく，本年度

表面の第 1 層から 4 層までが見分けられるよう異は，MD シミュレーションのコアプログラム内でなる色で表現している．図より，Xe 原子の衝突箇の動的な時間刻み配分，並列コアへの資源割り当所において原子配置が大きく乱されることが観察ての調整法の検討を当初の目標とした．しかし，7

学際大規模情報基盤共同利用・共同研究拠点 平成 23 年度共同研究 最終報告書 2012 年 5 月

リソース配分およびプリ/ポストプロセスへの応

ログラム，および周辺を支える解析ツール群の高用の観点から，Xcrypt を活用したジョブレベルで機能，高性能化への要求がより明確になった．平の並列化によるアプローチという予定とは異なる

成 21 年度より本事業で実施した，コアとなるシミ手法をとることとなった．Xcrypt という軽量言語

ュレーションプログラム、解析プロセスの効率化のフレームワークを用いることで，タスクの並行に加え、今後は特に高速かつ高可搬性を備えた巨

実行，および依存性のあるタスクの効率的な制御

大座標データを扱うためのソフトウェア基盤の構が極めて容易に実現できた．このことは，本年度

築に重点を置くことで，それぞれの基盤要素の連の目標とした「大規模な系の中から注目すべき領携を図り，原子衝突現象におけるシミュレーショ

域を抽出し，この部分の精度よく解析する」といン技術の展開に資していきたい．

う目的に合致したものであり，本年度の目標につ5.

いては一応の達成を見たといえる．

研究成果リスト(1) 学術論文

4.

T. Aoki, T. Seki, J. Matsuo, “Molecular dynamics

今後の展望
大規模分子動力学シミュレーションにより得ら

simulations of large fluorine cluster impact on silicon

れた多量の原子座標データに対する，効率的なポwith supersonic velocity” Nucl. Instr. and Meth. B 269

ストプロセシングの方法として，Xcrypt による電(2011) 1582-1585.

子顕微鏡イメージシミュレーションプログラムのT. Aoki, T. Seki, J. Matsuo, “Molecular dynamics

タスク並行実行を試みた．Xcrypt はローカルとリ

study of crater formation by core-shell structured

モート，それぞれのシステムにおけるジョブ実行

cluster impact” Instr. and Meth. B 269 in press

手順の違いを吸収する．そのため，これまで利用(available online on 2011).

してきたローカル環境でのデータ解析プログラム， (2) 国際会議発表
スクリプトをほとんど変更することなく，そのま

T. Aoki, T. Seki and J. Matsuo, “Molecular

ま遠隔のスーパーコンピュータ上にて分散並行処

dynamics study of species and structure effect on large

理を行うことができた．本手法は今回示した電子

cluster

顕微鏡イメージシミュレーションに限らず，各種(2011/5/9, Nice, France)

impact”

E-MRS

2011

SpringMeeting

データの大量一括変換，パラメータサーベイ等，

T. Aoki, T. Seki and J. Matsuo, “Large-scale

多様なプリ・ポストプロセスを迅速に分散並行処

molecular dynamics simulation of huge fluorine

理環境へ移行するための手段として，極めて有用

cluster impact on silicon” 21st MRS-Japan academicであると考える．

symposium (2011/12/19-21, Yokohama).

一方，巨大標的材料に対する衝突シミュレーシ

T. Aoki, T. Seki, J. Matsuo, T. Iwashita, T. Hiraishi,

ョンのように，シミュレーションタスクの実行に

S. Satake, T. Kenmotsu, “Concurrent job control of

前後の依存関係が存在する系についても，極めて

simulations for iterative radiation and post

シンプルなモジュールの導入により実現可能であprocessing”,

ることが示された．このことは個々のタスクがよComputer Simulations of Radiation Effects in Solids

り複雑に関連するシミュレーションシステムを設(2012/06/24-29, Santa Fe, USA)

計する場合においても，Xcrypt が柔軟に対応でき(3) 国内会議発表

ることを示すものであり，さらなる応用展開の可11th

International

Conferenceon青木学聡，瀬木利夫，松尾二郎「巨大クラスタ
ー衝突による表面加工プロセスの解析」第 1 回マ

能性を示すものといえる．
一方，ジョブレベルでの効率的な運用が可能と

ルチスケールマテリアルモデリングシンポジウム(2011/5/23-24, 大阪大学)なることにより，コアとなるシミュレーションプ
8