We combine Differential Evolution, a new technique, with the traditional grid based method for optimization of solar neutrino oscillation parameters $\Delta m^2$ and $\tan^{2}\theta$ for the case of two neutrinos. The Differential Evolution is a population based stochastic algorithm for optimization of real valued non-linear non-differentiable objective functions that has become very popular during the last decade. We calculate well known chi-square ($\chi^2$) function for neutrino oscillations for a grid of the parameters using total event rates of chlorine (Homestake), Gallax+GNO, SAGE, Superkamiokande and SNO detectors and theoretically calculated event rates. We find minimum $\chi^2$ values in different regions of the parameter space. We explore regions around these minima using Differential Evolution for the fine tuning of the parameters allowing even those values of the parameters which do not lie on any grid. We note as much as 4 times decrease in $\chi^2$ value in the SMA region and even better goodness-of-fit as compared to our grid-based results. All this indicates a way out of the impasse faced due to CPU limitations of the larger grid method.