We investigate the large time behavior of the hot spots of the solution to the Cauchy problemwhere ϕ ∈ L 2 (R N , e |x| 2 /4 dx) and V = V (r) decays quadratically as r → ∞. In this paper, based on the arguments in [K. Ishige and A. Mukai, preprint (arXiv:1709.00809)], we classify the large time behavior of the hot spots of u and reveal the relationship between the behavior of the hot spots and the harmonic functions for −∆ + V .