I wanted to know if the procedure for regularization of the Coulomb potential outlined in Celletti (2003): Basics of regularization theory could be generalized to arbitrary polynomial potentials. So instead of starting from an equation of motion with a force $F \sim x^{-2}$ I started with a generalized force $F \sim x^{n}$. So the equation of motion and the corresponding energy equation are: $$ \begin{align} \ddot{x} + K x^n = 0 && \frac{1}{2}\dot{x} + \frac{K}{n+1}x^{n+1} = \text{const.} := - E \end{align} $$ Now following Celletti I looked at the generalized transformation $$ \begin{align} x = u^{-n} && \frac{dt}{ds} = x = u^{-n} \end{align} $$ Inserting these transformations into the equation of motion and the energy equation yields with $\frac{du}{ds}=u'$ $$ \begin{align} u'' - u \left[ \left(\frac{u'}{u}\right)^2 + \frac{K}{n}u^{n(n+1)}\right] = 0 && \frac{1}{2}\left(\frac{u'}{u}\right)^2 + \frac{K}{n^2(n+1)}u^{n(n+1)} = -E \end{align} $$ Now the whole exercise of regularization is to transform the equation of motion to the simpler form of an harmonic oscillator $u'' + 2Eu = 0$ (Bartsch(2003): The Kustaanheimo-Stiefel transformation in geometric algebra). But we see from the two equations above, that the terms in the square bracket are only equal to the twice the energy, if $$ n = \frac{n^2(n+1)}{2} $$ or after solving the corresponding quadratic equation $n^2+n-2=0$ $$ \begin{align} n_1 = 1 \ &\text{and}\ n_2 = -2 \end{align} $$ So it seems that one can only regularize Coulomb forces $F \sim x^{-2}$ and harmonic forces $F \sim x$. This result intrigued me a bit, because it singles out the same forces/potentials as Bertrand's theorem. But what does the the property that all bound orbits are also closed orbits (Bertrand's theorem) have to do with regularization? What is so special about the Coulomb and the harmonic potential that they get singled out in both cases? Do these potentials have some common deeper symmetry, so that regularization and Bertrand's theorem can only be fulfilled by them? Or is this just a spurious relationship between Bertrand's theorem and regularization theory?

I wanted to know if the procedure for regularization of the Coulomb potential outlined in Celletti (2003): Basics of regularization theory could be generalized to arbitrary polynomial potentials. So instead of starting from an equation of motion with a force $F \sim x^{-2}$ I started with a generalized force $F \sim x^{n}$. So the equation of motion and the corresponding energy equation are: $$ \begin{align} \ddot{x} + K x^n = 0 && \frac{1}{2}\dot{x}^2 + \frac{K}{n+1}x^{n+1} = \text{const.} := - E \end{align} $$ Now following Celletti I looked at the generalized transformation $$ \begin{align} x = u^{-n} && \frac{dt}{ds} = x = u^{-n} \end{align} $$ Inserting these transformations into the equation of motion and the energy equation yields with $\frac{du}{ds}=u'$ $$ \begin{align} u'' - u \left[ \left(\frac{u'}{u}\right)^2 + \frac{K}{n}u^{n(n+1)}\right] = 0 && \frac{1}{2}\left(\frac{u'}{u}\right)^2 + \frac{K}{n^2(n+1)}u^{n(n+1)} = -E \end{align} $$ Now the whole exercise of regularization is to transform the equation of motion to the simpler form of an harmonic oscillator $u'' + 2Eu = 0$ (Bartsch(2003): The Kustaanheimo-Stiefel transformation in geometric algebra). But we see from the two equations above, that the terms in the square bracket are only equal to the twice the energy, if $$ n = \frac{n^2(n+1)}{2} $$ or after solving the corresponding quadratic equation $n^2+n-2=0$ $$ \begin{align} n_1 = 1 \ &\text{and}\ n_2 = -2 \end{align} $$ So it seems that one can only regularize Coulomb forces $F \sim x^{-2}$ and harmonic forces $F \sim x$. This result intrigued me a bit, because it singles out the same forces/potentials as Bertrand's theorem. But what does the the property that all bound orbits are also closed orbits (Bertrand's theorem) have to do with regularization? What is so special about the Coulomb and the harmonic potential that they get singled out in both cases? Do these potentials have some common deeper symmetry, so that regularization and Bertrand's theorem can only be fulfilled by them? Or is this just a spurious relationship between Bertrand's theorem and regularization theory?

I wanted to know if the procedure for regularization of the Coulomb potential outlined in Celletti (2003): Basics of regularization theory could be generalized to arbitrary polynomial potentials. So instead of starting from an equation of motion with a force $F \sim x^{-2}$ I started with a generalized force $F \sim x^{n}$. So the equation of motion and the corresponding energy equation are: $$ \begin{align} \ddot{x} + K x^n = 0 && \frac{1}{2}\dot{x} + \frac{K}{n+1}x^{n+1} = \text{const.} := - E \end{align} $$ Now following Celletti I looked at the generalized transformation $$ \begin{align} x = u^{-n} && \frac{dt}{ds} = x = u^{-n} \end{align} $$ Inserting these transformations into the equation of motion and the energy equation yields with $\frac{du}{ds}=u'$ $$ \begin{align} u'' - u \left[ \left(\frac{u'}{u}\right)^2 + \frac{K}{n}u^{n(n+1)}\right] = 0 && \frac{1}{2}\left(\frac{u'}{u}\right)^2 + \frac{K}{n^2(n+1)}u^{n(n+1)} = -E \end{align} $$ Now the whole exercise of regularization is to transform the equation of motion to the simpler form of an harmonic oszillator $u'' + 2Eu = 0$ (Bartsch(2003): The Kustaanheimo-Stiefel transformation in geometric algebra). But we see from the two equations above, that the terms in the square bracket are only equal to the twice the energy, if $$ n = \frac{n^2(n+1)}{2} $$ or after solving the corresponding quadratic equation $n^2+n-2=0$ $$ \begin{align} n_1 = 1 \ &\text{and}\ n_2 = -2 \end{align} $$ So it seems that one can only regularize Coulomb forces $F \sim x^{-2}$ and harmonic forces $F \sim x$. This result intrigued me a bit, because it singles out the same forces/potentials as Bertrands theorem. But what does the the property that all bound orbits are also closed orbits (Bertrands theorem) have to do with regularization? What is so special about the Coulomb and the harmonic potential that they get singled out in both cases? Do these potentials have some common deeper symmetry, so that regularization and Bertrands theorem can only be fulfilled by them? Or is this just a spurious relationship between Bertrands theorem and regularization theory?

I wanted to know if the procedure for regularization of the Coulomb potential outlined in Celletti (2003): Basics of regularization theory could be generalized to arbitrary polynomial potentials. So instead of starting from an equation of motion with a force $F \sim x^{-2}$ I started with a generalized force $F \sim x^{n}$. So the equation of motion and the corresponding energy equation are: $$ \begin{align} \ddot{x} + K x^n = 0 && \frac{1}{2}\dot{x} + \frac{K}{n+1}x^{n+1} = \text{const.} := - E \end{align} $$ Now following Celletti I looked at the generalized transformation $$ \begin{align} x = u^{-n} && \frac{dt}{ds} = x = u^{-n} \end{align} $$ Inserting these transformations into the equation of motion and the energy equation yields with $\frac{du}{ds}=u'$ $$ \begin{align} u'' - u \left[ \left(\frac{u'}{u}\right)^2 + \frac{K}{n}u^{n(n+1)}\right] = 0 && \frac{1}{2}\left(\frac{u'}{u}\right)^2 + \frac{K}{n^2(n+1)}u^{n(n+1)} = -E \end{align} $$ Now the whole exercise of regularization is to transform the equation of motion to the simpler form of an harmonic oscillator $u'' + 2Eu = 0$ (Bartsch(2003): The Kustaanheimo-Stiefel transformation in geometric algebra). But we see from the two equations above, that the terms in the square bracket are only equal to the twice the energy, if $$ n = \frac{n^2(n+1)}{2} $$ or after solving the corresponding quadratic equation $n^2+n-2=0$ $$ \begin{align} n_1 = 1 \ &\text{and}\ n_2 = -2 \end{align} $$ So it seems that one can only regularize Coulomb forces $F \sim x^{-2}$ and harmonic forces $F \sim x$. This result intrigued me a bit, because it singles out the same forces/potentials as Bertrand's theorem. But what does the the property that all bound orbits are also closed orbits (Bertrand's theorem) have to do with regularization? What is so special about the Coulomb and the harmonic potential that they get singled out in both cases? Do these potentials have some common deeper symmetry, so that regularization and Bertrand's theorem can only be fulfilled by them? Or is this just a spurious relationship between Bertrand's theorem and regularization theory?