Consider the equation $$\beta x\text{erf}(x)=\frac{1}{\sqrt{\pi}}\exp(-x^2),$$ where $\beta\ll1$ and $$\text{erf}(x)=\frac{2}{\sqrt{\pi}}\int_0^xe^{-t^2}\text{d}t$$ is the error function. I am looking for an approximate solution of the form $$x\approx x_1+\delta x_2,$$ where $\delta=\delta(\beta)$. For instance, for $\beta\gg1$, one can find the solutions $x\approx\pm(2\beta)^{-1/2}$, but I don't know where to start with the case where $\beta$ is very small.
Thanks in advance!