def METHOD_NAME():
    plot = figure(
        width=150, height=150,
        x_axis_type=None, y_axis_type=None,
        x_range = (-8,8), y_range = (-8,8),
        toolbar_location=None, outline_line_color=None,
        match_aspect=True,
    )
    groups = []
    for sc in df["electron shells"]:
        n = len(sc)
        group = RendererGroup(visible=False)
        groups.append(group)
        c0 = plot.circle(x=0, y=0, radius=list(range(1, n+1)), fill_color=None, line_color="black", visible=False)
        c0.group = group
        xs = np.array([])
        ys = np.array([])
        for i, c in enumerate(sc):
            da = 360/c
            r = i + 1
            A = np.radians(np.arange(0, 360, da))
            x = r*np.cos(A)
            y = r*np.sin(A)
            xs = np.append(xs, x)
            ys = np.append(ys, y)
            #plot.circle(x=x, y=y)
            #plot.polar.circle(r=i+1, phi=list(range(0, 360, da)), phi_units="deg")
        c1 = plot.circle(x=xs, y=ys, visible=False)
        c1.group = group
    return plot, groups