Research Papers

Meta-Atom Molecular Dynamics for Studying Material Property Dependent Deformation Mechanisms of Alloys

[+] Author and Article Information
Peng Wang

Institute of Applied Mechanics,
Zhejiang University,
Hangzhou 310027, China;
Institute of Advanced Engineering Structures
and Materials,
Zhejiang University,
Hangzhou 310027, China

Hongtao Wang

Institute of Applied Mechanics,
Zhejiang University,
Hangzhou 310027, China
e-mail: htw@zju.edu.cn

1Corresponding author.

Contributed by the Applied Mechanics Division of ASME for publication in the JOURNAL OF APPLIED MECHANICS. Manuscript received July 25, 2017; final manuscript received August 16, 2017; published online September 8, 2017. Editor: Yonggang Huang.

J. Appl. Mech 84(11), 111002 (Sep 08, 2017) (8 pages) Paper No: JAM-17-1404; doi: 10.1115/1.4037683 History: Received July 25, 2017; Revised August 16, 2017

Massively parallel molecular dynamics (MD) simulations have been performed to understand the plastic deformation of metals. However, the intricate interplay between the deformation mechanisms and the various material properties is largely unknown in alloy systems for the limited available interatomic potentials. We adopt the meta-atom method proposed by Wang et al., which unifies MD simulations of both pure metals and alloys in the framework of the embedded atom method (EAM). Owing to the universality of EAM for metallic systems, meta-atom potentials can fit properties of different classes of alloys. Meta-atom potentials for both aluminum bronzes and hypothetic face-centered-cubic (FCC) metals have been formulated to study the parametric dependence of deformation mechanisms, which captures the essence of competitions between dislocation motion and twinning or cleavage. Moreover, the solid-solution strengthening effect can be simply accounted by introducing a scaling factor in the meta-atom method. As the computational power enlarges, this method can extend the capability of massively parallel MD simulations in understanding the mechanical behaviors of alloys. The calculation of macroscopic measurable quantities for engineering oriented alloys is expected to be possible in this way, shedding light on constructing materials with specific mechanical properties.

Copyright © 2017 by ASME
Your Session has timed out. Please sign back in to continue.


Jacobsen, K. W. , and Schiotz, J. , 2002, “ Computational Materials Science: Nanoscale Plasticity,” Nat. Mater., 1(1), pp. 15–16. [CrossRef] [PubMed]
Yip, S. , 2004, “ Nanocrystalline Metals: Mapping Plasticity,” Nat. Mater., 3(1), pp. 11–12. [CrossRef] [PubMed]
Bragg, L. , and Nye, J. F. , 1947, “ A Dynamical Model of a Crystal Structure,” Proc. R. Soc. Lond. A. Math. Phys. Sci., 190(1023), pp. 474–481.
Georges, J. M. , Meille, G. , Loubet, J. L. , and Tolen, A. M. , 1986, “ Bubble Raft Model for Indentation With Adhesion,” Nature, 320(6060), pp. 342–344. [CrossRef]
Gouldstone, A. , Van Vliet, K. J. , and Suresh, S. , 2001, “ Nanoindentation: Simulation of Defect Nucleation in a Crystal,” Nature, 411(6838), p. 656. [CrossRef] [PubMed]
Daw, M. S. , and Baskes, M. I. , 1983, “ Semiempirical, Quantum Mechanical Calculation of Hydrogen Embrittlement in Metals,” Phys. Rev. Lett., 50(17), pp. 1285–1288. [CrossRef]
Daw, M. S. , and Baskes, M. I. , 1984, “ Embedded-Atom Method: Derivation and Application to Impurities, Surfaces, and Other Defects in Metals,” Phys. Rev. B, 29(12), pp. 6443–6453. [CrossRef]
Finnis, M. W. , and Sinclair, J. E. , 1984, “ A Simple Empirical N-Body Potential for Transition-Metals,” Philos. Mag. A, 50(1), pp. 45–55. [CrossRef]
Yamakov, V. , Wolf, D. , Phillpot, S. R. , Mukherjee, A. K. , and Gleiter, H. , 2004, “ Deformation-Mechanism Map for Nanocrystalline Metals by Molecular-Dynamics Simulation,” Nat. Mater., 3(1), pp. 43–47. [CrossRef] [PubMed]
Wang, P. , Xu, S. , Liu, J. , Li, X. , Wei, Y. , Wang, H. , Gao, H. , and Yang, W. , 2017, “ Atomistic Simulation for Deforming Complex Alloys With Application Toward TWIP Steel and Associated Physical Insights,” J. Mech. Phys. Solids, 98, pp. 290–308. [CrossRef]
Johnson, R. A. , 1989, “ Alloy Models With the Embedded-Atom Method,” Phys. Rev. B, 39(17), pp. 12554–12559. [CrossRef]
Daw, M. S. , Foiles, S. M. , and Baskes, M. I. , 1993, “ The Embedded-Atom Method: A Review of Theory and Applications,” Mater. Sci. Rep., 9(7–8), pp. 251–310. [CrossRef]
Mendelev, M. I. , Kramer, M. J. , Becker, C. A. , and Asta, M. , 2008, “ Analysis of Semi-Empirical Interatomic Potentials Appropriate for Simulation of Crystalline and Liquid Al and Cu,” Philos. Mag., 88(12), pp. 1723–1750. [CrossRef]
Rose, J. H. , Smith, J. R. , Guinea, F. , and Ferrante, J. , 1984, “ Universal Features of the Equation of State of Metals,” Phys. Rev. B, 29(6), pp. 2963–2969. [CrossRef]
Bain, E. C. , 1924, “ The Nature of Martensite,” Trans. Am. Inst. Min. Metall. Eng., 70(98), pp. 25–46.
Milstein, F. , Fang, H. E. , and Marschall, J. , 1994, “ Mechanics and Energetics of the Bain Transformation,” Philos. Mag. A, 70(4), pp. 621–639. [CrossRef]
Weixue, L. , and Tzuchiang, W. , 1998, “ Ab Initio Investigation of the Elasticity and Stability of Aluminium,” J. Phys.: Condens. Matter, 10(43), p. 9889. [CrossRef]
Plimpton, S. , 1995, “ Fast Parallel Algorithms for Short-Range Molecular Dynamics,” J. Comput. Phys., 117(1), pp. 1–9. [CrossRef]
Kittel, C. , 2004, Introduction to Solid State Physics, Wiley, Hoboken, NJ.
Pearson, W. B. , 1967, A Handbook of Lattice Spacings and Structures of Metals and Alloys, Pergamon Press, London.
Skriver, H. L. , and Rosengaard, N. M. , 1992, “ Surface Energy and Work Function of Elemental Metals,” Phys. Rev. B, 46(11), pp. 7157–7168. [CrossRef]
Rohatgi, A. , Vecchio, K. , and Gray , G., III , 2001, “ The Influence of Stacking Fault Energy on the Mechanical Behavior of Cu and Cu-Al Alloys: Deformation Twinning, Work Hardening, and Dynamic Recovery,” Metall. Mater. Trans. A, 32(1), pp. 135–145. [CrossRef]
Rohatgi, A. , Vecchio, K. S. , and Gray, I. G. T. , 2001, “ A Metallographic and Quantitative Analysis of the Influence of Stacking Fault Energy on Shock-Hardening in Cu and Cu–Al Alloys,” Acta Mater., 49(3), pp. 427–438. [CrossRef]
Zhang, Y. , Tao, N. R. , and Lu, K. , 2009, “ Effect of Stacking-Fault Energy on Deformation Twin Thickness in Cu–Al Alloys,” Scr. Mater., 60(4), pp. 211–213. [CrossRef]
Kibey, S. , Liu, J. B. , Johnson, D. D. , and Sehitoglu, H. , 2006, “ Generalized Planar Fault Energies and Twinning in Cu–Al Alloys,” Appl. Phys. Lett., 89(19), p. 191911. [CrossRef]
Liang, W. , Zhou, M. , and Ke, F. , 2005, “ Shape Memory Effect in Cu Nanowires,” Nano Lett., 5(10), pp. 2039–2043. [CrossRef] [PubMed]
Liang, W. , and Zhou, M. , 2006, “ Atomistic Simulations Reveal Shape Memory of FCC Metal Nanowires,” Phys. Rev. B, 73(11), p. 115409. [CrossRef]
Li, S. Z. , Ding, X. D. , Li, J. , Ren, X. B. , Sun, J. , and Ma, E. , 2010, “ High-Efficiency Mechanical Energy Storage and Retrieval Using Interfaces in Nanowires,” Nano Lett., 10(5), pp. 1774–1779. [CrossRef] [PubMed]
Wang, P. , Chou, W. , Nie, A. , Huang, Y. , Yao, H. , and Wang, H. , 2011, “ Molecular Dynamics Simulation on Deformation Mechanisms in Body-Centered-Cubic Molybdenum Nanowires,” J. Appl. Phys., 110(9), p. 093521. [CrossRef]
Van Swygenhoven, H. , Derlet, P. M. , and Froseth, A. G. , 2004, “ Stacking Fault Energies and Slip in Nanocrystalline Metals,” Nat. Mater., 3(6), pp. 399–403. [CrossRef] [PubMed]
Kibey, S. , Liu, J. B. , Johnson, D. D. , and Sehitoglu, H. , 2007, “ Energy Pathways and Directionality in Deformation Twinning,” Appl. Phys. Lett., 91(18), p. 181916. [CrossRef]
Bernstein, N. , and Tadmor, E. B. , 2004, “ Tight-Binding Calculations of Stacking Energies and Twinnability in FCC Metals,” Phys. Rev. B, 69(9), p. 094116. [CrossRef]
Rice, J. R. , 1992, “ Dislocation Nucleation From a Crack Tip—An Analysis Based on the Peierls Concept,” J. Mech. Phys. Solids, 40(2), pp. 239–271. [CrossRef]
Shimizu, F. , Ogata, S. , and Li, J. , 2007, “ Theory of Shear Banding in Metallic Glasses and Molecular Dynamics Calculations,” Mater. Trans., 48(11), pp. 2923–2927. [CrossRef]
Gorman, J. A. , Wood, D. S. , and Vreeland, T. , 1969, “ Mobility of Dislocations in Aluminum,” J. Appl. Phys., 40(2), pp. 833–841.
Chang, J. P. , Cai, W. , Bulatov, V. V. , and Yip, S. , 2001, “ Dislocation Motion in BCC Metals by Molecular Dynamics,” Mater. Sci. Eng. A, 309–310, pp. 160–163. [CrossRef]
Berry, J. , Grant, M. , and Elder, K. R. , 2006, “ Diffusive Atomistic Dynamics of Edge Dislocations in Two Dimensions,” Phys. Rev. E, 73(3), p. 031609. [CrossRef]


Grahic Jump Location
Fig. 1

(a) The fitting result to Rose’s equation from the developed meta-atom potential for Cu (refer to Sec. 3.1) and (b) the energy profile along the Bain’s transformation route (inset)

Grahic Jump Location
Fig. 2

(a) The stress–strain curve of a Cu NW upon uniaxial loading at 300 K. Insets are the typical configuration at different plastic deformation stages. (b) Reorientation of Cu NW. The red colored atoms are located on twin boundaries (see color figure online).

Grahic Jump Location
Fig. 10

(a) GPF curves coincides with each other for s = 1.0 and 1.1, (b) the dependence of Peierls stress on the scaling factor, and (c) the relation between the dislocation velocity and the applied shear stress. Details about the simulation and calculation can be found in the online supporting material.

Grahic Jump Location
Fig. 9

Atomic model of an extended dislocation subjected to an applied shear stress

Grahic Jump Location
Fig. 8

Scaling meta-atom model. (a) Equal number of α- and β-type meta-atoms randomly occupy the FCC lattice sites; (b) The spatial distribution of Von Mises strain for the meta-atom Cu-6 wt % Al model after full relaxation. The scaling factor is set to be s = 1.1. The atomic scale Von Mises strain is calculated according to the method in Ref. [34]. (c) The probability density of Von Mises strain for the scaling factor s = 1, 1.05, and 1.1 (symbols) is well fitted by the log-normal distribution (solid lines). (d) The dependence of the mean Von Mises strain and the standard deviation on the scaling factors. (e) The change of the lattice constant and the low-index surface energies with the scaling factors. (f) The change of elastic moduli with the scaling factors.

Grahic Jump Location
Fig. 7

(a) The generalized planar fault curves for meta-atom Cu and a brittle hypothetic FCC metal and (b) Brittle versus ductile response of constructed hypothetic FCC metals

Grahic Jump Location
Fig. 6

Snapshots of the cleavage process by MD simulation of a brittle hypothetic FCC metal: (a) ε = 0.037, (b) ε = 0.039, (c) ε = 0.048, and (d) ε = 0.054

Grahic Jump Location
Fig. 5

Snapshots of dislocation emission events by MD simulation of meta-atom Cu-6 wt % Al: (a) ε = 0.022, (b) ε = 0.025, (c) ε = 0.070, and (d) ε = 0.097

Grahic Jump Location
Fig. 4

Snapshots of dislocation emission events by MD simulation of meta-atom Cu: (a) ε = 0.010, (b) ε = 0.019, (c) ε = 0.053, (d) ε = 0.062, (e) ε = 0.073, and (f) ε = 0.1

Grahic Jump Location
Fig. 3

Meta-atom potentials for aluminum bronzes: (a) the embedding functions F(ρ), (b) the pair interaction functions ϕ(r), (c) the electron density function ρ(r), and (d) the corresponding generalized planar fault curves



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In