Next Article in Journal
Analysis of the Effects of Water Temperature on Water-Assisted Laser Trepanning in Superalloys
Previous Article in Journal
Effect of Material Extrusion Method on the Microstructure and Mechanical Properties of Copper Parts
Previous Article in Special Issue
Mechanical Properties of Interfaces between Mg and SiC: An Ab Initio Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Solute Segregation and Pinning Effect on Lateral Twin Boundary in Magnesium

by
Haoyan Zhang
1,
Qi Zhang
1,
Haowen Sun
1,
Mingyu Gong
1,*,
Jian Wang
2,* and
Yue Liu
1,*
1
State Key Lab of Metal Matrix Composites, School of Materials Science and Engineering, Shanghai Jiao Tong University, Shanghai 200240, China
2
Mechanical and Materials Engineering, University of Nebraska-Lincoln, Lincoln, NE 68588, USA
*
Authors to whom correspondence should be addressed.
Metals 2024, 14(8), 942; https://doi.org/10.3390/met14080942
Submission received: 2 July 2024 / Revised: 10 August 2024 / Accepted: 15 August 2024 / Published: 18 August 2024
(This article belongs to the Special Issue Multi-Scale Simulation of Metals and Alloys)

Abstract

:
Deformation twinning creates a three-dimensional twin domain via the migration of forward, normal and lateral twin boundaries (TBs) with respect to twin shear direction, normal to twin plane and twin lateral direction. Solute segregation and pinning effect on the forward and normal TBs have been experimentally observed and demonstrated via atomistic simulations. Here, we conducted a comprehensive study of solute segregation and the pinning effect on the lateral TBs in Mg. First-principles density functional theory calculations were used to obtain the segregation and formation energies of 19 alloying elements in coherent regions of lateral TBs. Alloying elements with greater difference in atomic radius from Mg generally show more negative segregation energy. Moreover, alloying elements with good solubility are selected to demonstrate the pinning effect on a coherent interface. Ge, Ga, Y, Gd, La and Ca show negative segregation energy and solubility energy, indicating that these elements can form stable segregation and have a strong pinning effect at the lateral boundary. Molecular dynamics simulations revealed that solutes in coherent regions are more effective in pinning lateral TBs than those in misfit regions. The results provide insight into the selection of solute atoms for tailoring twinning behavior.

1. Introduction

The easy activation of deformation twinning in Mg and Mg alloys results in strong plastic anisotropy and is responsible for limited room temperature formability [1,2]. Suppressing or retarding the nucleation and propagation of deformation twins can be achieved through alloy addition [3,4,5], grain refinement [6,7,8] and texture weakening [9,10]. A lot of experiments have demonstrated that alloying solute elements can modify the relative activity of different plastic deformation modes via the synergetic effect of solid-solution hardening on dislocations and solute segregation pinning on twin boundaries (TBs), improving the strength and ductility of Mg and Mg alloys [11]. The pinning effect of solute atoms on TBs has attracted attention recently [12,13], especially the segregation of solute atoms at coherent twin boundaries (CTBs) [14]. In practice, deformation twinning creates a three-dimensional twin domain embedded in a matrix. Corresponding to the crystallography of a twin, twin boundaries are categorized into forward, normal and lateral TBs (Figure 1a), with respect to η1 along twin shear direction, K1 normal to the twin plane, and λ along the twin lateral direction [15]. In addition to normal TB (usually referred to as a coherent twin boundary) which is related to twin thickening, twin propagation in the twinning direction and the lateral direction is accomplished through the migration of the forward and lateral TBs [15].
The segregation of alloying elements at TBs [14,16,17] and their pinning effect on TBs [12] have been revealed by extensive studies, most of which focus on behaviors of forward [17] and normal TBs [12,14,18,19]. In terms of crystallography, forward and normal TBs are tilt boundaries about the λ direction [20] and thus contain excess volume, which is defined as the Voronoi volume difference between the atoms close to the interface and in the matrix. The excess volume V i can be calculated by
V x = V i x V m
where V i x is the Voronoi volume of atoms close to the interface and V m is the Voronoi volume of atoms in the matrix. Solutes with a smaller/larger atomic radius than the matrix generally prefer the TBs’ location with negative/positive excess volume to lower the local stress field [14,21,22]. Atomistic simulations revealed that most of the alloying elements possess negative formation energy [21] and segregation energy [23], suggesting strong segregation along normal and forward TBs and a pinning effect on the migration of the two types of TBs. Meanwhile, lateral TBs are twist boundaries about the λ direction, possessing slight excess volume with respect to the matrix. Therefore, the segregation behaviors of lateral TBs cannot be easily predicted based solely on the excess volume. Based on a recent statistical analysis of electron back-scatter diffraction results [15] and atomistic simulation results [24], deformation twins in Mg propagate faster along the lateral direction than along the forward propagation direction as the twin is smaller than 40 μm during the twin growth. More details about the twin growth can be found in the Supplementary Materials. Therefore, at the beginning of twinning, it is more effective to suppress twin propagation/growth via the pinning migration of lateral TBs, further avoiding materials’ strong plastic anisotropy. In this sense, it is worth understanding solute segregation and the pinning effect on the lateral TBs.
Figure 1. Structural characteristics of lateral TBs. (a) Schematic showing a cube-shaped 3D twin domain. (b1) Plane view of a T-PP2 interface showing a coherent dichromatic complex (CDC) of two prismatic planes in twin (red empty symbols) and matrix (black solid symbols). In the unit hexagonal cells in twin (red rectangle) and matrix (black rectangle), the shadowed region indicates that the two atoms will occupy the same position after atomic shuffling. (b2) The disclination characteristics are revealed in the rotated coherent dichromatic pattern (RCDP). (b3) The coherent dichromatic pattern (CDP) indicates the 90° twist rotation of the T-PP2 interface. The disclination twist angle is φ = 3.71°, reprinted from Ref. [25]. (c1) The atomic structure of a twin domain enclosed by normal TBs (CTB) and lateral TBs (T-PP2), observed along η1 direction. (c2) The atomic structure of a T-PP2 interface observed along the λ direction.
Figure 1. Structural characteristics of lateral TBs. (a) Schematic showing a cube-shaped 3D twin domain. (b1) Plane view of a T-PP2 interface showing a coherent dichromatic complex (CDC) of two prismatic planes in twin (red empty symbols) and matrix (black solid symbols). In the unit hexagonal cells in twin (red rectangle) and matrix (black rectangle), the shadowed region indicates that the two atoms will occupy the same position after atomic shuffling. (b2) The disclination characteristics are revealed in the rotated coherent dichromatic pattern (RCDP). (b3) The coherent dichromatic pattern (CDP) indicates the 90° twist rotation of the T-PP2 interface. The disclination twist angle is φ = 3.71°, reprinted from Ref. [25]. (c1) The atomic structure of a twin domain enclosed by normal TBs (CTB) and lateral TBs (T-PP2), observed along η1 direction. (c2) The atomic structure of a T-PP2 interface observed along the λ direction.
Metals 14 00942 g001
In this work, we conducted first-principles density functional theory (DFT) calculations and molecular dynamics (MD) simulations to investigate solute segregation and the pinning effect on lateral TBs. The atomic structures of the lateral TBs in Mg were, for the first time, simulated and characterized by Liu et al. [25]. As shown in Figure 1(c1,c2), the lateral TBs of a { 1 - 012 } twin exhibit a typical semi-coherent twist boundary with misfit dislocations [25]. The coherent interface is parallel to prismatic planes (referred to as twist prismatic/prismatic (T-PP2) interface) in both the twin and matrix. The formation of the misfit dislocations cancels the φ = 3.71° twist disclination angle (Figure 1(b1–b3)) associated with the pile up of TDs on the T-PP2 interface. In the research [25], Y. Liu, N. Li, S. Shao, et al. used high-resolution transmission electron microscopy and atomic simulations to characterize the lateral boundary of the {10 1 ¯ 2} deformation twins in magnesium. They found that the structure of the lateral boundary is a semi-coherent T-PP2 boundary. In the research [26], J.P. Hirth, J. Wang and C.N. Tome used topological model, simulation and characterization methods to analyze the structure and defects of twin interfaces. They also demonstrated that the structure of the interface along the η direction contained screw dislocation. DFT calculations revealed that the coherent T-PP2 interface possesses slight excess volume but obvious bonding distortion and enabled to obtain the segregation and formation energies of 19 alloying elements in the coherent T-PP2 interface. Alloying elements with a greater difference in atomic radius to Mg generally show more negative segregation energy. Moreover, alloying elements with good solubility are selected to demonstrate the pinning effect on coherent T-PP2 interfaces. MD simulations revealed that solutes in coherent regions are more effective in pinning lateral TBs than those in misfit regions. The results provide an insight into the selection of solute atoms for tailoring twin response.

2. Methodology

The DFT calculations are performed with the Vienna Ab-initio Simulation Package (VASP) [27,28] to investigate the solubility and pinning effect of different alloying elements on a coherent region of lateral TBs. The generalized gradient approximation (GGA) with Perdew–Burke–Ernzerhof (PBE) parametrization [29,30] is used for the exchange and correlation functions. The core electrons are replaced by the projector augmented wave (PAW) pseudopotentials [31] with the valence states shown in Table 1. The plane-wave energy cut-off of 500 eV and a Methfessel–Paxton smearing of 0.01 eV are used. For all DFT calculations, the self-consistent iteration is stopped when the change in total energy is smaller than 10−5 eV. The convergence criterion of the geometry optimizations is that the forces acting on each atom are smaller than 0.01 eV/Å. The 80-atom 0.54 × 6.39 × 0.54 nm3 models containing coherent T-PP2 interfaces with/without a solute element (X) are constructed in the coordinates where the x-, y- and z-axes are along the [ 1 2 - 10 ]m||[ 1 2 - 10 ]t, [ 0001 ]m||[ 1 - 010 ]t and [ 1 - 010 ]m||[ 000 1 - ]t directions. With a 7 × 1 × 7 Monkhorst–Pack (M-P) K-point grid, all the models are relaxed.
The simulations of the migration of the T-PP2 interfaces, both in the absence and presence of solute segregation, are conducted using a bi-crystal model of the dimensions of 247.69 × 28.51 × 0.76 nm3 with a ( 1 - 012 ) twinning relationship in the coordinates where the x-, y- and z-axes are along the λ, K1 and η1 directions. This orientation ensures that the simulated interfaces accurately mimic the structural characteristics and migration behavior of real-world T-PP2 interfaces. The modified embedded atom method (MEAM) potential for the Mg-Ca system developed by Kim et al. [32] is employed in a Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) [33]. In the simulation, periodic boundary conditions are applied along the y- and z-directions. Additionally, a fixed 1 nm wide layer is introduced perpendicular to the x-direction. Prior to simulating the migration process, the model is subjected to a relaxation process to bring it to a near-equilibrium state. This relaxation is performed at a temperature of 10 K for 20 ps in order to obtain the most stable structure and allow for a more controlled investigation of the migration behavior. During the relaxation process, the atomic positions are allowed to adjust until the system reaches a stable configuration, characterized by minimal changes in energy and atomic positions.

3. Results and Discussion

Figure 2a manifests the atomic structure of a coherent T-PP2 interface. Although atoms near the interface have negligible excess volume ( V B : 0.138, V C : 0.2343, V D : 0.2342), the bonds between them and their neighboring atoms are significantly distorted. Figure 2(b1–b4) displays the twelve first nearest neighbor (1NN) bonds of atoms A-D; these are labeled in Figure 2a. The bond-i (i = 1, 2, …, 12) associated with atom-x (x = A, B, C, D) can be defined. Figure 2c compares the lengths of bond-i in Figure 2(b1–b4). All bonds of atom A that are far from the interface possess a 3.1~3.3 Å length. For the atoms B-D that are near the interface, their bond-7 is significantly extended; these measure 3.53, 3.87 and 3.85 Å. The length of the other bonds is still within 3.1~3.3 Å. Generally, bond strength becomes weaker when its length is larger than its equilibrium [34]. According to the research by Zhao, the type and length can significantly affect the elastic stiffness of the metallic glasses [35,36]. So, one can reasonably speculate that the region near the coherent T-PP2 interface is less stiff than the single crystal.
The segregation energy of a solute ( E s ), which is the energy difference between solutes on and far from the interface, evaluates the pinning effect, and is calculated by
E s = E tw Mg 79 X 1 C E tw Mg 79 X 1 A
where E tw Mg 79 X 1 C and E tw Mg 79 X 1 A are the system energies of the models with solutes X substituting atoms C and A, defined in Figure 2a. The segregation energy for 19 alloying elements is plotted in Figure 3a and summarized in Table 2. The 19 elements studied in reference [21] are selected to compare the segregation behavior between lateral TBs and normal/forward TBs. Most of the 19 elements are used to prepare the corresponding alloys in experiments [37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53]. It is found that the segregation energy becomes more negative with an increasing difference in the atomic radius between Mg (RMg) and solute X (RX), no matter whether RX is smaller or larger than RMg. This is different from the solute segregation behaviors on normal and forward TBs, on which a smaller/larger RX than RMg only results in negative segregation energy at a location with negative/positive excess volume to lower the associated elastic energy [21]. The unique solute segregation behaviors on lateral TB can still be explained from a mechanical point of view. Solutes with RX different from RMg prefer to locate in less stiff regions near the interface to lower the associated elastic energy [14]. It should be noted that Ce solute does not follow the trend because substituting Ce in a single crystal does not trigger an obvious strain field (Figure 3b).
In addition to negative segregation energy, ideal alloying elements that inhibit the migration of the coherent T-PP2 interface should have great solubility on the interface. The formation energy of a solute ( E f ) can evaluate the solubility, and is estimated by
E f = E tw Mg 79 X 1 C E tw Mg 80 E coh X + E coh Mg
where E tw Mg 80 is the system energy of the models without a solute on the interface. E coh X and E coh Mg are cohesive energies of bulk materials, as listed in Table 1. The formation energy with respect to the alloying element is summarized in Table 2. Ga, Cd, Pm and La solutes possess the highest solubility, while Ge, Li, Dy, Y, Gd and Ca solutes have a relatively lower solubility. The solubility energy associated with Zn, Al, Sc and Tb solutes is nearly zero and is positive when associated with the Zr, Na, K, Ce and Ba solutes. As shown in Figure 4, Ge, Ga, Y, Gd, Pm, La and Ca solutes with more negative formation and segregation energies have relatively greater solubility and can subsequently pin a coherent T-PP2 interface. Zr and Ce solutes with more positive formation and segregation energies should be less effective. Generally, when this method is used to measure the segregation behavior of solutes, the lower the segregation energy and the formation energy are, the more effective the segregation of solute atoms can be formed; that is, the solute can form a stable segregation and strong pinning effect on the interface. It is worth noting that this method is not only applicable to the T-PP2 interface, but can also be used to evaluate the segregation of solute atoms on any interfaces. The results may explain previously made experimental observations such as Y being more effective in reducing twin activity than Al [54], La addition significantly reducing the twin fraction in Mg-Ce alloys [55] and Gd being more effective than Al and Zn in delaying twinning [56]. Meanwhile, some numerical results also show that Y segregation can significantly reduce the twin volume fraction [57], while Gd segregation can hinder the motion of a twin tip and delay the propagation of a deformation twin [58].
The relaxed T-PP2 interface as shown in Figure 5a is semi-coherent. Misfit dislocations with the Burgers vector bmis = 0.47[ 1 - 01 1 - ] periodically locate on the interface [59]. The distance between two nearby misfit dislocations is 3.12 nm. The migration process of the T-PP2 interface without Ca segregation under σ13 = 1100 MPa at 300 K is shown in Figure 5(b1–b3). It should be noted that the high shear stress is applied to enable the migration of lateral TBs within the time-scale (~nano-second) of MD simulations. So, the results from MD simulations are only capable to show the trend in solute segregation along boundaries. At the beginning (10 ps, Figure 5b1), the coherent regions move leftward while the misfit dislocations remain steady. Then, the misfit dislocations can either move leftward (28 ps, Figure 5b2) or rightward (32 ps, Figure 5b3). The leftward motion of misfit dislocations is activated via the dragging of the migrating coherent regions, while the rightward motion is driven using the Peach–Koehler force [60] under external load. A rapid migration of the T-PP2 interface is observed after the depinning of the coherent regions from the misfit dislocations. As shown in Figure 5c, the T-PP2 interface without a Ca solute migrates 29.64 nm in 40 ps. Figure 5d,e illustrate the difference in the pinning effect between Ca solutes (0.16 atom/nm2 areal density) in coherent and misfit regions. The Ca solutes in misfit regions exhibit a slight pinning effect. The T-PP2 interface migrates 16.41 nm in 40 ps under the same loading condition (Figure 5d). Meanwhile, the Ca solutes in coherent regions strongly pin the migration of the T-PP2 interface because they significantly retard the motion of the coherent region that initiates the migration of the entire T-PP2 interface. Consequently, the T-PP2 interface only migrates 2.11 nm in 40 ps (Figure 5e). The MD results suggest that solutes in coherent regions are more effective in pinning lateral TBs than those in misfit regions. In this sense, alloying elements that can strongly segregate in and obstruct coherent regions of lateral TBs are desired for the inhibition of the lateral growth of deformation twins.

4. Conclusions

Deformation twins, as three-dimensional structural entities, are bounded by three kinds of interfaces, encompassing normal, forward, and lateral TBs. These boundaries play pivotal roles in governing the propagation and evolution of twinning. Among these, the rapidly migrating lateral TBs have emerged as critical factors influencing the growth of deformation twins. By effectively pinning the fast-moving lateral TBs, it becomes feasible to significantly impede the twin growth, thereby modulating the microstructural evolution and mechanical properties of the material. In this work, solute segregation and its pinning effect on lateral TBs (T-PP2 interface) have been investigated by MD modeling and DFT calculations. The results provide an insight into the selection of solute atoms for controlling twin activity.
  • The coherent T-PP2 interface of the lateral TBs with minor excess volume and obvious extended bonds should be less stiff than the region far from the interface. Therefore, alloying elements with a significant difference in atomic radius with regard to Mg possess more negative segregation energy.
  • Combined with the data of formation energy, DFT calculations enable the selection of alloying elements with a good solubility and pinning effect on a coherent T-PP2 interface. Ge, Ga, Y, Gd, Pm, La and Ca can form stable segregation on the T-PP2; meanwhile, they also have a strong pinning effect on the interface.
  • Along the semi-coherent T-PP2 interface, solute atoms segregated in coherent regions are more effective at pinning TB migration than those segregated in misfit regions.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/met14080942/s1, the supplementary information shows the analysis of three-dimensional deformation twin growth in pure Mg. Figures S1–S3: Growth of 3D twins under different loading stress at 100 K, 300 K, 500 K. Figure S4: Statistical graphs of growth processes of three-dimensional twin at different temperatures. Table S1: The dimensions of the deformation twin at 8 ps.

Author Contributions

Conceptualization, M.G. and J.W.; methodology, Q.Z.; validation, H.S.; formal analysis, H.Z. and M.G.; investigation, H.Z.; resources, J.W.; data curation, H.Z.; writing—original draft, H.Z.; writing—review and editing, Q.Z., H.S., M.G., Y.L. and J.W.; supervision, M.G., Y.L. and J.W.; funding acquisition, M.G. All authors have read and agreed to the published version of the manuscript.

Funding

The authors sincerely acknowledge the financial supports by the National Key Research and Development Program of China (No. 2022YFB3708900).

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors on request.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Koike, J.; Kobayashi, T.; Mukai, T.; Watanabe, H.; Suzuki, M.; Maruyama, K.; Higashi, K. The activity of non-basal slip systems and dynamic recovery at room temperature in fine-grained AZ31B magnesium alloys. Acta Mater. 2003, 51, 2055–2065. [Google Scholar] [CrossRef]
  2. Arul Kumar, M.; Beyerlein, I.J.; Lebensohn, R.A.; Tomé, C.N. Role of alloying elements on twin growth and twin transmission in magnesium alloys. Mater. Sci. Eng. A 2017, 706, 295–303. [Google Scholar] [CrossRef]
  3. Bian, M.Z.; Huang, X.; Chino, Y. Substantial improvement in cold formability of concentrated Mg–Al–Zn–Ca alloy sheets by high temperature final rolling. Acta Mater. 2021, 220, 117328. [Google Scholar] [CrossRef]
  4. Ghazisaeidi, M.; Hector, L.G.; Curtin, W.A. Solute strengthening of twinning dislocations in Mg alloys. Acta Mater. 2014, 80, 278–287. [Google Scholar] [CrossRef]
  5. Maghsoudi, M.; Ovri, H.; Krekeler, T.; Ziehmer, M.; Lilleodden, E.T. Effect of Gd solutes on the micromechanical response of twinning and detwinning in Mg. Acta Mater. 2023, 258, 119202. [Google Scholar] [CrossRef]
  6. Song, E.; Andani, M.T.; Misra, A. Investigation of grain size and geometrically necessary dislocation density dependence of flow stress in Mg-4Al by using nanoindentation. Acta Mater. 2024, 265, 119633. [Google Scholar] [CrossRef]
  7. Buha, J. Grain refinement and improved age hardening of Mg–Zn alloy by a trace amount of V. Acta Mater. 2008, 56, 3533–3542. [Google Scholar] [CrossRef]
  8. Wen, Y.; Guan, B.; Xin, Y.; Liu, C.; Wu, P.; Huang, G.; Liu, Q. Solute atom mediated Hall-Petch relations for magnesium binary alloys. Scr. Mater. 2022, 210, 114451. [Google Scholar] [CrossRef]
  9. Wu, J.; Jin, L.; Dong, J.; Wang, F.; Dong, S. The texture and its optimization in magnesium alloy. J. Mater. Sci. Technol. 2020, 42, 175–189. [Google Scholar] [CrossRef]
  10. Pei, R.; Zou, Y.; Zubair, M.; Wei, D.; Al-Samman, T. Synergistic effect of Y and Ca addition on the texture modification in AZ31B magnesium alloy. Acta Mater. 2022, 233, 117990. [Google Scholar] [CrossRef]
  11. Wang, J.Y.; Li, N.; Alizadeh, R.; Monclús, M.A.; Cui, Y.W.; Molina-Aldareguía, J.M.; Llorca, J. Effect of solute content and temperature on the deformation mechanisms and critical resolved shear stress in Mg-Al and Mg-Zn alloys. Acta Mater. 2019, 170, 155–165. [Google Scholar] [CrossRef]
  12. Li, Z.; Chen, H. Effects of solute segregation on the mobility and mechanical properties of {10 1 - 2} coherent twin boundary in magnesium. Comput. Mater. Sci. 2023, 227, 112297. [Google Scholar]
  13. Jiang, X.; Xi, C.; Liu, X.; Wang, L.; Kang, Y.; Huo, Q. Accelerating and decelerating effects of Gd content and peak-aging treatment on the grain growth of magnesium. Mater. Today Commun. 2023, 37, 107581. [Google Scholar] [CrossRef]
  14. Nie, J.F.; Zhu, Y.M.; Liu, J.Z.; Fang, X.Y. Periodic Segregation of Solute Atoms in Fully Coherent Twin Boundaries. Science 2013, 340, 957–960. [Google Scholar] [CrossRef]
  15. Liu, Y.; Tang, P.Z.; Gong, M.Y.; McCabe, R.J.; Wang, J.; Tomé, C.N. Three-dimensional character of the deformation twin in magnesium. Nat. Commun. 2019, 10, 3308. [Google Scholar] [CrossRef]
  16. Zhu, Y.M.; Xu, S.W.; Nie, J.F. {10 1 - 1} Twin boundary structures in a Mg–Gd alloy. Acta Mater. 2018, 143, 1–12. [Google Scholar]
  17. Wang, X.; Hu, Y.; Yu, K.; Mahajan, S.; Beyerlein, I.J.; Lavernia, E.J.; Rupert, T.J.; Schoenung, J.M. Room Temperature Deformation-induced Solute Segregation and its Impact on Twin Boundary Mobility in a Mg-Y Alloy. Scr. Mater. 2022, 209, 114375. [Google Scholar] [CrossRef]
  18. Tsuru, T.; Somekawa, H.; Chrzan, D.C. Interfacial segregation and fracture in Mg-based binary alloys: Experimental and first-principles perspective. Acta Mater. 2018, 151, 78–86. [Google Scholar] [CrossRef]
  19. He, C.; Li, Z.; Chen, H.; Wilson, N.; Nie, J.F. Unusual solute segregation phenomenon in coherent twin boundaries. Nat. Commun. 2021, 12, 722. [Google Scholar] [CrossRef]
  20. Wang, J.; Beyerlein, I.J. Atomic structures of symmetric tilt grain boundaries in hexagonal close packed (hcp) crystals. Model. Simul. Mater. Sci. Eng. 2012, 20, 024002. [Google Scholar] [CrossRef]
  21. Kumar, A.; Wang, J.; Tomé, C.N. First-principles study of energy and atomic solubility of twinning-associated boundaries in hexagonal metals. Acta Mater. 2015, 85, 144–154. [Google Scholar] [CrossRef]
  22. Nagarajan, D. Quantitative measurements of strain-induced twinning in Mg, Mg-Al and Mg-Zn solid solutions. Mater. Sci. Eng. A 2022, 860, 144292. [Google Scholar] [CrossRef]
  23. Huang, Z.; Nie, J.F. Interaction between hydrogen and solute atoms in {10 1 - 2} twin boundary and its impact on boundary cohesion in magnesium. Acta Mater. 2021, 214, 117009. [Google Scholar]
  24. Gong, M.Y.; Graham, J.; Taupin, V.; Capolungo, L. The effects of stress, temperature and facet structure on growth of {10 1 - 2} twins in Mg: A molecular dynamics and phase field study. Acta Mater. 2021, 208, 116603. [Google Scholar]
  25. Liu, Y.; Li, N.; Shao, S.; Gong, M.Y.; Wang, J.; McCabe, R.J.; Jiang, Y.; Tomé, C.N. Characterizing the boundary lateral to the shear direction of deformation twins in magnesium. Nat. Commun. 2016, 7, 11577. [Google Scholar] [CrossRef] [PubMed]
  26. Hirth, J.P.; Wang, J.; Tomé, C.N. Disconnections and other defects associated with twin interfaces. Prog. Mater. Sci. 2016, 83, 417–471. [Google Scholar] [CrossRef]
  27. Kresse, G.; Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metal–amorphous-semiconductor transition in germanium. Phys. Rev. B 1994, 49, 14251–14269. [Google Scholar] [CrossRef] [PubMed]
  28. Kresse, G.; Furthmuller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 1996, 54, 11169–11186. [Google Scholar] [CrossRef] [PubMed]
  29. Perdew, J.P.; Yue, W. Accurate and simple density functional for the electronic exchange energy: Generalized gradient approximation. Phys. Rev. B 1986, 33, 8800–8802. [Google Scholar] [CrossRef]
  30. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865–3868. [Google Scholar] [CrossRef]
  31. Blöchl, P.E. Projector augmented-wave method. Phys. Rev. B 1994, 50, 17953–17979. [Google Scholar] [CrossRef]
  32. Kim, K.H.; Jeon, J.B.; Lee, B.J. Modified embedded-atom method interatomic potentials for Mg–X (X = Y, Sn, Ca) binary systems. Calphad 2015, 48, 27–34. [Google Scholar] [CrossRef]
  33. Plimpton, S. Fast Parallel Algorithms for Short-Range Molecular Dynamics. J. Comput. Phys. 1995, 117, 1–19. [Google Scholar] [CrossRef]
  34. Morse, P.M. Diatomic Molecules According to the Wave Mechanics. II. Vibrational Levels. Phys. Rev. 1929, 34, 57–64. [Google Scholar] [CrossRef]
  35. Zhao, W.; Cheng, J.L.; Feng, S.D.; Li, G.; Liu, R.P. Intrinsic correlation between elastic modulus and atomic bond stiffness in metallic glasses. Mater. Lett. 2016, 175, 227–230. [Google Scholar] [CrossRef]
  36. Zhao, W.; Li, G.; Wang, Y.Y.; Feng, S.D.; Qi, L.; Ma, M.Z.; Liu, R.P. Atomic bond proportions and relations to physical properties in metallic glasses. Mater. Des. 2015, 65, 1048–1052. [Google Scholar] [CrossRef]
  37. Chen, X.; Venezuela, J.; Dargusch, M. The high corrosion-resistance of ultra-high purity Mg–Ge alloy and its discharge performance as anode for Mg–air battery. Electrochim. Acta 2023, 448, 142127. [Google Scholar] [CrossRef]
  38. Dai, Y.; Xiao, M.; Hu, Y.; Yang, Y.; Jiang, B.; Zheng, T.; Dong, L.; Yang, B.; Zheng, C. Effect of extrusion parameters on microstructure and mechanical properties of Mg-0.12Ca-0.08Ba alloy. J. Mater. Res. Technol. 2022, 20, 1570–1579. [Google Scholar] [CrossRef]
  39. Drozdenko, D.; Dobroň, P.; Yi, S.; Horváth, K.; Letzig, D.; Bohlen, J. Mobility of pinned twin boundaries during mechanical loading of extruded binary Mg-1Zn alloy. Mater. Charact. 2018, 139, 81–88. [Google Scholar] [CrossRef]
  40. Franke, P.; Neuschütz, D. Mg-Yb (Magnesium-Ytterbium). In Binary Systems. Part 5: Binary Systems Supplement 1: Phase Diagrams, Phase Transition Data, Integral and Partial Quantities of Alloys; Franke, P., Neuschütz, D., Eds.; Springer: Berlin/Heidelberg, Germany, 2007; pp. 1–3. [Google Scholar]
  41. Ge, B.; Yang, M.; Zu, Q.; Guo, J.; Tian, Y.; Peng, Q. Lithium Cluster Segregation in Coherent Contraction Twin Boundaries of Magnesium Alloy. Acta Mater. 2020, 201, 477–487. [Google Scholar] [CrossRef]
  42. Hanna, A.; Azzeddine, H.; Huang, Y.; Bradai, D.; Cabrera, J.M.; Langdon, T.G. An investigation of the thermal stability of an MgDy alloy after processing by high-pressure torsion. Mater. Charact. 2019, 151, 519–529. [Google Scholar] [CrossRef]
  43. Hui, W.; Sun, L.; Zhou, Y.; Gao, Y.; Zhai, W.; Dong, H.; Liu, Y.; Gao, Q.; Peng, J. Effect of Mg–Ga alloy intermetallic compounds on the properties of soluble aluminum alloy. Intermetallics 2022, 144, 107501. [Google Scholar] [CrossRef]
  44. Li, S.; Jin, Z.; Gong, X.; Huang, H.; Guan, S.; Yuan, G. Chemical ordering effects on martensitic transformations in Mg-Sc alloys. Acta Mater. 2023, 249, 118854. [Google Scholar] [CrossRef]
  45. Ma, B.; Tan, C.; Ouyang, L.; Shao, H.; Wang, N.; Zhu, M. Microstructure and discharge performance of Mg-La alloys as the anodes for primary magnesium-air batteries. J. Alloys Compd. 2022, 918, 165803. [Google Scholar] [CrossRef]
  46. Ma, X.; Jiao, Q.; Kecskes, L.J.; El-Awady, J.A.; Weihs, T.P. Effect of basal precipitates on extension twinning and pyramidal slip: A micro-mechanical and electron microscopy study of a Mg–Al binary alloy. Acta Mater. 2020, 189, 35–46. [Google Scholar] [CrossRef]
  47. Mishra, R.K.; Brahme, A.; Sabat, R.K.; Jin, L.; Inal, K. Twinning and texture randomization in Mg and Mg-Ce alloys. Int. J. Plast. 2019, 117, 157–172. [Google Scholar] [CrossRef]
  48. Molodov, K.D.; Al-Samman, T.; Molodov, D.A. Effect of gadolinium on the deformation and recrystallization behavior of magnesium crystals. Acta Mater. 2022, 240, 118312. [Google Scholar] [CrossRef]
  49. Pan, H.; Yang, C.; Yang, Y.; Dai, Y.; Zhou, D.; Chai, L.; Huang, Q.; Yang, Q.; Liu, S.; Ren, Y.; et al. Ultra-fine grain size and exceptionally high strength in dilute Mg–Ca alloys achieved by conventional one-step extrusion. Mater. Lett. 2019, 237, 65–68. [Google Scholar] [CrossRef]
  50. Gao, S.; Ren, S.; Liu, C.; Wu, Z.; Li, H. Effect of Cd Addition on Microstructure and Properties of Mg-Cd Binary Magnesium Alloy. Rare Met. Mater. Eng. 2015, 44, 2401–2404. [Google Scholar]
  51. Tan, X.F.; Kim, M.; Gu, Q.; Piraquive Pinzon, J.; Zeng, G.; McDonald, S.D.; Nogita, K. Na-modified cast hypo-eutectic Mg–Mg2Si alloys for solid-state hydrogen storage. J. Power Sources 2022, 538, 231538. [Google Scholar] [CrossRef]
  52. Tong, X.; Wu, G.; Easton, M.A.; Sun, M.; StJohn, D.H.; Jiang, R.; Qi, F. Exceptional grain refinement of Mg-Zr master alloy treated by tungsten inert gas arc re-melting with ultra-high frequency pulses. Scr. Mater. 2022, 215, 114700. [Google Scholar] [CrossRef]
  53. Zhang, Y.; Zhang, Z.; Nagaumi, H.; Yang, X. Study on the mechanism of dilute Mn addition in enhancing the creep resistance of hot-rolled Mg-1.5Y alloy sheet. Mat. Sci. Eng. A-Struct. 2023, 873, 145026. [Google Scholar] [CrossRef]
  54. Shi, Q.; Allison, J. An Experimental Study on Twinning Behavior in Mg Alloys with Different Solute Elements. In Magnesium Technology; Leonard, A., Barela, S., Neelameggham, N.R., Miller, V.M., Tolnai, D., Eds.; Springer Nature: Cham, Switzerland, 2024; pp. 69–71. [Google Scholar]
  55. Yang, W.; Tekumalla, S.; Gupta, M. Cumulative Effect of Strength Enhancer—Lanthanum and Ductility Enhancer—Cerium on Mechanical Response of Magnesium. Metals 2017, 7, 241. [Google Scholar] [CrossRef]
  56. Nagarajan, D.; Cáceres, C.H. The Friction Stress of the Hall–Petch Relationship of Pure Mg and Solid Solutions of Al, Zn, and Gd. Metall. Mater. Trans. A 2018, 49, 5288–5297. [Google Scholar] [CrossRef]
  57. Galido-Nava, E.I. Modelling twinning evolution during plastic deformation in hexagonal close-packed metals. Mater. Des. 2015, 83, 327–343. [Google Scholar] [CrossRef]
  58. Kristián, M.; Andrea, F.; Michal, K.; Andriy, O.; Gergely, F.; Gergely, N.; Petr, H.; Daria, D. The influence of gadolinium concentration on the twin propagation rate in magnesium alloys. J. Alloys Compd. 2023, 948, 169635. [Google Scholar]
  59. Gong, M.Y.; Hirth, J.P.; Liu, Y.; Shen, Y.; Wang, J. Interface structures and twinning mechanisms of twins in hexagonal metals. Mater. Res. Lett. 2017, 5, 449–464. [Google Scholar] [CrossRef]
  60. Peach, M.; Koehler, J.S. The Forces Exerted on Dislocations and the Stress Fields Produced by Them. Phys. Rev. 1950, 80, 436–439. [Google Scholar] [CrossRef]
Figure 2. Structural characteristics of a coherent T-PP2 interface. (a) The atomic configuration of the interface with coloring showing the atomic volume. The 1NN bonds of atoms (b1) far from and (b2b4) near to the interface. (c) A comparison of the length of the 1NN bonds shown in (b1b4).
Figure 2. Structural characteristics of a coherent T-PP2 interface. (a) The atomic configuration of the interface with coloring showing the atomic volume. The 1NN bonds of atoms (b1) far from and (b2b4) near to the interface. (c) A comparison of the length of the 1NN bonds shown in (b1b4).
Metals 14 00942 g002
Figure 3. Solute segregation behaviors on coherent T-PP2 interface. (a) Segregation energy of solutes on interface and (b) averaged change in bond length by solutes in single crystal with respect to alloying elements. (The horizontal dotted line represents 0 segregation energy in (a) which means that there is no extra pinning effect, and 0 change in bond length in (b), which means that there is no change in bond length. The vertical dotted line represents Mg in both (a) and (b), the radius of the elements on the left of the vertical dotted line are smaller than Mg, while the radius of the elements on the right of the vertical dotted line are larger than Mg.)
Figure 3. Solute segregation behaviors on coherent T-PP2 interface. (a) Segregation energy of solutes on interface and (b) averaged change in bond length by solutes in single crystal with respect to alloying elements. (The horizontal dotted line represents 0 segregation energy in (a) which means that there is no extra pinning effect, and 0 change in bond length in (b), which means that there is no change in bond length. The vertical dotted line represents Mg in both (a) and (b), the radius of the elements on the left of the vertical dotted line are smaller than Mg, while the radius of the elements on the right of the vertical dotted line are larger than Mg.)
Metals 14 00942 g003
Figure 4. Selection of alloying elements to inhibit migration of coherent T-PP2 interface based on formation and segregation energies data. The elements in the red circle have strong pinning effect and greater solubility on the T-PP2 interface.
Figure 4. Selection of alloying elements to inhibit migration of coherent T-PP2 interface based on formation and segregation energies data. The elements in the red circle have strong pinning effect and greater solubility on the T-PP2 interface.
Metals 14 00942 g004
Figure 5. Migration of lateral TBs with/without solutes. (a) Morphology of semi-coherent T-PP2 interface with misfit dislocations periodically located on coherent boundaries. Snapshots of migrating T-PP2 interface without Ca solute at (b1) 10 ps, (b2) 28 ps and (b3) 32 ps. Snapshots at 40 ps of T-PP2 interface (c) without Ca solute, (d) with Ca solute in misfit region and (e) with Ca solute in coherent region.
Figure 5. Migration of lateral TBs with/without solutes. (a) Morphology of semi-coherent T-PP2 interface with misfit dislocations periodically located on coherent boundaries. Snapshots of migrating T-PP2 interface without Ca solute at (b1) 10 ps, (b2) 28 ps and (b3) 32 ps. Snapshots at 40 ps of T-PP2 interface (c) without Ca solute, (d) with Ca solute in misfit region and (e) with Ca solute in coherent region.
Metals 14 00942 g005
Table 1. Details of first-principles calculations on bulk materials.
Table 1. Details of first-principles calculations on bulk materials.
ElementValencyK-MeshStructureLattice
(Å)
ECoh
(eV)
Zn1230 × 29 × 27Hexagonala = 2.61, c = 4.87−1.10
Ge147 × 12 × 23Cubica = 5.67−4.51
Ga139 × 9 × 29Orthorhombica = 4.43,
b = 7.60, c = 4.56
−2.89
Al37 × 12 × 23Cubica = 4.04−3.75
Cd1227 × 25 × 22Hexagonala = 2.96, c = 5.89−0.74
Li39 × 15 × 30Cubica = 4.33−1.91
Mg225 × 24 × 26Hexagonala = 3.17, c = 5.14−1.51
Zr1224 × 23 × 25Hexagonala = 3.24, c = 5.17−8.52
Sc119 × 14 × 28Cubica = 4.65−6.20
Na710 × 16 × 31Cubica = 4.21−1.31
Dy922 × 21 × 24Hexagonala = 3.64, c = 5.58−4.53
Y118 × 13 × 26Cubica = 5.10−6.40
Yb810 × 15 × 31Cubica = 4.24−1.46
Gd922 × 21 × 22Hexagonala = 3.64, c = 5.59−4.57
Pm1122 × 21 × 11Hexagonala = 3.68, c = 11.80−4.68
Ca107 × 12 × 24Cubica = 5.58−1.93
K98 × 12 × 24Cubica = 5.40−1.04
Ce129 × 14 × 28Cubica = 4.67−5.92
La1121 × 20 × 11Hexagonala = 3.80, c = 12.04−4.88
Ba109 × 8 × 18Hexagonala = 4.63,
b = 8.02, c = 7.29
−1.88
Table 2. Segregation and formation energies of different alloying elements on coherent T-PP2 interface.
Table 2. Segregation and formation energies of different alloying elements on coherent T-PP2 interface.
Solute E s (eV) E f (eV)
Zn−0.055−0.082
Ge−0.047−0.230
Ga−0.045−0.281
Al−0.0130.095
Cd−0.008−0.290
Li−0.031−0.177
Mg--
Zr0.0340.611
Sc−0.014−0.044
Na−0.0430.302
Dy−0.039−0.129
Y−0.048−0.224
Yb−0.064−0.032
Gd−0.043−0.197
Pm−0.061−0.259
Ca−0.074−0.127
K−0.1291.106
Ce0.0490.832
La−0.087−0.281
Ba−0.2120.256
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, H.; Zhang, Q.; Sun, H.; Gong, M.; Wang, J.; Liu, Y. Solute Segregation and Pinning Effect on Lateral Twin Boundary in Magnesium. Metals 2024, 14, 942. https://doi.org/10.3390/met14080942

AMA Style

Zhang H, Zhang Q, Sun H, Gong M, Wang J, Liu Y. Solute Segregation and Pinning Effect on Lateral Twin Boundary in Magnesium. Metals. 2024; 14(8):942. https://doi.org/10.3390/met14080942

Chicago/Turabian Style

Zhang, Haoyan, Qi Zhang, Haowen Sun, Mingyu Gong, Jian Wang, and Yue Liu. 2024. "Solute Segregation and Pinning Effect on Lateral Twin Boundary in Magnesium" Metals 14, no. 8: 942. https://doi.org/10.3390/met14080942

APA Style

Zhang, H., Zhang, Q., Sun, H., Gong, M., Wang, J., & Liu, Y. (2024). Solute Segregation and Pinning Effect on Lateral Twin Boundary in Magnesium. Metals, 14(8), 942. https://doi.org/10.3390/met14080942

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop