Phase reduction ($\dot{\theta }=\omega +ϵZ(\theta )\cdot p$)

Near a hyperbolic attracting limit cycle, a weak perturbation can be projected onto the neutral phase direction while transverse perturbations decay. Solving the periodic adjoint variational equation gives a phase response function. Averaging its inner product with weak forcing or coupling yields the first-order phase dynamics.

Linearizing weakly coupled locked phases gives $0=\xi -\Omega 1-kL\theta$, where $L$ is the directed graph Laplacian. Multiplying by a normalized left null vector ${\pi }^{\mathsf{T}}$ eliminates coupling and gives the displayed weighted frequency. In a balanced network, including symmetric reciprocal coupling, $\pi =1/N$ and the frequency is the arithmetic mean. A zero diagonal in the coupling matrix alone does not imply balance. For two directed weights $p,q>0$, the weights are $(q,p)/(p+q)$.

A periodic phase response function describes the first-order phase shift produced by a perturbation along a limit cycle. With period-one orbit $U$, choose the periodic adjoint solution normalized by ${\int }_0^{1}Y\cdot U^{\prime }ds=1$. The Fredholm solvability condition then projects forcing $p$ to ${\int }_0^{1}Y\cdot pds$. For voltage-difference coupling, the averaged interaction is $h(\varphi )={\int }_0^1{Y}^1(s)U^1(s+\varphi )ds$.