In Mathematical Methods of Classical Mechanics Arnold defines a central field as a vector field $\mathbf{F}\colon \mathbb{R}^n \to \mathbb{R}^n$ which is invariant under the group of galilean transformations which fix the origin (if I understand him correctly). From this he concludes that $\mathbf{F}(\mathbf{r}) = \Phi(r) \mathbf{e}_r$, where $r = \lVert \mathbf{r} \rVert,$ for some scalar function $\Phi\colon \mathbb{R}_{\geq 0} \to \mathbb{R}$.
This is fine. Next, he shows that any central field is conservative by showing that the line integral between two points $P_1$ and $P_2$ is independent of path. To this end he simply writes:
$$ \int_{P_1}^{P_2} \mathbf{F}(\mathbf{r}) \cdot \mathrm{d}\mathbf{S} = \int_{r(P_1)}^{r(P_2)} \Phi(r) \, \mathrm{d}r. $$
(Presumably $r(P)$ is the radial coordinate of the point $P$.) But I have a problem with this equality. Sure, it "follows" easily from calculating $\mathbf{F} \cdot \mathrm{d}\mathbf{S}$, but I have never seen any justification for this other than the usual physical one. I was wondering if there is any way to write out the definition of the line integral and show the above equality that way? I have tried choosing an arbitrary path and parametrisation, but I don't end up with anything constructive.
Incidentally, if there is any rigorous justification for simply calculating $\mathbf{F} \cdot \mathrm{d}\mathbf{S}$, I would love to hear that as well.
Not sure if can get to the core of your question, but I'll place a few key points here:
For any vector field $f:\mathbb{R}^3 \rightarrow \mathbb{R}^3$, the integral over a curve $\gamma$ is denoted by: $$ \int_\gamma f(\gamma)d\gamma $$ Which is something I also dislike.
Assuming it is a Parametrized curve $\gamma:\mathbb{R}\rightarrow \mathbb{R}^3 $ , then: $$ \gamma=\gamma(t) $$ $$ d\gamma =\gamma'(t)dt $$
Such that the previous integral becomes: $$ \int_\gamma f(\gamma)d\gamma = \int_{\gamma (t_1)}^{\gamma (t_2)} f(\gamma(t)). \gamma'(t) dt $$
Where the dot product was made explicit. Noting that: $$ g(t) = f(\gamma(t)). \gamma'(t) $$
We have $g:\mathbb{R}\rightarrow \mathbb{R}$. Assume there exists $F:\mathbb{R}^3\rightarrow \mathbb{R}$ such that: $$ \nabla F(x,y,z) = f(x,y,z) $$ Then: $$ \frac{dF(\gamma(t))}{dt} = \nabla F(\gamma(t)) . \gamma '(t) = f(\gamma(t)) . \gamma '(t) $$
So if there exists $G(t)$ such that $G'(t)=g(t)$, then $G(t) = F(\gamma(t)) $:
$$ \int_\gamma f(\gamma)d\gamma = \int_{\gamma (t_1)}^{\gamma (t_2)} f(\gamma(t)). \gamma'(t) dt = G(t_2)-G(t_1) = F(\gamma(t_2))-F(\gamma(t_1)) $$
So the line integral only depends on the points $\gamma (t_1)$ and $\gamma (t_2)$, not on which curve $\gamma$ was chosen, hence the path independence.
In the case the book presents $\gamma(r)=r\textbf{e}_r$ such that $\gamma'(r)=r'\textbf{e}_r+r\textbf{e}'_r=\textbf{e}_r+r\textbf{e}'_r$ ($r'=1$), thus:
$$ F(\textbf{r}) d\textbf{S} = (\Phi(r)\textbf{e}_r).(\textbf{e}_r+r\textbf{e}'_r) dr $$
But since $\textbf{e}_r$ is a unit vector, $\textbf{e}_r . \textbf{e}'_r = 0$. Therefore:
$$ F(\textbf{r}) d\textbf{S} = \Phi(r) dr $$
Does this help?