Diels-Alder (DÀ A) reaction shaped both the art and science of total synthesis to some degree, while the effect of temperature on DÀ A activity over polyoxometalates supported single-atom catalysts (SACs) has been infrequently studied and simulated using theoretical calculations. Herein, cycloaddition of 1,3butadiene (C 4 H 6 ) and ethylene (C 2 H 4 ) was employed as a model DÀ A reaction. The multitudes of SACs M 1 /PTA (M=Mn, Fe, Co, Ru, Rh, Pd, Os, Ir and Pt; PTA=[PW 12 O 40 ] 3À ) were examined by DFT-M06l computations to understand the reaction mechanism on a molecular level. The adsorption energies of reactant and product, and activation energy barriers for all the studied SAC systems have the same variation trends with the temperature variations. Considering that the adsorption for C 4 H 6 is always stronger than that of C 2 H 4 in all the studied systems, the initial adsorption configurations is the M 1 /PTA SACs adsorbed one C 4 H 6 molecule. Three SACs, namely the Co 1 /PTA and Rh 1 /PTA at 100 K, Rh 1 /PTA at 300 K were identified, which show predominant catalytic activity and the corresponding activation energy barriers are 4.21, 8.51 and 5.11 kcal mol À 1 , respectively. The bonding interaction between adsorbate C 4 H 6 and SACs arises from the occupied molecular orbitals (MOs) with a mixture of π orbitals of C 4 H 6 and d atomic orbitals of the metal single atom. These theoretical calculations give new guidelines to develop high catalytic activity and cost-effective SACs towards the DÀ A reaction.