Density functional theory (DFT) is used to explore the structure, stability, and bonding in CSiGaAl 2 À/0 and CGeGaAl 2 À/0 systems having planar tetracoordinate carbon (ptC). The neutral systems have 17 valence electrons and the mono-anionic systems have 18 valence electrons. The ab initio molecular dynamics simulations for 2000 fs time at two different temperatures (300 and 500 K) supported the kinetic stability of the systems. From the natural bond orbital (NBO) analysis it is shown that there is a strong electron donation from the ligand atoms to the ptC atom. We have used Li + ion for the neutralization of the mono-anionic systems and more interestingly it does not disrupt the planar structure. The most preferable site for binding of Li + ion is along the Al Al bond in both of the mono-anionic systems. All the systems in this work have both σ and π aromaticity which is predicted from the computations of nucleus independent chemical shift (NICS). Although the anionic species obey the 18 valence electronic rule, the neutral systems break the rule with 17 valence electrons. However, both sets of systems are stable in the planar form. The bonding analysis of the systems includes molecular orbital, adaptive natural density partitioning (AdNDP), quantum theory of atoms in molecules (QTAIM), electron localization function (ELF) basin, and aromaticity analyses. The energy decomposition analysis (EDA) determines the interaction of Li + ion with CSiGaAl 2 À and CGeGaAl 2 À in Li@SiGaAl 2 and Li@GeGaAl 2 , respectively.ab initio molecular dynamics simulations, energy decomposition analysis (EDA), planar tetracoordinate carbon, Wiberg bond index (WBI), σ/π aromaticity