Polynomial supersymmetry for matrix Hamiltonians

Abstract We study intertwining relations for matrix non-Hermitian Hamiltonians by matrix differential operators of arbitrary order. It is established that for any intertwining operator of minimal order there is operator that intertwines the same Hamiltonians in the opposite direction and such that the products of these operators are identical polynomials of the corresponding Hamiltonians. The related polynomial algebra of supersymmetry is constructed. The problems of minimization and reducibility of a matrix intertwining operator are considered and the criteria of minimizability and reducibility are presented. It is shown that there are absolutely irreducible matrix intertwining operators, in contrast to the scalar case.