ENERGETICS, STRUCTURE, AND COMPRESSIBILITY OF NAF DETERMINED BY THE PERIODIC HARTREE-FOCK METHOD

Basis set development methods for molecular cluster calculations were used to develop basis sets for periodic ab initio calculations of the energy-dependent properties of ionic crystals. Two extended, Gaussian basis sets for closed-shell sodium and fluorine ions were developed to compare with two minimal basis sets in calculating the crystal struct ures, enthalpies of formation, and high-pressure compression behavior of the B1 and B2 phases of NaF. Accurate calculation of the energy-dep endent properties to within experimental error required the use of the larger, extended basis sets. Enthalpies of formation for NaF-B1 at 0 K and 1 bar were calculated to be between - 921 and - 929 kJ/mol, whic h agrees well with an experimental value ( - 927 +/- 10 kJ/mol) estima ted from calorimetric data. For the B1 phase, the zero-pressure bulk m odulus was calculated to be between 0.485 and 0.499 Mbar, as compared to the measured value of 0.464 +/- 0.062 Mbar. The calculated molar vo lume at zero pressure agrees to within 0.1% of the experimental value. For the B2 phase, the zero-pressure bulk modulus was nearly the same as that for the B1 phase, ranging from 0.494 to 0.517 Mbar. However, t his calculated bulk modulus for the B2 phase was considerably smaller than the experimental value of 1.03 +/- 0.19 Mbar estimated from high- pressure compression measurements. Our calculations also suggest the V 02/V01 of the B2 phase is nearer to 0.871-0.895 than to the experiment al estimate of 0.81 1. The compression curves for both phases agreed, for the most part, with the experimental data to within the experiment al error. The calculated pressure for the NaF B1/B2 transition is betw een 0.26 and 0.29 Mbar, slightly higher than the averaged experimental value (0.23 +/- 0.03 Mbar) but within a pressure range where both pha ses metastably coexist along the compression curve. The calculated vol ume change accompanying the phase transition ( - 10.1%) is in good agr eement with experimental data (- 10.4%). The apparent discrepancies be tween the zero-pressure bulk modulus and the compression curves of the B2 phase are shown to be related to the range of molar volumes used i n deriving the zero-pressure properties from an equation of state.