I am just starting out with rietveld refinement. I have been trying to get a good fit for my data of the compounds Bi0.8La0.2Fe(1-x)Zrx using the software Fullprof. But in all my 4 compositions (x=0,0.02,0.05 and 0.1) there seem to be some intensity mismatch of the structural model and observed data. I am attaching the pcr file for the composition Bi0.8La0.2Fe0.98Zr0.02 and also the starting structural model. A picture showing the fit of the data is also attached. As can be seen, the peaks around 2theta = 22, 46 etc. show positive deviation (calculated value is higher) while some peaks have lower calculated value. My obtained chi^2=9.2. How can I improve the fitting ? I am thinking preferred orientation is not the cause since the same trend is seen in all four of my samples.