Part I – January 1969 - Papers - Monte Carlo Calculations of Configurational Entropies in Interstitial Solid Solutions

Monte Carlo methods have been used to compute the arrangements of interstitial atoms dissolved in tetrahedral sites in bcc lattices. It is assumed that the presence of an interstitial atom "blocks " a certain number of neighboring sites and prevents their occupancy. Sites "blocked" by more than one filled site are allowed for. The computed values of. the mean occupation number (defined as the ratio of the total number of sites blocked to the number of solute atoms are used to calculate the configurational entropies of the solutions. These entropies are compared with those resulting from previous theoretical studies of this problem and also with available experin~ental data for the p Zr-H, Nb-H, V-H, and Ta-H systems. Evidence is also given that the "blocking" explanation of low limiting compositions in these systems, rather than this being due to initial limitations on the number of sites available, is probably correct. THE ideal partial configurational entropy of mixing of an interstitial solute in a metal is given by: where p is the number of interstitial sites per metal atom and Xi is the atomic fraction of the interstitial. For the bcc lattice. which we shall be concerned with in this paper, the interstitial positions are shown in Fig. 1. It can be seen that for the tetrahedral sites, p=6. whereas for the octahedral sites, p = 3. Different emphasis has been placed on the relative importance of energy and entropy effects in determining deviations from ideality in interstitial solid solutions. In some cases the same system, e.g., Fe-C, has been described by the contradictory regular and athermal solution models indicating that the enthalpy and entropy functions, derived from equilibrium data, are frequently not accu.rate enough to differentiate between these treatments. However, for certain metal-hydrogen solutions the equilibrium data is available over sufficiently wide ranges of temperature and composition to permit a reasonably accurate determination of the compositional variation of the heats and entropies. Hoch&apos; has attempted to interpret the results of interstitial solid solutions in terms of a regular solution model. In the case of the Ta-H system where 13 = 6, this model entails fitting the experimental relative partial entropies of solution, asH, to the equation: where ASgs is the relative partial excess entropy of solution of hydrogen. Hoch found that the results of Mallett and Koeh1 could be fitted to this equation with an approximately constant value of AF up to XH = 0.25. However, it is apparent from the solubility isotherms in this system which become asymptotic to the composition TaH that, since (Xh /6 - ~Xh ) becomes infinite only at TaH6, it is necessary that AS<&apos; tends to infinity at TaH. In other words, the low saturation composition of TaH, instead of the anticipated TaH,, eliminates the possibility of applying regular solution theory to such systems. Rather large negative excess configurational entropies must exist at higher hydrogen concentrations in order to explain the lower saturation values. To account for these low limiting compositions and excess entropies two distinctly different approaches have been followed. Rees and many others1-l2 have assumed that not all interstitial sites are crystallographically equivalent with respect to the interstitial addition; that is, in Eq. [I] p is less than the value anticipated from geometrical considerations. To describe, say, a bcc metal-hydrogen system with a limiting composition of MH by this approach one would consider that p = 1 in the first instance instead of p = 6.&apos;j3 In some cases, nonintegral values of B have been taken in order to improve the fit with the experimental data over limited ranges of composition. The other approach which has been used to explain the low saturation compositions is to assume that, although all sites are available for occupancy, strong repulsive interactions exist between the neighboring interstitial atoms, and hence occupancy of any site excludes or blocks a certain number of neighboring sites from being occupied. Earliest treatments of this concept considered the exclusion of an integral number, of nearest-neighbor sites from being occupied at all concentrations. In this case, the partial configurational entropy is given by: These early treatments failed to allow for the overlap of the blocked sites which will arise at all but the very lowest concentrations. More recently attempts have been made to calculate the effect of this decrease in the number of blocked sites on the configurational entropy. Using the quasichemical treatment of interstitial solid solutions as given by Lacher and assuming that an infinite repulsive interaction energy existed between the solute atoms. atom obtained an approximate configurational entropy applicable to the blocking with overlap case: