Computational mechanochemistry studies use various techniques to deform systems of interest and explore the evolution of the system as an external force or displacement is increased. For example, SMD anchors atoms to an evolvable position, with restoring forces on the anchor site, using a simple potential, typically harmonic. CoGEF performs geometry optimizations while constraining the desired atoms. Evolving these constraints (e.g. increasing linear distance between two atoms) can simulate molecular deformation to incite mechanochemistry. Such approaches cannot readily capture complex many body molecular distortions observed in condensed matter under mechanical loads . To address this gap, we develop Many-Bodied Steered Molecular Dynamics (MB-SMD) in which we use four-body external biasing potential to deform molecular groups of interest along torsional and out-of-plane degrees of freedom. Four-body terms are motivated by our observation of the large molecular distortions observed under shock loading as well as their general applicability to explore the conformational space of a molecule. This external biasing potential is implemented in the form of four-body harmonic inter-atomic potential, 𝐸 = 𝐾(𝜙 𝑖𝑗𝑘𝑙 -𝜙 𝑜 ) 2