Using a recently developed effective field theory for the interactions of nucleons at non-relativistic energies, we calculate the rate for the fusion process p + p → d + e + + ν e to leading order in the momentum expansion. Coulomb effects are included non-perturbatively in a systematic way. The resulting rate is independent of specific models for the strong interactions at short distances and is in agreement with the standard result in the zero-range approximation.The first step in the different nuclear processes in the Sun which generate the observed luminosity is proton-proton fusion p + p → d + e + + ν e [1]. It was explained more than sixty years ago by Bethe and Critchfield[2] when nuclear physics was still in its infancy. When the field had more matured, it was reconsidered in the light of more modern developments by Salpeter[3] and later by Bahcall and May [4]. But in spite of the enormous progress in nuclear physics during this time, the methods and approximations made in these different calculations were essentially the same. The obtained accuracy in the obtained fusion rate was just a few percent. Including higher order electromagnetic and strong corrections the uncertainty in the rate is now around one percent [5][6]. This is very impressive for a strongly interacting process at low energies very ordinary perturbation theory cannot be used.In the light of the importance this fundamental process plays in connection with the solar neutrino problem and possible neutrino oscillations[1], it is natural to reconsider this process from the point of view of modern quantum field theory instead of the old potential models used previously. A first attempt in this direction was made by Ivanov et al. [7]. They obtained then a result which was significantly different from the standard result based upon potential models. Subsequently it was pointed out by Bahcall and Kamionowski[8] that their effective nuclear interaction was not consistent with what is known about protonproton scattering at low energies where Coulomb effects are important.The approach of Ivanov et al. [7] is based upon relativistic field theory and should in principle yield reliable results. But it is well known that in particular for bound states like the deuteron it is very difficult to use consistently a relativistic formulation. Also the uncertain nuclear physics part of the fusion process under consideration takes place at