Plane wave based ab initio DTF all valence electron spin polarized calculations are reported for the electronic structure and geometry of differently shaped graphene molecules. These polycyclic hydrocarbon molecules comprise four series: D 2h symmetric linear polyacenes C 4m+2 H 2m+4 (m = 2,,..,25); D 6h hexagulenes with zigzag edges C 6m**2 H 6m (m = 2,..., 10); D6h hexagulenes with crenelated edges C 6(3m**2-3m+1) H 6(2m-1) (m = 2,..., 6); D 3h zigzag edged triangulenes C m**2+4m+1 H 6m (m = 2,... 15). The systems variously display ground states that are spin paired singlet (S = 0), singlet anti-ferromagnetically ordered diradical and spin polarized S = 1/2(m - 1). Molecules with zigzag edges show evidence of electron delocalization along the perimeters with some bond alternation at the comers. In the acenes the spin paired singlet state of the short members switches to a singlet diradical at m ≈ 7-8 and this remains as the ground state for larger m. In contrast the triangulenes are magnetic and the atomic charge and spin density changes monotonically with distance from center to perimeter. In the hexagonal systems the development of a graphene core region, where the C-C bonds are 142pm, extends to within a few C-C bonds of the perimeter atoms. Zigzag hexagulenes have a spin paired singlet ground state for m ≤ 8 and a singlet diradical ground state for larger m. When hexagulenes are substitutionally doped with boron or nitrogen, acceptor or donor levels appear that track the valence or conduction band edges with increase in zigzag number. This result suggests the possibility of building several semiconductor device structures into the same graphene molecule.