Wha are $R^2$ values not available in fixed effects meta regressions

I'm currently conducting a meta-regression for my moderators from a meta-analysis my team recently conducted. In calculating the fixed effects meta-regressions, my objective is to identify how much heterogeneity is accounted for when the moderators are considered in the model. To do this, I used the rma() function from the metafor package in R:

Reg1 <- rma(yi = LogOdds, sei = SE, data = data, method = 'FE', mods = ~ mod1 + mod2)


The output provides out $$I^2$$ and $$H^2$$ values as well as the significance of our moderators, however, it does not provide an $$R^2$$ value for the amount of heterogeneity accounted for.

I did notice, however, that running an random effects meta regression dos provide the $$R^2$$.

My questions are: Does a meta-regression have to be random effects in order to garner the amount of heterogeneity accounted for? If I ran a fixed effects meta-analysis, must I also run a fixed effects meta-regression? And if not, how can I estimate the amount of heterogeneity accounted for in a fixed effects meta-regression?

The standard definition of the 'pseudo' $$R^2$$ statistic in meta-analysis is based on the proportional reduction in the estimate of $$\tau^2$$ (the amount of (residual) heterogeneity), that is: $$R^2 = \frac{\hat{\tau}_0^2 - \hat{\tau}_1^2}{\hat{\tau}_0^2},$$ where $$\hat{\tau}_0^2$$ is the estimate from a standard random-effects model (and hence this estimates the total amount of heterogeneity) and $$\hat{\tau}_1^2$$ is the estimate from the mixed-effects meta-regression model (and hence this estimates the residual amount of heterogeneity, that is, the amount not accounted for by the moderators included in the model).
For fixed-effects models, we do not estimate $$\tau^2$$ and hence this definition is no longer applicable. However, one could make the argument that an analogous version in this context is based on the standard definition of $$R^2$$ in regression models (although using weighted least squares estimation since we typically use inverse-variance weights in meta-analysis). I would suggest to compute the adjusted $$R^2$$ though, since $$k$$ is typically small in meta-analyses, so the adjustment is often relevant then.
Since you are using rma() from the metafor package: I have decided to implement this in the development version, so if you follow these instructions, you will now get the (adjusted) $$R^2$$ also for fixed-effects models.