Next Article in Journal
The Cone Beam Computed Tomography Evaluation of Cortical Bone Plate after Piezocision-Assisted Orthodontic Upper Arch Expansion: A Case Series
Next Article in Special Issue
Facile Synthesis of BiVO4@ZIF−8 Composite with Heterojunction Structure for Photocatalytic Wastewater Treatment
Previous Article in Journal
Correction: Poljak, M.; Matić, M. Metallization-Induced Quantum Limits of Contact Resistance in Graphene Nanoribbons with One-Dimensional Contacts. Materials 2021, 14, 3670
Previous Article in Special Issue
Multi-Scale Structure–Mechanical Property Relations of Graphene-Based Layer Materials
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of Mo Segregation at Grain Boundaries on the High Temperature Creep Behavior of Ni-Mo Alloys: An Atomistic Study

1
School of Science, North University of China, Taiyuan 030051, China
2
International Research Center for Computational Mechanics, State Key Laboratory of Structural Analysis for Industrial Equipment, Department of Engineering Mechanics, Faculty of Vehicle Engineering and Mechanics, Dalian University of Technology, Dalian 116024, China
3
Department of Physics and Astronomy, Clemson University, Clemson, SC 29634, USA
*
Author to whom correspondence should be addressed.
Materials 2021, 14(22), 6966; https://doi.org/10.3390/ma14226966
Submission received: 20 October 2021 / Revised: 11 November 2021 / Accepted: 15 November 2021 / Published: 18 November 2021

Abstract

:
Based on molecular dynamics simulations, the creep behaviors of nanocrystalline Ni before and after the segregation of Mo atoms at grain boundaries are comparatively investigated with the influences of external stress, grain size, temperature, and the concentration of Mo atoms taken into consideration. The results show that the creep strain rate of nanocrystalline Ni decreases significantly after the segregation of Mo atoms at grain boundaries due to the increase of the activation energy. The creep mechanisms corresponding to low, medium, and high stress states are respectively diffusion, grain boundary slip and dislocation activities based on the analysis of stress exponent and grain size exponent for both pure Ni and segregated Ni-Mo samples. Importantly, the influence of external stress and grain size on the creep strain rate of segregated Ni-Mo samples agrees well with the classical Bird-Dorn-Mukherjee model. The results also show that segregation has little effect on the creep process dominated by lattice diffusion. However, it can effectively reduce the strain rate of the creep deformation dominated by grain boundary behaviors and dislocation activities, where the creep rate decreases when increasing the concentration of Mo atoms at grain boundaries within a certain range.

1. Introduction

In the past few decades, nanocrystalline (NC) materials have attracted more and more attention because of their excellent properties such as ultra-high strength at room temperature [1,2,3,4,5,6,7]. However, when functioned at high temperatures and under continuous stresses, creep deformation, which may cause accidental deformation or even failure of materials and structures [8,9,10], is an inevitable important issue. Thus, extensive researches have been performed to study the creep mechanism of nanometallic materials, through both experiments and simulations [9,10,11,12,13,14,15,16,17,18].
Previous simulation results have clarified that the main creep mechanisms in NC materials are the lattice diffusion, grain boundary (GB) diffusion, GB sliding and dislocation activity. Different applied stresses may activate different creep mechanisms [9,19,20]. The well-known Brid-Dorn-Mukherjee classical equation [21], which comprehensively describes the effects of external stress, temperature and grain size, is widely used to analyze the creep behavior of NC materials. Its expression is as follows [21]:
ε ˙ = A D 0 G b k B T ( b d ) p ( σ G ) n exp ( Δ Q k B T )
where ε ˙ is the creep strain rate of the material in the steady-state creep stage, A is a dimensionless constant, D 0 is the diffusion coefficient, G is the shear modulus of the material, b is the magnitude of the Burgers vector, kB is the Boltzmann constant, T is the absolute temperature, d is the average grain size of the polycrystalline material, σ is the applied stress, Δ Q is the activation energy of a specific thermally activated creep mechanism, and p and n are the grain size exponent and stress exponent, respectively. Among them, the exponents n and p play an important role in determining the creep mechanism of materials. Specifically, different exponents correspond to different creep mechanisms. When n = 1, the creep deformation is always dominated by diffusion, and in this case, if p = 2, the vacancy diffusion occurs in the crystal, which is usually called the lattice diffusion (i.e., the Nabarro-Herring diffusion) [22]; if p is equal to 3, the vacancies flow along the GBs, which is called the GB diffusion, that is, the Coble diffusion [13]. When the stress exponent n = 2 or the grain size exponent p = 3, the creep process is dominated by the GB slip [12,19]. At the same time, it is worth noting that, the diffusion and slip of GBs always occurs simultaneously [19,23]. When the stress exponent n is larger than 4, the nucleation, slip, and climb of dislocations are the dominant creep deformation mechanisms [23], which is called as the dislocation creep.
One of the most intriguing properties of NC materials that are utilized at various applications is the ultra-high strength [24,25]. It is well accepted that its ultra-high strength is due to the existence of a large number of GBs and the related structures in the crystal. Generally speaking, the strength of NC materials increases with the decrease of grain size. However, when the grain size is refined to a certain extent, for example, about 10 nm to 30 nm, NC materials will soften, i.e., the strength will decrease with the further decrease of grain size [26]. This is mainly caused by the transition of plastic deformation mechanism from dislocation activities to GB behaviors [27]. At the same time, due to the large number of GB structures in NC materials, these materials are sensitive to high temperatures. The instability of GBs also greatly limits the application of NC materials in high temperature areas [28].
Accordingly, strategies against the softening behavior of NC materials caused by GB instability at extremely fine grain size are essential for the potential applications, which have attracted continuous efforts. The experimental results of Lu et al. showed that adding the second-phase atoms and inducing them to segregate towards GBs can effectively avoid the softening phenomenon [29]. Molecular dynamics simulations performed by Sellers et al. also showed that the segregation of Cu and Ag atoms at the GBs of β-Sn material affects the GB energy and shear stress [30]. The study conducted by Meyers et al. suggested that the segregated impurity atoms can pin the GBs, thus reducing the mobility of the GBs [31]. Based on these studies, it can be known that the segregation of the second-phase atoms to GBs has a significant effect on the GB behaviors. That is, the GB segregation effect can change the mechanical properties of NC materials by affecting the deformation process dominated by the GB behaviors.
The research based on molecular dynamics simulations that was conducted by Keblinski et al. indicated the Coble creep (the creep process dominated by GB diffusion) is one of the main creep mechanisms of ultrafine grained materials [32,33]. Therefore, the GB diffusivity, which was proved to be affected by the orientation of GBs, interface atoms diffusion, and the segregation of solute atoms at GBs, is a very important factor affecting the creep of NC materials [34,35,36]. Thus, is it possible to regulate the creep resistance of NC materials through the segregation of solute atoms at GBs? Some hints can be obtained from existing studies. For example, the simulation studies performed by Schafer et al. showed that the segregation of some certain types of solute atoms at GBs can improve the creep resistance of Cu materials significantly, and different sizes of the second-phase atoms have different influences [37]. The experiments based on Mg alloys also indicated that the segregation of solute atoms can improve the creep resistance of materials [38]. However, at present, there is still a lack of systematic research on the effect of segregation on the creep properties and mechanisms of NC materials. More importantly, understanding the creep mechanisms of ultrafine grained metal materials and their alloys at atomic scale is essential for the design of materials with stable mechanical and thermodynamic properties under various working conditions.
To study the creep properties of NC Ni with different grain sizes and concentrations of second-phase atoms at GBs under different external conditions systematically, a series of large-scale molecular dynamics simulations are carried out in this paper. The corresponding creep mechanisms are also revealed.

2. Simulation Methods and Parameters

All the three-dimensional polycrystalline models used in this paper are constructed by the ATOMSK software (ATOMSK-Beta 0.10.6, Pierre Hirel, Villeneuve d’Ascq, France) based on the Voronoi method [39]. Three types of models are constructed: the first one is the NC Ni sample without impurity atoms, the second type is the NC Ni sample with a certain proportion of Mo atoms randomly distributed in grains, and the third type is the NC Ni sample with a certain proportion of Mo atoms segregated at GBs, examining the effect of segregation on the creep behaviors of NC materials. The detailed simulation parameters are shown in Table 1. The dimensions of all the models are 35 nm × 35 nm × 35nm, containing about 4 million atoms. The average grain sizes range from 7.9 nm to 19.8 nm. Figure 1 gives the atomic configuration of the 3D polycrystalline Ni sample with the grain size of 7.9 nm, in which green and white atoms represent Ni atoms located in the grain interior and GBs, respectively, and black atoms represent the Mo atoms segregated at GBs. The influences of the addition and segregation of Mo atoms on the creep behaviors of NC Ni are studied through the cross-comparison of three types of models. The effects of applied stress, temperature, and grain size on the creep deformation are further investigated using the NC Ni sample with Mo atoms segregated at GBs.
All simulations in this work are performed using the LAMMPS package (LAMMPS-22 August 2018, Sandia National Laboratories, Albuquerque, NM, USA) [40]. The velocity-Verlet integration algorithm is used with the time step being set to 2 fs. The interaction between Ni and Mo atoms is described by the embedded atom potential parameterized by Zhou et al. [41], which can accurately predict some basic properties of materials, such as the lattice parameter, elastic constants, bulk modulus, vacancy formation energy, and so on. Besides, the GB energy of Ni calculated using this potential is in the range of the mean value of GB energy for the NC metals. Periodic boundary conditions are used in all three directions during the simulation. In order to accelerate the process of creep deformation with practical simulation time, a high temperature of 1100 K is adopted in all the simulations except for those used to study the effect of temperature.
During the simulation, each sample is first relaxed by the conjugated gradient method to get the equilibrium configuration at 0 K. Then, by using the Nosé-Hoover method [42], the sample is relaxed under an isothermal-isobaric ensemble for 200 ps to obtain the equilibrium configuration at the given temperature. Starting from the equilibrium, a constant stress is applied in the y direction for 200 ps to simulate creep load. The pressures in x and z directions are kept zero during the simulation. The microstructures are identified by the common neighbor analysis method [43], the face-centered-cubic, hexagonal-close-packed and non-structured atoms are colored green, red, and white, respectively.

3. Results and Discussion

3.1. Effect of Segregation on Creep Behavior

To investigate the effect of the segregation of Mo atoms on creep behaviors of NC Ni, creep simulations of the pure Ni samples and the Ni samples with Mo atoms segregated at GBs (referred as segregated Ni-Mo sample thereinafter) are carried out. Figure 2a,b exhibits the evolution of creep strain with respect to time in the 10.9 nm grained Ni and segregated Ni-Mo (3 at.%) samples under different constant stresses at 1100 K. According to the variation of the strain rate on each curve, a typical creep curve under a low stress can be divided into two stages: the initial creep stage and the steady creep stage. During the initial creep stage, the creep strain increases rapidly, the creep strain rate is very large, but gradually decreases with the increase of time. The creep strain rate in the initial stage reflects the instantaneous change of creep strain with time. During the steady creep stage, the strain increases almost linearly with time, implying that the creep rate is almost constant. Under high stresses, the third stage, namely the accelerated creep stage, could be observed, as shown in the curves corresponding to the stress over 1.4 GPa in Figure 2a and over 1.6 GPa in Figure 2b. In this stage, the creep strain rate increases gradually with time until the material breaks. It can also be found that in a certain time range, the critical stress for the third creep stage of the segregated sample is higher than that of the pure Ni sample, which indicates that the segregation of the second-phase atoms at GBs can improve the ability of NC materials to resist creep deformation.
By fitting the slope of creep curves in the steady state stage, the creep rates under different stresses are obtained. Then, the creep rates versus different applied stresses of pure Ni samples and segregated Ni-Mo samples are plotted in Figure 2c,d, respectively. Figure 2c,d shows that the creep rates of both the pure Ni samples and segregated Ni-Mo samples increase significantly with the increase of applied stress. According to the Bird-Dorn-Mukherjee equation, i.e., Equation (1), the stress exponent n can be obtained by fitting the slope of the dual logarithmic graph of creep rate and stress, which essentially reflects the strain rate sensitivity of the materials. Here, it can be seen from Figure 2c,d that in the low stress range, the values of n for both the pure Ni sample and segregated Ni-Mo sample are greater than 1 and less than 2, indicating that atoms diffusion dominates the creep deformation process within this stress range, but there is also a GB slip phenomenon during the deformation. In practical terms, GB diffusion and GB slip are generally coupled because the diffusion of the atoms at GBs will lead to the change of grain morphology, so the deformation compatibility at the GBs should be maintained by sliding between grains [19]. In the medium stress state, the values of n for pure Ni and segregated Ni-Mo samples are about 2, which implies that the GB slip is the main deformation mechanism within this stress range. In the high stress state, for pure Ni samples, the stress exponent approaches 4, suggesting that both the GB sliding and dislocation activities exist in the creep deformation process, but the dominated mechanism is the dislocation activities; for segregated Ni-Mo samples, the stress exponent is greater than 4, implying that the creep is governed by dislocation activities.
By comparing Figure 2c,d, it can be seen that although the stress exponent values of the segregated Ni-Mo sample are larger than those of pure Ni sample, the former is lower than the latter for a given stress. To further see the difference of creep rate between the two NC Ni samples more clearly, we put the creep rate of the two samples under different stress into the same figure, as shown in Figure 3. It can be found that the creep rate of polycrystalline Ni decreases significantly when there are Mo atoms segregated at the GBs, and the strengthening effect becomes more obvious with the increase of stress within a certain stress range.
To eliminate the influence of the addition of Mo atoms on the creep resistance of NC Ni, and further prove that the increase of creep resistance of NC Ni is due to the segregation effect at the GBs, we conduct a creep simulation of a group of the NC Ni samples with Mo atoms randomly distributed in the grains. The results show that the creep rate of the NC Ni sample with Mo atoms distributed in the grains is slightly higher than that of the pure Ni sample. Through linearly fitting the relationship between the creep rate and reciprocal temperature of these three sets of NC Ni models, the activation energies during the creep are obtained, and the results are shown in Figure 4. The results reveal that the activation energy of the segregated Ni-Mo sample is about 0.65 eV, which is significantly higher than that of the NC Ni sample without Mo atoms. However, the activation energy of the Ni sample with the Mo atoms distributed in the grains is about 0.34 eV, which is slightly lower than that of pure Ni sample. It should be noted that all the activation energies derived from molecular dynamics simulations are in the same order of magnitude as those observed in experiments [9]. The results show that the addition of Mo atoms randomly distributed in the grains reduces the activation energy of NC Ni, making creep easier to occur. This may be due to the lattice distortion caused by Mo atoms distributed in the grains, which increases the free energy of the material, and thereby reduces the external activation energy required for creep. The results also reveal that the segregation of Mo atoms at the GBs can visibly improve the activation energy of the material, which makes the creep difficult to occur. This also proves that the remarkable improvement of the creep resistance of NC Ni is due to the segregation effect of Mo atoms at the GBs.
Based on the analysis above, it can be seen that the NC Ni materials have better creep resistance when there are Mo atoms segregated at the GBs, which is consistent with some previous simulation and experimental results [37,38]. In the next section, the potential creep mechanism of NC materials under different applied stresses and the strengthening mechanism of segregation are discussed in detail.

3.2. Creep Mechanism Analysis

To better understand the inherent mechanisms of the effect of segregation at GBs on creep behaviors of NC materials, the microstructure evolution of NC pure Ni and segregated Ni-Mo samples under different stresses are investigated. Figure 5a exhibits the initial atomic configuration of the 10.9 nm grained pure Ni sample, and the deformed configurations of this sample after creep for 200 ps at the stress of 0.6 GPa, 1.0 GPa, and 1.6 GPa are shown in Figure 5b–d, respectively. In Figure 5b, the grain morphology changes obviously after creep at the stress of 0.6 GPa, compared with the initial atomic configuration before creep in Figure 5a. The grains with obvious morphological changes are marked by numbers in Figure 5a,b: the sizes of grains 2 and 3 increase obviously, while the size of grain 1 between them decreases; grain 4 gobbles up grain 5 gradually; grains 6 and 7 grow while grains 8 and 9 decrease in size and merge into one grain. These phenomena indicate that GB diffusion and GB slip occurred during the creep deformation. Besides, there is no dislocation found in the atomic configuration, implying that the dominant creep mechanisms under low stress states are the GB diffusion and GB slip, and this is consistent with the prediction of Figure 2c based on Equation (1). When the applied stress is 1.0 GPa, in addition to the obvious change of grain morphology, a small amount of dislocations are also observed in the atomic configuration, as indicated by the white arrows in Figure 5c, which illustrates that the prevalent mechanisms of the creep deformation are still GB behaviors, accompanied with a few dislocation activities. Figure 5d is the configuration after creep under the stress of 1.6 GPa, it is obvious to find that extensive partial dislocations nucleate at the GBs, and then pass through the grain interior, leaving behind stacking faults, as indicated by the white arrows. This indicates that the dislocation activities play a more important role during the creep deformation under high stress, which is also in line with the prediction of Figure 2c based on Equation (1).
Figure 6a–d shows the atomic configurations of segregated Ni-Mo (3 at.%) samples before creep and after creep for 200 ps under various applied stresses, respectively. Comparing Figure 6b,c with the initial atomic configuration in Figure 6a, it is found that the grain morphologies do not change obviously and there are almost no dislocation activities. Combined with the creep mechanisms predicted above according to Figure 2d, this implies that in the low and medium stress area, the creep deformation of segregated Ni-Mo sample is very small. Therefore, it can be found that the segregation of Mo atoms at GBs can effectively inhibit the creep process mediated by GB behaviors, thus improving the creep resistance of NC materials. When the stress is 1.6 GPa, a large number of stacking faults are observed in the grain interior, as indicated by the white arrows in Figure 6d. This means that the creep deformation is governed by dislocation activities when the applied stress is in the high stress regime, which is consistent with previous analysis [9,10].

3.3. Influences of Grain Size and Temperature

Generally speaking, as surface defects, the GBs in NC materials are in thermodynamically non-equilibrium state, and the creep is essentially a thermal activation process, so the GB and temperature both have important influence on the creep properties of materials. Therefore, this section focuses on the effects of grain size and temperature on the creep behaviors of segregated Ni-Mo samples. According to Equation (1), we can know that the creep strain rate will decrease with the increase of grain size for a fixed temperature and applied stress. The effect of the grain size on creep behaviors of the material is mainly characterized by the grain size exponent p. Figure 7 gives the relationship between the creep strain rate and the reciprocal of grain size of segregated Ni-Mo (3 at.%) sample under different applied stresses at 1100k. The double logarithmic plot in Figure 7 implies that the creep rate always increases with the decrease of grain size for all stress levels, and in a power-law relationship. From the figure, we can also see that the grain size exponents are about 1.92–2.36 for the applied stress within the low and medium stress range in Figure 2d. When the stress is in the range of 0.6–0.8 GPa, the values of stress exponent p are close to 2, which just corresponds to the Nabarro-Herring creep. When the stress continues to increase to 1.0 GPa, the stress exponent p is greater than 2 and less than 3, but closer to 2, which indicates that both the Nabarro-Herring and Coble diffusion or GB slip exist in the creep deformation process, but the dominant creep deformation mechanism is still the Nabarro-Herring diffusion. It can be seen that the mechanisms of lattice diffusion, GB diffusion and GB slip coexist under the low and medium stress state, but the segregation of Mo atoms at the GBs dramatically impedes the creep process mediated by GB behaviors. In addition, it should be noted that when the stress is 2.0 Gpa, within the high stress range in Figure 2d, the grain size exponent is 1.32, which means that the dislocation activities dominate the creep process.
Figure 8 exhibits the creep curves of segregated Ni-Mo (3 at.%) sample with grain size of 10.9 nm under an external stress of 1.0 GPa at different temperatures. It can be seen that for the segregated NC material, the creep rate increases with the increase of temperature, which is similar to that of pure NC materials [19,20,44]. This is because all creep mechanisms are essentially thermal activation processes.

3.4. Effect of Segregated Atoms Concentration at GBs

According to our previous quantitative research on the influence of segregation on GB energy of Ni [26], it can be known that, within a certain range, the GB energy decreases linearly with the increase of the concentration of Mo atoms at GBs, and the concentration has a significant effect on the plastic deformation process dominated by GB behaviors of NC Ni. The GB behaviors are also the important creep mechanisms of NC materials, so the concentration of Mo atoms segregated at the GBs also has an important influence on the creep behavior of NC Ni. To understand this influence more clearly, this section mainly investigates the creep behavior of polycrystalline Ni with different concentrations of Mo atoms segregated at GBs. Figure 9 shows the variation of creep strain rate with the concentration of Mo atoms segregated at GBs for the sample with an average grain size of 10.9 nm under different stresses. It can be seen that whether in the low stress regime or the medium and high stress regimes, as the Mo atoms concentration increases from 0 to 6 at.%, the creep strain rate of the material gradually decreases. This demonstrates that the creep resistance of NC materials can be optimized by adjusting the concentration of the second-phase atoms segregated at GBs. This may be useful for the design and fabrication of NC materials with excellent creep resistance.
In addition, it can also be found from Figure 9 that when the stress is between 1.0 GPa and 2.0 GPa, the influence of the concentration of Mo atoms on the creep strain rate is stronger than that when the applied stress is 0.6 GPa, which is due to the different creep mechanisms of materials under different stress states. From our previous work, we can know that the segregation of solute atoms at GBs can significantly improve the stability and reduce the mobility of GBs [26], which implies that the segregation of second-phase atoms at GBs has a greater influence on the creep process dominated by GB behaviors, such as when the stress is 1.0 GPa. From the analysis of grain size exponent above, it can be known that when the stress is 0.6 GPa, the deformation corresponds to the Nabarro-Herring creep, that is, the creep process is dominated by the lattice diffusion. Therefore, the segregation of Mo atoms at GBs has little effect on the creep rate, and with the increase of the concentration of Mo atoms, the creep rate quickly reaches a steady state and almost no longer decreases. Moreover, it is worth noting that the results in Figure 9 show that in the high stress region, i.e., 1.6 GPa and 2.0 GPa, the segregation of Mo atoms at GBs also has a significant effect on the creep process, while the analyses of stress and grain size exponents both indicate that the dominant mechanism of creep deformation is dislocation activities within this stress regime. This means that the segregation also has striking effect on creep deformation dominated by dislocation activities.
To explore the intrinsic mechanism of the influence of segregation on creep deformation dominated by dislocation activities, the dislocation densities of the 10.9 nm grained segregated Ni-Mo samples with different Mo concentrations after creep for 200 ps at 1100 K are shown in Figure 10. The corresponding data points at 0 GPa represent the dislocation density of equilibrium configuration before loading. It can be clearly seen from Figure 10 that when the stress is 0.6 GPa or 1.0 GPa, for all segregated Ni-Mo samples, the dislocation densities after creep 200 ps are basically the same as that in the corresponding equilibrium configurations. This also proves that there is almost no dislocation nucleation during the creep under low and medium stress states, and the creep deformation is dominated by the diffusion and GB sliding, in accordance with the previous analysis. When the applied stress increases to 1.6 GPa or 2.0 GPa, it can be found that the dislocation densities of all samples increase significantly compared with the equilibrium configurations. This trend also proves that the main creep mechanism under a high stress state is the dislocation activity. It should be noted that the dislocation densities of segregated samples with different Mo concentrations are much lower than that of the pure Ni sample, and the higher the proportion of Mo atoms, the lower the dislocation density in the segregated sample, which indicates that the segregation of Mo atoms at GBs inhibits the nucleation and proliferation of dislocation in NC Ni during the creep. When the proportion of Mo atoms is higher than 4 at.%, the dislocation density in the material under the applied stress of 1.6 GPa is even reduced to the same level of equilibrium configuration.
To show more evidences about the proposed dislocation activity as the dominant mechanism under high stress, we picked one grain among the 10.9 nm grained segregated Ni-Mo (6 at.%) sample during the steady state creep at 1.6 GPa, as shown in Figure 11. It can be obviously found that due to the limitation of grain size, the dominant creep mechanism is the dislocation nucleation from GBs rather than the collective dislocation dynamics inside grains, which has been confirmed by the previous studies on the deformation mechanism of NC materials [45]. By comparing the dislocation density of all samples at equilibrium state in Figure 10, it can be found that the dislocation density decreases with the increase of the concentration of Mo atoms, which implies that the segregation of Mo atoms at GBs destroys the dislocation structure at GBs. However, the dislocation structures at the GBs are the sources of dislocation nucleation during the creep, so the segregation of Mo atoms at GBs can also significantly reduce the creep rate of materials under high stress.

4. Conclusions

A train of large-scale molecular dynamics simulations are carried out to investigate the influence of the segregation of Mo atoms on the high-temperature creep behavior of NC Ni. The effects of the applied stress, grain size, temperature, and concentration of Mo atoms on the creep behaviors and mechanisms are systematically studied. The main conclusions are summarized as follows:
(1)
The segregated Ni-Mo sample corresponds to a higher activation energy than the pure Ni sample, which makes the creep of the NC Ni more difficult to occur and thus enhances the creep resistance of the material;
(2)
For both the pure Ni sample and segregated Ni-Mo sample, the creep mechanisms are the diffusion, GB slip, and dislocation activity in the low, medium, and high stress regimes, respectively;
(3)
The segregation of Mo atoms at GBs has little effect on the creep processes dominated by the lattice diffusion, but significantly slows down the creep processes dominated by the GB behavior and dislocation activity.

Author Contributions

Conceptualization, Y.Z.; Data curation, Q.L. and J.Z.; Formal analysis, Q.L. and Y.Z.; Funding acquisition, Y.Z. and H.Z.; Investigation, Q.L., J.Z., H.T., H.Z., H.Y. and Y.Z.; Methodology, Q.L., J.Z. and Y.Z.; Supervision, Y.Z.; Writing—original draft, Q.L.; Writing—review & editing, Y.Z., H.Y. and H.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the National Natural Science Foundation of China, grant numbers 12072062, 12072061 and 11772082, Liaoning Revitalization Talents Program, grant number XLYC1807193, Key Research and Development Project of Liaoning Province, grant number 2020JH2/10500003, and Fundamental Research Funds for the Central Universities, grant number DUT20LAB203.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The raw/processed data required to reproduce these findings cannot be shared at this time as the data also forms part of an ongoing study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Schiøtz, J.; Jacobsen, K.W. A Maximum in the strength of nanocrystalline copper. Science 2003, 301, 1357–1359. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Chen, J.; Lu, L.; Lu, K. Hardness and strain rate sensitivity of nanocrystalline Cu. Scr. Mater. 2006, 54, 1913–1918. [Google Scholar] [CrossRef]
  3. Zheng, Y.G.; Lu, J.; Zhang, H.W.; Chen, Z. Strengthening and toughening by interface-mediated slip transfer reaction in nanotwinned copper. Scr. Mater. 2009, 60, 508–511. [Google Scholar] [CrossRef]
  4. Zheng, Y.G.; Zhao, Y.T.; Ye, H.F.; Zhang, H.W. Size-dependent elastic moduli and vibrational properties of fivefold twinned copper nanowires. Nanotechnology 2014, 25, 315701. [Google Scholar] [CrossRef] [PubMed]
  5. Lu, L.; Chen, X.; Huang, X.; Lu, K. Revealing the maximum strength in nanotwinned copper. Science 2009, 323, 607–610. [Google Scholar] [CrossRef]
  6. Xiong, Z.; Zhong, L.; Wang, H.; Li, X. Structural defects, mechanical behaviors and properties of two-dimensional materials. Materials 2021, 14, 1192. [Google Scholar] [CrossRef]
  7. Li, X.; Lu, L.; Li, J.; Zhang, X.; Gao, H. Mechanical properties and deformation mechanisms of gradient nanostructured metals and alloys. Nat. Rev. Mater. 2020, 5, 706–723. [Google Scholar] [CrossRef]
  8. Yamamoto, Y.; Brady, M.P.; Lu, Z.P.; Maziasz, P.J.; Liu, C.T.; Pint, B.A.; More, K.L.; Meyer, H.M.; Payzant, E.A. Creep-resistant, Al2O3-forming austenitic stainless steels. Science 2007, 316, 433–436. [Google Scholar] [CrossRef]
  9. Yang, X.S.; Wang, Y.J.; Zhai, H.R.; Wang, G.Y.; Su, Y.J.; Dai, L.H.; Ogata, S.; Zhang, T.Y. Time-, stress-, and temperature-dependent deformation in nanostructured copper: Creep tests and simulations. Acta Mater. 2016, 94, 191–206. [Google Scholar] [CrossRef] [Green Version]
  10. Nie, K.; Wu, W.P.; Zhang, X.L.; Yang, S.M. Molecular dynamics study on the grain size, temperature, and stress dependence of creep behavior in nanocrystalline nickel. J. Mater. Sci. 2017, 52, 2180–2191. [Google Scholar] [CrossRef]
  11. Meraj, M.; Pal, S. The effect of temperature on creep behaviour of porous (1 at.%) nanocrystalline nickel. Trans. Indian Inst. Met. 2015, 69, 277–282. [Google Scholar] [CrossRef]
  12. Stevens, R.N. Grain-boundary sliding and diffusion creep in polycrystalline solids. Philos. Mag. 1971, 23, 265–283. [Google Scholar] [CrossRef]
  13. Coble, R.L. A model for boundary diffusion controlled creep in polycrystalline materials. J. Appl. Phys. 1963, 34, 1679–1682. [Google Scholar] [CrossRef]
  14. Nabarro, F.R.N. Steady-state diffusional creep. Philos. Mag. 1967, 16, 231–237. [Google Scholar] [CrossRef]
  15. Bhatia, M.A.; Mathaudhu, S.N.; Solanki, K.N. Atomic-scale investigation of creep behavior in nanocrystalline Mg and Mg-Y alloys. Acta Mater. 2015, 99, 382–391. [Google Scholar] [CrossRef] [Green Version]
  16. Millett, P.C.; Desai, T.; Yamakov, V.; Wolf, D. Atomistic simulations of diffusional creep in a nanocrystalline body-centered cubic material. Acta Mater. 2008, 56, 3688–3698. [Google Scholar] [CrossRef]
  17. Ford, J.M.; Wheeler, J.; Movchan, A. Computer simulation of grain-boundary diffusion creep. Acta Mater. 2002, 50, 3941–3955. [Google Scholar] [CrossRef]
  18. Yamakov, V.; Wolf, D.; Phillpot, S.R.; Gleiter, H. Grain-boundary diffusion creep in nanocrystalline palladium by molecular-dynamics simulation. Acta Mater. 2002, 50, 61–73. [Google Scholar] [CrossRef]
  19. Wang, Y.J.; Ishii, A.; Ogata, S. Transition of creep mechanism in nanocrystalline metals. Phys. Rev. B 2011, 84, 224102. [Google Scholar] [CrossRef]
  20. Zeng, Y.; Li, X. Atomistic simulations of high-temperature creep in nanotwinned TiAl alloys. Extrem. Mech. Lett. 2021, 44, 101253. [Google Scholar] [CrossRef]
  21. Mukherjee, A.K.; Bird, J.E.; Dorn, J.E. Experimental correlations for high-temperature creep. In Proceedings of the Detroit Materials Engineering Congress of American Society for Metals, Detroit, MI, USA, 15 October 1968. [Google Scholar]
  22. Herring, C. Diffusional viscosity of a polycrystalline solid. J. Appl. Phys. 1950, 21, 437–445. [Google Scholar] [CrossRef]
  23. Jiao, S.; Kulkarni, Y. Molecular dynamics study of creep mechanisms in nanotwinned metals. Comput. Mater. Sci. 2015, 110, 254–260. [Google Scholar] [CrossRef]
  24. Kelsall, R.W.; Hamley, I.W.; Geoghegan, M. Nanoscale Science and Technology; Wiley: New York, NY, USA, 2005. [Google Scholar]
  25. Murty, B.S.; Shankar, P.; Raj, B.; Al, E. Textbook of Nanoscience and Nanotechnology; Springer: Berlin, Germany, 2013. [Google Scholar]
  26. Li, Q.; Zhang, J.; Tang, H.; Ye, H.; Zheng, Y. Regulating the mechanical properties of nanocrystalline nickel via molybdenum segregation: An atomistic study. Nanotechnology 2019, 30, 275702. [Google Scholar] [CrossRef]
  27. Meyers, M.A.; Hahn, E.N. Grain-size dependent mechanical behavior of nanocrystalline metals. Mater. Sci. Eng. A 2015, 646, 101–134. [Google Scholar] [CrossRef] [Green Version]
  28. Padmanabhan, K.A.; Sripathi, S.; Hahn, H.; Gleiter, H. Inverse Hall–Petch effect in quasi- and nanocrystalline materials. Mater. Lett. 2014, 133, 151–154. [Google Scholar] [CrossRef]
  29. Hu, J.; Shi, Y.N.; Sauvage, X.; Sha, G.; Lu, K. Grain boundary stability governs hardening and softening in extremely fine nanograined metals. Science. 2017, 355, 1292–1296. [Google Scholar] [CrossRef]
  30. Sellers, M.S.; Schultz, A.J.; Kofke, D.A.; Basaran, C. Solute effects on β-Sn grain boundary energy and shear stress. J. Nanomechanics Micromechanics 2011, 1, 41–50. [Google Scholar] [CrossRef] [Green Version]
  31. Meyers, M.A.; Mishra, A.; Benson, D.J. Mechanical properties of nanocrystalline materials. Prog. Mater. Sci. 2006, 51, 427–556. [Google Scholar] [CrossRef]
  32. Keblinski, P.; Wolf, D.; Gleiter, H. Molecular-dynamics simulation of grain-boundary diffusion creep. Interface Sci. 1998, 6, 205–212. [Google Scholar] [CrossRef]
  33. Meraj, M.; Yedla, N.; Pal, S. The effect of porosity and void on creep behavior of ultra-fine grained nano crystalline nickel. Mater. Lett. 2016, 169, 265–268. [Google Scholar] [CrossRef]
  34. Lee, Y.; Basaran, C. Effect of Ni solute on grain boundary diffusivity and structure of β-Sn. Comput. Mater. Sci. 2014, 92, 1–7. [Google Scholar] [CrossRef]
  35. Michael, S.; Sellers, M.S.; Andrew, J.; Schultz, A.J.; Basaran, C.; David, A. β-Sn grain-boundary structure and self-diffusivity via molecular dynamics simulations. Phys. Rev. B 2010, 81, 134111. [Google Scholar] [CrossRef]
  36. Sellers, M.S.; Schultz, A.J.; Basaran, C.; Kofke, D.A. Atomistic modeling of β-Sn surface energies and adatom diffusivity. Appl. Surf. Sci. 2010, 256, 4402–4407. [Google Scholar] [CrossRef]
  37. Schäfer, J.; Ashkenazy, Y.; Albe, K.; Averback, R.S. Effect of solute segregation on thermal creep in dilute nanocyrstalline Cu alloys. Mater. Sci. Eng. A 2012, 546, 307–313. [Google Scholar] [CrossRef]
  38. He, S.M.; Zeng, X.Q.; Peng, L.M.; Gao, X.; Nie, J.F.; Ding, W.J. Microstructure and strengthening mechanism of high strength Mg–10Gd–2Y–0.5Zr alloy. J. Alloys Compd. 2007, 427, 316–323. [Google Scholar] [CrossRef]
  39. Hirel, P. Atomsk: A tool for manipulating and converting atomic data files. Comput. Phys. Commun. 2015, 197, 212–219. [Google Scholar] [CrossRef]
  40. Plimpton, S. Fast parallel algorithms for short-range molecular dynamics. J. Comput. Phys. 1995, 117, 1–19. [Google Scholar] [CrossRef] [Green Version]
  41. Zhou, X.W.; Johnson, R.A.; Wadley, H.N.G. Misfit-energy-increasing dislocations in vapor-deposited CoFe/NiFe multilayers. Phys. Rev. B 2004, 69, 144113. [Google Scholar] [CrossRef] [Green Version]
  42. Nosé, S. A unified formulation of the constant temperature molecular dynamics methods. J. Chem. Phys. 1984, 81, 511–519. [Google Scholar] [CrossRef] [Green Version]
  43. Stukowski, A. Visualization and analysis of atomistic simulation data with OVITO-the open visualization tool. Model. Simul. Mater. Sci. Eng. 2010, 18, 015012. [Google Scholar] [CrossRef]
  44. Karanjgaokar, N.; Chasiotis, I. Creep behavior of nanocrystalline Au films as a function of temperature. J. Mater. Sci. 2016, 51, 3701–3714. [Google Scholar] [CrossRef]
  45. Li, X.; Wei, Y.; Lu, L.; Lu, K.; Gao, H. Dislocation nucleation governed softening and maximum strength in nano-twinned metals. Nature 2010, 464, 877–880. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Atomic configuration of segregated polycrystalline Ni-Mo model with d = 7.9 nm.
Figure 1. Atomic configuration of segregated polycrystalline Ni-Mo model with d = 7.9 nm.
Materials 14 06966 g001
Figure 2. (a,b) Creep curves of the NC Ni and the segregated Ni-Mo (3 at.%) samples with d = 10.9 nm at 1100 K under various applied stresses, respectively; (c,d) Dual logarithmic plots of the creep strain rate versus stress of NC Ni sample and segregated Ni-Mo (3 at.%) samples with d = 10.9 nm at 1100 K, respectively.
Figure 2. (a,b) Creep curves of the NC Ni and the segregated Ni-Mo (3 at.%) samples with d = 10.9 nm at 1100 K under various applied stresses, respectively; (c,d) Dual logarithmic plots of the creep strain rate versus stress of NC Ni sample and segregated Ni-Mo (3 at.%) samples with d = 10.9 nm at 1100 K, respectively.
Materials 14 06966 g002
Figure 3. The comparison of the creep strain rate of two groups of NC Ni samples with d = 10.9 nm at 1100 K under different applied stresses.
Figure 3. The comparison of the creep strain rate of two groups of NC Ni samples with d = 10.9 nm at 1100 K under different applied stresses.
Materials 14 06966 g003
Figure 4. Arrhenius plot of the creep rate against the reciprocal temperature of NC Ni with d = 10.9 nm under an applied stress of 2 GPa.
Figure 4. Arrhenius plot of the creep rate against the reciprocal temperature of NC Ni with d = 10.9 nm under an applied stress of 2 GPa.
Materials 14 06966 g004
Figure 5. Snapshots of a pure Ni sample with d = 10.9 nm at 1100 K (a) before creep; (bd) after creep for 200 ps under an applied stress of 0.6 GPa, 1.0 GPa and 1.6 Gpa, respectively.
Figure 5. Snapshots of a pure Ni sample with d = 10.9 nm at 1100 K (a) before creep; (bd) after creep for 200 ps under an applied stress of 0.6 GPa, 1.0 GPa and 1.6 Gpa, respectively.
Materials 14 06966 g005
Figure 6. Snapshots of segregated Ni-Mo (3 at.%) samples with d = 10.9 nm at 1100 K (a) before creep; (bd) after creep for 200 ps under an applied stress of 0.6 GPa, 1.0 GPa and 1.6 GPa, respectively.
Figure 6. Snapshots of segregated Ni-Mo (3 at.%) samples with d = 10.9 nm at 1100 K (a) before creep; (bd) after creep for 200 ps under an applied stress of 0.6 GPa, 1.0 GPa and 1.6 GPa, respectively.
Materials 14 06966 g006
Figure 7. The relationship between the creep strain rate and the reciprocal of grain size d of the segregated Ni-Mo (3 at.%) sample at 1100 K under different applied stresses.
Figure 7. The relationship between the creep strain rate and the reciprocal of grain size d of the segregated Ni-Mo (3 at.%) sample at 1100 K under different applied stresses.
Materials 14 06966 g007
Figure 8. Creep curves of the segregated Ni-Mo (3 at.%) sample with d = 10.9 nm at different temperatures under an applied stress of 1.0 GPa.
Figure 8. Creep curves of the segregated Ni-Mo (3 at.%) sample with d = 10.9 nm at different temperatures under an applied stress of 1.0 GPa.
Materials 14 06966 g008
Figure 9. Variation in creep strain rate with the Mo concentration of segregated Ni-Mo samples with d = 10.9 nm under different stresses: (a) 0.6 GPa and 1.0 GPa, and (b) 1.6 GPa and 2.0 GPa.
Figure 9. Variation in creep strain rate with the Mo concentration of segregated Ni-Mo samples with d = 10.9 nm under different stresses: (a) 0.6 GPa and 1.0 GPa, and (b) 1.6 GPa and 2.0 GPa.
Materials 14 06966 g009
Figure 10. Relationship between the dislocation density and the applied stress in segregated Ni-Mo samples with different Mo concentrations after creep at 1100 K for 200 ps, the situation of stress at 0 GPa represents the dislocation density of the equilibrium configuration at 1100 K.
Figure 10. Relationship between the dislocation density and the applied stress in segregated Ni-Mo samples with different Mo concentrations after creep at 1100 K for 200 ps, the situation of stress at 0 GPa represents the dislocation density of the equilibrium configuration at 1100 K.
Materials 14 06966 g010
Figure 11. Dislocation nucleation from GBs and the left stacking faults during creep deformation under an applied stress of 1.6 GPa at 1100 K. The snapshot is from the segregated sample with 6 at.% Mo atoms and d = 10.9 nm, in which the perfect fcc atoms are not shown for clarity.
Figure 11. Dislocation nucleation from GBs and the left stacking faults during creep deformation under an applied stress of 1.6 GPa at 1100 K. The snapshot is from the segregated sample with 6 at.% Mo atoms and d = 10.9 nm, in which the perfect fcc atoms are not shown for clarity.
Materials 14 06966 g011
Table 1. Parameters in the creep simulation of NC Ni-Mo alloy.
Table 1. Parameters in the creep simulation of NC Ni-Mo alloy.
Group σ (GPa) d (nm) T (K) Mo (at.%)
10.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.1, 1.2, 1.3, 1.4, 1.6, 1.8, 2.010.911000.0, 3.0
20.6, 0.8, 1.0, 2.07.9, 10.9,
13.7, 19.8
11003.0
31.010.9700, 800, 900, 1000, 11003.0
40.6, 1.0, 1.6, 2.010.911000.0, 1.0, 2.0, 3.0, 4.0, 5.0, 6.0
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, Q.; Zhang, J.; Tang, H.; Zhang, H.; Ye, H.; Zheng, Y. Influence of Mo Segregation at Grain Boundaries on the High Temperature Creep Behavior of Ni-Mo Alloys: An Atomistic Study. Materials 2021, 14, 6966. https://doi.org/10.3390/ma14226966

AMA Style

Li Q, Zhang J, Tang H, Zhang H, Ye H, Zheng Y. Influence of Mo Segregation at Grain Boundaries on the High Temperature Creep Behavior of Ni-Mo Alloys: An Atomistic Study. Materials. 2021; 14(22):6966. https://doi.org/10.3390/ma14226966

Chicago/Turabian Style

Li, Qian, Jiayong Zhang, Huayuan Tang, Hongwu Zhang, Hongfei Ye, and Yonggang Zheng. 2021. "Influence of Mo Segregation at Grain Boundaries on the High Temperature Creep Behavior of Ni-Mo Alloys: An Atomistic Study" Materials 14, no. 22: 6966. https://doi.org/10.3390/ma14226966

APA Style

Li, Q., Zhang, J., Tang, H., Zhang, H., Ye, H., & Zheng, Y. (2021). Influence of Mo Segregation at Grain Boundaries on the High Temperature Creep Behavior of Ni-Mo Alloys: An Atomistic Study. Materials, 14(22), 6966. https://doi.org/10.3390/ma14226966

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