This research demonstrates a Contextual Subspace Variational Quantum Eigensolver on superconducting hardware that successfully models N₂ molecular bond-breaking with superior accuracy to single-reference wavefunction techniques, achieved through comprehensive error mitigation strategies and circuit optimization while requiring fewer quantum resources than comparable classical approaches.