What the standard error of estimated regression coefficients tell about their corresponding predictor variables? In linear regression modelling, we obtain the standard error of estimated coefficients, and we use them to draw inferences about the corresponding population-level coefficients, which are unknown. Besides this important role of standard errors, I was wondering what information they convey about the predictor's variability and its explanatory power. Although we assume predictor variables are fixed (not random), they practically exhibit variability, and sometimes we tend to prefer (for variable selection) more variable predictors than those with low variability. For this reason, I ran a simple simulation to correlate the variability of predictors and the standard error of their estimated coefficient. Please find the codes below to reproduce my results.
set(1345)
n <- 500 # sample size
p <- 30  # number of predictors (all continuous)
range.x <- 1:p  # The range of variables in uniform distribution


# simulating predictors
X <- sapply(1:p, function(j){
  runif(n, 0, j)
})

bet <- c(0, runif(p, 0.25, 0.3)) # true regression coefficients

# simulate the response (you know what I'm doing)
Y <- cbind(1, X) %*% as.matrix(bet) + rnorm(n, 0, 1)

fit.lm <- lm(Y~X) 
se.coef <- diag(vcov(fit.lm))[-1]


plot(range.x, se.coef, pch=19) 

In the simulation procedure, I simulated each predictor from a uniform distribution between 0 and p, where p id the index predictor p (p=,1, 2, ..., 30). That means each predictor has a different range, hence different variability between predictors. The rest of the simulation is very simple and does not require explanation. If not, please let me know.

The relationship looks fascinating, i.e. lower standard error means large variability of the predictor. And I wonder if I can interpret it as such? That is, a low standard error of estimated coefficients translates to a high explanatory power of the corresponding predictor for the response variable.

This is an extension of what is seen in single-predictor OLS. Under the assumption of normal error terms:

the estimator of the slope coefficient will itself be normally distributed with mean $\beta$ and variance $\sigma^2/\sum(x_i-\bar x)^2$, where $\sigma^2$ is the variance of the error terms

Other things being equal, a higher variance in $x$ values gives more precise estimates of the slope $\beta$. This makes some intuitive sense, as if you only have a narrow distribution of $x$ values you have a correspondingly narrow distribution of $y$ values and thus less information to get a good estimate of the slope, their ratio after subtracting their mean values.
