A theoretical Density Functional Theory study was performed for a benzene molecule in water cages. Two DFT functionals (B3LYP and BLYP) were employed. The optimized geometries of the studied clusters were used to calculate the aromaticity of benzene in a condensed phase using the aromaticity indices: HOMA, NICS, PDI, and H. The results were compared with aromaticity of a single benzene molecule in the gas phase and in the solvent environment provided by the PCM continuum model. It is argued that high aromaticity of benzene in the gas phase is retained in the water environment. (C) 2014 Elsevier B.V. All rights reserved.