How to specify this potential?

Questions
2

How is your current implementation looks like?
In general, the logp function (that you would wrap in a DensityDist) is alpha * fn_s(T), bounded by the constrain.