Si-E bondings in hexacoordinate silicon phthalocyanine were analyzed using bond order (BO), energy partition, atoms in molecules (AIM), electron localization function (ELF), and localized orbital locator (LOL). Bond models were proposed to explain differences between hexacoordinate and tetracoordinate Si-E bondings. Aromaticity of silicon phthalocyanine was investigated using nucleus-independent chemical shift (NICS), harmonic oscillator model of aromaticity (HOMA), conceptual density functional theory (DFT), ring critical point (RCP) descriptors, and delocalization index (DI). Structure, energy, bonding, and aromaticity of tetracoordinate silicon phthalocyanine were studied and compared with hexacoordinate one.