---
abstract: 'We modify Hamiltonian mechanics. We reformulate the law of conservation of energy.'
author:
- |
    S.E.Akrami[^1]\
    Departments of Mathematics and Physics,\
    Institute for Research in Fundamental Sciences,\
    P.O.Box:19395-5746, Tehran, Iran
title: Topics on Hamiltonian Mechanics
---

\[thm\][Lemma]{} \[thm\][Proposition]{} \[thm\][Corollary]{} \[thm\][Definition]{} \[thm\][Example]{}

23.6cm -.2in

**Acknowledgement** I thank Allah Rabbel-Alamin and Imam Zaman.

This paper is based on the following observation. Consider the following modified Hamiltonian equations [$$\frac{{\partial}Q}{{\partial}t_1}+\frac{{\partial}Q}{{\partial}t_2}=\frac{{\partial}H}{{\partial}P},~~~\frac{{\partial}P}{{\partial}t_1}+\frac{{\partial}P}{{\partial}t_2}=-\frac{{\partial}H}{{\partial}Q}\label{c-hamilton}$$]{} for unknown functions $Q=Q(t_1,t_2),P=P(t_1,t_2)$ and given Hamiltonian $H=H(Q,P)$.

If we set $q(t):=Q(t,t) ,p(t):=P(t,t)$ then since $\dot{q}(t)=\frac{{\partial}Q}{{\partial}t_1}(t,t)+\frac{{\partial}Q}{{\partial}t_2}(t,t)$ and $\dot{p}(t)=\frac{{\partial}P}{{\partial}t_1}(t,t)+\frac{{\partial}P}{{\partial}t_2}(t,t)$ we conclude that $q(t)$ and $p(t)$ satisfy the classical Hamilton’s equations. Thus if we computes the energy $H(q(t),p(t))$ then we get a constant $H(q(t),p(t))=H(q(0),p(0))$.

But instead of the above method, i.e. before putting $t_1=t_2=t$ if we expand $Q$ and $P$ in terms of powers of $t_2$ as [$$Q(t_1,t_2)=\sum_{n=0}^\infty
q_n(t_1)t_2^n,P(t_1,t_2)=\sum_{n=0}^\infty p_n(t_1)t_2^n,\label{powers}$$]{} and then substitute them in the equations (\[c-hamilton\]) and then by equating the coefficients of powers of $t_2$ in two sides of each equation, we get some recursive relations among the unknown coefficients $q_n(t_1),p_n(t_1)$ which can be solved by knowing the initial conditions $q_0(t_1),p_0(t_1)$ and then substituting these coefficients in (\[powers\]) we obtain $Q=Q(t_1,t_2)$ and $P=P(t_1,t_2)$ and at last by setting $t_1=t_2=t$ we obtain $q(t):=Q(t,t) ,p(t):=P(t,t)$. In all examples which we were able to solve the equations by this method, we observed that after calculating the energy $H(q(t),p(t))$ which as mentioned above that in general one expects to get a constant, in fact we got

[$$H(q(t),p(t))=\hat{h}(t),\label{qconservationlaw}$$]{} where $h(t):=H(q_0(t),p_0(t))$ and [$$\hat{h}(t):=\sum_{n=0}^{\infty}(-1)^n\frac{h^{(n)}(t)}{n!}t^n.\label{alter}$$]{}By $h^{(n)}$ we mean the $n$-th derivative of $h$.

Since we expect constant energy, the question arises that if for any function $f(t)$ one has $\hat{f}(t)=f(0)$ for all $t$? In fact if we differentiate the series term by term we get $$\begin{aligned}
\frac{d}{dt}\hat{f}(t)&=&
\sum_{n=0}^{\infty}(-1)^n\frac{f^{(n+1)}(t)}{n!}t^n+
\sum_{n=1}^{\infty}(-1)^n\frac{f^{(n)}(t)}{(n-1)!}t^{n-1}\nonumber\\&=&
\sum_{n=0}^{\infty}(-1)^n\frac{f^{(n+1)}(t)}{n!}t^n-\sum_{n=0}^{\infty}(-1)^n\frac{f^{(n+1)}(t)}{n!}t^n
\nonumber\\&=&0\nonumber\end{aligned}$$ Thus $\hat{f}(t)$ should be constant. But the point is that we are not allowed to differentiate term by term from a series even the series is uniformly convergent.

We were not able to prove this up to the present time. The only result is

If a function $f$ can be expanded around origin via power series $f(t)=\sum_{n=0}^{\infty}\frac{f^{(n)}(0)}{n!}t^n,t\in{\Bbb R}$, then $\hat{f}$ is absolutely and uniformly convergent to the constant $f(0)$. That is [$$\sum_{n=0}^{\infty}(-1)^n\frac{f^{(n)}(t)}{n!}t^n=f(0),~~~t\in{\Bbb R},\label{hat0}$$]{} for analytic functions.

Another evidence for strangeness of the series $\hat{f}(t)$ comes from the following argument. Suppose that $f$ is a smooth periodic function which by the Fourier analysis we know that it has Fourier expansion. That is we suppose $f(t)=\sum_{m=-\infty}^\infty
c_me^{im{\omega}t}.$ Then it is well known that we can differentiate to get $f^{(n)}(t)=\sum(im{\omega})^n c_me^{im{\omega}t}.$ Thus $$\begin{aligned}
\hat{f}(t)&=&\sum_n\sum_m\frac{(-1)^n(im{\omega})^{n}}{n!}c_me^{im{\omega}t}t^n\nonumber\\&=&
\sum_m\sum_n(\frac{(-1)^n(im{\omega})^{n}}{n!}t^n)c_me^{im{\omega}t}\nonumber\\&=& \sum_me^{-im{\omega}t}c_me^{im{\omega}t}\nonumber\\&=&\sum_{m=-\infty}^{\infty}c_m\nonumber\\&=&f(0).\nonumber\end{aligned}$$ But the above argument is again analytically ill, since we have exchanged the order of two infinite sums which from the theorems of mathematical analysis we are not allowed in general to do so. In order to be able to do so there exists a general theorem on double series which states that if for a double infinite series $\sum_{m,n}a_{mn}$ for each $n$ the series $\sum_m |a_{mn}|$ is convergent which we show its sum by $b_n$ and the series $\sum_n
b_n$ is also convergent then we are allowed to exchange the order of summation. That is we have $\sum_n\sum_ma_{mn}=\sum_m\sum_na_{mn}$. Now let us check if this criterion can be applied to the double series whose entries are $a_{mn}:=\frac{(-im)^n}{n!}c_me^{im t}t^n$. For simplicity we have assumed that ${\omega}=1$. We have $\sum_m
|a_{mn}|=2\frac{t^n}{n!}\sum_{m=1}^\infty m^n|c_m|$. But in Fourier analysis it is well known that the series $f^{(n)}(t)=\sum(im)^n
c_me^{im t}$ is absolutely convergent. That is the series $\sum_{m=1}^\infty m^n|c_m|$ is convergent which we show its sum by ${\alpha}_n$. Thus we have $\sum_m |a_{mn}|=2\frac{{\alpha}_nt^n}{n!}.$ Thus we should verify the convergence of the series $\sum_{n=0}^\infty\frac{{\alpha}_nt^n}{n!}$. But this is a power series whose convergence radius is given by $R={\alpha}^{-1}$ where ${\alpha}:=\limsup\sqrt[n]{\frac{{\alpha}_n}{n!}}=\limsup\sqrt[n]{\frac{{\alpha}_n}{n!}}$. The point is that we are not sure if $R\ne0$?

In summary the question of convergence of the series $\hat{f}$ and as well as the constancy of the sum of the series for a periodic or a general smooth function is remain open.

**Open Questions** Is it true that for all systems (\[qconservationlaw\]) holds?

Is there non-analytic function $f$ such that the series $\hat{f}$ is point-wise or uniformly convergent and among such functions if there is any, is there function such that the series converges to a non-constant function?

We conjecture that there are nonanalytic probably periodic functions $f$ such that the series $\hat{f}(t)$ is piecewise constant and takes discrete values. That is $\hat{f}(t)$ is constant only on some intervals and its values jump when one passes from an interval to its next interval. More details will appear soon in arXive.

Thus if $q_0(t),p_0(t)$ are simultaneously analytic, that is if both can be expanded as power series $$q_0(t)=\sum_{n=0}^{\infty}q_{0n}t^n,~~~p_0(t)=\sum_{n=0}^{\infty}p_{0n}t^n$$ then the function $h(t):=H(q_0(t),p_0(t))$ is also analytic and therefore we recover the classical conservation law $H(q(t),p(t))=h(0)=H(q_0(0),p_0(0))=H(q_{00},p_{00})$. Thus we observe that only the coefficients $q_{00}$ and $p_{00}$ of the expansion of $q_0(t)$ and $p_0(t)$ contribute in the energy. That is in the case in which $q_0(t),p_0(t)$ are simultaneously analytic, dependence of $q_0(t),p_0(t)$ to time does not effect the system and we can safely, as in the classical mechanics, assume that $q_0(t)=q_0,p_0(t)=p_0$ are constants. That is from the beginning we did not need two dimensional time $(t_1,t_2)$ and the modified equations (\[c-hamilton\]) are superfluous and the classical Hamilton’s equations are enough. But remember that we were able only prove the constancy of $\hat{h}(t)$ only for simultaneously analytic primitive conditions $q_0(t),p_0(t)$. Thus if one of them is not analytic then we may find new feature in our model. So we need to study more such a strange object as (\[alter\])!

The following proposition will be useful in computations.

(i)If the series $\hat{f}$ and $\hat{g}$ are convergent then the series $\widehat{af+bg}$ are also convergent and $\widehat{af+bg}=a\hat{f}+b\hat{g}$

\(ii) If moreover at least one of series $\hat{f}$ and $\hat{g}$ are absolutely convergent then the series $\widehat{fg}$ is also convergent and $\widehat{fg}=\hat{f}\hat{g}$.

Examples
========

$H(q,p)=ap+bq$, where $a$ and $b$ are some given constants.

[[**Proof **]{}]{}The recursive equations obtained from (\[c-hamilton\]) after expansion (\[powers\]) are $(n+1)q_{n+1}+\dot{q}_n=a,(n+1)p_{n+1}+\dot{p}_n=-b.$ One can easily show by induction that $q_n(t_1)=\frac{(-1)^n}{n!}q^{(n)}_0(t_1),
~~~p_n(t_1)=\frac{(-1)^n}{n!}p^{(n)}_0(t_1)$ for $n\ne1$ and $q_1(t_1)=a-\dot{q}_0(t_1),p_1(t_1)=\dot{p}_0(t_1)-b.$ Thus $q(t)=Q(t,t)=\sum_{n=0}^\infty q_n(t)t^n=at+\widehat{q_0}(t)$ and $p(t)=-bt+\widehat{p_0}(t)$. Finally $H(q(t),p(t))=-abt+a\widehat{p_0}(t)+bat+b\widehat{q_0}(t)=\widehat{ap_0+bq_0}(t)=\hat{h}(t)$ where $h(t)=ap_0(t)+bq_0(t)=H(q_0(t),p_0(t)).$ [$\blacksquare$]{}

The harmonic oscillator $H(q,p)=\frac{1}{2}p^2+\frac{{\omega}^2}{2}q^2$.

[[**Proof **]{}]{}The recursive equations obtained from (\[c-hamilton\]) after expansion (\[powers\]) are [$$(n+1)q_{n+1}+\dot{q}_n=p_n,~~~(n+1)p_{n+1}+\dot{p}_n=-{\omega}^2
q_n.\label{h1}$$]{} We show by induction that for $n>0$ [$$q_n=\frac{(-1)^n}{n!}(\sum_{i=0}^{[\frac{n}{2}]}(-1)^iC^n_{2i}{\omega}^{2i}q_0^{(n-2i)}-\sum_{i=0}^{
[\frac{n-1}{2}]}(-1)^{i}C^n_{2i+1}{\omega}^{2i}p_0^{(n-2i-1)})\label{h3}$$]{}and [$$p_n=\frac{(-1)^n}{n!}(\sum_{i=0}^{[\frac{n}{2}]}(-1)^iC^n_{2i}{\omega}^{2i}p_0^{(n-2i)}+{\omega}^{2}\sum_{i=0}^{
[\frac{n-1}{2}]}(-1)^{i}C^n_{2i+1}{\omega}^{2i}q_0^{(n-2i-1)})\label{h4}$$]{} where $C_i^j:=\frac{j!}{i!(j-i)!}$.

First step of induction $n=1$. We put $n=0$ in (\[h1\]) to get $q_1=-\dot{q}_0+p_0$ and $p_1=-\dot{p}_0-{\omega}^2q_0$. On the other hand if we put $n=1$ in the right hand sides of (\[h2\]) we get $q_1=-(\dot{q}_0-p_0)$ and $p_1=-(\dot{p}_0+{\omega}^2 q_0)$. Thus the first step of induction is verified. Let (\[h1\]) is true for $n$ we prove it for $n+1$. First suppose $n=2k+1$ is odd. We have $$\begin{aligned}
-\frac{(n+1)!}{(-1)^{n+1}}q_{n+1}&=&\frac{n!}{(-1)^{n}}(-\dot{q}_n+p_n)\nonumber\\&=&
-\sum_{i=0}^{k}(-1)^iC^n_{2i}{\omega}^{2i}q_0^{(n-2i+1)}+\sum_{i=0}^{
k}(-1)^{i}C^n_{2i+1}{\omega}^{2i}p_0^{(n-2i)}\nonumber\\&+&
\sum_{i=0}^{k}(-1)^iC^n_{2i}{\omega}^{2i}p_0^{(n-2i)}+{\omega}^{2}\sum_{i=0}^{
k}(-1)^{i}C^n_{2i+1}{\omega}^{2i}q_0^{(n-2i-1)}\nonumber\\&=&
-\sum_{i=0}^{k}(-1)^iC^n_{2i}{\omega}^{2i}q_0^{(n-2i+1)}+{\omega}^{2}\sum_{i=1}^{
k+1}(-1)^{i-1}C^n_{2(i-1)+1}{\omega}^{2(i-1)}q_0^{(n-2(i-1)-1)}\nonumber\\&+&\sum_{i=0}^{
k}(-1)^{i}C^n_{2i+1}{\omega}^{2i}p_0^{(n-2i)}+
\sum_{i=0}^{k}(-1)^iC^n_{2i}{\omega}^{2i}p_0^{(n-2i)}\nonumber\\&=&
-q_0^{(n+1)}+(-1)^k{\omega}^{2k+2}q_0-\sum_{i=1}^{k}(-1)^i(C^n_{2i}+C^n_{2i-1}){\omega}^{2i}q_0^{(n-2i+1)}\nonumber\\&+&\sum_{i=0}^{
k}(-1)^{i}(C^n_{2i}+C^n_{2i+1}){\omega}^{2i}p_0^{(n-2i)}\nonumber\\&=&
-q_0^{(n+1)}+(-1)^k{\omega}^{2k+2}q_0-\sum_{i=1}^{k}(-1)^iC^{n+1}_{2i}{\omega}^{2i}q_0^{(n-2i+1)}\nonumber\\&+&\sum_{i=0}^{
k}(-1)^{i}C^{n+1}_{2i+1}{\omega}^{2i}p_0^{(n-2i)}\nonumber\\&=&
-\sum_{i=0}^{k+1}(-1)^iC^{n+1}_{2i}{\omega}^{2i}q_0^{(n-2i+1)}+\sum_{i=0}^{
k}(-1)^{i}C^{n+1}_{2i+1}{\omega}^{2i}p_0^{(n-2i)}\nonumber\end{aligned}$$ Thus $q_{n+1}=\frac{(-1)^{n+1}}{(n+1)!}(\sum_{i=0}^{[\frac{n+1}{2}]}(-1)^iC^{n+1}_{2i}{\omega}^{2i}q_0^{(n-2i+1)}-\sum_{i=0}^{
[\frac{n}{2}]}(-1)^{i}C^{n+1}_{2i+1}{\omega}^{2i}p_0^{(n-2i)})$. This means that the part (\[h3\]) of assertion of induction holds for $n+1$ if assertion of induction holds for odd $n$. Similarly one can show that the part (\[h4\]) of assertion of induction holds for $n+1$ if assertion of induction holds for odd $n$. The case of even $n$ is similar.

Moreover we show that if the primitive conditions $q_0(t)$ and $p_0(t)$ are such that the series $\widehat{q_0}(t)$ and $\widehat{q_0}(t)$ converges absolutely then [$$q(t)=\widehat{q_0}(t)\cos{\omega}t+{\omega}^{-1}\widehat{p_0}(t)\sin{\omega}t\label{h5}$$]{} and [$$p(t)=\widehat{p_0}(t)\cos{\omega}t-{\omega}\widehat{q_0}(t)\sin{\omega}t.\label{h6}$$]{}

Since the series $\cos({\omega}t)=\sum_{n=0}^\infty\frac{(-1)^n}{(2n)!}{\omega}^{2n}t^{2n}$ and $\sin({\omega}t)=\sum_{n=0}^\infty\frac{(-1)^n}{(2n+1)!}{\omega}^{2n+1}t^{2n+1}$ are absolutely convergent, by the Mertens’ theorem [@Apo], we can multiply the following series by the Cauchy’s product rule $$\begin{aligned}
\widehat{q_0}(t)\cos{\omega}t+{\omega}^{-1}\widehat{p_0}(t)\sin{\omega}t&=&
\sum_{n=0}^{\infty}(-1)^n\frac{q_0^{(n)}(t)}{n!}t^n\sum_{n=0}^\infty\frac{(-1)^n}{(2n)!}{\omega}^{2n}t^{2n}\nonumber\\&+&
{\omega}^{-1}\sum_{n=0}^{\infty}(-1)^n\frac{p_0^{(n)}(t)}{n!}t^n\sum_{n=0}^\infty\frac{(-1)^n}{(2n+1)!}{\omega}^{2n+1}t^{2n+1}\nonumber\\&=&
\sum_{n=0}^{\infty}\sum_{i=0}^{[\frac{n}{2}]}(\frac{(-1)^{n-2i}}{(n-2i)!}q_0^{(n-2i)}\frac{(-1)^{i}}{(2i)!}{\omega}^{2i})t^{n-2i}t^{2i}
\nonumber\\&+&{\omega}^{-1}\sum_{n=0}^{\infty}\sum_{i=0}^{[\frac{n-1}{2}]}(\frac{(-1)^{n-2i-1}}{(n-2i-1)!}p_0^{(n-2i-1)}\frac{(-1)^{i}}{(2i+1)!}{\omega}^{2i+1})t^{n-2i-1}t^{2i+1}
\nonumber\\&=&\sum_{n=0}^{\infty}(\frac{(-1)^n}{n!}\sum_{i=0}^{[\frac{n}{2}]}(-1)^iC^n_{2i}{\omega}^{2i}q_0^{(n-2i)}\nonumber\\&-&\sum_{i=0}^{
[\frac{n-1}{2}]}(-1)^{i}C^n_{2i+1}{\omega}^{2i}p_0^{(n-2i-1)}))t^n\nonumber\\&=&\sum_{n=0}^{\infty}q_n(t)t^n\nonumber\\&=&Q(t,t)
\nonumber\\&=&q(t)\nonumber\end{aligned}$$ This proves (\[h5\]). Similarly one can prove (\[h6\]).

Thus $H(q(t),p(t))=\frac{1}{2}\widehat{p_0}^2+\frac{{\omega}^2}{2}\widehat{q_0}^2=\hat{h}(t)$, where $h(t):=\frac{1}{2}p_0^2(t)+\frac{{\omega}^2}{2}q^2_0(t)=H(q_0(t),p_0(t)).$ [$\blacksquare$]{}

$H(q,p)=pq^2.$

[[**Proof **]{}]{}The recursive equations obtained from (\[c-hamilton\]) after expansion (\[powers\]) are [$$(n+1)q_{n+1}+\dot{q}_n=\sum_{i=0}^nq_iq_{n-i},~~~(n+1)p_{n+1}+\dot{p}_n=-2\sum_{i=0}^nq_ip_{n-i}.\label{}$$]{} Solution of these equations are difficult. So instead of solving them, we just show that the few first terms of the series $q(t)=Q(t,t)=\sum_{n=0}^{\infty}q_n(t)t^n$ and $\frac{\hat{q}_0(t)}{1-t\hat{q}_0(t)}$ coincide and also the few first terms of the series $p(t)=P(t,t)=\sum_{n=0}^{\infty}p_n(t)t^n$ and $\hat{p}_0(t)(1-t\hat{q}_0(t))^2$ coincide. That is we conjecture that [$$q(t)=\frac{\hat{q}_0(t)}{1-t\hat{q}_0(t)},~~~p(t)=\hat{p}_0(t)(1-t\hat{q}_0(t))^2.\label{x}$$]{} By solving the first few terms of these recursive relations we get $q_1=q_0^2-q_0',q_2=q_0^3,-2q_0q_0'+\frac{1}{2}q_0'',q_3=q_0^4-3q_0^2q_0'+q_0q_0''+q_0'^2-\frac{1}{6}q_0''',{\ldots}$

On the other hand if we set $\sum_{n=0}^\infty
a_nt^n:=\frac{\dot{q}_0(t)}{1-t\dot{q}_0(t)}$. Then $$\begin{aligned}
\sum_{n=0}^\infty\frac{(-1)^n}{n!}q_0^{(n)}t^{n}&=&
(1-\sum_{n=0}^\infty\frac{(-1)^n}{n!}q_0^{(n)}t^{n+1})\sum_{n=0}^\infty
a_nt^n\nonumber\\&=&a_0+(a_1-a_0q_0)t+(a_2-a_1q_0+a_0q_0')t^2\nonumber\\&+&(a_3-a_2q_0+a_1q_0'-\frac{1}{2}a_0q_0'')t^3+{\cdots}\nonumber\end{aligned}$$Thus $a_0=q_0,a_1-a_0q_0=-q_0',a_2-a_1q_0+a_0q_0'=\frac{1}{2}q_0'',a_3-a_2q_0+a_1q_0'-\frac{1}{2}a_0q_0''=\frac{-1}{6}q_0''',{\ldots}$. Thus $a_0=q_0,a_1=q_0^2-q_0',a_2=q_0^3,-2q_0q_0'+\frac{1}{2}q_0'',a_3=q_0^4-3q_0^2q_0'+q_0q_0''+q_0'^2-\frac{1}{6}q_0''',{\ldots}$.Hence we see that $a_0=q_0,a_1=q_1,a_2=q_2,a_3=q_3,{\ldots}$. This proves the first relation of (\[x\]). The other part is proved similarly.

Up to the present time we have not been able to prove the conjecture (\[x\]). But suppose this conjecture is true. Then we have $H(q(t),p(t))=(\dot{q}_0(t))^2\dot{p}_0(t)=\hat{h}(t)$ where $h(t)=(q_0(t))^2p_0(t)=H(q_0(t),p_0(t))$. [$\blacksquare$]{}

Tom M. Apostol *Mathematical Analysis*, Addison-Wesley,1975

[^1]: Supported financially by the grant 83810319 from Institute for Research in Fundamental Sciences, P.O.Box:19395-5746, Tehran, Iran. E-mail:akramisa@ipm.ir
