Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Suzudo, Tomoaki; Ebihara, Kenichi; Tsuru, Tomohito; Mori, Hideki*
Journal of Applied Physics, 135(7), p.075102_1 - 075102_7, 2024/02
Times Cited Count:0 Percentile:0.00(Physics, Applied)Fracture of body centred cubic (bcc) metals and alloys below the ductile-to-brittle transition temperature is brittle. This is theoretically explained by the notion that the critical stress intensity factor of a given crack front for brittle fracture is smaller than that for plasticdeformation; hence, brittle fracture is chosen over plastic deformation. Although this view is true from a macroscopic point of view, such brittle fracture is always accompanied by small-scale plastic deformation in the vicinity of the crack tip, i.e. crack tip plasticity. This short paper investigates the origin of this plasticity using atomistic modeling with a recently developed machine-learning interatomic potential of -Fe. The computational results identified the precursor of crack tip plasticity, i.e. the group of activated atoms dynamically nucleated by fast crack propagation.
Suzudo, Tomoaki; Ebihara, Kenichi; Tsuru, Tomohito; Mori, Hideki*
Zairyo, 73(2), p.129 - 135, 2024/02
Body-centered-cubic transition metals, such as Fe and W, cleave along the {100} plane. To find out the mechanism of this response, atomistic simulations of curved crack-fronts of bcc Fe were conducted at 0 K using an interatomic potential created by an artificial neural network (ANN) technique. We discovered that dislocations can be emitted from the curved crack fronts along the {110} crack plane, and this phenomenon explains why the cleavage is observed only along the {100} plane. In addition, the cleavage simulations along {100} at the elevated temperature were found to be accompanied by plasticity; namely, they represented more realistic fracture.
Yabuuchi, Kiyohiro*; Suzudo, Tomoaki
Journal of Nuclear Materials, 574, p.154161_1 - 154161_6, 2023/02
Times Cited Count:2 Percentile:19.69(Materials Science, Multidisciplinary)In nuclear materials, irradiation defects cause degradation of mechanical properties. In these materials, the relationship between dislocations and voids is particularly important for mechanical strength. Although only spherical voids have been studied in the past, this study focuses on faceted voids, which are observed simultaneously with spherical voids. In the current study, molecular dynamics was used to analyze the effect of faceted voids in the irradiation hardening of pure iron. Specifically, we clarified the difference in obstacle strength and interaction processes between spherical voids and faceted voids, and that even faceted voids show differences in interaction depending on their crystallographic arrangement with dislocations.
Suzudo, Tomoaki; Ebihara, Kenichi; Tsuru, Tomohito; Mori, Hideki*
Scientific Reports (Internet), 12, p.19701_1 - 19701_10, 2022/11
Times Cited Count:7 Percentile:56.17(Multidisciplinary Sciences)Body-centered-cubic (bcc) transition metals, such as -Fe and W, cleave along the {100} plane, even though the surface energy is the lowest along the {110} plane. To unravel the mechanism of this odd response, large-scale atomistic simulations of curved cleavage cracks of -Fe were conducted in association with stress intensity factor analyses of straight crack fronts using an interatomic potential created by an artificial neural network technique. The study provides novel findings: Dislocations are emitted from the crack fronts along the {110} cleavage plane, and this phenomenon explains why the {100} plane can be the cleavage plane. However, the simple straight crack-front analyses did not yield the same conclusion. It is suggested that atomistic modeling, at sufficiently large scales to capture the inherent complexities of materials using highly accurate potentials, is necessary to correctly predict the mechanical strength. The method adopted in this study is generally applicable to the cleavage problem of bcc transition metals and alloys.
Chen, J.*; Yoshida, Kenta*; Suzudo, Tomoaki; Shimada, Yusuke*; Inoue, Koji*; Konno, Toyohiko*; Nagai, Yasuyoshi*
Materials Transactions, 63(4), p.468 - 474, 2022/04
Times Cited Count:1 Percentile:10.52(Materials Science, Multidisciplinary)In situ electron irradiation using high-resolution transmission electron microscopy (HRTEM) was performed to visualize the Frank loop evolution in aluminium-copper (Al-Cu) alloy with an atomic-scale spatial resolution of 0.12 nm. The HRTEM observation along the [110] direction of the FCC-Al lattice, Frank partial dislocation bounding an intrinsic stacking fault exhibited an asymmetrical climb along the 112 direction opposed to those in the reference pure Al under an electron irradiation, with a corresponding displacement-per-atom rate of 0.055-0.120 dpa/s. The asymmetrical climb of the partial dislocation was described as pinning effects due to Cu-Cu bonding in Guinier-Preston zones by a molecular dynamics simulation.
Ebihara, Kenichi; Suzudo, Tomoaki
Metals, 12(4), p.662_1 - 662_10, 2022/04
Times Cited Count:3 Percentile:31.34(Materials Science, Multidisciplinary)Phosphorus atoms in steels accumulate at grain boundaries via thermal and/or irradiation effects and induce grain boundary embrittlement. Quantitative prediction of phosphorus segregation at grain boundaries under various temperature and irradiation conditions is therefore essential for preventing embrittlement. To develop a model of grain boundary phosphorus segregation in -iron, we studied the migration of a phosphorus atom in two types of symmetrical tilt grain boundaries (3[1-10](111) and 5[100](0-13) grain boundaries) using molecular dynamics simulations with an embedded atom method potential. The results revealed that, in the 3 grain boundary, phosphorus atoms migrate three-dimensionally mainly in the form of interstitial atoms, whereas in the 5 grain boundary, these atoms migrate one-dimensionally mainly via vacancy-atom exchanges. Moreover, de-trapping of phosphorus atoms and vacancies was investigated.
Wang, J.*; Hatano, Yuji*; Toyama, Takeshi*; Suzudo, Tomoaki; Hinoki, Tatsuya*; Alimov, V. Kh.*; Schwarz-Selinger, T.*
Journal of Nuclear Materials, 559, p.153449_1 - 153449_7, 2022/02
Times Cited Count:3 Percentile:52.93(Materials Science, Multidisciplinary)To study the effect of the content of chromium (Cr) in the tungsten (W) matrix on the vacancy formation and retention of hydrogen isotopes, the samples of the W-0.3Cr alloy were irradiated with 6.4 MeV Fe ions in the temperature range of 523-1273 K. These displacement-damaged samples were exposed to D gas at a temperature of 673 K. The addition of 0.3% Cr into the W matrix resulted in a significant decrease in the retention of deuterium compared to pure W after irradiation especially at high temperature. Positron lifetime for W-0.3Cr alloy irradiated at 1073 K was almost similar to that for non-irradiated one. These facts indicate the suppression of the formation of vacancy-type defects by 0.3% Cr addition.
Toyama, Takeshi*; Suzudo, Tomoaki; Nagai, Yasuyoshi*; 9 of others*
Journal of Nuclear Materials, 556, p.153176_1 - 153176_7, 2021/12
Times Cited Count:4 Percentile:47.47(Materials Science, Multidisciplinary)We performed a high-precision investigation of radiation-enhanced diffusion (RED) using electron irradiation and three-dimensional atom probe (3D-AP). Cu-Fe diffusion pairs were created using high-purity Fe and Cu as base materials, and irradiated by 2 MeV electron. Cu diffusion into the Fe matrix was observed at the atomic level using 3D-AP, and the diffusion coefficient was obtained directly using Fick's law. RED was clearly observed, and the ratio of diffusion under irradiation to thermal diffusion was enhanced at low temperature. RED was quantitatively evaluated using the reaction kinetics model, and the model which consider only vacancies gave a good agreement. This gave experimental clarification that RED was dominated by irradiation-induced vacancies. In addition, the direct experimental results on the effect of irradiation on the solubility limits of Cu in Fe was obtained; solubility limits under irradiation were found to be lower than those under thermal aging.
Gilbert, M. R.*; Arakawa, Kazuto*; Suzudo, Tomoaki; Tsuru, Tomohito; 26 of others*
Journal of Nuclear Materials, 554, p.153113_1 - 153113_31, 2021/10
Times Cited Count:46 Percentile:90.27(Materials Science, Multidisciplinary)Modelling will continue to do so until the first generation of fusion power plants come on line and allow long-term behaviour to be observed. In the meantime, the modelling is supported by experiments that attempt to replicate some aspects of the eventual operational conditions. In 2019, a group of leading experts met under the umbrella of the IEA to discuss the current position and ongoing challenges in modelling of fusion materials and how advanced experimental characterisation is aiding model improvement. This review draws from the discussions held during that workshop. Topics covering modelling of irradiation-induced defect production and fundamental properties, gas behaviour, clustering and segregation, defect evolution and interactions are discussed, as well as new and novel multiscale simulation approaches, and the latest efforts to link modelling to experiments through advanced observation and characterisation techniques.
Zhao, C.*; Suzudo, Tomoaki; Toyama, Takeshi*; Nishitani, Shigeto*; Inoue, Koji*; Nagai, Yasuyoshi*
Materials Transactions, 62(7), p.929 - 934, 2021/07
Times Cited Count:3 Percentile:20.03(Materials Science, Multidisciplinary)We succeeded in measuring the diffusion coefficient of Cu in Fe in a low temperature range that had not been measured so far. Since the diffusion couple, which is a general method for measuring the diffusion coefficient, can be applied only at high temperature, atom probe tomography and Cu precipitation rate theory were used in this study. The estimated diffusion coefficient was found to be more reliable than that obtained in previous studies. Therefore, it is considered that the estimation by the atom probe provided higher accuracy. Furthermore, the kinetic Monte Carlo simulation revealed that the diffusion coefficient estimated by this method tends to be slightly overestimated as the temperature decreases.
Suzudo, Tomoaki; Tsuru, Tomohito
AIP Advances (Internet), 11(6), p.065012_1 - 065012_7, 2021/06
Times Cited Count:5 Percentile:29.17(Nanoscience & Nanotechnology)In the current study, we analyzed the self-interstitial atoms (SIAs) in BCC molybdenum (Mo) and tungsten (W) in comparison with other BCC transition metals utilizing first-principles method; particularly, we focused on uncommon dumbbells, whose direction are inclined from 111 toward 110 on the {110} plane. Such a direction is not stable neither in the group 5 BCC metals (i.e., vanadium, niobium, and tantalum) nor in -iron. Our first-principles relaxation simulations indicated that inclined dumbbells were more energetically-favored than common 111 dumbbells in Mo, while this is not necessarily the case for W. However, under a certain degree of lattice strain, such as shear or expansive strain, could make inclined dumbbells more favored also in W, suggesting that the lattice strain can substantially influence the migration barrier of SIAs in these metals because inclined dumbbells generally have a larger migration barrier than 111 dumbbells.
Suzudo, Tomoaki
Materials Science Forum, 1024, p.87 - 94, 2021/03
Tungsten (W) is suitable for solid targets of spallation neutron source due to its high neutron yield. The prediction of radiation effects of W is, therefore, of importance; especially, the influence of solute elements are complex and are not clearly known to date. We discuss here the solute effects using the first principles and kinetic Monte Carlo calculations and show that rhenium (Re) and osmium (Os), which are nuclear transmutation products of W, can largely change the stability and mobility of radiation defects. Such influences of the solute elements seem to explain the unsolved mechanism of the microstructural evolution of W-based materials under irradiation.
Ebihara, Kenichi; Sugiyama, Yuri*; Matsumoto, Ryosuke*; Takai, Kenichi*; Suzudo, Tomoaki
Metallurgical and Materials Transactions A, 52(1), p.257 - 269, 2021/01
Times Cited Count:9 Percentile:45.14(Materials Science, Multidisciplinary)We simulated the thermal desorption spectra of a small-size iron specimen to which was applied during charging with hydrogen atoms using a model incorporating the behavior of vacancies and vacancy clusters. The model considered up to vacancy clusters , which is composed of nine vacancies and employed the parameters based on atomistic calculations, including the H trapping energy of vacancies and vacancy clusters that we estimated using the molecular static calculation. As a result, we revealed that the model could, on the whole, reproduced the experimental spectra except two characteristic differences, and also the dependence of the spectra on the aging temperature. By examining the cause of the differences, the possibilities that the diffusion of clusters of and is slower than the model and that vacancy clusters are generated by applying strain and H charging concurrently were indicated.
Suzudo, Tomoaki; Ebihara, Kenichi; Tsuru, Tomohito
AIP Advances (Internet), 10(11), p.115209_1 - 115209_8, 2020/11
Times Cited Count:13 Percentile:58.14(Nanoscience & Nanotechnology)The mechanism of their brittle fracture of BCC metals is not fully understood. In this study, we conduct a series of three-dimensional molecular dynamics simulations of cleavage fracture of -iron. In particular, we focus on mode-I loading starting from curved crack fronts. In the simulations, brittle fractures are observed at cleavages on the {100} plane, while the initial cracks become blunted on other planes as a result of dislocation emissions. Our modeling results agreed with a common experimental observation, that is, {100} is the preferential cleavage plane in bcc transition metals.
Suzudo, Tomoaki; Takamizawa, Hisashi; Nishiyama, Yutaka; Caro, A.*; Toyama, Takeshi*; Nagai, Yasuyoshi*
Journal of Nuclear Materials, 540, p.152306_1 - 152306_10, 2020/11
Times Cited Count:11 Percentile:76.86(Materials Science, Multidisciplinary)Spinodal decomposition in thermally aged Fe-Cr alloys leads to significant hardening, which is the direct cause of the so-called 475C-embrittlement. To illustrate how spinodal decomposition induces hardening by atomistic interactions, we conducted a series of numerical simulations as well as reference experiments. The numerical results indicated that the hardness scales linearly with the short-range order (SRO) parameter, while the experimental result reproduced this relationship within statistical error. Both seemingly suggest that neighboring Cr-Cr atomic pairs essentially cause hardening, because SRO is by definition uniquely dependent on the appearance probability of such pairs. A further numerical investigation supported this notion, as it suggests that the dominant cause of hardening is the pinning effect of dislocations passing over such Cr-Cr pairs.
Ebihara, Kenichi; Suzudo, Tomoaki
Proceedings of Joint International Conference on Supercomputing in Nuclear Applications + Monte Carlo 2020 (SNA + MC 2020), p.65 - 69, 2020/10
Phosphorus (P) is known as one of the elements which cause the grain boundary (GB) embrittlement in steels and its GB segregation is promoted by the increase of vacancies and self-interstitial atoms due to irradiation. Thus we have been developing the rate-theory model for estimating GB P segregation under several temperatures and irradiation conditions. Because the model does not include the trapping and de-trapping processes properly, however, the model cannot calculate GB P coverage which is measured by experiments. As for the de-trapping process, so far, we have considered the migration of a P atom in the GB region of 3 symmetrical tilt GB using molecular dynamics (MD). In the current study, we also simulated the P migration in 5 GB using MD and compared the result with that of 3. As a result, at 800K, it was found that a P atom cannot migrate in 5 without vacancies while a P atom can migrate between iron atoms in 3.
Du, Y.*; Yoshida, Kenta*; Shimada, Yusuke*; Toyama, Takeshi*; Inoue, Koji*; Arakawa, Kazuto*; Suzudo, Tomoaki; Milan, K. J.*; Gerard, R.*; Onuki, Somei*; et al.
Materialia, 12, p.100778_1 - 100778_10, 2020/08
In order to ensure the integrity of the reactor pressure vessel in the long term, it is necessary to understand the effects of irradiation on the materials. In this study, irradiation-induced dislocation loops were observed in neutron-irradiated reactor pressure vessel specimens during annealing using our newly developed WB-STEM. It was confirmed that the proportion of loops increased with increasing annealing temperature. We also succeeded in observing the phenomenon that two loops collide into a loop. Moreover, a phenomenon in which dislocation loops decorate dislocations was also observed, and the mechanism was successfully explained by molecular dynamics simulation.
Ebihara, Kenichi; Suzudo, Tomoaki
TMS 2020; 149th Annual Meeting & Exhibition Supplemental Proceedings, p.995 - 1002, 2020/02
Times Cited Count:1 Percentile:55.88(Materials Science, Multidisciplinary)Phosphorus (P) is known as an element which causes grain boundary (GB) embrittlement in steels. In addition, GB P segregation is promoted by the increase of vacancies and self interstitial atoms due to irradiation. Thus, the diffusion rate theory model for estimating irradiation-induced GB P segregation has been developed based on the atomic processes. Since the present model does not include the trapping and de-trapping processes at GBs, however, it cannot calculate the value which is directly compared with experimental results. In this study, we simulated the migration of a P atom in the 3(111) symmetrical tilt GB. In addition, by tracking the migration of the P atom, the diffusion barrier energy was evaluated. As a result, the diffusion barrier energy was almost the same as the P segregation energy of an interstitial site in the GB, and it was found that P atoms migrate via interstitial sites in the GB.
Tsuru, Tomohito; Wakeda, Masato*; Suzudo, Tomoaki; Itakura, Mitsuhiro; Ogata, Shigenobu*
Journal of Applied Physics, 127(2), p.025101_1 - 025101_9, 2020/01
Times Cited Count:7 Percentile:38.96(Physics, Applied)We explored softening/strengthening behavior in various solute in W matrix by density functional theory (DFT) calculations combined with solid solution model. As a result of DFT calculations for various solute, a clear trend was observed in the interaction energy between a solute and a screw dislocation, which has predominant influence on solid solution behavior. The predictions based on the solid solution model with DFT can reasonably reproduce the complicated softening/strengthening as a function of temperature and solute concentration. We conclude that this specific balance is the origin of macroscopic solid solution softening.
Suzudo, Tomoaki; Onitsuka, Takashi*; Fukumoto, Kenichi*
Modelling and Simulation in Materials Science and Engineering, 27(6), p.064001_1 - 064001_15, 2019/08
Times Cited Count:17 Percentile:62.18(Materials Science, Multidisciplinary)Plasticity of body-centered-cubic (BCC) metals at low temperatures is determined by screw dislocation kinetics. Because the core of screw dislocation in these metals has non-planar structure, its motion is complex and unpredictable. For example, although density functional theory (DFT) predicts slip on a { 110 } plane, the actual slip plane at elevated temperatures departs from the prediction, its mechanism having been a mystery for decades. Here we conduct a series of molecular dynamics simulations to track the screw dislocation motion and successfully reproduced the transition of the slip plane. We then devised an algorithm to scrutinize the activation of dislocation jump over the Peierls barrier and discovered the possible origin of this unexpected phenomenon, i.e., a large fluctuation leads to the kink-pair nucleation for the cross-slip jump without transition of dislocation core structure.