Viral immune evasion by sequence variation is a significant barrier to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) vaccine design and coronavirus disease-2019 diffusion under lockdown are unpredictable with subsequent waves. Our group has developed a computational model rooted in physics to address this challenge, aiming to predict the fitness landscape of SARS-CoV-2 diffusion using a variant of the bidimensional Ising model (2DIMV) connected seasonally. The 2DIMV works in a closed system composed of limited interaction subjects and conditioned by only temperature changes. Markov chain Monte Carlo method shows that an increase in temperature implicates reduced virus diffusion and increased mobility, leading to increased virus diffusion.