For a molecule you can calculate the polarizability alpha as p=alpha·E, where p is the dipole moment of the molecule and E the external electric field. So in VASP, you can set an external electric field with the EFIELD_PEAD tag, and don't forget to include the dipole corrections with the LDIPOL and IDIPOL tags.