Further optimization and parallelization of our candidate searching algorithm will allow us to simulate larger systems over longer timescales. With the main limitation of this method being the timescales that can be accessed by MD, more sophisticated biasing techniques could be used to enhance dynamical sampling of bonding topologies and treatment of simulations in the rare-event regime.