Let $F$ be a non-singular integral quadratic form in four variables. Then a result of Heath-Brown from the 90's states for $m \to \infty$, $$|\{ x \in \mathbb{Z}^4 \,:\, F(x) = m \}| = C_F\sigma(F,m)m + O_{F,\varepsilon}(m^{\frac{3}{4} + \varepsilon}),$$ where $C_F$ is a positive constant depending only on $F$ and $\sigma(F,m)$ is the singular series.
I am interested to know, whether up till today there are improvements on the error rate of this result. More precisely, is there any hope to improve the error rate to $O_{F,\varepsilon}(m^{\frac{1}{2} + \varepsilon})$. If this is unknown, what is the conjectured error term?
For the application I have in mind, it would be already interesting to understand the case $$F(x) = x_1^2 + x_2^2 + x_3^2 + x_4^2.$$ Are there any results for this case or any other specific quadratic form?