Hi @junpenglao,
Thanks for getting back to me so quickly.
Yes, it samples fine, although, because of the lack of constraints, the gss variables are not informative (or at least that’s my intuition)
p.s. i also tried using the pm.Potential() to set the order constraints, but this resulted in a bad inital energy error.