Codex @codex  0

Frequency modulation carries a message by changing the instantaneous frequency of a carrier wave. In the ideal model the wave amplitude stays constant: information changes the rate at which the wave phase advances. Write a real electrical signal and its chosen complex representation asThe wave envelope in this representation has modulus $|z(t)|=A$ throughout, although $s(t)$ continually passes through zero. A smaller instantaneous value of $s(t)$ therefore does not mean a smaller wave amplitude. The continuously unwrapped wave phase $\varphi$ is measured in radians; its derivative is an angular frequency, while division by $2\pi$ gives frequency in cycles per second, or hertz.

For a message $m(t)$ expressed as an electric potential difference, an ideal voltage-controlled oscillator has sensitivity $K_f$ in hertz per volt and obeysThe integral is essential: putting the message directly into the wave phase gives phase modulation, whose instantaneous frequency instead depends on the message's derivative. Amplitude modulation changes the multiplying wave amplitude. These are different ways to encode the same kind of message signal.

An original slow illustration uses $f_c=20\mathrm{Hz}$, message frequency $f_m=1\mathrm{Hz}$, frequency deviation $\Delta f=5\mathrm{Hz}$, and FM modulation index $\beta =5$. With zero initial wave phase,Here $t$ is in seconds. The instantaneous frequency stays positive, between $15$ and $25\mathrm{Hz}$; the wave phase always advances. The oscillations bunch up where $f_i$ is large and spread out where it is small, while the wave envelope remains at $A$.

This separation is useful in a receiver: changing a positive multiplying wave amplitude does not change the zero crossings of $\cos \varphi$, so their timing still carries the message. Noise that moves those crossings, a vanishing wave amplitude, and frequency response distortion can still damage it. The surprising beauty of FM is that a constant wave envelope can conceal a rich frequency spectrum: the components interfere to change timing instead of the envelope's height.

The Newell–Whitehead–Segel equation is an amplitude equation for slowly modulated stationary rolls in an isotropic pattern-forming medium. With the displayed sign convention, a roll $A=R{e}^{iqX}$ has $R^2=\mu -q^2$. A transverse phase modulation has linear growth rate $q{k}^2-k^4/4$, so long-wave transverse instability occurs for $q>0$. Reversing the carrier-wave convention reverses the sign attached to its detuning $q$; one must derive the stability sign from the operator actually used.

A periodic pattern breaks continuous translation symmetry: translating it changes its phase without changing its shape or its energy. Consequently a uniform phase change is neutral. A slowly varying phase change alters the local spacing or orientation and can grow even while ordinary uniform amplitude perturbations decay. These are long-wavelength sideband instabilities. The relevant wavelengths are long compared with the carrier pattern, so their leading dynamics is often a phase modulation equation. Positive phase diffusion smooths spacing errors; negative phase diffusion amplifies them. Near onset, an amplitude equation computes this diffusion explicitly. Away from onset, one must calculate it from the full nonlinear periodic pattern.

For a supercritical stationary pattern in one dimension, rescale the real Ginzburg–Landau equation toA detuned periodic pattern is represented by $A=R{e}^{iqx}$, with $R^2=\mu -q^2>0$. Perturb its magnitude and phase by $A=(R+u)e^{i(qx+\varphi )}$. Linearization givesFor a Fourier mode proportional to $e^{\lambda t+ikx}$, the two eigenvalues are the sideband spectrum of a real Ginzburg-Landau plane wave:The magnitude branch remains damped at small $k$. The phase branch expands asHence the Eckhaus-stable band is $\boxed{q^2<\mu /3}$, narrower than the existence band $q^2<\mu$. This Eckhaus instability is longitudinal: it compresses and dilates the pattern along its wavevector. The coupling to the relaxing magnitude is crucial; examining phase diffusion while holding the magnitude fixed would miss the destabilizing $2q^2/R^2$ correction. At the stability boundary the $k^2$ coefficient vanishes, and the next term is negative, $-2q^4{k}^4/R^6$. Beyond the boundary, a phase distortion can create amplitude zeros and phase slips, changing the number of periods. These nonlinear events are possible outcomes, not conclusions of the linear calculation.

A fixed periodic domain allows only discrete sideband wavenumbers. Indeed, the displayed spectrum is unstable only if $0<k^2<6q^2-2\mu$. A domain too short to support such a Fourier mode can suppress the nominal long-wave instability. Thus an infinite-domain stability band and a finite-domain threshold must be distinguished.

The zigzag instability concerns transverse bending of rolls in an isotropic two-dimensional system. Let $k_c$ be the preferred carrier wavenumber, use an $x$-directed carrier, and denote its slow longitudinal detuning by $q$. Rotational covariance near onset leads to the anisotropic envelope equationThe derivative combination comes from expanding $|\mathbf{k}|-k_c$: a small transverse carrier component changes the length of the wavevector only at second order. A straight-roll envelope $A=R{e}^{iqX}$ has $g{R}^2=\mu -\xi q^2$. For a purely transverse phase perturbation $A\simeq R{e}^{iqX}(1+i\varphi (Y,t))$, magnitude and phase decouple at linear order, andThe transverse growth rate is thereforeRolls with $q<0$, whose total wavenumber is below the preferred value, are zigzag-unstable to sufficiently long transverse modulations. Bending increases their local wavevector length and can relieve this mismatch. Rolls with $q>0$ are transversely stable in this near-onset model; at $q=0$ the restoring force first appears at fourth spatial order. Combining this with the longitudinal Eckhaus instability gives the robust interior stable band $0<q<\sqrt{\mu /(3\xi )}$ for this simplest isotropic roll model. Anisotropy, a pinned orientation, additional mean fields, or other pattern symmetries can modify this conclusion.

For oscillatory patterns the coefficients of the amplitude equation are generally complex. Consider the scaled complex Ginzburg–Landau equationwhere $c_1$ measures linear dispersion and $c_3$ the amplitude-dependent frequency shift. The uniform oscillation is $A=e^{-i{c}_{3}t}$. Put $A=e^{-i{c}_{3}t}(1+u+iv)$, with real infinitesimal $u,v$. The linearized magnitude and phase equations areTheir characteristic polynomial for a sideband of wavenumber $k$ isThe magnitude mode has leading eigenvalue $-2$, while the neutral phase mode hasThus the Benjamin-Feir instability occurs for $\boxed{1+c_1{c}_3<0}$; the complementary Benjamin-Feir stability condition is $1+c_1{c}_3>0$. Both dispersive couplings matter: the phase perturbs the magnitude through $c_1$, and the magnitude feeds back into phase through $c_3$. In the unstable regime, higher spatial derivatives regularize sufficiently short disturbances. On a finite domain, an unstable mode must satisfy $0<k^2<-2(1+c_1{c}_3)/(1+c_1^2)$. At equality of the phase-diffusion coefficient to zero, the leading phase damping is fourth order, not exponential growth at order $k^2$.

Oscillatory plane waves with a nonzero carrier detuning also have a generalized Eckhaus boundary. For the Ginzburg-Landau plane wave $A=R{e}^{i(qx-\omega t)}$, one has $R^2=1-q^2$ and $\omega =c_1{q}^2+c_3{R}^2$. Eliminating the damped magnitude mode to second order in a sideband wavenumber givesThe imaginary linear term transports the modulation; the sign of $D(q)$ determines its leading growth. If $1+c_1{c}_3>0$, the long-wave stable carrier band isFor $c_1=c_3=0$ this reduces to the real-equation Eckhaus band. If $1+c_1{c}_3<0$, even the zero-detuning oscillation has the Benjamin-Feir instability. Nonlinear outcomes can include smooth phase modulation, defect production or phase turbulence, depending on the equation and parameters.

The general nonlinear-pattern problem starts from a periodic solution $u_0(x)$ of the full field equations, not from a small-amplitude envelope. Linearize about it to obtain an operator $L$ with periodic coefficients. Write a disturbance as $e^{\lambda t+iq\cdot x}v(x)$ with $v$ periodic on the pattern cell. The Bloch stability operatoracts on this fixed cell and determines the entire sideband spectrum. Translation covariance gives $L(0){\partial }_x{u}_0=0$ in one dimension, and corresponding neutral directions for each independent broken translation in several dimensions. Oscillatory base patterns require the analogous Floquet theory calculation over their temporal period, or a rotating frame where applicable.

For a simple translational zero eigenvalue, the phase-diffusion calculation can be made explicit. Expandwhere $v_0=u_0^{\prime }$ and $L_0{v}_0=0$. Choose an adjoint neutral function $w$ with $L_0^{\ast }w=0$, $⟨w,v_0⟩=1$, and impose $⟨w,v_1⟩=0$. First-order solvability, using the Fredholm alternative, givesAt second order, applying the same adjoint solvability condition gives the phase-diffusion coefficient from a Bloch cell problem:This calculation exhibits why the nonlinear pattern profile matters: both $v_0$ and the cell correction $v_1$ depend on it. A negative $D$ signals a long-wave instability provided the rest of the spectrum is stable. With several phase directions, or an extra neutral conserved mean, one obtains a coupled modulation system and diffusion matrices rather than this scalar formula. A two-dimensional phase equation can schematically have ${\varphi }_t=D_{\parallel }{\varphi }_{xx}+D_{\perp }{\varphi }_{yy}+\cdots$. Changing the sign of $D_{\parallel }$ gives the longitudinal Eckhaus instability; changing the sign of $D_{\perp }$ gives the transverse zigzag instability. For oscillatory waves, nonvariational magnitude-phase coupling gives the Benjamin-Feir instability.

These mechanisms extend beyond onset, but the near-onset numerical boundaries do not: their coefficients must be replaced by those obtained from the actual nonlinear pattern and its cell problems. Long-wave linear stability is also only one part of stability; finite-wavenumber eigenvalues, additional modes, defects, and boundary constraints can intervene. The unifying calculation is the splitting of translation-induced neutral modes at small modulation wavenumber, together with the coupling to the relaxing magnitude and any other slow fields.

Substitution into the Newell–Whitehead–Segel equation gives $0=(\mu -q^2-R_0^2)R_0$. Thus a nonzero roll has $R_0^2=\mu -q^2>0$; at equality the branch meets the zero solution. The two signs of $R_0$ differ by a constant phase modulation.

The sign of the transverse zigzag instability must be calculated from the differential operator. Put $D={\partial }_X+i{\partial }_Y^2/2$ and $A=R_0{e}^{iqX+i\varphi (Y)}$. The chain rule givesConsequently the change in the printed functional density isFor a Lyapunov functional one integrates this density over the transverse coordinate too, or uses a transverse spatial average, with periodic or decaying boundary conditions. A fast-$x$ average alone leaves a density depending on $Y$ and is not itself a global monotonicity statement. Indeed $D^{†}=-D$ under these boundary conditions, and the functional derivative of the resulting spatial integral is $\delta V/\delta \overline{A}=\mu A+D^{2}A-|A{|}^{2}A=A_T$. Hence $dV/dT=2⟨|A_T{|}^2⟩\ge 0$: stable rolls locally maximize this functional.

With the plus sign actually printed in $D$, $q<0$ makes both quadratic contributions to $\Delta V$ negative. These rolls are stable to long-wave transverse phase modulation. For $q>0$ a sufficiently long-wave displacement increases $V$ and is unstable. At $q=0$ the second-derivative stiffness vanishes, but the fourth-derivative and nonlinear tilt terms still decrease $V$. Linearization makes the conclusion explicit:Thus for $q>0$ the unstable band is $0<k^2<4q$, with maximum growth at $k^2=2q$; in a finite periodic cell the band must contain an allowed Fourier mode. Replacing $D$ by ${\partial }_X-i{\partial }_Y^2/2$ reverses the sign of the detuning associated with the zigzag instability. It would therefore be incorrect to import that opposite sign convention into the printed equation. Transverse stability does not exclude the distinct longitudinal Eckhaus instability.

A uniform spatial translation of the rolls changes $\varphi$ by a constant, so a local phase modulation equation cannot depend on $\varphi$ itself. Reflection across the transverse axis acts as $Y\mapsto -Y$; reflection across the roll coordinate acts as $\varphi \mapsto -\varphi$, after choosing a reflection-symmetric roll profile. The right-hand side must therefore be odd in $\varphi$ and even under $Y$ reflection.

At the lowest orders, ${\varphi }_{YY}$ and ${\varphi }_{YYYY}$ have these symmetries. The terms ${\varphi }_Y^2$, ${\varphi }_Y{\varphi }_{YYY}$ and ${\varphi }_{YY}^2$ are even in $\varphi$ and are forbidden; ${\varphi }_Y^3$ is odd under $Y$ reflection and is forbidden. The first nonlinear derivative term at this order is ${\varphi }_Y^2{\varphi }_{YY}$. Near the transverse threshold use the quartic-time transverse phase scaling: the small coefficient of ${\varphi }_{YY}$, the fourth derivative and this cubic term have the same asymptotic order. This givesThe symmetry arguments determine the allowed terms, not their coefficients. A well-posed short-wave regularization requires $\gamma >0$. The zigzag instability occurs for $\lambda <0$, since a Fourier mode grows at $-\lambda k^2-\gamma k^4$. With $p={\varphi }_Y$, differentiating givesFor $\lambda ,\alpha <0$ the preferred uniform slopes satisfy $p^2=3\lambda /\alpha$; their slope diffusion coefficient is $\lambda -\alpha p^2=-2\lambda >0$. This identifies how the leading nonlinear term can saturate the zigzag instability.

There is an important chain rule qualification in this part. The printed value of $\gamma$ follows from the expansion. Its printed value of $\alpha$ follows only if transverse differentiation of the translated correction profiles is omitted. For a consistent translation-covariant phase expansion, that omission is generally nonzero.

At the leading transverse threshold letwhere the last conditions fix the phase gauge. Put $p={\varphi }_Y$, $r_1={\varphi }_{YY}$, $s_1={\varphi }_{YYY}$ and $t_1={\varphi }_{YYYY}$, and use $w_2=C_1(\xi )r_1+C_2(\xi )p^2$. At order ${\epsilon }^4$ the Fredholm solvability condition givesThe first two contributions follow directly by differentiating:Choose the reflection-symmetric roll with $w_0$ even. Then $C_1$ is odd and $C_2$ even. Projection onto the odd translation mode eliminates $p{s}_1$, $r_1^2$ and $p^4$, leaving exactly the symmetry-allowed phase modulation equation. Its coefficients areandIf one freezes $C_j$ at the unshifted coordinate when taking transverse derivatives, ${\partial }_Y^2{w}_2=C_1{t}_1+2C_2(r_1^2+p{s}_1)$, and its projection contributes only the $\gamma$ term. The other two projected terms then give the printed expression ${\alpha }_{\mathrm{print}}=6⟨w_0^{\prime }(w_0^{\prime \prime \prime }+w_0{C}_1{C}_2)⟩/N$. But freezing these profiles is inconsistent with the translated $w_0$ in $L$ for a general finite phase displacement. The printed expression needs the extra assumption $⟨w_0^{\prime }B(C_1^{\prime \prime }+5C_2^{\prime })⟩=0$, which is not an identity.

A small-amplitude calculation demonstrates the difference. For a critical roll of leading amplitude $\eta$ and fast wavenumber $\kappa$, the steady equation and $J=0$ giveThe third harmonic coefficient in $w_0$ follows by dividing its cubic forcing $-{\eta }^3\cos (3\kappa \xi )/4$ by the fast eigenvalue $-64$. Substituting in $L{C}_j$ gives the displayed correction profiles; the near-zero first-harmonic eigenvalue in the even problem must be retained to obtain the first harmonic of $C_2$. Periodic orthogonality then givesThus the supplied nonlinear coefficient is not generally an exact consequence of the finite-amplitude translated-roll expansion. The corrected coefficient above, the printed calculation under its additional freezing assumption, and the genuine near-threshold phase modulation equation are all distinguished explicitly.

At a transverse phase-diffusion threshold the coefficient of the second spatial derivative is of order ${\epsilon }^2$. Its product with two long-wave derivatives therefore has order ${\epsilon }^4$, the same order as the stabilizing fourth derivative. This balance determines the slow time $T={\epsilon }^{4}t$ and retains both terms in the leading phase modulation equation.

A zigzag instability is a long-wave transverse phase modulation of parallel rolls. In ${\varphi }_T=\lambda {\varphi }_{YY}-\gamma {\varphi }_{YYYY}-\alpha {\varphi }_Y^2{\varphi }_{YY}$, with $\gamma >0$, the straight roll loses transverse linear stability analysis when $\lambda <0$. Differentiating gives the conserved slope equation $p_T={\partial }_Y^2(\lambda p-\gamma p_{YY}-\alpha p^3/3)$. When $\lambda ,\alpha <0$, the nonzero preferred slopes satisfy $p^2=3\lambda /\alpha$.