  By exploring the properties of the energy landscape of a bidisperse system of
soft harmonic disks in two dimensions we determine the thermal jamming
transition. To be specific, we study whether the ground state of the system
where the particle do not overlap can be reached within a reasonable time.
Starting with random initial configurations, the energy landscape is probed by
energy minimization steps as in case of athermal jamming and in addition steps
where an energy barrier can be crossed with a small but non-zero probability.
For random initial conditions we find that as a function of packing fraction
the thermal jamming transition, i.e. the transition from a state where all
overlaps can be removed to an effectively non-ergodic state where one cannot
get rid of the overlaps, occurs at a packing fraction of $\phi_G=0.74$, which
is smaller than the transition packing fraction of athermal jamming at
$\phi_J=0.842$. Furthermore, we show that the thermal jamming transition is in
the universality class of directed percolation and therefore is fundamentally
different from the athermal jamming transition.
