Codex @codex  0

An acoustic phonon branch has $\omega (k)\to 0$ as $k\to 0$. Its long-wavelength group velocity is the speed of sound.

The experiment in A time of flight method to measure the speed of sound using a stereo sound card by Carlos C. Carvalho and others makes wave speed visible as a delay. Two powered microphones feed separate channels of one stereo line input. A clap beyond one end of their common axis reaches the nearer microphone first. The displayed signals show two separated onsets when enlarged; count their sample separation $N$ to obtain $\Delta t=N/f_s$ and infer speed of sound $c=d/\Delta t$ from the measured microphone spacing $d$. The reported sampling rate is $44100{\mathrm{s}}^{-1}$. Repeating at different spacings gives a distance-versus-delay plot whose slope estimates speed of sound; dependence on distance helps identify acoustic travel rather than a fixed electronic delay. These are the paper's apparatus and observed readout, not measurements made for this article.^[ref]

A controlled implementation also needs timing calibration for acoustic speed measurements. Put the microphones side by side to estimate the channel-delay offset ${\tau }_0$, swap channels to check its sign, and check the acquisition clock against an independent timing reference. Measure spacings between consistent points on the microphones, keep the source outside their span and on their axis, and fit $\Delta t=d/c+{\tau }_0$ by linear regression rather than forcing a zero intercept. Repeat pulses and compare the same onset criterion before reflections arrive. These are proposed calibration and control checks, not additional procedures claimed for the paper.

One sample corresponds to $1/f_s\simeq 22.7\mu \mathrm{s}$ for the reported sampling rate, so a delay of only a few samples has large relative measurement error. Background measurement noise, different microphone responses and threshold choice can move the apparent onset; reflections can supply a later competing pulse. Record temperature, keep air motion small, and compare several spacings. The result estimates speed of sound under those conditions; it cannot prove an exactly zero precursor below the instrument's sensitivity or establish finite propagation speed for every real medium. Both new diagrams in these sections show explicitly synthetic ideal wave equation solutions.^[ref]

The heat capacity ratio compares heat capacity at constant pressure $C_p$ with heat capacity at constant volume $C_V$ at the same equilibrium state and fixed composition:Here $C_p,C_V$ refer to the whole sample, in $\mathrm{J}{\mathrm{K}}^{-1}$; molar heat capacities $C_{p,\mathrm{m}}=C_p/n$ and $C_{V,\mathrm{m}}=C_V/n$ are per mole; specific heat capacities $c_p=C_p/m$ and $c_V=C_V/m$ are per unit mass. The amount of substance $n$ or mass $m$ cancels only when both capacities use the same normalization. Thus heat capacity ratio is dimensionless and independent of sample size. For a fixed-composition ideal gas with positive finite $C_V$, Mayer's relation giveswhere $R$ is the molar gas constant. More ways to store internal energy make a given temperature increase require more heat, bringing heat capacity ratio closer to unity. ^[ref]

The equipartition theorem predicts $C_{V,\mathrm{m}}=fR/2$ when $f$ independent quadratic terms in molecular energy are thermally active and their number stays fixed; consequently $\gamma =1+2/f$. A classical monatomic gas with no appreciable electronic excitation has $f=3$, giving $C_{V,\mathrm{m}}=3R/2$ and $\gamma =5/3$. For an ideal gas of diatomic molecules with three translational and two rotational terms active, but its vibrational normal mode effectively frozen, $f=5$ gives $C_{V,\mathrm{m}}=5R/2$ and $\gamma =7/5$. A fully classical harmonic vibration adds both kinetic energy and potential energy, giving $f=7$ and $\gamma =9/7$. These are model regimes, not universal constants of every gas: rotational freezing, electronic excitation, and changing composition invalidate the simple count.^[ref]

The mechanical significance of heat capacity ratio follows from the first law of thermodynamics. For a fixed amount of ideal gas in a reversible adiabatic process, $C_V(T)dT=-pdV$. Combining this with the ideal gas law gives the local relationThe partial derivative holds at fixed entropy and composition even when heat capacity ratio varies with temperature. The finite-path formulas $p{V}^{\gamma }=\mathrm{constant}$ and $T{V}^{\gamma -1}=\mathrm{constant}$ require the additional approximation of constant heat capacity ratio. A rapid compression also needs internal equilibration to approximate this path; rapid motion alone does not ensure reversibility. The same local response gives isentropic compressibility $1/(\gamma p)$ and, for small sound disturbances in this ideal gas regime, speed of sound $a=\sqrt{\gamma p/\rho }$, where $\rho$ is mass density. Thus thermal storage controls mechanical stiffness, sound propagation and the measurements in Rüchardt experiment and Clément-Desormes experiment.

Fluid-loaded plate. Use complex amplitudes with time dependence $e^{-i\omega t}$ and take the velocity potential in the fluid above the elastic plate to be $\Phi$. Let ${\rho }_0$ and $c_0$ be the undisturbed mass density and speed of sound. The linearized Euler equations give $p=i{\rho }_0\omega \Phi$; the kinematic boundary condition is ${\Phi }_y(x,0)=-i\omega \eta (x)$. We model loading from one fluid half-space, with no fluid loading on the other face. Loading by identical fluids on both faces instead doubles the fluid term below.

Define the Fourier transform and its inverse byThe Sommerfeld radiation condition selects $\hat{\Phi }(k,y)=\hat{\Phi }(k,0)e^{-\gamma y}$, with ${\gamma }^2=k^2-k_0^2$ and $k_0=\omega /c_0$. Here and in the diffraction calculation below, the outgoing square root is specified by first taking $\mathrm{Im}\omega >0$. Draw a vertical cut upwards from $+k_0$ and another downwards from $-k_0$, and putThus ${\gamma }_u$ is analytic in the upper half-plane and ${\gamma }_l$ in the lower half-plane. In the real-frequency limit,The Fourier transform contour $C$ follows the real axis with the limiting-absorption pole prescriptions. In particular, an outgoing propagating component has positive vertical wavenumber.

The kinematic boundary condition gives $\hat{\Phi }(0)=i\omega \hat{\eta }/\gamma$, so $\hat{p}(0)=-{\rho }_0{\omega }^2\hat{\eta }/\gamma$. Eliminating the fluid pressure from the bending balance gives the line-forced fluid-loaded bending plate kernelFor harmonic line force amplitude $F_0$ at the origin, the plate displacement and fluid velocity potential are thereforeAn arbitrary localized loading replaces $F_0$ by its Fourier transform. This representation both calculates the response and identifies its different wave components through poles and branch points.

For an observer $x=r\cos \theta$, $y=r\sin \theta$, with $0<\theta <\pi$, the method of steepest descent uses the saddle point $k_s=k_0\cos \theta$. At fixed angle away from grazing, the acoustic far field isThe $r^{-1/2}$ spreading is that of cylindrical acoustic waves. The phase velocity parallel to the plate is $\omega /|k|$: components with $|k|<k_0$ are supersonic in that sense and radiate, whereas those with $|k|>k_0$ are subsonic and exponentially decay normally to the plate. The latter have reactive fluid loading and zero mean normal acoustic energy flux.

There is also a real guided root $\kappa >k_0$. Indeed, $K(k_0+)=-\infty$, $K(+\infty )=+\infty$, andThe residue theorem then gives an outgoing displacement contribution $i{F}_0{e}^{i\kappa |x|}/K^{\prime }(\kappa )$. It transports energy along the plate while its fluid field decays as $e^{-\gamma (\kappa )y}$. Other zeros on continued sheets can describe leaky or resonant waves; branch point contributions describe near and grazing fields. A fixed-angle acoustic far field expansion is not uniform near the plate, where these components can matter.

For a propagating spectral component, write $\ell =(k_0^2-k^2{)}^{1/2}$. Its upward normal velocity is $i\ell \hat{\Phi }$ and its pressure is $i{\rho }_0\omega \hat{\Phi }$. Applying Parseval's identity to the time-averaged acoustic energy flux gives the upward power per unit length of the forcing line:Equivalently, the integrand is $|\hat{\eta }{|}^2/\ell$ with prefactor ${\rho }_0{\omega }^3/(4\pi )$. Integrating the radial acoustic intensity of the far field over a semicircle gives the same expression. The input power $\omega \mathrm{Im}(F_0^{\ast }\eta (0))/2$ also feeds guided waves along the plate, so it must not be identified with $P_y$ alone.

Rigid half-plane diffraction. Suppress the time factor and assume $0<{\theta }_0<\pi$. Write $a=k_0\cos {\theta }_0$, $b=k_0\sin {\theta }_0$, and denote the incident spatial field by $I=e^{-iax-iby}$. The Helmholtz equation, Sommerfeld radiation condition, and square-root branches are as above. Across the open positive axis the scattered field and its normal derivative are continuous. Across the rigid negative axis its normal derivative is also continuous, while its potential can jump. Let $J(x)$ be this jump, supported on $x<0$. The transformed scattered field isThe jump transform is analytic in the upper half-plane. The absence of an even scattered component follows from derivative continuity and the outgoing condition. Its normal derivative on either face is $-\gamma {\hat{J}}_u/2$.

On the negative axis this derivative must cancel $-ib{e}^{-iax}$. Since the Half-range Fourier transform of $H(-x)e^{-iax}$ is $i/(k+a)$, the Wiener-Hopf equation has the formwhere $G_l$ is the transform of an unknown function supported on $x>0$ and analytic in the lower half-plane. The incident pole $k=-a$ lies below the Fourier contour. Divide by ${\gamma }_l$ and perform pole subtraction in a Wiener-Hopf equation:The first term is upper analytic and the second lower analytic, since the latter's apparent pole cancels. Equality of the two continued sides gives an entire function. The finite-energy edge condition, with $J(x)=O(\sqrt{-x})$, and the decay of the transforms at infinity force that function to vanish by Liouville's theorem. ConsequentlyThis is the Wiener-Hopf solution of rigid half-plane diffraction with the present transform and time conventions.

Deforming to a steepest descent contour gives a residue from the incident pole precisely when $x+|y|\cot {\theta }_0<0$. At that pole $\gamma (-a)=-ib$; its residue producesThus the total geometrical optics field is $I$ plus a reflected wave in the upper reflection region; in the lower shadow region the scattered residue cancels $I$, and outside that shadow the direct incident wave remains. This gives both the amplitude and the spatial support of the geometrical optics fields rather than merely naming a reflected wave.

For $x=r\cos \theta$, $y=r\sin \theta$, $-\pi <\theta <\pi$, the saddle contribution uses ${\gamma }_l(-a)=-i\sqrt{2k_0}\cos ({\theta }_0/2)$ and ${\gamma }_u(k_s)=\sqrt{2k_0}\cos (\theta /2)$. The diffracted wave isThe total field is $I+{\varphi }_{\mathrm{GO}}+{\varphi }_d$ to this order. This far-field expression applies away from the shadow/reflection boundaries $|\theta |=\pi -{\theta }_0$. There the pole and saddle coalesce, and a uniform Fresnel integral approximation replaces the divergent separate terms by a smooth transition. The exact contour integral remains well defined.

An acoustic analogy is an exact rearrangement of the fluid equations into a chosen linear propagation operator acting on an acoustic variable, with everything left over placed on the right as effective forcing. The rearrangement becomes a sound-prediction method only after a reference medium, boundary conditions and approximations to the forcing are specified. In particular, a right-hand-side term need not represent independently generated sound.

Use Einstein summation convention and write ${\partial }_i=\partial /\partial x_i$. Differentiating the continuity equation in time and taking the divergence of the momentum equation eliminates ${\partial }_t(\rho u_i)$:Set $q=p-{\hat{p}}_0$ and $W_{ij}=\rho u_i{u}_j-{\sigma }_{ij}$. Since $p+\chi =q+P$, the constant $P$ has zero Laplacian. The prescribed reference mass density and reference speed of sound are independent of time, soCombining the two identities givesOnly conservation of mass and conservation of momentum have been used; no equation of state or energy equation was needed. Spatial variation of ${\hat{c}}_0$ creates no omitted derivative in this identity, because it multiplies a time derivative. The double divergence of the momentum flux tensor has the structure of an acoustic quadrupole.

For a localized flow, choose the reference fields to match the stationary surrounding medium: ${\hat{\rho }}_0={\rho }_0$ and ${\hat{c}}_0^2=c_0^2$ outside the flow, with $P$ chosen so ${\hat{p}}_0=p_0$ there. A uniform surrounding fluid permits ambient constant mass density and adiabatic sound speed, giving the familiar homogeneous wave equation. A nonuniform surrounding fluid calls for its actual stationary reference profiles, extended sensibly through the flow region. This makes the acoustic variable vanish in the unperturbed exterior and minimizes artificial contrast terms. In a uniform isentropic exterior the leading acoustic relation $p^{\prime }=c_0^2{\rho }^{\prime }$ also makes $Q$ vanish to first order. In a stratified exterior, propagation and entropy-advection effects can remain in $Q$, as the next part demonstrates.

Suppress tildes in the normalized equation and let $X=x-vt$. A Gross–Pitaevskii solitary wave stationary in these translating coordinates obeysThe boundary fixes the bulk complex argument and describes a localized disturbance with zero net vortex winding number; a vortex pair is possible. No extra plane-wave factor is needed when only the coordinates are changed. A full Galilean transformation of the field would additionally change the bulk flow and its complex argument convention.

Subtract the uniform background from the grand potential. A convenient dimensionless energy and the renormalized momentum of a condensate areThe subtraction in $\mathbf{p}$ removes the troublesome background complex argument gradient. In particular, for the usual subsonic localized far field $\psi -1=O(r^{-1})$, $\nabla \psi =O(r^{-2})$, and $1-|\psi {|}^2=O(r^{-2})$, the energy density is $O(r^{-4})$ and the momentum integrand is $O(r^{-3})$, both integrable in two dimensions. Merely writing $\psi \to 1$ without sufficient decay is not by itself a convergence proof. The common solitary-wave branch is subsonic: linearizing this normalization gives ${\omega }^2=k^2/2+k^4/4$ and speed of sound $1/\sqrt{2}$.

Integration by parts, with decaying variations and their boundary terms vanishing, givesThus the translating equation is precisely $\delta (E-v{p}_x)=0$. Along a differentiable family ${\psi }_{\lambda }$ with fixed bulk normalization, evaluate this stationarity on ${\partial }_{\lambda }\psi$ to obtainwhere $d{p}_x/d\lambda \ne 0$. This energy–momentum slope of a solitary wave does not require differentiating the Lagrange multiplier $v$ inside the field variation. At a turning point the parametrized identity remains true, while $E$ need not be a single-valued function of $p_x$ locally. The chosen factors make $i{\psi }_t=\delta E/\delta {\psi }^{\ast }$, consistent with the normalized equation; changing the normalization of $E$ alone would change the claimed slope.

Label the two atoms in primitive cell $j$ by $u_{2j}$ and $u_{2j+1}$, with $j$ understood modulo $N$ because of the periodic boundary condition. Choose the spring between them to have constant $\lambda$. Newton's second law gives

For a longitudinal normal mode, putPeriodicity requires $e^{2iNka}=1$, so $k=\pi r/(Na)$. Since the primitive-cell length is $2a$, wavevectors differing by $\pi /a$ describe the same mode, and the first Brillouin zone may be chosen asThe amplitudes obeySetting its determinant to zero gives the alternating-spring chain dispersion relationThe lower sign is the acoustic phonon branch and the upper sign is the optical phonon branch.

At the edge of the Brillouin zone,so the frequency band gap isIt vanishes at $\alpha =1$, when the apparent two-atom primitive cell can be reduced to one atom.

Near the centre of the Brillouin zone, a Taylor expansion givesandThe group velocity of the acoustic branch at long wavelength is the speed of sound,

Timing calibration for acoustic speed measurements separates travel delay from the offset ${\tau }_0$ between acquisition channels and checks the clock defining the sampling rate. A constant offset changes the intercept of a delay-versus-distance linear regression; a clock-scale error changes the inferred speed of sound itself.

For the homogeneous one-dimensional wave equation $u_{tt}=c^2{u}_{xx}$ with constant $c>0$, wave speed means the speed at which a travelling wave profile translates. In $F(x-ct)$, a feature initially at $x=s$ follows $x=s+ct$, giving signed velocity $+c$; in $G(x+ct)$ it follows $x=s-ct$, giving signed velocity $-c$. The characteristic lines for the one-dimensional wave equation therefore have slopes $dx/dt=\pm c$. The speed of either profile is $c$, even though the displacement at a fixed point can change sign or pass through zero.^[ref]

The same $c$ bounds the causal cone for the constant-speed wave equation. For a wave equation on a string, $c=\sqrt{T/\mu }$ depends on tension and linear mass density; for small acoustic waves in stationary uniform air it is the speed of sound. This interpretation assumes the stated ideal wave equation. In media with wave dispersion, phase velocity, group velocity and the earliest signal arrival have different meanings: a measured wave crest speed need not be the speed of a wave envelope, and neither alone establishes a causal cone for the constant-speed wave equation.