Prediction of the relative stabilities and phase transition behavior of molecular crystalline polymorphs is highly coveted as distinct phases can possess different physical and chemical properties while having similar energies. Crystalline tetracyanoethylene (TCNE, C 6 N 4 ) is known to exhibit rich solid state phase behavior under different thermodynamic conditions, as demonstrated by a wealth of experimental studies on this system. Despite this fact, the role of temperature and kinetics on the phase diagram of TCNE remains poorly understood. Here, first-principles calculations and highresolution Fourier-transformed infrared (HR-FTIR) spectroscopy experiments are used to study the relative stabilities of the cubic and monoclinic phases of TCNE as a function of temperature. Specifically, density-functional theory with the van der Waals interactions method of Tkatchenko and Scheffler (DFT+vdW) is employed. The accuracy of this approach is demonstrated by the excellent agreement between the calculated and experimental structures. We find that the cubic phase is the most stable polymorph at 0 K, but becomes less favorable than the monoclinic phase at 160 K. This temperature-induced phase transition is explained on the basis of varying close contacts and vibrational entropies as a function of temperature. These findings are supported by a temperaturedependent HR-FTIR linewidth study of the CMN vibrons.