lme4: what is gained by scaling variables? I am performing a multilevel logistic regression with glmer from the lme4-package. Currently, it's a very simple model of the form: glmer(y ~ x + (1 | group), data = data, family = binomial). Running this in R returns the model fit accompanied by a warning message: "Model is nearly unidentifiable: very large eigenvalue". 
Therefore, I centered the variable x by using scale(x), as was mentioned for example here. This indeed lets the warning message disappear. Hooray! 
However, the output of the model is actually almost the same as before. The only things that changed are the estimate of the coefficient for y and its standard error, which is not surprising at all, the intercept and the correlation between the intercept and y. Everything else is the same as before scaling, e.g. the estimate of the random effects, the p-value of the coefficient of y, the AIC, BIC, deviance.
So, basically, all parameters that I am interested in are the same, regardless whether I center my variables or not. One thing that changes for sure, is the coefficient of my centered variable. However, this changes not in a good way, since I can't properly interpret the coefficient of a centered variable (or at least not as properly as the coefficient of the original variable).
Thus, my questions are: 


*

*Is my first model with the non-centered variable necessarily problematic? 

*If the answer to the first questions is "yes", is my second model with the centered variable really non-problematic, since it gives essentially the same output as the other model?

*Should I do more to improve my models?

 A: You mention that you cannot properly interpret the coefficient of a standardized variable, and I would disagree. Imagine you have a variable that is scaled from 0 to 100. A 1-unit difference in such a variable might not represent a very meaningful quantity. By using scale(x) you standardize that variable relative to a normal distribution - from each person's/observation's value of $x$, you subtract $\bar{x}$ and divide by sd($x$). Recall that regression coefficients are interpreted as the association between a 1-unit increase/decrease in x and the outcome (here the coefficient is either in log odds units or odds ratios). By standardizing x, a 1-unit increase in x is equivalent to 1 standard deviation increase in x. A lot of people like that interpretation. 
Regarding your question:

Is my first model with the non-centered variable necessarily
  problematic?

Only insofar is that you are getting warnings about reaching the limits of identification for the model. This is a somewhat common problem people encounter with glmer model fits. See this CV post for some guidelines. The linked thread deals with your other questions as well. Bottom, line you can increase the number of iterations, try using different optimizers, and check for singularity. 
You can also try fitting your model with a different package, such as GLMMadaptive to check if the estimates are similar and if you get similar warnings.
