The internal structure of stripes in the two dimensional Hubbard model is studied by going beyond the Hartree-Fock approximation. Partially filled stripes, consistent with experimental observations, are stabilized by quantum fluctuations, included through the Configuration Interaction method. Hopping of short regions of the stripes in the transverse direction is comparable to the bare hopping element. The integrated value of n k compares well with experimental results.By now it is well established that charged stripes are formed in a significant doping range of cuprate oxides [1][2][3]. The existence of these stripes was predicted, on the basis of mean field calculations, in advance of its observation, both in the one [4][5][6][7] and the three bands Hubbard model [8]. This is one of the scarce theoretical results in the field of high-T c superconductivity which was confirmed after its prediction. Interest in these calculations decreased, as it was generally understood that the Hartree-Fock approximation was unable to obtain the partially filled stripes observed experimentally. Unrelated calculations found stripes with different fillings in the t-J model [9], and, using different techniques, in the Hubbard model [10], although other numerical calculations show conflicting results [11,12]. Alternatively, it has been argued that stripes in doped Mott antiferromagnets arise from a tendency towards phase separation, frustrated by electrostatic interactions [13,14]. Stripes in the Hubbard model have also been analyzed within slaveboson techniques, which use the Hartree-Fock solutions as input [15].In the present work we show that the mean field calculations, initially used to demonstrate the existence of stripes, can be systematically improved in order to study the partially filled stripes observed experimentally. In addition, they provide significant insight into the internal structure of the stripes and their fluctuations in the transverse direction. The Hartree-Fock method can be considered a quasiclassical approximation to the spin and charge degrees of freedom. Their low amplitude quantum fluctuations can be incorporated by using the Random Phase Approximation [16]. In addition, one needs to consider quantum tunneling processes between degenerate, or nearly degenerate, Hartree-Fock solutions, when there are many. This is achieved with the Configuration Interaction method (CI), widely used in quantum chemistry [17,18]. The combination of the Hartree-Fock and CI methods gives reasonable results even when applied to one dimensional systems [19]. The CI method restores the symmetries broken by the Hartree-Fock approximation, and provides information on the quantum dynamics of the static solutions obtained in mean field, which can be broadly classified into spin polarons or stripes.Previous mean field studies have focused on filled stripes [7], which tend to be the solutions with the lowest energy per hole in this approximation, especially for the values of U/t ∼ 4 used in the initial studies [4][5][6][7]. There are,...