Phil Lucht Math & Physics Archive
Home / Math and Physics Files / Math / Special Functions / Legendre Functions

The search for an asymmetric q

DOCX · 81.3 KB
Open DOCX file

Working notes by Phil dated 2.21.10 on Legendre functions Q. He tries to rotate the cut (-∞,1) of Bateman forms (32), (43) and (37), including rotating the cut of the hypergeometric function F via its integral representation, to get a cut-free Q on (-1,1) like Q1(iζ) = ζ cot^-1(ζ) - 1. He concludes the attempt yields nothing new, mostly reproducing the symmetric q function, and gives up. It is motivated by an oblate hyperboloid Green's function problem.

AI-written summary; may contain errors.

Extracted text (machine-read; may contain errors)
The search for an asymmetric q PhL 2.21.10 Goal. Somehow I want to take a cut (-∞,+1) and rotate it down and to the right and thereby remove the cut in the region (-1,1) and by so doing, try to find a q function which solves the Legendre equation, but which is uncut in this region, and therefore generates functions like Q1(iζ)0R = ζ cot-1(ζ) – 1 which are asymmetric in ζ. Cut Structure for Q(32) 1 How to rotate a branch cut. 1 Apply this rotation to our Q (32) formula. 2 Comments on the new function Q0R 3 Cut Structure for Q(43) 3 Comments 4 Cut Structure of Q (37) 5 How do you rotate the cut of F(a,b,c,z) ? 5 Apply to Q (37) 6 Time To Review 8 Cut Structure for Q(32) Let's start with Q (32). If has this form Q = A (z-1)μ/2 (z+1)-μ/2 F(a,b,c, [1-z]/2) + B (z-1)-μ/2 (z+1)μ/2 F(a,b,c', [1-z]/2) What is the cut structure of this animal? To find for F, we plot its arg and we see that for z = (-∞,-1) this makes arg = (1,∞), so we conclude that (-∞,-1) is its cut. The factor (z+1)μ/2 has a branch point at z = -1, so we ignore it because we are interested in z = +1. Each term has a branch point at z = +1, and it is these guys that we want to "rotate" somehow. How to rotate a branch cut. Suppose we have f(z) = (z-1)α0L The idea of rotating the branch cut here was studied in some detail in "study of (1-z)^a and related.doc", Section 3. The function with the cut rotated is called (1-z)α0R and we found this relation: (z-1)α0R = (z-1)α0L Im(z) ≥ 0 (z-1)α0R = e+iα2π (z-1)α0L Im(z) ≤ 0 or (z-1)α0R = (z-1)α0L Im(z) 0 or (z-1)α0R = eiαπ e∓iαπ (z-1)α0L = e+iπ[1∓1]α (z-1)α0L Im(z) 0 Inverting we get (z-1)α0L = e-iπ[1∓1]α (z-1)α0R Apply this rotation to our Q (32) formula. From above we copy down and add the 0L notations Q0L = A (z-1)μ/20L (z+1)-μ/20L F(a,b,c, [1-z]/2) + B (z-1)-μ/20L (z+1)μ/20L F(a,b,c', [1-z]/2) Remember that the F cut is (-∞,-1) which we leave put. We also leave put the z+1 factors. So (z-1)μ/20L = e-iπμ/2 e±iπμ/2 (z-1)μ/20R = e-iπ[1∓1]μ/2 (z-1)μ/20R (z-1)-μ/20L = e+iπμ/2 e∓iπμ/2 (z-1)-μ/20R = e+iπ[1∓1]μ/2 (z-1)-μ/20R We install these to get Q0R = A (z-1)μ/20L (z+1)-μ/20L F(a,b,c, [1-z]/2) + B (z-1)-μ/20L (z+1)μ/20L F(a,b,c', [1-z]/2) = A e-iπμ/2[1∓1] (z-1)μ/20R (z+1)-μ/20L F(a,b,c, [1-z]/2) + B e+iπμ/2[1∓1] (z-1)-μ/20R (z+1)μ/20L F(a,b,c', [1-z]/2) I need to maintain the 0R subscript because this is not the normal function Maple uses! But we can drop the 0L subscript since that is a standard function. Let's now write the above out in full for (32) e-iπμ Qνμ(z)0R = (1/2)f(ν,μ)Γ(-μ) e-iπμ/2[1∓1] (z-1)μ/20R (z+1)-μ/2 F(-ν,1+ν,1+μ, [1-z]/2) + (1/2) Γ(+μ) e+iπμ/2[1∓1] (z-1)-μ/20R (z+1)μ/2 F(-ν,1+ν,1-μ, [1-z]/2) e-iπm Qnm(z)0R = (1/2)f(n,m)Γ(-m) e-iπm/2[1∓1] (z-1)m/20R (z+1)-m/2 F(-n,1+n,1+m, [1-z]/2) + (1/2) Γ(+m) e+iπm/2[1∓1] (z-1)-m/20R (z+1)m/2 F(-n,1+n,1-m, [1-z]/2) Comments on the new function Q0R (a) Since we have merely rotated a cut by 180 degrees, this function still solves the Legendre ODE. (b) the angle arg(z-1) lies in the range (0,2π), whereas angle of arg(z+1) is the usual (-π,π). (c) this formula is no good for μ = 0 (or any integer), so we cannot test it for Q1(z)0R . Cut Structure for Q(43) Q (43) has this form: ( see "zeros and poles..." for gi) e-iπmQnm(z) = 2m-1 g1(n,m) e∓iπ(n+1/2) (z+1)n/2 (z-1)n/2 F(a,b,c; z2/(z2-1)) + 2m g2(n,m) e∓iπ(n-1/2) z (z+1)(n-1)/2 (z-1) (n-1)/2 F(a',b',c'; z2/(z2-1)) As z runs ∞ to 1+ε, the argument in the F function runs 1+ε to ∞. That is the right side red curve below in this plot of the arg versus z. Therefore, the F function part of the functions has a cut (+1,+∞) and (-∞,-1). We can now treat this equation just as we did Q (32) above, using the following rule (z-1)α0L = e-iπ[1∓1]α (z-1)α0R => (z-1)n/20L = e-iπ[1∓1]n/2 (z-1)n/20R (z-1)(n-1)/20L = e-iπ[1∓1](n-1)/2 (z-1)(n-1)/20R Our result is then: e-iπmQnm(z) = 2m-1 g1(n,m) e∓iπ(n+1/2) (z+1)n/2 e-iπ[1∓1]n/2 (z-1)n/20R F(a,b,c; z2/(z2-1)) + 2m g2(n,m) e∓iπ(n-1/2) z (z+1)(n-1)/2 e-iπ[1∓1](n-1)/2 (z-1)(n-1)/20R F(a',b',c'; z2/(z2-1)) Now let's try some phase combining: first term: ∓ iπ(n+1/2) -iπ[1∓1]n/2 = ∓ iπ(2n+1)/2 -iπ[1∓1]n/2 = iπ/2{∓(2n+1)-1 ±n} = iπ/2{∓(n+1)-1} second term: ∓ iπ(n– 1/2) -iπ[1∓1](n-1)/2 = ∓ iπ(2n– 1)/2 -iπ[1∓1](n-1)/2 = iπ/2 { ∓(2n– 1) - [1∓1](n-1) } = iπ/2 { ∓(2n– 1) - (n-1) ±(n-1) } = iπ/2 { ∓(2n) - (n-1) ±(n) } = iπ/2 { ∓(n) - (n-1)} So summarize these two results as first term phase = iπ/2{∓(n+1)-1} = iπ/2{∓n∓1-1} second term phase = iπ/2{ ∓(n) - (n-1)} = iπ/2{∓n - n+1)} Our result is then: e-iπmQnm(z)0R = 2m-1 g1(n,m) eiπ/2{∓n ∓1-1} (z+1)n/2 (z-1)n/20R F(a,b,c; z2/(z2-1)) + 2m g2(n,m) eiπ/2{∓n-n+1) z (z+1)(n-1)/2 (z-1)(n-1)/20R F(a',b',c'; z2/(z2-1)) I was hoping to get something real here, but it is not working right! We still have a cut in (-1,1) which I am trying to get rid of. Comments I have been looking for a Bateman form which can replicate my Q1(iζ)0R = ζ cot-1(ζ) – 1 asymmetric function. But when I came up with this expression starting with Q1(z) = (z/2) ln[(z+1)/(z-1)] - 1, the branch cut I "rotated" was a log one, not a power one. None of the Bateman forms "exposes" this log branch point, it is "hidden inside the F function" in all the forms. So I think this says that I can NEVER arrive at a form like ζ cot-1(ζ) – 1 by trying to do cut rotations in the Bateman forms. Somehow we have to rotate a cut inside the F function. Cut Structure of Q (37) The F functions both have argument w = 2/(1-z) and therefore z = 1 - 2/w. Here are some pictures showing the cut of F(a,b,c,w) in w space and in z space. Here is the usual starting picture: Now if we could rotate the w-space cut to the left, here is what happens: So this would then clean up the (-1,1) region in Q(37). How do you rotate the cut of F(a,b,c,z) ? Consider the integral representation p 114 (2) F(a,b,c,z)0R = K ∫dt tb-1 (1-t)c-b-1 (1-zt)-a0R (1-zt) = t (1/t - z) = K ∫dt tb-1 (1-t)c-b-1 t-a(1/t-z)-a0R Now use the rule which says: (z-1)α0L = e±iπα (1-z)α0R (z-1/t)α0L = e±iπα (1/t-z)α0R => (1/t-z)α0R = e∓iπα(z-1/t)α0L (1/t-z)-a0R = e±iπa(z-1/t)-a0L Then we have F(a,b,c,z)0R = K ∫dt tb-1 (1-t)c-b-1 t-a(1/t-z)-a0R = K ∫dt tb-1 (1-t)c-b-1 t-a e±iπa(z-1/t)-a0L = e±iπa K ∫dt tb-1 (1-t)c-b-1 t-a (z-1/t)-a0L = e±iπa K ∫dt tb-1 (1-t)c-b-1 t-a (z-1/t)-a0L t (z-1/t) = (tz - 1) (*) = e±iπa K ∫dt tb-1 (1-t)c-b-1 (zt - 1)-a0L = e±iπa K ∫dt tb-1 (1-t)c-b-1 (zt - 1)-a We might then define this integral as F(a,b,c,z)0L and from the (*) form you see that it has a branch cut going from +1 to the left, so (-∞,1), which is what I was after: F(a,b,c,z)0L ≡ K ∫dt tb-1 (1-t)c-b-1 (zt - 1)-a = e∓iπa F(a,b,c,z)0R In retrospect, this seems pretty simple. We are just using this phase to flip remove the cut from the right and add it to the left. The question was to find out what phase does the trick. Since F(a,b,c,z)0R is symmetric in a and b, I guess F(a,b,c,z)0L no longer has this symmetry. Apply to Q (37) We showed above that for w = 2/(1-z) we could clear the (-1,1) region in z space if we could take the F cut off to the left. This is Bateman (37) for Q, so write that as e-iπμ Qνμ(z)0R = 2ν [ Γ(1+ν) Γ(1+ν+μ)/Γ(2+2ν)] ) (z+1)μ/2 (z-1)-ν-1-μ/2F(1+ν+μ,1+ν,2+2ν, 2/(1-z) )0R = K (z+1)μ/2 (z-1)-ν-1-μ/2F(1+ν+μ,1+ν,2+2ν, 2/(1-z) )0R Now if we want to "clear out (-1,1)" we have to do two things. First we have to rotate the power cut (z-1)α0R = = e+iπ[1∓1]α (z-1)α0L Im(z) 0 (z-1) -ν-1-μ/20R = e+iπ[1∓1]( -ν-1-μ/2) (z-1) -ν-1-μ/20L Im(z) 0 (z-1)-ν-1-μ/20L = e+iπ[1∓1]( ν+1+μ/2) (z-1)-ν-1-μ/20R Second, we have to flip the F function cut: F(1+ν+μ,1+ν,2+2ν, 2/(1-z) )0R = e±iπ(1+ν+μ) F(1+ν+μ,1+ν,2+2ν, 2/(1-z) )0L Our result is then: = K (z+1)μ/2 (z-1)-ν-1-μ/20R F(1+ν+μ,1+ν,2+2ν, 2/(1-z) )0L e+iπ[1∓1]( ν+1+μ/2) e±iπ(1+ν+μ) Now examine the phase: e+iπ[1∓1]( ν+1+μ/2) e±iπ(1+ν+μ) = eiπ(ν+1+μ/2) eiπ [∓( ν+1+μ/2) ±(1+ν+μ) = eiπ(ν+1+μ/2) eiπ [∓( μ/2) ±(μ)] = eiπ(ν+1+μ/2) e±iπμ/2) I was hoping to get no ± phase so I could claim we are clear on (-1,1). But taking as it comes, we have e-iπμ Qνμ(z)0R = eiπ(ν+1+μ/2) e±iπμ/2) K (z+1)μ/2 (z-1)-ν-1-μ/20R F(1+ν+μ,1+ν,2+2ν, 2/(1-z) )0L Let's just try this for μ = 0 and ν = 1 Q1(z)0R = eiπ(ν+1) K(z-1)-ν-1 F(2,2,4, 2/(1-z) )0L This does seem to have no cut. And we can say: F(2,2,4, 2/(1-z) )0L = F(2,2,4, 2/(1-z) )0R e∓iπ(1+ν) = F(2,2,4, 2/(1-z) )0R e∓iπ2 = F(2,2,4, 2/(1-z) )0R so we then have Q1(z)0R = 2[ Γ(2) Γ(2)/Γ(4)] (z-1)-2 F(2,2,4, 2/(1-z) )0R = (1/3) (z-1)-2 F(2,2,4, 2/(1-z) ) Maple tells us that: Time To Review 1. I wanted to solve the Green's function for an oblate hyperboloid. 2. I constructed a Smythian solution using P(iζ) and Q(iζ) and got a solution, but the solution was complex because these functions are complex. 3. I looked for two independent real solutions of the flipped Legendre equation in ζ. I found some and called them p(ζ) and q(ζ). Both were symmetrical in ζ which I was not happy about. I felt I could not construct an asymmetric solution for the bloid using symmetric functions. The solution has to be asymmetric since the Green's charge is one the upper half and not on the bottom half. 4. I then kept looking for some more general solutions that had asymmetric properties like Q1(iζ)0R = ζcot-1ζ–1. I thought I might find some by "rotating cuts" to the right to "clear out" the (-1,1) region. I then realized that you have to somehow rotate cuts of powers and in addition rotate cuts of F functions. Despite a lot of effort, I came up with nothing new. By considering Q (43), I merely reproduced my symmetric q function. So having done due diligence, I think I will now "give up" on trying to "find" symmetric q-like functions hidden in all those Bateman forms. It is not clear to me how I would have used an asymmetric q function, since it would blow up at least in one direction.