Molecular dynamics simulation has been performed to study the effect of the polarizabilities of model anions on the ionic solvation in water clusters. The primary focus is given to the surface versus interior solvation behavior of the anions. To this end, various combinations of polarizable/nonpolarizable water and anion models were considered. Using the nonpolarizable TIP4P water with polarizable Cl- and Br- models, the Cl- is fully solvated inside the (H2O)(60) cluster, whereas the Br- is partially solvated at the surface of the cluster. However, when the polarizability of the Br- is turned off, the "Br-" anion is fully solvated. Using the polarizable Dang-Chang water, both Cl- and Br- reside at the surface of (H2O)(60) as well as (H2O)(500) clusters, consistent with the finding of Stuart and Berne [J. Phys. Chem. 100, 11934 (1996)] based on the polarizable TIP4P-FQ water with the polarizable Drude halide model. When the polarizabilities of the halide anions are turned off, the smaller size "Cl-" anion is fully solvated in the interior of the Dang-Chang water cluster, whereas the larger "Br-" anion is still partially solvated at the surface of the cluster, indicating the importance of the anion-size effect. We have also calculated the free energy change for the Cl- moving from the center of a lamella water slab to the surface. The free-energy change is on the order of 1 kcal/mol, indicating that the Cl- can easily access the surface region of the Dang-Chang water slab. (C) 2003 American Institute of Physics.