# Create a conjugated carbon molecule
mol1 = CarbonMolecule("Double Honeycomb")
# Add atoms
c1_1 = mol1.add_atom(x=-0.25, y=1)
c1_2 = mol1.add_atom(x=0.25, y=1)
c1_3 = mol1.add_atom(x=-0.5, y=0.5)
c1_4 = mol1.add_atom(x=0.5, y=0.5)
c1_5 = mol1.add_atom(x=-0.25, y=0)
c1_6 = mol1.add_atom(x=0.25, y=0)
c1_7 = mol1.add_atom(x=-0.5, y=-0.5)
c1_8 = mol1.add_atom(x=0.5, y=-0.5)
c1_9 = mol1.add_atom(x=-0.25, y=-1)
c1_10 = mol1.add_atom(x=0.25, y=-1)
# Add bonds
b1_1 = mol1.add_bond(c1_1, c1_2, order="single")
b1_2 = mol1.add_bond(c1_2, c1_4, order="single")
b1_3 = mol1.add_bond(c1_1, c1_3, order="single")
b1_4 = mol1.add_bond(c1_3, c1_5, order="single")
b1_5 = mol1.add_bond(c1_4, c1_6, order="single")
b1_6 = mol1.add_bond(c1_5, c1_6, order="single")
b1_7 = mol1.add_bond(c1_5, c1_7, order="single")
b1_8 = mol1.add_bond(c1_6, c1_8, order="single")
b1_9 = mol1.add_bond(c1_7, c1_9, order="single")
b1_10 = mol1.add_bond(c1_8, c1_10, order="single")
b1_11 = mol1.add_bond(c1_9, c1_10, order="single")
# Draw with Schemdraw
mol1.draw()
# Create a conjugated carbon molecule
mol2 = CarbonMolecule("Letter Envelope")
# Add atoms
c2_1 = mol2.add_atom(x=-0.8, y=0.5)
c2_2 = mol2.add_atom(x=0.8, y=0.5)
c2_3 = mol2.add_atom(x=-0.3, y=0)
c2_4 = mol2.add_atom(x=0.3, y=0)
c2_5 = mol2.add_atom(x=-0.8, y=-0.5)
c2_6 = mol2.add_atom(x=0.8, y=-0.5)
# Add bonds
b2_1 = mol2.add_bond(c2_1, c2_2, order="single")
b2_2 = mol2.add_bond(c2_1, c2_5, order="single")
b2_3 = mol2.add_bond(c2_2, c2_6, order="single")
b2_4 = mol2.add_bond(c2_1, c2_3, order="single")
b2_5 = mol2.add_bond(c2_2, c2_4, order="single")
b2_6 = mol2.add_bond(c2_3, c2_4, order="single")
b2_7 = mol2.add_bond(c2_3, c2_5, order="single")
b2_8 = mol2.add_bond(c2_4, c2_6, order="single")
b2_9 = mol2.add_bond(c2_5, c2_6, order="single")
# Draw with Schemdraw
mol2.draw()
# Create a conjugated carbon molecule
mol3 = CarbonMolecule("Triple Honeycomb")
# Add atoms
c3_1 = mol3.add_atom(x=-0.25, y=1)
c3_2 = mol3.add_atom(x=0.25, y=1)
c3_3 = mol3.add_atom(x=-0.5, y=0.5)
c3_4 = mol3.add_atom(x=0.5, y=0.5)
c3_5 = mol3.add_atom(x=-0.25, y=0)
c3_6 = mol3.add_atom(x=0.25, y=0)
c3_7 = mol3.add_atom(x=-0.5, y=-0.5)
c3_8 = mol3.add_atom(x=0.5, y=-0.5)
c3_9 = mol3.add_atom(x=-0.25, y=-1)
c3_10 = mol3.add_atom(x=0.25, y=-1)
c3_11 = mol3.add_atom(x=1, y=0.5)
c3_12 = mol3.add_atom(x=1.25, y=0)
c3_13 = mol3.add_atom(x=1, y=-0.5)
# Add bonds
b3_1 = mol3.add_bond(c3_1, c3_2, order="single")
b3_2 = mol3.add_bond(c3_2, c3_4, order="single")
b3_3 = mol3.add_bond(c3_1, c3_3, order="single")
b3_4 = mol3.add_bond(c3_3, c3_5, order="single")
b3_5 = mol3.add_bond(c3_4, c3_6, order="single")
b3_6 = mol3.add_bond(c3_5, c3_6, order="single")
b3_7 = mol3.add_bond(c3_5, c3_7, order="single")
b3_8 = mol3.add_bond(c3_6, c3_8, order="single")
b3_9 = mol3.add_bond(c3_7, c3_9, order="single")
b3_10 = mol3.add_bond(c3_8, c3_10, order="single")
b3_11 = mol3.add_bond(c3_9, c3_10, order="single")
b3_12 = mol3.add_bond(c3_4, c3_11, order="single")
b3_13 = mol3.add_bond(c3_11, c3_12, order="single")
b3_14 = mol3.add_bond(c3_12, c3_13, order="single")
b3_15 = mol3.add_bond(c3_13, c3_8, order="single")
# Draw with Schemdraw
mol3.draw()