Calculation of solubility in titanium alloys from first-principles
We present an approach to calculate the atomic bulk solubility in binary alloys based on the statistical-thermodynamic theory of dilute lattice gas. The model considers all the appropriate ground states of the alloy and results in a simple Arrhenius-type temperature dependence determined by a “low-solubility formation enthalpy”. This quantity, directly obtainable from first-principle calculations, is defined as the composition derivative of the compound formation enthalpy with respect to nearby ground states. We apply the framework and calculate the solubility of the A specie in A-Ti alloys (A=Ag,Au,Cd,Co,Cr,Ir,W,Zn). In addition to determining unknown low-temperature ground states for the eight alloys, we find qualitative agreements with solubility experimental results. The presented formalism, correct in the low-solubility limit, should be considered as an appropriate starting point for determining if more computationally expensive formalisms are otherwise needed.
High-throughput ab initio methods, capable of predicting properties of an ample set of materials from quantum mechanics calculations, are becoming important tools for scientists working in rational materials development. These methods allow researchers to correlate between different systems and to observe trends converging toward the predictions of new materials (2); (3); (4); (5); (6); (7). The main difference between the “several-calculations” and the “high-throughput” philosophies is that the latter requires rapid estimations of materials properties so that the correlations between systems, even if roughly characterized, become the target information instead of the throughout and accurate description of a small subset. Clearly, to analyze the extensive amount of information and to extract correlations, ad-hoc algorithms and appropriate computer softwares have to be developed. Furthermore, once the space of the search is narrowed, a detailed study can be employed on the obtained reduced set of feasible candidate systems.
Several examples have appeared in literature in recent years, for instance the “data-mining of quantum calculations” method leading to the principle component analysis of the formation energies of many alloys in several configurations (5); (6), the evolutionary approach for determining hamiltonian (8), the “Pareto-optimal” alloys and catalysts (9); (10) the prediction of the lithium-boron superconductor (11), the “high-throughput Kohn-anomalies” search in ternary lithium-borides (12); (13), and the “multi-optimization” techniques used in studying high-temperature reactions in multicomponent hydrides (14); (15); (16).
This manuscript focuses on the high-throughput formalism for the calculations of solubility in binary alloys (solvus lines). The knowledge of solubility is crucial for designing new alloys with particular physical, chemical, and mechanical properties. For example, if we have to enhance an alloy property by adding extra specie as solute, it is necessary to know the equilibrium solubility to understand if it is possible to dissolve the candidate specie, and, if possible, to avoid supersaturation-precipitation and subsequent modification of the target property (aging effect). In superconducting materials research, the problem emerges frequently: often expensive and difficult experiments are undertaken to enhance the critical temperature (17); (18). Even in catalysis research, recent experiments and modeling have shown that the solubility of carbon is responsible of thermodynamic instabilities hindering the catalytic activity of very small Fe and Fe:Mo clusters (19); (20). The calculation of solubility of Zr in Al has already been addressed with success within the regular solution model (21) fit to ab initio calculations (22). Solubility can also be extracted from the knowledge of the phase diagram which, in the case of lattice-conserving alloys, can be generated within the the cluster expansion (23); (24) and Monte Carlo approaches(25); (26). However, a straightforward formalism leading to the estimation of equilibrium solubility for general alloys is still lacking.
In the present paper we devised a statistical-thermodynamic approach for the calculation of atomic solubility in alloys. The advantage of our approach consists in taking into account all available ground states rather than just the pure species configurations. To test the method, we present calculations for a number of binary titanium systems. The paper is organized as follows: In Sec. II we rewrite the equations governing solubility in the case of vacancies and substitutional impurities in binary alloys. Section III is devoted to the discussions of capabilities and limits of our formalism. Examples of phase diagrams and solubilities are addressed in Sec. IV for the following test Ti-A systems (A=Ag,Au,Cd,Co,Cr,Ir,W,Zn). Conclusions are given in Sec. V.
Ii Solubility formalism
Let us consider a disordered dilute solution of A-atoms and vacancies (V) in a pure B-solid with a given Bravais crystal lattice (labeled with the “dis”). Without taking into account the A-A, A-V, and V-V interactions and assuming that A and V concentrations are small, the approximate enthalpy of the considered alloy can be written as:
where , and are the enthalpy of the pure B solid per unit cell, the change in enthalpy of the solid upon substitution of one B with an A-atom, and the change in enthalpy upon removal of one B atom, respectively. In addition, and (=A,B,V) are the total numbers of crystal lattice sites and of atoms of -type:
and (=A,B) are energy and the volume per atom of the -pure solid, is the pressure, and represent the change of volume of the B-pure solid upon introduction of one A-atom or one vacancy. The framework introduced by Eq. (1) is similar to Wagner-Schottky model of a system of non-interacting particles.(27). The quantities with superscript “0” are considered to be temperature, pressure, and composition independent being calculated at zero temperature and pressure. Hence, the accuracy of the following results will be better in the case of limited pressures or for systems with very high bulk modulus, where the elastic energy fraction of is negligible.
The quantities in Eq. (2) can be easily approximated as differences of first-principles energies and volumes between large supercells (sc) with or without defects (substitutional A-atom or vacancy):
As the size of the supercell grows, the approximate quantities in Eqs. (3) approach their exact values. In literature(28) and in this paper, , and (=A,V) are called the “raw” (composition unpreserving) -defect formation energy, volume and enthalpy, respectively.
The enthalpy per atom is obtained from Eq. (1) as:
where (=A,B,V) are the atomic concentrations
Then, we follow the convention of using the formation enthalpy calculated with respect to the pure A- and B-solids (29):
where the quantities and , defined as
Equilibrium Gibbs Free Energy
The formation Gibbs free energy, , is defined as:
where is described by Eq. (7) and the formation entropy can be obtained within the mean-field approximation as
where , , (=A,B,V) are the temperature, the Boltzmann constant, and the site concentrations of atoms:
By changing variables from site concentrations to atomic concentrations :
we rewrite Eq. (10) as
In alloy with fixed atomic composition , the equilibrium concentration of vacancies is determined by minimizing the formation Gibbs free energy:
To conclude, the “true” vacancies formation enthalpy determines the equilibrium concentration of vacancies with an Arrhenius-type equation (see also Ref. (29)). In the next section, we show that resolves the solubility in the case of a phase-separating alloy having no intermediate ground states.
At a given temperature, the solubility of A-atoms in a B-solid phase, , is defined as the maximum homogenously achievable concentration of A without the formation of a new phase (Fig. 1(a)). The accurate calculation of requires considering of the nearby ground state (labeled as “gs”) with respect to the increase of , see the “red line” in Fig. 1(b). It is implemented by minimizing the formation Gibbs free energy of a mixture of (1) a disordered dilute solution of A-atoms and vacancies in a B-rich solid phase at composition (the “dis”-phase of the previous section), and (2) an on- or off-stoichiometric ground state “gs”-phase at composition . The lever rule gives the fractions of the two phases:
and the minimization is performed with respect to and ( is the overall composition of A in the two-phase mixture, ):
which approximates as an Arrhenius-type relation at low temperature:
The quantity is defined as
where is the formation enthalpy of the ground state “gs”. In the low-solubility limit, becomes:
where is the non-vacancy contribution (labeled as “nv”):
where the two contributions, non-vacancy and vacancy, are factorized. The last expression indicates that the presence of vacancies effectively increases the solubility by decreasing the number of host B-atoms in the solution.
Iii Interpretation of
For low-solubility calculations of non-interacting defects, the framework can be implemented through first-principles calculation of (Eq. (24)). It requires the knowledge of the enthalpy and composition of the ground-state “gs”, as well as the “true” defect formation enthalpies and (Eq. (8)).
To capture the physical meaning of , let us consider an arbitrary dilute solution “sc” at composition (see Fig. 2). The label “sc” indicates that the solution is generated as a supercell of the B-solid upon insertion of defects, randomly distributed but not too close to avoid interactions. The first part of Eq. (25) can be rewritten according to Eq. (9) as:
The second part of Eq. (25) becomes:
where is the formation enthalpy of the mixture of the pure B-solid with the ground state “gs” with overall composition (point (B+“gs”) in Fig. 2). Thus, we obtain:
The comparison of Eq. (27) with Eq. (29) leads to the conclusion that both and are derivatives of supercell formation energies with respect to A-composition. For , the supercell formation enthalpy is determined with respect to pure A and B phases (the distance between points “sc” and (B+A) in Fig. 2). For , the supercell formation enthalpy is determined with respect to B-pure and the ground state “gs” (the distance between points “sc” and (B+gs) in Fig. 2). To conclude, and are characterized by the angles and between the B-“sc”/B-A and B-“sc”/B-“gs” lines, respectively. In analogy with the definition in Eq. (8), can be considered as the “low-solubility formation enthalpy”.
The quantities and are identical only for phase-separating alloys having no intermediate ground states (i.e. “gs” A). In this case the low-solubility can be formally determined by minimizing (Eq. (14)) with respect to (e.g. Ref. (30)). However, in the general case, the existence of ordered ground-states must be verified so the appropriate formalism is used. Generally, in ordering alloys and differ, and might even have different signs.
Note that if the “sc” point is below (B+gs) in Fig. 2 () the solubility expressions (20-21) are not valid, and there must exist an undetected ground state (it might be “sc” itself) with concentration lower than . In this case, such undetected ground state should be used for the calculation of solubility rather than initial “gs” (31).
The expression for non-binary low-solubility within the regular solution model derived in Refs. [(21); (22)] coincides with our derivation in the case of binary alloys without vacancies and high-temperature contributions. This is because the regular solution model corresponds to our model for the free energy in case of dilute solution.
Iv Ground states of selected Titanium alloys
As an example of our formalism, we calculate the solubility of a set of metals in titanium. First, we explore the possible ground states of the A-Ti systems (A=Ag,Au,Cd,Co,Cr,Ir,W,Zn) and then we apply the construction described in the previous section.
The low temperature stability of A-Ti is performed by using our high-throughput quantum calculations framework (6); (5); (11); (13), based on first-principles energies obtained with the VASP software (32). We use projector augmented waves (PAW) pseudopotentials (33) and exchange-correlation functionals as parameterized by Perdew-Burke-Ernzerhof (34) for the generalized gradient approximation (GGA). Simulations are carried out with spin polarization except Ti-Co, at zero temperature and pressure (35), and without zero-point motion. All structures are fully relaxed (shape and volume of the cell and internal positions of the atoms). The effect of lattice vibrations is omitted. Numerical convergence to within about 1 meV/atom is ensured by enforcing a high energy cut-off (357 eV) and dense 6,000 k-point meshes.
The number of crystal structures considered for the calculations of each A-Ti system is 194. In addition to the 176 described in Ref. (6), we included the following prototypes A5, A6, A7, A11, CaGe, NbNi (PtTi), VZn, C36 and the whole complete set of hcp-superstructures with up to 4 atoms/cell. The whole process is performed in an automatic fashion through the software AFLOW which generates the prototypes, optimize the parameters, perform the calculations, correct possible errors, and calculate the phase diagrams (6); (36).
We propose degenerate results for: 1) CoTi, experimentally reported as C15, but with ab initio formation energies of -311.4 meV/at. and -304.0 meV/at. for C14 and C15, respectively; 2) CoTi, experimentally reported as NiTi, but with ab initio formation energies of -291.2 meV/at., -287.3 meV/at., and -285.7 meV/at. for C37, CuZr, and NiTi, respectively; 3) TiZn, experimentally reported as B2, but with ab initio formation energies of -195.4 meV/at. and 193.4 meV/at. for L1 and B2, respectively.
We propose the novel results for: 1) IrTi, experimentally reported as a two phases region above 500C, but with a possible ab inito low temperature ground state CaGe with formation energy of -373.8 meV/at. 2) IrTi, experimentally reported as a two phases region above 500C, but with a possible ab initio low temperature ground state C11, with formation energy of -716.0 meV/at. 3) TiW, experimentally reported as a two phases region above 500C, but with a possible ab initio low temperature ground state BCC, with formation energy of -82.7 meV/at. (see notation for the prototype in Ref. (6)). 4) TiW, experimentally reported as a two phases region above 500C, but with a possible ab initio low temperature ground state D1, with formation energy of -83.8 meV/at. 5) TiZn, experimentally not explored, but with a possible ab initio low temperature ground state A15, with formation energy of -120.0 meV/at.
In particular, the results indicate that in the Ir-Ti system the low temperature Ir-rich part of the known phase diagram is not complete, and that the Ti-W alloy has an ordering tendency at low temperature, in contrast to common belief (37); (38). In addition, in Ti-Zn there must exist a Ti-rich compound with Ti composition higher than the reported TiZn, (37); (38).
|(Refs. (38); (37))||result||(meV/at.)||group(39)|
V Results and discussion of solubility in Ti alloys
The “raw” formation enthalpies of Ti alloys with one substitutional atom or a vacancy were obtained through Eqs. (3) considering supercells of the hcp Ti. Supercell dimensions were chosen to limit the defect-defect interactions, being their distance at least three times larger than the nearest neighbor Ti-Ti bond. The results are presented in Table 2. Solubilities temperature dependencies are presented in Fig. 11, while Fig. 12 illustrates a comparison of experimental and theoretical data at T=700C.
Table 2 shows negative “true” formation enthalpies for A=Zn, Au, and Ir. This indicates that calculation of Zn, Au, and Ir solubilities in Ti is not faceable without considering nearby intermetallic ground states.
Figures 11-12 show that the highest theoretical solubilities in Ti occurs for Zn, Cd, Ag, and Au as consequence of their low solubility formation enthalpies (see Table 2). This high solubility has also been observed experimentally for Cd, Ag, and Au, whereas, to our best knowledge, solubility of Zn does not seem to have been studied. High solubility of late transition metals Zn, Cd, Ag, and Au in the early transition metal titanium is due to the substantial localized stability provided by the filling tendency of the -band of Ti (Ref. (40)).
Because of the very high formation enthalpy of vacancies in Ti (reported in Table 2), the vacancy equilibrium concentration was found to be very low at all considered temperatures (e.g. at C). Thus, the effect of vacancies on the solubilities of Ti alloys should be negligible.
Although, theoretical and experimental results follow similar trends, theoretical solubilities are considerably smaller for most of the considered alloys. A similar discrepancy was also observed for Al-Zr in Ref. (22). The discrepancy could be due to shortcomings of theory and/or experiment. The main approximations of our model consist of (a) neglecting the defect interaction, (b) neglecting the spatial defect correlation and (c) assuming low concentration of defects. However, as such approximations are somehow related, subsequent solubility calculations suggesting low values validate the assumptions (the mean-field approximation neglects the interatomic positional correlations but it should work well when the deviation from complete stoichiometric (pure Ti-solid in our case) is small - and so is the solubility in our case - see Sec. 19 in Ref. (41).) For all considered alloys except Ti-Zn, we found that at intermediate temperatures the calculated solubilities are small enough for our approximations to be valid, although smaller than experiment, as mentioned before. In Ref.(22), the authors added the defect interactions through the Cluster Expansion method but without increasing the solubility considerably. Actually, our formalism includes the solute-solvent ordering tendency by considering the real intermetallic ground state other then the pure Ti-solid. Thus, we conclude that our approximations (or those of the cluster expansion method) are not responsible for the theory-experiment solubility discrepancy.
In case of Ti-Zn, the formation enthalpy is very low (see Table 2). Correspondingly, the theoretically solubility at intermediate temperatures is very high, violating the assumptions of the model. Thus, for Zn in Ti, a more precise solute interaction and correlation parameterizations is required, which can be obtained by using, for example, Cluster Expansion and Monte Carlo simulations. Hence, the presented formalism should be considered as an appropriate starting point to determine if more computationally expensive formalisms are needed.
The other approximation of our model is the assumed independence of our model energy and volume parameters on temperature. The dependence can be caused by non configurational degrees of freedom, like vibrational or anharmonicity (magnetic ordering is not actual for considered alloys). However, theory-experiment solubility discrepancy is observed even at low enough temperatures (e.g. 374-385 K (42)), where the vibrational contribution to the free energy is not important. In fact, vibrational entropy is substantially smaller than the configurational contribution (43), so its inclusion can not modify solubility results much. In fact, even in Ref. (22), the authors added phonon contribution without increasing the solubility considerably.
On the other side, the experimental equilibrium solubility tends usually to be overestimated. In fact, the formation of metastable and/or unstable states which are subsequently frozen at low temperatures, can make solubility measurements very challenging. In such scenarios, the measured solubility may correspond to spinodal concentration rather then actual binodal concentration or simply characterize the frozen out of equilibrium solubility remaining from the initial specimen preparation at higher temperature. Besides, the segregation of defects into grain-boundaries, especially in multicrystalline samples prepared through non optimal cooling dramatically affect the amount of frozen defects and solutes.
Based on the statistical-thermodynamic theory of dilute lattice gas, we developed an approach to calculate the atomic bulk solubility in alloys. The advantage of our approach consists in considering all the appropriate ground states rather than the pure species. It was shown that the low-solubility follows simple Arrhenius-type temperature dependence determined by a “low-solubility formation enthalpy”. This quantity is defined as the composition derivative of the compound formation enthalpy with respect to nearby ground states. “Low-solubility formation enthalpy” coincides with the usual defect formation enthalpy only in the case of a phase-separating alloy having no intermediate ground states and vacancies. The key quantities of our model can be directly obtained by first-principles calculations of by fitting experimental temperature solubility dependence. Generalization of our model to intermediate phases and/or to multicomponent, multisublattice, interstitial-substitutional alloys is straightforward.
As examples, we applied the framework for a set ot eight Ti alloys A-Ti (A=Ag,Au,Cd,Co,Cr,Ir,W,Zn). We have found that the highest solubility for Zn, Cd, Ag, and Au is in qualitative agreement with available experimental data and band structure expectations. The quantitative differences between the theory and experiment observed in the present and other similar studies are discussed.
In conclusion, our formalism is correct in the limit of low-solubility and should be considered as an appropriate starting point for determining if more computationally expensive formalisms are otherwise needed.
We acknowledge Wahyu Setyawan, Mike Mehl, and Ohad Levy for fruitful discussions. This research was supported by ONR (Grant No. N00014-07-1-0878) and NSF (Grant No. DMR-0639822) We thank the Teragrid Partnership (Texas Advanced Computing Center, TACC) for computational support (MCA-07S005).
- Electronic address: email@example.com
- Ceder G, Chiang YM, Sadoway DR, Aydinol MK, Jang YI, B. Huang, Nature 1998;392:694.
- Johannesson GH, Bligaard T, Ruban AV, Skriver HL, Jacobsen KW, Nørskov JK, Phys. Rev. Lett. 2002;88:255506.
- Stucke DP, Crespi VH, Nano Lett. 2003;3:1183.
- Curtarolo S, Morgan D, Persson K, Rodgers J, Ceder G, Phys. Rev. Lett. 2003;91:135503.
- Curtarolo S, Morgan D, Ceder G, Calphad 2005;29:163.
- Fischer C, Tibbetts K, Morgan D, Ceder G, Nature Mater. 2006;5:641.
- Hart GL, Blum V, Walorski MJ, Zunger A, Nature Materials 2005;4:391.
- Bligaard T, Johannesson GH, Ruban AV, Skriver HL, Jacobsen KW, Nørskov JK, Appl. Phys. Lett. 2003;83:4527.
- Persson M, Bligaard T, Kustov A, Larsen KE, Greeley J, Johannessen T, Christensen CH, Nørskov JK, J. Catal. 2006;239:501.
- Kolmogorov AN, Curtarolo S, Phys. Rev. B 2006;73:180501(R).
- Calandra M, Kolmogorov AN, Curtarolo S, Phys. Rev. B 2007;75:144506.
- Kolmogorov AN, Calandra M, Curtarolo S, Phys. Rev. B 2008;78:094520.
- Wolverton C, Siegel DJ, Akbarzadeh AR, Ozoliņš V, J. Phys.: Condens. Matter 2008;20:064228.
- Siegel DJ, Wolverton C, Ozoliņš V, Phys. Rev. B 2007;76:134102.
- Akbarzadeh AR, Ozoliņš V, Wolverton C, Advanced Materials 2007;19:3233.
- Cava RJ, Zandbergen HW, Inumaru K, Physica C 2003;385:8.
- Singh RK, Shen Y, Gandikota R, Wright D, Carvalho C, Rowell JM, Newman N, Supercond. Sci. Technol. 2008;21:025012.
- Harutyunyan AR, Awasthi N, Mora E, Tokune T, Jiang A, Setyawan W, Bolton K, Curtarolo S, Phys. Rev. Lett. 2008;100:195502.
- Curtarolo S, Awasthi N, Setyawan W, Jiang A, Bolton K, Harutyunyan AR, Phys. Rev. B 2008;78:054105.
- Sigli C, Maenner L, Sztur C, Shahni R, in Proceedings of the 6th International Conference on Aluminum Alloys (ICAA-6), edited by T. Sato et al. (Japan Institute of Light Metals, 1998), Vol. 1, p. 87
- Clouet E, Sanchez JM, Sigli C, Phys. Rev. B 2002;65:094105.
- de Fontaine D, in Solid State Physics, edited by H. Ehrenreich and D. Turnbull (Academic Press, New York, 1994), Vol. 47, pp. 33â176.
- Ceder G, Comput. Mater. Sci. 1993;1:144.
- Xu Q, Van der Ven A, Phys. Rev. B 2007;76:064207.
- Van der Ven A, Ceder G, Phys. Rev. B 2005;71:054102.
- Wagner C, Schottky W, Z. Physik. Chem. B 1930;11:163.
- Mishin Y, Farkas D, Phil. Mag. A 1997;75:169. Y. Mishin C. Herzig, Acta Mater. 2000;48:589.
- Korzhavyi PA, Ruban AV, Lozovoi AY, Vekilov YK, Abrikosov IA, Johansson B, Phys. Rev. B 2000;61:6003.
- Matysina ZA, Milyan MI, Zaginaichenko SY, J. Phys. Chem. Solids 1988:49;737.
- The existence of such unknown ground state was predicited upon Be substitition of Mg in MgB in Chepulskii RV, Mazin II, Curtarolo S, First-principles search for potential high temperature superconductors in the Mg-B-A (A=alkali and alkaline earth metals) system with high boron content, submitted (2009).
- Kresse G, Hafner J, Phys. Rev. B 1993;47:558.
- Blochl PE, Phys. Rev. B 1994;50:17953.
- Perdew JP, Burke K, Ernzerhof M, Phys. Rev. Lett. 1996;77:3865.
- So that by calling enthalpy we actually mean energy in this section.
- Curtarolo S, Aflow: a computational software to calculate properties of materials in a high-throughput fashion. http://materials.pratt.duke.edu/aflow.html
- Villars P, Cenzual K, Daams JLC, Hulliger F, Massalski TB, Okamoto H, Osaki K, Prince A, Iwata S, Crystal Impact, Pauling File. Inorganic Materials Database and Design System, Binaries Edition, ASM International, Metal Park, OH (2003)
- Binary Alloy Phase Diagrams, edited by T. B. Massalski (ASM International, Metals Park, OH, 1992)
- International Tables for Crystallography, Vol. A, edited by T. Hahn (D. Reidel, Dordrecht, 1983).
- J. Friedel, in Physics of Metals, edited by J. M. Ziman (Cambridge University Press, Cambridge, 1969), p. 495.
- Krivoglaz MA, Smirnov AA, The Theory of Order-Disorder in Alloys (Macdonald, London, 1964).
- Chen Q, Sundman B, Acta Mater. 2001;49:947.
- van de Walle A, Ceder G, Rev. Mod. Phys. 2002;74:11.