Part VIII – August 1969 – Papers - On the Stability of Spherically Symmetric Compositional Inhomogeneities in Solid Solutions

A convolution procedure is described for obtaining an analytical expression for the composition profiles of solute concentration:; in binary solid solutions. The free energy of solutions containing such compositional inhomogeneities can then be evaluated and the stability of spherically symmetric clusters (G.P. zones, for short) us a function of their radii and central compositions can be studied accordingly. Two characteristic zone radii are defined: the minimum radius (RM) corresponding to the smallest zone for which growth is energetically favored over dissolution, and the critical radius RK) corresponding to zones associated with the minimum work of formation. The values of RM and RK coincide at the two classical limits of nucleation and growth and spinodal decomposition. As the spinodal is approached, however, RM and RK diverge increasingly and it is proposed that the initial] decomposition morphology is governed more closely by the minimum, rather than by the critical radius since the former possesses no singularity within the m eta stable phase boundary. Reversion is briefly discussed in the light of the foregoing and it is shown that a set of reversion temperatures located between the coherent miscibility gap and the coherent spinodal must correspond to a given distribution of low temperature zones. SINCE the early works of Guinier1 and preston,2 the kinetics and morphology of GP zones have been studied extensively, see Ref. 3, for example, for a recent review. Nevertheless, a completely satisfying theoretical treatment of the decomposition of meta-stable or unstable binary solid solutions is still lacking due to the almost unsurmountable mathematical difficulties that one encounters in attempting to solve the pertinent diffusion equation which, in principle, describes the whole process. This equation, first derived by Cahn for the initial stages of spinodal decomposition,4 is a nonlinear partial differential equation of the fourth order in anisotropic space, a general analytical solution of which is completely out of the question. From this equation, some prediction can be made, however, concerning the later stages of spinodal decomposition5 and a recent numerical solution6 of the equation in the case of one- and two-dimensional composition modulations of a solid solution inside the spinodal has led to the following considerations: the sinusoidal morphology7 is to be expected only for solutions of average composition c, located near the center of the coherent miscibility gap.8 For solutions with co closer to the spinodal decomposition, depleted shell zones9 should constitute the first recognizable decomposition product, later evolving to iso- lated, enriched preprecipitates in a uniformly depleted matrix as this must, in all cases, constitute the final metastable morphology, barring, of course, loss of coherency and change in crystal structure. In agreement with the qualitative description of Bonfiglioli and Guinier,l0 it thus appears that the various models discussed by previous authors are, in fact, idealized successive stages of the general decomposition process. Unfortunately, corresponding calculations for a solid solution outside the spinodal are not yet available because of two additional difficulties: it is essential that the diffusion equation be solved in three-dimensional space, and furthermore, one must include in the equation an appropriate "source term" describing the statistics of the creation of composition fluctuations required to trigger the nucleation events. At present, a general treatment of this problem does not appear to be tractable. Nevertheless, by using a modelistic approach, i.e., by selecting a priori certain composition variation profiles and by evaluating their associated free energies, one can draw certain conclusions relative to the stability of these isolated composition modulations or "zones", both outside and inside the spinodal. This is the object of the present investigation. THEORY We shall now describe a standard four-parameter spherically symmetric GP zone, evaluate its free energy according to the Cahn and Hilliard11 formula (Eq. [5], below), and plot the corresponding energy surfaces as a function of zone amplitude and radius as