The phase space between the reacting centers flattens due to the annihilation of the (m, s2a) pair through a fold upon further stretching of the B-N bond (Figure Panels c and d). Detailed information concerning the B-N critical values can be found in Ref. 29. Finally, two pairs of saddles coalesce not simultaneously; see Figure , Panels d-f. Regardless that any topographical change resulting from merging/arising saddles cannot be describe via any unfolding, the |det H| at all these CPs is provided for illustrating the assigning procedure; see Figure Panel g. Relying on Thom's catastrophe theory, no umbilic fingerprint is detected, as no critical point changes its index. This indicates that a function with one parameter and one state variable is the correct polynomial for characterizing this reaction from a topological standpoint. This is contrary to the report by Silvi and collaborators. Additionally, the number of eigenvalues of |det H| tending to zero is a critical variable often overlooked when classifying relevant chemical processes along a reaction pathway via unfoldings. From a pure topological perspective, distinguishing between unfoldings with one and two state variables can only be done considering the number of null eigenvalues. Hence, we propose that the BH3NH3 dissociation has to be characterized via a cuspoid (i.e., a one-state-variable function) and not in terms of an umbilic (i.e., a two-state-variable function), as only one eigenvalue of the CPs involved in the fold catastrophe tends to zero value. For this reaction, it is appropriate to consider the fold parameter controlling the changes in the ELF topography as the negative of the B-N enlargement since catastrophic bifurcations result from changes in this bond length.