The problem of d-dimensional Schrödinger equations with a position-dependent mass is analyzed in the framework of first-order intertwining operators. With the pair (H,H1) of intertwined Hamiltonians one can associate another pair of secondorder partial differential operators (R,R1), related to the same intertwining operator and such that H (resp. H1) commutes with R (resp. R1). This property is...