By using ab initio molecular dynamics calculations, we show that even where the graphene lattice constant contracts, as previously reported for freestanding graphene below room temperature, the average carbon-carbon distance increases with temperature, in both free and supported graphene. This results in a larger corrugation at higher temperature, which can affect the interaction between graphene and the supporting substrate. For a weakly interacting system as graphene=Irð111Þ, we confirm the results using an experimental approach which gives direct access to interatomic distances. DOI: 10.1103/PhysRevLett.106.135501 PACS numbers: 81.05.ue, 68.65.Àk, 65.80.Ck Graphene's (GR) mechanical [1] and thermal [2,3] properties are crucial for applications such as the cooling of electronic devices [4,5], but their microscopic mechanism is often not well understood [6], in contrast to that of many electronic phenomena [7][8][9][10]. A particular challenge lies in the fact that graphene, while being two dimensional, exists in a three-dimensional world, permitting low-lying vibrational excitations and the formation of large-scale ripples perpendicular to the plane [11,12]. This is also important because it introduces a difference between freestanding GR and the technologically more important supported material.In particular, the thermal expansion coefficient of GR and the underlying microscopic mechanism is currently being debated. For freestanding GR and below % 500-700 K, a state-of-the-art atomistic simulation [13], a nonequilibrium Green's function approach [14], and a harmonic density functional theory (HDFT) calculation [15,16] all show a decreasing in-plane lattice parameter a, i.e., a thermal contraction. Above 900 K or so, the two former techniques show a trend reversal with a thermal expansion, while in HDFT the contraction persists over the entire temperature range studied (up to 2500 K). Experiments clearly confirm the thermal contraction of a below room temperature [3,17], but no high temperature measurements have been reported.In order to clarify this issue and to gain insight into its microscopic origin, we present a study based on ab initio simulations and core level photoelectron spectroscopy. Temperature-dependent ab initio molecular dynamics (AIMD) calculations were performed on both freestanding and supported GR. We have used the VASP code [18] with the projector-augmented wave method [19,20], the Perdew-Burke-Ernzerhof exchange-correlation energy [21], and an efficient extrapolation for the charge density [22]. Single particle orbitals were expanded in plane waves with a cutoff of 400 eV. We used the NPT ensemble (constant particles number N, pressure P, and temperature T), as recently implemented in VASP [23,24]. For the present slab calculations, we only applied the constant pressure algorithm to the two lattice vectors parallel to the surface, leaving the third unchanged during the simulation.In the case of freestanding GR the calculations were done using unit cells of different sizes (8 Â 8, 10 Â...