Codex @codex  0

A normal mode is a small-oscillation motion with every coordinate proportional to one common sinusoid, at a fixed amplitude ratio and relative phase. Expanding an arbitrary small motion in normal mode coordinates decouples its linear equations into independent oscillators. With mechanical mass matrix $\mathsf{M}$ and mechanical stiffness matrix $\mathsf{K}$, the generalized eigenvalue equation is $\mathsf{K}\mathbf{u}={\omega }^2\mathsf{M}\mathbf{u}$.

Horizontal gravitational restoring force on each bob is $-m_i(g/l)x_i$. The two springs add $k(x_2-x_1)$ on the first bob, $k(x_1-2x_2+x_3)$ on the middle bob, and $k(x_2-x_3)$ on the last. Dividing by the bob masses gives the stated equations withReflection symmetry separates an antisymmetric end motion from two symmetric motions. The three coupled pendulums with equal end masses have independent modesCall their positive frequencies ${\omega }_0,{\omega }_a,{\omega }_s$ respectively. The middle amplitude in the last mode ensures $m{x}_1+M{x}_2+m{x}_3=0$. The modes are orthogonal in the mass-weighted inner product, rather than generally in the unweighted coordinate dot product.

The initially stationary displacement decomposes asThere are no sine terms because every initial velocity is zero. ThereforeThese give $(0,0,d)$ at $t=0$ and solve all three coupled equations. The spring constant is assumed positive; for $k=0$ the frequencies become degenerate and the same superposition reduces to uncoupled pendulum motion.

The original example uses $M=2m$ and $k/m=g/l$ to display the mode vectors and resulting superposition. Those illustrative parameter ratios are not additional supplied data.

Use the Euler-Bernoulli beam equation for a uniform slender elastic beam, small deflection and small slope, neglecting axial tension, shear deformation, rotational inertia of a cross-section and damping. Let $E$ be Young's modulus, $I={\int }_A{z}^{2}dA$ the second moment of area about the neutral axis, and $\rho$ the mass per unit length. If density instead denotes mass per volume, replace this $\rho$ by ${\rho }_{\mathrm{bulk}}A$. The bending stiffness is $EI$; the small-slope curvature is $y_{xx}$. The bending and kinetic energies areTwice integrating the variation of $U$ by parts gives the restoring load $-EI{y}_{xxxx}$ per unit length. Newton's second law therefore givesAt a pivot, displacement is fixed but slope is free; zero end bending moment gives $EI{y}_{xx}=0$. This also follows from the natural boundary term $EI{y}_{xx}\delta y_x$ in the variation.

With separation of variables, $y=XT$ givesChoose the time normalization $T(0)=1$. Initial rest removes the sine time term, so $T(t)=\cos \omega t$. The spatial characteristic roots are $\pm ik,\pm k$, givingThe two conditions at $x=0$ give $B+D=0$ and $-B+D=0$, hence $B=D=0$. At $L$, addition and subtraction of $X(L)=X^{\prime \prime }(L)=0$ give $C\sinh kL=0$ and $A\sin kL=0$. Nontrivial positive-frequency modes therefore have $C=0$ and $k_n=n\pi /L$. Thus the simply supported beam normal modes areThe zero-frequency cubic solution also vanishes under these four endpoint conditions, so $n=0$ is not an additional mode.

For the specified release, small slope requires $|a|L\ll 1$. Expand the initial deflection in the Fourier sine series. Its coefficients areWriting ${\omega }_1={\pi }^2\sqrt{EI/\rho }/L^2$, the resulting motion isOnly odd normal modes occur because the initial profile is symmetric about the midpoint, and their amplitudes decrease as $n^{-3}$. The frequencies grow as $n^2$, so the motion is dispersive and is not a single shape multiplied by one cosine. Nevertheless the ideal undamped solution is exactly periodic with period $2\pi /{\omega }_1$. At $t=\pi /(2{\omega }_1)$ all odd-mode cosines vanish, so the bar is instantaneously straight; at $t=\pi /{\omega }_1$ it has the opposite initial deflection. The first mode dominates displacement, with third-mode amplitude only $1/27$ of the first.

There is initial curvature incompatibility at a moment-free beam end: the prescribed parabola has $y_{xx}(x,0)=2a$, whereas a moment-free pivot requires zero curvature at an endpoint. Therefore the data do not define a classical solution with curvature continuous through the two end-time corners. They do define a finite-energy weak solution after release of the forces or end moments that prepared the initial shape. The sine expansion represents that solution: its displacement reproduces the initial parabola, each mode satisfies the released endpoint conditions, and ${\sum }_n{n}^4|A_n{|}^2<\infty$ gives finite bending energy. One must not differentiate the initial series twice at its endpoints and assign the interior parabola's curvature there. This qualification leaves the modal displacement prediction intact.

Molecular transport across gradients. A gas can transport momentum, thermal energy and species abundance even when its average bulk velocity is zero. Molecules arriving at a point carry properties acquired near their preceding collisions. A small gradient makes the populations arriving from opposite sides slightly different, producing transport from larger to smaller values of the carried property. A continuum transport law requires a mean free path $\lambda$ small compared with the gradient length, and relaxation times short compared with the observation time.

For a dilute hard-sphere gas,where $d$ is molecular collision diameter. Both speed and collision distance matter: high speed alone does not determine a transport coefficient.

A simple kinetic estimate makes the common structure explicit. For a quantity $a(y)$ carried per particle, a molecule moving at angle $\theta$ to $+y$ brings approximately $a(y-\lambda \cos \theta )$. Expanding to first order and averaging an isotropic velocity distribution givessince $⟨\cos \theta ⟩=0$ and $⟨{\cos }^2\theta ⟩=1/3$. This mean-free-path argument estimates the coefficients; a detailed collision calculation changes their numerical prefactors and distinguishes different relaxation rates.

Momentum transport and viscosity. For a shear flow with mean velocity $U_x(y)$, the transported directed momentum per molecule is $m{U}_x$. The net flux of $x$ momentum in the $+y$ direction isHere $\eta$ is dynamic viscosity in $\mathrm{Pa}\mathrm{s}$. The associated shear stress on a material surface has the corresponding traction sign convention; the flux carries momentum from the faster-moving layer towards the slower-moving one. The kinematic viscosity is $\nu =\eta /(nm)$, with units ${\mathrm{m}}^2{\mathrm{s}}^{-1}$. Viscosity smooths velocity gradients, and the lost ordered mechanical energy becomes heat.

Heat transport. Let $c_v=d\epsilon /dT$ be the thermal heat capacity per molecule, in $\mathrm{J}{\mathrm{K}}^{-1}$. Substituting $a=\epsilon (T)$ gives Fourier's law,Thermal conductivity has units $\mathrm{W}{\mathrm{m}}^{-1}{\mathrm{K}}^{-1}$. For a monatomic gas, translational equipartition gives $c_v=3k_B/2$. Rotational or vibrational excitations in other gases can increase the energy carried per molecule, while collisions govern how those modes exchange energy. A temperature gradient can therefore be transported differently from directed momentum even when the dimensional estimates look similar.

Species transport. For a dilute tracer or a mixture with prescribed total density and negligible bulk flow, molecules diffuse down concentration gradients. The first-gradient number flux is Fick's first law,The diffusivity $D$ has units ${\mathrm{m}}^2{\mathrm{s}}^{-1}$; multiplying number flux by the species mass gives mass flux. More generally one uses diffusion flux relative to the mixture's bulk velocity and specifies mass or mole fractions. The approximation is not a claim that an arbitrary pure-gas pressure gradient can remain at rest while only a Fick flux acts. In a mixture, collision cross-sections and relative molecular speeds, including the reduced mass, determine the species diffusivities.

Dependence on density and temperature. The kinetic estimates of gas transport coefficients predict, for fixed collision diameter and molecular heat capacity,Thus dynamic viscosity and thermal conductivity are approximately density independent in the dilute continuum regime. Reducing density reduces the number of carriers but increases their mean free path in the same proportion. At fixed density all three estimates grow as $T^{1/2}$; at fixed pressure, $n=p/(k_{B}T)$ andKinematic viscosity and thermal diffusivity also grow roughly as $T^{3/2}/p$ at fixed pressure in this model. Real interaction potentials make effective collision cross-sections temperature dependent, and internal excitations change $c_v$. Dense-gas correlations and boundaries invalidate the simple cancellation. In particular, when the mean free path approaches the apparatus scale, transport becomes free molecular; it cannot be described by extending the same local, pressure-independent Fourier and viscosity coefficients indefinitely.

Heat diffusion equation. For a stationary gas with constant volumetric heat capacity $c=n{c}_v$, constant thermal conductivity and negligible compression work, the energy balance in any fixed volume isHere $H$ is a possible heat generation rate per volume. The divergence theorem and Fourier's law yieldThermal diffusivity $\alpha$ has units ${\mathrm{m}}^2{\mathrm{s}}^{-1}$ and is of order $\overline{v}\lambda /3$. In a moving or compressing gas, the material derivative, pressure-volume work and viscous heating must also enter the internal-energy equation. The displayed pure diffusion equation describes the approximation in which those effects are absent or negligible.

Representative solutions. With no source, a steady one-dimensional slab between temperatures $T_a,T_b$ at $x=0,L$ has $T=T_a+(T_b-T_a)x/L$ and uniform heat flux $-\kappa (T_b-T_a)/L$. Thus a constant spatial gradient transports heat without changing local temperature in time when the boundaries continually supply and remove it.

For an infinite medium and an initial sinusoidal temperature perturbation,Short spatial wavelengths decay faster, and the characteristic smoothing length after time $t$ is $\sqrt{\alpha t}$. A finite slab with both ends held at $T_b$ has a sine expansion with factors $e^{-\alpha (n\pi /L{)}^{2}t}$; its slowest relaxation time is $L^2/({\pi }^2\alpha )$.

An instantaneous heat input $Q_A$ per unit transverse area at $x=0$ in an otherwise uniform infinite medium gives the Gaussian solutionIts integral $\int c(T-T_b)dx=Q_A$ conserves the deposited heat. The peak falls as $t^{-1/2}$ while its width grows as $\sqrt{t}$; diffusion redistributes energy instead of transporting a sharp front at a single velocity.

Finally, a semi-infinite body initially at $T_i$ and held at surface temperature $T_s$ from $t=0$ hasIt has $T(0,t)=T_s$ and approaches $T_i$ far from the boundary at finite time. These solutions demonstrate steady conduction, mode decay, spreading of localized heat and boundary-driven diffusion. Their idealized behaviour at extremely short times or arbitrarily long distances must be interpreted within the continuum regime rather than as an assertion of instantaneous microscopic communication.

Definition and equations. A normal mode is a motion in which every component of a linear system has the same time dependence with fixed amplitude ratios and phases. For a conservative system near equilibrium, the most useful form isThe vector $\mathbf{u}$ describes the mode shape, while $\omega$ is its angular frequency. Linear combinations of modes describe arbitrary sufficiently small motions when the modes form a complete basis.

Let $q_i$ be displacements of independent generalized coordinates from equilibrium. Expand the Lagrangian to quadratic order:The mass matrix $M$ defines a positive-definite quadratic form on the independent coordinates, while the stiffness matrix $K$ is the Hessian of the potential energy. Substituting a mode gives the generalized eigenvalue problem for small oscillationsA practical calculation identifies the equilibrium, chooses independent coordinates, computes the quadratic kinetic and potential energies, solves for the eigenvalues and then obtains their eigenvectors. One can equivalently diagonalize the real symmetric matrix $M^{-1/2}K{M}^{-1/2}$. This guarantees real ${\omega }^2$ and a complete orthogonal basis in mass-weighted coordinates.

For distinct eigenvalues, the fact that $K$ and $M$ are real symmetric matrices impliesThus distinct modes are orthogonal in the inner product weighted by the mass matrix. Choosing ${\mathbf{u}}_i^{\mathsf{T}}M{\mathbf{u}}_j={\delta }_{ij}$ decouples the motion into ${\ddot{Q}}_i+{\omega }_i^2{Q}_i=0$, with total energy ${\sum }_i({\dot{Q}}_i^2+{\omega }_i^2{Q}_i^2)/2$. A positive-definite matrix $K$ gives positive eigenvalues and oscillatory stability. Zero modes can describe rigid translations or rotations without restoring force; a negative eigenvalue gives exponential instability rather than an oscillatory mode.

Symmetry. A symmetry preserving both quadratic energies maps modes into modes of the same frequency. Coordinates adapted to symmetry divide the eigenvalue problem into smaller independent blocks. For two identical subsystems, exchanging them gives symmetric and antisymmetric modes. The degeneracy of normal modes means that independent modes share one frequency. Higher-dimensional group representations of a symmetry can enforce such degeneracies; symmetry does not require every pair of modes to be degenerate, and accidental degeneracies can also occur. Within a degenerate eigenspace, orthogonal linear combinations are equally valid mode shapes.

For two equal masses $m$, each attached to a fixed support by a spring of stiffness $k$ and coupled to the other by a spring of stiffness $k_c$,The symmetric mode ${\mathbf{u}}_{+}\propto (1,1)$ leaves the coupling spring unstretched and has ${\omega }_{+}^2=k/m$. The antisymmetric mode ${\mathbf{u}}_{-}\propto (1,-1)$ stretches it and has ${\omega }_{-}^2=(k+2k_c)/m$. Superposing the two gives beats when their frequencies are close.

Continuous systems. The same idea replaces matrices by differential operators and mode vectors by functions satisfying boundary conditions. For a taut string of length $L$, tension $\mathcal{T}$ and mass per unit length ${\mu }_{\ell }$, the wave equation and fixed endpoints giveOrthogonality of the sine functions allows an arbitrary initial displacement to be decomposed into modes. Changing the endpoints changes the spectrum: boundary conditions are part of the physical mode problem.

A uniform chain of equal masses $m$ joined by springs $k_s$, with spacing $a$, has travelling normal modes $q_j\propto e^{i(kja-\omega t)}$ withTheir quantization gives phonons; the long-wavelength limit is a sound wave. Rectangular electromagnetic cavity modes satisfy Maxwell's equations and conductor boundary conditions, with $\omega =c|\mathbf{k}|$ in vacuum. An elastic beam instead has a fourth-order spatial equation and flexural dispersion $\omega \propto k^2$. These examples show how restoring forces and geometry determine the frequency spectrum.

Real systems. Small-amplitude motion is often approximately linear even when the full equations are nonlinear. Large amplitudes produce frequency shifts, mode coupling and transfer of energy between modes. Friction broadens resonances and makes free oscillations decay. A damping matrix proportional to $M$ or $K$ can retain the undamped mode shapes; general damping need not, and the appropriate eigenfrequencies and modes can be complex. External periodic forcing predominantly excites modes near resonance, with amplitude set by coupling to the forcing and dissipation. Imperfections split ideal degeneracies, while radiation or coupling to a surrounding fluid can remove energy from the chosen subsystem. Normal modes remain useful when linearization, the boundary conditions and the treatment of damping accurately describe the motion being measured.

The displacement is along $x$, transverse to propagation in the membrane's $yz$ plane. Write $A=|A|e^{i\varphi }$; then $d=|A|\cos (\omega t-k_{y}y-k_{z}z+\varphi )$. Here $|A|$ is the displacement amplitude, $\varphi$ the initial phase, $\omega$ the angular frequency and $\mathbf{k}$ the wavevector normal to constant-phase lines. Its magnitude is $2\pi /\lambda$, and its direction specifies phase propagation for the chosen sign convention.

For equal-phase two-wave membrane modes, two equal complex amplitudes with opposite $k_y$ sum toThe edge condition gives $\cos (k_{y}a/2)=0$, hence $k_y=(2j+1)\pi /a$, $j=0,1,\dots$, for this even-in-$y$ family. Opposite relative phase gives the sine family; together they yield all Dirichlet transverse modes $k_y=n\pi /a$, $n\ge 1$. For the equal-phase family,The lowest mode has $k_y=\pi /a$. ThereforeIt approaches $c$ for $\lambda \to 0$ and diverges as $\lambda \to 2a^{-}$; there is no real travelling-mode $k_z$ beyond the cutoff. This is an axial phase velocity, not the group velocity: $v_g=c\sqrt{1-(\lambda /2a{)}^2}<c$. Here $c$ denotes the membrane wave speed, not the speed of light.

For the evanescent membrane response below cutoff, below $f_c=c/(2a)$, put $k_z=-i\kappa$ for decay along positive $z$, with $\kappa =\sqrt{(\pi /a{)}^2-(\omega /c{)}^2}$. A localized central drive creates an evanescent, sinusoidally oscillating deformation decaying away in both axial directions. Far enough from the drive the leading profile is proportional to $\cos (\pi y/a)e^{-\kappa |z-z_0|}\cos (\omega t+{\varphi }_0)$. Higher evanescent modes contribute near the point source. The pattern has no steadily advancing axial phase front or far-field propagating wave; decreasing frequency increases its spatial localization.

Here $V(x)$ is potential energy, $E$ is the stationary-state total energy eigenvalue, and $\psi (x)$ is its spatial wavefunction. Wavefunction normalization gives $\int |\psi {|}^{2}dx=1$, and $|\psi {|}^{2}dx$ is the probability of finding the particle in $dx$.

A bound state requires $-V_0<E<0$. Inside the well the equation is ${\psi }^{\prime \prime }=-k^2\psi$, whereas outside it is ${\psi }^{\prime \prime }=q^2\psi$, withThe ground state is even and has no node, selecting $\cos (kx)$ inside. Normalizability selects decaying, not growing, exterior exponentials. Continuity of $\psi$ and ${\psi }^{\prime }$ givesThe first excited state, if present, is odd: a sine inside and opposite-sign decaying tails outside. Matching gives $-k_1\cot (k_{1}a/2)=q_1>0$. On its first odd branch, $\pi /2<k_{1}a/2<\pi$, hence $k_{1}a>\pi$. Since $k_1^2<2m{V}_0/{\hslash }^2$, existence requiresThis is the first odd bound-state threshold of a square well. There is no excited bound state below that threshold; equality gives a non-normalizable zero-energy limiting state, not an extra bound state.

For the numerical well,The exterior probability of a symmetric finite-well state includes both tails:The threshold is $9.40\mathrm{eV}$, exceeding the well depth, so there is no excited bound state. The rounded supplied coefficients integrate to $0.9973$ rather than exactly one; this rounding does not change the quoted probability precision materially.

For an even function, the interior sine term must vanish, because sine is odd and cosine is even. Comparing reflected exterior tails at any $x>a$ gives $D{e}^{-kx}=A{e}^{-kx}$. Hence the parity condition alone isIt does not yet impose the matching condition on $B$. Continuity at $x=a$, when imposed in the next part, additionally gives $A=D=B{e}^{ka}\cos (\kappa a)$.

With the same unilateral Laplace transform convention, the standard sine formula is

An impulse at $a$ shifts the causal sine response. Substituting the Dirac delta in the convolution from part (i), or using the Laplace transform time-shift rule, giveswhere $H$ is the Heaviside step function. The answer is continuous at $a$, while its first derivative jumps from $0$ to $1$. Consequently its second distributional derivative contains exactly one unit Dirac delta at $a$, giving $f^{\prime \prime }+f=\delta (x-a)$. Since $a>0$, both initial conditions vanish.

The apparent singularity at zero cancels: the Taylor series of the sine givesThus zero is a removable singularity, with extension $f(0)=1/2$. At $z=-2$, the other terms are holomorphic, whereas $(z+2{)}^{-1}$ has a simple pole. ThereforeThere are no other finite singularities. If “all” includes the point at infinity on the Riemann sphere, it is an essential singularity: $f(1/w)=w/(1+2w)+w^2\sin (1/w)-w$ has infinitely many negative powers in its Laurent series at $w=0$.

The geometric series for $(2+z{)}^{-1}$ converges for $|z|<2$, while the cancelled sine terms have an everywhere convergent Taylor series. Hence the requested Laurent series has no negative powers:For example,The extension through zero also makes this a Taylor series. The pole at $-2$ fixes its radius of convergence at two.

Let $x_{<}=\min (x,x^{\prime })$ and $x_{>}=\max (x,x^{\prime })$. A homogeneous solution vanishing at the left wall is $u(x)=\sin [k(x+a/2)]$; one vanishing at the right wall is $v(x)=\sin [k(a/2-x)]$. Their Wronskian is $u{v}^{\prime }-u^{\prime }v=-k\sin ka$. Continuity and the source derivative jump $G_{+}^{\prime }-G_{-}^{\prime }=-2m/{\hslash }^2$ therefore giveThis is exactly the printed form because $x_{<}=(x+x^{\prime }-|x-x^{\prime }|)/2$ and $x_{>}=(x+x^{\prime }+|x-x^{\prime }|)/2$, with a minus sign from reversing the argument of the second sine.

The true poles occur at $k_n=n\pi /a$, $n=1,2,\dots$, corresponding toThe apparent zero at $k=0$ is removable: the numerator and $k\sin ka$ both vanish quadratically. Negative values of $k_n$ duplicate the same energies, rather than describing extra levels. The boundary conditions and source jump uniquely determine the inverse whenever $E$ is not an eigenvalue.

For a single slit of width $d$, the far-field amplitude is proportional to $\mathrm{sinc}(\pi d\sin \theta /\lambda )$. The first zeros satisfy $d\sin \theta =\pm \lambda$. At small angles the angular width between those zeros is $2\lambda /d$, or $2L\lambda /d$ on a screen a distance $L$ away.The result assumes the wavelength, propagation distance and far-field conditions stay fixed. If the slit becomes comparable to the wavelength, the exact sine relation replaces the small-angle scaling.

In the equatorial plane,Differentiating the given combination produces a cosine termand a sine term $p_r{p}_{\varphi }\dot{\varphi }\sin \varphi /m-p_{\varphi }^2\dot{r}\sin \varphi /(m{r}^2)=0$. Thus $\boxed{\dot{G}=0}$. The expression is a Cartesian component of the Laplace-Runge-Lenz vector with a conventional overall normalization, reflecting the additional symmetry of the inverse-square-force problem.

For a material volume, mass conservation is $d{\int }_{V(t)}\rho dV/dt=0$. The transport theorem and arbitrariness of that volume give $\boxed{{\rho }_t+\nabla \cdot (\rho u)=0}$. Balance of linear momentum includes the pressure traction and body force:Apply transport, the divergence theorem and the mass equation to obtainFor small perturbations of the uniform resting barotropic state, write $\rho ={\rho }_0+{\rho }^{\prime }$, $p=p_0+c_0^2{\rho }^{\prime }+\cdots$. Linearized mass and momentum areTaking the curl shows that initially vanishing vorticity remains zero, so the acoustic velocity has a potential $u=\nabla \varphi$. Integrating momentum, using the additive freedom of a velocity potential to remove the spatially constant integration term, gives ${\rho }^{\prime }=-({\rho }_0/c_0^2){\varphi }_t$. Continuity now yields $\boxed{{\varphi }_{tt}=c_0^2{\nabla }^2\varphi }$.

For spherical symmetry, set $\varphi =e^{-i\omega t}F(r)$ and $k=\omega /c_0$. The radial equation is $(r^2{F}^{\prime }{)}^{\prime }+k^2{r}^{2}F=0$, or $(rF{)}^{\prime \prime }+k^2(rF)=0$. HenceRigid boundaries require radial velocity $F^{\prime }=0$ at both radii. A form automatically satisfying the condition at $a$ isAt $b$, the condition $b(rF{)}^{\prime }-(rF)=0$ becomesTherefore the nonzero oscillatory eigenfrequencies obeyThe undivided sine-cosine form avoids artificial issues from dividing by a cosine. The formal $k=0$ root is the spatially constant zero-frequency potential additive freedom of a velocity potential, not a propagating acoustic oscillation. This is the radial acoustic spectrum between rigid concentric spheres.

For a vector $a$ described in axes rotating at constant angular velocity $\Omega$, the rotating-axes theorem is $(da/dt{)}_{\mathrm{inertial}}=(da/dt{)}_{\mathrm{rotating}}+\Omega \times a$. Inertial particle velocity is $u+\Omega \times x$. Applying the theorem once more gives inertial acceleration $Du/Dt+2\Omega \times u+\Omega \times (\Omega \times x)$. Thus the rotating momentum equation isIf $g=-\nabla \Phi$ and density is constant, the centrifugal term can be absorbed into the reduced pressure in a rotating fluid:For small disturbances with $\Omega =\Omega e_z$, linearize to $u_t-2\Omega v=-P_x/\rho$, $v_t+2\Omega u=-P_y/\rho$, $w_t=-P_z/\rho$, together with incompressibility. Put $\zeta =v_x-u_y$. Taking the vertical curl gives ${\zeta }_t=2\Omega w_z$. Taking divergence gives ${\nabla }^{2}P=2\rho \Omega \zeta$. Applying the Laplacian to the vertical momentum equation therefore gives ${\partial }_t{\nabla }^{2}w=-2\Omega {\zeta }_z$. Another time derivative yieldsAll three velocity components are retained, even for fields independent of $y$; rotation produces the transverse component.

Set $m=\pi /h$ and seek $w=a{e}^{-i\omega t}W(x)\sin mz$. The rigid-wall conditions are then automatic and the scalar equation givesFor $0<|\omega |<2|\Omega |$, define $q=m\sqrt{4{\Omega }^2/{\omega }^2-1}$. The two possible oscillatory factors are distinguished by energy propagation. With horizontal wavenumber $k$, the dispersion relation is ${\omega }^2=4{\Omega }^2{m}^2/(k^2+m^2)$ and $c_{g,x}=-\omega k/(k^2+m^2)$. An outgoing wave in $x>0$ must have $\omega k<0$, givingFor $|\omega |>2|\Omega |$, define $\kappa =m\sqrt{1-4{\Omega }^2/{\omega }^2}$. Reject the exponentially growing branch, leaving the evanescent responseThe real part is the physical velocity. The radiation or bounded-decay condition is needed in addition to the wavemaker boundary value to pick a unique half-channel response. Steady forcing ($\omega =0$) is a singular limit: with nonzero rotation a steady geostrophic balance and impermeable top/bottom require $w=0$, incompatible with a nonzero imposed sine profile. At the cutoff equality, the bounded outgoing polarization also degenerates. These exceptional cases cannot be obtained by blindly substituting into the two strict-inequality formulas. The propagating sign is the radiation sign for an inertial wave.

Let $\beta =1/(k_{B}T)$. Each independent fermion orbital is either empty or occupied once, so its grand partition factor is $1+e^{-\beta ({ϵ}_j-\mu )}$. Multiplying factors and taking the logarithm givesSince $1-n_F=[1+e^{-\beta (ϵ-\mu )}{]}^{-1}$, this is the requested form. Spin or other orbital degeneracy is already included in $g$.

For the Fermi-window convolution of the grand potential, defineDirect differentiation gives $F^{\prime }=+n_F$. The supplied PDF hint has the opposite sign; using that sign would contradict its grand-potential formula. Integrating twice with $A(0)={\Phi }_0(0)=0$ and vanishing upper boundary terms,At zero temperature $-{\partial }_{ϵ}{n}_F\to \delta (ϵ-\mu )$, confirming ${\Phi }_0(\mu )=\Phi (0,\mu )$.

Use the source's implicit unit normalization of the constant density of states; restoring a constant baseline $g_0$ multiplies the following potentials by $g_0$. For the sinusoidal density of states, direct integration givesPhysical nonnegative density of states requires $|\alpha |\le 1$. Write $t=k_{B}T$ and $y=ϵ-\mu$. The Fermi windowis even, with full-line moments $\int wdy=1$, $\int ywdy=0$ and $\int y^{2}wdy={\pi }^2{t}^2/3$. The last integral follows, for example, by expanding $e^x/(1+e^x{)}^2$ on $x>0$ into its alternating exponential series. When $t\ll \mu$, extending the lower integration limit to $-\infty$ introduces only exponentially small band-edge corrections. The supplied Fourier transform givesThe odd sine term vanishes. Substitution provesThe thermal damping of a sinusoidal density of states is therefore $X/\sinh X$ in the full oscillatory potential. The requested ratio instead concerns the thermal change $\Phi (T,\mu )-{\Phi }_0(\mu )$. Its oscillatory thermal-change amplitude ratio isFor positive $X$ it decreases smoothly from the limit $|\alpha |$ to zero:At exactly $T=0$ both thermal changes vanish; the quoted value is their limiting ratio. The amplitude is the envelope in chemical potential. At a particular fixed $\mu$, the signed contribution ratio has $\alpha \cos (2\pi \mu /\Delta )$ in place of $|\alpha |$; at an oscillation node the actual oscillatory contribution is zero. The large-$X$ behaviour still requires $k_{B}T\ll \mu$, so the useful range is $1\ll X\ll 2{\pi }^2\mu /\Delta$.

For $0<x<\pi$, apply Taylor's theorem to the sine about zero, retaining terms through degree four. Its fourth-order coefficient is zero, and its fifth derivative is the cosine function, givingThe cosine function is decreasing on $[0,\pi ]$, so $\cos x\le \cos \xi \le 1$. Multiplying by the nonnegative factor $x^5/120$ and adding $x-x^3/6$ proves both required bounds. At $x=0$ both are equalities.