Computational modeling of non-adiabatic chemistry at the conical intersection seam in polyatomic molecules