Institute of Metals Division - Thermodynamics of Interstitial Solid Solutions with Repulsive Solute-Solute Interactions

An exact statistical treatment of a one-dimensional model is used as a basis for evoluating the reliability of certain simplified expressions for the activity of the solute in interstitial solutions, including one obtained from the exact expression by setting the repulsive interaction equal to infinity. The latter approximation is found to be satisfactory at low and moderate concentration if the repulsive interaction is large, even though not infinite. A similar expression (identical if the co-odination number is two) is derived from the quasichemical expression of Lacher, and is recommended as the best available expression for the excess configurational entropy of interstitial solutions with excluded sites. Some reasonable models are discussed, and the nature of the saturated solutions is determined by inspection. Some of the models reduce to the one -dimensional case. An analysis is given of the excess partial entropy of hydrogen in V-H; Nb-H; and To-H solutions. MOST treatments of the statistical thermodynamics of interstitial solid solutions have followed the classic paper1 of Lacher in making the simplifying assumption that the configurational entropy of the solution is ideal. However, it is becoming increasingly apparent that there are many interstitial solutions with very large so lute-solute repulsions, and for these the assumption of ideal entropy is not valid or useful. It is important to realize that with substitutional solutions large repulsions between the component atoms must lead to phase separation, whereas in interstitial solutions the free energy of the solution is not drastically increased by large solute-solute repulsions until intrinsic saturation is reached at the concentration where further solute would be forced to enter a site in which it would experience the repulsive effect of one or more solute atoms already present. In the limiting case of an infinitely large repulsive interaction, the excess free energy would be attributable entirely to excess entropy, the enthalpy of mixing being zero. AS will be shown below, even if the repulsions are less than infinite, a treatment based on an assumption of infinite repulsions may be very satisfactory up to moderately high concentrations of the interstitial component. Often in solutions where large repulsive interactions exist, there are also small interactions — often attractive—between solute atoms in configurations other than that corresponding to the large repulsion. In such cases the excess free energy will consist of an excess entropy term attributable to the large repulsive interactions, and an enthalpy term corresponding to the other small interactions. Nomenclature to differentiate succinctly between important cases would be a convenience. In this paper the nomenclature shown in Table I will be used. In Table I, and in the preceeding discussion, excess quantities are defined in terms of standard states which are pure solid solvent and pure (possibly hypothetical) solid saturated phase of the structure in question. In practice, it is more convenient to choose the interstitial element as a component, and its conventional standard state. This will add a composition-independent term to the excess entropy and the enthalpy. The earliest paper known to the present author which treats the thermodynamics of athermal interstitial solutions was given by schei12 in 1951, but the statistical derivations in that paper are open to criticism. Speiser and Spretnak were the first to give a correct statistical treatment,3 limited, however, to concentrations sufficiently low that the number of empty sites excluded from occupancy by more than one filled site is negligible. The purpose of the present paper is to extend the statistical treatment to more concentrated solutions, and to examine the magnitude of the errors introduced by assuming that the repulsive interactions are infinite when in fact they must be finite. THE QUASICHEMICAL APPROXIMATION Fortunately, a standard method already exists for taking into account the effect of large interactions upon the entropy of mixing. This is the quasi-chemical method, in which the probability of existence of a given pair of solute atoms in a certain proximate configuration is assumed to be proportional to exp(-w/kT), where w is the energy increase of the solution when the two atoms are moved from isolated locations in the solution to the configuration in question. A quasichemical treatment of interstitial solutions was given in 1937 in a widely neglected paper by Lacher.4 The result comes out