switch from depth first to nx bfs algorithm

[M-R-Schaefer]

The home made depth first search in the unwrapping code could get stuck in infinite recursion for pathological structures such as fullerenes.