We have developed a perturbation theory for the calculation of Jastrow factors, where the perturbation energies converge to the variational energy of the Jastrow wave function. The perturbation methods are designed for linear as well as exponential Jastrow ansatz, which are closely related to many-body perturbation theory and coupled cluster methods. Test calculations are performed for homogeneous Fermi gas models with screened Coulomb interaction (Yukawa potential) e^{- μ r}/r and with periodic boundary conditions (supercell approach). The out-coming Jastrow factors are compared with the ones calculated by Fermi hyper-netted-chain method (FHNC). The correlation energy is compared with standard MP2 result, and is judged by variational Monte Carlo (VQMC) and pure diffusion Monte Carlo calculations (PDMC).