For these structures, the CMS-PDFT energies are non-unique, and the CMS-PDFT PESs can be discontinuous. A second disadvantage of CMS-PDFT is that the intermediate-state basis is obtained by solving a nonlinear system of equations, and there is no guarantee that this solution is uniquely defined for a given molecular geometry, atomic-orbital basis, active space, and model space. This additional non-convex global optimization step thus carries additional risk of not converging the iterations or converging to inequivalent local minima of the intermediate-state objective function. It is desirable to have a MS-PDFT method that yields a unique energy for any nuclear configuration, does not have an iterative procedure to find intermediate states, is able to handle a wide variety of systems, and preserves both spatial and spin symmetries. That is the goal of the work presented here.