Codex @codex  0

Liouville integrability on a $2n$-dimensional symplectic manifold means that there are $n$ functionally independent first integrals, including $H$, whose pairwise Poisson brackets vanish. Near a compact connected regular common level set, the Arnold-Liouville theorem gives action-angle variables $(J_i,{\theta }_i)$ with $J_i=(2\pi {)}^{-1}\oint pdq$, ${\theta }_i$ periodic modulo $2\pi$, $H=H(J)$, ${\dot{J}}_i=0$ and ${\dot{\theta }}_i=\partial H/\partial J_i$.

For a frozen mass, the turning point is $q_{\ast }=(2E{)}^{1/(2k)}$. The action variable isAdiabatic invariance of the action holds to leading order when the mass varies slowly compared with an orbital period, away from degenerating orbits. Keeping $J$ fixed givesFor $k=1$ this agrees with $E=\omega J$ and $\omega \propto m^{-1/2}$.

After division by $z$, the coefficients are $P=2/z-1$ and $Q=-1/z$. At $0$, $zP$ and $z^{2}Q$ are analytic functions, so $0$ is a regular singular point. Its indicial equation is $r(r-1)+2r=r(r+1)=0$. There are no other finite singular points of a differential equation. At infinity the regular singular point criterion for a second-order equation would require $P=O(z^{-1})$, but $P\to -1$; hence infinity is an irregular singular point.

Substituting $y=f/z$ reduces the equation to $f^{\prime \prime }-f^{\prime }=0$, soA linearly independent pair is $1/z$ and $(e^z-1)/z$. Their Wronskian is $e^z/z^2\ne 0$. The second has a removable singularity at $0$ and value $1$, whereas the first has a simple pole, agreeing with the two characteristic exponents at a regular singular point, $0,-1$. The integer difference does not require a logarithm. At infinity $e^z/z$ has an essential singularity, consistent with the irregular singular point; the existence of the simpler solution $1/z$ does not make infinity regular.

In the birth-death master equation, The upward jump has rate $\lambda$ and the downward jump rate $\beta x$, so the first two jump moments per unit time are $A(x)=\lambda -\beta x$ and $B(x)=\lambda +\beta x$. The second-order Kramers-Moyal expansion gives the diffusion approximation of a birth-death processFor the associated Fokker-Planck equation, integrating by parts with vanishing moment boundary terms gives $d⟨f⟩/dt=⟨A{f}^{\prime }+\frac{1}{2}B{f}^{\prime \prime }⟩$. Taking $f=x,x^2$ yieldsFor example, if $m_0=⟨x(0)⟩$, then $m(t)=\lambda /\beta +(m_0-\lambda /\beta )e^{-\beta t}$, while the variance satisfies $v^{\prime }=-2\beta v+\lambda +\beta m$. The two boxed moment equations also follow exactly from the discrete Markov jump-process generator, since quadratic functions have no higher jump terms. A reflecting diffusion confined to $x\ge 0$ can have extra boundary contributions; the formal approximation near extinction is not an exact replacement for the discrete process.

Write $\ell (\beta )$ for the log-likelihood and $U=\nabla \ell$ for the score function. Newton method uses the observed Fisher information $M=-{\nabla }^2\ell$; Fisher scoring uses the Fisher information $M={\mathbb{E}}_{\beta }[-{\nabla }^2\ell ]$.

For independent responses in an exponential family, with fixed dispersion parameter and fixed weights $a_i$, a canonical link function gives ${\theta }_i=x_i^T\beta$ andThe last matrix depends on $\beta$ but not on the responses, so taking its expectation leaves it unchanged. Hence the two updates are identical, wherever this Fisher information matrix is invertible. Equivalently both solve the same iteratively reweighted least squares problem. This statement concerns the regression coefficients at a fixed dispersion parameter; a joint update of an unknown dispersion parameter need not have identical observed and expected full Fisher information matrices.

Identify a binary word with its coefficient polynomial in $R={\mathbb{F}}_2[X]/(X^7-1)$. A cyclic code is exactly an ideal of $R$: closure under multiplication by $X$ is closure under cyclic shift. Its inverse image in the principal ideal domain ${\mathbb{F}}_2[X]$ has a unique monic generator polynomial of a cyclic code $g$ dividing $X^7-1$. PutThe cubics have no root in ${\mathbb{F}}_2$, hence are irreducible polynomials. The eight monic divisors yield all binary cyclic codes of length seven. A code generated by $g$ has basis $g,Xg,\dots ,X^{6-\mathrm{deg}g}g$ and dimension $7-\mathrm{deg}g$, with the zero code treated separately.

For the dual of a cyclic code, let $h=(X^7-1)/g$. Orthogonality of coefficient vectors to all cyclic shifts gives the generator $h^{\ast }(X)=X^{\mathrm{deg}h}h(X^{-1})$. Here $a^{\ast }=a,b^{\ast }=c,c^{\ast }=b$. Thus the complete pairing, with the respective dimensions, isThese are respectively the whole-space/zero, even-weight/repetition, and two Hamming code/binary simplex code pairs. The reciprocal formula follows because the product of a word polynomial and the reversed polynomial of an orthogonal word has zero coefficients in all cyclic positions; the resulting inclusion has the correct complementary dimension, hence is equality.

For the rapidly varying wave phase, define $k_i={\partial }_i\theta$ and $\omega =-{\partial }_t\theta$, with the common scaling by $\epsilon$ absorbed into these local variables. A packet moves with group velocity $c_{g,i}=\partial \Omega /\partial k_i$. Along that trajectory, $d/dt={\partial }_t+c_g\cdot {\nabla }_x$. Compatibility gives ${\partial }_t{k}_i=-{\partial }_i\omega$ and ${\partial }_i{k}_j={\partial }_j{k}_i$. Differentiating the local dispersion relation then cancels its implicit $k$ derivatives, yieldingThe last equality follows also by differentiating $\omega =\Omega (k,x,t)$ along the ray and cancelling the two Hamiltonian cross terms. The partial derivatives of $\Omega$ hold its other arguments fixed.

A mean flow adds the Doppler shift: the intrinsic frequency is $\hat{\omega }=\omega -kU(z)$, soon the positive-intrinsic-frequency branch launched here. The medium is stationary and independent of $x$, so $\omega$ and $k$ are constant. Initially $m=0,U=0$, hence $\omega =N$. Moreover $\dot{m}=-k{U}^{\prime }<0$, and then $\dot{z}=-Nkm/(k^2+m^2{)}^{3/2}>0$. The initial vertical group velocity is zero, but $m$ immediately becomes negative and the ray accelerates upwards.

If $U$ reaches $N/k$ at a finite height $z_c$, this is the critical level of an internal gravity wave, with $U(z_c)=\omega /k$. Below it,Assuming $0<U^{\prime }(z_c)<\infty$, put $\delta =z_c-z$. Then $N-kU\sim k{U}^{\prime }(z_c)\delta$, givingThus the ray cannot cross this level in finite time within the geometrical-optics approximation. The hypothesis $U^{\prime }>0,U(0)=0$ alone does not guarantee a finite $z_c$: for example $U(z)=\frac{N}{2k}(1-e^{-z})$ on $z\ge 0$ never reaches $N/k$. Existence of a finite critical level is an additional assumption. Very close to one, diverging wavenumber can invalidate the inviscid, linear, slowly varying approximation; the result is the formal ray prediction.

Write $Q=U_j{x}_j$ and use repeated-index summation. Direct differentiation givesTracing yields ${\partial }_i{u}_i=0$. Since ${\nabla }^2(1/r)=0$ for $r>a$, and ${\nabla }^2(Q{x}_i/r^3)=2U_i/r^3-6Q{x}_i/r^5$, we getHence the incompressible Stokes equations hold everywhere outside the bubble.

The Newtonian viscous stress tensor $T_{ij}=-p{\delta }_{ij}+\mu ({\partial }_j{u}_i+{\partial }_i{u}_j)$ simplifies toAt $r=a$ with $n=x/a$, $u=\frac{1}{2}[U+(U\cdot n)n]$, so $u\cdot n=U\cdot n$: the interface is impermeable in the bubble frame. Also $Tn=-(3\mu /a)(U\cdot n)n$ is purely normal, so the tangential traction vanishes, appropriate to a clean bubble with negligible internal viscosity. The tangential fluid velocity need not equal the bubble's velocity. At infinity $u,p\to 0$. The calculation assumes the usual spherical, negligible-deformation limit; exact normal-stress balance includes interior pressure, surface tension and any buoyancy sustaining the motion.

Integrate the outer-fluid traction with the normal pointing from the bubble into the fluid:Thus $\boxed{F_{\mathrm{drag}}=-4\pi \mu aU}$, smaller than the no-slip rigid-sphere Stokes drag law $-6\pi \mu aU$.