Vibrational and thermodynamic properties of Ni3S2 polymorphs from first-principles calculations

We have calculated the compressional, vibrational, and thermodynamic properties of Ni₃S₂ heazlewoodite and the high-pressure orthorhombic phase (with Cmcm symmetry) using the generalized gradient approximation to the density functional theory in conjunction with the quasi-harmonic approximation. The predicted Raman frequencies of heazlewoodite are in good agreement with room-temperature measurements. The calculated thermodynamic properties of heazlewoodite at room conditions agree very well with experiments, but at high temperatures (especially above 500 K) the heat capacity data from experiments are significantly larger than the quasi-harmonic results, indicating that heazlewoodite is anharmonic. On the other hand, the obtained vibrational density of states of the orthorhombic phase at 20 GPa reveals a group of low-frequency vibrational modes which are absent in heazlewoodite. These low-frequency modes contribute substantially to thermal expansivity, heat capacity, entropy, and Grüneisen parameter of the orthorhombic phase. The calculated phase boundary between heazlewoodite and the orthorhombic phase is consistent with high-pressure experiments; the predicted transition pressure is 17.9 GPa at 300 K with a negative Clapeyron slope of −8.5 MPa/K.