Modelling constant phase elements with ModellingToolkit

Yes, it’s not possible at present. Right now, directly use DelayDiffEq.jl