I have tried to do it quite a long. I have used around 50 different basis sets for each computation, but the results resemble quite scattered. I could not able to find any functional which can replicate the experimental result. I have scaled the value for water molecule ( experimental polarizability= 9.72) and only through the basis aug-cc-pvdz I have got the value for polarizability of water as 9.64, the closest among all values with an experiment. Now, I tried to solve L_Dopa molecule with the same functional, but the calculation did not converge.