Local structural evolution of Fe54C18Cr16Mo12 bulk metallic glass during tensile deformation and a temperature elevation process: A molecular dynamics study

Hui Lung Chen, Chia Hao Su, Shin Pon Ju*, Shih Hao Liu, Hsin Tsung Chen

*Corresponding author for this work

Research output: Contribution to journalJournal Article peer-review

4 Scopus citations


The mechanical and thermal properties of Fe54C18Cr16Mo12 bulk metallic glasses (BMGs) were investigated by a molecular dynamics simulation with the 2NN modified embedded-atom method (MEAM) potential. The fitting process of the cross-element parameters of 2NN MEAM (Fe-C, Fe-Cr, Fe-Mo, C-Cr, C-Mo, and Cr-Mo) was carried out first by the force matching method (FMM) on the basis of the reference data from density functional theory (DFT) calculations. With these fitted parameters, the structure of Fe54C18Cr16Mo12 BMG was constructed by the simulated-annealing basin-hopping (SABH) method, and the angle distribution range of the X-ray diffraction profile of the predicted Fe54C18Cr16Mo12 BMG closely matches that of the experiment profile, indicating the fitted 2NN MEAM parameters can accurately reflect the interatomic interactions of Fe54C18Cr16Mo12 BMG. The Honeycutt-Andersen (HA) index analysis results show a significant percentage of icosahedral-like structures within Fe54C18Cr16Mo12 BMG, which suggests an amorphous state. According to the tensile test results, the estimated Young's modulus of Fe54Cr16Mo12C18 bulk metallic glass is about 139 GPa and the large plastic region of the stress-strain curve shows that the Fe54C18Cr16Mo12 BMG possesses good ductility. Local strain distribution was used to analyze the deformation mechanism, and the results show that a shear band develops homogeneously with the tensile fracture angle (T) at about 50 degrees, in agreement with experimental results 45° < T < 90°. For the temperature elevation results, the discontinuity of the enthalpy-temperature profile indicates the melting point of Fe54Cr16Mo12C18 BMG is about 1310 K. The diffusion coefficients near the melting point were derived by the Einstein equation from the mean-square-displacement (MSD) profiles between 800-1400 K. On the basis of diffusion coefficients at different temperatures, the diffusion barriers of Fe54Cr16Mo12C18 can be determined by the Arrhenius equation. The diffusion barriers of total for Fe, Cr, Mo, C are 31.88, 24.68, 35.26, 22.50 and 31.79 kJ mol-1, respectively. The diffusion barriers of Fe and Cr atoms are relatively lower, indicating Fe and Cr atoms more easily diffuse with the increasing temperature.

Original languageEnglish
Pages (from-to)103925-103935
Number of pages11
JournalRSC Advances
Issue number126
StatePublished - 2015
Externally publishedYes

Bibliographical note

Publisher Copyright:
© The Royal Society of Chemistry 2015.


Dive into the research topics of 'Local structural evolution of Fe54C18Cr16Mo12 bulk metallic glass during tensile deformation and a temperature elevation process: A molecular dynamics study'. Together they form a unique fingerprint.

Cite this