A Monte Carlo replica exchange (MCRE) algorithm was used to compute the sodium cation distribution in bare faujasite zeolite with a number of cations per unit cell ranging from zero (Si:Al --> infinity) to 96 (Si:Al = 1). Eight independent realizations of the system were simulated simultaneously, in the temperature range 300-2325 K. The resulting distributions at room temperature were found to be in very good agreement with both available experiments and the analytical quasi-chemical model of Mortier and co-workers. A single canonical simulation at room temperature using site-to-site hopping yields identical results. The main advantage of the MCRE method is that no assumption is needed on the actual cation adsorption sites. One could thus, in principle, predict cation locations and distributions in a nanoporous solid for which the precise location of extraframework cations is not known.