Abstract

The present paper deals with the theoretical investigation of thermodynamical and structural properties like internal energy , entropy , Helmholtz free energy , isothermal compressibility , specific Heat , structure factor and long wave length limit of structure factor of 3d liquid transition metals. To describe electron-ion interaction we have used our newly constructed parameter free model potential. To perform this task, we have used different reference systems like Percus Yevick Hard Sphere (PYHS), One Component Plasma (OCP), and Charged Hard Sphere (CHS) reference systems. We have also seen the influence of different local field correction functions like Hartree (HR), Taylor (TR), and Sarkar et al. (SR) on thermodynamical properties of 3d liquid transition metals. Finally we conclude that the proper choice of the model potential along with reference system plays a vital role in the study of thermodynamical and structural properties of 3d liquid transition metals.

1. Introduction

The structure and thermodynamics of liquid metals have been widely studied with an increasing sophistication in the modelling of the interionic forces and in the classical statistical mechanics treatment of ionic correlations [148]. This study has drawn much theoretical attention both for its intrinsic interest and for the relevance to an understanding of electronic properties. The calculations of the atomic interactions in bonded metals have mainly been based on the density dependent pairwise potentials derived from electron-ion pseudopotential with the linear response theory and second order perturbation theory. The statistical mechanics side of the problem has been calculated with perturbation theory integral equations and computer simulations [110]. Although a lot of work has been done on both the structure and thermodynamics, still some questions await a definite answer. The limitations of linear response theory based interactions are well known. It is also possible to consider the interactions based on second order perturbation theory as effective pairwise potentials [11, 12]. However, in such case one would expect that a combined study of the structural and thermodynamic consequences of the interaction could illuminate in a quantitative way the relevance of such effective model interactions for the explanation of the physical and chemical properties of the liquid metals. We believe that only a combined analysis of structural and thermodynamic data should be able to assess the quality of any model for the interionic forces. However, only a few studies have been made addressing simultaneously these two problems.

In all these attempts the use of pseudopotential approach is found a remarkable success. But many existing studies have been limited to local model potentials with empirically determined parameters [1315]. Also it is found in our literature survey [150] that some abinitio pseudopotentials suitable for the perturbation theory of the ionic interactions were generated. In general, the pseudopotentials like Ashcroft empty core model (AS) [1619], Heine-Abarenkov model (HA) [17, 20, 21], harmonic model potential (HMP) [17, 20, 21], generalized nonlocal model potential (GNMP) [2224], energy independent nonlocal model potential (EINMP) [25, 26] and so forth. have been applied to study the structure and thermodynamics of liquids.

On the other side of the treatment of the atomic correlations, computer simulation has a precision limited only by numerical errors. However, when focusing on trends of thermodynamic quantities, approximate but reliable theories aimed at modelling the free energy are still very useful tools. The simplest of such method is Gibbs-Bogoliubov (GB) approach [24, 27, 28].

The general idea underlying the GB [24, 27, 28] inequality is that the true Helmholtz free energy of a liquid metal is bounded from above by the free energy of a suitably chosen reference system plus the difference in average potential energy between the actual and the reference system, calculated using the distribution functions of the reference systems. Thus the choice of the reference system is clearly very important. Ideally the basis for choosing the best reference system is that one expects it to give the lowest free energy and the expressions required for the thermodynamic and structure can be expressed possibly in an analytic form. People have applied popular reference system like Percus-Yevick Hard Sphere (PYHS) [23, 2938], One Component Plasma (OCP) [3744], Charged Hard Sphere (CHS) [29, 30, 4548], Hard Sphere Yukawa (HSY) [49, 50], Optimized Random Phase Approximation (ORPA) [24], and Soft Sphere (SS) systems [23].

Thus it is highly desirable to search a better reference system which gives good explanation of various properties when it is used with a particular model potential. Hence proper combination of a model potential and a reference system is one of the basic requirement in explaining various structure and thermodynamics of liquid metals. So, we thought it worthwhile to study the structure and thermodynamics of liquid transition metals using our own model potential using different reference systems.

In the last two decades, the considerable efforts have been made to understand the structure and thermodynamics of several simple liquid metals, liquid transition metals, liquid rare-earth metals, and their alloys [150]. The interatomic potentials of simple metals [24, 3739, 4549, 5155] have been fully investigated and their thermodynamic properties could be derived with sufficient accuracy, but in the case of transition metals the hybridization of d electron with s electron makes the things complex. Computer experiments are always intended to propose a plausible interpretation of experimental results in some cases to give the solution to an experimentally inaccessible problem. Despite the success of the theory in the solid state, results for the structure factor of liquid transition metals using molecular dynamics and other complicated liquid state theories have not been so reliable [110]. The reliability of the predicted values, however, entirely depends on the validity of a given interatomic potential and the model used. Compare to simple liquid metals only a few attempts are made on liquid transition metals [210, 50, 5558]. Therefore, the present paper deals with the computation of thermodynamical and structural properties like internal energy , entropy , Helmholtz free energy , isothermal compressibility , specific heat , structure factor , and long wave length limit of structure factor of 3d liquid transition metals using our newly constructed parameter free model potential [43, 59, 60] with different reference systems like Percus-Yevick Hard Sphere (PYHS) [23, 2938], One Component Plasma (OCP) [3744], and Charged Hard Sphere (CHS) [29, 30, 4548] systems.

An important application of pseudopotential is the calculation of structural properties of disordered materials such as liquid, amorphous, and their alloys. The problem with model pseudopotential is that of their transferability, because still there is no accurate method to obtain the form factor by which all the properties of liquid metals may be successfully investigated. It is found that a particular pseudopotential may be suitable for some properties of some metals and unsuitable for other properties of other metals. The usefulness of any model potential depends on how many number of parameters it involves. Generally, the potential involving less number of parameters is considered to be, comparatively a better one than that involving more number of parameters because it avoids more complexities in the calculation. It has been observed that a unique method of the determination of the potential parameter has not been pointed out so far. A number of investigators have used fitting procedure in which the potential parameters are fitted in such a way that a good agreement with experimental findings is to be obtained. Such a procedure will generally give good results for a certain property, but the same set of parameters will not give good results for other properties. Hence, we thought it worthwhile to construct a parameter free model potential which, by employing, explains the physical as well as chemical properties of condensed matter.

A pseudopotential method employing a simple model of a solid composed of atomic cores and valence electrons can predict the existence and properties of new solids and their properties [61]. By implementing this idea we have constructed a new model potential which is splitted into three regions [43, 59, 60],

Basically, this is the modified version of the Ashcroft’s empty core model [62]. Here we have considered actual core of an ion as an empty core rather than pseudizing it. The whole effective region is splitted into three parts. The effective weak potential felt by an electron is sandwiched between empty core and long range coulomb potential which is continuous at and is considered between the ionic radius () and atomic radius .

In the reciprocal space, the corresponding bare-ion form factor of the present model potential is given by [43, 59, 60]

Here , ,  , , , and are the valency, electronic charge, atomic volume, wave vector, ionic radius, and atomic radius, respectively. The local field correction functions like Hartree (HR) [63], Taylor (TR) [64], and Sarkar et al. (SR) [65] have been applied to see the influence of exchange and correlation effect on thermodynamical properties of liquid transition metals.

2. Theory

The Helmholtz free energy lies at the heart of the pseudopotential perturbation scheme to calculate the thermodynamical properties of liquid metals. The Helmholtz free energy of the system is defined by the use of the GB inequality as follow: where is the Helmholtz free energy (per ion) of the reference system, and denote the Hamiltonians of the real and reference system, respectively.

Again, with , , and are internal energy, entropy, and absolute temperature respectively.

The internal energy of reference system is given by is the first order electronic term arising from the average interaction of the valence electron with noncoulombic part of the bare ion pseudopotential. One has is the free energy of the electron gas. One has

Here is the density of the states per unit volume at Fermi energy. is the band structure energy and is given in terms of , the normalized energy wave number characteristic incorporating exchange and correlation corrections as,

Here is the structure factor of reference system. is the normalized energy wave number characteristic. It is expressed as

Here is called the modified dielectric function. It includes exchange and correlation effect in the screening through local field correction function . is the unscreened form factor. It depends on the core-ion pseudopotential and hence allows us to investigate its effect on the thermodynamics function. Here is the modified Hartree dielectric screening function which takes into account the conduction electron interaction

is the static dielectric function and is the correction factor for the exchange and correlated motion of the conduction electrons. is the Madelung energy and it can be written as follows. For PYHS systems For OCP system The constants , , , and are −0.89752, 0.94544, 0.17954, and −0.80049, respectively. is effective plasma parameter, , is the radius of the atomic sphere, . , where and are Boltzmann constant and absolute temperature, respectively.

For CHS systems where and .

The computation of the entropy is more straight forward and computed using with, where term represents the contribution from the different reference system. For PYHS system For OCP system For CHS system

For a given number density , packing fraction (HS, CHS, and SS reference systems) and plasma parameter (OCP reference system) are varied to obtain a minimum variational upper bound for the Helmholtz free energy .

The long wave length limit of structure factor is defined as follows.For PYHS For OCP system For CHS system

The isothermal compressibility and specific heat at constant volume can be expressed as

3. Results and Discussion

The thermodynamical variational properties of 3d series of liquid transition metals elements have been investigated without any adjustable parameter. Table 1 represents the input parameters used in the present computation. Self-consistent calculations by Moriarty [73] indicate that elements whose properties are affected by d bands, the valency lie in a narrow range . By implementing this concept, we have taken for 3d liquid transition metals, respectively. Therefore the number of d electrons, which do not participate to conduction, is fixed since is constant. We have displaced in Table 2 the values of the variational parameter with Hartree (HR) [63], Taylor (TR) [64] and Sarkar et al. (SR) [65] local field correction functions using Percus-Yevick Hard Sphere (PYHS) [23, 2938], One Component Plasma (OCP) [3744], and Charged Hard Sphere (CHS) [29, 30, 4548] reference systems for 3d liquid transition metals. The present results of variational parameter are quite good in agreement with experimental data [67] as well as other theoretical data [36, 45, 48, 50] whereas in the case of plasma parameter the experimental data is not available, but with other theoretical data [44, 45, 56] it is reasonable because we have taken actual plasma parameter rather than fitting it as done by others [45]. We have noticed that the effect of local field correction function on variational parameters ( and ) is very small. It is also found that the largest values of variational parameter ( and ) are obtained when the band is filled to a large extent, while lowest values are obtained when the band is nearly empty in the case of 3d liquid transition metal. This is clearly related to the strongest bonding which occurs when only the bonding but not the antibonding states are occupied. Figure 1 represents the variational parameters ( and ) with Hartree (HR) [63] local field correction for 3d liquid transition metals using different reference systems.

Figure 2 represents the Helmholtz free energy versus variational parameter of Cu liquid transition metals due to Sarkar et al. (SR) [65] local field correction function. The OCP system gives the lowest Helmholtz free energy compared to other reference systems. Figure 2 represents Helmholtz free energy versus variational parameter of liquid Cu for different local field correction function due to CHS system, which indicates that Taylor (TR) [64] local field correction function gives the lowest Helmholtz free energy, while Hartree (HR) [63] gives maximum Helmholtz free energy.

Figure 3 shows calculated and terms of 3d liquid transition metals. and terms are independent of the reference system and influence of local field correction function so that they remain constant for each system. Figure 4 represents Band structure energy and the effect of local field correction on term of 3d liquid transition metals. depends upon both reference system and local field correction functions. The OCP system gives the lowest value of Band structure energy. The maximum value of is obtained for 3d series for liquid Co. The deviation of local field correction on is going on increasing from Sc to Cu as the shell is going to fill up (or atomic number is going on increasing). The Taylor (TR) [64] shows less deviation of while Sarkar et al. (SR) [65] shows much deviation from Hartree (HR) [63] for CHS system. Figure 5 represents calculated value of Madelung energy term and effect of local field correction functions on 3d liquid transition metals for CHS system. The CHS system gives the lowest value of term while OCP system gives the highest values of term for liquid transition metals. From Figure 5, it is clear that term is not much sensitive to the choice of local field correction function. The calculated value of internal energy () of 3d liquid transition metals is tabulated in Table 3 along with other available data [10, 35] and also plotted in Figure 6 for Sarkar et al. (SR) [65] local field correction function. From these, it is clear that CHS system gives a lower internal energy () compares to other systems for all liquid transition metals. It is also noticed that the value of internal energy has not been much affected by the influence of local field correction functions. The and give negative contribution to internal energy while and give positive contribution to internal energy . term includes the electron-ion interaction through model potential. In this calculation major contribution comes from the term while contribution of is very small.

The calculated entropy of liquid transition metals is presented in Table 4. The minimum numerical values have been obtained due to the PYHS system while the maximum due to the OCP system. From the present results, we have observed that the OCP system gives excellent agreement with present results and experimental data [50, 68] for entropy compared to other reference system for 3d liquid transition metals (except Cu). A good agreement has been found with experimental data [50, 68] for liquid Ti, V, and Mn due to Sarkar et al. (SR) [65] local field correction function for OCP system while for liquid Fe, Co, and Ni due to Hartree (HR) [63] local field correction function for OCP system. The different contributions of entropy , are plotted in Figures 79 along with Sarkar et al. (SR) [65] local field correction function of 3d liquid transition metals. Among three contributions, only depends on structure and potential terms while and depend only on density and are independent of reference system and local field correction functions. The numerical value of is much larger than . The increases the absolute value of entropy whereas tends to decrease the total entropy of system. From Figures 79, it is clear that PYHS system gives the lowest value of term while OCP system gives the highest numerical value of term of liquid transition metals. It is also noticed that term is not much sensitive to the use of different local field correction function. We have compared our results with available theoretical data [10, 35]. The minimum numerical values of entropy have been obtained due to the PYHS system while the maximum due to the OCP system.

Table 5 represents Helmholtz free energy of 3d liquid transition metals. Figure 10 represents calculated Helmholtz free energy of 3d liquid transition metals due to Sarkar et al. (SR) [65] local field correction function. The OCP system gives minimum Helmholtz free energy compared to other reference system for liquid transition metals (except Co and Cu). For Co and Cu, the CHS system gives the minimum Helmholtz free energy. The values of Helmholtz free energy have not been much affected by the influence of local field correction functions. Hence we say that structure and potential are responsible to lower the free energy of the system. Therefore OCP and CHS systems show a lower Helmholtz free energy than PYHS. Among them OCP gives a lower Helmholtz free energy than other systems.

In our previous [60] attempt, we have thermodynamical properties like entropy , internal energy , and Helmholtz free energy of 3d liquid transition metals using variational principle based on the Gibbs-Bogolyubov inequality with Percus-Yevick hard sphere reference system. But, here we have used three different reference systems like Percus-Yevick Hard Sphere (PYHS) [23, 2938], One Component Plasma (OCP) [3744], and Charged Hard Sphere (CHS) [29, 30, 4548]. Our present results for thermodynamical properties are improved compared to our previous attempt [60]. Thus in present investigation, we have seen the influence of structure factor and different local field correction functions on thermodynamical and structural properties of liquid transition metals. It plays much more crucial role in the study of thermodynamical and structural properties of 3d liquid transition metals compare to our previous studies [60].

The structure factor obtained from the variational parameter ( and ) is illustrated in Figure 11. It is observed that in most of the cases overestimated the magnitude of the first peak, although in those cases the locations of the maxima and minima are reasonable predicted. We have explained the slight overestimation of our value when compared with experiment in cases where they occur, as being that in those cases where is greater than 0.45 (which is the value expected from liquid state theory). The implication of this is that effective hard sphere diameter is too large and by interference, the repulsive part of the pair potential in polyvalent metals extends too far [30]. CHS systems gives better results for structure factor compared to OCP system. OCP system gives much high peak in the structure factor compare to CHS system.

Table 6 represents calculated long wave length limit of structure factor of 3d liquid transition metals using our obtained variation parameters. The present results of long wave length limit is also compared with available experimental [71, 72] as well as other data [69, 70]. The present results are in a good agreement with available experimental data [71, 72] for Ti, V, Cr, Fe, and Cu liquid metals due to PYHS system, whereas for Sc and Mn due to CHS system and for Co and Ni liquid metals it is due to OCP system. Tables 7 and 8 show the calculated isothermal compressibility, and specific heat of 3d liquid transition metals, respectively, along with available experimental data [45, 71] and other theoretical data [45]. The good agreement with available experimental data [45, 71] has been obtained for isothermal compressibility due to CHS for Mn and Cu liquid metals, OCP for Sc, Co, and Ni liquid metals and PYHS for V and Cr liquid metals. It is found that good agreement with available experimental data [45] have been achieved for specific heat of Sc, V, Cr, Fe, Co, Ni liquid metals due to CHS system, whereas for Ti and Cu liquid metals it is due to PYHS and OCP systems respectively. So that overall CHS system is best suitable for calculating structural and thermodynamical properties of liquid transition metals along with our model potential compared to other reference system. This combination can be used to calculate other properties like electronic, magnetic, transport, and vibrational properties of liquid metals.

4. Summary and Conclusions

(i) The thermodynamical and structural properties of liquid transition metals have been successfully investigated along with our newly constructed parameter free model potential. (ii)These properties are very sensitive to the choice of reference systems. The good qualitative and quantitative data of the thermodynamical and structural properties of liquid transition metals have been obtained and would be beneficial to the other theoreticians as well as experimentalists working in the same field. (iii)These properties of liquid metals are not sensitive to the use of the different local field correction functions. (iv)Thus the application of our newly proposed model potential in the present study definitely adds new contribution to understand the thermodynamical and structural properties of transition metals in liquid state. (v)Thus the proper choice of the model potential along with the reference system plays a vital role in the study of the thermodynamical and structural properties of 3d liquid transition metals.

Acknowledgment

P. B. Thakor and Y. A. Sonvane acknowledge the financial support from the University Grants Commission, New Delhi, India, under a Major Research Project F. no. 33-26/2007 (SR).