An atomistic simulation of H-Pd(100) provided a phase diagram for the c2×2 H overlayer phase. The Embedded Atom Method (EAM) calculated energy of each configuration of atoms and the Metropolis Monte Carlo algorithm equilibrated the structure and generated configurations from which to sample the structure factor for the H overlayer. The procedure provided the expectation of the square of the structure factor modulus, < |S
2| >, as a function of temperature at three coverages. The inflection point of the < |S
2| > versus T curve estimated the critical temperature for disordering, T
c,, for one value of coverage, θ. The plot of T
c versus θ, the phase boundary for the c2×2 phase, lay about 125 K below the experimentally determined boundary. A comparison of the energies of ordered and disordered phases showed ΔE = 0.016 eV per hydrogen atom. Equating this unrealistically small energy difference to thermal kinetic energy (3/2)kBTc
at the critical temperature implies Tc
≈ 100 K. Obtaining – |S
2| > values relatively free of noise at such low temperatures required large numbers of Monte Carlo steps. The c2×2 phase is the experimentally determined stable low temperature phase, and was assumed to be the lowest-energy phase possible in this simulation. The very small ΔE indicates that some other ordered phase may be more stable than c2×2 in the EAM model.