We study a class of discrete focusing nonlinear Schrödinger equations (DNLS) with general nonlocal interactions. We prove the existence of onsite and offsite discrete solitary waves, which bifurcate from the trivial solution at the endpoint frequency of the continuous spectrum of linear dispersive waves. We also prove exponential smallness, in the frequency-distance to the bifurcation point, of the Peierls-Nabarro energy barrier (PNB), as measured by the difference in Hamiltonian or mass functionals evaluated on the onsite and offsite states. These results extend those of the authors for the case of nearest neighbor interactions to a large class of nonlocal short-range and long-range interactions. The appearance of distinct onsite and offsite states is a consequence of the breaking of continuous spatial translation invariance. The PNB plays a role in the dynamics of energy transport in such nonlinear Hamiltonian lattice systems.Our class of nonlocal interactions is defined in terms of coupling coefficients, Jm, where m P Z is the lattice site index, with Jm » m´1´2 s , s P r1, 8q and Jm " e´γ |m| , s " 8, γ ą 0, (Kac-Baker). For s ě 1, the bifurcation is seeded by solutions of the (effective / homogenized) cubic focusing nonlinear Schrödinger equation (NLS). However, for 1{4 ă s ă 1, the bifurcation is controlled by the fractional nonlinear Schrödinger equation, FNLS, with p´∆q s replacing´∆. The proof is based on a Lyapunov-Schmidt reduction strategy applied to a momentum space formulation. The PN barrier bounds require appropriate uniform decay estimates for the discrete Fourier transform of DNLS discrete solitary waves. A key role is also played by non-degeneracy of the ground state of FNLS, recently proved by Frank, Lenzmann & Silvestre.Key words. Discrete nonlocal nonlinear Schrödinger equation, onsite and offsite solitary waves, bifurcation from continuous spectrum, Peierls-Nabarro energy barrier, intrinsic localized modes Here, L is a linear and nonlocal interaction operator:(1.2)The range of the nonlocal interaction is characterized by the rate of decay of the coupling sequence, tJ m u. Nonlocal interactions arising in applications typically have coupling sequences, J m " J s m with polynomial decay J s m » m´1´2 s , 0 ă s ă 8, or exponential (Kac-Baker) decay: J 8 m » e´γ m , γ ą 0. We refer to the case s ą 1 as the short-range interaction case, and the case 0 ă s ă 1 as long-range interaction case. The case s " 1 is the critical or marginal range interaction. The case J 0 " 0, J˘1 " 1 and J m " 0, |m| ě 1, the operator L corresponds to the nearest-neighbor discrete Laplacian. 1 arXiv:1601.04598v2 [nlin.PS] 9 Jan 2017 J pQq for λ s pαq´1 ď |Q| ď π{α, (λ s pαq ą α and λ s pαq Ó 0 as α Ñ 0 ) in terms of those for 0 ď |Q| ď λ s pαq´1 (low frequency components of z E α,σ J). This yields a closed system for the low-frequency components, which we study perturbatively about the continuum FNLS limit using the implicit function theorem.