Application of the generalized gradient corrected functional within standard density-functional theory results in a dramatic failure for Au, leading to divergent thermodynamic properties well below the melting point. By combining the upsampled thermodynamic integration using Langevin dynamics technique with the random phase approximation, we show that inclusion of nonlocal many-body effects leads to a stabilization and to an excellent agreement with experiment. Most present-day ab initio databases contain results from T = 0 K density-functional-theory (DFT) calculations using standard exchange-correlation functionals. This is a limitation in two respects: (a) First, materials development and application takes place at finite temperatures at which entropic contributions become important. In particular, anharmonic phonon-phonon excitations have been revealed to strongly modify materials properties, e.g., for dynamically unstable systems [2] or for defect formation [3]. (b) Second, standard functionals are known to have intrinsic deficiencies, e.g., for elements with nearly full electron shells. Hybrid functionals provide an improvement by introducing a fraction of exact exchange [4][5][6], but the corresponding mixing parameter is not always well defined, especially for composite systems. Further improvement requires also the correlation energy to include nonlocal many-body effects. A possible route utilizes the random phase approximation (RPA) within the adiabatic-connection-fluctuation-dissipation theorem (ACFDT) [7]. RPA has been successfully applied to various systems [8][9][10] and recent developments aim at improving the computational efficiency [11,12]. However, despite the various developments and applications, RPA has been employed only at T = 0 K, whereas the impact of nonlocal many-body interactions at finite temperatures has remained unknown.In this Rapid Communication, we compute the material properties of Au-a prototype closed shell element-using the RPA within the ACFDT formalism up to the melting point. This is made possible by advancing our previously introduced upsampled thermodynamic integration using Langevin dynamics (UP-TILD) technique [13] towards a combination of molecular dynamics simulations performed at the standard DFT level in conjunction with upsampled snapshots computed within RPA. We shall demonstrate that the inclusion of nonlocal many-body effects remedies the deficiencies introduced by standard DFT.We start by investigating the thermodynamics of Au within the Perdew-Burke-Ernzerhof generalized gradient approximation ) in conjunction with the projectoraugmented-wave (PAW) formalism [15] as implemented in VASP [16,17] using the newly designed GW-PAW potentials [10]. We concentrate on the heat capacity, which is a representative thermodynamic quantity. The state-of-the-art approach to calculate its temperature dependence includes electronic excitations in conjunction with the quasiharmonic approximation. We find that electronic excitations give a negligible contribution [18] due to ...