It was already discussed in How does ModelingToolkit use bounds metadata? - #6 by ChrisRackauckas
And I opened an issue based on @ChrisRackauckas reply here Consider supporting bounds for `NonlinearProblem` resolution · Issue #183 · SciML/NonlinearSolve.jl · GitHub