I have to solve the following integral $$\Phi(\boldsymbol r) =- \int_{\boldsymbol r_0}^{\boldsymbol r}\boldsymbol{E}(\tilde{\boldsymbol r}) \cdot d\boldsymbol s$$ This comes from Electrostatics (I'm computing the electric potential of a charge Q at the origin, with reference point $\boldsymbol r$ being asymptotic infinity) We know that in general for a charge at the origin$$\boldsymbol E(\boldsymbol r) = \frac{1}{4\pi\epsilon_0}\frac{Q}{r^2}\hat{\boldsymbol e}_r$$
My trial
Clearly we use polar coordinates $\{r, \theta, \phi\}$ because of the spherical symmetry.Then the line element becomes $d\boldsymbol s = dr\hat{\boldsymbol e}_r+rd\theta\hat{\boldsymbol e}_{\theta}+r\sin{\theta}d\phi\hat{\boldsymbol e}_{\phi}$ so then $$-\int_{\boldsymbol r_0}^{\boldsymbol r}\boldsymbol{E}(\tilde{\boldsymbol r}) \cdot d\boldsymbol s = -\frac{Q}{4\pi\epsilon_0}\int_{\boldsymbol r_0}^{\boldsymbol r}\frac{1}{\tilde{r}^2}\hat{\boldsymbol e}_{\tilde r}\cdot d\boldsymbol s =-\frac{Q}{4\pi\epsilon_0}\int_{\boldsymbol r_0}^{\boldsymbol r}\frac{1}{\tilde{r}^2}d\tilde{r}$$
Because the spherical polar basis is an orthogonal basis, hence they are all perpendicular to each other and the result of the dot product is $1$ coming from the radial vector. But now I have a scalar integral and I will have to evaluate the result at a vector?
Can someone explain to me (possibly explaining each passage and with some strong mathematical rigor) how someone would finish this calculation?