II. METHOD III. RESULTS AND DISCUSSION function of frequency are presented. Finally the conclusions are given in Section IV. The independent electron susceptibility is defined as: where ϵ$_{nk}$ and f$_{nk}$ are the one-electron energies and the corresponding Fermi functions  and G   G ′ are reciprocal lattice vectors. In this work we accurately compute χ$_{0}$ by performing the summations in Eq. (1) using a random sampling over the Brillouin zone (BZ). We used ∼ 3000 independent k − points per band  chosen according to a stochastic algorithm which accumulates them around the FS for bands crossing E$_{F}$ . The final results are obtained by averaging over 40 runs each containing completely independent k -point set. This procedure nearly completely eliminates the numerical noise  and shows a good convergence both in terms of number of independent runs and number of k -points within a single run. We have included 65 bands in order to ensure convergence with respect to the number of empty bands. The energy bands and the matrix elements for LaOFeAs have been calculated  within the (spin-independent) local density approximation (LDA) 24 to the exchange-correlation functional. Our calculations have been done using the Full-potential Linearized Augmented Plane Wave (FPLAPW) method. This choice is necessary because of the extreme sensitivity of the electronic structure (in particular  of the bands close to E$_{F}$ ) to the method $^{25}$. No numerical approximation has been made in the evaluation of matrix elements in Eq. (1). Given the huge debate about the dependence of the results on the position of the As atom in the unit-cell $^{25 26}$  in the present work we show results both for the theoretically optimized ( z$_{As}$ = 0 . 638) and experimental 27 ( z$_{As}$ = 0 . 6513) position of the As atom. We first look at the band structure and the FS for the parent and the doped LaOFeAs. The ground state calculations are performed with the tetragonal unit cell  containing two Fe atoms  at the experimental lattice constants a = 4 . 03 ˚ A and c/a = 2 . 166. The energy bands 4 χ$_{0}$ ( q   G   G $^{′}$  ω ) = ∑ nn $^{′}$k f$_{nk}$ − f$_{n}$ $^{′}$k + q ϵ$_{nk}$ − ϵ$_{n}$ $^{′}$k + q + ¯ h ( ω + iη ) 〈 n $^{′}$k + q ∣ e i ( q + G ) · r ∣ n k 〉 〈 n k ∣ e − i ( q + G $^{′}$) · r ∣ n $^{′}$k + q 〉 (1)