We have investigated the solvent effects on $Delta log;K_s $(the difference of stability constant of binding) and the different free energies of binding of $Nd^{3+}$ and $Eu^{3+}$ ions to 18-crown-6, i.e., the selectivity of $Nd^{3+}$ and $Eu^{3+}$ ions to 18-crown-6 using a Monte Carlo simulation of statistical perturbation theory (SPT) in diverse solvents. The stability constant ($Delta log;K_s $) of binding of $Nd^{3+}$ and $Eu^{3+}$ ions to 18-crown-6, in $CH_3OH$ was calculated in this study as -1.06 agrees well with the different experimental results of -0.44~-0.6, respectively. We have reported here the quantitative solvent-polarity relationships (QSPR) studied on the solvent effects the relative free energies of binding of $Nd^{3+}$ and $Eu^{3+}$ ions to 18-crown-6. From the calculated coefficients of QSPR, we have noted that solvent polarity (ET) and Kamlet -Tafts solvatochromic parameters (b ) dominate the differences in relative solvation Gibbs free energies of $Nd^{3+}$ and $Eu^{3+}$ ions but basicity (Bj) dominates the negative values in differences in the stability constant ($Delta log;K_s $) as well as the relative free energies of binding of $Nd^{3+}$ and $Eu^{3+}$ ions to 18-crown-6 and acidity (Aj) dominates the positive values in differences in the stability constant ($Delta log;K_s $) as well as the relative free energies of binding of $Nd^{3+}$ and $Eu^{3+}$ ions to 18-crown-6.