Thermodynamic properties of graphene bilayers

[EN] Thermodynamic properties of graphene bilayers are studied by path-integral molecular dynamics (PIMD) simulations, considering quantization of vibrational modes and anharmonic effects. Bilayer graphene has been studied at temperatures between 12 and 1500 K for zero external stress, using the LCBOPII effective potential. We concentrate on the thermal expansion, in-plane and out-of-plane compressibility, and specific heat. Additional insight into the meaning of our results for bilayer graphene is obtained from a comparison with data obtained from PIMD simulations for monolayer graphene and graphite. They are also analyzed in view of experimental data for graphite. Zero-point and thermal effects on the in-plane and "real" area of bilayer graphene are studied. The thermal expansion coefficient αxy of the in-plane area is negative at low temperatures and positive for T≳800K. The minimum αxy is -6.6×10-6K-1 at T≈220K. Both in-plane (χxy) and out-of-plane (χz) compressibilities of graphene bilayers are found to increase for rising temperature, and turn out to be lower than that corresponding to monolayer graphene and higher than those found for graphite. At 300 K, we find for the bilayer χxy=9.5×10-2Å2/eV and χz=2.97×10-2GPa-1. Results for the specific heat obtained from the simulations are compared with those given by a harmonic approximation for the vibrational modes. This approach is noticeably accurate at temperatures lower than 200 K.The authors acknowledge the help of J. H. Los in the implementation of the LCBOPII potential. This work was supported by Ministerio de Ciencia, Innovación y Universidades (Spain) through Grants No. FIS2015-64222-C2 and No. PGC2018-096955-B-C44