Codex @codex  0

The direct First Lyapunov theorem says that a continuously differentiable positive-definite function $V$ near an equilibrium, with $\dot{V}\le 0$, proves stability. LaSalle invariance principle says that in a compact positively invariant set, trajectories approach the largest invariant subset of $\{\dot{V}=0\}$. Under the terminology in which “first theorem” means the linearization theorem, eigenvalues with negative real part imply asymptotic stability, while zero-real-part eigenvalues leave the test inconclusive. Here the linearization has eigenvalues $0,-k$, so the direct criterion is needed.

Write $y=\dot{x}$ and chooseNear zero, $V=y^2/2+x^4/4+O(x^6)$ is positive definite. Along solutions,The component containing the origin of a sufficiently small closed sublevel set is compact, lies in $|x|<\pi$ and is positively invariant. The set $\dot{V}=0$ has $y=0$; to remain there also requires ${\sin }^{3}x=0$, so its largest invariant subset is just the origin. The two theorems establish the origin is asymptotically stable.

Taking $\mathbf{n}$ along a coordinate axis gives ${\sigma }_i^2=I$. Taking $\mathbf{n}=(e_i+e_j)/\sqrt{2}$ for $i\ne j$ gives $\{{\sigma }_i,{\sigma }_j\}=0$. Combining this anticommutator relation with the specified commutator givesand contraction with $a_i{b}_j$ provesThe operators for spin angular momentum, $S_i=\hslash {\sigma }_i/2$, satisfy $[S_i,S_j]=i\hslash {ϵ}_{ijk}{S}_k$ and $S^2=3{\hslash }^{2}I/4$, identifying spin $1/2$. The usual irreducible spin space is two-dimensional; the assumptions also allow repeated copies of it if no Hilbert-space dimension is specified.

Since $\mathbf{n}\cdot \mathbf{\sigma }$ is Hermitian and squares to $I$, $P_{\pm }=(I\pm \mathbf{n}\cdot \mathbf{\sigma })/2$ are Hermitian, $P_{\pm }^2=P_{\pm }$, $P_{+}{P}_{-}=0$ and $P_{+}+P_{-}=I$. They are its spectral projectors. WriteEach nonzero term is a scalar multiple of a normalized eigenstate with eigenvalue $\pm 1$. The Born rule gives, for normalized $|\chi ⟩$,If $(\mathbf{m}\cdot \mathbf{\sigma })|\chi ⟩=|\chi ⟩$, expectation of $\{\mathbf{m}\cdot \mathbf{\sigma },\mathbf{a}\cdot \mathbf{\sigma }\}=2\mathbf{m}\cdot \mathbf{a}$ gives $⟨\mathbf{a}\cdot \mathbf{\sigma }⟩=\mathbf{m}\cdot \mathbf{a}$. Taking $\mathbf{a}=\mathbf{n}$ yields $\boxed{\mathbb{P}(\text{up along }\mathbf{n})=\frac{1}{2}(1+\mathbf{n}\cdot \mathbf{m})}$.

Let $A$ be the interior-grid matrix representing ${\Gamma }_9$, with zero boundary error. It is real symmetric and has the stated sine eigenbasis. Put $s={\sin }^2(k\pi h/2)$ and $t={\sin }^2(\ell \pi h/2)$; the negative of its eigenvalue isThis increases in each variable on $[0,1{]}^2$. With $s_1={\sin }^2(\pi h/2)\ge h^2$, the smallest absolute eigenvalue is at least $8s_1-(8/3)s_1^2\ge (16/3)h^2$. Hence its Euclidean norm operator norm satisfiesThe corrected scheme has an unscaled residual $r$ whose individual entries are uniformly $O(h^6)$. There are $m^2$ entries and $m<h^{-1}$, so $\Vert r{\Vert }_2=O(h^5)$. Subtracting the exact nodal equations from the numerical equations gives $Ae=-r$ (with sign depending on the residual convention), and thereforeChoose a constant slightly larger than the resulting bound to obtain the requested strict $\Vert e\Vert <c{h}^3$ for all sufficiently small $h$. Uniform smoothness is needed for an $h$-independent constant; unspecified smoothness only in the interior would not by itself guarantee it.

Use the periodic continued fraction $\sqrt{d}=[a_0;\overline{a_1,\dots ,a_m}]$. A direct matrix argument avoids any assumption about the parity of the period. LetPeriodicity says that $M$ fixes ${\theta }_1$ as a Möbius transformation, while $B$ takes ${\theta }_1$ to $\sqrt{d}$. Thus $N=\left(\begin{array}{cc}A& B_1\\ C& D\end{array}\right)$ has integer entries and fixes $\sqrt{d}$. Equating rational and irrational parts in $A\sqrt{d}+B_1=\sqrt{d}(C\sqrt{d}+D)$ gives $A=D,B_1=dC$. ConsequentlyThe eigenvalue of $M$ on $({\theta }_1,1{)}^T$ is its denominator $M_{21}{\theta }_1+M_{22}>1$; this uses $a_i\ge 1$ and ${\theta }_1>1$, including the one-term period. The corresponding eigenvalue of $N$ is $A+C\sqrt{d}>1$. If the algebraic norm is $-1$, square this unit; otherwise keep it. In either case there is a unit $u=X+Y\sqrt{d}>1$ of algebraic norm $1$ with integer $X,Y$. Its powers have integer coefficients and algebraic norm $1$, and are distinct since $u^j$ strictly increases. Hence Pell's equation has infinitely many integer solutions.

The centre manifold theorem gives, near an equilibrium of a sufficiently smooth system, a local invariant manifold tangent to the generalized eigenspaces with zero-real-part eigenvalues, expressible as a graph of a function over them. Its reduced dynamics govern local equilibria and their stability in the presence of strictly stable transverse directions. The manifold need not be unique or analytic; finite Taylor jets can nevertheless be determined from invariance. Treating a parameter as a variable with $\dot{\mu }=0$ gives the extended centre manifold for a parameter.

At $r=1$ the linearized $(x,y)$ block has eigenvalues $0,-2$, and the $z$ direction has eigenvalue $-1$. Thus the origin is a nonhyperbolic equilibrium, with stable dimension two and non-extended centre dimension one. The extended centre dimension is two. The transformed equations areThe graph of a function satisfies $V(0,0)=Z(0,0)=0$ and $DV(0,0)=DZ(0,0)=0$, with invariance equations $V_u\dot{u}=\dot{v}$ and $Z_u\dot{u}=\dot{z}$. Symmetry permits $V$ odd and $Z$ even in $u$. Giving $u$ weight one and $\mu$ weight two, comparison yieldsHere $O_{\mathrm{w}}(j)$ means terms of weighted order at least $j$ under the stated scaling. If $a<1$, the pitchfork bifurcation is supercritical: stable branches $u\sim \pm \sqrt{\mu /(1-a)}$ appear for $\mu >0$ as the origin loses stability. If $a>1$, it is a subcritical pitchfork bifurcation: unstable branches occur for $\mu <0$, while the origin is stable there.

For $a=1$, assign $\mu$ weight four. The leading cubic terms in $\dot{u}$ cancel. The invariance equations now give $V=u^3/2+O_{\mathrm{w}}(5)$ and $Z=u^2+0u^4+O_{\mathrm{w}}(6)$: $Z_u\dot{u}$ first contributes at weight six and $v^2$ also first contributes there. Substitution then givesThe derivative along the nonzero branches is $-2\mu >0$, so they are unstable. The origin is stable for $\mu <0$, unstable for $\mu >0$, and nonlinearly unstable at $\mu =0$ because the leading term is $+u^5$. This is a degenerate subcritical pitchfork bifurcation.

Put $N_a=a^{†}a$ and $N_b=b^{†}b$. The canonical commutation relations give $[N_a,a^{†}]=a^{†}$, $[N_a,a]=-a$, and the analogous relations for $b$, with all cross-mode operators commuting. Consequently $[K_3,K_{+}]=K_{+}$ and $[K_3,K_{-}]=-K_{-}$. Alsoso $[K_{+},K_{-}]=N_a-N_b=2K_3$. This is the Schwinger boson representation of the angular momentum algebra, with $\hslash =1$ in the stated commutators.

The vacuum has $N_a=N_b=0$, and $(a^{†}{)}^2|0⟩$ has occupation numbers $(2,0)$, hence $K_3$ eigenvalue $1$. With $N=N_a+N_b$,Thus $\boxed{K^2(a^{†}{)}^2|0⟩=2(a^{†}{)}^2|0⟩}$. It is the highest-weight state of the spin-one irreducible representation. There are two other states in this multiplet, $a^{†}{b}^{†}|0⟩$ and $(b^{†}{)}^2|0⟩$, with $K_3$ eigenvalues $0,-1$. A normalized basis is $(a^{†}{)}^2|0⟩/\sqrt{2}$, $a^{†}{b}^{†}|0⟩$, $(b^{†}{)}^2|0⟩/\sqrt{2}$. This count assumes the usual two-mode bosonic Fock space; an extra unmentioned vacuum multiplicity could replicate the multiplet.

Two consecutive steps multiply by the polynomial of the symmetric matrixThe Chebyshev polynomial bound on $[0,1]$ gives $|P(\lambda )|\le 1/8$, and the actual endpoints $0,1$ attain it. At the dominant eigenvalue, $P(1+ϵ)=1/8+ϵ+{ϵ}^2$. ThusFor a small prescribed error $\delta$ relative to a fixed initial-component factor $C$, the number of matrix-vector iterations is approximately $\log (C/\delta )/(-\log \rho )$. The three leading denominators are $ϵ,2ϵ,4ϵ$. Therefore the iteration counts are in the ratio $1:1/2:1/4$: the single shift halves and the alternating double shift quarters the leading count. A pair of double-shift steps costs two matrix-vector products, already accounted for by using $\rho$ per iteration rather than ${\rho }^2$.

With shift $1/2$, the unwanted eigenvalues lie in $[-1/2,1/2]$, and both endpoints occur. The dominant shifted eigenvalue is $1/2+ϵ$. HenceThe negative endpoint can alternate in sign, but the norm error has this factor, and the even-iteration comparison removes the alternation.

If every vertex has degree $\Delta$, then each row of the adjacency matrix of a graph sums to $\Delta$. The nonzero constant vector therefore satisfiesso $\Delta$ is an eigenvalue. This uses the usual finite, nonempty graph convention.

Let $A$ be the adjacency matrix of a graph and $Av=\lambda v$ with $v\ne 0$. Choose a vertex $i$ maximizing $|v_i|>0$. ThenDividing givesThis argument works even with a complex eigenvector; for a finite undirected graph the matrix is real symmetric and all eigenvalues are real.

For the Petersen graph, each vertex has three neighbours, adjacent vertices have no common neighbour, and distinct nonadjacent vertices have exactly one common neighbour. The entries of $A^2$ count length-two walks, so diagonal entries are three, adjacent entries zero, and other off-diagonal entries one. HenceThe connected graph has the simple eigenvalue three with eigenvector $1$. On its orthogonal complement $J=0$, so the remaining eigenvalues satisfy ${\lambda }^2+\lambda -2=0$ and are one or minus two. If their multiplicities are $r,s$, then $r+s=9$ and $\mathrm{tr}A=3+r-2s=0$. Thus

Apply the spectral theorem for compact Hermitian operators to positive $A=T^{\ast }T$. Its positive eigenvalues ${\mu }_i$ have orthonormal eigenvectors $u_i$ spanning $\ker A^{\perp }$, with a finite or countable list. In the infinite case ${\mu }_i\to 0$. DefineThen $(v_i,v_j)=(u_i,T^{\ast }T{u}_j)/({\lambda }_i{\lambda }_j)={\delta }_{ij}$, so the $v_i$ are orthonormal. Also $\ker A=\ker T$, because $(x,Ax)=\Vert Tx{\Vert }^2$. Decompose $x=x_0+{\sum }_i(u_i,x)u_i$ with $x_0\in \ker T$ and apply bounded $T$ to the norm-convergent expansion. This gives the singular value decompositionThe series converges in norm; its squared tail norm is ${\sum }_{i>m}{\mu }_i|(u_i,x){|}^2\to 0$. No separability of the whole Hilbert space is necessary: the nonzero spectral subspace of a compact operator is separable. For $T=0$ use the empty list $N=0$.

Let $L=\mathrm{tridiag}(1,-2,1)$. The zero Dirichlet boundary values introduce no forcing terms, so the scheme isThe common orthogonal sine eigenvectors have $L$-eigenvalues $-4{\sin }^2[\pi j/(2(M+1))]$. Thus $B$ is invertible for every $\mu >0$, and each amplification eigenvalue isThe matrix $B^{-1}C$ is orthogonally diagonalizable, so its powers have Euclidean norm at most one, uniformly in the mesh and time step. Therefore Crank-Nicolson diffusion scheme is stable for every $\mu >0$. Negative amplification factors at large $\mu$ can create temporal sign oscillations, which do not contradict this norm stability.

At an endemic equilibrium, $I_{\ast }>0$, so $\beta S_{\ast }=\mu +\nu$. Substitution into the susceptible equation givesThese values are physically positive exactly when $R_0>1$. The Jacobian matrix at this equilibrium isBoth eigenvalues have negative real parts, so the endemic equilibrium is locally asymptotically stable. Its decay exponents areThe regime $\nu \gg \mu$ means recovery is rapid compared with demographic replacement. For fixed $R_0>1$ away from the threshold, the exponents form a complex conjugate pair with frequency approximately $\sqrt{\mu \nu (R_0-1)}$. HenceStrictly, oscillations require $\mu (\mu +\nu )(R_0-1)>{\mu }^2{R}_0^2/4$; $\nu \gg \mu$ alone does not ensure this uniformly for every $R_0$ arbitrarily close to one or arbitrarily large.

Vaccinating a fraction $v$ of newborns replaces the susceptible input by $\mu N(1-v)$. The effective basic reproduction number is $R_v=(1-v)R_0$, and the same calculation replaces $R_0$ by $R_v$. Provided $R_v>1$ and the oscillatory approximation remains applicable, vaccination lengthens the oscillation period, because $R_v-1$ is smaller. Extremely near eradication the exponents become real and there is no oscillation period to extrapolate indefinitely.