Methodologies beyond the Born-Oppenheimer (BO) approximation are nowadays important to explain high precision spectroscopic measurements. Most previous evaluations of the BO correction are, however, focused on light-element molecules and based on a nonrelativistic Hamiltonian, so no information about the BO approximation (BOA) breakdown in heavy-element molecules is available. The present work is the first to investigate the BOA breakdown for the entire periodic table, by considering scalar relativistic effects in the Diagonal BO correction (DBOC). In closed shell atoms, the relativistic EDBOC scales as Z(1.25) and the nonrelativistic EDBOC scales as Z(1.17), where Z is the atomic number. Hence, we found that EDBOC becomes larger in heavy element atoms and molecules, and the relativistic EDBOC increases faster than nonrelativistic EDBOC. We have further investigated the DBOC effects on properties such as potential energy curves, spectroscopic parameters, and various energetic properties. The DBOC effects for these properties are mostly affected by the lightest atom in the molecule. Hence, in X2 or XAt molecule (X = H, Li, Na, K, Rb, and Cs) the effect of DBOC systematically decreases when X becomes heavier but in HX molecules, the effect of DBOC seems relatively similar among all the molecules.