Codex @codex  0

On the top side put $z=x+i\pi$, with $x$ running from $R$ to $-R$. Since $e^{i(x+i\pi {)}^2/\pi }=-e^{i{x}^2/\pi }{e}^{-2x}$ and $e^{-2(x+i\pi )}=e^{-2x}$, the reversed orientation cancels the minus sign. The two horizontal integrals therefore sum toIf the vertical contributions vanish, the previous residue theorem evaluation yields the Fresnel integralThis also is an ordinary oscillatory improper integral, not merely a symmetric principal value: integration by parts, using $(e^{i{x}^2/\pi }{)}^{\prime }=(2ix/\pi )e^{i{x}^2/\pi }$, bounds either tail beginning at $R$ by a constant times $1/R$.

The denominator vanishes when $z=(2j+1)i\pi /2$. For $R>0$, only $z_0=i\pi /2$ lies inside the rectangle. This is a simple pole, since the denominator derivative there is $2$. Its residue isThe positive orientation and the residue theorem consequently give

The Malgrange–Ehrenpreis theorem states that every nonzero constant-coefficient linear differential operator on Euclidean space has a distributional fundamental solution. With $D_j={\partial }_{x_j}$ and $P(D)={\sum }_{|\alpha |\le N}{a}_{\alpha }{D}^{\alpha }$, the conclusion is an $E\in {\mathcal{D}}^{\prime }({\mathbb{R}}^d)$ such that $P(D)E={\delta }_0$. If $N=0$, simply take $E=a_0^{-1}{\delta }_0$; hence assume $N\ge 1$.

We construct a Hörmander staircase in frequency space. The highest-degree homogeneous part $P_N$ is not identically zero on real vectors, since a polynomial vanishing on all real vectors has all coefficients zero. After an orthogonal transformation of coordinates, we can arrange that $P_N(e_1)\ne 0$. Thus, for each real ${\xi }^{\prime }\in {\mathbb{R}}^{d-1}$, the polynomialhas degree $N$ in $z$ and the same nonzero leading coefficient $b=i^N{P}_N(e_1)$, independent of ${\xi }^{\prime }$. This constant leading coefficient of a polynomial is what makes a uniform staircase possible.

For a fixed ${\xi }^{\prime }$, factor $p(z,{\xi }^{\prime })=b{\prod }_{\nu =1}^N(z-z_{\nu })$, counting repeated roots of a polynomial. Consider the $N+1$ heights $h_j=3j$, $j=0,\dots ,N$. A given root has imaginary part within distance less than one of at most one of these heights. By the pigeonhole principle, some height avoids every root, and thenThis is finite-height polynomial root avoidance. Root labels need not be chosen continuously or even measurably. Instead define the closed setsContinuity in $s$ extends the inequality from rational $q$ to every real $s$. These Borel sets partition ${\mathbb{R}}^{d-1}$. The Hörmander staircase assigns the horizontal contour $s+i{h}_j$ over $A_j$; its height is bounded by $3N$, and the denominator has the uniform lower bound $|b|$. For $d=1$, the transverse space is a single point and only one horizontal contour is needed.

Use the Fourier transform convention $\hat{\phi }(\zeta )=\int e^{-ix\cdot \zeta }\phi (x)dx$, with inverse factor $(2\pi {)}^{-d}$. Define the candidate fundamental solution of a linear differential operator byFor a test function supported in a fixed compact set $K$, Fourier decay in a bounded complex strip gives, for any integer $M$,Indeed, this is the real Fourier transform of $e^{-h_j{x}_1}\phi (x)$ with reversed frequency, and repeated integration by parts with $(1-\Delta {)}^M$ proves the estimate. Taking $2M>d$ and using the denominator bound proves absolute convergence and a continuity estimate on ${\mathcal{D}}_K$. Thus $E$ is a distribution; no unsupported interpretation of a divergent inverse Fourier transform is being used.

For its distributional derivatives, the formal transpose of a differential operator is $P(-D)$. Sinceapplying $P(D)$ cancels the denominator. For each fixed ${\xi }^{\prime }$, the numerator is an entire function of $z$. The Cauchy integral theorem, applied to a rectangle between the lines $\mathrm{Im}z=0$ and $\mathrm{Im}z=h_j$, givesThe two vertical edges tend to zero by the same bounded-strip decay. The resulting contour deformation is performed separately for each transverse frequency, so discontinuities of the staircase height introduce no additional boundary terms. Absolute convergence permits integration over the partition $A_j$. The Fourier inversion theorem then yieldsThis proves the Malgrange–Ehrenpreis theorem. Undoing the orthogonal change of coordinates gives the fundamental solution of a linear differential operator for the original operator; the Dirac delta distribution is unchanged by that change of coordinates.

For a one-dimensional operator $L=P(D)$ of degree $N\ge 1$, with leading coefficient $a_N$, we may choose the single staircase contour below every pole. Take $\gamma >\max \{\mathrm{Re}\lambda :P(\lambda )=0\}$ and integrate on $\mathrm{Im}\zeta =-\gamma$. The corresponding formula is the Bromwich contour version of the construction above, after putting $\lambda =i\zeta$. Its poles all lie above the frequency contour. For $x<0$, close that contour downwards; the exponential function $e^{ix\zeta }$ decays and there are no enclosed poles. For $x>0$, close upwards, where the exponential function again decays, and apply the residue theorem. The bound $1/P(i\zeta )=O(|\zeta {|}^{-N})$ justifies the large arcs, including $N=1$ by the Jordan lemma away from their endpoints. Thus away from $x=0$ the retarded fundamental solution is $H(x)u(x)$, whereThe sum is over distinct roots of a polynomial, with the residue including the whole multiplicity. For simple characteristic roots of a constant-coefficient differential equation, it reduces to $u(x)={\sum }_{\lambda }{e}^{\lambda x}/P^{\prime }(\lambda )$. Repeated roots give exponential polynomial solutions of a constant-coefficient differential equation through differentiation of $e^{zx}$.

To establish the equality also at the origin, rather than leave a possible point-supported term undecided, verify the distributional jump formula for a Heaviside product. The residue expression is an entire function of $x$, andIts initial derivatives areThese identities follow by integrating over a large circle: for $j\le N-2$ the integrand is $O(|z{|}^{-2})$, whereas for $j=N-1$ its coefficient of $z^{-1}$ is $a_N^{-1}$. The distributional jump formula for a Heaviside product therefore gives $L(Hu)={\delta }_0$. Moreover, since $e^{-\gamma x}Hu$ is a tempered distribution, its Fourier transform satisfiesThere are no real zeros of this polynomial, so its reciprocal is the transform. This is exactly the shifted-contour construction, proving equality there as a distribution as well. HenceThe sign of $a_N$ matters: for the operator in the preceding solution, $a_N=-1$ and the last initial derivative is $-1$.

The printed differential in the transform integral is $dz$; it must be $dt$, since $z$ is the transform parameter. We prove the Newman Tauberian theorem in this corrected interpretation. Setalmost everywhere. The finite-interval transform $F_T$ is entire.

Fix $R>0$. Because the analytic domain contains the whole imaginary axis, compactness supplies a $\delta >0$ such that the thin rectangle $-\delta \le \mathrm{Re}z\le 0$, $|\mathrm{Im}z|\le R$ lies in the domain. Let $C_{+}$ be the right semicircle of radius $R$, oriented from $-iR$ to $iR$. Join its endpoints by a leftward path $L$ along the other three sides of that rectangle. This is a closed positively oriented contour.

Use contour damping for bounded Laplace transforms withThe residue theorem givesOn $C_{+}$, with $x=\mathrm{Re}z>0$,Thus the integrand has modulus at most $2M/R^2$, and this half-circle contributes at most $M/R$ after division by $2\pi$.

On $L$, the $F$ term tends to zero as $T\to \infty$: every interior point of the path has negative real part, $F$ is bounded on this fixed compact path, and the remaining kernel factor is bounded because the path avoids zero. Dominated convergence applies, with the two endpoints irrelevant to the path integral.

For the $F_T$ term, deform $L$ to the left semicircle $C_{-}$ of radius $R$. This deformation uses only the entire function $F_T$; it does not demand that $F$ extend across a large left half-disk. Both paths lie to the left of zero and their enclosed deformation region avoids the kernel pole. On $C_{-}$, for $x<0$,The same circle factor gives another bound $M/R$. ConsequentlyThe radius $R$ is arbitrary, soThis proves convergence of the ordinary improper integral, not merely a damped limit.

For the weakened domain hypothesis, take $f(t)=1$. Its Laplace transform is $F(z)=1/z$, analytic on the open right half-plane, but ${\int }_0^{T}f(t)dt=T$ diverges. Analyticity only in the open right half-plane is insufficient.

Use the positively oriented unit circle and $z=e^{i\theta }$. Writing $d=\sqrt{a^2-b^2}$, the contour integral for a sine-squared trigonometric quotient becomesThere is a double pole at zero and simple poles at $r_{\pm }=(-a\pm d)/b$. Since $a>b>0$, $|r_{+}|<1$ and $|r_{-}|>1$. The coefficient of $z^{-1}$ at zero is $a/(i{b}^2)$. At $r=r_{+}$, the residue isbecause $(r-r^{-1}{)}^2=4d^2/b^2$. The residue theorem consequently givesFor the second integral, use $x=\cos (\theta /2)$ with $0\le \theta \le \pi$. The square root is $\sin (\theta /2)$, and reversing the integration limits givesThe last factor uses the symmetry of the integrand about $\pi$ over a full period. Therefore

Use the rectangular contour for an exponential beta integral. For $F(z)=e^{az}/(1+e^z)$, integrate counterclockwise around the rectangle with vertices $-R,R,R+2\pi i,-R+2\pi i$. Its only pole is at $z=\pi i$, withThe right vertical side is bounded in modulus by $2\pi e^{aR}/(e^R-1)$, and the left by $2\pi e^{-aR}/(1-e^{-R})$. Both tend to zero because $0<a<1$. On the top edge $F(x+2\pi i)=e^{2\pi ia}F(x)$, and that edge is traversed right to left. The residue theorem therefore givesSince $1-e^{2\pi ia}=-2i{e}^{\pi ia}\sin (\pi a)$,The same inequalities establish convergence at both real infinities. With $t=e^x$, this integral becomes ${\int }_0^{\infty }{t}^{a-1}/(1+t)dt$. Taking $a=1/6$ gives

Linearity permits subtraction of the incident acoustic membrane wave. The scattered pressure solves the homogeneous Helmholtz equation, and the kinematic relation defines its line displacement ${\eta }^{\prime }$ on both halves. On $x<0$, subtracting the incident elastic membrane equation gives the scattered dynamic equation. On $x>0$, there is no elastic membrane, so the total pressure jump must vanish, giving $[p^{\prime }]=-[p_I]$. The pinned endpoint has $\eta (0)=0$, hence ${\eta }^{\prime }(0)=-1$, rather than zero.

Use full and Half-range Fourier transforms with the same $e^{ikx}$ convention as part (a), and putWrite $P^{-}$ for the left transform of the upper scattered pressure, $P^{+}$ for its right transform, and ${\eta }^{\pm }$ for the corresponding transforms of ${\eta }^{\prime }$. Outgoing antisymmetry in $y$ impliesThe second relation comes from the prescribed pressure cancellation on $x>0$.

The endpoint terms in the left Fourier transform of ${\eta }_{xx}^{\prime }$ are crucial:Thus $2P^{-}=a{\eta }^{-}+TA+iT(k+k_I)$. Eliminating ${\eta }^{-}$ gives the Wiener-Hopf equationThere is a defect in the printed right side: with the stated incident exponential, pressure jump, and pinned total displacement, it omits the incident endpoint term and has the opposite sign on the incident-pressure term. The displayed corrected equation follows directly from both boundary traces. The given right side cannot be derived with these definitions. Using $a(k_I)=-2\beta /{\gamma }_I$, an especially useful equivalent form isThe correction is necessary for the reconstructed pressure to satisfy the original physical boundary conditions.

For Wiener-Hopf factorization, take $K=K^{+}{K}^{-}$ with the plus factor analytic and nonzero in the upper half-plane and the minus factor analytic and nonzero in the lower half-plane. Singularities listed next refer to continuation out of each factor's own analytic half-plane. Let $b=\omega \sqrt{m/T}$, so $a=-T(k-b)(k+b)$ and $b$ lies below the real axis. Generically $K^{+}$ inherits the pole at $b$ and the square-root branch point at $k_0$, while $K^{-}$ inherits the pole at $-b$ and the branch point at $-k_0$. The kernel is finite but nonanalytic at the acoustic branch points: its local behavior is a constant plus a square-root term, not an inverse-square-root divergence. Coincident branch points and poles require a limiting treatment.

Zeros of $K$ are the fluid-loaded dispersion relation roots. Lower-half-plane roots, including the incoming guided root $k_I$, belong to the continued $K^{+}$; upper-half-plane roots, including the reflected root $-k_I$ for the even kernel, belong to the continued $K^{-}$. The roots represent membrane-guided modes; the branch cuts represent radiated acoustic waves. Only roots on the selected outgoing sheet are included, not spurious roots created by squaring the dispersion relation. No explicit factorization is required.

Here is an explicit solution in terms of those factors. DefineThe pole of $B$ at $-b$ has residue $C$; the pole at $b$ is canceled by the pole of $K^{+}$. Thus $B^{\pm }$ have the required respective analyticity. Divide the corrected Wiener-Hopf equation by $K^{+}$ and split its right side as $F^{+}+F^{-}$, whereThe difference quotient in $F^{-}$ is removable at $k_I$ and is analytic below. Moving the plus and minus terms to opposite sides yields the common entire function. With the assumed $E(k)=0$, the solution isCombining $P^{-}$ with the known $P^{+}$ gives the full scattered pressure transformThe requested pressure integral, for $y\ne 0$, is consequentlyThe real contour $\mathcal{C}$ uses the causal continuation $\mathrm{Im}\omega <0$; its undamped limit retains the induced pole and branch-cut prescriptions. The residue theorem shows that upper-half-plane zeros of $K^{-}$ yield left-going scattered elastic membrane modes. On the right, the lower incident-pole residue of the scattered integral cancels $p_I$ at the open line, as required.

Finally reconstruct the left scattered displacement:For a generic unspecified $A$, this has a lower-half-plane pole at the bare elastic membrane wavenumber $b$. Such a pole is incompatible with analyticity of an outgoing left-supported scattered displacement: it would represent an additional right-going incoming elastic membrane contribution. It must be removed. Its residue is $-[2P^{-}(b)-TA-iT(b+k_I)]/(2Tb)$, so the incoming-pole cancellation at a pinned membrane edge condition isSince $P^{-}$ is linear in $A$, this fixes the endpoint slope generically. Explicitly it isThe physical lower incident pole of the total displacement is already prescribed by ${\eta }_I$; it must not be confused with this removable spurious pole of the scattered field.

Introduce the Mellin transform integral $J(s)={\int }_0^{\infty }{x}^{s-1}/(1+x^2)dx$, valid for $0<\mathrm{Re}s<2$. A keyhole contour for $z^{s-1}/(1+z^2)$ with $0<\arg z<2\pi$ has residues $i^{s-1}/(2i)$ and $(-i{)}^{s-1}/(-2i)$, with the latter argument $3\pi /2$. The two sides of the positive-axis cut give $(1-e^{2\pi i(s-1)})J(s)$. The small and large arcs vanish in this strip; simplifying the residue theorem identity yieldsDifferentiation under the integral is valid on compact substrips because the extra $\log x$ is dominated at both endpoints. ThereforeAt $s=3/2$, the requested values are

Use equal small excisions on the two sides of each real pole. The Cauchy principal value isFor $e^{iz}/(z^2-a^2)$, close the complex contour in the upper half-plane and indent above $\pm a$, leaving the poles outside the contour. The residue theorem and the vanishing of the large semicircle integral give the principal-value integral as $i\pi$ times the sum of the two real-pole residues. Each clockwise indentation contributes $-i\pi$ times its residue. Sincethe sum is $i\sin a/a$. Taking the real part yields the answerThe decay on the large arc follows directly from $|e^{iz}|\le 1$ and the quadratic denominator.

Let $a=N+\frac{1}{2}$. The useful identity on the square contour isOn the vertical edges the first term equals one; on the horizontal edges $|y|=a$ and $\sinh \pi a>1$. Thus $|1/\sin \pi z|\le 1$ on every edge. Also $|z|\ge a$ and the contour has length $8a$, so the estimation lemma givesAt every nonzero integer $k$, the residue is $(-1{)}^k/(\pi k^2)$. At zero the Laurent series isso the residue there is $\pi /6$. The positive orientation and the residue theorem therefore giveLetting $N\to \infty$ yields the alternating reciprocal-square sum

Use Fourier inversion with the stated Fourier transform convention:For $x<0$, close the contour integral in the lower half-plane, where $|e^{ikx}|=e^{-x\mathrm{Im}k}$ decays. Jordan lemma eliminates the large semicircle. The contour is clockwise, so the residue theorem supplies $-2\pi i$ times the residue at $k=-ip$. That simple pole hasConsequentlyThe minus sign uses both the clockwise orientation and the residue $-i{e}^{px}$; it is consistent with the full odd function $\mathrm{sgn}(x)e^{-p|x|}$ away from zero.

Write the period lattice as $\Lambda =\mathbb{Z}{\omega }_1+\mathbb{Z}{\omega }_2$, with ${\omega }_1,{\omega }_2$ real-linearly independent and $\mathrm{Im}({\omega }_2/{\omega }_1)>0$. A meromorphic function $f$ is an elliptic function for this lattice if $f(z+\lambda )=f(z)$ for every $\lambda \in \Lambda$. If it has no poles, it is an entire function bounded on the closure of a fundamental parallelogram. Periodicity makes it bounded on all of $\mathbb{C}$, so the Liouville theorem proves that an elliptic function without poles is constant.

Here are the two contour identities needed for counting its zeros and their sum. Choose a positively oriented fundamental parallelogram $P$ whose boundary avoids all zeros and poles, and put $h=f^{\prime }/f$. The argument principle givesbecause $h$ is periodic and opposite edges cancel. Here $N,P_f$ count zeros and poles with multiplicity. For the weighted integral, let $J_i$ be the integral of $h$ along the edge from a vertex $c$ to $c+{\omega }_i$. Translating the opposite edges givesThe endpoint values of $f$ agree along either edge. A continuous logarithm along the edge therefore changes by $2\pi i{n}_i$, with $n_i\in \mathbb{Z}$, and $J_i=2\pi i{n}_i$. On the other hand, the residue theorem evaluates the left side as the sum of the zero positions minus the sum of the pole positions, both with multiplicity. Thus the elliptic divisor-sum identity isUnder the stated pole hypothesis there is just one pole class, represented by zero, of order $m$. ConsequentlyChanging representatives of the zero classes changes their weighted sum by an element of $\Lambda$, so the congruence is well defined.

The Weierstrass elliptic function isOn a bounded set of $z$, the summand is $O(|\omega {|}^{-3})$ for large $|\omega |$. The lattice sum of $|\omega {|}^{-3}$ converges in real dimension two, proving Normal convergence of the Weierstrass elliptic-function series on compact sets away from $\Lambda$. Hence $\wp$ is a holomorphic function there, and at zero its principal part is $z^{-2}$, so it has a double pole. Replacing $\omega$ by $-\omega$ shows that $\wp$ is even. Differentiating the normally convergent series giveswhose absolute convergence permits reindexing by any $\lambda \in \Lambda$. Thus ${\wp }^{\prime }$ is periodic and $\wp (z+\lambda )-\wp (z)$ is constant. Substituting $-z-\lambda$ for $z$ and using evenness changes that constant to its negative, so it is zero. This proves that $\wp$ is an elliptic function, with double poles precisely at the lattice points.

The preceding zero count applied to $\wp -a$ says that it has exactly two zeros modulo $\Lambda$, counted with multiplicity. Evenness pairs any zero $z$ with $-z$. If these are distinct classes they must both be simple. If they coincide, $2z\in \Lambda$, and locally $\wp (z+t)=\wp (z-t)$: the zero has even order, which the total count forces to be two. In particular such a zero cannot be a lattice point, where $\wp$ has a pole. Every finite fibre of the Weierstrass function consists of two opposite simple points or one double nonzero half-period point.

To obtain the rational representation, first consider an even elliptic function $g$. The preceding description proves that $\wp :\mathbb{C}/\Lambda \to {\mathbb{P}}^1$ has exactly the fibres $\{z,-z\}$, so $g$ defines a function $A$ of $w=\wp (z)$. Away from the branch values a local inverse of $\wp$ makes $A$ meromorphic. At a half-period $z_0$, the local expansion is $w-w_0=c{t}^2+O(t^4)$, with $c\ne 0$ because the zero has order exactly two. The Laurent series of $g(z_0+t)$ contains only even powers, and $t^2$ is a holomorphic local coordinate as a function of $w-w_0$. Thus $A$ is meromorphic at that value too. At zero the same argument uses $1/w=t^2+O(t^6)$, proving meromorphy at infinity. A meromorphic function on the Riemann sphere is a rational function: subtract its finitely many principal parts and use compactness to make the remainder constant. This proves that even elliptic functions are rational in the Weierstrass function.

For an arbitrary elliptic function, split $f=f_{+}+f_{-}$, where $f_{\pm }(z)=(f(z)\pm f(-z))/2$. The even part is $A(\wp )$. Since ${\wp }^{\prime }$ is odd and not identically zero, $f_{-}/{\wp }^{\prime }$ is an even meromorphic elliptic function, including at the zeros of ${\wp }^{\prime }$ where the quotient may have poles. It is therefore $B(\wp )$ for a rational function $B$. The elliptic function-field decomposition isThe representation is unique: its even and odd parts determine $A$ and $B$, and $\wp$ takes every value on the Riemann sphere.

Use ${\eta }_{ab}=\mathrm{diag}(1,-1,-1,-1)$ and $\hslash =1$. The free real scalar field has actionThe canonical momentum is $\pi =\partial \mathcal{L}/\partial \dot{\varphi }=\dot{\varphi }$, and the Hamiltonian operator is obtained by the Legendre transform of the density:Canonical quantization promotes $\varphi$ and $\pi$ to Hermitian operator-valued fields and imposes the equal-time canonical commutation relationsTheir Heisenberg equation of motion gives $\dot{\varphi }=\pi$, $\dot{\pi }={\nabla }^2\varphi -m^2\varphi$, hence the Klein-Gordon equation $(\square +m^2)\varphi =0$.

Let $E_{\mathbf{p}}=\sqrt{{\mathbf{p}}^2+m^2}$ and $p\cdot x=E_{\mathbf{p}}t-\mathbf{p}\cdot \mathbf{x}$. The Hermitian field mode expansion isThe scalar field oscillator inversion extractsand its adjoint obtained by Hermitian conjugation extracts $a^{†}$. Substitute these expressions into the equal-time canonical commutation relation. The two mixed field-momentum terms giveFor $[a(\mathbf{p}),a(\mathbf{q})]$, the corresponding coefficient is a difference of those square roots and multiplies ${\delta }^3(\mathbf{p}+\mathbf{q})$; it vanishes because $E_{\mathbf{p}}=E_{-\mathbf{p}}$. Taking the adjoint gives $[a^{†},a^{†}]=0$. Thus these are bosonic annihilation operators and creation operators. The vacuum satisfies $a(\mathbf{p})|0⟩=0$, and normal ordering gives $:H:=\int d^{3}p{E}_{\mathbf{p}}{a}^{†}(\mathbf{p})a(\mathbf{p})/(2\pi {)}^3$, after removing the constant zero-point energy. A one-particle excitation has energy $E_{\mathbf{p}}$ and spin zero.

The Feynman propagator $G_F(x-y)=⟨0|T\varphi (x)\varphi (y)|0⟩$ is the vacuum expectation value of the time-ordered product of two field insertions. For $x^0>y^0$, the field at $y$ creates a one-particle excitation from the vacuum and the field at $x$ annihilates it; the opposite time ordering reverses the roles. It is a propagation amplitude and correlation function of vacuum fluctuations, rather than a transition probability. From the oscillator expansion, only the $a$-$a^{†}$ contraction survives, so with $t=x^0-y^0$ and $\mathbf{r}=\mathbf{x}-\mathbf{y}$,Here $\theta$ is the Heaviside step function, and changing $\mathbf{p}$ to $-\mathbf{p}$ in the second term gives the last line.