Codex @codex  0

First establish rational conjugation of finite-index modular subgroups without assuming that $\Gamma$ is a congruence subgroup. Multiply $\gamma$ by a positive integer to obtain an integral matrix $A$, and let $D=\det A>0$. Conjugation is unchanged by this scalar. If $h=I+DB\in \Gamma (D)$, thenThus the principal congruence subgroup $\Gamma (D)$ is contained in $H=S{L}_2(\mathbb{Z})\cap {\gamma }^{-1}S{L}_2(\mathbb{Z})\gamma$. It has finite index because reduction modulo $D$ has finite image. Inside $H$, pullback under conjugation of $\Gamma$ has relative index at most $[S{L}_2(\mathbb{Z}):\Gamma ]$. ConsequentlyThis argument does not assert that an arbitrary finite-index subgroup contains a principal congruence subgroup.

Use the determinant-normalized slash operatorThe positive real power of the determinant is used; on $S{L}_2(\mathbb{Z})$ this reduces to the usual slash operator for modular forms. The automorphy factor identity gives the right-action rule $(f{|}_k\gamma ){|}_k\eta =f{|}_k(\gamma \eta )$.

A modular form on a finite-index subgroup of integer weight $k$ is a holomorphic function on the complex upper half-plane, invariant under this weight-$k$ action of $\Gamma$, and holomorphic at a cusp at each of its cusps. A cusp of a modular group is a $\Gamma$ orbit in ${\mathbb{P}}^1(\mathbb{Q})$. If $\sigma \in S{L}_2(\mathbb{Z})$ carries infinity to its representative, choose a positive integer $h$ with $\sigma T^h{\sigma }^{-1}\in \Gamma$. Such $h$ exists by finite index. Then $f{|}_k\sigma$ is periodic and has a convergent expansion in $q_h=e^{2\pi iz/h}$ near zero; holomorphy means no negative exponents, and being a cusp form means zero constant term. Using an actual translation period avoids possible signs if a smaller width of a cusp is defined only modulo the center, particularly in odd weights.

For cusp holomorphy under rational slash operators, choose $\sigma \in S{L}_2(\mathbb{Z})$ with $\sigma \infty =\gamma \infty$, possible by completing a primitive integer pair to a determinant-one matrix. Then ${\sigma }^{-1}\gamma =\left(\begin{array}{cc}a& b\\ 0& d\end{array}\right)$, with $a/d>0$. Up to a nonzero constant factor,The imaginary part of the argument tends to infinity with that of $z$, so this remains bounded by the cusp expansion of $f{|}_k\sigma$. It tends to zero if $f$ is a cusp form. Moreover $f{|}_k\gamma$ is invariant under ${\Gamma }^{\prime }$: for $\eta \in {\Gamma }^{\prime }$, $\gamma \eta {\gamma }^{-1}\in \Gamma$ and the right-action rule applies. Finite index gives a translation period for ${\Gamma }^{\prime }$, so boundedness is a removable singularity at zero in that periodic parameter. This proves holomorphy at infinity. For every other cusp, apply the same argument to the rational matrix $\gamma \rho$, with $\rho \in S{L}_2(\mathbb{Z})$. Thus all cusp conditions hold, and

For the character twist by rational translations of a cusp form, put $B_j=\left(\begin{array}{cc}1& j/N\\ 0& 1\end{array}\right)$ and $K={\Gamma }_1(N)\cap {\Gamma }_0(N^2)$. For $\eta =\left(\begin{array}{cc}a& b\\ c& d\end{array}\right)\in K$, direct conjugation givesIndeed $N^2\mid c$ and $a\equiv d\equiv 1(\mathrm{mod}N)$. Every $f{|}_k{B}_j=f(z+j/N)$ is therefore $K$-invariant and vanishes at all its cusps by the preceding rational-translate argument. Their finite weighted sum is a cusp form, for every Dirichlet character:

There is, however, a missing primitivity hypothesis in the printed final expansion claim. The exact Fourier expansion of a modular form is alwaysValues of a Dirichlet character on units have modulus one, so $\chi (j{)}^{-1}=\overline{\chi (j)}$. For unit $n$, substitution gives $S_{\overline{\chi }}(n)=\chi (n)g(\overline{\chi })$, where $g$ is the Gauss sum of a Dirichlet character. For nonunit $n$, this vanishing formula requires a primitive Dirichlet character.

Here is its proof in that case. Choose a prime $p\mid (n,N)$. Primitivity supplies a unit $u\equiv 1(\mathrm{mod}N/p)$ with $\chi (u)\ne 1$: otherwise the character would factor through the surjective reduction to units modulo $N/p$. Surjectivity follows by lifting a unit and, if needed, adjusting the lift to avoid the additional prime $p$, using the Chinese remainder theorem. Multiplication by $u$ fixes $e^{2\pi inj/N}$ because $p\mid n$, but multiplies the character factor by a nontrivial constant. Hence $S_{\overline{\chi }}(n)=0$. The finite Fourier transform of a primitive Dirichlet character now gives the corrected formulaThe constant is nonzero: finite exponential orthogonality gives ${\sum }_{n\mathrm{mod}N}|S_{\overline{\chi }}(n){|}^2=N\phi (N)$, whereas the proved formula makes this $\phi (N)|g(\overline{\chi }){|}^2$. Thus $|g(\overline{\chi })|=\sqrt{N}$.

For a concrete counterexample to the printed unrestricted claim, take $N=2$, the principal Dirichlet character, and $f=\Delta$. The translation sum is $\Delta (z+1/2)={\sum }_n(-1{)}^n\tau (n)q^n$, whose $q^2$ coefficient is $-24$. Any constant multiple of the proposed odd-index-only series has $q^2$ coefficient zero. Thus the general modularity conclusion is proved, while the claimed simplification is false without the stated extra hypothesis.

Write $\Gamma =S{L}_2(\mathbb{Z})$ and $G=G{L}_2(\mathbb{Q}{)}^{+}$. The double-coset Hecke algebra consists of $\Gamma$-bi-invariant complex functions on $G$ supported on finitely many double cosets, with convolutionEach double coset has finitely many orbits under left multiplication by $\Gamma$, by rational conjugation of finite-index modular subgroups. Thus the sum is finite and independent of representatives. The characteristic functions of the double cosets form its basis, and its right action on invariant modular forms is $f\ast h={\sum }_{\Gamma x}h(x)f{|}_{k}x$, using the determinant-normalized slash operator.

For a positive integer $n$, let ${\mathcal{M}}_n$ be all integral two-by-two matrices of positive determinant $n$. It is $\Gamma$-bi-invariant. Define $T(n)$ to be its indicator function, equivalently the sum of the distinct double cosets it contains, each with coefficient one. This is the convention consistent with the requested formula. For composite $n$, it need not be the single double coset of $\mathrm{diag}(1,n)$: for instance $2I\in {\mathcal{M}}_4$ cannot lie in that double coset, since multiplying by unimodular integral matrices preserves the greatest common divisor of the entries.

Prove all the needed subgroup facts directly. For a finite-index subgroup $L\subseteq {\mathbb{Z}}^2$, take $a$ to be the least positive first coordinate appearing in $L$ and $d$ the least positive second coordinate on its intersection with the second axis. Euclidean division then shows that the first-coordinate projection is $a\mathbb{Z}$ and $L\cap (\{0\}\times \mathbb{Z})=\{0\}\times d\mathbb{Z}$. Positivity follows, for example, because the finite quotient group kills a nonzero multiple of each coordinate vector. Choose $(a,b)\in L$ and reduce $b$ modulo $d$ to $0\le b<d$. Every vector of $L$ has first coordinate a multiple of $a$, and subtracting that multiple of $(a,b)$ leaves a multiple of $(0,d)$. Therefore these two vectors form a basis of $L$. The parameters $a,d,b$ are unique. Reducing the first coordinate modulo $a$ and then the second modulo $d$ gives precisely $ad$ quotient representatives, so $[{\mathbb{Z}}^2:L]=ad$.

Apply this row Hermite normal form in rank two to the row lattice of an integral matrix $A$ with $\det A=n$. Its row lattice contains $n{\mathbb{Z}}^2$, since $\mathrm{adj}(A)A=nI$, so it has finite index. Its two rows and the displayed two rows are bases of the same row lattice. The two inverse change-of-basis matrices have integer entries, so their determinants are integers whose product is one. Thus the change-of-basis matrix has determinant $\pm 1$; since both orientations are positive, its determinant is one. Thus each orbit under left multiplication by $\Gamma$ has exactly one of the determinant-n matrix representatives for Hecke operatorsThere are ${\sum }_{d\mid n}d$ such representatives, proving finiteness as well as the formula. The subgroup argument is a proof of the relevant Hermite normal form, not an invocation of an unproved lattice classification.

Consequently the normalized Hecke operator isIt preserves $M_k(\Gamma )$: right multiplication by $\Gamma$ permutes the left-multiplication orbits in ${\mathcal{M}}_n$, giving invariance, and cusp holomorphy under rational slash operators gives the holomorphy of each term at every cusp.

For a triangular representative the determinant-normalized slash operator is $f{|}_{k}A=n^{k/2}{d}^{-k}f((az+b)/d)$. Summing the Fourier expansion of a modular form over $b$ kills every index not divisible by $d$, by finite exponential orthogonality. ThusThe Fourier coefficients of a composite-index Hecke operator are thereforeIn particular $a_1(T_{n}f)=a_n(f)$. If $T_{n}f=\lambda f$, comparison of the $q$ coefficients gives

Finally suppose $a_0(f)\ne 0$ and $f$ is a simultaneous eigenfunction. Constant coefficients force ${\lambda }_n={\sigma }_{k-1}(n)$, so $a_n(f)=a_1(f){\sigma }_{k-1}(n)$ for all $n\ge 1$. Weight two cannot occur, by vanishing of weight-two level-one modular forms. For even $k\ge 4$, use the normalized Eisenstein series with its Fourier expansion of a normalized Eisenstein serieswhere $B_k$ is the Bernoulli number. Then $h=f-a_0(f)E_k$ is a cusp form with coefficients $(a_1(f)-a_0(f)c_k){\sigma }_{k-1}(n)$. The Fourier coefficient bound for a cusp form bounds these by $O(n^{k/2})$, but at arbitrarily large primes $p$ their magnitude is $|a_1-a_0{c}_k|(1+p^{k-1})$. Since $k-1>k/2$, the coefficient factor must vanish. All coefficients of $h$ then vanish, including its constant coefficient, so $h=0$ by its cusp expansion and the identity theorem. This proves the noncuspidal level-one Hecke eigenform characterization:

The two Lie groups have the same local infinitesimal structure, but different global topology. This difference determines which Lie algebra representations integrate to representations of each group.

An element of the SU(2) group is a unitary matrix of determinant one. Orthogonality of its columns and its determinant give the unique formThus its group manifold is the unit three-sphere in ${\mathbb{C}}^2\simeq {\mathbb{R}}^4$: this is SU(2) as the three-sphere. In particular it is compact, connected and simply connected. In terms of the Pauli matrices one may also write $U=a_{0}I-i\mathbf{a}\cdot \mathbf{\sigma }$ with real coefficients satisfying $a_0^2+{\mathbf{a}}^2=1$.

Differentiate $U(t{)}^{†}U(t)=I$ and $\det U(t)=1$ at $U(0)=I$. The tangent space consists of traceless skew-Hermitian matrices:The Lie bracket of a Matrix Lie group is the matrix commutator. With $T_a=-i{\sigma }_a/2$, the Pauli matrix multiplication law givesThis derives the SU(2) Lie algebra as a three-dimensional real Lie algebra. The Hermitian physics generators $J_a={\sigma }_a/2$ instead obey $[J_a,J_b]=i{ϵ}_{abc}{J}_c$; they are $i$ times the skew-Hermitian tangent generators, so these are consistent conventions.

The SO(3) group consists of real orthogonal matrices with determinant one. Differentiating $R(t{)}^{T}R(t)=I$ at the identity givesThe determinant condition gives no additional infinitesimal constraint because a skew-symmetric matrix already has trace zero. Define $L_a\mathbf{v}={\mathbf{e}}_a\times \mathbf{v}$. The cross product identity impliesTherefore the SO(3) Lie algebra has the same structure constants and $T_a\mapsto L_a$ is a Lie algebra isomorphism.

The global relation is the Adjoint double cover from SU(2) to SO(3). For $V=\mathbf{v}\cdot \mathbf{\sigma }$, define $R(U)$ byConjugation preserves the real space of traceless Hermitian matrices and its inner product $\frac{1}{2}\mathrm{tr}(VW)=\mathbf{v}\cdot \mathbf{w}$. Hence $R(U)$ is orthogonal. Continuity and connectedness, together with $R(I)=I$, put it in $SO(3)$. Composition of conjugations makes $R$ a group homomorphism. If $R(U)=I$, then $U$ commutes with every Pauli matrix, hence is scalar; unitarity and determinant one leave precisely $U=\pm I$. The differential sends $T_a$ to $L_a$, so it is an isomorphism. More concretely,induces rotation through angle $\theta$ about $\mathbf{n}$, by the Rodrigues rotation formula. Every three-dimensional rotation has such an axis and angle, proving surjectivity. ConsequentlyThe matrices $U$ and $-U$ are antipodal points on the three-sphere, so the SO(3) group manifold is Real projective space ${\mathbb{RP}}^3$. Equivalently, the closed axis-angle ball $\Vert \theta \mathbf{n}\Vert \le \pi$ has opposite boundary points identified. The fundamental group is ${\pi }_1(SO(3))\simeq {\mathbb{Z}}_2$, whereas ${\pi }_1(SU(2))=0$. A $2\pi$ rotation lifts from $I$ to $-I$; a $4\pi$ rotation returns to $I$. Thus the covering is the universal cover and the groups are not globally isomorphic.

For representation theory, specify finite-dimensional complex continuous representations. Compactness permits an invariant Hermitian inner product, obtained by averaging against Haar measure, and therefore complete reducibility. Complexifying either real Lie algebra gives the sl2 Lie algebra. Its finite-dimensional irreducibles are indexed by $n\in {\mathbb{Z}}_{\ge 0}$, have highest weight $n$, and have dimension $n+1$. By integration of a Lie-algebra representation, since $SU(2)$ is simply connected, every such Lie algebra representation integrates uniquely. The resulting homogeneous polynomial representation of SU2 isIn the spin angular momentum notation $j=n/2$, its Hermitian $J_3$ eigenvalues are $j,j-1,\dots ,-j$. The central matrix $-I$ acts on the symmetric power by $(-1{)}^n$, so descent of an SU(2) representation to SO(3) occurs exactly when $n$ is even. HenceFor a reducible representation, every summand must satisfy the descent condition. The spin-one-half doublet is a genuine representation of the covering group but does not define a single-valued representation of $SO(3)$; the spin-one triplet does and is its vector representation. The distinction is topological, rather than a difference in their isomorphic Lie algebras.

Put $q=r_0^2$ and $\eta =d\psi +\frac{1}{2}\cos \theta d\varphi$. Expanding the metric tensor in the stationary direction gives the useful exact identitiesTaking the exterior derivative givesThe fibered form of the metric tensor has determinant $\det g=-r^6{\sin }^2\theta /16$. Its inverse components needed for the flux areThereforesince $h(1-\alpha \Omega )=1$. This calculation includes the rotational mixed term; dropping it prematurely would miss the exact cancellation.

With the stipulated negative coordinate orientation, ${ϵ}_{tr\psi \theta \varphi }=-r^3\sin \theta /4$. The pullback of the Hodge star operator to the constant-$t,r$ three-dimensional spacelike submanifold is consequentlyThe coordinate ranges and standard sphere identifications giveThus the five-dimensional Komar mass, in the normalization given and $G_5=1$, isThe flux is already independent of $r$ in the vacuum region, so its limit at infinity has the same value. The negative orientation in the question is crucial for the positive sign; reversing that orientation reverses the flux.

Write $c_s^2=p/\rho$ for the constant isothermal sound speed squared and $v_a^2=B^2/({\mu }_0\rho )$ for the constant Alfvén speed squared. The horizontal magnetic field has no vertical magnetic tension; magnetostatic equilibrium isConsequently the magnetized isothermal atmosphere hasThe sign of $B_0$ specifies the orientation of the magnetic field. The adiabatic sound speed is $v_s^2=\gamma p/\rho =\gamma c_s^2$, also constant.

For the proposed normal mode, let $A=\exp (i\mathbf{k}\cdot \mathbf{x}+z/(2H)-i\omega t)$, $a=1/(2H)$, and temporarily write $\mathbf{\xi }$ for its constant complex amplitude. The divergence amplitude is $\Delta =i{k}_x{\xi }_x+i{k}_y{\xi }_y+(i{k}_z+a){\xi }_z$. Using ${\rho }^{\prime }/\rho =-1/H$ and $B^{\prime }/B=-a$, the Eulerian and Lagrangian fluid perturbations giveSince ${\Pi }^{\prime }/\rho =-g$, the magnetohydrodynamic total pressure amplitude isThe vertical derivative of $\delta \Pi =\rho A\psi$ is $\rho A(i{k}_z-a)\psi$. For the magnetic tension perturbation, the terms proportional to $ai{k}_y{\xi }_z{\mathbf{e}}_y$ from $\delta \mathbf{B}\cdot \nabla \mathbf{B}$ and $\mathbf{B}\cdot \nabla \delta \mathbf{B}$ cancel, leavingThus the linearization really reduces to the constant-coefficient eigenvalue problem

For a static ideal magnetohydrodynamics equilibrium with conservative gravitational acceleration, the magnetohydrodynamic energy principle gives a Hermitian operator for the restoring force in the mass density weighted inner product, provided the surface terms vanish. The factor $e^{z/(2H)}$ removes the exponential mass density weight. This can be checked directly, without assuming stability: set $s=v_s^2$, $b=v_a^2$, and $q=s+b$. Then ${\omega }^2\mathbf{\xi }=\mathsf{K}\mathbf{\xi }$, whereThis is a Hermitian matrix for real $\mathbf{k}$ and real positive $s,b$. Therefore every squared frequency ${\omega }^2$ is real; negative values give exponential growth, rather than an oscillation with a complex squared frequency. In an infinite atmosphere these are generalized plane waves, not finite-total-energy global modes.

At ${\omega }^2=0$ and $k_y\ne 0$, the $y$ equation, after cancelling its magnetic terms, gives $s\Delta =g{\xi }_z$. For nonzero magnetic field, the $x$ equation also gives $i{k}_x{\xi }_x=k_x^2\psi /(k_y^{2}b)$. Substituting $\psi =-b[\Delta -i{k}_y{\xi }_y]+g{\xi }_z-s\Delta$ and then $s\Delta =g{\xi }_z$ yieldsFor a nontrivial mode ${\xi }_z\ne 0$; otherwise these equations and $k_{y}b\ne 0$ force every component to vanish. Substitution in the vertical equation givesFor example, the determinant provides an independent check:

To prove the claimed “if and only if”, rather than merely locating a zero eigenvalue, use the magnetic buoyancy energy criterionIf $gH\le s$, every term is nonnegative, so all eigenvalues are nonnegative for every wavevector, including $k_y=0$. If $gH>s$, choose $k_z=0$, any fixed $k_x\ne 0$, and ${\xi }_z=1$, ${\xi }_x=-a/(i{k}_x)$, ${\xi }_y=g/(si{k}_y)$. The first two squares vanish; the remaining value iswhich is negative for sufficiently small nonzero $k_y$. The Rayleigh quotient of this displacement is therefore negative, proving a negative eigenvalue and an unstable normal mode. This uses precisely the assumption that boundary conditions do not exclude the chosen wavelengths. Equality is a neutral threshold, not exponential growth.

Finally the plasma beta is $\beta =p/[B^2/(2{\mu }_0)]=2c_s^2/v_a^2$. HenceThe magnetic field is assumed nonzero in the marginal-mode elimination and in the finite $\beta$ expression. The positive-sound-speed, nonmagnetic limit remains covered by the energy identity; for the usual $\gamma >1$ it is stable. Parker instability occurs because magnetic pressure supports part of the atmosphere's weight while the gas can drain along gently bent magnetic field lines. The resulting buoyant displacement can release Newtonian gravitational potential energy; choosing long wavelengths along the magnetic field makes the stabilizing magnetic tension small. The completed-square calculation quantifies this competition, including the adiabatic process restoring force.

Reset $t=0$ when gravity is switched off. The subsequent overdamped relaxation of a clamped--free filament obeys ${\zeta }_{\perp }{h}_t=-B{h}_{xxxx}$ with the same homogeneous endpoint conditions. Separate a normal mode as $h=\varphi (x/L)e^{-\Gamma t}$. Writing $s=x/L$ givesThe two clamp conditions reduce the general hyperbolic/trigonometric solution toThe two free-end conditions have matrixIts determinant is $2(1+\cosh \alpha \cos \alpha )$. Hence the clamped--free bending mode values satisfyThere is no zero mode: at zero eigenvalue $\varphi$ is cubic, and the four clamp/free conditions force all coefficients to vanish. With the following normalization a mode isThe fourth-order operator is self-adjoint and positive under these endpoint conditions, since $\int f{g}^{\prime \prime \prime \prime }=\int f^{\prime \prime }{g}^{\prime \prime }$. Its orthogonal eigenfunction expansion givesFor the previously attained steady profile, the first coefficient is nonzero: $h_s$ is positive and the first mode can be chosen positive in the interior. Two applications of integration by parts, using both functions' endpoint conditions, also giveThus the dominant long-time answer, including its shape and amplitude, isThe next mode decays much faster. For other initial shapes with an accidentally zero first-mode projection, the lowest mode with $c_n\ne 0$ determines the asymptotic behaviour. The decay time has the correct units and exhibits the strong $L^4$ dependence characteristic of bending opposed by local viscous drag.

We prove the conjugate product of a quadratic ideal identity $\mathfrak{a}\overline{\mathfrak{a}}=(N(\mathfrak{a}))$, which applies to the given two-generator ideal without assuming its displayed generators form an integral basis.

The quadratic ring of integers of a number field has a basis $1,\omega$ over $\mathbb{Z}$: the vector 1 is primitive, since a rational algebraic integer must be an integer, and extends to an integral basis. Write ${\omega }^2-T\omega +U=0$, with $T,U\in \mathbb{Z}$. Let $B$ be the least positive integer in $\mathfrak{a}$ and let $s>0$ generate the set of coefficients of $\omega$ in elements of $\mathfrak{a}$. Its lattice basis can be chosen as $B,r+s\omega$, with $r\in \mathbb{Z}$. Since $B\omega \in \mathfrak{a}$, comparison of coefficients shows $s\mid B$ and $s\mid r$. Write $B=sd$, $r=se$ and $\theta =e+\omega$. ThenThe scaled ideal $\mathfrak{j}=s^{-1}\mathfrak{a}$ remains an ${\mathcal{O}}_L$-ideal. From ${\theta }^2-T^{\prime }\theta +U^{\prime }=0$ we get $U^{\prime }\in \mathfrak{j}\cap \mathbb{Z}=d\mathbb{Z}$, so $d\mid U^{\prime }$. The product $\mathfrak{j}\bar{\mathfrak{j}}$ is generated as an ideal by $d^2,d\theta ,d\bar{\theta },U^{\prime }$ and is contained in $(d)$.

Moreover $\gcd (d,T^{\prime },U^{\prime }/d)=1$. If a prime $\ell$ divided all three, then $\ell \mid T^{\prime }$ and ${\ell }^2\mid U^{\prime }$, so $\theta /\ell$ would satisfy the monic integer polynomial $z^2-(T^{\prime }/\ell )z+U^{\prime }/{\ell }^2$. It would lie in ${\mathcal{O}}_L$, impossible because its $\omega$ coefficient in the integral basis is $1/\ell$. Bezout identity now expresses one as an integer combination of $d,T^{\prime },U^{\prime }/d$. Multiplying by $d$, the terms are $d^2$, $d{T}^{\prime }=d\theta +d\bar{\theta }$, and $U^{\prime }$, all in the product. Hence $d$ belongs to it, proving $\mathfrak{j}\bar{\mathfrak{j}}=(d)$.

The determinant of the lattice basis $sd,s(e+\omega )$ gives the ideal norm $N(\mathfrak{a})=s^{2}d$. ConsequentlyIn particular $⟨b,\alpha ⟩⟨b,\bar{\alpha }⟩$ is principal. This does not require the original $b$ to be nonzero or the ideal itself to be principal.

A fraction $1-k$ of adults survives each time step; $k$ is the adult mortality fraction. Adults produce $\mu a_n$ larvae, which mature one step later with recruitment reduced by adult-dependent competition through $1/(1+a_n)$. Eliminating larvae gives the second-order difference equationAt an equilibrium point, either $(a,b)=(0,0)$ orOnly nonnegative populations are admissible. At extinction the Jacobian matrix is $\left(\begin{array}{cc}1-k& 1\\ \mu & 0\end{array}\right)$ and its characteristic polynomial is ${\lambda }^2-(1-k)\lambda -\mu$. Its positive root exceeds $1$ exactly when $\mu >k$. For $0<\mu <k$ both roots have modulus less than $1$; at $\mu =k$ the roots are $1,-k$, so linearization is marginal. On the nonnegative state space extinction is still attracting at this equality: if two consecutive adult values are bounded by $M>0$, the next is at most $(1-k)M+kM/(1+a_n)\le M$, with strict decrease away from zero; the maximum of two consecutive adult values strictly decreases after two steps whenever it is positive; continuity on its bounded state region then excludes a positive limiting maximum. Thus the extinction equilibrium point is unstable exactly when the positive equilibrium point exists, with equality understood through nonlinear rather than strict linear stability.

At the positive equilibrium point, the trace and determinant of the Jacobian matrix are $T=1-2k+k^2/\mu$ and $-k$. The Jury stability criterion for ${\lambda }^2-T\lambda -k$ requires $1-T-k>0$, $1+T-k>0$ and $1-k>0$. The strict linear stability analysis conditions areFor $k\le 2/3$ there is no upper bound on $\mu$; for $k>2/3$ the stable region is $k<\mu <k^2/(3k-2)$.

At the lower boundary $\mu =k$, the positive equilibrium point merges with extinction, with a multiplier $+1$ and slow recovery. At the upper boundary a multiplier reaches $-1$ (the other is $k$), giving alternating adult/larval fluctuations. The nonlinear map confirms a supercritical period-doubling bifurcation, as follows. For $k>2/3$, put $h=1-k$ and ${\mu }_c=k^2/(3k-2)$. A nonconstant period-two adult sequence $A,B$ satisfies $B=hA+\mu B/(1+A)$ and $A=hB+\mu A/(1+B)$. Subtracting and adding these equations givesThus positive unequal values emerge for $\mu >{\mu }_c$, with amplitude proportional to $\sqrt{\mu -{\mu }_c}$. The two larval phases are $\mu B,\mu A$. If $T_2,{\Delta }_2$ are the trace and determinant of the two-step Jacobian matrix, direct multiplication givesAt onset ${\Delta }_2=k^2<1$ and $T_2=1+k^2$, so the other two strict Jury stability criterion inequalities hold for sufficiently small positive $\mu -{\mu }_c$. Hence just beyond the upper boundary, a stable period-two population cycle replaces the stable equilibrium. This is a local conclusion, not stability for arbitrarily large reproduction rates.

At the upper boundary itself, the positive equilibrium point is still locally asymptotically stable, despite its multiplier $-1$. To decide this equality case, apply the centre manifold theorem for a discrete dynamical system. Put $u=a-a_{\ast }$ and write the centre graph as $b-b_{\ast }=H(u)$ with $H^{\prime }(0)=-{\mu }_c$. If the reduced map is $f(u)=-u+c_2{u}^2+c_3{u}^3+O(u^4)$, its invariance equation is $H(f(u))={\mu }_{c}u$. Expanding the original recruitment map and this equation through cubic order givesConsequently $f^2(u)=u-2(c_3+c_2^2)u^3+O(u^4)$. The stability at a nondegenerate flip bifurcation criterion shows algebraic attraction on the centre direction; the transverse multiplier $k$ has modulus less than one. Including this boundary, the complete local asymptotic-stability condition for a positive population is thereforeThe shaded figure shows the strict linear-stability region; its upper boundary adds this nonhyperbolic attracting case.

Demographic stochasticity near the lower boundary can lead to absorption at extinction; near the upper boundary it excites alternating fluctuations and can blur the deterministic period-doubling bifurcation. These qualitative predictions depend on the chosen stochastic recruitment and mortality rules.

Multiplicativity of the determinant gives $\delta (gh)=\delta (g)\delta (h)$ and $\delta (1)=1$, so $\delta$ is a linear character. Its finite image is a cyclic group of roots of unity. If it contains $-1$, its order is $2m$ and the character of a representation ${\delta }^m$ maps onto $\{\pm 1\}$. The kernel of a group homomorphism is the preimage of its identity, so $\ker ({\delta }^m)$ is a normal subgroup of index $2$. The subgroup $\ker \delta$ itself need not have index $2$.