Alloy Gene Sequence Project and System Science Philosophy

Youqing Xie^1,2,3, Xiaobo Li⁴, Xinbi Liu^1,2,3, Hongjian Peng⁵, Yaozhuang Nie⁶

¹School of Materials Science and Engineering, Central South University, Changsha, China.
²Powder Metallurgy Research Institute, Central South University, Changsha, China.
³State Key Laboratory of Powder Metallurgy, Central South University, Changsha, China.
⁴School of Materials Science and Engineering, Xiangtan University, Xiangtan, China.
⁵School of Chemistry and Chemical Engineering, Central South University, Changsha, China.
⁶School of Physical Science and Technique, Central South University, Changsha, China.

DOI: 10.4236/ijcns.2017.1012019                

Abstract

System Science Philosophy is a knowledge system constructed of universal principle and law sequences. Alloy gene is a characteristic atom existing in the center of coordination cluster and occupying the lattice point of a lattice cell, and carries holographic information about valence electron structure, physical and thermodynamic properties obtained by alloy gene theory. Alloy gene potential energy curve function has developed traditional atom pair interaction potential functions into many atoms’ interaction function associated with valence electron structure, bond length, bond energy, which makes alloy gene Gibbs energy function can be established. Alloy gene Gibbs energy transmissive function has developed traditional partition function. Its equilibrium and sub-equilibrium transmissive modes produce alloy holographic network positioning bank, which is operable platform to achieve transformation from “trial and error” method to the “whole obtained from a few part” law for getting advanced alloys. It has become possible to launch alloy gene sequence project.

Keywords

Alloy Gene Sequence, System Science Philosophy, Alloy Gene Theory, Alloy Gene Gibbs Energy Transmissive Function, Alloy Holographic Network Positioning Bank

1. Introduction

Materials are always one of the levers to propel the human society forward. There are 81 metal elements accounting for 79% in the Element Periodic Table, of which the 3240 kinds of binary alloy systems, 85,320 kinds of ternary alloy systems, and so on, can be composed. In these alloy systems there are a large number of potential advanced alloys to be discovered. Early in the 1937, Slater said that “the advance of physical metallurgy in the next few years is likely to be almost as dependent on the electron theory of metals as the advances of the last generation were dependent on thermodynamics” [1] . Nowadays, the biology has reached the level to be capable of transferring genes and “clone” creature, while in the same old and historical field of metal materials, the design of advanced alloys is still at the “trial and error” method stage. It is attributed to questionable thinking mode, theories and understanding of some experimental data in the QMAC (quantum mechanical ab initio calculations) electron theory [2] - [10] and CALPHAD (calculation of phase diagrams) thermodynamics [11] - [18] communities, and that there are very few scholars who regard seeking Alloy Gene (AG) as the first condition for realizing scientific design of advanced alloys. The System Science Philosophy (SSP) is a knowledge system constructed of universal principle and law sequences, which are integrated from theories and methods in biological, non-biological and societal systems. The SSP functions in the scientific practice are not only to be discovery, explanation and guide but also to strengthen scientific faith and courage of explorers. In order to quickly and efficiently discover, design and produce advanced alloys, we have discovered AG, established Systematic Metal Materials (SMM) science and created Alloy Holographic Network Positioning (AHNP) methods/techniques road by the SSP guide [19] - [24] . With these advances, it has become possible to launch the Alloy Gene Sequence Project (AGSP), of which the major goals are: Achieving transformation of the theoretic basis from Pauling’s chemical bond theory [25] [26] [27] and Lupis’s “central atoms” theory [28] to the SMM science, which drives further developments of the AG theory [29] - [35] , Alloy Phase (AP) theory [36] - [45] and Alloy Systems (AS) theory [46] [47] [48] , then provides theoretical support for the AGSP; Achieving transformation of technical law from “trial and error” method to “Whole Obtained from a few Parts (WOP)” law, which creates the AHNP techniques, including the methods of the AG holographic warehouse and AHNP warehouses of alloys, phases and alloy systems, as well as techniques of discovery, design, certification, manufacturing and deployment of advanced alloys. These techniques will significantly decrease drain on time, manpower, material, financial and energy resources, then provide technical support for the AGSP. The AGSP will create a new era of the SMM science and AG engineering, during the 21^st century. This paper includes common description of theories chain in Text, explanation of methods chain and equations chain in Supplementary Information (1) and (2), respectively.

2. SMM Science

Creative SSP principles will evolve real systematic science theories! It can be explained by the SMM science and taking fcc-based lattice Au-Cu system with Au₃Cu-, AuCu- and AuCu₃-type ordered and disordered phases as an example, which has been considered to be a classical paragon studied nearly by almost alloy theories for one century.

2.1. SMM Science Framework

According to the SSP-systematic globality principle, i.e., “The globality of a complex system is attributed to the multilevel of structures and properties and to the opening loop correlativity of structures, properties and environment.” The SMM science has been established (Figure 1). It includes three level theories: AG theory for determining AG-basic atom state (BAS) sequences and AG sequences, and establishing AG holographic warehouse; Equilibrium and Subequilibrium Alloy Phase (EAP, SAP) theories for determining EAP and SAP sequences, establishing EAP-EAHNP and SAP-SAHNP warehouses; Equilibrium and Subequilibrium Alloy System (EAS, SAS) theories for establishing EAS-EAHNP and SAS- SAHNP warehouses, based on EAP and SAP sequences, respectively.

The EAP and SAP theories are connected, respectively, with AG theory through the equilibrium Minimal Gibbs Energy Path (MGEP) and subequilibrium AG-Tracking Path (AG-TP) modes, based on AG-Gibbs Energy Transmissive (AG-GET) function, and associated with the EAS and SAS theories by Equilibrium Lever Numerical (ELN) and Subequilibrium Crosspoint Numerical

Figure 1. SMM science framework. It consists of alloy gene (AG) theory, equilibrium and subequilibrium alloy phase (EAP, SAP) theories, equilibrium and subequilibrium alloy system (EAS, SAS) theories, as well as AG-Gibbs energy transmissive (AG-GET) function, equilibrium minimal Gibbs energy path (MGEP) and subequilibrium alloy gene-tracking path (AG-TP) modes, and equilibrium lever numerical (ELN) and subequilibrium crosspoint numerical (SCN) methods.

(SCN) methods of isothermal Gibbs Energy (GE) curves, respectively. The SMM science is an integrated theory, in which the valence electron arranging structures and holographic properties in micro-AG theory, AG arranging structures and holographic properties in submicro-AP theory, AP arranging structures and holographic properties in macro-AS theory and environments have been linked up.

2.2. AG Sequences

The key step for establishing SMM science is to seek the AG sequences. According to the SSP-systematic genes (SG) universality principle, i.e., “Any system with appearing diversity of structures and properties possesses oneself of SG sequences.” We have discovered that, in the alloy system with gathering many atoms of some component elements, “a single AG is a characteristic atom existed in the center of the coordination cluster and occupied on the lattice point of a lattice cell.” In the fcc-based lattice Au-Cu system, the AG-sequences are the central $A_i^{Au}$ - and $A_i^{Cu}$ -sequences respectively in the basic coordination clusters $B_i^{Au}$ - and $B_i^{Cu}$ -sequences (Supplementary Information (1)), that are attributed to split of constituent atom states, under influence of coordinative configurations. The $B_i^{Au}=A_i^{Au}\cdot \left[\left(I-i\right)Au+iCu\right]$ and $B_i^{Cu}=A_i^{Cu}\cdot \left[\left(I-i\right)Au+iCu\right]$ coordinative clusters consist of the central $A_i^{Au}$ and $A_i^{Cu}$ characteristic atoms respectively, and the nearest coordinative configuration $\left[\left(I-i\right)Au+iCu\right]$ . Here, I is the coordinative number and equal to 12; i is the number of Cu-atoms and can change from 0 to 12; $\left(I-i\right)$ is the number of Au-atoms. The diversity of alloy phase structures is attributed to combination and arrangement of alloy genes in the AG sequences. The diversity of alloy phase properties is attributed to AG properties and their transmissive modes.

2.3. AG Theory

AG theory includes AG separated theory and AG Characteristic Crystal (AGC) theory consisting of AGC-valence bond theory and AGC-thermodynamics. According to separable ${\epsilon }_i$ and $v_i$ functions in the AG separated theory, the ${\epsilon }_i^{Au}$ , ${\epsilon }_i^{Cu}$ , $v_i^{Au}$ and $v_i^{Cu}$ sequences of the Au-Cu system have been separated from experimental enthalpies and volumes of few alloys by WOP law (Supplementary Information (1) 1.1).

In order to determine AG holographic information about AG valence electron structures, physical and thermodynamic properties, and according to SSP combining rule of imaginal thinking and mathematic thinking, we have proposed the AGC concept $C\left(T\right)$ , established AGC mixed model $C\left(x,T,\sigma \right)$ and additive law $Q\left(x,T,\sigma \right)$ of AGC extensive properties ( $q_i^A$ , $q_i^B$ ) of alloy phases as functions of composition x, temperature $T$ and order degree $\sigma$ :

$C\left(x,T,\sigma \right)={\displaystyle \underset{i=0}{\overset{I}{\sum }}{x}_i^A\left(x,T,\sigma \right)C_i^A}\left(T\right)+{\displaystyle \underset{i=0}{\overset{I}{\sum }}{x}_i^B\left(x,T,\sigma \right)C_i^B}\left(T\right)$ (1)

$Q\left(x,T,\sigma \right)={\displaystyle \underset{i=0}{\overset{I}{\sum }}{x}_i^A\left(x,T,\sigma \right)q_i^A\left(T\right)}+{\displaystyle \underset{i=0}{\overset{I}{\sum }}{x}_i^B\left(x,T,\sigma \right)q_i^B\left(T\right)}$ (2)

where, $x_i$ denotes AGC concentration. Their magical effects are reflected in the following areas:

1) The AGC is a fictitious concept, which consists of the same alloy genes. Therefore, there are $C_i^{Au}$ and $C_i^{Cu}$ sequences consisted with $A_i^{Au}$ and $A_i^{Cu}$ sequences in the fcc-based lattice Au-Cu system [29] [30] . Among them, the primary $C_0^{Au}$ and $C_{12}^{Cu}$ denote respectively the pure Au and Cu metals (Supplementary Information (1)). Then, the difficult problem to obtain AG holographic information may be resolved by the theory of single substances with known ${\epsilon }_i$ and $v_i$ sequences, i.e., AGC theory [31] [32] [33] [34] [35] .

2) The AGC mixed model may be imagined to be that when primary $C_0^{Au}$ and $C_{12}^{Cu}$ characteristic crystals form an alloy, they are split into the $C_i^{Au}$ and $C_i^{Cu}$ sequences at first, then to be mixed into the alloy [36] [37] [38] .

3) The AGC holographic information may be obtained by AGC theory. Each AG holographic information may be represented by corresponding AGC holographic information (Supplementary Information (1) 1.2-1.4).

2.4. AG ψ_i-Sequences

The diversity of AG valence electron structures ( ${\psi }_i$ ) is attributed to combination and arrangement of quantum valence electrons, that may be described by hybridization of some basic atom states in the basic atom state sequences (Supplementary Information (1) 1.2). In the AGC valence bound theory of the AGC theory, the AGC lattice constant equation $a_i\left(\psi ,r\right)$ and AGC potential energy equation $E_i\left(\psi ,r,\epsilon \right)$ as functions of valence electron structure $\psi$ , bond length $r$ and bond energy $\epsilon$ may be respectively used to translate the AG $v_i$ sequences into AGC $a_i\left(\psi ,r\right)$ sequences, and the AG ${\epsilon }_i$ sequences into AGC $E_i\left(\psi ,r,\epsilon \right)$ sequences. The AGC ${\psi }_i$ sequences may be calculated, through the crosspoint method of the equilattice constant $a_i\left(\psi ,r\right)$ -line and equipotential energy $E_i\left(\psi ,r,\epsilon \right)$ -line (Supplementary Information (1) 1.4) the essence of crystal symmetry and properties of AGC sequences may be explained by AGC ${\psi }_i$ sequences [39] [45] .

2.5. AG G_i-Sequences

In the AGC thermodynamics, AG Gibbs energy is the holographic energy state with multi-level thermodynamic properties. It may be divided into two parts [45] :

$G_i\left(\psi ,r,\epsilon ,T_f,M_f,E_f,F_f\right)=E_i\left(\psi ,r,\epsilon \right)+X_i\left(\psi ,r,\epsilon ,T_f,M_f,E_f,F_f\right)$ (3)

where, the AG interaction potential energy $E_i\left(\psi ,r,\epsilon \right)$ is responsible to make AG stabilization, the AG generalized action free energy $X_i\left(\psi ,r,\epsilon ,T_f,M_f,E_f,F_f\right)$ is responsible to make AG suit variations in environments, under influences of temperature $T_f$ , magnetic $M_f$ , electric $E_f$ and force $F_f$ fields. Under the influence of $T_f$ field, the $G_i^{Au}\left(T\right)$ and $G_i^{Cu}\left(T\right)$ sequences have been calculated by experimental enthalpies of few alloys (Supplementary Information (1) 1.3, 1.4).

2.6. AG W_i-Sequences

In the AGC valence bond theory, the AG potential energy curve function is of great significance: 1) It has developed traditional atom pair interaction potential functions, such as Mie potential and Lennard-Jones potential, into the MAI (Many Atoms Interaction)-AG $W\left(\psi ,r,\epsilon \right)$ function [31] , that makes AG theory have solid physical foundation and can establish the $G_i$ function worthy of the name. 2) It has developed elasticity: The force constant $K\text{'}(\psi ,r,\epsilon ,E,\theta )$ , bulk modulus $B\left(\psi ,r,\epsilon ,E,\theta \right)$ , Young’s modulus $Y\left(\psi ,r,\epsilon ,E,\theta \right)$ and shear modulus $\mu \left(\psi ,r,\epsilon ,E,\theta \right)$ equations associated with the $\psi ,r,\epsilon ,E$ and Debye temperature $\theta$ have been derived, that make $G_i$ function can respond to $F_f$ field. 3) It has developed Debye specific heat theory: The ${\theta }_i$ sequences may be calculated approximately by $\theta \left(E,V\right)$ function, that make $G_i$ function can respond to $T_f$ field, through Debye vibration free energy $X_i^D$ and generalized attaching action free energy $X_i^E$ . 4) It has developed Grüneisen’s thermal expansion theory: The Grüneisen’s two constants $K\left(\psi ,r,\epsilon ,E,\theta \right)$ and $Q\left(\psi ,r,\epsilon ,E,\theta \right)$ have been derived, which are useful to calculate AGC thermal expansion ${\alpha }_i$ sequences, that make $G_i$ function can respond to $F_f$ and $T_f$ fields. 5) It has developed magnetism and electricity: The AGC total potential energy $E_i\left(\psi ,r,\epsilon \right)$ includes contributions of bond energies of covalent, magnetic and free electrons, that makes $G_i$ function can respond to $T_f$ , $M_f$ , $E_f$ and $F_f$ fields. These achievements have been tested by studies of some pure metals [32] [33] [34] [35] (Supplementary Information (1) 1.3, 1.4).

2.7. AG-GET Function

According to the SSP-systematic opening circuit evolution principle, i.e., “The sources of system are systematic genes, which determine structures of systems; structures determine properties; properties suit environments; environments optimize genes; genes optimize systems; systems optimize environments.” we have established the AG-Gibbs energy transmissive AG-GET $\Omega \left(x,T,\sigma \right)$ function of alloy phases in $T_f$ field, which describes the systematic correlativity of the characteristic Gibbs energy $G^{*}\left(x,T,\sigma \right)$ (chemical potential of alloy phase), AG- probabilities ( $x_i^{Au}\left(x,T,\sigma \right)$ , $x_i^{Cu}\left(x,T,\sigma \right)$ ) occupied at the AG-GE ( $G_i^{Au}$ , $G_i^{Cu}$ ) levels and degeneracy g-factor of AG-probabilities occupied at the lattice points [37] [40] [45] [46] [47] [48] . It is the AG-GET function in the real sense of the term, which may be also called as the systematic opening evolution function (Supplementary Information (2) 4.1,4.2):

$\Omega \left(x,T,\sigma \right)=g\left(x_i^{Au}\left(x,T,\sigma \right),x_i^{Cu}\left(x,T,\sigma \right)\right)\times \exp \left[-G^{*}\left(x,T,\sigma \right)/\left(k_{B}T\right)\right]$ (4)

where, $k_B$ is Boltzmann’s constant, and $G^{*}\left(x,T,\sigma \right)$ may be obtained by the AG-GE additive law:

$G^{*}\left(x,T,\sigma \right)={\displaystyle \underset{i=0}{\overset{I=12}{\sum }}\left(x_i^{Au}\left(x,T,\sigma \right)\times G_i^{Au}\left(T\right)+x_i^{Cu}\left(x,T,\sigma \right)\times G_i^{Cu}\left(T\right)\right)}$ (5)

The first Gibbs energy function worthy of the name for alloy phases in the history has been derived:

$G\left(x,T,\sigma \right)=-k_{B}T\times \ln \left(\Omega \left(x,T,\sigma \right)\right)=G^{*}\left(x,T,\sigma \right)-T{S}^c\left(x,T,\sigma \right)$ (6)

The AG arranging configurational entropy $S^c\left(x,T,\sigma \right)$ -function is obtained from the degeneracy g-factor:

$S^c\left(x,T,\sigma \right)=k_B\ln \left(g\left(x_i^{Au}\left(x,T,\sigma \right),x_i^{Cu}\left(x,T,\sigma \right)\right)\right)$ (7)

2.8. AG-GET Modes

The essential definition of equilibrium evolution described by AG-minimal Gibbs energy path (MGEP) equilibrium transmissive mode is that “the AG-GE ( $G_i^{Au}$ , $G_i^{Cu}$ ) levels and AG-probabilities ( $x_i^{Au}$ , $x_i^{Cu}$ ) occupied at the lattice points/AG-GE levels can respond immediately and change synchronously with each small variation in temperature, and proceed along the minimal Gibbs energy $G_{\min }\left(x\right)-T$ path, supposing to have no obstacle to atom movement”. The $\Delta G_{\min }^m\left(x\right)-T$ , $x_i^{Au}\left(x\right)-T$ , $x_i^{Cu}\left(x\right)-T$ equilibrium path charts may be calculated by iso-order degree $\Delta G_{\sigma }^m\left(x\right)-T$ or isothermal $\Delta G_T^m\left(x\right)-\sigma$ path methods (Figure 2, Supplementary Information (1) 2.1). The MGEP mode leads to equilibrium EAP theory, which includes EAP thermodynamics [45] [46] [47] [48] , EAP valence bond theory [41] [42] [43] [44] and EAHNP phase diagrams of alloy phases (Figure 3). Then, the equilibrium EAS theory and EAHNP phase diagrams of Au-Cu system (Figure 4) have been established by Equilibrium Lever numerical NLN method of isothermal Gibbs energy $\Delta G_T^m-x$ curves.

The essential definition of subequilibrium evolution described by AG-tracking path (AG-TP) subequilibrium transmissive mode is that “the AG-Gibbs energy ( $G_i^{Au}$ , $G_i^{Cu}$ ) levels can change immediately with each small variation in temperature, but AG-probabilities ( $x_i^{Au}$ , $x_i^{Cu}$ ) occupied at the lattice point/AG-GE levels cannot change synchronously, due to exist potential well obstacle to atom

Figure 2. Path charts on disordering AuCu alloy. (a) Three dimension mixed Gibbs energy $\Delta G^m-T-\sigma$ chart with equilibrium $\Delta G_e^m-T$ path and subequilibrium $\Delta G_s^m-T$ path; (b) Two dimension equilibrium $\Delta G_e^m-T$ path and subequilibrium $\Delta G_s^m-T$ path.

Figure 3. EAHNP diagrams of AuCu type phase. (a) Two dimension equilibrium isothermal $\Delta G_T^m-x$ phase diagram with phase boundary curves $\Delta G_{PB}^m\left(x\right)$ ; (b) Two dimension equilibrium iso-order degree $T_{\sigma }-x$ phase diagrams with phase boundary curves $T_{PB}\left(x\right)$ and experimental jumping $T_j$ -temperatures [15] . The symbols “O”, “D”, “C” denote respectively ordered region, disordered region, stoechmetric compound.

Figure 4. Phase diagrams of Au-Cu system. (a) Equilibrium order degree EAHNP phase diagram with iso-order degree $T_{\sigma }-x$ curves and experimental jumping $T_j$ -temperatures [15] ; (b) Subequilibrium order degree SAHNP phase diagram with iso-order degree $T_{\sigma }-x$ curves and experimental jumping $T_j$ -temperatures [46] .

movement, and proceed along the higher $G_s\left(x\right)-T$ path than $G_{\min }\left(x\right)-T$ path.” The $G_s\left(x\right)-T$ , $x_{i,s}^{Au}\left(x\right)-T$ , $x_{i,s}^{Cu}\left(x\right)-T$ subequilibrium path charts may be calculated by the AG-Tracking Path AG-TP mode of the experimental mixed enthalpy $\Delta H^m-T$ path (Figure 2, Supplementary Information (1) 3). This AG-TP mode leads to subequilibrium SAP theory, which includes SAP thermodynamics, SAP valence bond theory and SAHNP phase diagrams of alloy phases. Then, the SAS theory and SAHNP phase diagrams of Au-Cu system (Figure 4) have been established by the Subequilibrium Crosspoint Numerical SCN method of isothermal Gibbs energy $\Delta G_T^m-x$ curves.

2.9. AHNP Bank of Alloy Systems

The AHNP information bank of alloy systems consists of AG holographic warehouse (Supplementary Information (1) 1), equilibrium EAHNP and subequilibrium SAHNP warehouses of alloys, phases and systems (Supplementary Information (1) 2, 3), in which all of path charts and phase diagrams are calculated by steps $\Delta x=0.\text{5at}\%$ , $\Delta T=1\text{K}$ , and $\Delta \sigma =0.0001$ . As clicking one network point, the holographic information of the clicked alloy may be obtained immediately, that includes the AG valence electron arranging structures and holographic properties of its alloy genome, AG arranging structures and holographic properties of its alloy phases and alloy phase arranging structures and holographic properties of its alloy organization.

The AHNP bank is operable platform for researchers to achieve social road from the discovery to deployment of advanced alloys: AHNP intellectual discovery → AG arranging holographic design → AG tracking holographic certification → AG tracking high effectiveness manufacturing → AG tracking safe deployment, then allows us to significantly decrease drain on time, manpower, material, financial and energy resources for getting advanced alloys and provides information support for the AGSP.

3. Discussions

The SSP is a universality knowledge system to study relationships among components-genes-structures-properties-environments of various systems. It consists of principles of the globality, SG universality, opening circuit evolution, unity of opposites and so on, and laws of WOP, AG-TP, combination of probability and determinacy, and so on.

1) Discussion on SG universality principle

In the single level system, the correctness of AG universality can be directly explained. The genes of the free atoms system are the quantum state electron sequence described by four quantum numbers. The diversity of electronic structures of free atoms is attributed to combination and arrangement of quantum state electrons, and diversity of their properties is attributed to response of diverse structures to environments. The genes of musical compositions system are the musical note sequence. The diversity of structures of compositions is attributed to combination and arrangement of musical notes, and diversity of their properties is attributed to response of their structures to performers.

It is very difficult to seek genes of biologic system with multi-levels. Since Mendel’s pioneering investigation in genetic characters of peas to formally began the Human Genome Project (HGP), it took 125 years. It is also very difficult to seek alloy genes of alloy system with three levels. Since 1978, we began to seek alloy genes, and up to now we propose AGSP, which has taken 40 years. During the same period, we have discovered that the human social genes (HSG) are the characteristic people sequences. A single HSG is a characteristic people existing in the center of the clan cluster and occupying the professional position in the professional cluster, as well as carrying biological and social genetic information about structures and properties from one generation to the next.

2) Discussion on systematic globality principle

In the SSP, each level has its own structure unit sequences and corresponding properties sequences. Properties are that the structure unit sequences make response to environment act on the system. Relationship of the structures and properties of each level with environments is described by corresponding level theory. Relationship between two level theories is described by transmissive equation. Therefore, the SMM science based on the QMAC-electron theory of free atoms may be described by theories-, structure unit sequences- and equations- chains besides information-chain (Figure 1):

The theories chain: QMAC-electron theory → AG theory → EAP and SAP theories → EAS and SAS theories.

The structure units sequences chain: quantum state electron sequences → basic atom state sequences with quantum state valence electron sequences → AG sequences → EAP and SAP sequences.

The equations chain: Schrödinger equation of free atoms → AG separable functions, MAI-AG potential energy curve function, AG Gibbs energy function → AG-GET function, equilibrium MGEP and subequilibrium AG-TP modes → equilibrium ELN and subequilibrium SCN methods.

These achievements allow us to recognize that the QMAC-electron theory of alloys is a questionable alloy theory lacking of theories-, structure unit sequences- and equations-chains and that the quantum state electron sequences and Schrödinger equation in the QMAC-electron theory of free atom system are unreasonably used to establish electron theory of alloys, violating the matched rule of structure unit sequence and structure level theory.

These achievement allow us to recognize that the QMAC-cluster variation thermodynamics and CALPHAD-thermodynamics of alloy phases are the questionable theories, in which the “so-called Gibbs energy functions” established unreasonably by “constituent atoms”, “constituent atom pairs” and “constituent atom clusters”, carrying no valence electron structures, physical properties and Gibbs energies worthy of the name, contain a number of unknown parameters to be adjusted until the Gibbs energy function is capable of representing the chosen information well. They don’t know that the constituent atom states have split in the alloys. These advances allow us also to recognize that Prigogine’s dissipative structures thermodynamics is a phenomenological theory without genes. It can not reveal essence of chaotic phases (Supplementary Information (1) 3.4).

These advances allow us also to recognize that the “Materials Genome Initiative (MGI)” [18] [49] [50] is a questionable plan without materials gene (MG) sequences, MG-theory and MG-techniques. It is unthinkable that the MGI has been launched as a strategic plan in the United States America.

The SMM science established by guide of SSP-principles and laws will prove stimulating to researchers working on biologic, non-biologic and social systems, and who may well find values in using them as philosophic thinking modes for establishing applied system sciences. The AGSP will become a paragon to create a new era for establishing various kinds of applied system sciences and systematic gene engineerings.

Supplementary Information (1): Methods Chain

Creative SSP laws/rules evolve AHNP methods/techniques! The SSP laws/ rules mainly include the WOP law, combining rule of imaginal logic thinking and mathematical logic thinking, matched rule of structure unit sequence and structure level, determined rule for SG-sequence, united rule of probability and determinacy, integrated rule of analysis and synthesis, combining rule of theory and practice, and so on. These SSP laws/rules have been demonstrated to be very useful for establishing the SMM science and AHNP bank. Among the SSP laws/rules the WOP law has significant scientific, useful and economic values.

The WOP law means that a system can be reproduced from its few parts in the biosystems. Such as a tree can be reproduced from a seed, a branch and even a leaf of this tree; and an animal can be reproduced from its a egg cell and even a body cell. The secrets of these processes are deeply hided in the structures and properties of biological gene sequences, as well as their transmissive mode, that is attributed to natural selection. The WOP law also means that the holographic information of a system can be produced from its few parts. Due to discover AG sequences (Figure 5), the AHNP bank of the fcc-based lattice Au-Cu system, which includes the AG holographic warehouse, equilibrium EAHNP and subequilibrium SAHNP warehouses of alloys, phases and Au-Cu system, has been established based only experimental enthalpies $H\left(x,T\right)$ and volumes $V\left(x,T\right)$ of a few alloys. The secrets of this process are also hided in the structures and properties of AG sequences, as well as their transmissive principles and modes, that have been revealed by the AHNP methods/techniques road of the WOP law in the following sections.

1) AG-WOP methods

The establishment of AG-holographic warehouse is the first stage for establishing EAHNP bank of alloy systems. It can be explained by taking the fcc-based lattice Au-Cu system, which has been considered to be classical paragon studied nearly by almost conventional and modern measurement methods for one century [51] - [77] . The AG-WOP methods include the separated and chosen methods of AG potential energies and volumes sequences, and the calculated methods of AG valence electron structures, physical and thermodynamic properties [78] - [86] (Figure 6).

a) WOP methods of AG ε_i and v_i sequences

Determining AG potential energy ${\epsilon }_i$ and volume $v_i$ sequences is the first step to establish AG-holographic warehouse by the WOP law. In the $q\left(x,T,\sigma \right)$ functions of alloys, each function contains only four unknown parameters of primary states ( $q_0^A$ , $q_I^B$ ) and terminal states ( $q_I^A$ , $q_0^B$ ) of AG-sequences ( $q=\epsilon$ or $v$ ,). We have separated out the ${\epsilon }_i^{Au}$ , ${\epsilon }_i^{Cu}$ , $v_i^{Au}$ and $v_i^{Cu}$ sequences from only experimental enthalpies and volumes of pure Au and Cu metals, as well as AuCu and AuCu₃ compounds, and chosen the 5^th $\Delta {\epsilon }^m\left(x,T,\sigma \right)$ and 6^th $\Delta v^m\left(x,T,\sigma \right)$ functions suitable for fcc-based lattice Au-Cu system by excess $\Delta {\epsilon }^{ex}\left(x,T,\sigma \right)$ and excess $\Delta v^{ex}\left(x,T,\sigma \right)$ of AuCu- and AuCu₃-type sublattice alloys relative to disordered Au_(1−x)Cu_x alloys [45] (Figure 7):

Figure 5. Three structural unit sequences. (a) Basic coordination cluster B-sequences and alloy gene A-sequences in fcc-based lattice Au-Cu system: $B_0^{Au}\cdots B_8^{Au}\cdots B_{12}^{Au}$ , $B_{12}^{Cu}\cdots B_4^{Cu}\cdots B_0^{Cu}$ and $A_0^{Au}\cdots A_8^{Au}\cdots A_{12}^{Au}$ , $A_{12}^{Cu}\cdots A_4^{Cu}\cdots A_0^{Cu}$ ; (b) AG-characteristic crystal C-sequences: $C_0^{Au}\cdots C_8^{Au}\cdots C_{12}^{Au}$ , $C_{12}^{Cu}\cdots C_4^{Cu}\cdots C_0^{Cu}$ .

Figure 6. Theoretic road of AG-WOP methods. (a) AG-theory based on the experimental techniques or the first-principles electron theory of alloys; (b) (c) AG potential energies and volumes-sequences calculated by separated theory of potential energies and volumes of alloy genes; (d) (e) AG potential energy curves and physical properties calculated by valence bond theory of characteristic crystals; (f) (g) AG Gibbs energy curves and thermodynamic properties calculated by thermodynamics of characteristic crystals.

Figure 7. Choice among E- and V-functions. (a) (b) Excess potential energy differences $\Delta E_{AuCu}^{ex}\left(x\right)$ and $\Delta E_{AuC{u}_3}^{ex}\left(x\right)$ relative to disordered Au_(1−x)Cu_x alloys; (c) (d) Excess volume differences $\Delta V_{AuCu}^{ex}\left(x\right)$ and $\Delta V_{AuC{u}_3}^{ex}\left(x\right)$ relative to disordered Au_(1−x)Cu_x alloys.

$\begin{array}{c}\Delta {\epsilon }^m\left(x,T,\sigma \right)={\displaystyle \underset{i=0}{\overset{I}{\sum }}\left[{\left(\frac{i}{I}\right)}^2{x}_i^{Au}\right]}\left({\epsilon }_I^{Au}\left(T\right)-{\epsilon }_0^{Au}\left(T\right)\right)\\ +{\displaystyle \underset{i=0}{\overset{I}{\sum }}\left[{\left(\frac{I-i}{I}\right)}^2{x}_i^{Cu}\right]\left({\epsilon }_0^{Cu}\left(T\right)-{\epsilon }_I^{Cu}\left(T\right)\right)}\end{array}$ (8)

$\begin{array}{c}\Delta v^m\left(x,T,\sigma \right)={\displaystyle \underset{i=0}{\overset{I}{\sum }}\left[{\left(\frac{i}{I}\right)}^2{x}_i^{Au}\right]}\left(v_I^{Au}\left(T\right)-v_0^{Au}\left(T\right)\right)\\ +{\displaystyle \underset{i=0}{\overset{I}{\sum }}\left[\frac{2I\left(I-i\right)-{\left(I-i\right)}^2}{I^2}{x}_i^{Cu}\right]\left(v_0^{Cu}\left(T\right)-v_I^{Cu}\left(T\right)\right)}\end{array}$ (9)

where, the $x_i^{Au}=x_i^{Au}\left(x,T,\sigma \right)$ , $x_i^{Cu}=x_i^{Cu}\left(x,T,\sigma \right)$ , $I=12$ .

We have also calculated the potential energies and volumes of Au, Cu, Au₃Cu, AuCu and AuCu₃ by three methods with local density approximation in the QMAC electron theory of alloys for obtaining $\Delta {\epsilon }^m\left(x,T,\sigma \right)$ and $\Delta v^m\left(x,T,\sigma \right)$ functions. It has been discovered that the accuracies of both the potential energies and volumes at the same instant cannot allow correctly to determine $\Delta \epsilon \left(x,T,\sigma \right)$ and $\Delta v\left(x,T,\sigma \right)$ function [84] [85] . Perhaps there is “a uncertainty principle on relation between potential energy and volume” in the QMAC electron theory of alloys, that is likely to be similar to the “Heisenberg uncertainty principle”.

b) WOP methods of AG ψ_i sequences

Determining AG-VE structure ${\psi }_i$ sequences is the second step to establish AG-holographic warehouse based on AG ${\epsilon }_i$ and $v_i$ sequences. In the AGC valence bond theory, the AGC/AG-VE structure is described by one-atom valence electron y-hybrized by some basic atom states in the basic atom state ${\phi }_k$ sequences: $\psi ={\displaystyle \underset{k}{\sum }{C}_k{\phi }_k}$ . Each ${\phi }_k$ state is described by free electron ( $s_f$ ), covalent electron ( $s_c$ , $p_c$ , $d_c$ , $f_c$ ), magnetic electron ( $d_m$ , $f_m$ ) and non-valence electron ( $p_n$ , $d_n$ , $f_n$ ) (Table 1, Table 2). Then, the bond length $r\left(\psi \right)$ and lattice constant $a\left(\psi ,r\right)$ equations, as well as bond energy $\epsilon \left(\psi ,r\right)$

and potential energy $E\left(\psi ,r,\epsilon \right)$ equations have been established. The presive solutions of valence electron structures of AG ${\psi }_i^{Au}$ and ${\psi }_i^{Cu}$ sequences have been calculated, through the crosspoint method of the equilattice constant $a\left(\psi ,r\right)$ -line and equipotential energy $E\left(\psi ,r,\epsilon \right)$ -line (Figure 8). The primary

Table 1. The basic atomic states (the number of free electrons (s_f), covalent electrons (s_c, d_c), magnetic electrons (d_m) and non-valent electrons (d_n)) and the characteristic properties of some fcc and bcc pseudo-crystals of Fe metal [33] .

Table 2. The basic atomic states (the number of free electrons (s_f, p_f), covalent electrons (s_c, d_c) and non-valent electrons (d_n)) and the characteristic properties (single bond radii (R), lattice constant (a) cohesive energy (E_c)) of some fcc pseudo-crystals of Au and Cu metals [86] .

Figure 8. AG-valence electron structures. (a) Hybridized compositions of three basic states and electronic structures for primary $C_0^{Au}$ characteristic crystals at 0 K; (b) Hybridized compositions of three basic states and electronic structures for primary $C_{12}^{Cu}$ characteristic crystals at 0 K.

( ${\psi }_0^{Au}$ , ${\psi }_{12}^{Cu}$ ) and terminal ( ${\psi }_{12}^{Au}$ , ${\psi }_0^{Cu}$ ) valence electron structure in the ${\psi }_i^{Au}$ and ${\psi }_i^{Cu}$ sequences (Table 3) are as follows [45] [86] :

${\psi }_0^{Au}=\left[Xe\right]5d_n{}^{3.4013}5d_c{}^{5.8506}6s_f{}^{1.7481}$ ; ${\psi }_{12}^{Au}=\left[Xe\right]5d_n{}^{3.4218}5d_c{}^{6.0936}6s_f{}^{1.4847}$ .

${\psi }_{12}^{Cu}=\left[Ar\right]3d_n{}^{5.0778}3d_c{}^{4.7254}4s_f{}^{1.1968}$ ; ${\psi }_0^{Cu}=\left[Ar\right]3d_n{}^{4.2877}3d_c{}^{5.04954}4s_f{}^{1.6628}$ .

where, $\left[Xe\right]$ ―Xenon electronic structure, $\left[Ar\right]$ ―Argon electronic structure.

According to potential surface forms and filling rules of $s_c$ , $p_c$ , $d_c$ and $f_c$ shells, the AG arranging crystallography of alloy phases may be established for revealing essence of crystalline structure categories of alloy phases.

According to energies, lengths and directions of $s_c$ , $p_c$ , $d_c$ and $f_c$ covalent bonds, the AG arranging covalent bond theory of alloy phases may be established for revealing essence of mechanical properties of alloy phases [41] [42] [43] [44] .

c) WOP methods of AG physical properties

Determining AG physical properties sequences is the third step to establish AG-holographic warehouse based on AG ${\epsilon }_i$ , $v_i$ and ${\psi }_i$ sequences. In the AGC-VB theory, the AG potential curve $W\left(\psi ,r,\epsilon ,E\right)$ equation has been established, besides bond length, lattice constant and potential energy equations. According to $W\left(\psi ,r,\epsilon ,E\right)$ equation, the potential energy curves (Figure 6), cohesive energies ( $E_c$ ), Debye temperatures ( $\theta$ ), bulk moduli ( $B$ ) and so on of AGC crystals have been calculated (Table 3). The $W\left(\psi ,r,\epsilon ,E\right)$ function has developed traditional atom pair interaction potential theory [87] [88] , Debye specific heat theory, Grüneisen’s thermal expansion theory, elasticity, magnetism and electricity, that makes $G_i$ function can respond to $T_f$ , $M_f$ , $E_f$ and $F_f$ fields [31] [32] [33] [34] [35] .

According to coupling energies and transmissive mode of magnetic and free electrons, the magnetic and electric conducting theories of alloy phases will be respectively established for revealing essence of magnetic and electric properties of alloy phases.

d) WOP methods of AG G_i-sequences

Determining AG Gibbs energy $G_i$ sequences is the fourth step to establish AG-holographic warehouse. Under the influence of temperature field, the AG $G\left(\psi ,r,\epsilon ,E,\theta ,T_f,M_f,E_f,F_f\right)$ may be simply rewritten as the $G_i\left(T\right)$ , which may be described by two parts: a temperature independent potential energy $E\left(0\right)$ and a temperature dependent generalized vibration/action free energy $X_i^v\left(T\right)$ :

$G_i\left(T\right)=E_i\left(0\right)+X_i^v\left(T\right)$ . (1)

The $X_i^v\left(T\right)$ includes Debye vibration free energy $X_i^D\left(T\right)$ and generalized attaching action free energy $X_i^E\left(T\right)$ :

$X_i^v\left(T\right)=X_i^D\left(T\right)+X_i^E\left(T\right)$ . (2)

The relationship of the $X_i^D\left(T\right)$ with Debye vibration energy $U_i^D\left(T\right)$ , Debye vibration entropy energy $T{S}_i^D\left(T\right)$ and Debye vibration specific heat $c_{pi}^D\left(T\right)$ is,

Table 3. The valence electron structures (the number of free electrons (s_f), covalent electrons (s_c, d_c) and non-valent electrons (d_n)), volumes(v), potential energies (ε) and single bond radii (R) of alloy genes, cohesive energy ( $E_c$ ), Debye temperatures ( $\theta$ ) and bulk moduli (B) of AG characteristic crystals in the fcc based Au-Cu system [45] .

$X_i^D\left(T\right)=U_i^D\left(T\right)-T{S}_i^D\left(T\right)={\displaystyle \int c_{pi}^D\left(T\right)\text{d}T}-T{\displaystyle \int \frac{c_{pi}^D\left(T\right)}{T}\text{d}T}$ (3)

The $c_{pi}^D\left(T\right)$ may be calculated by Debye theory.

The relationship of the $X_i^E\left(T\right)$ with attaching action energy $U_i^E\left(T\right)$ , attaching action entropy energy $T{S}_i^E\left(T\right)$ and attaching action specific heat $c_{pi}^E\left(T\right)$ is as following,

$X_i^E\left(T\right)=U_i^E\left(T\right)-T{S}_i^E\left(T\right)={\displaystyle \int c_{pi}^E\left(T\right)\text{d}T}-T{\displaystyle \int \frac{c_{pi}^E\left(T\right)}{T}\text{d}T}$ (4)

The $c_{pi}^E\left(T\right)$ may be obtained only by experimental enthalpies $H\left(x,T\right)$ of few alloys (Supplementary Information (2) 3), which includes contributions of electron excitation, expansion work ( $P\Delta V$ ) of volume, energy of formation of holes, variation of potential energy $\Delta E\left(T\right)$ with temperature, and so on. For fcc-based Au-Cu system, the AG $G\left(T\right)$ function with multi-level energetic functions in the real sense of the term have been established and shown in the Figure 9.

2) WOP methods of EAHNP warehouse

The establishment of equilibrium EAHNP path charts of single alloys and EAHNP phase diagrams of single phases and Au-Cu system is the second stage for establishing AHNP bank of the fcc-based lattice Au-Cu system only based on AG holographic warehouse, that needs no other experimental data.

a) WOP methods of EAHNP path of alloys

According to the $G_i^{Au}\left(T\right)$ and $G_i^{Cu}\left(T\right)$ sequences, mixed Gibbs energy $\Delta G^m\left(x,T,\sigma \right)$ function and essential definition of equilibrium order ƒ disorder transition, the equilibrium $\Delta G_e^m-T$ , ${\sigma }_e-T$ , $x_{i,e}^{Au}-T$ and $x_{i,e}^{Cu}-T$ paths on equilibrium order ƒ disorder transition as function of temperature for the stoichiometric Au₃Cu alloy have been calculated by the minimal Gibbs energy path MGEP method, which may be performed by the iso-order degree Gibbs energy $\Delta G_{\sigma }^m-T$ path method or isothermal Gibbs energy $\Delta G_T^m-\sigma$ path method (Figure 10), using calculated steps $\Delta T=1\text{K}$ and $\Delta \sigma =0.0001$ (Figure 11). Then, EAHNP $q_e-T$ path charts of the Au₃Cu alloy are calculated by $x_{i,e}^{Au}-T$ and $x_{i,e}^{Cu}-T$ paths. It has following general behaviors: The order → disorder transition upon heating and disorder → order transition upon cooling proceed along the same minimal mixed Gibbs energy $\Delta G_{\min }^m-T$ path. Namely, there exists no so-called hysteresis phenomenon between both transitions; The $\Delta G_{\min }^m-T$ path is continuous and homogenous, and has stability and linearity without jumping phenomenon. The $T_c$ critical temperature is determined by the crosspoint method of $\Delta G_{\min }^m-T$ and $\Delta G_{dis.}^m\left(x,\sigma =0\right)-T$ curves.

b) EAHNP phase diagrams of phases

By the $\Delta G_{\min }^m-T$ path MGEP method, the $\Delta G_e^m(x,T,\sigma )$ , ${\sigma }_e^{}(x,T)$ , $x_{i,e}^{Au}(x,T,\sigma )$ and $x_{i,e}^{Cu}(x,T,\sigma )$ path charts on equilibrium order ƒ disorder transition as function of composition, temperature and order degree for all alloys in the Au₃Cu-type sublattice phase may be calculated, using calculated steps $\Delta x=0.\text{5at}.\%$ , $\Delta T=1\text{K}$ , $\Delta \sigma =0.0001$ . Then, the three-dimension $q_e-x-T$ and two-dimension $q_{x,e}-T$ , $T_{q,e}-x$ and $q_{T,e}-x$ EAHNP phase diagrams may be constructed [48] .

By analogous method, the EAHNP phase diagrams of the AuCu- and AuCu₃-type sublattice ordered phases [45] [47] , as well as fcc-based lattice disordered phase may be established.

Figure 9. AG-thermodynamic properties. (a) (b) Generalized vibration heat capacity curves in Au-Cu system: $C_{p.0}^{Au.v}\cdots C_{pi}^{v.Au}\cdots C_{p12}^{v.Au}$ , $C_{p.0}^{Cu.v}\cdots C_{p.i}^{Cu.v}\cdots C_{p.12}^{Cu.v}$ ; (c) (d) Generalized vibration free energy curves: $X_0^{Au.v}\cdots X_i^{Au.v}\cdots X_{12}^{Au.v}$ , $X_0^{Cu.v}\cdots X_i^{Cu.v}\cdots X_{12}^{Cu.v}$ ; (e) (f) Enthalpy of formation curves: $H_0^{Au}\cdots H_i^{Au}\cdots H_{12}^{Au}$ , $H_0^{Cu}\cdots H_i^{Cu}\cdots H_{12}^{Cu}$ ; (g) (h) Gibbs energy curves: $G_0^{Au}\cdots G_i^{Au}\cdots G_{12}^{Au}$ , $G_{12}^{Cu}\cdots G_i^{Cu}\cdots G_0^{Cu}$ .

Figure 10. Minimal Gibbs energy path method. (a) The iso-order degree Gibbs energy $\Delta G_{\sigma }^m-T$ path method; (b) The isothermal Gibbs energy $\Delta G_T^m-\sigma$ path method; (c) $\Delta G_e^m\left(\Delta G_{\min }^m\right)-T$ path on equilibrium order ƒ disorder transition of stoichiometric Au₃Cu alloy; (d) $\sigma -T$ path on equilibrium order ƒ disorder transition of stoichiometric Au₃Cu alloy; (e) (f) Three-dimension $x_{i,e}^{Au}-T-i$ and $x_{i,e}^{Cu}-T-i$ path charts on disordering or ordering the stoichiometric Au₃Cu alloy.

Figure 11. AG-tracking path method. (a) Minimal mixed Gibbs energy $\Delta G_{\min }^m-T$ path on disordering AuCuI alloy obtained by the iso-order degree Gibbs energy equilibrium path method; (b) Experimental mixed enthalpy $\Delta H_s^m-T$ path method based on the network chart of the iso-order degree mixed enthalpies equilibrium paths, attaching equilibrium $\Delta H_e^m-T$ path; (c) $\Delta G_{\min }^m$ (OABC) and $\Delta G_s^m$ (Oabc) paths of the equilibrium and subequilibrium order → disorder transitions, where $\Delta G_h^m$ denotes superheated driving Gibbs energy; (d) ${\sigma }_e$ (OABCD) and ${\sigma }_s$ (OabcD) paths of the equilibrium and subequilibrium order → disorder transitions, where $\Delta {\sigma }_h$ denotes hysteresis order degree; (e) Three-dimension $x_i^{Au}-T-i$ SAHNP chart; (f) Three-dimension $x_i^{Cu}-T-i$ SHNP chart.

c) EAHNP phase diagrams of alloy system

The solid fcc-based lattice Au-Cu system is a competitive system of the Au₃Cu-, AuCu- and AuCu₃-type sublattice ordered phases, as well as fcc-based lattice disordered phase. The equilibrium three dimension $\Delta G^m-T-x$ phase diagram may be established by the equilibrium lever numerical method of isothermal Gibbs energy curves [23] . Then, the EAHNP phase diagrams of fcc-based lattice Au-Cu system have been established (Figure 4).

3) WOP methods of SAHNP warehouses

The subequilibrium SAHNP warehouse established by the AG-tracking path AG-TP method is the third stage for establishing AHNP bank of the fcc-based lattice Au-Cu system only based on AG holographic warehouse, EAHNP warehouse and experimental mixed enthalpy $\Delta H^m\left(x,T\right)$ or volume $\Delta V^m\left(x,T\right)$ paths of alloys! It includes SAHNP path charts of single alloys [46] , and SAHNP phase diagrams of single phases and Au-Cu system [23] .

a) WOP methods of SAHNP path of alloys

The steps of the AG-TP method for subequilibrium SAHNP path charts of alloys may be presented by taking ordered AuCuI(AuCu) alloy as an example, which is a classic paragon (Figure 11):

i) To establish equilibrium $\Delta G_e^m\left(\sigma \right)-T$ path chart by the iso-order degree Gibbs energy method at first, then, the $\Delta H_e^m-T$ path and iso-order degree enthalpies $\Delta H_{\sigma }^m-T$ path network chart may be drawn up (Figure 11(a), Figure 11(b)).

ii) To draw up the experimental mixed enthalpy $\Delta H_s^m-T$ path on the network chart of the $\Delta H_{\sigma }^m-T$ path, then, the subequilibrium, ${\sigma }_s-T$ , ${\displaystyle \sum x_{i,s}^{Au}}-T-i$ , ${\displaystyle \sum x_{i,s}^{Cu}}-T-i$ , $\Delta G_s^m-T$ (Figures 11(c)-(f)) and SAHNP path charts of the AuCu alloy may be calculated.

The common behaviors of the SAHNP charts are as follows: The order degree keeps unchanged in a broad temperature range from 0 K to 593 K by slow heating rate; There exists complex variation in order degree, during the temperature range from onset temperature ( $T_{onset}$ ) to jumping $T_j$ -temperature, that is only about 100 K interval; Upon increasing heating rate, the $T_{onset}$ -temperature increases, and the $T_j$ -temperature decreases, i.e., so called Retro-Effect [68] [70] ; The $\Delta G_s^m-T$ path is higher than the equilibrium $\Delta G_e^m-T$ path. Their difference ( $\Delta G_e^m\left(T\right)-\Delta G_s^m\left(T\right)$ ) is the superheated driving Gibbs energy of the subequilibrium alloy; There exists a thermal hysteresis Gibbs energy loop between the order → disorder and disorder → order transitions.

b) SAHNP phase diagrams of phases

By the AG-tracking path (AG-TP) method on experimental mixed enthalpy path, the $\Delta H_s^m\left(x,T,\sigma \right)$ , ${\sigma }_s\left(x,T\right)$ , $x_{i,s}^{Au}\left(x,T,\sigma \right)$ and $x_{i,s}^{Cu}\left(x,T,\sigma \right)$ path charts on subequilibrium order → disorder transition as function of composition, temperature and order degree for all alloys in the AuCu-type sublattice phase may be calculated, using calculated steps $\Delta x=0.\text{5at}.\%$ , $\Delta T=1\text{K}$ and $\Delta \sigma =0.0001$ . Then, the three-dimension $q_s-x-T$ and two-dimension $q_{x,s}-T$ , $T_{q,s}-x$ and $q_{T,s}-x$ SAHNP phase diagrams may be constructed [47] [48] .

By analogous method, the SAHNP phase diagrams of the Au₃Cu- and AuCu₃- type ordered phases may be established.

c) SAHNP phase diagrams of alloy systems

The subequilibrium three dimension $\Delta G^m-x-T$ phase boundries of solid fcc-based lattice Au-Cu system may be established by the crosspoint numerical method of isothermal Gibbs energy curves [23] . Then, the SAHNP phase diagrams of fcc-based lattice have been established (Figure 4).

d) Values of AG-TP method

The AG-TP method has special values. It can make us describe all real situations happened in evolutional progress of alloys, through tracking one property only. Then, we can reveal essence of chaotic phases, observe some new phenomena, discover some new mechanisms, explain some difficult questions left over by history and establish SAHNP path charts, that may be explained by the disordering path of AuCuI( $A_8^{Au}{A}_4^{Cu}$ ) compound [46] .

i) When $0\text{K}\le T<T_{onset}\left(593\text{K}\right)$ , it is the subequilibrium unchanged period of the AuCuI( $A_8^{Au}{A}_4^{Cu}$ ), of which the behaviors as fellows (Figure 11): There exists no variation in order degree until $T_{onset}=593\text{K}$ , $x_8^{Au}=x_4^{Cu}=\text{5}0\text{at}.\%$ , other AG concentrations $x_i^{Au}=x_i^{Cu}=0$ , $A_8^{Au}$ and $A_4^{Cu}$ still occupied at the $G_8^{Au}$ and $G_4^{Cu}$ energy levels, respectively; The superheated driving Gibbs energy ( $\Delta G_e^m\left(T\right)-\Delta G_s^m\left(T\right)$ ) between the equilibrium (OABC) and experiment (Oabc) paths increases with rising temperature, of which the difference (-373J/mol at 593 K) is very small (Figure 11(c)); The generalized action energies ( $U_8^{Au.v}$ and $U_4^{Cu.v}$ ) increase with rising temperature, but their values are much smaller than their potential well depths. At 593K the $\left|E_8^{Au}/U_8^{Au.v}\right|=27.78$ , $\left|E_4^{Cu}/U_4^{Cu.v}\right|=27.11$ (Table 4, Table 5); The new knowledge is that only depending on superheated driving Gibbs energy and generalized action energy, the $A_8^{Au}$ and $A_4^{Cu}$ “stem alloy genes” cannot change their lattice positions, if without a special mechanism.

ii) When $T_{onset}<T\le 620\text{K}$ and $1.000>{\sigma }_s\ge 0.\text{99}0$ (Figure 12(a), Figure 12(b)), it is the unequilibrium growing period of the high order degree AuCu(H) alloy. According to the peculiar orders in AG concentrations (Table 6, Table 7): $x_5^{Au}<x_4^{Au}>x_6^{Au}$ , $x_6^{Cu}<x_8^{Cu}>x_7^{Cu}$ , we have discovered the 3(RA-SA) mechanism, which is that three AG-pairs in a cell-scale region occur Resonance Activating (RA) to accumulate vibration energies and Synchro Alternating (SA) to change their lattice positions and energetic states by two manners (Figure 13(a), Figure 13(b)). The new knowledge is that there exists cell-scale heterogeneity in the AuCu(H) alloy, which includes early-3(RA-SA) cell region with more 3(RA-SA) cells, middle-3(RA-SA) cell region with less (RA-SA) cells and late-AuCuI ( $A_8^{Au}{A}_4^{Cu}$ ) cell region with no (RA-SA) cell, that lay a structural foundation for forming statistic periodic anti-direction (SPAD) region. The $A_8^{Au}$ and $A_4^{Cu}$ are named as the “signal alloy genes of 3 (RA-SA) mechanism”.

iii) When $\text{62}0\text{K}<T\le \text{65}0\text{K}$ and $0.990>{\sigma }_s\ge 0.925$ (Figure 12(c), Figure 12(d)), it is the unequilibrium growing period of the SPAD-AuCuII region. We

Table 4. The $E_8^{Au}$ , $U_8^{Au.v}$ , $c_p{}_8^{Au}$ , $\left(E_8^{Au}+U_8^{Au.v}\right)$ , $\left(T\times S_8^{Au.v}\right)$ , $G_8^{Au}$ , $\left|E_8^{Au}/U_8^{Au.v}\right|$ of the $A_8^{Au}$ -alloy gene [46] .

1) The generalized vibration energy includes contributions of Debye vibration energy, existing energy of electrons, forming energy of holes, expansion work of volume and variation in potential energy with temperature; 2) The melting of pure Au metal; 3) The melting of pure Cu metal.

Table 5. The $E_4^{Cu}$ , $U_4^{Cu.v}$ , $c_p{}_4^{Cu}$ , $\left(E_4^{Cu}+U_4^{Cu.v}\right)$ , $\left(T\times S_4^{Cu.v}\right)$ , $G_4^{Cu}$ , $\left|E_4^{Cu}/U_4^{Cu.v}\right|$ of the $A_4^{Cu}$ -alloy gene [46] .

1)The generalized vibration energy includes contributions of Debye vibration energy, existing energy of electrons, forming energy of holes, expansion work of volume and variation in potential energy with temperature; 2) The melting of pure Au metal; 3) The melting of pure Cu metal.

Figure 12. AG tracking structure path charts. (a) (b) As $0.990\le {\sigma }_s\le 1$ the $x_i^{Au}$ and $x_i^{Cu}$ concentration distributions on disordering AuCuI( $A_8^{Au}{A}_4^{Cu}$ ); (c) (d) As $0.925\le {\sigma }_s\le 0.990$ , the $x_i^{Au}$ and $x_i^{Cu}$ concentration distributions; (e) (f) As $0.807\le {\sigma }_s\le 0.925$ , the $x_i^{Au}$ and $x_i^{Cu}$ concentration distributions; (g) (h) As ${\sigma }_j=0.807\to {\sigma }_L=0.445$ , the $x_i^{Au}$ and $x_i^{Cu}$ concentration distributions.

Table 6. Systematic correlativity of the equilibrium transition temperature $T_e$ , subequilibrium transition temperature $T_s$ , order degree $\sigma$ and $A_i^{Au}$ -alloy gene concentrations ( $x_i^{Au}$ ) on disordering AuCuI( $A_8^{Au}{A}_4^{Cu}$ ) compound.

The $\sigma$ denotes ${\sigma }_e$ corresponding to $T_e$ or ${\sigma }_s$ corresponding to $T_s$ .

Table 7. Systematic correlativity of the equilibrium transition temperature $T_e$ , subequilibrium transition temperature $T_s$ , order degree $\sigma$ and $A_i^{Cu}$ -alloy gene concentrations ( $x_i^{Cu}$ ) on disordering AuCuI( $A_8^{Au}{A}_4^{Cu}$ ) compound.

The $\sigma$ denotes ${\sigma }_e$ corresponding to $T_e$ or ${\sigma }_s$ corresponding to $T_s$ .

have discovered that as $T_s=650\text{K}$ and ${\sigma }_s=0.925$ , the concentrations $x_9^{Au}$ and $x_3^{Cu}$ display the “emergent phenomenon”, which is defined as that their concentrations are larger than ones of the corresponding alloy genes in the disorder state: $x_9^{Au}>x_{9,dis}^{Au}$ , $x_3^{Cu}>x_{3,dis}^{Cu}$ . During this period, there exists “region- scale heterogeneity”: the original early-3(RA-SA) cell regions with more (RA-SA) cells grow into the early-(SPAD-AuCuII) regions, the original middle-(RA-SA) cell regions with less 3(RA-SA) cells grow into the middle-(RA-SA) cell regions with more 3(RA-SA) cells and the original late-AuCuI( ${\epsilon }_0^B$ ) regions with no (RA-SA) cell grow into the late-(RA-SA) cell regions with less (RA-SA) cells (Figure 13(c), Figure 13(d)). The new knowledge is that the AuCu(H) region and SPAP-AuCuII region are co-existed, two distinct x-ray diffraction patterns may be observed at the $v_0^B$ . This alloy with region-scale heterogeneity is called as the “Pseudo-Two-Phases (PTP)” AuCu alloy with the AuCu(H) regions and the SPAD-AuCuII regions. The use of the phrase “PTP” is intended to convey the impression that two heterogeneous subequilibrium regions have the same composition, different order degrees, and belong in the same ordered AuCu-type sublattice phase.

iv) When $650\text{K}<T<683\text{K}$ and $0.925>\sigma >0.807$ (Figure 12(e), Figure 12(f)), it is the unequilibrium mature period of the SPAD-AuCuII alloy. During

Figure 13. The 3(RA-SA) mechanism. (a) AG arranging structure after A-manner alternating in 3(RA-SA) mechanism of three $A_8^{Au}$ - $A_4^{Cu}$ pairs in a single cell on disordering AuCuI( $A_8^{Au}{A}_4^{Cu}$ ); (b) AG arranging structure after B-manner alternating in 3(RA-SA) mechanism of three $A_4^{Cu}$ - $A_8^{Au}$ pairs in a single cell; (c) AuCu(H) alloy containing three regions with cell-scale heterogeneity; (d) PTP-[AuCuI(H) + PAP-AuCuII] alloy containing three regions with region-scale heterogeneity; (e) AG arranging structure of a early SPAP-AuCuII cell.

this period, the middle-(RA-SA) cell regions and late-(RA-SA) cell regions are respectively transferred into the middle- and late-SPAD-AuCuII regions. The new knowledge is “the SPAD-AuCuII alloy containing early-, middle- and late-SPAD-AuCuII regions is a statistic periodic anti-direction structure stacking of incommensurate tetragonal cells containing more $A_8^{Au}$ and $A_4^{Cu}$ stem alloy genes and anti-direction boundary with (RA-SA) cells containing less $A_8^{Au}$ and $A_4^{Cu}$ stem alloy genes along b axis.” It has no strict long periodic cell, and may be called as Prigogine’s dissipative structures or chaotic phase [89] [90] .

v) When $683\text{K}\ge T_s$ and $0.807\ge {\sigma }_s\ge 0.786$ (Figure 12(g), Figure 12(h)), it is the unequilibrium jumping period of the SPAD-AuCuII alloy. We have discovered that as ${\sigma }_j=0.807$ , $A_9^{Au}$ and $A_3^{Cu}$ display the “maximum emergent phenomenon”: $x_9^{Au}=x_3^{Cu}=\text{6}.\text{515at}\%$ ; as $\sigma =0.786$ , $A_7^{Au}$ and $A_5^{Cu}$ display the “maximum emergent phenomenon”, $x_7^{Au}=x_5^{Cu}=13.458\text{at}\%$ . Therefore, the $A_9^{Au}$ , $A_3^{Cu}$ , ${\sigma }_{j,1}=0.807$ and $A_7^{Au}$ , $A_5^{Cu}$ , ${\sigma }_{j,2}=0.786$ are called as the first and second “jumping alloy genes”, “jumping order degrees”, respectively. As the ${\sigma }_s$ -order degree less than the ${\sigma }_j$ -order degree, the disordering begins to translate from a single splitting of the $A_8^{Au}$ and $A_4^{Cu}$ stem alloy genes to a “universal splitting” or “snowslide splitting” of the stem and jumping alloy genes, which is an essential condition to happen a discontinuous jump change from the high order degree SPAD-AuCuII to the Low order degree AuCu(L) with low order degree ( $\sigma <0.4545$ ). After the $T_j$ -temperature, the AuCu(L) alloy is continuously unequilibrium transformed into the equilibrium Disordered AuCu (D) alloy. The new knowledge is that the heterogeneous successive transitions on disordering stoichiometric AuCuI( $A_8^{Au}{A}_4^{Cu}$ ) by slow heating rate are as follows: AuCuI ( $A_8^{Au}{A}_4^{Cu}$ ) → AuCu (H) → PTP-AuCu → SPAD-AuCuII → AuCu (L) → AuCu (D).

The AG-TP method may lead to the AG tracking holographic certification, AG tracking high effectiveness manufacturing and AG tracking safe deployment, as well as AG subequilibrium and unequilibrium kinemics.

The fcc-based lattice Au-Cu system has been studied as a classical paragon with complete order-disorder transition for about one century. However, a equilibrium phase diagram of Au-Cu system in the real sense of the term has not been established, and numerous questions about order-disorder transition have been raised [73] [74] , because one of many reasons is that researchers have got used to recognizing experimental phenomena observed during very slow variation in temperature to be thermodynamic equilibrium phenomena: The experimental jumping $T_j$ -temperatures are erroneously considered as the $T_c$ -critical temperatures of equilibrium order-disorder transition of alloys, although the experimental jumping ${\sigma }_j$ -order degrees are about 0.8 ~ 0.6 and the experimental short range order degrees exist at the temperatures considerable above the $T_j$ -temperature [53] [65] ; The experimental composition-dependent $T_j-x$ curve is erroneously considered as the phase boundary $T_c-x$ curve between ordered and disordered phases. The researchers in the CALPHAD- and QMAC-communities took these miss understandings of experimental phenomena as the selected information, then adjusted parameters in so-called Gibbs energy functions and established so-called equilibrium phase diagrams to achieve the best representation of the selected information [18] . These phase diagrams have no scientific significance as well as useful value [15] [91] [92] [93] [94] (Figure 14).

Figure 14. Phase diagrams of Au-Cu system calculated by CALPHAD- and QMAC-thermodynamics. (a) Experimental jumping T_j temperatures, erroneously considered as equilibrium critical T_c temperatures and experimental limit composition ranges of long range ordered (LRO) Au₃Cu-, AuCu- and AuCu₃-type at room temperature [60] [61] [91] [92] ; (b) Calculated Au-Cu phase diagrams by CALPHAD-thermodynamics [92] ; (c) Calculated Au-Cu phase diagram by CALPHAD-thermodynamics [15] ; (d) Calculated Au-Cu phase diagram by CALPHAD-thermodynamics [15] ; (e) Calculated Au-Cu phase diagram by CALPHAD- thermodynamics [93] ; (f) Calculated Au-Cu phase diagram by QMAC-thermodynamics [94] .

Supplementary Information (2): Equations Chains

The logical thinking in images is a foundation of the logical thinking in mathematics; while it is deepened by the logical thinking in mathematics. A clear and perfect theory can be created by comprehensive application of the two ways of thinking. The SMM science is the framework consisting of the correlated three-level theories. Each level theory is quantitatively described by main equations. We have encountered some difficult mathematical problems during the mathematization process. It has taken a few years or even more than ten years of effort to solve each of the problems, such as the AG separable functions, MAI-AGC potential energy curve function and the AG Gibbs energy transmissive AG-GET function. We have deeply realized that the “mathematization is a heavy storm fortification of establishing SMM science”.

1) Equations in AG separable theory

a) AG additive law

$\{\begin{array}{l}q\left(x,T,\sigma \right)={\displaystyle \underset{i=0}{\overset{I}{\sum }}{x}_i^A\left(x,T,\sigma \right)q_i^A\left(T\right)}+{\displaystyle \underset{i=0}{\overset{I}{\sum }}{x}_i^B\left(x,T,\sigma \right)q_i^B\left(T\right)}\\ q_A\left(x,T,\sigma \right)={\displaystyle \underset{i=0}{\overset{I}{\sum }}{x}_i^A\left(x,T,\sigma \right)q_i^A\left(T\right)/x_A}\\ q_B\left(x,T,\sigma \right)={\displaystyle \underset{i=0}{\overset{I}{\sum }}{x}_i^B\left(x,T,\sigma \right)q_i^B\left(T\right)/x_B}\end{array}$ (1-1)

Where, the q denotes AG potential energy (E), volume (v) and other extensive properties, which are the functions of temperature (T) only; the $x_i^A\left(x,T,\sigma \right)$ and $x_i^B\left(x,T,\sigma \right)$ are AG concentrations; which are the functions of the composition (x), temperature (T) and order degree ( $\sigma$ ); $x_A$ and $x_B$ denote compositions of the A and B components, respectively.

b) AG separable functions

The potential energies and volumes are the most basic properties of AG sequences. In order to make the Equation (1-1) become the simple, applicable and separable $q\left(x,T,\sigma \right)$ function performed by WOP law, three type relations of $q_i^A$ and $q_i^B$ with i have been designed [37] :

Type I of linear relation,

$\{\begin{array}{l}{q}_i^A=q_0^A+\left(i/I\right)\left(q_I^A-q_0^A\right)\\ q_i^B=q_I^B+\left[\left(I-i\right)/I\right]\left(q_0^B-q_I^B\right)\end{array}$ (1-2)

Type II of concave parabola relation,

$\{\begin{array}{l}{q}_i^A=q_0^A+{\left(q/I\right)}^2\left(q_I^A-q_0^A\right)\\ q_i^B=q_I^B+{\left[\left(I-i\right)/I\right]}^2\left(q_0^B-q_I^B\right)\end{array}$ (1-3)

Type III of convex parabola relation,

$\{\begin{array}{l}{q}_i^A=q_0^A+\left(2i/I\right)\left(q_I^A-q_0^A\right)-{\left(i/I\right)}^2\left(q_I^A-q_0^A\right)\\ q_i^B=q_I^B+\left[2\left(I-i\right)/I\right]\left(q_0^B-q_I^B\right)-{\left[\left(I-i\right)/I\right]}^2\left(q_0^B-q_I^B\right)\end{array}$ (1-4)

where, the q denotes AG potential energy or volume; I denotes coordinative number; the $q_0^A$ and $q_I^B$ denote, respectively, potential energy ${\epsilon }_0^A$ and ${\epsilon }_I^B$ or volume $v_0^A$ and $v_I^B$ of the primary state alloy genes; the $q_I^A$ and $q_0^B$ denote, respectively, ${\epsilon }_I^A$ and ${\epsilon }_0^B$ or $v_I^A$ and $v_0^B$ of the terminal state alloy genes.

By combining Equations (1-2)-(1-4) and substituting them into Equation (1-1), nine $q\left(x,T,\sigma \right)$ functions can be obtained. In Table 8, the nine general $q\left(x,T,\sigma \right)$ functions of alloy phases can be used to compounds, ordered and disordered alloy phases. In Table 9, the nine $q\left(x,T\right)$ -functions can be used for the disordered alloy phases only.

c) Choices of E- and V-functions

The 5^th E-function in Table 8 has been chosen for describing Au-Cu system, according to the excess potential energies $\Delta E^{ex}\left(x,0,{\sigma }_{\max }\right)$ of AuCu- and AuCu₃-type sublattice ordered Au_(1-x)Cu_x alloys with maximum order degrees calculated by nine E-functions at room temperature [45] (Figure 7(a), Figure 7(b)).

The 6^th V-function in Table 8 has been chosen for describing Au-Cu system, according to the calculated values of the mixed volumes $\Delta v^m\left(x,T_r,{\sigma }_{\max }\right)$ and

Table 8. Nine general $q\left(x,T,\sigma \right)$ -functions of alloy phases at T K.^*

*: Here, the $x_i^A=x_i^A\left(x,T,\sigma \right)$ and $x_i^B=x_i^B\left(x,T,\sigma \right)$ ; $q=\epsilon$ or v; i is the number of B-atoms in the nearest coordinative configuration $\left[\left(I-i\right)\text{A}+i\text{B}\right]$ ; (I-i) is the number of A-atoms in the nearest coordinative configuration $\left[\left(I-i\right)\text{A}+i\text{B}\right]$ ; I is the coordinative number.

Table 9. $q\left(x,T\right)$ -functions of disordered alloy phases ( $\sigma =0$ ) at T K.^*

*: Here, the $x_A$ and $x_B$ denote compositions of the A and B components; $q=\epsilon$ or v; I is the coordinative number.

excess volumes $\Delta v^{ex}\left(x,T_r,{\sigma }_{\max }\right)$ of the AuCu- and AuCu₃-type sublattice ordered Au_(1−x)Cu_x alloys with maximum order degrees calculated by the nine V-functions at room temperature (Figure 7(c), Figure 7(d)).

In the Au-Cu system, the correlative function between AG potential energies,

$\{\begin{array}{l}{\epsilon }_i^{Au}={\epsilon }_0^{Au}+{\left(i/I\right)}^2\left({\epsilon }_I^{Au}-{\epsilon }_0^{Au}\right)\\ {\epsilon }_i^{Cu}={\epsilon }_I^{Cu}+{\left[\left(I-i\right)/I\right]}^2\left({\epsilon }_0^{Cu}-{\epsilon }_I^{Cu}\right).\end{array}$ (1-5)

The potential energy function of alloys is,

$\begin{array}{c}\epsilon \left(x,T,\sigma \right)=x_{Au}{\epsilon }_0^{Au}\left(T\right)+x_{Cu}{\epsilon }_I^{Cu}\left(T\right)+{\displaystyle \underset{i=0}{\overset{I}{\sum }}\left[{\left(\frac{i}{I}\right)}^2{x}_i^{Au}\right]}\left({\epsilon }_I^{Au}\left(T\right)-{\epsilon }_0^{Au}\left(T\right)\right)\\ +{\displaystyle \underset{i=0}{\overset{I}{\sum }}\left[{\left(\frac{I-i}{I}\right)}^2{x}_i^{Cu}\right]\left({\epsilon }_0^{Cu}\left(T\right)-{\epsilon }_I^{Cu}\left(T\right)\right)}\end{array}$ (1-6)

The correlative function between AG volumes:

$\{\begin{array}{l}{v}_i^{Au}=v_0^{Au}+{\left(i/I\right)}^2\left(v_I^{Au}-v_0^{Au}\right)\\ v_i^{Cu}=v_I^{Cu}+\left[2\left(I-i\right)/I\right]\left(v_0^{Cu}-v_I^{Cu}\right)-{\left[\left(I-i\right)/I\right]}^2\left(v_0^{Cu}-v_I^{Cu}\right).\end{array}$ (1-7)

The volume function of alloys is,

$\begin{array}{c}v\left(x,T,\sigma \right)=x_{Au}{v}_0^{Au}\left(T\right)+x_{Cu}{v}_I^{Cu}\left(T\right)+{\displaystyle \underset{i=0}{\overset{I}{\sum }}\left[{\left(\frac{i}{I}\right)}^2{x}_i^{Au}\right]}\left(v_I^{Au}\left(T\right)-v_0^{Au}\left(T\right)\right)\\ +{\displaystyle \underset{i=0}{\overset{I}{\sum }}\left[\frac{2I\left(I-i\right)-{\left(I-i\right)}^2}{I^2}{x}_i^B\right]\left(v_0^{Cu}\left(T\right)-v_I^{Cu}\left(T\right)\right)}\end{array}$ (1-8)

where, ${\epsilon }_0^{Au}=-3.8140\text{eV}$ , ${\epsilon }_I^{Au}=-4.0555\text{eV}$ , ${\epsilon }_0^{Cu}=-3.6488\text{eV}$ , ${\epsilon }_I^{Cu}=-3.4824\text{eV}$ ; $v_0^{Au}=16.9581\left({10}^{-3}{\text{nm}}^3\right)$ , $v_I^{Au}=15.5213\left({10}^{-3}{\text{nm}}^3\right)$ , $v_0^{Cu}=12.6828\left({10}^{-3}{\text{nm}}^3\right)$ , $v_I^{Cu}=11.8124\left({10}^{-3}{\text{nm}}^3\right)$ , $I=12$ . The results calculated by Equations ((1-6) and (1-8)) are in good agreement with experimental values. Their essential features are as follows:

i) The ${\epsilon }_i^{Au}$ decreases with increasing number i of Cu-atoms in the coordinative configuration: $\left({\epsilon }_i^{Au}-{\epsilon }_0^{Au}\right)={\left(i/I\right)}^2\left({\epsilon }_I^{Au}-{\epsilon }_0^{Au}\right)$ .

ii) The ${\epsilon }_i^{Cu}$ decreases with increasing number (I-i) of Au-atoms in the coordinative configuration: $\left({\epsilon }_i^{Cu}-{\epsilon }_I^{Cu}\right)={\left[\left(I-i\right)/I\right]}^2\left({\epsilon }_0^{Cu}-{\epsilon }_I^{Cu}\right)$ .

iii) The $v_i^{Au}$ decreases with increasing number i of Cu-atoms in the coordinative configuration: $\left(v_i^{Au}-v_0^{Au}\right)={\left(i/I\right)}^2\left(v_I^{Au}-v_0^{Au}\right)$ .

iv) The $v_i^{Cu}$ increase with increasing number (I-i) of Au-atoms in the coordinative configuration: $\left(v_i^{Cu}-v_I^{Cu}\right)=\left[2\left(I-i\right)/I\right]\left(v_0^{Cu}-v_I^{Cu}\right)-{\left[\left(I-i\right)/I\right]}^2\left(v_0^{Cu}-v_I^{Cu}\right)$ .

2) Equations of AGC valence bond theory

AGC valence bond theory of AGC theory includes the description of AG valence electron structures $\psi$ and three main functions: MAI-AGC lattice constant equation $a\left(\psi ,r\right)$ with many atoms interaction (MAI), MAI-AGC cohesive energy equation $E_c\left(\psi ,r,\epsilon \right)$ and MAI-AGC potential energy curve function $W\left(\psi ,r,\epsilon \right)$ , that have developed Pauling’s, Yu’s and Grigorovvich’s metallic bond theories [25] [26] [27] .

a) Description of AG ψ -sequences

An AG $\psi$ is represented by AGC $\psi$ of corresponding AGC crystal and formed by some basic atom states in the basic atom state ${\phi }_k$ sequences [32] [33] [34] [35] (Table 1, Table 2):

$\psi ={\displaystyle \underset{k=1}{\overset{n}{\sum }}{c}_k{\phi }_k}$ (2-1)

In our opinion the metallic bond is a mixed bond containing unsaturated covalent electron bond, free electron bond and magnetic electron bond, except non-valence electrons. The unsaturated covalent electron bond network is formed by covalent electrons ( $n_c$ ) of an atom with its neighbors. The number and symmetry of the overlapping hybrid orbitals directly determine the crystalline structure of the AGC crystal. These unsaturated covalent bonds are responsible for the most fraction of cohesive energy of the crystals. The free electron bond consists of near free electrons ( $n_f$ ), which uniformly fill the free space of the unsaturated covalent bond network and are responsible for the conductivity and plasticity. The magnetic electron bond consists of magnetic electrons ( $n_m$ ) responsible for the magnetism. Both the free electron bond and the magnetic electron bond have small contributions to cohesive energy of the crystal. The non-valence electrons ( $n_n$ ) can be treated as the “chemically inert” ion core.

The occupations of electrons in each basic atomic state obey Pauli Exclusion principles. If $s_k^c$ , $p_k^c$ and $d_k^c$ denote respectively the number of the covalent electrons of s, p and d shells in the k basic atomic state; $s_k^f$ , $d_k^m$ and $d_k^n$ denote respectively the numbers of free s electrons, magnetic d electrons and non-valence d electrons, then the state parameters of one-atom state of the AGC crystal can be obtained by the following equations:

$\{\begin{array}{l}{s}_c={\displaystyle \underset{k}{\sum }{c}_k{s}_k^c},p_c={\displaystyle \underset{k}{\sum }{c}_k{p}_k^c},d_c={\displaystyle \underset{k}{\sum }{c}_k{d}_k^c}\\ d_m={\displaystyle \underset{k}{\sum }{c}_k{d}_k^m},s_f=n_f={\displaystyle \underset{k}{\sum }{c}_k{s}_k^f},d_n={\displaystyle \underset{k}{\sum }{c}_k{d}_k^n}\\ n_c=s_c+p_c+d_c,n_v=n_c+n_f\\ R\left(1\right)={\displaystyle \underset{k}{\sum }{c}_k{R}_k\left(1\right)},{\displaystyle \underset{k}{\sum }{c}_k=1}\text{.}\end{array}$ (2-2)

where $R_k\left(1\right)$ is the single bond radius which can be obtained by Pauling’s equation slightly modified by the authors.

For Fe-AGC crystal,

$\{\begin{array}{l}{R}_k\left(1\right)=\left(0.1481-0.04{\delta }_k\right)\text{nm}\\ {\delta }_k=d_k^c/n_v\end{array}$ (2-3)

For Au-AGC crystal,

$\{\begin{array}{l}{R}_k\left(1\right)=\left(0.1520-0.0497{\delta }_k\right)\text{nm}\\ {\delta }_k=d_k^c/n_v\end{array}$ (2-4)

For Cu-AGC crystal,

$\{\begin{array}{l}{R}_k\left(1\right)=\left(0.1350-0.025{\delta }_k\right)\text{nm}\\ {\delta }_k=d_k^c/n_v\end{array}$ (2-5)

From the state parameters of the basic state atoms, the characteristic properties (lattice constant, cohesive energy, and atomic magnetic moment) of the relative pseudo-crystal consisting of atoms in each kind of basic atomic state can be calculated by a series of equations presented in the following sections (Table 1, Table 2).

b) The MAI lattice constant equation

An AG volume may be described by lattice constant equation $a\left(\psi ,r\right)$ with many atoms interaction (MAI) of corresponding AGC crystal, which is the function of AG valence electron structure $\psi$ and bond length $r_s$ . Pauling’s bond length equation should be suitable to the various bond lengths in the AGC crystal [29] (in nm):

$\{\begin{array}{l}{r}_1=2R\left(1\right)-\beta \lg \left(n_1\right)\\ r_2=2R\left(1\right)-\beta \lg \left(n_2\right)\\ ⋮\\ r_s=2R\left(1\right)-\beta \lg \left(n_s\right)\end{array}$ (2-6)

where $r_s$ denotes the length of the s neighboring bond, $R\left(1\right)$ denotes the single bond radius, $n_s$ denotes the number of pairs of covalent electrons on the relative bond, $\beta$ is taken as 0.060 determined by Pauling. From Equation (2-6), the following equations may be derived as

$\{\begin{array}{l}{n}_c={\displaystyle \underset{s}{\sum }{I}_s{n}_s}\\ n_1=\frac{n_c}{{\displaystyle \underset{s}{\sum }{I}_s}\times {10}^{\left(r_1-r_s\right)/\beta }}\\ n_s=n_1\times {10}^{\left(r_1-r_s\right)/\beta }\end{array}$ (2-7)

where $I_s$ represents the number of the same kind of the s bond. For a AGC crystal, the lattice constant can be expressed by its bond length

$a\left(\psi ,r\right)=r_s/G_s$ (2-8)

For a crystal with fcc structure, $G_1=1/\sqrt{2}$ , $G_2=1$ , and $G_3=\sqrt{6}/2$ .

For a crystal with bcc structure, $G_1=\sqrt{3}/2$ , $G_2=1$ , and $G_3=\sqrt{2}$ .

Combining Equations ((2-7) and (2-8)), we have (in nm)

$a\left(\psi ,r\right)=\left(1/G_1\right)\left\{2R\left(1\right)-\beta \lg \left[n_c/{\displaystyle \underset{s}{\sum }{I}_s}\times {10}^{\left(G_1-G_s\right)a/\beta }\right]\right\}$ (2-9)

It has been proved that the $a\left(\psi ,r\right)$ is the convergent transcendental equation; and that there is a certain lattice constant for a pseudo-crystal with known basic atom $\psi$ state. It may be obtained by repeated calculation methods [29] (Table 1, Table 2).

c) AGC cohesive energy equation

Considering contributions of bond energies of covalent electrons, near free electrons, and magnetic electrons to cohesive energy of AGC crystal, the MAI-AGC cohesive energy equation $E_c\left(\psi ,r,\epsilon \right)$ has been established [31] [32] [33] [34] [35] ,

$E_c\left(\psi ,r,\epsilon \right)=-E\left(\psi ,r,\epsilon \right)=N\left({\displaystyle \underset{s}{\sum }{\epsilon }_{c,s}}+{\epsilon }_f+{\epsilon }_m\right)$ (2-10)

where, N is Avogardro’s number.

In the AGC crystal, the bond energy of a kind of bonds should be directly proportional to the number $n_{c,s}$ of covalent electron pairs on each bond and the number $I_s$ of the same kind of bonds and inversely proportional to its bond length, i.e.,

${\epsilon }_{c,s}\propto \left(I_s{n}_{c,s}/r_{s,0}\right)\times f$ (2-11)

In Pauling’s and Yu’s opinion, the bonding capacities of s, p, d and f orbitals are $f_s=1$ , $f_p=\sqrt{3}$ , $f_d=\sqrt{5}$ and $f_f=\sqrt{7}$ , respectively. If the number of covalent electrons of s, p, d and f orbitals in the hybrid orbital are respectively $s_c$ , $p_c$ , $d_c$ and $f_c$ , the bonding capacity of the hybrid orbital should be

$f=\sqrt{\alpha }+\sqrt{3\beta }+\sqrt{5\gamma }+\sqrt{7\delta }$ (2-12)

where,

$\{\begin{array}{l}\alpha =s_c/n_T^v,\beta =p_c/n_T^v,\gamma =d_c/n_T^v,\delta =f_c/n_T^v\\ n_T^v=s_c+p_c+d_c+f_c+n_f+n_m\end{array}$ (2-13)

In the AGC crystal, there should be free electrons besides the covalent electrons. The free electrons are regarded as distributed through the whole crystal. If the average energy ${\epsilon }_f$ of the free electron bond is assumed directly proportional to the number $n_f$ of free electrons and inversely proportional to the effective average band length, i.e.

${\epsilon }_f\propto \left(n_f/{\bar{r}}_0\right)\times f^{\prime }$ (2-14)

where ${\bar{r}}_0$ is the average effective bond lengths of the AGC crystal in equilibrium,

${\stackrel{\rightharpoonup }{r}}_0=\frac{{\displaystyle \underset{s}{\sum }{I}_s{r}_{s,0}}}{{\displaystyle \underset{s}{\sum }{I}_s}}$ (2-15)

In Yü’s opinion, the bonding capacity of free electrons should be

$f^{\prime }=\sqrt{2{\alpha }^{\prime }}$ (2-16)

where,

${\alpha }^{\prime }=n_f/n_T^v$ (2-17)

For ferromagnetic AGC crystals, it also may be assumed that the average magnetic coupling energy is directly proportional to the number n_m of magnetic electrons,

${\epsilon }_m\propto B\left(n_m/r_{m,0}\right)$ or ${\epsilon }_m\propto B^{\prime }{n}_m$ (2-18)

where the $B^{\prime }$ denotes the bonding capacity of the magentic electrons.

Therefore, the MAI cohesive energy equation of AGC crystals may be written as (in KJ/mol) [31]

$E_c\left(\psi ,r_0,\epsilon \right)=A\left[{\displaystyle \sum \frac{I_s{n}_s}{r_{s,0}}}f+\frac{n_f}{{\bar{r}}_0}{f}^{\prime }+B^{\prime }{n}_m\right]$ (2-19)

Where, the constant may be determined by experimental cohesive energies of primary AGC crystals, i.e., pure elements. The $B^{\prime }$ value may be determined by the experimental specific heat peak due to magnetic transformation of primary AGC crystal.

d) AGC potential energy curve function

Under the condition of homogeneous deformation, various bond lengths may be expressed by the following equations:

$\{\begin{array}{l}{r}_{s,0}=R_s{r}_{1,0},r_s=R_s{r}_1\\ r_{m,0}=R_m{r}_{1,0},r_m=R_m{r}_1\\ {\stackrel{\rightharpoonup }{r}}_0={\displaystyle \underset{s}{\sum }{I}_s{R}_s{r}_{1,0}}/{\displaystyle \underset{s}{\sum }{I}_s},\bar{r}={\displaystyle \underset{s}{\sum }{I}_s{R}_s{r}_1}/{\displaystyle \underset{s}{\sum }{I}_s}\end{array}$ (2-20)

And their reduced bond lengths are identically equal:

$\frac{r_{s,0}}{r_s}=\frac{r_{m,0}}{r_m}=\frac{{\bar{r}}_0}{\bar{r}}=\frac{r_0}{r}$ (2-21)

where, the $r_0$ and r denote the shortest bond lengths of the AGC crystal in equilibrium and unequilibrium states, respectively.

Therefore, the MAI-AGC potential function with many atoms interaction (MAI) for one molar crystal may be simplified as [31]

$W\left(\psi ,r,\epsilon \right)=E_c\left[-n{\left(\frac{r_0}{r}\right)}^x+\left(n-1\right){\left(\frac{r_0}{r}\right)}^{nx/\left(n-1\right)}\right]$ (2-22)

In the equilibrium state, r = r₀, the absolute value of molar potential energy is approximately equal to the molar cohesive energy, i.e., $W\left(\psi ,r_0,\epsilon \right)=-E_c\left(\psi ,r_0,\epsilon \right)$ at 0K. The MAI-AGC potential energy curve function $W\left(\psi ,r,\epsilon \right)$ makes AG theory have solid physical foundation.

e) Equations of elasticities

i) The force constant

Suppose that the MAI-AGC crystal is subjected to a hydrostatic pressure, which produces a small deformation u. The potential function can be written as [31]

$W\left(\psi ,\left(r_0+u\right),\epsilon \right)=W\left(\psi ,r_0,\epsilon \right)+{\left(\frac{\text{d}W}{\text{d}r}\right)}_{r=r_0}u+\frac{1}{2}{\left(\frac{{\text{d}}^{2}W}{\text{d}{r}^2}\right)}_{r=r_0}{u}^2+\cdots$ (2-23)

According to the principle of mechanics and Debye model, the force constant k and exponents x in MAI-AGC potential energy curve function can be obtained

$k=4{\theta }^2{k}_B^2{r}_0^{2}m/{\hslash }^2{j}^2$ (2-24)

$x^2=4{\theta }^2{k}_B^2{r}_0^{2}m\left(n-1\right)/n{\hslash }^2{j}^2{E}_c$ (2-25)

where $\theta$ is Debye temperature; $k_B$ , Boltzman’s constant; $\hslash$ , Planck’s constant; m, atomic mass; j and n, which are respectively the multiple of the half cutoff wave length and another exponent (n) of the MAI-AGC potential function, which can be obtained from two experimental values of the experimental linear thermal expansion coefficients at two temperatures of primary AGC crystals [33] .

ii) Bulk modulus

Suppose that the crystal is subjected to a hydrostatic pressure, which produces the Hookeian volume change $\left(V-V_0\right)/V_0=-\epsilon$ . Then B, defined by [31]

$B=\left(\text{d}W/\text{d}V\right)/\left[\left(V-V_0\right)/V_0\right]$ (2-26)

is the bulk modulus of elasticity.

According to $W\left(\psi ,r,\epsilon \right)$ and $x^2$ functions,

$B=\frac{4{\theta }^2{k}_B^2{r}_0^{2}m}{9{\hslash }^2{V}_0}\cdot \frac{1}{j^2}\left(1-2\epsilon \right)$ (2-27)

iii) Young’s modulus, shear modulus and Poisson’s ratio

Suppose that the polycrystal undergoes a very small homogeneous extension along an arbitrary direction such as $r_0$ . Then, according to the definition of Young’s modulus, we have

$Y=\left(\text{d}W/\text{d}V\right)/\left[3\left(r-r_0\right)/r_0\right]$ (2-28)

The factor 1/3 is introduced due to dW/dV being the total differential and Young’s modulus being one-dimensional quantity.

Suppose the extension strain to be ${\epsilon }^{\prime }$ . Then, we have

${\epsilon }^{\prime }=\left(r-r_0\right)/r_0,{\epsilon }^{\prime }\ge 0$ (2-29)

According to $W\left(\psi ,r,\epsilon \right)$ and $x^2$ functions,

$Y=\frac{4{\theta }^2{k}_B^2{r}_0^{2}m}{9{\hslash }^2{V}_0}\cdot \frac{1}{j^2}\left(1-{\epsilon }^{\prime }/2\right)$ (2-30)

It is well known that there are the relationships among various moduli as follows

$\{\begin{array}{l}B=Y/3\left(1-2\sigma \right)\\ Y=2\mu \left(1+\sigma \right)\\ \mu =Y/\left[2\left(1+\sigma \right)\right]\\ \sigma =1/2-1/6\left(1-{\epsilon }^{\prime }/2+2\epsilon \right)\end{array}$ (2-31)

where, $\mu$ is the shear modulus, $\sigma$ , the Poisson’s ratio.

f) Equations of specific heats

The specific heat of magnetic AGC crystals includes two ports: the normal specific heat $C_{vn}$ caused by lattice vibration and magnetic specific heat $C_{vm}$ caused by magnetic transformation [33] [34] [35] . They can be calculated from the Debye theory and Brillouin function respectively

$\{\begin{array}{l}{C}_{vn}=3R\left[12{\left(\frac{T}{{\theta }_i}\right)}^3{\displaystyle {\int }_0^{{\theta }_i/T}\frac{y^3\text{d}y}{{\text{e}}^y-1}-}3\frac{{\theta }_i/T}{{\text{e}}^{{\theta }_i/T}-1}\right]\\ C_{vm}=-A\cdot B^{\prime }\cdot \frac{\text{d}{m}_s\left(T\right)}{\text{d}T}\\ {\theta }_i={\theta }_0\times \sqrt{E_i\left(0\right)/{\left(V_i\left(0\right)\right)}^{\left(2/3\right)}}/\sqrt{E_0\left(0\right)/{\left(V_0\left(0\right)\right)}^{\left(2/3\right)}}\end{array}$ (2-32)

where R is the gas constant; T is the temperature; ${\theta }_0$ , $E_0\left(0\right)$ and $V_0\left(0\right)$ are Debye temperature, potential energy and volume of primary AGC crystal, respectively.

When $\frac{T}{T_c}\le \frac{T_p}{T_c}$ , the temperature dependence of the atomic magnetic moment $m_s\left(T\right)$ can be expressed by Brillouin function

$\frac{m_s\left(T\right)}{m_s\left(0\right)}=\tanh \left[\left(\frac{m_s\left(T\right)}{m_s\left(0\right)}\right)/\left(\frac{T}{T_c}\right)\right]$ (2-33)

where $T_c$ is the Curie temperature, and $T_p$ is temperature at the experimental specific heat peak.

From this, it may be obtained that

$\frac{\text{d}{m}_s\left(T\right)}{\text{d}T}=\frac{m_s\left(0\right)\tanh \left[\left(\frac{m_s\left(T\right)}{m_s\left(0\right)}\right)/\left(\frac{T}{T_c}\right)\right]{\text{sech}}^2\left[\left(\frac{m_s\left(T\right)}{m_s\left(0\right)}\right)/\left(\frac{T}{T_c}\right)\right]}{T\left\{{\text{sech}}^2\left[\left(\frac{m_s\left(T\right)}{m_s\left(0\right)}\right)/\left(\frac{T}{T_c}\right)\right]-\left(\frac{T}{T_c}\right)\right\}}$ (2-34)

When $\frac{T}{T_c}\ge \frac{T_p}{T_c}$ , we have chosen the following equation to express the temperature dependence of $m_s(\; T\; )$

$\frac{m_s\left(T\right)}{m_s\left(0\right)}=\frac{1}{a+b{\text{e}}^{\left(T/T_c-T_p/T_c\right)}}$ (2-35)

From this, it may be obtained that

$\frac{\text{d}{m}_s\left(T\right)}{\text{d}T}=-\frac{m_s\left(0\right)\times b{\text{e}}^{\left(T/T_c-T_p/T_c\right)}}{T_c\times {\left[a+b{\text{e}}^{\left(T/T_c-T_p/T_c\right)}\right]}^2}$ (2-36)

The values of a and b constants can be obtained by the conditions of boundary at $T_p$ temperature

$\{\begin{array}{l}a=\frac{m_s\left(0\right)}{m_s\left(T\right)}-b\\ b=-\left(\frac{\text{d}{m}_s}{\text{d}T}\right)/\left(\frac{m_s^2\left(T\right)}{T_c{m}_s\left(0\right)}\right)\end{array}$ (2-37)

g) Equations of thermal expansion

The linear thermal expansion coefficient $\alpha \left(T\right)$ of magnetic AGC crystals includes two parts: the normal linear thermal expansion coefficient ${\alpha }_n\left(T\right)$ arisen from inharmonic lattice vibration and the magnetic linear thermal expansion coefficient ${\alpha }_m\left(T\right)$ arisen from magnetic transformation [31] [33] [34] [35] .

The ${\alpha }_n$ can be obtained by Grüneisen formula

${\alpha }_v\left(T\right)=\frac{C_v\left(T\right)}{3Q{\left[1-K\left(U/Q\right)\right]}^2}$ (2-38)

$U\left(T\right)={\displaystyle {\int }_0^T{C}_v\left(T\right)\text{d}T}$ (2-39)

The constants K and Q have been derived from $W\left(\psi ,r,\epsilon \right)$ function:

$K=\frac{2\text{π}\theta k_B{r}_0}{3hj}\left(2n-1\right)\sqrt{\frac{m}{n\left(n-1\right)E_c}}+\frac{1}{2}$ (2-40)

$Q=m{\left[\frac{4\text{π}\theta k_B{r}_0}{3hj}\right]}^2/\left[\frac{1}{3}\left[\frac{\theta k_B{r}_0}{hj}\cdot \sqrt{\frac{m}{E_c}\cdot \frac{n}{n-1}}+1\right]\right]$ (2-41)

Considering influence of magnetic electrons on the bond length, the Equation (2-9) should be modified as

$a\left(\psi ,r\right)=\left(\frac{1}{G_1}\right)\left\{2R\left(1\right)-\beta \lg \left[\frac{\left(n_c+A^{\prime }{m}_s\left(T\right)\right)}{\left({\displaystyle \underset{s}{\sum }{I}_s}\times {10}^{\left(G_1-G_s\right)a/\beta }\right)}\right]\right\}$ (2-42)

From this equation, the ${\alpha }_m\left(T\right)$ can be obtained by

${\alpha }_m\left(T\right)=-\frac{G_1}{a_{293}}\times \frac{\beta }{\ln 10}\times \frac{A^{\prime }}{n_c+A{m}_s\left(T\right)}\times \frac{\text{d}{m}_s\left(T\right)}{\text{d}T}$ (2-43)

where $m_s\left(T\right)$ is given by Equations ((2-33) and (2-35)); the $\frac{\text{d}{m}_s\left(T\right)}{\text{d}T}$ is given

by Equations ((2-34) and (2-36)). The value of $A^{\prime }$ constant can be obtained from an experimental ${\alpha }_m$ value at a temperature near the Curie point of primary AGC crystals.

The correctness of these equations in this chapter were demonstrated by agreement between theoretical and experimental results of the primary AGC Fe, Co, Ni, Cu, Ag, Au crystals, i.e., pure metals [32] [33] [34] [35] . These achievements lead to establish AG Gibbs energy $G_i\left(\psi ,r,\epsilon ,T_f,M_f,E_f,F_f\right)$ sequences and can make them respond to $T_f$ , $M_f$ , $E_f$ , and $F_f$ fields.

3) Equations of AGC Gibbs energy

Under the influence of temperature field, AG Gibbs energy $G_i\left(\psi ,r,\epsilon ,T_f,M_f,E_f,F_f\right)$ function may be rewritten simply as the $G_i\left(T\right)$ , which may be split into two parts: the AG temperature-independent potential energy $E_i\left(0\right)$ and temperature-dependent generalized vibration/action free energy $X_i^v\left(T\right)$ (Equation (3-1)). The AG enthalpy $H_i\left(T\right)$ of each characteristic crystal may be also split into two parts: the AG potential energy $E_i\left(0\right)$ and generalized vibration/action energy $U^v(T)$ (Equation (3-2)). The generalized vibration/action free energy $X_i^v\left(T\right)$ includes the generalized vibration energy $U_i^v\left(T\right)$ , which include Debye vibration energy $U_i^D\left(T\right)$ and attaching vibration energy $U_i^E\left(T\right)$ , and generalized vibration/action entropy energy $T{S}_i^v\left(T\right)$ , which include Debye vibration entropy energy $T{S}_i^D\left(T\right)$ and attaching action entropy energy $T{S}_i^E\left(T\right)$ . The generalized vibration/action heat capacity $C_{p,i}^v\left(T\right)$ includes Debye vibration heat capacity $C_{p,i}^D\left(T\right)$ and attaching action heat capacity $C_{p,i}^E\left(T\right)$ . The attaching $C_{p,i}^E\left(T\right)$ includes contributions of electron excitation, energy of formation of holes, variation of potential energy and volume expansion work with temperature, and so on. Therefore, the $G_i\left(T\right)$ is the AG holographic energy function with multi-level thermodynamic properties [45] [46] [47] [48] :

$\{\begin{array}{l}{G}_i^{Au}\left(T\right)=E_i^{Au}\left(0\right)+X_i^{Au.v}\left(T\right)\\ G_i^{Cu}\left(T\right)=E_i^{cu}\left(0\right)+X_i^{Cu.v}\left(T\right)\end{array}$ (3-1)

$\{\begin{array}{l}{H}_i^{Au}\left(T\right)=E_i^{Au}\left(0\right)+U_i^{Au.v}\left(T\right)\\ H_i^{Cu}\left(T\right)=E_i^{Cu}\left(0\right)+U_i^{Cu.v}\left(T\right)\end{array}$ (3-2)

$\{\begin{array}{l}{X}_i^{Au.v}\left(T\right)=U_i^{Au.v}\left(T\right)-T{S}_i^{Au.v}\left(T\right)\\ X_i^{Cu.v}\left(T\right)=U_i^{Cu.v}\left(T\right)-T{S}_i^{Cu.v}\left(T\right)\end{array}$ (3-3)

$\{\begin{array}{l}{X}_i^{Au.v}\left(T\right)=X_i^{Au.D}\left(T\right)+X_i^{Au.E}\left(T\right)\\ X_i^{Cu.v}\left(T\right)=X_i^{Cu.D}\left(T\right)+X_i^{Cu.E}\left(T\right)\end{array}$ (3-4)

$\{\begin{array}{l}{X}_i^{Au.D}\left(T\right)=U_i^{Au.D}\left(T\right)-T{S}_i^{Au.D}\left(T\right)\\ X_i^{Cu.D}\left(T\right)=U_i^{Cu.D}\left(T\right)-T{S}_i^{Cu.D}\left(T\right)\end{array}$ (3-5)

$\{\begin{array}{l}{X}_i^{Au.E}\left(T\right)=U_i^{Au.E}\left(T\right)-T{S}_i^{Au.E}\left(T\right)\\ X_i^{Cu.E}\left(T\right)=U_i^{Cu.E}\left(T\right)-T{S}_i^{Cu.E}\left(T\right)\end{array}$ (3-6)

$\{\begin{array}{l}{U}_i^{Au.v}\left(T\right)=U_i^{Au.D}\left(T\right)+U_i^{Au.E}\left(T\right)\\ U_i^{Cu.v}\left(T\right)=U_i^{Cu.D}\left(T\right)+U_i^{Cu.E}\left(T\right)\end{array}$ (3-7)

$\{\begin{array}{l}{S}_i^{Au.v}\left(T\right)=S_i^{Au.D}\left(T\right)+S_i^{Au.E}\left(T\right)\\ S_i^{Cu.v}\left(T\right)=S_i^{Cu.D}\left(T\right)+S_i^{Cu.E}\left(T\right)\end{array}$ (3-8)

$\{\begin{array}{l}{U}_i^{Au.D}\left(T\right)={\displaystyle \underset{0}{\overset{T}{\int }}{C}_{p.i}^{Au.D}\left(T\right)\text{d}T}\\ U_i^{Cu.D}\left(T\right)={\displaystyle \underset{0}{\overset{T}{\int }}{C}_{p.i}^{Cu.D}\left(T\right)\text{d}T}\end{array}$ (3-9)

$\{\begin{array}{l}{U}_i^{Au.E}\left(T\right)={\displaystyle \underset{0}{\overset{T}{\int }}{C}_{p.i}^{Au.E}\left(T\right)\text{d}T}\\ U_i^{Cu.E}\left(T\right)={\displaystyle \underset{0}{\overset{T}{\int }}{C}_{p.i}^{Cu.E}\left(T\right)\text{d}T}\end{array}$ (3-10)

$\{\begin{array}{l}{S}_i^{Au.D}\left(T\right)={\displaystyle \underset{0}{\overset{T}{\int }}\frac{C_{p.i}^{Au.D}\left(T\right)}{T}\text{d}T}\\ S_i^{Cu.D}\left(T\right)={\displaystyle \underset{0}{\overset{T}{\int }}\frac{C_{p.i}^{Cu.D}\left(T\right)}{T}\text{d}T}\end{array}$ (3-11)

$\{\begin{array}{l}{S}_i^{Au.E}\left(T\right)={\displaystyle \underset{0}{\overset{T}{\int }}\frac{C_{p.i}^{Au.E}\left(T\right)}{T}\text{d}T}\\ S_i^{Cu.E}\left(T\right)={\displaystyle \underset{0}{\overset{T}{\int }}\frac{C_{p.i}^{Cu.E}\left(T\right)}{T}\text{d}T}\end{array}$ (3-12)

$\{\begin{array}{l}{\theta }_i^{Au}={\theta }_0^{Au}\times \sqrt{E_i^{Au}\left(0\right)/{\left(V_i^{Au}\left(0\right)\right)}^{\left(2/3\right)}}/\sqrt{E_0^{Au}\left(0\right)/{\left(V_0^{Au}\left(0\right)\right)}^{\left(2/3\right)}}\\ {\theta }_i^{Cu}={\theta }_I^{Cu}\times \sqrt{E_i^{Cu}\left(0\right)/{\left(V_i^{Cu}\left(0\right)\right)}^{\left(2/3\right)}}/\sqrt{E_{12}^{Cu}\left(0\right)/{\left(V_{12}^{Cu}\left(0\right)\right)}^{\left(2/3\right)}}\end{array}$ (3-13)

$\{\begin{array}{l}{C}_{p.i}^{Au.D}\left(T\right)=9R{\left(\frac{T}{{\theta }_i^{Au}}\right)}^3{\displaystyle {\int }_0^{{\theta }_i^{Au}/T}\frac{{\text{e}}^x{x}^4}{{\left({\text{e}}^x-1\right)}^2}\text{d}x}\\ C_{p.i}^{Cu.D}\left(T\right)=9R{\left(\frac{T}{{\theta }_i^{Cu}}\right)}^3{\displaystyle {\int }_0^{{\theta }_i^{Cu}/T}\frac{{\text{e}}^x{x}^4}{{\left({\text{e}}^x-1\right)}^2}\text{d}x}\end{array}$ (3-14)

$\{\begin{array}{l}{C}_{p.0}^{Au.E}\left(T\right)=1.858993\times {10}^{-10}T+4.222584\times {10}^{-6}{T}^2\\ C_{p.12}^{Au.E}\left(T\right)=2.763134\times {10}^{-4}T+5.481502\times {10}^{-6}{T}^2\end{array}$ (3-15)

$\{\begin{array}{l}{C}_{p.12}^{Cu.E}\left(T\right)=0.00168035T+2.826912\times {10}^{-6}{T}^2\\ C_{p.0}^{Cu.E}\left(T\right)=0.00564614T+3.632353\times {10}^{-6}{T}^2\end{array}$ (3-16)

where, ${\theta }_0^{Au}=165\text{K}$ , ${\theta }_{12}^{Au}=175.21\text{K}$ , ${\theta }_{12}^{Cu}=343\text{K}$ , ${\theta }_0^{Cu}=343.07\text{K}$ .

It should be emphasized that the AG holographic warehouse of fcc-based lattice Au-Cu system has been established only based on experimental enthalpies and volumes of pure Au and Cu metals, as well as AuCu and AuCu₃ compounds (Figure 15)!

Figure 15. Calculated road of AG-WOP methods. AG-holographic information includes AG potential energies ( $E_i^{\alpha }$ ) and volumes ( $V_i^{\alpha }$ ) of basic properties, as well as AG valence electron structures ( ${\psi }_i^{\alpha }$ ) of basic structure; the AG physical property sequences (mechanical, magnetical, electric, optical and so on) obtained by AG-AGC valence bond theory; the AG Gibbs energy sequences with multi-level thermodynamic properties may be obtained by AG-AGC thermodynamics.

4) Equations of AG arranging thermodynamics

The AG arranging phase theory includes the AG arranging thermodynamics, which describes the AG arranging Gibbs energy band structures of alloy genes occupied at the Gibbs energy levels in the Gibbs energy space and thermodynamic properties of alloy phases; the AG arranging crystallography, which describes the AG arranging crystalline structures of alloy genes occupied at the lattice points in the crystalline space and crystalline symmetries of alloy phases; and AG arranging valence bond theory, which describes the AG arranging bond network structures of alloy genes occupied at the lattice points in the crystalline bond network space and physical and mechanical properties of alloy phases. Here, the main equations of the AG arranging thermodynamics are presented only, taking AuCu-type sublattice phase system as an example.

a) Equations chain of AG-GET function

The AG Gibbs energy transmissive AG-GET $\Omega \left(x,T,\sigma \right)$ function produces alloy phase holographic energy Gibbs energy $G\left(x,T,\sigma \right)$ function and configurational entropy $S^c\left(x,T,\sigma \right)$ function [45] :

$\Omega \left(x,T,\sigma \right)=g\left(x_i^{Au}\left(x,T,\sigma \right),x_i^{Cu}\left(x,T,\sigma \right)\right)\times \exp \left[-G^{\ast }\left(x,T,\sigma \right)/\left(kT\right)\right]$ (4-1)

$G\left(x,T,\sigma \right)=-k_{B}T\times \ln \left(\Omega \left(x,T,\sigma \right)\right)=G^{\ast }\left(x,T,\sigma \right)-T{S}^c\left(x,T,\sigma \right)$ (4-2)

$S^c\left(x,T,\sigma \right)=k_B\ln \left(g\left(x_i^{Au}\left(x,T,\sigma \right),x_i^{Cu}\left(x,T,\sigma \right)\right)\right)$ (4-3)

AG arranging characteristic (chemical potential) Gibbs energy $G^{*}\left(x,T\right)$ function of alloy phases:

$G^{\ast }\left(x,T,\sigma \right)={\displaystyle \underset{i=0}{\overset{12}{\sum }}\left(x_i^{Au}\left(x,T,\sigma \right)\times G_i^{Au}\left(T\right)+x_i^{Cu}\left(x,T,\sigma \right)\times G_i^{Cu}\left(T\right)\right)}$ (4-4)

AG arranging degeneracy function of alloy phases:

$\begin{array}{l}g\left(x_i^{Au}\left(x,T,\sigma \right),x_i^{Cu}\left(x,T,\sigma \right)\right)\\ =\frac{\left[N{\displaystyle \underset{i=0}{\overset{I}{\sum }}\left(x_i^{Au\left(1\right)}+x_i^{Cu\left(1\right)}\right)}\right]!}{\left({\displaystyle \underset{i=0}{\overset{I}{\sum }}N{x}_i^{Au\left(1\right)}}\right)!\left({\displaystyle \underset{i=0}{\overset{I}{\sum }}N{x}_i^{Cu\left(1\right)}}\right)!}\times \frac{\left[N{\displaystyle \underset{i=0}{\overset{I}{\sum }}\left(x_i^{Au\left(2\right)}+x_i^{Cu\left(2\right)}\right)}\right]!}{\left({\displaystyle \underset{i=0}{\overset{I}{\sum }}N{x}_i^{Au\left(2\right)}}\right)!\left({\displaystyle \underset{i=0}{\overset{I}{\sum }}N{x}_i^{Cu\left(2\right)}}\right)!}\end{array}$ (4-5)

Concentration equations of the AG $A_i^{Au}$ and $A_i^{Cu}$ at the (1) and (2) sublattice points:

$\{\begin{array}{l}{x}_i^{Au}\left(x,T,\sigma \right)=x_i^{Au\left(1\right)}\left(x,T,\sigma \right)+x_i^{Au\left(2\right)}\left(x,T,\sigma \right)\\ x_i^{Cu}\left(x,T,\sigma \right)=x_i^{Cu\left(1\right)}\left(x,T,\sigma \right)+x_i^{Cu\left(2\right)}\left(x,T,\sigma \right)\end{array}$ (4-6)

Concentration equations of the AG $A_i^{Au}$ and $A_i^{Cu}$ occupied at the (1) and (2) sublattice points:

$\{\begin{array}{l}{x}_i^{Au\left(1\right)}\left(x,T,\sigma \right)=v^{\left(1\right)}{P}_{Au}^{\left(1\right)}{\omega }_i^{\left(1\right)}\\ x_i^{Cu\left(1\right)}\left(x,T,\sigma \right)=v^{\left(1\right)}{P}_{Cu}^{\left(1\right)}{\omega }_i^{\left(1\right)}\\ x_i^{Au\left(2\right)}\left(x,T,\sigma \right)=v^{\left(2\right)}{P}_{Au}^{\left(2\right)}{\omega }_i^{\left(2\right)}\\ x_i^{Cu\left(2\right)}\left(x,T,\sigma \right)=v^{\left(2\right)}{P}_{Cu}^{\left(2\right)}{\omega }_i^{\left(2\right)}\end{array}$ (4-7)

Probability equations of Au and Cu atoms occupied at (1) and (2) sublattice points:

$\{\begin{array}{l}{P}_{Au}^{\left(1\right)}=x_{Au}+\left(1-v^{\left(1\right)}\right)\sigma ,\\ P_{Cu}^{\left(1\right)}=1-P_{Au}^{\left(1\right)}=1-\left[x_{Au}+\left(1-v^{\left(1\right)}\right)\sigma \right]\\ P_{Au}^{\left(2\right)}=x_{Au}-v^{\left(1\right)}\sigma ,\\ P_{Cu}^{\left(2\right)}=1-P_{Au}^{\left(2\right)}=1-\left(x_{Au}-v^{\left(1\right)}\sigma \right).\end{array}$ (4-8)

Fractions of the (1) and (2) sublattice points for AuCu-type sublattice phase:

$\{\begin{array}{l}{v}^{\left(1\right)}=N^{\left(1\right)}/N=1/2\\ v^{\left(2\right)}=N^{\left(2\right)}/N=1/2\end{array}$ (4-9)

Definition of order degree:

$\sigma =\frac{P_{Au}^{\left(1\right)}-x_{Au}}{1-v^{\left(1\right)}}=\frac{P_{Cu}^{\left(2\right)}-x_{Cu}}{1-v^{\left(2\right)}}=\frac{x_i^{Au\left(1\right)}}{v^{\left(1\right)}{v}^{\left(2\right)}{\omega }_i^{\left(1\right)}}-\frac{x_{Au}}{v^{\left(2\right)}}=\frac{x_i^{Cu(2)}}{v^{\left(1\right)}{v}^{\left(2\right)}{\omega }_i^{\left(2\right)}}-\frac{x_{Cu}}{v^{\left(1\right)}}$ (4-10)

The probability ${\omega }_i^{\left(1\right)}$ equation of the coordinative cluster $\left[\left(I-i\right)\text{Au},i\text{Cu}\right]$ surrounding the (1)-sublattice point:

${\omega }_i^{\left(1\right)}={\omega }_{i^{\prime }}^{\left(2\right)}\times {\omega }_{\left(i-i^{\prime }\right)}^{\left(1\right)}={\displaystyle \underset{i^{\prime }=k}{\overset{K}{\sum }}\left[C_8^{i^{\prime }}{\left(P_{Au}^{\left(2\right)}\right)}^{\left(8-i^{\prime }\right)}{\left(P_{Cu}^{\left(2\right)}\right)}^{\left(i^{\prime }\right)}\times C_4^{\left(i-i^{\prime }\right)}{\left(P_{Au}^{\left(1\right)}\right)}^{\left[4-\left(i-i^{\prime }\right)\right]}{\left(P_{Cu}^{\left(1\right)}\right)}^{\left(i-i^{\prime }\right)}\right]}$ (4-11)

In the ${\omega }_i^{\left(1\right)}$ equation, the probability ${\omega }_{i^{\prime }}^{\left(2\right)}$ of coordinative cluster $\left[\left(8-i^{\prime }\right)\text{Au},i^{\prime }\text{Cu}\right]$ occupied at the (2)-sublattice points, and the probability ${\omega }_{\left(i-i^{\prime }\right)}^{\left(1\right)}$ of coordinative cluster $\left[4-\left(i-i^{\prime }\right)\text{Au},\left(i-i^{\prime }\right)\text{Cu}\right]$ occupied at the (1)-sublattice points:

$\{\begin{array}{l}{\omega }_{i^{\prime }}^{\left(2\right)}={\displaystyle \underset{i^{\prime }=k}{\overset{K}{\sum }}\left[C_8^{i^{\prime }}{\left(P_{Au}^{\left(2\right)}\right)}^{\left(8-i^{\prime }\right)}{\left(P_{Cu}^{\left(2\right)}\right)}^{\left(i^{\prime }\right)}\right]}\\ {\omega }_{\left(i-i^{\prime }\right)}^{\left(1\right)}={\displaystyle \underset{i^{\prime }=k}{\overset{K}{\sum }}\left[C_4^{\left(i-i^{\prime }\right)}{\left(P_{Au}^{\left(1\right)}\right)}^{\left[4-\left(i-i^{\prime }\right)\right]}{\left(P_{Cu}^{\left(1\right)}\right)}^{\left(i-i^{\prime }\right)}\right]}\end{array}$ (4-12)

where, the i’ is the number of the Cu atoms occupied at the (2)-sublattice points.

The probability ${\omega }_i^{\left(2\right)}$ equation of coordinative cluster $\left[\left(I-i\right)\text{Au},i\text{Cu}\right]$ surrounding the (2)-sublattice point:

${\omega }_i^{\left(2\right)}={\omega }_{i^{\prime }}^{\left(1\right)}\times {\omega }_{\left(i-i^{\prime }\right)}^{\left(2\right)}={\displaystyle \underset{i^{\prime }=k}{\overset{K}{\sum }}\left[C_8^{i^{\prime }}{\left(P_{Au}^{\left(1\right)}\right)}^{\left(8-i^{\prime }\right)}{\left(P_{Cu}^{\left(1\right)}\right)}^{\left(i^{\prime }\right)}\times C_4^{\left(i-i^{\prime }\right)}{\left(P_{Au}^{\left(2\right)}\right)}^{\left[4-\left(i-i^{\prime }\right)\right]}{\left(P_{Cu}^{\left(2\right)}\right)}^{\left(i-i^{\prime }\right)}\right]}$ (4-13)

where, the i’ is the number of the Cu atoms occupied at the (1)-sublattice points

In the ${\omega }_i^{(2)}$ equation, the probability ${\omega }_{i\text{'}}^{(1)}$ of coordinative cluster $\left[\left(8-i^{\prime }\right)\text{Au},i^{\prime }\text{Cu}\right]$ occupied at the (1)-sublattice points, and the probability ${\omega }_{(i-i\text{'})}^{(2)}$ of coordinative cluster $\left[4-\left(i-i^{\prime }\right)\text{Au},\left(i-i^{\prime }\right)\text{Cu}\right]$ occupied at the (2)-sublattice points:

$\{\begin{array}{l}{\omega }_{i^{\prime }}^{\left(1\right)}={\displaystyle \underset{i^{\prime }=k}{\overset{K}{\sum }}\left[C_8^{i^{\prime }}{\left(P_{Au}^{\left(1\right)}\right)}^{\left(8-i^{\prime }\right)}{\left(P_{Cu}^{\left(1\right)}\right)}^{\left(i^{\prime }\right)}\right]}\\ {\omega }_{\left(i-i^{\prime }\right)}^{\left(2\right)}={\displaystyle \underset{i^{\prime }=k}{\overset{K}{\sum }}\left[C_4^{\left(i-i^{\prime }\right)}{\left(P_{Au}^{\left(2\right)}\right)}^{\left[4-\left(i-i^{\prime }\right)\right]}{\left(P_{Cu}^{\left(2\right)}\right)}^{\left(i-i^{\prime }\right)}\right]}\end{array}$ (4-14)

In the Equations ((4-9) to (4-11)), k and K take their values, according to

$\{\begin{array}{l}k=0,K=i,\text{if}i\le 4;\\ k=i-4,K=i,\text{if}4<i\le 8;\\ k=i-4,K=8,\text{if}8<i\le 12.\end{array}$ (4-15)

The maximum order degree ${\sigma }_{\max }$ as function of composition $x_{Au}$ :

For the stoichiometric alloy, due to $P_{Au}^{\left(1\right)}=1$ , ${\sigma }_{\max }=1$ ;

For the alloys with $x_{Au}<v^{\left(1\right)}$ , due to $P_{Au}^{\left(2\right)}=0$ , ${\sigma }_{\max }\left(x\right)=\frac{x_{Au}}{v^{\left(1\right)}}$ ;

For the alloys with $x_{Au}>v^{\left(1\right)}$ , due to $P_{Au}^{\left(1\right)}=1$ , ${\sigma }_{\max }\left(x\right)=\frac{\left(1-x_{Au}\right)}{\left(1-v^{\left(1\right)}\right)}$ .

The methods for calculating ${\omega }_i^{\left(1\right)}$ and ${\omega }_i^{\left(2\right)}$ are listed in Table 10, Table 11.

The configurational entropy equations of the alloy and components:

$\{\begin{array}{l}{S}^c\left(x,T,\sigma \right)=k_B\ln \left(g\left(x_i^{Au}\left(x,T,\sigma \right),x_i^{Cu}\left(x,T,\sigma \right)\right)\right)\\ =x_{Au}{S}^{Au.c}\left(x,T,\sigma \right)+x_{Cu}{S}^{Cu.c}\left(x,T,\sigma \right)\\ S^{Au.c}=-\frac{R}{x_{Au}}[\left({\displaystyle \underset{i=0}{\overset{I}{\sum }}{x}_i^{Au\left(1\right)}\left(x,T,\sigma \right)}\right)\ln \left({\displaystyle \underset{i=0}{\overset{I}{\sum }}\left(P_{Au}^{\left(1\right)}{\omega }_i^{\left(1\right)}\right)}\right)\\ +\left({\displaystyle \underset{i=0}{\overset{I}{\sum }}{x}_i^{Au\left(2\right)}\left(x,T,\sigma \right)}\right)\ln \left({\displaystyle \underset{i=0}{\overset{I}{\sum }}\left(P_{Au}^{\left(2\right)}{\omega }_i^{\left(2\right)}\right)}\right)]\\ S^{Cu.c}=-\frac{R}{x_{Cu}}[\left({\displaystyle \underset{i=0}{\overset{I}{\sum }}{x}_i^{Cu\left(1\right)}\left(x,T,\sigma \right)}\right)\ln \left({\displaystyle \underset{i=0}{\overset{I}{\sum }}\left(P_{Cu}^{\left(1\right)}{\omega }_i^{\left(1\right)}\right)}\right)\\ +\left({\displaystyle \underset{i=0}{\overset{I}{\sum }}{x}_i^{Cu\left(2\right)}\left(x,T,\sigma \right)}\right)\ln \left({\displaystyle \underset{i=0}{\overset{I}{\sum }}\left(P_{Cu}^{\left(2\right)}{\omega }_i^{\left(2\right)}\right)}\right)]\end{array}$ (4-16)

The mole potential energy equations of the alloy and components:

$\{\begin{array}{l}E\left(x,T,\sigma \right)=x_{Au}{E}^{Au}\left(x,T,\sigma \right)+x_{Cu}{E}^{Cu}\left(x,T,\sigma \right)\\ E^{Au}\left(x,T,\sigma \right)=\left(1/x_{Au}\right){\displaystyle \underset{i=0}{\overset{I=12}{\sum }}{x}_i^{Au}\left(x,T,\sigma \right)E_i^{Au}\left(0\right)}\\ E^{Cu}\left(x,T,\sigma \right)=\left(1/x_{Cu}\right){\displaystyle \underset{i=0}{\overset{I=12}{\sum }}{x}_i^{Cu}\left(x,T,\sigma \right)E_i^{Cu}\left(0\right)}\end{array}$ (4-17)

The Debye vibration energy equations of alloy and components:

$\{\begin{array}{l}{U}^D\left(x,T,\sigma \right)=x_{Au}{U}^{Au.D}\left(x,T,\sigma \right)+x_{Cu}{X}^{Cu.D}\left(x,T,\sigma \right)\\ U^{Au.D}\left(x,T,\sigma \right)=\left(1/x_{Au}\right){\displaystyle \underset{i=0}{\overset{I=12}{\sum }}{x}_i^{Au}\left(x,T,\sigma \right)}{\displaystyle {\int }_0^T{C}_{p,i}^{Au.D}\left(T\right)\text{d}T}\\ U^{Cu.D}\left(x,T,\sigma \right)=\left(1/x_{Cu}\right){\displaystyle \underset{i=0}{\overset{I=12}{\sum }}{x}_i^{Cu}\left(x,T,\sigma \right)}{\displaystyle {\int }_0^T{C}_{p,i}^{Cu.D}\left(T\right)\text{d}T}\end{array}$ (4-18)

The attaching action energy equations of alloy and components:

$\{\begin{array}{l}{U}^E\left(x,T,\sigma \right)=x_{Au}{U}^{Au.E}\left(x,T,\sigma \right)+x_{Cu}{X}^{Cu.E}\left(x,T,\sigma \right)\\ U^{Au.E}\left(x,T,\sigma \right)=\left(1/x_{Au}\right){\displaystyle \underset{i=0}{\overset{I=12}{\sum }}{x}_i^{Au}\left(x,T,\sigma \right)}{\displaystyle {\int }_0^T{C}_{p,i}^{Au.E}\left(T\right)\text{d}T}\\ U^{Cu.E}\left(x,T,\sigma \right)=\left(1/x_{Cu}\right){\displaystyle \underset{i=0}{\overset{I=12}{\sum }}{x}_i^{Cu}\left(x,T,\sigma \right)}{\displaystyle {\int }_0^T{C}_{p,i}^{Cu.E}\left(T\right)\text{d}T}\end{array}$ (4-19)

Table 10. The methods for calculating the probability ${\omega }_i^{\left(1\right)}$ of coordinative cluster $\left[\left(I-i\right)\text{Au},i\text{Cu}\right]$ surrounding the (1)-sublattice point in AuCu-type order alloys.

Table 11. The methods for calculating the probability ${\omega }_i^{\left(2\right)}$ of coordinative cluster $\left[\left(I-i\right)\text{Au},i\text{Cu}\right]$ surrounding the (2)-sublattice point in AuCu-type order alloys.

The generalized action energy equations of the alloy and components:

$\{\begin{array}{l}{U}^v\left(x,T,\sigma \right)=x_{Au}{U}^{Au.v}\left(x,T,\sigma \right)+x_{Cu}{U}^{Cu.v}\left(x,T,\sigma \right)\\ U^{Au.v}\left(x,T,\sigma \right)=\left(1/x_{Au}\right){\displaystyle \underset{i=0}{\overset{I=12}{\sum }}{x}_i^{Au}\left(x,T,\sigma \right)U_i^{Au.v}\left(T\right)}\\ U^{Cu.v}\left(x,T,\sigma \right)=\left(1/x_{Cu}\right){\displaystyle \underset{i=0}{\overset{I=12}{\sum }}{x}_i^{Cu}\left(x,T,\sigma \right)U_i^{Cu.v}\left(T\right)}\end{array}$ (4-20)

The Debye vibration entropy equations of the alloy and components:

$\{\begin{array}{l}{S}^D\left(x,T,\sigma \right)=x_{Au}{S}^{Au.D}\left(x,T,\sigma \right)+x_{Cu}{S}^{Cu.D}\left(x,T,\sigma \right)\\ S^{Au.D}\left(x,T,\sigma \right)=\left(1/x_{Au}\right){\displaystyle \underset{i=0}{\overset{I=12}{\sum }}{x}_i^{Au}\left(x,T,\sigma \right)}{\displaystyle {\int }_0^T\frac{C_{p,i}^{Au.D}\left(T\right)}{T}\text{d}T}\\ S^{Cu.D}\left(x,T,\sigma \right)=\left(1/x_{Cu}\right){\displaystyle \underset{i=0}{\overset{I=12}{\sum }}{x}_i^{Cu}\left(x,T,\sigma \right)}{\displaystyle {\int }_0^T\frac{C_{p,i}^{Cu.D}\left(T\right)}{T}\text{d}T}\end{array}$ (4-21)

The attaching action entropy equations of the alloy and components:

$\{\begin{array}{l}{S}^E\left(x,T,\sigma \right)=x_{Au}{S}^{Au.E}\left(x,T,\sigma \right)+x_{Cu}{S}^{Cu.E}\left(x,T,\sigma \right)\\ S^{Au.E}\left(x,T,\sigma \right)=\left(1/x_{Au}\right){\displaystyle \underset{i=0}{\overset{I=12}{\sum }}{x}_i^{Au}\left(x,T,\sigma \right)}{\displaystyle {\int }_0^T\frac{C_{p,i}^{Au.E}\left(T\right)}{T}\text{d}T}\\ S^{Cu.E}\left(x,T,\sigma \right)=\left(1/x_{Cu}\right){\displaystyle \underset{i=0}{\overset{I=12}{\sum }}{x}_i^{Cu}\left(x,T,\sigma \right)}{\displaystyle {\int }_0^T\frac{C_{p,i}^{Cu.E}\left(T\right)}{T}\text{d}T}\end{array}$ (4-22)

The generalized action entropy equations of the alloy and components:

$\{\begin{array}{l}{S}^v\left(x,T,\sigma \right)=x_{Au}{S}^{Au.v}\left(x,T,\sigma \right)+x_{Cu}{S}^{Cu.v}\left(x,T,\sigma \right)\\ S^{Au.v}\left(x,T,\sigma \right)=\left(1/x_{Au}\right){\displaystyle \underset{i=0}{\overset{I=12}{\sum }}{x}_i^{Au}\left(x,T,\sigma \right)S_i^{Au.v}\left(T\right)}\\ S^{Cu.v}\left(x,T,\sigma \right)=\left(1/x_{Cu}\right){\displaystyle \underset{i=0}{\overset{I=12}{\sum }}{x}_i^{Cu}\left(x,T,\sigma \right)S_i^{Cu.v}\left(T\right)}\end{array}$ (4-23)

The generalized action free energy equations of the alloy and components:

$\{\begin{array}{l}{X}^v\left(x,T,\sigma \right)=x_{Au}{X}^{Au.v}(x,T,\sigma )+x_{Cu}{X}^{Cu.v}\left(x,T,\sigma \right)\\ X^{Au.v}\left(x,T,\sigma \right)=\left(1/x_{Au}\right){\displaystyle \underset{i=0}{\overset{I=12}{\sum }}{x}_i^{Au}\left(x,T,\sigma \right)X_i^{Au.v}\left(T\right)}\\ X^{Cu.v}\left(x,T,\sigma \right)=\left(1/x_{Cu}\right){\displaystyle \underset{i=0}{\overset{I=12}{\sum }}{x}_i^{Cu}\left(x,T,\sigma \right)X_i^{Cu.v}\left(T\right)}\end{array}$ (4-24)

The enthalpy equations of the alloy and components:

$\{\begin{array}{l}H\left(x,T,\sigma \right)=x_{Au}{H}^{Au}\left(x,T,\sigma \right)+x_{Cu}{H}^{Cu}\left(x,T,\sigma \right)\\ H^{Au}\left(x,T,\sigma \right)=\left(1/x_{Au}\right){\displaystyle \underset{i=0}{\overset{I=12}{\sum }}{x}_i^{Au}\left(x,T,\sigma \right)}\left(E_i^{Au}\left(0\right)+U_i^{Au.v}\left(T\right)\right)\\ H^{Cu}\left(x,T,\sigma \right)=\left(1/x_{Cu}\right){\displaystyle \underset{i=0}{\overset{I=12}{\sum }}{x}_i^{Cu}\left(x,T,\sigma \right)}\left(E_i^{Cu}\left(0\right)+U_i^{Cu.v}\left(T\right)\right)\end{array}$ (4-25)

The characteristic Gibbs energy equations of the alloy and components:

$\{\begin{array}{l}{G}^{*}\left(x,T,\sigma \right)=x_{Au}{G}^{*Au}\left(x,T,\sigma \right)+x_{Cu}{G}^{*Cu}\left(x,T,\sigma \right)\\ G^{*Au}\left(x,T,\sigma \right)=\left(1/x_{Au}\right){\displaystyle \underset{i=0}{\overset{I=12}{\sum }}{x}_i^{Au}\left(x,T,\sigma \right)G_i^{Au}\left(T\right)}\\ G^{*Cu}\left(x,T,\sigma \right)=\left(1/x_{Cu}\right){\displaystyle \underset{i=0}{\overset{I=12}{\sum }}{x}_i^{Cu}\left(x,T,\sigma \right)G_i^{Cu}\left(T\right)}\end{array}$ (4-26)

The Gibbs energy equations of the alloy and components:

$\{\begin{array}{l}G\left(x,T,\sigma \right)=x_{Au}{G}^{Au}\left(x,T,\sigma \right)+x_{Cu}{G}^{Cu}\left(x,T,\sigma \right)\\ G^{Au}\left(x,T,\sigma \right)=\left(1/x_{Au}\right)\left(G^{*Au}\left(x,T,\sigma \right)-T{S}^{Au.c}\left(x,T,\sigma \right)\right)\\ G^{Cu}\left(x,T,\sigma \right)=\left(1/x_{Cu}\right)\left(G^{*Cu}\left(x,T,\sigma \right)-T{S}^{Cu.c}\left(x,T,\sigma \right)\right)\end{array}$ (4-27)

The mixed heat capacity equations of the alloy:

$\{\begin{array}{l}\Delta C_p^m\left(x,T,\sigma \right)=\Delta C_p^{m.v}\left(x,T,\sigma \right)+\Delta C_p^{m.s}\left(x,T,\sigma \right)\\ \Delta C_p^{m.v}\left(x,T,\sigma \right)={\displaystyle \underset{i=0}{\overset{I}{\sum }}{x}_i^{Au}\left(x,T,\sigma \right)\left(C_{p,i}^{Au.v}\left(T\right)-C_{p,0}^{Au.v}\left(T\right)\right)}\\ +{\displaystyle \underset{i=0}{\overset{I}{\sum }}{x}_i^{Cu}\left(x,T,\sigma \right)\left(C_{p,i}^{Cu.v}\left(T\right)-C_{p,12}^{Cu.v}\left(T\right)\right)}\\ \Delta C_p^{m.s}\left(x,T,\sigma \right)={\displaystyle \underset{i=0}{\overset{I}{\sum }}\frac{\partial x_i^{Au}\left(x,T,\sigma \right)}{\partial \sigma }\cdot \frac{\text{d}\sigma }{\text{d}T}\left(H_i^{Au}\left(T\right)-H_0^{Au}\left(T\right)\right)}\\ +{\displaystyle \underset{i=0}{\overset{I}{\sum }}\frac{\partial x_i^{Cu}\left(x,T,\sigma \right)}{\partial \sigma }\cdot \frac{\text{d}\sigma }{\text{d}T}\left(H_i^{Cu}\left(T\right)-H_{12}^{Cu}\left(T\right)\right)}\end{array}$ (4-28)

The activity equations of the components:

$\{\begin{array}{l}{a}_{Au}\left(x,T,\sigma \right)=x_{Au}\cdot \exp \left(\Delta G_{Au}^{*m}\left(x,T,\sigma \right)/RT\right)\\ a_{Cu}\left(x,T,\sigma \right)=x_{Cu}\cdot \exp \left(\Delta G_{Cu}^{*m}\left(x,T,\sigma \right)/RT\right)\\ \Delta G^{*Au}\left(x,T,\sigma \right)=\left(1/x_{Au}\right){\displaystyle \underset{i=0}{\overset{I=12}{\sum }}{x}_i^{Au}\left(x,T,\sigma \right)\left(G_i^{Au}\left(T\right)-G_0^{Au}\left(T\right)\right)}\\ \Delta G^{*Cu}\left(x,T,\sigma \right)=\left(1/x_{Cu}\right){\displaystyle \underset{i=0}{\overset{I=12}{\sum }}{x}_i^{Cu}\left(x,T,\sigma \right)}\left(G_i^{Cu}\left(T\right)-G_I^{Cu}\left(T\right)\right)\end{array}$ (4-29)

The volume equations of the alloy and components:

$\{\begin{array}{l}v\left(x,T,\sigma \right)=x_{Au}{v}^{Au}\left(x,T,\sigma \right)+x_{Cu}{v}^{Cu}\left(x,T,\sigma \right)\\ v^{Au}\left(x,T,\sigma \right)=\left(1/x_{Au}\right){\displaystyle \underset{i=0}{\overset{I=12}{\sum }}{x}_i^{Au}\left(x,T,\sigma \right)v_i^{Au}\left(T\right)}\\ v^{Cu}\left(x,T,\sigma \right)=\left(1/x_{Cu}\right){\displaystyle \underset{i=0}{\overset{I=12}{\sum }}{x}_i^{Cu}\left(x,T,\sigma \right)v_i^{Cu}\left(T\right)}\end{array}$ (4-30)

The mixed volume thermal expansion coefficient equations of the alloy:

$\{\begin{array}{l}\Delta {\beta }^m\left(x,T,\sigma \right)=\Delta {\beta }^{m.v}\left(x,T,\sigma \right)+\Delta {\beta }^{m.s}\left(x,T,\sigma \right)\\ \Delta {\beta }^{m.v}\left(x,T,\sigma \right)={\displaystyle \underset{i=0}{\overset{I}{\sum }}{x}_i^{Au}\left(x,T,\sigma \right)\left({\beta }_i^{Au.v}\left(T\right)-{\beta }_0^{Au.v}\left(T\right)\right)}\\ +{\displaystyle \underset{i=0}{\overset{I}{\sum }}{x}_i^{Cu}\left(x,T,\sigma \right)\left({\beta }_i^{Cu.v}\left(T\right)-{\beta }_{12}^{Cu.v}\left(T\right)\right)}\\ \Delta {\beta }^{m.s}\left(x,T,\sigma \right)={\displaystyle \underset{i=0}{\overset{I}{\sum }}\left(\frac{1}{v_i^{Au}\left(T\right)}\right)\frac{\partial x_i^{Au}\left(x,T,\sigma \right)}{\partial \sigma }\cdot \frac{\text{d}\sigma }{\text{d}T}\left(v_i^{Au}\left(T\right)-v_0^{Au}\left(T\right)\right)}\\ +{\displaystyle \underset{i=0}{\overset{I}{\sum }}\left(\frac{1}{v_i^{Cu}\left(T\right)}\right)\frac{\partial x_i^{Cu}\left(x,T,\sigma \right)}{\partial \sigma }\cdot \frac{\text{d}\sigma }{\text{d}T}\left(v_i^{Cu}\left(T\right)-v_{12}^{Cu}\left(T\right)\right)}\end{array}$ (4-31)

It should be emphasized that the new $\Omega \left(x,T,\sigma \right)$ -function can be used to describe both long-distance order and short-distance order states of alloy phases; and that, the thermodynamic properties of alloy components need no representatives by the partial molar properties. We have proved that the partial molar properties cannot represent the average properties of the corresponding components [21] . It should be emphasized again that the equilibrium EAHNP and subequilibrium SAHNP path charts on disordering AuCu( $A_8^{Au}{A}_4^{Cu}$ ) compound (Figures 16-27), the EAHNP phase diagrams of Au₃Cu-, AuCu- and AuCu₃-type sublattice systems (Figures 28-30) and the EAHNP and SAHNP phase diagrams of Au-Cu system (Figure 31, Figure 32) have been established only based on experimental enthalpies and volumes of pure Au and Cu metals, as well as AuCu and AuCu₃ compounds! The composition ranges of long range ordered phases are shown in Figure 33.

Figure 16. EAHNP and SAHNP charts on disordering AuCuI( $A_8^{Au}{A}_4^{Cu}$ ). (a) Three-dimensional mixed Gibbs energy $\Delta G^m-T-\sigma$ EAHNP chart with $\Delta G_e^m-T$ and $\Delta G_s^m-T$ paths; (b) $\Delta G_e^m-T$ and $\Delta G_s^m-T$ paths.

Figure 17. EAHNP and SAHNP charts on disordering AuCuI( $A_8^{Au}{A}_4^{Cu}$ ). (a) Three-dimensional mixed Gibbs energy $\Delta G^{\ast m}-T-\sigma$ EAHNP chart with $\Delta G_e^{\ast m}-T$ and $\Delta G_s^{\ast m}-T$ paths; (b) $\Delta G_e^{\ast m}-T$ and $\Delta G_s^{\ast m}-T$ paths.

Figure 18. EAHNP and SAHNP charts on disordering AuCuI( $A_8^{Au}{A}_4^{Cu}$ ). (a) Three-dimensional mixed enthalpy $\Delta H^m-T-\sigma$ EAHNP chart with $\Delta H_e^m-T$ and $\Delta H_s^m-T$ paths; (b) $\Delta H_e^m-T$ and $\Delta H_s^m-T$ paths.

Figure 19. EAHNP and SAHNP charts on disordering AuCuI( $A_8^{Au}{A}_4^{Cu}$ ). (a) Three-dimensional mixed potential energy $\Delta E^m-T-\sigma$ EAHNP chart with $\Delta E_e^m-T$ and $\Delta E_s^m-T$ paths; (b) $\Delta E_e^m-T$ and $\Delta E_s^m-T$ paths.

Figure 20. EAHNP and SAHNP charts on disordering AuCuI( $A_8^{Au}{A}_4^{Cu}$ ). (a) Three-dimensional mixed volume $\Delta V^m-T-\sigma$ EAHNP chart with $\Delta V_e^m-T$ and $\Delta V_s^m-T$ paths; (b) $\Delta V_e^m-T$ and $\Delta V_s^m-T$ paths.

Figure 21. EAHNP and SAHNP charts on disordering AuCuI( $A_8^{Au}{A}_4^{Cu}$ ). (a) Three-dimensional mixed potential energy $\Delta U^m-T-\sigma$ EAHNP chart with $\Delta U_e^m-T$ and $\Delta U_s^m-T$ paths; (b) $\Delta U_e^m-T$ and $\Delta U_s^m-T$ paths.

Figure 22. EAHNP and SAHNP charts on disordering AuCuI( $A_8^{Au}{A}_4^{Cu}$ ). (a) Three-dimensional configurational entropy energy $T\Delta S^c-T-\sigma$ EAHNP chart with $T\Delta S_e^c-T$ and $T\Delta S_s^c-T$ paths; (b) $T\Delta S_e^c-T$ and $T\Delta S_s^c-T$ paths.

Figure 23. EAHNP and SAHNP charts on disordering AuCuI( $A_8^{Au}{A}_4^{Cu}$ ). (a) Three-dimensional generalized vibration entropy energy $T\Delta S^v-T-\sigma$ EAHNP chart with $T\Delta S_e^v-T$ and $T\Delta S_s^v-T$ paths; (b) $T\Delta S_e^v-T$ and $T\Delta S_s^v-T$ paths.

Figure 24. EAHNP and SAHNP charts on disordering AuCuI( $A_8^{Au}{A}_4^{Cu}$ ). (a) Three-dimensional generalized vibration free energy $\Delta X^v-T-\sigma$ EAHNP chart with $\Delta X_e^v-T$ and $\Delta X_s^v-T$ paths; (b) $\Delta X_e^v-T$ and $\Delta X_s^v-T$ paths.

Figure 25. Heat capacities on disordering AuCuI( $A_8^{Au}{A}_4^{Cu}$ ). (a) $c_p^e-T$ and $c_p^s-T$ paths on $c_p-x-T$ EAHNP phase diagram. (b) $c_p^e-T$ and $c_p^s-T$ paths. (These charts are calculated by temperature step $\Delta T=1\text{K}$ and order degree step $\Delta \sigma =0.00001$ ).

Figure 26. Thermal expansion coefficient on disordering AuCuI( $A_8^{Au}{A}_4^{Cu}$ ). (a) ${\alpha }^e-T$ and ${\alpha }^s-T$ paths on $\alpha -x-T$ EAHNP phase diagram. (b) ${\alpha }^e-T$ and ${\alpha }^s-T$ paths. (These charts are calculated by temperature step $\Delta T=1\text{K}$ and order degree step $\Delta \sigma =0.00001$ ).

Figure 27. Activities on disordering AuCuI( $A_8^{Au}{A}_4^{Cu}$ ). (a) $a_{Au}^e-T$ and $a_{Au}^s-T$ paths on $a_{Au}-x-T$ EAHNP phase diagram; (b) $a_{Cu}^e-T$ and $a_{Cu}^s-T$ paths on $a_{Cu}-x-T$ EAHNP phase diagram; (c) $a_{Au}^e-T$ and $a_{Au}^s-T$ paths; (d) $a_{Cu}^e-T$ and $a_{Cu}^s-T$ paths.

b) Choice of configurational entropy functions

In order to understand how to calculate the configurational entropy from the AG-probabilities $\left(x_i^{Au},x_i^{Cu}\right)$ occupied at the AG-Gibbs energy levels, we have calculated the $S^c\left(x,T,\sigma \right)$ EHNP-diagrams of the Au₃Cu-, AuCu-, AuCu₃-type sublattice systems by Equations ((4-32) to (4-37)), which are shown in Figures 34-36. In our opinion, the AG-probabilities occupied at the AG-Gibbs energy levels in the degeneracy g-function should be degenerated by the constituent atom probabilities occupied at the lattice points for calculating constituent atom arranging entropy, based on following reasons:

i) The AG generalized action entropy $S^v\left(x,T,\sigma \right)$ has been accounted in the $G^{*}\left(x,T,\sigma \right)$ -characteristic Gibbs energy (see Equations (4-16) to (4-27)). Therefore, the configurational entropy $S^c\left(x,T,\sigma \right)$ arisen from the AG-probabilities occupied at the lattice points should be degenerated by the constituent atom probabilities occupied at the lattice points.

Figure 28. EAHNP phase diagrams of sublattice systems. (a) Three-dimension $\Delta G_e^m-T-x$ EAHNP diagrams of Au₃Cu-type sublattice system; (b) $T-x$ equilibrium phase diagram with iso-mixed Gibbs energy $T_{\Delta G^m}-x$ curves of Au₃Cu-type sublattice system; (c) Three-dimension $\Delta G_e^m-T-x$ EAHNP diagrams of AuCu-type sublattice system; (d) $T-x$ equilibrium phase diagram with iso-mixed Gibbs energy $T_{\Delta G^m}-x$ curves of AuCu-type sublattice system; (e) Three-dimension $\Delta G_e^m-T-x$ EAHNP diagrams of AuCu₃-type sublattice system; (f) $T-x$ equilibrium phase diagram with iso-mixed Gibbs energy $T_{\Delta G^m}-x$ curves of AuCu₃-type sublattice system.

Figure 29. EAHNP phase diagrams of sublattice systems. (a) Three-dimension ${\sigma }_e-x-T$ EAHNP diagrams of Au₃Cu-type sublattice system; (b) $T-x$ equilibrium phase diagram with iso-order degree $T_{\sigma }-x$ curves, top-points and experimental jumping $T_j$ -temperatures of Au₃Cu-type sublattice system; (c) Three-dimension ${\sigma }_e-x-T$ EAHNP diagrams of AuCu-type sublattice system; (d) $T-x$ equilibrium phase diagram with iso-order degree $T_{\sigma }-x$ curves, top-points and experimental jumping $T_j$ - temperatures of AuCu-type sublattice system; (e) Three-dimension ${\sigma }_e-x-T$ EAHNP diagrams of AuCu₃-type sublattice system; (f) $T-x$ equilibrium phase diagram with iso-order degree $T_{\sigma }-x$ curves, top-points and experimental jumping $T_j$ -tempera- tures of AuCu₃-type sublattice system.

Figure 30. EAHNP phase diagrams. (a) (b) Three-dimension $x_{i,e}^{Au}-x-T$ and $x_{i,e}^{Cu}-x-T$ diagrams of Au₃Cu-type alloy; (c), (d) Three-dimension $x_{i,e}^{Au}-x-T$ and $x_{i,e}^{Cu}-x-T$ diagrams of AuCu-type alloy; (e) (f) Three-dimension $x_{i,e}^{Au}-x-T$ and $x_{i,e}^{Cu}-x-T$ diagrams of AuCu₃-type alloy.

Figure 31. Equilibrium EAHNP phase diagrams of Au-Cu system. (a) Equilibrium EAHNP phase diagram with iso-mixed Gibbs energy $T_{\Delta G^m}-x$ curves of Au-Cu system; (b) Equilibrium EAHNP phase diagram with iso-order degree $T_{\sigma }-x$ curves of Au-Cu system. (These diagrams are calculated by temperature step $\Delta T=1\text{K}$ , composition step $\Delta x=0.5\%$ and ordered degree step $\Delta \sigma =0.001$ ).

Figure 32. Subequilibrium SAHNP phase diagrams of Au-Cu system. (a) Subequilibrium SAHNP phase diagram with iso-mixed Gibbs energy $T_{\Delta G^m}-x$ curves of Au-Cu system; (b) Subequilibrium SAHNP phase diagram with iso-order degree $T_{\sigma }-x$ curves of Au-Cu system. (These diagrams are calculated by temperature step $\Delta T=1\text{K}$ , composition step $\Delta x=0.5\%$ and ordered degree step $\Delta \sigma =0.001$ .)

(a)(b)

Figure 33. Ranges of long range ordered phases. (a) Composition ranges of ordered Au₃Cu-type, AuCu-type and AuCu₃-type in Au-Cu system calculated by their first jumping order degrees of Au₃Cu- and AuCu₃-type alloys and subequilibrium crosspint numerical method of isothermal Gibbs energy curves [47] [48] , together with experimental jumping points, which were erroneously considered as the equilibrium order-disorder transformation critical $T_c$ -temperatures from references given in [88] ; (b) Comparison of composition ranges of long range ordered phases. The alloys were annealed at appropriate temperature different periods (on six months and one year) and slowly cooled down to room temperature [60] [61] . For alloys on six months, the experimental composition range is B-D-G-H; For alloys on one year, the experimental composition range is A-E-F-H, which approachs to calculated limited composition range J₁-O₁-O₂-J₄. Each point is described by composition, mixed Gibbs energy and order degree: J₁ (18.875, −4220.47, 0.7540), A (21.00, −4807.10, 0.8400), C₁ (25.00, −5955.72, 1), B (27.50, −6094.93, 0.9667), D (35.50, −6411.03, 0.8600), O₁ (36.20, −6428.97, 0.8507), E (36.80, −6443.06, 0.8427), J₂ (43.375, −6518.13, 0.7553), L₁ (43.40, −6518.13, 0.7547), C₂ (50, −8746.23, 1), J₃ (56.625, 0.7547), F (62.50, −7094.10, 0.8333), O₂ (63.05, −7108.69, 0.8407), G (66.50, −7176.61, 0.8667), L₂ (70.35, −7203.16, 0.9380), C₃ (75, −7163.61, 1), H (79.50, −5726.05, 0.8200), J₄ (81.125, −5320.59, 0.7550).

Figure 34. Congfiguration entropies of the Au₃Cu-type sublattice system: (a) calculated by $S^c{\left(x,T,\sigma \right)}^{\left(1\right)}$ ; (b) calculated by $S^c{\left(x,T,\sigma \right)}^{\left(2\right)}$ ; (c) calculated by $S^c{\left(x,T,\sigma \right)}^{\left(3\right)}$ (the maximal order degree ${\sigma }_{\max }$ is the function of $x_{Cu}$ ).

Figure 35. Congfiguration entropies of the AuCu-type sublattice system: (a) calculated by $S^c{\left(x,T,\sigma \right)}^{\left(1\right)}$ ; (b) calculated by $S^c{\left(x,T,\sigma \right)}^{\left(2\right)}$ ; (c) calculated by $S^c{\left(x,T,\sigma \right)}^{\left(3\right)}$ (the maximal order degree ${\sigma }_{\max }$ is the function of $x_{Cu}$ ).

Figure 36. Congfiguration entropies of the AuCu₃-type sublattice system: (a) calculated by $S^c{\left(x,T,\sigma \right)}^{\left(1\right)}$ ; (b) calculated by $S^c{\left(x,T,\sigma \right)}^{\left(2\right)}$ ; (c) calculated by $S^c{\left(x,T,\sigma \right)}^{\left(3\right)}$ (the maximal order degree ${\sigma }_{\max }$ is the function of $x_{Cu}$ ).

ii) The configurational entropy of each ordered alloy calculated by a degeneracy function can change continually from the configurational entropy in the maximum order degree state to one in the perfectly disordered state.

iii) For binary alloy systems, the $S^c\left(x,T,\sigma \right)$ should lie in the range:

$0\le S^c\left(x,T,\sigma \right)\le \text{5}.\text{763}\left(\text{J}\cdot {\text{mol}}^{-\text{1}}\cdot {\text{K}}^{-\text{1}}\right)$

where, the upper limit is the value for a perfectly disordered equiatomic Au_0.5Cu_0.5 alloy.

The results show that the AG arranging degeneracy function $g{\left(x,T,\sigma \right)}^{\left(1\right)}$ and entropy function $S^c{\left(x,T,\sigma \right)}^{\left(1\right)}$ may be used [48] .

The AG arranging degeneracy function $g{\left(x,T,\sigma \right)}^{\left(1\right)}$ and entropy function $S^c{\left(x,T,\sigma \right)}^{\left(1\right)}$ :

$\begin{array}{c}g{\left(x,T,\sigma \right)}^{\left(1\right)}=\frac{\left[N{\displaystyle \underset{i=0}{\overset{I}{\sum }}\left(x_i^{Au\left(1\right)}\left(x,T,\sigma \right)+x_i^{Cu\left(1\right)}\right)\left(x,T,\sigma \right)}\right]!}{\left({\displaystyle \underset{i=0}{\overset{I}{\sum }}N{x}_i^{Au\left(1\right)}\left(x,T,\sigma \right)}\right)!\left({\displaystyle \underset{i=0}{\overset{I}{\sum }}N{x}_i^{Cu\left(1\right)}\left(x,T,\sigma \right)}\right)!}\\ \times \frac{\left[N{\displaystyle \underset{i=0}{\overset{I}{\sum }}\left(x_i^{Au\left(2\right)}\left(x,T,\sigma \right)+x_i^{Cu\left(2\right)}\right)\left(x,T,\sigma \right)}\right]!}{\left({\displaystyle \underset{i=0}{\overset{I}{\sum }}N{x}_i^{Au\left(2\right)}(x,T,\sigma )}\right)!\left({\displaystyle \underset{i=0}{\overset{I}{\sum }}N{x}_i^{Cu\left(2\right)}\left(x,T,\sigma \right)}\right)!}\end{array}$ (4-32)

$\begin{array}{c}{S}^c{\left(x,T,\sigma \right)}^{\left(1\right)}=k_B\ln \frac{\left[N{\displaystyle \underset{i=0}{\overset{I}{\sum }}\left(x_i^{Au\left(1\right)}\left(x,T,\sigma \right)+x_i^{Cu\left(1\right)}\left(x,T,\sigma \right)\right)}\right]!}{\left({\displaystyle \underset{i=0}{\overset{I}{\sum }}N{x}_i^{Au\left(1\right)}\left(x,T,\sigma \right)}\right)!\left({\displaystyle \underset{i=0}{\overset{I}{\sum }}N{x}_i^{Cu\left(1\right)}\left(x,T,\sigma \right)}\right)!}\\ \times \frac{\left[N{\displaystyle \underset{i=0}{\overset{I}{\sum }}\left(x_i^{Au\left(2\right)}\left(x,T,\sigma \right)+x_i^{Cu\left(2\right)}\left(x,T,\sigma \right)\right)}\right]!}{\left({\displaystyle \underset{i=0}{\overset{I}{\sum }}N{x}_i^{Au\left(2\right)}\left(x,T,\sigma \right)}\right)!\left({\displaystyle \underset{i=0}{\overset{I}{\sum }}N{x}_i^{Cu\left(2\right)}\left(x,T,\sigma \right)}\right)!}\end{array}$ (4-33)

The AG arranging degeneracy function $g{\left(x,T,\sigma \right)}^{\left(2\right)}$ and entropy function $S^c{\left(x,T,\sigma \right)}^{\left(2\right)}$ :

$\begin{array}{c}g{\left(x,T,\sigma \right)}^{\left(2\right)}=\frac{\left[N{\displaystyle \underset{i=0}{\overset{I}{\sum }}{x}_i^{Au\left(1\right)}\left(x,T,\sigma \right)}\right]!\left[{\displaystyle \underset{i=0}{\overset{I}{\sum }}N{x}_i^{Cu\left(1\right)}\left(x,T,\sigma \right)}\right]!}{\left({\displaystyle \underset{i=0}{\overset{I}{\prod }}N{x}_i^{Au\left(1\right)}\left(x,T,\sigma \right)}\right)!\left({\displaystyle \underset{i=0}{\overset{I}{\prod }}N{x}_i^{Cu\left(1\right)}\left(x,T,\sigma \right)}\right)!}\\ \times \frac{\left[N{\displaystyle \underset{i=0}{\overset{I}{\sum }}{x}_i^{Au\left(2\right)}\left(x,T,\sigma \right)}\right]!\left[{\displaystyle \underset{i=0}{\overset{I}{\sum }}N{x}_i^{Cu\left(2\right)}\left(x,T,\sigma \right)}\right]!}{\left({\displaystyle \underset{i=0}{\overset{I}{\prod }}N{x}_i^{Au\left(2\right)}\left(x,T,\sigma \right)}\right)!\left({\displaystyle \underset{i=0}{\overset{I}{\prod }}N{x}_i^{Cu\left(2\right)}\left(x,T,\sigma \right)}\right)!}\end{array}$ (4-34)

$\begin{array}{c}{S}^c{\left(x,T,\sigma \right)}^{\left(2\right)}=k_B\ln \frac{\left[N{\displaystyle \underset{i=0}{\overset{I}{\sum }}{x}_i^{Au\left(1\right)}\left(x,T,\sigma \right)}\right]!\left[{\displaystyle \underset{i=0}{\overset{I}{\sum }}N{x}_i^{Cu\left(1\right)}\left(x,T,\sigma \right)}\right]!}{\left({\displaystyle \underset{i=0}{\overset{I}{\prod }}N{x}_i^{Au\left(1\right)}\left(x,T,\sigma \right)}\right)!\left({\displaystyle \underset{i=0}{\overset{I}{\prod }}N{x}_i^{Cu\left(1\right)}\left(x,T,\sigma \right)}\right)!}\\ \times \frac{\left[N{\displaystyle \underset{i=0}{\overset{I}{\sum }}{x}_i^{Au\left(2\right)}\left(x,T,\sigma \right)}\right]!\left[{\displaystyle \underset{i=0}{\overset{I}{\sum }}N{x}_i^{Cu\left(2\right)}\left(x,T,\sigma \right)}\right]!}{\left({\displaystyle \underset{i=0}{\overset{I}{\prod }}N{x}_i^{Au\left(2\right)}\left(x,T,\sigma \right)}\right)!\left({\displaystyle \underset{i=0}{\overset{I}{\prod }}N{x}_i^{Cu\left(2\right)}\left(x,T,\sigma \right)}\right)!}\end{array}$ (4-35)

The AG arranging degeneracy function $g{\left(x,T,\sigma \right)}^{\left(3\right)}$ and entropy function $S^c{\left(x,T,\sigma \right)}^{\left(3\right)}$ :

$\begin{array}{c}g{\left(x,T,\sigma \right)}^{\left(3\right)}=\frac{\left[N{\displaystyle \underset{i=0}{\overset{I}{\sum }}\left(x_i^{Au\left(1\right)}\left(x,T,\sigma \right)+x_i^{Cu\left(1\right)}\left(x,T,\sigma \right)\right)}\right]!}{\left({\displaystyle \underset{i=0}{\overset{I}{\prod }}N{x}_i^{Au\left(1\right)}\left(x,T,\sigma \right)}\right)!\left({\displaystyle \underset{i=0}{\overset{I}{\prod }}N{x}_i^{Cu\left(1\right)}\left(x,T,\sigma \right)}\right)!}\\ \times \frac{\left[N{\displaystyle \underset{i=0}{\overset{I}{\sum }}\left(x_i^{Au\left(2\right)}\left(x,T,\sigma \right)+x_i^{Cu\left(2\right)}\left(x,T,\sigma \right)\right)}\right]!}{\left({\displaystyle \underset{i=0}{\overset{I}{\prod }}N{x}_i^{Au\left(2\right)}\left(x,T,\sigma \right)}\right)!\left({\displaystyle \underset{i=0}{\overset{I}{\prod }}N{x}_i^{Cu\left(2\right)}\left(x,T,\sigma \right)}\right)!}\end{array}$ (4-36)

$\begin{array}{c}{S}^c{\left(x,T,\sigma \right)}^{\left(3\right)}=k_B\ln \frac{\left[N{\displaystyle \underset{i=0}{\overset{I}{\sum }}\left(x_i^{Au\left(1\right)}\left(x,T,\sigma \right)+x_i^{Cu\left(1\right)}\left(x,T,\sigma \right)\right)}\right]!}{\left({\displaystyle \underset{i=0}{\overset{I}{\prod }}N{x}_i^{Au\left(1\right)}\left(x,T,\sigma \right)}\right)!\left({\displaystyle \underset{i=0}{\overset{I}{\prod }}N{x}_i^{Cu\left(1\right)}\left(x,T,\sigma \right)}\right)!}\\ \times \frac{\left[N{\displaystyle \underset{i=0}{\overset{I}{\sum }}\left(x_i^{Au\left(2\right)}\left(x,T,\sigma \right)+x_i^{Cu\left(2\right)}\left(x,T,\sigma \right)\right)}\right]!}{\left({\displaystyle \underset{i=0}{\overset{I}{\prod }}N{x}_i^{Au\left(2\right)}\left(x,T,\sigma \right)}\right)!\left({\displaystyle \underset{i=0}{\overset{I}{\prod }}N{x}_i^{Cu\left(2\right)}\left(x,T,\sigma \right)}\right)!}\end{array}$ (4-37)

The AHNP-social road is to turn the spiritual wealth of SMM science into the material wealth of advanced alloys. Therefore, the AGSP should add other three goals: We have to “surmount super era theory trap”, before we can start the AGSP. Because that the philosophic thinking mode and theoretical model of the SMM science are different from ones of QMAC electron theory and CALPHAD thermodynamics, specially due to social/interest factor, the SMM science has been inhibited; Achieving shift in social form to the AHNP-social organization, that drives deeper collaboration, vaster open, fuller utilization and fairer share, then provides social relations support for the AGSP; Achieving transformation of the human thinking mode from “reductionism” to the SSP, that drives well-qualified scientists, technicians and management personnel with SSP thinking mode, then provides human resources support for the AGSP.

Acknowledgements

This research was supported by the National Natural Science Foundation of China (51071181), the Natural Science Foundation of Hunan Province (2013FJ4043), Development and Reform Commission of Hunan Province (No.2016-01), China.

Conflicts of Interest

The authors declare no conflicts of interest.

References

[1]	Slater, J.C. (1937) Electronic Structure of Alloys. Journal of Applied Physics, 8, 385-390. https://doi.org/10.1063/1.1710311
[2]	Hume-Rothery, W. (1962) Atomic Theory for Students of Metallurgy. Institute of metals, London.
[3]	Hohenberg, P. and Kohn, W. (1964) Inhomogeneous Electron Gas. Physical Review B, 136, 864-871. https://doi.org/10.1103/PhysRev.136.B864
[4]	Kohn, W. and Sham, L. (1965) Self-Consistent Equations Including Exchange and Correlation Effects. Physical Review A, 140, 1133-1138. https://doi.org/10.1103/PhysRev.140.A1133
[5]	Kikuchi, R.A. (1951) Theory of Cooperative Phenomemena. Physical Review, 81, 988-1003. https://doi.org/10.1103/PhysRev.81.988
[6]	Asta, M., Wolverton, C., de Fontaine, D. and Dreyssé, H. (1991) Effective Interaction from Cluster-Variation Formalism I. Physical Review B, 44, 4907-4913. https://doi.org/10.1103/PhysRevB.44.4907
[7]	Wolverton, C., Asta, M., Dreyssé, H. and de Fontaine, D. (1991) Effective Interaction from Cluster-Variation Formalism II. Physical Review B 44, 4914-4924. https://doi.org/10.1103/PhysRevB.44.4914
[8]	Asta, M. and de Fontaine D. (1993) First-Principles Study of Phase Stability of Ti-Al Intermetallic Compounds. Journal of Materials Research, 8, 2554-2568. https://doi.org/10.1557/JMR.1993.2554
[9]	Ellis, D.E. and Perdwew, J.P. (1994) Density Function Theory of Molecules, Clusters and Solids. Kluwer Academic, Holland.
[10]	Pettifor, D.G. (2003) Electron Theory in Materials Modeling. Acta Materialia, 51, 5649-5673. https://doi.org/10.1016/S1359-6454(03)00466-X
[11]	Bragg, W.L. and Williams E.J. (1934) The Effect of Thermal Agitation on Atomic Arrangement in Alloys. Proceedings of the Royal Society of London A, 145, 699-730. https://doi.org/10.1098/rspa.1934.0132
[12]	Bragg, W.L. and Williams E.J. (1935) The Effect of Thermal Agitation on Atomic Arrangement in Alloys-II. Proceedings of the Royal Society of London A, 151, 540-566. https://doi.org/10.1098/rspa.1935.0165
[13]	Bethe, H.A. (1935) Statistical Theory of Superlattices. Proceedings of the Royal Society of London A, 150, 552-575. https://doi.org/10.1098/rspa.1935.0122
[14]	Kaufman L. and Bernstein H. (1970) Computer Calculation of Phase Diagrams. Academic Press, New York.
[15]	Sundman, B., Fries, S.G. and Oates, W.A. (1999) A Calphad Assessment of the Au-Cu System Using the Cluster Variatiation Method. Zeitschrift für Metallkunde, 90, 267-273.
[16]	Hillert, M. (2001) The Compound Energy Formalism. Journal of Alloys and Compounds, 320, 161-176. https://doi.org/10.1016/S0925-8388(00)01481-X
[17]	Oates, W.A. (2006) Configurational Entropies of Mixing in Solid Alloys. Journal of Phase Equilibria and Diffusion, 28, 79-89. https://doi.org/10.1007/s11669-006-9008-3
[18]	Kaufman, L. and Agren, J. (2014) CALPHAD, First and Second Generation-Birth of the Materials Genome. Scripta Materialia, 70, 3-6. https://doi.org/10.1016/j.scriptamat.2012.12.003
[19]	Xie, Y.Q. (1998) Systematic Science of Metallic Materials. Central South University of Technology Press, Changsha. (In Chinese)
[20]	Xie, Y.Q. (1998) Development Course of Metallic Materials Science and Evolution of Human Thinking Mode. Materials Review, 12, 6-10. (In Chinese)
[21]	Xie, Y.Q. (2011) Relationship between Partial and Average Atomic Volumes of Components in Au-Ni Alloys. Transactions of Nonferrous Metals Society of China, 21, 1801-1807. https://doi.org/10.1016/j.scriptamat.2012.12.003
[22]	Xie, Y.Q. (2012) Systematic Science of Alloys. Central South University Press, Changsha. (In Chinese)
[23]	Xie, Y.Q., Liu X.B., Li, X.B., Peng, H.J. and Nie, Y.Z. (2015) Holographic Alloy Positioning Design System and Holographic Network Phase Diagrams of Au-Cu System. Transactions of Nonferrous Metals Society of China, 25, 885-906. https://doi.org/10.1016/S1003-6326(15)63678-0
[24]	Xie, Y.Q., Li X.B., Liu, X.B., Peng, H.J. and Nie, Y.Z. (2015) Alloy Gene, Holographic Alloy Positioning Design System. The Chinese Journal of Nonferrous Metals, 25, 2798-2843.
[25]	Pauling, L. (1960) The Nature of the Chemical Bonding. Conell University Press, Ithaca.
[26]	Yu, R.H. (1978) Empirical Electron Theory of Solids and Molecules. Science Bulletin, 23, 217-224.
[27]	Grigorovvich, V.K. (1989) The Metallic Bond and Structure of Metals. Nova Science, New York.
[28]	Lupis, C.H.P. and Elliott, J.F. (1967) Prediction of Enthalpy and Entropy Interaction Coefficients by the “Central Atoms” Theory. Acta Materialia, 15, 265-276. https://doi.org/10.1016/0001-6160(67)90202-7
[29]	Xie, Y.Q. and Ma, L.Y. (1985) Theoretical Lattice Parameters of Valence Electron Structures of Crystals. Journal of Central South Institute of Mining and Metallurgy, 1, 1-8. (In Chinese)
[30]	Xie, Y.Q. and Ma, L.Y. (1985) Generalized Vegard’s Law and Generalized Yü’s Law. Journal of Central South Institute of Mining and Metallurgy, 3, 1-7. (In Chinese)
[31]	Xie, Y.Q. (1993) A New Potential Function with Many-Atom Interactions in Solid. Science in China (Series E), 36, 90-99.
[32]	Xie, Y.Q., Zhang, X.D., Zhao, L.Y. and Ma, X.L. (1993) Electronic Structure and Properties of Cu Metal. Science in China (Series E), 36, 487-496.
[33]	Xie, Y.Q. (1994) Electronic Structure and Properties of Pure Iron. Acta Metallurgica et Materialia, 42, 3705-3715. https://doi.org/10.1016/0956-7151(94)90436-7
[34]	Xie, Y.Q., Zhang, X.D. and Chen, J.Y. (1996) Electronic Structure and Properties of Pure Cobalt. Science in China (Series E), 39, 394-403.
[35]	Xie, Y.Q., Lü, M.S. and Zhao, L.Y. (1993) Electronic Structure and Properties of Ni Metal. Science in China (Series A), 36, 495-503.
[36]	Xie, Y.Q., Zhang, X.D., Zou, P. and Zhao, L.Y. (1993) Microstructure and Properties of Cu-Ni Alloys. Science in China, A, 36, 612-623.
[37]	Xie, Y.Q. (1998) Atomic Energies and Gibbs Energy Functions of Ag-Cu Alloys. Science in China, Series E: Technological Sciences, 41, 146-156. https://doi.org/10.1007/BF02919677
[38]	Xie, Y.Q. and Zhang, X.D. (1998) Atomic Volume and Volume Functions of Ag-Cu Alloys. Science in China, Series E: Technological Sciences, 41, 57-168. https://doi.org/10.1007/BF02919678
[39]	Xie, Y.Q. and Zhang, X.D. (1998) Electronic Structure of Ag-Cu Alloys. Science in China, Series E: Technological Sciences, 41, 225-236. https://doi.org/10.1007/BF02919432
[40]	Xie, Y.Q. and Zhang, X.D. (1998) Phase Diagram and Thermodynamic Properties of Ag-Cu Alloys. Science in China, Series E: Technological Sciences, 41, 348-356. https://doi.org/10.1007/BF02917006
[41]	Xie, Y.Q., Peng, K. and Liu, X.B. (2004) Influences of xTi/xAl Atomic States, Lattice Constants and Potential-Energy Planes of Ordered FCC TiAl Type Alloys. Physica B, 344, 5-20. https://doi.org/10.1016/j.physb.2003.06.002
[42]	Xie, Y.Q., Liu, X.B. and Peng, K. (2004) Atomic State, Potential Energies, Volumes, Stability and Brittleness of Ordered FCC TiAl3 Type Alloys. Physica B, 353, 15-33. https://doi.org/10.1016/j.physb.2004.08.022
[43]	Xie, Y.Q., Peng, H.J., Liu, X.B. and Peng, K. (2005) Potential Energies, Volumes, Stability and Brittleness of Ordered FCC Ti3Al Type Alloys. Physica B, 362, 1-17. https://doi.org/10.1016/j.physb.2005.01.197
[44]	Xie, Y.Q., Tao, H.J., Peng, H.J., Li, X.B., Liu, X.B. and Peng, K. (2005) Atomic State, Potential Energies, Volumes, Stability and Brittleness of Ordered FCC TiAl2 Type Alloys. Physica B, 366, 17-37. https://doi.org/10.1016/j.physb.2005.04.041
[45]	Xie, Y.Q., Li, X.B., Liu, X.B., Nie, Y.Z. and Peng, H.J. (2013) Alloy Gene Gibbs Energy Partition Function and Equilibrium Holographic Network Phase Diagrams of AuCu-Type Sublattice System. International Journal of Communications, Network and System Sciences, 6, 415-442. https://doi.org/10.4236/ijcns.2013.610045
[46]	Xie, Y.Q., Peng, H.J., Liu, X.B., Li, X.B. and Nie, Y.Z. (2014) New Atom Movement Mechanism for Tracking Path on Disordering AuCuI(A8AuA4Cu) Compound. Transactions of Nonferrous Metals Society of China , 24, 3221-3256. https://doi.org/10.1016/S1003-6326(14)63464-6
[47]	Xie, Y.Q., Li, X.B., Liu, X.B., Nie, Y.Z. and Peng, H.J. (2014) Alloy Gene Gibbs Energy Partition Function and Equilibrium Holographic Network Phase Diagrams of AuCu3-Type Sublattice System. Transactions of Nonferrous Metals Society of China, 24, 3585-3610. https://doi.org/10.1016/S1003-6326(14)63505-6
[48]	Xie, Y.Q., Nie, Y.Z., Li, X.B., Peng, H.J. and Liu, X.B. (2015) Alloy Gene Gibbs Energy Partition Function and Equilibrium Holographic Network Phase Diagrams of Au3Cu-Type Sublattice System. Transactions of Nonferrous Metals Society of China, 25, 211-240. https://doi.org/10.1016/S1003-6326(15)63598-1
[49]	OSTP (2011) Materials Genome Intiative Strantegic for Global Competiveness.: Office of Science and Technology Policy, OSTP, Washington DC. OSTP (2014) Materials Genome Intiative Strantegic Plan. Office of Science and Technology Policy, Washington DC.
[50]	Zhao, J.C. (2014) A Perspective on Materials Genome Initiative. Chinese Journal of Nature, 36, 89-104.
[51]	Johannsson, C.H. and Linde, J.O. (1936) Roentgenographic and Electrical Investigations of the Cu-Au System. Annals of Physics, 25, 1-48.
[52]	Nix, F.C. and Macnair, D. (1941) A Dilatometric Study of the Order-Disorder Transformation in Cu-Au Alloys. Physical Review, 60, 320-329. https://doi.org/10.1103/PhysRev.60.320
[53]	Sutcliffe, C.H. and Jaumot Jr., F.E. (1953) Order-Disorder in Au-Cu Alloy. I. Short-Range Order in an Alloy Containing Atomic Percent Au. Acta Metallurgica, 1, 725-730. https://doi.org/10.1016/0001-6160(53)90032-7
[54]	Frank, E., Jaumot, J.R. and Charles, H.S. (1954) Order-Disorder in Au-Cu alloy II. The Nature of the Order-Disorder Transformation and Long-Range Order. Acta Metallurgica, 2, 63-74. https://doi.org/10.1016/0001-6160(54)90095-4
[55]	Kuczynski, G.C., Doyama, M. and Fine, M.E. (1956) Transformations in Disordered Gold Copper Alloys. Journal of Applied Physics, 27, 651-655. https://doi.org/10.1063/1.1722450
[56]	O’brien, J.L. and Kuczynski, G.C. (1959) X-Ray Study of the Kinetics of Ordering in AuCu. Acta Metallurgica., 7, 803-806. https://doi.org/10.1016/0001-6160(59)90096-3
[57]	Scott, R.E. (1960) New Complex Phase in the Copper-Gold System. Journal of Applied Physics, 31, 2112-2117. https://doi.org/10.1063/1.1735509
[58]	Orr, R.L., Luciat-Labry, J. and Hultgren, R. (1960) Energy of the Order-Disorder Transformation in AuCu. Acta Metallurgica, 8, 431-434. https://doi.org/10.1016/0001-6160(60)90029-8
[59]	Marcinkowski, M.J. and Zwell, L. (1963) Transmission Electron Microcopy Study of the Off-Staichiometric Cu3Au Superlattices. Acta Metallurgica, 11, 373-390. https://doi.org/10.1016/0001-6160(63)90162-7
[60]	Lu, S.S. and Liang, C.K. (1966) The Superlattice Formation and Lattice Spacing Changes in Copper-Gold Alloys. Acta Physica Sincica, 22, 639-696 (in Chinese).
[61]	Lu, S.S. and Liang, C.K. (1966) The Existence of the Superlattice CuAu3 in the Cu-Au System. Kexue Tongbao, 17, 395-396.
[62]	Marques, M.C., Pinheiro, M.F. and Braga, M.E. (1977) Electrical Resistivity and Phase Change in CuAu. Scripta Metallurgica, 11, 197-200. https://doi.org/10.1016/0036-9748(77)90052-7
[63]	Claeson, T. and Boyce, J.B. (1984) Order-Disorder Transformation in Au-Cu Alloys Studied by Extended X-Ray-Absorption Fine Structure. Physical Review, 29, 1551-1557. https://doi.org/10.1103/PhysRevB.29.1551
[64]	Lander, G.H. and Brown, P.J. (1985) A Neutron Diffraction Study of AuCu3 near the Order-Disorder Temperature. Journal of Physics C: Solid State Physics, 18, 2017-2024. https://doi.org/10.1088/0022-3719/18/10/007
[65]	Butler, B.D. and Cohen, J.B. (1989) The Structure of Cu3Au above the Critical Temperature. Journal of Applied Physics, 65, 2214-2219. https://doi.org/10.1063/1.342833
[66]	Potekaev, A.I. (1995) Long-Period States of Ordered Metal Alloys. 1. Analysis of Structural Features. Russian Physics Journal., 38, 549-562. https://doi.org/10.1007/BF00559919
[67]	Potekaev, A.I. (1996) Long-Period States of Ordered Metal Alloys. 2. Physical Properties. Russian Physics Journal, 39, 521-533. https://doi.org/10.1007/BF02437017
[68]	Sprusil, B. and Pfeiler, W. (1997) The Retro-Effect in Stoichiometric CuAu: A Resistometric Study. Intermetallics, 5, 501-505. https://doi.org/10.1016/S0966-9795(97)00019-8
[69]	Bogdanoff, P.D. and Fultz, B. (1999) Vibrational Entropy of L12 Cu3Au Measured by Inelastic Neutron Scattering. Physical Review B, 36, 3976-3981. https://doi.org/10.1103/PhysRevB.60.3976
[70]	Sprusil, B. and Pfeiler, W. (2000) The Retro-Effect in CuAu and Its Kinetic Limitations. Intermetallics, 8, 81-83. https://doi.org/10.1016/S0966-9795(99)00076-X
[71]	Pfeiler, W. and Sprusil, B. (2002) Atomic Ordering in Alloys: Stable States and Kinetics. Materials Science and Engineering A, 324, 34-42. https://doi.org/10.1016/S0921-5093(01)01280-1
[72]	Kubiak, R. and Janczak, J. (1991) X-Ray Study of Ordered Phase Formation in Au31.6Cu68.4, Au50Cu50 and Au75Cu25. Journal of Alloys and Compounds, 176, 133-140. https://doi.org/10.1016/0925-8388(91)90018-Q
[73]	Bonneaux, J. and Guymont, M. (1999) Study of the Order-Disorder Transition Series in AuCu by In-Situ Temperature Electron Microscopy. Intermetallics, 7, 797-805. https://doi.org/10.1016/S0966-9795(98)00128-9
[74]	Feutelais, Y., Legendre, B. and Guymont, M. (1999) New Enthalpies Determination and In Situ X-Ray Diffraction Observations of Order/Disorder Transitions in Au0.5Cu0.5. Acta Materialia, 47, 2539-2551. https://doi.org/10.1016/S1359-6454(99)00089-0
[75]	Volker, E., Williams, F.J., Calvo, E.J., Jacob, T. and Schiffrin, D.J. (2012) O2 Induced Cu Surface Segregation in Au-Cu Alloys Studied by Angle Resolved XPS and DFT Modelling. Physical Chemistry Chemical Physics, 14, 7448-7455. https://doi.org/10.1039/c2cp40565b
[76]	Miranda, G., Silva, F.S. and Soares, D. (2013) Solid State Transformations and Equilibrium Crystal Structures of an Au-Cu Alloy with Shape Memory Effect. Materials Science Forum, 732, 859-864.
[77]	Wu, J., Zou, L., Wang, T., Chen, Y. and Cai, J. (2016) Spin Hall Angle and Spin Diffusion Length in Au-Cu Alloy. IEEE Transactions on Magnetics, 52, 1-4. https://doi.org/10.1109/TMAG.2016.2522938
[78]	Xie, Y.Q., Li, X.B., Liu, X.B., Peng, H.J. and Nie, Y.Z. (2009) Potential Energy Sequences of Characteristic Atoms on Basis of Heats of Formation of Disordered Au(1−x)Cux Alloys (Part 1). Journal of Materials Science and Engineering, 3, 51-68.
[79]	Xie, Y.Q., Li, X.B., Liu, X.B., Peng, H.J. and Nie, Y.Z. (2009) Potential Energy Sequences of Characteristic Atoms on Basis of Heats of Formation of Disordered Au(1−x)Cux Alloys (Part 2). Journal of Materials Science and Engineering, 3, 44-57.
[80]	Xie, Y.Q., Li, X.B., Liu, X.B., Peng, H.J. and Nie, Y.Z. (2009) Volume Sequences of Characteristic Atoms Separated from Experimental Volumes of Disordered Au(1−x)Cux Alloys. Journal of Materials Science and Engineering, 3, 54-75.
[81]	Xie, Y.Q., Liu, X.B., Li, X.B., Peng, H.J. and Nie, Y.Z. (2009) Potential Energies of Characteristic Atoms Separated from Experimental Heats of Formation of AuCu and AuCu3 Compounds. Transactions of Nonferrous Metals Society of China, 19, 1243-1256. https://doi.org/10.1016/S1003-6326(08)60436-7
[82]	Xie, Y.Q., Liu, X.B., Li, X.B., Peng, H.J. and Nie, Y.Z. (2009) Volume of Characteristic Atoms Separated from Experimental Volumes of AuCu and AuCu3 Compounds. Transactions of Nonferrous Metals Society of China, 19, 1599-1617. https://doi.org/10.1016/S1003-6326(09)60077-7
[83]	Xie, Y.Q., Nie, Y.Z., Liu, X.B., Li, X.B., Peng, H.J. and Li, Y.F. (2012) Potential Energies of Characteristic Atoms Separated from First-Principles Calculated Heats of Formation of AuCu and AuCu3 Compounds. In: Xie, Y.Q., Ed., Systematic Science of Alloys, Central South University Press, Changsha, 405-427.
[84]	Xie, Y.Q., Nie, Y.Z., Liu, X.B., Li, X.B., Peng, H.J. and Li, Y.F. (2012) Volume of Characteristic Atoms Separated from First-Principles Calculated Volume of L10-AuCu and L12-AuCu3 Compounds. S In: Xie, Y.Q., Ed., Systematic Science of Alloys, Central South University Press, Changsha, 428-460.
[85]	Xie, Y.Q., Nie, Y.Z., Li, X.B., Liu, X.B., Peng, H.J. and Li, Y.F. (2011) Characteristic Atoms Potential Energy Partition Function of Au3Cu Type Ordered Alloys in FP-SSA Framework. The Chinese Journal of Nonferrous Metals, 21, 2489-2501.
[86]	Xie, Y.Q., Li, Y.F., Liu, X.B., Li, X.B., Peng, H.J. and Nie, Y.Z. (2011) Characteristic Atoms Occupation Patterns of Au3Cu, AuCu3, AuCuI and AuCuII Based on Experimental Data of Disordered Alloys. Transactions of Nonferrous Metals Society of China, 21, 1092-1104. https://doi.org/10.1016/S1003-6326(11)60827-3
[87]	Mie, G. (1903) Zur kinetischen theorie der einatomigen korper. Annals of Physics 316, 657-697. https://doi.org/10.1002/andp.19033160802
[88]	Lennard-Jones, J.E. (1931) Cohesion. Proceedings of the Physical Society, 43, 461-482. https://doi.org/10.1088/0959-5309/43/5/301
[89]	Prigogine, I. (1969) Structure, Disspipation and Life, First Conference on Theoretical Physics and Biology. North-Holland Publishing Company, Amsterdam.
[90]	Prigogine, I. (1967) Structure, Introduction to Thermodynamics of Irreversible Processes. 3rd Edition, Wiley-Interscience, New York.
[91]	Lu, S.S. and Liang, C.K. (1966) Experimental Evidence of Order-Disorder Transitions as of the Second Degree. Kexue Tongbao, 17, 441-443.
[92]	Sundman, B., Fries, S.G. and Oates, W.A. (1998) A Thermodynamic Assessment of the Au-Cu System. CALPHAD, 22, 335-354. https://doi.org/10.1016/S0364-5916(98)00034-0
[93]	Cao, W., Chang, Y.A., Zhu, J., Chen, S. and Oates, W.A. (2007) Thermodynamic Modeling of the Cu-Ag-Au System Using the Cluster/Site Approximation. Intermetallics, 15, 1438-1446. https://doi.org/10.1016/j.intermet.2007.05.003
[94]	Wei, S.H., Mbaye, A.A., Ferrira, L.G. and Alex, Z. (1987) First-Principles Calculations of the Phase Diagrams of Noble Metals: Cu-Au, Cu-Ag and Ag-Au. Physical Review B, 36, 4163-4185. https://doi.org/10.1103/PhysRevB.36.4163