The Dynamical Cluster Approximation with Betts clusters is used to calculate the antiferromagnetic phase diagram of the 3D Hubbard model at half filling. Betts clusters are a set of periodic clusters which best reflect the properties of the lattice in the thermodynamic limit and provide an optimal finite-size scaling as a function of cluster size. Using a systematic finite-size scaling as a function of cluster space-time dimensions, we calculate the antiferromagnetic phase diagram. Our results are qualitatively consistent with the results of Staudt et al. [Eur. Phys. J. B 17 411 (2000)], but require the use of much smaller clusters: 48 compared to 1000.The accurate and efficient solution of lattice Hamiltonians such as the Hubbard model is a long standing challenge in the theoretical condensed matter community. These lattice models are routinely solved on a finite periodic lattice, for example with Monte Carlo, and the calculated properties extrapolated to the infinite limit. Due to the numerical expense in solving these models for large lattices, it is imperative to choose lattices that are efficient for the estimation and extrapolation of the physical properties of interest.In this paper we use the Dynamical Cluster Approximation (DCA) [1, 2, 3] (for a review see [4]) to explore the antiferromagnetic instability in the 3D Hubbard model at half filling, withwhere c ( †) iσ (creates) annihilates an electron with spin σ on site i, n iσ is the corresponding number operator, t the hopping amplitude between nearest neighbors i, j and U the on-site Coulomb repulsion. We solve this model on a series of finite clusters chosen according to the criteria proposed by Betts et al. [5,6]. We obtain converged results extrapolating from clusters of up to only 48 sites, which are in agreement with the calculations of Staudt et al. [7], who used conventional cubic lattices of up to 1000 sites and obtained the Néel temperature via the specific heat.To solve the Hamiltonian (1) we utilized the DCA [4]. For a 3D system the DCA maps the original lattice model onto a periodic cluster of size N c = L 3 c embedded in a selfconsistent host. Thus, correlations up to a range ξ < ∼ L c are treated directly, while the longer length scale physics is described at the mean-field level. With increasing cluster size, the DCA systematically interpolates between the single-site dynamical mean-field result and the exact result, while remaining in the thermodynamic limit. We solve the cluster problem using Quantum Monte Carlo (QMC) [8]. At half-filling there is no QMC sign problem; the only systematic error in the Monte Carlo is the time step error, which can be extrapolated away.In order to calculate the phase diagram of the system in the thermodynamic limit, we employ the scalingis the Néel temperature obtained from a DCA calculation with a cluster of linear cluster size L c . This form is justified if we envision the lattice as perfectly tiled by a periodic array of non-overlapping clusters. This system becomes ordered when the antiferromagnetic ...