Hi Ben,
I’m afraid my answer to that other question was misleading, because it results in estimating two parameters (in the current version of your model alpha_ar and d15n_ar), while what you want is to estimate a_r jointly from the two formulas.
I think this is currently not possible in brms but you could modify the generated Stan code to achieve it.