Fitting several regression models by changing only one independent variable within mutate()

Question

I suspect that this question might be a duplicate, however, I found nothing satisfactory. Imagine a simple dataset with a structure like this:

set.seed(123)
d