A theoretical investigation on the planarity of molecular structure of zinc phthalocyanine (ZnPc) and its aromaticity has been performed using B3LYP and M06-2X density functionals combined with selected Pople-type basis sets. The effect of the applied calculation method on the optimized structure of ZnPc and ZnPc•••H 2 O, both in the gas phase and in the polar solvent, was analyzed. To quantify the aromaticity of the ZnPc and ZnPc•••H 2 O complexes, both the geometric and magnetic criteria, i.e., Harmonic Oscillator Model of Aromaticity (HOMA) index and the nucleus-independent chemical shift (NICS) values at the centers or 1 Å above the centers of structural subunits, were calculated. The energies of highest energy occupied molecular orbital (HOMO) and lowest energy unoccupied molecular orbital (LUMO) and energy gaps were also estimated. The results show that the free ZnPC molecule is flat in the gas phase and nonplanar in the polar environments (DMSO and water). ZnPC•••H 2 O is nonpolar in the gas phase and polar solvent which is in agreement with recent X-ray reports. Both HOMA and NICS indexes indicate the presence of highly aromatic macrocycle and benzene rings while these parameters for pyrrolic ring are significantly smaller than in free pyrrole. The presence of polar solvents practically does not change aromaticity of the ring subunits of the studied compounds.