The search for the ground state structures through geometry optimizations onto the ground state with minimal energies was performed using DFT-based plane-wave Vienna Ab initio Simulation Package (VASP) . For the carbon atomic potential including valence states, the projector augmented wave (PAW) method was applied . Treating at the same level the exchange X and the correlation C, the exchange-correlation (XC) effects were considered using a generalized gradient approximation (GGA) . The relaxation of the atoms onto ground state geometry was done applying a conjugate-gradient algorithm . Blöchl tetrahedron method with corrections according to Methfessel-Paxton scheme was applied for geometry optimization and energy calculations. A special k-point sampling was applied for approximating the reciprocal space Brillouin-zone (BZ) integrals. For better reliability, the optimization of the structural parameters was carried out along with successive self-consistent cycles with increasing mesh until the forces on atoms were less than 0.02 eV/Å and the stress components below 0.003 eV/Å 3 . Besides the elastic constants calculated to infer the mechanical stabilities and hardness, further calculations of phonon dispersion curves were also carried out to verify the dynamic stability of the new carbon allotropes. In the present work, the phonon modes were computed considering the harmonic approximation via finite displacements of the atoms around their equilibrium positions to obtain the forces from the summation over the different configurations. The phonon dispersion curves along the direction of the Brillouin zone are subsequently obtained using "Phonopy" interface code based on Python language . The structure representations were obtained thanks to VESTA graphical program .