The DCA
+algorithm was recently introduced by Stear, Maier, and Schulthess [Phys. Rev. B 88, 115101 (2013)] to extend the dynamic cluster approximation (DCA) with a continuous lattice self-energy in order to achieve better convergence with cluster size. Here we extend the DCA
+algorithm to the calculation of two-particle correlation functions by introducing irreducible vertex functions with continuous momentum dependence consistent with the DCA
+self-energy. This enables a significantly more controlled and reliable study of phase transitions than with the DCA. We test the new method by calculating the superconducting transition temperature
Tc in the attractive Hubbard model and show that it reproduces previous determinantal quantum Monte Carlo results. We then calculate
Tc in the doped repulsive Hubbard model, for which previous DCA calculations could only access the weak-coupling (
U=4t) regime for large clusters. We show that the new algorithm provides access to much larger clusters and delivers asymptotically converged results for
Tc for both the weak (
U=4t) and intermediate (
U=7t) coupling regimes, and thereby enables the accurate determination of the exact infinite cluster size result.
Ground-state properties of the asymmetric Hubbard model are studied on one-dimensional clusters (rings) with N=6, 8, and 10 sites. The ground-state energy, correlation functions, and phase diagrams are determined for various sets of model parameters. Computations are performed both by exact diagonalization (for N=6) and by an approximate method by which correlation effects can be examined on larger clusters than exact diagonalization allows. In the limiting cases (one limit corresponds to the spinless Falicov-Kimball model and the other to the Hubbard model) our results agree quite well with those obtained analytically for infinite systems.