In principle I don't see why you couldn't estimate these models in mirt, though for now you'd have to switch to the Rsolnp optimizer to account for the reverse-sign equality constraints in the transformed loadings matrix (what they call $\breve{\Lambda}$). There may be a problem with high-dimensionality as well, but I'm not familiar enough with the models to confirm that.
Otherwise, if you wanted to use the default optimizers for unbounded parameters, I would have to add a 'signed' argument somewhere to indicate that there are equality constraints but that their signs are flipped (i.e., instead of a1_1 = a1_2, it would have to be a1_1 = -a1_2), or more generally a transformation matrix for the entire parameter set (which I already use, but it's generally built up from inputs rather than directly accessible by the user). The latter would be the most general, so I'll probably consider adding that sometime in the future anyway.
Phil