High-Temperature Heat Capacity of the Na2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\text {Na}}_{2}$$\end{document} Molecule

The heat capacity of the Na2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\text {Na}}_{2}$$\end{document} molecule is investigated in the 2000 K–6000 K temperature range. At these temperatures, the classical treatment is appropriate. The contributions from the ground state and all the excited states dissociating to the 3s + 3s and 3s + 3p limits are taken into account. The results are compared with other known data and are discussed under the influence of unbound rovibrational and electronic excited states.


Introduction
Difficulties in the calculation of thermophysical properties at high temperatures are widely known. The NIST-JANAF Table [1] have many inaccuracies, for example, due to neglecting electronic excited states as in the case of the BBr molecule [2]. Spectacular problems with the partition functions of some molecules obtained by the summation of rovibrational levels were reported by Babou et al. [3] Heat capacity is a thermodynamic quantity which is very sensitive to all kinds of approximations (because of the second derivative in its relation with the partition function), making it a challenge to obtain reliable results. The other quantities can often be reliably calculated even if serious approximations, such as neglecting the rovibrational coupling, are used. It is particularly easy to calculate entropy [2,[4][5][6].
The unbound states (both metastable and scattering) also significantly influence thermodynamic quantities; the heat capacity is again the most sensitive. The recent tendency is to use all the unbound states in thermodynamic calculations [7][8][9]. Fortunately, in the case of the classical approach, it is the simplest option available but at high temperatures one also needs to know the purely repulsive PECs (potential energy curves) to describe all the interactions between the atoms.  This article presents the heat capacity of the Na 2 molecule in the temperature range  of 2000-6000 K. This molecule was chosen because the results for this weakly bound  molecule are vulnerable to various factors and approximations and because the results  known so far do not agree. Formulas (and calculations) of the internal partition functions and the internal heat capacities are in atomic units and the resulting heat capacities are given in J·(molK) −1 . All the calculations were performed with a Wolfram Mathematica computing system.
The rest of the article is organized as follows: Sect. 2 presents the methods for the classical calculation of heat capacity and the PECs used in those calculations. Section 3 contains the results and the discussion, and finally Sect. 4 contains the concluding remarks.

Calculation of Heat Capacity
Heat capacity is expressed as the second derivative of the molecular partition function: where Q int is the internal partition function of the molecule, (5/2)R is the contribution to the heat capacity resulting from the translational motion, and β = 1/(k B T ) is the inverse temperature. The internal partition function for a given PEC (V ) is calculated with the bound states approach [10] where σ is given by the condition V (σ ) = 0, γ is the incomplete gamma function, and μ is the reduced mass of the molecule; the approach with all unbound states [2,7] and finally with the approach excluding the negative contribution to the partition function [11] Q N E rovib =

The PECs of the Na 2 Molecule
In the present calculations all states with the 3s + 3s and 3s + 3p dissociation limits were taken into account in an exact or approximate manner. The two states with the limit 3s + 3s are X 1 + g and 1 3 + u . The excited states with the 3s + 3p limit are 3 u , [12]. There are also higher lying excited states but at the temperatures considered here they can be neglected.
As discussed later the accuracy of the PECs at low and high inter-atomic distances is crucial so it is important to use ab initio PECs fitted to actually calculated energies as in Zhang et al. [13], not ones which are based only on spectroscopic constants (which describes best the vicinity of the PEC minimum and are often empirical) as in Biolsi [9] and Song et al. [4].
In the present study, the analytical PECs fitted to the ab initio energies were used. For the ground X 1 + g state the PEC from Zhang et al. [13] (CBS) was used. This publication also reports the first excited state (1 3 + u ) but it has the wrong short-range limit (− ∞ instead of ∞).
All the excited states were approximated with the use of the Konowalow data [14] fitted to the expression: and the parameters in this expression are given in Table 1. The results were calculated with the explicit 1 3 + u (statistical weight 3), 3 u (statistical weight 6), and 1 + u (statistical weight 1) states. The 3 + g , 1 u , 1 + g states were approximated with 3 + g (with statistical weight 6), and the 1 g , 3 g , 2 3 + u states were approximated with the 1 g state (with statistical weight 11). The energy shift between the states dissociating to the 3s + 3s and 3s + 3p limits is ΔE = 0.07255E h [14]. Finally, the expression for the internal partition function of the Na 2 molecule is The above expression in its full form or with the chosen electronic states is used to obtain the results in the following section.

Results
Heat capacity-ground and the first excited state Comparing different approaches one can see that for the ground state at high temperatures inclusion of the unbound states plays a very significant role ( Table 2, the upper entries-the results are based on the CBS PES, for the lower entries-the results are based on the IRM PES). At 2500 K, the Table 1 Parameters of the fits to Konowalow data with Eq. 5  difference between the bound and unbound heat capacities is 4.81 J · (molK) −1 ; Biolsi [9] also reports a significant discrepancy but in this case the method of comparison is not specified. The results in Table 2 were restricted to 3000 K to minimize the influence of the excited states.
The difference between the IRM [4] and CBS PECs is also notable (the lower and upper entries): for the bound version of the heat capacity it is 0.61 J · (molK) −1 at 2000 K and for the unbound the corresponding discrepancy is 1.16 J · (molK) −1 . It means that the whole range is important for the results, and this fact is a guideline for PEC construction-it is important to provide accurate results not only in the vicinity of the minimum but also at the high and low inter-atomic distances. Note that Musiał et al. [12] and Zhang et al. [13] calculate the energies for inter-atomic distances only down to around 2a 0 but for Konowalow's results the lowest distance is 3.8a 0 (and highest distance is 15a 0 while the integration had to be done to 20a 0 or more for fully converged results). This problem is also important for shallow excited states. The low lying 1 3 + u state of Na 2 will, especially, be affected: the correct partition function for this state is negative at higher temperatures so that, in case of inaccuracies, even a wrong sign of contribution is possible. This can cause a serious problem for heat capacity which is a very sensitive quantity for all kind of effects [2,8,15].
For comparison, Table 2 shows the results of Song et al. [4] (approximate values read out from the plot) based on the IRM ground state PEC which differs because of two reasons-neglecting the unbound contribution and the separation of rotations and vibrations.
Finally, note that heat capacity is a sum of the internal and translational contributions so that all inaccuracies in a total value are relatively more pronounced in the internal contribution. The internal part of the heat capacity is more important than the translational contribution at lower temperatures, whereas at higher temperatures the translational part dominates (Fig. 1).
To investigate the accuracy of the lowest excited state (1 3 + u ) PEC, the respective partition functions were calculated at 2000 K and 6000 K for both PECs. First, the analytical PEC was obtained by the fit to Konowalow data. Second, the CBS PEC of Zhang et al. [13], because of the incorrect behavior at small inter-molecular distances, was approximated by the infinite wall at distances smaller than 2.5a 0 . The results are presented in Table 3. The results were compared with those of Mies and Julienne [7] for the partition functions for this electronic state. The infinite wall (hard spheres) approximation fails for this shallow potential at both temper- Fig. 1 Ratio of the internal to translational contributions to heat capacity C p,int /C p,trans -dependence on temperature atures and it cannot be used, unlike for the deeper ground state [11]. The other results are similar only at 6000 K, showing the importance of the shallow electronic states and that there is not necessarily a bigger discrepancy at higher temperatures. The heat capacities for the ground and the first excited state were calculated and the results are compared with the use of Konowalow fit to study the sensitivity of the heat capacity to the differences in the lowest excited state. The bound (excluding the negative contribution) and unbound versions of heat capacity are compared in Table  4. The NIST [1] data and Biolsi's results for the ground and first excited state are also displayed [9]. Inclusion of the first excited state changes all the results especially those with all unbound states (C U p ) even at the lowest temperature of 2000 K which is actually not surprising because it is a low excited state dissociating to the same limit as the ground state. The results in the present study are significantly below Biolsi's results, whereas the NIST data are mostly much higher than all the others. It should be noted that the NIST data also contain some of the higher electronic states.
Heat capacity-ground and all excited states The results obtained with all excited states are summarized in Table 5. The values here are generally closer to the NIST ones, whereas Biolsi's results are much lower at most temperatures.
In order to look for a possible source of the discrepancy, the results of this study were recalculated removing some of the electronic states. This test was done at a temperature of 6000 K. The results for the ground and the first excited state followed the Biolsi pattern, so the 3 g state that is missing in Biolsi's calculations was suspected. Removal of the 3 g state gives a value of 58.04, removing all states which were approximated (i.e., leaving three excited states) yields 52.80. Removing the two lowest states dissociating to the 3s + 3p limit ( 3 u , 1 + u ) from the calculation yields 46.63. Another way to analyze the influence of the excited states is the convergence diagram. Figure 2 shows the influence of consecutive terms describing the excited states partition functions in Eq. 6 for the unbound version of the heat capacity. As is already known, the inclusion of the first excited state lowers the heat capacity and adding more and more states increases and keeps the heat capacity at higher levels again. It is tempting to attribute that increase in value to the second excited state ( 3 u ) but excluding only this state from the full calculation also gives a similar value (51.38 J · (molK) −1 ).
The final conclusion about this behavior is that, in general, the higher excited states (dissociating to the 3s + 3p limit) increase the heat capacity again. Note that it is not possible to attribute this to any given electronic states because heat capacity is not additive with respect to the partition functions of the electronic states (the reason for that is the logarithm in Eq. 1), meaning it is not possible to split the total heat capacity between the heat capacities of the electronic states. Because of the second derivative in the definition of heat capacity, it can also be concluded that all higher excited states have similar second derivatives. The specific source of the discrepancy for the low values at high temperatures in Biolsi's results must be then be attributed generally to the differences in the higher electronic states. The difference in the partition functions of electronic states, which can seem unimportant, can lead to a big difference in the second derivative.
The present results of heat capacity with all the excited electronic states were fitted to the function similar to the Shomate equation, as in the NIST database, with the result C p = 751.766 − 0.589972T + 0.000187928T 2 − 2.58304 · 10 −8 T 3 + 1.28916 · 10 −12 T 4 − 4.15595 · 10 8 /T 2 , where T is the temperature in Kelvin. The difference between the calculated points and the values of this equation is in the range of 0.01 % to 0.3 %.

Conclusions
Heat capacity at high temperatures is a thermophysical quantity which is very sensitive to the accuracy of the partition function of all the electronic states. The exact calculations require knowledge of the exact PECs of many electronic states, as in case of weakly bounded molecules as the Na 2 molecule. In light of the present results, it is crucial to have more exact results than the old Konowalow data and to use analytical representations which are more sophisticated than the Hulburt-Hirschfelder potential used by Biolsi which is known to be unsatisfactory at low inter-atomic distances [16]. Both the present results and those of Biolsi certainly suffer from inaccuracies at the low inter-atomic part of PECs so that it is important that the PECs are calculated at very small inter-atomic distances. If the results at even higher temperatures are to be calculated reliably, the electronic states dissociating to the 3s+4s limit will also have to be taken into account, but unfortunately analytical PECs are not available for these.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.