PARAEQUILIBRIUM
This is an option that is available with the ADVANCED_OPTIONS command with the POLY and ED-EXP modules.
This calculates a paraequilibrium between two specific phases in an alloy system with one or more interstitial component(s) as fast diffusion species. Under the paraequilibrium state, two partially-equilibrated phases have the same chemical potential (but different contents) for one or more interstitial components (such as C, N, O, S, etc., as individual or combined)], along varied temperature or along a composition variable (of the matrix or one substitutional component) which has already set as the stepping variable with SET_AXIS_VARIABLE.
To ensure a successful point calculation of paraequilibrium state between two specific phases in a defined alloy system, it is important that you first have made a single-point equilibrium calculation with an initial overall composition in the current system before performing this advanced-option calculation; however, it is unnecessary to obtain an equilibrium in which either one or both of the target phases is stable. The initial overall composition must have a reasonable setting for the desired paraequilibrium calculation for the two target phases. This is especially true for cases where there are more than one interstitial components to be considered in the paraequilibrium state, because different interstitial components (for example C and N combined) may have significant different behaviors as partitioning into different structured phases; otherwise, for one chosen interstitial component the initial overall composition is OK for the paraequilibrium calculation between the specified two phases, but for other chosen interstitial component(s) it might be impossible to calculate the paraequilibrium state.
Note the following:
- Always check if the chosen phases A and B have the exactly same definition of elements and if the chosen interstitial components are all in the vacancy sublattice sites of the two phases; otherwise the program cannot find the paraequilibrium state (as it is impossible to correctly calculate u-fractions).
- Always have a comprehensive understanding of the normal phase diagram for the currently investigated system, so that you make the appropriate choice of the phase pair and staring bulk composition for the system.
- Always set the status of the chosen interstitial components as SPECIAL using the POLY command: CHANGE_STATUS
Component <interstitial component> = SPECIAL
. By doing this, you get a clear picture of u-fractions of various substitutional and interstitial components, which are different from the overall composition in the system. The SPECIAL status means that specified component(s) are not included in summations for mole or mass fractions. Therefore, all the composition variables plotted from paraequilibrium calculations are u-fraction related quantities.

Name of first phase: <Phase A>
Name of second phase: <Phase B>
The names of the two target phases A and B, between which the paraequilibrium state is to be calculated, must be entered subsequently or on the same (first) line at once then separated by an empty space, e.g. FCC#1 BCC
or FCC#2 M23C6
.
You need to understand what you are dealing with in terms of calculating a paraequilibrium state between the two specified phases.
Specifically, there are four distinguished cases to understand: (1) both chosen phases must have similar interstitial/vacancy sublattices where the fast-diffusion interstitial component(s) occupy; (2) the choice on the target phase pair must be reasonable for the defined system and specified initial overall composition; (3) both target phases should have phase constitution definitions that cover all the defined substitutional and interstitial components of the current alloy system; or (4) it is impossible to calculate the paraequilibrium state between the target phase pairs with given interstitial component(s) in the currently defined system.
Fast diffusing component: /c/: <interstitial component(s)>
Fast diffusing component: /none/: <interstitial component(s)>
The name(s) of the fast-diffusing component(s) (C as the default single component) must be given at the above prompts subsequently or at the same (first) prompt. It is possible to specify more than one interstitial component as fast diffusion species.
Note the following:
Such specified interstitial component(s) must be appropriately defined according to the phase constitution definitions of the two selected phases: these must be located on the interstitial/vacancy sublattices in both chosen phases;
If there is only one fast-diffusing component which is carbon, press the <Enter> key to accept the default input at the first prompt; if the single fast-diffusing component is another element (e.g. N), type its name at the first prompt;
If there are two or more fast-diffusing components (e.g. C
and N
), type their names at the above prompts subsequently or at the same (first) prompt (separated by an empty space, such as C N
);
To finish the input of fast-diffusing elements, accept NONE at a repeated prompt, i.e. by pressing <Enter> key to start the paraequilibrium point calculation;
If NONE or a non-existing component name is typed at the first prompt, it means no back diffusion is to be considered, and the para-equilibrium calculation is thus canceled entirely.
If the paraequilibrium state between the two specified phases is successfully calculated, the messages displays e.g.
NP(FCC) = 0.3586 with U-fractions C = 2.71821E-02 N = 4.1548129E-03
NP(BCC) = 0.6414 with U-fractions C = 7.10061E-04 N = 2.3781027E-04
All other compositions the same in both phases
Note: LIST-EQUILIBRIUM is not relevant
The first and second lines list the phase amounts expressed in mole-percent [NP(phase)] and the contents of the interstitial components C and N in a specific phase expressed in the so-called u-fractions [u-f(phase,C) and u-f(phase,N)], for the phase A (in this case as FCC) and phase B (in this case as BCC), respectively. The third line states that the compositions of the matrix component and all the remaining compositions (regarding substitutional components) in both the target phase A and target phase B are the same at the current paraequilibrium state, while these are not shown on screen. The last line indicates that after this advanced-option calculation the LIST_EQUILIBRIUM command is irrelevant and does not list the paraequilibrium state for the system at the current condition.
However, if the single-point calculation of the paraequilibrium state between the two specified phases has failed, these messages display:
*** Error 4 in ns01ad
*** Numerical error
This implies that the chosen target phase pair may be unreasonable for the defined alloy system or for the defined initial overall composition, or one or both phases may have inappropriate phase constitution definitions regarding the specified interstitial component(s). Then, you must either modify the settings of initial overall composition or specify the reasonable target phase pair with an appropriate choice of the fast diffusion interstitials in the defined alloy system.