ZenGen, a tool to generate ordered configurations for systematic first-principles calculations: The Cr–Mo–Ni–Re system as a case study☆
Graphical abstract
Introduction
The field of thermodynamic modeling has been recently stimulated by the progress of techniques allowing the calculation of thermodynamic quantities from first-principles calculations, such as the Density Functional Theory (DFT) [1]. These methods allow the estimation of formation enthalpies of fully ordered compounds, taking into account their crystal structures. These calculations can be done not only for stable compounds, but also for metastable ones which play an important role in the description of these phases within the Compound Energy Formalism (CEF) [2], [3]. By using the CEF, any intermetallic phase could be described by a sublattice model for which every ordered configuration heat of formation has to be calculated. As an example, a binary phase with five crystal sites described in a 5-sublattice model generates different ordered configurations, a ternary a huge number, but which can be calculated with today's super-computers.
Technically, performing calculations on a large number of end-members may cause two types of problems: (i) a mistake in the distribution of atoms among all different sites; (ii) a too fast relaxation of crystal structure, thus losing the initial symmetry. To avoid these problematics, “ZenGen” code was created. This code is able to generate the necessary input files required to start DFT calculations of the corresponding ordered configurations for a given system. It has been tested on several phases, such as Laves phases (C14, C15, etc.), or other topologically close packed phases (, , D8b, P, δ, etc.). It can also be used to run Special Quasi-random Structures (SQS) calculations [4]. A basic introduction of ZenGen workflow is given is Section 2.
Thereafter, in order to illustrate ZenGen capacity, we have investigated the challenging quaternary Cr–Mo–Ni–Re system. Our aim was not to assess thermodynamically this system, but rather to show that systematic DFT calculations can be run contently in this very complex system, that they allow the calculation of a preliminary ab initio computed phase diagram, and that they can be used as an input for a traditional Calphad assessment. We have demonstrated this approach in our previous works [5], [6]. The results are presented in Section 3.
Section snippets
The ZenGen workflow
“ZenGen” is a free and open source code, governed by the CeCILL-B license under French law [7], which is officially recognized by Open Source Initiative (OSI). It can be downloaded from http://zengen.cnrs.fr. ZenGen can be installed on Unix–Linux machines and uses Bash, Perl and Python languages. In its present version, it is compatible with VASP program [8], [9] for DFT calculations; however it can be adapted to be compatible with other first-principles codes, as well.
φ, the conventional name
Known phase diagrams
The quaternary Cr–Mo–Ni–Re system is well suited for an automatic investigation of phase stability. Solid solutions with three different structures exist: (Ni), (Cr, Mo) and (Re), and two of elements are magnetic. The σ-phase is stable in the four ternaries. This phase is a hard brittle intermetallic compound and has a deleterious effect on the mechanical properties of many technologically important systems, like Cr–Mo–Ni–Re, a key system for Ni-based superalloys. The
Conclusions
ZenGen is a tool to automate the construction of the input files required for systematic DFT calculations of all the ordered configurations of a multicomponent phase. This approach can be used in parallel with a thermodynamic assessment by the Calphad method. It has been illustrated through the calculation of the quaternary Cr–Mo–Ni–Re system for which all the binary solid solutions and every end-member of the quaternary σ-phase have been calculated for the first time.
Many improvements will be
Acknowledgments
DFT calculations were performed using HPC resources from GENCI-CINES (Grant 2014-96175). Financial support, from the Agence Nationale de la Recherche (ANR), Project Armide ANR-10-BLAN-912-01 is acknowledged. N. Bourgeois is thankful to the French ANR agency for financial support from the National Program Investments for Future ANR-11-LABX-022-01. The authors wish to thank Nathalie Dupin for useful and fruitful discussions.
References (34)
- et al.
A regular solution model for phases with several components and sublattices, suitable for computer applications
J. Phys. Chem. Solids
(1981) - et al.
A compound-energy model of ordering in a phase with sites of different coordination numbers
Acta Metall.
(1986) - et al.
Phase stability in the ternary Re–W–Zr system
Acta Mater.
(2014) - et al.
The thermo-calc databank system
Calphad
(1985) Crystal structure and stability of complex precipitate phases in Al–Cu–Mg(Si) and Al–Zn–Mg alloys
Acta Mater.
(2001)- et al.
Determination of phase equilibria in the system Ni–V–Cr–Mo–Re at 1425 K using the graph method
J. Alloys Compd.
(1999) Crystal chemistry and calphad modeling of the σ phase
Prog. Mater. Sci.
(2008)- et al.
Ab initio ternary σ-phase diagram: the Cr–Mo–Re system
Calphad
(2010) - et al.
Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set
Comput. Mater. Sci.
(1996) - et al.
TCP phase predictions in Ni-based superalloysstructure maps revisited
Acta Mater.
(2011)
Ab initio and thermodynamic study of the Cr–Re system
Calphad
Modeling of Ni–Cr–Mo based alloysPart I—Phase stability
Calphad
Modeling of molar volume of the sigma phase involving transition elements
Comput. Mater. Sci.
Phase stability and mechanical properties of sigma phase in Co–Mo system by first principles calculations
Comput. Mater. Sci.
Experimental investigation of the Mo–Ni–Re system
J. Alloys Compd.
Nobel Lectures, Chemistry 1996–2000
Special quasirandom structures
Phys. Rev. Lett.
Cited by (43)
Thermodynamic assessment of the Ni–Co-Re system and diffusion study of its fcc phase
2024, Calphad: Computer Coupling of Phase Diagrams and ThermochemistryThermodynamic modeling of the Mo–Pt system based on the NACEF approach
2024, Calphad: Computer Coupling of Phase Diagrams and ThermochemistryAl–Ni–Ti thermodynamic database from first-principles calculations
2024, Calphad: Computer Coupling of Phase Diagrams and ThermochemistryApplying enhanced active learning to predict formation energy
2024, Computational Materials ScienceThermodynamic reassessment of Fe–Nb–V system
2023, Calphad: Computer Coupling of Phase Diagrams and Thermochemistry
- ☆
Fully documented manual and program are available on http://zengen.cnrs.fr