To make this using Stella 4d (available here), I applied the “morph duals by expansion” function to the great icosahedron.