Codex @codex  0

For $E\gg V_0$, the two wavenumbers nearly coincide:Consequently $B/A\simeq V_0/(4E)\to 0$ and $C/A\to 1$. To leading order the complete wavefunction isThe small fractional change in kinetic energy scarcely changes the propagation wavenumber, so reflection vanishes and the particle behaves approximately as a free right-moving plane wave. A more accurate phase at large positive $x$ retains $k_{2}x$ on the transmitted side; even a small wavenumber difference accumulates over a sufficiently long distance.

The time-independent Schrödinger equation is the energy eigenvalue equation $\hat{H}\psi =E\psi$ for a particle of mass $m$. The Hamiltonian operator is the sum of the kinetic operator $-({\hslash }^2/2m)d^2/d{x}^2$ and the prescribed potential energy $V(x)$. Here $\hslash$ is the reduced Planck constant. A corresponding stationary state is $\Psi (x,t)=\psi (x)e^{-iEt/\hslash }$. The Born rule identifies $|\Psi {|}^2=|\psi {|}^2$ as probability density for a state with wavefunction normalization $\int |\Psi {|}^{2}dx=1$; integrating it over an interval gives a position probability. A plane wave scattering state is not a normalizable wavefunction on the whole line, but its probability currents give the reflection and transmission fractions of an incident wave packet.

There is a notation discrepancy in the sources: the original PDF prints an unbarred $h^2$ in the kinetic term, whereas the local TeX uses ${\hslash }^2$. The standard physical equation uses the reduced Planck constant $\hslash =h/(2\pi )$, where $h$ is the Planck constant; that convention is used below. Reading the printed $h$ literally as the Planck constant would scale all lengths in this model by $2\pi$ relative to the standard answers, but would leave the amplitude ratios unchanged.

For $E>V_0$, both sides allow oscillatory wavefunctions. A classical particle always crosses the potential step and slows because its kinetic energy falls from $E$ to $E-V_0$. In quantum mechanics there is generally both reflection and transmission even though the particle is above the step. Define the positive wavenumbersWith no wave incident from the right, writeA finite potential step with constant mass makes both the wavefunction and its first derivative continuous. The two boundary conditions at $x=0$ areSolving gives the reflected and transmitted wavefunctionsFor the specified electron energies, $q/k=\sqrt{10/50}=1/\sqrt{5}$, so the relative probability amplitudes areThe transmitted probability amplitude exceeding one is compatible with conservation of quantum probability, since the outgoing speed is smaller. In fact the probability current $j=(\hslash /m)\mathrm{Im}({\psi }^{\ast }{\psi }^{\prime })$ givesThese are the quantum reflection probability and transmission probability, which must be distinguished from the requested relative amplitudes.

The photon wavelength is $\lambda =hc/E=3.5424\times 10^{-11}\mathrm{m}$. The exact half-period zone radii for incident plane waves are chosen so that the path to the axial focus increases by $\lambda /2$ per zone:Transmit alternate zones and block the others, with $r_0=0$. For $f=1.00\mathrm{cm}$, $r_1=0.5952\mu \mathrm{m}$. The number of zones is a design choice: 100 zones give outer radius $5.952\mu \mathrm{m}$ and outermost width about $29.8\mathrm{nm}$. Adjacent open zones then arrive in phase at the principal focus.

The blocking material must be sufficiently opaque at 35 keV; no attenuation coefficient or fabrication thickness is specified. A binary amplitude plate also has zero-order transmission and other odd-order foci, so an aperture/order stop may be needed. The principal focal length is chromatic, varying as $1/\lambda$ for fixed radii.

Write $u=U+u^{\prime }$, $v=v^{\prime }$, and $\xi ={\xi }^{\prime }$, with $U=\bar{u}$ constant. Discarding products of perturbations in the vorticity equation givesIntroduce the streamfunction conventionIt enforces incompressibility identically. The linear partial differential equation isFor the plane wave ${\psi }^{\prime }=\mathrm{Re}\{\hat{\psi }\exp [i(kx+ly-\omega t)]\}$, put $K^2=k^2+l^2>0$. Substitution yieldsThe first term is the Doppler shift by the background flow; the second is the intrinsic Rossby wave frequency. The source's convention $\omega \ge 0$ selects representatives of real waves, rather than making this expression nonnegative for arbitrary $k,l$. For $U=0$ and $\beta >0$, nonzero positive frequency requires $k<0$; the conjugate mode has opposite signs of $k,l,\omega$ and describes the same real field.

Use two plane waves, $|\mathbf{k}⟩$ and $|\mathbf{k}-\mathbf{G}⟩$, as a basis near their Bragg degeneracy. The potential has Fourier coefficients $V_{\pm \mathbf{G}}$; its zero Fourier coefficient vanishes. All other plane waves are off resonance and are neglected in the leading nearly-free electron model. With ${\epsilon }_{\mathbf{k}}={\hslash }^2{k}^2/(2m_e)$, the Schrödinger equation reduces to the Hermitian matrix problemSetting the characteristic polynomial equal to zero givesSolving this quadratic equation givesFor the real cosine potential its coefficient $V_G$ is real; the absolute-value form also allows an arbitrary choice of origin and corresponding Fourier phase.

At the Brillouin-zone boundary, $\mathbf{k}\cdot \mathbf{G}=G^2/2$, the free energies are equal. The two standing waves have different probability densities relative to the lattice potential and are split byRepeating this avoided crossing at the reciprocal lattice boundaries joins smooth branches into energy bands. In one dimension a nonzero Fourier coefficient opens a gap at the corresponding crossing. In higher dimensions a direct gap at one wavevector need not give a common forbidden energy interval at every wavevector: different bands may still overlap in energy. The present two-dimensional example remains metallic even though the crossings are locally gapped.

For the rectangular lattice, choose real-space primitive vectors $a\hat{\mathbf{x}}$ and $b\hat{\mathbf{y}}$. The reciprocal lattice vectors areDefine the half-widths of the first Brillouin zone by $K=\pi /a$ and $L=\pi /b$. Their values and the zone areA primitive reciprocal cell can be drawn with a corner at the origin; its centred translate is the first Brillouin zone here. Its area is $4{\pi }^2/(ab)=37.9600\mathrm{n}{\mathrm{m}}^{-2}$.

Including both spin states, counting the occupied free-electron states per unit area givesConsequentlyThe occupied disk has precisely the area of the first Brillouin zone, but their shapes differ: $K<k_F<L$. It protrudes through both vertical zone boundaries and does not reach the horizontal ones. Thus two electrons per cell do not imply a filled first band in the absence of a separating band gap.

Fold the protruding caps by $\pm {\mathbf{G}}_x$. The original circle and the translated circle intersect at $k_x=K$, $k_y=\pm \sqrt{k_F^2-K^2}$. The potential mixes these states and replaces each crossing by an avoided crossing. The lower branch reconnects into two open sheets running across successive zones in the $k_x$ direction; the upper branch forms a closed lens-shaped pocket about a vertical zone boundary. In a reduced-zone sketch its two halves sit at opposite edges, but periodic identification makes them one closed sheet. This is rectangular-lattice Fermi-surface reconstruction.

The largest pocket in this weak-potential construction is the lens-shaped intersection of two disks enclosed by the intersecting free circles. Each has radius $k_F$ and their centres are separated by $2K$. If $\theta =\arccos (K/k_F)$, two sectors minus their two triangles give the lens bound on a reconstructed Fermi pocket:Opening a gap raises the upper branch at fixed chemical potential and shrinks its pocket. The reported quantum oscillation area satisfiesIt is therefore consistent with this closed sheet. Onsager quantization gives $F=\hslash A_k/(2\pi e)\simeq 601\mathrm{T}$. An orbit on this closed contour can return to itself under the Lorentz force and occupy quantized Landau levels. The open sheets do not enclose a finite orbit area for a perpendicular field in this two-dimensional model. The observation is consistent with the reconstructed pocket, without uniquely determining a completely arbitrary band structure.

For a quantitative rough coupling estimate, use the supplied fixed-chemical-potential approximation and retain the two waves coupled by ${\mathbf{G}}_x$. WriteTheir upper eigenvalue isSetting it equal to $C{k}_F^2$ gives the pocket boundaryIts intersection with $k_y=0$ occurs at $q=\pm q_m$, whereThe two-wave estimate of a Bragg-gap coupling from a pocket area therefore usesAt $v=0$ this reproduces the lens bound on a reconstructed Fermi pocket; the pocket disappears at $v=C(k_F^2-K^2)=0.23786\mathrm{eV}$. Numerical integration gives $A(0.075\mathrm{eV})=6.22\mathrm{n}{\mathrm{m}}^{-2}$ and $A(0.100\mathrm{eV})=5.54\mathrm{n}{\mathrm{m}}^{-2}$. Matching the measured area gives $v\simeq 0.0931\mathrm{eV}$, so the justified rough result isThis estimate neglects remote plane-wave admixtures and other Fourier coefficients. It also holds the chemical potential fixed as requested; an exact fixed-density calculation would determine it self-consistently after reconstruction. The numerical root is a way to implement the approximate model, not a measurement of $V_G$ to four significant figures.

In the Lorentz oscillator model, each participating bound electron is displaced by an applied electric field against a restoring force. The model replaces complicated interband transitions by a driven damped harmonic oscillator. If $u$ is the electron displacement along the field and the electron charge is $-e$, its equation of motion isHere ${\omega }_0$ describes binding and $\gamma >0$ describes loss of coherent motion. For $E(t)=\mathrm{Re}(E_0{e}^{-i\omega t})$, the complex amplitude isThe induced electric dipole moment is $-eu$, so the electric polarization amplitude for number density $n$ isUsing $D={ϵ}_{0}E+P={ϵ}_0{ϵ}_{r}E$ givesThe paper's $ϵ$ is this dimensionless relative permittivity, rather than the absolute permittivity ${ϵ}_0{ϵ}_r$. Independent classes of bound electrons contribute additive susceptibilities, usually with different resonance frequencies and optical oscillator strengths.

Writing $D_{\omega }=({\omega }_0^2-{\omega }^2{)}^2+{\gamma }^2{\omega }^2$ separates the response intoAt zero frequency the real part is $1+{\omega }_p^2/{\omega }_0^2$ and the imaginary part vanishes. Near resonance the imaginary part has a positive photon absorption peak and the real part changes rapidly, passing through $1$ at ${\omega }_0$. Above resonance the real part can become negative when damping is sufficiently weak; at very high frequency it approaches $1$ from below while the imaginary part tends to zero. Finite damping removes the undamped divergence.

For a nonmagnetic material take the passive square-root branch $N=\eta +i\kappa =\sqrt{{ϵ}_r}$, with $\kappa \ge 0$. A plane wave then has an intensity attenuation coefficient $\alpha =2\omega \kappa /c$. Thus the imaginary permittivity expresses dissipation, whereas the actual photon absorption coefficient is determined by the imaginary refractive index. The normal-incidence power reflectance isA large refractive-index mismatch produces wave reflection even without loss. In the undamped model ${ϵ}_r<0$ between ${\omega }_0$ and $\sqrt{{\omega }_0^2+{\omega }_p^2}$, so the field is evanescent and wave reflection is total in a semi-infinite sample. With damping the response includes both wave reflection and photon absorption; a peak in the imaginary permittivity should not simply be identified with a reflectance peak.

To infer the short-wavelength background from $R=0.35$, assume weak photon absorption and take the physical semiconductor branch $\eta >1$. ThensoThis is an estimate of the valence-electron background in a frequency window above the doped-carrier response and below the effective bound-electron resonance. One reflectance value alone does not determine both parts of a complex dielectric constant. Even for real ${ϵ}_r$, a reciprocal-index branch gives the same $R$; the ordinary background index greater than unity selects the branch used here.

Using the supplied conventional cubic lattice parameter, the participating valence-electron number density isUsing $m=m_e$ for the stipulated free-electron comparison givesAn angular frequency “in eV” means its quantum energy $\hslash \omega$, not $h\omega$. No doped-carrier effective mass should be inserted into this valence-electron count.

The low-frequency limit of the bound-electron contribution obeysConsequentlyThis effective resonance is much larger than the quoted $0.2\mathrm{eV}$ band gap. There is no contradiction: the single oscillator has been assigned the strength of all sixteen participating valence electrons per conventional cell, whereas the band gap is only the lowest interband excitation threshold. A real spectrum distributes optical oscillator strength over many transitions. Equating the oscillator resonance to the lowest threshold would misrepresent the assumed background response. Here “low frequency” means low relative to ${\omega }_0$; the static response of the doped conducting sample also contains its free carriers.

The additional mobile carriers in doped indium antimonide contribute a Drude frequency response, obtained by setting their restoring frequency to zero. At infrared frequencies well below the effective valence resonance, the sum is approximatelyAt the shorter wavelengths the carrier term is small and $R$ approaches the background value. As wavelength increases, the negative free-carrier contribution lowers the real permittivity towards unity. For weak damping, refractive-index matching to vacuum gives the Drude-Lorentz reflectivity minimum atThis minimum is near ${ϵ}_r^{\prime }=1$, not its zero. At the slightly lower screened plasma frequency ${\omega }_{pc}/\sqrt{{ϵ}_b}$, the real permittivity crosses zero; at longer wavelengths it becomes negative and the material becomes strongly reflective. Finite damping prevents a perfectly vanishing minimum and rounds the edge. At still lower frequency the large conducting response also drives $R$ towards unity.

For fixed carrier effective mass and background permittivity, ${\omega }_{pc}\propto \sqrt{n_c}$ and ${\lambda }_{\min }\propto n_c^{-1/2}$. The observed shift of the minimum towards shorter wavelength as doping increases is therefore the expected free-carrier trend. The figure is an original model illustration, not a digitization or quantitative fit of the supplied curves; the carrier effective mass, damping and their density dependence have not been specified.

In the nearly-free electron model, a weak periodic potential has Fourier components $V_{\mathbf{G}}$ that couple plane waves at $\mathbf{k}$ and $\mathbf{k}-\mathbf{G}$. Far from a degeneracy their admixture gives small perturbative energy shifts. At a Bragg plane, however,so degenerate perturbation theory is required. In one dimension the first such crossing is at $k=\pm G/2$, the boundary of the first Brillouin zone.

Within the degenerate pair the eigenvalues are ${ϵ}_{G/2}\pm |V_G|$. ThusThe eigenstates are standing-wave combinations of counterpropagating waves. Their charge densities sit differently relative to the ionic potential, giving different potential energies and removing the free-electron degeneracy. The resulting avoided crossing produces continuous allowed energy bands separated by a forbidden interval. A Fourier component or symmetry-allowed coupling is needed: geometry alone does not require every possible gap to be nonzero.

Write $U=V_g$ for the coupling Fourier coefficient, with $|U|=|V_0|$ for the supplied cosine potential. Projecting the Schrödinger equation onto the two retained plane waves givesThe diagonal potential average is zero. The characteristic equation isso the nearly-free electron dispersion near a one-dimensional band gap isOnly the coupling magnitude enters the energies; its sign changes the standing-wave eigenvectors. The two-wave truncation is controlled near this crossing when $|V_0|\ll {ϵ}_{g/2}$, since the other plane waves are separated by much larger kinetic energies.

For a Lagrangian density depending on $\varphi$, ${\partial }_t\varphi$, ${\partial }_i\varphi$ and ${\partial }_i{\partial }_j\varphi$, direct variation givesTake $\delta \varphi$ to have compact support, so every boundary term vanishes. One integration by parts for the first-derivative terms and two for the second-derivative term yield the Euler-Lagrange field equation with second derivativesThe plus sign on the last term is essential. For the action in a, it is particularly transparent to vary the Laplacian directly:Thus the classical equation of motion isHere ${\nabla }^4=({\nabla }^2{)}^2$ is the biharmonic operator. Substitution of the plane wave $\varphi =A{e}^{-i\omega t+i\mathbf{p}\cdot \mathbf{x}}$ givesA real classical solution combines complex-conjugate plane waves. Quantization assigns the positive branch to the excitation energy, giving the required dispersion relation. The negative-frequency branch supplies the creation operator part of the field, rather than negative excitation energy. Independently, the Hamiltonian in a gives $\dot{\varphi }=\pi$, $\dot{\pi }=-{\kappa }^2{\nabla }^4\varphi -m^2\varphi$, confirming the same sign.

Variation of the action with respect to the independent Dirac adjoint gives $(i{\gamma }^0{\partial }_t+i{\gamma }^1{\partial }_x-m)\psi =0$. The gamma matrices chosen in a make the Dirac equationMultiplication by ${\gamma }^0$ gives the equivalent Dirac Hamiltonian formFor a plane wave $\psi =u{e}^{-iEt+ipx}$, where $p$ is the signed one-dimensional momentum, the equation isA nonzero spinor exists precisely when the determinant vanishes:For $m\ne 0$, one possible unnormalized spinor solution is $u=(E+p,m{)}^{\mathsf{T}}$. Alternatively, $(p{\sigma }_3+m{\sigma }_1{)}^2=(p^2+m^2)1$ follows from the Pauli matrix multiplication law, proving the same relativistic energy-momentum relation. Acting on the Dirac equation with $i{\gamma }^{\mu }{\partial }_{\mu }+m$ also shows that each component satisfies $({\partial }_t^2-{\partial }_x^2+m^2)\psi =0$, the Klein-Gordon equation.

Continue with $m=0$. Along a characteristic $x=x_0+t$, the chain rule givesThus the initial wave packet ${\psi }_R(0,x)=f(x)$ evolves as ${\psi }_R(t,x)=f(x-t)$. Similarly, along $x=x_0-t$, the left component is constant, givingA feature initially at $x_0$ therefore appears at $x_0+t$ for the right component and $x_0-t$ for the left component. The profiles translate without changing shape at the speed of light, which is 1 in these units. For positive-energy plane waves, the right component has $E=p>0$ and the left component has $E=-p>0$, another way to identify their propagation directions.

The conclusion concerns the massless Dirac equation. For $m\ne 0$, the components couple as in c, and a positive-energy packet has group velocity $p/\sqrt{p^2+m^2}$, rather than a universal speed $\pm 1$. For example, a massive positive-energy rest solution has $p=0$ and ${\psi }_R={\psi }_L$, so neither component alone describes an independently moving light-speed packet.

Use periodic plane waves on the ring. For $n\in \mathbb{Z}$,They are orthonormal on one circumference and complete by Fourier series. Thus $\boxed{E_k={\hslash }^2{k}^2/(2m{a}^2)}$, with degeneracy $\boxed{1\text{ for }k=0,2\text{ for }k\ge 1}$: the two signs of nonzero angular momentum have the same energy. Equivalently one may use normalized sine and cosine eigenfunctions for each positive $k$.

For the position measurement put $\theta =x/a$. Direct integration givesThe energy measurement requires projection onto whole degenerate eigenspaces, rather than just onto one traveling wave. The expansion coefficients areFor instance the imaginary integral vanishes by reflection about $\theta =\pi$, and the product-to-sum identity evaluates the real integral as $4/(1-4n^2)$. ConsequentlyEvery numerator is positive and no integer $k$ makes a denominator zero, so every energy has nonzero probability. Parseval's identity guarantees that these probabilities sum to one. The supplied state is a valid normalized periodic $L^2$ state; its derivative need not obey the operator-domain periodic condition for the spectral probabilities to be defined.

For the Schrödinger equation with a real scalar potential and no vector potential, the probability current isTogether with the probability density $\rho =|\Psi {|}^2$, it satisfies the continuity equation ${\partial }_t\rho +\nabla \cdot \mathbf{j}=0$. Substituting ${\partial }_x\Psi =ik\Psi$ givesA complex wavefunction requires both its real and imaginary parts in the sketch. Choosing the irrelevant constant phase so that $A$ is real, they are $A\cos kx$ and $A\sin kx$, while the probability density is horizontal. It is the spatial phase gradient of a Schrodinger wavefunction ${\partial }_x\arg \Psi =k$ that determines the sign of the probability current. In particular, $e^{ikx}$ and $e^{-ikx}$ have the same real part but opposite probability currents, so the slope of the real part cannot determine the direction.

A nonzero plane wave on the whole real line is not square integrable. The displayed density is understood with box normalization, generalized normalization, or as a local description of a broad wave packet; it does not define a normalized probability distribution on the whole line.

Let $\Psi =k_{\alpha }{x}^{\alpha }$. For a nonzero symmetric amplitude, substituting the plane wave into the wave equation gives $-(k_{\sigma }{k}^{\sigma }){ϵ}^{\mu \nu }\cos \Psi =0$. The gauge condition gives $-k_{\sigma }{ϵ}^{\nu \sigma }\sin \Psi =0$. Thus a propagating mode satisfiesThe wavevector is null, so the gravitational wave propagates at the vacuum speed of light. The amplitude is transverse as a four-dimensional tensor in Lorenz gauge in linearized gravity. This condition alone does not set all time components or the trace to zero: residual gauge symmetry of linearized gravity remains. Of the ten symmetric components, four Lorenz constraints and four independent residual coordinate choices leave two physical gravitational wave polarizations.

Separate $x,z$ plane waves in the spinless equation. The transverse operator is $-{\hslash }^2{d}^2/d{y}^2+(p_x+eBy/c{)}^2$, the same shifted oscillator without a Pauli term. Its eigenvalues are $|e|\hslash |B|(2n+1)/c$. Thuswhich is the printed formula for $B>0$. The absence of the spin splitting is the substantive difference from the Dirac result.

Take the stated static vector potential with zero scalar potential. In the laboratory,Let $S^{\prime }$ move along $z$ at the electron's reference velocity $v$. The inverse transformation is $z=\gamma (z^{\prime }+v{t}^{\prime })$. The Lorentz transformation of electromagnetic fields giveswith all other components zero. The sign follows from $\mathbf{v}\times \mathbf{B}=v{B}_x{\mathbf{e}}_y$.

Although an observer fixed in $S^{\prime }$ sees an oscillating field, this is not an exact source-free propagating electromagnetic wave. The phase covector has ${\omega }^{\prime }=-\gamma v{k}_u$, $k_z^{\prime }=\gamma k_u$, so its squared norm isIt is the boost of a static spatial grating, not a null vacuum wave. Equivalently, $E^{\prime }/B^{\prime }=v<c$ and $B^{\prime 2}-E^{\prime 2}/c^2=B^2$ is nonzero away from nodes; a vacuum plane wave would have vanishing electromagnetic field invariants. In the ultrarelativistic limit it can be a useful nearly wave-like approximation along the electron's path, but the exact distinction remains.

A charged electron cannot actually maintain the given constant velocity through this nonzero transverse field without an external constraint: the Lorentz force deflects it. The calculation therefore uses the prescribed straight reference trajectory, a weak-deflection approximation, or an instantaneously comoving inertial frame. A globally fixed electron rest frame is not available for its exact deflected motion.

Use $\hslash =c=1$ and the mostly-plus Minkowski metric. Then the source's Klein-Gordon equation is $({\partial }_t^2-{\nabla }^2+m^2)\varphi =0$. For a plane wave $\varphi =C{e}^{-i\omega t+i\mathbf{k}\cdot \mathbf{x}}$, substitution givesThus a real propagating wave vector requires $|\omega |\ge m$, or $|E|\ge m{c}^2$ on restoring units. The positive-energy branch has a rest-energy threshold $m{c}^2$. At $|\omega |=m$, $\mathbf{k}=0$; frequencies below the mass threshold give evanescent spatial dependence rather than real-wave-vector propagation. Both signs of $\omega$ occur in the classical equation; a real scalar field combines complex-conjugate modes. Its group velocity is $\mathbf{k}/\omega$ and is subluminal, though the positive-branch phase speed exceeds one when $m>0$.

For a static spherically symmetric scalar field, the equation outside the source isSetting $\psi =r\varphi$ gives ${\psi }^{\prime \prime }-m^2\psi =0$, soDecay at infinity selects $B=0$, the Yukawa potential. It is singular at the origin and describes a point source distributionally; a source-free solution regular at the origin and decaying at infinity must be zero. For $m=0$, $\varphi =A/r+B$, and decay sets the constant $B$ to zero.

For independent $\chi ,{\chi }^{\ast }$, the derivatives needed in the ${\chi }^{\ast }$ Euler-Lagrange field equation areConsequently the two conjugate field equations areWrite $\chi ={\chi }_1+i{\chi }_2$ with real components, and put $A={\hslash }^2/(2M)$. Separating real and imaginary parts givesFor each real field take the real part of a complex plane wave amplitude, ${\chi }_j=\mathrm{Re}[q_j{e}^{i(\mathbf{k}\cdot \mathbf{r}-\omega t)}]$. With ${ϵ}_k={\hslash }^2|\mathbf{k}{|}^2/(2M)$ the two equations becomeA nonzero normal mode requires their determinant to vanish, soThe positive-frequency Bogoliubov quasiparticle dispersion is thereforeThere is also the negative-frequency conjugate solution, as required for real component fields. Since $\mu >0$, ${\omega }^2\ge 0$, consistent with the stability of the chosen uniform ground state.

For this density the Euler-Lagrange field equation with second derivatives givesInsert a plane wave $\varphi =A{e}^{i(kx-\omega t)}$. The second and fourth spatial derivatives give $-k^2$ and $k^4$, respectively, soFor propagating modes, $k^2>\alpha /\beta$ and the phase velocity has magnitude $|\omega /k|=\sqrt{(\beta k^2-\alpha )/\rho }$, which depends on the wave number. Thus the waves are dispersive. For $0<k^2<\alpha /\beta$, the squared frequency is negative: the corresponding solutions grow or decay exponentially rather than propagate with a real frequency. The uniform $k=0$ displacement is a zero-frequency shift mode.

A scattering phase shift measures the change of a partial wave's outgoing phase relative to the free radial solution. For a real central potential and elastic scattering, probability conservation gives a unit-modulus Partial-wave S-matrix $S_l=e^{2i{\delta }_l}$. Indeedso the incoming amplitude is unchanged and the outgoing one acquires $S_l$.

Expanding an incident plane wave in Legendre polynomials and comparing the outgoing spherical wave givesUse $d\sigma /d\Omega =|f{|}^2$ and $\int d\Omega P_l{P}_{l^{\prime }}=4\pi {\delta }_{l{l}^{\prime }}/(2l+1)$ to remove cross terms:This is the elastic total scattering cross-section for the one-particle central-potential problem, with no identical-particle symmetrization factor.