We study the nonlinear Schrödinger equation with a periodic delta-function potential. This realizes a nonlinear Krönig-Penney model, with physical applications in the context of trapped Bose-Einstein condensate alkaly gases and in the transmission of signals in optical fibers. We find analytical solutions of zero-current Bloch states. Such wave-functions have the same periodicity of the potential, and, in the linear limit, reduce to the Bloch functions of the Krönig-Penney model. We also find new classes of solutions having a periodicity different from that of the external potential. We calculate the chemical potential of such states and compare it with the linear excitation spectrum.PACS numbers: 03.75. Kk, 03.75. Lm, 05.45.Yv Nonlinearity can deeply modify the Bloch theory of non-interacting atoms trapped in periodic potentials. Loop structures, energetic and dynamical instabilities, solitons and "generalized" Bloch states (i.e. states which do not share the same periodicity of the lattice), all arise in the context of a nonlinear Schrödinger (or GrossPitaevskii) equation. Applications span, for instance, the physics of dilute Bose-Einstein condensed gas trapped in optical lattices [1][2][3][4][5][6][7][8][9] or the propagation of signals in optical fibers [10] In this paper we study analitically the nonlinear Schrödinger equation with an external Krönig-Penney (KP) potential, given by a periodic array of deltafunctions. The linear Schrödinger equation with the same potential has been solved quite early in the '30, playing a distinguished role as a model in metal's theory [11,12].It is noticeable that also several properties of the nonlinear Schrödinger equation, with the same KP external potential, can be derived analitically. The most interesting results, however, are related with the emergence of new properties which do not have a counterpart in the linear case. Stationary solutions of the GPE which do not reduce to any of the eigenfunctions for a vanishing nonlinearity have been studied in [6] using a tight binding approximation, in and two-wells systems, in [13][14][15][16][17].The mean-field model of a quasi-1D BEC trapped in a KP potential is governed by the following nonlinear Schrödinger (or Gross-Pitaevski) equation:where µ is the chemical potential and g the nonlinear coupling constant. The KP external potential is given by equispaced delta-functions: V (x) = p ∞ n=−∞ δ(x−na), having a lattice constant a. Since the external potential has a step-like shape, it is useful to rewrite the GPE in hydrodynamic form. With ψ (x) = ρ (x) exp [−iΘ (x)] (and in dimensionless units), we have:, aρ → ρ and normalization n±1 n dx ρ(x) = 1. N 0 is the number of atoms in each well, and the integration constants α, β are fixed by the boundary conditions. In particolar, α has a simple physical meaning, being the current carried by the order parameter ψ(x):(where J a 2 h → J). Bloch state. We derive a class of stationary solutions of Eq.s(2) having (i) the same periodicity of the external potential (Bloch states), an...