Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
山本 耀二郎*; 早川 頌*; 沖田 泰良*; 板倉 充洋
Computational Materials Science, 229, p.112389_1 - 112389_9, 2023/10
被引用回数:2 パーセンタイル:23.35(Materials Science, Multidisciplinary)ヘリウムバブルは、核融合炉条件下で生成する特有の微細構造である。これらのバブルは、自らの移動によって接近・融合し、微細構造や材料特性に大きな影響を与える。しかし複数の金属原子の移動が含まれるこれらのプロセスは、分子動力学(MD)では時間スケールの制約から扱うことができない。この研究では、自己発展原子論的キネティックモンテカルロ(SEAKMC)法により時間スケールを拡張し、Fe中のバブルの融合過程を再現した。He原子の微小な振動やFe原子の短距離変位など、活性化エネルギーが極めて低い些細なイベントを避けるため、SEAKMCに2つのアルゴリズムを導入した。He原子の微小な振動を回避するために二段階のサドルポイント検索を行い、Fe原子の短距離変位を避けるためにFe原子の変位距離に対する閾値を設定した。さらに、活性化エネルギーの上限を設定する別のアルゴリズムを追加することで、非現実的に高い活性化エネルギーを持つイベントの選択を防ぎ、MDの時間尺度よりも8桁長い秒までのシミュレーション時間において、ダンベルから楕円形への構造の変化を再現することに成功した。開発された手法は、軽元素を含む金属材料の微細構造を分析するのに効果的であり、実験と比較可能な時間スケールに到達できる唯一の手法である。
早川 頌*; 山本 耀二郎*; 沖田 泰良*; 板倉 充洋; 鈴木 克幸*
Computational Materials Science, 218, p.111987_1 - 111987_10, 2023/02
被引用回数:1 パーセンタイル:6.14(Materials Science, Multidisciplinary)On-the-fly kinetic Monte Carlo (kMC), a computational technique for atomistic simulations, has attracted attention because it increases the simulation timescale beyond that of molecular dynamics (MD) simulations while maintaining atomistic fidelity. However, for most kMC methods, when events with high and low activation energies coexist in the event list, trivial events with extremely low activation energies that do not essentially affect the phenomena of interest, so-called flicker events, are frequently selected, making it challenging to observe the key dynamics. In this study, we use Self-Evolving Atomistic kMC (SEAKMC), one of the on-the-fly kMC methods, to model the unstable-to-stable transformations of irregular three-dimensional self-interstitial-atom (SIA) clusters in Cu generated through collision cascade. By setting an activation energy threshold once every five steps, transformations into stable configurations are enhanced. The algorithm renders the simulation timescales one or two orders of magnitude longer than those possible with MD simulations. Further, the probability of transformations into stable configurations is increased by 40 times compared to that of the original SEAKMC method. In addition, we find that the stable configurations obtained by the transformation of the SIA clusters are mostly Frank loops. In summary, this new algorithm for the SEAKMC method helps to resolve the inefficiency of kMC methods resulting from the selection of flicker events and will aid the study of meso-timescale atomistic dynamics.
津川 聖人*; 早川 頌*; 沖田 泰良*; 愛知 正温*; 板倉 充洋; 鈴木 克幸*
Computational Materials Science, 215, p.111806_1 - 111806_8, 2022/12
被引用回数:5 パーセンタイル:35.40(Materials Science, Multidisciplinary)Molecular dynamics simulations were conducted to evaluate the interactions between an edge dislocation and a rigid, impenetrable precipitate in Cu by changing the distance between the glide plane of the dislocation and the center of the precipitate (). In these calculations, the precipitate was introduced as a super particle that moved according to the total force exerted by the matrix atoms on the precipitate atoms. When the center of the precipitate was close to the glide plane, an Orowan loop was formed around the precipitate after the dislocation detached, and the critical resolved shear stress (CRSS) was similar to the value evaluated by the results at
. However, when the glide plane was far from the center of the precipitate, either a vacancy loop or loops generated through the Hirsch mechanism were formed, depending on whether the center of the precipitate was below or above the glide plane. The magnitude of the CRSS was not symmetric about
. This study confirmed that it is necessary to analyze the CRSS by changing
to construct a predictive model for the hardening caused by the formation of lattice defects, and that precipitate hardening appears to be smaller than the value estimated using the results at
.
津川 聖人*; 早川 頌*; 岩瀬 祐樹*; 沖田 泰良*; 鈴木 克幸*; 板倉 充洋; 愛知 正温*
Computational Materials Science, 210, p.111450_1 - 111450_9, 2022/07
被引用回数:16 パーセンタイル:75.72(Materials Science, Multidisciplinary)Precipitation strengthening has been utilized to improve the properties of metallic materials so far. Since interactions between precipitates and dislocations are micro-mechanisms responsible for this phenomenon, a molecular dynamics (MD) simulation is a powerful tool for quantifying this phenomenon. In this study, we introduced a method to simulate a rigid and impenetrable precipitate against a direct contact with a dislocation using a single interatomic potential representing the bulk material. The total force exerted on all atoms in the precipitate region was divided by the number of atoms in the region. This average force was then applied to each atom in the region to simulate one super particle that moved depending on the total force exerted by the matrix atoms on the precipitate atoms. We used MD simulations to quantify the interaction of a precipitate with an edge dislocation. After the dislocation overcame the precipitate, an Orowan loop was formed along the outer circumference of the precipitate. The energy of the loop was 2.1 0.1 eV/b, which was higher than that obtained using the elasticity. The hardening caused by the precipitate was larger than that caused by voids of the same size. The proposed method can be applied to simulate interactions of precipitates with dislocations in any type of metallic material, especially when a dislocation bypasses a precipitate without changing its structure, except when a strong repulsive force acts between them.
森 承宇*; 松田 那由多*; 沖田 泰良*; 愛知 正温*; 板倉 充洋; 鈴木 克幸*
Materialia, 21, p.101371_1 - 101371_6, 2022/03
The nonlinear ultrasonic (NLU) technique is a nondestructive method for detecting nanostructure in crystalline materials. In this study, a method was developed to quantify the changes in NLU signals associated with nanostructure using molecular dynamics (MD). A nonreflective boundary, which reduces the computational cost to the first power of the wavelength, was used to achieve this. This method is distinct from previous studies using a conventional MD, for which the computational cost is proportional to the square of the wavelength. The nonreflective boundary eliminates the influence of reflected waves at the detection position by setting a buffer region at the end of the simulation cell opposite from the wave source, and periodically resetting the displacements and velocities of all atoms in this region. This method allows the introduction of elastic waves with wavelengths longer than the cell size, and only an extension of time is required, according to the extension of the wavelength, without increasing the cell size. Hence, it is possible to extend the NLU wavelength by approximately four orders of magnitude, which approaches the wavelengths used for inspections and, thus, to use MD to simulate the changes in the NLU signals induced by nanostructure. The NLU signal values obtained by the two methods were in good agreement for a perfect Fe crystal and a Fe crystal containing 1% monovacancies. No significant frequency dependence of the acoustic nonlinearity parameter was found at 0 K. This method will contribute to the development of an inspection technique based on scientific principles.
沖田 泰良*; 寺山 怜志*; 津川 聖人*; 小林 恵太; 奥村 雅彦; 板倉 充洋; 鈴木 克幸*
Computational Materials Science, 202, p.110865_1 - 110865_9, 2022/02
被引用回数:7 パーセンタイル:34.77(Materials Science, Multidisciplinary)In this study, a Neural Network Potential (NNP) using an Artificial Neural Network (ANN) was developed for Zr, which is used as fuel cladding material in light water reactors. The reference data were obtained through first-principles calculations of various quantities, such as strained hexagonal-closed-packed (hcp) cells, strained face-centered cubic cells, cells containing a vacancy, several vacancies, and surfaces and -surface energy on all five slip planes in the hcp structures. These data were converted to training data for the ANN, which were invariant to the rotation and translation of the atoms and independent of the number of atoms in the cells. The ANN was defined as a three-layer structure and the number of the nodes was set to 26-12-18-1. The NNP reproduced the first-principles calculations, particularly for the shear deformation, vacancy formation energy, surface energy, and
-surface energy, with much higher accuracy than any of the existing potentials that have been developed for classical molecular dynamics simulations. The NNP was applied to identify the formation process of c-type dislocation loops in Zr, which is a key microstructure responsible for abrupt increases in hydrogen absorption. The formation process was determined by the balance of the vacancy formation energy, surface energy and the
-surface energy on the basal plane, both of which were precisely reproduced only by the NNP developed in this study. The formation process was identified based on the atomistic behavior of the NNP.
寺山 怜志*; 岩瀬 祐樹*; 早川 頌*; 沖田 泰良*; 板倉 充洋; 鈴木 克幸*
Computational Materials Science, 195, p.110479_1 - 110479_12, 2021/07
被引用回数:9 パーセンタイル:43.25(Materials Science, Multidisciplinary)Austenitic stainless steels, which are used as incore structural materials in light water reactors, are characterized by an extremely low stacking fault energy (SFE) among face-centered cubic (FCC) metals. To evaluate the effects of SFE on defect formation under high-energy particle irradiation, molecular dynamics simulations were performed using the interatomic potential sets for FCC metals with different SFEs and a primary knock-on atom energy (E) of 100 keV at 600 K. The results show that the number of residual defects is independent of the SFE. However, the characteristics of self-interstitial atom (SIA) clusters do depend on the SFE. For clusters smaller than a certain size, the ratio of glissile SIA clusters decreases as the SFE increases, which is similar to the trend observed at the low E
. However, for larger clusters, which can be detected only at a high E
, the ratio of glissile clusters increases. These results correspond to static energy calculations, in which the difference in the formation energy between a Frank loop and perfect loop (
E
) for the small clusters decreases as the SFE increases. In contrast, for the larger clusters, the SFE dependence of
E
changes due to the shape restrictions of stable perfect loops. At a high temperature of 600 K, large vacancy clusters with stacking faults can be detected at E
= 100 keV, resulting in the enhanced formation of these clusters at lower SFEs. Furthermore, several of these clusters were similar to perfect loops, with the edges split into two partial dislocations with stacking faults, although the largest clusters detected at low E
s were similar to stacking fault tetrahedrons.
早川 頌*; 土井原 康平*; 沖田 泰良*; 板倉 充洋; 愛知 正温*; 鈴木 克幸*
Journal of Materials Science, 54(17), p.11509 - 11525, 2019/09
被引用回数:18 パーセンタイル:53.56(Materials Science, Multidisciplinary)We performed molecular dynamics simulations to evaluate the effects of stacking fault energy (SFE) on interactions between a screw dislocation and spherical voids in face-centered cubic (fcc) metals. It was observed that the frequency of the cross-slips is a critical factor affecting the interaction, with primarily three different interaction morphologies being observed: (1) the two partial dislocations detach from the void independently with a time lag, (2) the two partial dislocations detach from the void almost simultaneously on a single slip plane, and (3) the two partial dislocations detach from the void almost simultaneously while involving more than one cross-slip and a jog formation. The magnitude of the critical resolved shear stress (CRSS) increases in the order mentioned above. The CRSS values for interaction morphology (2), which was observed most frequently in this study, were in good agreement with those predicted analytically by adjusting the parameters dependent on the SFE. Based on the obtained results, we discussed the applicability of the analytical model for void hardening in fcc metals. The results of this work contribute significantly to the modeling of mechanical property degradation in irradiated metals.
早川 頌*; 沖田 泰良*; 板倉 充洋; 川畑 友弥*; 鈴木 克幸*
Journal of Materials Science, 54(16), p.11096 - 11110, 2019/08
被引用回数:13 パーセンタイル:44.33(Materials Science, Multidisciplinary)We performed molecular dynamics simulations of displacement cascades in FCC metals under Poisson's deformation using interatomic potentials differing in stacking fault energy (SFE), in order to investigate the effect of tensile strain on the SFE dependence of defect formation processes. There was no clear SFE dependence of the number of residual defects and the size distribution of defect clusters under both no strain and the applied strain, while the strain enhanced the defect formation to a certain extent. We also observed that the strain affected the formations of self-interstitial atom (SIA) clusters depending on their size and the Burgers vector. These results were consistent with the analysis based on the defect formation energies. Meanwhile, the number of SIA perfect loops was higher at lower SFE under both no strain and the applied strain, leading to an increase in the ratio of glissile SIA clusters with a decrease in SFE. Further, the absolute number of SIA perfect loops was increased by the applied strain, while the SFE dependence of the number of SIA perfect loops was not affected. These findings were associated with the difference in formation energy between an SIA perfect loop and an SIA Frank loop. The insights extracted from this study significantly contribute to the modeling of microstructural evolution in nuclear materials under irradiation, especially for low SFE metals such as austenitic stainless steels.
中西 大貴*; 川畑 友弥*; 土井原 康平*; 沖田 泰良*; 板倉 充洋; 鈴木 克幸*
Philosophical Magazine, 98(33), p.3034 - 3047, 2018/09
被引用回数:11 パーセンタイル:42.12(Materials Science, Multidisciplinary)原子炉内部で使用されるオーステナイト鋼は積層欠陥エネルギーが低いことで知られ、そのため積層欠陥を伴う照射欠陥クラスタが大きくなることが予想される。この傾向を定量評価するため、FCC金属のモデルポテンシャルにおいて積層欠陥エネルギーだけを変化させたものを用い照射欠陥生成の分子動力学シミュレーションを様々な温度において実行し、生成した照射欠陥クラスタの種類と数を調べた。その結果移動可能な格子間原子クラスタの生成数が積層欠陥エネルギーとともに大きく変化することが見出された。
早川 頌*; 沖田 泰良*; 板倉 充洋; 愛知 正温*; 鈴木 克幸*
Philosophical Magazine, 98(25), p.2311 - 2325, 2018/06
被引用回数:8 パーセンタイル:35.77(Materials Science, Multidisciplinary)BCC鉄における格子間原子クラスタは容易方向へのすべり運動とは別に垂直方向への保存的上昇運動を行い、これにより三次元的な運動を行う。この拡散の三次元性は転位への吸収確率に大きく影響し、材料のスエリングなどの照射劣化挙動を左右する。この保存的上昇運動を分子動力学計算とモンテカルロ計算を組み合わせた手法で解析し、移動障壁エネルギーを評価したところ、従来の弾性理論により予想されていた値よりはるかに大きい値が得られた。
土井原 康平*; 沖田 泰良*; 板倉 充洋; 愛知 正温*; 鈴木 克幸*
Philosophical Magazine, 98(22), p.2061 - 2076, 2018/05
被引用回数:22 パーセンタイル:68.58(Materials Science, Multidisciplinary)原子炉材料、特に炉内構造材のオーステナイトの照射脆化を評価する際に、オーステナイトの特徴である低い積層欠陥エネルギーの影響を評価する必要がある。本研究では積層欠陥エネルギーだけを変化させたFCC金属のモデルポテンシャルを用いて転位とボイドの相互作用形態を評価し、ボイドの大きさと積層欠陥エネルギーの値によって様々な相互作用形態があることを見出した。
沖田 泰良*; 板倉 充洋
日本原子力学会誌ATOMO, 59(12), p.712 - 716, 2017/12
原子力材料の照射下特性変化を原子スケールの現象に基づいてモデル化する分子シミュレーションは、極短時間・極微小スケールの現象を定量化・解明・再現するため、幅広く用いられるようになってきた。本稿では、(1)分子動力学法を用いた高速中性子との弾性衝突に伴う結晶欠陥形成過程の定量化、(2)分子動力学法を用いた結晶欠陥集合体-転位接触による転位チャネリング形成過程の解明、(3)分子動力学-キネッティックモンテカルロ融合法を用いた結晶欠陥集合体の3次元的運動・転位へ吸収される過程の再現、に関して近年の研究進展を紹介する。
浅利 圭亮*; Hetland, O. S.*; 藤田 智*; 板倉 充洋; 沖田 泰良*
Journal of Nuclear Materials, 442(1), p.360 - 364, 2013/11
被引用回数:24 パーセンタイル:85.13(Materials Science, Multidisciplinary)炉内構造材料であるオーステナイトについて、その特徴である低い積層欠陥エネルギーが照射脆化にどう影響するかを調べるため、積層欠陥エネルギーの異なる分子動力学ポテンシャルを用いて照射欠陥による脆化への影響を調べた。その結果、積層欠陥エネルギーが低い場合は転位が広く分裂し二本の独立な転位として照射欠陥から抜けていく過程が観察され、またその過程に必要な外部応力が低下することが分かった。この結果はオーステナイトにおける脆化のモデル化に必須な情報となる。
板倉 充洋; 蕪木 英雄; 山口 正剛; 沖田 泰良*
Acta Materialia, 61(18), p.6857 - 6867, 2013/10
被引用回数:151 パーセンタイル:98.75(Materials Science, Multidisciplinary)原子力材料は長年の中性子照射によって硬化する。これは金属材料が折れることなく曲がる塑性変形が、照射による材料変化によって阻害される現象であるが、これを解明するには塑性のメカニズムを原子スケールで明らかにする必要がある。それには塑性変形を担う転位線と呼ばれる格子欠陥の動きがどのように照射損傷により阻害されるを知る必要があり、これは実験で直接観察できないので大規模な量子計算が必要になる。本発表では量子計算によって初めてこの転位の移動阻害メカニズムを明らかにしたので報告する。大規模な計算が必要という課題の解決にあたっては、新たに考案した多階層計算手法を用い少ない原子数で多数の原子での計算に相当する精度を出すことを可能にしたことが上げられる。これによって転位の運動阻害プロセスを定量的にモデル化することが可能となり、照射硬化をシミュレーションで定量評価するための道が開けたと言える。
板倉 充洋; 沖田 泰良*
no journal, ,
金属材料は析出物や格子欠陥により転位の運動が阻害され硬化する。特に原子炉材料では中性子照射による格子欠陥の生成と成長によりボイド, 析出物クラスタなど様々な欠陥が生成し硬化に寄与する。こうした寄与を定量評価するために線形弾性理論を用いた評価式が提唱され、特定の状況においては分子動力学計算の結果をよく再現することが知られている。しかし原子論的なモデリングによる補完が必要となる場合もある。例えばFCC金属では転位線が二本の部分転位に分裂するため、この分裂幅とボイドの直径の比によって部分転位が独立にボイドから抜けていく場合と一本の転位として同時に抜けていく場合があり、CRSSが異なる。またBCC金属では螺旋転位がパイエルス障壁のため直線状となり、ピン止めを脱する過程はこの螺旋転位の熱活性化過程による移動に支配されるため単純な応力のつり合い条件でCRSSを評価することはできない。本発表では分子動力学計算を用いてこれらピン止めのCRSS臨界を評価した結果を発表する。
沖田 泰良*; 板倉 充洋; 中西 大貴*; 川畑 友弥*
no journal, ,
Molecular dynamics simulations were conducted to evaluate the influence of the stacking fault energy on defect formation by collision cascades in FCC metals. The simulations were performed for PKA energies of 1-50 keV by using six sets of EAM potentials that differed only in SFEs. Neither the number of residual defects nor their clustering behavior is affected by SFE, except for the mean size of vacancy clusters at highest PKA energy. The mean size increases as the SFE decreases because of enhanced formation of large vacancy clusters, which prefer to comprise stacking fault structure inside them. The ratio of glissile SIA clusters decreases as the SFE increases. At higher SFEs, both the number of Frank loops and perfect loops tend to decrease; instead, three-dimensional clusters form with higher densities. The effect of SFE on the number of Frank loops becomes apparent only at the highest PKA energy, where comparably large SIA clusters can be formed with a higher density.
板倉 充洋; 沖田 泰良*
no journal, ,
現在核廃棄物を減らす新型原子炉の設計が行われており、その原子炉で使用する燃料被覆管の安全性について原子力機構, 東京大学, 日立製作所で共同研究を行っている。これまでより強い照射条件において被覆管の寿命の決定因子である水素化物析出までの時間がどのように変化するかを理論的なモデルにより予測し、適切な材料を開発することが目的である。本発表では照射欠陥クラスタへの水素偏析が駆動現象になると仮定し、その仮定を検証するため第一原理計算で純Zr中の空孔クラスタや積層欠陥への偏析エネルギー等を評価した結果を公表する。得られた結果は、最終的な材料開発に不可欠なモデルのパラメータを与えるものであり、新型炉の設計に大きく貢献するものである。
早川 頌*; 沖田 泰良*; 板倉 充洋; Haixuan, X.*; Osetsky, Y.*
no journal, ,
本研究は、照射下原子炉構造材に特有の結晶欠陥である自己格子間原子集合体に着目し、これら集合体の成長・合体の誘発メカニズムである集合体の三次元的挙動の定量化を目的とする。本発表では、シミュレーション技術開発室がこれまで行なってきた集合体三次元的挙動に関するon-the-fly kinetic Monte Carloによる計算結果と、それらから抽出できる知見を発表する。
森 承宇*; 松田 那由多*; 沖田 泰良*; 板倉 充洋
no journal, ,
非線形超音波応答(NLU)により、結晶性材料の材質劣化を検出する技術が様々な工学分野で用いられてきた。近年では、軽水炉圧力容器鋼の照射脆化を検出する技術への展開も図られているが、科学的根拠に基づいた検査技術構築のためには、非線形成分変化と微細組織変化の対応を原子挙動に基づいて解析することが求められる。特に分子動力学(MD)はその強力な手法であるが、計算の都合上導入可能な波長が実験で用いる波長より4桁以上短い。本研究では無反射境界を用いて計算速度を4000倍にまで高速化し、検査で使用される波長と同オーダーの波長でナノ構造に起因するNLU信号の変化を定量化することを目的とする。