We study Andreev states near atomic scale modulations in the pairing potential in both s-and d-wave superconductors with short coherence lengths. For a moderate reduction of the local gap, the states exist only close to the gap edge. If one allows for local sign changes of the order parameter, however, resonances can occur at energies close to the Fermi level. The local density of states (LDOS) around such pairing potential defects strongly resembles the patterns observed by tunneling measurements around Zn impurities in Bi2Sr2CaCu2O8+x (BSCCO). We discuss how this phase impurity model of the Zn LDOS pattern can be distinguished from other proposals experimentally.PACS numbers: 74.45.+c, 74.62.Dh Motivated by the experimental ability to determine the LDOS with high resolution in both energy and real space using scanning tunneling microscopy (STM), there has recently been a large interest in the perturbations caused by impurities in superconductors. This is because impurities disturb the underlying superconducting state and hence constitute a natural probe of the state in which they are embedded [1]. For magnetic adatoms on the surface of conventional s-wave superconductors such experiments were performed by Yazdani et al [2]. In the superconducting state of BSCCO, STM measurements have provided detailed information about the electronic structure near Ni and Zn impurities [3,4,5]. Near Zn it was found that each impurity generates a sharp conductance peak close to the Fermi level (ω B ∼ −1.5meV). By fixing the bias voltage between the tip and the sample at ω B , the spatial structure of this state can be mapped out: it is strongly localized near the Zn atom with maximum intensity on the impurity site and second largest intensity on the next-nearest neighbor sites.Although there exists a large number of theoretical treatments dealing with the resonant states generated by nonmagnetic impurities in d-wave superconductors, no consensus has been reached on their relevance to experiments [1]. Within the most straightforward scenario where Zn is modelled as a delta-function potential scatterer, low-energy resonant states are generated in the unitary limit where the potential is large compared to all other energy scales of the problem [6]. For realistic bands with particle-hole asymmetry, a particular value of the impurity potential V must be chosen to tune the resonance energy ω B to be close to the Fermi level. Thus, the potential scattering model can reproduce the energetics of the resonant state. However, as is well-known, it produces a low-energy spatial LDOS pattern with a severely suppressed amplitude on the impurity site, in contrast to experimental results on Zn [5]. Physically, this is clear since a large on-site potential V penalizes any substantial amplitude of the impurity-state wavefunction on the impurity site. In fact, within this model the largest intensity is found on the nearest neighbor site to the impurity.These properties of the delta-function potential scatterer are shown in Fig.1. The discrepanc...