Part IV – April 1968 - Communications - Discussion of "A Model for Concentrated Interstitial Solid Solutions; Its Application to Solutions of Carbon in Gamma Iron"*

On the basis of a statistical thermodynamic treatment of the data of smithz2 on the activity of carbon in austenite, a Darken and smith23 deduced that the interaction energy, wy, between carbon atoms occupying nearest-neighboring interstitial sites is repulsive. Aaronson, Domian, and pound 24 (ADP) have applied two different statistical thermodynamic analyses to the same data. One analysis, a composite of the treatments of Darken and Smith,23 Speiser and Spret-nak," and Kaufman, Radcliffe, and Cohen,26 primarily takes account of the positional entropic effects of wy in dilute solutions. The other, due to Lacher 27 and Fowler and Guggenheim," is also applicable to more concentrated solid solutions. The equations which both analyses produce for a,, can be rearranged so that, when In a,, is plotted against the appropriate function of the carbon content of the austenite, the plot will have a slope of unity when the correct value of wy is employed. Since no other constant need be simultaneously evaluated in order to make this determination, ascertaining the most probable value of wy from Smith's data on the bases of these analyses was an elementary problem in least-squares statistics. ADP demonstrated that the wy's obtained from both analyses vary with temperature in a statistically significant manner. McLellan, Gerrard, Horowitz, and sprague 29 (MGHS), on the other hand, have recently made a new statistical thermodynamic study of a,, and have concluded from this study that wy is independent of temperature. The equation which they derived for a,, is: where 0 = ratio of the mole fraction of carbon atoms to that of iron atoms, z = number of sites from which occupancy is excluded by the presence of a carbon atom and the repulsive character of its interaction with nearest-neighboring carbon atoms (z = zb in Ref. 24), and x = a geometric parameter taking account of overlapping of the exclusion shells of nearby carbon atoms. z is 1 when wy = 0 and 13 when wy = Q. For a given value of z there is a definite value of X . In the absence of an equation connecting z and X, MGHS computed geometrically the value of x at several integral values of z; x increases from 0.47 to 8.67 as z is increased from 2 to 5. MGHS fitted Eq. [I] to the a,, data of Ellis, Davidson, and Bodsworth30 at 925° and 1050°C and to the data of smithz2 at 1000°C. Confining z to integral values in the range 4 to 6, and also somewhat restricting AG, the values of z, X, and AG were sought which would minimize the average root mean square percent deviation, D, of the calculated a,, values from those obtained experimentally. Finding that D was almost always smallest when z = 5, they concluded that z is independent of temperature. Since:23, 26 where R = gas constant and T = absolute temperature, and w,, was considered to be similarly invariant, this statement cannot be strictly true. The temperature range employed, however, was small and the interval between the values of z tested was sufficiently large so that the results could still be consistent with a constant The statistical treatment of the a,, data applied by MGHS is clearly quite approximate. Since more than three values of a,, are available at each temperature, z, and X, and AG have actually been over deter mined, rather than underdetermined as their treatment implies. In order to make more complete and exact use of the a,, data, we shall employ the Legendre Method of Least Squares, as described by Whittaker and Robinson. This technique permits conversion of the overdetermined system of equations produced by the substitution into Eq. [I] of the available pairs of a,, and 0 data at a given temperature into an even determined system. Solution of the latter system yields the most probable value of each of the three constants. This procedure will allow a statistically rigorous test to be made of the conclusion drawn by MGHS in respect of the temperature dependence of In order to obtain linear iiequations of conditions" from Eq. [I], the exact values of the three unknown constants are replaced by G/RT = =/RT + y, where z, i, and aG are estimated values of the constants. This allows the conditional equations to be written in the form:'l where fi -f6, i, =/RT) = the right-hand side of Eq. [I] for a given pair of ayi and values. When applied to Eq. [I], these relationships become: The i equations of condition are reduced to three "normal equations" in the manner discussed by Whittaker and Robinson. The normal equations are then solved for a, ß, and y. Both of these operations are most conveniently performed by means of the computer technique of Winterbottom and Gjostein.32 This technique was used iteratively until the value of 0 (the most sensitive of the three correction factors) became 10.01. (Until they became small, a, 8, and y