We report on a numerical study of disorder effects in 2D d-wave BCS superconductors. We compare exact numerical solutions of the Bogoliubov-deGennes (BdG) equations for the density of states ρ(E) with the standard T-matrix approximation. Local suppression of the order parameter near impurity sites, which occurs in self-consistent solutions of the BdG equations, leads to apparent power law behavior ρ(E) ∼ |E| α with non-universal α over an energy scale comparable to the single impurity resonance energy Ω0. We show that the novel effects arise from static spatial correlations between the order parameter and the impurity distribution. 74.25.Bt,74.25.Jb,74.40.+k In spite of strong electronic correlations in the normal state, the superconducting state of high T c materials seems to be accurately described by a conventional BCS-like phenomenology. The debate over the k-space structure of the BCS order parameter ∆ k has now been resolved in favour of pairing states with d-wave symmetry. Since this symmmetry implies that ∆ k changes sign under rotation by π/2, there are necessarily points on the 2D Fermi surface at which ∆ k vanishes. The unique low-energy properties of high T c superconductors are determined by the quasiparticle excitations in the vicinity of these nodal points.For conventional s-wave superconductors, the density of states (DOS) ρ(E) has a well defined gap and is largely unaffected by non-magnetic disorder. In contrast, ρ(E) ∝ |E| in clean d-wave superconductors, and can be substantially altered by disorder. Much of the current understanding of disorder effects comes from perturbative theories, such as the widely-used self-consistent T-matrix approximation (SCTMA). In particular, the SCTMA is exact in the limit of a single impurity, and has been used in studies of the local DOS near an isolated scatterer [1]. For sufficiently strong scatterers, an isolated impurity introduces a pair of resonances at energies ±Ω 0 [2,3], where Ω 0 < ∆ is a function of of the impurity potential u 0 and of the band asymmetry. Analytic expressions for Ω 0 (u 0 ) have been given for a symmetric band [2], and in this special instance the unitary limit Ω 0 → 0 coincides with u 0 → ∞. For a realistic (asymmetric) band, the relationship is more complex [3].For a finite concentration of impurities n i , the SCTMA predicts that the impurity resonances broaden, with tails which overlap at the Fermi energy, leading to a finite residual DOS ρ(0) [3][4][5][6][7][8]. The region over which ρ(E) ≈ ρ(0), the "impurity band", and has a width comparable to the scattering rate γ at E = 0. In the Born limit, the ±Ω 0 resonances are widely separated in energy, and the overlap of their tails is exponentially small. In the strong-scattering limit, however, the overlap is substantial and γ ∼ √ Γ∆ d , where ∆ d is the magnitude of the d-wave gap, Γ = n i /πN 0 is the scattering rate in the normal state, and N 0 is the 2D normal-state DOS at the Fermi level. Several recent experiments [9][10][11][12] have studied quasiparticles in the impurity ...