Phase competition underlies many remarkable and technologically important phenomena in transition metal oxides. Vanadium dioxide (VO2) exhibits a first-order metal-insulator transition (MIT) near room temperature, where conductivity is suppressed and the lattice changes from tetragonal to monoclinic on cooling. Ongoing attempts to explain this coupled structural and electronic transition begin with two alternative starting points: a Peierls MIT driven by instabilities in electron-lattice dynamics and a Mott MIT where strong electron-electron correlations drive charge localization. A key missing piece of the VO2 puzzle is the role of lattice vibrations. Moreover, a comprehensive thermodynamic treatment must integrate both entropic and energetic aspects of the transition. Here we report that the entropy driving the MIT in VO2 is dominated by strongly anharmonic phonons rather than electronic contributions, and provide a direct determination of phonon dispersions. Our ab initio calculations identify softer bonding in the tetragonal phase, relative to the monoclinic phase, as the origin of the large vibrational entropy stabilizing the metallic rutile phase. They further reveal how a balance between higher entropy in the metal and orbital-driven lower energy in the insulator fully describes the thermodynamic forces controlling the MIT. Our study illustrates the critical role of anharmonic lattice dynamics in metal oxide phase competition, and provides guidance for the predictive design of new materials.

We use quantitative experimental and theoretical approaches to characterize the vibrational dynamics of the Fe atom in porphyrins designed to model heme protein active sites. Nuclear resonance vibrational spectroscopy (NRVS) yields frequencies, amplitudes, and directions for 57Fe vibrations in a series of ferrous nitrosyl porphyrins, which provide a benchmark for evaluation of quantum chemical vibrational calculations. Detailed normal mode predictions result from DFT calculations on ferrous nitrosyl tetraphenylporphyrin Fe(TPP)(NO), its cation [Fe(TPP)(NO)]+, and ferrous nitrosyl porphine Fe(P)(NO). Differing functionals lead to significant variability in the predicted Fe-NO bond length and frequency for Fe(TPP)(NO). Otherwise, quantitative comparison of calculated and measured Fe dynamics on an absolute scale reveals good overall agreement, suggesting that DFT calculations provide a reliable guide to the character of observed Fe vibrational modes. These include a series of modes involving Fe motion in the plane of the porphyrin, which are rarely identified using infrared and Raman spectroscopies. The NO binding geometry breaks the four-fold symmetry of the Fe environment, and the resulting frequency splittings of the in-plane modes predicted for Fe(TPP)(NO) agree with observations. In contrast to expectations of a simple three-body model, mode energy remains localized on the FeNO fragment for only two modes, an N-O stretch and a mode with mixed Fe-NO stretch and FeNO bend character. Bending of the FeNO unit also contributes to several of the in-plane modes, but no primary FeNO bending mode is identified for Fe(TPP)(NO). Vibrations associated with hindered rotation of the NO and heme doming are predicted at low frequencies, where Fe motion perpendicular to the heme is identified experimentally at 73 and 128 cm-1. Identification of the latter two modes is a crucial first step toward quantifying the reactive energetics of Fe porphyrins and heme proteins.