Heat capacity calculation of the CuAlS2 semiconductor compound using the Lambert W function

ABSTRACT In term of the Lambert W function, simple and reliable analytical expression for the estimation of the specific heat capacity of the CuAlS2 semiconductor is proposed by using the Lambert Boltzmann distribution. The obtained results have been compared with the corresponding experimental results in the 100-300 K temperature range. The estimated results are in good agreement with the experimental results over the entire temperature range. However, accurate values of the heat capacity cannot be calculated using the standard Boltzmann formula.


Introduction
Among with I-III-VI 2 semiconductors compounds, disulphide of copper and aluminum (CuAlS 2 ) is of special technological interest. It is known that these semiconductor compounds are suitable for many technical applications, such as light emitting devices [1], optical filters [2] and solar cells [3], they are also being studied as materials for Spintronics [4]. The CuAlS 2 semiconductor is also suitable for the use as window material for solar cells and light-emitting diodes since it has a wide direct gap of 3.5 eV.
For optimizing the synthesis of the CuAlS 2 semiconductor compound and estimating its thermodynamic characteristics, the evaluating data on the heat capacities of this semiconductor compound is of great importance. So far, the heat capacity of CuAlS 2 has been measured by the adiabatic calorimeter and were more refining using the least squares method [5]. However, no theoretical methods were used to calculate the heat capacity of CuAlS 2 .
The aim of the present study is to present a simple and reliable analytical expression for the estimation of the specific heat capacities of semiconductor compounds such as CuAlS 2 in terms of the Lambert W function using the Lambert Boltzmann distribution in the temperature range 100-300 K.

Heat capacities C v and C p expression and its application to the semiconductor compound CuAlS 2
The expression for internal energy U and heat capacity C v of indistinguishable particles are, respectively, The constant volume heat capacity is expressed as Here k is the Boltzmann constant, T the temperature, N the Avogadro number and Z is the partition function. This latter was determined for the Lambert-Boltzmann distribution which is the revised Boltzmann statistics for a finite number of particles established by Kakorin [6].
In the literature to derive the exponential Boltzmann distribution, the Stirling approximation, ln(N!) = N ln(N)-N is used. By applying a more exact version of Stirling's approximation and using the exact form log N!, respectively, Kakorin derives two more accurate versions of the Botlzmann distribution for a finite number of particles. The first of these involves the Lambert W function [7]. In particular, the Lambert partition function is [6], Here N W is the number of particles, is the Lambert maximum energy, β = 1/kT, ε = k θ E where θ E , k are, respectively, the Einstein temperature and the Boltzmann constant and W denotes the Lambert W function. This function solves the equation [7]. The Lambert W function allows the explicit solution of entire classes of differential equations, which actually only could be solved numerically and is experiencing today a renaissance in various fields of science and engineering [8][9][10][11][12][13][14][15][16][17][18].
The Lambert W function is a complex function with an infinite number of branches, only two of them having real values. If x is real, then for − 1 e ≤ x ≤ 0, there are two possible real values of W(x), as displayed in Figure 1.
Also, the C v versus C p relationship can be obtained from the Nernst-Lindemann equation [19,20] Implying that where A 0 = 5.1 × 10 −3 J −1 . K. mol and T m is the melting temperature. The square root in Equation (5) can be expanded conveniently into a Taylor series expansion [21], Inserting Equation (2) into Equation (1) yields and by replacing C v in Equation (6), we obtain the expression of the heat capacity Cp as: Numerical values of the principal Lambert W function can be computed using packages such as Mathematica and Maple.

Results and discussion
Firstly, the method was applied to calculate the values of the specific heat capacity at constant volume C v which we compared with experimental data [5].
The experimental and our calculated values of heat capacity C v are presented in Table 1. It is evident from this table that the result obtained is in excellent agreement with the experiment data [5] in the temperature range from 100 to 300 K with a large discrepancy of 0.7% at 120 K, 160 K and 250 K.
In an attempt to confirm the method of calculation proposed in this paper, the specific heat capacity at constant pressure C p for the semiconductor compound CuAlS 2 is calculated in the temperature range 100-300 K. In comparison with the experimental data available in the literature, the calculated result is shown   in Table 2. Our results are in excellent agreement with experimental data obtained by Korzun et al. [5] over the temperature range from 100 to 300 K with the largest discrepancy of 0.9% at 170 K. For further improvement for this method of calculation, Figures 2 and 3 show respectively the comparison of the value of heat capacity at constant pressure Cp and the heat capacity at constant volume C v calculated by the standard Boltzmann formula and the Lambert Boltzmann distribution with the experimental data [5]. The Boltzmann formula gives significantly larger discrepancies with the data and a large disagreement appears between both results. This confirms that the standard Boltzmann distribution is unable to accurately represent the specific heat.

Conclusion
A new method has been presented to accurately calculate the specific heat capacities of the semiconductor compound CuAlS 2 using the Lambert-Boltzmann distribution which is a generalization of the exponential Boltzmann distribution when the number of particles is finite. The calculated results of specific heat capacities relatively to the semiconductor compound CuAlS 2 have been compared with those reported in the literature and are found to be in good agreement. By comparison, the standard Boltzmann formula is systematically over-biased.

Disclosure statement
No potential conflict of interest was reported by the author.