A charged relativistic free particle has the Hamiltonian in general:

$$ \mathcal{H} = \sqrt{{\bf p}^2c^2+m^2c^4}.$$

I read somewhere that says, it is possible to go further and say that the EoM are Hamilton's equations. But it is not done as there is "less interest" in such a discussion.

Is there something deeper to this? Like another formalism is ''better''.

(My guess is a more trivial one though. That is, it is not useful because the equations just get very cluttered and ugly)

  • 1
    $\begingroup$ I'm...not sure what your question is. The e.o.m exist, but they get quite ugly due to the square root. We usually introduce the einbein to get a more convenient Hamiltonian formalism. What is the question about that? $\endgroup$ – ACuriousMind Apr 20 '16 at 17:54
  • $\begingroup$ Related: physics.stackexchange.com/q/194877/2451 and links therein. $\endgroup$ – Qmechanic Sep 3 '16 at 21:07
  • $\begingroup$ You can construct Hamiltonian which produces correct EOMs. For example consider $H(p,x)=g^{\mu \nu}(x) p_{\mu} p_{\nu}$. There are no square roots, but you lose reparametrization invariance. $\endgroup$ – Blazej Sep 17 '16 at 23:37

I) Here we will assume that OP is taking about a relativistic point particle with zero spin in a Minkowski spacetime with metric $\eta_{\mu\nu}$ of sign convention $(−,+,+,+)$. Also we put $c=1$ for simplicity. (OP mentioned that the particle has charge but since it is free that is irrelevant.)

Note that the relativistic point particle has world-line reparametrization invariance, which is a gauge symmetry/redundancy in the formulation. We are (to a large extent) free to parametrize the world-line of the point particle in any way we wish. Let us call the world-line parameter for $\tau$ (which does not have to be the proper time). This gauge freedom can be encoded in an einbein field $e=e(\tau)>0$. The resulting Hamiltonian Lagrangian is$^1$

$$ L_H~:=~ p_{\mu} \dot{x}^{\mu} - \underbrace{\frac{e}{2}(p^2+m^2)}_{\text{Hamiltonian}}, \tag{1} $$

cf. e.g. this Phys.SE post. Here dot means differentiation wrt. $\tau$. The square of the $4$-momentum is $$ p^2~:=~ \eta^{\mu\nu} p_{\mu} p_{\nu}~=~-(p^0)^2+{\bf p}^2 ~=~-2p^+p^- + {\bf p}_{\perp}^2, \tag{2}$$ where we have used light-cone coordinates in the last expression.

II) Static gauge $x^0=\tau$. If we integrate out $p^0$ and $e$, we get OP's square root model $$ \left. L_H\right|_{x^0=\tau} \quad\stackrel{p^0}{\longrightarrow}\quad {\bf p}\cdot \dot{\bf x}- \underbrace{\left(\frac{1}{2e} + \frac{e}{2}({\bf p}^2+m^2)\right)}_{\text{Hamiltonian}} \quad\stackrel{e}{\longrightarrow}\quad {\bf p}\cdot \dot{\bf x} - \underbrace{\sqrt{{\bf p}^2+m^2}}_{\text{Hamiltonian}} \tag{3} .$$

For sufficiently short$^2$ times $\Delta \tau=\tau_f-\tau_i$, the path integral becomes$^3$

$$ \langle {\bf x}_f,\tau_f \mid {\bf x}_i,\tau_i\rangle\qquad\qquad\qquad$$ $$~=~ \underbrace{\frac{i\Delta\tau}{2\pi\hbar}}_{\text{Norm. factor}} \int_{\mathbb{R}^3} \!\frac{\mathrm{d}^3{\bf p}}{(2\pi\hbar)^3} \int_{\mathbb{R_+}} \!\frac{\mathrm{d}e}{2}~ \underbrace{\sqrt{\frac{2\pi\hbar}{ie\Delta\tau}}}_{\text{Gauss. } p^0\text{-int.}} \exp\left[\frac{i}{\hbar}\left( {\bf p}\cdot \Delta {\bf x} -\underbrace{\left( \frac{1}{2e} + \frac{e}{2}({\bf p}^2+m^2)\right)}_{\text{Hamiltonian}}\Delta\tau\right) \right]$$ $$~\stackrel{(5)}{=}~ \int_{\mathbb{R}^3} \!\frac{\mathrm{d}^3{\bf p}}{(2\pi\hbar)^3} \frac{1}{2\sqrt{{\bf p}^2+m^2}} \exp\left[\frac{i}{\hbar}\left( {\bf p}\cdot \Delta {\bf x} - \Delta \tau \underbrace{\sqrt{{\bf p}^2+m^2}}_{\text{Hamiltonian}}\right)\right], \tag{4} $$ which is the standard on-shell scalar propagator $\langle 0 | \phi (x_f)\phi (x_i)|0\rangle $ in QFT, cf. e.g. Refs. 1 & 2. As is well-known, it is Lorentz covariant and falls off exponentially outside the light-cone. In eq. (4) we have used the integral $$ \int_{\mathbb{R}_+} \!\frac{\mathrm{d}e}{\sqrt{e}}\exp\left[-ae-\frac{b}{e}\right] ~=~\sqrt{\frac{\pi}{a}} \exp\left[-2\sqrt{ab}\right], \qquad {\rm Re}(a), {\rm Re}(b)~>~0.\tag{5} $$

III) Light-cone gauge $x^+=\tau$. If we integrate out $p^-$ and $e$, we get
$$ \left. L_H\right|_{x^+=\tau} \quad\stackrel{p^-,~e}{\longrightarrow}\quad -p^+\cdot \dot{x}^- +{\bf p}_{\perp}\cdot \dot{\bf x}_{\perp} - \underbrace{\frac{{\bf p}_{\perp}^2+m^2}{2p^+}}_{\text{Hamiltonian}} .\tag{6} $$

IV) We stress that the Euler-Lagrange equations for either of the Hamiltonian Lagrangians (1), (3), and (6) lead to Hamilton's equations. The point is now that physical quantities should not depend on the choice of gauge-fixing. We are free to use the most convenient gauge choice. Each formulation (1), (3), and (6) are valid, and have their pros and cons. The static gauge choice (3) is disfavored because of the square root.


  1. M.E. Peskin & D.V. Schroeder, An Intro to QFT; eq. (2.50).

  2. M.D. Schwartz, QFT and the Standard Model; eq. (6.25).


$^1$ Strictly speaking, there are also Faddeev-Popov ghost terms and gauge-fixing terms, which we have ignored for simplicity. These action terms are consistently generated in the BFV formulation, cf. e.g. my Phys.SE post here. The normalization factor in eq. (4) can be derived via Gaussian integration in the BFV formulation over the 2 bosonic variables $x^0$, $B$; and the 4 fermionic variables $\bar{C}$, $P$, $C$, $\bar{P}$.

$^2$ Here we just consider a single time slice for simplicity. The full path integral is the continuum limit of multiple time slice discretizations with insertion of corresponding completeness relations. It turns out that the result (4) for the free theory does not depend on the number of time slice discretizations.

$^3$ The Gaussian integration over $p^0_E=i p^0_M$ becomes damped after a Wick-rotation $\tau_E=i\tau_M$, $x^0_E=ix^0_M$ to Euclidean signature.

| cite | improve this answer | |
  • $\begingroup$ Notes for later: Integral (5) is a modified Bessel function $K_{1/2}$. $\endgroup$ – Qmechanic Jul 9 '19 at 16:25

Well, not really. We COULD write hamiltonian as square root - if we know, what is a square root of an operator. Of course we have simple approximation:

$$\sqrt{1+x}=1+\frac x2-\frac{x^2}{8}+O(x^3)$$

Using this we could write your hamiltonian as:

$$\mathcal H=mc^2\sqrt{1+\frac{p^2}{m^2c^2}}=mc^2+\frac{p^2}{2m}+O(p^4).$$

The problem is that this form of hamiltonian allows us for superluminal particle transport - evolving particle with this hamiltonian gives nonzero probability in large distances - very problematic for hamiltonian grown from relativistic theory.

The answer is to find other solution, which after getting square gives us the hamiltonian before. That's how you derive the Dirac equation - you assume some matrices in equation and you hope that you can get a solution:

$$\mathcal H_D=\alpha^ip_ic+\beta mc^2$$ $$\mathcal H_D^2=p^2c^2+m^2c^4$$

The following procedures are basics of every relativistic quantum theory script, for example: http://inspirehep.net/record/459479/files/Forshaw.pdf , page 9.

| cite | improve this answer | |
  • $\begingroup$ Thanks for the informative reply. So I saw from the link you gave me, they applied the Lorentz Transformations on the Dirac equation. But not on the original Hamiltonian of the free particle in my post. Does this suggest that the LT can't really be applied properly on the original Hamiltonian? $\endgroup$ – Candy Man Apr 21 '16 at 3:14
  • $\begingroup$ Good question - I don't know. We could suppose that if energy-momentum equation is conserved, then LT are applied. It just does not conserve the locality principle, which is something else than Lorentz Transform. $\endgroup$ – Cheshire Cat Apr 21 '16 at 22:24
  • $\begingroup$ The question was about classical physics though. $\endgroup$ – Blazej Sep 17 '16 at 23:31

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.