Jammed packings of repulsive elastic spheres have emerged as a rich model system within which elastic properties of disordered glassy materials may be elucidated. Most of the work on these packings have focused on the case of vanishing temperature. Here, we explore the elastic properties of the associated connectivity network for finite temperatures, ignoring the breaking of bonds and the formation of new ones. Using extensive Monte Carlo simulations, we find that, as the temperature is increased, the resulting spring network shrinks and exhibits a rapidly softening bulk modulus via a cusp. Moreover, the shear modulus stiffens in a fixed volume ensemble but not in a fixed pressure ensemble. These counter-intuitive behaviors may be understood from the characteristic spectrum of soft modes near isostaticity, which resembles the spectrum of a rod near its buckling instability. Our results suggest a generic mechanism for negative thermal expansion coefficients in marginal solids. We discuss the consequences of bond breaking and an apparent analogy between thermalization and shear.Packings of repulsive elastic particles have emerged as a rich model granular system with potential relevance to amorphous solids [1,2]. A number of remarkable features emerge, in particular, when the packing is close to the onset of rigidity. One of the key characteristics of jammed packings is an excess of low-frequency vibrations [3], the so-called boson peak. Numerous consequences can be derived from this peculiar vibrational spectrum, with regards to, for instance, elastic or transport properties [4][5][6]. Some of these features are shared with lattices close to isostaticity [7][8][9][10][11][12], which may be exploited to develop meta-materials with novel mechanical properties [13,14].While most studies have focused on the zerotemperature consequences of the vibrational spectrum, we here study the impact of thermal fluctuations. Specifically, we consider the harmonic connectivity network obtained from a jammed packing of repulsive, frictionless spheres close to isostaticity, and study its mechanical properties as we heat up the system to a low but finite temperature.Elastic properties of ordered and disordered networks of springs at finite temperatures have been studied previously [15][16][17]. Most recently, motivated by the attractive properties of highly responsive marginal solids for material science and biophysics, spring networks have been studied near the isostatic threshold [10][11][12]18]. These studies revealed, amongst others, interesting anomalies in the entropic elasticity. While these studies have focused on networks that have soft bulk and shear moduli, as in rigidity percolation [19], it is a characteristic of jammed networks, studied in this work, to have a finite * Electronic address: ohallats@berkeley.edu bulk modulus at isostaticity [8]. As we will see, this has major consequences for the impact of thermal fluctuations on the material properties of the network.Simulation approach. To prepare the initial conditions...