Transform your 0-1 variable so that it is asymptotically -Inf to +Inf, then do an unconstrained smooth, then invert the transformation.
Logit, or probit or something should help. Apologies for incompleteness, but you didn't supply a reproducible example and I've got things to do.