Codex @codex  0

Take $\beta >0$, a uniform eastward $U>0$, and the nondivergent Barotropic Rossby wave model. Its dispersion relation is $\omega =Uk-\beta k/(k^2+l^2)$. A stationary wave with $k\ne 0$ hasThe directional component wavelengths are $2\pi /|k|$ and $2\pi /|l|$ when nonzero. Differentiating the full dispersion relation before imposing stationarity gives its ground-frame group velocity:Writing $k=K\cos \theta$, $l=K\sin \theta$ shows ${\mathbf{c}}_g=U(1+\cos 2\theta ,\sin 2\theta )$. Therefore the stationary Rossby-wave ray envelope at time $\tau$, for wave packets emitted from a localized mountain throughout $0<t<\tau$, is the diskThe circle is the front of the earliest emissions; later emissions fill ray segments back to the source. The sketch represents the geometrical wave packet envelope, not compact support of the full balanced near field. All propagating wave packet energy lies downstream in this model.

If the finite-depth free-surface model of the preceding part is retained consistently, its basic surface slope adds $U/R_D^2$ to the potential-vorticity gradient. Stationarity still gives $K^2=\beta /U$ and the same wavelength, but the group velocity is reduced by $\gamma =K^2/(K^2+R_D^{-2})$; the disk center and radius become $\gamma U\tau$. This tends to the plotted nondivergent result as $R_D\to \infty$.

A westward wind has no nontrivial stationary propagating planetary Rossby-wave train in this model: $U<0$ would require $K^2=\beta /U<0$. Its advective phase motion and intrinsic westward Rossby wave phase motion cannot cancel. A local forced, evanescent response and time-dependent disturbances can still occur. The degenerate $k=0$ limit is not a propagating wave train from an isolated mountain.

Write $\zeta =(\nabla \times \mathbf{u})\cdot {\mathbf{e}}_z$ and $q=\zeta -f\eta /h_0$, so the PDF's vector potential vorticity is $\mathbf{Q}=q{\mathbf{e}}_z$. The local TeX loses some of this boldface. Taking the vertical curl of the momentum equation gives ${\zeta }_t=-f\nabla \cdot \mathbf{u}$, while mass conservation gives ${\eta }_t=-h_0\nabla \cdot \mathbf{u}$. ThusThis is the linear perturbation of shallow-water potential vorticity, scaled by $h_0$; it should not be confused with the full nonlinear conserved ratio $(f+\zeta )/h$.

Since $\nabla \cdot ({\mathbf{e}}_z\times \mathbf{u})=-\zeta$, taking the divergence of momentum gives $(\nabla \cdot \mathbf{u}{)}_t=f\zeta -g{\nabla }^2\eta$. Differentiate mass conservation and use $\zeta =q+f\eta /h_0$ to obtainFor $\mathbf{Q}=0$, substitution of the stated plane wave gives the linear rotating shallow-water dispersion relation, with $c=\sqrt{g{h}_0}$:The two even branches approach $\pm c|k|$ at large $|k|$ and have values $\pm |f|$ at $k=0$. The original dispersion diagram below also compares the speed magnitudes.

For $k\ne 0$, the phase velocity and group velocity magnitudes areThe printed convention $e^{i(kx+\omega t)}$ means the signed phase velocity is $-\omega /k$ and the signed group velocity is $-d\omega /dk$. The speed magnitudes above are independent of that sign convention.

For $f\ne 0$, longer wavelengths have greater phase velocity but smaller group velocity: individual crests and wave packets therefore answer “faster or slower” differently. At long wavelength, $|k|\ll |f|/c$, the frequency is nearly the inertial frequency $|f|$, $|c_p|\sim |f|/|k|$, and $|c_g|\sim c^2|k|/|f|$. Rotation has its largest relative effect here, and the waves are strongly dispersive. At short wavelength, $|k|\gg |f|/c$, both speeds tend to $c$ and rotation gives only a small, weakly dispersive correction. If $f=0$, the gravity waves are nondispersive, with speed $c$; at $k=0$ the phase velocity formula is undefined.

The wavefunction is a plane wave and a simultaneous formal eigenfunction of the momentum operator and the time-translation energy operator:Its momentum and energy are therefore definite, while its probability density is spatially uniform:For a particle of mass $m$, the probability current is $j=(\hslash /m)\mathrm{Im}(\overline{\varphi }{\partial }_x\varphi )=\hslash k{A}^2/m$, directed according to the sign of $k$.

On the whole real line, a nonzero constant-amplitude plane wave is not square-integrable and cannot be normalised to total probability one. It represents an idealised momentum eigenstate, interpreted through generalized eigenfunctions or as a limit of wave packets. Also, arbitrary $k,E$ do not automatically solve the Schrodinger equation for a specified potential energy: for a free particle, $E={\hslash }^2{k}^2/(2m)$, or this value plus the constant potential energy for a constant-potential region. If $A=0$, the zero function is not a physical quantum state.

The oscillator-generated covariant photon Fock space has an indefinite Hermitian form, not a positive Hilbert space inner product. For a normalizable one-photon wave packet of polarization $\lambda$, its squared norm is proportional to $-{\eta }^{\lambda \lambda }\int d^{3}p|f(\mathbf{p}){|}^2$. Thus $\lambda =0$ gives a negative-norm photon state. The divergent ${\delta }^3(0)$ of an unsmeared momentum eigenstate is a separate normalization issue, avoided by the wave packet.

Choose contravariant polarization vectors ${ϵ}^{0\mu }=(1,0)$, ${ϵ}^{3\mu }=(0,\hat{\mathbf{p}})$, and two ${ϵ}^{r\mu }=(0,{\mathbf{e}}_r)$ with ${\mathbf{e}}_r\cdot \mathbf{p}=0$ for $r=1,2$. The third spatial polarization is longitudinal polarization; the first two are transverse polarization. The timelike photon polarization is distinct from the longitudinal one.

The Gupta-Bleuler quantization condition sets the divergence of the positive-frequency part of a quantum field to zero on physical states:Restore the factors $e^{-i|\mathbf{p}|t}$ to the annihilation terms. Since $p_{\mu }{ϵ}^{0\mu }=|\mathbf{p}|$ and $p_{\mu }{ϵ}^{3\mu }=-|\mathbf{p}|$, Fourier transform gives the equivalent conditionIts sign depends on the chosen sign of the longitudinal polarization vector; the covariant condition does not.

To see its content for a general Fock state, temporarily discretize momentum and decompose one unphysical oscillator sector as $|\Psi ⟩={\sum }_{n,m\ge 0}|n,m{⟩}_{0,3}\otimes |{\chi }_{nm}{⟩}_{\perp }$, where $|n,m⟩=(a^{0†}{)}^n(a^{3†}{)}^m|0⟩/\sqrt{n!m!}$. The temporal oscillator obeys $a^0|n,m⟩=-\sqrt{n}|n-1,m⟩$, while $a^3|n,m⟩=\sqrt{m}|n,m-1⟩$. The condition therefore becomesEquivalently, the allowed finite-particle states use the transverse creation operators and only $b^{†}=a^{0†}-a^{3†}$ in the unphysical sector. Indeed the constraint acts on a polynomial of the two unphysical creation operators as $-{\partial }_{a^{0†}}-{\partial }_{a^{3†}}$, whose kernel consists of polynomials in their difference. Also $[b,b^{†}]=-1+1=0$, and a state containing $b^{†}$ is orthogonal to every constrained state because $b$ annihilates every such state. This argument applies mode by mode and extends by smearing to continuum momentum.

For one photon, $c_0{a}^{0†}|0⟩+c_3{a}^{3†}|0⟩$ is constrained only when $c_3=-c_0$, producing a null state. The Gupta-Bleuler null-state quotient removes these null directions. The condition excludes negative-norm physical states; quotienting its null states leaves the two positive-norm transverse photon polarizations. The condition alone gives a positive semidefinite Hermitian form, not yet a positive definite Hilbert space.

Assume a complete relativistically normalized basis of scalar momentum eigenstates transforms asThis states both reversal of spatial momentum and preservation of the normalization of basis quantum states. Antilinearity alone would not suffice: multiplying a conjugation operator by two gives an antilinear operator that multiplies squared norms by four.

Expand $|\varphi ⟩=\int d{\Pi }_p\varphi (\mathbf{p})|\mathbf{p}⟩$ and similarly for $|\psi ⟩$, with $d{\Pi }_p=d^{3}p/((2\pi {)}^{3}2E_p)$. This Lorentz-invariant phase-space measure is unchanged under $\mathbf{p}\mapsto -\mathbf{p}$. The antilinearity of the quantum time-reversal operator gives conjugated expansion coefficients. Orthogonality of the momentum eigenstates and cancellation of their unit-modulus phases yieldSince momentum reversal is a bijection of the complete basis, $\hat{T}$ is also onto. ThusThis proves antiunitarity, with the normalization hypothesis explicitly included. For scalar multiparticle quantum states the same argument uses the complete occupation-state basis and reverses all momenta; it is not restricted to a single-particle wave packet.

Work first in units $G=c=\hslash =k_B=1$. Black-hole thermodynamics starts with the laws of black-hole mechanics. For a connected stationary regular Killing horizon, the Zeroth law of black-hole mechanics makes its surface gravity $\kappa$ constant when the Einstein field equations and dominant energy condition hold. This parallels uniform equilibrium temperature. For neighboring stationary asymptotically flat four-dimensional Einstein-Maxwell black holes, the First law of black-hole mechanics iswhere ${\Omega }_H$ is the horizon angular velocity and ${\Phi }_H$ its electric potential relative to infinity. The last terms are rotational and electromagnetic work, analogous to the work terms in $\delta E=T\delta S+\text{work}$. The second law of black-hole mechanics, expressed by Hawking's area theorem, makes the total future horizon area nondecreasing for classical matter obeying the null energy condition and global assumptions such as strong asymptotic predictability. The third law of black-hole mechanics is an unattainability statement: subject to the usual regularity assumptions, the weak energy condition, and a bounded stress-energy tensor, no finite physical process reduces the surface gravity of a regular horizon to zero. It does not assert that extremal black holes have zero area or zero entropy.

Classical geometry alone fixes an analogy, not a nonzero physical temperature. Quantum field theory supplies Hawking temperature $T_H=\kappa /(2\pi )$. Comparing the area term of the first law with $T_H\delta S$ gives $\delta S=\delta A/4$ for nonextremal holes, and hence the standard Bekenstein-Hawking entropy. Restoring physical constants,Here restored $\kappa$ is the physical acceleration surface gravity and ${\ell }_P^2=G\hslash /c^3$ is the squared Planck length. The first-law comparison determines the area coefficient, leaving an additive entropy constant unspecified; the displayed formula is the usual convention. Its area scaling, rather than ordinary volume scaling, signals that the available thermodynamic degrees of freedom of a gravitating system are constrained by the horizon geometry.

Particle production explains the quantum input. For a real scalar field obeying the massless covariant wave equation, the Klein-Gordon inner product on solutions isThe wave equation makes its current divergence-free, so the product is independent of a Cauchy hypersurface when the boundary flux vanishes. Use suitably normalized wave packets or the appropriate continuum distributions. Asymptotic past and future Minkowski spacetimes give preferred positive-frequency solution mode bases $u_j^{\mathrm{in}}$ and $u_i^{\mathrm{out}}$, normalized by $(u_i,u_j)={\delta }_{ij}$, $(u_i^{\ast },u_j^{\ast })=-{\delta }_{ij}$ and $(u_i,u_j^{\ast })=0$. Between those regions, a nonstationary geometry generally has no preferred positive-frequency splitting.

The mode bases are related by a Bogoliubov transformation,The last identities follow from conservation of the Klein-Gordon inner product and are the canonical identities for a bosonic Bogoliubov transformation. Expanding the field in either complete basis and extracting its positive-frequency coefficient givesThus the state annihilated by every $a_j^{\mathrm{in}}$ hasThis particle number from Bogoliubov coefficients is nonzero when time evolution mixes positive and negative frequencies. It describes one state as vacuum in the early particle basis and populated in the late basis; it does not require an arbitrary choice of a vacuum at each intermediate time. For a finite implementable transformation the state remains a squeezed pure quantum state. In infinitely many modes, a common unitary bosonic Fock space implementation requires additional conditions, such as the beta map being a Hilbert-Schmidt operator; finite wave packet observables need not share every global divergence of an ideal continuum calculation.

For gravitational collapse, the future is not globally Minkowski: it contains a black hole. The relevant late out-basis includes modes reaching future null infinity together with modes entering the future event horizon. Modes on infinity alone are not a complete Cauchy basis. Nevertheless the same mode-mixing calculation determines the outgoing particle flux. Near a nonextremal horizon, the logarithm in the tortoise coordinate converts regular early null coordinates into the Hawking exponential ray mapwhere $U$ is the early affine null coordinate and $u$ is late retarded time at infinity. Backward propagation of a late outgoing mode $e^{-i\omega u}$ gives a factor $(U_H-U{)}^{i\omega /\kappa }$ for $U<U_H$. Its positive- and negative-frequency Fourier integrals differ by analytic continuation around the logarithmic branch. With a convergence regulator they reduce toTaking $ϵ\downarrow 0$ gives the thermal ratio of Hawking Bogoliubov coefficients, $|\beta {|}^2/|\alpha {|}^2=e^{-2\pi \omega /\kappa }$. Combining this with the canonical normalization yields the bosonic occupation $[e^{2\pi \omega /\kappa }-1{]}^{-1}$. Late-time wave packets turn formal continuum coefficients into a finite number flux.

A Schwarzschild black hole has $\kappa =1/(4M)$, so this spectrum has $T_H=1/(8\pi M)$. The exterior curvature potential partly reflects the outgoing modes. Its transmission probabilities are the greybody factors, giving, for one massless scalar species in the stationary late-time approximation,Thus the horizon temperature is universal, while the spectrum received at infinity is not a perfect featureless blackbody radiation spectrum. The derivation assumes the near-horizon quantum state inherited from a regular collapse vacuum and ignores rapid backreaction over the timescale of the packets. Exponentially blueshifted precursor frequencies indicate the usual short-distance assumption in the semiclassical calculation, not an independently demonstrated quantum-gravity description of the endpoint.

The outgoing positive energy flux drives black-hole evaporation. At leading order for a large isolated Schwarzschild black hole, area scales as $M^2$ and temperature as $M^{-1}$, giving $dM/du\sim -\gamma /M^2$ and a lifetime of order $M^3$, with the coefficient depending on species and transmission factors. Its negative heat capacity of a Schwarzschild black hole, $dM/d{T}_H=-8\pi M^2$, means that it heats up as it loses mass and cannot be in stable canonical equilibrium with an unlimited thermal reservoir. The approximation fails when quantum-gravitational scales are reached and does not fix whether the endpoint is complete evaporation, a remnant, or something else.

Area loss does not contradict the classical area theorem, because the quantum stress-energy tensor need not obey the classical null energy condition; negative horizon energy accompanies the positive outgoing flux. The appropriate thermodynamic statement is the generalized second law, involving $S_{\mathrm{BH}}+S_{\mathrm{outside}}$. Finally, Hawking radiation raises the black hole information paradox: if a pure quantum state completely evaporates into an exactly thermal mixed quantum state and its partners disappear, the result conflicts with ordinary unitary time evolution. Early outgoing radiation can be mixed simply because of its entanglement with interior modes while the complete state remains a pure quantum state; that fact alone is not information loss. Resolving the fate of all correlations through the evaporation endpoint requires more than the leading semiclassical flux calculation.

At the absolute wavenumber, the group velocity $d\omega /dk$ is zero. The corresponding accessible saddle therefore governs the disturbance seen at a fixed position. Its exponential growth rate isSaddle accessibility can be checked directly here. With $D=1+i{c}_d$, whose real part is positive, Fourier inversion gives the Green function of the linear complex Ginzburg-Landau equation,At fixed $x$ its exponential rate is $\mathrm{Im}{\omega }_0$, while along the packet centre $x=Ut$ it is $\mu$. For real wavenumbers of Fourier modes, the temporal growth rate is $\mathrm{Im}\omega (k)=\mu -k^2$, so the flow is temporally unstable exactly when $\mu >0$.

Consequently the classifications for a localized disturbance in this laboratory frame areFor convective hydrodynamic instability, a travelling wave packet amplifies but the response at each fixed position decays. For absolute hydrodynamic instability, that fixed-position response amplifies. The equality is the marginal absolute threshold: its exponential rate is zero and the impulse response has a $t^{-1/2}$ prefactor. For $\mu <0$ there is no temporal growth; $\mu =0$ is temporally marginal. When $U=0$ the convective window is empty. These conclusions use the infinite-line impulse problem; a general dispersion relation requires its own spatial-branch or saddle selection.

The normalized Bloch states $|k⟩=N^{-1/2}{\sum }_{n=1}^N{e}^{ikna}|n⟩$ satisfy periodicity when $k=2\pi j/(Na)$. Acting with the hopping terms shifts the coefficients by one site, givingThe first Brillouin zone is a fundamental cell of reciprocal space, for example $[-\pi /a,\pi /a)$; wavevectors differing by $2\pi /a$ label the same state. It contains exactly $N$ distinct allowed wavevectors and hence $N$ orbital states. The written Hamiltonian has no spin index; if electron spin degeneracy is included there are $2N$ states.

For a narrow wave packet, group velocity and band-curvature effective mass areFor the usual $J>0$, the mass is negative in the outer half of the zone, $\pi /(2a)<|k|\le \pi /a$, taking one endpoint representative. At $|k|=\pi /(2a)$ the inverse mass vanishes and the mass diverges. For arbitrary nonzero $J$, the sign criterion is $J\cos ka<0$; $J=0$ is a flat band.

For electron charge $-e$ with $e>0$ in a constant electric field $\mathcal{E}$, the single-band semiclassical equations are $\hslash \dot{k}=-e\mathcal{E}$, $\dot{x}=v(k)$. Thus $k(t)=k_0-e\mathcal{E}t/\hslash$, and for $\mathcal{E}\ne 0$,These are Bloch oscillations of period $2\pi \hslash /(|e\mathcal{E}|a)$. The zero-field limit is $x-x_0=v(k_0)t$. This assumes a narrow packet and neglects scattering and interband transitions.

For many electrons obeying the Pauli exclusion principle, a partially filled band has available nearby states and can carry current under a displaced occupation distribution: it describes a metal within band theory. A completely filled band has cancelling velocities over the zone, so its current is zero; if separated from empty higher bands by a nonzero gap, it describes an insulator at sufficiently low temperature and weak electric field. With spin degeneracy this isolated band is filled at two electrons per site. The single-electron Hamiltonian alone does not specify the filling or a gap to a next band, so those extra physical assumptions are required for classifying a material.

Without self-gravity and for $k>0$,The phase velocity and group velocity areso

With self-gravity restored,The product changes sign atFor $k>k_{\mathrm{crit}}$, crests and a localized wave packet travel in the same radial direction. For $k<k_{\mathrm{crit}}$, the group velocity and phase velocity have opposite signs: the envelope and its wave energy propagate opposite to the individual crests. At $k_{\mathrm{crit}}$ the group velocity vanishes and neighboring wave components cause the packet to spread.

Consider first a scalar constant-coefficient Cauchy problem for a partial differential equation on the whole line. A translation-invariant spatial discretization is a convolution operator, so the discrete Fourier transform turns it into multiplication by its Fourier symbol. This reduction is exact because Fourier modes are the simultaneous generalized eigenfunctions of every translation-invariant stencil.

For a fully discrete one-step methodwrite $H_h(\theta )$ for its amplification factor. The Parseval identity givesConsequently the exact necessary and sufficient condition for stability on every bounded time interval $0\le n\Delta t\le T$ iswith $C$ independent of the mesh. Sufficiency follows directly from Parseval; necessity follows by choosing transformed initial data concentrated where the multiplier is largest. For a contractive scheme one can take $C=0$, yielding the familiar Von Neumann stability analysis condition $|H_h(\theta )|\le 1$.

For a semidiscretization $u_t=A_{h}u$ with scalar symbol $a_h(\theta )$, the corresponding exact criterion isIndeed, the Fourier multiplier of the solution operator is $e^{t{a}_h(\theta )}$. For systems, eigenvalues alone cease to be sufficient when the matrix symbol is non-normal; the necessary and sufficient statement is the uniform bound $\Vert e^{t{A}_h(\theta )}\Vert \le C_T$.

On the whole lattice, a finite stencil defines a Laurent operator, whose operator norm is the essential supremum of its symbol. Restricting the same stencil to a half-line gives a Toeplitz operator. The Cauchy symbol condition remains necessary for a stable initial-boundary scheme because data localized far from the boundary behave like the whole-line problem for a finite time. It is not sufficient: the boundary closure can support growing modes or amplify incoming modes even when every interior Fourier mode is stable.

Boundary stability therefore requires a uniform estimate for the forced half-line recurrence. After a Laplace transform in time and a Fourier transform in tangential variables, one solves a normal-direction recurrence. The Uniform Kreiss--Lopatinskii condition requires the boundary equations to determine its decaying roots with a uniformly bounded inverse. In practical terms, one imposes one independent boundary condition for each incoming characteristic or numerical mode and none for outgoing modes. The group velocity $d\omega /dk$ identifies the direction in which a narrow wave packet carries energy, so its sign helps determine which boundary is inflow and explains why a numerically generated high-frequency branch can require a boundary condition different from that suggested by its phase velocity.

For a two-step method, the Fourier substitution produces an amplification polynomial of a multilevel finite difference scheme. Every root must lie in the closed unit disk, and unit-modulus roots must be simple, uniformly in the wavenumber. For example, leapfrog differencing of $u_t+c{u}_x=0$ givesand henceThe roots have unit modulus when $|\nu \sin \theta |\le 1$, giving the usual Courant restriction $|\nu |\le 1$. At the endpoint, a wavenumber with $|\nu \sin \theta |=1$ produces a repeated unit root and violates the uniform root condition; strict $|\nu |<1$ avoids this marginal linear growth. The second root is the familiar oscillatory computational mode, illustrating why a multilevel scheme requires its full amplification polynomial rather than a single multiplier.

A wave crest moves at the phase velocity, while the envelope of a narrow wave packet moves at the group velocity. Since $c_g>c_p$ for every nonzero wavenumber, the crests move more slowly than the packet.

For particle positions ${\mathbf{r}}_1,{\mathbf{r}}_2$, the two-particle Schrodinger equation isWrite $|a{⟩}_i$ for the localized wave packet ${\psi }_{ai}$ and $d_{ab}=|x_{a1}-x_{b2}|$. Neglecting packet spreading and branch overlap, the initial product state evolves branchwise asup to phases generated independently on the two particles. These branch-dependent phases generally cannot be separated into one phase depending only on $a$ and one depending only on $b$, so the Newtonian gravitational potential energy creates gravitationally induced entanglement.

If $d_{10}=|x_{11}-x_{02}|=d$ is much smaller than the other separations, remove their nearly common phase and retain onlyThe state is approximatelyIts concurrence is $|\sin (\varphi /2)|$, so it becomes maximally entangled first at $\varphi =\pi$. For $m_1=m_2=m$,Thus the near-maximal entanglement time is about $\boxed{10\mathrm{s}}$ within the stated approximation.

A single prescribed classical gravitational potential gives a Hamiltonian of the form $H_1[\Phi ]\otimes I+I\otimes H_2[\Phi ]$. Its evolution factorizes as $U_1\otimes U_2$ and preserves every initial product state, so it cannot generate this entanglement. A semiclassical mean field sourced only by expectation values likewise gives each particle a local one-body potential and does not provide a quantum mediator carrying branch correlations.

An entanglement witness is a Hermitian operator $W$ whose expectation is nonnegative on every separable state but negative on at least one entangled state. At $\varphi =\pi$, defineThe largest Schmidt coefficient of $|{\Psi }_{\ast }⟩$ is $1/\sqrt{2}$, so every product state $|u⟩|v⟩$ in the four-dimensional branch subspace satisfies $|⟨{\Psi }_{\ast }|u,v⟩{|}^2\le 1/2$. By closure under convex combinations, $\mathrm{Tr}(W{\rho }_{\mathrm{sep}})\ge 0$ for every separable mixture, whereasA negative measured value therefore certifies entanglement. Under the assumptions that the masses began unentangled and interacted only through gravity, such certification would show that the mediator can transmit quantum coherence; it would be evidence against a purely classical gravitational channel and for the quantum nature of gravity.

On a ray $x/t=V$, a localized wave packet is governed by a physically selected saddle point satisfying $d\omega /dk=V$ in its analytic dispersion relation. The exponent's real growth is $\mathrm{Im}({\omega }_{\star }-V{k}_{\star })$. The fixed-frame case $V=0$ gives the absolute growth rate. This ray formulation makes the frame dependence of convective hydrodynamic instability explicit.

Wave dispersion means that phase velocity depends on wavenumber or, equivalently, wavelength. Wave packets generally travel with group velocity, which need not equal phase velocity. The linear rotating shallow-water dispersion relation ${\omega }^2=f^2+c^2{k}^2$ has faster long-wave phase motion but slower long-wave packet motion when $f\ne 0$.