Alphabeta Math
LemmaStatement: AI-adaptedProof: AI-adaptedPipeline-generatedprecheck passjudge pass (gpt-6.1-sol)
How statement and proof provenance work

The first chip identifies the source of the statement or construction; the second identifies the source of its local proof or verification.

  • Literature-sourced: the exact statement appears in a cited source; only wording and notation differ.
  • AI-adapted: a semantically identical restatement of literature-sourced material, modulo indexing, notation, and boundary cases adopted by the library.
  • AI-generated: a genuinely novel statement formulated by AI, with no source for the claim itself.

These labels describe origin, not correctness: citations and verification chips remain separate evidence.

Spherical Hankel functions: explicit solutions, constant Wronskian and exact Dirichlet-to-Neumann bounds

Statement

For l∈N0 and z>0 set

hl(z)=(−i)l+1eizz∑j=0lij (l+j)!j! (l−j)! (2z)j.

Then hl is a complex-valued C2 function of z>0 solving

hl′′(z)+2zhl′(z)+[1−l(l+1)z2]hl(z)=0,

it never vanishes on (0,∞), and its Wronskian-type identity Im⁡(hl(z)‾ hl′(z))=z−2 holds. For k>0 and R>0 the spherical Dirichlet-to-Neumann coefficient

tl=k hl′(kR)hl(kR)

satisfies the exact bounds

Re⁡tl≤−1R,∣Re⁡tl∣≤l+1R,0<Im⁡tl≤k.

These are exact identities and inequalities for the displayed explicit functions, not asymptotic formulas; the argument uses only elementary differentiation and a constant Wronskian, so no large-l asymptotic input is required.

Facts & Assumptions

Given: An integer l∈N0, the explicit function hl of the statement on (0,∞), and reals k>0, R>0. Complex-valued functions of the real variable z are handled componentwise under the identification C≅R2 of C=R[x]/(x2+1) as the Euclidean plane and as a normed real algebra: what the identification preserves, with derivative and integral as in The derivative and the Riemann integral of a vector-valued function: an intrinsic derivative and a componentwise integral.

[L1]

Differentiation. The complex exponential is entire and its complex derivative is itself: the complex exponential is entire with exp⁡′=exp⁡. Linearity, product, reciprocal, and quotient rules for complex derivatives: sums, scalar multiples, products and (where the denominator is nonzero) quotients of complex differentiable functions are complex differentiable with the usual rules; in particular z↦eizP(z) and z↦eizP(z)/z are holomorphic on their domains for polynomial P, hence C2 on (0,∞) along the real axis. Sums, scalar multiples, products and quotients: (f+g)′(c)=f′(c)+g′(c), (αf)′(c)=αf′(c), (fg)′(c)=f′(c)g(c)+f(c)g′(c), and (f/g)′(c)=(f′(c)g(c)−f(c)g′(c))/g(c)2 when g(c)≠0: the same rules hold for real-valued functions, and (componentwise) for the real and imaginary parts of a complex-valued differentiable function: if y=u+iv with u,v real differentiable on an interval, then y′=u′+iv′, Re⁡y′=(Re⁡y)′, Im⁡y′=(Im⁡y)′.

[L2]

Constants from a vanishing derivative. A function continuous on an interval I whose derivative vanishes at every interior point of I is constant on I; consequently two such functions with the same derivative differ by a constant: a real-valued function continuous on an order-convex set, with derivative zero at every interior point, is constant; and two such functions with equal derivatives differ by a constant.

Proof

1.1L1algebra

Write P(z)=∑j=0lcjz−j with cj=ij(l+j)!j! (l−j)! 2j, so that hl=(−i)l+1eizP(z)/z, and put y=zhl=(−i)l+1eizP(z). Substituting y=eizP into the equation y′′+(1−l(l+1)z−2)y=0, which is equivalent to the stated radial equation because h=y/z gives h′′+2zh′+(1−l(l+1)z−2)h=1z(y′′+(1−l(l+1)z−2)y), produces by the product rule [L1] the equation P′′+2iP′−l(l+1)z−2P=0 for P. Inserting the finite expansion P=∑jcjz−j and matching the coefficient of z−m−2 gives [m(m+1)−l(l+1)]cm−2i(m+1)cm+1=0, that is cm+1=i(l−m)(l+m+1)2(m+1)cm, and the displayed factorials satisfy this recursion because cm+1cm=i(l−m)(l+m+1)2(m+1); the recursion terminates with cl+1=0, so P is exactly this polynomial and y, hence hl, satisfies the stated equation and is C2 on (0,∞). The leading term is (−i)l+1eiz/z with remainder O(z−2).

2.1step 1.1L1L2algebra

With y=u+iv, u=Re⁡y, v=Im⁡y, put W=uv′−vu′=Im⁡(y‾ y′). The componentwise real and imaginary parts of the equation of step 1.1 give u′′+(1−l(l+1)z−2)u=0 and v′′+(1−l(l+1)z−2)v=0; the product rule [L1] gives W′=uv′′−vu′′=u[−(1−l(l+1)z−2)v]−v[−(1−l(l+1)z−2)u]=0 on (0,∞), so W is constant by [L2]. Since P(z)→c0=1 and P′(z)→0 as z→∞, y=Cleiz(1+o(1)) and y′=iCleiz(1+o(1)), where Cl=(−i)l+1 has modulus one, so y‾ y′=∣y∣2(i+o(1)) and W→Im⁡i=1; hence W≡1. Consequently Im⁡(hl‾ hl′)=Im⁡((y‾/z)(y′/z−y/z2))=W/z2−∣y∣2Im⁡(1/z3)=z−2, which is nonzero for every z>0; therefore hl(z)≠0 on (0,∞).

3.1step 1.1step 2.1L1

Put F=∣y∣2 and c(z)=1−l(l+1)z−2. By step 2.1 the complex quantity A=y‾ y′ has Im⁡A=1, so A′=Re⁡(A)′=12F′′ is real (componentwise product rule [L1]); on the other hand A′=∣y′‾∣2+y‾ y′′=∣y′∣2−cF. Differentiating 12F′′=∣y′∣2−cF gives 12F′′′=2Re⁡(y′‾y′′)−c′F−cF′; since y′′=−cy, Re⁡(y′‾y′′)=−cRe⁡(y′‾y)=−cF′2, and with c′=2l(l+1)z−3 this yields the linear third-order equation F′′′+4(1−l(l+1)z−2)F′+4l(l+1)z−3F=0. Since P=∑jcjz−j with cj=ij(real), all cross terms of ∣P∣2=∑j,kcjck‾z−j−k with j−k odd cancel in the sum over (j,k) and (k,j), so F=∣P∣2 is a finite even polynomial in z−1, F=∑j≥0ajz−2j with real aj and a0=1; matching the coefficient of z−2m−3 in the third-order equation gives −4m(m+1)(2m+1)am+8ml(l+1)am+4l(l+1)am−8(m+1)am+1=0, i.e. am+1=(2m+1)[l(l+1)−m(m+1)]2(m+1)am, whose positive solution is aj=(2j)! (l+j)!22j(j!)2(l−j)! for 0≤j≤l and aj=0 for j>l, because the factor [l(l+1)−j(j+1)] vanishes at j=l. In particular aj>0 for j≤l, a0=1, and ∣hl(z)∣2=F(z)/z2=∑j=0lajz−2j−2≥z−2.

4.1step 2.1step 3.1algebra

On (0,∞) one has hl′hl=hl‾ hl′∣hl∣2 and, from step 3.1 and y‾ y′=F′2+i (step 2.1), Re⁡(hl‾ hl′)=Re⁡(y‾ y′z2−∣y∣2z3)=F′2z2−Fz3=−1z3∑j=0l(1+j)ajz−2j and Im⁡(hl‾ hl′)=z−2; dividing by ∣hl∣2=F(z)/z2 gives Re⁡hl′(z)hl(z)=−1z ∑j=0l(1+j)ajz−2j∑j=0lajz−2j and Im⁡hl′(z)hl(z)=1F(z). With z=kR and weights aj(kR)−2j≥0, the real part of tl=k hl′(kR)/hl(kR) is −1R times a weighted average of the numbers 1+j, j=0,…,l, so −(l+1)/R≤Re⁡tl≤−1/R; the imaginary part is k/F(kR), and since F(kR)=∑j≥0aj(kR)−2j≥a0=1 and F<∞, one has 0<Im⁡tl≤k. This proves all the stated bounds, and the bounds are exact identities for the displayed hl.

5.1step 1.1step 2.1step 3.1step 4.1∎

Collecting: step 1.1 gives the C2 solution property and the explicit leading term; step 2.1 gives the constant Wronskian Im⁡(hl‾hl′)=z−2 and nonvanishing; steps 3.1 and 4.1 give the positive coefficient representation ∣hl∣2=z−2∑jajz−2j≥z−2 and the exact bounds Re⁡tl≤−1/R, ∣Re⁡tl∣≤(l+1)/R, 0<Im⁡tl≤k for the Dirichlet-to-Neumann coefficient.

Depends on

Dependency tree · 2 levels

Sources