In this work, the CVHD method is implemented with Plumed, 32,33 which serves as a plug-in in our MD simulation of reaction systems. While biasing CVs that monotonously increase with distance quite naturally favors bond dissociations, bond formation processes are conceptually more dicult to accelerate using CVHD. A convenient solution to also bias bond formation in addition to bond dissociation is provided by the use of a switching function (SWF), 44 enabled by the CVHD implementation in Plumed. 33 A rational SWF is described in Eq. 2. It maps bond distances r from the interval of ]0, +∞[ to a variable s between ]0, 1[.