THERMOPHYSICAL PROPERTIES OF THORIUM COMPOUNDS FROM FIRST PRINCIPLES Vinayak Mishra a,* and Shashank Chaturvedi a a Computational Analysis Division, Bhabha Atomic Research Centre, Visakhapatnam 530012, India * Email of corresponding author: vinayakm@barc.gov.in Email: chats@barc.gov.in ABSTRACT Thorium carbide and nitride are potential candidates for their use as fuel materials in fast breeder reactors. Therefore, knowledge of their thermo-physical properties at high temperatures is necessary. In this paper, we present results of the first principles calculations of properties such as specific heat, volume and bulk modulus at high temperatures. The all-electron FPLAPW method has been combined with the quasi harmonic approximation for performing these calculations. Our results are in reasonably good agreement with published experimental results. Keywords: Thorium compounds, Specific heat, High temperature bulk modulus, FPLAPW, QHA. INTRODUCTION Thorium-based materials such as its binary compounds formed with light elements C and N are currently being investigated in relation with their potential utilization in fast breeder reactors, owing to their high melting point, high thermal conductivity, high density and good compatibility with the coolant (liquid Na) [1-2]. In order to use these materials as advanced fuels, it is very important to know their thermo-physical properties at high temperatures, which are required for modelling fuelpin behaviour. Only a few authors have attempted theoretical studies of thermophysical properties of these materials at high temperature. Their studies of these properties are either based on pseudo-potentials, or on molecular dynamics with empirical potentials. In the present work, we report results of a few useful properties of ThN and ThC, calculated from all electron FPLAPW + QHA method. COMPUTATIONAL DETAILS The calculations have been performed within the framework of density functional theory employing full potential linearized augmented plane wave (FPLAPW) method, as implemented in the WIEN 2k code [3]. Exchange correlation potentials based on LDA and GGA both have been used. Perdew - Wang parameterization [4] of LDA and Perdew - Burke - Ernzerhof parameterization [5] of GGA have been used. The basis function is expanded up to R MT K max = 9, R MT is the muffin tin radius and K max is the plane wave cut off parameter. The R MT values were selected as 2.2 for Th and 1.6 for C and N. The self consistent cycles were run until the energy convergence criterion of 10 5 Ry was reached. 5000 k - points were taken in the calculations. Thermal effects have been included using quasi-harmonic approximation. The (p, T ) equilibrium states are obtained by minimizing Gibbs free energy(g), with respect to V. Where, G can be expressed as G(V; p, T ) = E(V) + pv + A vib (Θ(V); T ) (1) The contribution A vib can be calculated by using quasi-harmonic Debye model. A vib (Θ(V); T ) = nkt [9Θ/8T + 3ln (1 e ( Θ/T) ) D (Θ/T )] (2) where n is the number of atoms per formula unit, k is the Boltzmann s constant, Θ is the Debye temperature and D(Θ/T) is the Debye integral. Once the equilibrium state for a given (p, T ) has
been obtained, then the required thermodynamic properties can be determined by using the corresponding volume in the appropriate thermodynamic equations. RESULTS AND DISCUSSION We first performed calculations of the total energy as a function of volume using the FPLAPW method. Next, we fitted the calculated E-V data with Birch-Murnaghan equation of state, in order to obtain lattice constant (a), bulk modulus (B) and pressure derivative of the bulk modulus (B P ) at ground state. In the case of ThN, the a, B and B P results, obtained from the BM fitting of our GGA (LDA) calculations are 5.178 (5.1055) A, 174.6 (200.2) GPa and 4.11 (3.32) respectively. The experimental results of these properties taken from the literature [6, 7] are: a = 5.16 A, B = 175 GPa and B P = 4.0. In the case of ThC, our GGA (LDA) calculations predict the a, B and B P as 5.355 (5.268) A, 131.9 (141.4) GPa and 3.0 (2.8) respectively. The published experimental results a and B of ThC are 5.344 A and 125 GPa respectively [1]. Thus, LDA tends to over-bind, as a result LDA lattice constants are too small and bulk moduli are too large. Though GGA overcorrects the lattice constant but the GGA results are in better agreement as compared to LDA. Therefore, in our subsequent calculations we have used GGA for the exchange and correlation. Fig. 1 : Comparison of total energies of B1 and B2 phases of ThC and ThN. B1 phase is stable near normal density and in expansion. B2 phase becomes stable at compression. The materials with B 1 phase can transform to B 2 phase on applying pressure. Such transformations may restrict the use of the materials as they lead to sudden change in the material properties. In order to test the stability of the B 1 phase of ThC and ThN in the compression and expansion, we have calculated the volume dependent total energies of their B 2 phase also. The total energies of these two structures have been compared in the Fig. 1. We can see that near normal volume and at expanded volumes, the B 1 phase of both materials is stable. However, at compressed volumes, 28.4 A 3 for ThC and 25.1 A 3 for ThN, the total energy of B 2 phase becomes smaller than the total energy of their B 1 phase. This indicates a phase transformation at compression. In order to find pressures at which this B 1 - B 2 phase transition takes place, we have compared their Gibbs free energies at various pressures. We have found that for ThC this transition takes place at 42 GPa and for ThN this transition takes place at 75.5 GPa. Hence these materials can be used without any problem up to these pressures. As these materials are used at very high temperatures, it is necessary to understand the high temperature behavior of the bulk modulus of these materials. Equilibrium volumes are obtained for each (p T ) by minimizing the Gibbs free energy as described in the section computational
details. Using the (p V) points for each T in the following formula one can calculate the bulk modulus at high temperatures. B T = V ( p V ) T (3) Fig. 2 : Variation of volume and bulk modulus of ThC and ThN with temperature. Top part of this figure represents variation of the volume per formula unit and bottom part represents the variation of bulk modulus. The evolution of the volume per formula unit of ThC and ThN with temperature is reported in Fig. 2 for zero pressure. We can observe in Fig. 2 that the volume per formula unit increases with temperature, as expected. The temperature evolution of the bulk modulus, have also been shown in the bottom part of the Fig. 2. Bulk modulus decreases with the temperature. We have fitted the temperature evolution of the bulk modulus with the Wachtmann formula [8]: B T = B 0 b 1 T e (-T0/T) (4) where B 0 is the bulk modulus at absolute zero, T is the temperature, and b 1 and T 0 are fitting parameters. At very high temperature the temperature dependence becomes linear according to this formula which is consistent with the obtained bulk modulus versus temperature data. We can see that the fitting is excellent for ThC. For ThN the fitting is slightly deviated. The fitting parameters are b 1 = 0.0154 GPa/K, T 0 = 118. 1 K for ThC and b 1 = 0.0363 Gpa/K, T 0 = 46.5 K for ThN. Temperature dependent bulk modulus is not available in the literature therefore our results can serve as prediction.
Heat capacity C p is important as it is one of the thermodynamic functions. Also, it is essential for the calculation of thermal conductivity from thermal diffusivity. The heat capacity (C p ) can be calculated by using following equation C p = C V (1 + αγt ) (5) Here, α is the thermal expansion coefficient, γ is the Gruneisen parameter and C V is the specific heat at constant volume which can be calculated by using following formula 3(Θ/T ) C V = 3 4 D(Θ/T ) nk[ e (Θ/T ) 1] (6) Figs. 3 shows the calculated specific heats at constant pressure and volume respectively for zero pressure. The temperature dependence of C V and C p follows a T 3 law at low temperatures. As is usual, C V tends to the 3R limit for higher values of T, where R is the gas constant. For the temperatures larger than 300K, Cp shows a nearly linear dependence with temperature. Our calculated values are compared with the experimental data from [9-11]. We can see reasonably good agreement between the calculated and experimental results. Fig. 3 : Cp, and CV results of ThC and ThN from present calculations, Experimental Results are from Samsonov et. al. [9], GTT -data [10] and Danan [11]. CONCLUSIONS In summary, we have presented results of our first principles calculations of the ground state properties such as a, B and B P ; of ThC and ThN. We have also tested the energetic stability of the B 1 phase, which is the ground state phase of both compounds, as compared to their B 2 phase and have found that at high pressure B 2 phase becomes stable. The phase transition pressures are 42 and 75.5 GPa for ThC and ThN respectively. We have also calculated the temperature dependent thermo-physical properties such as volume, bulk modulus and specific heats. Our results of specific heat C p are in reasonably good agreement with the available experimental data. REFERENCES [1] H. Kleykamp, Thorium carbides, Gmelin Handbook of Inorganic and Organometallic Chemistry, Eighth ed. Thorium Supplement, Vol. C3, Springer, Berlin (1987).
[2] R. Benz, A. Naoumidis, Thorium Compounds with Nitrogen, Gmelin Handbook of Inorganic Chemistry, eighth. ed. Thorium Supplement, vol. C6, Springer Berlin, (1992). [3] Blaha P, Schwarz K, Madsen GKH, Kvasnika D, Luitz J. WIEN 2k-userguide: An augmented plane wave plus local orbital program for calculating crystal properties, Viena:Viena University of Technology (Austria); [2014]. Available from http://www.wien2k.at/reg user/textbooks/usersguide.pdf/. [4] J. P. Perdew and Y. Wang, Accurate and simple analytic representation of the electron-gas correlation energy, Phys. Rev. B, 45 13244, (1992). [5] J. P. Perdew, S. Burke and M. Ernzerhof, Generalized Gradient Approximation Made Simple, Phys. Rev. Lett. 77, 3865, (1996). [6] F. A. Wedgwood, Actinide chalcogenides and pnictides. III. Optical-phonon frequency determination in UX and ThX compounds by neutron scattering, J. Phys. C: Solid State Phys. 7, 3203 (1974). [7] L. Gerward, J. S. Olsen, U. Benedict, J. P. Itié and J. C. Spirlet, The crystal structure and the equation of state of thorium nitride for pressures up to 47 GPa, J. Appl. Cryst. 18, 339, (1985). [8] J. B. Wachtman Jr, W. E. Tefft, D. G. Lam Jr, C. S. Apstein, Exponential temperature dependence of young s modulus for several oxides, Phys. Rev., 122, 1754, (1961). [9] G. V. Samsonov and I. M. Vanitskii, Handbook of refractory compounds, Plenum publishing corporation, New York (1982). [10] The SGTE pure substance and Solution Database, GTT-Data Services, (http://gtt.mch.rwthaachen.de/gtt-web/data/sgte-databases) (1996). [11] J. Danan, Chaleur specifique de 2 a 300 K du monocarbure de thorium (French), Specific heat of thorium monocarbide from 2 to 300 K (English translation), J. Nucl. Mater. 57, 280, (1975).