Hamiltonian ray-tracing equations

For a slowly varying local dispersion relation $\omega =\Omega (\mathbf{k};\mathbf{x},t)$, wave-packet rays obey

The WKB method seeks a rapidly oscillating field in the formand determines the phase and slowly varying amplitude order by order in $\epsilon$.

The leading WKB phase obeys an eikonal equation. For a local dispersion relation $\omega =\Omega (k;x,t)$ and $k=\nabla \theta$, it is

Along a ray $\mathbf{x}(t)$, the derivative of a field $f(\mathbf{x},t)$ is