Realized it was due to plot scaling, not a programming error that the harmonic bonds did not appear correct