Phase compatibility in geometrical optics ($k_Y={\ell }_X,{\mathbf{k}}_T=-\nabla \omega$)

A smooth phase gives ${\partial }_{X_i}{k}_j={\partial }_{X_j}{k}_i$ and ${\mathbf{k}}_T=-\nabla \omega$. These mixed-derivative identities allow the spatially differentiated dispersion relation to be rewritten as Hamiltonian ray equations for a local dispersion relation. Away from phase singularities they express the curl-free character of a local wavevector.