special relativity - Why don't people use Hamilton's equations for a relativistic free charged particle?

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)

Answer

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.

References:

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.