Research Article | Open Access
V. I. Razumovskiy, D. Scheiber, I. M. Razumovskii, V. N. Butrim, A. S. Trushnikova, S. B. Varlamova, A. G. Beresnev, "New Cr-Ni-Base Alloy for High-Temperature Applications Designed on the Basis of First Principles Calculations", Advances in Condensed Matter Physics, vol. 2018, Article ID 9383981, 8 pages, 2018. https://doi.org/10.1155/2018/9383981
New Cr-Ni-Base Alloy for High-Temperature Applications Designed on the Basis of First Principles Calculations
We use ab initio calculations to analyze the influence of 4d and 5d transition metal alloying elements on cohesive properties of the bulk and a representative grain boundary in Cr within the framework of the Rice–Thomson–Wang approach. The results obtained for Cr are combined with the analogous results for Ni to select Ta and Nb as promising alloying additions to dual-phase () Cr-Ni-base high-temperature alloys. Ta and Nb are added to the alloying system of an existing alloy I (Cr-Ni-W-V-Ti) in an attempt to design a chemical composition of a new alloy II (Cr-Ni-W-V-Ti) + (Ta-Nb). Investigation of the microstructure of the Ta-bearing Cr-Ni-alloy reveals a Ta enrichment of large -areas near GBs in -matrix that we consider as potency to increase the cohesive strength of GBs and the cohesive energy of the bulk in -phase. Mechanical testing of alloys I and II demonstrates that the alloy II has improved tensile strength and creep resistance at high temperatures.
Traditional materials for production of critical components of gas-turbine engines, where significant resistance to loading at high temperatures is required, are Ni-base superalloys [1, 2]. One of the most important performance characteristics of Ni-base superalloys is their creep resistance. In the polycrystalline critical parts, for example, turbine disks [3, 4], one of the weakest elements of microstructure is grain boundaries (GBs). GBs are normally characterized by enhanced diffusivity , which leads to increased creep at high temperatures. One can include into the alloying system elements that segregate to GBs and increase their work of separation in order to strengthen GBs and improve their creep resistance. The is a fundamental thermodynamic quantity that controls the mechanical strength of an interface and allows one to study the problem of GB embrittlement from first principles using the Rice–Thomson–Wang model [6, 7]. In accordance with the low alloying additions concept [8, 9], the “useful” alloying addictions have to increase simultaneously the cohesive strength of GBs and the cohesive energy of the bulk. By such an approach, Ta, Nb, Zr, and Hf have been predicted to be beneficial alloying addictions to Ni-base alloys and used to design a new polycrystalline Ni-base superalloy for powder metallurgy applications .
The main weakness of Ni-base superalloys for high-temperature applications is insufficient corrosion resistance. One of the most effective approaches to improve this is alloying with Cr. Cr-Ni-base high-temperature alloys represent their own type of superalloys , which is of great practical importance, for instance, for manufacturing combustion chambers of thermocatalytic engines intended for correction of the orbit and orientation of a spacecraft . Such alloys typically have the following chemical composition: Cr-(31–35)Ni-(1–3)W-(0.1–0.4)V-(0.05–0.3)Ti, wt.% (alloy I, ), and they consist of two major metallic phases: (solid solution based on Cr with bcc lattice) and (solid solution based on Ni with fcc lattice). A distinguishing feature of their microstructure is formation of large -phase areas near GBs in -matrix. The deformation behavior of the whole + microstructure in this type of alloys is controlled by -phase due to its high degree of ductility .
At first sight, it seems attractive to add the same elements that are beneficial for Ni-base alloys, which are Zr, Hf, Nb, and Ta, also into Cr-Ni-base alloys to obtain stronger GBs in the -phase. However, the influence of the aforementioned elements on the cohesive properties of -matrix remains unknown and data dealing with Zr, Hf, Nb, and Ta influence on mechanical properties of Cr-base alloys are controversial . That is why it is important to investigate an embrittling potency of Zr, Hf, Nb, and Ta in Cr-base solid solutions as well.
In this paper, we investigate first the effect of a number of transition metals including Zr, Hf, Nb, and Ta on the cohesive properties of the bulk -phase and a representative GB in Cr alloys using ab initio calculations. Second, we determine the embrittling potency of 4d and 5d transition metals in Cr alloys in the framework of the low alloying additions concept [8, 9]. Finally, we introduce the elements acceptable for both and phases into a typical alloy to fabricate a new Cr-Ni-base alloy and investigate its microstructure and mechanical properties.
2. Computational Details
2.1. Electronic Structure Calculations
Electronic structure calculations have been performed using the projector augmented wave method  as implemented in Vienna ab initio simulation package (VASP) [19, 20] with generalized gradient approximation (GGA) . Since the primary interest of the present study is in the cohesive properties of Cr-Ni-base alloys at elevated temperatures, all calculations have been performed at a relevant high-temperature bcc lattice constant of 2.92 Å . As the magnetic moment of the high-temperature paramagnetic phase of chromium has been shown to be zero within the framework of the disordered local moment model (see for instance ), we have performed only non-spin-polarized calculations. The convergence accuracy of the total energy has been chosen as 10-5 eV and 9 × 10−3 eV/A for the forces. The ionic relaxations have been included in the all types of calculations.
2.2. Semiempirical Calculations
The semiempirical simulations made use of the second nearest neighbor modified embedded atom method (2NN-MEAM) potentials  within the LAMMPS environment [25, 26]. These potentials were developed specifically for bcc metals, for which the contribution of the second nearest neighbors is not to be neglected. The simulation cells were set up with the equilibrium lattice constant of the potential; that is, 0 = 2.88 Å. For the relaxations, the Hessian-free truncated Newton algorithm was selected, and the chosen convergence criterion for the relaxation was relative energy changes between iterations of less than 10−11 and a 2-norm of the global force vector of less than 10-13 eV/Å.
2.3. Grain Boundary and Free Surface Modeling
A special coincident site lattice (CSL) model GB Σ5 (210) and a free surface (FS) (210) have been used in all DFT calculations as a representative GB and the corresponding surface created after the GB cleavage. The GB has been modeled by a supercell containing 40 atomic layers of Cr (one atom per layer) separated by 7 Å of vacuum as shown in Figure 1. The same supercell, but without GB, has been used for the free surface (FS) calculations. Only substitutional sites have been considered for the GB and surface segregation of 4d and 5d impurities. The 8 × 4 × 2 Monkhorst-Pack  -point mesh has been used for GB and FS calculations, whereas the 6 × 6 × 6 mesh has been employed in pure bulk calculations with a 54-atom bcc supercell. The structure of Σ5 (210) GB has been additionally relaxed by shifting two grains in the slab with respect to each other. The found minimum energy structure (shown in Figure 1) has been used in all GB slab calculations.
For the semiempirical GB simulations, the size of the simulation cells was increased due to the low computational cost. Therefore, the size of the grains normal to the GB was set to at least 30 Å for each grain separated, which were separated by 7 Å of vacuum. A large set of different GB structures was built using the atomic simulation environment (ASE) python package , and, analogous to the DFT simulations, the minimum energy structure was identified for each of these GBs by shifting the grains with respect to each other.
3.1. The Concept of Low Alloying Additions
In accordance with the Rice–Thomson–Wang approach [6, 7] the mechanical strength of a solid/solid interface (GB) is determined by the work of separation , which is easily accessible from first principles (ab initio) calculations based on density functional theory (DFT) . Here, we use within the concept of low alloying additions described in detail in [8, 9] to characterize the influence of 4d and 5d transition metals on the cohesive strength of selected alloys. As discussed in [8, 9], the alloying elements will be selected using the following three criteria:(i)a positive partial cohesive energy  where is the cohesive energy of an alloy and is the concentration of an th component. characterizes an influence of impurities on the cohesive energy of the bulk of alloy, and an increase in bulk cohesion is indicated by positive value of .(ii)a negative change in work of separation per impurity atom, that is, the strengthening energy η: where is the work of separation of a GB in pure Cr, is the work of separation of this GB with an impurity, is the impurity excess at the GB, is the impurity segregation energy to the GB, and is the impurity segregation energy to the FS. characterizes the influence of impurities on the GB strength with negative values corresponding to an increase in GB strength.(iii)a negative GB segregation energy defined for a substitutional impurity as where represents the total energy of the GB supercell containing host atoms and impurity atoms. The term is the energy difference between two bulk supercells before and after one of the host atoms is replaced by an impurity atom. In order to achieve a better cancellation of numerical errors originating from the -point sampling and the finite size of the plane-wave basis set, we must distinguish between the energy of the Cr atom obtained from the cubic 54-atom bcc supercell bulk and the same energy calculated in the slab supercell (by completely filling it up with layers of Cr, i.e., the so-called “space filling slab”), . Both and are the total energies of the respective cells divided by the number of atoms. We apply a similar formula for the case of the free surface (FS) substitutional segregation used in (3): where represents the total energy of the FS supercell containing host and impurity atoms.
All parameters entering the formulae above can be determined using ab initio calculations based on DFT. This procedure has already been applied to Ni [8, 9], which is why the focus of ab initio calculations in this paper lays primarily on the investigation of the influence of the series of 4d and 5d transition metals on the cohesive strength of Cr. The data from our previous study of Ni [8, 9] will be used along with the current data on Cr for the composition design of a new Cr-Ni alloy.
3.2. Choice of Representative Grain Boundary and Segregation Site
Parameters and for each impurity atom are calculated using a representative GB and a representative segregation site. To select a suitable GB structure, we have performed a systematic scan of a number of GB geometries in Cr using 2NN-MEAM calculations. The GB energies of these GBs are plotted in Figure 3 and demonstrate that the GB energy is highly anisotropic. While two Σ3 GBs exhibit GB energies below 0.5 J/m2, most GBs display higher GB energies within the range 1–1.5 J/m2. As can be seen in Figure 3, the Σ5 (210) GB already used in  represents a reasonable compromise due to its minimal structure size (40 atoms) and due to the fact that its GB energy falls in the midrange of special GBs formation energies in Cr. There is a large offset between DFT and 2NN-MEAM data in the case of the Σ5 (210) GB, but a similar deviation of 2NN-MEAM data from DFT has been observed in the case of Fe calculations, which can be traced back to a fitting of the 2NN-MEAM potentials to experimental surface energies that are in general below the values found with DFT . However, as it has been shown also in , MEAM and DFT provide qualitatively the same results for Fe. Therefore, 2NN-MEAM results may be used for a relative energy analysis.
The site preference within the selected Σ5 (210) GB has been checked by calculating the segregation energies of a Ta-impurity at G1–G9 GB sites in the vicinity of the GB (see Figure 1). The results show that G5 and G7 sites have almost identical values (within 0.01 eV) of with the absolute minimum at G5 site. In what follows, we will use this GB site (GB5) for and calculations. By doing so, we assume that since all 4d and 5d impurities are larger in size than Cr , they will exhibit the same site preference as one of the largest impurity atoms (Ta). Such trends have already been observed in other systems, for example, for 4d and 5d impurities in Mo and W  or selected solutes in Fe .
4. Theoretical Results
4.1. Pure Cr
In Table 1, we compare our results for the Σ5 (210) GB and the corresponding (210) free surface (FS) energies of pure Cr. A direct comparison to the DFT data from Ossowski et al.  shows good qualitative agreement between two sets of calculations. Ossowski et al. have investigated the same GB and FS geometries. We relate the differences in absolute values for both GB and FS energies to the differences in the exchange correlation functional, the assumed magnetic states, and the relaxation procedures employed in the calculations. Here, we have performed an additional minimum energy search in the GB plane by shifting two grains with respect to each other and scanning the energy landscape for the minimum energy configuration. The results for the FS energy are also compared to the DFT data from Vitos et al.  obtained for a set of different surface planes (excluding (210)) without local atomic relaxations, DFT data from Ossowski and Kiejna for the (210) FS , and available experimental data [37, 38] (not surface facet specific). Taking into account the scatter in FS energies expected for different FS facets from , the FS energy predicted in our calculations can be considered to be in good agreement with the available experimental data (representing an averaged FS energy) as well.
4.2. Impurity in Bulk
First, we analyze a possible cohesive strengthening effect of the impurity atoms in the bulk of Cr using the partial cohesive energy . The calculated values of are shown in Figure 2. The elements from the middle of the transition elements series have the largest strengthening effect. This trend is in agreement with an analogous trend for 4d and 5d impurities in Ni-matrix. The elements improving the bulk cohesive strength in Cr are Nb, Mo, Tc from the 4d and Ta, W, Re, Os, from the 5d series.
4.3. Impurity Segregation
Results for the segregation energies from first principles calculations of 4d and 5d transition elements to the FS (210) and to the Σ5 (210) GB are presented in Figure 2. The trend for the segregation energies with respect to the d-band filling displays some wave-like features, somewhat different from what has been observed for Ni [8, 9]; however, it resembles closely the trend found for d-band filling in Mo and W (site 1 in Figures and in ), which suggests that this trend is related to the d-band filling of the host element. The only experimental data point available for comparison is related to Mo segregation in a Cr-rich Cr-Mo alloy , where the segregation energy of Mo to the FS was found by means of the Auger spectroscopy to be virtually zero (−0.06 eV/at). This result is similar to the result of our calculations for Mo (210) surface where we also see an essentially zero segregation energy with a slight tendency towards antisegregation (see Figure 2).
The strengthening energy η as the difference between segregation to the GB and to the FS is also presented in Figure 2. It follows a trend similar to that of the partial cohesive energy , yet shifted somewhat to right (about one element). According to these results, the elements leading to an increased GB cohesion if present at the GB are Mo, Tc, Ru, Ta, W, Re, Os, and Ir.
Figure 2 shows a correlation between the influence of the alloying elements on the bulk and GB cohesion in Cr: the elements from the middle of the series have the most pronounced cohesive strengthening effect on the cohesion while the elements from the beginning and the end of the series have the most pronounced embrittling effect. Such behavior is similar to the result obtained in [8, 9] for Ni and seems to be a generic feature for alloying with transition elements in dilute limit in general .
4.4. Choice of Low Alloying Additions
We use the three main parameters of the low alloying additions concept , , and [8, 9] presented in Figure 2 to identify the cohesion strengthening elements in Cr. Elements with the most remarkable strengthening effect on both GB and bulk are located in the middle of each element series. If we exclude precious and life hazardous elements from our consideration, we identify Re, W, Mo, and Ta as the elements with a positive effect on both bulk and GB cohesion. Let us also remark that Nb (a beneficial addition to Ni-base alloys) only slightly decreases GB strength in Cr but still strengthens the bulk phase. In the Ni alloys, the most cohesion enhancing elements were identified as Zr, Hf, Nb, and Ta [8, 9]. Consequently, our choice fell on Ta and Nb, which are the only additions with a GB and bulk strengthening effect in γ-phase (Ni), practically neutral for GB in α-phase (Cr), and with some strengthening effect in the bulk of α-phase. So far as Ta and Nb have been absent from the alloying system of the typical Cr-Ni-base alloy (Cr, Ni, W, Ti, and V, alloy I), here, we propose to include them into the chemical composition of the Cr-Ni-base alloy as cohesion strengthening elements to design a new test alloy with improved mechanical properties (alloy II).
5. Experimental Procedure and Results
5.1. Materials and Methods
The ingots (~60 kg) of a test alloy II and alloy I have been melted in a vacuum-induction furnace ISV-0.16; electric slag has been remelted and annealed at 900°C for 16 hours. The obtained alloys have had the following chemical compositions (wt.%): alloy I (Cr-32.7Ni-1.92W-0.3Ti-0.3V); alloy II (Cr-32.2Ni-2.06W-0.08Ti-0.29V) + (0.15Ta, 0.16Nb) . For mechanical testing, the ingots of the alloys have been hot extruded from the rods 100 mm to the bars Ø28 mm (standard technology for production of the industrial alloy I). The X-ray diffraction (XRD) analysis (Siemens D500 X-ray with Cu-Kα radiation) and scanning electron microscopy (SEM, microscope JEOL JSM 6610LV) have been used to investigate the microstructure of the alloys. We have used the Instron 3382 machine for tensile testing of the samples at 20, 900, and 1080°С.
5.2. Microstructure of the Test Alloy II
The XRD data demonstrate that alloys I and II consist of the two main phases: α (bcc) and γ (fcc). Figure 4 shows a typical microstructure of the Cr-Ni-base alloys I and II after melting, electric slag remelting, and annealing. One can observe a formation of the relatively large light areas () that decorate GBs in -matrix and make GBs visible in the microstructure.
A not etched sample of the (Ta, Nb)-bearing alloy, Figure 5 (SEM), was used for investigation of the distribution of the alloying elements in the microstructure. Electron probe microanalysis maps, Figure 5, obtained from micrograph Figure 4 (SEM) demonstrate that the light areas near GBs are enriched with Ni and dark areas in the core of the grains are enriched with Cr. It is interesting to note that a distribution of Ta in the microstructure is fully correlated with a distribution of Ni; that means that Ta is predominantly located in the large γ-areas near GBs. W, Ti, and V are practically homogeneously distributed in the microstructure.
5.3. Mechanical Testing
The results of tensile mechanical testing at different temperatures of alloys I and II are given in Table 2. One can see that both alloys have practically the same tensile strength and elongation δ at room temperature. However, alloy II demonstrates a remarkable improvement of the strength at high temperature 1080°C; namely, the tensile strength increases from 89 to 137 MPa. In addition, we have been able to suppress unwanted superplastic behavior of the alloy I at 1080°C by reducing the elongation δ from 136 to 20%.
The results of creep testing of alloys I and II at temperature 800°C under loading 100 and 140 MPa are presented in Figure 6, which demonstrates a significant increase in the creep resistance of alloy II due to addition of small amounts of Ta and Nb into alloy I.
We have used DFT calculations to investigate the influence of 4d and 5d transition metal alloying elements on the bulk and GB cohesive strength in Cr within the low alloying additions concept based on the ideas of the Rice–Thomson–Wang approach. The results of this investigation have been combined with the analogous results for the Ni to select Ta and Nb as acceptable additions to dual-phase () Cr-Ni-base alloys. Both elements fit the criteria of the low alloying additions concept and strengthen GB and bulk of both γ and α phases best out of all considered 4d and 5d alloying elements.
These new alloying elements have been used for designing a new test Cr-Ni-base alloy II by means of adding (0,15Ta, 0,16Nb) to the known Cr-Ni-base alloy I (Cr-32,2Ni-2,06W-0,08Ti-0,29V) (wt.%). These alloys have been produced by means of melting in a vacuum-induction furnace ISV-0.16, electric slag remelting, and annealing at 900°C for 16 hours.
Investigation of the microstructure of Ta-bearing Cr-Ni-alloy has revealed enrichment of Ta in the large γ-phase areas near GBs in the α-matrix that we consider as a potency to increase simultaneously the cohesive strength of GBs and the cohesive energy of the bulk in γ-phase. Mechanical testing of both alloys has demonstrated that the test alloy II with Ta and Nb has had improved tensile strength and creep resistance at high temperatures with respect to the original alloy I.
The data on DFT segregation and strengthening energies calculations are included in the SEGROcalc database available in https://www.mcl.at/software/segrocalc/.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
V. I. Razumovskiy would like to thank Professor A. V. Ruban and Dr. A. Kartsev for useful discussions of some of the presented theoretical results. Computer resources at the National Supercomputer Center (NSC) in Linköping were provided by the Swedish National Infrastructure for Computing (SNIC). The financial support under the scope of the COMET program within the K2 Center “Integrated Computational Material, Process and Product Engineering (IC-MPPE)” (Project no. 859480) is gratefully acknowledged. This program is supported by the Austrian Federal Ministries for Transport, Innovation and Technology (BMVIT) and for Digital and Economic Affairs (BMDW), represented by the Austrian Research Funding Association (FFG), and the federal states of Styria, Upper Austria, and Tyrol.
- C. T. Sims, N. S. Stoloff, and W. C. Hagel, Superalloys II, Wiley, NY, USA, 1987.
- R. C. Reed, The Superalloys: Fundamentals and Applications, Cambridge University Press, Cambridge, UK, 2006.
- D. Chang, D. Krueger, and R. Sprague, “Superalloy powder processing, properties, and turbine disk applications,” in The Metallurgical Society of AIME, R. H. Bricknell, W. B. Kent, M. Gell, C. S. Kortovich, and J. F. Radavich, Eds., pp. 245–273, Superalloys.
- G. S. Garibov and A. V. Vostrikov, “Current trends of PM superalloys discs production technology for gas turbine engines,” in Proceedings of the 2005 International Conference Hot Isostatic Press, pp. 86–91, Paris, 2005.
- S. Z. Bokshtein, S. S. Ginzburg, S. T. Kishkin, I. M. Razumovskii, and G. B. Stroganov, Autoradiography of Grain and Interphase Boundaries and Structural Stability of Alloys , Moscow, 1987 (Russian).
- J. R. Rice and J.-S. Wang, “Embrittlement of interfaces by solute segregation,” Materials Science and Engineering: A Structural Materials: Properties, Microstructure and Processing, vol. 107, no. C, pp. 23–40, 1989.
- J. R. Rice and R. Thomson, “Ductile versus brittle behaviour of crystals,” Philosophical Magazine, vol. 29, no. 1, pp. 73–97, 1974.
- V. I. Razumovskiy, A. Y. Lozovoi, and I. M. Razumovskii, “First-principles-aided design of a new Ni-base superalloy: Influence of transition metal alloying elements on grain boundary and bulk cohesion,” Acta Materialia, vol. 82, pp. 369–377, 2015.
- V. I. Razumovskiy, A. Y. Lozovoi, and I. M. Razumovskii, “First-principles-aided design of a new Ni-base superalloy: Influence of transition metal alloying elements on grain boundary and bulk cohesion,” Acta Materialia, vol. 106, pp. 401-402, 2016.
- A. G. Beresnev, A. I. Logacheva, A. V Logachev, A. V Logunov, I. M. Razumovskii, and P. V Taran, “Heat-resistant Granular Ni-based Alloy Containing C,” Yt and Ni in Specific Weight Percentage, Patent RU2386714-C1, 2010.
- J. R. Davis, Nickel, Cobalt, and Their Alloys, ASM International, 2000.
- V. N. Butrim, “Technological Aspects of Manufacturing of High-Quality High-Chromium Alloy Semifinished Products for Production of Critical Spacecraft Thruster Components,” Light Alloy Technology, vol. 2, pp. 95–104, 2015 (Russian).
- G. M. Voronin, S. T. Kishkin, and I. O. Panasyuk, “Refractory alloys for aero- and space scpece applications,” in Aviation materials between the 20-21st centuries, G. M. Voronin, S. T. Kishkin, I. O. Panasyuk, V. N. Podyachev, and E. V. Sivakova, Eds., pp. 264–273, Nauka, VIAM, Moscow, 1994.
- A. V. Perepelkin, G. F. Sarzhan, S. A. Firstov, and G. G. Kurdyumova, “Investigation of the deformation and rupture mechanisms in two-phase Cr-Ni alloy VKh-4,” The Physics of Metals and Metallography, vol. 48, pp. 588–592, 1979.
- P. T. Dawson and S. A. Petrone, “Studies on Cr-Mo binary alloys II. Cross-over in surface segregation for chromium-rich alloys,” Surface Science, vol. 304, no. 1-2, pp. 119–130, 1994.
- D. Scheiber, R. Pippan, P. Puschnig, and L. Romaner, “Ab initio calculations of grain boundaries in bcc metals,” Modelling and Simulation in Materials Science and Engineering, vol. 24, no. 3, Article ID 035013, 2016.
- E. M. Savitsky and G. S. Burkhanov, Metallurgy of Refractory and Rare Metal Based Alloys, Moscow, 1971.
- P. E. Blöchl, “Projector augmented-wave method,” Physical Review B: Condensed Matter and Materials Physics, vol. 50, no. 24, pp. 17953–17979, 1994.
- G. Kresse and J. Furthmüller, “Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set,” Physical Review B: Condensed Matter and Materials Physics, vol. 54, no. 16, pp. 11169–11186, 1996.
- G. Kresse and J. Furthmüller, “Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set,” Computational Materials Science, vol. 6, no. 1, pp. 15–50, 1996.
- J. P. Perdew, K. Burke, and M. Ernzerhof, “Generalized gradient approximation made simple,” Physical Review Letters, vol. 78, no. 7, article no 1396, 1997.
- N. A. Dubrovinskaia, L. S. Dubrovinsky, S. K. Saxena, and B. Sundman, “Thermal expansion of chromium (Cr) to melting temperature,” Calphad, vol. 21, no. 4, pp. 497–508, 1997.
- P. A. Korzhavyi, A. V. Ruban, J. Odqvist, J.-O. Nilsson, and B. Johansson, “Electronic structure and effective chemical and magnetic exchange interactions in bcc Fe-Cr alloys,” Physical Review B: Condensed Matter and Materials Physics, vol. 79, no. 5, Article ID 054202, 2009.
- B. Lee, M. Baskes, H. Kim, and Y. Koo Cho, “Second nearest-neighbor modified embedded atom method potentials for bcc transition metals,” Physical Review B: Condensed Matter and Materials Physics, vol. 64, no. 18, Article ID 184102, 2001.
- S. Plimpton, “Fast parallel algorithms for short-range molecular dynamics,” Journal of Computational Physics, vol. 117, no. 1, pp. 1–19, 1995.
- S. Plimpton, “Fast parallel algorithms for short-range molecular dynamics,” Journal of Computational Physics, vol. 117, no. 1, pp. 1–19, 1995, http://lammps.sandia.gov.
- H. J. Monkhorst and J. D. Pack, “Special points for Brillouin-zone integrations,” Physical Review B: Condensed Matter and Materials Physics, vol. 13, no. 12, pp. 5188–5192, 1976.
- S. R. Bahn and K. W. Jacobsen, “An object-oriented scripting interface to a legacy electronic structure code,” Computing in Science & Engineering, vol. 4, no. 3, pp. 56–66, 2002.
- Y. Mishin, M. Asta, and J. Li, “Atomistic modeling of interfaces and their impact on microstructure and properties,” Acta Materialia, vol. 58, no. 4, pp. 1117–1151, 2010.
- I. M. Razumovskii, A. V. Ruban, V. I. Razumovskiy et al., “New generation of Ni-based superalloys designed on the basis of first-principles calculations,” Materials Science and Engineering: A Structural Materials: Properties, Microstructure and Processing, vol. 497, no. 1-2, pp. 18–24, 2008.
- T. Ossowski, E. Wachowicz, and A. Kiejna, “Effect of iron additions on intergranular cohesion in chromium,” Journal of Physics: Condensed Matter, vol. 21, no. 48, Article ID 485002, 2009.
- J. C. Slater, “Atomic radii in crystals,” The Journal of Chemical Physics, vol. 41, no. 10, pp. 3199–3204, 1964.
- D. Scheiber, R. Pippan, P. Puschnig, A. Ruban, and L. Romaner, “Ab-initio search for cohesion-enhancing solute elements at grain boundaries in molybdenum and tungsten,” International Journal of Refractory Metals and Hard Materials, vol. 60, pp. 75–81, 2016.
- H. Jin, I. Elfimov, and M. Militzer, “Study of the interaction of solutes with Σ5 (013) tilt grain boundaries in iron using density-functional theory,” Journal of Applied Physics, vol. 115, no. 9, Article ID 093506, pp. 1–8, 2014.
- L. Vitos, A. V. Ruban, H. L. Skriver, and J. Kollár, “The surface energy of metals,” Surface Science, vol. 411, no. 1-2, pp. 186–202, 1998.
- T. Ossowski and A. Kiejna, “Density functional study of surface properties of chromium,” Surface Science, vol. 602, no. 2, pp. 517–524, 2008.
- F. R. De Boer, Cohesion in Metals, North-Holland, Amsterdam, 1988.
- W. Tyson and W. Miller, “Surface free energies of solid metals: Estimation from liquid surface tension measurements,” Surface Science, vol. 62, no. 1, pp. 267–276, 1977.
- W. T. Geng, A. J. Freeman, and G. B. Olson, “Influence of alloying additions on grain boundary cohesion of transition metals: First-principles determination and its phenomenological extension,” Physical Review B: Condensed Matter and Materials Physics, vol. 63, no. 16, 2001.
- V. N. Butrim, I. M. Razumovskiy, V. N. Kashirtsev, and A. G. Beresnev, “Chromating Alloy and Method of Alloy Melting,” Patent RU2620405-C1, 2016.
Copyright © 2018 V. I. Razumovskiy et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.